fir_window.m

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

%程序 7.7
%凯泽窗生成程序
fpts=input('输入频带截止频率=');
mag=input('输入所需的幅度值=');
dev=input('输入每个频带的波纹=');
[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('增益,dB');
wvtool(w);