Base des télécommunications
TP 1 –Traitement numérique du signal (Échantillonnage)
L’objectif de ce travail est d’aborder la manipulation de signaux discrets en utilisant le
logiciel MATLAB.
Introduction
Ce TP est l’occasion de se familiariser avec le logiciel MATLAB au traitement du
signal via des exemples simples : représentation de signaux en temps et en fréquence,
application des théorèmes de base du traitement numérique du signal (échantillonnage de
Shannon).
I. Génération d’un signal numérique
Engendrer un signal numérique représentant N échantillons d’une cosinusoïde de fréquence
f 0 = 2 kHz.
1
T e=
Ce signal échantillonné est en fait une suite de N points tous les f e , que l’on peut donc
stocker dans un vecteur de taille N.
f0 = 2000; % La fréquence du signal
fe = …. % La fréquence d’échantillonnage
N = 2*fe/f0; %Le nombre d’échantillons dépend de ces fréquences
t = [Link]N-1]; %Nous générons un vecteur normé
T = t / fe; %Vecteur de temps [Te 2*Te 3*Te ... (N-1)*Te]
Y = cos(2*pi*f0*T); % Y est donc un vecteur de N points représentant un cosinus
échantillonné
plot(T, Y, '-+') % On trace la courbe
title('Echantillonnage d''un cosinus')
xlabel('Temps (s)')
ylabel('Amplitude')
E = ……… %Calcul de l’énergie : Somme des carrés des termes divisée par
le nombre d’échantillons
I.1. Tracer ce signal en faisant apparaître en abscisse le temps et en variant la
fréquence d’échantillonnage fe dans la plage suivante : 3, 4, 8 et 16 kHz.
I.2. Quelle est la valeur de E pour la plage de fe. Qu’est-ce que Vous remarquez ?
Remarque : Pour faire le calcul des sommes de carrés, nous utilisons la propriété du produit
Y .Y '=∑ y(k )²
scalaire d’un vecteur et de sa transposée : k
II. Théorème de Shannon
Tracer les spectres des signaux numériques sinusoïdaux de fréquence : 4 kHz et 7 kHz.
Justifier les résultats obtenus.
Nous rédigeons un script MATLAB qui va nous permettre de comparer les deux courbes :
f1 = 4000;
f2 = 7000;
fe = 12000;
N1 = 2*fe/f1
N2 = 2*fe/f2
t1 = [Link]N1-1];
T1 = t1 / fe;
t2 = [Link]N2-1];
T2 = t2 / fe;
Y1 = sin(2*%pi*f1*T1);
Y2 = sin(2*%pi*f2*T2);
figure
plot((t1)*fe/N1,abs(fft(Y1)))
title('Transformée de Fourrier d''une sinusoïde de fréquence 4kHz')
xlabel('Fréquence (Hz)')
ylabel('Amplitude')
figure
plot((t2)*fe/N2,abs(fft(Y2)), 'r')
title('Transformée de Fourrier d''une sinusoïde de fréquence 7kHz')
xlabel('Fréquence (Hz)')
ylabel('Amplitude')
Sur la figure 1 :
II.1 Tracer le module du spectre correspondant en fonction de la fréquence.
II.2 Observer et justifier le spectre obtenu. En particulier, préciser l’incrément
fréquentiel et montrer l’influence du nombre de points dans la déduction de la
fréquence du signal.
II.3 Remarquer et justifier l’amplitude des raies, la largueur de la raie ; déterminer
l’énergie de ce spectre. Conclure.
Sur la figure 2 :
II.4 Que se passe-t-il ? Comment justifier cette figure ? Comment éviter ce
phénomène ?
Chargé du module : Abdi et Benbakreti
