倭算数理研究所

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

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

特殊関数の公式を証明していくシリーズ(目次)。 今回は、以下の合流型超幾何微分方程式 (confluent hypergeometric differential equation)

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

(ただし { a,\,c } は定数)に関して、級数解を求めてみましょう。 まぁ、それは合流型超幾何関数になるんですが。

参考

導出

{ w(z) } を以下の冪級数の形で求めましょう:

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

この級数が合流型超幾何微分方程式 (*) を満たすという条件から展開係数 { w_n } の満たす漸化式を導き、{ w_n } の一般項を求めます。

{ w_n } の漸化式を求める
では、まず展開係数 { w_n } が満たす漸化式を求めましょう。 まぁ、級数の形の { w(z) }微分方程式 (*) にぶち込めばいいんですが、一気にやるとゴチャゴチャするのでステップ・バイ・ステップでやっていきましょう。 次の4つを { z } の冪級数で表します:

  { \displaystyle \begin{align*}
  w,\quad \frac{dw}{dz},\quad z\frac{dw}{dz}, \quad z\frac{d^2w}{dz^2}
\end{align*}}

後の計算のため、(途中はともかく)最終的には級数を以下のような形に整理することにします:

  • 和は { n=0 } から { n=\infty } まで
  • 全て { z }{ n } 乗({ n-1 } 乗とかではなく)で書く
  • 分母には { n! }{ (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=1}^\infty \frac{w_n}{(n-1)!}z^{n-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{d^2 w}{dz^2} &= \sum_{n=0}^\infty \frac{nw_{n+1}}{n!}z^n &\cdots(5) \qquad (\because (3) )
\end{align*}}

{ n } についての和で { n=0 } の項が実際には消えるところもありますが、和の範囲を揃えておく方が後が楽。 で、(1), (3), (4), (5) 式を合流型超幾何微分方程式 (*) の左辺に代入して変形していくと

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

これが { z } に関して恒等的に0となるとして

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

もう少し書き換えると(簡単のため { c } は非正の整数ではないとします。 つまり、以下の式で全ての { n } に対して分母が0にならないとします。)

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

おっ、綺麗な漸化式が出てきたっ。 まぁ、そうなるようにうまく微分方程式の方を作ってるんでしょうけどね。

{ w_n } の一般項を求める
上記で得られた { w_n } の一般項はそんなに難しくなく求まります。 初項 { w_0 } とりあえずそのままにして、{ w_n \, (n \ge 1) }{ w_0 } で表しましょう。

  { \displaystyle \begin{align*}
  w_n
      &= \frac{a+n-1}{c+n-1}w_{n-1} \\
      &= \frac{(a+n-1)(a+n-2)}{(c+n-1)(c+n-2)}w_{n-2} \\
          &\qquad \vdots \\ &= \frac{(a+n-1)(a+n-2)\cdots a}{(c+n-1)(c+n-2)\cdots c}w_0 \\
      &= \left(\prod_{k=0}^{n-1}\frac{a+k}{c+k}\right)w_0
\end{align*}}

ここでポッホハマー記号 (Pochhammer symbol) { (x)_n }

  { \displaystyle \begin{align*}
  (x)_n = \begin{cases}
    1 & (n = 0) \\[2mm]
    \displaystyle \prod_{m=0}^{n-1}(x+m) \quad \Big(= x(x+1)\cdots (x+n-1) \Big) & (n \geqq 1)
  \end{cases}
\end{align*}}

によって導入すると

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

となります。 これは { n=0 } のときも成り立ちます。 ちなみに、ポッホハマーの記号は以前の記事「下降階乗冪ってのがあったので上昇階乗冪ってのも試してミルカ」でやった、上昇階乗冪というのと同じです。

ちなみに { a } が非正の整数の場合、つまり { a = -m \; (m=0,\,1,\,2,\,\cdots) } と表されるとき、{ w_n } の漸化式より

  { \displaystyle \begin{align*}
  w_{m+1} &= \frac{-m+m}{c+m}w_m = 0
\end{align*}}

となり、また { n > m+1 } を満たす整数に対しても { w_n = 0 } となります。

{ w(z) }級数
上記で求めた { w_n }{ w(z) } の冪級数に代入すると

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

定数 { w_0 } を除いた級数

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

合流型超幾何関数 (confluent hypergeometric function)と呼びます。 前節最後に言及した事項より、{ a = -m \; (m = 0,\,1,\,2,\,\cdots) } のとき、{ F(a,\,c;\,z) }{ m } 次の多項式になります。

合流型超幾何微分方程式のもう一つの基本解

さて、合流型超幾何微分方程式は2階の線形常微分方程式ですが、こういう微分方程式には2つの基本解*1があるそうですね。 1つは上記の合流型超幾何関数 { F(a,\,c;\,z) } ですが、もう1つも合流型超幾何関数を使って

  { \displaystyle \begin{align*}
  z^{1-c}F(a-c+1,\,2-c;\,z)
\end{align*}}

と書けるそうです(ただし { c } が正の整数でない場合)。 これをちょっと確かめてみましょう。

解の形
上記の関数をそのまま微分方程式に代入して成り立っているのを確かめてもいいんですが、{ w(z) }

  { \displaystyle \begin{align*}
  w(z) &= z^\omega u(z) & \cdots (6)
\end{align*}}

という形に仮定して({ \omega } は後で決める定数)、合流型超幾何微分方程式{ u(z) }微分方程式に書き換えてその解の形を考えるという方法をとります。

{ u(z) } の満たすべき微分方程式
まず後の計算のために (5) 式の { w(z) }{ z }微分しておきましょう。

  { \displaystyle \begin{align*}
  \frac{dw}{dz} &= \omega z^{\omega-1} u + z^\omega \frac{du}{dz} \\
  \frac{d^2w}{dz^2}
    &= \omega(\omega-1)z^{\omega-2}u + 2\omega z^{\omega-1}\frac{du}{dz} + z^\omega \frac{d^2u}{dz^2}
\end{align*}}

さて、これらを合流型超幾何微分方程式の左辺に (5) 式を代入して変形していくと

  { \displaystyle \begin{align*}
  & z\frac{d^2w}{dz^2} - (c-z)\frac{dw}{dz} -a w \\
  &\quad = \omega(\omega-1)z^{\omega-1}u + 2\omega z^\omega\frac{du}{dz} + z^{\omega+1} \frac{d^2u}{dz^2} \\
  &\quad\qquad + c\omega z^{\omega-1} u + cz^\omega \frac{du}{dz} \\
  &\quad\qquad -\omega z^\omega u - z^{\omega+1} \frac{du}{dz} \\
  & \quad\qquad - a z^\omega u \\
  &\quad = z^\omega\left[z\frac{d^2u}{dz^2} + (2\omega+c-z)\frac{du}{dz}
      + \left\{\frac{\omega(\omega+c-1)}{z}-(a+\omega)\right\}u\right]
\end{align*}}

ここで { \omega = 1 - c } とおくと、微分方程式

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

となり、大括弧内 { [ \cdots ] } はパラメータを

  { \displaystyle \begin{align*}
  \begin{cases}
      a \longrightarrow a-c+1 \\
      c \longrightarrow 2-c
  \end{cases}
\end{align*}}

と置き換えた合流型超幾何微分方程式となります。 よって { u(z) } は合流型超幾何関数を使って書けて

  { \displaystyle \begin{align*}
  u(z) = F(a-c+1,\,2-c;\,z)
\end{align*}}

となります。 { \omega = 1 - c } も合わせて、

  { \displaystyle \begin{align*}
  z^{1-c}F(a-c+1,\,2-c;\,z)
\end{align*}}

が合流型超幾何微分方程式を満たすことが分かります。 これが { F(a,\,c;\,z) } と線型独立だとかいうのは他を当たってもらうことにしましょう。

同様のことを(合流型でない)超幾何微分方程式についてもやる予定ですが・・・

まとめ

合流型超幾何関数

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

の解は、 { c } が非正の整数でない場合、合流型超幾何関数

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

  { \displaystyle\begin{align*}
  z^{1-c}F(a-c+1,\,2-c;\,z)
\end{align*}}

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

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

です。

【修正】

  • 展開係数  { a_n } はパラメータの  { a } と混同する可能性がないでもないので、 { w_n } に変更しました。

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

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

*1:一般解をそれらの線型結合で書ける、って定義だっけ?