「条件付停留値問題」の版間の差分

2014年12月29日 (月) 01:48時点における版

メインページ>数学の部屋>ラグランジュの未定乗数法

条件付き停留値問題への応用

ここまでで、ラグランジュの未定乗数法の根本は直交補空間の直交補空間だと述べてきましたが、 これの条件付停留値問題への応用を示しましょう。

ベクトル　 の関数 が停留値を持つのは

( 1 )

の時ですが、

ベクトル に対し、K個の条件(拘束条件)

( 2 )

がある場合、関数 が停留値を取る条件は変わってきます。これを考えてみましょう。

の微小変化量を とすると、

{{eqn||3}

構文解析に失敗 (字句解析エラー): \delta G_l({\boldsymbol x}) = \mathrm{grad}G_l({\boldsymbol x})\cdot\delta{\boldsymbol x} = 0 (l=1～K) \label{eq4}

が成り立つようにしか は変化できません。

関数 が停留値を持つということは、 の微小変化量 　に対し が変化しないということなので

構文解析に失敗 (不明な関数「\label」): \mathrm{grad}f({\boldsymbol x}) = \left(\begin{array}{c} \frac{\partial f({\boldsymbol x})}{\partial x_1} \\ \frac{\partial f({\boldsymbol x})}{\partial x_2} \\ \vdots \\ \frac{\partial f({\boldsymbol x})}{\partial x_n} \\ \end{array}\right) \label{eq5}

構文解析に失敗 (不明な関数「\label」): \delta f({\boldsymbol x}) = \mathrm{grad}f({\boldsymbol x})\cdot\delta{\boldsymbol x} = 0 \label{eq6}

ラグランジュ未定乗数法の基本部分 で示しましたように、 は が所属する部分空間の直交補空間に属します。従って、適当な係数 を使って

構文解析に失敗 (不明な関数「\label」): \mathrm{grad}f({\boldsymbol x}) = -\lambda_1\mathrm{grad}G_1({\boldsymbol x}) - \lambda_2\mathrm{grad}G_2({\boldsymbol x})\cdots -\lambda_K\mathrm{grad}G_K({\boldsymbol x}) \label{eq7}

とかけます。この式の意味は停留値を取るような はこの式を満たすということであり、 が であらわされる部分空間全体になるわけではありません。

構文解析に失敗 (不明な関数「\label」): g({\boldsymbol x}) = f({\boldsymbol x}) + \lambda_1 G_1({\boldsymbol x}) + \lambda_2 G_2({\boldsymbol x})\cdots +\lambda_K G_K({\boldsymbol x}) \label{eq8}

とすると、

構文解析に失敗 (不明な関数「\label」): \mathrm{grad}g({\boldsymbol x}) = 0\label{eq9}

つまり \ref{eq2} と \ref{eq9} が停留値の条件となります。式 \ref{eq2} と \ref{eq9} は 個の方程式を含みますから、 全ての未知数を求められますので、これで十分です。実際、ラグランジュの未定乗数法というのは 式 \ref{eq2} と \ref{eq9} から停留点を求めるのが普通ですが、ここでもうひとひねりしてみましょう。

を関数の独立変数とすれば

構文解析に失敗 (不明な関数「\label」): h({\boldsymbol x}, \lambda_1, \lambda_2, \cdots , \lambda_K) = f({\boldsymbol x}) + \lambda_1 G_1({\boldsymbol x}) + \lambda_2 G_2({\boldsymbol x})\cdots +\lambda_K G_K({\boldsymbol x}) \label{eq10}

と書けます。 は拘束条件そのものなので、を にまで拡張すれば

構文解析に失敗 (不明な関数「\label」): \mathrm{grad}h({\boldsymbol x}, \lambda_1, \lambda_2, \cdots , \lambda_K) = 0\label{eq11}

が が停留値を持つ条件です。拘束条件は \ref{eq11} に取り込まれてしまった上、式 \ref{eq11} は独立変数を にまで拡大した停留点を求める方程式になっています。

つまり、ラグランジュの未定乗数法は、束縛条件付き停留値問題を束縛条件なしの停留値問題に巧妙に変換するのです。

まとめ

条件付き停留値問題を、条件なし停留値問題に変換する手順は以下の通りです。

1. 元の関数 と束縛条件 から式

を作ります。

2. 以下の式より、停留値での を求めます。