01何を探せばよいか

x¨+p x˙+q x=0(p, q は実数の定数)\ddot x + p\,\dot x + q\,x = 0\qquad (p,\ q \text{ は実数の定数})x¨+p x˙+q x=0(p, q は実数の定数)\begin{gathered}\ddot x + p\,\dot x + q\,x = 0\\[4pt](p,\ q \text{ は実数の定数})\end{gathered}

解を全部見つけるのは大変そうですが、この式は線形で右辺が 0 なので、独立な解を 2 つ見つければ、全部の解はその線形結合で書けます(線形微分方程式の重ね合わせ。結合の係数は複素数でもよい)。2 つが独立かは、Wronskian W=x1x˙2−x2x˙1W = x_1\dot x_2 - x_2\dot x_1(2×2 の行列式)が 0 でないことで確かめます。

02指数関数を試す

どんな関数を試すか。eλte^{\lambda t} は微分しても λ\lambda 倍になるだけなので、入れると全部の項が eλte^{\lambda t} でくくれます。

x¨+p x˙+q x=(λ2+p λ+q) eλt\ddot x + p\,\dot x + q\,x = \bigl(\lambda^{2} + p\,\lambda + q\bigr)\,e^{\lambda t}

eλte^{\lambda t} は 0 にならないので(λ\lambda が複素数でも ∣eλt∣=e(Re λ)t>0|e^{\lambda t}| = e^{(\mathrm{Re}\,\lambda) t} > 0。微分の公式 (eλt)′=λeλt(e^{\lambda t})' = \lambda e^{\lambda t} も複素数の λ\lambda で成り立つ。Euler の公式)、eλte^{\lambda t} が解になるのは λ\lambda が特性方程式の解のときで、そのときに限ります。微分方程式の解 x(t)x(t) と区別するため、特性方程式の解は根と呼びます。

特性方程式
λ2+p λ+q=0\lambda^{2} + p\,\lambda + q = 0
λ=−p±p2−4q2\lambda = \frac{-p \pm \sqrt{p^{2} - 4q}}{2}

特性方程式は λ\lambda の二次方程式なので、根は判別式 p2−4qp^2 - 4q の符号で 3 通り考えられます(二次方程式)。それぞれで解の形が決まります。

3 つの x(t)。上から、p² − 4q > 0 では 2 つの指数関数の和、p² − 4q = 0 では t をかけた指数関数、p² − 4q < 0 では揺れながら減る
根の種類ごとの解の形。どれも p > 0 の例

03異なる 2 つの実数の根 (p2−4q>0p^2 - 4q > 0)

2 つの根から、解 eλ1te^{\lambda_1 t} と eλ2te^{\lambda_2 t} が出ます。確かめるのは、この 2 つが独立かどうかです。λ1≠λ2\lambda_1 \neq \lambda_2 なら Wronskian は 0 にならず、独立です。

W=eλ1t⋅λ2eλ2t−eλ2t⋅λ1eλ1t=(λ2−λ1) e(λ1+λ2)t≠0\begin{aligned} W &= e^{\lambda_1 t}\cdot\lambda_2 e^{\lambda_2 t} - e^{\lambda_2 t}\cdot\lambda_1 e^{\lambda_1 t} \\ &= (\lambda_2 - \lambda_1)\,e^{(\lambda_1 + \lambda_2)t} \neq 0 \end{aligned}
x(t)=C1 eλ1t+C2 eλ2tx(t) = C_{1}\,e^{\lambda_{1}t} + C_{2}\,e^{\lambda_{2}t}
λ1,2=(−p±p2−4q)/2\lambda_{1,2} = \bigl(-p \pm \sqrt{p^2 - 4q}\bigr)/2

例題 1x¨+5x˙+6x=0\ddot x + 5\dot x + 6x = 0

λ2+5λ+6=(λ+2)(λ+3)=0,x=C1e−2t+C2e−3t\lambda^2 + 5\lambda + 6 = (\lambda + 2)(\lambda + 3) = 0,\qquad x = C_1 e^{-2t} + C_2 e^{-3t}λ2+5λ+6=(λ+2)(λ+3)=0,x=C1e−2t+C2e−3t\begin{gathered}\lambda^2 + 5\lambda + 6 = (\lambda + 2)(\lambda + 3) = 0,\\[4pt]x = C_1 e^{-2t} + C_2 e^{-3t}\end{gathered}

演習 1

x¨+3x˙+2x=0\ddot x + 3\dot x + 2x = 0 の一般解を求めよう。

答えを見る

特性方程式は

λ2+3λ+2=(λ+1)(λ+2)=0\lambda^{2} + 3\lambda + 2 = (\lambda + 1)(\lambda + 2) = 0

で、根は λ=−1, −2\lambda = -1,\ -2(異なる実数)。よって

x(t)=C1e−t+C2e−2t.x(t) = C_{1}e^{-t} + C_{2}e^{-2t}.

04重根 (p2−4q=0p^2 - 4q = 0)

根が λ0=−p/2\lambda_0 = -p/2 の 1 つだけだと、指数関数の解も 1 つしか出ず、2 つ目が足りません。そこで x=u(t) eλ0tx = u(t)\,e^{\lambda_0 t} とおいて、uu の式を調べます(どんな xx も u=x e−λ0tu = x\,e^{-\lambda_0 t} でこう書ける)。

x=u eλ0t とおくとx¨+p x˙+q x=[u¨+(2λ0+p) u˙+(λ02+p λ0+q) u] eλ0tx = u\,e^{\lambda_0 t}\ \text{とおくと}\qquad \begin{aligned} &\ddot x + p\,\dot x + q\,x \\ &= \Bigl[\ddot u + (2\lambda_0 + p)\,\dot u + (\lambda_0^{2} + p\,\lambda_0 + q)\,u\Bigr]\,e^{\lambda_0 t} \end{aligned}x=u eλ0t とおくとx¨+p x˙+q x=[u¨+(2λ0+p) u˙+(λ02+p λ0+q) u] eλ0t\begin{gathered}x = u\,e^{\lambda_0 t}\ \text{とおくと}\\[4pt]\begin{aligned} &\ddot x + p\,\dot x + q\,x \\ &= \Bigl[\ddot u + (2\lambda_0 + p)\,\dot u + (\lambda_0^{2} + p\,\lambda_0 + q)\,u\Bigr]\,e^{\lambda_0 t} \end{aligned}\end{gathered}

重根なので λ02+pλ0+q=0\lambda_0^2 + p\lambda_0 + q = 0 かつ 2λ0+p=02\lambda_0 + p = 0 で、残るのは u¨=0\ddot u = 0。2 回積分して u=C1+C2tu = C_1 + C_2 t なので

x(t)=(C1+C2 t) eλ0t,λ0=−p2x(t) = (C_{1} + C_{2}\,t)\,e^{\lambda_{0}t},\qquad \lambda_0 = -\frac{p}{2}x(t)=(C1+C2 t) eλ0t,λ0=−p2\begin{gathered}x(t) = (C_{1} + C_{2}\,t)\,e^{\lambda_{0}t},\\[4pt]\lambda_0 = -\frac{p}{2}\end{gathered}

すべての解についての置き換えなので、これで解は全部です。2 つ目の解 t eλ0tt\,e^{\lambda_0 t} は当てずっぽうではなく、こうして出てきます。

例題 2x¨+4x˙+4x=0\ddot x + 4\dot x + 4x = 0

λ2+4λ+4=(λ+2)2=0,x=(C1+C2t) e−2t\lambda^2 + 4\lambda + 4 = (\lambda + 2)^2 = 0,\qquad x = (C_1 + C_2 t)\,e^{-2t}λ2+4λ+4=(λ+2)2=0,x=(C1+C2t) e−2t\begin{gathered}\lambda^2 + 4\lambda + 4 = (\lambda + 2)^2 = 0,\\[4pt]x = (C_1 + C_2 t)\,e^{-2t}\end{gathered}

演習 2

x¨+2x˙+x=0\ddot x + 2\dot x + x = 0 の解で、x(0)=1, x˙(0)=0x(0) = 1,\ \dot x(0) = 0 を満たすものを求めよう。

答えを見る

特性方程式 λ2+2λ+1=(λ+1)2=0\lambda^{2} + 2\lambda + 1 = (\lambda + 1)^{2} = 0 は重根 λ0=−1\lambda_0 = -1 を持つので

x(t)=(C1+C2t) e−t,x˙(t)=(C2−C1−C2t) e−t.x(t) = (C_{1} + C_{2}t)\,e^{-t},\qquad \dot x(t) = (C_{2} - C_{1} - C_{2}t)\,e^{-t}.x(t)=(C1+C2t) e−t,x˙(t)=(C2−C1−C2t) e−t.\begin{gathered}x(t) = (C_{1} + C_{2}t)\,e^{-t},\\[4pt]\dot x(t) = (C_{2} - C_{1} - C_{2}t)\,e^{-t}.\end{gathered}

x(0)=C1=1x(0) = C_{1} = 1、x˙(0)=C2−C1=0\dot x(0) = C_{2} - C_{1} = 0 より C2=1C_{2} = 1。

x(t)=(1+t) e−t.x(t) = (1 + t)\,e^{-t}.

05複素数の根 (p2−4q<0p^2 - 4q < 0)

根が複素数だと、解 eλte^{\lambda t} も複素数になります。ほしいのは実数の解です。

λ=α±iβ,α=−p2,β=q−p24>0\lambda = \alpha \pm i\beta,\qquad \alpha = -\frac{p}{2},\quad \beta = \sqrt{q - \frac{p^{2}}{4}} > 0λ=α±iβ,α=−p2,β=q−p24>0\begin{gathered}\lambda = \alpha \pm i\beta,\\[4pt]\alpha = -\frac{p}{2},\quad \beta = \sqrt{q - \frac{p^{2}}{4}} > 0\end{gathered}

02 節のとおり、e(α±iβ)te^{(\alpha \pm i\beta)t} は複素数の値をとる解です。Euler の公式で書き直し、足して 2 で割る・引いて 2i2i で割ると、実数の解が 2 つ出ます(線形結合なので、やはり解)。

e(α±iβ)t=eαt(cos⁡βt±isin⁡βt)e^{(\alpha \pm i\beta)t} = e^{\alpha t}\bigl(\cos\beta t \pm i\sin\beta t\bigr)
e(α+iβ)t+e(α−iβ)t2=eαtcos⁡βt,e(α+iβ)t−e(α−iβ)t2i=eαtsin⁡βt\frac{e^{(\alpha+i\beta)t} + e^{(\alpha-i\beta)t}}{2} = e^{\alpha t}\cos\beta t,\qquad \frac{e^{(\alpha+i\beta)t} - e^{(\alpha-i\beta)t}}{2i} = e^{\alpha t}\sin\beta te(α+iβ)t+e(α−iβ)t2=eαtcos⁡βt,e(α+iβ)t−e(α−iβ)t2i=eαtsin⁡βt\begin{gathered}\frac{e^{(\alpha+i\beta)t} + e^{(\alpha-i\beta)t}}{2} = e^{\alpha t}\cos\beta t,\\[4pt]\frac{e^{(\alpha+i\beta)t} - e^{(\alpha-i\beta)t}}{2i} = e^{\alpha t}\sin\beta t\end{gathered}
W(eαtcos⁡βt, eαtsin⁡βt)=β e2αt≠0W(e^{\alpha t}\cos\beta t,\ e^{\alpha t}\sin\beta t) = \beta\,e^{2\alpha t} \neq 0
x(t)=eαt(C1cos⁡βt+C2sin⁡βt)x(t) = e^{\alpha t}\bigl(C_{1}\cos\beta t + C_{2}\sin\beta t\bigr)
λ=α±iβ\lambda = \alpha \pm i\betaα=−p/2\alpha = -p/2β=4q−p2/2\beta = \sqrt{4q - p^2}/2

x(0)x(0) と x˙(0)\dot x(0) が実数なら C1,C2C_1, C_2 も実数です。合成すると x=A eαtcos⁡(βt+ϕ)x = A\,e^{\alpha t}\cos(\beta t + \phi) とも書けます。

例題 3x¨+4x=0\ddot x + 4x = 0(単振動)

λ2+4=0,λ=±2i (α=0, β=2),x=C1cos⁡2t+C2sin⁡2t\lambda^2 + 4 = 0,\qquad \lambda = \pm 2i\ (\alpha = 0,\ \beta = 2),\qquad x = C_1\cos 2t + C_2\sin 2tλ2+4=0,λ=±2i (α=0, β=2),x=C1cos⁡2t+C2sin⁡2t\begin{gathered}\lambda^2 + 4 = 0,\\[4pt]\lambda = \pm 2i\ (\alpha = 0,\ \beta = 2),\\[4pt]x = C_1\cos 2t + C_2\sin 2t\end{gathered}

演習 3

x¨+2x˙+5x=0\ddot x + 2\dot x + 5x = 0 の一般解を、実数の関数で書こう。

答えを見る

特性方程式 λ2+2λ+5=0\lambda^{2} + 2\lambda + 5 = 0 の根は

λ=−2±4−202=−1±2i\lambda = \frac{-2 \pm \sqrt{4 - 20}}{2} = -1 \pm 2i

で、α=−1, β=2\alpha = -1,\ \beta = 2。よって

x(t)=e−t(C1cos⁡2t+C2sin⁡2t).x(t) = e^{-t}\bigl(C_{1}\cos 2t + C_{2}\sin 2t\bigr).