原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验报告3插值与拟合)
淮海工学院实验报告
班级 精算111 姓名 张琪 学号 2011122373
指导教师 张滦云 课程名称 数学建模与实验 成绩
实验室 数学建模实验室 实验项目 3、插值与拟合
同组实验者 实验日期 2013-3-22 一、实验目的
掌握用MATLAB计算三种插值（拉格朗日、分段线性、三次样条）的方法，掌握用MATLAB作线性最小二乘拟合的方法，并能运用插值与拟合方法解决实际问题。
二、实验所用软件及版本
MATLAB 7.1
三、实验内容
1、编制计算拉格朗日插值的M文件。
2、根据机床加工问题所提出的加工要求，以表1（P62）中所给的数据为节点，在x = 0到x = 15范围内求步长为0.1的插值（用三种插值方法）。通过数值和图形输出，将插值结果与精确值进行比较。
3、对录像机计数器问题，利用P42表1中的数据，拟合确定模型t = an2+bn中的参数a, b.
四、实验过程记录与结果报告（含基本步骤、主要程序、实验结果、异常情况等）
1. function y=lagr(x0,y0,x)
n=length(x0);m=length(x);
for i=1:m
z=x(i);
s=0;
for k=1:n
p=1;
for j=1:n
if j~=k
p=p*(z-x0(j))/(x0(k)-x0(j));
end
end
s=p*y0(k)+s;
end
y(i)=s;
end
2. x0=[0 3 5 7 9 11 12 13 14 15];
y0=[0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6];
x=0:0.1:15;
y1=lagr(x0,y0,x);
y2=interp1(x0,y0,x);
y3=interp1(x0,y0,x,spline);
[x,y1,y2,y3]
subplot(3,1,1)
plot(x0,y0,k+,x,y1,r)
grid,title(lagrange)
subplot(3,1,2)
plot(x0,y0,k+,x,y2,r)
grid,title(piecewise linear)
subplot(3,1,3)
plot(x0,y0,k+,x,y3,r)
grid,title(spline)
ans =
0 0 0 0
0.1000 -4.9437 0.0400 0.0499
0.2000 -8.8200 0.0800 0.0990
0.3000 -11.7726 0.1200 0.1474
0.4000 -13.9310 0.1600 0.1951
0.5000 -15.4117 0.2000 0.2421
0.6000 -16.3190 0.2400 0.2884
0.7000 -16.7457 0.2800 0.3340
0.8000 -16.7747 0.3200 0.3788
0.9000 -16.4792 0.3600 0.4230
1.0000 -15.9238 0.4000 0.4665
1.1000 -15.1651 0.4400 0.5094
1.2000 -14.2524 0.4800 0.5515
1.3000 -13.2287 0.5200 0.5930
1.4000 -12.1307 0.5600 0.6338
1.5000 -10.9898 0.6000 0.6739
1.6000 -9.8326 0.6400 0.7134
1.7000 -8.6812 0.6800 0.7523
1.8000 -7.5539 0.7200 0.7904
1.9000 -6.4654 0.7600 0.8280
2.0000 -5.4272 0.8000 0.8
文档评论（0）