关于不定线性方程组的若干预处理子的注记.pdfVIP

河南师范大学学报(自然科学版〉第 36 卷第 4 期 2008 年 7 月 }ou俨nalof Henan Normal Uηiversity (Natural Science) 文意编号: 1000-2367(2008)04一 0014…03 Vol.36 No.4 July.2008 关于不定线性方程组的若干预处理子的注记 孙丽英1 ，薛占熬2 (1.广东教育学院数学系，广州 510303 ， 2. 河南师范大学计算机与信息技术学院，河南新乡 453007) 摘 票:利用 Schur 分解，提出 KKT 现实不定线性系统的若干预处理子，讨论了这些预处理情形下的 Krylov 子空间方法收敛所需的法代步数，从而说明这些预处理方法是非常有效的. 关键词:不定矩阵 E预处理Tl极小多项式; Krylov 子空间方法 中固分提号 :0152.21 文献栋识码 :A 考虑如下 KKT 型实不定线性系统 ,A B T 、 ， x\ ， 1 、 (BTy +Ax = 1 u=(-_ \( )=(J \, NP!_ , (1) 飞 C 0 / \y) 飞 g/ \Cx g b 这里 A E RnXn ,B ,C E RmXn ， n 注 m. 当\]!来自约束变分或最优化问题时，常有 B =C;也有些问题要求 A 是 对称的，进而 ψ也是对称的.但不论 ψ对称与杏，通常\]!都是不定的(即其特征根的实部可正可负) . 形如(1)的线性方程组的选代解法有很多，特别是 Krylov 子空间方法(比如，应用于对称线性方程组的 MINRES 法，以及应用于非对称情形的广义极小残最 CGMRES) 怯).这略方法的共同特点是经过不姐过 η十m次选代可得到钱性系统(1)的精确解，但对于大型线性系统来说，这个选代次数的上限常常变得不实 际.这就需要预先处理(1)的系数矩阵，改善其谐性质，把解线性系统(1)所需的迭代总次数控制在尽可能小 的范围内. t也就是说，利用预处理子 P，对形如(1)的线性方程组进行预处理，减少收敢所需要的选代步数，同 时在每次选代时不1?著增加计算最. 为此，目前国内外许多学者致力于研究形式各界的预处理子，分析其收敛所需的选代步数、收缴濡度等 问题，这巳成为计算数学领域的一个热门课题.好的预处理子的必要条件是选代矩阵 T=P… 1\]! 有次数较低 的极小多项式，进而其相关的 Krylov 子空间的维数也就越小，该跚处理方法收敛所槽的选代步数就越少.这 个条件~T是非退化的或者其 lordan 标准型有小维数的 Jordan 块时，往往表现为 T有几个不阔的特征棍. 在 Schur 分解的基础上，考虑关于(1)的如下形式预处理子 z ,A 0 、 (A B T ) (A B T ) P 1 = PQ = (. _.. \ ,Pz = I _. . __1 P3 = 1 _ . __1 \0 CA…l J-. 10 士 CA…1 BT 1 . - 1 C 2CA …lBTI 证明了这些预处理子的迭代矩阵至多有 4 个不问特征根，即得到关于迭代姐阵 T=P叫 1\]! 至多是 4 次的极 小多项式.从而 Krylov 子空间方法(例如 GMRES 方法)应用于系数短阵是 T=P- 1 \]! 的预处理线性方程组 时，在 4 步戒更少步数选代内是收放的.即利用这蜘预处理子，对于 Krylov 子空间方法来说，烹多遮代 4 次 就会终止航线性系统(1)的精确解上，说明这蜡方法对于解决战性系统(1)是非常有效的. 1 定义与符号 Krylov 子空间方法是目前常用的最重要的选代方法之一，它把正交投影的思想应用到 Krylov 子空间 收稿日期 :2007-09-12 基金项目:广东省自然科学基金(7004344); 广东省高等学校自然科学研究建点项目 (052026) , 广东教育学院中青年学术 骨干培养项目 作者简介:孙丽英(1965 …) ，女，河南平顶山人，广东教育学院教授，研究方向 s数值代数、非负矩阵论. 第 4 期 孙商英等 z 关于不定线性方程组的若干预处理子的法记 15 ，，(矶的芸芸 span忡，胁 ，lf b ， …， ψ..-Ib ， …}上.其最常见的方法有共辄梯脏怯(CG) 、最小残量法(MINRES) 、 广义最小础盘法(GMRES)等(参辛苦文献[1.2J). 设fI E R(时酬〉圳，什时，满足 f(ψ) = 0 的次数最小的首一多项式称为 ψ 的极小多项式. 设fI E C(时时X(肿时，则存在四姐阵 U ，使得 U 怜 w S.该分解式称作 fI 的 Schur 分解.其中 S 角炬阵，ψ称为的 Schur 标准形;而且适当班取 U ，可使 S 的财角元紫按任意指定顺序排列.通