基于改进粒子群优化和Stackelberg博弈的武器部署

刘富樯; 刘中阳; 周伦; 皮阳军; 蒲华燕; 罗均

doi:10.16383/j.aas.c240257

基于改进粒子群优化和Stackelberg博弈的武器部署

doi: 10.16383/j.aas.c240257 cstr: 32138.14.j.aas.c240257

刘富樯^{1, 2,},
刘中阳^{1, 2,},
周伦^{1, 2,},
皮阳军^{1, 2,},
蒲华燕^{1, 2,},
罗均^{1, 2,}

1.
重庆大学机械与运载工程学院重庆 400044
2.
重庆大学高端装备机械传动全国重点实验室重庆 400044

基金项目: 国家自然科学基金(62033001), 重庆市技术创新与应用发展重点项目(CSTB2023TIAD-KPX0057) 资助

详细信息

作者简介:
刘富樯：重庆大学机械与运载工程学院副教授. 2015年获得西北工业大学博士学位. 主要研究方向为机器人系统设计, 博弈决策与容错控制. 本文通信作者. E-mail: liufq@cqu.edu.cn

刘中阳：重庆大学机械与运载工程学院硕士. 2021年获得吉林大学学士学位. 主要研究方向为智能设备决策规划. E-mail: liuzhongyang0323@163.com

周伦：重庆大学机械与运载工程学院博士研究生. 2020年获得石河子大学硕士学位. 主要研究方向为智能设备决策规划和深度强化学习. E-mail: Lunz0630@163.com

皮阳军：重庆大学机械与运载工程学院教授. 2010年获得浙江大学博士学位. 主要研究方向为分布式参数系统的控制, 智能无人系统和振动控制. E-mail: cqpp@cqu.edu.cn

蒲华燕：重庆大学机械与运载工程学院教授. 2011年获得华中科技大学博士学位. 主要研究方向为智能无人系统, 振动控制和机器人学. E-mail: phygood_2001@shu.edu.cn

罗均：重庆大学机械与运载工程学院教授. 2000年获得上海交通大学博士学位. 主要研究方向为人工智能, 传感技术, 智能无人系统和特种机器人. E-mail: luoj@cqu.edu.cn

计量
- 文章访问数: 42
- HTML全文浏览量: 34
- 被引次数: 0
出版历程
- 收稿日期: 2024-05-10
- 录用日期: 2025-02-08
- 网络出版日期: 2025-03-16

Weapon Deployment Based on Improved Particle Swarm Optimization and Stackelberg Game

LIU Fu-Qiang^{1, 2
,},
LIU Zhong-Yang^{1, 2
,},
ZHOU Lun^{1, 2
,},
PI Yang-Jun^{1, 2
,},
PU Hua-Yan^{1, 2
,},
LUO Jun^{1, 2
,}

1.
College of Mechanical and Vehicle Engineering, Chongqing University, Chongqing 400044
2.
State Key Laboratory of Mechanical Transmission for Advanced Equipment, Chongqing University, Chongqing 400044

Funds: Supported by National Natural Science Foundation of China (62033001) and Key Project of Technological Innovation and Application Development of Chongqing (CSTB2023TIAD-KPX0057)

More Information

Author Bio:
LIU Fu-Qiang　Associate professor at the College of Mechanical and Vehicle Engineering, Chongqing University. He received his Ph.D. degree from Northwestern Polytechnical University in 2015. His research interest covers robotic system design, game decision making and fault-tolerant control. Corresponding author of this paper

LIU Zhong-Yang　Master student at the College of Mechanical and Vehicle Engineering, Chongqing University. He received his bachelor degree from Jilin University in 2021. His main research interest is intelligent equipment decision planning

ZHOU Lun　Ph.D. candidate at the College of Mechanical and Vehicle Engineering, Chongqing University. He received his master degree from Shihezi University in 2023. His research interest covers intelligent equipment decision planning and deep reinforcement learning

PI Yang-Jun　Professor at the College of Mechanical and Vehicle Engineering, Chongqing University. He received his Ph.D. degree from Zhejiang University in 2010. His research interest covers control of distributed parameter systems, intelligent unmanned systems, and vibration control

PU Hua-Yan　Professor at the College of Mechanical and Vehicle Engineering, Chongqing University. She received her Ph.D. degree from Huazhong University of Science and Technology in 2011. Her research interest covers intelligent unmanned system, vibration controlling, and robotics

Luo Jun　Professor at the College of Mechanical and Vehicle Engineering, Chongqing University. He received his Ph.D. degree from Shanghai Jiao Tong University in 2000. His research interest covers artificial intelligence, sensing technology, intelligent unmanned systems, and special robotics

摘要

摘要: 为应对来袭目标的机动调整对防区防御能力的影响, 针对性设计全新的部署优化模型和求解算法. 首先, 从战术层面出发, 提出一种考虑攻防信息变化的新型武器部署模型, 该模型能够动态调整部署策略以提高防御系统的整体效能; 其次, 设计基于混沌映射机制和K均值聚类与重心法的算法初始化方案, 以应对资源紧缺和充足两种情况, 降低算法陷入局部最优的风险; 然后, 设计基于Metropolis准则的个体最优更新方法和基于Stackelberg博弈模型的全局最优更新方法用以指导种群的进化方向; 最后, 通过提供多规模场景仿真实验, 验证了新模型和所提算法的有效性, 对比实验结果表明新模型能够更准确地反映部署方案之间的差异, 所提算法在求解质量与收敛性方面均有显著提高.
- 武器部署 /
- Stackelberg博弈 /
- 粒子群算法 /
- K均值聚类 /
- 重心法
Abstract: In order to cope with the impact of incoming target＇s maneuvering adjustment on the defensive capability of the defense zone, a new deployment optimization model and a solution algorithm are designed. Firstly, at the tactical level, a new weapon deployment model is proposed that takes into account changes in offensive and defensive information, allowing for dynamic adjustment of the deployment strategy to improve the overall effectiveness of the defense system; Secondly, an initialization scheme for the algorithm, based on chaotic mapping mechanism and K-means clustering and center-of-gravity method, is designed to cope with both limited and ample resources, reducing the risk of the algorithm falling into local optima; Then, an individual optimal renewal method based on the Metropolis criterion and a global optimal renewal method based on the Stackelberg game model are designed to guide the evolutionary directions of the population; Finally, the effectiveness of the new model and the proposed algorithm is verified through multi-scale scenario simulation experiments. The results of the comparative experiments show that the new model reflects deployment scheme differences more accurately and the proposed algorithm significantly improves the solution quality and convergence.
- Weapon deployment /
- Stackelberg game /
- particle swarm optimization (PSO) /
- K-means clustering /
- centre-of-gravity method

HTML全文

近年来随着脑机交互、元宇宙概念热度的持续增加, 作为主要的脑机交互技术, 脑电技术备受关注. 脑电采集设备主要分为侵入式和非侵入式两类. 侵入式设备通过手术将微电极阵列植入大脑内部, 能够提供高精度、高信噪比的脑电信号, 但也伴随着手术风险、脑区限制、感染风险和高成本等问题^[1]. 相比之下, 非侵入式设备则通过将电极放置在头皮表面, 避免了手术带来的各种风险和不适, 因此在临床和科研中使用更多.

临床上相比于磁共振和功能近红外等设备, 非侵入式脑电采集设备具有操作简便、安全性高、成本较低、信号实时等优势, 不需要复杂的医疗操作经验就可以实现脑电信号的采集, 显著提高了用户的接受度和使用便利性. 国外脑电行业发展较早, 大部分厂商的非侵入脑电设备体积、重量及操作流程都无法照顾便携可穿戴, 尽管近年来本领域陆续有便携非侵入脑电设备推出, 但如何让受试者更加方便、快速佩戴设备, 同时兼顾脑电信号采集的质量一直是不断优化的方向.

本团队自主研制的非侵入穿戴脑电设备 (以下简称SignBrain设备, 型号P)在对实验人员培训后, 可在3分钟之内完成佩戴及使用 (设备快速佩戴的展示视频链接: https://pan.cstcloud.cn/s/uLpXYelIRbI). SignBrain设备采用爪形结构设计, 前端装配有电极的万向抱紧件 (图5)可根据不同人群头型的大小, 自动调节与头皮的接触角度, 减少佩戴设备调整接触阻抗的时间, 提高实验效率. SignBrain设备采用盐水海绵作为导电介质, 提高使用的便捷性, 佩戴前后无需洗头, 即戴即用. 设备采集的脑电数据通过蓝牙发送至上位机软件, 在线显示脑电波形及实时接触阻抗, 并可随时离线存储或上传至远端服务器进行查看.

SignBrain设备参照JJG1043-2008^[2] (脑电图机检定规程)和GB 9706.226-2021^[3] (医用电气设备第2-26部分: 脑电图机的基本安全和基本性能专用要求)等标准, 已通过国家食品药品监督管理局天津医疗器械质量监督检测中心性能指标检测和安全性检测 (图1), 信号噪声≤1 uV, 共模抑制比为111 dB, 设备能够有效抑制共模噪声和干扰, 提高信号的清晰度和稳定性.

图 1 SignBrain设备性能检测报告

Fig. 1 SignBrain device performance test report

下载: 全尺寸图片幻灯片

设备已在多家医院及相关单位完成小批量的实验测试, 并通过两个脑机交互范式, 高频视觉诱发和想象写字实验验证设备信号的高精度和实时性, 具体实验及数据将在下文中详细介绍.

1. SignBrain穿戴式脑电采集设备

1.1 结构介绍

SignBrain设备结构由主干、爪形、万向抱紧件、电极、海绵等配件构成, 选用合成的改性尼龙材料一体注塑加工而成. 整机重量约为160 g (如图7左图所示的佩戴效果).

1.1.1 主干结构

主干结构 (图2)是双侧爪形及盖板等配件的支撑结构, 提供设备在佩戴中的稳定性. 主干结构的弓形弯折角度, 决定了设备电极与头皮的贴合程度, 经过不断优化主干弯折角度, 进而调整万向抱紧件与头皮接触的压力, 使设备可在长时间佩戴使用时保持良好的贴合度及抗疲劳度. 主干左右侧分别设计有独立的舱室, 左侧舱室用于放置设备电路及射频模块, 右侧为电池舱室.

图 2 设备主干结构

Fig. 2 Device backbone structure

下载: 全尺寸图片幻灯片

1.1.2 爪形结构

SignBrain设备设计有双侧爪形, 共计有两个爪形支体. 单侧爪形设计有8个有效通道, 每一个通道占用爪形结构的一条分支. 总共有16个采集通道、1个参考通道、1个接地通道, 所有通道在爪形结构上设计的位置符合国际10-20标准导联排布^[4] (如图3所示).

图 3 16通道SignBrain设备导联位置排布图

Fig. 3 Layout of the electrodes for 16-channel SignBrain device

下载: 全尺寸图片幻灯片

爪形结构的优势在于佩戴时每一个分支通道可以穿过头发, 不同于脑电帽盖住头发的问题, 即便面对头发较多的女性, 稍微抬起所在分支, 使用棉签拨开头发即可快速降低阻抗, 减少被试在实验前的等待时间.

左侧爪形 (图4(a))包含O1、P7、T7、C3、AF3、F3、F7、FC5、GND, 其中GND电极为偏置电极, 设计在前额叶.

图 4 设备爪形结构

Fig. 4 Claw structure of the device

下载: 全尺寸图片幻灯片

右侧爪形 (图4(b))包含O2、P8、T8、C4、AF4、F4、F8、FC6、REF, 其中REF电极为中心参考电极, 设计在Cz的位置.

SignBrain设备16通道中前额叶脑区包含AF3, AF4, F3, F4, FC5, FC6, F7, F8, 可能涉及情绪调节、注意力和执行功能; 颞叶包含T7, T8可能涉及听觉处理和语言; 枕叶和顶叶有P7、P8、O1、O2, 可能涉及视觉处理.

1.1.3 万向抱紧件

SignBrain设备每一个爪形前端都设计有万向抱紧件, 该部件是设备爪形与头皮接触的重要组件

万向抱紧件 (图5)采用碗状结构, 内部结构为倒锥形, 底部可放置盐水电极及柱形海绵, 盐水电极外侧留有盐水补充口, 及时补充盐水能减少长时间实验带来的基线漂移. 顶部设计有与爪形分支末端锁紧的球形结构, 使万向抱紧件具有30°范围内的万向活动空间.

图 5 万向轴抱紧件

Fig. 5 Universal shaft

下载: 全尺寸图片幻灯片

万向抱紧件内部的盐水电极采用粉末烧结工艺制造^[5]而成, 该电极具有高导电性, 保证信号采集的准确性和稳定性. 粉末烧结材料具有良好的生物相容性, 能减少与皮肤接触时的过敏反应, 适合长期佩戴. 烧结过程可以增强电极的耐腐蚀性, 延长其使用寿命. 设备电极也可以使用常规的Ag/AgCl电极.

万向抱紧件是SignBrain设备在结构上的主要创新之一, 相比传统的脑电帽或爪形结构脑电产品, 电极通过海绵 (或其他导电载体)与头皮的接触角度是固定得, 对于不同大小、形状的头型可能会存在不完全适配. 脑电采集过程中, 影响信号质量的重要因素之一是电极与头皮的接触阻抗, 接触阻抗越大, 引入的噪声越多^[6]. 在本设备中, 由于采用柱形的盐水海绵作为头皮与电极的导电介质, 需与头皮垂直贴合. 因此在爪形结构末端设计有万向组件, 佩戴设备后, 可根据不同头型, 调整抱紧件的角度, 使其与头皮垂直, 保证海绵与头皮接触面积最大, 快速降低接触阻抗, 开展脑电实验.

1.1.4 海绵设计

考虑到不同性别、发量、头型的人群进行实验, 会存在个别通道抱紧件与头皮接触空位情况, 为此设计团队分别设计了短、中、长三种长度海绵 (图6), 设备默认组装短号海绵, 佩戴后若抱紧件与头皮的空位较大 (超过5 mm), 则建议更换中号长度海绵, 当抱紧件与头皮距离超过10 mm, 则建议更换长海绵. 更换海绵无需重新摘下设备, 只需抬起爪形分支, 重新插取海绵即可.

图 6 短中长三种海绵规格展示

Fig. 6 Three types of sponge: short, medium and long

下载: 全尺寸图片幻灯片

1.2 材料介绍

SignBrain设备在实际使用过程中, 需要对设备结构进行频繁的弯折操作 (包括佩戴过程、不同头型的挤压等), 在材料上需要考虑韧性, 没有足够韧性的材料, 主干结构无法夹住头部, 会在实验过程中滑落, 疲劳度较差的材料在几次实验后, 也将失去会弹性, 降低与头部的夹持力度. 我们在结构设计过程以及小批量生产之前进行了大量的材料选型比较, 包括尼龙、PE、ABS等材料, 综合比较发现改性尼龙可以获得好的韧性和回弹效果.

因此, 设备主体材料选用合成的改性尼龙材质, 该材质具有很好的韧性, 能够承受反复弯折而不易断裂或损坏, 大量反复弯折具有良好的抗疲劳性, 能够在长时间使用后保持性能稳定, 并且尼龙具有较好的耐磨损性能, 在一定温度范围内保持其机械性能, 对环境温度变化具有一定的适应性, 适合在摩擦和磨损较大的环境中使用. 尼龙还是一种轻质材料, 非常适合可穿戴^[7].

1.3 硬件设计方案

SignBrain设备可以采集16通道的脑电信号, 支持有线/无线蓝牙方式实时传输脑电数据 (图7), 支持在线阻抗检测及提醒, Marker同步记录, 可连续工作8个小时以上.

图 7 SignBrain设备效果图

Fig. 7 SignBrain device effect picture

下载: 全尺寸图片幻灯片

系统内部设计前级滤波、ADC、低功耗MCU等电路^[8−12]. 信号流程图如图8所示. 前端脑电信号经过前置滤波电路滤波, 由内部集成的两片ADS1299对其信号进行放大转换, 得到的数字量经过低功耗MCU压缩^[13], 通过低功耗蓝牙或USB传输至软件系统.MCU实时计算各通道阻抗, 超过软件阻抗阈值后, 在每个万向抱紧件上配置有警示灯实时提醒主试, 被试各通道电极的接触情况.

图 8 信号流程图

Fig. 8 Signal flow illustration

下载: 全尺寸图片幻灯片

SignBrain设备在开启蓝牙传输时, 最大采样率支持976 S/s, 采用USB传输最大采样率支持2 kS/s.

1.4 接口设计

SignBrain设备硬件上可通过Type-C连接外设控制器, 记录实验过程中不同trial的关键信息, 关键信息标签位与该时刻16通道脑电数据组合, 通过无线蓝牙传输至软件系统, 无线延时小于10 ms. 同样上位机软件系统也可通过USB或者UART等接口与外设控制器进行连接, 完成Marker记录. 对于软件系统中集成的实验范式, 可在软件中直接记录.

1.5 核心参数国内外对比

相比于32/64/128导联脑电设备, SignBrain爪形设备并不针对高通量脑电采集. 在设计之初, 设备主要考虑是能更加轻便、佩戴更加舒适, 使其适合长时间使用. SignBrain设备采用盐水作为导电介质, 实验前后无需洗头, 在临床门诊或高效采集中, 可穿戴、轻便、快速使用变得更为重要. 表1列举了美国Emotiv、美国CGX设备、奥地利g.tec g.Nautilus与SignBrain设备的对比. 其中美国CGX设备、美国Emotiv为类似可穿戴爪形结构设计的脑电设备, 奥地利g.tec g.Nautilus为脑电帽结构脑电设备.

表 1 SignBrain与其他厂商便携脑电设备的技术指标对比

Table 1 Comparison of technical specificities between SignBrain and the existing portable EEG devices

对比参数	SignBrain	Emotiv^[14]	美国 CGX^[15]	奥地利 g.tec^[16]
通道数量	16	14	20	16
采样频率 (Hz)	≥976	256	500	500
A/D位数	24	16	24	24
A/D分辨率 (μV)	0.53	0.51	—	0.53
阻抗	> 500 MΩ	—	—	> 100 MΩ
共模抑制比 (dB)	≥ 111	—	—	—
噪声 (uV)	≤ 1 μVpp	—	< 1 μVrms	< 0.6 μVrms
结构设计	爪形	爪形	爪形	脑电帽
重量	160 g	170 g	526 g	140 g

下载: 导出CSV

| 显示表格

可以发现, SignBrain的输入阻抗明显高于其他三个厂商, 对于微弱的脑电信号, 输入阻抗越高, 有助于更好地捕捉微弱信号, 数据采集更加准确, 并且可以减少外部噪声对信号的影响, 从而提高信号的信噪比. 同时SignBrain具有比其他厂商更高的采样率, 高采样率能够捕捉到快速变化的脑电信号, 如尖峰和快速振荡. 这对于分析脑电图的动态特征非常重要, 高采样率提供更好的时间分辨率, 在时间上相近的脑电活动可以被分开, 从而更好地分析不同事件相关的脑电活动. 佩戴方式上, 设备具有除脑电帽结构外, 比其他厂商更轻的重量, 对于快速应用的脑电场景中非常重要.

为进一步实际对比信号质量上差异, 在使用相同被试同一天进行静息态闭眼脑电信号采集, 采用Matlab对信号完成0.5 ~ 40 Hz带通滤波及FFT后, 得到频谱图发现SignBrain设备的底噪要好于Emotiv (图9). 采用相同处理方法, 将其与德国BP非穿戴设备的信号进行对比 (图10), 发现SignBrain设备的信号同样出色, 其中BP设备的具体型号为: BrainAmp DC^[17], 它是一款高性能的32导联脑电数据采集台式设备, 广泛应用于神经科学、心理学及临床研究. 表2为SignBrain设备与德国BP设备技术指标对比. 从指标上看, SignBrain设备的性能相比于德国BP高性能台式脑电设备并没有显著优势, 但通过频谱图上观察, SignBrain设备底噪的控制可以达到BP台式脑电设备的水平.

图 9 SignBrain与Emotiv的频谱图对比

Fig. 9 Spectrogram comparison between SignBrain and Emotiv

下载: 全尺寸图片幻灯片

图 10 SignBrain与BP的频谱图对比

Fig. 10 Spectrogram comparison between SignBrain and BP

下载: 全尺寸图片幻灯片

表 2 SignBrain设备与德国BP设备技术指标对比

Table 2 Comparison of technical specificities between SignBrain device and German BP device

对比参数	SignBrain	BrainAmp DC
通道数量	16	32
采样频率 (Hz)	976	5000
A/D分辨率 (μV)	0.53	0.1
输入动态范围μV(Vpp)	±17578	±16384
阻抗	实时阻抗	非实时阻抗
共模抑制比 (dB)	≥ 111	≥110
噪声 (uV)	≤ 1	< 1
电极材质	Ag/AgCL	Ag/AgCL

下载: 导出CSV

| 显示表格

2. SignBrain在临床医院验证

SignBrain设备已在各大医院及相关单位完成超过400例关于正常人、情绪检测、抑郁患者等数据的采集, 这些数据中包括静息态和任务态数据. 实际应用测试场景如图11所示.

图 11 实际应用测试图

Fig. 11 Practical application test

下载: 全尺寸图片幻灯片

前期主要包含两家医院, 北京大学第六医院和中国中医科学院广安门医院进行抑郁患者脑电数据的采集, 截至目前, 设备在两家医院采集的总人数为353例. 下面将简单介绍医院所采集脑电数据, 验证设备在实地不同人群的适用性及本设备所采集的脑电信号在抑郁症分类效果, 进而证明设备的易用性及应用价值.

图12为医院随机抽取5名男性被试实验照片截图, 图13为医院随机抽取5名女性被试实验照片截图. 通过视频记录发现10名被试在头发数量上也存在明显差异.

图 12 5名男性被试

Fig. 12 5 male participants

下载: 全尺寸图片幻灯片

图 13 5名女性被试

Fig. 13 5 female participants

下载: 全尺寸图片幻灯片

分别对男性女性被试相同通道的阻抗的均值进行对比, 以柱状图显示可以发现 (图14): (a) 10名被试中, 男性和女性佩戴设备后接触阻抗均小于18 kΩ; (b)从视频记录中看出女性的头发较长情况下, 接触阻抗仍小于18 kΩ. 同时我们对男性和女性被试所有通道阻抗值做双样本T检验, p值为0.93, 说明两个性别组无显著差异, 进一步证明设备接触阻抗不受性别影响.

图 14 不同性别各通道阻抗均值比较

Fig. 14 Comparison of the average impedance of each channel for different genders

下载: 全尺寸图片幻灯片

上述主要原因在于SignBrain设备采用爪形结构, 佩戴过程中, 每一个分支均可以很轻松穿过头发达到目标导联区域, 不会将整体头发压住, 这也是调节接触阻抗时间明显短于脑电帽的原因.

2.1 设备信号在抑郁分类上效果

通过抑郁自评量表SDS和具体得分人数统计, 发现353例医院的数据中有297人具有抑郁倾向, 无抑郁倾向56人.

同时, 对SignBrain设备采集的脑电数据进行整理、数据预处理, 采用深度学习DCAAN^[18][19]模型对处理后的数据集进行训练, 区分正常、轻度、中重度抑郁, 得到结果如表3所示, 发现其准确率达到70.61% (相关结果待发表), 在自有数据集实现了临床抑郁症的精准诊断.

表 3 不同模型在SignBrain数据集上进行抑郁分类效果

Table 3 The performance of different models in classifying depression on the SignBrain dataset

Model	Acc	Pre	Rec	F1-score
DAN (Long et al. , 2015)	74.36	58.92	60.19	58.91
DANN (Ganin et al. , 2016)	74.37	65.92	63.09	62.88
ADA (Haeusser et al. , 2017)	80.47	69.37	68.81	67.95
DCANN (Ours)	82.73	70.61	69.42	68.67

下载: 导出CSV

| 显示表格

上述数据集为SignBrain设备在临床医院两年内的数据积累, 其中大部分为门诊数据, 考虑到门诊问诊时间较短, 设备在研制设计初期, 初衷是可以快速完成可穿戴, 高质量完成脑电数据采集. 并且SignBrain设备结合上位机软件可以在数据采集后快速回看并上传至远端服务器, 这一特点大大减少脑电实验过程中需要科研人员协同医生完成脑电采集的时间, 有问题的实验数据也可快速提醒纠正. 实现了在临床医院, 患者无需过多等待, 医生无需过多操作, 科研人员远端回溯分析的特点.

3. 设备想象写字及高频视觉诱发验证

本节将简单介绍SignBrain设备在想象写字及高频视觉诱发实验范式上的完整测试和验证, 通过在线脑机交互实验, 验证设备在脑机交互领域上的可行性.

3.1 设备验证: 想象写字实验

3.1.1 想象写字方案设计

实验过程中, 参与者全程闭眼. 按照预设的26个字母的顺序, 在脑海中依次想象每个字母的书写轨迹. 为确保数据的精确性, 每个字母的想象被控制在约3秒的时间内. 在开始想象和结束想象这两个时间点, 参与者需要打Marker标记, 用于准确提取每个字母的脑电数据标签. 为了保持受试者的注意力和精力, 两次字母想象之间设置了1秒的休息时间. 这段时间允许大脑短暂的放松, 为下一次想象做好准备. 实验范式见图15.

图 15 想象写字范式图

Fig. 15 Experiment for imaginary writing recognition

下载: 全尺寸图片幻灯片

3.1.2 预处理及网络结构

对原始数据进行“滤波-降采样-重参考-标准化”等预处理流程. 首先提取0.5 ~ 40 Hz频段的信号, 去除50 Hz工频, 降采样到250 Hz. 其次对数据进行重参考, 消除电极位置的影响. 最终数据的维度是(16, 750, n). 维度16是信号通道的数量; 维度750 = 250*3是数据采样点的数量, 每一个字母是3 s的数据, 降采样到250 Hz; 维度n是预处理后提取了多少个频段的信号, 本次提取了n = 3个频段的信号, 分别为: Alpha波段、Beta波段以及包含多个波段的0.5 Hz ~ 40 Hz.

网络结构主要由时间维度卷积网络和域对抗全连接层构成. 首先将(16, 750, n = 3)维度的数据串联构造新维度(48, 750, 1)数据, 在此维度上完成不同视野卷积核在时间维度750上进行卷积, 其中特征提取器的损失函数见式 (1).

$$ {L}_{predict}=(1-\mathrm{e})H({\mathit{y}}_{\mathit{t}\mathit{u}\mathit{r}\mathit{e}}, \;{\mathit{y}}_{\mathit{p}\mathit{r}\mathit{e}\mathit{d}})+eH ( {\mathit{y}}_{\mathit{s}\mathit{o}\mathit{f}\mathit{t}}, \;{\mathit{y}}_{\mathit{p}\mathit{r}\mathit{e}\mathit{d}} )$$

(1)

式中$ {L}_{predict} $是最终损失, e为超参数 (强化模型对某些字母的关注), $ \left(1-\mathrm{e}\right)$为权重, $ H({\mathit{y}}_{\mathit{t}\mathit{u}\mathit{r}\mathit{e}}, \;{\mathit{y}}_{\mathit{p}\mathit{r}\mathit{e}\mathit{d}})$是真实标签与预测标签之间的交叉熵, $ H $($ {\mathit{y}}_{\mathit{s}\mathit{o}\mathit{f}\mathit{t}}, \;{\mathit{y}}_{\mathit{p}\mathit{r}\mathit{e}\mathit{d}} $)是标签平滑的向量与预测标签之间的交叉熵, $ {\mathit{y}}_{\mathit{s}\mathit{o}\mathit{f}\mathit{t}} $是进行标签平滑的向量.

加入域对抗模块^[19]将梯度反转层置于特征提取器的末端, 通过计算源域特征和目标域特征的最大均值差异, 模型可以更加关注不同域之间的特征分布差异问题, 提升模型对跨时间脑电信号的泛化性. 最大均值差异是基于再生核希尔伯特空间的距离来度量两个分布的差异. 具体损失函数见式 (2). 特征提取器朝着域不变特征的方向更新其权重, 提取同一类别不同时间的相似特征, 提升脑电信号跨时间分类.

$$ \begin{split}{L}_{MMD}=\;&D\left({M}^{s}, \;{M}^{t}\right)=\|\frac{1}{a}\sum _{i=1}^{a}H\left({M}_{i}^{s}\right)\\&-\frac{1}{b}\sum _{i=1}^{b}H\left({M}_{i}^{t}\right)\| \end{split}$$

(2)

式中H是希尔伯特空间变换; $ {M}^{s} $是源域的特征; $ {M}^{t} $是目标域的特征.

3.1.3 模型训练及结果评估

使用SignBrain设备采集800组数据, 模型以80%作为训练集, 10%作为验证集, 10%作为测试集的配比进行训练. 训练好的模型经过优化, 将其封装成独立想象测试软件, SignBrain设备实时采集的脑电想象信号蓝牙上传至软件, 经预处理, 模型计算后, 可在0.5s内输出被试者想象的字母, 在对26个字母进行整体测试发现平均准确度在86.5%, 具体结果见图16. 被试使用SignBrain设备真实想象字母照片如图17所示.

图 16 26字母分类测试结果

Fig. 16 26 letter recognition results

下载: 全尺寸图片幻灯片

图 17 被试使用设备想象字母

Fig. 17 Subjects used the device to imagine letters writing

下载: 全尺寸图片幻灯片

3.1.4 BP设备想象写字信号对比

为增加SignBrain与商用BP设备在想象写字信号结果上的对比, 采用相同范式、相同被试、相同模型训练, 使用Brain Products采集的900组数据和使用SignBrain设备采集的800组数据, 分别进行模型训练、结果如图18所示, 实验1为BP分类准确率为81.1%, 实验2为SignBrain设备的分类准确率为86.5%, 由实验1和2可见脑电信号的实际运用场景中, SignBrain设备和Brain Products产品均能达到较好的实验效果.

图 18 采用BP设备与SignBrain设备想象写字实验模型训练结果

Fig. 18 Training results of the imaginary writing recognition using BP device and SignBrain device

下载: 全尺寸图片幻灯片

3.2 设备验证: 高频视觉诱发实验

3.2.1 高频视觉诱发实验设计

高频视觉诱发实验要求测试人员佩戴好SignBrain设备, 所有电极的接触阻抗在10 kΩ左右, 屏幕上依次按照A到Z的顺序闪烁26个字母和一个删除键和一个空格键, 频率从8 Hz到26.9 Hz, 以0.7 Hz频率间隔为步长, 1.5 s标记一次, 实验3 s, 休息0.5 s. 该实验模式下测试两组. 所有测试数据经过预处理、分段后、使用eCCA^[20] (expanded Canonical Correlation Analysis)从诱发的信号片段中提取出干净的SSVEP成分, 抑制噪声和共模干扰, 最终将整个实验步骤方法集成到在线视觉诱发软件中进行在线测试.

3.2.2 高频视觉诱发实验结果

为展示设备在视觉通道上采集诱发信号的准确度, 选用O1采集的视觉诱发信号, 从8 Hz到26.9 Hz中随机选择10.8、13.6、15、20.6、22.7、24.1这些频率点进行FFT, 可发现设备完整保留了诱发信号的频率 (图19). 同时为对比视觉诱发频率与刺激频率的误差, 将上述5个频率点挑选出计算误差 (表4), 发现设备在被试视觉皮层O1通道诱发频率与刺激频率误差小于10%, 诱发频率的误差完全满足步长0.7 Hz的设计. 视觉诱发展示链接: https://pan.cstcloud.cn/s/dKHCh5IoQLU.

图 19 SignBrain设备6个频率点视觉诱发O1电极记录脑电信号FFT图

Fig. 19 Demonstrations for the six frequency evoked signals by SSVEP using SignBrain Device

下载: 全尺寸图片幻灯片

表 4 视觉皮层O1通道诱发频率与刺激频率的误差

Table 4 The error between the induced frequency of the O1 channel and the stimulus frequency

刺激频率 (Hz)	10.8	13.6	15	20.6	22.7	24.1
实际诱发频率 (Hz)	10.75	13.58	15.03	20.65	22.64	24.07
误差 (%)	7.1	2.8	4.3	7.1	8.6	4.3

下载: 导出CSV

| 显示表格

经实验测试, 整体26个字母在该频率区间进行视觉诱发均可以超过95%的准确率. 该实验范式无需在特定环境下进行, 被试佩戴好SignBrain设备即可开始视觉诱发实验, 脑电数据通过低延时的蓝牙发送至上位机软件并在极短时间内完成信号的比对分析, 最终呈现在显示器上. 由此说明, SignBrain设备在脑电信号采集中表现了较高的信噪比, 通信传输延时率低, 可长时间稳定采集脑电信号.

4. 讨论

SignBrain设备结构上设计为爪形, 易于佩戴, 在很多场景及范式上可以高效的完成脑电采集. 数据通过无线蓝牙发送至配套软件系统, 不受具体场地约束.

设备已在多家医院及单位累计采集超过400例脑电数据, 并对相应任务的数据集分析, 均达到预期实验结果. 本文介绍的26个字母想象范式是受试者在家中佩戴设备进行测试. 字母想象范式实验在侵入式脑电设备上应用较多, 由于微型电极直接接触神经元放电, 信号质量非常高, 可以很轻松取得较高的准确率. 但在非侵入设备上, 外部噪声及头皮外获取信号的影响, 会大大影响信号质量, 非侵入设备难度极大, 需要设备具有更高的分辨率, 信噪比及共模抑制比, 实时分辨的效果还取决于信号传输及处理过程的时间. SignBrain设备采集想象字母的数据, 对其分类, 26个字母中有22个字母的准确率达到80%以上, 有12个字母的准确率达到100%. 可见, SignBrain设备利用16通道电极即可完成想象写字, 具有较高的信噪比.

5. 结论

本文介绍了团队自主研发的SignBrain可穿戴脑电设备的特点和优势. SignBrain设备轻巧便携, 爪形结构设计确保电极始终贴合头皮, 信号稳定, 使用盐水电极便于在各种场景快速佩戴. 设备具有高信噪比和低噪声, 并在想象字母范式和视觉诱发实验中验证了信号质量, 为后续研究脑电活动信号变化提供了可靠技术支持. 该设备针对快速可穿戴用途设计, 同时脑电性能和安全性均通过医疗级检验, 可在临床疾病诊断、脑机交互、特殊场景的脑认知检测, 以及特殊移动场景下实现快速应用. 关于设备的开发和应用的讨论网站: www. SignBrain.cn, 欢迎访问交流.

致谢

SignBrain项目研发及应用工作的快速推进特别感谢四川省科技计划项目2022YFSY0014 (2022-2024)、中国科学院先导专项(Grant No. QYZDJ-SSW-SMC019)、科技部科技创新2030 - 脑科学与类脑智能项目(Grant No. 2021ZD0200200)支持与资助, 同时也感谢SignBrain项目团队及相关实验医院的协助.

图 1 投弹线

Fig. 1 Bombing lines

下载: 全尺寸图片幻灯片

图 2 火力重叠

Fig. 2 Firepower overlap

下载: 全尺寸图片幻灯片

图 3 初始化种群方法框架

Fig. 3 Framework for population initialization methods

下载: 全尺寸图片幻灯片

图 4 场景1进攻路径规划结果

Fig. 4 Path planning results in the scenario 1

下载: 全尺寸图片幻灯片

图 6 场景3进攻路径规划结果

Fig. 6 Path planning results in the scenario 3

下载: 全尺寸图片幻灯片

图 5 场景2进攻路径规划结果

Fig. 5 Path planning results in the scenario 2

下载: 全尺寸图片幻灯片

图 7 三种场景下防御部署方案适应度值收敛曲线

Fig. 7 Convergence curves of fitness values for defence deployment schemes in three scenarios

下载: 全尺寸图片幻灯片

图 8 两种部署模型评价指标值对比

Fig. 8 Comparison of the evaluation indicator values of two deployment models

下载: 全尺寸图片幻灯片

表 1 最优进攻路径代价对比

Table 1 Comparison of optimal attack path costs

仿真实例	对比指标	PSO	AGAPSO	SSGWO	IPSO
规模场景1	最大值	87.9157	90.1453	83.7577	81.1419
	最小值	84.7922	85.7749	79.3922	78.6324
	平均值	85.9452	87.0357	81.1712	79.4218
	均方差	6.3320	6.9340	7.0318	7.8558
规模场景2	最大值	84.1190	72.7094	72.0462	71.9816
	最小值	79.1262	70.2760	70.4387	70.1361
	平均值	82.4984	71.6797	71.1869	70.8235
	均方差	3.9597	3.6551	4.4898	3.9926
规模场景3	最大值	53.3404	52.7513	52.6991	52.6948
	最小值	51.2769	51.2769	51.2769	51.2769
	平均值	52.8909	51.9593	51.9407	51.9171
	均方差	1.3500	2.8143	2.2543	1.9502

下载: 导出CSV

表 2 各算法计算时间对比

Table 2 Comparison of computation time

算法	规模场景1	规模场景2	规模场景3
PSO	0.6864	0.8446	0.9469
AGAPSO	0.6314	0.8430	0.9308
SSGWO	0.6437	0.8236	0.9107
IPSO	0.6357	0.7909	0.8991

下载: 导出CSV

表 3 两种模型在对比模型评价指标中的差异

Table 3 Differences between two models in comparing model evaluation indicators

规模	组别	评价指标均值	平均值差	t	显著性
场景1	对比模型	0.5927	0.0288	−1.6265	0.1212
场景1	本文模型	0.6215	0.0288	−1.6265	0.1212
场景2	对比模型	0.6438	0.0088	0.3301	0.7451
场景2	本文模型	0.6350	0.0088	0.3301	0.7451
场景3	对比模型	0.7500	0.0036	−0.9334	0.3629
场景3	本文模型	0.7536	0.0036	−0.9334	0.3629

下载: 导出CSV

表 4 两种模型在本文模型评价指标中的差异

Table 4 Differences between two models in the model evaluation indicators of this paper

规模	组别	评价指标均值	平均值差	t	显著性
场景1	对比模型	75.5000	3.5126	−2.4348	0.0255*
场景1	本文模型	79.0126	3.5126	−2.4348	0.0255*
场景2	对比模型	73.4573	2.6694	2.2003	0.0411*
场景2	本文模型	76.1267	2.6694	2.2003	0.0411*
场景3	对比模型	63.1573	1.9786	−1.7206	0.1025
场景3	本文模型	65.1359	1.9786	−1.7206	0.1025

下载: 导出CSV

参考文献(22)

[1]	陈杰, 方浩, 辛斌, 邓方. 数字化陆用武器系统中的建模、优化与控制. 自动化学报, 2013, 39(7): 943−962 Chen Jie, Fang Hao, Xin Bin, Deng Fang. Modeling, optimization and control in ground-based digital weapon systems. Acta Automatica Sinica, 2013, 39(7): 943−962
[2]	周天颜, 冯小恩, 范云锋, 董诗音, 李玉庆, 金慧中. 避免防空火力过剩的地面兵力防御部署优化模型. 空天防御, 2022, 5(4): 19−23 doi: 10.3969/j.issn.2096-4641.2022.04.003 Zhou Tian-Yan, Feng Xiao-En, Fan Yun-Feng, Dong Shi-Yin, Li Yu-Qing, Jin Hui-Zhong. Optimization model of ground air defense force deployment to avoid excessive air defense firepower. Air & Space Defense, 2022, 5(4): 19−23 doi: 10.3969/j.issn.2096-4641.2022.04.003
[3]	曹伟伟. 高寒山地防空火力配系研究 [硕士学位论文], 国防科技大学, 中国, 2019. Cao Wei-Wei. Study on Air Defense Fire Deployment in Alpine Mountain [Master thesis], National University of Defense Technology, China, 2019.
[4]	王凤山, 杨志宏, 郭子曜. 基于三维实景地图的智能兵力部署研究. 舰船电子工程, 2021, 41(10): 107−111 doi: 10.3969/j.issn.1672-9730.2021.10.024 Wang Feng-Shan, Yang Zhi-Hong, Guo Zi-Yao. Research on intelligent force deployment based on 3D real map. Ship Electronic Engineering, 2021, 41(10): 107−111 doi: 10.3969/j.issn.1672-9730.2021.10.024
[5]	Bai X S, Fielbaum A, Kronmüller M, Knoedler L, Alonso-Mora J. Group-based distributed auction algorithms for multi-robot task assignment. IEEE Transactions on Automation Science and Engineering, 2023, 20(2): 1292−1303 doi: 10.1109/TASE.2022.3175040
[6]	Bai X S, Yan W S, Ge S S. Distributed task assignment for multiple robots under limited communication range. IEEE Transactions on Systems, Man, and Cybernetics: Systems, 2022, 52(7): 4259−4271 doi: 10.1109/TSMC.2021.3094190
[7]	Li Y, Zhang X, Zhao J, Yang X Z, Xi M Y. Position deployment optimization of maneuvering conventional missile based on improved whale optimization algorithm. International Journal of Aerospace Engineering, 2022, 2022: Article No. 4373879
[8]	温包谦, 王涛, 成坤, 张济众. 基于PSO-GA混合算法的末端防御兵力优化部署方法. 兵器装备工程学报, 2019, 40(11): 45−49 doi: 10.11809/bqzbgcxb2019.11.010 Wen Bao-Qian, Wang Tao, Cheng Kun, Zhang Ji-Zhong. End-defense force optimization deployment method based on PSO-GA hybrid algorithm. Journal of Ordnance Equipment Engineering, 2019, 40(11): 45−49 doi: 10.11809/bqzbgcxb2019.11.010
[9]	Sun W, Cao Z, Wang G, Song Y F, Guo X K. An optimized double-nested anti-missile force deployment based on the deep kuhn-munkres algorithm. Mathematics, 2022, 10(23): Article No. 4627 doi: 10.3390/math10234627
[10]	李珂铭, 李定主, 石昊. 基于改进NSGA-Ⅲ 的抗饱和攻击防空编队部署方法. 火力与指挥控制, 2023, 48(4): 153−160 doi: 10.3969/j.issn.1002-0640.2023.04.024 Li Ke-Ming, Li Ding-Zhu, Shi Hao. Deployment method of air defense formation against saturated attack based on improved NSGA-III. Fire Control & Command Control, 2023, 48(4): 153−160 doi: 10.3969/j.issn.1002-0640.2023.04.024
[11]	岳韶华, 何晟, 王刚, 刘伟, 陈晨. 基于改进遗传蜂群算法的武器系统优化部署问题研究. 兵器装备工程学报, 2022, 43(8): 80−86 doi: 10.11809/bqzbgcxb2022.08.012 Yue Shao-Hua, He Sheng, Wang Gang, Liu Wei, Chen Chen. Optimization of weapon system deployment based on improved genetic bee colony algorithm. Journal of Ordnance Equipment Engineering, 2022, 43(8): 80−86 doi: 10.11809/bqzbgcxb2022.08.012
[12]	张伟, 黄卫民. 基于种群分区的多策略自适应多目标粒子群算法. 自动化学报, 2022, 48(10): 2585−2599 Zhang Wei, Huang Wei-Min. Multi-strategy adaptive multi-objective particle swarm optimization algorithm based on swarm partition. Acta Automatica Sinica, 2022, 48(10): 2585−2599
[13]	刘建华, 刘国买, 杨荣华, 胡文瑜. 粒子群算法的交互性与随机性分析. 自动化学报, 2012, 38(9): 1471−1484 doi: 10.3724/SP.J.1004.2012.01471 Liu Jian-Hua, Liu Guo-Mai, Yang Rong-Hua, Hu Wen-Yu. Analysis of interactivity and randomness in particle swarm optimization. Acta Automatica Sinica, 2012, 38(9): 1471−1484 doi: 10.3724/SP.J.1004.2012.01471
[14]	刘景森, 吉宏远, 李煜. 基于改进蝙蝠算法和三次样条插值的机器人路径规划. 自动化学报, 2021, 47(7): 1710−1719 Liu Jing-Sen, Ji Hong-Yuan, Li Yu. Robot path planning based on improved bat algorithm and cubic spline interpolation. Acta Automatica Sinica, 2021, 47(7): 1710−1719
[15]	Shi K, Wu W D, Wu Z T, Jiang B P, Karimi H R. Coverage path planning for cleaning robot based on improved simulated annealing algorithm and ant colony algorithm. Signal Image and Video Processing, 2024, 18: 3275−3284 doi: 10.1007/s11760-023-02989-y
[16]	刘健. 地空导弹火力单元配置间距计算方法. 弹箭与制导学报, 2007, 27(1): 393−395 Liu Jian. A method to calculate the disposition interval of fire units under meeting shooting requests. Rockets, Missiles and Guidance, 2007, 27(1): 393−395
[17]	Xiong R, Wang S L, Huang Q, Yu C M, Fernandez C, Xiao W, et al. Improved cooperative competitive particle swarm optimization and nonlinear coefficient temperature decreasing simulated annealing-back propagation methods for state of health estimation of energy storage batteries. Energy, 2024, 292: Article No. 130594 doi: 10.1016/j.energy.2024.130594
[18]	Casorrán C, Fortz B, LabbéM, Ordóñez F. A study of general and security Stackelberg game formulations. European Journal of Operational Research, 2019, 278(3): 855−868 doi: 10.1016/j.ejor.2019.05.012
[19]	Bo L, Zhang Z, Liu Y, Yang S Q, Wang Y W, Wang Y Y, et al. Research on path planning method of solid backfilling and pushing mechanism based on adaptive genetic particle swarm optimization. Mathematics, 2024, 12(3): Article No. 479 doi: 10.3390/math12030479
[20]	Feng J X, Sun C L, Zhang J H, Du Y, Liu Z G, Ding, Y M. A UAV path planning method in three-dimensional space based on a hybrid gray wolf optimization algorithm. Electronics, 2024, 13(1): Article No. 68
[21]	Chen C, Chen J, Xin B. Hybrid optimization of dynamic deployment for networked fire control system. Journal of Systems Engineering and Electronics, 2013, 24(6): 954−3961 doi: 10.1109/JSEE.2013.00112
[22]	周宏宇, 王小刚, 单永志, 赵亚丽, 崔乃刚. 基于改进粒子群算法的飞行器协同轨迹规划. 自动化学报, 2022, 148(11): 2670−2676 Zhou Hong-Yu, Wang Xiao-Gang, Shan Yong-Zhi, Zhao Ya-Li, Cui Nai-Gang. Synergistic path planning for multiple vehicles based on an improved particle swarm optimization method. Acta Automatica Sinica, 2022, 148(11): 2670−2676