4运筹学——线性规划单纯形法.ppt

教案要点 文 件 名：062OR05.PPT；Simlex1.XLS 授课时间：2007年10月11日 授课班级：元培学院信计、自动化05级 授课内容：LP问题的单纯形法 预备知识：初等变换、凸集合、Excel 复习可行解、基可行解，基及非基变量。 难 点：单纯形法的步骤 重 点：单纯形法的思路；步骤：初始表，检验数，判优，进基、比值、出基、迭代，表上操作；利用Excel。 下节预习：教材《单纯形》，运输问题。 运筹学第五讲 元培学院05级 解LP问题单纯形法 LP问题的最优解 LP问题的标准型 标准型LP问题的解 线性代数问题 LP问题的代数解法 LP问题的代数解法 此线性规划问题图解法 LP问题的代数解法 实际上这是在可行域的顶点中搜索最优解： 从一个顶点开始，是否好？不好如何寻找出更好的一个？ 解LP的单纯形法 此法的证明不作要求。 先标准化,列表操作,借助Excel 实现。 定义7： X(1) , X(2) , … ,X(k) 是n维欧氏空间中的k点,若有一组数μ1,μ2 , …, μk 满足 0 ≤ μi ≤1 (i=1,… ,k) 定义8 顶点——凸集D, 点 X?D，若找不到0?1 的?和两个不同的点X(1) , X(2) ?D 使得: X = ?X(1) + (1-?)X(2), 则称X为 D的顶点。 定理1：LP问题的可行解域一定是凸集。 证明：设LP问题的可行解域为集合D, D={ X| AX=b,X ? 0 } 任取X(1),X(2)?D,则 X=? X(1) +(1- ? )X(2) ? 0 (0 ? ? ? 1) 又因为 AX(1) = b， AX(2) = b ∴ AX = A[? X(1) + (1-?) X(2) ] = ? b +(1-?)b = [? +(1-?)] b= b ∴ X?D， ∴D为凸集。 引理1：D为有界凸多面集,X?D，X必可为D的顶点的凸组合 。 只须证明： D的k个顶点X(1) , … ,X(k) ，有 0 ≤ μi ≤1， 使 X= μ1 X(1) + … + μk X(k) 证明可用归纳法（略） LP问题解的性质小结 若(LP)问题有可行解，则可行解集(可行域)是凸集(可能有界，也可能无界)，有有限个顶点。 单纯形解法 判断：基本解是否是可行基？ 判断：可行基是否是最优基？ 如何找第一个可行基？ 不是最优基的可行基如何改进？(有限次内要能结束。) 具体如何操作?单纯形表! 例1、Max Z=40x1+50x2 先把问题标准化。 单纯形法解LP问题 求解资源有约束求Max的LP问题单纯形法基本步骤 x1+2x2 +x3 = 30 3x1+2x2 +x4 = 60 2x2 +x5 = 24 x1 , x2 , x3 , x4 , x5 ?0 1 2 1 0 0 3 2 0 1 0 0 2 0 0 1 P1 P2 P3 P4 P5 A= 例1、 基变量 非基变量 x1 x2 x3 x4 x5 X= b= 30 60 24 x3=30-( x1+2 x2) x4=60-(3x1+2 x2) x5 =24 -2 x2 B=(P3 P4 P5)=I 可逆基N=(P1 P2) 令非基变量为0,解出基变量得基可行解 X= = = 是一个可行解 XN 0 XB B-1 b 0 0 30 60 24 1 2 1 又如Ｂ1 =(P1 P2 P3 )= 3 2 0 0 2 0 |B1|=6≠0, B1也可逆 令x1 = x2 =0, x3=30, x4=60, x5=24 x1= 12 - (1/3 x4 -1/3 x5) x2= 12 - (1/2 x5 ) x3 = -6 - (- 1/3 x4 -2/3 x5 ) Z = 40x1 +50x2 = 40[12-(1/3