[发明专利]一种基于单通道的脑电信号中肌电噪声消除方法

权利要求书

1.一种基于单通道的脑电信号中肌电噪声消除方法，其特征是按如下步骤进行：

步骤1：采用N通道脑电测量设备记录N个通道的脑电信号x_i(t)，1≤i≤N，构成脑电信号矩阵x(t)＝[x₁(t),x₂(t),.…,x_N(t)]^T，i赋值为1；

步骤2：利用总体平均经验模态分解将第i个通道的脑电信号x_i(t)分解成为P个本征模式分量l₁(t)、l₂(t)、…、l_P(t)，构成本征模式分量矩阵l(t)＝[l₁(t),l₂(t),.…,l_P(t)]^T，每个本征模式分量满足两个条件：

条件一：极值点的数量和过零点的数量相等或者相差一个；

条件二：由局部最大值和局部最小值分别定义的包络在任一时间点t上的均值为零；

步骤3：用多集典型相关分析对本征模式分量矩阵l(t)进行盲信号分离l(t)＝Ay(t)，得到混合矩阵A和P个典型变量y_j(t)，1≤j≤P，构成典型变量矩阵y(t)＝[y₁(t),y₂(t),.…,y_P(t)]^T；

步骤4：计算每个典型变量y_j(t)的自相关系数R_j，将自相关系数低于阈值θ的典型变量判定为肌电噪声；

步骤5：在典型变量矩阵中将判定为肌电噪声的典型变量置零，得到消除肌电噪声的典型变量矩阵用步骤3得到的混合矩阵A还原消除肌电噪声的本征模式分量矩阵将中的P个本征模式分量相加得到第i个通道的消除肌电噪声的脑电信号${\tilde{x}}_i\left(t\right)=\underset{o=1}{\overset{P}{\mathrm{\Σ}}}{\tilde{l}}_p\left(t\right);$

步骤6：若i<N，取i＝i+1，对下一个通道的脑电信号x_i(t)重复步骤2到步骤5，直到得到第N个通道的消除肌电噪声的脑电信号构成消除肌电噪声的脑电信号矩阵$\tilde{x}\left(t\right)={[{\tilde{x}}_1\left(t\right),{\tilde{x}}_2\left(t\right),....,{\tilde{x}}_N\left(t\right)]}^T.$