Lagrangeの未定乗数法

質問

例えば、これは熱力学のところで出てくる式ですが、
　　　　　ΣδN_j[ln(V_j/N_j)-1]=0　(※)
に対して、
δN=ΣδN_j=0 (※※)
という制限をつけたときに(※※)に未定乗数λをかけて(※)に加えますよね。
だったら、別にλじゃなくてもなんらかの数字例えば3をかけて加えてもいいんじゃないかと思ってしまうんですが、そうするとあとあと他の条件
　　　　　　　　 N=ΣN_j
からλを決定する必要が無くなりλ=3に勝手に決まってしまいおかしなことになってしまいます。これなら最初に
(※※)でかける数字によってλの値が一つに決められて何通りもの答えが出てきてしまいます。
　教えて頂きたいのは数学的な部分で何故、(※※)に未定乗数をかけて加えてよくて、何故数字をかけて加えてはいけないのかという理由です。そんなの当たり前だと思われるかも知れませんがよく分からないのでお願いします。

回答No.8

physicist_naka · ベストアンサー

No.4のものです。
> これは分かりました。しかし、これがλの代わりに決まった数をかけて
> はいけないという説明になっているというのが分からないのでよろしく
> お願いします。

ｈ＝ｆ－３ｇの極値を求めるにはどうしますか？
　　∂ｈ／∂ｘ＝０
　　∂ｈ／∂ｙ＝０
　　∂ｈ／∂ｚ＝０
　　・・・
としますよね。でもこれで決まるｘ、ｙ、ｚ・・・はｈ＝ｆ－３ｇの極値
を与えてもｇ＝０を満たすとは限りませんね。
ですからこのやり方ではまずいです。
それともこういう意味でしょうか。
λ＝３のとき、このやり方がまずいならλのままでもまずいのでは
ということでしょうか。これはつまりこういうことです。
　　∂ｈ／∂ｘ＝０
　　∂ｈ／∂ｙ＝０
　　∂ｈ／∂ｚ＝０
　　・・・
というやり方で、ｈ＝ｆ－λｇの極値を与えますが、このやり方で
出てくるｘ、ｙ、ｚ・・・（λの関数）がｇ＝０を満たすように
λを決めておけばいいですね。

お礼

ありがとうございます。分かりました！！そういう意味の説明だったんですね！
微分積分を含む式では普通の代数式と違って単純に何かをかけて加えるとその数字によって式の表す内容が変わってくるんですね。
つまり、例えば
x+y=1 (1)
x+2y=3 (2)
という連立方程式では(2)の式に例えば3をかけて加えてもxとyの数値に変化はないけれども、
このような微分を含む関数ではかける数字によって
x,yの数値が変わってくるというような事ですね？

回答No.2

nubou

ｆ（ｘ，ｙ，ｚ）がｇ（ｘ，ｙ，ｚ）＝０という条件付で（ａ，ｂ，ｃ）で極値をもつとき
∂ｆ（ａ，ｂ，ｃ）／∂ｘ＝λ・∂ｇ（ａ，ｂ，ｃ）／∂ｘ
∂ｆ（ａ，ｂ，ｃ）／∂ｙ＝λ・∂ｇ（ａ，ｂ，ｃ）／∂ｙ
∂ｆ（ａ，ｂ，ｃ）／∂ｚ＝λ・∂ｇ（ａ，ｂ，ｃ）／∂ｚ
である（ただし偏微分可能性等は暗黙の前提とする）
の場合について説明する

（ａ，ｂ，ｃ）から
∂ｇ（ａ，ｂ，ｃ）／∂ｓ＝０・・・（１）
を満たすように線分を引き線分の長さをｓとし線分の先を
（α（ｓ），β（ｓ），γ（ｓ））とし
ｐ＝ｄα（０）／ｄｓ，ｑ＝ｄβ（０）／ｄｓ，ｒ＝ｄγ（０）／ｄｓ
とする

ｆ（ｘ，ｙ，ｚ）がｇ（ｘ，ｙ，ｚ）＝０という条件付で（ａ，ｂ，ｃ）で極値を持つから
∂ｆ（ａ，ｂ，ｃ）／∂ｓ＝０・・・（２）

（１）より
ｐ・∂ｇ（ａ，ｂ，ｃ）／∂ｘ＋ｑ・∂ｇ（ａ，ｂ，ｃ）／∂ｙ＋ｒ・∂ｇ（ａ，ｂ，ｃ）／∂ｚ＝０・・・（３）

（２）より
ｐ・∂ｆ（ａ，ｂ，ｃ）／∂ｘ＋ｑ・∂ｆ（ａ，ｂ，ｃ）／∂ｙ＋ｒ・∂ｆ（ａ，ｂ，ｃ）／∂ｚ＝０・・・（４）

∂ｇ（ａ，ｂ，ｃ）／∂ｘ＝∂ｇ（ａ，ｂ，ｃ）／∂ｙ＝∂ｇ（ａ，ｂ，ｃ）／∂ｚ＝０
ならばｐ，ｑ，ｒは任意に取れるから
∂ｆ（ａ，ｂ，ｃ）／∂ｘ＝∂ｆ（ａ，ｂ，ｃ）／∂ｙ＝∂ｆ（ａ，ｂ，ｃ）／∂ｚ＝０
でなければならない

∂ｇ（ａ，ｂ，ｃ）／∂ｘ≠０ならば
任意のｑ，ｒにたいして（３）をみたすようなｐが存在し
（∂ｆ（ａ，ｂ，ｃ）／∂ｘ）／（∂ｇ（ａ，ｂ，ｃ）／∂ｘ）＝λとすると
（４）－（３）×λより
（∂ｆ（ａ，ｂ，ｃ）／∂ｙ－λ・∂ｇ（ａ，ｂ，ｃ）／∂ｙ）・ｑ＋
（∂ｆ（ａ，ｂ，ｃ）／∂ｚ－λ・∂ｇ（ａ，ｂ，ｃ）／∂ｚ）・ｒ＝０
ｑとｒは任意に設定できるから
∂ｆ（ａ，ｂ，ｃ）／∂ｙ－λ・∂ｇ（ａ，ｂ，ｃ）／∂ｙ＝０
∂ｆ（ａ，ｂ，ｃ）／∂ｚ－λ・∂ｇ（ａ，ｂ，ｃ）／∂ｚ＝０

いずれにしても
∂ｆ（ａ，ｂ，ｃ）／∂ｘ－λ・∂ｇ（ａ，ｂ，ｃ）／∂ｘ＝０
∂ｆ（ａ，ｂ，ｃ）／∂ｙ－λ・∂ｇ（ａ，ｂ，ｃ）／∂ｙ＝０
∂ｆ（ａ，ｂ，ｃ）／∂ｚ－λ・∂ｇ（ａ，ｂ，ｃ）／∂ｚ＝０

これはいずれも必要条件であって十分条件でない
極値の候補をすべて網羅できるが中には極値でないものもあるので
別の手段で検証しなければならない

なお
Ｆ（ｘ，ｙ，ｚ）＝ｆ（ｘ，ｙ，ｚ）－λ・ｇ（ｘ，ｙ，ｚ）
とすると
∂ｆ（ａ，ｂ，ｃ）／∂ｘ－λ・∂ｇ（ａ，ｂ，ｃ）／∂ｘ＝０
∂ｆ（ａ，ｂ，ｃ）／∂ｙ－λ・∂ｇ（ａ，ｂ，ｃ）／∂ｙ＝０
∂ｆ（ａ，ｂ，ｃ）／∂ｚ－λ・∂ｇ（ａ，ｂ，ｃ）／∂ｚ＝０
ｇ（ａ，ｂ，ｃ）＝０
は
∂Ｆ（ａ，ｂ，ｃ）／∂ｘ＝０
∂Ｆ（ａ，ｂ，ｃ）／∂ｙ＝０
∂Ｆ（ａ，ｂ，ｃ）／∂ｚ＝０
∂Ｆ（ａ，ｂ，ｃ）／∂λ＝０
とみなすことができる

多少煩雑になるが変数が４以上のとき条件式が２以上のときも上記から推察できる

お礼

ありがとうございます。
(１）よりｐ・∂ｇ（ａ，ｂ，ｃ）／∂ｘ＋ｑ・∂ｇ（ａ，ｂ，ｃ）／∂ｙ＋ｒ・∂ｇ（ａ，ｂ，ｃ）／∂ｚ＝０・・・（３）
と変形できるのがよく分からないです。

回答No.1

spinflip

f(x,y,z...)の極値を、制限g(x,y,z...)=0のもとで求める際に、
代わりに、f+λgという新しい関数を制限なしで、
極値を求めればOKというのがラグランジュの未定乗数法ですね。
この、新しい関数の極値を制限無しで求めるために偏微分すると、
(d/dx)f=-λ(d/dx)g
(d/dy)f=-λ(d/dy)g
(d/dz)f=-λ(d/dz)g
:
という一連の方程式が得られます。この左辺を
ベクトルだと思うと、これはfのグラジエントです。
fのグラジエントとは、fをポテンシャルだと思えば、
そのポテンシャルの中におかれた質点に働く力になります。
特に変数が三つの場合は、∇fですから、力学で教わった通り、
ポテンシャルによる力です(但し方向は逆ですが)。
変数が三つ以外の場合でも多次元空間での仮想的な力と
みなすことができます。

一方、右辺もベクトルだと思うと、これは、
曲線を表す方程式g=0の偏微分ですから、その曲線の
法線(曲線に直交する線)となります。(これはご存知ないかも
知れませんが、たとえば、x^2+y^2-1=0という曲線の
法線は、左辺を偏微分して並べたベクトル(2x,2y)になっています)。

すると、ラグランジュの未定定数法のビジュアルなイメージは、
猿が、ある曲線g=0にぶらさがって、すいすいとポテンシャルfの
中を動いている様子になります。ですから、猿が止まる点は、
猿に働く力が、曲線と直交する所、すなわち法線と平行に
なる所、になります。

このように、二つのベクトルが平行になれば良い(一致しなく
とも良い)のですから、その比例係数は「未定」なのです。

お礼

ありがとうございました。この一見無味乾燥そうな数式にこんな図形的なイメージを当てはめることが出来るんですね。

補足

＞一方、右辺もベクトルだと思うと、これは、曲線を表す方程式g=0の偏微分ですから、その曲線の法線(曲線に直交する線)となります。(これはご存知ないかも知れませんが、たとえば、x^2+y^2-1=0という曲線の法線は、左辺を偏微分して並べたベクトル(2x,2y)になっています)。
こんな事初めて知りました！でもどうしてそうなるんですか？質問とはずれますが説明可能ならよろしくお願いします。