ex0825.asv

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

% Chapter 08: Example 8.25 Design a filter
b=remez(21,[0.05 1],[1 1],'hilbert');
bb=remez(20,[0.05 0.95],[1 1],'hilbert');
[h,w]=freqz(b,1,512,2);
[hh,ww]=freqz(bb,1,512,2);
subplot(211); plot(w,abs(h)); grid on
title('Highpass Hilbert');
xlabel('Normalized Frequency'); ylabel('magnitude(dB)');
subplot(212); plot(ww,abs(hh)); grid on
title('Bandpass Hilbert');
xlabel('Normalized Frequency'); ylabel('magnitude(dB)');