[发明专利]基于变分模态分解和二阶盲辨识的脑电信号伪迹去除方法

权利要求书

1.基于变分模态分解和二阶盲辨识的脑电信号伪迹去除方法，其特征在于，包括以下步骤：

S1、采用变分模态分解VMD方法对被眼电伪迹EOG和肌电伪迹EMG同时污染的脑电混合信号y(t)进行分解，得到p个固有模态分量IMF，将固有模态分量IMF表示为固有模态分量矩阵u(t)，u(t)＝[u₁(t)，u₂(t)，...，u_p(t)]；

S2、将p个IMF作为盲源分离算法二阶盲辨识SOBI的输入信号，进行进一步分离，得到p维信号分量矩阵s(t)，s(t)＝[s₁(t)，s₂(t)，...，s_p(t)]；

S3、计算经过S1和S2步骤后，得到的s(t)中的每个信号分量的模糊熵值，根据模糊熵值识别出伪迹分量，并将识别出伪迹分量置零，得到不含伪迹成分的信号矩阵

S4、对信号分量矩阵s(t)中剩余的信号分量，即没有置零的信号分量，进行逆SOBI变换，得到不含伪迹成分的固有模态分量矩阵

S5、对固有模态分量矩阵求和，重构出干净脑电信号x(t)。

2.根据权利要求1所述的单通道脑电信号伪迹去除方法，其特征在于，所述步骤S1具体包括：

S101、对被眼电伪迹(EOG)和肌电伪迹(EMG)同时污染的脑电混合信号y(t)建立变分约束模型：

其中，k为模态分解个数，K为总的模态分解个数，u_k(t)为第k个模态分量，ω_k为u_k(t)的中心频率，是每个模态分量信号相应的解析信号，是带约束的变分模型，是模型约束，表示分解的各模态分量相加得到原始信号y(t)；

S102、在约束条件下重构出{u_k}和{ω_k}，在式(1)基础上引入二次惩罚因子α和拉格朗日乘子λ，将上述约束模型转化为非约束模型，表达式为：

S103、利用交替方向乘法器算法求解上述非约束模型得到拉格朗日函数的鞍点，从而得到窄带IMF分量，

1)初始化参数

2)n＝n+1，根据下式更新u_k

3)k＝k+1，重复步骤2，直到k＝K，终止上式，根据下式更新ω_k

4)k＝k+1，重复步骤3，直到k＝K，终止上式，

5)更新拉格朗日乘数λ

式中，γ为噪声容忍度，满足信号分解的保真度要求，n为迭代次数n＝{1，...，N}。分别对应y(t)，λ(t)的傅里叶变换，

6)重复步骤(2)-(5)，直到满足下式收敛条件，循环迭代的收敛条件以各IMF更新前后之差与更新前IMF的比值作为系数：

对于给定判断精度ε＞0，结束循环，输出各IMF的中心频率ω_k。

3.根据权利要求2所述的单通道脑电信号伪迹去除方法，其特征在于，所述步骤S2具体包括：

S201、首先对上述VMD分解得到的p维信号u(t)进行预白化处理，计算白化矩阵Q，白化后的信号记为z(t)，

z(t)＝Qu(t) (7)

S202、计算z(t)的多个时延τ∈{τ₁，τ₂，......τ_k}的采样协方差矩阵

R(τ)＝E[z(t+τ)z(t)^T]＝AR_z(τ)A^T (8)

S203、对上式计算的各个协方差矩阵R(τ_j)，进行联合近似对角化，计算正交矩阵V

V^TR(τ_j)V＝D_j (9)

式中，{D_j}是一组对角矩阵，

S204、估计混合矩阵A和源信号s(t)：

A＝Q^-1V (10)

s(t)＝A^-1u(t)＝V^TQu(t)#(11) 。