function?[Vint]=IDW（original_xvectororiginal_yvectororiginal_zvectorinterpolated_xvectorinterpolated_yvectordistance_weightmoderadius_lenght）

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%***************************************
%%??INVERSE?DISTANCE?WEIGHT????????????**
%%***************************************
%******INPUTS***********?
%original_xvector?=?stations?x?coordinates?（columns）?[vector]?
%original_yvector?=?stations?y?coordinates?（rows）?[vector]?
%original_zvector?=?variable?values?on?the?point?[original_xvector?original_yvector]??
%interpolated_xvector?=?interpolation?points??interpolated_xvector?coordinates?[vector]?
%interpolated_yvector?=?interpolation?points?interpolated_yvector?coordinates?[vector]?
%distance_weight?=?distance?weight?form?-2/3?to?-5/3?but?not?large?2?（very?important）
%?mode?should?be:?‘fr‘?:?fixed?radius?????or?????‘ng‘?:?neighbours
%radius_lenght?---?radius?lenght?
%%%?OUTPUTS?
%Vint?---?Matrix?[length（interpolated_yvector）length（interpolated_xvector）]?with?interpolated??variable?values?
%
%******************************EXAMPLES?
%%%????V_spa=IDW（x1y1v1xy-2‘ng‘length（x1））;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%???Yang?CS???--?ycsgps@163.com?
%???Copyright?2009?
%???$Date:?2009/08/02?$
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vint=zeros（length（interpolated_yvector）length（interpolated_xvector））;
original_xvector=reshape（original_xvector1length（original_xvector））;
original_yvector=reshape（original_yvector1length（original_yvector））;
original_zvector=reshape（original_zvector1length（original_zvector））;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if??strcmp（mode‘fr‘）
????for?i=1:length（interpolated_xvector）
????????for?j=1:length（interpolated_yvector）
????????????D=[];?V=[];?wV?=[];?vcc=[];
????????????D=?sqrt（（interpolated_xvector（i）-original_xvector）.^2?+（interpolated_yvector（j）-original_yvector）.^2）;
????????????D=D（D????????????V?=?vcc.*（D.^distance_weight）;
????????????wV?=?D.^distance_weight;
????????????if?isempty（D）
????????????????V=NaN;
????????????else
????????????????V=sum（V）/sum（wV）;
????????????end
????????????Vint（ji）=V;
????????end
????end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
else
????if?（radius_lenght?>?length（original_zvector））?||?（radius_lenght<1）?
????????disp（‘Error:?Number?of?neighbours?not?sufficient‘）?
????end?
????for?i=1:length（interpolated_xvector）
????????for?j=1:length（interpolated_yvector）
????????????D=[];?V=[];?wV?=[];vcc=[];
????????????D=?sqrt（（interpolated_xvector（i）-original_xvector）.^2?+（interpolated_yvector（j）-original_yvector）.^2）;
????????????[DI]=sort（D）;
????????????vcc=original_zvector（I）;
????????????V?=?vcc（1:radius_lenght）.*（D（1:radius_lenght）.^distance_weight）;
????????????wV?=?D（1:radius_lenght）.^distance_weight;
????????????V=sum（V）/sum（wV）;
????????????Vint（ji）=V;
?

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????3078??2009-08-02?22:37??IDW.m

-----------?---------??----------?-----??----

?????????????????3078????????????????????1