Distributed Fusion Estimation for Multisensor Multirate Systems with Stochastic Observation Multiplicative Noises

This paper studies the fusion estimation problem of a class of multisensor multirate systems with observation multiplicative noises. The dynamic system is sampled uniformly. Sampling period of each sensor is uniform and the integer multiple of the state update period. Moreover, different sensors have the different sampling rates and observations of sensors are subject to the stochastic uncertainties of multiplicative noises. At first, local filters at the observation sampling points are obtained based on the observations of each sensor. Further, local estimators at the state update points are obtained by predictions of local filters at the observation sampling points. They have the reduced computational cost and a good real-time property. Then, the cross-covariance matrices between any two local estimators are derived at the state update points. At last, using thematrix weighted optimal fusion estimation algorithm in the linear minimum variance sense, the distributed optimal fusion estimator is obtained based on the local estimators and the cross-covariance matrices. An example shows the effectiveness of the proposed algorithms.


Introduction
In networked systems or sensor networks, there often exist various uncertainties during the transmission process of signals due to the imperfection of the communication channels.It makes impossible to use linear model to describe some systems.The uncertainties can be approximated mathematically by an additive noise or a multiplicative noise [1][2][3][4][5][6].These systems are widely used in petroleum seismic exploration, target detection, speech processing, and other areas; thus, the research on systems with multiplicative noise has the important practical significance.In the early references [1], the optimal linear filters have been proposed for systems with uncertain observations described by the multiplicative noise.For more general case with stochastic parameters, the optimal linear estimation is designed in [2].References [3][4][5] study the polynomial filters; however, the proposed nonlinear filters have expensive computational cost.For networked systems with multiplicative noises and packet dropouts, optimal linear estimators including filter, predictor, and smoother have been proposed in [6].However, the above-mentioned literatures are all concerned with single sensor case but do not take multiple sensors into account.
As the sensor technology is widely used in military, civilian, scientific research, and many other fields, single sensor has failed to meet the performance requirements in many aspects.Moreover, as the development of electronics technologies, various sensors have been developed and applied to many practical fields such as target tracking since they can provide more information than any single sensor.Therefore, multisensor information fusion has received considerable research attention in recent years [7].For systems with a single sampling rate, the optimal state weighted fusion filter in the linear minimum variance sense [8] and the self-tuning fusion filter with unknown noise variances [9] have been presented.Recently, the multirate multisensor asynchronous fusion algorithms have been studied in [10][11][12].References [13,14] adopt the state augmentation approach to give the estimators with the expensive computational cost.Though [15,16] adopt the nonaugmented approach to design the filters, a modeling error is made by ignoring the process noise.Therefore, there is the accuracy loss.By considering the process noise to eliminate the modeling error, an optimal filter is presented to improve the estimation accuracy [17].Furthermore, the missing measurements are also taken into account in [16,18].In [19], a multirate multisensor distributed fusion estimator is proposed for two-sensor systems with one-step cross-covariance noises.However, most of the above-mentioned literatures do not take the multiplicative noises into account.In sensor networks, there often exist various sensors with different sampling rates and stochastic uncertainty of multiplicative noises.It is significant to use the nonaugmented approach to deal with the multirate multisensor systems.This motivates our work.
This paper studies the fusion estimation problem of a class of multisensor multirate systems with observation multiplicative noises.State is sampled uniformly at the finest rate.Different sensors have different sampling periods that are integer multiples of the state update period.Local estimators at the state sampling points are obtained based on the local filters at the observation sampling points by the filtering and prediction.By using the distributed optimal weighted fusion estimation algorithm in the linear minimum variance sense [8], a distributed optimal fusion estimator is obtained.It avoids the state and observation augmentation.It has a good reliability since it has the distributed parallel structure.Moreover, the estimation error cross-covariance matrices between any two local estimators are derived according to the different sampling cases.

Problem Formulation
Consider the following multisensor multirate system with observation multiplicative noises: where () where the symbol  is the mathematical expectation.
Remark 1.The sampling case of multisensor multirate systems can be described by Figure 1.The horizontal axis denotes time while the vertical axis denotes different sensors.Three sensors are shown in Figure 1.Black circle solid points represent the sampling time of different sensors.The sampling rate goes from the highest (sensor 1) to the lowest (sensor 3).As shown in Figure 1, the three sensors all sample uniformly.The first sensor has the same sampling rate as the state update rate; that is, the sampling period is .The sampling period of the second sensor is 2 and the third is 3.It is clear that the least common multiple of three sample periods is 6.This means that the samplings of different sensors are asynchronous in each data block of the length 6.
The objective of this paper is to find the distributed optimal fusion estimator x () of () based on the local estimators x () from different sensors.
To obtain the distributed fusion estimator by using the optimal weighted fusion estimation algorithm in the linear minimum variance sense [8], we need to compute the local estimators and variance matrices from each sensor and the cross-covariance matrices between any two local estimators.In the latter text, we will give the computation of local estimators and cross-covariance matrices.

Local Filters at the Observation Sampling Points
At first, we give the filter at the observation sampling points of each sensor.
From the iteration of (1), we have Mathematical Problems in Engineering ; we have the state space model at the observation sampling points for the th sensor as follows: with the noise statistical information Then, we have the filter at the observation sampling points of each sensor based on the above model.Lemma 2 (see [6]).For system (4), the local filters at the observation sampling points of the th sensor are computed by where

Local Estimators at the State Update Points
Based on the filters at the observation sampling points in Lemma 2, we have the following state estimation algorithms at the state update points.
Remark 4. The local estimators at the state update points have been obtained by filtering and prediction based on the filter at the observation sampling points.State augmentation is avoided.They are simple and have a good real-time property.Now, we have obtained the local estimators at the state update points based on the observations of each sensor.Next, we compute the cross-covariance matrices between any two local estimators.

Computation of Cross-Covariance Matrix
Theorem 5.The estimation error cross-covariance matrices between any two local estimators can be computed in the following three cases.
(i) If the th sensor and the th sensor both have observations at time , the estimation error cross-covariance matrix is given as (ii) If the th sensor has an observation and the th sensor does not have any observation at time , the estimation error cross-covariance matrix is given as (iii) If both the th sensor and the th sensor do not have observations at time , the estimation error cross-covariance matrix is given as The initial value is   (0) =  0 . Proof. where Then we obtain the estimation error equation as ( The cross-covariance matrix can be computed by where Substituting ( 17 Then, we have the estimation error equation as Similarly to the derivation of the case (i), (11)

Distributed Fusion Estimator
In the preceding sections, we have obtained the local estimators at the state update points and their covariance matrices.Applying the distributed matrix weighted optimal fusion estimation algorithm in the linear minimum variance sense [8], we can obtain the distributed fusion estimator as follows: The optimal weighted matrices are computed by where () = (  ()) is an  ×  matrix whose (, ) block is   () and  = [  , . . .,   ]  is an  ×  matrix.Then, the optimal fusion estimation error variance matrix is computed by Furthermore, we have   () ≤   ().
Remark 6.Compared to the centralized fusion estimator, the distributed fusion estimator has the flexibility, fault tolerance, and reliability since it has the distributed parallel structure [8].

Simulation
An uninterruptible power system (UPS) with three sensors subject to the multiplicative noises is taken as an example to demonstrate the effectiveness and applicability of the proposed method.We consider the UPS with 1KVA.The discrete-time model (1) can be obtained with sampling time 10 ms at half-load operating point as follows [20]: where (), V  (  ), and   (  ) are uncorrelated white noises with zero mean and variances   ,  V  , and    .In simulation, we take   = 0.64,   1 = 1.2,   2 = 0.7,   3 = 0.3,  V 1 = 3,  V 2 = 2,  V 3 = 1,  1 = 1,  2 = 2,  3 = 3,  11 = [8 12 6],  12 = [10 8 5],  13 = [8 6 10], the initial values (0) = 0, and  0 = 0.1 3 .The sampling case is shown in Figure 1. Figure 2 gives the distributed fusion estimator.We see that the fusion estimator has the effective estimation performance.Figure 3 gives the comparison of variances of the distributed fusion estimator and local estimators.We see that the proposed fusion estimator outperforms the local estimators.Figure 4 gives the comparison of variances of the distributed fusion estimator and centralized fusion estimator.We see that the distributed fusion estimator has the small accuracy loss.However, it is significant that the distributed fusion estimator has better reliability than the centralized fusion estimator since it is convenient to detect and isolate the faults of sensors from distributed structure.

Conclusion
A distributed fusion estimator has been designed for systems with multiple sensors of different sampling rates and observation multiplicative noises.Compared with the centralized

Figure 3 :
Figure 3: Comparison of variances of distributed fusion filter and local filters.

Figure 4 :
Figure 4: Comparison of variances of centralized and distributed fusion filters.
is the system state at time  and  is the state update period.Φ, Γ,  0 , and  1 are constant matrices.(  ) is the observation of the th sensor at time   ;   is the ratio of the observation sampling period and the state update period.is the number of sensors.() and V  (  ) are white noises with zero mean and variances   and  V  , respectively.The observation multiplicative noise   (  ) is scalar white noise with zero mean and variance    .(), V  (  ), and   (  ) are uncorrelated with each other.The initial state (0) is independent of (), V  (  ), and   (  ) and satisfies that x (   |   ) and x (   |    −   ) are the filter and predictor at the observation sampling points, respectively.  (   |   ) and   (   |    −   ) are the corresponding covariance matrices.  (  ) is the innovation sequence with the covariance matrix   (  ).  (  ) is the filtering gain matrix.  (  ) is the state secondorder moment.The initial values are x (0 | 0) =  0 ,   (0 | 0) =  0 , and   (0) =  0   0 +  0 .