求解双鞍点系统的迭代算法研究.pdfVIP

摘 要 作为具有特殊结构的线性系统, 鞍点问题一直是科学计算领域的研究热点. 事实上, 关于传统鞍点问题的研究已经基本趋向成熟. 本文在这些已有研究的基 础上, 主要探讨了具有特殊双鞍点结构的双鞍点系统的高效求解方法. 首先, 通过对双鞍点系统系数矩阵(双鞍点矩阵) 的合理划分, 分别建立了分 块上三角和分块对角预处理子. 详细分析了预处理后双鞍点矩阵的特征值性质并 给出了更具实际意义的非精确形式. 在数值实验部分, 本文用预处理的广义最小 残量法(PGMRES) 求解双鞍点问题, 同时将数值结果与已有研究中预处理子的表 现作对比, 进而体现这些新提出预处理子的有效性. 接着, 引入非零参数和对称正定子矩阵, 并对原双鞍点矩阵进行分裂, 进而直 接定义了一个基于该分裂的分块下三角预处理矩阵. 同样也分析了预处理后系数 矩阵的特征值分布情况 对于预处理子中的两个对称正定子矩阵 本文给出了几 . , 种不同的选择并记录了这几种情形下该分块下三角预处理子在数值例子中的表 现. 数值结果表明, 与已有研究相比, 该预处理子具有一定的竞争力. 最后, 本文提出新的对称三角分解算法. 该算法可以将原双鞍点矩阵分解为 分块三角矩阵与对称正定矩阵的乘积. 此分解算法的设计思路是利用分块递归 的思想, 不断降低被分解矩阵的维数直至一维. 相比于已有的对称三角分解, 新 算法在分解时间上占有较大的优势. 借助该分解算法, 可以成功地将非对称不定 的双鞍点系统等价转换为对称正定系统, 从而原本无法直接使用的共轭梯度法即 可被用于求解该对称正定系统. 以这样的方式所得到的数值结果相对于直接运用 PGMRES 求解的表现是有可比性的. 关键词: 双鞍点系统; 分块预处理子; 对称三角分解; 特征值分析; 分块递归 i ABSTRACT As a class of linear systems with special structures, saddle point problems have always been a research hotspot in the field of scientific computing. In fact, there are many research on traditional saddle point problems. Based on them, we mainly study the efficient method to solve double saddle point systems with special structure in this paper. First of all, by dividing the coefficient matrix of double saddle point system (i.e. double saddle point matrix) reasonably, we establish a block upper triangular precon- ditioner and a block diagonal preconditioner respectively. We analyze the eigenvalue properties of the preconditioned matrices in detail and give inexact forms for practical application. In numerical experiment part, we use PGMRES to solve the double saddle point problem. At the same time, we compare the numerical results with the perfor- mance of the preconditioners in the existing researches, which show the effect