EEGLAB使用流程

MATLAB版本 2016a; EEGLAB版本 v2020.0

數(shù)據(jù)的準(zhǔn)備

%matlab code clear all ;clc; filepath = strcat('200mM yuanshi\subject\subject3_',num2str(1),'.txt') sig=importdata(filepath); x=sig.data(2:1:end,2:2:end);%偶數(shù)項(xiàng)為腦電信號(hào) event = zeros(1,length(x(:,1))); kaishi = 2.5;%這是應(yīng)看到，5s一次刺激，后面改進(jìn)使用算法尋找事件標(biāo)記的開始 event(kaishi*256:kaishi*256 + 10) = 1; event((kaishi+5)*256:(kaishi+5)*256 + 10) = 1; event((kaishi+10)*256:(kaishi+10)*256 + 10) = 1; event((kaishi+15)*256:(kaishi+15)*256 + 10) = 1; event((kaishi+20)*256:(kaishi+20)*256 + 10) = 1; event((kaishi+25)*256:(kaishi+25)*256 + 10) = 1; x(:,end+1) = event;%事件標(biāo)記通道 x = x';

通道定位的文件

Standard-10-20-CCostmos.locs:

1 -18 0.51111 Fp1 2 18 0.51111 Fp2 3 -39 0.33333 F3 4 39 0.33333 F4 5 -90 0.25556 C3 6 90 0.25556 C4 7 -141 0.33333 P3 8 141 0.33333 P4 9 -162 0.51111 O1 10 162 0.51111 O2 11 -54 0.51111 F7 12 54 0.51111 F8 13 -90 0.51111 T3 14 90 0.51111 T4 15 -126 0.51111 T5 16 126 0.51111 T6 17 -43 0.5458 Sp1 18 43 0.5458 Sp2 19 0 0.25556 Fz 20 90 0 Cz 21 180 0.25556 Pz 22 180 0.51111 Oz 23 -90 0.63889 A1 24 90 0.63889 A2

數(shù)據(jù)的導(dǎo)入（MATLAB array）

本數(shù)據(jù)采集的時(shí)候沒(méi)有事件標(biāo)記通道，是將采集的被試者的腦電信號(hào)后，通過(guò)算法計(jì)算出來(lái)的事件標(biāo)記點(diǎn)，該數(shù)據(jù)一共有25個(gè)通道，24個(gè)通道是腦電通道，第25個(gè)通道是后期通過(guò)計(jì)算得出的事件標(biāo)記點(diǎn)。

File->Import data->Using EEGLAB functions and plugins->From ASCII/float or Matlab array.MATLAB工作空間變量的名字是 x,采樣率為256Hz(Data sampling rate),點(diǎn)擊OK即可,也可以其他數(shù)據(jù)格式的導(dǎo)入。

輸入數(shù)據(jù)的命名。點(diǎn)擊OK
可以通過(guò)Plot->Channel data(scroll)把數(shù)據(jù)EEG波形圖畫出來(lái)觀察一下。

現(xiàn)在點(diǎn)擊REJECT就會(huì)把選中的數(shù)據(jù)丟失掉，就可以觀察將不正常的數(shù)據(jù)剔除。只要是類似的界面，都可以實(shí)現(xiàn)數(shù)據(jù)的剔除

確定事件標(biāo)記通道

File->Import event info->From data chanel

因?yàn)榈?5個(gè)通道是事件標(biāo)記通道，所以輸入25，也可以采取其他事件標(biāo)記形式。
也可以在Edit -> Events values下面窗口進(jìn)行事件標(biāo)記的添加、命名和刪除。

通道定位

Edit->Channels locations，下拉列表選擇第二個(gè)參數(shù)，點(diǎn)擊OK即可
彈出下面窗口
點(diǎn)擊Read locations,選擇通道定位文件，我的文件是Standard-10-20-CCostmos.locs(不同實(shí)驗(yàn)定位文件不同，需要根據(jù)實(shí)驗(yàn)修改文件)，直接自動(dòng)監(jiān)測(cè)，點(diǎn)擊OK。

數(shù)據(jù)的濾波

只說(shuō)明最經(jīng)典的FIR濾波器,Tools->Filter the data->Basic FIR filter(new,dault)

第一個(gè)參數(shù)是高通濾波，第二個(gè)是低通濾波，當(dāng)兩個(gè)參數(shù)都有輸入的時(shí)候就是帶通濾波，但是EEGLAB官網(wǎng)不推薦帶通濾波，會(huì)出現(xiàn)錯(cuò)誤，所以要高通低通濾波兩次，這里參數(shù)選擇0.5Hz-45Hz。將濾波的數(shù)據(jù)命名為subject_FIR。

運(yùn)行ICA，去除偽跡

Tools->Decompose data by ICA，其他默認(rèn)即可，點(diǎn)擊OK,
ICA運(yùn)行的時(shí)候是可以打斷的，但是打斷之后再運(yùn)行是重新開始計(jì)算。
單純查看ICA分量：Plot->Component maps等
查看并剔除ICA分量：Tools->Inspect/label component by maps

點(diǎn)擊每個(gè)分量進(jìn)行查看判斷是否為偽跡

點(diǎn)擊ACCEPT將其轉(zhuǎn)化為REJECT就是將該分量舍棄。點(diǎn)擊OK。

繪制ERP

Tools->Extract epochs

第一個(gè)參數(shù)(bengin):事件標(biāo)記點(diǎn)
第二個(gè)參數(shù)(-0.5 1.5):數(shù)據(jù)選取事件標(biāo)記點(diǎn)前面0.5秒和事件標(biāo)記點(diǎn)1.5秒
第三個(gè)參數(shù)(subject_FIR_epochs):新數(shù)據(jù)組合的名字
點(diǎn)擊OK，

數(shù)據(jù)的基準(zhǔn)線是事件標(biāo)記點(diǎn)前500ms。
之后就可以在Plot菜單欄中繪制你想要的ERP波形了

Plot->Sum/Compare comp.ERPs可以繪制多個(gè)數(shù)據(jù)的ERP的均值方差對(duì)比

第一個(gè)參數(shù)是選用哪個(gè)數(shù)據(jù)，直接填寫數(shù)據(jù)的編號(hào)集

總結(jié)

以上是生活随笔為你收集整理的EEGLAB使用流程的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。