Register to reply

Eigen Values of Positive Definite Matrix

by vaibhavtewari
Tags: definite, eigen, matrix, positive, values
Share this thread:
vaibhavtewari
#1
Sep13-11, 03:05 PM
P: 67
Hi, i was wondering is there a way I can find the eigenvalues of a positive definite matrix. ?
Phys.Org News Partner Science news on Phys.org
Experts defend operational earthquake forecasting, counter critiques
EU urged to convert TV frequencies to mobile broadband
Sierra Nevada freshwater runoff could drop 26 percent by 2100
micromass
#2
Sep13-11, 03:32 PM
Mentor
micromass's Avatar
P: 18,333
Yes, find the characteristic polynomial and find its roots. The roots are the eigenvalues and should be positive.
vaibhavtewari
#3
Sep13-11, 04:03 PM
P: 67
That is true for any matrix, I was hoping to get an algorithm ..my matrix is large

micromass
#4
Sep13-11, 04:07 PM
Mentor
micromass's Avatar
P: 18,333
Eigen Values of Positive Definite Matrix

Google uses the power iteration (http://en.wikipedia.org/wiki/Power_iteration ). Maybe that'll help?
AlephZero
#5
Sep13-11, 05:15 PM
Engineering
Sci Advisor
HW Helper
Thanks
P: 7,177
How large is "large?" 100x100 or 1 millon x 1 million? Do you want all the eigenvalues, or just a few of them (either the largest or smallest ones?)

There are several different methods that will find all the eigenvalues of a moderate sized matrix (up to about 500 x 500 or 1000 x 1000 if you get lucky) efficiently. Google for the EISPACK library for computer codes. Trying to do this by hand is far too much work for anything bigger than about a 3x3 matrix.

For finding some eigenvalues of bigger postive definite matrices, probably the best method is the Lanczos algorithm, but I woudn't advise trying to write your own code to implement it (even though the maths appears to be quite straightforward) because there are some subtle numerical issues that you have to get right, or else it will fail to converge.
vaibhavtewari
#6
Sep13-11, 05:45 PM
P: 67
What are your thoughts regarding QR decomposition ?
vaibhavtewari
#7
Sep13-11, 05:46 PM
P: 67
my maximum matrix size is 100*100, I want to find all the eigen values.
DuncanM
#8
Sep13-11, 11:25 PM
P: 90
Just use one of the many routines or algorithms available. For example, MATLAB, Maple, and Mathematica are commercial software packages that could do it. Free packages are also available (for example, Octave is a free MATLAB clone). Or if you want an algorithm, use one of the many ones available on the NETLIB site. You could probably find one specific to your problem (for example, eigenvalues only, real entries, symmetric, etc.).
AlephZero
#9
Sep14-11, 07:13 PM
Engineering
Sci Advisor
HW Helper
Thanks
P: 7,177
Quote Quote by vaibhavtewari View Post
What are your thoughts regarding QR decomposition ?
Almost any method wiill work on a 100x100 matrix. On a modern PC the run time will be of the order of 1 second or less, so efficiency probably isn't important unless you want to solve thousands of matrices. QR decomposition should work fine.


Register to reply

Related Discussions
Positive Definite Matrix Linear & Abstract Algebra 8
Eigen Values of a 2x2 Matrix Calculus & Beyond Homework 4
Finding eigen values of a 2x2 matrix Calculus & Beyond Homework 2
Symmetric definite matrix is being pushed to indefinite in generalized eigen problem. Linear & Abstract Algebra 3
Ker Positive Definite Matrix Calculus & Beyond Homework 2