% Program 8_10
% Spectral Analysis Using the DFT
%
N = input('Type in the length of DFT = ');
T = input('Type in the sampling period = ');
freq = input('Type in the sinusoid's frequency = ');
k = 0:N-1;
f = sin(2*pi*freq*k*(1/T));
F = fft(f);
stem(k,abs(F));grid
xlabel('k');
ylabel('|X[k]|');