计算方法 源码 道立特 牛顿差值 最小二乘法拟合函数

//import java.util.*; public class small{ static int m; static int n=3; public static void main(String args[]){ small small1=new small(); Scanner cin=new Scanner(System.in); System.out.println("输入多项式节点的个数："); m=cin.nextInt(); int[][] a=new int[m][2]; for(int i=0;i<m;i++){ System.out.println("输入第"+i+"个节点和函数值:"); for(int j=0;j<2;j++){ a[i][j]=cin.nextInt(); } } int[][] b=new int[n][n]; int[] y=new int[n]; small1.ercheng(b,y,a); small1.finish(b,y); } public void ercheng(int[][] b,int[] y,int[][] a){ int[] s=new int[m]; for(int x=0;x<m;x++){ s[x]=1; } for(int i=0;i<n;i++){ b[0][0]=m; for(int j=0;j<n;j++){ for(int l=0;l<m;l++){ s[l]=a[l][0]; for(int k=1;k<i+j;k++){ s[l]*=a[l][0]; } b[i][j]+=s[l]; } } } //求解Y for(int i=0;i<n;i++){ for(int j=0;j<m;j++){ int e=1; for(int k=1;k<=i;k++){ e*=a[j][0]; } int t=e*a[j][1]; y[i]+=t; } } } public void finish(int[][] b,int[] y){ float t=0; float l[][]=new float[n][n]; float u[][]=new float[n][n]; float Y[]=new float[n]; float x[]=new float[n]; for(int i=0;i<n;i++){ l[i][i]=1; } //u first row 求出U的第一行 for(int j=0;j<n;j++){ u[0][j]=b[0][j]; } //l first line 求出L的第一列 for(int i=1;i<n;i++){ l[i][0]=b[i][0]/u[0][0]; } //set the matrix L and U 依次求出L和U矩阵 for(int r=1;r<n;r++){ for(int i=r;i<n;i++){ t=0; for(int k=0;k<r;k++){ if(r==1) t=l[r][k]*u[k][i]; else { t+=l[r][k]*u[k][i]; } } u[r][i]=b[r][i]-t; } t=0; for(int i=r+1;i<n;i++){ for(int k=0;k<=r-1;k++){ t+=l[i][k]*u[k][r]; } l[i][r]=(b[i][r]-t)/u[r][r]; } } float e=0; Y[0]=y[0]; for(int i=0;i<n;i++){ e=0; for(int j=0;j<=i-1;j++){ e+=l[i][j]*Y[j]; } Y[i]=y[i]-e; } float q=0; x[n-1]=y[n-1]/u[n-1][n-1]; for(int i=n-1;i>=0;i--){ q=0; for(int j=i+1;j<n;j++){ q+=u[i][j]*x[j]; } x[i]=(Y[i]-q)/u[i][i]; } System.out.println(" C is:"); for(int i=0;i<n;i++){ System.out.println("c"+(i+1)+"="+x[i]); } System.out.println(); System.out.println(" "+"y="+x[0]+"+"+x[1]+"x"+"+"+x[2]+"x*x"); } }