exa5_20.m

matlab应用实例。如：FIR、IIR数字滤波器设计等。

%--------------------------------------------------------------------------
% 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('幅度');