最近遇到一个题目,就是不使用Math方法去求解一个数的立方根和平方根,搜索了好多,现在记录下自己掌握的方法。
首先,在编码之前,我们得熟悉 什么是牛顿迭代法。下面我为大家找了两处链接。大家可以去看一下,然后我再说我的程序编写,如果只看程序,不熟悉这个算法的人是看不懂的。
下面 我来说说 我所理解的: 经过以上参考 以立方根举例 我们求a的立方根 首先列式 a=x^3 (a等于x的3次方) 。根据牛顿迭代法 我们把它转成 f(x) = x^3-a; 那么 f'(x) =3*x^2;
则 x1 = x0-f(x0)/f'(x0); (我们求次方就是:x1 = x0 -(x0^3-a)/(3*x0^2); 化简可得 :x1 = (2*x0+a/x0^2)/3。) x0怎么来的呢? 我们可以先假设一个x0 = a; 然后计算得出的x1-x0是否小到极致。例如 x1-x0<0.0000000001时,x1就是那个立方根。
我来写一个总的例子: 求6的平方根
1. 6 = x^3;
2. f(x) = x^3-6;
3.f'(x) = 3*x^2;
4.假设答案是x0=6;
5 . 求一波 x1 = x0-f(x0)/f'(x0); 下一步:x1 = x0 -(x0^3-a)/(3*x0^2);下一步:x1 = (2*x0+a/x0^2)/3;(这一步式子我用在了代码中)
6. 解得 x1 = 4.055555555555555;(小数位多吧,用计算机算又不是你算)
7.刚才猜的答案x0是6 ,x0-x1>0.0000000001 (精度自己设) 了, 不够准确啊这答案,怎么办 木有事
8. 让 x0=4.055555555555555; 再求一遍x1
9。此时 求得x1=2.825302502727911; x0-x1>0.00000001 还是不准确 继续求
10。。。。不断求 终于 x0-x1<0.00000001,结束 x1就是最终答案。
说了这么多 上代码:
import Java.util.Scanner;
/*
* 牛顿迭代法
*
* 求解立方根
*计算一个数字的立方根,不使用Math方法
*/
public class ClassTest09 {
public static void main(String[] args) {
System.out.println("请输入一个数字");
Scanner sc= new Scanner(System.in);
double a= sc.nextDouble();
double num =getLiFangGen(a);
System.out.println(a+"的立方根为:"+num);
double num1 = getPingFangGen(a);
System.out.println(a+"的平方根为:"+num1);
}
private static double getLiFangGen(double num){
if(num==0){
return 0;
}
double num1 = num;
double num2 = ((2*num1)+(num/(num1*num1)))/3;
System.out.println(num2);
while(num1-num2>0.00000001 || num1-num2<(-0.00000001)){
num1 = num2;
num2 = ((2*num1)+(num/(num1*num1)))/3;
System.out.println(num2);
}
return num2;
}
private static double getPingFangGen(double num){
if(num==0){
return 0;
}
if(num<0){
System.out.println(num+"没有平方根");
System.exit(0);
}
double num1 = num;
double num2 = (num1+(num/num1))/2;
while(num1-num2>0.00000001){
num1 = num2;
num2 = (num1+(num/num1))/2;
}
return num2;
}
}