矩量法 Method of Moment课件

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 微扰算子微扰算子 矩量法矩量法(简称简称MoM)，就其数值分析而言就是广义，就其数值分析而言就是广义Galerkin(伽略金伽略金)法。矩量法包括两个过程，法。矩量法包括两个过程，离散化过程离散化过程和选配和选配过程过程，从而把线性算子方程转化为矩阵方程

2、。这，从而把线性算子方程转化为矩阵方程。这里先举一个简单的例子。里先举一个简单的例子。例例1无限薄导体圆盘上的电荷分布问题。无限薄导体圆盘上的电荷分布问题。 试讨论半径为试讨论半径为a的无限薄理想导体圆盘，在中心线的无限薄理想导体圆盘，在中心线距离距离d处有一点电荷处有一点电荷 ，如图，如图5-17-1所示，求解导体圆盘所示，求解导体圆盘上的电荷分布。上的电荷分布。 解解 假设导体圆盘上电荷密度为假设导体圆盘上电荷密度为 ，根据电，根据电磁学的基本概念可知：磁学的基本概念可知：(1) 由外加电荷由外加电荷Q在导体圆盘上产生的电位在导体圆盘上产生的电位e 和导体圆和导体圆盘本身感应电荷密度盘本身

3、感应电荷密度所产生的电位所产生的电位i之和之和U 在盘上处在盘上处处相等，即保证导体圆盘是等位面。处相等，即保证导体圆盘是等位面。(2) 由于本问题中是感应电荷，因此总电荷由于本问题中是感应电荷，因此总电荷Qi0，其中，其中( ,)x y图图5-17-1导体圆盘上的电荷分布导体圆盘上的电荷分布 (5-17-1) (5-17-2) (5-17-3)22204eQxyd 0( ,)4isx ydSr ( ,)isQx y dS于是，问题可写为于是，问题可写为 (5-17-4)式中式中r= ，其中打撇的表示源点，不打撇的表，其中打撇的表示源点，不打撇的表示场点。示场点。 这个问题，采用电磁学经典解析

4、方法不能很好的解决，这个问题，采用电磁学经典解析方法不能很好的解决，因为未知量因为未知量 处于积分内部，是一个典型的积分方程。为此，处于积分内部，是一个典型的积分方程。为此，把圆盘分割成两部分：中心小圆和外部环带把圆盘分割成两部分：中心小圆和外部环带(如图如图5-17-1所示所示)，并并假定每一部分内的电荷密度假定每一部分内的电荷密度 (i=1,2)近似为常数近似为常数，于是，于是 (5-17-5)式中式中(5-17-6) 0()eiiUQ 约束条件22()()xxyyi211221( ,)()()ix yP SP S 1 ()0 iiiSSP SSS 称为脉冲函数，这时问题方程称为脉冲函数，

5、这时问题方程(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)内只有三个未知数内只有三个未知数 、 和和 ，于是，

6、于是方程组超定。方程组超定。()iP S21021 40 ieiiSiiidSUrS 12U 为了把为了把超定方程组超定方程组转化为转化为唯一解唯一解的方程组，可以采用很的方程组，可以采用很多办法。矩量法中，习惯用选配过程解决这个问题。简单说多办法。矩量法中，习惯用选配过程解决这个问题。简单说来，即在来，即在每个离散的单元每个离散的单元上上只选取一个场点只选取一个场点作为代表来建立作为代表来建立方程。例如，在方程。例如，在例例1中对于离散的中对于离散的 和和 分别取分别取 和和 两点做试验点，如图两点做试验点，如图5-17-2所示。具体写出方程组所示。具体写出方程组(5-17-9)其中其中11

7、1122121122221122 0 eellUllUSS 第1试验点第2试验点1S11(,)x y22(,)xy2S121122011122201121220221 (5-17-10)4()()1 (5-17-11)4()()1 4()()SSdSlxxyydSlxxyydSlxxyy122222022 (5-17-12)1 (5-17-13)4()()SSdSlxxyy图图5-17-2 圆盘上的试验点圆盘上的试验点111122121122221122 0 eellUllUSS 第1试验点第2试验点其中其中 表示表示 面元电荷面元电荷在在 处产生场的自作用单元；处产生场的自作用单元； 表示表

8、示 面元电荷在面元电荷在 处产生场的自作用单元；处产生场的自作用单元； 表示表示 面元电荷在面元电荷在 处产生场的互作用单元；处产生场的互作用单元； 表示表示 面元电荷在面元电荷在 处产生场的处产生场的互互作用单元。作用单元。1S1S2S2S11l22l12l21l11(,)x y11(,)x y22(,)xy22(,)xy121122011122201121220221 (5-17-10)4()()1 (5-17-11)4()()1 4()()SSdSlxxyydSlxxyydSlxxyy122222022 (5-17-12)1 (5-17-13)4()()SSdSlxxyy又有又有 (5-

9、17-14)经过经过离散化过程和选配离散化过程和选配过程，将积分方程组过程，将积分方程组(近似地近似地)转化为矩转化为矩阵方程阵方程 (5-17-15)由此得出电荷分布的解为由此得出电荷分布的解为 (5-17-16)122201122220221414eeQxydQxyd 111211212222121100eellllSSU 1111121221222121100eellllUSS图图 5-17-3 矩量法的一般过程矩量法的一般过程图图5-17-3所示的矩量法求解问题的一般过程。所示的矩量法求解问题的一般过程。讨论讨论 (1)矩量法的原问题并不限于积分方程，也可以是微分矩量法的原问题并不限于

10、积分方程，也可以是微分方程或其他方程。但必须能抽象成方程或其他方程。但必须能抽象成算子方程算子方程。从这一点而言，。从这一点而言，它是普遍的；另一方面，矩量法最终要转化为它是普遍的；另一方面，矩量法最终要转化为矩阵方程矩阵方程加以加以解决。因此，原问题必须属于线性算子范畴。例如，最速下解决。因此，原问题必须属于线性算子范畴。例如，最速下降线所构成的积分方程降线所构成的积分方程 不是线性泛函，所以无不是线性泛函，所以无法采用矩量法。法采用矩量法。(2) 电磁理论中计算的矩阵单元，一般均表示某个源在一个区电磁理论中计算的矩阵单元，一般均表示某个源在一个区域所产生的场，而实际产生的场往往都随着源的距

11、离增加而域所产生的场，而实际产生的场往往都随着源的距离增加而减少。换句话说，矩量法中矩阵一般是对角占优的：自作用减少。换句话说，矩量法中矩阵一般是对角占优的：自作用单元单元 比互作用单元比互作用单元 所起的作用要大。这一点在概念所起的作用要大。这一点在概念上十分重要。上十分重要。212yJdygynnl()mnlmn矩量法的研究对象是一般非齐次方程矩量法的研究对象是一般非齐次方程 (5-17-17)线性算子线性算子 的运算空间称为定义域，而的运算空间称为定义域，而 组成的空间称为值组成的空间称为值域。式域。式(5-17-17)中中 是已知的是已知的激励函数激励函数， 为为未知函数未知函数。令。

12、令 在在 的定义域内展开成的定义域内展开成 的组合，有的组合，有 (5-17-18)2.2 矩量法的一般过程矩量法的一般过程( )L ug( )L uguL12 ,nnuu uu即1NTnnnuuu其中其中1N为展开系数矩阵u 表示矩阵转置，应该注意到：表示矩阵转置，应该注意到：展开函数与基函数展开函数与基函数是有区是有区别的。一般来说，基函数是一无限展开。从完备基转化为近别的。一般来说，基函数是一无限展开。从完备基转化为近似有限截断基已经构成误差了，再从有限截断基转化为有限似有限截断基已经构成误差了，再从有限截断基转化为有限展开函数就很难保证展开函数就很难保证 能收敛于能收敛于 ，这也是矩量

13、法的，这也是矩量法的研究中需要深入研究的一个问题。这里且写出研究中需要深入研究的一个问题。这里且写出 (5-17-19)而而1Nuuu为展开函数1Nnnnuu1()NnnnL ug从算子方程从算子方程(5-17-17)到式到式(5-17-19)即构成离散化过程。它可以即构成离散化过程。它可以是是函数离散函数离散，也可以是，也可以是区域离散区域离散，或两者兼有。，或两者兼有。 T现在规定适当的内积现在规定适当的内积 。在算子。在算子L的值域内定义一类的值域内定义一类权函数权函数(或检验函数或检验函数) ,作用于式作用于式(5-17-19)两边，两边，且取内积，有且取内积，有 (5-17-20)

14、这就是所谓的选配过程或试验过程，矩量法的名称也由此这就是所谓的选配过程或试验过程，矩量法的名称也由此而来，即把激励矢量而来，即把激励矢量 和和 分别向权空间投影，取它分别向权空间投影，取它的矩，根据矩的大小确定展开系数。的矩，根据矩的大小确定展开系数。 如果展开函数的数目与权函数数目相等，则可把式如果展开函数的数目与权函数数目相等，则可把式(5-17-20)写成矩阵形式写成矩阵形式 (5-17-21)其中其中 (5-17-22)于是可以解出于是可以解出(5-17-23) ,g12,N 1, (),NnnnmmL ug1,2,mNg()nL umnnmlg, ()mnnmlL u1nmnmlg若

15、规定函数矩阵若规定函数矩阵 (5-17-24)于是待求的函数为于是待求的函数为 (5-17-25)矩量法的一般过程的数学表示如图矩量法的一般过程的数学表示如图5-17-4所示。所示。十分清楚，矩量法的结果优劣取决于：十分清楚，矩量法的结果优劣取决于：离散化程度；离散化程度； 和和 的选取；的选取；线性方程组的求解。在线性方程组的求解。在 = 的特的特殊情况下，可称为殊情况下，可称为Galerkin(伽略金伽略金)法法，于是矩量法也称于是矩量法也称为广义为广义Galerkin法。法。12,TnNuu uu1TTnnmnmuuu lgnunumm图图5-17-4 矩量法一般过程的数学表示矩量法一般

16、过程的数学表示 例例2研究研究 ，其中，其中解解 已经知道，此问题存在精确解已经知道，此问题存在精确解 本例采用矩量法求解，选择本例采用矩量法求解，选择 再选择权函数再选择权函数( )L ug222,14, (0)(1)0dLgx uudx 240511( )623uxxxx11,2,.(,)nnuxnxxN1muxxmm即采用即采用Galerkin法，内积定义为法，内积定义为 于是可给出一般计算结果于是可给出一般计算结果 10,( ) ( )gx g x dx 121120d, ()()()ddmnmnmnlL uxxxxxx 10(1)()d1nm nmnn nxxxmn 归纳起来有归纳起

17、来有 1120,()(14)dmmmggxxxx 11330(44)dmmxxxxx114(38)12242(2)(4)mmmmmm (38)(38),12(2)(4)mnmmmmmlgmnmm情况情况1：N=1111111111,33010lg于是有于是有21( )()10u xxx情况情况2: N=2121111323014725121241111521030601127233122323122312( )()()10330103u xxxxxxxx情况情况3: N=3123113113253014712512513917057363336191det( )1051010125328105

18、00l 12313111135705030233617105000705253012151111370503060 2424011511( )()()( )23623u xxxxxxxxux十分明显，十分明显，N=3时已得到了精确解时已得到了精确解 。 矩量解的曲矩量解的曲线如图线如图5-17-5所示。所示。 0( )ux( )u x图图5-17-5 u (x)矩量解矩量解第二章第二章 矩量法矩量法(Method of Moment)2.3 选配和离散过程选配和离散过程2.3.1 点选配点选配2.3.2 脉冲分域基脉冲分域基2.3.3 三角形函数分域基三角形函数分域基2.4 算子研究算子研究2

19、.4.1 近似算子近似算子2.4.2 扩展算子扩展算子2.4.3 微扰算子微扰算子2.3 选配和离散过程选配和离散过程 从上面的典型例子可知，矩量法的精华在于选配和从上面的典型例子可知，矩量法的精华在于选配和离散过程，值得单独进行研究。离散过程，值得单独进行研究。2.3.1点选配点选配 点选配是一种最简单而最典型的选配函数。因为矩阵点选配是一种最简单而最典型的选配函数。因为矩阵单元为单元为 ，一般说来，其中所含的积分计算，一般说来，其中所含的积分计算十分困难，这种情况下，最简单的办法是做某些点的投影，十分困难，这种情况下，最简单的办法是做某些点的投影，即所谓的点选配，实际上相当于把即所谓的点选

20、配，实际上相当于把权函数权函数取为取为 函数。函数。, ()mnmnlL uDirac 例例3 任研究任研究 。解解 设设 ，在这个例子中取在这个例子中取 函数为权函数即函数为权函数即 其中，其中， 是这个问题的选配点，于是有是这个问题的选配点，于是有222d( )14, (0)(1)0du xx uux 1( )nnuxxx(),1,2,.,1mmmxxmxmNNmx 例例3 任研究任研究 。解解 设设 ，可得到，可得到 在这个例子中取在这个例子中取 函数为权函数即函数为权函数即 其中，其中， 是这个问题的选配点，于是有是这个问题的选配点，于是有222d( )14, (0)(1)0du xx

21、 uux 1( )nnuxxx21221d()1 4dNnnnxxxx (),1,2,.,1mmmxxmxmNNmx121 (1)14Nnnnn nxx 1110, ()(1)()d(1)()1nnmnmnmmlL un nxxxxn nN 1220,(14) ()d14()1mmmmggxxxxN 归结起来，可写出归结起来，可写出1(1)()1nmnmln nN214()1mmgN 情况情况1: N=1111122,2( )u xlgxx 情况情况2: N=21213229242591213142191822262493 2323121312( )()()18318183u xxxxxxxx

22、情况情况3: N=312333522442332927132244可以得出可以得出3121413452921549427427321023631292)(312165)(31)(21)(04242xuxxxxxxxxu对于点选配情况对于点选配情况N=3，又一次回复到精确解。，又一次回复到精确解。讨论讨论 (1) 对于点选配的情况，对于点选配的情况，N+1阶矩阵中的阶矩阵中的N阶主子阵并不阶主子阵并不等于在等于在N时的系数矩阵时的系数矩阵(和和Galerkin情况不同情况不同)。因此当。因此当N逐逐渐变大时计算量无法节约。渐变大时计算量无法节约。(2) 点选配虽然看起来非常简单，然而其内在的道理

23、极其点选配虽然看起来非常简单，然而其内在的道理极其深刻。这一点可以从数值积分看出。研究表明任何数值积深刻。这一点可以从数值积分看出。研究表明任何数值积分方法，不论矩形、梯形、二次样条等，说到底都是选择分方法，不论矩形、梯形、二次样条等，说到底都是选择积分区域的点和区域点所对应的系数，由此产生积分区域的点和区域点所对应的系数，由此产生Gauss积积分的思想。所以在矩量法中，研究最佳点选配将是一个十分的思想。所以在矩量法中，研究最佳点选配将是一个十分有意义的课题分有意义的课题 。2.3.2脉冲分域基脉冲分域基 矩量法在离散化过程中用展开函数取代基函数，带来矩量法在离散化过程中用展开函数取代基函数，

24、带来了方便和自由。但是，随之而来的如何确保解的收敛性的了方便和自由。但是，随之而来的如何确保解的收敛性的问题却值得人们重视。问题却值得人们重视。 在尚未了解在尚未了解u(x)函数性态的条件下，采用有限个展开函数性态的条件下，采用有限个展开函数函数ui(x)，i=1，2，.，N时要确保解收敛显然在理论上时要确保解收敛显然在理论上存在不少困难，采用分域基函数可以说是比较稳妥的一种存在不少困难，采用分域基函数可以说是比较稳妥的一种解决方案。因为大多数良态函数解决方案。因为大多数良态函数(不做高速振荡不做高速振荡)均可以采均可以采用有限段直线或样条加以逼近，如图用有限段直线或样条加以逼近，如图5-17

25、-6所示。所示。图图5-17-6 分域基函数近似分域基函数近似 下面从最简单的脉冲函数着手展开讨论。下面从最简单的脉冲函数着手展开讨论。 一般的脉冲函数可以表述为一般的脉冲函数可以表述为112(1)2(1)112(1)2(1)(,)1()0(,)iiNNiiiNNxxxP xxxxx(5-17-26) 式式(5-17-26)表示以表示以i为中点，密度为为中点，密度为1/(N+1)的脉冲函数，的脉冲函数，在实际情况下，密度可以根据问题灵活改变，如图在实际情况下，密度可以根据问题灵活改变，如图5-17-7所示。所示。图图5-17-7 脉冲函数脉冲函数 图图5-17-8 三角形函数三角形函数 2.3

26、.3三角形函数分域基三角形函数分域基三角形函数也是常用的一种分域基，如图三角形函数也是常用的一种分域基，如图5-17-8所示。所示。 若采用三角形函数展开未知函数若采用三角形函数展开未知函数(x)，则有，则有 (5-17-27)所得的所得的解的合成相当于折线连接解的合成相当于折线连接，分段三角形函数所得的，分段三角形函数所得的折线包络如图折线包络如图5-17-9所示。所示。 为了研究具体例子，这里先给出三角形函数的导数概为了研究具体例子，这里先给出三角形函数的导数概念。引入如图念。引入如图5-17-10所示的阶梯函数所示的阶梯函数H(x-xi)，其定义为，其定义为)(01)(1, 11, 11

27、iiiixxxxixxxxxxxxTiii)()(1nNnnxxTxu图图5-17-9 分段三角形函数所得的折线包络分段三角形函数所得的折线包络图图5-17-10 H(x-xi) 函数函数 121()0iiiixxH xxxxxx (5-17-28) 再引入大家熟悉的再引入大家熟悉的Dirac-函数，也即脉冲函数，其定义函数，也即脉冲函数，其定义为为 (5-17-29)()0iiixxxxxx如图如图5-17-11所示。所示。图图5-17-11 (x-xi)函数函数 Dirac-函数有两个重要的性质：函数有两个重要的性质：1.归一性归一性 (5-17-30)2.选择性选择性 (5-17-31)

28、 这里不加证明的给出这里不加证明的给出Dirac-函数和阶梯函数之间的函数和阶梯函数之间的重要关系。重要关系。 (5-17-32)有了以上基础就可以把三角形函数的导数用阶梯函数有了以上基础就可以把三角形函数的导数用阶梯函数H表表示，具体为示，具体为 (5-17-33) ()()iidH xxxxdx1111111111()()()()()iiiiiiiiiiiiixxxxxxxxdT xxdxH xxH xxH xx()1ixx dx( ) ()( )iif xxx dxf x图图5-17-12给出形象的几何表示。给出形象的几何表示。图图5-17-12 三角形函数导数的几何表示三角形函数导数的

29、几何表示 例例4 重新研究重新研究Harrington(哈林登哈林登)问题，问题，L(u)=g，其中，其中L= ，g= ，边界条件为，边界条件为u(0)=u(1)=0。试用以三角函数作为展开函数，脉冲函数作为权函数的矩试用以三角函数作为展开函数，脉冲函数作为权函数的矩量法求解。量法求解。22ddx241x 解解 根据要求可写出根据要求可写出 于是有于是有 上式已计及上式已计及 选择权函数选择权函数于是矩阵单元于是矩阵单元 上式要分三种情况讨论。上式要分三种情况讨论。)()(1nNnnxxTxu2211()()(1)()2 ()()dnnnnndxL T xxT xxNxxxxxx 1)(Nmm

30、mmxxxP10()mnmnlL T xxdx1110()(1)()2 ()()mnnnP xxNxxxxxxdx11111iiiiNxxxx此外，激励单元为此外，激励单元为 结果可归纳为结果可归纳为12(1)12(1)1122200(14)()(14)(14)mNmNxmmmxgxdxP xxxdxxdx110) 1() 1(2nmnmnmNNlmn2132411(1)1mmNNg1mmNx10()mnmnlL T xxdx1110()(1)()2 ()()mnnnP xxNxxxxxxdx情况情况1：N=1 考虑到对比：考虑到对比： 则有则有 和和 的对比如图的对比如图5-17-13所示。

31、所示。114,L 12524g 250.26041696i 240511623uxxxx01130.2708333248u 1ux0( )u x图图 5-17-13 和 1ux0( )u x情况情况2：N=2l= g= 容易得到容易得到 同样对比有同样对比有 和和 的的 对比图如图对比图如图5-17-14所示。所示。633640127763271263401367621874680.213991715760.26237442187001530.2181073243650.267489722433uu 2ux0( )u x图图 5-17-14 和 2ux0( )u x情况情况3：N=3 于是有于

32、是有 同样对比有同样对比有 和和 的的 对比图如图对比图如图5-17-15所示所示840484048l 614819744815748g123483216611326432974915216324815785440.17382811131840.268229149152116160.2363281000141350.1757812112080.270833327681830.238281234uuu0( )ux 3ux图图 5-17-15 和 3( )u x0( )u x 讨论讨论 分域基在分域基在N不大的情况下与精确解的差距是明显的。不大的情况下与精确解的差距是明显的。但是它的相应矩阵是三条

33、带矩阵，可较明显地缩小计算量。但是它的相应矩阵是三条带矩阵，可较明显地缩小计算量。因此选择因此选择N不大的分域基并进行顶点拟合将会是一个比较不大的分域基并进行顶点拟合将会是一个比较好的方案。好的方案。2.4 算子研究算子研究 算子方程是矩量法建模的关键。它应该有两个方面算子方程是矩量法建模的关键。它应该有两个方面的要求：的要求： 一方面算子方程必须符合物理一方面算子方程必须符合物理(或工程或工程)问题的主要本问题的主要本质；质； 另一方面它又必须适合数值计算。另一方面它又必须适合数值计算。 这两个方面构成了算子研究的基础。这两个方面构成了算子研究的基础。2.4.1 近似算子近似算子 细心的读者

34、一定会提出这样一个问题，即细心的读者一定会提出这样一个问题，即例例4中为什么不采用中为什么不采用脉冲函数作为分域基展开？其实原因十分简单，因为脉冲函数的二脉冲函数作为分域基展开？其实原因十分简单，因为脉冲函数的二阶导数表示有很大困难。但是，倘若引进近似算子的概念，则可以阶导数表示有很大困难。但是，倘若引进近似算子的概念，则可以较好地解决这个问题。较好地解决这个问题。 算子近似含义相当广泛。作为例子，可采用有限差分代替微分。算子近似含义相当广泛。作为例子，可采用有限差分代替微分。 1()()22du xxxu xu xdxx 22211()()()2()22d u xxxu xu xu xxu

35、xu xxdxxx例例5 研究研究 的的Harrington问题，即问题，即 , ,试采用差分近似算子试采用差分近似算子 ,脉冲脉冲展开点选配的矩量法求解。展开点选配的矩量法求解。（做一般了解）（做一般了解）解解 为确保为确保 的边界条件，在两端各留出半的边界条件，在两端各留出半段为强制零段。因此当选择段为强制零段。因此当选择N个脉冲函数时，全部区域个脉冲函数时，全部区域(0,1)应分成应分成(N+1)段。即段。即 于是有于是有 L ug22,dLdx 214gx 010uuLL 010uu11xN 21112 ( )()11aLuNu xu xu xNN且做点选配有且做点选配有mmxx1mm

36、xN这样可以获得矩阵单元这样可以获得矩阵单元 的表示式的表示式,m nl12,n010P1112111am nmmlLdxNnnnP xP xP xxxdxNNN1120014mmmggdxxxxdx可以归纳为可以归纳为 22,211101m nNmnLNmnmn 2141mmgN 情况情况1：N=1 =8 ， 于是得到于是得到 对比对比 这里的这里的 和和 的对比如图的对比如图5-17-16所示所示 表面看来，与图表面看来，与图5-17-13类似，实际上脉冲函数和三类似，实际上脉冲函数和三角函数意义有很大不同，又注意到图角函数意义有很大不同，又注意到图5-17-16中中 和和 各强制置零半段

37、。各强制置零半段。1,1L12g 110.254a 01130.270833248u 1ux0( )ux104，31.04，图图5-17-16 和和 1ux0( )ux情况情况2: 于是有于是有 对比对比 和和 如图如图5-17-17所示。所示。1899 18l131259g121891319182521874590.209 876515670.259 259221870010.218 10730.267 489 723uu0( )ux 2ux图图5-17-17 和和 0( )ux 2ux2N 情况情况3: 于是有于是有 作为对比有作为对比有 和和 的对比如图的对比如图5-17-18所示。所示

38、。3N 3216016321601632l 518413g12376851225651512 10245128655362565127681311 2640.171 875117 4080.2656256553615 3600.234375000140.175 781210.270 833320.238281 234uuu 3ux0( )ux图图5-17-18 和和 3ux0( )ux2.4.2扩展算子扩展算子 算子包括定义域和运算域。如同数学上经常所做的那算子包括定义域和运算域。如同数学上经常所做的那样，可以采用扩展算子来增加展开函数或权函数选择的自样，可以采用扩展算子来增加展开函数或权函数

39、选择的自由度。由度。 原算子和扩展算子的逻辑关系如图原算子和扩展算子的逻辑关系如图5-17-19所示。很所示。很明显，扩展算子不改变原算子的运算。明显，扩展算子不改变原算子的运算。图图5-17-19 原算子和扩展算子的逻辑关系原算子和扩展算子的逻辑关系 例例6 希望希望Harrington问题问题 采用脉冲函数作为展开采用脉冲函数作为展开函数的并引入扩展算子概念。函数的并引入扩展算子概念。（做一般了解）（做一般了解）解解 从上面论述中已知从上面论述中已知 在原来的定义域中不存在。但在原来的定义域中不存在。但深入研究矩量法后发现，矩量法并不要求深入研究矩量法后发现，矩量法并不要求 有定义，而有定

40、义，而只要内积只要内积 有定义即可。有定义即可。 (5-17-36)于是可以放松要求为，所选择的权函数于是可以放松要求为，所选择的权函数 满足定义域，即满足定义域，即 (5-17-37)则可引入扩展算子则可引入扩展算子 (5-17-38) L ug L P L P ,L u 2112001100,dduL uu dxddxdxduddudxdxdxdx 010 10,dduL udxdxdx从而避免从而避免 问题，于是设问题，于是设 L u1Nnnnmmua P xxT xx则可知矩阵单元为则可知矩阵单元为1010,111211121212121mnmnmndT xxdP xxlL P xxd

41、xdxdxmmmNH xH xH xNNNnnxxdxNN以及激励单元以及激励单元 1120011331300013110m33111144443334111322441133mmmmmmmmmmNmNggdxT xxxdxT xxd xxT xxxxxxdT xxNxxP xxxP xxxdxNxxdxNxx m+111NmNdx归纳起来是归纳起来是 它和三角函数展开脉冲函数检验所得到的公式差距极其细它和三角函数展开脉冲函数检验所得到的公式差距极其细微。微。当当N增大时，增大时， 彼此相当接近。彼此相当接近。情况情况1: 这种情况与这种情况与 完全吻合。完全吻合。211101mnNmnlNmnmn 22241