| 查看: 688 | 回复: 8 | |||
| 当前主题已经存档。 | |||
calos818木虫 (著名写手)
新人菜鸟
|
[交流]
【求助】求助： 快速傅里叶变换 为何码位倒置？
|
||
| 在看Numerical Recipes 上的FFT源程序，采用2进制，设N=8,000-111，搞不懂为何要bit reverse，如100变成001，010变为010等，程序在下面，请各位指导下，万分感谢 |
回复此楼» 猜你喜欢
招博士
已经有6人回复
-
限项规定
已经有8人回复
-
国家基金申请书模板内插入图片不可调整大小?
已经有5人回复
-
交叉科学部支持青年基金，对三无青椒是个机会吗？
已经有3人回复
-
国家级人才课题组招收2026年入学博士
已经有5人回复
-
Fe3O4@SiO2合成
已经有6人回复
-
青年基金C终止
已经有4人回复
-
青椒八年已不青，大家都被折磨成啥样了？
已经有7人回复
-
为什么nbs上溴 没有产物点出现呢
已经有10人回复
-
救命帖
已经有11人回复
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
calos818
木虫 (著名写手)
新人菜鸟
- 应助: 11 (小学生)
- 金币: 3135.7
- 散金: 77
- 红花: 3
- 帖子: 1200
- 在线: 730.5小时
- 虫号: 698507
- 注册: 2009-02-08
- 专业: 凝聚态物性 II ：电子结构
3楼2009-08-09 09:49:32
calos818
木虫 (著名写手)
新人菜鸟
- 应助: 11 (小学生)
- 金币: 3135.7
- 散金: 77
- 红花: 3
- 帖子: 1200
- 在线: 730.5小时
- 虫号: 698507
- 注册: 2009-02-08
- 专业: 凝聚态物性 II ：电子结构
4楼2009-08-09 13:20:07
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
5楼2009-08-09 14:30:49
calos818
木虫 (著名写手)
新人菜鸟
- 应助: 11 (小学生)
- 金币: 3135.7
- 散金: 77
- 红花: 3
- 帖子: 1200
- 在线: 730.5小时
- 虫号: 698507
- 注册: 2009-02-08
- 专业: 凝聚态物性 II ：电子结构
6楼2009-08-09 14:53:07
calos818
木虫 (著名写手)
新人菜鸟
- 应助: 11 (小学生)
- 金币: 3135.7
- 散金: 77
- 红花: 3
- 帖子: 1200
- 在线: 730.5小时
- 虫号: 698507
- 注册: 2009-02-08
- 专业: 凝聚态物性 II ：电子结构
7楼2009-10-08 11:35:34
senlia
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 贵宾: 3.521
- 金币: 2178.1
- 红花: 1
- 帖子: 223
- 在线: 6.7小时
- 虫号: 849005
- 注册: 2009-09-16
- 性别: GG
- 专业: 计算数学与科学工程计算
8楼2009-10-10 09:34:59
jjdg
版主 (知名作家)
- 程序强帖: 9
- 应助: 223 (大学生)
- 贵宾: 3.926
- 金币: 51003.4
- 红花: 238
- 帖子: 7227
- 在线: 474.2小时
- 虫号: 27000
- 注册: 2003-10-21
- 专业: 药物代谢与药物动力学
- 管辖: 程序语言
9楼2009-10-11 12:08:28