基于na-memd和互信息的脑电imf分量选择方法
【技术领域】
[0001] 本发明属于脑-机接口领域，涉及一种基于NA-MEMD和互信息的脑电頂F分量选 择方法。
【背景技术】
[0002] 脑电信号（electroencephalogram,EEG)是大脑内部的神经细胞活动在大脑皮层 的综合反映，包含着与大脑状况，思维过程等方面的相关信息，因而对脑电信号进行处理并 提取有用信息已经成为脑科学研究的重要内容。
[0003] 噪声辅助的多变量经验模式分解（noise-assistedmultivariateempirical modedecomposition,NA-MEMD)方法，是一种非线性的时频域分析方法，可以自适应地将 多变量数据同时分解成对应多通道的具有不同特征尺度的多组本征模态函数（intrinsic modefunction，頂F)。但该方法存在分解结果包含大量虚假的頂F分量问题，这会对后续 的分类产生不良影响。因此，从分解得到的各阶IMF分量中选取包含有用信息的成分显得 尤为重要。但现有的选取方法依赖于先验信息，具有主观性。
[0004] 互信息方法广泛应用于量化两个随机变量间共有的线性和非线性信息以衡量变 量的相关程度，可以用于頂F的选取。但传统互信息方法用于脑电信号时，一次只能处理两 个通道的信号。当用于多通道处理时，运算复杂，不能充分利用多通道间的有效信息。
【发明内容】
[0005] 本发明的目的就是针对现有的选取有用分量的方法存在的依赖于先验信息，运算 复杂度高，有效信息利用不充分的问题，提出基于NA-MEMD和互信息的提取脑电信号有用 分量的方法。
[0006] 为了实现以上目的，本发明方法主要包括以下步骤：
[0007] 步骤（1)使用NA-MEMD同时对多通道的EEG和辅助噪声共同组成的信号进行分 解，得到数量和频率一致的IMF分量；
[0008] 步骤⑵分别计算出多通道的EEG与多通道的EEG的頂F分量、辅助噪声与辅助 噪声的MF分量、多通道的EEG的頂F分量与辅助噪声的頂F分量的相关性；
[0009] 步骤（3)根据相关性得到的敏感因子筛选出有用的頂F分量并重构，使用共空 间模式（commonspatialpattern,CSP)对重构信号进行特征提取，并使用支持向量机 (SupportVectorMachine，SVM)对特征向量进行分类。
[0010] 本发明具有如下有益效果：
[0011] 1、减少运算量，提高特征提取效率
[0012] 本发明对传统的互信息进行改进，将单通道脑电信号看成是一个向量，将多通道 的信号组合成矩阵，然后对矩阵进行二维互信息分析。可以一次性揭示多通道的相关性，降 低了计算复杂度。
[0013] 2、不依赖于先验信息，提高特征区分度
[0014] 本发明根据二维互信息计算相关性得到的敏感因子选取有用分量，且加入不相关 的噪声通道辅助选取，选取依据客观充分，且不需要预先明确脑电信号的有用频带信息。
[0015] 图1为本发明的实施流程图。
[0016] 下面结合附图详细描述本发明基于NA-MEMD和二维互信息的提取脑电信号有用 分量的方法。
[0017] 下面逐一对各步骤进行详细说明。
[0018] 步骤一：使用NA-MEMD同时对多通道的EEG和辅助噪声共同组成的信号进行分解， 得到数量和频率一致的IMF分量：
[0019] (1)引入s通道不相关高斯白噪声信号n(t)，其长度与V通道的待处理信号x(t) 等长，假设都为L;
[0020] (2)将生成的s通道噪声信号n(t)加入到V通道待处理信号x(t)中，构成q通道 的多元信号z(t)，其中q=s+v;
[0021] (3)在一个球面上，选择q_l个采样点集，即得到q维空间的方向向量
[0022] (4)计算信号z(t)沿着方向向量，的投影，记作
[0023] (5)找出投影信号{/1 的局部极值（包括极大值和极小值）所对应的瞬时时 间fj表示极值点位置，jG[1，L];
[0024] (6)对利用样条插值法获得多变量包络曲线二：
[0025](7)对于K个方向向量，包络曲线的均值m(t)为：
[0026] = KjJl/ (t) Cl)
[0027] (8)求取Ci (t):计算z(t)与m(t)之差，则Ci (t) =z(t)-m(t)，i表示第i阶IMF 分量。如果Ci (t)满足多变量IMF的停止条件，则!Ti (t) =z(t) -Ci (t);
[0028](9)对!Ti⑴重复上述⑷~⑶步进行w次，当rw(t)变成一个单调函数，则停 止筛分，否则对C1U)重复上述步骤。
[0029] 如此，便将多通道组合信号z(t)分解成对应的頂F分量，C1 (t)即为分解得到的 IMF分量，如下公式所示：
[0031] 步骤二：分别计算出多通道的EEG与多通道的EEG的頂F分量、辅助噪声与辅助噪 声的MF分量、多通道的EEG的頂F分量与辅助噪声的頂F分量的相关性：
[0032] (1)计算互信息Ami:即计算多通道的脑电信号x(t)的第i阶頂F分量与多通道的 脑电信号x(t)本身的相关度。
[0033]Ami=I(cxi(t) ；x(t)) (3)其中，CxiU) = {cx^t)，cx2i(t)，…，cxsi(t)}，i =1，2, ...，W0
[0034](2)计算互信息BMI:即计算辅助噪声信号n(t)的第i阶頂F分量与辅助噪声信号 n(t)本身的相关度。
[0035]Bmi= I(cn；(t) ；n(t)) (4)其中，Cni (t) =Icnli (t)，cn2i (t)，… ，cnti(t)}，i= 1，2,…，W0
[0036] (3)计算Emi:即计算多通道的脑电信号x(t)的第i阶頂F分量与辅助噪声信号 n(t)的对应阶IMF分量的相关度。
[0037] Emi= I (cx i(t)；cni(t)) (5)
[0038] (4)计算敏感因子R:
[0039] R= (Ami+Bmi)/2-Emi (6)
[0040] 步骤三：根据相关性得到的敏感因子筛选出有用的MF分量并重构，再对重构信 号进行特征提取和分类。
[0041] 下面以两类脑电信号分类为例进行说明。
[0042]a、重构分量并使用CSP提取特征向量
[0043] 比较各阶頂F分量对应的敏感因子R，其值越大，说明该阶的頂F分量越敏感，特征 区分度越大，包含的有用信息越多。根据R的大小选取包含有用信息的IMF分量。为了避 免直接使用选取的IMF分量进行计算导致高的特征维数，则先重构所选分量，然后对重构 后的信号进行特征提取，具体求解算法如下：
[0044] (1)将选取的包含有用信息的IMF分量累加得到重构后的多通道脑电信号为矩阵 X，X为vXt维，V表示脑电采集时的通道数，t表示脑电采集时每次任务各通道的采样点数。
[0045] (2)将脑电信号X归一化处理后可得到协方差矩阵：
[0047] 其中不类别d样本的脑电信号，dG{1，2}，XdT表不X满转置，trace(…）表 示矩阵的迹，即矩阵对角元素之和。
[0048] (3)合成的空间协方差矩阵分解：
[0049] COVd= COV !+COV2= U〇5：U01 (8)
[0050] 其中，S为特征值对角矩阵，U。为其对应的特征向量矩阵，在这一变换中，将2特 征值按降序排列，其相对应的U。也重新排列。
[0051] (4)求白化矩阵，定义白化矩阵p如下：
[0053](5)白化变换协方差矩阵COVl和C0V2 :
[0054] WM1 = PCOV y，WM2 = PCOV 2Pt (10)
[0055](6)主成分分解：
[0056]WM1= UJWM2= U2S2U2t^以取2 :和5： 2中最大的E个特征值对应的特征 向量组成投影矩阵u，U= [U'iU'2]，其中U'i、u'2分别为最大的个特征值对应的特征 向量。
[0057] b、使用SVM对特征向量分类
[0058] 为了凸显本发明的有效性，并方便与其它方法做比较，将所提取特征向量送入SVM 中分类：
[0059] (1)在使用SVM分类前，训练样本和测试样本按140 :60的固定比例随机选择数 据，并将训练数据和测试数据都归一化到[0，1]区间。
[0060] (2)SVM的核函数选择分类效果较稳定的RBF核函数。
[0061] (3)在训练时使用5倍交叉验证，并使用网格搜索（Girdsearch)法寻找最佳惩罚 系数C和RBF参数gamma。
[0062] (4)通过训练建立支持向量机模型，并使用该模型对测试集进行预测。
[0063] (5)重复上次分类过程100次，得出平均分类准确度。
【主权项】
1. 基于NA-MEMD和互信息的脑电頂F分量选择方法，其特征在于，该方法具体包括以下 步骤： 步骤一：使用NA-MEMD同时对多通道的EEG和辅助噪声共同组成的信号进行分解，得到 数量和频率一致的IMF分量； 步骤二：分别计算出多通道的EEG与多通道的EEG的頂F分量、辅助噪声与辅助噪声的 頂F分量、多通道的EEG的頂F分量与辅助噪声的頂F分量的相关性： ⑴计算互信息Ami:即计算多通道的脑电信号x(t)的第i阶頂F分量与多通道的脑电 信号x(t)本身的相关度； Ami= I (cx ； (t) ；x(t)) (I) 其中，Cxi (t) = Icxli (t)，cx2i (t)，…，cxsi (t)}，i = 1，2, · · ·，w ;其中 w 为 IMF 分量的 分解总阶数； (2) 计算互信息Bmi:即计算辅助噪声信号n (t)的第i阶頂F分量与辅助噪声信号n (t) 本身的相关度； Bmi= I (cn ； (t) ；n(t)) (2) 其中，Cni (t) = Icnli (t)，cn2i (t)，…，cnti (t)}，i = 1，2, · · ·，w ; (3) 计算Emi:即计算多通道的脑电信号x(t)的第i阶頂F分量与辅助噪声信号η⑴ 的对应阶IMF分量的相关度； Emi= I (cx ； (t) !Cni (t)) (3) (4) 计算敏感因子R : R = (Ami+Bmi) /2-EMi (4) 步骤三：根据相关性得到的敏感因子筛选出有用的MF分量并重构，再对重构信号进 行特征提取和分类。2. 根据权利要求1所述的基于NA-MEMD和互信息的脑电頂F分量选择方法，其特征在 于，步骤一具体包括以下步骤： 步骤（1):引入s通道不相关高斯白噪声信号n (t)，其长度与V通道的待处理信号X (t) 等长，假设都为L ; 步骤（2):将生成的s通道噪声信号n (t)加入到V通道待处理信号x(t)中，构成q通 道的多元信号z (t)，其中q = s+v ; 步骤（3):在一个球面上，选择q-1个采样点集，即得到q维空间的方向向量; 步骤⑷：计算信号z(t)沿着方向向量#的投影，记作 步骤（5):找出投影信号的局部极值（包括极大值和极小值）所对应的瞬时 时间$，j表示极值点位置，j e [1，L]; 步骤（6):对利用样条插值法获得多变量包络曲线; 步骤（7):对于K个方向向量，包络曲线的均值m(t)为：(5) 步骤（8):求取cjt):计算z(t)与m(t)之差，则CiU) = z(t)-m(t)，i表示第i阶 頂F分量；如果C1⑴满足多变量頂F的停止条件，则Γι (t) = z (t) -C1⑴； 步骤（9):对!Ti (t)重复上述步骤（4)~（8)步进行W次，irw(t)变成一个单调函数， 则停止筛分，否则对C1 (t)重复上述步骤； 如此，便将多通道组合信号z (t)分解成对应的頂F分量，C1 (t)即为分解得到的頂F分 量，如下公式所示：
【专利摘要】本发明涉及一种基于NA-MEMD和互信息的选择IMF分量的方法。在脑电信号处理中，NA-MEMD方法是一种有效的自适应性的非线性时频域分解方法，但该方法存在分解结果包含大量虚假的IMF分量问题，而现有的方法选取有用IMF分量多依赖于先验信息。利用二维互信息计算出各通道信号及其IMF分量、噪声及其IMF，分量、信号IMF分量与噪声IMF分量的相关性。然后，根据相关性得到的敏感因子筛选出有用的IMF分量并重构，使用CSP对重构信号进行特征提取。本发明方法降低了直接选用IMF分量作为特征值导致的高的特征维数，使用二维互信息减少了算法的处理时间，加强了特征区分度，是一种有效的选取有用信息的方法。
【IPC分类】G06K9/00
【公开号】CN105069440
【申请号】CN201510511801
【发明人】佘青山, 韩笑, 陈希豪, 田卓
【申请人】杭州电子科技大学
【公开日】2015年11月18日
【申请日】2015年8月19日