fig623.m

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

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

wp=2;Rp=2;

N1=4;
[B1,A1]=butter(N1,wp,'s');			%计算巴特沃斯低通模拟滤波器系统函数系数
wk=0:0.01:5;
Hk=freqs(B1,A1,wk);
subplot(2,2,1);plot(wk,abs(Hk));
line([0,1.72],[0.9,0.9]);line([0,2],[0.707,0.707])
line([1.72,1.72],[0,0.9]);line([2,2],[0,0.707])
line([0,3.5],[0.1,0.1]);line([3.5,3.5],[0,0.1])
xlabel('Ω');ylabel('|Ha(jΩ)|')
axis([0,5,0,1])
subplot(2,2,2);plot(wk,20*log10(abs(Hk)));
line([0,1.72],[20*log10(0.9),20*log10(0.9)]);line([0,2],[-3,-3])
line([1.72,1.72],[-35,20*log10(0.9)]);line([2,2],[-3,-35])
line([0,3.55],[-20,-20]);line([3.55,3.55],[-40,-20])
xlabel('Ω');ylabel('-A(Ω)')
axis([0,5,-25,0])