第10章 矩量法讲解

第十章 矩量法

解析方法仅适用于结构简单的散射体。如果散射目标结构复杂，必须选用数值方法。数值方法是对所求解的微分方程或积分方程实施离散，采用一组基函数表示电场、磁场或感应电流等未知量，然后将电磁场微分方程或积分方程转换为一组线性代数方程，即可按照标准的数值程序求解这些线性方程组。数值方法的优点在于容易处理结构复杂的散射体，而且通常可以获得高精度解。随着高性能计算机的飞速发展，数值方法已经成为解决实际问题的日益重要的工具。现今已有多种数值方法，各具特色，分别适用于求解不同的电磁问题。典型的数值方法是矩量法（MoM ）、时域有限差分法（FDTD ）和有限元法（FEM ）等。本章讨论矩量法，后两章将分别介绍时域有限差分法和有限元法。

矩量法是求解算子方程的有效方法，这些算子通常是微分算子、积分算子或者是两者的组合。20世纪60年代， R. F. Harrington 首先将矩量法用于电磁问题的求解[1]。目前已经广泛地用于天线分析、微波器件的设计以及复杂目标的雷达散射截面（RCS ）的计算。通常认为矩量法是精度最高的数值方法，因此引起更多的关注。如今很多商用软件的开发都基于矩量法。但是，矩量法需要求解稠密的矩阵方程。对于电大尺寸的散射体，它将十分消耗大量机时及内存。为了解决这个问题，人们作了很多努力，研发快速计算和有效的存储方法。因此发展了很多有关积分方程的快速求解算法，大力推动了矩量法的应用。 10-1一般步骤

典型的算子方程可以表示为下列形式

h Lf = （10-1-1）

式中L 为线性算子，可以是微分、积分或两者组合，h 为一个已知函数，f 为待求的未知函数。这些函数可以是矢量或标量，且定义域可为一维、二维或三维空间。因此，在电磁学中它们可以是空间及时间函数。矩量法的一般步骤是，首先将未知函数表示为一组基函数的线性组合，然后匹配算子方程，最后由离散的线性方程组求出展开系数。下面详述矩量法的具体步骤。

首先令N f f f ,,,21 为一组基函数，那么，未知函数)(x f 可以近似表示为

∑==+++≈N

n n n N N x f a x f a x f a x f a x f 1

2211)()()()()(

（10-1-2）

式中),,3,2,1(N n a n =为展开系数，它们是未知的。如果N 足够大，上述表示式将非常精确。将上式代入式（10-1-1），得

)()(1

x h x Lf a N

n n

n

≈∑=

（10-1-3）

下一步是选择一组权函数，N w w w ,,,21 ，以每个权函数与上式各项逐一相乘，并且在未知函数的定义域内求积，建立一组未知系数为n a 的线性代数方程。该组方程可以表示为

1

, 1,2,3,,N

mn n

m n Z

a b m N ===∑

（10-1-4）

该方程组的系数及右边项分别为

⎰=x Lf x w Z n m mn d )( （10-1-5） ⎰=x x h x w b m m d )()(

（10-1-6）

求出未知系数后，即可近似地决定未知函数，并由此求得其它场量。

上面简述了矩量法的求解过程，现在需要讨论几个问题。

首先是基函数的选择。对于基函数的两个基本要求是完备性和正交性。完备性是指选择的基函数可以精确地表示任何未知函数，且其精度随着基函数的数目增加而提高。正交性可以放宽为线性独立，即要求一组基函数中任何两个必须是线性独立的。众所周知，一组线性独立函数总可以应用所谓Gram-Smit 方法使其正交化。此外，表示式的有效性通常也是选择基函数的重要判椐。

如果基函数可使未知函数易于满足实际的边界条件，那么即是一种较好的基函数。具有实际应用的典型基函数有两种：其一称为全域基函数，另一个称为分域基函数或称为子域基函数。每一个全域基函数都在相同的域中定义，而每一个分域基函数的非零区域是在未知函数的部分域中定义。例如下列积分方程

2d )(),(1

1

='''⎰

-x x f x x g

（10-1-7）

其中未知函数)(x f 的定义域为]1,1[-，因此，基函数]1,1[ ),2/sin()(-∈=x nx x f n 是一组全域基函数。全域基函数的很大优点是各个基函数具有相同的表示精度。与分域基函数比较，采用全域基函数时通常待求的未知数的数目较少。因为使用全域基函数时无须网格剖分，数值计算也相对地易于实现。一些有用的全域基函数是多项式（例如 ,,,,,12n x x x 。），以及正弦和余弦函数。例如，对于区域],0[a ，可以选择 ),/3sin(),/2sin(),/sin(a x a x a x 作为一组基函数。全域基函数通常用于求解一维问题的线性积分方程，定义域为矩形的二维问题有时也可采用。全域基函数也可与分域基函数组合使用。一个典型的例子是，旋转体散射问题的求解。此时，每个基函数是一个随角度变化的全域基函数与一个轴向变量的分域基函数（例如方波函数或三角函数）的乘积。全域基函数的主要缺点是，它们仅可用于形状规则的定义域，例如一维导线和二维的方形或矩形域。对于边界形状复杂的区域，定义全域基函数是十分困难的。

一个分域基函数仅在部分函数域内定义有非零值，通常这部分域的尺寸远小于波长。除了方波函数（又称为脉冲基函数）以外，几乎全部分域基函数（例如三角函数）具有重叠的非零区。因此，为了定义分域基函数，通常需要将求解区域划分为很多小片的集合，每个小片称为一个网格单元或简称为一个单元。这样的单元集合构成目标的近似表示，因此称为网格。一些典型的网格单元形状如图10-1-1所示。典型的网格单元的形状，对于线状结构为线段；对于面状结构为三角形和方形；对于体状结构为立方体、四面体及棱柱体[2~5]。

实际上，广泛地选择矩形基函数和三角基函数作为分域基函数[2~5]。许多其它基函数是该两种基函数的变形或组合。一个矩形基函数在一个单元内定义为1，而在其余全部单元内定义为零。因此，任何两个矩形基函数的非零的子域不会重叠。三角基函数恰好相反，每个三角基函数在两个单元中具有非零区域。显然，基函数的选择有时与网格的形状有关。以一维为例，将未知函数定义在一个区域[0,1]，再将该区域划分为相同尺寸的N 个子区。第n 个子区的数学定义为

ΔΔ)1(n x n ≤≤-，这里N n ,,3,2,1 =，而N 1/Δ=。

一个矩形基函数)(x f n 定义为

⎩⎨

⎧≤≤-=其它区域

,0ΔΔ)1(

,1)(n x n x f n （10-1-8）

该函数如图10-1-2所示。

(a)

(b) (c) 图10-1-1 典型的网格单元：(a)线段，(b)平面三角形，(c)六面体