% Exercice sur matlab.
% F. Chambat 1998.

% Partie 1. Test.

x=[4;3;7;-9;1;0;0;0]
y=fft(x)

% Partie 2. TF de la somme de deux sinusoïdes.

t=0:0.001:0.6;				% Vecteur temps
x=2*cos(2*pi*50*t)+cos(2*pi*120*t);	% Signal : frequences 50 et 120 Hz
y=fft(x);
spy=y.*conj(y)/(600*600);		% Energie spectrale
plot(spy)
print -deps exam98-spy1			% Sortie PostScript
sauv=spy';
save exam98-spy1.dat sauv -ascii	% Sortie Ascii

% Représentation em fonction de la fréquence.
nyquist=1000/2 				% Fréquence de Nyquist
dnu=1/0.6				% Pas en fréquence
n=nyquist/dnu				% Nombre de points
f=0:dnu:(n-1)*dnu;			% Vecteur des fréquences
plot(f,spy(1:300))
print -deps exam98-spy2			% Sortie PostScript
sauvb(:,1)=f';				% Matrice a sauvegarder
sauvb(:,2)=spy(1:300)';
save exam98-spy2.dat sauvb -ascii	% Sortie Ascii