function AA
load inter2.mat;
xi=inter2xy(x);%计算插值节点（中点坐标x）
yi=inter2xy(y);%计算插值节点（中点坐标y）
Fl=TriScatteredInterp(x(:),y(:),z(:),'linear');%matlab 线性散点插值z
zil=Fl(xi,yi);
zi=inter2xy(z);%matlab 计算散点中点值z
mesh(x,y,z)
hold on
plot3(xi,yi,zi,'o')
plot3(xi,yi,zil,'*')

function xx1=inter2xy(x)
[m,n]=size(x);
xx=zeros(2*m-1,2*n-1);
for i=1:m
    for j=1:n
    xx(2*i-1,2*j-1)=x(i,j);
    end
end
xxo=xx;
for i=1:m
    for j=1:n-1
    xx(2*i-1,2*j)=1/2*(xx(2*i-1,2*j-1)+xx(2*i-1,2*j+1));
    end
end
for i=1:m-1
    for j=1:2*n-1
    xx(2*i,j)=1/2*(xx(2*i-1,j)+xx(2*i+1,j));
    end
end
xx1=xx-xxo;

3楼: Originally posted by 411wu at 2011-12-06 23:21:35:
我试了一下，系统提示 function AA Error: Function definitions are not permitted in this context. 2011-12-06 23:20:54貌似不可以啊，谢谢了