mintwo51.java.bak

数值分析之计算方法与实验5：曲线拟合 by java

class Mintwo51{
  static double sum1(int k,double x[][]){
    double sum=0.0;
    for(int j=0;j<=x.length;j++)
      sum=sum+Math.pow(x[0][j],k); 
    return sum;
  }
  static double sum2(int k,double x[][]){
    double sum=0;
    for(int j=0;j<=x.length;j++)
      sum=sum+Math.pow(x[0][j],k)*x[1][j];
    return sum;
  }        
  static double [][] zhenguix(int m,double x[][]){
    int n=x.length;
    double s[]=new double[2*m+1];
    double a[][]=new double[m+1][m+1];
    for(int k=0;k<=2*m;k++)
      s[k]=sum1(k,x);
    for(int i=0;i<=m;i++)
      for(int j=0;j<=m;j++)
         a[i][j]=s[i+j];
   
    double y[]=new double[m+1];
    for(int k=0;k<=m;k++)
      y[k]=sum2(k,x);
   
    
    double xx[][]=new double[m+2][m+2];
    for(int i=0;i<=m;i++)
      for(int j=0;j<=m;j++)
        xx[i][j]=a[i][j];
    for(int j=0;j<=m+1;j++)xx[m+1][j]=0;
    for(int i=0;i<m+1;i++)xx[i][m+1]=y[i];
    return xx;
  }            
     
  public static void main(String[]args){
    double x[][]={{1.0,1.5,2.0,2.5,3.0,3.5,4.0,4.5,5.0,5.5,6.0,6.5,7.0,7.5,8.0},
    {33.4,79.5,122.65,159.05,189.15,214.15,238.65,252.50,267.55,280.50,296.65,301.40,310.40,318.15,325.15}};
    double xx[][]=zhenguix(2,x);
    System.out.println("a,b,c分别为"+xx);
   
  }
}