EEGLAB使用流程

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

数据的准备

%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);%偶数项为脑电信号
event = zeros(1,length(x(:,1)));
kaishi = 2.5;%这是应看到，5s一次刺激，后面改进使用算法寻找事件标记的开始
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;%事件标记通道
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

数据的导入（MATLAB array）

本数据采集的时候没有事件标记通道，是将采集的被试者的脑电信号后，通过算法计算出来的事件标记点，该数据一共有25个通道，24个通道是脑电通道，第25个通道是后期通过计算得出的事件标记点。

File->Import data->Using EEGLAB functions and plugins->From ASCII/float or Matlab array.MATLAB工作空间变量的名字是 x,采样率为256Hz(Data sampling rate),点击OK即可,也可以其他数据格式的导入。

输入数据的命名。点击OK
可以通过Plot->Channel data(scroll)把数据EEG波形图画出来观察一下。

现在点击REJECT就会把选中的数据丢失掉，就可以观察将不正常的数据剔除。只要是类似的界面，都可以实现数据的剔除

确定事件标记通道

File->Import event info->From data chanel

因为第25个通道是事件标记通道，所以输入25，也可以采取其他事件标记形式。
也可以在Edit -> Events values下面窗口进行事件标记的添加、命名和删除。

通道定位

Edit->Channels locations，下拉列表选择第二个参数，点击OK即可
弹出下面窗口
点击Read locations,选择通道定位文件，我的文件是Standard-10-20-CCostmos.locs(不同实验定位文件不同，需要根据实验修改文件)，直接自动监测，点击OK。

数据的滤波

只说明最经典的FIR滤波器,Tools->Filter the data->Basic FIR filter(new,dault)

第一个参数是高通滤波，第二个是低通滤波，当两个参数都有输入的时候就是带通滤波，但是EEGLAB官网不推荐带通滤波，会出现错误，所以要高通低通滤波两次，这里参数选择0.5Hz-45Hz。将滤波的数据命名为subject_FIR。

运行ICA，去除伪迹

Tools->Decompose data by ICA，其他默认即可，点击OK,
ICA运行的时候是可以打断的，但是打断之后再运行是重新开始计算。
单纯查看ICA分量：Plot->Component maps等
查看并剔除ICA分量：Tools->Inspect/label component by maps

点击每个分量进行查看判断是否为伪迹

点击ACCEPT将其转化为REJECT就是将该分量舍弃。点击OK。

绘制ERP

Tools->Extract epochs

第一个参数(bengin):事件标记点
第二个参数(-0.5 1.5):数据选取事件标记点前面0.5秒和事件标记点1.5秒
第三个参数(subject_FIR_epochs):新数据组合的名字
点击OK，

数据的基准线是事件标记点前500ms。
之后就可以在Plot菜单栏中绘制你想要的ERP波形了

Plot->Sum/Compare comp.ERPs可以绘制多个数据的ERP的均值方差对比

第一个参数是选用哪个数据，直接填写数据的编号集