0 引言
如今，腦電信號廣泛應用于情感識別研究領域。 英國倫敦大學的 Eimer 和 Holmes 使用分別帶有驚恐和中性表情的人臉圖片作為視覺刺激并利用事件相關電位來找出驚恐情感和中性情感間的差異[1]. 美國休斯頓大學的 Zouridakis 等人通過過濾將原始腦電信號映射到常見的五個頻帶上，再分別求出五個頻帶所對應的頻帶能量，以此作為腦電特征進行情感識別[2]. Frantzidis 等人以圖片作為視覺刺激對測試者進行試驗，他們將所得到的數據分為心情愉悅和心情低落，高興奮度和低興奮度，并利用馬氏距離分類法進行情感分類，具有很高的分類精度[3]. 雖然對情感分析的方法種類有許多[4-8],但通過事件相關去同步化/同步化方法研究情感識別的研究還有待深入，并且本研究在實現了當左右腦對應區域均為活動狀態時，通過左右腦系數差值百分比判斷哪一個區域活動較為劇烈的功能。
1 預備知識
當大腦皮質某區域受到感官、動作指令或想象運動等刺激而開始激活時，該區域的代謝和血流增加，同時進行信息加工可以導致相應頻段的腦電信號振蕩的幅度減低或者阻滯，這種電生理現象稱作事件相關去同步化（ ERD）[9-10]. 大腦在靜息或惰性狀態下表現出明顯波幅增高的電活動，將其稱為事件相關同步化（ ERS）[11-12]. ERD/ ERS 具有頻段特異性： 大腦振蕩的活動頻率范圍波動較大，這些頻率常被分為以下頻段： delta（ < 4Hz） ,theta（ 4 ~7 Hz） ,alpha（ 8 ~12 Hz） ,beta（ 13 ~ 30 Hz） 及 gamma 頻段（ > 30 Hz） .
Alpha 和 beta 頻段的 ERD 可作為大腦激活的指標，而 ERS 的出現則表明大腦處于失活狀態，提示皮質區域恢復到靜息或惰性狀態[13]. 本研究采用事件相關去同步化和同步化方法研究過濾后的腦電信號。 現有常用的 ERD/ERS 分析方法主要有時頻圖譜法（TFS），功率譜密度（PSD）曲線法，ERD / ERS 系數法，腦地形圖法，功率譜嫡分析法，小波嫡分析法等[14-15]. 本文將主要介紹 ERD / ERS系數法。
2 方法
在介紹 ERD/ERS 系數法之前必須先引入另外兩個概念：短時傅里葉變換（Short Time FourierTransform,STFT） 和功率譜密度 （ Power SpectralDensity,PSD）。
2. 1 短時傅里葉變換
一般情況下，ERD/ERS 的發生是與特定時間段的特定頻段相關的，使用可以同時提取時域與頻域信息的時頻分析方法是更符合實際情況的選擇。 因此基于短時傅里葉變換的時頻圖譜和功率譜密度法常用在誘發腦電信號的分析中。 短時傅里葉變換的公式如下：
公式涵義：在時域用窗函數去截信號，對截下來的局部信號作傅里葉變換，即在 t 時刻該段信號的傅里葉變換，不斷地移動 t,即不斷地移動窗函數的中心位置，即可得到不同時刻的傅里葉變換，這些傅里葉變換的集合，既是 Sx（ω，t）。
2. 2 功率譜密度
定義：當波的頻譜密度乘以一個適當的系數后將得到每單位頻率波攜帶的功率，這被稱為信號的功率譜密度。 功率譜密度（PSD）同樣是基于短時傅里葉變換（STFT）其具體公式如下：
2. 3 ERD /ERS 系數
ERD / ERS 系數的計算公式如下：
其中，tR 為參考時間，tA 為測試時間，在本次研究實驗中用當前選定的時間段作為參考時間，而間隔 0. 1 s 的后一個時間段作為測試時間。
2. 4 研究方法
整個實驗 ERD/ERS 判斷方法（ERD/ERS 系數零一法）：將原始腦電數據過濾為腦電極為 Fp2（右前額）、Fp1（左前額）、C4（右腦中央區域）、C3（左腦中央區域）為四個電極，beta 頻帶（14,18,22,26 和 30 Hz） 的腦電波數據。 設定時間段的時間長度為 0. 4 s,時間間隔為 0. 1 s,依次提取十一個時間段的腦電數據，計算十個事件相關去同步/同步系數。 將計算得到的同步去同步系數大于 0的值（表示同步）均標記為 1,小于 0 的值（表示去同步）均標記為 0,將 10 個數相加，若所得結果大于 7 則表明同步現象明顯，小于 3 則表明去同步現象明顯。
當左右腦對應區域均為 ERD 時，可以通過比較左右腦參考時間和測試時間功率譜差值的絕對值大小比較活動強度。 絕對值較大的代表能量下降多即活動強度大。 比較 Fp2 和 Fp1,C4 和 C3 腦電極十個時間段參考時間和測試時間功率譜差值的絕對值大小，若 Fp2 值大于 Fp1,證明 Fp2 區域能量下降的多，即活動強烈，C4,C3 同理。
3 實驗
3. 1 測試者
所有測試者均為在校大學生，年齡在 20 ~25歲之間，25 名男性，10 名女性。
3. 2 實驗準備
安裝了如下電極：Fp1、Fpz、Fp2、C3、C4、F7、F8、Fz、C3、C4、T3、T4、Cz、P3、P4、T5、T6、Pz、O1、O2、Oz,奇數表示大腦左側，偶數表示大腦右側，Fpz、Fz、Cz、Pz、Oz 表示大腦中軸，兩個耳垂電極A1、A2 及 Z 作為參考電極，并要求有地線接地。腦電極的具體位置根據 Neurofax EEG - 9100K（NKC,Japan）的國際電極位置系統的標準。 采樣頻率設為 1 000 Hz.
3. 3 實驗步驟
第一步，讓實驗受試者保持放松閉眼，并記錄其腦電波，時間大約 120 s,此過程的最大目的是為了讓實驗參與者能處于相當比較平和的狀態；第二步，讓受試者睜眼，給予高興圖片視覺刺激，引起其情感變化，并記錄其腦電圖；第三步，讓受試者閉眼放松約 120 s,再給予不高興圖片視覺刺激，引起其情感變化，并記錄腦電圖。
4 結果
計算測試者分別在頻段為 14、18、22、26、30Hz 時，Fp2、Fp1、C4、C3 四個腦電極上，十個時間段的 ERD/ERS 系數。 其具體的數據和圖形情況，以某測試者為例，如表 1 和圖 1.
表 1 表示 beta 波 14 Hz,視覺刺激為高興時，某測試者在各電極各時間段的 ERD/ERS 系數。
每行表示 14 Hz 時，每個時間段的上各電極的ERD / ERS 系數，每一列表示特定電極上各時間段的 ERD/ERS 系數。
圖1 第一行表示視覺刺激為高興時 5 個 beta頻段上特定電極的 ERD/ERS 系數圖，每子圖的橫坐標表示時間段，縱坐標表示 ERD/ERS 系數最大上限和最小下限。 第二行表示視覺刺激為不高興時5 個 beta 頻段上特定電極的 ERD / ERS 系數圖。
根據 ERD/ERS 系數零一法，分別計算所有男性和女性在視覺刺激時，beta 頻段，各電極上的ERD / ERS 系數零一值和的平均值來判斷其大腦活動狀態，其具體內容如表 2 所示。
由表 2 可知，視覺刺激為高興時，男性電極Fp2,Fp1 和 C3 均變現為 ERD,C4 表現為 ERS;視覺刺激為不高興時，男性的 Fp2 和 Fp1,C3 和 C4均表現為 ERD,可以推斷出男性對高興和不高興的視覺刺激反應均比較明顯。 視覺刺激為高興時，女性電極 Fp2,Fp1 和 C3 均變現為 ERD,C4 表現為 ERS;視覺刺激為不高興時，女性的 Fp2,Fp1,C3 和 C4 均表現為 ERS,可以推斷出女性只對高興的視覺刺激明顯。
當左右腦對應電極均表現為 ERD 時，即男性：高興刺激時的 Fp2 和 Fp1 電極，不高興刺激時的 Fp2 和 Fp1 電極，C4 和 C3 電極；女性：高興刺激時的 Fp2 和 Fp1 電極。 比較功率譜差值絕對值，來判斷大腦活動的強度，得出結果：視覺刺激為高興時，男性 Fp1 上的譜密度差值較大，即左前額區域活動較為強烈；視覺刺激為不高興時，男性 Fp2和 C4 上譜密度差值較大，即右額區域和右腦中央區域活動較為強烈，而女性 Fp2 上的譜密度差值較大，即右前額區域活動較為強烈。
5 結論
為了研究由視覺刺激產生的情感變化時大腦的活動狀態，本文采用了 ERD/ERS 系數零一法通過特定時間段細化研究受刺激階段的腦電信號，保證了狀態檢測的準確信，充分發揮 ERD/ERS 系數的定量性的特點，并在此基礎上研究了左右腦均活動時的強烈程度。 本文結論對不同性別針對不同情感的情感識別有一定參考價值。
參考文獻
[1]EIMER MARTIN,HOLMES AMANDA. An ERP studyon the time course of emotional face processing. [J].NeuroReport,2002,13（4）：1-5.
[2]EIMER M,HOLMES A. Event-related brain potentialcorrelates of emotional face processing[J]. Neuropsy-chologia,2007,45（1）：15-31.
[3]ZOURIDAKIS G,PATIDAR U,PADHYE N S,et al.Spectral power of brain activity associated with emotion-apilot MEG study[C]/ /17th International Conference OnBiomagnetism Advances in Biomagnetism-BIOMAG2010,IFMBE Proceedings. Dubrovnik,Croatia:Springer,2010:354-357.
