Notes introductives Matlab

Le logiciel Matlab consiste en un langage interprt qui s'excute dans une fentre dite d'excution. L'intrt de Matlab tient, d'une part, sa simplicit d'utilisation : pas de compilation, dclaration implicite des variables utilises et, d'autre part, sa richesse fonctionnelle : arithmtique matriciel et nombreuses fonctions de haut niveau dans de nombreux domaines (analyse numrique, graphique, ...). La programmation sous Matlab consiste crire des scripts de commandes Matlab, excutables dans la fentre d'excution. En outre, grce aux diverses Toolboxes spcialiss (ensemble de scripts Matlab), Matlab s'enrichit au fur et mesure.

1. Lancement de Matlab
Lors de son lancement (via la commande Matlab par exemple) la fentre d'excution s'ouvre. Il est alors possible d'excuter diffrents types de commandes dans cette fentre, par exemple1 :

help HELP topics: toolbox\local matlab\datafun matlab\elfun matlab\elmat matlab\funfun matlab\general matlab\color matlab\graphics matlab\iofun matlab\lang matlab\matfun matlab\ops matlab\plotxy matlab\plotxyz matlab\polyfun matlab\sounds matlab\sparfun matlab\specfun matlab\specmat matlab\strfun matlab\dde matlab\demos toolbox\wintools Local function library. Data analysis and Fourier transform functions. Elementary math functions. Elementary matrices and matrix manipulation. Function functions - nonlinear numerical methods. General purpose commands. Color control and lighting model functions. General purpose graphics functions. Low-level file I/O functions. Language constructs and debugging. Matrix functions - numerical linear algebra. Operators and special characters. Two dimensional graphics. Three dimensional graphics. Polynomial and interpolation functions. Sound processing functions. Sparse matrix functions. Specialized math functions. Specialized matrices. Character string functions. DDE Toolbox. The MATLAB Expo and other demonstrations. GUI tools for MATLAB for MS Windows.

For more help on directory/topic, type "help topic".

La commande suivante permet d'affiner l'aide sur les fonctions mathmatiques lmentaires :

dans tous les exemples qui suivent, apparat aprs la commande le rsultat de cette commande 1

help elfun Elementary math functions. Trigonometric. sin sinh asin asinh cos cosh acos acosh tan tanh atan atan2 atanh sec sech asec asech csc csch acsc acsch cot coth acot acoth Exponential. exp log log10 sqrt Complex. abs angle conj imag real Numeric. fix floor ceil round rem sign Sine. Hyperbolic sine. Inverse sine. Inverse hyperbolic sine. Cosine. Hyperbolic cosine. Inverse cosine. Inverse hyperbolic cosine. Tangent. Hyperbolic tangent. Inverse tangent. Four quadrant inverse tangent. Inverse hyperbolic tangent. Secant. Hyperbolic secant. Inverse secant. Inverse hyperbolic secant. Cosecant. Hyperbolic cosecant. Inverse cosecant. Inverse hyperbolic cosecant. Cotangent. Hyperbolic cotangent. Inverse cotangent. Inverse hyperbolic cotangent. Exponential. Natural logarithm. Common logarithm. Square root. Absolute value. Phase angle. Complex conjugate. Complex imaginary part. Complex real part. Round towards zero. Round towards minus infinity. Round towards plus infinity. Round towards nearest integer. Remainder after division. Signum function.

permettant ainsi de voir toutes les fonctions mathmatiques lmentaires dont dispose Matlab. On peut maintenant prciser la recherche si l'on veut avoir une ide plus prcise de la fonction log par exemple :
help log LOG Natural logarithm. LOG(X) is the natural logarithm of the elements of X. Complex results are produced if X is not positive. See also LOG10, EXP, LOGM.

Quelques commandes d'environnement importantes

path : permet de savoir quels sont les dossiers que Matlab rfrence et de spcifier de nouveaux dossiers Unix o se trouvent des ressources personnelles, par exemple : path=(path,'~/monapplication/mesfichiersmatlab') cd : positionne Matlab dans un dossier Unix, par exemple : cd ~/monapplication/mesfichiersmatlab dir ou ls : permet de faire la liste des objets du dossier courant, par exemple la commande suivante :
ls . basename.exe dos4gw.exe .. binpatch.m fmex.bat basefnam.exe cmex.bat is_ext.exe matlab.exe matlab.hlp ml_bang.pif ml_dos.pif spr.exe

clear all : efface tous les objets en mmoire clf : dtruit les figures Si monscript.m est un script Matlab que vous avez crit et qui est accessible (soit via le path ou un commande cd) pour l'excuter, il suffira de saisir la commande : monscript

Par dfaut, Matlab inscrit le rsultat la suite de la commande. Si l'on ne dsire pas voir le rsultat d'une commande il suffit de terminer cette commande par ;. Pour certaine commande (ls,help,...) cel ne change rien!

2. Les variables sous Matlab

Matlab gre de faon automatique : les nombres entiers, rels, complexes de faon indiffrente, les chaines de caractres ainsi que les tableaux de nombre. En aucun cas, il n'est utile de dclarer le type de la variable que l'on manipule, y compris les tableaux. Ainsi les instructions suivantes, dclarent les variables lors de leur affectation :
a=1 a = 1 b=1.01 b = 1.0100 X=1.0e+05 X = 100000

nom=' mon nom' nom = mon nom

c=1+2i c = 1.0000 + 2.0000i

la constante i est le nombre imaginaire prdclar, de mme que certaines constantes (e,pi,...) On dclare un vecteur colonne de la faon suivante :
u=[1;3;-1] u = 1 3 -1

un vecteur ligne de la faon suivante :

v=[1,3,-1] v = 1 3

-1

et une matrice d'ordre 3x2 :

A=[1,2 ; -1, 3; 4, 0] A = 1 2 -1 3 4 0

la ',' sert sparer les lments d'une ligne et ';' les lments colonnes. En fait, on peut remplacer la ',' par un espace , ce qui amliore la lisibilit :
vb=[1 3 -1] vb = 1 3

-1

Pour spcifier un lment d'un vecteur, d'une matrice, on utilise la syntaxe suivante :
u(2) ans = 3 v(3) ans = -1 A(3,2) ans = 0

L'utilisation d'indice hors limite provoque une erreur, comme le montre cet exemple:
A(3,3) ??? Index exceeds matrix dimensions.

On peut utiliser des raccourcis bien utiles et plus efficaces pour remplir des vecteurs ou des tableaux . En voici quelques exemples :
u1=1:10 u1 =

(incrmentation automatique de 1 10 avec pas de 1)

1 v1=1:2:10 v1 = 1

10

(incrmentation automatique de 1 10 avec pas de 2)

3 5 7 9

Id3=eye(3) Id3 = 1 0 0 1 0 0 Un=ones(2) Un = 1 1 1 1 Z=zeros(2,3) Z = 0 0 0 0

(matrice identit d'ordre3)

0 0 1

(matrice constitue de 1 d'ordre 2)

(matrice nulle d'ordre 2x3)

0 0

De mme, il existe des syntaxes particulires permettant d'extraire des lignes, des colonnes de matrices :
A1=[11 12 13;21 A1 = 11 12 21 22 31 32 A1(:,1) ans = 11 21 31 A1(2,:) ans = 21 22 23;31 32 33] 13 23 33

(colonne 1 de la matrice A1)

(ligne 2 de la matrice A1)

22 23

Une des difficults que rencontre le dbutant est le mauvais maniement des matrices du une mauvaise connaissance des tailles des matrices. Pour vrifier ces tailles, on pourra utiliser la commande size:
size(A1) ans = 3 size(u) ans = 3

3. Oprations lmentaires sous Matlab

Les oprations sur les scalaires sont standards : addition +, soustraction -, multiplication *, division /, puissance ^. La racine carre s'obtient par la function sqrt.On dispose de toutes les fonctions usuelles sur les scalaires,

faire help elfun pour plus de dtails. Attention les fonctions peuvent renvoyer des complexes mme dans des situations anodines :
sqrt(-1) ans = 0 + 1.0000i acos(2) (function arc cosinus, ici on a affaire au prolongement dans le plan complexe de cette fonction!) ans = 0 + 1.3170i

En ce qui concerne les vecteurs et matrices ces oprateurs se prolongent au sens du calcul vectoriel et matriciel. En particulier, il faut veiller la compatibilit des tailles des objets! Voici quelques exemples :
u=[1 2 3] u = 1

v=[-1 1 1] v = -1 1 w=u+v w = 0 ut=u' ut = 1 2 3 3

(addition)
4

(transposition d'un vecteur ligne)

ut2=[ut ut] (concatnation de deux vecteurs colonnes donne une matrice 3x2) ut2 = 1 1 2 2 3 3 ps=v*ut ps = 4 M=ut*v M = -1 -2 -3

(Produit qui conduit au produit scalaire)

(Produit qui conduit une matrice)

1 2 3 1 2 3

L=M+2*eye(3) L = 1 1 -2 4 -3 3 Y=L\ut Y = 0.1667 0.3333 0.5000

1 2 5

(Rsolution du systme linaire L.Y=ut)

E=u/L' E = 0.1667

(Rsolution du systme linaire E.L' =u)

0.3333 0.5000

On prendra garde au sens de la division. Si les matrices ne sont pas inversibles, un message vous prvient. On peut effectuer des oprations tensorielles sur les vecteurs et matrices, par exemple effectuer le produit de deux vecteurs colonnes, composante par composante par l'adjonction d'un . l'oprande *. Par exemple :
ut.*Y ans = 0.1667 0.6667 1.5000 ut.^Y ans = 1.0000 1.2599 1.7321

De mme, Matlab autorise l'utilisation de toutes les fonctions scalaires dans un contexte vectoriel. Ainsi, si h est un vecteur de dimension n, sin(h) sera un vecteur de mme dimension :
h=0:pi/4:pi h = 0 sin(h) ans = 0

0.7854 0.7071

1.5708 1.0000

2.3562 0.7071

3.1416 0.0000

Pour les multiples oprations sur les matrices (inverse, puissance, trace, dterminant, factorisation, ...) faire help elmat et help matfun. Pour ce qui est des oprations sur les chanes de caractres, ces dernires tant considres comme des vecteurs ligne de caractres ascii, la concatnation de deux chanes s'effectuera de la faon suivante :
c1='texte' c1 = texte c2=' et suite de texte' c2 = et suite de texte c3=[c1 c2] c3 = texte et suite de texte

pour les autres oprations voir help strfun.

4. Contrle de flux
Oprateurs boolens

Avant de dcrire la syntaxe du test sous Matlab, indiquons les principaux oprateurs de relation ainsi que les oprateurs boolens qu'utilisent Matlab.

< <= > >= == ~= & | ~

strictement infrieur infrieur ou gal strictement suprieur suprieur ou gal gal diffrent de et logique (and) ou logique (or) non logique (not)

Le rsultat d'un test est un boolen, qui sous Matlab, prend la valeur 1 pour vrai et 0 pour faux. Par exemple, on a les rsultats suivants :
r=1<2 r = 1 r=~((1>2)|(0~=0)) r = 1

(traduction Matlab de l'expression logique : non (1>2 ou 00))

Il existe d'autres fonctions boolennes, par exemple xor, isfinite, isnan, isinf,..dont on trouvera la description en faisant help ops.

Syntaxe du test (if)

if expression boolenne instructions end if expression boolenne instructions else instructions end if expression boolenne instructions elseif expression boolenne instructions else instructions end

Syntaxe du branchement (switch) (expression est un scalaire ou une chane de caractres) switch expression case value1 instructions (instructions effectues si expression=value1) case value2 instructions ... otherwise instructions end Syntaxe de boucle (while et for)
while expression instructions end for indice=debut:pas:fin instructions end (si le pas n'est pas prcis, par dfaut il vaut 1)

5. Utilisation des fonctions

La notion de fonction existe sous Matlab. Sa syntaxe est la suivante : function [args1,args2,...] = nomfonction(arge1,arge2,...) instructions args1,args2,... sont les arguments de sortie de la fonction et peuvent tre de n'importe quel type arge1,arge2,... sont les arguments d'entre de la fonction et peuvent tre de n'importe que type instructions est un bloc d'instructions quelconque devant affecter les arguments de sortie args1,args2,... Lorsqu'il n'y a qu'un seul argument de sortie, on peut utiliser la syntaxe plus simple : function args = nomfonction(arge1,arge2,...) L'appel la fonction s'opre de la faon suivante : [vars1,vars2,...] = nomfonction(vare1,vare2,...) avec compatibilit des variables d'entres vare1,vare2,... avec les arguments d'entre arge1,arge2,... et compatibilit des variables de sorties vars1,vars2,..., si celles-ci ont dj t utilises, avec les arguments de sortie args1,args2,.... Remarque : il n'est pas obligatoire de fournir tous les arguments d'entres et de sortie lors de l'appel d'une fonction, mais ceux que l'on omet doivent tre les derniers des listes d'entre ou de sortie. Ainsi, supposons nomfonction soit une fonction 2 arguments d'entre et 2 arguments de sortie, on peut alors utiliser l'appel suivant : [vars1]= nomfonction(vare1) mais pas l'appel : [vars2]= nomfonction(vare2).

La limitation du nombre d'arguments de sortie est gre de faon automatique dans Matlab. Par contre, les variations du nombre d'arguments d'entre doivent tre gres par le concepteur de la fonction l'aide du paramtre nargin qui indique le nombre d'arguments en entre lors de l'appel de la fonction. Voici un exemple de calcul d'un produit scalaire ou d'une norme au carre, illustrant son utilisation : function r = psnorm2(a,b) if (nargin==1) r=a'*a elseif (nargin==2) r=a'*b end (a,b sont des vecteurs colonne)

psnorm2(u,v) renvoie le produit scalaire et psnorm2(u) renvoie la norme au carre de u. Remarquer que psnorm2(u,u) renvoit galement la norme au carre de u. Bien qu'il soit parfois dangereux d'utiliser cette possibilit, il est important de la connatre car de nombreuses fonctions natives de Matlab l'utilisent.

Porte d'une fonction

On peut dclarer des fonctions dans un script principal la suite des instructions du script principal. Mais il est souvent plus recommand, pour des raisons d'organisation, de les placer dans un autre script, voire plusieurs. Dans ce cas, le nom du script externe devra porter le nom de la fonction avec l'extension .m. Si dans ce script externe, on y place plusieurs fonctions, seule la fonction dont le script porte le nom sera accessible depuis le script d'appel. Les autres fonctions ne seront accessibles que dans le script o elles ont t dclares; on parle, dans ce cas, de sous-fonctions.

Egalement lie aux fonctions, la notion de variable globale permet de rendre visible des variables d'un script l'autre. Par dfaut, aucune variable n'est globale. Supposons que l'on ait un script principal et un script secondaire contenant la dclaration de la fonction nomfonction. Pour dclarer de faon globale, la variable de nom varg, on introduira l'instruction suivante dans le script principal et la fonction nomfonction : global varg

Fonction comme argument d'une fonction

Pour clore cette partie sur les fonctions, mentionnons l'existence de commandes permettant d'valuer une chane de caractres comme une commande Matlab. La command eval s'utilise ainsi :
comd='s1=sin(1)' comd = s1=sin(1) eval(comd) s1 = 0.8415

La commande feval permet quant elle de transfrer, en temps que chane de caractres, le nom d'une fonction que l'on veut valuer :
fon='sin' fon = sin s1=feval(fon,1) s1 = 0.8415

(derrire le nom de la fonction valuer, apparaissent les arguments)

6. Lecture et criture
Comme il l'a dj t mentionn, par dfaut toute commande excute produit un rsultat qui apparat dans la fentre d'excution la suite de la commande. On peut empcher l'affichage du rsultat en terminant la commande par ';'. Ainsi, on a :
H=[1 2 3] H = 1 2 H=[1 2 3];

(Rsultat affich) 3 (Aucun rsultat affich)

10

Dans les scripts d'excution, l'affichage d'un rsultat est exceptionnel. Il sert, essentiellement, dtecter les erreurs et afficher des rsultats finaux. On prendra donc garde ne pas oublier le ';' la fin de chaque ligne de commande.

Impression l'cran
La commande standard d'criture dans la fentre d'excution est fprintf qui a la structure gnrale suivante : fprintf(format,var1,var2,...) o format est une chane de caractres dcrivant le format d'criture des variables var1,var2,... que l'on souhaite afficher. Les principaux types de formats d'criture sont : %d %f %e %g %s entier rel exponentiel rel double prcision chane de caractres %5d : entier de longueur 5, par exemple 34562 %5.2f : rel de longueur 5 avec 2 chiffres aprs la virgule, par exemple 32.42 %10.8e : nombre de la forme -21.01e+05 mode automatique de dtection entre %e et %f

Par souci de simplicit, on peut se contenter d'utiliser les formats %d,%f,%e sans spcifier de longueur prcise. Par ailleurs on dispose de certains oprateurs de mise en forme , par exemple \n pour passer la ligne. Ainsi, on crira, par exemple : fprintf('\n Convergence en %d iterations ',it) o it dsigne une variable contenant un entier. Pour plus de dtails, faire help fprintf.

Impression dans un fichier

Il est galement possible d'crire les rsultats dans un fichier (et souhaitable lorsqu'il y en a beaucoup). Pour ce faire, on utilise encore la commande fprintf, mais en spcifiant un numro nfic associ un nom de fichier de rsultats, nomm ici ficres. On effectue, les oprations suivantes : nfic=fopen(ficres,'rw'); (ouvre le fichier ficres en mode lecture et criture) fprintf(fid,'\n Convergence en %d iterations ',it); (crit dans le fichier ficres) status=fclose(fid); (ferme le fichier ficres) L'opration d'ouverture de fichier par la commande fopen a chou si nfic vaut -1 et status renvoie 0 si l'opration de fermeture par la commande fclose est russie et -1 sinon.

Lecture de donnes
Afin de lire des donnes utiles l'excution, on peut procder de deux faons : soit en interrompant l'excution du programme et en demandant l'utilisateur d'indiquer les donnes, soit en lisant un fichier de donnes. Cette deuxime solution tant bien souvent prfrable. Pour interrompre l'excution et demander une valeur, on utilise la commande input, dont voici un exemple d'utilisation : data=input('Donnez votre valeur (par defaut 0)');

11

A l'issue de la rponse, la variable data contiendra la rponse envoye, qui peut prendre n'importe quel type et mme prendre la valeur vide [] si on a taper sur enter. C'est d'ailleurs par ce moyen que l'on gre les valeurs par dfaut : if (data ==[]) data=0; end Pour lire des fichiers de donne, on utilise la commande fscanf dont le principe de fonctionnement est voisin de la commande fprint. Pour lire la donne data dans le fichier ficdon, on utilisera la suite de commande suivante : nfic=fopen(ficdon,'r'); data=fscanf(fid,'%d'); status=fclose(fid); (ouvre le fichier ficdon en mode lecture) (lit dans le fichier ficdon) (ferme le fichier ficdon)

Pour plus d'informations, faire help fscanf, help fopen et help fclose. Il existe d'autres mthodes de lecture et d'criture sur fichier, faire help iofun pour de plus amples informations.

7. Graphique sous Matlab

Afin de terminer cette brve introduction Matlab, indiquons quelques fonctionnalits graphiques de Matlab. Donnons deux exemples : le trac d'une courbe et le trac d'un champ d'isovaleurs. Supposons que l'on veuille reprsenter graphiquement la courbe y=sin(x) sur l'intervalle [-,] avec 200 points. On excute alors les commandes suivantes :
x=-pi:pi/100:pi; y=sin(x); plot(x,y);

1 0.8 0.6 0.4 0.2 0 -0.2 -0.4 -0.6 -0.8 -1 -4

-3

-2

-1

Il existe de nombreuses options pour contrler l'affichage des courbes (couleur, axe, commentaire ...), faire help plot et plus gnralement help graph2d.

12

Maintenant, supposons que l'on veuille reprsenter la surface dfinie par la fonction z=sin(xy) sur le carr suivant [-/2,/2]x[-/2,/2] l'aide d'une grille de points 30x30. On utilise la squence de commandes suivantes :
[xi,yi]=meshgrid(-pi/2:pi/30:pi/2); zi=sin(xi.*yi); surf(xi,yi,zi,zi)

Pour les autres possibilits faire help graph3d.

13