「線形計画」の版間の差分

提供: ORWiki
ナビゲーションに移動 検索に移動
 
(3人の利用者による、間の7版が非表示)
1行目: 1行目:
'''【せんけいけいかく linear programming】'''
+
'''【せんけいけいかく (linear programming)】'''
  
 +
== 概要 ==
 +
最適化問題(数理計画問題)
 +
 +
 +
<table align="center">
 +
<tr>
 +
<td><math>\mbox{max.} \, </math></td>
 +
<td><math>f(x) \ ( \,</math>あるいは, <math>\min. \ f(x)) \,</math></td>
 +
</tr>
 +
<tr>
 +
<td><math>\mbox{s.t.} \, </math></td>
 +
<td><math>x = (x_1,x_2,\ldots,x_n) \in F,
 +
\,</math></td>
 +
</tr>
 +
</table>
 +
 +
 +
において, 目的関数 <math>f \,</math> が線形であり, かつ, 実行可能集合 <math>F \,</math> が線形等式と線形不等式を用いて表現されている問題.この問題への定式化, および, 解法を含めて線形計画と呼ぶ.
 +
 +
== 詳説 ==
 
 [[線形計画]](linear programming)(線形計画法, 線形計画問題)は, 複数の等式あるいは不等式で与えられる線形制約のもとで, [[目的関数]](objective function)と呼ばれる線形関数を最大化(または最小化)する問題である.  線形計画問題は通常以下の形式により表現される.
 
 [[線形計画]](linear programming)(線形計画法, 線形計画問題)は, 複数の等式あるいは不等式で与えられる線形制約のもとで, [[目的関数]](objective function)と呼ばれる線形関数を最大化(または最小化)する問題である.  線形計画問題は通常以下の形式により表現される.
  
 
+
<table align="center">
:<math>\mbox{(P)} \quad
+
<tr>
 +
<td>
 +
<math>\mbox{(P)} \quad
 
\begin{array}{lll}
 
\begin{array}{lll}
 
           & \mbox{max.} & {\displaystyle \sum_{j=1}^{n}c_j x_j} \\
 
           & \mbox{max.} & {\displaystyle \sum_{j=1}^{n}c_j x_j} \\
11行目: 33行目:
 
       x_1,x_2,\ldots ,x_n \geq 0.
 
       x_1,x_2,\ldots ,x_n \geq 0.
 
\end{array}</math>
 
\end{array}</math>
 +
</td>
 +
</tr>
 +
</table>
  
  
ここで, <math>c_j (j=1,2,\ldots,n), \, b_i (i=1,2,\ldots,m),\,a_{ij} (i=1,2,\ldots,m, \ j=1,2,\ldots,n)</math>は実数, '''''x''''' <math>=(x_1,x_2,\ldots, x_m)</math>は<math>n\,</math>個の変数からなる<math>n\,</math>次元ベクトルである.  すべての線形計画問題は, 簡単な変換により, この形式に帰着される.  問題(P)の制約条件を満たす'''''x'''''を[[実行可能解]](feasible solution)と呼び, 実行可能解のなかで目的関数を最大にするものを[[最適解]]と呼ぶ.
+
ここで, <math>c_j (j=1,2,\ldots,n), \, b_i (i=1,2,\ldots,m),\,a_{ij} (i=1,2,\ldots,m, \ j=1,2,\ldots,n)</math>は実数, <math>\boldsymbol{x}=(x_1,x_2,\ldots, x_m)</math>は<math>n\,</math>個の変数からなる<math>n\,</math>次元ベクトルである.  すべての線形計画問題は, 簡単な変換により, この形式に帰着される.  問題(P)の制約条件を満たす<math>\boldsymbol{x}</math>を[[実行可能解]](feasible solution)と呼び, 実行可能解のなかで目的関数を最大にするものを[[最適解]]と呼ぶ.
  
 
 実社会における多くの問題は線形計画問題として定式化できる.  線形計画の応用は, 初期の頃は, 軍事, 経済学やゲーム理論が中心であったが, 次第にその重点が産業の分野へと移行された [2, 3].  1947年にダンツィク(G. B. Dantzig)[2] によって提案された[[単体法]](simplex method)は, コンピュータの急速な発展と大規模な線形方程式を処理する技術の向上とあいまって, 線形計画問題に対する極めて実用的な解法となっている.  しかし, 単体法は Klee-Minty [4] により, 問題の入力サイズ(変数の個数と制約の本数)に関する多項式時間解法ではないことが指摘された.  カチヤン(L. G. Khachian)は[[楕円体法]]を提案し, 最初の多項式時間解法を与えた.  [[カーマーカー法]]およびその後開発された[[内点法]]は, 超大規模な線形計画問題に対し, 理論および実用の両面において, 単体法より優れた解法として認められている.
 
 実社会における多くの問題は線形計画問題として定式化できる.  線形計画の応用は, 初期の頃は, 軍事, 経済学やゲーム理論が中心であったが, 次第にその重点が産業の分野へと移行された [2, 3].  1947年にダンツィク(G. B. Dantzig)[2] によって提案された[[単体法]](simplex method)は, コンピュータの急速な発展と大規模な線形方程式を処理する技術の向上とあいまって, 線形計画問題に対する極めて実用的な解法となっている.  しかし, 単体法は Klee-Minty [4] により, 問題の入力サイズ(変数の個数と制約の本数)に関する多項式時間解法ではないことが指摘された.  カチヤン(L. G. Khachian)は[[楕円体法]]を提案し, 最初の多項式時間解法を与えた.  [[カーマーカー法]]およびその後開発された[[内点法]]は, 超大規模な線形計画問題に対し, 理論および実用の両面において, 単体法より優れた解法として認められている.
  
 以下では, 線形計画の理論で重要な役割を果たす[[双対問題]](dual problem), [[双対定理]](duality theorem)および[[相補性定理]](complementarity slackness theorem)について述べる.
+
 以下では, 線形計画の理論で重要な役割を果たす[[双対問題 (線形計画の)|双対問題]](dual problem), [[双対定理]](duality theorem)および[[相補性定理]](complementarity slackness theorem)について述べる.
  
 
 各々の線形計画問題(P)には, 双対問題と呼ばれる以下の線形計画問題(D)を対応させることができる:
 
 各々の線形計画問題(P)には, 双対問題と呼ばれる以下の線形計画問題(D)を対応させることができる:
  
  
:<math>\mbox{(D)} \quad
+
<table align="center">
 +
<tr>
 +
<td>
 +
<math>\mbox{(D)} \quad
 
\begin{array}{lll}
 
\begin{array}{lll}
 
           & \mbox{min.} & {\displaystyle \sum_{i=1}^{m}b_i y_i} \\
 
           & \mbox{min.} & {\displaystyle \sum_{i=1}^{m}b_i y_i} \\
28行目: 56行目:
 
                           \geq c_j \;  (j=1,2,\ldots,n),  
 
                           \geq c_j \;  (j=1,2,\ldots,n),  
 
           y_1,y_2,\ldots,y_m \geq 0.
 
           y_1,y_2,\ldots,y_m \geq 0.
\end{array}</math>
+
\end{array}</math></td>
 +
</tr>
 +
</table>
  
  
 
ここで, <math>c_j \, (j=1,2,\ldots,n), \, b_i \, (i=1,2,\ldots,m),\, a_{ij} \, (i=1,2,\ldots,m, \ j=1,2,\ldots,n)</math>は元の問題(P)で与えられたデータと同一である.  元の問題(P)は主問題, 問題(D)は主問題(P)に対する双対問題と呼ばれる.  双対問題の双対問題が主問題になることは, 問題(D)を問題(P)の形式に書き直すことにより容易に示せる.  主問題と双対問題の実行可能解について, 次のようなことが言える.
 
ここで, <math>c_j \, (j=1,2,\ldots,n), \, b_i \, (i=1,2,\ldots,m),\, a_{ij} \, (i=1,2,\ldots,m, \ j=1,2,\ldots,n)</math>は元の問題(P)で与えられたデータと同一である.  元の問題(P)は主問題, 問題(D)は主問題(P)に対する双対問題と呼ばれる.  双対問題の双対問題が主問題になることは, 問題(D)を問題(P)の形式に書き直すことにより容易に示せる.  主問題と双対問題の実行可能解について, 次のようなことが言える.
  
*主問題の実行可能解 '''''x''''' と双対問題の実行可能解 '''''y''''' に対し, <math>\sum_{j=1}^{n}c_jx_j\leq \sum_{i=1}^{m}b_i y_i</math>が成り立つ.
+
:*主問題の実行可能解<math>\boldsymbol{x}</math>と双対問題の実行可能解<math>\boldsymbol{y}</math>に対し, <math>\textstyle \sum_{j=1}^{n}c_jx_j\leq \sum_{i=1}^{m}b_i y_i</math>が成り立つ.
*主問題の実行可能解 '''''x''''' と双対問題の実行可能解 '''''y''''' に対し, それらの目的関数値が一致するならば, '''''x''''' '''''y''''' はそれぞれの問題の最適解である.
+
 
 +
:*主問題の実行可能解<math>\boldsymbol{x}</math>と双対問題の実行可能解<math>\boldsymbol{y}</math>に対し, それらの目的関数値が一致するならば, <math>\boldsymbol{x}</math><math>\boldsymbol{y}</math>はそれぞれの問題の最適解である.
  
 
さらに, 以下の双対定理と相補性定理が示すように, 主問題と双対問題のいずれか一方の最適解は, 他方の問題の最適解に関する情報を含む.
 
さらに, 以下の双対定理と相補性定理が示すように, 主問題と双対問題のいずれか一方の最適解は, 他方の問題の最適解に関する情報を含む.
40行目: 71行目:
 
'''双対定理:''' 主問題か双対問題のいずれか一方が最適解をもつならば, 他方もまた最適解をもち, それらの目的関数値は一致する.  また, いずれか一方の問題の実行可能解に対する目的関数値が有界でなければ, 他方の問題は実行可能解をもたない.
 
'''双対定理:''' 主問題か双対問題のいずれか一方が最適解をもつならば, 他方もまた最適解をもち, それらの目的関数値は一致する.  また, いずれか一方の問題の実行可能解に対する目的関数値が有界でなければ, 他方の問題は実行可能解をもたない.
  
'''相補性定理:''' 主問題の実行可能解 '''''x''''' と双対問題の実行可能解 '''''y''''' がそれぞれの問題の最適解であるための必要十分条件は, 以下の2条件(i),(ii)が成り立つことである.
+
'''相補性定理:''' 主問題の実行可能解<math>\boldsymbol{x}</math>と双対問題の実行可能解<math>\boldsymbol{y}</math>がそれぞれの問題の最適解であるための必要十分条件は, 以下の2条件(i),(ii)が成り立つことである.
 
 
  
:<math>\mbox{(i)} \quad(c_j-\sum_{i=1}^{m}a_{ij}y_i)x_j=0, \ j=1,2,\ldots,n</math>
 
  
:<math>\mbox{(ii)}\quad(\sum_{j=1}^{n}a_{ij}x_j-b_i)y_i =0, \ i=1,2,\ldots,m</math>
+
<table align="center">
 +
<tr>
 +
<td>
 +
<math>\begin{array}{lll}
 +
\mbox{(i)} & (c_j-\sum_{i=1}^{m}a_{ij}y_i)x_j=0,\ j=1,2,\ldots,n,\\
 +
\\
 +
\mbox{(ii)} & (\sum_{j=1}^{n}a_{ij}x_j-b_i)y_i =0,\ i=1,2,\ldots,m.
 +
\end{array}</math></td>
 +
</tr>
 +
</table>
  
  
52行目: 90行目:
  
 
----
 
----
参考文献
+
'''参考文献'''
<BR>[1]  V. Chv&aacute;tal, ''Linear Programming'', W. H. Freeman and Company, 1983, 阪田省二郎, 藤野和建 訳, 『線形計画法』上, 下, 啓学出版.
+
 
<BR>[2]  G. B. Dantzig, ''Linear Programming and Extensions'', Princeton University Press, 1963.
+
[1]  V. Chv&aacute;tal, ''Linear Programming'', W. H. Freeman and Company, 1983, 阪田省二郎, 藤野和建 訳, 『線形計画法』上, 下, 啓学出版.
<BR>[3]  S. I. Gass, ''Linear Programming and Extensions'', MaGram-Hill, 1975.
+
 
<BR>[4]  V. Klee and G. J. Minty "How good is the simplex algorithm?" in ''Inequalities-III'', O. Shisha, eds., Academic Press, 159-175, 1989.
+
[2]  G. B. Dantzig, ''Linear Programming and Extensions'', Princeton University Press, 1963.
 +
 
 +
[3]  S. I. Gass, ''Linear Programming and Extensions'', MaGram-Hill, 1975.
 +
 
 +
[4]  V. Klee and G. J. Minty "How good is the simplex algorithm?" in ''Inequalities-III'', O. Shisha, eds., Academic Press, 159-175, 1989.
 +
 
 +
[[Category:線形計画|せんけいけいかく]]

2008年8月6日 (水) 17:02時点における最新版

【せんけいけいかく (linear programming)】

概要

最適化問題(数理計画問題)


あるいは,


において, 目的関数 が線形であり, かつ, 実行可能集合 が線形等式と線形不等式を用いて表現されている問題.この問題への定式化, および, 解法を含めて線形計画と呼ぶ.

詳説

 線形計画(linear programming)(線形計画法, 線形計画問題)は, 複数の等式あるいは不等式で与えられる線形制約のもとで, 目的関数(objective function)と呼ばれる線形関数を最大化(または最小化)する問題である. 線形計画問題は通常以下の形式により表現される.


ここで, は実数, 個の変数からなる次元ベクトルである. すべての線形計画問題は, 簡単な変換により, この形式に帰着される. 問題(P)の制約条件を満たす実行可能解(feasible solution)と呼び, 実行可能解のなかで目的関数を最大にするものを最適解と呼ぶ.

 実社会における多くの問題は線形計画問題として定式化できる. 線形計画の応用は, 初期の頃は, 軍事, 経済学やゲーム理論が中心であったが, 次第にその重点が産業の分野へと移行された [2, 3]. 1947年にダンツィク(G. B. Dantzig)[2] によって提案された単体法(simplex method)は, コンピュータの急速な発展と大規模な線形方程式を処理する技術の向上とあいまって, 線形計画問題に対する極めて実用的な解法となっている. しかし, 単体法は Klee-Minty [4] により, 問題の入力サイズ(変数の個数と制約の本数)に関する多項式時間解法ではないことが指摘された. カチヤン(L. G. Khachian)は楕円体法を提案し, 最初の多項式時間解法を与えた. カーマーカー法およびその後開発された内点法は, 超大規模な線形計画問題に対し, 理論および実用の両面において, 単体法より優れた解法として認められている.

 以下では, 線形計画の理論で重要な役割を果たす双対問題(dual problem), 双対定理(duality theorem)および相補性定理(complementarity slackness theorem)について述べる.

 各々の線形計画問題(P)には, 双対問題と呼ばれる以下の線形計画問題(D)を対応させることができる:



ここで, は元の問題(P)で与えられたデータと同一である. 元の問題(P)は主問題, 問題(D)は主問題(P)に対する双対問題と呼ばれる. 双対問題の双対問題が主問題になることは, 問題(D)を問題(P)の形式に書き直すことにより容易に示せる. 主問題と双対問題の実行可能解について, 次のようなことが言える.

  • 主問題の実行可能解と双対問題の実行可能解に対し, が成り立つ.
  • 主問題の実行可能解と双対問題の実行可能解に対し, それらの目的関数値が一致するならば, はそれぞれの問題の最適解である.

さらに, 以下の双対定理と相補性定理が示すように, 主問題と双対問題のいずれか一方の最適解は, 他方の問題の最適解に関する情報を含む.

双対定理: 主問題か双対問題のいずれか一方が最適解をもつならば, 他方もまた最適解をもち, それらの目的関数値は一致する. また, いずれか一方の問題の実行可能解に対する目的関数値が有界でなければ, 他方の問題は実行可能解をもたない.

相補性定理: 主問題の実行可能解と双対問題の実行可能解がそれぞれの問題の最適解であるための必要十分条件は, 以下の2条件(i),(ii)が成り立つことである.



これらの最適性条件を次のように言い換えることができる. いずれか一方の問題の番目の制約式が不等号で成り立つならば, 他方の問題において対応する変数の値はゼロになる. また, いずれか一方の問題の番目の変数の値が正ならば, 他方の問題においてそれに対応する制約式は等号で成り立つ.



参考文献

[1] V. Chvátal, Linear Programming, W. H. Freeman and Company, 1983, 阪田省二郎, 藤野和建 訳, 『線形計画法』上, 下, 啓学出版.

[2] G. B. Dantzig, Linear Programming and Extensions, Princeton University Press, 1963.

[3] S. I. Gass, Linear Programming and Extensions, MaGram-Hill, 1975.

[4] V. Klee and G. J. Minty "How good is the simplex algorithm?" in Inequalities-III, O. Shisha, eds., Academic Press, 159-175, 1989.