线性规划

在數學中,線性規劃(英語:Linear Programming,簡稱LP)特指目標函數約束條件皆為線性最佳化問題。

线性规划

線性規劃是最優化問題中的一個重要領域。在作業研究中所面臨的許多實際問題都可以用線性規劃來處理,特別是某些特殊情況,例如:網路流、多商品流量等問題,都被認為非常重要。目前已有大量針對線性規劃演算法的研究。很多最佳化問題算法都可以分解為線性規劃子問題,然後逐一求解。在線性規劃的歷史發展過程中所衍伸出的諸多概念,建立了最優化理論的核心思維,例如「對偶」、「分解」、「凸集」的重要性及其一般化等。在微观经济学和商业管理领域中,线性规划亦被大量应用于例如降低生产过程的成本等手段,最終提升產值與營收。對線性規劃有早期貢獻的列昂尼德·维塔利耶维奇·康托罗维奇特亚林·科普曼斯於1975年共同獲得諾貝爾經濟學獎。

标准型

描述线性规划问题的常用和最直观形式是标准型。标准型包括以下三个部分:

  • 一个需要极大化的线性函数,例如
 
  • 以下形式的问题约束,例如:
 
 
 
  • 和非负变量,例如:
 
 

线性规划问题通常可以用矩阵形式表达成:

maximize  
subject to  

其他类型的问题,例如极小化问题,不同形式的约束问题,和有负变量的问题,都可以改写成其等价问题的标准型。

例子

以下是一個線性規劃的例子。假設一個農夫有一塊 平方千米的農地,打算種植小麥大麥,或是兩者依某一比例混合種植。該農夫只可以使用有限數量的肥料 農藥 ,而單位面積的小麥和大麥都需要不同數量的肥料和農藥,小麥以 表示,大麥以 表示。設小麥和大麥的售出價格分別為  ,則小麥與大麥的種植面積問題可以表示為以下的線性規劃問題:

  (最大化利潤 - 目標函數)
    (種植面積的限制)
  (肥料數量的限制)
  (農藥數量的限制)
  (不可以栽種負數的面積)

增广矩阵(松弛型)

在用单纯型法求解线性规划问题之前,必须先把线性规划问题转换成增广矩阵形式。增广矩阵形式引入非负松弛变量英语Slack variable将不等式约束变成等式约束。问题就可以写成以下形式:

Maximize   in:
 
 

这里 是新引入的松弛变量,  需要极大化的变量。

例子

以上例子的转换成增广矩阵:

maximize   (目标函数)
subject to   (限制式)
  (限制式)
  (限制式)
 

这里 ,是(非负)松弛变量。

写成矩阵形式:

Maximize   in:
 

对偶

每个线性规划问题,称为原问题,都可以变换为一个对偶问题。可将“原问题”表达成矩阵形式:

maximize  
subject to  

而相应的对偶问题就可以表达成以下矩阵形式:

minimize  
subject to  

这里用 来作为未知向量。

例子

上述线性规划例子的对偶问题:

假如有一个种植园主缺少肥料和农药,他希望同这个农夫谈判付给农夫肥料和农药的价格。可以构造一个数学模型来研究如何既使得农夫觉得有利可图肯把肥料和农药的资源卖给他,又使得自己支付的金额最少?

问题可以表述如下

假设 分别表示每单位肥料和农药的价格,则所支付租金最小的目标函数可以表示为

 
 
  (控制肥料與農藥的價格,使得农夫觉得比起拿那些肥料和農藥去種植小麥,賣給園主更有利可图)
  (與上相似,但改為大麥)
  (不可用负数单位金額购买)

理論

幾何上,線性約束條件的集合相當於一個凸包或凸集,叫做可行域。因為目標函數亦是線性的,所以其極值點會自動成為最值點。線性目標函數亦暗示其最優解只會出現在其可行域的邊界點中。

在兩種情況下線性規劃問題沒有最優解。其中一種是在約束條件相互矛盾的情況下(例如  ),其可行域將會變成空集,問題沒有解,因此亦沒有最優解。在這種情況下,該線性規劃問題會被稱之為「不可行」。

另一種情況是,約束條件的多面體可以在目標函數的方向無界(例如: ),目標函數可以取得任意大的數值,所以沒有最優解。

除了以上兩種病態的情況以外(問題通常都會受到資源的限制,如上面的例子),最優解永遠都能夠在多面體的頂點中取得。但最優解未必是唯一的:有可能出現一組最優解,覆蓋多面體的一條邊、一個面、甚至是整個多面體(最後一種情況會在目標函數只能等於0的情況下出現)。

演算法

 
兩個變量的線性規劃問題中,一組線性約束條件劃定了對兩個變量的解的可行域。可解的問題會有一個簡單多邊形的可行域。

單純形演算法利用多面體的頂點構造一個可能的解,然後沿著多面體的邊走到目標函數值更高的另一個頂點,直至到達最優解為止。雖然這個演算法在實際上很有效率,在小心處理可能出現的「迴圈」的情況下,可以保證找到最優解,但它的最壞情況可以很壞:可以構築一個線性規劃問題,單純形演算法需要問題大小的指數倍的運行時間才能將之解出。事實上,有一段時期內人們曾不能確定線性規劃問題是NP完全問題還是可以在多項式時間裏解出的問題。

第一個在最壞情況具有多項式時間複雜度的線性規劃算法在1979年由前蘇聯數學家列昂尼德·哈奇揚英语Leonid Khachiyan提出。這個算法建基於非線性規劃瑙姆·Z·索爾英语Naum Z. Shor發明的橢球法(ellip-soid method),該法又是阿爾卡迪·內米羅夫斯基(2003年約翰·馮·諾伊曼理論獎英语John von Neumann Theory Prize得主)和 D. Yudin凸集最優化橢球法的一般化。

理論上,「橢球法」在最惡劣的情況下所需要的計算量要比「單形法」增長的緩慢,有希望用之解決超大型線性規劃問題。但在實際應用上,Khachiyan的演算法令人失望:一般來說,單純形演算法比它更有效率。它的重要性在於鼓勵了對内点法的研究。內點演算法是針對單形法的「邊界趨近」觀念而改採「內部逼近」的路線,相對於只沿著可行域的邊沿進行移動的單純形演算法,內點演算法能夠在可行域內移動。

1984年,貝爾實驗室印度裔數學家納倫德拉·卡爾瑪卡爾英语Narendra Karmarkar提出了投影尺度法(又名卡爾瑪卡爾演算法英语Karmarkar's algorithm)。這是第一個在理論上和實際上都表現良好的算法:它的最壞情況僅為多項式時間,且在實際問題中它比單純形演算法有顯著的效率提升。自此之後,很多內點演算法被提出來並進行分析。一個常見的內點演算法為Mehrotra predictor-corrector method。儘管在理論上對它所知甚少,在實際應用中它卻表現出色。

單形法沿著邊界由一個頂點移動到「相鄰」的頂點,內點演算法每一步的移動考量較周詳,「跨過可行解集合的內部」去逼近最佳解。當今的觀點是:對於線性規劃的日常應用問題而言,如果演算法的實現良好,基於單純形法和內點法的演算法之間的效率沒有太大差別,只有在超大型線性規劃中,頂點幾成天文數字,內點法有機會領先單形法。

線性規劃的求解程式在各種各樣的工業最優化問題裡被廣泛使用,例如運輸網路的流量的最優化問題,其中很多都可以不太困難地被轉換成線性規劃問題。

线性规划理论中存在几个尚未解决的问题,这些开放问题的答案将会是数学运算中的根本突破,并且很可能是解决大规模线性规划问题的主要进展。

  • LP存在强多项式时间算法吗?
  • LP存在多项式时间算法以得到一个严格互补解吗?
  • LP在实数(单位成本)模型下存在多项式时间算法吗?

这些问题已经由斯蒂芬·斯梅尔二十一世纪十八个尚未解决的最伟大的问题中应用。用斯梅尔的话来说,“第三个问题是线性规划理论中最主要的尚未解决的问题”。然而,对于线性规划问题存在弱多项式时间算法,比如椭球法和内点法,尚未发现限制在约束条件个数和变量个数的强多项式时间算法,此算法的发展将会带来理论上重大意义,或者是解决大规模线性规划上的实际收益。

儘管赫爾希博士猜想近來被證明是錯誤的,但是它依舊留下下面的開放問題:

  • Are there pivot rules which lead to polynomial-time Simplex variants?
  • Do all polytopal graphs have polynomially-bounded diameter?

整數規劃

要求所有的未知量都為整數的線性規劃問題叫做整數規劃(integer programming, IP)或整數線性規劃(integer linear programming, ILP)問題。相對於即使在最壞情況下也能有效率地解出的線性規劃問題,整數規劃問題的最壞情況是不確定的,在某些實際情況中(有約束變量的那些)為NP困難問題

0-1整數規劃是整數規劃的特殊情況,所有的變量都要是0或1(而非任意整數)。這類問題亦被分類為NP困難問題

只要求當中某幾個未知數為整數的線性規劃問題叫做混合整數規劃(mixed integer programming, MIP)問題。這類問題通常亦被分類為NP困難問題

存在著幾類IP和MIP的子問題,它們可以被有效率地解出,最值得注意的一類是具有完全單位模約束矩陣,和約束條件的右邊全為整數的一類。

一個解決大型整數線性規劃問題的先進演算法為delayed column generation

参见

参考

外部連結

求解軟件包