Next Article in Journal
Different Estimation Methods for Type I Half-Logistic Topp–Leone Distribution
Next Article in Special Issue
Absence of Global Solutions for a Fractional in Time and Space Shallow-Water System
Previous Article in Journal
Normed Interval Space and Its Topological Structure
Previous Article in Special Issue
q-Analogue of Differential Subordinations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional Order Unknown Inputs Fuzzy Observer for Takagi–Sugeno Systems with Unmeasurable Premise Variables

by
Abdelghani Djeddi
1,†,
Djalel Dib
1,†,
Ahmad Taher Azar
2,3,*,† and
Salem Abdelmalek
4,†
1
Department of Electrical Engineering, Larbi Tebessi University, Tebessa 12002, Algeria
2
College of Engineering, Robotics and Internet-of-Things Lab (RIOTU), Prince Sultan University, Riyadh 12435, Saudi Arabia
3
Faculty of Computers and Artificial Intelligence, Benha University, Benha 13511, Egypt
4
Department of Mathematics, Larbi Tebessi University, Tebessa 12002, Algeria
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2019, 7(10), 984; https://doi.org/10.3390/math7100984
Submission received: 21 August 2019 / Revised: 9 October 2019 / Accepted: 11 October 2019 / Published: 16 October 2019

Abstract

:
This paper presents a new procedure for designing a fractional order unknown input observer (FOUIO) for nonlinear systems represented by a fractional-order Takagi–Sugeno (FOTS) model with unmeasurable premise variables (UPV). Most of the current research on fractional order systems considers models using measurable premise variables (MPV) and therefore cannot be utilized when premise variables are not measurable. The concept of the proposed is to model the FOTS with UPV into an uncertain FOTS model by presenting the estimated state in the model. First, the fractional-order extension of Lyapunov theory is used to investigate the convergence conditions of the FOUIO, and the linear matrix inequalities (LMIs) provide the stability condition. Secondly, performances of the proposed FOUIO are improved by the reduction of bounded external disturbances. Finally, an example is provided to clarify the proposed method. The obtained results show that a good convergence of the outputs and the state estimation errors were observed using the new proposed FOUIO.

1. Introduction

Recently, the interest in fractional derivatives and integral applications, as well as in theoretical and practical works, has grown immensely, see for example [1,2,3,4,5,6]. The main aspects, concept and several applications of fractional calculus are outlined, for example, in [7,8,9,10,11,12,13,14]. This is essentially due to the fact that various physical systems are well described by a fractional order state equation [15,16,17].
Growing applications have attracted interest in studying the state estimation of fractional differential equations in a linear case [18,19,20,21] and in a nonlinear case [22,23,24,25]. It is well known that the study of the problem of stabilization of the fractional order system is particularly important for the synthesis of the observer [26,27,28,29,30,31,32,33,34,35,36].
Takagi–Sugeno (TS) fuzzy models have also attracted attention in recent years. The main feature of this class of nonlinear models is to represent the local dynamics of each fuzzy implication (rule) by linear system models. It has been effectively employed in the implementation of nonlinear systems [37,38,39,40,41]. Takagi–Sugeno models have been broadly utilized to represent nonlinear integer-order systems. However, the fuzzy Takagi–Sugeno scheme remains very efficient for nonlinear fractional order systems (FOS) [42,43,44]. Therefore, the use of fractional-order Takagi–Sugeno (FOTS) models to represent nonlinear FOS will be introduced in this paper. Several approaches confirm that the validity functions of Takagi–Sugeno representation rely on measurable premise variables, whereas various applications, like diagnosis, consider that those variables rely on the input and state variables of the system that are usually immeasurable [45,46,47,48].
Takagi–Sugeno uses premise variables for computing weighting functions. Premise variables can be known (inputs or outputs of the system), or unknown variables taken as the state of the system to be estimated. State variables are usually unmeasurable, but they can be measured by the introduction of sensors, with an additional cost, but the right choice is to estimate the state variables in order to avoid the effects of sensor and shareholder faults that may have appeared on the inputs or outputs of the system considered. This justifies the research works on the state estimation of systems [48,49]. In order to use the state of the system as premise variables, then the states must be estimated, hence the need to synthesize an adequate observer able to estimate the state of the system despite the presence of unknown inputs and disturbances. Hence, it is motivating to deem the common state of unmeasurable variables such as system states. The problem appears especially in the structure of the state of the TS observer.
To implement a fuzzy observer for TS systems with unmeasurable premise variables (UPV), several methods have been evolved, comprising those which take account of analytical advances of an estimation error [50,51,52], and those which use the error description by a TS model with uncertainty or unorganized disruption [49]. The present work presents the Takagi–Sugeno unknown input fractional order observers design for FOTS models with UPV.
The main objective of the current paper is to found new stability and stabilization conditions using FOTS systems with UPV in the continuous case, to implement observers for nonlinear systems. The case where the weighting functions rely on premise variables depending on unmeasurable system states is considered. First, the representation of FOTS systems with UPV and their observers will be considered, which are given under the linear matrix inequalitie (LMI) formulation. Then, an analysis of the stability of the state estimation error studied by using the minimization of the L 2 norm of the transfer from bounded unknown exogenous disturbances to the state estimation error will be established. An application example is designed to demonstrate the performance of the suggested approach.
This paper is organized in the following way. The next section provides some background on the fractional calculus. The FOTS model is presented in Section 3. The main results of the paper, namely the synthesis of the fractional fuzzy observer based unmeasurable premise variables, are presented in Section 4. A new proposed method for unknown input estimation of the fractional order Takagi–Sugeno unknown input observer is given in Section 5. A numerical example is given in Section 6 to demonstrate the efficiency and validity of the proposed approach. Finally, the paper ends with concluding remarks and future perspectives in Section 7.

2. A Brief Introduction to Fractional Calculus

The fractional differo-integral operators are symbolized by a D t α f ( t ) , where a and t, are the bounds of the operation and α R is a generalization of the standard integration and differentiation to an arbitrary order, which can be rational, irrational or even complex. The basic continuous differo-integral operator is given by the following:
a D t α : = d α d t α , for α > 0 , 1 , for α = 0 , a t d τ α , for α < 0 .
In the literature, different definitions can be found concerning fractional order systems. The best most commonly used definitions of fractional order derivatives are:
The Riemann–Liouville (RL) definition [53]:
a D t α f ( t ) = 1 Γ ( m α ) d d t m a t f ( τ ) ( t τ ) 1 ( m α ) d τ .
The Grunwald–Letnikov (GL) definition [53]:
a D t α f ( t ) = lim h 0 1 h α k = 0 ( t a ) / h ( 1 ) k α k f ( t k h ) .
The Caputo definition of the fractional differ-integral operator for the function f ( t ) is adopted in this paper as the Caputo definition allows the initial values of classical integer-order derivatives with clear physical interpretation to be used as follows [51,53]:
a D t α f ( t ) = 1 Γ ( m α ) a t f m ( τ ) ( t τ ) 1 ( m α ) d τ .
In these expressions m 1 < α < m , and Γ ( . ) is the well-known Eulers gamma function:
Γ ( x ) = 0 e t t ( x 1 ) d t , x > 0 .

3. Fractional Order Takagi–Sugeno Model

Consider the following nonlinear system given by [54]:
a D t α x ( t ) = f x , u , y ( t ) = g x , u ,
where x ( t ) R n and α is the fractional order derivative. f and g are nonlinear functions.
Using the well-known transformation by nonlinear sector, the following TS fuzzy system is given [49]:
a D t α x ( t ) = i = 0 M h i ξ ( t ) A i x ( t ) + B i u ( t ) , y ( t ) = i = 0 M h i ξ ( t ) C i x ( t ) + D i u ( t ) ,
where u ( t ) R m is the input vector, x ( t ) R n is the state vector, and y ( t ) R p represents the output vector. A i R n × n , B i R n × n u , C i R n × n y and D i R n y × n u are known matrices. h i ξ ( t ) are the weighting functions relying on the premise variables ξ ( t ) which can be measurable (input or output of the system) or unmeasurable variables (state of the system). It can also be an external signal. These functions confirm the following so-called convex sum property:
i = 0 M h i ξ ( t ) = 1 , 0 h i ξ ( t ) 1 , i 1 , 2 , , M .
In this paper, the target is to model a fractional order fuzzy TS observer with UPV. Thus, the work is dedicated to the problem of state estimation for nonlinear fractional order systems characterized by continuous time TS models, with unknown input u ¯ ( t ) .
Under the hypothesis C 1 = C 2 = = C and D i = 0 , the FOTS model in the presence of unknown inputs and a measurable premise variable can be defined as follows:
a D t α x ( t ) = i = 0 M h i x ( t ) A i x ( t ) + B i u ( t ) + E i u ¯ ( t ) , y ( t ) = C x ( t ) + E u ¯ ( t ) ,
where u ¯ ( t ) R q ( q < p ) is the input vector, and E i R n × n u ¯ and E R n y × n u ¯ are known matrices. This can be rewritten as:
a D t α x ( t ) = i = 0 M h i x ^ ( t ) A i x ( t ) + B i u ( t ) + E i u ¯ ( t ) + h i x ( t ) h i x ^ ( t ) A i x ( t ) + B i u ( t ) + E i u ¯ ( t ) y ( t ) = C x ( t ) + E u ¯ ( t ) .
After rewriting the model (9), we obtain:
a D t α x ( t ) = i = 0 M h i x ^ ( t ) A i x ( t ) + B i u ( t ) + E i u ¯ ( t ) + ω ( t ) , y ( t ) = C x ( t ) + E u ¯ ( t ) .
This form corresponds to a perturbed FOTS model with measurable premise variables (estimated state of the system), where:
ω ( t ) = i = 1 M h i x ( t ) x ^ ( t ) A i x ( t ) + B i u ( t ) + E i u ¯ ( t ) .
This term is considered a global bounded and asymptotically vanishing perturbation.

4. Fractional Order Takagi–Sugeno Unknown Input Observer

The proposed fractional order Takagi–Sugeno fuzzy unknown input observer is given by the following equations:
a D t α x ( t ) = i = 0 M h i x ^ ( t ) N i z ( t ) + G i u ( t ) + L i y ( t ) , x ^ ( t ) = z ( t ) + H y ( t ) .
The state and the output estimation can be defined as:
x ˜ ( t ) = x ( t ) x ^ ( t ) , = x ( t ) z ( t ) + H C x ( t ) + H E u ¯ ( t ) , = P x ( t ) z ( t ) + H E u ¯ ( t ) ,
where
P = I + H C .
Hence, the dynamics of the state estimation error is
t 0 D t α x ^ ( t ) = P t 0 D t α x ( t ) t 0 D t α z ( t ) + H E t 0 D t α u ¯ ( t ) = i = 0 M h i x ^ ( t ) P A i x ( t ) + P B i u ( t ) + P E i u ¯ ( t ) + P ω ( t ) N i z ( t ) G i u ( t ) L i y ( t ) + H E t 0 D t α u ¯ ( t ) ,
replacing y ( t ) and z ( t ) by their respective expressions given by (11) and (13), the state error is given as follows:
t 0 D t α x ˜ ( t ) = i = 0 M h i x ^ ( t ) P A i N i K i C x ( t ) + P B i G i u ( t ) + P E i K i E u ¯ ( t ) + P ω ( t ) + N i e ( t ) + H E t 0 D t α u ¯ ( t ) ,
with K i = N i H + L i .
If the next conditions are satisfied:
H E = 0 ,
N i = P A i K i C ,
P B i = G i ,
P E i = K i E
and
L i = K i N i H .
Then, the dynamics of the state estimation error become:
t 0 D t α x ˜ ( t ) = i = 1 M h i x ^ ( t ) N i x ˜ ( t ) + P ω ( t ) ,
thus showing that the dynamics of the state estimation error is disturbed by ω ( t ) .
To synthesize the matrices of the observer (13), two techniques are proposed.

4.1. First Approach

It is assumed that the term ω ( t ) defined in (12) satisfies the following Lipschitz condition:
ω ( t ) δ x ˜ ( t ) ,
where δ is a positive constant.
Lemma 1
(see [55]). Let M and N be matrices of the appropriate sizes, then the following property holds:
M T N + N T M η M T M + η 1 N T N , η > 0 .
Theorem 1.
A fractional order unknown input observer (13) for system (11) exists if there exists a positive definite matrix X, matrices M i , S, positive scalars η and δ satisfying the following conditions for all i = 1 , , M :
Θ i X + S C X + S C T λ I < 0 ,
S E = 0 ,
X + S C E i = M i E ,
where
Θ i = A i T X + C T S + X + S C A i C T M i T M i C + η δ 2 I .
Then, the fractional order observer (13) is completely defined by:
H = X 1 S ,
K i = X 1 M i ,
N i = I + H C A i K i C ,
L i = K i N i H ,
and
G i = I + H C B i .
Proof of Theorem 1. 
In order to found the existence conditions of the fractional order observer in Theorem 1, Lemma 1 can be introduced:
Considering the following quadratic Lyapunov function:
V ( t ) = x ˜ ( t ) T X x ˜ ( t ) , X = X T > 0 ,
its derivative with regard to time is given by:
t 0 D t α V ( t ) t 0 D t α x ˜ ( t ) T X x ˜ ( t ) ( t ) + x ˜ ( t ) T X t 0 D t α x ˜ ( t ) .
By substituting (24), the dynamic of the quadratic Lyapunov function becomes:
t 0 D t α V ( t ) i = 1 M μ i x ^ ( t ) x ˜ ( t ) T N i T X + X N i x ˜ ( t ) + x ˜ ( t ) T X P ω ( t ) + ω ( t ) T P T X x ˜ ( t ) ,
and when using Lemma 1 and (25), this allows for the following:
x ˜ T X P ω + ω T P T X x ˜ η ω T ω + η 1 x ˜ T X P P T X x ˜ η γ 2 x ˜ T x ˜ + η 1 x ˜ T X P P T X x ˜ .
Substituting (38) in the fractional derivative of the Lyapunov function (36) yields:
t 0 D t α V = i = 1 M μ i x ^ x ˜ T N i T X + X N i + η γ 2 I + η 1 X P P T X x ˜ .
Since the activation functions satisfy condition (8), the fractional derivative of the Lyapunov function is negative if:
N i T X + X N i + η γ 2 I + η 1 X P P T X < 0 .
According to (19), Equation (40) becomes:
P A i K i C T X + X P A i K i C + η δ 2 I + η 1 X P P T X < 0 .
It is noted unfortunately that the matrix inequality (41) gives a disadvantage since it is nonlinear with respect to the variables K i , X and η (more precisely bilinear). A numerical procedure of resolution by linearization is obtained in the following section.
In order to convert these conditions into an LMI formulation, the following change of variables is considered:
M i = X K i
and by using the Schur complement [16], the linear matrix inequality is obtained:
A i T P T X + X P A i C T M i T + η δ 2 I X P P T X η I < 0 .
To satisfy condition (18), the equality can be solved:
X H E = 0 .
Using the change of variable S = X H , linear matrix equality is obtained:
S E = 0 .
The conditions (21) must be satisfied simultaneously, and using the change of variable (42) gives:
X + S C E i = M i E .
Since P = I + H C , replacing P in (43), the matrix inequality of Theorem 1 can be obtained.
The conditions (26)–(28) of Theorem 1 are thus demonstrated. □

4.2. Second Approach

In the case where hypothesis (24) is not satisfied, meaning that the information on its bounded δ is not available, the method established in the previous section cannot be applied.
In this section, another method based on the use of the L 2 approach is proposed.
Theorem 2.
A fractional order unknown input observer (13) for system (11) exists if there exists a positive definite matrix X, matrices M i , S and positive scalars δ ¯ satisfying the following conditions for all i = 1 , , M :
Θ i X + S C X + S C T γ ¯ I < 0 ,
S E = 0
and
X + S C E i = M i E ,
where
Θ i = A i T X + C T S + X + S C A i C T M i T M i C + I .
Then, the fractional order UI observer (13) is completely defined by (30)–(34).
Proof of Theorem 2. 
To prove Theorem 2, the real bounded Lemma 1 [56] is used.
The dynamics of the fractional state estimation error are given by:
t 0 D t α x ˜ ( t ) = i = 1 M h i x ^ ( t ) N i x ˜ ( t ) + P ω ( t ) .
Consider the following Lyapunov quadratic function:
V ( t ) = x ˜ ( t ) T X x ˜ ( t ) , X = X T > 0 .
Its derivative with regard to time is specified by:
t 0 D t α V ( t ) t 0 D t α x ˜ ( t ) T X x ˜ ( t ) + x ˜ ( t ) T X t 0 D t α x ˜ ( t ) .
By substituting (23), the dynamic of the quadratic Lyapunov function is obtained:
t 0 D t α V ( t ) i = 1 M μ i x ^ ( t ) x ˜ ( t ) T N i T X + X N i x ˜ ( t ) + x ˜ ( t ) T X P ω ( t ) + ω ( t ) T P T X x ˜ ( t ) .
In order to mitigate the impact of ω ( t ) on the state estimation error, the L 2 [52] will be used. It can guarantee:
x ˜ ( t ) 2 ω ( t ) 2 < δ , δ > 0 .
The system of the state estimation error is stable and the gain L 2 noted δ of the transfer from ω ( t ) to x ˜ ( t ) is bounded, if:
t 0 D t α V x ˜ ( t ) + x ˜ ( t ) T x ˜ ( t ) δ 2 ω ( t ) T ω ( t ) < 0 .
By substituting t 0 D t α V x ˜ ( t ) , the inequality (56) becomes:
i = 0 M μ i x ^ ( t ) x ˜ ( t ) T N i T X + X N i x ˜ ( t ) + x ˜ ( t ) T X P ω ( t ) + ω ( t ) T P T X x ˜ ( t ) + x ˜ ( t ) T x ˜ ( t ) δ 2 ω ( t ) T ω ( t ) .
The estimation error converges to zero and the gain L 2 of the transfer from ω ( t ) to x ˜ ( t ) is bounded by δ if the following inequality is verified:
i = 1 M μ i x ^ ( t ) N i T X + X N i + I X P P T X δ 2 I < 0 .
The convex sum property of the activation functions makes it possible to write the following sufficient condition:
N i T X + X N i + I X P P T X δ 2 I < 0 .
Using the expression (19) of N i and the changes of variables M i = X K i and δ ¯ = δ 2 (48) becomes:
Θ i X P P T X δ ¯ I < 0 , i = 1 , , M
where
Θ i = A i T P T X + X P A i C T M i T M i C + I .
To satisfy condition (18), the equality can be solved:
X H E = 0 .
Using the change of variable S = X H , the linear matrix equality is obtained:
S E = 0 .
The conditions (21) must be satisfied simultaneously, and by using the change of variable (42) gives:
X + S C E i = M i E .
Since P = I + H C , replacing P in (60), the matrix inequality of Theorem 2 is obtained.
The conditions (47)–(49) of Theorem 2 are thus demonstrated. □

5. Unknown Inputs Estimation

In system (11), the unknown input u ¯ ( t ) appears with the influence matrix:
Φ ( t ) = i = 0 M h i x ^ ( t ) E i E .
For the estimation of the unknown input, it is necessary that the rank of the matrix Φ ( t ) is verified at each time t for the following condition:
rank Φ ( t ) = q ,
where q is the dimension of u ¯ ( t ) . If this condition is satisfied, Φ ( t ) is full-rank column and its pseudo-inverse left Φ 1 ( t ) exists:
Φ ( t ) = Φ T ( t ) Φ ( t ) 1 Φ T ( t ) .
The unknown input can then be calculated according to the state estimated as follows:
u ¯ ^ ( t ) = Φ t 0 D t α x ^ ( t ) i = 0 M h i x ^ ( t ) A i x ^ ( t ) + B i u ( t ) y t C x ^ t ,
under condition (65) the asymptotic convergence from x ^ ( t ) to x ( t ) results in the asymptotic convergence of u ¯ ^ ( t ) to u ¯ ( t ) .

6. Example and Comparisons

To validate the advantages of the proposed fractional order unknown input observer, the system (68) represented by the FOTS model with UPV is considered with α = 0.8 . The state estimation is carried out by means of two fuzzy unknown input observers, the first with integer order and the second one with fractional order. The unknown inputs considered may be noise, faults or modeling uncertainties.

Example and Simulation Results

Consider the FOTS model (8), which is defined as follows:
t 0 D t α x ^ ( t ) = i = 0 M h i x ( t ) A i x ( t ) + B i u ( t ) + F i u ¯ ( t ) , y ( t ) = C x t + G u ¯ ( t ) ,
where: A 1 = 2 1 1 1 3 0 2 1 4 , A 2 = 3 2 2 5 3 0 0.5 0.5 4 , B 1 = 1 0.3 0.5 , B 2 = 0.5 1 0.25 , C = 1 1 1 1 0 1 , F 1 = 0.5 1 0.25 , F 2 = 1 0.52 1 , G = 0.9 0.9 .
The activation functions are chosen in the form:
h 1 x = 1 tanh x 1 2 , h 2 x = 1 h 1 x = 1 + tanh x 1 2 .
Two cases are considered for simulation, the first one in the absence of unknown inputs (unknown inputs are null), and the second one in the presence of unknown inputs. The unknown input considered is accompanied by an additive noise. To have a treatment close to reality, the initial values of the system are chosen non-null, but the initial values of the two unknown input observers are chosen equal to zero.
The outputs and the states of the FOTS system with their estimations given by the fuzzy integer and fractional order unknown input observers, and the unknown inputs and their estimates will be compared and analyzed.
Case 1: Absence of unknown input.
At first, the case of the absence of unknown inputs (unknown inputs are null) will be evaluated. Figure 1 shows the two outputs of the considered FOTS system ( y s 1 , y s 2 ) , the outputs estimated by the FOUIO ( y o 1 , y o 2 ) and the fractional order unknown input observer (FOUIO) ( y f o 1 , y f o 2 ) in the absence of unknown inputs. Figure 2 shows the outputs estimation error (a and b) in the absence of unknown inputs ( y s 1 y o 1 , y s 2 y o 2 ) and ( y s 1 y f o 1 , y s 2 y f o 2 ) . The two Figure 1 and Figure 2 show that the FOUIO gives better output estimation for the considered system. The decreased quality of the output estimation at the moment t = 0 is because of the choice of the initial values.
Figure 3 presents the states of the FOTS system ( x s 1 , x s 2 , x s 3 ) and their estimations given by the FUIO ( x o 1 , x o 2 , x o 3 ) and the FOUIO ( x f o 1 , x f o 2 , x f o 3 ) in the absence of unknown inputs. Figure 4 shows the state estimation errors (a and b) in the absence of unknown inputs ( x s 1 x o 1 , x s 2 x o 2 , x s 3 x o 3 ) and ( x s 1 x f o 1 , x s 2 x f o 2 , x s 3 x f o 3 ) . The two Figure 3 and Figure 4 show that the fuzzy observer with unknown inputs gives a better state estimation of the FOTS system. The decreased quality of the state estimation at the moment t = 0 is because of the choice of the initial values.
Now, the case of the presence of unknown inputs will be evaluated.
Case 2: Presence of unknown input and measurement noise simultaneously.
Figure 5 shows the outputs of the FOTS system ( y s 1 , y s 2 ) , and their estimations given by the FUIO ( y o 1 , y o 2 ) and the fuzzy FOUIO ( y f o 1 , y f o 2 ) in the presence of unknown inputs. Figure 6 shows the outputs estimation error (a and b) in the presence of unknown inputs ( y s 1 y o 1 , y s 2 y o 2 ) and ( y s 1 y f o 1 , y s 2 y f o 2 ) . The two Figure 5 and Figure 6 show that the FOUIO gives a better output estimation for the FOTS system. The decreased quality of the outputs estimation at the moment t = 0 is because of the choice of the initial values.
Figure 7 presents the states of the FOTS system ( x s 1 , x s 2 , x s 3 ) and their estimations given by the FUIO ( x o 1 , x o 2 , x o 3 ) and the FOUIO ( x f o 1 , x f o 2 , x f o 3 ) in the presence of unknown inputs. Figure 8 shows the state estimation errors (a and b) in the presence of unknown inputs ( x s 1 x o 1 , x s 2 x o 2 , x s 3 x o 3 ) and ( x s 1 x f o 1 , x s 2 x f o 2 , x s 3 x f o 3 ) . The two Figure 7 and Figure 8 show that the FOUIO gives a better state estimation for the FOTS system. The decreased quality of the state estimation at the moment t = 0 is because of the choice of the initial values.
Analyzing the convergence conditions of the proposed FOUIO, if the condition (23) on the term ω ( t ) is not satisfied or the value of the constant δ is very important (impossibility of finding a solution with Theorem 1, Theorem 2) offers the possibility of designing the observer with unknown input.
Figure 9 shows the considered unknown input with normal noise (ubar), their estimations given by the FUIO (ubar FUIO), FOUIO (ubar FOUIO) and the unknown input without noise (ubar without noise).
Figure 10 shows the unknown input estimation errors (ubar-ubar FUIO, ubar-ubar FOUIO).
The two Figure 9 and Figure 10 show that the FOUIO gives a better unknown input estimation of the FOTS system, but it cannot be decoupled from the noise. The decreased quality of the unknown input estimation at the moment t = 0 is because of the choice of the initial values.
In the presence of adding random measurement noises bounded by 0.01, the unknown input estimated based on the proposed observer is noisy. Indeed, the presence of measurement noise, at high frequency, decreases the quality of reconstruction of the unknown input.

7. Conclusions

In this paper, a new approach is proposed for designing a fractional order Takagi–Sugeno unknown input observer for a nonlinear system described by FOTS models with UPV. The first step is to rewrite the FOTS system in the form of a disturbed equivalent FOTS and with measurable premise variables. After that, two cases are considered; the first one uses the hypothesis that the perturbation, which appears after rewriting the FOTS model, verifies a Lipschitz condition, while the second one does not use this hypothesis. In this second case, another method is developed and based on an L 2 approach. The convergence conditions of the proposed observers are given in the form of linear matrix inequalities (LMI) that can be easily solved with conventional digital tools.
The obtained results show that a good convergence of the outputs and the state estimation errors is observed using the new proposed FOUIO. The state of the system can be estimated even in the presence of an unknown input varying rapidly since it is totally decoupled from the state. An improvement in the dynamics of the proposed observer is possible by placing the poles.
In future work, it would be interesting to study the decoupling of the noise and the estimation of the unknown inputs using the augmented systems.

Author Contributions

Conceptualization, A.D. and D.D.; methodology, D.D. and A.T.A.; software, A.D., D.D. and S.A.; validation, A.D. and A.T.A.; formal analysis, A.T.A.; investigation, S.A.; resources, S.A.; writing–original draft preparation, A.D. and S.A.; writing–review and editing, A.D., D.D. and A.A.; visualization, D.D. and S.A.; funding acquisition, A.T.A.

Funding

This research was funded by Prince Sultan University, Riyadh, Saudi Arabia.

Acknowledgments

The authors would like to thank Prince Sultan University, Riyadh, Saudi Arabia for supporting and funding this work. Special acknowledgment to Robotics and Internet-of-Things Lab (RIOTU) at Prince Sultan University, Riyadh, SA. Also, the authors wish to acknowledge the editor and anonymous reviewers for their insightful comments, which have improved the quality of this publication.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
FOTSFractional Order Takagi–Sugeno
FOSFractional Order Systems
FUIOFuzzy Unknown Input Observer
FOUIOFractional Order Unknown Input Observer
LMILinear Matrix Inequalities
UPVUnmeasurable Premise Variables
MPVMeasurable Premise Variables

References

  1. Machado, J.T.; Kiryakova, V.; Mainardi, F. Recent history of fractional calculus. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 1140–1153. [Google Scholar] [CrossRef] [Green Version]
  2. Oldham, K.B.; Spanier, J. The Fractional Calculus, Theory and Applications of Differentiation and Integration of Arbitrary Order; Academic Press: New York, NY, USA, 1974. [Google Scholar]
  3. Miller, K.S.; Ross, B. The Fractional Calculus, an Introduction to the Fractional Calculus and Fractional Deferential Equations; John Wiley & Sons Inc.: New York, NY, USA, 1993. [Google Scholar]
  4. George, A.; Argyros, I.K. Intelligent Numerical Methods: Applications to Fractional Calculus; Springer International Publishing: Cham, Switzerland, 2016; Volume 624. [Google Scholar]
  5. Li, W.; Zhao, H.M. Rational function approximation for fractional order differential and integral operators. Acta Autom. Sin. 2011, 37, 999–1005. [Google Scholar]
  6. Xe, D.Y.; Zhao, C.N.; Chen, Y.Q. Robust control for fractional order four-wing hyperchaotic system using LMI. In Proceedings of the IEEE Conference on Mechatronics and Automation, Luoyang, China, 25–28 June 2006; pp. 1043–1048. [Google Scholar]
  7. Li, C.L.; Su, K.L.; Zhang, J.; Wei, D.Q. Rational function approximation for fractional order differential and integral operators. Opt. Int. J. Light Electron Opt. 2013, 124, 5807–5810. [Google Scholar] [CrossRef]
  8. Yuan, J.; Shi, B.; Ji, W.Q. Adaptive sliding mode control of a novel class of fractional chaotic systems. Adv. Math. Phys. 2013, 2013, 6709. [Google Scholar] [CrossRef]
  9. Li, W.; Peng, C.; Wang, Y. Frequency domain subspace identification of commensurate fractional order input time delay systems. Int. J. Control. Autom. Syst. 2011, 9, 310–316. [Google Scholar]
  10. Vinagre, B.M.; Podlubny, I.; Dorcak, L.; Feliu, V. On fractional PID controllers: A frequency domain approach. In Proceedings of the IFAC Workshop on Digital Control: Past, Present and Future of PID Control, Terrasa, Spain, 5–7 April 2000; pp. 53–58. [Google Scholar]
  11. Aldair, A.A.; Wang, W.J. Design of fractional order controller based on evolutionary algorithm for a full vehicle nonlinear active suspension systems. Int. J. Control. Autom. 2010, 3, 33–46. [Google Scholar]
  12. Ostalczyk, P. Discrete Fractional Calculus: Applications in Control and Image Processing. In Series in Computer Vision; World Scientific Publishing Co.: Singapore, 2016; Volume 4. [Google Scholar]
  13. Mozyrska, D.; Wyrwas, M. The Z-transform method and delta type fractional difference operators. Discret. Dyn. Nat. Soc. 2015, 2–3, 1–12. [Google Scholar] [CrossRef]
  14. Das, S. Functional Fractional Calculus for System Identification and Controls; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  15. Ibrir, S. Robust state estimation with q-integral observers. In Proceedings of the American Control Conference, Boston, MA, USA, 30 June–2 July 2004; pp. 3466–3471. [Google Scholar]
  16. Farges, C.; Moze, M.; Sabatier, J. Pseudo-state feedback stabilization of commensurate fractional order systems. Automatica 2010, 46, 1730–1734. [Google Scholar] [CrossRef]
  17. Stanisławski, R.; Rydel, M.; Latawiec, K.J. Modeling of discrete-time fractional-order state space systems using the balanced truncation method. J. Frankl. Inst. 2017, 354, 3008–3020. [Google Scholar] [CrossRef]
  18. Doye, I.; Darouach, M.; Voos, H.; Zasadzinsk, M. Design of unknown input fractional-order observer for fractional-order systems. Int. J. Appl. Math. Comput. Sci. 2013, 23, 491–500. [Google Scholar] [CrossRef]
  19. Wei, Y.H.; Sun, Z.Y.; Hu, Y.S.; Wang, Y. On fractional order adaptive observer. Int. J. Autom. Comput. 2015, 12, 664–670. [Google Scholar] [CrossRef] [Green Version]
  20. Sabatier, J.; Farges, C.; Merveillaut, M.; Feneteau, L. On Observability and Pseudo State Estimation of Fractional Order Systems. Eur. J. Control. 2012, 18, 260–271. [Google Scholar] [CrossRef]
  21. Safarinejadian, B.; Asad, M.; Sha Sadeghi, M. Simultaneous state estimation and parameter identification in linear fractional order systems using colored measurement noise. Int. J. Control. 2016, 89, 2277–2296. [Google Scholar] [CrossRef]
  22. Li, F.; Wu, R.; Liang, S. Observer-based state estimation for non-linear fractional systems. Int. J. Dyn. Syst. Differ. Equ. 2015, 5, 322–335. [Google Scholar] [CrossRef]
  23. Fuli, Z.; Hui, L.; Shouming, Z. State estimation based on fractional order sliding mode observer method for a class of uncertain fractional-order nonlinear systems. Signal Process. 2016, 127, 168–184. [Google Scholar]
  24. Kong, S.; Saif, M.; Liu, B. Observer design for a class of nonlinear fractional-order systems with unknown input. J. Frankl. Inst. 2017, 354, 5503–5518. [Google Scholar] [CrossRef]
  25. Djeghali, N.; Djennoune, S.; Bettayeb, M.; Ghanes, M.; Barbot, J.P. Observation and sliding mode observer for nonlinear fractional-order system with unknown input. ISA Trans. 2016, 63, 1–10. [Google Scholar] [CrossRef]
  26. Ding, B.; Sun, H.; Yang, P. Further studies on LMI based relaxed stabilization conditions for nonlinear systems in Takagi–Sugeno’s form. Automatica 2006, 42, 503–508. [Google Scholar] [CrossRef]
  27. Kruszewski, A.; Wang, R.; Guerra, T.M. Nonquadratic stabilization conditions for a class of uncertain nonlinear discrete time TS fuzzy models: A new approach. IEEE Trans. Autom. Control. 2008, 53, 606–611. [Google Scholar] [CrossRef]
  28. N’Doye, I.; Darouach, M.; Zasadzinski, M.; Radhy, N.E. Robust stabilization of uncertain descriptor fractional-order systems. Automatica 2013, 49, 1907–1913. [Google Scholar] [CrossRef]
  29. Lu, J.G.; Chen, Y.Q. Robust stability and stabilization of fractional-order interval systems with the fractional order α: The case 0 < α < 1. IEEE Trans. Autom. Control. 2010, 55, 152–158. [Google Scholar]
  30. Trigeassou, J.; Maamri, N.; Sabatier, J.; Oustaloup, A. A Lyapunov approach to the stability of fractional differential equations. Signal Process. 2011, 91, 437–445. [Google Scholar] [CrossRef]
  31. Yu, W.; Li, Y.; Wen, G.; Yu, X.; Cao, J. Observer design for tracking consensus in second-order multi-agent systems: Fractional order less than two. IEEE Trans. Autom. Control. 2017, 62, 894–900. [Google Scholar] [CrossRef]
  32. Park, J.H.; Park, T.S.; Kim, S.H. Approximation-Free Output-Feedback Non-Backstepping Controller for Uncertain SISO Nonautonomous Nonlinear Pure-Feedback Systems. Mathematics 2019, 7, 456. [Google Scholar] [CrossRef]
  33. Faieghi, M.; Mashhadi, S.K.M.; Baleanu, D. Sampled-data nonlinear observer design for chaos synchronization: A Lyapunov-based approach. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 2444–2453. [Google Scholar] [CrossRef]
  34. Zhang, X.; Ding, F.; Xu, L.; Alsaedi, A.; Tasawar, H. A Hierarchical Approach for Joint Parameter and State Estimation of a Bilinear System with Autoregressive Noise. Mathematics 2019, 7, 356. [Google Scholar] [CrossRef]
  35. Ibrir, S.; Bettayeb, M. New sufficient conditions for observer-based control of fractional-order uncertain systems. Automatica 2015, 59, 216–223. [Google Scholar] [CrossRef]
  36. Song, C.; Fei, S.; Cao, J.; Huang, C. Robust Synchronization of Fractional-Order Uncertain Chaotic Systems Based on Output Feedback Sliding Mode Control. Mathematics 2019, 7, 599. [Google Scholar] [CrossRef]
  37. Liu, S.; Dong, X.; Zhang, Y. A New State of Charge Estimation Method for Lithium-Ion Battery Based on the Fractional Order Model. IEEE Access 2019, 7, 122949–122954. [Google Scholar] [CrossRef]
  38. Shi, S.L.; Li, J.X.; Fang, Y.M. Fractional-disturbance-observer-based Sliding Mode Control for Fractional Order System with Matched and Mismatched Disturbances. Int. J. Control. Autom. Syst. 2019, 17, 1184–1190. [Google Scholar] [CrossRef]
  39. Trinh, H.; Huong, D.C.; Nahavandi, S. Observer design for positive fractional-order interconnected time-delay systems. Trans. Inst. Meas. Control. Publ. 2018, 41, 378–391. [Google Scholar] [CrossRef]
  40. Dabiri, A.; Butcher, E. Optimal observer-based feedback control for linear fractional-order systems with periodic coefficients. J. Vib. Control. 2019, 25, 1379–1392. [Google Scholar] [CrossRef]
  41. Kong, S.; Saif, M.; Cui, G. Estimation and Fault Diagnosis of Lithium-Ion Batteries: A Fractional-Order System Approach. Math. Probl. Eng. 2018, 8705363, 1–12. [Google Scholar] [CrossRef]
  42. Yang, B.; Yu, T.; Shu, H.; Zhu, D.; Sang, Y.; Jiang, L. Passivity-based fractional-order sliding-mode control design and implementation of grid-connected photovoltaic systems. J. Renew. Sustain. Energy 2018, 10, 43701. [Google Scholar] [CrossRef]
  43. Chadli, M.; Karimi, H.R. Robust observer design for unknown inputs Takagi–Sugeno models. IEEE Trans. Fuzzy Syst. 2013, 21, 158–164. [Google Scholar] [CrossRef]
  44. Liu, S.; Li, X.; Wang, H.; Yan, J. Adaptive fault estimation for T-S fuzzy systems with unmeasurable premise variables. Adv. Differ. Equ. 2018, 105. [Google Scholar] [CrossRef]
  45. Takagi, T.; Sugeno, M. Fuzzy identification of systems and its applications to modeling and control. IEEE Trans. Syst. Man Cybern. Part B Cybern. 1985, SMC-15, 116–132. [Google Scholar] [CrossRef]
  46. Krokavec, D.; Filasova, A. On observer design methods for a class of Takagi Sugeno fuzzy systems. In Proceedings of the Third International Conference on Advanced Information Technologies & Applications, Dubai, UAE, 7–8 November 2014; pp. 279–290. [Google Scholar]
  47. Djeddi, A.; Harkat, M.F.; Soufi, Y. A New Approach for State Estimation of Uncertain Multiple model with Unknown Inputs. Application to Sensor Fault Diagnosis. Mediterr. J. Meas. Control. 2016, 12, 537–545. [Google Scholar]
  48. Chadli, M.; Guerra, T.M. LMI solution for robust static output feedback control of Takagi–Sugeno fuzzy models. IEEE Trans. Fuzzy Syst. 2012, 20, 1060–1065. [Google Scholar] [CrossRef]
  49. Oukacine, S.; Djamah, T.; Djennoune, S.; Mansouri, R.; Bettayeb, M. Multi-model identification of a fractional nonlinear system. IFAC Proc. 2013, 46, 48–53. [Google Scholar] [CrossRef]
  50. Junmin, L.; Yuting, L. Robust Stability and Stabilization of Fractional Order Systems Based on Uncertain Takagi–Sugeno Fuzzy Model With the Fractional Order 1 ≤ ν < 2. ASME J. Comput. Nonlinear Dyn. 2013, 8, 41005. [Google Scholar]
  51. Gao, Z.; Liao, X. Observer-based fuzzy control for nonlinear fractional-order systems via fuzzy T-S models: The 1 < α < 2 case. In Proceedings of the 19th World Congress, The International Federation of Automatic Control, Cape Town, South Africa, 24–29 August 2014. [Google Scholar]
  52. Ichalal, D.; Marx, B.; Mammar, S.; Maquin, D.; Ragot, J. How to cope with unmeasurable premise variables in Takagi–Sugeno observer design: Dynamic extension approach. Eng. Appl. Artif. Intell. 2018, 67, 430–435. [Google Scholar] [CrossRef]
  53. Shantanu, D. Functional Fractional Calculus, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  54. Petras, I. Fractional-Order Nonlinear Systems Modeling, Analysis and Simulation; Higher Education Press: Beijing, China; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  55. Akhenak, A.; Chadli, M.; Ragot, J.; Maquin, D. Estimation of state and unknown inputs of a nonlinear system represented by a multiple model. IFAC Proc. Vol. 2004, 37, 385–390. [Google Scholar] [CrossRef]
  56. Boyd, S.; El Ghaoui, L.; Feron, E.; Balakrishnan, V. Linear Matrix Inequalities in System and Control Theory; SIAM: Philadelphia, PA, USA, 1994. [Google Scholar]
Figure 1. Outputs of the fractional-order Takagi–Sugeno (FOTS) system and its estimation by FOUIO and fractional order unknown input observer (FOUIO) in a free of fault case.
Figure 1. Outputs of the fractional-order Takagi–Sugeno (FOTS) system and its estimation by FOUIO and fractional order unknown input observer (FOUIO) in a free of fault case.
Mathematics 07 00984 g001
Figure 2. Output estimation Error in a free of fault case. (a) Error output estimation of the FOTS system by FUIO in a free of fault case. (b) Error output estimation of the FOTS system by FOUIO in a free of fault case.
Figure 2. Output estimation Error in a free of fault case. (a) Error output estimation of the FOTS system by FUIO in a free of fault case. (b) Error output estimation of the FOTS system by FOUIO in a free of fault case.
Mathematics 07 00984 g002
Figure 3. State of the FOTS system and its estimation by FUIO and FOUIO in a free of fault case.
Figure 3. State of the FOTS system and its estimation by FUIO and FOUIO in a free of fault case.
Mathematics 07 00984 g003
Figure 4. State estimation error in free of fault case. (a) State estimation error between the FOTS system and FUIO in a free of fault case. (b) State estimation error between the FOTS system and FOUIO in a free of fault case.
Figure 4. State estimation error in free of fault case. (a) State estimation error between the FOTS system and FUIO in a free of fault case. (b) State estimation error between the FOTS system and FOUIO in a free of fault case.
Mathematics 07 00984 g004
Figure 5. Outputs for the FOTS system and its estimation by FUIO and FOUIO in a faulty case.
Figure 5. Outputs for the FOTS system and its estimation by FUIO and FOUIO in a faulty case.
Mathematics 07 00984 g005
Figure 6. Output for error estimation faulty case. (a) Outputs for error estimation between the FOTS system and FUIO in a faulty case. (b) Outputs error estimation between the FOTS system and FOUIO in a faulty case.
Figure 6. Output for error estimation faulty case. (a) Outputs for error estimation between the FOTS system and FUIO in a faulty case. (b) Outputs error estimation between the FOTS system and FOUIO in a faulty case.
Mathematics 07 00984 g006
Figure 7. State estimation for the FOTS system and its estimations by FUIO and FOUIO in a faulty case.
Figure 7. State estimation for the FOTS system and its estimations by FUIO and FOUIO in a faulty case.
Mathematics 07 00984 g007
Figure 8. Output estimation error in faulty case. (a) State estimation error between the FOTS system and FUIO in a faulty case. (b) State estimation error between the FOTS system and FOUIO in a faulty case.
Figure 8. Output estimation error in faulty case. (a) State estimation error between the FOTS system and FUIO in a faulty case. (b) State estimation error between the FOTS system and FOUIO in a faulty case.
Mathematics 07 00984 g008
Figure 9. Unknown input and their estimations.
Figure 9. Unknown input and their estimations.
Mathematics 07 00984 g009
Figure 10. Unknown input errors estimation.
Figure 10. Unknown input errors estimation.
Mathematics 07 00984 g010

Share and Cite

MDPI and ACS Style

Djeddi, A.; Dib, D.; Azar, A.T.; Abdelmalek, S. Fractional Order Unknown Inputs Fuzzy Observer for Takagi–Sugeno Systems with Unmeasurable Premise Variables. Mathematics 2019, 7, 984. https://doi.org/10.3390/math7100984

AMA Style

Djeddi A, Dib D, Azar AT, Abdelmalek S. Fractional Order Unknown Inputs Fuzzy Observer for Takagi–Sugeno Systems with Unmeasurable Premise Variables. Mathematics. 2019; 7(10):984. https://doi.org/10.3390/math7100984

Chicago/Turabian Style

Djeddi, Abdelghani, Djalel Dib, Ahmad Taher Azar, and Salem Abdelmalek. 2019. "Fractional Order Unknown Inputs Fuzzy Observer for Takagi–Sugeno Systems with Unmeasurable Premise Variables" Mathematics 7, no. 10: 984. https://doi.org/10.3390/math7100984

APA Style

Djeddi, A., Dib, D., Azar, A. T., & Abdelmalek, S. (2019). Fractional Order Unknown Inputs Fuzzy Observer for Takagi–Sugeno Systems with Unmeasurable Premise Variables. Mathematics, 7(10), 984. https://doi.org/10.3390/math7100984

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