Matlab浮点数累加误差

## 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
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值