
close all; clear all; clc
%% Resolução do exemplo 9.2 do livro do Lathi

No = 32;
n = 0 : No-1;

xn = [ones(1,5) zeros(1,23) ones(1,4)];

for r = 0 : No-1
    
    Xr(r+1) = sum(xn .* exp(-j*r*2*pi/No*n))/32;
    
end

r = n;

subplot(2,1,1)
stem(r,real(Xr),'k')
xlabel('r'); ylabel('X_r')
axis([0 31 -0.1 0.3])
legend('SFTD por cálculo direto',0)

Xr = fft(xn)/No;

subplot(2,1,2)
stem(r,real(Xr),'k')
xlabel('r'); ylabel('X_r')
axis([0 31 -0.1 0.3])
legend('SFTD pela FFT',0)
