Discussion :

[Sevrynn]

Bonjour,

Je dois réaliser sur matlab un programme permettant de résoudre Ax=b (avec A matrice et b vecteur) par la méthode de Gauss Jordan. J'ai réussi à mettre en place les 3/4 du programme seulement je n'arrive pas à mettre en place la partie concernant la recherche de la ligne de A contenant le pivot maximal afin de réaliser ensuite une inversion avec la première ligne. Je ne sais pas si je suis claire. En tout cas ce petit problème sur Matlab m'ennuie légèrement...

Merci de votre aide..

[rostomus]

Bonjour,

je te propose ça:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
d=diag(x);
b=(d==max(d));
h=x(b,:);
x(b,:)=x(1,:);
x(1,:)=h;

[Sevrynn]

Merci beaucoup... et si je veux remplacer la diagonale par la première colonne, il faudrait que je le rentre ainsi ?

p=x(:,1);
b=(p==max(p));
h=x(b,:);
x(b,:)=x(1,:);
x(1,:)=h;

[rostomus]

remplacer la diagonale par la premiere colonne:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
d=eye(size(x));
d=(d==1);
h=x(d);
x(d)=x(:,1);
x(:,1)=h;

[Sevrynn]

Merci beaucoup de votre aide

[Jerome Briot]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
d=eye(size(x));
d=(d==1);
h=x(d);

Bonjour,

ceci peut se simplifier avec la fonction DIAG :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

h=diag(x);

Pour le problème posé ici, on peut aussi jouer avec l'indexage (si la matrice est carrée) :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
x=rand(5)
 
nrows=size(x,1);
idx=[2:nrows nrows+2:nrows+1:numel(x)];
x(idx)=x([idx(nrows:end) idx(1:nrows-1)])

Il est possible de généraliser cette méthode pour remplacer n'importe quelle ligne/colonne par la diagonale.