原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验1拉格朗日插值.ppt
实验1：拉格朗日插值 1．实验原理 2．实验内容 function fi=Lagran_(x,y,xi) %根据给出的数据向量x和y，构 %造拉格朗日插值多项式，并计算在xi处的值。 fi=zeros(size(xi)); % 全零阵 np1=length(y); for i=1:np1 z=ones(size(xi)); %全1阵 for j=1:np1 if i~=j z=z.*(xi-x(j))/(x(i)-x(j)); %插值基函数 end end fi=fi+z*y(i); end return 3、例题：已知， 利用抛物插值计算 的值。 实验2：牛顿插值 经过 个点 ，构造一个n次多项式，形如： 其中 为各阶差商。 K阶差商的计算公式如下： 2．实验内容 function [C,D]=newpoly(x,y) %牛顿插值 n=length(x); D=zeros(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 3. 例题：函数表为 x 0.40 0.55 0.65 0.80 0.90 1.05 f(x) 0.4075 0.57815 0.69675 0.88811 1.02652 1.2538 利用牛顿插值多项式计算f(0.596)的值 3．实验内容 function a=ployfit(x,y,w) %最小二乘拟合 m=length(w); G=zeros(m); D=zeros(m,1); for i=1:m for j=1:m G(i,j)=sum(w.*(x.^(i+j-2))); end D(i,1)=sum(w.*y.*(x.^(i-1))); end G=G(1:2,1:2) D=D(1:2) a=G^(-1) *D * * 经过 个点， ，构造一个n次多项式，形如： 使得 成立。 其中 为插值基函数。 1．实验原理
文档评论（0）