program_10_5.m

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

% Program 10_5
% Lowpass Filter Design Using the Kaiser Window
%
fpts = input('Type in the bandedges = ');
mag = input('Type in the desired magnitude values = ');
dev = input('Type in the ripples in each band = ');
[N,Wn,beta,ftype] = kaiserord(fpts,mag,dev)
kw = kaiser(N+1,beta);
b = fir1(N,Wn, kw);
[h,omega] = freqz(b,1,512);
plot(omega/pi,20*log10(abs(h)));grid;
xlabel('\omega/\pi'); ylabel('Gain, dB');