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

1h"B-x  
本人是新手，刚接触FDTD，希望各位不吝指点 2-CK:)n/#  
>pYgF=J  
a&<<X:$Hy  
这是设置系数的部分 _8`|KY  
&,)tD62s  
%*********************************************************************** NdZ: 
7  
%    Updating coefficients for space region with nonpermeable media U9q*zP_jV  
%*********************************************************************** >9klh-f  
cP Y^Bf5)  
scfact=dt/muz/dx; 1hN! 2Y:  
AuCVpDH  
ca=(1.0-(dt*sig)/(2.0*epsz*eps))/(1.0+(dt*sig)/(2.0*epsz*eps)); [+1 i$d  
cb=scfact*(dt/epsz/eps/dx)/(1.0+(dt*sig)/(2.0*epsz*eps)); =wQ=`  
Cm&itG  
0'<S7?~|  
LSs={RD2+p  
qnlj~]NV  
这里只定义了CA和CB，没有CP和CQ，而且CB的定义比起文献中多了scfact（dt/muz/dx） g$97"d'  
.6hH}BM  
&0"`\~lA  
Z~R i%XG  
YUH/t
l  
而在计算磁场时 -_}EQ9Q  
%*********************************************************************** -Z@p   
%    Update magnetic fields WZr~Pb9  
%*********************************************************************** x{- caOH  
'SY&-<t(  
hy(1:ie)=hy(1:ie)+ez(2:ib)-ez(1:ie); ppYz~ {"r  
*!Dzst-J3  
^Ve^}|qPc  
F,bl>;{[{  
按照文献应该是hy(1:ie)=cp*hy(1:ie)+cq*(ez( .. (o|E@d  
A4^+p0@  

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

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

"ILWIzf.]  
sigh， 6>:~?gs  
程序这种东西，一般别人帮不了太多的。 2HXKz7da  
It4z9Gh  

引用第1楼zb5700于2008-03-26 11:20发表的  : X:lPWz!7{  
你加我qq吧，我详细告诉你，说来话长了 454839848