Non-normal matrix eigen-values inaccuracy
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
I am interested in the eigenvalues and eigenvectors of a non-normal (complex non-hermitian) matrix. I find that when I use eig the eigenvectors that were computed are not accurate i.e. when I plug them back into the matrix equation they do not satisfy the eigenvalue condition to great accuracy. The eigenvalues are all distinct and hence the matrix should be diagonalizable. The form of the matrix I am interested in is: M = [a,0,b,b;0,-conj(a),-conj(b),-conj(b);conj(b),b,4,0;-conj(b),-b,0,-4] with a = (1.3800 - 0.0364i)*10^4; and b = 0.0127 - 0.4820i for example. [v,d] = eig(M); u = M*v(:,1)-d(1,1)*ones(4,1); the vector u should ideally be very close to zero but it is not!
I think this may be a bug with Lapack but I wanted to first ask in the Matlab community before going to Lapack fora. Cheers Prasanna
0 Kommentare
Antworten (0)
Siehe auch
Kategorien
Mehr zu Linear Algebra finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!