[发明专利]一种基于非高斯时序模型的脑电特征提取方法

说明书

技术领域

本发明涉及脑电数据分析领域，具体涉及一种基于非高斯时序模型的脑电特征提取方法。

背景技术

头皮脑电信号数据包含了大量与大脑特征以及状态相关的信息，是进行大脑状态判断的重要工具。由于大脑的时变性和非线性性，其产生的信号也是时变的、非线性的，同时由于头皮脑电信号容易受到眼电、肌电以及心电信号的干扰，有效而鲁棒的脑电信号特征提取方法成为难题。

申请号为200910196746.3的发明公开了一种脑电波分析方法，应用了经典的时频域分析和主成份分析方法解决了脑电信号特征提取的难题，成功地提取了与人体紧张，疲劳与放松息息相关的时频域参数，并且将其映射到主成分空间中去，还运用了支持向量机在主成份空间中高效地分析非线性关系，提高了判读的精确性和有效性。但是，该发明能够处理的频率范围有限，不能满足更宽频段的需求，而且对于脑电信号中的噪声信号没有进行有效的处理。

申请号为200710163302.0的发明公开了一种脑电超慢波分析方法及其数据处理装置，它通过对0.5-50Hz范围内的脑电波进行系列分析，提取出其中超慢波，并分析超慢谱系功率，用以反映所对应的脑神经递质的功能状况，为建立一种无创伤脑神经递质功能检测方法提供了条件。但是，该发明对于脑电数据特征具有针对性，当特征值不一致时效果较差。

通过对脑电数据信号的监测和分析，可以得到大脑的状态，如深度睡眠，大脑疲劳，大脑警觉等状态，同时，也为人类与计算机的交互控制提供了新的手段。

对脑电数据信号进行特征提取，得到可识别的，可靠程度高的特征值是脑电数据分析的基础环节，具有重要价值，因此，需要提供一种脑电特征的提取方法，提取得到能够反应大脑状态变化信息的特征值。

发明内容

本发明提供了一种基于非高斯时序模型的脑电特征提取方法，结合脑电信号的时频特征、形态特征以及复杂度特征等多种特征，对大脑状态进行有效刻画，同时通过状态空间模型对大脑状态进行实时跟踪，有效去除伪迹，能够有效反应大脑状态的变化，敏感性高，抗噪性好。

一种基于非高斯时序模型的脑电特征提取方法，包括以下步骤：

（1）获取待处理脑电数据和两组训练脑电数据，去除待处理脑电数据和两组训练脑电数据中的伪迹，分别获得待处理脑电数据的有效频段和两组训练脑电数据的有效频段，再分别将待处理脑电数据的有效频段和每组训练脑电数据的有效频段分为若干数据段；每组训练脑电数据中均包括两种大脑状态下的脑电数据。

不同的大脑状态对应具有不同特征的脑电数据，例如深度睡眠时，脑电数据波动较为平缓；剧烈运动时，脑电数据波动很大。

作为优选，每组训练脑电数据中的两种大脑状态所对应的时长相同，且每种大脑状态所持续的时间不少于1分钟。

例如，需要区分的两种大脑状态分别为深度睡眠和剧烈运动时的大脑状态，深度睡眠和剧烈运动的时长分别为2分钟。

作为优选，所述步骤（1）采用带通滤波去除伪迹，去除伪迹后获得的有效频段频率为1.6～70Hz。带通滤波可以采用2阶巴特沃兹滤波器（Butterworth），有效频段的频率可以根据需要通过滤波参数进行选择。

步骤（1）中数据分段采用滑动时间窗的方法，时间窗的长度为1s，步长为0.5s。数据分段的目的在于确定数据处理的最小单元。滑动时间窗的长度以及滑动步长可以根据需要选取，时间窗的长度越短，滑动步长越短，划分的数据段越多，则提取得到的时频特征值和形态特征值越精确，但是，相应的计算量也大。

（2）提取步骤（1）中每个数据段的时频特征、形态特征和复杂度特征，得到相应的时频特征值、形态特征值和复杂度特征值，每个数据段的时频特征值、形态特征值和复杂度特征值构成一个特征向量。

步骤(2)中的时频特征提取采用经验模态分解，取前三个固有模态函数，利用下式计算得到每个数据段的时频特征值VoIMF：