单纯形法Python实现

最新推荐文章于 2025-07-29 19:45:28 发布

运筹小菜鸡

最新推荐文章于 2025-07-29 19:45:28 发布

阅读量6.9k

点赞数 10

CC 4.0 BY-SA版权

文章标签： python 算法

本文链接：https://blog.csdn.net/qq_43434469/article/details/123999470

本文基于《运筹学教程》第四版介绍了单纯形法的理论基础，并提供了Python代码实现，包括循环求解、最优性检验和结果输出。代码仅考虑模型有最优解的情况，未涵盖无最优解或无限最优解的处理。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档

前言

本文代码参考《运筹学教程》（第四版）中的第一章第四节中所介绍的单纯形法理论基础所实现，其中的各种符号表示与书中的符号对应。

其次，本文中只考虑模型有最优解的情况，不考虑模型无最优解或无限最优解。

一、单纯形法理论基础

单纯形法的理论可以在网上搜索到许多相关资料，在此不多赘述，直接进入正题。

二、Python代码实现

单纯形法的实现可以分为以为几个模块：循环、最优性检验、结果输出

1.循环求解

这里是单纯形法的主体部分，主要是分为以下几个步骤：求解检验数σ、确定换入基变量、确定换出基向量、对矩阵A进行高斯变换。以此循环。

代码如下：

def pivot(j,a,b):
    Cb=[c[i] for i in j]      #更新Cb

    # 对sigma进行求解
    sigma = np.array([0]*a.shape[1]).astype(float)   #求解sigma，用做后续的判定
    for i in range(a.shape[1]):

最低0.47元/天解锁文章

200万优质内容无限畅学