【运筹学】线性规划数学模型 ( 单纯形法 | 第二次迭代 | 方程组同解变换 | 生成新单纯形表 | 计算检验数 | 最优解判定 | 线性规划解个数分析 ) 后续博客 , 在上一篇博客中进行了 第二次迭代 , 使用中心元变换得到新的系数矩阵 , 计算检验数 , 验证最优解 , 计算入基变量 , 出基变量 , 本篇博客开始进行第三次迭代 ;
一、第二次迭代 : 进行矩阵变换
在上一篇博客中 , 第一次迭代后 , 找到 入基变量 x1 , 出基变量 x4 , 使用 x1 替换基变量中的 x4 位置 ;
新的单纯形表为 :
cj |
cj |
|
1 |
2 |
1 |
0 |
0 |
|
CB 基变量系数 (目标函数) |
基变量 |
常数 b
|
x1 |
x2 |
x3 |
x4 |
x5 |
θi |
0 ( 目标函数 x4 系数 c4 ) |
x4 |
15 |
2 |
−3 |
2 |
1 |
0 |
− (θ4) |
0 ( 目标函数 x5 系数 c5) |
x5 |
20 |
31 |
1 |
5 |
0 |
1 |
20 ( θ5 ) |
σj ( 检验数 ) |
|
|
1 ( σ1 ) |
2 ( σ2 ) |
1 ( σ3 ) |
0 |
0 |
|
第一次迭代 |
– |
– |
– |
– |
– |
– |
– |
– |
0 ( 目标函数 x4 系数 c4 ) |
x4 |
75 |
3 |
0 |
17 |
1 |
3 |
25 ( θ4 ) |
2 ( 目标函数 x2 系数 c2) |
x2 |
20 |
31 |
1 |
5 |
0 |
1 |
60 (θ2) |
σj ( 检验数 ) |
|
|
31 ( σ1 ) |
0 |
−9 ( σ3 ) |
0 |
−2 ( σ5 ) |
|
第二次迭代 |
– |
– |
– |
– |
– |
– |
– |
– |
1 ( 目标函数 x1 系数 c1 ) |
x1 |
? |
? |
? |
? |
? |
? |
? ( θ1 ) |
2 ( 目标函数 x2 系数 c2) |
x2 |
? |
? |
? |
? |
? |
? |
? (θ2) |
σj ( 检验数 ) |
|
|
0 |
0 |
? ( σ3 ) |
? ( σ4 ) |
? ( σ5 ) |
|
中心元 : 入基变量 x1 所在列 , 与出基变量 x4 所在行 , 相交的位置叫做中心元 ;
-
中心元系数 : 转换成 1 ;
-
中心元所在列另一个系数 : 转换成 0 ;
当前的线性规划标准形式等式方程组 : ⎩⎪⎪⎨⎪⎪⎧3x1+0x2+17x3+x4+3x5=7531x1+x2+5x3+0x4+x5=20
中心元转换为 1 : 将 3x1+0x2+17x3+x4+3x5=75 中的 x1 系数变成 1 , 只需要将方程等式两边都乘以 31 即可 ;
(3x1+0x2+17x3+x4+3x5)×31x1+0x2+317x3+31x4+x5==75×3125
与中心元同一列的 x1 的系数转换为 0 : 将 31x1+x2+5x3+0x4+x5=20 中的 x1 系数转为 0 :
- 将 x1+0x2+317x3+31x4+x5=25 方程 左右变量乘以 −31 ;
- 将上个步骤得到的等式与 31x1+x2+5x3+0x4+x5=20 相加 , 即可使 x1 系数为 0 ;
(x1+0x2+317x3+31x4+x5)×−31+(31x1+x2+5x3+0x4+x5)=25×−31+20(−x1+0x2−917x3−91x4−31x5)+(31x1+x2+5x3+0x4+x5)=3350x1+x2+928x3−91x4+32x5=335
新的单纯形表为 :
cj |
cj |
|
1 |
2 |
1 |
0 |
0 |
|
CB 基变量系数 (目标函数) |
基变量 |
常数 b
|
x1 |
x2 |
x3 |
x4 |
x5 |
θi |
0 ( 目标函数 x4 系数 c4 ) |
x4 |
15 |
2 |
−3 |
2 |
1 |
0 |
− (θ4) |
0 ( 目标函数 x5 系数 c5) |
x5 |
20 |
31 |
1 |
5 |
0 |
1 |
20 ( θ5 ) |
σj ( 检验数 ) |
|
|
1 ( σ1 ) |
2 ( σ2 ) |
1 ( σ3 ) |
0 |
0 |
|
第一次迭代 |
– |
– |
– |
– |
– |
– |
– |
– |
0 ( 目标函数 x4 系数 c4 ) |
x4 |
75 |
3 |
0 |
17 |
1 |
3 |
25 ( θ4 ) |
2 ( 目标函数 x2 系数 c2) |
x2 |
20 |
31 |
1 |
5 |
0 |
1 |
60 (θ2) |
σj ( 检验数 ) |
|
|
31 ( σ1 ) |
0 |
−9 ( σ3 ) |
0 |
−2 ( σ5 ) |
|
第二次迭代 |
– |
– |
– |
– |
– |
– |
– |
– |
1 ( 目标函数 x1 系数 c1 ) |
x1 |
25 |
1 |
0 |
317 |
31 |
1 |
? ( θ1 ) |
2 ( 目标函数 x2 系数 c2) |
x2 |
335 |
0 |
1 |
928 |
−91 |
32 |
? (θ2) |
σj ( 检验数 ) |
|
|
0 |
0 |
? ( σ3 ) |
? ( σ4 ) |
? ( σ5 ) |
|
二、第二次迭代 : 计算检验数
1 . 计算非基变量 x3 的检验数 σ3 :
σ3=1−(12)×⎝⎜⎜⎛317928⎠⎟⎟⎞=1−(1×317+2×928)=99−17×3+28×2=−998
2 . 计算非基变量 x4 的检验数 σ4 :
σ4=0−(12)×⎝⎜⎜⎛31−91⎠⎟⎟⎞=0−(1×31−2×91)=−91
3 . 计算非基变量 x5 的检验数 σ5 :
σ5=0−(12)×⎝⎜⎛132⎠⎟⎞=0−(1×1+2×32)=−37
新的单纯形表为 :
cj |
cj |
|
1 |
2 |
1 |
0 |
0 |
|
CB 基变量系数 (目标函数) |
基变量 |
常数 b
|
x1 |
x2 |
x3 |
x4 |
x5 |
θi |
0 ( 目标函数 x4 系数 c4 ) |
x4 |
15 |
2 |
−3 |
2 |
1 |
0 |
− (θ4) |
0 ( 目标函数 x5 系数 c5) |
x5 |
20 |
31 |
1 |
5 |
0 |
1 |
20 ( θ5 ) |
σj ( 检验数 ) |
|
|
1 ( σ1 ) |
2 ( σ2 ) |
1 ( σ3 ) |
0 |
0 |
|
第一次迭代 |
– |
– |
– |
– |
– |
– |
– |
– |
0 ( 目标函数 x4 系数 c4 ) |
x4 |
75 |
3 |
0 |
17 |
1 |
3 |
25 ( θ4 ) |
2 ( 目标函数 x2 系数 c2) |
x2 |
20 |
31 |
1 |
5 |
0 |
1 |
60 (θ2) |
σj ( 检验数 ) |
|
|
31 ( σ1 ) |
0 |
−9 ( σ3 ) |
0 |
−2 ( σ5 ) |
|
第二次迭代 |
– |
– |
– |
– |
– |
– |
– |
– |
1 ( 目标函数 x1 系数 c1 ) |
x1 |
25 |
1 |
0 |
317 |
31 |
1 |
? ( θ1 ) |
2 ( 目标函数 x2 系数 c2) |
x2 |
335 |
0 |
1 |
928 |
−91 |
32 |
? (θ2) |
σj ( 检验数 ) |
|
|
0 |
0 |
−998 ( σ3 ) |
−91 ( σ4 ) |
−37 ( σ5 ) |
|
三、第二次迭代 : 最优解判定
上述三个检验数 , ⎩⎪⎪⎪⎪⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎧σ3=−998σ4=−91σ5=−37 , 三个检验数都小于等于 0 , 该基可行解是最优解 ;