「よくわかる解析力学」（東京図書）サポート掲示板 †
よくわかる解析力学サポートページに戻る
- mathjaxを使って、TeX形式で数式を打てるようにしてあります。＄または＄＄（もちろんほんとは全角じゃなく半角の「ドル」です）で囲んで入力してください。
- spam避けに、httpを含む文章と、英字のみの文章は登録できなくしてあります。
「サポート掲示板3」が不具合で書き込めなくなりましたので作りました。こちらに書き込んでください。
6.2.4二重振り子 †
大学生? (2020-10-12 (月) 17:02:06)
6.59で固有ベクトルが分かりましたが、その後のθとの関係性がなぜそうなるのかよくわかりません。
- もともと、このベクトルは$\left(\begin{array}{c}\theta_1\\ \theta_2\end{array}\right)$を求めていたのですから、決まった固有ベクトルとθは比例します。 -- 前野?
- θのベクトルが、一つ目の固有ベクトルに定数を掛けた量と、量と、二つ目の量と、二つ目の固有ベクトルに定数を掛けた量になる、ということです。 -- 前野?
p321-循環座標について †
s? (2020-10-12 (月) 14:47:19)
ラグラジアンLがxiを含まない場合に、dtで微分したら、0になる理由がわかりません。xi含んでいるとなぜ0にならないのですか？
- そんなことは本には書いてないと思います。「$L(\xi,\dot \xi)$が$\xi$を含まないときには、$\xi$に対応する運動量である${\partial L\over\partial \dot \xi}$を$t$で微分したら0になる」ということは書いてあります。なぜそうなるかというとオイラー／ラグランジュ方程式からです。 -- 前野?
- すみません。サポートページを間違えました。無視してください。 -- s?
（第6刷） 7.3.1 自由に回転する剛体 †
川西? (2020-10-10 (土) 20:04:29)
p186脚注†¹⁸ の $\frac{\mathrm d \vec \omega}{dt}=\vec\omega\times\vec L$ はなぜ成り立つのでしょうか。
計算すると左辺は $\frac{\mathrm d \vec \omega}{\mathrm dt}=\dot\omega_X{\bf\vec e}_X+\dot\omega_Y{\bf\vec e}_Y+\dot\omega_Z{\bf\vec e}_Z$ になって右辺は $$\begin{eqnarray*} \vec\omega\times\vec L &=& \left(\omega_X{\bf\vec e}_X+\omega_Y{\bf\vec e}_Y+\omega_Z{\bf\vec e}_Z\right)\times\left(I_{XX}\omega_X{\bf\vec e}_X+I_{YY}\omega_Y{\bf\vec e}_Y+I_{ZZ}\omega_Z{\bf\vec e}_Z\right) \\ &=& \left(I_{ZZ}-I_{YY}\right)\omega_Y\omega_Z{\bf\vec e}_X+\left(I_{XX}-I_{ZZ}\right)\omega_Z\omega_X{\bf\vec e}_Y+\left(I_{YY}-I_{XX}\right)\omega_X\omega_Y{\bf\vec e}_Z \\ &=& -I_{XX}\dot\omega_X{\bf\vec e}_X-I_{YY}\dot\omega_Y{\bf\vec e}_Y-I_{ZZ}\dot\omega_Z{\bf\vec e}_Z \end{eqnarray*}$$ になるかと思いましたが、これだと方向が違うし、もともと両辺で次元が合っていないように見えます。ここはどう計算すればいいのでしょうか。
- すいません、この式は明らかに間違いです。 -- 前野?
- わかりました。ありがとうございます。 -- 川西?
演習問題5-3 †
数学苦手? (2020-10-07 (水) 19:43:50)
静止摩擦力が向心力となる理由がわかりません。
- 「理由がわかりません」と言うのは具体的にどういう部分ででしょう？？ 物体を円運動させるには向心力が必要で、その向心力は力であればなんでも言いわけで、静止摩擦力ではダメだとする理由はないと思うのですが。 -- 前野?
- 説明不足ですみません。向心力は中心に向かう力なので、静止摩擦力が向心力であることに違和感を感じました。 -- 数学苦手?
- 因果関係が逆です。最初から「中心」があるわけじゃありません。働いた力によって、どの点を中心に円運動するかが決まります。そのときに働いている力のうち、物体が中心周りにまわるようにするのに使われているものを「向心力」といいます（それは運動方程式からそのとき決まった「中心」を向く）。 -- 前野?
5.4章末問題 演習問題5-1 (2) †
大学生? (2020-10-07 (水) 17:58:00)
万有引力ポテンシャルの微分が位置と速度の内積で表されてますが、位置ベクトルと速度ベクトルが同じ方向でない場合、計算に違いが生じると思います。
- 計算は単純に$-{GMm\over\sqrt{x^2+y^2+z^2}}$を微分して${GMm\over(x^2+y^2+z^2)^{3\over2}}(x\dot x+y\dot y+z\dot z)$となる（$x\dot x+y\dot y+z\dot z$の部分が$\vec x$と$\vec v$の内積になる）だけなので、同じ方向でなくても成り立つと思いますが。 -- 前野?
- 納得しました。ありがとうございます -- 大学生?
（第6刷） 5.3.3 変数の消去 †
川西? (2020-10-03 (土) 10:20:28)
p136 (5.96)の左辺第2項の左側の偏微分が $\frac{\partial G_j}{\partial Q_k}\biggr|_{q_i,Q_\overline{k}}$ のように、 $q_i$ と $Q_\overline{k}$ を一定にすることになっていますが、この場合 $q_i$ 以外の $q_*$ が勝手に動いても影響ないのでしょうか。
以下のところまではわかったつもりです。
すなわち、拘束条件を解いた(5.88) $Q_i=Q_i(\{q_*\})$ を $G_j (\{q_*\}, \{Q_★\})$ に代入してから $q_i$ で偏微分したものは
$$
\frac{\partial G_j (\{q_*\}, \{Q_★(\{q_*\})\})}{\partial q_i}{\biggr|_{{\{q_\overline{i}}\}}}= \frac{\partial G_j (\{q_*\}, \{Q_★\})}{\partial q_i}{\biggr|_{\{{q_\overline{i}}\},\{Q_★\}}}~+ \sum_k\frac{\partial G_j (\{q_*\}, \{Q_★\})}{\partial Q_k}{\biggr|_{\{q_*\},{\{Q_\overline{k}}\}}}
\frac{\partial Q_k (\{q_*\})}{\partial q_i}{\biggr|_{{\{q_\overline{i}}\}}}
$$
のように表されるが、左辺は定数0の微分なので0。だから右辺も0。
（(5.96)と(5.95)では 𝑗 と 𝑘 の役割が入れ替わっている。(5.96)の $\{Q_*\}$ は(5.87)・(5.92)の $\{Q_★\}$ に対応。）
このように考えましたが、上式の右辺と(5.96)の左辺を比べると第2項の左側の偏微分に差異があって、 $q_i$ 以外の $q_*$ を一定にしないと何が起きるのかわかりませんでした。
- すいません、この微分の固定する変数は$\{q_*\},\{q_{\bar k}\}$が正しいです。 -- 前野?
- わかりました。ありがとうございます。TeXの修正もお手数をおかけしました。 -- 川西?
p.44 式(2.68) †
0? (2020-09-30 (水) 22:33:04)
左辺第一項の$\frac{dy}{dl}$は、
$$\frac{\partial L}{\partial x}=\frac{\partial}{\partial x} \left( x\frac{dy}{dl}-y\frac{dx}{dl} \right)$$
の計算の結果かと思いますが、$L$の右辺第二項の$x$微分がなぜゼロになるのでしょうか。
微小量$dx$の微小量を取るから、近似的に0ということでしょうか。
- Lは第１項に$x$を、第２項に${\mathrm dx\over\mathrm d\ell}$を含んでいます。$x$の微分と${\mathrm dx\over\mathrm d\ell}$の微分は別々に考えていて、(2.68)の第１項が$x$微分の方、第２項が${\mathrm dx\over\mathrm dl}$の方です。 -- 前野?
- 計算の手続きとしての話をすると、${\partial L\over\partial x}$は「${\mathrm dx\over\mathrm d\ell}$は一定として$x$で微分」の結果であり、第２項で出てくる${\partial L\over\partial \left({\mathrm dx\over\mathrm d\ell}\right)}$は「$x$を一定としての${\mathrm dx\over\mathrm d\ell}$微分」です。 -- 前野?
- $x$で偏微分する時は、他の文字と同様に$dx$も定数とみなして計算するということですね。理解できました、ありがとうございます。 -- 0?
- 関連してもう一つ、計算方法についての質問なのですが、p.45の式(2.72), (2.73)の左辺第2項の導出がいまいちわかりません。例えば式(2.72)について考えてみると、左辺第2項は式(2.70)の$-\frac{d}{d\ell} \left( 2\lambda \frac{dx}{d\ell} \right)$の計算結果ということでしょうか。式(2.72)の第3項が出てくるのはわかるのですが、第2項の導出がいまいちわかりません。ご教授いただけませんでしょうか。 -- 0?
- $\mathrm dx$は数ではないので「定数とみなす」のは${\mathrm dx\over\mathrm d\ell}$です。 -- 前野?
- ${\mathrm d\over\mathrm d\ell}(2\lambda{\mathrm dx\over\mathrm d\ell})$の結果が第２項と第３項です。微分のライプニッツ則${\mathrm d\over\mathrm dx}(AB)={\mathrm dA\over\mathrm dx}B+A{\mathrm dB\over\mathrm dx}$を使ってます。 -- 前野?
- 微小量は数では無く、微小量の比を取ることに意味があるという話は、前野先生の他の著書に書いてありましたね。復習しておきます。 -- 0?
- 微分のライプニッツ則を用いるということは、未定乗数$\lambda$はここでは$\ell$に関する関数($\lambda(\ell)$)であるということですか？ -- 0?
- もちろん、そうです。 -- 前野?
- このラグランジュ未定定数は、$\ell$に依存する拘束条件を導入するためのものなのだから、$\ell$の関数でないと困ります（$\lambda$で微分することで、各々の$\ell$に対する拘束条件が出てこないといけない）。 -- 前野?
- 下の方で$\lambda$が定数になっているから勘違いしたのかもしれませんが、これは運動方程式を解いたら定数であることがわかったのであって、最初は定数でないと考えなくては、欲しい拘束条件が出てきません。 -- 前野?
- なるほど、理解できました。ラグランジュ未定乗数についてきちんと理解できていませんでした。この本を通して正しい理解につなげたいと思います。 -- 0?
（第6刷） 3.6.3 一般座標におけるラプラシアン †
川西? (2020-09-26 (土) 10:15:46)
pp. 80-81 で変分原理を使って3次元極座標のラプラシアンを求める手順の最後で、(3.83)の2行目から3行目のようなを変形する、すなわち両辺を −2r² sin θ で割るべきだというのは何から判断するのでしょうか。
あらかじめ答え(3.72)を知っているなら、 −2r² sin θ で割ればちょうどいいことがわかるのですが。
ウェブサイトの「いろもの物理Tips集」の「極座標のラプラシアンの出し方いろいろ」も拝見しました。
そちらでは ∆𝑓 = 𝜌 というポアソン方程式が出てくる作用を使うやり方が載っていました。そのやり方なら、右辺が 𝜌 になるように両辺を何倍かして調整すれば、左辺が自動的に ∆𝑓 になって、ラプラシアンが求められることが理解できました。
一方、本のやり方はラプラス方程式 ∆𝑓 = 0 が出てくる作用なので、方程式は出てきましたが、右辺が0なのでラプラシアン自体を導出するために最後に両辺を何倍すればいいかわかりませんでした。
- 割る必要はありません。右辺が＝0なのですから、割っても割らなくても数式の内容は同じです。 -- 前野?
- じゃあなんで割ったのかというと、「知っている式に合わせたかった」というだけです。合わせる気がないなら割らなくてもいいです。 -- 前野?
- ありがとうございます。ラプラス「方程式」 ∆𝑓 = 0 を出したい場面なら、割っても割らなくても数式の内容が同じであることはわかります。ですが、∆ の中身が3次元極座標でどう表されるかを知らない（忘れた）のでそれを知りたいとき、つまりp.79の(3.72)を導出したいときに、どうすればいいのかが気になっていました。本に載っている方法は、そんなことは気にしなくてもいきなりラプラス方程式が出ますよ、という意味でしょうか。 -- 川西?
- つまりは直交座標の場合と（係数も含めて）一致する式を出したい、ということでしょうか。それでしたら、オイラー・ラグランジュの式を出す手前で計算をやめて、変分を$\int \delta f \triangle f \mathrm d^3\vec x$の形にまとめてから$\triangle f=0$という式を直交座標と極座標で出せば、係数も含めて同じ式（極座標の方は$r^2\sin\theta$で割った式）が出てきます。 -- 前野?
- ありがとうございます。できました。$I = \int \left(-\frac{1}{2} \vec\nabla f \cdot \vec\nabla f \right) \mathrm d^3\vec x$ の変分を計算したら $\delta I = \int \delta f \triangle f r^2 \sin \theta \mathrm dr \mathrm d\theta \mathrm d\phi$ になって、係数も含めて(3.72)と同じ ∆ が出てきました。−2r² sin θ の出所は(3.78)の先頭の −2 と $\mathrm d^3\vec x = r^2 \sin \theta \mathrm dr \mathrm d\theta \mathrm d\phi$に入っている r² sin θ の積だったのですね。よくわかりました。 -- 川西?
（第6刷） 3.6.2 懸垂線の方程式 †
川西? (2020-09-25 (金) 13:02:09)
横道にそれた疑問かもしれませんが、pp. 76-78 で(3.58)の変分から解(3.66)を求める物理的な意味について確認したいです。
まず前置きとして、長さ L （定数）の糸の両端を特定の位置 (x₁, y₁) と (x₂, y₂) に固定したときの糸の形を求めるために、【問い3-3】または【問い3-4】の答えの式を使えばよいことはわかりました。
例えば p.355 (D.47) の y = C cosh (C⁻¹x + D) + E を使えば、
(1) 曲線が (x₁, y₁) を通る。
(2) 曲線が (x₂, y₂) を通る。
(3) 曲線の長さが L である。（【問い3-3】の問題文1〜2行目の式を使う）
の3条件を用いて3つの任意定数 C, D, E の値が決まります、と理解しています。
それで本題ですが、定数項がない(3.66)のような式 y = B⁻¹ cosh (Bx + D) （面倒なので √C を B と置きました; Bは正）は任意定数が2つしかありません。
この場合、上記(1)(2)の条件（両端の座標）を用いれば B と D の値が決まり、それに伴って糸の長さ L も自動的に決まってしまう（ L は自由に選べない）のでしょうか。
つまり（(3.68)でなく）(3.58)の変分を考えて解(3.66)を求めた作業は、実はこういう問題を解いたことになるのでしょうか →「両端の座標を決めたとき、糸の長さ L をいくらにすれば糸全体の位置エネルギーが停留値になるか、そしてそのときの糸の形はどうなるか?」（ここで求まった解は、Lの値を変える実現可能な変化に対しても、懸垂線じゃない形に変える実現不可能な変化に対しても、どちらの方向にも停留する。）
そうだとすると、この問題は y 座標の原点（というか位置エネルギーの基準点）の位置に依存して解の L が（したがって糸の形も）変わりますよね。
簡単のため両端の位置を (-a, b), (a, b) とする（aは正）とします。
すると D=0 となり、 B については、
(i) b∕a がある値(*)より小さい場合、条件を満たす B はありません。この場合は停留値は存在せず、 L を大きくするほど糸全体の位置エネルギーは単調に減少し、いくらでも -∞ に近づきます。
(ii) b∕a がそれより大きい場合、条件を満たす B は2つあります。この場合は L を 2a から始めてだんだん長くしていくと糸全体の位置エネルギーが、最初は減少し、あるところで増加に転じ、またあるところで減少に転じた後いくらでも -∞ に近づきます。このとき極小値を取る L は、b/a が大きいほど小さくなります（例えばaを固定すると、bが大きい場所（高いところ）でやるほど L を小さくするべし、となる）。
のようになると思われます。
(*) X² (arsinh² X - 1) = 1 の正の解
この解釈であっているでしょうか、どこか間違っているでしょうか。
(3.58)と(3.68)の違いは、【問い3-3】や†³⁰では「yの平行移動」という簡単な説明になっているので、(3.58)の変分の意味が本当に上のような複雑な問題と考えられるのか不安になりました。
- まず最初に言っておきますが、条件が３つ、任意定数が３つあるというのが正しいやり方です。任意定数が２つしかないのは、なんらかの簡略化を行ったからです。つまり「正しいやり方」は任意定数が３つある方です。別の言い方をすると、定数項を０にするということはｙの最小値（糸の最下点）を$B^{-1}$に決めてしまったことに対応します。そういう意味では、任意定数が２個しかない計算は「間違った計算」です。 -- 前野?
- 簡略化した間違った計算だけど、その違いは「ｙの平行移動」でしかないので、まぁいいか、ということを確認したのが問い3-3だということになります。 -- 前野?
- というわけで、「ｙ座標の原点に依存して解のLが変わる」のは、本来省略してはいけない「L一定の条件」を外してしまったからだ、ということは言えそうですね。 -- 前野?
- この問題は今から思えば、最初からラグランジュ未定乗数を入れた式を書いておくべきだったかなぁ、と思うのですが、最初に書いたときは「いきなりラグランジュ未定乗数を使って、とやると読者も戸惑うかなぁ」と思ってそこを後回しにしてしまったのでした。 -- 前野?
- ラグランジュ未定乗数を入れない場合、「Lが変化する場合の変分も考えている」というのはそのとおりです。ただ、おっしゃるとおり解が正しい形では出てこないです。正しい計算をするとｙ方向に平行移動ので、ｂの値は原点ずらしで変化できるので、解はどんなLでも（短すぎると糸が渡せませんがそれ以上なら）大丈夫なんですが、そこがうまく出てこない式になってしまっているということになってます。 -- 前野?
- よくわかりました。素早いご回答ありがとうございました。 -- 川西?
2.2.1 反射の法則 †
川西? (2020-09-20 (日) 10:11:34)
初めて投稿します。私は大学生のときにサボってしまい今になって勉強している中年の者です。
第6刷を買って読み始めたところですが、わからないことがありますので教えていただければ幸いです。
p.30 の真ん中らへんに、点Rをずらしたときの入射光の経路の縮みが bR 、反射光の経路の伸びが ab と書かれています。
これについて私は、入射光の経路の縮みは (aR)sin(入射角) 、反射光の経路の伸びは (aR)sin(反射角) だと考えました。これらはそれぞれ bR 、 ab とは一致しないように思えます。
これはどのように考えればよいのでしょうか。
- すいません、確かに aR sin(角)と考える方が合理的です。bR、abとした考えは、入射光の「無限遠からaまで」と「無限遠からｂまで」を「平行光線だからここまでは同じと見てもいいだろう」という考えからですが、あまり良いやり方ではありません。 -- 前野?
- ありがとうございました。 -- 川西?
ラグランジアン の導出とダランベールの原理の関係 †
夏休み? (2020-09-17 (木) 00:30:24)
ラグランジアン を導出するには経路の変分を考えたとき、運動方程式が出てくるような物を求めるだけでラグランジアン は導出されると思うのですが、その導出とダランベールの原理の関係性がわかりません。ダランベールの原理なんて考えなくてもラグランジアン が導出できると思うのですが、この考えは間違ってますか？
- 「経路の変分を考える」という段階まで行ったら、もうダランベールの原理は要らない、と言ってもいいので、その考えは間違っていません。ダランベールの原理があるおかげで「運動方程式も変分原理で」という考えに至ることができる、と思った方がいいと思います。 -- 前野?
- ８７ページの真ん中の表で「仮想仕事の式」と同じ形で「ダランベールの原理による仮想仕事の式の拡張」ができて、その右の「？」は何だろう？という次の段階へと進むことができたわけです。 -- 前野?
- 段階を考えるのですね。理解できました。返答ありがとうございました！ -- 夏休み?
P.22ni †
FumaRu? (2020-09-11 (金) 01:45:13)
つりあいの位置や条件をポテンシャルの微分が0であることから求める際に、「力がつりあっている」のは物体が静止するための必要条件にすぎないので、つり合いの条件というのは静止または等速直線運動するための条件となると思うのですがこれは正しいでしょうか。
- 『よくわかる熱力学』の23ページを読んで上のように思ったのですが、間違えてこちらの掲示板に投稿してしまいました。お手数おかけしてすみません。 -- FumaRu?
- 『よくわかる熱力学』の23ページを読んで上のように思ったのですが、間違えてこちらの掲示板に投稿してしまいました。お手数おかけしてすみません。 -- FumaRu?
p.38 FAQについて †
Sh? (2020-08-14 (金) 21:17:29)
突然失礼いたします。
掲示板を遡ると表題の件で同様の質問をされていた方がおり、先生にご回答いただいておりましたが、
それを拝見した部外の私がちゃんと理解できているかお尋ねさせていただければと思います。質問は以下です。
(なにか)=0が結論できるのは、(なにか)の後ろについてくるもの、つまりδy(x)が独立な時だけ~
という主張の根拠はなんでしょうか。また、δy'が独立でないということもよくわからないです。(以下略)
この質問に対し、先生は
(前略)独立なのは$\delta y(x)$（$x$は任意の場所）です。$y(x+\Delta x)$を変化させると、$y'(x)$と$y'(x+\Delta x)$が連動して変わります。
そういう意味で独立ではないです。
と回答されました。
私なりに解釈したイメージなのですが、
$\delta y(x)$の形だと独立(任意?)な場所$x$で関数$y$に任意の微小変化$\delta y(x)$を与えられる⇒(なにか)は常に0でなければならない、ことに対し、
$\frac{d}{dx} \delta y(x)$では、例えばどこかの場所$x$で関数$y$に任意の微小変化$\delta y(x)$を与えたとすると、
その隣の場所$x$に対する$\delta y(x)$も変わってしまう。(こちらで勝手に決めることができない⇒独立でない)
⇒連鎖的に全ての場所$x$での$\delta y(x)$が決定してしまう。⇒うまくやれば(なにか)が常に0にならずとも式が成立する組み合わせが存在する。
以上のように認識しておりますが、何かマズい理解をしているところはあるでしょうか。
浅学、長文、冗長で申し訳ありませんが、ご回答いただければ幸いです。
- ちょっと返事遅れましたすみません。「こっちで勝手に決められない」→「独立でない」という考えはそれでしいです。 -- 前野?
- 独立でないと（なにか）＝０と主張できないということの根拠は本にも書いたとおりですが、（なにか）が０でなくとも足し算の結果が０になっている場合がありえるからです。後ろにかかっているものが独立なら、「Σ（なにか）δｙが、どんなδｙをとっても０になるということは（なにか）＝０」と結論できます。 -- 前野?
- 独立でなければ０になっても（なにか）が０でない場合の一番シンプルなのは「Ax+B×(-x)=0」で、xと(-x)は独立ではないのでA=0,B=0は言えません（A=Bは言えますが）。 -- 前野?
Liouvilleの定理と正準変換について †
FumaRu? (2020-08-11 (火) 02:18:08)
深夜遅く失礼します。題名の通りLiouvilleの定理についてです。
他の資料に「体積積分が正準変換で不変となっている」ことがLiouvilleの定理だと記述されていたのですが、本書では「時間発展に関して位相空間内の体積が変化しない」ことがLiouvilleの定理だと記述されています。これらは記述している内容は違っても本質的には同じことを指しているのでしょうか。
返信お待ちしております。
- 本書では先にLiouvilleの定理を説明してから次の章で正準変換を説明するという順番になっているので、「時間発展で不変」ということでLiouvilleの定理を導入し、後で時間発展は正準変換の一種（ハミルトニアンを正準変換の生成子にすれば時間発展になる）という順番で話しました。 -- 前野?
- 時間発展に限らず正準変換なら位相空間の体積は一定になるので、より一般的なそっちをリウビルの定理と呼ぶ場合もあるのかと思います。 -- 前野?
- なるほど。理解できました。ありがとうございます。 -- FumaRu?
p.114(5.10)について †
KYU? (2020-08-09 (日) 11:41:51)
(5.10)の一行目から二行目の変形で$d/dt(\partial L/\partial \dot{Q}_j \partial \dot{Q}_j/\partial \dot{q}_i)=d/dt(\partial L/\partial \dot{Q}_j)\partial \dot{Q}_j/\partial \dot{q}_i +d/dt(\partial \dot{Q}_j/\partial \dot{q}_i)\partial L/\partial \dot{Q}_j となるはずですがd/dt(\partial L/\partial \dot{Q}_j)\partial \dot{Q}_j/\partial \dot{q}_iではなくd/dt(\partial L/\partial \dot{Q}_j)\partial Q_j/\partial q_iとなっているのはなぜですか。
- すみません、$で囲むのを忘れていました。 --
- $d/dt(\partial L/\partial \dot{Q}_j \partial \dot{Q}_j/\partial \dot{q}_i)=d/dt(\partial L/\partial \dot{Q}_j)\partial \dot{Q}_j/\partial \dot{q}_i +d/dt(\partial \dot{Q}_j/\partial \dot{q}_i)\partial L/\partial \dot{Q}_j となるはずですがd/dt(\partial L/\partial \dot{Q}_j)\partial \dot{Q}_j/\partial \dot{q}_iではなくd/dt(\partial L/\partial \dot{Q}_j)\partial Q_j/\partial q_i$ これです。 --
- 一つ前のページの(5.7)を見てください。 -- 前野?
- 見落としていました、ありがとうございます。 -- KYU?
P.208(9.9)について †
FumaRu? (2020-08-06 (木) 18:40:13)
【FAQ】に(9.3)の両辺をpiで正しく微分すると(9.9)の式が出るとあります。私は(9.3)をqiで微分し、(9.6)の計算結果を使い、(9.9)と同様の式を出しました。計算結果は(9.3)と同じになりましたが、果たして正しい計算なのか不安になり、質問させていただきました。
- これだけではどういう計算をしたのか書いてないので、正しいのかどうか判定できません。 -- 前野?
- 「(9.3)を$q_i$で微分し」とありますが、何を独立変数として微分したのか？（$p_i,q_i$を独立変数として微分したのか、$q_i,p_i$を独立変数として微分したのか、それとも何か違う方法を用いたのか？） -- 前野?
- 言葉が足らず申し訳ありません。qとqドットを独立変数として微分しました。 -- FumaRu?
- ∂H/∂qi + Σ{(∂H/∂pj)(∂pj/∂qi)} = Σ(∂pj/∂qi)(dqj/dt) - ∂L/∂qi とqとqドットを独立変数として微分しました。 -- FumaRu?
- Texが使えないため、見づらくなってしまい申し訳ありません。右辺のΣのqドットは独立変数であると分かっていますが、見やすくするためにdq/dtと表記させていただきました。 -- FumaRu?
- その計算だと、$\sum_j {\partial H\over\partial p_j}{\partial p_j\over\partial q_i}=\sum_j{\partial p_j\over\partial q_j}\dot q_j$となりますね。${\partial p_j\over \partial q_j}$を行列と考えたときの逆行列があればいいですが、そうでなかったらこの式から${\partial H\over\partial p_j}=\dot q_j$は導けません。 -- 前野?
- 多くの場合、${\partial p_j\over\partial q_j}$は逆行列がありません（簡単な系では0です）。 -- 前野?
- 返信ありがとうございます。確かに∂pj/∂iが正則であること無意識に仮定して計算してしまっていました。理解出来ました。ありがとうございます。 -- FumaRu?
- 返信ありがとうございます。確かに∂pj/∂iが正則であること無意識に仮定して計算してしまっていました。理解出来ました。ありがとうございます。 -- FumaRu?
P.86(4.1)について †
FumaRu? (2020-08-04 (火) 17:06:12)
P.86(4.1)でつりあいの式⇄仮想仕事=0となっており、、P.61において→の証明をされています。同様に←の証明もできるのでしょうか？
- 追加の質問失礼します。P.109においてポテンシャルがV(x1-x2)になったことで∂V/∂x1=-∂V/∂x2になるとありますが、こらはどうやって導かれるのでしょうか。 -- FumaRu?
- 「つりあいの式←仮想仕事＝0」は仮想仕事の式に並進や回転など「起こり得る仮想変位」を代入してあげればいいと思います。 -- 前野?
- $V(\vec x_1-\vec x_2)$を実際に微分してみてください。$V$は$\vec yi =\vec x_1-\vec x_2)$の関数なので、まず$\vec y$で微分することになります。1次元の場合で書くと$V(y(x_1-x_2))$の微分は${dV\over dy}{dy\over dx_1}$になります。 -- 前野?
- $x_2$の微分の方は最後が${dy\over dx_2}$になるので符号が逆です。 -- 前野?
- 回答ありがとうございます。確かにおっしゃる通りになりました。 -- FumaRu?
５．１オイラーラグランジュ方程式と座標変換に関して †
田島? (2020-07-16 (木) 17:14:32)
ｐ１１２の（５．４）の形で書かれている変換に関して、オイラーラグランジュ方程式が共変になることの証明で、ｐ１１４の（５．１２）で∑(∂Ｑｊ/∂ｑi)（∂ｑi/∂Ｑｋ）＝δjkとありますが、この式はどうして成り立つのでしょうか？（（５．４）のヤコビ行列の逆行列が∂ｑi/∂Ｑｋであるのはどうやって確かめられるのでしょうか？）
- $Q$を$q$で表して、その$Q$の中の$q$を$Q$で表したとします。つまり、$Q_i(\{q_*(\{Q_\star\}))$という式です。この式は$Q_i$そのものなので、を$Q_j$で微分したら、$i=j$なら1、そうでなかったら0になります。 -- 前野?
- その式が${\partial Q_i\over\partial Q_j}=\sum_k {\partial Q_i\over\partial q_k}{\partial q_k\over\partial Q_j}$です（つまりは偏微分の式です）。 -- 前野?
- 回答ありがとうございます。いろいろうっかりしておりました -- 田島?
６．４連続な物体への極限に関して †
田島? (2020-07-03 (金) 20:13:13)
①Ｎ＋１本のばねがあり、ばね1本当たりのばね定数（Ｎ／ｘ）をｋ＝κ(N+1/L)としてありますが、κは単位長当たりのばね定数でしょうか？単位はどうなっているのでしょうか？
②そのあと、前の章でＮ個のモードの調和振動子に分解した解の角振動数に関してＮ→∞としたときの角振動数（６．８７）を出し、「この解の1個1個のモードを見ると、個々の点が振幅Ｃｐsin(pπ/Lx)角振動数ωｐで振動していると考えられる」とありますが、sin関数の式はどうやって導出されたものかがピンときませんでした。また、このsinの式は時間に依存しない（ｔを含まない）ように見えますが、それは私の勘違いでしょうか？教えてください。
- バネ定数は長さに反比例するので「単位長さあたりのバネ定数」というのはちょっと変な言葉ですが、κは「長さが1m（単位長さ）のときのバネ定数」と思ってください。長さが$L$から${N+1\over L}$になったのでばね定数が${N+1}$倍になったと考えてもいいです。 -- 前野?
- 単位は合うように決めます。kが[N/m]だからκは[N]です。 -- 前野?
- sinは(6.81)からいて、それの極限を取っていった結果が(6.87)ですが、(6.81)がわからないということでしょうか？？ -- 前野?
- それは(6.79)でバネ定数に対応するものの中に$\sin^2$がいたからですが、その部分がわからないということなのかな？？（ちなみにこれはばね定数なので、もちろん時間に依存しません）。 -- 前野?
- 回答ありがとうございます。６．４．１の冒頭でやっていることの意味を理解していないのかもしれません。 -- 田島?
- Ｎ→∞の極限をとることは長さＬの弦においてどのような状態を考えていることになるのでしょうか？また、Ｎ→∞の極限をとった後、各調和振動子の角振動数の式が（６．８７）で表されること、個々の振幅は弦の中での位置ｘに応じて異なることは理解しましたが、この振幅の式はどうやって導出されるのでしょうか？（６．８０）で表されたＮ個の調和振動子の運動方程式から出そうとしましたがうまくできません。 -- 田島?
- Nが有限のときは、実際には連続的な物質（原子レベルなら不連続でしょうけどそこは問わないことにして）である弦を、いったん「N個の不連続な質点」という架空の物体に置き換えているわけです。N→∞は、「元々考えたかったものに戻す」ということです。 -- 前野?
- 別の言い方をすれば「N→∞にしたものこそが弦（それまでは「弦もどき」）」ということになります。 -- 前野?
- (6.80)では$Y_i$という変数を使って書いていて、$Y_i$のそれぞれが調和振動子になってます。よって$Y_i=A_i \sin(\omega_j t + \alpha_i)$のような振動をすることがわかります。 -- 前野?
- $Y_i$の振幅である$A_i$は、運動方程式からは決まりません。初期条件で決まります。これは普通の単振動でも同じで、運動方程式からは振幅と初期位相は出てきません。 -- 前野?
- $Y_i$が求まれば、それに対応して$y_i$が(6.78)の関係を使ってわかります。 -- 前野?
- 初期条件により、振幅が変動することを忘れておりました。そこを含め、疑問点を解決できました。お答えいただいてありがとうございました。 -- 田島?
６．２．４2重振り子に関して †
田島? (2020-07-01 (水) 21:16:15)
ｐ１５２の（６．５３）の形のラグランジアンをｐ１５１の同時対角化の方法論を用いて変形して、最終的に（６．４９）の形つまり、（６．４２）と同形のラグランジアンを得ているわけですが、最後の分析で、θ１：θ２＝√ｍ：√Ｍでという条件が常に成り立つのがよくわかりません。その後ろの角速度ωをどうやって得たかは、√ばね定数/質量にあたるのが（６.49)のＬの対角成分に並ぶ固有値だからですよね？
- その上で固有ベクトルを$\left(\begin{array}{c} \sqrt{m}\\ \sqrt{M}\end{array}\right)$と求めています。 -- 前野?
- それはつまり、$\left(\begin{array}{c}\theta_1\\ \theta_2\end{array}\right)$がそれに比例するということです。 -- 前野?
- 回答ありがとうございます。固有ベクトルの意味を忘れていました。 -- 田島?
６．２．４の2重振り子のポテンシャルの平衡点に関して †
田島? (2020-07-01 (水) 18:15:29)
ｐ１５２の６．２．４2重振り子において、（６．５２）を得た後、ポテンシャルＵの平衡点を求め、ｐ１４１での「テイラー展開して、3次以上の項を無視★」という技法を用いていますが、今回の場合、Ｕ＝-Mgl（cosθ1＋cosθ２）という形で、第1項、第2項それぞれ独立して★がつかえますが、もしＵの項にθ1θ2のような項が入っていたら、どうするのでしょうか？
- その場合は６．２節の連成振動と同じで、３次以上を無視した後で行列の対角化を使って固有振動の重なりに書き換えていきます。 -- 前野?
- 回答ありがとうございます。連成振動を読み直します。 -- 田島?
６．１．２微小振動の単振動の微分方程式の解に関して †
田島? (2020-07-01 (水) 18:06:37)
ｐ１４１の６．１．２微小振動において、ポテンシャルＵ（ｘ）をテイラー展開し、釣り合い点において、3次以上の項を無視した運動方程式（６．７）を導出し、それを解いて、（６．８）を得ていますが、（６．７）の解は、ｘ－ｘ０αexp(+kt)+βexp(-kt)(ただし、ｋ＝√ーU``(x0)/m、α、βは初期条件により定まる任意定数)の形ですが、（６．８）の形は任意定数がＣのみの1つで表現されています。これは、どういうことなのでしょうか？
- （６．８）は複号$\pm$を含んでいて、＋とーの２項があります。省略して書いてますが、$C_+ e^{+\alpha}+C_-e^{-\alpha}$のように本当は２項あると思ってください。 -- 前野?
- もう1つ同じページに関して質問があるので、ここに書かせてください。「テイラー展開して、3次以上の項を無視」という近似の手法の説明で、途中Ｕ‘‘（ｘ０）＞０との条件を課していますが、これはなぜ必要なのでしょうか？Ｕ‘‘（ｘ０）＜０では、（６．７）のーＵ‘‘（ｘ０）が正になって、釣り合い点からどんどん離れる（飛び去る）から近似が成り立たなくなるという理解でよいでしょうか？ -- 田島?
- 1件目に関して理解しました。各項を強調するためということですね。 -- 田島?
- もう一つの質問についてはその通りです。不安定な釣り合い点の回りで展開しても意味はありません。 -- 前野?
梯子の釣り合いの条件式をポテンシャルを用いて導出する方法 †
田島? (2020-06-26 (金) 20:25:55)
ｐ６８で梯子が釣り合う条件をポテンシャルを用いて導出しています。
Ｆ（梯子に加える力）が保存力だとして、そのポテンシャルを出し、拘束条件（床と壁に梯子は常に接触する）を考慮して、全ポテンシャルが極値をとる条件を出していますが、ポテンシャルが極値も持つならば、どうして力が釣り合っているといえるのでしょうか？
- ポテンシャルが極値ということは仮想仕事が0というのと同じことですが、それはわかっているでしょうか。 -- 前野?
- 回答ありがとうございます。わかっておりませんでした。ポテンシャルが極値をとるならば、そこからδｒ＝（δｘ,δｙ）という微小変位をしてもＦ・δｒ＝０となるということですか。あまりピンときません。 -- 田島?
- ポテンシャルの定義は$F_x=-{\partial U\over\partial x}$のようになることなので、$\vec F\cdot\delta \vec x=0$は$U$の微分が０というのと同じことです。 -- 前野?
- ありがとうございます。納得しました。 -- 田島?
３．３仮想仕事の原理を使う例題に関して †
田島? (2020-06-26 (金) 16:06:33)
ｐ６５において、
①（３．２５）（２ｘ）＾２＋（２ｙ）＾２＝Ｌ＾２というはしごが床と壁に接触しながら動くという拘束条件から、xδｘ＋ｙδｙ＝０★という式を出していますが、これは数学的にはどのような操作を行って導出したものなのでしょうか？
②また、（３．２５）の少し上で、はしごに対して重力がする仕事はーｍｇδｙ（この状況ではδｙは負である）とありますが、δｙ＜０はどこからわかるのでしょうか？x,yはともに正であり、δｘ＞０に動かしているので、★よりδｙ＜０としましたが、式を経由せずともわかるような（自明な）事実なのでしょうか？
- $(2x)^2+(2y)^2=L^2$と$(2(x+\delta x))^2+(2(y+\delta y))^2 =L^2$が両立すると考えれば出ます。 -- 前野?
- そもそも実はδｙは正だろうが負だろうが構わないのですが、重力が正の仕事をするだろう、と考えればδｙが負だろうな、とわかります。 -- 前野?
- ありがとうございます。 -- 田島?
- 自分でも考え直してみたのですが、（２ｘ）＾２＋（２ｙ）＾２＝Ｌ＾２の左辺の関数をｆとして、∂ｆ/∂ｘ＋∂ｆ/∂ｙ＝ｄｆという全微分の式を考えておいて、ｆ＝Ｌ＾２で一定であるからｄｆ＝０としても同じことでしょうか？ -- 田島?
- それは同じ計算です。 -- 前野?
- ありがとうございます。 -- 田島?
２．３．２一般的な図形の等周問題について †
田島? (2020-06-25 (木) 19:05:29)
ｐ４５において、「長さｌ１を式の中で条件に入れていないからである。」とあり、別の言い方をすれば、「dl^1=dx^2+dy^2」を加えればよいとあるのですが、これが図形の全周＝ｌ１と同値（等価）な条件になることがよくわかりません。
私なりに考えましたが、ｄｌが微小長さに対応する条件が「dl^1=dx^2+dy^2」であり、ｐ４４の（２．６６）のようにｌ＝０～ｌ１まで積分してあることと併せて全周＝ｌ１を考慮したことになっていると思いました。
- 条件をつけてないときの$\ell$は「長さ」という意味を持っていません。$d\ell^2 =dx^2+dy^2$という条件をつけることで、$\ell$が長さという意味を持つわけです。というわけで、その考えでOKです。 -- 前野?
- 回答ありがとうございます。 -- 田島?
２．３．３最速降下線に関して †
田島? (2020-06-25 (木) 18:57:02)
ｐ４９の†２４「ｔ＝０においてdy/dx＝∞」とありますが、これはどこからわかるのでしょうか？どうもピンときません。
- $C=0$は$y=0$ということで、それは解にはならない（ということは$C$は０ではない量で、$y=0$でも$y{dy\over dx}$が０ではない）ということからわかります。 -- 前野?
- 回答ありがとうございます。納得しました。 -- 田島?
レヴィ・チビタ記号について †
理科大学生? (2020-06-23 (火) 06:13:53)
(i,j,k)が(1,2,3)の偶置換とは
$~\sigma = \left(\begin{array}{ccc} 1 & 2 & 3 \\ i & j & k \end{array} \right)$
における置換σが偶置換であるということでしょうか。
- Tex形式で入力しようとしましたが、自分の力不足で表現されておらず、すみません。 -- 理科大学生?
- そういう事です。TeX直しました。 -- 前野?
- なるほど、毎回ご返信ありがとうございます。 -- 理科大学生?
作用の並進不変性について †
理科大学生? (2020-06-20 (土) 14:47:09)
P.194で並進不変性を仮定した時に、経路が同等になりそうなのは分かりました。この時に、確かになりそうではあるのですが、ハミルトンの主関数の値が変化しないのは何故でしょうか。
- ハミルトン主関数は「経路」に沿っての作用の積分で、並進不変性があるなら、経路を全部その方向に移動させても「作用の積分」は同じになる、ということになります。 -- 前野?
- なるほど、x、t平面に対して、ハミルトニアンを垂直な軸にとった時に確かに積分を考えるとなりますね。ありがとうございます。 -- 理科大学生?
- なるほど、x、t平面に対して、ハミルトニアンを垂直な軸にとった時に確かに積分を考えるとなりますね。ありがとうございます。 -- 理科大学生?
誤植について †
理科大学生? (2020-06-19 (金) 09:20:05)
第5刷ですが、P.207の(9.4)式の下が「*(アスタリスク)→★(星)」だと思われます。
(6.60)式について †
理科大学生? (2020-06-18 (木) 00:01:53)
ωを(6.57)式のλの値を用いていますが、これは(T^t)MT=EになるようなTを固有ベクトルに既に(6.59)式で選択済みであるということなのでしょうか。
- (T^t)MT=Eになるようにλに依存せずに変化できるので、(6.49)式の最後はλが求まった瞬間に立式できそうですね。自己解決しました。 -- 理科大学生?
- (T^t)MT=Eになるようにλに依存せずに変化できるので、(6.49)式の最後はλが求まった瞬間に立式できそうですね。自己解決しました。 -- 理科大学生?
無題 †
理科大学生? (2020-06-12 (金) 04:21:41)
2016年に先生のご返答を確認して(5.96)式が成り立つことを理解しました。念のため確認ですが、微積分学的には陰関数定理ということでいいのでしょうか。
- 追記ですが、P.136です。Gがqの関数で、0の値から動かない定数関数(ただし、Gの中にQが残っているので注意しなければならない)になったのですね。 -- 理科大学生?
- 陰関数定理と考えていいです。 -- 前野?
P.92 †13について †
理科大学生? (2020-06-09 (火) 01:24:22)
(4.14)式の表面項自体は(4.12)式の右辺第２項の変分の表面項であるから、P.98の議論を根拠としてここを無視しても良いというのはどうかなと思うのですが、ここについて教えていただきたいです。
- そういう意味では、(4.12)の表面項も無視してます。我々の目標は運動方程式なので、運動方程式に効かない部分は一貫して無視してます。 -- 前野?
p.337_付録C.1_練習問題【問いC-1】について †
kf? (2020-05-26 (火) 16:52:38)
p.337_付録C.1_練習問題【問いC-1】の解答(p.369)について、(C.12)をレヴィ・チビタ記号で計算したものは掲載されていましたが、(C.13)についての解答はありませんでした。
サポートページのログ及びこちらの掲示板にもないようでしたので、こちらの問題の解答を教えていただきたいです。
万一既に掲載済みで私が気が付いていないだけでしたら申し訳ありません。
ご回答よろしくお願いいたします。
- (C.13)の方は、ヒントより、$(\vec A\times\vec B)\cdot (\vec C\times \vec D)=\sum_i \sum_{j,k}\epsilon_{ijk}A_jB_k\sum_{\ell,m}\epsilon_{i\ell m}C_\ell D_m$に$\sum_{i} \epsilon_{ijk}\epsilon_{i\ell m}=\delta_{j\ell}\delta_{km}-\delta_{jm}\delta_{k\ell}$を代入して、 $$ \sum_{j,k,\ell,m}\left( \delta_{j\ell}A_j C_\ell\delta_{km}B_k D_m-\delta_{jm}A_j D_m\delta_{k\ell}B_k C_\ell\right)=(\vec A\cdot\vec C)(\vec B\cdot\vec D)-(\vec A\cdot\vec D)(\vec B\cdot\vec C) $$ となる。 -- 前野?
- すいません、答えは↑のようになります。ヒントの式も少し間違ってました。 -- 前野?
- 昨日の今日でご返信くださり、忘れない内に答え合わせができました。ありがとうございます。 -- kf?
p.151(6.47)式について †
SP? (2020-04-21 (火) 00:12:32)
(6.47)式の下の説明に「T_1とT_2を定数倍して、T^tMT=Eとなるようにできる。」とありますがどのように定数倍すればよろしいのでしょうか。
浅学ゆえの質問ですが、回答宜しくお願い致します。
- シンプルな話です。たとえば$(\vec T_1)^t {\bf M}\vec T_1=a$だったとすれば、$\vec T_1$を${1\over\sqrt{a}}$倍すれば（$\vec T_1\to {1\over\sqrt{a}}\vec T_1$と置き換えれば）、$(\vec T_1)^t {\bf M}\vec T_1=1$になります。 -- 前野?
- 同様に$(\vec T_2)^t{\bf M}\vec T_2$も1にすれば、行列は単位行列になります。 -- 前野?
- 理解しました。素早い返信ありがとうございました。 -- SP?
p.151(6.47)について †
SP? (2020-04-18 (土) 23:33:42)
(T1)t M (T2)=0ならば(6.47)が成り立つ理由がわかりません。
また、行列内の（T1)1から(T2)2の意味するところがわからないので教えていただけますか。
回答よろしくお願いいたします。
- $(\vec T_1)_1$の意味がわからなかったら(6.47)が成り立つ理由もわかるわけはないので、そっちから説明すると、$(\vec T_1)_1$はベクトル$\vec T_1$の第１成分です。$\vec T_1=\left(\begin{array}{cc}(\vec T_1)_1\\(\vec T_1)_2\end{array}\right)$と書いてもよい。 -- 前野?
- という意味を知った上で(6.47)の左辺を計算してみてください。非対角な部分に出てくるのは$(\vec T_2)^t{\bf M}\vec T_1$とその転置になって、0だとわかります。 -- 前野?
- 理解できました。 -- SP?
- 丁寧なご回答ありがとうございました。 -- SP?
P254(10.38)について †
やまだ? (2020-04-15 (水) 18:36:36)
「n回繰り返し」の項についている係数についての質問です。Nを無限大にした時、
nがNに対して十分に小さい有限な値なら確かに係数は(10.38)にある通りに収束します。
ですが、nは0からNの範囲を取りうるはずです。nがNに十分に近づいても本当にこの係数は保たれますか？
Nを無限大にしない時、係数は
$\prod_{i=0}^{n-1} \frac{N-i}{N}\frac{1}{n!}$
です。ここでnをNに対して十分に大きく(例えばNのおおよそ半分)すれば、
$\prod_{i=0}^{n-1}\frac{N-i}{N}$
はNを無限大にすると0に収束するはずです。
返答よろしくお願いします。
- Nは→∞の極限をとるという話しなので、有限個であるnは「Nの半分」にはできないです（どんなnに対しても、それを「小さい」と思えるような大きなNが取れる、というのが極限の考え方です）。 -- 前野?
- あと、${N-i\over N}$は（$N$は$i$に比べいくらでも大きくなれるのだから）1に収束すると考えます。 -- 前野?
- TeXの変換まずかったところも直しました。 -- 前野?
- わかりました。ありがとうございます。 -- やまだ?
p.39練習問題(問2-4)について †
med? (2020-04-15 (水) 14:29:34)
p.39練習問題(問2-4)の解答についての質問です。
I の変分の第一項のルートの中にδy'が含まれていないのはなぜなのでしょうか。
回答よろしくお願い致します。
- ここでは経路の変更を「ｙを変えずにｘ方向に動かす」という形にしているからです。ｘを変えずにｙ方向に動かすという変分もできますが、それは同じ「直線」という結果になります。 -- 前野?
- 納得しました。ありがとうございました。 -- med?
p363 問い10-5の解答について †
小泉? (2020-03-20 (金) 09:23:45)
１）D.98でPを一定にしてQでPを微分した場合、
∂P/∂Q=0ではないでしょうか？
私は次のようにしました。qを一定と仮定,(D,12)を利用
∂P/∂p=(∂P/∂Q)×(∂Q/∂p)=(∂P/∂Q)×(-∂q/∂P)=-∂q/∂Q
となります。
∂P/∂p=-∂q/∂Q、∂p/∂Q=-∂P/∂qを使うと｛Q,P}=0となってしまいます。
２）D.98, D.99が成立した場合でも、D.100に於いて
(∂Q/∂q)*(∂q/∂Q)=1, (∂P/∂q)*(∂q/∂Q)=1
なので、｛Q,P}=1+1=2 となるのではないでしょうか？
よろしくお願い致します。
- こういうとき、「何を変数としてどの変数を固定してどの変数で微分しているか」をちゃんと把握してないと間違えます。まず本文の(D.98)の上に「$P=P(q(Q,P),Q)$をPを一定としてQで微分して」というのは「P,Qを変数としてPを一定としてQで微分」です。 -- 前野?
- 「Pを微分」というのは「Pの変化量を考える」ことなので、「Pを一定として」という条件をつけて微分したら０に決まってます。だから(D.98)の左辺が０なのはこれでいいわけです。 -- 前野?
- (D.98)の右辺でやっているのは、$P(q(Q,P),Q)$の中にはQが２箇所あるので、その両方微分して、足しましょうということで、もちろんこうやって微分したって答えは０になるよね、というのが(D.98)という等式です。 -- 前野?
- ２）の、${\partial Q\over\partial q}*{\partial q\over\partial Q}=1$というのは間違いです。偏微分ではこういう単純な「約分」はできません。 -- 前野?
- 正確に書いておくと、${\partial Q(p,q)\over\partial q}\biggr|_p{\partial q(Q,P)\over\partial Q}\biggr|_P$（ポアッソン括弧に出てくるのはこの式）は、１ではありません。 ${\partial Q(p,q)\over\partial q}\biggr|_p{\partial q(p,Q)\over\partial Q}\biggr|_p$なら１です（固定する変数が共通で、ｐ）。-- 前野?
- ご説明いただきよく理解できました。(D.99)はqを言ってウニ -- 小泉?
- ご説明いただきよく理解できました。(D.99)ではqは一定なので、∂q(Q,P)/∂Q=0 -- 小泉?
- 当方のPCの調子が悪く同じことを記載して申し訳りません。上記から続けます。従って、D.100の第一項は０となり、第二項はQを一定にしているので１となるといった理解でよろしいでしょうか？ -- 小泉?
- いいえ違います。１＋０＝１になるというのは大間違いです。 -- 前野?
- ポアッソン括弧を真面目に書くと${\partial Q(q,p)\over\partial q}\bigr|_q{\partial P(q,p)\over\partial p}\bigr|_q-{\partial P(q,p)\over\partial q}\bigr|_q{\partial Q(q,p)\over\partial p}\bigr|_q$です。 -- 前野?
- (D.100)までで証明したことは${\partial Q(q,p)\over\partial p}\bigr|_q=-{\partial q(Q,P)\over\partial P}\bigr|_Q$と、${\partial P(q,p)\over\partial p}\bigr|_q={\partial q(Q,P)\over\partial Q}\bigr|_P$です。 -- 前野?
- この二つを入れると、ポアッソン括弧が${\partial Q(q,p)\over\partial q}\bigr|_p {\partial q(Q,P)\over\partial Q}\bigr|_P+{\partial P(q,p)\over\partial q}\bigr|_p {\partial q(Q,P)\over\partial P}\bigr|_Q$になりますが、これは $q(Q(q,p),P(q,p))$を$q$で微分したものだから、１です。 -- 前野?
- このあたりは練習問題なのとスペースの関係で省略記法を使ってますが、最初のうちは引数や固定する変数をいちいち書きながら計算することを勧めます（慣れてくれば省略記法で計算できるようになりますが）。 -- 前野?
- 大変丁寧にご説明いただき、ありがとうございます。よく理解できました。引数、固定する変数を書きながら計算するように致します。 -- 小泉?
p219 (9.42)について †
小泉? (2020-03-01 (日) 13:59:55)
H=ΣP_i(dq_i/dt)-Lであるが、(9.42)ではラグランジアンLの扱いが省略されており、何故 H(q+ε･･,p-ε･･)=H(q,p)となるのか、理解できません。ご教示をお願いします。
- (9.42)の１行目は$f(x+\epsilon b,y)=f(x,y)+\epsilon b {\partial f(x,y)\over\partial x}$のような偏微分の計算をやってます（$b={\partial H\over\partial p}$とかになっている）。 これはラグランジアンの形とかは関係なくできる計算です。 -- 前野?
- 二つめの＝は、${\partial H\over\partial p}{\partial H\over\partial q}-{\partial H\over\partial q}{\partial H\over\partial p}=0$を使って消しているだけで、これもラグランジアンの形とは関係なく成り立つ式です。 -- 前野?
- よくわかりました。テイラー展開が使われていることに気が付きませんでした。ご教示ありがとうございます。 -- 小泉?
p.267(10.83)について †
tatsu? (2020-02-27 (木) 16:32:24)
$L=\frac{1}{2}m\left(\dot{Q}+gt\right)^2$
という加速系のラグランジアンから求めた運動量は$P=m\left(\dot{Q}＋gt\right)$であり、
$K=P\dot{Q}-\frac{P^2}{2m}=\frac{P^2}{2m}-Pgt$ (10.83)
がハミルトニアンである。
という説明がありますが、$\left(q,p\right)\rightarrow \left(Q,P\right)$において、$\left\{Q,P\right\}_{\left(q,p\right)}=1$を確認しておりません。その後の母関数を使った変換で導いたハミルトニアン$K$と(10.83)のハミルトニアン$K$は一致するのですが、"(10.28)の段階"でのハミルトニアン$K$は正準変換によるハミルトニアンと言えるのでしょうか？
お忙しいところ申し訳ありませんがよろしくお願いいたします。
- ここでは、加速度系の作用から直接ハミルトニアンを求めてます。このハミルトニアンは、「作用から通常の手続きで作った」という意味で、正しいハミルトニアンです。 -- 前野?
- 要は「慣性系のラグランジアン」と「それを座標変換したラグランジアン」があって、それぞれを元にハミルトニアンを二つ作り、その二つが果たして正準変換でつながるかどうかをその先で調べるわけです。(10.83)のKが正しいハミルトニアンで、Q,Pのポアッソン括弧が正しい形なのは、「作用から作ったから」で保証されてます。 -- 前野?
- 分りました。ありがとうございます。 -- tatsu?
p.255の説明の件 †
tatsu? (2020-02-24 (月) 06:31:30)
p.255に以下のような説明があります。
$\left(q,p\right)\rightarrow\left(Q,P\right)$という正準変換を行ったとき、作用も
$\int\left(p\dot{q}-H\left(q,p\right)\right)dt\rightarrow\int\left(P\dot{Q}-H\left(q\left(Q,P\right),p\left(P,Q\right)\right)\right)dt$
と変化する。それでも正準方程式が変化しない為には、どんな条件が必要だろうか。・・・・・・・・・
つまりこの場合、正準方程式が変わらずにラグラジアンが変化する。そうなるのは
$\int\left(p\dot{q}-H\left(q,p\right)\right)dt=\int\left(P\dot{Q}-H\left(q\left(Q,P\right),p\left(P,Q\right)\right)+\frac{dG}{dt}\right)dt$
のように「表面項」になる量が付け加わった場合である。
上記の説明だと、$p\dot{q}=P\dot{Q}+\frac{dG}{dt}$という条件は正準方程式が変化しない条件というより、作用が変化しないための条件だと読み取れるのですが、$p\dot{q}=P\dot{Q}+\frac{dG}{dt}$という条件が正準方程式が変化しない条件である理由を教えて頂けませんでしょうか？
宜しくお願いいたします。
- ${\mathrm dG\over\mathrm dt}$の項が運動方程式（正準方程式）に効かない、ということはいいですね？ だとすると、$(p,q)$系の正準方程式は$\int (p\dot q-H(p,q))\mathrm dt$の変分から、$(P,Q)$系の正準方程式は$\int (P\dot Q-H(P,Q))\mathrm dt$の変分から得られます。 -- 前野?
- するとどっちも$\dot p=-{\partial H\over\partial q},\dot q={\partial H\over\partial p}$、$\dot P=-{\partial H\over\partial Q},\dot Q={\partial H\over\partial P}$の形の正準方程式になります。$p\dot q=P\dot Q$にG以外の余分なのがつくと、変分から導かれる方程式がこの形になりません。 -- 前野?
- 理解できました。ありがとうございます。 -- tatsu?
- すみません。まだ良く分っていなかったです。$\frac{dG}{dt}$を加えて作用を等しくする理由は何故でしょうか？ -- tatsu?
- 作用を等しくしないと正準方程式は同じ形にならないです。${\mathrm dG\over\mathrm dt}$は加えたいから加えるのではなくて「この形になっていれば加えても支障はない」ということです。 -- 前野?
- 一般的には座標変換すれば、$p\dot q\to P\dot Q+(?)$になりますが、$(?)$の部分が「悪い形」だと、新しい変数は正準方程式を満たしません。$(?)$の「悪くない形」が${\mathrm dG\over\mathrm dt}$です。 -- 前野?
- 分りました。お忙しいところありがとうございました。 -- tatsu?
P.240の(9.40)式について †
tatsu? (2020-02-20 (木) 00:10:50)
P.240の(9.40)式に$\frac{d\left(p_i\delta q_i\right)}{dt}$という式がありますが、これは$\sum_i\frac{d\left(p_i\delta q_i\right)}{dt}$ではないでしょうか？宜しくお願い致します。
- Texが上手く変換されなかったので、再度質問を書きます。P.240の(9.40)式に -- tatsu?
- Texが上手く変換されなかったので、再度質問を書きます。P.240の(9.40)式に $\frac{d\left(p_i\delta q_i\right)}{dt}$ がありますが、これは $\sum_i\frac{d\left(p_i\delta q_i\right)}{dt}$ ではないでしょうか？ よろしくお願いいたします。 -- tatsu?
- 度々間違えてすみません。p.240の(9.40)式ではなく、P.214の(9.40)式の$\frac{d\left(p_i\delta q_i\right)}{dt}$が\$sum_i\frac{d\left(p_i\delta q_i\right)}{dt}$ではないか？という質問です。 -- tatsu?
- 度々度々間違えてすみません。p.240の(9.40)式ではなく、"P.214"の(9.40)式の$\frac{d\left(p_i\delta q_i\right)}{dt}$が$\sum_i \frac{d\left(p_i\delta q_i\right)}{dt}$ではないか？という質問です。 -- tatsu?
- コメント遅れましたすみません。この式は間違ってますが、むしろ$\sum_i$の位置が最初のカッコの後ろにあるのが間違いで、最初のカッコより前にあって、全部にかかっているべきですね。 -- 前野?
- 訂正としてはもちろん、${\mathrm d(p_i\delta q_i)\over \mathrm dt}$の前に$\sum_i$をつけても構いません。 -- 前野?
- $-{\partial H\over\partial p_i}$の後ろに二つある）のうち一つを、${\mathrm d(p_i\delta q_i)\over\mathrm dt}$の後ろに持っていく、という修正でもいいですね。 -- 前野?
- 分りました。ありがとうございました。 -- tatsu?
ネーターの定理が時間並進性のとき成り立つ条件について †
tatsu? (2020-02-18 (火) 16:52:05)
p.202のネーターの定理に「ある変数変換 $q_i \rightarrow q_i+\delta q_i$を行ったとき、
$L\left(\left\{ｑ_*\right\},\left\{\dot{q}_i\right\}\right) \rightarrow L\left(\left\{ｑ_*\right\},\left\{\dot{q}_i\right\}\right)+\frac{\text{d}J}{\text{d}t}$
とあった時」とありますが、単振動の時のラグラジアン
$L=\frac{1}{2}m\dot{x}^{2}-\frac{1}{2}kx^{2}$
は、$x\rightarrow x-\epsilon\dot{x}$の時、$L=\frac{１}{２}ｍ\dot{x}^2-\frac{１}{２}kx^2+\epsilon kx\dot{x}-\frac{１}{２}\epsilon^2 k\dot{x}^2$ ・・・・(1)
となりますが、
$L\left(x,\dot{x}\right) \rightarrow L\left(x,\dot{x}\right)+\frac{\text{d}J}{\text{d}t}=L\left(x,\dot{x}\right)+\frac{\text{d}}{\text{d}t}\left(-\epsilon L\right)$
の条件は、$L(x,\dot{x})+\frac{d}{dt}(−ϵL)＝\frac{1}{2}m\dot{x}^2−\frac{1}{2}kx^2ーϵ(m\dot{x}\ddot{x}−kx\dot{x})$ ・・・・(2)
となり、(1)式と(2)式が一致しません。
ここで、(1)式において$\epsilon^2$の項は無視して、更に(2)式において$m\ddot{x}=-kx$を代入すると、(1)式は
$L=\frac{１}{２}ｍ\dot{x}^2-\frac{１}{２}kx^2+\epsilon kx\dot{x}$ ・・・(3)
(2)式は、
$L=\frac{１}{２}ｍ\dot{x}^2-\frac{１}{２}kx^2+\epsilon2kx\dot{x}$ ・・・(4)
となり、(3)式と(4)式は$\epsilon kx\dot{x}$の分だけ一致しません。
これは、一次元単振動のラグラジアン$L$が時間並進性を持たないと理解して宜しいでしょうか？宜しくお願い致します。
- 最後の行は「時間並進性をもたない」ではなく「時間の並進不変性をもたない」でした。すみません。 -- tatsu?
- 最後の行は「時間並進性をもたない」ではなく「時間の並進不変性をもたない」でした。すみません。 -- tatsu?
- 時間の並進不変性 -- tatsu?
- すみません。途中で送信してしまいました。時間の並進不変性があるラグラジアンを教えていただけないでしょうか。よろしくお願いいたします。 -- tatsu?
- (1)を出すとき、$\dot x$も$\dot x-\epsilon\ddot x$とずらす必要があります。 -- 前野?
- 分りました。$\epsilon^2$の項は無視するという理解で宜しいでしょうか？ -- tatsu?
- $\epsilon^2$は無視です。無視したくないなら、そもそも$x-\epsilon\dot x$という展開を１次で止めてはいけません。 -- 前野?
- 分りました。お忙しいところお時間を割いて頂きありがとうございました。 -- tatsu?
p151 λ_1 = λ_2 の場合扱い †
小泉? (2020-01-25 (土) 09:41:25)
λ_1 = λ_2の場合、最終的には（６．４９）の１行目のLの式を、対角行列に変換することが目的で論理を組み立てる必要があると考えます。シュミットの直交化と同様の方法を使うと欄外に説明がありますが、シュミットの直交化は互いに直交するベクトルをつくることが目的。どのように利用するかヒントをいただきたいと思います。
- $\lambda_1\neq\lambda_2$なら自動的に$\vec T_1^t{\mathbf M}\vec T_2=0$になりますが、$\lambda_1=\lambda_2$ではそうはいかないのが今の困ったところで、要は$\vec T_1^t{\mathbf M}\vec T_2=0$になるようにしてやればよい。 -- 前野?
- そのためにシュミットの直交化と似た方法を使います。具体的には$\vec T_3=\vec T_1 + a\vec T_2$のようなベクトルを作って、$\vec T_1^t{\mathbf M}\vec T_3=0$になるように定数を決める。 -- 前野?
- $\vec T_2$を使うのをやめて新しく作った$\vec T_3$を使う（規格化が必要ならやりなおす）ことにすれば後は同じようにできます。 -- 前野?
- 素直にaを計算し、a=-(T_11*T_11*m_1+T_12*T_12*m_2)/(T_11*T_21*m_1+ -- 小泉?
- T_12*T_22*m_2)となりT_3を求めてあとは教科書通りすすめることができました -- 小泉?
- ただ、aの分母がゼロになる場合を考えるとT_1ベクトル、T_2ベクトルが平行 -- 小泉?
- になることまでは理解できましたが、その場合の阻害に関してはこれから考えようと思います。 -- 小泉?
- ご教示ありがとうございます。 -- 小泉?
- システム更改で業務多忙であったため検討が遅れて申し訳ありませんでした。 -- 小泉?
- T_1､_2ベクトルは固有値ベクトルですので独立である必要があり、平行ではないので分母はゼロになってはいけないことに気が付きました。 -- 小泉?
p66 3.3 仮想仕事の原理を使う例題 変位δθの扱い †
小泉? (2020-01-03 (金) 09:51:48)
角θがδθ増える場合、手の仕事がFδθd(Lcosθ)/dθとあります。
これまでの議論から推論するに、δX=δθd(Lcosθ)/dθと考えることができると思いますが。なで、この等号＝が成り立つか理解できないので、ご教授いただきたく思います。よろしくお願いいたします。
- $X=L\cos\theta$なので、このθがδθ変化したら？という計算をしています。 -- 前野?
- $\delta X=L\cos(\theta+\delta\theta)-L\cos\theta$で１次までの展開をしているだけです。 -- 前野?
- よくわかりました。ありがとうございます。 -- 小泉?
これより古い記事は
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
