最適性条件とは?

辞典・百科事典の検索サービス - Weblio辞書

初めての方へ

参加元一覧


用語解説|動画|文献|商品|全文検索
Weblio 辞書 > 学問 > OR事典 > 最適性条件の意味・解説 

OR事典

日本オペレーションズ・リサーチ学会日本オペレーションズ・リサーチ学会

最適性条件 (非線形計画における)

読み方さいてきせいじょうけん
【英】:optimality condition

概要

非線形計画問題において, 最適解満たすべき条件, あるいは最適解になることを保証する条件総称. 通常, それらは目的関数制約関数勾配ベクトルヘッセ行列用いて表現される. 最適性条件には, 1次の最適性条件, 2次の最適性必要条件, 2次の最適性十分条件等があるが, 最も基本的なのがカルーシュ・キューン・タッカー条件である.

詳説

非線形計画問題とは,有限個の不等式等式制約条件の下で,目的関数最小(最大)化する問題であり,次のように表される.


\begin{array}{ll} \min. & f(x) \\ \mbox{s. t.} & g(x) \le 0, \\ & h(x) = 0. \end{array}(1) \,


ここで x=(x_1,\dots,x_n)^{\top} (^{\top}転置記号), 関数 f(x)\,, g(x)=(g_1(x),\dots,g_m(x))^{\top}h(x)=(h_1(x),\dots,h_{\ell}(x))^{\top}微分可能で,導関数連続仮定する.制約条件を満たす\bar{x}\,実行可能解(以下,可能解という)と呼ぶ.任意の可能解 x\, に対して f(\bar{x})\leq f(x)成立するとき, \bar{x}\,最適解と呼ぶ.しかしながら一般に非線形計画問題最適解求めることは容易ではなく局所的最適な解である局所的最適解考えることが多い.即ち,可能解 \bar{x}\, で,適当な半径 \delta>0\, を選べば,\bar{x}\,中心とする半径 \delta\, 内にあるどの可能解 x\, に対してf(\bar{x})\leq f(x)成立するものを局所的最適解と呼ぶ.

局所的最適解求める際,局所的最適解満たすべき条件局所的最適解であることを保証する条件基本的役割を果たすが,それらを総称して最適性条件(optimality condition)と呼ぶ.通常,最適性条件は目的関数制約関数勾配ベクトル \nabla f(x):=(\partial f/\partial x_1,\dots,\partial f/\partial x_n)^{\top} や,(i,j)\, 成分\partial^2 f/\partial x_i\partial x_j で定義されるヘッセ行列 \nabla^2 f(x)用いて記述されるが,特に勾配ベクトルのみで表された最適性条件を1次の最適性条件(first-order optimality condition),ヘッセ行列利用して表されたものを2次の最適性条件と呼ぶ.

1次の最適性条件で,最もよく利用されるのが次の カルーシュ・キューン・タッカー条件(Karush-Kuhn-Tucker conditionKKT条件)である.

\bar{x}\,局所的最適解ならば,後述する制約想定の下で,適当なベクトル (\lambda,\mu)\in \mathbf{R}^{m+\ell}存在して,ラグランジュ関数L(x,\lambda,\mu):=f(x)+\sum_{j=1}^m\lambda_jg_j(x)+\sum_{k=1}^{\ell}\mu_kh_k(x) と定義するとき,次のKKT条件成立する.


\left\{\begin{array}{l} \nabla_x L(\bar{x},\lambda,\mu)= \nabla f(\bar{x})+\sum_{j=1}^m\lambda_j\nabla g_j(\bar{x})+ \sum_{k=1}^{\ell}\mu_k\nabla h_k(\bar{x})=0, \\ g_j(\bar{x})\leq 0,\ \lambda_j\geq 0,\ \lambda_jg_j(\bar{x})=0 \ (j=1,\dots,m),\ \ h(\bar{x})=0, \end{array}\right.


ここで,ベクトル (λ,μ)ラグランジュ乗数条件 \lambda_jg_j(\bar{x})=0相補性条件呼ばれる

一方, KKT条件次の等式不等式系が解 y\in \mathbf{R}^n を持たないことを,二者択一定理用いて言い換え条件でもある.


\nabla f(\bar{x})^{\top}y<0, \nabla h(\bar{x})^{\top}y=0,\ \ \nabla g_j(\bar{x})^{\top}y\leq 0\ (j\in I(\bar{x})).(2) \,


ただし,I(\bar{x}) := \{j:\, g_j(\bar{x})=0\}アクティブ不等式制約, すなわち点 \bar{x}\, において等号成り立つ不等式制約集合を表わす.条件(2)に現れる集合 \{y:\,\nabla g_j(\bar{x})^{\top}y\leq 0\ (j\in I(\bar{x})),\ \nabla h(\bar{x})^{\top}y=0\}実行可能集合 \{x:\,g(x)\leq 0,\ h(x)=0\} の,点 \bar{x}\, における近似集合考えられるが,実際に良い近似になるのは,何らかの仮定満たす場合に限られる.この種の仮定総称して,制約想定 (constraint qualificationCQ),あるいは正則条件と呼ぶ.キューン・タッカーの制約想定はじめとして数多く制約想定提案されているが,なかでも1次独立制約想定(LICQ)とマンガサリアン・フロモヴィッツ条件(MF 条件)が重要である.


LICQ: \nabla g_j(\bar{x})\ (j\in I(\bar{x})), \ \nabla h_k(\bar{x})\ (k=1,\dots,\ell)1次独立である.
MF条件:(i) \nabla h_k(\bar{x})\ (k=1,\dots,\ell)1次独立である.(ii) \nabla h(\bar{x})^{\top}z=0, \nabla g_j(\bar{x})^{\top}z<0 (j\in I(\bar{x}))満たす z\in \mathbf{R}^n存在する.


LICQ は MF 条件より強い仮定であり,LICQ が成立するとき,ラグランジュ乗数 (\lambda,\mu)\,一意に決まる.他方制約想定が満たされない場合でも,点 \bar{x}局所的最適解ならば,ゼロでないベクトル (\lambda_0,\lambda,\mu)\in \mathbf{R}^{1+m+\ell}存在して,次のフリッツ・ジョン条件成立する.


\left\{\begin{array}{l} \lambda_0\nabla f(\bar{x})+\sum_{j=1}^m\lambda_j\nabla g_j(\bar{x})+ \sum_{k=1}^{\ell}\mu_k\nabla h_k(\bar{x})=0,\ \ \lambda_0\geq 0 \\ g_j(\bar{x})\leq 0,\ \lambda_j\geq 0,\ \lambda_jg_j(\bar{x})=0\ (j=1,\dots,m), \ \ h(\bar{x})=0. \end{array}\right.


フリッツ・ジョン条件は,\lambda_0\neq 0 のとき(従って λ0 = 1 としてよい),KKT 条件一致する.これにより,\lambda_0\neq 0保証する条件制約想定と呼ぶことも多い.

凸計画問題場合(\bar{x},\lambda,\mu)KKT条件を満たせば,\bar{x}最適解になる.しかし一般には,KKT条件最適性保証するとは限らない最適性詳しく調べるには,2次の最適性条件が必要になる.以下において,各関数は2回連続微分可能であると仮定する.即ち,2回までの偏導関数全て連続とする.2次の最適性十分条件(second-order sufficient optimality condition)としては,次の定理がよく利用される.

今,(\bar{x},\lambda,\mu)KKT条件を満たすとする.さらに3つの条件 (a) \nabla g_j(\bar{x})^{\top}y=0\ \mbox{if}\ \lambda_j>0, (b) \nabla g_j(\bar{x})^{\top}y\leq 0\ \mbox{if}\ j\in I(\bar{x})\ \mbox{and}\ \lambda_j=0, (c) \nabla h(\bar{x})^{\top}y=0満たすゼロベクトルでない任意の y\in \mathbf{R}^n に対し,y^{\top} \nabla_x^2 L(\bar{x},\lambda,\mu)y>0成立するならば,\bar{x}\,孤立局所最適解になる.すなわち,適当な半径 \delta>0\, を選べば,\bar{x}\,中心とする半径 \delta\, 内にある任意の可能解 x\neq \bar{x}\, に対してf(\bar{x})<f(x)\,成立する.ここで,\nabla_x^2 L\,ラグランジュ関数x\, に関するヘッセ行列を表わす.

この十分条件と対をなす2次の最適性必要条件(second-order necessary optimality condition)としては,次の定理が知られている.

可能解 \bar{x}局所的最適解で,(\bar{x},\lambda,\mu)KKT条件を満たすとする.このとき,1次独立制約想定の下で,(a)-(c)満たす任意の y\in \mathbf{R}^n に対し,y^{\top} \nabla_x^2 L(\bar{x},\lambda,\mu)y\geq 0成立する.

2次の最適性十分条件は,最適性判定以外にも,安定性理論感度分析アルゴリズム収束議論において重要な役割を演じる.安定性理論(stability theory)と感度分析(sensitivity analysis)は,目的関数制約関数微小変化させたとき,最適解最適関数どのように変化するかを調べ問題であり,主としてそれらの連続性を論じるのが安定性理論微分可能性や変化率取り扱うのが感度分析である.通常,パラメトリックな問題


\begin{array}{ll} \min. & f(x,u) \\ \mbox{s. t.} & g(x,u) \le 0, \\ & h(x,u) = 0, \end{array}(3) \,


や,その一般化舞台に,最適関数 \phi (u)=\inf\{f(x,u): g(x,u)\leq 0, h(x,u)=0\}最適解集合 \Phi (u)=\{x:\, f(x,u)=\phi(u), g(x,u)\leq 0, h(x,u)=0\}挙動考察する.ただし,u\in \mathbf{R}^qパラメータであり,u=\bar{u} のとき(3)は非線形計画問題 (1) に一致するものとする

このとき, もし各関数(x, u)\, に関して2回連続微分可能で, (\bar{x}, \lambda, \mu)3つの条件

(i) 2次の最適性十分条件,
(ii) 1次独立制約想定,
(iii) 狭義相補性,

すなわち, 全ての j \in I(\bar{x})\, に対し\lambda_j > 0\,, を満たすならば, \bar{u}\, に十分近い任意の u\, に対して, 問題 (3) に対すKKT条件を満たす (x(u),\lambda(u),\mu(u))\,(\bar{x},\lambda,\mu)\,近くにただひとつ存在する.また,それらは2回連続微分可能で,条件 (i)-(iii) を満たし,x(u)\, は (3) の孤立局所最適解になる. さらに,それらの1回,2回の微分は (3) に対すKKT条件用いて計算できる.この結果ラグランジュ乗数経済学の重要な概念であるシャドープライスを表すことがわかる.



参考文献

[1] J.F. Bonnans and A. Shapiro, Perturbation Analysis of Optimization Problems, Springer, 2000.

[2] A.V. Fiacco, Introduction to Sensitivity and Stability Analysis in Nonlinear Programming, Academic Press, 1983.

[3] A.V. Fiacco and G.P. McCormick, Nonlinear Programming, SIAM, 1990.

[4] 福島雅夫,非線形最適化の基礎,朝倉書店2001

[5] 今野浩山下浩,非線形計画法日科技連,1978.

[6] 川崎英文極値問題横浜図書,2004.






最適性条件に関係した商品


最適性条件のページへのリンク
「最適性条件」の関連用語
最適性条件のお隣キーワード
モバイル
モバイル版のWeblioは、下記のURLからアクセスしてください。
http://m.weblio.jp/
» モバイルで「最適性条件」を見る
_ _   


最適性条件のページの著作権
Weblio 辞書情報提供元は参加元一覧にて確認できます。

  
日本オペレーションズ・リサーチ学会日本オペレーションズ・リサーチ学会
Copyright (C) 2012 (社)日本オペレーションズ・リサーチ学会 All rights reserved.

©2012 Weblio RSS