%
% Fig7_10.m - M-file for generating Figure 7-10
%
L = 128;
n = [0:L-1];
x1 = sin(2*pi*30.5*n/128); % generate sinewave of 30.5 Hz
X1 = abs(fft(x1));         % magnitude spectrum
beta = 8.96;
wn = (kaiser(L,beta))';    % Kaiser window
x2 = x1.*wn;               % windowing
X2 = abs(fft(x2));         % windowed magnitude spectrum
subplot(2,1,1),plot(n,X1),axis([0 64 0 70]),
title('(a) Sinewave at 30.5 Hz'), subplot(2,1,2),
plot(n,X2),axis([0 64 0 70]),title('(b) Using Kaiser window'),
xlabel('Frequency, Hz');