陳建國(guó)， 楊 徐， 李玉榕

(福州大學(xué)電氣工程與自動(dòng)化學(xué)院， 福建省醫(yī)療器械和醫(yī)藥技術(shù)重點(diǎn)實(shí)驗(yàn)室， 福建 福州 350116)

0 引言

近年來(lái)， 腦卒中(腦中風(fēng))或外傷引起的足下垂患者數(shù)量逐年增多. 足下垂(foot drop)是指患者踝關(guān)節(jié)無(wú)法產(chǎn)生背屈， 走路時(shí)無(wú)法使腳背彎起的癥狀. 足下垂患者在行走時(shí)腳趾拖地， 與正常人的下肢差異很大， 不僅行走不便， 而且容易摔倒. 如果在康復(fù)訓(xùn)練時(shí)摔倒更易導(dǎo)致二次中風(fēng)， 后果將更加嚴(yán)重. 因此, 如何利用新技術(shù)、 采用新方法使足下垂康復(fù)治療儀體積更小、 功能更健全、 使用更方便， 幫助患者恢復(fù)正常的行走能力， 減少其對(duì)家人和醫(yī)護(hù)人員的依賴， 對(duì)于足下垂患者本人和家庭而言具有重要的意義， 同時(shí)能產(chǎn)生良好的社會(huì)效益. 功能性電刺激是目前應(yīng)用最為廣泛的足下垂康復(fù)治療方法之一， 通過(guò)在患者行走的擺動(dòng)階段以低頻脈沖電流刺激脛骨前肌， 促進(jìn)脛骨前肌主動(dòng)收縮， 使足部產(chǎn)生背屈運(yùn)動(dòng)， 從而治療足下垂[1]. 相比其他康復(fù)治療方法， 功能性電刺激可以使患者行走時(shí)更加自信、 步態(tài)更穩(wěn)健且不易疲勞. 步態(tài)即為人體走路時(shí)的姿態(tài)， 包含了許多運(yùn)動(dòng)信息， 步態(tài)信息的分段在功能性電刺激中起著非常重要的作用. 目前， 步態(tài)分段方法存在識(shí)別率低、 實(shí)時(shí)性差、 過(guò)于復(fù)雜等缺點(diǎn). 如何利用采集到的步態(tài)數(shù)據(jù)進(jìn)行步態(tài)分段， 實(shí)時(shí)識(shí)別出步態(tài)的擺動(dòng)相和站立相， 是應(yīng)用功能性電刺激治療足下垂的關(guān)鍵因素.

荷蘭生物醫(yī)學(xué)技術(shù)研究所的Monaghan C C等[2]采用閾值分類法， 將采集到的角速度進(jìn)行積分得到角度， 通過(guò)設(shè)定閾值來(lái)進(jìn)行步態(tài)分段. 該方法的積分計(jì)算會(huì)帶來(lái)累積誤差， 容易導(dǎo)致誤判. 為解決積分帶來(lái)的誤差問(wèn)題， 英國(guó)牛津大學(xué)的Morri[3]采用一個(gè)步態(tài)周期的起始和結(jié)束時(shí)刻的角度相等且此時(shí)腳的速度為零的方法來(lái)進(jìn)行步態(tài)分段， 但其有效性和實(shí)時(shí)性需進(jìn)一步驗(yàn)證. 臺(tái)灣大學(xué)的SHIHY Y等[4]利用神經(jīng)網(wǎng)絡(luò)算法將傳感器測(cè)得的角度進(jìn)行訓(xùn)練， 識(shí)別出步態(tài)周期的各個(gè)階段， 并采用PID控制功能性電刺激的輸出. 此方法雖然能夠進(jìn)行實(shí)時(shí)步態(tài)分段， 但每個(gè)患者使用前都必須先用傳感器采集步態(tài)信息進(jìn)行訓(xùn)練， 所以使用起來(lái)比較復(fù)雜. 馬來(lái)西亞莫納什大學(xué)的Gouwand[5]利用角速度信號(hào)對(duì)步態(tài)周期各階段進(jìn)行識(shí)別， 大連理工大學(xué)的韓寧[6]將傳感器采集到的角速度用基于知識(shí)的步態(tài)分段算法進(jìn)行步態(tài)分段. 雖然兩者能夠比較準(zhǔn)確地識(shí)別出步態(tài)各個(gè)階段， 誤差小， 但是此算法中判斷前一個(gè)步態(tài)階段需要依賴于后一個(gè)步態(tài)階段的數(shù)據(jù)， 實(shí)時(shí)性差. 大連理工大學(xué)的Wang和Qiu等[7]將加速度和角速度信號(hào)采用滑動(dòng)方差法， 即分別對(duì)加速度和角速度求方差作為閾值， 對(duì)一個(gè)周期內(nèi)的步態(tài)各階段進(jìn)行識(shí)別， 分別是腳跟離地、 擺動(dòng)相、 腳跟著地以及完全站立. 該算法需要采集所有數(shù)據(jù)后才能求方差， 不能實(shí)現(xiàn)一邊采集數(shù)據(jù)一邊進(jìn)行步態(tài)識(shí)別， 因此不能應(yīng)用到實(shí)時(shí)功能性電刺激治療中. 丹麥奧爾堡衛(wèi)生科技部的Jan Rueterbories和Erika G Spaich[8]采用步態(tài)相位檢測(cè)方法， 雖然步態(tài)分段時(shí)有較好的實(shí)時(shí)性， 但是系統(tǒng)容易受到患者跛行步態(tài)影響， 導(dǎo)致較高的錯(cuò)誤率.

本研究提出的基于角速度信號(hào)的步態(tài)分段算法不需要積分運(yùn)算， 避免了累積誤差， 算法簡(jiǎn)單且準(zhǔn)確性較高， 最重要的是該算法實(shí)時(shí)性好， 能夠?yàn)檠芯看┐魇焦δ苄噪姶碳は到y(tǒng)提供基礎(chǔ)與保證.

1 系統(tǒng)總體設(shè)計(jì)方案

本研究為基于閾值步態(tài)分段算法的功能性電刺激系統(tǒng)， 系統(tǒng)總體方案如圖1所示. 將穿戴下肢的MPU6050傳感器采集到的數(shù)據(jù)通過(guò)I2C傳輸?shù)紺C2541單片機(jī)中， 然后采用步態(tài)分段算法進(jìn)行處理， 并控制功能性電刺激的輸出. 當(dāng)檢測(cè)到腳跟離地時(shí)， 控制功能性電刺激模塊對(duì)足下垂患者的脛前肌施加電刺激信號(hào)， 促進(jìn)踝關(guān)節(jié)產(chǎn)生背屈； 當(dāng)檢測(cè)到腳跟著地時(shí)， 停止電刺激輸出信號(hào)， 使患者能夠正常行走， 從而治療足下垂患者.

圖1 系統(tǒng)總體方案圖Fig.1 Overall scheme of the system

2 基于閾值的角速度步態(tài)分段算法

2.1 步態(tài)周期分析

圖2 正常人步態(tài)周期的時(shí)間分解Fig.2 Time breakdown of normal gait cycle

人的步行過(guò)程是一個(gè)有規(guī)律的模式， 對(duì)于規(guī)律的運(yùn)動(dòng)而言可以研究其周期性. 一個(gè)步態(tài)周期是指在行走過(guò)程中， 從一只腳完全站立在地面上到這只腳再次完全站立在地面上所經(jīng)歷的時(shí)間. 一個(gè)完整的步態(tài)周期可以分為三個(gè)階段： 站立階段、 推動(dòng)階段、 擺動(dòng)階段[9]， 如圖2所示. 文[10]研究表明， 正常人的步態(tài)周期約1.3 s， 其中站立階段的時(shí)間約620 ms， 從腳跟離地到腳尖離地有一個(gè)力推動(dòng)過(guò)程， 約200 ms， 擺動(dòng)階段的時(shí)間約470 ms.

腳跟離地(heel off， HO)是推動(dòng)階段的起點(diǎn)， 也是站立階段的終點(diǎn). 腳尖離地(toe off， TO)是擺動(dòng)階段的起點(diǎn)， 也是推動(dòng)階段的終點(diǎn). 腳跟著地(heel strike， HS)是站立階段的起始， 也即擺動(dòng)階段終點(diǎn)的標(biāo)志. 功能性電刺激通過(guò)在行走步態(tài)的擺動(dòng)階段刺激脛骨前肌從而治療足下垂， 也就是需要在HO時(shí)刻開(kāi)啟電刺激， 在HS時(shí)刻停止電刺激. 因此， 為了正確控制功能性電刺激輸出， 需要研究步態(tài)分段算法， 實(shí)時(shí)確定HO時(shí)刻和HS時(shí)刻.

在功能性電刺激時(shí)， 刺激的開(kāi)啟時(shí)刻尤為重要， 因?yàn)槿绻阆麓够颊吣_尖已經(jīng)離地， 而此時(shí)如果脛骨前肌仍舊無(wú)力的話， 會(huì)導(dǎo)致患者重心不穩(wěn)而跌倒. 相關(guān)實(shí)驗(yàn)研究表明[11-12]， 肌肉在外部電刺激作用下從零肌肉力輸出到最大肌肉力輸出需要一定時(shí)間， 包括電信號(hào)在肌肉的傳輸延遲和肌肉本身的反應(yīng)時(shí)間， 大約50～150 ms， 個(gè)體有差異， 如圖2中陰影部分所示. 因此， 實(shí)際應(yīng)用中， 應(yīng)當(dāng)選擇腳尖離地TO時(shí)刻的前50～150 ms的時(shí)刻作為電刺激的起點(diǎn). 考慮到尋找該特征點(diǎn)時(shí)， 算法本身會(huì)存在時(shí)延， 因此， 在選擇功能性電刺激起點(diǎn)時(shí)， 必須再加上算法的延時(shí). 可以選擇推動(dòng)階段的起點(diǎn)， 即腳跟離地HO時(shí)刻作為電刺激起點(diǎn). 而功能性電刺激的結(jié)束時(shí)刻的準(zhǔn)確性相比較而言可以放松些， 因?yàn)榧词鼓_跟著地后， 電刺激仍然存在， 患者只是感覺(jué)上有一點(diǎn)不舒服， 不至于發(fā)生危險(xiǎn).

2.2 步態(tài)數(shù)據(jù)分析

從文[13]可以得出， 平地自由行走模式下將傳感器穿戴于小腿和踝關(guān)節(jié)， 其矢狀面(通過(guò)人體沿垂軸和縱軸的平面及與其平行的所有平面都稱為矢狀面)角速度幅值大、 變化比較均勻、 規(guī)律性明顯、 易于尋找特征點(diǎn). 而且穿戴于小腿的不同位置， 角速度數(shù)據(jù)規(guī)律性一致， 故小腿以及踝關(guān)節(jié)可以認(rèn)為是傳感器的理想穿戴位置， 而且小腿穿戴位置的變化不至于影響結(jié)果， 這樣便于患者使用時(shí)可選擇在小腿上某個(gè)舒適的位置穿戴.

在本次步態(tài)數(shù)據(jù)采集實(shí)驗(yàn)的過(guò)程中， 選擇8名志愿者參與實(shí)驗(yàn)， 其中男女各四名. 為驗(yàn)證所設(shè)計(jì)的系統(tǒng)可以利用步態(tài)分段算法的識(shí)別結(jié)果實(shí)現(xiàn)功能性電刺激的正確輸出， 選擇了健康的志愿者， 后續(xù)階段的研究將針對(duì)不同類別的足下垂患者進(jìn)行實(shí)驗(yàn). 實(shí)驗(yàn)過(guò)程中志愿者們均穿著跑步鞋， 并在指定位置穿戴好步態(tài)測(cè)量的傳感器， 以在直線自由行走大約10 m和在跑步機(jī)上以不同速度行走兩種情況來(lái)進(jìn)行步態(tài)數(shù)據(jù)采集. 由于此系統(tǒng)主要用于足下垂患者康復(fù)訓(xùn)練， 所以要求志愿者們?cè)趯?shí)驗(yàn)中步行速度盡量放慢.

采集數(shù)據(jù)時(shí)， 為研究傳感器不同穿戴位置對(duì)步態(tài)分段結(jié)果的影響， 需將采集數(shù)據(jù)的傳感器穿戴在下肢的各個(gè)部位. 因此， 研究前期將系統(tǒng)穿戴在人體下肢， 一邊采集數(shù)據(jù)一邊將其儲(chǔ)存在單片機(jī)中. 步態(tài)實(shí)驗(yàn)結(jié)束后再將數(shù)據(jù)通過(guò)串口調(diào)試助手傳輸?shù)接?jì)算機(jī)， 用Matlab軟件進(jìn)行處理， 尋找最優(yōu)傳感器穿戴位置及進(jìn)行算法研究. 最終再將算法在單片機(jī)系統(tǒng)中實(shí)現(xiàn)， 以實(shí)現(xiàn)實(shí)時(shí)控制.

通過(guò)分析人體的行走特點(diǎn)可以發(fā)現(xiàn)， 人行走時(shí)的角速度信號(hào)規(guī)律比加速度信號(hào)規(guī)律更明顯， 利用角速度劃分三個(gè)階段的步態(tài)周期更簡(jiǎn)便[14]. 而且實(shí)驗(yàn)所利用的陀螺儀， 又稱之為角速度傳感器， 基本能夠消除重力對(duì)運(yùn)動(dòng)的影響， 并且噪聲很小[15]. 因此， 采用角速度信號(hào)作為步態(tài)分析的特征數(shù)據(jù).

2.3 步態(tài)分段算法的實(shí)現(xiàn)

圖3 一個(gè)步態(tài)周期矢狀面角速度變化曲線 Fig.3 Angular velocity of sagittal plane in one gait cycle

從圖2可以發(fā)現(xiàn)， 一個(gè)完整的步態(tài)周期可以分為三個(gè)階段： 站立階段、 推動(dòng)階段和擺動(dòng)階段. 推動(dòng)階段前的時(shí)刻為腳跟離地HO時(shí)刻， 站立階段前的時(shí)刻為腳跟著地HS時(shí)刻， 通過(guò)本章2.1小節(jié)的分析， 可以將腳跟離地HO和腳跟著地HS時(shí)刻分別作為功能性電刺激開(kāi)啟和關(guān)斷的時(shí)刻. 因此， 需要找到合適的算法確定出兩個(gè)關(guān)鍵的時(shí)刻： 腳跟離地HO和腳跟著地時(shí)刻HS. 實(shí)驗(yàn)采集的一個(gè)標(biāo)準(zhǔn)步態(tài)周期的小腿上矢狀面角速度如圖3所示, 分析在步態(tài)周期中關(guān)鍵特征點(diǎn).

1) 腳跟離地時(shí)刻的特征. 腳跟離地時(shí)刻為站立階段的終點(diǎn)， 理論上站立階段矢狀面的角速度為0， 當(dāng)矢狀面的角速度從0變成負(fù)時(shí)為腳跟離地時(shí)刻(按左手守則， 以逆時(shí)針?lè)较驗(yàn)檎?. 由于測(cè)量時(shí)陀螺儀靜止漂移等因素， 站立階段的角速度往往不為零， 所以認(rèn)為腳跟著地以后， 在一定時(shí)間之內(nèi)矢狀面角速度連續(xù)變小的點(diǎn)就是腳跟離地時(shí)刻.

2) 腳尖離地時(shí)刻的特征. 當(dāng)腳跟離地后矢狀面角速度達(dá)到負(fù)方向最大時(shí)(波谷)， 為腳尖離地時(shí)刻.

3)擺動(dòng)階段中點(diǎn)的特征. 擺動(dòng)階段中點(diǎn)有非常明顯的特征， 即當(dāng)腳尖完全離開(kāi)地面進(jìn)入擺動(dòng)階段時(shí)， 矢狀面角速度經(jīng)歷一個(gè)由負(fù)到正的變化， 直至達(dá)到最大角速度(波峰)， 即為擺動(dòng)階段中點(diǎn).

4) 腳跟著地時(shí)刻的特征. 從擺動(dòng)階段轉(zhuǎn)變?yōu)槟_跟著地時(shí)， 角速度會(huì)出現(xiàn)一個(gè)負(fù)的波谷， 所以在擺動(dòng)中點(diǎn)后的第一個(gè)波谷為腳跟著地點(diǎn)， 完成一個(gè)步態(tài)周期， 并準(zhǔn)備進(jìn)入下一個(gè)步態(tài)周期階段.

基于上述所示的步態(tài)周期矢狀面角速度曲線中的關(guān)鍵特征點(diǎn)， 且為簡(jiǎn)化算法保證其在單片機(jī)中實(shí)時(shí)實(shí)現(xiàn)， 采用基于閾值的角速度信號(hào)步態(tài)分段算法進(jìn)行步態(tài)周期分段. 其步態(tài)分段過(guò)程為：

1) 首先讀取數(shù)據(jù)， 當(dāng)N=0時(shí)， 系統(tǒng)將執(zhí)行判斷擺動(dòng)階段中點(diǎn)的子程序， 對(duì)于擺動(dòng)階段中點(diǎn)， 一開(kāi)始先判斷是否大于閾值λ1且為正的最大值(波峰)， 若滿足此條件則判斷為擺動(dòng)階段中點(diǎn)， 則N=1， 否則繼續(xù)讀取下一個(gè)數(shù)據(jù).

2) 當(dāng)N=1時(shí)， 系統(tǒng)將執(zhí)行判斷腳跟著地的子程序， 對(duì)于腳跟著地， 讀取數(shù)據(jù)后， 判斷此點(diǎn)是否為擺動(dòng)階段中點(diǎn)的第一個(gè)波谷， 若滿足條件則判斷為腳跟著地， 此時(shí)N=2， 否則繼續(xù)讀取數(shù)據(jù).

3) 當(dāng)N=2時(shí)， 系統(tǒng)將執(zhí)行判斷腳跟離地的子程序， 為了更容易識(shí)別這個(gè)點(diǎn)， 本算法中將腳跟著地后角速度絕對(duì)值開(kāi)始變大且大于一定閾值， 而且前后一定時(shí)間內(nèi)角速度值都遞減的時(shí)刻作為腳跟離地點(diǎn)， 因此讀取數(shù)據(jù)后先要判斷：

系統(tǒng)根據(jù)此算法一邊進(jìn)行數(shù)據(jù)采集， 一邊進(jìn)行步態(tài)分段， 既簡(jiǎn)單又具有較好的實(shí)時(shí)性， 且在單片機(jī)中易于實(shí)現(xiàn).λ1、λ2、 Δt1和Δt2這四個(gè)閾值的具體取值將在下文中進(jìn)一步討論.

本研究所設(shè)計(jì)的足下垂功能性電刺激系統(tǒng)是通過(guò)識(shí)別腳跟離地HO和腳跟著地HS來(lái)分別控制功能性電刺激開(kāi)啟和停止， 對(duì)實(shí)時(shí)性要求比較高， 因此， 算法延時(shí)的長(zhǎng)短將對(duì)電刺激的康復(fù)治療效果有很大的影響. 從上面分析可以發(fā)現(xiàn)， 該算法識(shí)別腳跟著地HS沒(méi)有滯后， 但識(shí)別腳跟離地HO存在Δt2的滯后. 本算法中， Δt2最長(zhǎng)時(shí)間為120 ms， 加上電刺激信號(hào)在肌肉的傳輸延遲和肌肉本身的反應(yīng)時(shí)間50～150 ms， 因此， 從算法檢測(cè)到腳跟離地， 到脛前肌實(shí)際輸出最大肌力的總時(shí)間為大約170～270 ms. 由于本系統(tǒng)主要設(shè)計(jì)用于足下垂患者康復(fù)治療， 足下垂患者從腳跟離地HO到腳尖離地TO的推動(dòng)階段的時(shí)間要大于正常人的200 ms， 因此， 算法對(duì)功能性電刺激的開(kāi)啟時(shí)刻和停止時(shí)刻是恰當(dāng)?shù)模?可以滿足足下垂患者功能性電刺激的要求.

3 算法識(shí)別結(jié)果

在Matlab軟件上用步態(tài)分段算法識(shí)別步態(tài)周期的各個(gè)階段， 即對(duì)穿戴于踝關(guān)節(jié)和小腿上的傳感器在平地自由行走模式和跑步機(jī)平臺(tái)下測(cè)得的角速度信號(hào)采用步態(tài)分段算法識(shí)別三個(gè)特征點(diǎn)： 腳跟離地(HO)、 擺動(dòng)階段中點(diǎn)(MS)、 腳跟著地(HS). 算法中λ1的取值范圍是35～45 (°)·s-1，λ2的取值范圍是-30～-20 (°)·s-1， 時(shí)間閾值Δt1的取值范圍是80～90 ms， Δt2的取值范圍是110～120 ms.

由于λ1的取值僅對(duì)擺動(dòng)階段中點(diǎn)的識(shí)別產(chǎn)生影響， 擺動(dòng)階段中點(diǎn)的識(shí)別有利于輔助腳跟著地和腳跟離地這兩個(gè)時(shí)刻的識(shí)別， 因此有必要對(duì)λ1的取值進(jìn)行分析， 找到最佳閾值. 對(duì)λ1取不同值時(shí)擺動(dòng)階段中點(diǎn)的識(shí)別率進(jìn)行實(shí)驗(yàn)和統(tǒng)計(jì)， 實(shí)驗(yàn)中共采集8個(gè)志愿者在平地直線自由行走大約10 m和跑步機(jī)平臺(tái)下行走15 s的步態(tài)數(shù)據(jù)， 每個(gè)志愿者在每個(gè)穿戴位置大概完成了36個(gè)步態(tài)周期， 統(tǒng)計(jì)結(jié)果如表1， 其中在小腿位置選取了三個(gè)不同的測(cè)試點(diǎn).

表1 平地行走模式下擺動(dòng)階段中點(diǎn)的識(shí)別率

從表1統(tǒng)計(jì)結(jié)果可知，λ1取40 (°)·s-1時(shí)擺動(dòng)階段中點(diǎn)的識(shí)別率最高， 基本達(dá)到99.9%； 當(dāng)λ1低于40 (°)·s-1時(shí)識(shí)別率相對(duì)較低， 這是由于腳跟離地后存在次波峰造成的誤判； 當(dāng)λ1高于40 (°)·s-1時(shí)的識(shí)別率也相對(duì)較低， 這是由于選擇的閾值過(guò)大， 當(dāng)行走的角速度幅值較小時(shí)系統(tǒng)識(shí)別不到造成的結(jié)果. 因此， 本算法選擇λ1=40 (°)·s-1為最佳閾值.

對(duì)于λ2、 Δt1和Δt2這三個(gè)參數(shù)只對(duì)腳跟離地的識(shí)別產(chǎn)生影響， 而腳跟離地識(shí)別的準(zhǔn)確性對(duì)功能性電刺激最佳刺激時(shí)相至關(guān)重要. 為了得到最佳識(shí)別率的閾值， 通過(guò)取8組不同λ2、 Δt1和Δt2取值組合， 將采集到的數(shù)據(jù)在Matlab軟件上利用步態(tài)分段算法進(jìn)行步態(tài)識(shí)別， 對(duì)平地自由行走模式下腳跟離地的識(shí)別結(jié)果進(jìn)行統(tǒng)計(jì)， 如表2.

表2 平地行走模式下腳跟離地的識(shí)別率

從表2統(tǒng)計(jì)結(jié)果可知， 在三個(gè)參數(shù)各自對(duì)應(yīng)的取值范圍內(nèi)，λ2的取值比Δt1、 Δt2的取值對(duì)腳跟離地識(shí)別率的影響更大， 當(dāng)Δt1和Δt2值一定，λ2= -20 (°)·s-1時(shí)， 腳跟離地的識(shí)別率普遍比較低， 不會(huì)超過(guò)70%； 當(dāng)λ2=-30 (°)·s-1時(shí)， 腳跟離地的識(shí)別率提高很多， 都超過(guò)了90%； 因此λ2值取-30 (°)·s-1時(shí)為最佳閾值. 當(dāng)λ2和Δt1值一定， Δt2取120 ms時(shí)， 腳跟離地的識(shí)別率相較于Δt2取110 ms時(shí)更高， 為最佳閾值； 當(dāng)λ2和Δt2值一定， Δt2取80 ms時(shí)， 腳跟離地的識(shí)別率相較于Δt2取90 ms時(shí)更高， 為最佳閾值. 綜上分析可得， 取閾值λ2= -30 (°)·s-1、 Δt1=80 ms和Δt2=120 ms時(shí)， 對(duì)腳跟離地的識(shí)別率最高， 系統(tǒng)擁有較好的魯棒性.

為進(jìn)一步驗(yàn)證閾值取λ1=40 (°)·s-1、λ2=-30 (°)·s-1、 Δt1=80 ms和Δt2=120 ms是否在其他條件下也具有較高的識(shí)別率， 在平地自由行走模式和各種速度模式下對(duì)三個(gè)特征點(diǎn)： 腳跟離地(HO)、 擺動(dòng)階段中點(diǎn)(MS)和腳跟著地(HS)進(jìn)行識(shí)別， 對(duì)腳跟離地和腳跟著地的識(shí)別率進(jìn)行統(tǒng)計(jì)， 結(jié)果如表3和表4所示.

表3 平地自由行走模式下步態(tài)識(shí)別結(jié)果

表4各種速度模式下步態(tài)識(shí)別結(jié)果

分析表3和表4統(tǒng)計(jì)結(jié)果可見(jiàn)， 在平地自由行走模式下， 對(duì)穿戴在踝關(guān)節(jié)和小腿上測(cè)得的角速度用步態(tài)分段算法識(shí)別腳跟離地和腳跟著地的識(shí)別率分別為97%和99%以上； 在各種速度模式下， 對(duì)穿戴在踝關(guān)節(jié)和小腿上測(cè)得的角速度用步態(tài)分段算法識(shí)別腳跟離地時(shí)刻和腳跟著地時(shí)刻的識(shí)別率分別為97%和98%以上. 結(jié)果表明， 該算法閾值取λ1=40 (°)·s-1、λ2=-30 (°)·s-1、 Δt1=80 ms和Δt2=120 ms時(shí)， 在平地自由行走模式和各種速度模式下均具有較高的識(shí)別率， 尤其對(duì)腳跟著地識(shí)別率幾乎100%， 而腳跟離地識(shí)別率相對(duì)低一點(diǎn). 這是由于腳要離地時(shí)容易晃動(dòng)， 有多個(gè)次波谷， 而且每個(gè)人幅值有差異， 導(dǎo)致誤判. 跑步機(jī)上設(shè)置最小速度采集時(shí)， 由于速度太慢， 實(shí)驗(yàn)者不易適應(yīng)， 造成步態(tài)有些紊亂， 所以識(shí)別率也會(huì)稍低一點(diǎn). 但總體而言本算法的識(shí)別率可達(dá)到97%～100%， 而文[8]中所提到的利用人步行時(shí)的徑向和切向加速度數(shù)據(jù)進(jìn)行離散分析的算法， 其識(shí)別率為84.2%～100%. 由此可看出, 本算法雖然簡(jiǎn)單， 但是依舊能夠保持較好的識(shí)別率， 可以準(zhǔn)確控制功能性電刺激輸出.

另外， 算法對(duì)傳感器穿戴在小腿和踝關(guān)節(jié)上采集到的角速度信號(hào)均能很好地識(shí)別， 所以小腿和踝關(guān)節(jié)均能作為傳感器穿戴位置. 但考慮到傳感器穿戴在踝關(guān)節(jié)時(shí)穿鞋容易受到限制， 而且對(duì)足下垂患者進(jìn)行電刺激治療的最佳刺激位置在小腿上， 將傳感器穿戴在小腿能夠使刺激位置和采集位置一體化， 穿戴更加方便和美觀， 因此， 傳感器穿戴在小腿上位置最佳.

4 結(jié)語(yǔ)

針對(duì)目前存在的步態(tài)分段方法識(shí)別率低、 實(shí)時(shí)性差、 過(guò)于復(fù)雜等缺點(diǎn)， 提出了基于角速度信號(hào)的功能性電刺激系統(tǒng)的步態(tài)分段算法， 用于識(shí)別步態(tài)周期的三個(gè)特殊點(diǎn)： 腳跟著地點(diǎn)、 腳跟離地點(diǎn)、 擺動(dòng)階段的中點(diǎn). 當(dāng)檢測(cè)到腳跟離地時(shí)， 輸出電刺激信號(hào)刺激患者脛前肌肉， 當(dāng)檢測(cè)到腳跟著地時(shí)停止電刺激， 以此來(lái)治療足下垂. 該算法識(shí)別率高、 實(shí)時(shí)性好、 容易實(shí)現(xiàn)， 在足下垂治療領(lǐng)域具有一定的應(yīng)用價(jià)值.

參考文獻(xiàn)：

[1] 許林海， 韓麗雅．功能強(qiáng)化訓(xùn)練結(jié)合肌電生物反饋對(duì)急性腦卒中患者上肢功能的影響[J]．中國(guó)康復(fù)， 2015(3)： 185-188．

[2] MONAGHAN C C, RIEL W J B M V, VELTINK P H. Control of triceps surae stimulation based on shank orientation using a uniaxial gyroscope during gait[J]. Medical & Biological Engineering & Computing, 2009, 47(11): 1181-1188.

[3] ZDRAGKAS G I, AVARITSIOTIS J N. Gait analysis and automatic gait event identification using accelerometers[C]//IEEE International Conference on Bioinformatics and Bioengineering. Athens: IEEE, 2008: 1-6.

[4] CHEN Y L, SHIH Y Y, CHEN W L,etal． Clincal evaluation of the tilt sensors feedback controlled FES for hemiplegia[C]//Proceedings of the 26th IEEE Annual International Conference on the EMBS. San Francisco: IEEE, 2004: 4737-4740．

[5] GOUWANDA D, SENANAYAKE S M N A. Application of hybrid multi-resolution wavelet decomposition method in detecting human walking gait events[C]//International Conference of Soft Computing and Pattern Recognition. Malacca: IEEE, 2009: 580-585.

[6] 韓寧. 基于體感網(wǎng)的步態(tài)分析算法的研究[D]. 大連： 大連理工大學(xué), 2013.

[7] WANG Z L, QIU S, CAO Z L,etal. Quantitative assessment of dual gait analysis based on inertial sensors with body sensor network[J]. Sensor Review, 2013, 33(1): 48-56.

[8] RUETERBORIES J, SPAICH E G, ANDERSEN O K. Gait event detection for use in FES rehabilitation by radial and tangential foot accelerations[J]. Medical Engineering & Physics, 2014, 36(4): 502-508.

[9] 王桂茂, 嚴(yán)雋陶, 劉玉超, 等. 腦卒中偏癱步態(tài)的時(shí)空參數(shù)與骨盆運(yùn)動(dòng)學(xué)分析[J]. 中國(guó)康復(fù)醫(yī)學(xué)雜志, 2010, 25(12): 1148-1151.

[10] KOTIADIS D, HERMENS H J, VELTINK P H. Inertial gait phase detection for control of a drop foot stimulator inertial sensing for gait phase detection[J]. Medical Engineering & Physics, 2010, 32(4): 287-297.

[11] 譚志梅, 姜文文, 燕鐵斌. 多通道功能性電刺激及其在腦卒中偏癱側(cè)肢體康復(fù)中的應(yīng)用[J]. 中華物理醫(yī)學(xué)與康復(fù)雜志, 2011, 33(6): 464-467.

[12] 馬長(zhǎng)波． 基于電機(jī)輔助功能性電刺激腳踏車康復(fù)訓(xùn)練系統(tǒng)的研究[D]． 哈爾濱： 哈爾濱工業(yè)大學(xué)， 2010．

[13] 黃金蘭, 李玉榕, 陳建國(guó). 基于微電子機(jī)械系統(tǒng)傳感器的功能性電刺激系統(tǒng)[J]. 生物醫(yī)學(xué)工程研究, 2015, 34(4):207-211.

[14] TIEN I, GLASER S D, BAJCSY R,etal. Results of using a wireless inertial measurirlg system to quantify gait motions in control subjects[J]. IEEE Transactions on Information Technology in Biomedicine, 2010, 14(4): 904-915.

[15] 侯季健. 人體步態(tài)信息采集系統(tǒng)設(shè)計(jì)[D]. 大連： 大連理工大學(xué), 2015.

[16] GOUWANDA D, GOPALAI A A. A robust real-time gait event detection using wireless gyroscope and its application on normal and altered gaits[J]. Medical Engineering & Physics, 2015, 37(2): 219-225.