[转载]关于Sullivan Dennis的Matlab程序两点疑问

W$E!}~Ro  
本人是新手，刚接触FDTD，希望各位不吝指点 a@>P?N~LA9  
^U[c:Rz  
8OYw72&  
这是设置系数的部分 P5vxQR_*lc  
v61'fQ1Qg!  
%*********************************************************************** `%PU_;Y5Q  
%    Updating coefficients for space region with nonpermeable media +7HM7cw  
%*********************************************************************** b+Ly%&  
!N, Oe<  
scfact=dt/muz/dx; DOk(5gR  
6xFvu7L_c;  
ca=(1.0-(dt*sig)/(2.0*epsz*eps))/(1.0+(dt*sig)/(2.0*epsz*eps)); O@$wU9D<  
cb=scfact*(dt/epsz/eps/dx)/(1.0+(dt*sig)/(2.0*epsz*eps)); Vh%=JL sK  
Tc{n]TV  
Gw\-e;,  
 X`REhvT  
w_O3];  
这里只定义了CA和CB，没有CP和CQ，而且CB的定义比起文献中多了scfact（dt/muz/dx） +\@) 1  
OPYl#3I  
82mKI+9&"  
/hNZ7\|P  
$)V_oQSqn  
而在计算磁场时 9&FFp*'3  
%*********************************************************************** v6Vd V.BI  
%    Update magnetic fields 9ZKB,  
%*********************************************************************** ;r**`O  
m^6& !`CD  
hy(1:ie)=hy(1:ie)+ez(2:ib)-ez(1:ie); `V(zz  
H<gC{:S  
nhbCk6Y5LZ  
Rds_Cd C  
按照文献应该是hy(1:ie)=cp*hy(1:ie)+cq*(ez( .. s]D&):  
!!jitFHzb  

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

你加我qq吧，我详细告诉你，说来话长了 454839848

!7m )QNV  
sigh， Ab"mX0n  
程序这种东西，一般别人帮不了太多的。 "Yk3K^`1T.  
%LL*V|  

引用第1楼zb5700于2008-03-26 11:20发表的  : @YH<H
c  
你加我qq吧，我详细告诉你，说来话长了 454839848