[求助]入射波如何加入？

这是个入射波源的程序： r'Wf4p^Xd  
***********************产生入射波*************************** P!Fykg  
}xC2~  
    do i=IncidentStart,IncidentEnd-1 mH\zSk  
      Hin(i)=Hin(i)-FH*(Ein(i+1)-Ein(i))  !1D FDTD for Hin. 7zCJ3p  
    end do Dx1(}D  
V<#KFm$>C  
    do i=IncidentStart+1,IncidentEnd-1 )1!<<;@0  
      Ein(i)=Ein(i)+FE*(Hin(i-1)-Hin(i))  !1D FDTD for Ein. nBp6uNK[  
    end do {!D(3~MI  
.i
y>N/u  
c     Set the unit source ~(v5p"]dj  
    Ein(Isource)=sin(OMIGA*N)      !入射波源 e1a\--  
    if(N.le.WL)Ein(Isource)=Ein(Isource)*.5*(1.-cos(OMIGA*N)) 5H.~pc
2y  
                                   !开关函数（升余弦函数） S>I` y]qlR  
    Ein(IncidentStart)=EBin(4)     !入射波吸收边界 D&F{0  
    EBin(4)=EBin(3) ~2\Sn-`  
    EBin(3)=Ein(IncidentStart+1) %:'G={G`QH  
    Ein(IncidentEnd)=EBin(2) !b_(|~7Lc  
    EBin(2)=EBin(1) [L2+k? *  
    EBin(1)=Ein(IncidentEnd-1) VIdKe&,  
lGr=I-=  
最后 .. 4ams~  
}O7sP^  

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

EBIN ~D3S01ecM  
只是临时变量 不用赋值  是用来放值的。 _?H3*!>3  
Ein(IncidentStart)=EBin(4)     !入射波吸收边界 g>h/|bw4  
    EBin(4)=EBin(3) @\w}p E  
    EBin(3)=Ein(IncidentStart+1) yn}Dj9(q  
    Ein(IncidentEnd)=EBin(2) :.ZWYze  
    EBin(2)=EBin(1) y@l&B+2ks  
    EBin(1)=Ein(IncidentEnd-1) )>{.t=#  
以上就是入射波的吸收边界啊

注意他的边界吸收条件。他的边界吸收条件，仅适用于背景是空气的1D计算。不是通用的吸收边界条件。

Ein(IncidentStart)=EBin(4)-------不是把EBin(4)的值赋给Ein(IncidentStart)的意思吗？可是它本身是什么啊？ JG^GEJ  
EBin(3)=Ein(IncidentStart+1) mQVlE__ub  
Ein(IncidentEnd)=EBin(2) ,1 H|{<  
EBin(2)=EBin(1) O+mEE>:w%  
EBin(1)=Ein(IncidentEnd-1)

还有用的是什么样的边界条件？我之前看的都是PML，这个不是吧？

只适用于1D真空环境下，完全吸收边界，推导一下就可以得到。

这个是辐射边界条件，葛德彪老师书上有现成的公式。 G:PcV_ihx  
其实，像这样的自由空间的入射波（一般尽可能把连接条件加在自由空间，故这种情况是最常见的），可以根据电磁波传播速度解析求出任意点的入射波（与入射源的相差由空间距离换算得到），这样就不需要做入射波的差分和边界条件了——更重要的是，入射方向与连接边界不平行时，这样要方便很多，但正如楼上所说：只适用于自由空间。

不好意思，我刚刚接触这个，还是不太清楚，在弱弱的问一句，下面这部分的介绍或推到在哪一章节里有？ [H2"z\\u  
Ein(IncidentStart)=EBin(4)     !入射波吸收边界 9"mcN3x:\e  
    EBin(4)=EBin(3) 3eS *U`_  
    EBin(3)=Ein(IncidentStart+1) 1Igo9rv  
    Ein(IncidentEnd)=EBin(2) ;b1*2-  
    EBin(2)=EBin(1) Lr`G. e  
    EBin(1)=Ein(IncidentEnd-1)

嘿嘿，那个书在办公室——我可背不下来——你看书中辐射边界条件那一节，最终可以用于编程序的公式就是你看到的这一段

楼猪引用的这段程序是入射波的产生代码 [TAW68f'  
m@hmu}qz-  
后续采用相关连接边界将这段程序计算好的入射波加如总场区就OK了啊   你这是断章取义