Different results for two equal functions

1 次查看（过去 30 天）
Beaya
Beaya 2013-10-27
评论： Beaya 2013-10-27
I am comparing these two functions:
f = sqrt(x.^2+1.0)-1.0
g = x.^2./(sqrt(x.^2+1.0)+1.0)
x = [8^-6 8^-7, 8^-8]
Even though f = g, I get slightly different values e.g. for x = 8^-7 I get f(x) = 1.136868377216 160 e-13 and g(x) = 1.136868377216 096 e-13.
I'm guessing this has something to do with the reduction of significant digits, however I don't know which results are more reliable. My guess is that function f = sqrt(x.^2+1.0)-1.0 is more precise because we only have one x in it but I am not sure...

请先登录，再回答此问题。

采纳的回答

Walter Roberson
Walter Roberson 2013-10-27
The two functions you show algebraically yield identical values only for x = 0.
  3 个评论
显示 1更早的评论
Walter Roberson
Walter Roberson 2013-10-27
Apologies, I had a typo when I did the calculations. The two do come out the same.
To more precision the 8^(-7) answer is
1.136868377216095673908443130184740364*10^(-13)
so your g(x) formula is more precise.
Beaya
Beaya 2013-10-27
Thank you but do you know what might be the reason to that? I mean why g(x) is more precise?

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Logical 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by