Compte Rendu TP 1 TS
Binômes : SELMI Manel / 181831063830
CHENDER Zineb / 171731045092
EXPERIMENTATION
A) clc
clear all
close all
f0=5; %la fréquence de signal sinus
T0=1/f0; %La période
Te=0.01; %le pas
tmax=4*T0; %le temps maximal
t=0:Te:tmax; %le temps
x=sin(2*pi*f0*t); %la fonction ( le signal )
plot(t,x); %la représentation graphique continue
hold on %pour tracer sur la même figure
stem(t,x) % la représentation graphique discrète
B) clc
clear all
close all
t=-[Link]; %le temps
x=[zeros(1,9),1,ones(1,11)]; %le vecteur
stem(t,x); % la représentation graphique discrète
axis([-10 10 -0.5 1.5]); %pour fixer les limites sur les axes des ordonnés et d’abscisses
title ('Impulsion unite'); %nommer le graphe
xlabel('n'); %pour attribuer un texte à l’axe des abscisses
ylabel('Amplitude'); % pour attribuer un texte à l’axe des ordonnées
Partie 2 : REALISATION
Manip 1 : Génération de quelques signaux
1. Représentation de l’impulsion de Dirac :
δ(t) .
% Génération d'une impulsion unité
% Génération d'un vecteur de -10 0 20
n=-10 :20 ;
% Génération de l’impulsion unité
u=[ zero(1,10) 1 zero(1,20) ]
%Tracer le signal généré
Stem(n,u) ;
xlabel (‘ Temps indexé en n ’) ;
ylabel (‘ Amplitude ‘) ;
title (‘ impulsion unité ‘) ;
axis ([ -10 20 0 1.2 ]) ;
3δ(t) .
% Génération d'une impulsion unité
% Génération d'un vecteur de -10 0 20
n=-10 :20 ;
% Génération de l’impulsion unité
u=[ zero(1,10) 3 zero(1,20) ]
%Tracer le signal généré
Stem(n,u) ;
xlabel (‘ Temps indexé en n ’) ;
ylabel (‘ Amplitude ‘) ;
title (‘ impulsion unité ‘) ;
axis ([ -10 20 0 4]) ;
2δ(t + 2) .
% Génération d'une impulsion unité
% Génération d'un vecteur de -10 0 20
n=-8 :18;
% Génération de l’impulsion unité
u=[ zero(1,10) 2 zero(1,20) ]
%Tracer le signal généré
Stem(n,u) ;
xlabel (‘ Temps indexé en n ’) ;
ylabel (‘ Amplitude ‘) ;
title (‘ impulsion unité ‘) ;
axis ([ -10 20 0 3]) ;
δ(t-1).
% Génération d'une impulsion unité
% Génération d'un vecteur de -10 0 20
n=-11 :21 ;
% Génération de l’impulsion unité
u=[ zero(1,10) 1 zero(1,20) ]
%Tracer le signal généré
Stem(n,u) ;
xlabel (‘ Temps indexé en n ’) ;
ylabel (‘ Amplitude ‘) ;
title (‘ impulsion unité ‘) ;
axis ([ -10 20 0 2 ]) ;
δ(t + 3) .
% Génération d'une impulsion unité
% Génération d'un vecteur de -10 0 20
n=-7 :23;
% Génération de l’impulsion unité
u=[ zero(1,10) 1 zero(1,20) ]
%Tracer le signal généré
Stem(n,u) ;
xlabel (‘ Temps indexé en n ’) ;
ylabel (‘ Amplitude ‘) ;
title (‘ impulsion unité ‘) ;
axis ([ -10 20 0 2 ]) ;
2. Génération d’une impulsion rectangulaire :
clc
clear all
close all
t=-1:0.0001:1;
x1=rectpuls(t,0.05);
plot(t,x1) ; axis([-0.1 0.1 -0.2 12]) ; grid ;
xlabel('temps(sec)');
ylabel('amplitude');
title('impulsion rectangulaire');
Modification pour construire une impulsion rectangulaire de durée T=0.08s
et l’amplitude A=2V
clc
clear all
close all
t=-1:0.0001:1;
x1=2*rectpuls(t,0.08);
plot(t,x1) ; axis([-0.1 0.1 -0.2 3]) ; grid ;
xlabel('temps(sec)');
ylabel('amplitude');
title('impulsion rectangulaire');
3.Génération du sinus cardinal : sinc(x)
Soit la fonction y(x)=Sinc(x)=sin(x)/x.
clc
clear all
close all
%traçage de la fonction sinus cardinal
%domaine des valeurs de la variable x
x=-4*pi:pi/100:4*pi;
%valeurs de la fonction
y=sinc(x)
%y=(x==0)+sin(x)./(x+(x==0));
%tracage de la fonction sinus cardinal
plot(x,y)
grid
title('sinus cardinal y=sin(x)/x')
Manip 2 : Produit de convolution de deux signaux analogiques
clc
clear all
close all
x=zeros(1,100);
for i=41:60,x(i)=2;end;
y=zeros(1,100);
for i=31:70,y(i)=3;end;
N=100;
a=-N/2:N/2-1;
b=-N:N-2;
subplot(2,2,1);plot(a,x);axis([-50,50,0,2.5]);grid;%subplot:pour afficher les différents fonctions
dans ces différents sous-figures
title('premier signal');xlabel('t');ylabel('x(t)');%title:pour attribuer un titre au graphique
subplot(2,2,2);plot(a,y);axis([-50,50,0,3.5]);grid;%axis:pour fixer les limites sur les axes des
abscisses et des ordonnées
title('deuxième signal');xlabel('t');ylabel('y(t)');%xlbel et ylabel:pour attribuer un texte à l'axe
des abscisses et des ordonnées
subplot(2,2,4);plot(b,conv(x,y));grid;%plot:pour tracer une représentation graphique continue
axis([-60 60 0 140]);title('convolution');xlabel('t');ylabel('conv(x,y)');
Changement de programme pour avoir le produit de convolution de deux
signaux rectangulaires identiques :
clc
clear all
close all
x=zeros(1,100);%
for i=41:60,x(i)=2;end;%
y=zeros(1,100);%
for i=31:70,y(i)=3;end;%
N=100;%
a=-N/2:N/2-1;%
b=-N:N-2;%
subplot(2,2,1);plot(a,x);axis([-50,50,0,2.5]);grid;%
title('premier signal');xlabel('t');ylabel('x(t)');%
subplot(2,2,2);plot(a,y);axis([-50,50,0,3.5]);grid;%
title('deuxième signal');xlabel('t');ylabel('y(t)');%
subplot(2,2,4);plot(b,conv(x,x));grid;%
axis([-60 60 0 140]);title('convolution');xlabel('t');ylabel('conv(x,x)');