matlab优化工具箱与lingo.ppt

第六讲 一、线性规划问题的求解 当规划中目标函数和约束条件都是决策变量的线性函数时，就是线性规划。 例如：min z=-2x1-x2+3x3-5x4 s.t. X1+2x2+4x3-x4=6 2x1+3x2-x3+x4=12 x1+x3+x4=4 x1,x2,x3,x4=0 1、利用matlab 优化工具箱中的linprog函数求解。 假设线性规划的数学模型为： 请同学们练习： max z=0.15x1+0.1x2+0.08x3+0.12x4 s.t. x1-x2-x3-x4=0 x2+x3-x4=0 x1+x2+x3+x4=1 x1,x2,x3,x4=0 思考：如果目标函数是max，约束是大于等于怎么办？ 将它转换为标准形式： min z=-0.15x1-0.1x2-0.08x3-0.12x4 s.t. x1-x2-x3-x4=0 -x2-x3+x4=0 x1+x2+x3+x4=1 x1,x2,x3,x4=0 求解代码： f=[-0.15,-0.1,-0.08,-0.12]; A=[1,-1,-1,-1;0,-1,-1,1]; b=[0;0]; Aeq=[1 1 1 1]; beq=[1]; lb=zeros(4,1); [x,fval]=linprog(f,A,b,Aeq,beq,lb) 2、利用lingo软件求解线性规划 对于线性规划：min z=-2x1-x2+3x3-5x4 s.t. X1+2x2+4x3-x4=6 2x1+3x2-x3+x4=12 x1+x3+x4=4 x1,x2,x3,x4=0 在lingo窗口中输入： min=-2*x1-x2+3*x3-5*x4; x1+2*x2+4*x3-x4=6; 2*x1+3*x2-x3+x4=12; x1+x3+x4=4; 注意lingo程序的输入规则！！！！！！！ 3、复杂一点的线性规划问题 假设某公司有6个货栈向8个销售商供应小装饰品，每一个货栈的供应量都是有限的，每一个销售商的需求量必须得到满足。该公司要决定如何调运货栈的装饰品满足销售商以使总运输成本最少。 数据： 货栈： 1 2 3 4 5 6 可供量： 60 55 51 43 41 52 销售商 1 2 3 4 5 6 7 8 需求量 35 37 22 32 41 32 43 38 运输成本： 1 2 3 4 5 6 7 8 1 6 2 6 7 4 2 5 9 2 4 9 5 3 8 5 8 2 3 5 2 1 9 7 4 3 3 4 7 6 7 3 9 2 7 1 5 2 3 9 5 7 2 6 5 6 5 5 2 2 8 1 4 3; 模型 Matlab求解此线性规划（代码）： f=[6;2;6;7;4;2;5;9;4;9;5;3;8;5;8;2;5;2;1;9;7;4;3;3;... 7;6;7;3;9;2;7;1;2;3;9;5;7;2;6;5;5;5;2;2;8;1;4;3]; A=[1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,0,0,0,0,0,0