最小二乘法曲线拟合c.docVIP

四．C++源代码#includeiostream#includemath.h#include iomanipusing namespace std;double max(double array[4]);double gt(double x[11],double y[11]);double GS(double a[4][4],double b[4],double e[],double s[]);void main(){ double a[11]={1,1,1,1,1,1,1,1,1,1,1},b[11],c[11],d[11],e[11],f[11],g[11]; cout请输入变量t值:; for(int i=0;i11;i++) { cinb[i]; c[i]=b[i]*b[i]; d[i]=b[i]*b[i]*b[i]; } cout请输入y值:; for(i=0;i11;i++) cine[i]; f[0]=gt(a,a);f[1]=gt(a,b);f[2]=gt(a,c);f[3]=gt(a,d);f[4]=gt(b,b);f[5]=gt(b,c);f[6]=gt(b,d); f[7]=gt(c,c);f[8]=gt(c,d);f[9]=gt(d,d);g[0]=gt(a,e);g[1]=gt(b,e);g[2]=gt(c,e);g[3]=gt(d,e); double array[4][4]={{f[0],f[1],f[2],f[3]},{f[1],f[4],f[5],f[6]},{f[2],f[5],f[7],f[8]},{f[3],f[6],f[8],f[9]}}; GS(array,g,e,b);}double gt(double x[11],double y[11]){ double sum=0; for(int i=0;i11;i++) sum=x[i]*y[i]+sum; return sum;}double GS(double a[4][4],double b[4],double e[11],double s[11]){ double y[11]={0}; double c[4]={0}; double x0[4]={0}; int i,k,j; double x[4]; double r,sum=0; double w; for(k=1;;k++) { for(i=0;i4;i++) { for(j=0;j4;j++) { sum=a[i][j]*x0[j]+sum; } x[i]=x0[i]+(b[i]-sum)/a[i][i]; c[i]=fabs(x[i]-x0[i]); x0[i]=x[i]; sum=0; } r=max(c); if(r0 { cout解得的三次拟合曲线系数为:\n; for(i=0;i4;i++) coutAi=x[i]endl; break; } } for(i=0;i11;i++) {y[i]=x[0]+x[1]*s[i]+x[2]*s[i]*s[i]+x[3]*s[i]*s[i]*s[i]; w=fabs(y[i]-e[i]); couty[i] e[i]endl; coutYi的误差为wsetprecision(15)endl;} return 0;}double max(double array[4]) { double a=array[0]; int i; for(i=0;i4;i++) {f(aarray[i]) a=array[i];} return a;}五．流程图输入t与y的值（在程序中以数组形式保存）输入t与y的值（在程序中以数组形式保存）利用权为1的向量进行正则三界方程组的构建。利用GS迭代法计算出拟合曲线系数A0，A1,A2,A3将系数代入拟合函数中求出各个t值处的函数值，并与测量值进行误差比较六．实验结果七．拟合曲线曲线拟合的最小二乘法一．问题提出从一组数据中找出其规律性，给出其数学模型的近似表达式问题。通常可以利用最小二乘法解决这样的问题。在某冶金的过程中，通过实验检测的含碳量与时间的关系的数据如下，试求含量y与时间t内在的关系拟合曲线。t 0 5 10 15 20 25 30 35 40 45 50y0 1.27 2.16 2.86 3.44