# Fault Diagnosis and Fault Tolerant Control for Non-Gaussian Singular Time-Delayed Stochastic Distribution Systems

Authors: Lina Yao; Long Feng; Bin Jiang
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

New integrated fault diagnosis (FD) and fault tolerant control (FTC) algorithm for non-Gaussian singular time-delayed stochastic distribution control (SDC) system is proposed in the paper. Different from general SDC systems, in singular time-delayed SDC systems, the relationship between the weights and the control input is expressed by a singular time-delayed state space model, which largely increases the difficulty in the FD and FTC design. An iterative learning observer (ILO) is designed to carry out the fault estimation. The influence of the time delay term is eliminated in the process of the proof of the observation error stability. The fault may be constant, slow varying, or fast varying. Combined with the switching control theory, based on the estimated fault information, the fault tolerant controller can be designed to make the postfault probability density function (PDF) still track the given distribution. Simulations are given to show the effectiveness of the proposed integrated FD and FTC algorithms.


---

## Body

## 1. Introduction

Reliability and stability are of paramount importance for practical processes. Fault detection, diagnosis (FDD), and fault tolerant control (FTC) have long been regarded as an important and integrated part in control system design. The research of FDD and FTC has been carried out for almost thirty years, mainly aiming at the deterministic systems. However, in practical system, there are various stochastic disturbances such as sensor noises, stochastic disturbance, fault, or random variation of parameters. Thus, the system representation should be based on various stochastic models as the starting point of FDD and FTC. To enhance the reliability of the practical stochastic system, the research of FDD and FTC of stochastic systems has long been one of the important areas of control theory and application [1–4]. Usually, two kinds of approaches can be used to deal with the FDD problems of stochastic systems. The first group of methods is originated from the statistic theory, where the ratio of likelihood and Bayesian methods are used to estimate the abrupt changes of the states or parameters of the concerned systems. In this context, the FDD algorithms are obtained by using some numerical computation methods such as the Monte Carlo methods or the particle filtering methods [5–8]. As for the second group of methods, observer and filtering design theory are employed where the min-max optimization techniques have been applied to the estimation error system in order to guarantee some of the required performances [9, 10]. Passive method and active method are used to deal with the FTC problem of stochastic systems. For the passive fault tolerant control of stochastic systems, the related literatures are mainly concerning the robust fault tolerant control of stochastic systems [11–13]. For the active fault tolerant control of stochastic systems, the fault estimation information is used to carry out the controller reconfiguration [14–16].There are many systems in practice where the output concerned is the probability density function (PDF) of the system output, rather than the actual output values. For such a group of systems, the output that can be used for the feedback control is the measured output PDFs. Such types of stochastic systems are called stochastic distribution systems which were defined in [17], for instance, the paper evenness control, chemical polymerization process, grinding process fineness distribution control, and the flame distribution control. These kinds of system are described by the relationship between the input and the output PDF instead of the normal relationship between the input and the output. Stochastic distribution systems can not only represent Gaussian systems, but also represent non-Gaussian systems, making the assumption of being subject to Gaussian distribution be deleted. The aim of FDD in SDC systems is to use the measured input and output PDFs to obtain the information of the fault aim [18–20]. An observer based algorithm was previously developed in [21] to detect the fault in the non-Gaussian SDC system. For most stochastic distribution systems, the research results focus on the FDD of such systems and little attention is paid to the FTC of such systems. When the objective probability density function is known, the purpose of FTC of the stochastic distribution systems is to design a fault tolerant controller to make the postfault PDF still track the given distribution [14–16]. When the objective PDF is not known in advance, the purpose of FTC of stochastic distribution systems is to design a fault tolerant controller to make the uncertainty of the system output still be minimized. In this case, the concept of entropy is introduced to the stochastic distribution system and the performance index about the entropy is minimized considering the effect of the fault, leading to the minimum entropy fault tolerant control of stochastic distribution systems.In the above-mentioned SDC systems, only dynamic links between the inputs and the weights are considered. However, in practice, some algebraic relations also exist between the input and the weights, leading to a singular state space model between the weights and the control input. Such systems are called singular stochastic distribution control systems. The final purpose of the grinding processes is to obtain ore pulp with certain fineness distribution by grinding the raw ore, which is a typical stochastic distribution system. The main control system of the grinding process, such as the amount of ore, amount of water, two times amount of water, sump level, cyclone feed concentration, pressure, and rate of flow, has direct influence on the probability density function of fineness distribution. Such control systems can be regarded as the input control variables which control the system distribution probability density function. It shows that a singular stochastic distribution control model between the probability density function and the multilevel control loops composed of the above-mentioned main control systems can be built. In practice, time delays may exist in system modelling, the data (or image) processing, and transmission phases [22]. This is true especially in chemical processes and paper making processes, where ignorance of time delays will make the related FDD and FTC approaches totally invalid. Thus, the singular stochastic distribution system model should include the time delay factor, leading to singular time-delayed stochastic distribution system model. Few literatures have been reported on the FDD and FTC of such singular time-delayed stochastic distribution system. This forms the main purpose of the work in this paper. In this context, an iterative learning (ILO) based fault diagnosis algorithm is proposed to estimate the change of fault. Using the fault estimation information and switching control theory, a fault tolerant controller is designed to make the postfault PDF still track the given distribution.The rest of this paper is organized as follows. Section2 presents model description. In Section 3, the iterative learning observer based fault diagnosis algorithm is proposed. Section 4 gives the design of the fault tolerant control. Simulation results of FD and FTC are presented in Section 5, followed by some concluding remarks in Section 6.


## 2. Model Description

Denoteγ
(
y
,
u
(
t
)
) as the probability density function of the system output y being defined on a known bounded interval [
a
,
b
]. The model of the continuous singular time-delayed stochastic distribution system is given as follows:
(1)
E
x
˙
(
t
)
=
A
x
(
t
)
+
A
d
x
(
t
-
τ
(
t
)
)
+
B
u
(
t
)
+
N
F
(
t
)
,
V
(
t
)
=
D
x
(
t
)
,
γ
(
y
,
u
(
t
)
)
=
C
(
y
)
V
(
t
)
+
T
(
y
)
,
where x
∈
R
n is the system state vector, u
(
t
)
∈
R
m is the system input, V
(
t
)
∈
R
n
-
1 is the weight vector, F
∈
R
m is the system fault vector, and τ
(
t
) is the time delay term. A
∈
R
n
×
n, A
d
∈
R
n
×
n, B
∈
R
n
×
m, D
∈
R
(
n
-
1
)
×
n, E
∈
R
n
×
n, and N
∈
R
n
×
m are system parameter matrices with rank
⁡
(
E
)
=
q
<
n (i.e., E is a singular matrix). C
(
y
), T
(
y
) are decided by the selected basis functions.Assumption 1.
The system is regular; that is,|
s
E
-
A
|
≠
0.Assumption 2.
The system has no impulse; that is,rank
⁡
E
=
deg
⁡
|
s
E
-
A
|.If the above-mentioned assumption conditions are satisfied, there exist two nonsingular matricesQ and P such that
(2)
Q
E
P
=
[
I
q
0
0
0
]
,
Q
A
P
=
[
A
1
0
0
I
n
-
q
,
]
,
Q
A
d
P
=
[
A
d
1
0
0
I
n
-
q
,
]
,
where Q
,
P
∈
R
n
×
n, A
1
∈
R
q
×
q, A
d
1
∈
R
q
×
q, and I
i is an i-order identity matrix:
(3)
Q
B
=
[
B
1
B
2
]
.
D
P
=
[
D
1
,
D
2
]
,
D
1
∈
R
(
n
-
1
)
×
q
,
D
2
∈
R
(
n
-
1
)
×
(
n
-
q
)
.By applying the following state coordinate transformation(4)
x
(
t
)
=
P
[
x
1
(
t
)
x
2
(
t
)
]
,
where x
1
(
t
)
∈
R
q
×
1, x
2
(
t
)
∈
R
(
n
-
q
)
×
1, it can be obtained from (1) that
(5)
x
˙
1
(
t
)
=
A
1
x
1
(
t
)
+
A
d
1
x
1
(
t
-
τ
(
t
)
)
+
B
1
u
(
t
)
+
N
1
F
(
t
)
,
x
2
(
t
)
=
-
B
2
u
(
t
)
-
N
2
F
(
t
)
,
V
(
t
)
=
D
1
x
1
(
t
)
+
D
2
x
2
(
t
)
,
γ
(
y
,
u
(
t
)
)
=
C
(
y
)
V
(
t
)
+
T
(
y
)
=
C
(
y
)
[
D
1
x
1
(
t
)
+
D
2
x
2
(
t
)
]
+
T
(
y
)
.


## 3. Fault Diagnosis

The iterative learning observer based method is used to carry out fault diagnosis for the singular time-delayed stochastic distribution system, using the adaptive tuning property of the iterative learning observer. The fault diagnosis is constructed for the transformed system (5) as follows:
(6)
x
^
˙
1
m
(
t
)
=
A
1
x
^
1
m
(
t
)
+
B
1
u
(
t
)
+
A
d
x
^
1
m
(
t
-
τ
(
t
)
)
+
L
ξ
m
+
v
(
t
)
,
x
^
˙
2
m
(
t
)
=
-
B
2
u
(
t
)
,
v
(
t
)
=
K
1
v
(
t
-
g
)
+
K
2
ξ
m
(
t
-
g
)
,
F
^
˙
=
W
v
(
t
)
,
γ
m
(
y
,
u
(
t
)
)
=
C
(
y
)
[
D
1
x
^
˙
1
m
(
t
)
+
D
2
x
^
2
m
(
t
)
]
+
T
(
y
)
,
where x
^
1
m and x
^
2
m are the state estimation, g is the time delay of the ILO, ξ
m
(
t
) is the measured residual signal, L and K
i
(
i
=
1,2
) are the gain matrices, and v
(
t
) is the input signal of the ILO. From the third equation of the fault diagnosis observer, it can be seen that v
(
t
) is updated by the previous residual signal and the input information. The concept of “iterative” shows that the repeated operation can be carried out.Denote the observation error vector as follows:(7)
e
m
(
t
)
=
x
1
(
t
)
-
x
^
1
m
(
t
)
.
Then the residual signal can be formulated to give
(8)
ξ
m
(
t
)
=
∫
a
b
‍
σ
(
y
)
[
γ
(
y
,
u
(
t
)
)
-
γ
m
(
y
,
u
(
t
)
)
]
d
y
=
∫
a
b
‍
σ
(
y
)
C
(
y
)
d
y
[
D
1
x
1
(
t
)
-
x
^
1
m
(
t
)
+
D
2
N
2
F
]
=
Σ
D
1
e
m
(
t
)
+
Σ
D
2
N
2
F
,
where Σ
=
∫
a
b
‍
σ
(
y
)
C
(
y
)
d
y and σ
(
y
) is the function of y defined on [
a
,
b
]. If σ
(
y
)
=
1, ε
m is not used as γ
(
y
,
u
(
t
)
) and γ
m
(
y
,
u
(
t
)
) are integrated into 1 on [
a
,
b
] simultaneously. From the transformed stochastic distribution system (5) and fault diagnosis observer (6), the observation error dynamic system can be obtained as
(9)
e
˙
m
(
t
)
=
x
˙
1
(
t
)
-
x
^
˙
m
(
t
)
=
A
1
e
m
(
t
)
+
N
1
F
-
L
ξ
m
-
v
(
t
)
+
A
d
e
m
(
t
-
τ
(
t
)
)
,
where G
=
N
1
-
L
Σ
D
2
N
2. It can be further formulated to give
(10)
e
˙
m
(
t
)
=
(
A
1
-
L
Σ
D
1
)
e
m
(
t
)
+
G
F
-
v
(
t
)
+
A
d
e
m
(
t
-
τ
(
t
)
)
.
It is assumed that (
A
1
,
Σ
D
1
) is observable.Lemma 3.
If the input of the iterative learning observer is defined by (6), the following inequality holds:
(11)
v
T
v
≤
3
v
T
(
t
-
g
)
K
1
T
K
1
v
(
t
-
g
)
+
3
e
m
T
(
t
-
g
)
(
K
2
Σ
D
1
)
T
(
K
2
Σ
D
1
)
e
m
(
t
-
g
)
+
3
(
K
2
Σ
D
2
N
2
F
)
T
(
K
2
Σ
D
2
N
2
F
)
.Assumption 4.
Assume that whent
≥
0, the input u
(
t
) and fault of the system are bounded:
(12)
∥
u
∥
≤
b
u
,
∥
F
∥
≤
f
d
.Theorem 5.
It is assumed that Assumption4 is satisfied and there exist positive definite matrices P
1, Q
1, and R
1 to make the following equation hold:
(13)
(
A
1
-
L
Σ
D
1
)
T
P
1
+
P
1
(
A
1
-
L
Σ
D
1
)
+
R
1
+
P
1
P
1
=
-
Q
1
.Then the observation error system (10) is stable and observation error e
(
t
) is bounded.Proof.
The Lyapunov function is selected as follows:(14)
π
=
e
m
T
P
1
e
m
+
∫
t
-
g
t
‍
e
m
T
(
θ
)
R
1
e
m
(
θ
)
d
θ
+
∫
t
-
g
t
‍
v
T
(
α
)
v
(
α
)
d
α
+
η
∫
t
-
τ
(
t
)
t
‍
e
m
T
(
s
)
e
m
(
s
)
d
s
.
The first order derivative of the Lyapunov function π is derived as follows:
(15)
π
˙
≤
e
m
T
(
(
A
1
-
L
Σ
D
1
)
T
P
1
+
P
1
(
A
1
-
L
Σ
D
1
)
+
R
1
+
P
1
P
1
)
e
m
+
2
f
d
∥
P
1
G
∥
∥
e
m
∥
+
V
T
(
t
-
g
)
×
(
(
6
+
3
σ
)
(
K
2
Σ
D
1
)
T
(
K
2
Σ
D
1
)
-
R
1
)
e
m
(
t
-
g
)
-
σ
v
T
v
+
(
6
+
3
σ
)
(
K
2
Σ
D
2
N
2
F
)
T
(
K
2
Σ
D
2
N
2
F
)
+
2
T
2
∥
V
∥
∥
e
m
∥
∥
P
1
L
Σ
∥
∥
D
∥
∥
e
m
∥
∥
E
∥
.
It is known that the following inequality holds:2
∥
e
m
T
P
1
∥
∥
v
(
t
)
∥
≤
e
m
T
P
1
P
1
e
m
+
v
T
v.
It can be further obtained that(16)
π
˙
≤
2
f
d
λ
max
⁡
(
P
1
G
)
∥
e
m
∥
-
(
λ
min
⁡
(
Q
1
)
-
2
T
2
∥
P
1
L
Σ
∥
∥
D
∥
∥
E
∥
)
∥
e
m
∥
2
+
(
6
+
3
σ
)
(
K
2
Σ
D
2
N
2
F
)
T
(
K
2
Σ
D
2
N
2
F
)
-
σ
v
T
v
.
When the condition λ
min
⁡
(
Q
1
)
≥
2
T
2
∥
P
1
L
Σ
∥
∥
D
∥
/
∥
E
∥ is satisfied, the observation error system is asymptotically stable.


## 4. Fault Tolerant Control

When the fault is diagnosed, the fault tolerant control should be considered. Using the fault estimation information, the fault tolerant controller is designed to compensate the system performance loss caused by the fault. In this section, the performance index and the switching rules of the switched singular time delay system are introduced. When the fault happened, the performance index is minimized combining with the system model and the fault estimation information, leading to fault tolerant control of the singular time-delayed stochastic distribution system.


### 4.1. Design of the Faultless Controller

In this section, the purpose of the faultless PDF tracking controller is to choose a suitable control input to make the output probability density function track the given distribution. Different from the ordinary stochastic distribution system,E before x
˙ is a singular matrix in the dynamic model.The performance index of the switching time-delayed singular system is given as follows:(17)
J
1
(
u
(
t
)
)
=
∫
a
b
‍
(
γ
(
y
1
,
u
(
t
)
)
-
g
(
y
)
)
2
d
y
+
u
T
(
t
)
R
2
u
(
t
)
=
∫
0
+
∞
‍
(
x
T
(
t
)
Q
2
x
(
t
)
+
u
σ
(
t
)
T
(
t
)
R
2
u
σ
(
t
)
(
t
)
)
d
t
,
where Q
2
=
Q
2
T
∈
R
m
×
m
>
0 and R
2
=
R
2
T
∈
R
m
×
m
>
0 is the positive definite weight matrix of restricting the control input u
1
(
t
). u
σ
(
t
)
=
G
t
x
(
t
) is the control rule of the switching time-delayed singular system.It is further formulated to give that(18)
J
=
∫
0
+
∞
‍
(
x
T
(
t
)
Q
2
x
(
t
)
+
u
σ
(
t
)
T
R
2
u
σ
(
t
)
(
t
)
)
d
t
=
∑
i
=
1
m
‍
∑
j
=
1
∞
‍
∫
t
i
j
t
i
j
+
1
‍
(
x
T
(
t
)
Q
2
x
(
t
)
+
u
i
T
(
t
)
R
2
u
i
(
t
)
)
d
t
=
∑
i
=
1
m
‍
∑
j
=
1
∞
‍
∫
t
i
j
t
i
j
+
1
‍
x
T
(
t
)
(
Q
2
+
G
i
T
R
2
G
i
)
x
(
t
)
d
t
<
-
∑
k
=
0
∞
‍
∫
t
k
t
k
+
1
‍
V
˙
(
x
t
)
d
t
=
V
(
x
(
t
0
)
)
=
Φ
T
(
0
)
P
2
T
E
Φ
(
0
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
=
Φ
1
T
(
0
)
P
2
Φ
1
(
0
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
,
where Z
∈
R
m
×
m is a symmetrically positive definite matrix.From the above-mentioned inequality, it is known that(19)
J
≤
Φ
1
T
(
0
)
P
2
Φ
1
(
0
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
,
where Φ
(
t
) is the continuous system initial function, Φ
(
t
)
=
[
Φ
1
(
t
)
,
Φ
2
(
t
)
], Φ
1
(
t
)
∈
R
r, Φ
2
(
t
)
∈
R
n
-
r, and P
2
=
P
2
T
>
0.DenoteJ
* as the upper bound of the performance function and it is shown that
(20)
J
*
=
Φ
1
T
(
0
)
P
2
Φ
1
(
0
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
.Using the switching rule(21)
σ
(
t
)
=
argmin
{
x
T
(
γ
(
t
)
+
P
2
T
A
σ
(
t
)
Z
-
1
A
σ
(
t
)
T
P
2
)
x
}
,
γ
t
=
P
2
T
A
σ
(
t
)
+
A
σ
(
t
)
T
P
2
+
Z
+
Q
2
+
G
t
T
R
2
G
t
.
The Lyapunov function is selected as follows:
(22)
V
(
t
)
=
X
T
(
t
)
P
2
T
E
x
(
t
)
+
∫
t
-
τ
t
‍
x
T
(
s
)
Z
x
(
s
)
d
s
=
Φ
1
T
(
0
)
P
2
Φ
1
(
0
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
=
J
*
.
The whole closed-loop system can be written as
(23)
E
x
˙
(
t
)
=
A
σ
(
t
)
x
(
t
)
+
A
γ
σ
(
t
)
x
(
t
-
τ
)
+
B
σ
(
t
)
u
σ
(
t
)
(
t
)
,
x
(
t
)
=
Φ
(
t
)
,
t
∈
[
-
τ
,
0
]
,
where A
γ
σ
(
t
), A
σ
(
t
), and B
σ
(
t
) are the constant matrices with appropriate dimensions under the designed switching rule σ
(
t
) of the system.The upper bound of the performance function is obtained as(24)
J
*
=
J
1
(
u
(
t
)
)
=
∫
a
b
‍
(
g
-
(
y
,
t
)
-
H
(
y
)
u
(
t
)
)
2
d
y
+
u
T
(
t
)
R
2
u
(
t
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
=
∫
a
b
‍
g
-
2
(
y
,
t
)
d
y
-
[
∫
a
b
‍
g
-
(
y
,
t
)
γ
(
y
)
d
y
]
u
(
t
)
-
u
T
(
t
)
∫
a
b
‍
H
T
(
y
)
g
-
(
y
,
t
)
d
t
+
u
T
(
t
)
∫
a
b
‍
H
T
(
y
)
H
(
y
)
d
y
u
(
t
)
+
u
T
(
t
)
R
2
u
(
t
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
,
where s
=
∫
a
b
‍
H
T
(
y
)
H
(
y
)
d
y, H
(
y
)
=
C
(
y
)
D
2
B
2, and g
-
(
y
,
t
)
=
C
(
y
)
D
1
x
1
(
t
)
+
T
(
y
)
-
g
(
y
).Then the controller can be formulated to give(25)
u
(
t
)
=
(
s
+
R
2
)
-
1
∫
a
b
‍
H
T
(
y
)
g
-
(
y
,
t
)
d
y
-
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
.


### 4.2. Design of the Fault Tolerant Controller after Fault Occurs

When the fault happens in the system, a new performance index function is given as(26)
J
2
(
u
′
(
t
)
)
=
∫
a
b
‍
(
γ
(
y
,
u
′
(
t
)
)
-
g
(
y
)
)
2
d
y
+
u
′
(
t
)
R
2
u
′
(
t
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
,
where u
′
(
t
) is the controller input after fault occurs. It can be obtained that
(27)
γ
(
y
,
u
′
(
t
)
)
-
g
(
y
)
=
C
(
y
)
D
1
x
1
(
t
)
+
C
(
y
)
D
2
x
2
(
t
)
+
T
(
y
)
-
g
(
y
)
=
C
(
y
)
D
1
x
1
(
t
)
+
C
(
y
)
D
2
N
2
F
(
t
)
-
C
(
y
)
D
2
B
2
u
′
(
t
)
+
T
(
y
)
-
g
(
y
)
=
g
˘
(
y
,
t
)
-
H
(
y
)
u
′
(
t
)
,
where g
˘
(
y
,
t
)
=
C
(
y
)
D
1
x
^
1
(
t
)
+
C
(
y
)
D
2
N
2
F
^
(
t
)
+
T
(
y
)
-
g
(
y
).The fault tolerant controlleru
′
(
t
) can be obtained by solving the following optimization problem using the defined performance index and the system model:
(28)
min
⁡
x
,
w
,
z
,
a
,
M
a
+
tr
⁡
(
M
)
s
.
t
.
{
∑
t
=
1
m
‍
a
1
[
*
]
<
0
[
-
a
Φ
1
T
(
0
)
Φ
1
(
0
)
-
x
1
]
<
0
[
-
M
N
′
T
N
′
-
Z
]
<
0
,
where ∫
-
τ
0
‍
Φ
(
s
)
Φ
T
(
s
)
d
s
=
N
′
N
′
T, N
′ and M are symmetrical positive definite matrix:
(29)
[
*
]
=
[
Ξ
i
x
T
x
T
Ψ
i
T
A
τ
i
Z
x
-
A
0
0
0
x
0
-
Q
-
1
0
0
Ψ
i
0
0
-
R
-
1
0
Z
A
τ
i
T
0
0
0
-
Z
]
<
0
.
Ψ
i is a known parameter matrix.It can be formulated that(30)
J
*
=
k
(
constant
)
.Similar tou
(
t
), the following fault tolerant controller can be obtained as
(31)
u
′
(
t
)
=
(
s
+
R
2
)
-
1
∫
a
b
‍
H
T
(
y
)
g
˘
(
y
,
t
)
d
y
+
(
s
+
R
2
)
-
1
∫
a
b
‍
H
T
(
y
)
C
(
y
)
D
1
d
y
+
(
s
+
R
2
)
-
1
∫
a
b
‍
H
T
(
y
)
C
(
y
)
D
2
N
2
d
y
F
^
(
t
)
+
(
s
+
R
2
)
-
1
∫
a
b
‍
H
T
(
y
)
[
T
(
y
)
-
g
(
y
)
]
d
y
-
∫
-
τ
0
‍
Φ
T
(
s
)
Z
-
1
Φ
(
s
)
d
s
.Remark 6.
Based on the non-Gaussian singular stochastic distribution model of literature [16], time delay factor is considered in this paper. Similar to the fault diagnosis algorithm of literature [16], the iterative learning observer based fault diagnosis algorithm is used to carry out fault diagnosis for non-Gaussian time delay singular stochastic distribution system. Combined with the switch control theory, the corresponding fault tolerant control law is obtained.


## 4.1. Design of the Faultless Controller

In this section, the purpose of the faultless PDF tracking controller is to choose a suitable control input to make the output probability density function track the given distribution. Different from the ordinary stochastic distribution system,E before x
˙ is a singular matrix in the dynamic model.The performance index of the switching time-delayed singular system is given as follows:(17)
J
1
(
u
(
t
)
)
=
∫
a
b
‍
(
γ
(
y
1
,
u
(
t
)
)
-
g
(
y
)
)
2
d
y
+
u
T
(
t
)
R
2
u
(
t
)
=
∫
0
+
∞
‍
(
x
T
(
t
)
Q
2
x
(
t
)
+
u
σ
(
t
)
T
(
t
)
R
2
u
σ
(
t
)
(
t
)
)
d
t
,
where Q
2
=
Q
2
T
∈
R
m
×
m
>
0 and R
2
=
R
2
T
∈
R
m
×
m
>
0 is the positive definite weight matrix of restricting the control input u
1
(
t
). u
σ
(
t
)
=
G
t
x
(
t
) is the control rule of the switching time-delayed singular system.It is further formulated to give that(18)
J
=
∫
0
+
∞
‍
(
x
T
(
t
)
Q
2
x
(
t
)
+
u
σ
(
t
)
T
R
2
u
σ
(
t
)
(
t
)
)
d
t
=
∑
i
=
1
m
‍
∑
j
=
1
∞
‍
∫
t
i
j
t
i
j
+
1
‍
(
x
T
(
t
)
Q
2
x
(
t
)
+
u
i
T
(
t
)
R
2
u
i
(
t
)
)
d
t
=
∑
i
=
1
m
‍
∑
j
=
1
∞
‍
∫
t
i
j
t
i
j
+
1
‍
x
T
(
t
)
(
Q
2
+
G
i
T
R
2
G
i
)
x
(
t
)
d
t
<
-
∑
k
=
0
∞
‍
∫
t
k
t
k
+
1
‍
V
˙
(
x
t
)
d
t
=
V
(
x
(
t
0
)
)
=
Φ
T
(
0
)
P
2
T
E
Φ
(
0
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
=
Φ
1
T
(
0
)
P
2
Φ
1
(
0
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
,
where Z
∈
R
m
×
m is a symmetrically positive definite matrix.From the above-mentioned inequality, it is known that(19)
J
≤
Φ
1
T
(
0
)
P
2
Φ
1
(
0
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
,
where Φ
(
t
) is the continuous system initial function, Φ
(
t
)
=
[
Φ
1
(
t
)
,
Φ
2
(
t
)
], Φ
1
(
t
)
∈
R
r, Φ
2
(
t
)
∈
R
n
-
r, and P
2
=
P
2
T
>
0.DenoteJ
* as the upper bound of the performance function and it is shown that
(20)
J
*
=
Φ
1
T
(
0
)
P
2
Φ
1
(
0
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
.Using the switching rule(21)
σ
(
t
)
=
argmin
{
x
T
(
γ
(
t
)
+
P
2
T
A
σ
(
t
)
Z
-
1
A
σ
(
t
)
T
P
2
)
x
}
,
γ
t
=
P
2
T
A
σ
(
t
)
+
A
σ
(
t
)
T
P
2
+
Z
+
Q
2
+
G
t
T
R
2
G
t
.
The Lyapunov function is selected as follows:
(22)
V
(
t
)
=
X
T
(
t
)
P
2
T
E
x
(
t
)
+
∫
t
-
τ
t
‍
x
T
(
s
)
Z
x
(
s
)
d
s
=
Φ
1
T
(
0
)
P
2
Φ
1
(
0
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
=
J
*
.
The whole closed-loop system can be written as
(23)
E
x
˙
(
t
)
=
A
σ
(
t
)
x
(
t
)
+
A
γ
σ
(
t
)
x
(
t
-
τ
)
+
B
σ
(
t
)
u
σ
(
t
)
(
t
)
,
x
(
t
)
=
Φ
(
t
)
,
t
∈
[
-
τ
,
0
]
,
where A
γ
σ
(
t
), A
σ
(
t
), and B
σ
(
t
) are the constant matrices with appropriate dimensions under the designed switching rule σ
(
t
) of the system.The upper bound of the performance function is obtained as(24)
J
*
=
J
1
(
u
(
t
)
)
=
∫
a
b
‍
(
g
-
(
y
,
t
)
-
H
(
y
)
u
(
t
)
)
2
d
y
+
u
T
(
t
)
R
2
u
(
t
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
=
∫
a
b
‍
g
-
2
(
y
,
t
)
d
y
-
[
∫
a
b
‍
g
-
(
y
,
t
)
γ
(
y
)
d
y
]
u
(
t
)
-
u
T
(
t
)
∫
a
b
‍
H
T
(
y
)
g
-
(
y
,
t
)
d
t
+
u
T
(
t
)
∫
a
b
‍
H
T
(
y
)
H
(
y
)
d
y
u
(
t
)
+
u
T
(
t
)
R
2
u
(
t
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
,
where s
=
∫
a
b
‍
H
T
(
y
)
H
(
y
)
d
y, H
(
y
)
=
C
(
y
)
D
2
B
2, and g
-
(
y
,
t
)
=
C
(
y
)
D
1
x
1
(
t
)
+
T
(
y
)
-
g
(
y
).Then the controller can be formulated to give(25)
u
(
t
)
=
(
s
+
R
2
)
-
1
∫
a
b
‍
H
T
(
y
)
g
-
(
y
,
t
)
d
y
-
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
.


## 4.2. Design of the Fault Tolerant Controller after Fault Occurs

When the fault happens in the system, a new performance index function is given as(26)
J
2
(
u
′
(
t
)
)
=
∫
a
b
‍
(
γ
(
y
,
u
′
(
t
)
)
-
g
(
y
)
)
2
d
y
+
u
′
(
t
)
R
2
u
′
(
t
)
+
∫
-
τ
0
‍
Φ
T
(
s
)
Z
Φ
(
s
)
d
s
,
where u
′
(
t
) is the controller input after fault occurs. It can be obtained that
(27)
γ
(
y
,
u
′
(
t
)
)
-
g
(
y
)
=
C
(
y
)
D
1
x
1
(
t
)
+
C
(
y
)
D
2
x
2
(
t
)
+
T
(
y
)
-
g
(
y
)
=
C
(
y
)
D
1
x
1
(
t
)
+
C
(
y
)
D
2
N
2
F
(
t
)
-
C
(
y
)
D
2
B
2
u
′
(
t
)
+
T
(
y
)
-
g
(
y
)
=
g
˘
(
y
,
t
)
-
H
(
y
)
u
′
(
t
)
,
where g
˘
(
y
,
t
)
=
C
(
y
)
D
1
x
^
1
(
t
)
+
C
(
y
)
D
2
N
2
F
^
(
t
)
+
T
(
y
)
-
g
(
y
).The fault tolerant controlleru
′
(
t
) can be obtained by solving the following optimization problem using the defined performance index and the system model:
(28)
min
⁡
x
,
w
,
z
,
a
,
M
a
+
tr
⁡
(
M
)
s
.
t
.
{
∑
t
=
1
m
‍
a
1
[
*
]
<
0
[
-
a
Φ
1
T
(
0
)
Φ
1
(
0
)
-
x
1
]
<
0
[
-
M
N
′
T
N
′
-
Z
]
<
0
,
where ∫
-
τ
0
‍
Φ
(
s
)
Φ
T
(
s
)
d
s
=
N
′
N
′
T, N
′ and M are symmetrical positive definite matrix:
(29)
[
*
]
=
[
Ξ
i
x
T
x
T
Ψ
i
T
A
τ
i
Z
x
-
A
0
0
0
x
0
-
Q
-
1
0
0
Ψ
i
0
0
-
R
-
1
0
Z
A
τ
i
T
0
0
0
-
Z
]
<
0
.
Ψ
i is a known parameter matrix.It can be formulated that(30)
J
*
=
k
(
constant
)
.Similar tou
(
t
), the following fault tolerant controller can be obtained as
(31)
u
′
(
t
)
=
(
s
+
R
2
)
-
1
∫
a
b
‍
H
T
(
y
)
g
˘
(
y
,
t
)
d
y
+
(
s
+
R
2
)
-
1
∫
a
b
‍
H
T
(
y
)
C
(
y
)
D
1
d
y
+
(
s
+
R
2
)
-
1
∫
a
b
‍
H
T
(
y
)
C
(
y
)
D
2
N
2
d
y
F
^
(
t
)
+
(
s
+
R
2
)
-
1
∫
a
b
‍
H
T
(
y
)
[
T
(
y
)
-
g
(
y
)
]
d
y
-
∫
-
τ
0
‍
Φ
T
(
s
)
Z
-
1
Φ
(
s
)
d
s
.Remark 6.
Based on the non-Gaussian singular stochastic distribution model of literature [16], time delay factor is considered in this paper. Similar to the fault diagnosis algorithm of literature [16], the iterative learning observer based fault diagnosis algorithm is used to carry out fault diagnosis for non-Gaussian time delay singular stochastic distribution system. Combined with the switch control theory, the corresponding fault tolerant control law is obtained.


## 5. A Computer Simulation Example

To demonstrate the effectiveness of the fault tolerant control algorithm, the output probability density functionγ
(
y
,
u
(
t
)
) is given as
(32)
γ
(
y
,
u
(
t
)
)
=
ω
1
ϕ
1
(
y
)
+
ω
2
ϕ
2
(
y
)
+
ω
3
ϕ
3
(
y
)
.The following B-spline basis functions are used to approximateγ
(
y
,
u
(
t
)
):
(33)
ϕ
1
(
y
)
=
(
1
6
y
3
+
3
2
y
2
+
9
2
y
+
9
2
)
f
1
hhhhhh
+
(
-
1
2
y
3
-
5
2
y
2
-
7
2
y
-
5
6
)
f
2
hhhhhh
+
(
1
2
y
3
+
1
2
y
2
-
1
2
y
+
1
6
)
f
3
hhhhhh
+
(
-
1
6
y
3
+
1
2
y
2
-
1
2
y
+
1
6
)
f
4
,
ϕ
2
(
y
)
=
(
1
6
y
3
+
y
2
+
2
y
+
4
3
)
f
2
hhhhhh
+
(
-
1
2
y
3
-
y
2
+
2
3
)
f
3
hhhhhh
+
(
1
2
y
3
-
y
2
+
2
3
)
f
4
hhhhhh
+
(
-
1
6
y
3
+
y
2
-
2
y
+
4
3
)
f
5
,
ϕ
3
(
y
)
=
(
1
6
y
3
+
1
2
y
2
+
1
2
y
+
1
6
)
f
3
hhhhhh
+
(
-
1
2
y
3
+
1
2
y
2
+
1
2
y
+
1
6
)
f
4
hhhhhh
+
(
1
2
y
3
-
5
2
y
2
+
7
2
y
-
5
6
)
f
5
hhhhhh
+
(
-
1
6
y
3
+
3
2
y
2
-
9
2
y
+
9
2
)
f
6
,
where f
i
(
i
=
1,2
,
3,4
,
5,6
) is defined as follows:
(34)
f
i
(
y
)
=
{
1
,
y
∈
[
i
-
4
,
i
-
3
]
,
0
,
otherwise
.
The static model of the probability density function is shown as
(35)
γ
(
y
,
u
(
t
)
)
=
C
(
y
)
V
(
t
)
+
T
(
y
)
.
The dynamic model of the singular time-delayed stochastic distribution system is given as
(36)
E
x
˙
(
t
)
=
A
x
(
t
)
+
A
d
x
(
t
-
τ
(
t
)
)
+
B
u
(
t
)
+
N
F
(
t
)
,
V
(
t
)
=
D
x
(
t
)
.
The main matrices of the dynamic model are given as
(37)
E
=
(
1.0108
0.5794
5.7745
3.8263
-
2.0132
13.7264
-
5.5036
6.8344
-
12.1285
)
,
D
=
(
0.4953
3.3627
-
0.7513
0.7674
0.8753
-
0.4944
)
,
A
=
(
-
1.5675
-
2.5243
0.5341
0.4014
1.4281
-
9.1249
1.9462
-
0.1041
4.8273
)
,
A
d
=
(
-
1.0054
-
2.7458
0
1.5782
0.1461
-
4.8423
0
-
0.5293
1.8249
)
,
B
=
(
-
3.5044
-
0.3042
-
8.2882
)
,
N
=
(
2.2470
6.4054
-
2.1299
)
.Select the following nonsingular transformation matrix:(38)
P
=
(
0.2019
-
0.1422
0.4830
0.1631
0.3211
0.2028
0.1976
0.1455
-
0.1049
)
,
Q
=
(
0.1020
0.1422
-
0.1690
0.3428
0.1860
0.4515
-
0.4155
0.2784
0.1172
)
.
The dynamic model of the time-delayed singular stochastic distribution system can be transformed as
(39)
x
˙
1
(
t
)
=
A
1
x
1
(
t
)
+
B
1
u
(
t
)
+
A
d
x
^
1
m
(
t
-
τ
(
t
)
)
+
L
ξ
m
+
V
(
t
)
,
x
^
˙
2
m
(
t
)
=
-
B
2
u
(
t
)
,
V
(
t
)
=
K
1
V
(
t
-
g
)
+
K
2
ξ
m
(
t
-
g
)
,
γ
m
(
y
,
u
(
t
)
)
=
C
(
y
)
[
D
1
x
1
(
t
)
+
D
2
x
2
(
t
)
]
+
T
(
y
)
.
The following parameter matrices and vectors can be obtained as
(40)
A
1
=
(
-
0.5
-
0.25
0.11
-
0.17
)
,
A
d
=
(
-
0.1
0
0
-
0.5
)
,
B
1
=
(
1
-
5
)
,
B
2
=
0.4
,
N
2
=
0.6
,
N
1
=
(
1.5
1
)
,
D
1
=
(
0.5
1.9
0.2
0.1
)
,
D
2
=
(
1
0.6
)
,
K
1
=
(
0.62
0.31
-
0.28
0.43
)
,
K
2
=
(
0.41
0.36
)
,
L
=
(
0.5
-
0.6
)
,
Φ
(
t
)
=
(
1
4
e
t
+
1
0
)
,
t
∈
[
-
1,0
]
,
Φ
(
0
)
=
(
1
4
e
0
)
≈
(
0.6796
0
)
,
M
=
(
1.5680
0
0
1.1459
)
,
Ψ
1
=
(
-
0.0362
-
0.3155
)
,
Ψ
2
=
(
-
0.0275
-
0.1961
)
.
In this paper, the fault diagnosis and fault tolerant control algorithms can apply to abrupt fault and gradual fault. To demonstrate the effectiveness of the algorithms, gradual-varying and abrupt fault are considered, respectively.(1) A Gradual-Varying Fault. A gradual-varying fault is considered as follows:
(41)
F
=
{
0
t
<
5
s
0.4
(
t
-
5
)
5
s
≤
t
<
15
s
5.2
-
0.26
(
t
-
10.3
)
t
≥
15
s
,
where the learning rate W is chosen as W
=
(
33
15
). Figure 1 gives the fault estimation result. It can be seen that the tracking effect is good enough to show the effectiveness of the fault diagnosis algorithm. Figure 2 represents the faultless PDF and Figure 3 gives the postfault output probability density function. It is shown that fault tolerant control algorithm is effective to the fast-varying fault.Figure 1
Fault estimation result.Figure 2
The faultless PDF.Figure 3
The postfault PDF.(2) An Abrupt Fault. It is assumed that an abrupt fault is given as follows:
(42)
F
=
{
0
t
<
15
s
0.8
t
≥
15
s
.The learning rateW is chosen as W
=
(
36
14
). Figure 4 gives the fault estimation result of the abrupt fault. Figure 5 gives the faultless output PDF and the postfault probability density function result is shown in Figure 6. The postfault PDF can still track the given distribution, leading to good fault tolerant control results to the abrupt fault.Figure 4
Fault estimation result.Figure 5
The faultless PDF.Figure 6
The postfault PDF.


## 6. Conclusions

For non-Gaussian singular stochastic distribution system, the time delay factor is considered. The iterative learning observer is used to carry out the fault diagnosis. The concept of switching system is introduced to the time delay singular stochastic distribution system. A new fault tolerant control algorithm is proposed combining switch with the time delay singular stochastic distribution system to make the postfault probability density function still track the given distribution. The active fault tolerant control of the time delay singular stochastic distribution system is realized. The computer simulation results show that the fault diagnosis and fault tolerant control algorithms can apply not only to the abrupt fault, but also to the gradual-varying fault.



---

*Source: 937583-2014-02-27.xml*