Discussion :

[Boule de coco]

Bonjour,

je dois écrire l'algorithme du gradient conjugué, est ce que quelqu'un pourrait me dire si mon programme est juste ?

le voici :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
function [x, k] = Gradient_conjugue(A, b, x, e)
k=0;
e=0;
x=x0;
r0=b-A*x0;
r=r0;
p0=r0;
while norm(r)<e*norm(r0)
    f=dot(r,r)/(A*p,p);
    x=x+f*p;
    r=r-f*A*p;
    B=(dot(r,r))^2/(dot(r,r))^2;
    p=r+B*p;
end

ou x^0 appartient à R^n, arbitraire, k=0,1,...
p^0=r^0=b-A*x^0
f^k=(r^k,r^k)/(A*p^k,p^k)
x^(k+1)=x^k+f^k*p^k
r^(k+1)=r^k-f^k*A*p^k
B^(k+1)= llr^(k+1)ll^2 / llr^kll^2
p^(k+1) = r^(k+1) + B^(k+1)*p^k

K(e) est le nb d'itération pour que llr^kll soit inférieur à epsilon(noté e) fois llr^0ll.


Merci bcp et bonne journée

[Jerome Briot]

Pour faire simple :

• ceci n'est pas un algorithme
• la variable p n'est définie (initialisée ?) nulle part
• la syntaxe (A*p,p) est illégale

[mr_samurai]

• la syntaxe (A*p,p) est illégale

dot(A*p,p) à la place , mais je ne sais pas si ton code fait bien ce que tu veux.

[mr_samurai]

Tu ne va jamais passer par ta boucle while :

Code :

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

1
2
3
e=0;
%...
while norm(r)<e*norm(r0))

++

[Boule de coco]

Vous voulez dire qu'en mettant e=0; ca fait que ma boucle while ne fonctionne pas, c'est ca?

[mr_samurai]

euuuh,

tu comprend au moins le code que t'as posté ?

e est un argument de ta fonction.... Si tu veux lui donner une valeur par default, il faut utiliser nargin.


PS :
pour répondre à ta question :

Vous voulez dire qu'en mettant e=0; ca fait que ma boucle while ne fonctionne pas, c'est ca?

Code :

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

while norm(r)<0

un norme est par definition >= 0 . donc la condition de passage est fausse.