ITPub博客

首页 > 应用开发 > IT综合 > 牛顿下山法(C++版)

牛顿下山法(C++版)

原创 IT综合 作者:chenjinglzz 时间:2007-10-07 15:25:28 0 删除 编辑

本程序为本人参考一些资料后编写,在visual c++6.0下调试通过。

[@more@]//////////////
////由于牛顿迭代法的收敛性依赖于初值的选取,故有下面的
////牛顿下山法:其具体原理为
////牛顿下山公式为:x_(k+1)=x_k-r*f(x_k)/f'(x_k),其中r的范围为(0,1],称为下山因子
////具体做法为:
////在集合{1,1/2,1/4,...,1/2^n,...}中挑选下山因子,使|f(x_k)|>|f(x_(k+1))|,否则重新选取初值
////相应的迭代函数为:g(x)=x-r*f(x)/f'(x),x=g(x)
////输出每一次迭代的结果,将迭代模块做全面
//////////////
#include
#include
using namespace std;
///////////
double fun(double x);
double dfun(double x);
void NewtonIterative(double(*pfun)(double x),double (*pdfun)(double x));
//////////
//////////
int main()
{
NewtonIterative(fun,dfun);
return 0;
}
//////////
//////////
double fun(double x)
{
return (x*x*x-x-1);
}
//////////
double dfun(double x)
{
return (3*x*x-1);
}
//////////
void NewtonIterative(double(*pfun)(double x),double (*pdfun)(double x))
{
int maxn,maxflag,k,n;
k=1;n=0;
double x0,x,temp;
double err;
cout< exit(1);
}
x=temp;
if(n==0)
cout<<" 第"< else
{
cout<<" 第"< cout<<" 下山因子为:r=1/2^"< }
n=0;
++k;
}while(fabs(x-x0)>err&&k if(k>=maxflag)
cout<

来自 “ ITPUB博客 ” ,链接:http://blog.itpub.net/11355887/viewspace-974947/,如需转载,请注明出处,否则将追究法律责任。

请登录后发表评论 登录
全部评论

注册时间:2007-12-10

  • 博文量
    17
  • 访问量
    334554