<?xml version="1.0"?>
<feed xmlns="http://www.w3.org/2005/Atom" xml:lang="ja">
	<id>https://orsj-ml.org/orwiki/wiki/api.php?action=feedcontributions&amp;feedformat=atom&amp;user=Tamura</id>
	<title>ORWiki - 利用者の投稿記録 [ja]</title>
	<link rel="self" type="application/atom+xml" href="https://orsj-ml.org/orwiki/wiki/api.php?action=feedcontributions&amp;feedformat=atom&amp;user=Tamura"/>
	<link rel="alternate" type="text/html" href="https://orsj-ml.org/orwiki/wiki/index.php?title=%E7%89%B9%E5%88%A5:%E6%8A%95%E7%A8%BF%E8%A8%98%E9%8C%B2/Tamura"/>
	<updated>2026-04-08T09:58:35Z</updated>
	<subtitle>利用者の投稿記録</subtitle>
	<generator>MediaWiki 1.35.3</generator>
	<entry>
		<id>https://orsj-ml.org/orwiki/wiki/index.php?title=%E5%86%85%E7%82%B9%E6%B3%95&amp;diff=9995</id>
		<title>内点法</title>
		<link rel="alternate" type="text/html" href="https://orsj-ml.org/orwiki/wiki/index.php?title=%E5%86%85%E7%82%B9%E6%B3%95&amp;diff=9995"/>
		<updated>2008-08-22T07:24:58Z</updated>

		<summary type="html">&lt;p&gt;Tamura: リンク法の修正と読書案内という見出しの追加&lt;/p&gt;
&lt;hr /&gt;
&lt;div&gt;'''【ないてんほう (interior point method)】'''&lt;br /&gt;
&lt;br /&gt;
=== 概要 ===&lt;br /&gt;
線形計画問題に対するカーマーカー法によって触発され,発展した制約付き最適化問題の反復解法の総称. 実行可能領域の内部を通って最適解へ近づくのでこう呼ばれる.特に大規模問題に有効とされ, [[凸計画問題]], [[半正定値計画問題]]等へも拡張されている. 主問題または双対問題のどちらか一方の実行可能領域に点列を生成する内点法を主内点法, 主問題および双対問題の双方の実行可能領域に点列を生成する内点法を主双対内点法と呼ぶ.&lt;br /&gt;
&lt;br /&gt;
=== 詳説 ===&lt;br /&gt;
　[[内点法]](interior point method) は[[最適化問題]]の制約領域の境界上ではなく内部に, 最適解に収束する点列を生成する逐次反復解法である. 特に[[線形計画|線形計画問題]]に対しては, 最適値に十分に近い近似解が得られれば, 変数次元の連立方程式を1回解く手間で最適解が得られることから, 内点法で最適解を得るまでに必要な計算の手間(算術演算回数)を算出することが可能であり, 多くの理論的な成果が報告されている. また実用面では, 特に大規模な問題に対して単体法よりも計算効率が 優れていることが確認されており, この解法を取り入れた商用コードも多く存在する. 線形計画問題以外への応用も積極的に行なわれており, 凸2次計画問題, [[相補性問題]], あるいは半正定値計画問題, [[2次錐計画|2次錐計画問題]]等の凸計画問題に対しても有効な解法であることが示されている.&lt;br /&gt;
&lt;br /&gt;
　1979年に提案された楕円体法は, 線形計画問題に対する初めての[[多項式時間解法]]であるが, 実用的な解法としては単体法に及ばなかった. これに対して1984年にATTベル研究所のカーマーカー(N. Karmarkar) によって提案された初めての内点法(interior point)である[[カーマーカー法]](Karmarkar's algorithm) [1] は, 多項式時間解法である上に, カーマーカー自身の計算機実験では 単体法を遥かに上回る結果であると報告された. この論文および報告を機として, 多くの最適化分野の研究者がこの解法に興味をもち, この結果多様なバリエーションが生み出された.&lt;br /&gt;
&lt;br /&gt;
　カーマーカー法は, 特殊な線形計画問題&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;\begin{array}{llllll}&lt;br /&gt;
\mbox{max.} &amp;amp; c^{\top}x  &amp;amp; \mbox{s. t.}  &amp;amp; Ax = 0,&lt;br /&gt;
&amp;amp; e^{\top}x = 1,  &amp;amp; x \geq 0, &lt;br /&gt;
\end{array}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
(&amp;lt;math&amp;gt;A\,&amp;lt;/math&amp;gt;は&amp;lt;math&amp;gt;m \times n\,&amp;lt;/math&amp;gt;行列, &amp;lt;math&amp;gt;c \in \mathbf{R}^n&amp;lt;/math&amp;gt;, &amp;lt;math&amp;gt;e \in \mathbf{R}^n&amp;lt;/math&amp;gt;は要素がすべて&amp;lt;math&amp;gt;1\,&amp;lt;/math&amp;gt;のベクトル)を対象とし, 対数障壁関数を用いた[[ポテンシャル関数 (内点法の)|ポテンシャル関数]](potential function) の値で最適解からの誤差を測り, [[射影変換]] (projective transformation) を用いて探索方向を決定するという特徴をもつ. 対数障壁関数を用いるという点で非線形計画法の一解法とも考えられるが, 着想, 解析共に従来の常識とは異なる斬新な解法であった. より一般的で記述が簡潔なアルゴリズムとして提案された解法が, [[アフィン変換法]](affine scaling method) である. その後1988年に, ロシアの数学者ディキン(I. Dikin) が20年以上前の1967年に同じ解法を提案していたことが明かになり, 話題となった. アフィン変換法は実用的であり, 大域的収束性が示されているが, 解法の多項式時間性については未だ不明である.&lt;br /&gt;
&lt;br /&gt;
　他方, 従来のSUMT法, ホモトピー法, [[ニュートン法]]といった枠組を用いて内点法を構築する試みも行なわれた. 中でも大きな役割を果たしたのは, ゾンネべンド(Gy. Sonnevend)が提唱した[[中心パス]](path of centers) であり, この概念を利用してレネガー(J. Renegar)はニュートン法を用いた初めての 多項式時間解法を提案した. さらに, 主問題と双対問題の双方を1つの問題とみなした場合の解析的 中心に関する研究も行なわれ, [[主双対内点法]](primal-dual interior point method) の提案へと結び付いた. 主双対内点法を改良した[[予測子修正子内点法]](predictor-corrector interior point method) は, 現在もっとも普及している内点法の1つである. これらの解法については [5] 等を参照されたい.&lt;br /&gt;
&lt;br /&gt;
　内点法における大きな課題の1つは, 条件のよい許容領域内部の初期点をいかに用意するかという問題である. この課題の解決策として [[非許容初期点内点法]] (infeasible interior point method) や, 理論的にはさらに優れた [[同次自己双対内点法]] (homogeneous self-dual interior point method) が提案されており, 実装の際にはどちらかが用いられることが多い.&lt;br /&gt;
&lt;br /&gt;
　以上の解法の中で, アフィン変換法以外の解法は多項式時間解法であることが示されている. ネステロフとネミロフスキーは, [[自己整合障壁関数]]の概念を導入した上で, 内点法という解法のクラスの中で, 多項式時間解法である解法からなる, 1つの大きなクラスを示した [2]. この研究は, 1990年代から活発に行なわれている, [[半正定値計画|半正定値計画問題]]や[[2次錐計画|2次錐計画問題]]等に内点法を適用する研究の重要な基礎となっている.&lt;br /&gt;
&lt;br /&gt;
=== 読書案内 ===&lt;br /&gt;
&lt;br /&gt;
　近年内点法が記述された教科書も多く出版されている. [12] は内点法を総合的に扱った和書である. 線形計画法のみならず, 半正定値計画問題, 非線形計画問題に対する内点法についも述べている. その他の洋書の多くは主双対内点法と非許容初期点内点法を基本とした内容であるが, カーマーカー法については [3], [9], アフィン変換法については [4], [6], 同次自己双対内点法については, [6], [10] 等にそれぞれ詳しく述べられている. [8], [9], [10] はそれぞれ線形計画法の教科書であり, 図表が多く入門的な内容であるが, [8] は単体法, [9] は内点法を包括する理論, [10] は実装 にそれぞれ重きを置いている点が特徴である. また, 内点法に関する論文の多くは, インターネットを通じて入手可能である [11].&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
----&lt;br /&gt;
'''参考文献'''&lt;br /&gt;
&lt;br /&gt;
[1] K. Karmarkar, &amp;quot;A New Polynomial-Time Algorithm for Linear Programming,&amp;quot; ''Combinatorica'', '''4''' (1984), 373-395.&lt;br /&gt;
&lt;br /&gt;
[2] Y. Nesterov and A. S. Nemirovskii, ''Interior Point Polynomial Algorithms in Convex Programming'', SIAM, 1994.&lt;br /&gt;
&lt;br /&gt;
[3] S.-C. Fang and S. Puthenpra, ''Linear Optimization and Extentions'' Prentice Hall, 1993.&lt;br /&gt;
&lt;br /&gt;
[4] R. Saigal, ''Linear Programming: A Modern Integrated Analysis'' Kluwer Academic Publishers, 1995.&lt;br /&gt;
&lt;br /&gt;
[5] J. Wright, ''Primal-Dual Interior-Point Methods'', SIAM, 1996.&lt;br /&gt;
&lt;br /&gt;
[6] T. Terlaky, eds., ''Interior Point Methods of Mathematical Programming'', Kluwer Academic Publishers, 1996.&lt;br /&gt;
&lt;br /&gt;
[7] Y. Ye, ''Interior Point Algorithms: Theory and Analysis'', Jhon Wiley &amp;amp; Sons, 1997&lt;br /&gt;
&lt;br /&gt;
[8] D. Bertsimas and J. N. Tsisiklis, ''Introduction to Linear Optimization'', Athena Scientific, 1997.&lt;br /&gt;
&lt;br /&gt;
[9] C. Roos, T. Terlaky and J.-Ph. Vial ''Theory and Algorithns for Linear Optimization'', John Wiley &amp;amp; Sons, 1997.&lt;br /&gt;
&lt;br /&gt;
[10] R. J. Vanderbei, ''Linear Programming: Foundations and Extensions'', Kluwer Academic Publishers, 1998.&lt;br /&gt;
&lt;br /&gt;
[11] http://www-unix.mcs.anl.gov/otc/InteriorPoint&lt;br /&gt;
&lt;br /&gt;
[12] 小島政和, 土谷隆, 水野眞治, 矢部博, 『内点法』, 朝倉書店, 2001.&lt;br /&gt;
&lt;br /&gt;
[[Category:線形計画|ないてんほう]]&lt;/div&gt;</summary>
		<author><name>Tamura</name></author>
	</entry>
	<entry>
		<id>https://orsj-ml.org/orwiki/wiki/index.php?title=%E5%86%85%E7%82%B9%E6%B3%95&amp;diff=9994</id>
		<title>内点法</title>
		<link rel="alternate" type="text/html" href="https://orsj-ml.org/orwiki/wiki/index.php?title=%E5%86%85%E7%82%B9%E6%B3%95&amp;diff=9994"/>
		<updated>2008-08-22T07:22:31Z</updated>

		<summary type="html">&lt;p&gt;Tamura: &lt;/p&gt;
&lt;hr /&gt;
&lt;div&gt;'''【ないてんほう (interior point method)】'''&lt;br /&gt;
&lt;br /&gt;
=== 概要 ===&lt;br /&gt;
線形計画問題に対するカーマーカー法によって触発され,発展した制約付き最適化問題の反復解法の総称. 実行可能領域の内部を通って最適解へ近づくのでこう呼ばれる.特に大規模問題に有効とされ, [[凸計画問題]], [[半正定値計画問題]]等へも拡張されている. 主問題または双対問題のどちらか一方の実行可能領域に点列を生成する内点法を主内点法, 主問題および双対問題の双方の実行可能領域に点列を生成する内点法を主双対内点法と呼ぶ.&lt;br /&gt;
&lt;br /&gt;
=== 詳説 ===&lt;br /&gt;
　[[内点法]](interior point method) は[[最適化問題]]の制約領域の境界上ではなく内部に, 最適解に収束する点列を生成する逐次反復解法である. 特に[[線形計画|線形計画問題]]に対しては, 最適値に十分に近い近似解が得られれば, 変数次元の連立方程式を1回解く手間で最適解が得られることから, 内点法で最適解を得るまでに必要な計算の手間(算術演算回数)を算出することが可能であり, 多くの理論的な成果が報告されている. また実用面では, 特に大規模な問題に対して単体法よりも計算効率が 優れていることが確認されており, この解法を取り入れた商用コードも多く存在する. 線形計画問題以外への応用も積極的に行なわれており, 凸2次計画問題, [[相補性問題]], あるいは半正定値計画問題, [[2次錐計画|2次錐計画問題]]等の凸計画問題に対しても有効な解法であることが示されている.&lt;br /&gt;
&lt;br /&gt;
　1979年に提案された楕円体法は, 線形計画問題に対する初めての[[多項式時間解法]]であるが, 実用的な解法としては単体法に及ばなかった. これに対して1984年にATTベル研究所のカーマーカー(N. Karmarkar) によって提案された初めての内点法(interior point)である[[カーマーカー法]](Karmarkar's algorithm) [1] は, 多項式時間解法である上に, カーマーカー自身の計算機実験では 単体法を遥かに上回る結果であると報告された. この論文および報告を機として, 多くの最適化分野の研究者がこの解法に興味をもち, この結果多様なバリエーションが生み出された.&lt;br /&gt;
&lt;br /&gt;
　カーマーカー法は, 特殊な線形計画問題&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;\begin{array}{llllll}&lt;br /&gt;
\mbox{max.} &amp;amp; c^{\top}x  &amp;amp; \mbox{s. t.}  &amp;amp; Ax = 0,&lt;br /&gt;
&amp;amp; e^{\top}x = 1,  &amp;amp; x \geq 0, &lt;br /&gt;
\end{array}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
(&amp;lt;math&amp;gt;A\,&amp;lt;/math&amp;gt;は&amp;lt;math&amp;gt;m \times n\,&amp;lt;/math&amp;gt;行列, &amp;lt;math&amp;gt;c \in \mathbf{R}^n&amp;lt;/math&amp;gt;, &amp;lt;math&amp;gt;e \in \mathbf{R}^n&amp;lt;/math&amp;gt;は要素がすべて&amp;lt;math&amp;gt;1\,&amp;lt;/math&amp;gt;のベクトル)を対象とし, 対数障壁関数を用いた[[ポテンシャル関数 (内点法の)|ポテンシャル関数]](potential function) の値で最適解からの誤差を測り, [[射影変換]] (projective transformation) を用いて探索方向を決定するという特徴をもつ. 対数障壁関数を用いるという点で非線形計画法の一解法とも考えられるが, 着想, 解析共に従来の常識とは異なる斬新な解法であった. より一般的で記述が簡潔なアルゴリズムとして提案された解法が, [[アフィン変換法]](affine scaling method) である. その後1988年に, ロシアの数学者ディキン(I. Dikin) が20年以上前の1967年に同じ解法を提案していたことが明かになり, 話題となった. アフィン変換法は実用的であり, 大域的収束性が示されているが, 解法の多項式時間性については未だ不明である.&lt;br /&gt;
&lt;br /&gt;
　他方, 従来のSUMT法, ホモトピー法, [[ニュートン法]]といった枠組を用いて内点法を構築する試みも行なわれた. 中でも大きな役割を果たしたのは, ゾンネべンド(Gy. Sonnevend)が提唱した[[中心パス]](path of centers) であり, この概念を利用してレネガー(J. Renegar)はニュートン法を用いた初めての 多項式時間解法を提案した. さらに, 主問題と双対問題の双方を1つの問題とみなした場合の解析的 中心に関する研究も行なわれ, [[主双対内点法]](primal-dual interior point method) の提案へと結び付いた. 主双対内点法を改良した[[予測子修正子内点法]](predictor-corrector interior point method) は, 現在もっとも普及している内点法の1つである. これらの解法については [5] 等を参照されたい.&lt;br /&gt;
&lt;br /&gt;
　内点法における大きな課題の1つは, 条件のよい許容領域内部の初期点をいかに用意するかという問題である. この課題の解決策として [[非許容初期点内点法]] (infeasible interior point method) や, 理論的にはさらに優れた [[同次自己双対内点法]] (homogeneous self-dual interior point method) が提案されており, 実装の際にはどちらかが用いられることが多い.&lt;br /&gt;
&lt;br /&gt;
　以上の解法の中で, アフィン変換法以外の解法は多項式時間解法であることが示されている. ネステロフとネミロフスキーは, [[自己整合障壁関数]]の概念を導入した上で, 内点法という解法のクラスの中で, 多項式時間解法である解法からなる, 1つの大きなクラスを示した [2]. この研究は, 1990年代から活発に行なわれている, [[半正定値計画]]問題や[[2次錐計画]]問題等に内点法を適用する研究の重要な基礎となっている.&lt;br /&gt;
&lt;br /&gt;
　近年内点法が記述された教科書も多く出版されている. [12] は内点法を総合的に扱った和書である. 線形計画法のみならず, 半正定値計画問題, 非線形計画問題に対する内点法についも述べている. その他の洋書の多くは主双対内点法と非許容初期点内点法を基本とした内容であるが, カーマーカー法については [3], [9], アフィン変換法については [4], [6], 同次自己双対内点法については, [6], [10] 等にそれぞれ詳しく述べられている. [8], [9], [10] はそれぞれ線形計画法の教科書であり, 図表が多く入門的な内容であるが, [8] は単体法, [9] は内点法を包括する理論, [10] は実装 にそれぞれ重きを置いている点が特徴である. また, 内点法に関する論文の多くは, インターネットを通じて入手可能である [11].&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
----&lt;br /&gt;
'''参考文献'''&lt;br /&gt;
&lt;br /&gt;
[1] K. Karmarkar, &amp;quot;A New Polynomial-Time Algorithm for Linear Programming,&amp;quot; ''Combinatorica'', '''4''' (1984), 373-395.&lt;br /&gt;
&lt;br /&gt;
[2] Y. Nesterov and A. S. Nemirovskii, ''Interior Point Polynomial Algorithms in Convex Programming'', SIAM, 1994.&lt;br /&gt;
&lt;br /&gt;
[3] S.-C. Fang and S. Puthenpra, ''Linear Optimization and Extentions'' Prentice Hall, 1993.&lt;br /&gt;
&lt;br /&gt;
[4] R. Saigal, ''Linear Programming: A Modern Integrated Analysis'' Kluwer Academic Publishers, 1995.&lt;br /&gt;
&lt;br /&gt;
[5] J. Wright, ''Primal-Dual Interior-Point Methods'', SIAM, 1996.&lt;br /&gt;
&lt;br /&gt;
[6] T. Terlaky, eds., ''Interior Point Methods of Mathematical Programming'', Kluwer Academic Publishers, 1996.&lt;br /&gt;
&lt;br /&gt;
[7] Y. Ye, ''Interior Point Algorithms: Theory and Analysis'', Jhon Wiley &amp;amp; Sons, 1997&lt;br /&gt;
&lt;br /&gt;
[8] D. Bertsimas and J. N. Tsisiklis, ''Introduction to Linear Optimization'', Athena Scientific, 1997.&lt;br /&gt;
&lt;br /&gt;
[9] C. Roos, T. Terlaky and J.-Ph. Vial ''Theory and Algorithns for Linear Optimization'', John Wiley &amp;amp; Sons, 1997.&lt;br /&gt;
&lt;br /&gt;
[10] R. J. Vanderbei, ''Linear Programming: Foundations and Extensions'', Kluwer Academic Publishers, 1998.&lt;br /&gt;
&lt;br /&gt;
[11] http://www-unix.mcs.anl.gov/otc/InteriorPoint&lt;br /&gt;
&lt;br /&gt;
[12] 小島政和, 土谷隆, 水野眞治, 矢部博, 『内点法』, 朝倉書店, 2001.&lt;br /&gt;
&lt;br /&gt;
[[Category:線形計画|ないてんほう]]&lt;/div&gt;</summary>
		<author><name>Tamura</name></author>
	</entry>
	<entry>
		<id>https://orsj-ml.org/orwiki/wiki/index.php?title=%E5%86%85%E7%82%B9%E6%B3%95&amp;diff=9993</id>
		<title>内点法</title>
		<link rel="alternate" type="text/html" href="https://orsj-ml.org/orwiki/wiki/index.php?title=%E5%86%85%E7%82%B9%E6%B3%95&amp;diff=9993"/>
		<updated>2008-08-22T07:20:22Z</updated>

		<summary type="html">&lt;p&gt;Tamura: &lt;/p&gt;
&lt;hr /&gt;
&lt;div&gt;'''【ないてんほう (interior point method)】'''&lt;br /&gt;
&lt;br /&gt;
=== 概要 ===&lt;br /&gt;
線形計画問題に対するカーマーカー法によって触発され,発展した制約付き最適化問題の反復解法の総称. 実行可能領域の内部を通って最適解へ近づくのでこう呼ばれる.特に大規模問題に有効とされ, [[凸計画問題]], [[半正定値計画問題]]等へも拡張されている. 主問題または双対問題のどちらか一方の実行可能領域に点列を生成する内点法を主内点法, 主問題および双対問題の双方の実行可能領域に点列を生成する内点法を主双対内点法と呼ぶ.&lt;br /&gt;
&lt;br /&gt;
=== 詳説 ===&lt;br /&gt;
　[[内点法]](interior point method) は[[最適化問題]]の制約領域の境界上ではなく内部に, 最適解に収束する点列を生成する逐次反復解法である. 特に[[線形計画|線形計画問題]]に対しては, 最適値に十分に近い近似解が得られれば, 変数次元の連立方程式を1回解く手間で最適解が得られることから, 内点法で最適解を得るまでに必要な計算の手間(算術演算回数)を算出することが可能であり, 多くの理論的な成果が報告されている. また実用面では, 特に大規模な問題に対して単体法よりも計算効率が 優れていることが確認されており, この解法を取り入れた商用コードも多く存在する. 線形計画問題以外への応用も積極的に行なわれており, 凸2次計画問題, [[相補性問題]], あるいは半正定値計画問題, [[2次錐計画|2次錐計画問題]]等の凸計画問題に対しても有効な解法であることが示されている.&lt;br /&gt;
&lt;br /&gt;
　1979年に提案された楕円体法は, 線形計画問題に対する初めての[[多項式時間解法]]であるが, 実用的な解法としては単体法に及ばなかった. これに対して1984年にATTベル研究所のカーマーカー(N. Karmarkar) によって提案された初めての内点法(interior point)である[[カーマーカー法]](Karmarkar's algorithm) [1] は, 多項式時間解法である上に, カーマーカー自身の計算機実験では 単体法を遥かに上回る結果であると報告された. この論文および報告を機として, 多くの最適化分野の研究者がこの解法に興味をもち, この結果多様なバリエーションが生み出された.&lt;br /&gt;
&lt;br /&gt;
　カーマーカー法は, 特殊な線形計画問題&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;\begin{array}{llllll}&lt;br /&gt;
\mbox{max.} &amp;amp; c^{\top}x  &amp;amp; \mbox{s. t.}  &amp;amp; Ax = 0,&lt;br /&gt;
&amp;amp; e^{\top}x = 1,  &amp;amp; x \geq 0, &lt;br /&gt;
\end{array}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
(&amp;lt;math&amp;gt;A\,&amp;lt;/math&amp;gt;は&amp;lt;math&amp;gt;m \times n\,&amp;lt;/math&amp;gt;行列, &amp;lt;math&amp;gt;c \in \mathbf{R}^n&amp;lt;/math&amp;gt;, &amp;lt;math&amp;gt;e \in \mathbf{R}^n&amp;lt;/math&amp;gt;は要素がすべて&amp;lt;math&amp;gt;1\,&amp;lt;/math&amp;gt;のベクトル)を対象とし, 対数障壁関数を用いた[[ポテンシャル関数 (内点法の)|ポテンシャル関数]](potential function) の値で最適解からの誤差を測り, [[射影変換]] (projective transformation) を用いて探索方向を決定するという特徴をもつ. 対数障壁関数を用いるという点で非線形計画法の一解法とも考えられるが, 着想, 解析共に従来の常識とは異なる斬新な解法であった. より一般的で記述が簡潔なアルゴリズムとして提案された解法が, [[アフィン変換法]](affine scaling method) である. その後1988年に, ロシアの数学者ディキン(I. Dikin) が20年以上前の1967年に同じ解法を提案していたことが明かになり, 話題となった. アフィン変換法は実用的であり, 大域的収束性が示されているが, 解法の多項式時間性については未だ不明である.&lt;br /&gt;
&lt;br /&gt;
　他方, 従来のSUMT法, ホモトピー法, ニュートン法といった枠組を用いて内点法を構築する試みも行なわれた. 中でも大きな役割を果たしたのは, ゾンネべンド(Gy. Sonnevend)が提唱した[[中心パス]](path of centers) であり, この概念を利用してレネガー(J. Renegar)はニュートン法を用いた初めての 多項式時間解法を提案した. さらに, 主問題と双対問題の双方を1つの問題とみなした場合の解析的 中心に関する研究も行なわれ, [[主双対内点法]](primal-dual interior point method) の提案へと結び付いた. 主双対内点法を改良した[[予測子修正子内点法]](predictor-corrector interior point method) は, 現在もっとも普及している内点法の1つである. これらの解法については [5] 等を参照されたい.&lt;br /&gt;
&lt;br /&gt;
　内点法における大きな課題の1つは, 条件のよい許容領域内部の初期点をいかに用意するかという問題である. この課題の解決策として [[非許容初期点内点法]] (infeasible interior point method) や, 理論的にはさらに優れた [[同次自己双対内点法]] (homogeneous self-dual interior point method) が提案されており, 実装の際にはどちらかが用いられることが多い.&lt;br /&gt;
&lt;br /&gt;
　以上の解法の中で, アフィン変換法以外の解法は多項式時間解法であることが示されている. ネステロフとネミロフスキーは, [[自己整合障壁関数]]の概念を導入した上で, 内点法という解法のクラスの中で, 多項式時間解法である解法からなる, 1つの大きなクラスを示した [2]. この研究は, 1990年代から活発に行なわれている, [[半正定値計画]]問題や[[2次錐計画]]問題等に内点法を適用する研究の重要な基礎となっている.&lt;br /&gt;
&lt;br /&gt;
　近年内点法が記述された教科書も多く出版されている. [12] は内点法を総合的に扱った和書である. 線形計画法のみならず, 半正定値計画問題, 非線形計画問題に対する内点法についも述べている. その他の洋書の多くは主双対内点法と非許容初期点内点法を基本とした内容であるが, カーマーカー法については [3], [9], アフィン変換法については [4], [6], 同次自己双対内点法については, [6], [10] 等にそれぞれ詳しく述べられている. [8], [9], [10] はそれぞれ線形計画法の教科書であり, 図表が多く入門的な内容であるが, [8] は単体法, [9] は内点法を包括する理論, [10] は実装 にそれぞれ重きを置いている点が特徴である. また, 内点法に関する論文の多くは, インターネットを通じて入手可能である [11].&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
----&lt;br /&gt;
'''参考文献'''&lt;br /&gt;
&lt;br /&gt;
[1] K. Karmarkar, &amp;quot;A New Polynomial-Time Algorithm for Linear Programming,&amp;quot; ''Combinatorica'', '''4''' (1984), 373-395.&lt;br /&gt;
&lt;br /&gt;
[2] Y. Nesterov and A. S. Nemirovskii, ''Interior Point Polynomial Algorithms in Convex Programming'', SIAM, 1994.&lt;br /&gt;
&lt;br /&gt;
[3] S.-C. Fang and S. Puthenpra, ''Linear Optimization and Extentions'' Prentice Hall, 1993.&lt;br /&gt;
&lt;br /&gt;
[4] R. Saigal, ''Linear Programming: A Modern Integrated Analysis'' Kluwer Academic Publishers, 1995.&lt;br /&gt;
&lt;br /&gt;
[5] J. Wright, ''Primal-Dual Interior-Point Methods'', SIAM, 1996.&lt;br /&gt;
&lt;br /&gt;
[6] T. Terlaky, eds., ''Interior Point Methods of Mathematical Programming'', Kluwer Academic Publishers, 1996.&lt;br /&gt;
&lt;br /&gt;
[7] Y. Ye, ''Interior Point Algorithms: Theory and Analysis'', Jhon Wiley &amp;amp; Sons, 1997&lt;br /&gt;
&lt;br /&gt;
[8] D. Bertsimas and J. N. Tsisiklis, ''Introduction to Linear Optimization'', Athena Scientific, 1997.&lt;br /&gt;
&lt;br /&gt;
[9] C. Roos, T. Terlaky and J.-Ph. Vial ''Theory and Algorithns for Linear Optimization'', John Wiley &amp;amp; Sons, 1997.&lt;br /&gt;
&lt;br /&gt;
[10] R. J. Vanderbei, ''Linear Programming: Foundations and Extensions'', Kluwer Academic Publishers, 1998.&lt;br /&gt;
&lt;br /&gt;
[11] http://www-unix.mcs.anl.gov/otc/InteriorPoint&lt;br /&gt;
&lt;br /&gt;
[12] 小島政和, 土谷隆, 水野眞治, 矢部博, 『内点法』, 朝倉書店, 2001.&lt;br /&gt;
&lt;br /&gt;
[[Category:線形計画|ないてんほう]]&lt;/div&gt;</summary>
		<author><name>Tamura</name></author>
	</entry>
	<entry>
		<id>https://orsj-ml.org/orwiki/wiki/index.php?title=%E5%86%85%E7%82%B9%E6%B3%95&amp;diff=9992</id>
		<title>内点法</title>
		<link rel="alternate" type="text/html" href="https://orsj-ml.org/orwiki/wiki/index.php?title=%E5%86%85%E7%82%B9%E6%B3%95&amp;diff=9992"/>
		<updated>2008-08-22T07:19:07Z</updated>

		<summary type="html">&lt;p&gt;Tamura: &lt;/p&gt;
&lt;hr /&gt;
&lt;div&gt;'''【ないてんほう (interior point method)】'''&lt;br /&gt;
&lt;br /&gt;
=== 概要 ===&lt;br /&gt;
線形計画問題に対するカーマーカー法によって触発され,発展した制約付き最適化問題の反復解法の総称. 実行可能領域の内部を通って最適解へ近づくのでこう呼ばれる.特に大規模問題に有効とされ, [[凸計画問題]], [[半正定値計画問題]]等へも拡張されている. 主問題または双対問題のどちらか一方の実行可能領域に点列を生成する内点法を主内点法, 主問題および双対問題の双方の実行可能領域に点列を生成する内点法を主双対内点法と呼ぶ.&lt;br /&gt;
&lt;br /&gt;
=== 詳説 ===&lt;br /&gt;
　[[内点法]](interior point method) は[[最適化問題]]の制約領域の境界上ではなく内部に, 最適解に収束する点列を生成する逐次反復解法である. 特に[[線形計画|線形計画問題]]に対しては, 最適値に十分に近い近似解が得られれば, 変数次元の連立方程式を1回解く手間で最適解が得られることから, 内点法で最適解を得るまでに必要な計算の手間(算術演算回数)を算出することが可能であり, 多くの理論的な成果が報告されている. また実用面では, 特に大規模な問題に対して単体法よりも計算効率が 優れていることが確認されており, この解法を取り入れた商用コードも多く存在する. 線形計画問題以外への応用も積極的に行なわれており, 凸2次計画問題, [[相補性問題]], あるいは半正定値計画問題, 2次錐計画問題等の凸計画問題に対しても有効な解法であることが示されている.&lt;br /&gt;
&lt;br /&gt;
　1979年に提案された楕円体法は, 線形計画問題に対する初めての[[多項式時間解法]]であるが, 実用的な解法としては単体法に及ばなかった. これに対して1984年にATTベル研究所のカーマーカー(N. Karmarkar) によって提案された初めての内点法(interior point)である[[カーマーカー法]](Karmarkar's algorithm) [1] は, 多項式時間解法である上に, カーマーカー自身の計算機実験では 単体法を遥かに上回る結果であると報告された. この論文および報告を機として, 多くの最適化分野の研究者がこの解法に興味をもち, この結果多様なバリエーションが生み出された.&lt;br /&gt;
&lt;br /&gt;
　カーマーカー法は, 特殊な線形計画問題&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;\begin{array}{llllll}&lt;br /&gt;
\mbox{max.} &amp;amp; c^{\top}x  &amp;amp; \mbox{s. t.}  &amp;amp; Ax = 0,&lt;br /&gt;
&amp;amp; e^{\top}x = 1,  &amp;amp; x \geq 0, &lt;br /&gt;
\end{array}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
(&amp;lt;math&amp;gt;A\,&amp;lt;/math&amp;gt;は&amp;lt;math&amp;gt;m \times n\,&amp;lt;/math&amp;gt;行列, &amp;lt;math&amp;gt;c \in \mathbf{R}^n&amp;lt;/math&amp;gt;, &amp;lt;math&amp;gt;e \in \mathbf{R}^n&amp;lt;/math&amp;gt;は要素がすべて&amp;lt;math&amp;gt;1\,&amp;lt;/math&amp;gt;のベクトル)を対象とし, 対数障壁関数を用いた[[ポテンシャル関数 (内点法の)|ポテンシャル関数]](potential function) の値で最適解からの誤差を測り, [[射影変換]] (projective transformation) を用いて探索方向を決定するという特徴をもつ. 対数障壁関数を用いるという点で非線形計画法の一解法とも考えられるが, 着想, 解析共に従来の常識とは異なる斬新な解法であった. より一般的で記述が簡潔なアルゴリズムとして提案された解法が, [[アフィン変換法]](affine scaling method) である. その後1988年に, ロシアの数学者ディキン(I. Dikin) が20年以上前の1967年に同じ解法を提案していたことが明かになり, 話題となった. アフィン変換法は実用的であり, 大域的収束性が示されているが, 解法の多項式時間性については未だ不明である.&lt;br /&gt;
&lt;br /&gt;
　他方, 従来のSUMT法, ホモトピー法, ニュートン法といった枠組を用いて内点法を構築する試みも行なわれた. 中でも大きな役割を果たしたのは, ゾンネべンド(Gy. Sonnevend)が提唱した[[中心パス]](path of centers) であり, この概念を利用してレネガー(J. Renegar)はニュートン法を用いた初めての 多項式時間解法を提案した. さらに, 主問題と双対問題の双方を1つの問題とみなした場合の解析的 中心に関する研究も行なわれ, [[主双対内点法]](primal-dual interior point method) の提案へと結び付いた. 主双対内点法を改良した[[予測子修正子内点法]](predictor-corrector interior point method) は, 現在もっとも普及している内点法の1つである. これらの解法については [5] 等を参照されたい.&lt;br /&gt;
&lt;br /&gt;
　内点法における大きな課題の1つは, 条件のよい許容領域内部の初期点をいかに用意するかという問題である. この課題の解決策として [[非許容初期点内点法]] (infeasible interior point method) や, 理論的にはさらに優れた [[同次自己双対内点法]] (homogeneous self-dual interior point method) が提案されており, 実装の際にはどちらかが用いられることが多い.&lt;br /&gt;
&lt;br /&gt;
　以上の解法の中で, アフィン変換法以外の解法は多項式時間解法であることが示されている. ネステロフとネミロフスキーは, [[自己整合障壁関数]]の概念を導入した上で, 内点法という解法のクラスの中で, 多項式時間解法である解法からなる, 1つの大きなクラスを示した [2]. この研究は, 1990年代から活発に行なわれている, [[半正定値計画]]問題や[[2次錐計画]]問題等に内点法を適用する研究の重要な基礎となっている.&lt;br /&gt;
&lt;br /&gt;
　近年内点法が記述された教科書も多く出版されている. [12] は内点法を総合的に扱った和書である. 線形計画法のみならず, 半正定値計画問題, 非線形計画問題に対する内点法についも述べている. その他の洋書の多くは主双対内点法と非許容初期点内点法を基本とした内容であるが, カーマーカー法については [3], [9], アフィン変換法については [4], [6], 同次自己双対内点法については, [6], [10] 等にそれぞれ詳しく述べられている. [8], [9], [10] はそれぞれ線形計画法の教科書であり, 図表が多く入門的な内容であるが, [8] は単体法, [9] は内点法を包括する理論, [10] は実装 にそれぞれ重きを置いている点が特徴である. また, 内点法に関する論文の多くは, インターネットを通じて入手可能である [11].&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
----&lt;br /&gt;
'''参考文献'''&lt;br /&gt;
&lt;br /&gt;
[1] K. Karmarkar, &amp;quot;A New Polynomial-Time Algorithm for Linear Programming,&amp;quot; ''Combinatorica'', '''4''' (1984), 373-395.&lt;br /&gt;
&lt;br /&gt;
[2] Y. Nesterov and A. S. Nemirovskii, ''Interior Point Polynomial Algorithms in Convex Programming'', SIAM, 1994.&lt;br /&gt;
&lt;br /&gt;
[3] S.-C. Fang and S. Puthenpra, ''Linear Optimization and Extentions'' Prentice Hall, 1993.&lt;br /&gt;
&lt;br /&gt;
[4] R. Saigal, ''Linear Programming: A Modern Integrated Analysis'' Kluwer Academic Publishers, 1995.&lt;br /&gt;
&lt;br /&gt;
[5] J. Wright, ''Primal-Dual Interior-Point Methods'', SIAM, 1996.&lt;br /&gt;
&lt;br /&gt;
[6] T. Terlaky, eds., ''Interior Point Methods of Mathematical Programming'', Kluwer Academic Publishers, 1996.&lt;br /&gt;
&lt;br /&gt;
[7] Y. Ye, ''Interior Point Algorithms: Theory and Analysis'', Jhon Wiley &amp;amp; Sons, 1997&lt;br /&gt;
&lt;br /&gt;
[8] D. Bertsimas and J. N. Tsisiklis, ''Introduction to Linear Optimization'', Athena Scientific, 1997.&lt;br /&gt;
&lt;br /&gt;
[9] C. Roos, T. Terlaky and J.-Ph. Vial ''Theory and Algorithns for Linear Optimization'', John Wiley &amp;amp; Sons, 1997.&lt;br /&gt;
&lt;br /&gt;
[10] R. J. Vanderbei, ''Linear Programming: Foundations and Extensions'', Kluwer Academic Publishers, 1998.&lt;br /&gt;
&lt;br /&gt;
[11] http://www-unix.mcs.anl.gov/otc/InteriorPoint&lt;br /&gt;
&lt;br /&gt;
[12] 小島政和, 土谷隆, 水野眞治, 矢部博, 『内点法』, 朝倉書店, 2001.&lt;br /&gt;
&lt;br /&gt;
[[Category:線形計画|ないてんほう]]&lt;/div&gt;</summary>
		<author><name>Tamura</name></author>
	</entry>
	<entry>
		<id>https://orsj-ml.org/orwiki/wiki/index.php?title=%E6%A5%95%E5%86%86%E4%BD%93%E6%B3%95&amp;diff=9991</id>
		<title>楕円体法</title>
		<link rel="alternate" type="text/html" href="https://orsj-ml.org/orwiki/wiki/index.php?title=%E6%A5%95%E5%86%86%E4%BD%93%E6%B3%95&amp;diff=9991"/>
		<updated>2008-08-22T07:14:40Z</updated>

		<summary type="html">&lt;p&gt;Tamura: /* 詳説 */&lt;/p&gt;
&lt;hr /&gt;
&lt;div&gt;'''【だえんたいほう (ellipsoid method)】'''&lt;br /&gt;
&lt;br /&gt;
=== 概要 ===&lt;br /&gt;
　[[楕円体法]](ellipsoid method) は, 微分不可能な[[凸計画問題]]に対する解法として, 1976年にユーディンとネミロフスキー(Yudin-Nemirovskii) によって提案された.  その後1979年にカチヤン(L. G. Khachiyan)は, 楕円体法が[[線形計画|線形計画問題]]に対する[[多項式時間解法]] に成り得ることを示した.  この結果は, 長年にわたって未解決のままであった,「線形計画問題を多項式時間で解くことが出来るか?」という計算複雑度に関する問題への肯定的な答を与えるものであり, 数理計画の分野に衝撃をもたらした.  線形計画法に対する解法としての楕円体法は, 実用の面からは[[単体法]]などに比べて効率が悪く, また理論的な計算量の面からも, 1984年以降に現れた[[内点法]]に劣っていることが分かっている. しかしながら, 線形, 非線形, 組合せなどの最適化問題に対し, 様々な理論的な結果を残している. 「&amp;lt;math&amp;gt;n&amp;lt;/math&amp;gt; 次元空間における楕円体を用いた2分探索」のような解法である.&lt;br /&gt;
&lt;br /&gt;
=== 詳説 ===&lt;br /&gt;
　楕円体法を説明するために, 与えられた多面体 &amp;lt;math&amp;gt;P = \{x \in \mathbf{R}^n \mid a_i^{\top} x \leq b_i\ (i = 1, 2, \cdots, m)\}&amp;lt;/math&amp;gt; が空であるか否かを判定し, また空でなければ &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; に含まれる点を求める, という非空性問題について考える. なお, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; の上での線形関数 &amp;lt;math&amp;gt;c^{\top} x&amp;lt;/math&amp;gt; の最大化問題を解くには, &amp;lt;math&amp;gt;\{x \in P \mid c^{\top} x \leq \gamma\}&amp;lt;/math&amp;gt; という多面体に対する非空性問題を繰り返し解き, &amp;lt;math&amp;gt;\gamma \in \mathbf{R}&amp;lt;/math&amp;gt; に関する2分探索を行えば良い.&lt;br /&gt;
&lt;br /&gt;
　この問題に対し, 以下の2点を仮定する.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
:&amp;lt;math&amp;gt;\mbox{(i)} \quad  P \subseteq \{x \in \mathbf{R}^n \mid \|x - x_0\| \leq R\}&amp;lt;/math&amp;gt;なる &amp;lt;math&amp;gt;x_0 \in \mathbf{R}^n&amp;lt;/math&amp;gt; 及び &amp;lt;math&amp;gt;R \in \mathbf{R}&amp;lt;/math&amp;gt; が既知である.&lt;br /&gt;
&lt;br /&gt;
:&amp;lt;math&amp;gt;\mbox{(ii)} \quad P\, &amp;lt;/math&amp;gt; が空でないならば, その体積は &amp;lt;math&amp;gt;\varepsilon (&amp;gt; 0)&amp;lt;/math&amp;gt; 以上である.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
　楕円体法では, いわば 「&amp;lt;math&amp;gt;n\,&amp;lt;/math&amp;gt;次元空間での [[楕円体]]を用いた2分探索」を行うことにより, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt;の点を求める. 各反復では, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt;を含む楕円体&amp;lt;math&amp;gt;E_k\,&amp;lt;/math&amp;gt; を常に保持する. 楕円体 &amp;lt;math&amp;gt;E_k\,&amp;lt;/math&amp;gt;は, 中心と呼ばれるベクトル &amp;lt;math&amp;gt;x_k \in \mathbf{R}^n&amp;lt;/math&amp;gt;及び&amp;lt;math&amp;gt;n \times n&amp;lt;/math&amp;gt;正定値対称行列&amp;lt;math&amp;gt;B_k\,&amp;lt;/math&amp;gt;を用いて,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&amp;lt;math&amp;gt;E_k = \{x \in \mathbf{R}^n \mid (x - x_k)^{\top} B_k^{-1}(x - x_k) \leq 1\}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
と表される. なお, 上記の仮定 (i) より, &amp;lt;math&amp;gt;B_0 = R^2 I&amp;lt;/math&amp;gt; と選ぶことができる.&lt;br /&gt;
&lt;br /&gt;
　各反復では, &amp;lt;math&amp;gt;x_k \in P\,&amp;lt;/math&amp;gt; が成り立つか否かを調べる.  成り立つ場合はアルゴリズムを停止する.  成り立たない場合は, &amp;lt;math&amp;gt;a_i^{\top} x_k &amp;gt; b_i&amp;lt;/math&amp;gt; なる &amp;lt;math&amp;gt;i\,&amp;lt;/math&amp;gt; が存在するので, そのような &amp;lt;math&amp;gt;i\,&amp;lt;/math&amp;gt; を見つける.  中心 &amp;lt;math&amp;gt;x_k\,&amp;lt;/math&amp;gt; を通る超平面 &amp;lt;math&amp;gt;a_i^{\top} x = a_i^{\top} x_k&amp;lt;/math&amp;gt; により,  &amp;lt;math&amp;gt;E_k\,&amp;lt;/math&amp;gt; を半分にして得られる集合 &amp;lt;math&amp;gt;E_k' = \{x \in E_k \mid a_i^{\top} x \leq a_i^{\top} x_k\}&amp;lt;/math&amp;gt;  に対し, これを含み, かつ体積が最小の楕円体を &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; とおく.  そのような楕円体 &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; を表す &amp;lt;math&amp;gt;x_{k+1}\,&amp;lt;/math&amp;gt; 及び &amp;lt;math&amp;gt;B_{k+1}\,&amp;lt;/math&amp;gt; は 次のようにして求められる:&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;x_{k+1} = x_k - \frac{1}{n+1}d,\qquad&lt;br /&gt;
B_{k+1} = \frac{n^2}{n^2 - 1}\left(B_k - \frac{2}{n+1}d d^{\top} \right).&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
ここで, &amp;lt;math&amp;gt;d = (1/\sqrt{a_i^{\top}B_k a_i})B_k a_i&amp;lt;/math&amp;gt;である.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
　上記のように &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; を定めると, 楕円体の体積の比は,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;\frac{{\rm vol}(E_{k+1})}{{\rm vol}(E_{k})}&lt;br /&gt;
 = \left(\left(\frac{n}{n+1}\right)^{n+1}&lt;br /&gt;
 \left(\frac{n}{n-1}\right)^{n-1}\right)^{1/2}&lt;br /&gt;
 &amp;lt; e^{-1/(2n)} &amp;lt; 1&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
となる.  したがって, 楕円体の体積は線形空間の次元 &amp;lt;math&amp;gt;n\,&amp;lt;/math&amp;gt; のみに依存する一定の比率で 減少していく.  反復回数 &amp;lt;math&amp;gt;k\,&amp;lt;/math&amp;gt; が十分大きくなり 楕円体の体積が &amp;lt;math&amp;gt;\varepsilon&amp;lt;/math&amp;gt; 未満になった とすると, 仮定 (ii) により &amp;lt;math&amp;gt;P = \emptyset&amp;lt;/math&amp;gt; であることが分かるので, アルゴリズムを停止する.  &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; が多面体であることから, 仮定 (i), (ii) での &amp;lt;math&amp;gt;x_0\,&amp;lt;/math&amp;gt;, &amp;lt;math&amp;gt;R\,&amp;lt;/math&amp;gt;, &amp;lt;math&amp;gt;\varepsilon&amp;lt;/math&amp;gt; をうまく選ぶことができ, それにより多項式回の反復で終了することが導かれる.&lt;br /&gt;
&lt;br /&gt;
　上記で述べた楕円体法の改良として,「深い」(deep)カットを用いた方法が 提案されている.  この方法では, 半楕円体 &amp;lt;math&amp;gt;E_k'\,&amp;lt;/math&amp;gt; の代わりに, より小さな集合 &amp;lt;math&amp;gt;\{x \in E_k \mid a_i^{\top} x \leq b_i\}&amp;lt;/math&amp;gt; に注目し, それを含む最小の楕円体 &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; へと更新していく.  この場合, 連続する2つの楕円体の比は, &amp;lt;math&amp;gt;\alpha = (a_i^{\top} x_k - b_i)/\sqrt{a_i^{\top}B_k a_i}&amp;lt;/math&amp;gt; とおくと,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;\frac{{\rm vol}(E_{k+1})}{{\rm vol}(E_{k})}&lt;br /&gt;
 = \left(\left(\frac{n(1 - \alpha)}{n+1}\right)^{n+1}&lt;br /&gt;
 \left(\frac{n (1 + \alpha)}{n-1}\right)^{n-1}\right)^{1/2}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
と表され, また &amp;lt;math&amp;gt;0 \leq \alpha \leq 1&amp;lt;/math&amp;gt; なので, 元の楕円体法に比べて反復回数が減少することがわかる.&lt;br /&gt;
&lt;br /&gt;
　さらに, 代用(surrogate)カット, 平行(parallel)カットを用いれば, 楕円体の体積比をより減少させることができる. また, 深いカットと全く逆の考え方により得られる「浅い」(shallow)カットを用いると, 連続する楕円体の体積比は逆に増加するが, 理論的には, 他のカットでは得られない興味深い結果を得ることができる. これらのカットについての詳細は参考文献を参照されたい.&lt;br /&gt;
&lt;br /&gt;
　線形計画法に対する楕円体法の適用は, 線形計画のみならず組合せ最適化においても重要な結果を残している. その1つが, 最適化問題と[[分離問題]]の多項式時間に関する等価性であろう.&lt;br /&gt;
&lt;br /&gt;
　&amp;lt;math&amp;gt;P \subseteq \mathbf{R}^n&amp;lt;/math&amp;gt; を多面体とする.  ただし, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; の不等式系による表現は必ずしも与えられていない, と仮定する.  分離問題とは, 任意のベクトル &amp;lt;math&amp;gt;x_* \in \mathbf{R}^n&amp;lt;/math&amp;gt; が与えられたとき,  &amp;lt;math&amp;gt;x_* \in P&amp;lt;/math&amp;gt; か否かを判定し, &amp;lt;math&amp;gt;x_* \not\in P&amp;lt;/math&amp;gt; ならば &amp;lt;math&amp;gt;a^{\top} x \leq b&amp;lt;/math&amp;gt; &amp;lt;math&amp;gt;(\forall x \in P)&amp;lt;/math&amp;gt; かつ &amp;lt;math&amp;gt;a^{\top} x_* &amp;gt; b&amp;lt;/math&amp;gt; なる &amp;lt;math&amp;gt;a \in \mathbf{R}^n&amp;lt;/math&amp;gt; 及び &amp;lt;math&amp;gt;b \in \mathbf{R}&amp;lt;/math&amp;gt; を求める問題のことである.  上記で述べた楕円体法のアルゴリズムを見ると分かるが, 必ずしも &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; の不等式表現が陽に与えられている必要はなく, 各反復において &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; と &amp;lt;math&amp;gt;x_k\,&amp;lt;/math&amp;gt; に関する分離問題を多項式時間で解けるならば, 最適化問題もまた多項式時間で解くことができる.  一方, この逆もまた成り立つことが Gr&amp;amp;ouml;tschel-Lov&amp;amp;aacute;sz-Schrijver (1981) 及び Padberg-Rao (1981) によって証明されている.&lt;br /&gt;
&lt;br /&gt;
　この他にも, [[劣モジュラ関数]]の最小化や[[パーフェクトグラフ]]での最大 安定集合問題などの組合せ最適化問題が多項式時間で解ける, という結果が, 楕円体法の適用により導かれている.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
----&lt;br /&gt;
'''参考文献'''&lt;br /&gt;
&lt;br /&gt;
[1] R. G. Bland, D. Goldfarb and M. J. Todd, &amp;quot;The ellipsoid method: a survey,&amp;quot; ''Operations Research'', '''29'''(1981) 1039-1091.&lt;br /&gt;
&lt;br /&gt;
[2] M. Gr&amp;amp;ouml;tschel, L. Lov&amp;amp;aacute;sz and A. Schrijver, ''Geometric Algorithms and Combinatorial Optimization'', Springer-Verlag, 1991.&lt;br /&gt;
&lt;br /&gt;
[3] 伊理正夫,『線形計画法』, 共立出版, 1986.&lt;br /&gt;
&lt;br /&gt;
[4] G. L. Nemhauser, A. H. G. Rinnooy Kan and M. J. Todd, eds., ''Optimization'', Handbooks in Operations Research and Management Science Vol. 1, North-Holland, 1989. 伊理正夫, 今野浩, 刀根薫 監訳,『最適化ハンドブック』,朝倉書店, 1995.&lt;br /&gt;
&lt;br /&gt;
[[Category:線形計画|だえんたいほう]]&lt;/div&gt;</summary>
		<author><name>Tamura</name></author>
	</entry>
	<entry>
		<id>https://orsj-ml.org/orwiki/wiki/index.php?title=%E6%A5%95%E5%86%86%E4%BD%93%E6%B3%95&amp;diff=9990</id>
		<title>楕円体法</title>
		<link rel="alternate" type="text/html" href="https://orsj-ml.org/orwiki/wiki/index.php?title=%E6%A5%95%E5%86%86%E4%BD%93%E6%B3%95&amp;diff=9990"/>
		<updated>2008-08-22T07:13:17Z</updated>

		<summary type="html">&lt;p&gt;Tamura: /* 詳説 */&lt;/p&gt;
&lt;hr /&gt;
&lt;div&gt;'''【だえんたいほう (ellipsoid method)】'''&lt;br /&gt;
&lt;br /&gt;
=== 概要 ===&lt;br /&gt;
　[[楕円体法]](ellipsoid method) は, 微分不可能な[[凸計画問題]]に対する解法として, 1976年にユーディンとネミロフスキー(Yudin-Nemirovskii) によって提案された.  その後1979年にカチヤン(L. G. Khachiyan)は, 楕円体法が[[線形計画|線形計画問題]]に対する[[多項式時間解法]] に成り得ることを示した.  この結果は, 長年にわたって未解決のままであった,「線形計画問題を多項式時間で解くことが出来るか?」という計算複雑度に関する問題への肯定的な答を与えるものであり, 数理計画の分野に衝撃をもたらした.  線形計画法に対する解法としての楕円体法は, 実用の面からは[[単体法]]などに比べて効率が悪く, また理論的な計算量の面からも, 1984年以降に現れた[[内点法]]に劣っていることが分かっている. しかしながら, 線形, 非線形, 組合せなどの最適化問題に対し, 様々な理論的な結果を残している. 「&amp;lt;math&amp;gt;n&amp;lt;/math&amp;gt; 次元空間における楕円体を用いた2分探索」のような解法である.&lt;br /&gt;
&lt;br /&gt;
=== 詳説 ===&lt;br /&gt;
　楕円体法を説明するために, 与えられた多面体 &amp;lt;math&amp;gt;P = \{x \in \mathbf{R}^n \mid a_i^{\top} x \leq b_i\ (i = 1, 2, \cdots, m)\}&amp;lt;/math&amp;gt; が空であるか否かを判定し, また空でなければ &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; に含まれる点を求める, という非空性問題について考える. なお, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; の上での線形関数 &amp;lt;math&amp;gt;c^{\top} x&amp;lt;/math&amp;gt; の最大化問題を解くには, &amp;lt;math&amp;gt;\{x \in P \mid c^{\top} x \leq \gamma\}&amp;lt;/math&amp;gt; という多面体に対する非空性問題を繰り返し解き, &amp;lt;math&amp;gt;\gamma \in \mathbf{R}&amp;lt;/math&amp;gt; に関する2分探索を行えば良い.&lt;br /&gt;
&lt;br /&gt;
　この問題に対し, 以下の2点を仮定する.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
:&amp;lt;math&amp;gt;\mbox{(i)} \quad  P \subseteq \{x \in \mathbf{R}^n \mid \|x - x_0\| \leq R\}&amp;lt;/math&amp;gt;なる &amp;lt;math&amp;gt;x_0 \in \mathbf{R}^n&amp;lt;/math&amp;gt; 及び &amp;lt;math&amp;gt;R \in \mathbf{R}&amp;lt;/math&amp;gt; が既知である.&lt;br /&gt;
&lt;br /&gt;
:&amp;lt;math&amp;gt;\mbox{(ii)} \quad P\, &amp;lt;/math&amp;gt; が空でないならば, その体積は &amp;lt;math&amp;gt;\varepsilon (&amp;gt; 0)&amp;lt;/math&amp;gt; 以上である.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
　楕円体法では, いわば 「&amp;lt;math&amp;gt;n\,&amp;lt;/math&amp;gt;次元空間での [[楕円体]]を用いた2分探索」を行うことにより, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt;の点を求める. 各反復では, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt;を含む楕円体&amp;lt;math&amp;gt;E_k\,&amp;lt;/math&amp;gt; を常に保持する. 楕円体 &amp;lt;math&amp;gt;E_k\,&amp;lt;/math&amp;gt;は, 中心と呼ばれるベクトル &amp;lt;math&amp;gt;x_k \in \mathbf{R}^n&amp;lt;/math&amp;gt;及び&amp;lt;math&amp;gt;n \times n&amp;lt;/math&amp;gt;正定値対称行列&amp;lt;math&amp;gt;B_k\,&amp;lt;/math&amp;gt;を用いて,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&amp;lt;math&amp;gt;E_k = \{x \in \mathbf{R}^n \mid (x - x_k)^{\top} B_k^{-1}(x - x_k) \leq 1\}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
と表される. なお, 上記の仮定 (i) より, &amp;lt;math&amp;gt;B_0 = R^2 I&amp;lt;/math&amp;gt; と選ぶことができる.&lt;br /&gt;
&lt;br /&gt;
　各反復では, &amp;lt;math&amp;gt;x_k \in P\,&amp;lt;/math&amp;gt; が成り立つか否かを調べる.  成り立つ場合はアルゴリズムを停止する.  成り立たない場合は, &amp;lt;math&amp;gt;a_i^{\top} x_k &amp;gt; b_i&amp;lt;/math&amp;gt; なる &amp;lt;math&amp;gt;i\,&amp;lt;/math&amp;gt; が存在するので, そのような &amp;lt;math&amp;gt;i\,&amp;lt;/math&amp;gt; を見つける.  中心 &amp;lt;math&amp;gt;x_k\,&amp;lt;/math&amp;gt; を通る超平面 &amp;lt;math&amp;gt;a_i^{\top} x = a_i^{\top} x_k&amp;lt;/math&amp;gt; により,  &amp;lt;math&amp;gt;E_k\,&amp;lt;/math&amp;gt; を半分にして得られる集合 &amp;lt;math&amp;gt;E_k' = \{x \in E_k \mid a_i^{\top} x \leq a_i^{\top} x_k\}&amp;lt;/math&amp;gt;  に対し, これを含み, かつ体積が最小の楕円体を &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; とおく.  そのような楕円体 &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; を表す &amp;lt;math&amp;gt;x_{k+1}\,&amp;lt;/math&amp;gt; 及び &amp;lt;math&amp;gt;B_{k+1}\,&amp;lt;/math&amp;gt; は 次のようにして求められる:&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;x_{k+1} = x_k - \frac{1}{n+1}d,\qquad&lt;br /&gt;
B_{k+1} = \frac{n^2}{n^2 - 1}\left(B_k - \frac{2}{n+1}d d^{\top} \right).&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
ここで, &amp;lt;math&amp;gt;d = (1/\sqrt{a_i^{\top}B_k a_i})B_k a_i&amp;lt;/math&amp;gt;である.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
　上記のように &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; を定めると, 楕円体の体積の比は,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;\frac{{\rm vol}(E_{k+1})}{{\rm vol}(E_{k})}&lt;br /&gt;
 = \left(\left(\frac{n}{n+1}\right)^{n+1}&lt;br /&gt;
 \left(\frac{n}{n-1}\right)^{n-1}\right)^{1/2}&lt;br /&gt;
 &amp;lt; e^{-1/(2n)} &amp;lt; 1&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
となる.  したがって, 楕円体の体積は線形空間の次元 &amp;lt;math&amp;gt;n\,&amp;lt;/math&amp;gt; のみに依存する一定の比率で 減少していく.  反復回数 &amp;lt;math&amp;gt;k\,&amp;lt;/math&amp;gt; が十分大きくなり 楕円体の体積が &amp;lt;math&amp;gt;\varepsilon&amp;lt;/math&amp;gt; 未満になった とすると, 仮定 (ii) により &amp;lt;math&amp;gt;P = \emptyset&amp;lt;/math&amp;gt; であることが分かるので, アルゴリズムを停止する.  &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; が多面体であることから, 仮定 (i), (ii) での &amp;lt;math&amp;gt;x_0\,&amp;lt;/math&amp;gt;, &amp;lt;math&amp;gt;R\,&amp;lt;/math&amp;gt;, &amp;lt;math&amp;gt;\varepsilon&amp;lt;/math&amp;gt; をうまく選ぶことができ, それにより多項式回の反復で終了することが導かれる.&lt;br /&gt;
&lt;br /&gt;
　上記で述べた楕円体法の改良として,「深い」(deep)カットを用いた方法が 提案されている.  この方法では, 半楕円体 &amp;lt;math&amp;gt;E_k'\,&amp;lt;/math&amp;gt; の代わりに, より小さな集合 &amp;lt;math&amp;gt;\{x \in E_k \mid a_i^{\top} x \leq b_i\}&amp;lt;/math&amp;gt; に注目し, それを含む最小の楕円体 &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; へと更新していく.  この場合, 連続する2つの楕円体の比は, &amp;lt;math&amp;gt;\alpha = (a_i^{\top} x_k - b_i)/\sqrt{a_i^{\top}B_k a_i}&amp;lt;/math&amp;gt; とおくと,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;\frac{{\rm vol}(E_{k+1})}{{\rm vol}(E_{k})}&lt;br /&gt;
 = \left(\left(\frac{n(1 - \alpha)}{n+1}\right)^{n+1}&lt;br /&gt;
 \left(\frac{n (1 + \alpha)}{n-1}\right)^{n-1}\right)^{1/2}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
と表され, また &amp;lt;math&amp;gt;0 \leq \alpha \leq 1&amp;lt;/math&amp;gt; なので, 元の楕円体法に比べて反復回数が減少することがわかる.&lt;br /&gt;
&lt;br /&gt;
　さらに, 代用(surrogate)カット, 平行(parallel)カットを用いれば, 楕円体の体積比をより減少させることができる. また, 深いカットと全く逆の考え方により得られる「浅い」(shallow)カットを用いると, 連続する楕円体の体積比は逆に増加するが, 理論的には, 他のカットでは得られない興味深い結果を得ることができる. これらのカットについての詳細は参考文献を参照されたい.&lt;br /&gt;
&lt;br /&gt;
　線形計画法に対する楕円体法の適用は, 線形計画のみならず組合せ最適化においても重要な結果を残している. その1つが, 最適化問題と[[分離問題]]の多項式時間に関する等価性であろう.&lt;br /&gt;
&lt;br /&gt;
　&amp;lt;math&amp;gt;P \subseteq \mathbf{R}^n&amp;lt;/math&amp;gt; を多面体とする.  ただし, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; の不等式系による表現は必ずしも与えられていない, と仮定する.  分離問題とは, 任意のベクトル &amp;lt;math&amp;gt;x_* \in \mathbf{R}^n&amp;lt;/math&amp;gt; が与えられたとき,  &amp;lt;math&amp;gt;x_* \in P&amp;lt;/math&amp;gt; か否かを判定し, &amp;lt;math&amp;gt;x_* \not\in P&amp;lt;/math&amp;gt; ならば &amp;lt;math&amp;gt;a^{\top} x \leq b&amp;lt;/math&amp;gt; &amp;lt;math&amp;gt;(\forall x \in P)&amp;lt;/math&amp;gt; かつ &amp;lt;math&amp;gt;a^{\top} x_* &amp;gt; b&amp;lt;/math&amp;gt; なる &amp;lt;math&amp;gt;a \in \mathbf{R}^n&amp;lt;/math&amp;gt; 及び &amp;lt;math&amp;gt;b \in \mathbf{R}&amp;lt;/math&amp;gt; を求める問題のことである.  上記で述べた楕円体法のアルゴリズムを見ると分かるが, 必ずしも &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; の不等式表現が陽に与えられている必要はなく, 各反復において &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; と &amp;lt;math&amp;gt;x_k\,&amp;lt;/math&amp;gt; に関する分離問題を多項式時間で解けるならば, 最適化問題もまた多項式時間で解くことができる.  一方, この逆もまた成り立つことが Gr&amp;amp;ouml;tschel-Lov&amp;amp;aacute;sz-Schrijver (1981) 及び Padberg-Rao (1981) によって証明されている.&lt;br /&gt;
&lt;br /&gt;
　この他にも, 劣モジュラ関数の最小化やパーフェクトグラフでの最大 安定集合問題などの組合せ最適化問題が多項式時間で解ける, という結果が, 楕円体法の適用により導かれている.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
----&lt;br /&gt;
'''参考文献'''&lt;br /&gt;
&lt;br /&gt;
[1] R. G. Bland, D. Goldfarb and M. J. Todd, &amp;quot;The ellipsoid method: a survey,&amp;quot; ''Operations Research'', '''29'''(1981) 1039-1091.&lt;br /&gt;
&lt;br /&gt;
[2] M. Gr&amp;amp;ouml;tschel, L. Lov&amp;amp;aacute;sz and A. Schrijver, ''Geometric Algorithms and Combinatorial Optimization'', Springer-Verlag, 1991.&lt;br /&gt;
&lt;br /&gt;
[3] 伊理正夫,『線形計画法』, 共立出版, 1986.&lt;br /&gt;
&lt;br /&gt;
[4] G. L. Nemhauser, A. H. G. Rinnooy Kan and M. J. Todd, eds., ''Optimization'', Handbooks in Operations Research and Management Science Vol. 1, North-Holland, 1989. 伊理正夫, 今野浩, 刀根薫 監訳,『最適化ハンドブック』,朝倉書店, 1995.&lt;br /&gt;
&lt;br /&gt;
[[Category:線形計画|だえんたいほう]]&lt;/div&gt;</summary>
		<author><name>Tamura</name></author>
	</entry>
	<entry>
		<id>https://orsj-ml.org/orwiki/wiki/index.php?title=%E6%A5%95%E5%86%86%E4%BD%93%E6%B3%95&amp;diff=9989</id>
		<title>楕円体法</title>
		<link rel="alternate" type="text/html" href="https://orsj-ml.org/orwiki/wiki/index.php?title=%E6%A5%95%E5%86%86%E4%BD%93%E6%B3%95&amp;diff=9989"/>
		<updated>2008-08-22T07:11:39Z</updated>

		<summary type="html">&lt;p&gt;Tamura: /* 概要 */&lt;/p&gt;
&lt;hr /&gt;
&lt;div&gt;'''【だえんたいほう (ellipsoid method)】'''&lt;br /&gt;
&lt;br /&gt;
=== 概要 ===&lt;br /&gt;
　[[楕円体法]](ellipsoid method) は, 微分不可能な[[凸計画問題]]に対する解法として, 1976年にユーディンとネミロフスキー(Yudin-Nemirovskii) によって提案された.  その後1979年にカチヤン(L. G. Khachiyan)は, 楕円体法が[[線形計画|線形計画問題]]に対する[[多項式時間解法]] に成り得ることを示した.  この結果は, 長年にわたって未解決のままであった,「線形計画問題を多項式時間で解くことが出来るか?」という計算複雑度に関する問題への肯定的な答を与えるものであり, 数理計画の分野に衝撃をもたらした.  線形計画法に対する解法としての楕円体法は, 実用の面からは[[単体法]]などに比べて効率が悪く, また理論的な計算量の面からも, 1984年以降に現れた[[内点法]]に劣っていることが分かっている. しかしながら, 線形, 非線形, 組合せなどの最適化問題に対し, 様々な理論的な結果を残している. 「&amp;lt;math&amp;gt;n&amp;lt;/math&amp;gt; 次元空間における楕円体を用いた2分探索」のような解法である.&lt;br /&gt;
&lt;br /&gt;
=== 詳説 ===&lt;br /&gt;
　楕円体法を説明するために, 与えられた多面体 &amp;lt;math&amp;gt;P = \{x \in \mathbf{R}^n \mid a_i^{\top} x \leq b_i\ (i = 1, 2, \cdots, m)\}&amp;lt;/math&amp;gt; が空であるか否かを判定し, また空でなければ &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; に含まれる点を求める, という非空性問題について考える. なお, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; の上での線形関数 &amp;lt;math&amp;gt;c^{\top} x&amp;lt;/math&amp;gt; の最大化問題を解くには, &amp;lt;math&amp;gt;\{x \in P \mid c^{\top} x \leq \gamma\}&amp;lt;/math&amp;gt; という多面体に対する非空性問題を繰り返し解き, &amp;lt;math&amp;gt;\gamma \in \mathbf{R}&amp;lt;/math&amp;gt; に関する2分探索を行えば良い.&lt;br /&gt;
&lt;br /&gt;
　この問題に対し, 以下の2点を仮定する.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
:&amp;lt;math&amp;gt;\mbox{(i)} \quad  P \subseteq \{x \in \mathbf{R}^n \mid \|x - x_0\| \leq R\}&amp;lt;/math&amp;gt;なる &amp;lt;math&amp;gt;x_0 \in \mathbf{R}^n&amp;lt;/math&amp;gt; 及び &amp;lt;math&amp;gt;R \in \mathbf{R}&amp;lt;/math&amp;gt; が既知である.&lt;br /&gt;
&lt;br /&gt;
:&amp;lt;math&amp;gt;\mbox{(ii)} \quad P\, &amp;lt;/math&amp;gt; が空でないならば, その体積は &amp;lt;math&amp;gt;\varepsilon (&amp;gt; 0)&amp;lt;/math&amp;gt; 以上である.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
　楕円体法では, いわば 「&amp;lt;math&amp;gt;n\,&amp;lt;/math&amp;gt;次元空間での [[楕円体]]を用いた2分探索」を行うことにより, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt;の点を求める. 各反復では, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt;を含む楕円体&amp;lt;math&amp;gt;E_k\,&amp;lt;/math&amp;gt; を常に保持する. 楕円体 &amp;lt;math&amp;gt;E_k\,&amp;lt;/math&amp;gt;は, 中心と呼ばれるベクトル &amp;lt;math&amp;gt;x_k \in \mathbf{R}^n&amp;lt;/math&amp;gt;及び&amp;lt;math&amp;gt;n \times n&amp;lt;/math&amp;gt;正定値対称行列&amp;lt;math&amp;gt;B_k\,&amp;lt;/math&amp;gt;を用いて,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&amp;lt;math&amp;gt;E_k = \{x \in \mathbf{R}^n \mid (x - x_k)^{\top} B_k^{-1}(x - x_k) \leq 1\}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
と表される. なお, 上記の仮定 (i) より, &amp;lt;math&amp;gt;B_0 = R^2 I&amp;lt;/math&amp;gt; と選ぶことができる.&lt;br /&gt;
&lt;br /&gt;
　各反復では, &amp;lt;math&amp;gt;x_k \in P\,&amp;lt;/math&amp;gt; が成り立つか否かを調べる.  成り立つ場合はアルゴリズムを停止する.  成り立たない場合は, &amp;lt;math&amp;gt;a_i^{\top} x &amp;gt; b_i&amp;lt;/math&amp;gt; なる &amp;lt;math&amp;gt;i\,&amp;lt;/math&amp;gt; が存在するので, そのような &amp;lt;math&amp;gt;i\,&amp;lt;/math&amp;gt; を見つける.  中心 &amp;lt;math&amp;gt;x_k\,&amp;lt;/math&amp;gt; を通る超平面 &amp;lt;math&amp;gt;a_i^{\top} x = a_i^{\top} x_k&amp;lt;/math&amp;gt; により,  &amp;lt;math&amp;gt;E_k\,&amp;lt;/math&amp;gt; を半分にして得られる集合 &amp;lt;math&amp;gt;E_k' = \{x \in E_k \mid a_i^{\top} x \leq a_i^{\top} x_k\}&amp;lt;/math&amp;gt;  に対し, これを含み, かつ体積が最小の楕円体を &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; とおく.  そのような楕円体 &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; を表す &amp;lt;math&amp;gt;x_{k+1}\,&amp;lt;/math&amp;gt; 及び &amp;lt;math&amp;gt;B_{k+1}\,&amp;lt;/math&amp;gt; は 次のようにして求められる:&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;x_{k+1} = x_k - \frac{1}{n+1}d,\qquad&lt;br /&gt;
B_{k+1} = \frac{n^2}{n^2 - 1}\left(B_k - \frac{2}{n+1}d d^{\top} \right).&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
ここで, &amp;lt;math&amp;gt;d = (1/\sqrt{a_i^{\top}B_k a_i})B_k a_i&amp;lt;/math&amp;gt;である.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
　上記のように &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; を定めると, 楕円体の体積の比は,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;\frac{{\rm vol}(E_{k+1})}{{\rm vol}(E_{k})}&lt;br /&gt;
 = \left(\left(\frac{n}{n+1}\right)^{n+1}&lt;br /&gt;
 \left(\frac{n}{n-1}\right)^{n-1}\right)^{1/2}&lt;br /&gt;
 &amp;lt; e^{-1/(2n)} &amp;lt; 1&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
となる.  したがって, 楕円体の体積は線形空間の次元 &amp;lt;math&amp;gt;n\,&amp;lt;/math&amp;gt; のみに依存する一定の比率で 減少していく.  反復回数 &amp;lt;math&amp;gt;k\,&amp;lt;/math&amp;gt; が十分大きくなり 楕円体の体積が &amp;lt;math&amp;gt;\varepsilon&amp;lt;/math&amp;gt; 未満になった とすると, 仮定 (ii) により &amp;lt;math&amp;gt;P = \emptyset&amp;lt;/math&amp;gt; であることが分かるので, アルゴリズムを停止する.  &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; が多面体であることから, 仮定 (i), (ii) での &amp;lt;math&amp;gt;x_0\,&amp;lt;/math&amp;gt;, &amp;lt;math&amp;gt;R\,&amp;lt;/math&amp;gt;, &amp;lt;math&amp;gt;\varepsilon&amp;lt;/math&amp;gt; をうまく選ぶことができ, それにより多項式回の反復で終了することが導かれる.&lt;br /&gt;
&lt;br /&gt;
　上記で述べた楕円体法の改良として,「深い」(deep)カットを用いた方法が 提案されている.  この方法では, 半楕円体 &amp;lt;math&amp;gt;E_k'\,&amp;lt;/math&amp;gt; の代わりに, より小さな集合 &amp;lt;math&amp;gt;\{x \in E_k \mid a_i^{\top} x \leq b_i\}&amp;lt;/math&amp;gt; に注目し, それを含む最小の楕円体 &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; へと更新していく.  この場合, 連続する2つの楕円体の比は, &amp;lt;math&amp;gt;\alpha = (a_i^{\top} x_k - b_i)/\sqrt{a_i^{\top}B_k a_i}&amp;lt;/math&amp;gt; とおくと,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;\frac{{\rm vol}(E_{k+1})}{{\rm vol}(E_{k})}&lt;br /&gt;
 = \left(\left(\frac{n(1 - \alpha)}{n+1}\right)^{n+1}&lt;br /&gt;
 \left(\frac{n (1 + \alpha)}{n-1}\right)^{n-1}\right)^{1/2}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
と表され, また &amp;lt;math&amp;gt;0 \leq \alpha \leq 1&amp;lt;/math&amp;gt; なので, 元の楕円体法に比べて反復回数が減少することがわかる.&lt;br /&gt;
&lt;br /&gt;
　さらに, 代用(surrogate)カット, 平行(parallel)カットを用いれば, 楕円体の体積比をより減少させることができる. また, 深いカットと全く逆の考え方により得られる「浅い」(shallow)カットを用いると, 連続する楕円体の体積比は逆に増加するが, 理論的には, 他のカットでは得られない興味深い結果を得ることができる. これらのカットについての詳細は参考文献を参照されたい.&lt;br /&gt;
&lt;br /&gt;
　線形計画法に対する楕円体法の適用は, 線形計画のみならず組合せ最適化においても重要な結果を残している. その1つが, 最適化問題と[[分離問題]]の多項式時間に関する等価性であろう.&lt;br /&gt;
&lt;br /&gt;
　&amp;lt;math&amp;gt;P \subseteq \mathbf{R}^n&amp;lt;/math&amp;gt; を多面体とする.  ただし, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; の不等式系による表現は必ずしも与えられていない, と仮定する.  分離問題とは, 任意のベクトル &amp;lt;math&amp;gt;x_* \in \mathbf{R}^n&amp;lt;/math&amp;gt; が与えられたとき,  &amp;lt;math&amp;gt;x_* \in P&amp;lt;/math&amp;gt; か否かを判定し, &amp;lt;math&amp;gt;x_* \not\in P&amp;lt;/math&amp;gt; ならば &amp;lt;math&amp;gt;a^{\top} x \leq b&amp;lt;/math&amp;gt; &amp;lt;math&amp;gt;(\forall x \in P)&amp;lt;/math&amp;gt; かつ &amp;lt;math&amp;gt;a^{\top} x_* &amp;gt; b&amp;lt;/math&amp;gt; なる &amp;lt;math&amp;gt;a \in \mathbf{R}^n&amp;lt;/math&amp;gt; 及び &amp;lt;math&amp;gt;b \in \mathbf{R}&amp;lt;/math&amp;gt; を求める問題のことである.  上記で述べた楕円体法のアルゴリズムを見ると分かるが, 必ずしも &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; の不等式表現が陽に与えられている必要はなく, 各反復において &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; と &amp;lt;math&amp;gt;x_k\,&amp;lt;/math&amp;gt; に関する分離問題を多項式時間で解けるならば, 最適化問題もまた多項式時間で解くことができる.  一方, この逆もまた成り立つことが Gr&amp;amp;ouml;tschel-Lov&amp;amp;aacute;sz-Schrijver (1981) 及び Padberg-Rao (1981) によって証明されている.&lt;br /&gt;
&lt;br /&gt;
　この他にも, 劣モジュラ関数の最小化やパーフェクトグラフでの最大 安定集合問題などの組合せ最適化問題が多項式時間で解ける, という結果が, 楕円体法の適用により導かれている.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
----&lt;br /&gt;
'''参考文献'''&lt;br /&gt;
&lt;br /&gt;
[1] R. G. Bland, D. Goldfarb and M. J. Todd, &amp;quot;The ellipsoid method: a survey,&amp;quot; ''Operations Research'', '''29'''(1981) 1039-1091.&lt;br /&gt;
&lt;br /&gt;
[2] M. Gr&amp;amp;ouml;tschel, L. Lov&amp;amp;aacute;sz and A. Schrijver, ''Geometric Algorithms and Combinatorial Optimization'', Springer-Verlag, 1991.&lt;br /&gt;
&lt;br /&gt;
[3] 伊理正夫,『線形計画法』, 共立出版, 1986.&lt;br /&gt;
&lt;br /&gt;
[4] G. L. Nemhauser, A. H. G. Rinnooy Kan and M. J. Todd, eds., ''Optimization'', Handbooks in Operations Research and Management Science Vol. 1, North-Holland, 1989. 伊理正夫, 今野浩, 刀根薫 監訳,『最適化ハンドブック』,朝倉書店, 1995.&lt;br /&gt;
&lt;br /&gt;
[[Category:線形計画|だえんたいほう]]&lt;/div&gt;</summary>
		<author><name>Tamura</name></author>
	</entry>
	<entry>
		<id>https://orsj-ml.org/orwiki/wiki/index.php?title=%E6%A5%95%E5%86%86%E4%BD%93%E6%B3%95&amp;diff=9988</id>
		<title>楕円体法</title>
		<link rel="alternate" type="text/html" href="https://orsj-ml.org/orwiki/wiki/index.php?title=%E6%A5%95%E5%86%86%E4%BD%93%E6%B3%95&amp;diff=9988"/>
		<updated>2008-08-22T07:08:05Z</updated>

		<summary type="html">&lt;p&gt;Tamura: 詳説の第１段落を概要に移した．旧概要は最後の一文以外削除．（ほぼ同様の文章を繰返していたため）&lt;/p&gt;
&lt;hr /&gt;
&lt;div&gt;'''【だえんたいほう (ellipsoid method)】'''&lt;br /&gt;
&lt;br /&gt;
=== 概要 ===&lt;br /&gt;
　[[楕円体法]](ellipsoid method) は, 微分不可能な凸計画問題に対する解法として, 1976年にユーディンとネミロフスキー(Yudin-Nemirovskii) によって提案された.  その後1979年にカチヤン(L. G. Khachiyan)は, 楕円体法が線形計画問題に対する[[多項式時間解法]] に成り得ることを示した.  この結果は, 長年にわたって未解決のままであった,「線形計画問題を多項式時間で解くことが出来るか?」という計算複雑度に関する問題への肯定的な答を与えるものであり, 数理計画の分野に衝撃をもたらした.  線形計画法に対する解法としての楕円体法は, 実用の面からは単体法などに比べて効率が悪く, また理論的な計算量の面からも, 1984年以降に現れた内点法に劣っていることが分かっている. しかしながら, 線形, 非線形, 組合せなどの最適化問題に対し, 様々な理論的な結果を残している. 「&amp;lt;math&amp;gt;n&amp;lt;/math&amp;gt; 次元空間における楕円体を用いた2分探索」のような解法である.&lt;br /&gt;
&lt;br /&gt;
=== 詳説 ===&lt;br /&gt;
　楕円体法を説明するために, 与えられた多面体 &amp;lt;math&amp;gt;P = \{x \in \mathbf{R}^n \mid a_i^{\top} x \leq b_i\ (i = 1, 2, \cdots, m)\}&amp;lt;/math&amp;gt; が空であるか否かを判定し, また空でなければ &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; に含まれる点を求める, という非空性問題について考える. なお, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; の上での線形関数 &amp;lt;math&amp;gt;c^{\top} x&amp;lt;/math&amp;gt; の最大化問題を解くには, &amp;lt;math&amp;gt;\{x \in P \mid c^{\top} x \leq \gamma\}&amp;lt;/math&amp;gt; という多面体に対する非空性問題を繰り返し解き, &amp;lt;math&amp;gt;\gamma \in \mathbf{R}&amp;lt;/math&amp;gt; に関する2分探索を行えば良い.&lt;br /&gt;
&lt;br /&gt;
　この問題に対し, 以下の2点を仮定する.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
:&amp;lt;math&amp;gt;\mbox{(i)} \quad  P \subseteq \{x \in \mathbf{R}^n \mid \|x - x_0\| \leq R\}&amp;lt;/math&amp;gt;なる &amp;lt;math&amp;gt;x_0 \in \mathbf{R}^n&amp;lt;/math&amp;gt; 及び &amp;lt;math&amp;gt;R \in \mathbf{R}&amp;lt;/math&amp;gt; が既知である.&lt;br /&gt;
&lt;br /&gt;
:&amp;lt;math&amp;gt;\mbox{(ii)} \quad P\, &amp;lt;/math&amp;gt; が空でないならば, その体積は &amp;lt;math&amp;gt;\varepsilon (&amp;gt; 0)&amp;lt;/math&amp;gt; 以上である.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
　楕円体法では, いわば 「&amp;lt;math&amp;gt;n\,&amp;lt;/math&amp;gt;次元空間での [[楕円体]]を用いた2分探索」を行うことにより, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt;の点を求める. 各反復では, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt;を含む楕円体&amp;lt;math&amp;gt;E_k\,&amp;lt;/math&amp;gt; を常に保持する. 楕円体 &amp;lt;math&amp;gt;E_k\,&amp;lt;/math&amp;gt;は, 中心と呼ばれるベクトル &amp;lt;math&amp;gt;x_k \in \mathbf{R}^n&amp;lt;/math&amp;gt;及び&amp;lt;math&amp;gt;n \times n&amp;lt;/math&amp;gt;正定値対称行列&amp;lt;math&amp;gt;B_k\,&amp;lt;/math&amp;gt;を用いて,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&amp;lt;math&amp;gt;E_k = \{x \in \mathbf{R}^n \mid (x - x_k)^{\top} B_k^{-1}(x - x_k) \leq 1\}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
と表される. なお, 上記の仮定 (i) より, &amp;lt;math&amp;gt;B_0 = R^2 I&amp;lt;/math&amp;gt; と選ぶことができる.&lt;br /&gt;
&lt;br /&gt;
　各反復では, &amp;lt;math&amp;gt;x_k \in P\,&amp;lt;/math&amp;gt; が成り立つか否かを調べる.  成り立つ場合はアルゴリズムを停止する.  成り立たない場合は, &amp;lt;math&amp;gt;a_i^{\top} x &amp;gt; b_i&amp;lt;/math&amp;gt; なる &amp;lt;math&amp;gt;i\,&amp;lt;/math&amp;gt; が存在するので, そのような &amp;lt;math&amp;gt;i\,&amp;lt;/math&amp;gt; を見つける.  中心 &amp;lt;math&amp;gt;x_k\,&amp;lt;/math&amp;gt; を通る超平面 &amp;lt;math&amp;gt;a_i^{\top} x = a_i^{\top} x_k&amp;lt;/math&amp;gt; により,  &amp;lt;math&amp;gt;E_k\,&amp;lt;/math&amp;gt; を半分にして得られる集合 &amp;lt;math&amp;gt;E_k' = \{x \in E_k \mid a_i^{\top} x \leq a_i^{\top} x_k\}&amp;lt;/math&amp;gt;  に対し, これを含み, かつ体積が最小の楕円体を &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; とおく.  そのような楕円体 &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; を表す &amp;lt;math&amp;gt;x_{k+1}\,&amp;lt;/math&amp;gt; 及び &amp;lt;math&amp;gt;B_{k+1}\,&amp;lt;/math&amp;gt; は 次のようにして求められる:&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;x_{k+1} = x_k - \frac{1}{n+1}d,\qquad&lt;br /&gt;
B_{k+1} = \frac{n^2}{n^2 - 1}\left(B_k - \frac{2}{n+1}d d^{\top} \right).&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
ここで, &amp;lt;math&amp;gt;d = (1/\sqrt{a_i^{\top}B_k a_i})B_k a_i&amp;lt;/math&amp;gt;である.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
　上記のように &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; を定めると, 楕円体の体積の比は,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;\frac{{\rm vol}(E_{k+1})}{{\rm vol}(E_{k})}&lt;br /&gt;
 = \left(\left(\frac{n}{n+1}\right)^{n+1}&lt;br /&gt;
 \left(\frac{n}{n-1}\right)^{n-1}\right)^{1/2}&lt;br /&gt;
 &amp;lt; e^{-1/(2n)} &amp;lt; 1&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
となる.  したがって, 楕円体の体積は線形空間の次元 &amp;lt;math&amp;gt;n\,&amp;lt;/math&amp;gt; のみに依存する一定の比率で 減少していく.  反復回数 &amp;lt;math&amp;gt;k\,&amp;lt;/math&amp;gt; が十分大きくなり 楕円体の体積が &amp;lt;math&amp;gt;\varepsilon&amp;lt;/math&amp;gt; 未満になった とすると, 仮定 (ii) により &amp;lt;math&amp;gt;P = \emptyset&amp;lt;/math&amp;gt; であることが分かるので, アルゴリズムを停止する.  &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; が多面体であることから, 仮定 (i), (ii) での &amp;lt;math&amp;gt;x_0\,&amp;lt;/math&amp;gt;, &amp;lt;math&amp;gt;R\,&amp;lt;/math&amp;gt;, &amp;lt;math&amp;gt;\varepsilon&amp;lt;/math&amp;gt; をうまく選ぶことができ, それにより多項式回の反復で終了することが導かれる.&lt;br /&gt;
&lt;br /&gt;
　上記で述べた楕円体法の改良として,「深い」(deep)カットを用いた方法が 提案されている.  この方法では, 半楕円体 &amp;lt;math&amp;gt;E_k'\,&amp;lt;/math&amp;gt; の代わりに, より小さな集合 &amp;lt;math&amp;gt;\{x \in E_k \mid a_i^{\top} x \leq b_i\}&amp;lt;/math&amp;gt; に注目し, それを含む最小の楕円体 &amp;lt;math&amp;gt;E_{k+1}\,&amp;lt;/math&amp;gt; へと更新していく.  この場合, 連続する2つの楕円体の比は, &amp;lt;math&amp;gt;\alpha = (a_i^{\top} x_k - b_i)/\sqrt{a_i^{\top}B_k a_i}&amp;lt;/math&amp;gt; とおくと,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;\frac{{\rm vol}(E_{k+1})}{{\rm vol}(E_{k})}&lt;br /&gt;
 = \left(\left(\frac{n(1 - \alpha)}{n+1}\right)^{n+1}&lt;br /&gt;
 \left(\frac{n (1 + \alpha)}{n-1}\right)^{n-1}\right)^{1/2}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
と表され, また &amp;lt;math&amp;gt;0 \leq \alpha \leq 1&amp;lt;/math&amp;gt; なので, 元の楕円体法に比べて反復回数が減少することがわかる.&lt;br /&gt;
&lt;br /&gt;
　さらに, 代用(surrogate)カット, 平行(parallel)カットを用いれば, 楕円体の体積比をより減少させることができる. また, 深いカットと全く逆の考え方により得られる「浅い」(shallow)カットを用いると, 連続する楕円体の体積比は逆に増加するが, 理論的には, 他のカットでは得られない興味深い結果を得ることができる. これらのカットについての詳細は参考文献を参照されたい.&lt;br /&gt;
&lt;br /&gt;
　線形計画法に対する楕円体法の適用は, 線形計画のみならず組合せ最適化においても重要な結果を残している. その1つが, 最適化問題と[[分離問題]]の多項式時間に関する等価性であろう.&lt;br /&gt;
&lt;br /&gt;
　&amp;lt;math&amp;gt;P \subseteq \mathbf{R}^n&amp;lt;/math&amp;gt; を多面体とする.  ただし, &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; の不等式系による表現は必ずしも与えられていない, と仮定する.  分離問題とは, 任意のベクトル &amp;lt;math&amp;gt;x_* \in \mathbf{R}^n&amp;lt;/math&amp;gt; が与えられたとき,  &amp;lt;math&amp;gt;x_* \in P&amp;lt;/math&amp;gt; か否かを判定し, &amp;lt;math&amp;gt;x_* \not\in P&amp;lt;/math&amp;gt; ならば &amp;lt;math&amp;gt;a^{\top} x \leq b&amp;lt;/math&amp;gt; &amp;lt;math&amp;gt;(\forall x \in P)&amp;lt;/math&amp;gt; かつ &amp;lt;math&amp;gt;a^{\top} x_* &amp;gt; b&amp;lt;/math&amp;gt; なる &amp;lt;math&amp;gt;a \in \mathbf{R}^n&amp;lt;/math&amp;gt; 及び &amp;lt;math&amp;gt;b \in \mathbf{R}&amp;lt;/math&amp;gt; を求める問題のことである.  上記で述べた楕円体法のアルゴリズムを見ると分かるが, 必ずしも &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; の不等式表現が陽に与えられている必要はなく, 各反復において &amp;lt;math&amp;gt;P\,&amp;lt;/math&amp;gt; と &amp;lt;math&amp;gt;x_k\,&amp;lt;/math&amp;gt; に関する分離問題を多項式時間で解けるならば, 最適化問題もまた多項式時間で解くことができる.  一方, この逆もまた成り立つことが Gr&amp;amp;ouml;tschel-Lov&amp;amp;aacute;sz-Schrijver (1981) 及び Padberg-Rao (1981) によって証明されている.&lt;br /&gt;
&lt;br /&gt;
　この他にも, 劣モジュラ関数の最小化やパーフェクトグラフでの最大 安定集合問題などの組合せ最適化問題が多項式時間で解ける, という結果が, 楕円体法の適用により導かれている.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
----&lt;br /&gt;
'''参考文献'''&lt;br /&gt;
&lt;br /&gt;
[1] R. G. Bland, D. Goldfarb and M. J. Todd, &amp;quot;The ellipsoid method: a survey,&amp;quot; ''Operations Research'', '''29'''(1981) 1039-1091.&lt;br /&gt;
&lt;br /&gt;
[2] M. Gr&amp;amp;ouml;tschel, L. Lov&amp;amp;aacute;sz and A. Schrijver, ''Geometric Algorithms and Combinatorial Optimization'', Springer-Verlag, 1991.&lt;br /&gt;
&lt;br /&gt;
[3] 伊理正夫,『線形計画法』, 共立出版, 1986.&lt;br /&gt;
&lt;br /&gt;
[4] G. L. Nemhauser, A. H. G. Rinnooy Kan and M. J. Todd, eds., ''Optimization'', Handbooks in Operations Research and Management Science Vol. 1, North-Holland, 1989. 伊理正夫, 今野浩, 刀根薫 監訳,『最適化ハンドブック』,朝倉書店, 1995.&lt;br /&gt;
&lt;br /&gt;
[[Category:線形計画|だえんたいほう]]&lt;/div&gt;</summary>
		<author><name>Tamura</name></author>
	</entry>
	<entry>
		<id>https://orsj-ml.org/orwiki/wiki/index.php?title=%E5%8D%98%E4%BD%93%E6%B3%95&amp;diff=9987</id>
		<title>単体法</title>
		<link rel="alternate" type="text/html" href="https://orsj-ml.org/orwiki/wiki/index.php?title=%E5%8D%98%E4%BD%93%E6%B3%95&amp;diff=9987"/>
		<updated>2008-08-22T07:02:13Z</updated>

		<summary type="html">&lt;p&gt;Tamura: 線形計画問題（線形計画にはる）にリンクをはり，文章の「てにをは」を修正&lt;/p&gt;
&lt;hr /&gt;
&lt;div&gt;'''【たんたいほう (simplex method)】'''&lt;br /&gt;
&lt;br /&gt;
=== 概要 ===&lt;br /&gt;
1947年, ダンツィク(G.B. Dantzig) によって提案された,[[線形計画|線形計画問題]]を解くための手法.理論的には有限回反復での収束性しか示されていないが,実用的には高速に最適解が得られることが知られている.また, 実用の際には, 単体法を改良した改訂単体法が良く用いられる.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
=== 詳説 ===&lt;br /&gt;
　ダンツイック（G. B. Dantzig）によって提案された[[単体法]]（simplex method）（[[シンプレックス法]]）は, [[カーマーカー法]]の出現する1984年までは線形計画問題を解くもっとも有効な解法とされていた.  単体法は問題の構造を巧妙に利用しており, いまでも, 線形計画, 整数計画, 組合せ最適化の多くの分野で使われている[2, 3, 4].  以下, 次のような基準形と呼ばれる線形計画問題を考えることにする.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;\mbox{(P)} \quad&lt;br /&gt;
\begin{array}{lll}&lt;br /&gt;
          &amp;amp; \mbox{max.}  &amp;amp; {\displaystyle \sum_{j=1}^{n}c_j x_j} \\&lt;br /&gt;
&amp;amp; \mbox{s. t.}  &amp;amp; {\displaystyle \sum_{j=1}^{n}a_{ij} x_j}&lt;br /&gt;
                           \leq b_i \; \;  (i=1,2,\ldots,m), &lt;br /&gt;
              x_1,x_2,\ldots ,x_n \geq 0.&lt;br /&gt;
\end{array}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
ここで, &amp;lt;math&amp;gt;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)&amp;lt;/math&amp;gt;は実数, &amp;lt;math&amp;gt;\boldsymbol{x}=(x_1,x_2,\ldots,x_m)&amp;lt;/math&amp;gt;は&amp;lt;math&amp;gt;n\,&amp;lt;/math&amp;gt;個の変数からなる&amp;lt;math&amp;gt;n\,&amp;lt;/math&amp;gt;次元ベクトルである.  すべての線形計問題は, 簡単な変換によって基準形に帰着できる.  問題（P）の制約条件を満たす&amp;lt;math&amp;gt;\boldsymbol{x}&amp;lt;/math&amp;gt;を[[実行可能解]]（feasible solution）, その集まりを[[実行可能多面体]]（feasible polyhedron）と呼ぶ.  実行可能解のなかで目的関数を最大にするものを[[最適解]]（optimal solution）と呼ぶ.&lt;br /&gt;
&lt;br /&gt;
　ここで, 問題（P）にスラック変数と呼ばれる非負の変数&amp;lt;math&amp;gt;x_{n+i} \, (i=1,\ldots,m)&amp;lt;/math&amp;gt;を導入する&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;\begin{array}{ll}&lt;br /&gt;
\mbox{max.}        &amp;amp; c_1 x_1+\cdots+ c_nx_n   \\&lt;br /&gt;
\mbox{s. t.} &amp;amp; a_{i1} x_1 +\cdots+a_{in} x_n + x_{n+i} = b_i&lt;br /&gt;
             \; \; (i=1, 2, \ldots, m), \\&lt;br /&gt;
            &amp;amp; x_1 \geq 0,\ldots,x_{n+m}\geq 0.&lt;br /&gt;
\end{array}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
さらに, 目的関数を &amp;lt;math&amp;gt;z=c_0+\sum_{j=1}^n c_j x_j \ (c_0=0)&amp;lt;/math&amp;gt; と書き, 制約式左辺のスラック変数以外の項を移行すると, 上の問題は以下のようになる.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&lt;br /&gt;
&amp;lt;math&amp;gt;\begin{array}{ll}&lt;br /&gt;
\mbox{max.}        &amp;amp; z  \\&lt;br /&gt;
\mbox{s. t.} &amp;amp; x_{n+i}= b_i-a_{i1}x_{1}-\cdots-a_{in}x_{n}&lt;br /&gt;
               \quad (i=1, 2, \ldots, m), \\&lt;br /&gt;
            &amp;amp; z=c_0+c_1 x_1+\cdots+c_n x_n, \\&lt;br /&gt;
            &amp;amp; x_1 \geq 0,\ldots, x_{n+m} \geq 0.&lt;br /&gt;
\end{array}&amp;lt;/math&amp;gt;&lt;br /&gt;
&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
制約式の左辺の変数を基底変数, 右辺の変数を非基底変数と呼ぶ.  ここで, 非基底変数の添字を&amp;lt;math&amp;gt;N(1)=1,\ldots,N(n)=n&amp;lt;/math&amp;gt;, 基底変数の添字を&amp;lt;math&amp;gt;B(1)=n+1,\ldots,B(m)=n+m&amp;lt;/math&amp;gt;とおき, 各変数の係数の添字も対応するものにする.　非基底変数の値を任意に固定すると基底変数の値は一意に定まるが,　特に非基底変数をすべて&amp;lt;math&amp;gt;0\,&amp;lt;/math&amp;gt;に固定して得られる解を[[基底解]]（basic solution）と呼ぶ.  また, 上の等式系を辞書と呼ぶ.  係数&amp;lt;math&amp;gt;b_i \, (i=1,\ldots,m)&amp;lt;/math&amp;gt;が非負のとき, 基底解は実行可能解であるが, このような辞書を実行可能辞書と呼ぶ.  線形計画問題の実行可能な基底解は実行可能多面体の頂点に対応する [3].  単体法は目的関数の値を最大化する実行可能な基底解を求めるものであり, 辞書を等価な辞書へと繰り返し変換することによって実現される.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
単体法&lt;br /&gt;
&lt;br /&gt;
入力：実行可能な辞書&lt;br /&gt;
&lt;br /&gt;
ステップ１（最適性判定）&lt;br /&gt;
&lt;br /&gt;
:(1) &amp;lt;math&amp;gt;c_{N(s)}&amp;gt;0\,&amp;lt;/math&amp;gt;となる添字&amp;lt;math&amp;gt;s (1\le s \le n)&amp;lt;/math&amp;gt;を1つ選ぶ.&lt;br /&gt;
&lt;br /&gt;
:(2) もし, そのような添字がなければ, 現在の基底解は最適となり, 終了する.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
ステップ２（有界性判定）&lt;br /&gt;
&lt;br /&gt;
:(1) &amp;lt;math&amp;gt;\textstyle \frac{b_r}{a_{rN(s)}}=\min\{\frac{b_i}{a_{iN(s)}} : a_{iN(s)}&amp;gt;0, 1 \le i \le m \}&amp;lt;/math&amp;gt;を計算し, 行番号&amp;lt;math&amp;gt;r\,&amp;lt;/math&amp;gt;を求める（この操作は現在の基底解から実行可能性を保ったまま非基底変数&amp;lt;math&amp;gt;x_{N(s)}\,&amp;lt;/math&amp;gt;だけをできる限り増加させる）.&lt;br /&gt;
&lt;br /&gt;
:(2) そのような番号&amp;lt;math&amp;gt;r\,&amp;lt;/math&amp;gt;がなければ（つまり, 変数&amp;lt;math&amp;gt;x_{N(s)}\,&amp;lt;/math&amp;gt;を限りなく増加させることができる）, 問題は非有界となり, 終了する.&lt;br /&gt;
&lt;br /&gt;
:(3) &amp;lt;math&amp;gt;r\,&amp;lt;/math&amp;gt;が存在する場合, ステップ３へ.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
ステップ３（（&amp;lt;math&amp;gt;r, s\,&amp;lt;/math&amp;gt;）を中心とするピボット演算を行う）&lt;br /&gt;
&lt;br /&gt;
:（1）辞書の&amp;lt;math&amp;gt;r\,&amp;lt;/math&amp;gt;番目の式について&amp;lt;math&amp;gt;x_{N(s)}\,&amp;lt;/math&amp;gt;の表現式を求める.&lt;br /&gt;
&lt;br /&gt;
:（2）求めた&amp;lt;math&amp;gt;x_{N(s)}\,&amp;lt;/math&amp;gt;の式を辞書の他の行 (&amp;lt;math&amp;gt;i\not = r\,&amp;lt;/math&amp;gt;) に代入する.&lt;br /&gt;
&lt;br /&gt;
:（3）添字&amp;lt;math&amp;gt;N(s)\,&amp;lt;/math&amp;gt;と&amp;lt;math&amp;gt;B(r)\,&amp;lt;/math&amp;gt;を交換し, ステップ１へ.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
実行可能解をもつ問題に対して単体法が有限回で終了すれば, その問題は最適解をもつか, または非有界である.  しかし, 辞書の右辺に&amp;lt;math&amp;gt;0\,&amp;lt;/math&amp;gt;の定数&amp;lt;math&amp;gt;b_i\,&amp;lt;/math&amp;gt;を少なくとも1つ含むとき, 辞書は退化しており, 単体法は有限回で終わるとは限らない.  このとき, ピボット演算を施しても, 基底解が変化しないことになり, 目的関数値も変化しない.  有限回収束を保証するための手段として, Blandのピボット選択規則 [1] が提案されており, これを用いると単体法は必ず有限回で収束する.&lt;br /&gt;
&lt;br /&gt;
　先の説明で述べたように, 単体法の入力として実行可能な辞書が必要とされる.  現在の辞書が実行可能でない場合, 単体法を用いて実行可能な辞書を求めることができる.  そのために, 新しい変数（人工変数）&amp;lt;math&amp;gt;x_a\,&amp;lt;/math&amp;gt;を用いて補助問題&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&amp;lt;math&amp;gt;\begin{array}{ll}&lt;br /&gt;
\mbox{max.}        &amp;amp; -x_a  \\&lt;br /&gt;
\mbox{s. t.} &amp;amp; x_{n+i}= b_i-a_{i1}x_{1}-\cdots-a_{in}x_{n}+x_a&lt;br /&gt;
               \; \;  (i=1, 2, \ldots, m),  \\&lt;br /&gt;
            &amp;amp; x_1 \geq 0,\ldots, x_{n+m} \geq 0, \; x_a \geq 0, &lt;br /&gt;
\end{array}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
を作る.  元問題が実行可能であることと補助問題の最適目的関数値が&amp;lt;math&amp;gt;0\,&amp;lt;/math&amp;gt;であることは等価なので, 補助問題を解くことにより, 元問題の実行可能解の有無を判定できる.&lt;br /&gt;
&lt;br /&gt;
　一般に, 基準形の線形計画問題を解く際には, 2段階単体法と呼ばれるものを用いる.  第1段階では, 補助問題を解き, 実行可能解の有無を判定し, 実行可能解が存在する場合は実行可能辞書を求める.  第2段階では, 求めた実行可能辞書を初期の辞書として, もう一度単体法を使い, 元問題を解く.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
----&lt;br /&gt;
'''参考文献'''&lt;br /&gt;
&lt;br /&gt;
[1]  R. G. Bland, &amp;quot;New finite pivot ruless for the simplex method&amp;quot;, ''Mathematics of Operations Research'' '''2'''（1977）, 103-107.&amp;lt;BR&amp;gt;&lt;br /&gt;
&lt;br /&gt;
[2]  V. Chv&amp;amp;aacute;tal, ''Linear Programming'', W. H. Freeman and Company, 1983.&amp;lt;BR&amp;gt;&lt;br /&gt;
&lt;br /&gt;
[3]  G. B. Dantzig, ''Linear Programming and Extensions'', Princeton Univesity Press, 1963.&amp;lt;BR&amp;gt;&lt;br /&gt;
&lt;br /&gt;
[4]  G. L. Nemhauser and L. A. Wolsey, ''Integer and Combinatorial Optimization'', John Wiley &amp;amp; Sons, 1988.&amp;lt;BR&amp;gt;&lt;br /&gt;
&lt;br /&gt;
[[Category:線形計画|たんたいほう]]&lt;/div&gt;</summary>
		<author><name>Tamura</name></author>
	</entry>
	<entry>
		<id>https://orsj-ml.org/orwiki/wiki/index.php?title=%E6%9C%80%E9%81%A9%E5%8C%96%E5%95%8F%E9%A1%8C&amp;diff=9986</id>
		<title>最適化問題</title>
		<link rel="alternate" type="text/html" href="https://orsj-ml.org/orwiki/wiki/index.php?title=%E6%9C%80%E9%81%A9%E5%8C%96%E5%95%8F%E9%A1%8C&amp;diff=9986"/>
		<updated>2008-08-22T06:55:28Z</updated>

		<summary type="html">&lt;p&gt;Tamura: 最小費用フロー問題，最大フロー問題（ネットワークフロー問題にはる），巡回セールスマン問題にリンクをはった&lt;/p&gt;
&lt;hr /&gt;
&lt;div&gt;'''【さいてきかもんだい (optimization problem)】'''&lt;br /&gt;
&lt;br /&gt;
=== 概要 ===&lt;br /&gt;
「与えられた制約条件の下で目的を最適に達成するための数理モデル」で数理計画問題(mathematical programming problem)ともいう. 数学的には,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&amp;lt;math&amp;gt;\mbox{max.} \, &amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&amp;lt;math&amp;gt;f(x) ( \,&amp;lt;/math&amp;gt;あるいは, &amp;lt;math&amp;gt;\mbox{min.} \quad f(x) ) \,&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&amp;lt;math&amp;gt;\mbox{s.t.}\, &amp;lt;/math&amp;gt; &amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&amp;lt;math&amp;gt;x = (x_1,x_2,\ldots,x_n) \in F,&lt;br /&gt;
\,&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
と表現される. ここで, &amp;lt;math&amp;gt;F \,&amp;lt;/math&amp;gt; は &amp;lt;math&amp;gt;n \,&amp;lt;/math&amp;gt; 次元ベクトル空間 &amp;lt;math&amp;gt;\mathbf{R}^n \,&amp;lt;/math&amp;gt; の部分集合(実行可能集合)で,  &amp;lt;math&amp;gt;f \,&amp;lt;/math&amp;gt; は &amp;lt;math&amp;gt;\mathbf{R}^n \,&amp;lt;/math&amp;gt; で定義された実数値関数(目的関数).&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
=== 詳説 ===&lt;br /&gt;
　[[最適化問題]] （optimization problem）（[[数理計画問題]]）は, 「与えられた制約条件の下でより良い目的を達成するための数理モデル」で, 数学的には,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&amp;lt;math&amp;gt;\mbox{max.}  \quad  f(x)   \mbox{(}&amp;lt;/math&amp;gt;あるいは, &amp;lt;math&amp;gt;\mbox{min.}\quad f(x))&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&amp;lt;math&amp;gt;\mbox{s.t.} \quad\quad x = (x_1,x_2,\ldots,x_n) \in F,&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
のように表現される. ここで, &amp;lt;math&amp;gt;F\, &amp;lt;/math&amp;gt; は &amp;lt;math&amp;gt;n\, &amp;lt;/math&amp;gt;次元ベクトル空間 &amp;lt;math&amp;gt;\mathbf {R}^n&amp;lt;/math&amp;gt;の部分集合で[[実行可能集合]]（feasible region）（[[許容集合]]）と呼ばれ,&amp;lt;math&amp;gt;f\, &amp;lt;/math&amp;gt; は &amp;lt;math&amp;gt;\mathbf {R}^n&amp;lt;/math&amp;gt; で定義された実数値関数で[[目的関数]]（objective function）と呼ばれる. また, &amp;lt;math&amp;gt;x \in F&amp;lt;/math&amp;gt; を[[実行可能解]]（[[許容解]]）, 実行可能解の中で目的関数の最大（あるいは, 最小）を達成する&amp;lt;math&amp;gt;x\, &amp;lt;/math&amp;gt;を[[最適解]]（optimal solution）と呼ぶ. 実行可能集合 &amp;lt;math&amp;gt;F\, &amp;lt;/math&amp;gt; は,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&amp;lt;table align=&amp;quot;center&amp;quot;&amp;gt;&lt;br /&gt;
&amp;lt;tr&amp;gt;&lt;br /&gt;
&amp;lt;td&amp;gt;&amp;lt;math&amp;gt;F = \{ x \in S : g_i(x) \leq 0 \ (i=1,2,\ldots,m) \}&amp;lt;/math&amp;gt;&amp;lt;/td&amp;gt;&lt;br /&gt;
&amp;lt;/tr&amp;gt;&lt;br /&gt;
&amp;lt;/table&amp;gt;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
のように表現されることが多い. ただし, &amp;lt;math&amp;gt;g_i\, &amp;lt;/math&amp;gt; は &amp;lt;math&amp;gt;\mathbf {R}^n&amp;lt;/math&amp;gt; で定義された実数値関数である. また, &amp;lt;math&amp;gt;S\, &amp;lt;/math&amp;gt; は対象とする最適化問題を記述するのに用いられる基礎となる空間と考えればよい. 基礎となる空間 &amp;lt;math&amp;gt;S\, &amp;lt;/math&amp;gt;, 実行可能領域 &amp;lt;math&amp;gt;F\, &amp;lt;/math&amp;gt; を表現するのに使われる関数 &amp;lt;math&amp;gt;g_i (i=1,2,\ldots,m),&amp;lt;/math&amp;gt; および, 目的関数 &amp;lt;math&amp;gt;f\, &amp;lt;/math&amp;gt; に種々の条件を課した多くのクラスの最適化問題が研究されている. 最適化問題は,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
:* &amp;lt;math&amp;gt;S = \mathbf {R}^n&amp;lt;/math&amp;gt;（より一般的には, &amp;lt;math&amp;gt;S\, &amp;lt;/math&amp;gt; は &amp;lt;math&amp;gt;\mathbf {R}^n&amp;lt;/math&amp;gt; の開集合）&lt;br /&gt;
:* 関数 &amp;lt;math&amp;gt;g_i (i=1,2,\ldots,m)&amp;lt;/math&amp;gt; は連続（多くの場合, 微分可能）&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
が仮定され, 変数ベクトル &amp;lt;math&amp;gt;x\, &amp;lt;/math&amp;gt; が連続的な値をとる[[連続最適化問題]]（continuous optimization problem）と,&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
:* &amp;lt;math&amp;gt;S\, &amp;lt;/math&amp;gt; は有限集合, または, 離散的な集合, たとえば, &amp;lt;math&amp;gt;S = \{ x \in \mathbf {R}^n : x_j = 0&amp;lt;/math&amp;gt; または &amp;lt;math&amp;gt;\mathrm{1}\, &amp;lt;/math&amp;gt;&amp;lt;math&amp;gt;\}\, &amp;lt;/math&amp;gt;（有限集合）, &amp;lt;math&amp;gt;S =\, &amp;lt;/math&amp;gt; あるグラフの部分グラフの集まり（有限集合）, &amp;lt;math&amp;gt;S = \{ x \in \mathbf {R}^n : x_j =&amp;lt;/math&amp;gt; 自然数&amp;lt;math&amp;gt;\}\, &amp;lt;/math&amp;gt;（離散無限集合）.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
が仮定され, 変数ベクトル &amp;lt;math&amp;gt;x\, &amp;lt;/math&amp;gt; が離散的な値をとる[[離散最適化問題]]（discrete optimization problem）に, 大きく分けることができる.  関数 &amp;lt;math&amp;gt;f, \ g_i (i=1,2,\ldots,m)&amp;lt;/math&amp;gt; に連続性（あるいは, 微分可能性）のみを仮定する非線形離散最適化問題のクラスも考えることができるが, そのような問題は非常に難しく, 関数 &amp;lt;math&amp;gt;f, \ g_i (i=1,2,\ldots,m)&amp;lt;/math&amp;gt; が線形（あるいは, 高々２次）関数である場合に限定することが多い. このように限定したとしても, 離散最適化問題には広範な応用がある. 個々の最適化問題は, それが定式化された元の（現実）問題と結びつけた名前で呼ばれることも多い. たとえば, [[最小費用フロー問題]], [[ネットワークフロー問題|最大フロー問題]], [[巡回セールスマン問題]], グラフ分割問題等である.  また, 上記の最適化問題に関する分類は厳密ではなく, 連続最適化と離散最適化の両方の特徴を共有する問題（たとえば, 施設配置問題, 線形混合整数計画問題）や, それらからはみ出た最適化問題も存在する.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
----&lt;br /&gt;
&lt;br /&gt;
'''参考文献'''&lt;br /&gt;
&lt;br /&gt;
[1] 福島雅夫, 『数理計画法入門』, 朝倉書店, 1996.&lt;br /&gt;
&lt;br /&gt;
[[Category:線形計画|さいてきかもんだい]]&lt;/div&gt;</summary>
		<author><name>Tamura</name></author>
	</entry>
	<entry>
		<id>https://orsj-ml.org/orwiki/wiki/index.php?title=%E5%9F%BA%E7%A4%8E%E7%B7%A8%EF%BC%9A%E9%A0%85%E7%9B%AE%E4%B8%80%E8%A6%A7&amp;diff=9985</id>
		<title>基礎編：項目一覧</title>
		<link rel="alternate" type="text/html" href="https://orsj-ml.org/orwiki/wiki/index.php?title=%E5%9F%BA%E7%A4%8E%E7%B7%A8%EF%BC%9A%E9%A0%85%E7%9B%AE%E4%B8%80%E8%A6%A7&amp;diff=9985"/>
		<updated>2008-08-22T06:44:51Z</updated>

		<summary type="html">&lt;p&gt;Tamura: /* 線形計画 */&lt;/p&gt;
&lt;hr /&gt;
&lt;div&gt;===線形計画===&lt;br /&gt;
&lt;br /&gt;
[[最適化問題]]　[[線形計画]]　[[単体法]]　[[楕円体法]]　[[内点法]]　[[半正定値計画]]&lt;br /&gt;
&lt;br /&gt;
=== 非線形計画 ===&lt;br /&gt;
&lt;br /&gt;
[[《非線形計画》]]　[[《最適性条件》]]　[[《双対性理論》]]　[[《制約なし最適化》]]　[[《制約付き最適化》]]　[[《大域的最適化》]]　[[《相補性問題》]]　[[《大規模問題の分解法》]]　[[《凸解析》]]　[[《高速微分法》]]　[[《多項式最適化問題》]]&lt;br /&gt;
&lt;br /&gt;
=== 組合せ最適化 ===&lt;br /&gt;
&lt;br /&gt;
[[《整数計画》]]　[[《組合せ最適化問題》]]　[[《多面体理論》]]　[[《アルゴリズム》]]　[[《データ構造》]]　[[《計算の複雑さ》]]　[[《パーフェクトグラフ》]]　[[《グレブナー基底》]]&lt;br /&gt;
&lt;br /&gt;
===グラフ・ネットワーク===&lt;br /&gt;
&lt;br /&gt;
[[《グラフ・ネットワーク》]]　[[《グラフの連結度》]]　[[《最短路問題》]]　[[《最小木問題》]]　[[《巡回セールスマン問題》]]　[[《ネットワーク・フロー問題》]]　[[《マッチング問題》]]　[[《マトロイド》]]　[[《劣モジュラ最適化》]]　[[《離散凸解析》]]　[[《複雑ネットワーク》]]&lt;br /&gt;
&lt;br /&gt;
===スケジューリング===&lt;br /&gt;
&lt;br /&gt;
[[《スケジューリング理論》]]　[[《スケジューリング問題》]]　[[《スケジューリングアルゴリズム》]]&lt;br /&gt;
&lt;br /&gt;
===計算幾何===&lt;br /&gt;
&lt;br /&gt;
[[《凸多面体》]]　[[《アレンジメント》]]　[[《ボロノイ図》]]　[[《三角形分割》]]　[[《幾何グラフ》]]　[[《バケット法》]]　[[《双対変換》]]　[[《木》]]　[[《ランダマイゼーション》]]　[[《ロバスト化技術》]]　[[《計算幾何学》]]&lt;br /&gt;
&lt;br /&gt;
===動的・確率・多目的計画===&lt;br /&gt;
&lt;br /&gt;
[[《動的計画》]]　[[《両的計画》]]　[[《多段確率決定樹表(ツリーテーブル)》]]　[[《不変埋没原理》]]　[[《多目的計画》]]　[[《最適停止》]]　[[《確率計画》]]　&lt;br /&gt;
&lt;br /&gt;
===近似・知能・感覚的手法===&lt;br /&gt;
&lt;br /&gt;
[[《近似アルゴリズム（ヒューリスティックアルゴリズム）》]]　[[《メタヒューリスティクス》]]　[[《ファジイ理論》]]　[[《ソフトコンピューティング》]]&lt;br /&gt;
[[《ラフ集合》]]　[[《ファジィランダム変数》]]&lt;br /&gt;
[[《ニューラルネットワーク》]]　[[《制約充足問題》]]　[[《人工知能》]]　[[《論理プログラミング》]]　[[《サポート・ベクター・マシン》]]&lt;br /&gt;
&lt;br /&gt;
===ゲーム理論===&lt;br /&gt;
&lt;br /&gt;
[[《ゲーム理論》]]　[[《非協力ゲーム理論》]]　[[《戦略形ゲーム》]]　[[《展開形ゲーム》]]　[[《進化と学習のゲーム理論》]]　[[《協力ゲーム理論》]]　[[《交渉ゲーム》]]　[[《提携形ゲーム》]]　[[《ゲームと実験》]]　[[《ゲーム理論の応用》]]　[[《ゲームの解の計算》]]　[[《生物学における進化ゲーム理論》]]&lt;br /&gt;
&lt;br /&gt;
===確率と確率過程===&lt;br /&gt;
&lt;br /&gt;
[[《確率論》]]　[[《確率過程》]]　[[《マルコフ連鎖》]]　[[《ポアソン過程と出生死滅過程》]]　[[《ランダム・ウォークとブラウン運動》]]　[[《マルコフ決定過程》]]　[[《マルコフ連鎖の数値解法》]]&lt;br /&gt;
&lt;br /&gt;
===統計===&lt;br /&gt;
&lt;br /&gt;
[[《回帰分析》]]　[[《クラスター分析》]]　[[《判別関数》]]　[[《多次元尺度構成法》]]　[[《数量化法》]]　[[《多変量解析》]]&lt;br /&gt;
&lt;br /&gt;
===予測===&lt;br /&gt;
&lt;br /&gt;
[[《予測》]]　[[《指数平滑法》]]　[[《季節調整法》]]　[[《自己回帰和分移動平均(ARIMA)モデル》]]　[[《カルマンフィルター》]]　[[《非集計行動モデル》]]　[[《生態学モデル》]]　[[《バス(Bass)モデル》]]&lt;br /&gt;
[[《複雑系による予測モデル》]]&lt;br /&gt;
&lt;br /&gt;
===シミュレーション===&lt;br /&gt;
&lt;br /&gt;
[[《シミュレーション》]]　[[《離散型シミュレーション》]]　[[《モンテカルロ法》]]　[[《一様乱数》]]　[[《非一様乱数》]]　[[《離散型シミュレーションの統計的側面》]]　[[《シミュレーションソフトウェア》]]　[[《シミュレーションモデルの検証》]]　[[《ペトリネット》]]&lt;br /&gt;
&lt;br /&gt;
===待ち行列===&lt;br /&gt;
&lt;br /&gt;
[[待ち行列モデル|《待ち行列》]]　[[《待ち行列モデルの標準形》]]　[[《待ち行列の各種モデル》]]　[[待ち行列モデル M/M/c]]　[[《待ち行列における関係式》]]　[[待ち行列モデル M/G/1]]　&lt;br /&gt;
[[《待ち行列に対するアルゴリズム的解法》]]&lt;br /&gt;
[[待ち行列のバケーションサーバモデル]]　[[《待ち行列における近似》]]　[[《待ち行列における希少事象の評価》]]　[[《待ち行列の極限モデル（流体近似と拡散近似）》]]&lt;br /&gt;
&lt;br /&gt;
===待ち行列ネットワーク===&lt;br /&gt;
&lt;br /&gt;
[[待ち行列ネットワーク]]　[[ジャクソンネットワーク|《待ち行列ネットワーク(ジャクソン型とその応用)》]]　[[BCMPネットワーク|《待ち行列ネットワーク(BCMP型とその応用)》]]　&lt;br /&gt;
[[《積形式解ネットワークとなるための条件》]]　[[《待ち行列ネットワークの近似解析》]]　[[《待ち行列ネットワークの安定性》]]&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
===待ち行列の応用===&lt;br /&gt;
&lt;br /&gt;
[[待ち行列の通信への応用]]　[[待ち行列のコンピュータへの応用]]　[[待ち行列の生産システムへの応用]]&lt;br /&gt;
&lt;br /&gt;
===信頼性・保全性===&lt;br /&gt;
&lt;br /&gt;
[[《信頼性》]]　[[《寿命分布》]]　[[《保全性》]]　[[《予防保全》]]　[[《システムの安全性》]]　[[《故障データ解析》]]　[[《ベイズ信頼性》]]　[[《システムの信頼性》]]　[[《フォールトトレランス》]]　[[《ソフトウェア信頼性》]]&lt;br /&gt;
&lt;br /&gt;
===探索理論===&lt;br /&gt;
&lt;br /&gt;
[[《探索理論》]]　[[《目標存在分布の推定》]]　[[《センサーの探知論》]]　[[《探索モデルと探索の運動学》]]&lt;br /&gt;
[[《静止目標物の最適探索》]]&lt;br /&gt;
[[《移動目標物の最適探索》]]&lt;br /&gt;
[[《探索ゲーム》]]&lt;br /&gt;
[[《ランデブー探索》]]&lt;br /&gt;
[[《探索理論の応用と実例》]]&lt;br /&gt;
&lt;br /&gt;
===経営・経済性工学===&lt;br /&gt;
&lt;br /&gt;
[[《経営戦略》]]　[[《経営モデル》]]　[[《分権管理》]]　[[《経営意思決定》]]　[[《利益計画》]]　[[《経営分析》]]　[[《間接費管理》]]　[[《経済計算》]]　[[《投資案件の評価》]]　[[《財務管理》]]　[[《企業価値評価》]]&lt;br /&gt;
&lt;br /&gt;
===マーケティング===&lt;br /&gt;
&lt;br /&gt;
[[《マーケティング概説》]]　[[《マーケティングモデル》]]　[[《マーケティング・リサーチ》]]&lt;br /&gt;
&lt;br /&gt;
===生産・在庫・ロジスティクス===&lt;br /&gt;
&lt;br /&gt;
[[《生産管理》]]　[[《JIT生産システム》]]　[[《ラインバランシング》]]　[[《在庫管理》]]　[[《経済発注量モデル(EOQモデル)》]]　[[《動的ロットサイズ決定問題》]]　[[《ロットスケジューリング》]]　[[《ロジスティクス》]]　[[《運搬経路問題(配送計画問題, トラック配送問題, 配送問題, 輸送経路問題)》]]　[[《施設配置問題》]]　[[《ロジスティクスネットワーク設計問題》]]&lt;br /&gt;
[[《APS》]]&lt;br /&gt;
&lt;br /&gt;
===企画・開発・プロジェクト・品質・ヒューマン===&lt;br /&gt;
&lt;br /&gt;
[[《製品企画開発》]]　[[《研究開発》]]　[[《プロジェクト管理》]]　[[《総合的品質管理》]]　[[《QC手法》]]　[[《人的資源管理》]]　[[《勤務スケジューリング》]]&lt;br /&gt;
&lt;br /&gt;
===ファイナンス===&lt;br /&gt;
&lt;br /&gt;
[[《モダンポートフォリオ理論(概論)》]]　[[《資産評価理論》]]　[[《企業財務》]]　[[《資産運用モデル》]]　[[《株価変動モデル》]]　[[《証券市場モデル》]]　[[《金利変動モデルと債券価格》]]　[[《金融派生証券(デリバティブ)(概論)》]]　[[《デリバティブ評価モデル》]]&lt;br /&gt;
[[《行動ファイナンス》]]&lt;br /&gt;
[[《証券化》]]&lt;br /&gt;
[[《倒産確率の推計》]]&lt;br /&gt;
[[《リアルオプション》]]&lt;br /&gt;
[[《CAPM》]]&lt;br /&gt;
[[《無裁定価格理論》]]&lt;br /&gt;
&lt;br /&gt;
===公共システム===&lt;br /&gt;
&lt;br /&gt;
[[《選挙制度》]]　[[《議員定数配分問題》]]　[[《投票理論》]]　[[《公共政策OR-I》]]　[[《公共政策OR-II》]]　[[《産業連関分析》]]　[[《エネルギー・環境政策》]]　[[《軍事モデル》]]&lt;br /&gt;
&lt;br /&gt;
===都市システム===&lt;br /&gt;
&lt;br /&gt;
[[《積分幾何学》]]　[[《都市構造分析》]]　[[《地理的最適化》]]　[[《ウォードロップの原理》]]　[[《地域間相互作用モデル》]]　[[《ODの調査》]]　[[《地理情報システム》]]&lt;br /&gt;
&lt;br /&gt;
===システム分析・意思決定支援・特許===&lt;br /&gt;
&lt;br /&gt;
[[《システム分析》]]　[[《リエンジニアリング》]]　[[《意思決定支援システム》]]　[[《効用関数》]]　[[《データマイニング》]]　[[《過程決定計画図》]]　[[《発想法》]]　[[《モデル管理》]]　[[《アルゴリズム特許》]]&lt;br /&gt;
[[《モデリング》]]&lt;br /&gt;
&lt;br /&gt;
===AHP(階層的意思決定法)===&lt;br /&gt;
&lt;br /&gt;
[[《AHP》]]　[[《AHP一対比較法》]]　[[《AHP重要度算出法》]]　[[《AHP整合性尺度》]]　[[《拡張型AHP》]]　[[《AHP重要度評価法》]]　[[《グループAHP》]]　[[《ANP》]]　[[《AHPの諸問題》]]　[[《大規模AHP》]]　[[《AHPの誤差》]]　[[《AHPの理論的解釈》]]&lt;br /&gt;
&lt;br /&gt;
===DEA(包絡分析法)===&lt;br /&gt;
&lt;br /&gt;
[[《DEA(包絡分析法)》]]　[[《CCRモデル》]]　[[《BCCモデル》]]　[[《効率性》]]　[[《規模の収穫》]]　[[《SBMモデル》]]&lt;/div&gt;</summary>
		<author><name>Tamura</name></author>
	</entry>
</feed>