电磁仿真原理——4. 矩量法（Moment Methods）

矩量法介绍

矩量法的经典例子

1. 如下图所示是带有电荷的一段有限直导线，导线上的电荷分布是不均匀的（因为电荷存在同性互斥的作用，导致两端电荷密度大，中间电荷密度小），空间中任一点电位定义为Φ_e (r )，r为场点的位置，而r^，定义为源点。

2. 将其分割成很多小区域，每一段的长度为Δx，那么原有的源点与场点的距离可表示为如下（为了分析简单，忽略z轴方向）：
   
3. 定义一个近似的解：首先假定电荷的密度在每一小段Δx上是均匀的，那么可以将电荷密度定义为脉冲函数的累加：
   
4. u_n(x)是基函数，在变分法的介绍中，我们定义的是全域的基函数，但是全域的基函数不好找（满足完备和边界两个条件），而且对于通用的算法的时候，我们是没有办法知道需要计算什么模型的，不能保证这个基函数适用于所有模型。
   所以我们需要定义一个子域的基函数，例如上述的例子，将其划分为不同的区域，每个区域中的边界条件，和场量变化是相对简单的。但这个基函数只适用于这个子域。
   
5. 定义电位为1（导线内部和外部电位是相等的），对公式进行变形，利用了脉冲函数的特性将累加符号移到外部：
   
   
6. 其中为了处理奇异点的问题（分母有可能为0），如果源点和场点位于不同的子段的时候分母肯定不为0，但是如果源点和场点位于同一个子段，则存在为0的情况，将半径固定为a，场点固定在导体表面，则距离改写为：
   
7. 将累加符号展开，然后因为脉冲函数的定义，所以积分项会有变化，方程左边是精确解，方程右边是近似解，两个一减就是残差。有了残差以后，再选择加权函数与其内积，另加权函数内积为0（也就是余数R=0），那么就可以让近似解