整数规划与matlab吴颖丹.ppt

数学建模培训;整数规划;精品资料; 你怎么称呼老师？ 如果老师最后没有总结一节课的重点的难点，你是否会认为老师的教学方法需要改进？ 你所经历的课堂，是讲座式还是讨论式？ 教师的教鞭 “不怕太阳晒，也不怕那风雨狂，只怕先生骂我笨，没有学问无颜见爹娘 ……” “太阳当空照，花儿对我笑，小鸟说早早早……”;例一、合理下料问题 设用某型号的圆钢下零件A1, A2,…,Am 的毛坯。在一根圆钢上下料的方式有B1,B2, … Bn 种，每种下料方式可以得到各种零件的毛坯数以及每种零件的需要量，??表所示。怎样安排下料方式，使得即满足需要，所用的原材料又最少？; 设：xj 表示用Bj (j=1.2…n) 种方式下料根数 模型：;例三、机床分配问题 设有m台同类机床，要加工n种零件。已知各种零件的加工时间分别为a1,a2,…an ，问如何分配，使各机床的总加工任务相等，或者说尽可能平衡。 ;;因此，求xij ，使得;整数规划一般形式;纯整数规划：所有决策变量要求取非负整数（这时引进的松弛变量和剩余变量可以不要求取整数）。 全整数规划：除所有决策变量要求取非负整数外，系数aij和常数bi也要求取整数（这时引进的松弛变量和剩余变量也必须是整数）。 混合整数规划：只有一部分的决策变量要求取非负整数，另一部分可以取非负实数。 0－1整数规划：所有决策变量只能取 0 或 1 两个整数。;从数学模型上看整数规划似乎是线性规划的一种特殊形式，求解只需在线性规划的基础上，通过舍入取整，寻求满足整数要求的解即可。但实际上两者却有很大的不同，通过舍入得到的解（整数）也不一定就是最优解，有时甚至不能保证所得倒的解是整数可行解。;例：设整数规划问题如下 ;用图解法求出最优解 x1＝3/2, x2 = 10/3 且有Z = 29/6;因此，可将集合内的整数点一一找出，其最大目标函数的值为最优解，此法为完全枚举法。 如上例：其中（2，2）（3，1）点为最大值，Z=4。 目前，常用的求解整数规划的方法有： 分支定界法和割平面法； 对于特别的0－1规划问题采用隐枚举法和匈牙利法。;考虑纯整数问题：;1、先不考虑整数约束，解( IP )的松弛问题( LP )，可能得到以下情况之一： ⑴.若( LP )没有可行解，则( IP )也没有可行解，停止计算。 ⑵.若( LP )有最优解，并符合( IP )的整数条件，则( LP )的最优解即为( IP )的最优解，停止计算。 ⑶.若( LP )有最优解，但不符合( IP )的整数条件，转入下一步。为讨论方便，设( LP )的最优解为： ; 2、定界： 记( IP )的目标函数最优值为Z* ,以Z(0) 作为Z* 的上界，记为 ＝Z(0) 。再用观察法找的一个整数可行解 X′，并以其相应的目标函数值 Z′作为Z* 的下界，记为Z＝ Z′，也可以令Z＝－∞，则有: Z ≤ Z* ≤; ; 为了更好地说明用分枝定界法求整数规划最优解的过程，我们选择只有两个变量的引例3.2.1 进行求解。 例3.3.1 用分枝定界法求解整数规划问题 (P)：max Z = 15x1 ? 20x2 s.t. 6x1 ? 4x2 ? 25 x1 ? 3x2 ? 10 xi ? {0, 1, 2, …} ;解：1. 求解如下的松弛问题 (P0)： (P0)：max Z = 15x1 ? 20x2 s.t. 6x1 ? 4x2 ? 25 x1 ? 3x2 ? 10 xi ? 0，i = 1, 2 得最优解 x1 = 2.5，x2 = 2.5，最优值 z = 87.5。由于原问题 (P) 目标函数的系数为整数，xi? {0, 1, 2, …}，故 z* ? 87，即最优值的上界 = 87。令最优值的下界 z = 0，则有 z = 0 z* ? 87 = 。我们将这些结果记录在树形图图 3.3.3 中。 ; 2. 因为此时两个变量都不是整数，我们从中选择一个变量进行分枝。假定选择 x1，在 (P0) 的约束之外，增加两个互相排斥的