Département du Génie Electrique

Module : Traitement numérique du signal et de

l’image

Cycle Ingénieur : SiCoM1

Semestre 2
Compte rendu : TP 2

Réalisé par : Encadré par :

Redouane Ayadi Pr. Farid Abdi

AU : 2019-2020
TP3 : Transformation de Fourier
 discrète , TZ et analyse
spectrale
1) Programmation de la TFD.

1-1) Définition théorique

La TFD d’un signal est définie par la relation suivante :
N−1 2π
−j kn
Xp [k ]= ∑ xp [n ]e N

n=0

2) Application
Nous avons commencé de générer une fonction TFD qui permet de calculer la TFD d’un
signal quelconque.

Après on veut appliquer cette fonction sur trois signaux différents et on visualise les résultats
trouvés.
On génère un signal composé de trois sinusoïdes pour calculer sa TFD, pour cela on
programme une fonction qui a pour entrée trois fréquences f1, f2, f3, et le nombre des
échantillons N.
>> generer_signal(10,15,15,20)

>> GenererSin(10,15,15,20)
Théoriquement :
N−1
3
TF ( f )= ∑ [δ ( f + K f 1 ) ¿−δ ( f −K f 1 ) ]+[δ ( f + K f 2 )−δ ( f −K f 2 ) ]+ ¿ ¿]
2 j k=0

Remarque : puisque on a calculé la TF en valeur absolue, on obtient que les valeurs positives.
Pour 2N échantillons : N=100, f1=10Hz, f2=15Hz, et f3=30 Hz, on obtient les résultats
suivants :

Impulsion de Dirac définie par d= [1, zéros(1, N-1)]

Avec N=8 on trouve l’allure suivante.

Pour analyser le résultat trouvé on fait appel à la théorie : δ ( n )=1 si n=0 et δ ( n )=0 si n ≠ 0
Pour calculer la TFD, on utilise la formule suivante :
k=N −1 n
−2 π j k
P ( n )= ∑ δ (k ) ⅇ N

k=0

On applique la TFD sur un signal auditif.

Cette partie est divisée sur trois sous parties :

1- On génère un signal auditif, on le génère sur un nombre N des échantillons, et on

visualise ses composants
On choisit audio= « SCHUBERT.WAV », et N=1024
Résultat

2- On calcule la fft et la TFD de ce signal et on compare les résultats :

Pour N=1024

On remarque qu’avec les deux commandes donne mêmes résultats

Pour N=2028

Conclusion
 Le résultat trouvé par les deux méthodes la TFD et la FFT sont équivalents.

3- On utilise la fonction FFT pour trouver la fréquence fondamentale de la note jouée par le violon,
et on déduit après la note équivalente.
 La transformée en Z
But du TP:

 Dans cette partie, nous essayons de manipuler la TZ et la TZI et les méthodes pour leur
détermination et leurs applications.
 Aussi on va voir la notion de filtrage, la présentation des filtres théoriques, la conception de
filtres et les filtres analogiques et numériques.

I. Détermination de la TZ et la TZI en utilisant le calcul

symbolique de Matlab :

a. Calcul théorique de la TZ

On considère les signaux suivants afin de déterminer leur transformée en Z théoriquement :

* h(n)=0.9 u ( n ) ;
*g(n)=u ( n ) −u ( n−10 ) ;

¿ f ( n )=cos (w 0 n)u(n);
Pour le signal H(n) :
D’après l’utilisation de la formule générale de la transformée en Z nous avons trouvé le
résultat suivant :
1 z
=
H(Z)= 1− 9 Z−1 z− 9
10 10
z 1
G(Z)= z−1 − 9
z (z−1)
2
z −Zcos(w 0)
F(Z)= 2
z −2 zcos ( w 0 )+1
 D’après l’utilisation de la formule théorique on passe maintenant à l’outil symbolique
Toolbox sous Matlab.

b. Calcul par Ztrans

 Résultat du calcul:

 On conclut que Les deux méthodes de calcul donnent le même résultat.

 Pour le signal H(z)=z/(z-9/10)

Pour le signal G(z)

Pour le signal cosinus

c. Calcul théorique de la TZI

Nous avons :

En utilisant la transformée en Z inverse pour retourner à la fonction de transfert en n pour cela

en utilisant la d
 Calcul de la TZI :
On a

A B
H ( Z )= −1
+ −1
1+0.5 Z 1−0.5 Z
−1
1+ Z −1
avec Z =−2 donc A=
−1
A=
1−0.5 Z
−1
2

−1
1+Z 3
avec Z =2 donc B=
−1
B=
1+0.5 Z
−1
2

−1
( )( ) ( )
n n
1 1 1 −1 −1 3 1
X ( z )= + , x ( n ) = +
2 1+ 0.5 Z −1
2 1−0.5 Z −1
2 2 2 2
 d. Calcul par iZtrans :

 Résultats de la TZI
 On remarque que avec les deux méthodes nous avons les mêmes résultats.
II. Détermination de la TZI en utilisant la décomposition en
fraction partielle

 Résultats obtenus:

♣ Pour cette partie on a vérifié le résultat de la fonction poly par la fonction roots, d’après

l’exécution on remarque que les pôles et les zéros sont justes .

La détermination des pôles et zéros par la fonction zplane.
Les pôles sont 0.5 et -0.5, et les zéros 0 et -1 d’après la présentation. On peut conclure que
notre système est stable car ses pôles sont à m’intérieur du cercle unité.

♣En utilisant la fonction residuez pour calculer les résidus et pôles :

Pour la fonction residuz il facilite le calcul pour trouver les coefficients du système.