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

引用第1楼cem-uestc于2009-10-09 09:36发表的  : O<$j}?2  
对应书上边界媒质处理的 %mD{rG9  

引用第2楼wufamei0022于2009-10-09 10:00发表的  :
谢谢你的指教，还有一个小疑问： !,rF(pz  
我个人按自己的方法推导了一下，最后形式能一样但FE1和FH1的符号与书中CA和CP相反，FE2和FH2除了符号相反外，pi/wl所处位置不对。按照个人所想，感觉是这样的：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./(Pi*((Media(1,i)+Media(1,j))/WL+.5*(Media(3,i)+Media(3,j)))   ?Iij[CbU  
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)   cM4{ e^  
  FH2(i,j)=1./(Pi*((Media(2,i)+Media(2,j))/WL+.5*(Media(4,i)+Media(4,j)))   8rFP*K9  
....... }n#$p{e$i