Next Article in Journal
Wireless Power Transfer Techniques for Implantable Medical Devices: A Review
Next Article in Special Issue
Fault Detection, Isolation, Identification and Recovery (FDIIR) Methods for Automotive Perception Sensors Including a Detailed Literature Survey for Lidar
Previous Article in Journal
Prediction of Body Weight of a Person Lying on a Smart Mat in Nonrestraint and Unconsciousness Conditions
Previous Article in Special Issue
A Reliable Fault Diagnosis Method for a Gearbox System with Varying Rotational Speeds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Towards Simultaneous Actuator and Sensor Faults Estimation for a Class of Takagi-Sugeno Fuzzy Systems: A Twin-Rotor System Application

Institute of Control and Computation Engineering, University of Zielona Góra, ul. Szafrana 2, 65-516 Zielona Góra, Poland
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(12), 3486; https://doi.org/10.3390/s20123486
Submission received: 6 May 2020 / Revised: 12 June 2020 / Accepted: 17 June 2020 / Published: 19 June 2020
(This article belongs to the Special Issue Sensors Fault Diagnosis Trends and Applications)

Abstract

:
The paper is devoted to the problem of estimating simultaneously states, as well as actuator and sensor faults for Takagi–Sugeno systems. The proposed scheme is intended to cope with multiple sensor and actuator faults. To achieve such a goal, the original Takagi–Sugeno system is transformed into a descriptor one containing all state and fault variables within an extended state vector. Moreover, to facilitate the overall design procedure an auxiliary fault vector is introduced. In comparison to the approaches proposed in the literature, a usual restrictive assumption concerning fixed fault rate of change is removed. Finally, the robust convergence of the whole observer is guaranteed by the so-called quadratic boundedness approach which assumes that process and measurement uncertainties are unknown but bounded within an ellipsoid. The last part of the paper portrays an exemplary application concerning a nonlinear twin-rotor system.

1. Introduction

Nowadays, owing to the interest in highly efficient systems, industrial companies permanently extend the number of sensors and actuators being used. Indeed, with the advent of IoT the number of components is systematically proliferating. This is mainly due to the fact of their relative low cost and wide accessibility. Irrespective of such an appealing effect, such a growth can increase the chance that actuator and sensor faults appear simultaneously. Moreover, the probability of multiple actuator and sensor faults is also increased. Thus, fault estimation is an important actor in modern Fault Diagnosis (FD) [1,2,3,4]. Indeed, it may give a knowledge about the presence, location and size of the fault. Such information is also a necessary ingredient of the active Fault-Tolerant Control (FTC) [5,6,7,8], which can accommodate the fault using appropriate control-oriented recovery actions.
There is no doubt that the problem of fault estimation was approached from many different angles. However, most of them focus on estimating either actuator or sensor faults. For a good representative example of recent developments in the area of actuator fault estimation the reader is referred to [9,10]. The case of sensor fault estimation can be realized similarly (cf. Pazera et al. [6] and the references therein).
The number of observers for simultaneous actuator and sensor faults is proliferating as well. However, most of them suffer due to the assumption of a limited fault rate, i.e., the fault derivative (or difference in consecutive discrete samples) is assumed to be close to zero.
Among the existing fault estimation strategies, the crucial attention is focused on the ones based on the sliding-mode observers [11,12], as well as the Kalman filter [13,14] and the related applications.
Taking into account the above discussion, it is also natural that the problem of simultaneous actuator and sensor fault estimation receives growing research attention [8,15,16,17,18]. However, the researches are developed for linear systems while the number of strategies capable of handling some classes of nonlinear systems is rather limited. Indeed, a recent state-of-the-art overview clearly indicates the trends for handling multiple simultaneous actuator and sensor fault estimation for nonlinear systems. Indeed, Gu et al. [19] converted an original Lipschitz system into a Linear Parameter Varying (LPV) one by a suitable reformulation of the Lipschitz property. The resulting design approach provides optimal fault estimates according to H criterion within the frequency domain. Another appealing strategy for Lipschitz class of systems was introduced in Abdollahi [20]. It transforms the system into two subsystems while each of them is affected by either a sensor or an actuator fault, respectively. Finally, two separate Sliding Mode Observers (SMOs) are designed. This scheme suffers from the fact that an asymptotically convergence is guaranteed without considering any robustness to modelling uncertainities and/or external disturbances. A similar SMO was also developed in Hmidi et al. [21] and the authors remove the above drawback by ensuring robustness to bounded disturbances.
An important group of strategies convert the original nonlinear system into an equivalent Takagi-Sugeno one. A representative example of such strategies is presented in Bounemeur et al. [22]. As a result, an adaptive fuzzy estimator is developed capable of handling numerous fault scenarios. Unfortunately, it is devoted to deterministic systems neglecting uncertainty factors. Another strategy is developed in Fu et al. [23] and tackles the systems with switching non-linearities. Unfortunately, it also does not provide appropriate means for settling the robustness issue. Finally, in Shaker [24] a novel multiple integral unknown input observer is developed that is capable of decoupling disturbances. The final group of approaches are the ones for polynomial systems [25,26]. In this case, the design methodology reduces to converting the system by augmenting the state vector with fault variables. This strategy can also cope with process disturbances.
The unappealing feature of the existing approaches to simultaneous actuator and sensor fault estimation is that they do not provide sufficient information about the estimation quality. To handle this issue, two kind of approaches can be utilized: the first one uses post-fault measurements [27] and the second predicts the fault based on its historical values and past measurements [28]. It should be noted that the latter is the only scheme which can be applied in the active FTC. Unfortunately, such schemes struggle with the problem of minimising the effect of fault prediction [9,28].
The scheme proposed in this paper can handle the above mentioned issues, which constitutes the novelty of the proposed approach:
  • the problem of one-step fault prediction and the related fault rate of change is removed by transforming the system into the descriptor Takagi-Sugeno one whose state vector contains both original states and the faults.
  • the original system with the actuator fault with a time-varying distribution matrix is transformed into an equivalent one with a constant distribution matrix and the so-called auxiliary fault vector. This strategy naturally reduces the design conservatives.
  • the effect of external disturbances is tackled with the so-called quadratic boundedness. As a result, the estimation quality can be assessed by the so-called uncertainty intervals of both states and faults. Thus, they can be perceived as possible worst cases of the unknown faults and states.
The paper is organised as follows. In Section 2, a simultaneous estimation of actuator and sensor faults problem is described along with suitable feasibility assumptions. Furthermore, a suitable convergence analysis is performed. Section 3 exhibits an illustrative example pertaining the twin-rotor system. Finally, Section 4 concludes the paper.

2. Observer Design

Let us start from the formulation of a Takagi-Sugeno (T-S) fuzzy system:
x k + 1 = A v k x k + B v k u k + B f v k f ¯ a , k + W 1 w 1 , k = i = 1 M h i v k A i x k + B i u k + B f i f ¯ a , k + W 1 w 1 , k ,
y k = C x k + C f f s , k + W 2 w 2 , k ,
with
h i ( v k ) 0 , i = 1 , , M m , i = 1 M m h i ( v k ) = 1 ,
where k stands for the discrete time and x k X R n , u k R r , y k R m are state, input and output vectors, respectively. Furthermore, f ¯ a , k F a R n a and f s , k F s R n s signify the actuator and sensor fault vectors, respectively, where n a and n s stand for the number of actuator and sensor faults, respectively. Matrices A , B and C are the known state, input and output matrices, respectively, while M m stands for the number of submodels. Thus, C f denotes the sensor fault distribution matrix, where rank ( C f ) = n s and rank ( B f ( v k ) ) = n a satisfy the inequality n a + n s m . This means that it is impossible to estimate more faults than there are measured outputs. Finally, W 1 and W 2 are the distribution matrices of w 1 , k and w 2 , k , which are exogenous disturbance vectors for the process and measurement uncertainties, respectively. The activation functions h i ( · ) depend on the vector of premise variables v k = v k 1 , v k 2 , , v k p T , which is assumed to depend on measurable variables, e.g., system outputs and known inputs [29]. However, an extension towards unmeasurable premise variables is possible via direct applications of the solution proposed, e.g., in Ichalal et al. [30].
Finally, in the remaining part of the paper, the following notation is used X v k = i = 1 M h i v k X i . First, let us start with transforming the state Equation (1) into an equivalent form
x k + 1 = i = 1 M h i v k A i x k + B i u k + B f f a , k + W 1 w 1 , k ,
with an auxiliary matrix B f satisfying rank ( B f ) = n a and an auxiliary actuator fault vector f a , k . Comparing Equation (1) and Equation (4), it can be observed that
B f ( v k ) f ¯ a , k = B f f a , k .
Thus, having f k , the original fault vector can be determined with
f ¯ a , k = ( B f ( v k ) ) B f f a , k .
where stands for the pseudo inverse operator. Thus, the objective of further deliberations is to propose a novel observer capable of providing x k , f a , k and f s , k estimates simultaneously. It should be also noted that the selection of B f is not critical, i.e., it can be set as B f = 1 M i = 1 M B f i , which clearly preserves rank ( B f ) = n a .
The proposed strategy starts with transforming Equation (4) and Equation (2) into an equivalent descriptor form with the following state variable
x ¯ k = x k T f a , k 1 T f s , k T T .
Using the above state vector Equation (7), the system Equation (4) and Equation (2) can be rewritten as:
E x ¯ k + 1 = A ¯ v k x ¯ k + B ¯ v k u k + W ¯ 1 w ¯ k ,
y k = C ¯ x k + W ¯ 2 w ¯ k ,
with:
E = I n B f 0 0 0 0 0 0 0 , A ¯ v k = A v k 0 0 0 0 0 0 0 0 , B ¯ v k = B v k 0 0 , W ¯ 1 = W 1 0 0 0 0 0 , W ¯ 2 = 0 W 2 , w ¯ k = w 1 , k T w 2 , k T T , C ¯ = C 0 C f .
This means that both faults as well as the state of the system are incorporated into a super state vector x ¯ k . Thus, for the purpose of the state estimation obeying Equation (8) and Equation (9), the following observer is proposed:
z k + 1 = N v k z k + M v k u k + L v k y k ,
x ¯ ^ k = z k + T 2 y k ,
where z k R n + n a + n s is the internal state of the estimator while x ¯ ^ k R n + n a + n s stands for the estimate of Equation (7).
It should be noted that the proposed approach eliminates the usual assumption concerning a bounded rate of change of actuator and sensor fault, which increase the conservatives level of the approaches proposed in the literature (see, e.g., Pazera et al. [28] and the references therein.)
Let us start with assuming that there exist matrices T 1 and T 2 such that
T 1 E + T 2 C ¯ = I ,
or
T 1 T 2 E C ¯ = I ,
which yields Equation (10) and Equation (11) design condition
rank E C ¯ = n + n a + n s .
Under the above assumption, it is possible to derive the error of state estimation with Equation (11), i.e.,
e k = x ¯ k x ¯ ^ k = x ¯ k z k T 2 y k = I T 2 C ¯ x ¯ k z k T 2 W ¯ 2 w ¯ k ,
which, by using Equation (12), boils down to
e k = T 1 E x ¯ k z k T 2 W ¯ 2 w ¯ k .
Thus, by substituting Equations (8)–(10), the dynamics of state estimation error obeys
e k + 1 = T 1 E x ¯ k + 1 z k + 1 T 2 W ¯ 2 w ¯ k + 1 = T 1 A ¯ v k x ¯ k + T 1 B ¯ v k u k + T 1 W ¯ 1 w ¯ k N v k z k M v k u k L v k C ¯ x ¯ k L v k W ¯ 2 w ¯ k T 2 W ¯ 2 w ¯ k + 1 = T 1 A ¯ v k L v k C ¯ v k T 1 E x ¯ k + T 1 B ¯ v k M v k u k + N v k e k + T 1 W ¯ 1 L v k W ¯ 2 + N v k T 2 W ¯ 2 w ¯ k T 2 W ¯ 2 w ¯ k + 1 .
From Equation (17), it is evident that the following relations should be satisfied:
T 1 A ¯ v k N v k T 1 E L v k C ¯ = 0 ,
T 1 B ¯ v k M v k = 0 .
Indeed, by satisfying Equation (18), the term related to x ¯ k vanished from Equation (17). Similarly, satisfying Equationv Equation (19) means that Equation (18) is no longer dependent on the system input u k . Applying Equation (12) to Equation (18) gives
T 1 A ¯ v k N v k I T 2 C ¯ L v k C ¯ = 0 ,
or
N v k = T 1 A ¯ v k L v k N v k T 2 C ¯ .
Finally, by defining
K v k = L v k N v k T 2 ,
equality Equation (21) boils down to
N v k = T 1 A ¯ v k K v k C ¯ ,
which makes it possible to transform Equation (17) into
e k + 1 = N v k e k + T 1 W ¯ 1 L v k W ¯ 2 + N v k T 2 W ¯ 2 w ¯ k T 2 W ¯ 2 w ¯ k + 1 = T 1 A ¯ v k K v k C ¯ e k + T 1 W ¯ 1 w ¯ k K v k W ¯ 2 w ¯ k N v k T 2 W ¯ 2 w ¯ k + N v k T 2 W ¯ 2 w ¯ k T 2 W ¯ 2 w ¯ k + 1 = T 1 A ¯ v k K v k C ¯ e k + T 1 W ¯ 1 w ¯ k K v k W ¯ 2 w ¯ k T 2 W ¯ 2 w ¯ k + 1 .
Subsequently, by defining the following super-vector
w ˜ k = w ¯ k w ¯ k + 1 ,
equality Equation (24) can be rewritten into a simpler form
e k + 1 = X v k e k + Z v k w ˜ k ,
with:
X v k = A ˜ v k K v k C ¯ , Z v k = W ˜ 1 K v k W ˜ 2 ,
where:
W ˜ 1 v k = T 1 W ¯ 1 T 2 W ¯ 2 , W ˜ 2 = W ¯ 2 0 , A ˜ v k = T 1 A ¯ v k .
For the purpose of further convergence analysis, let us start with reminding the Finsler’s Lemma [31]:
Lemma 1.
The following expressions are equivalent:
1. 
x ˜ k T Q x ˜ k < 0 , x ˜ { x ˜ R n x | x ˜ 0 , R x ˜ = 0 } ,
2. 
M ¯ R n + n s × m such that Q + M ¯ R + R T M ¯ T < 0 .
Let us also define the Lyapunov function
V k = e k T P e k ,
with P 0 . Furthermore, the estimation error Equation (26) can be rewritten in an alternative form
X v k e k + Z v k w ˜ k e k + 1 = 0 ,
which implies that the following extended-vector can be defined
x ˜ k = e k T w ˜ k T e k + 1 T T ,
and as a consequence the following statement can be formulated
R v k e k w ˜ k e k + 1 = 0 .
Thus, based on Equation (28), Equation (29) and Equation (30), it can be shown that
R v k = X v k Z v k I ,
where
R v k x ˜ k = 0 ,
in order to satisfy Equation (28). The convergence of the proposed observer is to be determined with the so-called Quadratic Boundedness (QB) approach [32]. This technique can be perceived as an extension of the usual Lyapunov approach towards the systems with external bounded disturbances. The usefulness of QB approach was proven in many papers while in Pazera et al. [28] it was proven that the standard H framework can be perceived as a special case of QB. To use the QB approach, it is necessary to assume that w ˜ k is bounded by the following ellipsoid:
E w = { w ˜ k : w ˜ k T Q w w ˜ k 1 } ,
with Q w 0 . Under such assumptions, the following definition can be recalled:
Definition 1.
The system Equation (26) is strictly quadratically bounded for all w ˜ k E w , k 0 , if V k > 1 V k + 1 V k < 0 for any w ˜ k E w .
As shown in Alessandri et al. [32] and Pazera et al. [28], the stability condition associated with
V k + 1 1 α V k α w ˜ k T Q w w ˜ k < 0 ,
with 0 < α < 1 . Based on the above considerations, the following theorem is established:
Theorem 1.
The observer-based system Equation (26) is strictly quadratically bounded for all w ˜ k E w if there exist matrices P 0 , U , N ¯ as well as α ( 0 , 1 ) such that the following holds:
α P P 0 A ˜ T v k U T C ¯ T N ¯ T v k 0 α Q w W ˜ 1 T U T W ˜ 2 T N ¯ T v k U A ˜ v k N ¯ v k C ¯ U W ˜ 1 N ¯ v k W ˜ 2 P U U T 0 ,
Proof. 
Using Equation (34) and setting
Q = α P P 0 0 0 α Q w 0 0 0 P ,
M ¯ = 0 T 0 T U T T .
along with Lemma 1 leads to
α P P 0 X T v k U T 0 α Q w Z T v k U T U X v k U Z v k P U U T 0 .
Setting
U X v k = U A ˜ v k K v k C ¯ = U A ˜ v k N ¯ v k C ¯ ,
U Z v k = U W ˜ 1 K v k W ˜ 2 = U W ˜ 1 N ¯ v k W ˜ 2 ,
into Equation (38) completes the proof. □
Irrespective of the incontestable appeal of the proposed solution Equation (35), its numerical tractability is feasible only if it is transformed to the set of linear matrix inequalities under fixed α :
α P P 0 A ˜ i T U T C ¯ T N ¯ i T 0 α Q w W ˜ 1 T U T W ˜ 2 T N ¯ i T U A ˜ i N ¯ i C ¯ U W ˜ 1 N ¯ i W ˜ 2 P U U T 0 , i = 1 , , M .
Finally, the design procedure boils down to:
Step 0: 
Determine T 1 and T 2 by solving Equation (12).
Step 1: 
Set α , 0 < α < 1 and determine N ¯ i and U by solving Equation (41).
Step 2: 
Calculate:
K i = U 1 N ¯ i ,
N i = T 1 A ˜ i K i C ¯ ,
M i = T 1 B ¯ i ,
L i = K i + N i T 2 .
While the application procedure leads to:
Step 0: 
Substitute k = 0 and set the initial conditions z 0 .
Step 1: 
Calculate z k + 1 and x ^ k with Equations (10) and (11).
Step 2: 
Set k = k + 1 and go to Step 1.
As demonstrated in Pazera et al. [28], the value of α has a direct influence of the convergence rate of the underlying observer. On the other hand, matrix P shapes the ellipsoidal bound of the estimation error:
e k T P e k ζ k , ζ k = 1 + ( 1 α ) k ( 1 e 0 T P e 0 ) ,
which can be directly used to determine the uncertainty intervals of the state variables [28]:
x ¯ ^ k , i ζ k c i T P 1 c i x ¯ k , i x ¯ ^ k , i + ζ k c i T P 1 c i
where c i signifies ith column of an identity matrix of an appropriate dimension.

3. Case Study: Twin-Rotor System

The aim of this section is to validate the performance and correctness of the novel observer. Accordingly, the proposed actuator and sensor fault estimation method has been implemented to the Twin-Rotor System (TRS) [28], illustrated in Figure 1. Note that all equations and a detailed nonlinear model of the TRS can be found in [28]. This model was further transformed into the Takagi-Sugeno form using the dedicated methodology proposed in Rotondo et al. [33]. The matrices shaping the model Equation (1) and Equation (2) are presented in the Appendix A.
The state vector of the system is given by
x = ω v Ω v θ v ω h Ω h θ h T ,
hence the input vector is described by
u = u v u h T ,
where ω v and Ω v stand for the rotational and angular velocities of the main rotor while θ v signifies the pitch angle of the beam. Furthermore, ω h and Ω h are the tail rotor rotational as well as angular velocities whilst θ h is the yaw angle of the beam. Finally, u v and u h stand for actuation due to the main and tail DC motors, respectively. It should be noted that both inputs u v and u h operate in range ( 1 , 1 ) . Moreover, all details and parameters can be found in manufacturer’s user manual [34] and in Pazera et al. [28].
The following fault scenarios have been considered:
f a , 1 , k = 0.1 · u 1 , k 3000 k 6500 0 otherwise ,
f a , 2 , k = 0.1 · u 2 , k 4500 k 8000 0 otherwise ,
f s , 1 , k = y k 2.5 5000 k 7000 0 otherwise ,
f s , 2 , k = y k + 1.6 6000 k 9500 0 otherwise .
According to the approach proposed along with Equation (6):
B f = 1 M m i = 1 M m B f i = 49.4084 · 10 6 66.2480 · 10 6 14.3375 · 10 3 13.2367 · 10 3 303.3781 0 6.9349 · 10 6 3.6497 · 10 6 1.3860 · 10 3 1.0854 · 10 3 0 58.0016 ,
which means that all actuators are considered as possibly faulty and their faults have to be estimated. Whilst the sensor fault distribution matrix is obtained from matrix C (see Appendix A, Equation (A25)) by extracting the rows corresponding to the first and second sensor, respectively and it is given as follows
C f = 0 0 1 0 0 0 0 0 1 0 T .
It can be clearly viewed that within this scenario, both actuator and sensor faults are partially at the same time. Moreover, the considered faults are constant biases, which are either positive or negative. From Equations (50)–(53), it can be clearly viewed that actuator and sensor faults appear partially at the same time. The considered actuator fault Equation (50) pertains a 10% performance decrease of the first actuator, while Equation (51) concerns a 10% performance increase of the second actuator. Moreover, the considered sensor faults Equation (52) and Equation (53) are represented by either positive or negative constant biases equal to 1.6 or 2.5 , respectively. Thus, they denote a significant sensor readings inaccuracies. However, the distribution matrix of sensor fault denotes that they influence the measurements of the angular velocity of the tail rotor Ω h as well as the angular velocity of the main rotor Ω v . Firstly, let us present a comparison between nonlinear estimation approach described in Pazera et al. [28] and the proposed Takagi-Sugeno scheme, which is shown in Figure 2a,b. As it can be observed, the Takagi-Sugeno model-based approach follows the response of the original nonlinear model-based one, and hence, it does not impair the quality of state and fault estimates significantly.
Figure 3a,b illustrate the real actuator faults of the main f a , 1 and tail f a , 2 rotor with blue dash-dotted lines along with their estimates given with red dashed lines. Additionally, the real and estimated values are overbounded by uncertainty intervals which are indicated with black dashed lines. Moreover, the real sensor faults f s , 1 and f s , 2 are presented in Figure 4a,b with blue dash-dotted lines as well as their estimates depicted by red dashed lines along with the uncertainty intervals indicated with black dashed lines. From these figures it is easily seen that the actuator and sensor faults were estimated with a very good accuracy under the measurement uncertainties. Consequently, the system states have been correctly reconstructed.
Figure 5a,b show the real rotational velocities of the main ω v and tail ω h rotor with blue dash-dotted lines. Their measured outputs are presented by light green solid lines along with their estimates given with red dashed lines as well as the uncertainty intervals indicated with black dashed lines. Moreover, Figure 6a,b present the angular velocities of the main Ω v and of the tail Ω h rotor with blue dash-dotted lines whilst their estimates are given with red dashed lines as well as the measured outputs depicted by light green solid lines. Additionally, the uncertainty intervals are indicated with black dashed lines. As can be observed in Figure 6a,b, irrespective of the intermittent fault (marked in light green) the state estimates are very close to the original states. The rotational velocities of the main and tail rotor and the pitch angle of the beam have been accurately estimated even if the actuator and sensor faults occurred simultaneously. This fact is illustrated in Figure 7a,b and Figure 8, which show the evolution of the state estimation error for the above variables. The figures clearly show that the states are properly reconstructed under the actuator glitch along with positive and negative incorrect readings of the sensor.

4. Concluding Remarks

An important aspect of the paper was to develop an observer for the Takagi-Sugeno systems which would be capable of estimating both state with both actuator and sensor. The first part of the proposed strategy was developed in such a way that no adaptive observer is required to obtain the sensor fault estimates. In other words, it is obtained directly by simply transforming an output equation. The second part of the observer concerns the adaptive actuator fault estimation and the state as well. In order to achieve robustness against process and measurement uncertainties and to ensure the stability performance a Quadratic Boundedness strategy is employed. For the purpose of obtaining the observer gain matrices, a set of LMIs is needed to be solved. The other part of the paper presents a performance validation. It has been achieved by implementing the approach to a nonlinear Twin-Rotor MIMO System. The obtained results plainly confirm that the desired properties of the observer have been accomplished. The state has been estimated properly even in the case of both the actuator and sensor fault case. While analysing an actuator fault, it can be concluded that its estimation error is at a very low level in both cases: fault and fault-free as well. Furthermore, the obtained sensor fault estimates quality strongly depends on the quality of the state estimation. The future work will be devoted to applying such an approach to a fault-tolerant control scheme for a T–S fuzzy systems.

Author Contributions

M.P. derived the design procedure for fault estimation. M.W. contributed to the development of the observer structure and determination of the uncertainty interval calculation rules. N.K. contributed to the experiments and described the experimental results. M.B. contributed to the development of the convergence conditions and review of the literature. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by the National Science Centre, Poland under Grant: UMO-2017/27/B/ST7/00620.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

A 1 = 1 9.9709 · 10 3 4.5486 · 10 7 5.1930 · 10 7 1.7300 · 10 9 1.5629 · 10 14 0 9.9418 · 10 1 8.5060 · 10 5 1.0376 · 10 4 5.1856 · 10 7 6.2171 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9988 · 10 1 9.9785 · 10 3 1.7793 · 10 7 0 0 0 2.4758 · 10 2 9.9566 · 10 1 3.4967 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
A 2 = 1 9.9709 · 10 3 4.5486 · 10 7 4.6737 · 10 7 1.5570 · 10 9 1.4066 · 10 14 0 9.9418 · 10 1 8.5060 · 10 5 9.3380 · 10 5 4.6671 · 10 7 5.5954 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9988 · 10 1 9.9785 · 10 3 1.7793 · 10 7 0 0 0 2.4758 · 10 2 9.9566 · 10 1 3.4967 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
A 3 = 1 9.9707 · 10 3 4.5690 · 10 7 5.2162 · 10 7 1.7378 · 10 9 1.5699 · 10 14 0 9.9415 · 10 1 8.5440 · 10 5 1.0422 · 10 4 5.2088 · 10 7 6.2449 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9988 · 10 1 9.9785 · 10 3 1.7793 · 10 7 0 0 0 2.4758 · 10 2 9.9566 · 10 1 3.4967 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
A 4 = 1 9.9707 · 10 3 4.5690 · 10 7 4.6946 · 10 7 1.5640 · 10 9 1.4129 · 10 14 0 9.9415 · 10 1 8.5440 · 10 5 9.3798 · 10 5 4.6879 · 10 7 5.6204 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9988 · 10 1 9.9785 · 10 3 1.7793 · 10 7 0 0 0 2.4758 · 10 2 9.9566 · 10 1 3.4967 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
A 5 = 1 9.9709 · 10 3 4.5486 · 10 7 5.1929 · 10 7 1.7300 · 10 9 1.5629 · 10 14 0 9.9418 · 10 1 8.5060 · 10 5 1.0375 · 10 4 5.1856 · 10 7 6.2171 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9981 · 10 1 9.9783 · 10 3 1.7793 · 10 7 0 0 0 3.7387 · 10 2 9.9559 · 10 1 3.4966 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
A 6 = 1 9.9709 · 10 3 4.5486 · 10 7 4.6736 · 10 7 1.5570 · 10 9 1.4066 · 10 14 0 9.9418 · 10 1 8.5060 · 10 5 9.3378 · 10 5 4.6670 · 10 7 5.5954 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9981 · 10 1 9.9783 · 10 3 1.7793 · 10 7 0 0 0 3.7387 · 10 2 9.9559 · 10 1 3.4966 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
A 7 = 1 9.9707 · 10 3 4.5690 · 10 7 5.2161 · 10 7 1.7377 · 10 9 1.5699 · 10 14 0 9.9415 · 10 1 8.5440 · 10 5 1.0422 · 10 4 5.2088 · 10 7 6.2449 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9981 · 10 1 9.9783 · 10 3 1.7793 · 10 7 0 0 0 3.7387 · 10 2 9.9559 · 10 1 3.4966 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
A 8 = 1 9.9707 · 10 3 4.5690 · 10 7 4.6945 · 10 7 1.5640 · 10 9 1.4129 · 10 14 0 9.9415 · 10 1 8.5440 · 10 5 9.3796 · 10 5 4.6879 · 10 7 5.6204 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9981 · 10 1 9.9783 · 10 3 1.7793 · 10 7 0 0 0 3.7387 · 10 2 9.9559 · 10 1 3.4966 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
A 9 = 1 9.9709 · 10 3 3.1840 · 10 7 5.1930 · 10 7 1.7300 · 10 9 1.5629 · 10 14 0 9.9418 · 10 1 5.9542 · 10 5 1.0376 · 10 4 5.1857 · 10 7 6.2172 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9995 · 10 1 9.9787 · 10 3 1.7793 · 10 7 0 0 0 9.2220 · 10 3 9.9573 · 10 1 3.4968 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
A 10 = 1 9.9709 · 10 3 3.1840 · 10 7 4.6737 · 10 7 1.5570 · 10 9 1.4066 · 10 14 0 9.9418 · 10 1 5.9542 · 10 5 9.3383 · 10 5 4.6671 · 10 7 5.5955 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9995 · 10 1 9.9787 · 10 3 1.7793 · 10 7 0 0 0 9.2220 · 10 3 9.9573 · 10 1 3.4968 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
A 11 = 1 9.9707 · 10 3 3.1983 · 10 7 5.2162 · 10 7 1.7378 · 10 9 1.5699 · 10 14 0 9.9415 · 10 1 5.9808 · 10 5 1.0422 · 10 4 5.2089 · 10 7 6.2450 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9995 · 10 1 9.9787 · 10 3 1.7793 · 10 7 0 0 0 9.2220 · 10 3 9.9573 · 10 1 3.4968 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
A 12 = 1 9.9707 · 10 3 3.1983 · 10 7 4.6946 · 10 7 1.5640 · 10 9 1.4129 · 10 14 0 9.9415 · 10 1 5.9808 · 10 5 9.3800 · 10 5 4.6880 · 10 7 5.6205 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9995 · 10 1 9.9787 · 10 3 1.7793 · 10 7 0 0 0 9.2220 · 10 3 9.9573 · 10 1 3.4968 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
A 13 = 1 9.9709 · 10 3 3.1840 · 10 7 5.1930 · 10 7 1.7300 · 10 9 1.5629 · 10 14 0 9.9418 · 10 1 5.9542 · 10 5 1.0376 · 10 4 5.1857 · 10 7 6.2171 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9989 · 10 1 9.9785 · 10 3 1.7793 · 10 7 0 0 0 2.1851 · 10 2 9.9567 · 10 1 3.4967 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
A 14 = 1 9.9709 · 10 3 3.1840 · 10 7 4.6737 · 10 7 1.5570 · 10 9 1.4066 · 10 14 0 9.9418 · 10 1 5.9542 · 10 5 9.3381 · 10 5 4.6671 · 10 7 5.5954 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9989 · 10 1 9.9785 · 10 3 1.7793 · 10 7 0 0 0 2.1851 · 10 2 9.9567 · 10 1 3.4967 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
A 15 = 1 9.9707 · 10 3 3.1983 · 10 7 5.2162 · 10 7 1.7378 · 10 9 1.5699 · 10 14 0 9.9415 · 10 1 5.9808 · 10 5 1.0422 · 10 4 5.2088 · 10 7 6.2449 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9989 · 10 1 9.9785 · 10 3 1.7793 · 10 7 0 0 0 2.1851 · 10 2 9.9567 · 10 1 3.4967 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
A 16 = 1 9.9707 · 10 3 3.1983 · 10 7 4.6946 · 10 7 1.5640 · 10 9 1.4129 · 10 14 0 9.9415 · 10 1 5.9808 · 10 5 9.3798 · 10 5 4.6880 · 10 7 5.6205 · 10 12 0 0 6.6291 · 10 1 0 0 0 0 0 0 9.9989 · 10 1 9.9785 · 10 3 1.7793 · 10 7 0 0 0 2.1851 · 10 2 9.9567 · 10 1 3.4967 · 10 5 0 0 0 0 0 9.0333 · 10 1 ,
B 1 = 5.7998 · 10 5 5.4643 · 10 5 1.6830 · 10 2 1.0918 · 10 2 3.0338 · 10 2 0 6.9349 · 10 6 3.6497 · 10 6 1.3860 · 10 3 1.0854 · 10 3 0 58.002 , B 2 = 5.7998 · 10 5 5.4643 · 10 5 1.6830 · 10 2 1.0918 · 10 2 3.0338 · 10 2 0 6.9349 · 10 6 3.6497 · 10 6 1.3860 · 10 3 1.0854 · 10 3 0 58.002 ,
B 3 = 5.8257 · 10 5 5.4887 · 10 5 1.6905 · 10 2 1.0967 · 10 2 3.0338 · 10 2 0 6.9349 · 10 6 3.6497 · 10 6 1.3860 · 10 3 1.0854 · 10 3 0 58.002 , B 4 = 5.8257 · 10 5 5.4887 · 10 5 1.6905 · 10 2 1.0967 · 10 2 3.0338 · 10 2 0 6.9349 · 10 6 3.6497 · 10 6 1.3860 · 10 3 1.0854 · 10 3 0 58.002 ,
B 5 = 5.7998 · 10 5 1.2162 · 10 4 1.6830 · 10 2 2.4301 · 10 2 3.0338 · 10 2 0 6.9348 · 10 6 3.6497 · 10 6 1.3859 · 10 3 1.0854 · 10 3 0 58.002 , B 6 = 5.7998 · 10 5 1.2162 · 10 4 1.6830 · 10 2 2.4301 · 10 2 3.0338 · 10 2 0 6.9348 · 10 6 3.6497 · 10 6 1.3859 · 10 3 1.0854 · 10 3 0 58.002 ,
B 7 = 5.8257 · 10 5 1.2217 · 10 4 1.6905 · 10 2 2.4410 · 10 2 3.0338 · 10 2 0 6.9348 · 10 6 3.6497 · 10 6 1.3859 · 10 3 1.0854 · 10 3 0 58.002 , B 8 = 5.8257 · 10 5 1.2217 · 10 4 1.6905 · 10 2 2.4410 · 10 2 3.0338 · 10 2 0 6.9348 · 10 6 3.6497 · 10 6 1.3859 · 10 3 1.0854 · 10 3 0 58.002 ,
B 9 = 4.0598 · 10 5 1.0576 · 10 5 1.1781 · 10 2 2.1132 · 10 3 3.0338 · 10 2 0 6.9350 · 10 6 3.6497 · 10 6 1.3860 · 10 3 1.0854 · 10 3 0 58.002 , B 10 = 4.0598 · 10 5 1.0576 · 10 5 1.1781 · 10 2 2.1132 · 10 3 3.0338 · 10 2 0 6.9350 · 10 6 3.6497 · 10 6 1.3860 · 10 3 1.0854 · 10 3 0 58.002 ,
B 11 = 4.0780 · 10 5 1.0623 · 10 5 1.1834 · 10 2 2.1226 · 10 3 3.0338 · 10 2 0 6.9350 · 10 6 3.6497 · 10 6 1.3860 · 10 3 1.0854 · 10 3 0 58.002 , B 12 = 4.0780 · 10 5 1.0623 · 10 5 1.1834 · 10 2 2.1226 · 10 3 3.0338 · 10 2 0 6.9350 · 10 6 3.6497 · 10 6 1.3860 · 10 3 1.0854 · 10 3 0 58.002 ,
B 13 = 4 0598 · 10 5 7.7558 · 10 5 1.1781 · 10 2 1.5496 · 10 2 3.0338 · 10 2 0 6.9349 · 10 6 3.6497 · 10 6 1.3860 · 10 3 1.0854 · 10 3 0 58.002 , B 14 = 4.0598 · 10 5 7.7558 · 10 5 1.1781 · 10 2 1.5496 · 10 2 3.0338 · 10 2 0 6.9349 · 10 6 3.6497 · 10 6 1.3860 · 10 3 1.0854 · 10 3 0 58.002 ,
B 15 = 4.0780 · 10 5 7.7904 · 10 5 1.1834 · 10 2 1.5566 · 10 2 3.0338 · 10 2 0 6.9349 · 10 6 3.6497 · 10 6 1.3860 · 10 3 1.0854 · 10 3 0 58.002 , B 16 = 4.0780 · 10 5 7.7904 · 10 5 1.1834 · 10 2 1.5566 · 10 2 3.0338 · 10 2 0 6.9349 · 10 6 3.6497 · 10 6 1.3860 · 10 3 1.0854 · 10 3 0 58.002 ,
W 1 = 0.01 I , W 2 = 0.05 I , C = 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 .

References

  1. Rodrigues, M.; Hamdi, H.; Theilliol, D.; Mechmeche, C.; BenHadj Braiek, N. Actuator fault estimation based adaptive polytopic observer for a class of LPV descriptor systems. Int. J. Robust Nonlinear Control 2015, 25, 673–688. [Google Scholar] [CrossRef] [Green Version]
  2. Witczak, M.; Mrugalski, M.; Korbicz, J. Towards robust neural-network-based sensor and actuator fault diagnosis: Application to a tunnel furnace. Neural Process. Lett. 2015, 42, 71–87. [Google Scholar] [CrossRef] [Green Version]
  3. Shao, H.; Jiang, H.; Zhang, X.; Niu, M. Rolling bearing fault diagnosis using an optimization deep belief network. Meas. Sci. Technol. 2015, 26, 115002. [Google Scholar] [CrossRef]
  4. Zhirabok, A.; Shumsky, A. Fault diagnosis in nonlinear hybrid systems. Int. J. Appl. Math. Comput. Sci. 2018, 28, 635–648. [Google Scholar] [CrossRef] [Green Version]
  5. Yao, L.; 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 Pt A, 538–543. [Google Scholar] [CrossRef]
  6. Pazera, M.; Buciakowski, M.; Witczak, M. Robust Multiple Sensor Fault–Tolerant Control For Dynamic Non–Linear Systems: Application To The Aerodynamical Twin–Rotor System. Int. J. Appl. Math. Comput. Sci. 2018, 28, 297–308. [Google Scholar] [CrossRef] [Green Version]
  7. Rotondo, D.; Ponsart, J.C.; Theilliol, D.; Nejjari, F.; Puig, V. A virtual actuator approach for the fault tolerant control of unstable linear systems subject to actuator saturation and fault isolation delay. Annu. Rev. Control 2015, 39, 68–80. [Google Scholar] [CrossRef] [Green Version]
  8. Liu, X.; Gao, Z.; Zhang, A. Robust Fault Tolerant Control for Discrete-Time Dynamic Systems With Applications to Aero Engineering Systems. IEEE Access 2018, 6, 18832–18847. [Google Scholar] [CrossRef]
  9. Lan, J.; Patton, R.J. A decoupling approach to integrated fault-tolerant control for linear systems with unmatched non-differentiable faults. Automatica 2018, 89, 290–299. [Google Scholar] [CrossRef]
  10. 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]
  11. Xia, J.; Guo, Y.; Dai, B.; Zhang, X. Sensor fault diagnosis and system reconfiguration approach for an electric traction PWM rectifier based on sliding mode observer. IEEE Trans. Ind. Appl. 2017, 53, 4768–4778. [Google Scholar] [CrossRef]
  12. Li, J.; Pan, K.; Su, Q. Sensor fault detection and estimation for switched power electronics systems based on sliding mode observer. Appl. Math. Comput. 2019, 353, 282–294. [Google Scholar] [CrossRef]
  13. Zhang, Q. Adaptive Kalman filter for actuator fault diagnosis. Automatica 2018, 93, 333–342. [Google Scholar] [CrossRef]
  14. Gou, L.; Zhou, Z.; Liang, A.; Wang, L.; Liu, Z. Dynamic Threshold Design Based on Kalman Filter in Multiple Fault Diagnosis. In Proceedings of the 2018 37th Chinese Control Conference (CCC), Wuhan, China, 25–27 July 2018; pp. 6105–6109. [Google Scholar]
  15. Youssef, T.; Chadli, M.; Karimi, H.; Wang, R. Actuator and sensor faults estimation based on proportional integral observer for TS fuzzy model. J. Frankl. Inst. 2017, 354, 2524–2542. [Google Scholar] [CrossRef]
  16. Lan, J.; Patton, R. A new strategy for integration of fault estimation within fault-tolerant control. Automatica 2016, 69, 48–59. [Google Scholar] [CrossRef]
  17. Chaves, E.; André, F.; Maitelli, A. Robust Observer-Based Actuator and Sensor Fault Estimation for Discrete-Time Systems. J. Control. Autom. Electr. Syst. 2019, 30, 160–169. [Google Scholar] [CrossRef]
  18. Liu, M.; Cao, X.; Shi, P. Fuzzy-model-based fault-tolerant design for nonlinear stochastic systems against simultaneous sensor and actuator faults. IEEE Trans. Fuzzy Syst. 2013, 21, 789–799. [Google Scholar] [CrossRef]
  19. Gu, Y.; Yang, G. Simultaneous actuator and sensor fault estimation for discrete-time Lipschitz nonlinear systems in finite-frequency domain. Optimal Control Appl. Methods 2018, 39, 410–423. [Google Scholar] [CrossRef]
  20. Abdollahi, M. Simultaneous sensor and actuator fault detection, isolation and estimation of nonlinear Euler-Lagrange systems using sliding mode observers. In Proceedings of the 2018 IEEE Conference on Control Technology and Applications (CCTA), Copenhagen, Denmark, 21–24 August 2018; pp. 392–397. [Google Scholar]
  21. Hmidi, R.; Brahim, A.; Hmida, F.; Sellami, A. Robust fault tolerant control for Lipschitz nonlinear systems with simultaneous actuator and sensor faults. In Proceedings of the 2018 International Conference on Advanced Systems and Electric Technologies (IC_ASET), Hammamet, Tunisia, 22–25 March 2018; pp. 277–283. [Google Scholar]
  22. Bounemeur, A.; Chemachema, M.; Essounbouli, N. Indirect adaptive fuzzy fault-tolerant tracking control for MIMO nonlinear systems with actuator and sensor failures. ISA Trans. 2018, 79, 45–61. [Google Scholar] [CrossRef]
  23. Fu, S.; Qiu, J.; Chen, L.; Mou, S. Adaptive fuzzy observer design for a class of switched nonlinear systems with actuator and sensor faults. IEEE Trans. Fuzzy Syst. 2018, 26, 3730–3742. [Google Scholar] [CrossRef]
  24. Shaker, M. Hybrid approach to design Takagi-Sugeno observer-based FTC for non-linear systems affected by simultaneous time varying actuator and sensor faults. IET Control Theory Appl. 2019, 13, 632–641. [Google Scholar] [CrossRef]
  25. Brizuela Mendoza, J.; Sorcia Vázquez, F.; Guzmán Valdivia, C.; Osorio Sánchez, R.; Martínez García, M. Observer design for sensor and actuator fault estimation applied to polynomial LPV systems: A riderless bicycle study case. Int. J. Syst. Sci. 2018, 49, 2996–3006. [Google Scholar] [CrossRef]
  26. Li, X.; Ahn, C.; Lu, D.; Guo, S. Robust simultaneous fault estimation and nonfragile output feedback fault-tolerant control for Markovian jump systems. IEEE Trans. Syst. Man, Cybern. Syst. 2018, 49, 1769–1776. [Google Scholar] [CrossRef]
  27. Gillijns, S.; De Moor, B. Unbiased minimum-variance input and state estimation for linear discrete-time systems. Automatica 2007, 43, 111–116. [Google Scholar] [CrossRef]
  28. Pazera, M.; Witczak, M. Towards robust simultaneous actuator and sensor fault estimation for a class of nonlinear systems: Design and comparison. IEEE Access 2019, 7, 97143–97158. [Google Scholar] [CrossRef]
  29. Takagi, T.; Sugeno, M. Fuzzy identification of systems and its application to modeling and control. IEEE Trans. Syst. Man Cybern. 1985, 15, 116–132. [Google Scholar] [CrossRef]
  30. 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] [Green Version]
  31. Skelton, R.; Iwasaki, T.; Grigoriadis, D. A Unified Algebraic Approach to Control Design; CRC Press: Boca Raton, FL, USA, 1997. [Google Scholar]
  32. Alessandri, A.; Baglietto, M.; Battistelli, G. Design of state estimators for uncertain linear systems using quadratic boundedness. Automatica 2006, 42, 497–502. [Google Scholar] [CrossRef]
  33. Rotondo, D.; Puig, V.; Nejjari, F.; Witczak, M. Automated generation and comparison of Takagi-Sugeno and polytopic quasi-LPV models. Fuzzy Sets Syst. 2015, 277, 44–64. [Google Scholar] [CrossRef] [Green Version]
  34. INTECO. Two Rotor Aerodynamical System, User’s Manual; INTECO: Kraków, Poland, 2007. [Google Scholar]
Figure 1. A nonlinear twin-rotor MIMO system.
Figure 1. A nonlinear twin-rotor MIMO system.
Sensors 20 03486 g001
Figure 2. Comparison between nonlinear and Takagi-Sugeno response of the system: angular velocity of the main rotor Ω h (a) and angular velocity of the tail rotor Ω v (b).
Figure 2. Comparison between nonlinear and Takagi-Sugeno response of the system: angular velocity of the main rotor Ω h (a) and angular velocity of the tail rotor Ω v (b).
Sensors 20 03486 g002
Figure 3. Actuator faults f a , 1 (a) and f a , 2 (b).
Figure 3. Actuator faults f a , 1 (a) and f a , 2 (b).
Sensors 20 03486 g003
Figure 4. Sensor faults f s , 1 (a) and f s , 2 (b).
Figure 4. Sensor faults f s , 1 (a) and f s , 2 (b).
Sensors 20 03486 g004
Figure 5. State variables ω v (a) and ω h (b).
Figure 5. State variables ω v (a) and ω h (b).
Sensors 20 03486 g005
Figure 6. State variables Ω v (a) and Ω h (b).
Figure 6. State variables Ω v (a) and Ω h (b).
Sensors 20 03486 g006
Figure 7. Evolution of the state estimation error for the rotational velocities of the main rotor (a) and of the tail rotor (b).
Figure 7. Evolution of the state estimation error for the rotational velocities of the main rotor (a) and of the tail rotor (b).
Sensors 20 03486 g007
Figure 8. Evolution of the state estimation error for the pitch angle of the beam.
Figure 8. Evolution of the state estimation error for the pitch angle of the beam.
Sensors 20 03486 g008

Share and Cite

MDPI and ACS Style

Pazera, M.; Witczak, M.; Kukurowski, N.; Buciakowski, M. Towards Simultaneous Actuator and Sensor Faults Estimation for a Class of Takagi-Sugeno Fuzzy Systems: A Twin-Rotor System Application. Sensors 2020, 20, 3486. https://doi.org/10.3390/s20123486

AMA Style

Pazera M, Witczak M, Kukurowski N, Buciakowski M. Towards Simultaneous Actuator and Sensor Faults Estimation for a Class of Takagi-Sugeno Fuzzy Systems: A Twin-Rotor System Application. Sensors. 2020; 20(12):3486. https://doi.org/10.3390/s20123486

Chicago/Turabian Style

Pazera, Marcin, Marcin Witczak, Norbert Kukurowski, and Mariusz Buciakowski. 2020. "Towards Simultaneous Actuator and Sensor Faults Estimation for a Class of Takagi-Sugeno Fuzzy Systems: A Twin-Rotor System Application" Sensors 20, no. 12: 3486. https://doi.org/10.3390/s20123486

APA Style

Pazera, M., Witczak, M., Kukurowski, N., & Buciakowski, M. (2020). Towards Simultaneous Actuator and Sensor Faults Estimation for a Class of Takagi-Sugeno Fuzzy Systems: A Twin-Rotor System Application. Sensors, 20(12), 3486. https://doi.org/10.3390/s20123486

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