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
评论0