EXTREMELY WEIRD arithmetic bug

2 次查看(过去 30 天)
Hi there,
Basically I'm adding 0.001 to a number over and over again. The series starts with zero so the resulting number should ALWAYS be a multiple of 0.01, but instead it eventually degenerates into numbers like 5.874000000000296 - obviously NOT a multiple of 0.001!
Can you explain why? Can it REALLY BE that Matlab just can't cope with numbers of that precision? It's a mathematics package!!!
Run the following code to reproduce the bug. I'm running 7.10.0(R2010a).
Just run log=testcase() and then inspect the latter half of log in the variable editor- when you double-click on a cell to expand it you will be able to see the lack of precision.
function [ log ] = testcase( )
%RUN LIKE THIS: log=testcase()
%Then inspect the log in the variable editor.
runTime=20;
t=0;
dx=0.001;
i=1;
log(i,1)=t;
it=2;
for i=2:runTime/dx
t=t+dx;
log(i,1)=t;
end
end

采纳的回答

Arnaud Miege
Arnaud Miege 2011-8-3
  1 个评论
louise
louise 2011-8-3
I'm being an idiot. It's the floating point problem. I know all about the floating point problem.
SIGH.
I was stupid enough to assume that Matlab was using some futuristic high-precision arithmetic. Nope.
As usual, the problem is BAD ASSUMPTIONS.
Thanks for your fast help.

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Matrix Indexing 的更多信息

产品

Community Treasure Hunt

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

Start Hunting!

Translated by