《最小二乘拟合法》由会员分享,可在线阅读,更多相关《最小二乘拟合法(3页珍藏版)》请在金锄头文库上搜索。
1、最小二乘拟合法Lagrange 插值法与最小二乘拟合法2.实验题目在某化学反应里,测得某物质的浓度y (单位:)随时间t (单位:min) 的变化数据如表 5-7所列。表 5-7 变化数据tytytyty1 4.002 6.403 8.004 8.805 9.226 9.507 9.708 9.869 10.0010 10.2011 10.3212 10.4213 10.5014 10.5515 10.5816 10.60理论上已知y与t的关系为y=aeb/t,其中a0和bvO为待定系数。上式两端取对数可得lny=lna+b/t。做变量替换z=lny, x=l/t,并记 A=lna, B=b,
2、则有z=A+Bx。根据所测数据,利用最小二乘直线拟合法先确定系数A和B,进而给出y与t 间的关系。3. 实验目的验证最小二乘拟合法对数据的拟合程度,掌握用C+程序来编程实现最小二 乘拟合法的算法并对其加深理解。4. 基础理论5 实验环境Visual C+ 语言6.实验过程#include#includeint main()int n=l6;double t,y,xl6,zl6;for(int j=0;jn;j+)cout当 t=j+1y;xj=1/(j+1);zj=log(y);doubledouble D,a,b;for(int i=0;in;i+)Sum_x2=Sum_x2+xi*xi;Sum_z=Sum_z+zi;Sum_xz=Sum_xz+xi*zi;coutendl;D=Sum_x2*n-Sum_x*Sum_x;b=(n*Sum_xz-Sum_x*Sum_z)/D; a=(Sum_x2*Sum_z-Sum_x*Sum_xz)/D a=exp(a);couty=aeA(b/t)endl;return 0;