草庐IT

c++ - e^x 的无限和收敛到一个不正确的值

coder 2024-02-19 原文

我正在尝试实现一种无需使用阶乘即可计算 e^x 的方法。我通过获取每两个连续项之间的比率并通过将该比率乘以最后一项来计算下一项来做到这一点。所有这些都将添加到结果总和中，直到项小到无关紧要为止。

这适用于较小的 x 值。对于 x 的 “极端” 值，由于某种原因，这会中断。我试过 -50，它应该打印 1.9287498e-22(固定格式)。 我得到的是 GNU GCC 中的 -56676.4235303065 和 VC++ 中的 2041.8329628977。两者都大错特错。这是我的代码:

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;

int main()
{
    double x;
    cout << fixed << setprecision(10) ;
    cout << "Enter x: " << endl;
    cin >> x;
    double sum = 0;
    long long i = 0;
    double term = 1;
    do
    {
        sum += term;
        term *= (x / (++i));
    } while (fabs(term) > 1e-10);
    cout << sum << endl;
}

这是我的问题:导致此错误仅出现的问题是什么 x 的大值，我该如何解决这个错误？，为什么 GNU GCC 和 VC++2017得到完全不同的答案？

最佳答案

此求和中的最大项比 double 的精度高约 10 个数量级: ~ 2.92e+20与 double 相比的 13 平方英尺或者。这些术语中的误差范围本身比您的求和结果大 ~30 个数量级。

因此，您的级数没有正确收敛也就不足为奇了，因为连续的项(相反的符号)可能不会抵消理论量。即使使用一些数值技巧，例如 Kahan-Neumaier 求和和在添加前对项进行排序，结果仍然只能减少到 ~6000。 .请注意，对于正 x 不会发生这种情况因为连续的任期不需要取消。

克服这个问题的一种方法是对 x 施加一个小的下限。 , 并使用平方取幂按比例放大到 x 的正确值.

更新:上述方法的实现:

// integer exponentiation by squaring (won't explain here)
double pow_square(double x, unsigned a)
{
   double r = 1.0;
   while (a > 0)
   {
      if (a % 2 == 1)
      {
         a--;
         r *= x;
      }
      a /= 2;
      x *= x;
   }
   return r;
}

// original method
double exp_original(double x, double e)
{
   double sum = 0.0;
   unsigned i = 0;
   double term = 1.0;
   do
   {
      sum += term;
      term *= (x / (++i));
   } while (fabs(term) > e);
   return sum;
}

// new adaptive method
double exp_new(double x, double e)
{
   static const double min_X = -3;

   // if within limit, simply use original function
   if (x >= min_X)
      return exp_original(x, e);

   // compute smallest possible scaling coefficient
   unsigned s = (unsigned)(x / (-min_X) + 0.5);
   double p = exp_original(x / s, e);
   return pow_square(p, s);
}

大范围测试 x值证实新方法可以更好地处理极值(负)情况:

x    | exp (C-library)    exp_original       exp_new
-------------------------------------------------------------
-10  | 4.53999297625e-05  4.53998989141e-05  4.53999299001e-05
-30  | 9.35762296884e-14  6.10299992426e-06  9.35762292245e-14
-50  | 1.92874984796e-22  2041.8329629       1.92874983803e-22
-60  | 8.7565107627e-27   722745700.93       8.75651067587e-27
-80  | 1.80485138785e-35  2.45082011705e+17  1.80485137011e-35
-100 | 3.72007597602e-44  8.1446527451e+25   3.72007589785e-44
-150 | 7.17509597316e-66  -9.14622659954e+47 7.1750957953e-66
-200 | 1.38389652674e-87  7.69097143891e+69  1.38389648613e-87

关于c++ - e^x 的无限和收敛到一个不正确的值，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/49987973/

不正ampdoublecodesectionc++math

有关c++ - e^x 的无限和收敛到一个不正确的值的更多相关文章

  1. ruby - 树顶语法无限循环 - 2

    我脑子里浮现出一些关于一种新编程语言的想法，所以我想我会尝试实现它。一位friend建议我尝试使用Treetop(Rubygem)来创建一个解析器。Treetop的文档很少，我以前从未做过这种事情。我的解析器表现得好像有一个无限循环，但没有堆栈跟踪；事实证明很难追踪到。有人可以指出入门级解析/AST指南的方向吗？我真的需要一些列出规则、常见用法等的东西来使用像Treetop这样的工具。我的语法分析器在GitHub上，以防有人希望帮助我改进它。class{initialize=lambda(name){receiver.name=name}greet=lambda{IO.puts("He

  2. ruby - 使用 Vim Rails，您可以创建一个新的迁移文件并一次性打开它吗？ - 2

    使用带有Rails插件的vim，您可以创建一个迁移文件，然后一次性打开该文件吗？textmate也可以这样吗？ 最佳答案 你可以使用rails.vim然后做类似的事情::Rgeneratemigratonadd_foo_to_bar插件将打开迁移生成的文件，这正是您想要的。我不能代表textmate。 关于ruby-使用VimRails，您可以创建一个新的迁移文件并一次性打开它吗？，我们在StackOverflow上找到一个类似的问题： https://sta

  3. ruby-on-rails - Rails - 一个 View 中的多个模型 - 2

    我需要从一个View访问多个模型。以前，我的links_controller仅用于提供以不同方式排序的链接资源。现在我想包括一个部分(我假设)显示按分数排序的顶级用户(@users=User.all.sort_by(&:score))我知道我可以将此代码插入每个链接操作并从View访问它，但这似乎不是“ruby方式”，我将需要在不久的将来访问更多模型。这可能会变得很脏，是否有针对这种情况的任何技术？注意事项:我认为我的应用程序正朝着单一格式和动态页面内容的方向发展，本质上是一个典型的网络应用程序。我知道before_filter但考虑到我希望应用程序进入的方向，这似乎很麻烦。最终从任何

  4. ruby-on-rails - 渲染另一个 Controller 的 View - 2

    我想要做的是有2个不同的Controller，client和test_client。客户端Controller已经构建，我想创建一个test_clientController，我可以使用它来玩弄客户端的UI并根据需要进行调整。我主要是想绕过我在客户端中内置的验证及其对加载数据的管理Controller的依赖。所以我希望test_clientController加载示例数据集，然后呈现客户端Controller的索引View，以便我可以调整客户端UI。就是这样。我在test_clients索引方法中试过这个:classTestClientdefindexrender:template=>

  5. ruby-on-rails - 如何优雅地重启 thin + nginx？ - 2

    我的瘦服务器配置了nginx，我的ROR应用程序正在它们上运行。在我发布代码更新时运行thinrestart会给我的应用程序带来一些停机时间。我试图弄清楚如何优雅地重启正在运行的Thin实例，但找不到好的解决方案。有没有人能做到这一点？ 最佳答案 #Restartjustthethinserverdescribedbythatconfigsudothin-C/etc/thin/mysite.ymlrestartNginx将继续运行并代理请求。如果您将Nginx设置为使用多个上游服务器，例如server{listen80;server

  6. ruby-on-rails - 如何使用 instance_variable_set 正确设置实例变量？ - 2

    我正在查看instance_variable_set的文档并看到给出的示例代码是这样做的:obj.instance_variable_set(:@instnc_var,"valuefortheinstancevariable")然后允许您在类的任何实例方法中以@instnc_var的形式访问该变量。我想知道为什么在@instnc_var之前需要一个冒号:。冒号有什么作用？ 最佳答案 我的第一直觉是告诉你不要使用instance_variable_set除非你真的知道你用它做什么。它本质上是一种元编程工具或绕过实例变量可见性的黑客攻击

  7. ruby-on-rails - 如果 Object::try 被发送到一个 nil 对象，为什么它会起作用？ - 2

    如果您尝试在Ruby中的nil对象上调用方法，则会出现NoMethodError异常并显示消息:"undefinedmethod‘...’fornil:NilClass"然而，有一个tryRails中的方法，如果它被发送到一个nil对象，它只返回nil:require'rubygems'require'active_support/all'nil.try(:nonexisting_method)#noNoMethodErrorexceptionanymore那么try如何在内部工作以防止该异常？ 最佳答案 像Ruby中的所有其他对象

  8. ruby - 为什么 SecureRandom.uuid 创建一个唯一的字符串？ - 2

    关闭。这个问题需要detailsorclarity.它目前不接受答案。想改进这个问题吗？通过editingthispost添加细节并澄清问题.关闭8年前。Improvethisquestion为什么SecureRandom.uuid创建一个唯一的字符串？SecureRandom.uuid#=>"35cb4e30-54e1-49f9-b5ce-4134799eb2c0"SecureRandom.uuid方法创建的字符串从不重复？

  9. ruby - 如果指定键的值在数组中相同，如何合并哈希 - 2

    我有一个这样的哈希数组:[{:foo=>2,:date=>Sat,01Sep2014},{:foo2=>2,:date=>Sat,02Sep2014},{:foo3=>3,:date=>Sat,01Sep2014},{:foo4=>4,:date=>Sat,03Sep2014},{:foo5=>5,:date=>Sat,02Sep2014}]如果:date相同，我想合并哈希值。我对上面数组的期望是:[{:foo=>2,:foo3=>3,:date=>Sat,01Sep2014},{:foo2=>2,:foo5=>5:date=>Sat,02Sep2014},{:foo4=>4,:dat

  10. ruby - 检查字符串是否包含散列中的任何键并返回它包含的键的值 - 2

    我有一个包含多个键的散列和一个字符串，该字符串不包含散列中的任何键或包含一个键。h={"k1"=>"v1","k2"=>"v2","k3"=>"v3"}s="thisisanexamplestringthatmightoccurwithakeysomewhereinthestringk1(withspecialcharacterslike(^&*$#@!^&&*))"检查s是否包含h中的任何键的最佳方法是什么，如果包含，则返回它包含的键的值？例如，对于上面的h和s的例子，输出应该是v1。编辑:只有字符串是用户定义的。哈希将始终相同。 最佳答案

随机推荐