Introduction Matlab

Ce poly constitue une introduction l'utilisation de la plate-forme Matlab.

Prrequis:
Savoir naviguer sur Internet, dans les dossiers grs par le systme d'exploitation Windows. Les tudiants
devront se munir d'une disquette personnelle.

Objectifs:
- acqurir un certain "vocabulaire"
- acqurir une certaine "grammaire"
- acqurir une "mthodologie" dans la conception des programmes

Contrle des Connaissances:

Un ou des exercices similaires ceux noncs dans la suite sera tir au sort et propos(s) chaque tudiant
dans le cadre des examen de sessions I et II.

note: ce poly a t crit en utilisant la version 7 de MATLAB.

***
SOMMAIRE
1. GENERALITES ET PRISE EN MAIN.................................................................................................. 4
1.1. DEMARRAGE, QUITTER ........................................................................................................................ 4
1.2. AIDE, DOCUMENTATION EN LIGNE, CONFIGURATION DU BUREAU ........................................................ 4
1.3. CALCULS ELEMENTAIRES, ANS, FORMAT ............................................................................................. 4
1.4. HISTORIQUE DES COMMANDES............................................................................................................. 4
1.5. FICHIERS DE COMMANDES (SCRIPTS) ................................................................................................... 5
1.5.1. Principe gnral .......................................................................................................................... 5
1.5.2. O doit se trouver mon fichier de commande ?.......................................................................... 5
1.5.3. Commentaires et auto-documentation ........................................................................................ 5
1.5.4. Suppression de l'affichage .......................................................................................................... 5
1.5.5. Pause dans l'excution ................................................................................................................ 5
1. Premier programme sous Matlab........................................................................................... 5
2. Input et format ........................................................................................................................ 6
3. Addpath (ajoute un chemin).................................................................................................... 6
2. VARIABLES ET FONCTIONS PREDEFINIES .................................................................................. 6
2.1. VARIABLES .......................................................................................................................................... 6
2.2. EFFACEMENT DES VARIABLES .............................................................................................................. 7
2.3. VARIABLES PREDEFINIES ..................................................................................................................... 7
2.4. FONCTIONS PREDEFINIES ..................................................................................................................... 7
4. Angle et abs............................................................................................................................. 7
5. Type......................................................................................................................................... 7
3. MATRICES ET TABLEAUX ................................................................................................................. 7
3.1. DEFINITION D'UN TABLEAU .................................................................................................................. 7
3.2. ACCES A UN ELEMENT D'UN TABLEAU.................................................................................................. 8
3.3. EXTRACTION DE SOUS-TABLEAUX ....................................................................................................... 9
3.4. CONSTRUCTION DE TABLEAUX PAR BLOCS........................................................................................... 9
6. Matrice et extraction............................................................................................................. 10
3.5. OPERATIONS SUR LES TABLEAUX ....................................................................................................... 10
1
 3.5.1. Addition et soustraction............................................................................................................ 10
Matrice d'lments alatoires............................................................................................... 10
7.
3.5.2. Multiplication, division et puissance terme terme ................................................................. 10
3.5.3. Multiplication, division et puissance au sens matriciel ............................................................ 10
8. Oprations sur une matrice .................................................................................................. 11
3.5.4. Transposition ............................................................................................................................ 11
3.6. LONGUEURS DE TABLEAU .................................................................................................................. 11
3.7. GENERATION RAPIDE DE TABLEAUX .................................................................................................. 11
3.7.1. Matrices classiques ................................................................................................................... 11
3.7.2. Listes de valeurs ....................................................................................................................... 11
3.8. MATRICES MULTIDIMENSIONNELLES ................................................................................................. 12
9. Rsolution d'un systme d'quations..................................................................................... 13
10. Gnration rapide de matrices ......................................................................................... 13
11. Gnration rapide de vecteurs.......................................................................................... 13
12. Manipulation de matrices ................................................................................................. 13
4. GRAPHIQUE 2D.................................................................................................................................... 14
4.1. L'INSTRUCTION PLOT ......................................................................................................................... 14
4.1.1. Tracer une courbe simple.......................................................................................................... 14
4.1.2. Superposer plusieurs courbes ................................................................................................... 15
4.1.3. Le pige classique..................................................................................................................... 15
4.1.4. Attributs de courbes.................................................................................................................. 15
4.2. DECORATION DES GRAPHIQUES.......................................................................................................... 16
4.2.1. Titre........................................................................................................................................... 16
4.2.2. Labels........................................................................................................................................ 16
4.2.3. Lgendes ................................................................................................................................... 16
4.2.4. Tracer un quadrillage................................................................................................................ 17
4.3. AFFICHER PLUSIEURS GRAPHIQUES (SUBPLOT) .................................................................................. 17
4.4. AXES ET ZOOM ................................................................................................................................... 17
4.5. INSTRUCTIONS GRAPHIQUES DIVERSES .............................................................................................. 18
4.5.1. Maintien du graphique.............................................................................................................. 18
4.5.2. Effacement de la fentre graphique .......................................................................................... 18
4.5.3. Saisie d'un point la souris....................................................................................................... 18
13. Subplot et ginput ............................................................................................................... 18
5. PROGRAMMATION MATLAB.......................................................................................................... 19
5.1. FONCTIONS ........................................................................................................................................ 19
5.1.1. Fonctions inline ........................................................................................................................ 19
5.1.2. Fonctions dfinies dans un fichier ............................................................................................ 20
14. Crer une fonction: vec2col.............................................................................................. 21
15. Equations du 1ier et du second degr ................................................................................ 21
16. Nargin, nargout ................................................................................................................ 21
5.1.3. Porte des variables .................................................................................................................. 21
17. Matrice de Vandemonde ................................................................................................... 22
18. Variable locale.................................................................................................................. 22
19. Variable globale ............................................................................................................... 22
5.2. STRUCTURES DE CONTROLE ............................................................................................................... 23
5.2.1. Oprateurs de comparaison et logiques .................................................................................... 23
20. Utilisation d'oprateurs logiques ..................................................................................... 24
5.2.2. La commande find .................................................................................................................... 24

2
 21. Reshape et find.................................................................................................................. 25
5.2.3. Instructions conditionnelles if................................................................................................... 25
22. Fonction avec if et is*....................................................................................................... 25
5.2.4. Boucles for................................................................................................................................ 26
23. Boucle conditionnelle, break ............................................................................................ 26
24. Tic et toc ........................................................................................................................... 26
5.2.5. Boucles while............................................................................................................................ 26
25. Boucle tant que ................................................................................................................. 26
5.2.6. return......................................................................................................................................... 26
5.2.7. Switch, case, otherwise............................................................................................................. 26
26. Switch................................................................................................................................ 26
27. Boucles implicites ............................................................................................................. 27
5.3. TYPES DE DONNEES (DATA TYPES)..................................................................................................... 27
5.3.1. Types numriques......................................................................................................................... 27
5.3.2. Type logique ............................................................................................................................. 27
5.3.3. Caractres et chanes................................................................................................................. 28
5.3.4. Structures .................................................................................................................................. 28
5.3.5. Cellules ..................................................................................................................................... 28
28. Fonctions signenombre et signematrice ........................................................................... 29
29. Ecrire une fonction dont l'en-tte est : ............................................................................. 29
30. Fonction estentier ............................................................................................................. 29
31. Types numriques ............................................................................................................. 29
32. Type logique...................................................................................................................... 29
33. Type caractres et chanes de caractres......................................................................... 29
34. Types structures et cellules............................................................................................... 30
6. CONSEILS DANS LA PROGRAMMATION..................................................................................... 30
6.1. CONVENTIONS D'ECRITURE ................................................................................................................ 30
6.2. FICHIERS, ORGANISATION, AMELIORATION DES PERFORMANCES ....................................................... 31
7. GRAPHIQUES 3D.................................................................................................................................. 32
7.1. COURBES EN 3D ................................................................................................................................ 32
7.2. SURFACES .......................................................................................................................................... 32
7.2.1. Gnration des points (meshgrid) ............................................................................................. 32
7.2.2. Trac de la surface .................................................................................................................... 34
7.2.3. Courbes de contour ................................................................................................................... 34
7.2.4. Contrle de l'angle de vue......................................................................................................... 34
35. Trac d'une surface en 3D ................................................................................................ 34
8. ECHANGES ENTRE MATLAB ET L'EXTERIEUR........................................................................ 35
8.1. SAUVEGARDE DE DONNEES ................................................................................................................ 35
8.2. IMPORTER DES TABLEAUX ................................................................................................................. 35
8.3. INCLURE DES COURBES MATLAB DANS UN DOCUMENT ................................................................... 35
9. CALCUL NUMERIQUE AVEC MATLAB ........................................................................................ 35
9.1. RECHERCHE DES ZEROS D'UNE FONCTION .......................................................................................... 36
9.2. INTERPOLATION ................................................................................................................................. 37
9.3. APPROXIMATION (ESTIMATION DE PARAMETRES OU "FITTING" OU "AJUSTEMENT") .......................... 37
9.3.1. Linaire ..................................................................................................................................... 38
36. Simulation d'une exponentielle dcroissante bruite ....................................................... 39

3
MATLAB (MATrix LABoratory) comprend de nombreuses fonctions graphiques, un systme puissant
d'oprateurs s'appliquant des matrices, des algorithmes numriques, ainsi qu'un langage de programmation
extrmement simple utiliser. L'aspect modulaire est l'un des plus grands atouts de MATLAB : l'utilisateur
peut lui-mme dfinir ses propres fonctions, en regroupant des instructions MATLAB dans un fichier portant
le suffixe ".m". La syntaxe est bien plus abordable que dans les langages classiques et devrait liminer les
rticences habituelles des programmeurs dbutants pour crire des fonctions. En termes de vitesse
d'excution, les performances sont infrieures celles obtenues avec un langage de programmation classique.
L'emploi de MATLAB devrait donc tre restreinte des problmes peu gourmands en temps calcul, mais
dans la plupart des cas, il prsente une solution lgante et rapide mettre en oeuvre.

1. Gnralits et prise en main

1.1. Dmarrage, quitter
Pour lancer le programme, double-cliquez l'icne matlab sur le bureau. Une fentre logo fait une brve
apparition, puis dans la fentre de commande, le symbole ">>" apparat : c'est l'invite de MATLAB qui
attend vos commandes. Vous pourrez quitter la session avec la commande quit.

1.2. Aide, documentation en ligne, configuration du bureau

L'aide en ligne peut tre obtenue directement dans la session en tapant help nom de commande.

1.3. Calculs lmentaires, ans, format

Commenons par les oprateurs les plus courants : +, -, *, /, ^. Le dernier signifie "puissance". Les
parenthses s'utilisent de manire classique. Nous avons tout pour effectuer un premier calcul : tapez une
expression mathmatique quelconque et appuyez sur <Entre>. Par exemple :
>> (3*2)/(5+3)
ans =
0.7500
Le rsultat est mis automatiquement dans une variable appele ans (answer). Celle-ci peut tre utilise pour
le calcul suivant, par exemple :
>> ans*2
ans =
1.5000
Ensuite, vous remarquerez que le rsultat est affich avec 5 chiffres significatifs, ce qui ne signifie pas que
les calculs soient faits avec aussi peu de prcision. Si vous voulez afficher les nombres avec plus de
prcision, tapez la commande format long. Pour revenir au comportement initial : format short.

1.4. Historique des commandes

Toutes les commandes que vous aurez tap sous MATLAB peuvent tre retrouves et dites grce aux
touches de direction. Appuyez sur pour remonter dans les commandes prcdentes, pour redescendre et
utilisez et pour diter une commande. Pour relancer une commande, inutile de remettre le curseur la
fin, vous appuyez directement sur la touche <Entre>. Encore plus fort : vous pouvez retrouver toutes les
commandes commenant par un groupe de lettres. Par exemple pour retrouver toutes les commandes
commenant par <plot>, tapez plot, puis appuyez plusieurs fois sur . Exprimenter cela ds maintenant.

4
 1.5. Fichiers de commandes (scripts)
1.5.1. Principe gnral
Le principe est simple : regrouper dans un fichier une srie de commandes MATLAB et les excuter en bloc.
Tout se passera comme si vous les tapiez au fur et mesure dans une session MATLAB. Il est fortement
conseill de procder comme ceci. Cela permet notamment de rcuprer facilement votre travail de la veille.
Les fichiers de commandes peuvent porter un nom quelconque, mais doivent finir par l'extension ".m".
Pour crer un fichier de commandes, taper des commandes MATLAB l'intrieur de l'diteur intgr, sauver
par exemple sous le nom "toto.m", puis sous MATLAB, tapez :
>> toto
Toutes les commandes du fichier seront excutes en bloc.

1.5.2. O doit se trouver mon fichier de commande ?

Crer votre rpertoire personnel sur le disque dur avec le navigateur Windows. Utiliser votre nom pour le
nom du rpertoire. Pour que Matlab utilise les scripts qui y seront placs, modifier le rpertoire courant de
Matlab pour l'assigner celui que vous venez de crer : barre des menus "current repertory". Tous les scripts
qui seront crits dans la suite devront tre sauvs dans ce rpertoire ou ventuellement sur une disquette avec
un nom votre convenance et suffisamment explicite.

1.5.3. Commentaires et auto-documentation

Les lignes blanches sont considres comme des commentaires. Tout ce qui se trouve aprs le symbole "%"
sera galement considr comme un commentaire.
Il est galement possible de documenter ses fichiers de commande. Ainsi, dfinissez une srie de lignes de
commentaires ininterrompues au dbut de votre fichier toto.m. Lorsque vous taperez :
>> help toto
ces lignes de commentaires apparatront l'cran. C'est ainsi que fonctionne l'aide en ligne de MATLAB.
C'est tellement simple que ce serait dommage de s'en priver.

1.5.4. Suppression de l'affichage

Jusqu' prsent, nous avons vu que toutes les commandes tapes sous MATLAB affichaient le rsultat. Pour
certaines commandes (cration de gros tableaux), cela peut s'avrer fastidieux. On peut donc placer le
caractre ";" la fin d'une ligne de commande pour indiquer MATLAB qu'il ne doit pas afficher le rsultat.

1.5.5. Pause dans l'excution

Si vous entrez la commande pause dans un fichier de commandes, le programme s'arrtera cette ligne tant
que vous n'avez pas tap < Entre>.
1. Premier programme sous Matlab
Ecrire l'exemple suivant, le sauver sous le nom 'premierprog' et l'excuter partir de la fentre de
commande.
% premier programme sous Matlab
disp('Ceci est mon premier programme')
disp(' Il affiche l''cran les 10 premiers entiers')
for i=1:10
i
end
% fin du programme
Modifier le programme en ajoutant un ';' aprs le 'i'. Que donne maintenant son execution?

5
2. Input et format
Ecrire un script qui invite l'utilisateur entrer un nombre, l'affiche en format court simple puis court avec
puissance de 10. Indices : input, format. La prcision dans un calcul avec chacun de ces nombres est-elle
affecte?
3. Addpath (ajoute un chemin)
Ecrire un programme simple et le sauver dans un rpertoire qui ne soit pas le rpertoire courant. Tenter de
l'excuter en l'appelant de la fentre de commandes. Utiliser la commande addpath pour rendre l'appel
effectif.

2. Variables et fonctions prdfinies

2.1. Variables
Le calcul effectu plus haut n'a gure d'intrt en soi. Il est bien sr possible de conserver un rsultat de
calcul et de le stocker dans des variables. Gros avantage sur les langages classiques : on ne dclare pas les
variables. Leur type (entier, rel, complexe) s'affectera automatiquement en fonction du calcul effectu. Pour
affecter une variable, on dit simplement quoi elle est gale. Exemple :
>> a=1.2
a=
1.2000
On peut maintenant inclure cette variable dans de nouvelles expressions mathmatiques, pour en dfinir une
nouvelle :
>> b = 5*a^2+a
b=
8.4000
et ensuite utiliser ces deux variables :
>> c = a^2 + b^3/2
c=
297.7920
J'ai maintenant trois variables a, b et c. Ces variables ne sont pas affiches en permanence l'cran. Mais
pour voir le contenu d'une variable, rien de plus simple, on tape son nom :
>> b
b=
8.4000
ou on double-click sur son nom dans le workspace.
On peut aussi faire des calculs en complexe 1 s'crit indiffremment i ou j, donc pour dfinir un
complexe :
>> a+ b*i
ans =
1.2000 + 8.4000i
Le symbole * peut tre omis si la partie imaginaire est une constante numrique. Tous les oprateurs
prcdents fonctionnent en complexe. Par exemple :
>> (a+b*i)^2
ans =
-69.1200 +20.1600i
Un dernier point sur les variables :
- MATLAB fait la diffrence entre les minuscules et les majuscules.
- Les noms de variables peuvent avoir une longueur quelconque.
- Les noms de variables doivent commencer par une lettre.
6
 2.2. Effacement des variables
La commande clear permet d'effacer une partie ou toutes les variables dfinies jusqu' prsent.
Syntaxe :
clear a b c . . .
Si aucune variable n'est spcifie, toutes les variables seront effaces.

2.3. Variables prdfinies

Il existe un certain nombre de variables pr-existantes. Nous avons dj vu ans qui contient le dernier rsultat
de calcul, ainsi que i et j qui reprsentent 1 .
Il existe aussi pi, qui reprsente , et quelques autres. Retenez que eps, nom que l'on a souvent tendance
utiliser est une variable prdfinie. ATTENTION : ces variables ne sont pas protges, donc si vous les
affectez, elles ne gardent pas leur valeur initiale. C'est souvent le problme pour i et j que l'on utilise souvent
spontanment comme indices de boucles, de telle sorte qu'on ne peut plus ensuite dfinir de complexe.
Quelle est la valeur de eps ?

2.4. Fonctions prdfinies

Toutes les fonctions courantes et beaucoup parmi les moins courantes existent. La plupart d'entre elles
fonctionnent en complexe. On retiendra que pour appliquer une fonction une valeur, il faut mettre cette
dernire entre parenthses. Exemple :
>> sin(pi/12)
ans =
0.16589613269342
Voici une liste non exhaustive :
- fonctions trigonomtriques et inverses : sin, cos, tan, asin, acos, atan
- fonctions hyperboliques (on rajoute "h") : sinh, cosh, tanh, asinh, acosh, atanh
- racine, logarithmes et exponentielles : sqrt, log, log10, exp
- fonctions erreur : erf, erfc
- fonctions de Bessel et Hankel : besselj, bessely, besseli, besselk, besselh.
4. Angle et abs
Chercher la syntaxe utilise pour les fonctions angle et abs. Dterminer alors, la phase et le module du
nombre complexe 1+4*j.
5. Type
Essayer type premierprog, puis type sin et encore type angle. Que fait la commande type?

3. Matrices et tableaux
En dpit du titre de cette section, MATLAB ne fait pas de diffrence entre les deux. Le concept de tableau est
important car il est la base du graphique : typiquement pour une courbe de n points, on dfinira un tableau
de n abscisses et un tableau de n ordonnes. Mais on peut aussi dfinir des tableaux rectangulaires deux
indices pour dfinir des matrices au sens mathmatique du terme, puis effectuer des oprations sur ces
matrices.

3.1. Dfinition d'un tableau

On utilise les crochets [ et ] pour dfinir le dbut et la fin de la matrice. Ainsi pour dfinir une variable M
1 2 3
11 12 13
contenant la matrice 21 22 23 , on crira :
7
>> M = [1,2,3;11,12,13;21,22,23]
M=
1 2 3
11 12 13
21 22 23
On utilise le symbole "," qui sert de sparateur de colonne et ";" de sparateur de ligne. On peut aussi dfinir
des vecteurs, ligne ou colonne, partir de cette syntaxe (un espace est interprt comme une virgule). Par
exemple :
>> U = [1 2 3]
U=
123
dfinit un vecteur ligne, alors que (une commande "entre" est interprte comme un poit-virgule):
>> V = [11
12
13]
V=
11
12
13
dfinit un vecteur colonne. Le passage d'une ligne de commande la suivante s'effectuant par la frappe de la
touche <Entre>.
On aurait pu aussi dfinir ce dernier par :
>> V=[11;12;13]

3.2. Accs un lment d'un tableau

Il suffit d'entrer le nom du tableau suivi entre parenthses du ou des indices dont on veut lire ou crire la
valeur. Exemple si je veux la valeur de M, ligne 3, colonne 2 :
>> M(3,2)
ans =
22
Pour modifier seulement un lment d'un tableau, on utilise le mme principe. Par exemple, je veux que M32
soit gal 32 au lieu de 22 :
>> M(3,2)=32
M=
1 2 3
11 12 13
21 32 23
Vous remarquerez que MATLAB raffiche du coup toute la matrice, en prenant en compte la modification.
On peut se demander ce qui se passe si on affecte la composante d'une matrice qui n'existe pas encore.
Exemple :
>> P(2,3) = 3
P=
000
003
Voil la rponse : MATLAB construit automatiquement un tableau suffisamment grand pour arriver
jusqu'aux indices spcifis, et met des zros partout sauf au terme considr. Vous remarquerez que
contrairement aux langages classiques, inutile de dimensionner les tableaux l'avance : ils se construisent au
fur et mesure.

8
 3.3. Extraction de sous-tableaux
Il est souvent utile d'extraire des blocs d'un tableau existant. Pour cela on utilise le caractre ":". Il faut
spcifier pour chaque indice la valeur de dbut et la valeur de fin. La syntaxe gnrale est donc la suivante
(pour un tableau deux indices) :
tableau(dbut : fin, dbut : fin)
2 3

Ainsi pour extraire le bloc 12 13
on tapera :
>> M(1:2,2:3)
ans =
2 3
12 13
On utilise le caractre ":" seul pour prendre tous les indices possibles. Exemple :
>> M(1:2,:)
ans =
1 2 3
11 12 13
C'est bien pratique pour extraire des lignes ou des colonnes d'une matrice. Par exemple pour obtenir la
deuxime ligne de M :
>> M(2,:)
ans =
11 12 13

3.4. Construction de tableaux par blocs

Vous connaissez ce principe en mathmatiques. Par exemple, partir des matrices et vecteurs prcdemment
dfinis, on veut dfinir la matrice
M V
U 0

qui est une matrice 4x4. Pour faire cela sous MATLAB, on fait comme si les blocs taient des scalaires, et on
crit tout simplement :
>> N=[M V
U 0]
N=
1 2 3 11
11 12 13 12
21 32 23 13
1 2 3 0

ou bien en utilisant les caractres "," et ";"

>> N=[M, V; U, 0]

Cette syntaxe est trs utilise pour allonger des vecteurs ou des matrices, par exemple si je veux ajouter une
colonne M, constitue par V :
>> M=[M V]
M=
1 2 3 11
1 12 13 12
21 32 23 13
9
Si je veux lui ajouter une ligne, constitue par U :
>> M = [M;U]
M=
1 2 3
11 12 13
21 32 23
1 2 3
6. Matrice et extraction
Ecrire un script qui dfinisse la matrice M suivante et affiche la dernire ligne:

1 2 3 4 5
9 10 11 12 13
-1 -2 -3 -4 -5

3.5. Oprations sur les tableaux

3.5.1. Addition et soustraction
Les deux oprateurs sont les mmes que pour les scalaires. A partir du moment o les deux tableaux
concerns ont la mme taille, le tableau rsultant est obtenu en ajoutant ou soustrayant les termes de chaque
tableau.
7. Matrice d'lments alatoires
Dterminer la syntaxe de la fonction rand, crer 2 matrices A et B de taille 5X6, valeurs alatoires
comprises entre -10 et 10. Afficher C= A+B

3.5.2. Multiplication, division et puissance terme terme

Ces oprateurs sont nots ".*", "./" et ".^" (attention ne pas oublier le point). Ils sont prvus pour effectuer
des oprations termes terme sur deux tableaux de mme taille. Ces symboles sont importants lorsque l'on
veut tracer des courbes.

3.5.3. Multiplication, division et puissance au sens matriciel

Puisque l'on veut manipuler des matrices, il parat intressant de disposer d'un multiplication matricielle.
Celle-ci se note simplement * et ne doit pas tre confondue avec la multiplication terme terme. Il va de soi
que si l'on crit A*B le nombre de colonnes de A doit tre gal au nombre de lignes de B pour que la
multiplication fonctionne.
La division a un sens vis- -vis des inverses de matrices. Ainsi A/B reprsente A multipli (au sens des
matrices) par la matrice inverse de B. Il existe aussi une division gauche qui se note \. Ainsi A\B signifie
l'inverse de A multipli par B. Ce symbole peut tre aussi utilis pour rsoudre des systmes linaires : si v
est un vecteur, A\v reprsente mathmatiquement A-1v, c'est- -dire la solution du systme linaire Ax = v.
La puissance n-ime d'une matrice reprsente cette matrice multiplie n fois au sens des matrices par elle-
mme.
Pour bien montrer la diffrence entre les oprateurs ".*" et "*", un petit exemple faisant intervenir la matrice
1 2

identit multiplie par la matrice 3 4 .
Voici la multiplication au sens des matrices :
>> [1 0; 0 1] * [1 2; 3 4]
ans =
12
34

10
et maintenant la multiplication terme terme :
>> [1 0; 0 1] .* [1 2; 3 4]
ans =
10
04

8. Oprations sur une matrice

a=[1,4;6,8]; Est-il possible d'appliquer la fonction sin la matrice a ? Quel est le rsultat ? Ecrire un script
qui affiche les rsultats successifs de sin(a).^n avec n variant de 1 10 et utilisant une boucle for...end.

3.5.4. Transposition
L'oprateur transposition est le caractre ' (prime) et est souvent utilis pour transformer des vecteurs lignes
en vecteurs colonnes et inversement.

3.6. Longueurs de tableau

La fonction size applique une matrice renvoie un tableau de deux entiers : le premier est le nombre de
lignes, le second le nombre de colonnes. La commande fonctionne aussi sur les vecteurs et renvoie 1 pour le
nombre de lignes (resp. colonnes) d'un vecteur ligne (resp colonne).
Pour les vecteurs, la commande length est plus pratique et renvoie le nombre de composantes du vecteur,
qu'il soit ligne ou colonne.

3.7. Gnration rapide de tableaux

3.7.1. Matrices classiques
On peut dfinir des matrices de taille donne ne contenant que des 0 avec la fonction zeros, ou ne contenant
que des 1 avec la fonction ones. Il faut spcifier le nombre de lignes et le nombre de colonnes. Voici deux
exemples :
>> ones(2,3)
ans =
111
111
>> zeros(1,3)
ans =
000
L'identit est obtenue avec eye. On spcifie seulement la dimension de la matrice (qui est carre. . .)
>> eye(3)
ans =
100
010
001
Il existe galement une fonction diag permettant de crer des matrices diagonales (voir l'aide intgre).

3.7.2. Listes de valeurs

Cette notion est capitale pour la construction de courbes. Il s'agit de gnrer dans un vecteur une liste de
valeurs quidistantes entre deux valeurs extrmes. La syntaxe gnrale est :
variable = valeur dbut: pas: valeur fin
Cette syntaxe cre toujours un vecteur ligne. Par exemple pour crer un vecteur x de valeurs quidistantes de
0.1 entre 0 et 1 :
>> x = 0:0.1:1
11
x=
Columns 1 through 7
0 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000
Columns 8 through 11
0.7000 0.8000 0.9000 1.0000
Il est conseill de mettre un point-virgule la fin de ce type d'instruction pour viter l'affichage fastidieux du
rsultat.
Autre exemple pour crer 101 valeurs qui-rparties sur l'intervalle [0; 2] :
>> x = 0: 2*pi/100 : 2*pi;
On peut aussi utiliser la fonction linspace :
>> x=linspace(0,2*pi,101);
Les valeurs peuvent tre rparties de manire logarithmique avec la fonction logspace.

3.8. Matrices multidimensionnelles

Ci-dessous, tableau tri-dimensionnel, pour accder l'lment situ sur la seconde ligne, troisime colonne,
quatrime page, utiliser les indices (2,3,4)

autre exemple, la matrice A :

12
9. Rsolution d'un systme d'quations
Soit le systme de 5 quations 5 inconnues :
1=x+2y+3z+4u+5v
2=x+2y+3z+4u+v
3=x+2y+3z+u+v
4=-x-2y+3z+4u+v
5=x+y+z+u+v
Ecrire le script rsolvant ce systme
10. Gnration rapide de matrices
Soit a=ones(2,3). Est-il possible de calculer a-1, a+a, 2*a, a*a, a.*a ? Analyser les rsultats.
11. Gnration rapide de vecteurs
En une ligne de commande, crer le vecteur colonne ayant comme valeurs (sans les crire explicitement) :
[10-60;10-59; ...; 1059; 1060]
indice : logspace
12. Manipulation de matrices
Crer la matrice suivante en utilisant la commande reshape :

1 11 21 31 41 51 61 71 81 91
2 12 22 32 42 52 62 72 82 92
3 13 23 33 43 53 63 73 83 93
4 14 24 34 44 54 64 74 84 94
5 15 25 35 45 55 65 75 85 95
6 16 26 36 46 56 66 76 86 96
7 17 27 37 47 57 67 77 87 97
8 18 28 38 48 58 68 78 88 98
9 19 29 39 49 59 69 79 89 99
10 20 30 40 50 60 70 80 90 100

Utiliser triu pour crer :

1 11 21 31 41 51 61 71 81 91
0 12 22 32 42 52 62 72 82 92
0 0 23 33 43 53 63 73 83 93
0 0 0 34 44 54 64 74 84 94
0 0 0 0 45 55 65 75 85 95
0 0 0 0 0 56 66 76 86 96
0 0 0 0 0 0 67 77 87 97
0 0 0 0 0 0 0 78 88 98
0 0 0 0 0 0 0 0 89 99
0 0 0 0 0 0 0 0 0 100

13
 4. Graphique 2D
Une courbe 2D est reprsent par une srie d'abscisses et une srie d'ordonnes. Le logiciel trace
gnralement des droites entre ces points. La fonction de MATLAB s'appelle plot.

4.1. L'instruction plot

4.1.1. Tracer une courbe simple
L'utilisation la plus simple de l'instruction plot est la suivante.
plot ( vecteur d'abscisses, vecteur d'ordonnes ) ou encore,
plot ([ x1 x2 ... xn ], [ y1 y2 ... yn ])
Les vecteurs peuvent tre indiffremment ligne ou colonne, pourvu qu'ils soient tous deux de mme type. En
gnral ils sont lignes car la gnration de listes de valeurs vue prcdemment fournit par dfaut des vecteurs
lignes.
Par exemple, si on veut tracer sin(x) sur l'intervalle [0; 2], on commence par dfinir une srie (de taille
raisonnable, disons 100) de valeurs quidistantes sur cet intervalle :
>> x = 0: 2*pi/100 : 2*pi;
puis, comme la fonction sin peut s'appliquer terme
terme un tableau :
>> plot(x, sin(x))
qui fournit le graphe suivant dans la fentre graphique :

On voit que les axes s'adaptent automatiquement aux

valeurs extrmales des abscisses et ordonnes.
On remarquera que tout ce que demande plot, c'est un
vecteur d'abscisses et un vecteur d'ordonnes. Les
abscisses peuvent donc tre une fonction de x plutt que x
lui-mme. En d'autres termes, il est donc possible de tracer
des courbes paramtres :
>> plot(cos(x), sin(x))

14
 4.1.2. Superposer plusieurs courbes
Il suffit de spcifier autant de couples (abscisses, ordonnes) qu'il y a de courbes tracer. Par exemple pour
superposer sin et cos :
>> plot(x, cos(x),x, sin(x))

Les deux courbes tant en ralit dans des couleurs diffrentes. Cette mthode fonctionne mme si les
abscisses des deux courbes ne sont pas les mmes.

4.1.3. Le pige classique

Vous tomberez rapidement dessus, et rassurez-vous, vous ne serez pas le premier. Essayez par exemple de
tracer le graphe de la fonction x ==> x*sin (x) sur [0; 2]. Fort des deux section prcdentes, vous y allez
franco :
>> x = 0:0.1:2*pi;
>> plot(x, x*sin(x))
??? Error using ==> *
Inner matrix dimensions must agree.
Pour MATLAB vous manipulez des tableaux : x en est un, et sin(x) en est un deuxime de mme taille, et
tous deux sont des vecteurs lignes. Est-ce que cela a un sens de multiplier deux vecteurs lignes au sens des
matrices ? Non. C'est pourtant ce que vous avez fait, puisque vous avez utilis le symbole * . Vous avez
compris : il fallait utiliser la multiplication terme terme .* soit :
>> plot(x, x.*sin(x))

4.1.4. Attributs de courbes

MATLAB attribue des couleurs par dfaut aux courbes. Il est possible de
Specifier Color
modifier la couleur, le style du trait et celui des points, en spcifiant aprs
chaque couple (abscisse, ordonne) une chane de caractres (entre primes) r red
pouvant contenir les codes suivants (obtenus en tapant help plot) : g green
b blue
Specifier Line Style
c cyan
- solid line (default)
m magenta
-- dashed line
y yellow
: dotted line
k black
-. dash-dot line
w white

15
 Specifier Marker Type
+ plus sign
o circle
* asterisk
. point
x cross
s square
d diamond
^ upward pointing triangle
v downward pointing triangle
> right pointing triangle
< left pointing triangle
p five-pointed star (pentagram)
h six-pointed star (hexagram)

Lorsque l'on utilise seulement un style de points, MATLAB ne trace plus de droites entre les points
successifs, mais seulement les points eux-mme. Ceci peut tre pratique par exemple pour prsenter des
rsultats exprimentaux.
Les codes peuvent tre combins entre eux. Par exemple :
>> plot(x,sin(x),':',x,cos(x),'r-.')

4.2. Dcoration des graphiques

4.2.1. Titre
C'est l'instruction title laquelle il faut fournir une chane de caractres. Le titre apparat en haut de la fentre
graphique :
>> plot(x,cos(x),x,sin(x))
>> title('Fonctions sin et cos')

4.2.2. Labels
Il s'agit d'afficher quelque chose sous les abscisses et cot de l'axe des ordonnes :
>> plot(x,cos(x))
>> xlabel('Abscisse')
>> ylabel('Ordonne')

4.2.3. Lgendes
Utilisation de l'instruction legend. Il faut lui communiquer autant de chanes de caractres que de courbes
traces l'cran. Un cadre est alors trac sur le graphique, qui affiche en face du style de chaque courbe, le
texte correspondant. Par exemple :
>> plot(x,cos(x),':',x,sin(x),'-.',x,sqrt(x),'--')
>> legend('cosinus','sinus','racine')
16
 4.2.4. Tracer un quadrillage
C'est l'instruction grid, qui utilis aprs une instruction plot affiche un quadrillage sur la courbe. Si on tape
nouveau grid, le quadrillage disparat.

4.3. Afficher plusieurs graphiques (subplot)

Voil une fonctionnalit trs utile pour prsenter sur une mme page graphique plusieurs rsultats. L'ide
gnrale est de dcouper la fentre graphique en pavs de mme taille, et d'afficher un graphe dans chaque
pav. On utilise l'instruction subplot en lui spcifiant le nombre de pavs sur la hauteur, le nombre de pavs
sur la largeur, et le numro du pav dans lequel on va tracer :
subplot (Nbre pavs sur hauteur, Nbre pavs sur largeur, Numro pav)
La virgule peut tre omise. Les pavs sont numrots dans le sens de la lecture d'un texte :
de gauche droite et de haut en bas :
Une fois que l'on a tap une commande subplot, toutes les commandes graphiques suivantes seront excutes
dans le pav spcifi. Ainsi, le graphique suivant a t obtenu partir de la suite d'instructions :
>> subplot(221)
>> plot(x,sin(x))
>> subplot(222)
>> plot(x,cos(x),x,sin(x),'-.')
>> subplot(223)
>> plot(cos(x),sin(x))
>> subplot(224)
>> plot(sin(2*x),sin(3*x))

4.4. Axes et zoom

Il y a deux manires de modifier les valeurs extrmes sur les axes, autrement dit de faire du zooming sur les
courbes. La plus simple est d'utiliser le bouton de la fentre graphique. Vous pouvez alors :
- encadrer une zone zoomer avec le bouton de gauche de la souris, ou
- cliquer sur un point avec le bouton de gauche. Le point cliqu sera le centre du zoom, ce dernier tant
effectu avec un facteur arbitraire
- cliquer sur un point avec le bouton de droite pour dzoomer
Pour faire des zooms plus prcis, utiliser la commande axis (voir l'aide intgre).

17
 4.5. Instructions graphiques diverses
4.5.1. Maintien du graphique
Par dfaut une instruction plot efface systmatiquement le graphique prcdent. Il est parfois utile de le
conserver et de venir le surcharger avec une nouvelle courbe. Pour cela on utilise la commande hold. Pour
dmarrer le mode surcharge, taper hold on, pour revenir en mode normal, hold off.
Il est conseill de ne pas abuser de cette commande.

4.5.2. Effacement de la fentre graphique

Tapez simplement clf. Cette commande annule galement toutes les commandes subplot et hold passes.

4.5.3. Saisie d'un point la souris

La commande ginput(N) permet de cliquer N points dans la fentre graphique. La commande renvoie alors
deux vecteurs, l'un contenant les abscisses, l'autre les ordonnes. Utilise sans le paramtre N, la commande
tourne en boucle jusqu ce que la touche "Entre" soit tape.

13. Subplot et ginput

Ecrire un script qui reprsente les fonctions cos(5t) (ligne bleue et des 'X') et tan(t) (ligne rouge et des '+')
dans l'intervalle t [0,10] sur 2 graphes d'une mme figure, placer titres et labels. Utiliser la fonction text
pour inviter l'utilisateur cliquer 2 fois dans l'un des graphes. Le script comportera la fonction ginput pour
rcuprer le rsultat des 2 clics et tracera une ligne entre les 2 points ainsi dfinis. Obtenir quelque chose
qui ressemble la figure ci-dessous:

la fonction cosinus
1

0.5
cos(5t)

-0.5

-1
0 1 2 3 4 5 6 7 8 9 10
t (s)

la fonction tangente
2

1
tan(2.5t)

0 Cliquer 2 fois pour tracer une ligne

-1

-2
0 1 2 3 4 5 6 7 8 9 10
t (s)

18
 5. Programmation MATLAB
5.1. Fonctions
Nous avons vu un certain nombre de fonctions prdfinies. Il est possible de dfinir ses propres fonctions. La
premire mthode permet de dfinir des fonctions simples sur une ligne de commande. La seconde, beaucoup
plus gnrale permet de dfinir des fonctions trs volues en la dfinissant dans un fichier.

5.1.1. Fonctions inline

Admettons que je veuille dfinir une nouvelle fonction que j'appelle sincos dfinie mathmatiquement par :
sincos(x) = sin (x) - x cos (x)
On crira :
>> sincos = inline('sin(x)-x*cos(x)', 'x')
sincos =
Inline function:
sincos(x) = sin(x)-x*cos(x)
C'est trs simple. N'oubliez pas les primes, qui dfinissent dans MATLAB des chanes de caractres. On peut
maintenant utiliser cette nouvelle fonction :
>> sincos(pi/12)
ans =
0.0059
Essayons maintenant d'appliquer cette fonction un tableau de valeurs :
>> sincos(0:pi/3:pi)
??? Error using ==> inline/subsref
Error in inline expression ==> sin(x)-x*cos(x)
??? Error using ==> *
Inner matrix dimensions must agree.
Cela ne marche pas. Vous avez compris pourquoi ? (c'est le mme problme que pour l'instruction
plot(x,x*sin(x)) vue prcdemment) MATLAB essaye de multiplier x vecteur ligne par cos(x) aussi vecteur
ligne au sens de la multiplication de matrice ! Par consquent il faut bien utiliser une multiplication terme
terme .* dans la dfinition de la fonction :
>> sincos = inline('sin(x)-x.*cos(x)')
sincos =
Inline function:
sincos(x) = sin(x)-x.*cos(x)

>> sincos(0:pi/3:pi)
ans =
0 0.3424 1.9132 3.1416
On peut donc noncer la rgle suivante :
Lorsque l'on dfinit une fonction, il est prfrable d'utiliser systmatiquement les oprateurs terme terme .*,
./ et .^ au lieu de * / et ^ si l'on veut que cette fonction puisse s'appliquer des tableaux.
On peut de mme dfinir des fonctions de plusieurs variables :
>> sincos = inline('sin(x)-y.*cos(x)', 'x', 'y')
sincos =
Inline function:
sincos(x,y) = sin(x)-y.*cos(x)
L'ordre des variables (affich l'cran) est celui dans lequel elles apparaissent dans la dfinition de la
fonction. Essayer :
>> sincos(1,2)

19
 5.1.2. Fonctions dfinies dans un fichier
Dans les mathmatiques on crit souvent les fonctions sous la forme :
y = f(x) ou z = g(x, y) ou ....
Nous verrons que c'est la mme chose sous MATLAB. Les noms de fonctions doivent tre plus explicites
que f ou g. Par exemple une fonction renvoyant le carr d'un argument x pourra tre nomme carre (on vite
les accents):
y = carre(x)
Dans le fichier dfinissant la fonction, cette syntaxe est aussi prcde du mot cl function :
Forme gnrale d'une fonction dans un fichier nomfonction.m :

function Resultat = nomfonction(Parametres_Passes_a_la_Fonction)

% lignes de commentaires (aide descriptive apparaissant dans la fentre de commandes % si on tape :
help nomfonction )
lignes de commandes
...
Resultat=...

Commenons par reprendre l'exemple prcdent (sincos). L'ordre des oprations est le suivant :
1. Editer un nouveau fichier appel sincos.m
2. Taper les lignes suivantes :

function s = sincos(x)
% sincos calcule sin(x)-x.*cos(x)
s = sin(x)-x.*cos(x);

3. Sauvegarder
Le rsultat est le mme que prcdemment :
>> sincos(pi/12)
ans =
0.0059
On remarquera plusieurs choses :
- l'utilisation de .* pour que la fonction soit applicable des tableaux.
- la variable s n'est l que pour spcifier la sortie de la fonction.
- l'emploi de ; la fin de la ligne de calcul, sans quoi MATLAB afficherait deux fois le rsultat de la
fonction.
Donc encore une rgle :
Lorsque l'on dfinit une fonction dans un fichier, il est prfrable de mettre un ; la fin de chaque commande
constituant la fonction. Attention cependant ne pas en mettre sur la premire ligne.
Un autre point important : le nom du fichier doit porter l'extension .m et le nom du fichier sans suffixe doit
tre exactement le nom de la fonction (apparaissant aprs le mot-cl function)
En ce qui concerne l'endroit de l'arborescence o le fichier doit tre plac, la rgle est la mme que pour les
scripts.

Voyons maintenant comment dfinir une fonction comportant plusieurs sorties. On veut raliser une fonction
appele cart2pol qui convertit des coordonnes cartsiennes (x; y) (entres de la fonction) en coordonnes
polaires (r; fi) (sorties de la fonction). Voil le contenu du fichier cart2pol.m :

function [r, theta] = cart2pol (x, y)

% cart2pol convertit des coordonnes cartsiennes (x; y) (entres de la fonction)
% en coordonnes polaires (r; fi) (sorties de la fonction)
20
% sqrt = racine carre (fonction interne de Matlab)
% atan = fonction rciproque de tangente (fonction interne de Matlab)
r = sqrt(x.^2 + y.^2);
theta = atan (y./x);

On remarque que les deux variables de sortie sont mises entre crochets, et spares par une virgule.
Pour utiliser cette fonction, on crira par exemple :
>> [module,phase] = cart2pol(1,1)
module =
1.4142
phase =
0.7854
On affecte donc simultanment deux variables module et phase avec les deux sorties de la fonction, en
mettant ces deux variables dans des crochets, et spars par une virgule.
Il est possible de ne rcuprer que la premire sortie de la fonction. MATLAB utilise souvent ce principe
pour dfinir des fonctions ayant une sortie principale et des sorties optionnelles.
Ainsi, pour notre fonction, si une seule variable de sortie est spcifie, seule la valeur du rayon polaire est
renvoye. Si la fonction est appele sans variable de sortie, c'est ans qui prend la valeur du rayon polaire :
>> cart2pol(1,1)
ans =
1.4142
14. Crer une fonction: vec2col
Ecrire la fonction vec2col.m qui transforme tout vecteur (ligne ou colonne) pass en argument en vecteur
colonne. Un traitement d'erreur (fonction error) testera que la variable passe la fonction est bien un
vecteur et non une matrice (utiliser size et l'instruction conditionnelle if dcrite plus bas).
15. Equations du 1ier et du second degr
Complter les 2 fonctions dont les 2 premires lignes sont respectivement :

function x=equ1(a,b)
% resout a*x+b=0
et
function x=equ2(a,b,c)
% resout a*x^2+b*x+c=0

qui rsolvent les quations du premier et du second degr. Utiliser la fonction nargin pour que l'appel
equ2(a,b, c) soit valide. Vrifier que les solutions sont correctes.
16. Nargin, nargout
Chercher avec l'aide la syntaxe et l'utilisation de ces commandes. Les appliquer equa1 et equa2.

5.1.3. Porte des variables

A l'intrieur des fonctions comme celle que nous venons d'crire, vous avez le droit de manipuler trois types
de variables :
- Les variables d'entre de la fonction. Vous ne pouvez pas modifier leurs valeurs.
- Les variables de sortie de la fonction. Vous devez leur affecter une valeur.
- Les variables locales. Ce sont des variables temporaires pour dcouper des calculs par exemple. Elles n'ont
de sens qu' l'intrieur de la fonction et ne seront pas "vues" de l'extrieur. Et C'EST TOUT!

Imaginons maintenant que vous criviez un fichier de commande, et une fonction (dans deux fichiers
diffrents bien sr), le fichier de commande se servant de la fonction. Si vous voulez passer une valeur A du
21
fichier de commandes la fonction, vous devez normalement passer par une entre de la fonction.
Cependant, on aimerait bien parfois que la variable A du fichier de commandes soit utilisable directement par
la fonction.
Pour cela on utilise la directive global. Prenons un exemple : vous disposez d'une relation donnant le Cp d'un
gaz en fonction de la temprature T :
Cp(T) = AT3 + BT2 + CT + D
et vous voulez en faire une fonction. On pourrait faire une fonction 5 entres, pour T, A, B, C, D, mais il est
plus naturel que cette fonction ait seulement T comme entre.
Comment passer A,B,C,D qui sont des constantes ? Rponse : on dclare ces variables en global tant dans le
fichier de commandes que dans la fonction, et on les affecte dans le fichier de commandes.
Fichier de commandes :
global A B C D
A = 9.9400e-8;
B = -4.02e-4;
C = 0.616;
D = -28.3;
T = 300:100:1000; % Temprature en Kelvins
plot(T, cp(T)) % On trace le CP en fonction de T

La fonction :
function y = cp(T)
global A B C D
y = A*T.^3 + B*T.^2 +C*T + D;

17. Matrice de Vandemonde

Ecrire la fonction v=vnd(x) qui partir du vecteur de composante (x0,...,xn) construit la matrice suivante dite
de Vandermonde :
xn L x2 x 1
0 0 0
xn M x2 x 1
v= 1 1 1
M M M M M

x n
n L x 2 x 1
n n

18. Variable locale

Ecrire un script de nom base contenant la ligne unique :
a=12345;
Excuter base, taper ensuite dans la fentre de commande :
a
Taper dans l'diteur la fonction de nom local contenant les 2 lignes :
function local
a=0
Excuter la fonction local, taper ensuite dans la fentre de commande :
a
Expliquer les rsultats affichs.

19. Variable globale

Ecrire une fonction de nom globa contenant les lignes :
22
function globa
% fonction illustrant la notion de variable globale
% voir aussi globb
global A
A=101010;
Ecrire une fonction de nom globb contenant les lignes :
function globb
% fonction illustrant la notion de variable globale
% voir aussi globa
global A
disp(A)
Executer globa, puis globb. Taper A dans la fentre de commandes. Taper maintenant global A dans la
fentre de commandes, puis A. Expliquer les rsultats affichs. Expliquer la diffrence entre les commandes
clear A et clear global A.

5.2. Structures de contrle

On va parler ici de tests et de boucles. Commenons par les oprateurs de comparaison et les oprateurs
logiques

5.2.1. Oprateurs de comparaison et logiques

Notons tout d'abord le point important suivant, inspir du langage C :
MATLAB reprsente la constante logique "FAUX" par 0 et la constante "VRAIE" par 1. Ceci est
particulirement utile, par exemple pour dfinir des fonctions par morceaux. Ensuite, les deux tableaux
suivants comportent l'essentiel de ce qu'il faut savoir :
Oprateur Syntaxe MATLAB

Operator Description
Operator Description
< Less than
& AND
<= Less than or equal to
> Greater than | OR

>= Greater than or equal to ~ NOT

== Equal to
~= Not equal to

On peut se demander ce qui se passe quand on applique un oprateur de comparaison entre deux tableaux, ou
entre un tableau et un scalaire. L'oprateur est appliqu terme terme. Ainsi :
>> A=[1 4 ; 3 2]
A=
14
32
>> A> 2
ans =
01
10
Les termes de A suprieur 2 donnent 1 (vrai), les autres 0 (faux).
23
Une application importante des oprateurs logiques concerne le traitement automatique des grandes matrices.
Crons une matrice (2X5, petite pour l'exemple):
>> A=[10, 12, 8, 11, 9; 10, 10, 128, 11, 9];
A=
10 12 8 11 9
10 10 128 11 9
L'lment 128 s'carte fortement de toutes les autres. Il peut s'agir par exemple d'un point aberrant apparu
lors d'une mesure exprimentale et on dsire remplacer les valeurs suprieure 20 par la valeur 10. Pour une
petite matrice telle que A, cela est facilement ralisable, il suffit de taper: A(2, 3)=10; Mais cette mthode est
irraliste si on a faire avec une matrice 1000X1000 ! Utilisons plutt:
>> A(A>20)=10;
Et le tour est jou:
A=
10 12 8 11 9
10 10 10 11 9
On peut aussi construire des fonctions par morceaux. Imaginons que l'on veuille dfinir la fonction suivante :
f(x) = sin(x) si x > 0 sinon f(x) = sin(2x)
Voil comment crire la fonction:
>> f = inline('sin(x).*(x>0) + sin(2*x).*(~(x>0))', 'x')
f=
Inline function:
f(x) = sin(x).*(x>0) + sin(2*x).*(~(x>0))
On ajoute les deux expressions sin x et sin 2x en les pondrant par la condition logique dfinissant leurs
domaines de validit. Simple mais efficace ! Reprsentons la fonction ainsi dfinie :
>> x=-2*pi:2*pi/100:2*pi;
>> plot(x,f(x))
donne la courbe suivante gauche:
1

0.8

0.6

0.4

0.2

-0.2

-0.4
Fonction carre
gnre partir d'un sinus et
-0.6 d'un
oprateur logique
-0.8

-1

0 2 4 6 8 10 12 14 16

20. Utilisation d'oprateurs logiques

A partir d'une fonction sinus, crer une fonction carre telle que celle, ci-dessus droite.

5.2.2. La commande find

La commande find est utile pour extraire simplement des lments d'un tableau selon un critre logique
donn.
>> k = find (x)
renvoie dans la variable k la liste des indices du tableau x dont les lments sont non nuls.

24
Sachant que "FAUX" et 0 sont identiques dans MATLAB, cela permet de trouver tous les indices d'un
tableau respectant un critre logique donn :
>> x = [ -1.2, 0, 3.1, 6.2, -3.3, -2.1]
x=
-1.2000 0 3.1000 6.2000 -3.3000 -2.1000
>> inds = find (x < 0)
inds =
1 5 6
On peut ensuite facilement extraire le sous-tableau ne contenant que les lments ngatifs de x en crivant
simplement :
>> y = x(inds)
y=
-1.2000 -3.3000 -2.1000

Retenez bien cet exemple. Ce type de squence s'avre trs utile dans la pratique car ici encore, il conomise
un test (if). Notons enfin que la commande find s'applique aussi aux tableaux 2D.
21. Reshape et find
Expliquer prcisment ce que fait la commande suivante :
>>a=reshape(fix((rand(1,400)-0.5)*100),20,20);
Trouver tous les lments de valeur gale 7.
Transformer la matrice a en un vecteur colonne b, trier b par ordre croissant. Combien d'lments ont pour
valeur 7?
Indices : find, sort.
ATTENTION : moins d'ajouter un commentaire explicite, ce genre d'expression est proscrire car il rend
la comprhension du code difficile.

5.2.3. Instructions conditionnelles if

La syntaxe est la suivante :
if condition logique
instructions
elseif condition logique
instructions
...
else
instructions
end
On remarquera qu'il n'y a pas besoin de parenthses autour de la condition logique. Notons qu'il est souvent
possible d'viter les blocs de ce type en exploitant directement les oprateurs de comparaison sur les tableaux
(voir la fonction carre et la fonction par morceaux dfinie dans la section prcdente), ainsi que la
commande find.
22. Fonction avec if et is*
Un certain nombre de fonctions dtecte l'tat des entits MATLAB. Crer une fonction parite(x) qui renvoie
la parit du nombre x ('x est pair' ou 'x est impair'), prcise s'il s'agit d'un nombre premier ('x est premier'),
traite le cas o x n'est pas un nombre et celui o x n'est pas entier (en utilisant la fonction estentier). La
fonction devra utiliser les fonctions isnumeric, isprime et les commandes de contrles if...elseif...end.

25
 5.2.4. Boucles for
for variable = valeur dbut: pas: valeur fin
instructions
end

Exemple:
for i = 1:m
for j = 1:n
H(i,j) = 1/(i+j);
end
end
23. Boucle conditionnelle, break
Combien de fois cette boucle est-elle value ? Quelles sont les valeurs de i affiches l'cran ?
for i=12:-2:1;
disp(i)
if i==8
break
end
end

24. Tic et toc

Ecrire un script utilisant les fonctions tic et toc et une (ou des) boucle(s) telle(s) que le temps coul soit de
l'ordre de la seconde. Le rsultat est-il toujours le mme ?

5.2.5. Boucles while

while condition logique
instructions
end
25. Boucle tant que
Que fait le script suivant ?
eps = 1;
while (1+eps) > 1
eps = eps/2;
end
eps = eps*2

5.2.6. return
Un retour la fonction appelante ou la fentre de commande, est provoqu lorsque MATLAB rencontre la
commande return.

5.2.7. Switch, case, otherwise

26. Switch
Que fait le script ci-dessous ?
clc
couleur = input('Entrer une couleur : ', 's');
couleur=lower(couleur);

switch couleur
26
 case 'rouge'
chaine='La couleur choisie est : rouge';
case 'vert'
chaine='La couleur choisie est : vert';
case 'bleu'
chaine='La couleur choisie est : bleu';
otherwise
chaine='Couleur inconnue';
end
disp(chaine)
27. Boucles implicites
Ecrire explicitement le vecteur u dfini par des boucles implicites :
u(1:2:10)= (1:2:10).^2;
u(2:2:10)= (2:2:10).^3;
Expliquer les rsultats.

5.3. Types de donnes (data types)

5.3.1. Types numriques
Les entiers sont cods sur 8, 16 32 ou
64 bits, il peut-tre avantageux de les utiliser
pour conomiser de la mmoire. C'est le cas
pour les images dont les niveaux d'intensit des
pixels sont reprsents en entiers non-signs sur
8 bits.

Essayer et interprter ceci:

r=2^7*rand(100,100,3);
r=uint8(r);
image(r);

Les nombres virgules flottantes double-prcision (cods sur 64 bits) sont les nombres utiliss par
dfaut par MATLAB. En simple-prcision, ils sont cods sur 32 bits.

MATLAB manipule les nombres complexes.

inf = infini et NaN = Not a Number sont 2 valeurs spciales utiles, reprsentant par exemple 1/0 et
0/0 (ou encore inf/inf) respectivement.

5.3.2. Type logique

Les variables logiques ont dj ts abordes plus haut. Elles ne prennent que 2 valeurs 0 ou 1, "faux" ou
"vrai" logique cod sur 1 octet. Essayer :
>> [30 40 50 60 70] > 40
ans =
0 0 1 1 1
Les oprateurs logiques suivants retournent un type logique :
& (and), | (or), ~ (not), xor, any, all
27
Ils peuvent tre utiliss lorsque l'on veut excuter une portion de code conditionnel , extraire des portions de
matrices (voir exercice correspondant).

5.3.3. Caractres et chanes

Les chanes sont traits comme des tableaux.

5.3.4. Structures
Les structures (bases de donnes) sont des tableaux avec des "conteneurs de donnes" appels "champs". Les
champs peuvent tre de n'importe quel type. Par exemple, un champ peut contenir une chane de caractre
reprsentant un nom, un autre contenir un nombre reprsentant l'ge, un troisime une matrice de notes, ...
Exemple :
etudiant.nom = 'Jean Durand';
etudiant.age = 21;
etudiant.notes = [12; 13; 10; 14; 14];
etudiant est une structure 3 champs ne comportant pour l'instant qu'un seul lment. Ajoutons en un second:

etudiant(2).nom = 'Gabriel Dupont';

etudiant(2).age = 22;
etudiant(2).notes = [19; 18; 18; 17; 19];

Les structures sont particulirement utiles pour passer en une seule variable un grand nombre de paramtres.

5.3.5. Cellules
Un tableau de cellules permet de stocker des donnes de type diffrents. La figure ci-dessous montre un
tableau A (2 X 3) de cellules contenant: un tableau d'entiers non signs (cell 1,1), un tableau de caractres
(cell 1,2), un tableau de nombres complexes (cell 1,3), un tableau de nombres virgules flottantes (cell 2,1),
un tableau d'entiers signs (cell 2,2), et un autre tableau de cellules (cell 2,3).

28
Pour accder aux donnes d'un tableau de cellules, utiliser la mme mthode d'indexation que pour les autres
matrices MATLAB avec des accolades { } la place des parenthses ( ). Par exemple A{1, 3} donne le
contenu de la cellule place sur la premire ligne, la troisime colonne, ici le tableau de nombre complexes.
Pour crer la cellule (1, 3), on peut crire:
A(1, 3)={0.25+3i, 8-16i; 34+5i, 7+0.92i}; ou
A{1, 3}=[0.25+3i, 8-16i; 34+5i, 7+0.92i];

Les tableaux de cellules peuvent tre une alternative aux structures.

28. Fonctions signenombre et signematrice
En utilisant la structure de contrle if...elseif...else...end, crire la fonction signenombre renvoyant la valeur
-1, 0 ou 1 suivant que le nombre pass en entre est ngatif, nul ou positif respectivement.
Ecrire maintenant une fonction signematrice s'appliquant aux matrices et renvoyant une matrice de -1, 0

29. Ecrire une fonction dont l'en-tte est :

function i=recherche(v,x)
recherchant dans un vecteur v la position du dernier lment gal une valeur x donne (ou 0 si il est
absent).

30. Fonction estentier

Ecrire une fonction { function r=estentier(x) } utilisant fix pour tester si un nombre x est entier (r=1) ou non
(r=0). La fonction doit pouvoir s'appliquer des matrices.
31. Types numriques
Dans la fentre de commandes, taper
a=65;
b=int8(65);
puis comparer la taille mmoire occupe par ces 2 variables, c'est- dire le nombre de bytes (octets) utiliss
pour les coder indiqus dans la fentre "Workspace".
32. Type logique
Que fait le script suivant :
a=input('Entrer un nombre entier ');
if ~isprime(a)
disp([num2str(a),' n''est pas premier'])
return
end
disp([num2str(a),' est premier'])
Celui-ci ?
A = 5:5:50
B = [1 3 6 7 10]
A(B)
Essayer et expliquer ceci :
A = fix(25*rand(5))
B = A>20
A(~B) = 0
33. Type caractres et chanes de caractres
Concatner les 2 chanes : a='abc' et z='uvwxyz' en ligne chligne=[a,z], puis en colonne chcolonne=[a;z].
La commande char, ajoute des espaces pour adapter la diffrence de longueurs des chanes. Essayer
char(abc, z).

29
Transformer chligne en entiers correspondant aux codes ASCII de ses caractres : asc=uint8(chligne).

Revenir aux caractres avec la commande char(ans)

Rechercher la position du caractre 'w' : position = findstr('w', chligne)

Que fait cette ligne de commande : dec2bin(a(findstr('w', a)))

34. Types structures et cellules

Cration d'une structure fictive etudiant, comportant les champs : nom, prenom, age, notes et moyenne.

1/ Ecrire une fonction enAlea=abaleatoire(a, b) qui renvoie un entier alatoire compris entre a et b.
Cette fonction servira dterminer de faon alatoire
- la longueur de chanes de caractres reprsentant les noms et les prnoms d'une liste de personnes fictives,
cette longueur devra tre comprise entre 5 et 12
- l'ge qui devra tre compris entre 16 et 21

2/ Ecrire une fonction chAlea=chainealeatoire(n) qui renvoie une chane de caractres alatoires de
longueur n. Les codes ASCII des lettres a, b, c, ...,z s'tendent de 97 122. Les noms et prnoms des
tudiants seront gnrs par cette fonction.

3/ Ecrire une fonction noAlea=notesaleatoires qui renvoie une matrice alatoire (5X1) de nombre compris
entre 0 et 20 et dfinis une dcimale prs.

4/ Ecrire la fonction principale (principale dans le sens o elle utilise les fonctions prcdentes)
etStruct=generationetudiant(n) qui gnre la structure etudiant comportant n lments (donc n tudiants).

5/ Ecrire la fonction afficheunetudiant(s) qui indique l'utilisateur le nombre d'lments de la structure s,

l'invite entrer un numro d'lment et affiche l'cran le nom et le contenu des champs de cet lment.

6. Conseils dans la programmation

Conventions largement adopts par les programmeurs, mme dans d'autres langages. Elles ne sont
videmment pas obligatoires mais permettent une meilleure lisibilit et assurent l'entretien des programmes.

6.1. Conventions d'criture

* Les noms de variables devraient documenter leur sens ou usage, commencer par une minuscule avec
ventuellement, des majuscules isoles pour assurer leur lisibilit (pas d'accentuation):
masse2
poidsVehicule

* La longueur d'un nom est li la porte de la variable. Des variables entires comme i, j, k, n, ... ou relles
comme x, y, u, ... ont une porte de quelques lignes de code.

* Le prfixe n devrait tre rserv pour indiquer des nombres :

nPoints
nFichiersOuverts

* Les itrateurs peuvent tre prfixs par i :

30
for iFichier=1:nFichier
:
end

* Les noms de structures devraient commencer par une majuscule:

Segment.longueur
Segment.epaisseur

* Les noms de fonctions devraient tre intgralement en minuscules avec ventuellement des "underscores":
calcule_moyenne(.)
triangule(.)
extraitcolonne(.)

* Les constantes devraient tre en majuscules:

BOLTZMANN=1.38e-23; % constante de Boltzmann
LUMIERE=3e8;% vitesse de la lumire

6.2. Fichiers, organisation, amlioration des performances

* Le meilleur moyen d'crire un gros programme est de l'assembler partir de petites pices, gnralement des
fonctions spcialises dans une tche spcifique.

* Une fonction interagit avec les programmes partir de ses entres, de ses sorties et des variables globales. Les
structures permettent d'viter de longues listes de variables d'entre/sortie.

* Vectoriser les boucles, c'est dire les convertir en oprations sur des vecteurs ou des matrices. L'exemple suivant
calcule le sinus de 1001 valeurs comprises entre 0 et 10:
i = 0;
for t = 0:.01:10
i = i + 1;
y(i) = sin(t);
end
La version vectorise du mme programme est:
t = 0:0.01:10;
y = sin(t);
Comparer le temps d'excution avec tic et toc.

* Pr-allouer de l'espace mmoire pour les matrices. Dans les boucles for ou while qui augmente la taille
d'une structure chaque boucle, MATLAB perd du temps chaque fois rechercher des blocs contigus de
mmoire de plus en plus grands. La rapidit du programme sera amlior si la matrice est initialise ds le
dbut sa taille maximale. Voici une version dispendieuse en temps d'excution:
x = 0;
for k = 2:1000
x(k) = x(k-1) + 5;
end
Changer la premire ligne en allouant un bloc de mmoire de 1 X 1000 pour x initialis zro produit un
srieux gain de temps machine:
x = zeros(1, 1000);
for k = 2:1000
x(k) = x(k-1) + 5;
end

31
 7. Graphiques 3D
7.1. Courbes en 3D
Une courbe en 2D est dfinie par une liste de doublets (x; y) et une courbe en 3D par une liste de triplets (x;
y; z). Puisque l'instruction plot attendait deux arguments, un pour les x, un pour les y, l'instruction plot3 en
attend trois : les x, les y et les z.
Voici un exemple de courbe paramtre :
>> t=-10:0.1:10;
>> plot3(exp(-t/10).*sin(t), exp(-t/10).*cos(t), exp(-t))
>> grid

7.2. Surfaces
7.2.1. Gnration des points (meshgrid)
Pour dfinir une surface, il faut un ensemble de triplets (x; y; z). En gnral les points (x; y) tracent dans le
plan un maillage rgulier mais ce n'est pas une obligation. La seule contrainte est que le nombre de points
soit le produit de deux entiers mn (on comprendra pourquoi trs bientt).
Si on a en tout mn points, cela signifie que l'on a mn valeurs de x, mn valeurs de y et mn valeurs de z. Il
apparat donc que abscisses, ordonnes et cotes des points de la surface peuvent tre stockes dans des
tableaux de taille mn.
Toutes les instructions de trac du surface, par exemple surf auront donc la syntaxe suivante:
surf ( tableau d'abscisses, tableau d'ordonnes, tableau de cotes)

x L x 1n y L y1n z L z1n
11 11 11
x 21 M x 2n y 21 M y 2n z 21 M z 2n

M M M M M M M M M

x m1 L x mn y m1 L y mn z m1 L z mn

Il reste maintenant construire ces tableaux. Prenons l'exemple de la surface dfinie par z = x2 + y2 dont on
veut tracer la surface reprsentative sur [-1; 1] [-2; 2].
Pour dfinir un quadrillage de ce rectangle, il faut d'abord dfinir une suite de valeurs x1; ... ; xm pour x et une
suite de valeurs y1; ... ; yn pour y, par exemple :
>> x = -1:0.2:1

32
x=
Columns 1 through 7
-1.0000 -0.8000 -0.6000 -0.4000 -0.2000 0 0.2000
Columns 8 through 11
0.4000 0.6000 0.8000 1.0000
>> y = -2:0.2:2
y=
Columns 1 through 7
-2.0000 -1.8000 -1.6000 -1.4000 -1.2000 -1.0000 -0.8000
Columns 8 through 14
-0.6000 -0.4000 -0.2000 0.2000 0.4000 0.6000
Columns 15 through 21
0.8000 1.0000 1.2000 1.4000 1.6000 1.8000 2.0000
En combinant toutes ces valeurs de x et y, on obtient mn points dans le plan (x; y). Il faut maintenant
construire deux tableaux, l'un contenant les 1121 abscisses de ces points l'autre les 1121 ordonnes, soit:

1 0.8 L 1 2 2 L 2

1 0.8 L 1 1.8 1.8 L 1.8
X = Y =

M M M M M M M M

1 0.8 L 1 2 2 L 2
Ainsi, si on considre par exemple la premire colonne de ces matrices, on obtient tous les couples (x, y)
pour x=-1: (-1, -2), (-1, -1.8), (-1, -1.6), (-1, -1.4), ..., (-1, 1.8) et (-1, 2). Les autres colonnes couvrent le reste
de la portion du plan xy. Ci-dessous, le maillage correspondant:

2
1.8
1.6
1.4
1.2
1
0.8
0.6
0.4
0.2
Y 0
-0.2
-0.4
-0.6
-0.8
-1
-1.2
-1.4
-1.6
-1.8
-2
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1

33
Rassurez-vous pas besoin de boucles, la fonction meshgrid se charge de gnrer ces 2 matrices:
>> [X,Y] = meshgrid(x, y);
Il reste maintenant calculer les z = f(x, y) correspondants. C'est l que les calculs terme terme sur les
matrices montrent encore leur efficacit : on applique directement la formule aux tableaux X et Y, sans
oublier de mettre un point devant les oprateurs *, / et ^
>> Z = X .^2 + Y.^2;

7.2.2. Trac de la surface

Ensuite on peut utiliser toutes les fonctions de trac de
surface, par exemple mesh :
>> mesh(X,Y,Z)

Les instructions les plus courantes sont

- mesh qui trace une srie de lignes entre les points de la surface
- meshc qui fonctionne comme mesh mais ajoute les courbes de niveau dans le plan (x; y)
- surf qui "peint" la surface avec une couleur fonction de la cote
- surfl qui "peint" la surface comme si elle tait claire.
- surfc qui fonctionne comme mesh mais ajoute les courbes de niveau dans le plan (x, y)

Notons enfin qu'il existe des fonctions de conversion entre les coordonnes cartsiennes, cylindriques et
sphriques, permettant de tracer facilement des courbes dfinies dans l'un de ces systmes de coordonnes.
On regardera par exemple la documentation de cart2pol.

7.2.3. Courbes de contour

On utilise la fonction contour. Elle s'utilise comme les
instructions prcdentes, mais fournit un graphe 2D dans le
plan (x; y). Plusieurs paramtres optionnels peuvent tre
spcifis, notamment le nombre de courbes de contours
afficher :
>> contour(X,Y,Z, 10)

7.2.4. Contrle de l'angle de vue

Il existe la commande view, mais le plus simple est de
travailler partir de la fentre graphique, cliquer sur .
Cliquez ensuite avec le bouton gauche de la souris pour faire tourner la figure. ATTENTION : tout clic de la
souris provoque un raffichage, et si votre surface contient beaucoup de points, elle met beaucoup de temps
se rafficher.
35. Trac d'une surface en 3D
Tracer la surface z=sin(x2 + y2) pour 0<x, y<4sqrt(pi).

34
 8. Echanges entre MATLAB et l'extrieur
8.1. Sauvegarde de donnes
Commandes save et load. La premire permet d'crire toute ou une partie des variables dans un fichier
binaire dont l'extension est ".mat". La seconde permet de recharger les variables contenues dans le fichier.
Syntaxe :
save nom fichier var1 var2 var3
Les variables sont simplement spares par des blancs. Si aucune variable n'est spcifie, toutes les variables
sont sauvegardes dans nom_fichier.mat. La syntaxe de load est identique. Les donnes sont dans un format
binaire.
Noter que des options de ces commandes permettent de sauver et de charger des fichiers de donnes dans le
format ASCII (format texte ==> voir l'aide).

8.2. Importer des tableaux

MATLAB est souvent utilis comme logiciel d'exploitation de rsultats. Exemple classique: on a un
programme FORTRAN ou C ou ... qui calcule des sries de valeurs, et on souhaite tracer une srie en
fonction de l'autre.
Le moyen le plus simple est que votre programme crive un fichier texte organis sous forme de n colonnes,
spares par des blancs (espaces), contenant chacune m lignes:
x1 y1 z1 : : :
x2 y2 z2 : : :
x3 y3 z3 : : :
:::
xm ym zm : : :
Cette structure est obligatoire pour que la mthode simple expose ici fonctionne correctement. Ce fichier
doit porter une extension, qui peut tre quelconque, mais diffrente de ".m" ou ".mat". Admettons qu'il soit
nomm toto.res, tapez dans MATLAB :
>> load toto.res
Ce faisant, vous crez un nouveau tableau MATLAB, qui porte le nom du fichier sans extension, soit toto.
Vous pouvez maintenant extraire les colonnes de ce tableau (voir section 3.3), et les mettre dans des vecteurs,
par exemple :
>> x = toto(:,1)
>> y = toto(:,2)
>> z = toto(:,3)
et vous n'avez plus qu' tracer l'un de ces vecteurs en fonction d'un autre.

8.3. Inclure des courbes MATLAB dans un document

Dans la fentre graphique, choisissez File -> Export. Vous pouvez alors crer toutes sortes de fichier
graphiques avec le contenu de la fentre graphique.
Pour les documents du pack office de Windows, vous pouvez aussi de faon assez simple, cliquer dans la
fentre graphique sur edit, copy figure, puis coller dans votre document comme cela a t ralis dans celui-
ci.

9. Calcul numrique avec MATLAB

L'objectif de ce chapitre est de prsenter brivement des routines de calcul numrique fournies dans la
version de base de MATLAB. La liste des problmes fournie dans ce polycopi est non exhaustive et on
trouvera des informations complmentaires dans les diffrentes documentations en ligne MATLAB. De
mme, les concepts thoriques sont supposs connus par le lecteur et il ne sera pas fait mention des mthodes
utilises.

35
 9.1. Recherche des zros d'une fonction
Problme : On cherche x0 tel que f(x0) = 0.
La fonction fzero permet de rsoudre ce problme. Il faut fournir d'une part la fonction f elle-mme, et
d'autre part une estimation (valeur initiale) de x0. L'efficacit de l'algorithme est comme toujours dpendante
de la valeur initiale choisie. La fonction f peut-tre dfinie par une directive inline ou bien crite dans une
fichier.
Par exemple on cherche le zro de :
f(x) = cos(x) - x
Une approche graphique permet souvent de trouver une estimation de x0. La figure suivante montre ainsi que
les fonctions x x et x cos(x) se coupent en un point sur [-; ]. Une valeur raisonnable pour l'estimation
de x0 est par exemple 0.

On peut donc crire, en utilisant inline :

>> f = inline('x-cos(x)')
>> fzero(f,0)
ans =
0.7391
On remarquera que la variable f envoye la fonction fzero est elle-mme une fonction. Toutes les routines
de calcul numrique de MATLAB ncessitant l'criture d'une fonction par l'utilisateur fonctionnent selon ce
principe.
On peut galement crire la fonction dans un fichier f.m :
function y = f(x)
y = x-cos(x);
et ensuite on crira :
>> fzero('f',0)
ans =
0.7391
ATTENTION! On remarquera ici que le symbole f est mis entre primes. Cela vient du fait qu'ici, la dfinition
de la fonction f est faite dans un fichier. Il est prfrable en gnral de dfinir les fonctions dans des fichiers
car cela offre une plus grande souplesse.
Un dernier point important : comment faire lorsque la dfinition de la fonction dpend en plus d'un paramtre
? Par exemple, on veut chercher le zro de la fonction f(x) = cos(mx) - x
o m est un paramtre susceptible de varier entre deux excutions. On ne peut pas rajouter un argument la
dfinition de notre fonction f car fzero impose que f ne dpende que d'une seule variable. Une solution serait
d'affecter m dans le corps de la fonction mais elle est mauvaise car :
- lorsque l'on veut changer la valeur de m il faut modifier la fonction,
- cette fonction sera peut-tre appele des dizaines voire des centaines de fois par fzero, et on rptera
chaque fois la mme instruction d'affectation de m.

36
La bonne solution est d'utiliser la directive global vue pcdemment. La fonction f s'crira donc :
function y = f(x)
global m
y = x-cos(m*x);
et on cherchera son zro en crivant :
>> global m
>> m = 1;
>> fzero(f,0)
Notons enfin que l'on aura tout intrt mettre les trois lignes ci-dessus dans un fichier de commandes, et
lancer ce fichier de commandes en bloc.

9.2. Interpolation
Le problme : connaissant des points tabuls (xi; yi), construire une fonction polynomiale passant par ces
points. En termes plus simples, c'est ce que fait la fonction plot quand elle "relie" les points que vous lui
donnez : elle fait passer un polynme de degr 1 entre deux points conscutifs.
La fonction interp1 rsout le problme. Il suffit de lui donner les valeurs tabules ainsi qu'une srie
d'abscisses o l'on souhaite connatre les ordonnes interpoles. Sur chaque intervalle on peut interpoler par
un polynme de degr 0 (option 'nearest'), de degr 1, qui est le comportement par dfaut (option 'linear') ou
de degr 3 (option 'cubic').
Les lignes suivantes fournissent un exemple gnrique: partir de 8 points (xi, yi), on cherche une
interpolation f(x) entre ces 8 points.
clear; close all;
xi = [-1 -1/2 0 1/2 1 2 3 4];
yi = [1.5 1 10 1 1/2 10 -10 0];
x = -2:0.1:5;
ycubic = interp1 (xi, yi, x, 'cubic');
plot(xi,yi,'ok', x,ycubic,'-k')
La figure suivante illustre l'interpolation (en traits pleins) cubiques
obtenues par interp1 sur ces 8 points (illustrs par des cercles).

9.3. Approximation (estimation de paramtres ou "fitting" ou "ajustement")

Problme : faire passer une courbe d'quation connue au milieu d'un nuage de points de telle sorte que la
courbe soit "la plus proche possible" de l'ensemble des points. Comme on cherche minimiser une distance
entre deux fonctions, on parle en gnral d'approximation aux moindres carrs ou encore de rgression.

37
 9.3.1. Linaire
Illustrons la mthode par un exemple: on cherche ajuster des points exprimentaux (xi; yi) par une fonction
de la forme :
a2 a
y = f (x) = a 1 + + 32
x x
Ou encore de faon plus gnrale, la forme suivante : y = a1 f1(x) + a2 f2(x) + ... +am fm(x).
Les 5 points exprimentaux sont donns par : xi = [0.25 1/2 1 2 3] et yi = [18 7 4 2 1.1].
On a donc un systme de 5 quations 3 inconnues (le systme est dit "surcontraint", c'est dire qu'il y a
plus d'quations que d'inconnues) :

1 1 1
x x2

y 1 1
1 1
x x2
y a
2 1 1 1
y =
1
x x2 a2
3
y4 1 1 a3

1
x
x2
y5
1 1 1
x x2

de la forme Y = F A. On trouve la matrice des paramtres par inversion : F-1 Y = A qui s'crit sous
MATLAB par F \ Y = A.
La squence d'instructions est la suivante :
>> xi = [1/4 1/2 1 2 3]; xi = xi'; % on dfinit les abscisses et ordonnes
>> yi = [18 7 4 2 1.1]; yi = yi'; % en colonnes
>> F = [ones(size(xi)) 1./xi 1./(xi.^2)];
>> A = F \ yi
>> x = 1/4:1/16:3; % On trace la courbe obtenue
>> plot(xi,yi,'o', x, A(1) + A(2)./x + A(3)./(x.^2))
Ce programme trace la courbe ajuste ainsi que les points exprimentaux :

38
Cette mthode marche ds que la fonction cherche dpend linairement des coefficients inconnus (le cas
oppos dpasse le cadre d'une simple introduction). Cela fonctionne en particulier pour faire de
l'approximation polynomiale, mais dans ce cas on peut utiliser directement la fonction polyfit, qui en plus
des coefficients du polynme renvoie des informations sur la qualit de l'approximation ralise. La fonction
polyfit attend en entre simplement les xi, les yi et l'ordre du polynme recherch. La fonction polyval peut -
tre utilise pour recalculer le polynme d'approximation en tout point.
L'exemple suivant calcule deux polynmes d'approximation, un d'ordre 3, l'autre d'ordre 4 et trace le rsultat.
>> xi = [-2 -1.5 -1 0 1 1.5 2];
>> yi = [3.5 2 1.5 0 1.3 2.5 3.9];
>> A4 = polyfit (xi, yi, 4)
>> A3 = polyfit (xi, yi, 3)
>> x = -2:0.1:2;
>> plot(xi, yi, 'o', x, polyval(A4,x), x, polyval(A3,x), '--' )

A4 =
-0.0652 0.0313 1.1411 -0.0022 0.1473

A3 =
0.0313 0.8570 -0.0022 0.3247

36. Simulation d'une exponentielle dcroissante bruite

On dsire simuler la mesure d'une grandeur exprimentale pseudo-priodique.
a/ Ecrire une fonction { function s=singau(param, t) } qui renvoie une sinusodale de la forme gnrale :
s(t) = a * exp(b * (t td)2) * cos(2 * * f * t + phi)
param sera un vecteur avec la correspondance :
param(1)= a
param(2)= b
param(3)= td
param(4) = f
param(5) = phi
t est un vecteur temps de N points ( t = 0, t, 2t, 3t, ..., (N-1)t ) cr pralablement, a l'amplitude, b un
coefficient d'amortissement, td un dcalage temporel, f la frquence et phi la phase initiale. La fonction sera
correctement documente.
b/ Ecrire une fonction { function yb=bruit(y, sig) } qui additionne un bruit normal d'cart type sig un
vecteur quelconque y. Ce bruit simulera celui toujours prsent lors des mesures exprimentales.

39
c/ Appliquer la fonction bruit s (Une valeur sig = 0.1*max(s) parat raisonnable et signifie qu'on commet
une erreur accidentelle de 10% du calibre (=max(s)) sur les mesures), comparer dans le mme graphe les 2
courbes s et bruit(s, 0.1*max(s)). Ci-dessous une partie de la figure qu'on peut obtenir.

40