小波分析程序问题

Panda°W

本人在做小波分析过程中遇到些麻烦，以下是程序，在运行时候显示数据读取错误，想请教下高手可否帮忙解决
! -------------------- wavelet program -----------------------
      Program WAVELET
      parameter (ns=1,na=41,nb=58)
      real subh(nb,ns),a(na),b(na),c,d
      real wf(na,nb),ff(nb),wlet(na,ns,nb)
! character*98 str,chr*10
      real A0
data A0/1.0/
      open(10,file='d:\gradstt\sta_ave\ave9mo_752sta.txt')
      do i=1,nb
        read(10,'(f8.4)')(subh(i,is),is=1,ns)
      end do
do is=1,ns
        pj=0.0
   do i=1,nb
          pj=pj+subh(i,is)
        enddo
   pj=pj/nb
   do i=1,nb
          ff(i)=subh(i,is)-pj
        enddo
        call wavetran(A0,ff,wf,na,nb)
        do i=1,nb
     do il=1,na
            wlet(il,is,i)=wf(il,i)
     end do
        enddo
      end do
      open(20,file='d:\gradstt\wave\ave9mo_752sta.dat'
     &       ,form='unformatted',access='direct',recl=na*ns)
      do i=1,nb
        write(20,rec=i)((wlet(il,is,i),il=1,na),is=1,ns)
      end do



数据读入没问题了，可是写出dat文件时候又出问题了

问题解决了，果然是数据文本格式问题，数据是无格式的，我把它按列对齐弄成一维的就好了。至于写出dat文件的问题把上面语句改成下面的就ok了，大概就是因为格式不对的原因吧，具体怎么回事也不大懂，只是尝试了很多改法，瞎蒙蒙对了{:soso_e143:}希望对有同样问题的朋友有所帮助，但愿没有帮倒忙哦{:soso_e121:}open(20,file='d:\gradstt\wave\ave9mo_752sta.dat'
     &       ,form='binary')
      do i=1,nb
        write(20)((wlet(il,is,i),il=1,na),is=1,ns)

river

读取错误的提示呢，它肯定大概定位了啊。这么长的程序你让别人从头看到尾啊，还有你用的什么资料也要简单说一下啊

Panda°W

river 发表于 2013-5-20 20:43

                               
登录/注册后可看大图

读取错误的提示呢，它肯定大概定位了啊。这么长的程序你让别人从头看到尾啊，还有你用的什么资料也要简单说 ...

file:///C:/Users/wangh/AppData/Roaming/Tencent/Users/823707966/QQ/WinTemp/RichOle/7R]A]Y9QH0LU90)C%7D5_P_(K.jpg

数据用的就是58年总降水量的txt文件，里面是58个一维数据

river

Panda°W 发表于 2013-5-20 21:47

                               
登录/注册后可看大图

数据用的就是58年总降水量的txt文件，里面是58个一维数据

错误提示呢？既然你自己都说是一维数据了，怎么还读成两维的呢  read(10,'(f8.4)')(subh(i,is),is=1,ns)

Panda°W

river 发表于 2013-5-20 21:53

                               
登录/注册后可看大图

错误提示呢？既然你自己都说是一维数据了，怎么还读成两维的呢  read(10,'(f8.4)')(subh(i,is),is=1,ns)

可是ns定义了 ns=1.这个程序也是我在站内搜到的，不是自己编的

river

Panda°W 发表于 2013-5-20 22:11

                               
登录/注册后可看大图

可是ns定义了 ns=1.这个程序也是我在站内搜到的，不是自己编的

那也要根据自己的资料修改啊，数据只有一维，你非要读成二维的，可不就提示读资料的时候断掉了么

蔚蓝的天

我也来顶一个

llq

楼主请问你用的是什么小波？小弟正在做这个，准备用python写这个，遇到问题了