Next Article in Journal
Prediction and Evaluation of Zero Order Entropy Changes in Grammar-Based Codes
Previous Article in Journal
Muscle Fatigue Analysis of the Deltoid during Three Head-Related Static Isometric Contraction Tasks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Minimum Entropy Active Fault Tolerant Control of the Non-Gaussian Stochastic Distribution System Subjected to Mean Constraint

School of Electrical Engineering, Zhengzhou University, Zhengzhou 450001, China
*
Author to whom correspondence should be addressed.
Entropy 2017, 19(5), 218; https://doi.org/10.3390/e19050218
Submission received: 23 March 2017 / Revised: 24 April 2017 / Accepted: 3 May 2017 / Published: 11 May 2017

Abstract

:
Stochastic distribution control (SDC) systems are a group of systems where the outputs considered is the measured probability density function (PDF) of the system output whilst subjected to a normal crisp input. The purpose of the active fault tolerant control of such systems is to use the fault estimation information and other measured information to make the output PDF still track the given distribution when the objective PDF is known. However, if the target PDF is unavailable, the PDF tracking operation will be impossible. Minimum entropy control of the system output can be considered as an alternative strategy. The mean represents the center location of the stochastic variable, and it is reasonable that the minimum entropy fault tolerant controller can be designed subjected to mean constraint. In this paper, using the rational square-root B-spline model for the shape control of the system output probability density function (PDF), a nonlinear adaptive observer based fault diagnosis algorithm is proposed to diagnose the fault. Through the controller reconfiguration, the system entropy subjected to mean restriction can still be minimized when fault occurs. An illustrative example is utilized to demonstrate the use of the minimum entropy fault tolerant control algorithms.

1. Introduction

Fault diagnosis (FD) and fault-tolerant control (FTC) has long been regarded as an important and integrated part in control systems. The research of this aspect has been carried out for nearly thirty years, where the major research is aimed at deterministic systems. However, most practical systems are subjected to random inputs from sensor noises [1], random disturbances [2] or random parameter changes. In order to improve the reliability of the control system, the research of fault diagnosis and fault-tolerant control for stochastic dynamic systems has been one of the liveliest research areas in control theory and practice [3]. For stochastic systems, the so-far obtained FDD (fault detection and diagnosis) approaches can be classified as
  • The system identification technique [4];
  • Observer or filter-based method [5];
  • The statistic approach based on the Bayesian theorem, Monte Carlo approach, likelihood method, and hypothesis test technique [6].
Moreover, in a safety-critical system, the role of a fault-tolerant control system is extremely important. One of its functions is to steer the process to a safe state whenever undesirable events occur. A fault-tolerant controller can be able to attenuate the effects caused by random variables.
For most control algorithms, it is supposed that system noise obeys Gaussian distribution. Traditional statistics such as mean and variance can be used to represent the randomness of the system output. However, for non-Gaussian systems, the variance cannot fully represent the randomness of the system output. A more general measure, rather than just the variance of the system output, needs to be considered for the randomness of the system output of non-Gaussian stochastic systems. Therefore, entropy concept [7] is proposed to measure the uncertainty of the system output of stochastic systems.
In existing FTC results of SDC systems, it is assumed that the target PDF is pre-specified. Nonlinear filters are designed to detect and diagnose the fault of the stochastic distribution system based on square-root approximation [8]. However, fault tolerant control is not considered. A nonlinear adaptive observer-based fault diagnosis algorithm is proposed to diagnose the fault in the dynamic part of the non-Gaussian SDC systems based on the rational square-root B-spline approximation model [9]. Through the controller reconfiguration, a good output PDF tracking can still be realized when fault occurs. An iterative learning observer based fault estimation algorithm is proposed [10] and the fault tolerant controller is designed so as to make the post-fault PDF still track the given distribution. However, in most cases, the objective PDF can’t be determined in advance and, as a result, the reasonable fault tolerant tracking controller cannot be designed. Therefore, minimum entropy control of the system output can be considered as an alternative strategy. A minimum entropy control method of closed-loop tracking errors is proposed [11]. A nonlinear observer is designed to obtain the pseudo-weights vector [12]. Then, the PDF tracking control and the minimum entropy control are developed. An integrated design of fault diagnosis and fault tolerant control (FTC) for nonlinear systems using Takagi–Sugeno (T–S) fuzzy models is proposed in [13]; meanwhile in the presence of the model uncertainty along with actuator/sensor faults and external disturbance. An augmented state unknown input observer is proposed to estimate the fault and the system state simultaneously, and using the fault estimation information, an FTC controller is developed to ensure the robust stability of the closed-loop system. The design scheme of thruster failure tolerant control for underwater vehicles is proposed [14], where each thruster of the underwater vehicle can rotate, offering a significant advantage to optimize its control. Active and passive fault-tolerant control systems (FTCSs) aer demystified by examining the similarities and differences between these two approaches from both philosophical and practical points of view in [15]. A new method is proposed [16] to minimize the closed loop randomness for general dynamic stochastic systems using the entropy concept. However, a fault item is not involved in this literature. This forms the point of this paper, where the entropy concept is used for design of the required fault tolerant controllers of stochastic systems.
In this paper, the rational square-root B-spline model [17] is used to represent the dynamics between the output PDF and input. The considered system in [17] is subjected to any arbitrary bounded random input and the purpose of the control input design is to make the output probability density function of the system output track a given distribution function as close as possible, that different from the objective PDF known in literature [17], the target PDF in this paper can not be determined in advance. Considering the existence of the fault, a nonlinear adaptive observer based fault diagnosis algorithm is proposed to diagnose the fault. Through the controller reconfiguration, the system entropy can still be minimized when fault occurs. The contribution of this paper is that the entropy concept is introduced to design of fault tolerant control for non-Gaussian stochastic distribution systems when the objective PDF can not be determined in advance.
The rest of this paper is organized as follows. Section 2 presents the model description, where a nonlinear dynamic relationship is considered. An observer is constructed in Section 3 for fault detection. In Section 4, the nonlinear adaptive observer based fault diagnosis algorithm is proposed, and the rational square-root B-spline model, which is different from that in [10], is used for the shape control of the system output PDF. Section 5 gives the design process of fault tolerant control. Simulation results of FD and FTC are presented in Section 6, followed by some concluding remarks in Section 7.

2. Model Description

Denote η [ a , b ] as a uniformly bounded stochastic process variable and assume that it represents the output of a stochastic system at sample time t. Denote u ( t ) as the control input that controls the distribution of η ( t ) . At any time, η ( t ) can be characterized by its probability density function γ ( y , u ( t ) ) , which is defined by
P ( a η ( t ) < ε | u ( t ) ) = a ε γ ( y , u ( t ) ) d y ,
where P ( a η ( t ) < ε | u ( t ) ) represents the probability of the output y ( t ) lying inside the interval [ a , ε ) when u ( t ) is applied to the system. It also means that the shape of the output probability density function, γ ( y , u ( t ) ) , is controlled by a crisp input u ( t ) . It is assumed that interval [ a , b ] is known and the output probability density functions γ ( y , u ( t ) ) is continuous and bounded. Using the rational square-root B-spline function approximation principle [17], the rational square-root B-spline expansion as follows,
γ ( y , u ( t ) ) = i = 1 m ω i ( u ( t ) ) ϕ i ( y ) i , j = 1 m ω i ( u ( t ) ) ω j ( u ( t ) ) a b ϕ i ( y ) ϕ j ( y ) d y ,
where ϕ i ( y ) ( i = 1 , · · · , m ) are the pre-specified basis functions defined on [ a , b ] and ω i ( i = 1 , · · · , m ) are the corresponding weights.
For most systems, a dynamic relationship between ω i and u is expressed. In this paper, a nonlinear dynamic relationship will be considered. This leads to the following model for the considered dynamic stochastic system:
x ˙ ( t ) = A x ( t ) + H g ( x ( t ) ) + B u ( t ) + G F ( t ) , V ( t ) = D x ( t ) ,
γ ( y , u ( t ) ) = C ( y ) V ( t ) V ( t ) T E V ( t ) ,
where x ( t ) R n is the state vector, V ( t ) R m is the output weight vector, u ( t ) R r is the control input vector and F ( t ) R n is the fault vector. A R n × n , B R n × r , D R m × n , E R m × m , H R n × n and G R n × n are system parameter matrices. Equation (2) represents the dynamic model of the weights vector. Equation (3) describes the static output PDF model using the rational square-root B-spline expression, where
V ( t ) = [ ω 1 ( u ( t ) ) , ω 2 ( u ( t ) ) , . . . , ω m ( u ( t ) ) ] T R m , V ( t ) 0 ,
C ( y ) = [ ϕ 1 ( y ) , ϕ 2 ( y ) , . . . , ϕ m ( y ) ] R 1 × m ,
E = a b C ( y ) T C ( y ) d y .
Assumption A1.
The nonlinear function g ( x ( t ) ) should satisfy the following Lipschitz condition:
g ( x i ) - g ( x j ) t m x x i - x j ,
where m x is a Lipschitz constant.

3. Fault Detection

The purpose of the fault detection is to use input vector in Equation (2) and output probability density function γ ( y , u ( t ) ) to detect the fault F. In order to detect such a fault, an observer is constructed [18] as follows:
x ^ ˙ ( t ) = A x ^ ( t ) + H g ( x ^ ( t ) ) + B u ( t ) + K ε d ( t ) , ε d ( t ) = a b σ ( y ) [ γ ^ ( y , u ( t ) ) - γ ( y , u ( t ) ) ] d y ,
γ ^ ( y , u ( t ) ) = C ( y ) V ^ ( t ) V ^ ( t ) T E V ^ ( t ) ,
where x ^ ( t ) R n is the estimated state, ε d ( t ) is the residual signal, σ ( y ) is the pre-specified weighting vector defined on [ a , b ] and K is the gain of the observer.
Denote the observation error vector as
e d ( t ) = x ^ ( t ) - x ( t ) .
The residual signal can be calculated as follows:
ε d ( t ) = a b σ ( y ) [ γ ^ ( y , u ( t ) ) - γ ( y , u ( t ) ) ] d y = a b σ ( y ) C ( y ) d y [ V ^ ( t ) V ^ ( t ) T E V ^ ( t ) - V ( t ) V ( t ) T E V ( t ) ] d y = Σ D e d ( t ) V ^ ( t ) T E V ^ ( t ) + Σ V ( t ) V ^ ( t ) T E V ^ ( t ) - Σ V ( t ) V ( t ) T E V ( t ) ,
where Σ = a b σ ( y ) C ( y ) d y .
Using Equations (2) and (4), the error dynamic system can be formulated that
e ˙ d ( t ) = x ^ ˙ ( t ) - x ˙ ( t ) = A e + H [ g ( x ^ ( t ) ) - g ( x ( t ) ) ] - G F ( t ) + K ε d ( t ) = A e + H [ g ( x ^ ( t ) ) - g ( x ( t ) ) ] - G F ( t ) + K Σ D e d ( t ) V ^ ( t ) T E V ^ ( t ) + K ( Σ V ( t ) V ^ ( t ) T E V ^ ( t ) - Σ V ( t ) V ( t ) T E V ( t ) ) .
Assume that the basis function ϕ i ( y ) and σ ( y ) are selected such that the pair ( A , Σ ) is observable. The gain L 1 can be selected such that matrix M = A + L 1 Σ D is H u r w i t z .
Lemma 1.
There exists a λ 1 ( T 1 = λ m i n ( E ) / λ m a x ( E ) , T 2 = λ m a x ( E ) / λ m i n ( E ) ) , T 1 | λ 1 | T 2 ) such that the following equation
V T E V - V ^ T E V ^ = λ 1 ( V T V - V ^ T V ^ ) = λ 1 ( V - V ^ )
holds.
Theorem 1.
It is supposed that there exist two positive definite matrices P and Q such that the following equalities
M T P + P M = - Q , P = P T > 0 , Q = Q T > 0
hold, and then the observation error system is asymptotically stable.
Proof. 
When no fault occurs, F = 0 . Denote K = L 1 V ^ ( t ) T E V ^ ( t ) , and L 1 is the gain matrix of the observer. Then, the error dynamic system can be further formulated as
e ˙ d ( t ) = ( A + L 1 Σ D ) e d ( t ) + H [ g ( x ^ ( t ) ) - g ( x ( t ) ) ] + L 1 Σ V ( t ) ( V ( t ) - V ^ ( t ) V ( t ) T E V ( t ) ) .
The Lyapunov function is selected as follows:
π = 1 2 e d T ( t ) P e d ( t ) .
It can be shown that
π ˙ = 1 2 e ˙ d T ( t ) P e d ( t ) + 1 2 e d T ( t ) P e ˙ d ( t ) = - 1 2 e d T ( t ) Q e d ( t ) + e d T ( t ) P H [ g ( x ^ ( t ) ) - g ( x ( t ) ) ] + e d T ( t ) P L 1 Σ V ( t ) ( λ 1 ( V ( t ) - V ^ ( t ) ) V ( t ) T E V ( t ) ) - 1 2 ( λ m i n ( Q ) - 2 T 2 P L 1 Σ D E - 2 m x P H ) e d ( t ) 2 ,
when
λ m i n ( Q ) 2 T 2 P L 1 Σ D E + 2 m x P H
hold, and then π ˙ < 0 . The proof is completed. ☐
The fault can be detected when ε d ( t ) > τ , where τ ( τ > 0 ) is a pre-specified threshold.

4. Fault Diagnosis Algorithm

Once the fault has been detected, the fault diagnosis needs to be carried out in order to estimate the size of the fault. A fault diagnosis observer is constructed [19] as follows:
x ^ ˙ m ( t ) = A x ^ m ( t ) + H g ( x ^ m ( t ) ) + B u ( t ) + K d ε m ( t ) + G F ^ ( t ) ε m ( t ) = a b σ ( y ) [ γ m ( y , u ( t ) ) - γ ( y , u ( t ) ) ] d y γ m ( y , u ( t ) ) = C ( y ) V m ( t ) V m ( t ) T E V m ( t ) ,
where F ^ is the estimation of F. It is assumed that e m = x ^ m - x , F ˜ = F ^ - F , and then
ε m ( t ) = a b σ ( y ) [ γ m ( y , u ( t ) ) - γ ( y , u ( t ) ) ] d y = Σ D e m V m ( t ) T E V m ( t ) + Σ V ( t ) V m ( t ) T E V m ( t ) - Σ V ( t ) V ( t ) T E V ( t ) .
Denote K d = L 2 V m ( t ) T E V m ( t ) , and the gain L 2 can be selected such that matrix ( A + L 2 Σ D ) is H u r w i t z . Then, the error dynamic system can be obtained as follows:
e ˙ m = A e m + H [ g ( x ^ m ( t ) ) - g ( x ( t ) ) ] + K d ε m ( t ) + G F ˜ ( t ) = ( A + L 2 Σ D ) e m + H [ g ( x ^ m ( t ) ) - g ( x ( t ) ) ] + G F ˜ + L 2 Σ V ( t ) ( λ 2 ( V ( t ) - V m ( t ) ) V ( t ) T E V ( t ) ) .
Assumption A2.
Suppose that F ( t ) N 2 , where N is a positive constant.
Theorem 2.
It is supposed that there exist constant Γ > 0 and two positive definite matrices P 1 , Q 1 such that the following equality holds:
( A + L 2 Σ D ) T P 1 + P 1 ( A + L 2 Σ D ) = - Q 1 ,
and the adaptive tuning law is selected as follows:
d F ^ d t = - Γ V m ( t ) T E V m ( t ) ε m ( F ^ N 2 ) , 0 ( F ^ N 2 ) ,
and then the error dynamic system (14) is stable.
Proof. 
The Lyapunov function is selected as follows:
π 1 = 1 2 e m T P e m + 1 2 F ˜ T F ˜ ,
and then the first derivative of π 1 is formulated as follows:
π ˙ 1 = 1 2 e ˙ m T P 1 e m + 1 2 e m T P 1 e ˙ m + F ˜ T F ˜ ˙ = 1 2 e m T [ ( A + L 2 Σ D ) T P 1 + P 1 ( A + L 2 Σ D ) ] e m + e m T P 1 H [ g ( x ^ m ( t ) ) - g ( x ( t ) ) ] + e m T P 1 G F ˜ ( t ) + F ˜ T F ˜ ˙ + e m T P 1 L 2 Σ V ( t ) λ 2 ( V ( t ) - V m ( t ) ) V ( t ) T E V ( t ) - 1 2 ( λ m i n ( Q 1 ) - 2 P 1 L 2 Σ D T 2 E - 2 m x P 1 H ) e m 2 + e m ( P 1 G - D T Σ T Γ + Γ Σ D T 2 E ) F ˜ = - δ 1 e m 2 + δ 2 e m F ˜ ,
where
δ 1 = 1 2 ( λ m i n ( Q 1 ) - 2 P 1 L 2 Σ D T 2 E - 2 m x P 1 H ) ,
δ 2 = ( P 1 G - D T Σ T Γ + Γ Σ D T 2 E ) .
From Assumption 2 and Theorem 2, it can be obtained that
π ˙ 1 - δ 1 ( e m - N δ 2 2 δ 1 ) 2 + N 2 ( δ 2 ) 2 4 ( δ 1 ) 2 .
When e m N δ 2 ( δ 1 + δ 1 ) 2 δ 1 , it can be obtained that π ˙ 1 0 . This means that when fault occurs in system (2), the observation error system is stable. ☐

5. Fault Tolerant Control

Once the fault has been diagnosed, the controller that is designed for the healthy system should be modified so as to compensate for the performance losses caused by the fault in the system.
When the objective PDF can not be determined in advance, or no objective PDF can be used, control target of the SDC system can be translated into controlling the entropy of the system output. It has been pointed out that the entropy has been used to characterize the uncertainty of the output variables for non-Gaussian stochastic systems. Thus, the uncertainty is expected to be as small as possible.
When the position of the output variable’s PDF is not determined, the minimum entropy of the output variables is generally difficult to be determined in practice. Therefore, control target subjected to mean constraint will be suitable for the solution of the minimum entropy controller.
Since the mean and entropy of a non-Gaussian stochastic variable are independent of each other, these two terms may be addressed separately. For this purpose, the following performance function is used
J ( u ( t ) ) = ( μ - μ g ) 2 - a b γ ( y , u ( t ) ) l n γ ( y , u ( t ) ) d y + 1 2 u T ( t ) R u ( t ) .
It can be seen that the first term in Equation (20) is the error between the mean and its target value μ g , where μ = a b y γ ( y , u ) d y . The second term is the Shannon entropy of the output variables, and the third term is a natural quadratic constraint for the control input, where R = R T > 0 .
The purpose of the controller design is to find the required optimal control input such that the performance index function is minimized. This purpose can be realized by selecting a control input such that the performance function does not increase. The performance index of the close-loop system can be seen as a Lyapunov function. In order to minimize J in Equation (20), the first order derivative of J can be readily formulated that
d J d t = ( μ - μ g ) 2 V ( t ) V ˙ ( t ) + u T ( t ) R u ˙ ( t ) - a b γ ( y , u ( t ) ) V ( t ) V ˙ ( t ) ( l n γ ( y , u ( t ) ) + 1 ) d y .
At this stage, by selecting the following control input
u T ( t ) R u ˙ ( t ) t = - λ | μ - μ g | - ( μ - μ g ) 2 V ( t ) V ˙ ( t ) + a b γ ( y , u ( t ) ) V ( t ) V ˙ ( t ) ( l n γ ( y , u ( t ) ) + 1 ) d y = ( a b γ ( y , u ( t ) ) V ( t ) ( l n γ ( y , u ( t ) ) + 1 ) d y - ( μ - μ g ) 2 V ( t ) ) V ˙ ( t ) - λ | μ - μ g | ,
where λ > 0 , it can then be formulated that
d J d t = - λ | μ - μ g | < 0 .
Therefore, the performance index is convergent.
From Equation (2), it can be obtained that
V ˙ ( t ) = D x ˙ ( t ) = D A x ( t ) + D H g ( x ( t ) ) + D B u ( t ) + D G F ( t ) .
From Equation (3), it can be obtained that
γ ( y , u ( t ) ) = ( V ( t ) T ( C ( y ) ) T ) ( C ( y ) V ( t ) ) V ( t ) T E V ( t ) μ = a b y γ ( y , u ( t ) ) d y = a b y [ V ( t ) T ( C ( y ) ) T ] [ C ( y ) V ( t ) ] V ( t ) T E V ( t ) d y = V ( t ) T ( a b y ( C ( y ) ) T C ( y ) d y ) V ( t ) V ( t ) T E V ( t ) ,
and then
γ ( y , u ( t ) ) V = 2 ( V ( t ) T C ( y ) T ) ( C ( y ) ( V ( t ) T E V ( t ) ) - ( C ( y ) V ( t ) ) ( V ( t ) T E ) ) ( V ( t ) T E V ( t ) ) 2 , μ ( t ) V ( t ) = 2 V ( t ) T a b y C ( y ) T C ( y ) d y ( V ( t ) T E V ( t ) - V ( t ) ( V ( t ) T E ) ) ( V ( t ) T E V ( t ) ) 2 .
When no fault occurs, F = 0 , using Equation (24), it can be obtained that
u T ( t ) R u ( t ) = - λ | a b y γ ( y , u ( t ) ) d y - μ g | + a b [ 2 ( V ( t ) T E V ( t ) ) [ V ( t ) T ( C ( y ) ) T ] C ( y ) ( V ( t ) T E V ( t ) ) 2 - 2 [ V ( t ) T ( C ( y ) ) T ] [ C ( y ) V ( t ) ] ( V ( t ) T E ) ( V ( t ) T E V ( t ) ) 2 ] ( ln γ ( y , u ( t ) ) + 1 ) d y × ( D A x ( t ) + D H g ( x ( t ) ) + D B u ( t ) ) - 2 ( μ - μ g ) × [ 2 ( V ( t ) T E V ( t ) ) V ( t ) T a b y ( C ( y ) ) T C ( y ) d y ( V ( t ) T E V ( t ) ) 2 - 2 V ( t ) T a b y ( C ( y ) ) T C ( y ) d y V ( t ) ( V ( t ) T E ) ( V ( t ) T E V ( t ) ) 2 ] × ( D A x ( t ) + D H g ( x ( t ) ) + D B u ( t ) ) .
When the fault occurs and has been diagnosed, it can then be rewritten as
u 1 T ( t ) R u 1 ( t ) = - λ | a b y γ ( y , u 1 ( t ) ) d y - μ g | + a b [ 2 ( V m ( t ) T E V m ( t ) ) [ V m ( t ) T ( C ( y ) ) T ] C ( y ) ( V m ( t ) T E V m ( t ) ) 2 - 2 [ V m ( t ) T ( C ( y ) ) T ] [ C ( y ) V m ( t ) ] ( V m ( t ) T E ) ( V m ( t ) T E V m ( t ) ) 2 ] ( ln γ ( y , u 1 ( t ) ) + 1 ) d y × ( D A x m ( t ) + D H g ( x m ( t ) ) + D B u 1 ( t ) + D G F ^ ( t ) ) - 2 ( μ - μ g ) [ 2 ( V m ( t ) T E V m ( t ) ) V m ( t ) T a b y ( C ( y ) ) T C ( y ) d y ( V m ( t ) T E V m ( t ) ) 2 - 2 V m ( t ) T a b y ( C ( y ) ) T C ( y ) d y V m ( t ) ( V m ( t ) T E ) ( V m ( t ) T E V m ( t ) ) 2 ] × ( D A x m ( t ) + D H g ( x m ( t ) ) + D B u 1 ( t ) + D G F ^ ( t ) ) ,
where u 1 is the reconfigured controller. The state x and fault F can be substituted by the state of the fault diagnosis observer state x m and the fault estimation F ^ , respectively.

6. Simulation Example

The following B-spline is selected as ϕ i ( y ) ( i = 1 , 2 , 3 )
ϕ 1 ( y ) = ( 1 6 y 3 + 3 2 y 2 + 9 2 y + 9 2 ) I 1 + ( - 1 2 y 3 - 5 2 y 2 - 7 2 y - 5 6 ) I 2 + ( 1 2 y 3 + 1 2 y 2 - 1 2 y + 1 6 ) I 3 + ( - 1 6 y 3 + 1 2 y 2 - 1 2 y + 1 6 ) I 4 , ϕ 2 ( y ) = ( 1 6 y 3 + y 2 + 2 y + 4 3 ) I 2 + ( - 1 2 y 3 - y 2 + 2 3 ) I 3 + ( 1 2 y 3 - y 2 + 2 3 ) I 4 + ( - 1 6 y 3 + y 2 - 2 y + 4 3 ) I 5 , ϕ 3 ( y ) = ( 1 6 y 3 + 1 2 y 2 + 1 2 y + 1 6 ) I 3 + ( - 1 2 y 3 + 1 2 y 2 + 1 2 y + 1 6 ) I 4 + ( 1 2 y 3 - 5 2 y 2 + 7 2 y - 5 6 ) I 5 + ( - 1 6 y 3 + 3 2 y 2 - 9 2 y + 9 2 ) I 6 ,
where
I i ( y ) = 1 y [ i - 4 , i - 3 ] 0 o t h e r w i s e ( i = 1 , 2 , 3 , 4 , 5 , 6 ) ,
Σ = a b σ ( y ) C ( y ) d y = 1 1 1 .
It is assumed that the system dynamics can be represented as
x ˙ ( t ) = A x ( t ) + H g ( x ( t ) ) + B u ( t ) + G F ( t ) V ( t ) = D x ( t ) γ ( y , u ( t ) ) = C ( y ) V ( t ) V ( t ) T E V ( t ) ,
where g ( x ( t ) ) = s i n ( 0 . 8 x ( t ) ) ,
A = - 3 . 1866 - 1 . 0208 - 0 . 2857 - 3 . 5050 , B = 2 . 6 1 1 - 2 . 3 ,
G = 1 1 , D = 0 . 1 0 . 9 0 . 2 0 . 1 1 0 . 6 , H = 1 . 9 0 0 - 0 . 505 .
Positive definite matrices P 1 and Q 1 are selected as follows:
P 1 = 0 . 3505 0 . 0958 0 . 0958 0 . 5123 , Q 1 = 1 0 0 1 .
It can be solved via linear matrix inequality (LMI) related to Equations (9) and (15) to obtain the detect observer gain L 1 = 0 . 3 - 3 . 15 diagnostic observer gain L 2 = 0 . 3036 1 . 1448 and the parameters of fault estimation adaptive tuning rule Γ = 7 . 95 .
To simulate the algorithm, it is assumed that the three types of faults are constructed as follows:
The step fault:
F 1 = 0 t < 5 s 1 t 5 s ,
The intermittent fault:
F 2 = 1 , 2 < t < 12 ; 25 < t < 35 ; 45 < t < 50 0 , e l s e ,
The incipient fault:
F 3 = 1 + 0 . 3 sin ( - 0 . 5 ( t - 40 ) ) ,
The ramp-like fault:
F 4 = 1 t < 10 1 + 0 . 4 ( t - 10 ) 10 t < 25 7 - 0 . 2 ( t - 25 ) 25 t 50 .
Case 1: When a step fault occurs, the fault detection signal is shown in Figure 1 and the fault diagnosis result is shown in Figure 2. From the two figures, it can be concluded that desired fault detection and diagnosis results have been obtained. If the system works normally, then we have F 1 = 0 . Figure 3 shows the result of the output entropy and Figure 4 shows the result of the output PDF. Figure 5 presents the response of the mean. In the simulation, fault F 1 is set to 1 after 5 s. Figure 6 shows the changes of the entropy when a fault has occurred in the system. It can be seen that the control input is driving the system towards the direction of less randomness. From Figure 7, it can be seen that the post-fault PDF can still follow the faultless PDF, leading to good fault tolerant control results. Figure 8 presents the response of the mean when a fault has occurred in the system, and it can be seen that the mean value levels off and approaches the objective value.
Case 2: When the intermittent fault, incipient fault and ramp-like fault occur, the fault diagnosis results are shown in Figure 9, Figure 10 and Figure 11. From the three figures, it can be concluded that diagnosis results have been obtained. Figure 12, Figure 13 and Figure 14 show the change of the entropy when the three faults have occurred in the system. It can be seen that the control input is driving the system towards the direction of less randomness.

7. Concluding Remarks

In this paper, when objective PDFs can not be determined in advance, the entropy concept is applied to the design of a fault tolerant controller for a non-Gaussian stochastic distribution system. The mean represents the center location of the stochastic variable, and it is reasonable that a minimum entropy fault tolerant controller can be designed under the restriction of the mean. A rational square-root B-spline function approximation model and the nonlinear dynamic system have been used to formulate the stochastic distribution system. A performance function that includes the entropy term and a mean constraint term is used to design the fault tolerant controller. The performance index is regarded as a Lyapunov function. Using the Lyapunov stability criterion, the stability of the whole system can be guaranteed. Furthermore, through the controller reconfiguration, the system entropy can still be minimized. The simulations have further confirmed the proposed fault diagnosis and fault tolerant control results.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (NSFC) under Grant 61374128, the Science and Technology Innovation Talents 14HASTIT040 in Colleges and Universities in Henan Province, China, and the Open Project of the Chinese State Key Laboratory of Synthetical Automation for Process Industries.

Author Contributions

Haokun Jin, Yacun Guan and Lina Yao conceived and designed the experiments; Haokun Jin and Yacun Guan performed the experiments; Haokun Jin analyzed the data; Lina Yao contributed analysis tools; Haokun Jin and Yacun Guan wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Juricic, D.; Zele, M. Robust detection of sensor faults by means of a statistical test. Automatica 2002, 38, 737–742. [Google Scholar] [CrossRef]
  2. Guo, L.; Wang, H. Bibliographical review on reconfigurable fault-tolerant control systems. Annu. Rev. Control 2006, 51, 695–700. [Google Scholar]
  3. Wang, H. Bounded Dynamic Stochastic Systems: Modelling and Control; Springer: London, UK, 2000. [Google Scholar]
  4. Isermann, R. Model-based fault-detection and diagnosis-status and applications. Annu. Rev. Control 2005, 29, 71–85. [Google Scholar] [CrossRef]
  5. Frank, P.M.; Ding, X. Survey of robust residual generation and evaluation methods in observer-based fault detection systems. J. Process Control 1997, 7, 403–424. [Google Scholar] [CrossRef]
  6. Li, P.; Kadirkamanathan, V. Fault detection and isolation in nonlinear stochastic systems a combined adaptive Monte Carlo filtering and likelihood ratio approach. Int. J. Control 2004, 77, 1101–1104. [Google Scholar] [CrossRef]
  7. Chen, R.H.; Mingori, D.L.; Speye, J.L. Optimal stochastic fault detection filter. Automatica 2003, 39, 377–390. [Google Scholar] [CrossRef]
  8. Guo, L.; Wang, H. Fault detection and diagnosis for general stochastic systems using B-spline expansions and nonlinear filter. IEEE Trans. Circuits Syst. I 2005, 52, 1644–1652. [Google Scholar]
  9. Yao, L.N.; Wang, A.P.; Wang, H. Fault detection, diagnosis and tolerant control for non-Gaussian stochastic distribution systems using a rational square-root approximation model. Int. J. Model. Identif. Control 2008, 3, 162–172. [Google Scholar] [CrossRef]
  10. Yao, L.N.; Qin, J.F.; Wang, H.; Jiang, B. Design of new fault diagnosis and fault tolerant control scheme for non-Gaussian singular stochastic distribution systems. Automatica 2012, 48, 2305–2313. [Google Scholar] [CrossRef]
  11. Yue, H.; Wang, H. Minimum entropy control of closed-loop tracking errors for dynamic stochastic systems. IEEE Trans. Autom. Control 2003, 48, 118–122. [Google Scholar]
  12. Zhou, J.L.; Yue, H.; Wang, H. Shaping of output PDF based on the rational square-root B-spline model. Acta Autom. Sin. 2005, 31, 343–351. [Google Scholar]
  13. Lan, J.L.; Patton, R.J. Integrated design of fault-tolerant control for nonlinear systems based on fault estimation and T-S fuzzy modelling. IEEE Trans. Fuzzy Syst. 2016, 99. [Google Scholar] [CrossRef]
  14. Ciabattoni, L.; Fasano, A.; Ferracuti, F.; Freddi, A.; Longhi, S.; Monteriu, A. A thruster failure tolerant control scheme for underwater vehicles. In Proceedings of the 12th IEEE/ASME International Conference on Mechatronic and Embedded Systems and Applications (MESA), Auckland, New Zealand, 29–31 August 2016; pp. 1–6. [Google Scholar]
  15. Jiang, J.; Yu, X. Fault-tolerant control systems: A comparative study between active and passive approaches. Annu. Rev. Control 2012, 36, 60–72. [Google Scholar] [CrossRef]
  16. Wang, H. Minimum entropy control of non-Gaussian dynamic stochastic systems. IEEE Trans. Autom. Control 2002, 47, 398–403. [Google Scholar] [CrossRef]
  17. Wang, H.; Zhang, J.H. Bounded Stochastic distributions control for pseudo armax stochastic systems. IEEE Trans. Autom. Control 2001, 46, 486–490. [Google Scholar] [CrossRef]
  18. Yao, L.N.; Peng, B. Fault diagnosis and fault tolerant control for the non-Gaussian time-delayed stochastic distribution control system. In Proceedings of the American Control Conference (ACC), Washington, DC, USA, 17–19 June 2013; pp. 5409–5414. [Google Scholar]
  19. Yao, L.N.; Feng, L. Fault diagnosis and fault tolerant tracking control for the non-Gaussian singular time-delayed stochastic distribution system with PDF approximation error. Neurocomputing 2016, 175, 538–543. [Google Scholar] [CrossRef]
Figure 1. Fault detection signal.
Figure 1. Fault detection signal.
Entropy 19 00218 g001
Figure 2. The fault and fault estimation.
Figure 2. The fault and fault estimation.
Entropy 19 00218 g002
Figure 3. The response of entropy when system is normal.
Figure 3. The response of entropy when system is normal.
Entropy 19 00218 g003
Figure 4. The response of the output PDF when F = 0.
Figure 4. The response of the output PDF when F = 0.
Entropy 19 00218 g004
Figure 5. The response of the mean when system is normal.
Figure 5. The response of the mean when system is normal.
Entropy 19 00218 g005
Figure 6. The response of entropy after fault-tolerant control.
Figure 6. The response of entropy after fault-tolerant control.
Entropy 19 00218 g006
Figure 7. The output PDF with fault-tolerant control.
Figure 7. The output PDF with fault-tolerant control.
Entropy 19 00218 g007
Figure 8. The response of the mean after fault-tolerant control.
Figure 8. The response of the mean after fault-tolerant control.
Entropy 19 00218 g008
Figure 9. The intermittent fault and fault estimation.
Figure 9. The intermittent fault and fault estimation.
Entropy 19 00218 g009
Figure 10. The incipient fault and fault estimation.
Figure 10. The incipient fault and fault estimation.
Entropy 19 00218 g010
Figure 11. The ramp-like fault and fault estimation.
Figure 11. The ramp-like fault and fault estimation.
Entropy 19 00218 g011
Figure 12. The response of entropy after fault-tolerant control for intermittent fault.
Figure 12. The response of entropy after fault-tolerant control for intermittent fault.
Entropy 19 00218 g012
Figure 13. The response of entropy after fault-tolerant control for incipient fault
Figure 13. The response of entropy after fault-tolerant control for incipient fault
Entropy 19 00218 g013
Figure 14. The response of entropy after fault-tolerant control for ramp-like fault.
Figure 14. The response of entropy after fault-tolerant control for ramp-like fault.
Entropy 19 00218 g014

Share and Cite

MDPI and ACS Style

Jin, H.; Guan, Y.; Yao, L. Minimum Entropy Active Fault Tolerant Control of the Non-Gaussian Stochastic Distribution System Subjected to Mean Constraint. Entropy 2017, 19, 218. https://doi.org/10.3390/e19050218

AMA Style

Jin H, Guan Y, Yao L. Minimum Entropy Active Fault Tolerant Control of the Non-Gaussian Stochastic Distribution System Subjected to Mean Constraint. Entropy. 2017; 19(5):218. https://doi.org/10.3390/e19050218

Chicago/Turabian Style

Jin, Haokun, Yacun Guan, and Lina Yao. 2017. "Minimum Entropy Active Fault Tolerant Control of the Non-Gaussian Stochastic Distribution System Subjected to Mean Constraint" Entropy 19, no. 5: 218. https://doi.org/10.3390/e19050218

APA Style

Jin, H., Guan, Y., & Yao, L. (2017). Minimum Entropy Active Fault Tolerant Control of the Non-Gaussian Stochastic Distribution System Subjected to Mean Constraint. Entropy, 19(5), 218. https://doi.org/10.3390/e19050218

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