8.3 矩量法与高频近似
8.3 矩量法与高频近似
CEM 的收官一站处理两类"前两节不方便"的问题：开放空间的辐射与散射（FDTD/FEM 需要把周囲大块真空都网格化，代价高昂）交给矩量法——只在导体表面离散；目标尺寸数百个波长的"电大尺寸"问题（所有网格法都内存爆炸）交给高频近似——光学思维回归电磁。
MoM：把体问题缩成面问题
天线与散射问题的本质是积分方程：入射场在导体表面激励电流 J，J 又辐射场、反过来满足边界条件（表面切向总场为零）。矩量法的操作：
- 只在表面剖分（线段/面片），把连续的电流分布 J 展开为基函数组合 J ≈ ΣIₙ·fₙ（RWG 基函数是三角面片上的标准选择）；
- 选一组权函数（伽辽金法取权=基）对方程做内积，得到 Z·I = V：Z 是 n×n 的阻抗矩阵（Z_mn 刻画第 n 个基函数的场在第 m 个检验函数处的反应——互阻抗，第 7 章互耦的直接产物），V 是激励向量；
- 解出表面电流 I，一切想要的量（方向图、增益、RCS）由电流积分而得。
好处看得见：未知量从"三维体积"降到"二维表面"，天线问题的规模骤减；辐射条件自动满足（积分方程的格林函数已把"波跑到无穷远"编码在内），不需要 PML。代价也看得见：Z 是稠密矩阵（每个基函数与所有其他基函数都耦合），存储 O(n²)、直接解 O(n³)——这正是第 7 章互耦"人人都相关"的矩阵化身。万级未知量就让直接解吃力，于是催生了一整支快速算法：**快速多极子方法（FMM/MLFMA）**把远区耦合分层聚合-转移-发散，把复杂度压到 O(n log n)，百万未知量的飞机级 RCS 计算由此可行。
演练：一个 10 段细线半波振子的 MoM 缩微版思路（完整实现见 Harrington 经典教材；此处给出骨架代码）：
# 半波振子 MoM 的脉冲基-点匹配骨架（概念演示） import numpy as np N = 20; L = 0.5; dz = L/N; k = 2*np.pi # λ=1 z = (np.arange(N) - N/2 + 0.5) * dz Z = np.zeros((N, N)) for m in range(N): for n in range(N): r = abs(z[m] - z[n]) Z[m, n] = (np.exp(-1j*k*r)/r).real * dz * 0.01 if r > 1e-9 else 0.5 # 近项正则化 V = np.zeros(N); V[N//2] = 1.0 # 中点馈电 I = np.linalg.solve(Z, V) print("电流分布（对称性检验）:", np.round(np.abs(I)/np.max(np.abs(I)), 2)) # 输出：中点最大、两端趋零，验证半波驻波分布（第7章图3的数值版）
打印出的电流分布呈中点峰、两端零的对称钟形——第 7 章"半波振子正弦驻波"的数值自证。
高频近似：光学思维的回归
目标尺寸比波长大几个量级时（飞机在 X 波段长达数百 λ），场的局部行为趋近几何光学：波长趋于零，麦克斯韦方程退化为光线追踪。两条主力技术：
物理光学（PO）：表面电流近似为切向入射场的两倍（J = 2n̂×H_inc，点亮面），对表面电流积分得散射场。计算便宜、精度在大平缓曲面处好；边缘绕射不准，需 GTD/PTD 修正。
几何绕射理论（GTD）：把光线模型扩展"绕射射线"——边缘、尖顶、曲面爬波各自有渐近绕射系数。雷达特征工程（进气道、缝隙、边缘的散射热点定位）的主力工具。
近似方法与全波方法的分工是频率的函数：同样的无人机，L 波段（波长 30 cm）用 MoM/FMM 全波算，毫米波（波长 8 mm）只能 PO/GTD——"模型层级随电尺寸切换"是 CEM 工程的第一判断。混合方法（PO 修大平面 + MoM 修边缘）在商业软件中已是常规选项。
验证与可信度：仿真的"最后一公里"
任何 CEM 结果在进入设计决策前要过三道关：解析极限对照（退化为已知解析解的情形：球、平板、波导尺寸公式）；收敛性研究（网格/基函数加密后目标量稳定）；实验对标（暗室测量、S 参数实测）。行业黑话"garbage in, garbage out"在 CEM 里格外贴切：材料参数错、网格质量差、边界条件不当，仿真给出的漂亮彩色图只是"数字艺术品"。第 9 章的 EMC 排障会反复引用本节的结论——仿真指引设计，测量裁决真相。
易错点
MoM 解介质体想当然：表面积分方程对介质需要双电流（等效电/磁流）与多区域 formulation，纯导体代码套介质问题结果全错。
PO 用在小目标：电尺寸小于十几个波长时绕射主导，PO 误差不可控；"近似方法省时间"只在它的频段成立。
忘了谐振结构对网格的特殊要求：高 Q 谐振（第 6 章腔体）的场在谐振频率附近剧变，扫频步长与网格密度都要加倍小心——错过了谐振峰，滤波器仿真就白做。
补充辨析：方法论的"最后一公里"
三道收尾是非题：第一道："MoM 解出表面电流后，远场方向图是精确的"——电流在基函数近似意义上精确，方向图的误差来自电流误差与远场积分截断；报告增益时应同时给出电流收敛指标。第二道："PO 算大目标又快又准"——快是实情，"准"要看部位：镜面区好、边缘与阴影区差；工程惯例是 PO 打底、边缘绕射修正（PTD）补强，重要目标仍抽样全波核对。第三道："机器学习正在取代 CEM 数值方法"——目前的角色是"代理模型"（用少量全波样本训练快速预测器），适合优化迭代与初筛；真正的设计定型仍回到全波与实测，因果机制（麦克斯韦方程）没有让位于统计拟合。本章三节合起来的总功课是：为给定问题选对工具组合，并诚实标注每一步的误差来源——这份"工具选型 + 误差预算"的双重素养，比精通任何单一软件都稀缺，也是从本章走向第 9 章实战前最值得打包带走的行李。
数字感受收尾：万级未知量的稠密矩阵存储约 800 MB（复双精度）——MoM 的"稠密税"一眼可见；MLFMA 把同规模降到百 MB、迭代求解分钟级。而一个 50 波长的飞机在 X 波段的网格法需求以亿计单元——高频近似的存在不是偷懒，而是物理尺寸逼出来的必然分层。数字看懂，选型表就变成了常识。