- 1、本文档共58页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
[理学]单纯形法
* 单纯形表 Cj C1 C2 Cj θ CB XB B-1b X1 X2 Xj P1 P2 pj σ(σ≤0?) A b 分别将系数C,矩阵A和资源限制b填入, 得到单纯形表的初表,计算检验数,进行最优性检验,若非优,就做初等行变换,形成新基的单纯形表 σ j=Cj-CBB-1Pj * 单纯形表的主要结构: 问题:第一张表的B-1=? ——单位阵I。 检验数的公式是什么? * 例5:用单纯形法求解例1 问题:标准模型的A中是否含I? ——松弛变量系数恰好构成I。 * [ ] [ ] 中表示进基列与出基行的交叉元,下一张表将实行以它为主元的初等行变换(称高斯消去)。方法是:先将主元消成1,再用此1将其所在列的其余元消成0。 * [ ] [ ] (请解释其实际意义) 此行元素除以主元10 30 0.3 1 0 0 0.1 此行元素乘以(-5)加入上行 此行元素除以主元2.5 * 练习:用单纯形法求解下面的线性规划 * [ ] 注:1. 表上每一列的含义: 2. 每张表上B-1的位置在哪?——对应于初表中I 的位置。 * [ ] [ ] 最优基解XB=(x3,x1,x2) 最优基B=(p3,p1,p2)在哪里? * [ ] 例6: 填表: 初表 终表 * 练习:用单纯形法求解下面的线性规划 * [ ] 问题:本题的单纯形终表检验数有何特点? —— 非基变量x2 的检验数等于零。 * 注:(1)解的几种情况在单纯形表上的体现(Max型): - 唯一最优解:终表非基变量检验数均小于零; - 多重最优解:终表非基变量检验数中有等于零的; - 无界解:任意表有正检验数相应的系数列均非正。 不同类型的解的判别: * (2) Min型单纯形表与Max型的区别仅在于:检验数反号,即 - 令负检验数中最小的对应的变量进基;minσ - 当检验数均大于等于零时为最优。 Min型的单纯形法 * 大M法 方法: 在一个标准的线性规划问题的约束条件中加进人工变量Xr并给定该人工变量在目标函数中的系数为(-M)(M为任意大的正数) 对于min型 有何不同? 最优性——最优基中不含人工变量 若基变量中不含有非零的人工变量,表示原问题有解。若对于MAX型,当 σ≤0 ,而还有人工变量(非零)时,则表示原问题无可行解。 设规划模型约束条件为AX ≥ b或者AX=b,此时需加入人工变量,以便在A矩阵中得到一个m×m的单位矩阵,使初始基可行。 * 例8 现有线性规划问题,试用大M法求解。 * 解 在上述问题的约束条件中加入松弛变量x4,剩余变量x5,人工变量x6,x7,得到 这里M是一个任意大的正数。 思考:能不能将第二个等式乘以-1? 不能。基解要可行,不能为负 * 因本例的目标函数是要求min,所以用所有cj-zj≥0来判别目标函数是否实现了最小化.用单纯形法进行计算时,初表如下: Cj -3 1 1 0 0 M M θ CB XB B-1b X1 X2 X3 X4 X5 X6 X7 P1 P2 P3 P4 P5 P6 P7 0 X4 11 1 -2 1 1 0 0 0 11 M X6 3 -4 1 2 0 -1 1 0 3/2 M X7 1 -2 0 [1] 0 0 0 1 1 σ(σj≥0?) -3+6M 1-M 1-3M 0 M 0 0 MINσ进基 * Cj -3 1 1 0 0 M M θ CB XB B-1b X1 X2 X3 X4 X5 X6 X7 P1 P2 P3 P4 P5 P6 P7 0 X4 10 3 -2 0 1 0 0 -1 -- M X6 1 0 [1] 0 0 -1 1 -2 1 1 X3 1 -2 0 1 0 0 0 1 -- σ(σj≥0?) -1 1-M 0 0 M 0 3M-1 0 X4 12 [3] 0 0 1 -2 2 -5 4 1 X2 1 0 1 0 0 -1 1 -2 -- 1 X3 1 -2 0 1 0 0 0 1 -- σ(σj≥0?) -1 0 0 0 1 M-1 M+1 * Cj -3 1 1 0 0 M M θ CB XB B-1b X1 X2 X3 X4 X5 X6 X7 P1 P2 P3 P4 P5 P6 P7 -3 X1 4 1 0 0 1/3 -2/3 2/3 -5/3 1 X2 1 0 1 0 0 -1 1 -2 1
文档评论(0)