mann-kendall檢驗(yàn)的程序

1、mannkendall檢驗(yàn)的程序分類:書籍軟件2007.1.28 12:07作者：步行者|評(píng)論：0 |閱讀：5233c.oThis is a program for testing climate jumpc Mann-Kendall testcparameter(iy=116,ihou = 1,is=1)dimension x1(500,ihou),x2(500,ihou),x3(500,ihou),iyear(500)dimension x4(iy),x5(iy),x6(iy)open(10,file=rumei.40,status=old)open(20,file=mk1.txt,sta

2、tus=unknown)cread(10,*)do 100 mxx=1,isdo m1 = 1,iyread(10,188)(x1(m1,m2),m2 = 1,ihou)iyear(m1) = m1 + 1884end doc188format(17x,3f6.0)188 format(f5.0)do m1 = 1,ihoudo m2 = 1,iyx4(m2)=x1(m2,m1)end docall MKtest(x4,iy,x5,x6)do m2 = 1,iyx2(m2,m1)=x5(m2)x3(m2,m1)=x6(m2)end doend dodo m1 = 1,ihoudo m2 = 1

3、,iywrite(20,288)iyear(m2),m1,x2(m2,m1),x3(m2,m1)cwrite(*,*)iyear(m2),m1,x2(m2,m1) enddoend do100continue288format(i4,i2,3f6.2)close(10)close(20)stopendsubroutine mktest(x1,m,u1,u3)dimension x1(m),d1(m),e1(m),var1(m)dimension u1(m),u2(m),d2(m),xp(m),u3(m)d1(1) = 0.0d2(1) = 0.0u1(1) = 0.0u2(1) = 0.0do

4、 m1 = 1,mxp(m1)=x1(m-m1 + 1)end dodo k1 = 2,mkk=0do k2 = 1,k1if(x1(k1).gt.x1(k2) thenkk=kk+1elseend ifend dod1(k1)=d1(k1-1)+float(kk)e1(k1)=float(k1*(k1-1)/4.var1(k1)=float(k1*(k1-1)*(2*k1+5)/72.u1(k1) = (d1(k1)-e1(k1)/sqrt(var1(k1)ckk=0do k2 = 1,k1if(xp(k1).gt.xp(k2) thenkk=kk+1elseend ifend dod2(k

5、1)=d2(k1-1)+float(kk)U2(k1) = (d2(k1)-e1(k1)/sqrt(var1(k1)end dou3(1) = -u2(m)do m1 = 2,mu3(m1) = -u2(m-m1 + 1)end doreturnend突變分為如下主要的幾種：1。均值突變（最常見）；2。方差突變；3。線性回歸突變（也稱趨勢(shì)突變）；4。概率突變；5。空間型突變;6。譜突變；7。模型參數(shù)突變。等等。每種形式的突變都有其特定的提法，不了解其固有提法對(duì)真正深入的討論是沒有幫助的，而且，實(shí) 際分析時(shí)也是無的放矢。這些突變都有若干種檢測(cè)方法。1。 均值突變的檢測(cè)方法有：（1）低通濾波法；（

6、2）t-檢驗(yàn)法；（3）Cramer法；（4）滑動(dòng)t-檢驗(yàn)法;（5）Yamaoto 法；（6）M-K 法；（7）最小二乘法；（8）連續(xù)滑動(dòng) t-檢驗(yàn)法；（9）Pttitt 法；（10） Lepage法；（11）局部比較法。等等。在均值突變分析中，除了突變點(diǎn)、突變個(gè)數(shù)的估計(jì)外，一般還要 分析躍變度等。2。 方差突變的檢測(cè)方法有：（1）F一檢驗(yàn)法；（2）滑動(dòng)F檢驗(yàn)法；（3）連續(xù)滑動(dòng)F一檢驗(yàn)法。等等。3。線性回歸突變（也稱趨勢(shì)突變）的檢測(cè)方法有：（1）最小二乘法；（2）局部比較法；滑動(dòng)參數(shù)法等c This is a program for testing climate jumpi by use of

7、 c Mann-Kendall test.cparameter(iy=100)dimension x(iy),u1(iy),u2(iy)cRead Datacdo i = 1,50 x(i) = 0.x(i+50) = 1.enddocMann-Kendall test methodcall MKtest(iy,x,u1,u2)cdo i = 1,iywrite(31,10)i,x(i),u1(i),u2(i)end do10 format(1x,i4,4f8.2)stopendcsubroutine MKtest(m,x,u1,u2)dimension x(m),u1(m),u2(m)dimension xr(m),ur(m) !work arraydo i = 1,mxr(i)=x(m + 1-i)end docall MKrank(m,x,u1)call MKrank(m,xr,ur)do i = 1,mu2(i) = -ur(m + 1-i)enddoreturnendc-cMann-Kendall Rank Statisticsubroutine MKrank(m,x,u)dimension x(m),u(m)u(1)=0.0d=0.0do i = 2,mnum=0do j = 1,iif(x(i