运筹学

陈建华

目录

  • 1 第一章    绪论
    • 1.1 第一节 运筹学的定义与发展简史
    • 1.2 第二节 运筹学的基本特点和工作步骤
    • 1.3 第三节 运筹学的主要分支
    • 1.4 第四节 运筹学的应用
  • 2 第二章  线性规划
    • 2.1 第一节 线性规划概述
    • 2.2 第二节 线性规划问题及其数学模型
    • 2.3 第三节 线性规划图解法及其几何意义
    • 2.4 第四节 线性规划单纯形法与单纯形表
    • 2.5 第五节 单纯形法的矩阵描述
    • 2.6 第六节 人造基下的单纯形法
    • 2.7 第七节 线性规划典型例题及应用
  • 3 第三章 运输问题
    • 3.1 第一节 运输问题的数学模型及其特征
    • 3.2 第二节 运输模型的求解---表上作业法
    • 3.3 第三节 运输问题的推广
  • 4 第四章 整数规划
    • 4.1 第一节 整数规划概念与特点
    • 4.2 第二节 分枝定界法
    • 4.3 第三节 割平面法
    • 4.4 第四节 0—1规划与隐枚举法
    • 4.5 第五节 指派问题与匈牙利法
    • 4.6 第六节 典型例题及应用
  • 5 第五章 图与网络
    • 5.1 第一节 图的基本概念
    • 5.2 第二节 树
    • 5.3 第三节 最短路问题
    • 5.4 第四节 网络最大流问题
    • 5.5 第五节 Euler图
    • 5.6 第六节 中国邮递员问题
  • 6 第六章 网络计划
    • 6.1 第一节 网络计划图
    • 6.2 第二节 网络计划图的时间参数
    • 6.3 第三节 网络计划的优化
  • 7 第七章 排队论
    • 7.1 第一节 排队论的基本概念
    • 7.2 第二节 排队系统常用分布
    • 7.3 第三节 单服务台模型
  • 8 第八章 存储论
    • 8.1 第一节 存储论基础
    • 8.2 第二节 确定性库存模型
    • 8.3 第三节 确定性库存模型的参数分析
    • 8.4 第四节 随机型存储模型
  • 9 第九章 决策论
    • 9.1 第一节 决策论基本问题
    • 9.2 第二节 完全不确定型决策
    • 9.3 第三节 风险型决策
    • 9.4 第四节 效用理论在决策中的应用
第三节 割平面法

第三节割平面法

割平面法是1958年由R.E.Gomory提出来的,因此又称为Gomory的割平面法。其基本思路是,先不考虑变量是整数的约束条件,解整数规划相应的松驰问题,所得的最优解如不符合整数条件,再增加新线性约束条件(其几何术语称为割平面),使得从原可行域中切割掉一部分,这部分只包含非整数解,但没有切割掉任何整数可行解。割平面法就是指出怎样找到适当的割平面(不见得一次就找到),使切割后最终得到的可行域中的一个整数坐标的极点,它恰好是问题的最优解。由此可见,割平面法的基础仍然是用解线性规划的方法去求解整数规划问题。

割平面法在求解纯整数规划问题的基本步骤:

(1)解松驰问题。先不考虑整数条件,求解整数规划(IP)对应的松弛问题(LP),可能得到以下情况之一:

(LP)没有可行解,则(IP)也没有可行解,停止计算。

(LP)有最优解,并符合(IP)的整数条件,则(LP)的最优解即为(IP)的最优解,停止计算。

(LP)有最优解,但不符合(IP)的整数条件,转入下一步。

(2)求切割方程。按以下方法求切割方程:

为松驰问题(LP)最优解中分为数值的一个基变量,由最优单纯形表可得

                                        (3-3)

都分解成整数部分N和非负真分数f之和,即:

           (3-4)

(3-3)代入(3-4)可得

                        (3-5)

当变量约束为整数时,(3-5)式左边看必须是整数,但右边因为有,所以不能为正,即

                                        (3-6)

(3-6)就是一个切割约束,引入松驰变量,得切割方程如下:

(3)解新的松驰问题。在原来松驰问题(LP)约束条件的基础上增加一个切割约束,构成新的线性规划问题(LP1),再求出(LP1)的最优解。再返回到步骤(2)

(4)重复以上步骤,直到求得最优解或判断原问题无可行解为止。

现举例说明割平面法在求解纯整数规划问题中的应用。

3.2】用割平面法求解下列整数规划。

解:用单纯形法求解该整数规划对应的松弛问题(LP),得其最优单纯形表3-1

3-1 最优单纯形表

                                                             

 

 
 

1

 
 

1

 
 

0

 
 

0

 
 

 
 

 
 

 
 

 
 

 
 

 
 

 
 

2

 
 

 
 

5/2

 
 

0

 
 

1

 
 

1/2

 
 

-1/2

 
 

3

 
 

 
 

13/4

 
 

1

 
 

0

 
 

-1/4

 
 

3/4

 
 

 
 

0

 
 

0

 
 

-1/4

 
 

-5/4

 

其中,为松驰变量。(LP)的最优解为,显然不符合整数条件。

从最优单纯形表中两个约束条件中提出一个如下:

将上式中系数和常数项分解成整数和非负真分数之和,得以下式子:

移项使整数部分在等号的左边,小数部分在等号的右边得:

都取非负整数时,由于上式的左边取整数,与之相等的右边也只能取整数,而右边不可能取得正整数,所以只能取负整数,即有不等式:

即:

引入松驰变量,得割平面方程:

                  (3-7)

将上式作为一个新的约束条件,增加到(LP),构成(LP1)。按如下方法求解(LP1):将(3-7)式添加在(LP)的最优单纯形表中,得到表3-2

3-2 割平面法计算过程

                                                                                       

 

 
 

1

 
 

1

 
 

0

 
 

0

 
 

0

 
 

 
 

 
 

 
 

 
 

 
 

 
 

 
 

 
 

2

 
 

 
 

5/2

 
 

0

 
 

1

 
 

1/2

 
 

-1/2

 
 

0

 
 

3

 
 

 
 

13/4

 
 

1

 
 

0

 
 

-1/4

 
 

3/4

 
 

0

 
 

0

 
 

 
 

-1/2

 
 

0

 
 

0

 
 

-1/2

 
 

-1/2

 
 

1

 
 

 
 

0

 
 

0

 
 

-1/4

 
 

-5/4

 
 

0

 

3-2中,检验数全部非负,但b列数存在负分量,用对偶单纯形法继续迭代,得最优单纯形表3-3

3-3割平面法计算过程

                                                                                       

 

 
 

1

 
 

1

 
 

0

 
 

0

 
 

0

 
 

 
 

 
 

 
 

 
 

 
 

 
 

 
 

 
 

2

 
 

 
 

2

 
 

0

 
 

1

 
 

0

 
 

-1

 
 

1

 
 

3

 
 

 
 

7/2

 
 

1

 
 

0

 
 

0

 
 

1

 
 

-1/2

 
 

0

 
 

 
 

1

 
 

0

 
 

0

 
 

1

 
 

1

 
 

-2

 
 

 
 

0

 
 

0

 
 

0

 
 

-1

 
 

-1/2

 

最优解为,不符合整数条件,需要再求一次切割方程。

从表3-3中提出第二个约束条件:

按照前述方法,求得切割约束为:

引入松驰变量,得割平面方程:                   (3-8)

(3-8)式添加在(LP1)的最优单纯形表3-3中,得到表3-4

3-4割平面法计算过程

                                                                                                                     

 

 
 

1

 
 

1

 
 

0

 
 

0

 
 

0

 
 

0

 
 

 
 

 
 

 
 

 
 

 
 

 
 

 
 

 
 

 
 

2

 
 

 
 

2

 
 

0

 
 

1

 
 

0

 
 

-1

 
 

1

 
 

0

 
 

3

 
 

 
 

7/2

 
 

1

 
 

0

 
 

0

 
 

1

 
 

-1/2

 
 

0

 
 

0

 
 

 
 

1

 
 

0

 
 

0

 
 

1

 
 

1

 
 

-2

 
 

0

 
 

0

 
 

 
 

-1/2

 
 

0

 
 

0

 
 

0

 
 

0

 
 

-1/2

 
 

1

 
 

 
 

0

 
 

0

 
 

0

 
 

-1

 
 

-1/2

 
 

0

 

用对偶单纯形法继续迭代,得最优单纯形表3-5

3-5割平面法计算过程

                                                                                                                     

 

 
 

1

 
 

1

 
 

0

 
 

0

 
 

0

 
 

0

 
 

 
 

 
 

 
 

 
 

 
 

 
 

 
 

 
 

 
 

2

 
 

 
 

1

 
 

0

 
 

1

 
 

0

 
 

-1

 
 

0

 
 

2

 
 

3

 
 

 
 

4

 
 

1

 
 

0

 
 

0

 
 

1

 
 

0

 
 

-1

 
 

0

 
 

 
 

3

 
 

0

 
 

0

 
 

1

 
 

1

 
 

0

 
 

-4

 
 

0

 
 

 
 

1

 
 

0

 
 

0

 
 

0

 
 

0

 
 

1

 
 

-2

 
 

 
 

0

 
 

0

 
 

0

 
 

-1

 
 

0

 
 

-1

 

最优解为,符合整数条件,对应的目标函数

所以原问题的最优解,对应的目标函数