矩量法求解电磁散射.docVIP

目录 第1章 格林函数 1 1.1 位函数 1 1.2 自由空间格林函数 2 1.2.1 二维情况 2 1.2.2 三维情况 3 1.3 半空间格林函数 4 第2章 无限长柱散射－CFIE 1 2.1 理想导体柱散射－TMz极化 1 2.2 理想导体柱散射－TEz极化 3 2.3 均匀介质柱散射－TMz极化 5 2.4 均匀介质柱散射－TEz极化 8 第3章 理想导体目标散射－CFIE 12 3.1 表面积分方程 12 3.2 线性方程组 12 3.2.1 RWG基函数和试函数 12 3.2.2 矩阵元素计算 13 3.2.3 奇异性处理 16 3.3 计算散射场 17 3.4 振荡核在平面多边形上的积分 19 3.5 数值算例 21 第4章 均匀介质目标散射－PMCHWT 23 4.1 表面积分方程 23 4.2 线性方程组 24 4.2.1 矩阵元素计算 24 4.2.2 奇异性处理 25 4.3 计算散射场 25 4.4 数值算例 26 第5章 非均匀介质目标散射－VIE 27 5.1 体积分方程 27 第6章 快速多极子 28 6.1 引言 28 6.2 快速多极子 29 6.3 树型结构 31 6.4 球谐函数 32 6.5 多层快速多极子 33  格林函数 位函数 对于时谐场，将麦克斯韦方程分为只有电型源和只有磁型源的两部分 ▽× Ee = ? jωμHe ▽× Em = ? Jm ? jωμHm ▽× He = Je + jωεEe ▽× Hm = jωεEm ▽? De = qe ▽? Dm = 0 ▽? Be = 0 ▽? Bm = qm 其中E表示电场，H表示磁场，D表示电通量，B表示磁通量，Je表示电流，Jm表示磁流。 由于▽? Be = 0，无散场可以用矢量的旋度表示，引入矢量位A，使得Be = ▽× A，可以得到▽× Ee = ? jωμHe = ? jωBe = ? jω▽× A，因此▽× (Ee + jωA) = 0。无旋场可以用标量的梯度表示，引入标量位φe，则Ee + jωA = ? ▽φe，因此 Ee = ? jωA ? ▽φe (1.1.1) 于是▽× He = Je + jωεEe = Je + jωε(? jωA ? ▽φe)，将▽× He = μ-1▽× (▽× A) = μ-1[▽(▽? A) ? ▽2A]代入，整理得 ▽2A + k2A = ? μJe + ▽(▽? A + jωμεφe) (1.1.2) 不妨令▽? A + jωμεφe = 0 (Lorentz Gauge Condition)，于是有 ▽2A + k2A = ? μJe Ee = ? jωA ? j(ωμε)-1▽(▽? A) He = μ-1▽× A 由于▽? Dm = 0，无散场可以用矢量的旋度表示，引入矢量位F，使得Dm = ? ▽× F，可以得到▽× Hm = jωεEm = jωDm = ? jω▽× F，因此▽× (Hm + jωF) = 0。无旋场可以用标量的梯度表示，引入标量位φm，则Hm + jωF = ? ▽φm，因此 Hm = ? jωF ? ▽φm (1.1.3) 于是▽× Em = ? Jm ? jωμHm = ? Jm ? jωμ(? jωF ? ▽φm)，将▽× Em = ? ε-1▽× (▽× F) = ? ε-1[▽(▽? F) ? ▽2F]代入，整理得 ▽2F + k2F = ? εJm + ▽(▽? F + jωμεφm) (1.1.4) 令▽? F + jωμεφm = 0 (Gauge Condition)，于是有 ▽2F + k2F = ? εJm Hm = ? jωF ? j(ωμε)-1▽(▽? F) Em = ? ε-1▽× F 由于时谐场是线性的，因此对于既有电型源也有磁型源的时谐场，可以将上述的结论叠加得到。于是有 ▽2A + k2A = ? μJe (1.1.5) ▽2F + k2F = ? εJm (1.1.6) E = ? jωA ? j(ωμε)-1▽(▽? A) ? ε-1▽× F (1.1.7) H = ? jωF ? j(ωμε)-1▽(▽? F) + μ-1▽× A (1.1.8) 在自由空间中，矢量位A满足方程▽2A + k2A = ? μJe，矢量位F满足方程▽2F + k2F = ? εJm。因此，可以将A和F分别理解为源Je和源Jm产生的场。不妨设自由空间中点源产生的场为g(r)，r为观察点到点源的距离，则g(r)满足方程 (▽2 + k2)g(r) = ? δ(r) (1.2.1