一种改进的单调增强单纯形算法.pdfVIP

第 28卷第 4期 徐 州 工 程 学 院 学 报 (自然 科 学 版 ) 2013年 12月 Vol_28No．4 JournalofXuzhou InstituteofTechnology (NaturalSciences Edition) Dec．2O13 一 种改进的单调增强单纯形算法 高培旺 (闽江学院，福 州 350121) 摘要 ：考察单调增强单纯形算法 的实际计算性能，并解析其计算效率较低 的原 因．该文提 出一 种改进方法，即从第一阶段算法开始 ，每旋 出一个人工变量，就使非负缩减费用系数 的个数得到单 调增加；在第二阶段算法中，放松对枢轴行 的选择要求，从而可使驱动变量尽快旋入基 中，产生一个 对偶可行解，然后再应用对偶单纯形算法获得 问题的最优解或无可行解的结论．大规模数值试验对 改进算法进行检验 的结果表 明，这种改进算法的计算效率优于经典单纯形算法，单调增强单纯形算 法理论具有实用价值． 关键词 ：线性规划；可行域；单纯形算法；单调增强单纯形算法；计算效率 中图分类号 ：O221．1 文献标志码 ：A 文章编号 ：1674—358X(2013)04—0005—06 单纯形法是求解线性规划实际问题非常有效的算法，至今仍然是人们感兴趣 的研究课题．由于单纯形法 在枢轴主元的选择上相 当灵活 ，因而产生 了许多变式，如各种定价准则_l1。]，原始一对偶方法 ]，单调增强算 法_7]．其中，单调增强算法 (简称 MBU单纯形算法)是由Anstreicher和 Terlaky开发的，具有独特的设计思 想 ，亦即保持非负缩减费用系数在迭代过程中仍然非负，并使非负缩减费用系数 的个数单调增加 ，直至所有 缩减费用系数成为非负数 ，此时，问题的最优解也随之产生．有趣的是 ，MBU单纯形算法也许会迭代到可行 域之外的极点上，但随着驱动变量进入基中，算法又将迭代拉 回到可行域的极点． 根据 MBU单纯形算法的思想 ，Bilen等 提 出了处理退化可行性 问题的有 限枢轴准则 ；本人结合其它 方法讨论具有不可行基的一般线性规划问题的求解_9 ．遗憾 的是，Anstreicher等在文献[7]中并没有报告 MBU单纯形算法的数值试验结果 ，致使该算法理论没有得到充分重视和发掘． 本研究对 MBU单纯形算法的实际计算性能进行考察 ，发现其计算效率不如经典单纯形算法，并通过迭 代计算过程 ，找出其计算效率较低 的原 因．为此 ，提出一种改进方法 ，从第一阶段算法开始 ，每旋出一个人工 变量 ，在保持非负缩减费用系数在迭代过程 中仍然非负的前提下，就使非负缩减费用系数的个数得到单调增 加 ；在第二阶段算法 中，放松对枢轴行的选择要求 ，这样可使驱动变量尽快旋人基中，产生一个对偶可行解 ， 然后再应用对偶单纯形算法获得 问题 的最优解或无可行解 的结论．最后，从线性规划标准测试库 NETLIB[1和整数规划标准测试库 MIPLIB_1中选取了22个典型算例 ，通过 MATLAB编程在计算机上实 现大规模数值试验．结果表明这种改进算法的计算效率优于经典单纯形算法 ，因而，单调增强单纯形算法理 论是有实用价值的． 1 MBU单纯形算法原理及计算性能 考虑如下形式的线性规划问题 ： (LP) max 厂一CX， S．t． A =b， ≥0， 这里 ，A∈R ( )，C、b分别是相应维数的行向量和列向量 ，且 6≥0． 不妨设 B与 』v分别是A 的基与非基子矩阵，C、CN分别是基与非基相应的费用系数向量 ，b是右手边 向 收稿 日期 ：2013—08—25 基金项 目：广西 自然科学基金项 目(0728260)；国家星火计划项 目(2013GA690426)；闽江学院人才引进基金项 目(MJU2012001) 作者简介 ：高培 旺(1964一)，男，湖南宁远人 ，教授 ，博士，硕士生导师 ，主要从事最优化理论及其应用研 究．