滤波器的设计（巴特沃斯滤波器）

fp=3000;
Rp=2;
fst=6000;
As=30;
OmegaP=2*pi*fp;
OmegaST=2*pi*fst;
[N,OmegaC]=buttord(OmegaP,OmegaST,Rp,As,'s');%求解滤波器的阶数N和3dB截止频率OmegaC
[b,a]=butter(N,OmegaC,'s');%求解N阶滤波器的具体参数b和a
fk=0:10:20000;%取点，从0-20000，每隔10取一个点
wk=2*pi*fk;%注意模拟滤波器用的频率都是模拟角频率，要乘上2pi的
Hk=freqs(b,a,wk);%对于取的每个点，求该处的频率响应大小
plot(fk/1000,20*log10(abs(Hk)),'linewidth',2);%横坐标单位是kHz,纵坐标单位是dB
grid;
title('巴特沃斯滤波器幅频响应')
xlabel('f/kHz');
ylabel('dB');
axis([0,16,-100,0]);%设置横纵坐标轴范围