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

引用第1楼cem-uestc于2009-10-09 09:36发表的  : se2+X>@>  
对应书上边界媒质处理的 iV(B0z  

引用第2楼wufamei0022于2009-10-09 10:00发表的  :
谢谢你的指教，还有一个小疑问： "|.>pD#0&  
我个人按自己的方法推导了一下，最后形式能一样但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)))   ZgZ}^x  
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)   e
l;^cMY  
  FH2(i,j)=1./(Pi*((Media(2,i)+Media(2,j))/WL+.5*(Media(4,i)+Media(4,j)))   mQ('X~l  
....... m:/nw,