ルジャンドル変換

関数 f ( x ) {\displaystyle \scriptstyle f(x)} のルジャンドル変換 f ( p ) = max x ( p x f ( x ) ) {\displaystyle \scriptstyle f^{\star }(p)=\max _{x}(px-f(x))} (ただし p = f ˙ ( x ) {\displaystyle \scriptstyle p={\dot {f}}(x)} )は、 x {\displaystyle \scriptstyle x} における f ( x ) {\displaystyle \scriptstyle f(x)} の接線の切片に対応している。

ルジャンドル変換(ルジャンドルへんかん、: Legendre transformation)とは、凸解析において、関数の変数をその微分に変えるために用いられる変換である。このとき実数関数 f(x) は微分可能でなくてもよいが連続関数だとする[1]

名前はフランス数学者アドリアン=マリ・ルジャンドルに因む。ルジャンドル変換は点と双対性、つまり凸な関数 y = f (x)(x, y) の点の集合によって表現できるが、それらの傾きと切片の値で指定される接線の集合によっても等しく充分に表現できることに基いている。

凸関数をルジャンドル変換する際、変換前の関数が保持している情報は、変換後の関数においても完全に保たれる[2]解析力学においてはこの性質を利用して、ラグランジアンからルジャンドル変換によってハミルトニアンが得られる(⇒#解析力学)。物理学等において他にも広く応用されており、熱力学における熱力学関数間の変換などにも用いられる(⇒#熱力学)。

ルジャンドル変換をより一般化したものはルジャンドル=フェンシェル変換と呼ばれる。

定義

関数 f(x) のルジャンドル変換 f*(p) は次のように定義される[3][4]

f ( p ) sup x ( p x f ( x ) ) {\displaystyle f^{*}(p)\equiv \sup _{x}(px-f(x))}

または

f ( p ) inf x ( f ( x ) p x ) . {\displaystyle f^{*}(p)\equiv -\inf _{x}(f(x)-px)\,.}

ここで supx および infx は変数 x を動かしたときの上限および下限を表す[注 1]

また、下に凸な関数に対して等価な、別の定義として以下がある[6]

f ( p ) p x f ( x )  at  x  such that  f ( x 0 ) p f ( x + 0 ) . {\displaystyle f^{*}(p)\equiv px-f(x)\quad {\text{ at }}x{\text{ such that }}f'(x-0)\leq p\leq f'(x+0)\,.}

at 以下は、「f′(x − 0) ≤ pf′(x + 0) を満たすような x を、一意でなければどれでもいいから代入せよ」という意味である。

変換される関数 fx 以外の変数を持っていてもよく、多変数関数に対しては各変数についてルジャンドル変換をすることができる。

性質

その定義から明らかなように、関数 g(x) = pxf(x) の上限が定まる場合のみ、ルジャンドル変換によって新しい関数 f*(p) を与えることができる[7]

それが定まる場合でも、もととなる関数 f(x)凸関数でない場合、新たに定義された関数 f*(p) は逆変換してももとの関数 f(x) へは戻らない[7]。したがって応用上は初めから関数 f(x) は凸であると仮定することが多い。

滑らかな関数に対する変換

もとの関数 f(x) の一階の導関数 f'(x)x について連続であり単調増加する場合、すなわち関数 f(x) が下に凸で滑らかな場合、関数 g(x) = pxf(x) の上限となる x は、g(x)x の一階の導関数 g'(x) が 0 になる点であるから、ルジャンドル変換は次のように書き直せる。

f ( p ) = p x ( p ) f ( x ( p ) ) . {\displaystyle f^{*}(p)=px^{*}(p)-f(x^{*}(p))\,.}

ここで関数 x*(p)f(x) の導関数 f'(x)逆関数である:

x ( p ) = f 1 ( p ) . {\displaystyle x^{*}(p)={f'}^{-1}(p)\,.}

これは方程式 f'(x) = p の解である。

また、f もそのルジャンドル変換 f * も2階微分可能なら、両者は逆数の関係にある[8]。すなわち

f ( x ) f ( p ) = 1 . {\displaystyle f''(x)f^{*}{''}(p)=1\,.}

ただし xpp = f′(x) を満たすとする。

逆変換

関数 f(x) のルジャンドル変換 f*(p) に対して再びルジャンドル変換を施した関数を f**(x) とする:

f ( x ) sup p ( x p f ( p ) ) . {\displaystyle f^{**}(x)\equiv \sup _{p}(xp-f^{*}(p))\,.}

f(x) が下に凸であれば、f**(x) はもとの関数 f(x) に等しい:

f ( x ) = f ( x ) . {\displaystyle f^{**}(x)=f(x)\,.}

つまりルジャンドル変換の逆変換はルジャンドル変換そのものとなる。

簡単な証明として、関数が滑らかな凸関数である場合についてこのことを示す。まず関数 f を 2 回ルジャンドル変換をすると以下のようになる。

f ( x ) = x p ( x ) f ( p ( x ) ) = x p ( x ) p ( x ) x ( p ( x ) ) + f ( x ( p ( x ) ) ) . {\displaystyle {\begin{aligned}f^{**}(x)&=xp^{*}(x)-f^{*}(p^{*}(x))\\&=xp^{*}(x)-p^{*}(x)x^{*}(p^{*}(x))+f(x^{*}(p^{*}(x)))\,.\end{aligned}}}

p*(x)d/dp f*(p) の逆関数であり、x*(p)f'(x) = d/dx f(x) の逆関数なので、

d f ( p ) d p = x ( p ) + p d x ( p ) d p d x ( p ) d p d f ( x ( p ) ) d x = x ( p ) . {\displaystyle {\begin{aligned}{\frac {\mathrm {d} f^{*}(p)}{\mathrm {d} p}}&=x^{*}(p)+p{\frac {\mathrm {d} x^{*}(p)}{\mathrm {d} p}}-{\frac {\mathrm {d} x^{*}(p)}{\mathrm {d} p}}{\frac {\mathrm {d} f(x^{*}(p))}{\mathrm {d} x}}\\&=x^{*}(p)\,.\end{aligned}}}

p*(x)x*(p) の逆関数でもあり、x*(p*(x)) = x が成り立つ。 このことから、f** はもとの関数 f に等しいことが示される。

f ( x ) = x p ( x ) p ( x ) x + f ( x ) = f ( x ) . {\displaystyle f^{**}(x)=xp^{*}(x)-p^{*}(x)x+f(x)=f(x)\,.}

ヤングの不等式

以下の不等式が成り立つ。この種の不等式はヤングの不等式と呼ばれる。

f ( p ) + f ( x ) p x . {\displaystyle f^{*}(p)+f(x)\geq px\,.}

ルジャンドル変換の定義より、

f ( p ) = sup x ( p x f ( x ) ) p x f ( x ) {\displaystyle f^{*}(p)=\sup _{x}(px-f(x))\geq px-f(x)}

であるから、f(x) を両辺に足せば上述の不等式が成り立つ。

その他

  • f (x) が下(上)に凸なら f *(p) も下(上)に凸。すなわちルジャンドル変換は凸性を保持する[8]
  • f (x)左微分と右微分が異なる( f (x) が折れ線となる)点は、f *(p)p に関して1次関数(直線)となる領域に対応する。逆に f (x)x に関して直線となる領域は f *(p) が折れ線となる点に対応する[9]。これは感覚的には、 f ( x ) {\displaystyle f''(x)\to \infty } f ( p ) 0 {\displaystyle f^{*}{''}(p)\to 0} が対応するとも解釈できる。
  • 象徴的に書けば、以下のように xp について対称な関係がある。
    p = d f d x , x = d f d p {\displaystyle p={\frac {\mathrm {d} f}{\mathrm {d} x}},\qquad x={\frac {\mathrm {d} f^{*}}{\mathrm {d} p}}}
  • f (x) の逆関数 f -1 (x) のルジャンドル変換は以下である(後述のヘルムホルツエネルギーマシュー関数(英語版)の関係などに応用がある)。
    ( f 1 ) ( p ) = p f ( 1 / p ) . {\displaystyle (f^{-1})^{*}(p)=-pf^{*}(1/p).}

多変数関数

多変数関数に対しては、その一部の変数に関してだけのルジャンドル変換を考えることができる。

2変数関数 f(x, y)x についてルジャンドル変換した関数を f*(p, y) とする。このとき、変換されない変数 y はスペクテータと呼ばれる[10]。スペクテータ y による偏微分はルジャンドル変換の影響を受けない。すなわち次式が成り立つ[10]

( f ( x , y ) y ) x = ( f ( p , y ) y ) p . {\displaystyle \left({\frac {\partial f(x,y)}{\partial y}}\right)_{x}=-\left({\frac {\partial f^{*}(p,y)}{\partial y}}\right)_{p}\,.}

応用例

熱力学

ルジャンドル変換で繋がっている熱力学関数とその変数のまとめ。

熱力学では、熱力学関数間の変換、すなわち内部エネルギー U(S, V)エンタルピー H(S, p)ヘルムホルツの自由エネルギー F(T, V) に、またそれらからギブスの自由エネルギー G(T, p) に変換する際にルジャンドル変換が用いられる。

H ( S , p ) = U ( S , V ) + p V , F ( T , V ) = U ( S , V ) T S , G ( T , p ) = H ( S , p ) T S = F ( T , V ) + p V = U ( S , V ) + p V T S . {\displaystyle {\begin{aligned}H(S,p)&=U(S,V)+pV,\\F(T,V)&=U(S,V)-TS,\\G(T,p)&=H(S,p)-TS\\&=F(T,V)+pV\\&=U(S,V)+pV-TS\,.\end{aligned}}}

ここで、V体積p圧力SエントロピーT温度である。US, V について下に凸である[11]ため、U(S, V), H(S, p), F(T, V), G(T, p) はルジャンドル変換を介して互いに等価である(同じ情報を持つ)。

ただしここではルジャンドル変換は

f ( p ) = f ( x ) x p  at  x  such that  f ( x 0 ) p f ( x + 0 ) {\displaystyle f^{\dagger }(p)=f(x)-xp\quad {\text{ at }}x{\text{ such that }}f'(x-0)\leq p\leq f'(x+0)}

という定義が用いられる[12]。この定義でも上に述べた性質はほぼ同様に成り立つが、符号や凸性の変化などがある。たとえば逆変換は

f ( x ) = f ( p ) + x p {\displaystyle f(x)=f^{\dagger }(p)+xp}

に変わるという不便さがある。しかし多変数関数をこの定義で変換した場合、凸性の上下については、変換した変数についてのみ逆転し残りの変数についてはもとのまま保持されるという簡便さがある。

熱力学では導関数の不連続性は相転移として現れる。

解析力学

解析力学におけるルジャンドル変換にThermodynamic square(英語版)を適用したときのオイラー=ラグランジュ方程式。
Thermodynamic square(英語版)を適用したときの正準方程式。

解析力学では、ラグランジアン Lハミルトニアン H に変換する際に、ルジャンドル変換が用いられる。座標を q としたときに正準運動量p = ∂L/∂·q として、ハミルトニアンは

H = q ˙ p L {\displaystyle H={\dot {q}}p-L}

定義される[13]。これによって、L(q, ·q) から H(q, p) になる。 実際これは以下の関係を満たす。

H q ˙ = p L q ˙ = 0 . {\displaystyle {\frac {\partial H}{\partial {\dot {q}}}}=p-{\frac {\partial L}{\partial {\dot {q}}}}=0\,.}

このハミルトニアンとオイラー=ラグランジュ方程式あるいは最小作用の原理を組み合わせることで正準方程式が導かれる[14]。 ハミルトニアンの全微分は、

d H = H p d p + H q d q + H t d t {\displaystyle \mathrm {d} H={\frac {\partial H}{\partial p}}\mathrm {d} p+{\frac {\partial H}{\partial q}}\mathrm {d} q+{\frac {\partial H}{\partial t}}\mathrm {d} t}

と書けるが、一方でハミルトニアンの定義より、

d H = p d q ˙ + q ˙ d p L q d q L q ˙ d q ˙ L t d t = p d q ˙ + q ˙ d p p ˙ d q p d q ˙ L t d t = q ˙ d p p ˙ d q L t d t {\displaystyle {\begin{aligned}\mathrm {d} H&=p\mathrm {d} {\dot {q}}+{\dot {q}}\mathrm {d} p-{\frac {\partial L}{\partial q}}\mathrm {d} q-{\frac {\partial L}{\partial {\dot {q}}}}\mathrm {d} {\dot {q}}-{\frac {\partial L}{\partial t}}\mathrm {d} t\\&=p\mathrm {d} {\dot {q}}+{\dot {q}}\mathrm {d} p-{\dot {p}}\mathrm {d} q-p\mathrm {d} {\dot {q}}-{\frac {\partial L}{\partial t}}\mathrm {d} t\\&={\dot {q}}\mathrm {d} p-{\dot {p}}\mathrm {d} q-{\frac {\partial L}{\partial t}}\mathrm {d} t\end{aligned}}}

となるので、ハミルトニアンの偏微分は以下の関係を満たす。この内、正準変数 p, q の偏微分に関する式をまとめて正準方程式 (canonical equations) と呼ぶ。

{ H p = q ˙ H q = p ˙ d H d t = L t {\displaystyle {\begin{cases}&{\dfrac {\partial H}{\partial p}}={\dot {q}}\\&{\dfrac {\partial H}{\partial q}}=-{\dot {p}}\\&{\dfrac {\mathrm {d} H}{\mathrm {d} t}}=-{\dfrac {\partial L}{\partial t}}\end{cases}}}

逆にハミルトニアンからラグランジアンを得る場合には、関数 L を以下のように定義し、

L = q ˙ p H {\displaystyle L={\dot {q}}p-H}

変数 p に対する偏微分が 0 になるようにする。すなわち、

L p = q ˙ H p = 0 . {\displaystyle {\frac {\partial L}{\partial p}}={\dot {q}}-{\frac {\partial H}{\partial p}}=0\,.}

結局このとき変数 ·q はハミルトニアンの運動量微分に等しくなる。

多変数の場合には、ラグランジアンのすべての一般化速度についてルジャンドル変換を施したものがハミルトニアンと呼ばれる。また部分的にルジャンドル変換をしたものはラウシアン(英語版) (Routhian) と呼ばれる[15]

正準変換

具体例

滑らかな関数の例

最も簡単な例として特異性のない関数を挙げる[16]。定数 s > 1 に対して

f ( x ) = x s s , x > 0 {\displaystyle f(x)={\frac {\;x^{s}}{s}}\,,\quad x>0}

とする。この関数 f は下に凸かつ十分滑らかである。関数 f のルジャンドル変換 f*

f ( p ) = p t t , p > 0 {\displaystyle f^{*}(p)={\frac {\;p^{t}}{t}}\,,\quad p>0}

となる。ただし t1/s + 1/t = 1 を満たす定数である。この例では f(x)f *(p) は対称な形となる。特に s = t = 2 の場合、f(x) = x2/2 はルジャンドル変換で形を変えない。

導関数に定数となる領域がある例

f (x) がある区間で1次関数(直線)となる例を挙げる[17]

f ( x ) = { x 3 + x 4 ( 0 < x < 1 ) , x 1 2 ( 1 x < 2 ) , x 2 2 x + 3 2 ( 2 x ) . {\displaystyle f(x)={\begin{cases}{\frac {x^{3}+x}{4}}&(0<x<1),\\x-{\frac {1}{2}}&(1\leq x<2),\\{\frac {x^{2}}{2}}-x+{\frac {3}{2}}&(2\leq x).\end{cases}}}

導関数 f'(x) が全領域で連続、かつ Template:Math2= で定数 f'(x) = 1 であることに注意する。この関数のルジャンドル変換は p = 1 で折れ線となり、

f ( p ) = { 1 2 ( 4 p 1 3 ) 3 / 2 ( 1 4 < p < 1 ) , 1 2 ( p 2 + 2 p 2 ) ( 1 p ) {\displaystyle f^{*}(p)={\begin{cases}{\frac {1}{2}}\left({\frac {4p-1}{3}}\right)^{3/2}&({\frac {1}{4}}<p<1),\\{\frac {1}{2}}(p^{2}+2p-2)&(1\leq p)\end{cases}}}

となる。

斉次関数

f(x)k 次の斉次関数であるなら、そのルジャンドル変換は

f x x f = ( k 1 ) f {\displaystyle {\frac {\partial f}{\partial x}}x-f=(k-1)f}

である[18]。特に k = 1 のとき、f(x) のルジャンドル変換は 0 となる。

注釈

  1. ^ より初等的な最大値や最小値を用いて
    f ( p ) = max x ( p x f ( x ) ) {\displaystyle f^{*}(p)=\max _{x}(px-f(x))}
    または
    f ( p ) = min x ( f ( x ) p x ) {\displaystyle f^{*}(p)=-\min _{x}(f(x)-px)}
    としている文献もある[5]。ただしこの場合最大値や最小値には無限大も含むこととする。

出典

  1. ^ 谷村省吾,. “ルジャンドル変換、物質情報学 1(解析力学)講義ノート 6” (pdf). 2022年8月27日閲覧。
  2. ^ 田崎 2000, p. 152.
  3. ^ 田崎 2000, p. 271, 脚注.
  4. ^ 清水 2007, p. 248.
  5. ^ 田崎 2000, p. 270.
  6. ^ 清水 2007, p. 223.
  7. ^ a b 田崎 2000, p. 272.
  8. ^ a b 清水 2007, p. 225.
  9. ^ 清水 2007, p. 240.
  10. ^ a b 二間瀬 & 綿村 2010, p. [要ページ番号].
  11. ^ 清水 2007, p. 99.
  12. ^ 清水 2007, p. 227.
  13. ^ 須藤 2010, p. 42.
  14. ^ 須藤 2010, p. 43.
  15. ^ 須藤 2010, p. 47.
  16. ^ 田崎 2000, p. 275.
  17. ^ 清水 2007, p. 234.
  18. ^ 菅野 1996, p. 152.

参考文献

  • 田崎, 晴明『熱力学 現代的な視点から』(初版)培風館〈新物理学シリーズ 32〉、2000年4月12日、302頁。ISBN 9784563024321。  付録 H. Legendre 変換、pp.270-278.
  • 須藤, 靖『解析力学・量子論』(第 2 版)東京大学出版会、2010年8月20日、274頁。ISBN 9784130626101。  5.2 ルジャンドル変換、pp.45-47.
  • 清水, 明『熱力学の基礎』東京大学出版会、2007年、221-254頁。ISBN 978-4-13-062609-5。 
  • 菅野, 礼司『微分形式による特殊相対論』丸善、1996年、152頁。ISBN 4-621-04262-9。 
  • 二間瀬, 敏史、綿村, 哲『解析力学と相対論』朝倉書店、2010年。ISBN 978-4-254-13772-9。 

関連項目

外部リンク

  • 表示
  • 編集