ex07018.m

来自「数字信号处理（机械工业出版社）的源码」· M 代码 · 共 8 行

% Chapter 7: Example 7.18 filter design
%
Wp = 40/500; Ws = 150/500;
Rp = 3; Rs = 60;
[n,Wn] = ellipord(Wp,Ws,Rp,Rs)
[b,a] = ellip(n,Rp,Rs,Wn);
freqz(b,a,512,1000); title('n=4 Elliptic Lowpass Filter')