24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

calos818

木虫 (著名写手)

新人菜鸟

应助: 11 (小学生)
金币: 3135.7
散金: 77
红花: 3
帖子: 1200
在线: 730.5小时
虫号: 698507
注册: 2009-02-08
专业: 凝聚态物性 II ：电子结构

[交流] 【求助】求助：快速傅里叶变换为何码位倒置？

在看Numerical Recipes 上的FFT源程序，采用2进制，设N=8,000-111，搞不懂为何要bit reverse，如100变成001，010变为010等，程序在下面，请各位指导下，万分感谢

回复此楼

» 猜你喜欢

1楼 2009-08-09 09:39:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

calos818

木虫 (著名写手)

新人菜鸟

应助: 11 (小学生)
金币: 3135.7
散金: 77
红花: 3
帖子: 1200
在线: 730.5小时
虫号: 698507
注册: 2009-02-08
专业: 凝聚态物性 II ：电子结构

还是看不懂这个FFT

赞一下

回复此楼

4楼2009-08-09 13:20:07

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 9 个回答

calos818

木虫 (著名写手)

新人菜鸟

应助: 11 (小学生)
金币: 3135.7
散金: 77
红花: 3
帖子: 1200
在线: 730.5小时
虫号: 698507
注册: 2009-02-08
专业: 凝聚态物性 II ：电子结构

源代码

SUBROUTINE four1(data,nn,isign)
   INTEGER isign,nn
   REAL data(2*nn)
   INTEGER i,istep,j,m,mmax,n
   REAL tempi,tempr
   DOUBLE PRECISION theta,wi,wpi,wpr,wr,wtemp
   n=2*nn
   j=1
   do 11 i=1,n,2
      if(j.gt.i)then
      tempr=data(j)
      tempi=data(j+1)
      data(j)=data(i)
      data(j+1)=data(i+1)
      data(i)=tempr
      data(i+1)=tempi
      endif
      m=n/2
1    if ((m.ge.2).and.(j.gt.m)) then
      j=j-m
      m=m/2
      goto 1
      endif
      j=j+m
11 continue
   mmax=2
2    if (n.gt.mmax) then
      istep=2*mmax
      theta=6.28318530717959d0/(isign*mmax)
      wpr=-2.d0*sin(0.5d0*theta)**2
      wpi=sin(theta)
      wr=1.d0
      wi=0.d0
      do 13 m=1,mmax,2
      do 12 i=m,n,istep
         j=i+mmax
         tempr=sngl(wr)*data(j)-sngl(wi)*data(j+1)
         tempi=sngl(wr)*data(j+1)+sngl(wi)*data(j)
         data(j)=data(i)-tempr
         data(j+1)=data(i+1)-tempi
         data(i)=data(i)+tempr
         data(i+1)=data(i+1)+tempi
12       continue
      wtemp=wr
      wr=wr*wpr-wi*wpi+wr
      wi=wi*wpr+wtemp*wpi+wi
13    continue
      mmax=istep
   goto 2
   endif
   return
   END
C  (C) Copr. 1986-92 Numerical Recipes Software ,4-#.

赞一下

回复此楼

2楼2009-08-09 09:40:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖