function?[g]=check（Xa1c1c01a2c2c02）
%?交叉驗(yàn)證
%?輸入:?X、Y方向的實(shí)驗(yàn)變差函數(shù)參數(shù)a1c1c01a2c2c02，原始數(shù)據(jù)X（N3）
%?輸出:?數(shù)組?g（）———g（:1）___x（）
%??????????????????????g（:2）___y（）
%??????????????????????g（:3）___估計(jì)值
%??????????????????????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））;??????????%?存放搜索到的最近的六個(gè)點(diǎn)
????????liu（k2）=y（bj（k））;
????????liu（k3）=z（bj（k））;
????end
??
????%克里格估值??矩陣?K
????for?d=1:6?%?設(shè)計(jì)可以由用戶(hù)更改
????????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;%??更改參與插值個(gè)數(shù)后記得更改

????%矩陣?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;
????%?計(jì)算點(diǎn)估值方差
????for?b=1:7
????????sum1（b）=lm（b1）.*m（b1）;
????end
????????sum2=sum（sum1）;
????????g（i5）=c+c0-sum2;
????%?估計(jì)偏差
????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‘）







????
????
????

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

?????文件???????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個(gè)文件信息