[发明专利]阵列式sEMG信号自动分解方法

说明书

本发明提供一种阵列式sEMG信号自动分解方法，首先对阵列式sEMG信号预处理，并且计算初始发放序列向量；其次利用时域减法确定运动单元个数，即聚类个数；然后根据聚类个数，利用最小距离分类器对阵列式sEMG波形聚类；最后重新计算新的运动单元发放序列向量，循环执行程序直至分解完成，并对所有发放序列归类整理，优化结果。该分解方法准确性高，计算快速、实现简单。

技术领域

本发明涉及一种阵列式sEMG信号自动分解方法。

背景技术

表面肌电信号( surface EMG, sEMG) 是利用表面电极从人体体表检测肌电信号，与针电极肌电信号(Needle EMG, NEMG)相比，它具有无创性、易于患者接受的特点，因此应用前景广阔。实验表明，利用阵列式sEMG能够提高运动单元(MU)的检测率，特别是提高小幅值运动单元活动电位(MUAP)的检出和识别效果。sEMG的起源是MUAP，活动电位由给定肌肉收缩过程中所激活的每一个运动单位所释放。在任何一个给定的募集模式，众多的运动单位以非同步的模式被激活，这些运动单位活动的总和构成了肌电信号的强度。阵列式sEMG信号实质上是多个运动单位活动电位的总和，其波幅典型地在1~5000uv之间，频率范围为10~400Hz。在临床上，通过阵列式sEMG可以较全面地了解神经肌肉的功能状态，鉴别神经源性和肌源性疾病，判断神经损伤的部位、程度及恢复状况，并且阵列式sEMG信号的检测分析对康复医学及运动医学也有具有重要意义。

阵列式sEMG分解实质上是对sEMG包含的运动单元发放序列进行分类，提取出运动单元各电位序列。目前，sEMG分类方法主要有：K均值聚类算法、模板匹配法、人工神经网络(ANN)算法、实时线性混叠盲信号分离算法、独立成分分折(ICA)、卷积核补偿算法等方法。K-均值聚类算法需要指定聚类的类别数，而在肌电信号缺少运动单元发放的先验知识，难以对类别进行精确的指定。模板匹配法由于模板获取困难，应用受限。ANN可以解决含更多迭加波形情形和在低信噪比时更好地消除绝对误差，然而，ANN方法一旦训练后，网络就固定不变，当模板的形状发生变化时，神经网络还需重新训练，所以强健性不好。ICA是一种盲信号分解技术, 它假设构成肌电信号的各MUAPT 相互独立，然后把信号分解成若干相互独立的成分。卷积核补偿算法方法是一种盲信号分解方法，该方法在波形叠加严重时效果仍然不理想。阵列式sEMG的信噪比较低，MUAP波形的变异性强且相互间的叠加程度较大，这是导致其分解困难的主要原因。很多研究者将插入式肌电信号的分解方法做一些修改应用到表面肌电信号的分解研究中，但是都不能达到插入式肌电信号的分解效果。总的来看，表面肌电信号分解研究还处于探索阶段，是肌电研究领域的难点之一。

发明内容

鉴于上述问题，本发明的目的在于提供一种阵列式sEMG信号自动分解方法。

为实现上述目的，针对阵列式表面电极肌电信号，提出基于最小距离器的聚类方法，对发放波形进行分类，并且在这个过程中，不求解矩阵，通过卷积核补偿方法得到肌肉运动单元的数目和发放序列，实现阵列式sEMG的分解。由于该方法同时考虑发放波形和发放时刻，相对于其它方法，该方法具有sEMG分解精度高的优点。

本发明公开了一种阵列式sEMG信号自动分解方法。其包括以下步骤：

步骤一：阵列式sEMG信号预处理：对sEMG信号滤波，剔除干扰；

步骤二：计算初始发放序列向量：利用阵列sEMG信号各通道信号的相关性，计算初始发放序列向量，作为提取初值；

步骤三：确定运动单元个数：利用时域减法，从sEMG中获得运动单元个数，作为聚类个数；

步骤四：对阵列式sEMG波形聚类：根据所述聚类个数，利用最小距离分类器，对阵列式sEMG波形聚类；

步骤五：计算新的运动单元发放序列向量：根据所述波形聚类结果，取含最多时刻的类，计算该类时刻波形均值，计算新的发放序列向量；