计算方法最小二乘拟合c语言程序.docVIP

计算方法最小二乘拟合c语言程序。 运行截图： 程序代码： #include stdio.h #include math.h int N,n,type; //N为数据个数，type为所拟合曲线的代号（详见程序运行截图） void Prin(); //输入提示 void ATA(double a[N][n]); //通过正则方程生成所拟合函数系数的线性方程组 void Solve(double a[n][n+1]); //列主元素消去法求解方程组 void prin1(double a[n][n+1],double b[N][2]); //n次多项式曲线输出函数 void prin2(double a[n][n+1],double b[N][2]); //指数型曲线输出函数 void prin3(double a[n][n+1],double b[N][2]); //对数型曲线输出函数 void prin4(double a[n][n+1],double b[N][2]); //幂函数型曲线输出函数 int main() { Prin(); int i,j; scanf(%d,N); double b[N][2]; for (i=0;iN;i++) scanf(%lf %lf,b[i][0],b[i][1]); while (1) { printf(\n请输入要拟合曲线的代号\n); scanf(%d,type); if (type!=101type!=102type!=103) n=type+1; else n=2; if (typeN) { double a[N][n+1]; for (i=0;iN;i++) a[i][1]=b[i][0],a[i][n]=b[i][1]; for (i=0;iN;i++) a[i][0]=1; for (j=1;jn;j++) for (i=0;iN;i++) a[i][j]=pow(a[i][1],j); ATA(a); Solve(a); prin1(a,b); } if (type==101) //zhishu { double a[N][n+1]; for (i=0;iN;i++) a[i][0]=1,a[i][1]=b[i][0],a[i][n]=log(b[i][1]); ATA(a); Solve(a); prin2(a,b); } if (type==102) { double a[N][n+1]; for (i=0;iN;i++) a[i][0]=1,a[i][1]=log(b[i][0]),a[i][n]=b[i][1]; ATA(a); Solve(a); prin3(a,b); } if (type==103) { double a[N][n+1]; for (i=0;iN;i++) a[i][0]=1,a[i][1]=log(b[i][0]),a[i][n]=log(b[i][1]); ATA(a); Solve(a); prin4(a,b); } } return 0; } void Prin() { printf(请在第一行输入数据对个数n，以下n行输入数据，); printf(最后一行输入要拟合的函数类型代号\n\n); printf((提示：\nn次多项式的函数类型代号为n\n); printf(指数型函数y=A*e^(Bx)的函数类型代号为101\n); printf(对数型函数y=A+B*lnx的函数类型代号为102\n); printf(幂函数型y=A*x^B的函数代号为103)\n\n); } void ATA(double a[N][n+1]) { int i,j,k; double b[n][n+1]; for (i=0;in;i++) { for (j=0;jn+1;j++) { double s=0; for (k=0;kN;k++) s+=a[k][i]*a[k][j]; b[i][j]=s; } } for (i=0;in;i++) for (j=0;jn+1;j++) a[i][j]=b[i][j]; } void Solve(double a[n][n+1]) { int i,j,k; for (i=0;in-1;i++) { double max=