Abstract
This paper investigates the finite-time fault estimation problem for linear time-delay systems, where the delay appears in both state and measurement equations. Firstly, the design of finite horizon fault estimation is converted into a minimum problem of certain quadratic form. Then we introduce a stochastic system in Krein space, and a sufficient and necessary condition for the minimum is derived by applying innovation analysis approach and projection theory. Finally, a solution to the fault estimation is obtained by recursively computing a partial difference Riccati equation, which has the same dimension as the original system. Compared with the conventional augmented approach, the solving of a high dimension Riccati equation is avoided.
1. Introduction
Krein-space theory has proven to be an effective tool in dealing with the indefinite quadratic control/filtering problems. Some recent researches on filtering have led to an interesting connection with Kalman filtering in Krein space [1, 2]. Comparing with the linear estimation approaches in Hilbert space, the Krein-space theory can lead to not only less conservative results but also computationally attractive algorithms. It has been shown in [1] that a finite horizon linear estimation problem can be cast into a problem of calculating the minimum point of a certain quadratic form. By applying linear estimation in Krein space, one can calculate recursively the minimum point via Riccati equation. In [2] the authors consider the prediction problem for time-varying continuous-time systems with delayed measurements in the finite horizon case. The necessary and sufficient condition for the existence of an predictor is obtained by applying a reorganized innovation approach in Krein space.
On the other hand, fault estimation is one of the most important issues. The paper [3] designs a fuzzy fault detection filter for T-S fuzzy systems with intermittent measurements, and all the results are formulated in the form of linear matrix inequalities. In [4], a sufficient condition for the existence of a fault filter is exploited in terms of certain linear matrix inequality. Reference [5] is concerned with the robust fault detection problem for a class of discrete-time networked systems with distributed sensors. The existence of the desired fault detection filter can be determined from the feasibility of a set of linear matrix inequalities. The paper [6] addresses the fault detection problem for discrete-time Markovian jump systems. The characterization of the gains of the desired fault detection filters is derived in terms of the solution to a convex optimization problem that can be easily solved by using the semidefinite program method. As for the fault estimation problem, the Krein-space approach has received much attention so far [7–12]. A Krein-space approach is presented in [7] to fault estimation for LDTV system, where the augmented approach [13] is also used. Different from [7], a more further result is obtained by a Krein-space approach and nonaugmented approach for the same problem in [8]. Recently, by applying Krein-space approach and reorganized innovation approach, [9] considers the finite-horizon fault estimation for linear discrete time-varying systems with delayed measurements [9]. Finite-horizon fault estimation for uncertain linear discrete time-varying systems with known inputs is considered in [10].
Recently, we note that time-delay systems have received much attention [14–19]. For optimal estimation problem, when the delay appears in state, the reorganized innovation approach is not suitable for estimation problem. Motivated by this point, we consider the finite-horizon problem fault estimation for linear discrete systems with time delay, where the delay appears in both state and measurement, which contain [9] as a special case. To the best of our knowledge, this problem has not yet been investigated, and this constitutes the primary motivation for our research. On the other hand, we desire to obtain the necessary and sufficient condition for the existence of an fault estimator. A natural idea is to use Krein space to deal with the finite-horizon fault estimation for linear time-delay systems, and this gives rise to another motivation of our work. The main contributions of the paper are highlighted as follows. (i) The necessary and sufficient condition will be derived for fault estimation problem with time delay. (ii) Compared with the augmentation approach [13], our result on estimation is given based on a partial difference Riccati equation, and hence the solving of an high dimension Riccati equation is avoided.
The organization of this paper is as follows. The problem statement is given in Section 2. Section 3 presents the fault estimator design in terms of a partial difference Riccati equation. A numerical example is given to demonstrate the effectiveness of the approach in Section 4, and the paper is concluded in Section 5.
Notation. Throughout this paper, a real symmetric matrix (≥0) denotes being a positive definite (or positive semidefinite) matrix. denotes an identity matrix of appropriate dimension. The superscripts “” and “” represent the inverse and transpose of a matrix. denotes the -dimensional Euclidean space. is the set of all real matrices. for and . For stochastic vectors and , inner product equals the covariance matrix of and . means . Matrices, if the dimensions are not explicitly stated, are assumed to have compatible dimensions for algebraic operations.
2. Problems Statement
Consider the following linear systems with time delay: where , , and are the state, the driving disturbance, and the fault to be estimated, respectively. Also, and are measurements and noises, respectively. Without loss of generality, the delays are assumed to be of an increasing order: , . Moreover, it is assumed that , , and belong to . For simplicity of presentation, we assume that , , , , and are constant matrices even though the later development and results can be easily adapted to the time-varying case.
Problem. Given the observation , seek a fault estimator () such that where is a given positive scalar, is a positive definite matrix, and Without loss of generality, the initial state estimator is assumed to be zero. The value of is assumed to be zero, where , , .
Define the fault estimation error between and as and introduce the following quadratic form: Obviously, the performance (3) is satisfied if and only if for all .
3. Main Result
We consider constructing an equivalent Krein-space problem to the minimum for . To do so we need to introduce the following stochastic systems in a Krein space: with Let ; then the linear space generated by the measurements in the Krein space up to time can be written as It is readily known that satisfies
In the sequel, we denote the Krein-space projection of onto by . Construct the innovations Defining , we further have Furthermore, we define the cross-covariance matrices of the state estimation error: By employing the Krein-space theory, a necessary and sufficient condition for the existence of the desired fault estimator is given in the following.
Lemma 1. Consider the stochastic systems (7)–(9). For a given , a fault estimator achieving the performance (3) exists if and only if Furthermore, if the above conditions are satisfied, the desired fault estimator is given by where , and the minimum of the quadratic form is
Proof. It can be seen from (12) and (15) that Equation (20) can be further written as where We can draw the conclusion from (17) and (21) that and have the same inertia. Therefore following the same line as in [1], the minimum value of can be obtained as follows: Thus the rest of the proof is clear.
In the following, we are devoted to the estimator design in terms of the solution to a partial difference Riccati equation.
Lemma 2. The state estimate is recursively calculated as where the initial values are and can be calculated by while is calculated by the partial difference Riccati equation as
Proof. Applying projection theory, we have
where is given as
Noting that , then based on (15) one has
Then one has
Next it follows from (7) and (32) that
Based on (36), we have the estimation error covariance matrices
which is (26). Combining (7) and (31), we obtain
Furthermore, one has
Finally (29) is straightforward by virtue of the definition of (16). Thus the proof is completed here.
Theorem 3. Consider the system (1)-(2). For a given , a fault estimator that achieves the performance index (3) exists if and only if and , where and are defined in Lemma 1. In this case, one possible finite-time fault estimator is given by where
Proof. According to [1], we can see that the fault estimation problem addressed for the deterministic systems (1) with (2) and (5) is partially equivalent to that for the stochastic systems (7) with (8) and (9) in a Krein space, and therefore the proof is readly and we omitted here.
Remark 4. In fact, the problem mentioned in this paper can be converted into the problem in [7, 8] by applying augmented approach. However, due to the existence of time delay, we need to solve a high dimension Riccati equation. Here the solutions to the fault estimator can be obtained by solving partial difference Riccati equations (26)–(28), which have the same dimension as the original system (1). Therefore solving a high dimension Riccati equation is avoided. Here we present a simple explanation. Because the multiplications and divisions cost much more in computation than additions, hence we only use the number of multiplications and divisions as the operation count. Denote and as the numbers of multiplications and divisions for augmentation method and our proposed approach in one step, respectively. According to [18], one can see that the order of in is 3, while the order of in is 2. Therefore if is large enough, .
Remark 5. Recently, by applying Krein-space approach and reorganized innovation approach, [9] has considered the finite-horizon fault estimation for linear discrete time-varying systems with two-channel single measurement delay. In this paper, we have investigated the finite-horizon problem fault estimation for linear discrete systems with time delay, where the delay appears in both state and measurement, which contain [9] as a special case.
4. Numerical Example
Consider the linear discrete-time system:
The finite time horizon concerned here is . The driving disturbance and measurement noise are selected as , . The fault to be estimated is assumed to be Set initial values as , , , , , . By using the result given in Theorem 3, the desired fault estimator is designed with . Figure 1 shows the fault and its estimate, which confirm that the designed estimator performs very well.
5. Conclusion
The finite-time fault estimation problem for linear time-delay systems has been investigated. The design of finite horizon fault estimation has been converted into a minimum problem of certain quadratic form. Then an stochastic system in Krein space has been proposed, and a sufficient and necessary condition for the minimum has been derived by applying innovation analysis approach and projection theory. Finally a solution to the fault estimation has been obtained by recursively computing a partial difference Riccati equation. Compared with the conventional augmented approach, the presented approach lessens the computational demand when the delay is large. In the further study, we will consider the fault estimation problem for linear time-delay systems with multiplicative noise.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported in part by the National Natural Science Foundation of China (61104126, 61304013, and 61170145), the Doctoral Fund of Ministry of Education of China (20113704120005), and the Excellent Young Scholars Research Fund of Shandong Normal University.