Java实现牛顿插值

喜欢ヅ旅行 2023-09-23 20:01 161阅读 0赞

![在这里插入图片描述][739af135a0d74f84b68d4d8a4fd0eb20.png]

public class Test {
        
        static double[][] format;
    
        public static void main(String[] args) {
        
            // TODO Auto-generated method stub
            Scanner in = new Scanner(System.in);
            System.out.print("输入插入总点数:");
            int n=in.nextInt();//插入点总数
            double[][] xy=new double[2][n];//二维数组来存储x和y
            for(int j=0;j<n;j++) {
        
                System.out.print("输入第"+(j+1)+"个点的x:");
                xy[0][j]=in.nextDouble();//第0行放各点的x值
                System.out.print("输入第"+(j+1)+"个点的y:");
                xy[1][j]=in.nextDouble();//第1行放各点的y值
            }
            creatformat(xy, n);
            System.out.println("输入计算几个点：");
            int v=in.nextInt();//所求点的总数
            for(int i=0;i<v;i++) {
        
                System.out.println("输入x= ");
                double x=in.nextDouble();
                System.out.print("f("+x+")≈"+"Nn("+x+")"+"=");
                System.out.printf("%.5f", Nn(xy, n, x));
                System.out.println();
            }
        }
    
    
        static void creatformat(double[][] xy,int n) {
        
            format=new double[n-1][n-1];
            for(int i=0;i<n-1;i++) {
        
                for(int j=0;j<n-1-i;j++) {
        
                    if(i==0) {
        
                        format[i][j]=(xy[1][j]-xy[1][j+1])/(xy[0][j]-xy[0][j+1]);
                    }
                    else{
        
                        format[i][j]=(format[i-1][j]-format[i-1][j+1])/(xy[0][j]-xy[0][j+i+1]);
                    }
                }
            }
        }
    
        static double Nn(double[][] xy,int n,double x) {
        
            //牛顿插值公式
            double sum=xy[1][0];//求和的初始值为f(X0)
            for(int i=0;i<n-1;i++) {
        
                double c=1;
                for(int j=0;j<i+1;j++) {
        
                    c=(x-xy[0][j])*c;
                }
                sum=sum+format[i][0]*c;
            }
            return sum;
        }
    }

运行结果得：  
![在这里插入图片描述][c72825c84728456391ec45d24e7de672.png]

[739af135a0d74f84b68d4d8a4fd0eb20.png]: https://img-blog.csdnimg.cn/739af135a0d74f84b68d4d8a4fd0eb20.png
[c72825c84728456391ec45d24e7de672.png]: https://img-blog.csdnimg.cn/c72825c84728456391ec45d24e7de672.png