Discrete Dynamics in Nature and Society

Volume 2015, Article ID 160683, 9 pages

http://dx.doi.org/10.1155/2015/160683

## Robust Distributed Filtering for Nonlinear Systems with Sensor Saturations and Fractional Uncertainties with Digital Simulation

State Grid Smart Grid Research Institute, North Zone of Future Sci-Tech City, Beiqijia, Changping District, Beijing 102211, China

Received 13 February 2015; Revised 30 March 2015; Accepted 5 April 2015

Academic Editor: Zidong Wang

Copyright © 2015 Dong Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This paper is concerned with the robust distributed filtering problem for nonlinear systems subject to sensor saturations and fractional parameter uncertainties. A sufficient condition is derived for the filtering error system to reach the required performance in terms of recursive linear matrix inequality method. An iterative algorithm is then proposed to obtain the filter parameters recursively by solving the corresponding linear matrix inequality. A numerical example is presented to show the effectiveness of the proposed method.

#### 1. Introduction

The sensor networks have gradually received more and more research interest in the past decades. A typical example is the power grids which have evolved over the past century from a series of small independent community based systems to large-scale and complex systems involving many kinds of intelligent electronic devices [1]. As is well known, the information collected by the sensor networks need to be further processed before used for diagnosis, control, optimization, and other transactions for the smart grids, which requires much more advanced methodologies that are more sensitive, reliable, and economic than those applied in the traditional centralized electricity networks. Due to its clear engineering insights, recently, sensor networks have gained increasing research interests during the past decades in various branches of theoretical research and industrial applications. In particular, the distributed filtering over sensor networks has been ongoing research issue that attracts special attention from researchers in the area. Different from traditional filtering techniques based on single or centralized structured/located sensors, the information available for the filter design on an individual node of the sensor network is not only from its own measurement but also from its neighboring sensors’ measurements according to the given topology [2]. As such, the objective of filtering based on a sensor network can be achieved in a distributed yet collaborative way. Such a problem is usually referred to as the distributed filtering problem. It is noticed that one of the main challenges in designing distributed filters lies in how to cope with the complicated couplings between one sensor and its neighboring sensors by reflecting such couplings in the filter structure specification.

It is worth mentioning that in real-world applications, due to a variety of reasons, such as abrupt changes of working conditions, internal/external disturbances, erosions, equipment aging, the sensor outputs are usually corrupted by these different disadvantages. Such phenomena are always referred to as incomplete information. Up to now, filtering problems with incomplete information have been widely studied and a lot of results have been reported in the literature; see [3–14] for some latest publications. As far as the distributed filtering problems subject to incomplete information are concerned, some recent representative work can be summarized as follows. In [15], the distributed filtering problem was solved when randomly occurring saturations and successive packet dropouts appear over sensor networks. In [16], distributed filtering problem was studied for a class of Markovian jump nonlinear time-delay systems over lossy sensor networks. In [17], distributed filtering for a class of time-varying systems was discussed over sensor networks with quantization errors and successive packet dropouts. The distributed state estimation with stochastic parameters and nonlinearities through sensor networks was solved in [18] over the finite-horizon. With randomly varying nonlinearities and missing measurements, [19] developed a distributed filtering method to achieve the performance requirements. In [20], distributed state estimation problem for uncertain Markov-type sensor networks with mode-dependent distributed delays was solved. Using linear matrix inequality approach, the distributed filters were given in [21, 22], respectively, for different types of stochastic systems. However, it is worth mentioning that, up to now, the distributed filtering problem has not been studied for nonlinear time-varying systems subject to randomly occurring sensor failures.

Motivated by the above discussion, in this paper, it is the objective to design a robust filter for a class of discrete time-varying nonlinear stochastic systems subject to sensor saturations and fractional uncertainties. The contribution of this paper is twofold: (i) for the first time, the filtering problem is discussed while taking both sensor saturations and fractional uncertainties into consideration. (ii) An iterative algorithm is developed to solve the robust filtering problem and then seek the desired filter parameters step by step. The rest of this paper is organized as follows: Section 2 formulates the distributed event-based filter design problem for nonlinear discrete time-varying stochastic system. The main results are presented in Section 3 where sufficient conditions for the existence of the desired filter are given in terms of recursive linear matrix inequalities. Section 4 gives a numerical example. Section 5 is the conclusion.

#### 2. Problem of Formulation

In this paper, it is assumed that the sensor network has sensor nodes which are distributed in the space according to a specific interconnection topology characterized by a directed graph , where denotes the set of sensor nodes, is the set of edges, and is the nonnegative adjacency matrix associated with the edges of the graph, that is, , which means that there is information transmission from sensor to sensor . If , then node is called one of the neighbors of node . For all , denote , which means that, in the sensor network, sensor node can receive the information from its neighboring nodes according to the given network topology.

Let us consider the discrete-time nonlinear stochastic system with sensors defined on :where , , , and represent the state, the measured output of the th sensor, estimated output, and disturbance belonging to , respectively. , , , , and are known real time-varying matrices with appropriate dimensions.

*Assumption 1. *The nonlinear function is assumed to obey the following constraint:where are known positive definite matrix sequence with appropriate dimensions describing the shape of the ellipsoids with being the center of the ellipsoids and is a known positive scalar.

The nonlinear function stands for the nonlinearity that is unknown, bounded, and deterministic but reside within an ellipsoidal set. Such a type of nonlinearity is usually occurring in practical engineering practice and is always a main origin for the degradation of the system performance.

*Assumption 2 (see [23]). * has the linear fractional form as follows:with and , where , , and are known constant matrices and denotes the unknown matrix functions with Lebesgue measurable elements.

*Definition 3. *A nonlinear function is said to satisfy the sector-bounded condition iffor some real matrices and , where is a symmetric positive definite matrix. In this case, we say that belongs to .

As discussed in [23], there exist matrices and such that ; the sensor saturation function can be written aswhere is a nonlinear vector-valued function satisfying the sector-bounded condition with and . In this case, can be represented as follows:where .

The objective of this paper is to design a filter with the following form for system (1):where is the state estimate, is the output to be estimated, and and are the filter parameters to be determined.

From system (1) and filter structure (7), one can obtainwhich can be rewritten in a more compact form as follows:where

Define ; it is the objective of this paper to design the filter parameters and to satisfyfor some given disturbance attenuation level and given positive matrix .

#### 3. Main Results

Before giving the main results, some useful lemmas are firstly introduced.

Lemma 4. *Let be quadratic functions of the variable η: with . If there exist scalars such that*

*then the following is true*

*Lemma 5 (Schur Complement Equivalence). Given constant matrices , , and , where and , thenif and only ifor, equivalently,*

*Lemma 6 (see [23]). Let , , and be real matrices of appropriate dimensions with satisfying . Then, for any scalar , we have*

*3.1. Performance Analysis*

*3.1. Performance Analysis**The following theorem gives a sufficient condition under which the required performance can be satisfied.*

*Theorem 7. Given a performance and the initial matrix . Let the filter parameters and be given. If there exist three sequences of positive scalars , , and a sequence of positive definite matrices satisfying such that the following matrix inequalitywhere is feasible, then performance can be satisfied.*

*Proof. *Definingthen the filtering error system can be expressed byDefiningone can get from (9) thatwhere Adding the zero term to both sides of (23), one can getwhere Next, by summing up (26) on both sides from to with respect to , one can havewhich is equivalent toFrom (29), since and the initial condition , one can see that the performance can be satisfied if is satisfied.

Next, one can see from the ellipsoidal nonlinear function (2) thatwhere with .

Moreover, it also can be obtained from the sensor fault constraint (6) thatwhere with .

By Schur Complement Equivalence Lemma, the matrix inequality (18) can be rewritten as follows:By Lemma 4, can be satisfied, which means that the required performance can be satisfied. The proof is complete.

*Remark 8. *By means of iterative matrix inequality method, Theorem 7 presents a sufficient condition for filtering error system to achieve the desired specification subject to the sensor saturations. Moreover, notice that the obtained condition is expressed by certain matrix inequalities including fractional parameter uncertainties ; we cannot solve the matrix inequality directly to obtain the explicit parametric expression of the required filter parameters and . As such, in the following, we will provide a theorem with which the addressed problem can be solved and the filter structure can then be obtained step by step.

*3.2. Filter Design*

*3.2. Filter Design*

*Theorem 9. Let the performance and the initial condition matrix be given. If there exist a sequence of matrices , a sequence of matrices , three sequences of positive scalars , , and , a sequence of positive definite matrices , and a sequence of positive definite matrices satisfying such thatwhere with and , and is the solution to the following optimization problemthen the addressed filtering problem is solvable and the desired filter parameters at each time can be obtained by solving the corresponding matrix inequality.*

*Proof. *The proof is quite straightforward based on Theorem 7. First, define . One can obtain that the matrix inequality (18) is equivalent to the following inequality:In what follows, we will try to prove that inequality (38) can be satisfied if inequality (35) is fulfilled. To this end, rewrite (38) as follows:where From Lemma 6, it follows thatFurthermore, one can obtain thatwhich, by Schur Complement Lemma, is equivalent to the linear matrix inequality (35). Thus, inequality (18) is satisfied if (35) is met, and, according to Theorem 7, the required filtering problem is solved.

*Remark 10. *Theorem 9 outlines the methodology with which the addressed robust distributed filtering problem is solved for nonlinear systems subject to sensor saturations and fractional uncertainties. An iterative linear matrix inequality approach has been developed to calculate the desired filter parameters step by step. It should be pointed out that it will inevitably lead to the larger computing burdens. However, the algorithm we proposed in this paper is capable of dealing with the time-varying case, and the increase in computing burden can be largely reduced with the fast development of modern computer technology. In what follows, we will present a computing algorithm with which the numerical values of filter parameters can be obtained in each time step.

*3.3. Computational Algorithm*

*3.3. Computational Algorithm*

*Algorithm 11 ( Filter Design Algorithm). *

*Step 1*. Set the required , the initial condition , the initial , and . Choose properly the values of , , and . Select the initial values for and satisfying the initial condition constraint.

*Step 2*. Set . Solve (35) for , , , and . Then the and can also be obtained.

*Step 3*. Set . Using the obtained and , solve (35) for , , , and .

*Step 4*. If , then stop. Else go to Step 3.

*The Filter Design Algorithm gives a recursive way to obtain the numerical values of the desired filter parameters at each time point . It should be pointed out that the existence of the filter is expressed as the feasibility of certain linear matrix inequalities that can be solved forward in time. The possible research topic in the future is to consider more performance indices and give the filtering schemes that are able to satisfy multiple requirements simultaneously.*

*4. Simulation Example*

*4. Simulation Example*

*In this section, an illustrative example is presented to show the effectiveness of the proposed Filtering Design Algorithm.*

*Let us consider the nonlinear system (1) with the parameters given below: Let and set . Let .*

*Choose , , and , and we solve the optimaization problem given in Theorem 9 using Matlab toolbox and obtain . Set , , , and . Set the initial value of the system state and its estimate by and , and the initial values of and are selected to be . We can see the effectiveness and applicability of the proposed filtering algorithm from the simulation results (Figures 1-2).*