Re: generalized eigenvalue problem



You mentioned earlier that the matrix A is singular. This does not
agree with your subsequent assertion that it is SPD. I'm not sure
what the inequality is that you put. This looks like something to do
with diagonal dominance.

Are your matrices full? If they are than you can try converting it to
a regular EVP. This is not a good idea if they are sparse, since the
resulting EVP will likely be full.

Infinite eigenvalues can happen:

A=[1,1;1,2]
B=[1,-1;-1,1]

eig(A,B) returns the eigenvalues 1/5 and Inf. I do know that it is
possible to get complex eigenvalues from a symmetric pencil though.

If you say something more about A,B and where they come from, I might
be able to add something more useful.

Edwin Ramayya wrote:


Dear Greg,
Thankyou so much for ur quick response. I did try what you
suggested!
but still I get some junk eigenvalues which are complex.

I checked for (SPD) Aii*Ajj > |Aij|, i ~= j, and checked for
Hermitian condition too.

A, B satisfy the condition!

anyother idea?

thanks,
Edwin.
.



Relevant Pages

  • inverse iteration for complex eigenvalues
    ... i have been using the 'inverse algorithm' to find eigenvalues of a real ... but it seems to only converge satisfactorily for real eigenvalues. ... for complex eigenvalues it seems to give spurious and widely converging ...
    (sci.math.num-analysis)
  • Re: Signs of eigenvalues after a positve row transformation
    ... >>1) the number of complex eigenvalues of the new matrix is greater than or equal ... >>3) the number of negative eigenvalues of the new matrix is not greater than the ... >It may depend on what you mean by "perturb". ...
    (sci.math)
  • Re: Signs of eigenvalues after a positve row transformation
    ... >1) the number of complex eigenvalues of the new matrix is greater than or equal ... >3) the number of negative eigenvalues of the new matrix is not greater than the ... It may depend on what you mean by "perturb". ...
    (sci.math)
  • Re: QR algorithm for eigenvalues
    ... I don't remember this but I think there are some modifications of the QR algorithm for computing complex eigenvalues. ... See e.g. Golub: Matrix Computations. ...
    (sci.math.num-analysis)