[讨论]matlab取整不整的一个问题

[austinlu]

在下面的程序中,第八行通过round或fix语句进行取整,但是发现在计算后,所得结果却并没有取整,减一个同大小的数字后也不是0,而是一个极小的数字.

clear %1
clc
p(1,:)=[6.84,125000]; %3
p(2,:)=[4.56,125000];
L=11.4; %5
l=0.57;
for j=1:length(p);
qoele=round(p(j,1)/l)+1 %8;
b=qoele*l-p(j,1)
a=l-(qoele*l-p(j,1))
end

运行结果为:
qoele =
13
b =
0.569999999999999
a =
5.55111512312578e-016
qoele =
9
b =
0.57
a =
-3.33066907387547e-016

其实这个程序用手算结果进行比较很明显的可以看出不是显示倍数的问题,我也用format等显示格式的函数进行了设置,但是还是找不到问题的所在.最后只能怀疑是round函数的问题.
我也用过fix函数取整,但是同样的问题还是存在.

[watcher]

这个不是round函数的问题，matlab做数值运算，a之所以不等于零是由于舍入误差引起的。简单的讲，matlab计算1/9*9得到的是个实数，而不是整数。
科学计算中恰当的误差是可以接受的，matlab有自己的误差控制方法，所以除非特殊的计算要求，否则不必过分注意舍入误差。

但是有些情况还是会带来困惑。比如，用matlab绘制一个带有奇点的函数，这时获得的结果往往和程序中使用的离散步长有关，即不同的离散距离对应不同的图形，这不是matlab计算有错，而是因为离散过程往往会遗漏奇点信息（事实上如果你强制计算奇点处的函数值，matlab会报错）。一种可选的解决方案是减小离散步长，但是不能从根本上解决问题。当我遇到这类问题时，通常是选择使用maple来得到函数图形。