MSOR迭代及MSOR波型松驰并分析.pdfVIP

．谚35678上 MSOR迭代及MSOR波型松弛 叶 宁 山东大学数学与系统科学学院应用数学系 摘 要 本文由七部分组成： 从第一部分到第四部分，我们考虑非奇异线形系统 Ax=b (1．1) 其中A∈C…和b∈C”，并且假设A是一个2X2的块矩阵形式： 4也：。12] 令 。=(台三卜[三。：]u=(：台卜㈡Z，] A相应的Gauss—Jacobian阵为： B=D。眇+三) MSOR迭代阵形式为： L。=(D—nz)。((，一n)D+flU) =(幺爿(卜静。譬诲。] 现在，我们寻求B的特征值满足何种要求时，能使三。迭代收敛因此得到 如下结论： 定理一：系数阵是一个二循环阵的MSOR方法收敛当且仅当相应的Gauss— Jacobian阵的特征值一和MSOR方法松弛因子q，m：满足： (i)103。-111∞：-1I1 (ii)所有的∥2位于椭圆内部： p 口bo≤o2厅 2=d+P14(acos8+ibsin口) 其中，中心点d：竺!±竺二! 国1∞2 卯】∞2 口：鹕姬丑回 a=掣 国l∞2 椭圆通过(1，0)点。 定理二．如果所有对应于Gauss—Jaeobian矩阵的谱半径平方／a2位于一个椭圆 内，其中心为d长半轴为a，短半轴为岛长半轴与∥2的实轴夹角为口且此椭 圆不包含 (1，O)。那么，在使MSOR方法收敛的最佳松弛因子∞。，国：的形式 为 ‰妒半兰篙等竽 1一d+、／(1一d)2一c2P丑。 并且MSOR迭代阵的谱半径表达式如下： 铲属I网2面杀舞司 从第五部分到第七部分主要研究了如下形式的方程： (2．1)对应的MSOR迭代方程进行付氏变换后，相应的迭代矩阵日。为 Ⅳ。：(D。一彬。，丘。，)_1【(，一哦。，)D。+蛾。，q。，】 2 日。=(D。一哦。，厶。，)-1№一哦。，)D。+哦。，U。。J 其中 Do=iD．／+D 得如下结果： L=D。。@+工) 应用一到四部分的结果，可得到J。在满足一定条件下H。对应的最佳松弛因 子： ‰吁丝等蔫裂笋 1一d(oj+√Q—qo)厂一c{o)。P。酬“’ 把其结果应用于波形方程： 竺一d：罂：， ’ a 瓠i 关键词：MSOR迭代、复松弛因子。 引 言 SOR方法是解决线形系统Ax=b的一种非常流行的数值化方法。把矩阵A 分解为A…DL(，，其中D是对角阵，￡和∥分别是严格的下三角和上三角 阵。SOR方法的迭代形式可表达为： 工㈨=(D一础)。1“1一国)D+roU)X[‘‘1]+国(D一础)’1b 则对应方的SOR迭代阵形式为： H。=CD—oz)。((1一∞)D+caU) Gauss—Jacploarl迭代阵形式为： B=D。1(L+u) 其中，对于SOR方法一个很重要的问题就是如何确定SOR方法中的松弛因 子co的值。对于松弛因子的研究已有很广泛的研究。 阵日。的特征值之间有一个严格的关系： 氕|+∞一1=H|∞