IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 63, NO. 4, APRIL 2018 1003 Kalman Filtering Over Gilbert–Elliott Channels: Stability Conditions and Critical Curve Junfeng Wu , Guodong Shi , Brian D. O. Anderson , Life Fellow, IEEE, and Karl Henrik Johansson , Fellow, IEEE
Abstract—This paper investigates the stability of Kalman tal question has inspired various significant results focusing on filtering over Gilbert–Elliott channels where random packet the interface of control and communication and has become drops follow a time-homogeneous two-state Markov chain a central theme in the study of networked sensor and control whose state transition is determined by a pair of failure and recovery rates. First of all, we establish a relaxed condition systems [2]Ð[4]. guaranteeing peak-covariance stability described by an in- Early works on networked control systems assumed that sen- equality in terms of the spectral radius of the system matrix sors, controllers, actuators, and estimators communicate with and transition probabilities of the Markov chain. We further each other over a finite-capacity digital channel, e.g., [2] and show that the condition can be interpreted using a linear [5]Ð[13], with the majority of contributions focused on one or matrix inequality feasibility problem. Next, we prove that the peak-covariance stability implies mean-square stability, both finding the minimum channel capacity or data rate needed if the system matrix has no defective eigenvalues on the for stabilizing the closed-loop system, and constructing optimal unit circle. This connection between the two stability no- encoderÐdecoder pairs to improve system performance. At the tions holds for any random packet drop process. We prove same time, motivated by the fact that packets are the funda- that there exists a critical curve in the failure-recovery rate mental information carrier in most modern data networks [3], plane, below which the Kalman filter is mean-square stable and no longer mean-square stable above. Finally, a lower many results on control or filtering with random packet dropouts bound for this critical failure rate is obtained making use appeared. of the relationship we establish between the two stability State estimation, based on collecting measurements of the criteria, based on an approximate relaxation of the system system output from sensors deployed in the field, is embedded in matrix. many networked control applications and is often implemented Index Terms—Estimation, Kalman filtering, Markov pro- recursively using a Kalman filter [14], [15]. Clearly, channel cesses, stability, stochastic systems. randomness leads to that the characterization of performance is not straightforward. A burst of interest in the problem of the stability of Kalman filtering with intermittent measurements has I. INTRODUCTION arisen after the pioneering work [16], where Sinopoli et al. mod- A. Background and Related Works eled the statistics of intermittent observations by an independent IRELESS communications are being widely used nowa- and identically distributed (i.i.d.) Bernoulli random process and W days in sensor networks and networked control systems. studied how packet losses affect the state estimation. Tremen- New challenges accompany the considerable advantages wire- dous research has been devoted to stability analysis of Kalman less communications offer in these applications, one of which is filtering or the closed-loop control systems over i.i.d. packet how channel fading and congestion influence the performance lossy packet networks in [17]Ð[22]. of estimation and control. In the past decade, this fundamen- To capture the temporal correlation of realistic commu- nication channels, the Markovian packet loss model has Manuscript received May 14, 2017; accepted July 14, 2017. Date of been introduced to partially address this problem. Since the publication July 27, 2017; date of current version March 27, 2018. This GilbertÐElliott channel model [23], [24], a classical two-state work was supported in part by the Knut and Alice Wallenberg Foundation, time-homogeneous Markov channel model, has been widely in part by the Swedish Research Council, and in part by the National Natural Science Foundation of China under Grant 61120106011. This applied to represent wireless channels and networks in indus- paper was presented in part at the 34th Chinese Control Conference, trial applications [25]Ð[27], the problem of networked control Hangzhou, China, July 2015. Recommended by Associate Editor S. over GilbertÐElliott channels has drawn considerable attention. Andersson. J. Wu and K. H. Johansson are with the ACCESS Linnaeus Center, Huang and Dey [28], [29] considered the stability of Kalman School of Electrical Engineering, Royal Institute of Technology, Stock- filtering with Markovian packet losses. To aid the analysis, holm SE-100 44, Sweden (e-mail: [email protected]; [email protected]). they introduced the notation of peak covariance, defined by G. Shi and B. D. O. Anderson are with the Research School of Engineering, The Australian National University, Canberra, ACT 0200, the expected prediction error covariance at the time instances Australia (e-mail: [email protected]; brian.anderson@anu. when the channel just recovers from failed transmissions, as edu.au). an evaluation of estimation performance deterioration. The Color versions of one or more of the figures in this paper are available online at http://ieeexplore.ieee.org. peak-covariance stability (PCS) can be studied by lifting the Digital Object Identifier 10.1109/TAC.2017.2732821 original systems at stopping times. Sufficient conditions for the
0018-9286 © 2017 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications standards/publications/rights/index.html for more information. 1004 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 63, NO. 4, APRIL 2018
PCS were proposed for general vector systems, and a necessary be recast as a linear matrix inequality (LMI) feasibility and sufficient condition for scalar systems. The existing problem. These conditions are theoretically and numer- literature [28]Ð[31] has however made restrictive assumptions ically shown to be less conservative than those in the on the plant dynamics and the communication channel. We literature. show by numerical examples that existing conditions for PCS 2) We prove that PCS implies MSS if the system matrix has only apply to relative reliable channels with low failure rate. no defective eigenvalues on the unit circle. Remarkably Moreover, existing results rely on calculating an infinite sum enough this implication holds for any random packet drop of matrix norms for checking stability conditions. process that allows PCS to be defined. This result bridges The PCS describes covariance stability at random stopping two stability criteria in the literature and offers a tool for times, whereas the mean-square stability (MSS) describes co- studying MSS of the Kalman filter through its PCS. Note variance stability at deterministic times. In the literature, the that MSS was previously studied using quite different main focus is on MSS rather than PCS. This is partly because methods such as analyzing the boundness of the expec- the definition of the former is more natural and practically useful tation of a kind of randomized observability Gramians than that of the latter. However, it is difficult to analyze MSS over a stationary random packet loss process to establish over GilbertÐElliott channels directly through the random Ric- the equivalence between stability in stopping times and cati equation, due to temporal correlation between the packet stability in sampling times [32], and characterizing the losses. If we could build clearer connections relating these two decay rate of the prediction covariance’s tail distribution stability notions, the MSS can be conveniently studied through for so-called nondegenerate systems [34]. PCS. The relationship between the MSS and the PCS was pre- 3) We prove that there is a critical p − q curve, with p be- liminarily discussed [29]. Improvements to these results can be ing the failure rate and q being the recovery rate of the found in [30] and [31]. Particularly in [32], by investigating the GilbertÐElliott channel, below which the expected pre- estimation error covariance matrices at each packet reception, diction error covariance matrices are uniformly bounded necessary and sufficient conditions for the MSS were derived for and unbounded above. The existence of a critical curve second-order systems and some certain classes of higher order for Markovian packet losses is an extension of that of systems. Although it was proved that with i.i.d. packet losses the critical packet loss rate subject to i.i.d. packet losses the PCS is equivalent to the MSS for scalar systems and sys- in [16]. However, the proof method of [16] does not ap- tems that are one-step observable [29], [31], for vector systems ply to Markovian packet losses. In this paper, the critical with more general packet drop processes, this relationship is curve is proved via a novel coupling argument, and to the unclear. best of our knowledge, this is the first time phase transi- There are some other works studying distribution of error co- tion is established for Kalman filtering over Markovian variance matrices. Essentially, the probabilistic characteristics channels. Finally, we present a lower bound for the crit- of the prediction error covariance are fully captured by its prob- ical failure rate, making use of the relationship between ability distribution function. Motivated by this, Shi et al. [33] the two stability criteria we established. This lower bound studied Kalman filtering with random packet losses from a prob- holds without relying on the restriction that the system abilistic perspective where the performance metric was defined matrix has no defective eigenvalues on the unit circle. In using the error covariance matrix distribution function, instead other words, we obtain an MSS condition for general lin- of the mean. Mo and Sinopoli [34] studied the decay rate of the ear time-invariant (LTI) systems under Markovian packet estimation error covariance matrix and derived the critical ar- drops. rival probability for nondegenerate systems based on the decay We believe these results add to the fundamental understanding rate. Weak convergence of Kalman filtering with packet losses, of Kalman filtering under random packet drops. i.e., that error covariance matrix converges to a limit distribu- The remainder of this paper is organized as follows. Section II tion, were investigated in [35]Ð[37] for i.i.d., semi-Markov, and presents the problem setup. Section III focuses on the PCS. Markov drop models, respectively. Section IV studies the relationship between the peak-covariance and mean-square stabilities, the critical p − q curve, and presents a sufficient condition for MSS of general LTI systems. B. Contributions and Paper Organization Section V demonstrates the effectiveness of our approach com- In this paper, we focus on the PCS and MSS of Kalman filter- pared with the literature using two numerical examples. Finally, ing with Markovian packet losses. We first derive relaxed and Section VI concludes this paper. N Sn explicit PCS conditions. Then, we establish a result indicating Notations: is the set of positive integers. + is the set of that PCS implies MSS under quite general settings. We even- n by n positive semidefinite matrices over the complex field. tually make use of these results to obtain MSS criteria. The For a matrix X, σ(X) denotes the spectrum of X and λX contributions of this paper are summarized as follows. denotes the eigenvalue of X that has the largest magnitude. 1) A relaxed condition guaranteeing PCS is obtained de- X∗, X, and X are the Hermitian conjugate, transpose, and scribed by an inequality in terms of the spectral ra- complex conjugate of X, respectively. Moreover, · means dius of the system matrix and transition probabilities of the 2-norm of a vector or the induced 2-norm of a matrix. ⊗ is the the Markov chain, rather than an infinite sum of matrix Kronecker product of two matrices. The indicator function of a norms as in [28]Ð[31]. We show that the condition can subset A⊂Ω is a function 1A :Ω→{0, 1}, where 1A(ω)=1 WU et al.: KALMAN FILTERING OVER GILBERT–ELLIOTT CHANNELS: STABILITY CONDITIONS AND CRITICAL CURVE 1005
The estimator computes xˆk|k , the minimum mean-squared error estimate, and xˆk+1|k , the one-step prediction, accord- ing to xˆk|k = E[xk |Fk ] and xˆk+1|k = E[xk+1|Fk ].LetPk|k and Pk+1|k be the corresponding estimation and prediction er- Fig. 1. Estimation over an erasure channel. ror covariance matrices, i.e., Pk|k = E[(xk − xˆk|k )(·) |Fk ] and Pk+1|k = E[(xk+1 − xˆk+1|k )(·) |Fk ]. They can be computed if ω ∈A, otherwise 1A(ω)=0. For random variables, σ(·) is recursively via a modified Kalman filter [16]. The recursions for the σ-algebra generated by the variables. xˆk|k and xˆk+1|k are omitted here. To study the Kalman filtering system’s stability, we focus on the prediction error covariance matrix P | , which is recursively computed as II. KALMAN FILTERING WITH MARKOVIAN PACKET LOSSES k+1 k Consider an LTI system Pk+1|k = APk|k−1 A + Q −1 xk+1 = Axk + wk (1a) −γk APk|k−1 C (CPk|k−1 C + R) CPk|k−1 A .
yk = Cxk + vk (1b) It can be seen that Pk+1|k inherits the randomness of {γt }1≤t≤k . n×n m ×n where A ∈ R is the system matrix and C ∈ R is We focus on characterizing the impact of {γk }k∈N on Pk+1|k . ∈ Rn the observation matrix, xk is the process state vec- To simplify notations in the sequel, let Pk+1 Pk+1|k , and ∈ Rm ∈ Rn k k Sn → Sn tor and yk is the observation vector, wk and define the functions h, g, h , and g : + + as follows m vk ∈ R are zero-mean Gaussian random vectors with auto- E ≥ E covariance [wk wj ]=δkjQ (Q 0), [vk vj ]=δkjR (R> h(X) AXA + Q (3) 0), E[wk vj ]=0∀j, k. Here, δkj is the Kronecker delta func- g(X) AXA + Q − AXC(CXC + R)−1 CXA tion with δkj =1if k = j and δkj =0otherwise. The initial state x is a zero-mean Gaussian random vector that is uncor- k ◦ ◦···◦ k ◦ ◦···◦ 0 h (X) h h h (X) and g (X) g g g (X) related with wk and vk and has covariance Σ0 ≥ 0. We assume ktimes that (C, A) is detectable and (A, Q1/2 ) is stabilizable. By apply- ktimes ing a similarity transformation, the unstable and stable modes (4) of the considered LTI system can be decoupled. An open-loop prediction of the stable mode always has a bounded estimation where ◦ denotes the function composition. The following sta- error covariance, therefore, this mode does not play any key role bility notion is standard. in the stability issues considered in this paper. Without loss of Definition 1: The Kalman filtering system with packet losses E ∞ generality, we assume that (A1) All the eigenvalues of A have is mean-square stable if supk∈N Pk < . magnitudes no less than 1. Certainly A is nonsingular. Define We consider an estimation scheme illustrated in Fig. 1, where τ min{k : k ∈ N,γ =0} the raw measurements of the sensor {yk }k∈N are transmitted to 1 k the estimator via an erasure communication channel over which β1 min{k : k>τ1 ,γk =1} packets may be dropped randomly. Denote by γ ∈{0, 1} the k . arrival of yk at time k: yk arrives errorfree at the estimator if . γk =1; otherwise γk =0. Whether γk takes value 0 or 1 is τ min{k : k>β− ,γ =0} assumed to be known by the receiver at time k. Define Fk as j j 1 k the filtration generated by all the measurements received by the βj min{k : k>τj ,γk =1}. (5) estimator up to time k, i.e., Fk σ (γt yt ,γt ;1≤ t ≤ k) and ∞ F = σ (∪ Fk ). We will use a triple (Ω, F, P) to denote the k=1 It is straightforward to verify that {τj }j∈N and {βj }j∈N are probability space capturing all the randomness in the model. two sequences of stopping times because both {τj ≤ k} and To describe the temporal correlation of realistic communica- {βj ≤ k} are Fk −measurable; see [38] for details. Due to the tion channels, we assume the GilbertÐElliott channel [23], [24], strong Markov property and the ergodic property of the Markov where the packet loss process is a time-homogeneous two-state chain defined by (2), the sequences {τ } ∈N and {β } ∈N have { } j j j j Markov chain. To be precise, γk k∈N is the state of the Markov finite values P-almost surely. Then, we can define the so- chain with initial condition, without loss of generality, γ1 =1. ∗ journ times at the state 1 and state 0, respectively, by τj and The transition probability matrix for the GilbertÐElliott channel β∗ ∀j ∈ N as is given by j ∗ − ∗ − 1 − qq τj τj βj−1 and βj βj τj P = (2) p 1 − p where β0 =1by convention. A result given by [29, Lemma 2] P | { ∗} { ∗} where p (γk+1 =0γk =1)is called the failure rate, and demonstrates that τk j∈N and βk j∈N are mutually indepen- q P(γk+1 =1|γk =0)the recovery rate. Assume that (A2) dent and have geometric distribution. Let us denote the predic- ∈ The failure and recovery rates satisfy p, q (0, 1). tion error covariance matrix at the stopping time βj by Pβj and 1006 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 63, NO. 4, APRIL 2018
Definition 3: The observability index Io is defined as the − smallest integer such that [C ,AC ,...,(AIo 1 ) C ] has rank n.IfIo =1, the system (C, A) is called one-step observable. L Sn → Sn We also introduce the operator K : + + defined as − I o 1 i−1 i (i) (i) ∗ i (i) (i) LK (X)=p (1−p) (A +K C ) ΦX (A + K C ) i=1 (6) where ΦX is the positive definite solution of the Lyapunov equa- 2 tion (1 − q)A ΦX A + qA XA =ΦX with |λA | (1 − q) < 1, − and K =[K(1),...,K(Io 1)] with each matrix K(i) having compatible dimensions. It can be easily shown that LK (X) is linear and nondecreasing in the positive semidefinite cone. We have the following result. 2 Theorem 1: Suppose |λA | (1 − q) < 1. If any of the follow- Fig. 2. Road map of this work. ing conditions hold. − i) ∃ K [K(1),...,K(Io 1)], where K(i)s are matrices with compatible dimensions, such that |λH (K )| < 1, 1 call it the peak covariance at βj . We introduce the notion of where PCS [29] as follows: − − 1 Definition 2: The Kalman filtering system with packet losses H(K)= qp (A ⊗ A) 1 − (1 − q)I is said to be peak-covariance stable if sup ∈N EP < ∞. j βj − I o 1 Remark 1: In [32], the authors defined stability in stopping · (Ai + K(i)C(i)) times as the stability of Pk at packet reception times. Note that i=1 {βj }j∈N , at which the peak covariance is defined, can also be treated as the stopping times defined on packet reception times. ⊗ (Ai + K(i)C(i))(1 − p)i−1 . (7) Clearly, in scalar systems, the covariance is at maximum when − ii) There exists K [K(1),...,K(Io 1)] with each ma- the channel just recovers from failed transmissions; therefore, trix K(i) having compatible dimensions such that peak covariance sequence gives an upper envelop of covariance Lk ∈ Sn limk→∞ K (X)=0for any X + . matrices at packet reception times. For higher order systems, − iii) There exist K [K(1),...,K(Io 1)] with each matrix the relation between them is still unclear. (i) The results of this paper are organized as follows (cf., Fig. 2). K having compatible dimensions and P>0 such that L First of all, we present a relaxed condition for PCS (Theorem 1). K (P )
⎡ √ ⎤ − I −1 (i) X p(A Y + C F ) ··· p(1 − p)Io 2 (A o ) Y +(C ) F − √ 1 Io 1 ⎢ ∗ ··· ⎥ ⎢ p(YA+ F1 C) Y 0 ⎥ Ψ ⎢ ⎥ ⎣ . . .. . ⎦ . . . . − − ∗ p(1 − p)Io 2 YAIo 1 + F C(i) 0 ··· Y Io −1 (9) WU et al.: KALMAN FILTERING OVER GILBERT–ELLIOTT CHANNELS: STABILITY CONDITIONS AND CRITICAL CURVE 1007
dynamic and channel characteristics on the one hand, and PCS In light of Lemma 1, the spectrum of pn (A ⊗ A) is on the other hand. These results cover the ones in [28], [29], given by σ (pn (A ⊗ A)) = {pn (λi λj ): λi, λj ∈ σ(A)}. Since and [31], as is evident using the subadditivity property of matrix A is a real matrix, its complex eigenvalues, if any, al- 2 norm, and the fact that the spectral radius is the infimum of all ways occur in conjugate pairs. Therefore, |λA | must possible matrix norms. To see this, one should notice that be an eigenvalue of A ⊗ A, and the spectral radius of p ⊗ |λ | n |λ |2i − n (A A) can be computed as pn (A⊗A) = i=1 A (1 |λ | − H (k) i−1 Io 1 |λ |2j − j−1 q) q j=1 A (1 p) p. It is evident that the sequence −1 2 −1 {|λp ⊗ |} ∈N is monotonically increasing. When |λ | (1 − ≤ qp (A ⊗ A) − (1 − q)I n (A A) n A q) < 1,wehave − I o 1 i−1 i (i) (i) lim pn (A ⊗ A)=H(0) (11) (1 − p) (Ai + K(i)C(i)) ⊗ (A + K C ) n→∞ i=1 and ∞ I −1 i−1 i i |λ |2 o ≤ qp (1 − q) A ⊗ A q A 2j j−1 lim |λp ⊗ | = |λ | (1 − p) p. n (A A) 2 A i=1 n→∞ 1 −|λ | (1 − q) A j=1 − I o 1 (12) i−1 i (i) (i) i (i) (i) (1 − p) (A + K C ) ⊗ (A + K C ) As |λX | is continuous with respect to X, (11) and (12) alto- |λ |2 − i=1 |λ | q A Io 1 |λ |2j − j−1 gether lead to = 2 A (1 p) p. H (0) 1−|λA | (1−q) j=1 ∞ − I o 1 K(i) =0∀1 ≤ i ≤ I − 1 − − Letting o , the condition provided in − i 1 i 2 − i 1 i (i) (i)2 2 = q (1 q) A p (1 p) A + K C Theorem 1 becomes: (i) |λA | (1 − q) < 1, and (ii) |λH (0)| < 1. i=1 i=1 Since the left side of (10) is positive, it imposes the positivity of 1 −|λ |2 (1 − q), whereby the conclusion follows. in which the first inequality follows from |λH (K )|≤H(K) A and the submultiplicative property of matrix norms, and Although computationally friendly, Proposition 1 only pro- i ⊗ vides a comparably rough criterion. It can be expected that, the last equality holds because, for a matrix X, X (i) i i 2 given the ability of searching for K s on the positive semidef- X = λ i = λ i ∗ i = X . Compar- (X (X )i )⊗(X i (X ∗)i ) X (X ) inite cone, Theorem 1 is less conservative than Proposition 1; ison with the related results in the literature is also demonstrated this is demonstrated by Example I in Section V. by Example I in Section V. Remark 3: The left side of (10) is strictly positive when In addition, in [28], [29], and [31], the criteria for PCS are Io ≥ 2, while it vanishes when Io =1. In the latter case, plus the 2 difficult to check since some constants related to the operator necessity as shown in [30], |λA | (1 − q) < 1 thereby becomes g are hard to explicitly compute. A thorough numerical search a necessary and sufficient condition for PCS. This observation may be computationally demanding. In contrast, the stability is consistent with the conclusion of [31, Corollary 2]. check of Theorem 1 uses an LMI feasibility problem, which can Remark 4: Proposition 1 covers the result of [30, Th. 3.1]. often be efficiently solved. To see this, notice that In the following proposition, we present another condition for I −1 q|λ |2 o PCS, which is, despite being conservative, easy to check. The A |λ |2j (1 − p)j−1 p new condition is obtained by making all K(i)s in Theorem 1 1 −|λ |2 (1 − q) A A j=1 take the value zero. ∞ − Proposition 1: If the following condition is satisfied I o 1 2i i−1 2j j−1 = |λA | (1 − q) q |ΛA | (1 − p) p I −1 o i=1 j=1 pq|λ |2 |λ |2i (1 − p)i−1 < 1 −|λ |2 (1 − q), (10) A A A − ∞ I o 1 i=1 ≤ Ai 2 (1 − q)i−1 q Aj 2 (1 − p)j−1 p, then the Kalman filtering system is peak-covariance stable. i=1 j=1 Proof: The proof requires the following lemma. which implies the less conservativity of [30, Th. 3.1]. Lemma 1 ([39, Th. 1.1.6]): Let p(·) be a given polynomial. If λ is an eigenvalue of a matrix A, then, p(λ) is an eigenvalue IV. MEAN-SQUARE STABILITY of the matrix p(A). Define a sequence of polynomials of the matrix A ⊗ A as In this section, we will discuss MSS of Kalman filtering with {pn (A ⊗ A)}n∈N , where Markovian packet losses. n i i−1 A. From PCS to MSS pn (A ⊗ A)= (A ⊗ A) (1 − q) q i=1 Note that the PCS characterizes the filtering system at − I o 1 stopping times defined by (5), whereas MSS characterizes the × (A ⊗ A)j (1 − p)j−1 p. property of stability at all sampling times. In the literature, j=1 the relationship between the two stability notations is still an 1008 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 63, NO. 4, APRIL 2018 open problem. In this section, we aim to establish a connection between PCS and MSS. First, we need the following definition for the defective eigenvalues of a matrix. Definition 4: For λ ∈ σ(A) where A is a matrix, if the alge- braic multiplicity and the geometric multiplicity of λ are equal, then λ is called a semisimple eigenvalue of A.Ifλ is not semisim- ple, λ is called a defective eigenvalue of A. We are now able to present the following theorem indicating that as long as A has no defective eigenvalues on the unit cir- cle, i.e., the corresponding Jordan block is 1 × 1, PCS always implies MSS. In fact, we are going to prove this connection for general random packet drop processes {γk }k∈N , instead of limiting to the GilbertÐElliott model. Theorem 2: Let {γk }k∈N be a random process over an un- derlying probability space (S , S,μ) with each γk taking its value in {0, 1}. Suppose {βj }j∈N take finite values μ−almost surely, and that A has no defective eigenvalues on the unit circle. Fig. 3. Relationships between the PCS and MSS over the space of Then, the PCS of the Kalman filter always implies MSS, i.e., the filtering systems under consideration. Theorem 2 indicates that the E ∞ E ∞ supk∈N Pk < whenever supj∈N Pβj < . intersection of the two sets, PCS and NoDEUC, is contained in the set MSS. Note that {βj }j∈N can be defined over any random packet loss processes, therefore, the PCS with packet losses that the filtering system is undergoing remains in accord with Definition 2. Theorem 2 bridges the two stability notions of Kalman fil- tering with random packet losses in the literature. Particu- larly, this connection covers most of the existing models for packet losses, e.g., i.i.d. model [16], bounded Markovian [40], GilbertÐElliott [28], and finite-state channel [25], [26]. Although E E supk∈N Pk and supj∈N Pβj are not equal in general, this connection is built upon a critical understanding that, no matter to which interarrival interval between two successive βj s the time k belongs, Pk is uniformly bounded from above by an affine function of the norm of the peak covariances at the starting and ending points thereof. This point holds regardless of the model of packet loss process. The proof of Theorem 2 was given in Appendix C. We also remark that there is difficulty in relaxing the assump- tion that A has no defective eigenvalues on the unit circle in Theorem 2. This is due to the fact that As defective eigenvalues Fig. 4. p − q plane is divided into MSS and non-MSS regions by the on the unit circle will influence both the PCS and MSS in a critical curve fc (p, q)=0. When p + q =1, the Markovian packet loss process is reduced to an i.i.d. process. As a result, the intersection point nontrivial manner (see Fig. 3). of the curves fc (p, q)=0and p + q =1gives the critical packet drop Remark 5: In [29], for a scalar model with i.i.d. packet probability established in [16]. losses, it has been shown that the PCS is equivalent to MSS, whereas for a vector system even with i.i.d. packet losses, the 2 relationship between the two is unclear. In [31], the equivalence Proposition 2: Let the recovery rate q satisfy |λA | (1 − q) between the two stability notions was established for systems < 1. Then, there exists a critical value pc ∈ (0, 1] for the failure that are one-step observable, again for the i.i.d. case. Theorem 2 rate in the sense that E ∞ ≥ now fills the gap for a large class of vector systems under general i) supk∈N Pk < for all Σ0 0 and 0
Theorem 3: There exists a critical curve defined by fc (p, q)=0, which reads two nondecreasing functions p = · −1 · 2 pc (q) and q = qc (p) with qc ( )=pc ( ), dividing (0, 1) into two disjoint regions such that ∈ E ∞ i) If (p, q) fc (p, q) > 0 , then, supk∈N Pk < for all Σ0 ≥ 0; ii) If (p, q) ∈ fc (p, q) < 0 , then, there exists Σ0 ≥ 0 un- E ∞ der which supk∈N Pk = . Remark 6: If the packet loss process is an i.i.d. process, where p + q =1 in the transition probability matrix defined Fig. 5. Sample path of P with p =0.99, q =0.65 in Example I. in (2), Proposition 2 and Theorem 3 recover the result of k [16, Th. 2]. It is worth pointing out that whether MSS holds or not exactly on the curve fc (p, q)=0is beyond the reach of the current analysis (even for the i.i.d. case with p + q =1): such an understanding relies on the compactness of the stability or nonstability regions.
C. MSS Conditions We can now make use of the PCS conditions we obtained in the last section, and the connection between PCS and MSS Fig. 6. Sample path of Pk with p =0.45, q =0.5 in Example II. indicated in Theorem 2, to establish MSS conditions for the con- sidered Kalman filter. It turns out that the assumption requiring p, [29] concludes that p<0.04 guarantees PCS; while Proposi- 2 no defective eigenvalues on the unit circle can be relaxed by an 1−|λA | (1−q) tion 1 requires p< 4 , which generates a less conser- |λA | q approximation method. We present the following result. vative condition p<0.22. Similarly, the example in [30] allows |λ |2 − Theorem 4: Let the recovery rate q satisfy A (1 q) < 1. p =0.1191 at most. We also note that it is rather convenient to ≤ Then, there holds p pc , where check the condition in Proposition 1 even with manual calcula- tion; in contrast, to check the conditions in [29] and [30] involves p sup p : ∃(K, P)s.t. L (P )
0 (13) K a considerable amount of numerical calculation. Additionally, we use the criterion established in Theorem 1 to check for the i.e., for all Σ0 ≥ 0 and 0