基于肌肉协同的上肢多关节同步比例肌电控制方法与系统与流程

本发明属于肌电信号处理及肌电控制领域，涉及肌肉协同分析方法和基于协同激活模型的多关节同步连续运动估计方法，特别涉及基于肌肉协同的上肢多关节同步比例肌电控制方法与系统。

背景技术：

表面肌电信号(surfaceelectromyography,semg)因其抗干扰性好易获取，且包含丰富的运动控制信息，被广泛应用于康复机器人、智能假肢和人机交互等领域。semg信号中包含的丰富信息能够反映肢体的运动意图和运动状态，为智能康复辅助设备提供多种不同模式的运动控制指令，与其它依靠按钮或者人工语音操控的辅助设备相比，肌电控制有着更好的控制本能性，患者使用起来更加灵活自然，也更能从心理上接纳。

当前的肌电控制接口主要采用动作模式分类的方法，这种方法虽然能够取得较好的分类准确率，但在这种控制策略下，设备的运行速度、位移等参数都是预先设定的，并且运动控制过程中多类动作之间的转换存在“模糊期”，导致智能辅助设备运行不流畅，这大大降低了该方法的实用性。为了实现对康复辅助设备连续流畅的控制，使控制过程更加灵活和直观，使控制量有更加明确的物理意义，采用表面肌电信号估计运动过程中关节的连续运动信息，来获取连续流畅的肌电控制指令逐渐受到关注。

在连续运动估计方面，国内外学者主要采用骨骼肌肉模型法和黑箱法这两类方法，均取得了不错的效果，获得了较高的估计准确率。但是这些研究大都是针对单关节单自由度的序列运动，而人体的自然运动往往是多关节多自由度的同步运动。为了解决多自由度运动问题，近年来一种新的基于肌肉协同分析的连续运动估计方法正得到广泛的研究和应用。

因此，本发明在准确有效的肌肉协同分析的基础上，建立协同激活模型实现上肢多关节的同步连续运动估计，并结合多自由度并行比例控制策略，搭建上肢多关节同步比例肌电控制系统，实现智能设备多自由度的同步连续运动控制，使其具备完整复现人体自然运动的潜能。

技术实现要素：

本发明一方面的目的在于提供一种基于肌肉协同的上肢多关节同步比例肌电控制方法，包括：

s1，对于上肢腕关节和肘关节的单独屈伸运动和同步屈伸运动过程中的表面肌电信号进行采集，并对采集的信号进行预处理；

s2，采用一种半监督式非负矩阵分解方法进行肌肉协同分析，具体包括：

s21，对预处理之后的表面肌电信号计算积分肌电值iemg作为肌肉激活度矩阵f(e)；

s22，根据肌肉协同理论构建肌肉协同模型f(e)n×t＝wn×k×h(t)k×t；

设定单个关节屈伸运动的协同元个数为2，对单独腕关节运动时的肌肉激活度进行nmf分解，得到两个肌肉协同元w1,w2；对单独肘关节运动时的肌肉激活度进行nmf分解，得到两个肌肉协同元w3,w4；其中，wi＝[w1i,w2i,…,wni],(i＝1,2,3,4)，wni(n＝1,2,...,n)代表第n块肌肉对第i个肌肉协同元wi的贡献度；

取多组肌电信号样本，多次计算wi，并分别对每个wi的多组计算结果取均值处理，最终得到的四个wi(mean)按列组合即可获得作为监督元的总体肌肉协同元矩阵w＝[w1(mean)w2(mean)w3(mean)w4(mean)]；

s23，提取出激活系数序列h(t)＝w⁺·f(e)，式中w⁺为协同元矩阵w的m-p逆；

s3，通过支持向量回归分别构建上肢腕关节和肘关节的激活系数序列和关节角度信息的协同激活模型，利用所述协同激活模型对上肢运动过程中腕关节和肘关节的角度进行同步估计。

在上述的方法，所述协同激活模型的训练方法包括：同步采集腕关节和肘关节运动过程中的表面肌电信号和角度信号并进行预处理；提取肌肉激活度矩阵f(e)，获取协同元矩阵w，提取出激活系数序列h(t)；将激活系数序列h(t)低通滤波预处理后归一化到[0,1]区间作为协同激活模型的输入，对相应关节的实测角度值进行相同区间的归一化处理，作为目标值进行模型的训练。

在上述的方法，利用sj＝αj·θj+βj(j＝1,2)多自由度并行比例肌电控制策略将估计出的所述角度信息转化为智能设备多个自由度的连续运行位移，其中，sj是智能设备第j个自由度的目标位移，αj是相应的比例系数，θj为人体上肢运动过程中估计出的两个关节角度，其中j＝1代表腕关节，j＝2代表肘关节，βj则表示修正系数。

在上述的方法，将所述位移数据输出至智能设备。

在上述的方法，采集表面肌电信号时肌电电极放置在桡侧腕长伸肌、桡侧腕屈肌、肱二头肌和肱三头肌，获取四通道肌电信号。

在上述的方法，对采集的表面肌电信号的预处理包括依次进行的高通滤波预处理和去基线处理。

在上述的方法，通过下式的乘性迭代规则来进行初级nmf分解获取w矩阵和h矩阵的初值：

通过下式的加性迭代规则进行二级nmf分解求得所需的协同元矩阵：

在上述的方法，协同激活模型采用头脑风暴优化算法进行模型参数的优化选择。

本发明另一方面的目的在于提供一种基于肌肉协同的上肢多关节同步比例肌电控制系统，包括处理器；存储器；以及计算机程序；其中，所述计算机程序被存储在所述存储器中，并且被配置为由所述处理器执行，所述计算机程序包括用于执行上述技术方案所述的方法的指令。

本发明针对上肢多关节运动解析问题，提供了一种半监督式非负矩阵分解协同分析方法，能够实现多关节运动有效解耦，并获取更加稳定的激活系数序列。本发明协同元表征的是腕关节、肘关节两个屈/伸自由度，每个自由度用两个协同元来表示，因此一共是4个协同元，但这4个协同元与两个关节屈/伸自由度上所有的运动都是相关的，相比于传统代表离散动作的协同元，更能体现连续运动的思想。采用头脑风暴算法优化的svr构建激活系数序列和关节角度信息的协同激活模型，能够实现上肢多关节同步连续运动估计，且运动估计效果更加精确稳定。结合多自由度并行比例肌电控制策略，搭建上肢多关节同步比例肌电控制系统，能够为智能设备提供更加流畅稳定的控制指令，且控制系统的交互性更强。

附图说明

图1为本发明基于肌肉协同的上肢多关节同步比例肌电控制方法的整体流程图。

图2为半监督式nmf协同分析方法流程图。

图3为基于肌肉协同的上肢多关节同步比例肌电控制系统结构框图。

具体实施方式

下面结合附图和具体实施方式对本发明作进一步详细的说明。

图1示出了基于肌肉协同的上肢多关节同步比例肌电控制方法的整体流程，利用肌肉协同收缩理论在解决多自由度运动问题的优势，提供一种半监督式的肌肉协同分析方法提取运动相关激活系数序列，并通过支持向量回归(supportvectorregression,svr)构建协同激活模型同步估计出上肢腕关节和肘关节的运动角度信息，最后结合多自由度并行比例控制策略搭建上肢多关节同步比例肌电控制系统，为智能设备提供流畅稳定的运动控制指令。

一种基于肌肉协同的上肢多关节同步比例肌电控制方法，包括：

s1.表面肌电信号的采集及预处理

对于上肢腕关节和肘关节的单独屈伸运动和同步屈伸运动过程中的表面肌电信号进行采集，并对采集的信号进行预处理。腕关节单独屈伸运动定义为dof-1，肘关节单独屈伸运动定义为dof-2，两个关节的同步屈伸运动定义为dof-1&2。表面肌电信号的采集过程中肌电电极放置在作用于腕关节和肘关节屈伸运动的四块肌肉：桡侧腕长伸肌、桡侧腕屈肌、肱二头肌和肱三头肌，获取四通道肌电信号。对采集的表面肌电信号进行高通滤波预处理，之后进行去基线处理，即每个通道的肌电信号减去手臂放松状态下对应通道的肌电信号初值。

s2.提供一种半监督式非负矩阵分解方法进行肌肉协同分析，具体包括以下子步骤：

s21.肌肉激活度提取对预处理之后的表面肌电信号计算积分肌电值(integratedemg,iemg)作为肌肉激活度矩阵f(e)，其中n表示的是分析窗口的窗口长度，xi(i＝1,2,...,n)表示分析窗口中第i个semg信号的采样值。计算得到的特征矩阵f(e)作为肌肉激活度用于后续nmf分解来得到肌肉协同元矩阵w和就激活系数矩阵h，计算肌肉激活度时采用非叠加的滑动时间窗，窗口长度可设定为100ms。

s22.协同元矩阵训练

根据肌肉协同理论构建肌肉协同模型f(e)n×t＝wn×k×h(t)k×t，式中f(e)n×t为运动过程中的肌肉激活度，wn×k为肌肉协同矩阵，hk×t为激活系数矩阵，n是肌肉通道数，t是信号的样本点数，k是协同元的个数。

利用级联的非负矩阵分解方法(non-negativematrixfactorization,nmf)求解协同元矩阵，通过下式的乘性迭代规则来进行初级nmf分解获取w矩阵和h矩阵的初值：

通过下式的加性迭代规则进行二级nmf分解求得所需的协同元矩阵：

本发明将肌肉协同模型和nmf算法结合，根据肌电信号特征f(e)与肌肉协同矩阵w，以及激活系数矩阵h之间的关系f(e)＝wh+e分解得到w和h，其中e代表的是分解误差。假设分解误差服从高斯分布，那么应用最大似然函数可得到目标函数为使目标函数最小，可通过梯度下降法得到w和h的乘性迭代法则，如前两个公式所示。假设分解误差服从泊松分布，目标函数为类似得到迭代法则如后两个公式所示。本发明取乘性迭代规则(前两公式)来进行初级nmf分解获取w矩阵和h矩阵的初值，选取加性迭代规则进行二级nmf分解求得所需的肌肉协同矩阵及其系数矩阵。

本发明研究对象为肘、腕两个关节的背屈/跖屈自由度，设定每个关节的屈伸自由度对应两个协同元，那么共需要4个肌肉协同元来参与研究，例如，对单独腕关节运动(dof-1)时的肌肉激活度进行nmf分解，得到两个肌肉协同元w1,w2；对单独肘关节运动(dof-2)时的肌肉激活度进行nmf分解，得到两个肌肉协同元w3,w4。wi＝[w1i,w2i,…,wni],(i＝1,2,3,4)，wni(n＝1,2,...,n)代表第n块肌肉对第i个肌肉协同元wi的贡献度，即这块肌肉在这个协同元中所占的比重。

本发明所强调的创新点在于对于一名受试者通过nmf分解得到的肌肉协同元的稳定性，同时要求这种半监督式的nmf具有较好的泛化能力，即对多名受试者都具有很好的效果，因此可以提供一个分解结果的“参考”，即监督元。所以，本发明选取多名受试者的肌电信号数据样本提取每个人的肌肉协同元，并且每个协同元wi都是一个列向量，这样按列组合，就能得到每名受试者的肌肉协同元矩阵[w1,w2,w3,w4]。再将结果得到的多组肌肉协同元矩阵取均值，将均值作为监督元，w＝[w1(mean),w2(mean),w3(mean),w4(mean)]。

s23.激活系数序列提取

对于腕关节或肘关节单独运动以及两个关节同步运动时的表面肌电信号，按照步骤s21获取其肌肉激活度f(e)，利用步骤s22中求得的总体肌肉协同元矩阵w，结合下式moore-penrose广义逆求解方法，即可提取出激活系数序列：

h(t)＝w⁺·f(e)

其中，w⁺为协同元矩阵w的m-p逆。

s3.构建协同激活模型实现上肢多关节同步连续运动估计

通过支持向量回归(supportvectorregression,svr)构建激活系数序列和关节角度信息的协同激活模型，模型类型选择为ε-svr，核函数选为rbf核函数，并采用头脑风暴优化算法(brainstormoptimization，bso)对模型参数进行优化选择。本发明对上肢两个关节(腕关节和肘关节)分别建立一个bso优化的svr协同激活模型，来同步估计两个关节的运动角度信息。

在模型训练阶段，同步采集关节运动过程中的semg信号和角度信号并进行预处理，其中semg信号的采集及预处理参照步骤s1，角度信号则采用关节角度传感器采集，进行滑动平均滤波之后，再将其降到10hz实现与肌电信息的同步。按照步骤s21提取肌肉激活度矩阵f(e)，按照步骤s22获取协同元矩阵w，按照步骤s23提取出激活系数序列h(t)。对激活系数序列进行1hz低通滤波预处理后将其归一化到[0,1]区间，作为svr协同激活模型的输入；对相应关节的实测角度值进行相同区间的归一化处理，作为目标值进行模型的训练。

在运动估计阶段，对于在线采集的semg信号，按照所述协同分析方法步骤s21提取肌肉激活度矩阵，直接使用模型训练阶段获取的协同元矩阵w，按照步骤s23提取出激活系数序列h(t)，将其进行与训练阶段相同的预处理和归一化，输入到训练好的svr协同激活模型中；按照训练阶段的实测角度归一化参数对输出值进行反归一化，即可得到估计角度值。

利用多自由度并行比例肌电控制策略将估计出的所述角度信息转化为智能设备多个自由度的连续运行位移，具体公式如下：

sj＝αj·θj+βj(j＝1,2)

其中，sj是智能设备第j个自由度的目标位移，αj是相应的比例系数，θj为人体上肢运动过程中估计出的两个关节角度，其中j＝1代表腕关节，j＝2代表肘关节，βj则表示修正系数。

然后可将位移数据输出(例如串口通信方式)至智能设备控制系统为智能设备提供流畅稳定的运动控制指令。

本实施例中的智能设备可以是康复机器人、智能假肢等。

一种基于肌肉协同的上肢多关节同步比例肌电控制系统，包括处理器；存储器；以及计算机程序；其中，所述计算机程序被存储在所述存储器中，并且被配置为由所述处理器执行，所述计算机程序包括用于执行如上述实施例所述的上肢多关节同步比例肌电控制方法的指令。

上述实施例中的基于肌肉协同的上肢多关节同步比例肌电控制方法及系统可以是集成在智能设备内的，也可以是独立的系统，当作为独立的系统时系统框架图可参考图3。