Nothing Special   »   [go: up one dir, main page]

Next Article in Journal
Robust Filtering Techniques for RTK Positioning in Harsh Propagation Environments
Next Article in Special Issue
A Robust Steered Response Power Localization Method for Wireless Acoustic Sensor Networks in an Outdoor Environment
Previous Article in Journal
FOD Detection Method Based on Iterative Adaptive Approach for Millimeter-Wave Radar
Previous Article in Special Issue
A Two-Phase Distributed Filtering Algorithm for Networked Uncertain Systems with Fading Measurements under Deception Attacks
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Event-Triggering State and Fault Estimation for a Class of Nonlinear Systems Subject to Sensor Saturations

1
College of Information Science and Technology, Donghua University, Shanghai 201620, China
2
Engineering Research Center of Digitalized Textile and Fashion Technology, Ministry of Education, Shanghai 201620, China
3
College of Electrical Engineering and Automation, Shandong University of Science and Technology, Qingdao 266590, China
4
Artificial Intelligence Energy Research Institute, Northeast Petroleum University, Daqing 163318, China
5
Heilongjiang Provincial Key Laboratory of Networking and Intelligent Control, Northeast Petroleum University, Daqing 163318, China
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(4), 1242; https://doi.org/10.3390/s21041242
Submission received: 27 December 2020 / Revised: 4 February 2021 / Accepted: 5 February 2021 / Published: 10 February 2021
(This article belongs to the Special Issue Sensor Fusion and Signal Processing)
Figure 1
<p>Structure of the state and fault estimation.</p> ">
Figure 2
<p>Trajectories <math display="inline"><semantics> <msubsup> <mi>x</mi> <mi>k</mi> <mn>1</mn> </msubsup> </semantics></math>, <math display="inline"><semantics> <msubsup> <mi>x</mi> <mi>k</mi> <mn>2</mn> </msubsup> </semantics></math> and their estimates.</p> ">
Figure 3
<p>Trajectories of fault <math display="inline"><semantics> <msub> <mi>f</mi> <mi>k</mi> </msub> </semantics></math> and its estimate.</p> ">
Figure 4
<p><math display="inline"><semantics> <mrow> <mrow> <mi>Log</mi> <mo>(</mo> </mrow> <msubsup> <mi>MSE</mi> <mi>k</mi> <mi>x</mi> </msubsup> <mo>)</mo> </mrow> </semantics></math> of <math display="inline"><semantics> <msubsup> <mi>x</mi> <mi>k</mi> <mn>1</mn> </msubsup> </semantics></math>, <math display="inline"><semantics> <msubsup> <mi>x</mi> <mi>k</mi> <mn>2</mn> </msubsup> </semantics></math> and the traces of their upper bounds.</p> ">
Figure 5
<p><math display="inline"><semantics> <mrow> <mrow> <mi>Log</mi> <mo>(</mo> </mrow> <msubsup> <mi>MSE</mi> <mi>k</mi> <mi>f</mi> </msubsup> <mo stretchy="false">)</mo> </mrow> </semantics></math> of <math display="inline"><semantics> <msub> <mi>f</mi> <mi>k</mi> </msub> </semantics></math> and the trace of its upper bound.</p> ">
Figure 6
<p>The position of the target <math display="inline"><semantics> <msub> <mi>x</mi> <mrow> <mn>1</mn> <mo>,</mo> <mi>k</mi> </mrow> </msub> </semantics></math> and its estimate.</p> ">
Figure 7
<p>The velocity of the target <math display="inline"><semantics> <msub> <mover accent="true"> <mi>x</mi> <mo>˙</mo> </mover> <mrow> <mn>1</mn> <mo>,</mo> <mi>k</mi> </mrow> </msub> </semantics></math> and its estimate.</p> ">
Figure 8
<p>The position of the target <math display="inline"><semantics> <msub> <mi>x</mi> <mrow> <mn>2</mn> <mo>,</mo> <mi>k</mi> </mrow> </msub> </semantics></math> and its estimate.</p> ">
Figure 9
<p>The velocity of the target <math display="inline"><semantics> <msub> <mover accent="true"> <mi>x</mi> <mo>˙</mo> </mover> <mrow> <mn>2</mn> <mo>,</mo> <mi>k</mi> </mrow> </msub> </semantics></math> and its estimate.</p> ">
Figure 10
<p>The actual fault and its estimate.</p> ">
Versions Notes

Abstract

:
This paper is concerned with the state and fault estimation issue for nonlinear systems with sensor saturations and fault signals. For the sake of avoiding the communication burden, an event-triggering protocol is utilized to govern the transmission frequency of the measurements from the sensor to its corresponding recursive estimator. Under the event-triggering mechanism (ETM), the current transmission is released only when the relative error of measurements is bigger than a prescribed threshold. The objective of this paper is to design an event-triggering recursive state and fault estimator such that the estimation error covariances for the state and fault are both guaranteed with upper bounds and subsequently derive the gain matrices minimizing such upper bounds, relying on the solutions to a set of difference equations. Finally, two experimental examples are given to validate the effectiveness of the designed algorithm.

1. Introduction

State estimation/filtering problems have always been one of the fundamental issues in the areas of target tracking, navigation and positioning, electric power systems, econometrics, biosystems, etc. Therefore, enormous research attention has been paid to the state estimation problems and some elegant work has been reported, see e.g., [1,2,3,4,5]. According to different performance indices, the current state estimation approaches include Kalman filtering (KF), extend Kalman filtering (EKF), H filtering and so on. To be specific, the famous KF approach has been proposed in [2] by providing optimal state estimates in the sense of minimal mean-squared error under the assumption that system parameters and noise statistics are precisely known. The H filtering method proposed [6,7] is capable to attenuate the influence from the exogenous disturbance to the filtering error. When it comes to the case that the system model is nonlinear or uncertain, the celebrated EKF approach has been shown to be a useful tool for the state estimation issues. For instance, in [8], the EKF approach has been developed to cope with the nonlinear systems subject to missing measurement. Moreover, in [9], the filtering approach has also been applied in the complex networks with incomplete measurements.
It is often the case that the faults are inevitable in practical applications because of a variety of reasons including component failure, ageing equipment, complex external environment, bandwidth limitation, etc. During the past few decades, fault detection and fault-tolerant control issues have gained considerable research enthusiasm due to the demand in reliability and security of the practical systems. It should be pointed out that the accurate information of the fault signals are hard to acquire, whilst the recently emerging state and fault estimation (SFE) approach provides a good solution to obtain the sufficient information of the state and fault signals simultaneously. By such a merit, increasing research attention has been paid on this aspect recently, and some inspiring work has been available in the literature [10,11,12]. For instance, an SFE algorithm relying on a recursive approach has been designed in [13] for uncertain systems with missing measurements and stochastic nonlinearities. H SFE problems have also been studied for various dynamic systems, such as fuzzy systems [12,14], nonlinear systems [15,16] and 2-D systems [17,18]. Nevertheless, the SFE problems have not been thoroughly investigated yet and still have been a research hotspot in control/filtering community.
In reality, sensors may not always provide signals with unlimited amplitudes owing to the physical constraints. If the sensor saturation is not properly handled, it will severely decrease the system performance. The main challenge of research on the saturations in control community is how to design a filtering/control algorithm that can effectively dealing with the nonlinearities brought by the sensor saturations. As a consequence, the filtering/control issues subject to sensor saturations have gained initial research focus, see, e.g., [19,20]. For example, a recursive filtering issue has been solved for uncertain systems with faults and sensor saturations in [21]. In [22,23], the H filtering issues have been settled for nonlinear systems with incomplete measurements and sensor saturations. In [20], mean-squared consensus control problem has been studied for stochastic multi-agent systems subject to sensor saturations where the desired controllers have been designed depending on the solutions to recursive matrix inequalities.
On another research front, the event-triggering mechanism (ETM) has become a research hotspot recently due mainly to its superiority of effectively reducing communication resources compared with the traditional time-triggering protocol. Under the ETM the current measurement will be transmitted only when the predefined triggering condition is met, and thereby the transmission numbers can be reduced largely. Based on this idea, various control and filtering issues under the ETM have been studied, see e.g., [6,24,25,26]. Very recently, considerable research attention has been paid on the event-triggering fault estimation issue owing to its vital role in the practical engineering. Accordingly, the event-triggering fault estimation problems have been investigated for various systems, such as nonlinear systems with missing measurements [27], stochastic systems subject to nonlinearities and packet dropouts [28], and stochastic systems with deception attacks [29]. However, to the best of the authors’ knowledge, the event-triggering state and fault estimation (ETSFE) problem for nonlinear systems with sensor saturations has not been fully studied, which constitutes the main motivation of this paper.
In terms of the methodologies, due to the effects brought by nonlinearities (including the saturation functions) and event-triggering protocol, it is almost impossible in the Kalman filtering framework to minimize the estimation error covariance through adjustment of the gain matrices. In [30], an alternative way has been proposed to handle the effects of norm bounded parameter uncertainties and a robust filter has been designed such that an upper bounded matrix of the estimation error covariance is minimized. Enlightened by this idea, such a filtering approach has been applied in various complex systems such as complex networks [9] and sensor networks [31]. However, it should be pointed out that, for the state and fault simultaneous estimation problem in the existing literature, the estimation error covariance minimization method has still been the main method which is incapable of dealing with more real complex phenomena. Therefore, it is the second motivation of this paper to develop the filtering approach proposed in [30] to handle the state and fault estimation problems with sensor saturations under the event-triggering strategy.
The novelties of this paper are emphasized as follows: (1) a novel ETSFE issue is, for the first time, addressed when the effects of sensor saturations, nonlinearities as well as ETM are simultaneously taken into consideration; (2) the state and fault estimator is designed such that the upper bounds on the error covariances of the state and fault estimation are respectively guaranteed at each time instant; and (3) the gain matrices are designed via two recursions which minimize the obtained upper bounds. Finally, two illustrative examples are utilized to verify the feasibility of the developed ETSFE algorithm.
The remaining part of this paper is organized as follows. In Section 2, the problem to be investigated is addressed. The main results are listed in Section 3 where the desired state and fault estimators is designed. In Section 4, two illustrative examples are given and the conclusion is drawn in Section 5.
Notations: In this paper, the notations mentioned are standard. R n and R m × n respectively denote the n-dimensional Euclidean space and m × n real matrix. I is the identity matrix, while diag { a 1 , a 2 , , a N } represents the block-diagonal matrix with matrices a 1 , a 2 , , a N . For symmetric matrices x and y, x y ( x > y ) means that x y is positive semi-define (positive definite) matrix. The superscript “T” and “ 1 ” refer to matrix transposition and inverse, respectively. R ( M ) is the rank of the matrix M. E { x } denotes the mathematical expectation of the stochastic variable x. tr { M } denotes the trace of the matrix M. · stands for the Euclidean norm.

2. Problem Formulation

The estimation structure under consideration is shown in Figure 1 and the dynamics of the plant is given by
x k + 1 = h k ( x k ) + B k f k + w k
where x k R n x is the system state vector, f k R n f represents a fault signal, w k R n x is the process noise, and B k is a given compatible matrix.
The measurements with sensor saturation are described by
y k = ϑ ( C k x k ) + D k f k + v k
where y k R n y represents the measurement vector at time instant k, v k R n y is the measurement noise, and C k , D k are both appropriate-dimensional matrices.
Assumption 1.
The matrix D k is full column rank, i.e., R ( D k ) = n f , n f n y .
The nonlinear function h k ( · ) satisfies the following condition
h k ( X ) h k ( Y ) ν k X Y , X , Y R n x
where ν k > 0 is a known matrix.
The noise signals w k and v k have the following statistical properties
E { w k } = 0 , E { w k w l T } = R k δ k l , E { v k } = 0 , E { v k v l T } = Q k δ k l
where R k > 0 and Q k > 0 are known appropriate-dimensional matrices, and δ k l represents the Kronecker function with
δ k l = 1 , k = l , 0 , k l .
The saturation function ϑ ( · ) : R n y R n y is defined by
ϑ ( s ) = ϑ ( s 1 ) ϑ ( s 2 ) ϑ ( s n y ) T
with ϑ ( s i ) = sign ( s i ) min { ϱ i , | s i | } ( i = 1 , 2 , , n y ) , where s i is the ith element of vector s, sign ( · ) represents the signum function, and ϱ i denotes the saturation level for ith element.
For the sake of reducing limited communication resource, the ETM is adopted to govern the transmission frequency between the sensor and the estimator. We denote the transmission instants by 0 = k 0 < k 1 < k 2 < < k l < , which is determined by
k l + 1 = min k N | k > k l , y k y k l > τ
where τ > 0 is a given scalar and y k l is the measurement transmitted at the latest time.
For the purpose of estimating the state and fault simultaneously, we construct the estimator as follows
x ^ k + 1 = h k ( x ^ k ) + B k f ^ k + G k y k l C k x ^ k D k f ^ k f ^ k = L k y k l C k x ^ k
where f ^ k , x ^ k represent the estimates of fault and state respectively and L k , G k are the estimator parameters respectively.
Let the state estimation error and fault estimation error be x ˜ k = x k x ^ k and f ˜ k = f k f ^ k , respectively.
By noting (1) and (7), one has
x ˜ k + 1 = h k ( x k ) h k ( x ^ k ) + ( B k G k D k ) f ˜ k + G k ε k G k ϑ ( C k x k ) G k v k + G k C k x ^ k + w k
and
f ˜ k = f k L k D k f k L k ( ϑ ( C k x k ) + v k ε k C k x ^ k )
where ε k = y k y k l .
Assuming that the constraint condition L k D k = I is met, we eventually derive
f ˜ k = L k ( ϑ ( C k x k ) + v k ε k C k x ^ k ) .
Remark 1.
The constraint condition L k D k = I plays a key role in the estimator design. It is obvious that the fault estimator (9) contains the fault vectors f k . Since the dynamics of the faults is generally unknown, the fault term should be eliminated. To this end, L k D k = I is introduced as an additional condition.
Then, we define the estimation error covariances of the state and fault as follows
P k x = E { x ˜ k x ˜ k T } , P k f = E { f ˜ k f ˜ k T } .
Our main objective of this paper is to develop an event-triggering state and fault estimator of the form (7) such that, for all nonlinearities as well as sensor saturations, the upper bounds ( Σ k x and Σ k f ) for the estimation error covariances of state and fault are respectively guaranteed, that is
P k x = E { x ˜ k x ˜ k T } Σ k x , P k f = E { f ˜ k f ˜ k T } Σ k f .
Moreover, the designed gain matrices G k and L k are expected to minimize the upper bound Σ k x and Σ k f simultaneously at each iteration.

3. Main Results

In this section, the upper bounds on the estimation error covariances of the state and fault are expressed by means of recursions. Then, the proper gain matrices G k and L k are designed to minimize the upper bounds on the estimation error covariances and fault error covariances, respectively. The following lemmas will be used for obtaining the results.
Lemma 1 ([4]).
For k [ 0 , N ] , let the matrix function be Ξ k ( · ) : R n × n R n × n , and arbitrary symmetric matrices x > 0 and y > 0 . If Ξ k ( x ) Ξ k ( y ) for all x y , then under the initial condition G 0 = H 0 , the solutions G k and H k to difference equations G k + 1 = Ξ k ( G k ) , and H k + 1 = Ξ k ( H k ) satisfy G k + 1 H k + 1 .
Lemma 2 ([8]).
The following relationship is true for arbitrary real vectors M and N
M N T + M N T ϵ M M T + ϵ 1 N N T
where ϵ > 0 is an arbitrary scalar.
Lemma 3.
Under the constraint condition L k D k = I , the fault error covariance P k f = E { f ˜ k f ˜ k T } can be derived as follows
P k f = L k Λ k L k T
where
Λ k = E ϑ ( C k x k ) ϑ T ( C k x k ) + C k x ^ k x ^ k T C k T + ε k ε k T + ε k x ^ k T C k T + C k x ^ k ε k T ϑ ( C k x k ) ε k T ε k ϑ T ( C k x k ) C k x ^ k ϑ T ( C k x k ) + Q k ϑ ( C k x k ) x ^ k T C k T ε k v k T v k ε k T .
Proof. 
The validation of (13) can be verified by noting (10) and (11) and the rest of proof is omitted. □
Lemma 4.
Under the condition L k D k = I , the state estimation error covariance P k x = E { x ˜ k x ˜ k T } is derived by
P k + 1 x = E ( h k ( x k ) h k ( x ^ k ) ) ( h k ( x k ) h k ( x ^ k ) ) T + ( B k G k D k ) f ˜ k f ˜ k T ( B k G k D k ) T + G k ϑ ( C k x k ) ϑ T ( C k x k ) G k T + G k ε k ε k T G k T + G k C k x ^ k x ^ k T C k T G k T + R k + G k Q k G k T k , 1 k , 1 T + k , 2 + k , 2 T + k , 3 + k , 3 T + k , 4 + k , 4 T + k , 5 + k , 5 T + k , 6 + k , 6 T k , 7 k , 7 T k , 8 k , 8 T k , 9 k , 9 T + k , 10 + k , 10 T k , 11 k , 11 T
where
k , 1 = E ( h k ( x k ) h k ( x ^ k ) ) ϑ T ( C k x k ) G k T , k , 2 = E ( h k ( x k ) h k ( x ^ k ) ) f ˜ k T ( B k G k D k ) T , k , 3 = E ( h k ( x k ) h k ( x ^ k ) ) x ^ k T C k T G k T , k , 4 = E ( h k ( x k ) h k ( x ^ k ) ) ε k T G k T , k , 5 = ( B k G k D k ) E f ˜ k x ^ k T C k T G k T , k , 6 = ( B k G k D k ) E f ˜ k ε k T G k T , k , 7 = ( B k G k D k ) E f ˜ k ϑ T ( C k x k ) G k T , k , 8 = G k E ϑ ( C k x k ) x ^ k T C k T G k T , k , 9 = G k E ϑ ( C k x k ) ε k T G k T , k , 10 = G k C k E x ^ k ε k T G k T , k , 11 = G k v k ε k T G k T .
Proof. 
It can be shown that (14) follows directly from (8) and (11), and the proof is omitted for conciseness. □

3.1. Fault Estimation

The following Theorem 1 provides the explicit form of the upper bound on the fault estimation error covariance P k f in terms of the recursion.
Theorem 1.
Consider the fault estimation error covariance in (13). Assume that the condition L k D k = I is satisfied. For any given positive scalars a k , b k , c k and d k , the upper bound on the fault estimation error covariance P k f is obtained by
Σ k f = L k Λ ¯ k L k T
where
Λ ¯ k = κ 1 τ 2 I + κ 2 i = 1 n y ϱ i 2 I + κ 3 C k x ^ k x ^ k T C k T + κ 4 Q k , κ 1 = 1 + a k 1 + b k + d k , κ 2 = 1 + a k + c k κ 3 = 1 + b k 1 + c k 1 , κ 4 = 1 + d k 1
Proof. 
In view of the triggering condition (6), one has
ε k ε k T ε k T ε k I τ 2 I .
Using Lemma 2, we obtain
ε k ϑ T ( C k x k ) ϑ ( C k x k ) ε k T a k ϑ ( C k x k ) ϑ T ( C k x k ) + a k 1 ε k ε k T ,
ε k x ^ k T C k T + C k x ^ k ε k T b k ε k ε k T + b k 1 C k x ^ k x ^ k T C k T ,
ϑ ( C k x k ) x ^ k T C k T C k x ^ k ϑ T ( C k x k ) c k ϑ ( C k x k ) ϑ T ( C k x k ) + c k 1 C k x ^ k x ^ k T C k T ,
and
ε k v k T v k ε k T d k ε k ε k T + d k 1 v k v k T .
Moreover, from the definition of the saturation function, we have
E { ϑ ( C k x k ) ϑ T ( C k x k ) } i = 1 n y ϱ i 2 I .
It then follows from (13), (17)–(22) that
Λ k Λ ¯ k .
Finally, considering (13), (15) and (23), we have P k f Σ k f , which ends this proof. □
By the results obtained in Theorem 1, the following theorem is going to design a gain matrix L k such that the upper bound on the fault estimation error variance is minimized at each iteration.
Theorem 2.
Under the constraint condition L k D k = I and supposing that a k , b k , c k and d k are given positive scalars, the upper bound Σ k f on the fault estimation error covariance is minimized, if the estimator gain L k is chosen as
L k * = ( D k T Λ ¯ k 1 D k ) 1 D k T Λ ¯ k 1 .
Meanwhile, the minimum upper bound is given by
Σ k f * = ( D k T Λ ¯ k 1 D k ) 1 .
Proof. 
This proof is substantially to solve the following constrained optimization problem
min L k { Σ k f } s . t . L k D k = I .
By means of Lagrange multiplier method, we introduce the following Lagrange function
Ξ ( L k , Υ k ) = L k Λ ¯ k L k T + ( I L k D k ) Υ k + Υ k T ( I L k D k ) T
where Υ k is the Lagrange factor.
The derivatives of Ξ ( L k , Υ k ) with respect to L k and Υ k can be written as
Ξ ( L k , Υ k ) L k = 2 L k Λ ¯ k 2 Υ k T D k T
and
Ξ ( L k , Υ k ) Υ k = 2 ( I L k D k ) .
Letting the above derivatives be zero, we have
L k * = Υ k T D k T Λ ¯ k 1
and
L k * D k = I .
Then, substituting (28) into (29), one has
Υ k = ( D k T Λ ¯ k 1 D k ) 1
from which we have L k * = ( D k T Λ ¯ k 1 D k ) 1 D k T Λ ¯ k 1 .
Combining (27), (28) and (30), we derive the minimized upper bound on the fault estimation error covariance as follows
Σ k f * = Ξ ( L k * , Υ k ) = ( D k T Λ ¯ k 1 D k ) 1 .
The proof is now complete. □

3.2. State Estimation

In the following theorem, an upper bound on the state estimation error covariance is derived by means of the recursion and then minimized by the designed gain parameter G k * .
Theorem 3.
Consider the state estimation error covariance obtained in (14). Assume that the condition L k D k = I is satisfied and π i > 0 ( i = 1 , 2 , , 11 ) are arbitrary positive scalars. If matrix Σ k satisfies the following difference equation
Σ k + 1 x = ϕ 1 ν k ν k T tr ( Σ k x ) + R k + ϕ 2 ( B k G k D k ) Σ k f ( B k G k D k ) T + ϕ 3 τ 2 G k G k T + ϕ 4 i = 1 n y ϱ i 2 G k G k T + ϕ 6 G k Q k G k T + ϕ 5 G k C k x ^ k x ^ k T C k T G k T ,
with the initial value P 0 x Σ 0 x , where
ϕ 1 = 1 + π 1 + π 2 + π 3 + π 4 , ϕ 2 = 1 + π 2 1 + π 5 + π 6 + π 7 , ϕ 3 = 1 + π 4 1 + π 6 1 + π 9 1 + π 10 1 + π 11 1 , ϕ 4 = 1 + π 1 1 + π 7 1 + π 8 + π 9 , ϕ 5 = 1 + π 3 1 + π 5 1 + π 8 1 + π 10 , ϕ 6 = 1 + π 11 ,
then Σ k + 1 x is the upper bound of P k + 1 x , i.e., P k + 1 x Σ k + 1 x .
Moreover, if the gain matrix G k is selected by
G k * = Φ k Ψ k 1
where
Φ k = ϕ 2 B k Σ k f D k T
and
Ψ k = ϕ 2 D k Σ k f D k T + ϕ 3 τ 2 I + ϕ 4 i = 1 n y ϱ i 2 I + ϕ 5 C k x ^ k x ^ k T C k T + ϕ 6 Q k ,
then the upper bound Σ k + 1 x is minimized and the minimum upper bound is given by
Σ k + 1 x * = ϕ 1 ν k ν k T tr ( Σ k x ) + ϕ 2 B k Σ k f B k T + R k ϕ 2 B k Σ k f D k T ϕ 2 D k Σ k f D k T + ϕ 3 τ 2 I + ϕ 4 i = 1 n y ϱ i 2 I + ϕ 5 C k x ^ k x ^ k T C k T + ϕ 6 Q k 1 ϕ 2 B k Σ k f D k T ϕ 2 D k Σ k f D k T + ϕ 3 τ 2 I + ϕ 4 i = 1 n y ϱ i 2 I + ϕ 5 C k x ^ k x ^ k T C k T + ϕ 6 Q k 1 T .
Proof. 
By noting (14) and using Lemma 2, one has
k , 1 k , 1 T π 1 E { ( h k ( x k ) h k ( x ^ k ) ) ( h k ( x k ) h k ( x ^ k ) ) T } + π 1 1 G k E { ϑ ( C k x k ) ϑ T ( C k x k ) } G k T , k , 2 + k , 2 T π 2 E { ( h k ( x k ) h k ( x ^ k ) ) ( h k ( x k ) h k ( x ^ k ) ) T } + π 2 1 ( B k G k D k ) P k f ( B k G k D k ) T , k , 3 + k , 3 T π 3 E { ( h k ( x k ) h k ( x ^ k ) ) ( h k ( x k ) h k ( x ^ k ) ) T } + π 3 1 G k C k x ^ k x ^ k T C k T G k T , k , 4 + k , 4 T π 4 E { ( h k ( x k ) h k ( x ^ k ) ) ( h k ( x k ) h k ( x ^ k ) ) T } + π 4 1 G k E ε k ε k T G k T , k , 5 + k , 5 T π 5 ( B k G k D k ) P k f ( B k G k D k ) T + π 5 1 G k C k x ^ k x ^ k T C k T G k T , k , 6 + k , 6 T π 6 ( B k G k D k ) P k f ( B k G k D k ) T + π 6 1 G k E ε k ε k T G k T , k , 7 k , 7 T π 7 ( B k G k D k ) P k f ( B k G k D k ) T + π 7 1 G k E { ϑ ( C k x k ) ϑ T ( C k x k ) } G k T , k , 8 k , 8 T π 8 G k E { ϑ ( C k x k ) ϑ T ( C k x k ) } G k T + π 8 1 G k C k x ^ k x ^ k T C k T G k T , k , 9 k , 9 T π 9 G k E { ϑ ( C k x k ) ϑ T ( C k x k ) } G k T + π 9 1 G k E ε k ε k T G k T , k , 10 + k , 10 T π 10 G k C k x ^ k x ^ k T C k T G k T + π 10 1 G k E ε k ε k T G k T , k , 11 k , 11 T π 11 G k Q k G k T + π 11 1 G k E ε k ε k T G k T .
Combining (14), (33) and (36), we have
P k + 1 x ϕ 1 E { ( h k ( x k ) h k ( x ^ k ) ) ( h k ( x k ) h k ( x ^ k ) ) T } + ϕ 2 ( B k G k D k ) P k f ( B k G k D k ) T + ϕ 3 G k E ε k ε k T G k T + ϕ 4 G k E { ϑ ( C k x k ) ϑ T ( C k x k ) } G k T + ϕ 5 G k C k x ^ k x ^ k T C k T G k T + ϕ 6 G k Q k G k T + R k .
From (3), (17), (22) and (37), we obtain
P k + 1 x ϕ 1 ν k ν k T tr ( P k x ) + R k + ϕ 2 ( B k G k D k ) P k f ( B k G k D k ) T + ϕ 3 τ 2 G k G k T + ϕ 4 i = 1 n y ϱ i 2 G k G k T + ϕ 5 G k C k x ^ k x ^ k T C k T G k T + ϕ 6 G k Q k G k T .
By using Lemma 1, we arrive at P k + 1 x Σ k + 1 x .
Next, the gain parameter G k given by (34) is ready to be optimal in the sense of minimizing the upper bound Σ k + 1 x .
The derivative of Σ k + 1 x with respect to G k is computed by
G k Σ k + 1 x = 2 ϕ 2 ( B k G k D k ) Σ k f D k T + 2 ϕ 3 τ 2 G k + 2 ϕ 4 i = 1 n y ϱ i 2 G k + 2 ϕ 5 G k C k x ^ k x ^ k T C k T + 2 ϕ 6 G k Q k .
Letting the derivative in (39) be zero, one has
G k * = ϕ 2 B k Σ k f D k T ϕ 2 D k Σ k f D k T + ϕ 3 τ 2 I + ϕ 4 i = 1 n y ϱ i 2 I + ϕ 5 C k x ^ k x ^ k T C k T + ϕ 6 Q k 1 ,
and the minimum upper bound is given as the form as (35). The proof is now complete. □
Based on the above results, the developed ETSFE algorithm is summarized as follows (Algorithm 1).
Algorithm 1: ETSFE algorithm
1. Let parameters a k , b k , c k , d k , π i ( i = 1 , 2 , , 11 ) be given. Set initial values x ^ 0 = x ¯ 0 and Σ 0 x = P 0 x , the length of time horizon N and k = 0 ;
2. Calculate the fault estimator gain matrix L k * according to (24), the upper bound of the fault estimation error covariance Σ k f * via (25), and the fault estimate f ^ k according to (7);
3. Calculate the state estimator gain matrix G k * according to (34), the upper bound of the state estimation error covariance Σ k x * via (35), and the state estimate x ^ k according to (7);
4. If k N , set k = k + 1 and go to step 2, else go to step 5;
5. Stop.
Remark 2.
Theorems 2 and 3 provide the expressions for the estimator gain matrices L k and G k , respectively, and the desired estimator has been designed. In the design of the estimator, three difficulties can be identified as follows: (1) how can we choose an appropriate standard to evaluate the estimation performance when the systems are subject to nonlinearities and sensor saturations under the event-triggering protocol? (2) how can we construct an appropriate estimator structure which is closely related to the selection of the estimation analysis method? (3) After the estimator structure is determined, how can we obtain the expressions for the estimator gain matrices by using the mathematical tool appropriately. In this paper, these questions have been well answered.
Remark 3.
To date, we have studied the ETSFE problems for nonlinear systems subject to sensor saturations. Relying on the matrix analysis technique, the estimation error covariances for the state and fault are both guaranteed by upper bounds, and then such upper bounds have been minimized by appropriately designing the gain parameters L k * and G k * . Comparing with the existing methods in literature [10,21], the estimation approach proposed in this paper is capable of dealing with state and fault simultaneous estimation problem when the systems are subject to the nonlinearities including the sensor saturation nonlinearity under the event-triggering protocols. Actually, the proposed estimation approach is applicable to all those situations where the complex phenomena could be modelled by a bound-limited variable. Moreover, the corresponding ETSFE algorithm proposed is of a simple and recursive form, which is suitable for online computation. The following section will provide two illustrative examples to validate the usefulness of the developed ETSFE algorithm.

4. Experimental Simulation

In this section, two simulation examples are utilized to validate the usefulness of the developed ETSFE algorithm.
Example 1.
The system under consideration in (1) has the following parameters
A k = 0.6 0.5 0.2 0.4 + 0.4 cos ( k ) , B k = [ 0.9 0.7 ] T , C k = 0.4 1 0.3 1.5 , D k = [ 0.9 0.9 ] T , ν k = 0.15 .
Let the nonlinear function h k ( x k ) be
h k ( x k ) = A k x k + h ¯ k ( x k )
where h ¯ k ( x k ) = 0.17 sin ( x k 1 ) 0.16 sin ( x k 2 ) T .
The process noise w k and measurement noise v k are Gaussian noises with zero mean and their covariances are set as R k = 0.3 0 0 0.3 and Q k = 0.2 0 0 0.2 , respectively. The saturation levels are ϱ 1 = ϱ 2 = 0.1 and the threshold of the triggering threshold is τ = 1.2 . In this example, the parameters π i ( i = 1 , 2 , , 11 ) , a k , b k , c k and d k are chosen as 1, and the initial state is chosen as x 0 = [ 0.5 0.5 ] T . Based on the above parameters, the gain matrices L k * and G k * can be iteratively computed by (24) and (34).
The designed algorithm is validated in MATLAB (R2016a), Intel Core CPU i5-8265. Figure 2, Figure 3, Figure 4 and Figure 5 show the simulation results based on the parameters mentioned above. The above picture in Figure 2 plots the state trajectory and its estimate for x k 1 and the picture below shows the state trajectory and the estimate for x k 2 . The fault signals as well as its estimates are shown in Figure 3. It is seen from Figure 2 and Figure 3 that the proposed estimator performs well. The mean square errors ( MSE k ) and the traces of their minimal upper bounds for state are shown in Figure 4. The corresponding results for fault are given in Figure 5. The simulation results shown in Figure 4 and Figure 5 concur with our theoretical analysis.
Remark 4.
Due to the effects of the nonlinearities and the event-triggering protocol, upper bounds of the estimation error covariances are employed as the performance index of state and fault estimation. Such a practice inevitably gives rise to conservatism which may affect the real estimation accuracy. Although the estimator gains are designed to minimize the upper bounds at each step, the minimized upper bounds are actually not tight. Note from Theorems 2 and 3 that the minimized upper bounds are closely related to the parameters a k , b k , c k , d k , π i > 0 ( i = 1 , 2 , , 11 ) . Therefore, in the experiment, these parameters should be selected prudently and the appropriate selection of these parameters may further improve the experimental results.
Example 2.
In this experiment, a ballistic object tracking example is employed to further validate the feasibility of the proposed ETSFE algorithm. The ballistic object tracking system is described by [21,32]
x k + 1 = h k ( x k ) + B k f k + w k y k = ϑ ( C k x k ) + D k f k + v k
where h k ( x k ) = A k x k + h ¯ k ( x k ) , h ¯ k ( x k ) = ν k ( g k ( x k ) + G ) , g k ( x k ) = g ψ ( x 2 , k ) 2 Υ x ˙ 1 , k 2 + x ˙ 2 , k 2 x ˙ 1 , k x ˙ 2 , k , ψ ( x 2 , k ) = κ 1 · exp ( κ 2 x 2 , k ) and
A k = 1 T 0 0 0 1 0 0 0 0 1 T 0 0 0 1 , B k = 0.01 0 0.01 0 , C k = 1 0 0 0 0 0 1 0 , D k = 0.01 0 0.01 0 , G = 0 g , ν k = T 2 2 0 T 0 0 T 2 2 0 T .
The parameters involved in the above system are stated as follows: x k = [ x 1 , k x ˙ 1 , k x 2 , k x ˙ 2 , k ] T where x 1 , k , x ˙ 1 , k represent respectively the position and the velocity of the target on x-axis and x 2 , k , x ˙ 2 , k represent the position and the velocity of the target on y-axis; T is the sampling period; g denotes the gravity acceleration; Υ represents the ballistic coefficient; the exponentially decaying function ψ ( · ) is the air density; and w k , v k are both Gaussian noises with zero mean and their covariances are R k = λ · diag { R ¯ k , R ¯ k } and Q k = 10 I 2 where R ¯ k = T 3 3 T 2 2 T 2 2 T and λ is a parameter related to the process noise.
In this example, the triggering thresholds are set as τ 1 = τ 2 = 200 and the saturation levels are selected by ϱ 1 = ϱ 2 = 6000 . The other parameters are chosen as g = 9.8 m/s 2 , Υ = 4 × 10 4 kg/ms 2 , λ = 0.1   m 2 / s 3 , T = 1 s, κ 1 = 1.227 , κ 2 = 1.094 × 10 4 , π i = 1 ( i = 1 , 2 , , 11 ) , a k = b k = c k = d k = 1 .
The simulation results are displayed in Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10. The position of the target on x-axis x 1 , k and its estimate are plotted in Figure 6. Figure 7 shows the velocity of the target on x-axis x ˙ 1 , k and its estimate. The corresponding results for the target on y-axis are depicted in Figure 8 and Figure 9. The actual fault and its estimate are shown in Figure 10. From the simulation, it is observed that the proposed ETSFE algorithm is indeed effective in the ballistic object tracking system.

5. Conclusions

In this paper, we have addressed the state and fault estimation issues for nonlinear systems subject to sensor saturations. The ETM based scheduling protocol has been adopted to manage the measurement transmission, and thus relieving the communication burden. Relying on the received measurements, the state and fault estimator has been designed where the recursions of the upper bounds on the estimation error covariances for the state and fault have been given, respectively. Then, the derived upper bounds have been minimized by properly designing the estimator gain matrices. Finally, we have utilized two simulation examples to show the validity of the designed ETSFE algorithm. Our future research topics would be to extend the main results in this paper to the sensor networks where multiple coupling sensors are involved [22] and to apply the main results in practical engineering such as power systems [3].

Author Contributions

Conceptualization, C.H. and B.S.; methodology, C.H. and B.S.; software, C.H. and B.S.; validation, C.H., B.S., L.Z. and Y.S.; formal analysis, C.H. and B.S.; investigation, C.H. and B.S.; writing—original draft preparation, C.H. and B.S.; writing—review and editing, C.H., B.S., L.Z. and Y.S.; All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the Program of Shanghai Academic/Technology Research Leader under Grant 20XD1420100, the Natural Science Foundation of Shanghai under Grant 18ZR1401500, the China Postdoctoral Science Foundation under Grant 2018T110702, the Fundamental Research Funds for the Central Universities and Graduate Student Innovation Fund of Donghua University under Grant CUSF-DH-D-2019092.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Han, Y.; Wang, B.; Deng, Z.; Fu, M. Matching algorithm based on the nonlinear filter and similarity transformation for gravity-aided underwater navigation. IEEE/ASME Trans. Mechatron. 2018, 23, 646–654. [Google Scholar] [CrossRef]
  2. Kalman, R.E. A new approach to linear filtering and prediction problems. J. Basic Eng. Trans. 1960, 82, 35–45. [Google Scholar] [CrossRef] [Green Version]
  3. Qu, B.; Li, N.; Liu, Y.; Alsaadi, F.E. Estimation for power quality disturbances with multiplicative noises and correlated noises: A recursive estimation approach. Int. J. Syst. Sci. 2020, 51, 1200–1217. [Google Scholar] [CrossRef]
  4. Shen, B.; Wang, Z.; Wang, D.; Luo, J.; Pu, H.; Peng, Y. Finite-horizon filtering for a class of nonlinear time-delayed systems with an energy harvesting sensor. Automatica 2019, 100, 144–152. [Google Scholar] [CrossRef]
  5. Tan, H.; Shen, B.; Peng, K.; Liu, H. Robust recursive filtering for uncertain stochastic systems with amplify-and-forward relays. Int. J. Syst. Sci. 2020, 51, 1188–1199. [Google Scholar] [CrossRef]
  6. Gu, Z.; Shi, P.; Yue, D.; Ding, Z. Decentralized adaptive event-triggered H filtering for a class of networked nonlinear interconnected systems. IEEE Trans. Cybern. 2019, 49, 1570–1579. [Google Scholar] [CrossRef] [Green Version]
  7. Liu, J.; Yang, M.; Xie, X.; Peng, C.; Yan, H. Finite-time H filtering for state-dependent uncertain systems with event-triggered mechanism and multiple attacks. IEEE Trans. Circuits Syst. I Regul. Pap. 2020, 67, 1021–1031. [Google Scholar] [CrossRef]
  8. Mao, J.; Ding, D.; Song, Y.; Liu, Y.; Alsaadi, F.E. Event-based recursive filtering for time-delayed stochastic nonlinear systems with missing measurements. Signal Process. 2017, 134, 158–165. [Google Scholar] [CrossRef]
  9. Hu, J.; Wang, Z.; Liu, S.; Gao, H. A variance-constrained approach to recursive state estimation for time-varying complex networks with missing measurements. Automatica 2016, 64, 155–162. [Google Scholar] [CrossRef] [Green Version]
  10. Chan, J.C.L.; Tan, C.P.; Trinh, H.; Kamal, M.A.S. State and fault estimation for a class of non-infinitely observable descriptor systems using two sliding mode observers in cascade. J. Frankl. Inst. 2019, 356, 3010–3029. [Google Scholar] [CrossRef]
  11. Hosseini, I.; Fiacchini, M.; Karimaghaee, P.; Khayatian, A. Optimal reset unknown input observer design for fault and state estimation in a class of nonlinear uncertain systems. J. Frankl. Inst. 2020, 357, 2978–2996. [Google Scholar] [CrossRef]
  12. Zhang, K.; Jiang, B.; Shi, P.; Xu, J. Fault estimation observer design for discrete-time Takagi-Sugeno fuzzy systems based on piecewise lyapunov functions. IEEE Trans. Fuzzy Syst. 2012, 20, 192–200. [Google Scholar] [CrossRef]
  13. Song, Y.; Hu, J.; Chen, D.; Ji, D.; Liu, F. Recursive approach to networked fault estimation with packet dropouts and randomly occurring uncertainties. Neurocomputing 2016, 214, 240–349. [Google Scholar] [CrossRef]
  14. Zhang, K.; Jiang, B.; Shi, P.; Xu, J. Analysis and design of robust H fault estimation observer with finite-frequency specifications for discrete-time fuzzy systems. IEEE Trans. Cybern. 2015, 45, 1225–1235. [Google Scholar] [CrossRef]
  15. Jiang, S.; Fang, H. Fault estimation for nonlinear networked systems with time-varying delay and random packet dropout. Asian J. Control 2014, 16, 126–137. [Google Scholar] [CrossRef]
  16. Li, Z.; Wang, Z.; Ding, D.; Shu, H. H fault estimation with randomly occurring uncertainties, quantization effects and successive packet dropouts: The finite-horizon case. Int. J. Robust Nonlinear Control 2013, 25, 2671–2686. [Google Scholar] [CrossRef] [Green Version]
  17. Ichalal, D.; Marx, B.; Ragot, J. State estimation and fault detection for 2-D discrete-time systems. In Proceedings of the 26th Mediterranean Conference on Control and Automation, Zadar, Croatia, 19–22 June 2018; pp. 703–708. [Google Scholar]
  18. Zhao, D.; Wang, Y.; Li, Y.; Ding, S.X. H fault estimation for 2-D linear discrete time-varying systems based on Krein space method. IEEE Trans. Syst. Man Cybern. Syst. 2018, 48, 2070–2079. [Google Scholar] [CrossRef]
  19. Li, Q.; Shen, B.; Liu, Y.; Huang, T. Event-triggered H state estimation for discrete-time neural networks with mixed time delays and sensor saturations. Neural Comput. Appl. 2017, 28, 3815–3825. [Google Scholar] [CrossRef]
  20. Ma, L.; Wang, Z.; Lam, H.K. Event-triggered mean-square consensus control for time-varying stochastic multi-agent system with sensor saturations. IEEE Trans. Autom. Control 2017, 62, 3524–3531. [Google Scholar] [CrossRef] [Green Version]
  21. Hu, J.; Wang, Z.; Gao, H. Joint state and fault estimation for time-varying nonlinear systems with randomly occurring faults and sensor saturations. Automatica 2018, 97, 150–160. [Google Scholar] [CrossRef]
  22. Dong, H.; Wang, Z.; Lam, J.; Gao, H. Distributed filtering in sensor networks with randomly occurring saturations and successive packet dropouts. Int. J. Robust Nonlinear Control 2014, 24, 1743–1759. [Google Scholar] [CrossRef]
  23. Wang, Z.; Shen, B.; Liu, X. H filtering with randomly occurring sensor saturations and missing measurements. Automatica 2012, 48, 556–562. [Google Scholar] [CrossRef] [Green Version]
  24. Battistelli, G.; Chisci, L.; Selvi, D.A. Distributed Kalman filter with event-triggered communication and guaranteed stability. Automatica 2018, 93, 75–82. [Google Scholar] [CrossRef]
  25. Huang, C.; Shen, B.; Chen, H.; Shu, H. A dynamically event-triggered approach to recursive filtering with censored measurements and parameter uncertainties. J. Frankl. Inst. 2019, 356, 8870–8889. [Google Scholar] [CrossRef]
  26. Li, W.; Jia, Y.; Du, J. Event-triggered Kalman consensus filter over sensor networks. IET Control Theory Appl. 2016, 10, 103–110. [Google Scholar] [CrossRef]
  27. Liu, Y.; Wang, Z.; He, X.; Zhou, D. Event-triggered fault estimation for nonlinear systems with missing measurements. In Proceedings of the 33rd Chinese Control Conference, Nanjing, China, 28–30 July 2014; pp. 5533–5538. [Google Scholar]
  28. Li, Y.; Peng, L. Event-triggered fault estimation for stochastic systems over multi-hop relay networks with randomly occurring sensor nonlinearities and packet dropouts. Sensors 2018, 18, 731. [Google Scholar] [CrossRef] [Green Version]
  29. Li, Y.; Wu, Q.; Peng, L. Simultaneous event-triggered fault detection and estimation for stochastic systems subject to deception attacks. Sensors 2018, 18, 321. [Google Scholar] [CrossRef] [Green Version]
  30. Theodor, Y.; Shaked, U. Robust discrete-time minimum-variance filtering. IEEE Trans. Signal Process. 1996, 44, 181–189. [Google Scholar] [CrossRef]
  31. Shen, B.; Wang, Z.; Wang, D.; Liu, H. Distributed state-saturated recursive filtering over sensor networks under round-robin protocol. IEEE Trans. Cybern. 2020, 50, 3605–3615. [Google Scholar] [CrossRef]
  32. Farina, A.; Ristic, B.; Benvenuti, D. Tracking a ballistic target: Comparison of several nonlinear filters. IEEE Trans. Aerosp. Electron. Syst. 2002, 38, 854–867. [Google Scholar] [CrossRef]
Figure 1. Structure of the state and fault estimation.
Figure 1. Structure of the state and fault estimation.
Sensors 21 01242 g001
Figure 2. Trajectories x k 1 , x k 2 and their estimates.
Figure 2. Trajectories x k 1 , x k 2 and their estimates.
Sensors 21 01242 g002
Figure 3. Trajectories of fault f k and its estimate.
Figure 3. Trajectories of fault f k and its estimate.
Sensors 21 01242 g003
Figure 4. Log ( MSE k x ) of x k 1 , x k 2 and the traces of their upper bounds.
Figure 4. Log ( MSE k x ) of x k 1 , x k 2 and the traces of their upper bounds.
Sensors 21 01242 g004
Figure 5. Log ( MSE k f ) of f k and the trace of its upper bound.
Figure 5. Log ( MSE k f ) of f k and the trace of its upper bound.
Sensors 21 01242 g005
Figure 6. The position of the target x 1 , k and its estimate.
Figure 6. The position of the target x 1 , k and its estimate.
Sensors 21 01242 g006
Figure 7. The velocity of the target x ˙ 1 , k and its estimate.
Figure 7. The velocity of the target x ˙ 1 , k and its estimate.
Sensors 21 01242 g007
Figure 8. The position of the target x 2 , k and its estimate.
Figure 8. The position of the target x 2 , k and its estimate.
Sensors 21 01242 g008
Figure 9. The velocity of the target x ˙ 2 , k and its estimate.
Figure 9. The velocity of the target x ˙ 2 , k and its estimate.
Sensors 21 01242 g009
Figure 10. The actual fault and its estimate.
Figure 10. The actual fault and its estimate.
Sensors 21 01242 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, C.; Shen, B.; Zou, L.; Shen, Y. Event-Triggering State and Fault Estimation for a Class of Nonlinear Systems Subject to Sensor Saturations. Sensors 2021, 21, 1242. https://doi.org/10.3390/s21041242

AMA Style

Huang C, Shen B, Zou L, Shen Y. Event-Triggering State and Fault Estimation for a Class of Nonlinear Systems Subject to Sensor Saturations. Sensors. 2021; 21(4):1242. https://doi.org/10.3390/s21041242

Chicago/Turabian Style

Huang, Cong, Bo Shen, Lei Zou, and Yuxuan Shen. 2021. "Event-Triggering State and Fault Estimation for a Class of Nonlinear Systems Subject to Sensor Saturations" Sensors 21, no. 4: 1242. https://doi.org/10.3390/s21041242

APA Style

Huang, C., Shen, B., Zou, L., & Shen, Y. (2021). Event-Triggering State and Fault Estimation for a Class of Nonlinear Systems Subject to Sensor Saturations. Sensors, 21(4), 1242. https://doi.org/10.3390/s21041242

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop