求解大型稀疏鞍点问题的对称超松弛方法.docVIP

求解大型稀疏鞍点问题的对称超松弛方法 　　摘要:为了快速有效地求解大型稀疏鞍点问题,在SOR-like迭代算法的基础上,结合SSOR分裂,构造了一种解鞍点问题的SSOR-like迭代算法,并研究了该算法的收敛性。数值例子证明:通过参数值的选择,SSOR-like算法比SOR-like算法具有更快的收敛速度和更小的迭代次数,选择了合适的参数值后,可以大大提高算法的收敛效率。 　　关键词:鞍点问题;迭代法;SOR-like 方法;SSOR-like方法 　　中图分类号:O241文献标识码:A文章编号:1009-3044(2010)21-5979-04 　　The SSOR Iterative Algorithms for Large Sparse Saddle Point Problems 　　PAN Chun-ping 　　(Zhejiang Industry Polytechnic College, Shaoxing 312000, China) 　　Abstract: In this paper, we mainly study the iterative algorithms for large sparse saddle point problems(SPP).It is of great interest to develop fast and efficient iterative methods for saddle point problems.Based on the SOR-like iterative scheme and the SSOR splitting, we present the SSOR-like iterative scheme. The convergence results are also given. The numerical results show that our new method is more efficient than SOR-like method. 　　Key words: saddle-point problems; SOR-like method; SSOR-like method; iterative method 　　大型稀疏鞍点问题的求解非常重要,在很多领域都应用过。如电磁学Maxwell方程的有限元离散和计算流体力学中的Stokes方程以及二阶椭圆方程问题的混合有限元方法,含有图像处理问题,限制条件的二次优化问题,最小二乘问题,线性弹性力学问题等等。本文主要考虑具有如下形式的鞍点问题: 　　 (1) 　　其中A∈Rm×m为对称正定矩阵(SPD),B∈Rm×m为列满秩矩阵,即 rank(B)=n ,向量x,p∈Rm,向量y,q∈Rn (这里我们假定m≥n)。BT是B的转置矩阵。P,q是已知向量。在这样的假设下,鞍点问题(2)有唯一解[1]。 　　目前,求解鞍点问题已经存在很多方法,其中包括直接法、Uzawa 算法、SOR-like 算法[10-11]、HSS 算法、GAOR 算法[12] 等。由于直接法在求解大型稀疏线性方程组时,可能会产生许多非零元,从而导致运算量大量增加。同时,对于条件数很大的系数矩阵,直接法往往会产生很大的误差。 因此在求解大型稀疏线性方程组时,我们通常会选择迭代法。考虑(1)的系数矩阵: 　　 (2) 　　这说明鞍点问题的系数矩阵有正和负的特征值。它的另一个特点是对角块中有零矩阵。综上所述,鞍点问题(1)是对称不定的、系数矩阵对角块中含有奇异矩阵的线性系统。它的这些特点使得一些著名的迭代算法,如PCG [3] 算法和SOR算法[4] 难以直接应用。Golub等人通过对系数矩阵进行适当的分解,构造出了一种SOR-like迭代方法[6]。为了表示上的方便,我们将方程组(1)等价表示为: 　　 (3) 　　将方程组(3)的系数矩阵 作如下分裂: 　　 (4) 　　其中: 　　 (5) 　　这里Q∈Rn×n是对称非奇异矩阵。 设ω0为松弛参数,考虑如下迭代格式,得到 　　SOR-like算法: 　　 (6) 　　1)SSOR-like迭代算法: 　　下面仍将方程组(3) 的系数矩阵A 作如下分裂: 　　 (7) 　　其中: 　　 (8) 　　这里Q∈Rn×n是对称非奇异矩阵。考虑如下的SSOR 迭代格式: 　　 (9) 　　 (10) 　　迭代格式(9)可以简记为: 　　 (11) 　　其中迭代矩阵: 　　 (12) 　　其中J=Q-1BTA-1B。由此我们可以得到以下的SSOR-like 迭代算法: 　　设Q∈Rn×n是对称非奇异矩阵, 给定初始向量x(0)∈R