IDW反距离加权插值法，源代码

IDW的源代码哦！ 要的速度下！ 还带了讲解！

%IDW（反距离加权插值法）
%其中xyz为已知坐标及其函数值XY为要插值的坐标
%xyzXY最高为二维的，不可为三维
%不考虑x，y中出现重复坐标的情况
function [Z]=IDW（xyzXY）
[m0n0]=size（x）;
[m1n1]=size（X）;
%生成距离矩阵r（m0*m1*n1n0）
for i=1:m1
    for j=1:n1
        r（m0*n1*（i-1）+m0*（j-1）+1:m0*n1*（i-1）+m0*（j）:）=sqrt（（X（ij）-x）.^2+（Y（ij）-y）.^2）;
    end
end
%定义插值函数
for i=1:m1
    for j=1:n1
        if find（r（m0*n1*（i-1）+m0*（j-1）+1:m0*n1*（i-1）+m0*（j）:）==0）
            [m2n2]=find（r（m0*n1*（i-1）+m0*（j-1）+1:m0*n1*（i-1）+m0*（j）:）==0）;
            Z（ij）=z（m2n2）;
        else
            numerator=sum（sum（z./r（m0*n1*（i-1）+m0*（j-1）+1:m0*n1*（i-1）+m0*（j）:）））;
            denominator=sum（sum（1./r（m0*n1*（i-1）+m0*（j-1）+1:m0*n1*（i-1）+m0*（j）:）））;
            Z（ij）=numerator/denominator;
        end
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         843  2010-06-10 21:47  IDW.m
     文件      161521  2009-06-29 14:51  基于插值的雨量预报评价模型.pdf
     文件        5171  2010-06-10 21:02  空间插值方法汇总.txt