fir_kaiserfilter.m

来自「FIR滤波器设计中的凯撒窗函数设计」· M 代码 · 共 15 行

%程序 7.34
%凯泽窗生成程序
fpts=input('输入频带截止频率=');
fs=input('输入采样频率=');
mag=input('输入所需的幅度值=');
dev=input('输入每个频带的波纹=');
[N,Wn,beta,ftype]=kaiserord(fpts,mag,dev,fs);
kw=kaiser(N+1,beta);
b=fir1(N,Wn,ftype,kw,'noscale');
%b=round(b*2^15);
[h,f]=freqz(b,1,512,fs);
plot(f,20*log10(abs(h)));grid;
xlabel('频率,Hz');ylabel('增益,dB');
wvtool(kw);
wvtool(b);