改进单纯形法的MATLAB实现

Simplex Algorithm

由George Dantzig发明的单纯形法（simplex
algorithm）在数学优化领域中常用于线性规划问题的数值求解。原单纯形法不是很经济的算法。1953年美国数学家G.B.丹齐克为了改进单纯形法每次迭代中积累起来的进位误差，提出改进单纯形法。其基本步骤和单纯形法大致相同，主要区别是在逐次迭代中不再以高斯消去法为基础，而是由旧基阵的逆去直接计算新基阵的逆，再由此确定检验数。这样做可以减少迭代中的累积误差，提高计算精度，同时也减少了在计算机上的存储量。
[1] “改进单纯形法_百度百科,” Baidu.com, 2018. [Online]. Available: https://baike.baidu.com/item/%E6%94%B9%E8%BF%9B%E5%8D%95%E7%BA%AF%E5%BD%A2%E6%B3%95/558282?fr=aladdin. [Accessed: 28-Mar-2020].

由高中学习的线性规划知识，我们知道线性规划问题的最优解在约束条件所形成的的凸多边形的顶点上取到。
而单纯形法的几何本质可以理解为沿凸多边形的边缘逐个尝试不同的顶点，直道找到最优解。

改进单纯形法？

改进在哪里？

原单纯形法都在上一轮迭代的基础上进行再迭代，在计算机内部由于浮点数的二进制近似表示，会在多次迭代后产生一定的累积误差，而改进单纯形法巧妙的利用了分块矩阵的性质避开了这个大坑。Let’s see how

下面我们列出单纯形法的标准形式和分块后的形式