37 矩 量 法.pptVIP

3.7 矩 量 法 一、矩量法的基本原理 首先以有 域内满足第一类边值条件的本征值问题为例说明矩量法的基本原理。本征值问题写成一般的形式为 所以矩量法利用基函数和权函数将最初的本征值问题（式（3.7-1））转换成了矩阵的本征值问题（式（3.7-6））,通过求解矩阵方程可到近似解。 选定基函数和权函分别为 其中 是微带片上的总电荷。设微带片的电荷密度为 ，微带片上的电压与电荷密度间满足积分方程： * * * * * * * * *第三章 静电场边值问题解法 * 矩量法是由哈林登（Harrington）1968年提出的，已成功地用于求解许多实际的电磁问题。本节介绍矩量法的基本原理，包括矩量法的解题过程，基函数和权函数的选择。 其中算子 可以是微分算子也可以是积分算子。在前面中，我们已经讨论了在横截面为 的柱状空心波导中电磁波的传播规律，即 （3.7-1） 第一步是将式（3.7-1）中的未知函数 近似表示成函数 的线性组合，即 其中 是沿柱状导轴向的传播常数，这是一典型的本征值问题,本征值 算子 ，本征函数 。 在波导壁上 矩量法的解体步骤: （3.7-2） 第二步是用权函数 （又称检验函数）对式（3.7-3）两边取内积，即有 其中函数 是已知的独立函数，称为基数; 是未知的待定系数，将式（3.7-2）代入式（3.7-1）得 （3.7-3） 式（3.7-4）可以重新写成 （3.7-4） 矩阵 是 阶矩阵， 是 阶矩阵， 是 阶矩阵。 （3.7-5） 将式（3.7-5）写成矩阵的形式 （3.7-6） 其中: 解方程（3.7-7）可求得 个本征值 ,对每一个本征值 ,由式（3.7-6）可求得本征矢量 为使矩阵方程（3.7-6） 有非零解,其系数矩阵 的行列式必须为零，即 （3.7-7） 最后求得相应的本征函数 例1．求不定解域为的本征值问题。 解 将本征函数近似表示成 将选定的基函数和权函数代入式（3.7-6） 其中: 为简单起见，选 ，则方程（3.7-6）变成 （3.7-7） 为使上述方程有非零解，其系数矩阵的行列式必须为零，即本征值心须满足方程 （3.7-8） 求得两个本征值分别为 。 其中 可任意选取，选择它满足 从而对应的 本体征函数为: 求得: 同理可求得: 因此本征值问题的近似解为: 其精确解是 图1中给出了 的精确解（图中实线）和用矩量法求得的近似解（图中虚线）的曲线，图中可以看出两者是非常接近的。 图1 本征值问题的矩量法解与精确解的比较 二.基函数和权函数的选择 在例1中选用的是整个定解域上的基函数，称为全域基函数。除了全域基函数外也可采用分域基函数，它只定义在函数定义域的子域上。常用的分域基函数有脉冲函数与三角形函数。 1.脉冲函数 脉冲函数定义为: 对于一维问题，如图2所示，假定函数的定义域为 ，将定义域分成 个宽度相同的子区间，每个子区间的宽度为 ，其中 系数 等于函数 在 处的值，但是对包含二阶导数 的算子不能选脉冲函数作为基函数，这是因为脉冲函数的二阶导数包含有对 函数的导数 则脉冲基函数为 待求的函数用脉函数的线组合近似，于是有 其中: 2.三角形函数 如图4所示，三角形基函数定义为: 函数 可近似地用三角形基函数的线性组合来表示，即有 前已述及，矩量法的第二步是用权函数 求内积是很困难的，但如果将权数选为 函数，即 其中: 3.权函数的选择 例2求表示在图5中的微带片状电容器的电容。 解 设地为电位参考点，加在微带片上的电压为 ，根据电容的定义，微带片的电容为: 其中： 它表示 上的电荷密度是均匀的，数值为 ， 采用点配法，权函数为： 采用矩量法求解积分方程,将微带片的宽 和 长 分别等分成 等分,即将微带片分成 个小矩形的宽为 ,长为 ，每个小矩形的面积为 。选用脉冲基函数，将板上的电荷密度 表示成 * * * * * * * * *第三章 静电场边值问题解法 * * *