polyfit函数的C语言实现

void PolyfitCf(int n_poly,int Nwin_length,int Npoly,double * ypoly,double **fitcoef)
{
 int i,j,m;
 int nwin_length=Nwin_length;
 int poly_n=n_poly;
 int npoly=Npoly;

 double *x=NULL;//[nwin_length];
 if (x==NULL)
 {
  x=new double[nwin_length];
 }

 for (int ix=0;ix<nwin_length;ix++)
 {
  x[ix]=ix+1;
 }

 double **y=NULL;
 int iypoly=0;
 y=new double*[nwin_length*npoly];
 for (int ifc=0;ifc<npoly;ifc++)
 {
  y[ifc]=new double[nwin_length];
 }
 for (int ifcx=0;ifcx<npoly;ifcx++)
 {
  for (int ifcy=0;ifcy<nwin_length;ifcy++)
  {
   y[ifcx][ifcy]=ypoly[iypoly];
   iypoly++;
  }
 }

 double apoly[3];
 for (int