%
% ksw.m - Kaiser window function and its magnitude response
%
L = 41; beta = 8; n = 0:40;
wn = kaiser(L,beta);
omega = -pi:2*pi/200:pi;      % frequency range
Hd = freqz(wn,1,omega);       % frequency response
mag = 20*log10(abs(Hd));      % in dB scale
subplot(2,1,1), plot(n,wn), title('Kaiser window, L=41, beta=8'),...
ylabel('Amplitude'); axis([0 (L-1) 0 1]);
subplot(2,1,2), plot(omega,mag), title('Magnitude Response'),...
    xlabel('Frequency'), ylabel('Magnitude(dB)'), grid;
axis([-3.15 3.15 -80 40]);