kriging.m

来自「Kriging 内插算法 通过已知点的测量值来内插获取未知点的预测」· M 代码 · 共 38 行

function z0 = kriging(x,y,z,xi,yi)
%KRIGING Summary of this function goes here
%   Detailed explanation goes here

% STEP 1 Coefficient determination
% calculate the Distance among known points
for ii=1:length(x)
    r(ii,:)=sqrt((x-x(ii)).*(x-x(ii))+(y-y(ii)).*(y-y(ii)));
    
end

% calculate the Kriging Matrix
% form the Kriging Matrix AA
bb=zeros(length(x),1)+1;
AA=[r,bb];
bb=[bb',0];
AA=[AA;bb];

% known measurement array z
z=[z;0];

% coeffieicnt calculation kese
kese=inv(AA)*z;

% STEP 2 Estimation based on known measurement
[kk,ii]=size(xi);[jj,kk]=size(yi);
z0(ii,jj)=0;
for iter_i=1:ii;
    for iter_j=1:jj;
        xxi=xi(1,iter_i);yyi=yi(iter_j,1);
        rij=sqrt((x-xxi).*(x-xxi)+(y-yyi).*(y-yyi));
        rrij=[rij;1];
        z0(iter_i,iter_j)=rrij'*kese;
    end
end

end