牛顿迭代法求解立方根

//计算一个数字的立方根,不使用库函数
//详细描述:
//接口说明
//原型:
//public static double getCubeRoot(double input)
//输入:double 待求解参数

//返回值:double  输入参数的立方根

牛顿迭代法原理:

对于求a的立方根,可以设f(x)=x^3-a,从而转换成求解f(x)=0,即求方程的根。

f(x)展开泰勒公式:f(x)=f(x0)+(x-x0)*f'(x0)=0,得出的x=x0-f(x0)/f'(x0)=g(x0),此时递归调用该式子可以逐步接近于最终结果。

为什么会接近于最终结果?当然,牛顿迭代并不是无条件收敛的。

首先,要保证f'(x0)!=0,这样f(x)=0将等价于x=x-f(x)/f'(x)。而x_k+1=xk-f(xk)/f'(xk)=g(xk),由迭代过程收敛性定理可得abs(g'(xk))<=L<1(具体证明请读者自行查找)。

g'(x*)=f(x*)*f''(x*)/[f'(x)]^2,x*为f(x)的一个根,所以g‘(x*)=0,g''(x*)=f''(x*)/f'(x*)!=0,只要f’‘(x*)!=0,则牛顿迭代法收敛。

#include <stdio.h>
#include <stdlib.h>

#define jingdu 0.0001
double newton_diedai(double a){
	double xn,xn_1;

	xn=1;
	xn_1 = xn-((xn*xn*xn-a)/(3*xn*xn));
	while(xn-xn_1>jingdu || xn-xn_1<-jingdu){
		xn = xn_1;
		xn_1 = xn-((xn*xn*xn-a)/(3*xn*xn));
	}
	return xn_1;
}

int main(void){
	double a;

	scanf("%lf", &a);
	a = newton_diedai(a);
	printf("%.1f", a);

	system("pause");
	return 0;
}


  • 0
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值