王宇丁，陳民鈾，張 莉

(重慶大學(xué) 輸配電裝備及系統(tǒng)安全與新技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，重慶400044)

0 引 言

自從發(fā)現(xiàn)腦電信號(hào)以來(lái)，科學(xué)家就一直努力通過(guò)腦電信號(hào)了解大腦的活動(dòng)規(guī)律。隨著人體解剖學(xué)、計(jì)算機(jī)技術(shù)、通信技術(shù)、電子技術(shù)、心理學(xué)等學(xué)科的發(fā)展，出現(xiàn)了一個(gè)新的研究領(lǐng)域，這就是腦-機(jī)接口(brain－computer interface，BCI)技術(shù)［1～4］。BCI 技術(shù)是利用一些有規(guī)律的腦電信號(hào)特征，實(shí)現(xiàn)大腦與計(jì)算機(jī)或者其他電子設(shè)備的通信［5～7］。作為一項(xiàng)跨學(xué)科的技術(shù)，BCI 的研究涉及的領(lǐng)域廣泛，包括微電子學(xué)、生物學(xué)、計(jì)算機(jī)、通信工程、哲學(xué)、心理學(xué)、臨床醫(yī)學(xué)、數(shù)學(xué)等各個(gè)學(xué)科。

α 波是自發(fā)腦電的主要成分之一，對(duì)應(yīng)了大腦視覺(jué)皮層的閑散節(jié)律，頻率范圍為 8 ～12 Hz，幅值一般小于100 μV。大腦皮層各個(gè)位置都會(huì)產(chǎn)生 α 波，但是，枕部的α 波阻斷現(xiàn)象最為明顯［8，9］。不同受試者的 α 波幅值雖然具有明顯個(gè)體差異，但同一受試者在睜眼、閉眼2 種狀態(tài)下的α 波幅值差異顯著，通過(guò)分析腦電數(shù)據(jù)容易區(qū)分出2 種不同狀態(tài)。

按照信號(hào)控制方式的不同，BCI 系統(tǒng)分為同步(synchronous)模式和異步(asynchronous)模式［10］，同步 BCI 系統(tǒng)要求用戶(hù)必須在預(yù)定的時(shí)間進(jìn)入一種特定的精神狀態(tài)，異步BCI 系統(tǒng)需要不斷分析腦電信號(hào)，用戶(hù)可以自由啟動(dòng)某種思維任務(wù)來(lái)完成特定的控制。本文建立的實(shí)時(shí)光標(biāo)控制系統(tǒng)屬于異步模式的BCI 系統(tǒng)。

1 系統(tǒng)設(shè)計(jì)

圖1 為BCI 系統(tǒng)結(jié)構(gòu)圖，包括傳感器系統(tǒng)和信號(hào)處理系統(tǒng)。傳感器系統(tǒng)從受試者頭皮表面采集微弱的腦電信號(hào)，經(jīng)過(guò)信號(hào)放大、A/D 轉(zhuǎn)換后，發(fā)送給信號(hào)處理系統(tǒng)處理;信號(hào)處理系統(tǒng)將接收到的EEG 信號(hào)轉(zhuǎn)換成反映受試者意圖的命令控制光標(biāo)移動(dòng)。

圖1 BCI 系統(tǒng)結(jié)構(gòu)設(shè)計(jì)Fig 1 Structural design of BCI system

1.1 傳感器系統(tǒng)

傳感器系統(tǒng)采用Ag/AgCl 電極獲取腦電信號(hào)，腦電信號(hào)采集電路對(duì)其濾除干擾信號(hào)、程控放大、消除電平漂移、帶通濾波、A/D 轉(zhuǎn)換，然后發(fā)送給信號(hào)處理系統(tǒng)。采集電路采用C8051F410 內(nèi)部的12 位ADC 采樣，由于該ADC 分辨率不高而眼電信號(hào)的幅度在μV 量級(jí)，且腦電信號(hào)中含有很強(qiáng)的干擾，無(wú)法直接轉(zhuǎn)換為數(shù)字信號(hào)，因此，必須對(duì)其進(jìn)行放大和去干擾。對(duì)腦電信號(hào)進(jìn)行采樣，放大倍數(shù)需設(shè)定在1000 ～32000 倍，需通過(guò)兩級(jí)或兩級(jí)以上的電路來(lái)完成，以保證信號(hào)的線性放大。由于采用多級(jí)放大電路，腦電信號(hào)的前置級(jí)電路的性能直接影響到整個(gè)系統(tǒng)的特性，尤其是信噪比。

圖2 為單通道的腦電信號(hào)采集電路框圖，系統(tǒng)分三級(jí)實(shí)現(xiàn)系統(tǒng)的放大增益，前置級(jí)放大電路采用AD8221 儀表放大器，可以去除極化電壓和共模信號(hào)，其放大倍數(shù)為10.15 倍。前置級(jí)放大輸出的信號(hào)仍然存在基線不穩(wěn)的情況，基線漂移嚴(yán)重時(shí)輸出會(huì)超過(guò)放大器動(dòng)態(tài)范圍，測(cè)不出信號(hào)，因此，在前置級(jí)電路后面加上去基線漂移電路，即限幅電路。經(jīng)過(guò)去基線漂移的腦電信號(hào)還必須通過(guò)帶通濾波才能被識(shí)別，帶通濾波電路同時(shí)可放大腦電信號(hào)，放大倍數(shù)為100 倍。電平遷移電路把信號(hào)的基線電平調(diào)整到1.1 V 左右，程控放大電路進(jìn)行可調(diào)增益的后級(jí)放大，放大倍數(shù)為1～32 倍，最后由 C8051F410 內(nèi)部 ADC 采樣。此外輸入的腦電信號(hào)可能出現(xiàn)幅度過(guò)大的干擾信號(hào)，存在損壞電路的風(fēng)險(xiǎn)，因此，在儀表放大電路前增加了輸入保護(hù)電路，將輸入腦電信號(hào)電壓限制在一定范圍內(nèi)，從而保證電路正常工作。

圖2 單通道的傳感器系統(tǒng)框圖Fig 2 Block diagram of single channel sensor system

1.2 信號(hào)處理系統(tǒng)

信號(hào)處理系統(tǒng)包括信號(hào)處理和實(shí)時(shí)光標(biāo)移動(dòng)控制部分，本文利用LabVIEW 語(yǔ)言實(shí)現(xiàn)。信號(hào)處理系統(tǒng)接收到EEG 信號(hào)后對(duì)其進(jìn)行8～12 Hz 的帶通濾波，得到 O1，O2 導(dǎo)聯(lián)的α 波，濾波器采用10 階 Butterworth 無(wú)限長(zhǎng)沖擊響應(yīng)(IIR)濾波器，然后對(duì)這導(dǎo)聯(lián)的α 波能量進(jìn)行600 ms 滑窗平均處理。設(shè)Xij為第i 導(dǎo)的第j 個(gè)采樣點(diǎn)腦電數(shù)據(jù)，則N 個(gè)采樣點(diǎn)對(duì)應(yīng)第i 導(dǎo)腦電數(shù)據(jù)的平均能量為

本文選取采樣率為 1 000 Hz，N =600 時(shí) O1，O2 兩導(dǎo)聯(lián)的α 波平均能量P1和P2作為特征。最后根據(jù)Fisher 線性判別分析(Fisher discriminant analysis，F(xiàn)DA)進(jìn)行模式分類(lèi)，其原理是試圖通過(guò)訓(xùn)練數(shù)據(jù)找到一組最佳的投影方向w，在這些投影方向投影可以將訓(xùn)練集中不同類(lèi)別的樣本最大程度地區(qū)分開(kāi)來(lái)［11］。定義Fisher 的準(zhǔn)則函數(shù)

用lagrange 乘子法求解w 的最優(yōu)解w*

式中 Sw為樣本類(lèi)內(nèi)離散度，Sb為樣本類(lèi)間離散度，m1，m2為2 類(lèi)樣本的均值。

圖3 為Fisher 法二分類(lèi)示意圖，其中實(shí)心圓與空心圓代表二維空間的2 類(lèi)不同樣本。從圖中可以看出原樣本無(wú)論在x 軸或是y 軸上的投影都是混雜的，因此，直接將它們投影到x 軸或是y 軸上都是不好分類(lèi)的，但是投影到方向w上就很容易區(qū)分。

圖3 Fisher 法分類(lèi)示意圖Fig 3 Diagram of classification by Fisher method

2 實(shí)驗(yàn)設(shè)計(jì)

受試者通過(guò)調(diào)節(jié)閉眼時(shí)間的長(zhǎng)短實(shí)現(xiàn)控制光標(biāo)向4 個(gè)方向移動(dòng)。實(shí)驗(yàn)方法如下:在LabVIEW 平臺(tái)下，根據(jù)式(1)得到采樣率為1 000，N =600 時(shí)受試者頭部O1，O2 導(dǎo)聯(lián)對(duì)應(yīng)的α 波平均能量P1和P2，通過(guò)線性判別分析的方法判定受試者當(dāng)前睜眼閉眼狀態(tài)并輸出控制信號(hào)Cr，當(dāng)Cr =1時(shí)系統(tǒng)發(fā)出一聲提示音，如受試者保持閉眼狀態(tài)，系統(tǒng)每隔600 ms 發(fā)出一次提示音，當(dāng)受試者僅聽(tīng)到1 聲提示音后睜眼，則表示選擇方向‘上’;聽(tīng)到2 聲提示音后睜眼，則表示選擇方向‘下’;聽(tīng)到3 聲提示音后睜眼，則表示選擇方向‘左’;聽(tīng)到4 聲提示音后睜眼，則表示選擇方向‘右’。輸出的控制信號(hào)與所對(duì)應(yīng)的命令如表1 所示。

表1 控制命令表Tab 1 Control command table

正常人閉眼后α 波能量上升輸出Cr =1 需要2.4s，保持閉眼狀態(tài)輸出上下左右4 個(gè)方向分別需要1.2，1.8，2.4，3.0 s，所以，系統(tǒng)輸出一個(gè)控制命令平均需要4.5 s。

參加實(shí)驗(yàn)的10 位受試者均為在校研究生，均未經(jīng)過(guò)先期訓(xùn)練。具體實(shí)驗(yàn)如下:受試者控制光標(biāo)分別向上下左右4 個(gè)方向移動(dòng)5 步，統(tǒng)計(jì)每位受試者完成任務(wù)所需時(shí)間和出錯(cuò)次數(shù)。如果光標(biāo)未按受試者目標(biāo)方向移動(dòng)，則記為一個(gè)錯(cuò)誤。在相同實(shí)驗(yàn)條件下，首先進(jìn)行1 組測(cè)試實(shí)驗(yàn)，以確定最佳投影方向w，然后每位受試者重復(fù)完成5 組實(shí)驗(yàn)，作為系統(tǒng)實(shí)驗(yàn)用以檢驗(yàn)系統(tǒng)的正確率。

3 實(shí)驗(yàn)結(jié)果分析

表2 為系統(tǒng)實(shí)驗(yàn)記錄，從表中可以看出:受試者通過(guò)幾次實(shí)驗(yàn)訓(xùn)練即可有效提高操作正確率。分析出現(xiàn)誤操作的主要原因有以下幾種:1)受試者閉眼期間α 波幅值下降導(dǎo)致錯(cuò)誤輸出‘0’;2)受試者睜眼后α 波幅值下降延遲而導(dǎo)致多輸出‘1’;3)受試者睜眼狀態(tài)下α 波幅值出現(xiàn)抖動(dòng)上升而導(dǎo)致錯(cuò)誤輸出‘1’。

表2 系統(tǒng)實(shí)驗(yàn)記錄Tab 2 Records of system experiments

從實(shí)驗(yàn)數(shù)據(jù)看，不同受試者在睜眼和閉眼時(shí)α 波能量差異明顯，進(jìn)一步驗(yàn)證了α 波阻斷現(xiàn)象的普適性。通過(guò)實(shí)驗(yàn)數(shù)據(jù)統(tǒng)計(jì)和分析可以得出結(jié)論:

1)最佳投影方向w 的設(shè)定因受試者不同而異，系統(tǒng)實(shí)驗(yàn)前需進(jìn)行測(cè)試實(shí)驗(yàn)，用以確定最佳投影方向w。

2)α 波平均處理時(shí)間常數(shù)的設(shè)定直接影響了系統(tǒng)的抗干擾能力和系統(tǒng)響應(yīng)速度，由此影響了實(shí)驗(yàn)的正確率。將濾波后的α 波信號(hào)進(jìn)行平均處理，適當(dāng)?shù)钠骄幚頃r(shí)間可以有效提高系統(tǒng)抗干擾能力，如眨眼，打哈欠等;過(guò)大的平均處理時(shí)間會(huì)降低系統(tǒng)的靈敏度;平均處理時(shí)間同時(shí)也是提示音間隔時(shí)間，需考慮受試者反應(yīng)速度，通過(guò)統(tǒng)計(jì)分析多次實(shí)驗(yàn)數(shù)據(jù)將其設(shè)定為600 ms。

文獻(xiàn)［12］建立的基于腦電α 波的BCI 控制系統(tǒng)需要受試者注視有5 個(gè)指示燈(代表5 種不同操作命令)的控制面板，指示燈按順時(shí)針?lè)较蛞来窝h(huán)閃爍點(diǎn)亮并持續(xù)一段時(shí)間，如果某個(gè)指示燈點(diǎn)亮?xí)r受試者閉眼產(chǎn)生α 波，則選中該指示燈代表的操作命令，其工作方式為同步模式，需要受試者注視有指示燈的控制面板，該系統(tǒng)輸出一個(gè)控制命令需要9 s;文獻(xiàn)［10］建立的基于α 波的異步BCI 系統(tǒng)輸出一個(gè)控制命令需要45 s;本文建立的實(shí)時(shí)光標(biāo)控制系統(tǒng)輸出一個(gè)控制命令僅需3.6～5.4 s，操作簡(jiǎn)單，響應(yīng)速度快，正確率高，實(shí)用價(jià)值高，與現(xiàn)有的基于腦電α 波的BCI 系統(tǒng)相比，具有更好的實(shí)時(shí)性，見(jiàn)表3。

表3 基于腦電α 波的BCI 系統(tǒng)對(duì)比Tab 3 Comparision of BCI system based on EEG α wave

4 結(jié)束語(yǔ)

本文利用在腦神經(jīng)功能正常的人中普遍存在的腦電α 波阻斷現(xiàn)象，設(shè)計(jì)了一種基于傳感器技術(shù)的實(shí)時(shí)光標(biāo)控制BCI 系統(tǒng)。10 名受試者的操作實(shí)驗(yàn)表明:未經(jīng)先期訓(xùn)練的受試者可以利用該系統(tǒng)穩(wěn)定可靠地實(shí)現(xiàn)既定實(shí)驗(yàn)任務(wù)，驗(yàn)證了此系統(tǒng)的易操作性、較高的正確率和較快的響應(yīng)速度。同現(xiàn)有的基于腦電α 波的BCI 系統(tǒng)相比，具有更好的實(shí)用價(jià)值。本文的研究工作為進(jìn)一步開(kāi)發(fā)便攜式EEG 采集設(shè)備和實(shí)時(shí)控制的BCI 系統(tǒng)提供了技術(shù)基礎(chǔ)。

［1］ Wolpaw J R，McFarland D J，Neat G W，et al.An EEG－based brain－computer interface for cursor control［J］.Electroencephalography and Clinical Neurophysiology，1991，78(3):252 －259.

［2］ Wolpaw J R，Birbaumer N，McFarland D J，et al.Brain－computer interfaces for communication and control［J］.Clinical Neurophysiology，2002，113(6):767 － 791.

［3］ Wolpaw J R，McFarland D J，Vaughan T M.Brain－computer interface research at the Wadsworth Center［J］.IEEE Transactions on Rehabilitation Engineering，2000，8(2):222 －226.

［4］ 韓 雪，李醒飛，鐘 瑩.腦-機(jī)接口用微傳感器設(shè)計(jì)與制造方法研究［J］.傳感器與微系統(tǒng)，2009，28(1):52 －55.

［5］ Pfurtscheller G，Muller－Putz G R，Schlogl A，et al.15 years of BCI research at Graz University of Technology:Current projects［J］.IEEE Transactions on Neural Systems and Rehabilitation Engineering，2006，14(2):205 －210.

［6］ Pfurtscheller G，Kalcher J，Neuper C，et al.On－line EEG classification during externally－paced hand movements using a neural network－based classifier［J］.Electroencephalography and Clinical Neurophysiology，1996，99(5):416 －425.

［7］ Leeb R，Pfurtscheller G.Walking through a virtual city by thought［C］∥2004 The 26th Annual International Conference of IEEE Engineering in Medicine and Biology Society，IEMBS’04，2004:4503 －4506.

［8］ Dewan E M.Occipital alpha rhythm eye position and lens accommodation［J］.Nature，1967，214(5092):975 － 977.

［9］ 黃遠(yuǎn)桂，吳聲伶.臨床腦電圖學(xué)［M］.北京:人民衛(wèi)生出版社，1998.

［10］ 趙海濱，王 宏，李春勝.基于Alpha 波的異步腦-機(jī)接口系統(tǒng)［J］.東北大學(xué)學(xué)報(bào):自然科學(xué)版，2009(10):1377 －1380.

［11］ 邊肇祺，張學(xué)工.模式識(shí)別［M］.北京:清華大學(xué)出版社，1988.

［12］ 萬(wàn)柏坤，綦宏志，趙 麗.基于腦電 Alpha 波的腦-機(jī)接口控制實(shí)驗(yàn)［J］.天津大學(xué)學(xué)報(bào)，2006(8):978 －984.