% DFT Semplificata (Seno) 


clear;
close all;

Fs=8000;
Ts=1/Fs;

N=80;  % dimensione in numero di campioni della finestra di segnale da analizzare
f0=1/(N*Ts); % frequenza fondamentale di analisi

% Test pattern (10, 20, 40 Hz)

for n=0:N-1
    x(n+1)=10+sin(2*pi*100*n*Ts)+0.5*sin(2*pi*400*n*Ts)+0.25*sin(2*pi*2000*n*Ts)+0.125*sin(2*pi*3500*n*Ts);
end

plot(x);
% DFT Seno

for k=1:(N/2)
X(k)=0;
for n=0:N-1
    X(k)=X(k)+x(n+1)*sin(2*pi*((k-1)*f0)*n*Ts);   % misura della componente a 100 Hz
end
end

%calcolo dell'ampiezza effettiva;
    
X=(X/N)*2;

X(1)=X(1)/2;    % Valor medio (Frequenza 0)

% Calibrazione asse frequenziale (multipli interi della fondamentale di analisi (f0=1/NTs) 

F=f0:f0:f0*N/2;
figure;
plot(F,X);

    






