Abstract
The average consensus problem of distributed inference in a wireless sensor network under Markovian communication topology of uncertain transition probability is studied. A sufficient condition for average consensus of linear distributed inference algorithm is presented. Based on linear matrix inequalities and numerical optimization, a design method of fast distributed inference is provided.
1. Introduction
During the past few decades, consensus problems of multiagent systems by information exchange have been extensively studied by many researchers, due to their widespread applications in autonomous spacecraft, unmanned air vehicles, mobile robots, and distributed sensor networks. OlfatiSaber and Murray introduced in [1, 2] a theoretical framework for solving consensus problems. In [3, 4], consensus problems of firstorder integrator systems were proposed based on algebra graph theory. In [5, 6], consensus problems of directed secondorder systems were presented. In [5], the authors provided necessary and sufficient condition for reaching mean square consensus of discretetime second order systems. Consensus conditions were studied in [6] of continuoustime second order systems by Linear Matrix Inequality (LMI) approach.
Among consensus problems, the average consensus problem is challenging which requires distributed computation of the average of the initial state of a network [1, 2]. For a strongly connected network, [1] proved that the average consensus problem is solvable if and only if the network is balanced. The discretetime average consensus plays a key role in distributed inference in sensor networks. In networks of fixed topology, [7] gave necessary and sufficient conditions for linear distributed inference to achieve average consensus. A design method was presented in [7] to implement linear distributed inference of fastest consensus. Because of noisy communication channels, link failures often occur in a real network. Therefore, it is meaningful to study distributed inference in networks of swing topology. Through a common Lyapunov function, a result of [1] stated that distributed inference reaches average consensus in a network of swing topology if the network holds strongly connected and balanced topology. Reference [8] modeled a network of swing topology using a Bernoulli process and established a necessary and sufficient condition for average consensus of distributed inference. The condition is related to a mean Laplacian matrix.
The Bernoulli process in [8] means that the network link failure events are temporally independent. From the viewpoint of engineering, it is more reasonable to consider network link failures of temporal independence. The most famous and most tractable stochastic process of temporal independence is Markovian chain in which any future event is independent of the past events and depends only on the present event. This motivates us to model a network of swing topology using a Markovian chain and hence to study distributed inference using Markovian jump linear system method [9–12]. In practice, transition probabilities of a Markovian chain are not known precisely a priori, and only estimated values of transition probabilities are available. Hence, this paper thinks of networks with Markovian communication of uncertain transition probability.
In fact, in the research on networked control systems, Markovian chain has been used by several authors to describe random communication in networks. Reference [13] provided packetloss model by Markovian chain in networked control. Under network communication of update times driven by Markovian chain, [14] gave stability conditions of modelbased networked control systems. Networked control systems with bounded packet losses and transmission delays are modeled through Markovian chain in [15]. The networked predictive control system in [16] adopted 2 Markovian chains to express date transmission in both the controlleractuator channel and the sensorcontroller channel.
In this paper, is used to denote the set of all nonnegative integers. The real identity matrix of is denoted by . Let 1 be the vector whose elements are all equal to 1. The Euclidean norm is denoted by . If a matrix is positive (negative) definite, it is denoted by (<0). The notation within a matrix represents the symmetric term of the matrix. The expected value is represented by .
The paper is organized as follows. Section 2 contains a description of the network and linear distributed inference. Section 3 presents average consensus conditions and a design method. Numerical simulation results are in Section 4. Finally, Section 5 draws conclusions.
2. Network Description
Consider distributed inference in a wireless sensor network consisting of sensor. Each sensor collects a local measurement about the situation of environment. It is assumed that these local measurements are independent and identically distributed random variables. The goal of inference is for all sensors to reach the global measurement such that the true situation of environment can be monitored convincingly.
This paper studies iterative distributed inference. Define which includes all realizable undirected links in the wireless sensor network. At the th iteration , the successful communication links in the wireless sensor network are described by the set A pair means that the sensors and communicate with each other at . A pair but means that there is no communication link between the sensors and at . Due to noisy communication channels and limited network power budget, is assumed to be random and to be modeled as follows. Given distinct subsets . Let be a stochastic process taking values in and driven by a Markov chain with a transition probability matrix , where , for all , for all . However, these s are not known precisely. Each is expressed as with a known and a unknown whose absolute value is less than a given positive constant . For any for all , and . This paper models The neighborhood of sensor at is denoted by and the element number of set is denoted by .
For sensor , set its initial state . At the th iteration, each sensor obtains its neighbors’ states and updates its state using the following linear iteration law: where is the weight parameter which is assigned by designers. Our study on the above distributed inference has two objectives: one is to derive a condition on the convergence of to in the sense of mean square; the other is how to find such that a fast convergence is achieved.
3. Average Consensus Analysis
3.1. Convergence Condition
Denote The system in Section 2 can be described as where is a matrix with entries . For , For , From (2)~(6), it is seen that if and only if , and hence
Furthermore, (10)~(12) imply that for all Then, we have which means that for all Denote . The iterative distributed inference is said to be average consensus in mean square sense if for any initial condition and .
Theorem 1. The linear distributed inference (7) reaches average consensus by choice of , if there exist m positive definite matrices such that for all with .
Proof. Assume at time step . From (9), (14), and (16), one has
We now consider the stochastic Lyapunov function
Then for all , we have
Denote
When the conditions in Theorem 1 are satisfied, we have
where denotes the minimal eigenvalue of and
Therefore, for all , for all , for all ,
This means .
3.2. Optimal Design
From the above proof, it is seen that the conditions in Theorem 1 result in not only but also decreasing . Moreover, implies ; that is, also converges to zero. Therefore, the decrease rate of can express the convergence speed of distributed inference. The following theorem is about the decrease rate of .
Theorem 2. Given and , if there exist m positive definite matrices such that for all ,then in linear distributed inference (7), for any nonzero ,
Proof. According to Schur complement [17], condition (25) can be rewritten as
From (27), for any nonzero , one has
for all and for all . From (28), (19), and (20), it is known that
and hence
Condition (25) in Theorem 2 is a LMI. We denote condition (25) as and for any define Using the LMI toolbox of MATLAB, can be computed by Algorithm 1.

For a with , it is known from Theorems 1 and 2 that linear distributed inference (7) is average consensus and that is a bound of convergence speed. Since a less value of gives a faster convergence speed, the fast distributed inference problem is addressed as which is an unconstrainted optimization problem of only one variable. Many existing numerical optimization methods [18] can be utilized to solve this problem efficiently. When , the optimal parameter provides a fast linear distributed inference which reaches average consensus.
4. Numerical Example
In this section, we present simulation results for average consensus of distributed inference in a simple sensor network. The network has 10 sensor nodes and switched in three possible communication situations. Figure 1 illustrates 3 communication situations. The estimated transition probabilities of is The estimate error satisfies Using the computation procedure in Section 3, the optimization problem (32) is solved. Graph of is displayed in Figure 2. The result is and . For the communication situation in Figure 1(1), we use the design method in [7] of minimizing asymptotic convergence factor and obtain . The design method in [7] is also applied to the other 2 situations in Figure 1 and get .
In order to compare our method with that in [7], the initial states of each sensor node is selected as
Thus from (1), we have . The real transition probability matrix is set as Figures 3, 4, and 5 show state curves of all sensor nodes under , , and , respectively. It can be seen that all sensor states convergence to and that our has faster convergence rate than or has. Achieving faster convergence is because our method considers the random switching among the 3 communication situations while [7]’s method considers only 1 communication situation.
5. Conclusion
The distributed average consensus problem in sensor networks has been studied under a Markovian switching communication topology of uncertain transition probabilities. Stochastic Lyapunov functions have been employed to investigate average consensus of linear distributed inference. A sufficient condition of average consensus has been proposed based on feasibility of a set of coupled LMIs. The design problem of fast distributed inference has been solved by numerical optimization techniques.
Acknowledgment
This work was supported by 973 Program of China (Grant 2009CB320603).