fig624.m

高西全2009年最新版<数字信号处理>源码程序

%《数字信号处理--原理、实现及应用》第六章 图6.1.3绘图程序:fig613.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月
% 不同阶数巴特沃斯低通模拟滤波器幅频响应

wp=1;Rp=1;

N1=2;
[B1,A1]=butter(N1,wp,'s');			%计算巴特沃斯低通模拟滤波器系统函数系数
wk=0:0.01:3;
Hk=freqs(B1,A1,wk);
subplot(2,2,1);plot(wk,abs(Hk));grid on
xlabel('\Omega');ylabel('幅度')
axis([0,3,0,1.1])
hold on

N1=4;
[B1,A1]=butter(N1,wp,'s');			%计算巴特沃斯低通模拟滤波器系统函数系数
Hk=freqs(B1,A1,wk);
plot(wk,abs(Hk));

N1=8;
[B1,A1]=butter(N1,wp,'s');			%计算巴特沃斯低通模拟滤波器系统函数系数
Hk=freqs(B1,A1,wk);
plot(wk,abs(Hk));