矩量法作业DOC

摘要矩量法是将连续方程离散为代数方程组的方法，此法对于求解微分方程和积分方程均适用。本文以直立的线天线为例，详细介绍了矩量法的基本思想及原理、基函数以及检验函数的选择，系统阐述了直立的天线的双位积分方程的建立过程，最后利用矩量法求解双位积分方程、使用matlab编程实现得到直立的天线上的电流分布。本文采用1.77λ的线天线的频率是0.3GHz，采用矩量法所选用的基函数是正弦函数，检验函数是函数。关键字：线天线；矩量法；双位积分方程引言计算电磁学被广泛地定义为一门内在和常规地应用数字计算机来获得电磁问题的数值结果的科学。主要有两个类别：数值方法和高频或渐进方法，通常数值方法用于天线或散射体的尺寸在一个波长到几十个波长量级的场合；矩量法数学本质是一种求解线性方程的方法，在电磁场数值计算方面得到了广泛的应用。上世纪60年代，哈林登（Harrington）首先将矩量法应用于电磁场领域。对于辐射问题，是将含有未知电流的积分方程化为矩阵方程，通过求解该矩阵方程得到天线上的电流分布。在现代电磁工程中对于边界不复杂的问题可用解析法得到精确解,较复杂的边值问题,用解析法不能得到解答,需用数值法,如:MOM、FEM(有限元法)、FDTD（时域有限差分法）等。矩量法是一种误差最小的方法且相对其他的算法来说是一种相对简单高效的方法。对于天线问题可以建立描述天线表面感应电流的积分方程,求出该电流即可进一步得到天线的辐射特性。在矩量法分析过程中，有多种不同的积分方程可供选择，如双位积分方程、Hallen积分方程、Pocklington积分方程、Schelkunoff积分方程、响应积分方程等，这些积分方程一般应用于细的线天线，即假设电流只沿天线轴线流动，忽略天线周向电流和端面径向电流。本文采用双位积分方程，应用矩量法直立的线天线进行深入的分析，并在计算过程中对积分奇异项进行处理，不但简化了数值分析的复杂性，而且大大地减小了计算量并提高其精度。矩量法的基本原理矩量法的基本思想矩量法（MethodofMoments,MoM）是一种将连续方程离散化为代数方程组的方法，对求解微分方程和积分方程均适用。矩量法是基于频域方程的，它的基本思想是将一个泛函方程化为矩阵方程并求解该方程。由于求解过程中需要计算广义矩量，故得名。矩量法包括如下三个基本过程：（1）离散化过程主要目的是将算子方程化为代数方程，具体步骤是：①在算子L定义域内适当的选择一组线性无关的基函数fn；②将待求函数f表示为该组基函数的线性组合；③利用算子的线性，将算子方程化为代数方程。（2）取样检测过程主要目的是将求解代数方程的问题转化为求解矩阵方程的问题。基本步骤是：①在算子L的值域内适当的选择一组线性无关的权函数Wm；②将Wm与代数方程取内积进行N次抽样检验；③利用算子的线性和内积的性质，将N次抽样检验的内积方程化为矩阵方程。（3）矩阵求逆过程。R.F.Harrington在《计算电磁场的矩量法》一书中对其原理及过程进行了详尽的介绍.它所做的工作是将积分方程化为差分方程,或将积分方程中积分化为有限求和,从而建立代数方程组,故它的主要工作量是用计算机求解代数方程组.所以,在矩量法求解代数方程组过程中,矩阵规模的大小涉及到占用内存的多少,在很大程度上影响了计算的速度.如何尽可能的减少矩阵存储量,成为加速矩量法计算的关键.矩量法原理根据线性空间的理论，N个线性方程的联立方程组、微分方程、差分方程、积分方程都属于希尔伯特空间中的算子方程，这类算子可化为矩阵方程求解。设有算子方程：()Lfg（2-1）式中L为算子，可以是微分方程、差分方程或积分方程；g是已知函数如激励源；f为未知函数如电流。假定上述方程的解存在且是唯一的，则有逆算子的存在，使)(1gLf成立。互为逆算子令：nnnfaf（2-2）式中na是系数。nf被称为展开函数或基函数。用L算子的线性便可以得到:()nnnaLfg（2-3）规定了一个适当的内积f，g，那么在L的值域内定义一个权函数或检验函数123,,…的集合，并对每个m取式(2-3)的内积，则,,nmnmnaLfg（2-4）式中m＝1，2，3…。此方程组可以写成如下的矩阵形式mnnmlag（2-5）式中111212122212,,,,,,,,,nmnnmmmnLfLfLflLfLfLfLfLfLf（2-6）12nnaaaa12,,,mmgggg（2-7）如果矩阵l是非奇异性的，于是f可写成：1TTnnnmnmffaflg（2-8）矩量法的求解步骤依以上原理所述，矩量法的求解步骤可以分为以下四步：（1）将未知量f展成由基函数nf构成的级数（2）选取与基函数内积的检验函数mw（3）由内积构成矩阵方程mnnnlag（4）解矩阵方程求得未知量na基函数和检验函数的选择矩量法求解算子方程的关键问题是基函数和权函数的选择。当选择基函数等于检验函数时,称为伽略金法。在特定的问题中,主要任务是选择基函数和检验函数,它们必须是线性无关的。基函数可以分为整域基和分域基;权函数一般有点匹配法,伽略金法,最小二乘法等。应用矩量法需注意:①误差分析;②方程收敛性;③积分奇异点处理等。下面分别介绍基函数和检验函数的选择：本文的基函数选择的是正弦函数=1(z')=sin[(-)]NnnInkLz,其中k为波数。选取这样的基函数是考虑到它可以满足细导体末端电流为零的边界条件。文中用点匹配法选择检验函数。下面将简要介绍点匹配法。若选取狄拉克函数为检验函数，即令()mmrr（2-9）式中，狄拉克函数的定义为0()()()mmmrrrrrr（2-10）由于()()()vfrrrdvfr，因而矩量法方程(2-6)(2-7)相应矩阵元素的计算结果为,()()mnmnmnnmvlLfrrLfdvLfr（2-11）,()()mmmmvggrrgdvgr（2-12）由此表明mnl和mg的计算归结为只需计算mr所在点处的对应值，因此称这种方法为点匹配法。3对称振子天线的矩量法求解假设对称振子天线放在Z轴线上，原点位于中点。如下图所示：图3.1对称振子天线示意图图3.1为2L=/2的对称振子天线，可做如下假设：（1）电流沿导线轴流动，体电流密度J可以用线电流I来近似，体电流密度用线电荷密度近似；（2）忽略天线端面的周向电流II和径向电流；（3）天线上电流仅为长度变量z'的函数，即I=I(z')，与，无关。3.1由海伦方程应用矩量法的推导分析过程当对称振子天线上在加上馈电信号0V时，或在接收电场iE作用时，其上将产生电流，其电流分布按基函数展开，在此，我们将选用正弦函数作为基函数。通过以上假设可以知道，半波振子天线的矢量磁位=zzzAeAAe只有分量，所以，而-LA(,z)=(z')G(zz')dz'4LZI，（3-1）其中G(zz')，为格林函数，通过麦克劳林级数展开，舍弃高次项，这样其稳定度也较高。G(zz')，的表达式为：2222exp[-jk+(z-z')]G(zz')=4+(z-z')aa，（3-2）2k=为波数。于是z分量的电场表达式为：222=+zzzAjEkAz（3-3）由边界条件可知：=-zizeEE（由边界条件可知，在导体表面切向电场为零），其中iE（z）为电源所产生的场，令0z=(z)zieEV（），0V为电压，(z)为狄拉克函数，将上述条件和（3-1），（3-2）式代入（3-3）于是便得到：2202-L1-(z)=(+k)(z')G(zz')dz'4LjVIz，（3-4）对上式进行求解为：0-L01(z')G(zz')dz'=Bcoskz-jsin42LVIkz，（3-5）公式中B为常数，0=120为自由空间的波阻抗。该方程为对称振子天线的海伦积分方程（Hallen）。海伦积分方程(3-5)简单变型为：0-L(z')G(zz')dz'=4cos-sin60LVIBkzjkz，（3-6）选择分布电流展开表达式：=1(z')=sin[nk(-)]NnnILz（3-7）=1I(z')Nnnnfn为待定系数（3-8）本文中取N=4，那么电流表达式变为：4=1(z')=sin[nk(-)]nnILz，代入海伦积分方程式（3-6），整理得到：401sin',''cossin60LnLnVankLzGzzdzCkzkzj（3-9）式中：4CB，rezzGjkr)',(，2122])'[(azzr式中L为天线的一般长度，a为导体半径，z为导体表面上场点的坐标，'z为导体表面上源点的坐标，在z轴线上。（3-9）式中有五个未知量1、2、3a、4a、C，因此应选5个加权函数作为五个方程式。课题中我们求解的是L=λ/4，λ为波长，为了求对称振子天线上电流分布I(z')，采用点配法，检验函数wm选择狄拉克函数，即mw=(z-z)m，其中m=1,2……N，其中选择的点为：(m-1)z=mLN变换，即：10z，214zL，324zL，434zL，5zL。运用选择的mz的值对每一mW对（3-9）式两边求内积，这样可以将（3-9）式转变成矩阵的形式:1111111222222203333333444444455555560ABCDEPABCDEPVABCDEaPjABCDEPaABCDEPC(3-10)式中2222,sin',''exp'sin''LmmLmLLmAzzkLzGzzdzjkzzakLzzza2222,sin2',''exp'sin2''LmmLmLLmBzzkLzGzzdzjkzzakLzzza2222,sin3',''exp'sin3''LmmLmLLmCzzkLzGzzdzjkzzakLzzza2222,sin4',''exp'sin4''LmmLmLLmDzzkLzGzzdzjkzzakLzzza,coscosmmmEzzkzkz00,sinsin6060mmmVVPzzkzkzjj其中除了1，2，3，4，C五个未知量，其余各个元素都可以通过计算机编程来求解，得出各矩阵形式，然后再利用矩阵求逆运算，求解出1，2，3，4，C，这样根据公式（3-7）就可以求得天线上的电流分布。本文是基于海伦积分方程的半波振子天线矩量法分析，选择的条件为：=2GfHz，波长为=0.15m，振子长度=0.250.0375Lm，310022.7a/λ，01VV，基函数采用三角正弦函数，检验函数采用函数，采用matlab计算出的电流系数矩阵如下所示：3.2Matlab编程程序流程图Matlab软件在矩阵处理方面有其独特的优势，故采用matlab语言编写程序仿真天线的性能。Matlab实现矩量法程序流程图如下图所示。3.3由海伦方程应用矩量法计算得到的结果图及其分析利用公式（3-7）可得电流分布表达式为：4=1(z')=sin[nk(-)]nnILz从而得到实虚部电流分布图如图3.2所示：开始输入天线初始参数选定基函数和检验函数矩阵元素Z11，Z12，......Zmn的计算计算电流分布（