Abstract

For the multisensor linear discrete time-invariant stochastic control systems with different measurement matrices and correlated noises, the centralized measurement fusion white noise estimators are presented by the linear minimum variance criterion under the condition that noise input matrix is full column rank. They have the expensive computing burden due to the high-dimension extended measurement matrix. To reduce the computing burden, the weighted measurement fusion white noise estimators are presented. It is proved that weighted measurement fusion white noise estimators have the same accuracy as the centralized measurement fusion white noise estimators, so it has global optimality. It can be applied to signal processing in oil seismic exploration. A simulation example for Bernoulli-Gaussian white noise deconvolution filter verifies the effectiveness.

1. Introduction

An important application background of white noise estimation is signal processing in oil seismic exploration. After the explosives buried underneath earth surface are exploded, analyzing the reflection coefficient formed by the reflections of every oil layer can be used to determine whether there is oil underground and the geometry shape of the oil field. The reflection coefficient can be described by Bernoulli-Gaussian white noise. So the white noise estimation problem can be used in the oil seismic exploration. This problem has been deeply researched by Mendel [13], but Mendel has not solved the system’s measurement white noise estimation problem and also has not referred the problem of multisensor information fusion white noise deconvolution estimation.

In order to improve the estimation accuracy of single sensor to the white noise, in [4, 5], when input noise and measurement noise are not correlated, the multisensor information fusion white noise deconvolution filters are put forward, respectively, by Kalman filtering method and modern time series analysis method. In [68], multisensor information fusion white noise optimal filter is presented for the systems with correlated noise by Kalman filtering method. In [9], modern time series analysis method and Gevers-Wouters algorithm are used to present the information fusion white noise deconvolution filter for multisensor systems with correlated noises. Thus, the solution of Riccati equation is avoided and the self-tuning filter with unknown model parameters can be designed. The shortcomings of methods presented in [39] are to require computing the high-dimension cross-covariance matrix and the fusion accuracy is global suboptimal.

Recently, the weighted measurement fusion (WMF) method has gained great attention. Its basic principle is to weigh local sensor measurements to obtain a low-dimensional measurement equation according to some fusion criterion and then use a single Kalman filter to obtain the final fused state estimation. It can, not only reduce the computing burden greatly, but also give a global optimal estimation. This can be explained as the accuracy of WMF and the centralized measurement fusion (CMF) filter [10] is the same. Therefore, it is globally optimal. Gan and Chris [10] put forward the WMF algorithm with the assumption that all the sensors have the same measurement matrix and the measurement noises of each sensor are uncorrelated. Using the Lagrange multiplier method, the WMF algorithm is presented when the measurement noises of each sensor are correlated [11]. Ran et al. [12, 13] put forward the WMF algorithm when the extended measurement matrix of all the sensors has full column rank or the measurement matrices of all the sensors have the maximal right factor. Self-tuning WMF Kalman filtering algorithm is presented in the work of Gao et al. [14] and Ran and Deng [15]. In [16, 17], full-rank decomposition and weighted least square theory is used under correlated noises and different sensor measurement matrices; the WMF algorithm is presented.

However, using WMF method to solve the white noise estimation value problem of multisensor system with correlated noise and with different measurement matrices in each sensor is always a difficult issue to be solved, since the present white noise estimation theory is not suitable for WMF method. In this paper, we use the WMF algorithm to solve the white noise fusion estimation problem of multisensor systems with correlated noises and different measurement matrices. Firstly, under the assumption that the noise input matrix is of full column-rank, we present the CMF white noise estimators by the extended measurements of all sensors. They have a large computing burden due to the high-dimension measurement matrix. Further, the WMF white noise deconvolution estimators are presented to reduce the computing burden, which have the global optimality.

The paper is structured as follows. The CMF white noise deconvolution estimators are presented in Section 2. The multisensor WMF white noise deconvolution estimators are designed in Section 3. A simulation example follows in Section 4. Some conclusions are given in the end.

2. Multisensor CMF White Noise Deconvolution Filter

Consider the discrete time-invariant linear stochastic control systems with sensors where is the state, , are the measurements, is the known control input, , , and are constant matrices with compatible dimensions. is the measurement matrix of the sensor . and are correlated white noises with zero means, and where the symbol denotes the expectation, , (), , the superscript denotes the transpose. Combining measurement equations of (2.2) yields the centralized measurement equation: where

The fusion measurement white noise has variance matrix . The correlated function of and is .

To convert the systems (2.1) and (2.4) into the uncorrelated system, (2.1) is equivalent to where is a pending matrix. (2.6) can be converted into where , , . as output feedback becomes a part of the control item. Then primary system formulae (2.1) and (2.4) are equivalent to the system formed by formulae (2.4) and (2.7). To make , introduce which ensures that and are not correlated. Then variance matrix of is yielded as .

Lemma 2.1 (see [18, 19]). Multisensor systems (2.1) and (2.4) with correlated noise have CMF global optimal input white noise deconvolution estimators and the error variance matrices as satisfies Riccati equation: Smoothing gain has two computing methods:

Theorem 2.2. Supposing is full column-rank matrix, for systems (2.1) and (2.4), the optimal white noise deconvolution estimator of CMF input white noise is given by And the variance matrices of estimator errors are given as defining where is the pseudo-inverse of , that is,

Proof. Notice when N 0, , then (2.18) holds obviously. For (2.1), according to the projection theory [19], we have From (2.11) and (2.24), we have which can be simplified to Then, from the definition of (2.23), we easily obtain (2.19).
Subtracting (2.25) from (2.1) yields the error relation Notice ; there is a prediction error relation Substituting (2.29) into (2.28) yields And notice , then we have So, (2.21) holds. From (2.18), we have (2.20). The proof is completed.

Theorem 2.3. For (2.1) and (2.4), when is a full column-rank matrix, one has optimal white noise deconvolution smoothers of CMF input white noise as The error variance matrices are computed by Define

Proof. For (2.1), the projection theory is used, and we have From (2.15), we have From the recursive relation of (2.36), we have then we obtain (2.32). In fact, from (2.17), we have which is simplified to Define as Then we obtain (2.34). From the definition of (2.35), (2.32) is proved. Using minus both sides of (2.32), and from , (2.33) is proved. This completes the proof.

Corollary 2.4. For multisensor systems with correlated noise (2.1) and (2.4), the nonrecursive white noise smoothers are given by The error variances satisfy

Proof. (2.32) and (2.33) are iterated by times; (2.42) and (2.43) are obtained. This completes the proof.

3. Multisensor WMF White Noise Deconvolution Estimator

From [20], we know that any nonzero matrix has full-rank decomposition where is a full column-rank matrix with the rank and is a full row-rank matrix with the rank r, then measurement model (2.4) can be represented as Given that is a full column-rank matrix, it follows that is nonsingular. Then the weighted least squares (WLS) [21, 22] method is used and the Gauss-Markov estimate of is yielded as

then substituting (2.4) into (3.3) yields

The variance matrix of is given by .

For systems (2.7) and (3.5) using standard Kalman filtering algorithm [19], we can obtain WMF Kalman estimators , and its variance matrices , innovation , . It is proved in [16] that the WMF Kalman filter has the global optimality; that is, it is numerically identical to the CMF Kalman filter if they have the same initial values.

The above WMF method can obviously reduce the computing burden since the dimension of the measurement vector for the CMF is , , while that for the WMF is , and is much larger than generally.

Theorem 3.1. For (2.7) and (3.5), when is a full column-rank matrix, one has WMF optimal nonrecursive smoothers of input white noise The error variance matrices satisfy where If is satisfied, then WMF is numerically equivalent to CMF, that is, where white noise filter is

Proof. From [16], when , we have So, we have . (2.19), (2.42), and (2.43) are compared with (3.7), (3.8), and (3.12), then (3.11) are obtained.

Corollary 3.2. WMF input white noise recursive smoothers are The error variance matrices satisfy then, one has the following relation:

Proof. From Theorem 2.3 and (3.13), we have (3.14)–(3.17). This completes the proof.

4. Simulation Example

Consider the multisensor discrete linear stochastic ARMA signal system where is the signal, , are the measurement signals. , are Gaussian white noises with zero mean and variance matrix . And , where is Bernoulli white noise satisfying if , and if , where denotes probability. is independent of , then the variance matrix of is . is the back shift operator.

Our goal is to find the optimal white noise deconvolution estimators , , and , , and the corresponding error variance matrices , , and , to test the estimation result values of these three algorithms equal and to compare the computing burden.

The system (4.1) is converted to state space models

where we define

In the simulation, we set

So

It is not full column rank, so Hermite standard model [20] can be used to compute the full-rank decomposition of , , where

The simulation results are presented by Figures 1, 2, 3 and 4, Tables 1 and 2. In Figures 1, 2, 3 and 4, the -coordinates of endpoints on real lines represent real values , -coordinates of solid round points’ centers represent , ,1,2,3, -coordinates of hollow round points’ centers represent , . The -coordinates of hollow squares’ centers represent , . It can be seen that the results are completely unanimous.

In Figure 1, it can be seen that at some specified time, the results of these three algorithms are also completely the same. Figure 2 represents the sameness of estimation error variances of the three fusion methods. The fusion accuracy is better than local estimation accuracy of every single sensor.

These results all represent that WMF is completely equivalent to CMF, and the fusion accuracy is higher than local estimation accuracy of each single sensor.

In the other aspect, it can be seen that when the filtering algorithm of CMF is used, the dimension of is , the Kalman filtering needs to compute the inverse matrix of a matrix at every time, while when the WMF Kalman filtering algorithm is used, the dimension of fusion measurement is , then the matrix inversion computation is converted into scalar division, so the WMF method can obviously reduce the computational burden.

5. Conclusions

White noise deconvolution problems have great applications background in oil seismic exploration. Under the condition that the white noise input matrix is full column-rank, the centralized measurement fusion and weighted measurement fusion white noise estimators are presented based on the projection theory, respectively. Their function equivalence has been proved. Furthermore, the proposed weighted measurement fusion white noise estimators can obviously reduce the computational burden.

Acknowledgments

This work is supported by Science and Technology Research Foundation of Heilongjiang Province, Education Department under Grant 12521420.