网站大量收购闲置独家精品文档,联系QQ:2885784924

上海交通大学计算方法作业答案.pdfVIP

  1. 1、本文档共6页,可阅读全部内容。
  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文档。上传文档
查看更多

P50–1

%%牛顿插值多项式

function[c,d]=newpoly(x,y)

%这里x为n个节点的横坐标所组成的向量,y为纵坐标所组成的向量。

%c为所求的牛顿插值多项式的系数构成的向量。

n=length(x);

d=zeros(n,n);

d(:,1)=y;

forj=2:n

fork=j:n

d(k,j)=(d(k,j-1)-d(k-1,j-1))/(x(k)-x(k-j+1));

end

end

c=d(n,n);

fork=(n-1):-1:1

c=conv(c,poly(x(k)));

m=length(c);

c(m)=c(m)+d(k,k);

end

x=0.2:0.2:1;

y=[0.98,0.92,0.81,0.64,0.38];

c=newpoly(x,y)

c=-0.52080.8333-1.10420.19170.9800

%%三次样条插值

x=[0.2,0.4,0.6,0.8,1.0];

y=[0.98,0.92,0.81,0.64,0.38];

x0=[0.2,0.28,1.0,1.08];

pp=csape(x,y,variational);

%%三次样条函数表达式

disp(pp.coefs);

-1.3393-0.0000-0.24640.9800

0.4464-0.8036-0.40710.9200

-1.6964-0.5357-0.67500.8100

2.5893-1.5536-1.09290.6400

绘制曲线图

x2=0:0.01:1.2;

y11=polyval(c,x2);

y22=ppval(pp,x2);

x0=[0.2,0.28,1.0,1.08];

y110=polyval(c,x0);

y220=ppval(pp,x0);

plot(x2,y11,r,x0,y110,^,x2,y22,g,x0,y220,h)

legend(牛顿插值,牛顿插值样点,三次样条插值,三次样条插值样点)

P50-3

(1)x=[0,1,4,9,16,25,36,49,64];

y=0:8;

x1=0:0.1:64;

x2=0:0.01:1;

f=lagrange(x,y)

%%得到多项式函数表达式

L(x)=-3.28063e-10*x^8+6.71268e-8*x^7-0.00000542921*x^6+0.000222972*x^5-

0x^4+0.0604294*x^3-0.38141*x^2+1.32574*x

y1=lagrange(x,y,x1);

y2=lagrange(x,y,x2);

(2)(2)x=[0,1,4,9,16,25,36,49,64];

y=0:8;

x1=0:0.1:64;

x2=0:0.01:1;

%%得到三次样条差值函数表达式

pp=csape(x,y,not-a-knot);

disp(pp.coefs);

0.0266-0.29981.27320

0.0266-0.21990.75341.0000

-0.00210.01970.15292.0000

0.0005-0.01120.19553.0000

-0.0000-0.00010.11604.0000

0.0000-0.00140.10265.0000

文档评论(0)

172****7268 + 关注
实名认证
文档贡献者

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

1亿VIP精品文档

相关文档