Accurate estimation of complex eigenvalues

Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
1 reply · 2K views
soikez
Messages
5
Reaction score
0
Hello,
I use Arnoldi iterative algorithm in order to compute the eigenvalues of a matrix. I know that the eigenvalues are of the form [tex]\lambda(1+j/c)[/tex] and I can totally estimate them. The problem that occurs is that both the range of [tex]\lambda_0[/tex] and c is for example [100,1000]. That means that there is a significant difference in the order of real and imaginary part e.g. [tex](10^4,10^6)[/tex], so the algorithm I use defines with more accuracy the real part of the eigenvalue. To be more specific let me give an example. If I estimate an eigenvalue, using python or matlab, as [tex]\lambda=100+0.017j[/tex], the analytically computed eigenvalue is [tex]\lambda=100.1+0.012j[/tex]. This deviation in the imaginary part causes a lot of problems.
Is there any way of normalization in order to exceed this problem of accuracy?

Thanks in advance!
 
Physics news on Phys.org
You could try to describe it as a purely real problem by using the vector space isomorphism ##\mathbb{C}\cong \mathbb{R}\oplus \mathbb{R}##.