基于模糊近似熵与辛几何的肌肉疲劳分析.docVIP

基于模糊近似熵与辛几何的肌肉疲劳分析 　　摘 要 为探索人体肌肉在运动后疲劳程度的变化，本文运用非线性动力学中的模糊近似熵和辛几何2种方法，对12名健康受试者在腕关节负重60%最大自主性施力、直至力竭的情况下，采集静态表面肌电信号作为样本，进行处理分析。模糊近似熵是近似熵算法的改进，可以更好地衡量一个序列的复杂程度；辛几何则运用相空间重构原理，提取主分量谱。试验结果初步表明，应用上述2种方法可相互补充，对受试者的肌肉疲劳程度做出有效、定性的评估，而且2种方法对肌肉信号处理结果的一致性有助于进一步描述肌肉系统的生物力学特征。 关键词 表面肌电信号 肌肉疲劳 模糊近似熵 辛几何 香农熵 引 言 肌肉在持续收缩过程中， 会逐渐进入疲劳状态，由此记录下来的肌电信号就是疲劳肌电信号[1]。由于神经肌电活动方式将兴奋传递至肌肉引起肌电变化并伴随肌肉收缩， 探索运动时肌肉疲劳出现过程中的肌电活动变化规律及特点， 是监测运动性肌肉疲劳和研究疲劳机理必不可少的一种方法[2]。运动性肌肉疲劳是指运动引起肌肉产生最大收缩力量或者最大输出功率暂时性下降的生理现象[3]。而表面肌电信号是肌肉收缩时伴随的电信号，它蕴涵许多肌肉运动信息，是在体表无创检测肌肉活动的重要方法[4]。目前，国内外相关文献指出，当肌肉疲劳时，肌肉能够产生的力量将降低；而在疲劳性练习时，肌肉的活动和生理、生化特性将发生改变，从而导致sEMG信号发生相应的变化[5]。因此，通过对表面肌电信号的分析，来分析肌肉疲劳状态，在肌肉功能的评价中具有重要意义[6]。 对于疲劳肌电信号的研究，人们一直在探索着。过去，人们普遍利用肌电（EMG）信号幅值和频谱变化作为肌肉疲劳度指标，其中应用最为广泛的是中值频率法。另外也有人采用时频分析法、反射响应法、波形分析法和线性建模法等来研究肌肉的疲劳程度， 但这些方法由于受使用条件限制， 存在许多不稳定性， 因而没有得到推广应用[7]。而近年来，随着人们对神经肌肉系统非线性性质认识的深入，有人开始尝试采用非线性数学方法分析肌电信号[8]，从而使该领域研究步入一个崭新阶段。本文就将尝试运用非线性方法中的模糊近似熵（fApEn）和辛几何谱（SG）2种方法来分析肌电信号，从而为肌肉疲劳的界定和研究提供帮助和理论依据，并力图对肌电信号的物理解释做出一些探索性工作。 1 模糊近似熵和辛几何方法 1.1 模糊近似熵 近似熵（ApEn）是用一个非负数来表示一个时间序列的复杂性，越复杂的时间序列对应的近似熵越大[9]。它并不是企图完全重构吸引子，而是用一种有效的统计方法——边缘概率的分布来区分各种过程。近似熵表现出主要的特点是只需要比较短的数据就能估计出比较稳定的统计值[10]；具有较好的抗干扰能力；并且适用于随机信号，这样就适用于生理性时间序列的分析，特别是对肌电信号分析是非常有利的。 1.2 辛几何 辛几何（SG） 是一种相空间的几何学，是针对哈密顿体系提出的，而一切守恒真实物理过程都可以表示为哈密顿体系[12]。在对信号处理手段上，辛几何方法中需通过一系列变换得到“主分量谱”，通常对非线性序列求取嵌入维数时，常常会受到数据长度、延迟时间、噪声等主观因素的影响，而SG方法则受其干扰很小，而且它还能保持时间序列的基本性质不变，在非线性短时间序列处理上有很强的应用性[13]。 由于表面肌电信号是非线性、非平稳的[15]，在对其做近似处理时，如通常的分段截取来获得平稳信号等，必然要丢失一定的信息，为避免人为近似所带来有用信息的损失和冗余信息的增加，就可以采用目前国际上新出现的具有良好发展前景的信号处理方法：辛几何方法。 2 试验结果分析 2.1 对象与信号提取 所分析的肌电信号来自于12名健康受试者（8男4女，年龄：30.2±4.9岁）的肱二头肌在静态自愿等长收缩时记录下来的信号。受试者均没有神经肌肉疾病史，并且在每次实验前签署知情同意书。实验中使用一对表面肌电自粘式导电凝胶电极（直径25mm），先在肱二头肌最厚点处摩擦、清洁，然后立即纵向放置电极。肌电参考电极放置在尺骨近端的头部。实验开始时，要求受试者进行肘关节负重要达到他/她最大自主收缩（MVC）的60%，并且通过读取屏幕上的力矩来保持这一百分比，直至力竭。EMG的增益信号为1000，其带宽为10～400Hz。由肌电电极采集这时的信号，其采样频率为1kHz，并存储在计算机中以供进一步分析。实验共采集到12组数据，每名受试者一组，每组9个疲劳负荷肌电信号样本，为恢复期分别进行2s，4s，6、8s，10s，20s，30s，60s，120s的信号取样。 2.2 模糊近似熵处理结果 对于疲劳肌电的模糊近似熵fApEn（m，r，N），选择嵌入维数为m 2，而数据集N和阈值参数r根据疲劳肌电实验数据选择。以下r的取值均为r*SD（S