Distributed State Estimation for Dynamic Positioning Systems with Uncertain Disturbances and Transmission Time Delays

Key Laboratory of Cyber-physical System and Intelligent Control in Universities of Shandong, School of Information Science and Electrical Engineering, Ludong University, 264025 Yantai, Shandong, China Yantai Research Institute of New Generation Information Technology, Southwest Jiaotong University, 264025 Yantai, China School of Mathematics and Statistics, Ludong University, 264025 Yantai, Shandong, China


Introduction
For the complex and large-scale systems in the field of energy, transportation, logistics, and so on, due to the emerging problems such as strongly nonlinear and highly coupling, they significantly challenge the current computational tools. To solve science and engineering problems based on complex biological mechanisms and living phenomena, bio-inspired computing (BIC) method is used to design and develop the computer algorithms and models, such as neural networks [1], genetic algorithms [2], and artificial immune systems [3].
Taking the neural network as an example in the field of BIC, Tang et al. proposed a feedforward migration neural network to predict the battery aging trajectories, which arose from the highly nonlinear dynamics of battery aging [4]. In engineering application, the dynamic positioning system of unmanned underwater vehicles (UUVs) is a complex problem mainly due to the nonlinear dynamics, uncertainty in model parameters, and external disturbances. Using the adaptive method based on BIC computing is a promising solution to address these issues. Nevertheless, its full potential is yet to be realized, especially in underwater vehicle applications, mainly due to the highly nonlinear and timevarying nature of vehicle dynamics, unpredictable external disturbances, and the practical difficulty in accurately modeling the hydrodynamic effects [5]. In order to make the processing signal less sensitive to measurement noise, it suffers from scenarios such as the measurement information is inaccurate [6]. However, it is difficult to know the precise noise model in practice. Because of the absence of the full battery degradation model and the inevitable local aging fluctuations in the uncontrolled environments, a base model-oriented gradient-correction particle filter (GC-PF) [7] is proposed to predict aging trajectories of lithium-ion batteries. erefore, as an estimation method that does not need to know the precise model of noise, such as H ∞ filter, which considers the noise as random finite signals, has been developed in engineering applications like fault identification and estimation [6]. It should be noticed that both Kalman filter and H ∞ filter are not suitable for nonlinear and non-Gaussian systems. In nonlinear systems, there are uncertainties caused by modeling errors and external disturbances.
Parameter uncertainties are inevitable in practical systems, and they may lead to severe performance degradation or even instability of the closed-loop system. erefore, they should be fully considered in designing an estimator [8]. Liu et al. employed the Gaussian process regression (GPR) technique to capture the underlying mapping among capacity, storage temperature, and state-of-charge [9]. Moreover, two related data-driven models are developed based on GPR technique [10]. e uncertainties stem from uncertain hydrodynamic parameters, modeling errors, and unknown forces due to the ocean currents in an underwater environment. Specifically, the disturbance observer is developed to estimate the lumped disturbance composed of parametric model uncertainties, modeling errors, and unknown environmental forces. e command governor is formulated as a quadratically constrained quadratic programming problem [11]. Moreover, as opposed to the onestep prediction, multistep prediction is used herein, and a better performance can be achieved [12,13].
As shown in Figure 1, the ship operating in the ocean inevitably suffers the disturbances due to waves, wind, and currents. In practice, the ship dynamics is related to the ship's own characteristics and the operating conditions; thus, there evidently exist parameter uncertainties in the ship motion mathematical model [14]. From Figure 1, it is believed that the concept of hybrid control can provide a scalable and stringent methodology for the design of real industrial control applications. A lot of research studies use the relevant BIC method to deal with several control objectives, changing environmental as well as operational conditions [15]. A ship's speed measurement before GPS is important to the navigation system. In real-time, closedloop control applications, there are problems such as noisy measurements and derivative processes. A variety of filtering methods are available to reduce the noise in the literature [16]. e Kalman filter is applied not only in signal processing but also in the estimation of velocity. e derivation of position signal data may be required not having another sensor for measuring velocity. A general derivative process is the Euler method, but this method may not give accurate results [16,17].
Dynamic networks are interconnected dynamic systems with measured node signals and dynamic modules reflecting the links between the nodes. Various approaches have been developed for identification of dynamic networks, roughly divided into three categories. e first approach considers the identification of a single module in the dynamic network in the situation that the interconnection structure, or topology, of the network is known. e second approach focuses on identification of the full network dynamics for a given topology, and the last category deals with the identification of the topology (and dynamics) of the network [18]. Due to the influence of underwater observation conditions, positioning gross error will often appear; a method of adding the conditional constraints and confidence assessment to EKF was put forward to filter the positioning value, which can make the filtering result more robust and smooth [19]. Sato and Toda [20] proposed adaptive algorithms of gain tuning for Kalman filters and switching Kalman and H ∞ filters for discrete systems. Both gain tuning and switching rely on square means of innovations. Accurate estimation of the horizontal component of the crane force is essential for the feedforward solution, but it is extremely difficult to obtain it from a measurement due to its time-varying nature and the nonlinear system dynamics stemming from the ship-environment interaction [21].
Considering the nonlinear and uncertain dynamic characteristics, a fuzzy neural network backstepping controller is designed to deal with the uncertainties, including estimation and compensation [22]. Brouwer et al. [23] presented an estimator to quantify the random uncertainty of the mean (RUM) from a stationary single time series, without performing actual repeat tests. To develop a robust ship dynamic modeling approach by making efforts on determining a common ship dynamic model and designing a robust identification method [24], the ship dynamic model should be suitable for cases which mean the sufficient balance between the complexity and the accuracy of the model.
Filtering and state estimation are important features of a dynamic positioning system, and estimates of the velocities must be computed from the noisy position and heading measurements through a state observer [25]. To solve the problem of high-precision measurement and positioning, the space measurement and positioning systems with relevant measurement methods have been developed rapidly. With the increasing of the network scale, the scope of the 2 Complexity measurement and positioning systems is more and more large; meanwhile, the coverage mode is more and more flexible. Considering the complexity of the communication links, the energy and computing power of each node are limited. e configuration structure uses the strategy of all nodes, shares, coordinates, and communicates with the signal processing workloads. Accordingly, the network performance does not increase with the size of the network attenuation. erefore, the use of the distributed network structure can have the ability to increase the reliability of information and improve the precision of collaborative calculation. Moreover, the distributed network structure reduces the burden of communication [26]. e distributed estimation strategy is applied to target tracking and positioning, fault diagnosis, and so on. e objective is that state estimation accuracy after fusion is higher than that of each local estimation. Research on filtering is the basis of the information perception and fusion; according to different performance indicators, relevant distributed filtering has aroused more and more attention, including distributed particle filtering [27], distributed H ∞ filtering, and distributed Kalman filtering. Furthermore, distributed Kalman filtering is used for realizing the coordination and fusion of error estimation, which is an important scheme to solve the state estimation for large-scale systems [28]. A data-driven approach combining the RBF neural network (NN) and the extended Kalman filter (EKF) was proposed to estimate the internal temperature for lithium-ion battery thermal management [29,30].
Due to the restriction of communication facilities and scale for networked systems, the measurement information is inevitably affected by the unpredictable interference. Note that the fundamental Kalman filtering method cannot satisfy the performance requirement such that the robust Kalman filtering method is designed considering various noise sequences. Costa and Benites [31] investigated robust pattern independent filtering for discrete-time Markov jump systems with state-correlated noises. Considering the correlation for noise, Feng et al. and Yan et al. [32,33] presented the distributed Kalman filtering fusion method with crosscorrelation between measurement noise and process noise. In addition, for a class of noisy sequences with uncertain variance, Keshavarz-Mohammadiyan and Khaloozadeh [34] used the method of solving the least upper bound of all admissible uncertainties for its true filtering error covariance matrix. e noise variance is set to be suitable in the current system state during the iterative process, and the estimation precision is improved. To deal with multiple time-delay systems and reduce the computational complexity, the received information may be asynchronous in the presence of transmission time delays during network communication.
e method of measurement transformation developed the Kalman filter using the recombination of the innovation sequence. Due to restructuring the measurement sequence, the system with time delay is transformed into the correlative one with free time-delay [35].
In the engineering of the ship navigation applications, due to the complex issue of transmission over networks and external disturbances, the dynamic spatial positioning is difficult with the requirements of high-precision positioning and growing measurement errors. Motivated by the above analysis, based on the principle of spatial positioning of linear charge-coupled devices (CCD), the information measured by the three-dimensional (3D) photoelectric sensor is considered to be influenced by the uncertain factors, such as transmission time delays and cross-correlation noise, and the accuracy of state estimation is improved. In order to improve the robustness for the state estimation, the robust finite-horizon filtering is transformed into Kalman prediction to reduce the computational complexity using the method of reorganizing and transforming the measurement sequence as well as the innovation sequence. e proposed weighted fusion approach of filtering error covariance is used for coordinating and exchanging the measurement information between two subsystems. erefore, based on the optimal information fusion criterion, the distributed weighted fusion estimation using an adaptive selection scheme is the reoptimization of the local estimation, which is able to obtain higher estimation precision than the local one. e remainder of the paper is organized as follows. e dynamic spatial positioning system is modelled in Section 2. Section 3 presents the distributed weighted fusion method using an adaptive signal selection scheme. Simulation results and analysis for the dynamic spatial positioning system are presented in Section 4, and the concluding remarks are given in Section 5.

Dynamic Positioning System
Model e dynamic positioning system is designed by the 3D photoelectric sensor, and the principle of spatial positioning is shown in this section.

Principle of the 3D Photoelectric Sensor.
Firstly, the Metris Krypton measurement approach based on the threecoordinate position measurement is developed. e 3D spatial position of each measuring point is obtained using the optical tracker and the special light emission device (LED). e measured target is made of a special LED, and its 3D coordinates are read from three sets of charge-coupled devices (CCD) on the optical tracker. Special LED placement can be moved arbitrarily, and it can be placed on any moving parts. And then, the processor is used for calculating the position change and direction from the LED at a certain time interval. As shown in Figure 2(a), the 3D photoelectric sensor is composed of three linear CCDs, which are embedded three fixed-focus cylindrical optical mirrors. e three linear CCDs have a "Y" shape distribution on the same plane.
e positioning principle of the 3D photoelectric sensor is shown in Figure 2(b). e measured target is composed of the special LED. And then, each cylindrical optical mirror through light only projects to its corresponding plane under the CCD. A spatial point light source is projected into each CCD plane by a cylindrical optical mirror, which is straight line light intersecting with the CCD and is able to perceive its Complexity 3 position information on the CCD. If the targets are in the measurement range of the 3D photoelectric sensor, it will be decomposed into three planes through the three CCD lenses. erefore, the processor is used to compute the intersection point for the three planes, that is, the coordinates of the light source for the space points (X, Y, Z).

System Model Based on the Adaptive Signal Selection
Scheme. Furthermore, considering the uncertainty of the distributed systems, the system model is established from the system state and measurement information, which is obtained by the sensor during the transmission process with uncertain disturbances and random transmission time delays.
More points in the space are measured by optical spatial positioning systems. It is necessary to use the fusion data from the coordination for multisensors. For the tracking system of the space motion target, the distributed configuration scheme is established, and the optimal estimation is designed to obtain more accurate measurement information. e considered system model for each sensor is described by the following linear discrete-time system [36,37]: Note that x(k) ∈ R r represents the system state, and are known time-varying matrices, and F k are the timevarying parameter uncertainties. w k ∈ R and v i k ∈ R m i denote the process noise and measurement noise of the i-th sensor, which are zero-mean white noise with covariance matrices Q k and R i k , respectively. e initial state x(0) with mean μ 0 and covariance P 0 is assumed to be uncorrelated with other noise signals. Note that the uncertainties F k satisfy F k F T k ≤ I. Due to the limited bandwidth, the transmission from sensors to processors inevitably generates the phenomenon of network congestion. To actively drop the disordered packets from the random transmission time delays, an adaptive signal selection scheme is designed (ASS).
Assumed that the measurement is sampled with a constant sampling period T, the sampling time instant is represented as t ∈ kT, k ∈ N { } for any subsystem. Set the maximum of the transmission delay as N, and it does not exceed the current time instant, i.e., τ(k) ≤ k. To further describe the networked system with random transmission time delays and lost and disordered packets, Figure 3 illustrates a typical scenario. e current time instant k receives a signal with timestamp k 1 , and timestamp k 2 is expressed as the most recent signal before transmitting. From the sensor to the processor, the corresponding transmission time delays are recorded as τ(k 1 ) ∈ N and τ(k 2 ) ∈ N, which denote the transmission time delays from the received packages using the ASS scheme or not at the sampling time k, which satisfies the following relationships: Since the latest data packet during transmission is closer to the current data, a variable of delay β(k) ≥ 0 is developed to express the relationship between k 1 and k 2 , that is, From equations (3) and (4), When receiving a valid data packet with timestamp using the ASS, in order to reduce the computational burden, the stored signal is reorganized as Note that the data packet y ASS (k) is composed of the timestamp and transmission delay at each sampling time [38].

Uncertain Disturbances.
To simplify the modeling of the wave and current for the UUV, considering unknown model parameters and uncertain time-varying disturbances in engineering applications, it is assumed that the bias model and the wave model are driven by zero-mean Gaussian noise. e offset of the vessel in surge, sway, and yaw is small as compared to the size of the ship. Suppose that the distributed systems involve the cross-correlation noise during networked transmission [36,37,39]. Meanwhile, the process for the j-th subsystem are cross-related. e statistical attributes satisfy the following relationships: where the covariance satisfies Note that δ k,l is the Kronecker function. Since the stored measurement sequence contains random transmission time delays and lost and disordered packets, the distributed estimation problem is transformed into solving the optimal state estimation x ASS (k | k) using the ASS scheme, which fuses and compensates each local esti-

ASS-Based Estimation
For the considered distributed uncertain system in equations (1) and (2), a distributed state estimation approach is studied in this section based on robust finite-horizon filtering.

Reorganized Filter and Prediction Vector of the State.
First of all, assume that the stored packet is y i ASS (k) for the ith subsystem at the current sampling time k. Depending on the role of the ASS, the received valid packet is z i (k − τ(k 1 )) with delay τ(k 1 ) before transmission. From equations (2) and (5), the reorganized measurement is denoted as Based on the stored packet sequence y i ASS (0), . . . , y i ASS (k − 1), y i ASS (k) , the optimal state estimation x ASS (k | k) with random transmission time delays can be obtained by the following local recursive method: in which x i ASS (t | t) and x i ASS (t + 1 | t) represent the filter and the predictor, respectively, at the time instant t. Meanwhile, C i ASS,t , K i ASS,t , A i ASS,t , and L i ASS,t are filtering parameters.
In order to solve the upper boundary for the filter and prediction covariance matrix, the augmented vectors of the state are defined as follows: where e i In addition, combining with the upper   Complexity 5 boundary, the augmented system modelled from equations (1) and (8)-(10) is rewritten as and e symbols are defined as Next, the covariance matrices are defined as Θ i (11) and (12). At the same time, the estimated error covariance matrices are derived from equations (11)- (15), and they are as follows: Finally, to find the upper boundary for the error covariance matrices, the filtering covariance is defined from the given equation (10): Note that the corresponding definitions of covariance

Upper Boundary for Error
Covariance. Based on the augmented vector of the state, the corresponding filter and prediction error covariance matrices are designed by eorem 1 to detect the appropriate filtering parameters.  + 1). erefore, the upper boundaries for the error covariance matrices can be calculated by the following recursive formulas: Proof. e proof is similar to the derivation process in [13,38].
And then, to probe the appropriate filtering parameters, the upper boundaries of the error covariance matrices are derived separately: en, the optimal estimation based on the robust finitehorizon filtering in equations (8) and (9) is derived by the following eorem 2.

□ Theorem 2. Suppose that the received valid measurement is y i ASS (k) with τ(k 1 )-step transmission time delays at the current time k, and variable α t is a positive scalar. Matrices P(t) and Σ
i ASS (t) are positive solutions, and they are derived from the following iteration formulas: Therefore, the results of the filtering parameters are solved as follows: From equation (22), we define the following symbols as Proof. e upper boundary evolution is obtained by minimizing the estimation error covariance matrix, which is similar to solving the filtering parameters in [12,13]. □

State Estimation Based on Weighted Fusion Criteria.
For each subsystem, taking into account the network-induced random transmission time delays and lost and disordered packets, the current estimated state is calculated by eorem 2. In order to further improve the accuracy of state estimation, distributed weighted fusion estimation is studied in this section. Based on the weighted fusion criteria, the minimized estimation error cross-covariance matrix is used to perform the information exchange between any two subsystems.
Proof. Its derivation process is similar to the solution of eorem 2.
e optimal weighted fusion estimation can be used to solve the trace of the estimation error by minimizing the cross-covariance matrix [33,40]. Set x i ASS (k | t) as a local estimation with the r-dimensional random vector x(k), which contains the network-induced factor. e distributed 8 Complexity fusion estimator x ASS (k | k) is designed as a r-dimensional vector: It is worth noting that solving x i ASS (k | k) from x i ASS (k | t) develops the one-step prediction method to guarantee the precision of the estimation [35,41].
us, the distributed weighted fusion estimation performance is optimal only if the weighted matrix given by equation (25) meets the following criteria [37,41], i.e., In equation (26), Meanwhile, the corresponding optimal information fusion estimation of the cross-covariance is calculated by So far, the solution process of the distributed weighted fusion estimation method based on the robust finite-horizon filtering using the ASS scheme is shown in Algorithm 1. □ Remark 1. e distributed uncertain systems propose the ASS scheme to deal with the network-induced random transmission time delays and lost and disordered packet phenomena. As mentioned before, for the random transmission delay, the one-step prediction method is used to compensate the missing data packets. On the contrary, in order to realize the information exchange, the distributed state estimation based on the weighted fusion criteria is studied using the linear minimum covariance matrix to obtain the consistent expression of the system state. As a theoretical basis, the fusion scheme has higher estimation accuracy than each local subsystem.

Simulation of Spatial Positioning.
e spatial position structure of the 3D photoelectric sensor is shown in Figure 2(a). Set the length of the three equivalent linear CCDs as 8.8 mm, which are embedded in the 3D photoelectric sensor. e focal length of the column mirror is f � 20 mm, and the distance from the inner end of each CCD to the center of the coordinates is set as 10 mm. e simulation results are demonstrated as follows.

Each Intersection of the Projection Line on the CCD.
It is assumed that a spatial luminous point P0(x0, y0, z0) is used as the spatial positioning point, and the measurable range is 4000 * 4000 * 6000 mm 3 by the 3D photoelectric sensor. P0(x0, y0, z0) passes through each focal point. i.e., f1, f2, f3. And then, the light is divergent from the spatial luminous point. Note that the light is dispersed into a straight line through a cylindrical prism based on the refraction of light principle. e three CCDs project onto the X-Y plane, and it is vertical between each straight line and the corresponding CCD. Since the straight line over the CCD point can be perceived, there are four points which are coplanar for each CCD. at is, the spatial point light, the CCD through the lens of the light point, the lens focus, and the projection point on the X-Y plane. As a result, the three CCDs are extended into three planes, where only the spatial point light is communal. Note that the spatial point light positioning information is the intersection of the three planes. e positioning principle is shown in Figure 2(b).

Measurable Range on the Plane.
e dynamic positioning system is verified using the following three cases: Case 1: according to the principle of the spatial positioning system, every point on the z � 100 mm plane, and the range of 81 × 81 mm 2 , the values are obtained on the spatial plane, and the points perceived by the sensor are the spatial points in the three CCD measurable ranges. As shown in Figure 4, the red area represents the measurable point, and the cyan area represents the undetectable point. , Z � 6000, the space plane is set from Z � 100 mm to Z � 6000 mm, defining the measurable range of each spatial point be 500 mm as the additional steps. Figure 6(a) shows the measurable range for each spatial layer. Similarly, Figure 6(b) describes the results due to reducing the addition steps as 10 mm, i.e., the measurable range from Z � 90 mm to Z � 300 mm on the X-Y plane. Figure 6 shows that the measurable range of spatial points is presented from point, line to triangle, until the polygon, which is dynamically changing.
Similarly, it is obvious that the more precise the pace on the CCD is, the more perceptible the number of points could be sensed. erefore, the spatial positioning system possesses a wider measurement range. From the current simulation results, the range of measurable points in each spatial layer is still Y-axis symmetry, and the distribution profile of each layer point is from one point to the triangle, until the polygon measurement interval. Note that the schematic diagram of the dynamic measurement device is symmetry on the Y-axis as shown in Figure 2(a). According to the principle of the 3D photoelectric sensor, the projection from the spatial point source is also symmetry on each CCD plane.
where the sampling period is set as T � 0.1 s. e maximum value of the transmission time delay is set as 5 sampling periods. e time-varying parameter uncertainty is satisfied as F k � sin(0.6k). System state x(k) � s k _ s k € s k T is used to represent the values of position, velocity, and acceleration for the target at kT. Set parameters as F k � 0.1 0.1 0.1 T , E k � 0.02 0.02 0.02 , C 1 k � 0.6 0.8 1 , C 2 k � 1 0.8 0.5 , C 3 k � 0.3 1 0.7 , and H 1 k � H 2 k � H 3 k � 0.8. Define η k ∈ R as zero-mean white noise with the variance σ 2 η � 0.09. Considering the cross-correlation between the process and measurement noise, ζ i in equation (28) is set as ζ 1 � 2, ζ 2 � 0.8, and ζ 3 � 1, respectively. e covariance of the process noise in equation (28) is w k with the covariance j k Q k , and the mutual covariance is denoted as S i k � ζ i k Q k in equation (28). e initial values are set as x(0 | 0) � μ 0 � E(x(0)) � 1 1 1 T and P(0 | 0) � 0.01I 3 . For the numerical examples, 300 sampling points are selected, and the experimental results are verified based on 300 Monte Carlo simulation calculations. Comparing the traces of the estimated error covariance matrix, the results are shown in Figure 7. e estimated state x(k | k) is validated by the improved robust finite-horizon Kalman filtering (IRFHKF) in [38], as well as the proposed distributed weighted fusion state estimation using an adaptive signal selection scheme (DWFASS). First of all, taking into account the transmission with delay-free from Figure 7(a), the upper boundary of the estimation error Input: e initial state x(0), measurement z 1 (0), . . . , z L (0), estimated state x 0|0 , and error covariance P 0|0 . Output: e local state estimation x 1 ASS (k | k), . . . , x L ASS (k | k) and the fused state estimation x(k | k). (1) Perceive the measurements z 1 (k), . . . , z L (k) from sensor data packets with timestamp at the current time instant k, which are transmitted over the networked systems. (2) Judge the received valid signals using the ASS scheme. e reorganized measurement is denoted as y i (t) for the i-th subsystem.
(3) Probe the filtering parameters C i ASS,t , K i ASS,t , A i ASS,t , and L i ASS,t according to eorem 2. (4) Solve the filter x i ASS (t | t) and predictor x i ASS (t + 1 | t) in equations (8) and (9). (5) Obtain the weighted fusion estimation x(t | t) from equation (25). (6) Update k � k + 1 at the next time instant. Return to Step 2. (8) and (9), respectively. Θ     Complexity covariance matrix is small compared with the IRFHKF method. It is worth noting that the upper boundary of the minimum value for the fused cross-covariance matrix is smaller than the upper boundary for the other covariance matrices, since the weighted fusion criterion is the optimized again in the sense of the minimal variance. Comparing the two methods shown in Figures 7(b)-7(d), the proposed method is more suitable for detecting the optimal filtering parameters. erefore, the designed estimator is appropriate for the uncertainty in model parameters as well as external disturbances, such as random transmission time delays and lost and disordered packets.
In addition, to further analyze the performance of the proposed DWFASS method, the maximum absolute errors (MAE) [42] are obtained by On the one hand, the root mean square error (RMSE) [43] between the measured and estimated values is defined as where M represents the number of the measured points and x(k) and x i ASS (k | k) represent the measured and estimated states at the current sampling points k, respectively. Figure 7(e) shows the MAEs comparing DWFASS and IRFHKF methods; meanwhile, the RMSEs are demonstrated in Table 1. Considering the cross-correlation for the noise, the dynamic tracking trajectory is closer to the actual state value if the noise is more correlated.
Next, for the distributed networked systems, the processing effect of each subsystem is different, which is determined by the cross-correlation for the system parameters and uncertain disturbances. Note that if the process noise is more cross-correlated with the measurement noise, the estimation performance is more accurate.
Moreover, to further illustrate the effectiveness of the DWFASS method, the distributed weighted fusion results for state estimation are shown in Figure 8. e results of state estimation are calculated from eorems 2 and 3, and the method of the robust finite-horizon filtering is studied using the ASS scheme, as well as the distributed weighted fusion estimation criteria. Figure 8 shows that the ASS scheme is able to actively discard data packets owing to the random transmission time delays and lost and disordered packets. At the same time, Table 2 enumerates the RMSEs of each sensor. erefore, the distributed weighted fusion estimation uses the one-step prediction compensation for better performance and computational efficiency in target tracking.

Numerical Simulation Analysis.
According to the simulation results, the dynamic tracking results are compared by the method of recombination error covariance. e comparative results show that the proposed estimation instrument has better estimation accuracy, especially in the case of communication constraints with multistep random delays and lost and disordered network induction packets. Because the true estimation error covariance is lower than the upper boundary of the probe, the proposed distributed fusion estimation method has better performance in the phenomenon of network induction, including the ability to closely track the state of the system.

Conclusion
is paper, facing the spatial positioning principle of the linear CCD, introduces the workflow of the 3D photoelectric sensor dynamic positioning system and simulates the spatial positioning effect of the proposed DWFASS method. According to the measured spatial target, the trajectory of its motion is tracked and estimated, and a kind of distributed uncertain system with network-inducing factors is studied. In order to process the lost and disordered data packets by the random transmission time delays, the system is modelled based on the adaptive signal selection scheme. In order to overcome the influence of limited communication capacity and reduce the computational burden, this paper adopted the onestep prediction compensation strategy and proposed the optimal state estimation method to compensate for the missing packets. At the same time, the finite-horizontal filtering strategy is designed, and the real estimated error variance obtained by the DWFASS method is less than its upper boundary. en, the distributed estimation method based on filtering is introduced, which introduces the uniform expression of the weighted fusion standard. erefore, the distributed fusion estimation had higher estimation accuracy than that of each local subsystem. e numerical simulation is demonstrated using the spatial target positioning system containing three sensor devices, and the results show that the proposed distributed weighted fusion estimation has better tracking performance. To better process the data required for positioning and improve the accuracy of positioning, this paper has optimized the performance of positioning algorithms. e proposed method can be used in various motion control scenarios such as target tracking, trajectory tracking, path following, and formation control of AUVs subject to uncertainties and constraints. In addition, the future work will further focus on following some known distributions for the time delays or packet loss to further enhance the state estimation, such as the Gaussian or binary distribution based on the known bandwidth of the network.

Data Availability
e data used to support the findings of this study are included within the proposed information fusion estimation approach in the article. e data from numerical simulations

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.