1-1-1. 速さと速度

Speed and Velocity

速さ

速さ

定義

時間 Δt\Delta t あたりの移動距離が Δx\Delta x であるような等速直線運動をする物体の速さ vv は次のように表す。

v=ΔxΔtv = \dfrac{\Delta x}{\Delta t}
説明

「速さ」とは「単位時間あたりの移動距離」のことである。 時間 Δt\Delta t の間に距離 Δx\Delta x だけ移動したとき,単位時間あたりに進む距離,すなわち速さ vv を求めることを考える。

このとき,次の比の式が成り立つ。

Δt:Δx=1:v\Delta t : \Delta x = 1 : v

数学公式により内項と外項の積は等しいから,

v⋅Δt=Δxv \cdot \Delta t = \Delta x

両辺を Δt\Delta t で割ることで,定義通りの式が導かれる。

v=ΔxΔtv = \dfrac{\Delta x}{\Delta t}

平均の速さ

公式

時刻 t1t_1 のとき移動距離が x1x_1,時刻 t2t_2 のとき x2x_2 である物体の平均の速さ v‾\overline{v} は次のように求める。

v‾=∣x2−x1t2−t1∣\overline{v} = \left| \dfrac{x_2 - x_1}{t_2 - t_1} \right|
導出

移動距離 Δx=∣x2−x1∣\Delta x=|x_2-x_1|,経過時間 Δt=∣t2−t1∣\Delta t=|t_2-t_1| より,速さの式に代入して

v‾=ΔxΔt=∣x2−x1t2−t1∣\overline{v} = \dfrac{\Delta x}{\Delta t} = \left| \dfrac{x_2 - x_1}{t_2 - t_1} \right|

瞬間の速さ

公式

物体の移動距離が xx のときの瞬間の速さ vv は,そのときの時刻 tt を用いて次のように表す。

v=∣dxdt∣v = \left| \dfrac{dx}{dt} \right|
導出

時刻 tt における物体の移動距離を tt の関数として x(t)x(t) と表すとき,時刻 t+h  (h>0)t+h \; (h>0) における移動距離は x(t+h)x(t+h) である。時刻 tt から t+ht+h までの物体の平均の速さ v‾\overline{v} は

v‾=∣x(t+h)−x(t)(t+h)−t∣\overline{v} = \left| \dfrac{x(t+h) - x(t)}{(t+h) - t} \right|

ここで時刻 tt,距離 xx における瞬間の速さを求めるには hh を限りなく 00 に近づけることを考えればよい。これは導関数の定義そのものに他ならない。したがって

v=∣lim⁡h→0x(t+h)−x(t)h∣=∣dxdt∣v = \left| \lim_{h \to 0} \dfrac{x(t+h) - x(t)}{h} \right| = \left| \dfrac{dx}{dt} \right|

運動の独立性

定義

多次元空間における物体の運動において,ある成分方向の運動状態は,その方向の力および初速度のみによって決定され,それと垂直な方向の運動の影響を受けない。

説明

この原理の根拠は,後述するニュートンの運動方程式 F=ma\boldsymbol{F} = m\boldsymbol{a} が「ベクトル方程式」であるという点にある。 例として2次元平面(xyxy 平面)での運動を考えると,ベクトル形式の運動方程式は次のように成分ごとに分解して記述できる。

{Fx=maxFy=may\begin{cases} F_x = m a_x \\ F_y = m a_y \end{cases}

ここで重要なのは,xx 方向の加速度 axa_x を決定するのは xx 方向の力 FxF_x のみであり,yy 方向の力 FyF_y がいかに変化しようとも axa_x には一切影響を与えないという数学的事実である。

この原理を象徴する有名な思考実験に「弾丸の同時落下」がある。水平に高速で発射された弾丸と,その瞬間に真下に自由落下させた弾丸は,空気抵抗を無視すれば,水平方向の速度の有無に関わらず,全く同時に地面に到達する。これは,鉛直方向の運動(落下)が水平方向の運動(前進)から完全に独立していることを示している。

物理学において,この独立性を見出し,問題を扱いやすい単位に分解して考える手法は,電磁気学における電場と磁場の重ね合わせや,量子力学における状態の重ね合わせなどへ進むための基礎的な考え方の枠組みとなっている。

速度ベクトル

公式

物体の位置ベクトルが r\boldsymbol{r} のときの速度ベクトル v\boldsymbol{v} は,そのときの時刻 tt を用いて次のように表す。

v=drdt\boldsymbol{v} = \dfrac{d\boldsymbol{r}}{dt}
導出

時間 Δt\Delta t だけ経過して物体の変位が Δr\Delta\boldsymbol{r} だったとき,速度ベクトルは ΔrΔt\dfrac{\Delta\boldsymbol{r}}{\Delta t} で求められ,時刻 tt における瞬間の速度ベクトルは Δt\Delta t を限りなく 00 に近づけることを考えればよい。これは導関数の定義そのものだから

v=lim⁡Δt→0ΔrΔt=drdt\boldsymbol{v} = \lim_{\Delta t \to 0} \dfrac{\Delta\boldsymbol{r}}{\Delta t} = \dfrac{d\boldsymbol{r}}{dt}

極座標系

2次元極座標系

公式

原点を O\mathrm{O} とする xyxy 平面において, xx 軸,yy 軸と平行な単位ベクトルをそれぞれ ex\boldsymbol{e}_x および ey\boldsymbol{e}_y とする。

また,xyxy 平面上の位置 (x,y)(x, y) にある物体の位置ベクトルを r\boldsymbol{r} ,速度ベクトルを v\boldsymbol{v} ,原点 O\mathrm{O} からの距離を rr ,物体と原点を結んだ線 と xx 軸がなす角を θ\theta とし,rr と θ\theta のうちどちらか一方のみが増す方向の単位ベクトルをそれぞれ er\boldsymbol{e}_r および eθ\boldsymbol{e}_\theta とする。

座標系変換

er\boldsymbol{e}_r および eθ\boldsymbol{e}_\theta は,xx 軸方向の単位ベクトル ex\boldsymbol{e}_x と yy 軸方向の単位ベクトル ey\boldsymbol{e}_y およびそのときの偏角 θ\theta を用いて次のように表す。

{er=cos⁡θex+sin⁡θeyeθ=−sin⁡θex+cos⁡θey\begin{cases} \boldsymbol{e}_r = \cos\theta \boldsymbol{e}_x + \sin\theta \boldsymbol{e}_y \\ \boldsymbol{e}_\theta = -\sin\theta \boldsymbol{e}_x + \cos\theta \boldsymbol{e}_y \end{cases}
速度ベクトル

時刻 tt で原点からの距離が rr である物体の速度ベクトル v\boldsymbol{v} は,そのときの偏角 θ\theta,中心から遠ざかる向きを表す単位ベクトル er\boldsymbol{e}_r と,中心のまわりを回転するときの向きを表す横断方向の単位ベクトル eθ\boldsymbol{e}_\theta を用いて次のように求める。

v=drdter+rdθdteθ\boldsymbol{v} = \dfrac{dr}{dt} \boldsymbol{e}_r + r \dfrac{d\theta}{dt} \boldsymbol{e}_\theta
導出
座標系変換

導出1

まず,原点から物体へ向かう向きを正とする径方向の単位ベクトル er\boldsymbol{e}_r を考える。このベクトルは xx 軸の正の向きと角度 θ\theta をなす長さ 11 のベクトルだから,その xx 成分および yy 成分は三角関数の定義に基づき,次のように表せる。

er=cos⁡θex+sin⁡θey⋯①\boldsymbol{e}_r = \cos\theta \boldsymbol{e}_x + \sin\theta \boldsymbol{e}_y \quad \cdots \text{①}

次に,er\boldsymbol{e}_r と垂直で,角度 θ\theta が増加する向きを正とする回転方向の単位ベクトル eθ\boldsymbol{e}_\theta を考える。このベクトルは er\boldsymbol{e}_r を反時計回りに π2\dfrac{\pi}{2} だけ回転させたものに相当するため,①式の角度 θ\theta を θ+π2\theta + \dfrac{\pi}{2} に置き換えることで得られ,

eθ=cos⁡(θ+π2)ex+sin⁡(θ+π2)ey\boldsymbol{e}_\theta = \cos\left(\theta + \dfrac{\pi}{2}\right) \boldsymbol{e}_x + \sin\left(\theta + \dfrac{\pi}{2}\right) \boldsymbol{e}_y

ここで,三角関数の性質(あるいは加法定理)である cos⁡(θ+π2)=−sin⁡θ\cos\left(\theta + \dfrac{\pi}{2}\right) = -\sin\theta であり,また sin⁡(θ+π2)=cos⁡θ\sin\left(\theta + \dfrac{\pi}{2}\right) = \cos\theta を用いて

eθ=−sin⁡θex+cos⁡θey\boldsymbol{e}_\theta = -\sin\theta \boldsymbol{e}_x + \cos\theta \boldsymbol{e}_y
導出2

位置ベクトル r\boldsymbol{r} を,位置 (x,y)(x, y) ,単位ベクトル ex\boldsymbol{e}_x ,ey\boldsymbol{e}_y および偏角 θ\theta を用いて表すと,長さ rr の xx 成分および yy 成分は三角関数の定義よりそれぞれ rcos⁡θr\cos\theta および rsin⁡θr\sin\theta だから

r=xex+yey=rcos⁡θex+rsin⁡θey\boldsymbol{r} = x \boldsymbol{e}_x + y \boldsymbol{e}_y = r\cos\theta \boldsymbol{e}_x + r\sin\theta \boldsymbol{e}_y

まず er\boldsymbol{e}_r について考えると,位置ベクトル r\boldsymbol{r} を rr で偏微分,すなわち距離 rr のみを微小に増減させて位置ベクトル r\boldsymbol{r} の変化率を求めることにより物体が動く方向を調べると

∂r∂r=cos⁡θex+sin⁡θey\dfrac{\partial \boldsymbol{r}}{\partial r} = \cos\theta \boldsymbol{e}_x + \sin\theta \boldsymbol{e}_y

単位ベクトル ex\boldsymbol{e}_x および ey\boldsymbol{e}_y の大きさはともに 11 であり,cos⁡2θ+sin⁡2θ=1\sqrt{\cos^2\theta + \sin^2\theta} = 1 であるため,右辺は大きさが 11 の rr のみを増す単位ベクトルすなわち er\boldsymbol{e}_r そのものである。

cos⁡θex+sin⁡θey=er\cos\theta \boldsymbol{e}_x + \sin\theta \boldsymbol{e}_y = \boldsymbol{e}_r

eθ\boldsymbol{e}_\theta についても同様に,位置ベクトル r\boldsymbol{r} を θ\theta で偏微分,すなわち偏角 θ\theta のみを微小に増減させて位置ベクトル r\boldsymbol{r} の変化率を求めることにより物体が動く方向を調べると

∂r∂θ=−rsin⁡θex+rcos⁡θey\dfrac{\partial \boldsymbol{r}}{\partial \theta} = -r\sin\theta \boldsymbol{e}_x + r\cos\theta \boldsymbol{e}_y

同じく単位ベクトル ex\boldsymbol{e}_x および ey\boldsymbol{e}_y の大きさはともに 11 であり, (−rsin⁡θ)2+(rcos⁡θ)2=r\sqrt{(-r\sin\theta)^2 + (r\cos\theta)^2} = r であるから,大きさが 11 の θ\theta のみを増す単位ベクトルすなわち eθ\boldsymbol{e}_\theta は

eθ=1r∂r∂θ=−sin⁡θex+cos⁡θey\boldsymbol{e}_\theta = \dfrac{1}{r} \dfrac{\partial \boldsymbol{r}}{\partial \theta} = -\sin\theta \boldsymbol{e}_x + \cos\theta \boldsymbol{e}_y
速度ベクトル

物体の位置ベクトルが r\boldsymbol{r} のときの速度ベクトル v\boldsymbol{v} は,そのときの時刻 tt を用いて 速度ベクトル

v=drdt\boldsymbol{v} = \dfrac{d\boldsymbol{r}}{dt}
導出1

位置ベクトル r\boldsymbol{r} を,位置 (x,y)(x, y) ,単位ベクトル ex\boldsymbol{e}_x ,ey\boldsymbol{e}_y および偏角 θ\theta を用いて表すと,長さ rr の xx 成分および yy 成分は三角関数の定義よりそれぞれ rcos⁡θr\cos\theta および rsin⁡θr\sin\theta だから

r=xex+yey=rcos⁡θex+rsin⁡θey\boldsymbol{r} = x \boldsymbol{e}_x + y \boldsymbol{e}_y = r\cos\theta \boldsymbol{e}_x + r\sin\theta \boldsymbol{e}_y

距離 rr と偏角 θ\theta はともに時間 tt とともに変化する関数であることに注意すると,積の微分に関する数学公式 {f(x)g(x)}′=f′(x)g(x)+f(x)g′(x)\{f(x)g(x)\}' = f'(x)g(x) + f(x)g'(x) および合成関数の微分に関する数学公式 {f(g(x))}′=f′(g(x))⋅g′(x)\{ f(g(x)) \}' = f'(g(x)) \cdot g'(x) を用いて

v=drdt=(drdtcos⁡θ+r(−sin⁡θ)dθdt)ex+(drdtsin⁡θ+rcos⁡θdθdt)ey=drdtcos⁡θex−rdθdtsin⁡θex+drdtsin⁡θey+rdθdtcos⁡θey=drdt(cos⁡θex+sin⁡θey)+rdθdt(−sin⁡θex+cos⁡θey)\begin{aligned} \boldsymbol{v} &= \dfrac{d\boldsymbol{r}}{dt} \\ &= \left( \dfrac{dr}{dt}\cos\theta + r(-\sin\theta)\dfrac{d\theta}{dt} \right)\boldsymbol{e}_x + \left( \dfrac{dr}{dt}\sin\theta + r\cos\theta\dfrac{d\theta}{dt} \right)\boldsymbol{e}_y \\ &= \dfrac{dr}{dt}\cos\theta \boldsymbol{e}_x - r\dfrac{d\theta}{dt}\sin\theta \boldsymbol{e}_x + \dfrac{dr}{dt}\sin\theta \boldsymbol{e}_y + r\dfrac{d\theta}{dt}\cos\theta \boldsymbol{e}_y\\ &= \dfrac{dr}{dt} ( \cos\theta \boldsymbol{e}_x + \sin\theta \boldsymbol{e}_y ) + r\dfrac{d\theta}{dt} ( -\sin\theta \boldsymbol{e}_x + \cos\theta \boldsymbol{e}_y ) \end{aligned}

ここで er=cos⁡θex+sin⁡θey\boldsymbol{e}_r = \cos\theta \boldsymbol{e}_x + \sin\theta \boldsymbol{e}_y および eθ=−sin⁡θex+cos⁡θey\boldsymbol{e}_\theta = -\sin\theta \boldsymbol{e}_x + \cos\theta \boldsymbol{e}_y より

v=drdter+rdθdteθ\boldsymbol{v} = \dfrac{dr}{dt} \boldsymbol{e}_r + r \dfrac{d\theta}{dt} \boldsymbol{e}_\theta
導出2

位置ベクトル r\boldsymbol{r} を大きさと単位ベクトルに分けて記述すると

r=rer\boldsymbol{r} = r \boldsymbol{e}_r

ここで,積の微分に関する数学公式 {f(x)g(x)}′=f′(x)g(x)+f(x)g′(x)\{f(x)g(x)\}' = f'(x)g(x) + f(x)g'(x) を用いて 両辺を tt で微分すると,デカルト座標の単位ベクトル ex\boldsymbol{e}_x と ey\boldsymbol{e}_y は時間に依存せず一定であっても極座標の単位ベクトル er\boldsymbol{e}_r は物体の位置(偏角 θ\theta)によって向きが変わるため,時間 tt の関数として微分の対象となることに注意して

v=drdt=drdter+rderdt⋯②\boldsymbol{v} = \dfrac{d\boldsymbol{r}}{dt} = \dfrac{dr}{dt} \boldsymbol{e}_r + r \dfrac{d\boldsymbol{e}_r}{dt}\quad \cdots \text{②}

ここで,第 22 項に含まれる単位ベクトルの時間微分 derdt\dfrac{d\boldsymbol{e}_r}{dt} を求める。er=cos⁡θex+sin⁡θey\boldsymbol{e}_r = \cos\theta \boldsymbol{e}_x + \sin\theta \boldsymbol{e}_y であり,θ\theta は tt の関数であるから,合成関数の微分法を用いると,

derdt=ddt(cos⁡θex+sin⁡θey)=−sin⁡θdθdtex+cos⁡θdθdtey=dθdt(−sin⁡θex+cos⁡θey)\begin{aligned} \dfrac{d\boldsymbol{e}_r}{dt} &= \dfrac{d}{dt} ( \cos\theta \boldsymbol{e}_x + \sin\theta \boldsymbol{e}_y ) \\ &= -\sin\theta \dfrac{d\theta}{dt} \boldsymbol{e}_x + \cos\theta \dfrac{d\theta}{dt} \boldsymbol{e}_y \\ &= \dfrac{d\theta}{dt} ( -\sin\theta \boldsymbol{e}_x + \cos\theta \boldsymbol{e}_y ) \end{aligned}

eθ=−sin⁡θex+cos⁡θey\boldsymbol{e}_\theta = -\sin\theta \boldsymbol{e}_x + \cos\theta \boldsymbol{e}_y より

derdt=dθdteθ\dfrac{d\boldsymbol{e}_r}{dt} = \dfrac{d\theta}{dt} \boldsymbol{e}_\theta

②式に代入して

v=drdter+rdθdteθ\boldsymbol{v} = \dfrac{dr}{dt} \boldsymbol{e}_r + r \dfrac{d\theta}{dt} \boldsymbol{e}_\theta
説明

一般に,物体の位置を特定するために必要な変数が最低 22 つ必要な世界を二次元といい,デカルト座標では xx と yy ,極座標では rr(距離)と θ\theta(角度)がその役割を担う。主として極座標系は円運動の単元で便利だが,単なる円運動のための特殊な道具ではなく,デカルト座標系( x,yx, y )と並んで平面上の運動を記述する普遍的な手法である。

なお,

速度ベクトル
の
導出2
で現れる derdt=dθdteθ\dfrac{d\boldsymbol{e}_r}{dt} = \dfrac{d\theta}{dt} \boldsymbol{e}_\theta という式は,径方向の単位ベクトル er\boldsymbol{e}_r が変化(微分)したとき,その変化の方向は元のベクトルに垂直な回転方向 eθ\boldsymbol{e}_\theta を向き,その変化の激しさは回転の速さ dθdt\dfrac{d\theta}{dt} に比例することを意味し,極座標の本質をよく表す。座標軸そのものが回転するという極座標の動的な性質は,今後円運動や惑星の運動に限らず,物体の運動をどの視点(座標系)で捉えるのが最も自然かという視点を持つために重要である。

3次元極座標系

公式

原点を O\mathrm{O} とする xyzxyz 空間において,xx ,yy ,zz 軸と平行な単位ベクトルをそれぞれ ex\boldsymbol{e}_x ,ey\boldsymbol{e}_y ,ez\boldsymbol{e}_z とする。

また,xyzxyz 空間上の位置 (x,y,z)(x, y, z) にある物体の位置ベクトルを r\boldsymbol{r},速度ベクトルを v\boldsymbol{v},原点 O\mathrm{O} からの距離を rr,原点と物体を結んだ線分と zz 軸がなす角(天頂角)を θ\theta,物体から xyxy 平面に下ろした垂線の足と原点を結んだ線分と xx 軸がなす角(方位角)を ϕ\phi とし,rr ,θ\theta, ϕ\phi のうちいずれか一方向のみが増す方向の単位ベクトルをそれぞれ er\boldsymbol{e}_r, eθ\boldsymbol{e}_\theta, eϕ\boldsymbol{e}_\phi とする。

座標系変換

er\boldsymbol{e}_r,eθ\boldsymbol{e}_\theta,eϕ\boldsymbol{e}_\phi は,xx軸,yy軸,zz軸それぞれに平行な単位ベクトル ex\boldsymbol{e}_x,ey\boldsymbol{e}_y,ez\boldsymbol{e}_z および天頂角 θ\theta,方位角 ϕ\phi を用いて次のように表す。

{er=sin⁡θcos⁡ϕex+sin⁡θsin⁡ϕey+cos⁡θezeθ=cos⁡θcos⁡ϕex+cos⁡θsin⁡ϕey−sin⁡θezeϕ=−sin⁡ϕex+cos⁡ϕey\begin{cases} \boldsymbol{e}_r = \sin\theta\cos\phi \boldsymbol{e}_x + \sin\theta\sin\phi \boldsymbol{e}_y + \cos\theta \boldsymbol{e}_z \\ \boldsymbol{e}_\theta = \cos\theta\cos\phi \boldsymbol{e}_x + \cos\theta\sin\phi \boldsymbol{e}_y - \sin\theta \boldsymbol{e}_z \\ \boldsymbol{e}_\phi = -\sin\phi \boldsymbol{e}_x + \cos\phi \boldsymbol{e}_y \end{cases}
速度ベクトル

時刻 tt における物体の速度ベクトル v\boldsymbol{v} は,中心から遠ざかる向きの単位ベクトル er\boldsymbol{e}_r,天頂角が増す向きの単位ベクトル eθ\boldsymbol{e}_\theta,方位角が増す向きの単位ベクトル eϕ\boldsymbol{e}_\phi を用いて次のように求める。

v=drdter+rdθdteθ+rsin⁡θdϕdteϕ\boldsymbol{v} = \dfrac{dr}{dt} \boldsymbol{e}_r + r \dfrac{d\theta}{dt} \boldsymbol{e}_\theta + r \sin\theta \dfrac{d\phi}{dt} \boldsymbol{e}_\phi
導出
座標系変換

導出1

まず,原点から物体へ向かう向きを正とする径方向の単位ベクトル er\boldsymbol{e}_r を考える。er\boldsymbol{e}_r を zz 軸および xyxy 平面へ射影すると,三角関数の定義に基づきその成分は z=cos⁡θz = \cos\theta,および xyxy 平面上の長さ sin⁡θ\sin\theta となり,後者を xx 軸,yy 軸へ射影すると,それぞれ sin⁡θcos⁡ϕ\sin\theta\cos\phi および sin⁡θsin⁡ϕ\sin\theta\sin\phi となるから

er=sin⁡θcos⁡ϕex+sin⁡θsin⁡ϕey+cos⁡θez⋯①\boldsymbol{e}_r = \sin\theta\cos\phi \boldsymbol{e}_x + \sin\theta\sin\phi \boldsymbol{e}_y + \cos\theta \boldsymbol{e}_z \quad \cdots \text{①}
ドラッグして視点を回転できます

径方向の単位ベクトル は赤色で示している。これを青色で示した 軸および 平面へ射影することにより、その成分は 、および 平面上の長さ となる。

この 平面上の成分をさらに 軸、 軸へ射影すると、それぞれ および となるため、

次に,θ\theta 方向の単位ベクトル eθ\boldsymbol{e}_\theta を考える。これは er\boldsymbol{e}_r と同じ垂直面内にあり,θ\theta を θ+π2\theta + \dfrac{\pi}{2} に置き換えたものに相当するため,①式の θ\theta を置換し,三角関数の性質を用いることで得られる。

eθ=sin⁡(θ+π2)cos⁡ϕex+sin⁡(θ+π2)sin⁡ϕey+cos⁡(θ+π2)ez=cos⁡θcos⁡ϕex+cos⁡θsin⁡ϕey−sin⁡θez\begin{aligned} \boldsymbol{e}_\theta &= \sin\left(\theta + \dfrac{\pi}{2}\right)\cos\phi \boldsymbol{e}_x + \sin\left(\theta + \dfrac{\pi}{2}\right)\sin\phi \boldsymbol{e}_y + \cos\left(\theta + \dfrac{\pi}{2}\right) \boldsymbol{e}_z \\ &= \cos\theta\cos\phi \boldsymbol{e}_x + \cos\theta\sin\phi \boldsymbol{e}_y - \sin\theta \boldsymbol{e}_z \end{aligned}
ドラッグして視点を回転できます

は、 と同じ垂直面内にあり、 方向へさらに直角に倒れたベクトルである。

最後に,ϕ\phi 方向の単位ベクトル eϕ\boldsymbol{e}_\phi は,xyxy 平面内において ϕ\phi が増加する向き(反時計回り)の単位ベクトルである。これは zz 成分を持たず,2次元極座標の eθ\boldsymbol{e}_\theta と同様の形式となる。

eϕ=−sin⁡ϕex+cos⁡ϕey\boldsymbol{e}_\phi = -\sin\phi \boldsymbol{e}_x + \cos\phi \boldsymbol{e}_y
ドラッグして視点を回転できます

は、点 P を通る 平面に平行な面内に存在し、水平面の動径方向に対して が増加する向きへ直角に回転したベクトルである。

導出2

位置ベクトル r\boldsymbol{r} をデカルト座標の成分表示で表すと,幾何学的な関係より

r=rsin⁡θcos⁡ϕex+rsin⁡θsin⁡ϕey+rcos⁡θez\boldsymbol{r} = r\sin\theta\cos\phi \boldsymbol{e}_x + r\sin\theta\sin\phi \boldsymbol{e}_y + r\cos\theta \boldsymbol{e}_z

この位置ベクトル r\boldsymbol{r} を,各変数 r,θ,ϕr, \theta, \phi で偏微分し,それぞれの変化率を調べる。

∂r∂r=sin⁡θcos⁡ϕex+sin⁡θsin⁡ϕey+cos⁡θez∂r∂θ=rcos⁡θcos⁡ϕex+rcos⁡θsin⁡ϕey−rsin⁡θez∂r∂ϕ=−rsin⁡θsin⁡ϕex+rsin⁡θcos⁡ϕey\begin{aligned} \dfrac{\partial \boldsymbol{r}}{\partial r} &= \sin\theta\cos\phi \boldsymbol{e}_x + \sin\theta\sin\phi \boldsymbol{e}_y + \cos\theta \boldsymbol{e}_z \\ \dfrac{\partial \boldsymbol{r}}{\partial \theta} &= r\cos\theta\cos\phi \boldsymbol{e}_x + r\cos\theta\sin\phi \boldsymbol{e}_y - r\sin\theta \boldsymbol{e}_z \\ \dfrac{\partial \boldsymbol{r}}{\partial \phi} &= -r\sin\theta\sin\phi \boldsymbol{e}_x + r\sin\theta\cos\phi \boldsymbol{e}_y \end{aligned}

これらのベクトルの大きさを計算すると,それぞれ 1,r,rsin⁡θ1, r, r\sin\theta となる。各偏微分ベクトルをそれぞれの大きさで割って正規化することで,各方向の単位ベクトルが導かれる。

er=∂r∂r=sin⁡θcos⁡ϕex+sin⁡θsin⁡ϕey+cos⁡θezeθ=1r∂r∂θ=cos⁡θcos⁡ϕex+cos⁡θsin⁡ϕey−sin⁡θezeϕ=1rsin⁡θ∂r∂ϕ=−sin⁡ϕex+cos⁡ϕey\begin{aligned} \boldsymbol{e}_r &= \dfrac{\partial \boldsymbol{r}}{\partial r} = \sin\theta\cos\phi \boldsymbol{e}_x + \sin\theta\sin\phi \boldsymbol{e}_y + \cos\theta \boldsymbol{e}_z \\ \boldsymbol{e}_\theta &= \dfrac{1}{r} \dfrac{\partial \boldsymbol{r}}{\partial \theta} = \cos\theta\cos\phi \boldsymbol{e}_x + \cos\theta\sin\phi \boldsymbol{e}_y - \sin\theta \boldsymbol{e}_z\\ \boldsymbol{e}_\phi &= \dfrac{1}{r\sin\theta}\dfrac{\partial \boldsymbol{r}}{\partial \phi} = -\sin\phi \boldsymbol{e}_x + \cos\phi \boldsymbol{e}_y \end{aligned}
速度ベクトル
v=drdt\boldsymbol{v} = \dfrac{d\boldsymbol{r}}{dt}
導出1

位置ベクトル r\boldsymbol{r} のデカルト成分 x=rsin⁡θcos⁡ϕ,y=rsin⁡θsin⁡ϕ,z=rcos⁡θx = r\sin\theta\cos\phi, y = r\sin\theta\sin\phi, z = r\cos\theta を,積の微分法および合成関数の微分法を用いてそれぞれ時間 tt で微分する。

vx=drdtsin⁡θcos⁡ϕ+rdθdtcos⁡θcos⁡ϕ−rdϕdtsin⁡θsin⁡ϕvy=drdtsin⁡θsin⁡ϕ+rdθdtcos⁡θsin⁡ϕ+rdϕdtsin⁡θcos⁡ϕvz=drdtcos⁡θ−rdθdtsin⁡θ\begin{aligned} v_x &= \dfrac{dr}{dt}\sin\theta\cos\phi + r\dfrac{d\theta}{dt}\cos\theta\cos\phi - r\dfrac{d\phi}{dt}\sin\theta\sin\phi \\ v_y &= \dfrac{dr}{dt}\sin\theta\sin\phi + r\dfrac{d\theta}{dt}\cos\theta\sin\phi + r\dfrac{d\phi}{dt}\sin\theta\cos\phi \\ v_z &= \dfrac{dr}{dt}\cos\theta - r\dfrac{d\theta}{dt}\sin\theta \end{aligned}

これらを v=vxex+vyey+vzez\boldsymbol{v} = v_x \boldsymbol{e}_x + v_y \boldsymbol{e}_y + v_z \boldsymbol{e}_z としてまとめ,drdt,rdθdt,rdϕdtsin⁡θ\dfrac{dr}{dt}, r\dfrac{d\theta}{dt}, r\dfrac{d\phi}{dt}\sin\theta の各項で括り直すと,

v=drdt(sin⁡θcos⁡ϕex+sin⁡θsin⁡ϕey+cos⁡θez)+rdθdt(cos⁡θcos⁡ϕex+cos⁡θsin⁡ϕey−sin⁡θez)+rdϕdtsin⁡θ(−sin⁡ϕex+cos⁡ϕey)\begin{aligned} \boldsymbol{v} &= \dfrac{dr}{dt} (\sin\theta\cos\phi \boldsymbol{e}_x + \sin\theta\sin\phi \boldsymbol{e}_y + \cos\theta \boldsymbol{e}_z) \\ &+ r\dfrac{d\theta}{dt} (\cos\theta\cos\phi \boldsymbol{e}_x + \cos\theta\sin\phi \boldsymbol{e}_y - \sin\theta \boldsymbol{e}_z) \\ &+ r\dfrac{d\phi}{dt}\sin\theta (-\sin\phi \boldsymbol{e}_x + \cos\phi \boldsymbol{e}_y) \end{aligned}

括弧内がそれぞれ er,eθ,eϕ\boldsymbol{e}_r, \boldsymbol{e}_\theta, \boldsymbol{e}_\phi と一致するため,目的の速度ベクトルが導かれる。

導出2

位置ベクトルを r=rer\boldsymbol{r} = r \boldsymbol{e}_r とし,両辺を時間 tt で微分する。

v=drdter+rderdt⋯②\boldsymbol{v} = \dfrac{dr}{dt} \boldsymbol{e}_r + r \dfrac{d\boldsymbol{e}_r}{dt} \quad \cdots \text{②}

ここで er\boldsymbol{e}_r の時間微分を求める。er\boldsymbol{e}_r は θ\theta と ϕ\phi の関数であるから,連鎖律(全微分)より

derdt=∂er∂θdθdt+∂er∂ϕdϕdt\dfrac{d\boldsymbol{e}_r}{dt} = \dfrac{\partial \boldsymbol{e}_r}{\partial \theta} \dfrac{d\theta}{dt} + \dfrac{\partial \boldsymbol{e}_r}{\partial \phi} \dfrac{d\phi}{dt}

位置ベクトルの導出2での計算結果より,∂er∂θ=eθ\dfrac{\partial \boldsymbol{e}_r}{\partial \theta} = \boldsymbol{e}_\theta,および ∂er∂ϕ=sin⁡θeϕ\dfrac{\partial \boldsymbol{e}_r}{\partial \phi} = \sin\theta \boldsymbol{e}_\phi であるから,

derdt=dθdteθ+sin⁡θdϕdteϕ\dfrac{d\boldsymbol{e}_r}{dt} = \dfrac{d\theta}{dt} \boldsymbol{e}_\theta + \sin\theta \dfrac{d\phi}{dt} \boldsymbol{e}_\phi

これを②式に代入すると,空間における速度ベクトルの式が導出される。

v=drdter+rdθdteθ+rsin⁡θdϕdteϕ\boldsymbol{v} = \dfrac{dr}{dt} \boldsymbol{e}_r + r \dfrac{d\theta}{dt} \boldsymbol{e}_\theta + r \sin\theta \dfrac{d\phi}{dt} \boldsymbol{e}_\phi
説明

三次元空間において,物体の位置を特定するには最低 33 つの変数が必要であり,デカルト座標( x,y,zx, y, z )に対し,球座標では r,θ,ϕr, \theta, \phi がその自由度を担う。

球座標における速度ベクトルにおいて,rsin⁡θdϕdtr\sin\theta \dfrac{d\phi}{dt} という項が現れるのは,方位角 ϕ\phi の回転半径が zz 軸からの距離 rsin⁡θr\sin\theta に依存するという幾何学的な事実を反映している。また,

速度ベクトル
の
導出2
で示した通り,er\boldsymbol{e}_r の時間微分が eθ\boldsymbol{e}_\theta と eϕ\boldsymbol{e}_\phi の和で表されることは,径方向の単位ベクトルの「向きの変化」が,天頂角方向の回転と方位角方向の回転の重ね合わせとして記述できることを示している。このような「動く座標系」の視点は,人工衛星の軌道解析や,量子力学における水素原子の電子軌道の記述など,中心力場が支配する広大な空間の物理を理解する上での必須の基盤となる。

基底の動特性と物理学の広がり デカルト座標系(x,y,zx, y, z)の最大の特徴は,どこに移動しても単位ベクトル ex,ey,ez\boldsymbol{e}_x, \boldsymbol{e}_y, \boldsymbol{e}_z の向きが変わらないことである。そのため,微分しても単位ベクトル自体は消えてしまい,成分の微分だけを考えればよかった。

しかし,極座標や球座標では「物体の位置によって単位ベクトルの向き自体が変わる」という性質がある。速度の中に現れる rdθdtr \dfrac{d\theta}{dt} や rsin⁡θdϕdtr \sin\theta \dfrac{d\phi}{dt} といった項は,この**「基底の回転」**から生じるものである。

この考え方は,人工衛星の軌道計算や,電子の原子核まわりの運動(量子力学)を記述する際に不可欠となる。特に3次元の速度式に含まれる rsin⁡θdϕdtr \sin\theta \dfrac{d\phi}{dt} という項は,物体が経線方向に動く際に「緯度(θ\theta)によって円の半径が異なる」という幾何学的な事実を美しく数式に反映している。

相対速度

相対速度

定義

速度 vAv_{\mathrm{A}} の物体Aから,速度 vBv_{\mathrm{B}} の物体Bを見たときの物体Bの相対速度 vABv_{\mathrm{AB}} は次のように求める。

vAB=vB−vAv_{\mathrm{AB}} = v_{\mathrm{B}} - v_{\mathrm{A}}
60
100

相対速度ベクトル

公式

速度 vA\boldsymbol{v}_{\mathrm{A}} の物体Aから,速度 vB\boldsymbol{v}_{\mathrm{B}} の物体Bを見たときの物体Bの相対速度ベクトル vAB\boldsymbol{v}_{\mathrm{AB}} は次のように求める。

vAB=vB−vA\boldsymbol{v}_{\mathrm{AB}} = \boldsymbol{v}_{\mathrm{B}} - \boldsymbol{v}_{\mathrm{A}}
物体Aの速度ベクトル
物体Bの速度ベクトル
導出

xyzxyz 空間において,xx,yy,zz 軸方向それぞれの相対速度を考えると,vA=(vAx,vAy,vAz)\boldsymbol{v}_{\mathrm{A}} = (v_{\mathrm{A}x}, v_{\mathrm{A}y}, v_{\mathrm{A}z}),vB=(vBx,vBy,vBz)\boldsymbol{v}_{\mathrm{B}} = (v_{\mathrm{B}x}, v_{\mathrm{B}y}, v_{\mathrm{B}z}) として

vB−vA=(vBxvByvBz)−(vAxvAyvAz)=(vABxvAByvABz)=vAB\boldsymbol{v}_{\mathrm{B}} - \boldsymbol{v}_{\mathrm{A}} = \begin{pmatrix} v_{\mathrm{B}x} \\ v_{\mathrm{B}y} \\ v_{\mathrm{B}z} \end{pmatrix} - \begin{pmatrix} v_{\mathrm{A}x} \\ v_{\mathrm{A}y} \\ v_{\mathrm{A}z} \end{pmatrix} = \begin{pmatrix} v_{\mathrm{AB}x} \\ v_{\mathrm{AB}y} \\ v_{\mathrm{AB}z} \end{pmatrix} = \boldsymbol{v}_{\mathrm{AB}}