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òËƲ»¿ÉÒÔ°¡£¬Ð»Ð»ÁË