《计算电磁学》演示文稿.pptVIP

计 算 电 磁 学 概 论 二. 矩量法 1. 基本原理 对天线问题，积分方程的基本形式为 方程的核（已知） 激励函数（已知） 对激励的响应（未知） 写成算子方程 激励函数 待求函数 算子 二. 矩量法 在 L 的定义域内将待求函数 f 展开为一组线性无关的已知函数 fn 的集合 基函数 待定系数 （1）基函数展开 当 时， ; 当n为有限项时，f 为近似值。 代入到算子方程中，得 二. 矩量法 （2）加权求内积 在 L 的定义域内选取另一组已知函数系列 wm(m=1,2,…,N)，将每个wm对上式取内积，得到含N个未知量的线性方程组 其中，(m, n=1,2,…,N) 权函数 方程两边与权函数求内积的目的：使上式两边之差——余量的加权平均值最小，以确定未知系数 。 二. 矩量法 写成矩阵形式 其中， 二. 矩量法 （3）求解矩阵方程 若矩阵[l]是非奇异的，即[l]－1存在： 将求得的 代入到 ，便得到 f 的近似解。 二. 矩量法 2. 基函数和权函数的选择 基函数和权函数的选择是否恰当，关系到 解的精确度 矩阵元素计算的难易 能够求逆的矩阵大小 良态矩阵的可实现性 基函数和权函数的组合也有多种方式。 二. 矩量法 （1） 基函数 全域基函数：在 L的定义域内都有定义的基函数，应满足边界条件且彼此线性无关。 傅立叶级数： 幂级数： 切比雪夫多项式 勒让德多项式 二. 矩量法 分域基函数：在定义域中的各分域上才有定义的基函数。 脉冲基函数（分段单位函数）： 三角基函数（分段线性函数）： 0 1 0 1 x x 二. 矩量法 （2） 权函数 加略金法：选取权函数等于基函数，即 点匹配法：选择 函数为权函数，使未知函数的近似解只在某些离散点上满足算子方程式，即 利用 函数的挑选性: ，有 表示只需计算基函数 fn和激励函数 g 在 点处的值。取消了加权积分，简化了运算。 二. 矩量法 小结： 优点：只需很少几项展开函数的线性组合可较好地逼近未知函数；收敛很快 全域基 缺点：选择合适的全域基函数很困难 优点：简单灵活，使用方便，应用广泛 分域基 缺点：收敛较慢，计算量大 优点：实施容易，计算量小 加略金法 缺点：计算精度较低 点匹配法 优点：计算精度高 缺点：实施较繁，计算量大 二. 矩量法 例1 在区间[0,1]上，位场满足泊松方程和边界条件为 求解 解1 选择全域基的幂级数作为基函数，它满足边界条件，且在[0,1]全域存在 其精确解可由直接积分法求得 二. 矩量法 则近似解表示为 代入到算子方程，得 采用加略金法，选择权函数 ，与上式作内积 写成矩量方程为 二. 矩量法 其中矩阵元素为 二. 矩量法 当N＝1时， 则 当N＝2时，矩阵方程为 解得 则 二. 矩量法 当N＝3时，矩阵方程为 解得 则 近似解即为精确解。 二. 矩量法 矩量法解 二. 矩量法 解2 仍然选择幂级数为基函数，有 选择点匹配法， 函数为权函数 将区间[0,1]等分为N＋1个区间，匹配点的位置为 二. 矩量法 代入到矩量方程中作内积，各矩阵中元素 二. 矩量法 当N＝1时， 则 当N＝2时，矩阵方程为 解得 则 二. 矩量法 当N＝3时，矩阵方程为 解得 则 近似解即为精确解。 二. 矩量法 例2 Pocklington方程的矩量解 解1 脉冲基函数和点匹配法 z=-L/2 z=L/2 z=0 z 偶极子线天线 k Pocklington方程为： 将天线等分为N段，h=L/N。用脉冲函数 Pn(z’)展开天线上的电流 其中， 二. 矩量法 将 函数定义在天线表面的匹配点上，匹配点选在每分段中点zm Pocklington方程变为： 写成矩阵形式： 广义阻抗矩阵 广义电流矩阵 广义电压矩阵 二. 矩量法 解2 全域基函数和点匹配法 用全域基余弦函数展开天线上的电流 代入到Pocklington方程，有 选择 函数为权函数，定义在天线表面的匹配点zm上，则 * 一. 引 言 二. 矩量法 三. 有限元法 四. 时域有限差分法 邵 维 参 考 文 献 R.F.哈林登,计算电磁场的矩量法,国防工业出版社,1981。 李世智,电磁散射与辐射问题的矩量法,电子工业出版社,1985。 金建铭,电磁场有限元方法,西安电子科技大学出版社,1998。 王秉中,计算电磁学,科学出版社,2002。 葛德彪,闫玉波,电磁波时域有限差分法(第二版),西安电子科技大学出版社,2005。 盛新庆,计算电磁