[發明專利]一種多尺度補償傳遞熵的皮層肌肉功能耦合方法

在審

• 申請號: 202011383022.2
• 申請日: 2020-12-01
• 公開（公告）號: CN112617858A
• 公開（公告）日: 2021-04-09
• 發明（設計）人: 佘青山;張敏;孟明;張啟忠
• 申請（專利權）人: 杭州電子科技大學
• 主分類號: A61B5/369
• 分類號: A61B5/369;A61B5/389;A61B5/397;A61B5/372;A61B5/374;A61B5/00
• 代理公司: 浙江千克知識產權代理有限公司 33246
• 代理人: 周希良

1.一種多尺度補償傳遞熵的皮層肌肉功能耦合方法，其特征在于：該方法包括以下主要步驟：

步驟(1)，多通道腦肌電信號同步采集；

具體為：同時采集受試者者左右手在不同抓握動作下的多通道腦肌電信號；

步驟(2)，對同步采集的腦肌電信號采用APIT-MEMD方法進行時頻尺度化；

具體為：首先構建兩組腦電信號X＝{x₁,x₂,…x_i,…,x_M}和肌電信號Y＝{y₁,y₂,…y_i,…,y_M}時間序列，M為時間序列的長度，其次采用APIT-MEMD方法對兩組時間序列進行分解，各得到K個窄帶IMF分量，與時頻尺度相關的兩組腦肌電時間序列分別為和

步驟(3)，計算腦肌電信號在不同尺度上的McTE值；

具體為：兩個表示物理系統x、y隨時間變化的隨機過程X、Y，令x_t和y_t是描述X和Y在時間t處訪問狀態的隨機變量，而x_n:t是組成向量從時間n到時間t的所有x樣本，則補償傳遞熵cTE的定義為：

cTE_X→Y＝H(y_t|y_1:t-1,x_t)-H(y_t|y_1:t-1,x_1:t) (1)

構造信號到信號的多尺度補償傳遞熵McTE_EEG→EMG：

其中McTE_EEG→EMG表示EEG的k₁分量到EMG的k₂分量的多尺度補償傳遞熵值

同理可得：

其中，McTE_EMG→EEG表示EMG的k₂分量到EEG的k₁分量的多尺度補償傳遞熵值；

步驟(4)，采用步驟(2)和(3)的計算結果進行不同耦合方向、不同頻段間的皮層肌肉耦合分析。