最优化方法线性规划单纯形法PPT学习教案

1、会计学1最优化方法线性规划单纯形法最优化方法线性规划单纯形法线性规划线性规划：目标函数是线性的，约束条件是目标函数是线性的，约束条件是线性等式或不等式线性等式或不等式线性规划线性规划第1页/共77页第2页/共77页第3页/共77页第4页/共77页第5页/共77页 问题问题：确定食品数量，满足营养需求，花费最小？：确定食品数量，满足营养需求，花费最小？ 变量：变量：n种食品，种食品，m种营养成份；第种营养成份；第 j 种食品的单价种食品的单价每单位第每单位第 j 种食品所含第种食品所含第 i 种营养的数量种营养的数量食用第食用第 j 种食品的数量种食品的数量为了健康，每天必须食用第为了健康，每天

2、必须食用第i 种营养的数量种营养的数量 模型：模型：第6页/共77页产销产销平衡平衡/不平衡不平衡的运输问题的运输问题第7页/共77页第8页/共77页线性规划的一般形式线性规划的一般形式第9页/共77页线性规划的标准形线性规划的标准形(分析、算法分析、算法)标准形的特征：标准形的特征：极小化极小化、等式约束等式约束、变量非负变量非负向量表示：向量表示：第10页/共77页一般形式一般形式 标准形标准形转化转化称称 松弛松弛(slack)/盈余盈余(surplus)变量；变量； 自由自由变量变量第11页/共77页等价表示为等价表示为第12页/共77页定义：定义： 给定含有给定含有n个变量，个变量，

3、m个方程的线性方程组个方程的线性方程组Ax=b，设设B是由是由A 的列组成的任一非奇异的列组成的任一非奇异mm子阵，则如果置子阵，则如果置x的所有与的所有与B无关的无关的n-m个分量为零后，所得方程组的解是个分量为零后，所得方程组的解是Ax=b关于基关于基B的的基本解基本解(basic solution) ，称，称x中与基中与基B对应对应的分量为的分量为基变量基变量(basic variables)退化基本解：退化基本解：基本解中如果有一个或多个基变量的值为基本解中如果有一个或多个基变量的值为零零基本解与基变量基本解与基变量其中其中满秩假定：满秩假定： m n矩阵矩阵A满足满足mn，且，且A的

4、行向量线性无关的行向量线性无关 在满秩假定下，方程组在满秩假定下，方程组Ax=b总有解，且至少有一个基总有解，且至少有一个基本本 解解第13页/共77页基本可行解基本可行解定义定义 称称 的非负基本解是的非负基本解是标准形标准形的的基本可行解基本可行解( (basic feasible solution)；约束系统约束系统 第14页/共77页线性规划的基本定理线性规划的基本定理i) 若标准型有可行解，则若标准型有可行解，则必有必有基本可行解；基本可行解；ii) 若标准型有最优解，则若标准型有最优解，则必有必有最优最优基本可行解。基本可行解。 考虑线性规划标准形，其中考虑线性规划标准形，其中A是

5、秩为是秩为m的的mn阶阶矩阵，则以下结论成立：矩阵，则以下结论成立：基本可行解的个数基本可行解的个数不超过不超过第15页/共77页与凸性的关系与凸性的关系线性规划的基本定理线性规划的基本定理( (标准形标准形) )基本可行基本可行解解线性方程线性方程组的基本组的基本性质性质代数理论代数理论(与与表述形式有关表述形式有关)设计算法设计算法极点极点凸集理论凸集理论几何理论几何理论( (与表述形式与表述形式无关无关) )直观理解直观理解第16页/共77页凸性凸性( (凸集及性凸集及性质质) )几何解释几何解释：连接集合中任两点的线段仍含在该集合中：连接集合中任两点的线段仍含在该集合中性质性质 定义定

6、义 是凸集是凸集(convex set)，如果对，如果对S中任意中任意 两两 点点 x , y 和和(0,1)中的任一数中的任一数 满足满足 第17页/共77页一些重要的凸集一些重要的凸集有限个闭半空间的交集有限个闭半空间的交集多面集多面集(polyhedral convex set):推推广广平面上：多边形平面上：多边形注：注：任一线性规划的可行集是任一线性规划的可行集是多面集多面集！超平面超平面(hyperplane):正正/ /负闭半空间：负闭半空间：第18页/共77页极点极点几何上几何上：极点即不能位于连接该集合中其它两点：极点即不能位于连接该集合中其它两点的开线段上的点的开线段上的点

7、定义定义 称凸集称凸集C中的点中的点 x 是是C的极点，如果存在的极点，如果存在 C 中的点中的点 y, z 和某和某 ，有，有则必有则必有 y=z.第19页/共77页极点与基本可行解的等价性定理极点与基本可行解的等价性定理推论：推论：i) 若若K非空，则至少有一个极点非空，则至少有一个极点.ii) 若线性规划有最优解，则必有一个极点是最优解若线性规划有最优解，则必有一个极点是最优解.iii) Ax=b对应的约束集对应的约束集K最多有有限个极点最多有有限个极点.考虑线性规划标准形，其中考虑线性规划标准形，其中A是秩为是秩为m的的mn矩阵，令矩阵，令则则x是是 K 的极点，的极点，当且仅当当且仅

8、当x是线性规划的基本可行解是线性规划的基本可行解.（线性规划基本定理的几何形式）（线性规划基本定理的几何形式）第20页/共77页例例2.2. K有有2个极点个极点有有3个基本解，个基本解，2个个可可行行K 有有3个极点个极点有有3个基本解，个基本解，均可行均可行例例1. 1. 第21页/共77页例例3. 3. Subject to5个极点个极点极点极点第22页/共77页线性规划线性规划解的解的几何特征几何特征唯一唯一 解解(顶点顶点)！第23页/共77页线性规划解的线性规划解的几何特征几何特征无解无解可行集：可行集：多边形多边形( (二维二维) )多边集多边集( (高维空间高维空间) )给出给

9、出有效的代数刻画有效的代数刻画和和严谨的几何描述严谨的几何描述，从理论上证，从理论上证实上述几何特征，并实上述几何特征，并寻求有效算法寻求有效算法 有解：有解：唯一解唯一解/ /多个解多个解( (整条边、面、甚至整条边、面、甚至整个可行集整个可行集) ) 有顶点解有顶点解第24页/共77页顶点顶点一条边一条边无无(下下)界界第25页/共77页第26页/共77页 初始化：初始化：如何找到一个如何找到一个BFS？ 判断准则：判断准则：何时最优？何时无界？何时最优？何时无界？ 迭代规则：迭代规则：如何从一个极点如何从一个极点/BFS迭代到相迭代到相邻极点邻极点/BFS？第27页/共77页1.转轴转轴

10、(基本解基本解相邻相邻基本解基本解)满秩假定：满秩假定：A是行满秩是行满秩的的第28页/共77页规范形规范形(canonical form)基变量基变量基本解基本解非基变量非基变量等价变形等价变形不妨设不妨设 线性无关线性无关 一般地一般地：次序可以打乱！次序可以打乱！只要有只要有m个单位列个单位列第29页/共77页规范形的转换问题规范形的转换问题 什么时候可以替换？什么时候可以替换？ 替换后替换后新新规范形规范形是什么？是什么？ 替换问题替换问题假设在上述规范形中，想假设在上述规范形中，想用用第30页/共77页转轴转轴(pivot) 当且仅当当且仅当 ，可以替换，可以替换 替换后，新规范形的

11、系数替换后，新规范形的系数转轴公式转轴公式转轴元转轴元(pivot element)第31页/共77页转轴例例1. 求下列方程组以求下列方程组以 为基变量的基本为基变量的基本解解第32页/共77页转轴转轴转轴转轴转轴转轴x=(0,0,0,4,2,1)第33页/共77页2.BFS相邻相邻BFS(极点极点相邻相邻极点极点)问题：问题：确定出基变量，使转轴后确定出基变量，使转轴后新规范形新规范形对应对应BFS？设设x是是BFS, 且规范形如前，且且规范形如前，且假设假设 aq 进基进基因因为为令令可否选取合适的可否选取合适的 使得使得 是是BFS ？所所以以第34页/共77页确定离基变量确定离基变量

12、至少有一个正元至少有一个正元第35页/共77页例例3. 考虑线性方程组考虑线性方程组a4进基进基转轴转轴B=(a1,a2,a3)X=(4,3,1,0,0,0)x=(0,1,3,2,0,0)第36页/共77页3. BFS目标值减小的相邻目标值减小的相邻BFS 将将Ax=b的任一解的任一解用非基变量用非基变量表示；表示；设设x是是BFS，且规范形如前，则有，且规范形如前，则有问题：问题：确定进基变量，转轴后使确定进基变量，转轴后使新新BFS的目标值的目标值变小变小？用非基变量表示用非基变量表示. 选取进基变量的依据选取进基变量的依据 将将目标函数目标函数第37页/共77页相对相对/既约费用系数既约

13、费用系数(relative/reduced cost coefficients)将将 Ax=b 的任一解的任一解 x 用非基变量表示为用非基变量表示为度量变量相对于给定基的费用度量变量相对于给定基的费用第38页/共77页确定进基变量确定进基变量最优性定理最优性定理定理定理(BFS的提高的提高)相对费用系数的相对费用系数的经济解释经济解释！(合成价格、相对价格合成价格、相对价格) 给定目标值为给定目标值为z0的的非退化非退化基本可行解，且假定存基本可行解，且假定存在在 j 使得使得 rj 0，无无可行解！可行解！ z*= 0，有有可行解！可行解！ 基变量中基变量中无无人工变量人工变量x 是是BF

14、S，B 是对应的基是对应的基 基变量中基变量中有有人工变量人工变量驱赶人工变量出基驱赶人工变量出基假设第假设第 i 个基变量是人工变量，且当前单纯形表个基变量是人工变量，且当前单纯形表第第 i 行的前行的前n个数据是个数据是第第 i 个约束冗余；个约束冗余；删除单纯形表的第删除单纯形表的第 i 行数行数据据以以任一非零元任一非零元为转轴元转轴为转轴元转轴得辅助问题的一个新最优得辅助问题的一个新最优BFS，且基变量中少，且基变量中少1个人工变量！个人工变量！第55页/共77页例例1. 给出下面系统的一个基本可行解，或者说明其无解给出下面系统的一个基本可行解，或者说明其无解引入引入人工人工变量变量

15、目标目标：辅助问题辅助问题的的初始表初始表格格！BFS第56页/共77页第一张第一张单纯形表单纯形表第二张第二张单纯形表单纯形表注意基变量整列包括末行z在内除了基变量其他元素都是0第57页/共77页辅助问题的辅助问题的最优值是最优值是0 0. . 原问题的原问题的BFS：第58页/共77页两阶段法两阶段法可求任一可求任一线性规划问题线性规划问题第59页/共77页例例2. 利用两阶段法求解下面的问题利用两阶段法求解下面的问题辅助问题辅助问题第第I阶段：阶段：先构造辅助向量z=x4+x5第60页/共77页辅助问题的辅助问题的最后一张单纯形表最后一张单纯形表原问题的原问题的初始初始表格：表格：得到辅

16、助问题的最后一张单纯形表后，去掉辅助变量，将原始问题的z带入表格，启动单纯形法第61页/共77页原问题的原问题的最优解最优解：第62页/共77页6. 修正单纯形法修正单纯形法(Revised simplex method) 重要事实：重要事实： 通常仅有少数列发生转轴通常仅有少数列发生转轴(2m-3m) 核心问题：核心问题：如何更新如何更新当前基的逆当前基的逆新基的逆新基的逆理论上的表现理论上的表现表格实现表格实现 仅需原始数据仅需原始数据(c, A, b)和基和基 B 的逆矩阵的逆矩阵第63页/共77页7. 单纯形法的矩阵形式单纯形法的矩阵形式给定基给定基 B 及对应及对应BFS (xB,

17、0), 其中其中xB=B-1b用用非基非基变量表示变量表示目标函数目标函数：用用非基非基变量表示变量表示基基变量：变量：相对相对费用向量费用向量第64页/共77页初始表格单纯形表初始表格单纯形表初始表格初始表格通常不是单纯形表！通常不是单纯形表！与基矩阵与基矩阵 B 对应的对应的单纯形表单纯形表第65页/共77页修正单纯形法的计算步骤修正单纯形法的计算步骤步步2 选取选取 q 满足满足步步3 计算计算 yq=B-1aq；若若步步1 计算计算 。如果。如果 停；得最优解停；得最优解.步步0 给定给定BFS及对应的及对应的B-1。计算。计算核心核心计算：计算：B-1涉及涉及到的计算：到的计算： ， 停，停，问题问题无界无界；否则，选；否则，选 p 满足满足步步4 更新更新 B-1, B-1b和和 ，返步，返步1.第66页/共77页基的转换定理基的转换定理左乘左乘该矩阵等价于对矩阵进行该矩阵等价于对矩阵进行初等行变换初等行变换！定理定理 不妨设不妨设B= . 则则 aq 进基，进基，ap出基后所得新基出基后所得新基 的逆的逆这里这里 ei 表示表示n 维单位向量，向量维单位向量，向量 v 定义定义为为第67页/共77页相关数据的更新相关数据的更新初等行变换初等行变换设设转轴元转轴元是，即是，即 aq