Can any one suggest me how to make a large matrics positive definate.
2 views (last 30 days)
Walter Roberson on 27 May 2015
The trick is to use
pdmat = (YourMatrix + YourMatrix.')/2;
Or is the question to find out if a given matrix is positive definite? If so then note that the positive definite tests are very sensitive to the exact bit values of floating point numbers, so a matrix that looks positive definite might not be because of a single bit difference between two numbers. The above trick is routinely used on matrices that "should" be positive definite in order to force the last bits to be the same in the symmetric positions.
John D'Errico on 27 May 2015
Your question is not very clear.
You do not need for a matrix to be positive definite to compute the eigenvalues of that matrix. Just use EIG. SVD is not necessary.
Knowing what you will do with those eigenvalues might help us to help you more.
Alfonso Nieto-Castanon on 27 May 2015
the matrix B defined as:
B = (A+A')/2 + lambda*speye(size(A,1));
will have eigenvalues eig(B) equal to real(eig(A))+lambda, so for a lambda high enough the resulting matrix B will be positive definite.
How that would help you, though, is another question (I am not sure I understand what you are trying to do exactly)