matlab中polyfit(x,y,2)函数的c语言的源代码

/**********************************************
*Author :wacs5
*DATE :20090408(YYYMMDD)
*Functtion :多项式拟合polyfit
**********************************************/
#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
#include <math.h>
main()
{
int i,j,m,n=7,poly_n=2;
double x[7]={1,2,3,4,6,7,8},y[7]={2,3,6,7,5,3,2};
double a[3];
void polyfit(int n,double *x,double *y,int poly_n,double a[]);

system("cls");
polyfit(n,x,y,poly_n,a);

for (i=0;i<poly_n+1;i++)/*这里是升序排列，Matlab是降序排列*/
printf("a[%d]=%g\n",i,a[i]);
getch();
}


/*==================polyfit(n,x,y,poly_n,a)===================*/
/*=======拟合y=a0+a1*x+a2*x^2+……+apoly_n*x^poly_n========*/
/*=====n是数据个数 xy是数据值 poly_n是多项式的项数======*/
/*===返回a0,a1,a2,……a[poly_n]，系数比项数多一（常数项）=====*/
void polyfit(int n,double x[],double y[],int poly_n,double a[])
{