畳み込み級数

曖昧さ回避 畳み込み」とは異なります。

数学において、畳み込み級数(たたみこみきゅうすう、: telescoping series; 望遠鏡級数)は、各項からその近くの後続または先行する項と打ち消しあう部分をとりだして、次々に項が消えていくことで和が求まるような級数である。[1][2]。こうやって項を打ち消しあって和を求める方法は差分法 (method of differences) や和分法としても知られる。

たとえば、級数

n = 1 1 n ( n + 1 ) {\displaystyle \sum _{n=1}^{\infty }{\frac {1}{n(n+1)}}}

は、以下のように簡単になる:

n = 1 1 n ( n + 1 ) = n = 1 ( 1 n 1 n + 1 ) = ( 1 1 2 ) + ( 1 2 1 3 ) + = 1 + ( 1 2 + 1 2 ) + ( 1 3 + 1 3 ) + = 1. {\displaystyle {\begin{aligned}\sum _{n=1}^{\infty }{\frac {1}{n(n+1)}}&=\sum _{n=1}^{\infty }\left({\frac {1}{n}}-{\frac {1}{n+1}}\right)\\&=\left(1-{\frac {1}{2}}\right)+\left({\frac {1}{2}}-{\frac {1}{3}}\right)+\cdots \\&=1+\left(-{\frac {1}{2}}+{\frac {1}{2}}\right)+\left(-{\frac {1}{3}}+{\frac {1}{3}}\right)+\cdots =1.\end{aligned}}}

注意点

差分法のテクニックは便利だが、落とし穴もある。

0 = n = 1 0 = n = 1 ( 1 1 ) = 1 + n = 1 ( 1 + 1 ) = 1 {\displaystyle 0=\sum _{n=1}^{\infty }0=\sum _{n=1}^{\infty }(1-1)=1+\sum _{n=1}^{\infty }(-1+1)=1}

とするのは正しくない、個々の項が 0 に収斂でもしない限り、項を括りなおすのは有効ではないからである(グランディ級数の項を見よ)。こういった間違いを避けるには、まず N 項までの和を求めておいてから、そのあと N を無限大に飛ばした極限を計算する。

n = 1 N 1 n ( n + 1 ) = n = 1 N ( 1 n 1 n + 1 ) = ( 1 1 2 ) + ( 1 2 1 3 ) + + ( 1 N 1 N + 1 ) = 1 + ( 1 2 + 1 2 ) + ( 1 3 + 1 3 ) + + ( 1 N + 1 N ) 1 N + 1 = 1 1 N + 1 1   a s   N . {\displaystyle {\begin{aligned}\sum _{n=1}^{N}{\frac {1}{n(n+1)}}&{}=\sum _{n=1}^{N}\left({\frac {1}{n}}-{\frac {1}{n+1}}\right)\\&{}=\left(1-{\frac {1}{2}}\right)+\left({\frac {1}{2}}-{\frac {1}{3}}\right)+\cdots +\left({\frac {1}{N}}-{\frac {1}{N+1}}\right)\\&{}=1+\left(-{\frac {1}{2}}+{\frac {1}{2}}\right)+\left(-{\frac {1}{3}}+{\frac {1}{3}}\right)+\cdots +\left(-{\frac {1}{N}}+{\frac {1}{N}}\right)-{\frac {1}{N+1}}\\&{}=1-{\frac {1}{N+1}}\to 1\ \mathrm {as} \ N\to \infty .\end{aligned}}}

といった具合である。

込み入った例

  • 多くの三角関数は、連続する項の間で畳み込みを起こさせる、階差としての表示をも持つ。たとえば
    n = 1 N sin ( n ) = n = 1 N 1 2 csc ( 1 2 ) ( 2 sin ( 1 2 ) sin ( n ) ) = 1 2 csc ( 1 2 ) n = 1 N ( cos ( 2 n 1 2 ) cos ( 2 n + 1 2 ) ) = 1 2 csc ( 1 2 ) ( cos ( 1 2 ) cos ( 2 N + 1 2 ) ) {\displaystyle {\begin{aligned}\sum _{n=1}^{N}\sin \left(n\right)&{}=\sum _{n=1}^{N}{\frac {1}{2}}\csc \left({\frac {1}{2}}\right)\left(2\sin \left({\frac {1}{2}}\right)\sin \left(n\right)\right)\\&{}={\frac {1}{2}}\csc \left({\frac {1}{2}}\right)\sum _{n=1}^{N}\left(\cos \left({\frac {2n-1}{2}}\right)-\cos \left({\frac {2n+1}{2}}\right)\right)\\&{}={\frac {1}{2}}\csc \left({\frac {1}{2}}\right)\left(\cos \left({\frac {1}{2}}\right)-\cos \left({\frac {2N+1}{2}}\right)\right)\end{aligned}}}
    といったような具合である。
  • fg多項式でそれらの商が部分分数に分解されるような形の和
    n = 1 N f ( n ) g ( n ) {\displaystyle \sum _{n=1}^{N}{f(n) \over g(n)}}
    でも、今考えたい差分による和をとる方法が有効でないこともある。特に、
    n = 0 2 n + 3 ( n + 1 ) ( n + 2 ) = n = 0 ( 1 n + 1 + 1 n + 2 ) = ( 1 1 + 1 2 ) + ( 1 2 + 1 3 ) + ( 1 3 + 1 4 ) + + ( 1 n 1 + 1 n ) + ( 1 n + 1 n + 1 ) + ( 1 n + 1 + 1 n + 2 ) + = {\displaystyle {\begin{aligned}\sum _{n=0}^{\infty }{\frac {2n+3}{(n+1)(n+2)}}&{}=\sum _{n=0}^{\infty }\left({\frac {1}{n+1}}+{\frac {1}{n+2}}\right)\\&{}=\left({\frac {1}{1}}+{\frac {1}{2}}\right)+\left({\frac {1}{2}}+{\frac {1}{3}}\right)+\left({\frac {1}{3}}+{\frac {1}{4}}\right)+\cdots \\&{}\qquad \cdots +\left({\frac {1}{n-1}}+{\frac {1}{n}}\right)+\left({\frac {1}{n}}+{\frac {1}{n+1}}\right)+\left({\frac {1}{n+1}}+{\frac {1}{n+2}}\right)+\cdots \\&{}=\infty \end{aligned}}}
    のように項がぜんぜん打ち消しあわない場合には無力である。
  • k を正の整数とすれば
    n = 1 1 n ( n + k ) = H k k {\displaystyle \sum _{n=1}^{\infty }{\frac {1}{n(n+k)}}={\frac {H_{k}}{k}}}
    となる(1/(k − 1) 以降の項がすべて打ち消されて消える)。ここで Hkk-次調和数である。

確率論における応用

確率論において、ポアソン過程(英語版)はもっとも単純な場合には「発生」がランダムな時間で引き起こされるような確率過程で、次の発生までの待ち時間が無記憶(英語版)指数分布に従い、任意の時間間隔における発生回数は(期待値が時間間隔の長さに比例するという)ポアソン分布にしたがう。さて、Xt を時刻 t までの発生回数とし、Txx-番目の発生までの待ち時間として、確率変数 Tx確率密度関数を求めよう。ポアソン分布の確率質量関数を用いれば

Pr ( X t = x ) = ( λ t ) x e λ t x ! {\displaystyle \Pr(X_{t}=x)={\frac {(\lambda t)^{x}e^{-\lambda t}}{x!}}}

がわかる。ここで λ は任意の時間間隔 1 での発生回数の平均値である。 [Xt ≥ x] なる事象は [Txt] なる事象と同じことであり、したがって同一の確率を持つ。したがって、求める密度関数は

f ( t ) = d d t Pr ( T x t ) = d d t Pr ( X t x ) = d d t ( 1 Pr ( X t x 1 ) ) = d d t ( 1 u = 0 x 1 Pr ( X t = u ) ) = d d t ( 1 u = 0 x 1 ( λ t ) u e λ t u ! ) = λ e λ t e λ t u = 1 x 1 ( λ u t u 1 ( u 1 ) ! λ u + 1 t u u ! ) {\displaystyle {\begin{aligned}f(t)&{}={\frac {d}{dt}}\Pr(T_{x}\leq t)={\frac {d}{dt}}\Pr(X_{t}\geq x)={\frac {d}{dt}}(1-\Pr(X_{t}\leq x-1))\\[6pt]&{}={\frac {d}{dt}}\left(1-\sum _{u=0}^{x-1}\Pr(X_{t}=u)\right)={\frac {d}{dt}}\left(1-\sum _{u=0}^{x-1}{\frac {(\lambda t)^{u}e^{-\lambda t}}{u!}}\right)\\[6pt]&{}=\lambda e^{-\lambda t}-e^{-\lambda t}\sum _{u=1}^{x-1}\left({\frac {\lambda ^{u}t^{u-1}}{(u-1)!}}-{\frac {\lambda ^{u+1}t^{u}}{u!}}\right)\end{aligned}}}

であり、この和はほとんどが打ち消しあって、

f ( t ) = λ x t x 1 e λ t x ! {\displaystyle f(t)={\frac {\lambda ^{x}t^{x-1}e^{-\lambda t}}{x!}}}

だけが残る。

その他の応用

参考文献

  1. ^ Tom M. Apostol, Calculus, Volume 1, Blaisdell Publishing Company, 1962, pages 422–3
  2. ^ Brian S. Thomson and Andrew M. Bruckner, Elementary Real Analysis, Second Edition, CreateSpace, 2008, page 85