ex0818.m

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

% Chapter 08: Example 8.18 Design a filter
f = [0 0.3 0.4 0.6 0.7 1]; a = [0 0 1 1 0 0];
b = remez(17,f,a);
[h,w] = freqz(b,1,512)
plot(f,a,w/pi,abs(h))
legend('Ideal','remez Design')