ex07017.m

数字信号处理（机械工业出版社）的源码

% Chapter 7: Example 7.17 filter design
%
n = 10; Rp = 0.5; Rs = 20;
Wn = [100 200]/500;
[b,a] = ellip(n,Rp,Rs,Wn); 
freqz(b,a,128,1000); title('n=5 Chebyshev Type II Bandpass Filter')
figure
[y,t] = impz(b,a,101); stem(t,y)
title('Impulse Response of n=10 Elliptic Filter')