資料紹介
ラグランジュ方程式の利点
なぜこんなに有り難がるのか
計算のための断り書き
まず、今回の話に出てくる具体的な計算をすべて 2 次元で行うことを許してもらいたい。 3 次元で行うと式が非常に面倒になることはすぐ後で分かるだろう。 議論の本質が変わってしまうことはないのであまり気にしてはいけない。
またこれまで x の時間微分 dx/dt を v という記号で表していたが、これでは r や θ の時間微分などをどんな記号で表したら良いかという問題が出てきてしまう。 毎回ていねいに微分の記号を書くと式が非常に複雑になってしまい、理解しにくくなるのでこれは避けたい。
そこでそれぞれの変数の頭の上に点（ドット）をつけることで時間の 1 階微分を表現してやることにする。 もし 2 個の点がつけばそれは時間の 2 階微分を意味していることにする。
この表現法はとても助かる。
極座標における運動方程式
では本題に入ろう。 前回はニュートンの運動方程式をラグランジュの方程式に変形することを行った。 そのためにわざわざラグランジアンなどという量も定義した。 なぜニュートンの運動方程式をこのような形式に書き直す必要があるのだろうか？ それにはちゃんと理由があり、それを知ったらちょっと感動するはずである。
デカルト座標で成り立つニュートンの運動方程式は以下のようになっている。
ここではベクトル表記で書いたが、x, y の成分に分けると次のようになる。 ついでに力 F をポテンシャルエネルギー V の微分で表しておいた。
非常にすっきりした分かりやすい式である。 では、これを別の座標系を使って表した場合にはどのように書けるのだろうか？ 例えば極座標 ( r, θ ) の場合を考えて、座標変換してやろう。
まず左辺の変形から試みる。 極座標とデカルト座標の関係式は
であるので、これを左辺の x, y に代入してやって、時間で 2 階微分してやればいい。 r も θ も時間に依存しているので積の微分公式を使ってやらねばならず、ややこしいことになる。
次に右辺にある V の偏微分についてだが、これは第 1 部で説明した 微分の座標変換 を使ってやればいい。 それで結局、全体では
のようになる。 もはや一つの式の中に r と θ が乱れ飛んで収拾がつかない。 もっと分かりやすくまとめられないだろうか？
せめて右辺だけでも整理してやろうということで、
という計算をしてやれば、意外にも左辺もかなりきれいになり、
となる。 それでもまだ一つの式の中に r と θ の両方が存在している事には変わりない。 これがあの美しかったニュートンの運動方程式の極座標系における姿である。 この非対称の醜い姿をよくまぶたに焼き付けておくといい。
ラグランジュ方程式の利点
次にラグランジュの方程式を使ってやろう。 前回求めたラグランジュの方程式は、
と表される。 まだこの表記に慣れないだろうからわざわざ成分に分けて書けば、
ということである。 これを極座標 ( r, θ ) で表現してやると次のようになる。
・・・。 えーと、何かな？
・・・。 式変形の続きを期待されても困る。 これで終わりである。
あっけなかっただろうか？ 「記号を書き換えただけではないか」と思うかも知れない。 その通り。 つまりラグランジュ方程式は、 別の座標に変換しても同じ形の式がそのまま成り立っているのである。 ニュートンの方程式のように醜い形に化けたりはしない。
タグ
All rights reserved.
資料の原本内容
ラグランジュ方程式の利点
なぜこんなに有り難がるのか
計算のための断り書き
まず、今回の話に出てくる具体的な計算をすべて 2 次元で行うことを許してもらいたい。 3 次元で行うと式が非常に面倒になることはすぐ後で分かるだろう。 議論の本質が変わってしまうことはないのであまり気にしてはいけない。
またこれまで x の時間微分 dx/dt を v という記号で表していたが、これでは r や θ の時間微分などをどんな記号で表したら良いかという問題が出てきてしまう。 毎回ていねいに微分の記号を書くと式が非常に複雑になってしまい、理解しにくくなるのでこれは避けたい。
そこでそれぞれの変数の頭の上に点（ドット）をつけることで時間の 1 階微分を表現してやることにする。 もし 2 個の点がつけばそれは時間の 2 階微分を意味していることにする。
この表現法はとても助かる。
極座標における運動方程式
では本題に入ろう。 前回はニュートンの運動方程式をラグランジュの方程式に変形することを行った。 そのためにわざわざラグランジアンなどという量も定義した。 なぜニュートンの運動方程式をこのような形式に書き直す必要があるのだろうか？ それにはちゃんと理由があり、それを知ったらちょっと感動するはずである。
デカルト座標で成り立つニュートンの運動方程式は以下のようになっている。
ここではベクトル表記で書いたが、x, y の成分に分けると次のようになる。 ついでに力 F をポテンシャルエネルギー V の微分で表しておいた。
非常にすっきりした分かりやすい式である。 では、これを別の座標系を使って表した場合にはどのように書けるのだろうか？ 例えば極座標 ( r, θ ) の場合を考えて、座標変換してやろう。
まず左辺の変形から試みる。 極座標とデカルト座標の関係式は
であるので、これを左辺の x, y に代入してやって、時間で 2 階微分してやればいい。 r も θ も時間に依存しているので積の微分公式を使ってやらねばならず、ややこしいことになる。
次に右辺にある V の偏微分についてだが、これは第 1 部で説明した 微分の座標変換 を使ってやればいい。 それで結局、全体では
のようになる。 もはや一つの式の中に r と θ が乱れ飛んで収拾がつかない。 もっと分かりやすくまとめられないだろうか？
せめて右辺だけでも整理してやろうということで、
という計算をしてやれば、意外にも左辺もかなりきれいになり、
となる。 それでもまだ一つの式の中に r と θ の両方が存在している事には変わりない。 これがあの美しかったニュートンの運動方程式の極座標系における姿である。 この非対称の醜い姿をよくまぶたに焼き付けておくといい。
ラグランジュ方程式の利点
次にラグランジュの方程式を使ってやろう。 前回求めたラグランジュの方程式は、
と表される。 まだこの表記に慣れないだろうからわざわざ成分に分けて書けば、
ということである。 これを極座標 ( r, θ ) で表現してやると次のようになる。
・・・。 えーと、何かな？
・・・。 式変形の続きを期待されても困る。 これで終わりである。
あっけなかっただろうか？ 「記号を書き換えただけではないか」と思うかも知れない。 その通り。 つまりラグランジュ方程式は、 別の座標に変換しても同じ形の式がそのまま成り立っているのである。 ニュートンの方程式のように醜い形に化けたりはしない。
いやいや、ただ話を聞いただけで信じてはいけないのだ。 本当かどうかちゃんと確かめてやる必要がある。 やってみよう。 言っておくが以下で行うのはただの確認ではなく、その過程で学べるようにと思ってわざわざ行うのである。
本当にこんなことが成り立つのか
検算してやろう。 デカルト座標でのラグランジアンは
であった。 これを極座標で表すには、ちゃんと微分の計算もしてやらなくてはならないが、比較的簡単に出来て、
となる。 これをラグランジュ方程式に当てはめて計算してやろう。 まず初めの式は簡単に
となって確かに先ほど見たのと同じ結果になる。 2 番目の式も次のように計算できる。
おお！ こうして前に苦労して求めたのと同じ結果がこんな単純な手続きで出て来てしまうのである。 すげーよ。
種明かし
なぜこのようなことになっているのかをちゃんと説明しておこう。 理由さえ分かってしまえば不思議でもなんでもない。
例えばこれから、変数 という組から という組への座標変換をするとしよう。 具体的に書けば ( x, y, z ) → ( r, θ, ) などの変換をしようと言っているに過ぎないわけだが、デカルト座標とか極座標にこだわりたくないので、このような や といった表記を使って議論することにする。
変換するということはつまり、 は のような関数として表せるということだ。 このようにしてわざわざ変数を並べるのは面倒なので、このことを次のように表す。
これは単純なことを表しているに過ぎないが、こういう表現に慣れて欲しいのでわざわざ書いてみたまでである。 さて、これからやりたい事は、 を使った場合に
が成り立っているとしたら、変換後にも同じ形の
が成り立っていることを数学的に示したいのである。 今回のように予め結果を確信している場合には、絶対に旨く説明できるはずだと信じて大胆に行動した方が分かり易かろう。
変換後の L は L( , , t ) という形の関数になっているはずである。 これを (3) 式に当てはめてみよう。 とは言うものの、(3) 式では L を や で偏微分しているので、 L が , の関数であって、さらに , が の関数であることを考慮して、ちょっと面倒な形で表さないといけない。 まず (3) 式の第 1 項に含まれる偏微分だが、
となるだろう。 この計算の中になぜ×印が書いてあるかと言うと、初めに言ったように、 は の関数であって には依存していないから、本来なら初めから書かないでおいた方がいいくらいのものだからである。 しかし初心者に「何をしてはいけないか」が良く分かるようにわざわざ書いてみた。 他方、 は に依存している可能性があるので残った項についてはこのままにしておこう。
次に第 2 項の偏微分だが、
となるだろう。 これらを使って (3) 式を書き直してやると、
となる。 この式から絶対に (4) 式が導かれてくるはずなのだ・・・と信じてみよう。 うーん、何て強引な論理展開。（笑） そう言われてみれば、(4) 式を作るための部品が見えているではないか！ 残したい部分に丸を付けてやると分かり易いだろう。 それと、邪魔な項は四角で囲っておいてやろう。
これを実現するためには、次のような関係が成り立ってさえいれば良さそうだ。
もしこれらが正しければ、 2 つの邪魔な項については互いに打ち消しあって消えてくれるし、残したい 2 つの項はまとめることが出来て、
となるであろう。 ここで をどの で偏微分しようとも 0 となるなんてことは、空間から空間、平面から平面に対応させるような普通の健全な座標変換では起こらないので、これで (4) 式が再現できることになる。
残るは (5) 式、(6) 式が正しい事を示してやることである。 まず (5) 式から確かめてやろう。 を t で微分してやろうとすると、 は の関数なので、
と書ける。 (7) 式を で偏微分してやると、右辺は k = i となる項の係数のみが残る事になる。 つまり、
である。 これは (5) 式そのものである。 これで (5) 式は証明された。
次に (6) 式が正しいことを示そう。 (6) 式の左辺は次のように変形できる。
1 行目の変形は (5) 式を使っただけであり、 2 行目への変形は偏微分 を t で微分するにあたって、この偏微分が の関数になっていることを使ったのである。 なぜそんなことが言えるのかって？ が のみの関数なのだから、 を で偏微分したものだって の関数以外になりようがないではないか。 3 行目ではそれを整理して書き直しただけである。
これはこのまま置いておいて、(6) 式の右辺を変形してみたらどうだろうか。
1 行目の変形は (7) 式を当てはめただけである。 そして 1 行目の右辺のカッコの中身を で偏微分するわけだが、偏微分というのは、その関数の中にあからさまに含まれている変数で微分するという意味だから、 を で偏微分してみたところで 0 となるだけで意味がない。 それで 2 行目のようになるのである。 何と言うことか！ これは先ほどの (6) 式の左辺の変形結果と同じである！ こうして (6) 式も正しいという事が証明できた。
無事に目的を果たす事が出来た。 初めにやった具体的な例ではデカルト座標から極座標への変換を示したわけだが、どんな座標系に変換しようともラグランジュ方程式は形を変えることがないことが言えたのである。
本当の仕掛け
数学的なトリックに目を奪われてしまったかも知れない。 ひょっとして「過去の天才たちはどうしてこんな巧妙な仕掛けに気付く事が出来たのだろうか」と溜息を漏らしてしまっているだろうか。 しかし、このラグランジュ方程式は割と自然な成り行きで発見されたのである。 それは第 3 部で説明することになる。 ...
