Accurate estimation of complex eigenvalues

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
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