03 整数规划文档（土木工程系统分析）.ppt

第三章整数规划算例：水泥预制件运输最优方案建模：令x1，x2分别为预制梁、预制板的运输数量3.5IP模型求解第三章整数规划求解思路1：求得对应LP问题的最优解后取整不考虑x1，x2的整数限制→松弛问题3.5IP模型求解012x17564312x243第三章整数规划求解思路1：求得对应LP问题的最优解后取整不考虑x1，x2的整数限制→松弛问题结论：不可行IP问题→相应的LP问题LP问题的最优解→取整→IP问题的最优解3.5IP模型求解第三章整数规划求解思路2：穷举法若问题可行域有界→→穷举法x1：0、1、2、3、4；x2：0、1、23*5＝15可行、有效指派问题：n！ 4！＝24 10！＝3,628,80020！＝2432902008176640000不可行3.5IP模型求解第三章整数规划求解方法：分枝定界法割平面法3.5IP模型求解x1x2x1x2原问题问题B1问题B2x1原问题x2x1变化后问题x2第三章整数规划分枝定界法(BranchandBoundMethod)20世纪60年代初LanddoigDakin灵活、便于计算机求解可求解纯IP和混合IP问题，求解IP的重要方法基本原理通过增加成对的不等式右端为整数的约束（x1≥5，x1≤4），将原问题的可行域划分为两个子区域，得到两个问题的最优解，不断划分直到得到的最优解为整数解为止。3.5IP模型求解第三章整数规划分枝定界法——计算步骤松弛——分枝——剪枝3.5IP模型求解第三章整数规划分枝定界法——计算实例3.5IP模型求解x1＝4.81x2＝1.82z0＝356x1＝0x2＝0z0＝0x1x2123456780123456789109x1+7x2＝567x1+20x2＝70z＝40x1+90x2x1＝4.81x2＝1.82z0＝356x1x212345678012345678910BB2B3B1问题BX1＝4.81X2＝1.82Z0＝356问题B5X1＝5.44X2＝1.00Z0＝308问题B6无可行解问题B4X1＝1.42X2＝3.00Z0＝327问题B3X1＝4.00X2＝2.00Z0＝340问题B2X1＝5X2＝1.57Z0＝341问题B1X1＝4.00X2＝2.10Z0＝349x1x212345678012345678910分枝定界求解过程第三章整数规划0-1规划的隐枚举法(ImplicitEnumeration)3.5IP模型求解第三章整数规划0-1规划的隐枚举法(ImplicitEnumeration)3.5IP模型求解解条件满足条件Z值Θ1234（0,0,0）0×（0,0,1）5-1101√5（0,1,0）-2×（0,1,1）315×（1,0,0）31110√3（0,0,1）80211√8（1,1,0）1×（1,1,1）626×第三章整数规划0-1规划的隐枚举法(ImplicitEnumeration)3.5IP模型求解第三章整数规划割平面法（CuttingPlaneAlgorithm）R.E.Gomory提出，又称Gomory割平面法基本原理：通过增加线性约束（几何术语—割平面）将原可行域切掉一部分，这部分只包含非整数解，但没有切掉任何整数可行解。算法要点：如何构造割平面，使得切割后的可行域的一个整数坐标的极点刚好是问题的最优解。计算步骤：求解松弛LP问题的最优解→构造割平面→切割后再求解→直到LP问题的最优解刚好是整数解。3.5IP模型求解第三章整数规划割平面法（CuttingPlaneAlgorithm）计算实例3.5IP模型求解x1x2C（1,1）X1=3/4,x2=7/4,maxz=10/4Ci1100biXBX1X2X3X41X2013/4?7/41X110-1/41/43/400-1/2-1