iir10.m

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

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);

%function [db,mag,pha,Omega]=freqs_m(b,a,Omega_max);
Omega=[0:1:500]*Omega_max/500;
H=freqs(b,a,Omega);
mag=abs(H);
dB=20*log10((mag+eps)/max(mag));
pha=angle(H);