eig() shows positive eigenvalues while chol() failed

Question

Allen Lau 2015-5-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/218341-eig-shows-positive-eigenvalues-while-chol-failed

评论： John D'Errico 2016-6-18

Hi guys,

I was LU-decomposing a supposed-to-be positive-definite matrix using "chol" in a loop. However, occasionally it failed with the warning

 " Error using chol
   Matrix must be positive definite. "

eig() shows that the smallest eigenvalue is 1.8*10^(-4) which I think should be well above zero. BTW, the largest one is 7.7*10^9.

How can I get around this dilemma? Your suggestions are greatly appreciated. Many thanks.

6 个评论
显示 4更早的评论隐藏 4更早的评论

John D'Errico 2015-5-24

在 MATLAB Online 中打开

First of all, it is very often the case that a matrix that LOOKS symmetric is actually not so. You can test that. Try this:

norm(A - A')

or

max(abs(A - A'))

Do these tests return exact zeros? Very often the computations that are done don't create an exactly symmetric matrix. This happens because of floating point trash, that when you change the order of operations, you don't end up with exactly the same result. And in a matrix multiply, it is the BLAS that can do things like that to you.

As far as "twisting" the matrix, the difference in this case will be in the least significant bits of those numbers. If it was failing the chol test, that means your matrix is not SPD as far as MATLAB is concerned. Wanting that to be different is not enough.

Your matrix, if the spread of the eigenvalues is only 1e9 down to 1e-4 (so only 13 powers of 10) would not seem to be numerically singular. That dynamic range would normally need to be on the order of 1e16 for that to be a problem. So my guess is that your matrix is not indeed symmetric, even though you think it is so.

In any case, I have a hard time seeing why you feel it to be a problem to tweak the matrix to be usable. It was not SPD as you created it, so it was NOT truly a valid matrix for your purposes.

Ricardo Almeida 2016-6-18

I´m having the same problem. Without going into peculiarities of decomposition methods, I think it might be some technical issue. I'm running chol function in two different computers, both Windows 7 64bits and matlab 2015a. One flags a positive definite matrix and other don't (Maybe it's a coincidence but always return the number of columns).

John D'Errico 2016-6-18

Read my comment to the other question where you made the same comment.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

John D'Errico 2015-5-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/218341-eig-shows-positive-eigenvalues-while-chol-failed#answer_180220

编辑：John D'Errico 2015-5-24

I wrote nearestSPD (download it from the file exchange) for this purpose. It will tweak the values of your matrix slightly so that chol should always succeed. Of course, if the required tweak is necessarily large, it will still make a minimal change so that the result is as close as possible, yet will be positive definite.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

eig() shows positive eigenvalues while chol() failed

6 个评论
显示 4更早的评论隐藏 4更早的评论

采纳的回答

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

eig() shows positive eigenvalues while chol() failed

6 个评论 显示 4更早的评论隐藏 4更早的评论

采纳的回答

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

6 个评论
显示 4更早的评论隐藏 4更早的评论

0 个评论
显示 -2更早的评论隐藏 -2更早的评论