matlab les fonctions

bonjour
voila j'ai écris deux programme
je calcul la solution d'un système de 4 équations avec 4 variables avec deux méthodes
1ere méthode celle de cramer dont le code

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19

function d=determinant(a,b,n) e=a; s=b; for i=1:4 e(i,n)=s(i); end d=det(e); end function v=cramer(a,b) if det(a)~=0 for i=1:4 m(i)=determinant(a,b,i)/det(a); end else erreur('division par zero'); end v=m; end

---

2ème méthode methode de gauss pivot dont le code

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

function v=methpiv(a,b) if a(1,1)=0 for i=2:4 if a(i,1)~=0 for j=1:4 va=a(i,j); a(i,j)=a(1,j); a(1,j)=va end i=5; end end end e=zeros(4); r=a; f=b; for i=1:4 e(1,i)=r(1,i); end for i=1:4 e(2,i)=r(2,i)-(r(2,1)/r(1,1))*r(1,i); e(3,i)=r(3,i)-(r(3,1)/r(1,1))*r(1,i); e(4,i)=r(4,i)-(r(4,1)/r(1,1))*r(1,i); end f(2)=f(2)-(r(2,1)/r(1,1)*f(1)); f(3)=f(3)-(r(3,1)/r(1,1)*f(1)); f(4)=f(4)-(r(4,1)/r(1,1)*f(1)); h=zeros(4); for i=1:2 for j=1:4 h(i,j)=e(i,j); end end for i=2:4 h(3,i)=e(3,i)-(e(3,2)/e(2,2))*e(2,i); h(4,i)=e(4,i)-(e(4,2)/e(2,2))*e(2,i); end g(1)=f(1); g(2)=f(2); g(3)=f(3)-(e(3,2)/e(2,2))*f(2); g(4)=f(4)-(e(4,2)/e(2,2))*f(2); s=zeros(4); for i=1:3 for j=1:4 s(i,j)=h(i,j); end end for i=3:4 s(4,i)=h(4,i)-(h(4,3)/h(3,3))*h(3,i); end p(1)=g(1); p(2)=g(2); p(3)=g(3); p(4)=g(4)-(h(4,3)/h(3,3))*g(3); v(4)=p(4)/s(4,4); v(3)=(p(3)-s(3,4)*v(4))/s(3,3); v(2)=(p(2)-s(2,3)*v(3)-s(2,4)*v(4))/s(2,2); v(1)=(p(1)-s(1,2)*v(2)-s(1,3)*v(3)-s(1,4)*v(4))/s(1,1); end

---

voila je vais ajouter le fait que l'utilisateur s'il entre 1 le système se résous par la méthode de cramer ,2 se résous par la méthode gauss
j'ai pensé a ajouter une variable dans les arguments de chaque fonction mais comment faire
je voudrais que vous me vérifier si le code en rouge qui vérifie si le pivot est nul et cherche le pivot non nul dans la premiere colonne puis échange avec cette ligne est juste ou non

merci d'avance

Bonjour,

Pour le choix, je ne comprends pas bien : tu as deux fonctions différentes, donc il suffit de faire un test sur l'entrée de l'utilisateur et d'appeler la bonne fonction?

Pour ton code en rouge, il est juste, mais tu peux l'améliorer :
tout d'abord, lorsque tu trouves un terme non nul, pour arrêter la boucle for il faut éviter de modifier la variable d'incrémentation. Pour cela, tu peux utiliser un break (cela permet de mettre fin à la boucle for) :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12

if a(1,1)=0 for i=2:4 if a(i,1)~=0 for j=1:4 va=a(i,j); a(i,j)=a(1,j); a(1,j)=va end break end end end

---

La recherche du premier terme non nul et la permutation des lignes peut aussi être simplifié :

Code:

---

1 2 3 4 5 6

if a(1,1)=0 i = find(a(:,1),1); if ~isempty(i) a([1 i],:) = a([i 1],:); end end

---

J'ai mis un test car je ne sais pas si le cas ou la première colonne a des éléments tous nuls peut arriver.

De la même façon, la plupart des boucles de ton code peuvent être supprimées en prenant en compte le fait que matlab fait très bien des opérations sur les tableaux.
Ex :

Code:

---

1 2 3

for i=1:4 e(i,n)=s(i); end

---

s'écrit :

Code:

---

e(:,n)=s(i);

---

bonsoir
est ce que je peux faire ça
pour la fonction utilisant la méthode cramer

Code:

---

function v=cramer(a,b)

---

pour la fonction utilisant la méthode pivot

Code:

---

function v=methpiv(a,b)

---

la fonction pour le choix de l'utilisateur

Code:

---

1 2 3 4 5 6 7 8 9

function x=systemeL(a,b,n) if n==1 e=cramer(a,b); end if n==2 e=methodepivote(a,b); end x=e; end

---

svp ou je dois mettre break
merci

Oui, par exemple, tu peux faire ainsi.

J'avais copié le mauvais code dans mon précédent post, je l'ai édité en rajoutant le break.

merci beaucoup magelan