C++ 线性拟合程序

#include<iostream>
#include<cmath>
#include"nihe1.h"
using namespace std;
float linefit(Point l_point[],int n_point,float *R,float *a1,float *a0)  //通常的函数（非指针型）调用结束后只能返回一个变量或者对象
                          //要想函数结束时把大量数据返回给主调函数，可以使用指针型函数。
{
float av_x,av_y;
float L_xx,L_yy,L_xy;
float a,b;
    av_x=0;
av_y=0;
L_xx=0;
L_yy=0;
L_xy=0;
for (int i=0;i<n_point;i++)
{
      av_x+=l_point[i].X/n_point;
  av_y+=l_point[i].Y/n_point;
}
for(i=0;i<n_point;i++)
{
      L_xx+=(l_point[i].X-av_x)*(l_point[i].X-av_x);
      L_yy+=(l_point[i].Y-av_y)*(l_point[i].Y-av_y);
      L_xy+=(l_point[i].X-av_x)*(l_point[i].Y-av_y);
    }
a=L_xy/L_xx;
b=av_y-L_xy*av_x/L_xx;
//cout<<"This line can be fitted by y=ax+b."<<endl;
//cout<<"a="<<a<<endl;