msor迭代及msor波型松驰

msor迭代及msor波型松驰

ID:27075975

大小:4.45 MB

页数:6页

时间:2018-12-01

msor迭代及msor波型松驰_第1页
msor迭代及msor波型松驰_第2页
msor迭代及msor波型松驰_第3页
msor迭代及msor波型松驰_第4页
msor迭代及msor波型松驰_第5页
资源描述:

《msor迭代及msor波型松驰》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库

1、.谚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方法收敛当且仅当相应的Ga

2、uss—Jacobian阵的特征值一和MSOR方法松弛因子q,m:满足:(i)103。-111∞:-1I<1(ii)所有的∥2位于椭圆内部:p2=d+P14(acos8+ibsin口)口>bo≤o<2厅其中,中心点d:竺!±竺二!国1∞2卯】∞2口:鹕姬丑回国l∞2椭圆通过(1,0)点。定理二.如果所有对应于Gauss—Jaeobian矩阵的谱半径平方/a2位于一个椭圆内,其中心为d长半轴为a,短半轴为岛长半轴与∥2的实轴夹角为口且此椭圆不包含(1,O)。那么,在使MSOR方法收敛的最佳松弛因子∞。,国:的形式为‰妒半兰篙等竽并且MSOR迭代阵的谱半径表达式如下:

3、从第五部分到第七部分主要研究了如下形式的方程:(2.1)对应的MSOR迭代方程进行付氏变换后,相应的迭代矩阵日。为Ⅳ。:(D。一彬。,丘。,)_1【(,一哦。,)D。+蛾。,q。,】2a=掣铲属I网2面杀舞司1一d+、/(1一d)2一c2P丑。日。=(D。一哦。,厶。,)-1№一哦。,)D。+哦。,U。。J其中Do=iD./+D(2.1)式对应的Gauss—Jacobian迭代,进行付利叶变换后,相应的迭代矩阵可得如下结果:L=D。。@+工)应用一到四部分的结果,可得到J。在满足一定条件下H。对应的最佳松弛因子:1一d(oj+√Q—qo)厂一c{o)。P。酬“’把其结果应用

4、于波形方程:罂’关键词:MSOR迭代、复松弛因子。‰吁丝等蔫裂笋竺一d::,a瓠i引言SOR方法是解决线形系统Ax=b的一种非常流行的数值化方法。把矩阵AL(,,其中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的值。对于松弛因子的研究已有很广泛的研究。对于具

5、有矩阵A特征的二循环阵,Gauss—Jacobian迭代阵B和SOR迭代阵日。的特征值之间有一个严格的关系:氕|+∞一1=H|∞0AJ∥,和五,分别是B和Ⅳo的特征值a如果Gauss—Jacobian所有的特征值位于(一1,i),那么,对应于SOR方法的最佳松弛因子为:2国印f=——1==2=亍1十√1一∥o。而对于依赖于时间t的∞…的CSOR迭代阵为:(d+D)州旷三戊¨㈨h№=(鲁+。)x卜qt,+[u一(罢+。)]J:》(『)x。阽_1(f—rpr经付氏变换后,SOR迭代阵为:JH㈣=【D。q侧£r1【(1,删)D。呐删U】分解为A…D相应的Gauss-Jacobi

6、an阵为:B㈧=Dn-,@+£),如果假设B的特征值位于如下椭圆内:∥=P’“[acose+ibsin口】a≥b0≤口<2,2-则H:的最佳松弛因子表达式为2=l+√1一c2P以上是前人经过研究已有的部分结果。古典的SOR和MSOR方法的最佳松弛因子是实数。因此,收敛域是严格受限的~个区域。在本文中,我们将在复数域上考虑YSOR方法的最佳松弛因子和收敛区域。同时,发现6auss—Jacobian阵的特征值分布的收敛域要比实数最佳松弛因子收敛域大的多。MSOR方法要优于SOR方法。而MSOR方法的收敛域要比SOR方法的收敛域大的多。且MSOR方法的收敛速度也优于SOR方

7、法。6山=————{======2Ⅲ一感谢您试用AnyBizSoftPDFtoWord。试用版仅能转换5页文档。要转换全部文档,免费获取注册码请访问http://www.anypdftools.com/pdf-to-word-cn.html

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。