program_10_2.m

《数字信号处理：基于计算机的方法（第三版）》资源下载第三部分。 著者： （美）米特拉著 作译者： 孙洪等译 ISBN号： 7-121-02652-X 出版日期： 2006-06 出版

% Program 10_2
% Design of Equiripple Linear-Phase FIR Filters
%
format long
fedge = input('Band edges in Hz = ');
mval = input('Desired magnitude values in each band = ');
dev = input('Desired ripple in each band =');
FT = input('Sampling frequency in Hz = ');
[N,fpts,mag,wt] = remezord(fedge,mval,dev,FT);
b = remez(N,fpts,mag,wt);
disp('FIR Filter Coefficients'); disp(b)
[h,w] = freqz(b,1,256);
plot(w/pi,20*log10(abs(h)));grid;
xlabel('\omega/\pi'); ylabel('Gain, dB');