预条件IMGS迭代方法的比较定理.pdfVIP

数学物理学报 http：／／actams．wipm．ac．cn 预条件 IMGS迭代方法的比较定理 李继成 蒋耀林 (西安交通大学理学院 西安 710049) 摘要：该文通过比较几种迭代矩阵的谱半径，首次证明了经预条件子 ，+＆ 的 IMGS迭代 方法 L1-21比经典的AOR迭代方法收敛速度快；其次证明了预条件 SOR方法比经典的SOR 方法收敛速度快．最后，证明了有关迭代矩阵的谱半径关于参数的单调性，改进了一些已知的 结论 ． 关键词：预条件子；Gauss—Seidel迭代；AOR迭代；M 一矩阵；线性系统． MR(2000)主题分类：65F10；65F15；15A09 中图分类号：O241．6 文献标识码：A 文章编号：1003—3998(2011)04—880—07 1 引言 本文考虑如下线性系统 Ax=b， X，b∈R ． (1．1) A=(aij)∈R ”为一个非奇异矩阵．为了提高求解系统 (1．1)的迭代方法的收敛速度，通 常将系统 (1．1)等价地转化成预条件系统 P =Pb， (1．2) 非奇异矩阵P称为一个预条件子．最近，通过采用不同的预条件子 P，求解预条件系统 (1．2) 的多种迭代方法被提出 [1--3,6--11]．具体来讲，对于预条件子 P=I+ ，其中 t啦㈡一1’……’ 3， OL。∈(0，1]1i=1，2，… ，n一1．文献 [1]中作者从理论上证明：若 是非奇异M 一矩阵，则 IMGS方法 [1-2]不但 比经典的 SOR方法收敛速度快，而且比经典的Gauss—Seide方法收敛 速度也快．本文从理论上证明：对于预条件系统 (1．2)，IMGS方法比经典的AOR方法收敛 速度也快，而且，由于经典的SOR方法是AOR方法在 =r时的特殊情况，所以，本文推 广了一些已有结论．此外，本文也证明了预条件SOR方法比经典的SOR方法收敛速度快． 本文第二部分给出几个概念和引理，第三部分证明了几个主要的比较定理和关于松弛 因子u的单调性定理，第四部分给出数值例子来支持本文主要理论结果． 收稿 日期：2009一Ol一18；修订 日期：2010—03—25 E—mail：jcli@mail．xjtu．edu．cn 基金项 目：国家自然科学基金 、陕西省自然科学基金 (2007A16)和中央高校基本科研业务费专项 资金资助 No．4 李继成等：预条件 IMGS迭代方法的比较定理 881 2 预备知识 对于矩阵A∈R ，如果 A=sI—B，B 0，sp(B)，其中p(B)表示矩阵B的谱半 径，则称 为非奇异M 一矩阵． 如果M 是非奇异M 一矩阵，N 0，则称A=M —N为 一分裂．如果 -1 0且 M -1N 0， 则称A：M —N为弱正则分裂．显然 M 一分裂是弱正则分裂． 全文假设 A = (aij) 是非奇异不可约 一矩阵，而且满足 0 ai,i+la冲l，i 1， i=1，2，… ，n一1．考虑分裂 A…I L U， (2．1) 『为单位矩阵，一 和 一 分别表示 的严格下三角部分和严格上三角部分做成的相应阶 数的矩阵．求解系统 (1．1)的经典的SOR方法的迭代矩阵为 L =(I～coL)～ 【(1～ )+ ] (2．2) 为非 0的实数．经典的AOR方法的迭代矩