Université Amar Telidji de Laghouat

Département d’électronique

Niveau : 3ème Année Lic

Module : TP Traitement du Signal

TP 00 : Introduction à MATLAB

1.1 Qu’es ce que Matlab ?

MATLAB est un logiciel de calcul technique à haut-niveau équipé d'une interface facile à utiliser. Son nom
provient des mots MATrix et LABoratory comme il est basé sur l'utilisation des matrices. MATLAB est a
le produit de Mathworks, et la première version de MATLAB été sortie en 1984. Il a été adopté la première fois
par la communautés des ingénieurs d’automatique et des mathématiques appliquées, mais comme les années ont
passé de nombreuses fonctions (classé dans des boîtes à outils toolbox) ont été incluses ce qui a rendu MATLAB
un outil extrêmement puissant utile pour les scientifiques et les ingénieurs de diverses disciplines. Par exemple,
MATLAB peut être employé dans un éventail d’applications, telles que les télécommunications, traitement de
signal et d'image, commande, mathématiques, modélisation financières, biologie, aéronautique, et beaucoup plus.
Il est également très utiles et populaires dans l'éducation parce qu'un étudiant peut facilement dériver des résultats
maniables sans faire face à beaucoup de difficultés, comme langage de programmation la syntaxe et la structure de
Matlab est très simple et l'analyse de données et la visualisation est facile. Avec l'aide de MATLAB dans ce
manuel on s'attend à ce que le l’étudiant comprenne de manière pratique les concepts difficile du signal et de la
théorie des systèmes.
1.2 Prise en main :
Dans cette partie, on va illustrer comment effectuer quelques opérations de base sur MATLAB. La présentation
est faite sous forme de tables. Dans la partie gauche on donne l'instruction de MATLAB, alors que dans le côté
droit, les résultats et leurs explications.
1.2.1 Opérations Arithmétiques Simples :
Les quatre opérations {+, -, *, et / } .sont employés comme d'habitude pour additionner, soustraient, multiplier, et
diviser deux nombres. Dans MATLAB un opérateur gauche de division {\} est également disponible, et la
puissance d'un nombre est calculé par l'opérateur {^}.

Commandes Résultat / Commentaires

3+5 ans = 8 %Addition
5-2 ans = 3 %Soustraction
6*7 ans = 42 %Multiplication
2/4 ans = 0.5 %Division
2\4 ans = 2 %Division à Gauche
2^3 ans = 8 %Power
1.2.2 Commentaires
Le symbole % est utilisé pour insérer des commentaires. Quoi que ce soit écrit à droite de % est vu comme
commentaire et n'est pas pris en considération par MATLAB.
1.3 Gestion de la Mémoire
Les variables sont sauvegardées dans une région de mémoire dite « workspace ». Pour voir les variables
sauvegardées la commande (>> workspace) est utilisée. Alternativement, on peut employer les commandes
(>>who) et (>>whos).
who Your variable are: a b s % Les variables sauvegardées
whos Name Size Bytes Class
a 1x1 8 double array
b 1x1 8 double array
s 1x5 10 char array
1.4 Vecteurs
MATLAB est un outil de programmation spécialisé pour travailler avec des matrices (ou des tableaux). Même les
scalaires sont considéré comme des matrices de taille 1×1. Dans cette partie, on présente un cas spécial de
matrices, les vecteurs.

M. REGGAB Page 1
2015/2016
 1.4.1 Vecteurs Ligne
Pour définir un vecteur ligne, les éléments du vecteur sont donnés entre crochets [... ]. Des espaces ou des
virgules sont insérés entre ces éléments. Un élément de vecteur est indiqué par son index. L’index du premier
élément d'un vecteur est 1. Des parenthèses sont utilisées afin de se rapporter à une position dans un vecteur.

Commandes Résultat / Commentaires

A=[1 2 3 4 5] A=12345 % vecteur ligne A de 5 éléments.
B=[-6,-7,-8,-9] B = -6 -7 -8 -9 % vecteur ligne B de 4 éléments.
C=[0.3,8.5] C = 0.3 8.5 % vecteur ligne C de 2 éléments décimaux.
B(3) ans = -8 % Le 3ieme élément du vecteur B
B(6) ???index exceeds matrix dimensions
%le Vecteur B a 4 éléments. donc, se rapporter à sa 6eme position causera une erreur.
1.4.2 Addition/Soustraction
Les opérations entre les vecteurs sont généralement effectuées entre les éléments correspondants des vecteurs, c.-
à-d. on a des opérations par élément. Une condition nécessaire pour l'opération entre deux vecteurs est que les
vecteurs doivent être de la même taille.
 Si on suppose que a = [a1,a2,…an] et b = [ b1,b2,…,bn]. Alors a+b = [a1+ b1,a2+b2,…an+ bn] et
 a-b = [a1- b1,a2-b2,…an- bn]
1.4.3 Multiplication, Division, et puissance
Pour ce type d'opérations, on doit faire attention au fait que c'est une opération élément par élément. Pour exécuter
la multiplication et la division entre deux vecteurs ou pour calculer la puissance d'un vecteur l'opérateur point "."
doit être inséré avant l'opérateur.
 Si on suppose que a = [a1,a2,…an] et b = [ b1,b2,…,bn]. Alors

Commandes Résultat / Commentaires

a=1 :2 :11 a = 1 3 5 7 9 11
b=1 :6 b=123456 %Les vecteurs a et b.
a.*b ans = 1 6 15 28 45 66 %Multiplication élément par élément.
a./b ans = 1.00 1.50 1.67 1.75 1.80 1.833 %Division élément par élément.
a.\b ans = 1.00 0.67 0.60 0.57 0.56 0.54 % Division à gauche élément par élément.
a.^2 ans = 1 9 25 49 81 121 % la puissance carrée de chaque élément de a
2.^a ans = 2 8 32 128 512 2048 % 2 a la puissance de chaque élément de a
a.^b ans = 1 9 125 2401 59049 1771516 %Chaque élément de a a la puissance de l’élément
correspondent de b

1.5 Vecteurs de Colonne

Jusqu'ici les vecteurs considérés étaient des vectrices lignes. La manière de définir un vecteur de colonne est tout à
fait semblable. Dans ce cas, les éléments également sont insérés entre crochets. Cependant, les éléments ne sont
pas séparés par des virgule mais par des points-virgules ";".d’autres manières pour avoir un vecteur de colonne sont
démontrées dans le tableau suivent :
Commandes Résultat / Commentaires
a = [1;2;3;4] a=1 %méthode 1 les points-virgules sont utilisés
2 %pour aller a la ligne suivante.
3
4
b = [5 b=5 % Méthode 2 utilisé entrée
6 6 %pour aller a la ligne suivante

M. REGGAB Page 2
2015/2016
 7 7
8] 8
c = [9 10 11 12]
d = c’ d = 9 % méthode 3 utilisé l’operateur de transposition pour avoir
10 %un vecteur colonne a partir d’un vecteur ligne
11
12
z = [1:0.1:100]’; % créer un vecteur colonne qui commence par 1 et se termine par 100 avec un pas de 0.1
size(z) ans = 991 1 %la taille du vecteur z
length(z) ans = 991 % l’instruction length donne le nombre d’éléments du vecteur

1.5.1 produit scalaire de deux vecteur :

Il est nécessaire que la taille des deux vecteurs soit la même, soit a = [a 1,a2,…an] et b = [ b1,b2,…,bn]. le produit

scalaire de a et b est donné par il existe deux méthodes pour calculer le produit scalaire de deux
vecteurs (voir le tableau)
Commandes Résultat / Commentaires
a=1 :4 a=1 2 3 4 % Définition des vecteurs a et b
b=2 :5 b=2 3 4 5
ps1=dot(a,b) ps1=40 %méthode 1 utilisation de la commande dot
ps2=sum(a.*b) ps2=40 % méthode 2 calcule par la définition
1.6 Commandes utiles :
Commandes Résultat / Commentaires
a = [4 2 7 0 6] a=42706 %Définition du vecteur a
s = sum(a) S=19 %la somme des élément de a
c = cumsum(a) c = 4 6 13 13 19 %La somme cumulé des éléments de a
p = prod(a) p=0 %Le produit des éléments de a
d = diff(a) d = -2 5 -7 6 %La différence entre deux élément successives
[m,i] = max(a) m=7 %la valeur maximale de a et stocké dans m
i=4 %et son indice est stocké dans i
[m,i] = min(a) m=0 %la valeur minimale de a et stocké dans m
i=3 %et son indice est stocké dans i
max(a) ans=7 % si on ne donne de variable de sortie on obtient seulement la valeur max de a
mean(a) ans = 3.80 % donne la moyenne des élément de a
sort(a) ans = 0 2 4 6 7 %trier les éléments de a dans l’ordre croissant
sort(a,’descend’) ans = 7 6 4 2 0 %trier les éléments de a dans l’ordre décroissant
1.7 Les Matrices :
Il existe trois méthodes pour créer une matrice sous Matlab : (voir le tableau)
Commandes Résultat / Commentaires
A = [1 2 3; 4 5 6; 7 8 9] A=123
456
7 8 9 %Méthode 1
A = [1 2 3 A=123
456 456
7 8 9] 7 8 9 %Méthode 2
D = [0:.1:.5;0:-.1:-.5;0:pi/5:pi] D = 0 0.1 0.2 0.3 0.4 0.5
0 -0.1 -0.2 -0.3 -0.4 -0.5
0 0.63 1.26 1.88 2.51 3.14 %Méthode 3
1.8 Travailler sur des matrices :
On utilise deux indices pour désigner un élément d’une matrice, le premier désigne la ligne de la matrice et le
second sa colonne .on utilise l’operateur « : » pour désigner des sous-matrices.
Commandes Résultat / Commentaires
A = [1 2 3 ; 4 5 6 ; 7 8 9] A = 1 2 3 %La Matrice A(3,3)
456
789
A(: , :) ans = 1 2 3 %Affichage de tout les lignes et toutes les colonnes de A
M. REGGAB Page 3
2015/2016
 456
789
A(2 , :) ans = 4 5 6 %Affichage de la 2eme ligne et toutes les colonnes de A
A(: , 3) Ans =3 %Affichage de la 3eme colonne et toutes les lignes de A
6
9
A(1 : 2 , 1 : 2) ans = 1 2 %affichage de la sous matrice de A composée de la 1ere et 2eme
4 5 %ligne et la 1ere et 2eme colonne de A
A(1 : 2 : 3 , 1 : 2 : 3) ans = 1 3 %affichage de la sous matrice de A composée de la 1ere et 3eme
7 9 %ligne et la 1ere et 3eme colonne de A
A(1 : 2 : 3, 1 : 2) ans = 1 2 %affichage de la sous matrice de A composée de la 1ere et 3eme
7 8 %ligne et la 1ere et 2eme colonne de A
A( 1: 2 :3 , :) ans = 1 2 3 %affichage de la sous matrice de A composée de la 1ere et 3eme
7 8 9 % ligne et toutes les colonnes de A
1.8.1 Addition/Soustraction de deux matrices :
Les matrices doivent avoir la même taille pour pouvoir faire ces opérations.
Commandes Résultat / Commentaires
A = [1 2 3 ; 4 5 6 ; 7 8 9] A=123
456
789
B = [1 1 1 ; 3 3 3 ; 5 5 5] B=111
333
555
C = A+B C = 2 3 4 %addition
789
12 13 14
D = A-B D = 0 1 2 %soustraction
123
234
1.8.2 Multiplication de deux matrices :
Il nécessaire que le nombre de colonnes de la 1ere matrice soit egale au nombre de lignes de la 2eme matrice pour

pouvoir faire la multiplication c.-à-d. A(N ,M)*B(M,K)=C(N,K) Avec et i=1..N et j=1..K

Commandes Résultat / Commentaires
A = [1 2 3 ; 4 5 6] A=123
456 %A(2,3)
B = [1 1 1 1; 2 2 2 2; 3 3 3 3] B=1111
2222
3 3 3 3 %B(3,4)
C = A*B C = 14 14 14 14
32 32 32 32 %C(2,4)
D = B*A ???Error using =>* Inner matrix dimensions must agree.
% Le produit matriciel n’est pas commutatif

1.8.3 Inverse d’une matrice

Commandes Résultat / Commentaires
A = [2 4 6 ; 1 1 1 ; 3 4 1] A=246
1 1 1 %Matrice Carrée A
341
B = inv(A) B = -0.3750 2.5000 -0.2500
0.2500 -2.0000 0.5000
0.1250 0.5000 -0.2500 %matrice inverse de A
C = A^-1 C = -0.3750 2.5000 -0.2500
0.2500 -2.0000 0.5000
0.1250 0.5000 -0.2500 %2eme Methode A-1
A*B ans = 1.0000 0.0000 0.0000
0.0000 1.0000 0.0000 %A*B=I
0.0000 0.0000 1.0000
M. REGGAB Page 4
2015/2016
 B*A ans = 1.0000 0.0000 0.0000
0.0000 1.0000 0.0000 %B*A=I
0.0000 0.0000 1.0000
1.8.4 Déterminant d’une matrice :
Pour calculer l’inverse d’une matrice elle doit être carrée et son déterminant non nul, sous matlab on peut calculer
le déterminant d’une matrice
Commandes Résultat / Commentaires
A = [1 2 3 ; 4 5 6] A=123
456 %A(2,3)
Inv(A) ???Error using=>inv Matrix must be square.
% La matrice A doit être carrée
A = [2 4 6 ; 1 1 1 ; 3 4 1] A=246
111 %A(3,3)
341
det(A) ans = 8 %le déterminant pas nul A est invertible
inv(A) ans = -0.3750 2.5000 -0.2500
0.2500 -2.0000 0.5000 %L’inverse de A
0.1250 0.5000 -0.2500
B = [1 1 1; 2 2 2; 3 3 3] B=111
222 %B(3,3)
333
det(B) ans=0 %le déterminant est nul B n’est pas invertible
inv(B) ans = Inf Inf Inf
Inf Inf Inf
Inf Inf Inf
1.8.5 Matrices Spéciales :
Commandes Résultat / Commentaires
ones(2,3) ans = 1 1 1
111 %Matrice de taille 2x3 remplis de 1.
zeros(1,4) ans = 0 0 0 0 %Matrice de taille 1x4(ou vecteur) remplis de 0.
rand(3) ans = 0.8600 0.9000 0.4600
0.6000 0.1000 -0.3000 %Matrice de taille 3x3 remplis de valeurs aléatoires.
0.4540 0.6023 0.2700 %entre 0 et 1
eye(4,2) ans = 1 0
01 %Matrice de taille 4x2 avec des 1 dans la diagonale parincipale
00
00
eye(3) Ans = 1 0 0
010 %Matrice identité de taille 3x3
001
A = magic(3) A=816
3 5 7 %Matrice magique
492

Travail à faire :
1. créer un vecteur a=[0,0.1,0.2,…,10] et un vecteur b=[cos(0),cos(0.2),cos(0.4),….,cos(20)] puis calculer les valeurs
suivantes.
a) c=a/b
b) d=a4
c) le produit de a et b
d) le produit scalaire de a et b
2. créer une matrice A(3,3) dont les element sont des nombre aléatoire uniformément distribué dans l’intervalle [0,10]
et calculer son
a) inverse
b) transposé
c) déterminant
d) taille

M. REGGAB Page 5
2015/2016
3. soit les deux vecteurs add=a+b=[13,12,11,10,10] et sub=a-b=[-11,-6,-14,8] trouver les vecteurs a et b

M. REGGAB Page 6
2015/2016