习题一 P46 (a)
x2 4x12x244 3 2 1 0 1 2 3 4x16x26x1
1的所有x1,x2,此时目标函数值2该问题有无穷多最优解,即满足4x16x26且0x2z3。 (b)
x2 3 2 0 1 4 x1
用图解法找不到满足所有约束条件的公共范围,所以该问题无可行解。
(a) 约束方程组的系数矩阵
1236300A814020
300001基 p1 p2 p3 p1 p2 p4 p1 p2 p5 基解 x1 x2 x3 x4 x5 x6 1670 - 0 0 0 36是否基可行解 否 是 是 目标函数值 0 10 0 7 0 0 70 3 0 0 0 210 3 p1 p2 p6 p1 p3 p4 p1 p3 p5 p1 p3 p6 p1 p4 p5 p1 p4 p6 721 4 0 0 0 4450 0 8 0 0 230 0 0 8 0 211 0 0 0 3 2否 否 是 否 是 否 3 0 0 0 0 3 5 0 515 0 0 2 0 44最优解x0,10,0,7,0,0T。 (b) 约束方程组的系数矩阵
1234A2212 基 p1 p2 p1 p3 p1 p4 p2 p3 p2 p4 p3 p4 基解 x1 x2 x3 x4 114 0 0 2是否基可行解 否 是 否 是 否 是 目标函数值 211 0 0 55111 0 0 3610 2 0 210 0 2 243 5 5 0 0 1 1 5 211最优解x,0,,0。
55T
(a)
(1) 图解法
x2 4 3 2 1 0 1 2 3 x1
最优解即为3x14x29353的解x1,,最大值z
25x2x8221(2)单纯形法
首先在各约束条件上添加松弛变量,将问题转化为标准形式 max z10x15x20x30x43x4x2x39s.t. 15x12x2x48
则P3,P4组成一个基。令x1x20
得基可行解x0,0,9,8,由此列出初始单纯形表 cj cB 基 b 10 5 0 0 x1 x2 x3 x4 0 x3 9 0 x4 8 3 4 1 0 [5] 2 0 1 10 5 0 0 cjzj 12。min, cj cB 基 b 89538 510 5 0 0 x1 x2 x3 x4 3140 1 55 210 x3 5810 x1 52 11 0 55cjzj 0 1 0 2
20,min新的单纯形表为 cj cB 基 b 35 x2 22183, 142210 5 0 0 x1 x2 x3 x4 530 1 141410 x1 1 1 21 0 770 0 525 14143235 2cjzj
1,20,表明已找到问题最优解x11, x2 , x30 , x40。最大值 z*(b) (1) 图解法 6x12x224x2 12 9 x1x256 3 0 3 6 9 x1
最优解即为6x12x2241773的解x,,最大值z
222x1x25(2) 单纯形法
首先在各约束条件上添加松弛变量,将问题转化为标准形式
max z2x1x20x30x40x55x2x315s.t. 6x12x2x424xxx5125
则P3,P4,P5组成一个基。令x1x20
得基可行解x0,0,15,24,5,由此列出初始单纯形表 cj cB 基 b 2 1 0 0 0 x1 x2 x3 x4 x5 0 5 1 0 0 [6] 2 0 1 0 1 1 0 0 1 2 1 0 0 0 0 x3 15 0 x4 24 0 x5 5 cjzj
12。min,cj cB 基 b 245,4 612 1 0 0 0 x1 x2 x3 x4 x5 0 5 1 0 0 1 0 x3 15 2 x4 4 0 x5 1 cjzj 11 0 0 360 0 1 6321110 0 0 333315,24,
22520,min新的单纯形表为 cj cB 基 b 2 1 0 0 0 x1 x2 x3 x4 x5 0 0 1 15x0 3 2 515 427x2 4 2 30 x5 2cjzj 11 42310 1 0 24110 0 0 421 0 0 1,20,表明已找到问题最优解x11, x2值 z
*715,x40,x50。最大 ,x32217 2''''''x2 x20,x20(a) 在约束条件中添加松弛变量或剩余变量,且令x2x2,
'x3x3, z'z
该问题转化为
''''max z'3x1x2x22x30x40x5''''2x13x23x24x3x412''''4x1x2x22x3x58s.t. ''''3x1x2x23x36''''x1,x2,x2,x3,x4,x50
其约束系数矩阵为
233410A411201
311300在A中人为地添加两列单位向量P7,P8 2334100041120110 31130001''''x22x30x40x5Mx6Mx7 令max z'3x1x2得初始单纯形表
cj cB 基 b 3 1 1 2 0 0 M M ''''x1 x2 x2 x3 x4 x5 x6 x7 0 x4 12 M x6 8 M x7 6 2 3 3 4 1 0 0 0 4 1 1 -2 0 -1 1 0 3 1 1 -3 0 0 0 1 cjzj 37M 1 1 25M 0 M 0 0 '''(b) 在约束条件中添加松弛变量或剩余变量,且令x3x3x3 x30,x30该问题转化为
'''max z'3x15x2x3x30x40x5'''x12x2x3x3x46'''2x1x23x33x3x516s.t. '''xx5x5x101233'''x1,x2,x3,x3,x4,x50''', z'z
其约束系数矩阵为
121110
A213301115500在A中人为地添加两列单位向量P7,P8
1211101021330100 11550001令max z'3x15x2x3x30x40x5Mx6Mx7 得初始单纯形表 cj cB 基 b '''3 -5 1 -1 0 0 -M M x3 x4 x5 x6 x7 x1 x2 x3M x6 6 0 x5 16 M x7 10 cjzj 1 2 1 -1 -1 0 1 0 2 1 3 -3 0 1 0 0 1 1 5 -5 0 0 0 1 32M 53M 1+6M -1-6M -M 0 0 0
(a)解1:大M法
在上述线性规划问题中分别减去剩余变量x4,x6,x8,再加上人工变量x5,x7,x9,得
maxz2x1x22x30x4Mx50x6Mx70x8Mx9
x1x2x3x4x562x1x3x6x72 s,t,2x2x3x8x90x1,x2,x3,x4,x5,x6,x7,x8,x90
其中M是一个任意大的正数。据此可列出单纯形表 cj 2x11x22x30x4Mx50x6Mx70x8M x9 i cb基 b Mx56Mx72 Mx90cjzj 1202M10[2]111100M1000010M0100001M00 10 6 0 3M12MMx56Mx72 1x20cjzj 1202M00103/2[1]1/25M322100M1000010M01001/21/200 42 1/21/2M113M 2222Mx532x32 1x21cjzj [4]210010101001003/23/2111/21/21/21/200 1/21/23/4 4M500M03M35M3M113M 22222x13/42x37/2 1x27/4cjzj 1000001001001/41/21/41/41/21/43/83/81/41/81/41/81/81/41/4 3/89M 8 1/83/839M885/4M53/84 由单纯形表计算结果可以看出,40且ai40(i1,2,3),所以该线性规划问题有无界解 解2:两阶段法。
现在上述线性规划问题的约束条件中分别减去剩余变量x4,x6,x8,再加上人工变量
x5,x7,x9,得第一阶段的数学模型
据此可列出单纯形表 cj 0x10x20x30x41x50x61x70x81 x9 i cb基 b 1x561x72 1x90cjzj 120110[2]311111001100001010100001100 10 6 0 1x561x72 0x20cjzj 120100103/2[1]1/25/210011000010101001/21/20124 03 21/21/22 10x53x21cjzj [4]210010101001003/23/2111/21/21/21/200 1/21/23/4 0x32 000010101 2x13/42x37/2 1x27/4cjzj 1000*001001001/41/21/40T1/41/21/413/83/81/41/801/81/81/41/4 3/81 1/41/83/81*0 第一阶段求得的最优解X(,,,0,0,0,0,0,0),目标函数的最优值0。
*因人工变量x5x7x90,所以X(,377442377,,0,0,0,0,0,0)T是原线性规划问题的基可442行解。于是可以进行第二阶段运算。将第一阶段的最终表中的人工变量取消,并填入原问题的目标函数的系数,进行第二阶段的运算,见下表。 cjzj 2x11x22x30x40x60x8 0 i cb基 b 2x13/42x37/2 1x27/41000010101/41/21/43/81/41/81/4 1/83/8cjzj 0005/43/89/8 由表中计算结果可以看出,40且ai40(i1,2,3),所以原线性规划问题有无界解。
(b)解1:大M法
在上述线性规划问题中分别减去剩余变量x4,x6,x8,再加上人工变量x5,x7,x9,得
minz2x13x2x30x40x5Mx6Mx7
x14x22x3x4x683x12x2x5x76 s,t,x1,x2,x3,x4,x5,x6,x7,x8,x90
其中M是一个任意大的正数。据此可列出单纯形表 cj 2x11x22x30x4Mx50x6Mx7 0M i cb基 b MMx6813[4]22010011001x76 cjzj 23 24M1/4[5/2]36M1012M1/21M1/41/2M0101/41/20 0 13x2Mx722 84/5 cjzj 55M4201010M1231M423/101/5M3M3240 3x22x19/54/5 3/52/51/103/101/10 2/51/52/5cjzj 000*1/21/2M1/2M1/2 由单纯形表计算结果可以看出,最优解X(,49,0,0,0,0,0)T,目标函数的最优解值55z*24937。X存在非基变量检验数30,故该线性规划问题有无穷多最优解。 55解2:两阶段法。
现在上述线性规划问题的约束条件中分别减去剩余变量x4,x5,再加上人工变量x6,x7,得第一阶段的数学模型minx6x7
x14x22x3x4x683x12x2x5x76 s,t,x1,x2,x3,x4,x5,x6,x7,x8,x90据此可列出单纯形表 cj 0x10x20x30x40x51x61 x7 i cb基 b 1x61x786 13[4]22010011001 23 cjzj 41/4[5/2]61021/2111/41/210101/41/20 0 10x21x722 84/5 cjzj 5/2 01013/52/51/213/23/101/50 1/10 2/500x2x19/54/5013/101/101/52/5cjzj 0*000011 49,0,0,0,0,0)T,目标函数的最优值*0。 5549T因人工变量x6x70,所以(,,0,0,0,0,0)是原线性规划问题的基可行解。于是可
55第一阶段求得的最优解X(,以进行第二阶段运算。将第一阶段的最终表中的人工变量取消,并填入原问题的目标函数的系数,进行第二阶段的运算,见下表。
cjzj 2x13x21x30x40 x5 i cb基 b 3x22x19/54/5 01103/52/53/101/51/102/5 cjzj 000*1/21/2 由单纯形表计算结果可以看出,最优解X(,49,0,0,0,0,0)T,目标函数的最优解值55z*249由于存在非基变量检验数30,故该线性规划问题有无穷多最优37。55解。
表1-23 x1 x2 x3 x4 x5 x4 6 x5 1 2 4 -2 1 0 1 3 2 0 1 3 1 2 0 0 cjzj 表1-24 x1 x2 x3 x4 x5 x1 3 x5 1 1 2 1 12 0 0 5 1 12 1 0 7 5 32 0 cjzj
5 x2 83 0 x5 143 0 x6 293 3 5 4 0 0 0 x1 x2 x3 x4 x5 x6 23 1 0 13 0 0 43 0 5 23 1 0 53 0 4 23 0 1 cjzj 13 0 4 53 0 0 5 x2 83 4 x3 1415 x1 x2 x3 x4 x5 x6 23 1 0 13 0 0 415 0 1 215 15 0 0 x6 8915 4115 0 0 215 45 1 cjzj 1115 0 0 1715 45 0 5 x2 5041 4 x3 6241 3 x1 8941 x1 x2 x3 x4 x5 x6 0 1 0 1541 841 1041 0 0 1 641 541 441 1 0 0 241 1241 1541 cjzj 0 0 0 4541 2441 1141 最后一个表为所求。
习题二 P76 写出对偶问题 (a)
min z2x12x24x3max w2y13y25y3y12y2y32x13x24x322xx3xy3 对偶问题为:3yy4y2
1234123s.t. s.t. 4y13y23y34x14x23x35y10,y20,y3无约束x1,x20,x3无约束(b)
max z5x16x23x3min w5y13y28y3y1y24y35x12x22x35x5xx3 对偶问题为: 2y5y7y6 123123s.t. s.t. 4x17x23x382y1y23y33x1无约束,x20,x30y1无约束,y20,y30
(a)错误。原问题存在可行解,对偶问题可能存在可行解,也可能无可行解。
(b)错误。线性规划的对偶问题无可行解,则原问题可能无可行解,也可能为无界解。 (c)错误。 (d)正确。
对偶单纯形法 (a)
min z4x112x218x3x1 3x33 s.t. 2x22x35x,x,x0123解:先将问题改写为求目标函数极大化,并化为标准形式
max z'4x112x218x30x40x5x1 3x3x4 3s.t. 2x22x3 x55x0i1,,5i
列单纯形表,用对偶单纯形法求解,步骤如下 cj cB 基 b 4 12 18 0 0 x1 x2 x3 x4 x5 0 x4 3 0 x5 5 1 0 3 1 0 0 2 2 0 1 cjzj 0 x4 3 12 x2 52 4 12 18 0 0 1 0 3 1 0 0 1 1 0 12 cjzj 4 0 6 0 6 18 x3 1 12 x2 32 11 0 1 033 111 1 0 332 2 0 0 2 6 Tcjzj 3最优解为x0,1,, 目标值z39。
2(b)
min z5x12x24x33x1x2 2x34 s.t. 6x13x25x310x,x,x0123解:先将问题改写为求目标函数极大化,并化为标准形式
max z'5x12x24x30x40x53x1 x2 2x3x4 4 s.t. 6x13x25x3 x510x0i1,,5i列单纯形表,用对偶单纯形法求解 cj 5 2 4 0 0 cB 基 b x1 x2 x3 x4 x5 0 x4 4 0 x5 10 3 1 2 1 0 6 3 5 0 1 cjzj 5 2 4 0 0 23 0 x4 111 0 1 33 102 x2 3 cjzj512 1 0 33 221 0 0 33 4 x3 2 2 x2 0 3 0 1 3 1 3 1 0 5 2 1 0 0 2 0 cjzj 最优解为x0,0,2T, 目标值z8。 将该问题化为标准形式:
max z2x1x2x30x40x5x1x2x3x46s.t. x12x2x54x0i1,5i
用单纯形表求解 cj cB 基 b 2 1 1 0 0 x1 x2 x3 x4 x5 0 x4 6 0 x5 4 [1 ] 1 1 1 0 1 2 0 0 1 2 1 1 0 0 cjzj 6 cB 基 b x1 x2 x3 x4 x5 2 x1 6 0 x5 10 1 1 1 1 0 0 3 1 1 1 0 3 -1 2 0 cjzj 由于j0,所以已找到最优解X*6,0,0,0,10,目标函数值z*12 (a) 令目标函数
max z(21)x1(-1+2)x2(1+3)x3
(1)令230,将1反映到最终单纯形表中 cj cB 基 b 21 1 1 0 0 x1 x2 x3 x4 x5 21 x4 6 0 x5 10 1 1 1 1 0 0 3 1 1 1 0 -3-1 - 1 -1 2-1 0 cjzj 表中解为最优的条件:-3-10,- 1 -10,2-10,从而11 (2)令130,将2反映到最终单纯形表中 cj cB 基 b 2 12 1 0 0 x1 x2 x3 x4 x5 2 x1 6 0 x5 10 1 1 1 1 0 0 3 1 1 1 0 2-3 - 1 2 0 cjzj 表中解为最优的条件:2-3 0, 从而23 (3) 令120,将3反映到最终单纯形表中 cj cB 基 b 2 1 13 0 0 x1 x2 x3 x4 x5 2 x1 6 0 x5 10 1 1 1 1 0 0 3 1 1 1 0 -3 3 - 1 2 0 cjzj 表中解为最优的条件:3-10, 从而31
(b) 令线性规划问题为
max z2x1x2x3x1x2x364s.t. x12x245x0i1,3 i(1)先分析的变化
1011bB1b110
1使问题最优基不变的条件是bb100,从而16
1(2)同理有610,从而210 1026(c) 由于x(6,0,0,0,10)代入x12x362,所以将约束条件减去剩余变量后的方程x12x3x62直接反映到最终单纯形表中 cj cB 基 b 2 -1 1 0 0 0 x1 x2 x3 x4 x5 x6 1 1 1 1 0 0 0 3 1 1 1 0 1 0 -2 0 0 1 2 x1 6 0 x5 10 0 x6 -2 cjzj 0 -3 -1 -2 0 0 对表中系数矩阵进行初等变换,得 cj cB 基 b 2 -1 1 0 0 0 x1 x2 x3 x4 x5 x6 2 x1 6 0 x5 10 1 1 1 1 0 0 0 3 1 1 1 0 0 x6 -8 cjzj 0 -1 [-3] -1 0 1 0 -3 -1 -2 0 0 cj cB 基 b 2 -1 1 0 0 0 x1 x2 x3 x4 x5 x6 1 23 0 23 0
2 x1 10 31 33 0 20 x5 22 0 833 1
1 30 x6 8 310 1 1 1 0 33cjzj 38510 0 0 333因此增加约束条件后,新的最优解为
x11082228,x3,x5,最优值为 3333
(a) 线性规划问题
max z3x1x24x36x13x25x345s.t. 3x14x25x330x1,x2,x30
单纯形法求解 cB 基 b x1 x2 x3 x4 x5 0 x4 45 0 x5 30 6 3 5 1 0 3 4 5 0 1 cjzj 0 x4 15 4 x3 6 3 1 4 0 0 3 1 0 1 1 341 1 0 555cjzj 3 x1 5 4 x3 3 3114 0 0 5551111 0 333120 1 1 55cjzj 130 2 0 55最优解为x1,x2,x35,0,3 ,目标值z27。 (a) 设产品A的利润为3,线性规划问题变为
max z3x1x24x36x13x25x345s.t. 3x14x25x330x1,x2,x30
单纯形法求解 基 b x4 45 x5 30 x1 x2 x3 x4 x5 6 3 5 1 0 3 4 5 0 1 cjzj x4 15 x3 6 3 1 4 0 0 3 1 0 1 1 341 1 0 5553114 0 0 5551111 0 333120 1 1 55cjzj x1 5 x3 3 cjzj 130 2 0 35353113139,,都小于等于0,解得。 3535355
为保持最优计划不变,应使2(b) 线性规划问题变为
max z3x1x24x33x46x13x25x38x445s.t. 3x14x25x32x430x1,x2,x3,x40
单纯形法求解
cB 基 b x1 x2 x3 x4 x5 x6 0 x5 45 0 x6 30 6 3 5 8 1 0 3 4 5 2 0 1 cjzj 0 x4 15 4 x3 6 3 1 4 3 0 0 3 1 0 6 1 1 3421 1 0 555531174 0 0 55551111 0 2 3334120 1 1 555cjzj 3 x1 5 4 x3 3 cjzj 50 x4 21130 2 0 5551111 0 1 266621318 1 0 5151515129717 0 0 101530304 x3 5 cjzj
此时最优解为x1,x2,x30,0,5,目标值z20,小于原最优值,因此该种产品不值得生产。
(c) 设购买材料数量为y,则规划问题变为
max z3x1x24x30.4y6x13x25x345s.t. 3x14x25x3y30x1,x2,x3,y0
单纯形法求解 cB 基 b x1 x2 x3 y x4 x5 0 x5 45 0 x6 30 6 3 5 0 1 0 3 4 5 1 0 1 cjzj 0 x4 15 23 1 4 0 0 53 1 0 1 1 1 4 x3 6 3411 1 0 555531124 0 0 555511111 0 3333cjzj 3 x1 5 4 x3 3 2120 1 1 555cjzj 1130 2 0 5550 y 15 4 x3 9 3 1 0 1 1 1 631 1 0 0 5553922 0 0 5555cjzj 此时最优解为x1,x2,x3,y0,0,9,15,目标值z30,大于原最优值,因此应该购进原材料扩大生产,以购材料15单位为宜。
第三章
表
产地 销地 A1 A2 A3 A4 销量 B1 B2 B3 B4 9 8 12 13 10 10 12 14 8 9 11 12 10 10 11 12 6 14 35 5 产量 18 24 6 12
用vogel法求解得
B A A1 A2 A3 A4 14 24 6 0 11 4 B1 B2 B3 B4
用位势法检验,把上表中有数字的地方换成运价
B A A1 A2 A3 A4 Vj B1 B2 B3 B4 8 13 12 8 11 11 12 1 0 4 5 Ui 8 8 7 7
令v1=1
则u1+v2=8 所以u3=7 u1+v4=13 v3=4 u2+v3=12 u4=7 u3+v1=8 v5=8 u3+v3=11 u2=8 u4+v3=11 v2=0 u4+v4=12
得检验数表
B A A1 A2 A3 A4 0 0 1 2 1 2 0 2 3 B1 B2 B3 B4 表中所有的数字均大于等于零,故所求方案为最优方案 解:(a)用运价代替表中有数字的地方,求出位势和检验数
B A A1 A2 A3 Vj 1 11 12 k 9 2 1 k-11 -2 k-1 12-k 11 1 B1 B2 B3 B4 Ui 令v1=1则u1+v2=1 故v3=-2
u1+v4=11 u2=11 u2+v1=12 v2=k-11 u2+v2=k u1=12-k u2+v3=9 v4=k-1 u3+v1=2 u3=1
得检验数表
B B1 B2 B3 B4 A A1 A2 A3 k-3 10-k 30-k 24-k 15 18-k 令表中所有的检验数均大于等于零,得3≤k≤10
(b)由表和表计算出位势和检验数,令C24=M
位势表 B B1 B2 B3 B4 Ui A A1 A2 A3 Vj
检验数表
B A A1 A2 A3 4 17 M-17 17 17 11 B1 B2 B3 B4 1 11 12 7 9 2 1 -4 -2 6 5 11 1
当存在某非基变量的检验数大于等于零的时候有无穷多最优解 则M-17=0 所以M=17 故运价C24=17 销地 A1 A2 A3 供应量 产地 B1 500 540 580 2 B1’ 570 610 650 3 B2 M 600 640 4 B2’ M 670 710 2 B3 M M 550 1 B3’ M M 620 3 S 40 80 120 2 销量
3 3 4 由于产大于销,设有一个理想的销地A4,则 销地 产地 B1 B1’ B2 B2’ B3 B3’ S 销量
A1 A2 A3 A4 500 540 580 0 570 610 650 0 M 600 640 0 M 670 710 0 M M 550 0 M M 620 0 40 80 120 0 3 3 4 7 供应量 2 3 4 2 1 3 2 第四章
略
第五章
解:
Min z=P1 d1- +P2(d2-+d2+) . 11x1+3x2≥25
100 x1+50x2 + d1- - d1+ =1900 10x1+16x2 + d2- - d2+ =200 x1,x2, d1- ,d1+ , d2- ,d2+ ≥0
解:
目标规划模型为
Min z = P1 d1- +P2(d2-+d2+) +P2(d3-+d3+) +P3(d4-+d4+) +P3d5+ . 300 x1+450x2 + d1- - d1+ =10000 x1+ d2- - d2+ =10
x2 + d3- - d3+ =15
4x1+6x2 + d4- - d4+ =150
3x1+2x2 + d5- - d5+ =75
x1,x2, d1- ,di+ , di- ≥0 ,i=1,…,5
因篇幅问题不能全部显示,请点此查看更多更全内容