Université Sidi Mohammed Ben Abdellah

Faculté des Sciences et Techniques –Fès

Système mécatronique et ingénierie automobile

Rapport sur les travaux

Pratiques de traitement
analogique du Signal

Réalisé par : Encadré par :

 AMAL ETTALHAOUI  Hassane kabbaj
 YOUNNES ELOUARRATI
FEUILLE DE
TP1
Exercice n°1 :
1) Soit le signal x1(t) suivant :

x1(t) = e ∝t e jwt ; ᾳ= -3 ; w= 20 rd/s

programme :

Figure(1) : réel de x1(t) en fonction du temps

  Ajoutons sur le même programme la deuxième figure:

 Et on obtient le résultat suivant :

x1(t) en fonction du temps

2) On va essayer l’instruction stem au lieu de plot :
On va réecrire le programme précédent ,
en remplaçant l’instruction plot par stem :
 Les résultats sont représentés dans les deux figures suivantes :

Figure(1) : Partie réelle de x1(t) en fonction du temps

 la fonction stem donne des courbe discrets.
3) On reprend la question 1 en rajoutant un déphasage de pi/4 à l’origine en supposant le
signal x2(t) tel que :
x2(t) = e jπ / 4 e(α + jw)t

Avec : alpha= -3 et w= 20 rd/s

Programme :

 On obtient les résulats suivants :

Figure(3) : La partie réelle de x2(t) en fonction du temps

 Figure(4) : La partie imaginaire de x2(t) en fonction du temps

4) On va représenter les parties réelles et imaginaires des signaux x1(t)et x2(t) en une seule
figure , pour se faire on va utiliser l’instruction subplot :

Programme :

 On obtient la figure suivante :

Exercice n°2 :
1) On trace la fonction de transfert de second ordre suivante :
225
H(p) = p ²+3 p+ 225

Avec p=j*w

On fait varier la pulsation w d’une façon logarithmique entre 1 et 1000 rd/s :

Programme :
  On obtient la figure suivante :

2) On doir mesurer expérimentalement :

-Le gain statique du système
-La fréquence de coupure wc et le coefficient d’amortissement m :

D’après la figure , on remarque que le signal commence de 0dB ,Le gain statique H0 =1
La fréquence de coupure wo est définie par:

H(wo) = H max – 3dB

3) On compare les resultants trouvés aux valeurs théoriques :

Calcul du gain statique théoriquement :

On sait que le gain statique du système s’écrit sous la forme :

H 0=lim H ( p)
p→0

225
= lim p ²+3 p+225
p→0

225
= 225 = 1

 H0 =1

Calcul de fréquence de coupure :

On a une fonction de transfert du second ordre , sa forme est exprimée comme suit:
1
H(p) = 1+ 2 ε ∗p+ 1 ∗p ²
ω0 ω 0²

1
Or on a : H(p)= 1 ∗p ² + 3 ∗p +1
225 225

1 1
Par identification ,on trouve : ω 0²
= 225

D’ou : ω 0=√ 225 ¿15 rad/s.

FEUILLE DE TP2 :
« Représentation de
signaux discrets :
signaux de base »
Pour un index temporel n variant linéairement de -20 à 20 avec un pas de 1, on trace les
signaux suivants :

1) Signal delta ou Dirac numérique :

Programme :

On obtient le tracé suivant :

Figure(1) : Signal delta ou Dirac numérique

2) Signal échelon unite:
Programme:

Figure(1) :Signal échelon unité:

3) Signal exposant n ( A positive, , 0<alpha<1) :
Programme :

Figure(1) : Signal exposant n(A positive,0<alpha<1) :

4) Signal exposant n (A positif, -1<alpha<0) :

Programme :
Figure(1) : Signal
exposant n (A positif, -1<alpha<0):

5) Signal exposant n (A positif ,|alpha|>1 :

Programme :
 Figure(1) : Signal exposant n (A positif, |alpha|>1):

6) Signal sinusoïdal :
Programme :

 On obtient le tracé suivant :

Figure(1) : Signal sinusoïdal :

7) Signal sinusoïdal avec déphasage :
Programme :

 On obtient le Tracé suivant :

Figure(1) : Signal sinusoïdal avec déphasage :

8) Signal triangulaire centré (pair) de pentes 1 et -1 :
Programme :

 On obtient le tracé suivant :

Figure(1) : Signal triangulaire centré (pair)de pentes 1 et -1:

FEUILLE DE TP3 :
Première partie :
Transformée de Laplace
Exercice n°1 :
Soit le signal x(t) suivant :

x(t) = e−a t

avec : a est une constante

Calculons à la main la transformée de Laplace de x(t) :

+∞

X(f) = TF {x1(t)} ¿ ∫ e−(a + p) t dt

−∞

−(a + p )t
−e
¿[ ]
a+ p

1
= p+ a

On calcule la transformée de Laplace de x(t) à l’aide de Matlab :

Programme :

 Résultat :

On remarque qu’on a obtenu les mêmes résultats.

 Remarque:
- L’ instruction syms permet de déclarer les variables comme étant symboliques.
- L’instruction laplace : est une fonction qui permet d’obtenir la transformée de
Laplace d’un signal.
Exercice n°2 :
Soit le signal x(t) suivant :
−t /T
x (t )=1−e

Calculons à la main la transformée de laplace de x(t) :

+∞

X(f) = TF {x1(t)}=∫ (e ¿ ¿(− j 2 πft )¿ −e(− j 2 πft )∗e−t / T )dt ¿ ¿

0

=
[ e (− j2 πft ) ]
−¿ ¿
j2 πf

1
1
=p - p+
1
T

Avec w=2 πf et p =j*w

On calcule la transformée de Laplace de x(t) à l’aide de Matlab :

Programme :

Résultat :

On o a obtenu le même résultat.

Exercice n°3 :
Soit F1(s) la fonction dans le domaine de Laplace donnée par :
3 s+2
F1(s) = s ²+3 s +2

On doit retrouver f1(t) la fonction image de F1(s), pour cela on peut utiliser deux méthodes :

 Premiére méthode en utilisant l’instruction ilaplace : (Méthode directe) :

Programme :

Résultat :

 Deuxiéme méthode en utilisant l’instruction residue : (Méthode indirecte )

Programme :

Résultat :

r1 r2 1
Or on a : F(s)= + On sait que : −1
TF ([ ¿)=e−at
( s−p 1) ( s−p 2) s +a

On remplace r1 par 4 , r2= -1,p1=-2 et p2=- En effectuant la transformée de Laplace inverse , on trouve :
f(t)=4e
−2
- e−1
 4 1
1,on trouve : F(s)= -
( s+2) ( s+1)

Deuxième partie:
Transformée de
Fourier
Exercice n°4 :
Soit le signal :
−t ²
x (t )=e

Programme

Résultat :

Exercice n°5 :
Soit le signal : x (t )=t e−t

En utilisant Matlab, on donne l’expression de la transformée de fourier de F (ω):

Programme :

Résultat :
Exercice n° 6:
Soit le signal discret x(n) défini par : x(n) = exp(-0.2*n²),pour n variant de 0 à 20

En utilisant l’instruction fft,on doit tracer le signal x(n) ,le module et la phase de sa
transformée de fourier dans un même espace en utilisant l’instruction subplot.

Remarque : L’instruction fft est la transformée de fourier rapide pour un signal discret, par
contre l’instruction ft est la transformée de fourier pour un signal continu.
Programme :

Représentation des figures :

 FEUILLE DE TP4 :
Première partie : Notions
élémentaires de Simulink
  Avant filtrage :
Dans cette première partie du Tp ,on va apprendre à manipuler le toolBox Simulink de
Matlab ,pour cela on ouvre Matlab ,on clique sur le symbole coloré « Simulink »,on clique
sur File  New Model (Dans ce modèle on va créer notre projet Simulink) , et on enregistre
notre [Link]

 On visualise deux signaux sinusoïdes :

On active Source  Sine Wave
Composants d’affichage Sinks Scope

 Remarques :

Au départ, le bloc Scope ne comportait qu’un seul axe d’affichage. Pour en ajouter un
second, il faut double-cliquer sur le composant, cliquer sur l’icône des paramètres, puis
modifier le champ "Nombre d’axes" à 2, et valider avec OK.

Concernant la configuration du premier signal sinusoïdal, un double-clic sur le bloc

Sine Wave permet d'accéder aux paramètres suivants :
 Amplitude : 1
 Polarisation (bias) : 0
 Fréquence : 1 rad/s
 Phase : 0 rad
 Temps d’échantillonnage (sample time) : 0
De la même manière, pour le deuxième signal, un double-clic sur le bloc Sine Wave1
permet de saisir les valeurs suivantes :
 Amplitude : 2
 Polarisation (bias) : 0
 Fréquence : 2 rad/s
 Phase : 0 rad
 Temps d’échantillonnage (sample time) : 0
 Après la simulation, une figure apparait :

Figure : Signaux avant filtrage :

N.B : Pour changer l’amplitude du signal, on clique que sur le bouton droit

 Autoscale
  Après filtrage :
On va appliquer un filtre à un signal sinusoïdal :

On active Source  Sine Wave

Composants d’affichage Sinks Scope
Signal Processing blockSet  Filtering  Filter Design
Simulink  Signal Routing  Mux
Remarque : On a utilisé un multiplexeur pour visualiser les deux signaux à la fois

Pour wc=5 rd/s et w=10 rd/s :

On va cliquer deux fois sur le filtre : fréquence de coupure, wc= 5 rd/s

De la même façon, on clique deux fois sur la sinusoïde : pulsation w =10 rd/s

Après la simulation on obtient la figure suivante :

Figure : Signaux après le filtrage

Remarque : On remarque qu’il ya un décalage, qui est du au passage du signal dans le filtre
(filtrage).

Commentaire : On en déduit que le filtre est de type passe-bas, car sa fréquence de coupure
ωc\omega_cωc est inférieure à la fréquence du signal ω\omegaω (ωc<ω\omega_c < \
omegaωc<ω). Dans ce cas, le filtre ne laisse passer que les composantes fréquentielles
inférieures à ωc\omega_cωc, ce qui conduit à une atténuation du signal sinusoïdal.

Pour éviter cette atténuation, il est nécessaire de choisir une fréquence de coupure supérieure
à celle du signal, soit ωc>ω\omega_c > \omegaωc>ω. Ce principe est appliqué dans le cas
suivant :

Pour ωc=10 rad/s\omega_c = 10 \, \text{rad/s}ωc=10rad/s et ω=5 rad/s\omega = 5 \, \

text{rad/s}ω=5rad/s :

 On double-clique sur le bloc du filtre pour modifier la fréquence de coupure à 10

rad/s.

 Ensuite, on ajuste la fréquence de la sinusoïde à 5 rad/s.

Le filtre utilisé dans ce cas est un filtre Butterworth, connu pour sa réponse
fréquentielle sans ondulations dans la bande passante.

Après la simulation, on observe les deux signaux après le filtrage :

 Remarque : Pour changer le pas, on clique sur Simulink Configuration Parameters

 Stop time = 5s (par exemple)  Apply Ok
 Commentaire : On contate que le signal n’a pas changé,et qu’ il n’y a pas
d’atténuation du signal cela du au fait de choisir wc>w
 On change le filtrage :

On va doucle-cliquer sur le filtre Design Method Choisir Cheby Shev1  Ok

Avec : la fréquence de coupure = 10 rd/s, et la sinusoïde a pour pulsation

w =5 rd/s.

On obtient la figure suivante :

 Commentaire : On remarque qu’on a obtenu la même figure qu’auparavant , d’où le

changement de filtre n’a aucun impact sur l’allure du signal.
Deuxième partie :
Réponse indicielle
d’un système du
ème
4 ordre
Dans le Simulink library browser, on va choisir « Continuous », puis on va choisir la boite
« transfer function »et on la dépose sur notre nouveau document Simulink ,on va choisir
aussi une source d’entrée Step :

On va utiliser la fonction de tranfert suivante :

1
H(s)= 4 3
s +1.197 s +1.717 s ²+1.025 s+0.379

Pour entrer ses coefficients, on double cliquer sur la fonction de transfert, une fenêtre des
paramètres apparait :

On entre :

Numerator coefficicents :

N = [1]

Denominator coefficients :

D = [1 1.197 1.717 1.025 0 .379]

Ce qui nous permet d’obtenir la courbe suivante :
 Dernière
partie : TP N°5
But de ce TP : visualisation de différents signaux dans le domaine spectral par la méthode
FFT de Matlab ; En effet elle représente la transformée de Fourier rapide d’un signal discret
avec une période d’échantillonnage Ts.

 En ajoutant un filtre :