近年来，由于工业生产以及航空航天等领域对于控制系统安全性和可靠性要求的不断提高，故障检测技术得到了快速的发展，并且涌现出了大量的故障检测方法^[1-2].故障检测滤波器方法是最早发展的一类基于模型的检测方法，其主要原理是利用系统解析表达式对系统的输出做出估计，然后与实际测量的输出值相比较得到残差，再经过残差评价以及阈值比较的过程，最终判断故障是否发生^[3].但在实际的故障检测过程中，未知的扰动信号以及不确定的建模误差的存在也会对残差产生影响，因此，需要设计鲁棒故障检测滤波器.目前使用较多的是基于鲁棒优化的方法，大致可分为两类:1)使用H_∞范数作为衡量故障敏感性和鲁棒性的性能指标，将鲁棒故障检测滤波器的设计问题转化为传递函数的H_∞范数最小值求解问题^[4-5]；2)使用混合H_{_}/H_∞性能指标来研究鲁棒故障检测滤波器的设计问题，将其转化为多目标优化的求解问题^[6-7].

与此同时，切换系统作为一种典型的混杂系统，正在获得越来越多的关注.其结构是由若干个描述连续动态的子系统以及一个表示离散动态的切换信号组成.切换信号的作用是决定子系统的切换时刻和次序，并且保证在某个时刻，整个系统中只有一个子系统激活.目前，有关切换系统的研究比较活跃，涌现出了很多关于切换系统的研究成果^[8-10].并且，随着切换系统的应用领域日益增多，其故障检测问题已经引起了研究者的兴趣，并且取得了一些研究结果.文献^[11]针对带范数有界的模型不确定性和状态时滞的连续时间切换系统，将其故障检测问题转化为H_∞滤波问题.文献^[12]使用H_∞滤波器方法实现了具有区间时变时滞的离散时间切换系统的故障检测.文献^[13]则针对离散时间切换系统，研究了有限频率范围内的故障检测问题.除了使用H_∞滤波器的方法进行故障检测之外，文献^[14]考虑了在混合H_{_}/H_∞性能指标框架下，离散时间不确定切换系统鲁棒故障检测滤波器的设计.文献^[15]则考虑了混合H_{_}/H_∞性能指标下，连续时间切换系统鲁棒故障检测滤波器的设计,但在文献^[14-15]中，由于Lyapunov矩阵和系统矩阵存在耦合，得到的结果具有一定的保守性.

基于此，本文将考虑在混合H_{_}/H_∞框架下，切换系统鲁棒故障检测滤波器的设计问题.首先，利用平均驻留时间技术和多Lyapunov函数方法，得到了表示鲁棒性的加权H_∞性能指标和表示敏感性的H_{_}性能指标的线性矩阵不等式.其次，使用加入松弛变量的方法，解除了Lyapunov矩阵和系统矩阵的耦合，降低了结果的保守性.

1 问题形成

考虑离散时间线性切换系统，模型表示如下：

$\begin{align} & x\left( k+1 \right)={{A}_{\sigma \left( k \right)}}x\left( k \right)+{{B}_{\sigma \left( k \right)1}}u\left( k \right)+{{D}_{\sigma \left( k \right)1}}d\left( k \right)+{{F}_{\sigma \left( k \right)1}}f\left( k \right), \\ & y\left( k \right)={{C}_{\sigma \left( k \right)}}x\left( k \right)+{{B}_{\sigma \left( k \right)2}}u\left( k \right)+{{D}_{\sigma \left( k \right)2}}d\left( k \right)+{{F}_{\sigma \left( k \right)2}}f\left( k \right). \\ \end{align}$

式中：x(k)∈Rⁿ为系统的状态向量；y(k)∈R^m为系统的测量输出向量；d(k)∈R^p为干扰向量；f(k)∈R^q为故障，假设d(k)和f(k)都是l₂范数有界的信号.σ(k)为切换系统的切换信号，是一个分段常值函数，取值范围是自然数集合N={1,2,…,K}.(A_σ(k),B_σ(k)1,B_σ(k)2,C_σ(k),D_σ(k)1,D_σ(k)2,F_σ(k)1,F_σ(k)2)为每个子系统的系数矩阵，是具有适当维数且数值已知的常数矩阵.当切换信号σ(k)=i时，说明第i个子系统被激活，则相应的系统模型转化为

$\begin{align} & x\left( k+1 \right)={{A}_{i}}x\left( k \right)+{{B}_{i1}}u\left( k \right)+{{D}_{i1}}d\left( k \right)+{{F}_{i1}}f\left( k \right), \\ & y\left( k \right)={{C}_{i}}x\left( k \right)+{{B}_{i2}}u\left( k \right)+{{D}_{i2}}d\left( k \right)+{{F}_{i2}}f\left( k \right). \\ \end{align}$

(1)

系统(1)设计故障检测滤波器为

$\begin{align} & \hat{x}\left( k+1 \right)={{A}_{i}}\hat{x}\left( k \right)+{{B}_{i1}}u\left( k \right)+{{L}_{i}}\left( y\left( k \right)-\hat{y}\left( k \right) \right), \\ & \hat{y}\left( k \right)={{C}_{i}}\hat{x}\left( k \right)+{{B}_{i2}}u\left( k \right), \\ & r\left( k \right)={{V}_{i}}\left( y\left( k \right)-\hat{y}\left( k \right) \right). \\ \end{align}$

(2)

式中：r(k)∈R^m为残差向量；L_i∈R^n×m为滤波器的增益矩阵；V_i∈R^m×m为残差的权值矩阵；L_i、V_i分别为需要设计的滤波器参数.

根据式(1)、(2)可以得到估计误差系统：

$\begin{align} & e\left( k+1 \right)=\left( {{A}_{i}}-{{L}_{i}}{{C}_{i}} \right)e\left( k \right)+\left( {{D}_{i1}}-{{L}_{i}}{{D}_{i2}} \right)d\left( k \right)+ \\ & \left( {{F}_{i1}}-{{L}_{i}}{{D}_{i2}} \right)f\left( k \right), \\ & r\left( k \right)={{V}_{i}}{{C}_{i}}e\left( k \right)+{{V}_{i}}{{C}_{i2}}d\left( k \right)+{{V}_{i}}{{F}_{i2}}f\left( k \right), \\ \end{align}$

其中e(k)=x(k)－${\hat{x}}$(k)，再进一步得到：

$\begin{align} & e\left( k+1 \right)={{{\tilde{A}}}_{i}}e\left( k \right)+{{{\tilde{D}}}_{i1}}d\left( k \right)+{{{\tilde{F}}}_{i1}}f\left( k \right), \\ & r\left( k \right)={{{\tilde{C}}}_{i}}e\left( k \right)+{{{\tilde{D}}}_{i2}}d\left( k \right)+{{{\tilde{F}}}_{i2}}f\left( k \right). \\ \end{align}$

(3)

式中：${{{\tilde{A}}}_{i}}={{A}_{i}}-{{L}_{i}}{{C}_{i}},{{{\tilde{D}}}_{i1}}={{D}_{i1}}-{{L}_{i}}{{D}_{i2}},{{{\tilde{F}}}_{i1}}={{F}_{i1}}-{{L}_{i}}{{D}_{i2}},$${{{\tilde{C}}}_{i}}={{V}_{i}}{{C}_{i}},{{{\tilde{D}}}_{i2}}={{V}_{i}}{{D}_{i2}},{{{\tilde{F}}}_{i2}}={{V}_{i}}{{F}_{i2}}$.

可知误差系统也为切换系统，并假设其与原系统同步，残差信号r(k)与干扰d(k)和故障f(k)有关，为了得到设计鲁棒故障检测滤波器，需要设计滤波器参数满足以下3个条件：

1) 误差系统是全局指数稳定的；

2) 残差r(k)对干扰d(k)具有鲁棒性，即满足加权H_∞性能指标；

3) 残差r(k)对故障f(k)具有敏感性，即满足H_{_}性能指标.

下面用定义来描述上述的3个条件.

定义1 如果存在常数K>0,0＜β＜1，使得系统(3)的解e(k)满足‖e(k)‖≤Kβ^(k－k₀)‖e(k₀)‖，∀k≥k₀，则在切换信号σ(k)下，系统(3)是全局指数稳定的.

定义2 给定常数0＜α＜1和γ>0，当满足零初始条件和f(k)=0时，对于所有的非零d(k)∈l₂[0,∞),若以下不等式成立

$\sum\limits_{k=0}^{\infty }{{{\left( 1-\alpha \right)}^{k}}{{r}^{T}}\left( k \right)r\left( k \right)}\le {{\gamma }^{2}}\sum\limits_{k=0}^{\infty }{{{d}^{T}}\left( k \right)d\left( k \right)},$

(4)

则称系统(3)的残差r(k)对于干扰d(k)具有加权的H_∞性能指标，其中α为系统的衰减度.

定义3 给定常数β>0，当满足零初始条件和d(k)=0时，对于所有的非零f(k)∈l₂［0,∞),若以下不等式成立

$\sum\limits_{k=0}^{\infty }{{{r}^{T}}\left( k \right)r\left( k \right)}\ge {{\beta }^{2}}\sum\limits_{k=0}^{\infty }{{{f}^{T}}\left( k \right)f\left( k \right)},$

(5)

则称系统(3)的残差r(k)对于故障f(k)具有H_{_}性能指标.

另外，考虑切换系统的切换信号满足平均驻留时间，定义如下.

定义4^［9] 对于满足k₀＜k_s＜k_v的任意时刻，令N_σ(k_s,k_v)表示σ(k)在时间区间(k_s,k_v)内的切换次数，如果对于给定的τ_α>0，N₀≥0，可以得到：

${{N}_{\sigma }}\left( {{k}_{s}},{{k}_{v}} \right)\le {{N}_{0}}+\frac{{{k}_{v}}-{{k}_{s}}}{{{\tau }_{\alpha }}},$

成立，则τ_α称为切换信号σ(k)的平均驻留时间(average dwell time，ADT)，N₀称为抖振边界，通常情况下，令N₀=0.

2 鲁棒故障检测滤波器的设计

考虑在H_{_}/H_∞框架下，鲁棒故障检测滤波器的设计问题.首先将通过两个引理，得到初步的鲁棒故障检测滤波器的设计方法，但由于Lyapunov矩阵与系统矩阵之间存在耦合，所得的结果具有一定的保守性；其次，将考虑使用松弛变量的方法，实现Lyapunov矩阵与系统矩阵的解耦，进而得到具有较低保守性的设计结果.

引理1^[9] 给定常数0＜α＜1,μ>1，对于误差系统(3)，定义多Lyapunov函数V(e(k))={V_i(e(k)),i∈N}，如果其满足以下两个性质：

1) V_i(e(k+1))－V_i(e(k))＜－αV_i(e(k)),

2) V_i(e(k_l))＜μV_j(e(k_l)),

并且平均驻留时间满足：

${{\tau }_{\alpha }}\ge \tau _{\alpha }^{*}=\text{ceil}\left[ -\frac{\ln \mu }{\ln \left( 1-\alpha \right)} \right],$

(6)

其中，ceil(υ)表示取整函数，则误差系统(3)是全局指数稳定的.

注1 引理1使用多Lyapunov函数和平均驻留时间方法，研究了在平均驻留时间切换信号下，离散时间切换系统稳定的条件，它可以保证条件1)满足.基于引理1，可以得到如下的引理2.

引理2 给定常数0＜α＜1，μ>1，γ>0，β>0,当存在正定对称矩阵P_i和Q_i(i∈N)，使得矩阵不等式：

$\left[ \begin{matrix} -{{P}_{i}} & {{P}_{i}}{{{\tilde{A}}}_{i}} & {{P}_{i}}{{{\tilde{D}}}_{i1}} \\ * & -\left( 1-\alpha \right){{P}_{i}}+\tilde{C}_{i}^{\text{T}}{{C}_{i}} & \tilde{C}_{i}^{\text{T}}{{{\tilde{D}}}_{i2}} \\ * & * & -{{\gamma }^{2}}I+\tilde{D}_{i2}^{\text{T}}{{{\tilde{D}}}_{i2}} \\ \end{matrix} \right]＜0,$

(7)

${{P}_{i}}＜\mu {{P}_{j}},\left( i\in \text{N, }j\in \text{N,}i\ne j \right)$

(8)

$\left[ \begin{matrix} -{{Q}_{i}} & {{Q}_{i}}{{{\tilde{A}}}_{i}} & {{Q}_{i}}{{{\tilde{F}}}_{i1}} \\ * & -\left( 1-\alpha \right){{Q}_{i}}+\tilde{C}_{i}^{\text{T}}{{C}_{i}} & -\tilde{C}_{i}^{\text{T}}{{{\tilde{F}}}_{i2}} \\ * & * & {{\beta }^{2}}I-\tilde{F}_{i2}^{\text{T}}{{{\tilde{F}}}_{i2}} \\ \end{matrix} \right]＜0,$

(9)

${{Q}_{i}}＜\mu {{Q}_{j}},\left( i\in \text{N, }j\in \text{N,}i\ne j \right)$

(10)

成立，同时切换信号满足平均驻留时间式(6)，则误差系统(3)全局指数稳定，并且满足加权H_∞性能(4)和H_{_}性能指标(5).

证明 证明过程与文献^[14]中定理1类似，故此省略.

注2 由于引理2中的式(7)和式(9)中含有非线性项，为了将其转化为线性矩阵不等式，需要进行变量替换：

${{Z}_{i}}={{P}_{i}}{{L}_{i}},{{Y}_{i}}={{Q}_{i}}{{L}_{i}},$

则有

${{L}_{i}}=P_{i}^{-1}{{Z}_{i}},{{L}_{i}}=Q_{i}^{-1}{{Y}_{i}}.$

可知滤波器增益矩阵L_i的数值由两个公式确定，为解决该问题，则需要令P_i=Q_i.这个附加条件的引入，相当于用一个Lyapunov函数同时求解加权H_∞性能和H_{_}性能.显而易见，这样得到的结果具有一定的保守性.分析原因，不难发现是由于Lyapunov矩阵与系统矩阵的耦合，才导致问题的产生.因此，本文将通过引入松弛变量的方法，解除二者的耦合，进而得到保守性较低的结果.

定理1 给定常数0＜α＜1，μ>1，γ>0，β>0,当存在正定对称矩阵P_i，Q_i，M_i和矩阵X_i，G_i(i∈N)，使得线性矩阵不等式：

$\left[ \begin{matrix} {{P}_{i}}-{{G}_{i}}-G_{i}^{\text{T}} & G_{i}^{\text{T}}{{A}_{i}}-{{X}_{i}}{{C}_{i}} & G_{i}^{\text{T}}{{D}_{i1}}-{{X}_{i}}{{D}_{i2}} \\ * & -\left( 1-\alpha \right){{P}_{i}}+G_{i}^{\text{T}}{{M}_{i}}{{C}_{i}} & G_{i}^{\text{T}}{{M}_{i}}{{D}_{i2}} \\ * & * & -{{\gamma }^{2}}I+G_{i2}^{\text{T}}{{M}_{i}}{{D}_{i2}} \\ \end{matrix} \right]＜0,$

(11)

$\left[ \begin{matrix} {{Q}_{i}}-{{G}_{i}}-G_{i}^{\text{T}} & G_{i}^{\text{T}}{{A}_{i}}-{{X}_{i}}{{C}_{i}} & G_{i}^{\text{T}}{{F}_{i1}}-{{X}_{i}}{{F}_{i2}} \\ * & -\left( 1-\alpha \right){{Q}_{i}}-G_{i}^{\text{T}}{{M}_{i}}{{C}_{i}} & -G_{i}^{\text{T}}{{M}_{i}}{{F}_{i2}} \\ * & * & {{\beta }^{2}}I-F_{i2}^{\text{T}}{{M}_{i}}{{F}_{i2}} \\ \end{matrix} \right]＜0.$

(12)

及式(8)，(10)成立，且切换信号的平均驻留时间满足式(6)，则误差系统(3)全局指数稳定，并且满足加权H_∞性能(4)和H_{_}性能指标(5)，同时鲁棒故障检测滤波器的参数可根据

${{L}_{i}}={{\left( G_{i}^{\text{T}} \right)}^{-1}}{{X}_{i}},{{V}_{i}}={{\left( {{M}_{i}} \right)}^{1/2}}.$

(13)

求解得到.

证明 若式(11)成立，则P_i－G_i－G_i^T＜0，可得到G_i+G_i^T>P_i>0，说明G_i是可逆矩阵.因为P_i为正定对称矩阵，有

${{\left( {{P}_{i}}-{{G}_{i}} \right)}^{\text{T}}}P_{i}^{-1}\left( {{P}_{i}}-{{G}_{i}} \right)\ge 0,$

(14)

成立，将其展开可得:

$G_{i}^{\text{T}}P_{i}^{-1}{{G}_{i}}\ge {{G}_{i}}+G_{i}^{\text{T}}-{{P}_{i}}.$

(15)

在式(15)左右两端同时乘以－I，则有

$-G_{i}^{\text{T}}P_{i}^{-1}{{G}_{i}}\le {{P}_{i}}-G_{i}^{\text{T}}-{{G}_{i}},$

(16)

则可知矩阵不等式

$\left[ \begin{matrix} -G_{i}^{\text{T}}P_{i}^{-1}{{G}_{i}} & G_{i}^{\text{T}}{{A}_{i}}-{{X}_{i}}{{C}_{i}} & G_{i}^{\text{T}}{{D}_{i1}}-{{X}_{i}}{{D}_{i2}} \\ * & -\left( 1-\alpha \right){{P}_{i}}+G_{i}^{\text{T}}{{M}_{i}}{{C}_{i}} & C_{i}^{\text{T}}{{M}_{i}}{{D}_{i2}} \\ * & * & -{{\gamma }^{2}}I+D_{i2}^{\text{T}}{{M}_{i}}{{D}_{i2}} \\ \end{matrix} \right]＜0,$

(17)

成立，再将X_i=G_i^TL_i带入，又由于G_i可逆，则在不等式(17)两端分别左乘T:=diag{P_i^T(G_i^－1)^T,I,I}和右乘矩阵T^T，得

$\left[ \begin{matrix} -{{P}_{i}} & {{P}_{i}}\left( {{A}_{i}}-{{L}_{i}}{{C}_{i}} \right) & {{P}_{i}}\left( {{D}_{i1}}-{{L}_{i}}{{D}_{i2}} \right) \\ * & -\left( 1-\alpha \right){{P}_{i}}+C_{i}^{\text{T}}{{M}_{i}}{{C}_{i}} & C_{i}^{\text{T}}{{M}_{i}}{{D}_{i2}} \\ * & * & -{{\gamma }^{2}}I+D_{i2}^{\text{T}}{{M}_{i}}{{D}_{i2}} \\ \end{matrix} \right]＜0.$

(18)

根据误差系统参数矩阵的定义，可由式(18)得到式(7)成立.

同理，对式(12)进行相似的变换可得到式(9)成立.

再根据式(8)、(10)，以及切换信号满足的条件(6)，则根据引理2，可以证明误差系统(3)全局指数稳定，并且满足加权H_∞性能(4)和H_{_}性能指标(5)，鲁棒故障检测滤波器的参数可由式(13)计算.至此，定理可证.证毕.

注3 定理1通过引入松弛变量G_i实现了Lyapunov矩阵与系统矩阵的解耦，可以在P_i≠Q_i时求解滤波器增益矩阵L_i和残差的权值矩阵V_i，降低了结果的保守性.

通过定理1可以得到鲁棒残差生成器，但由于干扰信号的作用，即使没有故障出现，残差也不为零，因此需要对残差进行评估，以判断是否出现了故障.

首先，选取残差的范数作为评价函数：

${{J}_{L}}\left( r\left( k \right) \right)=\left( \sum\limits_{k=0}^{L}{{{r}^{\text{T}}}\left( k \right)r\left( k \right)} \right)_{d\left( k \right)\ne 0,f\left( k \right)=0.}^{1/2}$

其次，设置阈值为无故障时残差范数的最大值，即

${{J}_{L}}\left( r\left( k \right) \right)={{\left( \sup {{\left\| r\left( k \right) \right\|}_{2}} \right)}_{d\left( k \right)\ne 0,f\left( k \right)=0.}}$

将残差评价函数的值和阈值比较，当J_L(r(k))≥J_th(r(k))，说明系统存在故障，当J_L(r(k))＜J_th(r(k))，则系统无故障.

3 仿真举例

本文考虑具有两个子系统的离散时间线性切换系统，每个子系统的模型如式(1)所示，参数矩阵给定如下：

$\begin{align} & {{A}_{1}}=\left[ \begin{matrix} -0.9 & -1 & 0 \\ 2 & 0 & 1 \\ 1 & 1 & 0 \\ \end{matrix} \right],{{A}_{2}}=\left[ \begin{matrix} -0.8 & -1 & 0 \\ 2 & 0 & 1 \\ 1.01 & 1 & 0 \\ \end{matrix} \right]. \\ & {{B}_{11}}={{F}_{11}}={{\left[ \begin{matrix} 1 & 1 & 5 \\ \end{matrix} \right]}^{\text{T}}},{{B}_{21}}={{F}_{21}}={{\left[ \begin{matrix} 1 & 1 & -8 \\ \end{matrix} \right]}^{\text{T}}}. \\ & {{D}_{11}}={{\left[ \begin{matrix} -0.9 & -1 & 3 \\ \end{matrix} \right]}^{\text{T}}},{{D}_{21}}={{\left[ \begin{matrix} -1 & -1 & 4 \\ \end{matrix} \right]}^{\text{T}}}. \\ & {{C}_{1}}=\left[ \begin{matrix} -1 & 0.8 & 2 \\ 1 & -1.2 & 0 \\ \end{matrix} \right],{{C}_{2}}=\left[ \begin{matrix} -0.9 & -1 & 2 \\ 1 & -1 & 0 \\ \end{matrix} \right]. \\ & {{B}_{12}}={{F}_{12}}={{\left[ \begin{matrix} -10 & -100 \\ \end{matrix} \right]}^{\text{T}}}, \\ & {{B}_{22}}={{F}_{22}}={{\left[ \begin{matrix} -20 & -100 \\ \end{matrix} \right]}^{\text{T}}}, \\ & {{D}_{12}}={{D}_{22}}={{\left[ \begin{matrix} 1 & 1 \\ \end{matrix} \right]}^{\text{T}}}. \\ \end{align}$

给定α=0.2,μ=1.2，则切换信号需要满足的最小平均驻留时间计算为

$\tau _{\alpha }^{*}=\text{ceil}\left[ -\frac{\ln \mu }{\ln \left( 1-\alpha \right)} \right]=\text{ceil}\left[ -\frac{\ln 1.2}{\ln 0.8} \right]=1.$

故选定平均驻留时间为τ_α=4的切换信号，如图 1所示.输入信号u(k)=1,k∈[0, 300].干扰信号d(k)=0.1sin(0.1k),k∈[0, 300]，如图 2所示.

故障信号描述如下：

$f\left( k \right)=\left\{ \begin{matrix} 0, & k\in \left[ 1,150 \right]; \\ 0.5, & k\in \left[ 151,176 \right]; \\ 0, & k\in \left[ 177,300 \right]. \\ \end{matrix} \right.$

给定γ=0.2,β=3，使用定理1可以得到G_i，P_i和Q_i(i=1,2)为

$\begin{align} & {{G}_{1}}=\left[ \begin{matrix} 0.0132 & 0.0007 & 0.0116 \\ 0.0008 & 0.0027 & 0.0005 \\ 0.0119 & 0.0004 & 0.0115 \\ \end{matrix} \right], \\ & {{G}_{2}}=\left[ \begin{matrix} 0.0092 & -0.008 & 0.0070 \\ -0.0009 & 0.0017 & 0 \\ 0.0073 & 0 & 0.0064 \\ \end{matrix} \right]. \\ & {{P}_{1}}=\left[ \begin{matrix} 0.0092 & 0.0001 & 0.0072 \\ 0.0001 & 0.0026 & 0.0001 \\ 0.0072 & 0.0001 & 0.0070 \\ \end{matrix} \right], \\ & {{P}_{2}}=\left[ \begin{matrix} 0.0088 & -0.0005 & 0.0064 \\ 0.0064 & 0.0025 & 0.0003 \\ 0.0064 & 0.0003 & 0.0060 \\ \end{matrix} \right]. \\ & {{Q}_{1}}=\left[ \begin{matrix} 0.0113 & -0.0011 & 0.0099 \\ -0.0011 & 0.0027 & -0.0010 \\ 0.0099 & -0.0010 & 0.0100 \\ \end{matrix} \right], \\ & {{Q}_{2}}=\left[ \begin{matrix} 0.0111 & -0.0018 & 0.0090 \\ -0.0018 & 0.0027 & -0.0007 \\ 0.0090 & -0.0007 & 0.0091 \\ \end{matrix} \right]. \\ \end{align}$

可以看出P_i≠Q_i(i=1,2)，同时可计算得到滤波器增益矩阵L₁、L₂和残差的权值矩阵V₁、V₂为

$\begin{align} & {{L}_{1}}=\left[ \begin{matrix} 0.1739 & 0.0650 \\ -0.1910 & 0.0950 \\ -0.1519 & 0.0411 \\ \end{matrix} \right],{{L}_{2}}=\left[ \begin{matrix} 0.3215 & 0.0757 \\ -0.197.2 & 0.2290 \\ -0.3366 & 0.1693 \\ \end{matrix} \right]; \\ & {{V}_{1}}=\left[ \begin{matrix} 0.0176 & 0.0266 \\ 0.0266 & 0.0433 \\ \end{matrix} \right],{{V}_{2}}=\left[ \begin{matrix} 0.0075 & 0.0144 \\ 0.0144 & 0.0428 \\ \end{matrix} \right]. \\ \end{align}$

计算残差的评价函数和阈值为

$\begin{align} & {{J}_{L}}\left( r\left( k \right) \right)=\sqrt{\sum\limits_{k=1}^{161}{{{r}^{\text{T}}}\left( k \right)r\left( k \right)}}=15.9957, \\ & {{J}_{\text{th}}}={{\left( \sup {{\left\| r\left( k \right) \right\|}_{2}} \right)}_{f\left( k \right)=0}}=\sqrt{\sum\limits_{k=1}^{300}{{{r}^{\text{T}}}\left( k \right)r\left( k \right)}}=15.6597. \\ \end{align}$

可以看出在k=161时，J_L(r(k))>J_th(r(k))，检测到了故障，残差评价函数的变化曲线如图 3所示.

在图 3中，虚线表示没有故障的残差评价函数变化曲线，实线表示有故障的残差变化曲线.在故障发生前，两个曲线是重合的；当故障发生后，含有故障的残差评价曲线发生变化，明显地超出了无故障的曲线，说明检测出故障.

4 结 论

1) 将离散时间线性切换系统的故障检测，转化为了满足混合H_{_}/H_∞性能指标的故障检测滤波器设计问题，利用平均驻留时间和多Lyapunov函数方法设计了滤波器的参数.

2) 通过松弛变量方法解除了Lyapunov矩阵与系统矩阵的耦合问题，得到保守性较低的结果.

3)利用一个例子进行仿真验证，仿真结果表明本文设计的故障检测滤波器可以实现故障检测的功能.