TP 1 : Initiation à Matlab et Génération des signaux
FIA : S5
Enseignant : Y. BEN YOUSSEF
AU 2022/2023
1 Objectifs
Le but de ce TP est de vous familiariser avec le logiciel Matlab qui sera utilisé pour tous
les TP de traitement de signal et la génération des signaux : échelon, signe, rectangulaire,
impulsion, triangulaire, rampe, sinusoïdal et sinus cardinal.
Matlab (MATrix LABoratory) est un logiciel de calcul numérique, qui permet de traiter des
problèmes complexes dans de nombreux domaines scientifiques, d’exploiter et de visualiser
des donnèes en 1D/2D/3D. Le langage MATLAB a été conçu par Cleve Moler à la fin des
années 1970 à partir des bibliothèques Fortran, LINPACK et EISPACK2. Il est utilisé par les
scientifiques et les ingénieurs, pour réaliser les différentes tâches de traitement du signal.
De nombreuses fonctions sont disponibles dans la boîte à outils Signal Processing (Sptools).
Cette boîte à outils permet donc d’analyser et de comparer des signaux dans les domaines
temporel, fréquentiel et temps-fréquence, identifier des motifs et des tendances, extraire des
caractéristiques, et développer et valider des algorithmes personnalisés afin d’obtenir des
informations sur des données.
2 Initiation sur Matlab
2.1 Prise en main de Matlab
— Ouvrir Matlab et créer un dossier et nomer le TPtraitementsignal
— Ouvrir le script nomer le TP1 l’estension est .m et le sauvgrader dans votre dossier.
— Taper la commande suivante
u=ones(1,5)
u est un vecteur de valeur 1 : c’est une matrice 1 ligne et 5 colones
— Taper v=zeros(1,5)
v est un vecteur de valeur 0 : c’est une matrice 1 ligne et 5 colones
— Taper le produit
u ∗ v, u. ∗ v, u ∗ v ′ où v ′ est le transposé du vecteur v qui le transforme en un vecteur 5
lignes et 1 colone
3 Visualisation des signaux numériques
Nous allons voir les fonctions qui permettent de générer les signaux numériques usuels.
Lancez Matlab. Placez-vous sous votre répertoire de travail. Créez et nomez un répertoire
puis placez-vous sous ce répertoire.
AU 2022/2023 Page-1/4- TP Traitement du signal sous Matlab
3.1 Signal sinusoidal
t=-2 :0,01 :2 ;
f=50 ;
A1=1 x1 = A1 ∗ cos(2 ∗ pi ∗ f ∗ t);
plot(t,x1,’o’)
title(’Signal Sinusoidal’)
xlabel(’t en secondes’)
ylabel(’V en volts’)
— Ajouter des signaux de différentes amplitudes mais de fréquences 2f et 3f,....
3.2 Signal sinuscardinal
t=0 :pi/10 :2*pi ;
y=sinc(t) ;
plot(t,y,’o’)
title(’Mon premier graphique sous Matlab’)
xlabel(’t en secondes’)
ylabel(’signal sinus Cardinal’)
Remplacer plot par stem.
3.3 L’impulsion unité de Dirac
Soit le programme suivant :
t=-10 :10 ;
x=[zeros(1,10),1,zeros(1,10)] ;
stem(t,x) ;
axis([-10 10 -0.5 1.5]) ;
title() ;
xlabel(’n’) ;
ylabel(’Amplitude’) ;
— Taper et tester le programme.
— Définir le type de signal après l’exécution de programme.
— Comment doit proceder pour modifier stem par plot ?
— Faire les changements nécessaires au programme ci-dessus pour représenter les im-
pulsions suivantes : 5δ(t), 4δ(t + 4), et δ(t − 3).
3.4 Echelon unité
Echelon unité
t=-10 :10 ;
x=[zeros(1,10),ones(1,11)] ;
stem(t,x) ;
axis([-10 10 -0.5 1.5]) ;
AU 2022/2023 Page-2/4- TP Traitement du signal sous Matlab
title(’Echelon unite’) ;
xlabel(’n’) ;
ylabel(’Amplitude’) ;
Modifier stem par plot et noter la différence
3.5 Signal Signe
t=-10 :.001 :10 ;
sgn=sign(t) ;
plot(t,sgn), axis([-10 10,-2 2]), grid title(’Signal signe’) ;
xlabel(’t’) ;
ylabel(’signal signe’) ;
3.6 Signal Rampe
t=-5 : 0.1 : 5 ;
for j=1 :length(t) ;
if t(j)>0
r(j)=t(j) ;
else
r(j)=0 ;
end
end
figure (2)
plot(t,r) ;
axis([-5 5 -0.04 5]) ;
xlabel(’Temps(sec)’) ;
ylabel(’Amplitude’) ;
3.7 Signal triangle
t=-1 : 0.0001 : 1 ;
w=0.2 ;
x=tripuls(t,w) ;
plot(t,x) ; d=0 :1/3 :1 ;
y=pulstran(t,d,’tripuls’,0.1,-1) ;
plot(t,y) ;
3.8 Signal exponentiel décroissant
t=-10 :10 ;
u=[zeros(1,10),ones(1,11)] ;
x=exp(-0.5*t).*u ;
plot(t,x,’o’) ;
axis([-10 10 -1.5 1.5]) ;
AU 2022/2023 Page-3/4- TP Traitement du signal sous Matlab
title(’Exponentielle retardée’) ;
xlabel(’t’) ;
ylabel(’signal’) ;
3.9 Signal carré
fe=1000 ;
te=1/fe ;
% Définition du signal carré
subplot(2,1,1) ;
t=-5 :0.001 :5 ;
x=square(2*pi*t) ;
plot(t,x) ;
axis([-5,5,-2,2]) ;
xlabel(’temps’) ;
ylabel(’x(t)’) ;
title(’signal carre’) ;
AU 2022/2023 Page-4/4- TP Traitement du signal sous Matlab