[讨论]三维瞬态外推远场

do i=Iomin,Iomax `5[d9z/6  
            do j=Jomin,Jomax (>x4X@b  
               !下输出面     !79^M  
               tmp=i*sinto*cospo+j*sinto*sinpo+(Komin)*costo $T
\z  
               tmp0=rmax-tmp 1@)8E`u  
#cqia0.H  
               nn=int(n+0.5+2*tmp0) $/aZ/O)F  
               f=n+0.5+2*tmp0-nn -;TqdL@  
               4eym$UWw  
               !电场Ey m ?a&XZ  
               tmp=-0.25*(Ey(i,j,Komin-1)+Ey(i,j+1,Komin-1)+Ey(i,j,Komin)+Ey(i,j+1,Komin)) -,Q !:  
               ux(nn-1)=ux(nn-1)+(1.0-f)*tmp [<Wo7G1s  
               ux(nn)=ux(nn)+(2*f-1.0)*tmp M ]dS>W%U  
               ux(nn+1)=ux(nn+1)-f*tmp *v[WJ"8@  
1Z?en  
krPwFp2[*  
这是我瞬态外推的一点程序，我感觉没什么问题，结果 .. +-=w`  
c|lo%[]R!  

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

电场的值取在面元的中心，你的电场平均值不应该这样取，这是取磁场的方法。我和你用了同样的推导公式，我也遇到了问题，想请教你：我的后向散射交叉极化图在时间步的最后一小段会出现较大的抖动，哪怕我把时间步增加了，前面的图形都很好，可是在最后还是会出现同样的问题，不知你遇到了没？