clear;
close all
figure(5)
[N1,wc1,beta1,ftype1]=kaiserord([0.2,0.3],[0,1],[0.01,0.1]);
b1=fir1(N1,wc1,ftype1,kaiser(N1+1,beta1));
w=kaiser(N1+1,beta1);
n=0:N1;
plot(n,w)
title('凱塞窗函數(shù)')

figure(6)
freqz(b1,1,1024);