je cherche un algo ou un code source qui prend en entrée une matrices (de double) et qui donne TOUTES les valeurs propres ainsi que les vecteurs propresassociés à cette matrice.
non, j'ai que dalle
je suis laché en plein désert ....
y'a apparement la méthode des puissances ittérées (que je ne connais pas mais qui donne la plus grande val propre et son vect propre associé)
pour les autres ....
ciao
moi j'ai une méthode que j'ai apprise en cours (pas celles des puissances itérées), mais elle doit être difficile à implémenter (les valeurs propres de la matrice M sont les racines du polynome det(M-XId))
le déterminant ce n'est pas difficile, il existe des méthodes systématiques. Ce qui pose problème c'est comment représenter le 'X' (problème de programmation), et surtout résoudre le polynome final (la c'est un problème plus général)
bah c'est un peu chiant à expliquer comme ca.
déja je peux te dire que je l'ai fait dans une de mes sources (). En fait j'utilise ce qu'on appelle le développement par rapport à une ligne ou à un colonne:
*si la matrice est de taille 2*2
a|b
c|d
le déterminant vaut a*d-c*b
*sinon tu choisis un ligne
pour chaque élément de cette ligne tu calcule le déterminant de la matrice obtenue en enlevant la ligne et la colonne de l'élément et tu multiplie par -1 une fois sur deux (celle d'en haut à gauche vaut 1).
1|-1|1...
-1|1|-1...
...
le déterminant de la matrice est le somme de tous ces déterminants
voila j'avais prévenue, je veux bien croire que ce n'est pas très clair
La méthode réucrsive des cofacteurs est la pire des méthodes qui soit pour calculer un déterminant, vu que c'est qq chose comme du N^3, donc au-delà d'une matrice 10x10, c'est prohibitif.
Pour ça, il vaut mieux trigonaliser la matrice, et le déterminant est alors le produit des éléments diagonaux.
Pour trigonaliser, on peut utiliser la méthode de Gauss, Gauss-Seidel par exemple (chercher Gauss sur le site).
moi aussi ça m'intéresse ! Et les bibliothèques genre CLAPACK, tnt, jama, gmp, ça suxx, alors si vous avez un truc intéressant pour calculer les valeurs propres et les vecteurs propres, mettez un lien svp
--------------------------------------------------------------------
"Et là, ça compile... Est-ce que ça fait s'qu'on veut ?..."