7. 定数係数2階非斉次線形微分方程式 例題集

$Q1$.
次の微分方程式の一般解を求めなさい。

$\dfrac{d^2x}{dt^2} + 6\dfrac{dx}{dt} + 8x = -3t$
解答・解説を見る
$x = C_1e^{-4t}+C_2e^{-2t} - \dfrac{3}{8}t + \dfrac{9}{32}~~$ ($C_1$, $C_2$ は任意定数)

まず, 斉次線形微分方程式

$\dfrac{d^2x}{dt^2} + 6\dfrac{dx}{dt} + 8x = 0$

を解いていきます。特性方程式

$\lambda^2 + 6\lambda + 8=0$

を考えると

$\lambda^2 + 6\lambda + 8=(\lambda + 4)(\lambda + 2)=0$

$\lambda = -4,~-2$ より, 斉次方程式の一般解は $x = C_1e^{-4t}+ C_2e^{-2t}$ になります。

次に, 非斉次方程式の右辺が $t$ の $1$ 次式であることから, 特殊解を

$x = At + B$

として非斉次方程式に代入すると

$6A + 8(At+B) = 8At + (6A+8B) = -3t$

係数を比べると $A = -\dfrac{3}{8}$, $B = \dfrac{9}{32}$ となるので

$x = -\dfrac{3}{8}t + \dfrac{9}{32}$

が特殊解であることがわかります。

非斉次方程式の一般解は特殊解と斉次方程式の一般解の和になるので, 一般解は

$x = C_1e^{-4t}+C_2e^{-2t} - \dfrac{3}{8}t + \dfrac{9}{32} ~~$ ($C_1$, $C_2$ は任意定数)

となります。

※別解(特殊解の求め方)

特殊解の公式

$\displaystyle x = -x_1\int \dfrac{x_2R(t)}{W(x_1,x_2)} ~dt + x_2\int \dfrac{x_1R(t)}{W(x_1,x_2)}~dt$

を用いると $x_1 = e^{-4t}$, $x_2 =e^{-2t}$, $R(t) = -3t$ であり

$W(e^{-4t},e^{-2t}) = \left| e4te2t4e4t2e2t \right| = 2e^{-6t}$

であるから不定積分の部分を計算すると

$x2R(t)W(x1,x2) dt=3te2t2e6t dt=32te4t dt=32(14te4t14e4t dt)=332e4t38te4t$

同様に

$x1R(t)W(x1,x2) dt=3te4t2e6t dt=32te2t dt=32(12te2t12e2t dt)=38e2t34te2t$

よって特殊解は

$x=x1x2R(t)W(x1,x2) dt+x2x1R(t)W(x1,x2) dt=e4t(332e4t38te4t)+e2t(38e2t34te2t)=(33238t)+(3834t)=(3834)t+(38332)=38t+932$

となります。

$Q2$.
次の微分方程式の一般解を求めなさい。

$\dfrac{d^2x}{dt^2} + \dfrac{dx}{dt} = 2t$
解答・解説を見る
$x = C_1 + C_2e^{-t} + t^2 -2t ~~$ ($C_1$, $C_2$ は任意定数)

$\dfrac{d^2x}{dt^2} + \dfrac{dx}{dt} = 0$ とすると, 特性方程式は

$\lambda^2 + \lambda = \lambda(\lambda + 1)=0$

$\lambda = 0,-1$ より, 斉次方程式の一般解は $x = C_1 + C_2e^{-t}$ となります。

次に, 非斉次方程式の右辺が $t$ の $1$ 次式であり, 左辺に $x$ の項が含まれないことから特殊解が

$x = At^2 + Bt$

という形をしていると仮定します。

ここで, 定数項は斉次方程式の一般解になっているので, 特殊解では考える必要がないことに注意しましょう。

非斉次方程式に代入すると

$2A + (2At + B) = 2At + (2A + B) = 2t$

係数を比べると $A = 1$, $B = -2$ となるので

$x = t^2 - 2t$

が特殊解となることがわかります。

よって, 非斉次方程式の一般解は

$x = C_1 + C_2e^{-t} + t^2 -2t ~~$ ($C_1$, $C_2$ は任意定数)

となります。

※注意

ここで特殊解を

$x = At + B$

と仮定してしまうと, 非斉次方程式に代入した時に

$A = 2t$

となってしまい, 係数を求めることができません。

※別解(特殊解の求め方)

特殊解の公式を使うと $x_1 = 1$, $x_2 = e^{-t}$, $R(t) = 2t$ であり

$W(1,e^{-t}) = \left| 1et0et \right| = -e^{-t}$

であるから不定積分の部分を計算すると

$x2R(t)W(x1,x2) dt=2tetet dt=2t dt=t2$

同様に

$x1R(t)W(x1,x2) dt=2tet dt=2tet dt=2(tetet dt)=2et2tet$

よって特殊解は

$x=x1x2R(t)W(x1,x2) dt+x2x1R(t)W(x1,x2) dt=(t2)+et(2et2tet)=t22t+2$

となります。定数項の $2$ は一般解に含まれることにも注意しましょう。

$Q3$.
次の微分方程式の一般解を求めなさい。

$\dfrac{d^2x}{dt^2} + 5\dfrac{dx}{dt} + 6x = 4e^{3t}$
解答・解説を見る
$x = C_1e^{-3t} + C_2e^{-2t} + \dfrac{2}{15}e^{3t}~~$ ($C_1$, $C_2$ は任意定数)

$\dfrac{d^2x}{dt^2} + 5\dfrac{dx}{dt} + 6x=0$ とすると, 特性方程式は

$\lambda^2 + 5\lambda + 6 = (\lambda + 3)(\lambda + 2) =0$

$\lambda = -3,-2$ より, 斉次方程式の一般解は $x = C_1e^{-3t}+ C_2e^{-2t}$ となります。

次に, 非斉次方程式の右辺が指数関数なので, 特殊解を

$x = Ae^{3t}$

と仮定して, 非斉次方程式に代入すると

$(9A + 15A + 6A)e^{3t} = 4e^{3t}$

係数を比べると $30A = 4$ より $A= \dfrac{2}{15}$ となります。よって

$x = \dfrac{2}{15}e^{3t}$

が特殊解であることがわかります。以上から, 非斉次方程式の一般解は

$x = C_1e^{-3t} + C_2e^{-2t} + \dfrac{2}{15}e^{3t}~~$ ($C_1$, $C_2$ は任意定数)

となります。

※別解(特殊解の求め方)

ロンスキアンを計算すると

$W(e^{-3t},e^{-2t}) = \left| e3te2t3e3t2e2t \right| = e^{-5t}$

不定積分の部分を計算すると

$x2R(t)W(x1,x2) dt=4ete5t dt=4e6t dt=23e6t$

同様に

$x1R(t)W(x1,x2) dt=4e5t dt=4e5t=45e5t$

よって特殊解は

$x=x1x2R(t)W(x1,x2) dt+x2x1R(t)W(x1,x2) dt=23e3t+45e3t=215e3t$

となります。

$Q4$.
次の微分方程式の一般解を求めなさい。

$\dfrac{d^2x}{dt^2} - \dfrac{dx}{dt} - 2x = e^{2t}$
解答・解説を見る
$x = C_1e^{-t} + C_2e^{2t} + \dfrac{1}{3}te^{2t}~~$ ($C_1$, $C_2$ は任意定数)

$\dfrac{d^2x}{dt^2} - \dfrac{dx}{dt} - 2x=0$ とすると, 特性方程式は

$\lambda^2 - \lambda - 2 = (\lambda + 1)(\lambda - 2) =0$

$\lambda = -1,~2$ より, 斉次方程式の一般解は $x = C_1e^{-t}+ C_2e^{2t}$ となります。

次に, 非斉次方程式の右辺が指数関数で, $x = e^{2t}$ が一般解に含まれていることから, 特殊解を

$x = Ate^{2t}$

と仮定して, 非斉次方程式に代入すると

$(4Ae^{2t} + 4Ate^{2t}) -(Ae^{2t} + 2Ate^{2t}) -2Ate^{2t} = e^{2t}$

整理すると $3Ae^{2t} = e^{2t}$ となるので $A=\dfrac{1}{3}$ となります。よって

$x = \dfrac{1}{3}te^{2t}$

が特殊解であることがわかります。以上から, 非斉次方程式の一般解は

$x = C_1e^{-t} + C_2e^{2t} + \dfrac{1}{3}te^{2t}~~$ ($C_1$, $C_2$ は任意定数)

となります。

※注意

ここで特殊解を $x = Ae^{2t}$ と仮定してしまうと非斉次方程式は

$4Ae^{2t} -2Ae^{2t} -2Ae^{2t} = 0$

となってしまうので, 係数 $A$ を求めることができません。

※別解(特殊解の求め方)

ロンスキアンを計算すると

$W(e^{-t},e^{2t}) = \left| ete2tet2e2t \right| = 3e^t$

不定積分の部分を計算すると

$x2R(t)W(x1,x2) dt=e4t3et dt=13e3t dt=19e3t$

同様に

$x1R(t)W(x1,x2) dt=et3et dt=13t$

よって特殊解は

$x=x1x2R(t)W(x1,x2) dt+x2x1R(t)W(x1,x2) dt=19e2t+13te2t$

となります。$-\dfrac{1}{9}e^{2t}$ は一般解に含まれていることに注意しましょう。

$Q5$.
次の微分方程式の一般解を求めなさい。

$\dfrac{d^2x}{dt^2} + 5\dfrac{dx}{dt} + 6x = 5\sin 2t$
解答・解説を見る
$x = C_1e^{-3t} + C_2e^{-2t} + \dfrac{5}{52}\sin 2t - \dfrac{25}{52}\cos 2t~~$ ($C_1$, $C_2$ は任意定数)

$\dfrac{d^2x}{dt^2} + 5\dfrac{dx}{dt} + 6x=0$ とすると, 特性方程式は

$\lambda^2 + 5\lambda + 6 = (\lambda + 3)(\lambda + 2) =0$

$\lambda = -3,~-2$ より, 斉次方程式の一般解は $x = C_1e^{-3t}+ C_2e^{-2t}$ となります。

次に, 非斉次方程式の右辺が三角関数なので, 特殊解を

$x = A\sin 2t + B\cos 2t$

と仮定すると

$\dfrac{dx}{dt} = 2A\cos 2t - 2B\sin 2t$

$\dfrac{d^2x}{dt^2} = -4A\sin 2t - 4B\cos 2t$

非斉次方程式に代入して整理すると

$(2A -10B)\sin 2t + (10A +2B)\cos 2t = 5\sin 2t$

係数を比べると

${2A10B=510A+2B=0$

これを解くと $A =\dfrac{5}{52}$, $B = -\dfrac{25}{52}$ となります。よって

$x = \dfrac{5}{52}\sin 2t - \dfrac{25}{52}\cos 2t$

が特殊解であることがわかります。以上から, 非斉次方程式の一般解は

$x = C_1e^{-3t} + C_2e^{-2t} + \dfrac{5}{52}\sin 2t - \dfrac{25}{52}\cos 2t~~$ ($C_1$, $C_2$ は任意定数)

となります。