古典力学/イントロダクション/練習問題のソースを表示

==「平均の速度と瞬間の速度」の例題 ==
[[古典力学/イントロダクション#平均の速度と瞬間の速度]]
=== 速度一定の運動 ===
<math>x(t) = Vt + x(0)</math> で表される物体の運動について、物体の位置 <math>x(t)</math> を時刻 <math>t</math> について微分し、物体の各時刻における瞬間の速度 <math>v(t)</math> を求めよ。ただし <math>V</math> および <math>x(0)</math> は定数とする。
;【解】
位置 <math>x(t)</math> の時間微分は次のように計算できる。
:<math>\begin{align}
\frac{dx(t)}{dt} &= \lim_{h \to 0} \frac{x(t+h) - x(t)}{(t+h) - t} \\
                 &= \lim_{h \to 0} \frac{\left(V(t+h) + x(0)\right) - \left(Vt + x(0)\right)}{h} \\
                 &= \lim_{h \to 0} \frac{Vh}{h} \\
                 &= \lim_{h \to 0} V \\
                 &= V.
\end{align}</math>
位置の時間微分はいま物体の瞬間の速度に等しいのだから、
:<math>v(t) = \frac{dx(t)}{dt}</math>
より瞬間の速度は
:<math>v(t) = V</math>
と求まる。ここで <math>V = 0</math> とすれば物体は <math>x(t) = x(0)</math> で（ある座標系に対して）静止していることになる。このとき直ちに位置の時間微分は 0 となることが示される。つまり、変数 <math>t</math> に対して定数として振る舞う項の <math>t</math> 微分は 0 である。

=== 加速度一定の運動 ===
<math>x(t) = \frac{1}{2}at^2 + Vt + x(0)</math> で表される物体の運動について、物体の位置 <math>x(t)</math> を時刻 <math>t</math> について微分し、物体の各時刻における瞬間の速度 <math>v(t)</math> を求めよ。ただし <math>a</math>、<math>V</math> および <math>x(0)</math> は定数とする。
;【解】
位置 <math>x(t)</math> の時間微分は次のように計算できる。
:<math>\begin{align}
\frac{dx(t)}{dt} &= \lim_{h \to 0} \frac{
                    \left(\frac{1}{2}a\left(t + h\right)^2 + V(t+h) + x(0)\right) 
                  - \left(\frac{1}{2}at^2 + Vt + x(0)\right)}{h} \\
                 &= \lim_{h \to 0} \frac{ ath + \frac{1}{2}h^2 + Vh}{h} \\
                 &= \lim_{h \to 0} \left(V + at + \frac{1}{2}h\right) \\
                 &= V + at.
\end{align}</math>
従って物体の瞬間の速度は
:<math>v(t) = V + at</math>
と求まる。特に <math>t=0</math> の場合について <math>v(0) = V</math> が成り立つ。時刻 <math>t=0</math> における物体の速度 <math>v(0)</math> はその物体の'''初速''' (initial velocity) と呼ばれる。この問題では、<math>V</math> は物体の初速 <math>v(0)</math> に等しいので、こちらも初速と呼ばれる。

物体の速度 <math>v(t)</math> の時間変化率を考えると、<math>v(t)</math> の時間微分は <math>\dot{v}(t) = a</math> となることが分かる。速度の時間微分は'''[[w:加速度|加速度]]''' (acceleration) と呼ばれる。先ほど示した関係は、この問題における物体の加速度が一定であることを示している。加速度が一定であるような状況はたとえば、物体を地表付近で投射したときの物体の運動について、空気抵抗などを無視できるような場合などが当てはまる。あるいは傾斜が一定の斜面を物体が滑る場合についても同様に、物体の受ける加速度が一定であると見なすことができる。

日常的な意味で「加速」とは速度が増すことを言い、速度が遅くなることは「減速」と言って区別するが、物理学においてはどちらの結果が生じる場合でも、速度の時間変化率は「加速度」と呼ばれる<ref>ただし分かりやすさを求める上で、速度を減じるような加速度を特別に'''減速度''' (deceleration) と呼ぶことはときどきある。</ref>。また、「負の加速度」が必ずしも物体の速度を遅くするわけではなく、同様に「正の加速度」が加わっても物体の速度が増すとは限らないことに注意。たとえば「負の速度」を持つ物体に「負の加速度」を加えれば、物体より速く負の方向へ動いていくことになる。このことは一次元的な運動から離れて平面上の自由な物体の運動を考えてみると分かりやすい。平面上の運動では、正負の二方向に限らずあらゆる方向に速度を持つことができ、同様に加速度も平面上の様々な方向へ加えられる。このとき、速度や加速度が正であったり負であったりということの区別はそれほど意味を持たない。

=== 時間のベキで表される運動 ===
<math>x(t) = L\left(\frac{t}{\tau}\right)^n</math> で表される物体の運動について、物体の位置 <math>x(t)</math> を時刻 <math>t</math> について微分し、物体の各時刻における瞬間の速度 <math>v(t)</math> を求めよ。ただし <math>L</math>、<math>\tau</math> は定数であり、<math>n</math> は整数とする。
;【解】
指数 <math>n</math> によって場合分けをする。<math>n=0</math> の場合、<math>a^0=1</math> より <math>x(t)=L</math> だから <math>v(t) = 0</math>（ただし[[w:冪乗|ベキ]]の計算について <math>0^0 = 1</math> と定義する。あるいは <math>n=0</math> の場合に <math>x(t)=L</math> となるように <math>x(t)</math> を定義する）。
<math>n=1,2</math> の場合は例題 (a), (b) の結果を参考にすればよい。<math>n \ge 1</math> の場合について、
:<math>\alpha^n - \beta^n = \left(\alpha - \beta\right)\left(\alpha^{n-1} + \beta\alpha^{n-2} + \beta^2\alpha^{n-3} + \dots + \beta^{n-2}\alpha  + \beta^{n-1} \right) = \left(\alpha - \beta\right)\sum_{k=0}^{n-1}\alpha^{n-1-k}\beta^k</math>
と[[w:因数分解|因数分解]]できることに注意すれば、位置 <math>x(t)</math> の時間微分は次のように計算できる。
:<math>\begin{align}
\frac{dx(t)}{dt} &= \lim_{h \to 0} \frac{L\left(\frac{t+h}{\tau}\right)^n - L \left(\frac{t}{\tau}\right)^n}{h} \\
                 &= \frac{L}{\tau^n} \lim_{h \to 0} \frac{\left(t + h\right)^n - t^n}{h} \\
                 &= \frac{L}{\tau^n} \lim_{h \to 0} \frac{1}{h}
                    \left\{\left((t + h) - t\right) \left((t + h)^{n-1} + t(t + h)^{n-2} + \dots + t^{n-1}\right) \right\} \\
                 &= \frac{L}{\tau^n} \lim_{h \to 0} \left\{\sum_{k=0}^{n-1}(t + h)^{n-1-k}t^k\right\} \\
                 &= n\frac{L}{\tau^n}t^{n-1}  \quad (n \ge 1).
\end{align}</math>
この関係が指数 <math>n</math> が非負の整数であるとき成り立つことはすぐに確かめられる。<math>n</math> が負である場合について <math>t^n = 1/t^{-n}</math> だから、
:<math>\begin{align}
\frac{dx(t)}{dt} &= \lim_{h \to 0} \frac{L\left(\frac{t+h}{\tau}\right)^n - L \left(\frac{t}{\tau}\right)^n}{h} \\
                 &= \frac{L}{\tau^n} \lim_{h \to 0} \frac{1}{h}\left(\frac{1}{\left(t + h\right)^{-n}} - \frac{1}{t^{-n}}\right) \\
                 &= \frac{L}{\tau^n} \lim_{h \to 0} \frac{1}{h}\frac{t^{-n} - \left(t + h\right)^{-n}}{t^{-n}\left(t + h\right)^{-n}} \\
                 &= \frac{L}{\tau^n} \lim_{h \to 0} \frac{-\sum_{k=0}^{-n-1}(t + h)^{-n-1-k}t^k}{t^{-n}\left(t + h\right)^{-n}} \\
                 &= -\frac{L}{\tau^n} \lim_{h \to 0} \left\{\sum_{k=0}^{-n-1}(t + h)^{-1-k}t^{k+n}\right\} \\
                 &= n\frac{L}{\tau^n}t^{n-1} \quad (n < 0)
\end{align}</math>
となる。指数 <math>n</math> の正負の違いを除けば先に求めたものと全く同じ形になっている。つまり、任意の整数 <math>n</math> について、
:<math>\frac{d}{dt}\left(L\left(\frac{t}{\tau}\right)^n\right) = n\frac{L}{\tau^n}t^{n-1}</math>
である。

=== 時間のベキの和で表される運動 ===
<math>x(t) = L\sum_{k=0}^{n} \left\{f(k)\left(\frac{t}{\tau}\right)^k\right\}</math> で表される物体の運動について、物体の位置 <math>x(t)</math> を時刻 <math>t</math> について微分し、物体の各時刻における瞬間の速度 <math>v(t)</math> を求めよ。ただし <math>f(k)</math> は変数 <math>k</math> の関数であり、また <math>n</math> は正の整数である。
;【解】
(c) で得られた結果をそのまま各項に対して適用すればよい。<math>t</math> に関して微分をする際、<math>t</math> によらない関数は定数と見ることができるから、
:<math>\frac{d}{dt}\left(f(k)\left(\frac{t}{\tau}\right)^k\right) = f(k)\frac{d}{dt}\left(\frac{t}{\tau}\right)^k</math>
である。また各項の和について、
:<math>\begin{align}
\frac{d}{dt}\left(f(k)\left(\frac{t}{\tau}\right)^k + f(r)\left(\frac{t}{\tau}\right)^r\right) 
 &= \lim_{h \to 0} \frac{1}{h}\left\{
      f(k)\left(\frac{t+h}{\tau}\right)^k + f(r)\left(\frac{t+h}{\tau}\right)^r 
    - f(k)\left(\frac{t}{\tau}\right)^k - f(r)\left(\frac{t}{\tau}\right)^r
    \right\} \\
 &= \lim_{h \to 0} \left\{
      \frac{f(k)\tau^{-k}\left( \left(t+h\right)^k - t^k \right)}{h} 
    + \frac{f(r)\tau^{-r}\left( \left(t+h\right)^r - t^r \right)}{h}
    \right\} \\
 &= kf(k)\tau^{-k}t^{k-1} + rf(r)\tau^{-r}t^{r-1}
\end{align}</math>
が成り立つ。よく見ればこれはそれぞれの項の微分になっているから、すべての項の和についての微分は、各項についての微分の和を計算すればよいことになる。以上のことから位置 <math>x(t)</math> の時間微分より物体の瞬間の速度が求まる。
:<math>\begin{align}
v(t) = \frac{dx(t)}{dt} &= \frac{d}{dt}\sum_{k=0}^{n} \left\{f(k)\left(\frac{t}{\tau}\right)^k\right\} \\
                        &= \sum_{k=0}^{n} \frac{d}{dt}\left\{f(k)\left(\frac{t}{\tau}\right)^k\right\} \\
                        &= \sum_{k=0}^{n} \left\{kf(k)\tau^{-k}t^{k-1}\right\}.
\end{align}</math>

=== ある位置へ物体が到達する時刻 ===
<math>x(t) = \frac{1}{2}at^2 + Vt + x(0)</math> で表される物体の運動について、物体の位置 <math>x(t)</math> が <math>x(t) = X</math> となる時刻 <math>t(X)</math> および物体の位置に関する微分 <math>\frac{dt(X)}{dX}</math> を求めよ。ただし <math>a</math>、<math>V</math> および <math>x(0)</math> は定数とする。
;【解】
<math>t = t(X)</math> のとき <math>x(t) = X</math> である。すなわち、
:<math>x(t(X)) = X.</math>
このとき、物体の位置を具体的に書き下せば、
:<math>\frac{1}{2}at^2(X) + Vt(X) + x(0) = X</math>
と書けるので、これを <math>t(X)</math> について解けば、<math>a</math> が 0 でない場合について、
:<math>t(X) = -\frac{V}{a} \pm \frac{\sqrt{V^2 + 2a\left(X - x(0)\right)}}{a}</math>
と求まる。物体の位置に関する微分 <math>\frac{dt(X)}{dX}</math> について、
:<math>\begin{align}
\frac{dt(X)}{dX} &= \lim_{h \to 0} \left\{\frac{1}{h}\left(
                    \left(-\frac{V}{a} \pm \frac{\sqrt{V^2 + 2a\left(X + h - x(0)\right)}}{a}\right)
                  - \left(-\frac{V}{a} \pm \frac{\sqrt{V^2 + 2a\left(X - x(0)\right)}}{a}\right)
                    \right)\right\} \\
                 &= \pm \lim_{h \to 0} \left\{\frac{1}{ah}\left(
                    \sqrt{V^2 + 2a\left(X + h - x(0)\right)}
                  - \sqrt{V^2 + 2a\left(X - x(0)\right)}
                    \right)\right\} \\
                 &= \pm \frac{2}{\sqrt{V^2 + 2a\left(X - x(0)\right)}} 
                    \lim_{h \to 0} \left\{ \frac{V^2 + 2a\left(X - x(0)\right)}{2ah}\left(
                    \sqrt{1 + \frac{2ah}{V^2 + 2a\left(X - x(0)\right)}} - 1
                    \right)\right\}
\end{align}</math>
と表すことができる。ここで
:<math>\varepsilon^2 = 1 + \frac{2ah}{\sqrt{V^2 + 2a\left(X - x(0)\right)}}</math> 
と置き換えれば、
:<math>\begin{align}
\frac{dt(X)}{dX} &= \pm \frac{2}{\sqrt{V^2 + 2a\left(X - x(0)\right)}} 
                    \lim_{\varepsilon \to 1} \left\{
                    \frac{1}{\varepsilon^2 - 1}\left(
                    \varepsilon - 1
                    \right)\right\} \\
                 &= \pm \frac{2}{\sqrt{V^2 + 2a\left(X - x(0)\right)}} 
                    \lim_{\varepsilon \to 1} 
                    \frac{1}{\varepsilon + 1} \\
                 &= \pm \frac{1}{\sqrt{V^2 + 2a\left(X - x(0)\right)}} 
\end{align}</math>
が得られる。ところで <math>X = x(t)</math> と再び置き換えれば、<math>2a\left(x(t) - x(0)\right) = a^2t^2 + 2aVt</math> だから、
:<math>\left.\frac{dt(X)}{dX}\right|_{X = x(t)} = \pm \frac{1}{|V + at|} = \pm\frac{1}{|v(t)|}</math>
となる。ここで例題 (b) の結果を使った。右辺は物体の瞬間速度の[[w:絶対値|絶対値]] <math>|v(t)|</math> について[[w:逆数|逆数]]をとったものになっている。物体の瞬間速度 <math>v(t)</math> が正のとき、位置の変化量 <math>dX</math> は正であり、他方、物体の瞬間速度 <math>v(t)</math> が負のとき、位置の変化量 <math>dX</math> は負である。時間の向きを正に取るなら、瞬間速度 <math>v(t)</math> の向きが正負のいずれの場合でも時間の変化量 <math>dt</math> は正である。従って、速度の絶対値を取らず
:<math>\left.\frac{dt(X)}{dX}\right|_{X = x(t)} = \frac{1}{v(t)}</math>
と表せる。瞬間速度が位置の時間微分であったことを思い出せば、上記の関係は以下のように整理することができる。
:<math>\left.\frac{dt(X)}{dX}\right|_{X = x(t)}\frac{dx(t)}{dt} = 1.</math>
ここまで <math>a \ne 0</math> の場合についてを解いたが、<math>a = 0</math> の場合はより単純で、
:<math>X = Vt(X) + x(0)</math>
より
:<math>t(X) = \frac{X - x(0)}{V}</math>
となる。また、物体の位置 <math>X</math> に関する微分は
:<math>\frac{dt(X)}{dX} = \frac{1}{V}</math>
となる。<math>a \ne 0</math> の場合と同様に、時刻を位置に関して微分したものは物体の瞬間速度の逆数として与えられる。

----
<references />

==「座標系の変換」の例題==
[[古典力学/イントロダクション#座標系の変換]]

=== 距離とデカルト座標系の性質 ===
平面上の物体の位置 <math>\boldsymbol{x}</math> を以下のように表す。
:<math>\boldsymbol{x} = x_1\mathbf{e}^{(1)} + x_2\mathbf{e}^{(2)}.</math>
任意の場合について、物体の原点からの距離を
:<math>|\boldsymbol{x}|^2 = \boldsymbol{x}\cdot\boldsymbol{x} = (x_1)^2 + (x_2)^2</math>
と位置の成分を用いて表すためには、[[w:基底 (線型代数学)|基底ベクトル]]の間にどのような関係が成り立てばよいか示せ。

;【解】
物体の原点までの距離を基底ベクトルの[[w:内積|内積]]を用いて表せば、
:<math>|\boldsymbol{x}|^2 =  x_1x_1\mathbf{e}^{(1)}\cdot\mathbf{e}^{(1)}
                          + 2x_1x_2\mathbf{e}^{(1)}\cdot\mathbf{e}^{(2)}
                          +  x_2x_2\mathbf{e}^{(2)}\cdot\mathbf{e}^{(2)}
</math>
となる。これが任意の場合について
:<math>(x_1)^2 + (x_2)^2 = (x_1)^2\mathbf{e}^{(1)}\cdot\mathbf{e}^{(1)}
                         + 2x_1x_2\mathbf{e}^{(1)}\cdot\mathbf{e}^{(2)}
                         + (x_2)^2\mathbf{e}^{(2)}\cdot\mathbf{e}^{(2)}
</math>
という関係が成り立つためには、基底の内積が以下の関係を満たす必要がある。
:<math>
\mathbf{e}^{(1)}\cdot\mathbf{e}^{(1)} = 1,\,
\mathbf{e}^{(1)}\cdot\mathbf{e}^{(2)} = 0,\,
\mathbf{e}^{(2)}\cdot\mathbf{e}^{(2)} = 1.
</math>
これらの関係は[[w:クロネッカーのデルタ|クロネッカーのデルタ]] (Kronecker delta) 
:<math>\delta_{ij} = 
\begin{cases}
  1, \quad i = j \\
  0, \quad \textrm{otherwise}
\end{cases}
</math>
を用いて次のように表現できる。
:<math>
\mathbf{e}^{(i)}\cdot\mathbf{e}^{(j)} = \delta_{ij}
</math>
このように内積がクロネッカーのデルタとなるような基底の組を[[w:正規直交基底|正規直交基底]] (orthonormal basis) という。直交基底 (orthogonal basis) であることは、他の基底ベクトルとの内積が 0 となること、言い換えればすべての基底ベクトルが互いに[[w:直交|直交]]していること、正規 (normal) であることは基底ベクトルの自分自身との内積が 1 に[[w:正規化|正規化]]（規格化）されていることを指す。

=== 基底ベクトルの回転 ===
平面上の物体の位置 <math>\boldsymbol{x}</math> を以下のように表す。
:<math>\boldsymbol{x} = x_1\mathbf{e}^{(1)} + x_2\mathbf{e}^{(2)}.</math>
また、[[w:行列|行列]] <math>\textrm{R}_\theta</math> を[[w:回転行列|回転行列]]とする。
座標系の回転 <math>\textrm{R}_\theta</math> によって位置の成分 <math>\boldsymbol{\xi}^\mathrm{T} = (x_1,\, x_2)</math> が
:<math>\boldsymbol{\xi}' = \textrm{R}_\theta\boldsymbol{\xi}</math>
と変換される場合について、
:<math>\mathbf{e}' = \textrm{P}\mathbf{e}</math>
と表されるような、座標系の基底 <math>\{\mathbf{e}^{(i)}\}_{i=\{1,2\}}</math> に対する変換行列 <math>\textrm{P}</math> を求めよ。
ただし、<math>\mathbf{e} = (\mathbf{e}^{(1)},\,\mathbf{e}^{(2)})^\mathrm{T}</math> は基底ベクトルの組を表し、<math>\mathbf{e},\,\mathbf{e}'</math> はそれぞれ回転変換前と変換後の基底ベクトルの組であるとする。

;【解】
座標系の回転に対して位置ベクトルは不変だから、
:<math>
\boldsymbol{x} = \boldsymbol{\xi}^\mathrm{T}\mathbf{e} 
               = \boldsymbol{\xi}'{}^\mathrm{T}\mathbf{e}'
               = (\{\textrm{P}^{-1}\}^\mathrm{T}\boldsymbol{\xi})^\mathrm{T}\textrm{P}\mathbf{e}
</math>
が成り立つ。ここで、成分ベクトルに対する変換が回転行列で表されることから
:<math>\textrm{R}_\theta\boldsymbol{\xi} = \{\textrm{P}^{-1}\}^\mathrm{T}\boldsymbol{\xi}</math>
より次の関係が成り立つ。
:<math>\textrm{P}^{-1} = {\textrm{R}_\theta}^\mathrm{T}.</math>
回転行列が[[w:直交行列|直交行列]]であることを利用すれば、
:<math>
{\textrm{R}_\theta}^\mathrm{T}\textrm{R}_\theta = \textrm{I} 
~\Rightarrow~ {\textrm{R}_\theta}^\mathrm{T} = {\textrm{R}_\theta}^{-1}
</math>
より、基底の組に対する変換行列は
:<math>\textrm{P} = \textrm{R}_\theta</math>
となることが示される。

----
<references />

{{DEFAULTSORT:こてんりきかく いんとろたくしよん れんしゆうもんたい}}
[[Category:古典力学|いんとろたくしよん れんしゆうもんたい]]