单纯形法求解步骤：一个简单例子

1. 题目

使用单纯形法求解
$$\max z=50x_1+100x_2$$

$$\text{s.t.}\{\begin{array}{l}{x}_1+x_2\le 300\\ 2x_1+x_2\le 400\\ x_2\le 250\\ x_i\ge 0,i=1,2\end{array}$$

2. 求解过程

• 首先引入松弛变量，将线性规划问题化成标准型。
  $$\max z=50x_1+100x_2+0x_3+0x_4+0x_5$$

$$\text{s.t.}\{\begin{array}{l}{x}_1+x_2+x_3=300\\ 2x_1+x_2+x_4=400\\ x_2+x_5=250\\ x_i\ge 0,i=1,2,\dots ,5\end{array}$$

2.1. 第一次迭代

1. 首先画出如下的表格；（PS：表格第一行有 3 个 $C$、$\sigma$，应该合并，但是我不知道如何在Markdown中合并表格）

$C$	$C$	$C$
$C_B$	基	$b$	$x_1$	$x_2$	$x_3$	$x_4$	$x_5$
$\sigma$	$\sigma$	$\sigma$

2. 在第一行中，$C_j$ 表示目标函数的系数矩阵（向量），$C_j$ 右侧是目标函数的系数矩阵（向量）具体值，在题目中目标函数为 $\max z=50x_1+100x_2+0x_3+0x_4+0x_5$，因此按照 $x_1,x_2,\dots ,x_5$ 的顺序写出它们的系数，
   $$C=\left[\begin{array}{ccccc}50& 100& 0& 0& 0\end{array}\right]$$

$C$	$C$	$C$	$50$	$100$	$0$	$0$	$0$
$C_B$	基	$b$	$x_1$	$x_2$	$x_3$	$x_4$	$x_5$
$\sigma$	$\sigma$	$\sigma$

3. 在 $x_1,x_2,x_3,x_4,x_5$ 下面写上约束方程组的系数矩阵 $A$；需要注意的是：在方程中没有出现的变量，其系数为 0，必须写上；
   $$A=(a_{ij})=\left[\begin{array}{ccccc}1& 1& 1& 0& 0\\ 2& 1& 0& 1& 0\\ 0& 1& 0& 0& 1\end{array}\right]$$

$C$	$C$	$C$	50	100	0	0	0
$C_B$	基	$b$	$x_1$	$x_2$	$x_3$	$x_4$	$x_5$
			$1$	$1$	$1$	$0$	$0$
			$2$	$1$	$0$	$1$	$0$
			$0$	$1$	$0$	$0$	$1$
$\sigma$	$\sigma$	$\sigma$

4. 在 $b$ 列依次写上约束矩阵右侧的常数，分别是
   $$b=\left[\begin{array}{c}300\\ 400\\ 250\end{array}\right]$$

$C$	$C$	$C$	50	100	0	0	0
$C_B$	基	$b$	$x_1$	$x_2$	$x_3$	$x_4$	$x_5$
		$300$	1	1	1	0	0
		$400$	2	1	0	1	0
		$250$	0	1	0	0	1
$\sigma$	$\sigma$	$\sigma$

5. 在系数矩阵中找到一个单位矩阵，其对应的变量作为初始基变量，因此基变量为
   $$\text{基}=\left[\begin{array}{c}{x}_3\\ x_4\\ x_5\end{array}\right]$$
   将其填入基列；

$C$	$C$	$C$	50	100	0	0	0
$C_B$	基	$b$	$x_1$	$x_2$	$x_3$	$x_4$	$x_5$
	$x_3$	300	1	1	1	0	0
	$x_4$	400	2	1	0	1	0
	$x_5$	250	0	1	0	0	1
$\sigma$	$\sigma$	$\sigma$

6. 在 $C_B$ 列中填入基向量在目标函数的系数矩阵（向量）中值
   $$C_B=\left[\begin{array}{c}0\\ 0\\ 0\end{array}\right]$$

$C$	$C$	$C$	50	100	0	0	0
$C_B$	基	$b$	$x_1$	$x_2$	$x_3$	$x_4$	$x_5$
$0$	$x_3$	300	1	1	1	0	0
$0$	$x_4$	400	2	1	0	1	0
$0$	$x_5$	250	0	1	0	0	1
$\sigma$	$\sigma$	$\sigma$

7. 至此，迭代1中的初始单纯形表填入完毕，接下来计算检验和

$${\sigma }_j=C_j-\sum _i{C}_B\ast a_{ij}$$

因此，

$${\sigma }_1=C_1-\sum _{i=1}^3{C}_B\ast a_{i1}=50-{\left[\begin{array}{c}0\\ 0\\ 0\end{array}\right]}^T\ast \left[\begin{array}{c}1\\ 2\\ 0\end{array}\right]=50$$

$${\sigma }_2=C_2-\sum _{i=1}^3{C}_B\ast a_{i2}=100-{\left[\begin{array}{c}0\\ 0\\ 0\end{array}\right]}^T\ast \left[\begin{array}{c}1\\ 1\\ 1\end{array}\right]=100$$

$${\sigma }_3=C_3-\sum _{i=1}^3{C}_B\ast a_{i3}=0-{\left[\begin{array}{c}0\\ 0\\ 0\end{array}\right]}^T\ast \left[\begin{array}{c}1\\ 0\\ 0\end{array}\right]=0$$

$${\sigma }_4=C_4-\sum _{i=1}^3{C}_B\ast a_{i4}=0-{\left[\begin{array}{c}0\\ 0\\ 0\end{array}\right]}^T\ast \left[\begin{array}{c}0\\ 1\\ 0\end{array}\right]=0$$

$${\sigma }_5=C_5-\sum _{i=1}^3{C}_B\ast a_{i5}=0-{\left[\begin{array}{c}0\\ 0\\ 0\end{array}\right]}^T\ast \left[\begin{array}{c}0\\ 0\\ 1\end{array}\right]=0$$

因此在$\sigma$ 行填入
$$\sigma =\left[\begin{array}{ccccc}50& 100& 0& 0& 0\end{array}\right]$$

$C$	$C$	$C$	50	100	0	0	0
$C_B$	基	$b$	$x_1$	$x_2$	$x_3$	$x_4$	$x_5$
0	$x_3$	300	1	1	1	0	0
0	$x_4$	400	2	1	0	1	0
0	$x_5$	250	0	1	0	0	1
$\sigma$	$\sigma$	$\sigma$	$50$	$100$	$0$	$0$	$0$

8. 因为 $\exists j$，使得 ${\sigma }_j>0$ 成立，所以还需要继续进行迭代；选择使得 ${\sigma }_j$ 取得最大值的变量，作为入基变量；

$$k=\mathrm{argmax}{\sigma }_j,{\sigma }_j>0$$

$$\sigma =\left[\begin{array}{ccccc}50& 100& 0& 0& 0\end{array}\right]$$
中最大值为100，其对应的变量为 $x_2$（位置），所以 $x_2$ 作为入基变量，$k=2$；

9. 计算 ${\theta }_j$ 确认出基变量，

$${\theta }_j=b_j/a_{jk}$$

求使得 ${\theta }_j$ 取得最小值的变量，作为出基变量；

$${\theta }_1=b_1/a_{12}=300/1=300,(x_3)$$

$${\theta }_2=b_2/a_{22}=400/1=400,(x_4)$$

$${\theta }_3=b_3/a_{32}=250/1=250,(x_5)$$

所以，$\min {\theta }_j=250$，对应的变量为 $x_5$ 作为出基变量；

2.2. 第二次迭代

1. 入基变量 $x_2$，出基变量 $x_5$；

$C$	$C$	$C$	50	100	0	0	0
$C_B$	基	$b$	$x_1$	$x_2$	$x_3$	$x_4$	$x_5$
0	$x_3$	300	1	1	1	0	0
0	$x_4$	400	2	1	0	1	0
$100$	$x_2$	250	0	1	0	0	1
$\sigma$	$\sigma$	$\sigma$

2. 进行初等行变换，将基变量 $x_2$ 所在列中，除了基变量 $x_2$ 所在的行为1，其余全部消去；

$$(-1)\ast \left[\begin{array}{cccccc}250& 0& 1& 0& 0& 1\end{array}\right]+1\ast \left[\begin{array}{cccccc}300& 1& 1& 1& 0& 0\end{array}\right]=\left[\begin{array}{cccccc}50& 1& 0& 1& 0& -1\end{array}\right]$$

$$(-1)\ast \left[\begin{array}{cccccc}250& 0& 1& 0& 0& 1\end{array}\right]+1\ast \left[\begin{array}{cccccc}400& 2& 1& 0& 1& 0\end{array}\right]=\left[\begin{array}{cccccc}150& 2& 0& 0& 1& -1\end{array}\right]$$

$C$	$C$	$C$	50	100	0	0	0
$C_B$	基	$b$	$x_1$	$x_2$	$x_3$	$x_4$	$x_5$
0	$x_3$	$50$	$1$	$0$	$1$	$0$	$-1$
0	$x_4$	$150$	$2$	$0$	$0$	$1$	$-1$
100	$x_2$	250	0	1	0	0	1
$\sigma$	$\sigma$	$\sigma$

3. 计算检验和

$${\sigma }_j=C_j-\sum _i{C}_B\ast a_{ij}$$

$${\sigma }_1=C_1-C_B\ast a_{i1}=50-{\left[\begin{array}{c}0\\ 0\\ 100\end{array}\right]}^T\ast \left[\begin{array}{c}1\\ 2\\ 0\end{array}\right]=50,$$

$${\sigma }_2=C_2-C_B\ast a_{i2}=100-{\left[\begin{array}{c}0\\ 0\\ 100\end{array}\right]}^T\ast \left[\begin{array}{c}0\\ 0\\ 1\end{array}\right]=0,$$

$${\sigma }_3=C_3-C_B\ast a_{i3}=0-{\left[\begin{array}{c}0\\ 0\\ 100\end{array}\right]}^T\ast \left[\begin{array}{c}1\\ 0\\ 0\end{array}\right]=0,$$

$${\sigma }_4=C_4-C_B\ast a_{i4}=0-{\left[\begin{array}{c}0\\ 0\\ 100\end{array}\right]}^T\ast \left[\begin{array}{c}0\\ 1\\ 0\end{array}\right]=0,$$

$${\sigma }_5=C_5-C_B\ast a_{i5}=0-{\left[\begin{array}{c}0\\ 0\\ 100\end{array}\right]}^T\ast \left[\begin{array}{c}-1\\ -1\\ 1\end{array}\right]=-100,$$

$C$	$C$	$C$	50	100	0	0	0
$C_B$	基	$b$	$x_1$	$x_2$	$x_3$	$x_4$	$x_5$
0	$x_3$	50	1	0	1	0	-1
0	$x_4$	150	2	0	0	1	-1
100	$x_2$	250	0	1	0	0	1
$\sigma$	$\sigma$	$\sigma$	$50$	$0$	$0$	$0$	$-100$

4. 因为 $\exists j$，使得 ${\sigma }_j>0$ 成立，所以还需要继续进行迭代；选择使得 ${\sigma }_j$ 取得最大值的变量，作为入基变量，

$$k=\mathrm{argmax}{\sigma }_j,{\sigma }_j>0$$

$$\sigma =\left[\begin{array}{ccccc}50& 0& 0& 0& -100\end{array}\right]$$
中最大值为50，其对应的变量为 $x_1$（位置），所以 $x_1$ 作为入基变量，$k=1$；

5. 计算 ${\theta }_j$ 确认出基变量，

$${\theta }_j=b_j/a_{jk}$$

求使得 ${\theta }_j$ 取得最小值的变量，作为出基变量；

$${\theta }_1=b_1/a_{11}=50/1=50,(x_3)$$

$${\theta }_2=b_2/a_{21}=150/2=75,(x_4)$$

$${\theta }_3=b_3/a_{31}=250/0=?,(x_5)$$

所以，$\min {\theta }_j=50$，对应的变量为 $x_3$ 作为出基变量；

2.3. 第三次迭代

1. 入基变量 $x_1$，出基变量 $x_3$；

$C$	$C$	$C$	50	100	0	0	0
$C_B$	基	$b$	$x_1$	$x_2$	$x_3$	$x_4$	$x_5$
$50$	$x_1$	50	1	0	1	0	-1
0	$x_4$	150	2	0	0	1	-1
$100$	$x_2$	250	0	1	0	0	1
$\sigma$	$\sigma$	$\sigma$

2. 进行初等行变换，将基变量 $x_1$ 所在列中，除了基变量 $x_1$ 所在的行为1，其余全部消去；

$$(-2)\ast \left[\begin{array}{cccccc}50& 1& 0& 1& 0& -1\end{array}\right]+1\ast \left[\begin{array}{cccccc}150& 2& 0& 0& 1& -1\end{array}\right]=\left[\begin{array}{cccccc}50& 0& 0& -2& 1& 1\end{array}\right]$$

$C$	$C$	$C$	50	100	0	0	0
$C_B$	基	$b$	$x_1$	$x_2$	$x_3$	$x_4$	$x_5$
50	$x_1$	50	1	0	1	0	-1
0	$x_4$	$50$	$0$	$0$	$-2$	$1$	$1$
100	$x_2$	250	0	1	0	0	1
$\sigma$	$\sigma$	$\sigma$

3. 计算检验和
   $${\sigma }_j=C_j-\sum _i{C}_B\ast a_{ij}$$

$${\sigma }_1=C_1-C_B\ast a_{i1}=50-{\left[\begin{array}{c}50\\ 0\\ 100\end{array}\right]}^T\ast \left[\begin{array}{c}1\\ 0\\ 0\end{array}\right]=0,$$

$${\sigma }_2=C_2-C_B\ast a_{i2}=100-{\left[\begin{array}{c}50\\ 0\\ 100\end{array}\right]}^T\ast \left[\begin{array}{c}0\\ 0\\ 1\end{array}\right]=0,$$

$${\sigma }_3=C_3-C_B\ast a_{i3}=0-{\left[\begin{array}{c}50\\ 0\\ 100\end{array}\right]}^T\ast \left[\begin{array}{c}1\\ -2\\ 0\end{array}\right]=-50,$$

$${\sigma }_4=C_4-C_B\ast a_{i4}=0-{\left[\begin{array}{c}50\\ 0\\ 100\end{array}\right]}^T\ast \left[\begin{array}{c}0\\ 1\\ 0\end{array}\right]=0,$$

$${\sigma }_5=C_5-C_B\ast a_{i5}=0-{\left[\begin{array}{c}50\\ 0\\ 100\end{array}\right]}^T\ast \left[\begin{array}{c}-1\\ 1\\ 1\end{array}\right]=-50,$$

$C$	$C$	$C$	50	100	0	0	0
$C_B$	基	$b$	$x_1$	$x_2$	$x_3$	$x_4$	$x_5$
50	$x_1$	50	1	0	1	0	-1
0	$x_4$	50	0	0	-2	1	1
100	$x_2$	250	0	1	0	0	1
$\sigma$	$\sigma$	$\sigma$	$0$	$0$	$-50$	$0$	$-50$

4. 因为 $\forall j$，使得 ${\sigma }_j\le 0$ 成立，所以迭代终止，算法结束；

5. 最终解

$C$	$C$	$C$	50	100	0	0	0
$C_B$	基	$b$	$x_1$	$x_2$	$x_3$	$x_4$	$x_5$
50	$x_1$	$50$	1	0	1	0	-1
0	$x_4$	$50$	0	0	-2	1	1
100	$x_2$	$250$	0	1	0	0	1
$\sigma$	$\sigma$	$\sigma$	0	0	-50	0	-50

解为：
$$\{\begin{array}{l}{x}_1=50;\\ x_2=250;\\ x_3=0;\\ x_4=50;\\ x_5=0.\end{array}$$

$$\max z=50x_1+100x_2=50\ast 50+100\ast 250=27500.$$

3. 算法

1、 标准化（引入松弛变量）；
2、构造初始单纯形表；
3、确定基变量和非基变量；
4、计算检验和；
5、确定入基变量和出基变量；
6、初等行变换；
7、重复4-6步，直到所有的检验和小于等于0；
8、$b$ 值就是对应的基变量的解，非基变量的解为0。

4. 手稿

5. 参考

1. cruelty_angel：单纯形法
2. 妙恋非恋：单纯形法的计算步骤
3. mohu：常用数学符号的 LaTeX 表示方法
4. CodeAllen的博客：Markdown文字添加颜色方法总结（珍藏）