最小二乘法的编程实现..doc

最小二乘法： （用方法计算逆矩阵） #include stdio.h #include math.h #include stdlib.h #includestring.h #include iostream.h #define N 200 #define n 9 void Getdata(double sun[N])//从txt文档中读取数据（小数） { char data; char sunpot[10]={0000000000};//为防止结果出现‘烫’字 int i=0,j=0; double d; FILE *fp=fopen(新建文本文档.txt,r); if(!fp) { printf(cant open file\n); } while(!feof(fp)) { data=fgetc(fp); if(data!=\n) { sunpot[i]=data; i++; } else if(data==\n) { sunpot[i]=\0;//给定结束符 d=atof(sunpot);//将字符串转换成浮点数 sun[j]=d; j++; i=0;//将i复位 } } } void Normal(double sun[N],double sun1[N])//将数据进行标准化 { double mean,temp=0,variance=0; int i; for(i=0;iN;i++) temp=temp+sun[i]; mean=temp/N; for(i=0;iN;i++) { sun1[i]=sun[i]-mean; } for(i=0;iN;i++) { variance=variance+sun1[i]*sun1[i]; } variance=variance/N; variance=sqrt(variance); for(i=0;iN;i++) { sun1[i]=sun[i]/variance;//减去均值除以均方差进行归一化 } } void Matrix(double sun[N],double matrixl[n][n],double matrixr[n])//组建方程的左右两个矩阵 { double matrix1[N-n][n];//方程左边系数矩阵 double matrix_trans[n][N-n]; double matrix2[N-n];//方程右边系数矩阵 int i,j,k; double temp; for(i=0;iN-n;i++)//组建N-nxn维方程矩阵 { for(j=0;jn;j++) { matrix1[i][j]=sun[n+i-j-1];//此处与matlab不同，matlab是从1开始的 } } for(i=0;iN-n;i++) matrix2[i]=sun[i+n]; for(i=0;in;i++)//将N-nxn维矩阵转置 { for(j=0;jN-n;j++) { matrix_trans[i][j]=matrix1[j][i]; } } for(i=0;in;i++)//组建nxn维方阵(A*A) { for(k=0;kn;k++) { temp=0; for(j=0;jN-n;j++) { temp+=matrix_trans[i][j]*matrix1[j][k]; } matrixl[i][k]=temp; } } for(i=0;in;i++) { temp=0; for(j=0;jN-n;j++) { temp+=matrix_trans[i][j]*matrix2[j]; } matrixr[i]=temp; } } void ExchangeLine(double a[n][n], int i, int j, int p)//交换矩阵第i行和第j行 { int k; double temp; for (k=0;k=p-1;k++) { temp=a[j][k]; a[j][k]=a[i][k]; a[i][k]=temp; } } void Gauss(double a[n][n],int k,int p)//高斯消元,用第k行消去其他行第一个非零元素 { int i,j; double m; for(i=