LP问题进阶 Part 2 | 对偶理论

这一部分对应书上第七章(P45-P51),难度和第四章差不多。由于引用了一些非数学的概念,所以学习的时候难免会遇到一些看似无所根据的概念或者假设。建议大家不要在各种花里胡哨的假设上浪费时间,要把精力留给数学。
虽然线性规划的对偶理论在线性代数和数学模型中都没怎么被提到,但其理论本身有趣而且有用。对偶理论将原LP问题转化为其对偶LP问题,从而使与原问题相关的一些复杂的性质(如解是否为最优)变为对偶问题中较为简单的性质(对应地,解是否可行)
第五章和第六章暂时没有讲。 第五章主要是线性代数的内容,在此默认大家都学过。需要看一眼的是最后提到的用矩阵表示LP问题的方法。第六章是敏感性分析,我打算留到最后和算法复杂度一起讲。

为方便查阅,再link一下教材


基本概念

为方便理解我们引入一个例子:

例1: 现有一工厂可以生产小人偶和小火车这两种玩具。生产一个小人偶需要1单位木材和2单位颜料,产生3单位收益。生产一辆小火车需要1单位木材和1单位颜料,产生2单位收益。工厂现有80单位木材和100单位颜料,问如何分配火车与人偶的生产量可以使得收益最大?(见下图左侧)

为方便理解对偶性而引入的一个简单的例子

先补充几个前几章出现过但我没有提到的简单概念。由于没有了解过中文版教材所以翻译可能不太准确。

  • 对象(Item)是LP问题中的“名词”,即收益材料。在此问题中,对象有三个,即:收益、木材、颜料。对象对应着LP问题中的行。在此问题中,收益、木材与颜料分别对应着第一到第三行。
  • 行为(Activity)是LP问题中的“动词”,对应着标准表示中的列。通俗地讲,行为就是将材料转化为产品(收益)的过程。在本问题中,制造小火车与制造小人偶就是两种行为。上图中,x_{1}所在的列对应着制造小人偶的行为,x_{2}所在的列对应着制造小火车的行为。

回忆上一次在基本概念中讲到的基本解与字典的概念,结合线性代数的知识,我们可以将LP问题写为:

LP问题的矩阵表示

不要被上图中复杂的形式迷惑,其实推导过程非常简单。式中各个符号的定义见教材第32页。需要强调的两点是:

  • x_N是由非基本变量构成的向量,在具体取值的时候即为0向量。详见之前对非基本变量的定义。注意不要将非基本变量理解为0变量。非基本变量只是在考虑对应基本解时取值为0。
  • 此处假设B是可逆的,但我无法证明B一定是可逆的。但我相信大多数情况下B都是可逆的,因为在字典中x_{B}可有x_{N}表示,而联系二者的只有等式Ax = b。如果B不满秩的话会丢失x_{B}的信息进而使得x_{B}无法被x_{N}表示。(直观感受,未证明)
  • 基本解对应的目标函数的取值为c_{B}B^{-1}b,当c_{N} - c_{B}B^{-1}N的各项系数小于0时取最优解。此时基本变量的取值为B^{-1}b

目标函数的值(c_{B}B^{-1}b)和判断最优解的条件(c_{N} - c_{B}B^{-1}N < 0)这两个矩阵表示非常重要,建议各位同学自己推导一遍并记下来。

现在假设在市场上有以价格y_{1}y_{2}售卖的木材与颜料,如上图右侧所示。市场是一个很有趣的概念,可惜我不怎么懂,只知道本章许多内容都与之多多少少地相关。接下来介绍一个比较重要的概念。

  • 影子价格(百度百科): 某种材料的影子价格通俗地讲就是该种材料增加1单位所带来的额外的收益
    当某种材料过剩时,其影子价格即为0,因为该材料本来就多于所需,所以它的增加不带来额外收益。 关于影子价格,需要注意以下两点:
    • 影子价格的定义严格地讲是收益关于某种材料的变化率,但由于我们考虑的是线性问题所以干脆说成是“由于增添1单位某材料所带来的额外收益”。影子价格是一个局部的概念,不要钻牛角尖。
    • 书上没有任何解释地给出了影子价格的矩阵表示:c_{B}B^{-1}。其实这就是一个简单的求导运算。对比目标函数的值:c_{B}B^{-1}b可以发现某种材料的影子价格恰好就是目标函数关于该材料的导数。需要注意的是此时由于材料可以自由买卖,所以b变成了一个变向量。在某种意义上,目标函数其实本来就是材料的函数,只不过材料的多少会影响函数的形式(B^{-1}也会随材料改变。)

影子价格c_{B}B^{-1}是一个向量,其每个元素对应一种材料。

基本概念大概就是这些,接下来就是神奇的部分了。


对偶性理论

这一部分有许多条件并不是一般的,比如我们只讨论了最大值问题而没研究最小值问题。不过方法的本质都是一样的。

1. 基本假设与直观理解

对偶(Dual)不论是单词还是翻译都和泛函分析中函数空间的那种对偶一样。事实上对偶是数学中一种非常常见的思想方法。通过配合地研究原问题(Primal)对偶问题(Dual)我们可以得到很多性质颇好的结论。首先,我们引入一个假设:

假设1: 生产一件产品所需的材料的价格不小于出售该件产品所带来的收益。

之后我们简称生产某件产品所需的材料的价格为该产品的材料成本。对于该假设我们可以作如下理解:

  • 首先,我们引入的所有假设的目的都是为了解原LP问题,故引入新的变量不能改变原LP问题的最大值。只对数学模型感兴趣的同学可以直接跳过之后的三点。
  • 假设在此问题中,工厂本来就有一定数量的各种材料。生产一件产品的成本不包括材料成本
  • 在此假设下,购入更多的原材料只会使总收益变少,所以原LP问题的最大值不变。
  • 如果该假设不成立,那么该工厂就可以从市场中无限地买入原材料进而获得无限的利益。由于市场存在竞争,这种情况是不允许出现的。

其实在思考本假设的意义时我也有很多疑问,当然可能我现在对于此问题的理解仍然是错误的。我们不妨先默认这些看似不合理直观的假设,推得我们想要的理论,之后再严谨地用数学方法证明理论的合理性即可。当然,感兴趣的同学也可以参考经济学领域的相关书籍。
不论如何,我们不应在直观理解上浪费太多时间。我们现在引出原问题对偶问题的形式:

原问题与对偶问题

左边为原问题,右边为对偶问题。首先可以注意到的是对偶问题的限制条件即假设1。此时有显然的关系式:3x_{1}+2x_{2} \leq (y_{1} + 2 y_{2})x_{1} + (y_{1} + y_{2})x_{2} = y_{1}(x_{1} + x_{2}) + y_{2}(2x_{1}+x_{2}) \leq 80y_{1} + 100y_{2}左边的不等号对应于对偶问题的限制条件,右边的不等号则对应原问题的对偶条件。因此两边目标函数存在关系:Max \ 3x_{1}+2x_{2} \leq Min \ 80y_{1} + 100y_{2}是显然的。这其实就是待会要提到的弱对偶性

原问题与对偶问题的矩阵表示如下图所示:

原问题与对偶问题的矩阵表示

由此可知对偶问题的对偶问题就是原问题。即两问题互为对偶

2. 弱对偶性与强对偶性定理

关于对偶性,最关键的理论就是以下两个重要的定理:

定理1 (弱对偶性) 假设xy分别是原问题与对偶问题的可行解,则:c^{T}x \leq b^{T}y证明:由限制条件可知:c^{T}x \leq (A^{T}y)^{T}x = (y^{T}A)x = y^{T}(Ax) \leq y^{T}b = b^{T}y证明完毕。

由弱对偶性可知:

  • 若原问题的目标函数无界,则对偶问题必无可行解;同理若对偶问题的目标函数无界,则原问题必无可行解。注意前者的无界是无上界而后者的无界是无下界。若产品的收益都为正,则目标函数无界等价于可行域无界。
  • 书上说原问题和对偶问题都无可行解的情况是存在的。不过我还没有想到这样的例子。

定理2 (强对偶性) 假设xy分别是原问题与对偶问题的最优解,则:c^{T}x = b^{T}y 且满足y^{T}(b-Ax) = x^{T}(A^{T}y - c) = 0其中b-AxA^{T}y - c分别对应着原问题与对偶问题中的松弛变量。

我们现在还无法证明这个定理。

3. 原问题与对偶问题的性质

为了证明强对偶性定理并得到一些我们需要的性质,在此先叙述一些相对简单的概念与结论。由于原问题与对偶问题实际上互为对偶,故不失一般性地可以只叙述原问题的性质。

性质1: 原问题中每一个等号-限制条件都对应一个无符号限制的对偶变量。

此处“等号-限制条件”的含义应该是清晰的。原问题中的每一条限制条件都对应着一个对偶变量,这一点大家也需要牢记。强调以下几点:

  • 如果原问题中采用的是\leq-限制条件,则无符号限制的对应对偶变量会导致弱对偶性定理失效。因为负号会使不等式转向。对于\geq-限制条件也一样。现考虑例1中木材的限制条件对应的对偶变量(即木材的价格y_{1}):3x_{1}+2x_{2} \leq (y_{1} + 2 y_{2})x_{1} + (y_{1} + y_{2})x_{2} = y_{1}(x_{1} + x_{2}) + y_{2}(2x_{1}+x_{2})如果不要求y_1 \geq 0的话显然无法推出3x_{1}+2x_{2} \leq 80y_{1} + 100y_{2}
  • 如果原问题中的某条限制条件是等号-限制条件,则称该限制条件为紧(tight)的。这个概念之后还会用到。
  • 若原问题中某条限制是等号-限制条件则不论对应的对偶变量如何取值,弱对偶性仍然成立。依然考虑例1中的木材,但此时将其限制条件改为等号-限制条件,则有:3x_{1}+2x_{2} \leq (y_{1} + 2 y_{2})x_{1} + (y_{1} + y_{2})x_{2} = y_{1}(x_{1} + x_{2}) + y_{2}(2x_{1}+x_{2}) = 80 y_{1} + y_{2}(2x_{1}+x_{2}) \leq 80 y_{1} + 100 y_{2}即仍满足弱对偶性。
  • 类比地可以总结出规律:对于一个求最大值的LP问题,\leq-限制条件对应非负对偶变量;\geq-限制条件对应非正对偶变量;等号-限制条件对应无符号限制的对偶变量

接下来介绍互补松弛的概念。

定义1:称向量x = (x_{1}, \cdots , x_{n})y = (y_{1}, \cdots , y_{n})互补的(complementary),若:y^{T}(b-Ax) = x^{T}(A^{T}y - c) = 0

关于互补,有如下性质:

  • 强对偶性说明原问题的最优解与对偶问题的解是互补的。
  • 对于互补的向量xy,若y_{i} > 0则原问题的第i个限制条件为紧的。同理若x_{i} > 0则对偶问题的第i个限制条件是紧的。
  • 对于任意基本解x,由其定义的影子价格\pi与之互补。证明如下:

命题1:对于LP问题的任意基本解x,由其定义的影子价格\pi与之互补
证明:明天再说!

关于互补松弛性,有以下几点重要的性质。先不加证明地列出:

x是原问题的最优解,则

  • y是对偶问题的最优解,则xy互补。(强对偶性)
  • y是对偶问题的可行解且xy互补,则y是对偶问题的最优解。
  • 存在对偶问题的可行解y使得xy互补。

结合命题1可知:

x是原问题的基本可行解,\pi是其影子价格,则x为最优解当且仅当\pi是可行解。

如本文开头所言,判断最优性的问题被巧妙地转化为了判断可行性地问题。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 216,372评论 6 498
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,368评论 3 392
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 162,415评论 0 353
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,157评论 1 292
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,171评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,125评论 1 297
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,028评论 3 417
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,887评论 0 274
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,310评论 1 310
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,533评论 2 332
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,690评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,411评论 5 343
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,004评论 3 325
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,659评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,812评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,693评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,577评论 2 353

推荐阅读更多精彩内容

  • 这一部分对应书上第四章(P16-P26),其难度比前三章大得多。个人建议先粗读了解这个方法的思路,再精读把握具体技...
    StRygwyr阅读 2,353评论 0 1
  • ​一般来说凸优化(Convex Optimization, CO)中最一般的是锥规划 (Cone Programm...
    史春奇阅读 5,074评论 1 6
  • 人生的大多数时刻是要浑浑噩噩过的 但小部分时间也需要睁眼看看世界
    不知道起什么比较好呀阅读 109评论 0 0
  • 唐浩,湖南湘江新区发展集团有限公司,299期努力一组同修,湖南国学践行研修班第22期仁组学员,23期良组、24期和...
    大肚哥_be8c阅读 89评论 0 0
  • 今天失眠了,大半夜睡不着,索性干脆,更新日记吧。 昨天,看了,我们是真正的朋友。里面小s t被信说。她没有很认真的...
    翼宁小宁子阅读 137评论 0 0