TP2 Analyse de Fourier
TP2 Analyse de Fourier
TP2 Analyse de Fourier
Analyse de Fourier
I- Séries de Fourier :
Introduire la fonction suivante sous Matlab. Celle ci permet de calculer une approximation
d’un signal carré à partir d’un nombre limité de coefficients de Fourier.
function sf(n)
%
% Calcul de la série de Fourier correspondant au signal carré de période T,
% et échantillonné avec une période Ts,
% La somme de la série jusqu'à l'harmonique de rang n est calculée et tracée
% N : nombre de points du signal
%
m=floor(n/2);
Ts=0.22e-3;
T=24.e-3;
N=500;
p=(-m-1):m ;
cn=((-1).^p)./(2*p+1)/pi;
cc=(2*pi*i*Ts)*(2*p+1)/T;
for k=0:N-1
vn=exp(cc*k);
s(k+1)=cn*vn';
end
t=(0:N-1)*Ts;
f=(0:N-1)/N/Ts;
subplot(211); plot(t,real(s)); grid; title('spectre')
subplot(212); plot(f,abs(fft(s)/N)); grid; title('signal’)
1- Exécuter cette fonction sous Matlab pour différentes valeurs du nombre n d’harmoniques.
Retrouver le signal considéré et calculer les coefficients de sa série de Fourier. Comparer
alors avec les calculs mis en œuvre dans cette fonction. Quelles sont vos conclusions.
2- Modifier cette fonction pour considérer un autre signal périodique que vous avez déjà
étudié théoriquement , de période et de fréquence d’échantillonnage fixée.
II-Transformée de Fourier
La transformée de Fourier est calculée par Matlab par un algorithme rapide, par la fonction
fft( ), celle ci nécessite que le signal échantillonné soit de longueur égale à une puissance de 2.
Sinon le signal considéré est prolongé par des zéros pour respecter cette condition.
2- Réaliser le signal y(t) produit de convolution du signal x(t) avec lui même
L’instruction de ce produit de convolution est conv2(x,x)
>> y=conv2(x,x) ;
Représenter le signal y(t) et comparer avec le résultat obtenu par calcul.
3- Générez un signal sinusoïdal x(t) de fréquence fx=100Hz échantillonné à fe=8kHz sur 256
points. Modulez ce signal par un signal cosinus m(t) de fréquence fm=2500Hz appelée
fréquence porteuse (modulation : multiplication terme à terme des deux signaux).