matlab程序p186.txt

实现butterworth滤波器,能够输出幅频特性曲线

close all;
%把数字滤波器的频率特征转换成模拟滤波器的频率特征
Wp=32000000*2*pi; %截止频率
Ws=66000000*2*pi; %阻带起始频率
Rp=2; %通带最大衰减2dB
Rs=31; %阻带最小衰减30db
[n,Wn]=buttord(Wp,Ws,Rp,Rs,'s'); %选择滤波器的最小阶数n
[z,p,k]=buttap(n);     %创建Butterworth低通滤波器原型
[b0,a0]=zp2tf(z,p,k);  %零极点增益模型转换为传递函数模型
[b,a]=lp2lp(b0,a0,Wn); %实现低通向低通的转变

[h,w]=freqs(b,a); %求频率响应
plot(w/(2*pi),20*log10(abs(h))); %绘出频率响应曲线
grid on;
XLabel('频率(Hz)');
yLabel('幅频响应(db)');
xLim([2000000,80000000]);
ylim([-45,0]);