我最近记录的脑电信号的采样率为256Hz.The信号将在4-64赫兹频段通过。我需要一个代码来滤除脑电data.Is。
发布于 2014-01-12 13:09:34
您可以应用50或60 Hz的陷波滤波器
眼球运动产生的伪影通常有2-5 Hz的频率范围,所以你可以在那里应用高通滤波器。
小波具有使用小波包分解过滤噪声(硬阈值和软阈值)的阈值机制。
发布于 2014-05-01 20:47:09
在Fs =1000 Hz的情况下,使用此代码能够过滤信号并提取特征频带(alpha,beta,...)
S = "your EEG-Data-Row";
waveletFunction = 'db8' OR 'sym8' ;
[C,L] = wavedec(S,8,waveletFunction);
%% Calculation The Coificients Vectors
cD1 = detcoef(C,L,1); %NOISY
cD2 = detcoef(C,L,2); %NOISY
cD3 = detcoef(C,L,3); %NOISY
cD4 = detcoef(C,L,4); %NOISY
cD5 = detcoef(C,L,5); %GAMA
cD6 = detcoef(C,L,6); %BETA
cD7 = detcoef(C,L,7); %ALPHA
cD8 = detcoef(C,L,8); %THETA
cA8 = appcoef(C,L,waveletFunction,8); %DELTA
%%%% Calculation the Details Vectors
D1 = wrcoef('d',C,L,waveletFunction,1); %NOISY
D2 = wrcoef('d',C,L,waveletFunction,2); %NOISY
D3 = wrcoef('d',C,L,waveletFunction,3); %NOISY
D4 = wrcoef('d',C,L,waveletFunction,4); %NOISY
D5 = wrcoef('d',C,L,waveletFunction,5); %GAMMA
D6 = wrcoef('d',C,L,waveletFunction,6); %BETA
D7 = wrcoef('d',C,L,waveletFunction,7); %ALPHA
D8 = wrcoef('d',C,L,waveletFunction,8); %THETA
A8 = wrcoef('a',C,L,waveletFunction,8); %DELTA
希望这能有所帮助。
https://stackoverflow.com/questions/13801270
复制相似问题