Introduction au calcul scientifique avec Matlab
TP2
Correction
Exercice 1.
Observer l'effet des instructions suivantes.
>> x=rand(1,5) => crée un vecteur ligne de 5 termes de valeurs aléatoires comprises entre 0
et 1
>> mean(x) => calcule la valeur moyenne des termes du vecteur x
5
>> std(x) => calcule l’écart type des termes du vecteur x: std(x) =
√ ∑ ( xi−mean( x ))2
i=1
N
>> median(x) => retourne la valeur médiane des termes de x: si x1<x2<x3<x4<x5 alors x3 est la
valeur médiane
>> sort(x) => trie les termes de x du plus petit au plus grand
>> A=rand(3) => crée une matrice 3x3 par des valeurs aléatoires comprises entre 0 et 1
>> sort(A) => trie chaque colonne de la matrice A de manière ascendante (du plus petit au
plus grand)
>> [B, I]=sort(A) => retourne dans la matrice B la matrice A triée, et la matrice I contient des
indices qui spécifie comment le terme ij de A a été réarrangé dans B : for j = 1:n, B(:,j) =
A(I(:,j),j)
>> sort(A') => voir explication sort(A)
>> max(A) => retourne un vecteur qui contient la valeur maximum de chaque colonne de la
matrice A
>> max(A') => retourne un vecteur qui contient la valeur maximum de chaque ligne de la
matrice A
>> max(max(A)) => retourne la valeur maximum de tous les termes de la matrice A
>> sum(A) => retourne un vecteur contenant la somme des termes de chaque colonne de la
matrice A
>> cumsum(A) => retourne une matrice de la même taille que A où chaque colonne est une
somme cumulative de la colonne de A : si B= cumsum(A), alors pour chaque colonne j ;
B(i+1)j = A(i+1)j + Aij pour i=1, 2
>> prod(A) => rerourne un vecteur contenant le produit des termes de chaque colonne de A
FEKAK Fatima-Ezzahra, ENSA, Fès
CP1, 2019-2020
>> diff(A) => retourne une matrice 2x3 don’t chaque colonne correspond au vecteur
difference de la colonne de A: [A(2:n,:) - A(1:n-1,:)]
>> D=A([1,2],1:3) => D contient la sous matrice de A contenant les lignes 1 et 2 et les 3
colonnes
>> sum(D,1) => retourne un vecteur colonne dont les termes sont la somme des termes de
chaque colonne de D
>> sum(D,2) => retourne un vecteur ligne dont les termes sont la somme des termes de
chaque ligne de D
Exercice 2.
Soit X une matrice 2xn contenant les coordonnées de n points du plan. Comment faire pour
obtenir une matrice où les points sont ordonnés par ordre croissant des abscisses ? (sans
utiliser les boucles)
Réponse :
>> B = sortrows(X')' => la matrice B à la même dimension que A et contient les coordonnées
des n points avec des abscisses croisants
En effet, la fonction sortarrows permet de trier les colonnes d’une matrice comme un groupe.
Vu qu’on souhaite trier les lignes de A alors en fait >> sortrows(A') le résultat est une matrice
nx2. Enfin, transpose le résultat pour obtenir une matrice 2xn
Exemple pour n=10 :
>> X = rand(2,10)
X=
0.0971 0.6948 0.9502 0.4387 0.7655 0.1869 0.4456 0.7094 0.2760 0.6551
0.8235 0.3171 0.0344 0.3816 0.7952 0.4898 0.6463 0.7547 0.6797 0.1626
>> sortrows(X')
ans =
0.0971 0.8235
0.1869 0.4898
0.2760 0.6797
0.4387 0.3816
0.4456 0.6463
0.6551 0.1626
0.6948 0.3171
0.7094 0.7547
0.7655 0.7952
0.9502 0.0344
FEKAK Fatima-Ezzahra, ENSA, Fès
CP1, 2019-2020
>> sortrows(X')'
ans =
0.0971 0.1869 0.2760 0.4387 0.4456 0.6551 0.6948 0.7094 0.7655 0.9502
0.8235 0.4898 0.6797 0.3816 0.6463 0.1626 0.3171 0.7547 0.7952 0.0344
Remarque : Pour trouver la fonction « sortarrows » il faut regarder les fonctions qui ont une
relation avec la fonction principale qui est « sort ». Ces fonctions sont toujours proposer par
le help de la fonction principale
Exercice 3.
1. Soit le vecteur de dimension 8 de composantes : 3.2; 4.8; 3.3; 3.2; 3.1; 4.2; 3.2; 3.3
Construisez le vecteur y = (yi)i=1;…;8 correspondant.
2. Construisez, à l'aide des fonctions vectorielles (sans utiliser les boucles), une fonction
moyenne.m qui a comme argument d’entrée n et qui permet de calculer la suite des
moyennes mobiles:
3. Donner la valeur de
Réponse :
1. >> y= [3.2; 4.8;3.3; 3.2; 3.1; 4.2; 3.2; 3.3]
2. - Créer un fichier moyenne.m
- Y mettre le code suivant :
-------------------------------------------------------------------------------------------------------------
function [yn] = moyenne(n)
%MOYENNE calcule le terme n de la suite des moyennes mobiles du vecteur y
% vecteur y
y= [3.2; 4.8;3.3; 3.2; 3.1; 4.2; 3.2; 3.3];
yn = sum(y(1:n))/n;
end
-------------------------------------------------------------------------------------------------------------
3. y 8 = sum(y(1:8)) = 3.5375
Remarque : étant donné que y contient 8 terme alors y 8 est finalement la moyenne de y,
donc on peut aussi retrouver le même résultat en faisant >> y 8=mean(y)
Exercice 4.
Soit la matrice A :
FEKAK Fatima-Ezzahra, ENSA, Fès
CP1, 2019-2020
Construisez, à l’aide des opérateurs logiques, la fonction find.m qui a pour argument d’entrée
un scalaire k et qui permet de retourner les termes égaux à k ainsi que leur position.
Réponse :
Etant donné de « find » est déjà une fonction de matlab, et afin d’éviter tout conflit, on crée
une fonction appelée « find_A »
- Créer un fichier find_A.m
- Y mettre le code suivant :
-------------------------------------------------------------------------------------------------------------
function [I,J] = find_A(k)
%FIND_A retourne les positions ij des termes égaux à k de la matrice A
% I est le vecteur des indices i et J est le vecteur des indices j
% la matrice A
A=[1 2 3; 2 3 1 ; 3 1 2];
%initialisation
I=zeros(3,1);
J=zeros(3,1);
e=1;
l=1;
for i=1:3
for j=1:3
if A(i,j)==k
I(e)= i;
J(l)= j;
e=e+1;
l=l+1;
end
end
end
-------------------------------------------------------------------------------------------------------------
Pour vérifier que ma fonction est juste, je vais comparer son résultat au résultat de la
fonction « find » de Matlab.
Exemple pour k= 2:
>> [I,J]=find_A(2)
I=
1
2
3
J=
2
1
FEKAK Fatima-Ezzahra, ENSA, Fès
CP1, 2019-2020
3
>> [J,I]=find(A==2)
J=
2
1
3
I=
1
2
3
On retrouve bien les mêmes positions des termes égaux à 2, la seule différence est que le
premier argument de ma fonction find_A est le vecteur des indices i et le deuxième argument
et celui des indices j. Contrairement à la fonction « find » de Matlab pour laquelle le premier
argument est le vecteur des indices j et le deuxième est celui des indices i.
Remarque : Il faut toujours bien lire les notices des fonctions Matlab pour mieux comprendre
et bien utiliser les arguments d’entrée et de sortie de ces fonctions.
Exercice 5.
Reprenez le vecteur y de l’exercice 3.
1. Ecrire une fonction existe.m qui a comme argument un vecteur quelconque z et un
scalaire k et qui permet de retourner la réponse à la question suivante : existe-t-il une
coordonnée du vecteur z inférieur à k.
2. Utiliser la fonction existe.m pour répondre à la question suivante : existe-t-il une
coordonnée du vecteur y inférieure à 3.3
3. Ecrire un script qui permet de construire un vecteur logique z tel que la ième
coordonnée de ce vecteur sera 1 si la ième coordonnée du vecteur y est à l'extérieur
de l'intervalle :
Où σ8 est la racine carrée de la variance d'échantillonnage.
Réponse :
1. - Créer un fichier existe.m
- Y mettre le code suivant :
-------------------------------------------------------------------------------------------------------------
function [r] = existe(z,k)
%EXISTE permet de répondre à la question existe-t-il une coordonnée du vecteur z inférieur
àk
n=size(z);
r='non';
for i=1:n
FEKAK Fatima-Ezzahra, ENSA, Fès
CP1, 2019-2020
if z(i) < k
r = 'oui';
end
end
end
-------------------------------------------------------------------------------------------------------------
2. Exemple pour z=y, k=3.3 et k=1 :
>> r=existe(y,3.3)
r=
'oui'
>> r=existe(y,1)
r=
'non'
3. –Créer un fichier script appelé ici logic.m
- Y mettre le code suivant :
-------------------------------------------------------------------------------------------------------------
%initialisation des variables
y= [3.2; 4.8;3.3; 3.2; 3.1; 4.2; 3.2; 3.3];
n= size(y);
z = zeros(n);
% la racine carré de la variance est égale à l'écart type
sigma = std(y);
% la moyenne de y
ym = mean(y) ;
% les valeurs de l'intervale
inf = ym-sigma;
sup = ym + sigma;
for i=1:n
if (inf <= y(i)) && (y(i) <= sup)
z(i) = 1;
end
end
-------------------------------------------------------------------------------------------------------------
En exécutant logic dans la console Matlab et en affichant z, on obtient :
>> logic
>> z
z=
1
0
1
1
1
0
1
1
FEKAK Fatima-Ezzahra, ENSA, Fès
CP1, 2019-2020
