【已解决】小弟mtt，滑动t检验求助

jmcao

C        THIS IS A PROGRAM FOR DETECTING ABRUPT CLIMATIC CHANGE
C        BY USING MOVING t-TEST TECHNIQUE
        PROGRAM MTT
        DIMENSION Y(1000),YN1(1000),YN2(1000),YC1(1000),YC2(1000),
     &        S(1000),T(1000),NY1(1000),NY2(1000),NNY1(1000),NNY2(1000)
        WRITE(*,10)
  10        FORMAT(5X,'N=70,IH=5,NYEAR=1930')
        READ(*,*)N,IH,NYEAR
C        **************************************************
C        * N:    SAMPLE SIZE                              *
C        * IH:   LENGTH OF SUB-SERIES                         *
C        * NYEAR: FIRST YEAR OF THE TIME SERIES           *
C        * Y(N):  ORIGINAL TIME SERIES                    *
C        * ************************************************
        OPEN(2,FILE='d:\2.DAT')
        READ(2,*)(Y(I),I=1,N)
        N1=N-IH+1
        N2=N-2*IH+1
        C1=0.0
        C2=0.0
        DO 20 I=1,IH
  20        C1=C1+Y(I)
        DO 30 I=IH+1,2*IH
  30        C2=C2+Y(I)
        DO 40 I=1,N-IH
        D1=C1-Y(I)+Y(I+IH)
        YN1(I)=C1/IH
  40        C1=D1
        DO 50 I=IH+1,N-IH
        D2=C2-Y(I)+Y(I+IH)
        YN2(I-IH)=C2/IH
  50        C2=D2
        YN1(N1)=C1/IH
        YN2(N2)=C2/IH
        DO 60 I=1,N1
        YC1(I)=0.0
        DO 70 J=I,IH+I-1
  70        YC1(I)=YC1(I)+(Y(J)-YN1(I))*(Y(J)-YN1(I))
  60        CONTINUE
        DO 80 I=1,N2
        YC2(I)=0.0
        DO 90 J=IH+I,2*IH+I-1
  90        YC2(I)=YC2(I)+(Y(J)-YN2(I))*(Y(J)-YN2(I))
  80        CONTINUE
        DO 100 I=1,N2
  100        S(I)=SQRT((YC1(I)+YC2(I))/(IH+IH-2))
        DO 110 I=1,N2
  110        T(I)=(YN1(I)-YN2(I))/(S(I)*SQRT(2.0/IH))
        DO 120 I=1,N2
        NY1(I)=NYEAR+I-1
        NNY1(I)=NY1(I)+IH-1
        NY2(I)=NYEAR+I+IH-1
  120        NNY2(I)=NY2(I)+IH-1
        IF(IH.EQ.5)THEN
        A=4.0
        B=-4.0
        ELSE IF(IH.EQ.10)THEN
        A=3.2
        B=-3.2
        ELSE
        A=3.0
        B=-3.0
        END IF
        OPEN(3,FILE='d:\3.DAT',STATUS='NEW')
        WRITE(3,130)
  130        FORMAT(30X,'ABRUPT CLIMATIC CHANGE ANALYSIS'/)
        WRITE(3,140)
  140        FORMAT(25X,'TM',4X,'t-TEST(0.01)'/)
        DO 150 I=1,N2
        WRITE(3,160)NY1(I),NNY1(I),NY2(I),NNY2(I),T(I),A,B
  160        FORMAT(1X,I4,'-',I4,'--',I4,'-',I4,1X,3F8.2)
  150        CONTINUE
        OPEN(4,FILE='d:\4.DAT',STATUS='NEW')
        DO 170 I=1,N2
        WRITE(4,180)NNY1(I),T(I),A,B
  180        FORMAT(1X,I4,3F8.2)
  170        CONTINUE
        STOP
        END
我想求的是，1960-2013，气温序列，N=54，IH=5，
数据是，2.txt，希望得到的是，mtt.txt
大一新生，网上看到t检验的f程序，求各位大侠，哪些地方需要改动的，帮忙指出来，
OPEN(3,FILE='d:\3.DAT',STATUS='NEW')这个意思不明白，不是就一个读入的数据文件2.txt吗
谢谢

zhouling

大一新生就搞这个？

言深深

首先你要装一个fortran，并且掌握fortran程序的操作；然后按照：
1、你已经有的文件数据放到2.dat文件里面，和fortran文件装在同一个英文路径的文件夹下面；
2、运行程序之后，依次输入，”54,5,1960“，按回车运行程序，结束之后输出来的4.dat就是结果。
结果分析什么的，自己查书吧！

最后，我不能容忍的是，这个是MTT，不是mmt

Kongxiangye

言深深 发表于 2014-7-18 10:55
首先你要装一个fortran，并且掌握fortran程序的操作；然后按照：
1、你已经有的文件数据放到2.dat文件里面 ...

这个版主很可爱

jmcao

额，程序里面哪部分是需要改动的？这个编程是拷贝的，所以有的看不明白，一定.dat吗，可不可以.txt？

lqouc

jmcao 发表于 2014-7-18 14:17
额，程序里面哪部分是需要改动的？这个编程是拷贝的，所以有的看不明白，一定.dat吗，可不可以.txt？

改动的部分就是注释里面提到的参数，英文总不能看不懂吧。
读入数据的扩展名随便什么都行，只要数据排列方式对就好。
建议楼主从基础学起，不要上来就搞这种程序，否则大家怎么解释你也看不明白。

言深深

jmcao 发表于 2014-7-18 14:17
额，程序里面哪部分是需要改动的？这个编程是拷贝的，所以有的看不明白，一定.dat吗，可不可以.txt？

再补充一句：论坛的回复一定要点击我那边的回复我才能看到，你这个回复你自己的帖子我是看不到的·······

jmcao

分别是，IH=10，IH=5的t曲线，mk的交点有3个，1993，1994,1995，到底哪年突变啊

jmcao

图1过线的是
1988-1997--1998-2007        -3.59
1992-2001--2002-2011        -3.42
所以，到底哪年突变啊？

jmcao

言深深 发表于 2014-7-18 15:44
再补充一句：论坛的回复一定要点击我那边的回复我才能看到，你这个回复你自己的帖子我是看不到的···· ...

没事版主~好了谢谢~再帮忙看看帖子呗