Research on Single-trial EEG Decoding-based Class Bootstrap Method for Lie Prediction
-
摘要: 基于脑电(Electroencephalogram, EEG)的谎言检测技术依赖于对事件相关电位(Event-related potential, ERP)的有效解码, 当前主要采用手工设计特征进行脑电分析. 近年来, 单试次脑电分类方法取得了长足进步, 其中端到端的脑电分类方法能够实现对脑电的自动特征提取和分类, 但在谎言检测中缺乏研究和应用, 同时存在无法在测谎场景下直接应用的问题. 本研究设计基于复合反应范式(Complex trial protocol, CTP)进行自我面孔信息识别任务的实验, 采集了18 名被试的脑电数据. 研究了不同端到端的单试次ERP分类方法在谎言检测中的应用, 同时针对单试次脑电解码方法无法直接实际应用的问题, 提出了一种类自举算法. 算法基于数据分布假设, 通过对比各类刺激图像被视为探针刺激时所训练模型的性能, 来推断真正的探针刺激. 实验结果表明, 在基于自我面孔信息的CTP的谎言预测中, 所提出的类自举法性能优于传统探针预测方法, 在仅使用少量脑电数据情况下, 可实现准确的谎言预测.Abstract: Lie detection techniques based on electroencephalogram (EEG) rely on the effective decoding of event-related potential (ERP). At present, manual design features are mainly used for EEG analysis. In recent years, the single-trial EEG classification method has made progress. End-to-end EEG classification methods can realize automatically extract features from EEG and classify, which lacks research and application in lie detection, also those methods cannot be directly applied in lie detection. In this study, we designed the autobiographical-based face recognition task based on a complex trial protocol (CTP) and the EEG of 18 subjects was collected. The application of different single-trial ERP classification methods in lie detection are studied. A class bootstrap method is proposed to solve the problem that the single-trial EEG decoding method cannot be applied to practice directly. The class bootstrap method was based on the assumption of data distribution, the probe stimulus was deduced by comparing the classification performance of classifiers that were trained when each category of stimulus images was set as probe stimuli. The experimental results show that the proposed class bootstrap method outperforms the traditional lie detection method and can accurately predict lies when only a small amount of EEG data is used.
-
近年来, 随着计算机、通信技术的快速发展, 多智能体理论成为研究热点, 并且相关控制方法被广泛应用于无人机集群、无人船编队、车辆编队等领域[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. 结束语
本文研究基于事件触发通信机制的多智能体系统可靠控制问题. 通过设计新的自适应事件触发通信条件, 可以有效地减少智能体之间的通信频率. 基于该触发通信机制, 所设计的控制器避免了对于邻居实时状态信息的依赖. 此外, 该控制器具有通信间隔内时变的特点, 能保证在执行器失效故障和非连续通信条件下多智能体系统的一致性和稳定性. 最后, 通过仿真验证所提出的方法可以满足设计要求. 由于考虑的故障类型为线性多智能体系统执行器失效故障, 具有一定的局限性. 未来将面向执行器卡死、偏移故障和异构多智能体系统开展基于事件触发通信的容错控制研究.
-
图 2 类自举法的分布假设示意图
Fig. 2 Schematic diagram of distribution hypothesis of the class bootstrap method
图 4 类自举法中不同脑电标签训练解码模型的特征可视化
Fig. 4 Feature visualization of decoding models trained with different EEG labels in class bootstrap method
表 1 单试次脑电解码实验主要参数
Table 1 Main parameters of single-trial EEG decoding experiment
方法 主要参数 代码来源 时域 + SVM 时域特征: Pz导联数据
SVM分类器: C = 100; kernel = “rbf”利用python工具包sklearn实现 小波 + SVM 连续小波变换; wavelet = “morl”;
SVM分类器: C = 100; kernel = “rbf”利用python工具包pywt、sklearn实现 CSP + LDA 空域特征: 全脑64 通道 利用python工具包pyRiemann、sklearn实现 HDCA 时间窗: 100 ms 根据文献[22]复现 MDRM xDAWN-Covariances + MDM python工具包pyRiemann复现 OCLNN Batch size = 8; optimizer = “adam”
lr = 0.0002; betas = (0.9, 0.999)
weight_decay = 0.001根据文献[24]利用Pytorch包复现 EEGNet Batch size = 64; lr = 0.001
optimizer = “adam”; betas = (0.9, 0.99)https://github.com/vlawhern/arleegmodels PLNet Batch size = 16; lr = 0.001
optimizer = “adam”; betas = (0.9, 0.99) early stop根据文献[26]利用tensorflow.keras库复现 
下载: 导出CSV
表 2 不同方法在不同训练数据量下的分类均衡精度(均值±标准差) (%)
Table 2 Balanced accuracy of different methods under different training data (mean±standard deviation) (%)
方法 模型训练中采用的训练数据量(组) 1 2 3 4 5 时域 + SVM 57.31±5.21** 59.42±5.25*** 60.51±6.03*** 61.12±6.66*** 61.49±7.08*** 小波 + SVM 58.06±5.79*** 58.80±4.97*** 59.40±5.77*** 59.12±5.53*** 59.31±6.10*** CSP + LDA 59.16±7.04* 63.06±9.14** 65.56±9.11*** 67.66±10.91*** 68.73±10.77*** HDCA 56.80±5.41** 60.89±5.82** 62.69±7.09*** 64.68±7.34*** 67.42±7.88*** MDRM 52.52±2.60*** 57.99±6.44*** 62.47±9.04*** 64.91±8.34*** 67.69±8.93*** OCLNN 55.25±4.88*** 61.32±5.64*** 65.83±7.16*** 68.28±7.40*** 70.88±8.55*** EEGNet 61.79±7.57 64.75±6.79* 68.72±7.39** 69.59±7.73** 71.71±7.80** PLNet 60.57±6.54 68.38±7.22 73.42±8.44 74.88±8.69 76.73±8.94
下载: 导出CSV
表 3 不同方法在不同数据量下的探针预测正确率 (%)
Table 3 Probe prediction accuracy of different methods under different data volume (%)
方法 使用的总数据量（组） 2 3 4 5 6 BAD B-P 61.11 66.67 72.22 72.22 77.78 P-P 61.11 72.22 77.78 77.78 77.78 类自举法 HDCA 27.78 38.89 61.11 77.78 83.33 MDRM 55.56 50.00 66.67 88.89 88.89 OCLNN 88.89 72.22 83.33 88.89 88.89 EEGNet 83.33 72.22 83.33 94.44 100.00 PLNet 88.89 72.22 83.33 94.44 100.00
下载: 导出CSV
-
[1] 李贺, 傅元, 梁静, 傅小兰. 欺骗识别的准确率:影响因素与提高途径. 心理科学, 2020, 43(5): 1204-1210 doi: 10.16719/j.cnki.1671-6981.20200525Li He, Fu Yuan, Liang Jing, Fu Xiao-Lan. The accuracy of deceptin detection: Influence factors and improving approaches. Psychological Science, 2020, 43(5): 1204-1210 doi: 10.16719/j.cnki.1671-6981.20200525 [2] Santhiya P, Chitrakala S. A survey on emotion recognition from EEG signals: Approaches, techniques and challenges. In: Proceedings of the International Conference on Vision Towards Emerging Trends in Communication and Networking (ViTECoN). Vellore, India: IEEE, 2019. 1−6 [3] Honts C R, Thurber S D, Handler M D. A comprehensive meta-analysis of the comparison question polygraph test. Applied Cognitive Psychology, 2021, 35: 411-427 doi: 10.1002/acp.3779 [4] Ramele R, Villar A J, Santos J M. EEG waveform analysis of P300 ERP with applications to brain computer interfaces. Brain Sciences, 2018, 8(11):199 doi: 10.3390/brainsci8110199 [5] Steiner G Z, Barry R J, Gonsalvez C J. Can working memory predict target-to-target interval effects in the P300? International Journal of Psychophysiology: Official Journal of the International Organization of Psychophysiology, 2013, 89(3): 399−408 [6] Rosenfeld J P, Labkovsky E, Winograd M, Lui M, Vandenboom C, Chedid E. The Complex Trial Protocol (CTP): a new, countermeasure-resistant, accurate, P300-based method for detection of concealed information. Psychophysiology, 2008, 45(6): 906-919 doi: 10.1111/j.1469-8986.2008.00708.x [7] Chang W W, Wang H, Lu Z G, Liu C. A concealed information test system based on functional brain connectivity and signal entropy of audio–visual ERP. IEEE Transactions on Cognitive and Developmental Systems, 2020, 12: 361-370 doi: 10.1109/TCDS.2020.2991359 [8] Wang H, Qi Y, Yu H, Wang Y, Liu C, Hu G, et al. RCIT: An RSVP-based concealed information test framework using EEG signals. IEEE Transactions on Cognitive and Developmental Systems, 2022, 14: 541-551 doi: 10.1109/TCDS.2021.3053455 [9] 彭丝雨, 周到, 张家琦, 王宇, 高军峰.基于互信息的脑网络及测谎研究.电子学报, 2019, 47(7): 1551-1556 doi: 10.3969/j.issn.0372-2112.2019.07.021Peng Si-Yu, Zhou Dao, Zhang Jia-Qi, Wang Yu, Gao Jun-Feng. Research on mutual information-based brain network and lie detection. Acta Electronica Sinica, 2019, 47(7): 1551-1556. doi: 10.3969/j.issn.0372-2112.2019.07.021 [10] Dodia S, Edla D R, Bablani A, Cheruku R. Lie detection using extreme learning machine: A concealed information test based on short-time Fourier transform and binary bat optimization using a novel fitness function. Computational Intelligence, 2020, 36: 637-658, DOI: 10.1111/coin.12256 [11] Baghel N, Singh D, Dutta M K, Burget R, Myska V. Truth identification from EEG signal by using convolution neural network: Lie detection. In: Proceedings of the 43rd International Conference on Telecommunications and Signal Processing (TSP). Milan, Italy: IEEE, 2020. 550−553 [12] Bablani A, Edla D R, Kupilli V, Dharavath R. Lie detection using fuzzy ensemble approach with novel defuzzification method for classification of EEG signals. IEEE Transactions on Instrumentation and Measurement, 2021, 70: 1-13 [13] Javaid H, Dilawari A, Khan U G, Wajid B. EEG guided multimodal lie detection with audio-visual cues. In: Proceedings of the 2nd International Conference on Artificial Intelligence (ICAI). Islamabad, Pakistan: IEEE, 2022. 71−78 [14] Zou Z X, Shi Z W, Guo Y H, Ye J P. Object detection in 20 years: A survey [Online], available: https: //doi.org/10.48550/arXiv.1905.05055, November 15, 2022 [15] Qi F, Wu W, Yu Z L, Gu Z, Wen Z, Yu T, et al. Spatiotemporal-filtering-based channel selection for single-trial EEG classification. IEEE Transactions on Cybernetics, 2021, 51: 558-567 doi: 10.1109/TCYB.2019.2963709 [16] 柳长源, 李文强, 毕晓君.基于RCNN-LSTM的脑电情感识别研究.自动化学报, 2022, 48(3): 917-925 doi: 10.16383/j.aas.c190357Liu C Y, Li W Q, Bi X J. Research on EEG emotion recognition based on RCNN-LSTM. Acta Automatica Sinica, 2022, 48(3): 917-925 doi: 10.16383/j.aas.c190357 [17] Santamaría V E, Martínez C V, Vaquerizo V F, Hornero R. EEG-inception: A novel deep convolutional neural network for assistive ERP-based brain-computer interfaces. IEEE Transactions on Neural Systems and Rehabilitation Engineering, 2020, 28: 2773-2782 doi: 10.1109/TNSRE.2020.3048106 [18] Luo J, Shi W, Lu N, Wang J, Chen H, Wang Y, et al. Improving the performance of multisubject motor imagery-based BCIs using twin cascaded softmax CNNs. Journal of Neural Engineering, 2021, 18(3), DOI: 10.1088/1741-2552/abe357 [19] Fu L, Chong W, Yang L, Hao W, Boxun Fu, Ji Y S, et al. Phase preservation neural network for electroencephalography classification in rapid serial visual presentation task. IEEE Transactions on Biomedical Engineering, 2022, 69(6):1931-1942, DOI: 10.1109/TBME.2021.3130917 [20] Martínez C R, Delorme A, Truong D, Artoni F, Kreutz D K, Sivagnanam S, et al. The open EEGLAB portal interface: High-performance computing with EEGLAB. NeuroImage, 2021, 224(6): Article No. 116778 [21] Zhang C, Eskandarian A. A survey and tutorial of EEG-based brain monitoring for driver state analysis. IEEE/CAA Journal of Automatica Sinica, 2021, 8: 1222-1242 doi: 10.1109/JAS.2020.1003450 [22] Gerson A D, Parra L C, Sajda P. Cortically coupled computer vision for rapid image search. IEEE Transactions on Neural Systems and Rehabilitation Engineering, 2006, 14(2): 174-179 doi: 10.1109/TNSRE.2006.875550 [23] Barachant A, Congedo M. A plug and play P300 BCI using information geometry [Online], available: https://arxiv.org/abs/1409.0107, October 10, 2022 [24] Lawhern V J, Solon A J, Waytowich N R. EEGNet: A compact convolutional neural network for EEG-based brain-computer interfaces. Journal of Neural Engineering, 2018, 15(5): Article No. 056013 [25] Shan H, Liu Y, Stefanov T P. A simple convolutional neural network for accurate P300 detection and character spelling in brain computer interface. In: Proceedings of the 27th International Joint Conference on Artificial Intelligence (IJCAI). Stockholm, Sweden: 2018. 1604−1610 [26] Zang B, Lin Y, Liu Z. A deep learning method for single-trial EEG classification in RSVP task based on spatiotemporal features of ERPs. Journal of Neural Engineering, 2021, 18(4).DOI: 10.1088/1741-2552/ac1610 [27] Haider S K, Jiang A, Jamshed M A, Pervaiz H B, Mumtaz S. Performance enhancement in P300 ERP single trial by machine learning adaptive denoising mechanism. IEEE Networking Letters, 2019, 1: 26-29 doi: 10.1109/LNET.2018.2883859 [28] 何章锦, 陈晨, 张金, 黎峰, 张家琦, 蓝文威等.多通道EEG信号小波相干分析的测谎研究.电子测量与仪器学报, 2020, 34(10): 156-162 doi: 10.13382/j.jemi.B2003017He Zhang-Jin, Chen Chen, Zhang Jin, Li Feng, Zhang Jia-Qi, Lan Wen-Wei, et al. Lie detection study based on wavelet coherence analysis on multi-channel EEG signals. Journal of Electronic Measurement and Instrument, 2020, 34(10): 156-162. doi: 10.13382/j.jemi.B2003017 [29] 邓晓红, 李婷, 薛超, Rosenfeld J P, 卢洋, 王莹等.基于自我参照编码的复合反应范式:有效鉴别有罪者和无罪知情者. 心理学报, 2021, 53(10): 1105-1119 doi: 10.3724/SP.J.1041.2021.01105Deng Xiao-Hong, Li Ting, Xue Chao, Rosenfeld J P, Lu Yang, Wang Ying, et al. The complex trial protocol based on self-referential encoding: Discriminating the guilty from the knowledgeable innocent. Acta Psychologica Sinica, 2021, 53(10): 1105-1119 doi: 10.3724/SP.J.1041.2021.01105 -
下载:
计量
- 文章访问数: 689
- HTML全文浏览量: 161
- PDF下载量: 143
- 被引次数: 0
