| 查看: 1159 | 回复: 3 | ||
[求助]
帮忙解释一下这段fortran程序
|
|
INTEGER(2) STATUS DO I=1, N DO J=1, N PH(I,J)=0.0 FT(I,J)=0.0 END DO END DO NHTM=2**NUHTM ! 2^NUHTM为计算子孔径光斑的付氏变换点阵数 K=NHTM/2+1 M=INT(NSA/2) DO 200 L=1, NH I0=IC(L) J0=JC(L) DO 100 I=-M, M I1=I0+I DO 100 J=-M, M J1=J0+J FT(K+I,K+J)=CMPLX(XR(I1,J1),XI(I1,J1)) 100 CONTINUE NU0=NUHTM CALL FFTDB(NU0) I0=IH(L)-K !IH(M),JH(M): H-S子孔径中心座标 J0=JH(L)-K DO I=1, NHTM I1=I0+I DO J=1, NHTM J1=J0+J X=FT(I,J) PH(I1,J1)=PH(I1,J1)+X*CONJG(X) END DO END DO do i=1, nhtm do j=1, nhtm ft(i,j)=0.0 end do end do 200 CONTINUE R=PI/FLOAT(NHTM) M=K-1 DO 500 L=1, NH I0=IH(L) J0=JH(L) X1=0.0 Y1=0.0 U=0.0 DO I=-M+1, M I1=I0+I DO J=-M+1, M J1=J0+J V=PH(I1,J1) U=U+V X1=X1+I*V Y1=Y1+J*V END DO END DO V=R/U IF(U.NE.0.0) THEN TX(L)=V*X1-TX0(L) TY(L)=V*Y1-TY0(L) END IF 500 continue IF(IWS.EQ.0) THEN X1=0.0 Y1=0.0 DO 600 L=1, NH X1=X1+TX(L) Y1=Y1+TY(L) 600 CONTINUE X1=X1/NH Y1=Y1/NH DO 800 J=1, NH TX(J)=TX(J)-X1 800 TY(J)=TY(J)-Y1 TXM=X1-TXM0 TYM=Y1-TYM0 END IF call graph(2) STATUS = SETCOLOR(12) is=nccd/ns !ns,光束直径上的子孔径数;NCCD为传感器点阵数 js=is/2 DO I=1, NH ii=ih(i)-js+356 jj=jh(i)-js+310 STATUS = RECTANGLE($GBORDER,II,JJ,II+IS,JJ+IS) end do return END |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有7人回复
-
这样的filecode谁见过
已经有18人回复
-
关于代码变化问题，想知道的进来
已经有17人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求高人帮忙看一下我这个Fortran程序的错误出在哪里？
已经有12人回复
- 一个FORTRAN90程序，越界了，希望高手帮助
已经有5人回复
- 有会看fortran程序的，能帮忙看看这个程序吗？
已经有10人回复
- 哪位帮忙写份Fortran程序？
已经有10人回复
- 有关fortran的一次而问题，希望大家能帮帮忙，谢谢
已经有4人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】读文本数据出现段错误[fortran]
已经有14人回复
- 【求助】用fortran怎么描述下面这样的情况，我是新手，请大家帮忙
已经有8人回复
- 【求助】下面这个公式怎么用FORTRAN实现
已经有7人回复
- 【求助】请问一下fortran的可视化编程
已经有7人回复
- 【求助】求助个Fortran计算，应该不难，半天搞不定，请高手帮忙
已经有11人回复
janecool
铁杆木虫 (著名写手)
- 物理EPI: 1
- 应助: 54 (初中生)
- 贵宾: 0.006
- 金币: 5996.2
- 散金: 342
- 红花: 12
- 帖子: 2070
- 在线: 626.6小时
- 虫号: 830148
- 注册: 2009-08-18
- 性别: GG
- 专业: 粒子物理学和场论
2楼2012-04-20 21:16:52
yiluwangxi
木虫 (著名写手)
- 应助: 147 (高中生)
- 金币: 2919.1
- 散金: 386
- 红花: 22
- 帖子: 1939
- 在线: 1323.7小时
- 虫号: 829297
- 注册: 2009-08-17
- 性别: GG
- 专业: 声学
3楼2012-04-20 23:19:16
未注册用户
木虫 (职业作家)
- 应助: 50 (小学生)
- 金币: 3550.1
- 散金: 928
- 红花: 8
- 帖子: 3340
- 在线: 546.2小时
- 虫号: 1227538
- 注册: 2011-03-10
- 性别: GG
- 专业: 光学
4楼2012-04-21 09:25:16
