mann-kendall突变检验.doc

《mann-kendall突变检验.doc》由会员分享，可在线阅读，更多相关《mann-kendall突变检验.doc（9页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、Four short words sum up what has lifted most successful individuals above the crowd: a little bit more.-author-datemann-kendall突变检验SPSS中的MK检验SPSS中的MK检验program mainimplicit nonec This is a program for testing climate jumpi by use ofc Mann-Kendall test.c-integer,parameter : iy=100integer ireal x(iy),u

2、1(iy),u2(iy)c-Read Datac-open(31,file=d:mk.txt,form=formatted)do i=1,50x(i)=0.5x(i+50)=-0.5enddoc-Mann-Kendall test methodcall MKtest(iy,x,u1,u2)c-do i=1,iywrite(31,10)i,x(i),u1(i),u2(i)end do10 format(1x,i4,4f8.2)stopend!c-subroutine MKtest(m,x,u1,u2)dimension x(m),u1(m),u2(m)dimension xr(m),ur(m)

3、!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)enddoreturnend!c-!c-Mann-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).gt.x(j)thennum=num+1endifenddod=d+float(num)e=i*(i-1.)/4.var=i*(i-

4、1.)*(2.*i+5.)/72.u(i)=(d-e)/sqrt(var)enddoreturnendMK(MANN-KENDALL)是气象学/气候学中经常用来进行突变检验的一种方法，具体方法如下：设原始时间序列为y1，y2，yn，mi表示第i个样本yi大于yj（1ji）的累积数，定义统计量：.在原序列随机独立等假设下，dk的均值和方差分别为：将上面公式的dk标准化，得：UFk组成一条UF曲线，通过信度检验可得出其是否有明显的变化趋势。5把此方法引用到反序列中，计算得到另一条曲线UB，则两条曲线在置信区间内的交点确定为突变点。给定显著性水平=0.05，则统计量UF和UB的临界值为1.96。UF0，表示序列呈上升趋势；反之，表明呈下降趋势，大于或小于1.96，表示上升或下降趋势明显。.-