基于小波變換與Fisher準則的P300特征提取方法
【技術領域】
[0001] 本發(fā)明屬于認知神經(jīng)科學領域與信息技術領域的結(jié)合應用，涉及一種事件相關電 位P300特征提取方法，具體是基于小波變換與Fisher準則的P300特征提取方法。
【背景技術】
[0002] 腦機接口是一種為運動功能缺失而腦功能完好的患者，提供與外界交流的新路 徑。P300Speller是腦機接口的一種方式，其功能是通過分析用戶的腦電信號，識別出其 所希望輸出的字符，從而幫助用戶與外界交流。目前在P300Speller中使用的特征提取方 法為：對從各通道提取出的EEG數(shù)據(jù)段做下采樣，得到的即作為特征。這種特征提取方法存 在特征不明確的問題，從而影響到分類器效果，為了保證其準確率，通常需要較多次刺激重 復，一般為15次，存在的問題是，一方面被試易感到疲勞，另一方面輸出一個字符需要較多 時間，從而降低了字符傳輸速率。
[0003] 由于P300成分主要存在于低頻部分，目前已發(fā)表的利用小波變換的一些論文， 一般是對腦電數(shù)據(jù)進行多層小波變換。例如在論文《Off-lineanalysisoftheP300 event-relatedpotentialusingdiscretewavelettransform〉〉中對腦電數(shù)據(jù)進行 6 層 離散小波變換，將得到小波系數(shù)的近似部分作為特征?；蚴菍δX電信號去噪再重構(gòu)信號，提 取特征。但是兩者的實用性都不高。

【發(fā)明內(nèi)容】

[0004] 為解決了現(xiàn)有技術中腦電信號特征提取不明確、實用性不高的缺陷，本發(fā)明提供 一種基于小波變換與Fisher準則相結(jié)合的P300特征提取方法，該方法通過提取恰當特征， 使分類器的效果得到提升，在滿足準確率要求的前提下，能夠達到減少刺激重復次數(shù)、提高 字符傳輸速率的目的。
[0005] 本發(fā)明采用以下技術方案實現(xiàn)：一種基于小波變換與Fisher準則的P300特征提 取方法，其特征在于，包括以下步驟：步驟S1 :對EEG數(shù)據(jù)經(jīng)過預處理后，根據(jù)用戶設定的傳 輸通道，提取每個刺激后長度為4的數(shù)據(jù)段，記為向量把給定小波所對應的低通濾波器 記為向量h，其長度為12 ;設置相關參數(shù)：小波分解層數(shù)L、單通道特征個數(shù)r;步驟S2 :根據(jù) L、&、12和h確定wxB矩陣W;步驟S3 :把所有的e延拓為n維向量，仍用遺示；步驟S4 : 根據(jù)Fisher準則，由W和r為每個通道分別確定一個rxi!矩陣M ;步驟S5 :對所有EEG數(shù) 據(jù)段e按通道求特征向量：
，^為r維向量；步驟S6 :將各傳輸通道得到的特征向量 拼接構(gòu)成總特征向量。
[0006] 進一步的，所述步驟S2包括以下具體步驟：步驟S21 :k=l,


;步驟S22 :若k等于L，則輸出W;否則進入下一步；步驟S23 :
，m=p;步驟S24:把W延拓 為gXB矩陣，仍記為W;步驟S25 :生成矩陣A:
，
;步驟S26 :用AxTST 更新W，k=k+l，返回步驟S22。
[0007] 在本發(fā)明一實施例中，步驟S4包括以下具體步驟：步驟S41:把該通道的所有 為目標組G+和非目標組G、其中G+中的N+個e為目標刺激對應的數(shù)據(jù)段，G沖的『個e為非目標刺激對應的數(shù)據(jù)段；步驟S42 :利用公式*1=1%把所有的e轉(zhuǎn)化為4相應地，〇分 屬于G+和G、步驟S43 :計算

為m維向量；步驟S44 :將
按絕對值由大到小排序，得到
，表示w中 各元素的原位置；步驟S45:根據(jù)設定的單通道特征個數(shù)r，取出W的第q，q，___，cv行， 構(gòu)成rxn特征提取矩陣M，其中ISrim。
[0008] 與現(xiàn)有技術相比，本發(fā)明具有以下優(yōu)點：結(jié)合了小波變換與Fisher準則兩種方法 來確定特征。通過利用小波變換，將原始EEG數(shù)據(jù)映射到小波域上，在小波域上應用Fisher 準則尋找能將兩類數(shù)據(jù)很好地區(qū)分開的投影軸，并只取兩類差異最大的r維投影空間，作 為最終特征空間，從而實現(xiàn)了對高維特征空間的降維，此外所得到的特征空間，是能明顯表 現(xiàn)出兩類數(shù)據(jù)的差別，有利于分類器性能的提高，減少刺激重復次數(shù)，提高字符傳輸速率。 并且通過小波變換矩陣來為每個傳輸通道構(gòu)造對應特征提取矩陣，采用特征提取矩陣求各 通道特征，具有在線運算速度快的優(yōu)點，使本發(fā)明能夠適于實際應用。
【附圖說明】
[0009] 圖1為本發(fā)明的主要步驟流程圖。
【具體實施方式】
[0010] 下面結(jié)合附圖和【具體實施方式】對本發(fā)明做進一步說明。
[0011] 本發(fā)明提供一種基于小波變換與Fisher準則的P300特征提取方法，用于P300 Speller特征提取，對每個刺激對應的EEG數(shù)據(jù)段提取恰當特征，提高準確率，減少刺激重 復次數(shù)，提高字符傳輸速率。
[0012] 本發(fā)明的主要步驟流程圖參見圖1。該方法包括以下步驟： 步驟S1 :對EEG數(shù)據(jù)經(jīng)過預處理后，根據(jù)用戶設定的傳輸通道，提取每個刺激后長度為li的數(shù)據(jù)段，記為向量把給定小波所對應的低通濾波器記為向量h，其長度為;設置相 關參數(shù)：小波分解層數(shù)L、單通道特征個數(shù)r; 步驟S2 :根據(jù) 矩陣w; 步驟S3 :把所有的e延拓為n維向量，仍用e表示； 步驟S4 :根據(jù)Fisher準則，由W和r為每個通道分別確定一個rxn矩陣M; 步驟S5 :對所有EEG數(shù)據(jù)段e按通道求特征向量：v=Afe，r為r維向量； 步驟S6 :將各傳輸通道得到的特征向量拼接構(gòu)成總特征向量。
[0013] 進一步的，步驟S2包括以下具體步驟： 步驟S21:k=l，
；初始化W
步驟S22 :若k等于L，則輸出W;否則進入下一步； 步驟S23 :
，m=p; 步驟S24:把W延拓為矩陣，仍記為W; 步驟S25:生成矩陣A
步驟S26 :用Ax曹更新W，k=k+l，返回步驟S22。
[0014] 在本發(fā)明一實施例中步驟S4包括以下具體步驟： 步驟S41:把該通道的所有為目標組G+和非目標組(T，其中G+中的N+個e為目標 刺激對應的數(shù)據(jù)段，G沖的N，e為非目標刺激對應的數(shù)據(jù)段； 步驟S42 :利用公式il=IFe把所有的禱?；癁?相應地，〇分屬于G+和(T; 步驟S43 :計_

為m維向量； 步驟S44:將
按絕對值由大到小排序，得到
，其中
，表示AJ中各元素的原位置； 步驟S45:根據(jù)設定的單通道特征個數(shù)r，取出W的第<%，4，_"，行，構(gòu)成rxH特 征提取矩陣M，其中ISrSiii。
[0015] 在本發(fā)明一實施例中，選擇小波為"db4"小波，即可確定對應的低通濾波器h及其 長度12,對參數(shù)L、&、r可按如下設定：L=4, & =800,r=15。
[0016] 以上僅為本發(fā)明的較佳實施方式，不能因此而理解為對本發(fā)明專利范圍的限制。 應當指出的是，對于本領域的普通技術人員來說，在不脫離本發(fā)明構(gòu)思的前提下，還可以做 出若干變形和改進，這些都屬于本發(fā)明的保護范圍。因此，本發(fā)明專利的保護范圍應以所附 權利要求為準。
【主權項】
1. 一種基于小波變換與Fisher準則的P300特征提取方法，其特征在于，包括以下步 驟： 步驟Sl :對EEG數(shù)據(jù)經(jīng)過預處理后，根據(jù)用戶設定的傳輸通道，提取每個刺激后長度為 的數(shù)據(jù)段，記為向量把給定小波所對應的低通濾波器記為向量h，其長度為?2 ;設置相 關參數(shù)：小波分解層數(shù)U單通道特征個數(shù)r ; 步驟S2 :根據(jù)L、I1、?3和h確定 JRXlS 矩陣w; 步驟S3 :把所有的e延拓為η維向量，仍用e表示； 步驟S4 :根據(jù)Fisher準則，由W和r為每個通道分別確定一個m矩陣M ; 步驟S5 :對所有EEG數(shù)據(jù)段e按通道求特征向量：V = Me，r為r維向量； 步驟S6 :將各傳輸通道得到的特征向量拼接構(gòu)成總特征向量。2. 根據(jù)權利要求1所述的基于小波變換與Fisher準則的P300特征提取方法，其特征 在于：所述步驟S2包括以下具體步驟： 步驟S21:k=l，；初始化1:，先令〇 ? I =1,2,---,?* ? J=W---,η， 再令= ?，?=:!，】,--、》!*，; 步驟S22 :若k等于L，則輸出W ;否則進入下一步； 步驟 S23 :, m=p ； 步驟S24:把W延拓為矩陣，仍記為W;步驟S25 :生成矩陣A: ，先令= 〇，，-， J 二 12，凡 再令 =乓，?: = - - ,|1，? = Xt - - ,?2 ; 步驟S26 :用AxW更新W，k=k+l，返回步驟S22。3. 根據(jù)權利要求1所述的基于小波變換與Fisher準則的P300特征提取方法，其特征 在于：步驟S4包括以下具體步驟： 步驟S41 :把該通道的所有為目標組G +和非目標組0-，其中G+中的N +個e為目標 刺激對應的數(shù)據(jù)段，G沖的N，e為非目標刺激對應的數(shù)據(jù)段； 步驟S42 :利用公式rf = Iife把所有的e轉(zhuǎn)化為忒相應地，〇分屬于G +和0-; 步驟S43 :計;S，Hf為m維向量； 步驟S44 :將》二按絕對值由大到小排序，得到[氣》氣氣】，其中 表示講中各元素的原位置； 步驟S45:根據(jù)設定的單通道特征個數(shù)r，取出W的第C1, C2，_"，c,行，構(gòu)成rxn特 征提取矩陣M，其中IirSm。
【專利摘要】本發(fā)明以小波變換與Fisher準則相結(jié)合為基礎提供了一種針對P300 Speller腦機接口的腦電信號特征提取方法。主要包括以下步驟：根據(jù)用戶指定的小波函數(shù)及小波分解層數(shù)，為特定長度的EEG數(shù)據(jù)段構(gòu)造小波變換矩陣；利用小波變換矩陣，將EEG數(shù)據(jù)段映射到小波域，在小波域中利用Fisher準則尋找使得不同類別數(shù)據(jù)區(qū)分度最大的投影軸；利用投影軸從小波變換矩陣中抽取若干行構(gòu)成特征提取矩陣；每個通道的特征提取矩陣把對應的EEG數(shù)據(jù)段映射為特征向量；每次刺激在各通道上的特征向量被拼接為一個特征向量。該方法為每次刺激對應的EEG數(shù)據(jù)段計算特征向量，在滿足準確率要求的前提下，能減少刺激重復次數(shù)，從而提高字符傳輸速率。
【IPC分類】G06K9/00, G06K9/62
【公開號】CN104899573
【申請?zhí)枴緾N201510339891
【發(fā)明人】黃志華, 郭順英, 林蘇云, 文宇坤
【申請人】福州大學
【公開日】2015年9月9日
【申請日】2015年6月18日