[求助]麻烦大家看看这个一维FDTD中参数的设置

scfact=dt/muz/dx; ~UJ_Rr54  
 ca=(1.0-(dt*sig)/(2.0*epsz*eps))/(1.0+(dt*sig)/(2.0*epsz*eps)); c8jq.y v  
cb=scfact*(dt/epsz/eps/dx)/(1.0+(dt*sig)/(2.0*epsz*eps)); Au/n|15->C  
ca的设置是按照书上公式写的，cb有点不一样呀，看不大懂，望 .. )Hy|K1  
nEm+cHHo?  

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

看不懂 啊  啊 啊啊

完整的程序在WORD文档里面。

据程序中 ez(2:ie)=ca*ez(2:ie)+cb*(hy(2:ie)-hy(1:ie-1));   e`zx#v  
                       hy(1:ie)=hy(1:ie)+ez(2:ib)-ez(1:ie);   [Hw  
O*xC}$OOn  
求解hy 时没有用到系数，可以猜想一下 求cb时，应该是对Ez 与 hy 做了一些处理，就像用 Eta0=sqrt( mu0/eps0) 对电场和磁场进行归一化处理。。 具体是怎么样做的，可能还是要靠自己用笔去推算一下了。。。。应该难不到楼主。。呵呵