Floating point numeric problem

Question

Qian Feng 2016-12-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/316183-floating-point-numeric-problem

评论： Star Strider 2016-12-17

采纳的回答： Star Strider

在 MATLAB Online 中打开

I encounter a problem which I cannot understand.

a = 1.228269:0.000001:1.22828;

b = 1.228265:0.000001:1.22828;

ismember(a,b)

>> ismember(a,b)

ans =

1×12 logical array
   0   1   0   1   1   1   1   1   1   1   1   1

It is clearly to see that the logical array is not compatible with the values in vectors a and b. An expert has suggested that this has to do with the floating point numerics, do anyone can explain what happened to this example and how to we circumvent it ? Thanks a lot !

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2016-12-9

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/316183-floating-point-numeric-problem#answer_246592

What happened is best explained in: Why is 0.3 - 0.2 - 0.1 (or similar) not equal to zero?

You can circumvent it using the ismembertol function (in R2015a and later versions).

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

Guillaume 2016-12-9

在 MATLAB Online 中打开

And Star's answer is not about ismember either. The problem is indeed one of floating point accuracy.

The number 0.000001 cannot be represented exactly in binary (the same way that you can't write 1/3 = 0.333... exactly in decimal), so every time you're adding 0.000001 you're accumulating some error (the same way if you did 0.333 + 0.333 + 0.333 = 0.999 in decimal). Because you're not starting at the same point (1.228269 vs 1.228265) the number at which the error accumulation is enough for it to matter differs for each sequence.

In your case, the simplest way to avoid this error accumulation is to multiply all your numbers by 1/0.000001 so you have integer boundaries and steps. All integers (up to flintmax) can be represented exactly as floating point, so you can't get any accumulation error. You can then do comparison on these integer values divided by 100000.

ia = 1228269:1228280
ib = 1228265:1228280
a = ia / 1e6;
b = ib / 1e6;
ismember(a, b)

Because the floating point value is obtained the same way for a and b (integer with no error / 1000000), the error is exactly the same for both numbers and they'll always compare equal for equal integers.

Qian Feng 2016-12-17