H-矩阵的预条件AOR迭代法及其收敛性.pdfVIP

第24卷第6期 齐 齐 哈 尔 大 学 学 报 Vo1．24．No．6 2008年 11月 JournalofQiqiharUniversity Nov．，2008 H一矩阵的预条件AOR迭代法及其收敛性 张保祥 、 (淮阴师范学院 数学系，江苏 淮安 223300) 摘要：给出了H一矩阵的预条件AOR迭代法及其收敛性，并给出了松驰因子 与加速因子 的选取对收敛速度的 影响，同时通过数值实例验证了主要结果。 关键词：谱半径；AOR迭代法；预条件；H一矩阵 中图分类号：0151．21 文献标识码：A 文章编号：1007—984X(2008)06—0069—03 1 预备知识 冈为非奇异H一矩阵主对角元非零，所以本文总假定所涉及矩阵主对角元非零，并设A=(口，，)∈R 为 阶实方阵。 不失一般性，令A= 一 — ，其中一 ，一 为A的严格下三角与严格上j：角矩阵，则基本Gauss—Seidel 迭代法的迭代矩阵为 = (，一 ) 基本AOR迭代法的迭代矩阵为 L =(Jr—aL)叫[(1一 Jr+改 】 基本AOR迭代法的迭代矩阵为 L =(1一rt) [(1一 J『+(一 三+吐 ] 定义 1I1lA=()∈R 称为是z一矩阵，若 ≤0，f≠J，f，J=1，2，…，。 定义2 若A=()∈R舳是z一矩阵且有分解式Z=sl—B，使得 p(B)，则称A为M一矩阵。 定义3t”称矩阵A=(m)∈R脚为A的比较矩阵，其中mⅣ=l口，I，m =一I I，f≠J，f，=l，2， … 以 。 ， 定义4l【1令A=(a，，)∈R ，若A的比较矩阵A是M一矩阵，则称 是H一矩阵。 定义5” 矩阵分裂A=M—N称为 1)M一分裂，如果 是M一矩阵且N 0； 2)弱正则分裂，如果 存在且 ．0， N 0； 3)正则分裂，如果 存在且M 0，N 0。 引理 若A：M—N为弱正则分裂，那么p(M Ⅳ)l甘 A 0。 引理2l【J设 是 H一矩阵，则I 一‘lM 一。 引理3。 若A O且为不可约刀X刀矩阵，则有 1) 有一正特征值等于谱半径p(A)； 2)相应于 ()存在一特征向量 0。 引理4I”般 A 0。则 1)若存择 0， ≠0满足Ax≥锨，则 ()≥ ，进一步地，若Ax锨，则p(A) ； 2)若存在 0， ≠O满足Ax ，则p(A)≤ ，进一步地，若Ax ，则p(A) 。 收稿日期：2008-07-07 作者简介：张保祥 (1965一)．男，江苏沭阳人，讲师．硕士．主要从事数值代数方而的研究，Emaihhabao@163．tom。 · 7O· 齐 齐 哈 尔 大 学 学 报 2008钲 2 主要结论 取预条件算子：P= =I+S对 =(a)∈R 是H一矩阵作用 0 一 lal2 0 … O 0 0 。 o~2a23 … 0 = 0 0 0 … 一 1口 1