23._生成函数的唯一性和收敛性

生成函数的唯一性和收敛性

根据序列 {an}n=0\left\{a_n\right\}_{n=0}^{\infty} 的具体形式,生成函数 Ga(s)G_a(s) 可能对所有的 ss 均存在,也可能只对某些 ss 存在,还可能仅当 s=0s=0 时存在.(由于 Gs(0)=a0G_s(0)=a_0 ,所以这没有太多实际意义!)

考虑下面的例子. (1)最简单的情况是 a0=1a_0=1 且其他项都满足 an=0a_n=0 ,此时会得到 Ga(s)=1G_a(s)=1 .更一般的情况是,除了有限多个 nn 之外,其余 ana_n 全为 0 ,那么 Ga(s)G_a(s) 就是一个多项式。 (2)如果对所有的 nn 均有 an=1a_n=1 ,那么由几何级数公式可得,Ga(s)=n=0sn=G_a(s)=\sum_{n=0}^{\infty} s^n= 11s\frac{1}{1-s} .当然,为了使用几何级数公式,这里必须满足 s<1|s|<1 .对于较大的 ss ,上述级数不收玫。 (3)如果 an=1/na_n=1 / n !,那么 Ga(s)=n=0sn/n!G_a(s)=\sum_{n=0}^{\infty} s^n / n! .这是 ese ^s 的定义,所以 Ga(s)G_a(s) 对所有的 ss 均存在. (4)如果 an=2na_n=2^n ,那么 Ga(s)=n=02nsn=n=0(2s)nG_a(s)=\sum_{n=0}^{\infty} 2^n s^n=\sum_{n=0}^{\infty}(2 s)^n .这是比值为 2s2 s 的几何级数.当 2s<1|2 s|<1 时级数收玫,当 2s>1|2 s|>1 时级数发散.因此,如果 s<1/2|s|<1 / 2 ,那么 Ga(s)=(12s)1G_a(s)=(1-2 s)^{-1} . (5)如果 an=na_n=n !,那么不难看出,对于任意的 s>0,Ga(s)|s|>0, G_a(s) 均发散。最容易看出该级数发散的方法是,对于任意一个固定的 s0s \neq 0 ,只要 nn 足够大,就有 n!sn>1n!|s|^n>1 .因为级数中的项不趋向于 0 ,所以级数不能收玫.利用斯特林公式(参见第 18 章),我们可以估算 nn 必须取多大,才能保证 n!sn>1n!|s|^n>1 .由斯特林公式可知,n!(n/e)n2πnn!\sim(n / e )^n \sqrt{2 \pi n} ,从而有 n!sn>(ns/e)n.n!snn!|s|^n>(n|s| / e )^n . n!|s|^n 显然不会趋向于 0 ,因为当 n>e/sn> e /|s| 时,我们有 n!sn>1n!|s|^n>1

如果有序列 {an}n=0\left\{a_n\right\}_{n=0}^{\infty} ,那么显然可以得到它的生成函数(写出 Ga(s)G_a(s) 的解析表达式可能并不容易,但它确实有一个公式).反之亦然:如果有一个生成函数 Ga(s)G_a(s)(存在一个 δ\delta ,使得当 s<δ|s|<\delta 时该函数收玫),那么我们就可以得到原来的序列.当 Ga(s)G_a(s)存在任意阶微分时,我们能轻松地得到这个序列,因为此时有 an=1n!dnGa(s)dsna_n=\frac{1}{n!} \frac{ d ^n G_a(s)}{ d s^n} .这个结果极其重要,我们以后会经常用到它,所以现在有必要把它作为定理单独列出来.

定理 19.3.1(序列生成函数的唯一性)设 {an}n=0\left\{a_n\right\}_{n=0}^{\infty}{bn}n=0\left\{b_n\right\}_{n=0}^{\infty} 是生成函数分别为 Ga(s)G_a(s)Gb(s)G_b(s) 的两个数列。当 s<δ|s|<\delta 时,Ga(s)G_a(s)Gb(s)G_b(s) 均收敛。那么,这两个序列相等(即对于所有的 ii ,均有 ai=bia_i=b_i ),当且仅当对于所有的 s<δ|s|<\delta 均有 Ga(s)=Gb(s)G_a(s)=G_b(s) .通过对生成函数求微分,我们可以重新得到序列:an=1n!dnGa(s)dsna_n=\frac{1}{n!} \frac{ d ^n G_a(s)}{ d s^n}

证明:显然,如果 ai=bia_i=b_i ,那么 Ga(s)=Gb(s)G_a(s)=G_b(s) .对于另一个方向,如果可以对生成函数求任意阶微分,那么 ai=1i!diGa(s)dsia_i=\frac{1}{i!} \frac{ d ^i G_a(s)}{ d s^i}bi=1i!diGb(s)dsib_i=\frac{1}{i!} \frac{ d ^i G_b(s)}{ d s^i} .因为 Ga(s)=Gb(s)G_a(s)=G_b(s) ,所以这两个函数的导数也是相等的,从而有 ai=bia_i=b_i

注 19.3.2 除以 nn !看着有些不舒服,稍后我们会看到另一个不包含这个因子的生成函数。如果不想求微分,仍然可以利用生成函数来确定系数。显然,令 s=0s=0 ,我们就得到了 a0a_0 。接下来,观察 (Ga(s)a0)/s\left(G_a(s)-a_0\right) / s 并让表达式中的 ss 等于 0 ,这样就能求出 a1a_1 .按照这种方式继续进行下去,我们就可以求出任意的 ama_m .当然,注意这与微分方法是多么相似!

最后是一个简短的提醒:虽然我们写出了生成函数,但这并不意味着它是有意义的!遗憾的是,不管 ss 取什么值(当然,除了 s=0s=0 之外,此时平凡地收敛),最终的和可能都不收敛。幸运的是,与概率相关的生成函数通常是(但不总是)收玫的,至少对某些 ss 收敛,稍后我们将详细讨论这个问题。判别级数是收敛还是发散的方法有很多,B. 3 节总结了四种更流行且更强大的判别法(比值判别法,根值判别法,比较判别法以及积分判别法).

利用生成函数,最根本的是简化运算

离散型随机变量

定义19.4.1(序列的卷积)已知两个序列 {am}m=0\left\{a_m\right\}_{m=0}^{\infty}{bn}n=0\left\{b_n\right\}_{n=0}^{\infty} .它们的卷积被定义为下面这个新序列 {ck}k=0\left\{c_k\right\}_{k=0}^{\infty}

ck=a0bk+a1bk1++ak1b1+akb0=l=0kalbkl.c_k=a_0 b_k+a_1 b_{k-1}+\cdots+a_{k-1} b_1+a_k b_0=\sum_{l=0}^k a_l b_{k-l} .

通常情况下,我们把它记作 c=abc=a * b . 这个定义来自于多项式乘法.如果 f(x)=m=0amxmf(x)=\sum_{m=0}^{\infty} a_m x^mg(x)=n=0bnxng(x)=\sum_{n=0}^{\infty} b_n x^n ,并假设级数均收敛,那么

h(x)=f(x)g(x)=k=0ckxkh(x)=f(x) g(x)=\sum_{k=0}^{\infty} c_k x^k

其中 c=abc=a * b .例如,若 f(x)=2+3x4x2f(x)=2+3 x-4 x^2g(x)=5x+x3g(x)=5-x+x^3 ,那么 f(x)g(x)=f(x) g(x)= 10+13x23x2+6x3+3x44x510+13 x-23 x^2+6 x^3+3 x^4-4 x^5 .由定义可知,c2c_2 应该等于:

a0b2+a1b1+a2b0=20+3(1)+(4)5=23a_0 b_2+a_1 b_1+a_2 b_0=2 \cdot 0+3 \cdot(-1)+(-4) \cdot 5=-23

这恰好是 f(x)f(x)g(x)g(x) 相乘后得到的结果. 引理 19.4.2 设 Ga(s)G_a(s){am}m=0\left\{a_m\right\}_{m=0}^{\infty} 的生成函数,Gb(s)G_b(s){bn}n=0\left\{b_n\right\}_{n=0}^{\infty} 的生成函数,那么 c=abc=a * b 的生成函数是 Gc(s)=Ga(s)Gb(s)G_c(s)=G_a(s) G_b(s)

连续型随机变量

定义 19.5.1(概率生成函数)设 XX 是一个连续型随机变量,其概率密度函数为 ff ,那么

GX(s)=sxf(x)dxG_X(s)=\int_{-\infty}^{\infty} s^x f(x) d x

XX 的概率生成函数.

定义 19.5.2(函数卷积)两个函数 f1f_1f2f_2 的卷积,即 f1f2f_1 * f_2 被定义为

(f1f2)(x)=f1(t)f2(xt)dt.\left(f_1 * f_2\right)(x)=\int_{-\infty}^{\infty} f_1(t) f_2(x-t) d t .

如果 fif_i 都是概率密度函数,那么上述积分收敛.