www.gusucode.com > 地统计学中的克里格插值法matlab源码程序 > devico.m
function [Cs]=devico(key) keybar=mean(key); s=0; for i=1:length(key) s=s+(key(i)-keybar)^3*length(find(key==key(i))); end Cs=s/(length(key)*std(key)^3*keybar);