分数階微積分学

分数階微分積分学（ぶんすうかいびぶんせきぶんがく、英: fractional calculus）は解析学（特に微分積分学）の一分野で、微分作用素 D および積分作用素 J が実数冪あるいは複素数冪をとる可能性について研究する学問である。 この文脈における「冪」の語は作用素の合成を繰り返し行うという意味で用いており、それに従えばたとえば f2(x) = f(f(x)) ということになる。さてたとえば、微分作用素 D の平方根（あるいは微分を半分だけ作用させる）という意味での式 D = D 1 / 2 {\displaystyle {\sqrt {D}}=D^{1/2}} に何か意味のある解釈をつけられるかということを考える。この式は、つまりある作用素を「二度」作用させて、微分作用素 D と同じ効果を得られるということを意味しているのであり、あるいはもっと一般に、実数 s に対して微分作用素の冪 D s {\displaystyle D^{s}} にあたるものを決定できるかという問をも考えることができるだろう。このとき、s が整数 n を値にとるならば、n > 0 のときこの冪は通常の意味での n-階微分作用素となり、n < 0 のときは積分作用素 J の (−n)-乗となるように定義されるものでなければならない。 このようなことを考える理由はいくつかある。ひとつはそれによって「離散」的な変数 n で添字付けられる微分作用素の族 Dn 全体が作る半群を実数 s を径数とする「連続」的な半群のなかにあるとして考えられるようになることである。連続的半群というものは数学のさまざまなところに現われ、豊かな理論を備えている。分数階微分積分学では、冪として必ずしも有理数冪に限らず実数冪や複素数冪を一般に扱うため「分数階」という名称で呼ぶのは少々紛らわしいが、慣習的に「分数階微分積分学」の名称が使われている。 このような理論の存在については、1832年からのリウヴィルの論文にその素地を見ることができる。函数の階数 a の分数階微分は今日ではしばしばフーリエ変換あるいはメリン変換といった積分変換の意味で定義される。重要なことは、点 x における分数階微分というものが「局所的」な概念であるのは、a が整数値をとる場合に限られるという点である。つまり、非整数階の場合には、函数 f の点 x における分数階微分が x の極近くでの f のグラフのみに依存して決まるということができない（整数階微分であればこれが言える）。然るに、分数階微分作用素の理論においてはある種の境界条件や函数についてのさらなる情報が関わってくることが想定される。喩えるならば、分数階微分はある種の周辺視野を要求するのである。 まず相応に自然な疑問は、半微分 (英: half-derivative) と呼ばれるべき、作用素 H で H 2 f ( x ) = D f ( x ) = d d x f ( x ) = f ′ ( x ) {\displaystyle H^{2}f(x)=Df(x)={\frac {d}{dx}}f(x)=f'(x)} を満たすものは存在するかということであろう。そのような作用素は存在する。実際には任意の実数値 a > 0 に対して ( P a f ) ( x ) = D f ( x ) = d d x f ( x ) = f ′ ( x ) {\displaystyle (P^{a}f)(x)=Df(x)={\frac {d}{dx}}f(x)=f'(x)} を満たす作用素 P が存在することがいえる。言い方を変えれば、n-階微分 dny⁄dxn の n を任意の実数値に拡張することができるのである。 もう少し詳しく述べるに、階乗の非整数値への拡張としてのガンマ関数 Γ から始める。ガンマ関数は n ! = Γ ( n + 1 ) {\displaystyle n!=\Gamma (n+1)} を満たしていることを利用する。さて函数 f(x) は x > 0 で矛盾なく定義され、0 から x までの定積分 ( J f ) ( x ) = d − 1 d x − 1 f ( x ) = ∫ 0 x f ( t ) d t {\displaystyle (Jf)(x)={d^{-1} \over dx^{-1}}f(x)=\int _{0}^{x}f(t)\;dt} . ができるものと仮定する。これを繰り返して ( J 2 f ) ( x ) = d − 2 d x − 2 f ( x ) = ∫ 0 x ( J f ) ( t ) d t = ∫ 0 x ( ∫ 0 t f ( s ) d s ) d t {\displaystyle (J^{2}f)(x)={d^{-2} \over dx^{-2}}f(x)=\int _{0}^{x}(Jf)(t)\;dt=\int _{0}^{x}\left(\int _{0}^{t}f(s)\;ds\right)dt} や任意の自然数冪 Jnf に拡張することができるが、反復積分に関するコーシーの公式によれば ( J n f ) ( x ) = d − n d x − n f ( x ) = 1 ( n − 1 ) ! ∫ 0 x ( x − t ) n − 1 f ( t ) d t {\displaystyle (J^{n}f)(x)={d^{-n} \over dx^{-n}}f(x)={1 \over (n-1)!}\int _{0}^{x}(x-t)^{n-1}f(t)\;dt} である。また、階乗函数を用いる代わりに（Γ(n + 1) = n! あるいは同じことだが Γ(n) = (n − 1)! の関係にある）ガンマ関数に置き換えれば、 　　 ( J n f ) ( x ) = d − n d x − n f ( x ) = 1 Γ ( n ) ∫ 0 x ( x − t ) n − 1 f ( t ) d t {\displaystyle (J^{n}f)(x)={d^{-n} \over dx^{-n}}f(x)={1 \over \Gamma (n)}\int _{0}^{x}(x-t)^{n-1}f(t)\;dt} とも表せる。これを用いれば直接に n が実数だけでなく複素数である場合にまで一般化することができる。すなわち、ガンマ関数を複素数の範囲まで広げることにより、積分作用素を「分数（負の整数を除く複素数）階適用する」作用素の自然な候補として ( J α f ) ( x ) = d − α d x − α f ( x ) = 1 Γ ( α ) ∫ 0 x ( x − t ) α − 1 f ( t ) d t {\displaystyle (J^{\alpha }f)(x)={d^{-\alpha } \over dx^{-\alpha }}f(x)={1 \over \Gamma (\alpha )}\int _{0}^{x}(x-t)^{\alpha -1}f(t)\;dt} を与えることができる。これは実際に作用素として矛盾なく定まる。 作用素 J は可換かつ加法的である。つまり、 ( J α ) ( J β ) f = ( J β ) ( J α ) f = ( J α + β ) f = 1 Γ ( α + β ) ∫ 0 x ( x − t ) α + β − 1 f ( t ) d t {\displaystyle (J^{\alpha })(J^{\beta })f=(J^{\beta })(J^{\alpha })f=(J^{\alpha +\beta })f={1 \over \Gamma (\alpha +\beta )}\int _{0}^{x}(x-t)^{\alpha +\beta -1}f(t)\;dt} が成立する。この性質は、分数階微積分作用素の半群性と呼ばれる。残念ながら微分作用素 D に関しての同様の議論はこれよりももっと著しく複雑なものになってしまうが、それでも D が一般には可換にも加法的にもならないことは示すことができる。 ここで関数 f(x) として f ( x ) = x k {\displaystyle f(x)=x^{k}} という形の単項式を考える。この一階導函数は周知の如く f ′ ( x ) = d d x f ( x ) = k x k − 1 {\displaystyle f'(x)={\frac {d}{dx}}f(x)=kx^{k-1}} で与えられる。また、この二階導函数は f ″ ( x ) = d 2 d x 2 f ( x ) = k ( k − 1 ) x k − 2 {\displaystyle f''(x)={\frac {d^{2}}{dx^{2}}}f(x)=k(k-1)x^{k-2}} で与えられる。微分を繰り返せば一般に f ( a ) ( x ) = d a d x a x k = k P a x k − a = k ! ( k − a ) ! x k − a {\displaystyle f^{(a)}(x)={\frac {d^{a}}{dx^{a}}}x^{k}={_{k}P_{a}}x^{k-a}={\frac {k!}{(k-a)!}}x^{k-a}} を得る。ここで、階乗をガンマ関数に置き換えることにより d a d x a x k = k ! ( k − a ) ! x k − a = Γ ( k + 1 ) Γ ( k − a + 1 ) x k − a {\displaystyle {\frac {d^{a}}{dx^{a}}}x^{k}={\frac {k!}{(k-a)!}}x^{k-a}={\frac {\Gamma (k+1)}{\Gamma (k-a+1)}}x^{k-a}} が成り立つものと考えることができる。そのため、たとえば x の半微分（二分の一階導函数）は d 1 / 2 d x 1 / 2 x = Γ ( 1 + 1 ) Γ ( 1 − 1 / 2 + 1 ) x 1 − 1 / 2 = Γ ( 2 ) Γ ( 3 / 2 ) x 1 / 2 = 2 π 1 / 2 x − 1 / 2 = 2 π − 1 / 2 x 1 / 2 = 2 x π {\displaystyle {\frac {d^{1/2}}{dx^{1/2}}}x={\frac {\Gamma (1+1)}{\Gamma (1-1/2+1)}}x^{1-1/2}={\frac {\Gamma (2)}{\Gamma (3/2)}}x^{1/2}={\frac {2}{\pi ^{1/2}}}x^{-1/2}=2\pi ^{-1/2}x^{1/2}=2{\sqrt {\frac {x}{\pi }}}} で与えられる。これをもう一回行うと、 d 1 / 2 d x 1 / 2 2 π 1 / 2 x 1 / 2 = 2 π 1 / 2 Γ ( 1 / 2 + 1 ) Γ ( 1 / 2 − 1 / 2 + 1 ) x 1 / 2 − 1 / 2 = 2 π 1 / 2 Γ ( 3 / 2 ) Γ ( 1 ) x 0 = 2 π 1 / 2 π 1 / 2 2 1 = 1 {\displaystyle {\frac {d^{1/2}}{dx^{1/2}}}{\frac {2}{\pi ^{1/2}}}x^{1/2}={\frac {2}{\pi ^{1/2}}}{\frac {\Gamma (1/2+1)}{\Gamma (1/2-1/2+1)}}x^{1/2-1/2}={\frac {2}{\pi ^{1/2}}}{\frac {\Gamma (3/2)}{\Gamma (1)}}x^{0}={\frac {2}{\pi ^{1/2}}}{\frac {\frac {\pi ^{1/2}}{2}}{1}}=1} が得られる。これはすなわち、そもそも成り立って欲しかった性質である ( d 1 / 2 d x 1 / 2 d 1 / 2 d x 1 / 2 ) x = d d x x = 1 {\displaystyle \left({\frac {d^{1/2}}{dx^{1/2}}}{\frac {d^{1/2}}{dx^{1/2}}}\right)x={d \over dx}x=1} がきちんと満たされていることを意味している。ここで、上述のような微分作用素の拡張は、なにも実数冪のみに縛られるものではない。例えば (1 − i)-階導函数の (1 + i)-階導函数は二階微分を与えるものである。もちろん a が負の整数以外の負の値をとるならば適当な積分が与えられる(Γ(x+yi)の定義域はx>0)。 ラプラス変換に関する話として分数階微積分の問題を考えることもできる。 L { J f } ( s ) = L { ∫ 0 t f ( τ ) d τ } ( s ) = 1 s ( L { f } ) ( s ) , L { J 2 f } = 1 s ( L { J f } ) ( s ) = 1 s 2 ( L { f } ) ( s ) {\displaystyle {\begin{aligned}{\mathcal {L}}\left\{Jf\right\}(s)={\mathcal {L}}\left\{\int _{0}^{t}f(\tau )\,d\tau \right\}(s)={\frac {1}{s}}({\mathcal {L}}\left\{f\right\})(s),\\{\mathcal {L}}\left\{J^{2}f\right\}={\frac {1}{s}}({\mathcal {L}}\left\{Jf\right\})(s)={\frac {1}{s^{2}}}({\mathcal {L}}\left\{f\right\})(s)\end{aligned}}} などが成立することに注意して、ここでは J α f = L − 1 { s − α ( L { f } ) ( s ) } {\displaystyle J^{\alpha }f={\mathcal {L}}^{-1}\left\{s^{-\alpha }({\mathcal {L}}\{f\})(s)\right\}} が成り立つものと考える。例えば、 J α ( t k ) = L − 1 { Γ ( k + 1 ) s α + k + 1 } = Γ ( k + 1 ) Γ ( α + k + 1 ) t α + k {\displaystyle J^{\alpha }\left(t^{k}\right)={\mathcal {L}}^{-1}\left\{{\Gamma (k+1) \over s^{\alpha +k+1}}\right\}={\Gamma (k+1) \over \Gamma (\alpha +k+1)}t^{\alpha +k}} となることが期待される。実際に、畳み込み公式 L { f ∗ g } = ( L { f } ) ( L { g } ) {\displaystyle {\mathcal {L}}\{f*g\}=({\mathcal {L}}\{f\})({\mathcal {L}}\{g\})} が与えられれば（p(x) = xα−1 として） ( J α f ) ( t ) = 1 Γ ( α ) L − 1 { ( L { p } ) ( L { f } ) } = 1 Γ ( α ) ( p ∗ f ) = 1 Γ ( α ) ∫ 0 t p ( t − τ ) f ( τ ) d τ = 1 Γ ( α ) ∫ 0 t ( t − τ ) α − 1 f ( τ ) d τ {\displaystyle {\begin{aligned}(J^{\alpha }f)(t)&={\frac {1}{\Gamma (\alpha )}}{\mathcal {L}}^{-1}\left\{\left({\mathcal {L}}\{p\}\right)({\mathcal {L}}\{f\})\right\}\\&={\frac {1}{\Gamma (\alpha )}}(p*f)\\&={\frac {1}{\Gamma (\alpha )}}\int _{0}^{t}p(t-\tau )f(\tau )\,d\tau \\&={\frac {1}{\Gamma (\alpha )}}\int _{0}^{t}(t-\tau )^{\alpha -1}f(\tau )\,d\tau \end{aligned}}} となることが示される。これは上述のコーシーの公式に他ならない。 ラプラス変換が上手く計算できる函数は比較的少ないが、しかし分数階微分方程式を解くにあたってはしばしば有用である。 古典的な形での分数階微分積分学は、リーマン-リウヴィル積分によって与えられるもので、これは本質的には上で述べたような内容のものである。また、一定周期ごとに繰り返すという「境界条件」を課せば、周期関数に対する理論であるワイル微積分（英語版）が考えられる。これはフーリエ級数に対して定義され、一定のフーリエ係数が消えている（したがって単位円上の積分して 0になるような函数に適用できる）ことを要請する。 a D t − α f ( t ) = a I t α f ( t ) = 1 Γ ( α ) ∫ a t ( t − τ ) α − 1 f ( τ ) d τ {\displaystyle _{a}D_{t}^{-\alpha }f(t)={}_{a}I_{t}^{\alpha }f(t)={\frac {1}{\Gamma (\alpha )}}\int _{a}^{t}(t-\tau )^{\alpha -1}f(\tau )d\tau } 対して、グリュンバルト-レトニコフ微分（英語版） は積分の代わりに微分から始める理論である。 アダマール分数階積分は J. Hadamard によって導入され、次の式で与えられる。 a D t − α f ( t ) = 1 Γ ( α ) ∫ a t ( log ⁡ t τ ) α − 1 f ( τ ) d τ τ , t > a . {\displaystyle _{a}\mathbf {D} _{t}^{-\alpha }f(t)={\frac {1}{\Gamma (\alpha )}}\int _{a}^{t}\left(\log {\frac {t}{\tau }}\right)^{\alpha -1}f(\tau ){\frac {d\tau }{\tau }},\qquad t>a.} 関数解析学の文脈では、冪のみならずもっと一般に函数 f に対する作用素 f(D) についてスペクトル論の汎函数計算における研究がなされる。擬微分作用素の理論においても D の冪について考えることができる。この作用素は特異積分作用素（英語版） の例として得られる。また、古典理論の高次元への一般化はリースポテンシャルの理論と呼ばれる。したがって「分数階微分積分学」について論じるのに、いくつもの現代的理論を利用することができる。 特殊函数論で重要なErdélyi-Kober operatorも参照。