(完整版)数值分析插值法.docVIP

下载本文档

26
0
约6.75千字
约 10页
2021-10-16 发布于境外
举报
版权申诉

(完整版)数值分析插值法.doc

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

第二章插值法 2•在区间[-1,1]上分别取n=10,20用两组等距节点对龙哥函数 f(x)=1/(1+25*xA2)做多项式插值及三次样条插值，对每个 n值，分别画出插值函数及 f ( x)的图形。 ( 1 )多项式插值 ① 先建立一个多项式插值的 M-file ；输入如下的命令(如牛顿插值公式) ： function [C,D]=newpoly(X,Y) n=length(X); D=zeros(n,n) D(:,1)=Y for j=2:n for k=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); for k=(n-1):-1:1 C=conv(C,poly(X(k))) m=length(C); C(m)= C(m)+D(k,k); end ② 当n=10时，我们在命令窗口中输入以下的命令： clear,clf,hold on; X=-1:0.2:1; Y=1丿(1+25*X.A2); [C,D]=newpoly(X,Y); x=-1:0.01:1; y=polyval(C,x); plot(x,y,X ,Y ■)； grid on; xp=-1:0.2:1; z=1./(1+25*xp.A2); plot(xp,z,r) 得到插值函数和f( x)图形： ③ 当n=20时，我们在命令窗口中输入以下的命令: clear,clf,hold on; X=-1:0.1:1; Y=1./(1+25*X.A2); [C,D]=newpoly(X,Y); x=-1:0.01:1; y=polyval(C,x); plot(x,y,X ,Y .)； grid on; xp=-1:0.1:1; z=1./(1+25*xp.A2); plot(xp,z,r) 得到插值函数和f (x)图形： (2)三次样条插值 ①先建立一个多项式插值的 M-file ；输入如下的命令 : function S=csfit(X,Y,dx0,dxn) N=length(X)-1; H=diff(X); D=diff(Y)./H; A=H(2:N-1); B=2*(H(1:N-1)+H(2:N)); C=H(2:N); U=6*diff(D); B(1)=B(1)-H(1)/2; U(1)=U(1)-3*(D(1)); B(N-1)=B(N-1)-H(N)/2; U(N-1)=U(N-1)-3*(-D(N)); for k=2:N-1 temp=A(k-1)/B(k-1); B(k)=B(k)-temp*C(k-1); U(k)=U(k)-temp*U(k-1); end M(N)=U(N-1)/B(N-1); for k=N-2:-1:1 M(k+1)=(U(k)-C(k)*M(k+2))/B(k); end M(1)=3*(D(1)-dx0)/H(1)-M(2)/ 2; M(N+1)=3*(dxn-D(N))/H(N)-M(N)/ 2; for k=0:N-1 S(k+1,1)=(M(k+2)-M(k+1))/(6*H(k+1)); S(k+1,2)=M(k+1)/2; S(k+1,3)=D(k+1)-H(k+1)*(2*M(k+1)+M(k+2))/ 6; S(k+1,4)=Y(k+1); end ②当 n=10 时，我们在命令窗口中输入以下的命令： clear,clc X=-1:0.2:1; Y=1./(25*X92+1); dx0= 0.0739644970414201;dxn= -0.0739644970414201; S=csfit(X,Y,dx0,dxn) x1=-1:0.01:-0.5;y1=polyval(S(1,:),x1-X(1)); x2=-0.5:0.01:0;y2=polyval(S(2,:),x2-X(2)); x3=0:0.01:0.5; y3=polyval(S(3,:),x3-X(3)); x4=0.5:0.01:1;y4=polyval(S(4,:),x4-X(4)); plot(x1,y1,x2,y2,x3,y3,x4,y4, X,Y,.) 结果如图： ②当n=20时，我们在命令窗口中输入以下的命令： clear,clc X=-1:0.1:1; Y=1./(25*X92+1); dxO= 0.0739644970414201;dx n= -0.0739644970414201; S=csfit(X,YdxO,dx n) x1=-1:0.01:-0.5;y仁polyval(S(1,:),x1-X(1)); x2=-0.5:0.01:0;y2=polyval(S(2,:),x2-X (2)); x3=0:0.01:0.5