colon operator rounding problem
显示 更早的评论
Hello everyone!
I encountered a problem in one of my code, using the semi colon operator like this:
a = 0:0.1:120;
will not give me exactly what I want, it does return 0, 0.1, 0.2 etc, but with a small imprecision (equal to eps actually)
the following code :
a = 0:0.1:120;
disp(a(20));
disp(a(20)-1.9)
isequal(a(20),1.9)
is returning:
1.9
2.22044604925031e-16
ans =
0
Any help ? I really need this isequal(a(20),1.9) to return 1...
thanks !
2 个评论
Ambroise
2015-7-10
The problem is the real value of 0.1...
If you really want to see what the floating point values really are, try this:
采纳的回答
bio lim
2015-7-10
Try it like this.
a=(0:1200)/10;
disp(a(20));
disp(a(20)-1.9);
isequal(a(20),1.9)
Your previous code was not returning 1 because of rounding errors when doing finite precision arithmetic.
3 个评论
Ambroise
2015-7-10
Ambroise
2015-7-10
If you multiply it by decimals, such as 0.1, again you get rounding errors. That is why, I specifically wrote 10 in the first place. If you actually check your initial code, you can see that until a(4), it returns 0 but starts getting rounding error from that point.
1/0.1 returns 10 because, 1 is defined and when you divide you are going to get 10 without a rounding error.
更多回答（2 个）
Use
a = linspace(0, 120, 1201);
But in general don't use
a(20) == 1.9
but
abs(a(20) - 1.9) <= eps
If you don't want to do it this way, just define
a(20) = 1.9;
Steven Lord
2015-7-10
0 个投票
See question 1 in the Mathematics section of the FAQ for a more detailed explanation of this behavior.
类别
在 帮助中心 和 File Exchange 中查找有关 Matrices and Arrays 的更多信息
产品
Translated by 
