A Recursive Non-Uniform Sampling Estimator for Asynchronous Nonlinear Systems

In this paper, we consider the problem of asynchronous estimation in the presence of packet losses for the randomly sampling nonlinear system. Packet losses occur at the control input and at the measurement side. Firstly, the synchronization of the asynchronous sampling system is realized by weighting the state of the adjacent state update points. Secondly, the projection theorem is used to estimate the system state at the sampling time. Due to modeling errors and unmodeled dynamics, obtaining an accurate dynamic model is challenging. Therefore, observation inference based on interpolation techniques is proposed to solve the asynchronous estimation problem. Furthermore, the algorithm is extended to multi-sensor systems to obtain a distributed fusion estimator. Finally, simulation experiments are conducted to validate the effectiveness of the algorithm.


Introduction
Compared to the uniform sampling mechanism, the non-uniform sampling mechanism can better ensure the controllability and observability of the discretization system [1], which is widely used in battery health management [2], target tracking [3] and power systems [4].Therefore, it is of great theoretical and practical significance to study the estimation problem of networked asynchronous sampling systems.The difficulty in solving this problem lies in synchronizing asynchronous sampling systems.After decades of development, the current synchronization methods for asynchronous non-uniform sampling systems are mainly divided into data block [5] and model transformation [6].
The introduction of data block theory provides an intuitive approach to asynchronous estimation problems.Representative work in this area is as follows.For the estimation problem of randomly sampling a neural network with an unknown sampling interval, this paper designs a density-identifiable estimator based on data reception [7].The non-uniform sampling mechanism also poses a challenge for parameter identification in nonlinear systems.The decomposition technique effectively reduces the computational complexity of the algorithm [8].In networked control systems, there are many uncertain factors, such as time delays, packet losses, fading measurements, random parameter uncertainty, correlated noise, and more.Using innovation analysis and the covariance intersection (CI) technique, an asynchronous fusion estimator that accounts for fading measurements is proposed [9].Event-triggered mechanisms are crucial for achieving time synchronization in asynchronous sampling systems.A dynamic event-triggered mechanism [10] is used to develop the time-varying moving horizon asynchronous estimator.An event-triggered estimator is developed by incorporating packet losses and heavy-tailed noise into the system model.Based on covariance intersection, a fast CI event-triggered fusion estimator is proposed, and a sufficient condition for error bounded [11] is derived.
However, if the data block is large, this augmentation method imposes a large computational burden.Therefore, the method based on state weighting has gradually attracted the attention of researchers [12].A high-gain state estimator is proposed for nonlinear asynchronous continuous-discrete sampling systems.The convergence and stability of the estimator are proven by analyzing the Riccati equation [13].As the research progresses, multiplicative noise and missing measurements have also attracted attention.The estimation problem of asynchronous randomly sampling systems with packet losses is studied, and a centralized fusion algorithm based on measurement data reordering is proposed [14].Using a neural network compensation model, the optimal fusion estimator based on fault detector is proposed [15].In recent years, the study of event-triggered mechanisms has become a hot topic due to the need to conserve communication resources [16].The term "optimal estimation" refers to the best possible estimation of the state or parameter based on a specific estimation criterion when the measurements are known.Based on the minimum entropy error criterion, a sequential fusion estimator is proposed to solve the asynchronous sampling problem with communication delay, and the convergence proof is provided [17].Unlike traditional estimators designed to handle network-induced phenomena, an estimator that proactively designs communication protocols to mitigate potential non-ideal situations has garnered attention.The estimation problem for asynchronous sampling systems with repetitive processes is considered [18].
Since packet losses are inevitable due to network congestion and connection disruptions, it is necessary to consider packet losses in the estimation problem of asynchronous sampling systems.In practice, nonlinear systems offer more advantages than linear systems in describing the dynamic evolution process of the system.Therefore, the nonlinear asynchronous sampling fusion estimation problem holds greater research value.Compared with the augmentation method, measurement fusion has lower computational complexity, and its accuracy is close to that of the centralized method [19].Considering the characteristics of asynchronous sampling in network systems, sequential methods are more suitable for dealing with asynchronous fusion problems [20].In addition, the unmodeled dynamics, such as highly maneuvering systems [21] and anti-modeling [22], are easily overlooked when establishing state models, and the established models are also susceptible to modeling errors.Therefore, an observation inference method based on interpolation is proposed to solve the asynchronous estimation problem.In summary, the main contributions of this paper are as follows: 1.
For the state estimation problem of a nonlinear asynchronous sampling system, the synchronization of the state space model is realized by weighting two adjacent state update points at the measurement sampling point.

2.
Based on the concept of probability density distribution of an approximate nonlinear function, a deterministic sampling Gaussian estimation framework is proposed, which is considered superior to the algorithm framework with an approximate nonlinear function.

3.
Considering the potential modeling errors and unmodeled dynamics in the state equation, a non-recursive estimator based on interpolation operations that solely depends on measurement data is proposed, whose advantage lies in the minimal computational load and rapid calculation speed.
The paper is organized as follows.The system model, the basic theory, and the interested problem are presented in Section 2. The state estimator for an asynchronous sampling system is proposed, and an observation inference scheme independent of the state model is presented in Section 3. The distributed fusion estimator adapted to the multi sensor system is given in Section 4. Section 5 shows the performance analysis of the estimator.Section 6 provides simulation experiments.Section 7 concludes this paper.
Notation R n means the n-dimensional Euclidean space.A T and A −1 stand for the trace and determinant of a matrix A. P (•) and E(•) represent the probability and the expectation.

System Architecture
The asynchronous non-uniform sampling stochastic system with packet dropout is characterized by.
where x k ∈ R n denotes the state, u k ∈ R p means the control input, y k i ∈ R m is the measurement, and x k i denotes the state at sampling points y k i .The statistical properties of the noise w k ∈ R r and v k i ∈ R m are provided.Φ, B and Γ are matrices with appropriate dimensions.
The measurement function h(•) is the nonlinear smooth function.Random variables following a Bernoulli distribution are used to describe packet dropout.Specifically, ξ k and γ k i represent the data losses of control input and measurement end.
Figure 1 illustrates the general process of information flow in the networked control system.The initial state x 0 with mean µ 0 and variance P 0 is uncorrelated with the other random variables.Such a linear dynamic nonlinear observation system is common in the field of target tracking [23] and robot localization.The measurement data are usually obtained in the polar coordinate system or spherical coordinate system.For example, the distance and angle of the target are measured in the polar coordinate system, while the state equation is often established in the rectangular coordinate system.The former and the latter are undoubtedly nonlinear relationships [24].The non-uniform sampling mechanism is illustrated in Figure 2, where the state is updated uniformly, and measurements are randomly sampled.
the expectation.

System Architecture
The asynchronous non-uniform sampling stochastic system with packet dropout is characterized by.
The measurement function ( )   Figure 1 illustrates the general process of information flow in the networked control system.The initial state 0 x with mean 0  and variance 0 P is uncorrelated with the other random variables.Such a linear dynamic nonlinear observation system is common in the field of target tracking [23] and robot localization.The measurement data are usually obtained in the polar coordinate system or spherical coordinate system.For example, the distance and angle of the target are measured in the polar coordinate system, while the state equation is often established in the rectangular coordinate system.The former and the latter are undoubtedly nonlinear relationships [24].The non-uniform sampling mechanism is illustrated in Figure 2, where the state is updated uniformly, and measurements are randomly sampled.The difference between non-uniform sampling and uniform sampling but random packet loss lies mainly in two aspects: the sampling strategy and the effects of packet loss.Firstly, in terms of the sampling strategy, non-uniform sampling means that the data points are not collected according to a fixed time interval or frequency, but rather based on some non-uniform or random pattern.This can lead to an uneven distribution of data, with some areas being dense and others sparse.In contrast, uniform sampling but random The difference between non-uniform sampling and uniform sampling but random packet loss lies mainly in two aspects: the sampling strategy and the effects of packet loss.Firstly, in terms of the sampling strategy, non-uniform sampling means that the data points are not collected according to a fixed time interval or frequency, but rather based on some non-uniform or random pattern.This can lead to an uneven distribution of data, with some areas being dense and others sparse.In contrast, uniform sampling but random packet loss involves sampling at a fixed time interval or frequency, but the sampled packets may be randomly lost during transmission or processing.Secondly, regarding the impact of packet loss, although there is random packet loss in both cases, the effect on data quality will vary due to different sampling strategies.In the case of non-uniform sampling, the uneven distribution of data can be exacerbated by packet loss, leading to loss or difficulty in analyzing important information.In the case of uniform sampling, even though the data itself are uniformly distributed, random packet loss can still lead to discontinuity or missing information.

Actuator Plant Sensor
Several measurements are distinguished here for the convenience of deriving subsequent algorithms.The measurement under the non-uniform sampling mechanism is The non-uniform sampling measurement considering the effect of packet loss is The successfully received measurement at the estimator is defined as Next, we begin with the measurement equation and propose an estimation algorithm based on pure measurement data by employing interpolation techniques and inversion operations.

Gaussian Filter at Measurement Sampling Points
The pseudo-state equation is calculated using the state weighting.The recursive estimator at the random sampling point can be derived using the projection theorem.
Let k i and k i−1 be the sampling times of the i-th and (i − 1)-th successfully received measurements of the intervals (k − 1, k] and (l − 1, l], respectively.Define The state at a time instant k − α k i can be described This formula describes two scenarios for the position of the sensor output sampling points.In the first case, the sensor output is sampled as soon as the state is updated; therefore, the system state corresponds to the state measured by the sensor.It denotes that α k i = 0 and x k−α k i = x k .In the second case, the sensor output is sampled between the state update period, which means that 0 ≤ α k i < 1 and Similarly, the state at the sampling time l − α l i−1 can be described.
Lemma 1.The state equation at the random sampling point can be obtained through state weighting where The covariance matrix of the system noise is calculated by Remark 1. Unlike existing modeling methods, the approach used in this paper does not increase the dimension of the measurement Equation ( 12), hence it is referred to as a non-augmented state space model.State weighting refers to the weighted sum of two adjacent state update points to construct the new state, as shown in (5).After synchronization, the state equation can be obtained by iteratively updating the original state equation to maintain consistency with the measurement equation over time 14.However, in reality, there is no state update at the random sampling point, so the state equation is referred to as a pseudo-state equation.
After obtaining the synchronization model, the state estimators at the sampling points are determined.

Theorem 1.
For the state space model ( 7) and ( 8), the recursive asynchronous estimator is given as follows.
The state predictor and its covariance matrix are computed The state and noise gain matrices are calculated State and noise estimators are calculated.
where innovation is computed.
Sensors 2024, 24, 2882 6 of 20 The covariance matrices for state and noise are given.
Proof of Theorem 1. First, we need to prove the time updating part.The prediction error equation can be obtained by Substituting (22) into leads to (14).
Calculating the gain matrix K x k−α k i by definition (15) requires first calculating the covariance matrices Similarly, since the innovation covariance matrix requires computing only the covariance matrix where the cross-covariance matrix for noise and state is computed by Next, the proof of the measurement update part is presented.Based on the definition of the filtering error equation, The prediction error covariance matrix (20) can be obtained by The projection theorem demonstrates that the noise smoother Sensors 2024, 24, 2882 7 of 20 possesses a structure similar to that of the state estimator, allowing the derivation of the smoothing error equation.
Similarly, ( 21) and ( 33) can be obtained through the cross-covariance matrix P w The corresponding nonlinear estimator can be obtained through various deterministic sampling methods.For instance, the Unscented Kalman Filter (UKF) can be derived through unscented transformation, while the Cubature Kalman Filter (CKF) can be obtained by using the spherical radial cubature rule.The non-augmented model-based UKF and CKF algorithms proposed in this section are referred to as NM-UKF and NM-CKF algorithms.□ Remark 3.However, in reality, obtaining an accurate state model is challenging due to the presence of unmodeled dynamics.This, in turn, leads to a performance deterioration of conventional filters that rely on the accuracy of the model.Therefore, a measurement inference scheme based on interpolation is also proposed.
First, measurement data at the state update point can be obtained using the piecewise cubic Hermite interpolation algorithm.(It is worth noting that, compared with cubic spline interpolation [25], the modified Akima piecewise cubic interpolation (MAPCI) produces fewer fluctuations and is more suitable for handling rapid changes between platform regions.Compared to the shape-preserving piecewise cubic interpolation [26], MAPCI is less sharply flattened, which enables it to effectively handle oscillating data.)The next step is to estimate the state using measurements that are synchronized with the state updates.In other words, the state is inferred directly from the measurements, which are corrupted by noise [27].(The first scenario involves an invertible observation function, but the observation noise is unknown.In this case, the observation noise may be colored noise, correlated noise, non-Gaussian noise, or other complex situations.Since it cannot be accurately described, we choose to ignore the noise to simplify the problem.The second scenario occurs when the observation function is invertible, and the observation noise is known [28].If the system model adheres to the assumptions of the classical Kalman filter, then the inversion operation can be implemented directly.If the observation function is nonlinear and the observation noise is complex [29], it may be more reasonable to sample the observation noise and use the mean and variance of a set of samples to characterize the state estimate.The third scenario occurs when the observation function is not invertible.Facing this situation in practical engineering, a common solution is to enhance observability by increasing the number of sensors [30].The addition of sensors to address the irreversibility of the observation function mainly depends on the number of observation variables).A summary of the pure measurement methods used to solve the multirate estimation problem can be found in Appendix A.
The block diagram of our proposed Gaussian estimator with the O 2 inference based on Gaussian Hermite interpolation can be seen in Figure 3.The concept behind both methods Sensors 2024, 24, 2882 8 of 20 is to convert the asynchronous sampling system into a synchronous sampling system and then use an appropriate synchronous system scheme to derive the state estimation.The synchronization concept of the Gaussian estimator relies on the timing of measurement sampling to achieve synchronization through state weighting.Subsequently, a Gaussian estimation framework is introduced based on the concept of a probability density function that approximates the nonlinear function.The synchronization concept of O 2 inference based on Gaussian Hermite interpolation involves aligning the state update time through Gaussian Hermite interpolation, followed by deriving the state estimate using O 2 inference.observation variables).A summary of the pure measurement methods used to solve the multirate estimation problem can be found in Appendix A.
The block diagram of our proposed Gaussian estimator with the O2 inference based on Gaussian Hermite interpolation can be seen in Figure 3.The concept behind both methods is to convert the asynchronous sampling system into a synchronous sampling system and then use an appropriate synchronous system scheme to derive the state estimation.The synchronization concept of the Gaussian estimator relies on the timing of measurement sampling to achieve synchronization through state weighting.Subsequently, a Gaussian estimation framework is introduced based on the concept of a probability density function that approximates the nonlinear function.The synchronization concept of O2 inference based on Gaussian Hermite interpolation involves aligning the state update time through Gaussian Hermite interpolation, followed by deriving the state estimate using O2 inference.

Distributed Fusion Estimator
In the previous section, we proposed an asynchronous sampling state estimation algorithm.In this section, we provide an asynchronous fusion estimator that is suitable for multi-source systems.
For multi-sensor systems, the state equation is (1).The measurement equation of the multi-sensor random sampling system can be given as follows.
( ) ( ) where represents the ith measurement of sensor j, and When considering the multi-sensor system, an intuitive approach is to rearrange all measurements.This approach has the advantage of treating the system to some extent as a single sensor system.In this case, the proposed estimation algorithm in this paper can be used to solve the problem of asynchronous fusion estimation.However, the drawback of this scheme is that it requires all sensors to function properly, which means that the reliability of this scheme is not optimal.In other words, if any sensor fails, the estimator

Distributed Fusion Estimator
In the previous section, we proposed an asynchronous sampling state estimation algorithm.In this section, we provide an asynchronous fusion estimator that is suitable for multi-source systems.
For multi-sensor systems, the state equation is (1).The measurement equation of the multi-sensor random sampling system can be given as follows.
When considering the multi-sensor system, an intuitive approach is to rearrange all measurements.This approach has the advantage of treating the system to some extent as a single sensor system.In this case, the proposed estimation algorithm in this paper can be used to solve the problem of asynchronous fusion estimation.However, the drawback of this scheme is that it requires all sensors to function properly, which means that the reliability of this scheme is not optimal.In other words, if any sensor fails, the estimator diverges.To solve the issue of asynchronous sampling in fusion estimation, we propose distributed fusion estimation algorithms.
To avoid calculating the cross-covariance matrix [31], the CI fusion algorithm is adopted, which gives an upper bound on the error variance matrix by solving the op-timization problem online.The state estimator obtained by the CI fusion algorithm can be expressed as where P CI k|k is the error covariance matrix of the fusion estimator.ω l k 0 ≤ ω l k ≤ 1 is the weighting coefficient, which is obtained by solving the following optimization problem min The distributed fusion estimator presented in this paper can be adapted to non-ideal situations because some of its operations can be performed within the sensor.Meanwhile, it has good fault tolerance, flexibility, and reliability due to its low requirement for channel capacity.

Performance Analysis
We first analyze the estimation performance when there are no packet losses.Lemma 2. Assuming that the estimator xk−1|k−1 and variance matrix P k−1|k−1 of state x k−1 have been obtained, xk|k and xk|k represent the estimator of state x k at time k obtained by using N k or N k data in the period (k − 1, k] , respectively.The corresponding variance P k|k and P k|k .So, if N k ≥ N k ≥ 0 , then P k|k ≤ P k|k .Lemma 2 states that the more fully the measurement data are used over time, the more accurate the estimation result.It also shows that when the local estimation error is bounded, the fused estimation error is also bounded.Therefore, in the following, we analyze the performance of the local estimator Definition 1 (Stochastic Boundedness [32]).
Theorem 2. For matrix pair Φ, √ Q to be stabilizable, define The error covariance matrix is given by where It follows from (38), that we have Substituting ( 42) into ( 40) and ( 41) yields Theorem 2. □

Simulations 6.1. Target Tracking System
State estimation algorithms are widely used in target tracking, industrial control, wireless communication, and other fields.Among these systems, the target tracking system is often used to validate the effectiveness of estimation algorithms.The motion model is a crucial component of a target tracking algorithm.Most current target tracking algorithms rely on the motion model.In the research into target tracking algorithms, several simple and easily implementable motion models have been the focus of development work.The constant acceleration model is one of the most common basic motion models.The rangebearing model of active radar is one of the most common measurement models, where measurements include the range and bearing, which are affected by noise.Therefore, this paper uses the following discrete-time asynchronous random sampling system as the verification example in the simulation section.
where The range-bearing measurement equation is given by.
where the active sensor is located at S x , S y .Set ξ = 0.8, γ 1 = 0.8, γ 2 = 0.9.Similar to system noise, the variances in measurement noise v 1 Figure 4 illustrates the position tracking effect of the estimators for sensor 1.The accumulative mean square error (AMSE) of the estimation algorithms for sensor 1 is plotted in Figure 5 and the performance for different estimators is given in Table 1.The core of the Extended Kalman Filter (EKF) lies in the first-order linearization of the nonlinear function.Therefore, for plants with strong nonlinearity, the filtering effect is not ideal, or even unusable.Both the UKF and CKF adhere to the idea that "approximating the probability density function distribution of a nonlinear function is easier than approximating the nonlinear function itself," and they belong to the Gaussian filtering framework.The UKF algorithm obtains the sigma points and corresponding weights through unscented transformation.The weights are generally negative in high-dimensional systems, which introduces highorder truncation error terms and reduces the accuracy of the algorithm [33].The CKF algorithm acquires cubature points and propagates them using nonlinear equations.As a result, the weights are always positive, leading to a significant reduction in error.Therefore, in the simulation results, the CKF algorithm demonstrates higher numerical stability and filtering accuracy compared to the UKF algorithm.
However, obtaining the exact state equation is very challenging due to unmodeled dynamics in practical engineering.Therefore, this paper proposes a measurement inference method that does not rely on the state equation.First, timestamp registration is achieved through interpolation.Second, the system state is inferred directly through the O 2 inference.This process avoids the estimation error caused by unmodeled dynamics in the state equation.According to the simulation results, the MAPCI-O 2 algorithm and NM-CKF algorithm demonstrate equivalent tracking performance.Moreover, since the MAPCI-O 2 algorithm relies solely on measurement data, the running time and computational burden of the algorithm are significantly reduced compared to traditional filtering algorithms.
Figure 6 illustrates the feasibility of the proposed estimation algorithm.The AMSEs of local and fusion estimation algorithms are shown in Figure 7 and a performance comparison for local and fusion estimators is given in Table 2.The accuracy of the CI fusion estimation algorithm based on multi-sensor measurement data is higher than that based on singlesensor data.6 illustrates the feasibility of the proposed estimation algorithm.The AMSEs of local and fusion estimation algorithms are shown in Figure 7 and a performance comparison for local and fusion estimators is given in Table 2.The accuracy of the CI fusion estimation algorithm based on multi-sensor measurement data is higher than that based on single-sensor data.
Figure 8 shows the error curve of the estimator with different N for sensor 1, where 1 N / = N means that the sensor output is sampled once in the N state update cycle.
Simulation results indicate that increasing the number of sampling times during the same state update period leads to a more accurate estimation effect.The recursive estimation algorithm proposed in Section IV is a general framework, which means that most nonlinear filtering algorithms based on deterministic sampling can be applied within this framework.The UKF and CKF are typical representatives of deterministic sampling used in simulation experiments.Furthermore, the estimation error    6 illustrates the feasibility of the proposed estimation algorithm.The AMSEs of local and fusion estimation algorithms are shown in Figure 7 and a performance comparison for local and fusion estimators is given in Table 2.The accuracy of the CI fusion estimation algorithm based on multi-sensor measurement data is higher than that based on single-sensor data.
Figure 8 shows the error curve of the estimator with different N for sensor 1, where 1 N / = N means that the sensor output is sampled once in the N state update cycle.
Simulation results indicate that increasing the number of sampling times during the same state update period leads to a more accurate estimation effect.The recursive estimation algorithm proposed in Section IV is a general framework, which means that most nonlinear filtering algorithms based on deterministic sampling can be applied within this framework.The UKF and CKF are typical representatives of deterministic sampling used in simulation experiments.Furthermore, the estimation error  Figure 8 shows the error curve of the estimator with different N for sensor 1, where N = 1/N means that the sensor output is sampled once in the N state update cycle.Simulation results indicate that increasing the number of sampling times during the same state update period leads to a more accurate estimation effect.
The recursive estimation algorithm proposed in Section IV is a general framework, which means that most nonlinear filtering algorithms based on deterministic sampling can be applied within this framework.The UKF and CKF are typical representatives of deterministic sampling used in simulation experiments.Furthermore, the estimation error is used as the performance indicator.Simulation examples demonstrate that the NM-CKF algorithm exhibits the best performance, while the NM-UKF algorithm is slightly inferior, but both algorithms are capable of tracking the system state.However, the NM-EKF algorithm is unable to track the system state.This highlights the necessity and effectiveness of designing the proposed algorithm.When the state equation contains significant uncertainty or is difficult to obtain, the MAPCI-O 2 algorithm can be used as an effective complement to the traditional state estimator.In addition, in Section V, we extend the state estimation algorithm to include asynchronous sampling with packet losses in the multi-sensor system.Simulation results further confirm the effectiveness of the fusion algorithm.
is used as the performance indicator.Simulation examples demonstrate that the NM-CKF algorithm exhibits the best performance, while the NM-UKF algorithm is slightly inferior, but both algorithms are capable of tracking the system state.However, the NM-EKF algorithm is unable to track the system state.This highlights the necessity and effectiveness of designing the proposed algorithm.When the state equation contains significant uncertainty or is difficult to obtain, the MAPCI-O2 algorithm can be used as an effective complement to the traditional state estimator.In addition, in Section V, we extend the state estimation algorithm to include asynchronous sampling with packet losses in the multi-sensor system.Simulation results further confirm the effectiveness of the fusion algorithm.is used as the performance indicator.Simulation examples demonstrate that the NM-CKF algorithm exhibits the best performance, while the NM-UKF algorithm is slightly inferior, but both algorithms are capable of tracking the system state.However, the NM-EKF algorithm is unable to track the system state.This highlights the necessity and effectiveness of designing the proposed algorithm.When the state equation contains significant uncertainty or is difficult to obtain, the MAPCI-O2 algorithm can be used as an effective complement to the traditional state estimator.In addition, in Section V, we extend the state estimation algorithm to include asynchronous sampling with packet losses in the multi-sensor system.Simulation results further confirm the effectiveness of the fusion algorithm.

Spring Mass System
The spring mass system, which is used as a practical example in many studies, is modeled as

Spring Mass System
The spring mass system, which is used as a practical example in many studies, is modeled as where x k = [x 1,k , x 2,k ] T , x 1,k and x 2,k are the position and velocity of mass.κ is the spring constant.c is the damping coefficient.α is the sensor scale factor.The system parameters are set as follows: M = 1, κ = 1, c = 0.5 and α = 1.The packet losses phenomenon, which is considered a common network induced phenomenon in networked control systems, is described as the Bernoulli distribution.Set ξ = 0.8, and Q w = diag( [1,1]).The initial state is x 0 = [0, 0] T .Set γ = 0.8, and Figures 9 and 10 illustrate the tracking results of different estimators, while Table 3 displays the performance of different estimators.Similar to the conclusion for the target tracking system, the three estimators can approximately track the target.Specifically, the performance of UKF and CKF is similar, while the MAPCI-O 2 algorithm is slightly outperforms the Gaussian estimator, which may be attributed to the fact that the disturbance induced by measurement noise is partially mitigated by random sampling in the simulation.
To illustrate the impact of the packet loss phenomenon on the estimator, it is necessary to conduct a packet loss rate experiment, and it should be noted that this paper considers both the packet loss ξ k from the controller to the actuator and the packet loss γ k from the sensor to the estimator (SE).However, if ξ k changes, it will cause the system model to change, which means that the performance comparison of the estimator loses the benchmark at that time.To this end, we mainly focus on the packet loss phenomenon on the SE side and use NM-CKF as the base algorithm, Figure 11 shows the estimation performance of the estimator under different packet loss rates.It can be seen that the larger γ is, the smaller packet losses probability is, and the estimation performance of the filter is better.This means that the phenomenon of packet losses does affect the estimation performance of the filter.The primary focus of this paper is to solve the state estimation problem in asynchronous non-uniform sampling systems.The paper does not employ a compensation strategy to rectify the accuracy losses resulting from packet losses.The existing packet loss compensation strategies mainly include zero compensation, hold compensation, and prediction compensation [34].Among these, the prediction compensation strategy has the best effect.The main reason is that the predicted value of the measurement based on the state prediction is closer to the current time measurement than the zero and the last time measurement [35].In addition to filters, there are also observers that can help solve the state estimation problem, among which the intermediate observer is particularly advantageous for resolving packet loss issues [36].This also motivates us to focus on the fusion estimation problem for asynchronous non-uniform sampling with data packet loss in our future work.Specifically, we aim to investigate more effective and general packet loss compensation strategies and synchronization methods.
In order to verify the applicability of the covariance crossover idea, the distributed fusion estimator experiment is carried out next, where the simulation parameters are γ 1 = 0.9 for sensor 1 and γ 2 = 0.8 for sensor 2, and the rest of the parameters are the same as above Figure 12 displays a comparison diagram of estimations between the local estimator and the fusion estimator.It can be seen that the accuracy of the fusion estimator is better than that of the local estimator, indicating that the proposed algorithm remains effective in the context of multi-sensor systems.From the local estimator, it can be seen that a smaller packet loss rate leads to less information loss and higher accuracy of the estimator, as also depicted in Figure 13.
include zero compensation, hold compensation, and prediction compensation [34].Among these, the prediction compensation strategy has the best effect.The main reason is that the predicted value of the measurement based on the state prediction is closer to the current time measurement than the zero and the last time measurement [35].In addition to filters, there are also observers that can help solve the state estimation problem, among which the intermediate observer is particularly advantageous for resolving packet loss issues [36].This also motivates us to focus on the fusion estimation problem for asynchronous non-uniform sampling with data packet loss in our future work.Specifically, we aim to investigate more effective and general packet loss compensation strategies and synchronization methods.
In order to verify the applicability of the covariance crossover idea, the distributed fusion estimator experiment is carried out next, where the simulation parameters are 1 0.9 γ = for sensor 1 and 2 0.8 γ = for sensor 2, and the rest of the parameters are the same as above Figure 12 displays a comparison diagram of estimations between the local estimator and the fusion estimator.It can be seen that the accuracy of the fusion estimator is better than that of the local estimator, indicating that the proposed algorithm remains effective in the context of multi-sensor systems.From the local estimator, it can be seen that a smaller packet loss rate leads to less information loss and higher accuracy of the estimator, as also depicted in Figure 13.In general, multi-sensor fusion can usually be divided into three architectures: centralized [37], decentralized [38], and distributed [39].In centralized fusion, the observation equations of all sensors are combined into a high-dimensional observation equation using the dimension expansion method, and then Kalman filtering is applied simultaneously with the state equation.Although centralized fusion can theoretically achieve the global optimal fusion estimate, it comes with the drawbacks of high computational burden and poor fault tolerance [40].Decentralized fusion uses a weighted local Kalman filter to obtain the fusion estimate, so it is also referred to as weighted fusion.It has the advantages of low computational burden, easy fault diagnosis, and separation.However, such architectures typically necessitate highly reliable sensors, which are often In general, multi-sensor fusion can usually be divided into three architectures: centralized [37], decentralized [38], and distributed [39].In centralized fusion, the observation equations of all sensors are combined into a high-dimensional observation equation using the dimension expansion method, and then Kalman filtering is applied simultaneously with the state equation.Although centralized fusion can theoretically achieve the global optimal fusion estimate, it comes with the drawbacks of high computational burden and poor fault tolerance [40].Decentralized fusion uses a weighted local Kalman filter to obtain the fusion estimate, so it is also referred to as weighted fusion.It has the advantages of low computational burden, easy fault diagnosis, and separation.However, such architectures typically necessitate highly reliable sensors, which are often costly and not easily scalable.In addition, due to physical constraints, such as communication delays, the communication bandwidth is limited, and the fusion center cannot effectively communicate with all sensors in a large-scale sensor network [41].Each sensor node in the distributed architecture uses only the information obtained from its locally connected neighbors for fusion.This feature can offer improved built-in redundancy compared to the other two types of architecture, thereby enhancing robustness against sensor failures.The distributed fusion architecture can also reduce the communication burden because the function of the fusion center is amortized over the nodes [42].Therefore, this can improve flexibility.The main indicators used in the algorithm evaluation include global optimality, local consistency, fully distributed and communication burden [43], and the selection of different performance indicators results in different advantages and disadvantages of the algorithm, which also encourages us to conduct more comprehensive and detailed research on the fusion architecture in the future.

Conclusions
We investigate the asynchronous multirate estimator with packet dropout.The state equation at the sampling point is determined based on state weighting.The estimators at the sampling point of the measurement and the state update point are determined using the orthogonal projection theorem and innovation analysis technique.In addition, if the state equation cannot be obtained, a state estimation algorithm based on measurement inference is proposed.The fusion form of the proposed algorithm suitable for a multi-sensor system is provided, and the feasibility of the algorithm is confirmed with simulation experiments.

4.
Utilizing measurements to obtain the interpolation function.

R
are provided. , B and  are matrices with appropriate dimensions.
h is the nonlinear smooth function.Random variables following a Bernoulli distribution are used to describe packet dropout.Specifically,  k and  i k represent the data losses of control input and measurement end.

Figure 1 .
Figure 1.Networked control system under transmission control protocol.Figure 1. Networked control system under transmission control protocol.

Remark 2 .
If the sensor output is sampled as soon as the state is updated, the filter xk−α k i |k−α k i becomes the estimator xk|k at time k.If the measurement data are not sampled at time k but within the period (k − 1, k], the predictor xk|k−α k i based on the filter xk−α k i |k−α k i is used as the estimator xk|k at time k.If the sensor output is not sampled in time after state update, the predictor xk|k−1 is used as the estimator xk|k at time k.

Figure 3 .
Figure 3.Comparison of Gaussian estimator based on measurement sampling time and O2 inference based on state update time.

Figure 3 .
Figure 3.Comparison of Gaussian estimator based on measurement sampling time and O 2 inference based on state update time.

Proof of Theorem 2 .
It follows from the Bayesian formula γ and velocity, respectively.u k = [u 1,k , u 2,k ] T and w k = [w 1,k , w 2,k ] T denote the control input and system noise.The variances of the random variables w 1,k and w 2,k are Q w 1 = 0.09 and Q w 2 = 0.25.The state update period ∆ = 1s.Let N be the sensor output sampling times in the state update cycle.

Figure 6 .Figure 7 .
Figure 6.Estimation effect of local and CI fusion estimation algorithms.(a)

Figure 6 .
Figure 6.Estimation effect of local and CI fusion estimation algorithms.(a) p x k .(b) p y k .

Figure 6 .Figure 7 .
Figure 6.Estimation effect of local and CI fusion estimation algorithms.(a)

Figure 8 .
Figure 8.Comparison of AMSE of local estimator 1 with different N.

Figure 11 .Figure 10 .
Figure 11.Estimation performance of filters with different packet loss rates.

Figure 11 .Figure 11 .
Figure 11.Estimation performance of filters with different packet loss rates.

Figure 12 .
Figure 12.Tracking performance of different estimators.

Figure 12 .
Figure 12.Tracking performance of different estimators.

Figure 13 .
Figure 13.Estimation performance of local and fusion estimators.

Figure 13 .
Figure 13.Estimation performance of local and fusion estimators.

Table 1 .
Performance of different estimators.

Table 1 .
Performance of different estimators.

Table 1 .
Performance of different estimators.

Table 2 .
Performance of fusion estimator and local estimators.

Table 2 .
Performance of fusion estimator and local estimators.

Table 2 .
Performance of fusion estimator and local estimators.

Table 3 .
Performance comparison for different estimators.
P (x) = 3hs 2 − 2s 3 h 3 y k+1 + h 3 − 3hs 2 + 2s 3 h 3 y k + s 2 (s − h) h 2 d k+1 + s(s − h) 2 h 2 d k (A4) where h = x k+1 − x k , s = x − x k .5.Constructing the measurement data y k at the state update point x k involves utilizing the interpolation function, thus completing asynchronous sampling synchronization.6.If the observation function is invertible.7.If the observation noise statistics are known.The state estimation χk can be obtained through χi k = h −1 y k , υ i is one of the samples sets sampled from the observation noise v k .Otherwise, the observation function is not invertible The number of sensors has been increased to enhance the observability of the plant, enabling the acquisition of state estimation χk y 1,k = h 1,k (x k , v 1,k ) y 2,k = h 2,k (x k , v 2,k ). ..y n,k = h n,k (x k , v n,k ) (A8)