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

本发明专利技术公开了基于肌肉协同的上肢多关节同步比例肌电控制方法与系统，所述方法包括：首先对关节运动相关肌群的表面肌电信号进行采集并预处理，根据肌肉协同收缩模型，提供一种半监督式的非负矩阵分解协同分析方法对肌电信号进行解耦，有效提取关节运动相关的肌肉协同元及其激活系数序列；其次通过支持向量回归构建激活系数和关节角度的协同激活模型，并采用头脑风暴算法对模型参数进行优化选择，实现上肢多个关节运动角度信息的同步估计；最后结合多自由度并行比例肌电控制策略，构建基于肌肉协同的上肢多关节同步比例肌电控制系统，将估计的多关节运动角度信息转化为康复辅助设备多自由度的运行位移，为康复辅助设备提供流畅稳定的运动控制指令。

【技术实现步骤摘要】
基于肌肉协同的上肢多关节同步比例肌电控制方法与系统
本专利技术属于肌电信号处理及肌电控制领域，涉及肌肉协同分析方法和基于协同激活模型的多关节同步连续运动估计方法，特别涉及基于肌肉协同的上肢多关节同步比例肌电控制方法与系统。
技术介绍
表面肌电信号(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。表面肌电信号的采集过程中肌电电极放置在作用于腕关节和肘关节屈伸运动的四块肌肉：桡侧腕长伸肌、桡侧腕屈肌、肱二头肌和肱三头肌，获取四通道肌电信号。对采集的表面肌本文档来自技高网...

【技术保护点】
1.一种基于肌肉协同的上肢多关节同步比例肌电控制方法，其特征在于，包括：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，通过支持向量回归分别构建上肢腕关节和肘关节的激活系数序列和关节角度信息的协同激活模型，利用所述协同激活模型对上肢运动过程中腕关节和肘关节的角度进行同步估计。...

【技术特征摘要】
1.一种基于肌肉协同的上肢多关节同步比例肌电控制方法，其特征在于，包括：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，通过支持向量回归分别构建上肢腕关节和肘关节的激活系数序列和关节角度信息的协同激活模型，利用所述协同激活模型对上肢运动过程中腕关节和肘关节的角度进行同步估计。2.根据权利要求1所述的方法，其特征在于，所述协同激活模型的训练方法包括：同步采集腕关节和肘关节运动过程中的表面肌电信号和角度信号并进行预处理；提取肌肉激...

【专利技术属性】
技术研发人员：艾青松，张亚男，刘泉，陈昆，孟伟，左洁，
申请(专利权)人：武汉理工大学，
类型：发明
国别省市：湖北,42