Document not found! Please try again

Fault Diagnosis and Fault Tolerant Control for Non-Gaussian Singular

0 downloads 0 Views 2MB Size Report
Jul 28, 2016 - Non-Gaussian Singular Time-Delayed Stochastic Distribution ... the feedback of distribution tracking error, fault, and disturbance estimation to ...
Hindawi Publishing Corporation Journal of Control Science and Engineering Volume 2016, Article ID 9563481, 11 pages http://dx.doi.org/10.1155/2016/9563481

Research Article Fault Diagnosis and Fault Tolerant Control for Non-Gaussian Singular Time-Delayed Stochastic Distribution Systems with Disturbance Based on the Rational Square-Root Model Yuancheng Sun and Zhanhong Liang School of Electrical Engineering, Zhengzhou University, Zhengzhou, Henan 450001, China Correspondence should be addressed to Zhanhong Liang; [email protected] Received 23 March 2016; Revised 14 July 2016; Accepted 28 July 2016 Academic Editor: Manuel Pineda-Sanchez Copyright © 2016 Y. Sun and Z. Liang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. For the non-Gaussian singular time-delayed stochastic distribution control (SDC) system with unknown external disturbance where the output probability density function (PDF) is approximated by the rational square-root B-spline basis function, a robust fault diagnosis and fault tolerant control algorithm is presented. A full-order observer is constructed to estimate the exogenous disturbance and an adaptive observer is used to estimate the fault size. A fault tolerant tracking controller is designed using the feedback of distribution tracking error, fault, and disturbance estimation to let the postfault output PDF still track desired distribution. Finally, a simulation example is included to illustrate the effectiveness of the proposed algorithms and encouraging results have been obtained.

1. Introduction In order to improve the reliability of practical stochastic systems, fault diagnosis and fault tolerant control for stochastic dynamic systems has long been one of the important areas of control theory and application [1–4]. Stochastic distribution control (SDC) system is a new branch of stochastic system control in which the output is the non-Gaussian probability density function (PDF) of the system output. The equations of these systems describe the relationship between the input PDF and output PDF of systems rather than the traditional relationship between input and output. SDC theory was proposed by professor Wang [5], which has been applied on some actual processes, for instance, the paper evenness control in the process of paper making, high polymer polymerization process of chemical industry, and flame distribution control. In the framework of nonGaussian SDC systems, it has very important theoretical significances and deep application prospect to develop fault diagnosis and fault tolerant control technology for complex industrial processes in which the control product quality and the distribution of indirect indicators need to be controlled.

With the study of fault diagnosis for non-Gaussian SDC systems, some fault diagnosis algorithms have been proposed. For fault diagnosis of SDC systems, observer or filter-based methods are mainly used so far [6–10], in which the information of system output PDF and other measured information generate residuals in order to analyze and estimate the change of fault. A stable filter-based residual generator is constructed such that the fault can be detected and diagnosed for general stochastic systems in [6]. In [7], a nonlinear neural network observer is designed for fault diagnosis in which the adaptive tuning rule for network parameters is determined by the Lyapunov stability theorem. In [8] a fault diagnosis algorithm is proposed based on iterative learning observer for SDC system. Otherwise, in [11], a novel fault-estimation observer is designed for Takagi-Sugeno (T-S) fuzzy systems with actuator faults, and the problem of fault tolerant control is addressed. For fault tolerant control, the active fault tolerant control is mainly used so far. Combined with the controller design method without fault, such as optimal control, PI control, sliding mode control, and model reference adaptive control, the controller can be reconfigured or reconstructed when

2 fault occurred. When the target PDF is known, the purpose of fault tolerant control is to make the output PDF of the systems still track a given PDF as close as possible after the fault happened. It has been shown in [8] that an optimal control strategy is designed to reconfigure the controller to compensate the influence of the fault on the system performance. A fault tolerant controller based on PI tracking control is proposed for a time-delayed SDC system in [12] and for singular time-delayed SDC system in [13]. In [14], a fuzzy fault tolerant control scheme is developed to guarantee the closed-loop system to be exponentially stable in mean square. When the target PDF is unknown, introduce the concept of entropy to fault tolerant control of non-Gaussian SDC systems, in which the purpose is to minimize the uncertainty of the system output after fault happened [15, 16]. Although the model descriptions of SDC systems are always the dynamic link relations between the input and the weights of output PDFs, the conventional dynamic links do not meet the practical requirements because of the existence of some algebraic constraint conditions between some state variables. Thus, it is necessary to study the SDC systems in the framework of singular systems [8, 13, 16–18]. Due to the material delivering by conveyor belt, system modeling, and data operation and transmission, time-delay widely exists in actual industrial systems. Time-delay can make system performance degradation and even make systems unstable; at the same time, time-delay will largely reduce the effectiveness of the fault diagnosis and fault tolerant control. The results of fault diagnosis and fault tolerant control for singular timedelayed SDC systems are focused on the model approximated by the linear or square-root B-spline model [13, 18]. The rational square-root B-spline model can guarantee that the weights of feedback control are positive at the same time independent of each other compared with other B-spline models [19]. External disturbance widely exists in various industrial processes and this situation will turn to be very complicated for fault diagnosis [9, 20], but controller design to eliminate the influence of disturbance is not considered. Thus, it is significant to study fault diagnosis and fault tolerant control for non-Gaussian singular time-delayed SDC system with external disturbance based on the rational square-root model approximation. In this paper, a robust fault diagnosis and fault tolerant control approach is proposed for non-Gaussian singular time-delayed stochastic distribution system with external disturbance based on the rational square-root approximation. The external disturbance is considered and supposed to be generated by a linear exogenous system, and a full-order observer is designed to estimate the disturbance. Then an adaptive observer is constructed to estimate fault information. The gain matrices can be determined by solving the corresponding linear matrix inequalities (LMIs). In order to eliminate the influence of fault and disturbance, a fault tolerant tracking controller is designed to track the desired PDF. An augmentation control input is defined to contain the feedback of output PDF tracking error and the estimation of fault and disturbance, so that the control input of the SDC system can eliminate the influence of fault and disturbance on the system performance, leading to fault tolerant tracking

Journal of Control Science and Engineering control. Finally, the computer simulation results show the effectiveness of the algorithm. The main contributions of this paper can be summarized as follows. (1) Comparing with most of the existing studies of SDC system with disturbance, integrated fault diagnosis and fault tolerant control is achieved; rather only fault diagnosis is considered. Besides, the influence of disturbance to the system performance can be rejected by the fault tolerant controller. (2) The rational square-root model is used to approximate the output PDF of the singular time-delayed SDC system, as a contrast, linear, or square-root model is used in the existing studies of singular time-delayed SDC system.

2. Model Description Denoting 𝛾(𝑦, 𝑢(𝑡)) as the PDF of the system output with 𝑦 being defined on a known bounded interval [𝑎, 𝑏], the continuous singular time-delayed SDC system can be expressed as follows: 𝐸𝑥̇ (𝑡) = 𝐴𝑥 (𝑡) + 𝐴 𝑑 𝑥 (𝑡 − 𝜏) + 𝐵𝑢 (𝑡) + 𝑁𝐹 (𝑡) + 𝐵𝑑 𝑑 (𝑡) ,

(1)

𝑉 (𝑡) = 𝐷𝑥 (𝑡) , 𝑥 (𝑡) = 𝜑 (𝑡) , 𝑡 ∈ [−𝛿, 0] , √𝛾 (𝑦, 𝑢 (𝑡)) =

𝐶 (𝑦) 𝑉 (𝑡) √𝑉𝑇 (𝑡) Σ1 𝑉 (𝑡)

, (2)

𝑏

Σ1 = ∫ 𝐶𝑇 (𝑦) 𝐶 (𝑦) 𝑑𝑦, 𝑎

where 𝑥(𝑡) ∈ 𝑅𝑛 is the state vector, 𝑢(𝑡) ∈ 𝑅𝑚 is the control input vector, 𝑉(𝑡) ∈ 𝑅𝑛 is the weight vector and 𝐹 ∈ 𝑅𝑚 is the fault vector, and 𝜏 is the time-delay term. 𝜑(𝑡) is a real valued continuous function. 𝐴 ∈ 𝑅𝑛×𝑛 , 𝐴 𝑑 ∈ 𝑅𝑛×𝑛 , 𝐵 ∈ 𝑅𝑛×𝑚 , 𝐷 ∈ 𝑅𝑛×𝑛 , 𝐵𝑑 ∈ 𝑅𝑛×𝑚 , 𝐸 ∈ 𝑅𝑛×𝑛 , and 𝑁 ∈ 𝑅𝑛×𝑚 are system parameter matrices with rank(𝐸) = 𝑞 < 𝑛. Equation (2) represents the static model of the output PDF approximated by the square-root B-spline model. It is denoted that 𝐶 (𝑦) = [𝜙1 (𝑦) , 𝜙2 (𝑦) , . . . , 𝜙𝑛 (𝑦)] , 𝑉 = [𝜔1 , 𝜔2 , . . . , 𝜔𝑛 ]

𝑇

(𝑉 ≠ 0) ,

(3)

where 𝜙𝑖 (𝑦) (𝑖 = 1, 2, . . . , 𝑛, 𝑛 ≥ 2) is the prespecified basis function, 𝜔𝑖 (𝑖 = 1, 2, . . . , 𝑛) is the approximation weight which is only related to 𝑢(𝑡), and 𝑛 is the number of basis functions. 𝑑(𝑡) is the unknown external disturbance which can be supposed to be generated by a linear exogenous system described by 𝜔̇ (𝑡) = 𝑊𝜔 (𝑡) , 𝑑 (𝑡) = 𝑇𝜔 (𝑡) .

(4)

Remark 1. From literatures [21], many kinds of disturbances in engineering can be described by this model, for example, unknown constant and harmonics with unknown phase and

Journal of Control Science and Engineering

3

magnitude. In most existing results, the disturbances are restricted to be bounded exogenous signals. The state variable of the exogenous system 𝜔(𝑡) is the target of estimation, 𝑑(𝑡) is the form of disturbance, and 𝑊 and 𝑇 are the known parameter matrices with suitable dimensions. The following assumptions are used throughout this paper. Assumption 2. (𝐴, 𝐷) is observable, and (𝑊, 𝐵𝑑 𝑇) is observable. Assumption 3. The fault and disturbance occurring in SDC system are bounded; that is, ‖𝐹‖ ≤ 𝑀𝑓 and ‖𝑑(𝑡)‖ ≤ 𝑀𝑑 , where 𝑀𝑓 and 𝑀𝑑 are two positive constants.

3. Fault Diagnosis In order to reject the disturbance and eliminate the influence to the fault diagnosis, a full-order observer is constructed to 1 (𝑡) estimate the external disturbance. Denoting 𝑠(𝑡) = [ 𝑥𝜔(𝑡) ], the full-order system is constructed by augmenting state (1) with exogenous system (4) as follows: 𝑠̇ (𝑡) = 𝐴 0 𝑠 (𝑡) + 𝐴 𝑑0 𝑠 (𝑡 − 𝜏) + 𝐵0 𝑢 (𝑡) + 𝑁0 𝐹 (𝑡) , 𝑥2 (𝑡) = −𝐴 𝑑2 𝑥2 (𝑡 − 𝜏) − 𝐵2 𝑢 (𝑡) − 𝑁2 𝐹 (𝑡) − 𝐵𝑑2 𝑇𝜔 (𝑡) , 𝑉 (𝑡) = 𝐷1 𝑥1 (𝑡) + 𝐷2 𝑥2 (𝑡) = 𝐷0 𝑠 (𝑡) + 𝐷2 𝑥2 (𝑡) , where

Assumption 4. The system (1) is regular and impulse-free; that is, |𝑠𝐸 − 𝐴| ≠ 0 and rank(𝐸) = deg(|𝑠𝐸 − 𝐴|).

𝐴0 = [

With the assumptions, there exist two matrices (𝐿 1 , 𝐿 2 ) such that the following equation

𝐴 𝑑0 = [

𝐿 1 𝐸𝐿 2 = [

𝐼𝑞 0 ], 0 0

𝐴1 0 𝐿 1 𝐴𝐿 2 = [ ] 0 𝐼𝑛−𝑞

𝑁0 = [

̂𝑠̇ (𝑡) = 𝐴 0̂𝑠 (𝑡) + 𝐴 𝑑0̂𝑠 (𝑡 − 𝜏) + 𝐵0 𝑢 (𝑡)

− 𝐵𝑑2 𝑇̂ 𝜔 (𝑡) , ̂ (𝑡) = 𝐷0̂𝑠 (𝑡) + 𝐷2 𝑥 ̂ 2 (𝑡) , 𝑉 √𝛾̂ (𝑦, 𝑢 (𝑡)) =

where 𝐵1 , 𝐵𝑑1 , 𝑁1 ∈ 𝑅𝑞×𝑚 , 𝐵2 , 𝐵𝑑2 , 𝑁2 ∈ 𝑅(𝑛−𝑞)×𝑚 , 𝐷1 ∈ 𝑅(𝑛−1)×𝑞 , and 𝐷2 ∈ 𝑅(𝑛−1)×(𝑛−𝑞) , which can be determined as follows: 𝐵1

𝐿 2𝐵 = [ ] , 𝐵2 𝐵𝑑2

],

The full-order observer is designed as follows:

(6)

𝜔̇ (𝑡) = 𝑊𝜔 (𝑡) ,

𝐿 2 𝐵𝑑 = [

0

̂ 2 (𝑡 − 𝜏) − 𝐵2 𝑢 (𝑡) − 𝑁2 𝐹 (𝑡) ̂ 2 (𝑡) = −𝐴 𝑑2 𝑥 𝑥

𝑉 (𝑡) = 𝐷1 𝑥1 (𝑡) + 𝐷2 𝑥2 (𝑡) ,

𝐵𝑑1

𝑁1

(9)

+ 𝐾1 𝜀 (𝑡) ,

+ 𝐵𝑑1 𝑇𝜔 (𝑡) , − 𝐵𝑑2 𝑇𝜔 (𝑡) ,

𝐴 𝑑1 0 ], 0 0

𝐷0 = [𝐷1 0] .

𝑥̇ 1 (𝑡) = 𝐴 1 𝑥1 (𝑡) + 𝐴 𝑑1 𝑥1 (𝑡 − 𝜏) + 𝐵1 𝑢 (𝑡) + 𝑁1 𝐹 (𝑡) 𝑥2 (𝑡) = −𝐴 𝑑2 𝑥2 (𝑡 − 𝜏) − 𝐵2 𝑢 (𝑡) − 𝑁2 𝐹 (𝑡)

𝐴 1 𝐵𝑑1 𝑇 ], 0 𝑊

𝐵1 𝐵0 = [ ] , 0

(5)

holds. It is assumed that 𝐿 1 𝐴 𝑑 𝐿 2 = [ 𝐴0𝑑1 𝐴0𝑑2 ] holds simultaneously. Then the SDC system (1) can be transformed as

],

𝐷𝐿 1 = [𝐷1 𝐷2 ] , 𝑁1 𝐿 2𝑁 = [ ] . 𝑁2

(8)

𝑇 √̂ ̂ (𝑡) 𝑉 (𝑡) Σ1 𝑉

,

𝑏

𝜀 (𝑡) = ∫ (√𝛾 − √𝛾̂) 𝑑𝑦, 𝑎

where 𝐾1 is observer gain matrix, 𝜀(𝑡) is the residual signal, and 𝜀 (𝑡) =

(7)

̂ (𝑡) 𝐶 (𝑦) 𝑉

(10)

Σ2 (𝐷0 𝑠 (𝑡) + 𝐷2 𝑥2 (𝑡)) √𝑉𝑇 (𝑡) Σ1 𝑉 (𝑡) −

+

̂ 2 (𝑡)) Σ2 (𝐷0̂𝑠 (𝑡) + 𝐷2 𝑥 𝑇 √̂ ̂ (𝑡) 𝑉 (𝑡) Σ1 𝑉

Σ2 (𝐷0 𝑠 (𝑡) + 𝐷2 𝑥2 (𝑡)) 𝑇 √̂ ̂ (𝑡) 𝑉 (𝑡) Σ1 𝑉

4

Journal of Control Science and Engineering −

=

Σ2 (𝐷0 𝑠 (𝑡) + 𝐷2 𝑥2 (𝑡))

The purpose of fault diagnosis is to estimate the size of fault. The fault diagnosis observer is constructed as follows:

𝑇

√̂ ̂ (𝑡) 𝑉 (𝑡) Σ1 𝑉

̂ 1 (𝑡) + 𝐴 𝑑1 𝑥 ̂ 1 (𝑡 − 𝜏) + 𝐵1 𝑢 (𝑡) ̂̇ 1 (𝑡) = 𝐴 1 𝑥 𝑥

Σ2 (𝐷0 𝑒𝑠 (𝑡) + 𝐷2 𝑒2 (𝑡)) 𝑇 √̂ ̂ (𝑡) 𝑉 (𝑡) Σ1 𝑉

+

̂ (𝑡) + 𝐵𝑑1 𝑇̂ + 𝑁1 𝐹 𝜔 (𝑡) + 𝐾2 𝜀 (𝑡) ,

Σ2 (𝐷0 𝑠 (𝑡) + 𝐷2 𝑥2 (𝑡))

̂ (𝑡) ̂ 2 (𝑡 − 𝜏) − 𝐵2 𝑢 (𝑡) − 𝑁2 𝐹 ̂ 2 (𝑡) = −𝐴 𝑑2 𝑥 𝑥

𝑇

√̂ ̂ (𝑡) 𝑉 (𝑡) Σ1 𝑉

⋅(

− 𝐵𝑑2 𝑇̂ 𝜔 (𝑡) ,

𝑇 √̂ ̂ (𝑡) − √𝑉𝑇 (𝑡) Σ1 𝑉 (𝑡) 𝑉 (𝑡) Σ1 𝑉

√𝑉𝑇 (𝑡) Σ1 𝑉 (𝑡)

),

√𝛾̂ (𝑦, 𝑢 (𝑡)) = (11)

Lemma 5 (see [19]). There exists a constant 𝜆 such that the following equation

󵄩̂ 󵄩󵄩 󵄩󵄩 𝑇 󵄩󵄩 = 𝜆 (󵄩󵄩󵄩󵄩𝑉 (𝑡)󵄩󵄩󵄩 − 󵄩󵄩󵄩𝑉 (𝑡)󵄩󵄩󵄩)

̂ (𝑡) 𝐶 (𝑦) 𝑉 𝑇 √̂ ̂ (𝑡) 𝑉 (𝑡) Σ1 𝑉

̂ 1 (𝑡). Then the observation error where 𝑒1 (𝑡) = 𝑥1 (𝑡) − 𝑥 dynamic system is obtained as follows: ̃ (𝑡) 𝑒1̇ (𝑡) = 𝐴 1 𝑒1 (𝑡) + 𝐴 𝑑1 𝑒1 (𝑡 − 𝜏) + 𝑁1 𝐹 𝜔 (𝑡) − 𝐿 4 Σ2 𝐷1 𝑒1 − 𝐿 4 Σ2 𝐷2 𝑒2 + 𝐵𝑑1 𝑇̃

(12)

+ 𝐿 4 Σ2 𝑉 (

holds, where 𝜆 min (Σ1 )/𝜆 max (Σ1 ) ≤ 𝜆 ≤ 𝜆 max (Σ1 )/𝜆 min (Σ1 ) and 𝜆 max (Σ1 ) and 𝜆 min (Σ1 ) are the maximum and minimum eigenvalues of matrix Σ1 , respectively.

󵄩̂󵄩󵄩 󵄩󵄩) 𝜆 2 (‖𝑉‖ − 󵄩󵄩󵄩󵄩𝑉 󵄩 √ 𝑉 𝑇 Σ1 𝑉

− 𝐿 5 Σ2 𝑉 (

𝜀 (𝑡)

󵄩̂󵄩󵄩 󵄩󵄩) 𝜆 2 (‖𝑉‖ − 󵄩󵄩󵄩󵄩𝑉 󵄩 √ 𝑉 𝑇 Σ1 𝑉

Σ2 (𝐷0 𝑒𝑠 (𝑡) + 𝐷2 𝑒2 (𝑡)) 𝑇 √̂ ̂ (𝑡) 𝑉 (𝑡) Σ1 𝑉

(13)

󵄩󵄩̂ 󵄩󵄩 󵄩󵄩 𝑇 󵄩󵄩 (𝑡)󵄩󵄩󵄩 − 󵄩󵄩󵄩𝑉 (𝑡)󵄩󵄩󵄩) Σ2 (𝐷0 𝑠 (𝑡) + 𝐷2 𝑥2 (𝑡)) 𝜆 1 (󵄩󵄩󵄩𝑉 + . 𝑇 √̂ √𝑉𝑇 (𝑡) Σ1 𝑉 (𝑡) ̂ (𝑡) 𝑉 (𝑡) Σ1 𝑉

𝑒𝑠̇ (𝑡) = 𝑠̇ (𝑡) − ̂𝑠̇ (𝑡) = 𝐴 0 𝑒𝑠 (𝑡) + 𝐴 𝑑0 𝑒𝑠 (𝑡 − 𝜏) + 𝑁0 𝐹

+ 𝐿 3 Σ2 𝑉 (

󵄩̂󵄩󵄩 󵄩󵄩) 𝜆 1 (‖𝑉‖ − 󵄩󵄩󵄩󵄩𝑉 󵄩 √ 𝑉 𝑇 Σ1 𝑉

),

),

(17)

𝑇

̃ = 𝐹 − 𝐹, ̂ 𝜔 ̂ ̃ = 𝜔−𝜔 ̂ . 𝐿 4 = 𝐾2 /√̂ where 𝐹 𝑉 (𝑡)Σ1 𝑉(𝑡), 𝐿5 = 𝑇

The full-order observation error dynamic system can be formulated as follows:

− 𝐿 3 Σ2 𝐷0 𝑒𝑠 (𝑡) − 𝐿 3 Σ2 𝐷2 𝑒2 (𝑡)

(16)

̂ + 𝐿 5 Σ2 𝐷1 𝑒1 + 𝐿 5 Σ2 𝐷2 𝑒2 ̂̇ = −Γ1 𝐹 𝐹

Then, it can be obtained that

=

,

̂̇ = −Γ1 𝐹 ̂ + Γ2 𝜀 (𝑡) , 𝐹

𝑏

where 𝑒𝑠 (𝑡) = 𝑠(𝑡) − ̂𝑠(𝑡), Σ2 = ∫𝑎 𝐶(𝑦)𝑑𝑦, and 𝑒2 (𝑡) = 𝑥2 (𝑡) − ̂ 2 (𝑡). 𝑥

√𝑉 ̂𝑇 (𝑡) Σ1 𝑉 ̂ (𝑡) − √𝑉𝑇 (𝑡) Σ1 𝑉 (𝑡)

(15)

̂ (𝑡) = 𝐷1 𝑥 ̂ 1 (𝑡) + 𝐷2 𝑥 ̂ 2 (𝑡) , 𝑉

(14)

̂ 𝑉 (𝑡)Σ1 𝑉(𝑡), and Γ1 are gain matrices to be determined Γ2 /√̂ later. Combing (14), (16), and (17), the following augmentation error dynamic system can be obtained as 𝑒𝑠̇ (𝑡)

𝑒𝑠 (𝑡) 𝐴 𝑑0 𝑒𝑠 (𝑡 − 𝜏) ̇ 𝑒 (𝑡) ( 1 ) = 𝐴 (𝑒1 (𝑡)) + (𝐴 𝑑1 𝑒1 (𝑡 − 𝜏)) ̂ (𝑡) ̂̇ (𝑡) 0 𝐹 𝐹 𝑁0 𝐹 (𝑡)

−𝐿 3 Σ2 𝐷2 𝑒2

+ (𝑁1 𝐹 (𝑡)) + (−𝐿 4 Σ2 𝐷2 𝑒2 ) ),

𝑇 ̂ is the gain matrix to be deterwhere 𝐿 3 = 𝐾1 /√̂ 𝑉 (𝑡)Σ1 𝑉(𝑡) mined later.

0 𝐿 3 Σ2 ℎ1 + ( 𝐿 4 Σ2 ℎ2 ) , −𝐿 5 Σ2 ℎ2

𝐿 5 Σ2 𝐷2 𝑒2

(18)

Journal of Control Science and Engineering

5 For the constants 𝜇1 > 0 and 𝜇2 > 0, the generalized 𝐻∞ performance is denoted as follows:

where

[ 𝐴=[ [

𝐴 0 − 𝐿 3 Σ2 𝐷0 𝑇 0

[

0

󵄩2 󵄩 󵄩 󵄩2 𝐽∞ = 󵄩󵄩󵄩𝑠∞ 󵄩󵄩󵄩 − 𝜇12 ‖𝐹 (𝑡)‖2 − 𝜇22 󵄩󵄩󵄩𝑒2 (𝑡)󵄩󵄩󵄩

0

] 𝐴 1 − 𝐿 4 Σ2 𝐷1 −𝑁1 ] ], 𝐿5 −Γ1 ]

− 𝛿 (𝑄1 , 𝑄2 , 𝑄3 , 𝑄4 ) , 𝛿 (𝑄1 , 𝑄2 , 𝑄3 , 𝑄4 ) = 𝑒𝑠𝑇 (0) 𝑄1 𝑒𝑠 (0) + 𝑒1𝑇 (0) 𝑄2 𝑒1 (0)

𝑇 = [0 𝐵𝑑1 𝑇] , (19) ℎ𝑖 = 𝑉 (

󵄩̂󵄩󵄩 󵄩󵄩) 𝜆 𝑖 (‖𝑉‖ − 󵄩󵄩󵄩󵄩𝑉 󵄩 √ 𝑉 𝑇 Σ1 𝑉

+∫ ),

0

−𝜏

𝑖 = 1, 2.

(21)

𝜑1𝑇 𝑄3 𝜑1 𝑑𝛼

0

+ ∫ 𝜑2𝑇 𝑄4 𝜑2 𝑑𝛼. −𝜏

Theorem 6. For the matrices 𝐶𝑖 > 0 (𝑖 = 1, 2, 3, 4, 5) and 𝑄𝑖 > 0 (𝑖 = 1, 2, 3, 4) and constants 𝜇𝑖 > 0 (𝑖 = 1, 2), there exist positive definite symmetric matrices 𝑃1 and 𝑃2 , positive definite matrices 𝑅1 and 𝑅2 , and gain matrices 𝐿 3 , 𝐿 4 , 𝐿 5 , and Γ1 satisfying 𝑃1 ≤ 𝑄1 , 𝑃2 ≤ 𝑄2 , 𝑅1 ≤ 𝑄3 , 𝑅2 ≤ 𝑄4 , and the following LMI, and then the augmenting error dynamic system (18) is stable and satisfies 𝐽∞ < 0

Select the reference output as follows: ̃ 𝑠∞ = 𝐶1 𝑒𝑠 (𝑡) + 𝐶2 𝑒1 (𝑡) + 𝐶3 𝑒1 (𝑡 − 𝜏) + 𝐶4 𝐹 + 𝐶5 𝑒2 (𝑡) . 𝑇

𝜉11 𝑃1 𝐴 𝑑0 𝑇 𝑃2𝑇

[ [∗ [ [ [∗ [ [ [∗ [ [ Φ1 = [ ∗ [ [ [∗ [ [ [∗ [ [∗ [ [∗

(20)

0

𝑃1 𝐿 3 Σ2

0

𝑃1 𝑁0

0

−𝑃1 𝐿 3 Σ2 𝐷2

0

0

0

0

−𝑅1

0

0

0



𝜉33

𝑃2 𝐴 𝑑1

0





−𝑅2

0

0

0







−𝐼

0

0

0

𝑃2 𝐿 4 Σ2 𝑃2 𝑁1

𝜉38 0









−𝐼

0

Σ𝑇2 𝐿 5











−𝜇12

0













−2Γ1















] ] ] ] −𝑃2 𝐿 4 Σ2 𝐷2 ] ] ] ] 0 ] ] 0 ] ] ] 0 ] ] ] 0 ] ] 𝐿 5 Σ2 𝐷2 ] ] −𝜇22

(22)

]

𝑇

+ [𝐶1 0 𝐶2 𝐶3 0 0 𝐶4 −𝐶4 𝐶5 ] [𝐶1 0 𝐶2 𝐶3 0 0 𝐶4 −𝐶4 𝐶5 ] < 0, 𝑡

where

𝜋2 = 𝑒1𝑇 𝑃2 𝑒1 + ∫ [𝑚2 𝑒1𝑇 (𝛼) 𝑒1 (𝛼) − ℎ2𝑇 ℎ2 ] 𝑑𝛼 0

𝑇

𝜉11 = (𝐴 0 − 𝐿 3 Σ2 𝐷0 ) 𝑃1 + 𝑃1 (𝐴 0 − 𝐿 3 Σ2 𝐷0 ) + 𝑅1

𝑡

+∫

+ 𝑚1 𝐼,

𝑡−𝜏

𝑇

𝜉33 = (𝐴 1 − 𝐿 4 Σ2 𝐷1 ) 𝑃2 + 𝑃2 (𝐴 1 − 𝐿 4 Σ2 𝐷1 ) + 𝑅2

(23)

𝑒1𝑇 (𝛼) 𝑅2 𝑒1 (𝛼) 𝑑𝛼,

̂ ̂ 𝑇 𝐹. 𝜋3 = 𝐹

+ 𝑚2 𝐼,

(24)

𝜉38 = 𝐷1𝑇 Σ𝑇2 𝐿𝑇5 − 𝑃2 𝑁1 . Proof. Select the following Lyapunov functions 𝜋1 , 𝜋2 , and 𝜋3 as follows: 𝑡

𝜋1 = 𝑒𝑠𝑇 𝑃1 𝑒𝑠 + ∫ [𝑚1 𝑒𝑠𝑇 (𝛼) 𝑒𝑠 (𝛼) − ℎ1𝑇 ℎ1 ] 𝑑𝛼 0

+∫

𝑡

𝑡−𝜏

𝑒𝑠𝑇 (𝛼) 𝑅1 𝑒𝑠

(𝛼) 𝑑𝛼,

It can be obtained that the first-order derivatives of 𝜋1 , 𝜋2 , and 𝜋3 are given as follows: 𝑇

𝜋̇ 1 = 𝑒𝑠𝑇 [(𝐴 0 − 𝐿 3 Σ2 𝐷0 ) 𝑃1 + 𝑃1 (𝐴 0 − 𝐿 3 Σ2 𝐷0 ) + 𝑅1 + 𝑚1 𝐼] 𝑒𝑠 + 2𝑒𝑠𝑇 𝑃1 𝐴 𝑑0 𝑒𝑠 (𝑡 − 𝜏) + 2𝑒𝑠𝑇 𝑃1 𝑁0 𝐹 − 2𝑒𝑠𝑇 𝑃1 𝐿 3 Σ2 𝐷2 𝑒2 + 2𝑒𝑠𝑇 𝑃1 𝐿 3 Σ2 ℎ1 − ℎ1𝑇 ℎ1 − 𝑒𝑠𝑇 (𝑡 − 𝜏) 𝑅1 𝑒𝑠 (𝑡 − 𝜏) ,

6

Journal of Control Science and Engineering 𝑇

𝜋̇ 2 = 𝑒1𝑇 [(𝐴 1 − 𝐿 4 Σ2 𝐷1 ) 𝑃2 + 𝑃2 (𝐴 1 − 𝐿 4 Σ2 𝐷1 )

Denote 𝜋 = 𝜋1 + 𝜋2 + 𝜋3 , and consider an auxiliary function as the performance index

+ 𝑅2 + 𝑚2 𝐼] 𝑒1 + 2𝑒1𝑇 𝑃2 𝑇𝑒𝑠 + 2𝑒1𝑇 𝑃2 𝑁1 𝐹 + 2𝑒1𝑇 𝑃2 𝐴 𝑑1 𝑒1 (𝑡 − 𝜏) − 2𝑒1𝑇 𝑃2 𝐿 4 Σ2 𝐷2 𝑒2

𝑡

󵄩2 󵄩 󵄩 󵄩2 𝐽aux = ∫ (󵄩󵄩󵄩𝑠∞ 󵄩󵄩󵄩 − 𝜇12 ‖𝐹 (𝛼)‖2 − 𝜇22 󵄩󵄩󵄩𝑒2 (𝛼)󵄩󵄩󵄩 0

̂ + 2𝑒𝑇 𝑃2 𝐿 4 Σ2 ℎ2 − ℎ𝑇 ℎ2 − 𝑒𝑇 (𝑡 − 𝜏) − 2𝑒1𝑇 𝑃2 𝑁1 𝐹 1 2 1

+ 𝜋̇ (𝛼)) 𝑑𝛼.

⋅ 𝑅2 𝑒1 (𝑡 − 𝜏) , ̂ 𝑇𝐹 ̂ + 2𝑒𝑇 𝐷𝑇 Σ𝑇 𝐿𝑇 𝐹 ̂ ̂𝑇 𝜋̇ 3 = −2Γ1 𝐹 1 1 2 5 + 2𝐹 𝐿 5 Σ2 𝐷2 𝑒2 ̂ 𝑇 𝐿 5 Σ2 ℎ2 . + 2𝐹 (25)

Then it can be obtained that ‖𝑠∞ ‖2 − 𝜇12 ‖𝐹(𝑡)‖2 − 𝜇22 ‖𝑒2 (𝑡)‖2 + 𝜋̇ = 𝑞𝑇 Φ1 𝑞, where 𝑇

̂ 𝑇 (𝑡) 𝑒𝑇 (𝑡)] . 𝑞𝑇 = [𝑒𝑠𝑇 (𝑡) 𝑒𝑠𝑇 (𝑡 − 𝜏) 𝑒1𝑇 (𝑡) 𝑒1𝑇 (𝑡 − 𝜏) ℎ1𝑇 (𝑡) ℎ2𝑇 (𝑡) 𝐹𝑇 (𝑡) 𝐹 2

From the Schur complement lemma, Φ1 equals inequality (22), leading to 𝐽aux < 0. Then 𝑡

√𝛾𝑔 (𝑦) =

󵄩2 󵄩 󵄩 󵄩2 + 𝜋 (𝑡) − 𝜋 (0) ≤ 󵄩󵄩󵄩𝑠∞ 󵄩󵄩󵄩 − 𝜇12 ‖𝐹 (𝑡)‖2 − 𝜇22 󵄩󵄩󵄩𝑒2 (𝑡)󵄩󵄩󵄩

0

0

−𝜏

−𝜏

(28)

−𝜏



[𝑒𝑠𝑇 (0) 𝑄1 𝑒𝑠

+∫

(0) +

−𝜏

𝑒1𝑇 (0) 𝑄2 𝑒1

0

0

−𝜏

−𝜏

(31)

Then the tracking error dynamic system can be expressed as follows:

+ 𝐵𝑑 𝑑 (𝑡) + (𝐴 + 𝐴 𝑑 ) 𝑥𝑔 .

= [𝑒𝑠𝑇 (0) 𝑃1 𝑒𝑠 (0) + 𝑒1𝑇 (0) 𝑃2 𝑒1 (0) +∫

(30)

= 𝐸 (𝑥̇ (𝑡) − 𝑥̇ 𝑔 ) .

̇ (𝑡) = 𝐴𝑒𝑚 (𝑡) + 𝐴 𝑑 𝑒𝑚 (𝑡 − 𝜏) + 𝐵𝑢 (𝑡) + 𝑁𝐹 (𝑡) 𝐸𝑒𝑚

𝛿 (𝑃1 − 𝑄1 , 𝑃2 − 𝑄2 , 𝑅1 − 𝑄3 , 𝑅2 − 𝑄4 )

𝜙2𝑇 𝑅2 𝜙2 𝑑𝛼]

, ∀𝑦 ∈ [𝑎, 𝑏] ,

̇ (𝑡) 𝐸𝐷−1 𝑒2̇ (𝑡) = 𝐸𝐷−1 (𝑉̇ (𝑡) − 𝑉̇ 𝑔 ) = 𝐸𝑒𝑚

𝐽∞ < 𝐽aux + 𝛿 (𝑃1 − 𝑄1 , 𝑃2 − 𝑄2 , 𝑅1 − 𝑄3 , 𝑅2 − 𝑄4 ) ,

0

√𝑉𝑔𝑇 Σ1 𝑉𝑔

𝑒2 (𝑡) = 𝑉 (𝑡) − 𝑉𝑔 ,

̂ 1 = 0, 𝑠 = 𝜑1 (𝑡), and 𝑥1 = where when 𝑡 ∈ [−𝜏, 0], ̂𝑠 = 0, 𝑥 𝜑2 (𝑡). It can be obtained that

𝜙1𝑇 𝑅1 𝜙1 𝑑𝛼

𝐶 (𝑦) 𝑉𝑔

where 𝑉𝑔 is the expected weight vector of desired PDF 𝛾𝑔 (𝑦). Let

+ ∫ 𝜑1𝑇 𝑅1 𝜑1 𝑑𝛼 + ∫ 𝜑2𝑇 𝑅2 𝜑2 𝑑𝛼] ,

0

(27)

PDF still track the desired PDF. A given desired PDF can be described as follows:

󵄩2 󵄩 󵄩 󵄩2 𝐽aux = ∫ (󵄩󵄩󵄩𝑠∞ 󵄩󵄩󵄩 − 𝜇12 ‖𝐹 (𝑡)‖2 − 𝜇22 󵄩󵄩󵄩𝑒2 (𝑡)󵄩󵄩󵄩 ) 𝑑𝜏 0

− [𝑒𝑠𝑇 (0) 𝑃1 𝑒𝑠 (0) + 𝑒1𝑇 (0) 𝑃2 𝑒1 (0)

(26)

(32)

Denote (29)

(0)

+ ∫ 𝜙1𝑇 𝑄3 𝜙1 𝑑𝛼 + ∫ 𝜙2𝑇 𝑄4 𝜙2 𝑑𝛼] . Since 𝛿(𝑃1 − 𝑄1 , 𝑃2 − 𝑄2 , 𝑅1 − 𝑄3 , 𝑅2 − 𝑄4 ) < 0, 𝐽∞ < 0 is proved. This completes the proof.

4. Fault Tolerant Control After estimating the fault and disturbance, it is necessary to design a fault tolerant controller to make the postfault output

𝐸 = 𝐿 1 𝐸𝐿 2 = [ 𝐿−1 2 𝑒𝑚 (𝑡) = [

𝜉1 (𝑡) 𝜉2 (𝑡)

𝐼𝑞 0 0 0

], (33)

] = 𝜉 (𝑡) .

Then an augmentation control input containing the system input and the information of desired PDF is constructed as 𝑈(𝑡) = 𝑢(𝑡) + (𝐵𝑇 𝐵)−1 𝐵𝑇 (𝐴 + 𝐴 𝑑 )𝑥𝑔 , where the former item is the reconstructed control input and the latter is the extended information of desired PDF. The tracking error dynamic system can be rewritten as follows: 𝐸𝜉̇ (𝑡) = 𝐿 1 𝐴𝐿 2 𝜉 (𝑡) + 𝐿 1 𝐴 𝑑 𝐿 2 𝜉 (𝑡 − 𝜏) + 𝐿 1 𝐵𝑈 (𝑡) + 𝐿 1 𝑁𝐹 (𝑡) + 𝐿 1 𝐵𝑑 𝑑 (𝑡) .

(34)

Journal of Control Science and Engineering

7

Assume that 𝑈(𝑡) is the feedback of the tracking error, fault, and disturbance, then it can be obtained that

−1

𝜋̇ 4 = 𝜉𝑇 [𝑃 𝐿 1 𝐴𝐿 2 + 𝐿𝑇2 𝐴𝑇 𝐿𝑇1 𝑃

𝑏

𝑈 (𝑡) = Γ3 ∫ (√𝛾 (𝑦, 𝑢 (𝑡)) − √𝑔 (𝑦)) 𝑑𝑦 + Γ4 𝐹 (𝑡) 𝑎

Γ3 Σ2 𝐷𝐿 2 𝜉 (𝑡) √𝑉𝑔𝑇 Σ1 𝑉𝑔 ⋅

+

−1

−1

+ 2𝜉𝑇 𝑃 𝐿 1 𝐴 𝑑 𝐿 2 𝜉 (𝑡 − 𝜏) − ℎ3𝑇 ℎ3 + 2𝜉𝑇 (𝑃 𝐿 1 𝑁

√𝑉𝑔𝑇 Σ1 𝑉𝑔

√𝑉𝑇 (𝑡) Σ1 𝑉 (𝑡)

(35)

−1

−1

+ 𝑃 𝐿 1 𝐵Γ4 ) 𝐹 + 2𝜉𝑇 𝑃 𝐿 1 𝐵𝑑 𝑑 (𝑡) −1

+ 2𝜉𝑇 𝑃 𝐿 1 𝐵Γ5 𝑑 (𝑡) − 𝜉𝑇 (𝑡 − 𝜏) 𝑅1 𝜉 (𝑡 − 𝜏)

+ Γ4 𝐹 (𝑡) + Γ5 𝑑 (𝑡)

−1

−1

+ 2𝜉𝑇 𝑃 𝐿 1 𝐵𝐿 6 Σ2 ℎ3 ≤ 𝜉𝑇 [𝑃 𝐿 1 𝐴𝐿 2

= 𝐿 6 Σ2 𝐷𝐿 2 𝜉 (𝑡) + 𝐿 6 Σ2 ℎ3 + Γ4 𝐹 (𝑡) + Γ5 𝑑 (𝑡) ,

+ 𝐿𝑇2 𝐴𝑇 𝐿𝑇1 𝑃 =

Γ3 /√𝑉𝑔𝑇 Σ1 𝑉𝑔 and ℎ3

=

(𝜆 3 (‖𝑉𝑔 ‖ −

‖𝑉(𝑡)‖)/√𝑉𝑇 (𝑡)Σ1 𝑉(𝑡))𝑉(𝑡).

𝑇 𝑀4 𝐴 𝑑 𝑃𝑇 𝐵𝐿 6 Σ2 𝐿−1 2 𝑃

−𝑅1

0



−𝐼









−1

+ 𝑅1 + 2𝑃 𝐿 1 𝐵𝐿 6 Σ2 𝐷𝐿 2 + 𝛼3 ] 𝜉 −1

(38)

+ 𝜂42 𝑑𝑇 (𝑡) 𝑑 (𝑡) + 2𝜉𝑇 𝑃 𝐿 1 𝐴 𝑑 𝐿 2 𝜉 (𝑡 − 𝜏) − 𝜉𝑇 (𝑡 − 𝜏) 𝑅1 𝜉 (𝑡 − 𝜏) + 2𝜉𝑇 𝑃 𝐿 1 𝐵𝐿 6 Σ2 ℎ3 − ℎ3𝑇 ℎ3 + +

𝑇 −1 1 𝑇 −1 𝜉 𝑃 𝐿 1 𝑁 (𝑃 𝐿 1 𝑁) 𝜉 + 𝜂12 𝐹𝑇 𝐹 2 𝜂1

𝑇 −1 1 1 𝑇 −1 𝜉 𝑃 𝐿 1 𝐵Γ4 (𝑃 𝐿 1 𝐵Γ4 ) 𝜉 + 𝜂22 𝐹𝑇 𝐹 + 2 2 𝜂2 𝜂3 −1

1 1 𝐵Γ4 𝐵Γ 𝜂2 𝜂4 5 ] ] 0 0 ] ] ] 0 0 ] ] −𝐼 0 ] ] ∗ −𝐼 ]

𝑇

−1

⋅ 𝜉𝑇 𝑃 𝐿 1 𝐵𝑑 (𝑃 𝐿 1 𝐵𝑑 ) 𝜉 + 𝜂32 𝑑𝑇 (𝑡) 𝑑 (𝑡) + −1

1 𝜂42

𝑇

−1

⋅ 𝜉𝑇 𝑃 𝐿 1 𝐵Γ5 (𝑃 𝐿 1 𝐵Γ5 ) 𝜉 = 𝑞𝑇 Φ1 𝑞 + (𝜂12 + 𝜂22 ) 𝐹𝑇 𝐹 + (𝜂32 + 𝜂42 ) 𝑑𝑇 (𝑡) 𝑑 (𝑡) , (36)

< 0, −1 𝑇 𝑇 𝑀4 = 𝐴𝑃𝑇 + 𝑃𝐴𝑇 + 𝑃𝐿−𝑇 2 𝑅1 𝐿 2 𝑃 + 2𝐵𝐿 6 Σ2 𝐷𝑃 −1 𝑇 + 𝑃𝐿−𝑇 2 𝛼3 𝐿 2 𝑃 +

−𝑇

−1

Theorem 7. For the given constant 𝜂𝑖 (𝑖 = 1, . . . , 4) and a small positive constant 𝜆, suppose that there exist positive definite symmetric matrices 𝑃 and matrices Γ3 , Γ4 , and Γ5 such that the following LMI holds, and then the tracking error dynamic system (34) is stable.

[ [ [∗ [ Φ4 = [ [∗ [ [∗ [ [∗

+ 𝑅1

−1

Γ3 Σ2 𝑉 (𝑡)

󵄩 󵄩 𝜆 3 (󵄩󵄩󵄩󵄩𝑉𝑔 󵄩󵄩󵄩󵄩 − ‖𝑉 (𝑡)‖)

where 𝐿 6

−𝑇

+ 2𝑃 𝐿 1 𝐵𝐿 6 Σ2 𝐷𝐿 2 ] 𝜉 + 𝛼3 𝜉𝑇 𝜉

+ Γ5 𝑑 (𝑡) =

denote 𝛼3 = (𝜆 3 ‖𝐷𝐿 2 ‖/√‖Σ1 ‖)2 ; then 𝜋4 ≥ 0 holds. The firstorder derivative of 𝜋4 can be obtained as follows:

where 𝑞𝑇 = [𝜉𝑇 (𝑡) 𝜉𝑇 (𝑡 − 𝜏) ℎ3𝑇 ] and −1

−1

𝑀 𝑃 𝐿 1 𝐴 𝑑 𝐿 2 𝑃 𝐿 1 𝐵𝐿 6 Σ2 ] [ 1 ], Φ1 = [ −𝑅1 0 ] [∗ ∗ −𝐼 ] [∗

1 1 𝑁𝑁𝑇 + 2 𝐵𝑑 𝐵𝑑𝑇 + 𝜆𝐼, 𝜂12 𝜂3

−1

Γ3 = 𝐿 6 √𝑉𝑔𝑇 Σ1 𝑉𝑔 .

𝑀1 = 𝑃 𝐿 1 𝐴𝐿 2 + 𝐿𝑇2 𝐴𝑇 𝐿𝑇1 𝑃

−𝑇

+ 𝑅1 + 𝛼3

−1

Proof. Select the following Lyapunov function as follows: 𝑇 −𝑇

𝜋4 = 𝜉𝑇 𝐸 𝑃 𝜉 + ∫

𝑡

𝑡−𝜏

𝑡

𝜉𝑇 (𝑠) 𝑅1 𝜉 (𝑠) 𝑑𝑠 (37)

+ ∫ (𝛼3 𝜉𝑇 𝜉 − ℎ3𝑇 ℎ3 ) 𝑑𝑠, 0

where ℎ3 = (𝜆 3 (‖𝑉𝑔 ‖ − ‖𝑉(𝑡)‖)/√𝑉𝑇 (𝑡)Σ1 𝑉(𝑡))𝑉(𝑡). It can

be obtained that ℎ3𝑇 ℎ3 ≤ (𝜆 3 ‖𝐷𝐿 2 ‖/√‖Σ1 ‖)2 𝜉𝑇 (𝑡)𝜉(𝑡), and

+ 2𝑃 𝐿 1 𝐵𝐿 6 Σ2 𝐷𝐿 2 +

𝑇 −1 1 −1 𝑃 𝐿 1 𝑁 (𝑃 𝐿 1 𝑁) 2 𝜂1

+

𝑇 −1 1 −1 𝑃 𝐿 1 𝐵Γ4 (𝑃 𝐿 1 𝐵Γ4 ) 2 𝜂2

+

𝑇 −1 1 −1 𝑃 𝐿 1 𝐵𝑑 (𝑃 𝐿 1 𝐵𝑑 ) 2 𝜂3

+

𝑇 −1 1 −1 𝑃 𝐿 1 𝐵Γ5 (𝑃 𝐿 1 𝐵Γ5 ) . 2 𝜂4

(39)

8

Journal of Control Science and Engineering −1

= 𝐿𝑇2 𝑃−1 𝐿−1 Substituting 𝑃 = 𝐿 1 𝑃𝐿−𝑇 2 and 𝑃 1 into −𝑇 Φ1 and premultiplying diag(𝑃𝐿 2 𝐼 𝐼) and postmultiplying 𝑇 𝐼 𝐼) to the left and the right side of Φ1 , it can be diag(𝐿−1 2 𝑃 obtained that 𝑇 𝑀2 𝐴 𝑑 𝑃𝑇 𝐵𝐿 6 Σ2 𝐿−1 2 𝑃 ] [ ], Φ2 = [ 0 ] [ ∗ −𝑅1 ∗ ∗ −𝐼 ] [

+ +

2

2

𝜙2 (𝑦) = 0.5 (𝑦 − 3) 𝐼2 + (−𝑦2 + 9𝑦 − 19.5) 𝐼3 2

(40)

1 1 + 2 𝑁𝑁𝑇 + 2 𝐵Γ4 Γ4𝑇 𝐵𝑇 𝜂1 𝜂2

[ Φ3 = [ [∗

where 𝐼𝑖 (𝑖 = 1, 2, . . . , 5) is an interval function defined as

𝑇 𝐵𝐿 6 Σ2 𝐿−1 2 𝑃

−𝑅1

0

] ], ]



−𝐼

]

[∗

{1, 𝑦 ∈ [𝑖 + 1, 𝑖 + 2] 𝐼𝑖 = { 0, otherwise, {

𝜆

(42)

(43)

and the tracking error dynamic system (34) is stable. ̂ into (35), the practical fault ̂ and 𝑑(𝑡) By substituting 𝐹(𝑡) tolerant tracking controller can be formulated as follows:

[0 0 0] −3.2 0.9 −0.2 [ 0.3 −3.51 −0.1] 𝐴=[ ], 0

(44)

−1

̂𝑔. − (𝐵𝑇 𝐵) 𝐵𝑇 (𝐴 + 𝐴 𝑑 ) 𝑥

−1 ]

−0.1 −0.5 0.12 [ 0.2 −0.25 0.06] 𝐴𝑑 = [ ], [ 0

0

0.6 ]

0.1 [0.3] 𝐻 = [ ], [0.3] 1 [ 0 𝐵=[

−1

̂𝑔 𝑢 (𝑡) = 𝑈 (𝑡) − (𝐵𝑇 𝐵) 𝐵𝑇 (𝐴 + 𝐴 𝑑 ) 𝑥 ̂ (𝑡) ̂ (𝑡) + Γ5 𝑑 = 𝐿 6 Σ2 𝐷𝐿 2 𝜉 (𝑡) + 𝐿 6 Σ2 ℎ3 + Γ4 𝐹

1 0 0 [0 1 0] 𝐸=[ ],

[ 0

is satisfied, and then 𝜋̇ 4 < 0, leading to 2 2 2 2 2 2 󵄩󵄩 󵄩󵄩2 𝑀𝑓 (𝜂1 + 𝜂2 ) + 𝑀𝑑 (𝜂3 + 𝜂4 ) , 󵄩󵄩𝜉󵄩󵄩 ≤ 𝜆

(46)

The system parameter matrices of the SDC system (1) are given as follows:

When Φ3 < 0, it can be acquired that 𝜋̇ 4 ≤ −𝜆𝜉𝑇 𝜉 + (𝜂12 + 𝜂22 )𝐹𝑇 𝐹 + (𝜂32 + 𝜂42 )𝑑𝑇 (𝑡)𝑑(𝑡). According to the Schur complement lemma, Φ3 < 0 is equivalent to Φ4 < 0. Thus, when Φ4 < 0 holds, 𝑀𝑓2 (𝜂12 + 𝜂22 ) + 𝑀𝑑2 (𝜂32 + 𝜂42 )

𝑖 = 1, 2, . . . , 5.

(41)

𝑀3 = 𝑀2 + 𝜆𝐼.

󵄩󵄩 󵄩󵄩2 󵄩󵄩𝜉󵄩󵄩 >

2

𝜙3 (𝑦) = 0.5 (𝑦 − 4) 𝐼3 + (−𝑦2 + 11𝑦 − 29.5) 𝐼4 2

Adding 𝜆𝐼 to 𝑀2 , then 𝑀3 𝐴 𝑑 𝑃

(45)

+ 0.5 (𝑦 − 6) 𝐼4 ,

+ 0.5 (𝑦 − 7) 𝐼5 ,

1 1 𝐵 𝐵𝑇 + 𝐵Γ Γ𝑇 𝐵𝑇 . 𝜂32 𝑑 𝑑 𝜂42 5 5

𝑇

2

𝜙1 (𝑦) = 0.5 (𝑦 − 2) 𝐼1 + (−𝑦2 + 7𝑦 − 11.5) 𝐼2 + 0.5 (𝑦 − 5) 𝐼3 ,

−1 𝑇 𝑇 𝑀2 = 𝐴𝑃𝑇 + 𝑃𝐴𝑇 + 𝑃𝐿−𝑇 2 𝑅1 𝐿 2 𝑃 + 2𝐵𝐿 6 Σ2 𝐷𝑃 −1 𝑇 𝑃𝐿−𝑇 2 𝛼3 𝐿 2 𝑃

following B-spline functions 𝜙𝑖 (𝑦), 𝑖 = 1, 2, 3, is considered as

1 2

[0.002 0.005

−0.01

−0.02] ], 0.1 ]

1 0 0 [0 1 0] 𝐷=[ ], [0 0 1]

5. A Simulation Example To illustrate the effectiveness of the proposed algorithms, a SDC system whose output PDF can be approximated by the

0.1 [ 0 𝐵𝑑 = [

0.1 0.2

[0.002 0.005

−0.01

−0.02] ], 0.1 ]

Journal of Control Science and Engineering

9

1.0 0 −0.2 [ 0 1.0 −0.1] 𝐿1 = [ ], 0

0.1

1.0 ]

1 0 0 [0 1 0] 𝐿2 = [ ]. [0 0 1] (47) The exogenous disturbance can be described by (4), and the corresponding parameter matrices are given as follows: 2 0 0 [0 0 0 ] 𝑇=[ ],

Disturbance and its estimation

[0

0.15

0.05 0 −0.05 −0.1 −0.15 −0.2

0

2

4

6

8

10

12

14

16

18

20

t (s)

[0 0 0 ]

(48)

0 4 0 [−4 1 0] 𝑊=[ ].

Disturbance Disturbance estimation

Figure 1: Disturbance and its estimation.

[ 0 0 2] The time-delay 𝜏 is selected as 𝜏 = 2 s and the sampling time is chosen as 0.1 s. It is assumed that the fault is constructed as follows:

𝑅2 = [

(49)

𝑇

𝑇

𝐿 4 = [0.2435 −0.0223] , Γ1 = 0.01,

𝐶1 = [0.1 0.1 0.1 0.1 0.1] 𝐶2 = 𝐶3 = [0.2 0.2] ,

(50)

𝐶4 = 𝐶5 = 0.1. It can be calculated from the LMIs of Theorems 6 and 7 that

Γ2 = −0.7, 0.1201 0.0134 [0.0134 0.1206 𝑃=[ [

𝑃1 [ [ 0.1950 [ [ = [−0.5180 [ [ 0.5660 [ [ 0.0002 𝑃2 = [

0.002

] 0.003 ] ] −0.0773 6.8373 −0.0841 −0.0026] ], ] 0.1075 −0.0841 6.9108 0.0017 ] ] 0.0003 −0.0026 0.0017 0.0956 ] 2.0251 −0.0773 0.1075

1.160 0.0599 0.0599 2.1566

0

0

0 0

] ],

0.3652]

Γ3 = [−0.0169 −0.2712 −3.1352] , Γ4 = [0.5301 −0.1530 0.1035] , 0.2905 −0.0639 0.0098 [−0.0639 0.0669 0.0206] Γ5 = [ ]. [ 0.0098

0.0206 0.8149] (51)

],

𝑅1 7.0146 [ [−0.0797 [ [ = [−1.4354 [ [ 1.8075 [ [−0.0037

],

𝐿 5 = 0.1038,

The reference output is denoted such that

0.1950 −0.5180 0.5660

−0.5060 7.5992

𝐿 3 = [0.0090 −0.0316 0.0109 −0.0143 0] ,

𝑡 < 5s {0, 𝐹={ 0.5, 𝑡 > 5 s. {

1.8716

2.6380 −0.5060

−0.0797 −1.4354 1.8075 −0.0037

] 0.1264 −0.0035] ] 0.0394 1.3237 −0.3724 0.0004 ] ], ] 0.1264 −0.3724 1.4959 −0.0011] ] −0.0035 0.0004 −0.0011 1.1477 ] 7.5148

0.0394

The disturbance and its estimation have been shown in Figure 1. The fault diagnosis result has been presented in Figure 2. It can be seen that the desired fault diagnosis result has been obtained and the influence of disturbance has been rejected. The initial PDF, the desired PDF, and the final PDF with fault tolerant control are given in Figure 3. In Figures 4 and 5, the three-dimensional mesh plot shows the changes of the output PDF without and with fault tolerant control, respectively. Before fault occurs, the fault tolerant controller has not been reconstructed, so the disturbance cannot be

Journal of Control Science and Engineering

0.5

0.8

0.4

0.6 𝛾(y, u(t))

Fault and fault estimation

10

0.3

0.4 0.2

0.2

0 8

0.1

20

6 y

0

0

2

4

6

8

10 12 t (s)

14

16

18

4 2

20

0

5

10 t (s)

15

Figure 4: The output PDF of the whole process without fault tolerant control.

Fault estimation Fault

Figure 2: Fault and its estimation. 0.8

0.7

0.6 𝛾(y, u(t))

0.6 0.5

0.4 0.2

0.4 𝛾

0 8

0.3

6 y

0.2

4 2

0.1 0

0

5

10 t (s)

15

20

Figure 5: The output PDF of the whole process with fault tolerant control. 2

2.5

3

Initial PDF Final PDF

3.5

4

4.5 y

5

5.5

6

6.5

7

Expected PDF

Figure 3: The output PDF with fault tolerant control when fault occurs.

rejected. The influence of disturbance is shown in Figure 4. By comparing Figures 4 and 5, the effectiveness of the fault tolerant control can be seen.

make the postfault PDF still track the desired PDF. The Lyapunov stability theorem and LMI method are used to analyze the stability of the augmentation observation error dynamic system and tracking error dynamic system, and 𝐻∞ performance of fault diagnosis is guaranteed. The simulations have further confirmed the efficiency of the proposed fault diagnosis and fault tolerant control algorithm.

Competing Interests

6. Conclusions

The authors declare that they have no competing interests.

In this paper, a fault diagnosis and fault tolerant control algorithm is given for the non-Gaussian singular timedelayed SDC system based on the rational square-root Bspline approximation model. The external disturbance is taken into consideration. A full-order observer is designed to estimate the disturbance, and then an adaptive observer is constructed to estimate the fault. Using the feedback of output PDF tracking error and the estimation of fault and disturbance, the augmentation control input is designed to

Acknowledgments This work was supported by the National Natural Science Foundation of China (No. 61374128), State Key Laboratory of Synthetical Automation for Process Industries, the Science and Technology Innovation Talents 14HASTIT040 in Colleges and Universities in Henan Province, China, and Excellent Young Scientist Development Foundation 1421319086 of Zhengzhou University, China.

Journal of Control Science and Engineering

References [1] E. A. Garc´ıa and P. M. Frank, “Deterministic nonlinear observer-based approaches to fault diagnosis: a survey,” Control Engineering Practice, vol. 5, no. 5, pp. 663–670, 1997. [2] J. Chen and R. J. Pantton, Robust Model-Based Fault Diagnosis for Dynamic Systems, Kluwer Academic, Boston, Mass, USA, 1999. [3] R. Isermann, “Model-based fault-detection and diagnosisstatus and applications,” Annual Reviews in Control, vol. 29, no. 1, pp. 71–85, 2005. [4] E. Martin, J. Morris, and S. Lane, “Monitoring process manufacturing performance,” IEEE Control Systems Magazine, vol. 22, no. 5, pp. 26–39, 2002. [5] H. Wang, Bounded Dynamic Stochastic Systems, Advances in Industrial Control, Springer, London, UK, 2000. [6] L. Guo and H. Wang, “Fault detection and diagnosis for general stochastic systems using B-spline expansions and nonlinear filters,” IEEE Transactions on Circuits and Systems I: Regular Papers, vol. 52, no. 8, pp. 1644–1652, 2005. [7] L. Yao, J. Qin, A. Wang, and H. Wang, “Fault diagnosis and faulttolerant control for non-Gaussian non-linear stochastic systems using a rational square-root approximation model,” IET Control Theory & Applications, vol. 7, no. 1, pp. 116–124, 2013. [8] L. N. Yao, J. F. Qin, H. Wang, and B. Jiang, “Design of new fault diagnosis and fault tolerant control scheme for non-Gaussian singular stochastic distribution systems,” Automatica, vol. 48, no. 9, pp. 2305–2313, 2012. [9] L. Guo and Y. Zhang, “Generalized robust 𝐻∞ fault diagnosis filtering based on conditional stochastic distributions of system outputs,” Proceedings of the Institution of Mechanical Engineers. Part I: Journal of Systems & Control Engineering, vol. 221, no. 6, pp. 857–864, 2007. [10] L. Yao, Z. Yin, and H. Wang, “Robust 𝐻∞ fault diagnosis for stochastic distribution systems with disturbance rejection performance,” International Journal of Modelling, Identification and Control, vol. 21, no. 3, pp. 288–294, 2014. [11] Q. K. Shen, B. Jiang, and V. Cocquempot, “Fault-tolerant control for T–S fuzzy systems with application to near-space hypersonic vehicle with actuator faults,” IEEE Transactions on Fuzzy Systems, vol. 20, no. 4, pp. 652–665, 2012. [12] L. N. Yao and B. Peng, “Fault diagnosis and fault tolerant control for the non-Gaussian time-delayed stochastic distribution control system,” Journal of the Franklin Institute. Engineering and Applied Mathematics, vol. 351, no. 3, pp. 1577–1595, 2014. [13] L. Yao and L. Feng, “Fault diagnosis and fault tolerant control for non-Gaussian time-delayed singular stochastic distribution systems,” International Journal of Control, Automation and Systems, vol. 14, no. 2, pp. 235–442, 2016. [14] M. Liu, X. Cao, and P. Shi, “Fault estimation and tolerant control for fuzzy stochastic systems,” IEEE Transactions on Fuzzy Systems, vol. 21, no. 2, pp. 221–229, 2013. [15] H. Wang, “Minimum entropy control of non-Gaussian dynamic stochastic systems,” IEEE Transactions on Automatic Control, vol. 47, no. 2, pp. 398–403, 2002. [16] L. Yao, Y. Guan, and A. Wang, “Fault diagnosis and minimum entropy fault tolerant control for non-Gaussian singular stochastic distribution systems using square-root approximation,” International Journal of Modelling, Identification and Control, vol. 24, no. 3, pp. 206–215, 2015.

11 [17] J. Zhou, G. Li, and H. Wang, “Robust tracking controller design for non-Gaussian singular uncertainty stochastic distribution systems,” Automatica, vol. 50, no. 4, pp. 1296–1303, 2014. [18] L. N. Yao and L. Feng, “Fault diagnosis and fault tolerant tracking control for the non-Gaussian singular time-delayed stochastic distribution system with PDF approximation error,” Neurocomputing, vol. 175, pp. 538–543, 2016. [19] L.-N. Yao and H. Wang, “Fault diagnosis and tolerant control for non-Gaussian stochastic distribution control systems based on the rational square-root approximation model,” Control Theory and Applications, vol. 23, no. 4, pp. 561–568, 2006. [20] S. Cao, Y. Yi, and L. Guo, “Anti-disturbance fault diagnosis for non-Gaussian stochastic distribution systems with multiple disturbances,” Neurocomputing, vol. 136, pp. 315–320, 2014. [21] L. Guo and W.-H. Chen, “Disturbance attenuation and rejection for systems with nonlinearity via DOBC approach,” International Journal of Robust and Nonlinear Control, vol. 15, no. 3, pp. 109–125, 2005.

International Journal of

Rotating Machinery

Engineering Journal of

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

The Scientific World Journal Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

International Journal of

Distributed Sensor Networks

Journal of

Sensors Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Journal of

Control Science and Engineering

Advances in

Civil Engineering Hindawi Publishing Corporation http://www.hindawi.com

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014

Submit your manuscripts at http://www.hindawi.com Journal of

Journal of

Electrical and Computer Engineering

Robotics Hindawi Publishing Corporation http://www.hindawi.com

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014

VLSI Design Advances in OptoElectronics

International Journal of

Navigation and Observation Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Hindawi Publishing Corporation http://www.hindawi.com

Chemical Engineering Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014

Active and Passive Electronic Components

Antennas and Propagation Hindawi Publishing Corporation http://www.hindawi.com

Aerospace Engineering

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014

International Journal of

International Journal of

International Journal of

Modelling & Simulation in Engineering

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Shock and Vibration Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Advances in

Acoustics and Vibration Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014