Next Article in Journal
Timepix3: Temperature Influence on Radiation Energy Measurement with Si Sensor
Previous Article in Journal
Adversarial Auxiliary Weighted Subdomain Adaptation for Open-Set Deep Transfer Bridge Damage Diagnosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Event-Triggered Kalman Filter and Its Performance Analysis

1
School of Electronic Engineering, Heilongjiang University, Harbin 150080, China
2
Key Laboratory of Information Fusion Estimation and Detection, Harbin 150080, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(4), 2202; https://doi.org/10.3390/s23042202
Submission received: 16 January 2023 / Revised: 12 February 2023 / Accepted: 13 February 2023 / Published: 15 February 2023
(This article belongs to the Section Physical Sensors)

Abstract

:
In estimation of linear systems, an efficient event-triggered Kalman filter algorithm is proposed. Based on the hypothesis test of Gaussian distribution, the significance of the event-triggered threshold is given. Based on the threshold, the actual trigger frequency of the estimated system can be accurately set. Combining the threshold and the proposed event-triggered mechanism, an event-triggered Kalman filter is proposed and the approximate estimation accuracy can also be calculated. Whether it is a steady system or a time-varying system, the proposed algorithm can reasonably set the threshold according to the required accuracy in advance. The proposed event-triggered estimator not only effectively reduces the communication cost, but also has high accuracy. Finally, simulation examples verify the correctness and effectiveness of the proposed algorithm.

1. Introduction

The growing demands of communications, navigation positioning, fault detection, and environmental monitoring have spawned wireless sensor network systems [1,2,3,4,5]. In wireless sensor networks, due to the limitation of bandwidth and energy, frequent data transmission is unfavorable, which not only increases the communication costs, but also induces some negative network phenomena, such as packet loss, network delay, etc. In order to reduce the waste of resources in the operation of the equipment and ensure the quality of communication, a better transmission mechanism is worth seeking. Unlike traditional time-triggered mechanisms, an event-triggered one transmits information only when needed. It effectively alleviates the problems of network congestion and waste of resources, and has attracted widespread attention from scholars [6,7,8,9,10].
As early as 1983, Ho et al. first introduced the concept of event trigger into discrete systems [11]. Later, in the development of event-triggered filter estimation, a trigger mechanism named send-on-delta (SOD) emerges. The trigger mechanism is based on differences in the measured signals, and when the signal deviates from the delta, the communication is triggered. It is one of the classic solutions to solve redundant data transmission in wireless sensor networks [12]. Subsequently, in [13], using estimation variance as a trigger condition, Trimpe et al. correlated trigger decision with estimator performance, and estimation algorithms for distributed sensor nodes that can be applied in two cases were designed, and another popular trigger mechanism was proposed. Based on the innovation and its covariance information, Wang N.et al. proposed an event-triggered sequence fusion estimator. It is an improvement of the SOD mechanism, which solves the problem of event-triggered estimators of observation noise under the linear minimum variance criterion [14]. Zhang et al. introduced the event trigger mechanism into the  H  filtering framework. The trigger conditions depend on the current measurement and the difference between the latest transmission, which significantly saves communication resources, and the filtering error has proved to be asymptotically stable [15]. In [16], the event-triggered mechanism under malicious denial-of-service attacks is studied, and estimating performance and network communication are analyzed.
The trigger mechanisms used in the above works are all static; that is, a constant is used as the trigger threshold [17]. With the in-depth study of event triggers, the setting of the threshold is no longer conservative and dynamic event trigger mechanisms appeared [18]. In [19], taking into account the unknown input of the network control system, an adaptive threshold including system state vector and neural network weight estimation is designed, and the stability of the system is rigorously proved by using the Lyapunov stability theory. In [20], for the multi-agent system, event-triggered mechanisms in the centralized and distributed cases are developed, and the dynamic threshold of the distributed system is set by introducing internal variables and exponential functions, which greatly reduces the transmission burden. For uncertain systems with transmission delays, a distributed transmission scheme based on a dynamic event trigger is proposed in [21]. Different from [20], the trigger time series of the system is determined by the dynamic event variable rules. In [22], a novel decentralized trigger condition is designed, whose trigger threshold is constructed by the rate of change, and the Zeno behavior in the trigger process is avoided by adding a corresponding constant term.
In fact, the selection of the event-triggered threshold depends on the preset accuracy and trigger rate, but there are few works on the accuracy range and trigger rate at a certain threshold. Compared with the traditional trigger mechanism, the one proposed in this paper can give a quantitative relationship between trigger threshold, communication rate, and estimation error. Clearly, it is important for wireless sensors to clarify the relationship between the communication frequency and the corresponding estimation accuracy range. At the same time, the trigger mechanism with better performance at the same trigger rate is sought. For these reasons, the motivation for our study came into being. The main contributions of this paper are summarized as follows:
  • For wireless sensors, an event-triggered mechanism is proposed. The proposed event-triggered mechanism is based on a normal distribution constructed from the predicted and filtered differences. When the constructed normal distribution trigger function exceeds a tolerable threshold, the filtering results are transmitted.
  • The theoretical trigger probability and estimation accuracy under the proposed trigger mechanism are derived. For linear time-invariant systems, the Riccati equation can be used to obtain the estimation error variance of the event-triggered estimator. For linear time-varying systems without steady-state estimation, an approximate estimation accuracy can be obtained. Therefore, their trigger threshold can be set according to the precision requirements.
  • The simulation verifies the correctness of the proposed theorems and inferences. Compared with several types of event trigger mechanisms in the existing literature [14,22], the trigger mechanism proposed in this paper has higher estimation accuracy and better performance under the same trigger rate.
Notation:  R n  denotes the  n -dimensional Euclidean space,  R n × m  denotes the set of  n  by  m  real-valued matrices, ‘E’ denotes the mathematical expectation,  M T  and  M 1  are the transpose and inverse of matrix  M , respectively,  chol ( P )  is the Cholesky decomposition of matrix  P tr ( M )  is the trace of matrix  M I n  is the  n -dimensional identity matrix,  N ( μ , δ 2 )  denotes the Gaussian distribution with mean  μ  and variance  δ 2 δ k t  is the Kronecker delta function  ( δ tt = 1 , δ t k = 0   ( t k ) ) x ˜ k | k τ = x k x ^ k | k τ  is the state estimation error, and  z ˜ k | k τ = z k z ^ k | k τ ( τ = 0 , 1 , , k 1 )  is the measurement prediction error.

2. Problem Formulation

Consider the following linear system:
x k + 1 = Φ x k + Γ w k
z k = H x k + v k
where  k  is discrete time,  x k + 1 R n  is the state vector,  z k R m  is the measurement vector,  Φ R n × n  is the state transition matrix,  Γ R n × m  is the noise transfer matrix,  H R m × n  is the measurement matrix,  w k  is the process noise, and  v k  is the measurement noise.
Assumption 1.
w k  and  v k  are uncorrelated white Gaussian noise processes with zero means and variances  Q k  and  R k , respectively,
E [ ( w k v k ) ( w t T   v t T ) ] = [ Q k   0 0 R k ] δ k t
Assumption 2.
The initial state  x 0  is dependent on  w k  and  v k , and satisfies:
E [ x 0 ] = μ 0 ,   E [ ( x 0 μ 0 ) ( x 0 μ 0 ) T ] = P 0

2.1. Threshold Selection of Event-Triggered Mechanism

Taking into account the characteristics of intelligence and the low power consumption of current sensors, as well as the advantages of event-triggered mechanisms based on state estimation, an event-triggered Kalman filter based on the state estimation is designed. The smart sensor node is responsible for continuous filtering and judgment. If the trigger condition is not satisfied, the prediction actions will be performed by the center; otherwise, the filter estimation will be sent by the smart sensor node and the target states will be updated by the receiving center. This structure can effectively reduce communication redundancy.
Theorem 1.
For the system in Equations (1) and (2), the test statistics obey the following distribution:
γ k = ( chol ( P k , τ x ) ) 1 ( x ^ k | k τ x ^ k | k ) ~ N ( 0 , I n )  
where
x ^ k | k = x ^ k | k 1 + K k z ˜ k | k 1
K k = P k | k 1 H T [ H P k | k 1 H T + R k ] 1
z ˜ k | k 1 = z k H x ^ k | k 1
x ^ k | k 1 = Φ x ^ k 1 | k 1
P k | k 1 = Φ P k 1 | k 1 Φ T + Γ Q k 1 ( Γ ) T
x ^ k | k τ  is the   τ s t e p  Kalman prediction:
x ^ k | k τ = Φ τ 1 x ^ k τ + 1 | k τ , τ 1  
and the second moment  P k , τ x  of   ( x ^ k | k τ x ^ k | k )  can be written as:
P k , τ x = E { ( x ^ k | k τ x ^ k | k ) ( x ^ k | k τ x ^ k | k ) T }               = P k | k P k | k τ P k | k τ ( 1 ) P k | k τ ( 2 ) P k | k τ ( τ ) ( P k | k τ ( 1 ) ) T ( P k | k τ ( 2 ) ) T ( P k | k τ ( τ ) ) T
where the filtering error variance   P k | k is:
P k | k = [ I n K k H ] P k | k 1  
the prediction error variance   P k | k τ  in Equation (12) is calculated by:
P k | k τ = Φ τ 1 P k τ + 1 | k τ ( Φ ( τ 1 ) ) T + j = 2 τ Φ τ j Γ Q k τ + j 1 Γ T Φ ( τ j ) T , τ 2  
  P k | k τ ( m ) ( m = 1 , , τ )  in Equation (12) is calculated by:
P k | k τ ( m ) = Φ τ 1 P k τ + 1 | k τ ( A m ) T H T K k τ + m T ( Φ τ m ) T + B  
and the solution steps for   P k | k τ ( m ) ( m = 1 , , τ )  in Algorithm 1 is:
Algorithm 1: Solution steps for  P k | k τ ( m )
Initialize:   A 1 = I n , C 0 = I n .
Iterate:
for   m = 1 : τ
A m + 1 = Φ A m Φ K k τ + m H A m
end
if  τ = 1 , B = 0  
  else  τ = 2 , m = 2 , B = Φ τ m Γ Q k τ + m + 1 Γ T I n T H T K k τ + m T ( Φ τ m ) T
  else  B m = 0
for  m = 3 : τ  
  h = m 2 ,   p = m 2  
for  ξ = 1 : m 2  
for  = 0 : h 1  
C + 1 = Φ C Φ K k τ + m h + H C  
end
  B m ψ = Φ τ m + p Γ Q k τ + m p 1 Γ T C h T H T K k τ + m T ( Φ τ m ) T  
h = h 1 ,   p = p 1  
end
B m = B m + B m ψ  
end
B = B m Φ τ m Γ Q k τ + m + 1 Γ T I n T H T K k τ + m T ( Φ τ m ) T  
end
                      P k | k τ ( m ) = Φ τ 1 P k τ + 1 | k τ ( A m ) T H T K k τ + m T ( Φ τ m ) T + B
Proof. 
Due to  x ^ k | k τ L ( z 1 , , z k τ )  ( L  is the linear manifold of  z 1 , , z k τ ) and  x ^ k | k L ( z 1 , , z k )  in Equation (5),  γ k L ( z 1 , , z k )  obeys the Gaussian distribution.
(1)
The proof of the mean value of  γ k . Since  x ^ k | k τ  and  x ^ k | k  are unbiased,
E { γ k } = E { ( chol ( P k , τ x ) ) 1 ( x ^ k | k τ x ^ k | k ) } = E { ( chol ( P k , τ x ) ) 1 } ( E { x ^ k | k τ } E { x ^ k | k } ) = 0
(2)
The proof of the  P k , τ x  in Equation (12). Based on the Kalman filter:
x ^ k | k = x ^ k | k 1 + K k z ˜ k | k 1 = Φ x ^ k 1 | k 1 + K k z ˜ k | k 1 = Φ 2 x ^ k 2 | k 2 + Φ K k 1 z ˜ k 1 | k 2 + K k z ˜ k | k 1 = Φ τ   x ^ k τ | k τ + Φ τ 1 K k τ + 1 z ˜ k τ + 1 | k τ + Φ τ 2 K k τ + 2 z ˜ k τ + 2 | k τ + 1 + + K k z ˜ k | k 1 = x ^ k | k τ + i = 1 τ Φ τ i K k τ + i z ˜ k τ + i | k τ + i 1
estimated error variance matrix  P k , τ x  is computed by:
P k , τ x = E { ( x ^ k | k τ x ^ k | k ) ( x ^ k | k τ x ^ k | k ) T } = E { [ ( x ^ k | k τ x k ) + ( x k x ^ k | k ) ] [ ( x ^ k | k τ x k ) + ( x k x ^ k | k ) ] T } = P k | k τ + P k | k E { ( x ^ k | k τ x k ) ( x ^ k | k x k ) T } E { ( x ^ k | k x k ) ( x ^ k | k τ x k ) T }
From Equations (17) and (18),  E { ( x ^ k | k τ x k ) ( x ^ k | k x k ) T }  can be written as:
E { ( x ^ k | k τ x k ) ( x ^ k | k x k ) T } = E { ( x ^ k | k τ x k ) ( x ^ k | k τ x k + i = 1 τ Φ τ i K k τ + i z ˜ k τ + i | k τ + i 1 ) T } = P k | k τ + E { i = 1 τ ( x ^ k | k τ x k ) [ Φ τ i K k τ + i ( H x ˜ k τ + i | k τ + i 1 + v k τ + i ) ] T } + E { ( x ^ k | k τ x k ) [ Φ τ 2 K k τ + 2 ( H x ˜ k τ + 2 | k τ + 1 + v k τ + 2 ) ] T } + + E { ( x ^ k | k τ x k ) [ K k ( H x ˜ k | k 1 + v k ) ] T } = P k | k τ + P k | k τ ( 1 ) + P k | k τ ( 2 ) + + P k | k τ ( τ )
since  w k i  is white noise and independent of  v k τ + i P k | k τ ( 1 ) P k | k τ ( 2 ) P k | k τ ( 3 )  in Equation (19) are respectively calculated as:
P k | k τ ( 1 ) = E { ( x ^ k | k τ x k ) [ Φ τ 1 K k τ + 1 ( H x ˜ k τ + 1 | k τ + v k τ + 1 ) ] T } = E { [ Φ τ 1 ( x ^ k τ + 1 | k τ x k τ + 1 ) L ( w k τ + 1 w k 1 ) ] [ Φ τ 1 K k τ + 1 ( H x ˜ k τ + 1 | k τ + v k τ + 1 ) ] T } = Φ τ 1 P k τ + 1 | k τ H T K k τ + 1 T ( Φ τ 1 ) T
P k | k τ ( 2 ) = E { ( x ^ k | k τ x k ) [ Φ τ 2 K k τ + 2 ( H x ˜ k τ + 2 | k τ + 1 + v k τ + 2 ) ] T } = E { [ Φ τ 1 ( x ^ k τ + 1 | k τ x k τ + 1 ) L ( w k τ + 1 w k 1 ) ] [ Φ τ 2 K k τ + 2 ( H x ˜ k τ + 2 | k τ + 1 + v k τ + 2 ) ] T } = Φ τ 1 P k τ + 1 | k τ ( I n K k τ + 1 H ) T Φ T H T K k τ + 2 T ( Φ τ 2 ) T Φ τ 2 Γ Q k τ + 1 Γ T H T K k τ + 2 T ( Φ τ 2 ) T
P k | k τ ( 3 ) = E { ( x ^ k | k τ x k ) [ Φ τ 3 K k τ + 3 ( H x ˜ k τ + 3 | k τ + 2 + v k τ + 3 ) ] T } = E { [ Φ τ 1 ( x ^ k τ + 1 | k τ x k τ + 1 ) L ( w k τ + 1 w k 1 ) ] [ Φ τ 3 K k τ + 3 ( H x ˜ k τ + 3 | k τ + 2 + v k τ + 3 ) ] T } = Φ τ 1 P k τ + 1 | k τ [ Φ 2 ( I n K k τ + 1 H ) Φ K k τ + 2 H Φ ( I n K k τ + 1 H ) ] T H T K k τ + 3 T ( Φ τ 3 ) T Φ τ 2 Γ Q k τ + 1 Γ T ( I n K k τ + 2 H ) T Φ T H T K k τ + 3 T ( Φ τ 3 ) T Φ τ 3 Γ Q k τ + 2 Γ T H T K k τ + 3 T ( Φ τ 3 ) T = Φ τ 1 P k τ + 1 | k τ ( A 3 ) T H T K k τ + 3 T ( Φ τ 3 ) T Φ τ 2 Γ Q k τ + 1 Γ T ( I n K k τ + 2 H ) T Φ T H T K k τ + 3 T ( Φ τ 3 ) T Φ τ 3 Γ Q k τ + 2 Γ T H T K k τ + 3 T ( Φ τ 3 ) T
From Equations (20)–(22), the solution rules can be summarized as Equation (15). Similarly,  E { ( x ^ k | k x k ) ( x ^ k | k τ x k ) T }  is defined by:
E { ( x ^ k | k x k ) ( x ^ k | k τ x k ) T } = P k | k τ + ( P k | k τ ( 1 ) ) T + ( P k | k τ ( 2 ) ) T + + ( P k | k τ ( τ ) ) T
In summary, the error variance matrix in Equation (18) can be calculated as Equation (12).
(3)
The proof of  P k | k τ  in Equation (14), from Equation (1):
x k = Φ τ 1 x k τ + 1 + i = 2 τ Φ τ i Γ w k τ + i 1 , τ 2
since  w k  is white noise and independent of  v k , the Equation (11) is yielded. Based on  P k | k τ = E { ( x k x ^ k | k τ ) ( x k x ^ k | k τ ) T } , Equation (11) and Equation (24), the Equation (14) is yielded. The remaining equations are easily obtained by the Kalman filter and prediction. □
Considering that the test statistic  γ k  obeys the standard Gaussian distribution, under the selected level of significance  α , the trigger threshold will be given by using the hypothesis test that satisfies the Gaussian distribution [23,24].
Corollary 1.
Under the selected significance level  α , set original hypothesis  H 0 x ^ k | k τ x ^ k | k , and antithetic hypothesis  H 1 x ^ k | k τ = x ^ k | k γ k , i  is the  i th  element of  γ k , if the observed value satisfies
| γ k , i | Z α / 2 , i = 1 , 2 , , n  
then reject  H 0  (non-trigger communication), otherwise accept   H 0  (trigger communication).  Z α / 2  is the  α / 2  quantile of the standard Gaussian distribution.
Remark1: 
Selecting the original hypothesis  H 0 : x ^ k | k τ x ^ k | k  as the trigger event, the purpose is to protect the trigger event  H 0 ; thus, priority is given to ensuring the estimation accuracy of the receiving center. On the contrary, if the original hypothesis is the event  H 1 x ^ k | k τ = x ^ k | k , this is to protect the event  H 1 , so as to give priority to saving the communication of smart nodes. Usually,  1 α  is very small, which means that the probability of the following error events  A ¯ k , i  ( A k , i  is the trigger event caused by the  i th  element of the state vector at time  k , i = 1 n ) is low:
P { A ¯ k , i } = P { | γ k , i | Z α / 2 | x ^ k | k τ x ^ k | k } = 1 α

2.2. Kalman Filter Algorithm Based on Event-Triggered Mechanism

Since indiscriminate filter calculations and trigger judgments of smart nodes are described in Section 2.1, only the algorithm of the receiving center is introduced in this section. In the receiving center, the Kalman filter algorithm based on the event-triggered mechanism (ET-KF) can be expressed as:
x ^ k = s k x ^ k | k + ( 1 s k ) x ^ k | k τ k , k = 1 , 2 ,
or recursively calculated as:
x ^ k = s k x ^ k | k + ( 1 s k ) Φ x ^ k 1 , k = 1 , 2 ,  
where  x ^ k  is the final estimate of ET-KF in the receiving center at time  k k τ k  is the last trigger time before time  k s k  is the trigger state of smart sensor nodes at time  k ; i.e.,  s k = 1  means that the transmission event is triggered and the receiving center receives the filter  x ^ k | k  of smart sensor nodes, and  s k = 0  means that the transmission event is not triggered and the receiving center will predict, that is
s k = { 0 , | γ k , i | Z α / 2 , i = 1 , , n 1 , otherwise  
Therefore, the flow chart of ET-KF is shown in Figure 1.

3. Performance Analysis of the ET-KF

From Equation (26), the probability  P { A k , i }  of the trigger event is  α , and the event  A k , i  is determined by  | γ k | i , which is the linear combination of measurements  z 0 ~ k . When  rank ( z k ) = r n , it means that there are  n r  events that are completely related to other events. Note that the event-triggered statistic  γ k  obeys the standard Gaussian distribution (for events subject to normal distribution, linear independence is equivalent to mutual independence), that is,  E { γ k , i γ k , j } = 0 , i j  and  E { γ k , i γ k , i } = 1 , then the probability  P { A k , i A k , j } = 0 , i j . Selecting  r  linearly unrelated events  A k , 1 , A k , 2 A k , r , according to the trigger mechanism in (27)–(29), the event-triggered probability can be calculated by:
p = P { A k , 1 A k , 2 A k , r } = 1 P { A ¯ k , 1 A ¯ k , 2 A ¯ k , r } = 1 P { A ¯ k , 1 } P { A ¯ k , 2 } P { A ¯ k , r } = 1 ( 1 α ) r
Therefore,  P { s k = 1 } = p P { s k = 0 } = 1 p . The estimated error variance  P k | k  of the ET-KF in the receiving center is given by the following theorem.
Theorem 2.
The estimation error variance   P k  of the ET-KF in the receiving center can be calculated as:
(1) If the system in Equations (1) and (2) has a steady state, that is,  P k  has a steady state value  P ¯ , then  P k  can be calculated by the following Riccati equation:
= p P ¯ + ( 1 p ) Φ Φ T + ( 1 p ) Φ Γ Q k 1 Γ T + ( 1 p ) Γ Q k 1 Γ T Φ T + ( 1 p ) Γ Q k 1 Γ T
where
lim k P k =
(2) If the system is a time-varying system,   P k  can be approximately calculated as:
P k p P k | k + ( 1 p ) p P k | k 1
Proof. 
If the system has a steady state estimation and  P k | k  has a stationary value  P ¯ , then based on Equation (28),  P k  can be calculated as:
P k = E { x ˜ k x ˜ k T } = E { ( x k x ^ k ) ( ) T } = E { ( s k ( x k x ^ k | k ) + ( 1 s k ) ( x k Φ x ^ k 1 ) ) ( ) T } = E { ( s k ( x k x ^ k | k ) + ( 1 s k ) ( Φ x k 1 + Γ w k 1 Φ x ^ k 1 ) ) ( ) T } = E { ( s k x ˜ k | k + ( 1 s k ) ( Φ x ˜ k 1 + Γ w k 1 ) ) ( ) T }
where ‘ ( ) ’ means the same as the previous formula. Since  s k  obeys the Bernoulli distribution,
E { s k ( 1 s k ) x ˜ k | k ( Φ x ˜ k 1 ) T } = 0
Considering that  w k  is white noise, from Equation (35) and  E { s k 2 = 1 } = p E { ( 1 s k ) 2 } = 1 p , there are:
P k = E { s k 2 x ˜ k | k x ˜ k | k T } + E { ( 1 s k ) 2 Φ x ˜ k 1 x ˜ k 1 T Φ T } + E { ( 1 s k ) 2 Φ Γ w k 1 w k 1 T Γ T } + E { ( 1 s k ) 2 Γ w k 1 w k 1 T Γ T Φ T } + E { ( 1 s k ) 2 Γ w k 1 w k 1 T Γ T } = p P k | k + ( 1 p ) Φ P k 1 Φ T + ( 1 p ) Φ Γ Q k 1 Γ T + ( 1 p ) Γ Q k 1 Γ T Φ T + ( 1 p ) Γ Q k 1 Γ T
the limit value exists in Riccati Equation (36):
lim k P k =
where   is the solution of the following steady-state Riccati Equation:
= p P ¯ + ( 1 p ) Φ Φ T + ( 1 p ) Φ Γ Q k 1 Γ T + ( 1 p ) Γ Q k 1 Γ T Φ T + ( 1 p ) Γ Q k 1 Γ T
If the system is a time-varying system, from Equation (28):
x ^ k = s k x ^ k | k + ( 1 s k ) Φ [ s k 1 x ^ k 1 | k 1 + ( 1 s k 1 ) Φ x ^ k 2 ] = s k x ^ k | k + ( 1 s k ) s k 1 x ^ k | k 1 + ( 1 s k ) ( 1 s k 1 ) Φ 2 x ^ k 2
the error at time  k  is:
x ˜ k = x k x ^ k = s k ( x k x ^ k | k ) + ( 1 s k ) s k 1 ( x k x ^ k | k 1 ) + ( 1 s k ) ( 1 s k 1 ) ( x k Φ 2 x ^ k 2 ) = s k x ˜ k | k + ( 1 s k ) s k 1 x ˜ k | k 1 + ( 1 s k ) ( 1 s k 1 ) ( x k Φ 2 x ^ k 2 )
since
E { s k ( 1 s k ) s k 1 x ˜ k | k ( x ˜ k | k 1 ) T } = 0 E { ( 1 s k ) 2 s k 1 ( 1 s k 1 ) x ˜ k | k 1 ( x k Φ 2 x ^ k 2 ) T } = 0 E { s k ( 1 s k ) ( 1 s k 1 ) x ˜ k | k x ˜ k | k ( x k Φ 2 x ^ k 2 ) T } = 0
  P k  can written as:
P k = E { x ˜ k x ˜ k T } = E { s k 2 x ˜ k | k x ˜ k | k T } + E { ( 1 s k ) 2 s k 1 2 x ˜ k | k 1 x ˜ k | k 1 T } + E { ( 1 s k ) 2 ( 1 s k 1 ) 2 ( x k Φ 2 x ^ k 2 ) ( ) T } = p P k | k + ( 1 p ) p P k | k 1 + ( 1 p ) 2 E { ( x k Φ 2 x ^ k 2 ) ( ) T }
since  ( 1 p )  is usually very small ( p = 1 ( 1 α ) r ), the third term of Equation (42) is omitted, and there are:
P k p P k | k + ( 1 p ) p P k | k 1
Remark 2.
The estimation error variance  P k  in the receiving center under the event-triggered threshold  Z α / 2  is given in Theorem 2. That is, under the condition that the system has a steady-state filter,  P k  can be computed by Equation (31) and Equation (32);on the contrary, the approximate calculation is (33). In addition, for a system with steady-state estimation, we can set the threshold reasonably according to the desired theoretical accuracy  P k . It can be seen from Equations (6) to (13) and (31) to (33) that the final estimation accuracy  P k  is not only related to the trigger rate, but also directly proportional to the system noise statistics  Q k  and measurement noise statistics  R k .

4. Simulation Examples

  • Scenario 1.(Event-triggered steady-state Kalman filter.) Consider the following planar tracking system:
x k = [ 1 T 0 0 0 1 0 0 0 0 1 T 0 0 0 1 ] x k + [ 0.5 T 2 0 T 0 0 0.5 T 2 0 T ] w k , k = 1 , , l
where  x k = [ x k x ˙ k y k y ˙ k ] T  is the state of target movement,  x k y k x ˙ k  and  y ˙ k  are the positions and velocities of the target in  x - a x i s  and  y - a x i s , respectively. The measurement equation of the system is:
z k = [ 1 0 0 0 0 0 1 0 ] x k + v k , k = 1 , , l
The estimation performance is mean square error (MSE) [25,26]:
MSE k = 1 k t = 0 k 1 N j = 1 N ( x t j x ^ t j ) 2 , k = 1 , , l
where  x t j  and  x ^ t j  are the true states and their estimations of the  j th  Monte Carlo experiment at time  t .
In the simulation, the working period is  T = 0.5   s  and the length of work is  l = 200  steps, the variance of processing noise is  Q = d i a g ( 0.2   m 2 / s 4 , 0.2   m 2 / s 4 ) , and the variance of measurement noise is  R = d i a g ( 0.4   m 2 / s 4 , 0.4   m 2 / s 4 ) , the initial state is  x 0 = [ 10   m ,   1   m / s ,   10   m ,   1   m / s ] T . A total of 200 Monte Carlo experiments are carried out. Because the measurement Equation (45) is a two-dimensional system and only the positions on the  x - a x i s  and  y - a x i s  are measured, according to Equation (26), the theoretical trigger probability is  p = 1 ( 1 α ) 2 . Here,  α = 0.98 ,   0.8 ,   0.6 ,   0.4  are chosen for comparison; that is, the thresholds are 0.05, 0.25, 0.52, and 0.84; and the theoretical trigger rates are 0.9996, 0.96, 0.84, and 0.64. For the threshold  Z α / 2 , the theoretical trigger probability and the actual trigger frequency are shown in Table 1. The theoretical trigger probability is consistent with the result of the actual trigger frequency, which verifies the correctness of Theorem 1 and Corollary 1. It shows that, when the proposed trigger mechanism is used, the actual trigger frequency can be predicted, and communication resources can be adjusted in advance to reduce the possibility of network congestion.
The error covariance  P k  calculated by Equation (31) and  MSE k  obtained by Equation (46) are shown in Table 2. It can be seen that the theoretical accuracy (diagonal elements of error covariance matrix  P k ) is close to  MSE k  obtained in actual work, which verifies the correctness of theorem 2 under the steady-state situation. When adopting the proposed trigger mechanism, the proposed ET-KF accuracy can be pre-calculated, and then the most suitable trigger frequency can be obtained by a reasonable trigger threshold.
  • P k , i ( i = 1 , , n )  and  MSE 200 , i  is the  i , i  diagonal element in the error covariance matrix  P k  of the ET-KF, and  MSE 200 , i  is the  i  element of  MSE 200  for the ET-KF.
  • Scenario 2. (Event-triggered time-varying Kalman filter.) Consider the planar tracking system in Equation (44) and Equation (45), where  x 0 = [ 10   m   1   m / s   10   m   1   m / s ]  is the initial state,  Q k = diag ( 0.2 + 0.2 sin ( 2 π k / 100 ) m 2 / s 4 ,     0.2 + 0.2 sin ( 2 π k / 100 ) m 2 / s 4 )  is processing noise variance,  R k = diag ( 0.3 + 0.5 sin ( 2 π k / 100 ) m 2 / s 4   , 0.3 + 0.5 sin ( 2 π k / 100 ) m 2 / s 4 )  is measurement noise variance. The performance of the system is measured by the accumulated mean square error (AMSE) [27,28,29]:
AMSE k = t = 1 k 1 N j = 1 N ( x t j x ^ t j ) T ( x t j x ^ t j )
where  x t j  and  x ^ t j  are the true states and their estimations of the  j th  Monte Carlo experiment at time  t . Corresponding to AMSE, the theoretical estimation accuracy is measured by  Σ tr P k :
Σ tr P k = t = 1 k 1 N j = 1 N tr ( P t j )
where  tr ( P t j )  is the trace of the estimated error covariance matrix at time  t  of the  j th  Monte Carlo experiment.
In Table 3, it is shown that for time-varying systems, the theoretical trigger frequency is close to the actual trigger frequency, which verifies the correctness of Theorem 1 and Corollary 1.
In Figure 2, with a different  α , the comparison curves of AMSE of ET-KF are shown. It can be seen that the smaller  α  is, the greater the threshold  Z α / 2  is, the lower the trigger rate is, and the worse the accuracy AMSE is. Comparison curves between the theoretical approximate estimation accuracy  Σ tr P k  (the sum of trace of  P k ) calculated by (33) and the AMSE obtained by the experiment with  α = 0.60  are shown in Figure 3. It can be seen that the theoretical estimation accuracy is close to the AMSE obtained by the experiment, which verifies the correctness of Theorem 2 under the time-varying situation. The AMSEs in Figure 2 and Figure 3 rise in a straight line, which means that the corresponding MSEs are stable.
In Table 4 and Figure 4, ET-KF is compared with the event trigger of innovation standardization (IS-KF) [14], a dynamic event trigger constructed by the linear change in the rate of change (LC-KF) [22]. From Table 4, under the same trigger rate (80%), the  AMSE 200  (AMSE at step 200) of the proposed trigger mechanism is lower than the other two trigger mechanisms. Figure 4 shows the comparison curve of AMSE at different steps, as the number of steps increases, the AMSE value gap based on different trigger mechanisms becomes larger. Among them, the AMSE value of the ETKF mechanism is the smallest and the estimation accuracy is the highest, which verifies the correctness of the conclusion.

5. Conclusions

In estimation, for linear systems, an event-triggered Kalman filter (ET-KF) is proposed in this paper. The main work is as follows:
  • The event-triggered statistic is constructed, which proves that the statistic obeys the standard Gaussian distribution, according to the event-triggered statistic and hypothesis test of the Gaussian distribution, the significance of the event-triggered threshold is given, and then an event-triggered estimation mechanism is designed.
  • Based on the event-triggered threshold and mechanism proposed in this paper, the theoretical trigger frequency under different thresholds and the estimation accuracy of event-triggered systems are analyzed. The proposed ET-KF can accurately set the event trigger frequency in advance. For linear systems with steady-state estimation, the estimation accuracy can be obtained by the Riccati equation accurately, and the trigger threshold can be set according to the accuracy. For linear time-varying systems without steady-state estimation, the approximate estimation accuracy can be obtained.
  • The proposed trigger mechanism has higher estimation accuracy at the same trigger rate, and the trigger setting is reasonable.

Author Contributions

Conceptualization and methodology: X.L. and G.H.; software and validation: X.L.; formal analysis and writing—original draft preparation: G.H.; writing—review and editing supervision: X.L. and G.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work is funded by the National Natural Science Foundation (NNSF) of China under Grant No.61503127, by the University Basic Research Foundation of Heilongjiang Province 2021-KYYWF-0027, by Outstanding Youth Foundation of Heilongjiang University, by Key Laboratory of Information Fusion Estimation and Detection, Heilongjiang Province.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the Key Laboratory of Information Fusion Estimation and Detection, Heilongjiang Province, for providing the condition used in this research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wei, W.; Song, H.; Li, W.; Shen, P.; Vasilakos, A. Gradient-driven parking navigation using a continuous information potential field based on wireless sensor network. Inf. Sci. 2017, 408, 100–114. [Google Scholar] [CrossRef]
  2. Elamin, A.; Abdelaziz, N.; El-Rabbany, A. A GNSS/INS/LiDAR Integration Scheme for UAV-Based Navigation in GNSS-Challen. Sensors 2022, 22, 9908. [Google Scholar] [CrossRef]
  3. Zhang, C.; Qin, J.; Li, H.; Wang, Y.; Wang, S.; Zheng, W.X. Consensus-based distributed two-target tracking over wireless sensor networks. Automatica 2022, 146, 110593. [Google Scholar] [CrossRef]
  4. Lu, Y.; Li, B.; Karimi, H.R.; Zhang, N. Measurement outlier-Resistant target tracking in wireless sensor networks with energy harvesting constraints. J. Frankl. Inst. 2022; in press. [Google Scholar] [CrossRef]
  5. Marsh, B.; Sadka, A.H.; Bahai, H. A Critical Review of Deep Learning-Based Multi-Sensor Fusion Techniques. Sensors 2022, 22, 9364. [Google Scholar] [CrossRef]
  6. GHorbani, S.; Safavi, A.A.; Naghavi, S.V. Decentralized event-triggered robust MPC for large-scale networked Lipchitz non-linear control systems. IET Control Theory Appl. 2021, 15, 2316–2333. [Google Scholar] [CrossRef]
  7. Ma, C.; Zheng, S.; Xu, T.; Ji, Y. Finite-Time Asynchronous Event-Triggered Formation of UAVs with Semi-Markov-Type Topologies. Sensors 2022, 22, 4529. [Google Scholar] [CrossRef]
  8. Yang, X.; Wang, H.; Zhu, Q. Event-triggered predictive control of nonlinear stochastic systems with output delay. Automatica 2022, 140, 110230. [Google Scholar] [CrossRef]
  9. Chen, W.; Shi, D.; Wang, J.; Shi, L. Event-triggered state estimation: Experimental performance assessment and comparative study. IEEE Trans. Control. Syst. Technol. 2016, 25, 1865–1872. [Google Scholar] [CrossRef]
  10. Shoukry, Y.; Tabuada, P. Event-triggered state observers for sparse sensor noise/attacks. IEEE Trans. Autom. Control 2015, 61, 2079–2091. [Google Scholar] [CrossRef]
  11. Ho, Y.C.; Cao, X.; Cassandras, C. Infinitesimal and finite perturbation analysis for queueing networks. Automatica 1983, 19, 439–445. [Google Scholar] [CrossRef]
  12. Miskowicz, M. Send-on-delta concept: An event-based data reporting strategy. Sensors 2006, 6, 49–63. [Google Scholar] [CrossRef] [Green Version]
  13. Trimpe, S.; D’Andrea, R. Event-based state estimation with variance-based trigger. IEEE Trans. Autom. Control 2014, 59, 3266–3281. [Google Scholar] [CrossRef]
  14. Wang, N.; Sun, S. Event-triggered sequential fusion filters based on estimators of observation noises for multi-sensor systems with correlated noises. Digit. Signal Process. 2021, 111, 102960. [Google Scholar] [CrossRef]
  15. Zhang, X.M.; Han, Q.L. Event-based H∞ filtering for sampled-data systems. Automatica 2015, 51, 55–69. [Google Scholar] [CrossRef]
  16. Sun, Y.C.; Yang, G.H. Event-triggered remote state estimation for cyber-physical systems under malicious DoS attacks. Inf. Sci. 2022, 602, 43–56. [Google Scholar] [CrossRef]
  17. Liu, C.; Sun, S. Event-triggered optimal and suboptimal distributed Kalman consensus filters for sensor networks. J. Franklin. Inst. 2021, 358, 5163–5183. [Google Scholar] [CrossRef]
  18. Zhong, X.; He, H. An event-triggered ADP control approach for continuous-time system with unknown internal states. IEEE Trans. Cybern. 2016, 47, 683–694. [Google Scholar] [CrossRef]
  19. Chen, Z.; Niu, B.; Zhao, X.; Zhang, L.; Xu, N. Model-based adaptive event-triggered control of nonlinear continuous-time systems. Appl. Math. Comput. 2021, 408, 126330. [Google Scholar] [CrossRef]
  20. Hu, S.; Qiu, J.; Chen, X.; Zhao, F.; Jiang, X. Dynamic event-triggered control for leader-following consensus of multiagent systems with the estimator. IET Control Theory Appl. 2022, 16, 475–484. [Google Scholar] [CrossRef]
  21. Zhu, X.; Zhang, X.; Wei, J.; Du, H.; Xi, J. Dynamic event-triggered control of networked control systems with uncertainty and transmission delay. J. Franklin. Inst. 2022, 359, 477–491. [Google Scholar] [CrossRef]
  22. Long, J.; Wang, W.; Wen, C.; Huang, J.; Lü, J. Output feedback based adaptive consensus tracking for uncertain heterogeneous multi-agent systems with event-triggered communication. Automatica 2022, 136, 110049. [Google Scholar] [CrossRef]
  23. Gao, G.; Gao, B.; Gao, S.; Hu, G.; Zhong, Y. A Hypothesis Test-Constrained Robust Kalman Filter for INS/GNSS Integration with Abnormal Measurement. IEEE Trans. Veh. Technol. 2022; early access. [Google Scholar] [CrossRef]
  24. Gao, G.; Zhong, Y.; Gao, S.; Gao, B. Double-channel sequential probability ratio test for failure detection in multisensor integrated systems. IEEE Trans. Instrum. Meas. 2021, 70, 1–14. [Google Scholar] [CrossRef]
  25. Zhong, M.; Ding, S.X.; Han, Q.L.; He, X.; Zhou, D. A Krein space-based approach to event-triggered H∞ filtering for linear discrete time-varying systems. Automatica 2022, 135, 110001. [Google Scholar] [CrossRef]
  26. Cacace, F.; Conte, F.; d’Angelo, M.; Germani, A.; Palombo, G. Filtering linear systems with large time-varying measurement delays. Automatica 2022, 136, 110084. [Google Scholar] [CrossRef]
  27. Zhao, H.; Tian, B.; Chen, B. Robust stable iterated unscented Kalman filter based on maximum correntropy criterion. Automatica 2022, 142, 110410. [Google Scholar] [CrossRef]
  28. Gao, S.; Zhong, Y.; Li, W. Robust adaptive filtering method for SINS/SAR integrated navigation system. Aerosp. Sci. Technol. 2011, 15, 425–430. [Google Scholar] [CrossRef] [Green Version]
  29. Zhu, X.; Gao, B.; Zhong, Y.; Gu, C.; Choi, K.S. Extended Kalman filter for online soft tissue characterization based on Hunt-Crossley contact model. J. Mech. Behav. Biomed. Mater. 2021, 123, 104667. [Google Scholar] [CrossRef]
Figure 1. The flow chart of ET-KF.
Figure 1. The flow chart of ET-KF.
Sensors 23 02202 g001
Figure 2. AMSE of the ET-KF algorithm with different  α .
Figure 2. AMSE of the ET-KF algorithm with different  α .
Sensors 23 02202 g002
Figure 3. AMSE of ET-KF algorithm and  Σ tr P k  with  α = 0.60 .
Figure 3. AMSE of ET-KF algorithm and  Σ tr P k  with  α = 0.60 .
Sensors 23 02202 g003
Figure 4. AMSE of different trigger threshold algorithms.
Figure 4. AMSE of different trigger threshold algorithms.
Sensors 23 02202 g004
Table 1. Comparison of trigger frequency in the steady-state system for different  α .
Table 1. Comparison of trigger frequency in the steady-state system for different  α .
  α 0.980.80.60.4
Threshold  Z α / 2  0.050.250.520.84
Theoretical trigger probability0.99960.960.840.64
Actual trigger frequency99.5%96%83.5%63%
Table 2. The error covariance of the steady-state ET-KF algorithm and  MSE i .
Table 2. The error covariance of the steady-state ET-KF algorithm and  MSE i .
  P k , 1   MSE 200 , 1   P k , 2   MSE 200 , 2   P k , 3   MSE 200 , 3   P k , 4   MSE 200 , 4
  α = 0.98 0.17950.17900.14500.15440.17950.18160.14500.1345
  α = 0.80 0.18480.17920.14500.15480.18480.18070.14500.1342
  α = 0.60 0.20840.18350.14500.15940.20840.18490.14500.1365
  α = 0.40 0.26140.19880.14500.16620.26140.18500.14500.1362
Table 3. Comparison of trigger frequency in the time-varying system under different  α  values.
Table 3. Comparison of trigger frequency in the time-varying system under different  α  values.
  α 0.980.80.60.4
Threshold  Z α / 2  0.050.250.520.84
Theoretical trigger probability0.99960.960.840.64
Actual trigger frequency99.5%95.5%82%65%
Table 4. The performance comparison of different trigger threshold algorithms.
Table 4. The performance comparison of different trigger threshold algorithms.
Different Trigger Threshold Algorithms   AMSE 200 Trigger Frequency
ET-KF19880%
IS-KF22580%
LC-KF21080%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, X.; Hao, G. Event-Triggered Kalman Filter and Its Performance Analysis. Sensors 2023, 23, 2202. https://doi.org/10.3390/s23042202

AMA Style

Li X, Hao G. Event-Triggered Kalman Filter and Its Performance Analysis. Sensors. 2023; 23(4):2202. https://doi.org/10.3390/s23042202

Chicago/Turabian Style

Li, Xiaona, and Gang Hao. 2023. "Event-Triggered Kalman Filter and Its Performance Analysis" Sensors 23, no. 4: 2202. https://doi.org/10.3390/s23042202

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop