exa5_20.m

来自「新版Matlab 7辅助信号处理技术与应用」· M 代码 · 共 17 行

%--------------------------------------------------------------------------
% exa060603_kaiserord.m,   for example 6.6.3;
% to test kaiserord and to design bandpass FIR filter;
%-------------------------------------------------------------------------
fsamp = 8000;
fcuts = [1000 1300 2210 2410];
mags = [0 1 0];
devs = [0.01 0.05 0.01];
[n,Wn,beta,ftype] = kaiserord(fcuts,mags,devs,fsamp);
n = n + rem(n,2);
hh = fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');
[H,f] = freqz(hh,1,1024,fsamp);
plot(f,abs(H),'k')
grid on
title('凯塞窗函数设计的带通滤波器幅频特性');
xlabel('频率（f）');
ylabel('幅度');