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

权利要求书

1.一种脑电信号中眼电伪迹的去除方法，将回归方法同主成分分析方法相结合，包括以下步骤：

步骤1：对头表采集得到的数据进行主成分分解后，根据所得到的各个成分的时域波形及空间分布确定最可能为眼电伪迹的成分，其时域波形记为P，空间分布记为F；

步骤2：对步骤1确定的时域波形P按照眼电信号的特点进行1～8Hz的带通滤波，得到带通滤波后的眼电伪迹时域波形PF；

步骤3：利用各个头表电极采集的信号y_i与步骤2得到的带通滤波后的眼电伪迹时域波形PF之间的相关系数r_i计算得到各个头表电极的衰减因子B_i；

所述相关系数ri=Σj=1n(PFj-PF‾)(yij-yi‾)Σj=1n(PFj-PF‾)2·Σj=1n(yij-yi‾)2;]]>

所述衰减因子Bi=ri×sdyisdp,]]>其中：sd_yi是y_i的标准差，且sdyi=Σj=1n(yij-yi‾)2n;]]>sd_p是PF的标准差，且sdp=Σj=1n(PFj-PF‾)2n;]]>

所述n为头表电极的采样点个数；

步骤4：将步骤3得到的各个头表电极的衰减因子B_i，取其绝对值，然后求对应头表电极的衰减因子B_i的绝对值与步骤1所得的眼电伪迹的空间分布F的乘积，记为FA；

步骤5：利用FA及所有头表电极采集的脑电信号U_EEG，将FA的广义逆(FA)⁺作为空间滤波器，对所有头表电极采集的脑电信号U_EEG进行滤波，得到校正后的眼电伪迹成分的时间分布PA＝(FA)⁺×U_EEG；

步骤6：利用校正后的眼电伪迹成分的时间分布PA求该成分在各个头表电极处的产生的时域波形U_EOG＝FA×PA，最后从所有头表电极采集的脑电信号U_EEG中减去U_EOG即可得到去除眼电伪迹后的脑电信号U_TEEG。