基于表面肌電的肌肉疲勞檢測(cè)系統(tǒng)研究

糜超 陳陽(yáng) 鄒凌

摘 要： 為了利用表面肌電信號(hào)自動(dòng)檢測(cè)肌肉疲勞狀態(tài)，設(shè)計(jì)一套8導(dǎo)聯(lián)表面肌電采集系統(tǒng)。該系統(tǒng)主要包括模擬前端放大電路、主控模塊和電源模塊三部分，并研究了肌肉疲勞與正常狀態(tài)下表面肌電信號(hào)的特征提取及分類(lèi)方法。實(shí)驗(yàn)表明：自行設(shè)計(jì)的放大器分辨率為38.1 μV，系統(tǒng)噪聲小于15.3 μV，采樣率達(dá)到1 000 Hz，功耗約為30.7 mW；使用該放大器采集肌肉疲勞與正常狀態(tài)下肌電信號(hào)，支持向量機(jī)分類(lèi)的識(shí)別率高達(dá)98%。該系統(tǒng)具有體積小、精度高、功耗低及操作簡(jiǎn)便等優(yōu)點(diǎn)，可用于家庭健康監(jiān)護(hù)領(lǐng)域。

關(guān)鍵詞： 表面肌電信號(hào)； 肌肉疲勞； 特征提取； 支持向量機(jī)； 放大電路； 健康監(jiān)護(hù)

中圖分類(lèi)號(hào)： TN911.23?34； TP368.1； R318.6 文獻(xiàn)標(biāo)識(shí)碼： A 文章編號(hào)： 1004?373X（2018）20?0078?05

Abstract： An 8?channel sEMG acquisition system is designed to automatically detect the muscle fatigue status by using sEMG signals. The system mainly consists of the analog front?end amplifying circuit， main control module and power supply module. The feature extraction and classification method of the sEMG signal at the muscle fatigue status and normal status are studied. The experimental results show that the designed amplifier has resolution rate of 38.1 μV， system noise of smaller than 15.3 μV， sampling rate of as high as 1 000 Hz， and power consumption of about 30.7 mW； when the amplifier is used to collect the sEMG signal at the muscle fatigue status and normal status， the recognition rate of support vector machine （SVM） classification is as high as 98%； the system has the advantages of small size， high precision， low power consumption and easy operation， and can be used in the field of family health monitoring.

Keywords： sEMG signal； muscle fatigue； feature extraction； SVM； amplifying circuit； health monitoring

0 引 言

肌肉疲勞是當(dāng)肌肉不能滿(mǎn)足某一力量要求時(shí)的肌肉狀態(tài)。連續(xù)強(qiáng)烈的肌肉收縮是導(dǎo)致肌肉疲勞的主要原因[1]，如帕金森疾病、內(nèi)分泌紊亂、癌癥、營(yíng)養(yǎng)不良，以及長(zhǎng)時(shí)間使用肌肉電刺激儀等都會(huì)引起肌肉疲勞。

表面肌電信號(hào)（surface Electromyography，sEMG）是人體肌肉發(fā)生收縮活動(dòng)時(shí)產(chǎn)生的一種生理電信號(hào)，經(jīng)常被用來(lái)評(píng)估肌肉活動(dòng)狀態(tài)。1890年Etienne Jules Marey第一次記錄了肌肉活動(dòng)的電信號(hào)。肌電信號(hào)是一種隨機(jī)生理電信號(hào)，具有以下特點(diǎn)[2?3]：

1） 微弱性。肌電信號(hào)幅度范圍一般在100～5 000 mV，均方根在0～1.5 mV；相比于針電極檢測(cè)的肌電信號(hào)，表面肌電信號(hào)的幅度更小。

2） 低頻特性。表面肌電信號(hào)能量主要集中在1 000 Hz以下，頻譜主要集中在20～500 Hz之間，功率譜通常處在30～300 Hz之間。

3） 高阻抗性。當(dāng)皮膚干燥、無(wú)損的情況下，人體的電阻可高達(dá)40～400 kW。在肌肉處于疲勞狀態(tài)下時(shí)，功率譜密度會(huì)向低頻段漂移及在峰值上產(chǎn)生變化[4]。比如，在進(jìn)行最大力負(fù)重練習(xí)后，肌電信號(hào)的峰值會(huì)變大，并且向低頻段漂移。

肌電是目前肌肉疲勞檢測(cè)中較為準(zhǔn)確方便的測(cè)量方法。作為一種簡(jiǎn)單、無(wú)創(chuàng)傷、可定量的研究方法，利用肌電圖研究肌肉疲勞是勞動(dòng)生理學(xué)中常用的方法。表面肌電放大器的設(shè)計(jì)和基于表面肌電的肌肉疲勞檢測(cè)已經(jīng)取得了一定的進(jìn)展，如周兵等人設(shè)計(jì)了基于ARM 7的表面肌電信號(hào)采集儀[5]；Koutsos E等采集表面肌電信號(hào)后，進(jìn)行功率譜上的分析[6]。目前的檢測(cè)大都是基于觀測(cè)功率譜圖這一方法，需要一定的專(zhuān)業(yè)知識(shí)才能鑒別出肌肉狀態(tài)，不便于家庭化的推廣使用。

本文設(shè)計(jì)了一個(gè)高精度、便攜式、8導(dǎo)聯(lián)的表面肌電采集系統(tǒng)，利用該系統(tǒng)分別采集不同被試肌肉正常與疲勞狀態(tài)下的肌電信號(hào)。提取信號(hào)的有效特征，進(jìn)行肌肉正常與疲勞狀態(tài)分類(lèi)，實(shí)現(xiàn)肌肉疲勞的檢測(cè)。

1 肌電采集系統(tǒng)硬件設(shè)計(jì)

多通道肌電采集放大電路主要由模擬采集放大部分、主控芯片部分、電源部分三大塊組成。模擬前端主要對(duì)微弱的肌電信號(hào)進(jìn)行兩次放大以及帶通濾波，拾取有效的肌電信號(hào)；經(jīng)過(guò)放大濾波的模擬肌電信號(hào)傳送給24位高精度的ADS1278模/數(shù)轉(zhuǎn)換芯片（Analog?to?Digital converter，ADC），得到穩(wěn)定的數(shù)字信號(hào)；主控芯片將數(shù)字信號(hào)進(jìn)行存儲(chǔ)，同時(shí)發(fā)送給上位機(jī)，以便后期研究。

1.1 模擬前端電路設(shè)計(jì)

模擬前端需要對(duì)信號(hào)進(jìn)行調(diào)理，將微弱的肌電信號(hào)放大至ADC可識(shí)別的電壓范圍之內(nèi)，同時(shí)最大程度地去除干擾噪聲，保留肌電信號(hào)有效成分。在設(shè)計(jì)模擬前端時(shí)，需要根據(jù)表面肌電信號(hào)的特性，完成共模抑制比、放大倍數(shù)等參數(shù)的計(jì)算。

假設(shè)肌電電壓的分辨率ΔVEMG為500 nV，當(dāng)滿(mǎn)足式（10）時(shí)，最小增益[GAFEmin]為4.2。其中，LSB為ADC精度。

以上計(jì)算基于以下假設(shè)：

1） 放大器的量程FSR為5 V；

2） 電路中做了去直流處理；

3） 放大器具有大于75 dB的共模抑制比（50 Hz處）；

4） 結(jié)合ADS1278芯片手冊(cè)，假設(shè)信噪比失真參數(shù)SINAD為103 dB。綜上所述，前端模擬電路需要滿(mǎn)足如下要求：

綜合以上計(jì)算分析，本系統(tǒng)使用差分輸入的儀表放大器INA333進(jìn)行5倍的前級(jí)放大，并且同時(shí)對(duì)信號(hào)高通濾波，截止頻率為6 Hz。為了有效去除直流成分，使用OPA333運(yùn)算放大器進(jìn)行去直流操作。使用差分輸出放大器THS4521對(duì)信號(hào)進(jìn)行約40倍的后級(jí)放大，同時(shí)進(jìn)行截止頻率為500 Hz的低通濾波。模擬前端電路設(shè)計(jì)參數(shù)及原理圖如表1、圖1所示。

1.2 電源及電平轉(zhuǎn)換電路設(shè)計(jì)

肌電信號(hào)采集系統(tǒng)需要低噪聲、高精度的供電模塊，本系統(tǒng)電源設(shè)計(jì)如圖2所示。ADC芯片的參考電壓直接決定了工作的精度和穩(wěn)定性。本系統(tǒng)中選用ISL21007?25（U3）穩(wěn)壓芯片來(lái)實(shí)現(xiàn)2.5 V的參考電壓值。該芯片具有低噪音（峰峰值4.5 mV），高精度（±0.5 mV），低漂移（3 ppm/℃）的特點(diǎn)。前端的模擬采集放大部分的供電使用±2.5 V的雙電源供電。2.5 V的電源使用具有低壓差線(xiàn)性穩(wěn)壓器的TPS79625（U2）實(shí)現(xiàn)；-2.5 V的電源首先使用高頻（2 MHz）逆向電源芯片MAX889T（U4）得到-3 V的電壓，再使用具有低噪音、高電源供電抑制比的TPS72325（U5）芯片得到-2.5 V的電壓。前端濾波和ADC內(nèi)部供電使用5 V的單電源供電，使用TPS79650（U1）芯片實(shí)現(xiàn)。數(shù)字信號(hào)電源部分和ADC的I/O口3.3 V供電，以及ADC的數(shù)字內(nèi)核1.8 V供電使用TPS79645（U6）實(shí)現(xiàn)。

1.3 ADC選擇及主控模塊設(shè)計(jì)

ADC芯片的選擇需要滿(mǎn)足如下要求：具有能夠連續(xù)采樣的多通道；支持使用SPI口與主控進(jìn)行通信；較低的功耗。

基于以上幾點(diǎn)，本系統(tǒng)采用TI公司的ADS1278[8]芯片。該芯片能夠達(dá)到144 KSPS的采樣率，111 dB的信噪比，具有8個(gè)通道。

主控模塊采用STM32實(shí)現(xiàn)，主要功能是利用SPI接收來(lái)自ADC芯片的數(shù)字信號(hào)，將接收到的數(shù)據(jù)封裝成包，數(shù)據(jù)包大小為32 bit，包含一個(gè)8 bit的通道號(hào)和24 bit的肌電數(shù)據(jù)，并且將其通過(guò)USB口發(fā)送給上位機(jī)。主控程序的主要邏輯分為SPI總線(xiàn)控制、FIFO數(shù)據(jù)保存、USB數(shù)據(jù)發(fā)送三部分，如圖3所示。主控芯片通過(guò)SPI控制ADC芯片，設(shè)置定時(shí)器查詢(xún)ADC是否完成數(shù)據(jù)轉(zhuǎn)換，完成則循環(huán)讀取8次DOUT口數(shù)據(jù)，每次讀取1 bit，將STRB標(biāo)志置為1。當(dāng)主程序檢測(cè)到STRB置為1后，檢查FIFO是否已滿(mǎn)，若滿(mǎn)則將TH_FULL標(biāo)志置為1；未滿(mǎn)則將已讀取的8 bit數(shù)據(jù)保存至FIFO。主程序檢測(cè)到TH_FULL為1后，則將FIFO里的數(shù)據(jù)通過(guò)USB口發(fā)送給與之連接的PC端，同時(shí)清空FIFO。

1.4 放大器主要性能指標(biāo)

本文設(shè)計(jì)實(shí)現(xiàn)的表面肌電放大器實(shí)物尺寸為20 cm×10 cm×2 cm，體積較小。設(shè)備使用鋰電池進(jìn)行供電，電極使用可插拔式的DIN2.0接口，數(shù)據(jù)利用USB傳輸。系統(tǒng)測(cè)量參數(shù)如表2所示。結(jié)果表明，自行設(shè)計(jì)的8導(dǎo)聯(lián)表面肌電放大器符合表2中的設(shè)計(jì)參數(shù)，具有較低功耗、較高共模抑制比和較高精度的特點(diǎn)。與國(guó)標(biāo)相比[9]，共模抑制比、有效帶寬達(dá)到了要求；系統(tǒng)噪聲未能達(dá)標(biāo)，但已經(jīng)滿(mǎn)足肌肉疲勞檢測(cè)的要求。

2 特征提取及分類(lèi)

2.1 特征提取

提取有效的特征是疲勞檢測(cè)的基礎(chǔ)。肌電信號(hào)的特征提取主要分為三類(lèi)方法：時(shí)域法，頻域法，時(shí)域?頻域結(jié)合法。表面肌電信號(hào)的時(shí)域特征與頻域特征或時(shí)?頻域特征相比，具有計(jì)算量簡(jiǎn)單、獲取迅速等優(yōu)點(diǎn)，可以獲得相對(duì)較好的分類(lèi)特性。本文提取平均絕對(duì)值 （Integrated Electromyogram，IEMG）、過(guò)零率（Zero Crossing，ZC）、符號(hào)改變斜率（Slope Sign Change，SSC）、波長(zhǎng)（Waveform Length，WL）、4個(gè)時(shí)序上的肌電信號(hào)特征。

1） 平均絕對(duì)值

2） 過(guò)零率

3） 符號(hào)改變斜率

4） 波 長(zhǎng)

2.2 分 類(lèi)

本文使用4個(gè)分類(lèi)算法來(lái)檢測(cè)肌肉的疲勞和非疲勞兩種狀態(tài)。其中，K近鄰（KNN）和支持向量機(jī)（SVM）是非參數(shù)方法；線(xiàn)性判別（LDA）和樸素貝葉斯（NB）為參數(shù)化方法。所有方法使用Python實(shí)現(xiàn)。

2.3 實(shí)驗(yàn)結(jié)果與分析

本實(shí)驗(yàn)采用自主設(shè)計(jì)的8導(dǎo)聯(lián)肌電采集系統(tǒng)，1 kHz的采樣率，設(shè)置放大倍數(shù)為200倍。擬選取平均年齡在25歲，平均身高在175 cm，平均體重在60 kg的被試10名。被試無(wú)神經(jīng)肌肉方面的疾病，以及沒(méi)有受過(guò)特殊的力量訓(xùn)練。實(shí)驗(yàn)要求被試在休息一段時(shí)間以后，將一次性心電電極放置在肱三頭肌附近，每個(gè)電極之間間隔10 mm。要求被試在站立狀態(tài)下，以2 s每次的速度舉啞鈴數(shù)次。實(shí)驗(yàn)中，首先采集被試非疲勞狀態(tài)下的肌電信號(hào)，每次持續(xù)3 s。在被試負(fù)重訓(xùn)練，手臂感到疲憊后，采集被試靜止?fàn)顟B(tài)下的肌電信號(hào)，每次持續(xù)3 s。將采集到的肌電信號(hào)首先進(jìn)行20～500 Hz帶通濾波，去除50 Hz工頻干擾[10]。肌肉正常狀態(tài)下和疲勞狀態(tài)下的肌電圖分別如圖4a）、圖4b）所示，功率譜對(duì)比如圖5所示。當(dāng)肌肉處于疲勞狀態(tài)下時(shí)，表面肌電功率譜圖的峰值變大，向低頻段漂移，與已有的相關(guān)文獻(xiàn)報(bào)道相符[4]。

提取第2.1節(jié)中的平均絕對(duì)值、過(guò)零率、符號(hào)改變率、波長(zhǎng)4個(gè)特征，使用K近鄰、支持向量機(jī)、線(xiàn)性判別和樸素貝葉斯分類(lèi)器分類(lèi)，分類(lèi)結(jié)果如表3所示。

SVM分類(lèi)時(shí)正確率達(dá)到了98%，能夠有效區(qū)分肌肉疲勞與正常兩種狀態(tài)。

3 結(jié) 語(yǔ)

本文研制了一套8導(dǎo)聯(lián)表面肌電采集系統(tǒng)。對(duì)10位被試肌肉疲勞與正常狀態(tài)下的表面肌電信號(hào)采集與分析，利用機(jī)器學(xué)習(xí)方法分類(lèi)。實(shí)驗(yàn)表明，本系統(tǒng)能很好地完成對(duì)肌肉疲勞的檢測(cè)任務(wù)，適用于家庭健康監(jiān)護(hù)等場(chǎng)景。為進(jìn)一步提高系統(tǒng)性能，將在信號(hào)采集與分類(lèi)算法上進(jìn)行改進(jìn)，使系統(tǒng)實(shí)現(xiàn)低功耗與小型化。

注：本文通訊作者為鄒凌。

參考文獻(xiàn)

[1] BOYAS S， GU?VEL A. Neuromuscular fatigue in healthy muscle： underlying factors and adaptation mechanisms [J]. Annals of physical & rehabilitation medicine， 2011， 54（2）： 88?108.

[2] NA Y， CHOI C， LEE H D， et al. A study on estimation of joint force through isometric index finger abduction with the help of SEMG peaks for biomedical applications [J]. IEEE transactions on cybernetics， 2016， 46（1）： 2?8.

[3] DASOG M， KOIRALA K， LIU P， et al. Electromyogram bandwidth requirements when the signal is whitened [J]. IEEE transactions on neural systems and rehabilitation engineering， 2014， 22（3）： 664?670.

[4] ROGERS D R， MACISAAC D T. EMG?based muscle fatigue assessment during dynamic contractions using principal component analysis [J]. Journal of electromyography and kinesiology， 2011， 21（5）： 811?818.

[5] 周兵，紀(jì)曉亮，張榮，等.基于嵌入式技術(shù)的表面肌電信號(hào)采集儀設(shè)計(jì)[J].現(xiàn)代電子技術(shù)，2010，33（6）：55?57.

ZHOU Bing， JI Xiaoliang， ZHANG Rong， et al. Design of surface EMG acquisition instrument based on embedded technology [J]. Modern electronics technique， 2010， 33（6）： 55?57.

[6] KOUTSOS E， GEORGIOU P. An analogue instantaneous median frequency tracker for EMG fatigue monitoring [C]// Proceedings of IEEE International Symposium on Circuits and Systems. Melbourne： IEEE， 2014： 1388?1391.

[7] VAN RIJN A C M， PEPER A， GRIMBERGEN C A. High?quality recording of bioelectric events [J]. Medical & biological engineering & computing， 1990， 28（5）： 389?397.

[8] Texas Instruments. ADS1278?HT： octal， simultaneous?sampling 24?bit analog?to?digital converter [EB/OL]. [2013?04?01]. http：//www.ti.com/lit/ds/symlink/ads1278?ht.pdf.

[9] Tianjin Quality Supervision and Inspection Center of Medical Device， CFDA. Myoelectric biofeedback equipment： YY/T 1095?2015 [S/OL]. [2015?03?02]. http：//www.168i.cn/soft/show.asp？id=218720.

[10] 湯烈，葛良全.一種濾除表面肌電信號(hào)中工頻干擾的新型方法[J].現(xiàn)代電子技術(shù)，2015，38（4）：30?34.

TANG Lie， GE Liangquan. A novel method of filtering the power?line interference in surface EMG signal [J]. Modern electronics technique， 2015， 38（4）： 30?34.