怎样用小波分析的程序计算周期分析？

殷2333

ccccc小波的模方差贡献ccccc
        open(120,file='f:\xb\mo.grd',
     &        form='binary')
!        do i=1,mn
!     write(*,10) (w2(i,j),J=1,mm)
!        enddo

        do j=1,MM
                da(j)=0.0
                do i=1,MN
                        da(j)=da(j)+w1(i,j)
                enddo
        enddo
        do i=1,MM
        write(120) da(i)
        enddo
CCCCCCCCCCCCCCCCCCC写小波分析的不同时间周期的功率谱
!        open(122,file='g:\xb\xb\pu.txt',form='formatted')
       
        do j=1,MM
!        if((j.eq.5).or.(j.eq.9).or.(j.eq.15))then
       
!        write(122,'(47f6.2)') (iw(i,j),i=1,MN)
       
!     endif
        enddo


ccccc小波的模ccccc
      open(119,file='f:\xb\pattern1.grd',form='binary')
        write(119) ((w2(i,j),I=1,mn),J=1,mm)
cccccc小波的实部ccccc  
        open(188,file='f:\xb\real-j1.grd',form='binary')
        write(188) ((rw1(i,j),J=1,mm),I=1,mn)
cccccc小波的位相ccccc
      open(189,file='f:\xb\phase1.grd',form='binary')
      write(189) ((wx1(i,j),I=1,mn),j=1,mm)
cccccc小波的虚部ccccc     
        open(121,file='f:\xb\imaginary1.grd',form='binary')
        write(121) ((iw(i,j),I=1,mn),J=1,mm)

        close(119)
        close(188)
        close(189)
        close(121)
!---------------------------------------------------
      END
!---------------------------------------------
        SUBROUTINE RMORLET(T,A,B,PA,RMO)
       DOUBLE PRECISION Q,TT,PA,T,A,B,RMO
       TT=(T-B)/A
   !    Q=DBLE(1.0)
       Q=-((ABS(TT))**2)
       Q=Q/2           
       RMO=DCOS(2*PA*TT)*DEXP(Q)
        RETURN
       END
!---------------------------------------------
         SUBROUTINE  RIMORLET(T,A,B,PA,RIMO)
       DOUBLE PRECISION Q,TT,T,A,B,pa,RIMO
       TT=(T-B)/A
   !    Q=DBLE(1.0)
       Q=-((ABS(TT))**2)
       Q=Q/2.0
       RIMO=DSIN(2*PA*TT)*DEXP(Q)
        RETURN
       END
!---------------------------------------------
       SUBROUTINE NORMAL(N,X,XX)
       DOUBLE PRECISION X(N),XX(N),PX,S
       PX=0.d0
       DO I=1,N
        PX=PX+X(I)
       ENDDO
         PX=PX/FLOAT(N)
       S=0.d0
       DO I=1,N
        S=S+(X(I)-PX)**2
         ENDDO
       S=DSQRT(S/FLOAT(N))
       DO I=1,N
        XX(I)=(X(I)-PX)/S
       ENDDO
       END
!--------------------------------------------
      SUBROUTINE FILTER(AA,ITX,DT,PERC,PERL)
      DOUBLE PRECISION AA(ITX),WK1(ITX),WK2(ITX)
        DOUBLE PRECISION PI2,WO,WOS,W1,W2,DW,DWS,B1,B2,B3
        PI2=2.0*3.1415926d0
        WO=PI2/PERC
        WOS=WO*WO
        W1=PI2/PERL
        W2=WOS/W1
        DW=2.0*ABS((SIN(W1*DT))/(1.0+COS(W1*DT))-(SIN(W2*DT))/(1.0+
     &   COS(W2*DT)))
      DWS=(4.0*SIN(W1*DT)*SIN(W2*DT))/((1.0+COS(W1*DT))*(1.0+
     &    COS(W2*DT)))
      B3=(2.0*DW)/(4.0+2.0*DW+DWS)
        B2=(4.0-2.0*DW+DWS)/(4.0+2.0*DW+DWS)
        B1=(2.0*(DWS-4.0))/(4.0+2.0*DW+DWS)
      ITXM=ITX-1
        DO IT=1,ITX
       WK1(IT)=0.0
        ENDDO
        DO IT=3,ITX
       WK1(IT)=B3*(AA(IT)-AA(IT-2))-B1*WK1(IT-1)-B2*WK1(IT-2)
        ENDDO
        WK2(ITX-1)=WK1(ITX-1)
        WK2(ITX)=WK1(ITX)
        WK1(1)=AA(1)
        WK1(2)=AA(2)
        DO IT=2,ITXM
       II=ITX-IT
       WK2(II)=B3*(WK1(II)-WK1(II+2))-B1*WK2(II+1)-B2*WK2(II+2)
        ENDDO
        DO IT=1,ITX
       AA(IT)=WK2(IT)
        ENDDO
        END
这打开的都是什么鬼？

殷2333

模、位相、实部、虚部都是什么？

殷2333

有没有大神在？

lxf1896

MM是如何计算出来的，是窗口长度吗？