毕业论文

打赏
当前位置: 毕业论文 > 电子通信 >

有限元方法在波导计算的应用仿真(3)

时间:2017-02-07 21:38来源:毕业论文
当然,我们希望尽可能地用解析方法求解边值问题。而许多实际重要的工程问题都没有解析解。为了克服这种困难,人们已发展了各种近似方法,其中应用


当然,我们希望尽可能地用解析方法求解边值问题。而许多实际重要的工程问题都没有解析解。为了克服这种困难,人们已发展了各种近似方法,其中应用最广泛的是里兹方法和伽辽金方法。
1.3.2里兹(Ritz)变分法[1]
一种用变分表达式(泛函)来得到近似解的方法。这种方法将边值问题用泛函表示,泛函的极小值对应于给定边界条件下的控制方程的解。为了说明这种过程首先定义内积
              =                               ( 2 )
星号表示复数共轭。在这种内积定义下,如果有
                                        ( 3 )
则表示(1)中的算符£是自伴的,如果有
                                       ( 4 )
则(1)式中运算符£是正定的。可以证明。如果(1)式中运算符£既自伴又正定,那么,(1)式的解可通过求下式泛函对  的极小值得到
F( )=(<£ , >/2)-(< , >)/2-(< , >)/2         ( 5 )
 表示试探函数。
一旦确定了泛函,即可用下述步骤来解决,假定(5)式中的 近似展开为
                           ( 6 )
其中, 是定义在全域上的展开函数,  是待定的展开系数,{.}表示列向量。上标T表示向量的转置。(6)带入(5)
F=                 ( 7 )
为了求F( )极小,我们令其对 的偏导数为0
                                                ( 8 )
可将其写成下列矩阵方程
              [S]{c}={b}                             ( 9 )
   [s]的元素为
        =(1/2)                   ( 10)
     [b]的元素为
                                       ( 11)
显然,[S]是对称矩阵。引用运算符£的自伴性质,  可写成
                                         ( 12)
求解矩阵方程(9)式,即可得到(1)式的解。
1.3.3 伽辽金法[1]
假定 是控制方程( 1)的近似解。将 代入式( 1),由于与真实解之间存在误差,因此会得到一个非零余量
                                    ( 13) 有限元方法在波导计算的应用仿真(3):http://www.youerw.com/tongxin/lunwen_2635.html
------分隔线----------------------------
推荐内容