Networked Estimation for Event-Based Sampling Systems with Packet Dropouts

This paper is concerned with a networked estimation problem in which sensor data are transmitted over the network. In the event-based sampling scheme known as level-crossing or send-on-delta (SOD), sensor data are transmitted to the estimator node if the difference between the current sensor value and the last transmitted one is greater than a given threshold. Event-based sampling has been shown to be more efficient than the time-triggered one in some situations, especially in network bandwidth improvement. However, it cannot detect packet dropout situations because data transmission and reception do not use a periodical time-stamp mechanism as found in time-triggered sampling systems. Motivated by this issue, we propose a modified event-based sampling scheme called modified SOD in which sensor data are sent when either the change of sensor output exceeds a given threshold or the time elapses more than a given interval. Through simulation results, we show that the proposed modified SOD sampling significantly improves estimation performance when packet dropouts happen.


Introduction
Recent works have discussed event-driven alternatives to traditional time-triggered sampling schemes. It has been shown to be more efficient than time-triggered one in some situations, especially in network bandwidth improvement. In [1][2][3][4][5][6][7], event-based sampling scheme was applied by adjusting OPEN ACCESS the threshold value at each sensor node, data transmission rate is reduced so that the network can be used for other traffic.
However, analysis and simulation in the the works on event-driven sampling scheme were performed under ideal communication network conditions: no delays or packet dropouts are assumed, but in realistic applications, network induced delays and packet losses do happen.
The issues of network delays and packet dropouts in time-triggered systems have been addressed and solved by researchers in [8][9][10][11][12][13][14]. In [8] the stability of the Kalman filter in relation to the data arrival rate is investigated. It is shown that there exists a critical data arrival rate for an unstable system so that the mean filtering error covariance will be bounded for any initial condition. In a very recent study [13], the optimal H 2 filtering problems associated respectively with possible delay of one sampling period, uncertain observations and multiple packet dropouts are studied under a unified framework. The H 2 -norm of systems with stochastic parameters is defined and computed via a Lyapunov equation and a steady-state filter is designed via an LMI approach. In [14], the authors adopt a model similar to that of [13] for multiple packet dropouts to investigate finite-horizon optimal linear filtering, prediction and smoothing problems.
In conventional event-based sampling systems, also called send-on-delta (SOD) sampling [5][6][7], the issues of network delay and packet loss are difficult to solve because data transmission and reception do not use a periodical time-stamp mechanism as in the time-triggered sampling systems. Motivated by those issues, in this paper, we introduce a modified SOD sampling scheme in which the event-driven sampling is combined with a time-triggered sampling scheme to detect packet dropouts. Then, a networked estimator based on a Kalman filter is formulated to estimate states of the system periodically even when the sensor nodes do not transmit data. The proposed SOD sampling scheme has properties inherited from the conventional SOD sampling: so the benefits from event-driven sampling are still hold. Through theoretical analysis and simulation results, we show that the proposed SOD sampling scheme gives better estimation performance than the conventional SOD one when packet loss happens.

Modified SOD Sampling Scheme
Consider a networked control system described by the linear continuous-time model: is the state of the plant, u is the deterministic input signal, ( ) is the measurement output which is sent to the estimator node by the sensor nodes. ( ) w t is the process noise with covariance Q , and ( ) v t is the measurement noise with covariance R . We assume that ( ) w t and ( ) v t are uncorrelated, zero mean white Gaussian random processes.
The modified SOD sampling scheme illustrated in Figure 1b is stated as follows: be the last transmitted value of the i -th sensor output at instant , last i t . A new sensor value will be sent to the estimator node if one of two following conditions is satisfied:  Using the modified SOD sampling scheme above we will obtain some benefits. Firstly, the estimator can detect signal oscillations or steady-state error if the difference of output value remains within the threshold range during a long time. Secondly, the estimator can detect multiple packet dropouts if it does not receive sensor data within the interval , (0, ) t i d . Thirdly, theoretical analysis for SOD sampling is still applied for the modified SOD sampling.
However, this scheme has one disadvantage that sensor data transmission rate will be increased due to condition (2b). If , t i d is small, the estimator detects packet dropouts fast but data transmission rate is increased. If , t i d is large, transmission rate is small but the estimator detects packet dropouts slowly. Therefore, an optimal , t i d value is necessary to compromise these constraints.

Multiple packet dropouts detection
The estimator node detects packet dropouts of i-th sensor data by checking the instant i-th sensor data arrive. If there is no i-th sensor data arriving, the estimator node for the time We state the general case for multiple packet dropouts as follows: If the estimator node does not receive i-th sensor data for time ( ) . This is an inevitable flaw of the modified SOD sampling scheme. We can constraint this flaw by reducing the , t i d value, but sensor data transmission rate will be increased. Therefore, an optimal , t i d value is necessary to compromise between the two constraints.

State Estimation with Modified SOD Transmission Method
The networked estimation problem applying modified SOD transmission method can be described as follows: are sampled at the period T but their data are only sent to the estimator node when (2a) or (2b) is satisfied. 2. For simplicity in the problem formulation, transmission delay from the sensor nodes to the estimator node is ignored. 3. The estimator node estimates states of the plant regularly at the period T regardless of whether or not sensor data arrive. If there is no i -th sensor data received for ( ) , the estimator node considers that the measurement value of the i -th sensor output ( ) i y t is still equal to , last i y but the measurement noise increases from ( ) Note that if 0 i d = then there is no packet dropout, the estimator acts like a conventional SOD filter [5]. To formulate a state estimation problem, the boundry of needs to be determined as 0 i d ¹ (packet dropouts happen). In the next section, we will compute the covariance of , ( ) Note that (3) is also applied to the case of no packet dropout [5] Therefore, if there is no i -th sensor data received for

State estimation
A modified Kalman filter for state estimation ˆk x at step k, where there is a change in the measurement update part of the discrete Kalman filter algorithm [15], is given as in the Figure 4. We use the discretized system model sampled at period T : In the modified Kalman filter in Figure 4, the states of the plant are estimated regularly at every period T, regardless of whether or not sensor data arrive. If i -th sensor data arrive then as measurement noise covariance for state estimation. As stated in [8], if the system (1) is unstable and a packet loss rate is high, the proposed filter could diverge. For example, if all packets are lost, d i will increase and thus i R will become infinite. Thus P in Figure 4 could become infinite.

Optimal δ t,i Computing Problem
As mentioned in Section 3, δ t,i is a trade-off parameter between sensor data transmission rate and the response of packet dropouts detection. The response of packet dropout detection guarantees estimation performance. Because SOD sampling is more efficient than the time-triggered one in network bandwidth improvement, we should choose δ t,i such that sensor data transmission rate is reduced to promote ability of SOD sampling. In the next section, we will investigate the relation of δ t,i with transmission rate and the effect of δ t,i on estimation performance. Then an optimization problem is formulated to find the optimal δ t,i value according to the given estimation performance.

Sensor data transmission rate by condition (2b)
The total sensor data transmission rate caused by condition (2b) in a time unit: where p is the number of sensor output The average number of packet dropouts in the conventional SOD sampling per a time unit:

Estimation error covariance due to packet dropouts
In the proposed SOD sampling, the average number of packet dropouts within the time interval , t i d will be: We know from Section 4.1 that the larger number of consecutive packet dropouts is, the larger measurement noise covariance is. Measurement noise covariance is largest if i d packets are consecutively lost. Following the idea in (4), if there is i d packet loss, the measurement noise covariance should be increased as follows: ( ) The estimation performance in this case can be computed from the following discrete algebraic Riccati equation: Note that (10) does not provide the actual estimation error covariance of the filter. The main purpose of (10) is to evaluate how δ t,i affects the estimation performance. We can see that if δ t,i is large, the estimation error covariance P increases.
The solution of (10) is denoted by ( ) , t i P d . In the following optimization algorithm to find δ i , we try to reduce the sensor transmission rate caused by condition (2b) subject to the given estimation performance constraint: ). 0 P is also the estimation performance of the conventional SOD. m is the ratio to the estimation performance of conventional SOD filter in case of no packet dropout. If m is large, the , t i d optimization problem (11) is done with weaker estimation performance constraints.

Simulation
To verify the proposed filter, we consider an example of the second-order system with step input where the output is sampled by the SOD and modified SOD sampling: where the system parameters for performance evaluation are given by 30, a 5, b 1 M = = = (underdamped system) . The simulation process is implemented for 50 seconds. Choose 5 m = for the optimization problem (11 Table 1 shows the estimation error in two filters (SOD filter and modified SOD filter) as where i x is the reference state, ˆi x is the estimated state, and N = 5,000.
In Table 1, we see that when applying the modified SOD filter, the estimation error is significantly improved. For instance, in the case 1 2 0.05 x x = = , the total number of sensor data transmissions in the modified SOD (# 137) is just slightly greater than that in conventional SOD (# 126) but the estimation error is reduced so much ((e 1 = 0.0075, e 2 = 0.0096) compared to (e 1 = 0.0383, e 2 = 0.0167)).   The boundry of 1 e in the modified SOD filter (SODa) is much smaller than that in the conventional SOD filter. Figure 8 shows the instants the sensor node transmits data to the estimator node due to condition (2b). We see that the number of sensor data transmissions caused by condition (2b) is very small in comparison with the total number of sensor data transmissions [(n 1 = 7, n 2 = 7) compared to (n 1 = 101, n 2 = 36)]. When the modified SOD sampling is applied, the total number of sensor data transmissions is slightly increased, but the estimation error is significantly reduced. Therefore, the modified SOD sampling significantly improves estimation performance with only a little increase in the data transmission rate. Notice that if we just consider the transmission condition (Equation 2a), estimation error of the proposed method is worse for systems that the output varies slowly. However, an issue of conventional event-based sampling is that it can not detect signal oscillations or steady-state error if the difference of output value remains within the threshold range (because the output varies slowly). This fact causes estimation error to be increased. Whereas, the proposed method uses the transmission condition (Equation 2b) not only to detect packet dropouts but to reduce the error in case the output changes slowly.
As illustrated in Figures 7 and 8, where the estimation error of the proposed method (top-right graph of Figure 7) and of the conventional method (top-left graph of Figure 7) are shown according to the output y1 (top-left graph of Figure 8). We see when y1 varies slowly (time interval from 20s to 50s), the proposed method gives much smaller estimation error than the conventional one.
In case the output changes fast, it is obvious that ignoring packet dropout will introduce extremely incorrect result because we still use the wrong old measurement noise value even when we do not know how much the output value changes.

Conclusions
In this paper, the state estimation problem with modified SOD transmission method over networks, in which an event-based sampling is combined with a time-triggered sampling to detect packet loss situations, has been considered. We have shown that when using the proposed modified SOD filter, estimation performance is significantly improved with a small increase in sensor data transmission. If multiple packet dropouts happen, the estimator node will detect and compensate for them with an amount of additive measurement noise to improve estimation performance. This method is very useful for networks where data transmission is unreliable due to noise.