離散計算 では、 不定和 演算子( 反差分 演算子とも呼ばれる )は、 または [1] [2] で表され、 線形演算子 であり 、 前進差分演算子の逆である。不定積分が 導関数 と関係しているの と同様に、前進差分演算子と 不定積分 は関係している。したがって
、
∑
x
{\textstyle \sum _{x}}
Δ
−
1
{\displaystyle \Delta ^{-1}}
Δ
{\displaystyle \Delta }
Δ
∑
x
ふ
(
x
)
=
ふ
(
x
)
。
{\displaystyle \Delta \sum _{x}f(x)=f(x)\,.}
より明確に言えば、 ならば、
∑
x
ふ
(
x
)
=
ふ
(
x
)
{\textstyle \sum _{x}f(x)=F(x)}
ふ
(
x
+
1
)
−
ふ
(
x
)
=
ふ
(
x
)
。
{\displaystyle F(x+1)-F(x)=f(x)\,.}
F ( x ) が与えられた f ( x )に対するこの関数方程式の解である 場合 、周期 1 の 任意の周期関数 C ( x ) に対して F ( x )+ C ( x ) も解になります。したがって、各不定和は実際には関数の族を表します。ただし、 カールソンの定理 により、その ニュートン級数 展開に等しい解は、加法定数 C を除いて一意です。この一意の解は、 反差分演算子の形式的なべき 級数 形式で表すことができます 。
Δ
−
1
=
1
e
だ
−
1
{\displaystyle \Delta ^{-1}={\frac {1}{e^{D}-1}}}
離散積分学の基本定理
不定和は次の式で定和を計算するために使用できる: [3]
∑
け
=
1つの
b
ふ
(
け
)
=
Δ
−
1
ふ
(
b
+
1
)
−
Δ
−
1
ふ
(
1つの
)
{\displaystyle \sum _{k=a}^{b}f(k)=\Delta ^{-1}f(b+1)-\Delta ^{-1}f(a)}
定義
ラプラス の総和公式は、 不定和を 不定積分と 差分演算子の 反復から得られる補正項の和として表すことを可能にするが、もともとは積分を不定和と補正項の和として表す逆のプロセスのために開発された。不定和や不定積分の場合と同様に、 積分定数 を任意に選択するまで有効である 。 演算子代数 を使用すると、演算対象となる関数を繰り返しコピーして式が乱雑になるのを防ぐことができる。 [4]
∑
x
=
∫
+
1
2
−
1
12
Δ
+
1
24
Δ
2
−
19
720
Δ
3
+
3
160
Δ
4
−
⋯
{\displaystyle \sum _{x}=\int {}+{\frac {1}{2}}-{\frac {1}{12}}\Delta +{\frac {1}{24}}\Delta ^{2}-{\frac {19}{720}}\Delta ^{3}+{\frac {3}{160}}\Delta ^{4}-\cdots }
例えば、この式では、項は 与えられた関数を2で割る演算子を表します。 この式に現れる係数、、などはグレゴリー係数であり 、 ラプラス 数とも呼ばれます。項の係数は [4] です。
1
2
{\displaystyle {\tfrac {1}{2}}}
+
1
2
{\displaystyle +{\tfrac {1}{2}}}
−
1
12
{\displaystyle -{\tfrac {1}{12}}}
Δ
n
−
1
{\displaystyle \Delta ^{n-1}}
C
n
n
!
=
∫
0
1
(
x
n
)
d
x
{\displaystyle {\frac {{\mathcal {C}}_{n}}{n!}}=\int _{0}^{1}{\binom {x}{n}}\,dx}
ここで、左辺の 分子は第一種コーシー数と呼ばれるが、この名前はグレゴリー係数自体にも適用されることがある。 [4]
C
n
{\displaystyle {\mathcal {C}}_{n}}
∑
x
f
(
x
)
=
∑
k
=
1
∞
(
x
k
)
Δ
k
−
1
[
f
]
(
0
)
+
C
=
∑
k
=
1
∞
Δ
k
−
1
[
f
]
(
0
)
k
!
(
x
)
k
+
C
{\displaystyle \sum _{x}f(x)=\sum _{k=1}^{\infty }{\binom {x}{k}}\Delta ^{k-1}[f]\left(0\right)+C=\sum _{k=1}^{\infty }{\frac {\Delta ^{k-1}[f](0)}{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\,,}
ファウルハーバーの公式は、方程式の右辺が収束することを規定しています。
もし そうなら [5]
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 }\left(f(n)-f(n+x)\right)+C.}
∑
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}
定数項の選択
多くの場合、不定和の
定数 C は次の条件から決定されます。
させて
F
(
x
)
=
∑
x
f
(
x
)
+
C
{\displaystyle F(x)=\sum _{x}f(x)+C}
定数 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}
あるいは、ラマヌジャンの和を使用することもできます。
∑
x
≥
1
ℜ
f
(
x
)
=
−
f
(
0
)
−
F
(
0
)
{\displaystyle \sum _{x\geq 1}^{\Re }f(x)=-f(0)-F(0)}
または1
∑
x
≥
1
ℜ
f
(
x
)
=
−
F
(
1
)
{\displaystyle \sum _{x\geq 1}^{\Re }f(x)=-F(1)}
それぞれ [6] [7]
部分ごとの合計
部分ごとの不定合計:
∑
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)}
期間ルール
関数の周期で あれ
ば
T
{\displaystyle T}
f
(
x
)
{\displaystyle f(x)}
∑
x
f
(
T
x
)
=
x
f
(
T
x
)
+
C
{\displaystyle \sum _{x}f(Tx)=xf(Tx)+C}
が関数の反周期で ある
場合 、
T
{\displaystyle T}
f
(
x
)
{\displaystyle f(x)}
f
(
x
+
T
)
=
−
f
(
x
)
{\displaystyle 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}
代替使用法
上限の数値が指定されていない合計を説明するために、「不定の合計」というフレーズを使用する著者もいます。
∑
k
=
1
n
f
(
k
)
.
{\displaystyle \sum _{k=1}^{n}f(k).}
この場合、 和の
閉じた形式表現 F ( k )は、
F
(
x
+
1
)
−
F
(
x
)
=
f
(
x
+
1
)
{\displaystyle F(x+1)-F(x)=f(x+1)}
これはテレスコーピング方程式と呼ばれます。 [8]これは 後方差分 演算子の逆です 。これは、前述の離散計算の基本定理を使用して前方反差分演算子に関連付けられます。
∇
{\displaystyle \nabla }
不確定額の一覧
これはさまざまな関数の不定和のリストです。すべての関数が基本関数で表現できる不定和を持つわけではありません。
有理関数の反差分
∑
x
a
=
a
x
+
C
{\displaystyle \sum _{x}a=ax+C}
ここから 因数分解すると 1 が残り、別の形式 となります 。そこから次の式が得られます。
a
{\displaystyle a}
x
0
{\displaystyle x^{0}}
∑
x
x
0
=
x
{\displaystyle \sum _{x}x^{0}=\ x}
以下の合計については、
x
=
x
1
{\displaystyle x=x^{1}}
∑
x
x
=
x
(
x
+
1
)
2
+
C
{\displaystyle \sum _{x}x={\frac {x(x+1)}{2}}+C}
正の整数指数には ファウルハーバーの公式 が使える。負の整数指数には
∑
x
1
x
a
=
(
−
1
)
a
+
1
ψ
(
a
+
1
)
(
x
)
a
!
+
C
,
a
∈
Z
{\displaystyle \sum _{x}{\frac {1}{x^{a}}}={\frac {(-1)^{a+1}\psi ^{(a+1)}(x)}{a!}}+C,\,a\in \mathbb {Z} }
ここで ポリガンマ 関数 を使用できます。
ψ
(
n
)
(
x
)
{\displaystyle \psi ^{(n)}(x)}
より一般的には、
∑
x
x
a
=
{
−
ζ
(
−
a
,
x
+
1
)
+
C
1
,
if
a
≠
−
1
ψ
(
x
+
1
)
+
C
2
,
if
a
=
−
1
{\displaystyle \sum _{x}x^{a}={\begin{cases}-\zeta (-a,x+1)+C_{1},&{\text{if }}a\neq -1\\\psi (x+1)+C_{2},&{\text{if }}a=-1\end{cases}}}
ここで、 は フルヴィッツゼータ関数 、は ディガンマ関数 です 。 および は通常それぞれ (ここで 、 は リーマンゼータ関数 ) および オイラー・マスケローニ定数 に設定される定数です 。 変数 を に置き換えると、これは 一般化調和数 になります。 フルヴィッツゼータ関数 と ポリガンマ 関数の関係については、 バランスポリガンマ関数 と フルヴィッツゼータ関数#特殊なケースと一般化 を参照してください 。
ζ
(
s
,
a
)
{\displaystyle \zeta (s,a)}
ψ
(
z
)
{\displaystyle \psi (z)}
C
1
{\displaystyle C_{1}}
C
2
{\displaystyle C_{2}}
ζ
(
−
a
)
{\displaystyle \zeta (-a)}
ζ
(
s
)
{\displaystyle \zeta (s)}
a
{\displaystyle a}
−
a
{\displaystyle -a}
これから、 を使用すると 、別の形式が得られます。
∂
∂
a
ζ
(
s
,
a
)
=
−
s
ζ
(
s
+
1
,
a
)
{\displaystyle {\frac {\partial }{\partial a}}\zeta (s,a)=-s\zeta (s+1,a)}
∑
x
x
a
=
∫
0
x
−
a
ζ
(
1
−
a
,
u
+
1
)
d
u
+
C
,
if
a
≠
−
1
{\displaystyle \sum _{x}x^{a}=\int _{0}^{x}-a\zeta (1-a,u+1)du+C,{\text{ if }}a\neq -1}
∑
x
B
a
(
x
)
=
(
x
−
1
)
B
a
(
x
)
−
a
a
+
1
B
a
+
1
(
x
)
+
C
{\displaystyle \sum _{x}B_{a}(x)=(x-1)B_{a}(x)-{\frac {a}{a+1}}B_{a+1}(x)+C}
指数関数の反差分
∑
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}(x!)+C}
∑
x
log
b
a
x
=
log
b
(
x
!
a
x
)
+
C
{\displaystyle \sum _{x}\log _{b}ax=\log _{b}(x!a^{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
csch
(
a
2
)
sinh
(
a
x
−
a
2
)
+
C
{\displaystyle \sum _{x}\cosh ax={\frac {1}{2}}\operatorname {csch} \left({\frac {a}{2}}\right)\sinh \left(ax-{\frac {a}{2}}\right)+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-ディガンマ 関数です 。
ψ
q
(
x
)
{\displaystyle \psi _{q}(x)}
三角関数の反差
∑
x
sin
a
x
=
−
1
2
csc
(
a
2
)
cos
(
a
2
−
a
x
)
+
C
,
a
≠
2
n
π
{\displaystyle \sum _{x}\sin ax=-{\frac {1}{2}}\csc \left({\frac {a}{2}}\right)\cos \left({\frac {a}{2}}-ax\right)+C\,,\,\,a\neq 2n\pi }
∑
x
cos
a
x
=
1
2
csc
(
a
2
)
sin
(
a
x
−
a
2
)
+
C
,
a
≠
2
n
π
{\displaystyle \sum _{x}\cos ax={\frac {1}{2}}\csc \left({\frac {a}{2}}\right)\sin \left(ax-{\frac {a}{2}}\right)+C\,,\,\,a\neq 2n\pi }
∑
x
sin
2
a
x
=
x
2
+
1
4
csc
(
a
)
sin
(
a
−
2
a
x
)
+
C
,
a
≠
n
π
{\displaystyle \sum _{x}\sin ^{2}ax={\frac {x}{2}}+{\frac {1}{4}}\csc(a)\sin(a-2ax)+C\,\,,\,\,a\neq n\pi }
∑
x
cos
2
a
x
=
x
2
−
1
4
csc
(
a
)
sin
(
a
−
2
a
x
)
+
C
,
a
≠
n
π
{\displaystyle \sum _{x}\cos ^{2}ax={\frac {x}{2}}-{\frac {1}{4}}\csc(a)\sin(a-2ax)+C\,\,,\,\,a\neq n\pi }
∑
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\,,\,\,a\neq {\frac {n\pi }{2}}}
ここで q-ディガンマ 関数です 。
ψ
q
(
x
)
{\displaystyle \psi _{q}(x)}
∑
x
tan
x
=
i
x
−
ψ
e
2
i
(
x
+
π
2
)
+
C
=
−
∑
k
=
1
∞
(
ψ
(
k
π
−
π
2
+
1
−
x
)
+
ψ
(
k
π
−
π
2
+
x
)
−
ψ
(
k
π
−
π
2
+
1
)
−
ψ
(
k
π
−
π
2
)
)
+
C
{\displaystyle \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-x\right)+\psi \left(k\pi -{\frac {\pi }{2}}+x\right)-\psi \left(k\pi -{\frac {\pi }{2}}+1\right)-\psi \left(k\pi -{\frac {\pi }{2}}\right)\right)+C}
∑
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\,,\,\,a\neq {\frac {n\pi }{2}}}
∑
x
sinc
x
=
sinc
(
x
−
1
)
(
1
2
+
(
x
−
1
)
(
ln
(
2
)
+
ψ
(
x
−
1
2
)
+
ψ
(
1
−
x
2
)
2
−
ψ
(
x
−
1
)
+
ψ
(
1
−
x
)
2
)
)
+
C
{\displaystyle \sum _{x}\operatorname {sinc} x=\operatorname {sinc} (x-1)\left({\frac {1}{2}}+(x-1)\left(\ln(2)+{\frac {\psi ({\frac {x-1}{2}})+\psi ({\frac {1-x}{2}})}{2}}-{\frac {\psi (x-1)+\psi (1-x)}{2}}\right)\right)+C}
ここで、正規化された sinc関数 です 。
sinc
(
x
)
{\displaystyle \operatorname {sinc} (x)}
逆双曲関数の反差分
∑
x
artanh
a
x
=
1
2
ln
(
Γ
(
x
+
1
a
)
Γ
(
x
−
1
a
)
)
+
C
{\displaystyle \sum _{x}\operatorname {artanh} \,ax={\frac {1}{2}}\ln \left({\frac {\Gamma \left(x+{\frac {1}{a}}\right)}{\Gamma \left(x-{\frac {1}{a}}\right)}}\right)+C}
逆三角関数の逆差分
∑
x
arctan
a
x
=
i
2
ln
(
Γ
(
x
+
i
a
)
Γ
(
x
−
i
a
)
)
+
C
{\displaystyle \sum _{x}\arctan ax={\frac {i}{2}}\ln \left({\frac {\Gamma (x+{\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
)
{\displaystyle \Gamma (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
{\displaystyle (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}
( 超指数関数を 参照)
参照
参考文献
^ Man, Yiu-Kwong (1993)、「不定和の閉形式の計算について」、 Journal of Symbolic Computation 、 16 (4): 355–376、 doi :10.1006/jsco.1993.1053、 MR 1263873
^ ゴールドバーグ、サミュエル(1958)、経済学、心理学、社会学の実例付き差分方程式入門、ワイリー、ニューヨーク、チャップマン&ホール、ロンドン、p. 41、 ISBN 978-0-486-65084-5 , MR 0094249, が関数 の最初の差分が関数 である 場合 、 は の 不定和と呼ばれ 、 と表記されます。
Y
{\displaystyle Y}
y
{\displaystyle y}
Y
{\displaystyle Y}
y
{\displaystyle y}
Δ
−
1
y
{\displaystyle \Delta ^{-1}y}
; 1986年にドーバーブックスから再版
^ 「離散数学と組合せ数学のハンドブック」、ケネス・H・ローゼン、ジョン・G・マイケルズ、CRC プレス、1999 年、 ISBN 0-8493-0149-1
^ abc メルリーニ、ドナテッラ;スプルニョーリ、レンツォ。 Verri, M. Cecilia (2006)、「The Cauchy Numbers」、 離散数学 、 306 (16): 1906–1920、 doi :10.1016/j.disc.2006.03.065、 MR 2251571
^ Markus Müller. 非整数項の加算方法と、異常な無限和を生成する方法 2011-06-17 に Wayback Machine でアーカイブ(彼の研究では分数和の定義が少し異なっており、つまり逆差分の逆数であるため、彼の公式の下限は 1 となっていることに注意してください)
^ Bruce C. Berndt、「Ramanujan's Notebooks Archived 2006-10-12 at the Wayback Machine 」、 「Ramanujan's Theory of Divergent Series」 、第6章、Springer-Verlag (ed.)、(1939)、pp. 133–149。
^ Éric Delabaere、『ラマヌジャンの総和、 アルゴリズム セミナー 2001 ~ 2002』 、F. Chyzak (編)、INRIA、(2003)、83 ~ 88 ページ。
^ 非線形高階差分方程式のアルゴリズム、マヌエル・カウアーズ
さらに読む
「差分方程式: 応用入門」、Walter G. Kelley、Allan C. Peterson、Academic Press、2001 年、 ISBN 0-12-403330-X
マルクス・ミュラー。非整数項の加算方法と、異常な無限和の生成方法
マルクス・ミュラー、ディーク・シュライヒャー。分数和とオイラーのような恒等式
SP Polyakov. 有理関数の不定和と加算可能部分の最小化の追加。Programmirovanie、2008 年、第 34 巻、第 2 号。
「有限差分方程式とシミュレーション」、フランシス・B・ヒルデブランド、プレンクティス・ホール、1968年