线性互补问题的并行多分裂松弛迭代算法.docVIP

线性互补问题的并行多分裂松弛迭代算法 段班祥 李郴良 徐安农 （桂林电子工业学院 计算科学与应用物理系 广西桂林 541004） 摘要：运用矩阵多重分裂理论，同时考虑并行计算与松弛迭代法，得到一类求解线性互补问题的高效数值算法。当问题的系数矩阵为对角元为正的H-矩阵或对称半正定矩阵时，证明了算法的全局收敛性。该算法与已有算法相比，具有计算量小、计算速度快等特点，因而特别适于求解大规模问题。数值试验的结果说明了算法的有效性。 关键词：线性互补问题；多重分裂；松弛迭代 中图分类号：O241.6 AMS(1991)主题分类: 90C33 文献标识码： A 1 引言 互补问题是一类重要的优化问题，在线性规划、凸二次规划、双矩阵对策及经济与交通平衡等科学与工程中有着广泛的应用，因此研究互补问题的高效数值算法是非常必要的。矩阵的多重分裂方法最初是O’Leary 和 White在研究线性方程组时提出的概念，Frommer A. 和 Mayer G. 发展了这类算法，提出了多分裂松弛并行迭代算法。这类方法将线性问题的系数矩阵进行不同的分裂，在并行计算中每个处理器对相应的某种分裂所得到的线性子问题进行计算，再将每个结果按适当的方式组合构成整个问题的迭代解。这种方法的优点在于凭借适当的组合方式使得并行计算中相应的某些分量不必计算，从而大大提高计算效率。 本文考虑如下线性互补问题：求使得 ， ， （1.1）其中，.本文将[2]中求解线性方程组的多分裂松弛迭代法进行推广，得到一类求解线性互补问题的快速算法。在第二节我们给出必要的预备知识；第三节给出松弛迭代算法；第四节当为H-矩阵类时对算法进行收敛性分析；第五节当为对称矩阵类时对算法进行收敛性分析；最后一节给出算法的数值实例，以说明算法的有效性。 预备知识 定义2.1 如果，有，称矩阵为Z-矩阵；如果是非奇异的Z-矩阵，且，称为非奇异M-矩阵。设，其中 作者简介：段班祥（1971-），男，汉族，湖南武冈人，桂林电子工业学院计算科学与应用物理系硕士研究生，主要从事数值计算、应用软件开发等方面的研究。 联系地址：桂林电子工业学院计算科学与应用物理系，邮编：541004. 联系电话：0773-5601272. E-mail: duanbanx@ ，，称为的比较矩阵；如果的比较矩阵 是非奇异的M-矩阵，则称为非奇异的H-矩阵，简称为H-矩阵。 定义2.2 设，如果①是非奇异的且，②且，称为矩阵的正则分裂；如果且，则称为矩阵的弱正则分裂；若对，线性互补问题（1.1）总有解，则称为Q-矩阵；一个分裂称为Q-分裂，是指为Q-矩阵。 定义2.3 设为非奇异的实矩阵，，，满足下列条件：①，，②非奇异，，③（单位阵），则称三元素集，为矩阵的多重分裂。 引理2.4 设满足，那么. 引理2.5 设是H-矩阵，，，那么（1）非奇异；（2）；（3）是非奇异的且. 引理2.6 如果为对角元为正的H-矩阵，则对任意给定的，线性互补问题（1.1）有唯一解。 算法的建立 基于第二节定义2.3中的矩阵的多重分裂定义，我们首先给出求解问题（1.1）的多分裂松弛迭代算法。 算法3.1 （多分裂松弛迭代算法） （1）任意给定初始值，置； （2）对的多重分裂：（），并行求解： （3.1） （3）（为松弛因子）； （4），转第（2）步，直之收敛。 注：将划分为一些子集的并：， ，其中各子集可以是重叠的，也可以是非重叠的；定义，如下：，且 .因此，在计算中，只须计算对应于的分量，这样以来，计算量就减小了许多，又由于本算法适于并行计算，计算效率自然就提高了。 为H-矩阵类时收敛性分析 本节我们给出当为H-矩阵时算法3.1的收敛性定理。 引理4.1 设（）为的多重分裂，、分别为对角元为正的H-矩阵，为问题（1.1）的唯一解，则对任意的向量及任意的初始向量，由算法3.1产生的序列及满足 . （4.1） 证明：由于为对角元为正的H-矩阵，由引理2.6知：子问题（3.1）有唯一解，即是唯一定义的。对不等式（4.1），我们逐个分量的来加以证明。 对任意的，假设.在上述假设下，如果，则由于（4.1）左边向量的第个分量非正而右边向量的分量始终为非负的，（4.1）式显然成立。下面假设，则由算法3.1可得.另一方面我们有，两式相减可得：，从而 . （因为的对角元为正，）. 类似可证当时（4.1）也成立。 这就证明了引理4.1. 下面给出算法3.1的收敛性定理： 设，，，若非奇异，则令.为方便起见令，若为H-矩阵，则由引理