指数関数と微分演算子

指数関数に微分演算子を代入すると、適当な条件の下で

\[f(x+1) = \exp \left( \frac{d}{dx} \right) f(x) \]

となることを確認し、差分方程式の解法への応用を考える。(厳密な議論はしない。)

指数関数に微分演算子を代入する

指数関数のテイラー展開に微分演算子を代入して

\[\exp \left( \frac{d}{dx} \right ) =1 + \frac{d}{dx} + \frac{1}{2!} \left( \frac{d}{dx} \right) ^2 + \frac{1}{3!} \left( \frac{d}{dx} \right) ^3 +\cdots \]

について考える。これを関数に作用させると

\[\exp \left( \frac{d}{dx} \right ) f(x) = f(x) + f'(x) + \frac{f''(x)}{2!} + \cdots \] となるが、右辺は $f(x+1)$ と適当な条件下で等しくなる。つまり、

\[f(x+1) = f(x) + f'(x) + \frac{f''(x)}{2!} + \cdots \]

が適当な条件下で成立する。なぜならば

$f(w)$の$a$におけるテイラー展開を考えると \[f(w) = f(a) + f'(a)(w-a) + \frac{f''(a)}{2!}(w-a)^2 + \cdots\]

となるが、ここで、$a=x,\, w=x+1$を代入することで \[f(x+1) = f(x) + f'(x) + \frac{f''(x)}{2!} + \cdots \] が得られるからである。

例1 $ f(x) = x^2 $ の場合.

\[\exp \left( \frac{d}{dx} \right ) f(x) = \exp \left( \frac{d}{dx} \right ) x^2 =\left(1 + \frac{d}{dx} + \frac{1}{2!} \left( \frac{d}{dx} \right) ^2 +\cdots \right) x^2 = x^2 +2x +1 = (x+1)^2 \]

例2 $f(x) = e^x$ の場合.

\[\exp \left( \frac{d}{dx} \right ) f(x) = \exp \left( \frac{d}{dx} \right ) e^x =\left(1 + \frac{d}{dx} + \frac{1}{2!} \left( \frac{d}{dx} \right) ^2 +\cdots \right) e^x = (1+1+\frac{1}{2!}+\cdots ) e^x = e^{x+1} \]

差分方程式への応用

\[f(x+1)-f(x) = g(x)\] の形の差分方程式の解法について考える。

これを先ほど得た関係式を用いて書きかえると

\[ \left ( \exp \left( \frac{d}{dx} \right ) -1 \right) f(x) = g(x)\]

となる。これを$f(x)$について解くと

\[ f(x) = \frac{1}{ \exp \left( \frac{d}{dx} \right ) -1} g(x) \]

が得られ、形式的に差分方程式の解が得られる。右辺は

\[ \frac{1}{ \exp \left( \frac{d}{dx} \right ) -1} g(x)=\sum_{n=0}^\infty \frac{B_n}{n!} \left( \frac{d}{dx} \right)^{n-1} g(x)\]

と、$(e^x-1)^{-1}$のローラン展開を用いて計算する。ただし、

\[\left( \frac{d}{dx} \right) ^{-1} \] は積分と解釈する。

例. べき乗和の公式 ($g(x)=(x+1)^2$のケース.)

$f(x) = 1^2 + 2^2 + 3^2 + \cdots + x^2 $ とおくと

$f(x+1) - f(x) = (x+1)^2$ となる。したがって

\begin{align} f(x) &= \frac{1}{ \exp \left( \frac{d}{dx} \right ) -1} (x+1)^2 \\ &= \left( \left( \frac{d}{dx} \right)^{-1} + B_1 + \frac{B_2}{2}\frac{d}{dx} + \frac{B_3}{3!}\frac{d^2}{dx^2}+\cdots \right) (x+1)^2 \\ &= \left( \left( \frac{d}{dx} \right)^{-1} -\frac{1}{2} + \frac{1}{12}\frac{d}{dx} + \cdots \right) (x+1)^2 \\ &= \int (x+1)^2 dx -\frac{1}{2}(x+1)^2 + \frac{1}{12}\frac{d}{dx}(x+1)^2 \\ &= \frac{(x+1)^3}{3} -\frac{(x+1)^2}{2} + \frac{(x+1)}{6}+C \\ &= \frac{x\left( x+1\right) \left( 2 x+1\right) }{6} +C\end{align}

が得られるが $f(0)=0$ とすれば $C=0$ となるので、結局 \[f(x) = \frac{x\left( x+1\right) \left( 2 x+1\right) }{6} \] が導出できる。これはべき乗和の公式に他ならない。