単振動 : 微分方程式の解法 (solution of differential equation)

角振動数 ω の単振動の従う微分方程式

d 2 x d t 2 + ω 2 x = 0     - - - (1)

の一般解を求める:  解法1  解法2  (初期値問題は ⇒)


解法1

与式は定数係数の2階同次線形微分方程式であるので,解を x = e λ t とおくと, d x / d t = λ e λ t , d 2 x / d t 2 = λ 2 e λ t より

d 2 x d t 2 + ω 2 x = λ 2 e λ t + ω 2 e λ t = ( λ 2 + ω 2 ) e λ t = 0

となり, e λ t ≠ 0 から,特性方程式

λ 2 + ω 2 = 0     - - - (2)

を得る.この特性方程式の解は λ = ± i ω となり,それぞれ λ 1 = i ω , λ 2 = − i ω とすると,式(1)の一般解は2つの独立な解 e λ 1 t , e λ 2 t  の線形結合

x = c 1 e λ 1 t + c 2 e λ 2 t = c 1 e i ω t + c 2 e − i ω t    ( c 1   ,   c 2 : 任意定数(複素数))

として求まる.オイラーの公式 e ± i θ = cos θ ± i sin θ を用いると,一般解は

x = c 1 ( cos ω t + i sin ω t ) + c 2 ( cos ω t − i sin ω t ) = ( c 1 + c 2 ) cos ω t + i ( c 1 − c 2 ) sin ω t

と書ける.物理量 x が実数であることを考えると, c 1 = ( A 1 + i A 2 ) / 2 , c 2 = ( A 1 − i A 2 ) / 2 とおいて

x = A 1 cos ω t + A 2 sin ω t    ( A 1   ,   A 2 : 任意定数(実数))     - - - (3)

が求まる(2つの独立な解として, cos ω t , sin ω t  を選び,それらの線形結合をとることに対応).

また, A = A 1 2 + A 2 2 , cos α = A 1 / A , sin α = − A 2 / A とおいて,加法定理を用いると,一般解は

x = A ( A 1 A cos ω t + A 2 A sin ω t ) = A ( cos ω t cos α − sin ω t sin α ) = A cos ( ω t + α )     - - - (4)

と書ける.


解法2  ページトップ

与式に d x / d t をかけて整理すると,

d x d t ( d 2 x d t 2 + ω 2 x ) = d x d t d 2 x d t 2 + ω 2 x d x d t = 1 2 d d t {   ( d x d t ) 2   } + 1 2 ω 2 d d t { x 2 } = 0

となるので,両辺を時刻 t で積分すると

1 2 ∫ d d t {   ( d x d t ) 2   } d t + 1 2 ω 2 ∫ d d t { x 2 } d t = 1 2 ( d x d t ) 2 + 1 2 ω 2 x 2 = C   (定数)

が得られる(この式は力学的エネルギー保存則を表している).

左辺はすべて正の項なので右辺の定数も正であり, C = ( 1 / 2 ) ω 2 A 2 とおく( A :定数).よって,

1 2 ( d x d t ) 2 + 1 2 ω 2 x 2 = 1 2 ω 2 A 2     ⇒     ( d x d t ) 2 = ω 2 ( A 2 − x 2 )     ⇒     d x d t = ± ω A 2 − x 2

となり,変数分離して両辺を積分すると

± d x A 2 − x 2 = ω d t     ⇒     ∫ ± d x A 2 − x 2 = ∫ ω d t     ⇒     ∓ cos − 1 x A = ω t + α   ( α :定数)

が得られ,両辺の cos をとって A をかければ,与式の一般解として

x = A cos ( ω t + α )     - - - (5)

が求まる.ここで, cos ( − θ ) = cos θ を用いた.



初期値問題  ページトップ

初期条件 x ( 0 ) = x 0 , v ( 0 ) = v 0 を満たす特殊解を求める( v ( t ) = d x / d t ).


一般解を  x = A 1 cos ω t + A 2 sin ω t  とした場合

v = d x d t = − A 1 ω sin ω t + A 2 ω cos ω t

初期条件より

x ( 0 ) = A 1 cos 0 + A 2 sin 0 = A 1 = x 0

v ( 0 ) = − A 1 ω sin 0 + A 2 ω cos 0 = A 2 ω = v 0     ⇒     A 2 = v 0 ω

なので,初期条件を満たす解は次式となる.

x = x 0 cos ω t + v 0 ω sin ω t     - - - (6)


一般解を  x = A cos ( ω t + α )  とした場合

v = d x d t = − A ω sin ( ω t + α )

初期条件より

x ( 0 ) = A cos ( 0 + α ) = A cos α = x 0

v ( 0 ) = − A ω sin ( 0 + α ) = − A ω sin α = v 0     ⇒     A sin α = − v 0 ω

なので,任意定数 A , α は

A 2 = ( A cos α ) 2 + ( A sin α ) 2 = x 0 2 + ( − v 0 ω ) 2 = x 0 2 + v 0 2 ω 2     - - - (7)

tan α = sin α cos α = − v 0 ω x 0 = − v 0 x 0 ω      ( x 0 ≠ 0 )     - - - (8)

を満たすように決定する( x 0 = 0 のときは, cos α = 0 を満たす α を考えればよい ).

この場合, A については正負の, α については n π ( n :整数)の任意性が残っているが, A > 0 と制限すると,

A = x 0 2 + v 0 2 ω 2     - - - (9)

であり, α については,次の2式

cos α = x 0 A = x 0 ω ( x 0 ω ) 2 + v 0 2   ,   sin α = − v 0 A ω = − v 0 ( x 0 ω ) 2 + v 0 2     - - - (10)

を同時に満たすように − π < α ≤ π の範囲内で考えれば,一意的に決まる.


ホーム>>カテゴリー分類>>力学>>質点の力学>>単振動>>微分方程式の解法

最終更新日:2026年9月25日