numerical oddity with decimal multiplication
显示 更早的评论
Running R2022b, I noticed the following:
58*0.1-5.8
ans =
8.8818e-16
However,
58/(1/.1)-5.8
ans =
0
Why is that, and how can it be avoided? It messes up subsequent equality operations.
Some further testing shows that
58*0.1-5.8-4*eps
ans =
0
57*0.1-5.7
ans =
0
59*0.1-5.9
ans =
0
and in fact it's 0 for most numbers´I randomly tested, however
28*0.1-2.8-2*eps
ans =
0
Thanks in advance to anyone able and willing to clarify :)
1 个评论
"Why is that"
Because binary floating point numbers have finite precision, which has the effect that arithmetic operations on them necessarily accumulate floating point error. You (the programmer) need to understand that, and take that into account in your algorithms.
"and how can it be avoided?"
It can't be avoided (if you want to store data using binary floating point numbers).
This is worth reading as well:
回答（1 个）
the cyclist
2023-8-9
2 个投票
Possibly the most "popular" question on this forum. Read the Accuracy of Floating-Point Data section of the Floating-Point Numbers documentation.
(You will also get many other answers here, but they will all boil down to that same information.)
类别
在 帮助中心 和 File Exchange 中查找有关 Logical 的更多信息
2023-8-9
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
