[发明专利]一种自适应脑电信号中眼电伪迹的自动去除方法

权利要求书

1.一种自适应脑电信号中眼电伪迹的自动去除方法，其特征在于，包括以下步骤：

1)脑电信号的经验模态分解

脑电信号X(t)是非平稳信号，进行经验模态分解EMD，得到各IMF分量c_i(t)和余项r(t)；

其中，i为分解得到的IMF分量的个数，i的值完全由X(t)决定；

2)对脑电信号EMD分解所得的所有IMF分量c_i(t)进行希尔伯特变换，求得瞬时频率f_i(t)，基于瞬时频率，利用以下公式(Ⅰ)进行带通滤波，提取出脑电信号频带内的信号g_i(t)；

gi(t)=0iffi(t)>fh0iffi(t)<flci(t)ifohter---(I);]]>这里f_l为1Hz，f_h为30Hz；

3)利用公式阈值函数公式(Ⅱ)和阈值滤波函数公式(Ⅳ)，进行阈值滤波，去除眼电信号，得到信号w_i(t)；

τ_i＝mean|M_i-m_i|+std|M_i-m_i|.............................(Ⅱ)

其中，M_i为第i个IMF分量g_i(t)的各个极值点的值；mean表示求均值，std表示求标准差；m_i为第i个IMF分量g_i(t)的时间均值；

假设时间序列g_i(t)的长度为N，则

mi=1NΣt=1Ngi(t)...(III)]]>

阈值函数τ_i，表征了g_i(t)信号以均值m_i为基准线的振幅上限；在此基础上，针对眼电伪迹幅度大的特点，设计了阈值滤波函数公式(Ⅳ)；

wi(t)=miif|gi(t)-mi|>τigi(t)if|gi(t)-mi|≤τi...(IV)]]>

根据公式(Ⅴ)利用滤波后的全部模态分量w_i(t)进行数据重构，得到不含眼电的比较纯净的脑电信号Y(t)；

Y(t)=Σi=1nwi(t)...(V).]]>