累积误差

Question

我有一个非常直接的问题。以下代码打印出摄氏度和华氏度。 我的问题是关于它迭代的次数。对于较小的数字，例如从 0 开始，到 10 停止，步长为 1.1。循环完成后，它将打印出正确的迭代次数。

但对于 0-11000000 的大数，使用步骤 1.1 会打印出错误的迭代次数。为什么会发生这种情况？由于 1100000/1.1 应该约为 1000001，但我得到 990293。

#include <iostream>
#include <iomanip>

using namespace std;

int main()
{

   float start, stop, step;
   int count = 0;

   cout << "start temperature: ";
   cin >> start;
   cout << "stop temperature: ";
   cin >> stop;
   cout << "step temperature: ";
   cin >> step;

   cout << setw(10) << "celsius" << setw(15) << "fahrenheit" << endl;
   cout << setw(25) << "celsius" << setw(15) << "fahrenheit" << endl;

   while(start <= stop)
   {
      count++;
      float c, f;
      c = (5.0/9)*(start-32);
      f = 32+(9.0/5)*start;
      cout << setw(10) << fixed << setprecision(2) << c << setw(15) << start << setw(15) << fixed << setprecision(2) << f  << " count: " << count << endl;

      start = start + step;
   }
   cout << "The program loop made " << count << " iterations." << endl;

   return 0;
}

Answer 1

浮点舍入错误。从本质上讲，浮点数并不是 100% 准确的表示，你所做的每一次计算都会存在错误，并且当你反复添加它们时，你会添加越来越多的错误。您应该做的是计算一次步数，将其存储在一个整数中，然后循环多次。

Answer 2

根据记录，清理后的版本如下所示：

#include <iostream>
#include <iomanip>

using namespace std;

int main()
{

   double start, stop, step;

   cout << "start temperature: ";
   cin >> start;
   cout << "stop temperature: ";
   cin >> stop;
   cout << "step temperature: ";
   cin >> step;

   cout << setw(10) << "celsius" << setw(15) << "fahrenheit" << endl;

   unsigned steps = (stop - start) / step;

   for(unsigned i = 0; i < steps; ++i)
   {
      double temp = start + i * step;
      double c = (5.0 / 9.0) * (temp - 32.0);
      double f = 32.0 + (9.0 / 5.0) * temp;
      // This is a real good example of why people hate <iostream> formatting.
      // If you want formatting that's quite different from the default, it
      // gets too verbose too fast. Using C stdio: 
      //printf("%10.2f%15.2f\n", c, f);
      cout << setw(10) << fixed << setprecision(2) << c
           << setw(15) << fixed << setprecision(2) << f << endl;
   }

   cout << "The program loop made " << steps << " iterations." << endl;

   return 0;
}

这种循环风格的主要好处是每次迭代（输出除外）与顺序无关，因此可以展开并（理论上）并行化。