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

这是个入射波源的程序： ex|kD*=  
***********************产生入射波*************************** +p[~hM6?  
6 %=BYDF  
    do i=IncidentStart,IncidentEnd-1 >u4e:/5]  
      Hin(i)=Hin(i)-FH*(Ein(i+1)-Ein(i))  !1D FDTD for Hin. $?s^HKF~  
    end do :rj78_
e9  
7'8O*EoB'  
    do i=IncidentStart+1,IncidentEnd-1 :D,YR(])  
      Ein(i)=Ein(i)+FE*(Hin(i-1)-Hin(i))  !1D FDTD for Ein. ew"Fr1UGYZ  
    end do 9F-k:hD |  
UP`q6]P  
c     Set the unit source {c1qC zM4  
    Ein(Isource)=sin(OMIGA*N)      !入射波源 ~"wD4Ue  
    if(N.le.WL)Ein(Isource)=Ein(Isource)*.5*(1.-cos(OMIGA*N)) [a`i{(!  
                                   !开关函数（升余弦函数） E??%)q  
    Ein(IncidentStart)=EBin(4)     !入射波吸收边界 p']AXJ`Z  
    EBin(4)=EBin(3) %*gO<U4L]  
    EBin(3)=Ein(IncidentStart+1) F#PJ+W*h  
    Ein(IncidentEnd)=EBin(2) D!P?sq_5r  
    EBin(2)=EBin(1) 4sAshrUf  
    EBin(1)=Ein(IncidentEnd-1) |")x1'M  
t#q>U%!  
最后 .. nx4aGS"F:  
K$ &wO.  

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

EBIN !BX62j\?  
只是临时变量 不用赋值  是用来放值的。 P&GZe/6Y  
Ein(IncidentStart)=EBin(4)     !入射波吸收边界 V!P3CNK  
    EBin(4)=EBin(3) `#f=&S?k  
    EBin(3)=Ein(IncidentStart+1) AH`15k_i  
    Ein(IncidentEnd)=EBin(2) Z C93C7lJ  
    EBin(2)=EBin(1) |FFC8R%@]u  
    EBin(1)=Ein(IncidentEnd-1) 2{CSH_"Z7  
以上就是入射波的吸收边界啊

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

Ein(IncidentStart)=EBin(4)-------不是把EBin(4)的值赋给Ein(IncidentStart)的意思吗？可是它本身是什么啊？ ?Rt1CDu  
EBin(3)=Ein(IncidentStart+1) ^A11h6I  
Ein(IncidentEnd)=EBin(2) ^p"4)6p-W  
EBin(2)=EBin(1) J-, H6u  
EBin(1)=Ein(IncidentEnd-1)

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

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

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

不好意思，我刚刚接触这个，还是不太清楚，在弱弱的问一句，下面这部分的介绍或推到在哪一章节里有？ kus}WJ  
Ein(IncidentStart)=EBin(4)     !入射波吸收边界 _k2w(ew?  
    EBin(4)=EBin(3) y3IWfiz>/d  
    EBin(3)=Ein(IncidentStart+1) 6)[<)?A.[  
    Ein(IncidentEnd)=EBin(2) @6&JR<g*t  
    EBin(2)=EBin(1) qn"K
9k  
    EBin(1)=Ein(IncidentEnd-1)

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

楼猪引用的这段程序是入射波的产生代码 qL 5>o>J  
p/olCmHD)  
后续采用相关连接边界将这段程序计算好的入射波加如总场区就OK了啊   你这是断章取义