eg10_12.m

matlab7.0从入门到精通(光盘内容)

wp=0.2*pi;
ws=0.3*pi;
rp=0.5;
rs=20;
[n,wn]=ellipord(wp,ws,rp,rs,'s')
[z,p,k] = ellipap(n,rp,rs);
w = logspace(-1,1,1000);
h = freqs(k*poly(z),poly(p),w);
semilogx(w,abs(h));
grid