一种求解LP问题两阶段基点迭代转移方法.docVIP

一种求解LP问题两阶段基点迭代转移方法 　　摘要：利用线性规划的线性、几何平面这一两面性结构特点，定义了LP问题的一种特殊基点转移矩阵及其转移运算，并建立了单纯形基点的定向迭代转移模型，从而提出了一种求解LP问题的两阶段基点定向转移搜索方法.另外，借助新提出的可行域局部ε正则化方法，将退化基点迭代转移转化为非退化基点迭代转移，彻底消除了基点退化对极点转移搜索过程的不利影响. 关键词：线性规划； 基点转移矩阵；退化的；局部正则化；算法 中图分类号：O221.1 文献标识码：A 当今，尽管学术界将LP算法划分成三大类型[1-4]：单纯形类算法、椭球类算法和内点类算法，但它们具有一个共同特征，即均属于迭代算法，根本区别在于迭代方式不同，这也导致了它们在搜索效率上的差异.单纯形类算法的优点在于迭代过程通过 “换基”实现，所以迭代运算均为线性的.缺点是在迭代过程中目标函数值的非严格单调性，当然，出现目标函数值的非严格单调性的内在因素在于退化基点的“一解多基”现象，这一缺陷可能导致迭代过程中的基循环问题.而椭球类算法与内点类算法却恰好相反，它们的优点是在迭代过程中目标函数值的严格单调性，缺点是主要的迭代运算是非线性的，每次迭代的计算量巨大.也正因为这一缺陷，虽然它们是多项式时间的，但实际使用效果并非十分理想，甚至搜索效率还不如单纯形类算法.那么，是否存在迭代运算为线性的，而迭代过程的目标函数值严格单调的LP算法？回答是肯定的，摄动单纯形法[1]就是这样的LP算法. 所谓摄动单纯形法，实为一种先将普通的LP问题转化为无退化现象的LP问题，再利用单纯形法求解的LP算法.从理论上看，摄动单纯形法是一种理想的LP算法，它同时具备上述提到的两个优良特性.然而，在实践中，因为摄动项的添加，相当于要解一个双倍于原LP问题规模的新问题，导致迭代过程的计算量呈爆炸性增长.因此，这种算法的实用性大大降低了. 线性、几何平面特性为线性规划的两大基本结构特性.椭圆类算法、内点类算法都是从非线性规划中移植过来的，属于非线性化方法.自然地，这些算法无法充分反映与利用好线性规划的以上两大结构特性.而传统的单纯形类算法（包括摄动单纯形法）的运算平台是单纯形表，该平台的线性特征十分明显，但几何平面特征明显不足，这也导致了传统单纯形类算法的功能缺陷. 鉴于以上情况，本文欲实现的主要研究目标如下：1）利用线性规划的结构特性，构建更能反映线性规划线性、几何平面这一两面性结构特点的LP新解算模型；2）从分析退化基点的转移机理入手，找到退化基点的转移性缺陷，并制定有效的应对之策；3）在以上两方面研究成果的基础上，提出迭代运算为线性、目标函数值严格单调的LP迭代算法，这也是本文欲达成的最终之研究目标. 1单纯形基点的定向迭代转移模型 至此，通过构建单纯形（包括LP问题的可行域、线性不等式组的解空间）的基点定向转移矩阵，在单纯形的基点与数字矩阵之间建立起了一种对应关系.在矩阵的各功能块中，基解列代表基点，而其方向块的各列代表基点的迭代转移方向（即该超多面体顶点的极方向），强迫性价值系数列代表该基点转移的目标参考方向.本研究的目标之一在于，通过对基点定向转移矩阵的负旋转迭代（等价于基点的迭代转移），最终找到单纯形的一个优化极点. 实际上，依据上述性质1，仅解决了单纯形非退化基点的转移问题.要想利用基点定向转移矩阵的负旋转迭代运算来搜索单纯形的优化极点，还必须解决单纯形退化基点的转移问题.为此，引入下列定义及有关命题： 定义4若单纯形基点定向转移矩阵的基解列中所含零元素的个数大于决策变量数与剩余变量数之和，则称该基点定向转移矩阵对应的基点是退化的，将方向块中非单位行对应的基解列的零元素称为该基点的退化零元素，并称单纯形在该退化基点处具有局部非正则性. 定义5 将“用无穷小量正参数ε替代退化基点的所有退化零元素”的操作称为一次单纯形局部小量正参数ε正则化；若单纯形的极点都是非退化的，则称该单纯形具有正则性（也称LP问题是正则的LP问题）. 性质2只要小量正参数ε足够小，可行域的局部小量正参数ε正则化就不会改变LP问题的最优基. 性质3只要小量正参数ε足够小，经过局部小量正参数ε正则化的基点定向转移矩阵，不论通过多少次负旋转迭代运算，所得矩阵的ε零化矩阵（即，令其中的ε值为零而得到的矩阵）仍然为原LP问题的人工强迫性LP问题可行域基点的转移矩阵. 实际上，之所以要对单纯形进行局部小量正参数ε正则化处理，目的就是让可行域退化的基点非退化化，将退化基点的转移问题转化为非退化极点的转移问题，彻底消除基点退化对迭代搜索的不利影响. 2两阶段基点迭代转移算法 2.1基本思路与构想 先利用单纯形基点的定向