応用数学 第8回 (2) 級数解法

級数解法の例 その1

 微分方程式の解が整級数の形であると仮定して、 係数比較によって解を求める解法を級数解法と呼びます。
Ex.2  $y'=y^2$,  $y(0)=1$
 $y=y(x)$ が整級数 $(\ast)$ の形であると仮定します。 まず
$c_0=y(0)=1$
から
$y = 1 + c_1\, x + c_2\, x^2 + c_3\, x^3 + \cdots $
であり、
$y^2 = 1 + 2\,c_1\, x + (c_1^2+2\,c_2)\, x^2 + 2\,(c_1c_2+c_3)\,x^3 + \cdots $
となります。これと
$y' = c_1 + 2\,c_2\, x + 3\, c_3\, x^2 + 4\, c_4 \, x^3 + \cdots $
を係数比較して \begin{align} c_1 & = 1 \\ 2\, c_2 &= 2 \, c_1 = 2    \therefore   c_2 = 1 \\ 3\, c_3 &= c_1^2 + 2 \, c_2 = 3    \therefore   c_3 = 1 \\ 4\, c_4 &= 2\,(c_1c_2 + c_3) = 4    \therefore   c_4 = 1 \\ \end{align} どうやら
$y = 1 + x + x^2 + x^3 + x^4 + \cdots = \dps{\frac{1}{1-x}}$
のようです。

検算 $y(0) = 1$ は良くて、
$y' = \left((1-x)^{-1}\right)' = (-1)(1-x)^{-2}(1-x)' = (1-x)^{-2} = y^2$.
正解です。

級数解法の例 その2

Ex.3  $y'' - xy' - y = 0$
 今度は $\sum$ を使った式で計算しましょう。
$\dps{y = \sum_{n=0}^{\infty} c_n x^n}$
とすると
$\dps{y' = \sum_{n=0}^{\infty} n \, c_n x^{n-1}}$
より
$\dps{xy' = \sum_{n=0}^{\infty} n \, c_n x^n}$
また
$\dps{y'' = \sum_{n=2}^{\infty} n(n-1) \, c_n\, x^{n-2} = \sum_{n=0}^{\infty} (n+2)(n+1) \, c_{n+2}\, x^n}$
よって
左辺の $x^n$ の係数 $= (n+2)(n+1)\,c_{n+2} - n\,c_n - c_n = (n+1)\big\{(n+2)\,c_{n+2}-c_n\big\}$
これが $0$ なので
$\dps{c_{n+2} = \frac{1}{n+2}\,c_n},$   $\forall\,n$
すると $n$ が偶数のとき $c_n$ は $c_0$ で書けて
$\dps{c_2 = \frac{1}{2}c_0}$,   $\dps{c_4 = \frac{1}{4}c_2 = \frac{1}{4\times 2}c_0}$,   $\dps{c_6 = \frac{1}{6}c_4 = \frac{1}{6\times 4\times 2}c_0},$   $\cdots$
また $n$ が奇数のとき $c_n$ は $c_1$ で書けて
$\dps{c_3 = \frac{1}{3}c_1}$,   $\dps{c_5 = \frac{1}{5}c_3 = \frac{1}{5\times 3}c_1}$,   $\dps{c_7 = \frac{1}{7}c_5 = \frac{1}{7\times 5\times 3}c_1},$   $\cdots$
よって \begin{align} y & = c_0 + c_1 x + \frac{1}{4\times 2}c_0\,x^2 + \frac{1}{5\times 3}c_1\,x^3 + \cdots \\ & = \dps{ c_0 \left( 1 + \frac{1}{2}x^2 + \frac{1}{4\times 2}x^4 + \frac{1}{6\times 4\times 2}x^6 + \cdots \right) } \\ & + \dps{ c_1 \left( x + \frac{1}{3}x^3 + \frac{1}{5\times 3}x^5 + \frac{1}{7\times 5\times 3}x^7 + \cdots \right)} \\ \end{align} すなわち \begin{align} f & = \dps{1 + \frac{1}{2}x^2 + \frac{1}{4\times 2}x^4 + \frac{1}{6\times 4\times 2}x^6 + \cdots} \\ g & = \dps{x + \frac{1}{3}x^3 + \frac{1}{5\times 3}x^5 + \frac{1}{7\times 5\times 3}x^7 + \cdots} \\ \end{align} と置けば、
$y=c_0\,f + c_1 \, g$   ( $c_0$, $c_1$ は任意定数 )
と書け、2階線形微分方程式の一般解(2次元)の基底が $f$ と $g$ ということになります。