eig return complex values
조회 수: 17 (최근 30일)
이전 댓글 표시
Hello,
I'm trying to find the eigenvalues and eigenvectors of an invertible matrix. The eig function returns me complex values.
But the matrix is invertible: I invert it on Pascal.
How to explain and especially how to solve this problem please?
The matrix I am trying to invert is the inv(C)*A matrix, from the attached files.
Thanks,
Michael
댓글 수: 5
채택된 답변
Torsten
2022년 1월 22일
편집: Torsten
2022년 1월 22일
Use
E = eig(A,C)
instead of
E = eig(inv(C)*A)
or
E = eig(C\A)
댓글 수: 4
Torsten
2022년 1월 22일
편집: Torsten
2022년 1월 22일
Although negligible, eig(A,C) produces no imaginary parts.
E = eig(A,C) solves for the lambda-values that satisfy
A*x = lambda*C*x (*)
for a vector x~=0.
If C is invertible, these are the eigenvalues of inv(C)*A (as you can see by multiplying (*) with
inv(C) ).
추가 답변 (1개)
Matt J
2022년 1월 22일
편집: Matt J
2022년 1월 22일
It turns out that B=C\A does have real eigenvalues in this particular case, but floating point errors approximations produce a small imaginary part that can be ignored.
load matrices
E=eig(C\A);
I=norm(imag(E))/norm(real(E))
So just discard the imaginary values,
E=real(E);
댓글 수: 2
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!