实验4最小二乘拟合法.docVIP

实验名称 实验4 最小二乘拟合法 实验题目：在某化学反应里，测得某物的浓度y（单位：%）随时间t(单位：min)的变化数据如下： t y t y t y t y 1 4.00 5 9.22 9 10.00 13 10.50 2 6.40 6 9.50 10 10.20 14 10.55 3 8.00 7 9.70 11 10.32 15 10.58 4 8.80 8 9.86 12 10.42 16 10.60 理论上已知y与r间的关系为 y=aeb/t, 其中a0和b0为待定系数。上式两端取对数可得ln y=ln a+b/t.做变量替换z=ln y,x=1/t,并记A=ln a,B=b,则有 z=A+Bx. 根据所测数据，利用最小二乘直线拟合法先确定系数A和B，进而给出y与t的关系。 实验目的：掌握最小二乘法，?熟练掌握求函数的近似最小二乘拟合多项式. 基础理论：在科学实验数据处理中，往往要根据一组给定的实验数据求出自变量x与因变量y的函数关系，这是ai为待定参数, 由于观测数据总有误差，且待定参数ai的数量比给定数据点的数量少(即n＜m)，因此它不同于插值问题.这类问题不要求通过点，而只要求在给定点xi上的误差,的平方和最小。 。这就是最小二乘逼近，得到的拟合曲线为y=s(x)，这种方法称为曲线拟合的最小二乘法。 实验环境：操作系统：Windows XP； 实验平台：Visual C++ 6.0 实验过程： 附录： #includestdio.h int main() { double a0,a1,sumTi=0,sumTi2=0,sumRi=0,TiRi=0; double T[100],R[100]; int n,m,i,j; printf(输入数据的个数n:); scanf(%d,n); printf(输入各组数据T[i]与R[i]:\n); m=n; i=0;j=0; while(n--) scanf(%lf%lf,T[i],R[i++]); printf(%d\n,i); for(j=0;ji;j++) { sumTi+=T[j]; sumRi+=R[j]; sumTi2+=T[j]*T[j]; TiRi+=T[j]*R[j]; } a0=-(sumTi*TiRi-sumTi2*sumRi)/(m*sumTi2-sumTi*sumTi); a1=(m*TiRi-sumTi*sumRi)/(m*sumTi2-sumTi*sumTi); printf(a0=%.3lf a1=%.3lf\n,a0,a1); printf(最小二乘拟合曲线为:R=%.3lf+%.3lfT\n,a0,a1); return 0; }