关于求解线性互补问题不可行内点算法地若干新地整修过程.docxVIP

膏旦力学廊壬学岔爸卫 性质1．3方阵吖为sufficient阵的充分必要条件是存在某个茁≥0，使得M 为只(x)阵。 注：该性质在参考文献【11]中有详细的证明。 本文将对sufficient型线性互补问题(sufficient-matrix LCP)以及LP型线性 互补问题(LP LCP)、quasi．LP型线性互补问题(quaSi-LP LCP)。下面介绍这几 种线性互补问题。 对线性互补问题(LCP)： fw=Mz+q {wrz：0 1w≥o，z≥o 当系数矩阵M为sufficient阵时，LCP为sufficient型线性互补问题(记为 suffieient．matrix LCP)。 LCP中系数矩阵肘为只(K)阵时，该LCP被称为只(K)型线性互补问题(记 为只(K)LCP)。 求解线性规划问题(LP) rain JJ．瓜=6．x≥0 及其对偶问题(DP) mox b 7y sj． 47y+s=c，s≥0 (其中A∈R”“，C∈R”，b∈R”)的在线性互补规划中等价于求解一个LCP： 』ww，=；M：z。+g， 1w≥o，z≥0 其中M ?)．(?’)=(一≥跏一，=(?]，?Rm+n。此种娓 M=(一：r：] 的线性互补问题被称为线性规划型线性互补问题(记为LP LCP)。显然，LP LCP 复量大学硕士学位论文 (i)若存在f∈，={1，2， ，打}，m。≠0，则在m。处进行转轴，保持变量下标不 变，原问题变为w’；M。z’+g’，其中∽，z；)=z，，Ⅵ)，而在其它／(≠f)处， 0j，zj)=一，，z』)。称这种转轴为一阶主转轴； (ii)若存在i,E1={l，2， ，n}，m。=0，且存在／(≠f)∈，，有mF≠0(由性质 1．2．1知，mp≠0)，则在m口、m』处进行转轴，保持变量下标不变，原问题变为 w+=M’z’+g’，其中(woz；)=(z，，wf)、¨．zj)=(z，，_)，而在其它下标七(≠f，，)处， (w：，z：)：Wk．2k)．称这种转轴为二阶主转轴； 睡： 罚=Mi+日 茹7三：0 茹≥0．三≥0 则系数矩阵廊仍保持sufficient性质。 性质1．5对LP LCP的系数矩阵肘，经一阶或二阶主转轴后仍该LCP为LP LCP。 对非其次线性方程组爿z=b(A c；R椭)，不失一般性，设rd础o)：m，且彳的 前m列线性无关t 于是Ax=b等价于爿口如十AozD=6， 其中 复曼大学硕士学位论文 B={l，2， ，m)，D={m+1，m+2， ，n}。令孑=岛1AD，i=磊1b，于是有下面的定义： 定义1．5对上述非其次线性方程组Ax=6，我们称方程％=AxD+b为 定义1．6对非其次线性方程组的tableau表示形式：w=Mz+q(其中 肌∥“，令州12 一)的某个解向量(爿中的第川个元她，令向量e满 趴，吐旷啪Ⅲ《]=(麓。)闭呲糊称由(州≥。] 的操作为对z，进行向上驱动操作(“ve毛慨ingly)；当九。时，我们称(爿 到[≥]的操作为对z，进行向下驱动操作c“vc毛aecre髂tns时，。 下面介绍线性互补问题的极大互补解与严格互补解，为了叙述的方便我们引 陵i， 圆={i I若存在该LCP的一个解(wJz)，使得wfo，z，：0) Ⅳ_{i 1若存在该LCP的一个解(Ⅵz)，使得wf=o，z，o) 仁{i l对该LCP的任意一个解(wJz)，均满足w『：乃：0} 记f=(圆，M∞，则称f为该LCP的最优剖分。 定义1·8对LCP的一个互补解(w，z)，若z，o§f∈毋，且w，o§，∈Ⅳ 复垦大学硕士学位论文 定义1．9对LCP的一个互补解(w，2)，若BUⅣ兰{1,2， ，月)，即W、z满足： Wi0§zl=0，。』O§W1=0 则称(W，z)为该LCP的一个严格互补解。 性质1．6若LP LCP或quasi·LP LCP有解，则问题一定存在一个严格互补解。 内点算法因初始点不同可分为可行内点算法与不可行内点算法。 现行的许多算法都是要求以一个严格的可行内点作为初始点。所谓“严格的可行内点”，对于线性规划而言，就是满足所有约束条件，尤其是要严格满足不 等式约束：对于线性互补规划而言，就是满足条件w=Mz+吼W0，z0的点。 可行内点算法的初始点和以后的所有迭代点都保持在问题的可行域内。但是，如何找到一个严格的内点作为算法的初始点在实践中并不是一件容易的事。理论上，我们可以通过引入人工变量，把原问题的可行集的维数进行扩充，即对LCP而言，在等式约束的右边增加一条严格大于零的列向量，从而获得一个符合要求的初始点。然而，所增加的列向量很可能破坏了原问题系数矩阵M的某些好性质，这给算法的进行带来了一定的难度。在八十年代末，兴起了另一种内点算法的构造思路：不可行内点算法。 不可行内点算