bonjour
je croyais avoir déjà créer la discussion suivante mais je ne la retrouve pas
alors la voilà
La matrice A est carrée (taille n) et B vecteur de Rn. A inversible
méthode itérative x(0) dans Rn
x(k+1) =x(k) +d (b-A x(k) ) (1) pour k>=1
si la méthode converge la limite est l=A-1 b
soit Xd=I-dA et c= db permet décrire (1) sous la forme
x(k+1) = Xdx(k) +c
la mérhode converge si le rayon specttral de Xd <1
comme A est symétrique , Xd est aussi symétrique donc
son rayon spectral est égale à sa norme 2
la condition rayon specttral de Xd <1
équivaut à llI-dAll<1 soit dllAll>0 (c'est là où je ne suis pas sûre de mes calculs) or llAll>0 donc la condition est d>0
Enfin si A est symétrique définie positive
la condition sur d pour la convergence de la méthode devient
soit L valeur propre de A , on sait que L>0
soit F valeur propre de (I-dA)
(1-F)/d est valeur propre de A
réciproquement
pour L, on obtient -dL+1 est valeur propre de (I-dA)
pour la convergence , il faut que -1<-dL+1<1 et comme L>0
on obtient d doit appartenir à ]0, 2/L[
Ensuite, il faut trouver d* qui donne la convergence la plus rapide????
je ne sais pas faire, on dit qu'il faut que le rayon spectrale de
Xd*<=à celui de Xd
enfin il faut calculer le rayon spectral Xd*: comment faire?
merci par avance pour votre aide
fifrelette
-----