ex0817.asv

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

% Chapter 08: Example 8.17 Design a filter
n = 23
F = [0 0.3 0.4 0.6 0.7 0.9];
A = [0 1 0 0 0.5 0.5]
b = firls(24,F,A,'hilbert');
freqz(b)
    for i=1:2:6,
        plot([F(i) F(i+1)],[A(i) A(i+1)],'--'), hold on
    end
[H,f] = freqz(b,1,512,2);
plot(f,abs(H)), grid on, hold off