Discussion :

[tpscience]

Bonjour à tous,

J'ai un petit problème avec la fonction

Code :

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

eig

j'ai l'impression.
Je demande à mon programme de me sortir les vecteurs et valeurs propres d'une matrice, ceux-ci pouvant être réels ou complexes, suivant différentes valeurs d'une variable définissant ma matrice. Or, pour mes dernières valeurs de cette variable, matlab m'inverse le classement de mes 2 valeurs propres. Par défaut je pense qu'il classe d'abord la solution "positive" (dans le sens où je prend le signe devant le discriminant dans mes solutions) puis la solution "négative". Mais il m'inverse les 2 au bout d'un moment...

J'ai tenté de forcer le classement avec la fonction

Code :

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

eigs(A,k,'lr')

mais cela n'y fait rien...?!

Auriez-vous une idée ?

Merci.

[le fab]

pas d'idée sinon que eigs est sensée être utilisée sur des matrices sparses
après moi je ne m'en suis servi que pour vérifier que toutes les valeurs propres étaient de même signe, je n'ai jamais fait attention à l'ordre

[tpscience]

En fait l'ordre m'intéresse puisque je trace le module et l'argument de chacune des 2, donc si matlab me les inverse à un moment donné, cela n'a plus de sens !
C'est d'ailleurs étrange tout ça, car à une position donnée dans ma matrice diagonale devrait correspondre toujours le même type de solution...

[FLB]

Salut,
Matlab n'assure pas d'ordre en sortie de eig, mais tu peux appliquer la méthode sort à tes valeurs propres - tant que les valeurs propres que tu vois évoluer ne se croisent pas.
La seconde sortie de sort te permettra aussi de trier tes vecteurs propres dans l'ordre trié des valeurs propres.

[tpscience]

Merci pour vos réponses, en fait ma matrice est diagonale, donc

Code :

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

sort

ne suffit pas, puisque

Code :

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

sort

classe soit par colonne, soit par ligne. J'ai fais en fait une succession de calcul d'extrema ajoutée à des boucles conditionnelles. Donc c'est un peu en force brute, mais ça marche, donc...

Merci à vous.

[FLB]

Pour récupérer la diagonale de D, et la trier :

Code :

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

sort(diag(D))

[tachmou]

Bonjour,

Je trouve cette question très intéressante. L'interprétation des valeurs propres est - de mon point de vue - loin d'être si évidente que ça. Tu peux peut-être trouver une piste de réponse en regardant du côté des vecteurs propres associés.

En effet, tes valeurs propres s'expriment de pair avec le vecteur propre qui lui est associé. Matlab, de son côté, quand il te fourni des valeurs propres et des vecteurs propres, il t'assure juste qu'ils sont fournis dans le même ordre - l'un par rapport à l'autre.

L'ordonnancement des valeurs propres seules, en soit, n'a pas vraiment de sens établi. D'autant que le calcul de ces valeurs est effectuée numériquement.

Maintenant, ce qu'il serait intéressant de savoir pour te répondre de manière un peu plus adéquate, c'est pour quelle application tu étudies ces structures ?

Personnellement, "l'inversion" de deux valeurs propres peut paraître surprenante, mais possède une interprétation claire avec les vecteurs propres. En particulier, si tu cherches à faire de l'identification avec la seule aide des valeurs propres... excuse moi mais - selon mon point de vue - tu es mal barré. Il faut que tu prennes en compte les vecteurs propres aussi. Mais tout dépend de l'application je pense, d'où ma question.

A suivre...

Cordialement,