Fault Estimation for Discrete Time-Delay Systems

and Applied Analysis 3 It is readily known that y r (k) satisfies y r (k) = d


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 timevarying 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][8][9][10][11][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].Finitehorizon  ∞ 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][15][16][17][18][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 (≥0) denotes  being a positive definite (or positive semidefinite) matrix. denotes an identity matrix of appropriate dimension.The superscripts "−1" and "" represent the inverse and transpose of a matrix.R  denotes the dimensional Euclidean space.R × is the set of all  ×  real matrices.  = 0 for  ̸ =  and   = 1.For stochastic vectors  and , inner product ⟨, ⟩ equals the covariance matrix of  and .() ∈  2 [0, ] means ∑  =0   ()() < ∞.Matrices, if the dimensions are not explicitly stated, are assumed to have compatible dimensions for algebraic operations.

Problems Statement
Consider the following linear systems with time delay: where () ∈ R  , () ∈ R  , and () ∈ R  are the state, the driving disturbance, and the fault to be estimated, respectively.Also, () ∈ R  and V() ∈ R  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 V() belong to  2 [0, ].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.

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 y  () = [y  () r  ()]

𝑇
; then the linear space generated by the measurements in the Krein space up to time  can be written as It is readily known that y  () satisfies In the sequel, we denote the Krein-space projection of m() onto L{{y  ()}  =0 } by m( | ).Construct the innovations Defining x( | ) = x() − x( | ), 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  > 0, a fault estimator r() achieving the performance (3) exists if and only if Furthermore, if the above conditions are satisfied, the desired  ∞ fault estimator is given by where ỹ() = y() − ŷ( |  − 1), 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  k  () have the same inertia.Therefore following the same line as in [1], the minimum value of   can be obtained as follows: min 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 x(− | −1) = 0 and   () can be calculated by where   () is given as 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  ( + 1,  − , ) 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  > 0, a fault estimator () that achieves the performance index (3) exists if and only if Θ() > 0 and Ξ() < 0, where Θ() and Ξ() are defined in Lemma 1.In this case, one possible finitetime  ∞ fault estimator is given by (41) 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  aug and  new 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  aug is 3, while the order of ℎ  in  new is 2. Therefore if ℎ  is large enough,  aug ≫  new .
Remark 5. Recently, by applying Krein-space approach and reorganized innovation approach, [9] has considered the finite-horizon  ∞ fault estimation for linear discrete timevarying 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.

Conclusion
The finite-time  ∞ fault estimation problem for linear timedelay 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.