货机装运模型.docxVIP

PAGE PAGE 10 货机装运模型 问题重述：一架货机有三个货舱：前舱、中舱和后舱。三个货舱所能装载的货物的最大重量 和体积有限制如下表所示。并且为了飞机的平衡，三个货舱共装载的货物重量必须与其最大 的容许量成比例。 前舱 中舱 后舱 重量限制（吨） 10 16 8 体积限制（立方米） 6800 8700 5300 现有四类货物共该货机本次飞行装运，货物的规格以及装运后获得的利润如下表 重量（吨） 空间（立方米/吨） 利润（元/吨） 货物 1 18 480 3100 货物 2 15 650 3800 货物 3 23 580 3500 货物 4 12 390 2850 应如何安排装运，使得货机本次飞行获利最大？ 模型假设： （1） 每种货物可以无限细分； （2） 每种货物可以分布在一个或者多个货舱内； （3） 不同的货物可以放在同一个货舱内，并且可以保证不留空隙。 模型建立： 决策变量：每种货物放在每个货舱内的重量。用 xij 表示第 i 种货物放在第 j 个货舱内的重量， i = 1,2,3,4 分别表示货物 1，货物 2，货物 3 和货物 4。 j = 1,2,3 分别表示前舱、中舱和后 舱。 前舱（10 吨） 中舱（16 吨） 后舱（8 吨） （6800 立方米） （8700 立方米） （5300 立方米） x11，x21，x31， x41 x12，x22，x32， x42 x13，x23，x33， x43 货物 1 货物 2 货物 3 货物 4 x11+x12+x13 x21+x22+x23 x31+x32+x33 x41+x42+x43 决策目标：总利润的最大化，目标函数为 3100( x11 + x12 + x13 ) + 3800( x21 + x22 + x23 ) + 3500( x31 + x32 + x33 ) + 2850( x41 + x42 + x43 ) 约束条件： （1） 供装载的四种货物的总重量约束， ?? x11 + x12 + x13 ≤ 18 ? ?x21 + x22 + x23 ≤ 15 ? ?x31 + x32 + x33 ≤ 23 x41 + x42 + x43 ≤ 12 （2） 三个货舱的空间限制 ?480x11 + 650x21 + 580x31 + 390x41 ≤ 6800 ? ??480x12 + 650x22 + 580x32 + 390x42 ≤ 8700 ? ?480x13 + 650x23 + 580x33 + 390x43 ≤ 5300 （3） 三个货舱的重量限制 ?x11 + x21 + x31 + x41 ≤ 10 ? ??x12 + x22 + x32 + x42 ≤ 16 ? ?x13 + x23 + x33 + x43 ≤ 8 （4） 三个货舱装入重量的平衡约束 x11 + x21 + x31 + x41 = x12 + x22 + x32 + x42 = x13 + x23 + x33 + x43 10 16 8 模型求解： 使用计算软件求解（在 MATLAB 中，可以使用 linprog 命令求解） 求解结果为： ( x1 ; x2 ; x3 ; x4 ) = (0,0,0;10,0,5;0,12.947,3, ;0,3.053,0) MATLAB 实现线性规划的运算 为了避免这种形式多样性带来的不便，Matlab 中规定线性规划的标准形式为 min  c T x  such that Ax ≤ b Aeq ? x = beq lb ≤ x ≤ ub 其中 c 和 x 为 n 维列向量， A 、 Aeq 为适当维数的矩阵， b 、 beq 为适当维数的列向量。 例如线性规划 的 Matlab 标准型为  max cT x x  such that  Ax ≥ b min x ? cT x  such that ? Ax ≤ ?b 基本函数形式为 linprog(c,A,b)，它的返回值是向量 x 的值。还有其它的一些函数调用形式（在 Matlab 指令窗运行 help linprog 可以看到所有的函数调用形式），如： [x,fval]=linprog(c,A,b,Aeq,beq,LB,UB,X0,OPTIONS) 这里 fval 返回目标函数的值，LB 和 UB 分别是变量 x 的下界和上界， x0 是 x 的初始值， OPTIONS 是控制参数。 注： linprog 的格式 X=