运筹学 割平面法.pptVIP

Cj 1 1 0 0 0 CB XB b x1 x2 x3 x4 s1 1 x1 5/3 1 0 5/6 －1/6 0 1 x2 8/3 0 1 －2/3 1/3 0 0 s1 －2/3 0 0 －1/3 －1/3 1 σj －13/3 0 0 －1/6 －1/6 0 Cj 1 1 0 0 0 CB XB b x1 x2 x3 x4 s1 1 x1 5/3 1 0 5/6 －1/6 0 1 x2 8/3 0 1 －2/3 1/3 0 0 s1 －2/3 0 0 －1/3 －1/3 1 σj －13/3 0 0 －1/6 －1/6 0 Cj 1 1 0 0 0 CB XB b x1 x2 x3 x4 s1 1 x1 0 1 0 0 －1 0 1 x2 4 0 1 0 1 －2 0 x3 2 0 0 1 1 －3 σj －4 0 0 0 0 －1/2 ﹡ Cj 1 1 0 0 0 CB XB b x1 x2 x3 x4 s1 1 x1 0 1 0 0 －1 0 1 x2 4 0 1 0 1 －2 0 x3 2 0 0 1 1 －3 σj －4 0 0 0 0 －1/2 ﹡ 得到整数最优解，即为整数规划的最优解，而且此整数规划有两个最优解： X＊= (0, 4), Z = 4, 或 X＊= (2, 2), Z = 4。 CB XB b x1 x2 x3 x4 x5 0 x3 4 0 0 1 -1/3 4/3 4 x2 4/3 0 1 0 2/9 -5/9 11 x1 8/3 1 0 0 1/9 2/9 σj 0 0 0 -19/9 -2/9 CB XB b x1 x2 x3 x4 x5 s1 0 x3 0 0 0 1 -1 0 6 4 x2 3 0 1 0 1/2 0 -5/2 11 x1 2 1 0 0 0 0 1 0 x5 3 0 0 0 1/2 1 -9/2 σj 0 0 0 -2 0 -1 （2 ，3） （一）、计算步骤： 1、用单纯形法求解( IP )对应的松弛问题( LP )： ⑴.若( LP )没有可行解，则( IP )也没有可行解，停止计算。 ⑵.若( LP )有最优解，并符合( IP )的整数条件，则( LP )的最优解即为( IP )的最优解，停止计算。 ⑶.若( LP )有最优解，但不符合( IP )的整数条件，转入下一步。 第二节 割平面法 2、从(LP)的最优解中，任选一个不为整数的分量xr,,将最优单纯形表中该行的系数 和 分解为整数部分和小数部分之和，并以该行为源行，按下式作割平面方程： 3、将所得的割平面方程作为一个新的约束条件置于最优单纯形表中（同时增加一个单位列向量），用对偶单纯形法求出新的最优解，返回1。 的小数部分 的小数部分 例一：用割平面法求解整数规划问题 解：增加松弛变量x3和x4 ，得到(LP)的初始单纯形表和最优单纯形表： Cj 0 1 0 0 CB XB b x1 x2 x3 x4 0 x3 6 3 2 1 0 0 x4 0 -3 2 0 1 σj 0 0 1 0 0 Cj 0 1 0 0 CB XB b x1 x2 x3 x4 0 x1 1 1 0 1/6 -1/6 1 x2 3/2 0 1 1/4 1/4 σj -3/2 0 0 -1/4 -1/4 此题的最优解为：X＊ (1 , 3/2) Z = 3/2 但不是整数最优解，引入割平面。以x2 为源行生成割平面，由于 1/4=0+1/4, 3/2=1+1/2, 我们已将所需要的数分解为整数和分数，所以，生成割平面的条件为: Cj 0 1 0 0 CB XB b x1 x2 x3 x4 0 x1 1 1 0 1/6 -1/6 1 x2 3/2 0 1 1/4 1/4 σj -3/2 0 0 -1/4 -1/4 现将生成的割平面条件加入松弛变量，然后加到表中： Cj 0 1 0 0 0 CB XB b x1 x2 x3 x4 s1 0 x1 1 1 0 1/6 -1/6 0 1 x2 3/2 0 1 1/4 1/4 0 0 s1 -1/2 0 0 -1/4 -1/4 1 σj -3/2 0 0 -1/4 -1/4 0 CB XB b x1 x2 x3 x4 s1 0 x1 2/3 1 0 0 -1/3 2/3 1 x2 1 0 1 0 0 1 0 x3 2 0 0 1 1 -4 σj -1 0 0 0 0 -1 此时，X1 ＝(2/3, 1), Z=1,仍不是整数解。继续以x1为源行生成割平面，其条件为： CB XB b x1 x2 x3 x4 s1 0 x1 2/3 1 0 0 -1/3 2/