递推法c语言程序,最小二乘法,递推最小二乘法C语言源程序

outFile.close;

cout << "Over!" /inFile.open("yk.txt");

fai = new double *[row_y];

for(i = 0; i < row_y; i++)

fai[i] =new double [col_y];

double *pyk = yk;

double *puk = uk;

double *pfai = *fai;

pyk = &yk[1];

puk = &uk[1];

for(int j = 0; j < row_y; j++)

{

pfai = &fai[j][0];

*pfai = -(*(pyk+1));

*(pfai+1) = -(*pyk);

*(pfai+2) = *(puk+1);

*(pfai+3) = *puk;

pyk++;

puk++;

}

}

//iterative least square method

#define rowPn 4

#define colPn 9

void iterative_least_square_method(void)

{

double **Pn,*En_1,*Pn_En_1,*En_Pn_1,**PE,**Pn_1,**Pn_En_0;

double En_Pn_En_1