solve函数c语言_polyfit函数的C语言兑现

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

{

x[ix]=ix+1;

}

double **y=NULL;

int iypoly=0;

y=new double*[nwin_length*npoly];

for (int ifc=0;ifc

{

y[ifc]=new double[nwin_length];

}

for (int ifcx=0;ifcx

{

for (int ifcy=0;ifcy

{

y[ifcx][ifcy]=ypoly[iypoly];

iypoly++;

}

}

double apoly[3];

for (int ixy=0;ixy

{