小波分析的程序有点错误算不出结果，帮忙看下

wysyhd

程序
!C THIS IS A PROGRAM OF WAVELET ANALYSIS BY USING
!       GAUSS REAL OF MCXICAN HAT
        PROGRAM WA
         PARAMETER (N=61,M=20,NMAX=N*3,NS=0,KS=1)
         DIMENSION F(M,N),X(N),XM(NMAX),IT(M)
         REAL D,Y,Z
! *****************************************************
!        * N:  SAMPLE SIZE                            *
!        * M:  NUMBER OF TIME SCALE                   *
! * NS: STARTING SCALE                                *
! * KS: INTERVAL OF SCALE                             *
! * X:  PRIMARY TIME SERIES                           *
! * F:  WAVELET TRANSFORMATION COEFFICIENTS           *
! * IT: EXPANDED AND CONTRACTED  SCALES               *
! *****************************************************
        WAVE(Y)=(1.0-16.0*Y*Y)*EXP(-8.0*Y*Y)
        OPEN(6,FILE='c:\lw\pr_qypj.grd',form='binary')
         READ(6,*)(X(I),I=1,N)
         DO  I=1,N
         XM(I)=X(N+1-I)
         XM(N+I)=X(I)
         XM(2*N+I)=X(N+1-I)
   enddo
         DO  I=1,M
         IT(I)=I*KS+NS
   enddo
         DO  J=1,N
         F(I,J)=0.0
   enddo
        DO  I=1,M
        D=FLOAT(IT(I))
        II=IT(I)
        DO  J=N+1,2*N
        DO  K=J-II,J+II
        Z=(FLOAT(K)-FLOAT(J))/D
        F(I,J-N)=F(I,J-N)+WAVE(Z)*XM(K)/SQRT(D)
   enddo
   enddo
   enddo
        OPEN(3,FILE='c:\lw\Wavelet.txt',STATUS='NEW')
        WRITE(3,70)
   70 FORMAT(20X,'*****WAVELET ANALYSIS*****'/)
        WRITE(3,80)
   80 FORMAT(2X,'NO.',3X,'SCALE',4X,'COE.')
         DO  I=1,M
         DO  J=1,N
         WRITE(3,110)J,IT(I),F(I,J)
  110 FORMAT(2X,I3,3X,I3,3X,F8.2)
enddo
enddo
         STOP
         END
!!-------------------------------END---------------------------------

最后运行会出下面的问题，求大神帮忙，谢谢，，，


forrtl: severe (257): formatted I/O to unit open for unformatted transfers, unit
6, file E:\lunwen\qypj.grd

river

错误提示说的很明白了，读取文件的时候出错，你那个是二进制文件，是无格式有数据类型的文件，而你的READ(6,*)(X(I),I=1,N)里面的*是用来读有格式无数据类型的文件的，所以需要改成READ(6)(X(I),I=1,N)

曾记否

READ(6,*)(X(I),I=1,N)读取格式错误，去掉星号

wysyhd

谢谢，两位大神，

wysyhd

river 发表于 2013-5-17 22:07

                               
登录/注册后可看大图

错误提示说的很明白了，读取文件的时候出错，你那个是二进制文件，是无格式有数据类型的文件，而你的READ(6 ...

谢谢啊，还要说声对不起了，扣了你的贡献

river

wysyhd 发表于 2013-5-18 00:14

                               
登录/注册后可看大图

谢谢啊，还要说声对不起了，扣了你的贡献

没关系，那个是论坛随机的，还挺好玩儿的。

冬日的阳光

        DO  J=1,N
         F(I,J)=0.0
   enddo
运行程序的时候，有没有显示程序运行出错？越界

空格

M=20是怎么确定的呢

zhouxiaosky

多谢分享~~~