[求助]葛德彪教授书中附录二的程序中的fdtd迭代公式系数为什么是这种形式？

葛德彪教授一书中的fdtd迭代公式系数为什么是这种形式？   qbT].,?!U  
c\RDa|B,  
  FE1(i,j)=((Media(1,i)+Media(1,j))+Pi*.5*(Media(3,i)+Media(3,j))/WL)/((Media(1,i)+Media(1,j))-Pi*.5*(Media(3,i)+Media(3,j))/WL)     FE2(i,j)=1./((Media(1,i)+Media(1,j))-Pi*.5*(Media(3,i)+Media(3,j))/WL)     1miTE4;?  
?;i6eg17<  
FH1(i,j)=((Media(2,i)+Media(2,j))+Pi*.5*(Media(4,i)+Media(4,j))/WL)/((Media(2,i)+Media(2,j))-Pi*.5*(Media(4,i)+Media(4,j))/WL)     FH2(i,j)=1./((Media(2,i)+Media(2,j))-Pi*.5*(Media(4,i)+Media(4,j))/WL)  gZvl D  
OSq"q-Q  
谢谢各位 .. b,{?+
8  
Y<kvJb&1*  

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

葛德彪教授一书中的fdtd迭代公式系数为什么是这种形式？   |OhNQoTY  
u m(A3uQ  
  FE1(i,j)=((Media(1,i)+Media(1,j))+Pi*.5*(Media(3,i)+Media(3,j))/WL)/((Media(1,i)+Media(1,j))-Pi*.5*(Media(3,i)+Media(3,j))/WL)     FE2(i,j)=1./((Media(1,i)+Media(1,j))-Pi*.5*(Media(3,i)+Media(3,j))/WL)     uFL~^vz  
BTgL:  
FH1(i,j)=((Media(2,i)+Media(2,j))+Pi*.5*(Media(4,i)+Media(4,j))/WL)/((Media(2,i)+Media(2,j))-Pi*.5*(Media(4,i)+Media(4,j))/WL)     FH2(i,j)=1./((Media(2,i)+Media(2,j))-Pi*.5*(Media(4,i)+Media(4,j))/WL) %Mz(G-I.\  
?VO*s-G:J  
谢谢各位大侠赐教！ k9'%8(7M:  

2.5节中有，介质界面参数的选取。 APsd^J  

坚持 马上就能成功了

这里是不是有用归一化处理，其中第二式和第四式的系数与书中给的不是很相符，相差了一个1/2c，这里我取了deltaT=deltaX/2c，其它的问题已解决。

两点： 5 S$*YRp  
1、介质界面参数的选取，根据Maxwell方程组积分方程的第一、二个方程 S^]i  
2、电场和磁场的归一化处理。 ?{mFQ
  

我觉得是有些问题，尤其是程序中的K，z均是用的自由空间中的参数。 h lD0^8S  
在铝金属中，即，media（1）中迭代系数FE2分子应是sqrt(相对介电常数*相对磁导率)的。觉得这部分语句不严谨，虽然在铝中两者都为1，但这样程序就不具有通用性的 48nZ H=(Eh  

稍微纠正一下我上边的帖子，因为程序中对金属区和自由空间都采用的dt=ds/（2C0）C0是自由空间中的波速。归一化后是程序中的公式，归一化用的z0也是自由空间中的波阻抗。

写计算程序，还是自己根据公式写最好。因为看、理解别人的程序远比自己写痛苦 S)`%clN}J  

终于弄明白了啊