事件相关电位的呈现是在特定频带下的时域波形,从这个角度出发的话,其处理就涉及到滤波;此外由于脑电信号易受噪声干扰,所以单独测出的每两次的数据肯定是不同的,甚至会差别很大,即使是在同一条件下,这就需要去除噪声(这里不做详细介绍)、去除基线和叠加平均。
简单的处理流程:
- 滤波
- 分段&加标签
- 去基线
- 叠加平均
滤波(filt raw): 滤波的作用有两个,第一,取出我们关注的ERP特征最明显的频段的数据,比如8~13Hz;第二,滤除低频漂移(线性趋势),有的时候观察原始波形,发现会有一个线性上升或者下降,经过低频滤波便可去除,比如过滤掉3Hz以下的数据。
一个典型的MATLAB滤波函数如下:
% 滤波iir
function r=filtraw(raw,fs,lowcut,highcut,order)
sizeraw=size(raw);
udfilter=designfilt('bandpassiir','FilterOrder',order, ...
'HalfPowerFrequency1',lowcut,'HalfPowerFrequency2',highcut, ...
'SampleRate',fs);
for i=1