[已解决]平面波激励源到底怎么加啊？

以下是 matlab编的 一部分代码 t gHXIr}3  
n3JSEu;J  
但是我是看得一头雾水 m$bDWxm#e  
请教各位高手， 在总场/散射场连接边界 设置入射波的的话  在三维情况下是不是对六个界面都要做修正啊 k2ZMDU  
q;H5S<]/  
还有 下面这段代码  加的到底是什么波源？ 怎么会有一维波的传播？？  怎么体现的入射波的方向？？ uw9w{3]0f  
多谢!!! =!P$[pN2  
tY6QhhuS:  
for j=2:jb-1 TWTh!  
    ez_inc(j)=ez_inc(j)+(dt/(epsz*ds)).*(hx_inc(j-1)-hx_inc(j));                 `& }C*i"  
end VfJdCg_  
oAxCI/  
for j=1:je                                                                             ^<QF*!  
    hx_inc(j)=hx_inc(j)+(dt/(muz*ds)).*(ez_inc(j)-ez_inc(j+1)); %Y<3v\`_  
end U~YjTjbd  
*K2fp=Ns  
ez_inc(1)=ey_low_m2; $' >|r]  
ey_low_m2=ey_low_m1;                             $BWA=2$  
ey_low_m1=ez_inc(2);                               +}.S:w_xQ  
53)*i\9&  
ez_inc(jb)=ey_high_m2; 8VR! Y0`e  
ey_high_m2=ey_high_m1;                                                    4dI=  
ey_high_m1=ez_inc(jb-1); _cD-E.E%  
x]<0Kq9K  
ez_inc(3)=sin(omega*dt*n); K
{[N.dX(  
#FsoK*F  
*4RL  
for i=is1:is2                                                                               fGMuml?[ e  
      for j=js1:js2-1                                                            $%0A#&DVh  
ey(i,j,ks1)=ey(i,j,ks1)-(dt/(epsz*ds)).*hx_inc(j); 8PwPI
%Pb  
ey(i,j,ks2+1)=ey(i,j,ks2+1)+(dt/(epsz*ds)).*hx_inc(j); WTQd}f  
      end %~^:[@xa*  
end -JfO} DRI  
G.sf>
.[  
,NoWAmv  
for i=is1:is2                                                                   anxwK47  
    for k=ks1:ks2 V( SRw  
hx(i,js1-1,k)=hx(i,js1-1,k)+(dt/(muz*ds))*ez_inc(js1); gaxxB]8  
hx(i,js2,k)=hx(i,js2,k)-(dt/(muz*ds))*ez_inc(js2); TM^.y Y  
    end (`FY{]Wz!  
end eCX
w8  
(G`O[JF  
for j=js1:js2                     &nb .. pm k;5 d  
NJ$c0CNy  

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

总程序是什么，是否可以上传呢？

这个应该是总场/散射场加入射平面波的方法 pQ yH`  
4&+lc*  
程序中应该是加的垂直入射的平面波

程序时 etynqk老兄上传的那个 三维总场/散射场matlab程序  论坛内有 我权限不够还上传不了 sigh~ W+D{4:  
Nvj0MD{ X  
这个问题我大体已经想明白 此处加的是一个沿着y方向传播的平面波 i.G"21M  
'(!
U5j  
但程序有几个细节上的错误  设置散射球时  修改的迭代参数应该为 caez cbex caey cbey caez cbez。 再就是球半径应该是 格点数乘以ds

恩 对 加的是一个沿y方向传播的 平面波 4]mAV\1  
<n{-&;>  
但是我还请教一下    加平面波时 是不是要对 三维的六个连接面都要做 field correction?? (H-}z`sy/@  
4+'d">+|  
再就是 这个平面波怎么会是 只有 ez 和hx 呢？ 在 如y=js1界面上的切向分量不是 ex和ez吗？ 这里怎么只对ez和hx做了 correction？？