Fault isolation filter for networked control system with event-triggered sampling scheme.

In this paper, the sensor data is transmitted only when the absolute value of difference between the current sensor value and the previously transmitted one is greater than the given threshold value. Based on this send-on-delta scheme which is one of the event-triggered sampling strategies, a modified fault isolation filter for a discrete-time networked control system with multiple faults is then implemented by a particular form of the Kalman filter. The proposed fault isolation filter improves the resource utilization with graceful fault estimation performance degradation. An illustrative example is given to show the efficiency of the proposed method.

Introduction

Over recent years, fault diagnosis for networked control system (NCS) using the mode-based analytical redundancy method have received significant attention. In [1,2], the overviews of main ideas and results on fault diagnosis of NCS are given, including the fundamentals of fault diagnosis for NCS with information scheduling, fault diagnosis approaches based on the simplified time-delayed system models and the quasi T-S fuzzy model, and fault diagnosis for linear and nonlinear NCS with long time-delay. However, most of the available results make use of time-triggered state estimation techniques by sampling the output of plant at an essentially equidistant time instant. Because the sampling period is determined according to the worst case operation conditions that rarely occur, the time-triggered sampling leads to a conservative usage of the communication bandwidth.

On the other hand, recent advances in computing and communication technologies enable the wireless networks (e.g., Bluetooth, wirelessHART and ZigBee) to rapidly replace wired networks in many applications, including industrial control and monitoring, home automation and consumer electronics, security and military sensing, and health monitoring [3,4]. Though the wireless channels are easier and cheaper to deploy and avoid cumbersome cabling, they also pose serious resource constraints. Therefore, applying the time-triggered sampling method to the wireless NCS may have some negative effects on the estimation and control performance of system, such as wasting the scarce communication resource and further shortening the lifetime of overall system. Since the event-triggered sampling strategies present a number of potential advantages for NCS, such as clock-free operation, less traffic requirement, and better resource utilization, they have been regarded as the possible and important alternatives to the time-triggered sampling.

Until now, numerous event-triggered sampling concepts have been proposed in the literature, such as send-on-delta sampling [5,6], level-crossing sampling [7], deadband sampling [8], Lebesgue sampling [9], send-on-area sampling [10], error energy sampling [11], self-triggered sampling [12], etc. Although these schemes have different terminologies, the same attribute is that the signal is sampled only when an a priori defined events occurs in the data monitored by sensors. For instance, the studies in [5-9] are concerned with the same sampling criterion where the event is defined as that the difference Δ between the current sensor value and the last transmitted one is greater than a given threshold. While in [10] and [11], the event is that the integral and energy of Δ is greater than a given threshold, respectively. Because of inherent distinctive benefits, the so-called event-triggered state estimation and event-based control for NCS with event-triggered sampling schemes have gained increasing attention. In this paper, however, we focus our attention on the event-triggered state estimation for the purpose of fault diagnosis. In relation to a parallel line of research on the event-based control, we refer the readers to the literature, e.g., [12,13].

In the context of state estimation, although the time-triggered state estimation over networks with network-induced effects taken into account have made great progress (see, e.g., [14-17]), research on the event-triggered state estimation is relatively lacking apart from several works [6,10,18-26]. It is well known that utilizing more sensors can potentially improve the performance of the estimation algorithms. However, using too many sensors can in turn create bottlenecks in the communication resource when these sensors compete for bandwidth. As a result, the studies in [18-21] explore the tradeoff between communication and estimation performance. Rather than sending every raw measurement to the remote estimator via network, a so-called controlled communication policy was adapted, which firstly obtain the local estimate x̃ from the raw sensor measurements and then compare x̃ with the remote estimate to decide whether or not it is worth sending data x̃. Also, Reference [21] proposes an optimal communication policy by dynamic programming and value iteration to minimize a long-term average cost function, which is related to the difference between the local and remote estimate. Based on the send-on-delta method, Reference [6] proposes a modified Kalman filter where computed output with increased measurement noise covariance is used when there is no sensor data transmission. The authors also discuss how to choose the threshold which is a trade-off parameter between the sensor data transmission rate and the estimation performance. Reference [22] extends the previous work [6] to address how to determine the measurement value at a sensor node if it does not send data. To avoid the inability of send-on-delta method in detecting the signal oscillations or steady-state error, Reference [10] proposes a novel scheme called send-on-area and then formulates a networked estimator based on Kalman filter to estimate the states of the system. More recently, Reference [23] proposes a networked estimator for event-triggered sampling systems with packet dropouts. Reference [24] develops an event-triggered estimator which is updated both when an event occurs with a received measurement sample, as well as at sampling instants synchronous in time without receiving a measurement sample. However, to the authors’ knowledge, fault diagnosis of networked control systems making use of the event-triggered state estimation method has not been addressed, which motivates the current study of this paper.

In this paper, we show our attention on the implementation problem of a modified fault isolation filter (FIF) for NCS with event-triggered sampling and multiple faults. By the send-on-delta scheme which is one of the event-triggered sampling strategies, it means that the sensor data is transmitted only when the absolute value of difference between the current sensor value and the previously transmitted one is greater than the given threshold value. Based on this scheme, a modified FIF for a discrete-time NCS with multiple faults is then implemented by a particular form of the Kalman filter. The rest of this paper is organized as follows. The modelling of NCS with event-triggered sampling scheme is presented in Section 2. A modified FIF is proposed in Section 3. An illustrative example is presented in Section 4 to show the effectiveness of the result. The paper is concluded in Section 5.

Notations: In what follows, if not explicitly stated, matrices are assumed to have compatible dimensions. Z+ denotes the set of nonnegative integer numbers. ℛ and ℛ are, respectively, the n-dimensional Euclidean space and the set of n × m real matrices. A denotes the transpose matrix or vector A. A−1 and A+ represent the inverse and pseudo-inverse of A, respectively. diag(a1,…,a) refers to an n × n diagonal matrix with a as its ith diagonal entry. rank(A) stands for the rank operator of matrix A. ℰ(x) represents the mathematical expectation of random variable x. x ∼ 𝒩(μ, Σ) means that the random vector satisfies the normal distribution with mean value μ and covariance matrix Σ. 𝒫(x|y) means the conditional probability distribution of x given y. Sign function is defined as

NCS with Event-Triggered Sampling Scheme

The architecture of NCS with event-sampling discussed in this paper is shown in Figure 1, where the closed-loop system consists of a plant with smart sensors and actuators, a remote FIF and a wireless network channel. The controller, FIF and actuator are assumed to be logically integrated. Thus, control commands do not need to experience any wireless transmission. This configuration represents a system, e.g., wireless sensor/actuator system, where actuation is inexpensive but sensor measurements are transmitted to the controller or FIF by sensors with a limited energy.

Figure 1.

The architecture of NCS with event-triggered sampling.

Since the event generator, the controller and the FIF have to be implemented on smart sensors and actuators by means of digital hardware, a discrete-time plant model is considered as the alterative to the continuous plant together with a zero-order hold and a sampler. The state evolution and sensor measurement equation are given as follows, respectively: where x ∈ ℛ is the state vector, u is the known manipulated input, is the fault distribution matrix, n ∈ ℛ is the fault vector, and the output observation vector of sensors. We assume that the initial state vector x0, process noise w and measurement noise v are uncorrelated, zero mean white Gaussian random processes with x0 ˜ 𝒩 (μ0, Σ0), w(0, W) v(0, R), where Σ0, W and R are symmetric, positive definite matrices.

As shown in Figure 1, the smart sensor numbered i has a sampler which regularly samples the sensor measurement with period h and an event generator to decide whether or not to send new sensor measurement through the network. The event generator therein, also known as send-on-delta scheme, is illustrated in Figure 2 where the sensor data is transmitted only when the absolute value of difference between the current sensor value and the previously transmitted one is greater than the given threshold value δ, namely where δ is simply determined according to x-fold of the maximal amplitude of .

Figure 2.

The architecture of event generator.

Obviously, all the transmitted measurements are the event-triggered samplers which are subsequences of the raw measurement . For instance, if the time instant when the previous event occurs is denoted as k ∈ Z+, the event-triggered sampling condition (2) is further formulated as

However, the time instant k (j = 1, 2....) can not be precisely determined because of the sampler, which is significantly different from the continuous one introduced in e.g., [5,6]. As shown in Figure 3, the dash line and real line represent the sensor measurement in the continuous time and discrete time, respectively. In the continuous time case, the time t when an event occurs can be exactly known since the measurement of plant is continuously updated by the sensor. While in the discrete time, the events have to be generated at the subsequent discrete time steps, e.g., k1 equivalent to kh, k2 equivalent to (k + 4)h, since the measurement is only updated at some discrete time instants and remains constant in the inter-sampling interval. Although producing the unsent measurements, e.g., the sample at the time instant (k + 1)h, wastes a bit of computing resource, the unsent measurements in turn save a lot of communication resources. In a sense of improving the resource utilization, applying the event-triggered sampling schemes to discrete-time plant is also meaningful. The reason is that wireless communication consumes more energy than information processing. As noted in [5], a sensor node can execute 3,000 instructions for the same energy cost of sending a single bit at the distance of 100 meters by radio.

Figure 3.

Send-on-delta sampling in discrete-time case.

In the sequel, we further assume that all sensor measurements are time-stamped and the network is communication link without packet losses and time delays. For the purpose of reducing communication, only parts of the raw measurements will be communicated to the remote FIF by the send-on-delta scheme. Through the communication link, all the event-triggered samplers are then stored in an infinite buffer. If sensor measurement does not yet arrive at the buffer at time instant k, it means that the current value has not significantly changed in contrast to the previous event-triggered sampler . In this case, the previous buffer value whose value is equivalent to will be stored in the k-slot of the buffer, as illustrated in Figure 4.

Figure 4.

The configuration of the buffer.

Furthermore, the arrival of the measurement at time k is defined as a binary variable , namely when arrives at the buffer at time instant k, otherwise . Rather than considering the design problem where only the probabilities of at each time instant is known and the research interest lies in studying the effect of loss and delay probabilities, we address the implementation problem where the value of at each time instant is known in advance. The last received value of i-th sensor output at time instant k is denoted as . If there is no sensor data received for k > k, the estimator node considers that the measurement value of the i-th sensor output is still equal to , but the measurement noise is increased from to where satisfies .

From the existing literature [6,22], the assumption that has a uniform distribution with zero mean and a variance is valid only if the measurement covariance R > δ2. Otherwise, the mean and variance of Δ(k, k) is and δ/12 respectively. Thus, the measurements which the FIF will use for fault diagnosis are formulated as and the output noise of the smart sensor numbered i is defined as: if . Otherwise, the output noise is defined as:

Moreover, the following selector is designed to flexibly determine whether (5) or (6) is applied to the FIF: where ε > 0 is a sufficiently small threshold.

By the selector (7) and compensating some values for the unsent measurements, we can also regularly implement some existing fault isolation filter algorithms, e.g., [27] to NCS with event-triggered sampling even though the measurements are transmitted irregularly.

Modified Fault Isolation Filter

Fault isolation filter, a special dynamic observer which generates the directional residuals in response to a particular fault, is an attractive way for enhancing the fault isolability [28,29]. It was first developed by Beard [28] and [29] and later revisited by Massoumnia [30] in the geometric framework and by White and Speyer [31] and Park and Rizzoni [32] in the context of eigenstructure assignment. Further improvements were suggested by Liu and Si [33], and Keller [27]. For linear continuous time-invariant system, Liu and Si [33] have proposed a fault isolation filter such that faults can be asymptotically detected and isolated. To guarantee that the ith component of the output residual is decoupled from all but the ith fault, the columns of the fault detectability matrix are assigned as the eigenvectors of the filter’s transition matrix with a set of fixed eigenvalues. Keller extended this approach to discrete-time stochastic linear systems. A new fault isolation filter has been developed to isolate q faults with at least q output measurements. More recently, Reference [34] addressed the fault detection problem for a class of linear networked control systems by extending the FIF proposed in [27].

In this section, we further construct a modified Keller’s FIF to detect and isolate the multiple faults in NCS with event-triggered sampling.

Recalling the definitions of fault detectability indexes and matrice introduced in [27]:

Definition 1. The linear stochastic system (1) is said to have fault detectability indexes ρ = {ρ1, ρ2, …, ρ} if ρ = min{ν: CAν−1 f ≠ 0, ν = 1,2,…}

Definition 2. If the linear stochastic system (1) has finite fault detectability indexes, the fault detectability matrix D is defined as: with

Now, the following filter is presented as the residual generator of discrete-time plant (1): where x̂ is the state of the filter, α the residual generator or the fault indicator. Filter gain K ∈ ℛ and projector L ∈ ℛ are unknown matrices to be found for the solution of the fault detection and isolation problem.

From Equations (1) and (9), the state estimation error e = x and the output of the filter α propagate as

Let G(z) be the transfer function from n to the output residual α. Then the following theorem is presented to design K and L such that which ensures the isolation of multiple faults.

Theorem 1. Under the condition rank(D) = q, the solutions of (11) can be parameterized as K = ωΠ+K̄Σ, L = Π, with Σ = β(I − DΠ), Π = D+, ω = AΨ and D = CΨ, where K̄ ∈ ℛ rank(Σ) = m − q

From Theorem 1, the fault isolation filter (9) is rewritten from the free parameter K̄ as where α is a deadbeat filter of fault n and given by: The fault of detectability index ρ directly affects the reduced output residual α with a time delay equals to its detectability index. ᾰ is the fault indicator signal without faults and propagates from the fault-free state estimation error ē = x̃ − x̂ as where Ā = A − ωΠC, C̄ = ΣC and x̃ is the fault-free state.

From Equation (14), the following theorem is then proposed to design the free parameter K̄ which minimizes the trace of the fault-free state estimation error covariance matrix .

Theorem 2. The proposed fault isolation filter described by the following relations: with V̄ = ΣRΣ, W̄ = W + ωΠRΠω.

Based on Theorem 2 and the measurement noise shown in Equations (5) and (6), the modified FIF for NCS with event-triggered sampling is proposed as follows:

set x0, P0, δ, ε

if ith sensor data are received

else if is a sufficiently small sclar

end if

Computing with V̄ = ΣR̄Σ, W̄ = W + ωΠR̄Πω.

Illustrative Example

In this section, we will present an example to illustrate the implementation approach proposed in this paper. The modified example is borrowed from [27] described by Equation (1), where the parameters are as follows: The faults associated with the first and second column of the fault distribution matrix F occur at time instant r1 = 10 with and r2 = 30 with , respectively.

By Theorem 1, we have D = CΨ and Ψ = F since rank(CF) = q (ρ1 = 1, ρ2 = 1). The parametrization of the FIF’s gain K and projector L is then given by

Furthermore, we choose -fold of the maximal amplitude of (i = 1,2,3). Then, the threshold values in Equation (3) are determined as δ1 = 4.8406, δ2 = 2.0382, δ3 = 1.0009. The sufficiently small threshold ε (i = 1,2,3) in Equation (7) is given as 10−4. Then, Figure 5(a)–7(a) show all the transmitted measurements of sensors with event-triggered scheme. Figures 5(b)–7(b) indicate all the transmitted measurements of sensors with time-triggered scheme. By comparison, the sensors with event-triggered scheme transmit only 66%, 67%, and 58% of samples produced by time-triggered scheme, respectively. In other words, the resource utilization by the event-triggered scheme can be obtained 34%, 33% and 42% improvement, respectively.

Figure 5.

Transmitted measurements of sensor by (a) event-triggered and (b) time-triggered.

Figure 7.

Transmitted measurements of sensor by (a) event-triggered and (b) time-triggered.

The measurements used in the simulation are indicated in Figure 8, where (i = 1,2,3) and (i = 1,2,3) represent the measurements with compensating and without compensating some values for the unsent ones, respectively. By Algorithm 1 and the measurements , a modified FIF is implemented to the NCS with event-triggered sampling. Figure 9 shows the innovation sequences of residual , and , which correspond to the residual obtained by the modified FIF described in Algorithm 1, the Keller’s FIF in [27] using the time-triggered samples and the Keller’s FIF using the measurements , (i = 1,2,3), respectively.

Figure 8.

Measurements (a) with compensation and (b) without compensation.

Algorithm 1.

Initialization

set x0, P0, δi, εi if ith sensor data are received (θki=1) R¯(i,i)=R(i,i),zki=yki,ηki=0 else if |yk_ji−yk_j−1i|≤ε_i% ε_i is a sufficiently small sclar R¯(i,i)=R(i,i)+δi2/3,ηki=0 else R¯(i,i)=R(i,i)+δi2/12,ηki=sign (yk_ji−yk_j−1i)*δi2/2 end if z_k=[zk1…zki…zkm]^Ty_k=[yk1…yki…ykm]^Tη_k=[ηk1…ηki…ηkm]^Tz_k=y_k+η_k Computing (16){x^_k+1=Ax^_k+Bu_k+ωα_k+K¯_kγ_kα_k=Π(z_k−Cx^_k)γ_k=Σ(z_k−Cx^_k)P¯_k+1=(A¯−K¯_kC¯)P¯_k(A¯−K¯_kC¯)^T+K¯_kV¯K¯kT+W¯+ωΠR¯Σ^TK¯kT+K¯_kΣR¯Π^Tω^TK¯_k=(A¯P¯_kC¯^T−ωΠR¯Σ^T)(C¯P¯_kC¯^T+V¯)^−1with V̄ = ΣR̄ΣT, W̄ = W + ωΠR̄ΠTωT.

Figure 9.

Innovation sequences of residual α, and .

In order to compare the simulation results obtained by different schemes, we propose the root mean square (RMS) of the fault estimation errors as a performance index. This error for the scalar variable n with respect to its estimate α for N simulation steps is defined as where is the value of the variable n in the jth step, is the estimate of and i = 1,2. Table 1 indicates the RMS of the fault estimation errors by different methods.

Table 1.

RMS of the fault estimation errors.

α	(α1)k	(α1#)k	(α1★)k	(α2)k	(α2#)k	(α2★)k
RMS	1.6258	0.9368	18.0693	1.7998	1.6235	14.2398

From Figure 9 and Table 1, the performance of the modified FIF using less sensor data transmission experiences graceful degradation in contrast to the Keller’s FIF using time-triggered samples. By graceful degradation, it means that a system degenerates in such a manner that it continues to operate, but provides a reduced level of service rather than causing total breakdown. On the other hand, the Keller’s FIF using the measurements (i = 1,2,3) fails to work completely, while the number of the total sensor data transmission is same as the modified FIF used. From these results, we can draw a conclusion that the modified FIF has an advantage in the tradeoff between communication cost and fault estimation performance.

In order to establish the relationships between the number of the sensor data transmissions, the RMS of the fault estimation errors by different methods and the threshold values δ (i = 1,2,3), some further simulations are done by adjusting the threshold values δ (i = 1,2,3) from 1/20 to 1/80 of the maximal amplitude of (i = 1,2,3). Figure 10 shows the percentage of transmitted samples to total samples in relation to δ by event-triggered scheme, wherein the real line, dash line and dash-dotted line represent the percentage for sensor , and , respectively. From Figure 10, it can be seen that the sensor data transmission rate is inversely proportional to δ, namely the communication cost reduced by the event-triggered scheme will increase when δ increases, and vice versa.

Figure 10.

Percentage of transmitted samples to total samples in relation to δ by event-triggered scheme.

In Figure 11, the real line, dash line and dash-dotted line represent the fault estimation performance of the modified FIF, the Keller’s FIF with time-triggered sampling and the FIF with the measurements , respectively. It can be seen that the fault estimation performance of the modified FIF is improved and eventually approaches the performance of Keller’s FIF with time-triggered sampling as δ decreases. However, the performance of the FIF with the measurement (i = 1,2,3) is very poor even by flexibly adjusting the threshold δ.

Figure 11.

RMS of fault detection errors in relation to δ.

Conclusions

This paper is concerned with the implementation problem of fault isolation filter for networked control system with the send-on-delta scheme which is one of the event-triggered sampling strategies. By send-on-delta, the sensor data is transmitted only when the absolute value of difference between the current sensor value and the previously transmitted one is greater than the given threshold value. Based on this scheme, a modified fault isolation filter for a discrete-time networked control system with multiple faults is then implemented by a particular form of the Kalman filter. In contrast to the Keller’s fault isolation filter using time-triggered samples, the proposed fault isolation filter improves the resource utilization with graceful fault estimation performance degradation. Also, we can improve the performance of the modified FIF by flexibly adjusting the threshold values with taking resource utilization into account.

Throughout the paper, no network-induced packet losses are taken into account in the model of the networked control system. Further study of the fault isolation filter for networked control system with event-triggered sampling and packet losses is encouraged.