4-LP对偶单纯形法12.pptVIP

表上作业法 初始基本可行解表上求解：西北角法 需求量 供给量 B1 B2 B3 B4 ai A1 6 4 5 9 10 A2 8 3 6 4 12 A3 5 8 7 8 8 bj 5 9 6 10 规则：优先安排标号小的产地和销地间的运输任务。 需求量 供给量 B1 B2 B3 B4 ai A1 5 6 5 4 X 5 X 9 10,5,0 A2 X 8 4 3 6 6 2 4 12,8,2,0 A3 X 5 X 8 X 7 8 8 8,0 bj 5 0 9 4 0 6 0 10 8 0 需求量 供给量 B1 B2 B3 B4 ai A1 6 4 5 9 10 A2 8 3 6 4 12 A3 5 8 7 8 8 bj 5 5 6 10 注意： 求解过程中，可能出现剩余供给和剩余需求 同时为0，这时只能划掉西北角元素所在的行列之一， 不能同时划掉，相应的基变量以0值出现在运输表上， 属于退化情形。 求 供 B1 B2 B3 B4 ai A1 5 6 5 4 X 5 X 9 10,5,0 A2 X 8 0 3 6 6 6 4 12,6,0 A3 X 5 X 8 X 7 4 8 8,0 bj 5 0 5 0 6 0 10 4 0 初始基本可行解: 求 供 B1 B2 B3 B4 ai A1 4 5 6 5 20 A2 7 10 5 6 20 A3 8 9 12 7 50 bj 15 25 20 30 练习：用西北角法求如下运输问题的初始基本可行解。 最优基本可行解表上求解法(单纯形算法应用) 运输问题的数学模型中，有一个约束是多余的，可去掉任一个！ 这个方程组行满秩，必有解。因此可求各个变量的判别数。 求 供 B1 B2 B3 B4 ai A1 5 8 3 4 10 A2 6 7 4 8 12 A3 7 6 8 5 9 bj 5 5 11 10 例 产销不平衡的运输问题 1.产大于销 产销不平衡的运输问题 2.产小于销 求 供 B1 B2 B3 B4 ai A1 6 4 3 7 9 A2 9 8 10 5 12 A3 4 7 6 10 14 bj 8 9 10 11 课堂练习 1.转化成产销平衡运输问题。 2.用西北角法求一基本可行解，并由此出发求最优解，是 总运费最小。 用MATLAB优化工具箱解线性规划 min z=cX 1、模型： 命令：x=linprog（c，A，b） 2、模型：min z=cX 命令：x=linprog（c，A，b，Aeq,beq） 注意：若没有不等式： 存在，则令A=[ ]，b=[ ]. 3、模型：min z=cX 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. 解 编写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) 解: 编写M文件如下： c=[-7 -5]; A=[3 2; 4 6; 0 7]; b=[90;200;210]; Aeq=[]; beq=[]; vlb=[0,0]; vub=[inf,inf]; [x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub) 问题2解答 实用优化方法 对偶单纯形法 设 是一个基， 满足约束 的对偶向量 称为对偶基