Discussion :

[jery190]

Bonsoir
quelqu'un(e) m'aide à compléter le programme de la résolution de AX=b par la méthode itérative richardson(gradient conjugué) pour une matrice quelconque inversible, mais je n'arrive pas comment placer alpha dans le prrogramme.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
A
b
M
a=det(M)
X0=ones(1)
itermax=100;
tol=0.0001;
R0=b-(A*X0)
if (a==0)
    la matrice  nest pas inversible"
else   a~=0
    for k=1:itermax
        while(norm(X-X0))/(norm(X0))<=0.0001
            R=M*Z
            X=X+alpha*Z
            R=R-alpha*A*Z
        end
        X
        R
    end
end

le programme s'arréte et fait l'error

Code :

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

1
2
3
4
5
6
 
??? Error using ==> unknown
Matrix dimensions must agree.
 
Error in ==> richardson at 8
R0=b-(A*X0)

Je veux SVP que qlq donne aussi le programme de gradient conjugué pour le cas dynamique???
Merci d'avance

[Jerome Briot]

Le message d'erreur t'indique que tu fais des opérations entre des matrices dont les dimensions ne concordent pas.

Vérifie les dimensions des variables b, A et X0 à l'aide de la fonction SIZE ou de la fenêtre Variable Editor

Je veux SVP que qlq donne aussi le programme de gradient conjugué pour le cas dynamique???

Nous ne sommes pas la pour faire le travail à ta place
Tu dois donc nous montrer ce que tu as commencé à coder... même si c'est faux.