ep642.m

来自「高西全2009年最新版<数字信号处理>源码程序」· M 代码 · 共 32 行

%《数字信号处理》第六章 例6.4.2 设计程序：ep642.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月
% 用双线性变换法设计DF

T=1;Fs=1/T;
wpz=0.2;wsz=0.3;
wp=2*tan(wpz*pi/2);ws=2*tan(wsz*pi/2);rp=1;rs=1 %预畸变校正转换指标
[N,wc]=buttord(wp,ws,rp,rs,'s');    %设计过渡模拟滤波器
[B,A]=butter(N,wc,'s');
%绘制模拟滤波器的损耗函数曲线
fk=0:1/512:1;wk=2*pi*fk;
Hk=freqs(B,A,wk);
subplot(2,2,1);
plot(fk,20*log10(abs(Hk)));grid on
xlabel('f(Hz)');ylabel('幅度(dB)')
axis([0,0.9,-100,5]);title('(a)         ');
[Bz,Az]=bilinear(B,A,Fs);      %用双线性变换法转换成数字滤波器
wk=0:pi/512:pi;
Hz=freqz(Bz,Az,wk);
subplot(2,2,2);
plot(wk/pi,20*log10(abs(Hz)));grid on;
xlabel('\omega/\pi');ylabel('幅度(dB)')
axis([0,1,-100,5]);title('(b)           ');
[N,wc]=buttord(wpz,wsz,rp,rs);      %调用buttord和butter直接设计数字滤波器   
[Bz,Az]=butter(N,wc);
wk=0:pi/512:pi;
Hz=freqz(Bz,Az,wk);
subplot(2,2,4);
plot(wk/pi,20*log10(abs(Hz)));grid on;
xlabel('\omega/\pi');ylabel('幅度(dB)')
axis([0,1,-100,5]);title('(b)           ');