基于条件经验模式分解和串并行CNN的脑电信号识别

唐贤伦; 李伟; 马伟昌; 孔德松; 马艺玮

doi:10.11999/JEIT190124

基于条件经验模式分解和串并行CNN的脑电信号识别

doi: 10.11999/JEIT190124

重庆邮电大学自动化学院重庆 400065

基金项目: 国家自然科学基金(61673079, 61703068)，重庆市基础研究与前沿探索项目(cstc2018jcyjAX0160)

详细信息

作者简介:
唐贤伦：男，1977年生，教授，博士，研究方向为智能系统与机器人，模式识别理论与应用

李伟：男，1996年生，硕士生，研究方向为深度学习、脑电信号识别

马伟昌：男，1995年生，硕士生，研究方向为机器人控制

孔德松：男，1993年生，硕士生，研究方向为深度学习

马艺玮：女，1980年生，副教授，博士，研究方向为智能信息处理

通讯作者:
李伟　cqyddxliwei@foxmail.com

中图分类号: TP391.4
计量
- 文章访问数: 2586
- HTML全文浏览量: 1429
- PDF下载量: 119
- 被引次数: 14
出版历程
- 收稿日期: 2019-03-01
- 修回日期: 2019-11-22
- 网络出版日期: 2019-12-14
- 刊出日期: 2020-06-04

Conditional Empirical Mode Decomposition and Serial Parallel CNN for ElectroEncephaloGram Signal Recognition

School of Automation, Chongqing University of Posts and Telecommunications, Chongqing 400065, China

Funds: The National Natural Science Foundation of China (61673079, 61703068), The Basic Research and Frontier Exploration Project of Chongqing (cstc2018jcyjAX0160)

摘要

摘要:
针对运动想象脑电信号(EEG)的非线性、非平稳特点，该文提出一种结合条件经验模式分解(CEMD)和串并行卷积神经网络(SPCNN)的脑电信号识别方法。在CEMD过程中，采用各阶固有模式分量(IMF)与原始信号的相关性系数作为第1个IMF筛选条件，在此基础上，提出各阶IMF之间的相对能量占有率作为第2个IMF筛选条件。此外，为了考虑脑电信号各个通道之间的特征和突出每个通道内的特征，该文提出SPCNN网络模型对进行CEMD过程后的脑电信号进行分类。实验结果表明，在自行采集的脑电数据集上平均识别率达到94.58%。在公开数据集BCI competition IV 2b上平均识别率达到82.13%，比卷积神经网络提高了3.85%。最后，在自行设计的智能轮椅脑电控制平台上进行了轮椅前进、左转和右转在线控制实验，验证了该文算法对脑电信号识别的有效性。
- 脑电信号识别 /
- 经验模式分解 /
- 卷积神经网络 /
- 特征提取 /
- 智能轮椅
Abstract:
For the non-linear and non-stationary characteristics of motor imagery ElectroEncephaloGram (EEG) signals, an EEG signal recognition method based on Conditional Empirical Mode Decomposition (CEMD) and Serial Parallel Convolutional Neural Network (SPCNN) is proposed. In the CEMD process, the correlation coefficient between the Intrinsic Mode Functions (IMFs) and the original signal is used as the first condition to select IMFs. Based on this, the relative energy occupancy rates between the IMFs are proposed as the second condition to select IMFs. Further, to consider the characteristics between the EEG signal channels and highlight the features in each EEG signal channel, a SPCNN model is proposed to classify the processed EEG signals. The experimental results show that the average recognition rate reaches 94.58% on the dataset collected by ourselves. And the average recognition rate reaches 82.13% on the BCI competition IV 2b dataset, which is 3.85% higher than the average recognition rate of convolutional neural network. Finally, the online control experiments are carried out on the designed intelligent wheelchair platform, which proves the effectiveness of the proposed algorithm for EEG signals recognition.
- Electro Encephalo Gram(EEG) recognition /
- Empirical Mode Decomposition (EMD) /
- Convolutional Neural Network (CNN) /
- Feature extraction /
- Intelligent wheelchair

HTML全文

1. 引言

在外辐射源雷达低空小目标探测中，其背景杂波受传播环境影响而复杂多变，尤其是功率突变的杂波会在距离多普勒谱上形成点状、块状的杂波边缘，极大地影响了低空小目标的恒虚警检测性能。

当前常用的恒虚警算法根据杂波在空间与时序上的分布特性与均匀性可划分为：空域类(也称为滑窗类)和时序类。

空域类恒虚警算法中具有代表性的是单元平均(Cell Averaging)CA-CFAR算法^[1,2]，其采用临近待检测单元的杂波均值来计算检测阈值，在均匀环境下效果较好，但是在多目标及杂波边缘环境下无法消除临近目标与能量突变干扰，使得检测性能恶化。在该算法的基础上，Trunk^[3]提出了最小选择(Smallest Of)SO-CFAR以解决多目标干扰问题，Hansen等人^[4]提出最大选择(Greatest Of)GO-CFAR以解决杂波边缘环境下的能量突变问题，但是这两种算法都在均匀环境下存在一定性能损失，且并不适用于另一种非均匀环境。Smith等人^[5]在以上算法的基础上提出了变指数(Varibility Index)VI-CFAR。该算法先判决当前杂波环境的类型，然后选择与其匹配的其他恒虚警算法^[6,7]，其本质是多种算法的组合。Rohling等人^[8-11]提出有序统计(Order Statistics)OS-CFAR，虽然该算法在多目标、杂波边缘环境下性能得到部分改善，但在均匀环境下性能也有一定损失。

时序类恒虚警算法^[12-16]中，杂波图(Clutter Map)CM-CFAR是该类算法的典型代表，CM-CFAR将杂波空间划分为杂波图单元用以存储杂波数据，每个杂波图单元依靠当前数据和历史数据进行加权迭代更新。由于这类算法在迭代初期采样的数据量少，所以检测性能较差；随着迭代次数增加，检测性能会得到较大的提升。但是这类算法也具有因迭代计算而带来的记忆效应这一不可忽视的缺点。

在外辐射源雷达低空小目标探测中，由于目标回波弱、背景杂波复杂多变，空域类算法在该类应用中性能损失严重、虚警概率恶化，传统时序类算法迭代初期性能较差并且还有记忆效应等问题，为此本文提出了动态有序矩阵(Dynamic Order Matrix)DOM-CFAR算法，其将杂波空间划分成有序矩阵，利用雷达初次扫描数据的滑窗中值对有序矩阵中的有序数组进行初始化用以弥补传统时序类算法在迭代初期的性能劣势，并通过雷达扫描数据对有序数组进行动态极值替换用以获取杂波估计中值并弥补记忆效应造成的影响，最后利用杂波估计中值来计算检测阈值。

2. 检测模型

2.1 杂波与目标

在判决雷达回波 $x(t)$ 中是否存在目标时，根据二元假设检测可知

$\left. \begin{aligned} & x(t) = n(t),\qquad\;\;\;{{\rm{H}}_0} \\ & x(t) = n(t) + s(t),{{\rm{H}}_1} \end{aligned} \right\}$

(1)

其中， $s(t)$ 为目标信号， $n(t)$ 为杂波和噪声信号， ${{\rm{H}}_{\rm{0}}}$ 表示只有杂波， ${{\rm{H}}_{\rm{1}}}$ 表示存在目标。

杂波分布模型按照其概率密度函数参数个数可以分为单参和多参分布模型^[17-20]。例如瑞利分布、指数分布等只由1个参数决定其幅度分布规律的为单参分布模型；韦布尔分布、对数正态分布、莱斯分布等由多个参数决定其幅度分布规律的为多参分布模型。本文以杂波数据经平方率检波后服从指数分布为背景，其杂波概率密度函数(Probability Density Function, PDF)为

$f(x) = \frac{1}{\lambda }\exp \left( - \frac{x}{\lambda }\right),0 \le x$

(2)

其中， $\lambda$ 为每个杂波单元总功率。在 ${{\rm{H}}_{\rm{0}}}$ 假设下， $\lambda$ 为杂波和噪声的功率 $\mu$ ；在 ${{\rm{H}}_{\rm{1}}}$ 假设下， $\lambda$ 为杂波、噪声和目标的总功率 $(1 + s)\mu$ , $s$ 表示目标与杂波噪声功率 $\mu$ 的比值。由式(2)可得检测概率(Probability of Detection, PD)与虚警概率(Probability of False Alarm, PFA)为

$\begin{split} {\rm{PD}} =\,& {\rm{P}} [x > {x_{\rm{t}}}|{{\rm{H}}_1}] = \int_{{x_{\rm{t}}}}^\infty {f(x){\rm{d}} x} |{{\rm{H}}_1}\ \ \\ =\,& \exp \left( { - \frac{{{x_{\rm{t}}}}}{{(1 + s)\mu }}} \right) \end{split}$

(3)

$\begin{split} {\rm{PFA}} =\,& {\rm{P}} [x > {x_{\rm{t}}}|{{\rm{H}}_{\rm{0}}}] = \int_{{x_{\rm{t}}}}^\infty {f(x){\rm{d}} x} |{{\rm{H}}_{\rm{0}}} \\ =\,& \exp \left( { - \frac{{{x_{\rm{t}}}}}{\mu }} \right) \end{split}$

(4)

在已知参数 $\mu$ 和预设虚警概率 ${\rm{PFA}}$ 时，可由式(4)求得检测阈值T为

$T = - \mu \ln ({\rm{PFA}})$

(5)

2.2 有序矩阵

对于雷达回波的RD谱数据 ${{X}}$ ，每帧数据既有距离信息也有多普勒信息，其表达式为

${{X}} = \left[ {\begin{array}{*{20}{c}} {{x_{1,1}}}&{{x_{1,2}}}& ··· &{{x_{1,d}}}& ··· &{{x_{1,D}}}\\ {{x_{2,1}}}&{{x_{2,2}}}& ··· &{{x_{2,d}}}& ··· &{{x_{2,D}}}\\ \vdots & \vdots &{\ddots}& \vdots &{\ddots}& \vdots \\ {{x_{r,1}}}&{{x_{r,2}}}& ··· &{{x_{r,d}}}& ··· &{{x_{r,D}}}\\ \vdots & \vdots &{\ddots}& \vdots &{\ddots}& \vdots \\ {{x_{R,1}}}&{{x_{R,2}}}& ··· &{{x_{R,d}}}& ··· &{{x_{R,D}}} \end{array}} \right]$

(6)

其中， ${x_{r,d}}$ 表示RD谱 ${{X}}$ 中待检测的第 $(r,d)$ 单元数据，其中 $r$ 为第 $r$ 个距离元、 $d$ 为第 $d$ 个多普勒元。由于本文所提算法利用待检测杂波的历史数据来估计中值 ${y_{{\rm{mid}}}}$ 求解检测阈值T，为了准确地估算估计中值 ${y_{{\rm{mid}}}}$ ，这里利用有序矩阵 ${{M}}$ 来存储相关历史数据，有序矩阵 ${{M}}$ 的表达式为

${{M}} = \left[ {\begin{array}{*{20}{c}} {{{{Y}}_{1,1}}}&{{{{Y}}_{1,2}}}& ··· &{{{{Y}}_{1,d}}}& ··· &{{{{Y}}_{1,D}}}\\ {{{{Y}}_{2,1}}}&{{{{Y}}_{2,2}}}& ··· &{{{{Y}}_{2,d}}}& ··· &{{{{Y}}_{2,D}}}\\ \vdots & \vdots &{\ddots}& \vdots &{\ddots}& \vdots \\ {{{{Y}}_{r,1}}}&{{{{Y}}_{r,2}}}& ··· &{{{{Y}}_{r,d}}}& ··· &{{{{Y}}_{r,D}}}\\ \vdots & \vdots &{\ddots}& \vdots &{\ddots}& \vdots \\ {{{{Y}}_{R,1}}}&{{{{Y}}_{R,2}}}& ··· &{{{{Y}}_{R,d}}}& ··· &{{{{Y}}_{R,D}}} \end{array}} \right]$

(7)

其中， ${{{Y}}_{r,d}}$ 为有序矩阵 ${{M}}$ 中第 $(r,d)$ 单元的有序数组，其长度为 $2K + 1$ ，用于存储 ${x_{r,d}}$ 历史数据，表达式为

${{{Y}}_{r,d}} = \left[ {\begin{array}{*{20}{c}} {{y_1}}\;\;\;{{y_2}}\;\;\; ···\;\;\; {{y_k}}\;\;\; ···\;\;\; {{y_{2K + 1}}} \end{array} } \right]$

(8)

其中， ${y_k}$ 为有序数组 ${{{Y}}_{r,d}}$ 中的第 $k$ 大的数据， $K$ 为本算法的半参考窗值，K值越大，所需的存储空间越多。待检测单元 ${x_{r,d}}$ 的杂波估计中值 ${y_{{\rm{mid}}}}$ 为 ${{{Y}}_{r,d}}$ 中的 ${y_{K + 1}}$ 。

3. DOM-CFAR原理

DOM-CFAR算法实现流程如图1所示，分为初始化、迭代优化与检测3个步骤。

图 1 DOM-CFAR算法实现流程图

下载: 全尺寸图片幻灯片

3.1 初始化

初始化时，以RD谱数据 ${{X}}$ 中待检测单元 ${x_{r,d}}$ 为中心的长度为 $2K + 1$ 的距离多普勒域滑窗进行排序整理，用滑窗中值 $x{'_{r,d}}$ 对有序数组 ${{{Y}}_{r,d}}$ 中的所有数据赋初值，如式(9)所示

${y_k} = x{'_{r,d}},k \in [1,2K + 1]$

(9)

初始化后，有序数组 ${{{Y}}_{r,d}}$ 中的数据 ${y_k}$ 服从 $N({u_0} + \tau ,{\sigma ^2})$ 的高斯分布，其均值与方差为

$\qquad {u_0} + \tau = \frac{{\displaystyle\sum\limits_{k = 1}^{k = 2K + 1} {{y_k}} }}{{2K + 1}}$

(10)

$\qquad {\sigma ^2} = \frac{{\displaystyle\sum\limits_{k = 1}^{k = 2K + 1} {{{[{u_0} + \tau - {y_k}]}^2}} }}{{2K + 1}}$

(11)

其中， ${u_0}$ 为杂波理论中值， $\tau$ 为 ${y_k}$ 与 ${u_0}$ 的差值， ${\sigma ^2}$ 表示 ${y_k}$ 的方差。

由式(9)可知，初始化后有序数组 ${{{Y}}_{r,d}}$ 中 ${y_k}$ 的数值都为 $x{'_{r,d}}$ ，方差值为0，并且当K值越大，采样滑窗所包含的杂波样本越多， $x{'_{r,d}}$ 越趋近 ${u_0}$ ， $\left| \tau \right|$ 数值越小。

由于本算法初始化参考OS-CFAR的采样策略，所以在迭代初期可获得与其接近的算法性能。

3.2 迭代优化

算法迭代通过将有序数组 ${{{Y}}_{r,d}}$ 中偏离理论中值 ${u_0}$ 的数据替换，保留接近 ${u_0}$ 的数据，从而获取更加准确的估计中值 ${y_{{\rm{mid}}}}$ 。

在数值替换中对极小值 ${y_1}$ 和极大值 ${y_{2K + 1}}$ 进行动态轮询替换，以雷达扫描次数t的奇偶性为标志，当雷达扫描次数t为奇数时，将待检测单元 ${x_{r,d}}$ 赋值给 ${{{Y}}_{r,d}}$ 中的 ${y_1}$ ；当t为偶数时，将 ${x_{r,d}}$ 赋值给 ${y_{2K + 1}}$ 。重新对 ${{{Y}}_{r,d}}$ 进行排序以保持有序性，以被替换极值点为起点，未被替换极值点为终点进行单向冒泡排序。K值越大，排序耗时越多。

由于初始化时 ${{{Y}}_{r,d}}$ 中的数据都相同，所以在算法迭代中需要经过2K+1轮迭代才能将 ${{{Y}}_{r,d}}$ 中的初始数据全部替换，之后 ${{{Y}}_{r,d}}$ 中的 ${y_{{\rm{mid}}}}$ 开始变化，初始数据被全部替换也标志着算法渡过2K+1轮迭代的初期阶段。

3.3 检测阈值计算

通过杂波估计中值 ${y_{{\rm{mid}}}}$ 逆推求解杂波模型参数，即令 $\displaystyle\int_0^{{y_{{\rm{mid}}}}} {f(x){\rm{d}} x} {\rm{ = }}\frac{1}{2}$ 可得杂波参数 $\mu = \dfrac{{{y_{{\rm{mid}}}}}}{{\ln 2}}$ 。

由于预设虚警概率已知，可通过式(4)求解标称因子 $\alpha$ 和检测阈值T为

$\qquad \alpha = - \frac{{\ln ({\rm{PFA}})}}{{\ln 2}}$

(12)

$\qquad T = \alpha {y_{{\rm{mid}}}}{\rm{ = }} - \frac{{{y_{{\rm{mid}}}}\ln ({\rm{PFA}})}}{{\ln 2}}$

(13)

使用本算法进行检测处理时，若杂波数据服从单参分布模型，可直接按原流程处理；若杂波数据服从多参分布模型，则需要先保留其中1个参数为未知状态，将其它参数由参数估计方法(最大似然估计、矩估计等)估算，之后按照原流程估算这个未知参数，计算检测阈值进行检测处理。

3.4 半参考窗值K的取值

半参考窗数值K需要综合考虑内存需求、运算速度、迭代速度、期望检测性能这4个因素进行设置。为便于分析K值对迭代速度和期望检测性能的影响，这里设计如下2个评价指标

${m_{t,K}} = \frac{1}{N}\sum\limits_{i = 1}^N {|{y_{t,K,i}} - {{\bar y}_{{\rm{mid}}}}|}$

(14)

$\sigma _{t,K}^2 = \frac{1}{N}\sum\limits_{i = 1}^N {{{(|{y_{t,K,i}} - {{\bar y}_{{\rm{mid}}}}| - {m_{t,K}})}^2}}$

(15)

其中， ${y_{t,K,i}}$ 为杂波估计中值， ${\bar y_{{\rm{mid}}}}$ 为杂波理论中值， ${m_{t,K}}$ 为 ${y_{t,K,i}}$ 与 ${\bar y_{{\rm{mid}}}}$ 之差绝对值的均值， $\sigma _{t,K}^2$ 为 ${y_{t,K,i}}$ 与 ${\bar y_{{\rm{mid}}}}$ 之差绝对值的方差， $t$ 为迭代次数， $K$ 为半参考窗值， $i$ 为蒙特卡罗实验次数， $N$ 为蒙特卡罗实验总次数。均值 ${m_{t,K}}$ 越小，表明 ${y_{t,K,i}}$ 越接近 ${\bar y_{{\rm{mid}}}}$ 、准确性越高；方差 $\sigma _{t,K}^2$ 越小，表明稳定性越高；准确性与稳定性越高计算出的检测阈值更接近理论值，可获得更优的期望检测性能。

利用杂波功率 $\lambda$ 为0 dBm的指数分布杂波数据进行2×10⁴组蒙特卡罗实验，每组实验进行2000次迭代、K取值范围为[1, 100]，对获得的 ${m_{t,K}}$ 和 $\sigma _{t,K}^2$ 取自然对数绘制等高线图谱如图2、图3所示。

图 2 均值等高线

下载: 全尺寸图片幻灯片

图 3 方差等高线

下载: 全尺寸图片幻灯片

由图2、图3可知，两图左上角出现三角形异常区域，这是由算法迭代初期杂波估计中值不变暂稳态导致。K取值在[10,20]范围内时，两图的评价指标处于等高线山谷区域，均值和方差较低(异常区域数据除外)，准确性和稳定性较高，所以应尽量将K取值在[10,20]范围内。

4. 仿真分析

将DOM-CFAR与最优检测器(opt), CA-CFAR, SO-CFAR, GO-CFAR, OS-CFAR和CM-CFAR进行4×10³组蒙特卡罗仿真实验，其中背景杂波功率 $\lambda$ 为20 dBm, DOM-CFAR的半参考窗值K为15, CA-CFAR, SO-CFAR, GO-CFAR和OS-CFAR的保护单元为4、滑窗长度为30, CM-CFAR的遗忘因子为1/128，预设虚警概率PFA为10^–4。

4.1 均匀环境下的性能分析

对均匀环境下单目标分别进行10次、200次、1000次迭代优化，其中信杂比取值范围为[0 dB, 30 dB]，其检测结果如图4、图5、图6所示。

图 4 均匀杂波下迭代10次

下载: 全尺寸图片幻灯片

图 5 均匀杂波下迭代200次

下载: 全尺寸图片幻灯片

图 6 均匀杂波下迭代1000次

下载: 全尺寸图片幻灯片

由图4可知，DOM迭代初期检测性能接近OS，优于SO，较CA, GO略差。

由图5可知，经过200次迭代DOM性能已优于参照对比算法。

由图6可知，经过1000次迭代DOM性能趋近于最优检测器性能。

4.2 多目标环境下的性能分析

在目标参考前窗设置相同信杂比的第2个目标形成拥有多目标的背景杂波环境进行10次、200次迭代优化，其中信杂比取值范围为[0 dB, 30 dB]，检测结果如图7、图8所示。

图 7 多目标下迭代10次

下载: 全尺寸图片幻灯片

图 8 多目标下迭代200次

下载: 全尺寸图片幻灯片

由图7可知，DOM迭代初期拥有与OS相近的多目标检测优势。

由图8可知，DOM检测性能不受邻近目标的影响，且经过200次迭代后，其检测性能已经优于参照对比算法。

4.3 杂波边缘环境下的性能分析

设置前500个距离单元杂波功率为20 dBm，后500个距离单元杂波功率为30 dBm，形成杂波边缘环境，分别进行10次、200次迭代检测，其结果如图9、图10。

图 9 杂波边缘下迭代10次

下载: 全尺寸图片幻灯片

图 10 杂波边缘下迭代200次

下载: 全尺寸图片幻灯片

由图9可知，DOM迭代初期的虚警概率在杂波边缘的低功率区相较于CA, GO和CM更接近预设值，在杂波边缘的高功率区相较于CA, SO更接近预设值。

由图10可知，DOM经过200次迭代，虚警概率在杂波边缘的低功率区更接近预设值，在杂波边缘的高功率区相较于CA, SO和OS更接近预设值。

4.4 记忆效应分析

DOM与CM都是利用杂波历史数据计算检测阈值的时序类恒虚警算法，算法处理时每帧数据(无论是否有目标能量信息)都参与迭代并影响后续检测，从而导致算法的记忆效应。

由于DOM的参考数据为估计中值 ${x_{{\rm{DOM}}}}(t)$ , CM的参考数据为估计均值 ${x_{{\rm{CM}}}}(t)$ ，为了统一评价指标进行如下处理

$\quad {s_{{\rm{DOM}}}}(t) = \frac{{|{x_{{\rm{DOM}}}}(t) - {P_{{\rm{dB}}}} \cdot {{\bar x}_{{\rm{DOM}}}}|}}{{{P_{{\rm{dB}}}} \cdot {{\bar x}_{{\rm{DOM}}}}}}$

(16)

$\quad {s_{{\rm{CM}}}}(t) = \frac{{|{x_{{\rm{CM}}}}(t) - {P_{{\rm{dB}}}} \cdot {{\bar x}_{{\rm{CM}}}}|}}{{{P_{{\rm{dB}}}} \cdot {{\bar x}_{{\rm{CM}}}}}}$

(17)

其中， ${P_{{\rm{dB}}}}$ 为杂波功率值， ${\bar x_{{\rm{DOM}}}}$ 为理论中值， ${\bar x_{{\rm{CM}}}}$ 为理论均值。评价指标 $s$ 的数值越小，参考值越接近理论值，估算的检测阈值更准确。

先在均匀杂波无目标的环境下进行2000次迭代使DOM和CM的性能稳定，从第2001次迭代开始，连续引入50次单目标回波，仿真目标能量参与迭代对性能的影响，其中目标能量信杂比分别为4 dB, 8 dB, 16 dB，结果如图11所示。

图 11 目标干扰对性能的影响

下载: 全尺寸图片幻灯片

由图11可知，当目标干扰出现并参与迭代，对DOM的影响有一定延时，当目标干扰消失时，DOM会以更快的速度重新优化，故其拥有更优的抗目标干扰能力。

4.5 算法复杂度与耗时

由于算法中每个待检测单元都需要一组长度为2K+1的有序数组用于存储其历史数据，因此空间复杂度(使用大O表示法)为O(MK)，其中M由待检测单元总个数决定。在迭代优化中动态极值替换的逻辑计算量很小几乎可以忽略，对运算速度起决定影响的是单向冒泡排序，因此其运算复杂度为O(K)。

以大O表示法为评价指标对各个算法的空间、运算复杂度以及在硬件为i7-7700HQ CPU条件下处理10帧实测RD谱数据的串行、并行计算耗时进行分析，结果如表1所示，其中N由滑窗大小决定、K由半参考窗值决定、M由待检测单元总个数决定。

表 1 算法复杂度及耗时比较

算法名称	空间复杂度	运算复杂度	串行耗时(s)	并行耗时(s)
CA	1	MN	51.337	2.776
SO	1	MN	73.588	3.260
GO	1	MN	75.176	3.289
OS	1	MNlgN	449.651	7.632
CM	M	M	0.087	0.073
DOM	MK	MK	88.490	5.560

下载: 导出CSV

| 显示表格

由表1及前文中各算法在均匀、多目标和杂波边缘环境下的性能对比可知，本算法与现有大部分算法相比，虽然存储空间需求较大、处理速度优势不突出，但在性能上得到了较大提高。同时根据已进行的雷达低空小目标探测实验可知，其算法复杂度对于目前的GPU硬件条件来说是可行的。

5. 实测数据验证

使用武汉大学电波传播实验室的外辐射源雷达系统^[21]在有风力发电机组环境下的实测数据进行验证分析，探测条件如图12所示。

图 12 雷达探测环境

下载: 全尺寸图片幻灯片

实测数据经多径杂波和多普勒扩展杂波抑制后的RD谱中，风机叶片转动的多普勒扩展杂波残余仍然较强，形成杂波边缘环境，如图13所示。

图 13 实测RD谱

下载: 全尺寸图片幻灯片

在实测RD谱上添加100个信杂比为15 dB的随机目标则形成拥有多目标的杂波边缘环境。将DOM-CFAR与CA-CFAR, SO-CFAR, GO-CFAR, OS-CFAR和CM-CFAR进行实测验证，其中DOM-CFAR的半参考窗值K为15, CA-CFAR, SO-CFAR, GO-CFAR和OS-CFAR的保护单元为4、滑窗长度为30, CM-CFAR的遗忘因子为1/128，迭代总次数为400，预设虚警概率为10^–6，结果如表2所示。

由表2数据可知，在迭代初期，DOM相较于参照对比算法拥有较优的检测性能，而虚警性能相对较差。但是DOM通过迭代优化后，检测概率会稳步提升，虚警概率会迅速下降，从而获得更优的检测性能与虚警性能。

表 2 实测的检测概率与虚警概率

CFAR	CA	SO	GO	OS	CM迭代次数				DOM迭代次数
CFAR	CA	SO	GO	OS	10	100	200	400	10	100	200	400
检测概率(%)	72.4	51.3	78.8	86.8	1.3	62.1	76.4	93.0	88.0	98.3	98.1	98.2
虚警概率(10^–5)	2.68	9.58	1.93	22.06	0.13	2.93	4.25	4.45	20.93	5.08	2.85	1.28

下载: 导出CSV

| 显示表格

6. 结束语

本文提出了一种基于动态有序矩阵的外辐射源恒虚警检测算法，利用历史杂波数据来动态替换有序矩阵中的极值，使得杂波估计中值向理论中值收敛，通过估计中值求解检测阈值。仿真与实测结果表明，该算法在迭代初期拥有较优的检测性能，通过多次迭代可以快速达到检测概率与虚警概率的最优性能，在均匀杂波、多目标和杂波边缘的复杂环境中检测性能稳定，并且相较于同属时序类算法的CM-CFAR，拥有更优的抗目标干扰能力。

图 1 串并行卷积神经网络结构图

下载: 全尺寸图片幻灯片

图 2 Emotiv脑电采集仪

下载: 全尺寸图片幻灯片

图 3 Emotiv脑电采集仪电极安放位置

下载: 全尺寸图片幻灯片

图 4 单次脑电信号采集过程

下载: 全尺寸图片幻灯片

图 5 设定不同的阈值 $\alpha$ 时识别率情况

下载: 全尺寸图片幻灯片

图 6 设定不同的阈值 $\beta$ 时识别率情况

下载: 全尺寸图片幻灯片

图 7 采用不同处理方法的识别准确率对比

下载: 全尺寸图片幻灯片

图 8 智能轮椅系统结构图

下载: 全尺寸图片幻灯片

表 1 不同算法对5受试者脑电信号的识别准确率(%)

算法	CSP	ACSP	DBN	CNN	STFT-CNN	SPCNN	本文CEMD-SPCNN
S01	65.00	77.50	87.08	86.25	88.75	90.42	93.33
S02	81.67	82.92	87.50	87.92	89.17	91.25	94.17
S03	98.33	97.08	95.83	95.83	96.67	97.08	99.16
S04	76.25	78.33	83.33	85.42	85.42	86.25	89.58
S05	95.42	96.25	93.75	91.67	92.50	94.17	96.67
均值	83.33	86.41	89.50	89.42	90.50	91.83	94.58
方差	190.01	91.88	26.51	18.61	18.18	16.62	13.02

下载: 导出CSV

表 2 不同算法对BCI competition IV 2b数据集的识别准确率(%)

算法	Chin	Gan	Coyle	CSP	ACSP	DBN	CNN	STFT-CNN	SPCNN	本文CEMD-SPCNN
B01	70.00	71.00	60.00	66.56	67.50	66.56	72.22	75.00	76.39	80.56
B02	61.00	61.00	56.00	57.81	55.31	62.50	61.03	61.76	63.24	64.71
B03	61.00	57.00	56.00	61.25	62.19	60.00	61.11	62.50	62.50	64.58
B04	98.00	97.00	89.00	94.06	94.69	96.87	98.65	98.65	99.32	99.32
B05	93.00	86.00	79.00	80.63	76.88	82.19	86.48	87.16	87.84	88.51
B06	81.00	81.00	75.00	75.00	75.94	77.50	79.17	80.56	81.25	83.33
B07	78.00	81.00	69.00	72.50	71.25	76.56	78.47	77.08	79.17	81.25
B08	93.00	92.00	93.00	89.38	89.38	88.75	86.18	86.18	86.84	90.13
B09	87.00	89.00	73.00	85.63	81.25	85.94	81.25	82.64	84.03	86.81
均值	80.22	79.44	72.22	75.86	74.93	77.43	78.28	79.06	80.06	82.13
方差	192.19	190.03	181.69	158.75	157.50	155.85	147.92	138.63	137.52	129.78

下载: 导出CSV

表 3 各类操作在线识别准确率(%)

操作	直行	左转	右转
S01	96	84	86
S02	94	90	82
S03	98	88	92

下载: 导出CSV

参考文献(21)

RAMADAN R A and VASILAKOS A V. Brain computer interface: Control signals review[J]. Neurocomputing, 2017, 223: 26–44. doi: 10.1016/j.neucom.2016.10.024

佘青山, 陈希豪, 高发荣, 等. 基于感兴趣脑区LASSO-Granger因果关系的脑电特征提取算法[J]. 电子与信息学报, 2016, 38(5): 1266–1270. doi: 10.11999/JEIT150851

SHE Qingshan, CHEN Xihao, GAO Farong, et al. Feature extraction of electroencephalography based on LASSO-Granger causality between brain region of interest[J]. Journal of Electronics &Information Technology, 2016, 38(5): 1266–1270. doi: 10.11999/JEIT150851

YU Tianyou, XIAO Jun, WANG Fangyi, et al. Enhanced motor imagery training using a hybrid BCI with feedback[J]. IEEE Transactions on Biomedical Engineering, 2015, 62(7): 1706–1717. doi: 10.1109/TBME.2015.2402283

KAUFMANN T, HERWEG A, and KÜBLER A. Toward brain-computer interface based wheelchair control utilizing tactually-evoked event-related potentials[J]. Journal of NeuroEngineering and Rehabilitation, 2014, 11: No.7. doi: 10.1186/1743-0003-11-7

WANG Haofei, DONG Xujiong, CHEN Zhaokang, et al. Hybrid gaze/EEG brain computer interface for robot arm control on a pick and place task[C]. The 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Milan, Italy, 2015: 1476–1479. doi: 10.1109/EMBC.2015.7318649.

陈强, 陈勋, 余凤琼. 基于独立向量分析的脑电信号中肌电伪迹的去除方法[J]. 电子与信息学报, 2016, 38(11): 2840–2847. doi: 10.11999/JEIT160209

CHEN Qiang, CHEN Xun, and YU Fengqiong. Removal of muscle artifact from EEG data based on independent vector analysis[J]. Journal of Electronics &Information Technology, 2016, 38(11): 2840–2847. doi: 10.11999/JEIT160209

HSU W Y. EEG-based motor imagery classification using neuro-fuzzy prediction and wavelet fractal features[J]. Journal of Neuroscience Methods, 2010, 189(2): 295–302. doi: 10.1016/j.jneumeth.2010.03.030

WU Shanglin, WU Chunwei, PAL N R, et al. Common spatial pattern and linear discriminant analysis for motor imagery classification[C]. 2013 IEEE Symposium on Computational Intelligence, Cognitive Algorithms, Mind, and Brain, Singapore, 2013: 146–151. doi: 10.1109/CCMB.2013.6609178.

ZHANG Yu, WANG Yu, ZHOU Guoxu, et al. Multi-kernel extreme learning machine for EEG classification in brain-computer interfaces[J]. Expert Systems with Applications, 2018, 96: 302–310. doi: 10.1016/j.eswa.2017.12.015

KEVRIC J and SUBASI A. Comparison of signal decomposition methods in classification of EEG signals for motor-imagery BCI system[J]. Biomedical Signal Processing and Control, 2017, 31: 398–406. doi: 10.1016/j.bspc.2016.09.007

PARK C, LOONEY D, REHMAN N U, et al. Classification of motor imagery BCI using multivariate empirical mode decomposition[J]. IEEE Transactions on Neural Systems and Rehabilitation Engineering, 2013, 21(1): 10–22. doi: 10.1109/TNSRE.2012.2229296

DOSE H, MØLLER J S, IVERSEN H K, et al. An end-to-end deep learning approach to MI-EEG signal classification for BCIs[J]. Expert Systems with Applications, 2018, 114: 532–542. doi: 10.1016/j.eswa.2018.08.031

ZHANG Jin, YAN Chungang, and GONG Xiaoliang. Deep convolutional neural network for decoding motor imagery based brain computer interface[C]. 2017 IEEE International Conference on Signal Processing, Communications and Computing, Xiamen, China, 2017: 1–5. doi: 10.1109/ICSPCC.2017.8242581.

AN Xiu, KUANG Deping, GUO Xiaojiao, et al. A deep learning method for classification of EEG data based on motor imagery[C]. The 10th International Conference on Intelligent Computing, Taiyuan, China, 2014: 203–210. doi: 10.1007/978-3-319-09330-7_25.

LI Shufang, ZHOU Weidong, YUAN Qi, et al. Feature extraction and recognition of ictal EEG using EMD and SVM[J]. Computers in Biology and Medicine, 2013, 43(7): 807–816. doi: 10.1016/j.compbiomed.2013.04.002

TARAN S, BAJAJ V, SHARMA D, et al. Features based on analytic IMF for classifying motor imagery EEG signals in BCI applications[J]. Measurement, 2018, 116: 68–76. doi: 10.1016/j.measurement.2017.10.067

MU Yashuang, LIU Xiaodong, and WANG Lidong. A Pearson’s correlation coefficient based decision tree and its parallel implementation[J]. Information Sciences, 2018, 435: 40–58. doi: 10.1016/j.ins.2017.12.059

LECUN Y, BOTTOU L, BENGIO Y, et al. Gradient-based learning applied to document recognition[J]. Proceedings of the IEEE, 1998, 86(11): 2278–2324. doi: 10.1109/5.726791

SAKHAVI S, GUAN Cuntai, and YAN Shuicheng. Learning temporal information for brain-computer interface using convolutional neural networks[J]. IEEE Transactions on Neural Networks and Learning Systems, 2018, 29(11): 5619–5629. doi: 10.1109/TNNLS.2018.2789927

SUN Shiliang and ZHOU Jin. A review of adaptive feature extraction and classification methods for EEG-based brain-computer interfaces[C]. 2014 International Joint Conference on Neural Networks, Beijing, China, 2014: 1746–1753. doi: 10.1109/IJCNN.2014.6889525.

LU Na, LI Tengfei, REN Xiaodong, et al. A deep learning scheme for motor imagery classification based on restricted Boltzmann machines[J]. IEEE Transactions on Neural Systems and Rehabilitation Engineering, 2017, 25(6): 566–576. doi: 10.1109/TNSRE.2016.2601240

施引文献

期刊类型引用(5)

1.	刘言，刘宁波，黄勇，王中训. 利用相位特征筛选参考单元的改进CFAR方法. 烟台大学学报(自然科学与工程版). 2023(03): 371-378 . 百度学术
2.	杜改丽，封治华. 一种改进现有机载雷达探测性能的多普勒处理新技术. 现代雷达. 2022(02): 75-83 . 百度学术
3.	邹俊杰，程丰，万显荣. 外源雷达空时联合恒虚警检测分析与实验. 雷达科学与技术. 2022(04): 415-420+428 . 百度学术
4.	周子铂，王彬彬，张朝伟，刘建卫，徐颖鑫，王志会. 基于迭代对消的外辐射源雷达目标检测方法. 雷达科学与技术. 2022(05): 555-564 . 百度学术
5.	杨威，崔恒荣. 77GHz近程雷达CFAR算法研究. 通信技术. 2021(09): 2125-2131 . 百度学术