《矩量法 of MomentPPT》由会员分享,可在线阅读,更多相关《矩量法 of MomentPPT(62页珍藏版)》请在金锄头文库上搜索。
1、1,第二章 矩量法(Method of Moment) 2.1 引言 2.2 矩量法的一般过程 2.3 选配和离散过程 2.3.1 点选配 2.3.2 脉冲分域基 2.3.3 三角形函数分域基 2.4 算子研究 2.4.1 近似算子 2.4.2 扩展算子 2.4.3 微扰算子,2,矩量法(简称MoM),就其数值分析而言就是广义Galerkin(伽略金)法。矩量法包括两个过程,离散化过程和选配过程,从而把线性算子方程转化为矩阵方程。这里先举一个简单的例子。,3,例1无限薄导体圆盘上的电荷分布问题。 试讨论半径为a的无限薄理想导体圆盘,在中心线距离d处有一点电荷 ,如图5-17-1所示,求解导体圆
2、盘上的电荷分布。 解 假设导体圆盘上电荷密度为 ,根据电磁学的基本概念可知: (1) 由外加电荷Q在导体圆盘上产生的电位e 和导体圆盘本身感应电荷密度所产生的电位i之和U 在盘上处处相等,即保证导体圆盘是等位面。 (2) 由于本问题中是感应电荷,因此总电荷Qi0,其中,4,图5-17-1导体圆盘上的电荷分布 (5-17-1) (5-17-2) (5-17-3),5,于是,问题可写为 (5-17-4) 式中r= ,其中打撇的表示源点,不打撇的表示场点。 这个问题,采用电磁学经典解析方法不能很好的解决,因为未知量 处于积分内部,是一个典型的积分方程。为此,把圆盘分割成两部分:中心小圆和外部环带(如
3、图5-17-1所示),并假定每一部分内的电荷密度 (i=1,2)近似为常数,于是 (5-17-5) 式中 (5-17-6),6,称为脉冲函数,这时问题方程(5-17-4)成为 (5-17-7) (5-17-8) 把问题方程(5-17-4)近似的转化为式(5-17-7)和式(5-17-8)的过程称为离散化过程。但是,必须注意到方程(5-17-7)中,场点r表示圆盘上的任意点(x,y),换句话它们是不定的,因而式(5-17-7)中包含着无限个方程。另一方面,离散后的方程组(5-17-7)和方程组(5-17-8)内只有三个未知数 、 和 ,于是方程组超定。,7,为了把超定方程组转化为唯一解的方程组,
4、可以采用很多办法。矩量法中,习惯用选配过程解决这个问题。简单说来,即在每个离散的单元上只选取一个场点作为代表来建立方程。例如,在例1中对于离散的 和 分别取 和 两点做试验点,如图5-17-2所示。具体写出方程组 (5-17-9) 其中,8,图5-17-2 圆盘上的试验点,9,其中 表示 面元电荷在 处产生场的自作用单元; 表示 面元电荷在 处产生场的自作用单元; 表示 面元电荷在 处产生场的互作用单元; 表示 面元电荷在 处产生场的互作用单元。,10,又有 (5-17-14) 经过离散化过程和选配过程,将积分方程组(近似地)转化为矩阵方程 (5-17-15) 由此得出电荷分布的解为 (5-1
5、7-16),11,图 5-17-3 矩量法的一般过程 图5-17-3所示的矩量法求解问题的一般过程。 讨论 (1)矩量法的原问题并不限于积分方程,也可以是微分方程或其他方程。但必须能抽象成算子方程。从这一点而言,它是普遍的;另一方面,矩量法最终要转化为矩阵方程加以解决。因此,原问题必须属于线性算子范畴。例如,最速下降线所构成的积分方程 不是线性泛函,所以无法采用矩量法。 (2) 电磁理论中计算的矩阵单元,一般均表示某个源在一个区域所产生的场,而实际产生的场往往都随着源的距离增加而减少。换句话说,矩量法中矩阵一般是对角占优的:自作用单元 比互作用单元 所起的作用要大。这一点在概念上十分重要。,1
6、2,矩量法的研究对象是一般非齐次方程 (5-17-17) 线性算子 的运算空间称为定义域,而 组成的空间称为值域。式(5-17-17)中 是已知的激励函数, 为未知函数。令 在 的定义域内展开成 的组合,有 (5-17-18),2.2 矩量法的一般过程,其中,13,表示矩阵转置,应该注意到:展开函数与基函数是有区别的。一般来说,基函数是一无限展开。从完备基转化为近似有限截断基已经构成误差了,再从有限截断基转化为有限展开函数就很难保证 能收敛于 ,这也是矩量法的研究中需要深入研究的一个问题。这里且写出 (5-17-19),而,从算子方程(5-17-17)到式(5-17-19)即构成离散化过程。它
7、可以是函数离散,也可以是区域离散,或两者兼有。,14,现在规定适当的内积 。在算子L的值域内定义一类权函数(或检验函数) ,作用于式(5-17-19)两边,且取内积,有 (5-17-20) 这就是所谓的选配过程或试验过程,矩量法的名称也由此而来,即把激励矢量 和 分别向权空间投影,取它的矩,根据矩的大小确定展开系数。 如果展开函数的数目与权函数数目相等,则可把式(5-17-20)写成矩阵形式 (5-17-21) 其中 (5-17-22) 于是可以解出 (5-17-23),15,若规定函数矩阵 (5-17-24) 于是待求的函数为 (5-17-25) 矩量法的一般过程的数学表示如图5-17-4所
8、示。 十分清楚,矩量法的结果优劣取决于:离散化程度; 和 的选取;线性方程组的求解。在 = 的特殊情况下,可称为Galerkin(伽略金)法,于是矩量法也称为广义Galerkin法。,16,图5-17-4 矩量法一般过程的数学表示,17,例2研究 ,其中 解 已经知道,此问题存在精确解 本例采用矩量法求解,选择 再选择权函数,18,即采用Galerkin法,内积定义为 于是可给出一般计算结果,19,归纳起来有,20,情况1:N=1,于是有,情况2: N=2,21,情况3: N=3,22,十分明显,N=3时已得到了精确解 。 矩量解的曲线如图5-17-5所示。,图5-17-5 u (x)矩量解,
9、23,第二章 矩量法(Method of Moment) 2.3 选配和离散过程 2.3.1 点选配 2.3.2 脉冲分域基 2.3.3 三角形函数分域基 2.4 算子研究 2.4.1 近似算子 2.4.2 扩展算子 2.4.3 微扰算子,24,2.3 选配和离散过程 从上面的典型例子可知,矩量法的精华在于选配和离散过程,值得单独进行研究。 2.3.1点选配 点选配是一种最简单而最典型的选配函数。因为矩阵单元为 ,一般说来,其中所含的积分计算十分困难,这种情况下,最简单的办法是做某些点的投影,即所谓的点选配,实际上相当于把权函数取为 函数。,25,例3 任研究 。 解 设 , 在这个例子中取
10、函数为权函数即 其中, 是这个问题的选配点,于是有,26,例3 任研究 。 解 设 ,可得到 在这个例子中取 函数为权函数即 其中, 是这个问题的选配点,于是有,27,归结起来,可写出,28,情况1: N=1,情况2: N=2,29,情况3: N=3,可以得出,对于点选配情况N=3,又一次回复到精确解。,30,讨论 (1) 对于点选配的情况,N+1阶矩阵中的N阶主子阵并不等于在N时的系数矩阵(和Galerkin情况不同)。因此当N逐渐变大时计算量无法节约。 (2) 点选配虽然看起来非常简单,然而其内在的道理极其深刻。这一点可以从数值积分看出。研究表明任何数值积分方法,不论矩形、梯形、二次样条等
11、,说到底都是选择积分区域的点和区域点所对应的系数,由此产生Gauss积分的思想。所以在矩量法中,研究最佳点选配将是一个十分有意义的课题 。,31,2.3.2脉冲分域基 矩量法在离散化过程中用展开函数取代基函数,带来了方便和自由。但是,随之而来的如何确保解的收敛性的问题却值得人们重视。 在尚未了解u(x)函数性态的条件下,采用有限个展开函数ui(x),i=1,2,.,N时要确保解收敛显然在理论上存在不少困难,采用分域基函数可以说是比较稳妥的一种解决方案。因为大多数良态函数(不做高速振荡)均可以采用有限段直线或样条加以逼近,如图5-17-6所示。,32,图5-17-6 分域基函数近似,33,下面从
12、最简单的脉冲函数着手展开讨论。 一般的脉冲函数可以表述为,(5-17-26),式(5-17-26)表示以i为中点,密度为1/(N+1)的脉冲函数,在实际情况下,密度可以根据问题灵活改变,如图5-17-7所示。,图5-17-7 脉冲函数 图5-17-8 三角形函数,34,2.3.3三角形函数分域基 三角形函数也是常用的一种分域基,如图5-17-8所示。 若采用三角形函数展开未知函数(x),则有 (5-17-27) 所得的解的合成相当于折线连接,分段三角形函数所得的折线包络如图5-17-9所示。 为了研究具体例子,这里先给出三角形函数的导数概念。引入如图5-17-10所示的阶梯函数H(x-xi),
13、其定义为,35,图5-17-9 分段三角形函数所得的折线包络,图5-17-10 H(x-xi) 函数,(5-17-28),36,再引入大家熟悉的Dirac-函数,也即脉冲函数,其定义为 (5-17-29),如图5-17-11所示。,图5-17-11 (x-xi)函数,37,Dirac-函数有两个重要的性质: 1.归一性 (5-17-30) 2.选择性 (5-17-31) 这里不加证明的给出Dirac-函数和阶梯函数之间的重要关系。 (5-17-32) 有了以上基础就可以把三角形函数的导数用阶梯函数H表示,具体为 (5-17-33),38,图5-17-12给出形象的几何表示。,图5-17-12
14、三角形函数导数的几何表示,例4 重新研究Harrington(哈林登)问题,L(u)=g,其中 L= ,g= ,边界条件为u(0)=u(1)=0。 试用以三角函数作为展开函数,脉冲函数作为权函数的矩量法求解。,39,解 根据要求可写出 于是有 上式已计及 选择权函数 于是矩阵单元 上式要分三种情况讨论。,40,此外,激励单元为 结果可归纳为,41,情况1:N=1 考虑到对比: 则有 和 的对比如图5-17-13所示。,图 5-17-13 和,42,情况2:N=2 l= g= 容易得到 同样对比有 和 的 对比图如图5-17-14所示。,图 5-17-14 和,43,情况3:N=3 于是有 同样
15、对比有 和 的 对比图如图5-17-15所示,44,图 5-17-15 和,讨论 分域基在N不大的情况下与精确解的差距是明显的。但是它的相应矩阵是三条带矩阵,可较明显地缩小计算量。因此选择N不大的分域基并进行顶点拟合将会是一个比较好的方案。,45,2.4 算子研究 算子方程是矩量法建模的关键。它应该有两个方面的要求: 一方面算子方程必须符合物理(或工程)问题的主要本质; 另一方面它又必须适合数值计算。 这两个方面构成了算子研究的基础。,46,2.4.1 近似算子 细心的读者一定会提出这样一个问题,即例4中为什么不采用脉冲函数作为分域基展开?其实原因十分简单,因为脉冲函数的二阶导数表示有很大困难。但是,倘若引进近似算子的概念,则可以较好地解决这个问题。 算子近似含义相当广泛。作为例子,可采用有限差分代替微分。,47,例5 研究 的Harrington问题,即 , ,试采用差分近似算子 ,脉冲展开点选配的矩量法求解。(做一般了解) 解 为确保 的边界条件