TP2 Analyse de Fourier

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 2

Licence GTR , TP Traitement du Signal :

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.

1- Construire le signal porte x(t) tel que:


>> a=zeros(1,8) ;
>> b=ones(1,16);
>> x=[a b a];
Fournir la densité spectrale d’énergie du signal x(t)

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.

Obtenir la représentation de la densité spectrale d’énergie du signal y(t) à partir de sa


transformée de Fourier, puis en n’utilisant que la transformée de Fourier du signal x(t).
Utiliser pour cela la propriété de la Transformée de Fourier d’un produit de convolution de
deux signaux.

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).

Comparez les spectres des signaux. Expliquer le résultat.

Vous aimerez peut-être aussi