[求助]从一个网站抄来的简单程序(新手看不懂)大家给解释一下,谢谢

考虑一维情况下，设TEM波沿z方向传播, Maxwell方程为： %)@(Tye -  
                   yeh adm\  
如果以图1所示的方式将计算空间进行离散可以得到    qSs^}eN  

FDTD差分公式为:

    式中电场E和磁场H在相隔1/2的时间步取样。而在实际程序中可以使用上次计算的结果表示相隔1/2时间步。于是根据（2）式，我们开始动手编制第一个FDTD计算程序。

如图1所示的划分方法，将计算空间划分为100个空间步，电场E和磁场H各相距半个空间步，电场位于元胞边界，磁场位于元胞中心，这样我们编制的程序如下：

real,allocatable::  Ex(:),Hy(:) G"R>aw  
real c0,Eps0,Mu0 }u$aPS<$!  
  KPvYq?F>4  
Imax=100 YSUH*i/%  
Imin=-100 6je%LHhL  
  t 1'or  
~\(>m=|C:H  
 allocate (Ex(Imin:Imax)) ~k_zMU-1  
allocate (Hy(Imin:Imax-1)) dg@/HLZ  
      61;5Yo  
c0=3.e8 [Z&s0f1Qb  
Eps0=8.854e-12 m&vYZ3vK[  
Mu0=1.256637e-6 mI*>7?  
{xx;zjt%}}  
derta=3e-1 wOl-iN=  
dt=derta/c0 h 7P?n.K  
      ;DWp>jgy  
t=5 RIpq/^Th  
t0=10 *j;r|P;g  
Ex=0. E;e2{@SX2K  
Hy=0. =1B&d[3;  
Y:Tt$EQ  
open(1,file='ex') ].+G-<.:  
open(2,file='z_cell') bI0+J)  
  Bc,z]  
 do n=1,150 Xg)FIaw]eT  
   do  i=Imin+1,Imax-1 wAwH8xLU  
    if (i.eq.-20) then o sH,(\4_  
      Ex(i)=exp(-((n-1)*1.-t0)**2/t**2) "f$A0RL  
    else tTxo:+xg  
          Ex(i)=Ex(i)+(Hy(i-1)-Hy(i))*dt/derta/Eps0 ,#FH8%Yf  
     endif            B=xZkc  
  enddo ^bw~$*"j
#  
  ju6_L<  
  do  i=Imin,Imax-1 -5*;J&.  
    if (i.ge.Imin.and.i.le.Imax-1) then cB'4{R@e  
X-^Oz@.>  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

real,allocatable::  Ex(:),Hy(:)  %定义变量 RT9@&5>il  
real c0,Eps0,Mu0 @x4Dt&:"  
$r_gFv  
Imax=100   %设置网格数目 4[wP$  
Imin=-100 {F[
Xe_=#"  
QI<3N  
sP7(1)\  
allocate (Ex(Imin:Imax)) %分配变量空间 :Tl?yGF  
allocate (Hy(Imin:Imax-1)) =)LpMTz  
     64@s|m*  
c0=3.e8 %常系数定义 >i<-rO>kN  
Eps0=8.854e-12 c^%k1pae(  
Mu0=1.256637e-6 RPQ)0.O7  
>kT~X ,o  
derta=3e-1 %空间和时间步长定义 m+ YgfR  
dt=derta/c0 9IV WbJ  
     <hdCO< 0(  
t=5%激励波形参数设置 b*i+uV?  
t0=10 Ckp=d  
Ex=0. NST6pu\,U  
Hy=0. ^DOcw@Z6HC  
/0(KKZ)  
open(1,file='ex') zKr(Gt8  
open(2,file='z_cell') bHM .&4G  
ualtIHXK)  
do n=1,150 %总共时间步数 /zKuVaC  
   do  i=Imin+1,Imax-1 hW&UG#PY>  
    if (i.eq.-20) then &$f?XdZ7  
      Ex(i)=exp(-((n-1)*1.-t0)**2/t**2)  %加激励源 h'+ swPh  
    else 4{pa`o3  
          Ex(i)=Ex(i)+(Hy(i-1)-Hy(i))*dt/derta/Eps0 %更新电场值 DEhA8.v  
     endif           \8ZNXCP  
  enddo I 3PnyNZ  
)2S\:&x  
  do  i=Imin,Imax-1 }z#8vE;  
    if (i.ge.Imin.and.i.le.Imax-1) then :W"ITY(  
      Hy(i)=Hy(i)-(Ex(i+1)-Ex(i))*dt/derta/Mu0%更新磁场 ' r/1+.  
    endif 3[4]
G@  
  enddo he#iWD'  
enddo     cCIEG e6  
write(1,*) (Ex(i),i=Imin,Imax) u1l#k60  
write(2,*) (i,i=Imin,Imax) 0Og =H79<  
end

不过我不觉得你看这个程序对你有太大的帮助，还不如在本论坛找找一些现有的源代码，感觉比这个规范多了。

哦,谢谢,刚开始学,还不知道怎么叫规范,我目前现在的阶段,程序在我眼里都一样,看不出区别~~版主不妨给我推荐几个三维的程序,谢谢

论坛上有三维的FDTD程序

http://bbs.rfeda.cn/read-htm-tid-16745-fpage-3.html Bm1yBKjO  
KD=T04v  
这是我以前贴的Susan的代码，你看一下。这是比较规范的。

引用第5楼gwzhao于2008-11-17 11:00发表的 : Fn yA;,*  
http://bbs.rfeda.cn/read-htm-tid-16745-fpage-3.html
这是我以前贴的Susan的代码，你看一下。这是比较规范的。

P.=&:ay7?  
谢谢斑竹,目前正在看葛老师的课本,看的头都大了,顺便问一下,我算电磁波经缝隙与电缆耦合的话用MATLAB好,还是用FORTRAN好,给点建议,谢谢

这要看你对编程语言的爱好、选择了 .R#p<"$I  
算法的实现与编程语言无关

引用第7楼cem-uestc于2008-11-20 19:15发表的 : H[wJ; l  
这要看你对编程语言的爱好、选择了 Cg |_) _w  
算法的实现与编程语言无关

cpF\^[D  
WZ"NG|  
这个吧,我学过C,有点基础,要用fortran得从头学,所以想偷点懒用MATLAB做,不过好像程序也看不懂啊 tO~DA>R  
7[rn ,8@  
有什么MATLAB与FDTD相关的编程资料吗?给推荐一下,谢谢(我们学校没有这方面的书)

matlab书籍很多很多的啊，网上都可以找到很多电子书的。