拋物線法、牛頓法、弦截法求根實例

計算方法第四次計算實習題

​

7.用下列方法求$f(x)=x^3-3x-1=0在x_0=2附近的根，根的準確值x^{*}=1.87938524..,要求計算結果準確到四位有效數字$

（1）用牛頓法

（2）用弦截法，取$x_0=2,x_1=1.9$

（3）用拋物線法，取$x_0=1,x_1=3,x_2=2$

解題思路：按部就班，套公式編寫程序即可注意控制精度，要求準確到四位有效數字，即要求準確解和所得近似解誤差不超過$0.5*10^{-4}$，同時要注意迭代時的變量關係，以下是源代碼：

(1)牛頓法：

/**
 * @Title: newton.java
 * @Desc: TODO
 * @Package: root
 * @author: glm233
 * @date: 2020年5月25日 下午10:53:16
 * @version 1.0
 */

package root;

/**
 * @ClassName: newton
 * @Desc: TODO
 * @author: glm233
 * @date: 2020年5月25日 下午10:53:16
 * @version 1.0
 */
import java.util.Scanner;
public class newton {
	 public static void main(String[] args) {
	        Scanner scanner = new Scanner(System.in);
	        System.out.print("請輸入起始點：");
	        double x = scanner.nextDouble();
	        System.out.print("請輸入誤差允許範圍：");
	        double e = scanner.nextDouble();
	        System.out.print("請輸入最大迭代次數：");
	        int N = scanner.nextInt();
	        scanner.close();
	        double res = getEistimate(x,e,N);
	        System.out.println("牛頓法得到的解爲：" + res);
	    }

	    private static double getEistimate(double x,double e,int N) {
	        double result,x1;
	        for (int k = 1; true ; k++ ) {
	            if(f(x) == 0){
	                throw new RuntimeException("出現奇異情況....");
	            }
	            x1 = x - F(x)/f(x);
	            if(Math.abs(x1 - x) < e ){
	                System.out.println("迭代" + k + "次得到結果");
	                return x1;
	            }
	            if(k == N){
	                throw new RuntimeException("迭代失敗，已達最大迭代次數N....");
	            }
	            x = x1;
	        }
	    }

	    private static double F(double x) {
	        return x*x*x -3*x-1;
	    }

	    private static double f(double x) {
	        return 3*x*x -3;
	    }

}

運行結果：

（2）用弦截法，取$x_0=2,x_1=1.9$

/**
 * @Title: secant.java

 * @Desc: TODO
 * @Package: root
 * @author: glm233
 * @date: 2020年5月25日 下午11:05:56
 * @version 1.0
 */

package root;
import java.util.Scanner;
/**
 * @ClassName: secant
 * @Desc: TODO
 * @author: glm233
 * @date: 2020年5月25日 下午11:05:56
 * @version 1.0
 */
public class secant {
	   public static void main(String[] args) {
	        Scanner scanner = new Scanner(System.in);
	        System.out.print("請輸入初值1：");
	        double x0 = scanner.nextDouble();
	        System.out.print("請輸入初值2：");
	        double x1 = scanner.nextDouble();
	        System.out.print("請輸入誤差允許範圍：");
	        double e = scanner.nextDouble();
	        System.out.print("請輸入最大迭代次數：");
	        int N = scanner.nextInt();
	        scanner.close();
	        double res = getEistimate(x0,x1,e,N);
	        System.out.println("快速弦截法得到的解爲：" + res);
	    }

	    private static double getEistimate(double x0,double x1,double e,int N) {
	        double x2;
	        for (int k = 0; true ; k++ ) {
	            if(f(x0,x1) == 0){
	                throw new RuntimeException("出現奇異情況....");
	            }
	            x2 = x1 - F(x1)/f(x0,x1);
	            if(Math.abs(x2 - x1) < e){
	                System.out.println("迭代" + k + "次得到結果");
	                return x2;
	            }
	            if(k == N){
	                throw new RuntimeException("迭代失敗，已達最大迭代次數N....");
	            }
	            x0 = x1;
	            x1 = x2;
	        }
	    }

	    private static double f(double x0, double x1) {
	        return (F(x1) - F(x0))/(x1 - x0);
	    }

	    private static double F(double x) {
	    	return x*x*x -3*x-1;
	    }
}

運行結果：

]

（3）用拋物線法，取$x_0=1,x_1=3,x_2=2$

/**
 * @Title: parabolic.java

 * @Desc: TODO
 * @Package: root
 * @author: glm233
 * @date: 2020年5月25日 下午11:10:01
 * @version 1.0
 */

package root;

import java.util.Scanner;

/**
 * @ClassName: parabolic
 * @Desc: TODO
 * @author: glm233
 * @date: 2020年5月25日 下午11:10:01
 * @version 1.0
 */
public class parabolic {
	
	public final static double check=1.87938524;
	public static void main(String[] args) {
		
        Scanner scanner = new Scanner(System.in);
        System.out.print("請輸入初值1：");
        double x0 = scanner.nextDouble();
        System.out.print("請輸入初值2：");
        double x1 = scanner.nextDouble();
        System.out.print("請輸入初值3：");
        double x2 = scanner.nextDouble();
        System.out.print("請輸入誤差允許範圍：");
        double e = scanner.nextDouble();
        System.out.print("請輸入最大迭代次數：");
        int N = scanner.nextInt();
        scanner.close();
        double res = getEistimate(x0,x1,x2,e,N);
        System.out.println("精確解爲：" + check);
        System.out.println("拋物線法得到的解爲：" + res);
    }

    private static double getEistimate(double x0,double x1,double x2,double e,int N) {
        double x3,w=omiga(x0,x1,x2);
        for (int k = 0; true ; k++ ) {
            if(w*w<4*F(x2)*f(x0,x1,x2)){
                throw new RuntimeException("出現奇異情況....");
            }
            if(w>0) {
            	x3 = x2 - 2*F(x2)/(w+Math.sqrt(w*w-4*F(x2)*f(x0,x1,x2)));
            }
            else	x3 = x2 - 2*F(x2)/(w-Math.sqrt(w*w-4*F(x2)*f(x0,x1,x2)));
            if(Math.abs(x3 - check) < e){
                System.out.println("迭代" + k + "次得到結果");
                //System.out.println(x3-check);
                return x2;
            }
            if(k == N){
                throw new RuntimeException("迭代失敗，已達最大迭代次數N....");
            }
            //System.out.println(x3);
            x0 = x1;
            x1 = x2;
            x2 = x3;
        }
    }

    private static double omiga(double x0,double x1,double x2){
    	return f(x0,x1,x2)*(x2-x1)+f(x2,x1);
    	
    	
    }
    private static double f(double x0, double x1) {
        return (F(x1) - F(x0))/(x1 - x0);
    }
    
    private static double f(double x0, double x1,double x2) {
        return (f(x1,x2)-f(x1,x0))/(x2-x0);
    }

    private static double F(double x) {
    	return x*x*x -3*x-1;
    }
}

運行結果：