Can any one suggest me how to make a large matrics positive definate.
1 次查看(过去 30 天)
显示 更早的评论
I have matrices whose positive eigen values have to be calculated. In order to find them the matrices have to be positive definite. Kindly anyone suggest.
Can SDV be used for it? If yes then how?
0 个评论
回答(3 个)
Walter Roberson
2015-5-27
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.
1 个评论
John D'Errico
2015-5-27
编辑:John D'Errico
2015-5-27
The above trick of averaging the matrix with its transpose will indeed suffice to make it symmetric. But a symmetric matrix can still have negative eigenvalues. There will be no complex eigenvalues.
A = -eye(5);
A is perfectly symmetric, yet not at all positive definite, and trying to make it so by symmetrizing it will not suffice.
John D'Errico
2015-5-27
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.
0 个评论
Alfonso Nieto-Castanon
2015-5-27
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)
0 个评论
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!