1 pièce(s) jointe(s)
compléter le programme d'une méthode itérative
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:
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:
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
