《MannKendall法(非参数检验方法)程序》由会员分享，可在线阅读，更多相关《MannKendall法(非参数检验方法)程序（2页珍藏版）》请在装配图网上搜索。
1、Mann-Kendall 法（非参数检验方法）程序统计知识 2009-02-0213:47:51 阅读 986 评论 1 字号：大中小订阅Mann-Kendall 法(非参数检验方法)用于气候突变检测 programmainimplicitnonecThisisaprogramfortestingclimatejumpibyuseofcMann-Kendalltest.cinteger,parameter:iy=100integerirealx(iy),u1(iy),u2(iy)cReadDatacopen(31,file=d:mk.txt,form=formatted)doi=1,50 x(
2、i)=0.5x(i+50)=-0.5enddocMann-KendalltestmethodcallMKtest(iy,x,u1,u2)cdoi=1,iywrite(31,10)i,x(i),u1(i),u2(i)enddo10format(1x,i4,4f8.2)stopend!csubroutineMKtest(m,x,u1,u2)dimensionx(m),u1(m),u2(m)dimensionxr(m),ur(m)!workarraydoi=1,mxr(i)=x(m+1-i)enddocallMKrank(m,x,u1)callMKrank(m,xr,ur)doi=1,mu2(i)=-ur(m+1-i)enddoreturnend!cMann-KendallRankStatisticsubroutineMKrank(m,x,u)dimensionx(m),u(m)u(1)=0.0d=0.0doi=2,mnum=0doj=1,iif(x(i).gt.x(j)thennum=num+1endifenddod=d+float(num)e=i*(i-1.)/4.var=i*(i-1.)*(2.*i+5.)/72.u(i)=(d-e)/sqrt(var)enddoreturnend
MannKendall法(非参数检验方法)程序
