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 'Determine whether ##125## is a unit in ##\mathbb{Z_471}##'

Thread 'Determine whether ##125## is a unit in ##\mathbb{Z_471}##'

This is the question, I understand the concept, in ##\mathbb{Z_n}## an element is a is a unit if and only if gcd( a,n) =1. My understanding of backwards substitution, ... i have using Euclidean algorithm, ##471 = 3⋅121 + 108## ##121 = 1⋅108 + 13## ##108 =8⋅13+4## ##13=3⋅4+1## ##4=4⋅1+0## using back-substitution, ##1=13-3⋅4## ##=(121-1⋅108)-3(108-8⋅13)## ... ##= 121-(471-3⋅121)-3⋅471+9⋅121+24⋅121-24(471-3⋅121## ##=121-471+3⋅121-3⋅471+9⋅121+24⋅121-24⋅471+72⋅121##...
View full post »

Thread 'About the existence of Hamel basis for vector spaces'

It is well known that a vector space always admits an algebraic (Hamel) basis. This is a theorem that follows from Zorn's lemma based on the Axiom of Choice (AC). Now consider any specific instance of vector space. Since the AC axiom may or may not be included in the underlying set theory, might there be examples of vector spaces in which an Hamel basis actually doesn't exist ?
View full post »

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 »

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
694
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