%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