数学建模整数规划分析.ppt

第四章 整数线性规划 4.1 汽车生产计划问题 4.2 分派问题 4.3 选课策略 y 问题二： 某厂每日8小时的产量不低于1800件。为了进行质量控制，计划聘请两种不同水平的检验员。一级检验员的标准为：速度25件/小时，正确率98%，计时工资4元/小时；二级检验员的标准为：速度15小时/件，正确率95%，计时工资3元/小时。检验员每错检一次，工厂要损失2元。为使总检验费用最省，该工厂应聘一级、二级检验员各几名？ 解 设需要一级和二级检验员的人数分别为x1、x2人, 则应付检验员的工资为： 因检验员错检而造成的损失为： 故目标函数为： 约束条件为： 且x1、x2均为整数 编写M文件xxgh4.m如下： c = [40;36]; A=[-5 -3]; b=[-45]; Aeq=[]; beq=[]; vlb = zeros(2,1); vub=[9;15]; %调用linprog函数： [x,fval] = linprog(c,A,b,Aeq,beq,vlb,vub) To Matlab (xxgh4) 结果为： x = 9.0000 0.0000 fval =360 即只需聘用9个一级检验员。 注：本问题应还有一个约束条件：x1、x2取整数。故它是一个整数线性规划问题。这里把它当成一个线性规划来解，求得其最优解刚好是整数：x1=9，x2=0，故它就是该整数规划的最优解。若用线性规划解法求得的最优解不是整数，将其取整后不一定是相应整数规划的最优解，这样的整数规划应用专门的方法求解。 例1 汽车厂生产计划 汽车厂生产三种类型的汽车，已知各类型每辆车对钢材、劳动时间的需求，利润及工厂每月的现有量。 小型 中型 大型 现有量 钢材（吨） 1.5 3 5 600 劳动时间（小时） 280 250 400 60000 利润（万元） 2 3 4 制订月生产计划，使工厂的利润最大。 设每月生产小、中、大型汽车的数量分别为x1, x2, x3 汽车厂生产计划 模型建立 小型 中型 大型 现有量 钢材 1.5 3 5 600 时间 280 250 400 60000 利润 2 3 4 线性规划模型(LP) 模型求解 3） 模型中增加条件：x1, x2, x3 均为整数，重新求解。 x = 64.5161 167.7419 0.0000 f = 632.2581 结果为小数，怎么办？ 1）舍去小数：取x1=64，x2=167，算出目标函数值z=629，与LP最优值632.2581相差不大。 2）试探：如取x1=65，x2=167；x1=64，x2=168等，计算函数值z，通过比较可能得到更优的解。 但必须检验它们是否满足约束条件 qiche1.m 注意：IP不能用Matlab软件直接求解 整数规划(Integer Programming,简记IP) 模型求解 怎么办？ 算法：1、割平面法 2、分支定界法 3、全枚举法 （应用较广泛） 全枚举法一般步骤： 1、确定决策变量数量和变化范围； 2、用for-end语句作决策变量的整型参数变化的循环，若有多个决策变量则要实现多重循环； 3、用if-end语句作不等式约束和等式约束条件是否满足的判断； 4、符合约束条件的一组决策变量，则进行目标函数计算，并进行存储，否则跳过 5、用函数max或min语句，搜索目标函数的最大值或最小值及相应的决策变量。 注意：当决策变量变化范围较大时，运行速度较慢 qiche2.m qiche2.c IP 的最优解x1=64，x2=168，x3=0，最优值z=632 注意：对于本题来说，决策变量的变化范围很大，故程序运行速度较慢。 分派问题 若干项任务分给一些候选人来完成，每人的专长不同，完成每项任务取得的效益或需要的资源