用MATLAB优化工具箱解线性规划.docxVIP

用 MATLAB 优化工具箱解线性规划 1、模型： min z=cX s.t. AX ? b 命令：x=linprog（c，A，b） 2、模型： min z ? cX s.t. AX ? b AeqX ? beq 命令：x=linprog（c，A，b，Aeq,beq） 注意：若没有不等式： AX ? b 存在，则令 A=[ ]，b=[ ]. 若没有等式约束, 则令 Aeq=[ ], beq=[ ]. 3、模型： min z ? cX s.t. AX ? b AeqX ? beq VLB ? X ? VUB 命令：[1] x=linprog（c，A，b，Aeq,beq, VLB，VUB） [2] x=linprog（c，A，b，Aeq,beq, VLB，VUB, X0） 注意：[1] 若没有等式约束, 则令 Aeq=[ ], beq=[ ]. [2]其中 X0 表示初始点4、命令：[x,fval]=linprog(…) 返回最优解ｘ及ｘ处的目标函数值fval. 例 1 max z ? 0.4x 1 ? 0.28x 2 ? 0.32x 3 ? 0.72x 4 ? 0.64x 5 ? 0.6x 6 s.t. 0.01x 1 ? 0.01x 2 ? 0.01x 3 ? 0.03x 4 ? 0.03x 5 ? 0.03x 6 ? 850 0.02x 1 ? 0.05x 4 ? 700 x 2 x ? 0.05x 5 ? 0.08x ? 100 ? 900 3 6 x j ? 0 j ? 1,2, 6 解 编写 M 文件小 xxgh1.m 如下： c=[-0.4 -0.28 -0.32 -0.72 -0.64 -0.6]; A=[0.01 0.01 0.01 0.03 0.03 0.03;0.02 0 0 0.05 0 0;0 0.02 0 0 0.05 0;0 0 0.03 0 0 0.08]; b=[850;700;100;900]; Aeq=[]; beq=[]; vlb=[0;0;0;0;0;0]; vub=[]; [x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub) 例 2 min z ? 6x 1 s.t. x 3x 2 x 4x 3 x  ? 120 1 2 3 x ? 30 1 0 ? x 2 ? 50 x ? 20 3 解: 编写 M 文件 xxgh2.m 如下： c=[6 3 4]; A=[0 1 0]; b=[50]; Aeq=[1 1 1]; beq=[120]; vlb=[30,0,20]; vub=[]; [x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub 例 3 （任务分配问题）某车间有甲、乙两台机床，可用于加工三种工件。 假定这两台车床的可用台时数分别为 800 和 900，三种工件的数量分别为 400、600 和 500，且已知用三种不同车床加工单位数量不同工件所需的台时数和加工费用如下表。问怎样分配车床的加工任务，才能既满足加工工件的要求，又使 加工费用最低？ 甲 0.4 1.1 1.0 13 9 10 800 乙 0.5 1.2 1.3 11 12 8 900 车床类 型单位工件所需加工台时数单位工件的加工费用工件 1工件 2工件 3工件 1工件 2工件 3可用台时数解 设在甲车床上加工工件 1、2、3 的数量分别为 x1、x2、 车床 类 型 单位工件所需加工台时数 单位工件的加工费用 工件 1 工件 2 工件 3 工件 1 工件 2 工件 3 可用台 时数 min z ? 13x 1 9x 2 ? 10x 3 ? 11x 4 ? 12x 5 8x 6 ??x ? x ? x1 4 x ?x ? 2 5 ? 4 0 0 ? 6 0 0 s.t. ?x ? x ? 5 0 0 ?3 6 ? 0.4x ? 1.1x ? x ? 8 0 0 ? 1 2 3 ?0.5x ? 4 ? 1.2x 5 ? 1.3x 6 ? 9 0 0 ??x ? i ? 0, i ? 1,2, ,6 编写 M 文件 xxgh3.m 如下: f = [13 9 10 11 12 8]; A = [0.4 1.1 1 0 0 0 0 0 0 0.5 1.2 1.3]; b = [800; 900]; Aeq=[1 0 0 1 0 0 0 1 0 0 1 0 0 0 1 0 0 1]; beq=[400 600 500]; vlb = zeros(6,1); vub=[]; [x,fval] = linprog(f,A,b,Aeq,beq,vlb,vub) 例 4．某厂每日 8 小时的产量不低于 1800 件。为了