【运筹学】线性规划 单纯形法 阶段总结 ( 初始基可行解 | 判定最优解 | 迭代 | 得到最优解 | 全流程详细解析 ) ★

单纯形法 参考博客 :

1 . 查找初始基可行解 :

• 【运筹学】线性规划数学模型 ( 单纯形法原理 | 单纯形法流程 | 查找初始基可行解 )

2 . 最优解判定 :

• 【运筹学】线性规划数学模型 ( 单纯形法 | 最优解判定原则 | 可行解表示 | 目标函数推导 | 目标函数最大值分析 )
• 【运筹学】线性规划数学模型 ( 单纯形法 | 最优解判定原则 | 单纯形表 | 系数计算方法 | 根据系数是否小于等于 0 判定最优解 )
• 【运筹学】线性规划数学模型 ( 单纯形法 | 最优解判定原则 | 线性规划求解示例 )

3 . 迭代原则 :

• 【运筹学】线性规划数学模型 ( 单纯形法 | 迭代原则 | 入基 | 出基 | 线性规划求解示例 )
• 【运筹学】线性规划数学模型 ( 单纯形法 | 第一次迭代 | 方程组同解变换 | 计算新单纯形表 | 计算检验数 | 入基变量选择 | 出基变量选择 )
• 【运筹学】线性规划数学模型 ( 单纯形法 | 第二次迭代 | 方程组同解变换 | 生成新单纯形表 | 计算检验数 | 最优解判定 | 线性规划解个数分析 )

4 . 单纯形法阶段总结 :

• 【运筹学】线性规划 单纯形法 阶段总结 ( 初始基可行解 | 判定最优解 | 迭代 | 得到最优解 | 全流程详细解析 ) ★ 推荐

一、线性规划示例

---

使用单纯形法求解线性规划最优解 :

$$\begin{array}{l}maxZ=3x_1+4x_2\\ \\ \{\begin{array}{ll}2x_1+x_2\le 40& \\ & \\ x_1+3x_2\le 30& \\ & \\ x_j\ge 0& (j=1,2)\end{array}\end{array}$$

二、转化标准形式

---

首先将该线性规划转为标准形式 :

参考 【运筹学】线性规划数学模型标准形式 ( 标准形式 | 目标函数转化 | 决策变量转化 | 约束方程转化 | 固定转化顺序 | 标准形式转化实例 ) 线性规划 普通形式 -> 标准形式 转化顺序说明 博客 , 先处理变量约束 , 再将不等式转为等式 , 最后更新目标函数 ;

① 变量约束 : 首先查看变量约束 , 两个变量都是 $\ge 0$ 的 , 符合线性规划标准形式要求 ;

② 不等式转换 : 两个等式都是 小于等于不等式 , 需要 在不等式左侧加入松弛变量 , 将其转为等式 ;

• $2x_1+x_2\le 40$ , 左侧加入松弛变量 $x_3$ , 变为 $2x_1+x_2+x_3=40$
• $x_1+3x_2\le 30$ , 左侧加入松弛变量 $x_4$ , 变为 $x_1+3x_2+x_4=30$

③ 更新目标函数 : 将 $x_3$ 和 $x_4$ 加入到目标函数中 , 得到新的目标函数 $maxZ=3x_1+4x_2+0x_3+0x_4$ ;

此时线性规划标准形式为 :

$$\begin{array}{l}maxZ=3x_1+4x_2+0x_3+0x_4\\ \\ \{\begin{array}{ll}2x_1+x_2+x_3+0x_4=40& \\ & \\ x_1+3x_2+0x_3+x_4=30& \\ & \\ x_j\ge 0& (j=1,2,3,4)\end{array}\end{array}$$

三、查找初始基可行解

---

找基矩阵 :

上述线性规划标准形式的系数矩阵为 $\left[\begin{array}{cccccc}& 2& 1& 1& 0& \\ & & & & & \\ & 1& 3& 0& 1& \end{array}\right]$ , 其中子矩阵中有 $\left[\begin{array}{cccc}& 1& 0& \\ & & & \\ & 0& 1& \end{array}\right]$ 单位阵 $I$ ;

使用该单位阵 $I$ 作为基矩阵 , 单位阵肯定是可逆的 , 其对应的基解 , 解出后的值就是右侧的常数值 , 肯定大于等于 $0$ , 是基可行解 ;

列出单纯形表 :

$c_j$	$c_j$		$3$	$4$	$0$	$0$
$C_B$ 基变量系数 (目标函数)	基变量	常数 $b$	$x_1$	$x_2$	$x_3$	$x_4$	${\theta }_i$
$0$ ( 目标函数 $x_3$ 系数 $c_3$ )	$x_3$	$40$	$2$	$1$	$1$	$0$
$0$ ( 目标函数 $x_4$ 系数 $c_4$)	$x_4$	$30$	$1$	$3$	$0$	$1$
		${\sigma }_j$	$3$	$4$	$0$	$0$

基变量是 $x_3$ 和 $x_4$ , 基变量在约束条件中的系数矩阵 $\left[\begin{array}{cccc}& 1& 0& \\ & & & \\ & 0& 1& \end{array}\right]$ 就是基矩阵 , 这是个单位阵 ;

基变量是 $x_3$ 和 $x_4$ 在目标函数中的系数是 $\left(\begin{array}{c}00\end{array}\right)$ ;

此时的基解是 $\left(\begin{array}{c}0\\ 0\\ 40\\ 30\end{array}\right)$ , 该解是初始解 , 下面判定该解是否是最优解 ;

四、初始基可行解的最优解判定

---

使用 检验数矩阵 $(C_N^T-C_B^T{B}^{-1}N)$ 判断上述解 , 是否是最优解 , 该矩阵计算结果中所有的数 , 都是检验数 $\sigma$ , 如果 所有的数都小于等于 $0$ , 说明该解就是最优解 ;

这里只求非基变量的检验数 , 即 $x_1$ , $x_2$ 的检验数 ;

列出目标函数非基变量系数 $(C_N^T-C_B^T{B}^{-1}N)$ 矩阵 :

• 非基变量在目标函数中的系数矩阵 : $C_N^T=\left(\begin{array}{c}34\end{array}\right)$

• 基变量在目标函数中的叙述矩阵 : $C_B^T=\left(\begin{array}{c}00\end{array}\right)$

• $B^{-1}N$ 是系数矩阵中经过矩阵变换后 , 基变量系数是单位阵 $I$ , 非基变量系数是 $B^{-1}N$ : $B^{-1}N=\left[\begin{array}{cccc}& 2& 1& \\ & & & \\ & 1& 3& \end{array}\right]$

$$(C_N^T-C_B^T{B}^{-1}N)=C_N^T=\left(\begin{array}{c}34\end{array}\right)-\left(\begin{array}{c}00\end{array}\right)\times \left[\begin{array}{cccc}& 2& 1& \\ & & & \\ & 1& 3& \end{array}\right]$$

$$=\left(\begin{array}{c}{\sigma }_1{\sigma }_2\end{array}\right)$$

其中 ${\sigma }_1$ 是目标函数中 $x_1$ 的系数 , ${\sigma }_2$ 是目标函数中的 $x_2$ 的系数 ;

如果上述两个系数都小于等于 $0$ , 那么当 非基变量 $X_N=\left(\begin{array}{c}{x}_1\\ x_2\end{array}\right)$ 取值为 $0$ 时 , 目标函数取值最大 ;

系数的计算公式为 : ${\sigma }_j=c_j-\sum c_i{a}_{ij}$ , 其中 $c_j$ 对应的是非基变量在目标函数系数 , $c_i$ 是基变量在目标函数中的系数 , $a_{ij}$ 是 $B^{-1}N$ 中的矩阵向量 , 代表一列 ;

${\sigma }_1=c_1-(c_3{a}_{11}+c_4{a}_{12})$

${\sigma }_1=3-(0\times 2)-(0\times 1)=3$ , 是从下面的单纯形表中的如下位置提取的数值 ;

${\sigma }_2=c_2-(c_3{a}_{21}+c_4{a}_{22})$

${\sigma }_2=4-(0\times 1)-(0\times 3)=4$ , 是从下面的单纯形表中的如下位置提取的数值 ;

如果这两个系数 , 如果都小于等于 $0$ , 该 基可行解$\left(\begin{array}{c}0\\ 0\\ 40\\ 30\end{array}\right)$ 才是最优解 , 这两个系数都大于 $0$ , 因此不是最优解 ;

五、第一次迭代 : 入基与出基变量选择

---

入基变量选择 : 具体哪个变量入基 , 是由检验数决定的 , 检验数 ${\sigma }_j$ 较大的入基 ; $x_2$ 的检验数 ${\sigma }_2$ 是 $4$ , 大于 ${\sigma }_1=3$ , 因此这里选择 $x_2$ 作为入基变量 ;

出基变量选择 : 系数矩阵中 , 常数列 $b=\left(\begin{array}{c}40\\ 30\end{array}\right)$ , 分别除以除以入基变量大于 $0$ 的系数列 $\left(\begin{array}{c}1\\ 3\end{array}\right)$ , 得出结果是 $\left(\begin{array}{c}40\\ 10\end{array}\right)$ , 然后选择一个最小值 $10$ , 查看该最小值对应的变量是 $x_4$ , 选择该变量作为出基变量 ;

这里将出基变量与入基变量选择好了 , $x_2$ 的检验数较大 , 选择 $x_2$ 作为入基变量 , $x_4$ 的 ${\theta }_4$ 较小 , 选择 $x_4$ 作为出基变量 ;

入基出基操作完成后 , 基变量变成了 $x_3,x_2$ ;

六、第一次迭代 : 方程组同解变换

---

方程组做同解变换 :

线性规划原始方程组为 $\{\begin{array}{l}2x_1+x_2+x_3+0x_4=40\\ \\ x_1+3x_2+0x_3+x_4=30\end{array}$ , 需要将 $x_2$ 的系数变为 $\left(\begin{array}{c}0\\ 1\end{array}\right)$ , $x_3$ 的系数保持 $\left(\begin{array}{c}1\\ 0\end{array}\right)$ 不变 ;

方程 $2$ 同解变换 : 在 $x_1+3x_2+0x_3+x_4=30$ 中 , 需要将 $x_2$ 的系数变成 $1$ , 在方程两端乘以 $\frac{1}{3}$ , 此时方程变成 $\frac{1}{3}{x}_1+x_2+0x_3+\frac{1}{3}{x}_4=10$ ;

方程 $1$ 同解变换 : 将上述方程 $2$ 作同解变换后 , 方程组变成 $\{\begin{array}{l}2x_1+x_2+x_3+0x_4=40\\ \\ \frac{1}{3}{x}_1+x_2+0x_3+\frac{1}{3}{x}_4=10\end{array}$ , 目前的需求是将方程 $1$ 的 $x_2$ 系数变为 $0$ , 使用方程 $1$ 减去 方程 $2$ 即可得到要求的矩阵 :

$$\begin{array}{lcl}(2-\frac{1}{3})x_1+0x_2+x_3+(0-\frac{1}{3})x_4& =& 40-10\\ & & \\ \frac{5}{3}{x}_1+0x_2+x_3-\frac{1}{3}{x}_4& =& 30\end{array}$$

最终方程 $1$ 转化为 $\frac{5}{3}{x}_1+0x_2+x_3-\frac{1}{3}{x}_4=30$ ;

同解变换完成后的方程组为 $\{\begin{array}{l}\frac{5}{3}{x}_1+0x_2+x_3-\frac{1}{3}{x}_4=30\\ \\ \frac{1}{3}{x}_1+x_2+0x_3+\frac{1}{3}{x}_4=10\end{array}$

七、第一次迭代 : 生成新的单纯形表

---

单纯形表变成如下形式 : 下面的单纯形表中 , 上面部分是初始基可行解对应的单纯形表 , 下面的部分是本次迭代后生成的新的单纯形表 ;

将同解变换后的方程组中的 系数矩阵 , 和 常数 , 填入新的单纯形表中 ;

$c_j$	$c_j$		$3$	$4$	$0$	$0$
$C_B$ 基变量系数 (目标函数)	基变量	常数 $b$	$x_1$	$x_2$	$x_3$	$x_4$	${\theta }_i$
$0$ ( 目标函数 $x_3$ 系数 $c_3$ )	$x_3$	$40$	$2$	$1$	$1$	$0$	$40$ ( ${\theta }_3$ )
$0$ ( 目标函数 $x_4$ 系数 $c_4$)	$x_4$	$30$	$1$	$3$	$0$	$1$	$10$ ( ${\theta }_4$ )
${\sigma }_j$			$3$ ( ${\sigma }_1$ )	$4$ ( ${\sigma }_2$ )	$0$	$0$
–	–	–	–	–	–	–	–
$0$ ( 目标函数 $x_3$ 系数 $c_3$ )	$x_3$	$30$	$\frac{5}{3}$	$0$	$1$	$-\frac{1}{3}$	$?$ ( ${\theta }_3$ )
$4$ ( 目标函数 $x_2$ 系数 $c_2$)	$x_2$	$10$	$\frac{1}{3}$	$1$	$0$	$\frac{1}{3}$	$?$ ( ${\theta }_2$ )
${\sigma }_j$ ( 检验数 )			$\frac{5}{3}$ ( ${\sigma }_1$ )	$0$	$0$	$-\frac{4}{3}$ ( ${\sigma }_4$ )

八、第一次迭代 : 解出基可行解

---

新的 基变量是 $x_3,x_2$ , 对应的基矩阵是 $\left(\begin{array}{c}10\\ 01\end{array}\right)$ , 非基变量是 $x_1,x_4$ , 对应的非基矩阵是 $\left(\begin{array}{c}\frac{5}{3}-\frac{1}{3}\\ \frac{1}{3}\frac{1}{3}\end{array}\right)$ , 将非基变量设置为 $0$ , 方程组为 $\{\begin{array}{l}\frac{5}{3}\times 0+0x_2+x_3-\frac{1}{3}\times 0=30\\ \\ \frac{1}{3}\times 0+x_2+0x_3+\frac{1}{3}\times 0=10\end{array}$ , 解出基变量为 $\{\begin{array}{l}{x}_3=30\\ \\ x_2=10\end{array}$ , 基可行解 为 $\left(\begin{array}{c}0\\ 10\\ 30\\ 0\end{array}\right)$

九、第一次迭代 : 计算检验数 ${\sigma }_j$ 判定最优解 并选择入基变量

---

根据 【运筹学】线性规划数学模型 ( 单纯形法 | 最优解判定原则 | 单纯形表 | 系数计算方法 | 根据系数是否小于等于 0 判定最优解 ) 博客中分析 , 检验数计算公式为 :

• 矩阵形式 : $C_N^T-C_B^T{B}^{-1}N$
• 单个检验数计算公式 : ${\sigma }_j=c_j-\sum c_i{a}_{ij}$

基变量的检验数是 $0$ , 主要是求非基变量的检验数 ${\sigma }_1,{\sigma }_4$ ;

${\sigma }_1=c_1-(c_3{a}_{11}+c_2{a}_{12})$

${\sigma }_1=3-(0\times \frac{5}{3})-(4\times \frac{1}{3})=\frac{5}{3}$ , 是从下面的单纯形表中的如下位置提取的数值 ;

${\sigma }_4=c_4-(c_3{a}_{41}+c_2{a}_{42})$

${\sigma }_4=0-(0\times -\frac{1}{3})-(4\times \frac{1}{3})=-\frac{4}{3}$ , 是从下面的单纯形表中的如下位置提取的数值 ;

检验数 $\{\begin{array}{l}{\sigma }_1=3-(0\times \frac{5}{3})-(4\times \frac{1}{3})=\frac{5}{3}\\ \\ {\sigma }_4=0-(0\times -\frac{1}{3})-(4\times \frac{1}{3})=-\frac{4}{3}\end{array}$ , ${\sigma }_1$ 是大于 $0$ 的 , 两个检验数必须都小于等于 $0$ , 该基可行解才算作是最优解 , 因此 该基可行解不是最优解 ;

根据检验数选择入基变量 : 继续迭代 , 选择检验数较大的非基变量 , 作为入基变量 , 这里入基变量是 $x_1$ ;

十、第一次迭代 : 根据入基变量计算并选择出基变量

---

参考博客 【运筹学】线性规划数学模型 ( 单纯形法 | 迭代原则 | 入基 | 出基 | 线性规划求解示例 ) 五、出基与入基变量选择

入基变量 根据检验数 $\sigma$ 选择的是 $x_1$ ;

出基变量是根据 $\theta$ 值来选择的 , 选择 $\theta$ 值较小的值对应的基变量作为出基变量 ;

$\theta$ 值计算 : 常数列 $b=\left(\begin{array}{c}30\\ 10\end{array}\right)$ , 分别除以除以入基变量 $x_1$ 大于 $0$ 的系数列 $\left(\begin{array}{c}\frac{5}{3}\\ \\ \frac{1}{3}\end{array}\right)$ , 计算过程如下 $\left(\begin{array}{c}\frac{30}{\frac{5}{3}}\\ \\ \frac{10}{\frac{1}{3}}\end{array}\right)$ , 得出结果是 $\left(\begin{array}{c}18\\ \\ 30\end{array}\right)$ , 然后选择一个最小值 $18$ , 查看该最小值对应的变量是 $x_3$ , 选择该变量作为出基变量 ;

$x_1$ 作入基变量 , $x_3$ 作出基变量 ; 使用 $x_1$ 替代基变量中 $x_3$ 的位置 ;

迭代后的基变量为 $x_1,x_2$ ;

更新一下单纯形表 :

$c_j$	$c_j$		$3$	$4$	$0$	$0$
$C_B$ 基变量系数 (目标函数)	基变量	常数 $b$	$x_1$	$x_2$	$x_3$	$x_4$	${\theta }_i$
$0$ ( 目标函数 $x_3$ 系数 $c_3$ )	$x_3$	$40$	$2$	$1$	$1$	$0$	$40$ ( ${\theta }_3$ )
$0$ ( 目标函数 $x_4$ 系数 $c_4$)	$x_4$	$30$	$1$	$3$	$0$	$1$	$10$ ( ${\theta }_4$ )
${\sigma }_j$			$3$ ( ${\sigma }_1$ )	$4$ ( ${\sigma }_2$ )	$0$	$0$
–	–	–	–	–	–	–	–
$0$ ( 目标函数 $x_3$ 系数 $c_3$ )	$x_3$	$30$	$\frac{5}{3}$	$0$	$1$	$-\frac{1}{3}$	$18$ ( ${\theta }_3$ )
$4$ ( 目标函数 $x_2$ 系数 $c_2$)	$x_2$	$10$	$\frac{1}{3}$	$1$	$0$	$\frac{1}{3}$	$30$ ( ${\theta }_2$ )
${\sigma }_j$ ( 检验数 )			$\frac{5}{3}$ ( ${\sigma }_1$ )	$0$	$0$	$-\frac{4}{3}$ ( ${\sigma }_4$ )

十一、第二次迭代 : 方程组同解变换

---

当前的方程组为 $\{\begin{array}{l}\frac{5}{3}{x}_1+0x_2+x_3-\frac{1}{3}{x}_4=30\\ \\ \frac{1}{3}{x}_1+x_2+0x_3+\frac{1}{3}{x}_4=10\end{array}$ , 选择 $x_1,x_2$ 作为基变量 , 基矩阵为 $\left(\begin{array}{c}\frac{5}{3}0\\ \\ \frac{1}{3}1\end{array}\right)$ , 进行同解变换 , 将基矩阵转为单位阵 ;

方程 $1$ 同解变换 :

将 $\frac{5}{3}{x}_1+0x_2+x_3-\frac{1}{3}{x}_4=30$ 方程中的 $x_1$ 的系数变为 $1$ , $x_2$ 的系数为 $0$ 保持不变 ;

方程的左右变量乘以 $\frac{3}{5}$ :

$$\begin{array}{lcl}\frac{5}{3}{x}_1+0x_2+x_3-\frac{1}{3}{x}_4& =& 30\\ & & \\ (\frac{5}{3}{x}_1+0x_2+x_3-\frac{1}{3}{x}_4)\times \frac{3}{5}& =& 30\times \frac{3}{5}\\ & & \\ x_1+0x_2+\frac{3}{5}{x}_3-\frac{1}{5}{x}_4& =& 18\end{array}$$

当前方程组变成 $\{\begin{array}{l}{x}_1+0x_2+\frac{3}{5}{x}_3-\frac{1}{5}{x}_4=18\\ \\ \frac{1}{3}{x}_1+x_2+0x_3+\frac{1}{3}{x}_4=10\end{array}$

方程 $2$ 同解变换 : 将方程 $1$ 乘以 $-\frac{1}{3}$ , 与方程 $2$ 相加 ;

① 方程 $1$ 乘以 $-\frac{1}{3}$ :

$$\begin{array}{lcl}(x_1+0x_2+\frac{3}{5}{x}_3-\frac{1}{5}{x}_4)\times -\frac{1}{3}& =& 18\times -\frac{1}{3}\\ & & \\ -\frac{1}{3}{x}_1+0x_2+-\frac{1}{5}{x}_3+\frac{1}{15}{x}_4& =& -6\end{array}$$

② 与方程 $2$ 相加 :

$$\begin{array}{lcl}(-\frac{1}{3}{x}_1+0x_2+-\frac{1}{5}{x}_3+\frac{1}{15}{x}_4)+(\frac{1}{3}{x}_1+x_2+0x_3+\frac{1}{3}{x}_4)& =& -6+10\\ & & \\ 0x_1+x_2-\frac{1}{5}{x}_3+\frac{2}{5}{x}_4& =& 4\end{array}$$

当前方程组变成 $\{\begin{array}{l}{x}_1+0x_2+\frac{3}{5}{x}_3-\frac{1}{5}{x}_4=18\\ \\ 0x_1+x_2-\frac{1}{5}{x}_3+\frac{6}{15}{x}_4=4\end{array}$

十二、第二次迭代 : 生成新的单纯形表

---

更新一下单纯形表 : 将第三次迭代的矩阵填入下面的单纯形表中 ;

$c_j$	$c_j$		$3$	$4$	$0$	$0$
$C_B$ 基变量系数 (目标函数)	基变量	常数 $b$	$x_1$	$x_2$	$x_3$	$x_4$	${\theta }_i$
$0$ ( 目标函数 $x_3$ 系数 $c_3$ )	$x_3$	$40$	$2$	$1$	$1$	$0$	$40$ ( ${\theta }_3$ )
$0$ ( 目标函数 $x_4$ 系数 $c_4$)	$x_4$	$30$	$1$	$3$	$0$	$1$	$10$ ( ${\theta }_4$ )
${\sigma }_j$			$3$ ( ${\sigma }_1$ )	$4$ ( ${\sigma }_2$ )	$0$	$0$
第一次迭代	–	–	–	–	–	–	–
$0$ ( 目标函数 $x_3$ 系数 $c_3$ )	$x_3$	$30$	$\frac{5}{3}$	$0$	$1$	$-\frac{1}{3}$	$18$ ( ${\theta }_3$ )
$4$ ( 目标函数 $x_2$ 系数 $c_2$)	$x_2$	$10$	$\frac{1}{3}$	$1$	$0$	$\frac{1}{3}$	$30$ ( ${\theta }_2$ )
${\sigma }_j$ ( 检验数 )			$\frac{5}{3}$ ( ${\sigma }_1$ )	$0$	$0$	$-\frac{4}{3}$ ( ${\sigma }_4$ )
第二次迭代	–	–	–	–	–	–	–
$3$ ( 目标函数 $x_1$ 系数 $c_1$ )	$x_1$	$18$	$1$	$0$	$\frac{3}{5}$	$-\frac{1}{5}$	$?$ ( ${\theta }_3$ )
$4$ ( 目标函数 $x_2$ 系数 $c_2$)	$x_2$	$4$	$0$	$1$	$-\frac{1}{5}$	$\frac{2}{5}$	$?$ ( ${\theta }_2$ )
${\sigma }_j$ ( 检验数 )			$0$	$0$	$?$ ( ${\sigma }_3$ )	$?$ ( ${\sigma }_4$ )

十三、第二次迭代 : 计算检验数、最优解判定

---

计算检验数 $\sigma$ :

${\sigma }_3=0-3\times \frac{3}{5}-4\times (-\frac{1}{5})=-\frac{9}{5}+\frac{4}{5}=-1$

${\sigma }_4=0-3\times (-\frac{1}{5})-4\times \frac{2}{5}=\frac{3}{5}-\frac{8}{5}=-1$

两个非基变量的检验数都是小于等于 $0$ 的 , 因此该基可行解 $\left(\begin{array}{c}18\\ 4\\ 0\\ 0\end{array}\right)$是最优解 ;