ep632.m

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

%《数字信号处理--原理、实现及应用》第六章 例6.2.2 设计程序：ep622.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月
% 用脉冲响应不变法设计DF

% T=1s
T=1;Fs=1/T;
wp=0.2*pi/T;ws=0.35*pi/T;rp=1;rs=10;
[N1,wc1]=buttord(wp,ws,rp,rs,'s');
[B1,A1]=butter(N1,wc1,'s');
%绘制带通滤波器的损耗函数曲线
fk=0:Fs/2/512:Fs/2;wk=2*pi*fk;
Hk=freqs(B1,A1,wk);
subplot(2,2,1);
plot(fk,20*log10(abs(Hk)));grid on
xlabel('f(Hz)');ylabel('幅度(dB)')
axis([0,0.5,-50,5]);title('(a) 模拟滤波器（T=1s）    ');
[Bz1,Az1]=impinvar(B1,A1,Fs);
Hz=freqz(Bz1,Az1,wk);
subplot(2,2,3);
plot(wk/pi,20*log10(abs(Hz)));grid on;
xlabel('\omega/\pi');ylabel('幅度(dB)')
axis([0,1,-50,5]);title('(b) 数字滤波器（T=1s）    ');
% T=0.1s
T=0.1;Fs=1/T;
wp2=wp/T;ws2=ws/T;rp=1;rs=10;
[N2,wc2]=buttord(wp2,ws2,rp,rs,'s');
[B2,A2]=butter(N2,wc2,'s');
%绘制带通滤波器的损耗函数曲线
fk=0:Fs/2/512:Fs/2;wk=2*pi*fk;
Hk=freqs(B2,A2,wk);
subplot(2,2,2);
plot(fk,20*log10(abs(Hk)));grid on
xlabel('f(Hz)');ylabel('幅度(dB)')
axis([0,5,-50,5]);title('(c) 模拟滤波器（T=0.1s）    ');
[Bz2,Az2]=impinvar(B,A,Fs);
Hz=freqz(Bz2,Az2,wk);
subplot(2,2,4);
plot(wk/pi,20*log10(abs(Hz)));grid on;
xlabel('\omega/\pi');ylabel('幅度(dB)')
axis([0,1,-50,5]);title('(d) 数字滤波器（T=0.1s）   ');