电磁数值方法的数学基础.docVIP

第二章 矩量法的基本原理 §2.1. 矩量法 矩量法（method of moment）大约在20世纪60年代开始引进电磁学，从此矩量法在各种各样电磁学研究中得到了广泛的应用。包括天线问题、天线设计、微波网络、散射及辐射效应、微带结构分析、电磁兼容等各种问题。在矩量法的发展中，许多电磁学科学家做出了贡献，如Richmond，Mittra，Wilton等，特别是Harrington。（MOM之父？） Roger F. Harrington K.K. Mei（梅冠香） 矩量法在电磁场分析中有着广泛的应用，其基本思想是：将泛函方程化为矩阵方程，并求解该矩阵方程；利用线形空间和算子来加以表达。 §2.1.1 如果非齐次方程为 （2-1） 式中是线性算子，为已知函数，为未知函数。令在的定义域中被展开为的组合，如 （2-2） 式中是系数。被称为展开函数或基函数。对于精确解，式（2-2）通常是无穷项之和，而形成一个基函数的完备集。对于近似解，则通常是有限项之和。式（2-2）代入式（2-1），再应用算子的线性便可以得到 （2-3） 对此问题若定义一个适当的内积，在的值域内定义一个权函数或检验函数的集合，并对每个取式（2-3）的内积，则 （2-4） 式中。此方程组可以写成如下的矩阵形式 （2-5） 式中 （2-6） （2-7） 如果矩阵是非奇异性的，其逆矩阵存在，则便由下式给出 （2-8） 的选取，不要求与相同，具有灵活性。当选择这种特殊情况时，通常称为伽略金法。如果自伴，则，成为对称矩阵。 在任何一个特定的问题中，主要任务是选择和。必须是线性无关的，并且使得它们的某种叠加形式能够很好地逼近。也应该是线性无关的，并且也应该使得内积取决于的相对独立性。 影响选择的其它因素： ① 精度取决于的选取； ② 计算矩阵单元的难易； ③ 能够求逆的矩阵大小； ④ 良态矩阵的可实现性。 §2.1.2 点选配 在实际问题中，式（2-6）的中的积分计算通常是很困难的。求近似解的一个简单方法是在所关心的区域内，要求在一些离散点上满足方程式（2-3），这种方法称之为点选配。就矩量法而言，这相当于用狄拉克函数作为检验函数，使得内积的计算简单化。 取，则 （2-9） 局限：选取的点数比较大。 分域基 另一种实用的近似方法是分域基函法。它采用的各个基函数只是在定义域的各个分域上存在，于是，展开式的每个只是在所关心区域的各个分域上才影响的近似。这种方法通常可以简化计算以及简化形成的矩阵。 脉冲基（积分方程较适合，不能用于微分算子） 图2-1 脉冲基函数 ： （2-10） 经典情况：整域基—分域基 三角形基（可扩展为屋脊基函数） （2-11） 图2-2 三角基函数 比脉冲基有所改进，可求一次微分、二次微分。 （2-12） （2-13） 分段三角函数基 （2-14） 与全域基函数相比，分段基函数适用范围更广泛。三角函数基在段内具有较好的连续性，得到了广泛的应用。 矢量三角（RWG）基函数 矢量三角形基函数的定义如下式： （2-15） 式中各参量的意义见下图: 图2-3 矢量三角形基函数 基函数由正、负两个有一条公共边的三角形组成。当源点在正三角形中时，矢径从正三角形顶点指向源点；当源点在负三角形中时，矢径从源点指向负三角形的顶点。源点处的电流大小与矢径成正比，且方向与矢径的方向一致。 矢量三角形基函数具有如下特点： 1. 在基函数的外边界上，即两个三角形的四个外边上，电流没有法向分量，因而在边界上没有线电荷的存在。 2. 在两个三角形的公共边界上，电流的法向分量连续且为常数1。这一点可以从图2-3中看出，的法向分量的大小等于相应三角形的高，乘以系数，正好为1。这也说明，不仅两三角形的外边上没有线电荷，而且在交界的边上也没有线电荷。同时，其法向分量为1，可以使我们认识到基函数对应的电流系数的意义就是通过公共边的法向电流大小。 3. 矢量三角形基函数的表面电荷密度，正比于表面电流的散度。它在两个三角形内为常数，而且整个三角形上的面电荷总量为0。 图2-4 矢量三角形基函数公共边上的电流 （2-16） 这样，对于给定的表面，首先对其进行三角形划分。然后可以用这种矢量三角形基函数进行模拟表面电流: （2-17） 这里，N是矢量三角形基函数的数目。是通过该边的法向电流的大小。 §2.2. 矩量法与变分法 伽略金法与瑞利-里兹变分法等效，而一般的矩量法也是一种变分方法。利用线性空间的概念来解释矩量法。设表示的值域，表示由张成的空间，表示由张成的空间，于是矩量法意味着精确值在上的投影等于近似值在上的投影。由于获得投影的方法使误差化为最小，所以矩量法是一种使误差为最小的方法。而由变分法可以得到与此相同的