Resolution des systemes lineaire et skyline

maestro -  
khunlun Messages postés 123 Statut Membre -
Bonjour les ami(e)s,
je voudrais un code matlab qui :
- Stoke une matrice, definie positive et symetrique, sous forme skyline
- Fait la decomposition de la matrice, A=LDL^t (decomposition de cholesky), lorsque A est stocke sous la forme skyline
- Resoud le systeme Ax=B pour un vecteur B donne en utilisant la decomposition de A ci-dessus.
merci d'avance

1 réponse

khunlun Messages postés 123 Statut Membre 27
 
Bonjour,

Je ne connais pas skyline.
Mais je connais les commandes suivantes qui pourrait aider.

- chol(M) pour la factorisation de Cholesky.
- Ax = b peut se résoudre en utilisant : x = A\b.

Bonne chance !
0
maestro
 
merci khunlun,
en fait jai une matrice creuse et de grande taille et skyline est un mode de stockage ki consiste qu'a ne stocker que les termes non nuls de la matrice (profil de la matrice).
la difficulte consite donc a manipuler le profil de la matrice en lieu et place de la matrice.
0
khunlun Messages postés 123 Statut Membre 27
 
Re bonjour,

Il y a la fonction sparse qui permet de stocker des matrices en ne spécifiant que la position est la valeur des éléments non nul.

Exemple de l'aide :
sparse(1:n,1:n,1) donne la matrice identité.

J'espère que ca aidera un peu :p

Bonne chance
0