Matlab+cpp矩量法代码演示

0前言

在上一篇博客中，我们详细介绍了矩量法（MoM）的原理及其数值求解过程。在这篇文章中，我们将以电场积分方程（EFIE）为例，介绍基于RWG函数的矩量法C++代码。
在正式解析代码之前，要先提一下关于矢量运算库函数的问题。
在求解EFIE方程的MoM方法中，涉及到了大量的矢量和矩阵运算过程，为了提高我们的编程效率，推荐选择一款成熟的矢量运算库。
这里，给大家推荐armadillo C++。这个运算库的使用方式非常类似Matlab，有matlab经验的人非常容易上手，我这里就不再过多介绍了。另外，我在个人网盘中上传了armadillo的配置和安装教程。感兴趣的朋友可以自行下载，链接如下：https://pan.baidu.com/s/1M0pZF_d56CNyhdF5mK9wZQ。

1三维目标几何剖分与网格信息处理

从上一张的介绍中不难发现，利用MoM求解EFIE，关键难点在于计算并填充阻抗矩阵方程。说白了，就是计算不同RWG函数之间基于电场积分算子联系起来的相互作用效应。那么问题来了，我们怎么将针对一个三维目标定义与其相对应的RWG函数呢？
根据RWG函数的定义，它的本质是依托一对三角面片的公共边定义的。因此，拿到一个三维目标，我们要做的第一步就是对其进行表面三角剖分，将其分解为一系列三角面片的组合。通常，这样的表面三角剖分工作依托一些现有的成熟工具（如FEKO、Ansys ICEM CFD）等就可以开展。
得到目标剖分后的三角面片网格信息之后，我们要对其进行一系列的处理，如定义三角