矩量法（MoM）— 积分方程基础的电磁场解析手法
矩量法（MoM）的理论基础
概述 — 什么是MoM
矩量法与FEM和FDTD有什么区别呢？这些都是求解电磁场的方法吧？
好问题。MoM是基于积分方程的方法，仅以导体表面的电流分布为未知量。无需体积网格，对开放域问题有压倒性优势。电线天线和RCS（雷达截面积）计算中，MoM比FEM效率高得多。
无需体积网格，这样做岂不是太省力了！
正是如此。FEM必须对结构周围的空气进行3D体积网格划分。而且处理开放空间还需要吸收边界条件（PML等）。MoM只需离散导体表面，格林函数自动满足无穷远处的辐射条件。NEC2是一个免费的MoM代码，业余无线电爱好者也广泛使用。
NEC2就是那个业余无线电爱好者常用的吗？专业现场也能用吗？
NEC2是专门针对电线天线的薄线近似MoM代码。偶极子天线和八木天线的设计仍是首选。但它不支持面元素，所以对于贴片天线或包含筐体的复杂形状，需要用Altair Feko或Ansys HFSS的MoM求解器。本文将从MoM的理论基础到实现方法，系统地讲解。
与FEM、FDTD的本质区别
电磁场数值解析主要有三种方法。理解各自的定式化和适用领域是选择合适方法的关键。
| 方法 | 定式化 | 离散化对象 | 矩阵性质 | 优势领域 |
|---|---|---|---|---|
| MoM | 积分方程 | 仅导体表面 | 密矩阵（小规模） | 天线辐射、RCS、开放域 |
| FEM | 微分方程（弱形式） | 整个体积 | 稀疏矩阵（大规模） | 复杂形状、非均质媒质、谐振器 |
| FDTD | 微分方程（差分） | 体积全部（正交网格） | 无矩阵（显式） | 宽带瞬态响应、非线性 |
控制方程 — EFIE（电场积分方程）
想看MoM的数学式子。基础方程是什么样的？
基本方程是EFIE（电场积分方程）。它由完全导体表面上切向电场为零的边界条件推导。
完全导体（PEC）表面 $S$ 上，入射电场 $\mathbf{E}^{inc}$ 与表面电流 $\mathbf{J}(\mathbf{r}')$ 的关系由以下EFIE描述：
其中格林函数 $G(\mathbf{r}, \mathbf{r}')$ 表示自由空间中点源的响应：
看起来相当复杂呢……简单说是什么意思？
简单说就是「导体上任意位置的电流通过格林函数这个'传播权重'相互影响」。左边是入射波（外来的电磁波），右边是电流产生的散射场。MoM的任务是求解使两者在导体上平衡的电流分布 $\mathbf{J}$。
Galerkin法的离散化
这个积分方程怎样在计算机上求解呢？
未知的表面电流 $\mathbf{J}$ 用基函数（basis function）$\mathbf{f}_n$ 的线性组合近似。这是MoM的核心。
将未知电流用 $N$ 个基函数展开：
$I_n$ 是各基函数的系数（未知数）。将此展开式代入EFIE，用测试函数 $\mathbf{f}_m$ 进行内积（测试）。Galerkin法中测试函数与基函数相同：
这里 $\mathcal{L}$ 是EFIE的积分算子，$\langle \cdot, \cdot \rangle$ 表示内积。
这不就是FEM的Galerkin法吗！
正确。从加权残差法来说，与FEM框架一样。区别在于算子是积分而不是微分。所以MoM有时也叫"加权残差法应用于积分方程"。这是Harrington在1968年的著作中系统化的。
阻抗矩阵 [Z][I]=[V]
Galerkin法的测试后，连续的积分方程转化为以下矩阵方程：
各分量定义如下：
$[Z]$ 是阻抗矩阵，$[I]$ 是电流系数向量，$[V]$ 是激励向量是吧。与FEM的 $[K]\{u\} = \{F\}$ 结构很像！
完全同意。但有一个决定性的区别。FEM的刚度矩阵 $[K]$ 是稀疏矩阵（只有相邻单元之间有耦合），而MoM的阻抗矩阵 $[Z]$ 是密矩阵（完全矩阵）。因为格林函数使所有基函数互相耦合。
那就是说，如果有 $N$ 个未知数，就需要 $O(N^2)$ 的内存吗？
聪慧的观察。直接法求解需要 $O(N^2)$ 的内存和 $O(N^3)$ 的计算时间。所以MoM长期以来被看作"规模小但精度高"的方法。MLFMM（多层级快速多极法）打破了这个局限，可以将计算量削减到 $O(N \log N)$。数值解法章节会详细说明。
基函数的选择
基函数 $\mathbf{f}_n$ 具体是什么样的形状？
代表性的基函数整理在表中。根据问题形状选用不同的函数。
| 基函数 | 形状 | 适用对象 | 特性 |
|---|---|---|---|
| 脉冲函数 | 区间内常数 | 电线（1D） | 最简单。NEC2采用 |
| 三角函数 | 区间内线性补间 | 电线（1D） | 保证电流连续性 |
| RWG函数 | 三角形对上的向量 | 面元素（2D） | Rao-Wilton-Glisson。面MoM标准 |
| 高阶基函数 | 多项式 | 面元素（2D） | WIPL-D采用。用较少单元达到高精度 |
RWG函数是面MoM的标准呢。这是不是类似FEM的边单元？
很好的比喻。边单元保证切向分量的连续性一样，RWG函数通过相邻三角形的共享边保证法向电流连续性。自1982年Rao、Wilton、Glisson提出以来，一直是面MoM的事实标准。例如Altair Feko也以RWG基函数为基础。
Harrington开创的积分方程世界
矩量法体系是由Roger Harrington（雪城大学）在1968年出版的《Field Computation by Moment Methods》一书中确立的。将表面电流用基函数展开，通过测试函数的内积构建阻抗矩阵的方法，使得之前无法解析的任意形状天线的辐射、散射问题有了数值解的可能。这本书仅约110页，却改变了电磁场计算的历史，至今仍被广泛引用。Harrington常说"学生理解不了的理论，说明理论本身有问题"。
矩量法（MoM）的数值计算手法
密矩阵问题及其处理
刚才说密矩阵需要 $O(N^2)$ 内存，具体规模限制有多大呢？
粗略估计，$N = 10{,}000$ 个未知数时，$[Z]$ 是 $10^4 \times 10^4 = 10^8$ 个复数单元，内存约1.6 GB。$N = 100{,}000$ 时160 GB。通过LU分解的直接法需要 $O(N^3)$ 运算量，$N$ 增加时很快就无法计算。
| 未知数 $N$ | 矩阵规模 | 内存（复数双精度） | LU分解时间（估计） |
|---|---|---|---|
| 1,000 | $10^6$ | 16 MB | 数秒 |
| 10,000 | $10^8$ | 1.6 GB | 数分钟 |
| 100,000 | $10^{10}$ | 160 GB | 数天 |
| 1,000,000 | $10^{12}$ | 16 TB | 直接法不可行 |
16 TB……普通工作站根本不可能。大问题怎么解呢？
这就是MLFMM的用武之地。不显式构造矩阵，而是用 $O(N \log N)$ 快速近似矩阵向量乘积。与迭代法（GMRES等）结合，百万级未知数的问题也能在普通工作站上求解。
MLFMM（多层级快速多极法）
MLFMM（多层级快速多极法）是一种不显式建立矩阵，而通过多极展开近似远场基函数间相互作用的技术。
多极展开……有直观的图像吗？
比方说查询远方城市的人口时，不逐个计算，而是"东京约1400万人"这样汇总处理。近处的人（邻接网格）精确计算，远处的群体用汇总数据近似。这种"近精确、远近似"的分层结构就是MLFMM的本质。
| 求解方法 | 内存 | 矩阵向量乘积 | 总计算量 |
|---|---|---|---|
| 直接法（LU分解） | $O(N^2)$ | — | $O(N^3)$ |
| 迭代法（显式保存矩阵） | $O(N^2)$ | $O(N^2)$ | $O(N^2 \cdot n_{iter})$ |
| MLFMM + 迭代法 | $O(N \log N)$ | $O(N \log N)$ | $O(N \log N \cdot n_{iter})$ |
EFIE、MFIE、CFIE的应用选择
EFIE之外还有MFIE、CFIE，有什么区别呢？
三个积分方程各有特点：
- EFIE（电场IE）：可用于开放面和电线结构。但闭合导体在内部共振频率处解不唯一
- MFIE（磁场IE）：仅用于闭合导体。也存在内部共振问题
- CFIE（组合场IE）：EFIE和MFIE的线性组合 $\alpha \cdot \text{EFIE} + (1-\alpha) \cdot \text{MFIE}$。完全消除内部共振，与MLFMM组合时迭代收敛也优异。通常 $\alpha = 0.5$
闭合导体的RCS计算应该毫不犹豫地用CFIE吧。
完全同意。Altair Feko也默认采用CFIE。但电线天线这样的开放结构无法使用MFIE，只能用EFIE。Feko和NEC2会自动判断结构类型选择合适的方程。
奇异积分的处理
格林函数在 $|\mathbf{r} - \mathbf{r}'| \to 0$ 时发散。自己单元的积分怎么办？
好眼光！这就是"奇异积分"问题，是MoM实现中最困难的部分之一。主要处理方法有：
- 奇异部分提取法：把 $1/R$ 的奇异部分解析积分，其余的光滑部分用数值积分
- Duffy变换：用坐标变换将三角形单元的自项非奇异化
- 自适应高斯积分：根据观测点与积分域的距离增加或减少积分点数
实务中，近场单元（$R < 0.5\lambda$）采用高精度积分，远场单元用低阶高斯积分或MLFMM处理。
矩量法（MoM）的实务应用
MoM解析的实务流程
实际用MoM解析时要走哪些步骤呢？
典型的MoM解析流程如下：
- 形状建模：电线结构分段，面结构从CAD生成三角形网格
- 励振设置：馈电点（电压缝隙）、入射平面波、近场源等
- 网格合理性检验：确认段长、三角形边长不超过 $\lambda/10$
- 求解器设置：选择EFIE/CFIE、是否使用MLFMM、频率扫描范围
- 解析执行：阻抗矩阵构建 → 线性方程求解 → 电流分布获取
- 后处理：远场图案、增益、输入阻抗、RCS的计算
表面网格的最佳实践
知道 $\lambda/10$ 法则了，还要注意什么呢？
MoM的网格与FEM有不同的要点：
- 边长上限：最高解析频率处的 $\lambda/10$ 为目标。可放宽到 $\lambda/7$，但要注意精度
- 边角部：导体边缘电流集中，应在其周围细密化网格
- 长宽比：三角形长宽比推荐5以下。细长三角形会降低RWG函数精度
- 电线与面的连接：电线段端点应与面网格节点重合
- 馈电间隙：NEC2中电压源的馈电段长度影响输入阻抗。推荐 $\lambda/1000$ 左右
电线建模的注意事项
用NEC2分析偶极子天线时，常见的失败是什么？
NEC2的电线建模基于"薄线近似"，要求电线半径 $a$ 远小于段长 $\Delta$。具体来说：
- $a / \Delta < 1$：违反此条件会破坏薄线近似。理想情况 $a / \Delta < 0.01$
- 邻近电线间距：两条电线接近时，间距应超过电线半径的5倍
- 段长均分：电线通常等间距分段。电流变化大的部分也可不等分细化
- 接合角度：电线接合角度小于90度时容易产生电流不连续
原来FEM网格质量的指标与完全不同呢。电线半径与段长的比例容易忽视。
常见失败及其对策
| 现象 | 原因 | 对策 |
|---|---|---|
| 输入阻抗异常 | 馈电段太长 / 网格过粗 | 馈电区附近细密化网格 |
| 内部共振导致发散 | 闭合导体仅用EFIE | 改用CFIE（$\alpha = 0.5$） |
| MLFMM迭代不收敛 | 多极展开项数不足 | 增加展开次数、改进预处理 |
| 图案对称性破坏 | 网格不对称 | 重新生成对称网格 |
| 低频条件数恶化 | EFIE低频分解 | 使用环路-星形基函数 |
天线与筐体的影响 — MoM-PO混合法
设计车载或飞机天线时，不仅要考虑天线本体，还要考虑与大型导体筐体（车身、机体）结合后的辐射图案。筐体的衍射与反射会大幅改变天线增益。MoM（特别是MoM-PO混合法）精细处理天线周围细节，用物理光学近似（PO）处理远方大型筐体，实现了巨型模型在合理计算时间内的解析。Altair Feko的HYBRID求解器就是这一方法的代表实现，被广泛应用。
矩量法（MoM）的软件比较
NEC2 — 免费的MoM代码
先来详细了解NEC2。什么时候开发的？
NEC（数值电磁学代码）由Lawrence Livermore国立实验室在NASA资助下于1970年代开发。NEC-2在1981年发布并进入公有领域。虽然是40多年前的代码，但电线天线解析的精度至今仍被公认。
- 定式化：EFIE（电线）+ MFIE（面）混合
- 基函数：电线上的三角函数（正弦、余弦组合）
- 优势领域：偶极子、八木、对数周期、螺旋线等电线天线
- 限制：面贴片仅支持平面（不支持曲面）、不支持介质
- 许可证：NEC-2为公有领域。NEC-4需要LLNL许可
业余无线电爱好者用的GUI工具吧？
是的。代表作是4nec2（Windows、免费）和EZNEC（付费、更直观的UI）。不需要懂NEC卡片组指令（GW/EX/FR等）就能用GUI设计天线。业余无线电爱好者日常优化天线的指向性和驻波比都用这些工具。
Altair Feko
专业现场用什么工具？
专业MoM求解器中应用最广的是Altair Feko。源于南非斯坦伦博什大学，通过EM Software & Systems公司开发，2014年被Altair收购。
- MoM求解器：基于RWG基函数。支持EFIE/MFIE/CFIE
- MLFMM：大规模问题用的高速求解器内置
- 混合手法：MoM/MLFMM/PO/RL-GO/UTD/FEM自动混用。混合手法最多
- 优势领域：车载天线布局、飞机RCS、EMC/EMI、反射天线
- 特色：支持CMA（特征模分析）。可以物理理解天线基本模
其他支持MoM的工具
| 工具名称 | 开发者 | MoM实现 | 特性 |
|---|---|---|---|
| WIPL-D | WIPL-D公司（塞尔维亚） | 高阶基函数MoM | 单元数大幅减少，效率优异 |
| Ansys HFSS（IE-Solver） | Ansys Inc. | MoM（积分方程求解器） | 易与FEM切换、耦合 |
| CST Studio Suite | Dassault Systemes SIMULIA | MoM（积分方程求解器） | FDTD/FEM/MoM切换灵活 |
| COMSOL（RF模块） | COMSOL AB | BEM（MoM相当） | 多物理场耦合强 |
| 概念（CONCEPT-II） | 汉堡工科大 | MoM | EMC解析专用、学术用 |
功能比较矩阵
MoM解析到底选哪个工具？
按用途整理一下。
| 用途 | 推荐工具 | 理由 |
|---|---|---|
| 电线天线设计 | NEC2（4nec2） | 免费、高精度、经验丰富 |
| 车载/飞机天线布局 | Altair Feko | 大型筐体混合解析 |
| RCS计算 | Altair Feko / WIPL-D | CFIE + MLFMM 成熟度高 |
| 需要FEM联动 | Ansys HFSS（IE-Solver） | 同一环境中FEM↔MoM切换 |
| EMC解析（全波+电缆） | CST Studio Suite | 电缆线束MoM解析 |
| 学术研究、自定义实现 | 自写代码 / NEC2改造 | 算法验证的自由度大 |
MoM方法用处这么广，工具选择完全不一样呢。
没错。自称"MoM强"的工具，与Feko这样的MoM专业工具在大规模问题上的性能和稳定性有很大差异。建议在选型前用基准问题做比较。
NEC2与业余无线电社群
NEC2进入公有领域后，在全球业余无线电爱好者中迅速扩大使用范围。CQ杂志和美国电台中继联盟（ARRL）的文章中充满了"用NEC2优化的八木天线"设计案例。特别是4nec2这个Windows前端出现后，摆脱了文本卡片组的束缚，用户群体爆炸式增长。专业天线工程师中也有很多人表示"第一个天线设计就是用NEC2"。一个40多年前的代码至今仍在现役，这是Harrington理论和MoM普遍性的证明。
矩量法（MoM）的前沿研究
混合手法（MoM-PO、MoM-FEM）
MoM的前沿研究朝什么方向发展？
最实用且成效显著的是混合手法。结合不同解法各自的优势，处理单一方法无法应对的问题。
- MoM-PO：天线周围精细部分用MoM，大型筐体用物理光学（PO）近似。飞机RCS和散射解析中是标准
- MoM-FEM：导体部分用MoM，介质、吸收体用FEM。罩型天线的常用方法
- MoM-UTD/GTD：远方棱角衍射用几何光学（GTD/UTD）处理。进一步削减计算成本
Feko的HYBRID求解器就是这么做的吧！
完全同意。Feko在单个模型中混合MoM/MLFMM/PO/RL-GO/UTD/FEM。用户指定各区域的求解方法，软件自动在边界处理连接条件。飞机天线解析中，天线本体（MoM）+ 机翼表面（MLFMM）+ 机体远场（PO）三层混合是典型配置。
宽带MoM与模型阶数降低
频率扫描很慢的问题怎么解决？
有几个方法：
- AWE（渐近波形评价）：从某频点的阻抗矩阵及其导数，用Pade近似外推邻近频率的响应。频率点数可削减10~100倍
- CBFM（特征基函数法）：大型结构分解为子区域，各部分的特征基函数预计算。频率扫描计算量大幅削减
- CMA（特征模分析）：结构固有模式频率非依赖计算。用于物理理解和最优设计
与机器学习的融合
最近AI/ML热潮对MoM也有影响吗？
出现了几个受关注的研究方向：
- 代理模型：用神经网络学习MoM的输入输出关系。加速参数优化的反复计算
- 前处理学习：为迭代解法学习最优的预处理矩阵
- 网格优化：基于电流分布预测的自适应网格细化
不过，MoM基于物理规律的精度不能被ML代替。主要还是用ML来"加速MoM"这样的方向。
矩量法（MoM）的故障排查
内部共振问题
听说闭合导体的MoM解析有时会出问题…
现象：在特定频率，阻抗矩阵的条件数急剧恶化，求解发散
原因：闭合导体的EFIE（或MFIE）在导体内部共振频率处产生非唯一性。这是数学定式的人为产物，非物理共振。
对策：
- 改用CFIE（$\alpha \cdot \text{EFIE} + (1-\alpha) \cdot \text{MFIE}$，$\alpha = 0.5$ 推荐）
- CFIE能完全消除内部共振
- 开放结构（电线）不会产生内部共振
MLFMM收敛不良
MLFMM迭代不收敛时，查什么？
现象：GMRES/BiCGStab等迭代法在指定次数内不收敛
可能原因与对策：
- 多极展开项数不足 → Feko中增加 "Number of multipoles per box"
- 用了EFIE → 改用CFIE。CFIE的矩阵谱特性好，收敛快
- 预处理不当 → 试试块对角预处理或SAI（稀疏近似逆）
- 网格质量差 → 检查是否有极小或极大的三角形
- 迭代次数限制太小 → 增加上限（500~1000）
内存不足（密矩阵）
老师，MoM出了"Out of Memory"…
现象：阻抗矩阵分配时内存不足
对策（优先顺序）：
- 启用MLFMM：矩阵不显式构造，内存从 $O(N^2) → O(N \log N)$ 削减
- 粗化网格：$\lambda/10$ → $\lambda/7$（精度下降要注意）
- 混合手法：远场部分改用PO近似，MoM未知数削减
- 核外求解：矩阵退避到磁盘（速度下降）
- 并行计算：MPI分布式内存，矩阵分散到多个节点
MoM的全貌理解了！虽然基于积分方程而与FEM完全不同，最后还是归结为矩阵方程这点是一样的呢。先从NEC2的电线天线试试。
好主意！从半波长偶极子开始，与理论值（输入阻抗 73+j42.5 Ω）对比，再改变段数观察收敛性。有不明白的地方随时问。
"解析不匹配"时
- 先深呼吸 — 不要急着乱改设置，否则问题更复杂
- 做最小重现 — 比如半波长偶极子这样已知问题，验证同样设置是否正确
- 只改一个参数 — 网格密度、积分方程种类、频率三选一变更，分离影响
- 回到物理 — 用可视化检查电流分布是否有非物理的尖刺。对称结构的电流应对称
更多
错误