求高手帮助滑动t检验

学长吉他

我想做60年数据的滑动t检验，步长选择 5和10，下面是我的程序请问高手需要改一下哪里？另外我的数据存放是60个数据一列这样对吗？
PROGRAM N_MTT
PARAMETER (N=60 ,K=5 ,NYEAR=1)
REAL *8 X(N),X1(N),X2(N),T(K:N-K+1)
REAL *8 V1,V2,EQ1,EQ2
OPEN (1,FILE="C:\Users\James\Desktop\001.TXT")
DO I=1,N
READ (1,*) X(I)
ENDDO
CLOSE(1)

IF (K.EQ.5) THEN
A=2.31
ELSE IF (K.EQ.10) THEN
A=2.10
ELSE IF (K.EQ.12) THEN
A=2.05
ENDIF
DO L=K,N-K+1
DO I=1,K
X1(I)=X(L-K+I)
X2(I)=X(L+I-1)
ENDDO
CALL SVAR(X1,V1,K,N)
CALL SVAR(X2,V2,K,N)
CALL SAVAG(X1,EQ1,K,N)
CALL SAVAG(X2,EQ2,K,N)
T(L)=(EQ1-EQ2)/SQRT((2.0/K)*((K-1)*(V1+V2)/(K+K-2)))
PRINT *,NYEAR+L-1
ENDDO
OPEN (2,FILE="C:\Users\James\Desktop\001.TXT")
DO L=K,N-K+1
WRITE (2,'(I4,3F18.2)') NYEAR+L-1,T(L),A,-A
ENDDO
CLOSE(2)
END
SUBROUTINE SVAR(X,V,K,N)
REAL *8 X(K),V
SUM=0
DO I=1,K
SUM=SUM+X(I)
ENDDO
AVAG=SUM/K
VAR=0
DO I=1,K
VAR=VAR+(X(I)-AVAG)**2
ENDDO
V=VAR/K
END
SUBROUTINE SAVAG(X,EQ,K,N)
REAL *8 X(K),EQ
SUMX=0
DO I=1,K
SUMX=SUMX+X(I)
ENDDO
EQ=SUMX/K
END

lqouc

学长吉他 发表于 2013-3-27 10:11

                               
登录/注册后可看大图

谢谢啊，呵呵，我的数据出来了，现在剩下画图了，你做完了吗？

最后考虑再三觉得太主观了，于是放弃这个方法了，依然用mk检验。其实滑动t检验检验也不是万能，排除步长的主观因素不说，感觉从原理上面只能判断均值的突变，不能很好的判断振幅和周期的突变（虽然可以通过步长调整解决一部分）。所以的所以......现在就这样了，以后再研究研究别的方法好了。

学长吉他

有人帮帮我吗

kongfeng0824

帮顶！！！！！！！！！

学长吉他

kongfeng0824 发表于 2013-3-25 22:22

                               
登录/注册后可看大图

帮顶！！！！！！！！！

谢谢啊，可是还是没人解答啊

lqouc

学长吉他 发表于 2013-3-26 10:32

                               
登录/注册后可看大图

谢谢啊，可是还是没人解答啊

你这个程序没有试试么？难道出错了么？六十个数一列是没问题的，要说问题的话最好把开始时间改一下，不要写1，改好了得到的结果就可以直接画图了，结果是对应了时间的。
我也是今天要用这个临时看了一眼，不一定对，多交流吧。

学长吉他

lqouc 发表于 2013-3-26 14:17

                               
登录/注册后可看大图

你这个程序没有试试么？难道出错了么？六十个数一列是没问题的，要说问题的话最好把开始时间改一下，不要 ...

我试了啊，出来的窗口是这样的不懂什么意思，你帮我再看看

lqouc

学长吉他 发表于 2013-3-26 16:18

                               
登录/注册后可看大图

我试了啊，出来的窗口是这样的不懂什么意思，你帮我再看看

那个意思是说你的文件读到末尾了，也就是说你的文件没有那么多，你确定是一行一个数，一共60行么？实在不行你就加一个iostat的判断，但是这个治标不治本，估计后面也还会出错。还是不行的话可以试试魏凤英的版本，或者直接换mk好了，因为这个的步长你还要根据分析的震荡周期自己调，主观性太强了。

学长吉他

lqouc 发表于 2013-3-26 16:30

                               
登录/注册后可看大图

那个意思是说你的文件读到末尾了，也就是说你的文件没有那么多，你确定是一行一个数，一共60行么？实在不 ...

谢谢啊，呵呵，我的数据出来了，现在剩下画图了，你做完了吗？

学长吉他

学长吉他 发表于 2013-3-27 10:11

                               
登录/注册后可看大图

谢谢啊，呵呵，我的数据出来了，现在剩下画图了，你做完了吗？

你QQ多少？我加你