关于线性互补问题的模系矩阵分裂迭代方法.pdfVIP

2012年11月 计算数学 第34卷第4期 NOV．，2012 MATHEMATICANUMERICASINICA V01．34，No．4 关于线性互补问题的模系矩阵分裂迭代方法术 张丽丽 (河南财经政法大学数学与信息科学系，郑州450002) 摘 要 模系矩阵分裂迭代方法是求解大型稀疏线性互补问题的有效方法之一．本文的目标是归纳总 结模系矩阵分裂迭代方法的最新发展和已有成果，主要内容包括相应的多分裂迭代方法，二级多 分裂迭代方法和两步多分裂迭代方法，以及这些方法的收敛理论． 关键词：线性互补问题；模系方法；矩阵多分裂；二级多分裂 MR(2000)主题分类：65F10，65K05，68W10，90C33 1．引 言 E E 对于给定的矩阵A∈瓞似”和向量gR”，线性互补问题就是要找到一对向量Z，rR“， 使得 71：=Az+q之0，z芝0，ZT(Az+q)=0(1．1) 成立．线性互补问题(1．1)常简记为LCP(q，A)． 线性互补问题在经济和工程中有着广泛的应用．早期的线性互补问题主要来源于线性规 划，二次规划和双矩阵对策的研究；详见【12，21，23]及其参考文献．随着线性互补问题在基本 理论和数值算法上的不断发展和日益完善，其实际应用也逐渐扩展到经济，金融，控制，交通 和力学等诸多领域，例如市场供需平衡，期权定价，交通网络平衡，弹性接触，自由边界问题等； 详见[12：14]- 线性互补问题的数值解法主要有直接方法和迭代方法两大类．直接方法常用于求解中小 型线性互补问题，具有求解速度快且解的精度高等优点。但是对于工程应用中出现的大型稀 疏线性互补问题，直接方法最大的缺点在于破坏矩阵稀疏结构，其计算复杂度和计算内存随着 问题规模的增加而剧烈增加，甚至可能会因舍入误差的影响或者内存的限制而导致求解过程 中断．在这种情况下，迭代方法相对于直接方法有着无可比拟的优势，主要体现在求解过程中 能够保持稀疏结构，数据存储量小，对舍入误差不敏感，求解速度快等方面．本文将主要介绍 模系矩阵分裂迭代方法，这是最近提出的用于求解线性互补问题的一种迭代方法，在实际应用 中非常有效且易于实现15J， 给定方阵A=(aij)∈R似”，其比较矩阵(A)=((aij))E酞“”定义为 ，、 f 10巧1，i=J， 《o巧)2 1 ．，． L—Iaijl，o尹J， 如果方阵A对于所有的i≠J都有aij≤0，则称其为Z一矩阵．给定非奇异方阵A，如果 +2012年7月3日收到 万方数据 374 计算数学 为日一矩阵[1l,2s]．特别地，如果一个H一矩阵的对角元素都是正的，则记为三0一矩阵【3】．设 向量q，LCP(q，A)都存在惟一解M引．在本文中我们主要考虑A为且}一矩阵时的收敛理论， 这时LCP(q，A)存在惟一解，不妨记为孙 2．模迭代方法 求解线性互补问题的模迭代方法是基于性质：对于任意的向量z∈腿“，…+z和蚓一。都 是非负向量且各分量之间具有互补性质．van 和r：=IX卜。，将线性互补问题(1．1)等价地转化为一个不动点方程组 (j+A)x=(I—A)lxl—q (2．1) 并由此构造了求解线性互补问题的模迭代方法： 方法2．1．(模迭代方法) (1)选取初始向量z(o)∈R”，置m：=0； (2)对于已得到的z(…，求解线性方程组 (J+A)z(“+1)=(I—