线性规划的一个扩展型Gay算法及其复杂性分析.pdfVIP

第1墙 1 工 程 数 学 学 报 v。1．15No_l 。。年 JOURNAlOFENGINEERING MATHEMATICS 附 ·l99a 线性规划的一个扩展型Gay算法 及其复杂性分析。 艾文宝 梁昔明 (西安交通大学理学院，西安 710049) 摘 要 给出了一十求线性规划问题初蛄 目标下界的内点算法，将它和Gay算法[_结合起来．可以去掉Gay算 法中需已知一个初始 目标函数下界的=笨件．并且证明了新算法的选代次数和原算法相比并没有增加．仍然 是 4(n+ 1)L／． 关键词 线性规划，内点算法，Gay算法 ．投影变换 分类号 AMS(1991)90C05- CCLO221 1 引 言 自1947年Dantzig首次提出解线性规划问题的单纯形法以来，它一直是一个非常有效 的算法，但随着线性规划问题的规模越来越大，人们开始关注单纯形法的计算复杂性问题 ， 并证明了单纯形法不是一个多项式时问算法”．1984年Karmarkar发表了线性规划的一个 多项式时H1『内点算法——Karmarkar算法0：．自那以后，有关线性规划的各种内点算法及复 杂性分析的文章如雨后春笋般出现，但归结起来不外乎三类：投影尺度类、仿射尺度类和原 始——对侣类，在投影尺度类算法中最有代表性的算法是Gay算法和Todd算法 ，它们都 需要事先给定一个初始 目标函数下界，这对一般线性规划问题来说是无法做到的．本文在第 2节中给出了一个求线性规划目标函数下界的内点算法，在第3节中将它和Gay算法结合成 一 个新的内点算法，并且证明了新算法也是一个多项式时间算法．尤其重要和有趣的是：虽 然新算法比原Gay算法增加了一个求初始 目标函数下界的算法，但新算法总的迭代次数却 和原Gay算法一样多．仍为4(n+1)L／y． 一 需要指出的是；同样的结合和证明也可以在Todd算法上进行． -本文1997年8月10日收判 墨 艾文宝等一线性规蜘的一十扩展型Gay算法及其复杂性分析 85 2 求 目标函数下界的算法 考虑线性规划问题 (LP) men{cTxl ∈S)， S= {lAx= 6． ≥0) 其中A∈R… 且rank(A)= ．c和 ∈R，b∈R。，最小值为 一cTx’． 对任取的，∈S，，0．作 =diag(}，P=(1，1’．．·．1)∈R一，定义投影变换 丁 ( 一 R一 ) _；m)=(H+l 以 ㈦= 厄 定义 万一[A，一]， ；一 ：『}， = ∈R… 1 = 0，eta-一o) 一 {1 一 0， 一H+ 1，j≥ o) 五= {l，i=n+ 1， ≥ o} {+ 乃)一 {l—e∈力) (LP)问题经过变换r变成下面的问题 (TLP) rain{；i／_J+l∈ ) 并且有 e— T()，cT—I1P垒 广，j‘一T(x)，子空问 上的投影矩阵为 P= f一 (XX)一 一 1／(H+1) ．对任取的实数 ≤广，定义^+= (0，o，…，o，1) ∈R一十l，哼() 一 P(c—ze．+1)垒c，一zeP，()= 脚(2)ll，d。一 lIc，llId一d,e，，d3= lle，f1，则显然有d2= PJ+= +I，d；=P ．+l；e+l，其中 。l和e+分别表示 和e的第n+ 1个分量． 引理 1 若 满足关系式广 一 = (n+ 1)d()，则 ≤ ‘． 证 明 若 (2)=0，则 =广．又因为 (；一fe )( 一P)一 (，’)一 (’一#)一 0 (’一P)；0．所以 ’一尸 一 ；P一广 = 一广一o．得 。=广．故有 = ’． 若rl(z)≠0．定义球体 口(，n+ 1)一 {x…i—ell≤