基于新息约束的大规模系统的自适应事件触发故障估计

Adaptive-Event-Triggered Fault Estimation for Large-Scale Systems via Innovation Constraints

摘要: 本文研究了具有测量输出量化的大规模系统的故障估计问题，其中子估计器之间的通信拓扑与子系统的物理耦合结构相同。针对信道带宽有限的问题，提出了一种新的自适应事件触发方案来控制子估计器之间的数据传输。利用李雅普诺夫分析方法，得到了保证估计误差的动力学在均方意义下指数最终有界的充分条件。利用矩阵不等式技术，得到了期望滤波器增益的解析表达式和滤波误差的最终界。最后，通过一个数值例子验证了该设计方法的有效性。

Abstract: This paper focuses on the fault estimation issues of large-scale systems with measurements output, where the communication topology among sub-estimators is the same with the physical coupling structure of subsystems. Considering the limited channel bandwidth, a novel adaptive event-triggered scheme is proposed for governing the data transmission among sub-estimators. With the help of Lyapunov analysis approaches, sufficient conditions are derived to ensure that the dynamics of estimation errors is exponentially ultimately bounded in mean-square sense. Furthermore, an analytic formula of the desired filter gain and the ultimate bound of filtering errors are obtained through the utilization of matrix inequality techniques. Finally, one numerical example is given to testify the effectiveness of the proposed design approach.

1. 引言

大规模系统在工程实践中被广泛运用，包括配电系统、道路交通网络和传感器、执行器网络 [1]。在现有的事件触发方案中，通常采用一个恒定的阈值，该协议不能通过动态调整通信频率来提高系统性能。因此，自适应方案得到了越来越多的研究关注 [2] [3]。在过去的几十年里，由于现代工业对安全性和可靠性的要求越来越高，故障检测技术一直受到研究关注，已有很多文献进行了大量研究 [4]，故障估计技术是准确跟踪故障大小、形状等信息的关键环节 [5]。应该注意的是，量化方法一般分为均匀量化和非均匀量化。在 [6] 中，一个扇区边界方案被提出来处理反馈控制系统中的对数量化效应，对于含执行器饱和的系统的控制与滤波问题也已经有了一定的研究 [7]，而设计滤波器增益的波动可能对所开发的滤波器的动态行为有重大影响 [8] [9]。

在现有的关于大规模系统故障估计的研究中，每个子系统的子估计器需要在每个时刻与所有邻居进行通信，这将消耗大量的通信资源，而在如今的大规模网络化通讯系统中网络资源尤为珍贵，频繁的数据交互势必会出现网络的负载饱和，甚至造成网络拥塞，有限的网络资源不可避免地会导致网络诱导现象，从而使系统性能低下。事件触发协议的一个基本特性是，只有当某个事件被激活时，信息才被优先传输。与其他协议相比，自适应的事件触发协议易于在应用层中设计和执行。引入的自适应规则通常是单调降低的，因此事件触发阈值越小，信息传递就越频繁。因此，当系统动态趋于稳定时，信息传输的频率不能降低，故应该设计一个自适应方案来消除现有方案的单调性。此外，考虑到固有的物理耦合，为大规模网络系统开发一个可扩展的算法也具有挑战性。因此，本文的主要动机之一是提出一种新的自适应规则来服务于大规模网络系统的故障估计问题。同时信号量化得到了广泛的关注，并有大量的研究成果被报道。实际上，执行器/传感器的饱和可能是实际控制系统中遇到的最常见的非线性，如果在控制器/滤波器设计中忽略这种非线性，可能会降低系统的性能，甚至导致不稳定。滤波技术由于在控制和信号处理中提供一些关键信息的优良能力，在电力系统、全球定位系统、工业过程监控、目标跟踪等领域得到了越来越多的应用。在实际工程中，即使滤波策略已设计良好，设计滤波器增益的波动可能在其实现过程中发生。这些波动是由于数字信号处理中数值计算的舍入误差、滤波器中电子元件的参数偏差以及电路中的老化元件造成的。所有这些波动，尽管很小，但不应被忽略，我们应该设计对增益的波动或变化不敏感的滤波策略。因此，本文的研究很有意义。

本文的研究目标是利用自适应事件触发机制解决基于新息饱和的大规模系统的故障估计问题。本文的主要贡献有三个方面：1) 针对具有测量输出量化的大规模系统，提出了一种具有自适应机制的事件触发方案来控制子估计器之间的信息传输；2) 结合估计误差的界，得到了估计误差动力学在均方意义下指数最终有界的充分条件；3) 利用李雅普诺夫理论，建立了基于线性矩阵不等式的设计框架，以获得期望的增益。

2. 问题描述

我们考虑由N个耦合子系统组成的离散时间线性大规模系统：

$x_{i,k+1}=A_{ii}{x}_{i,k}+{\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{A}_{ij}{x}_{j,k}}+D_i{\omega }_{i,k}+{\alpha }_k{F}_i{f}_{i,k}$ (1)

输出为

$y_{i,k}=C_i{x}_{i,k}$ (2)

其中 $x_{i,k}\in {ℝ}^n$ 和 $y_{i,k}\in {ℝ}^m$ 分别是子系统i的状态向量和测量输出。 $A_{ii},A_{ij},D_i,F_i$ 和 $C_i$ 是具有适当维数的已知的常数矩阵。耦合结构由一个图表ϝ $=\left(\nu ,\tau \right)$ 来刻画，其中 $\nu =\left\{1,2,\cdots ,N\right\}$ 是有限子系统集， $\tau \subseteq \nu \times \nu$ 是代表内在物理联系的边界集。 ${\mathbb{S}}_i=\left\{j:\left(j,i\right)\in \tau ,i\in \nu \right\}$ 表示子系统的邻居集， ${\omega }_{i,k}\in {ℝ}^n$ 表示满足 $\mathbb{E}\left\{{\omega }_{i,k}\right\}=0$ 和 $\mathbb{E}\left\{{\Vert {\omega }_{i,k}\Vert }_2\right\}\le \bar{\omega }$ 的有界过程噪音， $f_{i,k}\in {ℝ}^{n_f}$ 是一个未知的故障信号。如果 $A_{ij}\ne 0$，子系统j被称作是子系统i的邻居子系统。此外，对于子系统i，它的邻居子系统的数量记为 $m_i$，它包含了子系统i。假设可以表示为如下形式的故障动力学：

$\begin{array}{r}\hfill f_{i,k+1}=A_{f_i}{f}_{i,k}\end{array}$ (3)

其中 $A_{f_i}$ 是一个合适维数的已知矩阵。容易看到当 $A_{f_i}\equiv I$ 时，故障是一个常数。随机变量 ${\alpha }_k\in \left\{0,1\right\}$，描述了随机发生故障的现象，满足伯努利分布，取值0或1

$Prob\left({\alpha }_k=1\right)=\mathbb{E}\left({\alpha }_k\right)=\bar{\alpha }Prob\left({\alpha }_k=0\right)=1-\bar{\alpha }$ (4)

其中 $\bar{\alpha }\in \left[0,1\right]$ 是一个已知标量。

在本文中，每个子系统通过使用对数静态和时不变量化器对测量输出 $y_{i,k}$ 进行量化。表示对称量化器为 $g(\cdot )={\left[g_1(\cdot )g_2(\cdot )\cdots g_m(\cdot )\right]}^{\text{T}}$，也就是说 $g_j\left(-v\right)=-g_j\left(v\right),j=1,\cdots ,m$。量化过程的映射为

${\bar{y}}_{i,k}=g\left(y_{i,k}\right)={\left[g_1\left(y_{i,k}^{\left(1\right)}\right)g_2\left(y_{i,k}^{\left(2\right)}\right)\cdots g_m\left(y_{i,k}^{\left(m\right)}\right)\right]}^{\text{T}}$ (5)

对每个 $g_j(\cdot )\left(1\le j\le m\right)$，量化等级的集合描述为

$U_j=\left\{\pm {\hat{\mu }}_l^{\left(j\right)},{\hat{\mu }}_l^{\left(j\right)}={\chi }_j^l{\hat{\mu }}_0^{\left(j\right)},l=0,\pm 1,\pm 2,\cdots \right\}\cup \left\{0\right\},0<{\chi }_j<1,{\hat{\mu }}_0^{\left(j\right)}>0$ (6)

其中 ${\chi }_j\left(j=1,\cdots ,m\right)$ 称为量化密度。每一个量化层对应一个部分，这样量化器将整个部分映射到这个量化层。给出了如下的对数量化器

$g_j\left(y_{i,k}^{\left(j\right)}\right)=\{\begin{array}{l}{\hat{\mu }}_l^{\left(j\right)},\frac{1}{1+{\beta }_j}{\hat{\mu }}_l^{\left(j\right)}\le y_{i,k}^{\left(j\right)}\le \frac{1}{1-{\beta }_j}{\hat{\mu }}_l^{\left(j\right)}\\ 0,y_{i,k}^{\left(j\right)}=0\\ -g_j\left(-y_{i,k}^{\left(j\right)}\right),y_{i,k}^{\left(j\right)}<0\end{array}$ (7)

其中 ${\beta }_j=\left(1-{\chi }_j\right)/\left(1+{\chi }_j\right)$。从上面的定义可以很容易地看出 $g_j\left(y_{i,k}^{\left(j\right)}\right)=\left(1+{\Delta }_k^{\left(j\right)}\right)y_{i,k}^{\left(j\right)}$，其中 $\left|{\Delta }_k^{\left(j\right)}\right|\le {\beta }_j$。根据上面讨论的变换，量化效应可以转化为扇形有界不确定性。定义 ${\Delta }_k=diag\left\{{\Delta }_k^{\left(1\right)},\cdots ,{\Delta }_k^{\left(m\right)}\right\}$，具有量化效应的测量可以表示为

$\begin{array}{r}\hfill {\bar{y}}_{i,k}=\left(I+{\Delta }_k\right)y_{i,k}=\left(I+{\Delta }_k\right)C_i{x}_{i,k}\end{array}$ (8)

通过定义   $\bar{\Delta }=diag\left\{{\beta }_1,\cdots ,{\beta }_m\right\}$， $E_k={\Delta }_k{\bar{\Delta }}^{-1}$。我们可以得到一个满足 $E_k{E}_k^{\text{T}}\le I$ 的未知实值时变矩阵。

让我们看看通信网络。对于大规模系统的故障估计问题，子估计器的估计误差必须被传输到它们的耦合邻居，因此，如果在每个瞬间都执行数据传输，通信负担将大大增加。为了降低传输频率，减轻通信负担，提出了一种自适应的事件触发方案，事件触发发生器函数 ${\zeta }_i\left(\cdot ,\cdot \right):{ℝ}^{n_f}\times ℝ\to ℝ$ 定义如下：

$\begin{array}{r}\hfill {\zeta }_i\left({\phi }_{i,k},{\delta }_{i,k}\right)={\phi }_{i,k}^{\text{T}}{\Psi }_i{\phi }_{i,k}-{\delta }_{i,k}\end{array}$ (9)

其中 ${\phi }_{i,k}=y_{i,k}^t-y_{i,k}$， $y_{i,k}^t$ 是最新的传输测量， $y_{i,k}$ 为当前测量值， ${\Psi }_i$ 是一个预先给定的正定对称权矩阵， ${\delta }_{i,k}$ 为待设计的时变正标量。只要满足不等式

$\begin{array}{r}\hfill {\zeta }_i\left({\phi }_{i,k},{\delta }_{i,k}\right)>0\end{array}$ (10)

就触发执行。

因此，事件触发的瞬间序列 $0\le k_0^i<k_1^i<\cdots <k_d^i<\cdots$ 可以递归地确定

$k_{d+1}^i=\min \left\{k|k>k_d^i,\begin{array}{r}\hfill {\zeta }_i\left({\phi }_{i,k},{\delta }_{i,k}\right)>0\end{array}\right\}$ (11)

为了实现自适应调节，设计了标量 ${\delta }_{i,k}$ 的自适应规则：

${\delta }_{i,k+1}=\min \left\{\max \left\{{\delta }_m,{\delta }_{i,k}\left(1-\frac{2a}{\pi }\arctan \left[\Vert {\phi }_{i,k}\Vert -b_i\right]\right)\right\},{\delta }_M\right\}$ (12)

其中 $\arctan (\cdot )$ 是反正切函数。 ${\delta }_M$ 和 ${\delta }_m$ (满足 ${\delta }_m\le {\delta }_M$ )是两个提前给定的标量， $0<a<1$ 和 $b_i>0$ 是两个用以调整 $\arctan (\cdot )$ 输出值的给定常数。此外，令 ${\delta }_{i,0}={\delta }_m$。根据条件(12)，我们可以很容易地得到下面的引理。

引理1：对于有给定初始条件 ${\delta }_{i,0}={\delta }_m>0$ 的自适应事件触发条件(11)，标量 ${\delta }_{i,k}$ 对所有 $k\in \mathbb{N}$ 都满足 ${\delta }_m\le {\delta }_{i,k}\le {\delta }_M$，条件 ${\phi }_{i,k}^{\text{T}}{\Psi }_i{\phi }_{i,k}\le {\delta }_M$ 总是成立。

设 ${\bar{x}}_{i,k}={\left[x_{i,k}^{\text{T}}{f}_{i,k}^{\text{T}}\right]}^{\text{T}}$，根据(1)~(3)，我们有

$\begin{array}{r}\hfill \{\begin{array}{l}{\bar{x}}_{i,k+1}=\left(\bar{A}+{\tilde{\alpha }}_k{\bar{F}}_i\right){\bar{x}}_{i,k}+{\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\bar{A}}_{i,j}{\bar{x}}_{j,k}}+{\bar{D}}_i{\omega }_{i,k}\\ {\bar{y}}_{i,k}=\left(I+{\Delta }_k\right)\left({\bar{C}}_i{\bar{x}}_{i,k}+{\phi }_{i,k}\right)\end{array}\end{array}$ (13)

其中

$\begin{array}{l}{\bar{A}}_i=\left[\begin{array}{cc}{A}_{ii}& \bar{\alpha }{F}_i\\ 0& A_{f_i}\end{array}\right],{\bar{A}}_{i,j}=\left[\begin{array}{cc}{A}_{ij}& 0\\ 0& 0\end{array}\right],{\bar{F}}_i=\left[\begin{array}{cc}0& F_i\\ 0& 0\end{array}\right]\\ {\bar{D}}_i=\left[\begin{array}{c}{D}_i\\ 0\end{array}\right],{\bar{C}}_i=\left[\begin{array}{cc}{C}_i& 0\end{array}\right],{\tilde{\alpha }}_k={\alpha }_k-\bar{\alpha }\end{array}$

在本文中，我们设计了如下的故障估计器：

${\hat{x}}_{i,k+1}={\bar{A}}_i{\hat{x}}_{i,k}+{\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\bar{A}}_{i,j}{\hat{x}}_{j,k}}+K_{ii}Sat\left({\bar{y}}_{i,k}-{\bar{C}}_i{\hat{x}}_{i,k}\right)$ (14)

其中 ${\hat{x}}_{i,k}$ 表示滤波器状态向量， $K_{ii}$ 是待确定的滤波器参数。为了实现性能分析和增益设计，饱和函数 $Sat(\cdot ):{ℝ}^{n_y}\to {ℝ}^{n_y}$ 可按预定的工作范围分解为一个线性部分和一个非线性部分之和。具体地说，我们有

$Sat\left(\rho \right)=G_1\rho +h\left(\rho \right)$ (15)

其中已知矩阵 $G_1$ 满足 $G_1\le I$，非线性向量值函数 $h\left(\rho \right)$ 具有下列扇区有界条件：

$h^{\text{T}}\left(\rho \right)\left(h\left(\rho \right)-G_2\rho \right)\le 0$ (16)

其中 $G_2\ge 0$ 是一个已知矩阵。定义局部估计误差 $e_{i,k}={\bar{x}}_{i,k}-{\hat{x}}_{i,k}$，可以得到

$\begin{array}{c}{e}_{i,k+1}={\bar{A}}_i{e}_{i,k}+{\tilde{\alpha }}_k{\bar{F}}_i{\bar{x}}_{i,k}+{\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\bar{A}}_{i,j}{e}_{j,k}}+{\bar{D}}_i{\omega }_{i,k}-K_{ii}\{G_1[{\bar{C}}_i{e}_{i,k}+{\Delta }_k{\bar{C}}_i{\bar{x}}_{i,k}\\ +\left(I+{\Delta }_k\right){\phi }_{i,k}]+h\left({\bar{C}}_i{e}_{i,k}+{\Delta }_k{\bar{C}}_i{\bar{x}}_{i,k}+\left(I+{\Delta }_k\right){\phi }_{i,k}\right)\}\\ =\left({\bar{A}}_i-K_{ii}{G}_1{\bar{C}}_i\right)e_{i,k}+{\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\bar{A}}_{i,j}{e}_{j,k}}+{\bar{D}}_i{\omega }_{i,k}+\left({\tilde{\alpha }}_k{\bar{F}}_i-K_{ii}{G}_1{\Delta }_k{\bar{C}}_i\right){\bar{x}}_{i,k}\\ -K_{ii}{G}_1\left(I+{\Delta }_k\right){\phi }_{i,k}-K_{ii}h\left({\varsigma }_{i,k}\right)\end{array}$ (17)

其中 ${\varsigma }_{i,k}\triangleq {\bar{C}}_i{e}_{i,k}+{\Delta }_k{\bar{C}}_i{\bar{x}}_{i,k}+\left(I+{\Delta }_k\right){\phi }_{i,k}$。

令 ${\chi }_{i,k}\triangleq {\left[{\bar{x}}_{i,k}^{\text{T}}{e}_{i,k}^{\text{T}}\right]}^{\text{T}}$，可以得到以下的动态系统：

${\chi }_{i,k+1}=\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right){\chi }_{i,k}+{\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\mathcal{A}}_{i,j}{\chi }_{j,k}}+\tilde{D}{\omega }_{i,k}+{\bar{B}}_1{\phi }_{i,k}+{\bar{K}}_{ii}h\left({\text{ς}}_{i,k}\right)$ (18)

其中

$\begin{array}{l}{\mathcal{A}}_{i,k}=\left[\begin{array}{cc}{\bar{A}}_i+{\tilde{\alpha }}_k{\bar{F}}_i& 0\\ {\tilde{\alpha }}_k{\bar{F}}_i& {\bar{A}}_i\end{array}\right],{\mathcal{B}}_{i,k}=\left[\begin{array}{cc}0& 0\\ -K_{ii}{G}_1{\Delta }_k{\bar{C}}_i& -K_{ii}{G}_1{\bar{C}}_i\end{array}\right],{\mathcal{A}}_{i,j}=\left[\begin{array}{cc}{\bar{A}}_{i,j}& 0\\ 0& {\bar{A}}_{i,j}\end{array}\right]\\ {\tilde{D}}_i=\left[\begin{array}{c}{\bar{D}}_i\\ {\bar{D}}_i\end{array}\right],{\bar{B}}_i=\left[\begin{array}{c}0\\ -K_{ii}{G}_1\left(I+{\Delta }_k\right)\end{array}\right],{\bar{K}}_{ii}=\left[\begin{array}{c}0\\ -K_{ii}\end{array}\right]\end{array}$

定义1：系统(18)在均方意义上是指数最终有界的，如果存在一个函数 ${\beta }_1$ 和一个正常数 ${\theta }_2$ 使

$\mathbb{E}\left\{{\displaystyle {\sum }_{i=1}^N{\Vert {\chi }_{i,k}\Vert }^2}\right\}\le {\beta }_1\left(k,\mathbb{E}\left\{{\displaystyle {\sum }_{i=1}^N{\Vert {\chi }_{i,0}\Vert }^2}\right\}\right)+{\theta }_2$ (19)

本文的目标是为具有测量输出量化和事件触发通信协议的离散时间大规模系统(1)设计一组估计器(14)。具体地说，我们关注于寻找合适的估计参数 $K_{ii}$，使得估计误差的动力学(17)在均方意义上是指数最终有界的。

引理2 (S-过程)：H，M，N和L是具有适当维数的实矩阵，若L和M满足 $L=L^{\text{T}}$ 和 $M{M}^{\text{T}}\le I$。然后有 $L+HMN+N^{\text{T}}{M}^{\text{T}}{H}^{\text{T}}<0$，当且仅当存在正标量 $\epsilon >0$ 使得 $L+{\epsilon }^{-1}H{H}^{\text{T}}+\epsilon N^{\text{T}}N<0$ 或等价地

$\left[\begin{array}{ccc}L& H& \epsilon N^{\text{T}}\\ H^{\text{T}}& -\epsilon I& 0\\ \epsilon N& 0& -\epsilon I\end{array}\right]<0$

3. 主要结果

定理1：对于给定的正标量 ${\delta }_m$， ${\delta }_M$，a， $b_i$， ${\epsilon }_s>0\left(s=0,1,2,\cdots ,6\right)$，还有矩阵 $K_{ii}$， ${\Delta }_k$，有自适应事件触发传输方案(9)~(12)的估计误差动力学(18)在均方意义上是指数最终有界的，如果存在正定矩阵 $P_i$，对于一切 $i,j\in \nu$，正标量 ${\lambda }_i\left(i\in \nu \right)$ 和一个正标量 $ϵ>0$，使得下面的不等式成立

$\begin{array}{r}\hfill \Gamma \triangleq \left[\begin{array}{ccc}{L}_{1i}& \ast & \ast \\ H_{1i}^{\text{T}}& -ϵI& \ast \\ ϵ{N^{\prime }}_i& 0& -ϵI\end{array}\right]<0\end{array}$ (20)

其中

$L_{1i}=\left[\begin{array}{cccc}{\kappa }_{2i}-P_i+{\epsilon }_{0}I& \ast & \ast & \ast \\ 0& -{\lambda }_1{\Psi }_i& \ast & \ast \\ {\lambda }_2{G}_2{C^{\prime }}_i& {\lambda }_2{G}_{2}I& -2{\lambda }_{2}I& \ast \\ \sqrt{o_1}{\mathcal{A}}_{i,k}& 0& \sqrt{o_4}{\bar{K}}_{ii}& -P_i^{-1}\end{array}\right]$

$\begin{array}{r}\hfill H_{1i}={\left[\begin{array}{cccc}0& 0& {\lambda }_2{G}_2^{\text{T}}& 0\\ 0& 0& 0& \sqrt{o_3}{\left({\bar{K}}_{ii}{\bar{G}}_1\right)}^{\text{T}}\\ 0& 0& {\lambda }_2{G}_2^{\text{T}}& 0\\ 0& 0& 0& \sqrt{o_1}{\left({\bar{K}}_{ii}{G^{\prime }}_1\right)}^{\text{T}}\end{array}\right]}^{\text{T}}\end{array}$

${N^{\prime }}_i=\left[\begin{array}{cccc}\bar{\Delta }{C^{″}}_i& 0& 0& 0\\ 0& {\bar{\Delta }}^{″}& 0& 0\\ 0& \bar{\Delta }& 0& 0\\ N_i& 0& 0& 0\end{array}\right]$

此外，如果上述不等式成立，则估计误差动力学(18)的界为

$\theta =\frac{{\rho }_0}{\left({\rho }_0-1\right){\min }_{i\in N}{\lambda }_{\min }\left(P_i\right)}{\theta }_2$ (21)

其中

${\theta }_1={\displaystyle {\sum }_{i=1}^N\mathbb{E}}\left\{tr\left\{K_{3i}\right\}\right\},{\theta }_2={\displaystyle {\sum }_{i=1}^N\mathbb{E}}\left\{{\theta }_1+{\lambda }_1{\delta }_{i,k}\right\}$

其中 ${\rho }_0>1$，有

${\Psi }_1\left(\rho \right)=-{\epsilon }_0\rho +\left(\rho -1\right){\min }_{i\in N}{\lambda }_{\max }\left(P_i\right)$ (22)

证明：构造下面的李雅普诺夫函数

$V_k={\displaystyle {\sum }_{i=1}^N{\chi }_{i,k}^{\text{T}}{P}_i{\chi }_{i,k}}$ (23)

然后沿着轨迹(18)计算 $V_k$ 的差分，并取数学期望

$\begin{array}{c}\mathbb{E}\left\{\Delta V_k\right\}=\mathbb{E}\left\{V_{k+1}-V_k\right\}={\displaystyle {\sum }_{i=1}^N\mathbb{E}}\left\{{\chi }_{i,k+1}^{\text{T}}{P}_i{\chi }_{i,k+1}-{\chi }_{i,k}^{\text{T}}{P}_i{\chi }_{i,k}\right\}\\ ={\displaystyle {\sum }_{i=1}^N\mathbb{E}}\{{\left[\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right){\chi }_{i,k}+{\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\mathcal{A}}_{i,j}{\chi }_{j,k}}+\tilde{D}{\omega }_{i,k}+{\bar{B}}_i{\phi }_{i,k}+\bar{K}h\left({\varsigma }_{i,k}\right)\right]}^{\text{T}}{P}_i[\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right){\chi }_{i,k}\\ +{\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\mathcal{A}}_{i,j}{\chi }_{j,k}}+\tilde{D}{\omega }_{i,k}+{\bar{B}}_i{\phi }_{i,k}+\bar{K}h\left({\varsigma }_{i,k}\right)]-{\chi }_{i,k}^{\text{T}}{P}_i{\chi }_{i,k}\}\\ \le {\displaystyle {\sum }_{i=1}^N\mathbb{E}}\{{\chi }_{i,k}^{\text{T}}{\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right)}^{\text{T}}{P}_i\left[\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right){\chi }_{i,k}+2{\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\mathcal{A}}_{i,j}{\chi }_{j,k}}+2{\bar{B}}_i{\phi }_{i,k}+2\bar{K}h\left({\varsigma }_{i,k}\right)\right]\\ +{\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\chi }_{j,k}^{\text{T}}{\mathcal{A}}_{i,j}^{\text{T}}{P}_i\left[{\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\mathcal{A}}_{i,j}{\chi }_{j,k}}+2{\bar{B}}_i{\phi }_{i,k}+2\bar{K}h\left({\varsigma }_{i,k}\right)\right]}+{\phi }_{i,k}^{\text{T}}{\bar{B}}_i^{\text{T}}{P}_i\left[{\bar{B}}_i{\phi }_{i,k}+2\bar{K}h\left({\varsigma }_{i,k}\right)\right]\\ +h^{\text{T}}\left({\varsigma }_{i,k}\right){\bar{K}}^{\text{T}}{P}_i\bar{K}h\left({\varsigma }_{i,k}\right)-{\chi }_{i,k}^{\text{T}}{P}_i{\chi }_{i,k}+tr\left\{{\tilde{D}}^{\text{T}}{P}_i\tilde{D}\bar{\omega }\right\}\}\end{array}$ (24)

基于基本矩阵不等式，我们有

$\begin{array}{l}2{\chi }_{i,k}^{\text{T}}{\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right)}^{\text{T}}{P}_i{\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\mathcal{A}}_{i,j}{\chi }_{j,k}}\\ \le {\epsilon }_1{\chi }_{i,k}^{\text{T}}{\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right)}^{\text{T}}{P}_i\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right){\chi }_{i,k}+{\epsilon }_1^{-1}\left({\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\chi }_{j,k}^{\text{T}}{\mathcal{A}}_{i,j}^{\text{T}}}\right)P_i\left({\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\mathcal{A}}_{i,j}{\chi }_{j,k}}\right),\end{array}$ (25)

$2{\chi }_{i,k}^{\text{T}}{\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right)}^{\text{T}}{P}_i{\bar{B}}_i{\phi }_{i,k}\le {\epsilon }_2{\chi }_{i,k}^{\text{T}}{\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right)}^{\text{T}}{P}_i\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right){\chi }_{i,k}+{\epsilon }_2^{-1}{\phi }_{i,k}^{\text{T}}{\bar{B}}_i^{\text{T}}{P}_i{\bar{B}}_i{\phi }_{i,k},$ (26)

$2{\chi }_{i,k}^{\text{T}}{\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right)}^{\text{T}}{P}_i{\bar{K}}_{ii}h\left({\varsigma }_{i,k}\right)\le {\epsilon }_3{\chi }_{i,k}^{\text{T}}{\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right)}^{\text{T}}{P}_i\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right){\chi }_{i,k}+{\epsilon }_3^{-1}{h}^{\text{T}}\left({\varsigma }_{i,k}\right){\bar{K}}_{ii}^{\text{T}}{P}_i{\bar{K}}_{ii}h\left({\varsigma }_{i,k}\right),$ (27)

$2\left({\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\chi }_{j,k}^{\text{T}}{\mathcal{A}}_{i,j}^{\text{T}}}\right)P_i{\bar{B}}_i{\phi }_{i,k}\le {\epsilon }_4\left({\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\chi }_{j,k}^{\text{T}}{\mathcal{A}}_{i,j}^{\text{T}}}\right)P_i\left({\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\mathcal{A}}_{i,j}{\chi }_{j,k}}\right)+{\epsilon }_4^{-1}{\phi }_{i,k}^{\text{T}}{\bar{B}}_i^{\text{T}}{P}_i{\bar{B}}_i{\phi }_{i,k},$ (28)

$2\left({\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\chi }_{j,k}^{\text{T}}{\mathcal{A}}_{i,j}^{\text{T}}}\right)P_i{\bar{K}}_{ii}h\left({\varsigma }_{i,k}\right)\le {\epsilon }_5\left({\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\chi }_{j,k}^{\text{T}}{\mathcal{A}}_{i,j}^{\text{T}}}\right)P_i\left({\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\mathcal{A}}_{i,j}{\chi }_{j,k}}\right)+{\epsilon }_5^{-1}{h}^{\text{T}}\left({\varsigma }_{i,k}\right){\bar{K}}_{ii}^{\text{T}}{P}_i{\bar{K}}_{ii}h\left({\varsigma }_{i,k}\right)$, (29)

$2{\phi }_{i,k}^{\text{T}}{\bar{B}}_i^{\text{T}}{P}_i{\bar{K}}_{ii}h\left({\varsigma }_{i,k}\right)\le {\epsilon }_6{\phi }_{i,k}^{\text{T}}{\bar{B}}_i^{\text{T}}{P}_i{\bar{B}}_i{\phi }_{i,k}+{\epsilon }_6^{-1}{h}^{\text{T}}\left({\varsigma }_{i,k}\right){\bar{K}}_{ii}^{\text{T}}{P}_i{\bar{K}}_{ii}h\left({\varsigma }_{i,k}\right)$. (30)

现在，把(25)~(30)代入(24)可以得到

$\begin{array}{c}\mathbb{E}\left\{\Delta V_k\right\}\le {\displaystyle {\sum }_{i=1}^N\mathbb{E}}\{o_1{\left[\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right){\chi }_{i,k}\right]}^{\text{T}}{P}_i\left[\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right){\chi }_{i,k}\right]+o_2{\left({\displaystyle {\sum }_{j\in {\Gamma }_i}{\mathcal{A}}_{i,j}{\chi }_{j,k}}\right)}^{\text{T}}{P}_i\left({\displaystyle {\sum }_{j\in {\Gamma }_i}{\mathcal{A}}_{i,j}{\chi }_{j,k}}\right)\\ +o_3{\left({\bar{B}}_i{\phi }_{i,k}\right)}^{\text{T}}{P}_i\left({\bar{B}}_i{\phi }_{i,k}\right)+o_4{\left({\bar{K}}_{ii}h\left({\varsigma }_{i,k}\right)\right)}^{\text{T}}{P}_i\left({\bar{K}}_{ii}h\left({\varsigma }_{i,k}\right)\right)-{\chi }_{i,k}^{\text{T}}{P}_i{\chi }_{i,k}+tr\left\{{\tilde{D}}^{\text{T}}{P}_i\tilde{D}\bar{\omega }\right\}\}\end{array}$ (31)

其中

$\begin{array}{l}{o}_1=1+{\epsilon }_1+{\epsilon }_2+{\epsilon }_3,o_2=1+{\epsilon }_1^{-1}+{\epsilon }_4+{\epsilon }_5\\ o_3=1+{\epsilon }_2^{-1}+{\epsilon }_4^{-1}+{\epsilon }_6,o_4=1+{\epsilon }_3^{-1}+{\epsilon }_5^{-1}+{\epsilon }_6^{-1}\end{array}$

很容易得到

$\begin{array}{l}{\displaystyle {\sum }_{i=1}^N\mathbb{E}}\left\{{\left({\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\mathcal{A}}_{i,j}{\chi }_{j,k}}\right)}^{\text{T}}{P}_i\left({\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\mathcal{A}}_{i,j}{\chi }_{j,k}}\right)\right\}\\ \le {\displaystyle {\sum }_{i=1}^N\mathbb{E}}\left\{m_i{\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\chi }_{j,k}^{\text{T}}{\mathcal{A}}_{i,j}^{\text{T}}{P}_i{\mathcal{A}}_{i,j}{\chi }_{j,k}}\right\}\\ ={\displaystyle {\sum }_{i=1}^N\mathbb{E}}\left\{{\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\chi }_{i,k}^{\text{T}}{\mathcal{A}}_{j,i}^{\text{T}}{m}_j{P}_i{\mathcal{A}}_{j,i}{\chi }_{i,k}}\right\}\end{array}$ (32)

因此，令 ${\xi }_{i,k}={\left[{\chi }_{i,k}^{\text{T}}{\phi }_{i,k}^{\text{T}}{h}^{\text{T}}\left({\varsigma }_{i,k}\right)\right]}^{\text{T}}$，从(31)~(32)可以得出

$\begin{array}{c}\mathbb{E}\left\{\Delta V_k\right\}\le {\displaystyle {\sum }_{i=1}^N\mathbb{E}}\left\{{\chi }_{i,k}^{\text{T}}\left({\kappa }_{1i}+{\kappa }_{2i}-P_i\right){\chi }_{i,k}+o_3{\phi }_{i,k}^{\text{T}}{\bar{B}}_i^{\text{T}}{P}_i{\bar{B}}_i{\phi }_{i,k}+o_4{h}^{\text{T}}\left({\varsigma }_{i,k}\right){\bar{K}}_{ii}^{\text{T}}{P}_i{\bar{K}}_{ii}h\left({\varsigma }_{i,k}\right)+tr\left\{{\kappa }_{3i}\right\}\right\}\\ \le {\displaystyle {\sum }_{i=1}^N\mathbb{E}}\left\{{\xi }_{i,k}^{\text{T}}{\Gamma }_1{\xi }_{i,k}+{\theta }_1\right\}\end{array}$ (33)

其中

${\Gamma }_1=\left[\begin{array}{ccc}{\kappa }_{1i}+{\kappa }_{2i}-P_i& 0& 0\\ 0& {\kappa }_{5i}& 0\\ 0& 0& {\kappa }_{4i}\end{array}\right]$ (34)

$\begin{array}{l}{\kappa }_{1i}=o_1{\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right)}^{\text{T}}{P}_i\left({\mathcal{A}}_{i,k}+{\mathcal{B}}_{i,k}\right),{\kappa }_{2i}=o_2{\displaystyle {\sum }_{j\in {\mathbb{S}}_i}{\mathcal{A}}_{j,i}^{\text{T}}{m}_j{P}_i{\mathcal{A}}_{j,i}}\\ {\kappa }_{3i}={\tilde{D}}^{\text{T}}{P}_i\tilde{D}\bar{\omega },{\kappa }_{4i}=o_4{\bar{K}}_{ii}^{\text{T}}{P}_i{\bar{K}}_{ii},{\kappa }_{5i}=o_3{\bar{B}}_i^{\text{T}}{P}_i{\bar{B}}_i,{\theta }_1={\displaystyle {\sum }_{i=1}^N\mathbb{E}\left\{tr\left\{{\kappa }_{3i}\right\}\right\}}\end{array}$

针对自适应事件触发机制和扇形有界条件(16)，引入两个正标量 ${\lambda }_1,{\lambda }_2$，有

$\begin{array}{c}\mathbb{E}\left\{\Delta V_k\right\}\le {\displaystyle {\sum }_{i=1}^N\mathbb{E}}\left\{{\xi }_{i,k}^{\text{T}}{\Gamma }_1{\xi }_{i,k}-{\lambda }_1\left({\phi }_{i,k}^{\text{T}}{\Psi }_i{\phi }_{i,k}-{\delta }_{i,k}\right)-2{\lambda }_2{h}^{\text{T}}\left({\text{ς}}_{i,k}\right)\left(h\left({\varsigma }_{i,k}\right)-G_2{\varsigma }_{i,k}\right)+{\theta }_1\right\}\\ \le {\displaystyle {\sum }_{i=1}^N\mathbb{E}\left\{{\xi }_{i,k}^{\text{T}}{\Gamma }_2{\xi }_{i,k}\right\}}+{\theta }_2\end{array}$ (35)

其中

$\begin{array}{l}{\Gamma }_2=\left[\begin{array}{ccc}{\kappa }_{1i}+{\kappa }_{2i}-P_i& \ast & \ast \\ 0& {\kappa }_{5i}-{\lambda }_1{\Psi }_i& \ast \\ {\lambda }_2{G}_2\left({C^{\prime }}_i+{\Delta }_k{C^{″}}_i\right)& {\lambda }_2{G}_2\left(I+{\Delta }_k\right)& {\kappa }_{4i}-2{\lambda }_{2}I\end{array}\right]=\Gamma -diag\left\{{\epsilon }_{0}I,0,0\right\}\\ {\theta }_2={\displaystyle {\sum }_{i=1}^N\mathbb{E}}\left\{{\theta }_1+{\lambda }_1{\delta }_{i,k}\right\},{C^{″}}_i=\left[\begin{array}{cc}{\bar{C}}_i& 0\end{array}\right]\end{array}$

换句话说，

$\Gamma =\left[\begin{array}{ccc}{\kappa }_{1i}+{\kappa }_{2i}-P_i+{\epsilon }_{0}I& \ast & \ast \\ 0& {\kappa }_{5i}-{\lambda }_1{\Psi }_i& \ast \\ {\lambda }_2{G}_2\left({C^{\prime }}_i+{\Delta }_k{C^{″}}_i\right)& {\lambda }_2{G}_2\left(I+{\Delta }_k\right)& {\kappa }_{4i}-2{\lambda }_{2}I\end{array}\right]<0$ (36)

不难得到

${\mathcal{B}}_{i,k}={\bar{K}}_{ii}{G^{\prime }}_1{M}_1{N}_i,{\bar{B}}_i={\bar{K}}_{ii}{\bar{G}}_1{M}_1{\bar{\Delta }}^{″},{\bar{\Delta }}_k={\bar{\Delta }}^{\prime }{M}_{1}I$

其中

${G^{\prime }}_1=\left[\begin{array}{cc}{G}_1& G_1{\bar{C}}_i\end{array}\right],M_1=diag\left\{E_k,I\right\},N_i=\left[\begin{array}{cc}\overline{\stackrel{¯}{\Delta }}{\bar{C}}_i& 0\\ 0& I\end{array}\right]$

${\bar{G}}_1=\left[\begin{array}{cc}{G}_1& G_1\end{array}\right],{\bar{\Delta }}^{″}=\left[\begin{array}{c}\bar{\Delta }\\ I\end{array}\right]$

利用Schur引理，很容易知道 $\Gamma <0$ 等同于

$L_{1i}+H_{1i}{M}_2{N^{\prime }}_i+{\left(H_{1i}{M}_2{N^{\prime }}_i\right)}^{\text{T}}<0$ (37)

其中

$L_{1i}=\left[\begin{array}{cccc}{\kappa }_{2i}-P_i+{\epsilon }_{0}I& \ast & \ast & \ast \\ 0& -{\lambda }_1{\Psi }_i& \ast & \ast \\ {\lambda }_2{G}_2{C^{\prime }}_i& {\lambda }_2{G}_{2}I& -2{\lambda }_{2}I& \ast \\ \sqrt{o_1}{\mathcal{A}}_{i,k}& 0& \sqrt{o_4}{\bar{K}}_{ii}& -P_i^{-1}\end{array}\right],H_{1i}={\left[\begin{array}{cccc}0& 0& {\lambda }_2{G}_2^{\text{T}}& 0\\ 0& 0& 0& \sqrt{o_3}{\left({\bar{K}}_{ii}{\bar{G}}_1\right)}^{\text{T}}\\ {\lambda }_2{G}_2& 0& {\lambda }_2{G}_2^{\text{T}}& 0\\ 0& 0& 0& \sqrt{o_1}{\left({\bar{K}}_{ii}{G^{\prime }}_1\right)}^{\text{T}}\end{array}\right]}^{\text{T}}$

${N^{\prime }}_i=\left[\begin{array}{cccc}\bar{\Delta }{C^{″}}_i& 0& 0& 0\\ 0& {\bar{\Delta }}^{″}& 0& 0\\ 0& \bar{\Delta }& 0& 0\\ N_i& 0& 0& 0\end{array}\right],M_2=diag\left\{E_k,M_1,E_k,M_1\right\}$

通过引理2，由(37)可以得到

$L_{1i}+{ϵ}^{-1}{H}_{1i}{H}_{1i}^{\text{T}}+ϵ{\left({N^{\prime }}_i\right)}^{\text{T}}{N^{\prime }}_i<0$

也就是说

$\begin{array}{r}\hfill \Gamma =\left[\begin{array}{ccc}{L}_{1i}& \ast & \ast \\ H_{1i}^{\text{T}}& -ϵI& \ast \\ ϵ{N^{\prime }}_i& 0& -ϵI\end{array}\right]<0\end{array}$

现在，我们来揭示估计误差动力学的界。因为 $\Gamma <0$，不难由(20)和(36)得到下面的不等式成立：

$\mathbb{E}\left\{\Delta V_k\right\}\le -{\epsilon }_0\mathbb{E}\left\{{\displaystyle {\sum }_{i=1}^N{\Vert {\chi }_{i,k}\Vert }^2}\right\}+{\theta }_2$ (38)

然后，很容易得到

${\min }_{i\in N}{\lambda }_{\min }\left(P_i\right)\mathbb{E}\left\{{\displaystyle {\sum }_{i=1}^N{\Vert {\chi }_{i,k}\Vert }^2}\right\}\le \mathbb{E}\left\{V_k\right\}\le {\max }_{i\in N}{\lambda }_{\max }\left(P_i\right)\mathbb{E}\left\{{\displaystyle {\sum }_{i=1}^N{\Vert {\chi }_{i,k}\Vert }^2}\right\}$ (39)

在下面引入一个标量 $\varrho >1$，并结合李亚普诺夫函数(23)，我们可以计算

$\begin{array}{c}\mathbb{E}\left\{{\varrho }^{k+1}{V}_{k+1}-{\varrho }^k{V}_k\right\}={\varrho }^{k+1}\mathbb{E}\left\{V_{k+1}-V_k\right\}+{\varrho }^k\left(\varrho -1\right)\mathbb{E}\left\{V_k\right\}\\ \le -{\epsilon }_0{\varrho }^{k+1}\mathbb{E}\left\{{\displaystyle {\sum }_{i=1}^N{\Vert {\chi }_{i,k}\Vert }^2}\right\}+{\varrho }^{k+1}{\theta }_2+{\varrho }^k\left(\varrho -1\right){\max }_{i\in N}{\lambda }_{\max }\left(P_i\right)\mathbb{E}\left\{{\displaystyle {\sum }_{i=1}^N{\Vert {\chi }_{i,k}\Vert }^2}\right\}\\ \le {\psi }_1\left(\varrho \right){\varrho }^k\mathbb{E}\left\{{\displaystyle {\sum }_{i=1}^N{\Vert {\chi }_{i,k}\Vert }^2}\right\}+{\varrho }^{k+1}{\theta }_2\end{array}$ (40)

其中

${\psi }_1\left(\varrho \right)=-{\epsilon }_0\varrho +\left(\varrho -1\right){\max }_{i\in N}{\lambda }_{\max }(\; P\; i\; )$

把(39)的两边从0到s-1对k求和得到

$\mathbb{E}\left\{{\varrho }^s{V}_s\right\}-\mathbb{E}\left\{V_0\right\}\le {\psi }_1\left(\varrho \right){\displaystyle {\sum }_{k=0}^{s-1}{\varrho }^k\mathbb{E}\left\{{\displaystyle {\sum }_{i=1}^N{\Vert {\chi }_{i,k}\Vert }^2}\right\}}+\frac{\varrho \left(1-{\varrho }^s\right)}{1-\varrho }{\theta }_2$ (41)

我们可以看到存在一个标量 ${\varrho }_0=\frac{{\max }_{i\in N}{\lambda }_{\max }\left(P_i\right)}{{\max }_{i\in N}{\lambda }_{\max }\left(P_i\right)-{\epsilon }_0}$，使得 ${\psi }_1\left({\varrho }_0\right)=0$。那么，由(39)可知：

$\mathbb{E}\left\{{\varrho }^s{V}_s\right\}-\mathbb{E}\left\{V_0\right\}\le \frac{\varrho \left(1-{\varrho }^s\right)}{1-\varrho }{\theta }_2$ (42)

另一方面，很容易发现

$\mathbb{E}\left\{V_0\right\}\le {\max }_{i\in N}{\lambda }_{\max }\left(P_i\right)\mathbb{E}\left\{{\displaystyle {\sum }_{i=1}^N{\Vert {\chi }_{i,k}\Vert }^2}\right\}$ (43)

和

$\left\{{\varrho }^k{V}_k\right\}\ge {\varrho }^k{\min }_{i\in N}{\lambda }_{\min }\left(P_i\right)\mathbb{E}\left\{{\displaystyle {\sum }_{i=1}^N{\Vert {\chi }_{i,k}\Vert }^2}\right\}$ (44)

把这两个条件考虑进去，可以得到

$\mathbb{E}\left\{{\displaystyle {\sum }_{i=1}^N{\Vert {\chi }_{i,k}\Vert }^2}\right\}\le \frac{{\max }_{i\in N}{\lambda }_{\max }\left(P_i\right)}{{\varrho }_0^k{\min }_{i\in N}{\lambda }_{\min }\left(P_i\right)}\mathbb{E}\left\{{\displaystyle {\sum }_{i=1}^N{\Vert {\chi }_{i,0}\Vert }^2}\right\}+\frac{{\varrho }_0^k-1}{{\varrho }_0^{k-1}\left({\varrho }_0-1\right){\min }_{i\in N}{\lambda }_{\min }\left(P_i\right)}{\theta }_2$

根据定义1，我们可以得出系统(18)是均方意义上指数最终有界的，其最终界为：

$\theta ={\lim }_{k\to +\infty }\frac{{\varrho }_0^k-1}{{\varrho }_0^{k-1}\left({\varrho }_0-1\right){\min }_{i\in N}{\lambda }_{\min }\left(P_i\right)}{\theta }_2=\frac{{\varrho }_0}{\left({\varrho }_0-1\right){\min }_{i\in N}{\lambda }_{\min }\left(P_i\right)}{\theta }_2$ (45)

这个定理的证明就完成了。

定理1关注估计器的性能分析，为有测量输出量化和事件触发通信协议的离散时间大规模系统(1)开发所需的估计器设计方案奠定了基础。因此，我们可以考虑有测量输出量化和事件触发通信协议的离散时间大规模系统(1)的估计器设计问题。    

现在我们准备提出滤波器增益 $K_{ii}$ 的一个可行设计方案。

定理2：自适应事件触发传输方案(9)~(12)下的估计误差动力学(18)是均方意义下指数最终有界的，如果存在正定矩阵 $P_i$，对所有 $i,j\in \nu$ 存在正标量 ${\lambda }_i\left(i=1,2\right)$ 和一个正标量 $ϵ>0$，满足下面的不等式

$\begin{array}{r}\hfill {\Gamma }^{\text{*}}=\left[\begin{array}{cccc}{L}_{2i}& \ast & \ast & \ast \\ H_{2i}^{\text{T}}& -ϵI& \ast & \ast \\ ϵ{N^{\prime }}_i& 0& -ϵI& \ast \\ {U^{\prime }}_{1i}& 0& 0& {U^{\prime }}_{2i}\end{array}\right]<0\end{array}$ (46)

其中

$\begin{array}{l}{L}_{2i}=\left[\begin{array}{cccc}-P_i+{\epsilon }_{0}I& \ast & \ast & \ast \\ 0& -{\lambda }_1{\Psi }_i& \ast & \ast \\ {\lambda }_2{P}_i{G}_2{C^{\prime }}_i& {\lambda }_2{G}_{2}I& -2{\lambda }_{2}I& \ast \\ \sqrt{o_1}{P}_i{\mathcal{A}}_{i,k}& 0& \sqrt{o_4}{W}_{ii}& -P_i^{\text{T}}\end{array}\right],\begin{array}{r}\hfill H_{2i}^{\text{T}}=\left[\begin{array}{cccc}0& 0& {\lambda }_2{G}_2^{\text{T}}& 0\\ 0& 0& 0& \sqrt{o_3}{\left(W_{ii}{\bar{G}}_1\right)}^{\text{T}}\\ 0& 0& {\lambda }_2{G}_2^{\text{T}}& 0\\ 0& 0& 0& \sqrt{o_1}{\left(W_{ii}{G^{\prime }}_1\right)}^{\text{T}}\end{array}\right]\end{array}\\ {U^{\prime }}_{1i}=\left[\begin{array}{cccc}\sqrt{o_2}{P}_i{A}_{1,i}& 0& \cdots & 0\\ \sqrt{o_2}{P}_i{A}_{2,i}& 0& \cdots & 0\\ ⋮& ⋮& \ddots & ⋮\\ \sqrt{o_2}{P}_i{A}_{N,i}& 0& \cdots & 0\end{array}\right],{U^{\prime }}_{2i}=\left[\begin{array}{cccc}-m_1^{-1}{P}_i^{\text{T}}& \ast & \ast & \ast \\ 0& -m_2^{-1}{P}_i^{\text{T}}& \ast & \ast \\ ⋮& ⋮& \ddots & \ast \\ 0& 0& \cdots & -m_N^{-1}{P}_i^{\text{T}}\end{array}\right]\end{array}$

其中 $W_{ii}=P_i{\bar{K}}_{ii}$。

此外，增益参数可以通过 $\begin{array}{r}\hfill {\bar{K}}_{ii}=P_i^{-1}{W}_{ii}\end{array}$ 得到。进一步，估计误差动力学(18)的上界可以由(21)得到。

证明：根据Schur引理，由(20)可得

$\begin{array}{r}\hfill \Gamma \triangleq \left[\begin{array}{cccc}{L}_1& \ast & \ast & \ast \\ H_1^{\text{T}}& -ϵI& \ast & \ast \\ ϵN_2& 0& -ϵI& \ast \\ U_1& 0& 0& U_2\end{array}\right]<0\end{array}$ (47)

其中 $U_{1i}=\left[\begin{array}{cccc}\sqrt{o_2}{A}_{1,i}& 0& \cdots & 0\\ \sqrt{o_2}{A}_{2,i}& 0& \cdots & 0\\ ⋮& ⋮& \ddots & ⋮\\ \sqrt{o_2}{A}_{N,i}& 0& \cdots & 0\end{array}\right]$， $U_{2i}=\left[\begin{array}{cccc}-{\left(m_1{P}_i\right)}^{-1}& \ast & \ast & \ast \\ 0& -{\left(m_2{P}_i\right)}^{-1}& \ast & \ast \\ ⋮& ⋮& \ddots & \ast \\ 0& 0& \cdots & -{\left(m_N{P}_i\right)}^{-1}\end{array}\right]$。

对(47)分别左乘右乘 $diag\left\{I,I,I,P_i,I,I,I,I,I,I,I,I,P_i,\cdots ,P_i\right\}$ 和其转置，可得

$\begin{array}{r}\hfill {\Gamma }^{\text{*}}=\left[\begin{array}{cccc}{L}_2& \ast & \ast & \ast \\ H_2^{\text{T}}& -ϵI& \ast & \ast \\ ϵN_2& 0& -ϵI& \ast \\ {U^{\prime }}_1& 0& 0& {U^{\prime }}_2\end{array}\right]<0\end{array}$ (48)

滤波器增益可由 ${\bar{K}}_{ii}=P_i^{-1}{W}_{ii}$ 得到，这就完成了证明。

4. 数值仿真实例

在本节中，给出一个数值例子来验证所提出方法的有效性。取目标系统的参数分别为：

$A_{11}=\left[\begin{array}{cc}0.45& 1\\ 0& -0.84\end{array}\right],A_{22}=\left[\begin{array}{cc}-0.45& 0.9\\ 0& -0.66\end{array}\right],A_{12}=\left[\begin{array}{cc}0.1& 0.2\\ 0& 0.2\end{array}\right],A_{21}=\left[\begin{array}{cc}0.1& 0.2\\ 0& -0.4\end{array}\right]$

$D_1=\left[\begin{array}{cc}0.01& -0.04\\ 0.02& -0.01\end{array}\right],D_2=\left[\begin{array}{cc}0.24& -0.36\\ 0.15& -0.5\end{array}\right],C_1=\left[\begin{array}{cc}1& 0.25\end{array}\right],C_2=\left[\begin{array}{cc}1& 0.2\end{array}\right],$ $F_1=\left[\begin{array}{c}0.8\\ 0.3\end{array}\right],F_2=\left[\begin{array}{c}0.4\\ 0.7\end{array}\right]$

对于每个子系统，状态初始值为 $x_{1,0}={\left[1.28,1.39\right]}^{\text{T}}$， $x_{2,0}={\left[1.09,1.19\right]}^{\text{T}}$， ${\hat{x}}_{1,0}={\hat{x}}_{2,0}=0$。人为饱和约束水平设为0.3，扇形有界条件的参数设为 $G_1=I$， $G_2=diag\left\{0.2,0.2\right\}$，自适应事件触发阈值和可调参数为 ${\delta }_m=0.15$， ${\delta }_M=0.3$， $a=0.35$， $b_1=0.66$， $b_2=0.69$， $b_3=0.83$， $b_4=1.23$。此外，事件触发的加权矩阵为 ${\Psi }_i=0.35I$， ${\Psi }_2=0.24I$， ${\Psi }_3=0.28I$， ${\Psi }_4=0.37I$。预先给定的参数为 ${\epsilon }_1=0.1$， ${\epsilon }_2=0.021$， ${\epsilon }_3=0.23$， ${\epsilon }_4={\epsilon }_5={\epsilon }_6=0.01$， $\mu =0.001$。

根据定理2，使用带有LMI工具箱的Matlab软件，可以得到期望估计增益为

${\bar{K}}_{11}=\left[\begin{array}{c}\begin{array}{c}\begin{array}{c}-0.00007\\ 0.00005\end{array}\\ -0.00001\end{array}\\ 0.000048\end{array}\right],{\bar{K}}_{22}=\left[\begin{array}{c}\begin{array}{c}\begin{array}{c}-0.02046\\ 0.04617\end{array}\\ -0.00511\end{array}\\ 0.03737\end{array}\right]$

实际故障及其故障估计曲线见图1和图2，估计误差图见图3，子系统事件触发图见图4。

Figure 1. Actual fault and its estimation

图1. 实际故障及其估计

Figure 2. Actual fault and its estimation

图2. 实际故障及其估计

Figure 3. The evolution of errors

图3. 估计误差

Figure 4. Subsystem event-triggered graph

图4. 子系统事件触发图

5. 总结

本文研究了带有测量输出量化和新息饱和的离散时间大规模系统的故障估计问题。对于大规模系统，本文利用自适应事件触发协议减少了子系统间不必要的信息传输。假定过程噪声是有界的，随机变量 ${\alpha }_k$ 描述了随机发生的故障且满足伯努利分布，针对饱和引入一个滤波器来估计故障的情况，利用李雅普诺夫函数和线性矩阵不等式方法获得了估计器增益的参数。最后通过一个算例验证了该方法的有效性。此外，在本文研究的基础上，下一步将考虑基于新息约束、欺骗攻击的大规模系统的事件触发故障检测问题的研究。