最小二乘法 .txt

# include <stdio.h>
# include <math.h>
# define num 50
# define EPS 0.000001
# define wi 1

  int n;              /*存放所提供函数 含有自变量x的项数*/
  int m;              /*存放所提供数据的 组数 - 1*/
  float x[num];             /*存放给定数据组的x值*/
  float y[num];             /*存放给定数据组的y值*/
  float result_x[num][num]; /*返回xx()的结果*/
  float result_y[num];      /*返回yy()的结果*/
  float result[num];        /*存放高斯求解后的一组结果*/

void main()
{ 
  void xx();
  void yy(float y[]);
  float f(int x,int i);
  void Gauss(float x_x[][num], float y_y[]);
  void last_y(float ax);

  float ax;          /* 求出目标函数后，求y值对应的x值*/
  int i;  

/* x[0]=0; x[1]=0.2; x[2]=0.4; x[3]=0.6; x[4]=0.8; x[5]=1.0; x[6]=1.2; 
 y[0]=0.9; y[1]=1.9; y[2]=2.8; y[3]=3.3; y[4]=4.0; y[5]=5.7; y[6]=6.5;  */
  
  printf("请输入所给函数含有自变量x的项数：n = ");
  scanf("%d",&n);

  printf("\n请输入提供数据的(组数-1)：m = ");
  scanf("%d",&m);

  printf("请输入要求函数值对应的自变量x：x = ");
  scanf("%f",&ax);

  printf("\n请输入提供数据所有x值：\n");
  for(i=0; i<=m; i++)
     scanf("%f",&x[i]);

  printf("\n请输入提供数据所有y值：\n");
  for(i=0; i<=m; i++)
     scanf("%f",&y[i]);

  printf("所有数据输入完成!\n\n");   

  xx();                 /* 得到高斯消去法左边的二维矩阵*/
  yy(y);                /* 得到高斯消去法又边的一维数组*/

  Gauss(result_x,result_y);  /*求得给定函数的系数数组*/
  last_y(ax);           /* 由高斯方程求得的系数解和给定的自变量值求最终结果*/
}


 /*...............公用模块.......................*/
   float f(int h, int i) 
  {  
     int s;
     float d;
     d = 1;