第四篇 数学规划模型.jsp.pptVIP

model: sets: pattern/1..7/:residu,number; type/1..3/:required; link(pattern,type):aa; endsets data: residu=3 1 3 3 1 1 3; required=50 20 15; aa=4 0 0 3 1 0 2 0 1 1 2 0 1 1 1 0 3 0 0 0 2; enddata !目标函数(总余量最小) min=@sum(pattern:residu*number); !工地需求; @for(type(i): @sum(pattern(j): aa(j,i)*number(j))required(i)); !各变量整数约束; @for(pattern:@gin(number)); end 钢管切割问题的Lingo程序 当余料没有用处时，通常以总根数最少为目标 目标2（总根数） 钢管下料问题1 约束条件不变 最优解：x2=15, x5=5, x7=5, 其余为0； 最优值：25。 xi 为整数 按模式2切割15根，按模式5切割5根，按模式7切割5根，共25根，余料35米 虽余料增加8米，但减少了2根 与目标1的结果“共切割27根，余料27米” 相比 板材规格2： 长方形， 32?28cm， 2万张。 例2 易拉罐下料 每周工作40小时，每只易拉罐利润0.10元，原料余料损失0.001元 / cm2（不能装配的罐身、盖、底也是余料） 模式1：1.5秒 模式2：2秒 模式3：1秒 模式4：3秒 上盖 下底 罐 身 罐身高10cm，上盖、下底直径均5cm。 板材规格1： 正方形，边长24cm，5万张。 如何安排每周生产？ ? 罐身个数 底、盖 个数 余料损失 （cm2） 冲压时间（秒） 模式1 1 10 222.6 1.5 模式2 2 4 183.3 2 模式3 0 16 261.8 1 模式4 4 5 169.5 3 模式1: 正方形 边长24cm 问题分析 计算各种模式下的余料损失 上、下底直径d=5cm，罐身高h=10cm。 模式1 余料损失 242-10??d2/4 - ?dh=222.6 cm2 问题分析 目标:易拉罐利润扣除原料余料损失后的净利润最大 约束：每周工作时间不超过40小时； 原料数量：规格1（模式1 ~3）5万张， 规格2（模式4）2万张； 罐身和底、盖的配套组装 。 注意：不能装配的罐身、上下底也是余料 决策变量 xi ~ 按照第i 种模式的生产张数(i=1,2,3,4)； y1 ~ 一周生产的易拉罐个数； y2 ~ 不配套的罐身个数； y3 ~ 不配套的底、盖个数。 模型建立 目标 约束条件 时间约束 原料约束 产量 余料 时间 x1 222.6 1.5 x2 183.3 2 x3 261.8 1 x4 169.5 3 模型建立 y1 ~ 易拉罐个数；y2 ~ 不配套的罐身； y3 ~ 不配套的底、盖。 每只易拉罐利润0.10元，余料损失0.001元 / cm2 罐身面积?dh=157.1 cm2 底盖面积?d2/4=19.6 cm2 (40小时) 约束条件 配套约束 y1 ~ 易拉罐个数；y2 ~ 不配套的罐身； y3 ~ 不配套的底、盖。 罐身 底、盖 1 10 2 4 0 16 4 5 产量 x1 x2 x3 x4 虽然xi和y1，y2，y3应是整数，但是因生产量很大，可以把它们看成实数，从而用线性规划模型处理 。 将所有决策变量扩大10000倍（xi ~万张，yi ~万件） LINDO发出警告信息：“数据之间的数量级差别太大，建议进行预处理，缩小数据之间的差别” 模式2生产40125张， 模式3生产3750张， 模式4生产20000张， 共产易拉罐160250个 (罐身和底、盖无剩余)， 净利润为4298元 模型求解 OBJECTIVE FUNCTION VALUE 1) 0.4298337 VARIABLE VALUE REDUCED COST Y1 16.025000 0.000000 X1 0.000000 0.000050 X2