1次元の調和振動子 - 倭算数理研究所

古典力学のいろいろな系で運動方程式を解いていくシリーズ（目次）。　今回は1次元の調和振動子 (harmonic oscillator) の系を解いていきます。　1次元の調和振動子は ${ k }$ を正の定数（ばね定数）として

　　 ${ \displaystyle\begin{align*} F = -kx \end{align*}}$

で与えられる復元力によって運動します。　このときニュートンの運動方程式は

　　 ${ \displaystyle\begin{align*} m\frac{d^2x}{dt^2} = -kx \qquad\cdots (*) \end{align*}}$

となります。

パラメータの書き換え

まずは運動方程式を扱いやすくするために、パラメータ（ ${ m }$ , ${ k }$ ）を書き換えておきましょう。　運動方程式 (*) の両辺を ${ m }$ で割って ${ \omega = \sqrt{\frac{k}{m}} }$ とおくと

　　 ${ \displaystyle\begin{align*} &\frac{d^2x}{dt^2} = -\omega^2 x \qquad \left(\omega = \sqrt{\frac{k}{m}}\right) \\ \therefore \; &\frac{d^2x}{dt^2} + \omega^2 x = 0 \qquad\cdots (1) \end{align*}}$

となります。

特殊解を求める

次は (1) 式に対して

　　 ${ \displaystyle\begin{align*} x(t) = e^{\lambda t} \end{align*}}$

の形の特殊解を求めましょう。　この ${ x }$ を (1) 式に代入して

　　 ${ \displaystyle\begin{align*} &(\lambda^2 + \omega^2) x = 0 \\ \therefore \; &\lambda = \pm i\omega \end{align*}}$

を得ます。　よって特殊解は

　　 ${ \displaystyle\begin{align*} x(t) = e^{\pm i\omega t} \end{align*}}$

となります。

一般解と初期条件を課した解

斉次線形微分方程式の一般解は特殊解の線型結合で書けたんだったけな？

　　 ${ \displaystyle\begin{align*} x(t) = C_1\,e^{i\omega t} + C_2\,e^{-i\omega t} \end{align*}}$

ただし ${ C_1,\,C_2 }$ は（積分）定数です。

初期条件として

　　 ${ \displaystyle\begin{align*} x(0) &= x_0, & \frac{dx}{dt}(0) &= v_0 \end{align*}}$

を課すと

　　 ${ \displaystyle\begin{align*} C_1 + C_2 &= x_0 ,& i\omega (C_1 - C_2) &= v_0 \end{align*}}$

より

　　 ${ \displaystyle\begin{align*} C_1 = \frac{1}{2}\left(x_0 + \frac{v_0}{i\omega}\right),\quad C_2 = \frac{1}{2}\left(x_0 - \frac{v_0}{i\omega}\right) \end{align*}}$

よって初期条件を課した解は最終的に

　　 ${ \displaystyle\begin{align*} x(t) &= \frac{1}{2} \left(x_0 + \frac{v_0}{i\omega}\right)e^{i\omega t} + \frac{1}{2} \left(x_0 - \frac{v_0}{i\omega}\right)e^{-i\omega t} \\ &= \frac{x_0}{2} \left(e^{i\omega t} + e^{-i\omega t}\right) + \frac{v_0}{2i\omega}\left(e^{i\omega t} - e^{-i\omega t}\right) \\ &= x_0\cos \omega t + \frac{v_0}{\omega}\sin\omega t \end{align*}}$

となります。　これを時間微分で微分すると速度 ${ v(t) }$ も計算できて

　　 ${ \displaystyle\begin{align*} v(t) &= \frac{dx}{dt}(t) \\ &= v_0\cos \omega t - x_0\omega\sin\omega t \end{align*}}$

以上をまとめると

　　 ${ \displaystyle\begin{align*} x(t) &= x_0\cos \omega t + \frac{v_0}{\omega}\sin\omega t \\ v(t) &= v_0\cos \omega t - x_0\omega\sin\omega t \end{align*}}$