Next Article in Journal
Synthesis of the Hydrophobic Cationic Polyacrylamide (PADD) Initiated by Ultrasonic and its Flocculation and Treatment of Coal Mine Wastewater
Next Article in Special Issue
Improved Genetic Algorithm Tuning Controller Design for Autonomous Hovercraft
Previous Article in Journal / Special Issue
Economic Reliability-Aware MPC-LPV for Operational Management of Flow-Based Water Networks Including Chance-Constraints Programming
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Actuator and System Faults Via an Unknown Input Interval Observer for Takagi–Sugeno Systems

by
Citlaly Martínez-García
1,
Vicenç Puig
2,
Carlos-M. Astorga-Zaragoza
1,*,
Guadalupe Madrigal-Espinosa
1,3 and
Juan Reyes-Reyes
1
1
Department of Electronic Engineering, Tecnológico Nacional de México/CENIDET, Cuernavaca, Morelos 62490, Mexico
2
Advanced Control Systems (SAC) Research Group at Institut de Robótica i Informática Industrial (IRI), Universitat Politécnica de Catalunya-BarcelonaTech (UPC), 08028 Barcelona, Spain
3
National Institute of Electricity and Clean Energy, Cuernavaca, Morelos 62490, Mexico
*
Author to whom correspondence should be addressed.
Processes 2020, 8(1), 61; https://doi.org/10.3390/pr8010061
Submission received: 4 November 2019 / Revised: 20 December 2019 / Accepted: 24 December 2019 / Published: 2 January 2020
(This article belongs to the Special Issue Optimization for Control, Observation and Safety)

Abstract

:
This paper presents a simultaneous state variables and system and actuator fault estimation, based on an unknown input interval observer design for a discrete-time parametric uncertain Takagi–Sugeno system under actuator fault, with disturbances in the process and measurement noise. The observer design is synthesized by considering unknown but bounded process disturbances, output noise, as well as bounded parametric uncertainties. By taking into account these considerations, the upper and lower bounds of the considered faults are estimated. The gain of the unknown input interval observer is computed through a linear matrix inequalities (LMIs) approach using the robust H criteria in order to ensure attenuation of process disturbances and output noise. The interval observer scheme is experimentally evaluated by estimating the upper and lower bounds of a torque load perturbation, a friction parameter and a fault in the input voltage of a permanent magnet DC motor.

1. Introduction

Typically, an observer is an scheme for state estimation through the system input and output measurements. For instance, in [1] a nonlinear observer is applied to estimate the degree of polymerization in a series of polycondensation reactors. However, an observer can be designed for parameter estimation [2], unknown input estimation [3,4] or fault estimation [5,6] among other important applications where it is important to precisely know the actual value of the states, signals or parameters for multiple purposes.
Sometimes there are many technical difficulties in performing an exact estimation of the state, signals or parameters to be estimated. For instance: (i) Model uncertainties, (ii) simplifying assumptions of physical phenomena for modeling, and (iii) complexity reduction of models or the unmeasured disturbances, represent an important source of mismatch between a real process and a mathematical model. In these cases, an approximation of the estimated values can be performed. These approximations can be very useful in many applications where there is not necessary to know the exact value of a variable.
An alternative to estimate unknown variables in processes with uncertain models, interval observers can be used. These observers provide an interval estimation providing a lower and upper bound of the unknown estimated variables. The actual value of the corresponding unmeasured variable located inside the interval defined by these bounds assuming that the uncertainty bounds are known.
Although it is not possible to estimate the exact value of a variable, the information provided by an interval observer can be very useful for several applications. For instance, the authors in [7] propose an interval observer to estimate the lower and upper bounds of vehicle dynamics regardless of the presence of unknown inputs whose bounded interval is also estimated. The authors in [8] design an interval sliding mode observer for sensor fault detection and applied it to an electrical traction device. Another interesting application of interval observers is given in [9], where a trajectory control based on an interval observer is designed for a quadrotor. The interval observer is synthesized by using an uncertain model where all the uncertainties (parameters, disturbance, noise) are unknown but bounded with known bounds.
The main limitation of recent works regarding interval observers is that in most cases, the interval observer design considers linear systems, or a very particular structure of nonlinear systems which sometimes are transformed into linear ones. For instance, the observer in [7] has been designed for switched systems; therefore, its use is limited. In other cases of interval observer designs such as [9], no faults are considered to be estimated or there is a lack of procedure to detect actuator faults [8].
The objective of this paper is to design an interval observer for a wider variety of nonlinear processes by using the Takagi–Sugeno (T–S) approach. Most of the nonlinear models can be adequately transformed into a T–S model (e.g., [10,11]) by using two different methods [12]:
  • The nonlinear sector method, in this case the nonlinear model and its equivalent T–S model have exactly the same behavior. For this reason, this is the method used in this work.
  • The linearization method, in which the equivalent T–S model can be dynamically approximated to the original nonlinear model with a certain accuracy, depending on the design requirements.
Besides, many advantageous opportunities arise when interval observers are designed for processes modeled in T–S form: (i) Pole placement via linear matrix inequalities (LMI) regions is considered to compute the observer gains, in contrast with many nonlinear approaches where the observer gains are heuristically tuned; (ii) a standard methodological procedure can be used to compute the observer gains; (iii) many approaches originally conceived for linear systems can be easily extended to T–S systems. For these reasons, the design of interval observers for T–S systems is a recent and interesting research topic. For example in [13], the authors propose an interval observer for the state estimation of systems modeled in T–S form with parametric uncertainty, disturbances, and measurement noise. However, the work is limited to estimate the unmeasured states. The authors in [14] treat the problem of fault diagnosis of proton exchange membrane (PEM) fuel cells. However, this paper deals with only the case of sensor faults by means of a bank of observers. In [15] a robust fault detection procedure for vehicle lateral dynamics using a switched T–S interval observer is presented. The proposed method is conceived to detect but not to estimate faults.
The main contribution of this paper consists in the design of an interval observer that performs a simultaneous estimation of unmeasured states, actuator and system faults for processes modeled in T–S form with uncertainties. The conditions for the existence of such observers are given. Such conditions guarantee the observer stability and they are proved through a Lyapunov analysis combined with a LMI formulation. The interval observer scheme is experimentally evaluated by estimating the upper and lower bounds of a torque load perturbation, a friction parameter and a fault in the input voltage of a permanent magnet direct-current (DC) motor. These cases are typical faults that, if not detected in time, can become catastrophic failures such as short-circuits or machinery damages due to damaged bearings.

2. Problem Formulation and Preliminaries

Consider the following discrete-time T–S system:
x ( k + 1 ) = i = 1 m ξ i ( ρ ( k ) ) [ ( A i + Δ A i ) x ( k ) + B i u ( k ) ] + E f f ( k ) + G θ ( k ) + E w w ( k ) , y ( k ) = C x ( k ) + E v v ( k ) ,
where x ( k ) R n x , u ( k ) R n u , f ( k ) R n f , θ ( k ) R n θ , w ( k ) R n w and v ( k ) R n v represent the state variable, the input, the actuators fault vector, the unknown parameter, the disturbance and the output noise vector. A i , Δ A i , B i , G and C are matrices of appropriate dimensions. E f , E w and E v are matrices of the coupling distribution. k denotes the k th discrete time instant.
The term ξ i ( ρ ( k ) ) represents the i-th membership function, which is a weighting of the rule i, where i = 1 , 2 , , m . The membership functions are normalized, i.e., they satisfy the following conditions [12,16]:
i = 1 m ξ i ( ρ ( k ) ) = 1 0 ξ i ( ρ ( k ) ) 1 , i = 1 , 2 , , m .
To obtained a simultaneous estimation of parameters and faults, the system (1) is rewritten as follows
x ( k + 1 ) = i = 1 m ξ i ( ρ ( k ) ) ( A i + Δ A i ) x ( k ) + B i u ( k ) + E ¯ f E ¯ ( k ) + E w w ( k ) y ( k ) = C x ( k ) + E v v ( k ) ,
where the vector f E ¯ ( k ) is an augmented one, which is defined by the actuator fault vector f ( k ) and the unknown parameter vector θ ( k ) ; and consequently, the matrix E ¯ contains the fault coupling distribution matrix E f and the parameter matrix G, i.e.,:
E ¯ = E f G , f E ¯ ( k ) = f ( k ) θ ( k ) .
The following considerations are taken into account for the T–S system of the Equation (3):
  • The augmented fault vector f E ¯ ( k ) is defined as:
    f E ¯ ( k + 1 ) = f E ¯ ( k ) + w f E ¯ ( k ) ,
    where w f E ¯ ( k ) is considered as a variation of the actuator fault. Therefore, the estimation of f E ¯ ( k ) is equivalent to the estimation of f ^ ( k ) and θ ^ ( k ) .
  • The perturbation vector w ( k ) is considered unknown but bounded as follows:
    w ¯ ( k ) w ( k ) w ̲ ( k ) .
  • The noise vector v ( k ) is also considered as an unknown but bounded signal, i.e.,:
    | v ( k ) | V ( k ) .
  • The uncertain matrix Δ A i is considered bounded as follows,
    Δ A i ¯ Δ A i Δ A i ̲ .
  • Based on previous assumptions, the estimates to be obtained will be as follows
    x ^ ¯ ( k ) x ( k ) x ^ ̲ ( k ) ,
    f ^ ¯ E ¯ ( k ) f E ¯ ( k ) f ̲ ^ E ¯ ( k ) .
This means that we would get two estimates, i.e., the upper and lower limit of each variable. For that, we consider the following design based on a T–S interval observer.

3. Observer Design

In this section a similar procedure as that in [17] (where no parametric uncertainties nor noise nor disturbances were considered) is presented for the observer design. For this design, first it is considered the output vector at time instant ( k + 1 ) , i.e.,
y ( k + 1 ) = C x ( k + 1 ) + E v v ( k + 1 ) .
Substituting the state equation from system (3), it yields to:
y ( k + 1 ) = C i = 1 m ξ i ( ρ ( k ) ) [ ( A i + Δ A i ) x ( k ) + B i u ( k ) ] + E ¯ f E ¯ ( k ) + E w w ( k ) + E v v ( k + 1 ) .
Next, the following equation can be derived after the pertinent operations
C E ¯ f E ¯ ( k ) = y ( k + 1 ) C i = 1 m ξ i ( ρ ( k ) ) ( A i + Δ A i ) x ( k ) C i = 1 m ξ i ( ρ ( k ) ) B i u ( k ) C E w w ( k ) E v v ( k + 1 ) ,
where it is possible to obtain the fault vector f ( k ) as follows:
f ( k ) = O y ( k + 1 ) C i = 1 m ξ i ( ρ ( k ) ) [ ( A i + Δ A i ) x ( k ) B i u ( k ) ] E w w ( k ) ] E v v ( k + 1 ) ) ,
such that O comes from the following condition, which furthermore must be satisfied for the observer to exist [18]:
r a n k ( C E ¯ f E ¯ ) = r a n k ( E ¯ f E ¯ ) = n θ + n f .
The decoupling is achieved by computing
O = ( C E ¯ f E ¯ ) + ,
such that ( C E ¯ f E ¯ ) + ( C E ¯ f E ¯ ) = I n f is satisfied. Whereas the value of O is obtained as:
O = ( C E ¯ f E ¯ ) T C E ¯ f E ¯ 1 ( C E ¯ f E ¯ ) T .
Replacing fault vector Equation (13) in system Equation (3), the new T–S discrete-time system is obtained:
x ( k + 1 ) = i = 1 m ξ i ( ρ ( k ) ) [ ( A i + Δ A i ) x ( k ) + B i u ( k ) ] + E w w ( k ) + E ¯ O y ( k + 1 ) E ¯ O E v v ( k + 1 ) , y ( k ) = C x ( k ) + E v v ( k )
with
A i = ( I E f O C ) A i , Δ A i = ( I E f O C ) Δ A i ,
B i = ( I E f O C ) B i , E w = ( I E f O C ) E w .
Now, based on (17), the unknown input T–S interval observer structure can be written as follows [19]:
x ^ ̲ ( k + 1 ) = i = 1 m ξ i ( ρ ( k ) ) ( I E ¯ O C ) [ ( A i L ̲ i C ) x ^ ̲ ( k ) + B i u ( k ) + Δ A i x ( k ) + L ̲ i y ( k ) | L ̲ i | E v V ( k ) + E w w ̲ ( k ) + E ¯ O y ( k + 1 ) | L ̲ i | E v V ( k + 1 ) ] , x ^ ¯ ( k + 1 ) = i = 1 m ξ i ( ρ ( k ) ) ( I E ¯ O C ) [ ( A i L ¯ i C ) x ^ ̲ ( k ) + B i u ( k ) + Δ A i x ( k ) + L ¯ i y ( k ) + | L ¯ i | E v V ( k ) + E w w ¯ ( k ) + E ¯ O y ( k + 1 ) + | L ¯ i | E v V ( k + 1 ) ] , f ^ ̲ ( k ) = O [ y ( k + 1 ) C i = 1 m ξ i ( ρ ( k ) ) ( A i x ^ ̲ ( k ) + B i u ( k ) + Δ A i x ( k ) ) C E w w ̲ ( k ) E v V ( k + 1 ) w f ̲ ( k ) ] , f ^ ¯ ( k ) = O [ y ( k + 1 ) C i = 1 m ξ i ( ρ ( k ) ) ( A i x ^ ¯ ( k ) + B i u ( k ) + Δ A i x ( k ) ) C E w w ¯ ( k ) + E v V ( k + 1 ) w f ¯ ( k ) ] ,
with
Δ A i x ( k ) = A ̲ i + x ̲ + A ¯ i + x ̲ A ̲ i x ¯ + + A ¯ i x ¯ ,
Δ A i x ( k ) = A ¯ i + x ¯ + A ̲ i + x ¯ A ¯ i x ̲ + + A ̲ i x ̲ ,
where x ^ ̲ ( k ) and x ^ ¯ ( k ) R x n are the interval estimations of x ( k ) , f ^ ̲ ( k ) and f ^ ¯ ( k ) R s are the interval estimations of f E ¯ ( k ) . L ̲ i and L ¯ i are the observer gains used to compute the upper and lower bounds of the estimated states, faults and parameters, respectively.
The unknown input interval observer can be designer considering (18) in a way that ensures the simultaneous estimation of Equations (8) and (9). The following theorem is introduced to secure the stability analysis and robustness in the presence of unknown entries.
Theorem 1.
Consider the system given by (18) as an interval observer for system (17) for fault and parameter estimation. The observer (18) is stable and robust against the effects of unknown inputs such as bounded disturbances or noise if there exists a symmetric matrix P = P T > 0 , a matrix Q > 0 and the scalars ϵ 1 > 0 , γ > 0 and β > 0 such that:
Q G i , j W i Γ > 0 ,
ϕ i , i < 0 ,
ϕ i , j = I P + γ η 2 I 0 0 0 ( Q G i , j W i Γ ) T 0 γ I ϵ 1 P P H i P Φ 0 0 H i P β 2 I 0 H i T Q T + Φ T Q T 0 Φ P 0 β 2 I Q T Φ ( ) ( ) ( ) ( ) P Q Q T ,
2 m 1 ϕ i , i + ϕ i , j + ϕ j , i < 0 ,
for i , j = 1 , 2 , , m , 1 i j m , i.e., for all subsystems. The observer gains are given by
L ̲ i = Q 1 W ̲ i ,
L ¯ i = Q 1 W ¯ i .
Proof. 
For the stability analysis the following estimation error equations are considered:
e ̲ ( k ) = x ( k ) x ^ ̲ ( k ) ,
e ¯ ( k ) = x ^ ¯ ( k ) x ( k ) .
Substituting the state equation (17) and the estimate state equations (18), (24) and (25) it leads to:
e ̲ ( k + 1 ) = i = 1 m ξ i ( ρ ( k ) ) [ A i + Δ A i x ( k ) + B i u ( k ) ] + E w w ( k ) + E ¯ O y ( k + 1 ) E ¯ O E v v ( k + 1 ) i = 1 m ξ i ( ρ ( k ) ) ( I E ¯ O C ) [ ( A i L ̲ i C ) x ^ ̲ ( k ) + B i u ( k ) + Δ A i x ( k ) + L ̲ i y ( k ) | L ̲ i | E v V ( k ) + E w w ̲ ( k ) + E ¯ O y ( k + 1 ) | L ̲ i | E v V ( k + 1 ) ] ) ,
e ^ ¯ ( k + 1 ) = i = 1 m ξ i ( ρ ( k ) ) ( I E ¯ O C ) [ ( A i L ¯ i C ) x ^ ̲ ( k ) + B i u ( k ) + Δ A i x ( k ) + L ¯ i y ( k ) + | L ¯ i | E v V ( k ) + E w w ¯ ( k ) + E ¯ O y ( k + 1 ) + | L ¯ i | E v V ( k + 1 ) ] i = 1 m ξ i ( ρ ( k ) ) ( A i + Δ A i ) x ( k ) + B i u ( k ) + E w w ( k ) + E ¯ O y ( k + 1 ) E ¯ O E v v ( k + 1 ) ) ,
such that the resulting error equations are the following:
e ̲ ( k + 1 ) = i = 1 m ξ i ( ρ ( k ) ) ( A i L ̲ i C e ̲ ( k ) + Δ A i x ( k ) Δ A i x ( k ) + E w ( w ( k ) w ̲ ( k ) ) + | L ̲ i | E v V ( k ) L ̲ i E v v ( k ) E ¯ O E v v ( k + 1 ) + L ̲ i E v V ( k + 1 ) ] ,
e ¯ ( k + 1 ) = i = 1 m ξ i ( ρ ( k ) ) [ A i L ¯ i C e ¯ ( k ) + Δ A i x ( k ) Δ A i x ( k ) + E w ( w ¯ ( k ) w ( k ) ) + L ¯ i E v v ( k ) + | L ¯ i | E v V ( k ) E ¯ O E v v ( k + 1 ) + L ̲ i E v V ( k + 1 ) ] .
By convenience, the estimation error given by equations (28) and (29) are rewritten as follows
ε ( k + 1 ) = i = 1 m ξ i ( ρ ( k ) ) [ G i ε ( k ) + Θ Δ A + H i δ ( k ) ] + Φ δ ( k + 1 ) ,
ε ( k ) = e ̲ ( k ) e ¯ ( k ) , G i = A i L ̲ i C 0 0 A i L ¯ i C , Θ Δ A = Δ A i x ( k ) Δ A i x ( k ) Δ A i x ( k ) Δ A i x ( k ) ,
H i = E w L ̲ i E v | L ̲ i | E v 0 0 E w L ¯ i E v | L ¯ i | E v ,
Φ = 0 E ¯ O E v E v 0 0 0 E ¯ O E v E v δ ( k ) = w ( k ) w ̲ ( k ) v ( k ) V ( k ) w ¯ ( k ) w ( k ) v ( k ) V ( k ) T .
To show that the observer is stable and robust, the following Lyapunov quadratic function for stability analysis is proposed:
V 1 ( ε ( k ) ) = ε ( k ) T P ε ( k ) > 0 with P = P T > 0 ,
whose increment function corresponds to
Δ V 1 ( ε ( k ) ) = V 1 ( ε ( k + 1 ) ) V 1 ( ε ( k ) ) Δ V 1 ( ε ( k ) ) = ε ( k + 1 ) T P ε ( k + 1 ) ε ( k ) T P ε ( k ) ,
Thus, the the stability condition requires Δ V 1 ( ε ( k ) ) 0 , i.e.,
Δ V 1 ( ε ( k ) ) = i = 1 m j = 1 m ξ i ( ρ ( k ) ) ξ j ( ρ ( k ) ) ( [ G i ε ( k ) + Θ Δ A + H i δ ( k ) ] + Φ δ ( k + 1 ) ) T P , ( [ G i ε ( k ) + Θ Δ A + H i δ ( k ) ] + Φ δ ( k + 1 ) ) ε ( k ) T P ε ( k ) 0 .
If each function is substituted, Equation (33) can be expressed as:
i = 1 m j = 1 m ξ i ( ρ ( k ) ) ξ j ( ρ ( k ) ) ε ( k ) T ( G i T P G j P ) ε ( k ) + Θ Δ A T P Θ Δ A + δ ( k ) T H i T P H j δ ( k ) + δ ( k + 1 ) T Φ T P Φ δ ( k + 1 ) + 2 ε ( k ) T ( G i , j T P Θ Δ A + G i , j T P H j δ ( k ) + G i , j T P Φ δ ( k + 1 ) ) + 2 Θ Δ A T ( P H j δ ( k ) + P Φ δ ( k + 1 ) ) + 2 δ ( k ) T ( H i , j T P Φ δ ( k + 1 ) ) 0 .
Furthermore, for the unknown input T–S interval observer design, the criterion H for the robust estimation problem of T–S system is considered to minimize the effects of noise and disturbance signals:
lim k ε ( k ) = 0 for δ ( k ) = 0 k ,
ε ( k ) 2 < β δ ( k ) 2 for δ ( k ) 0 , ξ ( 0 ) = 0 ,
where β = ψ α correspond to a vector for minimizing the disturbance and noise. The criterion H corresponds to the following function:
ε ( k ) T ε ( k ) β 2 δ ( k ) T δ ( k ) β 2 δ ( k + 1 ) T δ ( k + 1 ) 0 ,
such that the increment of the Lyapunov function results in
V 1 ( ε ( k + 1 ) ) V 1 ( ε ( k ) ) + ε ( k ) T ε ( k ) β 2 δ ( k ) T δ ( k ) β 2 δ ( k + 1 ) T δ ( k + 1 ) 0 .
In addition to considering the stability analysis and robustness, the next condition is considered for the estimation speed Δ V ( ε ( k ) ) ( ϵ 1 P γ ) Δ A i ( k ) for all trajectory, equivalent to
i = 1 m j = 1 m ξ i ( ρ ( k ) ) ξ j ( ρ ( k ) ) ε ( k ) T ( G i T P G j P + I ) ε ( k ) + Θ Δ A T P Θ Δ A + δ ( k ) T H i T P H j δ ( k ) + δ ( k + 1 ) T Φ T P Φ δ ( k + 1 ) + 2 ε ( k ) T ( G i , j T P Θ Δ A + G i , j T P H j δ ( k ) + G i , j T P Φ δ ( k + 1 ) ) + 2 Θ Δ A T ( P H j δ ( k ) + P Φ δ ( k + 1 ) ) + 2 δ ( k ) T ( H i , j T P Φ δ ( k + 1 ) ) β 2 δ ( k ) T δ ( k ) β 2 δ ( k + 1 ) T δ ( k + 1 ) + γ Θ Δ A T Θ Δ A ϵ 1 Θ Δ A T P Θ Δ A 0 ,
whereas in Equation (39) it can be seen that Θ Δ A T P Θ Δ A is a global Lipschitz function such that
f ̲ ( x ̲ , x ¯ ) = ( A ̲ i + A ¯ i + ) x ̲ + A ̲ i x ¯ + + A ¯ i x ¯ ,
f ¯ ( x ̲ , x ¯ ) = ( A ¯ i + A ̲ i + ) x ¯ A ¯ i x ̲ + + A ̲ i x ̲ ,
| f ̲ ( x ̲ , x ¯ ) | Δ A ̲ i + Δ A ¯ i + 2 | x ̲ | + ( Δ A ̲ i 2 + Δ A ¯ i 2 ) | x ¯ | ,
| f ¯ ( x ̲ , x ¯ ) | Δ A ¯ i + Δ A ̲ i + 2 | x ¯ | + ( Δ A ¯ i 2 + Δ A ̲ i 2 ) | x ̲ | ,
and the resulting functions are given by
η = 2 ( Δ A ̲ i + Δ A ¯ i + 2 | + Δ A ̲ i 2 + Δ A ¯ i 2 ) .
Consequently, the resulting incremental Lyapunov function can be rewritten as follows
i = 1 m j = 1 m ξ i ( ρ ( k ) ) ξ j ( ρ ( k ) ) ε ( k ) T ( G i T P G j P + I + γ η 2 I ) ε ( k ) + δ ( k ) T H i T P H j δ ( k ) + δ ( k + 1 ) T Φ T P Φ δ ( k + 1 ) + 2 ε ( k ) T ( G i , j T P Θ Δ A + G i , j T P H j δ ( k ) + G i , j T P Φ δ ( k + 1 ) ) + 2 Θ Δ A T ( P H j δ ( k ) + P Φ δ ( k + 1 ) ) + 2 δ ( k ) T ( H i , j T P Φ δ ( k + 1 ) ) β 2 δ ( k ) T δ ( k ) + γ Θ Δ A T Θ Δ A ϵ 1 Θ Δ A T P Θ Δ A β 2 δ ( k + 1 ) T δ ( k + 1 ) 0 ,
and can be expressed in the following form:
ε ( k ) Θ Δ A δ ( k ) δ ( k + 1 ) T G i T P G j P + η 2 I + I G i , j T P G i T P H i G i T P Φ P G i , j γ I ϵ 1 P P H i P Φ H i P G i H i P H i T P H i β 2 I H i T P Φ Φ P G i Φ P Φ P H i Φ P Φ β 2 I ε ( k ) Θ Δ A δ ( k ) δ ( k + 1 ) 0
To relax the conservatism of (46), the following theorem is considered.
Theorem 2.
There exists a symmetric matrix P > 0 such that [20]
A T P A P < 0 ,
and a matrix G such that the following inequality implies (47)
P A T G T G A P G G T < 0 ,
Consequently, by applying this theorem, inequality (46) is equivalent to
I P + η 2 I 0 0 0 G i , j T Q T 0 γ I ϵ 1 P P H i P Φ 0 0 H i T P β 2 I 0 H i T Q T + Φ T Q T 0 Φ T P 0 β 2 I Φ T Q T Q G i , j 0 Q H i + Q Φ Q Φ P Q Q T 0 ,
such that denoting the inequality (49) as ϕ i , j , it follows
i = 1 m j = 1 m ξ i ( ρ ( k ) ) ξ j ( ρ ( k ) ) ε ( k ) Θ Δ A δ ( k ) δ ( k + 1 ) T ϕ i , j ε ( k ) Θ Δ A δ ( k ) δ ( k + 1 ) 0 .
In the inequality, (50) a bilinearity between the G Q matrices appears as can be been in
A i , j + Δ A ̲ i , j + 0 0 A i , j + Δ A ¯ i , j + G i , j L ̲ i 0 0 L ¯ C 0 0 C Γ Q 0 0 Q Q 0 .
To eliminate the bilinearity that there exists with L ̲ i , L ¯ i and Q matrices, it is possible to use the following change of variables W ̲ i = Q L ̲ i and W ¯ i = Q L ¯ i . Consequently, the following linear inequality is obtained
Q G i , j W i Γ < 0 ,
where W i correspond to W i = W ̲ i 0 0 W ¯ i . Finally, the inequality (21) is the result of using [21], which relaxes the double sum problem. □

4. Simulation Results

4.1. Case Study

A DC motor will be used to illustrate the fault estimation proposed in this paper. The following nonlinear mathematical model represents the dynamics of DC motor [22]:
i ˙ a ( t ) = R a L i a ( t ) K e L v m ( t ) + 1 L u ( t ) , v ˙ m ( t ) = K T J 1 i a ( t ) f r f p v m ( t ) J 1 v m ( t ) T 0 ( t ) T 2 ( t ) J 1 .
where i a ( t ) and v m ( t ) are the armature current and the rotational speed, u ( t ) is the input voltage, T 2 ( t ) and T 0 ( t ) correspond to the load and non-load torque. Table 1 summarizes the model parameter values.
L correspond to the inductance, R a is the armature resistance, K T is the torque-current coefficient, f p is the friction coefficient (due to aerodynamics), K e is the back-emf coefficient, f r is the friction coefficient (due to the bearing lubrication condition) and J 1 is the normalized inertial moment of the rotor.
The nonlinear model (53) can be transformed first into a continuous T–S representation (3) considering the following assumptions:
Assumption 1.
The torque T 0 and T 2 are considered to be unknown. Therefore, it is necessary to decouple their effect.
Assumption 2.
The rotational speed v m is a measurement and is considered as the scheduling parameter.
Assumption 3.
The armature current i a is the measured output.
Consequently, by considering that the rotational speed is scheduling variable ρ ( k ) = v m ( k ) = x 2 ( k ) varying in the interval ρ ( k ) [ ρ ̲ ρ ¯ ] , being ρ ̲ = 100 and ρ ¯ = 300 ] the minimal and maximal rotational speeds. The results T–S representation (3) has the following matrices:
A 1 = R a L K e L K T J 1 f r + f p ρ ¯ J 1 , A 2 = R a L K e L K T J 1 f r + f p ρ ̲ J 1 , B = 1 L 0 ,
E ¯ = 1 L 0 0 0 1 J 1 1 J 1 , E w = 1 L 0 , E v = 0.98 , C = 1 0 0 1 ,
Δ A i = 0.01 A i ; and f E ¯ = f ( k ) T 0 ( k ) T 2 ( k ) T .
The previous continuous-time T–S model can be expressed in discrete time with a sampling time T s = 1 . The resulting matrices are
A 1 = 0.2471 0.0088 8.3671 0.9178 , A 2 = 0.2480 0.0086 8.1747 0.8820 , B = 0.6591 5.9246 , Δ A 1 + = 0.0024 0.00008 0.0836 0.0091 ,
Δ A 2 + = 0.0024 0.00008 0.0817 0.0088 , Δ A 1 = 0.0024 0.00008 0.0836 0.0091 , Δ A 2 = 0.0024 0.00008 0.0817 0.0088 ,
Δ A ̲ 1 = 0.0024 0 0.0836 0.0091 , Δ A ̲ 2 = 0.0024 0 0.0817 0.0088 , Δ A ̲ 1 + = 0 0.8823 0 0 × 10 4 ,
Δ A ̲ 2 + = 0 0.8620 0 0 × 10 4 , Δ A ¯ 1 = 0 0.8823 0 0 × 10 4 , Δ A ¯ 2 = 0 0.8620 0 0 × 10 4 ,
Δ A ¯ 1 + = 0.0024 0 0.0836 0.0091 × 10 4 , Δ A ¯ 2 + = 0.0024 0 0.0817 0.0088 × 10 4 ,
E ¯ = 0.0117 0 0 0 1.4949 1.4949 × 10 2 , E w = 1.1764 0 , E v = 0.08 0.08 .
The solution of LMIs (19)–(21) of Theorem 1 (considering γ = 15 , ϵ 1 = 36.46 , η = 0.337 and β = 5.6214 4.1365 ) lead to the following solution
P = 3.6658 0.8521 0.8521 6.1022 , Q = 43.6972 14.9802 14.9802 19.6044 ,
L ̲ 1 = 0.0099 0.0111 0.1453 0.2328 , L ̲ 2 = 0.0085 0.0245 0.1847 0.3049 ,
L ¯ 1 = 0.1301 0.1171 0.8238 0.8665 × 10 13 , L ¯ 2 = 0.07489 0.2300 0.0382 0.2643 × 10 13 ,
The initial conditions for the T–S unknown input interval observer are x ^ ̲ ( 0 ) = 18 250 T , x ¯ ^ ( 0 ) = 5 150 T , f ^ ̲ ( 0 ) = 0.01 0.01 0.001 , f ¯ ^ ( 0 ) = 0.001 0.02 0 . Additionally, the system disturbance system and output noise is considered to be bounded with the following bounds: 0.98 w ( k ) 0.98 and | v ( k ) | 0.8 .

4.2. Experimental Tests

Two scenarios are considered for the evaluation of the interval observer. The armature current i a ( t ) , measurable via an oscilloscope, and the rotational speed v m ( t ) of the motor, measurable via an incremental encoder associated with an FPGA myRIO-1900 board of National Instruments is used for implementing the proposed approach.
In the evaluation tests, the laboratory prototype shown in Figure 1 is used. This prototype consists of a DC motor available at the TecNM/CENIDET in Mexico (1) coupled to a bearing train (2), and an incremental encoder (3) through a band, whose mathematical model is presented in Equation (53). The results show the good performance of the interval observer in the event of an actuator fault.
In the first evaluation test, the DC motor is powered with 14 V at time instant 390 s, an abrupt fault, almost instantaneous, is introduced in the motor supply voltage via a programmable testing power source. The fault in the motor input produces a decrease of 3.5 V.
Figure 2 shows the measurement of the armature current and the limits (upper and lower) estimated by the interval observer. It can be seen in the figure that the current and limits slightly change their value in the presence of the fault.
Figure 3 shows that the motor speed signal and the limits (upper and lower), estimated by the interval observer (18), present a fairly close dynamic behavior and the speed is always kept within the limits. When the fault disappears, the speed signal recovers its nominal value in approximately 120 s, with the dynamics of the motor coupled to a bearing train.
Figure 4 shows the estimated limits for the input fault, which corresponds to a change in the motor supply voltage. The limits are kept at a value of zero in the absence of failure and change their value when the fault is present.
Figure 5 and Figure 6 show the estimated values of parameters T 0 and T 2 , of the parameter vector θ ( k ) . It can be observed that these parameters remain relatively constant (around 0 and 0.5, respectively) and in the presence of the fault their values are modified. When the fault disappears, they converge again to their initial values.
Figure 7 shows the dynamic behavior of the membership functions, which meet the conditions described in Equation (33).
In the second evaluation test, the DC motor is powered with 15 V at time instant 420 s. An intermittent fault occurs in the supply voltage to the DC motor, caused by interruptions in the connection of the power supply.
Figure 8 shows the dynamic behavior of the armature current signal. Figure 9 shows the variations of the motor rotational speed signal. The current signal and the speed signal, both measurable, are maintained within their respective estimated intervals, in the presence of a fault.
Figure 10 shows the evolution of the estimated bounds for the input fault, which corresponds to change in the the voltage of the motor power supply. The limits are kept at a value centered around zero in the absence of fault and change their value when the fault is present.
Figure 11 and Figure 12 show the estimated values of parameters T 0 and T 2 , of the parameter vector θ ( k ) .
Figure 13 shows the dynamic behavior of the membership functions, which meet the conditions described in Equation (2).

5. Conclusions

A discrete-time unknown-input interval observer is proposed for a system modeled in T–S form with uncertainties. This observer allows the simultaneous estimation of unmeasured states, actuator and system faults despite disturbances and measurement noise. The structure of the proposed discrete-time T–S model has four additional terms: Three terms in the dynamic structure corresponding to the fault, disturbance and parametric uncertainty, and an additive noise term in the output (measurement noise). The conditions for the existence of the observer are formally given to guarantee the observer stability. Such conditions are derived through a Lyapunov analysis combined with a LMI formulation. The proposed discrete-time interval observer approach is experimentally evaluated by estimating the upper and lower bounds of a torque load perturbation, a friction parameter and a fault in the input voltage, in a permanent magnet DC motor.
The main advantage of the proposed T–S interval observer with respect to Kalman or Luenberger-like observers is that a great amount of nonlinear models can be transformed into the Takagi–Sugeno form, with a consequent benefit of preserving the model dynamics. This feature allows us to use this observer for a great number of nonlinear systems, in contrast with Kalman or Luenberger-like observes which requires linear or linearized systems to be implemented.

Author Contributions

All the authors have equally contributed. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been partially funded by the Spanish State Research Agency (AEI) and the European Regional Development Fund (ERFD) through the projects SCAV (ref. MINECO DPI2017-88403-R) and also by EU INTERREG POCTEFA (2014-2020) EFA 153/16 SMART.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ling, C.; Kravaris, C. State observer design for monitoring the degree of polymerization in a series of melt polycondensation reactors. Processes 2016, 4, 4. [Google Scholar] [CrossRef] [Green Version]
  2. Nagy-Kiss, A.M.; Schutz, G.; Ragot, J. Parameter estimation for uncertain systems based on fault diagnosis using Takagi–Sugeno model. ISA Trans. 2015, 56, 65–74. [Google Scholar] [CrossRef] [PubMed]
  3. Youssef, T.; Chadli, M.; Karimi, H.R.; Zelmat, M. Design of unknown inputs proportional integral observers for TS fuzzy models. Neurocomputing 2014, 123, 156–165. [Google Scholar] [CrossRef]
  4. Peng, C.-C. Nonlinear Integral Type Observer Design for State Estimation and Unknown Input Reconstruction. Appl. Sci. 2017, 7, 67. [Google Scholar] [CrossRef] [Green Version]
  5. Zhang, K.; Jiang, B.; Shi, P.; Xu, J. Fault estimation observer design for discrete-time systems in finite-frequency domain. Int. J. Robust Nonlinear Control. 2015, 25, 1379–1398. [Google Scholar] [CrossRef]
  6. Van Nguyen, T.; Ha, C. Sensor Fault-Tolerant Control Design for Mini Motion Package Electro-Hydraulic Actuator. Processes 2019, 7, 89. [Google Scholar] [CrossRef] [Green Version]
  7. Ifqir, S.; Ichalal, D.; Oufroukh, N.A.; Mammar, S. Robust interval observer for switched systems with unknown inputs: Application to vehicle dynamics estimation. Eur. J. Control 2018, 44, 3–14. [Google Scholar] [CrossRef]
  8. Zhang, K.; Jiang, B.; Yan, X.-G.; Shen, J. Interval Sliding Mode Observer Based Incipient Sensor Fault Detection With Application to a Traction Device in China Railway High-Speed. IEEE Trans. Veh. Technol. 2019, 68, 2585–2597. [Google Scholar] [CrossRef]
  9. Abadi, A.; El Amraoui, A.; Mekki, H.; Ramdani, N. Guaranteed trajectory tracking control based on interval observer for quadrotors. Int. J. Control 2019, 1–17. [Google Scholar] [CrossRef]
  10. Chang, Y.-C.; Tsai, C.-T.; Lu, Y.-L. Current Control of the Permanent-Magnet Synchronous Generator Using Interval Type-2 TS Fuzzy Systems. Energies 2019, 12, 2953. [Google Scholar] [CrossRef] [Green Version]
  11. Liu, F.; Li, R.; Dreglea, A. Wind Speed and Power Ultra Short-Term Robust Forecasting Based on Takagi–Sugeno Fuzzy Model. Energies 2019, 12, 3551. [Google Scholar] [CrossRef] [Green Version]
  12. Lendek, Z.; Guerra, T.M.; Babuska, R.; De Schutter, B. Stability Analysis and Nonlinear Observer Design Using Takagi-Sugeno Fuzzy Models; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  13. Li, J.; Wang, Z.; Shen, Y.; Wang, Y. Interval Observer Design for Discrete-Time Uncertain Takagi—Sugeno Fuzzy Systems. IEEE Trans. Fuzzy Syst. 2019, 27, 816–823. [Google Scholar] [CrossRef]
  14. Rotondo, D.; Fernandez-Canti, R.M.; Tornil-Sin, S.; Blesa, J.; Puig, V. Robust fault diagnosis of proton exchange membrane fuel cells using a Takagi-Sugeno interval observer approach. Int. J. Hydrogen Energy 2016, 41, 2875–2886. [Google Scholar] [CrossRef] [Green Version]
  15. Ifqir, S.; Ichalal, D.; Oufroukh, N.A.; Mammar, S. Adaptive Threshold Generation for Vehicle Fault Detection using Switched TS Interval observers. IEEE Trans. Ind. Electron. 2019. [Google Scholar] [CrossRef]
  16. Ohtake, H.; Tanaka, K.; Wang, H.O. Fuzzy modeling via sector nonlinearity concept. Integr. Comput. Aided Eng. 2003, 10, 333–341. [Google Scholar] [CrossRef]
  17. Rotondo, D.; Witczak, M.; Puig, V.; Nejjari, F.; Pazera, M. Robust unknown input observer for state and fault estimation in discrete-time Takagi–Sugeno systems. Int. J. Syst. Sci. 2016, 47, 3409–3424. [Google Scholar] [CrossRef] [Green Version]
  18. Hui, S.; Żak, S. Observer design for systems with unknown inputs. Int. J. Appl. Math. Comput. Sci. 2005, 15, 431–446. [Google Scholar]
  19. Efimov, D.; Raïssi, T.; Perruquetti, W.; Zolghadri, A. Estimation and control of discrete-time LPV systems using interval observers. In Proceedings of the 52nd IEEE Conference on Decision and Control, Florence, Italy, 10–13 December 2013; pp. 5036–5041. [Google Scholar]
  20. De Oliveira, M.C.; Bernussou, J.; Geromel, J.C. A new discrete-time robust stability condition. Syst. Control Lett. 1999, 37, 261–265. [Google Scholar] [CrossRef]
  21. Tuan, H.D.; Apkarian, P.; Narikiyo, T.; Yamamoto, Y. Parameterized linear matrix inequality techniques in fuzzy control system design. IEEE Trans. Fuzzy Syst. 2001, 9, 324–332. [Google Scholar] [CrossRef] [Green Version]
  22. Liu, X.-Q.; Zhang, H.-Y.; Liu, J.; Yang, J. Fault detection and diagnosis of permanent-magnet DC motor based on parameter estimation and neural network. IEEE Trans. Ind. Electron. 2000, 47, 1021–1030. [Google Scholar]
Figure 1. Laboratory prototype.
Figure 1. Laboratory prototype.
Processes 08 00061 g001
Figure 2. Measurement of armature current x 1 ( k ) and estimation of interval bounds.
Figure 2. Measurement of armature current x 1 ( k ) and estimation of interval bounds.
Processes 08 00061 g002
Figure 3. Measurement of rotational speed x 2 ( k ) and estimation of interval bounds.
Figure 3. Measurement of rotational speed x 2 ( k ) and estimation of interval bounds.
Processes 08 00061 g003
Figure 4. Estimation of the bounds for the input voltage fault.
Figure 4. Estimation of the bounds for the input voltage fault.
Processes 08 00061 g004
Figure 5. Estimation of T 0 ( k ) .
Figure 5. Estimation of T 0 ( k ) .
Processes 08 00061 g005
Figure 6. Estimation of T 2 ( k ) .
Figure 6. Estimation of T 2 ( k ) .
Processes 08 00061 g006
Figure 7. Membership functions.
Figure 7. Membership functions.
Processes 08 00061 g007
Figure 8. Measurement of armature current x 1 ( k ) and estimation of interval bounds.
Figure 8. Measurement of armature current x 1 ( k ) and estimation of interval bounds.
Processes 08 00061 g008
Figure 9. Measurement of rotational speed x 2 ( k ) and estimation of interval bounds.
Figure 9. Measurement of rotational speed x 2 ( k ) and estimation of interval bounds.
Processes 08 00061 g009
Figure 10. Estimation of the limits of input voltage fault.
Figure 10. Estimation of the limits of input voltage fault.
Processes 08 00061 g010
Figure 11. Upper and lower bound estimations of T 0 ( k ) .
Figure 11. Upper and lower bound estimations of T 0 ( k ) .
Processes 08 00061 g011
Figure 12. Upper and lower bound estimations of T 2 ( k ) .
Figure 12. Upper and lower bound estimations of T 2 ( k ) .
Processes 08 00061 g012
Figure 13. Membership functions.
Figure 13. Membership functions.
Processes 08 00061 g013
Table 1. Parameters of a DC motor.
Table 1. Parameters of a DC motor.
ParameterValue
L 850 × 10 3 H
R a 1.02 Ω
K T 0.1 N·m/A
f p 0.000000075 N/rpm 2
f r 0.0000035 N/rpm
K e 0.0134 V/rpm
J 1 0.00668933 N·m·s

Share and Cite

MDPI and ACS Style

Martínez-García, C.; Puig, V.; Astorga-Zaragoza, C.-M.; Madrigal-Espinosa, G.; Reyes-Reyes, J. Estimation of Actuator and System Faults Via an Unknown Input Interval Observer for Takagi–Sugeno Systems. Processes 2020, 8, 61. https://doi.org/10.3390/pr8010061

AMA Style

Martínez-García C, Puig V, Astorga-Zaragoza C-M, Madrigal-Espinosa G, Reyes-Reyes J. Estimation of Actuator and System Faults Via an Unknown Input Interval Observer for Takagi–Sugeno Systems. Processes. 2020; 8(1):61. https://doi.org/10.3390/pr8010061

Chicago/Turabian Style

Martínez-García, Citlaly, Vicenç Puig, Carlos-M. Astorga-Zaragoza, Guadalupe Madrigal-Espinosa, and Juan Reyes-Reyes. 2020. "Estimation of Actuator and System Faults Via an Unknown Input Interval Observer for Takagi–Sugeno Systems" Processes 8, no. 1: 61. https://doi.org/10.3390/pr8010061

APA Style

Martínez-García, C., Puig, V., Astorga-Zaragoza, C. -M., Madrigal-Espinosa, G., & Reyes-Reyes, J. (2020). Estimation of Actuator and System Faults Via an Unknown Input Interval Observer for Takagi–Sugeno Systems. Processes, 8(1), 61. https://doi.org/10.3390/pr8010061

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