运筹学-大M法或两阶段法的上机实验
- 格式:doc
- 大小:179.02 KB
- 文档页数:11
1.4节中关于大M 法和两阶段法的课堂例题讨论课堂讨论题 习题1.5(2)P.51 用大M 法和两阶段法求解,34233min 2121212121≥=+≤+≥++=x x x x x x x x x x z标准化后为,,,34233min 43212142132121≥=+=++=-++=x x x x x x x x x x x x x x z采用大M 法数学模型为,,,3423)(3min 62162142153216521≥=++=++=+-++++=x x x x x x x x x x x x x x x M x x z Λ单纯形表法求解如下:因此的最优解如下TX ]1030[*=,最优目标函数min z=3。
采用两阶段法。
第一阶段的数学模型是:,,,3423min 621621421532165≥=++=++=+-++=x x x x x x x x x x x x x x x Λω单纯形表法的求解如下第二阶段的数学模型基标准化了的模型。
通过第一阶段的计算,我们求得了一个可行基,即[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-==001011101342P P P B ,解出相应的基变量后,可的单纯形表如下上述结果和将第一阶段的最后一张表的人工变量列删除,并且将原问题的价值系数换上再计算检验数的结果一样。
x入基,上述问题是退化问题,若按Bland法则,第一次换基时应让1相应的过程如下:已得第一阶段最优解,第二阶段的可行基是[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-==011021111312P P P B初始单纯形表如下的最优解,与前述结果一样,但由于选择的换基次序不同(即寻优的方向不同),在第一阶段中多了一步迭代,同时第一阶段的最优解不是原问题(第二阶段)的最优解,又计算了一次迭代才得到最优解。
欢迎您的下载,资料仅供参考!致力为企业和个人提供合同协议,策划案计划书,学习资料等等打造全网一站式需求。
运筹学作业1请分别用大M 法和两阶段法解题12312312312312310151253956151525,,0Max z x x x x x x x x x x x x x x x =++++≤⎧⎪-++≤⎪⎨++≥⎪⎪≥⎩解: 1:大M 法将该线性规划转换为标准型为:Max z=10x1+15x2+12x3+0x4+0x5+0x6-Mx75x1+3x2+x3+x4=9 -5x1+6x2+15x3+x5=15 2x1+x2+x3-x6+x7=5 x1,x2,x3,x4,x5,x6,x7≧0由该单纯性表可知,所有的检验数σ=j j z c -<=0,而人工变量x7没有被迭代出去,x7=1/2≠0,所以原线性规划问题无可行解。
2:阶段法: Min ω=X7 5x1+3x2+x3+x4=9 -5x1+6x2+15x3+x5=15 2x1+x2+x3-x6+x7=5 X1,x2.x3,x4,x5,x6,x7≧0由表中最后一行的检验数可知:所有的非基变量检验数σ=j j z c -=0,而人工变量x7=1/2≠0,并没有被迭代出去,故原线性规划问题无可行解。
运筹学作业22.某职工从武汉调至上海工作,拟将行李用集装箱装箱,然后用汽车或轮船运往上海,目的是为了节省运费但只能任选其一种运输方式。
已知有关数据见下表:试引入整数变量和0-1决策变量等建立该问题的整数规划模型解:设X1,X2分别是甲乙货物的箱数,同时引入一个充分大的数M 以及0-1变量y,令:0,当采取车运方式y=1,当采取船运方式则该整数规划可用如下模型表示:Min z = (1-y)(100x1+150x2)+y(200x1+300x2)7x1+8x2≤35+yM20x1+15x2≤500+yM80x1+75x2≤85+(1-y)M35x1+42x2≤250+(1-y)MX1,X2整数。
§4.2 两阶段法与大M 法————初始可行基的求法求解线性规划的步骤是: 1) 已知一个初始基本可行解 2) 从初始基本可行解出发,写出单纯型表,求出进基离基变量,做主元消去法,求出一个新的基本可行解且使目标函数值得到改善。
3) 判断当前基本可行解是否是最优解 那末,当观察不出来初始基本可行解时,怎么办?下面介绍的方法是几种求初始基本可行解的方法4.2.1 两 阶 段 法mincxt s .b Ax =x ≥0其中A 是nm ⨯矩阵,b≥0。
若A 中有m 阶单位矩阵,则初始基本可行解立即得到。
比如,[]N I A m ,=,那么⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡=0b x x x N B就是一个基本可行解。
若A 中不包含m阶单位矩阵,就需要用某种方法求出一个基本可行解。
介绍两阶段法之前,先引入人工变量的概念。
设A 中不包含m阶单位矩阵,为使约束方程的系数矩阵中含有m阶单位矩阵,把每个方程增加一个非负变量,令b x Ax a =+ (4.2.2)x ≥0 ,a x ≥0即bx x I A a m =⎥⎦⎤⎢⎣⎡),( (4.2.3)x ≥0 ,ax≥0显然,⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡b x x a 0是4:2:3的一个基本可行解.向量x ®¸0是人为引入的,它的每个分量成为人工变量。
人变量与前面介绍过的松弛变量是两个不同的概念。
松弛变量的作用是把不等式约束改写成等式约束,改写前后的两个问题是等价的。
因此,松弛变量是“合法”的变量。
而人工变量的引入,改变了原来的约束条件从这个意义上讲,它们是“不合法”的变量。
第一阶段是用单纯形方法消去人工变量(如果可能的话):min a Tx es:t A x +x ®=b (4:2:1) x ¸0;x ®¸0其中e =(1;1;1;¢¢¢;1)T 是分量全是1的m 维列向量,x ®=(x n+1;¢¢¢;x n+m )T 是人工变量构成的m 维列向量。
实验报告实验课程名称运筹学
实验项目名称大M法或两阶段法的上机实验年级
专业
学生姓名
学号
00 学院
实验时间:年月日
实验内容(包括实验具体内容、算法分析、源代码等等):
1.书上P97页第6题:用大M 法和两阶段法求解下列线性规划问题。
max z=5;3213x x x ++ 约束条件:102x 4x x 321≥++, 16.x 2x -x 321≤+ A :大M 法
图1.1
图1.2
δ,得出目标函数的最优解x1=16,x2=0,由上面的结果可知,满足所求出的0
≤
j
x3=0,sx4=16,Rx5=0,sx=0,最优值是80。
当把M的值改为100000后,值还是一样的,这样就可以得出当M为100时,已经得出有效解。
B:两阶段法
图1.3
由图1.3可知,先进行线性规划的第一阶段,满足0≤j δ,且z 值为零,即说明存在一个可行解使得所有的人工变量都为零,此时x2=2.5,sx6=21,其余为0得出z=0。
接下来进行第二阶段,令z=5x1+x2+3x3-0sx4+0Rx5+0sx6,和大M 的分析方法一样,最终将得到满足0≤j δ时达到最优解:当x1=16,x2=0,x3=0,sx4=6,Rx5=0,sx6=0,最优值为80。
2.书上P97页第7题(4)大M 法和两阶段法求解下列线性规划问题 。
max z=;321x x 2x ++ 约束条件:,42x 2x 4x 321≥++ ,204x 2x 21≤+ ,162x 8x 4x 321≤++ A :大M 法
图2.1
图2.2
由上面的图 2.1可知,首先先输入数据即线性规划的系数如图 2.1所示令max z=321x x 2x ++-0sx4+0sx6+0sx7-MRx5; 进行下一次迭代,以同样的方法一直下去,直到所求出的为止0≤j δ,就可以得出目标函数的最优解:x1=4,sx4=12,sx6=12,其余为0时,最优值为8。
当把M 的值改为100000后,值还是一样的,这样就可以得出当M 为100时,已经得出有效解。
B:两阶段法
图2.3
10。