TP TNS M1 2023
TP TNS M1 2023
TP TNS M1 2023
la commande stem, sur une première figure, représenter le signal en étalonnant
TP1 : ETUDE DES SIGNAUX DISCRETS & ANALYSE SPECTRALE (FFT) correctement l’axe des abscisses. Quelles indications pouvez vous en retirer ?
En plus d’appliquer les notions essentielles du traitement du signal, ce TP est également une initiation à Matlab : Première analyse spectrale :
TRAVAIL PREPARATOIRE : Avec la commande fft, calculer la TFD Xn du signal xk.
1°) Rappeler la transformée de Fourier de rect( t ). En déduire celle de rect( t/ ) puis celle de x(t) Compte tenu de l’étude précédente, construire un signal xk1 constitué du signal xk auquel on
= rect((t‐t0)/ ). On posera TF[ x(t) ] = X(f).
concatènera 4N zéros.
2°) L’échantillonnage de x(t) dure un temps T > et fourni N échantillons : T = N.Ts. A partir des N
Calculer Xn1 la TFD de xk1.
échantillons {xk} nous calculons la transformée de Fourier discrète {Xn}. Combien d’échantillons comporte
Sur une deuxième figure, représenter simultanément (commande hold on) les modules de Xn et
{Xn} ? Quelle est la relation entre Xn et X(f) ?
Xn1 en étalonnant correctement l’axe des abscisses (utiliser des couleurs différentes et
Travail à réaliser : commencer par représenter Xn1 qui, ayant plus de points, risque de masquer Xn).
Chargement des données : Pouvez‐vous en déduire une forme analytique du signal étudié ?
Le signal étudié est dans le fichier de donnée sig1_rect.
Le chargement est réalisé par la commande load sig1_rect. Effet d’une fenêtre de pondération :
Dans l’espace de travail (workspace) sont générées deux variables : Fs la fréquence Avec la commande Hamming représenter graphiquement wham(n) ainsi que son spectre
d’échantillonnage du signal et xk le tableau de points résultant de l’échantillonnage. Tapez help hamning et étudiez les autres types de fenêtres, utilisez wintool, et sptool
Visualisation du signal :
Deuxième analyse spectrale
Avec la commande stem, représenter graphiquement le signal en étalonnant correctement l’axe
Créer le signal xkh en pondérant xk par une fenêtre de hanning ou autres (commande hanning)
des abscisses. En déduire T, Ts et . La condition de Shannon a‐t‐elle été respectée ?
Reprendre l’analyse précédente avec ce nouveau fichier soit :
Transformée de Fourier discrète : 1. Avec la commande fft, calculer la TFD Xnh du signal xkh.
Avec la commande fft, calculer Xn la TFD de xk. 2. Construire un signal xkh1 constitué du signal xkh auquel on concatènera 4N zéros.
Représenter graphiquement son module en étalonnant correctement l’axe des abscisses en 3. Calculer Xnh1 la TFD de xkh1.
fréquences. Nous pouvons nous contenter de la moitié des points : pourquoi ? 4. Sur une troisième figure, représenter simultanément les modules de Xnh et Xnh1 en étalonnant
Comparer à la valeur théorique. correctement l’axe des abscisses
Interpolation de la TFD : 5. Pouvez‐vous en déduire une forme analytique du signal étudié ?
Quelles sont vos conclusions sur l’analyse spectrale à partir 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 :
B) Influence de la période d’échantillonnage
Constituer un nouveau signal xk1 à partir de xk auquel on ajoute 4.N zéros. Le signal xk précédent est sous échantillonné pour fournir un signal xk2. La nouvelle période
Calculer la TFD Xn1 de ce nouveau signal. d’échantillonnage sera Fs1=Fs/4.
La représenter sur la même figure que Xn. Troisième analyse spectrale :
Comparer les résultats et interpréter la manipulation effectuée.
Créer le signal xk2 en prélevant dans xk un échantillon sur 4 (commande for).
A) ANALYSE spectrale d’un signal par la TFD : Créer le signal xkh2 en pondérant xk2 par une fenêtre de hanning (commande hanning).
Le signal étudié est contenu dans le fichier sig2_quidonc. Construire un signal xkh3 constitué du signal xkh2 auquel on concatènera 4N zéros.
Reprendre l’analyse avec ce nouveau fichier soit :
Travail à réaliser :
1. Avec la commande fft, calculer la TFD Xnh3 du signal xkh3.
Chargement des données : 2. Sur une quatrième figure, représenter le module de Xnh3 en étalonnant correctement
Le signal étudié est dans le fichier de donnée sig2_quidonc. l’axe des abscisses
Le chargement est réalisé par la commande load sig2_quidonc. 3. Comparer au spectre de la troisième figure et interpréter le résultat.
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 : 2
1