带通滤波程序提取周期的方法

cjjj

请问以下这个滤波程序具体的提取周期的方法（这个程序也是在本网站下载的，稍有改动，在此感谢程序提供者）：
C-------------PROGRAM MAIN---------------C
c  Array hgt(nx,ny,nt) shows the HGT data
c  'nx' is the number of grids in x-direction
c  'ny' is the number of grids in y-direction
c  'nt' is day
c  'mw1' is bottom days
c  'mw2' is top days
c----------------------------------------C
      parameter(nx=145,ny=73,nt=365,mw1=30,mw2=60)
      dimension hgt(nx,ny,nt),hgt1(nx,ny,nt)
      dimension cx(nt),xy(nt),bx(nt),tt(nt)
c-----------input data------------------

      open(11,file='f:\uwnd.sig995.2003.nc',form='unformatted')
      read(11) (((hgt(i,j,it),i=1,nx),j=1,ny),it=1,nt)
      close(11)
       
       
       
      do 70 i=1,nx
      do 70 j=1,ny
      do 50 it=1,nt
  50  tt(it)=hgt(i,j,it)
      call filter(nt,tt,bx,cx,xy,mw1,mw2)
      do 60 it=1,nt
  60  hgt1(i,j,it)=bx(it)
  70  continue

c----writing hgt1(nx,ny,nt) to file 'hgt.2000lb.grd'-----
      open(12,file='F:\uwnd\uwnd.sig995.2003.grd',form='binary')

        do 80 it=1,nt

  80  write(12) ((hgt1(i,j,it),i=1,nx),j=1,ny)
      close(12)
       
      end

c
      subroutine filter(n,x,bx,cx,xy,mw1,mw2)
      dimension x(n),bx(n),cx(n),xy(n)
      w1=2*3.1415926/float(mw1)
      w2=2*3.1415926/float(mw2)
      w0=sqrt(w1*w2)
      dt=1.0
      ww=2*abs(sin(w1*dt)/(1+cos(w1*dt))
     $ -sin(w2*dt)/(1+cos(w2*dt)))
      ww0=4*sin(w1*dt)*sin(w2*dt)/((1+cos(w1*dt))
     $ *(1+cos(w2*dt)))
      a0=2*ww/(4+2*ww+ww0)
      b1=2*(ww0-4)/(4+2*ww+ww0)
      b2=(4-2*ww+ww0)/(4+2*ww+ww0)
      do 89 i0=1,n-2
      i=i0+2
      bx(1)=0.0
      bx(2)=0.0
      cx(i)=a0*(x(i)-x(i-2))-b1*bx(i-1)-b2*bx(i-2)
  89  bx(i)=cx(i)
      do 91 i=1,n
  91  x(i)=bx(n+1-i)
      do 92 i0=1,n-2
      i=i0+2
      xy(1)=x(1)
      xy(2)=x(2)
      cx(i)=a0*(x(i)-x(i-2))-b1*xy(i-1)-b2*xy(i-2)
  92  xy(i)=cx(i)
      do 93 i=1,n
  93  bx(i)=xy(n+1-i)
      end

!这是fortran的一个带通滤波程序，是提取介于n1-n2之间波段。  

hellofriend

提取周期的话用小波分析最好吧

cjjj

hellofriend 发表于 2013-8-26 16:06

                               
登录/注册后可看大图

提取周期的话用小波分析最好吧

抱歉~好久不上了~可能是小波分析比较好吧~可惜学艺不精 不会额。。。

李中昊

提取周期常用的还是小波分析，这个程序是用来滤波的

cjjj

李中昊 发表于 2015-5-14 14:57
提取周期常用的还是小波分析，这个程序是用来滤波的

后面也是用来滤波了

是打发到高

额不是很清楚呢

清柳雨菡

谢谢楼主！

biglittle金

李中昊 发表于 2015-5-14 14:57
提取周期常用的还是小波分析，这个程序是用来滤波的

提取周期变化曲线也是用小波分析么

biglittle金

李中昊 发表于 2015-5-14 14:57
提取周期常用的还是小波分析，这个程序是用来滤波的

提取周期变化曲线也是用小波分析么

cgccq

谢谢，学习了！