结构矩阵线性互补问题的模系矩阵分裂迭代方法

Linear complementarity problems of structured matrices come from science and engineering computing in many fields, such as finding Nash equilibrium points in convex optimization problem of a rigid body in two, unilateral constrained motion of the electrical network, the ideal diode, widely applied to inequality and optimal control problems in fluid mechanics the free boundary problem. How to effectively solve linear complementarity problems of structured matrices began to become a hot research topic in computational mathematics circles. This project is mainly on the research of modulus-based matrix splitting iterative methods for linear complementarity problems of structured matrices. The contents include: various types of splits based on the system matrix, the establishment of a series of modulus matrix splitting iterative methods and discussing their convergence properties and how to select the parameters of the problem. Doing some perturbation theory analysis of the new iterative method, studying the system coefficient matrix to small perturbations or small errors of the numerical solution will produce what kind of effect, namely the numerical solution will produce what kind of error. At the same time, stability and sensitivity of the new iterative algorithm are also analyzed and measured. This research project aims to promote modulus matrix splitting iterative methods of the linear complementarity problem of structured matrices, and provides the effective methods and theories for solving linear complementarity problems of structured matrices. Therefore, it is of important theoretical and practical significance.

结构矩阵线性互补问题来源于科学与工程计算的许多领域, 如在凸二次优化问题中寻找纳什均衡点,运动的刚体单边约束,不平等的最优控制问题,流体力学中的自由边界问题等等问题中都有广泛的应用。如何有效地求解结构矩阵线性互补问题开始成为计算数学界的一个研究热点。本项目主要研究结构矩阵线性互补问题的模系矩阵分裂迭代方法。内容包括: 研究结构系统矩阵不同类型的分裂；基于这些分裂建立模系矩阵分裂迭代方法；讨论模系矩阵分裂迭代方法的收敛性质以及迭代公式中参数选取的问题。对模系矩阵分裂迭代方法做扰动分析的研究,研究系统系数矩阵发生小的扰动或小的误差对数值解会产生怎样的影响，即数值解会产生怎样的误差。同时，对模系矩阵分裂迭代算法的稳定性和敏感性方面进行分析。本项目旨在促进结构矩阵线性互补问题的模系矩阵分裂迭代方法研究, 为求解结构矩阵线性互补问题提供有效的方法和理论，有着一定的理论和实际意义。

在项目执行期间，我们主要围绕寻找和发展结构矩阵线性互补问题的模系矩阵分裂迭代方法展开了研究。根据实际问题出来的特殊矩阵结构，设计出更为有效的算法，并利用结构矩阵的相关工具，对算法的稳定性、解的敏感性等进行了必要的理论分析和研究，得到了如下研究结果：将一些结构矩阵线性互补问题LCP(q,A)转化为隐式定点方程, 如，H矩阵的线性互补问题，对称矩阵的线性互补问题等，通过对模系矩阵分裂迭代方法引进有效的参数矩阵，从而得到不同的模系矩阵分裂迭代方法。分析模系矩阵分裂迭代方法中迭代矩阵的谱半径，研究模系矩阵分裂迭代方法的收敛性；研究参数的选取使得迭代方法的收敛速度加快，给出一系列求解结构矩阵线性互补问题的模系矩阵分裂迭代方法并给出模系矩阵分裂迭代方法中迭代矩阵的谱分析，从而对迭代算法的收敛性进行分析。研究模系矩阵分裂迭代方法的扰动分析, 分析结构系统矩阵发生小的扰动或小的误差对数值解会产生怎样的影响，即数值解会产生什么样的误差。研究算法的向前(后)误差，对新算法的稳定性和敏感性进行分析。 .总的来说，课题组采用了理论分析和数值实验有机结合的研究方法，取得了比较理想的研究结果，很好的完成了该课题的研究任务。