- 1、本文档共6页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
清华杨顶辉数学实验三作业
第 10题
用分段线性插值
x=[0,3:2:9,11:1:15];y1=[0,1.8,2.2,2.7,3.0,3.1,2.9,2.5,2.0,1.6];y2=[0,1.2,1.7,2.0,2.1,2.0,1.8,1.2,1.
0,1.6];
xi=[0:0.1:15];
y1i=interp1(x,y1,xi);y2i=interp1(x,y2,xi);
plot(xi,y1i,xi,y2i)
s1=trapz(xi,y1i);s2=trapz(xi,y2i);s=s1-s2
s=
10.7500
用三次样条插值:
y1i=interp1(x,y1,xi,spline);y2i=interp1(x,y2,xi,spline);
plot(xi,y1i,xi,y2i)
s1=trapz(xi,y1i);s2=trapz(xi,y2i);s=s1-s2
s=
11.3444,用这种方法得到的结果和课后答案一样
【比较】在本题中,三次样条插值的图像比分段线性插值的图像平滑,更符合实际情况。
第11题
三次样条插值:
clear
x=[7.0,10.5,13.0,17.5,34.0,40.5,44.5,48.0,56.0,61.0,68.5,76.5,80.5,91.0,96.0,101.0,104.0,106.5,
111.5,118.0,123.5,136.5,142.0,146.0,150.0,157.0,158.0];
y1=[44,45,47,50,50,38,30,30,34,36,34,41,45,46,43,37,33,28,32,65,55,54,52,50,66,66,68];
y2=[44,59,70,72,93,100,110,110,110,117,118,116,118,118,121,124,121,121,121,122,116,83,81,
82,86,85,68];
xi=[7.0:0.1:158.0];
y1i=spline(x,y1,xi);y2i=spline(x,y2,xi);
plot(xi,y1i,xi,y2i)
s1=trapz(xi,y1i);s2=trapz(xi,y2i);s=s2-s1;
s_fact=s*(40*10^3/(18*10^-3))^2*10^-12
s_fact=
4.2468e+004
分段线性插值:
y1i=interp1(x,y1,xi);y2i=interp1(x,y2,xi);
plot(xi,y1i,xi,y2i)
s1=trapz(xi,y1i);s2=trapz(xi,y2i);s=s2-s1;
s_fact=s*(40*10^3/(18*10^-3))^2*10^-12
s_fact=
4.2414e+004
【比较】本题的三次样条插值效果不好,在[150,160]阶段,形成了一个突出的尖峰,明显偏
离书上的实际图案。而分段线性插值则和书上图案比较贴近。
第 12题
分段线性插值:
clear
x=[0,2,4,5,6,7,8,9,10.5,11.5,12.5,14,16:1:24];
y=[2,2,0,2,5,8,25,12,5,10,12,7,9,28,22,10,9,11,8,9,3];
xi=[0:1/60:24];
yi=interp1(x,y,xi);
plot(xi,yi)
s=fix(sum(yi))
s=
12992
三次样条插值:
yi=spline(x,y,xi);
plot(xi,yi)
s=fix(sum(yi))
s=
12670
文档评论(0)