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

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

牛顿插值matlab程序代码
clc;
formatlong;%显示15位
%输入初始数据
x0=[0.40.50.70.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记录差商表的一列
fori=1:n-1%构造差商表%%这里n为节点个数
dx0=dx;
forj=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