第章整数规划割平面法.docx

割平面法 求解整数规划问题： Max Z=3xi+2x2 2xi+3x2 14 4xi+2x2 18 xi,x 2 0,且为整数 解：首先，将原问题的数学模型标准化，这里标准化 有两层含义：（1）将不等式转化为等式约束，（2）将 整数规划中所有非整数系数全部转化为整数，以便于 构造切割平面。从而有： Max Z=3xi+2x2 2xi+3x2+X3=14 2xi+X2+X4=9 x1,x 2 0,且为整数 利用单纯形法求解，得到最优单纯形表，见表 1: 表1 CB X b 3 2 0 0 X 茨 X3 X4 2 X 5/2 0 1 1/2 -1/2 3 X 13/4 1 0 -1/4 3/4 j 59/4 0 0 1/4 5/4 最优解为：X1=13/4, x 2=5/2, Z=59/4 根据上表，写出非整数规划的约束方程，如: x2+1/2x 3-1/2x 4=5/2 (1) 将该方程中所有变量的系数及右端常数项均改写成 “整数与非负真分数之和”的形式，即： (1+0)x 2+(0+1/2)x 3+(-1+1/2)x 4=2+1/2 把整数及带有整数系数的变量移到方程左边，分数及 带有分数系数的变量称到方程右边，得： x2 - x 4-2 =1/2-(1/2x 3+1/2x 4) (2) 由于原数学模型已经 “标准化”，因此， 在整数最优解 中，X2和X4也必须取整数值，所以(2)式左端必为整数 或零，因而其右端也必须是整数。又因为 x3,x 4 0，所 以必有： 1/2-(1/2x 3+1/2x 4)1 由于 (2) 式右端必为整数，于是有： TOC \o 1-5 \h \z 1/2-(1/2x 3+1/2x 4) 0 (3) 或 x3+x4 1 (4) 这就是考虑整数约束的一个割平面约束方程，它是用 非基变量表示的，如果用基变量来表示割平面约束方 程，则有： 2x1+2x2 11 (5) 从图 1 中可以看出， (5) 式所表示的割平面约束仅割去 线性规划可行域中不包含整数可行解的部分区域，使 点 E(3.5 ， 2) 成为可行域的一个极点。 图1 在(3)式中加入松弛变量X5，得： -1/2x 3-1/2x 4+X5=-1/2 (6) 将⑹ 式增添到问题的约束条件中， 得到新的整数规划 问题： Max Z=3xi+2x2 2xi+3x2+X3=14 2xi+X2+X4=9 -1/2x 3-1/2x 4+X5=-1/2 Xi 0,且为整数，i=1,2,…,5 该问题的求解可以在表 1中加入(6)式，然后运用对偶 单纯形法求出最优解。具体计算过程见表 2: 表2 CB Xb b 3 2 0 0 0 X1 X2 % X4 X 2 X2 5/2 0 1 1/2 -1/2 C 3 X1 13/4 1 0 -1/4 3/4 0 0 X5 -1/2 0 0 [-1/2] -1 2 1 j 59/4 0 0 1/4 5/4 0 2 X2 2 0 1 0 -1 1 3 X1 7/2 1 0 0 1 -1/2 0 X3 1 0 0 1 1 -2 j 58/4 0 0 0 1 1/2 由此得最优解为：X1=7/2, x 2=2, z=58/4 该最优解仍不满足整数约束条件，因而需进行第二次 切割。为此，从表2中抄下非整数解xi的约束方程为： Xi+X4-1/2x 5 = 7/2 按整数、分数归并原则写成： TOC \o 1-5 \h \z xi +X4-X 5-3 = 1/2-1/2x 5 0 (7) 这就是一个新的割平面方程，用基变量来表示，得： X1+X2 5 (8) 在(7)中加入松弛变量X6,得： -1/2x 5+X6=-1/2 (9) 将(9)式增添到前一个问题的约束条件中去，得到又一 个新的整数规划问题，对它求解可以在表2中加入(7) 式，然后运用对偶单纯形法求出最优解。具体计算过 程见表3： 表3 CB X b 3 2 0 0 0 0 X1 X X X4 X5 X6 2 % 2 0 1 0 -1 1 0 3 X 7/2 1 0 0 1 -1/2 0 0 % 1 0 0 1 1 -2 0 0 X -1/2 0 0 0 0 [-1/2] 1 j 58/4 0 0 0 1 1/2 0 2 X2 1 0 1 0 -1 0 2 3 Xi 4 1 0 0 1 0 -1 0 X 3 0 0 1 1 0 -4 0 X 1 0 0 0 0 1 -2 14 0 0 0 1 0 由此得最优解为：Xi=4, x 2=1,z=14。该最优解符合整 数条件，因此也是原整数规划问题的最优解。 从图1中可以看出，由(8)式表示的割平面约束， 不 仅割去线性规划可行域中剩下的不含整数解域，而且 使最优整数解xi=4, X2=1 (即图2中的G点)，成为新