不定和分

数学における不定和分(ふていわぶん、: indefinite sumx または逆差分(ぎゃくさぶん、: antidifference; 反差分)Δ−1 [1][2][3] は、微分に対する不定積分(反微分)の離散版で、前進差分 Δ の逆演算となる線型作用素である。[注 1]

文献によっては "indefinite sum" の語を、例えば

k = 1 n f ( k ) {\displaystyle \sum _{k=1}^{n}f(k)}

のような和において、上の限界となる値 (この例では n) をとくに固定せずに考える場合を指すのに用いることもある。この場合、この和を表す閉じた式 F(n) は函数方程式(畳み込み方程式

F ( x + 1 ) F ( x ) = f ( x + 1 ) {\displaystyle F(x+1)-F(x)=f(x+1)}

の解[4]であり、これは後退差分作用素 の逆である。この後退和分作用素と先の(前進)和分作用素との間には後述の和分差分学の基本定理を通じて関係がある。

定義

与えられた函数 f(x) に対し F(x) := ∑xf(x)f(x)不定和分であるとは、F(x)函数方程式

Δ x f ( x ) = f ( x ) , {\displaystyle \Delta \sum _{x}f(x)=f(x),}

つまりより直接的に述べれば

F ( x + 1 ) F ( x ) = f ( x ) {\displaystyle F(x+1)-F(x)=f(x)}

の解であることを総称して言う。函数 F(x) が与えられた f(x) に対するこの函数方程式の解ならば、周期 1 を持つ任意の周期函数 C(x) に対して f(x) + C(x) もまた同方程式の解である[注 2]から、各不定和分とは実際にはそのような(互いに周期 1 函数だけ異なる)函数の族を表すものと理解される。ただし、解のうちで自身のニュートン級数(英語版)展開と一致するようなものは任意定数 C (和分定数)を加える違いを除いて一意に定まる。

注意

和分定数 C の選び方について、

F ( x ) = x f ( x ) + C {\displaystyle F(x)=\sum _{x}f(x)+C}

と置くとき、

0 1 F ( x ) d x = 0 {\displaystyle \int _{0}^{1}F(x)dx=0}

あるいは

1 2 F ( x ) d x = 0 {\displaystyle \int _{1}^{2}F(x)dx=0}

を満たすように C を固定することがしばしばある。ラマヌジャン和を用いて書けば、それぞれ

x 1 f ( x ) = f ( 0 ) F ( 0 ) {\displaystyle \sum _{x\geq 1}^{\Re }f(x)=-f(0)-F(0)}

あるいは

x 1 f ( x ) = F ( 1 ) {\displaystyle \sum _{x\geq 1}^{\Re }f(x)=-F(1)}

である[5][6]

性質

和分差分学の基本定理

微分積分学の基本定理の離散版として、不定和分を用いて定和分の計算ができる[7]。即ち

k = a b f ( k ) = Δ 1 f ( b + 1 ) Δ 1 f ( a ) {\displaystyle \sum _{k=a}^{b}f(k)=\Delta ^{-1}f(b+1)-\Delta ^{-1}f(a)}

が成り立つ。

部分和分法

詳細は「部分和分」を参照

部分積分法の離散版として、以下のように部分和分の公式が成り立つ。

不定和分に関する部分和分
x f ( x ) Δ g ( x ) = f ( x ) g ( x ) x ( g ( x ) + Δ g ( x ) ) Δ f ( x ) , {\displaystyle \sum _{x}f(x)\Delta g(x)=f(x)g(x)-\sum _{x}(g(x)+\Delta g(x))\Delta f(x),}
x f ( x ) Δ g ( x ) + x g ( x ) Δ f ( x ) = f ( x ) g ( x ) x Δ f ( x ) Δ g ( x ) . {\displaystyle \sum _{x}f(x)\Delta g(x)+\sum _{x}g(x)\Delta f(x)=f(x)g(x)-\sum _{x}\Delta f(x)\Delta g(x).}
定和分に関する部分和分
i = a b f ( i ) Δ g ( i ) = f ( b + 1 ) g ( b + 1 ) f ( a ) g ( a ) i = a b g ( i + 1 ) Δ f ( i ) . {\displaystyle \sum _{i=a}^{b}f(i)\Delta g(i)=f(b+1)g(b+1)-f(a)g(a)-\sum _{i=a}^{b}g(i+1)\Delta f(i).}

周期法則

周期函数 f(x) の周期が T に対し、

x f ( T x ) = x f ( T x ) + C {\displaystyle \sum _{x}f(Tx)=xf(Tx)+C}

が成り立つ。また T が函数 f(x) の反周期、即ち f(x + T) = −f(x) のとき、

x f ( T x ) = 1 2 f ( T x ) + C {\displaystyle \sum _{x}f(Tx)=-{\frac {1}{2}}f(Tx)+C}

が成り立つ。

ラプラス和公式

x f ( x ) = 0 x f ( t ) d t + k = 1 c k Δ k 1 f ( x ) k ! + C {\displaystyle \sum _{x}f(x)=\int _{0}^{x}f(t)dt+\sum _{k=1}^{\infty }{\frac {c_{k}\Delta ^{k-1}f(x)}{k!}}+C}

ただし、

c k = 0 1 Γ ( x + 1 ) Γ ( x k + 1 ) d x {\displaystyle c_{k}=\int _{0}^{1}{\frac {\Gamma (x+1)}{\Gamma (x-k+1)}}dx}

第二種ベルヌイ数(英語版)[8]である。

ニュートンの公式

x f ( x ) = k = 1 Δ k 1 f ( x ) k ! ( x ) k + C {\displaystyle \sum _{x}f(x)=-\sum _{k=1}^{\infty }{\frac {\Delta ^{k-1}f(x)}{k!}}(-x)_{k}+C}

ただし

( x ) k = Γ ( x + 1 ) Γ ( x k + 1 ) {\displaystyle (x)_{k}={\frac {\Gamma (x+1)}{\Gamma (x-k+1)}}}

下降階乗冪である。

ファウルハーバーの公式

x f ( x ) = n = 1 f ( n 1 ) ( 0 ) n ! B n ( x ) + C . {\displaystyle \sum _{x}f(x)=\sum _{n=1}^{\infty }{\frac {f^{(n-1)}(0)}{n!}}B_{n}(x)+C.}

ただし、右辺が存在する場合に限る。

ミューラーの公式

lim x + f ( x ) = 0 {\displaystyle \lim _{x\to {+\infty }}f(x)=0} のとき

x f ( x ) = n = 0 ( f ( n ) f ( n + x ) ) + C {\displaystyle \sum _{x}f(x)=\sum _{n=0}^{\infty }(f(n)-f(n+x))+C}

が成り立つ[9]

オイラー・マクローリンの公式

x f ( x ) = 0 x f ( t ) d t 1 2 f ( x ) + k = 1 B 2 k ( 2 k ) ! f ( 2 k 1 ) ( x ) + C . {\displaystyle \sum _{x}f(x)=\int _{0}^{x}f(t)dt-{\frac {1}{2}}f(x)+\sum _{k=1}^{\infty }{\frac {B_{2k}}{(2k)!}}f^{(2k-1)}(x)+C.}

以下、いくつかの函数の不定和分を記す。初等函数の不定和分でも必ずしも初等函数で書けないことに注意。

初等函数

  • x a = a x + C , {\displaystyle \sum _{x}a=ax+C,}
  • x x = x 2 2 x 2 + C , {\displaystyle \sum _{x}x={\frac {x^{2}}{2}}-{\frac {x}{2}}+C,}
  • x x a = B a + 1 ( x ) a + 1 + C , ( a Z ) , {\displaystyle \sum _{x}x^{a}={\frac {B_{a+1}(x)}{a+1}}+C,\quad (a\notin \mathbb {Z} ^{-}),}
    ただし、Ba(x) = −aζ(−a+1,x) は次数を実数に一般化したベルヌイ多項式.
  • x x a = ( 1 ) a 1 ψ ( a 1 ) ( x ) Γ ( a ) + C , a Z {\displaystyle \sum _{x}x^{a}={\frac {(-1)^{a-1}\psi ^{(-a-1)}(x)}{\Gamma (-a)}}+C,\,a\in \mathbb {Z} ^{-}}
    ただし、ψ(n)(x)ポリガンマ函数.
  • x 1 x = ψ ( x ) + C , {\displaystyle \sum _{x}{\frac {1}{x}}=\psi (x)+C,}
    ただし、ψ(x)ディガンマ函数.
  • x a x = a x a 1 + C , {\displaystyle \sum _{x}a^{x}={\frac {a^{x}}{a-1}}+C,}
    特に x 2 x = 2 x + C . {\displaystyle \sum _{x}2^{x}=2^{x}+C.}
  • x log b x = log b Γ ( x ) + C , {\displaystyle \sum _{x}\log _{b}x=\log _{b}\Gamma (x)+C,}
  • x log b a x = log b ( a x 1 Γ ( x ) ) + C . {\displaystyle \sum _{x}\log _{b}ax=\log _{b}(a^{x-1}\Gamma (x))+C.}

三角函数・双曲線函数

  • x sinh a x = 1 2 csch ( a 2 ) cosh ( a 2 a x ) + C , {\displaystyle \sum _{x}\sinh ax={\frac {1}{2}}\operatorname {csch} \left({\frac {a}{2}}\right)\cosh \left({\frac {a}{2}}-ax\right)+C,}
  • x cosh a x = 1 2 coth ( a 2 ) sinh a x 1 2 cosh a x + C , {\displaystyle \sum _{x}\cosh ax={\frac {1}{2}}\coth \left({\frac {a}{2}}\right)\sinh ax-{\frac {1}{2}}\cosh ax+C,}
  • x tanh a x = 1 a ψ e a ( x i π 2 a ) + 1 a ψ e a ( x + i π 2 a ) x + C , {\displaystyle \sum _{x}\tanh ax={\frac {1}{a}}\psi _{e^{a}}\left(x-{\frac {i\pi }{2a}}\right)+{\frac {1}{a}}\psi _{e^{a}}\left(x+{\frac {i\pi }{2a}}\right)-x+C,}
    ただし、ψq(x)q-ディガンマ函数.
  • x sin a x = 1 2 csc ( a 2 ) cos ( a 2 a x ) + C ( a n π ) , {\displaystyle \sum _{x}\sin ax=-{\frac {1}{2}}\csc \left({\frac {a}{2}}\right)\cos \left({\frac {a}{2}}-ax\right)+C\quad (a\neq n\pi ),}
  • x cos a x = 1 2 cot ( a 2 ) sin a x 1 2 cos a x + C ( a n π ) , {\displaystyle \sum _{x}\cos ax={\frac {1}{2}}\cot \left({\frac {a}{2}}\right)\sin ax-{\frac {1}{2}}\cos ax+C\quad (a\neq n\pi ),}
  • x sin 2 a x = x 2 + 1 4 csc ( a ) sin ( a 2 a x ) + C ( a n π 2 ) , {\displaystyle \sum _{x}\sin ^{2}ax={\frac {x}{2}}+{\frac {1}{4}}\csc(a)\sin(a-2ax)+C\quad (a\neq {\frac {n\pi }{2}}),}
  • x cos 2 a x = x 2 1 4 csc ( a ) sin ( a 2 a x ) + C ( a n π 2 ) , {\displaystyle \sum _{x}\cos ^{2}ax={\frac {x}{2}}-{\frac {1}{4}}\csc(a)\sin(a-2ax)+C\quad (a\neq {\frac {n\pi }{2}}),}
  • x tan a x = i x 1 a ψ e 2 i a ( x π 2 a ) + C ( a n π 2 ) , {\displaystyle \sum _{x}\tan ax=ix-{\frac {1}{a}}\psi _{e^{2ia}}\left(x-{\frac {\pi }{2a}}\right)+C\quad (a\neq {\frac {n\pi }{2}}),}
    ただし、ψq(x)q-ディガンマ函数.
  • x tan x = i x ψ e 2 i ( x + π 2 ) + C = k = 1 ( ψ ( k π π 2 + 1 z ) + ψ ( k π π 2 + z ) ψ ( k π π 2 + 1 ) ψ ( k π π 2 ) ) + C , {\displaystyle {\begin{aligned}\sum _{x}\tan x&=ix-\psi _{e^{2i}}\left(x+{\frac {\pi }{2}}\right)+C\\&=-\sum _{k=1}^{\infty }\left(\psi \left(k\pi -{\frac {\pi }{2}}+1-z\right)+\psi \left(k\pi -{\frac {\pi }{2}}+z\right)-\psi \left(k\pi -{\frac {\pi }{2}}+1\right)-\psi \left(k\pi -{\frac {\pi }{2}}\right)\right)+C\end{aligned}},}
  • x cot a x = i x i ψ e 2 i a ( x ) a + C ( a n π 2 ) , {\displaystyle \sum _{x}\cot ax=-ix-{\frac {i\psi _{e^{2ia}}(x)}{a}}+C\quad (a\neq {\frac {n\pi }{2}}),}
  • x artanh a x = 1 2 ln ( ( 1 ) x Γ ( 1 a ) Γ ( x + 1 a ) Γ ( 1 a ) Γ ( x 1 a ) ) + C , {\displaystyle \sum _{x}\operatorname {artanh} \,ax={\frac {1}{2}}\ln \left({\frac {(-1)^{x}\Gamma \left(-{\frac {1}{a}}\right)\Gamma \left(x+{\frac {1}{a}}\right)}{\Gamma \left({\frac {1}{a}}\right)\Gamma \left(x-{\frac {1}{a}}\right)}}\right)+C,}
  • x arctan a x = i 2 ln ( ( 1 ) x Γ ( i a ) Γ ( x + i a ) Γ ( i a ) Γ ( x i a ) ) + C . {\displaystyle \sum _{x}\arctan ax={\frac {i}{2}}\ln \left({\frac {(-1)^{x}\Gamma ({\frac {-i}{a}})\Gamma (x+{\frac {i}{a}})}{\Gamma ({\frac {i}{a}})\Gamma (x-{\frac {i}{a}})}}\right)+C.}

特殊函数

  • x ψ ( x ) = ( x 1 ) ψ ( x ) x + C , {\displaystyle \sum _{x}\psi (x)=(x-1)\psi (x)-x+C,}
  • x Γ ( x ) = ( 1 ) x + 1 Γ ( x ) Γ ( 1 x , 1 ) e + C , {\displaystyle \sum _{x}\Gamma (x)=(-1)^{x+1}\Gamma (x){\frac {\Gamma (1-x,-1)}{e}}+C,}
    ただし、Γ(s,x)不完全ガンマ函数.
  • x ( x ) a = ( x ) a + 1 a + 1 + C , {\displaystyle \sum _{x}(x)_{a}={\frac {(x)_{a+1}}{a+1}}+C,}
    ただし、(x)a下降階乗冪.
  • x sexp a ( x ) = ln a ( sexp a ( x ) ) ( ln a ) x + C . {\displaystyle \sum _{x}\operatorname {sexp} _{a}(x)=\ln _{a}{\frac {(\operatorname {sexp} _{a}(x))'}{(\ln a)^{x}}}+C.}
(ただし、sexp超指数函数(英語版))

関連項目

  1. ^ 一つのパラメータ h を導入して、歩み h の差分 Δhf(x) := f(x + h) − f(x) あるいは差分商 Δhf(x)Δhx = f(x+h)−f(x)h の逆演算として、歩み h の不定和分 f(xhx を考えることもある。h = 1 が本項における場合であり、また h→0 の極限で Δhf(x)Δhxdf(x)dx は微分商、f(xhx → ∫f(x)dx は不定積分となる。
  2. ^ 従って特に、函数 f(x) として定義域が整数全体となるようなもの、即ち数列 (an) を取るならば、周期 1 の周期函数とは定数列に他ならない。

参考文献

  1. ^ Indefinite Sum - PlanetMath.(英語)
  2. ^ On Computing Closed Forms for Indefinite Summations. Yiu-Kwong Man. J. Symbolic Computation (1993), 16, 355-376
  3. ^ "If Y is a function whose first difference is the function y, then Y is called an indefinite sum of y and denoted Δ−1y" Introduction to Difference Equations, Samuel Goldberg
  4. ^ Algorithms for Nonlinear Higher Order Difference Equations, Manuel Kauers
  5. ^ Bruce C. Berndt, Ramanujan's Notebooks Archived 2006年10月12日, at the Wayback Machine., Ramanujan's Theory of Divergent Series, Chapter 6, Springer-Verlag (ed.), (1939), pp. 133–149.
  6. ^ Éric Delabaere, Ramanujan's Summation, Algorithms Seminar 2001–2002, F. Chyzak (ed.), INRIA, (2003), pp. 83–88.
  7. ^ "Handbook of discrete and combinatorial mathematics", Kenneth H. Rosen, John G. Michaels, CRC Press, 1999, ISBN 0-8493-0149-1
  8. ^ Bernoulli numbers of the second kind on Mathworld
  9. ^ Markus Müller. How to Add a Non-Integer Number of Terms, and How to Produce Unusual Infinite Summations (note that he uses a slightly alternative definition of fractional sum in his work, i.e. inverse to backwards difference, hence 1 as the lower limit in his formula)

関連文献

  • "Difference Equations: An Introduction with Applications", Walter G. Kelley, Allan C. Peterson, Academic Press, 2001, ISBN 0-12-403330-X
  • Markus Müller. How to Add a Non-Integer Number of Terms, and How to Produce Unusual Infinite Summations
  • Markus Mueller, Dierk Schleicher. Fractional Sums and Euler-like Identities
  • S. P. Polyakov. Indefinite summation of rational functions with additional minimization of the summable part. Programmirovanie, 2008, Vol. 34, No. 2.
  • "Finite-Difference Equations And Simulations", Francis B. Hildebrand, Prenctice-Hall, 1968