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

do i=Iomin,Iomax H7{)"P]{f  
            do j=Jomin,Jomax kYbqb?  
               !下输出面     k }amSsE  
               tmp=i*sinto*cospo+j*sinto*sinpo+(Komin)*costo /e/%mo  
               tmp0=rmax-tmp 3K;V3pJ].  
#`Af  
               nn=int(n+0.5+2*tmp0) (*~'#k  
               f=n+0.5+2*tmp0-nn 'e&L53n  
               SvTd#>ke  
               !电场Ey b
suGZ  
               tmp=-0.25*(Ey(i,j,Komin-1)+Ey(i,j+1,Komin-1)+Ey(i,j,Komin)+Ey(i,j+1,Komin)) !wbO:py[8>  
               ux(nn-1)=ux(nn-1)+(1.0-f)*tmp I/2{I  
               ux(nn)=ux(nn)+(2*f-1.0)*tmp rra|}l4Y  
               ux(nn+1)=ux(nn+1)-f*tmp k4N_Pa$}\  
F,dPmR  
6d7E@}<  
这是我瞬态外推的一点程序，我感觉没什么问题，结果 .. ]rNM3@bVy  
p{@jM  

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

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