最小二乘法拟合曲线

/*
曲线拟合
int CurveFitter(int x_data_addr, float* y_data);
y_data--------输入--y数据
返回值-------最大值下标
*/
#define rank_ 3
float atemp[2 * (rank_ + 1)] = { 0 };
float b[rank_ + 1] = { 0 };
float a[rank_ + 1][rank_ + 1];
float fit_data[60];
int CurveFitter(float* y_data)
{
		atemp[1] =0;
		atemp[2] =0;
		atemp[3] =0;
		atemp[4] =0;
		atemp[5] =0;
		atemp[6] =0;
		b[0] =0;
		b[1] =0;
		b[2] =0;
		b[3] =0;
    for (int i = 0; i < 60; i++)
    {
        atemp[1] =atemp[1]+i;
        atemp[2] =atemp[2]+i*i;  //atemp[2] =atemp[2]+ pow(i,2);
        atemp[3] =atemp[3]+i*i*i;//atemp[2] =atemp[2]+ pow(i,3);
        atemp[4] =atemp[4]+i*i*i*i;//atemp[2] =atemp[2]+ pow(i,4);
        atemp[5] =atemp[5]+i*i*i*i*i;//atemp[2] =atemp[2]+ pow(i,5);
        atemp[6] =atemp[6]+i*i*i*i*i*i;//atemp[2] =atemp[2]+ pow(i,6);
        b[0] +