Suboptimal Filtering of Networked Discrete-Time Systems with Random Observation Losses

This paper studies the remote filtering problemover a packet-droppingnetwork.A generalmultiple-input-multiple-output (MIMO) discrete-time system is considered. The multiple measurements are sent over different communication channels every time step, and the packet loss phenomenon in every communication channel is described by an independent and identically distributed (i.i.d) Bernoulli process. A suboptimal filter is obtained which can minimize the mean squared estimation error. The convergence properties of the estimation error covariance are studied, andmean square stability of the suboptimal filter is proved under standard assumptions. A simulation example is exploited to demonstrate the effectiveness of the results.


Introduction
Filtering systems which transmit data packets through communication networks are called network-based filtering systems [1].The introduction of networks brings many attractive advantages, such as low cost, fast deployment, and flexible installation.However, communication networks are usually unreliable and may give rise to packet losses and networkinduced delays due to inherent limited bandwidth.The packet losses and network-induced delays can degrade the performance or even cause instability of the systems.Hence, it is not surprising that, in the past few years, the study of the state estimation problem for network-based filtering systems with packet losses and time-delays has been an active research area; see [2][3][4][5][6][7], to name a few.This paper is concerned with the design of the filter for networked discrete-time systems with random observation losses.
In the literature, there have been commonly two approaches for modeling the packet loss phenomenon in the network-based filtering systems.The first approach is to use a system with Markovian jumping parameter to represent random packet loss model [1].Note that such class of systems is a special class of Markovian jumping systems; hence, some results of control synthesis [8,9] and filtering methods [10] for Markovian jumping systems may be extended to these systems.The second approach is to view the packet loss as an independent Bernoulli process.And such approach has been used to deal with the estimation problems for network-based filtering systems with missing or intermittent observations [11][12][13].
As is well known, Kalman filtering [14] is one of the most popular and useful approaches to the filtering problem.In the literature, a few results have been reported on the Kalman filtering problem with observation losses [11][12][13][15][16][17].The early studies on the Kalman filtering with uncertain observation can be traced back to [15], where the linear minimum mean squared error (LMMSE) estimation method is considered.More recently, the LMMSE filter is obtained for systems with multiple packet dropouts in [16], where the number of possible consecutive packet dropouts is limited by a known bound.Besides, by the state augmentation, the LMMSE optimal filter, predictor, and smoother are designed for systems with finite consecutive packet dropouts [17].The LMMSE filtering only uses the statistics of the unobserved uncertainty sequence.In fact, the filter can get the information whether a packet has been delivered or not for networked filtering systems.In recent years, the paper [11] proposes a new filtering method which is called Kalman filtering with intermittent observations.The filter proposed in this paper exploits additional information regarding the packet arrival indicator sequence.As a consequence, the filter in [11] can give better performances.Nonetheless, the analysis of the filter is only limited to the boundary analysis due to the complex discussions [11,12].Motivated by the above analysis, the paper [13] proposes a new suboptimal estimator under a new performance index, which improves the performance of the LMMSE Kalman filter and possesses better properties of convergence and stability than the Kalman filter with intermittent.
Our paper extends the results in [13] to a more general case.In [13], the traditional assumption that all the measurements are encoded together and transmitted to the remote filter via a common communication channel is made.Thus the measurements are either received in full or lost completely.However, in the practical networked filtering systems, the measurements usually cannot be encapsulated into one data packet and multiple measurements must be transmitted through different communication channels.Moreover, the packet loss processes in different channels are often distinct.This is the motivation of the present paper.Our paper investigates the suboptimal filtering problem for the discretetime systems with multichannel transmission mechanism.The convergence of the estimation error covariance and the mean square stability of the filter are proved.It should be pointed out that the presented results can also be used for the systems with all measurements sent via one common communication channel.
This remainder of the paper is organized as follows.Section 2 formulates the problem and makes some preliminaries.The main results of this paper are presented in Section 3. The suboptimal filter is derived, and the convergence and stability of the suboptimal filter are proved under standard assumptions.A simulation example is given to demonstrate the effectiveness of the approach in Section 4. Section 5 concludes the paper.

Problem Statements and Preliminaries
Consider the following network-based filtering system: where () ∈   is the system state, () ∈   is the measurement, and () ∈   and V() ∈   are, respectively, the system noise and measurement noise with zero mean and covariance matrices {()  ()} =  , and {V()V  ()} =  , , in which  > 0,  > 0, and  , is the Kronecker delta function.The initial state (0) is also a random vector of mean  0 and covariance  0 . and  are constant matrices of appropriate dimensions.The scenario under consideration is illustrated in Figure 1, where the measurement () is partitioned into  parts; that is, The measurement components are transmitted to the remote filter through  different communication channels.The independent and identically distributed (i.i.  () ( ̸ = ) are independent for every  and ; then we have

Remote filter Plant
From the above assumption and analysis, the multiplicative noise matrix () can be expressed by the following diagonal binary random matrix with entries of either 1 or 0 in the diagonal: It is noted that, as shown in Figure 1,  measurements are transmitted to the filter via  channels.In fact, some measurements may be encoded together and sent over the network in a single packet; that is to say, the number of communication channels can be smaller than that of the measurements.The model proposed in this paper can be easily adjusted to describe the above case.For example, assume the number of communication channels is 4. The first two measurements are encoded together and transmitted by a common channel, and the last two measurements are transmitted by another common channel.Thus () can be written as Throughout this paper, without loss of generality, the following assumptions are made for technical convenience.
Assumption 1.The random processes (), V(), () for all  and the initial state (0) are mutually independent.Assumption 2. The remote filter can obtain the information regarding () by employing the time-stamp technique.
Similar to [13], we introduce the following innovation sequence {()}: where and () is the gain matrix of the filter, which can be chosen such that the following is minimized: According to [13, Lemma 1], it is easy to obtain that () is mutually uncorrelated noise with zero mean.The filtering problem considered in this paper is to find the state estimation x( + 1 | ) in (6), in which () is chosen such that ( 7) is minimized.Remark 3. Our paper extends the results in [13] where all the measurements are transmitted to the filter via a common communication channel.Assume that the packet loss processes in all the channels are identical; the filter proposed in this paper is equivalent to that in [13], which means that the suboptimal filter in [13] can be regarded as a special case of our filter.
Remark 4. Our estimation problem is different from [12] where the expectation is only taken over the system noise () and the measurement noise V().In this paper, the expectation is taken over not only () and V() but also the multiplicative noise matrix () with entries   () in the diagonal.Moreover, it is assumed that the measurements are sent to the filter via two communication channels for the sake of simplicity in [12], and while we address the filtering problem in the more general case, the observation processes can be sent by multiple communication channels.
Remark 5.The filter proposed in this paper is also different from the LMMSE optimal Kalman filter with multiplicative noise [18] where the innovation sequence {  ()} is defined as in which  is the mean of the random matrix () and  = diag( 1 ,  2 , . . .,   ).And the state estimation in [18] can be written as The aim of [18] is to find   () such that (7) is minimized.
Obviously, this filter only uses the statistics of the multiplicative noise matrix (), while the filter proposed in this paper exploits additional information regarding the packet arrival indicator sequence.Hence, our filter may give better performances.

Suboptimal Filter Design.
The following theorem gives the recursive equations of the suboptimal filter defined in (6).
Theorem 6.Consider the systems (1) and (2); the suboptimal filter defined in (6) is given as follows: where the gain matrix () is calculated as and the state estimation error covariance () is given by in which () is defined as in (15) below. where Then we can rewrite (14) as where If we choose () =  * (), (18) can be minimized.This completes the proof.

Convergence and Stability of the Suboptimal Filter.
In this subsection, the convergence and stability of the proposed filter are studied.The following preliminary lemmas are introduced before presenting the main objectives of this paper.
Lemma 7. Consider the following operators: where Assume  is symmetric and positive semidefinite; then the following facts are true.
(ii) The proof for fact (ii) is somewhat more technical since   in the operator   (, ) is the implicit expression of the vector .In order to facilitate the proof of fact (ii), we will transform (21) to the explicit expression of .First, we can rewrite (21) as follows: where Then, let us define It is easy to obtain that Note ( 23) and then we claim that (28) can be described by where   ( = 1, 2, . . ., ) are the diagonal matrices and they have all diagonal elements being zeros except for the th element being √  (1 −   ).The above analysis leads to the following expression of (24): It is clear that (30) is the explicit expression of vectors  and .We are now ready to derive fact (ii).
(i) For all  ≥ 0, lim  → ∞    () = 0. (ii) Define the linear system where  ≥ 0.Then, the sequence () is bounded.The proof can be derived similarly to [11,Lemma 3]; therefore, it is omitted.Theorem 9. Suppose there exists a matrix K and a positive definite matrix P > 0 such that P >   ( K, P); then for any initial condition  0 , (12) converges to a unique positive semidefinite matrix .That is to say, lim  → ∞ () =  ≥ 0.  Note that  = + K K > 0; we can easily obtain that   ( K, P) > ψ ( P).Then, in view of the assumption P >   ( K, P), we have P > ψ ( P).Therefore, ψ meets the condition of Lemma 8.
In the following, we show that () initialized at (0) =  0 ≥  also converges to the same limit .
It follows from Lemma 7 (iii) that Since   meets the condition of Lemma 8, from the above analysis, we have Thus we claim that () initialized at (0) =  0 ≥  also converges to the same limit .
We now establish that the sequence () converges to  for all initial conditions.Note that Since  0 () and   0 () converge to the same limit , lim  → ∞   0 () =  holds.This completes the proof.
Taking limits in (10), we can rewrite the estimator (10) as Next, we will present the result that (52) is mean square stable.
Theorem 10.Suppose there exists a matrix K and a positive definite matrix P > 0 such that P >   ( K, P); the filter (52) is mean square stable if the limit  of (12) is exactly positive definite.
Proof.First it is obvious that the mean square stability of the following system is equivalent to that of the filter (52): x ( + 1 | ) = ( −  () )  x ( |  − 1) .
From Lyapunov inequality in [19], we conclude that if we can find a positive definite matrix  1 > 0 satisfying From (45) and the assumption  > 0, it is obvious that there exists  1 =  such that (55) holds.Hence the proof is complete.

Figure 1 :
Figure 1: Diagram of networked filtering system with random packet losses.

Figure 2 :Figure 3 :
Figure 2: Tracking performance of the proposed filter.(a) The first component of ().(b) The second component of ().