1.一種基于多元經(jīng)驗?zāi)J椒纸獾募‰娦盘柾教幚矸椒?,主要包括以下步驟:
步驟1:通過p,p≥3個通道肌電信號采集儀器同步拾取人體在特定運動動作條件下典型肌肉群處的肌電信號,獲取肌電信號樣本數(shù)據(jù),構(gòu)成一個p維原始信號s(t),s(t)={s1(t),s2(t),…,sp(t)},其中t是時刻值;
步驟2:通過低差異的哈默斯利序列,將步驟1所獲取的s(t)映射到一個p-1維球體上形成新的數(shù)據(jù)點集up={u1,u2,…,up};
步驟3:將步驟2得到的點集,沿著方向向量集其中,p-1維球體上角度點集θv={θv1,θv2,…,θp-1|θv∈Rp},方向向量v=1,2,…,V分別進行投影并得到一組投影集
步驟4:計算步驟3得到的投影集的極值及對應(yīng)的時刻
步驟5:將步驟4得到的極值矩陣插值,得到包絡(luò)曲線集
步驟6:將步驟5得到的包絡(luò)曲線集求平均,得到平均包絡(luò)曲線,
步驟7:通過公式(5)提取成分d(t):
此時,若d(t)滿足停止條件,則將d(t)被定義為固有模態(tài)分量c(t),同時將s(t)-d(t)視為s(t),重復(fù)以上步驟1至7;若d(t)不滿足停止條件,則將d(t)視為s(t),重復(fù)以上步驟1至7得到殘余分量r(t),停止循環(huán)計算;
通過以上步驟重復(fù)篩選后,原始信號s(t)將被分解為N個固有模態(tài)分量,{cN(t)},N=1,2,…,N和殘余分量r(t)之和,如公式(6)所示:
其中,篩選停止條件是其中是平均包絡(luò)曲線的變化值,a(t)是包絡(luò)平均曲線的幅值,而σ是調(diào)節(jié)參數(shù),其值根據(jù)實際情況決定且σ∈(0,1)。
2.如權(quán)利要求1所述的一種基于多元經(jīng)驗?zāi)J椒纸獾募‰娦盘柾教幚矸椒?,其特征在于所述步驟2的獲得低差異的哈默斯利序列的具體方法為:
首先將獲得的多通道肌電信號數(shù)據(jù)構(gòu)建Halton序列,即假設(shè)整個質(zhì)數(shù)集由小到大排列后,前n個質(zhì)數(shù)為{x1,x2,…,xn},則某一維度上的Halton序列的第i個樣本為:
其中,基數(shù)x為整數(shù),aj(j=1,2,…,s)是樣本索引i按x進制展開后的j位上的數(shù),s表示樣本索引i按x進制展開后的總個數(shù),因此樣本索引i通過基數(shù)x的表示方法如下:
i=a0+a1×x+a2×x2+…+as×xs (2)
將質(zhì)數(shù)分別帶入公式(2),得到n個樣本的Halton序列如下:
由Halton序列得到n個樣本的哈默斯利序列如下:
3.如權(quán)利要求1所述的一種基于多元經(jīng)驗?zāi)J椒纸獾募‰娦盘柾教幚矸椒ǎ涮卣髟谟谒霾襟E3的具體方法為:
將哈默斯利序列采樣得到的(p-1)維球體上形成新的數(shù)據(jù)點集up={u1,u2,…,up}沿著(p-1)維球體的所有方向向量分別投影到球體上,得到一組投影集此時,投影集長度為單位長度,不具有物理意義,投影集方向向量值代表了up中各個元素。
4.如權(quán)利要求1所述的一種基于多元經(jīng)驗?zāi)J椒纸獾募‰娦盘柾教幚矸椒?,其特征在于所述步驟4的方法為:
通過設(shè)定檢索數(shù)據(jù)長度,并計算過零點個數(shù),得出投影集對應(yīng)的極值及時刻
5.如權(quán)利要求1所述的一種基于多元經(jīng)驗?zāi)J椒纸獾募‰娦盘柾教幚矸椒?,其特征在于所述本發(fā)明肌電信號包括通過介入和非介入兩種方式拾取的肌電電信號。