当线性规划与算法相遇:揭秘单纯形法(Simplex)的独特魅力

传统的解决线性规划问题的方法是图形法、代数法求解,但是图形法解题有极大的局限性,因为一旦变量超过3个,基本上就无法通过图形解决,而代数法虽然可以解题,但对于复杂的问题可能效果较差甚至无法求解!
相比图形法和代数法,单纯形法解决线性规划问题具有以下优势:

  1. 理论基础强:单纯形法是基于线性规划的基本理论,通过系统的迭代过程逐步逼近最优解。它是一种可行的、确定性的算法,能够找到问题的最优解或者确定问题是无界或无解的。
  2. 高效性:在实践中,单纯形法通常能够在合理的时间内找到线性规划问题的最优解。尤其对于具有稀疏性质的问题,单纯形法的性能更为出色。此外,单纯形法的计算复杂度与问题规模的增长呈多项式关系。
  3. 灵活性:单纯形法适用于各种类型的线性规划问题,包括有约束的和无约束的问题。它可以处理多目标函数、等式约束、不等式约束、非线性约束等多种情况。
  4. 可以进行优化:单纯形法可以通过一些优化策略来提高算法的效率,例如早期停止条件、对偶单纯形法等。这些优化措施可以在实际应用中进一步加速算法的执行。
  5. 可解释性强:单纯形法的迭代过程很容易理解和解释。每个迭代步骤都代表着一种改进,可以直观地解释为什么选择某个变量作为进基变量或出基变量,从而得到更优解。

尽管单纯形法具有以上优势,但对于大规模问题或非常稀疏的问题,单纯形法可能会遇到性能瓶颈。在这种情况下,可以考虑使用其他更高效的线性规划算法,例如内点法、启发式算法或者列生成法等。
具体使用单纯形法解题的步骤如下所示:

  1. 化一般型为标准型,求初始基本可信解,建立初始单纯形表;
  2. 求检验数并判断,若已得到最优解,结束计算;否则转入下一步;
  3. 进行基变换,构建新的单纯形表进行迭代;
  4. 重复步骤二、三,直到得出最优解、重复解或无最优解等。

1.化标准型

我们常见的线性规划模型的一般型为:
目标函数: m a x ( m i n ) Z = ∑ j = 1 n c j x j 约束条件: { ∑ j = 1 n a i j x j ≥ ( ⩽ ) b i , i = 1 , 2 , ⋯ , m x j ≥ 0 , j = 1 , 2 , ⋯ , n \begin{aligned} 目标函数:& max(min) Z=\sum_{j=1}^{n}c_{j}x_{j} \\ 约束条件:& \left\{\begin{matrix} \sum_{j=1}^{n}a_{ij}x_{j} \geq (\leqslant )b_{i} ,i=1,2,\cdots ,m\\ \\ x_{j}\geq 0,j=1,2,\cdots ,n \end{matrix}\right. \end{aligned} 目标函数:约束条件:max(min)Z=j=1ncjxj j=1naijxj()bi,i=1,2,,mxj0,j=1,2,,n
单纯形法是一种求解线性规划问题的常用方法,其第一步是将线性规划问题转化为标准型,主要是为了方便后续的计算和迭代。
将线性规划问题转化为标准型的目的有以下几个方面:

  1. 约束条件的统一表示:标准型可以将线性规划问题的约束条件统一表示为等式形式,即将不等式约束和非负约束都表示为等式约束。这样可以简化计算过程和算法的设计。
  2. 约束条件的非负性:标准型要求所有变量的取值都非负,这样可以确保问题的可行解存在。通过引入松弛变量或人工变量,将不等式约束转化为等式约束,并引入非负约束,确保问题的可行性。
  3. 目标函数的最大化或最小化:标准型要求将目标函数转化为最小化形式。对于最大化问题,可以通过将目标函数乘以-1来转化为最小化问题,并利用单纯形法求解。
    通过化为标准型,可以将线性规划问题转化为一个更加结构化和规范化的形式,方便应用单纯形法进行迭代计算。标准型的形式更加适合使用单纯形表格来表示和计算,使得单纯形法的步骤更加清晰和易于理解。
    需要注意的是,并非所有的线性规划问题都能够直接转化为标准型,有些问题需要经过一些额外的转化步骤才能达到标准型的形式。但是,一旦将问题转化为标准型,就可以直接应用单纯形法进行求解。
    具体线性规划模型的标准型为:
    目标函数: m a x Z = ∑ j = 1 n c j x j 约束条件: { ∑ j = 1 n a i j x j = b i , i = 1 , 2 , ⋯ , m x j ≥ 0 , j = 1 , 2 , ⋯ , n \begin{aligned} 目标函数:& maxZ=\sum_{j=1}^{n}c_{j}x_{j} \\ 约束条件:& \left\{\begin{matrix} \sum_{j=1}^{n}a_{ij}x_{j} =b_{i} ,i=1,2,\cdots ,m\\ \\ x_{j}\geq 0,j=1,2,\cdots ,n \end{matrix}\right. \end{aligned} 目标函数:约束条件:maxZ=j=1ncjxj j=1naijxj=bi,i=1,2,,mxj0,j=1,2,,n
    标准型的要求主要为:
    (1)目标函数为求最大值
    (2)约束条件均为等式方程
    (3)变量 x j x_{j} xj为非负
    (4)常数 b i b_{i} bi都大于等于零
    下面引入一个简单的例子并将其化为标准型:
    其线性规划问题的一般型为:
    m a x z = 6 x 1 − 2 x 2 + x 3 { 2 x 1 − x 2 + 2 x 3 ⩽ 2 x 1 + 4 x 3 ⩽ 4 x 1 , x 2 , x 3 ⩾ 0 \begin{aligned} max \ z=6x_{1}-2x_{2}+x_{3} \\ \left\{\begin{matrix} 2x_{1}-x_{2}+2x_{3}\leqslant 2\\ x_{1}+4x_{3}\leqslant 4\\ x_{1},x_{2},x_{3}\geqslant 0 \end{matrix}\right. \end{aligned} max z=6x12x2+x3 2x1x2+2x32x1+4x34x1,x2,x30
    将其转换为标准型为:
    m a x z = 6 x 1 − 2 x 2 + x 3 + 0 x 4 + 0 x 5 { 2 x 1 − 1 x 2 + 2 x 3 + 1 x 4 + 0 x 5 = 2 1 x 1 + 0 x 2 + 4 x 3 + 0 x 4 + 1 x 5 = 4 x 1 , x 2 , x 3 , x 4 , x 5 ≥ 0 \begin{aligned} max \ z=6x_{1}-2x_{2}+x_{3}+0x_{4}+0x_{5}\\ \left\{\begin{matrix} 2x_{1}-1x_{2}+2x_{3}+1x_{4}+0x_{5}=2\\ 1x_{1}+0x_{2}+4x_{3}+0x_{4}+1x_{5}=4\\ x_{1},x_{2},x_{3},x_{4},x_{5}\geq 0 \end{matrix}\right. \end{aligned} max z=6x12x2+x3+0x4+0x5 2x11x2+2x3+1x4+0x5=21x1+0x2+4x3+0x4+1x5=4x1,x2,x3,x4,x50

2.求检验数并判断最优解

在单纯形法中,一旦将线性规划问题转化为标准型,可以通过以下步骤求解检验数,并判断是否达到最优解:

  1. 制作初始单纯形表格:根据标准型的形式,构造初始单纯形表格,包括目标函数的系数、约束条件的系数矩阵、右侧常数项等。
  2. 计算检验数:在初始单纯形表格中,计算每个变量的检验数。检验数表示在目标函数中增加或减少一个单位变量的值时,目标函数值的变化。检验数的计算公式为:检验数 = 目标函数系数 - 系数矩阵中对应列的系数与目标函数系数的乘积。
  3. 判断是否达到最优解:若所有的检验数都为非负数,则当前解为最优解。因为如果存在负的检验数,将会导致目标函数值继续改善,因此需要进行下一步的迭代。
    具体初始单纯形表如下所示:
    在这里插入图片描述
    注意,这里对参数进行解释: c k c_{k} ck表示各个变量的价值系数,这里可以从目标函数 m a x z = 6 x 1 − 2 x 2 + x 3 + 0 x 4 + 0 x 5 max \ z=6x_{1}-2x_{2}+x_{3}+0x_{4}+0x_{5} max z=6x12x2+x3+0x4+0x5可以看出各个变量( x 1 、 x 2 、 x 3 、 x 4 、 x 5 x_{1}、x_{2}、x_{3}、x_{4}、x_{5} x1x2x3x4x5)下面的数值,我们只需要看约束条件即可。
    { 2 x 1 − 1 x 2 + 2 x 3 + 1 x 4 + 0 x 5 = 2 1 x 1 + 0 x 2 + 4 x 3 + 0 x 4 + 1 x 5 = 4 \begin{aligned} \left\{\begin{matrix} 2x_{1}-1x_{2}+2x_{3}+1x_{4}+0x_{5}=2\\ 1x_{1}+0x_{2}+4x_{3}+0x_{4}+1x_{5}=4\\ \end{matrix}\right. \end{aligned} {2x11x2+2x3+1x4+0x5=21x1+0x2+4x3+0x4+1x5=4
    X B X_{B} XB代表基变量,具体基变量的找法,我们只需要在系数矩阵中找到对应的单位矩阵,单位矩阵所对应的变量即为基变量,因此,从初始单纯形表可以看出, x 4 、 x 5 x_{4}、x_{5} x4x5即为基变量。
    c B c_{B} cB代表右侧基变量所对应的价值系数,初始单纯形表的基变量对应的价值系数分别是0、0
    b b b这一列称之为资源限量,填写的时候,只需要看系数矩阵中右侧的数字即可
    求基本可行解时,只需要设置全部非基变量为0,即令 x 1 、 x 2 、 x 3 = 0 x{1}、x{2}、x_{3}=0 x1x2x3=0,这样可以求出基变量 x 4 = 2 , x 5 = 4 x_{4}=2,x_{5}=4 x4=2,x5=4,所以基本可行解为 ( 0 , 0 , 0 , 2 , 4 ) T (0,0,0,2,4)^T (0,0,0,2,4)T
    接下来,我们需要通过计算检验数来判断该解是否为最优解,即分别计算 ( x 1 、 x 2 、 x 3 、 x 4 、 x 5 ) (x_{1}、x_{2}、x_{3}、x_{4}、x_{5}) (x1x2x3x4x5)所对应的检验数 σ j \sigma_{j} σj,其计算方式为 c j − c B x j c_{j}-c_{B}x_{j} cjcBxj
    若当前计算的所有 σ j \sigma_{j} σj都小于等于0,即表示当前基础可行解为最优解,否则还需要进行基变换来进一步求得最优解!

3.基变换

基变换的作用就会帮我们找到下一个可行解,简单来说就是用当前一个非基变量来替换基变量,也就是让非基变量入基,让基变量出基。
在确定哪个非基变量入基的时候,我们只需要看检验数,当前最大的检验数对应的变量,就是需要入基的非基变量,当前是 x 1 x_{1} x1,它对应的检验数为6,是当前所有检验数的最大值。
接下来就需要确定出基变量,那么首先就需要计算 θ \theta θ,其计算方法就是b这一列和确定入基 x 1 x_{1} x1这一列相除得到,即为 ( 1 , 4 ) T (1,4)^T (1,4)T,计算好 θ \theta θ后,我们只需要找 θ \theta θ值中的最小值,其最小值对应的变量 x 4 x_{4} x4就是所对应的出基变量。
紧接着就需要把入基变量 x 1 x_{1} x1和出基变量 x 4 x_{4} x4中相交的数字经过行列运算变换为1,其对应的同列元素全部变换为0,计算完毕后就对应着下一个单纯形表。
具体下一个单纯形表如下所示:
在这里插入图片描述
这里需要注意,你可能看到 θ \theta θ对对应的第一个数填的是-,那是因为当前入基 x 2 x_{2} x2对应的数字为负数,因此不需要计算。
由此可以看出,并不是所有的检验数都为小于等于零,因此,当前的基础可行解仍然不是最优解,所以还需要进一步进行基变换,下面就不对基变换的过程进行详解,仅仅展示最优求得最优解所对应的单纯形表。
在这里插入图片描述
至此,从图表中可以看出,全部变量对应的检验数都小于等于0,因此,此时多对应的解为最优解,最优解为 ( 4 , 6 , 0 , 0 , 0 ) T (4,6,0,0,0)^T (4,6,0,0,0)T

4.结果判定方法

具体解的判定方法如下所示:
在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/171981.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

python DevOps

在云原生中,python扮演的角色是什么? 在云原生环境中,Python 作为一种高级编程语言,在多个方面扮演着重要角色。云原生是指利用云计算的各种优势(如弹性、可扩展性和自动化),构建和运行应用程序…

Transformer英语-法语机器翻译实例

依照Transformer结构来实例化编码器-解码器模型。在这里,指定Transformer编码器和解码器都是2层,都使用4头注意力。为了进行序列到序列的学习,我们在英语-法语机器翻译数据集上训练Transformer模型,如图11.2所示。 da…

【Linux】部署单体项目以及前后端分离项目(项目部署)

一、简介 以下就是Linux部署单机项目和前后端分离项目的优缺点,希望对你有所帮助。 1、Linux部署单机项目: 优点: 简化了系统管理:由于所有服务都在同一台机器上运行,因此可以简化系统管理和维护。提高了性能&#x…

Spring Boot集成RESTful API

在Spring Boot中集成一个RESTful API是我们在实际开发中较为常见的一种开发任务,以下通过一个小的案例来展示在Spring Boot中创建RESTful API来编写一个单元测试。 本节使用到的注解: Controller:修饰class,用来创建处理http请求的…

(免费领源码) Asp.Net#SQL Server校园在线投票系统10557-计算机毕业设计项目选题推荐

摘 要 随着互联网大趋势的到来,社会的方方面面,各行各业都在考虑利用互联网作为媒介将自己的信息更及时有效地推广出去,而其中最好的方式就是建立网络管理系统,并对其进行信息管理。由于现在网络的发达,校园投票通过网…

吃瓜教程3|决策树

ID3算法 假定当前样本集合D中第k类样本所占比例为pk,则样本集合D的信息熵定义为 信息增益 C4.5算法 ID3算法存在一个问题,就是偏向于取值数目较多的属性,因此C4.5算法使用了“增益率”(gain ratio)来选择划分属性 CA…

财务数字化转型是什么?_光点科技

财务数字化转型是当今企业发展中的一项关键策略,旨在借助先进的数字技术,重新塑造和优化财务管理体系,以适应迅速变化的商业环境。这一转型不仅仅是技术的升级,更是对企业财务理念和流程的全面升级和改革。 财务数字化转型的核心在…

2.23每日一题(反常积分收敛性的判断)

解法一:用定义(当被积函数的原函数比较好找时): 积分结果为存在则收敛,不存在则发散。 解法二:通过p积分的比较法判断敛散性: 即被积函数与p积分相比较,使得两者同敛散;再…

Centos安装RabbitMQ,JavaSpring发送RabbitMQ延迟延时消息,JavaSpring消费RabbitMQ消息

1,版本说明 erlang 和 rabbitmq 版本说明 https://www.rabbitmq.com/which-erlang.html 确认需要安装的mq版本以及对应的erlang版本。 2,下载安装文件 RabbitMQ下载地址: https://packagecloud.io/rabbitmq/rabbitmq-server Erlang下载地…

如何公网远程访问OpenWRT软路由web界面

文章目录 1.openWRT安装cpolar2.配置远程访问地址3.固定公网地址 简单几步实现在公网环境下远程访问openWRT web 管理界面,使用cpolar内网穿透创建安全隧道映射openWRT web 界面面板443端口,无需公网IP,无需设置路由器。 1.openWRT安装cpola…

嵌入式PID算法理论+实践分析

1.1 概述 比例(Proportion)积分(Integral)微分(Differential)控制器(PID控制器或三项控制器)是一种采用反馈的控制回路机制,广泛应用于工业控制系统和需要连续调制控制的…

10款轻量型的嵌入式GUI库分享

LVGL LittlevGL是一个免费的开源图形库,提供了创建嵌入式GUI所需的一切,具有易于使用的图形元素、漂亮的视觉效果和低内存占用。 特点: 强大的构建模组 按钮、图表、列表、滑块、图像等 ​先进的图形 动画、反锯齿、半透明、平滑滚动 多样…

山西电力市场日前价格预测【2023-10-29】

日前价格预测 预测说明: 如上图所示,预测明日(2023-10-29)山西电力市场全天平均日前电价为318.01元/MWh。其中,最高日前电价为537.50元/MWh,预计出现在18:15。最低日前电价为0.00元/MWh,预计出…

C#WinformListView实现缺陷图片浏览器

C#&Winform&ListView实现缺陷图片浏览器 功能需求图像浏览行间距调整悬浮提示 功能需求 机器视觉检测系统中特别是缺陷检测系统,通常需要进行对已经检出的缺陷图片进行浏览查阅。主要是通过条件筛选查询出所需要的数据,进行分页再展示到界面中。…

反向传播神经网络(BPNN)的实现(Python,附源码及数据集)

文章目录 一、理论基础1、前向传播2、反向传播3、激活函数4、神经网络结构 二、BP神经网络的实现1、训练过程(BPNN.py)2、测试过程(test.py)3、测试结果4、参考源码及实验数据集 一、理论基础 反向传播神经网络(BPNN&…

微服务-统一网关Gateway

网关的作用 对用户请求做身份认证、权限校验将用户请求路由到微服务,并实现负载均衡对用户请求做限流 搭建网关服务 创建新module,命名为Gateway,引入依赖(1.SpringCloudGateway依赖;2.Eureka客户端依赖或者nacos的服…

【开源】基于SpringBoot的城市桥梁道路管理系统的设计和实现

目录 一、摘要1.1 项目介绍1.2 项目录屏 二、功能模块三、系统展示四、核心代码4.1 查询城市桥梁4.2 新增城市桥梁4.3 编辑城市桥梁4.4 删除城市桥梁4.5 查询单个城市桥梁 五、免责说明 一、摘要 1.1 项目介绍 基于VueSpringBootMySQL的城市桥梁道路管理系统,支持…

40.查找练习题(王道2023数据结构第7章)

试题1&#xff08;王道7.2.4节综合练习5&#xff09;&#xff1a; 写出折半查找的递归算法。 #include<stdio.h> #include<stdlib.h> #include<string.h>#define MAXSIZE 10 #define ElemType int #define Status inttypedef struct{int data[MAXSIZE]; /…

Python---break关键字对for...else结构的影响

for循环中添加else结构 循环可以和else配合使用&#xff0c; else下方缩进的代码指的是当循环正常结束之后要执行的代码。 强调&#xff1a; 循环 正常结束&#xff0c;else之后要执行的代码。 非正常结束&#xff0c;其else中的代码是不会执行的。&#xff08;如遇到br…

[计算机提升] Windows系统各种开机启动方式介绍

1.14 开机启动 在Windows系统中&#xff0c;开机启动是指开启电脑后&#xff0c;自动运行指定的程序或服务的技术。一些程序或服务需要在开机后自动启动&#xff0c;以便及时响应用户操作&#xff0c;比如防安防软件、即时通信工具、文件同步软件等。 同时&#xff0c;一些系统…