二维插值、牛顿插值的matlab程序代码(共4页).doc

《二维插值、牛顿插值的matlab程序代码(共4页).doc》由会员分享，可在线阅读，更多相关《二维插值、牛顿插值的matlab程序代码(共4页).doc（4页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上二维插值的matlab程序代码x=129.0140.0103.588.0185.5195.0105.5157.5107.577.081.0162.0162.0117.5;y=7.5141.523.0147.022.5137.585.560.5121.03.056.5116.584.043.5;plot(x,y,o)Z=48686889988949;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

2、中采用二维插值进行三维绘图之后，如何把图中的各个数据点导出？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); enden

3、dc=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. -0. -0. -0.;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专心-专注-专业