2

零数列による級数の加速

75
0

あいさつ

んちゃ!

目次
  1. 零集合
  2. Euler-Maclaurin展開
  3. 零数列の構成
  4. Stirlingの公式
  5. 応用1
  6. 超幾何関数に関する公式
  7. 応用2
  8. 最後に

表記
  1. N0:={0}N
  2. Zero({an}nN0):={{wn}nN0C|nN0wnan=0}
  3. limNSN({an}nN0)=S({an}nN0)

零集合

零集合

数列{an}nN0Cに対して定まる以下の様な数列の集合を零集合と呼ぶ。またその元を零数列と呼ぶ。
Zero({an}nN0):={{wn}nN0C|nN0wnan=0}

零集合Zero({an}nN0)に次の様な和、スカラー積を導入すると零集合はベクトル空間を成す。
{{wn}nN0,{w}nN0Zero({an}nN0):{wn}nN0+{w}nN0={wn+wn}nN0cC:{wn}nN0Zero({an}nN0):c{wn}nN0={cwn}nN0

[0]和、スカラー倍は閉じている。
{wn}nN0,{w}nN0Zero({an}nN0):n=0(wn+wn)an=n=0wnan+n=0wnan=0+0=0
また
cC:{wn}nN0Zero({an}nN0):n=0(cwn)an=cn=0wnan=c0=0
[1]結合法則
{wn}nN0,{w}nN0,{w}nN0Zero({an}nN0):n=0(wn+(wn+wn))an=n=0((wn+wn)+wn)an
[2]単位元:{0}Zero({an}nN0)
{wn}nN0Zero({an}nN0):n=0(wn+0)an=n=0(0+wn)an=n=0wnan
[3]逆元:{wn}nN0Zero({an}nN0)
{wn}nN0Zero({an}nN0):n=0(wn+(wn))an=n=0(wn+wn)an=n=00an=0
[4]可換性
{wn}nN0,{w}nN0Zero({an}nN0):n=0(wn+wn)an=n=0(wn+wn)an
[5]ベクトルに関する分配法則
cC:{wn}nN0,{w}nN0Zero({an}nN0):n=0(c(wn+wn))an=n=0(cwn+cwn)an
[6]スカラーに関する分配法則
c,dC:{wn}nN0Zero({an}nN0):n=0((c+d)wn)an=n=0(cwn+dwn)an
[7]スカラー積に関する結合法則
c,dC:{wn}nN0Zero({an}nN0):n=0((cd)wn)an=n=0(c(dwn))an
[8]スカラー積の単位元:1C
{wn}nN0Zero({an}nN0):n=0(1wn)an=n=0wnan

数列{an}nN0Cに対して定まる総和S({an}nN0):=n=0anは下記の性質を持つ。
{wn}nN0Zero({an}nN0):S({an}nN0)=n=0(1+wn)an

そのまま計算するだけ。
n=0(1+wn)an=n=0an+n=0wnan=S({an}nN0)+0=S({an}nN0)

D次の加速度を持つ零数列

数列{an}nN0に対して定まる総和SN({an}nN0):=n=0N1anに対して次の性質を満たす数列{wn}nN0Zero({an}nN0)D次の加速度を持つ零数列と呼ぶ。
AR+ s.t. |SN({an}nN0)n=0N1wnanS({an}nN0)|<AND

漸近表示

数列{an}nN0に対して定まる総和SN({an}nN0):=n=0N1anに対して定まる次の性質を満たすNに関する多項式が存在する時、D次の漸近表示と呼ぶ。
A0,A1,...,AD1 s.t. SN({an}nN0)=S({an}nN0)+n=0D1AnNn+O(ND)

Euler-Maclaurin展開

Euler-Maclaulinの和公式導出
やりたい事:Smn=k=mnf(k)を積分で近似したい。
Euler-Maclaurinの公式 Euler-Maclaurinの公式
[1]まずはとりあえず、総和と積分の差分でも計算してみましょう。
ϵmn=Smnf(n)mnf(x)dx=k=mn1f(k)k=mn1kk+1f(x)dx=k=mn1{f(k)01f(x+k)dx}
[2]区間[0,1]で定義された微分可能な関数ψ(x)を考えて次の条件が成り立つ場合を考えてみる。
01ψ(x)f(x+k)dx=ψ(1)f(k+1)ψ(0)f(k)01ψ(x)f(x+k)dx(ψ(0)=ψ(1)=12,ψ(x)=1ψ(x)=x12)=12{f(k)+f(k+1)}01f(x+k)dxこの事を用いると。望遠鏡和を使う事で下記の式を得る。
{ψ(x)=x12f(n)f(m)2+k=mn101ψ(x)f(x+k)dx=Smnf(n)mnf(x)dx
[3]次に以下の様な記号を定めてみよう。
{RmnN=k=mn101BN(x)f(N)(x+k)(N=1,2,...)BN(x)=BN1(x)(N=2,3,...)BN(0)=BN(1)=bN1(N=2,3,...)B1(x)=ψ(x)
A.Rmn(N1),RmnNとの間に成り立つ関係式
RmnN=k=mn101BN(x)f(N)(x+k)=k=mn1[BN(x)f(N1)(x+k)]0101BN(x)f(N1)(x)dx=bN1{f(N1)(n)f(N1)(m)}01BN1(x)f(N1)(x)dx=bN1{f(N1)(n)f(N1)(m)}Rmn(N1)=bN1{f(N1)(n)f(N1)(m)}bN2{f(N2)(n)f(N2)(m)}++(1)N2b1{f(1)(n)f(1)(m)}+(1)N1Rmn1以上の計算から
Rmn1=(1)NbN1{f(N1)(n)f(N1)(m)}+(1)N1bN2{f(N2)(n)f(N2)(m)}++b1{f(1)(n)f(1)(m)}+(1)N1RmnN=k=1N1(1)k+1bk{fk(n)f(k)(m)}+(1)N1RmnN
B.BN(x)を具体的に求める。
(i)まず次のような母関数
f(x,y)=1+N=1BN(x)yNを定める。
(ii)両辺をxで偏微分する。
xf(x,y)=N=1BN(x)yN=y{1+N=2BN1(x)yN1}=y{1+N=1BN(x)yn}=yf(x,y)この計算から、あるyに関する関数を用いて以下の結果を得る。
f(x,y)=g(y)eyx
(iii)さらに
{f(0,y)=112y+N=2bN1yNf(1,y)=1+12y+N=2bN1yNを用いると次式が得られる。
g(y)+12y=g(y)ey12yこの式を用いる事で下記の式を得る。
g(y)=yey1
(iv)以上の結果をまとめる事で
f(x,y)=yexyey1
(v)また以上の計算により{bN}の母関数は下記の様に書ける事が分かる。
f(1,y)=y1ey=1+12y+N=2bN1yN
[4]{bN}の満たす漸化式を求める。
y=(1+12y+N=2bN1yN)(1ey)=(M=1(1)MM!yM)(1+12y+N=2bN1yN)(1+12y+N=2bN1yN=N=0cNyN)=M=1N=0(1)M1M!cNyM+N(M+N=L)=L=1M=1L(1)M1M!cLMyL係数比較により下記の漸化式を得る。
M=1N(1)M1M!cNM=(1)N1N!+(1)N22(N1)!+M=1N2(1)M1M!bNM1=M=1N2(1)NM{N(M+1)}!bM+(1)N22(N1)!+(1)N1N!=0
[5]b2n=0(nN)を確認する。
f(1,y)12y=y1+ey2(1ey)=ycoshy2sinhy2この関数は遇関数なのでy2n+1の係数は0なので[5]の結果が得られる。
[6]b1,b3,b5を試しに計算してみよう。
(i)b1(N=3)について
b114+16=b112(1213)=b1112=0b1=112
(ii)b3(N=5)について
b3+16b1148+1120=b3+112(1614+110)=b3+112(112+110)=b3+124(16+15)=b3+124130=b3+1720=0b3=1720(iii)b5(N=7)について
b5+16b3+1120b111440+15040=b514320+1144011440+15040=b514320+15040=b5+1720(16+17)=b5172042=b5130240=0b5=130240
を得る。
[6]上記をまとめると
k=mnf(k)=mnf(x)dx+f(m)+f(n)2+f(1)(n)f(1)(m)12f(3)(n)f(3)(m)720+f(5)(n)f(5)(m)30240+Rmn7

零数列の構成

零数列の構成法
[1]数列{an}nN0に対して定まる総和SN({an}nN0)のEuler-Maclaurin展開を求める。
SN({an}nN0)=S({an}nN0)+n=1D1AnNn+O(1ND)[2]二重級数{fnm}n,mN0から以下の様な数列{wn}nN0を構成する。
wn=m=0D2fn,mBm[3]次に、n=0N1wnan=δ+n=1D1AnNn+O(1ND)を満たすように漸近展開を行う。
[4]最終的にan0を満たす最小のn0を求めてwn0wn0δan0とすれば{wn}nN0D次の加速度を持つ零数列になる。

Stirlingの公式

Wallisの公式

π=limn22n(n!)4n{(2n)!}2

[0]In:=0π2sinnxdxとおく。
[1]するとInnに関して単調減少。なぜなら、x(0,1):sinx(0,1)なのでnN0:x(0,1):0<sinn+1x<sinnx
ゆえに、0<In+1<In
[2]下記の様に計算ができる
In=0π2sinnxdx=sinn1xcosx|0π2+(n1)0π2sinn2xcos2xdx=(n1)In2(n1)In
In=n1nIn2
[3][2]を用いると以下の結果を得る。
{I2n=2n12n2n32n212π2=(2n)!22n(n!)2π2I2n+1=2n2n+12n22n123=22n(n!)2(2n+1)!
[4][1]よりI2n+1<I2n<I2n1とする。すると下記の様に計算できる。
1<I2nI2n1=n{(2n)!22n(n!)2}2π<11+12n
[5]
π=limn22n(n!)4n{(2n)!}2

Stirlingの公式

N!2πNNNeN+eO(1N)

[1]Euler-Maclaurin展開より下記の様に計算できる。
logN!=n=0N1log(n+1)=0N1log(x+1)dx+logN2+1N112+=logC+NlogNN+logN2+O(1N)=log(CNN+12eN)+O(1N)
[2]上記の計算により
C=limnN!eNNN+12
[3]定数CについてはWallisの公式より下記の様に計算できる。
π=limn22n(n!)2n12{(2n)!}=12limnn!ennn+12n!ennn+12(2n)2n+12(2n)!=C2
[4]ゆえにC=2πが得られるので、最終的に以下の式が証明された。
N!=2πNNNeN+eO(1N)

中央二項係数に関するStirlingの公式

(2nn)4nπn

(2nn)=(2n)!n!n!=4πn(2n)2ne2nen2πnnnen2πnnn=4nπn

応用1

f(x)=1(x+1)2とおくと、nN0:f(n)(x)=(1)n(n+1)!(x+1)n+2とするとEuler-Maclaurin展開により
n=0N11(n+1)2=π261N+12N2+O(1N3)
と書けます。この事を用いて3次の加速度を持つ零数列を求めてください。

[1]まず以下の様に二重数列を置いてみます。
fnm=256nAm(n+1)m(n+1)4(2(n+1)n+1)4
[2]
wn=256n+1A2(n+1)2+A1(n+1)+A0(n+1)4(2(n+1)n+1)4
[3]Stirlingの公式を用いると十分大きなNに対して下記の様に書ける。
n=0N1wn(n+1)2=n=0N1256n+1A2(n+1)2+A1(n+1)+A0(n+1)6(2(n+1)n+1)4n=0N1256n+1A2(n+1)2+A1(n+1)+A0(n+1)6(4n+1π(n+1))4=π2n=0N1A2(n+1)2+A1(n+1)+A0(n+1)4=A2π2n=0N11(n+1)2+A1π2n=0N11(n+1)3+A0π2n=0N11(n+1)4
[4]そこで、以下の様に記号を定めるfk(x)=1(x+1)4k
f(n)(x)=(1)n(4k)n(x+1)4k+n
ゆえに
n=0N11(n+1)4k=ζ(4k)1(3k)N3k+12N4k4k12N5k+O(1N7k)(k=0,1,2)
[5][3][4]より次の様に計算できる。
n=0N1wn(n+1)2=A2π2ζ(2)+A1π2ζ(3)+A0π2ζ(4)A2π2N+(A2π22A1π22)1N2+O(1N3)
[6][5]より
{A2=1π2A1=0A0=ζ(2)π2ζ(4)
ゆえに4次の加速度を持つ数列は下記の様に書ける。
n=0256n+1π2(n+1)215π4(n+1)6(2(n+1)n+1)4=δ=0.243004
👆n=0N1256n+1π2(n+1)215π4(n+1)6(2(n+1)n+1)4+1/N12N2=δ+O(1N3)により計算しています。
この問題を完全に解くには上記級数の真の値を求めないといけません。
{w0=16π2240π4δwn=256n+1π2(n+1)215π4(n+1)4(2(n+1)n+1)4

δR s.t. 116π2240π4+δ+n=2N{1256nπ2n215π4n4(2nn)4}1n2=ζ(2)+O(1N3)
因みにδは次の様に書ける。
δ:=n=1256nπ2n215π4n6(2nn)4

超幾何関数に関する公式

中央二項係数の逆数のべきを含む和による級数ってかっこいいですよね?
この様な級数をガンマ関数や対数関数などを使ってビシッと答えれてたらキターンって感じです!
ではその様な計算が出来るよう公式を作成しましょう。

下記の様なp,qN0(pq)によって定まる級数S(p,q):=n=11np(2nn)qは超幾何関数により次の様に書ける。
S(p,q)=12qq+1Fq({1}p+1,{2}qp;{32}q;14q)

[1]tn:=1np(2nn)qとおく。
[2]tnは超幾何項。
tn+1tn=np(2nn)q(n+1)p(2n+2n+1)q=np(n+1)2q(n+1)p(2n+2)q(2n+1)q=np(n+1)qp2q(2n+1)q=np(n+1)qp22q+1(n+12)q
[3]ゆえに
tn=(n1)pnqp22q+1(n12)qtn1=(1)n1p(2)n1qp22q(n1)(32)n1qt1=(1)n1p(2)n1qp2q(2n1)(32)n1q
[4]以上の事を代入すると
S(p,q)=n=1tn=n=0(1)np(2)nqp2q(2n+1)(32)nq=12qn=0(1)np+1(2)nqp(32)nq1n!4qn=12qq+1Fq({1}p+1,{2}qp;{32}q;14q)

nN0:aC:{(a)2n=22n(a2)n(a+12)n(a)2n+1=22na(a+12)n(a2+1)n

実際に計算すればよい。
(a)2n=a(a+1)(a+2)(a+2n2)(a+2n1)=a(a+2)(a+4)(a+2n2)(a+1)(a+3)(a+2n1)=22na2(a2+1)(a2+2)(a2+n1)(a2+12)(a2+12+1)(a2+12+n1)=22n(a2)n(a+12)n
(a)2n+1=22na(a+12)n(a2+1)n

nN0:aC:{(a)3n=33n(a3)n(a+13)n(a+23)n(a)3n+1=3na(a+13)n(a+23)n(a3+1)n(a)3n+2=3na(a+1)(a+23)n(a3+1)n(a+13+1)n

これも実際に計算すればいい。
(a)3n=a(a+3)(a+6)(a+3n3)(a+1)(a+4)(a+3n2)(a+2)(a+5)(a+3n1)=33na3(a3+1)(a3+2)(a3+n1)a+13(a+13+1)(a+1n+n1)a+23(a+23+1)(a+23+n1)=33n(a3)n(a+13)n(a+23)n
(a)3n+1=3na(a+13)n(a+23)n(a3+1)n
(a)3n+2=3na(a+1)(a+23)n(a3+1)n(a+13+1)n

{pFq(a1,a2,...,ap;b1,b2,...,bq;z)+pFq(a1,a2,...,ap;b1,b2,...,bq;z)=22pF2q+1(a12,a1+12,a22,a2+12,...,ap2,ap+12;b12,b1+12,b22,b2+12,...,bq2,bq+12,12;(2pq1z)2)pFq(a1,a2,...,ap;b1,b2,...,bq;z)pFq(a1,a2,...,ap;b1,b2,...,bq;z)=2a1a2apb1b2bqz2pF2q+1(a1+12,a1+22,a2+12,a2+22,...,ap+12,ap+22;b1+12,b1+22,b2+12,b2+22,...,bq+12,bq+22,32;(2pq1z)2)

[1]
pFq+1(a1,a2,...,ap;b1,b2,...,bq;z)+pFq(a1,a2,...,ap;b1,b2,...,bq;z)=2n=0(a1)2n(a2)2n(ap)2n(b1)2n(b2)2n(bq)2nz2n(2n)!=2n=022(pq)(a12)n(a1+12)n(a22)n(a2+12)n(ap2)n(ap+12)n(b12)n(b1+12)n(b22)n(ab+12)n(bq2)n(bq+12)nz2n22n(12)nn!=22pF2q+1(a12,a1+12,a22,a2+12,...,ap2,ap+12;b12,b1+12,b22,b2+12,...,bq2,bq+12,12;(2pq1z)2)
[2]
pFq(a1,a2,...,ap;b1,b2,...,bq;z)pFq(a1,a2,...,ap;b1,b2,...,bq;z)=2n=0(a1)2n+1(a2)2n+1(ap)2n+1(b1)2n+1(b2)2n+1(bq)2n+1z2n+1(2n+1)!=2zn=0(a1)2n+1(a2)2n+1(ap)2n+1(b1)2n+1(b2)2n+1(bq)2n+1z2n22n(32)nn!=2zn=022(pq)na1a2ap(a1+12)n(a12+1)n(a2+12)n(a22+1)n(ap+12)n(ap2+1)nb1b2bq(b1+12)n(b12+1)n(b2+12)n(b22+1)n(bq+12)n(bq2+1)n(32)n1n!(z2)2n=2a1a2apb1b2bqzn=0(a1+12)n(a12+1)n(a2+12)n(a22+1)n(ap+12)n(ap2+1)n(b1+12)n(b12+1)n(b2+12)n(b22+1)n(bq+12)n(bq2+1)n(32)n1n!(2pq1z)2n=2a1a2apb1b2bqz2pF2q+1(a1+12,a1+22,a2+12,a2+22,...,ap+12,ap+22;b1+12,b1+22,b2+12,b2+22,...,bq+12,bq+22,32;(2pq1z)2)

楕一種円積分

k[0,1]:K(k):=0π2dθ1k2sin2θ

K(k)=π22F1(12,12;1;k2)

[1]f(x)=11xとおくとf(n)(x)=(12)n(1x)2n+12なので
f(x)=n=0(12)nn!xn
[2]ゆえにこれを用いると
K(k)=0π2n=0(12)nn!k2nsin2nθdθ=n=0(12)nn!k2n0π2sin2nθdθ=n=0(12)nn!k2n(12)nn!π2=π22F1(12,12;1;k2)

{K(k)+K(ik)=π4F3(14,14,34,34;12,1,12;k4)K(k)K(ik)=π4k24F3(34,34,54,54;1,3232;k4)

K(12)=Γ(14)24π

K(12)=0π2dθ112sin2θ=20π2dθ2sin2θ(sinθ=x)=201dx(1x2)(2x2)(x2=y)=1201dyy(1y)(2y)(1y=z)=1201dzz(1z)(1+z)=1201dzz(1z2)=12201dww34(1w)12=122B(14,12)=Γ(14)Γ(12)22Γ(34)=Γ(14)2π22Γ(14)Γ(114)=Γ(14)24π

K(i12)=π62F1(12,12;1;13)

K(i12)=0π2dθ1+12sin2θ=20π2dθ2+sin2θ(sinθ=x)=201dx(1x2)(2+x2)(x2=y)=1201dyy(1y)(2+y)(y=1z)=1201dzz(1z)(3z)=1601z121(1z)121(113z)12dz=16Γ(12)Γ(12)Γ(1)2F1(12,12;1;13)=π62F1(12,12;1;13)

{4F3(14,14,34,34;12,1,12;14)=Γ(14)24ππ+162F1(12,12;1;13)4F3(34,34,54,54;1,32,32;14)=2Γ(14)2ππ4232F1(12,12;1;13)

log(1x)=x2F1(1,1;2;x)

[1]
dndxnlog(1x)=(n1)!(1x)n
[2]
log(1x)=n=11nxn=n=1(n1)!(n1)!n!xn(n1)!=xn=0n!n!(n+1)!xnn!=x2F1(1,1;2;x)

ψ(α,β,γ,z)=01xα1(1x)β1(1zx)γdxとすると以下の様に書ける。
ψ(α,β,γ,z)=Γ(α)Γ(β)Γ(α+β)2F1(α,γ;α+β;z)

ψ(α,β,γ,z)=01xα1(1x)β1(1zx)γdx=n=0(γ)nn!zn01xα+n1(1x)β1dx=n=0(γ)nn!znΓ(α+n)Γ(β)Γ(α+β+n)=Γ(α)Γ(β)Γ(α+β)n=0(α)n(γ)n(α+β)nznn!=Γ(α)Γ(β)Γ(α+β)2F1(α,γ;α+β;z)

nN0:NN:(1N)n(2N)n(N1N)n=(nN)!n!NnN

(1N)n(2N)n(N1N)n=1N2NN1NN+1NN+2N2N1N(n1)N+1N(n1)N+2NnN1N=(nN)!n!NnN

NN:Q,RN0(0<R<N) s.t. (QN+RN)n(QN+R+1N)n(QN+N1N)n((Q+1)N+1N)n((Q+1)N+R1N)n=Q!{(Q+n)N+R1}!(Q+n)!(QN+R1)!NnN

[1]R=1とした場合。
(QN+1N)n(QN+2N)n(QN+N1N)n=QN+1NQN+2NQN+N1NQN+N+1NQN+N+2NQN+N+N1N(Q+n1)N+1N(Q+n1)N+2N(n+Q)N1N=Q!{(Q+n)N}!(Q+n)!(QN)!NnN
[2]R1とした場合。
(QN+RN)n(QN+R+1N)n(QN+N1N)n((Q+1)N+1N)n((Q+1)N+R1N)n=QN+RNQN+R+1NQN+N1N(Q+1)N+1N(Q+1)N+2N(Q+1)N+R1N(Q+1)N+RN(Q+1)N+R+1N(Q+1)N+N1N(Q+2)N+1N(Q+2)N+R1N(Q+n1)N+RN(Q+n1)N+R+1N(Q+n1)N+N1N(Q+n)N+1N(Q+n)N+R1N=Q!{(Q+n)N+R1}!(Q+n)!(QN+R1)!NnN

  1. (12)n=(2n)!n!22n
  2. (32)n={2(n+1)}!(n+1)!n22n
  3. (14)n(24)n(34)n=(4n)!n!44n
  4. (24)n(34)n(54)n=(4n+1)!n!44n
上記の式を用いると
(14)n(14)n(34)n(34)n(12)n(12)n={(4n)!}2(n!)248n(n!)428n{(2n)!}4={(4n)!}2(n!)244n{(2n)!}4={(4n)!}2(n!)244n{(2n)!}4={(4n)!}2(n!)2{(2n)!}41256n
ゆえに、
4F3(14,14,34,34;12,1,12;14)=n=0{(4n)!}2{(2n)!}411024n=Γ(14)24ππ+16n=0{(2n)!}224n(n!)313n=Γ(14)24ππ+16n=0{(2n)!}2(n!)3148n
同様に
(34)n(34)n(54)n(54)n(32)n(32)n={(4n+1)!}2(n!)248n(n!)224n{(2n)!}2{(n+1)!}2n24n[{2(n+1)}!]2=(4n+1)2{(4n)!}2(n!)2n44n+1(2n+1)2{(2n)!}4
4F3(34,34,54,54;1,32,32;14)=14n=0(4n+1)2{(4n)!}2n(2n+1)2{(2n)!}411024n=2Γ(14)2ππ423n=0{(2n)!}2(n!)3148n

応用2

Watsonの補題

積分可能な関数f:R+Cが次の性質を満たすとする。
(1)a>0 s.t. |f(t)|=O(eat)
(2)f(t)m=0cmtamただし、{am}mN0は単調増加な数列
この時、次の積分:In=0entf(t)dtを定めると、次の様に書ける。
Inm=0cmΓ(am+1)nam+1

(1)の条件により、|a|<nの様にとる。すると
CR+ s.t. |f(t)|Ceat<Cent
ゆえに
|0f(x)enxdx|0|f(x)|enxdx<C0eaxenxdx=C0e(na)xdx=Cna<+
絶対収束する。
次にその様な自然数nを取ります。
In=0f(x)enxdx0m=0cmxamenxdx(nx=y)=m=0cmnam+10xam+11exdx=m=0cmΓ(am+1)nam+1

Abelの総和公式1

数列{an}nN,{bn}nNに対して次の様な式が成り立つ。
{An:=m=1namn=1Nanbn=ANbN+n=1N1An(bnbn+1)

下の図を参照すればよい。

Abelの総和法 Abelの総和法

Abelの総和公式2

数列{an}nNと微分可能な関数fC1(0,x)に対して以下の式が成り立つ。
{A(x):=0nxan0nxanf(n)=A(x)f(x)0xA(y)f(y)dy

[1]右辺の積分について少し考察すれば分かる。
0xA(y)f(y)dy=k=1xk1kA(y)f(y)dy=k=1xlimϵ0+0k1kϵA(y)f(y)dy=k=1xA(k1)limϵ0+0k1kϵf(y)dy=k=1xA(k1)limϵ0+0{f(kϵ)f(k1)}=k=1xA(k1){f(k)f(k1)}
[2]f(k)=bkとおくと、Abelの総和と一致する事が分かるので証明終了。

積分可能な関数f:R+Cが次の性質を満たすとする。
a>0 s.t. |f(t)|=O(eat)
すると次の様な積分変換を定義することが出来る。
a<Re(s):L[f](s)=0f(x)esxdx
すると次の積分変換により元に戻る。
L1[F](x)=12πicic+iF(s)esxds

逆Laplace変換 逆Laplace変換

解析的表示f(z)=n=0anznが存在する場合のみ示す。
[1]仮定より、Laplace変換は常に定義出来なおかつ、それは絶対収束しますので次の様な計算が出来る。
L[f](s)=0f(z)eszdz=0n=0anzneszdz=n=0an0zneszdz=n=0ann!sn+1
[2]この得られた級数について以下の計算を行います。
L1[f](z)=12πicic+in=0ann!sn+1eszds+Γn=0ann!sn+1eszds=n=012πicic+iann!sn+1m=0smzmm!ds=m=0n=0ann!zmm!12πicic+i1snm+1ds=m=0n=0ann!zmm!δn,m=n=0anzn=f(z)

Perronの公式

数列{an}nN0により定まる二つの級数を考える。
{f(s):=n=1annsA(x):=1nxan
すると以下の式が成り立つ。
A(logx)=12πicic+if(s)xssds

[1]Abelの総和法により、任意の微分可能な関数φC1(1,x)に対して以下の式が成り立つ。
{A(x)=1nxan1nxanφ(n)=A(x)φ(x)1xA(t)φ(t)dt
[2]φ(x):=xsとおく。すると下記の様に書き直せる。
1nxanns=A(x)xs+s1xA(t)ts+1dt
[3]xとする。すると以下の式が成り立つ。
f(s)=s1A(t)ts+1dt(t=eτ)=s0A(τ)esτdτ
[4]ゆえに
f(s)s=0A(τ)esτdτL[A(τ)](s)
[5][4]の右辺はLaplace変換なのでLaplace変換の逆変換を行う事で
A(x)=12πicic+if(s)sesxds
最後にx=logyの様に置き換えて
A(logy)=12πicic+if(s)sysds

応用1では適当に数列を選んで三次の加速度を持つ零数列を構成出来ましたが...
しかし、n=0wnanの値が不明であり実用に足らなかったです。
ですが、今なら超幾何関数に関する公式を幾つか考え、特に楕円積分の特殊値を求めました。これらの値は既に分かっているものですから、この様な数列の一次結合を用いて級数の加速を行う事が出来るかもしれません。
僕はその様な期待を胸に計算を進めてきたのですがここで残念なお知らせ。
少なくとも、今回求めた楕円関数の特殊値、log2,2の一次結合により構成したwnでは3次の加速度を持つ零数列を構成出来ないようです。
ただ、構成できない事を証明するのも案外有益なのではないかと思い一つ例をの載せます。
良かったらこの例にチャレンジしてみてください。

f(x)=1(x+1)2とおくと、nN0:f(n)(x)=(1)n(n+1)!(x+1)n+2とするとEuler-Maclaurin展開により
n=0N11(n+1)2=π261N+12N2+O(1N3)
と書けます。この事を用いて3次の加速度を持つ零数列{wn}nN0を作りたいとします。
しかし、下記の数列のC係数の一次結合をどのようにとっても3次の加速度を持つ零数列を構成できない事を証明してください。
{fn0={{2(n+1)}!}2(n+1)2{(n+1)!}332n+1fn1=(1)n(n+1)fn2=(2(n+1)n+1)(n+1)222(n+1)(n+1)!12n+1

[1]まず以下の様に適当な係数A,B,Cを用いて下記の様な数列を作成する。
wn=A{{2(n+1)}!}2(n+1)2{(n+1)!}332n+1+B(1)n(n+1)+C(2(n+1)n+1)(n+1)222(n+1)(n+1)!12n+1
[2]以下の式が成り立つ様A,B,Cを定める。
{n=0N1wn(n+1)2=δ1N+12N2+O(1N3)n=0wn(n+1)2=AΓ(14)22ππ+Blog2+C2
漸近表示を求める。
[3]
{(2n)!}2(n!)332n=2π(2n)2n+1e4ne3n2π2πn3n+3232n=1e2π22nnn+1232n=1e2π1nn+128n
(2nn)23nn!=4n23nπnen2πnn+12=1π2ennn+12n
n=1Nwn1n2A1e2πn=1N1nn+128n+Bn=1N(1)n1n+C1π2n=1Nennn+12n
[4]そこで以下の様な記号を定める。
{f(p,q;x):=1xx+pqxψ(x):=x(logx+logq)+plogxS(p,q;f)(N):=n=1Nf(p,q;n)F(p,q;s)=n=1f(p,q;n)ns
[5]Perronの公式を用いる。
S(p,q;f)(logN)=12πicic+iF(p,q;s)Nssds
s=0の留数について
F(p,q;s)Nss=1sF(p,q;s)eslogN
{F(p,q;s)=F(p,q;0)+sF(p,q;0)+O(s2)eslogN=1+slogN+O(s2)
F(p,q;s)Nss=1s{F(p,q;0)+sF(p,q;0)+O(s2)}{1+slogN+O(s2)}=F(p,q;0)s+F(p,q;0)logN+F(p,q;0)+O(s)
以上の計算により
Ress=0F(p,q;s)Nss=F(0)=n=1f(p,q;n)
すなわち
RR+:S(p,q,f)(logN)=n=11nn+pqn+O(NR)
[6]つまり、与えられた数列fn0,fn1,fn2C係数の一次結合をどのようにとっても3次の加速度を持つ零数列を構成できない事が証明できました。

最後に

今回の記事は零数列を足す事で、級数の収束速度を上げる方法について書きました。
途中超幾何関数の公式について書いたり脱線もしましたが、そういった部分も楽しんで頂けていたら嬉しいです。
零数列に関しても思いつき次第どんどん書いていきますので引き続きよろしくお願いいたします。

投稿日:28日前
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. あいさつ
  2. 零集合
  3. Euler-Maclaurin展開
  4. 零数列の構成
  5. Stirlingの公式
  6. 応用1
  7. 超幾何関数に関する公式
  8. 応用2
  9. 最後に