[发明专利]基于肌电模型和无迹卡尔曼滤波的关节运动估计方法

说明书

本发明涉及到一种基于肌电模型和无迹卡尔曼滤波的关节运动估计方法，首先采集膝关节在连续运动状态下股二头肌、股四头肌、股外侧肌、股内侧肌、半腱肌、股薄肌的肌电信号和实时角度，对其进行带通滤波处理，并提取小波系数和均方根特征，然后使用一种结合了肌肉动力学、关节动力学、骨骼动力学和相关肌电特征的状态空间肌电模型，通过无迹卡尔曼滤波算法，得出Sigma采样集χ_i和权重W_i，然后进行进一步的预测，计算出系统状态变量和协方差矩阵P(k+1|k)，迭代循环后，实现对膝关节连续运动的估计。该方法与传统的角度估计方法相比，减小了系统误差、累积误差和外部干扰的影响，精度高，稳定性好，对目标机动反应快速，有了明显的改进。

技术领域

本发明属于模式识别领域，涉及一种肌电信号模式识别方法，特别涉及一种基于状态空间肌电模型和无迹卡尔曼滤波的关节连续运动估计方法。

背景技术

表面肌电信号(surface Electromyography,sEMG)是现今比较热门的前沿科学技术人机交互的输入信号源，它是一种非平稳的微弱信号，是肌肉兴奋时与其相关的运动单元共同产生的一组动作电位序列，它的特征区分明显，蕴含的信息丰富，且采集简单无创伤，成为目前人机交互技术中热门的研究领域。对于表面肌电信号的研究主要集中在特征提取和模式识别这两个过程。相应的研究成果也比较成熟，可以识别多个离散的动作类别。但在康复医疗机器人等领域更多时候需要的是对患者连续运动变量的预测，以实现对康复机器人的平滑柔顺控制。

传统的关节连续运动估计方法有提取肌电特征，然后通过神经网络建立sEMG与连续运动的回归模型。生理肌肉模型也是估计连续关节运动的一种方法。Buchanan等人提出了一种基于肌电信号的正向动力学模型，该模型由希尔肌肉模型(Hill muscle model,HMM)、肌肉激活动力学和关节正向动力学组成。该模型涉及多个生理参数，计算困难，实际应用受到限制。HMM是估计连续关节运动最常用的肌肉模型，但存在两个问题：一是HMM涉及到许多难以识别的复杂的生理参数，计算负担也比较大；二是HMM可以直接从sEMG信号计算关节力矩，但如果需要连续的关节运动估计，则还需要从力矩中计算运动状态。这通常会带来累计误差，降低预测精度。

而将HMM与关节正向动力学结合起来，并对模型参数进行简化替代的方法能够有效解决上述问题，该方法不需要计算关节力矩，而是可以直接从sEMG信号中计算关节运动。同时使用肌电特征组成测量方程作为反馈，并使用闭环预测算法，可较为准确的估计关节连续运动。

发明内容

本发明涉及到一种用于关节角度估计的状态空间肌电模型和无迹卡尔曼滤波方法，首先采集膝关节在连续运动状态下股二头肌、股四头肌、股外侧肌、股内侧肌、半腱肌、股薄肌的肌电信号和实时角度，对其进行带通滤波处理，并提取小波系数和均方根特征，然后使用一种结合了肌肉动力学、关节动力学、骨骼动力学和相关肌电特征的状态空间肌电模型，通过无迹卡尔曼滤波算法，得出Sigma采样集χ_i和权重W_i，然后进行进一步的预测，计算出系统状态变量和协方差矩阵P(k+1|k)，迭代循环后，实现对膝关节连续运动的估计。

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

步骤一，采集关节连续运动时相关肌肉的肌电信号，即通过肌电信号采集仪采集关节运动时相关肌肉的肌电信号，然后采用带通滤波的方法对原始信号进行预处理。

步骤二，根据希尔肌肉模型和关节动力学求出状态空间肌电模型的非线性表达式；该状态空间肌电模型首先对希尔肌肉模型进行参数代替及简化处理，简化处理后的离散时间的预测模型如下：

T_s是采样时间，是k时刻的关节角速度，θ_k是k时刻的关节位置，s_i是代替参数，均为常数。