倭算数理研究所

科学・数学・学習関連の記事を、「倭マン日記」とは別に書いていくのだ!

実録! 超幾何微分方程式の級数解を求めてみる。 もちろんそれは超幾何関数になる。

特殊関数の公式を証明していくシリーズ(目次)。 前回の記事『合流型超幾何微分方程式の級数解を求めてみる。 もちろんそれは合流型超幾何関数になる。』で、合流型超幾何微分方程式級数解を構築していくと合流型超幾何関数が得られることを示しました。 今回はそれを超幾何微分方程式について行います。 ちなみに、「合流型 (confluent)」という小難しい修飾子が付いている方が実際には簡単な微分方程式、関数になります。 名前の由来は、超幾何微分方程式特異点(3つあるうちの2つ)が合流することから来ています。 特異点の数が減るので「合流型」の方が簡単になる、という寸法です。

さて、今回の本題の超幾何微分方程式 (hypergeometric function) とは、以下のような微分方程式です:

  { \displaystyle\begin{align*}
    z(1-z)\frac{d^2w}{dz^2} + \left[c - (a+b+1)z\right]\frac{dw}{dz} - ab w &= 0 & \cdots (*)
\end{align*}}

 { a,\,b,\,c } は定数です。

参考

級数解の導出

上記の微分方程式 (*) の級数解を前回と同様の手順で求めていきましょう。 級数解は以下の形で定義します:

  { \displaystyle\begin{align*}
    w(z) &= \sum_{n=0}^\infty \frac{w_n}{n!}z^n & \cdots (1)
\end{align*}}

手順は

  1. 微分方程式 (*) から、{ w_n } が満たす漸化式を導く
  2. { w_n } の一般項を求める

です。

{ w_n } の漸化式
まず、前回と同様に (1) 式の微分等を計算しましょう:

  { \displaystyle\begin{align*}
    \frac{dw}{dz}
        &= \sum_{n=0}^\infty \frac{nw_n}{n!}z^{n-1} & \cdots (2) \qquad (\because (1)) \\
        &= \sum_{n=0}^\infty \frac{w_{n+1}}{n!}z^n & \cdots (3) \\
    z\frac{dw}{dz} &= \sum_{n=0}^\infty \frac{nw_n}{n!}z^n &\cdots (4) \qquad (\because (2)) \\
    z\frac{dw^2}{dz^2} &= \sum_{n=0}^\infty \frac{nw_{n+1}}{n!}z^n & \cdots (5) \qquad (\because (4)) \\
    z^2\frac{dw^2}{dz^2} &= \sum_{n=0}^\infty \frac{n(n-1)w_n}{n!}z^n & \cdots (6) \qquad (\because (2))
\end{align*}}

これを微分方程式 (*) の左辺に代入して変形していくと

  { \displaystyle\begin{align*}
  &z(1-z)\frac{d^2w}{dz^2} + \left[c - (a+b+1)z\right]\frac{dw}{dz} - ab w \\
  &\quad = \sum_{n=0}^\infty \Big\{nw_{n+1} - n(n-1)w_n + c w_{n+1} - (a+b+1) nw_n - ab w_n\Big\}\frac{z^n}{n!}
\end{align*}}

この級数{ z } に関して恒等的に0になるとして

  { \displaystyle\begin{align*}
    & nw_{n+1} - n(n-1)w_n + c w_{n+1} - (a+b+1) nw_n - ab w_n = 0 \\
    \therefore \;& (c+n)w_{n+1} = \left\{n(n-1) + (a+b+1)n + ab\right\}w_n
\end{align*}}

右辺の中括弧内 {...} をもう少し整理すると

  { \displaystyle\begin{align*}
    n(n-1) + (a+b+1)n + ab
        &= (b+n)a + n(n-1) + (b+1)n \\
        &= (b+n)a + n(b+n) \\
        &= (a+n)(b+n)
\end{align*}}

というふうに綺麗にまとまります。 もとの漸化式に戻ると

  { \displaystyle\begin{align*}
    (c+n)w_{n+1} = (a+n)(b+n)w_n
\end{align*}}

よって { w_n } についての漸化式は(簡単のため { c } は非正の整数ではないとして)

  { \displaystyle\begin{align*}
    w_{n+1} &= \frac{(a+n)(b+n)}{c+n} w_n & (\textrm{for}\;n \ge 0)
\end{align*}}

となります。

{ w_n } の一般項
上記の漸化式から、 { w_n } の一般項は

  { \displaystyle\begin{align*}
  w_n &= \left(\prod_{m=0}^{n-1}\frac{(a+m)(b+m)}{c+m}\right)w_0
\end{align*}}

となることが分かります。 前回導入したポッホハマー記号 { (x)_n }

  { \displaystyle\begin{align*}
  (x_n) = 
    \begin{cases}
      1 & (n = 0) \\[2mm]
      \displaystyle \prod_{k=0}^{n-1}(x + k) & (n \geqq 1)
  \end{cases}
\end{align*}}

を使うと、{ w_n } の一般項は

  { \displaystyle\begin{align*}
    w_n = \frac{(a)_n(b)_n}{(c)_n}w_0
\end{align*}}

となります。

超幾何関数
前節で得られた { w_n } を (1) 式に代入すると

  { \displaystyle\begin{align*}
    w(z) = w_0\sum_{n=0}^\infty \frac{(a)_n(b)_n}{(c)_n}\frac{z^n}{n!}
\end{align*}}

定数 { w_0 } を除いた級数

  { \displaystyle\begin{align*}
    F(a,\,b,\,c;\,z) = \sum_{n=0}^\infty \frac{(a)_n(b)_n}{(c)_n}\frac{z^n}{n!}
\end{align*}}

を(Gauss の)超幾何関数 (hypergeometric function) と呼びます。 合流型超幾何微分方程式の場合と同様に、超幾何微分方程式の他の基本解も超幾何関数を使って書けるようですが、拙者の能力の許容量を超えてるのでスルー。

合流型超幾何微分方程式と超幾何微分方程式

超幾何微分方程式で、変数を { z \rightarrow z/b } に変換して { b \rightarrow \infty } の極限をとると合流型超幾何微分方程式になります。 これは超幾何微分方程式特異点 { z = 0,\,1,\,\infty } のうち、{ z = 1 } のものを { z = 0 } と「合流」させることに対応します({ z = 0,\,\infty }特異点はそのままに)。 これと同様の操作を超幾何関数に対して行うと合流型超幾何関数が得られます。 実際にやってみましょう。

超幾何微分方程式
まずは微分方程式の方。 変数変換 { z \rightarrow z/b } の下で微分の変換は

  { \displaystyle\begin{align*}
    \frac{d}{dz} \longrightarrow \frac{d}{d(z/b)} = b\frac{d}{dz}
\end{align*}}

となるので微分方程式 (*) の左辺は

  { \displaystyle\begin{align*}
    & z(1-z)\frac{d^2w}{dz^2} + \left[c - (a+b+1)z\right]\frac{dw}{dz} - abw \\
    &\quad\longrightarrow \frac{z}{b}\left(1-\frac{z}{b}\right)b^2\frac{d^2w}{dz^2}
        + \left[c - (a+b+1)\frac{z}{b}\right]b\frac{dw}{dz} - abw \\
    &\qquad = b\left\{z\left(1-\frac{z}{b}\right)\frac{d^2w}{dz^2}
        + \left[c - \left(1 + \frac{a+1}{b}\right)z\right]\frac{dw}{dz} - aw\right\}
\end{align*}}

と変換されます。 よって微分方程式 (*) は(両辺を { b } で割って)

  { \displaystyle\begin{align*}
    z\left(1-\frac{z}{b}\right)\frac{d^2w}{dz^2} + \left[c - \left(1 + \frac{a+1}{b}\right)z\right]\frac{dw}{dz} - aw = 0
\end{align*}}

となります。 ここで { b \rightarrow \infty } の極限をとると

  { \displaystyle\begin{align*}
    z\frac{d^2w}{dz^2} + (c - z)\frac{dw}{dz} - aw = 0
\end{align*}}

となり、合流型超幾何微分方程式が得られます。

超幾何関数
微分方程式の場合と同じ変換&極限を、超幾何関数 { F(a,\,b,\,c;\,z) } にも施してみましょう。

  { \displaystyle\begin{align*}
    F(a,\,b,\,c;\,z)
        &= \sum_{n=0}^\infty \frac{(a)_n(b)_n}{(c)_n}\frac{z^n}{n!} \\
        &\rightarrow \sum_{n=0}^\infty \frac{(a)_n(b)_n}{(c)_n}\frac{1}{n!}\left(\frac{z}{b}\right)^n \\
        &= \sum_{n=0}^\infty \frac{(a)_n}{(c)_n}\frac{(b)_n}{b^n}\frac{z^n}{n!}
\end{align*}}

ここでポッホハマー記号 { (b)_n } の定義より

  { \displaystyle\begin{align*}
    \frac{(b)_n}{b^n}
        &= \frac{1}{b^n}\prod_{k=0}^{n-1}(b+k)
        = \prod_{k=0}^{n-1}\left(1 + \frac{k}{b}\right) \\
        &\longrightarrow 1 & (b \rightarrow \infty)
\end{align*}}

なので

  { \displaystyle\begin{align*}
    F(a,\,b,\,c;\,z)
        &\longrightarrow \sum_{n=0}^\infty \frac{(a)_n}{(c)_n}\frac{z^n}{n!} & (b\rightarrow \infty)\\
        &\quad= F(a,\,c;\,z)
\end{align*}}

となって、{ b \rightarrow \infty } の極限で超幾何関数は合流型超幾何関数になることが分かりました。

まとめ

超幾何微分方程式

  { \displaystyle\begin{align*}
    z(1-z)\frac{d^2w}{dz^2} + \left[c - (a+b+1)z\right]\frac{dw}{dz} - ab w &= 0
\end{align*}}

級数解は

  { \displaystyle\begin{align*}
  F(a,\,b,\,c;\,z)
    &= \sum_{n=0}^\infty \left(\prod_{m=0}^{n-1}\frac{(a+m)(b+m)}{c+m}\right)\frac{z^n}{n!} \\[2mm]
    &= \sum_{n=0}^\infty \frac{(a)_n(b)_n}{(c)_n}\frac{z^n}{n!}
\end{align*}}

で与えられます。 ただし { (x)_n } はポッホハマー記号

  { \displaystyle\begin{align*}
  (x_n) = 
    \begin{cases}
      1 & (n = 0) \\[2mm]
      \displaystyle \prod_{k=0}^{n-1}(x + k) & (n \geqq 1)
  \end{cases}
\end{align*}}

です。

【修正】

  • ポッホハマーの公式の定義で符号が間違っていたので修正しました。
  • 展開係数として使っていた  { a_n } はパラメータの  { a } と混同のおそれがないでもないので、展開係数を  { w_n } に変更しました。

基礎物理数学第4版Vol.3 特殊関数 (KS理工学専門書)

基礎物理数学第4版Vol.3 特殊関数 (KS理工学専門書)