广义鞍点问题基于PSS的约束预条件子.pdfVIP

2012年5月 计算数学 第34卷第2期 Mav．2012 MATHEMATICANUMERICASINICA V01．34．No．2 广义鞍点问题基于PSS的约束预条件子q) 曹 阳·，z牛 强。 蒋美群， (1．苏州大学数学科学学院，江苏苏州215006) (2．南通大学交通学院，江苏南通226019) (3．西交利物浦大学数理中心，江苏苏州215123) 摘 要 对于(1，1)块为非Hermitian阵的广义鞍点问题，本文给出了一种基于正定和反对称分裂 and definiteskew—Hermitian (Positive splitting，PSS)的约束预条件子．该预条件子的(1，1)块 由求解非Hermitian正定线性方程组时的PSS迭代法所构造得到．文中分析了PSS约束预条件 子的一些性质并证明了预处理迭代法的收敛性．最后用数值算例验证了该预条件子的有效性， 关键词：广义鞍点问题；正定反对称分裂；迭代法；约束预条件子 MR(2000)主题分类：65F10 1．引言 考虑如下大型稀疏广义鞍点问题 脏A；洲小”6 间．形如(1．1)的方程组广泛来源于许多实际问题，如二次优化问题，最小二乘问题，椭圆型偏 微分方程的混合有限元离散，计算流体力学，径向基函数插值，无网格方法等等，详细的讨论 可参考【9]．在大多问题中，A，B，C为大型稀疏阵且广义鞍点矩阵4的条件数很大，因而用预 处理迭代法求解(1．1)较为实际，这也成为当前求解广义鞍点问题的主要方法之一． 众所周知，预处理方法是一种求解大型线性方程组非常有效的方法，构造一个好的预条 件子既能加快迭代速度又能减少每一步迭代的工作量．经过多年的发展，针对鞍点问题已 提出很多种预处理方法，如块对角预条件子【21]，块三角预条件子[5，11，17，”]，约束预处理方 法16，10,12-14,16，2引，HSS预处理方法[3，8·19|，限定的预处理共轭梯度法[4，7】等等，可参考综述 性文献[9]．这些预处理方法中的预条件子可通过代数的方法直接构造得到，也可针对具体问 题构造得到．本文主要考察约束预处理方法并用代数的方法构造了一类基于PSS的约束预 and 条件子，它将[1]中给出的一种基于正定和反对称分裂(Positive—definite +2011年9月5日收到． 1)基金项目： 苏州大学国家自然科学基金预研基金(SDY2011801)，江苏省普通高校研究生创新项目 万方数据 184 计算数学 2012侄 半正定阵，B行满秩，C为Hermitian半正定阵)通过一类基本的迭代格式得到 g仳‘+1=(g—A)乱2+6 (1．2) 其中 ， l G 炉【一B (1．3) 迭代法将G取为 G：=Gskew=二(，+aLs)(，+Q玩)， Ot ，为单位阵：Q为正常数，分析了该预处理方法的收敛性，并验证了该约束预条件子的有效性． and splitting，HSS G：=G日s