Research Article | Open Access
Distributed Robust Kalman Filtering with Unknown and Noisy Parameters in Sensor Networks
This paper investigates the distributed filtering for discrete-time-invariant systems in sensor networks where each sensor’s measuring system may not be observable, and each sensor can just obtain partial system parameters with unknown coefficients which are modeled by Gaussian white noises. A fully distributed robust Kalman filtering algorithm consisting of two parts is proposed. One is a consensus Kalman filter to estimate the system parameters. It is proved that the mean square estimation errors for the system parameters converge to zero if and only if, for any one system parameter, its accessible node subset is globally reachable. The other is a consensus robust Kalman filter to estimate the system state based on the system matrix estimations and covariances. It is proved that the mean square estimation error of each sensor is upper-bounded by the trace of its covariance. An explicit sufficient stability condition of the algorithm is further provided. A numerical simulation is given to illustrate the results.
Distributed estimation is a fundamental problem in sensor networks and has attracted broad attention of the researchers. Earlier works on distributed estimation often assume that the system model is precisely known to all sensors. However, in most practical applications, an exact model of the system may not be available to sensors and the robust performance of the distributed filter should be investigated.
Consensus is a simple and feasible strategy for distributed estimation over networks. There have been many consensus estimation algorithms. Olfati-Saber () proposed consensus Kalman filtering algorithms based on the consensus protocol and standard Kalman filter. Stanković et al. () gave an algorithm which was composed of decentralized overlapping estimators and a consensus scheme for sensor networks where each sensor had prior knowledge of submatrices of the system matrices. Li et al. () developed optimal and suboptimal filters and gave sufficient conditions for the stochastic stability of the suboptimal filter. Cattivelli et al. () studied the networks in which each sensor had observable aggregated measurement matrix and proposed diffusion Kalman filtering algorithms. Yang et al. () studied sensor networks where each sensor could measure the target and was activated with a certain probability. Ugrinovskii ([6, 7]) proposed a suboptimal consensus-based estimation algorithm for continuous-time system with deterministic disturbance. The above-mentioned algorithms are called one time-scale method; i.e., the consensus iteration time interval is equal to the estimation/filtering one. References [8–13] designed consensus Kalman filtering algorithms with two time-scale strategy; i.e., during one estimation/filtering interval, multiple consensus iterations were processed to fuse the prior and novel information. The two time-scale algorithms are fully distributed and their stability just requires the collective detectability of the network. However, in this algorithm vast communication costs are required. In most existing works on the one time-scale algorithm, to design the observer gains or provide stability conditions, global knowledge of the network topology is needed. References [2–7, 14–18] gave stability conditions based on LMIs or the spectral radius of the whole system matrix, which requires global topology information and large computations.
Earlier works on distributed estimation in sensor networks often suppose that an exact model of the target system is known. In most practical problems, it is difficult for all sensors to exactly know the model of the system, and each sensor may just have access to partial model information. The robust performance of the filter with respect to the prior unknown or noisy system parameters is an important issue. Although there have been many works solving the robust state estimation problem ([19–21]), there have been few works designing the consensus-based robust estimation algorithm and investigating the robust performance for sensor networks. Recently, Zhang et al. () proposed a distributed Kalman/ filter for sensor networks with partial unknown system parameters. Han et al. () proposed a distributed consensus filtering problem for a class of discrete-time-varying systems with stochastic nonlinearities and gave sufficient stability conditions based on the recursive linear matrix inequalities. However, the linear matrix inequalities of sensors are correlated and their verification requires global knowledge of topology and large computations.
In this paper, we mainly study the distributed filtering for discrete-time-invariant system whose system matrix parameters are previously unknown or not precisely known to the sensors in the network. Each sensor may not be able to obtain the measurement. Even if one sensor can obtain the measurement, its measurement may not be observable. We design a distributed robust Kalman filtering algorithm to estimate the system parameters and system state. To estimate the system parameters, due to noises in the available parameter information, we apply a consensus Kalman filtering algorithm. We prove that the estimation errors of the system parameters converge to zero in mean square sense, if and only if for any one of the system parameters its corresponding node subset in which each node can obtain this parameter information is globally reachable. To estimate the system state, due to parameter estimation errors, we apply a consensus robust Kalman filtering algorithm where the consensus weights are designed based on the covariances of the sensors. We provide a sufficient condition guaranteeing the boundedness of the mean square estimation errors of the sensors. The contributions of this paper include the following: (1) to propose a fully distributed robust Kalman filtering algorithm to deal with the noises in available system parameters; (2) to provide a sufficient stability condition based on some uncorrelated LMIs, whose computation is small and does not require global topology information.
Notation. In this paper, denotes a unit matrix of size . , , and represent the trace, the trace norm, and the 2-norm of a matrix, respectively. denotes the 2-norm of vector . denotes a block-diagonal matrix with diagonal blocks . Matrix is positively definite and is written as . denotes the expectation of . For a finite node set , denotes the number of nodes in this set.
2. Problem Formulation
Consider a target system with discrete-time linear dynamicswhere denotes the state of the system and is a system Gaussian noise with zero mean and covariance matrix . It is supposed that is nonsingular.
Suppose target (1) is monitored by a network of sensors with the following measuring system:where, for sensor , is the measurement vector at sensor , is the measuring matrix, and is the Gaussian measurement white noise with zero mean and covariance matrix . The measurement noises of different sensors and system noise are uncorrelated. Here may not be observable.
In the network, some sensors may not be able to obtain the measurements (2) due to limited sensing range and energy saving. Let denote whether the measurement is available. If sensor can obtain its measurement (2), then ; otherwise . Let () represent the sensing topology.
In applications, it is a really strong assumption that all sensors in the network precisely know the system parameters of the target. In many cases, if one sensor can just obtain the measurement outputs with respect to partial system state, it is likely to obtain the corresponding partial parameters in system matrix . This paper discusses the case in which the parameters may not be previously known to some sensors and sensors’ available system information contains noises. For (), we apply to denote whether the rank, column element of matrix is known to sensor . If the information is unknown to , ; otherwise, and sensor ’s corresponding available system information iswhere and is a zero mean Gaussian white noise with variance , which is independent of , , and for any one different , , or ().
In this network, even if one sensor can obtain the measurement of the target, it does not completely know the target’s system matrix and its measuring system may be also not observable, which makes it infeasible to estimate the state of the target without cooperation. The objective is to construct an estimation algorithm on the basis of the local system parameter information, local measurements (when available), and data received from all adjacent neighbors to estimate the state of the target.
3. Distributed Filtering Algorithm
To explain the principle of graph theory in distributed estimation, preliminaries about graph theory are firstly introduced.
In the cooperative estimation in the sensor network, each sensor is treated as a node and the nodes communicate according to the communication graph. For a communication digraph , denotes the node set, denotes the edge set, and denotes the adjacency matrix. If node can receive information from node , then there is a corresponding edge in ; i.e., . The adjacency elements associated with the edges of the graph are defined as , . Since each node can always get its own information, for all , . Neighbor set of sensor is denoted as .
In , a simple path of length from to is such that there exists a sequence of nodes with subsequent edges . For node and a node subset , there exists at least one path from to the node set if there exists at least one node such that there is a path from node to node . A node subset is said to be globally reachable in the communication topology if, for any node in , there exists at least one path from to the node set .
Consensus protocol is a simple and effective strategy in cooperation of sensors. In this distributed filtering algorithm, two consensus processes, aimed at estimating the system parameters and system state, respectively, are included. Let us assume that, at time , each node in the sensor network can transmit its system matrix information, including parameter estimation error covariances and parameter estimations , and its system state estimation information, including measurement , measurement matrix , measurement noise covariance (when available), state estimation error covariance , state covariance , and estimated state , to its neighbors. Then the consensus-based distributed robust Kalman filtering algorithm is summarized by Algorithm 1.
For node , when new data is received, it firstly locally computes the weighted average of its neighbors’ estimations and covariances for the parameters in the system matrix:where and denote the fused estimation and covariance for system matrix, respectively, is the weight of edge and can be any positive value satisfying , and denotes the neighbor set of sensor . And then a local Kalman filtering is constructed to estimate the system parameters.
In the second part, sensor locally aggregates the measurement vectors, measurement matrices, and measurement noise covariance of its neighbors , respectively, and obtains , , . Then it computes the weights of in-edges () and obtains the weighted average of its neighbors’ covariance matrices and estimated states
In the third part, a local robust Kalman filter is constructed based on the aggregated measurements and weighted average estimate state for each sensor node to estimate the state of the target. For sensor , where and are the filter gains to be determined. Assume .
The weights in the consensus estimation of the system matrix can be any positive values with row sum being equal to 1, while in the consensus estimation of the system state, due to the dynamical variation of the state, the weights cannot be arbitrary. When the weights are not properly designed, the estimation error of the overall network may be divergent. In [6, 7, 14–18] the weights are constants and concerning conditions are implied in the LMIs-based stability conditions, which is difficult to be computed. Here we apply a fully distributed weight design approach based on the traces of the estimation covariances of the sensors.
At time , for one node , if is observable, then the weight of edge is designed as and for node , ; otherwisewhere denotes the trace of covariance matrix .
Remark 1. The weighted average consensus-based Kalman filter in Algorithm 1 is a little similar to that proposed in  where Matei et al. studied sensor networks with completely known parameters and applied a distributed suboptimal filter with weighted average consensus protocol and standard Kalman filter. For sensor ,andBy Markovian jump linear system, a sufficient stability condition is given based on the feasibility of the following LMIs:If the system matrix is not precisely known to sensors and just its estimation can be used, when estimating the system state, standard Kalman filter is not feasible and a robust Kalman filtering algorithm is designed in this paper.
Remark 2. In the network, when one sensor can obtain the measurement outputs of the process, it just obtains partial parameters of the system matrix with noises. Thus, Algorithm 1 contains two consensus-based filtering processes. One is a consensus Kalman filtering for system parameters. The other is a distributed filtering for system state. Moreover, being different from existing algorithms ([6, 7, 14–18]), the weights in the consensus filtering for system state are adaptively designed. The correction step for filtering is done at sensor levels, and thus the filtering algorithm is completely distributed.
4. Stability Analysis
In this section, we analyze the stability properties of the distributed robust Kalman filtering algorithm in Algorithm 1. In this paper, we assume the network is time-invariant.
Firstly, we analyze the estimation of the system parameters. For sensor and , from Algorithm 1 we have thatwithwhere (). If sensor has an access to the parameter with noise, ; otherwise .
Proposition 3. For any and any node , holds if and only if in the communication topology there exists at least one path from to the node set , where .
(Sufficiency). If in the graph there exists one node , , and node is globally reachable from any other node, then we partition the graph into strongly connected subgraphs (, . Renumber the nodes such that is an upper block-diagonal matrix and the node subset corresponding to contains node . Define as a vector stacking the values , , , and define as a vector stacking the values , . Then,For irreducible stochastic matrix , there exists a vector with positive elements satisfying that . Thus, . Since the subgraph is strongly connected, if for one node , , then from (13), we have that, for any other node , holds. When , since , we have that and then , which infers that , and thus for , .
For matrix , there exists a diagonal matrix , in which the diagonal elements are all larger than 1 and there exists at least one diagonal element being strictly larger than 1, such that is an irreducible stochastic matrix. Define with positive elements satisfying that , and define ; then for , , and there exists at least one ; holds. Thus, if for all , , From the above analysis, for , ; thus, for , also holds.
Similarly, for , . Therefore, if in the graph there exists one node , , and node is globally reachable from any other node, then holds for all .
If there is no globally reachable node with , then the graph can be partitioned into subgraphs; in each subgraph there exists one node denoted as , , and node is globally reachable from any other node in the subgraph. From the above analysis, it is not difficult to prove that holds for all . The sufficiency part has been proved.
(Necessity). If there exists a node having no path from itself to the node set , then there exists a node subset in which each node has no neighbor in other node subsets and is not be able to obtain the information of . Obviously, the prediction error covariances of the nodes in this subset cannot converge to 0. The necessity part has been proved.
In the following, we analyze the estimation for the system state. To begin with, we give a lemma about robust Kalman filtering for a system with unknown random uncertainty modeled by white noise.
Lemma 4. Consider a target with discrete-time system:where is the state; is the measured output with observable measurement matrix ; and are uncorrelated zero mean Gaussian white noises with variance matrix and , respectively. The observer cannot obtain the precise system matrix and just obtains the value with noise; i.e., , where represents the unknown parameter noise and satisfies , for any or , . If there exist such thatwhere , satisfy the following two Riccati equationsandthen, the filterwith satisfiesif .
Proof. Define , ; the augmented system of (19) and (24) can be expressed bywhereFrom the properties of the parameter noises, we have . Define ; it is obvious that .
Define , then from (26),whereand satisfy .
From , if there exists such that , thenTherefore,Define an equationwith . Obviously, holds for all . Thus, implies .
If is assumed to have the form , from (32) we have , whereWhen , there holds , and thus is a block-diagonal positive definite matrix with diagonal blocks.In order to find the optimal filter gains minimizing , we take the first variation to (34) with respect to and equal it to zero.From (36) we obtain the optimal gains . This completes the proof.
In the following we analyze the stability of the distributed robust Kalman filtering algorithm.
For sensor , , it can only use its own estimated system matrix , and then the target dynamics can be depicted byHere and satisfies for any or , .