Discrete Dynamics in Nature and Society

Volume 2016, Article ID 1967875, 7 pages

http://dx.doi.org/10.1155/2016/1967875

## Positive Filtering with -Gain for Discrete-Time Positive Systems

College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China

Received 7 March 2016; Revised 30 March 2016; Accepted 7 April 2016

Academic Editor: Carmen Coll

Copyright © 2016 Xiaoming Chen. 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 positive filtering problem for discrete-time positive systems under the -induced performance. We aim to propose a pair of positive filters with error-bounding features to estimate the output of positive systems. A novel characterization is first constructed so that the filtering error system is asymptotically stable with a prescribed -induced performance. Then, necessary and sufficient conditions for the existence of required filters are presented, and the obtained results are expressed as linear programming problems. Moreover, it is pointed out that the results can be easily checked by standard software. In addition, a numerical example is given to show the effectiveness of the proposed design procedures.

#### 1. Introduction

In real world, many dynamical systems involve variables which are always confined to the positive orthant. This special category of systems is generally referred to as positive systems in the literature. Positive systems arise in different application fields such as physics, engineering, and social sciences [1, 2]. Since positive systems possess many unique features and have special structures, a lot of methods established for general systems cannot be used for positive systems. Due to their numerous applications and unique features, positive systems have received ever-increasing research interest in recent years [3–12]. For instance, the problem of reachability and controllability has been studied for positive systems in [13–15]. In [16, 17], the state-feedback controller synthesis results have been expressed as linear matrix inequality problems and linear programming problems. In [18, 19], a stability analysis method for compartmental dynamic systems has been proposed. An interesting result on the positive observer design problem has been given for positive systems in [20]. The positive realization problem has been discussed in [21]. In [22], a solution has been proposed to the model reduction problem for positive systems. Furthermore, the analysis and synthesis problems have been addressed for special classes of positive systems such as 2-D positive systems [4, 23] and time-delay positive systems [24–29].

It is remarkable that since many previous approaches used for the filtering problem of general systems fail to ensure the positivity of the filter, existing approaches cannot be directly applied for positive systems. Therefore, it is necessary to develop new techniques for positive systems. Moreover, differently from most existing results on stability and stabilizability of positive systems which were derived by resorting to the quadratic Lyapunov functions, the applications of the linear copositive Lyapunov functions led to many novel results in recent years [1, 16, 30–36]. In addition, linear copositive Lyapunov functions stimulate the use of -gain as a performance index for positive systems. In some situations, one might be interested in the sum of quantities for positive systems from a practical viewpoint. Therefore, some frequently used performance measures induced by signals such as norm are not very natural to describe the features of practical positive systems. On the other hand, -norm is more appropriate for positive systems because it represents the sum of the values of the components. For example, it is usually desirable to analyze the total mass of material in all the compartments for compartmental networks.

In this paper, the problem of -induced filtering is studied for positive systems with the positivity preserved in the filters. The main contributions of this paper are as follows. First, we present the -induced performance index for positive systems and characterize it analytically. Then, we propose a pair of -induced positive filters to estimate the output of positive systems at all times. Finally, we establish necessary and sufficient conditions to obtain the desired positive filters in terms of linear programming.

The layout of the paper is as follows. In Section 2, the problem addressed in this paper is formulated and the positive filters are introduced. The positive filter design procedure for positive systems is proposed in Section 3, followed by a numerical example in Section 4 to show the application of the theoretical results. Finally, we draw conclusions in Section 5.

*Notation 1. *All the matrices, if the dimensions are not explicitly stated, are assumed to have compatible dimensions for algebraic operations. Let denote the set of real numbers; is the -column real vectors; stands for the set of real matrices of dimension For a vector , the -norm is defined as . Let denote the nonnegative orthant of . For a matrix , the element located at the th row and the th column is denoted by ; and denote the th row and the th column, respectively. The -norm of an infinite sequence is defined as . (resp., ) means that, for all and , (resp., ). The induced 1-norm of a matrix is denoted by . The notation means that the matrix and denotes . The Euclidean norm for vectors can be expressed as . The space of all vector-valued functions defined on with finite norm is denoted by . stands for a column vector with each entry equal to .

#### 2. Problem Formulation

In this section, the -induced performance is first introduced for discrete-time positive systems. Moreover, the -induced filtering problem for a stable positive system is formulated.

Consider the following discrete-time linear system:where , , , and denote the state vector, disturbance signal, measurement, and the signal to be determined, respectively.

Next, the following definition is presented below, which will be used in the sequel.

*Definition 1. *System (1) is positive if , , and always hold for all and all input .

Before moving on, some useful results are introduced and the following lemmas are needed.

Lemma 2 (see [37]). *The discrete-time system (1) is positive if and only if*

Lemma 3. *System (1) with input is asymptotically stable if and only if there exists a vector (or ) satisfying*

Next, the definition of -induced performance is introduced. A stable positive system (1) is said to have -induced performance at the level if, under zero initial conditions,where is a given scalar.

Now we are in the position to introduce the following result which serves as a characterization on the asymptotic stability of system (1) with the -induced performance in (4) via linear programming. The performance characterization result is a theoretical basis for further development.

Lemma 4 (see [32]). *The positive system (1) is asymptotically stable and satisfies if and only if there exists a vector satisfying*

It is easy to see that the transient output cannot be estimated via conventional filters, which can only give an estimate of the output asymptotically. In order to design a filter which can be used to give the information of the transient output at all times, we intend to find a lower-bounding estimate and an upper-bounding one . More specifically, the signal can be encapsulated between the two estimates at all times. A pair of filters is proposed as follows:where , , , and . , and are filtering parameters.

First, the lower-bounding case is considered. Set , , and . From (1) and (6), the augmented system is described as follows:where

Since, the lower-bounding filter (6) is designed to approximate with , it is natural to require that the estimate is also positive, which means that the filter (6) should be a positive system. From Lemma 2, it can be seen that , , and are required. Based on the above discussion, the positive lower-bounding filtering (PLBF) problem is formulated as folows.

*Positive Lower-Bounding Filtering (PLBF).* Given a stable positive system (1), find a positive filter (6) with , , and that ensures that the filtering error system (8) is positive, asymptotically stable, with performance under zero initial conditions.

Similarly, the second filtering error system can be obtained by defining , , and , and it can be formulated as follows:where Then, the positive upper-bounding filtering (PUBF) problem is proposed below.

*Positive Upper-Bounding Filtering (PUBF).* Given a stable positive system (1), find a positive filter (7) with , , and that ensures that the filtering error system (10) is positive, asymptotically stable, with performance under zero initial conditions.

#### 3. Main Results

In this section, a pair of positive error-bounding filters is obtained which can bound the signal at all times with the -induced performance satisfied. Necessary and sufficient conditions in terms of linear programming are presented to design a lower-bounding filter, which is followed by parallel results obtained for the upper-bounding case.

Now, we are in a position to provide conditions to design the desired lower-bounding filter for system (1) in the following theorem.

Theorem 5. *Given a stable discrete-time positive system (1), a lower-bounding filter (6) exists such that the filtering error system (8) is positive, asymptotically stable with if and only if there exist vectors , and matrices , , such that the following LMIs are feasible:where ; . Then, a possible choice for and is given by*

*Proof . **Sufficiency.* From , (12), (13), and (20), we have that and , which together with guarantees that the lower-bounding filter (6) is a positive system.

Together with (20) and , (15)-(16) implyand we have that From (10), , , , and (14), it can be easily seen that the filtering error system (8) is positive.

Moreover, from (20), we obtain thatInequalities (17)–(19) together with (23) imply thatand we havewhere .

By Lemma 4, the filtering error system (8) is asymptotically stable with .*Necessity*. By Lemma 4, we can conclude that there exists such that inequality (25) holds. Set . We can deduce that the following inequalities hold:From and , it follows that the change of variablesmakes the problem linear and yields (17)–(20).

Next, the lower-bounding filter (6) is positive and we have , , and . Because of (27) and , and imply and .

Moreover, if the filtering error system (8) is positive, the following inequalities hold:and this implies (14)–(16). The whole proof is completed.

The parallel result is presented in the following for the upper-bounding case. We propose the following theorem to design the upper-bounding filter for positive systems. The proof is similar to the lower-bounding case and thus is omitted here.

Theorem 6. *Given a stable discrete-time positive system (1), an upper-bounding filter (7) exists such that the closed-loop system (10) is positive, asymptotically stable with if and only if there exist vectors , and matrices , , such that the following LMIs are feasible:where ; . Then, a possible choice for and is given by*

*Remark 7. *Conditions obtained in Theorems 5 and 6 can be easily solved by linear programming. Moreover, it is noted that the -induced error-bounding filters can be designed by combining Theorems 5 and 6.

#### 4. Illustrative Example

In this section, an illustrative example is given to illustrate the effectiveness of the theoretical results.

Consider system (1) with system matrices given by

For , by implementing the linear program in Theorem 5, we obtain the feasible solutionwhich yields

Moreover, for , by implementing the linear program described in Theorem 6, we obtain the feasible solutionwhich further yields

Next, the following disturbance is used in this example: Figure 1 shows the lower estimate , the output , and the upper estimate . From Figure 1, we see that the signal can be encapsulated at all times with the lower-bounding estimate and the upper-bounding one . In other words, by the filters designed through Theorem 6 and Theorem 5, the output of the original positive system can be estimated at all times.