数值分析第五版计算实习题.pdfVIP

  1. 1、本文档共8页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多

数值分析第五版计算实习题--第1页

.

数值分析计算实习题

第二章

2-1

程序:

clear;clc;

x1=[0.20.40.60.81.0];

y1=[0.980.920.810.640.38];

n=length(y1);

c=y1(:);

forj=2:n%求差商

fori=n:-1:j

c(i)=(c(i)-c(i-1))/(x1(i)-x1(i-j+1));

end

end

symsxdfd;

df(1)=1;d(1)=y1(1);

fori=2:n%求牛顿差值多项式

df(i)=df(i-1)*(x-x1(i-1));

d(i)=c(i)*df(i);

end

disp(4次牛顿插值多项式);

P4=vpa(collect((sum(d))),5)%P4即为4次牛顿插值多项式,并保留小数点后5位数

pp=csape(x1,y1,variational);%调用三次样条函数

q=pp.coefs;

disp(三次样条函数);

fori=1:4

S=q(i,:)*[(x-x1(i))^3;(x-x1(i))^2;(x-x1(i));1];

S=vpa(collect(S),5)

end

x2=0.2:0.08:1.08;

dot=[121112];

figure

ezplot(P4,[0.2,1.08]);

holdon

y2=fnval(pp,x2);

x=x2(dot);

y3=eval(P4);

y4=fnval(pp,x2(dot));

plot(x2,y2,r,x2(dot),y3,b*,x2(dot),y4,co);

title(4次牛顿插值与三次样条);

结果如下:

4次牛顿插值多项式

P4=-0.52083*x^4+0.83333*x^3-1.1042*x^2+0.19167*x+0.98

三次样条函数

..

数值分析第五版计算实习题--第1页

数值分析第五版计算实习题--第2页

.

x∈[0.2,0.4]时,S=-1.3393*x^3+0.80357*x^2-0.40714*x+1.04

x∈[0.4,0.6]时,S=0.44643*x^3-1.3393*x^2+0.45*x+0.92571

x∈[0.6,0.8]时,S=-1.6964*x^3+2.5179*x^2-1.8643*x+1.3886

x∈[0.8,1.0]时,S=2.5893*x^3-7.7679*x^2+6.3643*x-0.80571

输出图如下

文档评论(0)

180****0535 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档