Is this a new bug?

A'*A, its eigenvalues' calculated results has a negetive value:
C =
0.0215 -0.0014 -0.2572
-0.0014 0.0006 0.1029
-0.2572 0.1029 18.0000
eig(C'*C)
ans =
-1.22819777351878e-17
0.000316768541555439
324.153592202518

3 个评论

which release and OS are you using?
1.31720860254795e-10
0.000317607460000373
324.153629422408
Before asking for release make sure you guys exchange data not from screen display.
num2hex can show us the exact values you used

回答(1 个)

Bruno Luong
Bruno Luong 2020-7-21
编辑:Bruno Luong 2020-7-21
Numerical inaccuracy limitation. Nothing new in this front and no one consider them as BUG.
Sometime symmetrize the matrix help
D = C'*C;
eig(0.5*(D + D'))
Further reading if the above is not enough to solve the "bug": nearest SPD

此问题已关闭。

标签

提问:

2020-7-21

关闭:

2021-8-20

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by