eg10_12.m

来自「matlab7.0光盘内容。都是一些具体的例子。」· M 代码 · 共 11 行

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