Article
Fault detection of networked control systems in stochastic environments
Transactions of the Institute of Measurement and Control 35(4) 521–530 Ó The Author(s) 2012 Reprints and permissions: sagepub.co.uk/journalsPermissions.nav DOI: 10.1177/0142331212459540 tim.sagepub.com
Junqi Yang1,2 and Fanglai Zhu1
Abstract This paper deals with the problem of fault detection for networked control systems in the stochastic environments. The packet dropout, networkinduced delay and out-of-order packets due to long delay are assumed to exist stochastically in the sensor-to-controller link and the controller-toactuator link. By assuming that the stochastic environment could be described as a homogeneous Markovian chain, the networked control system is modelled as a discrete-time, non-uniform sampling stochastic parameter-varying Markov jump system. Based on the model, the residual generator is developed and the problem of fault detection is formatted as a filtering problem. By using the theory of a Markovian jump linear system, a fault detection filter which makes the residual generation system stochastically stable is considered, and a prescribed disturbance attenuation level is satisfied. A sufficient condition to solve this problem is given in terms of the linear matrix inequalities. A numerical example shows that the proposed fault detection filter is sensitive to the fault but robust to exogenous disturbance.
Keywords Networked control systems, fault detection, linear matrix inequalities, packet dropout, stochastic processes
Introduction With the rapid development of communication networks, conventional control system architectures have been evolving to modern networked control, and a great amount of effort has been devoted to the problems of networked control systems (NCSs). NCSs are control systems in which controller and plant are connected via a communication channel. A NCS integrates information, communications, and control into a single system in which control loops are closed over computer networks. Compared with conventional control systems, the networked control has a number of advantages, e.g. reduced cost, high resource utilization, simple installation and maintenance, increased system agility and reduced system wiring. Thus, it has received increasing interest in recent years. An introduction to networked control systems can be found (Tian and Levy, 2008; Yue et al., 2009). However, owing to the limitation of the network resources, the introduction of a communication network also brings some new problems and challenges, such as network-induced delay, data packet dropout, network scheduling and quantization problems, which all will inevitably degrade the performance of the NCSs and even cause system instability, and make analysis and synthesis of NCSs complex. The researches on NCSs have become important as international investigation focused on the automatic control domain, such as the control problem (Peng et al., 2011; Wang and Yang, 2011), the stability problem and stabilization problem (Zhang and Fang, 2011) of NCSs, and so on. Among these researches focusing on NCSs, the fault diagnosis and fault-tolerant control of NCSs have received
extensive attention. A review of work in the fault diagnosis domain of NCSs can be seen (Fang et al., 2007; Aubrun et al., 2008). A fault detection (FD) approach to NCSs with network-induced and unknown input based on eigendecomposition, adaptive evaluation and adaptive thresholds is proposed (Wang et al., 2007). A new FD scheme for networked control systems subject to uncertain time-varying delay is discussed (Wang et al., 2008). The FD problems of NCSs with random packet dropout are developed (Wang et al., 2009a, 2009b). Peng et al. (2010) addressed the problem of the FD for linear time-invariant systems over data networks with limited network quality of service (QoS) and the FD error dynamic systems are transferred to Markov jumping systems. A set of Kronecker delta functions is employed to characterize the random measurement delays and the stochastic data missing phenomenon, and the FD filter is designed in solving a certain set of linear matrix inequalities (LMIs) (He et al., 2008, 2010). Li and Tang (2010) presented a new design approach to the fault diagnostor for NCSs by constructing a novel reduced-order state observer
1
College of Electronics and Information Engineering, Tongji University, Shanghai China 2 College of Electrical Engineering and Automation, Henan Polytechnic University, Jiaozuo China Corresponding author: Junqi Yang, College of Electronics and Information Engineering, Tongji University, Shanghai 201804, China. Email:
[email protected]
522 of the augmented system, which can diagnose faults without using residual to embody faults. To the best of the authors’ knowledge, much literature has been reported about FD problem for NCSs, but most of published results on fault detection of NCSs considers either delay or packet dropouts by employing robust FD methods (Ding, 2008). One category of existing work assumes that the statistics of delay or packet dropouts are known (Zheng et al, 2006; He et al., 2008; Wang et al., 2009b). Another category transforms the influence of delay or packet dropouts into model uncertainties and unknown inputs (Wang et al., 2008, 2009b). Only a few results consider random network environments with network delay, packet dropouts and out-of-order packets simultaneously. This motivates us to study this interesting and challenging problem, which has great potential in practical applications. In this paper, the problems of modeling and fault detection for NCSs in a stochastic environment are considered, and the effect of the disturbance on the NCSs is also considered. Then the FD scheme is proposed and makes the residual generation system stochastically stable, and satisfies a prescribed disturbance attenuation level. The rest of this paper is organized as follows. In the next section we give the system description and propose a mathematical model to describe the stochastic NCS. A FD filter is then considered. Next we give some main results and simulation results, which are followed by the conclusions.
Notation The notation used throughout the paper is standard. Here R and Rn3m denote n-dimensional Euclidean space and set of all n 3 m real matrices, respectively. For matrices or vectors, superscript T indicates transposition. For symmetric matrices, X . 0 (\ 0) indicates that X is positive definite (negative definite). In some partitioned symmetric matrices, the symbol generically denotes each of its symmetric blocks. We use l2 to denote the space of square summable sequences and E( ) to denote the expectation of ( ).
Transactions of the Institute of Measurement and Control 35(4) (
x_ (t) = Ax(t) + Bu(t) + Dd(t) + Ff (t) y(t) = Cx(t)
ð1Þ
where x(t) 2 Rn , u(t) 2 Rs , y(t) 2 Rh are the state vector, the control input vector and the measurement output vector, respectively. The unknown input d(t) 2 Rp denotes modelling errors, or exogenous disturbance input or random noise, f (t) 2 Rq is the directly immeasurable fault signal vector. Without loss of generality, the l2 norms of u(t), d(t) and f (t) are assumed to be existing and bounded and A, B, C, D, F are real constant matrices of appropriate dimensions. Assumption 1. The sensor and actuator are clock driven and synchronized, the sampling period is T, and all measurement data packets are time-stamped. The actuator has a logic zeroorder hold (ZOH). Assumption 2. The controller is event driven and triggered by the arrival of a measurement data packet. Once the control data become available, it sends a control packet to the actuator of the control loop. The control packet carries the timestamps of the measurement and control data. Assumption 3. For any control loop, the measurement data of the plant is transmitted with a single packet and so is the control data. Assumption 4. No new control command is generated if a data packet is dropped out since the controller is event driven. In this paper, the signal adopted by the actuator will keep the previous input if no new control commands are updated at the actuator.
NCS model and problem statement
Remark 1. Under the aforementioned assumptions, here is no distinction between packet dropouts that occur in the sensorto-controller connection and the controller-to-actuator connection in the network. Indeed, for packet dropouts between the sensor and the controller (as shown the interval [(l + 7)T, (l + 8)T] in Figure 1) no new control update is computed by the controller, thus no new control input is sent to the actuator. In the case of packet dropouts between the controller and the actuator (as shown the interval [(l + 3)T, (l + 4)T] in Figure 1), the packet cannot be received by the actuator either.
Consider a linear continuous time-invariant system as the controlled plant
Remark 2. Employing the time-stamping technique described in Assumptions 1 and 2, the logic ZOH at the actuator stores
Figure 1 Time diagram for data transmission.
Yang and Zhu
523
the latest control packet. This implies that the logic ZOH discards all control packets but the most recent valid one (as shown in Figure 1). The actuator keeps its control signal unchanged until the output of the logic ZOH gets updated to a new value (Peng et al., 2011). An example of the timing diagram of data transmission of the considered NCS in a stochastic environment is shown in Figure 1, in which the two control signals shown in dashed lines are lost, and some signals are not used by the actuator due to the long delay and being out of order (as shown by (l + 1) and (l + 5)). It can be seen from the timing diagram that the control inputs acting on the actuator are different from one sampling interval to another, and thus the system models of the NCS vary from one sampling interval to another as the packet dropout situations change. In Figure 1, the tk (k = 1, 2, 3, . . . ) are some integers such that ft1 , t2 , t3 , . . .g f1, 2, 3, . . .g, and denote the time constant of the updating packet in the actuator side. Suppose that the time instants corresponding to two successively updating data packet are tk T and tk+1 T , respectively, with tk+1 T tk T = uk T (1 uk N ), then the results caused by delay, packet dropouts and being out of order can be regarded as non-uniform sampling with sampling intervals fuk T g at the actuator (Suplin et al., 2007; Wang et al., 2009b), and uk can be obtained from the time-stamp of the received control packet. Similar to the work of Wang et al. (2009b), the dynamics of the plant (1) at the time tk can be written as the following discrete-time Markov jump linear model (
x(k + 1) = A(uk )x(k) + B(uk )u(k) + D(uk )d(k) + F(uk )f (k) y(k) = Cx(k) ð2Þ kP 1
where t0 = 0, tk = uj T (k 1), j=0 R R uk T At u T , B(uk ) = 0 e dtB D(uk ) = 0 k eAt dtD, A(uk ) = e R uk T At F(uk ) = 0 e dtF, x(k) = x(tk ), u(k) = u(tk ), y(k) = y(tk ), fuk g is a discrete homogeneous Markov chain taking values in a finite state space Y = f1, 2, . . . , N g with transition probability matrix P = ½lij i, j2Y , and lij is defined as for
k=0,1,2,.,
denotes the packet dropout of the sensor–controller channel of Wang et al. (2009b).
Fault detection filter design A typical FD system consists of a residual generator and a residual evaluation stage including an evaluation function and a threshold. For the purpose of residual generation, the following fault detection filter is used
where
N P
lij=1 and A(uk ), B(uk ), D(uk ), F(uk ) are known real
j=1
constant matrices for all uk = i 2 Y. The matrices Ai := A(uk = i), Bi := B(uk = i), Di := D(uk = i) and Fi := F(uk = i), are known real constant matrices of appropriate dimensions. Remark 3. Although Equation (2) is obtained in a similar manner to the work of Wang et al. (2009b), the meaning of uk T is different. In Equation (2), uk T denotes the interval between two successfully received control packets at the actuator, and describes random character between the sensor–controller channel and controller–actuator channel, while uk T is assumed as the interval between two successive successfully received measurements at the controller and just
ð3Þ
where ^x(k) 2 Rm is the state vector of the filter, yf (k) and rk are the input and output of the filter, respectively, the rk is the residual vector and Afi , Bfi , Cfi and Dfi are appropriately dimensioned filter matrices to be determined. Remark 4. The design idea of the FD filter is similar to the work of Yao et al. (2011). The filter of Yao et al. (2011) is designed to detect the faults of singular systems, while our paper deals with the problem of modelling and FD for networked control systems in stochastic environments. So the results of Yao et al. (2011) and our paper are different. Data packet dropout happens unavoidably in the network environment, so the measurements y(k) will drop randomly. In this paper, it is assumed that the data packet dropout is described by a Bernoulli stochastic variable ak with ak = 0 denoting the measurements dropout condition, ak = 1 the correct receiving condition. Here ak can be defined as Prfak = 1g = Efak g = a, Prfak = 0g = 1 Efak g = 1 a where a 2 R is a known positive scalar. The relationship between y(k) and yf (k) can be established, i.e.
A(uk T )
lij = Prfuk+1 = jjuk = ig
^x(k + 1) = Afi^x(k) + Bfi yf (k) rk = Cfi^x(k) + Dfi yf (k)
yf (k) = ak y(k)
ð4Þ
By the second expression of Equations (2), (3) and (4), the FD filter is to be the following form
^x(k + 1) = Afi^x(k) + ak Bfi Cx(k) rk = Cfi^x(k) + ak Dfi Cx(k)
ð5Þ
The objective of FD is to identify the fault f (k) when it appears. Similar to the cases presented by Gao et al. (2008), Zhong et al. (2005) and Yao et al. (2010), the introduction of a suitable weighting matrix G(s) is used to limit the frequency interval, in which the fault should be identified, and the system performance could be improved. Here we have the equation ~f (s) = G(s)f (s), where f (s) and ~f (s) denote the Laplace transforms of f (t) and ~f (t), respectively. One state space realization of ~f (s) = G(s)f (s) can be
~x(k + 1) = Aw~x(k) + Bw f (k) ~f (k) = Cw ~x(k) + Dw f (k)
ð6Þ
where ~x(k) 2 Rr is the state vector, and Aw , Bw , Cw and Dw are chosen correspondingly.
524
Transactions of the Institute of Measurement and Control 35(4)
T We use the notation h(k) = xT (k) ~xT (k) ^xT (k) , T w(k) = uT (k) d T (k) f T (k) , e(k) = rk ~f (k) and a ~k = ak a, here a (0 a 1) is a positive constant. Based on (2), (5) and (6), the residual system is given by (
and when w(k) = 0 the system is stochastically stable, if there exist matrices Pi . 0, i 2 Y such that the following LMIs 2
~ 1i + a ~ 2i )h(k) + B ~ i w(k) h(k + 1) = (A ~k A ~ 1i + a ~ 2i )h(k) + D ~ i w(k) e(k) = (C ~k C
Pi 6 0 6 6 C 6 ~ 1i 6 ~ 6 uC 2i 6 4 uP ~ 2i iA ~ 1i iA P
ð7Þ
where 2
3 0 0 Ai 7 ~ 1i = 6 Aw 0 5 , A 4 0 aBfi C 0 Afi 2 3 Bi Di Fi ~ i = 4 0 0 Bw 5 B 0 0 0 ~ 1i = ½ aDfi C C
Cw
~i = ½ 0 0 D
Dw
2
3 0 0 7 0 0 5, 0 0
0 ~ 2i = 6 A 4 0 Bfi C
I 0 0 0
I 0 0
3 7 7 7 7 7\0 7 7 5 i P
i P 0
ð13Þ
hold for i 2 Y, where N X
lij Pj , u =
pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi a(1 a),
j=1
sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi (sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi) ‘ ‘ X X T e (k)e(k) , kwk2 = w(k)T w(k) kek2, E = E
0 0 ,
k=0
In the residual evaluation stage, an evaluation function and a threshold should be provided for the purpose of FD. Here, as in most contributions, we adopt the following residual evaluation function ( k X
)12
Proof. An index is introduced as JD = EfhT (k + 1)P(uk+1 )h(k + 1)jh(k), uk g hT (k)P(uk ) h(k) + EfeT (k)e(k)g g2 w(k)T w(k)
ð8Þ
JD = EfhT (k + 1)P(uk + 1 )h(k + 1)jh(k), uk g hT (k)P(uk )h(k) + EfeT (k)e(k)g g2 w(k)T w(k) ~ 1i + a ~ 1i + a ~ 2i )h(k) + B ~ 2i )h(k) ~ i w(k))T Pi ((A ~k A ~k A = Ef((A
The threshold is selected as sup 06¼d(k), u(k)2l2 , f = 0
k=0
Then, along the residual system (7),
riT ri
i=1
Jth =
g 2 I ~i D 0 0 iB ~i P
i = P
~ 2i = ½ Dfi C Cfi , C
J (k) =
ð12Þ
kek2, E \ gkwk2
sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi r X riT ri
~ i w(k))g h(k)T Pi h(k) +B ~ 1i + a ~ 1i + a ~ 2i )h(k) + D ~ 2i )h(k) ~ i w(k))T ((C ~k C ~k C + Ef((C
ð9Þ
i=1
~ i w(k))g g 2 w(k)T w(k) +D ~ T Pi A ~ T Pi A ~ T Pi B ~ 1i h(k) + a ~ 2i h(k) + hT (k)A ~ i w(k) ~ k hT (k)A = EfhT (k)A 1i 1i 1i
where r denotes the maximum time step of the evaluation function. The decision logic can be defined as J (k) . Jth ,
fault and alarm for fault
J (k) Jth ,
no fault
~ T Pi A ~ T Pi A ~ 1i h(k) + a ~ 2i h(k) ~ 2k hT (k)A +a ~ k hT (k)A 2i 2i T T ~ Pi B ~ i w(k) +a ~ k h (k)A 2i
ð10Þ
~ 1i h(k) + a ~ 2i h(k) + wT (k)B ~ T Pi A ~ T Pi A ~ T Pi B ~ i w(k)g ~ k wT (k)B + wT (k)B i i i T T T~ T T~ ~ ~ ~ k h (k)C C 2i h(k) h(k) Pi h(k) + Efh (k)C C 1i h(k) + a 1i
1i
~TD ~ + hT (k)C 1i i w(k) T ~ ~T C ~TD ~ 1i h(k) + a ~ 2i h(k) + a ~ i w(k) +a ~ k h (k)C T C ~ 2 hT (k)C ~ k hT (k)C
Main results
2i
For formulating some practically computable criteria to obtain the filter parameters described by (3), the following definition is useful in deriving the criteria. Definition 1 (Yao et al., 2011). When w(k) = 0, system (7) is said to be stochastically Markovian jump stable, if there exists a piecewise quadratic Lyapunov function T
with Pi = P(uk = i) . 0, i 2 Y, and DV = EfV (h(k + 1), uk + 1 )g V (h(k), uk ) decreasing.
its is
difference negative
Theorem 1. Consider the system (7) with zero initial conditions, let g . 0 be a given scalar, then the residual error e(k) satisfies the following performance index
2i
2i
g2 w(k)T w(k)
ð14Þ Note the fact that Ef~ a2k g = a(1 a), we have
Ef~ ak g = Efak ag = 0
h(k) JD = w(k)
ð11Þ
V (h(k), uk ) = h (k)P(uk )h(k)
k
T ~ ~ ~ TC ~ TC ~T ~ ~ k wT (k)D + wT (k)D i 1i h(k) + a i 2i h(k) + w (k)Di Di w(k)g
T
h(k) Li w(k)
and
where Li = 2 4
~T C ~T ~ ~T ~ ~T ~ ~ Pi + C 1i 1i + A1i Pi A1i + a(1 a)A2i Pi A2i + a(1 a)C2i C 2i
~TD ~ T Pi B ~i + C ~ A 1i 1i i
~ 1i + D ~ ~ T Pi A ~ TC B i i 1i
~ ~ TD ~T ~ g 2 I + D i i + Bi Pi Bi
3 5
Yang and Zhu
525
Applying Schur complement to (13), we have Li \ 0, and which leads to EfhT (k + 1)P(uk+1 )h(k + 1)jh(k), uk g hT (k)P(uk )h(k) + EfeT (k)e(k)g g 2 w(k)T w(k) \ 0
ð15Þ
Summing up both sides of the above inequality from 0 to ‘, and considering zero initial condition and EfhT (k + 1)P(uk+1 )h(k + 1)gjk!‘ . 0 we have ( E
‘ X
) T
e (k)e(k)
k=0
g2
‘ X
w(k)T w(k) \ 0
k=0
so the condition (12) is satisfied. Next, we proved that system (7) is stochastically Markovian jump stable. Given the Lyapunov function (11), when w(k) = 0, we can obtain from the deducing process of Equation (14)
2
Qi 6 0 g 2 I 6 6 6 aD C C ~i D I C fi 6 fi w 6 6 0 0 C 0 uD fi 6 6 " # 6 C 0 uB fi 6 0 0 6 6 uBfi C 0 6 6" # " # 6 Ui A + aBfi C Afi Ui Bi i 4 0 Vi Bi Vi Ai + aBfi C Afi
~ T Pi A ~ T Pi A ~ 1i + a(1 a)A ~ 2i \ 0 Pi + A 1i 2i
Lemma 1 (Zhong et al., 2005). Consider system (7) and let g . 0 be a given scalar, then LMIs (13) with Pi . 0, i 2 Y are feasible, if and only if there exist matrices Pi . 0 and Wi such that the following LMIs
0 0 ~i Wi B
\0
i Wi WiT + P 0
I
0
W i W i + Qi
0
0
W i W i + Qi
T
T
7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 5
1 Bfi G = i Dfi 0
Afi Cfi
0 I
Afi C fi
Bfi Dfi
ð18Þ
3 7 7 7 7 7 7 7 7 5 i Wi WiT + P
ð16Þ
hold for i 2 Y. Theorem 2 Consider the residual system (7). Let m = n + r. For a given positive constant g, there exist Qi = Q1i Q2i Ui G i . 0 and , matrices Afi , Bfi , C fi W = i QT2i Q3i Vi Gi and Dfi satisfying the following LMIs for i 2 Y
P1i PT2i
P2i W1i , WiT = P3i W4i
W2i W3i
ð19Þ
where P1i 2 R(n+r) 3 (n+r) , P2i 2 R(n+r) 3 m , P3i 2 Rm 3 m , W1i 2 R(n+r) 3 (n+r) , W2i 2 R(n+r) 3 m , W3i 2 Rm 3 m and W4i 2 Rm 3 (n+r) . The W3i and W4i are non-singular. Define the following matrices
so DV \ 0. We can obtain that system (7) is stochastically Markovian jump stable. The proof is completed.u
I 0 0
N P 0 A where , Qi = lij Qj , j 2 Y, Ai = i Bi = 0 Aw j=1 B i Di F i , C = ½ C 0 and C w = ½ 0 Cw . There exists 0 0 Bw a FD filter in the form of (3) such that residual system (7) is stochastically stable with the guaranteed performance index g. Moreover, if the above conditions are feasible, then the matrices for a desired filter in the form of (3) are given
from (13), and applying the Schur complement to (13), we have
I 0 0 0
ð17Þ
Pi =
g 2 I ~i D
\0
~ T Pi A ~ T Pi A ~ 1i + a(1 a)A ~ 2i )h(k) = hT (k)( Pi + A 1i 2i
Pi 6 0 6 6 ~ 6 C 1i 6 6 uC ~ 2i 6 6 ~ 2i 4 uWi A ~ 1i Wi A
Proof. Inspired by Gao et al. (2008) and Yao et al. (2011), let the matrices Pi and WiT be partitioned as
DV = EfV (h(k + 1), uk + 1 )g V (h(k), uk )
2
3
Ji =
I 0 , 1 0 W3i W4i
Qi = JiT Pi Ji =
Q1i QT2i
Q2i Q3i
Thus Ji is non-singular. By pre- and post-multiplying the resulting expression with matrices diagfJiT , I, I, I, JiT , JiT g and diagfJi , I, I, I, Ji , Ji g to (16) respectively, we have 2
JiT Pi Ji 6 0 6 ~ 1i Ji 6 C 6 6 uC ~ 2i Ji 6 4 uJ T Wi A ~ 2i Ji i ~ 1i Ji JiT Wi A
g2 I ~ Di 0 0 T ~i Ji W i B
I 0 0 0
I 0 0
T i )Ji Ji ( Wi WiT + P 0
3 7 7 7 7 7 7 5 T T Ji ( Wi Wi + Pi )Ji
ð20Þ
\0 thus, the following equations can be obtained
~ 1i Ji = aDfi C C ~ 2i Ji = u½ Dfi C 0 , C Cfi W3i1 W4i , uC w T W1iT Bi ~ 2i Ji = u W4i Bfi C 0 , J T Wi B ~i = uJiT Wi A i W4iT Bfi C 0 W4iT W3iT W2iT Bi ~ 1i Ji = JiT Wi A
W1iT Ai + aW4iT Bfi C T T T W4i W3i W2i Ai + aW4iT Bfi C
W4iT Afi W3i1 W4i , W4iT Afi W3i1 W4i
526
Transactions of the Institute of Measurement and Control 35(4)
let
design of the FD filter (3) is formulated as an H‘ -filtering problem, which can be solved by the following optimization problem
W i = JiT Wi Ji =
=
I 0 1 0 W3i W4i W1iT T T W4i W3i W2iT
Qi = JiT Pi Ji =
N X
T
W1iT W2iT
W4iT W3iT
W4iT W3i1 W4i W4iT W3i1 W4i lij Qj ,
I 0 1 0 W3i W4i
min l,
(i = 1, . . . , N )
j=1
then (20) can be rewritten as
P1 P3
\0 P2
ð21Þ
here 2
3 Qi 2 5 0 P1 = 4 g I ~i D I aDfi C C w Cfi W3i1 W4i 2 3 I T 5 P2 = 4 0 W i W i + Qi T 0 0 W i W i + Qi
6 6 6 6 6 4
3 ½ uDfi C 0 0 0 7 uW4iT Bfi C 0 0 07 7 uW4iT Bfi C 0 7 7 T T T 1 T 5 W1i Bi W4i Afi W3i W4i W1i Ai + aW4i Bfi C 0 W4iT W3iT W2iT Ai + aW4iT Bfi C W4iT Afi W3i1 W4i W4iT W3iT W2iT Bi
define
Ui G i W4iT W3i1 W4i W1iT = Vi Gi W4iT W3iT W2iT W4iT W3i1 W4i T 1 Bfi W4i 0 Afi Bfi W3i W4i 0 = Cfi Dfi 0 I Dfi 0 I T 1 T W4i Afi W3i W4i W4i Bfi ð22Þ = Cfi W3i1 W4i Dfi
Wi = Afi C fi
from (21), one readily obtains (17). From Equation (22), the FD filter matrices Afi , Bfi , Cfi and Dfi in Equation (3) can be written as
Afi Cfi
#" # T " 1 Bfi W4i 0 Afi Bfi (W3i1 W4i ) 0 = Dfi 0 I C fi Dfi 0 I #" " #" 1 1 1 1 1 Afi Bfi (W3i W4i ) (W3i W4i ) Gi 0 = C fi Dfi 0 I 0
Remark 5. Based on the filter in the present paper, it is easy to resolve the fault-tolerant control problem by means of feedback control. When some faults are detected, we can reconfigure the control law to compensate for the faults based on the Markovian jump system model. Meanwhile, by using the guaranteed cost control approach according to the performance index (12) and online controller switching based on the system model, the fault-tolerant control strategies based on the present filter in our paper can be designed to ensure stability of the closed-loop networked control system at all times. The issue of fault-tolerant control may become the topic of our future work.
In this section, we illustrate the effectiveness of our proposed methods with a numerical example. We consider a linear timeinvariant system (1) with the following parameters described by Peng et al. (2010)
ð23Þ
Numerical example
P3 = 2
s:t: (17)
A=
0 3
1 2 0:1 2 ,B= ,D= ,F = ,C=½1 1 4 1 0:05 1
Assume that the sampling period of the random NCS is T = 0:01, the finite state space of the Markov chain is Y = f1, 2, 3g. By some calculation, the Markovian jump system model is obtained as (2) with the following parameters in different modes
0:9999 0:0098 0:0200 , B1 = , 0:0294 :9606 0:0095 0:0010 0:0200 D1 = , F1 = 0:0005 0:0095
A1 =
0:9994 0:0192 0:0400 , B2 = , 0:0577 0:9225 0:0200 0:0020 0:0400 D2 = , F2 = 0:0200 0:0010
A2 =
0 I
#
which implies that (W3i1 W4i )1 can be viewed as a similarity transformation on the state space realization of the filter and has no effect on the filter mapping form yf (k) to rk . Without loss of generality, we may set (W3i1 W4i )1 = I, thus leading to (18). The proof is complete. u From Theorem 2, the parameters of the FD filter (3) can be obtained. Based on the performance index (12), the
0:9987 0:0283 0:0598 , B3 = , 0:0848 0:8857 0:0316 0:0030 0:0598 D3 = , F3 = 0:0016 0:0316
A3 =
The transition probability matrix is given as follows (Peng et al., 2010; He et al., 2008) 2
3 0:5 0:4 0:1 P = 4 0:2 0:6 0:2 5 0:2 0:3 0:5
Yang and Zhu
527
Figure 2 Random switching mode.
Figure 3 Disturbance signal.
Figure 2 indicates the switching mode with random environment satisfying the Markov transition probability matrix P. The unknown input or disturbance signal is supposed to be randomly uniformly distributed over ½0:5 0:5, and Figure 3 shows the signal pattern. The fault signal f (k) is given as
f (k) =
2 + 0:5 cos (4pk) for k = 100, 101, . . . , 200 0 others
Figure 4 gives its corresponding wave shape. With the given parameters and based on the optimal filter design problem (23), Table 1 shows the minimum guaranteed
528
Transactions of the Institute of Measurement and Control 35(4)
Figure 4 Fault signal of the NCS.
Figure 5 Residual signal of the NCS.
performances g in the terms of the feasibility of (23) for different values of a, from which one can see that the smaller the value of the a, the larger the value of g. This is reasonable, as smaller a implies a higher chance of measurements missing, and thus worse disturbance attenuation performance g.
From Table 1, we can see that g = 0:1283 when a = 0:9, under this situation, Figure 5 shows the generated residual signal rk , and the evolution in (8) is presented in Figure 6. We sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 5P 00 select the threshold as Jth = sup riT ri , after the 500 times f =0
i=1
Yang and Zhu
529
Table 1. Minimum g for different values of a a
0.6
0.65
0.7
0.75
0.8
0.85
0.9
0.95
0.99
g
0.1337
0.1333
0.1327
0.1320
0.1312
0.1300
0.1283
0.1250
0.1152
Figure 6. Evaluation of J(k).
simulations, the value Jth = 0:4251 is obtained. The simulation result in Figure 6 shows that J (101) = 0:3734 and J (102) = 0:4823, and J (102) . Jth . Thus, the appeared fault can be detected after two time steps.
Conclusions In this paper, the discussion has been based on NCSs under stochastic environments. A FD approach for NCSs has been presented. Considering the stochastic characters of NCSs with packet dropout, network-induced delay, and out-of-order packets, we deal with them by modelling the NCSs as a discrete-time, non-uniform sampling stochastic parametervarying Markov jump system. Then by using a design method for the FD filter, which makes the residual generation system stochastically stable and a sufficient condition is obtained by solving this problem in terms of LMI. At last, a numerical simulation example is given to illustrate our results, and which can detect the occurrence of a fault in a very short time. There are some further research topics. Similar to linear systems, new fault diagnosis theory for nonlinear NCSs should also be developed.
Funding This work was supported by the National Natural Science Foundation of China (project number 61074009) and supported by Shanghai Leading Academic Discipline Project (project number B004).
References Aubrun C, Sauter D and Yame´ J (2008) Fault diagnosis of networked control systems. International Journal of Applied Mathematics and Computer Science 18: 525–537. Ding SX (2008) Model-based Fault Diagnosis Techniques: Design Schemes, Algorithms, and Tools. Berlin: Springer-Verlag. Fang H, Ye H and Zhong M (2007) Fault diagnosis of networked control systems. Annual Reviews in Control 31: 55–68. Gao H, Chen T and Wang L (2008) Robust fault detection with missing measurements. International Journal of Control 81: 804–819. He X, Wang Z, Ji YD and Zhou D (2008) Network-based fault detection for discrete-time state-delay systems: A new measurement model. International Journal of Adaptive Control and Signal Processing 22: 510–528. He X, Wang Z, Ji YD and Zhou D (2010) Fault detection for discrete-time systems in a networked environment. International Journal of Systems Science 41: 937–945.
530
Transactions of the Institute of Measurement and Control 35(4)
Li J and Tang G (2010) Fault diagnosis for networked control systems with delayed measurements and inputs. IET Control Theory and Applications 4: 1047–1054. Peng C, Tian Y and Yue D (2011) Output feedback control of discrete-time systems in networked environments. IEEE Transactions on Systems, Man, and Cybernetics – Part A: Systems and Humans 41: 185–190. Peng C, Yue D, Tian E and Gu Z (2010) Observer-based fault detection for networked control systems with network Quality of Services. Applied Mathematical Modelling 34: 1653–1661. Suplin V, Fridman E and Shaked U (2007) Sampled-data H‘ control and filtering: nonuniform uncertain sampling. Automatica 43: 1072–1083. Tian Y and Levy D (2008) Dealing with network complexity in realtime networked control. International Journal of Computer Mathematics 85: 1235–1253. Wang Y, Ding SX, Ye H and Wang G (2008) A new fault detection scheme for networked control systems subject to uncertain timevarying delay. IEEE Transactions on Signal Processing 56: 5258– 5268. Wang J and Yang H (2011) H‘ control of a class of networked control systems with time delay and packet dropout. Applied Mathematics and Computation 217: 7469–7477. Wang Y, Ye H, Ding SX and Wang G (2009a) fault detection of networked control systems subject to access constraints and random packet dropout. Acta Automatica Sinica 35: 1230–1234. Wang Y, Ye H, Ding SX, Wang G and Zhou D (2009b) Residual generation and evaluation of networked control
systems subject to random packet dropout. Automatica 45: 2427–2434. Wang Y, Ye H and Wang G (2007) Fault detection of NCS based on eigendecomposition, adaptive evaluation and adaptive threshold. International Journal of Control 80: 1903–1911. Yao X, Feng Y, Wu L and Zheng W (2010) Filter design for discretetime Markovian jump singular systems with its application to fault detection. In The 3rd International Symposium on Systems and Control in Aeronautics and Astronautics, Harbin, pp. 321–326. Yao X, Wu L and Zheng W (2011) Fault detection filter design for Markovian jump singular systems with intermittent measurements. IEEE Transactions on Signal Processing 59: 3099–3109. Yue D, Tian EG, Wang ZD and Lam J (2009) Stabilization of systems with probabilistic interval input delays and its applications to networked control systems. IEEE Transactions on systems, Man, and Cybernetics – Part A: Systems and Humans 39: 939–945. Zhang Y and Fang H (2011) Stabilization of nonlinear networked systems with sensor random packet dropout and time-varying delay. Applied Mathematical Modelling 35: 2253–2264. Zheng Y, Fang H and Wang HO (2006) Takagi–Sugeno fuzzy-modelbased fault detection for networked control systems with Markov delays. IEEE Transactions on Systems, Man and Cybernetics – Part B: Cybernetics 36: 924–929. Zhong M, Ye H, Shi P and Wang G (2005) Fault detection for Markovian jump systems. IEE Proceedings of Control Theory and Applications 152: 397–402.