Discussion :

[NELLLY]

Salut, j'ai à chercher les valeurs propres et les vecteurs propres de la matrice suivante:

Code :

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

1
2
3
4
5
6
 
P=[0 0 0 1
7/8 3/4 0 0
1/16 1/8 1/2 0
1/16 1/8 1/2 0];
>> [v,D]=eig(P);

Jusqu'ici tout va bien. J'ai besoin ensuite de chercher l'élément de la diagonale qui contient la valeur 1. Pour celà j'ai utilisé la commande suivante:

Code :

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

1
2
 
>> find(diag(D)==1)

Mais au lieu d'avoir comme réponse "1" j'ai la réponse "empty-matrix". Je crois que c'est dû au fait qu'il y a des valeurs imaginaires dans la diagonale D. Comment faire?
Merci

[Jerome Briot]

Je crois que c'est dû au fait qu'il y a des valeurs imaginaires dans la diagonale D.

Non.

L'explication :

Code :

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

1
2
3
4
5
6
7
8
>> sprintf('%.18f\n',diag(D))
 
ans =
 
1.000000000000000400
0.124999999999999930
0.124999999999999930
0.000000000000000164

Comment faire?

En consultant la FAQ : Pourquoi 0.3-0.2-0.1 est-il différent de 0 ?

[NELLLY]

Salut
j'ai essayé l'instruction suivante

Code :

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

1
2
3
 
d=diag(D)
find(abs(d-1)<eps)

J'ai obtenu aussi la réponse"empty matrix"

[Jerome Briot]

Il ne faut pas toujours utiliser eps mais un ordre de grandeur propre à ton problème.

Par exemple :

Code :

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

find(abs(d-1)<1E-3)

ou

Code :

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

find(abs(d-1)<1E-6)

Nous allons corriger la FAQ.

[NELLLY]

Merci Dut. ça fonctionne à merveille maintenant.