Séance de TP 2

Séance de TP 2
Exercice 1
%Le code MATLAB suivant permet de générer et de tracer les signaux discrets basiques.
subplot(3,3,1);
stem([1;zeros(49,1)]);
title('impulsion de Dirac')
subplot(3,3,2); stem(ones(50,1));
title('signal échelon unitaire')
subplot(3,3,3);
stem([ones(1,5),zeros(1,3)])
title('signal Rectangulaire')
subplot(3,3,4);
stem(sin(2*pi/8*(0:15)))
title('signal Sinusoidal')
subplot(3,3,5); stem(sinc(0:0.25:8));
title('"Sinc" signal')
subplot(3,3,6); stem(exp(-(0:15)));
title('e^-n signal')
subplot(3,3,7);
stem(pow2(-0.5*(0:15)))
title('2^-(0.5n) signal')
subplot(3,3,8); stem(3.^(0:15));
title('3^n signal')
subplot(3,3,9); stem(randn(1,16));
title('Gaussian random signal')
Exercice 2 Génération de signaux
Rappel sur les signaux numériques
Un signal numérique est défini par un nombre d'échantillons N

relevés à une fréquence d'échantillonnage Fe.
On s’intéresse au domaines temporel et fréquentiel
les axes temporel et fréquentiel sont définis comme suit:
L'axe temporel est un vecteur de N points espacés de 1/Fe.

Sous Matlab, la syntaxe est : axe_temps = (0:N-1)/Fe;
L'axe fréquentiel est un vecteur de N points compris entre 0

et Fe. La syntaxe est : axe_freq = (0:N-1)*Fe/N;
% tracé d’un signal sinusoïdal continu et discret
Fs=16e3; t=0:1/Fs:5e-3;
n=0:length(t)-1;
length(t)%81
Lastt=5e-3*Fs%80
subplot(211);
plot(t,10*sin(2*pi*1200*t+pi/4));
xlabel('continuous time');
ylabel('amplitude')
title('Continuous-time sinusoidal signal')
subplot(212);
stem(10*sin(2*pi*1200/16000*n+pi/4))
xlabel('discrete time');
ylabel('amplitude')
title('Discrete-time sinusoidal signal')
Exercice 3
Tracé d’un signal sinusoïdal à fréquence variable
cos(2*pi* f*t)
dont la fréquence f varie avec le temps selon la loi
f = 0.01 t.
On effectuera le tracé sur une centaine de points

code
freq = 0.01 * [0:100];
theta = 2*pi*freq.*[0:100];
s = cos (theta);
plot ([0:100],s)
Exercice 4 % signal modulé AM:
n = linspace(0,1200,300); m = 0.7;
fH = 0.8 ; fL = 0.003;
xH = sin(2*pi*fH*n);
xL = sin(2*pi*fL*n);
y = (1+m*xL).*xH;
subplot(211)
stem(n,y); grid ;
xlabel('Time index n');
ylabel('Amplitude');
subplot(212)
plot(n,y); grid ;
Exercice 5
Modulation d’amplitude:
On prend un signal m(t) égal à la somme de 3 sinusoïdes de
Fréquences respectives 2130Hz, 3750Hz et 4960Hz
Et d’amplitudes respectives 1, 1.8 et 0.9
On choisit fo=50kHz. On prend k=0.25(indice de modulation)
écrire un programme qui engendre x(t) sur une durée de 2 ms.
La freq. d’échantillonnage fe est 500KHz

fe=500000; T = 2/1000;
N= T*fe;%T=N*Te=N/fe
n=(0:N); t = n/fe; k=1/4
f0=50000;f0r=f0/fe
fm=[2130;3750;4960];fmr=fm/fe;
Am=[1,1.8,0.9]
m=Am*cos(2*pi*fmr*n);
p=cos(2*pi*f0r*n);
x=(1+k*m).*p;
subplot(311); plot(t,m)
subplot(312); plot (t,x)
Lf=2048; f=(0:Lf)/Lf;
g=fft(x,Lf);
mod=abs(g);
subplot(313); plot(f(1:Lf/2),mod(1:Lf/2))
Exercice 6
Modulation de fréquence par une rampe
𝒕
s(t) = A cos theta 𝒔 ( 𝒕 )= 𝑨𝒄𝒐𝒔 (𝟐 𝝅 𝒇 𝟎 𝒕+𝟐 𝝅 𝒌 𝒇 ∫ 𝒎 𝒇 ( 𝒖 ) 𝒅𝒖)

𝟎
On prend pour signal modulant m(t) = t, on en déduit

theta(t) = 2 p (F0t+kf . t2/2)
On choisit F0 = 2000 Hz, une durée de de signal de 2s et une fréquence

d’échantillonnage fe=8000Hz.
Ecrire un programme qui demande la valeur de kf et qui Permet d’écouter par

la fonction sound le signal modulé obtenu.
Prendre kf = 1000 puis kf = 1500
Qu’observe-t-on?
fe=8000;F0=2000; T=2;
%disp('frequence d''echantillonnage : 8000 Hz’);
%disp('freq.instantanee:fi(t)=2000+kf * t (Hz)’);
%disp('duree de modulation : 2s');
kf=input('choisir kf (Hz/s) : ');
%tt=sprintf('Excursion en freq. de%iHz a %iHz', F0,
F0+kf*T);
%disp(tt)
t=(0:fe*T-1)/fe; theta=2*pi*(F0*t+0.5*kf*t.^2);
x=cos(theta);axis('auto'); sound(x,8000)
fi(t) = 2000+kf * t t varie de 0 à 2s
- Pour Kf =1000, la fréquence instantanée varie de 2000 Hz

(grave) à 4000 Hz (aigu).
- Pour Kf = 1500, fi(t) varie de 2000 Hz à 5000 Hz,
Mais comme la fréquence d’échantillonnage est de 8000 Hz, on a

du repliement de spectre :
les fréquences entre 4000 Hz et 5000 Hz sont perçues comme
allant de 4000 Hz à 3000 Hz.
Ce qui explique la variation de l’aigu vers le grave en fin

d’écoute.

Séance de TP 2

Transféré par

Droits d'auteur :

Formats disponibles

Séance de TP 2

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Séance de TP 2

Transféré par

Droits d'auteur :

Formats disponibles

Séance de TP 2

Rappel sur les signaux numériques

Un signal numérique est défini par un nombre d'échantillons N

L'axe temporel est un vecteur de N points espacés de 1/Fe.

L'axe fréquentiel est un vecteur de N points compris entre 0

dont la fréquence f varie avec le temps selon la loi

On effectuera le tracé sur une centaine de points

freq = 0.01 * [0:100];

On prend un signal m(t) égal à la somme de 3 sinusoïdes de

Fréquences respectives 2130Hz, 3750Hz et 4960Hz

Et d’amplitudes respectives 1, 1.8 et 0.9

On choisit fo=50kHz. On prend k=0.25(indice de modulation)

écrire un programme qui engendre x(t) sur une durée de 2 ms.

La freq. d’échantillonnage fe est 500KHz

s(t) = A cos theta 𝒔 ( 𝒕 )= 𝑨𝒄𝒐𝒔 (𝟐 𝝅 𝒇 𝟎 𝒕+𝟐 𝝅 𝒌 𝒇 ∫ 𝒎 𝒇 ( 𝒖 ) 𝒅𝒖)

On prend pour signal modulant m(t) = t, on en déduit

On choisit F0 = 2000 Hz, une durée de de signal de 2s et une fréquence

Ecrire un programme qui demande la valeur de kf et qui Permet d’écouter par

- Pour Kf =1000, la fréquence instantanée varie de 2000 Hz

Mais comme la fréquence d’échantillonnage est de 8000 Hz, on a

Ce qui explique la variation de l’aigu vers le grave en fin

Vous aimerez peut-être aussi