例题
求limn→∞(2πn−∑k=0∞(n−k)!nkn!)
利用Gamma函数的性质
k!=∫0∞xke−xdx
利用二项式定理
∑k=0∞(n−k)!nkn!=∑k=0nnkCnk∫0∞xke−xdx=∫0∞(1+nx)ne−xdx
这里因为求和是有限的,所以可以直接换序.
作换元x→nx,得到
n∫0∞[(1+x)e−x]ndx
为了研究积分的渐近行为,把底数放到指数上
[(1+x)e−x]n=enln[(1+x)e−x]
Laplace方法专门用于处理形如∫e−nf(y)dy的积分,故引入y:
−y=ln[(1+x)e−x]=ln(1+x)−x⇒y=x−ln(1+x)
计算微分
dy=1+xxdx
从而原积分化为
n∫0∞e−ny(1+x1)dy
因为e−ny衰减得很快,积分的值由y→0+影响很大,故考虑1+x1的行为.
对于
y=x−ln(1+x)
要反解出x(y),可以利用Lagrange反演定理.
Lagrange反演
设
w=f(z)=w0+ak(z−z0)k+ak+1(z−z0)k+1+⋯,ak=0,k∈N
则
z=z0+∑n=1∞n!1[dζn−1dn−1((f(ζ)−w0)k1ζ−z0)]ζ=z0(w−w0)kn
这里给出的并不是严格叙述的版本.
在定理中取k=2,有
x=ζ→0+limζ−ln(ζ+1)ζ⋅y+21ζ→0+lim(ζ−ln(ζ+1)ζ2)′⋅y=2y+32y+o(x)
从而
1+x1=2y2+32+o(1),y→0+
进而Laplace方法可得
n∫0∞e−ny(1+x1)dy=n∫0∞e−ny(2y2+32)dy+o(1)=22πn+32+o(1),n→∞
所以
limn→∞(2πn−∑k=0∞(n−k)!nkn!)=−32