求mann-kendall趋势检验的fortran程序

柠檬/db

要学的东东好多 啊，好急需，谢谢楼主

phenix2010

楼主试试这个，
C        THIS IS A PROGRAM FOR DETECTING ABRUPT CLIMATIC CHANGE
C        BY USING MANN-KENDALL TECHNIQUE
        PROGRAM MK
        DIMENSION Y(1000),YY(1000),U(1000),UF(1000),UB(1000),
     &        M(1000),MD(1000)
        WRITE(*,10)
  10        FORMAT(2X,'N=?,NYEAR=?')
        READ(*,*)N,NYEAR
C        ***************************************************
C        * N:     SAMPLE SIZE                              *
C        * NYEAR: FIRST YEAR OF THE TIME SERIES            *
C        * Y(N):  ORIGINAL TIME SERIES                     *
C        * UF(N): ORIGINAL SERIES OF U(LN)                 *
C        * UB(N): COUNTER SERIES OF U(LN)                  *
C        * A,B:   CRITICAL VALUE 1.96 AND -1.96            *
C        ***************************************************          
        OPEN(2,FILE=' ')
        READ(2,*)(Y(I),I=1,N)
        CALL SMK(Y,M,MD,UF,N)
        DO 20 I=1,N
  20        YY(I)=Y(N+1-I)
        CALL SMK(YY,M,MD,U,N)
        DO 30 I=1,N
  30        UB(I)=-U(N+1-I)
        OPEN(3,FILE=' ',STATUS='NEW')
        A=1.96
        B=-1.96
        DO 40 I=1,N
        WRITE(3,50)NYEAR+I-1,UF(I),UB(I),A,B
  50        FORMAT(1X,I4,4F8.2)
  40        CONTINUE
        CLOSE(3)
        STOP
        END
C***********************************************************
        SUBROUTINE SMK(Y,M,MD,U,N)
        DIMENSION Y(N),M(N),MD(N),U(N)
        M(1)=0
        DO 10 I=2,N
        M(I)=0
        MD(I)=0
        DO 20 J=1,I-1
        IF(Y(I).LT.Y(J))GOTO 20
        M(I)=M(I)+1
  20        CONTINUE
        MD(I)=MD(I-1)+M(I)
  10        CONTINUE
        U(1)=0.0
        DO 30 I=2,N
        E=I*(I-1)/4.00
        VAR=I*(I-1)*(2*I+5)/72.00
        U(I)=(MD(I)-E)/SQRT(VAR)
  30        CONTINUE
        RETURN
        END

mmpyb

挺难学的东西

lacie翼

jane1124 发表于 2011-8-9 17:52

                               
登录/注册后可看大图

不知道你能否看到图片，就是要求这个bita值，学得不好，不会编程。

lz的这篇文章是那篇？

西门杨

谢谢，长学问了

meredithxf

我非常想要这个matlab的程序哇

Lorenz混沌

这程序是突变分析的，不是趋势

zhiaicaizi

phenix2010 发表于 2011-8-18 10:30

                               
登录/注册后可看大图

楼主试试这个，
C        THIS IS A PROGRAM FOR DETECTING ABRUPT CLIMATIC CHANGE
C        BY USING MANN-KENDALL TE ...

比如说某地区40年（1900-1940）数据，想知道气温的突变，程序中N代表什么？是40年的样本数据，即N=40，那么NYEAR代表是1900年吗？尽管程序中说代表时间序列的第一年，我是初学者，没有明白什么意思，求指教？

三江魂

看不到图片

birdbuaa

也需要趋势MK程序