鞍点问题的PSS预处理技术探索.doc

鞍点问题的PSS预处理技术 学生姓名：王小二 学 号：1209**** 专业班级：信息与计算科学12-*班 2013年7月8日 线性表的设计和实现PSS预处理技术，并首次研究了在参数很小的情PSS预处理矩阵特征值的分布情况，用Navier-Stokes方程模型验证了所得理论结果的有效性。 关键词：鞍点问题；特征值估计；预处理子；埃尔米特和反埃尔米特分裂 ABSTRACT In this paper, we provide a concise overview of some preconditioning techniques for saddle point problems, and focus on PSS preconditioning technique. When parameter is sufficiently small, we firstly analysis the eigenvalue distribution of the preconditioned saddle point matrices. Navier-Stokes equation model is used to illustrate the effectiveness of the presented theoretical results. Keywords: saddle point problem; eigenvalue estimate; preconditioner; Hermitian and skew-Hermitian splitting 此处插入目录 第1章 前 言 1.1 鞍点问题背景 鞍点问题经常出现在许多科学和工程的应用中，如约束最优化、约束与加权最小二乘问题、经济学、电磁学、图像重构、图像定位和识别等（文献[2,10,12,15,18,21]）。事实上大多数带约束条件的问题都可产生鞍点系统。 许多应用问题可以用最小化约束问题阐述。通常一些问题是无限维的且非线性的， 通过离散化，可以产生大规模有限维问题。这些问题通常被替换成一系列拥有线性等价约束条件的二次最小化约束问题： （1-1） （1-2） 其中对称半正定， 和 是给定的向量，一阶最优化条件是通过下面的线性系统给出的： （1-3） 其中，是拉格朗日乘子矢量。（1-3）这种类型的线性系统被称为鞍点问题，因为其解 是下列拉格朗日函数的一个鞍点： （1-4） 拥有鞍点的线性系统也出现在一些离散的物理模型里，例如力学结构和RCL电路等。 更一般的，认为拥有下列形式的线性系统是鞍点问题，其中和跟之前的形式一样，为对称半正定矩阵。 （1-5） 此系统产生于稳定的混合有限元法和描述微压缩流体，固体学的方程的离散化，加权最小二乘问题和优化中的内点法等，并且相对于其它块来说有较小的模。这里，是大型稀疏矩阵，因此一般情况下（1-5）要用迭代法解（文献[5,6,22,27,30,31]），但不幸的是对鞍点系统直接使用Krylov法时收敛速度很慢，这就需要对（1-5）进行预处理再使用Krylov子空间法以获取更快的收敛速度。在过去的几年里，为了找到有效地预处理子，国内外学者做了很多的工作（文献[1,33,14,15,17]）。本文将介绍一些比较常用的预处理子，并对PSS预处理子进行谱分析。 1.2 鞍点矩阵的性质 本节介绍鞍点问题系数矩阵（鞍点矩阵）的有关性质，请见文[9]。 称下面的高阶，稀疏和非奇异的线性系统为广义鞍点问题： （1-6） 其中半正定，即对称半正定，满秩，对称半正定，且。 改变（6）中第二行块的符号可以产生一个等价的鞍点问题： （1-7） 记为的特征值集合，不管是否对称，下面定理给出了不对称鞍点矩阵的一些特定性质： 定理1.1 不对称鞍点矩阵都有一些特定的性质： 1.非奇异。 2.是半正定的，即对任意的有。 3.是正半稳定的，也就是说的特征值的实部都是非负的，即对于所有，。 4.如果是正定的，则是正稳定的：对于所有，。 【证】证明第1个：构造，这里，则可以得到 （1-8） 从还可得到 因为加号前后的两部分都是非负的，所以必定有 且 但是 因为是对称半正定的，所以（看文献[19]第400页）。类似的，，由对称半正定可得，因此（用（1-8）的第二个式子）。由可得。但是如果从（1-8）的第一个式子我们可以得到，由满秩可得。因此唯一的解是零解，所以非奇异。 证明第2个：对于任意的我们有，这里 是的对称部分。很明显是半正定的。 证明