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

葛德彪教授一书中的fdtd迭代公式系数为什么是这种形式？   VIP7OHJh  
-~eJn'W  
  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)     mcz+P |  
2^t#6XBk/  
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) pA{ 5V9
 
pE{yv1Yg  
谢谢各位 .. < B_Vc:Q  
bifS 2>c  

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

葛德彪教授一书中的fdtd迭代公式系数为什么是这种形式？   n$m"]inX  
]43[6Im  
  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)     #s5 pz8v  
^GN8V-X4y  
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) QbYc[8-[  
5p?
!ni9  
谢谢各位大侠赐教！ 1|,Pq9  

2.5节中有，介质界面参数的选取。 "h$R ]~eG  

坚持 马上就能成功了

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

两点： X\A]"su  
1、介质界面参数的选取，根据Maxwell方程组积分方程的第一、二个方程 JieU9lA^&B  
2、电场和磁场的归一化处理。 )mD\d|7f  

我觉得是有些问题，尤其是程序中的K，z均是用的自由空间中的参数。 9<~,n1b>x  
在铝金属中，即，media（1）中迭代系数FE2分子应是sqrt(相对介电常数*相对磁导率)的。觉得这部分语句不严谨，虽然在铝中两者都为1，但这样程序就不具有通用性的 {Ot[WF  

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

写计算程序，还是自己根据公式写最好。因为看、理解别人的程序远比自己写痛苦 4?P%M"\Iv  

终于弄明白了啊