广义鞍点问题的松弛维数分解预条件子.pdfVIP

2012年11月 计算数学 第34卷第4期 NOV．，2012 MATHEMATICANUMERICASINICA V01．34，No．4 广义鞍点问题的松弛维数分解预条件子q) 曹 阳 (南通大学交通学院，江苏南通226019； 苏州大学数学科学学院，江苏苏州215006) 谈为伟 蒋美群 (苏州大学数学科学学院江苏苏州215006) 摘要 dimensional 本文将Benzi等提出的松弛维数分解(Relaxedfactorization，RDF)预条件子进一 步推广到广义鞍点问题上，并称为GRDF(GeneralizedRDF)预条件子．该预条件子可看做是用维数 dimensional 分裂迭代法求解广义鞍点问题而导出的改进维数分裂(Modified split，MDS)预条件子 的松弛形式，它相比MDS预条件子更接近于系数矩阵，因而结合Krylov子空间方法(如GMRES) 有更快的收敛速度．文中分析了GRDF预处理矩阵特征值的一些性质，并用数值算例验证了新预 条件子的有效性． 关键词：广义鞍点问题；维数分裂迭代法；预处理；特征值 MR(2000)主题分类：65F10 1．引 言 我们考虑下面大型稀疏广义鞍点问题： 负“兰(二：署‘)(；)=(乏)三占． c·．·， 其中A∈Rn”是非对称正定矩阵，B∈Rm煳为行满秩的，C∈撼一撕为对称半正定的矩阵，x， 化问题，带等式约束的最小二乘问题，计算流体力学，无网格方法等等，有关这类线性方程组的 来源可参考综述性文献[8】． 预处理迭代法是当前求解大型稀疏线性方程组的研究热点之一，同样是求解广义鞍点问 题(1．1)的主要方法之一．预处理迭代法要预先构造出一个矩阵(称为预条件子)尹，并将(1．1) 等价地转化成与之有相同解的线性系统： 尹一1负“=尹一1占， (1．2) 然后再用迭代法(如一般的分裂迭代法，Krylov子空间迭代法等)进行求解?值得说明的是，将 (1．1)等价地转化为(1．2)的方法称为左预处理方法．此外，还有右预处理方法，分裂预处理方法 等．有关这些预处理方法更详细的论述，可参考[18]．为了求解(1．2)时有更快的收敛速度和较 少的计算复杂度，在构造预条件子尹时，有一些具体的要求，如尹是非奇异的、尹需尽量是贸 +2012年4月23日收到． 11基金项目：苏州大学国家自然科学基金预研基金(SDY2011801) 万方数据 352 计算数学 2012芷 的一个充分逼近、尹的逆不太难求等．预条件子尹的构造方法也是多种多样的，它可以从矩阵 分裂构造迭代法时得到，也可以由矩阵分解得到，还可以根据具体结构通过代数的方法构造得 束预条件子‘3，13】，HSS预条件子[4-7,17]等． 维数分裂(DS)预条件子．该预条件子是由维数分裂迭代法所导出的．在[9】中作者证明了由维 数分裂迭代法求解这类特殊结构的鞍点问题具有无条件收敛性．接着在此基础上，Benzi等人 数矩阵，因而在用迭代法求解时具有更好的收敛性．在这类鞍点问题中，系数矩阵用具有如下 结构： ， f A， o 曰；1 贸=10 A2 磅I，