一种基于微控制器的肌电信号处理及特征提取方法技术

本发明专利技术公开了一种基于微控制器的肌电信号处理及特征提取方法，包括信号采集放大处理、窗函数移动平均处理、阈值处理和波形长度特征提取四个步骤；其中，信号采集放大处理过程由硬件完成，采用三点式差动输入电极采集表面肌电信号；窗函数移动平均处理过程在极小运算量下对信号做平滑处理；阈值处理过程在极小的运算量下从肌电信号中提取出肌肉的收缩-舒张状态；波形特征提取过程在极小的运算量下获得波形变化幅度特征。本发明专利技术方法可以应用于运算资源与存储资源极为有限的微控制器上，利用较低成本采集处理肌电信号，获得肌肉的收缩-舒张状态与波形变化幅度特征。

A method of EMG signal processing and feature extraction based on Microcontroller

The invention discloses an EMG signal processing and feature extraction method based on the micro controller, including signal acquisition, amplification processing window moving average processing, threshold processing and waveform length feature extraction of four steps; wherein, the amplified signal acquisition processing by the hardware is completed by three point differential input electrode surface EMG signal acquisition window function; moving average process in minimal computation of signal smoothing; threshold processing in little calculation quantity from EMG signals extracted from muscle contraction diastolic state; waveform feature extraction process to obtain the amplitude characteristic waveform changes in the amount of computation under the minimum. The method of the invention can be applied to the micro controller with extremely limited computing resources and storage resources, and can obtain the characteristics of the muscle contraction relaxation state and the waveform change amplitude by using the low cost acquisition and processing of the EMG signal.

【技术实现步骤摘要】

本专利技术属于计算机协同医疗
，具体涉及一种基于微控制器的肌电信号处理及特征提取方法。
技术介绍
表面肌电信号(SurfaceElectromyogram，SEMG)是一种重要的人体生物信号，是通过人体表面电极，在人体皮肤表面采集人体肌肉运动时所产生的电位信号。它的来源是人体自主运动时神经肌肉活动发放的生物电信号，这些电信号沿肌纤维传播并经由皮肤、脂肪构成的容积导体滤波后在皮肤表面检测电极处的时间和空间上综合叠加，最终形成表面肌电信号。由于不同的动作引发不用的肌肉群作用，产生的肌电信号就不同，所以肌电信号在识别人体动作方面具有独到的优势，被广泛应用于假肢控制、康复训练、临床医学和体育科学等诸多领域中。表面肌电信号的频段主要分布在20～500Hz的区间内，主要能量集中在50～150Hz的区间内，基本服从正态分布。可以发现，人体肌电信号具有以下两方面特征：一是信号弱。直接从人体中检测出来的电信号的幅值往往较弱；在振幅上，由于肌电信号是一种肌肉产生而在皮肤表面提取的电信号，其电压信号非常微弱，峰值一般在0～6mV。因此，在处理肌电信号之前，需要使用硬件对信号进行放大处理。二是频率范围较低。经频谱分析可知，人体肌电信号的频率范围一般较小，在10Hz级，这使得在采样时对AD转换器的性能、MCU处理能力和通信模块传输能力没有较高的要求，远低于声音信号、图像信号的采集要求，对操作系统的实时性要求也较低。但是，微控制器的计算能力普遍较弱，无法使用频域变换、小波变换等需要较多计算和存储资源的方法来处理采集得到的原始肌电信号，也无法对肌电信号进行相应的特征提取。
技术实现思路
为了在计算和存储能力较弱的微控制器上，以较低成本处理肌电信号并提取特征，本专利技术提供了一种基于微控制器的肌电信号处理及特征提取方法，其采用的技术方案如下：首先，采用三点式差动输入电极采集表面肌电信号，并对该肌电信号进行放大以及整流平滑处理后由微控制器的ADC(模数转换)模块进行离散化采样；然后，采用窗函数对离散化后的肌电采样信号进行移动平均处理；最后，基于窗函数移动平均处理后的肌电采样信号进行特征提取，获得肌肉的收缩及舒张状态以及肌电信号波形变化幅度。本专利技术以P个采样周期大小的窗口对离散化后的肌电采样信号进行窗函数移动平均处理，其采用的窗函数表达式如下：signalMAV(i)=1PΣj=i-P+1isignal(j)]]>其中：signal(j)为离散化后的肌电采样信号中第j时刻的信号值，signalMAV(i)为经窗函数移动平均处理后的肌电采样信号中第i时刻的信号值，P为大于1的自然数，i为自然数，j为自然数且i-P+1≤j≤i。本专利技术通过对窗函数移动平均处理后的肌电采样信号进行阈值处理，获得肌肉的收缩及舒张状态。根据以下表达式对窗函数移动平均处理后的肌电采样信号进行阈值处理：signalthreshold(i)=1signalMAV(i-k)≥TH,k=0~T0signalMAV(i-k)<TH,k=0~Tnullothers]]>其中：signalthreshold(i)为经阈值处理后的肌电采样信号中第i时刻的信号值，signalMAV(i-k)为经窗函数移动平均处理后的肌电采样信号中第i-k时刻的信号值，TH为设定的阈值，i为自然数，k为自然数且0≤k≤T，T为设定的保持周期数且T为大于1的自然数，signalthreshold(i)＝1表示第i时刻肌肉处于收缩状态，signalthreshold(i)＝0表示第i时刻肌肉处于舒张状态。本专利技术通过对窗函数移动平均处理后的肌电采样信号进行波形长度特征提取，获得肌电信号波形变化幅度。根据以下表达式对窗函数移动平均处理后的肌电采样信号进行波形长度特征提取：signalWL(i)=1NΣp=i-N+1i(signalMAV(p-1)-signalMAV(p))2+12]]>其中：signalWL(i)为第i时刻肌电信号波形的线段长度和，signalMAV(p-1)和signalMAV(p)分别为经窗函数移动平均处理后的肌电采样信号中第p-1时刻和第p时刻的信号值，N为统计波形长度的采样周期数且N为大于1的自然数，i为自然数，p为自然数且i-N+1≤p≤i。由此可见，本专利技术基于上述技术方案可以应用于运算资源与存储资源极为有限的微控制器上，以较低的硬件成本采集处理肌电信号，获得肌肉的收缩-舒张状态与波形变化幅度特征。附图说明图1为本专利技术肌电信号处理及特征提取的流程示意图。图2为三点式差动输入电极的贴放位置示意图。具体实施方式为了更为具体地描述本专利技术，下面结合附图及具体实施方式对本专利技术的技术方案进行详细说明。本实施例采用FreescaleKL02作为微控制器，以C语言编写控制程序，采用CyborganOSCore作为微控制器实时操作系统，以肌电信号特征最明显的肱二头肌作为被测肌肉，采用三点式差动输入电极对肌电信号进行采集，使用AdvancerTechnology的MuscleSensorPlatinumv3.3对信号做硬件处理，如图1所示，通过躯干电信号处理服务对肌电信号依次进行信号采集、信号放大处理、信号整流平滑处理、AD采样离散化处理、窗函数移动平均处理、阈值处理以及波形长度特征提取，最后通过蓝牙将数据包发往CyborganOSBrain。其中信号采集按图2所示方式将MuscleSensorPlatinumv3.3的三点式差动输入电极贴放在肱二头肌对应皮肤表面上，将MuscleSensorPlatinumv3.3的SIG信号接入到FreescaleKL02的ADC0_SE0引脚上，将MuscleSensorPlatinumv3.3与FreescaleKL02开发板共地。使用CyborganOS配置工具配置一个定时器OsAlarm_Besp，该定时器自动触发，每1ms触发一次，用来周期性调用任务OsTask_Besp。任务OsTask_Besp首先启动ADC，完成一个数模转换，将获取到的10位电压信号存放于一个uint16类型的变量signal内，并以此执行窗函数移动平均处理、阈值处理、波形长度特征提取三个处理过程，并将结果用通信服务传递给上位机，具体过程如下：(1)窗函数移动平均处理。本实施方式使用窗函数对信号进行平均处理，公式如下，表示在P个采样周期窗口内对信号瞬本文档来自技高网...

【技术保护点】
一种基于微控制器的肌电信号处理及特征提取方法，其特征在于：首先，采用三点式差动输入电极采集表面肌电信号，并对该肌电信号进行放大以及整流平滑处理后由微控制器的ADC模块进行离散化采样；然后，采用窗函数对离散化后的肌电采样信号进行移动平均处理；最后，基于窗函数移动平均处理后的肌电采样信号进行特征提取，获得肌肉的收缩及舒张状态以及肌电信号波形变化幅度。

【技术特征摘要】
1.一种基于微控制器的肌电信号处理及特征提取方法，其特征在于：
首先，采用三点式差动输入电极采集表面肌电信号，并对该肌电信号进行
放大以及整流平滑处理后由微控制器的ADC模块进行离散化采样；
然后，采用窗函数对离散化后的肌电采样信号进行移动平均处理；
最后，基于窗函数移动平均处理后的肌电采样信号进行特征提取，获得肌
肉的收缩及舒张状态以及肌电信号波形变化幅度。
2.根据权利要求1所述的肌电信号处理及特征提取方法，其特征在于：以
P个采样周期大小的窗口对离散化后的肌电采样信号进行窗函数移动平均处理，
其采用的窗函数表达式如下：
signalMAV(i)=1PΣj=i-P+1isignal(j)]]>其中：signal(j)为离散化后的肌电采样信号中第j时刻的信号值，signalMAV(i)
为经窗函数移动平均处理后的肌电采样信号中第i时刻的信号值，P为大于1的
自然数，i为自然数，j为自然数且i-P+1≤j≤i。
3.根据权利要求1所述的肌电信号处理及特征提取方法，其特征在于：通
过对窗函数移动平均处理后的肌电采样信号进行阈值处理，获得肌肉的收缩及
舒张状态。
4.根据权利要求3所述的肌电信号处理及特征提取方法，其特征在于：根
据以下表达式对窗函数移动平均处理后的肌电采样信号进行阈值处理：
signalthreshold(i)=1signalMAV(i-k)≥TH,k=0~T...

【专利技术属性】
技术研发人员：李红，邵开来，王杰，杨国青，吴朝晖，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江;33