張小棟,郭晉,李睿,陸竹風(fēng)

(1.西安交通大學(xué)現(xiàn)代設(shè)計及轉(zhuǎn)子軸承系統(tǒng)教育部重點實驗室,710049,西安;2.西安交通大學(xué)機(jī)械工程學(xué)院,710049,西安;3.西安測繪總站裝備技術(shù)室,710054,西安)

?

表情驅(qū)動下腦電信號的建模仿真及分類識別

張小棟1,2,郭晉2,3,李睿1,2,陸竹風(fēng)2

(1.西安交通大學(xué)現(xiàn)代設(shè)計及轉(zhuǎn)子軸承系統(tǒng)教育部重點實驗室,710049,西安;2.西安交通大學(xué)機(jī)械工程學(xué)院,710049,西安;3.西安測繪總站裝備技術(shù)室,710054,西安)

針對自發(fā)型腦電信號識別率低、個體差異度大等問題,提出了一種新的基于表情驅(qū)動腦電信號的腦機(jī)接口方式,并進(jìn)行了建模仿真及實驗驗證。利用神經(jīng)元集群模型對表情驅(qū)動腦電信號進(jìn)行機(jī)理建模與仿真分析,得到自發(fā)表情的相關(guān)腦區(qū)及表情驅(qū)動腦電信號的頻率分布特性;提出了一種基于小波變換和人工神經(jīng)網(wǎng)絡(luò)模型映射的表情驅(qū)動下腦電信號分析識別方法,有效提高了表情驅(qū)動下腦電信號的識別率。從神經(jīng)生理學(xué)角度驗證了表情驅(qū)動腦電信號的特征來源是受大腦前額葉皮層和邊緣系統(tǒng)相互協(xié)調(diào)共同控制的,并通過實驗驗證了所提腦電信號分類識別算法切實有效,其最高分類準(zhǔn)確率可達(dá)85%。

表情驅(qū)動;腦電信號;機(jī)理建模;小波變換;神經(jīng)網(wǎng)絡(luò)

自1924年第一次在人腦中檢測到有規(guī)律的電活動開始,腦-機(jī)接口(brain-computer interface, BCI)技術(shù)成為了各國學(xué)者的研究熱點。然而,目前關(guān)于腦-機(jī)接口的應(yīng)用研究,都不能達(dá)到快速、靈活、準(zhǔn)確的控制要求,只能實現(xiàn)一些簡單的功能和進(jìn)行一些基礎(chǔ)性的研究工作;對腦電信號的識別利用,都只是從信號分析的層面入手,很少有人從大腦的生理機(jī)能角度出發(fā),深入研究大腦不同區(qū)域在活動狀態(tài)下的機(jī)理特性,通過構(gòu)建腦電信號(electroencephalogram, EEG)的機(jī)理模型提高其分類識別的效果[1]。

本文提出了一種基于表情驅(qū)動下的腦電信號分類識別方法。在表情驅(qū)動下的腦電信號采集過程中,受試者完成不同的表情動作,例如提眉、皺眉、撇嘴等,即可采集相應(yīng)動作下的腦電信號。表情動作對受試者要求低,而且在人們的日常生活中經(jīng)常利用不同表情傳遞信息,因此具有非常廣泛的適用性。

本文工作是針對表情驅(qū)動下腦-機(jī)接口范式的一項基礎(chǔ)性探索研究,即通過對表情驅(qū)動下的腦電信號進(jìn)行機(jī)理建模與仿真分析,得知表情驅(qū)動下腦電信號產(chǎn)生的相關(guān)大腦區(qū)域及頻率分布特性。以仿真分析結(jié)果為指導(dǎo),采集相應(yīng)大腦區(qū)域的腦電信號,同時利用小波特征提取方法和人工神經(jīng)網(wǎng)絡(luò)識別算法實現(xiàn)對腦電信號特征集中的頻率段信息進(jìn)行分類識別。

1 表情驅(qū)動下腦電信號的建模仿真

1.1 相關(guān)機(jī)理分析

2010年P(guān)rice等通過對情緒及表情相關(guān)神經(jīng)環(huán)路的研究認(rèn)為,表情與紋狀體和海馬回等相關(guān)區(qū)域為核心的邊緣系統(tǒng)有關(guān),隨著神經(jīng)解剖學(xué)研究的進(jìn)一步探索發(fā)現(xiàn),表情的表達(dá)與大腦邊緣系統(tǒng)和前額葉皮質(zhì)及其他一些環(huán)路相關(guān),尤其與內(nèi)側(cè)前額皮質(zhì)、杏仁核、內(nèi)側(cè)紋狀體、內(nèi)側(cè)丘腦、下丘腦及腦干的一些部分關(guān)系最緊密[2],如圖1所示。

圖1 大腦皮質(zhì)層及邊緣系統(tǒng)結(jié)構(gòu)圖

另外,研究者通過分析正常腦電信號頻率分類得出,θ波與α波是一種具有高度暗示性的狀態(tài)。當(dāng)人做出不同表情動作時,通常暗示大腦處于某種精神狀態(tài),例如提眉時瞳孔放大,更多視覺信息進(jìn)入眼部供大腦快速處理分析,皺眉時人的視覺往往會將焦點聚集到某處,使人的注意力集中。以上研究表明,表情的加工與產(chǎn)生主要是由大腦內(nèi)一個分布式神經(jīng)網(wǎng)絡(luò)所控制,其中前額葉皮質(zhì)和邊緣系統(tǒng)(杏仁核、海馬回、紋狀體等)主要負(fù)責(zé)協(xié)調(diào)和控制表情,而θ波和α波則與表情動作產(chǎn)生時的暗示性狀態(tài)有關(guān)。因此,可以對大腦前額葉皮質(zhì)層和邊緣系統(tǒng)這兩個區(qū)域以及θ波和α波這兩個頻段進(jìn)行建模仿真分析,以驗證相關(guān)神經(jīng)生理學(xué)研究的推斷,并對表情驅(qū)動下腦電信號的特征提取及分類識別研究提供依據(jù)。

1.2 表情驅(qū)動下腦電信號建模仿真

目前人們主要研究的腦神經(jīng)模型有兩類,一類是從微觀角度出發(fā),關(guān)注每個神經(jīng)元活動的模型,例如Hodgkin-Huxley模型[3]、Chay模型[4]、Traub模型[5]等。顯然,僅通過分析單個神經(jīng)元的活動對理解大腦的復(fù)雜行為是遠(yuǎn)遠(yuǎn)不夠的。另一類模型是從宏觀角度出發(fā)的集總參數(shù)模型,例如神經(jīng)元集群模型(neural mass model,NMM)[6],其思想是“平均區(qū)域近似”原理,就是大腦內(nèi)部具有相似功能的同一區(qū)域具有相同的集總參數(shù)。NMM具有非常明顯的生理意義,神經(jīng)元群之間的耦合作用可以反映不同腦區(qū)的相互作用,可以從宏觀角度上仿真大尺度的腦電信號。

1.2.1 基本神經(jīng)元集群的模型建立 如圖2a所示,一個基本神經(jīng)元群由興奮性和抑制性兩個子細(xì)胞群構(gòu)成,它們相互作用,產(chǎn)生神經(jīng)振蕩。其中:粗框內(nèi)的部分表示子群1,主要由錐體細(xì)胞組成,負(fù)責(zé)神經(jīng)沖動的傳出;粗框外的部分表示子群2,主要包括其他非錐體細(xì)胞,負(fù)責(zé)神經(jīng)沖動的增益和反饋,例如藍(lán)細(xì)胞和星形細(xì)胞。

興奮性和抑制性線性變換函數(shù)的作用是將突觸前膜動作電位平均脈沖密度轉(zhuǎn)變成突觸后膜的興奮或抑制膜電壓,其單位沖激響應(yīng)分別為

he(t)=u(t)He/τete-t/τe

(1)

hi(t)=u(t)Hi/τite-t/τi

(2)

(a)基本神經(jīng)元集群模型框圖

(b)多動態(tài)基本神經(jīng)元集群模型

(c)雙區(qū)域耦合神經(jīng)元集群模型圖2 神經(jīng)元集群模型

式中:u(t)為Heaviside函數(shù);參數(shù)He和Hi為興奮性和抑制性神經(jīng)元突觸增益;τe和τi分別為興奮性和抑制性神經(jīng)元樹突平均延時常數(shù)與細(xì)胞膜平均延時常數(shù)的和。線性函數(shù)he(t)和hi(t)可由如下形式的一階常微分方程表示

(3)

式中:G代表興奮和抑制條件下的增益參數(shù)He和Hi;g代表興奮和抑制條件下的延時參數(shù)τe、τi;x(t)為線性函數(shù)的輸入。對于子群1,x(t)=S(y1(t)-y2(t)),對于子群2中的興奮線性函數(shù)he(t),x(t)=p(t)+C2S(C1y0(t)),對于子群2中的抑制線性函數(shù)hi(t),x(t)=C4S(C3y0(t))。其中,錐體細(xì)胞群和其他非錐體細(xì)胞群之間的平均連接常數(shù)用C1、C2、C3、C4表示。z(t)為線性函數(shù)的輸出,在圖2a中,分別由y0、y1、y2表示。

靜態(tài)非線性函數(shù)的作用是將平均膜電壓轉(zhuǎn)換成動作電位的平均發(fā)放率,并作為線性函數(shù)he(t)和hi(t)的輸入

S(v)=2e0/(1+er(v0-v))

(4)

式中:2e0為最大發(fā)放率;v0為相對于發(fā)放率e0的突觸后膜電位;r表示S(v)函數(shù)的彎曲程度;v表示突觸前膜電壓。

子群1接收來自子群2的興奮和抑制的反饋信號。子群1與子群2的建模思想相同,但子群2只有興奮性信號輸入,所以只需要一個興奮性線性函數(shù)。模型仿真時來自大腦其他不定區(qū)域的信號由高斯信號p(t)表示。

綜上所述,基本神經(jīng)元集群模型可由以下微分方程組表示

(5)

式中:y3、y4、y5分別為線性函數(shù)輸出量y0、y1、y2的一階導(dǎo)數(shù)?；旧窠?jīng)元集群模型的輸出為興奮和抑制性細(xì)胞膜后電位的差,即y(t)=y1(t)-y2(t)。

按照不同頻率特性,腦電波可以分為δ、θ、α、β和γ波5個頻段?？紤]到δ波僅存在于人的深度睡眠中,γ波段屬于腦電波中的高頻段,其作用尚在研究中,故在多動態(tài)基礎(chǔ)神經(jīng)元集群模型中選擇腦電波的θ、α和β波頻段建立模型。表1給出了產(chǎn)生θ、α和β波特征波形的參數(shù)取值[7]。

通過調(diào)整線性轉(zhuǎn)換函數(shù)h(t)中的參數(shù)He、τe和Hi、τi,即可調(diào)整基本神經(jīng)元集群模型內(nèi)部的動態(tài)平衡,從而仿真多種節(jié)律腦電信號。圖3為基本神經(jīng)元集群模型模擬不同頻率的窄帶信號,通過調(diào)節(jié)興奮和抑制性的突觸增益和延時時間常數(shù),可以模擬大腦產(chǎn)生近似的θ、α和β波。

1.2.2 多動態(tài)神經(jīng)元集群的模型建立 利用基本神經(jīng)元集群模型模擬單一的節(jié)律波只能反映腦電信號的組成結(jié)構(gòu),還不能達(dá)到模擬真實采集得到的腦電信號的目標(biāo),因此還需要建立多動態(tài)基本神經(jīng)元集群模型,其結(jié)構(gòu)如圖2b所示。

表1 模型生理學(xué)參數(shù)的典型值

(a)θ波仿真節(jié)律波與歸一化頻譜

(b)α波仿真節(jié)律波與歸一化頻譜

(c)β波仿真節(jié)律波與歸一化頻譜圖3 仿真θ、α和β節(jié)律波及歸一化頻譜

1.2.3 表情驅(qū)動腦電信號的產(chǎn)生機(jī)理與建模仿真 由前述表情的加工和產(chǎn)生機(jī)理可知,大腦用于控制表情的區(qū)域集中位于前額葉皮質(zhì)層和邊緣系統(tǒng),因此本文利用神經(jīng)集群模型的多動態(tài)特性對這兩個區(qū)域進(jìn)行建模仿真,并試圖驗證不同表情下腦電信號包含的不同頻帶的分布特性,如圖2c所示。

令區(qū)域1和區(qū)域2為兩個多動態(tài)神經(jīng)元集群,分別表示前額葉皮層和邊緣系統(tǒng),每個區(qū)域包括興奮性和抑制性細(xì)胞子群,每個子群由3個線性轉(zhuǎn)換函數(shù)h(t)分別表示θ、α和β節(jié)律波,區(qū)域1與區(qū)域2之間的耦合強(qiáng)度分別用qkj表示,其中k,j∈(1,2)且k≠j,同時每個區(qū)域還接收來自其他不確定區(qū)域的信號,用p(t)表示。

兩個區(qū)域間傳遞的耦合信號設(shè)為Gkj,可表示為

(6)

式(6)中R(x)=x-mean(x)為求均值,用于保持各個區(qū)域信號之間的平衡。雙區(qū)域耦合神經(jīng)元集群模型的微分方程如下式所示

(7)

由于缺少相關(guān)的神經(jīng)解剖學(xué)實驗經(jīng)驗及研究資料,本文通過實驗反推法確定權(quán)重系數(shù)W。根據(jù)表情驅(qū)動下腦電信號的產(chǎn)生機(jī)理,采集受試者在不同表情驅(qū)動下的前額葉皮質(zhì)區(qū)域FC5(左)、FC6(右)及邊緣系統(tǒng)區(qū)域的F7(左)和F8(右)4通道腦電信號并進(jìn)行快速傅里葉變換,得到其功率譜密度函數(shù),其后將功率譜密度函數(shù)分成5個子頻帶,分別為0～4 Hz、4～8 Hz、8～16 Hz,16～32 Hz,32～64 Hz。

對腦電信號的5個子頻帶求對數(shù)頻帶能量Ei,即可得每個子頻帶在總頻帶中的能量分布

(8)

(9)

表2為所得的4～8 Hz、8～16 Hz及16～32 Hz頻帶能量分布比值。設(shè)定提眉動作仿真模型中W=[0.13,0.60,0.19],右撇嘴仿真模型中W=[0.19,0.40,0.26],即可得到真實表情動作下腦電信號的擬合仿真信號,如圖4所示。

表2 表情刺激下腦電信號不同頻帶分布

(a)提眉仿真信號和歸一化頻譜圖

(b)提眉真實信號和歸一化頻譜圖

(c)右撇嘴仿真信號和歸一化頻譜圖

(d)右撇嘴真實信號和歸一化頻譜圖圖4 F7通道下提眉與右撇嘴表情驅(qū)動產(chǎn)生的腦電信號與相應(yīng)的仿真結(jié)果對比圖

由圖4中的仿真結(jié)果可知,提眉動作下的腦電信號能量主要集中在低頻10 Hz以內(nèi),而右撇嘴動作下的腦電信號能量分布除了以低頻部分為主以外,還夾雜著一些中頻和高頻成分,這說明提眉動作和右撇嘴動作下的腦電信號頻域能量主要集中在低頻部分即θ波,以及一部分中頻部分即α波,而δ、β和γ波段中的信息較少。通過前額葉皮層區(qū)和邊緣系統(tǒng)的耦合建模,說明了人對表情的控制是通過這兩個區(qū)域相互協(xié)調(diào)來進(jìn)行的。

仿真結(jié)果驗證了之前關(guān)于表情驅(qū)動下腦電信號相關(guān)神經(jīng)生理學(xué)的機(jī)理分析結(jié)論,因此通過綜合考慮,本文將在前額葉皮層處和大腦邊緣系統(tǒng)附近采集腦電信號,并對表情驅(qū)動下腦電信號的θ和α波兩個子頻帶進(jìn)行特征提取和模式識別工作。

2 基于小波變換的表情驅(qū)動腦電信號多維特征提取研究

由于表情驅(qū)動下的腦電信號是一種非平穩(wěn)、非線性、時變的信號,單純使用時域或頻域分析方法都無法完整地展示腦電信號的特征,而小波分析作為應(yīng)用最廣泛的時頻域分析方法,具有多尺度、多分辨率分析的優(yōu)勢[9]。本文采用基于D-6小波基的小波變換Mallat算法提取表情驅(qū)動下腦電信號的特征。

本文研究中腦電信號的采樣頻率選為128 Hz,有效頻率范圍為0～64 Hz。進(jìn)行4級小波分解后,各分解頻率范圍分別與δ、θ、α、β和γ波的頻率范圍基本相同。在實際應(yīng)用中,小波分解后構(gòu)成的特征量往往取各頻段的小波系數(shù)、小波系數(shù)的統(tǒng)計信息或小波系數(shù)的模極大值等。本文研究中采用小波系數(shù)模均值作為特征量。

根據(jù)前文中對表情驅(qū)動下腦電信號建模仿真可知,表情驅(qū)動下腦電信號的能量分布主要集中于θ與α波內(nèi)。因此,對前額葉皮層附近和邊緣系統(tǒng)附近的F7、FC5、FC6和F8這4個通道的θ和α波兩個頻帶的數(shù)據(jù)進(jìn)行特征提取,由模均值組成1×8的特征向量如下

Xk=[xk1xk2xk3xk4xk5xk6xk8xk8]

(10)

式中:xk1～xk4分別為4個通道腦電信號中的θ波小波系數(shù)模均值;xk5～xk8分別為4個通道腦電信號中的α波小波系數(shù)模均值。

將特征向量Xk歸一化處理后,輸入神經(jīng)網(wǎng)絡(luò)系統(tǒng)中進(jìn)行模式識別。本文采用的歸一化方法為

x′=(x-xmin)/(xmax-xmin)

(11)

那么,經(jīng)過歸一化后的小波特征矩陣為

Tk=[tk1tk2tk3tk4tk5tk6tk7tk8]

(12)

本文在實驗中每種表情選取10組特征數(shù)據(jù),組成8×40的特征矩陣,作為神經(jīng)網(wǎng)絡(luò)模式識別系統(tǒng)的輸入,如圖5所示。

(a)提眉動作

(b)皺眉動作

(c)左撇嘴動作

(d)右撇嘴動作圖5 4種表情特征的小波系數(shù)模均值

由圖5可以看出,每種表情在重復(fù)實驗的情況下特征變化比較平穩(wěn),4種表情的10次樣本重復(fù)效果比較理想,盡管由于實驗過程中受試者的疲勞反應(yīng)、精神狀態(tài)不同導(dǎo)致每次動作下的特征值略有不同,但是每次樣本間的差異并不大,其波動范圍均在0～0.1之間,不會對模式識別結(jié)果造成太大影響。因此可以說明,通過對F7、FC5、FC6和F8這4個通道的θ和α波兩個頻帶進(jìn)行4級小波變換,提取的表情驅(qū)動下的腦電信號特征值具有一定的可分性,可以作為模式識別網(wǎng)絡(luò)的輸入量進(jìn)行分類。

3 基于人工神經(jīng)網(wǎng)絡(luò)的表情驅(qū)動腦電信號模式識別方法研究

人工神經(jīng)網(wǎng)絡(luò)算法具有較強(qiáng)的容錯能力、允許對問題的了解較少、可實現(xiàn)特征空間中較復(fù)雜的劃分等優(yōu)勢。誤差后向傳播神經(jīng)網(wǎng)絡(luò)(back-propagation neural network,BPNN)是目前人工神經(jīng)網(wǎng)絡(luò)中研究最深入、應(yīng)用最廣泛的一種模型[10]。因此,本文將研究其在表情驅(qū)動腦電信號模式識別中的具體應(yīng)用。

在本研究中,從小波變換得出的特征量為8個元素,因此輸入層為8個節(jié)點。目標(biāo)為實現(xiàn)4種表情動作(提眉、皺眉、左撇嘴、右撇嘴)腦電信號的模式識別,因此在輸出層只需y1和y2這2個神經(jīng)元即可滿足需求。經(jīng)過訓(xùn)練期望得到的輸出層目標(biāo)結(jié)果與相應(yīng)的表情動作見表3。

表3 輸出層目標(biāo)結(jié)果及相應(yīng)的表情動作

一般而言,一個3層的BP網(wǎng)絡(luò)可以完成任意從n維到m維的映射。本研究網(wǎng)絡(luò)屬于中小型結(jié)構(gòu),所以優(yōu)先選用一個隱含層。隱節(jié)點的個數(shù)選擇在理論上還沒有一個明確的規(guī)定,目前大多以經(jīng)驗為主,本文參考公式為

h=(N+M)0.5+a

(13)

式中:M為輸出單元數(shù);N為輸入單元數(shù);a為1～10之間的常數(shù)。在保證高穩(wěn)定識別率的前提下,遵循盡量采用較少的隱節(jié)點的原則,最終確定采用8-12-2的BPNN結(jié)構(gòu)用于表情驅(qū)動EEG的模式識別,其結(jié)果如圖6所示,其中圖6a、圖6b分別為神經(jīng)網(wǎng)絡(luò)輸出層y1、y2的網(wǎng)絡(luò)輸出與目標(biāo)結(jié)果的對比結(jié)果。

(a)y1的識別結(jié)果

(b)y2的識別結(jié)果圖6 BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練結(jié)果

4 實驗研究與討論

本實驗采用由美國加州舊金山的神經(jīng)科技公司Emotiv Systems開發(fā)研制的Emotiv-Epoc腦電采集系統(tǒng)作為腦電信號采集設(shè)備。采樣頻率為128 Hz,通道數(shù)為16。Emotive腦電帽與PC終端通過藍(lán)牙傳輸數(shù)據(jù),表情驅(qū)動腦電信號分類實驗平臺如圖7所示。

圖7 表情驅(qū)動腦電信號分類實驗場景圖

受試者為3名在校男生,年齡均在22～30歲,均無神經(jīng)疾病史。本實驗研究對象為4種不同的面部表情,即提眉、皺眉、左撇嘴、右撇嘴。受試者在實驗過程中要完成這4種表情動作,如圖8所示。

圖8 4種不同表情動作

實驗過程中,室內(nèi)光線勿強(qiáng),讓受試者背靠椅背,身體盡量放松,安靜、不眨眼、不轉(zhuǎn)動眼球、勿入睡,只做規(guī)定的4種表情動作,身體其他部位盡量不活動。實驗時間不超過30 min。實驗中應(yīng)確保各種表情動作具有良好的一致性和可重復(fù)性。

為了分析表情驅(qū)動腦電信號對不同受試者的差異性,本文基于上述表情驅(qū)動腦電信號的產(chǎn)生機(jī)理分析,選取大腦皮層中前額葉皮質(zhì)區(qū)域和邊緣系統(tǒng)區(qū)域的FC5、FC6、F7和F8這4個通道的腦電信號作為分析對象。實驗中要求受試者每種表情重復(fù)5次為一組,4種表情分別做20組,即每名受試者采集80組EEG數(shù)據(jù)。將80組數(shù)據(jù)中的60組作為訓(xùn)練樣本,20組作為測試樣本。表4為3名受試者在3次實驗中的單次最高離線識別率,其中最高離線識別率為受試者3的85%。

表4 不同表情的單次最高離線識別率

通過實驗結(jié)果可知,基于表情驅(qū)動下的腦電信號的特征提取與模式識別方法是切實可行的。但是,由于腦電信號的個體差異較大,因此有必要對不同受試者與智能假肢的人機(jī)智能系統(tǒng)的協(xié)同策略與方法展開研究,探討人與假肢在該系統(tǒng)中各自的職能與其交互工作過程的流程原理,以提高該人機(jī)系統(tǒng)的執(zhí)行效率、安全性和魯棒性,進(jìn)而提高整個腦控假手系統(tǒng)的精密控制性能。

5 結(jié) 論

綜上所述,本文通過表情驅(qū)動下腦電信號的建模仿真與分類識別研究,可以得出以下3個主要結(jié)論:

(1)利用神經(jīng)集群模型建立的表情驅(qū)動腦電信號仿真模型,能夠從神經(jīng)生理學(xué)角度驗證表情驅(qū)動下腦電信號的產(chǎn)生和控制是在大腦前額葉皮層和邊緣系統(tǒng)相互協(xié)調(diào)共同作用下完成的,對分析并利用復(fù)雜腦電信號中蘊(yùn)含的特征信息提供了依據(jù);

(2)針對表情驅(qū)動腦電信號提出的小波變換特征提取方法,能從多個不同尺度空間更為全面地反映原始腦電信號的特征,從而使提取的特征量具有良好的可分性;

(3)針對表情驅(qū)動腦電信號分類識別問題提出的BP神經(jīng)網(wǎng)絡(luò)模式識別策略,能夠?qū)⒉煌砬橄碌哪X電信號進(jìn)行識別分類,最高離線識別率可達(dá)85%。

但是,目前腦電信號采集系統(tǒng)對硬件的精確度要求較高,一般在腦電信號采集過程中會產(chǎn)生少量偽差,需要在數(shù)據(jù)預(yù)處理階段進(jìn)行必要的識別并剔除,否則會在一定程度上影響識別的效果。加之腦電信號的個體差異較大,因此有必要研究適應(yīng)具體受試者的多維特征選擇算法,從而提高不同表情驅(qū)動下腦電信號模式識別系統(tǒng)的針對性和泛化性。

另外,對表情驅(qū)動腦電信號的仿真建模,僅定性地分析了腦電波中包含的頻率成分,其結(jié)論是由真實采集到的腦電信號反推而來的,能夠在一定程度上幫助人們理解表情驅(qū)動下腦電信號的產(chǎn)生機(jī)理,要想得到準(zhǔn)確的腦電信號機(jī)理模型,還需要神經(jīng)解剖學(xué)、生理學(xué)、醫(yī)學(xué)等領(lǐng)域的進(jìn)一步協(xié)同研究與發(fā)展。

[1] 張小棟, 李睿, 李耀楠. 腦控技術(shù)的研究與展望 [J]. 振動、測試與診斷, 2014, 34(2): 205-211. ZHANG Xiaodong, LI Rui, LI Yaonan. Research and forecast on brain-computer interface [J]. Journal of Vibration, Measurement and Diagnosis, 2014, 34(2): 205-211.

[2] PRICE J L, DREVETS W C. Neurocircuitry of mood disorders [J]. Neuropsychopharmacology, 2010, 35(1): 192-216.

[3] HODGKIN A L, HUXLEY A F. A quantitative description of membrane current and its application to conduction and excitation in nerve [J]. The Journal of Physiology, 1952, 117(4): 500-544.

[4] CHAY T R. Chaos in a three-variable model of an excitable cell [J]. Physica: D Nonlinear Phenomena, 1985, 16(2): 233-242.

[5] TRAUB R D, CONTRERAS D, CUNNINGHAM M O, et al. Single-column thalamocortical network model exhibiting gamma oscillations, sleep spindles, and epileptogenic bursts [J]. Journal of Neurophysiology, 2010, 93(4): 2194-2232.

[6] LOPES F H, HOEKS A, SMITS H, et al. Model of brain rhythmic activity [J]. Kybernetik, 1974, 15(1): 27-37.

[7] WENDLING F, BELLANGER J J, BARTOLOMEI F, et al. Relevance of nonlinear lumped-parameter models in the analysis of depth-EEG epileptic signals [J]. Biological Cybernetics, 2000, 83: 367-378.

[8] JANSEN B H, RIT V G. Electroencephalogram and visual evoked potential generation in a mathematical model of coupled cortical columns [J]. Biological Cybernetics, 1995, 73(4): 357-366.

[9] 王登, 苗奪謙, 王睿智. 一種新的基于小波包分解的EEG特征抽取與識別方法研究 [J]. 電子學(xué)報, 2013, 41(1): 193-198. WANG Deng, MIAO Duoqian, WANG Ruizhi. A new method of EEG classification with feature extraction based on wavelet packet decomposition [J]. Acta Electronica Sinica, 2013, 41(1): 193-198.

[10]張宇博, 舒紅平. 改進(jìn)的BP神經(jīng)網(wǎng)絡(luò)算法及其應(yīng)用 [J]. 網(wǎng)絡(luò)安全技術(shù)與應(yīng)用, 2007, 27(12): 172-175. ZHANG Yubo, SHU Hongping. Optimized BP neural networks algorithm and its application [J]. Computer Simulation, 2015, 27(12): 172-175.

[本刊相關(guān)文獻(xiàn)鏈接]

劉蓮艾,徐進(jìn),徐燕芝,等.早期血管性癡呆患者的效應(yīng)連接研究.2016,50(3):146-150.[doi:10.7652/xjtuxb201603023]

徐光華,張鋒,謝俊,等.穩(wěn)態(tài)視覺誘發(fā)電位的腦機(jī)接口范式及其信號處理方法研究.2015,49(6):1-7.[doi:10.7652/xjtuxb201506001]

鄭楊,王剛,徐進(jìn),等.大腦皮層-肌肉交互作用與注意力關(guān)系的實驗研究.2014,48(6):134-138.[doi:10.7652/xjtuxb2014 06023]

董淵哲,李曉玲,姚磊,等.面向神經(jīng)工效學(xué)的觸控顯示系統(tǒng)人機(jī)交互評價方法.2014,48(5):123-128.[doi:10.7652/xjtuxb201405022]

趙夢琰,李擴(kuò),鮑剛,等.利用小波包的腦電源定位算法仿真研究.2013,47(12):130-136.[doi:10.7652/xjtuxb201312 022]

盛恒松,田洪君,鄭崇勛,等.去除腦電信號眨眼偽差的獨立成分時域相關(guān)算法.2013,47(10):127-131.[doi:10.7652/xjtuxb201310022]

婁武濤,徐進(jìn),盛恒松,等.血管性癡呆患者腦電信號的多通道線性描述符特征研究.2011,45(4):121-125.[doi:10.7652/xjtuxb201104022]

高軍峰,王沛,鄭崇勛.基于P300和機(jī)器學(xué)習(xí)的測謊方法研究.2010,44(10):120-124.[doi:10.7652/xjtuxb201010023]

趙敏,鄭崇勛,趙春臨,等.利用Fisher判別式和事件相關(guān)電位的心理意識真實性識別.2010,44(8):132-136.[doi:10.7652/xjtuxb201008026]

高軍峰,鄭崇勛,王沛.腦電信號中肌電偽差的實時去除方法研究.2010,44(4):114-118.[doi:10.7652/xjtuxb201004024]

(編輯 武紅江)

A Simulation Model and Pattern Recognition Method of Electroencephalogram Driven by Expression

ZHANG Xiaodong1,2,GUO Jin2,3,LI Rui1,2,LU Zhufeng2

(1. Key Laboratory of Education Ministry for Modern Design and Rotor-Bearing System, Xi’an Jiaotong University,Xi’an 710049, China; 2. School of Mechanical Engineering, Xi’an Jiaotong University, Xi’an 710049, China;3. Equipment and Technique Department, Xi’an Surveying and Mapping Station, Xi’an 710054, China)

An expression driven brain-computer interface system is proposed to solve the low recognition rate and huge individual differences of spontaneous electroencephalogram (EEG) signal, and simulation and experimental verification are conducted. A neural mass model is used in mechanism modeling and simulation of expressions driven EEG signal, and the results show the related brain areas of facial expression and the frequency distribution characteristics of expression driven EEG signals. A classification and recognition method based on the wavelet transform and an artificial neural network is put forward, to improve the recognition rate of expression driven EEG signal. This research elaborates that it is the prefrontal cortex and the limbic system to cooperate on the generation and control of expression driven EEG signals, and verifies the feasibility of the classification and recognition method whose highest off recognition rate is 85%.

expression control; electroencephalogram; mechanism modeling; wavelet transform; artificial neural network

2016-10-15。 作者簡介:張小棟(1967—),男,教授,博士生導(dǎo)師。 基金項目:國家自然科學(xué)基金資助項目(51275388)。

時間:2016-04-15

10.7652/xjtuxb201606001

TP911.7

A

0253-987X(2016)06-0001-08

網(wǎng)絡(luò)出版地址:http:∥www.cnki.net/kcms/detail/61.1069.T.20160415.1612.010.html