| 查看: 3192 | 回复: 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
» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有10人回复
-
2025年遐想
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
-
自然科学基金委宣布启动申请书“瘦身提质”行动
已经有4人回复
-
求个博导看看
已经有18人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 聚焦生物仿制药系列
已经有6人回复
- 蛋白质的共价修饰的测定的试验方法----凌波丽 个人总结
已经有19人回复
- matlab快速切换图
已经有6人回复
- 差值傅里叶图怎么看？
已经有3人回复
- Matlab二维海浪谱的绘制方法
已经有7人回复
- 请教一个傅里叶变换坐标顺序问题
已经有5人回复
- 已知一随机信号的功率谱密度函数，如何利用Matlab产生其单次时域样本
已经有6人回复
- 傅立叶变换的物理意义
已经有9人回复
- 【求助】氢原子波函数什么时候用指数函数，什么时候用三角函数？
已经有7人回复
- 【讨论】多普勒信号频率的测量方法
已经有6人回复
- 【交流】有关布里渊区、K空间的问题
已经有11人回复
panjunxing
木虫 (著名写手)
小学生
- 应助: 0 (幼儿园)
- 金币: 2994.8
- 散金: 948
- 红花: 2
- 帖子: 1103
- 在线: 69.8小时
- 虫号: 579223
- 注册: 2008-07-04
- 性别: MM
- 专业: 高分子组装与超分子结构
您能帮我一下吗？
3楼2014-03-13 09:10:20