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

高校应用数学学报 2012，27(4)：456—464 关于鞍点问题的预处理HSS—SOR 交替分裂迭代方法 潘春平 (浙江工业职业技术学院人文社科部，浙江绍兴 312000) 摘 要：为了高效地求解大型稀疏鞍点问题，在白中治，Golub和潘建瑜提出的预处理 对称／反对称分裂(PHss)迭代法的基础上，通过结合SOR—like迭代格式对原有迭代算 法进行加速，提出了一种预处t~HSS—SOR交替分裂迭代方法，并研究了该算法的收敛 性．数值例子表明：通过参数值的选择，新算法~hSOR—like和PHSS算法都具有更快的 收敛速度和更少的迭代次数，选择了合适的参数值后，可以提高算法的收敛效率． 关键词：鞍点问题；交替迭代；PHSS方法：SOR—like方法 中图分类号：0241．6 文献标识码：A 文章编号：1000—4424(2012)04—0456—009 §1 引 言 考虑具有如下形式的鞍点问题： ( )=) (1) 其中A∈Rm 为对称正定(SPD)矩阵，B∈R 为列满秩矩阵，~1]rank(B)=咒，向量 ，P∈ R ，向量Y，q∈ (这里假定m 几)．B ’是B的转置矩阵．P，g是已知向量．在这样的假设下， 鞍点问题f11有唯一解l1J_ 鞍点问题是一类对称不定的、系数矩阵对角块中含有奇异矩阵的线性系统．它出现在流 体力学、二次优化、线性弹性力学、电磁学、图像处理、最小二乘问题等计算科学与工程学 领域 (见 1『—31及其中的参考文献)．它的来源和应用相当广泛．比如在不可压缩流体力学中，著名 的Navier—Stokes方程是含约束条件的偏微分方程，这些约束条件代表了一些基本的守恒定律f动 量守恒定律、质量守恒定律等)，将这些偏微分方程离散化就产生了鞍点问题．电磁学Maxwell方 程的有限元离散以及二阶椭圆方程问题的混合有限元方法也会产生鞍点问题．在优化系统中，鞍 点问题也被叫做 K“KTSystem”． 收稿 日期：2012—03—26 修回日期：201209—29 基金项目：全国教育信息技术研究课题(126240641)；浙江省高职研究会课题(YBl115)；浙江工业职业技术学院科 技计划项 目(1023401092012014) 潘春平：关于鞍点问题的预处理HSS—SOR交替分裂迭代方法 457 由于这类问题经过某些方法离散后的系数矩阵通常是大型稀疏的，直接求解方法往往需要 大量的内存和计算量而不切实际，而迭代法每次迭代所需的运算量相对较少，所 以在研究大型 稀疏鞍点问题的数值解时，迭代法得到了更广泛的关注，并逐渐成为近年来数值代数方向的研 究热点之一，然而它的特点使得一些著名的迭代算法，~nPcG算法和s0R算法难以直接应用．因 此，目前有许多学者针对求解大型稀疏鞍点问题的迭代法进行了研究和探讨，并构造出了一些 较好的算法．~Uzawa方法、不精确Uzawa方法及其变形f3—51、S0R—like算法以及它的推广和 加速算法(见文献 1o])、HSS算法、PHSS算法及其推广和加速算法(见文献[11—17])、Krylov子 空间方法等．其中Uzawa算法是一类较早被提出用于求解鞍点问题的经典迭代算法．在此基础 上Golub等人通过对系数矩阵进行适当的分解，构造出了--~S0R—like迭代方法f61，Uzawa算法 可以看作是它的一个特例，HSS方法是一类基于矩阵的对称与反对称分裂的两步迭代方法[111． 它于2003年由白中治，GeneH．Golub$~MichaelK．Ng提出、之后受到了许多学者的深切关注，深 入研究和广泛应用． 本文在回顾SOR—like方法及预处理对称与反对称分裂迭代法的基础上，通过对两种方法的 综合利用，提出了一种新的预处理HSS—S0R交替分裂 (PHSS—SOR)迭代方法，并对新算法的收敛性 进行了研究，大量的数值实验表明，通过选择合适的参数，可以提高算法的收敛效率． §2预处NH