原创力文档- 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章
预备：编制计算拉格朗日插值的m文件lagr1.m:
function y=lagr1(x0,y0,x)
n=length(x0);m=length(x);
for i=1:m
z=x(i);
s=0.0;
for k=1:n
p=1.0;
for j=1:n
if j~=k
p=p*(z-x0(j))/(x0(k)-x0(j));
end
end
s=s+p*y0(k);
end
y(i)=s;
end
预备：对数值函数，用辛普森公式计算定积分的程序simp.m:
function s=simp(y,h,m)
s=0;
for k=1:m
s=s+4*y(2*k);
end
for k=1:(m-1)
s=s+2*y(2*k+1);
end
s=(s+y(1)+y(2*m+1))*h/3;
3.1 对函数, 在n 个节点上（ n 不要太大，如5至11）用拉格朗日、分段线性、三次样条三种插值方法，计算m 个插值点的函数值（m 要适中，如50至100）。通过数值和图形输出，将三种插值结果与精确值进行比较。适当增加n ，再作比较，由此作初步分析。
解：
函数插值比较程序如下：
clear;
n=5;%在n个节点上进行插值
m=75;%产生m个插值点，计算函数在插值点处的精确值，将来进行对比
x=-2:4/(m-1):2;
y=exp(-x.^2);
z=0*x;
x0=-2:4/(n-1):2;
y0=exp(-x0.^2);
y1=lagr1(x0,y0,x);% y1为拉格朗日插值
y2=interp1(x0,y0,x);% y2为分段线性插值
y3=spline(x0,y0,x);% y3为三次样条插值
[x y y1 y2 y3]
plot(x,z,k,x,y,r:,x,y1,g-.,x,y2,b,x,y3,y--)
gtext(Lagr.), gtext(Pieces. linear), gtext(Spline),
gtext(y=exp(-x.^2))
hold off;%比较插值所得结果与函数在插值点处的精确值
s = x y y1 y2 y3
[x y y1 y2 y3]
【MATLAB 计算结果】
n=5时，得到结果如下：
数值比较如下：
x y y1 y2 y3
－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－
-2.0000 0.0183 0.0183 0.0183 0.0183
-1.9467 0.0226 -0.0328 0.0370 -0.0082
-1.8933 0.0277 -0.0717 0.0556 -0.0276
-1.8400 0.0339 -0.0990 0.0742 -0.0404
-1.7867 0.0411 -0.1158 0.0929 -0.0468
-1.7333 0.0496 -0.1229 0.1115 -0.0472
-1.6800 0.0595 -0.1211 0.1302 -0.0419
-1.6267 0.0709 -0.1112 0.1488 -0.0314
-1.5733 0.0841 -0.0940 0.1675 -0.0159
-1.5200 0.0992 -0.0702 0.1861 0.0041
-1.4667 0.1164 -0.0406 0.2047 0.0284
-1.4133 0.1357 -0.0058 0.2234 0.0566
-1.3600 0.1573 0.0334 0.2420 0.0883
-1.3067 0.1813 0.0764 0.2607 0.1232
-1.2533 0.2079 0.1226 0.2793 0.1609
-1.2000 0.2369 0.1714 0.2980 0.2011
-1.1467 0.2685 0.2222 0.316
文档评论（0）