[发明专利]基于二代小波和ICA的肌电信号降噪与去混迭方法

说明书

技术领域

本发明属于多通道信号去混迭领域，涉及一种应用于肌电信号的消噪与去混迭方法。

背景技术

肌电信号包含了丰富的肌肉运动的信息，它能够反映肢体的运动模式。目前，肌电信号已经被广泛应用于临床诊断、康复工程、运动医学等领域。表面肌电信号(sEMG)是一种在皮肤表面募集到的，伴随着肢体运动导致肌肉收缩而产生的生物电信号。基于肌电信号的人机接口通过对表面肌电信号的处理和模式识别可识别出人体不同的动作模式，控制外部环境设备工作，如基于肌电信号的仿生假手的控制，基于肌电信号的跌倒辨识。

基于肌电信号的人机接口研究中，要通过肌电信号对多个动作实现模式识别，往往需要用多个电极同时进行多路肌电信号的采集。由于人体本身存在天然电流，各神经元之间又相互连接，以及采样电极相邻通道之间存在耦合串扰等因素，使得多通道肌电信号采集到的信号存在相邻肌电信号的混迭成分。为了得到纯粹的肌电信号，进行进一步研究分析，我们需要一种有效的方法对肌电信号进行去混迭处理。

罗志增等针对脑电信号中的噪声及夹杂的眼电、心电等伪迹，采用小波分解重构消噪和基于熵估计的RADICAL算法进行消除。对3路脑电观测信号进行小波消噪和白化处理后，通过雅可比旋转矩阵分别对其中的两两组合信号用RADICAL算法进行分离，得出最优分离矩阵完成盲源分离。张洪渊等利用随机变量概率密度函数非参数估计的核函数法,对混合信号的概率密度函数及其导数进行估计,可以直接应用于所有以非线性函数代替评价函数的信号源盲分离算法。以上方法都是不错的信号盲分离算法，但是基于sEMG是一种非线性非平稳信号，以及人体活动时各肌肉的协同工作的特点，使得以上方法运用于sEMG时显现出一些不足。

sEMG的有用能量主要分布在10-500Hz之间，非常微弱(其幅值仅为μV级)，极易受到噪声的干扰，在对sEMG作去混迭处理前必须滤除其中的噪声。以改善肌电信号的盲源分离效果。

传统的小波降噪方法对信号进行多尺度分解，采用一个小波基函数去逼近不同尺度上的信号，在小波基函数与逼近信号的差异处会产生细节信息，这些细节信号在阈值处理时会被当作噪声被滤除，使降噪后的信号丢失部分有用信息。这样，使得常规的小波降噪方法无法准确的刻画非平稳信号。

综上所述，已有的信号消噪与去混迭方法在运用于肌电信号时都存在不足之处。本发明提出了一种基于二代小波变换和独立分量分析（ICA）相结合的肌电信号降噪与去混迭方法。

发明内容

本发明的目的就是针对现有技术的不足，提出一种基于二代小波变换和独立分量分析（ICA）相结合的肌电信号降噪与去混迭方法。第二代小波变换继承了传统小波的多分辨率的特性；不依赖傅立叶变换；小波变换后的系数是整数；基于多项式内插的思想，所有的运算都在空间域进行，从而摆脱了对频域的依赖。利用二代小波进行信号降噪能滤除高频噪声并保留有用信息，还能较好地保留sEMG的边缘特性。

独立分量分析（ICA）是根据从传感器监测到的观测信号来估计未知的源信号和混合矩阵的。它从信号高阶统计特性出发，主要目的按是为了提取信号的独立成分，与已有的主分量分析(PCA)方法相比，它不仅可去除分量间的低阶相关性，而且还能够去除分量间的高阶相关性，使之相互独立。因此，ICA是一种有效的盲源分离方法，相对其他方法能更全面地体现数据间的本质结构。

所以，二代小波变换和独立分量分析（ICA）相结合的方法可以有效的降低肌电信号中的噪声并去除相邻通道信号的混迭成分。

为了实现以上目的，本发明方法主要包括以下步骤：

步骤一，对仪器采集到的含噪声和线性混迭的M路观测信号，X＝[x₁,x₂,…,x_M]^Τ进行二代小波分解。

对每一路信号x[n]进行二代小波分解的具体步骤如下：

①分裂：将信号序列x[n]分成两个互不相交的子集，通常按奇样本x_odd[n]、偶样本x_even[n]分成两个等长的序列。

②预测：采用一个与数据结构无关的预测算子P，保持偶样本不变，由偶样本估计奇样本x_odd[n]＝P(x_even[n])。如果信号具有局部相关结构，偶样本集合和奇样本集合一定是高度相关的，可以以一定的精度用一个子集估计另一个子集。原值与预测值之差作为小波系数d[n]，表示为d[n]＝x_odd[n]-P(x_even[n])。