L'éditeur Matlab est l'interface qui permet d'écrire les fonctions. Il peut être lancé en tapant "edit NomDeFichier" dans le command window ou en cliquant sur File>>New>>M-file.


Vous pouvez réaliser soit des scripts soit des fonctions:

• Les scripts sont des programmes sans argument dont le contenu s'exécute comme un copier-coller dans le command window. Les scripts ont accès aux variables stockées dans le workspace. Aucune entête n'est nécessaire.
  Pour l'exécuter, il suffit de taper le nom du fichier dans le command window
• Les fonctions sont des programmes qui reçoivent des arguments en entrée et qui renvoie un ou plusieurs résultats en sortie. Pour réaliser une fonction, il faut que la première ligne du fichier soit :
  function [output1,output2, ...]=NomDeFonction(Input1,Input2,...)
  Le fichier contenant la fonction doit avoir le même nom que celle-ci. Il est cependant possible d'insérer plusieurs fonctions dans un fichier mais celles-ci ne seront accessible qu'à l'intérieur même du fichier.

Les commentaires insérer dans les premières lignes (Signalés avec un "%") seront afficher lorsque vous taperez "help NomDeFichier". Voici un exemple de fonction réalisant la somme et la différence de 2 nombres :

function [s,d]=sumdif(a,b) % Effectue la somme et la différence de deux nombres % % Syntaxe : % [s,d]=sumdif(a,b); % % Sorties : % s somme % d différence s=a+b; d=a-b;

Un mauvais programmeur a codé la fonction suivante, permettant de résoudre un système d’équation de type Ax=b par la méthode de Gauss (Rappel ici). Le bon programmeur aurait vérifié chacune des étapes au fur et à mesure mais celui là a écrit le programme sans réfléchir et ne comprend pas pourquoi ça ne marche pas.
En utilisant le debugger vous devez corriger toutes les erreurs de ce programme. (On ne s’occupera pas des cas particuliers où A(k,k)=0).

function x=gaussres(A,b) n=size(A); for k=1:n for i=k:n A(i,:)=A(i,:)-A(i,k)/A(k,k)*A(k,:); b(i)=b(i)-A(i,k)/A(k,k)*b(k); end end x=zeros(n,1); for i=1:n x(i)=(b(i)-sum(A([i,i+1:end])*x(i:end)))/A(i,i); end

Un petit exemple de résolution que vous devez obtenir :


Correction : (Sélectionner pour voir la solution)

function x=gaussres(A,b) n=size(A,1); for k=1:n for i=k+1:n c=A(i,k)/A(k,k); A(i,:)=A(i,:)-c*A(k,:); b(i)=b(i)-c*b(k); end end x=zeros(n,1); for i=n:-1:1 x(i)=(b(i)-sum(A(i,i+1:end)*x(i+1:end)))/A(i,i); end