program_10_4.m

来自「《数字信号处理：基于计算机的方法（第三版）》资源下载第三部分。 著者： （美）」· M 代码 · 共 11 行

% Program 10_4
% Kaiser Window Generation
%
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)
w = kaiser(N+1,beta); w = w/sum(w);
[h,omega] = freqz(w,1,256);
plot(omega/pi,20*log10(abs(h)));grid;
xlabel('\omega/\pi'); ylabel('Gain, dB');