| 查看: 3258 | 回复: 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
» 猜你喜欢
335求调剂
已经有4人回复
-
各位老师好，我的一志愿为北京科技大学085601材料专硕
已经有12人回复
-
359求调剂
已经有7人回复
-
070300化学专业279调剂
已经有10人回复
-
环境工程297分求调剂一志愿杭高院
已经有9人回复
-
一志愿武理材料工程302调剂环化或化工
已经有9人回复
-
309分085801求调剂
已经有5人回复
-
一志愿 南京航空航天大学 ，080500材料科学与工程学硕
已经有5人回复
-
材料调剂
已经有4人回复
-
材料调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 聚焦生物仿制药系列
已经有6人回复
- 蛋白质的共价修饰的测定的试验方法----凌波丽 个人总结
已经有19人回复
- matlab快速切换图
已经有6人回复
- 差值傅里叶图怎么看？
已经有3人回复
- Matlab二维海浪谱的绘制方法
已经有7人回复
- 请教一个傅里叶变换坐标顺序问题
已经有5人回复
- 已知一随机信号的功率谱密度函数，如何利用Matlab产生其单次时域样本
已经有6人回复
- 傅立叶变换的物理意义
已经有9人回复
- 【求助】氢原子波函数什么时候用指数函数，什么时候用三角函数？
已经有7人回复
- 【讨论】多普勒信号频率的测量方法
已经有6人回复
- 【交流】有关布里渊区、K空间的问题
已经有11人回复
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
baobiao007
主管区长
中国特色
- 应助: 201 (大学生)
- 金币: 6482.7
- 散金: 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
- 专业: 模式识别
5楼2014-03-13 09:38:36
