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

以下是 matlab编的 一部分代码 >_jT.d  

$5jQm,V$K  
但是我是看得一头雾水 |`1lCyV\tE  
请教各位高手， 在总场/散射场连接边界 设置入射波的的话  在三维情况下是不是对六个界面都要做修正啊 G<#9`  
casva;  
还有 下面这段代码  加的到底是什么波源？ 怎么会有一维波的传播？？  怎么体现的入射波的方向？？ mFw`LvH?*  
多谢!!! SS[jk  
`|Pfa  
for j=2:jb-1 C}<j8a?  
    ez_inc(j)=ez_inc(j)+(dt/(epsz*ds)).*(hx_inc(j-1)-hx_inc(j));                 (, /`*GC  
end 7+0Kg'^+n  
7IB<0  
for j=1:je                                                                             (/Lo44wT  
    hx_inc(j)=hx_inc(j)+(dt/(muz*ds)).*(ez_inc(j)-ez_inc(j+1)); +3dWnBg?  
end yl/-!  
{`2! 3= "  
ez_inc(1)=ey_low_m2; brx 7hI  
ey_low_m2=ey_low_m1;                             y[:\kI  
ey_low_m1=ez_inc(2);                               I`^ 7Bk.r  
N DZ :`D  
ez_inc(jb)=ey_high_m2;  $A]2Iw!&  
ey_high_m2=ey_high_m1;                                                    HT0VdvLw  
ey_high_m1=ez_inc(jb-1);  S<#>g s4  
dQT A^m  
ez_inc(3)=sin(omega*dt*n); {VtmQU?cJ  
kU{+@MA;  
Mrysy)x  
for i=is1:is2                                                                               [-5%[ty9X  
      for j=js1:js2-1                                                            OC34@YUj[  
ey(i,j,ks1)=ey(i,j,ks1)-(dt/(epsz*ds)).*hx_inc(j); AiKja>Fl<  
ey(i,j,ks2+1)=ey(i,j,ks2+1)+(dt/(epsz*ds)).*hx_inc(j); v._Q XcE  
      end Ezsb'cUa(  
end )c!7V)z  
3eT5~Lbs  
]CHO5'%,$  
for i=is1:is2                                                                   wCw-EGLR  
    for k=ks1:ks2 }%Mj`Bh  
hx(i,js1-1,k)=hx(i,js1-1,k)+(dt/(muz*ds))*ez_inc(js1); tIn dve  
hx(i,js2,k)=hx(i,js2,k)-(dt/(muz*ds))*ez_inc(js2); & mwQj<Z  
    end R+{QZ'K.qg  
end hEUS&`K  
)r2$!(NQ  
for j=js1:js2                     &nb .. h,%`*Qg6  
@`R#t3)8JP  

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

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

这个应该是总场/散射场加入射平面波的方法 v7-z<'?s~  
{7d(B1[1  
程序中应该是加的垂直入射的平面波

程序时 etynqk老兄上传的那个 三维总场/散射场matlab程序  论坛内有 我权限不够还上传不了 sigh~
MS5X#B  
,I|3.4z  
这个问题我大体已经想明白 此处加的是一个沿着y方向传播的平面波 ,"en7  
6(pa2  
但程序有几个细节上的错误  设置散射球时  修改的迭代参数应该为 caez cbex caey cbey caez cbez。 再就是球半径应该是 格点数乘以ds

恩 对 加的是一个沿y方向传播的 平面波 JnE\z*
NB  
3g79/w  
但是我还请教一下    加平面波时 是不是要对 三维的六个连接面都要做 field correction?? Wf!u?nH.5  
J;DTh ]z?:  
再就是 这个平面波怎么会是 只有 ez 和hx 呢？ 在 如y=js1界面上的切向分量不是 ex和ez吗？ 这里怎么只对ez和hx做了 correction？？