今日の講義で微分方程式の特性根と一般解の話をしたら,重根をもつ場合について,高校で習う隣接3項漸化式で重根を持つ場合の解き方と同じ考え方で微分方程式も説明できることが分かったと講義の後で教えてくれました.以下にまとめてみます.
隣接3項漸化式で特性方程式が重根をもつ場合
[latex]y(k+2)=2ay(k+1)-a^2y(k)[/latex]
を変形して,差が等比数列の形になるようにする.
[latex]y(k+2)-ay(k+1)=a(y(k+1)-ay(k))[/latex]
右辺は
[latex]a^{k+1}(y(1)-ay(0))[/latex]
となる.
両辺を[latex]a^{k+1}[/latex]で割る.
[latex]
\displaystyle\frac{y(k+2)}{a^{k+1}}-\frac{y(k+1)}{a^k}=y(1)-ay(0)
[/latex]
となり,[latex]\displaystyle\frac{y(k+1)}{a^k}[/latex]は,初項[latex]y(0)/a^{-1}[/latex],公差[latex]y(1)-ay(0)[/latex]の等差数列になる.
したがって
[latex]
\displaystyle\frac{y(k)}{a^{k-1}}=\frac{y(0)}{a^{-1}} + (k-1)(y(1)-ay(0))
[/latex]
となる.
両辺を[latex]a^{k-1}[/latex]倍して
[latex]
y(k)=a^ky(0)+a^{k-1}(k-1)(y(1)-ay(0))
[/latex]
になる.
隣接3項漸化式では,[latex]y(k)[/latex]を求めることになるが,線形差分方程式とみた場合には,基底に[latex]a^k[/latex]と[latex]ka^k[/latex]を用いればよいことが分かる.
微分方程式
微分方程式の特性根が重根の場合,
[latex]
\displaystyle\frac{d^2y(t)}{dt^2} -2a\frac{dy(t)}{dt} + a^2y(t)=0
[/latex]
を変形して
[latex]
\displaystyle\frac{d^2y(t)}{dt^2} – a\frac{dy(t)}{dt} = a\left(\frac{dy(t)}{dt} – ay(t)\right)
[/latex]
となる.
これは,
[latex]
\displaystyle\frac{d}{dt}\left(\frac{dy(t)}{dt}-ay(t)\right)= a\left(\frac{dy(t)}{dt} – ay(t)\right)
[/latex]
なので,
[latex]
\displaystyle\frac{dy(t)}{dt}-ay(t) = e^{at}\left(\frac{dy(0)}{dt} – ay(0)\right)
[/latex]
になる.
[latex]e^{at}[/latex]で両辺を割ると,
[latex]
e^{-at}\displaystyle\frac{dy(t)}{dt}-ae^{-at}y(t) = \frac{dy(0)}{dt} -ay(0)
[/latex]
なので,左辺を次のように書き換えて
[latex]
\displaystyle\frac{d}{dt}\left(e^{at}y(t)\right) = \frac{dy(0)}{dt} -ay(0)
[/latex]
になる.
積分して,
[latex]
e^{-at}y(t) = \left(\displaystyle\frac{dy(0)}{dt} -ay(0)\right)t + c
[/latex]
になり([latex]c[/latex]は積分定数),
[latex]
y(t)=\left(\displaystyle\frac{dy(0)}{dt} -ay(0)\right)te^{at} + ce^{at}
[/latex]
になる,
つまり,基底は[latex]e^{at}[/latex]と[latex]te^{at}[/latex]になる.
実際問題として,重根のときの基底を導出する時間はないでしょうから,重根のときの基底を覚えておくとよいでしょう.また,講義で説明した定数変化法を理解しておくとよいと思います.