[發(fā)明專利]一種腦電圖信號(hào)分類方法

無效

• 申請(qǐng)?zhí)枺?/td>: 201210344890.9
• 申請(qǐng)日: 2012-09-17
• 公開（公告）號(hào): CN102824173A
• 公開（公告）日: 2012-12-19
• 發(fā)明（設(shè)計(jì)）人: 李志強(qiáng);郝新紅;栗蘋;于成大;閆曉鵬;梁營
• 申請(qǐng)（專利權(quán)）人: 北京理工大學(xué)
• 主分類號(hào): A61B5/0476
• 分類號(hào): A61B5/0476
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種腦電圖信號(hào)分類方法，其特征在于，包括以下步驟：

一、輸入腦電圖信號(hào)，并對(duì)輸入的腦電圖信號(hào)進(jìn)行經(jīng)驗(yàn)?zāi)B(tài)分解，把腦電圖信號(hào)分解成本征模函數(shù)的和；

二、對(duì)每個(gè)本征模函數(shù)c(t)=a(t)cosθ(t)進(jìn)行經(jīng)驗(yàn)調(diào)幅-調(diào)頻分解，得到經(jīng)驗(yàn)調(diào)頻分量F(t)=cosθ(t)；

三、判斷得到的經(jīng)驗(yàn)調(diào)頻分量中是否含有騎行波；

四、如果含有騎行波，進(jìn)行去騎行波處理；

五、通過計(jì)算經(jīng)驗(yàn)調(diào)幅分量；

六、計(jì)算經(jīng)驗(yàn)調(diào)頻分量的正交分量，為了得到四個(gè)象限的相位值以便解卷繞求解瞬時(shí)相位，對(duì)的符號(hào)根據(jù)三角函數(shù)特點(diǎn)做以下處理：如果F(t)單調(diào)遞增則取否則取這樣相位便可對(duì)應(yīng)四個(gè)象限的值；

七、計(jì)算瞬時(shí)相位：由經(jīng)驗(yàn)調(diào)頻分量及其正交分量反正切解卷繞得到瞬時(shí)相位θ(t)=unwrap[arctansinθ(t)F(t)];]]>

八、由步驟五得到的調(diào)幅分量和步驟七得到的瞬時(shí)相位，按照Cohen帶寬公式計(jì)算調(diào)幅帶寬和調(diào)頻帶寬：

<ω>=1E∫-∞+∞θ′(t)a2(t)]]>

B2AM=1E∫-∞+∞a′2(t)dt]]>

B2FM=1E∫-∞+∞(θ′(t)-<ω>)2a2(t)dt]]>

九、將步驟八計(jì)算得到的調(diào)幅帶寬和調(diào)頻帶寬作為支持向量機(jī)的輸入，對(duì)腦電圖信號(hào)進(jìn)行分類。

2.根據(jù)權(quán)利要求1所述的一種腦電圖信號(hào)分類方法，其特征在于，步驟四所述去騎行波處理包括以下步驟：

1)首先判斷其是否為孤立騎行波，判斷方法：對(duì)于上凸或下凸的騎行波，如果其右端點(diǎn)為下一個(gè)上凸或下凸的騎行波的左端點(diǎn)，則該騎行波為首尾相連非孤立的騎行波，否則為孤立的騎行波；

2)如果是孤立騎行波則直接進(jìn)行騎行波翻轉(zhuǎn)，即沿著連接騎行波的兩個(gè)端點(diǎn)的直線翻轉(zhuǎn)騎行波，假設(shè)s(t)為一段騎行波，l(t)為連接兩端的直線，騎行波s(t)經(jīng)翻轉(zhuǎn)后為2l(t)-s(t)；然后找到被翻轉(zhuǎn)的騎行波的極值，用該極值的絕對(duì)值去除騎行波左側(cè)第一個(gè)過零點(diǎn)到騎行波右側(cè)第一個(gè)過零點(diǎn)之間的波形數(shù)據(jù)，經(jīng)這一過程處理后不僅該孤立騎行波被去除掉而且局部的數(shù)據(jù)被歸一化；

3)如果步驟1)判斷不是孤立的騎行波而是連續(xù)出現(xiàn)幾個(gè)騎行波首尾相連，首先翻轉(zhuǎn)每個(gè)騎行波；然后選擇最左端騎行波的左端點(diǎn)和最右端騎行波的右端點(diǎn)再加上兩端點(diǎn)間的另一個(gè)點(diǎn)共三個(gè)點(diǎn)進(jìn)行拋物線插值擬合，擬合出來的拋物線作為母波的一部分；之后用拋物線的極值的絕對(duì)值去除最左端騎行波左側(cè)第一個(gè)過零點(diǎn)和最右側(cè)騎行波右側(cè)第一個(gè)過零點(diǎn)之間的波形數(shù)據(jù)，經(jīng)這一過程處理后不僅該首尾相連的騎行波被去除掉而且局部的數(shù)據(jù)被歸一化；

4)最后檢查是否還有騎行波，如果有則返回步驟1)重復(fù)執(zhí)行，否則退出。