二维插值、牛顿插值的matlab程序代码.docVIP

下载本文档

31
0
约1.19千字
约 4页
2018-02-24 发布于河北
举报
版权申诉

二维插值、牛顿插值的matlab程序代码.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文档。上传文档

二维插值、牛顿插值的matlab程序代码.doc

二维插值的matlab程序代码 ?x=[129.0?140.0?103.5?88.0?185.5?195.0?105.5?157.5?107.5?77.0?81.0?162.0?162.0?117.5]; ?y=[7.5?141.5?23.0?147.0?22.5?137.5?85.5?60.5?121.0?3.0?56.5?116.5?84.0?43.5]; ?plot(x,y,o) ?Z=[4?8?6?8?6?8?8?9?9?8?8?9?4?9]; ?h=-Z; ?xi=70:3:200;yi=0:3:150; ?[X,Y]=meshgrid(xi,yi); ?H=griddata(x,y,h,X,Y,cubic); ?mesh(X,Y,H); ?view(-60,30); 图形如下: matlab中采用二维插值进行三维绘图之后，如何把图中的各个数据点导出？[X,Y] = meshgrid(-3:.25:3); Z = peaks(X,Y); [XI,YI] = meshgrid(-3:.125:3); ZI = interp2(X,Y,Z,XI,YI); mesh(X,Y,Z) mesh(XI,YI,ZI)%这句话就导出了插入的点function [c,d]=newpoly(x,y) n=length(x); 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 牛顿插值matlab程序代码 clc; format long; %显示15位 %输入初始数据 x0=[0.4 0.5 0.7 0.8]; y0=[-0.916291 -0.693147 -0.356675 -0.223144]; x=0.6;%插值点 n=max(size(x0)); y=y0(1);%迭代初始值 disp(y); s=1; dx=y0;?? %每次循环用dx记录差商表的一列 for i=1:n-1%构造差商表? %%这里n为节点个数 ??? dx0=dx; ??? for j=1:n-i ??????? dx(j)=(dx0(j+1)-dx0(j))/(x0(i+j)-x0(j)); ??? end ??? df=dx(1); ??? s=s*(x-x0(i)); ??? y=y+s*df;%计算 ??? disp(y); end