劉 磊 楊 鵬 劉作軍 宋寅卯

(鄭州輕工業(yè)大學(xué)建筑環(huán)境工程學(xué)院 鄭州450002)

0 引言

下肢外骨骼助力機(jī)器人是一種可為人體下肢行走提供關(guān)節(jié)助力的裝置,其特點(diǎn)在于可以識(shí)別各種路況與步速[1-3]。由于外骨骼機(jī)器人實(shí)際運(yùn)動(dòng)意圖與期望運(yùn)動(dòng)意圖存在一定差異(運(yùn)動(dòng)不相容),穿戴者容易摔倒。外骨骼機(jī)器人的關(guān)節(jié)機(jī)構(gòu)與控制方法是目前研究的熱點(diǎn)問(wèn)題,但對(duì)外骨骼機(jī)器人的運(yùn)動(dòng)相容性識(shí)別研究不多。當(dāng)外骨骼機(jī)器人的運(yùn)動(dòng)軌跡與人體期望運(yùn)動(dòng)軌跡不一致時(shí),為實(shí)現(xiàn)外骨骼助行機(jī)器人可自適應(yīng)調(diào)整,本文對(duì)步幅過(guò)大、步幅過(guò)小及步幅相容3 種人機(jī)運(yùn)動(dòng)相容性進(jìn)行識(shí)別。

表面肌電信號(hào)(surface electromyography,sEMG)是人體的一種生理信號(hào),能夠準(zhǔn)確地反映各個(gè)肌肉的實(shí)時(shí)狀態(tài)[4]。因此本文選擇sEMG 作為信號(hào)源,特征提取和模式識(shí)別是其中最重要的兩個(gè)環(huán)節(jié)。在眾多的sEMG 特征參數(shù)提取方法中,時(shí)域特征、頻域特征、時(shí)頻域特征和非線性特征被證明具有良好的識(shí)別能力。

分類器的設(shè)計(jì)是模式識(shí)別領(lǐng)域中非常重要的環(huán)節(jié),隱馬爾可夫模型、BP(back propagation)神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)(support vector machine,SVM)、極限學(xué)習(xí)機(jī)(extreme learning machine,ELM)[5]是目前常用的分類器,基于ELM 的分類模型,因其實(shí)施簡(jiǎn)單、訓(xùn)練時(shí)間短,已在助行機(jī)器人運(yùn)動(dòng)相容性識(shí)別領(lǐng)域取得了較好的識(shí)別結(jié)果,但是單核特征表達(dá)的單一性,限制了模型分類精度的進(jìn)一步提高。為了解決此問(wèn)題,本文首先從時(shí)域、頻域、時(shí)頻域角度提取了表面肌電信號(hào)特征,然后利用灰狼算法優(yōu)化極限學(xué)習(xí)機(jī)核函數(shù)參數(shù),最后用多核極限學(xué)習(xí)機(jī)理論,獲得最優(yōu)的分類模型。實(shí)驗(yàn)結(jié)果表明,基于多核學(xué)習(xí)ELM 的助行機(jī)器人運(yùn)動(dòng)相容性識(shí)別準(zhǔn)確率較單核極限學(xué)習(xí)機(jī)有明顯提高。

1 極限學(xué)習(xí)機(jī)

極限學(xué)習(xí)機(jī)可以看作是單隱含層網(wǎng)絡(luò),內(nèi)部為全連接。其中輸入層、隱含層和輸出層的神經(jīng)元數(shù)目分別為n和m,其結(jié)構(gòu)如圖1 所示。

圖1 極限學(xué)習(xí)機(jī)網(wǎng)絡(luò)結(jié)構(gòu)

設(shè)輸入層與隱含層、隱含層與輸出層之間的連接權(quán)值分別如下:

隱含層閾值為

激活函數(shù)用g(x) 表示,輸入樣本數(shù)目為N,則輸出Y∈Rm×N可以表示為

其中H為隱含層輸出矩陣。

當(dāng)g(x) 無(wú)限可微時(shí),ELM 的參數(shù)可在訓(xùn)練前任意指定,且在訓(xùn)練過(guò)程中保持不變。輸出層權(quán)值v由求解下式獲得。

式(6)的解為

其中H+為隱含層輸出矩陣E的廣義逆矩陣。

ELM 將神經(jīng)網(wǎng)絡(luò)訓(xùn)練過(guò)程簡(jiǎn)化為矩陣求逆問(wèn)題,學(xué)習(xí)速度明顯提高。

2 多核學(xué)習(xí)極限學(xué)習(xí)機(jī)

假定n個(gè)樣本,采用m個(gè)核映射函數(shù)gi映射得到m個(gè)特征值,βi表示對(duì)應(yīng)核函數(shù)gi的ELM 輸出權(quán)重,多核ELM[6-8]優(yōu)化公式為

其中,β=β1,β2,…,βp」 表示多個(gè)核函數(shù)融合矩陣的輸出權(quán)重,ξp表示第i個(gè)樣本的訓(xùn)練誤差,C表示平衡因子,γp表示第p個(gè)核函數(shù)對(duì)應(yīng)的特征值的融合系數(shù),yi表示第i個(gè)樣本標(biāo)簽。

式(8)的拉格朗日乘子方程為

其中α和τ表示拉格朗日乘子,最優(yōu)結(jié)構(gòu)參數(shù)α和核融合系數(shù)τ是多核學(xué)習(xí)ELM 的重要參數(shù),通過(guò)式(9)可求解多核ELM 最優(yōu)結(jié)構(gòu)參數(shù)α和核融合系數(shù)τ。

其中,K(..γ) 為多個(gè)核函數(shù)通過(guò)融合系數(shù)τ融合后的核矩陣,最后計(jì)算輸出權(quán)重βp。

3 基于多核學(xué)習(xí)ELM 的助行機(jī)器人運(yùn)動(dòng)相容性識(shí)別

3.1 表面肌電信號(hào)采集

根據(jù)文獻(xiàn)[9,10]的研究結(jié)果,闊筋膜張肌、股內(nèi)側(cè)肌、長(zhǎng)收肌和股直肌針對(duì)下肢運(yùn)動(dòng)識(shí)別效果較好,能夠得到清晰的信號(hào),因此選擇這4 塊肌肉采集sEMG 信號(hào)。步幅過(guò)大、過(guò)小和相容3 種運(yùn)動(dòng)模式的步幅分別為700 mm、350 mm、500 mm。實(shí)驗(yàn)對(duì)象共6 名,實(shí)驗(yàn)對(duì)象均自愿參加實(shí)驗(yàn)并簽署了知情同意書,無(wú)下肢神經(jīng)及肌肉骨骼病史,在實(shí)驗(yàn)前被告知實(shí)驗(yàn)的過(guò)程及要求。每位實(shí)驗(yàn)對(duì)象3 種情況各采集100 組數(shù)據(jù)。表面肌電信號(hào)采用高精度、高性能的Trigno TMWireless sEMG 無(wú)線肌電采集系統(tǒng),圖2 為sEMG 傳感器位置圖,圖3 為sEMG 信號(hào)采集與顯示。

圖2 傳感器位置圖

圖3 信號(hào)采集與顯示

3.2 信號(hào)預(yù)處理與特征提取

本文選擇巴特沃斯濾波器對(duì)信號(hào)預(yù)處理。肌電信號(hào)特征值的提取是運(yùn)動(dòng)相容性識(shí)別的關(guān)鍵環(huán)節(jié),本文采用時(shí)域特征提取法、頻域特征提取法、肌電信號(hào)非線性特征提取法及小波包特征提取法[11],每名受試者的3 種運(yùn)動(dòng)狀態(tài)都截取512 個(gè)數(shù)據(jù)點(diǎn)分析。小波包特征提取方法如下:參考移動(dòng)窗采樣方法,移動(dòng)窗大小為64 點(diǎn),每采集256 點(diǎn)后進(jìn)行一次小波包分解,將4 塊肌肉表面肌電信號(hào)作4 層db2 小波包分解,提取24個(gè)頻帶的信號(hào)并計(jì)算各頻帶信號(hào)能量作為特征值。排列組合熵主要是用來(lái)度量sEMG 的復(fù)雜度和非線性,通過(guò)使用排列組合熵可以評(píng)估和分析信號(hào)時(shí)間序列所呈現(xiàn)的一定規(guī)律性。本文選擇排列組合熵表示sEMG 非線性特征,計(jì)算過(guò)程參考文獻(xiàn)[12]的方法。表1 給出了特征值的基本意義,表2 給出了特征值的具體內(nèi)容。

表1 特征屬性的基本意義

表2 肌電信號(hào)特征

3.3 核函數(shù)的選擇

ELM 分類器需要選擇核函數(shù),線性核函數(shù)、多項(xiàng)式核函數(shù)、高斯核函數(shù)和Sigmoid 核函數(shù)是目前廣泛使用的核函數(shù)。

線性核函數(shù):

多項(xiàng)式核函數(shù):

高斯核函數(shù):

Sigmod 核函數(shù):

本文將對(duì)表2 中時(shí)域特征、頻域特征、時(shí)頻域等特征選用不同的核函數(shù)得到識(shí)別結(jié)果,通過(guò)分析對(duì)比來(lái)選擇最優(yōu)核函數(shù)。

3.4 核參數(shù)選擇

多核學(xué)習(xí)ELM 參數(shù)確定的過(guò)程相對(duì)簡(jiǎn)單,核函數(shù)參數(shù)是影響識(shí)別結(jié)果的重要參數(shù),本文采用灰狼優(yōu)化算法對(duì)核參數(shù)進(jìn)行尋優(yōu)。

灰狼算法是基于對(duì)自然界灰狼群的社會(huì)等級(jí)制度和圍獵行為的模擬[13-14]。在D維搜索空間中,定義種群X={X1,X2,…,Xn},其中n表示種群中灰狼的個(gè)體數(shù)。定義Xi=表示第i只灰狼的位置,其中表示第i只灰狼在第d維上的位置。將在遍歷種群的當(dāng)前最優(yōu)解記為α,將適應(yīng)度當(dāng)前排名次優(yōu)的解記為β,將適應(yīng)度當(dāng)前排名第3 的解記為δ,種群中的其他個(gè)體記為ω,獵物的位置即為優(yōu)化問(wèn)題的全局最優(yōu)解。

在圍獵過(guò)程中,狼群根據(jù)下式對(duì)獵物進(jìn)行包圍:

式中,Xp(t) 表示種群迭代至第t代時(shí)獵物的位置;X(t) 表示第t代灰狼個(gè)體的位置;常數(shù)C表示擺動(dòng)因子。

在我國(guó)現(xiàn)階段的市場(chǎng)經(jīng)濟(jì)條件下，社會(huì)對(duì)于人才的需求具有自發(fā)性以及多樣性的特點(diǎn)，在允許畢業(yè)生自主就業(yè)的基礎(chǔ)下，政府減少了對(duì)于高校的行政約束，允許高校根據(jù)社會(huì)市場(chǎng)的需求進(jìn)行自主設(shè)置專業(yè)和自主招生，這樣有利于充分發(fā)揮市場(chǎng)人才資源配置的靈活性和高效性，有利于實(shí)現(xiàn)人才供需的有效性結(jié)合。政府需要加強(qiáng)對(duì)于人才資源配置的宏觀調(diào)控，有利于保障人才資源和社會(huì)經(jīng)濟(jì)的協(xié)調(diào)發(fā)展，完善人才統(tǒng)計(jì)指標(biāo)體系等，促進(jìn)高校對(duì)于教育體系的深化改革，轉(zhuǎn)變傳統(tǒng)的教育模式，加強(qiáng)高校就業(yè)工作，維護(hù)社會(huì)穩(wěn)定持續(xù)發(fā)展。

其中,r1=rand(0,1)。

根據(jù)下式更新灰狼個(gè)體的位置:

其中,A表示收斂因子。

其中,α隨著迭代次數(shù)的增加,從2 遞減至0,r2=rand(0,1)。

在狼群中,α、β、δ最靠近獵物,灰狼優(yōu)化算法利用3 個(gè)個(gè)體來(lái)確定是否為最優(yōu)解。在灰狼群中,個(gè)體圍殺獵物行為的數(shù)學(xué)描述如下。

計(jì)算在第t代時(shí),α、β、δ對(duì)種群內(nèi)灰狼個(gè)體位置的影響X1、X2、X3,如下:

式中,Xα(t)、Xβ(t)、Xδ(t) 表示種群迭代至第t代時(shí)α、β、δ的位置,X(t) 表示第t代灰狼個(gè)體的位置,C1、C2、C3和A1、A2、A3表示不同的擺動(dòng)因子與收斂因子。

3.5 基于多核學(xué)習(xí)ELM 運(yùn)動(dòng)相容性識(shí)別流程

利用多核學(xué)習(xí)ELM 進(jìn)行運(yùn)動(dòng)相容性識(shí)別具體步驟如下。

(1)采集6 名實(shí)驗(yàn)對(duì)象3 種運(yùn)動(dòng)模式下的sEMG。

(2)利用巴特沃斯帶通濾波器對(duì)sEMG 濾波處理,并將處理后數(shù)據(jù)的50%選為訓(xùn)練樣本、50%選為測(cè)試樣本。

(3)采用3.2 節(jié)方法對(duì)預(yù)處理后的sEMG 從不同角度進(jìn)行特征提取。

(4)構(gòu)建多核學(xué)習(xí)ELM 分類模型,利用灰狼算法尋找分類器中最優(yōu)核函數(shù)參數(shù)。

(5)利用訓(xùn)練樣本訓(xùn)練多核學(xué)習(xí)ELM 分類器的性能,利用測(cè)試樣本檢驗(yàn)訓(xùn)練效果,得出識(shí)別結(jié)果。

識(shí)別流程如圖4 所示。

圖4 算法流程圖

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

本文按照3.1 節(jié)方法采集sEMG,利用巴特沃斯帶通濾波消除因零點(diǎn)漂移所產(chǎn)生的基線。圖5 為步幅過(guò)大、過(guò)小、相容時(shí)長(zhǎng)收肌信號(hào)圖。圖6 為步幅過(guò)大不相容運(yùn)動(dòng)模式中股內(nèi)側(cè)肌原始sEMG、巴特沃斯濾波sEMG、整流sEMG 波形對(duì)比。

圖5 步幅過(guò)大、過(guò)小、相容長(zhǎng)收肌信號(hào)圖

圖6 整流前后股內(nèi)側(cè)肌波形對(duì)比

對(duì)每組數(shù)據(jù)分別提取4 種sEMG 特征參數(shù),分別是時(shí)域特征、頻域特征、灰度矩特征和小波域特征,按照式(28)對(duì)特征參數(shù)進(jìn)行歸一化處理。

其中,x、xnew分別表示歸一化前后的數(shù)據(jù),xmax和xmin分別表示數(shù)據(jù)的最大值和最小值。

由表3 可見(jiàn),時(shí)域特征的最佳核函數(shù)是高斯核函數(shù);時(shí)頻域特征的最佳核函數(shù)是高斯核函數(shù);頻域特征的最佳核函數(shù)為Sigmod 核函數(shù);非線性特征的最佳核函數(shù)為多項(xiàng)式核函數(shù)。利用最優(yōu)核函數(shù)組合多核核函數(shù),表4 給出了最優(yōu)核組合。

表3 單核ELM 識(shí)別結(jié)果

表4 最優(yōu)核函數(shù)組合

根據(jù)表4 的最優(yōu)組合構(gòu)造多核ELM 分類模型,并對(duì)參數(shù)利用灰狼算法尋找最優(yōu)核函數(shù)參數(shù)。灰狼算法參數(shù)設(shè)置如下,初始種群規(guī)模為30,最大迭代次數(shù)為300,維度為20,運(yùn)行次數(shù)為20,灰狼算法中的反向?qū)W習(xí)優(yōu)化選擇概率p=0.5。采用5 折交叉驗(yàn)證的識(shí)別準(zhǔn)確率作為適應(yīng)度函數(shù)。經(jīng)過(guò)20 次計(jì)算求均值,表5 給出了最優(yōu)核函數(shù)參數(shù)。將參數(shù)優(yōu)化之后的多核學(xué)習(xí)ELM 分類器應(yīng)用于測(cè)試樣本,結(jié)果如圖7 所示?？梢钥闯?4 種單核ELM 用于不同特征參數(shù)的運(yùn)動(dòng)相容性模式識(shí)別時(shí),識(shí)別準(zhǔn)確率差異很大,即單一核函數(shù)應(yīng)用于不同特征空間時(shí),無(wú)法保證穩(wěn)定的高識(shí)別率;而多核學(xué)習(xí)ELM 不但保證了在每個(gè)特征空間提供最優(yōu)識(shí)別率,而且在時(shí)頻域特征空間中把識(shí)別正確率從92.6%提高94.7%。

表5 最優(yōu)核函數(shù)參數(shù)

圖7 不同ELM 核函數(shù)與多核學(xué)習(xí)ELM 對(duì)比實(shí)驗(yàn)結(jié)果

本文運(yùn)用SVM、多核學(xué)習(xí)SVM、ELM、多核學(xué)習(xí)ELM 4 種識(shí)別方法進(jìn)行運(yùn)動(dòng)相容性識(shí)別。為了更好地比較4 種方法的識(shí)別結(jié)果,選擇了Matlab 中的方差分析(analysis of variance,ANOVA)的統(tǒng)計(jì)學(xué)方法進(jìn)行分析,顯著性系數(shù)選擇均為P=0.05。

首先對(duì)6 位實(shí)驗(yàn)對(duì)象基于多核學(xué)習(xí)SVM、SVM、ELM、多核學(xué)習(xí)ELM 4 種識(shí)別方法的識(shí)別結(jié)果進(jìn)行比較。圖8 給出基于SVM、多核學(xué)習(xí)SVM、ELM、多核學(xué)習(xí)ELM 平均識(shí)別準(zhǔn)確率分別為87.00% ±3.66%、89.66% ±3.58%、91.37% ±4.31%、94.37%±4.31%。經(jīng)ANOVA 分析發(fā)現(xiàn),無(wú)論是SVM 還是多核SVM 均與多核學(xué)習(xí)ELM 方法存在顯著性差異(P＜0.01),可以看出多核學(xué)習(xí)ELM 模式識(shí)別方法比SVM、多核學(xué)習(xí)SVM 模式識(shí)別方法更適合相容性識(shí)別。對(duì)于ELM 與多核學(xué)習(xí)ELM 的識(shí)別結(jié)果進(jìn)行ANOVA 分析發(fā)現(xiàn)并無(wú)顯著性差別(P=0.9517),從實(shí)驗(yàn)結(jié)果來(lái)看,多核學(xué)習(xí)ELM 要略優(yōu)于ELM 方法。

圖8 4 種方法對(duì)比實(shí)驗(yàn)結(jié)果

為了更清晰地分析不同方法的準(zhǔn)確率差異性,將6 位實(shí)驗(yàn)對(duì)象的每一類運(yùn)動(dòng)的識(shí)別準(zhǔn)確率進(jìn)行統(tǒng)計(jì)平均,得到結(jié)果如表6 所示。從表6 中可以看出,僅對(duì)步幅過(guò)大的情況基于ELM 識(shí)別方法的準(zhǔn)確率略高于多核學(xué)習(xí)ELM 方法,對(duì)于步幅過(guò)小、步幅相容情況多核學(xué)習(xí)ELM 方法均有明顯的優(yōu)勢(shì),對(duì)于助行機(jī)器人運(yùn)動(dòng)相容性識(shí)別來(lái)說(shuō),多核學(xué)習(xí)ELM 方法性能更好。

表6 不同方法每一類運(yùn)動(dòng)的識(shí)別準(zhǔn)確率

4 結(jié)論

本文研究了一種基于多核學(xué)習(xí)ELM 的助行機(jī)器人運(yùn)動(dòng)相容性識(shí)別方法,把采集到的4 個(gè)特征空間的特征數(shù)據(jù)作為分析對(duì)象進(jìn)行助行機(jī)器人運(yùn)動(dòng)相容性識(shí)別研究,得出以下結(jié)論。

(1)多核學(xué)習(xí)ELM 融合了多種sEMG 特征,能夠較為全面地描述運(yùn)動(dòng)信息。

(2)多核學(xué)習(xí)ELM 分類器針對(duì)每個(gè)運(yùn)動(dòng)特征空間分別選取最優(yōu)核函數(shù),并利用灰狼算法實(shí)現(xiàn)了核參數(shù)的優(yōu)化配置。

(3)實(shí)驗(yàn)結(jié)果表明,多核學(xué)習(xí)ELM 分類器可以通過(guò)不同核函數(shù)融合任意多組運(yùn)動(dòng)特征信息,識(shí)別準(zhǔn)確率優(yōu)于ELM 和單核SVM,具有較好的實(shí)用性和推廣性。