[发明专利]一种基于对称导极的脑电信号分析方法

说明书

技术领域

本发明属于生物医学工程领域，涉及脑电信号的分析方法，具体为一种基于对称导极的脑电信号分析方法。

背景技术

脑电是神经系统的电生理活动在大脑皮质或头皮表现出的电现象，可以反映神经系统的电生理活动。如今，脑电信号的采集已十分方便，且对人体不造成任何损害，因此，可以在较短时间（如10min）内获取大量有用的脑电信号。对这些信号进行特征提取和分析，对后续研究脑电信号与大脑的生理活动之间的关联性起着关键性的作用。

发明内容

在脑电信号的分析中，因为脑区的不同，相应的活跃程度也不同，获得的脑电信号其特征值必然存在极大的差异，所以不同脑区导极之间的特征参数不存在可比性，因此提出对称导极的分析方法来解决以上问题。

本发明方法具体为：

步骤1.采集脑电信号，采用头表式贴片电极，采集被试者16导联的信号。

步骤2.去除采集信号时的外部干扰。

步骤3.利用EEGLAB中的 IIR滤波器去除50Hz的工频干扰。

步骤4.提取脑电信号的特征值，例如：频谱、功率谱、近似熵等。

步骤5.利用对称导极分析方法，将被试者左侧对应导极所得脑电信号的特征值与对应右侧导极所得信号的特征值做比值。

步骤6.比较所得比值的差异性，采用T检验、单因素方差分析或者多重比较的方法来得到不同位置导极之间特征值比值的差异性。

本发明的有益效果：利用脑结构的对称性，将大脑左右侧对应位置导极上获取的特征参数做比值，然后以该比值作为新的特征参数，从而来比较不同位置上该参数是否存在差异性。因为脑结构的对称性，那么大脑左右侧对应位置上导极获取的脑电信号具有极大的相似性。以比值来量化脑电信号的特征，由此获得的数据才存在可比性，对应的结果才具有良好的可靠性。

附图说明

图1为脑电信号采集导联图。

具体实施方式

下面结合附图对本发明的实施例作详细说明。

本发明方法具体包括以下步骤：

步骤1.采集脑电信号，采用头表式贴片电极，采集被试者16导联的信号。

步骤2.去除采集信号时的外部干扰，如采集过程中人体四肢动作造成的肌电干扰、眨眼造成的眼电干扰。

步骤3.利用EEGLAB中的 IIR滤波器去除50Hz的工频干扰。

步骤4.提取脑电信号的特征值，例如：频谱、功率谱、近似熵等。

步骤5.利用对称导极分析方法，将被试者左侧对应导极所得脑电信号的特征值与对应右侧导极所得信号的特征值做比值；如                                                、（如图1所示）。 

步骤6.比较所得比值的差异性，采用T检验、单因素方差分析或者多重比较的方法来得到不同位置导极之间特征值比值的差异性。

以下给出具体实施例：

在上述理论的基础上，我们做了以下的实验和数据分析来验证本方法的可靠性。在本实验中，我们采用慢波系数和近似熵这两个特征参数来量化脑电信号，最后的差异性检验采用单因素方差分析的方法（采用T检验、多重比较的方法同样适用于本实施例），具体步骤如下：

1.脑电信号的采集

采用上海诺承公司提供的脑电采集仪，按照国际10/20标准系统安放头皮电极，电极阻抗，采样频率为256Hz，采集导联为16导，以双耳垂、为参考电极，前额中部电极为参考接地电极，以导联，，，，，，，获取脑电信号，如图1所示。

2.数据处理

(1)预处理：去除采集信号时的外部干扰，如采集过程中人体四肢动作造成的肌电干扰、眨眼造成的眼电干扰。

(2)利用EEGLAB脑电信号处理工具箱中的IIR滤波器去除50Hz的工频干扰。

(3)截取被试者安静状态下的平稳数据进行快速傅里叶变换(FFT)。

(4)特征提取：

①慢波系数

慢波系数来源于脑电信号的频谱特征。脑电信号的功率谱按频率分为六个频带，分别是，，，，，。因此，在计算出各个频带绝对功率值的基础上，我们定义慢波系数：

                       (1)

即以低频信号功率谱与高频信号功率谱的比值来定义慢波系数。由定义可以看出慢波系数能够真实反映脑电信号中低频部分所携带的能量信息。

②近似熵