matlab滤波器设计

在进行脑电信号处理时用到带通滤波和高通滤波，记录以便后续方便查找

1、带通滤波

 进行f1-f2Hz的带通滤波，采样率为fs，滤波器阶数为2阶

Wc1=2*f1/fs;                                          %下截止频率 f1Hz
Wc2=2*f2/fs;                                          %上截止频率 f2Hz
[b,a]=butter(2,[Wc1,Wc2],'bandpass');  % 二阶的巴特沃斯带通滤波
data=filter(b,a,data);

2、高通滤波 

 进行f1Hz的高通滤波，采样率为fs，滤波器阶数为2阶

Wc1=2*f1/250;                                          %截止频率为f1Hz
[b,a]=butter(2,Wc1,'high');                            % 二阶的巴特沃斯高通滤波
data=filter(b,a,data);

3、低通滤波

进行f1Hz的低通滤波，采样率为fs，滤波器阶数为2阶

 Wc1=2*f1/250;                                          %截止频率 f1Hz
 [b,a]=butter(2,Wc1,'low');                          % 二阶的巴特沃斯带通滤波
 data=filter(b,a,data);