Improved Minimum Entropy Filtering for Continuous Nonlinear Non-Gaussian Systems Using a Generalized Density Evolution Equation

This paper investigates the filtering problem for multivariate continuous nonlinear non-Gaussian systems based on an improved minimum error entropy (MEE) criterion. The system is described by a set of nonlinear continuous equations with non-Gaussian system noises and measurement noises. The recently developed generalized density evolution equation is utilized to formulate the joint probability density function (PDF) of the estimation errors. Combining the entropy of the estimation error with the mean squared error, a novel performance index is constructed to ensure the estimation error not only has small uncertainty but also approaches to zero. According to the conjugate gradient method, the optimal filter gain matrix is then obtained by minimizing the improved minimum error entropy criterion. In addition, the condition is proposed to guarantee that the estimation error dynamics is exponentially bounded in the mean square sense. Finally, the comparative simulation results are presented to show that the proposed MEE filter is superior to nonlinear unscented Kalman filter (UKF).


Introduction
State estimation theory has been regarded as an important research area in modern control systems.In particularly, the appearance of the Kalman filtering theory in the last century had a profound influence on modern optimal control [1][2][3].Kalman filter algorithm is based on the minimum variance estimation for the linear Gaussian systems.However, the noises in real systems may be non-Gaussian; and even for a nonlinear system with Gaussian noises, the system output could be a non-Gaussian variable.In these cases, Kalman filter may lead to poor estimation.Some attempts have been made at studying filtering algorithms for nonlinear systems with non-Gaussian noises.The existed methodologies to design filters for nonlinear non-Gaussian systems can be classified into three kinds: analytical approach, simulation-based approach and adaptive approach.
Analytical approaches to filtering for nonlinear non-Gaussian systems were investigated in [4][5][6][7].In [4,5], nonlinear non-Gaussian systems were described by combining an improved square-root B-spline model with a further nonlinear dynamic model.Once B-spline expansions have been made for probability density functions (PDFs), further modeling was still needed to reveal the relationship between the input and the weights related to the PDFs, a nonlinear filter was then constructed by minimizing the error between the measured output PDF and estimated output PDF.Although the proposed filter is suitable to nonlinear non-Gaussian systems, the output PDF should be measureable.Moreover, it is not easy to build the state space expression of the weights related to the PDFs in practical systems.The filter for a class of multivariate dynamic stochastic systems with non-Gaussian stochastic input and nonlinear output was investigated in [6] and [7] respectively.In [6], a new formulation of the residual PDF was made to link the residual PDF to the gain matrix of the filter, and the optimal filtering gain matrix was then solved by minimizing the entropy of the residual.The minimum entropy filter in [6] presented a good performance in reducing the randomness of the filter residual, and was more general and suited for non-Gaussian systems.However, minimum entropy criterion may not guarantee that the estimation errors approach to zero.In [7], following the minimum information divergence criterion, a hybrid characteristic function of the conditional estimation error was introduced to construct the performance index of the tracking filter.An analytical solution of the filter gain matrix was then obtained so that the PDFs of the filtering error can follow a target distribution shape.Nevertheless, it is a little complicated to calculate the analytical solution.
There are two types of simulation-based approaches to filtering for nonlinear non-Gaussian systems: numerical integration and sequential Monte Carlo.Although the filtering problem can be tackled using numerical integration [8,9], it is difficult to implement when the dimension of the state vector is higher.Sequential Monte Carlo simulation [10][11][12][13][14][15], which is also named particle filtering strategy, has shown its great advantages to deal with filtering problems for nonlinear non-Gaussian systems.Nevertheless, there are still many issues to be solved, for example, (1) random sampling may bring the accumulation of Monte Carlo error and even lead to filter divergence; (2) a large number of particles are needed to avoid degradation and to improve the estimation accuracy, which makes the calculation a sharp increase.
Adaptive approaches have been also investigated for solving filtering problem in nonlinear non-Gaussian systems in last decades.Since minimum error entropy criterion ensures that the estimation error has small uncertainty, it was used for supervised training of nonlinear stochastic system in [16,17].However, entropy does not change with the mean of the distribution, the algorithm may not yield zero-mean error.Therefore, the result may be corrected by properly modifying the bias of the output processing element of the neural networks.Maximum mutual information criterion was proposed for adaptive filtering in [18], this criterion is robust to measure distortions.Nevertheless, the maximum mutual information criterion leads to non-unique optimum solution.It is necessary to use a priori information about the unknown system in order to obtain unique solution.
In general, the filtering problem for nonlinear non-Gaussian systems calls for further investigation.The filtering problem addressed here is solved by combining the improved performance index with optimal design method.
Entropy is a natural extension beyond mean square error because it is a function of PDF.One of the most important problems for minimum entropy filtering is the formulation of the PDF of estimation error.For continuous nonlinear systems, the classical Liouville equation and its Dostupov-Pugachev extension have been used to obtain the PDF of the concerned stochastic variable.However, it is very difficult to have their analytical/numerical solutions due to high-dimensional partial differential equations.From the viewpoint of the random event description of the principle of preservation of probability, the generalized density evolution equation was developed in [19].Different from the traditional density evolution equations, the dimensions of the generalized density evolution equation just depend on the necessary physical quantities' dimensions not the whole original systems itself, which makes the dimensions of the partial differential equations much lower.Some application results (see e.g., [20,21]) illustrated the efficiency and conveniences of the formulated density evolution equation.
The contribution of this paper is to develop a new filtering strategy for multivariable nonlinear systems with non-Gaussian disturbances by utilizing a novel performance index which contains the entropy of estimation error, square error and constraints on gain matrix of the filter.A novel approach, which uses the principle of preservation of probability, is presented here to formulate the joint PDF evolution equation of the estimation errors.The PDF evolution equation explicitly reveals the relationship among the estimation errors of the filter, filter gain and random inputs.In addition, the entropy of estimation errors in the performance index is replaced by its information potential presented in [16,17] so as to simplify the calculation of the entropy.The filter gain matrix for nonlinear non-Gaussian systems is then designed by minimizing the proposed performance index.This filtering algorithm yields to the estimation errors which not only have small uncertainty but also approach to zero.Finally, the exponentially boundedness in the mean square sense is analyzed for the estimation error dynamics.
This paper is organized as follows: in section 2, the state-space model of a nonlinear non-Gaussian system and the structure of filter are built to formulate the filtering problem.The joint PDF evolution of estimation errors is formulated in Section 3 so as to calculate the entropy of estimation errors.The filtering algorithm is proposed in Section 4 by minimizing the improved entropy criterion and the exponentially mean-square boundedness condition for error systems is provided.Comparative simulation results are given to illustrate the efficiency and validity of the proposed method in Section 5. Conclusions are drawn in Section 6. Notation.n  and n m   denote the n-dimensional Euclidean space and the set of all n m  real matrices respectively.The superscript "T " denotes the transpose.
  min   represent the largest and smallest eigenvalue of A respectively.   Ε stands for the mathematical expectation of random variables.

System Model
Consider a nonlinear stochastic continuous system described by: where x   is the state vector, m y  is the measured output vector.
p    and q    are system noises and measurement noises respectively.( ) A t and ( ) D t are two known time-varying system matrices.The following assumptions are made for simplicity, which can be satisfied by many practical systems.
Assumption 1: ( ) F  is a known Borel measurable and smooth vector-value nonlinear function of its arguments.Assumption 2:  and  are bounded, mutually independent random vectors with known joint PDFs   and   .

Filter Dynamics
For the nonlinear stochastic system (1), the filter dynamics can be formulated as follows: The resulting estimation error ê x x   satisfies: (3) where is the gain matrix to be determined and can be denoted as is a stretched column vector.The purpose of filtering algorithm for nonlinear non-Gaussian systems is to ensure that the estimation errors achieve small dispersion and approach to zero.Entropy provides a measure of the average information contained in a given PDF.When entropy of estimation errors is minimized, all moments of the error PDF are constrained.Hence, the following quadratic Renyi's entropy of estimation error is still contained in the performance index to measure the dispersion of the estimation error: (4) where ( ) e p z is the joint PDF of estimation error and ( ) V e the information potential [14].It can be seen that it is necessary to provide the joint PDF of errors firstly.In the next section, we will deduce the expression of ( ) e p z according to the generalized density evolution equation.

Formulation for the Joint PDF of Error
In most cases, Equation ( 3) is a well-posed equation, and the error vector ( ) e t can be determined uniquely.It may be assumed to be: ( At the present stage, the explicit expression of ( ) H  is not requisite, and the sufficient condition just needs to know its existence and uniqueness.The derivative of ( ) e t can be assumed to take the form: It is observed that all the randomness involved in this error dynamics comes from noises  and  , thus, the augmented system ( ( ), , ) e t   is probability preserved.In other words, from the random event description of the principle of preservation of probability, it leads to: (7) where e  ,   and   are the distribution domains of e ,  and  , respectively; , , ( , , , )   e p e t Combining Equation (7) with Equation ( 8) and considering the arbitrariness of e      , it yields: (9) where Equation ( 6) and Equation ( 9) are the generalized density evolution equation (GDEE) for ( ( ), , ) e t   .The corresponding instantaneous PDF of ( ) e t can be obtained by solving a family of partial differential equations with the following given initial condition: (10) where ( )   is the Dirac-Delta function; and 0 0 0 ê x x   is deterministic initial value of ( ) e t .Then, we have: (11) where the joint PDF ( , , , ) is the solution of (8), which can be obtained according to the method presented in [22].

Performance Index
The following performance index for a minimum entropy filter was presented in [6]: where 1 R and 2 R are weighing matrices.The first term is the entropy of the estimation errors and the second term is used to penalize the elements of the gain.Although the MEE criteria can minimize both the probabilistic uncertainty and dispersion of the estimation error, it cannot guarantee the minimum error.In this work , an improved performance index J is considered: (13) where 2 ( ( )) H e t is the quadratic Renyi's entropy given in Equation ( 4). 3 0 R  is the weight corresponding to the mean squared error.The third term on the right side of the equation is utilized to make the estimation errors approach to zero.Since Renyi's entropy is a monotonic increasing function of the negative information potential, minimization of Renyi's entropy is equivalent to minimizing the inverse of the quadratic information potential 1 , so the performance index can be rewritten as follows: (14) Remark 1: The entropy of estimation error is replaced by quadratic information potential of error in order to simplify the calculation.Nonparametric estimation of the PDF of estimation errors is another alternative way to estimate the PDF except the method introduced in [22].Nonparametric entropy (or information potential) estimation of random error can be found in [16,17].

Optimal Filter Gain Matrix
After constructing the performance index J , the optimal filter gain can be carried out by minimizing the performance index Equation ( 14), which can be regarded as an unconstrained nonlinear programming problem.The classical conjugate gradient approach is employed here to obtain an iterative solution, which could avoid the difficulty of solving the extremum condition ( 0 ).The elements of optimal filter gain can be solved and summarized in Theorem 1.
Theorem 1: For a given accuracy 0   , the gain matrix of the suboptimal improved MEE filter is given by: ( where ( ) ( ) and the Hessian matrix ( ) Proof: Choose an initial value 0 l , and denote : 0 k  .Search the next value 1 k l  from the given point k l along the conjugate gradient direction k P .Repeat the above procedure until reach the given accuracy  .
Consequently, we can obtain the optimal filtering algorithm shown in Equation (15).Based on the conjugate gradient method, the steps to solve the optimal filter gain are summarized as follows: Step (1) Give an accuracy 0   and initialize 0 l , denote k = 0 Step (2) Set 0 0 ( ) P J l   .

Remark 2:
The conjugate gradient method here could overcome the slow convergence of the steepest descent approach.It is one of the most effective algorithms to deal with large-scale nonlinear optimization problems.

Exponentially Bounded in the Mean Square
In this section, the Lyapunov method is utilized to give the condition under which the filter error dynamics is exponentially ultimately bounded in the mean square.Firstly, in order to simplify the problem, we give the following assumption:

Assumption 3: The nonlinear vector function ( , )
F   is assumed to satisfy (0, 0) 0 F  [23]and: (17) where For estimation error dynamics (3), the exponentially mean-square boundedness is defined as follows: Definition 1 [24]: For all initial conditions 0 n e  , the dynamics of the estimation error (i.e., the solution of the system (3)) is exponentially ultimately bounded in the mean square if there exist constants 0 Moreover, filter ( 2) is said to be exponential if, for every 0 n e  , system (3) is exponentially ultimately bounded in the mean square.
Theorem 2: Let the filter parameter L be given.If there exist positive scalars 1  , 2  and 3  such that the following matrix inequality: has a positive definite solution 0 P  , then system (3) is exponentially ultimately bounded in the mean square.
where 0 P  is the solution to the matrix inequality Equation (19).For notational convenience, we give the following definitions: (21) and then it follows from Equation (3) that: (22) Hence, the derivative of V along a given trajectory Equation ( 22) is obtained as:  and 3  be positive scalars.Then the following matrix inequalities hold: Noticing Assumptions 3 and Equation ( 17), we have: Substituting Equations ( 24)-( 27) into Equation ( 23) results in:  .Since 0 T  is arbitrary, the definition of exponential ultimate boundedness in ( 18) is then satisfied, and this completes the proof of Theorem 2.

Simulation Results
In order to show the applicability of the proposed filtering algorithm, we consider the following nonlinear system represented as : R  , 2 2 R  and 3 10 R  , respectively.The simulation results based on the MEE filter are shown in Figures 2-6, 7a and 8a).And the comparative results between MEE filter and UKF are shown in Figures 7 and 8. Figure 2 demonstrates that the performance index decreases monotonically with time.In Figures 3 and 4, both the range and PDFs of errors are given, it can be seen that the shapes of PDFs of the tracking errors become narrower and sharper along with the increasing time, which illustrates the dispersions of estimation errors can be reduced.In order to clarify the improvements, the initial and final PDFs are shown in Figures 5 and 6.It can be shown from Figures 3-6 that the proposed improved MEE filter can decrease the uncertainties of the estimation errors and drive the estimation errors approaching to zero. Figure 7 shows that both the MEE filter and UKF make the estimation errors approach to zero.But the estimation errors in Figure7a are more closer to zero with smaller randomness.The variances of the filter gains in two cases are presented in Figure 8.It can be seen that the gains of UKF are not convergent.
From the above analyses, it is obvious that the proposed strategy has better performance than UKF.Therefore, the proposed MEE filter is more suitable for nonlinear stochastic systems with non-Gaussian noises.

Conclusions
In this paper, based on the improved MEE criterion, optimal filter design problem is studied for multivariate nonlinear continuous stochastic systems with non-Gaussian noises.A novel performance index, including Renyi's entropy of the estimation error, mean value of the squared error and constrains on filter gain, is formulated.By using the generalized density evolution equation, the relationship among the PDFs of estimation error, random disturbances and filter gain matrix is obtained.A recursive optimal filtering algorithm is obtained by minimizing the improved MEE criterion.The exponentially mean-square boundedness condition of the estimation error systems is established.In this work, the proposed MEE filter and UKF are both applied to a continuous nonlinear and non-Gaussian stochastic system.Comparative simulation results demonstrate the superiority of the presented control algorithm.
both sides from 0 to 0 T  and taking the expectation lead to: random disturbances  and  obey non-Gaussian, and their distributions are shown in Figure1.The weights in Equation (14) are