%  Calcul de la transformée de Fourier
%  discrète
%  Licence de Physique, Sciences physique
%     1995, 1996
% -------------------------------------

clear all;
clg

% Initialisation
% ---------------
 
N = ?;          % Nombre de points
freq = ?;      % Fréquence
phase = ?;



tau =?;     % pas de temps de l'échantillonnage

% calcul de la fonction à analyser
% --------------------------------


for j = 0:N-1; 
    t(j+1) = tau*j;
    y(j+1) = ?;
end

%
%
%   pour tracer la fonction avec plus de points
%
for j = 0:10*N-1; 
    tt(j+1) = j*tau/10;
    yy(j+1) = ?;
end




%%%%  Exercice 5 %%%% 
%load mauna.dat;
%N = length(mauna);
%for i = 1:N
%y(i) = mauna(i,1);
%t(i) = i;
%end

% detrending


%p = polyfit(t,y,1)
%for i = 1:N
%lin(i) = p(1)*t(i) + p(2);
%y(i) = y(i) - lin(i);
%end
%%%%%            %%%%%

% Calcul de la tranformée de Fourier
% -----------------------------------

deuxpiN = 2*pi*sqrt(-1)/N;

for k=0:N-1
  somme = 0;
  for j=0:N-1
    somme = somme + ????;
  end
  yt(k+1) = somme;
  s(k+1) = ????;
end



% Calcul des fréquences
% ---------------------


for k=1:N/2
    f(k) = ???;
end
for k=N/2+1:N
    f(k) = ????;
end


% Partie graphique
% ----------------

subplot(211)
plot(tt,yy,t,y,'r+');

% pour ex 5:
%plot(t,y);

%title('Données filtrées');
  ylabel('Amplitude');
  xlabel('Temps')

subplot(212)
%plot(f,real(yt),'-',f,imag(yt),'--');
%title('Réelle (--) et Imag (- -)');
plot(f,s,'-');
%title('spectre d energie (--)');
  ylabel('Transformée')
  xlabel('Fréquence')
