在这个问题中,我正在解决Matlab中的数值计算问题,并希望获得经验,以日后如何避免这种问题/错误
例如,让我们考虑遵循简单的代码
t = 0.4 + 0.1 - 0.5
t =
0
它工作正常,但是
u = 0.4 - 0.5 + 0.1
u =
2.7756e-17
当然要记住它也是0,但是为什么在第一次计算中没有得到相同的结果?或者有什么不同?也请看一下
v = (sin(2*pi) = = sin(4*pi))
v = (sin(2*pi)==sin(4*pi))
v =
0
它表明正弦函数不是周期性的,因此在这种情况下的一般建议是引入一些ε?
V=((sin(2*pi)-sin(4*pi))<eps)
V =
0
或者
EPS=0.000000000000001
EPS =
1.0000e-15
>> V=((sin(2*pi)-sin(4*pi))<EPS)
V =
1
请帮我
得到这些结果是正常的,因为Matlab中的浮点相对精度是
eps('double')
ans =
2.2204e-16
对于V=((sin(2*pi)-sin(4*pi))<eps)
,因为
sin(2*pi)-sin(4*pi)
ans =
2.4493e-16
大于eps('double')
,因此其结果将是V=0
。
而且V=((sin(2*pi)-sin(4*pi))<EPS)
因为EPS>2.4493e-16
,这样的结果将会是V=1
。
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句