教育房产时事环球科技商业
投稿投诉
商业财经
热点动态
科技数码
软件应用
国际环球
晨报科学
新闻时事
信息智能
汽车房产
办公手机
教育体育
生活生物

运筹学单纯形法(单纯形法各个步骤详解)

  运筹学单纯形法(单纯形法各个步骤详解)
  R帷幄』原创
  作者:臧永森
  作者:臧永森,清华大学工业工程系在读博士,研究方向:运筹优化算法的设计与应用、数据统计分析、大数据技术与应用,戚铭尧老师团队
  编者按
  此文属于电子书线性规划专题第三章单纯形法的内容。在前面的文章中,我们为引入单纯形法介绍了可行域、最优解、可行解、基解、基可行解等基础概念,也阐述了它们之间的关系(具体可见文章《在单纯形法之前》)。在明确了这些基本概念之后,这一节我们来探讨单纯形法的思想逻辑和求解步骤。
  我们已经知道,优化问题的最优解一定是基可行解,那么如何找到最优的基可行解就是最优化问题的求解思路。因此,单纯形法在求解过程,就是不断地寻求变量出入基的循环迭代过程,每次迭代都达到降低目标函数值(或增大目标函数值)的目的,最终得到最优解。那么在迭代过程中,如何使解在改善过程中向着最优解的方向尽快地收敛呢?我们下面用比较直观的方式来解析这个过程。
  单纯形法的基本思想与逻辑
  本文采用的思路参考Dimitris Bertsimas和 John N. Tsitsiklis在 Introduction to Linear Optimization一书中提出的方法[1]。考虑如下标准线性规划问题:
  我们将矩阵A拆分为n个列元素:A1, A2, A3,, An,那么我们可以将问题看成是满足非负约束(4)、凸约束(3)以及约束(5)的最小化问题。
  结合式(3)和(5)我们可以看出,原优化问题转化为求解能够构造出(b, z)的使得z值最小的关于(Ai, ci)的凸组合。为了更好地理解它们之间的几何关系,我们将一个平面视作包含A的一个m维空间,将与ci相关的成本项看作是一维垂直数轴,这时每一个点(Ai, ci)都可以唯一在该三维坐标系中表示出来,如图1所示:
  图1 线性规划问题1—4的"列几何"图示
  我们将(b, z)同样视为一条垂直线表示在图1中,这条垂直线叫做需求线,其与平面的交点是(b, 0)点。需求线与(Ai, ci)的凸组合在几何上有一定的关系,它们或相交或相离,这取决于我们对(Ai, ci)凸组合的选取,选取的凸组合不一样,几何关系就不同。很容易能理解,如果需求线和凸组合相交,说明(b, z)可以用相应的凸组合表示出来,也就表明这个凸组合就是原问题的一个可行解;而如果相离,则说明这个凸组合不满足能够表达(b, z)的条件,也就不是原问题的可行解。所有的凸组合构成了一个凸包,如果需求线能够与凸包相交,那么原问题就存在可行解,如果需求线不能与凸包相交,说明原问题无解。进一步将图1抽象,得到图2,从图中我们可以看出,点I、H、G就是三个不同的凸组合与需求线的交点,也就是原问题的三个可行解。
  图2 可行解的"列几何"图示
  经过上面的分析我们得知,要找到最优解,就是找到与需求线相交的使得z值最小的凸组合。那么如何找这样的凸组合呢?首先引入两个定义:
  如果向量
  是线性独立的,那么向量
  被称为Rn空间中的仿射独立或者仿射无关,其中k<=n。
  在Rn空间中由k+1个仿射无关向量组成的凸包被称为k维单纯形。
  对模型(1—4)来说,总共有m+1个等式约束,假定约束系数矩阵是满秩的,那么一个基可行解将对应m+1个线性独立的列向量,也就意味着有m+1个基点,根据上述定义,由基点之间的差向量线性独立可以得到其仿射独立,由此可以知道它们组成的凸包是m维单纯形。
  假设m维单纯形与需求线相交于点(b, z),由(5)知用来表示(b, z)的线性组合的权重向量是xi ,该向量就是一个基可行解,也就对应我们上节所分析的基变量的内容,当然z就是相应的目标函数值。我们用图2做一个解释,阴影区域的三角形CDF,就是一个2维单纯形,其与需求线的交点H点就是基可行解,点C、D、F是基点。
  我们对二维单纯形CDF做一些改变,会发现相应的z值(与需求线的交点)也会变化,比如我们令基点B取代基点F,单纯形变为BCD,这时可行解变为I点,相应的z值较之前有所增长。类似地,若点E取代点C成为基点,单纯形由CDF变为EDF,可行解就出现在G点,此时z值有所减小。从这些变化中我们找出这样一个规律,当且仅当新加入基的点在当前单纯形平面上方(下方)时,所得的交点(即可行解)对应的z值会增大(减小)。
  如果我们更加形象地描述这个基点变化的过程,就如同用手抓住单纯形CDF的基点C,保持D点和F点固定不变,用力向上拉(向下拉),将C点拉到B点(E点),也就产生了新的单纯形BCD(EDF)。单纯形法的旋转迭代过程,就是不断找到基点向上拉(向下拉)到新基点形成新单纯形的过程。
  单纯形法的求解过程
  简单总结一下单纯形法的求解原理。先找到一个基可行解,然后从非基解中找一个比较有前途的点入基,替换掉基可行解中有待改善的基点,从而达到改善目标函数的目的,如此重复迭代,直至无法找到可以入基的点。
  下面我们用一个例题来演示单纯形法的求解过程。用单纯形法求解如下LP问题:
  第一步:将上述LP转化为标准形式,目的是能够在初始单纯形表中很容易地获得初始基可行解。
  第二步,将标准LP列入第0个单纯形表,如表1:
  表1 单纯形表0
  上述单纯形表中可以看出初始基变量是(s1,s2,s3),从表中找一个能够入基的变量,要求该变量入基后能够使得目标函数值增大量最大。决策变量在第0行的系数看成是这个变量的缩减成本,就是当这个变量增加1时,目标函数z的值将减少的量。比如x1的系数是-2,就说明当x1每增加1,z值将减少-2,也就是增加2。因此如果我们要选择能够使目标函数增加量最大的量入基,应该选择第0行中系数最小的负值(读者可以考虑下为什么必须是负值)。因此这里选择x2入基。
  那如何选择出基变量呢?这里我们采用比值法,用右端项的值(即rhs列)除以出基变量对应的列系数(红色线框标注),从中选择最小的比值对应的基变量出基。如果不选择最小比值对应的基变量出基,将会导致后面的迭代过程出现负的右端项,相应行的基变量将为负值,这与LP标准型的变量非负约束相违背,因此这种操作是不被允许的。所以,表1中的比值优胜者是3,因此s3出基(蓝色线框标注)。
  第三步:通常我们会在x2所在列与s3所在行交汇点圈一个圈,也就是元素4。这表示这一点是我们的转轴点,通过初等变化,将该元素所在的行与列的其他元素变为0,该元素本身变为1,得到下一个单纯形表,如表2所示:
  表2 单纯形表1
  第四步:继续在第0行找负系数对应的入基变量,发现x1对应的系数是-2,可以入基。同时比值运算发现s1对应的变量需要出基,因此第一行、x1列对应的元素1是转轴点,圈一个圈,并进行列运算,得表3:
  表3 单纯形表2
  第五步:继续上述计算,注意这里因为入基变量s3对应的列有负值,在比值运算时直接赋值为空,因为比值只看正值,如果将负值也考虑进来取最小比值,同样将导致负的右端项。通过入基变量选取和比值测试,对元素2圈圈,做行列变换,得表4:
  表4 单纯形表3
  第六步:最新表中发现第0行的所有元素均为正值,此时选取任何变量入基,都会使得z值因为正的缩减成本而降低,很显然这对于最大化问题来说是不利的。因此,上表已经达到最优状态,单纯形法迭代结束。
  综上,原问题最优解就是
  本文主要介绍了单纯形法的基本逻辑思路,以及具体的求解过程,接下来我们将继续带领大家探索单纯形法求解过程中可能出现的几种解,以及单纯形法的变形求解方法。,希望大家继续关注【优化】板块,电子书线性规划专题的科普文。参考文献:
  [1] Dimitris Bertsimas, John N. Tsitsiklis, Introduction to Linear Optimization. Athena Scientific, Belmont, Massachusetts. P

魔兽世界转服多少钱(魔兽世界官方网站)魔兽世界转服多少钱(魔兽世界官方网站)服的小伙伴你们好,我是细细说,很高兴又跟大家见面了。曾经我不止一次分享过魔兽世界怀旧服转服赚金币的相关心得,但总感觉大部人听不太明白,所以今天杭州公交卡种类(杭州公交卡怎么开发票)杭州公交卡种类(杭州公交卡怎么开发票)杭区融媒体中心消息,在杭州,说到公交地铁出行,那就不得不提各类交通卡了。杭州的交通卡有几种类型?怎么办理?怎么充值?怎么挂失?月票是啥?小编送支付宝提现什么意思(支付宝提现手续费多少)支付宝提现什么意思(支付宝提现手续费多少)终于!继微信3月开始提现收费之后,支付宝这一天也到来了。9月12日,支付宝发布公布称,因为经营成本上升,为了给用户提供更多的优质服务,我们支付宝提现什么意思(支付宝提现手续费多少)支付宝提现什么意思(支付宝提现手续费多少)终于!继微信3月开始提现收费之后,支付宝这一天也到来了。9月12日,支付宝发布公布称,因为经营成本上升,为了给用户提供更多的优质服务,我们支付宝提现是什么意思(支付宝提现不收费技巧)支付宝提现是什么意思(支付宝提现不收费技巧)3月30日,支付宝开放平台发布公告,支付宝收钱码提现免费服务再延长3年,至2024年3月31日,且提现免费不设单笔上限和单日上限。支付宝余额宝提现要手续费吗(支付宝提现不收费技巧)余额宝提现要手续费吗(支付宝提现不收费技巧)现在的生活支付越来越离不开支付宝和微信了,即使不往里面存钱,也避免不了朋友通过微信或支付宝给你转账发红包。然而,如果需要将微信和支付宝里支付宝提现要多久(支付宝半夜提现要多久)支付宝提现要多久(支付宝半夜提现要多久)大家介绍一下(我)本人常用的免费提现到银行卡的方法方法一这个方法需要满足两个条件(需要用到两部手机),同一部手机切换账号系统会检测到。注册两支付宝怎样买东西(手机支付宝怎样网上买东西)支付宝怎样买东西(手机支付宝怎样网上买东西)南网如今在生活里消费购物的时候,使用花呗支付的人越来越多,花呗拥有着越来越多的用户,受到了大众的广泛喜爱。许多人都开始使用起支付宝的蚂蚁交话费送手机活动(2021移动充话费送手机)交话费送手机活动(2021移动充话费送手机)来源都市现场综合重庆新闻频道未经授权不得转载充话费就能送手机,而且每个月还有话费返还,这样算下来,不但手机白送还能赚两千多的手机话费。但移动充话费送手机活动(话费套餐送手机活动)移动充话费送手机活动(话费套餐送手机活动)我们经常在各种媒体上听见有什么商业模式的创新,我基本上都不会去看,不是因为我懂多少的商业,也不是因为有多么的自负,而是在商业上面要有什么创本地连接断开(本地连接已断开媒体连接)本地连接断开(本地连接已断开媒体连接)1。如果您安装过某个软件,可以先尝试卸载该软件,部分软件有可能导致此类问题的发生。2。若依然无法解决,需要重新手动开启服务NetworkSto
淘宝无法购买(淘宝无法购买东西怎么办)淘宝无法购买(淘宝无法购买东西怎么办)网友在某淘宝店不能下单,原因竟是之前在其他店面给过差评。另一名网友此前在淘宝一家店铺下单之后,因对方客服一直不在线,想给差评却没有找到评价页面哪里看油菜花(拍油菜花的姿势大全)哪里看油菜花(拍油菜花的姿势大全)青海在不经意间成为了油菜花观赏胜地,尤其是门源当然不仅仅是门源,青海湖周边广袤的高原上,到处都是油菜花。每到夏天青海湖边还有各个角落都开满了大片的油菜花是什么意思(油菜花的特点)油菜花是什么意思(油菜花的特点)油菜又叫芸苔,花期大约有三十天。受经纬度高度等影响,从每年的十二月底一直持续到来年的四月底都会有油菜花开。油菜花具有重要的经济价值,又有观赏价值,是卢本伟牛b是什么梗(卢本伟是什么梗全体起立)卢本伟牛b是什么梗(卢本伟是什么梗全体起立)什么梗?卢本伟牛逼出自55开经常开直播玩h1z1,只要杀到人了就说自己牛逼。卢本伟牛逼就在这样的情况下诞生了,现在在很多直播间都可以看到引号的作用有哪些(引号的作用有哪些举例)引号的作用有哪些(引号的作用有哪些举例)AndereLnder,andereSitten。这是当代大学德语课本中的一个标题,意为不同的国家有不同的习俗。而在这篇文章里,这句话也可以p型半导体和n型半导体(怎样区分N型和P型半导体)p型半导体和n型半导体(怎样区分N型和P型半导体)P型半导体即空穴浓度远大于自由电子浓度的杂质半导体N型半导体即自由电子浓度远大于空穴浓度的杂质半导体。P型和N型半导体都属于杂质半物理知识科普折射率公式三个公式关于到现在折射率公式三个公式这个话题,相信很多小伙伴都是非常有兴趣了解的吧,因为这个话题也是近期非常火热的,那么既然现在大家都想要知道折射率公式三个公式,小编也是到网上收集了一些与相对误差(相对误差计算公式)相对误差(相对误差计算公式)在测量过程中,由于测量器具和测量条件的限制,测量误差总是不可避免的。测量误差是指测得值与被测真值在数值上相差的程度,可用绝对误差和相对误差来表示。绝对误公顷与亩的换算(计量单位大全表)公顷与亩的换算(计量单位大全表)击右上方红色字体,关注农爱邦,获取农村最新热点资讯!作为一个在农村长大的孩子,说来惭愧,确实是不知道一亩地是多少平方米。那么,一亩地是多少平方米呢?公分和厘米的换算(厘米的换算单位和公式)公分和厘米的换算(厘米的换算单位和公式)参考答案890厘米(89)分米85千米(85000)米570厘米(57)分米460毫米(46)厘米900毫米(90)厘米81分米(810)厘数学知识负a分之b是求斜率吗关于到现在负a分之b是求斜率吗这个话题相信很多小伙伴都是非常有兴趣了解的吧因为这个话题也是近期非常火热的那么既然现在大家都想要知道负a分之b是求斜率吗小编也是到网上收集了一些与负a