TP TRAITEMENT AVANCÉ DU SIGNAL SOUS MATLAB UEM 1.1 RESPONSABLE A.

ALLICHE 2016/2017 TP TRAITEMENT AVANCÉ DU SIGNAL SOUS MATLAB UEM 1.1 RESPONSABLE A.ALLICHE 2016/2017
TP02 : ETUDE DES SYSTEMES A TEMPS DISCRETS
TRAVAIL PREPARATOIRE : Soit la table de TZ 1°) Développement en Partie Fractionnaire sous MATLAB
No. Signal Transformée z- Région de La fonction MATLAB ‘residuez’ peut être utilisée pour déterminer le développement en partie fractionnaire (DPF)
Convergence de la TZ. La syntaxe est donnée par :
[R, P, K] = residuez(B, A)
1 Les vecteurs B et A sont les coefficients du numérateur et du dénominateur des polynômes et
d , respectivement
2 1 Plan z
3
Exemple 1
4
Donner le développement en partie fractionnaire des TZ suivantes:
5 (a)
(b)
6
(c)
7
Solution
8
Partie(a) : Dans ce cas on a :
9
10 En multipliant chacun des numérateur et du dénominateur par on obtient,
11
Ceci peut être ecrit aussi :
12
13
En utilisant la fonction de MATLAB ‘conv’, on aurra
14 >> conv([1 -1],[1 -0.5])
ans =
où et 1.0000 -1.5000 0.5000
MATLAB
sont des constantes complexes définies
15 par Dans ce cas la correspondance est : >> conv([1 -1],[1 -0.5])
, ans =
1.0000 -1.5000 0.5000
1 Avec la commande de MATLAB ‘residuez’ on détermine les coefficients 2
Master Reseaux & Télécommunications Université A.MIRA Bejaia Faculté de Technologie Département d’Electronique Master Reseaux & Télécommunications Université A.MIRA Bejaia Faculté de Technologie Département d’Electronique
 TP TRAITEMENT AVANCÉ DU SIGNAL SOUS MATLAB UEM 1.1 RESPONSABLE A.ALLICHE 2016/2017 TP TRAITEMENT AVANCÉ DU SIGNAL SOUS MATLAB UEM 1.1 RESPONSABLE A.ALLICHE 2016/2017
Ou
MATLAB
>> [R, P, K] = residuez([1 0],[1 -1.5 0.5])
R= Part(c) :
2
-1
P=
Alors, 1.0000
>> conv(conv([1 -1],[1 -0.5]),[1 -0.5])
0.5000
ans =
K=
Ou 1.0000 -2.0000 1.2500 -0.2500 MATLAB
[]
>> conv([1 -1],[1 -1 0.5])
ans =
Partie(b) : 1.0000 -2.0000 1.2500 -0.2500
MATLAB
>> [R, P, K] = residue([1 0],[1 -2 1.25 -0.25])
>> conv([1 -1],[1 -1 0.5])
R=
ans =
MATLAB
1.0000 -2.0000 1.5000 -0.5000 4.0000 -4.0000 -1.0000
>> conv([1 -1],[1 -1 0.5])
P=
ans =
1.0000 0.5000 0.5000
1.0000 -2.0000 1.5000 -0.5000
K=
[]
MATLAB
>> [R, P, K] = residuez([1 1 0],[1 -2 1.5 -0.5]) 2°) Travail à réaliser :
R= A) Lire help les fichiers de Matlab fonction residuez ( ), roots, zplane et filter ( ).
Bien lire, et comprendre la procédure de calcul de la TZ inverse par la méthode précédente
4.0000 -1.5000 - 0.5000i -1.5000 + 0.5000i
P = 1.0000 0.5000 + 0.5000i 0.5000 - 0.5000i
K=
[]
Calculez et représenter la réponse impulsionnelle x(n) en utilisant le DPF et en utilisant la commande
filter, donner la représentation dans le plan z ,
3 4
Master Reseaux & Télécommunications Université A.MIRA Bejaia Faculté de Technologie Département d’Electronique Master Reseaux & Télécommunications Université A.MIRA Bejaia Faculté de Technologie Département d’Electronique
 TP TRAITEMENT AVANCÉ DU SIGNAL SOUS MATLAB UEM 1.1 RESPONSABLE A.ALLICHE 2016/2017 TP TRAITEMENT AVANCÉ DU SIGNAL SOUS MATLAB UEM 1.1 RESPONSABLE A.ALLICHE 2016/2017
B) Soit l’équation aux différences d’’un système LI à Temps discret : TP1 : ETUDE DES SIGNAUX DISCRETS & ANALYSE SPECTRALE (FFT)
En plus d’appliquer les notions essentielles du traitement du signal, ce TP est également une initiation à Matlab :
TRAVAIL PREPARATOIRE :
1) Déterminer manuellement la fonction de transfert H(z) 1°) Rappeler la transformée de Fourier de rect( t ). En déduire celle de rect( t/T ) puis celle de x(t)
2) Tracer la réponse impulsionnelle en utilisant les deux méthodes = rect((t-t0)/T ). On posera TF[ x(t) ] = X(f).
3) On considère l’entrée du système x(n) =U(n)- U(n-100) , donnez la représentation temporelle de x(n)
2°) Calculez la TFtD et la TFD de rect(n)|N donnez la représentation graphique pour N=6
et de y(n) pour n=[0 :200], en se basant sur la sortie est-ce un filtre passe Haut ou passe Bas
4) Appliquez au fichier music.wav, écouter la sortie du signal filtré et quelle est la nature du filtre ? Travail à réaliser :
A) Ecrire un programme sous Matlab permettant le calcul de la TFTD et de la TFD et inverse d’un
signal à temps discret, appliquez x= [1 1 1 1} comparer avec la commande fft et ifft de Matlab ?
B) Chargement des données :
x Le signal étudié est dans le fichier de donnée sig1_rect.
x Le chargement est réalisé par la commande load sig1_rect.
x Dans l’espace de travail (workspace) sont générées deux variables : Fs la fréquence
d’échantillonnage du signal et xk le tableau de points résultant de l’échantillonnage.
Visualisation du signal :
Avec la commande stem, représenter graphiquement le signal en étalonnant correctement l’axe
des abscisses. En déduire T, Ts et T. La condition de Shannon a-t-elle été respectée ?
Transformée de Fourier discrète :
x Avec la commande fft, calculer Xn la TFD de xk.
x Représenter graphiquement son module en étalonnant correctement l’axe des abscisses en
fréquences. Nous pouvons nous contenter de la moitié des points : pourquoi ?
x Comparer à la valeur théorique.
Interpolation de la TFD :
Nous ne disposons pas de suffisamment de points dans l’estimation de X(f). Pour cela nous allons
procéder à la manipulation suivante :
x Constituer un nouveau signal xk1 à partir de xk auquel on ajoute 4.N zéros.
x Calculer la TFD Xn1 de ce nouveau signal.
x La représenter sur la même figure que Xn.
x Comparer les résultats et interpréter la manipulation effectuée.
A) ANALYSE spectrale d’un signal par la TFD :
Le signal étudié est contenu dans le fichier sig2_quidonc.
Travail à réaliser :
Chargement des données :
x Le signal étudié est dans le fichier de donnée sig2_quidonc.
x Le chargement est réalisé par la commande load sig2_quidonc.
x Dans l’espace de travail (workspace) sont générées deux variables : Fs la fréquence
d’échantillonnage du signal et xk le tableau de points résultant de l’échantillonnage.
5 Visualisation du signal : 1
Master Reseaux & Télécommunications Université A.MIRA Bejaia Faculté de Technologie Département d’Electronique Master Reseaux & Télécommunications Université A.MIRA Bejaia Faculté de Technologie Département d’Electronique
 TP TRAITEMENT AVANCÉ DU SIGNAL SOUS MATLAB UEM 1.1 RESPONSABLE A.ALLICHE 2016/2017
Avec la commande stem, sur une première figure, représenter le signal en étalonnant
correctement l’axe des abscisses. Quelles indications pouvez vous en retirer ?
Première analyse spectrale :
x Avec la commande fft, calculer la TFD Xn du signal xk.
x Compte tenu de l’étude précédente, construire un signal xk1 constitué du signal xk auquel on
concatènera 4N zéros.
x Calculer Xn1 la TFD de xk1.
x Sur une deuxième figure, représenter simultanément (commande hold on) les modules de Xn et
Xn1 en étalonnant correctement l’axe des abscisses (utiliser des couleurs différentes et
commencer par représenter Xn1 qui, ayant plus de points, risque de masquer Xn).
x Pouvez-vous en déduire une forme analytique du signal étudié ?
Effet d’une fenêtre de pondération :
Avec la commande Hamming représenter graphiquement wham(n) ainsi que son spectre
Tapez help hamning et étudiez les autres types de fenêtres, utilisez wintool, et sptool
Deuxième analyse spectrale
x Créer le signal xkh en pondérant xk par une fenêtre de hanning ou autres (commande hanning)
x Reprendre l’analyse précédente avec ce nouveau fichier soit :
1. Avec la commande fft, calculer la TFD Xnh du signal xkh.
2. Construire un signal xkh1 constitué du signal xkh auquel on concatènera 4N zéros.
3. Calculer Xnh1 la TFD de xkh1.
4. Sur une troisième figure, représenter simultanément les modules de Xnh et Xnh1 en étalonnant
correctement l’axe des abscisses
5. Pouvez-vous en déduire une forme analytique du signal étudié ?
Quelles sont vos conclusions sur l’analyse spectrale à partir de la TFD ?
B) Influence de la période d’échantillonnage
Le signal xk précédent est sous échantillonné pour fournir un signal xk2. La nouvelle période
d’échantillonnage sera Fs1=Fs/4.
Troisième analyse spectrale :
x Créer le signal xk2 en prélevant dans xk un échantillon sur 4 (commande for).
x Créer le signal xkh2 en pondérant xk2 par une fenêtre de hanning (commande hanning).
x Construire un signal xkh3 constitué du signal xkh2 auquel on concatènera 4N zéros.
x Reprendre l’analyse avec ce nouveau fichier soit :
1. Avec la commande fft, calculer la TFD Xnh3 du signal xkh3.
2. Sur une quatrième figure, représenter le module de Xnh3 en étalonnant correctement
l’axe des abscisses
3. Comparer au spectre de la troisième figure et interpréter le résultat.
2
Master Reseaux & Télécommunications Université A.MIRA Bejaia Faculté de Technologie Département d’Electronique