## Matlab浮点数累加误差引起的判断错误
第一次发文,记录一下,用词不当,还望指正。
在使用Matlab的过程中,进行了如下的一次循环,(Q_site=[0,2],sub_Q初始化为1),如果使用(Q_site[1,sub_Q]+0.1==i)直接判断二者是否相等,在i循环至0.5时,由于浮点数累加导致的误差,此时的if值为假,无法进入if语句,从而无法得到后续的正确结果。在浏览了相关博客之后,采用abs(i-(Q_site(1,sub_Q)+0.1))<eps,通过eps判断当两者的误差在这个精度之内即为相等,但当i继续累加至1左右,累计的误差便超过了eps,由于我这里对精度要求在0.1,因此我选择将eps100来作为控制精度。同时在第二个if语句中,我的本意是当Q_site(1,sub_Q)的值比Q_site(1,sub_Q+1)的值小0.1时,便停止递增,但由于累加引起的损失误差,会导致这个差值比0.1略大,从而会在正确的累加基础上再累加一次,通过对第二个if语句进行计数可以验证,因此在判断条件中添加了一个控制精度,即100eps来确保正确的累加次数。可以在命令行窗口中输入“eps”获取eps的值。
for i=0:0.1:sum_span
if (abs(i-(Q_site(1,sub_Q)+0.1))<100*eps)&&(i>0)
value_of_shear(1,index)=value_of_shear(1,index-1)-0.1*Q_value(1,(sub_Q+1)/2);
if (Q_site(1,sub_Q)+0.1)+100*eps<Q_site(1,sub_Q+1)
Q_site(1,sub_Q)=Q_site(1,sub_Q)+0.1;
digit=digit+1;
else
if sub_Q<(num_Q-1)
sub_Q=sub_Q+2;
end
end
signal=1;
sign=1;
end
end