分数階微積分学

分数階微分積分学（ぶんすうかいびぶんせきぶんがく、英: 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も参照。 [脚注の使い方] ^ ここで積分作用素の J は integration の頭文字 I を用いるところ、I は恒等写像など他の意味に使われたり、I に似た字形の記号・文字がいろいろと使われたりすることによる混同を避けるためにしばしば使われる。 ^ Liouville, Joseph (1832), “Mémoire sur quelques questions de géométrie et de mécanique, et sur un nouveau genre de calcul pour résoudre ces questions”, Journal de l'École Polytechnique (Paris) 13: 1-69, https://gallica.bnf.fr/ark:/12148/bpt6k4336778/f2.item.r=Joseph%20Liouville . ^ Liouville, Joseph (1832), “Mémoire sur le calcul des différentielles à indices quelconques”, Journal de l'École Polytechnique (Paris) 13: 71-162, https://gallica.bnf.fr/ark:/12148/bpt6k4336778/f72.image . ^ この主題の歴史については、以下の修士論文（フランス語）Stéphane Dugowson, Les différentielles métaphysiques (histoire et philosophie de la généralisation de l'ordre de dérivation), Thèse, Université Paris Nord (1994) を参照。 ^ Hadamard, J. (1892), “Essai sur l'étude des fonctions données par leur développement de Taylor”, Journal of pure and applied mathematics 4 (8): 101–186, https://eudml.org/doc/233965  Fractional Integrals and Derivatives: Theory and Applications, by Samko, S.; Kilbas, A.A.; and Marichev, O. Hardcover: 1006 pages. Publisher: Taylor & Francis Books. ISBN 2-88124-864-0 Theory and Applications of Fractional Differential Equations, by Kilbas, A. A.; Srivastava, H. M.; and Trujillo, J. J. Amsterdam, Netherlands, Elsevier, February 2006. ISBN 0-444-51832-0 (http://www.elsevier.com/wps/find/bookdescription.cws_home/707212/description#description) An Introduction to the Fractional Calculus and Fractional Differential Equations, by Kenneth S. Miller, Bertram Ross (Editor). Hardcover: 384 pages. Publisher: John Wiley & Sons; 1 edition (May 19, 1993). ISBN 0-471-58884-9 The Fractional Calculus; Theory and Applications of Differentiation and Integration to Arbitrary Order (Mathematics in Science and Engineering, V), by Keith B. Oldham, Jerome Spanier. Hardcover. Publisher: Academic Press; (November 1974). ISBN 0-12-525550-0 Fractional Differential Equations. An Introduction to Fractional Derivatives, Fractional Differential Equations, Some Methods of Their Solution and Some of Their Applications., (Mathematics in Science and Engineering, vol. 198), by Igor Podlubny. Hardcover. Publisher: Academic Press; (October 1998) ISBN 0-12-558840-2 Fractional Calculus - An Introduction for Physicists by R. Herrmann, World Scientific, Singapore 2014. 500 pages. Fractals and quantum mechanics, by N. Laskin. Chaos Vol.10, pp.780-790 (2000). (http://link.aip.org/link/?CHAOEH/10/780/1) Fractals and Fractional Calculus in Continuum Mechanics, by A. Carpinteri (Editor), F. Mainardi (Editor). Paperback: 348 pages. Publisher: Springer-Verlag Telos; (January 1998). ISBN 3-211-82913-X Physics of Fractal Operators, by Bruce J. West, Mauro Bologna, Paolo Grigolini. Hardcover: 368 pages. Publisher: Springer Verlag; (January 14, 2003). ISBN 0-387-95554-2 Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models. by F. Mainardi, Imperial College Press, 2010. 368 pages. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media. by V.E. Tarasov, Springer, 2010. 450 pages. Fractional Derivatives for Physicists and Engineers by V.V. Uchaikin, Springer, Higher Education Press, 2012, 385 pages. Fractional Calculus and the Taylor-Riemann Series, Rose-Hulman Undergrad. J. Math. Vol.6(1) (2005). Operator of fractional derivative in the complex plane, by Petr Zavada, Commun.Math.Phys.192, pp. 261-285,1998. doi:10.1007/s002200050299 (available online or as the arXiv preprint) Relativistic wave equations with fractional derivatives and pseudodifferential operators, by Petr Zavada, Journal of Applied Mathematics, vol. 2, no. 4, pp. 163-197, 2002. doi:10.1155/S1110757X02110102 (available online or as the arXiv preprint) Fractional differentiation by neocortical pyramidal neurons, by Brian N Lundstrom, Matthew H Higgs, William J Spain & Adrienne L Fairhall, Nature Neuroscience, vol. 11 (11), pp. 1335 - 1342, 2008. doi:10.1038/nn.2212 (abstract) Equilibrium points, stability and numerical solutions of fractional-order predator-prey and rabies models, by Ahmed E., A.M.A. El-Sayed, H.A.A. El-Saka. 2007. Jour. Math. Anal. Appl. 325,452. 微分積分学 微積分作用素（英語版） 分数階微分方程式 函数的平方根 プロジェクト 数学 ポータル 数学 整数階でない微積分法について 非整数階の微積分 - ウェイバックマシン（2018年3月29日アーカイブ分） Weisstein, Eric W. "Fractional Calculus". mathworld.wolfram.com (英語). Weisstein, Eric W. "Fractional Derivative". mathworld.wolfram.com (英語). Weisstein, Eric W. "Fractional Differential Equation". mathworld.wolfram.com (英語). Fractional Calculus at MathPages Specialized journal: Fractional Calculus and Applied Analysis Initialized Fractional Calculus FRACTIONAL CALCULUS PROJECT Igor Podlubny's collection of related books, articles, links, software, etc. Les Différentielles métaphysiques――Histoire et philosophie de la généralisation de l'ordre de dérivation History, Definitions, and Applications for the Engineer (PDF) , by Adam Loverro, University of Notre Dame Fractional Calculus Modelling Introductory Notes on Fractional Calculus Pseudodifferential operators and diffusive representation in modeling, control and signal