Discrete Dynamics in Nature and Society
Volume 2018 (2018), Article ID 1516028, 12 pages
https://doi.org/10.1155/2018/1516028

Research Article

Kalman Filtering Algorithm for Systems with Stochastic Nonlinearity Functions, Finite-Step Correlated Noises, and Missing Measurements

Yonghui He,¹ Jibin Jiang,1,2 Hischuan Huang,³ Shufang Zhuo,¹ and Yanfeng Wu¹

¹Department of Mechatronics Engineering, Fujian Polytechnic of Information Technology, Fuzhou 350003, China
²School of Mechanical and Automotive Engineering, Fujian University of Technology, Fuzhou 350108, China
³Department of Automation Engineering, Chienkuo Technology University, Changhua 50094, Taiwan

Correspondence should be addressed to Yonghui He; auto_mitu@163.com

Received 8 August 2017; Revised 15 March 2018; Accepted 8 April 2018; Published 15 May 2018

Academic Editor: Manuel De la Sen

Copyright © 2018 Yonghui He et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The locally optimal filter is designed for a class of discrete-time systems subject to stochastic nonlinearity functions, finite-step correlated noises, and missing measurements. The multiplicative noises are employed to describe the random disturbances in the system model. The phenomena of missing measurements occur in a random way and the missing probability is characterized by Bernoulli distributed random variables with known conditional probabilities. Based on the projection theory, a class of Kalman-type locally optimal filter is constructed and the filtering error covariance matrix is minimized in the sense of minimum mean square error principle. Also, by solving the recursive matrix equation, we can obtain the filter gain. Finally, two examples are provided: one is a numerical example to illustrate the feasibility and effectiveness of the proposed filtering scheme; the other is to solve the problem of target estimation for a tracking system considering networked phenomena.

1. Introduction

The filtering theory is a kind of optimal approach to estimate the state of the target plant based on the measurement output of the observation signals and the certain filtering criteria. Kalman filtering theory is a modern filtering theory which was proposed by Kalman in [1]. Kalman filter is an optimal estimation technique, and it is also a recursive minimum variance estimation method. From the view of geometry, the estimation of the Kalman filtering can be regarded as the projection of the state variable on the linear space generated by the observation sequence. For the linear system with the Gauss distribution noises, the recursive least mean square estimation of the system state can be obtained according to the Kalman filtering method. The state space of modern control theory is introduced into the optimal filtering theory by Kalman filtering where the dynamic model of the system is described by the equation of state and the observation equation is used to describe the observation model. In addition, the Kalman filtering can deal with the problems of time-varying systems, nonstationary signals, and multidimensional signals. Subsequently, the recursive estimation method was developed by employing the minimum mean square error principle and the projection theory. Accordingly, a considerable number of research results were reported to design the optimal filters in many practical systems [2–7].

Due to the influence of certain factors, the phenomena of correlated noises and missing measurements are inevitably encountered and the accuracy of system state estimation is degraded. Hence, the problems of correlated noises and missing measurements attracted the attention of many scholars and a large number of results were published in [8–13]. Recently, in [14], the problem of state estimation was solved for a class of linear discrete-time stochastic systems subject to missing data and correlated noises, where the estimators were unbiased and the estimation error covariances were minimized. Based on the state augmentation approach and the projection theory, in [15], the optimal linear estimator was designed in the minimum variance sense for linear uncertain systems with correlated noises and incomplete measurements and the sufficient condition on the existence of steady-state estimators was shown. In [16], the distributed and centralized fusion filtering problems were discussed for the sensor networked systems with stochastic uncertainties and correlated random transmission delays, where the recursive algorithms for the optimal linear distributed and centralized filters were derived by the innovation sequence analysis method. Subsequently, in [17], according to the minimax robust estimation principle, the robust centralized fusion and weighted measurement fusion time-varying Kalman estimators were presented in a unified framework for uncertain multisensor systems subject to linearly correlated white noises.

In many practical systems, the stochastic nonlinear functions are commonly encountered due to the unreliability of the communication network [18–23]. Hence, it is necessary to deal with the stochastic nonlinear functions to improve the accuracy of system state estimation. To mention but a few, in [24], a finite-horizon recursive filter was designed for discrete-time-varying nonlinear system with stochastic nonlinearities and randomly occurring measurement delays which made the filtering error minimum at each sampling time. The probability-constrained filtering problem was investigated in [25] for a class of time-varying nonlinear stochastic systems subject to estimation error variance constraint and a new filter was constructed to guarantee a minimized upper-bound on the estimation error variance. In [26], by solving two Riccati-like difference equations, an easy-implemented recursive filtering algorithm was investigated for a type of time-delayed stochastic nonlinear systems subject to event-based communication protocols and missing measurements. Based on the linear matrix inequality approach, in [25], a filter was designed to solve the problem for a class of nonlinear systems with stochastic sensor saturations and Markovian measurement transmission delays and the sufficient conditions were established to guarantee the asymptomatic stability of the filtering process. However, to the best of the authors’ knowledge, these researches do not pay much attention to the problem of the optimal filtering for a class of discrete-time systems with stochastic nonlinearity functions, finite-step correlated noises, and missing measurements.

Based on the above discussions, in our paper, the purpose is to discuss the locally optimal filtering problem for a class of discrete-time systems with stochastic nonlinearity functions, finite-step correlated noises, and missing measurements. Due to the unreliable factors, the multiplicative noises are employed to describe the random disturbances in the system model. The considered phenomenon of missing measurements is modeled by a random variable satisfying Bernoulli distribution with known conditional probabilities. A class of Kalman-type locally optimal filter is designed and the filtering error covariance matrix is minimized in the sense of minimum mean square error principle by using the projection theory. By solving the recursive matrix equation, we can obtain the filter gain. Finally, two examples are given to illustrate the feasibility and effectiveness of the proposed filtering scheme. The contribution of this paper: (1) The system model is considered to add the multiplicative noises, stochastic nonlinearity, correlated noises, and missing measurements. (2) We make first attempt to design a locally optimal filter for the systems with complex network phenomena (multiplicative noises, stochastic nonlinearity functions, correlated noises, and missing measurements). (3) A new recursive algorithm is established to obtain the locally optimal filter which is suitable for online applications.

The rest of the paper is arranged as follows. In the next section, the models are proposed and a brief description of the systems is presented. In Section 3, a class of Kalman-type locally optimal filter is designed for discrete-time systems with stochastic nonlinearity functions, finite-step correlated noises, and missing measurements. In Section 4, an illustrative example is used to show the effectiveness of the proposed filtering method. Finally, the conclusions are shown in Section 5.

Notations 1. The symbols used in the paper are standard. denotes the -dimensional Euclidean space. represents the transpose of matrix . stands for the expectation of the variable . and represent the identity matrix and the zero matrix with appropriate dimensions, respectively. stands for a diagonal matrix with elements in the diagonal. is the Kronecker delta function. If , then ; else, . If the dimensions of the matrices are not definitely stated, they are considered to be well-matched for algebraic operations.

2. Problem Formulation and Preliminaries

Consider the following discrete-time stochastic systems with multiplicative noises and stochastic nonlinearity functions:

where is the state vector; represents the measured output. , , , and are zero-mean Gaussian white noises with unity covariances and , are employed to describe the multiplicative noises, which are uncorrelated with other noise signals. The matrices , , , , and are known matrices with appropriate dimensions. The random variable , which describes the phenomenon of missing measurements, obeys the Bernoulli distribution and has the following statistical properties:

where is a known scalar and is uncorrelated with other noise signals.

In this paper, we consider the process noise and observation noise being correlated Gaussian white noises and satisfy the following equations:

where , , , and is the number of the correlated steps. and are uncorrelated with other signals.

The functions and denote the stochastic nonlinearities with , and satisfy the following expectations and variances for all :

where is a known nonnegative integer; , , and are known matrices with appropriate dimensions.

The purpose of this paper is, based on the projection theory and the minimum mean square error principle, to design the Kalman-type recursive filter as the following form:

where and are the state estimation and the one-step prediction at time instant . is the filter gain matrix to be determined.

Remark 2. In model (2), if , , which represents the fact that the sensor receives the data at time instant successfully and we can estimate the state by using (7) and the observation values; if , , which stands for the fact that the sensor only receives the values of the stochastic nonlinearity function and the noises of the time instant . In this case, the sensor occurs with the missing measurements and the state is estimated only by using the values of the stochastic nonlinearity function and the noises. Equation (7) of the state estimation becomes with the following form:

Remark 3. Due to the presence of complex network phenomena (multiplicative noises, stochastic nonlinearity functions, correlated noises, and missing measurements) in the target system model, it is difficult to explain whether the designed filter is globally optimal. The general algorithms are focused on the local solution to reduce the amount of computation and the complexity of the algorithm for complex problems. Hence, a recursive filtering algorithm is proposed to design a locally optimal filter (6)-(7) for the target plant (1)-(2) based on the observation sequence .

3. Main Results

In this part, we aim to derive the filter gain matrix to obtain the class of Kalman-type locally optimal filter by applying the projection theory and the filtering error covariance matrix is minimized in the sense of minimum mean square error principle. Therefore, the following lemmas and definitions are introduced for promoting the subsequent derivation process.

Lemma 4. Define and , where . Then, the following equations can be calculated:

Proof. For , based on the definition of , it can be held as follows:

Since

then, .
By using the same method in (10), we have

Similarly, can be derived.
Thus, from the above proof, (9) and (12) are true. The proof of this lemma is complete.

Lemma 5. Letting the state covariance be , then one has the following equation:

Proof. By the definition of and Lemma 4, we can establish the following equation:

Hence, (13) can be held. The proof of this lemma is complete.

Lemma 6. Define the matrices and . Then, the following equations are given:

Proof. Based on (6) and (7) and Lemma 4, we have the following recursive equation:

Here, we assume that is -step correlated with the noise; then, . Thus,

Similarly, can be derived. The proof of this lemma is complete.

Lemma 7. Define the one-step prediction error and the filtering error . Then, the covariance matrix is computed as follows:

where

is calculated by Lemma 6 and is the filtering error covariance.

Proof. Substituting (1) and (6) into the definition of yields

By utilizing (7) and Lemmas 4 and 6, we have

Therefore, from the above proof, we can conclude that (20) and (21) are true. The proof of this lemma is complete.

Based on Lemmas 4–7, we are ready to solve the parameter problem for the Kalman-type filter (6)-(7), which makes the filtering error covariance in the sense of minimum mean square error principle at every sampling time instant.

Theorem 8. The filter gain matrix of the Kalman-type filter (6)-(7) is shown as follows:

and the covariance matrix of the filtering error is computed recursively by the following equation:

where

and the parameters , , , and are calculated by Lemmas 4–7, respectively.

Proof. Based on similar method of the projection theory, we have

Substituting (2) into , it has

Then,

where

Since the values of the parameters to are equal to zero, thus,

By using Lemmas 4 and 6, can be computed as follows:

From (29) and (30), one has

On the other hand, since the fact that and , the following equation can be established:

By (31) and (32), we can conclude that (22) and (24) are true.
Subsequently, the following derivations are shown to compute the covariance matrix of the filtering error . From the definition of and (25), we have

The proof of this theorem is complete.

Remark 9. In summary, the purpose is to obtain the optimal state estimator by the following five equations:

where the parameters can be computed by Lemmas 4–7. In our paper, due to considering the stochastic nonlinearity functions and finite-step correlated noises, we need to compute the values of the parameters , , and , , , and . In addition, because the phenomenon of missing measurements is taken into account, there are many cross terms in the parameter . Now, the algorithm steps of the optimal recursive filter are shown as follows.

Step 1. We can compute by Lemma 4. Then, substituting and into (13), we can obtain .

Step 2. Calculate by substituting and into (21). Then, substituting and into (20), we have .

Step 3. The parameter can be computed by substituting , , , , and into (24).

Step 4. Substituting , , and into (22), we get .

Step 5. The covariance matrix of the filtering error is computed by substituting , , and into (23).

Step 6. The optimal filter is calculated by substituting (6) and (22) into (7).

Remark 10. Note that, compared with the previous results, we consider the complex network phenomena (multiplicative noises, stochastic nonlinearity functions, correlated noises, and missing measurements) and design the locally optimal filter to solve the state estimation problem for complex network systems based on the projection theory. Due to the presence of complex network phenomena in the target system model and to further improve the accuracy of filtering estimation for the systems, we compute the parameters , , , and to deal with the finite-step correlated noises and compute the parameter to obtain the filter gain matrix and the covariance matrix . Hence, the proposed filter is suitable for online applications and we make the first attempt to propose the locally optimal filter for discrete-time stochastic systems with multiplicative noises, stochastic nonlinearity functions, correlated noises, and missing measurements. Also, the filtering algorithm accuracy can be improved because we have made great efforts to compensate the effects from multiplicative noises, stochastic nonlinearity functions, correlated noises, and missing measurements.

4. An Illustrative Example

In this section, a numerical example is proposed to show the feasibility and effectiveness of the developed main results and a target tracking system example is given to illustrate that the proposed filter is suitable for the target tracking systems object to the complex network phenomena (multiplicative noises, stochastic nonlinearity functions, correlated noises, and missing measurements).

Example 1. Consider the following systems with multiplicative noises and stochastic nonlinearity functions:

where

; and are uncorrelated zero-mean Gaussian white noises with unity covariances.

Here, we assumed that the noises are two-step correlated, that is, . Then, the process noises and the observation noises can expressed as follows:

where , , and and , , and ; is uncorrelated zero-mean Gaussian white noises with unity covariances. Then, the process noises satisfy the following equation:

Similarly, we can obtain the following results:

We chose the stochastic nonlinear functions and as follows:

where is the th element of the system state and , are zero-mean Gaussian white noises with unity covariances and uncorrelated with other signals. It is easy to see that the stochastic nonlinear functions satisfy the following equations:

where .

Let

and MSE denote the mean square error for the estimation of , that is, , where and is the number of simulation test.

According to Theorem 8, the locally optimal filter can be constructed by applying the projection theory and the minimum mean square error estimation principle. Based on the given Kalman filtering algorithm and Matlab software, the filter gain matrices and the error covariance matrices at every time step can be recursively computed and the values are shown as in Tables 1 and 2. The simulations are shown in Figures 1–6. Among them, Figures 1 and 2 plot the of the proposed filtering algorithm. The state of the actual system and the new designed estimation are plotted in Figures 3 and 4. Figures 5 and 6 plot the comparison between the filter in our paper and the filter without the multiplicative noises and missing measurements. From the simulations, it is easily seen that the proposed filtering algorithm can estimate the system state well irrespective of the correlated noise and the occurrence of the randomly missing measurements. When deriving the locally optimal filter, the additional efforts have been made to compensate the effects from stochastic nonlinearity functions, finite-step correlated noises, and missing measurements.

Table 1: Filter gain matrices .

Table 2: The error covariance matrices .

Figure 1: .

Figure 2: .

Figure 3: The trajectories of and .

Figure 4: The trajectories of and .

Figure 5: The trajectories of .

Figure 6: The trajectories of .

Example 2. Consider the following target tracking systems with multiplicative noises and stochastic nonlinearity functions:

where is the sampling time, . , is the location, and is the velocity at time instant . The other parameters are the same as those in Example 1.

The following simulation results in Figures 7 and 8 are obtained and we can see that the proposed filter, which can effectively estimate the position of the target, is suitable for the target tracking systems object to the complex network phenomena (multiplicative noises, stochastic nonlinearity functions, correlated noises, and missing measurements).

Figure 7: .

Figure 8: The trajectories of position.

5. Conclusion

The problem of locally optimal filtering has been investigated for a class of discrete-time systems with stochastic nonlinearity functions, finite-step correlated noises, and missing measurements. The multiplicative noises have been employed to describe the random disturbances in the system model. Due to the unreliable factors, the phenomena of the finite-step correlated noises and missing measurements have been considered in the target plant. Based on the projection theory, a class of Kalman-type locally optimal filter is designed and the filtering error covariance matrix is minimized in the sense of minimum mean square error principle. By solving the recursive matrix equation, we can obtain the filter gain. Finally, two examples are given to illustrate the feasibility and effectiveness of the proposed filtering scheme. One of the future research topics would be the stability analysis of the proposed algorithm.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China (NSFC) under Grant 51575110 and Fujian Young and Middle-Aged Teachers’ Science and Technology Research Project under Grant JAT160725.

References

R. E. Kalman, “A new approach to linear filtering and prediction problems,” Journal of Fluids Engineering, vol. 82, no. 1, pp. 35–45, 1960.

H. Hoely, K. J. H. Lawz, and R. Temponex, “Multilevel ensemble Kalman filtering,” SIAM Journal on Numerical Analysis, vol. 54, no. 3, pp. 1813–1839, 2016.

Z. Cao, R. Zhang, Y. Yang, J. Lu, and F. Gao, “Discrete-Time Robust Iterative Learning Kalman Filtering for Repetitive Processes,” IEEE Transactions on Automatic Control, vol. 61, no. 1, pp. 270–275, 2016.

V. Stojanovic and N. Nedic, “Robust Kalman filtering for nonlinear multivariable stochastic systems in the presence of non-Gaussian noise,” International Journal of Robust and Nonlinear Control, vol. 26, no. 3, pp. 445–460, 2016.

T. Soubdhan, J. Ndong, H. Ould-Baba, and M.-T. Do, “A robust forecasting framework based on the Kalman filtering approach with a twofold parameter tuning procedure: Application to solar and photovoltaic prediction,” Solar Energy, vol. 131, pp. 246–259, 2016.

S. Liu, G. Wei, Y. Song, and Y. Liu, “Extended Kalman filtering for stochastic nonlinear systems with randomly occurring cyber attacks,” Neurocomputing, 2016.

X. Gao, L. Mo, Z. Xiao, X. Chen, and S. Katayama, “Seam tracking based on kalman filtering of micro-gap weld using magneto-optical image,” The International Journal of Advanced Manufacturing Technology, vol. 83, no. 1-4, pp. 21–32, 2016.

W. Li, Y. Jia, and J. Du, “Tobit Kalman filter with time-correlated multiplicative measurement noise,” IET Control Theory & Applications, vol. 11, no. 1, pp. 122–128, 2017.

T. Tian, S. Sun, and N. Li, “Multi-sensor information fusion estimators for stochastic uncertain systems with correlated noises,” Information Fusion, vol. 27, pp. 126–137, 2016.

D. Chen, L. Xu, and J. Du, “Optimal filtering for systems with finite-step autocorrelated process noises, random one-step sensor delay and missing measurements,” Communications in Nonlinear Science and Numerical Simulation, vol. 32, pp. 211–224, 2016.

Z. Xing, Y. Xia, L. Yan, K. Lu, and Q. Gong, “Multisensor Distributed Weighted Kalman Filter Fusion With Network Delays, Stochastic Uncertainties, Autocorrelated, and Cross-Correlated Noises,” IEEE Transactions on Systems, Man, and Cybernetics: Systems, no. 99, pp. 1–11, 2016.

W. Liu, “Optimal Filtering for Discrete-Time Linear Systems with Time-Correlated Multiplicative Measurement Noises,” IEEE Transactions on Automatic Control, vol. 61, no. 7, pp. 1972–1978, 2016.

M. J. García-Ligero, A. Hermoso-Carazo, and J. Linares-Pérez, “Estimation from a multisensor environment for systems with multiple packet dropouts and correlated measurement noises,” Applied Mathematical Modelling, vol. 45, pp. 802–812, 2017.

H. Shu, S. Zhang, B. Shen, and Y. Liu, “Unknown input and state estimation for linear discrete-time systems with missing measurements and correlated noises,” International Journal of General Systems, vol. 45, no. 5, pp. 648–661, 2016.

S. Wang, H. Fang, and X. Tian, “Minimum variance estimation for linear uncertain systems with one-step correlated noises and incomplete measurements,” Digital Signal Processing, vol. 49, pp. 126–136, 2016.

R. Caballero-Águila, A. Hermoso-Carazo, and J. Linares-Pérez, “Networked fusion filtering from outputs with stochastic uncertainties and correlated random transmission delays,” Sensors, vol. 16, no. 6, 2016.

W.-Q. Liu, X.-M. Wang, and Z.-L. Deng, “Robust centralized and weighted measurement fusion Kalman estimators for uncertain multisensor systems with linearly correlated white noises,” Information Fusion, vol. 35, pp. 11–25, 2017.

N. Pletschen and K. J. Diepold, “Nonlinear state estimation for suspension control applications: a Takagi-Sugeno Kalman filtering approach,” Control Engineering Practice, vol. 61, pp. 292–306, 2017.

D. Ding, Y. Shen, Y. Song, and Y. Wang, “Recursive state estimation for discrete time-varying stochastic nonlinear systems with randomly occurring deception attacks,” International Journal of General Systems, vol. 45, no. 5, pp. 548–560, 2016.

Y. Sun, X. Wu, J. Cao, Z. Wei, and G. Sun, “Fractional extended Kalman filtering for non-linear fractional system with Lévy noises,” IET Control Theory & Applications, vol. 11, no. 3, pp. 349–358, 2017.

L. Li, D. Yu, Y. Xia, and H. Yang, “Stochastic stability of a modified unscented Kalman filter with stochastic nonlinearities and multiple fading measurements,” Journal of The Franklin Institute, vol. 354, no. 2, pp. 650–667, 2017.

J. Devakumar and K. Srinivasan, “Improved state estimation for stochastic nonlinear chemical reactor using particle filter based on unscented transformation,” Journal of Control Engineering and Applied Informatics, vol. 18, no. 4, pp. 650–667, 2016.

H. Yan, F. Qian, F. Yang, and H. Shi, “H∞ filtering for nonlinear networked systems with randomly occurring distributed delays, missing measurements and sensor saturation,” Information Sciences, vol. 370-371, pp. 772–782, 2016.

S. Zhang, J. Zhao, Y. Zhao, and G. Li, “Gain-Constrained Extended Kalman Filtering with Stochastic Nonlinearities and Randomly Occurring Measurement Delays,” Circuits, Systems and Signal Processing, vol. 35, no. 11, pp. 3957–3980, 2016.

Y. Liu, Z. Wang, X. He, and D. Hua, “H∞ filtering for non-linear systems with stochastic sensor saturations and Markov time delays: the asymptotic stability in probability,” IET Control Theory & Applications, vol. 10, no. 4, pp. 1706–1715, 2016.

J. Mao, D. Ding, Y. Song, Y. Liu, and F. E. Alsaadi, “Event-based recursive filtering for time-delayed stochastic nonlinear systems with missing measurements,” Signal Processing, vol. 134, pp. 158–165, 2017.