A Compressed Sensing Algorithm Based on Multi-scale Residual Reconstruction Network
-
摘要: 目前压缩感知系统利用少量测量值使用迭代优化算法重构图像.在重构过程中，迭代重构算法需要进行复杂的迭代运算和较长的重构时间.本文提出了多尺度残差网络结构，利用测量值通过网络重构出图像.网络中引入多尺度扩张卷积层用来提取图像中不同尺度的特征，利用这些特征信息重构高质量图像.最后，将网络的输出与测量值进行优化，使得重构图像在测量矩阵上的投影与测量值更加接近.实验结果表明，本文算法在重构质量和重构时间上均有明显优势.Abstract: In recent years, a small number of measurements and iterative optimization algorithms were exploited in compressed sensing to reconstruct images. In the process of reconstruction, most algorithms based on iteration for compressed sensing image reconstruction suffer from the complicatedly iterative computation and time-consuming. In this paper, we propose a novel multi-scale residual reconstruction network (MSRNet), and exploit the measurements to reconstruct images through the network. The multi-scale dilate convolution layer is introduced in the network to extract the feature of different scales in the image, and the feature information could improve the quality of reconstructed image. Finally, we exploit the output of the network and measurements to optimize our algorithm, so as to make the projection of the reconstructed image closer to the measurements. The experimental results show that the MSRNet requires less running time and has better performance in reconstruction quality than other compressed sensing algorithms.
-
Key words:
- Compressive sensing /
- convolution neural network /
- multi-scale convolution /
- dilate convolution
-
近年来, 随着计算机、通信技术的快速发展, 多智能体理论成为研究热点, 并且相关控制方法被广泛应用于无人机集群、无人船编队、车辆编队等领域[1-3]. 作为多智能体领域的基础研究, 多智能体系统的一致性问题一直备受关注. 多智能体系统的一致性是指智能体间通过信息交互实现系统状态或输出趋同. 在过去十几年间, 基于一致性理论的多智能体系统协同控制得到了广泛研究[4-6].
对于实际系统来说, 故障是在所难免的, 其中执行机构故障由于其普遍性和破坏性而被广泛研究[7-8]. 为了保证多智能体系统的安全可靠运行, 容错控制问题一直都是重要研究方向. 文献[9-10]所设计的容错控制器旨在保证执行器发生故障时系统的稳定. 为实现执行器故障补偿控制, 文献[11-12]通过构造观测器来估计故障参数, 并基于估计结果设计容错控制器. 文献[12]研究基于分布式观测器的领导−跟随多智能体系统容错控制问题. 随着自适应技术的发展, 基于自适应技术的容错控制得到了广泛应用. 通过在控制器中引入自适应参数增益, 系统具有良好的容错性. 文献[13]在控制器中引入自适应参数, 降低了执行器故障产生的负面影响, 保证了系统的稳定性. 文献[14]研究一类离散多智能体系统的容错跟踪控制问题, 基于自适应技术来补偿死区和执行器故障对控制性能的影响. 综上所述, 多智能体系统的容错控制已经得到广泛研究. 考虑到实际多智能体系统中计算、通信等硬件资源的有限性, 研究网络通信受限情况下多智能体系统的可靠控制问题显然更有意义.
众所周知, 数字控制系统主要考虑的是非连续信息传输情况, 网络通信环境下智能体通常也不是每时每刻获取邻居的状态信息. 近年来, 随着事件触发技术在控制领域的广泛应用和不断发展, 科研人员开始研究基于事件触发控制的多智能体系统一致性问题, 以减轻通信资源的占用[15-16]. 在早期的多智能体系统事件触发控制研究中, 相邻智能体仍然需要持续通信, 这无疑导致了通信带宽的浪费. 为解决这一问题, 文献[17-19]提出一种固定采样周期的通信机制, 智能体之间只在特定时刻进行通信. 但是固定时间周期采样通信有很大的局限性, 仍然会产生不必要的采样. 为克服上述问题, 文献[20]基于周期采样机制研究线性多智能体系统的事件触发控制问题, 该触发机制要求智能体之间的信息传递仅在离散采样时刻进行.
值得一提的是, 在实际应用过程中智能体的状态可能是未知的. 文献[21-22]研究基于观测器的容错控制问题, 控制器所需的状态信息从状态观测器获得, 避免对系统状态的依赖. 文献[23]基于事件触发技术设计全分布式观测器和控制器, 该控制器仅需获取邻居的状态以及自身状态的观测信息. 在此基础上, 一种更激进的事件触发通信协议被提出, 该协议只需要获取自身信息, 当满足触发条件时才获取相邻智能体的观测值. 在文献[24]中, 智能体的通信触发判定条件仅需要测量误差, 因此克服了连续通信的限制, 同时减少了控制器的更新次数. 文献[25]所设计的控制器不仅在自身达到事件触发阈值时更新, 而且在其邻居达到事件触发阈值时更新. 在此基础上, 文献[26]使用观测器对领导者的未知参数进行预测. 文献[27]针对有向生成树拓扑结构的多智能体系统, 提出一种基于自触发通信的自适应容错控制方案. 为进一步扩大触发间隔, 文献[28-29]提出动态事件触发的概念, 通过在触发函数中引入动态变量或自适应参数, 提高触发阈值[30]. 近年来, 动态事件触发方法被大量应用于多智能体系统[31-33], 并取得了很好的效果. 文献[33]在触发函数中加入一个正常量, 可以避免芝诺现象, 同时可以延长触发间隔, 但多智能体系统仅可以达成有限共识.
为使多智能体系统具有良好的容错性, 同时降低智能体之间的通信频率, 文献[34-36]研究基于事件触发的容错控制问题. 文献[34]研究基于自适应事件触发技术的多智能体系统可靠控制问题. 文献[35]基于观测器研究多智能体系统自适应事件触发容错同步问题. 文献[36]基于事件触发机制设计开环观测器, 减少了智能体间的通信次数, 同时使得多智能体系统具有良好的容错性. 文献[37]设计一种新的分布式鲁棒容错控制器, 解决了不连续通信条件下的多智能体系统容错控制问题. 文献[38]研究具有执行器失效故障的异构线性多智能体系统的事件触发容错一致性问题. 文献[39]针对严格反馈离散多智能体系统, 提出一种自适应多梯度递归强化学习事件触发容错控制方案. 尽管基于事件触发的多智能体的容错控制问题已经得到广泛研究, 但采用事件触发机制的目标是降低智能体的控制器更新频率, 忽略了对多智能体系统通信资源的优化. 此外, 基于事件触发机制的多智能体系统故障类型及其处理方式研究还不够深入. 考虑到已有研究的不足, 本文研究基于事件触发通信协议的多智能体系统可靠共识问题, 提出一种新的事件触发通信机制, 有效地减少智能体之间的通信次数, 证明所设计的事件触发通信机制可以有效地避免芝诺现象. 本文的创新点如下:
1) 本文设计的事件触发通信机制是用来决定智能体之间是否通信, 控制器的更新只是多智能体系统通信的结果. 此外, 为延长触发间隔, 同时保持理想的动态特性, 所设计的控制器在触发间隔内是时变的.
2) 与现有的动态事件触发机制相比, 本文提出一种新的触发函数, 可以在不影响多智能体系统最终一致性的情况下有效地降低智能体之间的通信频率.
3) 不同于已有的多智能体系统动态事件触发自适应可靠控制方法, 本文考虑单个智能体的每个执行器都可能发生不同的失效故障.
1. 预备知识
1.1 图论
图$ {\cal{G}} = ({\cal{V}},\;{\cal{E}}) $由顶点集$ {\cal{V}} $和边集$ {\cal{E}} $组成, 其中$ {\cal{V}} = \{1,\;2,\;\cdots ,\;N\} $表示节点, $ {\cal{E}}\subseteq {\cal{V}}\times {\cal{V}} $表示节点间的信息流. $ {\cal{A}} = [a_{im}]\in {\bf{R}}^{N\times N} $表示$ {\cal{G}} $的加权邻接矩阵. 如果$ a_{im} = 0 $, 意味着智能体$ i $无法从$ m $获取信息; 如果$ a_{im}>0 $, 表示智能体$ i $可以接收到$ m $的状态信息. $ L = [l_{im}]\in {\bf{R}}^{N\times N} $表示$ {\cal{G}} $的拉普拉斯矩阵. 若$ i = m $, 则$ l_{ii} = \sum_{m = 1,\;m\neq i}^N a_{im} $; 否则, $ l_{im} = -a_{im} $. 如果图之间任意两个相邻节点均可以相互交流信息, 则该图是无向的; 如果任意两个节点之间存在一条连接彼此的路径, 则认为该图是全连通的.
假设 1. 图$ {\cal{G}} $是全连通且无向的.
引理 1. 如果假设1成立, 拉普拉斯矩阵$ L $的最小特征值为0且其余的特征值非负, $ L $矩阵为半正定矩阵[34].
定义 1. 矩阵的谱范数或向量的欧几里得范数表示为$ \vert \vert \cdot \vert\vert $, $ \otimes $表示克罗内克积, $ \text{diag}\{\cdots \} $表示对角矩阵, $ (\cdots ) $表示列向量.
定义 2. 本文中的一致意味着任意两个智能体的状态$ \lim_{t\to\infty} \Vert x_i(t) -x_m(t)\Vert = 0 $.
1.2 问题描述
本文考虑具有$ N $个智能体的线性多智能体系统, 第$ i $个智能体的状态表示为
$$ \begin{split} \dot{x}_i(t) = Ax_i(t)+Bu_i(t)\\ \end{split} $$ (1) 其中, $ x_i(t) \in {\bf{R}}^z $表示智能体$ i $的系统状态, $ u_i(t) \in {\bf{R}}^n $表示智能体$ i $的控制输入, $ A $和$ B $为已知的具有适当维数的矩阵.
当系统发生执行器故障时, 智能体$ i $的控制输入可以改写成如下形式
$$ \begin{equation} u_i ^F(t) = ({\boldsymbol{I}}_{n}-\Lambda_i) u_i(t) \end{equation} $$ (2) 本文考虑到智能体控制输入的每一个执行器均有可能发生故障, 其中$ {\boldsymbol{I}}_{n} $表示$ n $维单位矩阵, $ \Lambda_i = \text{diag}\{\rho_{i1},\;\rho_{i2},\;\cdots ,\;\rho_{in}\} $, $ \rho_{ij} $为智能体$ i $的第$ j $个执行器的故障因子. 定义
$$ \begin{split} &\Lambda = \text{diag}\{ \Lambda_1,\;\Lambda_2,\; \cdots ,\;\Lambda_N\}\\& \hat{\Lambda}_i = \text{diag}\{\hat{\rho}_{i1},\;\hat{\rho}_{i2},\; \cdots ,\;\hat{\rho}_{in}\}\\& \hat{\Lambda} = \text{diag}\{\hat{\Lambda}_1,\;\hat{\Lambda}_2,\; \cdots ,\;\hat{\Lambda}_N\}\\& \bar{\Lambda}_i = \text{diag} \{\bar{\rho}_{i1},\;\bar{\rho}_{i2},\; \cdots ,\;\bar{\rho}_{in}\}\\& \bar{\Lambda} = \text{diag}\{\bar{\Lambda}_1,\;\bar{\Lambda}_2,\; \cdots ,\;\bar{\Lambda}_N\} \end{split} $$ $$ \begin{split} & \underline{\Lambda}_i = \text{diag}\{\underline{\rho}_{i1},\;\underline{\rho}_{i2},\; \cdots ,\;\underline{\rho}_{in}\}\\& \underline{\Lambda} = \text{diag}\{\underline{\Lambda}_1,\;\underline{\Lambda}_2,\; \cdots ,\;\underline{\Lambda}_N\}\\& \bar{\hat{\Lambda}}_i = \text{diag}\{\bar{\hat{\rho}}_{i1},\;\bar{\hat{\rho}}_{i2},\; \cdots ,\;\bar{\hat{\rho}}_{in}\}\\& \bar{\hat{\Lambda}} = \text{diag}\{\bar{\hat{\Lambda}}_1,\;\bar{\hat{\Lambda}}_2,\; \cdots ,\;\bar{\hat{\Lambda}}_N\}\\& \underline{\hat{\Lambda}}_i = \text{diag}\{\underline{\hat{\rho}}_{i1},\;\underline{\hat{\rho}}_{i2},\; \cdots ,\;\underline{\hat{\rho}}_{in}\}\\& \nonumber \underline{\hat{\Lambda}} = \text{diag}\{\underline{\hat{\Lambda}}_1,\;\underline{\hat{\Lambda}}_2,\; \cdots ,\;\underline{\hat{\Lambda}}_N\} \end{split} $$ 其中, $ \hat{\rho}_{ij} $为待设计的自适应参数; $ \bar{\rho}_{ij} $和$ \underline{\rho}_{ij} $分别表示$ \rho_{ij} $的上界和下界, $ \bar{\hat{\rho}}_{ij} $和$ \underline{\hat{\rho}}_{ij} $分别表示$ \hat{\rho}_{ij} $的上界和下界.
此外, 故障存在上界和下界, 即$ 0 < \underline{\hat{\rho}}_{ij} \leq \hat{\rho}_{ij} \leq $ $ \bar{\hat{\rho}}_{ij} <1 $, $ \bar{\rho}_{ij} = \bar{\hat{\rho}}_{ij} $, $ \underline{\rho}_{ij} = \underline{\hat{\rho}}_{ij} $, $ \bar{\Lambda}_i = \bar{\hat{\Lambda}}_i $, $ \underline{\Lambda}_i = \underline{\hat{\Lambda}}_i $.
假设 2. 系统状态矩阵$ A $和$ B $是可控的.
引理 2. 芭芭拉引理[34], 如果$ y(t) $和$ \dot{y}(t) $有界, 并且$ \int_{0}^{\infty}y^{\rm{T}}(t)y(t)\text{d}t<\infty $, 则$ \lim_{t\to\infty} y(t) = 0 $.
引理 3. 杨氏不等式[38], 如果$ x>1 $, $ y>1 $, 并且${1}/{x}+{1}/{y} = 1 $, 则$||a^{\rm{T}}b||\leq{||a^x||}/x+||b^y||/{y}$.
2. 主要结果
2.1 通信协议和容错控制器设计
本节提出一个新的事件触发通信协议和容错控制器$ u_i(t_k^i) $. 智能体$ i $的通信触发时刻为$ t_k^i $, 其中$ k = 0,\;1,\cdots $. 智能体$ i $在$ [0,\;t_{k+1}^i] $期间接收到智能体$ m $状态信息的时刻为$ {t_{k^m}^i} $. 多智能体系统的一致性误差为
$$ \begin{equation} \xi_i(t) = \sum\limits_{m = 1}^{N} a_{im}(x_i(t)-x_m(t)) \nonumber \end{equation} $$ $$ \begin{equation} \xi_i(t_k^i) = \sum\limits_{m = 1}^{N} a_{im}(x_i(t_k^i)-x_m(t_{k^m}^i)) \nonumber \end{equation} $$ 事件触发状态误差$ e_{i}(t) $为$ e_{i}(t) = x_i(t_k^i) - x_i(t) $, 事件触发一致性误差$ e_{\xi_i}(t) $定义为
$$ \begin{split} e_{\xi_i}(t)& = \xi_i(t_k^i)-\xi_i(t)\\ \end{split} $$ (3) 定义
$$ \begin{split} &\xi^{\text{T}}(t) = (\xi_1^{\text{T}}(t),\;\xi_2^{\text{T}}(t),\;\cdots ,\;\xi_N^{\text{T}}(t))\\&e^{\text{T}}_{\xi}(t) = (e^{\text{T}}_{\xi_1}(t),\;e^{\text{T}}_{\xi_2}(t),\;\cdots ,\; e^{\text{T}}_{\xi_N}(t))\\&\xi^{\text{T}}(t_k) = (\xi^{\text{T}}_1(t_k^1),\;\xi^{\text{T}}_2(t_k^2),\;\cdots ,\;\xi^{\text{T}}_N(t_k^N)) \nonumber \end{split} $$ 自适应参数$ \hat{\rho}_{ij} $的引入可以提升控制器的容错能力, 它被设计为如下形式
$$ \dot{\hat{\rho}}_{ij} = \text{Proj}_{[\bar{\rho}_{ij},\;\underline{\rho}_{ij}]}(\hat{\rho}_{ij}) = \left\{ \begin{aligned} & 0 , \left\{ \begin{aligned} &\hat{\rho}_{ij}(t) = \bar{\rho}_{ij}\; \text{且}\; \theta_{ij}> 0 \\& \text{或者} \\ &\hat{\rho}_{ij}(t) = \underline{\rho}_{ij} \;\text{且}\; \theta_{ij} < 0 \end{aligned} \right.\\ & \tau_{ij} \theta_{ij} ,\quad \quad \quad \quad \quad \quad \;\text{否则} \end{aligned} \right. $$ (4) 其中, $\text{Proj}_{[\bar{\rho}_{ij},\;\underline{\rho}_{ij}]}(\hat{\rho}_{ij})$表示投影运算, $ \hat{\rho}_{ij} $被投影到区间$[\bar{\rho}_{ij},\;\underline{\rho}_{ij}]$内, $ \tau_{ij} > 0 $, $ \theta_{ij} $被定义为
$$ \theta_{ij} = (1+\hat{\rho}_{ij}(t))\Vert \xi_i^{\rm{T}}(t_k^i)[PB]_{j}\Vert^2-\sigma_{ij}(t)\hat{\rho}_{ij}(t) $$ (5) 其中, $ \sigma_{ij}(t) $满足$\int_{0}^{ +\infty} \sigma_{ij}(t) {\rm{d}}t = 0$, $ P $是一个待求解的正实对称矩阵, $ [PB]_{j} $为矩阵$ PB $的第$ j $列.
控制器$ u_i(t_k^i) $被设计为
$$ \begin{split} u_i(t_k^i) = &- ({\boldsymbol{I}}_{n}+\hat {\Lambda}_i(t) )B^{\rm{T}}P\xi_i(t_k^i) \end{split} $$ (6) 注 1. 不同于文献[34]中所设计的容错控制器, 本文控制器$ u_i(t_k^i) $中的$ \xi_i(t_k^i) $不仅在自身达到触发条件时更新, 其邻居达到触发条件时也会更新.
注 2. 与文献[35-36, 40]的工作不同, 控制器$ u_i(t_k^i) $在触发间隔内是时变的. 事件触发通信协议仅决定智能体间是否通信, 通信协议的设计初衷是避免相邻智能体的实时或者频繁通信, 基于通信结果完成控制器更新. 因此, 控制器$ u_i(t_k^i) $在触发间隔内时变也是合理的. 基于自适应技术设计的自适应参数$ \hat{\rho}_{ij} $可以补偿执行器失效带来的负面影响. 此外, 得益于触发间隔内时变的控制器$ u_i(t_k^i) $, 系统具有更好的动态响应.
智能体$ i $的系统动态为
$$ \dot{x}_i(t) = Ax_i(t)+Bu_i(t_k^i)$$ (7) 对于多智能体系统而言, 需要针对每一个智能体给出其事件触发通信条件. 动态事件触发协议为
$$ t_{k+1}^i= \text{inf} \{t> t_{k}^i :f_i(e_i(t),\;\xi_i(t_k^i),\;\eta_i(t))\leq 0 \} $$ (8) 其中, $ \text{inf} $表示集合的下确界, 动态事件触发函数$ f_i(e_i(t), \;\xi_i(t_k^i),\;\eta_i(t)) $被设计为如下形式
$$ \begin{split} &f_i(e_i(t),\xi_i(t_k^i),\eta_i(t))=\\& \qquad(\omega_i\xi^{\rm{T}}_i(t_k^i) PB({\boldsymbol{I}}_n+\hat{\Lambda}_i(t)) B^{\rm{T}}P\xi_i(t_k^i)\;-\\&\qquad\kappa \gamma e^{\rm{T}}_i(t)PBB^{\rm{T}}Pe_i(t))\varphi_i+\eta_i(t) \end{split} $$ (9) 其中, $ \varphi_i $为正实数, 动态变量$ \eta_i(t) $满足如下条件
$$ \dot{\eta}_i(t) = \left\{ \begin{aligned} & -\lambda \eta_i(t)+\vartheta_i \text{e}^{-\digamma_it}, & & \Gamma > 0 \\ & -\lambda \eta_i(t) - \Gamma +\vartheta_i \text{e}^{-\digamma_it}, & & \text{否则} \end{aligned} \right. $$ (10) 其中, $ \Gamma $为
$$ \begin{split} \Gamma=\;& \omega_i\xi^{\rm{T}}_i(t_k^i) PB({\boldsymbol{I}}_{n}+\hat{\Lambda}_i(t)) B^{\rm{T}}P\xi_i(t_k^i)\;-\\ &\kappa \gamma e^{\rm{T}}_i(t)PBB^{\rm{T}}Pe_i(t) \end{split} $$ (11) 其中, $ \omega_i $, $ \kappa $, $ \gamma $, $ \lambda $, $ \vartheta_i $和$ \digamma_i $均为已知正实数.
采用事件触发通信协议是为了减少智能体间不必要的通信, 智能体$ i $在未通信时无法获取邻居的实时信息, 事件触发函数$ f_i(e_i(t), \;\xi_i(t_k^i),\; \eta_i(t)) $仅需要智能体自身的状态信息及其邻居在上一触发时刻的状态信息.
在触发函数中引入动态变量$ \eta_i(t) $可以提高触发阈值, 延长触发间隔$ \{t_{k+1}^i-t_{k}^i \} $. 文献[33]通过在事件触发函数中引入正常量来进一步延长触发间隔, 但会对多智能体系统的一致性产生负面影响. 本文设计的触发机制(9)不仅可以延长触发间隔, 还可以满足多智能体系统的最终一致性要求.
多智能体系统可以表示为
$$ \begin{split} \dot{x}(t) = \;&({\boldsymbol{I}}_N\otimes A)x(t)-({\boldsymbol{I}}_N\otimes B)({\boldsymbol{I}}_{Nn}-\Lambda\;+\\ &({\boldsymbol{I}}_{Nn}-\Lambda) \hat {\Lambda} ) ({\boldsymbol{I}}_N\otimes B^{\rm{T}}P)(e_{\xi}(t)+\xi(t)) \end{split} $$ 状态一致性误差$ \varepsilon $表示为
$$ \begin{split} \varepsilon(t) = (M\otimes {\boldsymbol{I}}_z)x(t) \end{split} $$ (12) 其中, $ M = {\boldsymbol{I}}_N-(1/N) {\bf{1}}_N\otimes{\bf{1}}_N^{\text{T}} $, $ {\bf{1}}_N $表示元素值全为$ 1 $的列向量. 容错控制器(6)的设计目的是使多智能体系统状态一致性误差$ \varepsilon $收敛至0
$$ \begin{split} \dot{\varepsilon}(t)=\;& (M\otimes A)x(t)-(M\otimes {\boldsymbol{I}}_z)({\boldsymbol{I}}_N\otimes B)({\boldsymbol{I}}_{Nn}-\Lambda)\;\times\\&({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t)) ({\boldsymbol{I}}_N\otimes B^{\rm{T}}P)\xi(t_k) \end{split} $$ 传统事件触发控制的设计目的是减少控制器$ u_i(t_k^i) $的更新次数, 但在多智能体系统中由于通信带宽的限制, 智能体间不能频繁通信. 文中采用事件触发机制(8)来决定智能体是否与邻居通信, 在触发通信间隔内控制器$ u_i(t_k^i) $只使用自身的状态信息和邻居在触发时刻的已知信息, 减少了不必要的的通信.
注 3. 自适应参数$ \hat{\Lambda}_i(t) $不仅可以补偿执行器故障, 还可以延长触发间隔$\{ t_{k+1}^i-t_k^i \}$.
2.2 一致性证明
针对多智能体系统事件触发可靠一致性控制问题, 本文给出如下定理.
定理 1. 如果假设1和假设2成立, 控制器(6)和事件触发机制(8)可以保证多智能体系统在发生执行器故障时的一致性.
证明. 构造李雅普诺夫函数$ V(t) = V_1(t)+ V_2(t)\;+ V_3(t) $, 其中
$$ \begin{split} V_1(t)& = \varepsilon^{\text{T}}(t)(L\otimes P)\varepsilon(t)\\ V_2(t)& = \sum\limits_{i = 1}^{N}\sum\limits_{j = 1}^{n}\frac{(\hat{\rho}_{ij}-\rho_{ij})^2}{4\tau_{ij}}\\ V_3(t)& = \sum\limits_{i = 1}^{N} \eta_i(t) \end{split} $$ 已知$ LM = L $, 等式可以改写为
$$ \begin{split} \dot{ V}_1(t) =\;& \varepsilon^{\text{T}}(t)(L\otimes (A^{\text{T}}P+PA))\varepsilon(t)\;-\\ &2\varepsilon^{\text{T}}(t)(L\otimes PB)({\boldsymbol{I}}_{Nn}-\Lambda\;+\\ &({\boldsymbol{I}}_{Nn}-\Lambda)\hat{\Lambda}(t) )({\boldsymbol{I}}_N\otimes B^{\text{T}}P)\xi(t_k) \end{split} $$ (13) 易得$ \xi(t) = (L\otimes {\boldsymbol{I}}_z)\varepsilon(t) $, 则下式成立
$$ \begin{split} \dot{ V}_1(t) =\;& \varepsilon^{\text{T}}(t)(L\otimes (A^{\text{T}}P+PA))\varepsilon(t)\;-\\ &2\xi^{\text{T}}(t)({\boldsymbol{I}}_N\otimes PB)({\boldsymbol{I}}_{Nn}-\Lambda\;+\\ &({\boldsymbol{I}}_{Nn}-\Lambda)\hat{\Lambda}(t) )({\boldsymbol{I}}_N\otimes B^{\text{T}}P)\xi(t_k) \end{split} $$ (14) 已知$ e_{\xi_i}(t) = \xi_i(t_k^i)-\xi_i(t) $, 可得
$$ \begin{split} \dot{ V}_1(t) =\;& \varepsilon^{\text{T}}(t)(L\otimes (A^{\text{T}}P+PA))\varepsilon(t)\;-\\&2\xi^{\text{T}}(t)({\boldsymbol{I}}_N\otimes PB)\Omega ({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t)) ({\boldsymbol{I}}_N\;\otimes\\ &B^{\text{T}}P)e_{\xi}(t)-2\xi^{\text{T}}(t)({\boldsymbol{I}}_N\otimes PB)\Omega({\boldsymbol{I}}_{Nn}\;+\\ &\hat{\Lambda}(t))({\boldsymbol{I}}_N\otimes B^{\text{T}}P)\xi(t) \\[-1pt]\end{split} $$ (15) 其中, $ {\boldsymbol{I}}_{Nn}-\Lambda = \Omega $.
根据$ e_{\xi_i}(t) = \xi_i(t_k^i)-\xi_i(t) $, 式(15)可写为如下形式
$$ \begin{split} \dot{ V}_1(t)=\;& \varepsilon^{\text{T}}(t)(L\otimes (A^{\text{T}}P+PA))\varepsilon(t)\;-\\&2\xi^{\text{T}}(t_k)({\boldsymbol{I}}_N\otimes PB)\Omega({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t))({\boldsymbol{I}}_N\;\otimes\\& B^{\text{T}}P)e_{\xi}(t)+2e_{\xi}^{\text{T}}(t)({\boldsymbol{I}}_N\otimes PB)\Omega({\boldsymbol{I}}_{Nn}\;+\\&\hat{\Lambda}(t))({\boldsymbol{I}}_N\otimes B^{\text{T}}P)e_{\xi}(t)-2\xi^{\text{T}}(t)({\boldsymbol{I}}_N\;\otimes\\& PB)\Omega({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t))({\boldsymbol{I}}_N\otimes B^{\text{T}}P)\xi(t)\\[-1pt] \end{split} $$ (16) 已知$ e_{\xi_i}(t) = \xi_i(t_k^i)-\xi_i(t) $, 式(16)可改写为
$$ \begin{split} \dot{ V}_1(t)=\;& \varepsilon^{\text{T}}(t)(L\otimes (A^{\text{T}}P+PA))\varepsilon(t)\;-\\&2\xi^{\text{T}}(t_k)({\boldsymbol{I}}_N\otimes PB)\Omega({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t)) \nonumber ({\boldsymbol{I}}_N\;\otimes\\& B^{\text{T}}P)\xi(t_k)+2\xi^{\text{T}}(t_k)({\boldsymbol{I}}_N\otimes PB)\Omega({\boldsymbol{I}}_{Nn}\;+\\&\hat{\Lambda}(t))({\boldsymbol{I}}_N\otimes B^{\text{T}}P)\xi(t)+2e_{\xi}^{\text{T}}(t)({\boldsymbol{I}}_N\;\otimes\\& PB)\Omega({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t))({\boldsymbol{I}}_N\otimes B^{\text{T}}P)e_{\xi}(t)\;-\\&2\xi^{\text{T}}(t)({\boldsymbol{I}}_N\otimes PB)\Omega({\boldsymbol{I}}_{Nn}\;+\\&\hat{\Lambda}(t))({\boldsymbol{I}}_N\otimes B^{\text{T}}P)\xi(t) \end{split} $$ 进一步, 根据引理3可得
$$ \begin{split} & 2\xi^{\text{T}}(t_k)({\boldsymbol{I}}_N\otimes PB)\Omega({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t))({\boldsymbol{I}}_N\otimes B^{\text{T}}P)\xi(t) \leq \\&\qquad \xi^{\text{T}}(t_k)({\boldsymbol{I}}_N\otimes PB)\Omega({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t))({\boldsymbol{I}}_N\;\otimes \\&\qquad B^{\text{T}}P)\xi(t _k)+\xi^{\text{T}}(t)({\boldsymbol{I}}_N\otimes PB)\Omega({\boldsymbol{I}}_{Nn}\;+\\&\qquad\hat{\Lambda}(t))({\boldsymbol{I}}_N\otimes B^{\text{T}}P)\xi(t)\\[-1pt] \end{split} $$ (17) 根据不等式(17), 可得
$$ \begin{split} \dot{ V}_1(t) \leq\;& \varepsilon^{\text{T}}(t)(L\otimes (A^{\text{T}}P+PA))\varepsilon(t)\;-\\&\xi^{\rm{T}}(t_k)({\boldsymbol{I}}_N\otimes PB)\Omega({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t))({\boldsymbol{I}}_N\;\otimes\\& B^{\text{T}}P)\xi(t_k)+2e_{\xi}^{\text{T}}(t)({\boldsymbol{I}}_N\otimes PB)\Omega({\boldsymbol{I}}_{Nn}\;+\\&\hat{\Lambda}(t))({\boldsymbol{I}}_N\otimes B^{\text{T}}P)e_{\xi}(t)-\xi^{\text{T}}(t)({\boldsymbol{I}}_N\;\otimes\\& PB)\Omega({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t))({\boldsymbol{I}}_N\otimes B^{\text{T}}P)\xi(t)\\[-1pt] \end{split} $$ (18) 如果$ \hat{\rho}_{ij}(t) = \bar{\rho}_{ij} $且$ \theta_{id} > 0 $, 或者$ \hat{\rho}_{ij}(t) = $$\underline{\rho}_{ij} $且$ \theta_{id} < 0 $, 可以得到$ 2(\hat{\rho}_{ij}-\rho_{ij}) \theta_{id} \geq 0 $.
$ \dot{V}_2(t) $满足如下不等式
$$ \begin{split} \dot{V}_2(t)\leq\;& \sum\limits_{i = 1}^{N}\sum\limits_{j = 1}^{n}\frac{1}{2}(\hat{\rho}_{ij}-\rho_{ij})((1\;+\\ &\hat{\rho}_{ij}(t))\Vert \xi_i^{\text{T}}(t_k^i)[PB]_{j}\Vert^2-\sigma_{ij}(t)\hat{\rho}_{ij}(t))\end{split} $$ (19) 已知$ \hat{\rho}_{ij}(t)\leq \bar{\hat{\rho}}_{ij} $, 则下式成立
$$ \begin{split} \dot{V}_2(t)\leq\;& \sum\limits_{i = 1}^{N}\sum\limits_{j = 1}^{n}\frac{1}{2}( \bar{\hat{\rho}}_{ij}-\rho_{ij})(1\;+\\ &\hat{\rho}_{ij}(t))\Vert \xi_i^{\text{T}}(t_k^i)[PB]_{j}\Vert^2\;-\\&\sum\limits_{i = 1}^{N}\sum\limits_{j = 1}^{n}\frac{1}{2}(\hat{\rho}_{ij}-\rho_{ij})\sigma_{ij}(t)\hat{\rho}_{ij}(t) \end{split} $$ (20) 基于等价变换, 可得
$$ \begin{split} \dot{V}_2(t)\leq\;&\frac{1}{2} \xi^{\text{T}}(t_k)({\boldsymbol{I}}_N\otimes PB)(\bar{\Lambda}-\Lambda)({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t))\;\times\\&({\boldsymbol{I}}_N\otimes B^{\text{T}}P)\xi(t_k)\;-\\ &\frac{1}{2}\sum\limits_{i = 1}^{N}\sum\limits_{j = 1}^{n}(\hat{\rho}_{ij}-\rho_{ij})\sigma_{ij}(t)\hat{\rho}_{ij}(t) \\[-1pt]\end{split} $$ (21) 已知$ \xi(t_k) = e_{\xi}(t)+\xi(t) $, 可得
$$ \begin{split} \dot{V}_2(t)\leq\;& \xi^{\text{T}}(t)({\boldsymbol{I}}_N\otimes PB)(\bar{\Lambda}-\Lambda)({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t))({\boldsymbol{I}}_N\;\otimes\\& B^{\text{T}}P)\xi(t)+e^{\text{T}}_{\xi}(t)({\boldsymbol{I}}_N\otimes PB)(\bar{\Lambda}-\Lambda)({\boldsymbol{I}}_{Nn}\;+\\ &\hat{\Lambda}(t))({\boldsymbol{I}}_N\otimes B^{\text{T}}P)e_{\xi}(t)\;-\\ &\sum\limits_{i = 1}^{N}\sum\limits_{j = 1}^{n}\frac{1}{2}(\hat{\rho}_{ij}-\rho_{ij})\sigma_{ij}(t)\hat{\rho}_{ij}(t) \\[-1pt]\end{split} $$ (22) 根据触发函数(8), 可得
$$ \begin{equation} \Gamma>-\frac{1}{\varphi_i}\eta_i(t) \end{equation} $$ (23) 关于动态变量$ \eta_i(t) $和$ \dot{V}_3 (t) $, 如下不等式成立
$$ \begin{split}& -\lambda \eta_i(t)+\vartheta_i \text{e}^{-\digamma_it}\leq \dot{\eta}_i(t)<\\ &\qquad-\left(\lambda- \frac{1} {\varphi_i}\right) \eta_i(t)+\vartheta_i \text{e}^{-\digamma_it} \end{split} $$ (24) $$ \begin{split} &\sum\limits_{i = 1}^{N}-\lambda \eta_i(t)+\vartheta_i \text{e}^{-\digamma_it}\leq\dot{ V}_3(t)<\\ &\qquad\sum\limits_{i = 1}^{N} -\left(\lambda- \frac{1} {\varphi_i}\right) \eta_i(t)+\vartheta_i \text{e}^{-\digamma_it} \end{split} $$ (25) 综合以上结果, $ \dot V(t) $满足如下不等式
$$ \begin{split} \dot V(t)=\;& \dot V_1(t)+\dot V_2(t)+\dot V_3(t)<\\&\varepsilon^{\text{T}}(t)(L\otimes (A^{\rm{T}}P+PA))\varepsilon(t)\;-\\&\xi^{\text{T}}(t_k)({\boldsymbol{I}}_N\otimes PB)\Omega({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t)) ({\boldsymbol{I}}_N\;\otimes\\& B^{\text{T}}P)\xi(t_k)+2e_{\xi}^{\text{T}}(t)({\boldsymbol{I}}_N\otimes PB)\Omega({\boldsymbol{I}}_{Nn}\;+\\&\hat{\Lambda}(t))({\boldsymbol{I}}_N\otimes B^{\text{T}}P)e_{\xi}(t)-\xi^{\text{T}}(t)({\boldsymbol{I}}_N\;\otimes\\& PB)\Omega({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t))({\boldsymbol{I}}_N\otimes B^{\text{T}}P)\xi(t)\;+\\ & \xi^{\text{T}}(t)({\boldsymbol{I}}_N\otimes PB)(\bar{\Lambda}-\Lambda)({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t))({\boldsymbol{I}}_N\;\otimes\\& B^{\text{T}}P)\xi(t)+e^{\text{T}}_{\xi}(t)({\boldsymbol{I}}_N\otimes PB)(\bar{\Lambda}-\Lambda)({\boldsymbol{I}}_{Nn}\;+\\&\hat{\Lambda}(t))({\boldsymbol{I}}_N\otimes B^{\text{T}}P)e_{\xi}(t)\;-\\&\sum\limits_{i = 1}^{N}\sum\limits_{j = 1}^{n}\frac{1}{2}(\hat{\rho}_{ij}-\rho_{ij})\sigma_{ij}(t)\hat{\rho}_{ij}(t)\;- \\&\sum\limits_{i = 1}^{N} \left(\left(\lambda- \frac{1} {\varphi_i}\right)\eta_i(t)-\vartheta_i \text{e}^{-\digamma_it}\right) \end{split} $$ 合并同类项, 可得
$$ \begin{split} \dot V(t) <&\;\varepsilon^{\text{T}}(t)(L\otimes (A^{\text{T}}P+PA))\varepsilon(t)\;-\\&\xi^{\text{T}}(t_k)({\boldsymbol{I}}_N\otimes PB)\Omega({\boldsymbol{I}}_{Nn}+\hat{\Lambda})({\boldsymbol{I}}_N\;\otimes\\& B^{\text{T}}P)\xi(t_k)-\xi^{\text{T}}(t)({\boldsymbol{I}}_N\otimes PB) \Pi ({\boldsymbol{I}}_N\;\otimes\\& B^{\text{T}}P)\xi(t)+e^{\text{T}}_{\xi}(t)({\boldsymbol{I}}_N\otimes PB)\Psi ({\boldsymbol{I}}_{Nn}\;+\\&\hat{\Lambda}(t))({\boldsymbol{I}}_N\otimes B^{\text{T}}P)e_{\xi}(t)\;-\\&\sum\limits_{i = 1}^{N}\sum\limits_{j = 1}^{n}\frac{1}{2}(\hat{\rho}_{ij}-\rho_{ij})\sigma_{ij}(t)\hat{\rho}_{ij}(t)\;-\\&\sum\limits_{i = 1}^{N} \left(\left(\lambda- \frac{1} {\varphi_i}\right)\eta_i(t)-\vartheta_i \text{e}^{-\digamma_it}\right) \end{split} $$ (26) 其中, $ \Pi = {\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t)-\bar{\Lambda}-\bar{\Lambda}\hat{\Lambda}(t) $, $ \Psi = 2{\boldsymbol{I}}_{Nn}+ \bar{\Lambda}\;-$ $3\Lambda $.
$$ \begin{split} &\sum\limits_{i = 1}^{N}e^{\text{T}}_{\xi _i}(t) PBB^{\text{T}}P e_{\xi _i}(t) = e^{\text{T}}(t)(L^2\otimes PBB^{\text{T}}P)e(t) \leq\\&\qquad \gamma \sum\limits_{i = 1}^{N}e^{\text{T}}_i(t)PBB^{\text{T}}Pe_i(t) \\[-1pt]\end{split} $$ (27) 其中, $ \gamma \;=\; \lambda_{\max}(L^2\otimes {\boldsymbol{I}}_n) \;=\; \lambda_{\max}(L^2)\lambda_{\max}( {\boldsymbol{I}}_n)\; =$ $ \lambda_{\max}(L^2) $, 定义$ \kappa = \max \nolimits_{i,\;j} ((\bar{\hat{\rho}}_{ij}+1)(2+\bar{\hat{\rho}}_{ij}-3\underline{\rho}_{ij})) $, $ \omega_i = \min \nolimits_{i} (1-\bar \rho _{ij}) $.
对矩阵进行变换, 可得
$$ \begin{split} \dot V(t) <\;&\varepsilon^{\text{T}}(t)(L\otimes (A^{\text{T}}P+PA))\varepsilon(t)\;-\\&\xi^{\text{T}}(t)({\boldsymbol{I}}_N\otimes PB) \Pi ({\boldsymbol{I}}_N\otimes B^{\text{T}}P)\xi(t)\;-\\&\sum\limits_{i = 1}^{N} \omega_i\xi^{\text{T}}_i(t_k^i) PB({\boldsymbol{I}}_n+\hat{\Lambda}_i(t)) B^{\text{T}}P\xi_i(t_k^i)\; +\\ &\sum\limits_{i = 1}^{N}\kappa e^{\text{T}}_{\xi_i}(t) PBB^{\text{T}}Pe_{\xi_ i}(t)\;-\\&\sum\limits_{i = 1}^{N}\sum\limits_{j = 1}^{n}\frac{1}{2}(\hat{\rho}_{ij}-\rho_{ij})\sigma_{ij}(t)\hat{\rho}_{ij}(t)\;-\\& \sum\limits_{i = 1}^{N} \left(\left(\lambda- \frac{1} {\varphi_i}\right)\eta_i(t)-\vartheta_i \text{e}^{-\digamma_it}\right)\\[-1pt] \end{split} $$ (28) 根据不等式(27), 可得
$$ \begin{split} \dot V(t) <\;& \varepsilon^{\text{T}}(t)(L\otimes (A^{\text{T}}P+PA))\varepsilon(t)\;-\\&\xi^{\text{T}}(t)({\boldsymbol{I}}_N\otimes PB) \Pi ({\boldsymbol{I}}_N\otimes B^{\text{T}}P)\xi(t)\;-\\&\sum\limits_{i = 1}^{N}\omega_i\xi^{\text{T}}_i(t_k^i) PB({\boldsymbol{I}}_n+\hat{\Lambda}_i(t)) B^{\text{T}}P\xi_i(t_k^i)\;+\\&\kappa \gamma \sum\limits_{i = 1}^{N}e^{\text{T}}_i(t)PBB^{\text{T}}Pe_i(t)\;-\\& \sum\limits_{i = 1}^{N}\sum\limits_{j = 1}^{n}\frac{1}{2}(\hat{\rho}_{ij}-\rho_{ij})\sigma_{ij}(t)\hat{\rho}_{ij}(t)\;- \\& \sum\limits_{i = 1}^{N} \left(\left(\lambda- \frac{1} {\varphi_i}\right) \eta_i(t)-\vartheta_i \text{e}^{-\digamma_it}\right) \end{split} $$ 基于所设计的触发函数(8), 上述不等式可以改写为
$$ \begin{split} \dot V(t) <\;&\varepsilon^{\text{T}}(t)(L\otimes (A^{\text{T}}P+PA))\varepsilon(t)\;-\\&\sum\limits_{i = 1}^{N}\sum\limits_{j = 1}^{n}\frac{1}{2}(\hat{\rho}_{ij}-\rho_{ij})\sigma_{ij}(t)\hat{\rho}_{ij}(t)\;-\\&\xi^{\text{T}}(t)({\boldsymbol{I}}_N\otimes PB) \Pi({\boldsymbol{I}}_N\otimes B^{\text{T}}P)\xi(t)\;+\\& \sum\limits_{i = 1}^{N} \frac{1} {\varphi_i}\eta_i(t) - \sum\limits_{i = 1}^{N} \left( \left(\lambda- \frac{1} {\varphi_i}\right) \eta_i(t)-\vartheta_i \text{e}^{-\digamma_it}\right) \end{split} $$ (29) 对不等式(29)合并同类项, 可得
$$ \begin{split} \dot V(t) <\;& \varepsilon^{\text{T}}(t)(L\otimes (A^{\text{T}}P+PA))\varepsilon(t)+\sum\limits_{i = 1}^{N} \vartheta_i \text{e}^{-\digamma_it}\;-\\&\varepsilon^{\text{T}}(t)(L\otimes PB) \Pi (L\otimes B^{\text{T}}P)\varepsilon(t) \;-\\&\sum\limits_{i = 1}^{N}\sum\limits_{j = 1}^{n}\frac{1}{2}(\hat{\rho}_{ij}-\rho_{ij})\sigma_{ij}(t)\hat{\rho}_{ij}(t)\;-\\ &\sum\limits_{i = 1}^{N} \left(\lambda- \frac{2} {\varphi_i}\right) \eta_i(t) \\[-1pt]\end{split} $$ (30) 其中, 需要确保以下不等式成立
$$ \begin{split} L\otimes (A^{\text{T}}P+PA)-(L\otimes PB) \Pi (L\otimes B^{\text{T}}P) < 0 \end{split} $$ (31) 在不等式(31)左右两边同时乘以$ {\boldsymbol{I}}_N\otimes P^{-1} $, 易得
$$ \begin{equation*} \begin{split} L\otimes (P^{-1}A^{\text{T}}+AP^{-1}) -(L\otimes B) \Pi (L\otimes B^{\text{T}}) <0 \end{split} \end{equation*} $$ 已知$ \Pi \geq ({\boldsymbol{I}}_{Nn}+\underline{\hat{\Lambda}})({\boldsymbol{I}}_{Nn}-\bar{\Lambda}) $, 可得
$$ \begin{split} &\varepsilon^{\text{T}}(t)(L\otimes (P^{-1}A^{\text{T}}+AP^{-1})-(L\otimes B) \Pi (L\;\otimes \\ &\qquad B^{\text{T}})) \varepsilon(t) \leq \varepsilon^{\text{T}}(t)(L\otimes (P^{-1}A^{\text{T}}+AP^{-1})\;-\\&\qquad(L\otimes B) ({\boldsymbol{I}}_{Nn}+\underline{\Lambda})({\boldsymbol{I}}_{Nn}-\bar{\Lambda}) (L\otimes B^{\text{T}})) \varepsilon(t) \leq \\&\qquad\varepsilon^{\text{T}}(t)(L\otimes (P^{-1}A^{\text{T}}+AP^{-1})-\delta(L^2\otimes B B^{\text{T}})) \varepsilon(t) \end{split} $$ (32) 其中, $ \delta $是正实常量, $ \delta = \min \nolimits_{i,\;j} ((\underline{\rho}_{ij}+1)(1- \bar{\hat{\rho}}_{ij})) $.
基于引理1, 必将存在一个酉矩阵$ U = [({\bf{1}}/ \sqrt{N}) Y] $, 使得$ U^{\text{T}}LU $ = $ \text{diag}\{0,\lambda_2,\cdots ,\lambda_N\} $. 令$\bar{ \varepsilon} = (U^{\text{T}}\;\otimes {\boldsymbol{I}}_z )\varepsilon$, $ \lambda_2 $表示拉普拉斯矩阵$ L $的最小特征根.
基于假设1, 不等式(32)可作如下变换
$$ \begin{split} & \varepsilon^{\text{T}}(t)(L\otimes (P^{-1}A^{\text{T}}+AP^{-1})-\delta(L^2\otimes B B^{\text{T}})) \varepsilon(t) \leq \\& \sum\limits_{i = 2}^{N}\lambda_i( \bar{ \varepsilon}_i^{\text{T}}(t) (P^{-1}A^{\text{T}}+AP^{-1}-\lambda_2\delta B B^{\text{T}}) \bar{ \varepsilon}_i(t) )<0 \end{split} $$ (33) 引理 4. 基于文献[41]中的黎卡提不等式$ PA\;+ A^{\text{T}}P-\varpi PBB^{\text{T}}P+Q<0 $, 在系统状态矩阵$ (A,\;B) $可控的情况下, $ P $可通过如下线性矩阵不等式求解. 其中$ P $和$ Q $均为正定实对称矩阵, $ \varpi $为正实数
$$ \left[ {\begin{array}{cc} P^{-1}A^{\text{T}}+ AP^{-1}-\varpi BB^{\text{T}} & P^{-1} \\ P^{-1} & -Q^{-1} \end{array} } \right] <0$$ (34) 注 4. 根据引理4, 不等式(33)一定是可解的, 并且文中所提不等式(31)具有较小的保守性. 当拉普拉斯矩阵$ L $已知时, 可以通过求解不等式(31)获得矩阵$ P $, 进而完成控制器设计. 若仅已知拉普拉斯矩阵$ L $的最小和最大特征值, 可以通过求解不等式(34) 来完成控制策略设计.
对不等式(30)两边从$ 0 $至$ t $开始积分
$$ \begin{split} {V}(t)\;- &\;\int_{0}^{t}\varepsilon^{\text{T}}(t)(L\otimes (A^{\text{T}}P+PA)-(L\;\otimes\\& PB) ({\boldsymbol{I}}_{Nn}+\hat{\Lambda}(t))({\boldsymbol{I}}_{Nn}-\bar{\Lambda}) (L\;\otimes\\& B^{\text{T}}P)) \varepsilon(t) {\text{d}}t < {V}(0)\;-\\&\sum\limits_{i = 1}^{N}\int_{0}^{t} \left(\lambda-\frac{2}{\varphi_i}\right) \eta_i(t) {\text{d}}t+\sum\limits_{i = 1}^{N}\int_{0}^{t} \vartheta_i \text{e}^{-\digamma_it} {\text{d}}t\;-\\&\sum\limits_{i = 1}^{N}\sum\limits_{j = 1}^{n}\int_{0}^{t}\frac{1}{2}(\hat{\rho}_{ij}-\rho_{ij})\sigma_{ij}(t)\hat{\rho}_{ij}(t) {\text{d}}t\\[-1pt] \end{split} $$ (35) 根据芭芭拉引理[34] (引理2)可以得到
$$ \begin{split} \lim_{t\to\infty} \Vert \varepsilon(t) \Vert = 0 \end{split} $$ (36) □ 以上完成了多智能体系统一致性的证明, 在第 2.3节中将证明所设计的事件触发通信机制(8)可以排除芝诺现象.
注 5. 根据不等式(35), 可以看出放松触发条件会对多智能体系统的一致性产生负面影响. 本文设计的事件触发函数$ f_i(e_i(t),\;\xi_i(t_k^i),\;\eta_i(t)) $可以延长触发间隔且多智能体系统最终可以实现一致性.
注 6. 本文设计事件触发机制的初衷是智能体分布式按需通信, 也就是智能体不需要实时获得邻居的信息. 文中所设计的触发通信机制(8)和控制器(6)仅需要获取自身状态信息$ x_i(t) $和触发时刻的邻居信息$ x_m(t_{k^m}^i) $.
2.3 芝诺现象的排除
芝诺现象的出现将会对多智能体系统的稳定性带来影响, 当事件触发通信机制在很短的时间内频繁触发, 智能体间频繁的通信势必会对多智能体系统通信带宽和计算资源带来挑战. 本节将证明所设计的事件触发通信协议(8)不存在芝诺现象.
定理 2. 如果假设1和假设2成立, 事件触发通信机制(8)可以避免芝诺现象的发生.
证明. 在触发间隔$ [t_k^i,\;t_{k+1}^i] $内, $ K e_i(t) $的右上迪尼导数$ \text{D}^+ \Vert K e_i(t)\Vert $可以表示为
$$ \begin{split} \text{D}^+ \Vert K e_i(t) \Vert& \leq \Vert K \dot e_i(t) \Vert = \Vert K \dot x_i(t) \Vert\leq \\& \Vert KA x_i(t)+KB u_i ^F(t) \Vert\leq\\ & \Vert A \Vert \Vert K e_i(t) \Vert+\Vert K A x_i(t_k^i)\Vert\; +\\&\Vert KB u_i ^F(t)\Vert \end{split} $$ (37) 其中, $ K = B^{\text{T}}P $, 令$ \Upsilon = \Vert K Ax_i(t_k^i)\Vert +\Vert KBu_i ^F(t)\Vert $
$$ \begin{split} \Vert K e_i(t) \Vert \leq \frac{\Upsilon}{\Vert A \Vert}(\text{e}^{\Vert A \Vert(t-t_k^i)}-1) \end{split} $$ (38) 根据事件触发条件(8), 可以得到
$$ \begin{split} &\omega_i\xi^{\text{T}}_i(t_k^i) PB\Xi B^{\text{T}}P\xi_i(t_k^i)\;-\\&\qquad\kappa \gamma e^{\text{T}}_i(t)PBB^{\text{T}}Pe_i(t)+\frac{1}{\varphi_i}\eta_i(t)> 0 \end{split} $$ (39) 其中, $ \Xi = {\boldsymbol{I}}_n+\hat{\Lambda}_i(t) $.
当事件触发条件$ f_i(e_i(t),\xi_i(t_k^i),\eta_i(t))\leq0 $满足时, 以下条件成立
$$ \begin{split} &\omega_i\xi^{\text{T}}_i(t_k^i) PB\Xi B^{\text{T}}P\xi_i(t_k^i)\;-\\&\qquad\kappa \gamma e^{\text{T}}_i(t)PBB^{\text{T}}Pe_i(t)+\frac{1}{\varphi_i}\eta_i(t) \leq 0 \end{split} $$ (40) $$ \begin{split} &\omega_i\xi^{\text{T}}_i(t_k^i) PB\Xi B^{\text{T}}P\xi_i(t_k^i)+\frac{1}{\varphi_i}\eta_i(t) \leq \\ &\qquad\kappa \gamma e^{\text{T}}_i(t)PBB^{\text{T}}Pe_i(t) \end{split} $$ (41) $$ \begin{split} \Vert K e_i(t) \Vert \geq \sqrt{\frac{ \omega_i\xi^{\text{T}}_i(t_k^i) PB\Xi B^{\text{T}}P\xi_i(t_k^i)+\frac{1}{\varphi_i}\eta_i(t)}{\kappa \gamma}} \end{split} $$ (42) 因此, 可以得到
$$ \begin{split} &\frac{\Upsilon}{\Vert A \Vert}(\text{e}^{\Vert A \Vert(t-t_k^i)}-1) \geq\\&\qquad \sqrt{\frac{ \omega_i\xi^{\text{T}}_i(t_k^i) PB\Xi B^{\text{T}}P\xi_i(t_k^i)+\frac{1}{\varphi_i}\eta_i(t)}{\kappa \gamma}} \end{split} $$ (43) $$ \begin{split} &\text{e}^{\Vert A \Vert(t-t_k^i)}-1\geq \\&\qquad\frac{\Vert A \Vert}{\Upsilon} \sqrt{\frac{ \omega_i\xi^{\text{T}}_i(t_k^i) PB\Xi B^{\text{T}}P\xi_i(t_k^i)+\frac{1}{\varphi_i}\eta_i(t)}{\kappa \gamma}} \end{split} $$ (44) $$ \begin{split} & t-t_k^i\geq \frac {1}{\Vert A \Vert} \;\times\\&\ln \left( \frac{\Vert A \Vert}{\Upsilon} {\sqrt{\frac{ \omega_i\xi^{\text{T}}_i(t_k^i) PB\Xi B^{\text{T}}P\xi_i(t_k^i)+\frac{1}{\varphi_i}\eta_i(t)}{\kappa \gamma}}+1}\right) \end{split} $$ (45) □ 至此, 证明了所设计的事件触发通信机制可以排除芝诺现象. 在第3节中, 将在无人船集群系统中验证所提方法的有效性.
注 7. 相较于文献[41], 通过引入$ \Gamma $和$ \vartheta_i \text{e}^{-\digamma_it} $, 所设计的动态函数$ \eta_i(t) $衰减速率更慢, 可进一步延长触发间隔, 并且设计的动态函数$ \eta_i(t) $以指数形式收敛至零.
注 8. 值得一提的是, 当智能体发生较为严重的执行器故障时, $ \kappa $将变大且$ \omega_i $将变小, 触发条件$ f_i(e_i(t),\;\xi_i(t_k^i),\;\eta_i(t))\leq0 $会更容易满足, 智能体间的通信会更频繁.
3. 仿真结果
在本节中, 将通过仿真验证所提方法可以满足预期要求. 随着无人船的大量应用, 无人船集群控制得到了广泛关注. 本文将文献[42]的无人船模型进行线性化, 该模型由六个状态组成: 纵荡位置、横荡位置、偏航角、纵荡速度、横荡速度和偏航速率, 该多智能体系统由五个相同的智能体构成. 系统矩阵$ A $和$ B $如下式
$$ A = \left[ { \begin{array}{cc} A1 & A2 \\ A3 & A4 \end{array} } \right] $$ $$ \begin{split} B=\left[ {\begin{array}{ccc} 0 & 0 & 0\\ 0 & 0 & 0\\ 0 & 0& 0\\ 0.250\,0 & 0 & 0\\ 0 & 0.339\,6 & -0.377\,4 \\ 0 & 0.377\,4 &\; \;\;2.641\,5 \end{array} } \right] \end{split} $$ (46) 其中
$$ \begin{split} A1 = \left[ { \begin{array}{ccc} 0 & 0 & 0 \\ 0 & 0 & 0\\ 0 & 0 & 0 \end{array} } \right],\;\; A2 = \left[ { \begin{array}{ccc} 1 &0& 0\\ 0&1&0\\ 0&0&1 \end{array} } \right] \end{split} $$ $$ \begin{split} A3 = \left[ { \begin{array}{ccc} -0.020\,0 & 0 & 0 \\ 0 &-0.034\,0&\;\;\;0.037\,7\\ 0&\;\;\;0.037\,7&-0.264\,2 \end{array} } \right] \end{split}$$ $$ \begin{split} A4 = \left[ { \begin{array}{ccc} -0.125\,0&0&0 \\ 0& -0.135\,8&\;\;\;0.169\,8\\ 0&\;\;\;0.150\,9&-1.188\,7 \end{array} } \right] \end{split} $$ 考虑如图1所示的由五个智能体组成的多智能体系统, 系统的拉普拉斯矩阵$ L $如下
$$ L = \left[ \begin{array}{*{20}{r}} 3 & -1 & 0 & -1& -1 \\ -1 & 1 & 0 & 0& 0 \\ 0 & 0 & 1 & -1& 0\\ -1& 0 & -1 &3 & -1\\ -1& 0 & 0 & -1 &2 \end{array} \right]$$ (47) 多智能体系统的故障参数$ \rho_{ij} $为
$$ \begin{split} \Lambda_1 = \text{diag}\{ 0.66,\;0.68,\;0.74 \}\\ \nonumber \Lambda_2 = \text{diag}\{0.62,\;0.76,\;0.69 \}\\ \Lambda_3 = \text{diag}\{ 0.73,\;0.56 ,\;0.72 \}\\ \Lambda_4 = \text{diag}\{0.67,\;0.33,\;0.64 \}\\ \Lambda_5 = \text{diag}\{0.55,\;0.45,\;0.67 \} \end{split} $$ 动态变量$ \eta_i(t) $的初值为
$$ \begin{split} \eta_i(t) = (0.055,\;0.043,\;0.034,\;0.047 ,\;0.051) \nonumber \end{split} $$ 通过求解线性矩阵不等式(31), 可以得到如下矩阵
$$ \begin{split} P = \left[ { \begin{array}{cc} P1 & P2 \\ P3 & P4 \end{array} } \right] \end{split} $$ (48) 其中
$$ \begin{split} & P1 = \left[ {\begin{array}{ccc} 0.003\,5 & 0 & 0 \\ 0 & \;\;\;0.003\,7 & -0.001\,9 \\ 0 & -0.001\,9 &\;\;\; 0.015\,8 \\ \end{array} } \right] \\ & P2 = \left[ {\begin{array}{ccc} 0.003\,7 & 0 & 0\\ 0 & \;\;\;0.002\,9 & -0.002\,0\\ 0 & -0.001\,2 &\;\;\; 0.015\,2\\ \end{array} } \right] \\ & P3 = \left[ {\begin{array}{ccc} 0.003\,7 & 0 & 0 \\ 0 &\;\;\; 0.002\,9 & -0.001\,2 \\ 0 & -0.002\,0 &\;\;\; 0.015\,2 \end{array} } \right] \\ & P4 = \left[ {\begin{array}{ccc} 0.146\,2 & 0 & 0\\ 0 & \;\;\;0.110\,7 & -0.000\,4\\ 0 & -0.000\,4 & \;\;\;0.122\,9 \end{array} } \right] \end{split} $$ 如图2 ~ 图7所示, 当多智能体系统发生执行器故障时, 自适应容错控制器可以保证多智能体系统的稳定性. 图8表示编队中无人船的纵荡位置和横荡位置, 仿真结果表明本文的容错控制方法可以使无人船集群很好地实现一致性.
从图9和图10可以看出, 如果智能体之间的状态差异过大, 智能体间的通信会更加频繁. 因此在实现一致性的过程中, 智能体之间需要频繁的通信. 当多智能体系统逐渐达成一致时, 智能体之间的通信频率不断降低, 所设计的事件触发通信机制减少了不必要的通信. 为说明本文所提事件触发机制具有更好的应用效果, 将文献[25]所提事件触发控制方法应用于本例中的无人船编队系统, 相应的事件触发时刻与触发间隔仿真效果如图10所示. 对比图9和图10可以看出, 本文所提方法能够更好地节省通信资源.
图 9 基于触发通信协议(8)的触发时刻和触发间隔, (a)、(b)、(c)、(d) 和(e)分别表示智能体1、智能体2、智能体3、智能体4和智能体5的触发时刻和触发间隔Fig. 9 Triggering moment and triggering interval for each agent by communication protocol (8), where (a), (b), (c), (d) and (e) indicate the triggering moment and triggering interval of the agent 1, agent 2, agent 3, agent 4 and agent 5 respectively
图 10 基于文献[25]的触发时刻和触发间隔, (a)、(b)、(c)、(d) 和(e)分别表示智能体1、智能体2、智能体3、智能体4和智能体5的触发时刻和触发间隔Fig. 10 Triggering moment and triggering interval for each agent in reference [25], where (a), (b), (c), (d) and (e) indicate the triggering moment and triggering interval of the agent 1, agent 2, agent 3, agent 4 and agent 5 respectively自适应参数$ \hat{\rho}_{ij} $和动态变量$ \eta_i(t) $如图11和图12所示. 当多智能体系统趋于一致时, 自适应参数$ \hat{\rho}_{ij} $趋于稳定, 动态变量$ \eta_i(t) $渐进收敛. 当智能体发生执行器故障时, 自适应参数$ \hat{\rho}_{ij} $可以补偿故障, 因此多智能体系统具有更好的控制性能.
仿真结果表明, 在执行器发生故障的情况下, 控制器(6)可以保证多智能体系统的稳定性和一致性. 此外, 触发通信协议(8)可以有效地降低智能体之间的通信频率.
4. 结束语
本文研究基于事件触发通信机制的多智能体系统可靠控制问题. 通过设计新的自适应事件触发通信条件, 可以有效地减少智能体之间的通信频率. 基于该触发通信机制, 所设计的控制器避免了对于邻居实时状态信息的依赖. 此外, 该控制器具有通信间隔内时变的特点, 能保证在执行器失效故障和非连续通信条件下多智能体系统的一致性和稳定性. 最后, 通过仿真验证所提出的方法可以满足设计要求. 由于考虑的故障类型为线性多智能体系统执行器失效故障, 具有一定的局限性. 未来将面向执行器卡死、偏移故障和异构多智能体系统开展基于事件触发通信的容错控制研究.
-
图 1 多尺度残差重构网络(MSRNet), s-Dconv表示扩张卷积, $s$ = 1, 2, 4
Fig. 1 Mult-scale residuce recontruction network, s-Dconv denotes s-dilate convolution, here $s$ = 1, 2 and 4
图 4 重构图像块修正前后误差的比较(Barbara图像)
Fig. 4 The comparison of the error of the reconstructed image block before and after refined (Barbara)
图 7 比较几种算法的重构性能(第1行到第3行采样率MR = 0.25, 0.10, 0.04)
Fig. 7 Comparison of reconstruction performance of various algotithms (MR = 0.25, 0.10, 0.04)
图 8 修正前后的重构图像对比(MR = 0.25)
Fig. 8 The comparison of reconstructed images before and after refined (MR = 0.25)
表 1 6幅测试图像在不同算法不同采样率下的PSNR
Table 1 PSNR values in dB for six testing images by different algorithms at different measurement rates
图像 算法 MR = 0.25 MR = 0.10 MR = 0.04 MR = 0.01 w/o BM3D w/BM3D w/o BM3D w/BM3D w/o BM3D w/BM3D w/o BM3D w/BM3D Barbara TVAL3 24.19 24.20 21.88 22.21 18.98 18.98 11.94 11.96 NLR-CS 28.01 28.00 14.80 14.84 11.08 11.56 5.50 5.86 D-AMP 25.08 25.96 21.23 21.23 16.37 16.37 5.48 5.48 ReconNet 23.25 23.52 21.89 22.50 20.38 21.02 18.61 19.08 DR$ ^2 $-Net 25.77 25.99 22.69 22.82 20.70 21.30 18.65 19.10 MSRNet 26.69 26.91 23.04 23.06 21.01 21.28 18.60 18.90 Boats TVAL3 28.81 28.81 23.86 23.86 19.20 19.20 11.86 11.88 NLR-CS 29.11 29.27 14.82 14.86 10.76 11.21 5.38 5.72 D-AMP 29.26 29.26 21.95 21.95 16.01 16.01 5.34 5.34 ReconNet 27.30 27.35 24.15 24.10 21.36 21.62 18.49 18.83 DR$ ^2 $-Net 30.09 30.30 25.58 25.90 22.11 22.50 18.67 18.95 MSRNet 30.74 30.93 26.32 26.50 22.58 22.79 18.65 18.88 Flinstones TVAL3 24.05 24.07 18.88 18.92 14.88 14.91 9.75 9.77 NLR-CS 22.43 22.56 12.18 12.21 8.96 9.29 4.45 4.77 D-AMP 25.02 24.45 16.94 16.82 12.93 13.09 4.33 4.34 ReconNet 22.45 22.59 18.92 19.18 16.30 16.56 13.96 14.08 DR$ ^2 $-Net 26.19 26.77 21.09 21.46 16.93 17.05 14.01 14.18 MSRNet 26.67 26.89 21.72 21.81 17.28 17.40 13.83 14.10 Lena TVAL3 28.67 28.71 24.16 24.18 19.46 19.47 11.87 11.89 NLR-CS 29.39 29.67 15.30 15.33 11.61 11.99 5.95 6.27 D-AMP 28.00 27.41 22.51 22.47 16.52 16.86 5.73 5.96 ReconNet 26.54 26.53 23.83 24.47 21.28 21.82 17.87 18.05 DR$ ^2 $-Net 29.42 29.63 25.39 25.77 22.13 22.73 17.97 18.40 MSRNet 30.21 30.37 26.28 26.41 22.76 23.06 18.06 18.35 Monarch TVAL3 27.77 27.77 21.16 21.16 16.73 16.73 11.09 11.11 NLR-CS 25.91 26.06 14.59 14.67 11.62 11.97 6.38 6.71 D-AMP 26.39 26.55 19.00 19.00 14.57 14.57 6.20 6.20 ReconNet 24.31 25.06 21.10 21.51 18.19 18.32 15.39 15.49 DR$ ^2 $-Net 27.95 28.31 23.10 23.56 18.93 19.23 15.33 15.50 MSRNet 28.90 29.04 23.98 24.17 19.26 19.48 15.41 15.61 Peppers TVAL3 29.62 29.65 22.64 22.65 18.21 18.22 11.35 11.36 NLR-CS 28.89 29.25 14.93 14.99 11.39 11.80 5.77 6.10 D-AMP 29.84 28.58 21.39 21.37 16.13 16.46 5.79 5.85 ReconNet 24.77 25.16 22.15 22.67 19.56 20.00 16.82 16.96 DR$ ^2 $-Net 28.49 29.10 23.73 24.28 20.32 20.78 16.90 17.11 MSRNet 29.51 29.86 24.91 25.18 20.90 21.16 17.10 17.33 平均PSNR TVAL3 27.84 27.87 22.84 22.86 18.39 18.40 11.31 11.34 NLR-CS 28.05 28.19 14.19 14.22 10.58 10.98 5.30 5.62 D-AMP 28.17 27.67 21.14 21.09 15.49 15.67 5.19 5.23 ReconNet 25.54 25.92 22.68 23.23 19.99 20.44 17.27 17.55 DR$ ^2 $-Net 28.66 29.06 24.32 24.71 20.80 21.29 17.44 17.80 MSRNet 29.48 29.67 25.16 25.38 21.41 21.68 17.54 17.82 
下载: 导出CSV
表 2 不同算法下11幅测试图像平均SSIM
Table 2 Mean SSIM values for 11 testing images by different algorithms
算法 MR = 0.01 MR = 0.04 MR = 0.10 MR = 0.25 ReconNet 0.4083 0.5266 0.6416 0.7579 DR$ ^2 $-Net 0.4291 0.5804 0.7174 0.8431 MSRNet 0.4535 0.6167 0.7598 0.8698
下载: 导出CSV
表 3 MSRNet重构图像修正后11幅测试图像的PSNR (dB)和SSIM
Table 3 The PSNR (dB) and SSIM of 11 test images of refined MSRNet reconstruction
图像 MR = 0.25 MR = 0.10 MR = 0.04 PSNR SSIM PSNR SSIM PSNR SSIM Monarch 29.74 0.9189 24.40 0.8078 19.62 0.6400 Parrots 30.13 0.9054 25.15 0.8240 22.06 0.7311 Barbara 27.53 0.8553 23.28 0.6630 21.39 0.5524 Boats 31.63 0.8999 26.73 0.7753 22.86 0.6355 C-man 27.17 0.8433 23.33 0.7400 20.51 0.6378 Fingerprint 28.75 0.9280 23.18 0.7777 18.81 0.5212 Flinstones 27.77 0.8529 22.14 0.7114 17.46 0.4824 Foreman 35.85 0.9297 31.71 0.8740 26.97 0.7939 House 34.15 0.8891 29.55 0.8196 25.60 0.7443 Lena 30.95 0.9019 26.68 0.7965 23.16 0.6882 Peppers 30.67 0.8898 25.43 0.7811 21.28 0.6382 平均值 30.39 0.8922 25.59 0.7791 21.79 0.6423
下载: 导出CSV
表 4 不同卷积方式在图 5的测试集中重构图像的平均PSNR (dB)
Table 4 Mean PSNR in dB for testing set in Fig. 5 by different convolution
卷积形式 MR = 0.01 MR = 0.04 MR = 0.10 MR = 0.25 普通卷积 17.50 21.23 24.79 29.05 扩张卷积 17.54 21.41 25.16 29.48
下载: 导出CSV
表 5 不同卷积方式在BSD500测试集中重构图像平均PSNR (dB)
Table 5 Mean PSNR in dB for BSD500 testing set by different convolution
算法 MR = 0.01 MR = 0.04 MR = 0.10 MR = 0.25 普通卷积 19.34 22.14 24.48 27.78 扩张卷积 19.35 22.25 24.73 27.93
下载: 导出CSV
表 6 重构一幅$ 256 \times 256 $图像的运行时间(s)
Table 6 Time (in seconds) for reconstruction a single $ 256 \times 256 $ image
算法 MR = 0.01 (CPU/GPU) MR = 0.04 (CPU/GPU) MR = 0.10 (CPU/GPU) MR = 0.25 (CPU/GPU) ReconNet 0.5363/0.0107 0.5369/0.0100 0.5366/0.0101 0.5361/0.0105 DR$ ^2 $-Net 1.2039/0.0317 1.2064/0.0317 1.2096/0.0314 1.2176/0.0326 MSRNet 0.4884/0.0121 0.5172/0.0124 0.5152/0.0117 0.5206/0.0126
下载: 导出CSV
表 7 不同算法在BSD500测试集的平均PSNR (dB)和平均SSIM
Table 7 Mean PSNR in dB and SSIM values for BSD500 testing images by different algorithms
模型 MR = 0.01 MR = 0.04 MR = 0.10 MR = 0.25 PSNR SSIM PSNR SSIM PSNR SSIM PSNR SSIM ReconNet 19.17 0.4247 21.40 0.5149 23.28 0.6121 25.48 0.7241 DR$ ^2 $-Net 19.34 0.4514 21.86 0.5501 24.26 0.6603 27.56 0.7961 MSRNet 19.35 0.4541 22.25 0.5696 24.73 0.6837 27.93 0.8121
下载: 导出CSV
表 8 比较ReconNet、DR$ ^2 $-Net和MSRNet三种算法对高斯噪声的鲁棒性(图 5中11幅测试图像)
Table 8 Comparison of robustness to Gaussian noise among of ReconNet, DR$ ^2 $-Net, MSRNet (11 testing images in Fig. 5)
模型 MR = 0.25 MR = 0.10 $ \sigma $ = 0.01 $ \sigma $ = 0.05 $ \sigma $ = 0.10 $ \sigma $ = 0.25 $ \sigma $ = 0.01 $ \sigma $ = 0.05 $ \sigma $ = 0.10 $ \sigma $ = 0.25 ReconNet 25.44 23.81 20.81 14.15 22.63 21.64 19.54 14.17 DR$ ^2 $-Net 28.49 25.63 21.45 14.32 24.17 22.70 20.04 14.54 MSRNet 29.28 26.50 22.63 18.46 25.06 23.56 21.11 15.46
下载: 导出CSV
表 9 比较ReconNet、DR$ ^2 $-Net和MSRNet三种算法对高斯噪声的鲁棒性(BSD500数据集)
Table 9 Comparison of robustness to Gaussian noise among of ReconNet, DR$ ^2 $-Net, MSRNet (BSD500 dataset)
模型 MR = 0.25 MR = 0.10 $ \sigma $ = 0.01 $ \sigma $ = 0.05 $ \sigma $ = 0.10 $ \sigma $ = 0.25 $ \sigma $ = 0.01 $ \sigma $ = 0.05 $ \sigma $ = 0.10 $ \sigma $ = 0.25 ReconNet 25.38 22.03 20.72 14.03 23.22 22.06 19.85 14.51 DR$ ^2 $-Net 27.40 24.99 21.32 14.47 24.17 22.74 20.26 14.85 MSRNet 27.78 25.53 22.37 18.38 24.67 23.34 21.22 16.37
下载: 导出CSV
-
[1] Donoho D L. Compressed sensing. IEEE Transactions on Information Theory, 2006, 52(4):1289-1306 doi: 10.1109/TIT.2006.871582 [2] Candes E J, Justin R J, Tao T. Robust uncertainty principles:exact signal reconstruction from highly incomplete frequency information. IEEE Transactions on Information Theory, 2006, 52(2):489-509 doi: 10.1109/TIT.2005.862083 [3] Candes E J, Wakin M B. An introduction to compressive sampling. IEEE Signal Processing Magazine, 2008, 25(2):21-30 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=f038caf0ee508db10a5b5fd193679bc2 [4] 任越美, 张艳宁, 李映.压缩感知及其图像处理应用研究进展与展望.自动化学报, 2014, 40(8):1563-1575 http://www.aas.net.cn/CN/abstract/abstract18426.shtmlRen Yue-Mei, Zhang Yan-Ning, Li Ying. Advances and perspective on compressed sensing and application on image processing. Acta Automatica Sinica, 2014, 40(8):1563-1575 http://www.aas.net.cn/CN/abstract/abstract18426.shtml [5] Tropp J A, Gilbert A C. Signal recovery from random measurements via orthogonal matching pursuit. IEEE Transactions on Information Theory, 2007, 53(12):4655-4666 doi: 10.1109/TIT.2007.909108 [6] Blumensath T, Davies M E. Iterative hard thresholding for compressed sensing. Applied and Computational Harmonic Analysis, 2009, 27(3):265-274 http://d.old.wanfangdata.com.cn/OAPaper/oai_arXiv.org_0805.0510 [7] Xiao Y H, Yang J F, Yuan X M. Alternating algorithms for total variation image reconstruction from random projections. Inverse Problems and Imaging, 2012, 6(3):547-563 doi: 10.3934/ipi.2012.6.547 [8] Dong W S, Shi G M, Li X, Zhang L, Wu X L. Image reconstruction with locally adaptive sparsity and nonlocal robust regularization. Signal Processing:Image Communication, 2012, 27(10):1109-1122 doi: 10.1016/j.image.2012.09.003 [9] Eldar Y C, Kuppinger P, Bolcskei H. Block-sparse signals:uncertainty relations and efficient recovery. IEEE Transactions on Signal Processing, 2010, 58(6):3042-3054 doi: 10.1109/TSP.2010.2044837 [10] 练秋生, 陈书贞.基于混合基稀疏图像表示的压缩传感图像重构.自动化学报, 2010, 36(3):385-391 http://www.aas.net.cn/CN/abstract/abstract13678.shtmlLian Qiu-Sheng, Chen Shu-Zhen. Image reconstruction for compressed sensing based on the combined sparse image representation. Acta Automatica Sinica, 2010, 36(3):385-391 http://www.aas.net.cn/CN/abstract/abstract13678.shtml [11] 沈燕飞, 李锦涛, 朱珍民, 张勇东, 代锋.基于非局部相似模型的压缩感知图像恢复算法.自动化学报, 2015, 41(2):261-272 http://www.aas.net.cn/CN/abstract/abstract18605.shtmlShen Yan-Fei, Li Jin-Tao, Zhu Zhen-Min, Zhang Yong-Dong, Dai Feng. Image reconstruction algorithm of compressed sensing based on nonlocal similarity model. Acta Automatica Sinica, 2015, 41(2):261-272 http://www.aas.net.cn/CN/abstract/abstract18605.shtml [12] Dong C, Loy C C, He K M, Tang X O. Image super-resolution using deep convolutional networks. IEEE Transactions on Pattern Analysis & Machine Intelligence, 2016, 38(2):295-307 http://d.old.wanfangdata.com.cn/Periodical/jsjfzsjytxxxb201709007 [13] Long J, Shelhamer E, Darrell T. Fully convolutional networks for semantic segmentation. In:Proceedings of the 2015 IEEE Conference on Computer Vision and Pattern Recognition. Boston, MA, USA:IEEE, 2015. 3431-3440 [14] Zhang K, Zuo W M, Chen Y J, Meng D Y, Zhang L. Beyond a gaussian denoiser:residual learning of deep CNN for image denoising. IEEE Transactions on Image Processing, 2017, 26(7):3142-3155 doi: 10.1109/TIP.2017.2662206 [15] Mousavi A, Patel A B, Baraniuk R G. A deep learning approach to structured signal recovery. In:Proceedings of the 53rd Annual Allerton Conference on Communication, Control, and Computing. Monticello, IL, USA:IEEE, 2015. 1336-1343 [16] Kulkarni K, Lohit S, Turaga P, Kerviche R, Ashok A. ReconNet:non-iterative reconstruction of images from compressively sensed measurements. In:Proceedings of the 2016 IEEE Conference on Computer Vision and Pattern Recognition. Las Vegas, NV, USA:IEEE, 2016. 449-458 https://ieeexplore.ieee.org/document/7780424/ [17] He K M, Zhang X Y, Ren S Q, Sun J. Deep residual learning for image recognition. In:Proceedings of the 2016 IEEE Conference on Computer Vision and Pattern Recognition. Las Vegas, NV, USA:IEEE, 2016. 770-778 http://www.oalib.com/paper/4016438 [18] Yao H T, Dai F, Zhang D M, Ma Y K, Zhang S L, Zhang Y D, et al. DR2-Net:deep residual reconstruction network for image compressive sensing[Online], available:https://arxiv.org/pdf/1702.05743.pdf, July 6, 2017. [19] Yu F, Koltun V. Multi-scale context aggregation by dilated convolutions[Online], available:https://arxiv.org/abs/1511.07122, April 30, 2016. [20] Gan L. Block compressed sensing of natural images. In:Proceedings of the 15th International Conference on Digital Signal Processing. Cardiff, UK:IEEE, 2007. 403-406 [21] Kingma D P, Ba J L. Adam:a method for stochastic optimization. In:Proceedings of the 3rd International Conference for Learning Representations. San Diego, USA:Computer Science, 2014. [22] Krizhevsky A, Sutskever I, Hinton G E. ImageNet classification with deep convolutional neural networks. Communications of the ACM, 2017, 60(6):84-90 doi: 10.1145/3065386 [23] Dabov K, Foi A, Katkovnik V, Egiazarian K. Image denoising by sparse 3-D transform-domain collaborative filtering. IEEE Transactions on Image Processing, 2007, 16(8):2080-2095 doi: 10.1109/TIP.2007.901238 [24] Wu J S, Jiang L Y, Han X, Senhadji L, Shu H Z. Performance evaluation of wavelet scattering network in image texture classification in various color spaces. Journal of Southeast University, 2015, 31(1):46-50 http://cn.bing.com/academic/profile?id=c73d072061e32a8b8d31062f4abad153&encoded=0&v=paper_preview&mkt=zh-cn [25] Abadi M, Agarwal A, Barham P, Brevdo E, Chen Z F, Citro C, et al. TensorFlow:large-scale machine learning on heterogeneous distributed systems[Online], available:http://tensorflow.org, September 19, 2017. [26] Glorot X, Bengio Y. Understanding the difficulty of training deep feedforward neural networks. In:Proceedings of the 13th International Conference on Artificial Intelligence and Statistics. Sardinia, Italy:JMLR, 2010. 249-256 [27] Li C B, Yin W T, Jiang H, Zhang Y. An efficient augmented lagrangian method with applications to total variation minimization. Computational Optimization & Applications, 2013, 56(3):507-530 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=252968b97ae188abfb677b3e976d5e15 [28] Dong W S, Shi G M, Li X, Ma Y, Huang F. Compressive sensing via nonlocal low-rank regularization. IEEE Transactions on Image Processing, 2014, 23(8):3618-3632 doi: 10.1109/TIP.2014.2329449 [29] Metzler C A, Maleki A, Baraniuk R G. From denoising to compressed sensing. IEEE Transactions on Information Theory, 2016, 62(9):5117-5144 doi: 10.1109/TIT.2016.2556683 期刊类型引用(6)
1. 傅卫沁，刘彬，徐长瑜. 基于多传感器的工业机械臂精细化操作远程控制方法. 机械与电子. 2023(05): 62-66+71 . 
百度学术2. 张玉泽，姚立忠，罗海军. 基于量子加权最小门限单元网络的出水COD预测. 环境工程技术学报. 2023(05): 1857-1864 . 百度学术3. 李文静，李治港，乔俊飞. 基于突触巩固机制的前馈小世界神经网络设计. 自动化学报. 2023(10): 2145-2158 . 本站查看4. 陈雨点. 模糊神经网络理论下的外贸出口量预测模型. 现代商业. 2022(17): 60-62 . 百度学术5. 丁海旭 ，汤健 ，夏恒 ，乔俊飞 . 基于TS-FNN的城市固废焚烧过程MIMO被控对象建模. 控制理论与应用. 2022(08): 1529-1540 . 百度学术6. 李浩楠，刘勇. 模糊神经网络的优化及其应用. 哈尔滨理工大学学报. 2020(06): 142-149 . 百度学术其他类型引用(12)
-
下载:
计量
- 文章访问数: 2987
- HTML全文浏览量: 929
- PDF下载量: 249
- 被引次数: 18
