Newton 牛顿迭代法求根公式及其改进算法

最新推荐文章于 2022-12-26 18:25:28 发布

油醋三椒

最新推荐文章于 2022-12-26 18:25:28 发布

阅读量4k

点赞数 3

分类专栏：计算方法文章标签：算法 Newton 牛顿迭代法求根

本文链接：https://blog.csdn.net/weixin_42664622/article/details/103682278

版权

计算方法专栏收录该内容

13 篇文章 16 订阅

订阅专栏

简介

在这里插入图片描述

代码

#include<bits/stdc++.h>
using namespace std;

double f(double x)
{
//	cout<<"f(x):  "<<x*x*x+2*x*x+10*x-20<<endl;
	return x*x*x+2*x*x+10*x-20;
 } 

double f_(double x){
//	cout<<"f_x:  "<<x*x+4*x+10<<endl;
	return x*x*3+4*x+10;
} 

int main()
{
	double x;
	cout<<"输入计算的初值点"<<endl;
	cin>>x;
	int maxx;
	cout<<"输入最大迭代次数"<<endl;
	cin>>maxx;
	for(int i=0;i<maxx;i++)
	{
		cout<<"i: "<<i<<" x: "<<x<<endl;
//		cout<<x<<" "<<f(x)<<" "<<f_(x)<<endl; 
		x = x - f(x)/f_(x); 
	} 
	return 0;
}

结果（默认函数为：x * x * x+2 * x * x+10 * x - 20 ）：
在这里插入图片描述

Newton下山法

在这里插入图片描述

代码

#include<bits/stdc++.h>
using namespace std;

double f(double x)
{
//	cout<<"f(x):  "<<x*x*x+2*x*x+10*x-20<<endl;
	return x*x*x+2*x*x+10*x-20;
 } 

double f_(double x){
//	cout<<"f_x:  "<<x*x+4*x+10<<endl;
	return x*x*3+4*x+10;
} 

int main()
{
	double x;
	cout<<"输入计算的初值点"<<endl;
	cin>>x;
	int maxx;
	cout<<"输入最大迭代次数"<<endl;
	cin>>maxx;
	double k ;
	cout<< "输入精度要求" <<endl;
	cin>>k ;
	cout<<"i: "<<0<<" x: "<<x<<endl;
	for(int i=0;i<maxx;i++)
	{
//		cout<<x<<" "<<f(x)<<" "<<f_(x)<<endl; 
		double z = 1.0,y; 
		for(int j=0;;j++)
		{
			y = x - z*f(x)/f_(x); 
			z = z*0.5;
			if(fabs(f(y))<fabs(f(x))) break;
		}
		if(fabs(y-x)<k ) { x = y;cout<<"i: "<<i<<" x: "<<x<<endl; break;		}
		x = y;
		cout<<"i: "<<i<<" x: "<<x<<endl;
	} 
	return 0;
}