iir.m

来自「FIR滤波器设计」· M 代码 · 共 29 行

clc
wp=0.2*pi;wr=0.3*pi;Ap=1;Ar=15;
T=1;Omegap=(2/T)*tan(wp/2);Omegar=(2/T)*tan(wr/2);
[b,a]=afd_butt(Omegap,Omegar,Ap,Ar)
[C,B,A]=sdir2cas(b,a)
[dB,mag,pha,Omega]=freqs_m(b,a,pi)
[hn,x,t]=impulse(b,a);
subplot(224);  
plot(t,hn)
  title('冲激响应');ylabel('hn(t)');
  xlabel('时间（秒)');

[b,a]=bilinear(b,a,T);
[C,B,A]=dir2cas(b,a);
[db,mag,pha,w]=freqz_m(b,a);
subplot(221);plot(w/pi,mag)
title('数字滤波器幅度响应|H|(ej/Omegar)'); ylabel('|H|');
subplot(222);plot(w/pi,db)
title('数字滤波器幅度响应(db)'); ylabel('分贝数');
subplot(223);plot(w/pi,pha/pi)
title('数字滤波器相位响应'); ylabel('弧度');
  xlabel('以Pi为单位频率');
delta_w=2*pi/1000;
Ap=-(min(db(1:1:wp/delta_w+1)))
Ar=-round(max(db(delta_w+1:1:501)))