c++ - C/C++ - 牛顿法的迭代计数器

Question

我创建了一个运行牛顿法的函数，用于逼近函数(定义为 f)的解。我的函数返回更好的根近似值就好了，但是它不会正确显示函数中执行的迭代次数。

这是我的代码:

#include <stdio.h> 
#include <math.h> 
#include <cstdlib>
#include <iostream>

double newton(double x_0, double newtonaccuracy);

double f(double x);

double f_prime(double x);

int main() 
{
   double x_0;  

   double newtonaccuracy;  

   int converged;  

   int iter;

   printf("Enter the initial estimate for x : ");

   scanf("%lf", &x_0);

   _flushall();

   printf("\n\nEnter the accuracy required : ");

   scanf("%lf", &newtonaccuracy);

   _flushall();


   if (converged == 1) 
      {
        printf("\n\nNewton's Method required %d iterations for accuracy to %lf.\n", iter, newtonaccuracy);

        printf("\n\nThe root using Newton's Method is x = %.16lf\n", newton(x_0, newtonaccuracy));
      } 

   else 
      {
        printf("Newton algorithm didn't converge after %d steps.\n", iter);
      }



      system("PAUSE");
} 


double newton(double x_0, double newtonaccuracy) 
{
   double x = x_0;

   double x_prev;

   int iter = 0;


   do 
   {
      iter++;


      x_prev = x;

      x = x_prev - f(x_prev)/f_prime(x_prev);

   } 
   while (fabs(x - x_prev) > newtonaccuracy && iter < 100);

   if (fabs(x - x_prev) <= newtonaccuracy)
   {
      int converged = 1;
   }  
   else
   {
      int converged = 0; 
   }   




    return x;
}  


double f(double x) {
       return ( cos(2*x) - x );
}  

double f_prime(double x) 
{
   return ( -2*sin(2*x)-1 ); 
}  

要尽可能具体，它是行:

printf("\n\nNewton's Method required %d iterations for accuracy to %lf.\n", iter, newtonaccuracy);

这给我带来了麻烦。每次我运行这个程序时，它都会说“牛顿法需要 2686764 次迭代......”但是这不可能是真的，前提是我已经正确编码(我的代码允许的最大迭代次数是 100)。

Answer 1

在 main 中使用的变量 iter 没有初始化，也没有在 newton 函数中使用，这里你使用了一个局部变量 iter 。您需要通过引用将 iter 传递给 newton，或者想办法从函数中返回它。

这是一个通过引用获取一些参数并修改它们的函数示例:

double foo(double& initial_value, int& iterations)
{
  initial_value *= 3.14159;
  iterations = 42;
  return initial_value/2.;
}

调用方:

double x + 12345.;
int iter = 0;
double y = foo(initial_value, iter);