ラグランジュの未定係数法 の履歴(No.5)
更新解きたい問題†
を、
個の拘束条件
の下で最大化・最小化したい 。
(ただし
とする。さもないと拘束条件が多すぎて、それだけで点が定まってしまう。)
実際には拘束条件の下で の 停留点 を探すことになる。
このような問題は以下に説明する ラグランジュの未定係数法 と呼ばれる手法を使うと簡単に解ける。
キモ†
「拘束条件下での停留点」とは、
を 次元空間のベクトルにおいて 拘束条件を満たす点として、
「 を拘束条件を破らない方向に取る限り」
の一次の変位量がゼロとなる:
という意味である。
拘束条件を破るような方向へ動かしたときに となっても構わないところが 拘束条件付き停留点探しのキモである。
ラグランジュの未定係数法†
未定係数 を用いて
という関数を構成し、
のすべての条件式を満たす点 およびその点における係数 を見つければ、 その点が停留点となる。
また逆に、全ての停留点に対して上記の条件式を満足する係数 が存在する。
すなわち、上の条件式はその点が停留点であるための必要十分条件になっている。
条件式の意味†
での微分からは元の拘束条件が現れるのみである。
一方、 での微分は、
となる。
そこで、 本の条件をすべてまとめてベクトル形式とすれば、
と書ける。これを変形すると、
となり、すなわち両者を合わせれば、
その点ですべての拘束条件が満たされており、
なおかつ
が
の一次結合で表せるような点が停留点である
と読める。
停留点の十分条件となっていること†
ラグランジュの未定係数法の条件式を満たす点 が必ず停留点となることは、 以下のように簡単に理解できる。
は拘束条件を満たすから、 を「すべての の値を変化させない方向」に取った時のみ、変位後の点も拘束条件を満たすことになる。
そのような に対しては、すべての に対して が成り立つ。
このことと条件式より、
&math( \Delta f=\bm\nabla f\cdot\Delta\bm x=\sum_i \lambda_i\underbrace{\bm\nabla g_i\cdot\Delta\bm x}_{=\,0}=0 );
となり、条件を満たす点が必ず停留点となることを確認できる。
停留点の必要条件となっていること†
逆に、すべての停留点に対して条件式を満たす が存在するだろうか?
ある点 が拘束条件下での停留点であるとすれば、
- は拘束条件を満たす
- を拘束条件を満たす方向へ動かしたときに が変化しない
が成り立つが、この 2. は、
- がすべての に垂直なら、 にも垂直である
と読み替えられる。この条件と、
- は の線形結合で表せる
が同値であることを言えば良いのだが・・・
直交補空間の概念に通じているなら
2. を
・ すべての が張る空間を とし、
・ その直交補空間を とすれば、
・ 任意の が に垂直であると言い換えることができて、これは が の直交補空間の元であることを示すから、 は を張る の線形結合で表せて、その係数が である。
ということで証明が終わるのであるが、、、前提知識の少なくて済む説明をするなら:
すべての が張る空間を とし、 そこに正規直交基底 を取る。
が の元ではない可能性を考えて、
と書く。ただし、 すなわち任意の に対して とする。
すると はすべての に垂直であるから、 と取れば、
となって、 でない限り仮定に反する。
すなわち であり、これは が の元であり、 の線形結合で表せることを意味する。