巴特沃斯滤波器matlab代码,MATLAB实现巴特沃斯数字滤波器

MATLAB实现巴特沃斯数字滤波器

前因：因为要准备保研面试，今年暑假就重新把烂尾的项目捡起来了。

为了提取采集到的脑电信号中有用的部分，想用数字带通滤波器实现，浏览了很多帖子。要不是只有代码，没有注释；要不就是只有理论，没有代码。索性自己写一篇，方便回顾。

1. 用↓观察频谱

f=fftshift(fft(b)); %b表示信号值data

w=linspace(-512/2,512/2,length(b)); %根据奈奎斯特采样定理，512/2为最大频率

plot(w,abs(f)); %Hz为单位

k、Hz等纵坐标如何判断(5.同理)

2. 频率转化

Fs=512

fp=0.5Hz - 50Hz

fs=0.25Hz - 55Hz

///

*Q：↑为何如此取值

A：为防止频谱泄露，滤波器并非完全垂直截止，需过渡衰减

回归正题，然后将单位为‘Hz’的模拟频率转化为单位为‘rad’的数字角频率

基础知识见本连接

(wp，ws分别为数