| 查看: 3479 | 回复: 5 | ||
panjunxing木虫 (著名写手)
小学生
|
[求助]
三维傅里叶变换 已有1人参与
|
|
将一组16*32*64的三维数组用fortran程序进行傅立叶变换，主程序该如何编写？ 下面是傅立叶变换源程序，可惜我不会用  ，另外附上自己的程序，请各位大侠多多指教！SUBROUTINE rlft3(data,speq,nn1,nn2,nn3,isign) INTEGER isign,nn1,nn2,nn3 COMPLEX data(nn1/2,nn2,nn3),speq(nn2,nn3) INTEGER i1,i2,i3,j1,j2,j3,nn(3) DOUBLE PRECISION theta,wi,wpi,wpr,wr,wtemp COMPLEX c1,c2,h1,h2,w c1=cmplx(0.5,0.0) c2=cmplx(0.0,-0.5*isign) theta=6.28318530717959d0/dble(isign*nn1) wpr=-2.0d0*sin(0.5d0*theta)**2 wpi=sin(theta) nn(1)=nn1/2 nn(2)=nn2 nn(3)=nn3 if(isign.eq.1)THEN call fourn(data,nn,3,isign) do i3=1,nn3 do i2=1,nn2 speq(i2,i3)=data(1,i2,i3) end do end do endif do i3=1,nn3 j3=1 if (i3.ne.1) j3=nn3-i3+2 wr=1.0d0 wi=0.0d0 do i1=1,nn1/4+1 j1=nn1/2-i1+2 do i2=1,nn2 j2=1 if (i2.ne.1) j2=nn2-i2+2 if(i1.eq.1)then h1=c1*(data(1,i2,i3)+conjg(speq(j2,j3))) h2=c2*(data(1,i2,i3)-conjg(speq(j2,j3))) data(1,i2,i3)=h1+h2 speq(j2,j3)=conjg(h1-h2) else h1=c1*(data(i1,i2,i3)+conjg(data(j1,j2,j3))) h2=c2*(data(i1,i2,i3)-conjg(data(j1,j2,j3))) data(i1,i2,i3)=h1+w*h2 data(j1,j2,j3)=conjg(h1-w*h2) endif end do wtemp=wr wr=wr*wpr-wi*wpi+wr wi=wi*wpr+wtemp*wpi+wi w=cmplx(sngl(wr),sngl(wi)) end do end do if(isign.eq.-1)then call fourn(data,nn,3,isign) endif return END |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 1.f90
2014-03-12 17:53:13, 9.72 K
» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 聚焦生物仿制药系列
已经有6人回复
- 蛋白质的共价修饰的测定的试验方法----凌波丽 个人总结
已经有19人回复
- matlab快速切换图
已经有6人回复
- 差值傅里叶图怎么看？
已经有3人回复
- Matlab二维海浪谱的绘制方法
已经有7人回复
- 请教一个傅里叶变换坐标顺序问题
已经有5人回复
- 已知一随机信号的功率谱密度函数，如何利用Matlab产生其单次时域样本
已经有6人回复
- 傅立叶变换的物理意义
已经有9人回复
- 【求助】氢原子波函数什么时候用指数函数，什么时候用三角函数？
已经有7人回复
- 【讨论】多普勒信号频率的测量方法
已经有6人回复
- 【交流】有关布里渊区、K空间的问题
已经有11人回复
baobiao007
木虫 (职业作家)
中国特色
- 应助: 201 (大学生)
- 金币: 6484.2
- 散金: 557
- 红花: 40
- 帖子: 3050
- 在线: 1009.9小时
- 虫号: 505962
- 注册: 2008-02-18
- 专业: 应用地球物理学
2楼2014-03-12 17:53:47
panjunxing
木虫 (著名写手)
小学生
- 应助: 0 (幼儿园)
- 金币: 2994.8
- 散金: 948
- 红花: 2
- 帖子: 1103
- 在线: 69.8小时
- 虫号: 579223
- 注册: 2008-07-04
- 性别: MM
- 专业: 高分子组装与超分子结构
您能帮我一下吗？
3楼2014-03-13 09:10:20
cooooldog
铁杆木虫 (著名写手)
ส็็็
- 应助: 237 (大学生)
- 金币: 6101.9
- 散金: 1114
- 红花: 39
- 帖子: 1380
- 在线: 553.8小时
- 虫号: 506699
- 注册: 2008-02-18
- 专业: 模式识别
【答案】应助回帖
感谢参与，应助指数 +1
|
楼主不如用这里的代码: http://fftwpp.sourceforge.net/ 可以从python, C++, Fortran调用,而且都有例子, 可以从例子里面学习用法 |
4楼2014-03-13 09:37:25
cooooldog
铁杆木虫 (著名写手)
ส็็็
- 应助: 237 (大学生)
- 金币: 6101.9
- 散金: 1114
- 红花: 39
- 帖子: 1380
- 在线: 553.8小时
- 虫号: 506699
- 注册: 2008-02-18
- 专业: 模式识别
5楼2014-03-13 09:38:36
lixinwen000
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 4182.6
- 散金: 274
- 红花: 3
- 帖子: 717
- 在线: 287.8小时
- 虫号: 1540526
- 注册: 2011-12-17
- 性别: GG
- 专业: 海岸工程
6楼2014-03-25 15:42:19
