matlab反距离权重插值

用于matlab的反距离权重插值函数，是比较方便的插值方法之一

function [g]=check（Xa1c1c01a2c2c02）
% 交叉验证
% 输入: X、Y方向的实验变差函数参数a1c1c01a2c2c02，原始数据X（N3）
% 输出: 数组 g（）———g（:1）___x（）
%                      g（:2）___y（）
%                      g（:3）___估计值
%                      g（:4）___采样值
%                      g（:5）___方差
%                      g（:6）___偏差
c=（c1+c2）./2;
c0=（c01+c02）./2;
if a1>a2
    kk=a1./a2;
else   % a1    q=a1;
    a1=a2;
    a2=q;
    kk=a1./a2
end
a=a1;
[ntp]=size（X）;
for i=1:n
    x（i）=X（i1）;
    y（i）=X（i2）;
    z（i）=X（i3）;
    g（i1）=X（i1）;
    g（i2）=X（i2）;
    g（i4）=X（i3）;
end
for i=1:n   
    t=0;
    for j=1:n     
        if i~=j
            t=t+1;
         dx=abs（X（j1）-X（i1））;
         dy=abs（X（j2）-X（i2））;
         dxy（t）=sqrt（dx.*dx+dy.*dy.*kk.*kk）;
        end
    end

    [dijbj]=sort（dxy）;
    for k=1:6%可以输入
       %di（k1）=dij（k）;
        liu（k1）=x（bj（k））;          % 存放搜索到的最近的六个点
        liu（k2）=y（bj（k））;
        liu（k3）=z（bj（k））;
    end
  
    %克里格估值  矩阵 K
    for d=1:6 % 设计可以由用户更改
        for b=1:6
            hx=abs（liu（d1）-liu（b1））;
            hy=abs（liu（d2）-liu（b2））;
            hxy=sqrt（hx.*hx+hy.*hy.*kk.*kk）;
               if hxy==0
                   rhh=0;
               elseif hxy>a
                   rhh=c0+c;
               else
                   rhh=c0+c.*（3.*hxy./（2.*a））-（hxy.*hxy.*hxy）./（2.*a.*a.*a）;
               end        
            k（db）=c0+c-rhh;
            if d==6
                k（d+1b）=1;
            end
        end
            k（db+1）=1;
    end
    k（77）=0;%  更改参与插值个数后记得更改

    %矩阵 M
    for d=1:6
        hx=abs（liu（d1）-X（i1））;
        hy=abs（liu（d2）-X（i2））;
        hxy=sqrt（hx.*hx+hy.*hy.*kk.*kk）;
           if hxy==0
              rhh=0;
           elseif hxy>a
               rhh=c0+c;
           else
               rhh=c0+c.*（3.*hxy./（2.*a））-（hxy.*hxy.*hxy）./（2.*a.*a.*a）;
           end  
     
        m（d1）=c0+c-rhh;
        if d==6
            m（d+11）=1;
        end
    end

    k1=inv（k）;
    lm=k1*m;

    for d=1:6
        lm1（d）=lm（d1）;
        lm2（d1）=liu（d3）;
    end

    g（i3）=lm1*lm2;
    % 计算点估值方差
    for b=1:7
        sum1（b）=lm（b1）.*m（b1）;
    end
        sum2=sum（sum1）;
        g（i5）=c+c0-sum2;
    % 估计偏差
    g（i6）=g（i4）-g（i3）;       % 偏差
end
me=mean（g（i6））;
for i=1:n
    ct（i）=（g（i6）-me）.*（g（i6）-me）;
end
sume=sum（ct）;
cte=sqrt（sume./（n-1））;
cte      %方差
me       %平均值 

%for i=1:n
 %   P（i1）=X（i1）;
  %  P（i2）=X（i2）;
   % P（i3）=g（i2）;
%end
dlmwrite（‘outdatakrigcheck.txt‘ g ‘delimiter‘‘ ‘‘precision‘ ‘%10.2f‘）







    
    
    

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2610  2010-06-15 09:07  mkriging\check.asv

     文件       2808  2010-06-15 09:10  mkriging\check.m

     文件     124020  2010-06-14 22:17  mkriging\data.txt

     文件      50938  2010-06-03 14:11  mkriging\data1.txt

     文件         49  2010-06-14 11:31  mkriging\duishu.m

     文件       1157  2010-06-12 14:08  mkriging\fan.m

     文件       1178  2010-06-12 14:22  mkriging\fan3.m

     文件        988  2010-06-14 23:06  mkriging\fancheck.asv

     文件       1191  2010-06-15 09:27  mkriging\fancheck.m

     文件    2074956  2010-06-08 09:55  mkriging\fin.txt

     文件       3423  2010-06-15 09:17  mkriging\guzhi.asv

     文件       2905  2010-06-15 09:19  mkriging\guzhi.m

     文件        484  2010-06-12 10:25  mkriging\krig3.txt

     文件         66  2010-06-12 10:25  mkriging\krig5.txt

     文件      13992  2010-06-12 10:25  mkriging\krig7.txt

     文件      72600  2010-06-12 10:25  mkriging\krig8.txt

     文件     157410  2010-06-14 22:27  mkriging\outdatafan.txt

     文件     262350  2010-06-14 23:09  mkriging\outdatafancheck.txt

     文件         66  2010-06-14 22:38  mkriging\outdatakrigcanshu.txt

     文件     314820  2010-06-14 22:38  mkriging\outdatakrigcheck.txt

     文件     783552  2010-06-14 22:39  mkriging\outdatakrigguzhi.txt

     文件       2145  2010-06-07 23:57  mkriging\output10.txt

     文件      96525  2010-06-06 15:29  mkriging\output3.txt

     文件     128700  2010-06-11 12:43  mkriging\output9.txt

     文件       3872  2010-06-13 21:42  mkriging\sdx.txt

     文件       3872  2010-06-13 23:13  mkriging\sdy.txt

     文件       1298  2010-06-14 23:51  mkriging\wang.asv

     文件       1293  2010-06-14 23:51  mkriging\wang.m

     文件       3880  2010-06-11 12:42  mkriging\xfit.asv

     文件       3822  2010-06-15 09:02  mkriging\xfit.m

............此处省略44个文件信息