T P1
T P3
Séries de Fourier
RAPPEL: Nous rappelons qu'un signal périodique x (t ) de fréquence f0=1/T se décompose en
une infinité de sinusoides de fréquences 0, f 0 , 2f0, 3f0, …,nf0, … qu'on appelle série de Fourier
dont on donne les amplitudes par les formules suivantes:
{
t +T t +T
{
+∞ +∞ 1 1
2π 2π 1 2 2π
x(t )=a0+∑ [Link](n. t )+∑ [Link](n. t )¿ a0= ∫ x(t )dt ;an= ∫ x(t ).cos(n. t )dt ;¿ ¿¿¿¿
n=1 T n=1 T T t 1
Tt T
1
(1)
Cette décomposition s'appelle la décomposition en série de Fourier trigonométique; il existe
une autre décomposiion en série de Fourier exponentielle:
n=+ ∞ j . 2 π . n. t
x ( t )= ∑ X n.e T
n=−∞
t +T t
1
1 − j . 2 π . n.
Avec : X n= .
T
∫ x ( t ). e T
. dt
t1
(2)
[Link] du signal train d'impulsions rectangulaires
Nous voulons dans un premier exemple étudier la décomposition en série de Fourier d'un
train d'impulsions rectangulaires de fréquence f 0 et d'amplitude 2a0, nous tracerons d'abord
une courbe donnant une somme de sinusoïdes puis nous calculerons la décomposition en SF
du signal considéré, Dans le deuxième programme on prendra un nombre d'harmoniques de
plus en plus élevé et à chaque fois on tracera la courbe de la somme considérée et on
remarquera les signaux résultants:
On rappelle que pour le signal suivant
A
t
1T D/2 −D/2
e
La décomposition en série de Fourier exponentielle , et trigonometrique donne
n=∞ n=∞ n=∞
x (t )= ∑ A . D . f 0 . sinc ( nD f e ) e j 2 πnf t= ∑ 2. A . D . f 0 . sinc ( nD f 0 ) . co s ( 2 πn f 0 t )= ∑ a n . co s ( 2 πn f 0 t )
e
n=−∞ n =0 n=0
Travail à effectuer:
1. Taper le programme suivant: Prendre une certaine valeur de f 0, de A et de D
calculer les coefficients a n pour n=0;n=1;n=2 et puis taper le suivant
programme
%Traçé d'une somme de sinusoides
t=0:0.05:10;
f0=input('input frequence alementfond: ');
a0=input('input amplitude a0: ');
a1=input('input amplitude a1: ');
a2=input('input amplitude a2: ');
x=a0*cos(2*pi*f0*t)+a1*cos(2*pi*2*f0*t)+a2*cos(2*pi*3*f0*t);
plot(t,x);
2. Que fait ce programme
.
3. Taper des programmes similaires en prennant un nombre d'harmoniques de plus
en plus élevé(on calculera d'abord les coefficients a n pour n=3;n=4;n=5… ). Que
remarquer vous.
%Pour le train d'impulsions rectangulaires nous allons tracer la
somme de ses 2N composantes exponentielles
%x(t)=Somme de n=-N à n=N de(A*D*fe*sinc(n*D*fe)exp(j*pi*n*fe*t)
t=-4:0.005:4;
N=input('number of harmonics freq+ = ');
A=0.5;fe=0.5;D=1/fe/2
xN=A*D*fe*ones(1,length(t));
for n=[Link]N,
xN=xN+A*D*fe*sinc(n*D*fe)*2*cos(n*2*pi*fe*t);
end;
plot(t,xN)
Que fait ce programme;prendre plusieurs exemples avec des
valeurs d'harmoniques de plus en plus élevé.Que remarquez
vous
4. Tracé du spectre:Taper le programme suivant en prenant plusieurs exemples
avec différents rapports D/Te,que peut on dire du spectre:
clear
2
>> n=-[Link];
>>Te=2;
>> D=Te/2;
>> SYMS t
>> X=abs(int(exp(-i*2*pi*n*t/Te),-D/2,D/2)*1/Te);
>> stem (X )
[Link] du signal triangulaire
pour le signal suivant
2π
t
T
On la décomposition en série de Fourier
∞
−8 1 2π
x (t )= 2 ∑ 2
cos (n t)
π k=0 (2 k +1) T
Taper un programme qui permet de tracer la somme de quelques harmoniques du signal
triangulaire considéré ci-dessus, que remarquer vous en prenant un nombre d'harmoniques de
plus en plus élevé.
Tracer le spectre de fréquence(amplitude seulement).
3
4