一个例子搞懂单纯形法大M法和两阶段法

news/2025/5/31 16:36:52

文章目录

  • 1. 题目
  • 2. 添加松弛变量
  • 3. 大M法
  • 4. 两阶段法


1. 题目

目标函数:
min⁡z=4x1+x2\min z = 4x_1 + x_2 minz=4x1+x2

约束条件:
s.t.{3x1+x2=34x1+3x2≥6x1+2x2≤4x1,x2≥0\text{s.t.} \begin{cases} 3x_1 + x_2 = 3 \\ 4x_1 + 3x_2 \geq 6 \\ x_1 + 2x_2 \leq 4 \\ x_1, x_2 \geq 0 \end{cases} s.t.3x1+x2=34x1+3x26x1+2x24x1,x20

2. 添加松弛变量

min⁡z=4x1+x2+0x3+0x4\min z = 4x_1 + x_2 + 0x_3 + 0x_4 minz=4x1+x2+0x3+0x4

s.t.{3x1+x2=34x1+3x2−x3=6x1+2x2+x4=4xi≥0,i=1,2,…,4\text{s.t.} \begin{cases} 3x_1 + x_2 = 3 \\ 4x_1 + 3x_2 - x_3 = 6 \\ x_1 + 2x_2 + x_4 = 4 \\ x_i \geq 0, i = 1,2,\dots, 4 \end{cases} s.t.3x1+x2=34x1+3x2x3=6x1+2x2+x4=4xi0,i=1,2,,4

3. 大M法

max⁡−z=−4x1−x2+0x3+0x4−Mx5−Mx6\max -z = -4x_1 - x_2 + 0x_3 + 0x_4 - Mx_5 - Mx_6 maxz=4x1x2+0x3+0x4Mx5Mx6

s.t.{3x1+x2+x5=34x1+3x2−x3+x6=6x1+2x2+x4=4xi≥0,i=1,2,…,6\text{s.t.} \begin{cases} 3x_1 + x_2 + x_5 = 3 \\ 4x_1 + 3x_2 - x_3 + x_6 = 6 \\ x_1 + 2x_2 + x_4 = 4 \\ x_i \geq 0, i = 1,2,\dots, 6 \end{cases} s.t.3x1+x2+x5=34x1+3x2x3+x6=6x1+2x2+x4=4xi0,i=1,2,,6

单纯形表

.CCC.-4-100-M-M
CBC_BCBbbbx1↓x_1 \downarrowx1x2x_2x2x3x_3x3x4x_4x4x5x_5x5x6x_6x6
-M←x5\leftarrow x_5x53[3]10010
-Mx6x_6x6643-1001
0x4x_4x44120100
.σ\sigmaσ.7M-44M-1-M000

入基变量 x1x_1x1,出基变量 x5x_5x5

.CCC.-4-100-M-M
CBC_BCBbbbx1x_1x1x2↓x_2 \downarrowx2x3x_3x3x4x_4x4x5x_5x5x6x_6x6
-4x1x_1x11113\frac{1}{3}310013\frac{1}{3}310
-M←x6\leftarrow x_6x620[53\frac{5}{3}35]-10-43\frac{4}{3}341
0x4x_4x43053\frac{5}{3}3501-13\frac{1}{3}310
.σ\sigmaσ.053\frac{5}{3}35M+13\frac{1}{3}31-M0-73\frac{7}{3}37M+43\frac{4}{3}340

入基变量 x2x_2x2,出基变量 x6x_6x6

.CCC.-4-100-M-M
CBC_BCBbbbx1x_1x1x2x_2x2x3↓x_3 \downarrowx3x4x_4x4x5x_5x5x6x_6x6
-4x1x_1x135\frac{3}{5}531015\frac{1}{5}510115\frac{1}{15}151-15\frac{1}{5}51
-1x2x_2x265\frac{6}{5}5601−35-\frac{3}{5}530−45-\frac{4}{5}5435\frac{3}{5}53
0←x4\leftarrow x_4x4100[1]11-1
.σ\sigmaσ.0015\frac{1}{5}510-M-815\frac{8}{15}158-M-15\frac{1}{5}51

入基变量 x3x_3x3,出基变量 x4x_4x4

.CCC.-4-100-M-M
CBC_BCBbbbx1x_1x1x2x_2x2x3x_3x3x4x_4x4x5x_5x5x6x_6x6
-4x1x_1x125\frac{2}{5}52100-15\frac{1}{5}51-215\frac{2}{15}1520
-1x2x_2x295\frac{9}{5}5901035\frac{3}{5}53−15-\frac{1}{5}510
0x3x_3x3100111-1
.σ\sigmaσ.000−15-\frac{1}{5}51-M-1115\frac{11}{15}1511-M

所有的 σj≤0\sigma_j \leq 0σj0,且所有的人工变量都是非基变量。
得到最优解:
x∗=(25,95,1,0)T,x^* = (\frac{2}{5}, \frac{9}{5}, 1, 0)^T, x=(52,59,1,0)T

z∗=4x1+x2=4∗25+95=175.z* = 4x_1 + x_2 = 4*\frac{2}{5} + \frac{9}{5} = \frac{17}{5}. z=4x1+x2=452+59=517.

4. 两阶段法

  • 第一阶段:
    min⁡w=0x1+0x2+0x3+0x4+x5+x6\min w = 0x_1 + 0x_2 + 0x_3 + 0x_4 + x_5 + x_6 minw=0x1+0x2+0x3+0x4+x5+x6

⇔\Leftrightarrow

max⁡−w=0x1+0x2+0x3+0x4−x5−x6\max -w = 0x_1 + 0x_2 + 0x_3 + 0x_4 - x_5 - x_6 maxw=0x1+0x2+0x3+0x4x5x6

s.t.{3x1+x2+x5=34x1+3x2−x3+x6=6x1+2x2+x4=4xi≥0,i=1,2,…,6\text{s.t.} \begin{cases} 3x_1 + x_2 + x_5 = 3 \\ 4x_1 + 3x_2 - x_3 + x_6 = 6 \\ x_1 + 2x_2 + x_4 = 4 \\ x_i \geq 0, i = 1,2,\dots, 6 \end{cases} s.t.3x1+x2+x5=34x1+3x2x3+x6=6x1+2x2+x4=4xi0,i=1,2,,6

单纯形表

.CCC.0000-1-1
CBC_BCBbbbx1↓x_1 \downarrowx1x2x_2x2x3x_3x3x4x_4x4x5x_5x5x6x_6x6
-1←x5\leftarrow x_5x53[3]10010
-1x6x_6x6643-1001
0x4x_4x44120100
.σ\sigmaσ.74-1000

入基变量 x1x_1x1,出基变量 x5x_5x5

.CCC.0000-1-1
CBC_BCBbbbx1x_1x1x2↓x_2 \downarrowx2x3x_3x3x4x_4x4x5x_5x5x6x_6x6
0x1x_1x11113\frac{1}{3}310013\frac{1}{3}310
-1←x6\leftarrow x_6x620[53\frac{5}{3}35]-10-43\frac{4}{3}341
0x4x_4x43053\frac{5}{3}3501-13\frac{1}{3}310
.σ\sigmaσ.053\frac{5}{3}35-10-73\frac{7}{3}370

入基变量 x2x_2x2,出基变量 x6x_6x6

.CCC.0000-1-1
CBC_BCBbbbx1x_1x1x2x_2x2x3x_3x3x4x_4x4x5x_5x5x6x_6x6
0x1x_1x135\frac{3}{5}531015\frac{1}{5}510115\frac{1}{15}151-15\frac{1}{5}51
0x2x_2x265\frac{6}{5}5601−35-\frac{3}{5}530−45-\frac{4}{5}5435\frac{3}{5}53
0x4x_4x4100111-1
.σ\sigmaσ.000000

第一阶段结束。

  • 第二阶段

max⁡−z=−4x1−x2+0x3+0x4\max -z = -4x_1 - x_2 + 0x_3 + 0x_4 maxz=4x1x2+0x3+0x4

.CCC.-4-100
CBC_BCBbbbx1x_1x1x2x_2x2x3↓x_3 \downarrowx3x4x_4x4
-4x1x_1x135\frac{3}{5}531015\frac{1}{5}510
-1x2x_2x265\frac{6}{5}5601−35-\frac{3}{5}530
0←x4\leftarrow x_4x4100[1]1
.σ\sigmaσ.0015\frac{1}{5}510

入基变量 x3x_3x3,出基变量 x4x_4x4

.CCC.-4-100
CBC_BCBbbbx1x_1x1x2x_2x2x3x_3x3x4x_4x4
-4x1x_1x125\frac{2}{5}52100-15\frac{1}{5}51
-1x2x_2x295\frac{9}{5}5901035\frac{3}{5}53
0x3x_3x310011
.σ\sigmaσ.000−15-\frac{1}{5}51

第二阶段结束。

最优解
x∗=(25,95,1,0)T,x^* = (\frac{2}{5}, \frac{9}{5}, 1, 0)^T, x=(52,59,1,0)T

z∗=4x1+x2=4∗25+95=175.z^* = 4x_1 + x_2 = 4*\frac{2}{5} + \frac{9}{5} = \frac{17}{5}. z=4x1+x2=452+59=517.

图片来自网络

文章来源:https://blog.csdn.net/PursueLuo/article/details/108740660
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:https://dhexx.cn/news/show-3036381.html

相关文章

计算半径从1到10的圆面积,当面积大于100时结束循环,break使用范例

/*计算半径从1到10的圆面积&#xff0c;当面积大于100时结束循环*/#include<stdio.h>#include<stdlib.h>#define PI 3.1415/*定义圆周率为3.1415*/ int main(void){ int r;/*定义半径为整型变量*/ float area;/*定义面积为浮点数*/ printf("%s\t%s\n",&…

论文阅读报告:Taxonomy-aware feature engineering for microbiome classification,Mai Oudah and Andreas Hen

文章目录1. HFE1.1. Feature engineering phase1.2. Correlation-based filtering phase1.3. Information Gain (IGIGIG) based filtering phase1.4. IGIGIG-based leaf filtering phase2. DOI1. HFE Hierarchical Feature Engineering&#xff0c;简写 HFE&#xff0c;包含四…

一个神奇的数字货币,终结了南非小哥每天步行20公里的烦恼

今天BCH基金会的BCH4Bicycles项目迎来了第一位受益人。贫困导致这位南非小哥不得不每日步行20公里上班。今天BCH终结了他的烦恼&#xff0c;他领取了一辆由BCH基金会捐赠的自行车和一台手机&#xff0c;并通过CoinText接收到了自己第一笔BCH。在不久的将来&#xff0c;世界各地…

pip install python-weka-wrapper3 error

文章目录1. Pip 安装 python-weka-wrapper3 错误2. 我的解决过程2.1. JDK2.2. Microsoft VC3. 解决问题1. Pip 安装 python-weka-wrapper3 错误 错误内容 (base) C:\Users\Administrator>pip install python-weka-wrapper3 Collecting python-weka-wrapper3 Downloading py…

linux下tengine2.2.0编译安装、开机启动、反向代理配置及健康检查

tengine2.2.0编译安装、开机启动、反向代理配置及健康检查 tengine是由淘宝发起的一个基于nginx的开源项目&#xff0c;nginx的吞吐量比较高、快速、稳定&#xff0c;而且反向代理和负载均衡使用nginx&#xff0c;也是最常见的。本文介绍在Linux&#xff08;centos&#xff09;…

yii2邮箱发送

yii2 邮件发送 163邮箱 1.在配置文件main-local.php components>[]里面配置 mailer > [ class > yii\swiftmailer\Mailer, // send all mails to a file by default. You have to set // useFileTransport to false and configure a transport // for the mailer …

Java:运算符优先级表

文章目录1. Java 运算符优先级表1. Java 运算符优先级表 运算符结合性[ ] . ( ) &#xff08;方法调用&#xff09;从左向右! ~ -- &#xff08;一元运算符&#xff09;- &#xff08;一元运算符&#xff09; ( )&#xff08;强制类型转换&#xff09;new从右向左\color{red}…

JavaScript-流程语句

1.计算2的n次幂&#xff0c;n可输入&#xff0c;n为自然数。 2. 计算n的阶乘&#xff0c;n可输入&#xff0c;n为自然数。 3. 输入a&#xff0c;b&#xff0c;c&#xff0c;不一样的3个数&#xff0c;打印出最大的。 4. 打印 1-100 的质数。 阅读更多

前端实现app引导页面动画效果

插件描述&#xff1a;jQuery引导插件TourTip 交互式可视化指南网页上的元素。使用方法 步骤1&#xff1a; 将以下标记添加到您的文档的<head> 你还需要复制旁边插件的css文件夹和下载的IMG文件夹。你就会有你需要启动一个引导网页的一切。 / *附加的Tourtip CSS文件* / &…

lcg_magic算法笔记:快速排序

文章目录1. 思想2. 示例3. 代码1. 思想 快速排序是一种常用的排序方法。 快速排序的思想是&#xff1a; 首先在数组中选定一个参考值。 这个参考值的作用是&#xff1a;将整个数组分成两个部分。小于这个参考值的所有值都在参考值的左边&#xff0c;大于这个参考值的所有值都…