孟慶豐，陶慶，2，來全寶，胡玉舸

（1.新疆大學(xué)機(jī)械工程學(xué)院，新疆 烏魯木齊 830047；2.新疆大學(xué)機(jī)械工程博士后研究中心，新疆 烏魯木齊 830047）

1 引言

肌肉疲勞是人們?nèi)粘Ｉ钪幸环N常見的生理現(xiàn)象，一般是指由長(zhǎng)時(shí)間運(yùn)動(dòng)或發(fā)力引起的肌無力或酸痛的感覺[1]。日常生活中我們或多或少都會(huì)遇到肌肉疲勞的現(xiàn)象，但是如果不顧肌肉的疲勞而過度運(yùn)動(dòng)或長(zhǎng)時(shí)間活動(dòng)的話，很可能會(huì)造成肌體的永久性損傷，更甚者會(huì)帶來生命危險(xiǎn)[2]，因此對(duì)于肌肉疲勞評(píng)價(jià)的研究足以引起我們的高度重視。對(duì)于肌肉疲勞評(píng)價(jià)方法的研究與分析，無論是在體育訓(xùn)練還是在康復(fù)醫(yī)學(xué)、運(yùn)動(dòng)醫(yī)學(xué)等領(lǐng)域都具有重要意義。

如今，已知通過記錄表面肌電信號(hào)計(jì)算疲勞指數(shù)的方法可以分為兩個(gè)主要類別[3]。第一類由經(jīng)典的基于信號(hào)可量化頻譜變化的算法組成。例如經(jīng)典的平均頻率以及頻譜矩比[4]和基于小波的方法[5]都屬于此類。第二類考慮信號(hào)的非線性成分。例如考慮信號(hào)重復(fù)模式或分形測(cè)量的遞歸量化分析和基于熵的方法的模糊近似熵都屬于此類[6]。

由于上肢的運(yùn)動(dòng)功能遠(yuǎn)比下肢復(fù)雜，上肢在人體日常的工作生活扮演著重要的角色，而肘關(guān)節(jié)的屈伸動(dòng)作甚至出現(xiàn)在人體運(yùn)動(dòng)的每個(gè)瞬間，并且上肢中肘關(guān)節(jié)的運(yùn)動(dòng)減少50%將導(dǎo)致上肢運(yùn)動(dòng)功能降低近80%[7]，肘關(guān)節(jié)的正常運(yùn)動(dòng)對(duì)人體的工作生活是非常重要的，因此探究更有效的肘關(guān)節(jié)肌肉疲勞評(píng)價(jià)方法迫在眉睫。

為此，國(guó)內(nèi)外學(xué)者已經(jīng)針對(duì)肘關(guān)節(jié)肌肉疲勞評(píng)價(jià)方法展開了大量的研究，不同方法得出的評(píng)價(jià)指標(biāo)獲得的肌肉疲勞評(píng)價(jià)效果不同。文獻(xiàn)[8]利用AEMG與MPF等線性和Lempel?Ziv 復(fù)雜度與DET等非線性指標(biāo)分析不同負(fù)載下肱二頭肌疲勞過程中表面肌電信號(hào)的變化，結(jié)果顯示疲勞指標(biāo)的變化具有非疲勞特異性。

文獻(xiàn)[9]采集不同狀態(tài)的上肢兩側(cè)肌電信號(hào)，將時(shí)域的積分肌電值、均方根振幅和頻域的平均功率頻率、中位頻率作為疲勞評(píng)價(jià)指標(biāo)，研究得出上肢一側(cè)近端肌肉疲勞會(huì)使另一側(cè)的運(yùn)動(dòng)信號(hào)產(chǎn)生伴隨性變化，即產(chǎn)生中樞共驅(qū)動(dòng)現(xiàn)象。

文獻(xiàn)[10]分析時(shí)域的積分肌電值、均方根振幅和頻域的平均功率頻率、中位頻率等指標(biāo)隨著疲勞程度加深的變化趨勢(shì)，提出一種肌肉疲勞出現(xiàn)時(shí)刻的估計(jì)方法。文獻(xiàn)[11]將時(shí)域的iEMG、RMS與頻域的MPF、Fc 分析評(píng)價(jià)指標(biāo)，利用SPSS 統(tǒng)計(jì)分析時(shí)域與頻域之間的參量關(guān)系以及上肢肌肉疲勞間的關(guān)系，結(jié)果與實(shí)際情況更貼合。

以上研究主要是對(duì)上肢疲勞時(shí)域或頻域指標(biāo)變化規(guī)律的探索，缺乏對(duì)各疲勞指標(biāo)的評(píng)價(jià)效果進(jìn)行分析。為此，研究了肘關(guān)節(jié)肌肉疲勞評(píng)價(jià)方法的評(píng)價(jià)效果，記錄測(cè)試對(duì)象在不同負(fù)載下肘關(guān)節(jié)做屈曲運(yùn)動(dòng)的sEMG 信號(hào)，計(jì)算MNF、SMR、WIRMI1551、fApEn與RQA%DET的指標(biāo)，并進(jìn)一步比較5種方法對(duì)肘關(guān)節(jié)肌肉疲勞的評(píng)價(jià)效果。

從而獲得更有效的肘關(guān)節(jié)肌肉疲勞評(píng)價(jià)方法。將從肌肉疲勞評(píng)價(jià)方法的抗干擾性與區(qū)分疲勞程度的能力兩個(gè)方面對(duì)5種疲勞評(píng)價(jià)方法的效果進(jìn)行比較，分析所得結(jié)果，為肘關(guān)節(jié)肌肉疲勞提供一個(gè)更加準(zhǔn)確的人因工程定量評(píng)價(jià)方法。

2 研究方法

2.1 肌肉疲勞評(píng)價(jià)方法

2.1.1 平均頻率

研究表明，隨著肌肉疲勞的累積，sEMG 信號(hào)功率譜的主頻率由高轉(zhuǎn)向低，平均功率值逐漸變小，并且疲勞程度越大，下降趨勢(shì)越強(qiáng)?；趯?duì)功率譜密度（PSD）的估計(jì)，計(jì)算了譜疲勞評(píng)價(jià)方法指標(biāo)。PSD 是利用Welch 方法估算的（其中k=15），并使用Hamming窗函數(shù)進(jìn)行50%的重疊[12]?；赑SD估計(jì)的MNF：

式中：f—頻率；fs—采樣率。

2.1.2 頻譜矩

作為平均頻率的替代，頻譜矩（SMR）方法被Dimitrov等人描述為新的光譜指數(shù)，它是對(duì)H/L比率方法的改進(jìn)。文獻(xiàn)[13]用平滑函數(shù)加權(quán)PSD估計(jì)，他們建議使用頻率的倒數(shù)作為低頻帶的加權(quán)函數(shù)，高頻帶的加權(quán)函數(shù)是頻率平方，立方，四次或五次方?；陔A數(shù)k的譜矩Mk：

頻譜矩為：

式中：f1=5Hz和f2=500Hz—積分的上下限，階數(shù)k=5。

2.1.3 基于離散小波的方法

量化頻率壓縮的其他增強(qiáng)疲勞評(píng)價(jià)方法是Gonzalez?Izal等人報(bào)道的基于小波的方法[14]。它們是基于固定小波變換（SWT）的頻譜估計(jì)而不是PSD 估計(jì)。WIRM1551 比值是從Gonzalez?Izal等人建議的小波方法的變體中選擇的。它被定義為：

式中：D1—使用sym5小波的SWT的第一尺度的功率譜；D5—SWT的第五尺度的功率譜。

正如建議的那樣，f1=8Hz和f2=500Hz為積分的上下限。

2.1.4 模糊近似熵

熵估計(jì)是用于疲勞評(píng)價(jià)的另一種非線性方法。樣本熵（Sam?pEn）是一種計(jì)算時(shí)間序列熵的方法，為條件概率的負(fù)自然對(duì)數(shù)，該條件概率是從一個(gè)紀(jì)元開始的兩個(gè)序列在m維向量相似，而在m+1 維向量仍然相似，當(dāng)向量的長(zhǎng)度m=2，閾值r=0.3 為最佳搭配。為了疲勞檢測(cè)的目的，文獻(xiàn)[15]報(bào)告了一種改進(jìn)的算法，稱為模糊近似熵（fApEn）。與樣本熵相反，如果兩個(gè)序列相差任意恒定偏移量d，則模糊近似熵會(huì)將它們視為相似。以m=2，r=0.6為最佳組合計(jì)算模糊函數(shù)來獲得模糊近似熵的近似度。模糊函數(shù)為：

2.1.5 遞歸量化分析

遞歸量化分析（RQA）是一種試圖量化分析信號(hào)中重復(fù)模式的非線性疲勞評(píng)價(jià)方法[6]。首先構(gòu)造距離矩陣，從時(shí)間信號(hào)的角度導(dǎo)出長(zhǎng)度為D（嵌入維）的向量，其中時(shí)間延遲間隔λ=5，嵌入尺寸D=15。然后矩陣通過其平均值進(jìn)行歸一化處理，研究表明，閾值r=0.75時(shí)矩陣會(huì)產(chǎn)生僅包含零和一的二進(jìn)制遞歸映射。%DET是用于評(píng)估復(fù)發(fā)圖的，其值為形成一定鏈長(zhǎng)的對(duì)角線與矩陣中所有對(duì)角線的數(shù)量之比。研究發(fā)現(xiàn)，在疲勞任務(wù)下，RQA%DET值的上升與運(yùn)動(dòng)單元的更長(zhǎng)時(shí)標(biāo)同步非常吻合。

2.2 疲勞評(píng)價(jià)方法的定量化分析

由于疲勞評(píng)價(jià)方法難以基于測(cè)試對(duì)象的主觀感受來校準(zhǔn)，因此本研究使用定量化分析方法比較疲勞評(píng)價(jià)方法的差異性。定量化分析方法主要包括同一疲勞條件下評(píng)價(jià)方法的抗干擾性及不同疲勞條件下區(qū)分疲勞程度的能力。相關(guān)研究表明，肌肉疲勞評(píng)估指數(shù)呈現(xiàn)具有恒定斜率的線[16]，本研究利用最小二乘法對(duì)五種疲勞評(píng)估方法進(jìn)行了線性回歸分析，并求解出每個(gè)方法的確定系數(shù)R2與回歸方程的斜率k。

2.2.1 抗干擾性分析

疲勞評(píng)價(jià)方法預(yù)期在疲勞情況下產(chǎn)生平滑變化的疲勞指數(shù)。偏離這個(gè)平滑過程被認(rèn)為是干擾。疲勞指數(shù)的干擾可以通過確定系數(shù)R2來量化[17]。R2表示具有恒定斜率的線反映信號(hào)數(shù)據(jù)的相關(guān)性，R2的數(shù)值范圍從0到1。R2的值為0意味著信號(hào)中沒有線性趨勢(shì)。R2的值為1則意味著信號(hào)恰好是具有非零斜率且沒有額外干擾的直線。因此R2的值越大表明疲勞評(píng)估算法的抗干擾性越強(qiáng)。

2.2.2 區(qū)分疲勞程度能力分析

人體肘關(guān)節(jié)肌肉在不同負(fù)載下會(huì)產(chǎn)生不同的肌肉疲勞程度。柯爾莫哥洛夫?斯米洛夫檢驗(yàn)（Kolmogorov?Smirnov test，即K-S檢驗(yàn)）是一種基于概率統(tǒng)計(jì)理論的方法，它能夠計(jì)算兩個(gè)數(shù)據(jù)經(jīng)驗(yàn)分布函數(shù)間的最大垂直距離Lmax，Lmax越大說明兩個(gè)數(shù)據(jù)相似性越小。

本研究利用K-S檢驗(yàn)分析五種評(píng)價(jià)方法區(qū)分疲勞程度的能力[18]。即對(duì)五種評(píng)價(jià)方法在不同負(fù)載下的線性回歸方程斜率k進(jìn)行K-S檢驗(yàn)得到的Lmax值，Lmax越大表明疲勞評(píng)估方法對(duì)不同負(fù)載下區(qū)分疲勞程度的能力越強(qiáng)。

3 肌電信號(hào)采集與處理

3.1 測(cè)試對(duì)象與試驗(yàn)準(zhǔn)備

3.1.1 測(cè)試對(duì)象

為了采集數(shù)據(jù)分析比較不同疲勞評(píng)價(jià)方法，本研究采集了12名測(cè)試對(duì)象（6名女性和6名男性）肱二頭肌的表面肌電信號(hào)。測(cè)試對(duì)象的年齡在（19～26）歲之間，無皮膚過敏史且肱二頭肌皮膚表面無創(chuàng)傷和疤痕，上肢運(yùn)動(dòng)能力處于正常狀態(tài)。在試驗(yàn)開始之前，所有測(cè)試對(duì)象都書面同意志愿參加本次試驗(yàn)。

3.1.2 試驗(yàn)準(zhǔn)備

為了施加恒定的力，要求每個(gè)測(cè)試對(duì)象以同一固定姿勢(shì)托舉啞鈴。實(shí)驗(yàn)過程中要求測(cè)試對(duì)象保持直立坐姿，肘部處于放松狀態(tài)后抓舉啞鈴進(jìn)行實(shí)驗(yàn)。有關(guān)動(dòng)作設(shè)置，如圖1所示。

圖1 實(shí)驗(yàn)設(shè)定Fig.1 Setup of the Experiment

數(shù)據(jù)采集過程中，首先引導(dǎo)測(cè)試對(duì)象以最大可實(shí)現(xiàn)的用力在短時(shí)間內(nèi)抓舉啞鈴以單獨(dú)測(cè)取每個(gè)測(cè)試對(duì)象的最大自主收縮力（MVC）。然后在隨后的三次試運(yùn)行期間，根據(jù)每個(gè)測(cè)試對(duì)象的MVC調(diào)整他們要托舉的重量。在三次試驗(yàn)中，要求測(cè)試對(duì)象分別托舉他們MVC 的20%，40%和60%的重量。在20%和40%MVC 負(fù)荷的情況下，所有測(cè)試對(duì)象持續(xù)3min 的指定時(shí)間，在60%MVC負(fù)荷水平下持續(xù)135s達(dá)到疲勞。每次試驗(yàn)后，測(cè)試對(duì)象休息30min。

3.2 試驗(yàn)設(shè)備與數(shù)據(jù)采集

3.2.1 試驗(yàn)設(shè)備

葡萄牙Plux公司生產(chǎn)ErgoPlux的8通道無線表面肌電測(cè)量系統(tǒng)、上海勵(lì)圖醫(yī)療器材有限公司生產(chǎn)的LT?7型理療電極片。

3.2.2 數(shù)據(jù)采集

測(cè)試肌肉位置，如圖2所示。先用75%的醫(yī)用酒精擦拭測(cè)試對(duì)象右臂肱二頭肌的皮膚表面，清理皮膚表面污漬與角質(zhì)，再使用干燥的棉球擦凈皮膚表面殘留的酒精后將八組電極片以高度密集的形式放置在肱二頭肌上方。即電極以彼此相鄰的兩列連接，且彼此之間沒有任何間隔。根據(jù)其直徑，相鄰放置的電極的中心間隔為26mm。測(cè)試電極沿著所測(cè)試肌肉的肌纖維方向，貼在肱二頭肌最隆處，接地電極貼在肘關(guān)節(jié)處。試驗(yàn)過程中，每名測(cè)試對(duì)象依次托舉他們最大負(fù)載（MVC）的20%、40%與60%重量，利用8 通道的表面肌電測(cè)量系統(tǒng)同步記錄測(cè)試對(duì)象的sEMG信號(hào)。

圖2 高密度電極片的貼點(diǎn)位置圖Fig.2 Location of the High?Density Electrode Sheet in the Biceps

3.3 數(shù)據(jù)處理

經(jīng)多次試驗(yàn)，測(cè)試對(duì)象在測(cè)試結(jié)束后均在不同負(fù)載下感受到不同程度的肌肉酸痛僵硬，即肌肉處于疲勞狀態(tài)。為保證所分析的sEMG信號(hào)數(shù)據(jù)均在疲勞狀態(tài)下，研究采用測(cè)試過程中最后的1min進(jìn)行數(shù)據(jù)分析。首先同步采集8個(gè)通道的sEMG 會(huì)產(chǎn)生基線偏移現(xiàn)象，利用sgolay濾波器消除趨勢(shì)化，然后sEMG信號(hào)采集過程中會(huì)受到50Hz工頻及其他信號(hào)干擾，利用50Hz陷波與巴特沃斯帶通濾波器對(duì)sEMG信號(hào)去噪[18]，最后將處理后的8通道信號(hào)取平均值。處理后的示例信號(hào)，如圖3所示。

圖3 處理后的示例信號(hào)Fig.3 Exemplary Signals After Processing

根據(jù)五種疲勞評(píng)價(jià)方法的定義，使用Matlab R2014b分別編寫五種疲勞評(píng)價(jià)方法的算法，計(jì)算出每個(gè)測(cè)試對(duì)象在不同負(fù)載下各疲勞評(píng)價(jià)方法的指標(biāo)值，然后對(duì)各疲勞評(píng)價(jià)方法進(jìn)行定量化分析。定量化分析分為兩個(gè)部分，一是利用最小二乘法對(duì)各方法的指標(biāo)進(jìn)行線性回歸分析（圖4 展示了1 名測(cè)試對(duì)象在不同負(fù)載下最后一分鐘的MNF最小二乘法線性回歸分析結(jié)果，其他4 種方法與MNF 做相同的處理），比較線性回歸方程的確定系數(shù)R2的均值，利用SPSS25對(duì)R2進(jìn)行單因素方差分析，進(jìn)一步評(píng)價(jià)疲勞評(píng)價(jià)方法的抗干擾性；二是對(duì)回歸方程的斜率k進(jìn)行K-S檢驗(yàn)得到不同負(fù)載下的Lmax值，并對(duì)Lmax進(jìn)行單因素方差分析，評(píng)價(jià)5種疲勞評(píng)價(jià)方法區(qū)分疲勞程度的能力。圖4中=0.792，正如預(yù)期的那樣，隨著負(fù)載水平的增加，R2值會(huì)增加，因?yàn)樨?fù)載水平越高，MNF的線性趨勢(shì)越強(qiáng)。

圖4 MNF的線性回歸分析結(jié)果Fig.4 Linear Regression Analysis Results of MNF

4 結(jié)果

4.1 抗干擾評(píng)價(jià)結(jié)果

為量化不同疲勞算法的抗干擾性，本研究計(jì)算了12個(gè)測(cè)試對(duì)象在不同負(fù)載下各疲勞評(píng)價(jià)方法線性回歸方程的確定系數(shù)R2平均值與標(biāo)準(zhǔn)差，并對(duì)R2進(jìn)行單因素方差分析結(jié)果，如圖5所示。

圖5（a）顯示在20%MVC情況下，SMR的R2平均值為5種疲勞評(píng)價(jià)方法的最大值，最大平均值=0.237，且與MNF（p=8.28e?04）和RQA%DET（p=5.45e?05）方法存在顯著性差異。表明SMR在20%MVC下具有最好的抗干擾性。WIRM1551具有較大的R2平均值=0.223，說明WIRM1551在20%MVC下具有較好的抗干擾能力。RQA%DET的R2平均值最小，表明RQA%DET的抗干擾能力最差。

圖5（b）顯示在40% MVC 情況下，SMR的R2平均值仍然最大=0.328，且與MNF（p=1.38e?05）和RQA%DET（p=3.69e?06）在評(píng)價(jià)肘關(guān)節(jié)肌肉疲勞的抗干擾性的差異具有統(tǒng)計(jì)學(xué)意義。表明SMR 在40%MVC 下具有最強(qiáng)的抗干擾能力。MNF、WIRM1551與fApEn的R2平均值相差很小，分別為=0.289，=0.288，說明這3種方法在40%MVC下具有相同的抗干擾能力。而RQA%DET的抗干擾性最差。

圖5（c）顯示MNF、WIRM1551與fApEn在60%MVC具有較大的R2平均值，但SMR仍然為最大值=0.528，且與RQA%DET（p=3.08e?07）的差異具有統(tǒng)計(jì)學(xué)意義。說明MNF、WIRM1551與fApEn在60%MVC下具有較好的抗干擾性，但SMR的抗干擾能力是最強(qiáng)的。RQA%DET依舊是5種方法中抗干擾性最差的。

圖5 五種肌肉疲勞評(píng)價(jià)方法在不同負(fù)載下的R2Fig.5 Five Muscle Fatigue Detection Methods for R2 Under Different Loads

綜上所述，SMR在不同負(fù)載下的R2平均值均為最大值，擁有最好的抗干擾性。WIRM1551 在不同負(fù)載下的R2平均值僅比SMR 小一些，同樣具有較好的抗干擾性。MNF 與fApEn 僅在大負(fù)載下?lián)碛休^大的R2平均值，但是在小負(fù)載，低疲勞情況下的R2平均值相比SMR還是有較大差距的，因此，MNF與fApEn僅在大負(fù)載下?lián)碛休^強(qiáng)的抗干擾能力。而RQA%DET在不同負(fù)載下的抗干擾能力是最弱的。

4.2 疲勞程度區(qū)分能力評(píng)價(jià)結(jié)果

為量化分析5種疲勞評(píng)價(jià)方法區(qū)分疲勞程度的能力，對(duì)線性回歸方程的斜率k進(jìn)行K-S檢驗(yàn)，得出不同負(fù)載下肘關(guān)節(jié)肌肉疲勞指標(biāo)經(jīng)驗(yàn)概率分布函數(shù)的最大垂直距離Lmax，并對(duì)Lmax進(jìn)行單因素方差分析。全局結(jié)果，如圖6所示。局部Lmax量化分析結(jié)果，如圖7所示。

圖6 全局Lmax量化分析Fig.6 Gl obalLmaxQuantitativeAnalysis

不同負(fù)載下全局Lmax量化分析結(jié)果，如圖6所示。結(jié)果顯示在疲勞程度較大（40%vs60%）的情況下，MNF、WIRM1551與fAp?En的Lmax平均值較大，SMR的Lmax平均值為最大值；在疲勞差異大（20%vs60%）的情況下，SMR與WIRM1551的Lmax平均值為最大值；在疲勞程度小（20%vs40%）的情況下，SMR的Lmax平均值呈現(xiàn)明顯的尖點(diǎn)。而RQA%DET的Lmax平均值在三種情況均為最小值。全局Lmax量化分析結(jié)果表明在5種方法中SMR擁有最強(qiáng)的區(qū)分肘關(guān)節(jié)肌肉疲勞程度的能力，而RQA%DET是最弱的。

20%與60%MVC的Lmax量化分析結(jié)果，如圖7（a）所示。結(jié)果顯示SMR 與WIRM1551的Lmax為最大值，均為0.883，表明SMR與WIRM1551 在區(qū)分20%與60%MVC 的疲勞程度的能力最出色。而RQA%DET 是5 種方法中評(píng)價(jià)疲勞程度效果最差的。SMR與RQA%DET的差異具有統(tǒng)計(jì)學(xué)意義（p=1.09e?02）。表明SMR與WIRM1551在疲勞程度差異較大時(shí)有最強(qiáng)的區(qū)分疲勞程度能力。40%與60%MVC的Lmax量化分析結(jié)果，如圖7（b）所示。結(jié)果表明MNF、SMR、WIRM1551 與fApEn 均呈現(xiàn)較大的Lmax平均值，說明4 種方法在評(píng)價(jià)大負(fù)載引起的疲勞程度均有很出色的能力。

其中SMR 是最出色的，Lmax=0.766，且與RQA%DET（p=8.93e?03）存在顯著性差異。而RQA%DET區(qū)分疲勞程度的能力最差。說明在疲勞程度較大時(shí)，MNF、SMR、WIRM1551與fApEn均具有較好的區(qū)分能力，且SMR 是最好的。20%vs40%MVC的Lmax量化分析結(jié)果，如圖7（c）所示。

結(jié)果表明SMR的Lmax平均值最大，最大值Lmax=0.622，且與RQA%DET（p=1.6e?03）存在顯著性差異。RQA%DET的Lmax平均值仍是最小的。說明在疲勞程度較小的情況下，SMR仍然擁有最強(qiáng)的區(qū)分肌肉疲勞程度的能力。

綜上所述，SMR相較于另外4種方法在不同負(fù)載下都具有最大的Lmax平均值，說明SMR在區(qū)分不同負(fù)載引起的肘關(guān)節(jié)肌肉疲勞程度的能力最出色。WIRM1551 雖沒有SMR 那么出色，但在不同負(fù)載下均具有較好的區(qū)分疲勞程度的能力。MNF與fApEn僅在疲勞程度大或疲勞差異較大時(shí)具有較好的區(qū)分能力。而RQA%DET 在不同負(fù)載下具有最小的Lmax平均值，說明RQA%DET在不同負(fù)載下引起的肌肉疲勞程度區(qū)分性最差。

5 結(jié)論

通過采集肱二頭肌的表面肌電信號(hào)計(jì)算了5種常見肌肉疲勞評(píng)價(jià)算法，從評(píng)價(jià)算法的抗干擾性與區(qū)分疲勞程度的能力兩個(gè)方面量化其疲勞評(píng)價(jià)性能，以R2大小體現(xiàn)抗干擾能力的強(qiáng)弱，以Lmax展現(xiàn)區(qū)分疲勞程度的能力。實(shí)驗(yàn)結(jié)果顯示，頻譜距SMR同時(shí)擁有最大的R2和Lmax值，即頻譜距SMR在評(píng)價(jià)肘關(guān)節(jié)肌肉疲勞時(shí)具有最好抗干擾和區(qū)分疲勞程度的能力。故頻譜距SMR能夠更加及時(shí)有效的檢測(cè)出肘關(guān)節(jié)肌肉疲勞狀態(tài)。因此在今后研究肘關(guān)節(jié)屈曲運(yùn)動(dòng)的肌肉疲勞時(shí)，建議使用頻譜矩（SMR）進(jìn)行實(shí)時(shí)監(jiān)測(cè)。