本文基于动态事件触发研究了受欺骗性攻击的非线性网络控制系统的故障检测问题。利用一个服从伯努利分布的随机变量来描述欺骗性攻击,设计了动态事件触发机制以节省网络资源,构造了相应的Lyapunov函数,得到了滤波误差系统均方渐近稳定且满足混合 H∞和无源性能的充分条件。故障检测滤波器参数的设计以线性矩阵不等式的形式给出。最后通过一个数值算例,比较了动态事件触发方案与静态事件触发方案,说明了所提方法的有效性。 This paper investigates the problem of fault detection in nonlinear networked control systems subject to deception attacks based on dynamic event-triggered mechanism. A random variable obeying the Bernoulli distribution is used to describe the deception attack. In order to save network resources, a dynamic event-triggered mechanism is designed, the corresponding Lyapunov function is constructed and sufficient conditions are obtained for the filtering error system to be mean-square asymptotically stable and to satisfy the mixed H∞and passive property. The design parameters of the fault detection filter are given in the form of linear matrix inequalities. Finally, the effectiveness of the proposed method is illustrated by comparing the dynamic event-triggered scheme with the static event-triggered scheme through a numerical arithmetic example.
本文基于动态事件触发研究了受欺骗性攻击的非线性网络控制系统的故障检测问题。利用一个服从伯努利分布的随机变量来描述欺骗性攻击,设计了动态事件触发机制以节省网络资源,构造了相应的Lyapunov函数,得到了滤波误差系统均方渐近稳定且满足混合H∞和无源性能的充分条件。故障检测滤波器参数的设计以线性矩阵不等式的形式给出。最后通过一个数值算例,比较了动态事件触发方案与静态事件触发方案,说明了所提方法的有效性。
网络控制系统,动态事件触发机制,欺骗性攻击,故障检测,混合H∞和无源性
Xuefen Hua, Weiqun Wang
School of Mathematics and Statistics, Nanjing University of Science and Technology, Nanjing Jiangsu
Received: Mar. 20th, 2024; accepted: Apr. 10th, 2024; published: Apr. 18th, 2024
This paper investigates the problem of fault detection in nonlinear networked control systems subject to deception attacks based on dynamic event-triggered mechanism. A random variable obeying the Bernoulli distribution is used to describe the deception attack. In order to save network resources, a dynamic event-triggered mechanism is designed, the corresponding Lyapunov function is constructed and sufficient conditions are obtained for the filtering error system to be mean-square asymptotically stable and to satisfy the mixed H∞and passive property. The design parameters of the fault detection filter are given in the form of linear matrix inequalities. Finally, the effectiveness of the proposed method is illustrated by comparing the dynamic event-triggered scheme with the static event-triggered scheme through a numerical arithmetic example.
Keywords:Networked Control Systems, Dynamic Event-Triggered Mechanism, Deception Attacks, Fault Detection, Mixed H∞and Passive Performance
Copyright © 2024 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
网络控制系统由于其成本低、易于安装、易于维护等优点,在过去的几十年里受到了广泛的关注,目前在多智能体系统、智能电网、无人机、工业自动化等工业有许多应用 [
由于网络的开放性,通信信道不可避免地会受到恶意攻击。欺骗性攻击是一种典型的网络攻击,研究者们也提出了一些有效的方法来对抗欺骗性攻击。文献 [
为了节约网络资源,事件触发机制已被广泛应用于网络控制系统。例如,文献 [
基于上述分析,考虑动态事件触发机制,研究欺骗性攻击下网络控制系统的混合H∞和无源故障检测是有实际价值的。本文的主要贡献总结如下:1) 比较了动态与静态事件触发方案,设计的动态事件触发方案可以更有效地利用网络资源。2) 获得了滤波误差系统在欺骗性攻击下均方渐近稳定的充分条件,并设计了混合H∞和无源性故障检测滤波器,能够快速检测故障。
考虑一类非线性网络控制系统:
{ x ˙ ( t ) = A x ( t ) + B v ( t ) + E f ( t ) + F g ( t , x ( t ) ) , y ( t ) = C x ( t ) , (1)
其中, x ( t ) ∈ R n 是状态向量, y ( t ) ∈ R m 是测量输出, f ( t ) ∈ R q 是待检测故障, v ( t ) 是满足集合 L d = { v ( t ) : 0 < ∫ 0 ∞ v T ( t ) v ( t ) d t ≤ d } 的外源性干扰信号, g ( t , x ( t ) ) 是一个非线性向量函数且满足:
g T ( t , x ( t ) ) g ( t , x ( t ) ) ≤ ϵ 2 x T ( t ) N T N x ( t ) , (2)
其中, ϵ > 0 为已知常数, A , B , C , E , F , N 为具有适当维数的已知矩阵。
当 f ( t ) ≠ 0 时,系统出现故障,为了检测故障,设计如下的故障检测滤波器:
{ x ˙ f ( t ) = A f x f ( t ) + B f y ^ ( t ) , r ( t ) = C f x f ( t ) , (3)
其中, x f ( t ) ∈ R n 为 x ( t ) 的状态估计, r ( t ) ∈ R p 为残差信号, y ^ ( t ) 表示故障检测滤波器的输入。矩阵 A f , B f , C f 是待确定的具有适当维数的系数矩阵。
为节省更多的网络资源,采用如下的动态事件触发方案:
e T ( t ) Φ e ( t ) ≥ δ ( t ) y T ( t k h ) Φ y ( t k h ) , (4)
其中, e ( t ) = y ( t ) − y ( t k h ) 测量误差, t k h 是最新的事件触发时刻,假设第一个事件发生在时刻 t 0 = 0 。 δ ( t ) 是动态事件触发参数, Φ 是一个待设计的正定权重矩阵。一旦满足了事件触发条件,测量的输出数据 y ( t k h ) 将被更新,如图1所示。下一个触发时刻 t k + 1 h 由以下触发条件决定:
图1. 基于事件触发的网络控制系统结构图
t k + 1 h = t k h + inf { l h | e T ( t ) Φ e ( t ) ≥ δ ( t ) y T ( t k h ) Φ y ( t k h ) } . (5)
动态事件触发参数 δ ( t ) 是时变的 [
δ ˙ ( t ) = u δ ( t ) [ 1 δ ( t ) − δ 0 ] e T ( t ) Φ e ( t ) , (6)
其中, u , δ 0 是给定的正常数。如果 u = 0 ,则上述动态事件触发方案将退化为具有固定阈值的静态事件触发方案。
当没有欺骗性攻击时,根据事件触发条件(5),滤波器的输入 y ^ ( t ) 为 y ( t k h ) ,并且保持不变直到下一个新的触发信号到达。用测量误差表示事件触发采样状态 y ( t k h ) ,即 y ( t k h ) = y ( t ) − e ( t ) 。
当存在欺骗性攻击时,攻击者通过在真实信号中注入一定的欺骗信号来破坏数据包。滤波器的真实输入描述如下:
y ^ ( t ) = [ 1 − α ( t ) ] q ( y ( t ) ) + α ( t ) y ( t k h ) , (7)
其中, q ( t ) 是欺骗信号, q ( t ) 满足以下约束条件:
‖ q ( t ) ‖ 2 ≤ ‖ Q y ( t ) ‖ 2 , (8)
Q > 0 是一个常数矩阵,用于表示欺骗性攻击函数的上界。
α ( t ) 是一个伯努利分布变量,满足以下分布:
{ P r o b { α ( t ) = 1 } = E { α ( t ) } = α 1 , P r o b { α ( t ) = 0 } = 1 − α 1 = α ¯ 1 , V a r { α ( t ) } = E { ( α ( t ) − α 1 ) 2 } = α 1 α ¯ 1 ,
其中, α 1 ∈ [ 0 , 1 ] 是已知常数。当 α 1 = 0 时,表示真实的输入被欺骗信号 q ( t ) 所取代。
记 ξ ( t ) = [ x T ( t ) , x f T ( t ) ] T , ω ( t ) = [ v T ( t ) , f T ( t ) ] T , r e ( t ) = r ( t ) − f ( t ) 。由方程(1)和(3),可以得到下列滤波误差系统:
{ ξ ˙ ( t ) = A ¯ ξ ( t ) + α ( t ) A 1 ξ ( t ) + B ¯ ω ( t ) − α ( t ) E ¯ e ( t ) + ( 1 − α ( t ) ) E ¯ q ( t ) + F ¯ g ( t , H ξ ( t ) ) , r e ( t ) = C ¯ ξ ( t ) − D ¯ ω ( t ) , (9)
其中,
A ¯ = [ A 0 0 A f ] , A 1 = [ 0 0 B f C 0 ] , B ¯ = [ B E 0 0 ] , E ¯ = [ 0 B f ] , F ¯ = [ F 0 ] , C ¯ = [ 0 C f ] , D ¯ = [ 0 I ] , H = [ I 0 ] .
根据残差信号,设计如下的残差估计函数:
F ( t ) = ( 1 t ∫ 0 t r T ( s ) r ( s ) d s ) 1 2 . (10)
残差估计阈值 F t h 为:
F t h = sup v ( t ) ∈ L d , f ( t ) = 0 F ( t ) . (11)
故障检测判据是:
{ F ( t ) > F t h ⇒ 报 警 , F ( t ) ≤ F t h ⇒ 不 报 警 . (12)
定义1. [
lim t → ∞ E { ‖ ξ ( t ) ‖ 2 } = 0 ,
则称系统(9)在均方意义上是渐近稳定的。
定义2. [
1) 当 ω ( t ) = 0 ,系统(9)是均方渐近稳定的。
2) 当 ω ( t ) ≠ 0 ,在零初始条件下有下式成立:
E { ∫ 0 ∞ ( − ρ r e T ( t ) r e ( t ) + 2 γ ( 1 − ρ ) r e T ( t ) ω ( t ) ) d t } ≥ E { − γ 2 ∫ 0 ∞ ω T ( t ) ω ( t ) d t } . (13)
当 ρ = 0 ,式(13)表示无源性能指标;当 ρ = 1 ,式(13)表示 H ∞ 性能指标;当 ρ ∈ ( 0 , 1 ) ,式(13)表示混合 H ∞ 和无源性能指标。
本文的目标是给出保证系统(9)在所设计的动态事件触发条件下是均方渐近稳定且满足混合 H ∞ 和无源性能的充分条件,并获得混合 H ∞ 和无源故障检测滤波器(3)的参数,使得出现故障时系统能快速报警。
下面的定理给出了保证滤波误差系统(9)在动态事件触发条件下是均方渐近稳定的且满足混合 H ∞ 和无源性能的充分条件。
定理1. 给定正数 α 1 , u , δ 0 , γ , ρ ,如果存在适当维数的对称正定矩阵 P , Φ ,使
Π = [ Π 11 * * * * * * * α 1 ( P E ¯ ) T − u δ 0 Φ * * * * * * Π 31 0 Π 33 * * * * * ( P F ¯ ) T 0 0 − I * * * * Π 51 0 0 0 − I * * * ρ C ¯ 0 − ρ D ¯ 0 0 − I * * u C ^ Φ − u Φ 0 0 0 0 − Φ * ϵ N ^ 0 0 0 0 0 0 − I ] < 0 , (14)
其中,
Π 11 = ( A ¯ + α 1 A 1 ) T P + P ( A ¯ + α 1 A 1 ) + Q ¯ T Q ¯ , Π 51 = α ¯ 1 ( P E ¯ ) T , Π 31 = ( P B ¯ ) T − γ ( 1 − ρ ) C ¯ , Π 33 = − γ 2 I + 2 γ ( 1 − ρ ) D ¯ T , C ^ = [ C 0 ] , N ^ = [ N 0 ] , Q ¯ = Q C ^ ,
则系统(9)是均方渐近稳定的且满足混合 H ∞ 和无源性能指标。
证明. 取Lyapunov 函数 V ( t ) ,
V ( t ) = ξ T ( t ) P ξ ( t ) + 1 2 δ 2 ( t ) , (15)
对 V ( t ) 沿着系统(9)的轨迹求导,可得:
V ˙ ( t ) = 2 ξ T ( t ) P ξ ˙ ( t ) + δ ˙ ( t ) δ ( t ) , (16)
要保证系统(9)是均方渐近稳定的,需要证明 E { V ˙ ( t ) } < 0 成立 [
E { V ˙ ( t ) + J ( t ) } ≤ 0 , J ( t ) = ρ r e T ( t ) r e ( t ) − 2 γ ( 1 − ρ ) r e T ( t ) ω ( t ) − γ 2 ω T ( t ) ω ( t ) . (17)
改写非线性项约束(2),有:
ϵ 2 ξ T ( t ) N ^ T N ^ ξ ( t ) − g T ( t , H ξ ( t ) ) g ( t , H ξ ( t ) ) ≥ 0 , (18)
改写欺骗信号约束(8),有:
ξ T ( t ) Q ¯ T Q ¯ ξ ( t ) − q T ( t ) q ( t ) ≥ 0 , (19)
根据事件触发条件(5),当
δ ˙ ( t ) δ ( t ) = u [ 1 δ ( t ) − δ 0 ] e T ( t ) Φ e ( t ) ≤ u y T ( t k h ) Φ y ( t k h ) − u δ 0 e T ( t ) Φ e ( t ) = u [ C ^ ξ ( t ) − e ( t ) ] T Φ [ C ^ ξ ( t ) − e ( t ) ] − u δ 0 e T ( t ) Φ e ( t ) , (20)
由此,可得:
E { V ˙ ( t ) + J ( t ) } ≤ E { V ˙ ( t ) + J ( t ) + ϵ 2 ξ T ( t ) N ^ T N ^ ξ ( t ) − g T ( t , H ξ ( t ) ) g ( t , H ξ ( t ) ) } ≤ 2 ξ T ( t ) P [ A ¯ ξ ( t ) + α 1 A 1 ξ ( t ) + B ¯ ω ( t ) − α 1 E ¯ e ( t ) + α ¯ 1 E ¯ q ( t ) + F ¯ g ( t , H ξ ( t ) ) ] + ρ [ C ¯ ξ ( t ) − D ¯ ω ( t ) ] T [ C ¯ ξ ( t ) − D ¯ ω ( t ) ] − 2 γ ( 1 − ρ ) [ C ¯ ξ ( t ) − D ¯ ω ( t ) ] T ω ( t )
− γ 2 ω T ( t ) ω ( t ) + u [ C ^ ξ ( t ) − e ( t ) ] T Φ [ C ^ ξ ( t ) − e ( t ) ] − u δ 0 e T ( t ) Φ e ( t ) + ϵ 2 ξ T ( t ) N ^ T N ^ ξ ( t ) − g T ( t , H ξ ( t ) ) g ( t , H ξ ( t ) ) + ξ T ( t ) Q ¯ T Q ¯ ξ ( t ) − q T ( t ) q ( t ) = ζ T ( t ) Ξ ζ ( t ) , (21)
其中,
ζ ( t ) = [ ξ T ( t ) , e T ( t ) , ω T ( t ) , g T ( t , H ξ ( t ) ) , q T ( t ) ] T ,
Ξ = [ Ξ 11 * * * * Ξ 21 − u δ 0 Φ * * * Ξ 31 0 Ξ 33 * * ( P F ¯ ) T 0 0 − I * Π 51 0 0 0 − I ] + ρ Ξ 1 T Ξ 1 + u Ξ 2 T Φ Ξ 2 ,
Ξ 11 = ( A ¯ + α 1 A 1 ) T P + P ( A ¯ + α 1 A 1 ) + ϵ 2 N ^ T N ^ + Q ¯ T Q ¯ , Ξ 21 = α 1 ( P E ¯ ) T , Ξ 31 = ( P B ¯ ) T − γ ( 1 − ρ ) C ¯ , Ξ 33 = − γ 2 I + 2 γ ( 1 − ρ ) G ¯ T , Π 51 = α ¯ 1 ( P E ¯ ) T , Ξ 1 = [ C ¯ 0 − D ¯ 0 0 ] , Ξ 2 = [ C ^ − I 0 0 0 ] .
对条件(14)中 Π < 0 应用Schur补定理 [
E { V ˙ ( t ) + J ( t ) } ≤ ζ T ( t ) Ξ ζ ( t ) < 0. (22)
当 ω ( t ) = 0 , J ( t ) = ρ r e T ( t ) r e ( t ) > 0 ,
E { V ˙ ( t ) } = E { J ( t ) + V ˙ ( t ) − J ( t ) } < 0. (23)
因此,由李亚普诺夫稳定性理论 [
当 ω ( t ) ≠ 0 ,在零初始条件下 E { V ( 0 ) } = 0 , E { V ( ∞ ) } ≥ 0 ,则:
E { ∫ 0 ∞ J ( t ) d t } ≤ − E { ∫ 0 ∞ V ˙ ( t ) d t } ≤ 0. (24)
因此,根据定义2,系统(9)是均方渐近稳定且满足混合 H ∞ 和无源性能指标的。证明完毕。
定理2. 对于给定正数 α 1 , u , δ 0 , γ , ρ ,如果存在适当维数的对称正定矩阵 X , Y , Φ 和矩阵 A F , B F , C F ,使得下列矩阵不等式成立:
Y − X < 0 , (25)
[ Π 11 * * * * * * * * * Π 21 Π 22 * * * * * * * * 0 α 1 B F − u δ 0 Φ * * * * * * * Π 41 B T X 0 − γ 2 I * * * * * * Π 51 E T X 0 Π 54 Π 55 * * * * * F T Y F T X 0 0 0 − I * * * * 0 α ¯ 1 B F 0 0 0 0 − I * * * Π 81 0 0 − ρ D ¯ 0 0 0 − I * * u C Φ u C Φ − u Φ 0 0 0 0 0 − Φ * ϵ N ϵ N 0 0 0 0 0 0 0 − I ] < 0 , (26)
其中,
Π 11 = Y A + A T Y + Q ˜ T Q ˜ , Π 21 = X A + A T Y + α 1 B F C + A F + Q ˜ T Q ˜ , Π 41 = B T Y − γ ( 1 − ρ ) C F , Π 22 = X A + A T X + α 1 B F C + α 1 C T B F T + Q ˜ T Q ˜ , Π 51 = E T Y − γ ( 1 − ρ ) C F , Π 54 = γ ( 1 − ρ ) D ¯ , Π 55 = 2 γ ( 1 − ρ ) D ¯ − γ 2 I , Π 81 = ρ C F , Q ˜ = Q C ,
则称系统(9)是均方渐近稳定且满足混合 H ∞ 和无源性能指标的,相应滤波器(3)的参数为:
A f = ( Y − X ) − 1 A F , B f = ( Y − X ) − 1 B F , C f = C F .
证明. 将定理1中P分块:
P = [ X 1 * X 2 X 3 ] , P − 1 = [ Y 1 * Y 2 Y 3 ] .
定义
M 1 = [ Y 1 I Y 2 0 ] , M 2 = [ I X 1 0 X 2 ] , J = d i a g { M 1 T , I , I , I , I , I , I , I } .
则有 P M 1 = M 2 ,矩阵 P > 0 成立当且仅当 M 1 T P M 1 > 0 ,即 X 1 − Y 1 − 1 > 0 。根据 P P − 1 = I ,有 X 1 Y 1 + X 2 T Y 2 = I ,故 X 2 , Y 2 为可逆矩阵。对条件(14)中每项左右两端分别乘以 J 和 J T ,之后再每项左右两边同乘以 d i a g { Y 1 − 1 , I , I , I , I , I , I , I , I , I } ,令:
Y = Y 1 − 1 , X = X 1 , A F = X 2 T A f Y 2 Y 1 − 1 , B F = X 2 T B f , C F = C f Y 2 Y 1 − 1 ,
可得式(26)成立。根据定理1,条件(25)和(26)能保证系统(9)是是均方渐近稳定且满足混合 H ∞ 和无源性能指标的,证明完毕。
本节给出一个例子来验证所设计的事件触发方案和故障检测滤波器的有效性。
例:系统(1)的系数矩阵如下:
A = [ 0 1 − 6 − 2 ] , B = [ 0.5 1 ] , C = [ − 0.1 1 ] , F = [ 1 0 0 1 ] , E = [ 1 1 ] .
取初始状态为 x ( 0 ) = [ − 0.64 , 0.51 ] T ,设外部扰动 v ( t ) = 0.2 e − 2 t ,非线性项 g T ( t , x ( t ) ) = [ 0.2 sin ( t ) , 0 ] ,令 γ = 2.267 , ρ = 0.5 , δ 0 = 0.6 ,故障信号为:
f ( t ) = { 1 , 3.5 ≤ t ≤ 5.5 , 0 , 其 他 .
利用MATLAB求解定理2中条件(25)和(26),得到事件触发参数矩阵为 Φ = 0.4489 ,相应混合 H ∞ 和无源故障检测滤波器参数如下:
A f = [ − 3.1037 − 1.6094 5.1596 − 1.4416 ] , B f = [ 0.4416 − 0.2280 ] , C f = [ − 0.4588 0.1302 ] .
图2(a)为残差信号,可以看出在有故障和无故障情况下残差信号最终均为零,说明设计的故障检测滤波器可以很好地估计故障。图2(b)为事件触发时刻,图3(b)为 E { ξ T ( t ) R ξ ( t ) } 的运动轨迹,虽然系统在3.5 s到5.5 s受到故障影响,但从轨迹图可以看出系统(9)是渐近稳定的。图3(a)为有故障和无故障情况下的残差估计函数 F ( t ) 与估计阈值 F t h ,当有故障发生时,残差估计函数 F ( t ) 在 t = 3.6 s 会超过估计阈值,系统会发生报警,说明设计的滤波器性能良好,可及时检测到故障。
对系统(9),当 u = 0 时,式(5)即为静态事件触发方案,动态与静态事件触发参数 δ 与事件触发次数之间的关系如表1所示。事件触发参数 δ 越大,事件触发次数越少,传输的数据越少。因此,与静态事件触发方案相比,动态事件触发方案(5)可以在保持系统性能的同时,减少信号的传输数量,节约了网络资源。
δ | 静态事件触发 | 动态事件触发 |
---|---|---|
δ = 0.1 | 68 | 36 |
δ = 0.3 | 47 | 34 |
δ = 0.5 | 39 | 30 |
表1. 不同 δ 下静态与动态事件触发次数
图2. 残差信号和事件触发时刻
图3. 残差估计函数与阈值、 E { ξ T ( t ) R ξ ( t ) } 的轨迹
本文研究了欺骗性攻击下的非线性网络控制系统的故障检测问题。利用伯努利分布的变量描述欺骗攻击,采用动态事件触发条件,构造了含动态事件触发参数的Lyapunov函数,给出了保证滤波误差系统均方渐近稳定且满足混合 和无源性的充分条件,设计了故障检测滤波器和事件触发参数。本文所设计的故障检测滤波器性能好,能及时检测到故障;与常用的静态事件触发条件比较,本文所设计的动态触发条件在保持系统性能的同时,减少了事件触发次数,节约了宝贵的网络资源。
在此对王为群教授对本篇论文提供的指导和帮助表示衷心的感谢。
花雪芬,王为群. 欺骗性攻击下基于动态事件触发的网络控制系统的混合H∞和无源故障检测Mixed H∞and Passive Fault Detection for Networked Control Systems Based on Dynamic Event-Triggered Mechanism under Deception Attacks[J]. 动力系统与控制, 2024, 13(02): 81-90. https://doi.org/10.12677/dsc.2024.132008
https://doi.org/10.1109/TIE.2015.2504351
https://doi.org/10.1016/j.rser.2017.05.092
https://doi.org/10.1109/TAC.2008.923688
https://doi.org/10.1109/TAC.2018.2856899
https://doi.org/10.1109/ACCESS.2019.2928306
https://doi.org/10.1016/j.sigpro.2014.12.014
https://doi.org/10.23919/ChiCC.2017.8028537
https://doi.org/10.1016/j.jfranklin.2020.12.018
https://doi.org/10.1016/j.amc.2017.09.007
https://doi.org/10.1016/j.jfranklin.2019.09.002
https://doi.org/10.1049/iet-cta.2016.0135
https://doi.org/10.1016/j.jfranklin.2018.01.048
https://doi.org/10.1016/j.jfranklin.2022.12.034
https://doi.org/10.1002/acs.3409
https://doi.org/10.1109/TCYB.2016.2631903
https://doi.org/10.1016/j.automatica.2010.03.007
https://doi.org/10.1109/TCSI.2013.2246213