ex0823.m

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

% Chapter 08: Example 8.23 Design a filter
n = 50;
f = [0 0.1 0.15 0.25 0.3 0.4 0.45 0.55 0.6 0.7 0.75 0.85 0.9 1]; 
m = [1 1 0 0 1 1 0 0 1 1 0 0 1 1];
b = remez(n,f,m);
[h,w] = freqz(b,1,512,2);
plot(w,abs(h)); grid on
xlabel('Normalized Frequency'); ylabel('magnitude(dB)');