线性规划及其对偶问题课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《线性规划及其对偶问题课件.ppt》由用户(晟晟文业)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 线性规划 及其 对偶 问题 课件
- 资源描述:
-
1、第第2章章 线性规划及其对偶问题线性规划及其对偶问题(Linear Programming&its Dual Problem)2.1线性规划线性规划2.2 对偶理论对偶理论2.3 灵敏度分析灵敏度分析*2.4 LINGO软件求解线软件求解线 性规划模型的方法性规划模型的方法 2.1 线线 性性 规规 划划(Linear Programming)2.1.1 线性规划问题的数学模型线性规划问题的数学模型2.1.2 线性规划问题解的概念线性规划问题解的概念2.1.3 求解线性规划问题的图解法求解线性规划问题的图解法2.1.4 求解线性规划问题的单纯形法求解线性规划问题的单纯形法2.1.5 单纯形法的
2、进一步讨论单纯形法的进一步讨论2.1.6 线性规划模型的应用线性规划模型的应用为了完成一项任务或达到一定的目的,怎样用最少的人力、物力去完成或为了完成一项任务或达到一定的目的,怎样用最少的人力、物力去完成或者用最少的资源去完成较多的任务或达到一定的目的,这个过程就是规划。者用最少的资源去完成较多的任务或达到一定的目的,这个过程就是规划。例、有一正方形铁皮,如何截取例、有一正方形铁皮,如何截取 x 使容积为最大?使容积为最大?xa 22xxav 此为无约束极值问题此为无约束极值问题22(2)(2)(2)0axxax 6ax 0 dxdv2.1.1 线性规划问题的数学模型线性规划问题的数学模型(一
3、)、问题的提出(一)、问题的提出 设设 备备产产 品品 A B C D利润(元)利润(元)2 1 4 0 2 2 2 0 4 3 有有 效效 台台 时时 12 8 16 12例、已知资料如下表例、已知资料如下表所示,问如何安排生所示,问如何安排生产才能使利润最大?产才能使利润最大?或如何考虑利润大,或如何考虑利润大,产品好销。产品好销。模模 型型max Z=2x1+3x2 x1 0,x2 0s.t.2x1+2x2 12 x1+2x2 8 4x1 16 4x2 12此为带约束的极值问题此为带约束的极值问题问题中总有未知的变量,需要我们去解决问题中总有未知的变量,需要我们去解决要求:有目标函数及约
4、束条件,一般有非负条件存在,由此组成规划要求:有目标函数及约束条件,一般有非负条件存在,由此组成规划数学模型。数学模型。如果在规划问题的数学模型中,变量是连续的(数值取实数)其目标如果在规划问题的数学模型中,变量是连续的(数值取实数)其目标函数是有关线性函数(一次方),约束条件是有关变量的线性等式或函数是有关线性函数(一次方),约束条件是有关变量的线性等式或不等式,这样,规划问题的数学模型是线性的。反之,就是非线性的不等式,这样,规划问题的数学模型是线性的。反之,就是非线性的规划问题(其中一个条件符合即可)。规划问题(其中一个条件符合即可)。(二)、数学模型(二)、数学模型 1、1 12 21
5、1 112 2111 12 21max(min)(,)(,)0,0n nn nmmmn nmnZcxc xc xa xa xa xba xa xa xbxx 目标函数:目标函数:约束条件:约束条件:2、线性规划数学模型的一般形式、线性规划数学模型的一般形式也可以记为如下形式也可以记为如下形式:11 max(min)Z (,)(i1,2,)0 (j 1,2,)njjjnijjijjc xa xbmxn 目标函数:目标函数:约束条件:约束条件:如将上例用表格表示如下:如将上例用表格表示如下:设变量设变量 (1,2,)jxjn 产产 品品 j 设设 备备 i 有效台时有效台时 利润利润 mnmijn
6、aaaaa 1111n 2 1m 2 1 mbbb 21nccc 21jcib向向 量量 形形 式:式:12(,)nCccc nxxX1 mjjjaap1 mbbb1max (min)(,).0 jjZCXp xbs tX 矩阵形式:矩阵形式:mnmnaaaaA1111m ax (m in)(,)0 ZC XA XbX 规规划划确定型确定型随机型随机型静态规划静态规划动态规划动态规划线性规划线性规划非线性规划非线性规划 整数规划整数规划 非整数规划非整数规划整数规划整数规划非整数规划非整数规划3、规划类型、规划类型一一 般般 有有两种方法两种方法图图 解解 法法单纯形法单纯形法两个变量、直角坐
7、标两个变量、直角坐标三个变量、立体坐标三个变量、立体坐标适用于任意变量、但需将适用于任意变量、但需将一般形式变成标准形式一般形式变成标准形式2.1.2 线性规划问题解的概念线性规划问题解的概念(一)、求解方法(一)、求解方法1、解的概念、解的概念 可行解:满足约束条件可行解:满足约束条件、的解为可行解。所有解的集的解为可行解。所有解的集合为可行解的集或可行域。合为可行解的集或可行域。最优解:使目标函数达到最大值的可行解。最优解:使目标函数达到最大值的可行解。基:基:B是矩阵是矩阵A中中mm阶非奇异子矩阵(阶非奇异子矩阵(B 0),则),则B是一个基。是一个基。111121(,)mmmm maa
8、Bp ppaa则称则称 Pj(j=1,2,m)为基向量。为基向量。Xj 为基变量,否则为非基变量。为基变量,否则为非基变量。(二)、线性规划问题的解(二)、线性规划问题的解 基本解基本解(基解基解):令所有非基变量的值全为:令所有非基变量的值全为0,求得的基变,求得的基变量解与这些量解与这些0非基变量所组成的解非基变量所组成的解 X=(x1,x2,xm,0,0)T称为称为LP的基本解,最多为的基本解,最多为 个。个。Cmn 基本可行解:满足非负约束条件的基本解,简称基可行解基本可行解:满足非负约束条件的基本解,简称基可行解 可行基:对应于基可行解的基称为可行基。可行基:对应于基可行解的基称为可
9、行基。非可行解非可行解可可行行解解基解基解基可行解基可行解2、解的基本定理、解的基本定理 若线性规划问题存在可行解,则问题的可行域是凸集若线性规划问题存在可行解,则问题的可行域是凸集(凸多边形)(凸多边形)凸集凸集凸集凸集不是凸集不是凸集顶顶 点点 最优解一定是在凸集的某一顶点实现(顶点数目不超最优解一定是在凸集的某一顶点实现(顶点数目不超过过 个)个)Cmn 先找一个基本可行解,与周围顶点比较,如不是最大,先找一个基本可行解,与周围顶点比较,如不是最大,继续比较,直到找出最大为止。继续比较,直到找出最大为止。3、解的情况、解的情况唯唯 一一 解解无无 穷穷 解解无无 界界 解解无可行解无可行
10、解有最优解有最优解无最优解无最优解建立直角坐标建立直角坐标 ,图中阴影部分及边界上的点均为,图中阴影部分及边界上的点均为其解,是由约束条件来反映的。其解,是由约束条件来反映的。例例2.3)0 ,(21xx1212121212m a x23221 228.4 1 6 41 20,0Zxxxxxxs txxxx2.1.3 求解线性规划问题的图解法求解线性规划问题的图解法01 2 3 4 5 6 7 8 1 2 3 4 5 6 作作 图图 有唯一最有唯一最 优优 解:解:x1=4,x2=2最优目标函数值最优目标函数值 Z*=14x2 x1(4 2)00124 16 48212223221212121
11、21xxxxxxxxxxZ,max 例、例、例、例、121212212max2263212.20,0Zxxxxxxs txxx无穷多最优解无穷多最优解12121212max22.1,0 Zxxxxs txxxx 无界解无界解x1x1x2 x2 1121212min32 1.236,0 Zxxxxs txxxxx1x2 无可行解无可行解例、例、(一)基本思想(一)基本思想将模型的一般形式变成标准形式,再根据标准型模型,从可行域中找将模型的一般形式变成标准形式,再根据标准型模型,从可行域中找一个基本可行解,并判断是否是最优。如果是,获得最优解;如果不一个基本可行解,并判断是否是最优。如果是,获得最
12、优解;如果不是,转换到另一个基本可行解,当目标函数达到最大时,得到最优解。是,转换到另一个基本可行解,当目标函数达到最大时,得到最优解。(二)线性规划模型的标准形式(二)线性规划模型的标准形式max (i1,2,m).0 (j1,2,n)jjijjijZc xa xbstx1、标准形式、标准形式2.1.4 求解线性规划问题的单纯形法求解线性规划问题的单纯形法 2、特征:、特征:.目标函数为求极大值,也可以用求极小值;目标函数为求极大值,也可以用求极小值;.所有约束条件(非负条件除外)都是等式,右端常数项所有约束条件(非负条件除外)都是等式,右端常数项为非负;为非负;.变量为非负。变量为非负。3
13、、转换方式:、转换方式:.目标函数的转换目标函数的转换 如果是求极小值即如果是求极小值即 ,则可将目标函数乘以,则可将目标函数乘以(1),可化为求极大值问题。),可化为求极大值问题。jjxcZmin jjxcZZmax也就是:令也就是:令 ,可得到上式。,可得到上式。ZZ 即即.约束方程的转换:由不等式转换为等式约束方程的转换:由不等式转换为等式 ijijbxa0 iniinjijxbxxa称为松弛变量称为松弛变量 ijijbxa0 iniinjijxbxxa称为剩余变量称为剩余变量.变量的变换变量的变换 若存在取值无约束的变量若存在取值无约束的变量 ,可令,可令 其中:其中:jxjjjxxx
14、 0,jjxx例例2.2 将下列线性规划问题化为标准形式将下列线性规划问题化为标准形式123123123123123m in235 7 42.325,0,Zxxxxxxxxxs txxxxxx 为无约束(无非负限制)为无约束(无非负限制)解解:用用 替换替换 ,且,且 ,54xx 3x0,54 xx将第将第3个约束方程两边乘以个约束方程两边乘以(1)将极小值问题反号,变为求极大值将极小值问题反号,变为求极大值标准形式如下:标准形式如下:12456712456124571245124567max23()005()7()2.52()5,0 Zxxxxxxxxxxxxxxxxstxxxxx xxxx
15、x76,xx引入变量引入变量例、将下列线性规划问题化为标准型例、将下列线性规划问题化为标准型12121212m in2385.340,Zxxxxs txxxx 为无约束为无约束解:解:1341345134613456max2()38()5 .3()4 ,0 Zxxxxxxxs txxxxxxxxx(三)、单纯形法(三)、单纯形法例例2.31212121212m a x23221 228.4 1 6 41 2,0 Zxxxxxxs txxxx变成标准型变成标准型12345612312415max23000022 12 2 8 .4 16 Zxxxxxxxxxxxxs txx26123456 4
16、12 ,0 xxxxxxxx约束方程的系数矩阵约束方程的系数矩阵 654321100040010004001021000122ppppppA IppppB100001000010000165436543 xxxx,21xx,可作为为基变量可作为为基变量为非基变量为非基变量I 为单位矩阵且线性独立为单位矩阵且线性独立令:令:12x 16x 8x 12x 0654321 xx则:则:基本可行解为基本可行解为(0,0,12,8,16,12)此时,此时,Z=0然后,找另一个基本可行解。即将非基变量换入基变然后,找另一个基本可行解。即将非基变量换入基变量中,但保证其余的非负。如此循环下去,直到找到量中,
17、但保证其余的非负。如此循环下去,直到找到最优解为止。最优解为止。注意:为尽快找到最优解,在换入变量时有一定的要求。注意:为尽快找到最优解,在换入变量时有一定的要求。如将目标系数大的先换入等。如将目标系数大的先换入等。找出一个初始可行解找出一个初始可行解是否最优是否最优转移到另一个目标函数转移到另一个目标函数(找更大的基本可行解)(找更大的基本可行解)最优解最优解是是否否循循环环直到找出为止,核心是:变量迭代直到找出为止,核心是:变量迭代结束结束其步骤总结如下:其步骤总结如下:jjcnmmcccc11BcBXbmcc 1mxx 1mbb 1nmmxxxx 11im 1mnmmnmaaaa1,11
18、,110010 0 ijijjacc min(0)ilikiikbaa(四)、单纯形表(四)、单纯形表判定标准:判定标准:0 0:min0 0:maxz j jjjjjjjjijijjjczzcczzcaczc 若求若求 或或12345612312415max23000022 12 2 8 .4 16 Zxxxxxxxxxxxxs txx26123456 4 12 ,0 xxx xxxxx例例 题:题:cj 2 3 0 0 0 0cBxBb x1 x2 x3 x4 x5 x60000 x3x4x5x61281612 2 2 1 0 0 0 1 2 0 1 0 0 4 0 0 0 1 0 0 4
19、 0 0 0 1j i2 3 0 0 0 012/28/212/4cj 2 3 0 0 0 0cBxBb x1 x2 x3 x4 x5 x6000 x3x4x516 4 0 0 0 1 0 ji3x23010001/42620100-1/210010 0-1/2cj 2 3 0 0 0 0cBxBb x1 x2 x3 x4 x5 x60003x3x4x5x262163 2 0 1 0 0 -1/2 1 0 0 1 0 -1/2 4 0 0 0 1 0 0 1 0 0 0 1/4j i2 0 0 0 0 -3/46/2216/4cj 2 3 0 0 0 0cBxBb x1 x2 x3 x4 x5
20、 x60203x3 x1x5x22283 0 0 1 -2 0 1/2 1 0 0 1 0 -1/2 0 0 0 -4 1 2 0 1 0 0 0 1/44412j 0 0 0 -2 0 1/4icj 2 3 0 0 0 0cBxBb x1 x2 x3 x4 x5 x60203x6 x1x5x2 4402 0 0 2 -4 0 1 1 0 1 -1 0 0 0 0 -4 4 1 0 0 1 -1/2 1 0 0j 0 0 -1/2 -1 0 0i 0 0 0 -2 0 1/44412 0 0 1 -2 0 1/2 1 0 0 1 0 -1/2 0 0 0 -4 1 2 0 1 0 0 0 1/
21、42283x3 x1x5x20203 x1 x2 x3 x4 x5 x6bxBcB 2 3 0 0 0 0cjijcj 2 3 0 0 0 0cBxBb x1 x2 x3 x4 x5 x60203x3 x1x6 x2 0442 0 0 1 -1 -1/4 0 1 0 0 0 1/4 0 0 0 0 -2 1/2 1 0 1 0 1/2 -1/8 0j 0 0 0 -3/2 -1/8 0i 0 0 0 -2 0 1/44412 0 0 1 -2 0 1/2 1 0 0 1 0 -1/2 0 0 0 -4 1 2 0 1 0 0 0 1/42283x3 x1x5x20203 x1 x2 x3 x4
22、 x5 x6bxBcB 2 3 0 0 0 0cjij练习练习1212121223515.6224 ,0maxZxxxxs txxxx121231241234235 15.62 24 ,0maxZxxxxxs txxxxxxx 0 0 -1/12 -7/24 x2x112 x1 x2 x3 x4bxBcB 2 1 0 0cji15/43/4011/4-1/810-1/125/24j(一)、模型情况(一)、模型情况 变变 量:量:xj0 xj0 xj无约束无约束 结结1、组成、组成 约束条件:约束条件:=b 目标函数:目标函数:max min 果果2、变量、变量 xj0 令令 xj=-xj,xj
23、0 xj0 不处理不处理 xj 无约束无约束 令令xj=xj xj,xj0,xj0 唯一最优唯一最优无穷最优无穷最优无界解无界解无可行解无可行解2.1.5 单纯形法的进一步讨论单纯形法的进一步讨论3 3、约束、约束 条件:条件:bxpbxpbxpjjjjjj 加入松弛变量加入松弛变量加入人工变量加入人工变量先减去先减去 再加上再加上axaxsxsx例例2.512312312313123m in3 21 1423.2 1,0 Zxxxxxxxxxs txxxxx 1231234 12356137123min3 2 11 42 3 .2 1 ,0 Zxxxxxxxxxxxxs txxxxxx 4
24、4、目标函数:、目标函数:max,min 设规划模型约束条件为设规划模型约束条件为 ,需加入人工变量,需加入人工变量而得到一个而得到一个mm的单位矩阵,即基变量组合。因人工变量的单位矩阵,即基变量组合。因人工变量为虚拟变量,且存在于初始基本可行解中,需要将它们从基为虚拟变量,且存在于初始基本可行解中,需要将它们从基变量中替换出来。若基变量中不含有非零的人工变量,表示变量中替换出来。若基变量中不含有非零的人工变量,表示原问题有解。若当原问题有解。若当 ,而还有人工变量(非零)时,而还有人工变量(非零)时,则表示原问题无可行解。则表示原问题无可行解。axjjp xb0j加入人工变量后,目的是找到一
25、个单位向量,叫人工基。其加入人工变量后,目的是找到一个单位向量,叫人工基。其目标价值系数要确定,但不能影响目标函数的取值。一般可目标价值系数要确定,但不能影响目标函数的取值。一般可采用两种方法处理:大采用两种方法处理:大M法和两阶段法。法和两阶段法。.大大M法:法:即假定人工变量在目标函数中的系数为即假定人工变量在目标函数中的系数为M(任(任意大正数),如果是求极大值,需加意大正数),如果是求极大值,需加-M;如果是求极小值,;如果是求极小值,需加需加M。如基变量中还存在。如基变量中还存在M,就不能实现极值。,就不能实现极值。1234567 min300ZxxxxxMxMx 如上例:0 j计算
展开阅读全文