complex eigenvalues
显示 更早的评论
Hi, I have a square symmetric matrix (5,5) with complex entries,the output eigenvalues when I use eig(T) are all complex .I want to determine the smallest negative eigenvalue.I don't know how ,any one can help.
2 个评论
Sean de Wolski
2011-12-5
smallest as in the closest to infinity or closest to zero?
zayed
2011-12-5
采纳的回答
Walter Roberson
2011-12-5
"smallest" is not defined for complex numbers. "negative" is not defined for complex numbers either.
You can compare real parts, or you can compare imaginary parts, or you can compare magnitudes.
[vals, idx] = min(real(E));
E(idx)
or
[vals, idx] = min(imag(E));
E(idx)
or
[vals, idx] = min(abs(E));
E(idx)
8 个评论
zayed
2011-12-5
Walter Roberson
2011-12-6
Euler's Theorem http://en.wikipedia.org/wiki/Euler%27s_theorem is completely irrelevant to this question.
MATLAB stores complex numbers as a pair of arrays, one for the real part and one for the imaginary part. real() and imag() copy the appropriate part _exactly_ with no loss of precision at all.
If you are concerned that abs() might perhaps not have sufficient precision, then you could instead use hypot(real(E), imag(E)) as hypot takes special care to avoid overflow and underflow; see http://www.mathworks.com/help/techdoc/ref/hypot.html
zayed
2011-12-6
Walter Roberson
2011-12-6
You don't get how to use _what_?
I already posted complete code for selecting the eigenvalue with minimum real component, with minimum complex component, or with minimum magnitude.
If you are asking how to use hypot as a replacement for abs() in order to take in to account possible overflow or underflow (which is not very likely in your situation, but just-in-case):
[vals, idx] = min(hypot(real(E), imag(E)));
E(idx)
which just replaces abs(E) with the hypot() code I showed above.
In previous threads, for your purposes you defined the minimum eigenvalue as being the one with the most negative real component. Your code to implement that looked acceptable -- not as efficient as what I wrote above, but workable. Was there a problem when you used that code?
zayed
2011-12-6
zayed
2011-12-6
Walter Roberson
2011-12-6
hypot has been in versions since sometime in 2008 or before; I have not traced it further.
When I read Loren's blog about hypot, I see in the comments that abs() is also implemented robustly, so there would be no advantage to using hypot() over using abs(), so you might as well not bother.
http://blogs.mathworks.com/loren/2008/02/07/why-hypot/
If you are looking for the eigenvalue with the smallest magnitude (such as min(abs(E)) would find), then you could instead use
gamma = eigs(T,1,'sm');
which will find just the one eigenvalue. Smallest magnitude could be positive or negative for the real or imaginary components, though -- the eigenvalue closest to 0. There is unfortunately no way with eigs to pick out just the complex eigenvalue with the real component or imaginary component closest to negative infinity: you will have to use one of the above min() forms for that.
zayed
2011-12-6
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Linear Algebra 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
