矩量法作业.docVIP

摘要 矩量法是将连续方程离散为代数方程组的方法，此法对于求解微分方程和积分方程均适用。本文以直立的线天线为例，详细介绍了矩量法的基本思想及原理、基函数以及检验函数的选择，系统阐述了直立的天线的双位积分方程的建立过程，最后利用矩量法求解双位积分方程、使用matlab编程实现得到直立的天线上的电流分布。本文采用1.77λ的线天线的频率是0.3GHz，采用矩量法所选用的基函数是正弦函数，检验函数是函数。 关键字：线天线；矩量法；双位积分方程 引言 计算电磁学被广泛地定义为一门内在和常规地应用数字计算机来获得问题的数值结果的科学。主要有类别数值方法和高频或渐进方法，通常数值方法用于天线或散射体的尺寸一个到几十个波长量级的场合矩量法数学本质是一种求解的方法(有限元法)、 FDTD（时域有限差分法）线基本思想矩量法（Method of Moments, MoM）是一种将离散化为代数方程组的方法，对求解和均适用基于频域方程的， 根据线性空间的理论，N个线性方程的联立方程组、微分方程、差分方程、积分方程都属于希尔伯特空间中的算子方程，这类算子可化为矩阵方程求解。 设有算子方程： 　　　　　　　　　 　　　　　　　　　　　　　　　　　　（2-1）　 式中L为算子，可以是微分方程、差分方程或积分方程；g是已知函数如激励源；f为未知函数如电流。 假定上述方程的解存在且是唯一的，则有逆算子 的存在，使 成立。 互为逆算子 令：　　　　　　　　 　　　　　　　　　　　　　　　　（2-2） 式中是系数。被称为展开函数或基函数。用算子的线性便可以得到: 　　　　　　　　　　　　　　（2-3） 规定了一个适当的内积，，那么在的值域内定义一个权函数或检验函数…的集合，并对每个取式(2-3)的内积，则 　　　　　　　　　（2-4） 式中＝1，2，3…。此方程组可以写成如下的矩阵形式 　　　　　　　　 　　　　　　　　　　　　　　　　（2-5） 式中 　　　　　　（2-6） 　　　　　　　　　 （2-7） 如果矩阵是非奇异性的，于是可写成： 　　　　　　　　　　　　　　　　　 （2-8） 矩量法的求解步骤 依以上原理所述，矩量法的求解步骤可以分为以下四步： （1）将未知量展成由基函数构成的级数 （2）选取与基函数内积的检验函数 （3）由内积构成矩阵方程 （4）解矩阵方程求得未知量 基函数和检验函数的选择 矩量法求解算子方程的关键问题是基函数和权函数的选择。当选择基函数等于检验函数时,称为伽略金法。在特定的问题中,主要任务是选择基函数和检验函数,它们必须是线性无关的。基函数可以分为整域基和分域基;权函数一般有点匹配法,伽略金法,最小二乘法等。应用矩量法需注意:①误差分析;②方程收敛性;③积分奇异点处理等。下面分别介绍基函数和检验函数的选择： 本文的基函数选择的是正弦函数,其中k为波数。选取这样的基函数是考虑到它可以满足细导体末端电流为零的边界条件。文中用点匹配法选择检验函数。下面将简要介绍点匹配法。 若选取狄拉克函数为检验函数，即令 　　　　　　　　　　　　　　（2-9） 式中，狄拉克函数的定义为 　　　　　　　　　　　　　　　　　　　　　　　（2-10） 由于，因而矩量法方程(2-6)( 2-7)相应矩阵元素的计算结果为 　　　　　　　　　　　　　（2-11） 　　　　　　　　　　　　　　 　（2-12） 由此表明和的计算归结为只需计算所在点处的对应值，因此称这种方法为点匹配法。 3 对称振子天线的矩量法求解 假设对称振子天线放在Z轴线上，原点位于中点。如下图所示： 图3.1 对称振子天线示意图 图3.1为的对称振子天线，可做如下假设：（1）电流沿导线轴流动，体电流密度可以用线电流I来近似，体电流密度用线电荷密度近似；（2）忽略天线端面的周向电流径向电流；（3）天线上电流仅为长度变量的函数，即 3.1 由海伦方程应用矩量法的推导分析过程 当对称振子天线上在加上馈电信号时，或在接收电场作用时，其上将产生电流，其电流分布按基函数展开，在此，我们将选用正弦函数作为基函数。 通过以上假设可以知道，半波振子天线的矢量磁位，而 （3-1） 其中为格林函数，通过麦克劳林级数展开，舍弃高次项，这样其稳定度也较高。的表达式为： （3-2） 为波数。 于是z分量的电