原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
《数值计算方法》实验报告
专业： 姓名： 学号： 班级： 成绩：
实验名称
实验5 最小二乘拟合法
实验题目
在某化学反应里，测得某物质的浓度 y （单位：％）随时间t
（单位：min）的变化数据如表5—7所列
t
y
t
y
t
y
t
y
/
5
9
13
2
6
10
14
3
7
11
15
4
8
12
16
理论上已知y与t间的关系为
b/t
y ae ,
其中a0和b0为待定系数。上式两端取对数可得In y=ln a+b/t.
做变量替换
z=ln y,x=/t,并记 A=ln a,B=b则有
z=A+Bx.
根据所测数据，利用最小二乘直线拟合法先确定系数 A和B,
进而给出y与t间的关系。
3.实验目的
熟练使用最小二乘拟合法
基础理论
最小二乘拟合法
实验环境
Microsoft Visual C++
实验过程
结果分析
A^l -0567,B=2.427?
y=n^1 ,427ftZt^
Pi*e ss ^Lih y k^：y Lu c； unt lriu@
本次试验令我更加熟悉最小二乘拟法；
附录：程序清单
#in clude
#in clude
void mai n(X
int i=0;
double z[16],x[16],D,a,b;
double t[16] ={1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16；
double y[16]={,,,,,,,,,,,,,,};
double sum_x=0,sum_x2=0,sum_y=0,sum_xy=0; for(i=0;i16;i++){
x[i]=1/t[i];
z[i]=log(y[i]);
}
for(i=0;i16;i++){
sum_x=sum_x+x[i];
sum_x2=sum_x2+x[i]*x[i];
sum_y=sum_y+z[i];
sum_xy=sum_xy+x[i] * z[i];
} 一 一
D=sum_x2* 16-sum_x^sum_x;
a=(16*sum_xy- sum_x^ sum_y)/ D;
b=(sum_x0 sum_y-sum_x* sum_xy)/D;
prin tf(A=%,B=%\n\n,a,b);
prin tf( y=eA%*eA(%t) ,a,b);
}
文档评论（0）