How do I estimate complex eigenvalues?

photis
Messages
7
Reaction score
0
Let A be a matrix with real elements. The problem is to estimate eigenvalues of A, real and complex. QR algorithm is fine for real eigenvalues, but obviously fails to converge on complex eigenvalues... So, I'm looking for an alternative that could provide an estimate for complex eigenvalues of A. Can anybody help?

Thanks,
photis
 
Physics news on Phys.org
Go to http://www.netlib.org/lapack/double/. The routine dgeev.f (a fortran subprogramme) is good and does not require complex arithmetic because you have a real matrix, so that all complex eigenvalues (if there are any) come in pairs of complex conjugates. I have used this routine successfully many many times in the past but I don't remember the details of the numerical method. Usually, the matrix has to be reduced to Hessenberg form first; then, I think there may be a variant of the QR algorithm that works. I found the discussions in "Numerical Recipes" by Press et al. and "Numerical methods that work" by F.S. Acton very imformative. Good luck---the world of non-symmetric matrices is not a pleasant place!

For more about the numerical procedure, look in
http://www.netlib.org/lapack/lug/node50.html
 
Last edited:
Thanks, your help was really valueable:!) . Following the links, I found this: http://www.acm.caltech.edu/~mlatini/research/qr_alg-feb04.pdf"

As eigenvalues come in conjugate pairs, QR apparently fails (no dominant eigenvalue exists). However, instead of generating a single eigenvalue estimation, QR produces a 2x2 matrix "containing" the conjugate pair. One can either (a) calculate the eigenvalues of the 2x2 matrix directly and proceed with next eigenvalue(s) or (b) use "Wilkinson shift" to move QR on the complex plain.

(a) may affect estimates of the remaining eigenvalues, but (b) introduces complex arithmetic, so (a) seems preferable. After all, I use QR to get reasonable initial approximations for inverse power algorithm.

Does anybody know if there is an alternative method (not QR) to ger Schur quasitriangular form of a real matrix?
 
Last edited by a moderator:
i have a matrix that its elements are complex ,ineed its eignvalue , did any onehave any subroutine for that?

(not double percision)
tnx
 

Thread 'Trouble understanding an online solution to an exercise in Dummit & Foote'

##\textbf{Exercise 10}:## I came across the following solution online: Questions: 1. When the author states in "that ring (not sure if he is referring to ##R## or ##R/\mathfrak{p}##, but I am guessing the later) ##x_n x_{n+1}=0## for all odd $n$ and ##x_{n+1}## is invertible, so that ##x_n=0##" 2. How does ##x_nx_{n+1}=0## implies that ##x_{n+1}## is invertible and ##x_n=0##. I mean if the quotient ring ##R/\mathfrak{p}## is an integral domain, and ##x_{n+1}## is invertible then...
View full post »

Thread 'Questions about non existence of GCDs vs (coimages, cokernels)'

The following are taken from the two sources, 1) from this online page and the book An Introduction to Module Theory by: Ibrahim Assem, Flavio U. Coelho. In the Abelian Categories chapter in the module theory text on page 157, right after presenting IV.2.21 Definition, the authors states "Image and coimage may or may not exist, but if they do, then they are unique up to isomorphism (because so are kernels and cokernels). Also in the reference url page above, the authors present two...
View full post »

Thread 'Decomposition into irreps of compact Lie group'

When decomposing a representation ##\rho## of a finite group ##G## into irreducible representations, we can find the number of times the representation contains a particular irrep ##\rho_0## through the character inner product $$ \langle \chi, \chi_0\rangle = \frac{1}{|G|} \sum_{g\in G} \chi(g) \chi_0(g)^*$$ where ##\chi## and ##\chi_0## are the characters of ##\rho## and ##\rho_0##, respectively. Since all group elements in the same conjugacy class have the same characters, this may be...
View full post »

Similar threads

I Getting eigenvalues of an arbitrary matrix with programming
Replies
1
Views
2K
I Can one find a matrix that's 'unique' to a collection of eigenvectors?
Replies
33
Views
676
A Eigenvalues of block matrix/Related non-linear eigenvalue problem
Replies
1
Views
2K
MHB Eigenvalues are real numbers and satisfy inequality
Replies
1
Views
2K
I Types of complex matrices, why only 3?
Replies
5
Views
2K
Finding the eigenvalues of a complex matrix
Replies
1
Views
3K
MHB Eigenvector of 3x3 matrix with complex eigenvalues
Replies
1
Views
4K
I Quintic Polynomial Tschirnhaus Transform: Possible complex Bring Radicals
Replies
0
Views
2K
F=ℝ: Normal matrix with real eigenvalues but not diagonalizable
Replies
1
Views
2K
I Question about an eigenvalue problem: range space
Replies
1
Views
2K

Hot Threads

Recent Insights

Back
Top