[理学]数学建模 线性动态规划问题 教程.ppt

NLP虽然可用现成的数学软件求解(如LINGO, MATLAB)，但是其结果常依赖于初值的选择。 方法3：化为非线性规划 非线性规划（Non- Linear Programming，简记NLP） 实践表明，本例仅当初值非常接近上面方法算出的最优解时，才能得到正确的结果。 若生产某类汽车，则至少生产80辆，求生产计划。 x1=0 或 ?80 x2=0 或 ?80 x3=0 或 ?80 应如何安排原油的采购和加工 ？ 例2 原油采购与加工 市场上可买到不超过1500吨的原油A： 购买量不超过500吨时的单价为10000元/吨； 购买量超过500吨但不超过1000吨时，超过500吨的 部分8000元/吨； 购买量超过1000吨时，超过1000吨的部分6000元/吨。 售价4800元/吨 售价5600元/吨 库存500吨 库存1000吨 汽油甲(A?50%) 原油A 原油B 汽油乙 (A?60%) 决策变量 目标函数 问题分析 利润：销售汽油的收入 - 购买原油A的支出 难点：原油A的购价与购买量的关系较复杂 甲(A?50%) A B 乙(A?60%) 购买x? x11 x12 x21 x22 4.8千元/吨 5.6千元/吨 原油A的购买量,原油A, B生产汽油甲,乙的数量 c(x) ~ 购买原油A的支出 利润(千元) c(x)如何表述？ 原油供应 约束条件 x? 500吨单价为10千元/吨； 500吨? x? 1000吨，超过500吨的8千元/吨； 1000吨? x? 1500吨，超过1000吨的6千元/吨。 目标函数 购买x? A B x11 x12 x21 x22 库存500吨 库存1000吨 目标函数中c(x)不是线性函数，是非线性规划； 对于用分段函数定义的c(x)，一般的非线性规划软件也难以输入和求解； 想办法将模型化简，用现成的软件求解。 汽油含原油A的比例限制 约束条件 甲(A?50%) A B 乙(A?60%) x11 x12 x21 x22 x1 , x2 , x3 ~以价格10, 8, 6(千元/吨)采购A的吨数 目标函数 只有当以10千元/吨的价格购买x1=500(吨)时，才能以8千元/吨的价格购买x2 方法1 非线性规划模型，可以用LINGO求解 模型求解 x= x1+x2+x3, c(x) = 10x1+8x2+6x3 500吨? x? 1000吨，超过500吨的8千元/吨 增加约束 x= x1+x2+x3, c(x) = 10x1+8x2+6x3 方法1：LINGO求解 Model: Max= 4.8*x11 + 4.8*x21 + 5.6*x12 + 5.6*x22 - 10*x1 - 8*x2 - 6*x3; x11+x12 x + 500; x21+x22 1000; x11 - x21 0; 2*x12 - 3*x22 0; x=x1+x2+x3; (x1 - 500) * x2=0; (x2 - 500) * x3=0; x1 500; x2 500; x3 500; x 0; x11 0; x12 0; x21 0; x22 0; x1 0; x2 0; x3 0; end Objective value: 4800.000 Variable Value Reduced Cost X11 500.0000 0.0000000E+00 X21 500.0000 0.0000000E+00 X12 0.0000000E+00 0.0000000E+00 X22 0.0000000E+00 0.0000000E+00 X1 0.1021405E-13 10.00000 X2 0.0000000E+00 8.000000 X3 0.0000000E+00 6.000000 X 0.0000000E+00 0.0000000E+00 LINGO得到的是局部最优解，还能得到更好的解吗？ 用库存的500吨原油A、500吨原油B生产汽油甲，不购买新的原油A，利润为4,800千元。 y1, y2 , y3=1 ~以价格10, 8, 6(千元/吨)采购A 增加约束 方法2 0-1线性规划模型，可用LINDO求解 y1,y2,y3 =0或1 OBJECTIVE FUNCTION VALUE 1)