一种基于卷积盲源分离的高密度表面肌电信号分解方法与流程

本发明涉及表面肌电信号分解领域，具体是一种基于卷积盲源分离的高密度表面肌电信号分解方法。

背景技术：

新发展起来的高密度表面肌电信号(high-densitysurfaceemg，hd-semg)技术利用粘贴在人体皮肤上的表面电极栅格采集皮下肌肉电信号，具有无创性、易于患者接受的特点。肌电信号是在脊髓中枢发出收缩指令或受到外界施加的刺激时，多个参与发放的运动单元(motorunit，mu)兴奋而产生的运动单元动作电位(motorunitactionpotential,muap)在检测电极处的时间和空间上综合叠加的结果。高密度表面肌电信号分解是指从观测到的多通道高密度肌电信号中识别出单个运动单元的运动元动作电位并将其归类，从而揭示构成肌电信号中主体运动单元动作电位序列的过程，进而获得运动单位的募集和发放信息、muap的波形特征及其发放率等信息，可应用于神经源性和肌源性疾病的辅助诊断和神经肌肉系统的正常、病理及疲劳过程的机理研究，对康复医学及运动医学也有具有重要意义。

由于人体组织具有低通滤波作用，表面肌电信号的信噪比较低，而muap波形的变异性强且叠加程度高，导致表面肌电信号分解的研究很困难。目前，hd-semg信号分解的研究方法主要有：模板匹配法、人工智能算法、盲源分离算法、独立成分分折(ica)、卷积核补偿算法等方法。目前研究方法在肌肉收缩程度较大时的分解效果不理想，主要原因是募集的运动单元较多时，各个运动单位相互之间的影响加大，检测到的肌电信号muap波形的变异和叠加程度随之增大。总之，高密度表面肌电信号分解的研究还处于起步阶段，是肌电信号处理研究领域的难点之一。

技术实现要素：

本发明目的在于提供了一种基于卷积盲源分离的高密度表面肌电信号分解方法。本发明有效提高分解的准确性，能应用于神经源性和肌源性疾病的辅助诊断和神经肌肉系统的正常、病理及疲劳过程的机理研究，对康复医学及运动医学也有具有重要意义。

为了达到上述目的，本发明的技术方案是：

一种基于卷积盲源分离的高密度表面肌电信号分解方法，其特征在于：包括如下步骤：

1)利用粘贴在人体皮肤上的表面电极栅格采集皮下肌肉电信号，即为恒力输出下的高密度表面肌电信号，纠正不同电极间的时延，并进行时间对齐；然后采用迭代紧缩方法进行卷积盲源分离过程，提取运动单元信号源；

2)重复卷积盲源分离过程，提取更多个运动单元信号源，直到重复迭代3次后仍未检测出新的运动单元信号源为止；

3)对所有提取出的运动单元信号源的发放序列归类整理，剔除重复以及不合理的发放时刻向量，再次审核运动单元信号源；

4)利用审核通过的运动单元信号源进行神经源性和肌源性疾病的辅助诊断、神经肌肉系统的正常、病理及疲劳过程的机理研究。

所述的步骤1)中对高密度表面肌电信号做时延纠正和时间对齐的方法具体是：利用非线性能量算子(nonlinearenergyoperator,neo)对测量信号滤波，然后利用巴特利特(bartlett)窗函数与滤波结果卷积以去除由于噪声产生的虚假波峰，通过计算neo滤波信号间的互相关性时延纠正，实现时间对齐。

所述的步骤1)中卷积盲源分离过程包括如下步骤：

(1)利用峭度对比函数最大化算法提取出一个运动单元源信号；

(2)对提取的运动单元源信号进行去噪处理，以保证在每个假定发放时刻的波形相似性；

(3)利用自适应阈值化算法检测波峰，以波峰为中心识别出此运动单元对应的所有活动段，提取运动单元发放序列；

(4)估算本次提取出的运动单元源在各个电极产生的信号贡献，并将其从测量信号中减除，得到更新的测量信号。

所述的步骤(1)中，首先利用公式(5)和公式(6)迭代计算以最大化公式(4)给出的峭度对比函数，提取出一个运动单元信号源；

所述的公式(4)为：

所述的公式(5)为：y^(p)(n)＝(b^(p)*x)(n)；

所述的公式(6)为:

其中x为测量所得信号；y(n)为运动单元信号源的估计向量；p为迭代系数；kurt(y)为y(n)的峭度；μ^(p)为步长，为梯度；b(k)为长度为k的有限冲击响应；j(b(k))为峭度对比函数。

所述的步骤(3)具体为：首先利用自适应阈值化算法检测波峰，其阈值选取为信号均方根的5倍，以波峰为中心识别出此运动单元对应的所有活动段，并设定时间窗宽度为10.5ms，记录所有活动段的时间、波形特征信息，利用欧氏距离计算所有活动段的相似度，并计算所有相似性较高的活动段的平均波形生成运动单元动作电位波形模板，保存此运动单元发放序列的相关数据。

所述的步骤3)中信号贡献计算方法为：利用线性最小二乘法，从所有测量信号中减除来自本次提取的运动单元源的贡献量，得到了(n-1)个源的紧缩矩阵；也就是说从降维矩阵中分离出剩余的(n-1)个信号源；具体方法为：从测量信号x^(q-1)(n)中减去按公式(8)计算的来自第(q-1)个源yq-1(n)的贡献z^(q-1)(n)，得到按公式(7)更新测量信号；

所述的公式(7)为：x^(q)(n)＝x^(q-1)(n)-z^(q-1)(n)；

所述的公式(8)为：

其中，(t^(q)(k))k∈z是一个滤波器的脉冲响应；

通过最小化公式(9)得到更新的测量信号；

所述的公式(9)为：ε(t^(q))＝e{|x^(q)(n)|²}，

ε(t^(q))为更新后的测量信号的能量；

所述的步骤3)中平均每秒多于5个波峰则认为是生理上不合理的运动单元信号源。

本发明的有益效果：本发明提出的基于卷积盲源分离的高密度表面肌电信号分解方法，将高密度表面肌电信号建模为线性卷积矩阵，并采用紧缩方法逐步迭代估计各个运动单元源，利用峭度对比函数最大化算法提取运动单元源信号，通过后续去噪、波形相似性比较以及再审核等处理获得生理上合理的运动单元发放序列，可以有效提高分解的准确性。本发明应用于神经源性和肌源性疾病的辅助诊断和神经肌肉系统的正常、病理及疲劳过程的机理研究，对康复医学及运动医学也有具有重要意义。

附图说明

图1是本发明的流程图。

具体实施方式

实施例1

以下结合附图实例对本发明作进一步详细描述。

1.高密度表面肌电信号的mimo-lti系统模型

hd-semg信号建模为多输入多输出线性时不变(multiple-input-multiple-outputlineartimeinvariant，mimo-lti)系统在n维运动单元信号源s作为输入产生的m维输出信号x。此mimo-lti系统为一个有限冲击响应滤波器，a是其长度为l的单位脉冲响应，为m×n矩阵。那么，测量所得信号x可以表示为：

为了实现盲源信号分离，通常假定：1)信号源都为零均值单位方差平稳随机信号，且相互之间是统计不相关的；2)运动单元源的数目不超过测量信号通道数，即n≤m。从测量信号中分离出源信号的问题，就转化为求解分离算子b，通过公式(2)计算可以得到其输出y(n)，y(n)即为运动单元信号源s(n)的估计向量。

2.基于紧缩迭代的卷积盲源分离算法

利用迭代紧缩方法，利用公式(3)通过估计分离算子b相应行逐步计算提取出y(n)，其中bi(k)是长度为k的有限冲击响应。

迭代紧缩方法分为以下几步：

1)最大化公式(4)所示的峭度对比函数，其中kurt(y)为y(n)的峭度。

采用梯度算法，在第p次迭代中，按公式(5)和公式(6)计算，计算运动单元信号源s(n)的估计y(n)，并更新分离算子：

y^(p)(n)＝(b^(p)*x)(n)

(5)

其中，μ^(p)为步长，为梯度。

2)通过上述步骤提取出一个运动单元信号源的估计后，利用线性最小二乘法，从所有测量信号中减除来自此源的贡献量，得到了(n-1)个源的紧缩矩阵。问题就转化为，从降维矩阵中分离出剩余的(n-1)个信号源。具体操作为：从测量信号x^(q-1)(n)中减去按公式(8)计算的来自第(q-1)个源yq-1(n)的贡献z^(q-1)(n)，得到按公式(7)更新测量信号：

x^(q)(n)＝x^(q-1)(n)-z^(q-1)(n)

(7)

其中，(t^(q)(k))k∈z是一个滤波器的脉冲响应，通过最小化公式(9)得到：

ε(t^(q))＝e{|x^(q)(n)|²}

(9)

3)重复迭代上述两步，直到分离出所有n个信号源为止。

基于上述mimo-lti系统模型和卷积盲源分离算法，本实施例的一种基于卷积盲源分离的高密度表面肌电信号分解方法，其流程图如图1所示，包括以下步骤。

1)利用粘贴在人体皮肤上的表面电极栅格采集皮下肌肉电信号，即为恒力输出下的高密度表面肌电信号，纠正不同电极间的时延，并进行时间对齐；然后采用迭代紧缩方法进行卷积盲源分离过程，提取运动单元信号源；

所述的步骤1)中对高密度表面肌电信号做时延纠正和时间对齐的方法具体是：利用非线性能量算子(nonlinearenergyoperator,neo)对测量信号滤波，然后利用巴特利特窗函数与滤波结果卷积以去除由于噪声产生的虚假波峰，通过计算neo滤波信号间的互相关性时延纠正，实现时间对齐。

所述的步骤1)中卷积盲源分离过程包括如下步骤：

(1)利用峭度对比函数最大化算法提取出一个运动单元源信号；

所述的步骤(1)中，首先利用公式(5)和公式(6)迭代计算以最大化公式(4)给出的峭度对比函数，提取出一个运动单元信号源；

所述的公式(4)为：

所述的公式(5)为：y^(p)(n)＝(b^(p)*x)(n)；

所述的公式(6)为:

其中x为测量所得信号；y(n)为运动单元信号源的估计向量；p为迭代系数；kurt(y)为y(n)的峭度；μ^(p)为步长，为梯度；b(k)为长度为k的有限冲击响应；j(b(k))为峭度对比函数。

(2)对提取的运动单元源信号进行去噪处理，以保证在每个假定发放时刻的波形相似性；

(3)利用自适应阈值化算法检测波峰，以波峰为中心识别出此运动单元对应的所有活动段，提取运动单元发放序列；

所述的步骤(3)具体为：首先利用自适应阈值化算法检测波峰，其阈值选取为信号均方根的5倍，以波峰为中心识别出此运动单元对应的所有活动段，并设定时间窗宽度为10.5ms，记录所有活动段的时间、波形特征信息，利用欧氏距离计算所有活动段的相似度，并计算所有相似性较高的活动段的平均波形生成运动单元动作电位波形模板，保存此运动单元发放序列的相关数据。所述的步骤3)中平均每秒多于5个波峰则认为是生理上不合理的运动单元信号源。

(4)估算本次提取的运动单元源在各个电极产生的信号贡献，并将其从测量信号中减除，得到更新的测量信号。

2)重复卷积盲源分离过程，提取更多个运动单元信号源，直到重复迭代3次后仍未检测出新的运动单元信号源为止；

3)对所有提取出的运动单元信号源的发放序列归类整理，剔除重复以及不合理的发放时刻向量，再次审核运动单元信号源；

所述的步骤3)中信号贡献方法为：利用线性最小二乘法，从所有测量信号中减除来自本次提取的运动单元源的贡献量，得到了(n-1)个源的紧缩矩阵；也就是说从降维矩阵中分离出剩余的(n-1)个信号源；具体方法为：从测量信号x^(q-1)(n)中减去按公式(8)计算的来自第(q-1)个源yq-1(n)的贡献z^(q-1)(n)，得到按公式(7)更新测量信号；

所述的公式(7)为：x^(q)(n)＝x^(q-1)(n)-z^(q-1)(n)；

所述的公式(8)为：

其中，(t^(q)(k))k∈z是一个滤波器的脉冲响应；

通过最小化公式(9)得到更新的测量信号；

所述的公式(9)为：ε(t^(q))＝e{|x^(q)(n)|²}；，

ε(t^(q))为更新后的测量信号的能量。

4)利用审核通过的运动单元信号源进行神经源性和肌源性疾病的辅助诊断、神经肌肉系统的正常、病理及疲劳过程的机理研究。

本实施例的一种基于卷积盲源分离的高密度表面肌电信号分解方法，将高密度表面肌电信号建模为线性卷积矩阵，并采用紧缩方法逐步迭代估计各个运动单元源，利用峭度对比函数最大化算法提取运动单元源信号，通过后续去噪、波形相似性比较以及再审核等处理获得生理上合理的运动单元发放序列，可以有效提高分解的准确性。本发明应用于神经源性和肌源性疾病的辅助诊断和神经肌肉系统的正常、病理及疲劳过程的机理研究，对康复医学及运动医学也有具有重要意义。