腦電信號的小波變換和樣本熵特征提取方法

張毅，羅明偉，羅元

(1.重慶郵電大學 自動化學院，重慶 400065;2.重慶郵電大學光電工程學院，重慶 400065)

腦-機接口(brain-computer interface，BCI)是不依賴于大腦外周神經(jīng)與肌肉系統(tǒng)，在人腦和計算機或外部設備之間建立起來的一種通信系統(tǒng).它能夠為肢體殘疾但思維意識正常的患者提供一種新型的對外信息交流手段，并在殘疾人康復、正常人輔助控制、娛樂等領域有著廣泛的應用前景［1-4］.

特征提取是BCI研究的關鍵技術之一，常用的方法有 FFT［5］(fast Fourier transform)、AR［6-7］(autoregressive)、AAR［8］(adaptive auto-regressive)、小波變換［9-11］、樣本熵［12］等方法.FFT、AR 和 AAR 通過把幅度隨時間變化的腦電信號變換為腦電功率隨頻率變化的譜圖，從而提取出腦電信號的頻域特征.然而這些方法都只適合分析平穩(wěn)信號，對于非平穩(wěn)的腦電信號，具有很大的局限性.小波變換是一種典型的時頻分析法，具有多分辨率特性，在時域和頻域均具有良好的分辨率，適用于非平穩(wěn)信號處理.因此，小波變換非常分析合適腦電信號.然而，大腦是一個典型的非線性系統(tǒng)，腦電信號可以看作是它的輸出，小波變換卻不能反映出腦電信號的非線性特征.而樣本熵是一種非線性分析法，能通過度量腦電信號的復雜度來反映它的非線性特征，但不能反映出腦電信號的時頻特征.為了盡量充分地提取出腦電信號的特征，提高其正確的識別率，本文提出了一種結合小波變換和樣本熵的特征提取方法.該方法通過結合小波變換提取出的時頻特征和樣本熵分析提取出的非線性特征，得到了能反映出腦電信號時域、頻域和非線性特征的特征向量.利用支持向量機(support vector machine，SVM)對左、右手運動想象腦電信號進行分類，得到的最高正確識別率為91.43%，明顯高于僅采用小波變換、樣本熵以及其他傳統(tǒng)方法［6-7］的正確識別率.

1 基于Emotiv的腦電采集系統(tǒng)

本文所使用的腦電采集儀是如圖1所示的Emotiv，其采樣頻率為128 Hz.電極是按照國際10-20標準電極安放法安放，其安放位置如圖2所示.

圖1 Emotiv腦電采集儀Fig.1 Emotiv EEG acquisition instrument

圖2 Emotiv電極安放位置Fig.2 Electrode positions on Emotiv

圖中“CMS”和“DRL”為參考電極.選取4個健康受試者(年齡為22～24歲)，分別采集左、右手運動想象腦電信號.以一個受試者進行左手運動想象為例，對整個數(shù)據(jù)采集過程進行說明:實驗開始時，讓受試者靜坐并保持放松狀態(tài)，當受試者聽到口令時，開始想象左手運動，大約1 min后停止想象，休息2 s后再次進行同樣的操作，總共進行5次實驗.按上述實驗步驟，分別采集4個受試者進行左、右手運動想象的腦電信號.以持續(xù)2 s的腦電信號作為一組樣本，一個受試者進行一種運動想象可獲得150組樣本，每組樣本有256個數(shù)據(jù).

圖3是受試者進行左、右手運動想象的FC5通道腦電信號的波形圖.

圖3 想象左、右手運動的FC5通道腦電信號的波形Fig.3 The EEG signals of FC5channel for imagery left-right hands movement

2 結合小波變換和樣本熵的腦電信號特征提取

2.1 小波變換特征提取

人在想象單側手運動時，其對側相應初級感覺運動皮層區(qū)的腦電μ節(jié)律(8～12 Hz)和β節(jié)律(14～30 Hz)節(jié)律幅值降低，這種現(xiàn)象稱為事件相關去同步(event-related desynchronization，ERD);而同側腦電μ節(jié)律和β節(jié)律幅度升高，稱為事件相關同步(event-related synchronization，ERS).根據(jù)這一特征，可使用μ節(jié)律和β節(jié)律來分析左、右手運動想象腦電信號.而小波變換能把信號的整個頻帶劃分為多個子頻帶，因此可使用小波變換來分析左、右運動想象腦電信號.為了減少特征向量的維數(shù)，本文僅分析β節(jié)律.設x(n)表示實驗采集的EEG離散信號，則x(n)的離散小波變換定義為

式中:ψ(n)為小波基函數(shù)，j、k分別代表頻率分辨率和時間平移量.采用Mallat算法，對信號進行有限層分解，即

表1 小波分解的各層頻帶范圍Table 1 Frequency band range of each level of wavelet decomposition

小波系數(shù)能表達信號在時域和頻域的能量分布，因此利用小波系數(shù)的能量能反映出腦電信號的時域和頻域特征.由表1可知D2(16～32 Hz)在腦電信號的β節(jié)律頻帶范圍附近，因此，可提取對應于D2頻帶的小波系數(shù)的能量均值作為特征量.同時，為了進一步突出想象單側手運動引起的FC5、FC6通道腦電信號的幅值差異，還提取了小波系數(shù)的能量均值差PS作為特征量，即

式中:PFC5為FC5通道的能量均值，PFC6為FC6通道的能量均值.

對左、右手運動想象任務，各選取100組樣本數(shù)據(jù)計算，得到小波系數(shù)的能量均值和能量均值差分別如圖4和圖5所示.從圖4可以看出:想象左手運動時，F(xiàn)C5通道的能量均值PFC5大于FC6通道的能量均值PFC6;想象右手運動時，F(xiàn)C5通道的能量均值PFC5小于FC6通道的能量均值PFC6.因此，腦電信號經(jīng)小波分解后，其對應于β節(jié)律頻帶的小波系數(shù)的能量均能作為區(qū)分左、右手運動想象腦電信號的特征量.從圖5可知，利用小波系數(shù)的能量值均值差提取得到的特征很好地反映出了左、右手運動想象腦電信號的差別.

圖4 小波系數(shù)的能量均值Fig.4 The average energy of wavelet coefficient

圖5 小波系數(shù)的能量均值差Fig.5 The average energy difference of wavelet coefficient

2.2 結合小波變換和樣本熵的特征提取

小波變換能反映出腦電信號的時頻特征，但不能反映出它的非線性特征.而樣本熵是一種非線性分析法，能通過度量腦電信號的復雜度來反映它的非線性特征，但卻不能反映出它的時頻特征.針對這2種方法的不足，本文將結合小波變換和樣本熵來提取腦電信號的特征，進而提高腦電信號的正確識別率.

樣本熵是在近似熵的基礎上提出的一種改進的復雜度測量方法［13］.它不但具備近似熵所有的優(yōu)點，而且避免了近似熵中統(tǒng)計量不一致的問題.

設一維時間序列為{X(i)}，i=1，2，…，N，其樣本熵可通過如下計算得到:

1)將序列{X(i)}按順序組成m維矢量，即2)對每一個i計算矢量Ym(i)與矢量Ym(j)之間的距離，即

3)給定閾值 r(r＞0)，對每一個 i值統(tǒng)計d［Ym(i)，Ym(j)］小于r的數(shù)目，記為Ami(r)，并記此數(shù)目與總的距離數(shù)目N-m的比值為Bmi(r)，即

5)將序列{X(i)}按順序組成m+1維矢量，再重復1)～4)，得到Bm+1i(r)和Bm+1(r).

6)序列{X(i)}的樣本熵為

在實際計算中，由于序列長度有限，因此最后得到序列長度為N時的樣本熵的估計值

顯然，樣本熵的值與m、r和N有關.通常取m=2，r=0.1～0.25 SD(SD 表示原始時間序列的標準差).在本文中取 m=2，r=0.2 SD，N=256.

利用樣本熵提取腦電信號的非線性特征的具體步驟為:首先對腦電信號加滑動時間窗，以1 s長度為滑動時間窗(即128點)，計算腦電信號的樣本熵，窗口每次移動一個采樣點，并計算下1 s時間窗的腦電信號的樣本熵，直到計算出最后1 s時間窗的腦電信號的樣本熵為止，從而獲得該樣本數(shù)據(jù)中腦電信號樣本熵的時間序列，然后將這一組樣本熵序列疊加平均，即獲得一組樣本數(shù)據(jù)的樣本熵.

對左右手運動想象的FC5通道腦電信號，各選取100組樣本數(shù)據(jù)計算，得到FC5通道腦電信號的樣本熵如圖6所示.從圖6知，想象左、右手運動的腦電信號的樣本熵雖有部分重疊，但仍具有一定的區(qū)分度.由此可見，樣本熵在一定程度上能反映左、右手運動想象腦電信號的特征.

圖6 FC5通道腦電信號的樣本熵Fig.6 The sample entropy of EEG signals of FC5channel

對采集的FC5、FC6通道的左、右手運動想象腦電信號，用樣本熵特征提取方法可得到二維特征向量，用小波變換可得到三維特征向量，結合2種方法可得到五維特征向量，即為本文的特征提取向量.

2.3 SVM 分類

SVM是一種基于統(tǒng)計學習理論的模式識別方法，能有效地避免傳統(tǒng)分類方法所存在的過學習、維數(shù)災難、局部極小等問題，并且在小樣本條件下仍然具有良好的泛化能力.本文將利用SVM對左、右手運動想象腦電信號進行分類.首先選取SVM的核函數(shù)，然后利用訓練樣本計算出SVM的參數(shù)(如拉格朗日乘子a、最優(yōu)超平面的法向量w和偏置值 b等)，并確定判別函數(shù)f(x)，最后輸入測試樣本，根據(jù)sgn(f(x))的值，輸出類別.如果sgn(f(x))為1，則判別結果為想象左手運動的腦電信號;如果sgn(f(x))為-1，則判別結果為想象右手運動的腦電信號.

3 實驗結果分析

在采集的150組樣本數(shù)據(jù)中，選取75組作為訓練樣本，剩下75組作為測試樣本.把用小波變換提取的三維特征向量、用樣本熵提取的二維特征向量以及兩者相結合的五維特征向量分別輸入SVM分類器中，得到不同受試者在不同的特征提取方法下的正確識別率見表2，其中SVM的核函數(shù)采用的是線性核函數(shù).

表2 不同受試者在不同的特征提取方法下的正確識別率Table 2 The right recognition rates of different subjects using different feature extraction methods %

從表2可以看出，在對左、右手運動想象腦電信號的識別中，使用小波變換和樣本熵相結合的特征提取算法得到的最高正確識別率為91.43%，僅使用小波變換和樣本熵得到的最高正確識別率分別為88.00%和 75.33%，而在文獻［6-7］中，僅使用 AR特征提取方法得到的最高正確識別率分別85.00%和75%.因此，使用小波變換和樣本熵相結合的特征提取算法的正確識別率明顯高于僅采用小波變換或樣本熵的特征提取算法的正確識別率，同時也高于傳統(tǒng)的使用AR特征提取方法的正確識別率.

4 結束語

本文提出了一種結合小波變換和樣本熵的特征提取方法.分別用小波變換和樣本熵對左、右手運動想象腦電信號進行了分析，并抽取出相應的特征組成特征向量，輸入SVM分類器中，實現(xiàn)了左、右手運動想象腦電信號的分類.研究結果表明，結合小波變換和樣本熵的特征提取算法對左右手運動想象腦電信號的正確識別率明顯高于僅采用小波變換、樣本熵以及其他傳統(tǒng)特征提取算法的正確識別率，其最高正確識別率達到91.43%.因此，該方法在BCI中具有應用價值.

［1］MUGLERR E，BENSCH M，HALDER S，et al.Control of an internet browser using the P300 event-related potential［J］.Int J Bioelectromagn，2008，10(1):56-63.

［2］王斐，張育中，寧廷會，等.腦—機接口研究進展［J］.智能系統(tǒng)學報，2011，6(3):189-199.

WANG Fei，ZHANG Yuzhong，NING Tinghui，et al.Research progress in a brain-computer interface［J］.CAAI Transactions on Intelligent Systems，2011，6(3):189-199.

［3］TANAKA K，MATSUNAGA K，WANG H.Electroencephalogram-based control of an electric wheelchair［J］.IEEE Trans Robotics，2005，21(4):762-766.

［4］REBSAMEN B，GUAN G T，ZHANG H H，et al.A brain controlled wheelchair to navigate in familiar environments［J］.IEEE Trans Neural Syst Rehabil Eng，2010，18(6):590-598.

［5］POLAT K，GüNES S.Artificial immune recognition system with fuzzy resource allocation mechanism classifier，principal component analysis and FFT method based new hybrid automated identification system for classification of EEG signals［J］.Expert System with Applications，2008，34(3):2039-2048.

［6］徐寶國，宋愛國.單次運動想象腦電的特征提取和分類［J］.東南大學學報:自然科學版，2007，37(4):629-633.

XU Baoguo，SONG Aiguo.Feature extraction and classification of single trial motor imagery EEG［J］.Journal of Southeast University:Natural Science Edition，2007，37(4):629-633.

［7］張毅，楊柳，李敏，等.基于AR和SVM的運動想象腦電信號識別［J］.華中科技大學學報:自然科學版，2011，39(suppl.2):103-106.

ZHANG Yi，YANG Liu，LI Min，et al.Recognition of motor imagery EEG based on AR and SVM［J］.Journal of Huazhong University of Science＆Technology:Natural Science Edition，2011，39(suppl.2):103-106.

［8］PFURSTCHELLER G，NEUPER C.Motor imagery and direct brain-computer communication［J］.Proceedings of the IEEE，2001，89(7):1123-1134.

［9］李明愛，王蕊，郝冬梅.想象左右手運動的腦電特征提取及分類研究［J］.中國生物醫(yī)學工程學報，2009，28(2):166-170.

LI Ming’ai，WANG Rui，HAO Dongmei.Feature extraction and classification of EEG for imagery left-right hands movement［J］.Chinese Journal for Biomedical Engineering，2009，28(2):166-170.

［10］RAFIEE J，RAFIEE M A，PRAUSE N，et al.Wavelet basis functions in biomedical signal processing［J］.Expert Systems with Applications，2011，38(5):6190-6201.

［11］徐寶國，宋愛國，費樹岷.在線腦機接口中腦電信號的特征提取與分類方法［J］.電子學報，2011，39(5):1025-1030.

XU Baoguo，SONG Aiguo，F(xiàn)EI Shumin.Feature extraction and classification of EEG in online brain-computer interface［J］.Acta Electronica Sinica，2011，39(5):1025-1030.

［12］周鵬，葛家怡，曹紅寶，等.基于樣本熵的運動想象分類研究［J］.信息與控制，2008，37(2):191-196.

ZHOU Peng，GE Jiayi，CAO Hongbao，et al.Classifica-tion of motor imagery based on sample entropy［J］.Information and Control，2008，37(2):191-196.

［13］RICHMAN J S，MOORMAN J R.Physiological time-series analysis using approximate entropy and sample entropy［J］.American Journal of Physiology—Heart and Circulatory Physiology，2000，278(6):H2039-H2049.