关于鞍点问题广义预处理HSS—SOR交替分裂迭代方法.pdfVIP

2013年 11月 计 算 数 学 第 35卷第 4期 NOV．，2013 MATHEMATICA NUMERICA SINICA Vo1．35，No．4 关于鞍点问题的广义预处理 HSS—SOR交替 分裂迭代方法 ) 潘春平 (浙江工业职业技术学院人文社科部，浙江绍兴312000) 摘 要 本文研究了鞍点问题的迭代法．在白中治，Golub和潘建瑜提出的预处理对称／反对称分裂 (PHSS)迭代法的基础上，通过结合 GSOR迭代格式，利用两个参数加速，提出了一种广义预处 理 HSS—SOR交替分裂迭代法，并研究了该方法的收敛性．数值结果表明本文所给方法是有效的， 关键词：鞍点问题；交替迭代；PHSS方法；CSOR方法 MR (2000)主题分类：65F10，65F50 1．引 言 考虑鞍点问题 (BAT北)=) 其中A∈Rm 为对称正定 (SPD)矩阵，B∈R 为列满秩矩阵，即rank(B)=n，x，P∈Rm， Y，q∈R (这里我们假定m n)．B 是B的转置矩阵．P，q是已知向量．在此假设下，鞍点问 题 (1)有唯一解 [18]． 鞍点问题是一类对称不定的、系数矩阵对角块中含有奇异矩阵的线性方程组，它出现在 流体力学、二次优化、线性弹性力学、电磁学、图像处理、最小二乘问题等计算科学与工程学 领域 (见 [18—20]及其参考文献)．它的来源和应用相当广泛，比如在不可压缩流体力学中，著 名的Navier．Stokes方程是含约束条件的偏微分方程，这些约束条件代表了一些基本的守恒定 律 (动量守恒定律、质量守恒定律等)，将这些偏微分方程离散化就产生了鞍点问题；电磁学 Maxwell方程 的有限元离散以及二阶椭圆方程 问题的混合有限元方法也会产生鞍点问题．优 化中鞍点问题也被叫做 K“KT方程组”． 目前有许多学者针对求解大型稀疏鞍点问题的迭代法进行了研究和探讨，并构造出了一 些很好的方法，如Uzawa方法、不精确 Uzawa方法及其变形 [6,20,21]、SOR-like方法以及它的 推广和加速 (见 [7，22—24】)、HSS方法、PHSS方法及其推广和加速 (见 [1-5，8，22])以及预处理 Krylov子空间方法等 [9,101．Uzawa方法是一类较早被用于求解鞍点问题的经典迭代法．在此 基础上，Golub等人通过对系数矩阵进行适当的分解，构造出了SOR-like迭代法 2【2】，白中治， Parlett和王增琦用两个参数对其加速，构造了GSOR迭代法 _7J．SOR-like迭代法和 Uzawa 2012年 11月 27日收到． )基金项 目：全国教育信息技术研究课题 (126240641)；浙江省高职研究会课题 (YBl115)；浙江工业职业技 术学院科技计划项 目(1o23401092012014)． 354 计 算 数 学 方法都可以看作 GSOR方法的特例．HSS方法是一类基于矩阵的对称与反对称部分分裂的 两步迭代法 1【]．它于 2003年由白中治，Golub和Ng提出，之后得到了许多学者的密切关注、 深入研究和广泛应用 (见 (2-5，8])．PHSS和 GPHSS方法均是HSS方法在鞍点问题中的具体 应用，它们的优点是收敛性对迭代参数的要求不高，只要参数为正即可． 矩阵多重分裂方法最初 由O’Leafy和 White在研究大规模稀疏线性方程组 Ax：b时 提出 ．国内外许多学者进一步发展了这类方法 [12--17]．对于分裂 A ：M —N =P～ Q(M，P非奇异)，迭代格式 I(+1)M -1Nx(k)+ M-17， 2，… 2() 、