前言
第1章 凸集和凸函數(shù) 1
1.1 最優(yōu)化問題和方法 1
1.2 凸集及相關性質(zhì) 4
1.3 凸集的分離 5
1.3.1 點和凸集的分離 5
1.3.2 凸集和凸集的分離 11
1.4 凸函數(shù) 13
1.4.1 凸函數(shù)及相關性質(zhì) 13
1.4.2 凸函數(shù)的判別 15
1.5 凸規(guī)劃問題 20
1.6 習題 21
第2章 線性規(guī)劃的基本性質(zhì) 23
2.1 線性規(guī)劃的形式 23
2.1.1 線性規(guī)劃的三種基本形式 23
2.1.2 三種形式相互等價 29
2.2 可行域和頂點 30
2.3 解線性規(guī)劃的圖解法 32
2.4 基本解和基本可行解 33
2.5 線性規(guī)劃基本定理 41
2.6 習題 43
第3章 單純形算法 45
3.1 單純形算法的基本思想 45
3.2 幾何形式的單純形算法 46
3.3 代數(shù)化的單純形算法 46
3.3.1 基本思想 46
3.3.2 代數(shù)化的單純形算法示例 46
3.4 一般的單純形算法 51
3.4.1 檢驗數(shù)向量 51
3.4.2 目標函數(shù)值和檢驗數(shù)向量的值 51
3.4.3 單純形算法 53
3.5 表格化的單純形算法 54
3.5.1 單純形表 54
3.5.2 旋轉(zhuǎn) 55
3.6 使用數(shù)學軟件解線性規(guī)劃 62
3.6.1 使用MATLAB解線性規(guī)劃 63
3.6.2 使用CPLEX解線性規(guī)劃 64
3.7 單純形算法的分析 66
3.8 退化問題的處理 71
3.9 兩階段法 72
3.9.1 兩階段法的基本思想 72
3.9.2 解輔助線性規(guī)劃 73
3.9.3 兩階段單純形算法 75
3.10 矩陣的全單位模性質(zhì) 85
3.11 再議解線性規(guī)劃 87
3.11.1 單純形算法的復雜性 87
3.11.2 解線性規(guī)劃的多項式時間算法 88
3.11.3 單純形算法的平滑分析 88
3.12 習題 89
第4章 線性規(guī)劃對偶理論 92
4.1 線性規(guī)劃的對偶 92
4.2 對偶定理 98
4.3 對偶單純形算法 107
4.4 關于單純形表檢驗數(shù)行和右端項的討論 112
4.5 原始對偶算法 113
4.5.1 最短路問題的整數(shù)規(guī)劃 114
4.5.2 原始對偶算法 115
4.5.3 算法分析 117
4.6 習題 121
第5章 整數(shù)規(guī)劃 124
5.1 整數(shù)規(guī)劃問題 124
5.1.1 背包問題 124
5.1.2 最小生成樹問題 125
5.1.3 旅行售貨商問題 126
5.1.4 整數(shù)線性規(guī)劃 127
5.2 割平面法 129
5.2.1 割平面法的基本思想 129
5.2.2 割平面的生成方法 129
5.3 分枝定界法 135
5.3.1 分枝定界法的基本思想 135
5.3.2 分枝定界法解整數(shù)規(guī)劃 136
5.4 習題 143
第6章 動態(tài)規(guī)劃 144
6.1 動態(tài)規(guī)劃的原理 144
6.1.1 多階段決策問題 144
6.1.2 最優(yōu)化原理 149
6.1.3 前向優(yōu)化和后向優(yōu)化 151
6.2 問題舉例 152
6.2.1 最長公共子序列問題 152
6.2.2 背包問題 154
6.2.3 從背包問題談時間復雜度 157
6.2.4 旅行售貨商問題 160
6.2.5 一般圖上的最短s-t路問題 164
6.3 習題 167
第7章 圖與網(wǎng)絡算法 169
7.1 最大流問題 169
7.1.1 最大流的增廣路算法 169
7.1.2 最大流和最小割 178
7.1.3 對偶理論的觀點 179
7.2 最小費用流問題 183
7.3 匹配問題概述 193
7.4 二分圖上不帶權(quán)重的最大匹配問題 194
7.4.1 使用最大流算法求解 194
7.4.2 增廣路算法 197
7.5 二分圖上帶權(quán)重的最大匹配問題 204
7.5.1 歸約到最小費用流問題的解法 204
7.5.2 匈牙利算法 206
7.6 一般圖上的最大匹配問題 211
7.7 習題 211
第8章 無約束優(yōu)化的基本概念 214
8.1 一元函數(shù)的極小化問題 214
8.1.1 黃金分割法 214
8.1.2 函數(shù)逼近法 218
8.2 下降方向 218
8.3 一維搜索的基本概念 219
8.4 習題 220
第9章 使用導數(shù)的無約束優(yōu)化方法 221
9.1 無約束優(yōu)化問題的一階極值條件 221
9.2 下降算法的一般形式 222
9.3 最速下降法 223
9.3.1 算法 223
9.3.2 搜索為什么要沿負梯度方向進行 223
9.3.3 最速下降法的鋸齒現(xiàn)象 227
9.4 牛頓法 228
9.4.1 一元優(yōu)化問題的牛頓法 228
9.4.2 多元優(yōu)化問題的牛頓法 230
9.4.3 阻尼牛頓法 234
9.4.4 牛頓法的進一步修正 237
9.5 共軛梯度法 238
9.5.1 基本概念 238
9.5.2 共軛方向的幾何意義 239
9.5.3 共軛梯度算法 241
9.5.4 一般無約束優(yōu)化問題的共軛梯度法 250
9.6 無約束優(yōu)化問題的二階極值條件 251
9.7 擬牛頓法 253
9.7.1 擬牛頓方程 253
9.7.2 DFP算法 254
9.7.3 BFGS算法 256
9.8 習題 256
第10章 約束優(yōu)化問題的基本概念和性質(zhì) 258
10.1 問題舉例 258
10.2 可行方向 260
10.3 不等式約束問題的一階最優(yōu)性條件 261
10.3.1 必要條件 262
10.3.2 充分條件 272
10.4 一般約束問題的一階最優(yōu)性條件 273
10.4.1 必要條件 273
10.4.2 充分條件 275
10.5 約束優(yōu)化問題的對偶理論 279
10.5.1 對偶問題 279
10.5.2 凸規(guī)劃的對偶 283
10.6 習題 288
第11章 約束優(yōu)化問題的解法 290
11.1 二次規(guī)劃的解法 290
11.1.1 等式約束二次規(guī)劃的直接消元法 290
11.1.2 等式約束二次規(guī)劃的拉格朗日方法 292
11.1.3 一種凸二次規(guī)劃的有效集方法 295
11.2 簡約梯度法 304
11.2.1 簡約梯度 304
11.2.2 構(gòu)造搜索方向 308
11.2.3 算法設計 312
11.3 罰函數(shù)法 321
11.3.1 外點罰函數(shù)法 321
11.3.2 內(nèi)點罰函數(shù)法 327
11.4 習題 330
第12章 若干基本的數(shù)學概念和定理 332
12.1 n維空間中的點與集合 332
12.2 連續(xù)和一致連續(xù) 333
12.3 無窮小量與無窮小量的階 333
12.4 一元函數(shù)可微和微分 334
12.5 二元函數(shù)可微和全微分 335
12.6 泰勒公式 336
12.6.1 帶佩亞諾余項的泰勒公式 336
12.6.2 帶拉格朗日余項的泰勒公式 336
12.7 方向?qū)?shù) 337
參考文獻 339
索引 343
第1章凸集和凸函數(shù)
本章首先簡要介紹最優(yōu)化方法這門學科,然后著重介紹凸集和凸函數(shù)這兩個基本概念及其相關的性質(zhì).凸集和凸函數(shù)在最優(yōu)化方法中廣泛出現(xiàn)在線性規(guī)劃、無約束優(yōu)化以及約束優(yōu)化等問題的描述和求解方法中.
1.1最優(yōu)化問題和方法
最優(yōu)化方法學科研究優(yōu)化問題的求解方法.現(xiàn)實生活中,優(yōu)化問題多種多樣,最優(yōu)化方法主要研究那些可以用數(shù)學模型刻畫的優(yōu)化問題.這些問題從總體上可分為組合優(yōu)化問題和連續(xù)優(yōu)化問題兩類,它們是按照數(shù)學模型中所使用的變量的不同類型來劃分的.
在歷史上,優(yōu)化技術(shù)出現(xiàn)得很早.17世紀時,在微積分技術(shù)的發(fā)展過程中,牛頓就研究過極值問題.后來又出現(xiàn)了處理帶約束的目標函數(shù)極值問題的拉格朗日(Joseph-Louis Lagrange,法國,1736-1813)乘數(shù)法 1847年,柯西(Augustin-LouisCauchy,法國,1789-1857)研究了函數(shù)值沿什么方向下降最快的問題,提出了最速下降法.1939年,康托羅維奇(Kantorovich,蘇聯(lián),1912-1986)提出了下料問題和運輸問題的求解方法.直至20世紀30年代,最優(yōu)化這個古老的課題尚未形成獨立的學科.
20世紀40年代以來,由于工業(yè)生產(chǎn)和科學研究的迅猛發(fā)展,特別是計算機的日益廣泛應用,最優(yōu)化理論和算法迅速發(fā)展起來,成為運籌學的主體內(nèi)容.運籌學(最優(yōu)化方法)中的組合優(yōu)化部分和計算機科學中的算法設計與分析部分是部分重合的,由于組合優(yōu)化技術(shù)更強調(diào)了數(shù)學方法的運用,因此組合優(yōu)化可以看作是傳統(tǒng)的算法設計與分析的延伸.運籌學中的連續(xù)優(yōu)化部分是在數(shù)學分析中的求導數(shù)和求極值的基礎上發(fā)展起來的.現(xiàn)在,機器學習中損失函數(shù)的求解正好是連續(xù)優(yōu)化問題.因此,運籌學的連續(xù)優(yōu)化部分也構(gòu)成了應用計算機科學解決實際問題的基礎.
最優(yōu)化方法的組合優(yōu)化技術(shù)主要包括整數(shù)線性規(guī)劃、動態(tài)規(guī)劃以及面向問題的各種啟發(fā)式方法等.最優(yōu)化方法中的連續(xù)優(yōu)化技術(shù)可分為線性規(guī)劃和非線性規(guī)劃,也可分為無約束優(yōu)化和約束優(yōu)化,這是分別按照不同的劃分標準劃分的.值得指出的是,線性規(guī)劃技術(shù),按照變量的類型劃分屬于連續(xù)優(yōu)化,但線性規(guī)劃技術(shù)也可以用于求解組合優(yōu)化問題.例如在近似算法中,線性規(guī)劃是獲得廣泛成功的一種算法設計方法.使用線性規(guī)劃技術(shù)設計組合優(yōu)化問題的近似算法時,人們首先用線性規(guī)劃描述組合優(yōu)化問題的松弛,然后求解線性規(guī)劃得到最優(yōu)解,再通過舍入等方法將線性規(guī)劃最優(yōu)解轉(zhuǎn)化為組合優(yōu)化問題的可行解.
下面我們看幾個典型的優(yōu)化問題.
例1.1.1利潤最大化問題.某工廠用種原料生產(chǎn)種產(chǎn)品.第種原料的數(shù)量為每單位的第種產(chǎn)品需要第種原料的數(shù)量為,可獲利潤為.問如何安排生產(chǎn),才能使總利潤最大?
例1.1.1中的利潤最大化問題是經(jīng)濟領域中經(jīng)常遇到的一個非;镜膯栴}.這個問題可以用線性規(guī)劃或整數(shù)線性規(guī)劃進行描述和求解.
在這里,我們定義為第種產(chǎn)品的產(chǎn)量.則所有種產(chǎn)品的總利潤就是.問題的目標是最大化這個總利潤.當然,利潤不可能無限大,因為安排生產(chǎn)受限于原材料的供給.當?shù)诜N產(chǎn)品的產(chǎn)量為時,生產(chǎn)這種產(chǎn)品對原料的需求為,顯然應有,因為第種原料總共有那么多.另外,每種產(chǎn)品的產(chǎn)量也不能為負數(shù).綜合起來,我們就得到
(1.1.1)
這就是例1.1.1的數(shù)學模型,這是一個線性規(guī)劃,因為在這個規(guī)劃中目標函數(shù)和約束條件都只包含變量的一次項,沒有變量的二次項以及更高次項,更沒有變量的冪以外的其他函數(shù).
值得指出的是,線性規(guī)劃(1.1.1)適用于描述那些產(chǎn)品“無限可分”的利潤最大化問題,如產(chǎn)品為面粉、食用油等.換言之,產(chǎn)品數(shù)量取值為實數(shù)應該是有意義的.
如果在利潤最大化問題中,產(chǎn)品數(shù)量僅能取值為整數(shù),例如產(chǎn)品為家具、門窗等,則我們需要使用整數(shù)線性規(guī)劃為利潤最大化問題建模,如線性規(guī)劃(1.1.2)所示.
(1.1.2)
1.1最優(yōu)化問題和方法
定義1.1.2設施選址問題(FacilityLocation Problem)
實例:有個設施和個客戶,設施和客戶之間,以及設施之間、客戶之間都有距離,這些距離滿足三角不等式.打開第個設施有一個費用方.每個客戶都需要連接到一個打開的設施上.
目標:打開若干設施,滿足客戶的連通需求,使得設施的打開費用和客戶的連接費用之和最小.
設施選址問題是組合優(yōu)化領域中一個著名的最小優(yōu)化問題.它描述了這樣一種場景:比如有一個工廠,它有分布在多處地點的零售點.現(xiàn)在這個工廠要建立若干倉庫,以滿足向零售點配送貨物的需要.建設倉庫就有一個建設費用,倉庫建好之后,零售點就會選擇距離最近的倉庫供貨.如何在若干的候選位置建設倉庫,以滿足零售點的供貨需求,就表達為設施選址問題.在這里,倉庫的候選位置就是設施位置,零售點就是客戶.設施的打開費用就是建設倉庫的費用,客戶的連接費用就是所有零售點到其最近的倉庫的距離之和.
在定義1.1.2中,一個問題是按照實例和目標兩部分敘述的.這是在計算機科學中對問題描述的一種常見的形式.定義1.1.2中的“實例”也可以叫做“輸入”,“目標”也可以叫做“輸出當問題是判定問題時(即,回答“是”和“否”的問題),“目標”也可以寫作“詢問.
關于設施選址問題,“從m個設施中選擇若干設施打開”和“從m個設施位置中選擇若干位置建立設施”這兩種說法是等價的.
在這里,我們繼續(xù)用線性規(guī)劃對設施選址問題進行建模,以領略線性規(guī)劃強大的表達能力.
定義變量表示是否打開設施,若,表示不打開設施,若,則表示打開設施定義變量表示客戶是否連接到設施上,它的值也是或者.這樣,問題的總費用就可以表達為里表示頂點和之間的距離.客戶要連接到設施上,必須設施打開的條件下才行,這可以用來表達.另外,每個客戶一定要連接到某個設施上,這可表達為.因此,設施選址問題的線性規(guī)劃模型為.
這個線性規(guī)劃是一個整數(shù)線性規(guī)劃,因為它的每個變量只能取整數(shù)值.
設施選擇問題是一個基本的NP困難問題.設施選址問題還有很多的變形.對于設施選址問題及其變形,學者們設計了很多的近似算法,其中有若干算法是非常有名的.感興趣的讀者可進一步參考等著作及其引用的論文.
例1.1.3機器學習中的函數(shù)優(yōu)化問題.神經(jīng)網(wǎng)絡是機器學習中的一種模型,它和優(yōu)化問題密切相關.多層感知機(Multi-Layer Perceptron,MLP)是一種由若干模擬的神經(jīng)元排列成矩陣結(jié)構(gòu)的前饋神經(jīng)網(wǎng)絡,它所完成的基本任務是進行分類.可以使用一個函數(shù)來表達多層感知機所完成的任務.這里,向量是輸入數(shù)據(jù),向量:是多層感知機使用的參數(shù),是多層感知機的輸出,其含義是把輸入數(shù)據(jù)在參數(shù)最的控制下分類成所表達的那種類型.例如,在實際應用中,可將圖像數(shù)據(jù)輸入給神經(jīng)網(wǎng)絡,讓它將圖像分成若干種類型.
為了使神經(jīng)網(wǎng)絡能夠以盡量高的準確率完成分類,需要設定參數(shù)最的合適的值,這是在訓練數(shù)據(jù)的幫助下完成的,即人們通常說的“訓練一個神經(jīng)網(wǎng)絡訓練集可表達為,即,已經(jīng)知道數(shù)據(jù),的類型為,也稱為數(shù)據(jù)的標簽.訓練一個神經(jīng)網(wǎng)絡,就是調(diào)整參數(shù),使得某種損失函數(shù)的值盡量地小.若選擇平方誤差為損失函數(shù),則多層感知機的優(yōu)化模型為.
其中,表示范數(shù),其定義為按照最優(yōu)化方法的觀點,這是一個無約束優(yōu)化問題.在這里,為了簡單說明問題,我們將優(yōu)化模型簡化了.實際使用的優(yōu)化模型還要加上一個正則項來進行修正.這些細節(jié)就略去了.
在足夠強大的訓練集的支持下將一個神經(jīng)網(wǎng)絡訓練好(即,選擇好了參數(shù)a:的值),就可以在應用中使用這個神經(jīng)網(wǎng)絡完成數(shù)據(jù)分類任務了.
1.2凸集及相關性質(zhì)
定義1.2.1
兩個點,的凸組合還是中的一個點.它們的所有凸組合構(gòu)成以這兩個點為端點的一條線段.
定義I.2.2
定義1.2.3
例如,歐氏空間中三角形區(qū)域、矩形區(qū)域、圓形區(qū)域都是凸集.在幾何直觀上,凸集就是邊界“向外鼓”的集合.
例1.2.4
例1.2.5
定義1.2.6
在幾何直觀上,三角形區(qū)域的頂點、矩形區(qū)域的頂點都是符合定義1.2.6的頂點.讀者可領略定義1.2.6的精妙.注意,圓形區(qū)域是凸集,它有無數(shù)個頂點.
凸集具有良好的性質(zhì),人們對凸集有比較深入的把握.當一個數(shù)學規(guī)劃的可行解的集合是一個凸集時,意味著這個數(shù)學規(guī)劃可能會比較容易求解.例如,線性規(guī)劃的可行域是凸集,當一個線性規(guī)劃有解時,它的最優(yōu)解可在頂點上取得.這些知識將在第2章中詳細介紹.
1.3凸集的分離
凸集是點的集合,不在一個凸集中的點,在直觀上,位于這個凸集的“外面我們?nèi)绾问褂脭?shù)學的語言來精確描述這種直觀?進而,兩個凸集若不相交,它們是“分離”開的,這樣的直觀用數(shù)學語言如何表達?令人驚奇的是,這種看上去非;镜闹庇^概念,卻蘊含著深刻的道理.例如,由點和凸集的分離定理1.3.3,可以推導出著名的Farkas引理(引理1.3.4)(G.Farkas,匈牙利,1847-1930),而Farkas引理可推導出線性規(guī)劃的強對偶定理(定理4.2.5).再如,由凸集和凸集的分離定理1.3.9可推導出戈丹引理(引理1.3.10)(P.A.Gordan,德國,1837—1912),而戈丹引理又可推導出關于約束優(yōu)化的著名的K-T定理(定理10.3.7).
1.3.1點和凸集的分離
在介紹定理1.3.2之前,先來回顧一下下確界符號inf.
定義1.3.1令SgR是一個集合,若3m,使得Vze民都有;則稱m是S的一個下界.記L是S的所有下界所組成的集合.則L中的最大元素a稱為S的下確界,記為
定理1.3.2
證明
則f(x)是連續(xù)函數(shù).
若S是有界集,則是非空、有界且閉的,所以函數(shù)在上可以達到最小值,即,使得
即,
現(xiàn)在假設不是有界集.在中任取一點以為球心,為半徑做一個球,如圖1.3.1所示.
……