w_kaiser.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 19 行

%凯泽窗，w_Kaiser
M=50;n=[0:49];
beta=5.6533;wk=Kaiser(M,beta);
subplot(2,2,1);stem(n,wk);
axis([0 50 0 1]);
title('凯泽窗 beta=5.6533') 
[db,mag,pha,grd,w]=freqz_m(wk,[1]);
subplot(2,2,3);plot(w/pi,db);
axis([0 1 -60 1]);
title('振幅响应') 
%
beta=6.6533;wk=Kaiser(M,beta);
subplot(2,2,2);stem(n,wk);
axis([0 50 0 1]);
title('凯泽窗 beta=6.6533') 
[db,mag,pha,grd,w]=freqz_m(wk,[1]);
subplot(2,2,4);plot(w/pi,db);
axis([0 1 -60 1]);
title('振幅响应')