随机试验与随机事件

随机试验

在自然界和人类社会中出现的各种现象大致可以分为两类:必然现象和随机现象。在一定条件下必然出现的现象叫做必然现象;在相同的条件下,可能出现不同的结果,而在试验或观测之前不能预知确切的结果的现象叫做随机现象

为了研究和揭示随机现象的统计规律性,我们需要对随机现象进行大量重复的观察、测量或试验。具有以下特点的试验称为随机试验

  1. 可重复性:可以在相同条件下重复进行;
  2. 可观测性:每次试验的所有可能结果都是明确的,可以观测的,并且试验结果有两个或两个以上;
  3. 随机性:每次试验结果不确定。

随机试验简称试验,通常用字母 EE 表示。随机试验 EE 的基本结果称为样本点,用 ω\omega 表示。称随机试验 EE 的所有基本结果的集合为样本空间,用 Ω=ω\varOmega = |\omega| 表示。

随机事件及其运算

我们把随机试验 EE 的样本空间 Ω=ω\varOmega = |\omega| 的子集称为随机试验 EE随机事件,简称为事件,用大写字母 A,B,CA, B, C 等表示。设 AΩA \subseteq \varOmega,如果试验结果 ωA\omega \in A,则称在这次试验中事件 AA 发生;如果 ωA\omega \notin A,则称事件 AA 不发生。

下面是一些基本概念:

  1. 随机事件:样本空间的子集,是由若干样本点构成的集合,记为 A,B,CA, B, C 等。
  2. 事件发生:当事件 AA 所包含的样本点有一个出现,就说事件 AA 发生了,否则就说事件 AA 未发生。
  3. 基本事件:由一个样本点组成的事件。
  4. 必然事件:包含样本空间 Ω\varOmega 的所有样本点,每次试验必然发生。
  5. 不可能事件:空集 \empty,不包含任何样本点,每次试验都不可能发生。

随机事件的关系与运算:

名称 表达式 含义
包含 ABA \subset B A 发生必然导致 B 发生
相等 A=BA = B\\ ABA \subset BBAB \subset A A 与 B 为同一事件
并(和) ABA \cup B A 与 B 至少有一个发生
交(积) AB=ABA \cap B = AB A 与 B 同时发生
ABA - B A 发生,B 不发生
互不相容(互斥) AB=AB = \empty A 与 B 不能同时发生
互逆(对立) A=B\overline{A} = B\\ AB=ΩA \cup B = \varOmegaAB=AB = \empty A 与 B 中有且只有一个发生

基本运算律:

运算顺序:逆交并差,括号优先。

  1. 交换律:AB=BAA \cup B = B \cup AAB=BAAB = BA
  2. 结合律:(AB)C=A(BC)(A \cup B) \cup C = A \cup (B \cup C)(AB)C=A(BC)(AB)C = A(BC)
  3. 分配律:(AB)C=(AC)(BC)(A \cup B)C = (AC) \cup (BC)A(BC)=(AB)(AC)A \cup (BC) = (A \cup B)(A \cup C)
  4. 对偶律:AB=A B\overline{A \cup B} = \overline{A}\ \overline{B}AB=AB\overline{AB} = \overline{A} \cup \overline{B}

i=1nAi=i=1nAii=1nAi=i=1nAi\overline{\bigcup_{i = 1}^n A_i} = \bigcap_{i = 1}^n \overline{A_i} \qquad \overline{\bigcap_{i = 1}^n A_i} = \bigcup_{i = 1}^n \overline{A_i}

  1. 吸收律:若 ABA \subset B,则 AB=BA \cup B = BAB=AAB = A.
  2. 差化积:AB=AB=AABA - B = A \overline{B} = A - AB
  3. 逆:A=A\overline{\overline{A}} = A
  4. AB=    AB    BAAB = \empty \iff A \subset \overline{B} \iff B \subset \overline{A}

随机事件的概率

频率与概率

定义: 设在相同的条件下进行的 nn 次试验中,事件 AA 发生了 nAn_A 次,则称 nAn_A 为事件 AA 发生的频数,称比值 nAn\dfrac{n_A}{n} 为事件 AA 发生的频率,记作 fn(A)f_n(A),即

fn(A)=nAnf_n(A) = \dfrac{n_A}{n}

频率具有如下基本性质:

  1. 非负性:对于任意事件 AA,有 fn(A)0f_n(A) \geq 0
  2. 规范性:对于必然事件 Ω\varOmega,有 fn(Ω)=1f_n(\varOmega) = 1
  3. 优先可加性: 对于两两互不相容的事件 A1,A2,,AmA_1, A_2, \cdots, A_m(即当 iji \neq j 时,有 AiAj=A_i A_j = \emptyi,j=1,2,,mi, j = 1, 2, \cdots, m),有

fn(i=1mAi)=i=1mfn(Ai)f_n \left( \bigcup_{i = 1}^m A_i \right) = \sum_{i = 1}^m f_n \left( A_i \right)

大量的试验表明,在相同的条件下重复进行 nn 次试验,当 nn 增大时,事件 AA 发生的频率 fn(A)f_n(A) 呈现出稳定性,逐渐稳定于某一常数 pp,我们用这一常数 pp 表示事件 AA 发生的可能性大小,称为事件 AA概率,记为 P(A)=pP(A) = p.

定义: 设随机试验 EE 的样本空间为 Ω\varOmega,如果对于 EE 的每一个事件 AA,有唯一的实数 P(A)P(A) 和它对应,并且这个事件的函数 P(A)P(A) 满足以下条件:

  1. 非负性:对于任一事件 AA,有 P(A)0P(A) \geq 0
  2. 规范性:对于必然事件 Ω\varOmega,有 P(Ω)=1P(\varOmega) = 1
  3. 可列可加性:对于两两互不相容的事件 A1,A2,A_1, A_2, \cdots,有

P(i=1Ai)=i=1P(Ai)P \left( \bigcup_{i = 1}^\infty A_i \right) = \sum_{i = 1}^\infty P \left( A_i \right)

则称 P(A)P(A) 为事件 AA 的概率。

根据上述定义,可以得到概率的一些性质:

性质 1: 对于不可能事件 \empty,有 P()=0P(\empty) = 0.

性质 2(有限可加性): 对于两两互不相容的事件 A1,A2,,AnA_1, A_2, \cdots, A_n,有

P(i=1nAi)=i=1nP(Ai)P \left( \bigcup_{i = 1}^n A_i \right) = \sum_{i = 1}^n P \left( A_i \right)

性质 3: 对于任一事件 AA,有 P(A)=1P(A)P(\overline{A}) = 1 - P(A).

性质 4: 如果 ABA \subseteq B,则有

P(BA)=P(B)P(A),P(A)P(B)P(B - A) = P(B) - P(A), \quad P(A) \leq P(B)

性质 5: 对于任一事件 AA,有 P(A)1P(A) \leq 1.

性质 6(减法公式): 对于任意两个事件 AABB,有

P(BA)=P(B)P(AB)P(B - A) = P(B) - P(AB)

性质 7(加法公式): 对于任意两个事件 AABB,有

P(AB)=P(A)+P(B)P(AB),P(AB)P(A)+P(B)P(A \cup B) = P(A) + P(B) - P(AB),\\ P(A \cup B) \leq P(A) + P(B)

推广:概率的一般加法公式为

P(i=1nAi)=i=1nP(Ai)1i<jnP(AiAj)+1i<j<knP(AiAjAk)++(1)n1P(A1A2An)\begin{aligned} P \left( \bigcup_{i = 1}^n A_i \right) &= \sum_{i = 1}^n P \left( A_i \right) - \sum_{1 \leq i < j \leq n} P \left( A_i A_j \right) + \sum_{1 \leq i < j < k \leq n} P \left( A_i A_j A_k \right)\\ &+ \cdots + (-1)^{n - 1} P \left( A_1 A_2 \cdots A_n \right) \end{aligned}

古典概型

如果随机试验具有以下两个特点:

  1. 试验的样本空间 Ω\varOmega 只包含有限个样本点;
  2. 在试验中每个基本事件发生的可能性相同,

则称这种试验为等可能概型古典概型

古典概型中事件 AA 的概率计算方式:

P(A)=A包含的基本事件个数Ω包含的基本事件总数P(A) = \frac{A \text{包含的基本事件个数}}{\varOmega \text{包含的基本事件总数}}

由这一公式可知,要计算古典概型中事件 AA 的概率,只需算出样本空间 Ω\varOmega 包含的基本事件总数及事件 AA 包含的基本事件个数。这时常常用到加法原理、乘法原理以及排列组合公式。

特点: 有限性和等可能性

排列组合相关知识:

  1. 加法原理:完成一件事有 nn 类方法,第 ii 类方法中有 mim_i 种具体的方法,则完成这件事共有 i=1nmi\displaystyle\sum_{i = 1}^n m_i 种不同的方法。
  2. 乘法原理:完成一件事有 nn 个步骤,第 ii 个步骤中有 mim_i 种具体的方法,则完成这件事共有 i=1nmi\displaystyle\prod_{i = 1}^n m_i 种不同的方法。
  3. 组合:从 nn 个不同的元素中(不放回地)取出 mm 个组成一组,不同的分法共有

Cnm=n!m!(nm)!C_n^m = \frac{n!}{m!(n - m)!}

  1. 排列:从 nn 个不同的元素中(不放回地)取出 mm 个按一定的次序排成一排,不同的排法共有

Anm=n(n1)(n2)(nm+1)=n!(nm)!A_n^m = n(n - 1)(n - 2) \cdots (n - m + 1) = \frac{n!}{(n - m)!}

  1. 全排列:Ann=n!A_n^n = n!
  2. 可重复排列:从 nn 个不同的元素中可重复地取出 mm 个排成一排,不同的排法共有 nmn^m.
  3. 分类分组nn 个元素分 mm 类,第 ii 类中有 kik_i 个相同地元素,且 k1+k2++km=nk_1 + k_2 + \cdots + k_m = n,不同的排法共有

Cnk1Cnk1k2Ckmkm=n!k1!k2!km!C_n^{k_1} C_{n - k_1}^{k_2} \cdots C_{k_m}^{k_m} = \frac{n!}{k_1!k_2! \cdots k_m!}

摸球问题: 袋中有 aa 只白球,bb 只红球,每次任取一只,取后不放回,进行 k+1 (k+1a+b)k + 1\ (k + 1 \leq a + b) 次,则最后一次取到白球的概率为

P(A)=Aa1Aa+b1kAa+bk+1=aa+bP(A) = \frac{A_a^1 A_{a + b - 1}^k}{A_{a + b}^{k + 1}} = \frac{a}{a + b}

随机取数问题:

例:从 0, 1, 2, …, 9 这十个数字中任意取出 4 个,求能排成一个四位偶数的概率。

解:从 0 至 9 中任取 4 个数并进行排列得到的数字的个数为 C104A44=5040C_{10}^4 A_4^4 = 5040.

  1. 当个位为 0 时,排列得到的四位偶数的个数为 C93A33=504C_9^3 A_3^3 = 504
  2. 当个位不为0时,四位偶数的个数为 C41C81C82A22=1792C_4^1 C_8^1 C_8^2 A_2^2 = 1792.

综上所述,所求概率为 504+17925040=4190\dfrac{504 + 1792}{5040} = \dfrac{41}{90}.

分房问题:nn 个人随机地分配到 N (nN)N\ (n \leq N) 个房间,则每个房间至多有一人的概率为

P(A)=ANnNnP(A) = \frac{A_N^n}{N^n}

某指定的 nn 个房间各有一人的概率为

P(B)=n!NnP(B) = \frac{n!}{N^n}

几何概型

在古典概型中利用等可能性计算了一类问题的概率,但古典概型要求基本事件的总数必须是有限个。人们希望把这种做法推广到无限个基本事件,而这些基本事件又有某种等可能性的情形。

如果随机试验是将一个点随机地投到某一区域 Ω\varOmegaΩ\varOmega 可以是直线上的某一区间,也可以是平面或空间内的某一区域)内,而这个随机点落在 Ω\varOmega 中任意两个度量相等的子区域内的可能性是一样的,则称这样的试验属于几何概型。对于任何有度量的子区域 AΩA \subseteq \varOmega,我们定义事件 $A = $ “随机点落在区域 AA 内” 的概率为

P(A)=A的度量Ω的度量P(A) = \frac{A \text{的度量}}{\varOmega \text{的度量}}

特点: 无限性和等可能性

蒲丰投针问题: 在平面上画有等距离的平行线,平行线间的距离为 2a (a>0)2a\ (a > 0). 向该平面任意投掷一枚长为 2l (l<a)2l\ (l < a) 的圆柱形的针,求此针与任一平行线相交的概率。

解: 针投在该平面上,以 xx 表示针的中点 MM 到最近的一条平行线的距离,以 φ\varphi 表示针与此直线的交角,则有

0xa,0φπ0 \leq x \leq a, \quad 0 \leq \varphi \leq \pi

针与最近的一条平行线相交的充分必要条件是

xlsinφx \leq l\sin{\varphi}

Ω={(φ,x)0φπ,0xa}A={(φ,x)0φπ,0xlsinφ}\begin{aligned} \varOmega &= \{ (\varphi, x) \mid 0 \leq \varphi \leq \pi, 0 \leq x \leq a \}\\ A &= \{ (\varphi, x) \mid 0 \leq \varphi \leq \pi, 0 \leq x \leq l\sin{\varphi} \} \end{aligned}

则所求概率为

p=A的面积Ω的面积=0πlsinφdφπa=2lπap = \frac{A \text{的面积}}{\varOmega \text{的面积}} = \frac{\int_0^\pi l\sin{\varphi} \mathrm{d}\varphi}{\pi a} = \frac{2l}{\pi a}

条件概率

条件概率与乘法公式

定义:AABB 是试验 EE 的两个事件,且 P(A)>0P(A) > 0,称 P(AB)/P(A)P(AB) / P(A) 为在事件 AA 已经发生的条件下,事件 BB 发生的条件概率,记为 P(BA)P(B \mid A),即

P(BA)=P(AB)P(A)P(B \mid A) = \frac{P(AB)}{P(A)}

由这个定义可知,对于任意两个事件 AABB,如果 P(A)>0P(A) > 0,则有

P(AB)=P(A)P(BA)P(AB) = P(A) P(B \mid A)

称上式为概率的乘法公式

条件概率具有如下性质:

  1. 非负性:对任意事件 BB,有 P(BA)0P(B \mid A) \geq 0
  2. 规范性:对于必然事件 Ω\varOmega,有 P(ΩA)=1P(\varOmega \mid A) = 1
  3. 可列可加性:对于两两互不相容的事件 B1,B2,B_1, B_2, \cdots, 有

P((i=1Bi)A)=i=1P(BiA)P\left( \left( \bigcup_{i = 1}^\infty B_i\right) \mid A \right) = \sum_{i = 1}^\infty P(B_i \mid A)

可由条件概率的三个基本性质推导出其他一些性质:

  1. P(A)=0P(\empty \mid A) = 0
  2. P(BA)=1P(BA)P(\overline B \mid A) = 1 - P(B \mid A)
  3. P((B1B2)A)=P(B1A)P(B1B2A)P((B_1 - B_2) \mid A) = P(B_1 \mid A) - P(B_1 B_2 \mid A)
  4. P((B1B2)A)=P(B1A)+P(B2A)P((B1B2)A)P((B_1 \cup B_2) \mid A) = P(B_1 \mid A) + P(B_2 \mid A) - P((B_1 B_2) \mid A).

乘法公式的推广:

  1. P(AB)>0P(AB) > 0 时,P(ABC)=P(A)P(BA)P(CAB)P(ABC) = P(A) P(B \mid A) P(C \mid AB)
  2. P(A1A2An1)>0P(A_1 A_2 \cdots A_{n - 1}) > 0 时,P(A1A2An)=P(A1)P(A2A1)P(AnA1A2An1)P(A_1 A_2 \cdots A_n) = P(A_1) P(A_2 \mid A_1) \cdots P(A_n \mid A_1 A_2 \cdots A_{n - 1}).

例:设 M 件产品中有 m 件是不合格品,从中任取两件,在所取产品中有一件是不合格品的条件下,求另一件也是不合格品的概率。

解:设 A 表示事件 “所取产品中至少有一个是不合格品”,B 表示事件 “所取产品都是不合格品”,则

P(A)=Cm2+Cm1CMm1CM2,P(B)=Cm2CM2P(A) = \frac{C_m^2 + C_m^1 C_{M - m}^1}{C_M^2}, \quad P(B) = \frac{C_m^2}{C_M^2}

又因为 P(AB)=P(B)P(AB) = P(B),所以所求概率为

P(BA)=P(AB)P(A)=P(B)P(A)=m12Mm1P(B \mid A) = \frac{P(AB)}{P(A)} = \frac{P(B)}{P(A)} = \frac{m - 1}{2M - m - 1}

结论:A,BA, B 为随机事件,若 0<P(A)<1,0<P(B)<10 < P(A) < 1, 0 < P(B) < 1,则 P(AB)>P(AB)P(A \mid B) > P(A \mid \overline{B}) 的充要条件是 P(BA)>P(BA)P(B \mid A) > P(B \mid \overline{A}).

证明:

P(AB)>P(AB)P(A \mid B) > P(A \mid \overline{B}) 等价于 P(AB)P(B)>P(AB)P(B)=P(A)P(AB)1P(B)\dfrac{P(AB)}{P(B)} > \dfrac{P(A \overline{B})}{P(\overline{B})} = \dfrac{P(A) - P(AB)}{1 - P(B)},即 P(AB)>P(A)P(B)P(AB) > P(A) P(B).

P(BA)>P(BA)P(B \mid A) > P(B \mid \overline{A}) 等价于 P(AB)P(A)>P(B)P(AB)1P(A)\dfrac{P(AB)}{P(A)} > \dfrac{P(B) - P(AB)}{1 - P(A)},即 P(AB)>P(A)P(B)P(AB) > P(A) P(B).

全概率公式

A1,A2,,AnA_1, A_2, \cdots, A_nΩ\varOmega 中的事件,满足

  1. 事件 A1,A2,,AnA_1, A_2, \cdots, A_n 互不相容(即 AiAj=, ij, i,j=1,2,A_i A_j = \empty, \ i \neq j, \ i, j = 1, 2, \cdots);
  2. iAi=Ω\bigcup_i A_i = \varOmega.

则称 A1,A2,,AnA_1, A_2, \cdots, A_nΩ\varOmega 的一个 (有限)分割完全事件组

如果 P(Ai)>0 (i=1,2,,n)P(A_i) > 0 \ (i = 1, 2, \cdots, n),则

B=i=1nAiB,(AiB)(AjB)=B = \bigcup_{i = 1}^n A_i B, \quad (A_i B)(A_j B) = \empty

根据乘法公式得

P(B)=i=1nP(Ai)P(BAi)P(B) = \sum_{i = 1}^n P(A_i) P(B \mid A_i)

这个公式称为全概率公式,它是概率论的一个基本公式。

例:某车间有四个班组生产同一种产品,其产量分别占总产量的 15%、20%、30%、35%,次品率分别为 0.05、0.04、0.03、0.02,现从全部产品中任取一件,问恰好取到次品的概率是多少?

解:设 AiA_i 表示事件 “取到 ii 组的产品”,i=1,2,3,4i = 1, 2, 3, 4BB 表示事件 “恰好取到次品”,

P(B)=i=14P(Ai)P(BAi)=0.15×0.05+0.2×0.04+0.3×0.03+0.35×0.02=0.0315\begin{aligned} P(B) &= \sum_{i = 1}^4 P(A_i) P(B \mid A_i)\\ &= 0.15 \times 0.05 + 0.2 \times 0.04 + 0.3 \times 0.03 + 0.35 \times 0.02\\ &= 0.0315 \end{aligned}

贝叶斯(Bayes)公式

Ω\varOmega 为试验 EE 的基本空间,BB 为任意事件,A1,A2,,AnA_1, A_2, \cdots, A_nΩ\varOmega 的一个划分,且 P(Ai)>0 (i=1,2,,n)P(A_i) > 0 \ (i = 1, 2, \cdots, n),如果 P(B)>0P(B) > 0,则

P(AjB)=P(Aj)P(BAj)i=1nP(Ai)P(BAi), j=1,2,,nP(A_j \mid B) = \frac{P(A_j) P(B \mid A_j)}{\sum\limits_{i = 1}^n P(A_i) P(B \mid A_i)}, \ j= 1, 2, \cdots, n

这个公式称为贝叶斯公式

其中 P(Ai)P(A_i) 称为先验概率,它是由以往的经验得到的,是事件 B 的原因;条件概率 P(AiB)P(A_i \mid B) 称为后验概率,它是事件 B 发生后,再对导致 B 发生的原因的可能性大小重新加以评估得到的。

例:学生在做一道有 4 个选项的选择题时,如果他不知道问题的正确答案,就做随机猜测。假设学生认为自己知道正确答案的概率为 0.2,试求学生确实知道正确答案的概率。

解:设事件 AA 为 “题答对了”,事件 BB 为 “知道正确答案”,则由题意得

P(AB)=1,P(AB)=0.25P(A \mid B) = 1, \quad P(A \mid \overline{B}) = 0.25

所以由贝叶斯公式得

P(BA)=P(B)P(AB)P(B)P(AB)+P(B)P(AB)=0.2×10.2×1+0.8×0.25=0.5\begin{aligned} P(B \mid A) &= \frac{P(B) P(A \mid B)}{P(B) P(A \mid B) + P(\overline{B}) P(A \mid \overline{B})}\\ &= \frac{0.2 \times 1}{0.2 \times 1 + 0.8 \times 0.25} = 0.5 \end{aligned}

所求概率即为 0.5。

事件的独立性

定义:A,BA, B 是试验 EE 的两个随机事件,如果

P(AB)=P(A)P(B)P(AB) = P(A) P(B)

则称事件 AA 和事件 BB 相互独立

容易得到如下结论:对于同一试验 EE 的两个事件 AABB,如果 P(A)>0P(A) > 0,则 AABB 相互独立的充分必要条件为 P(BA)=P(B)P(B \mid A) = P(B);如果 P(B)>0P(B) > 0,则 AABB 相互独立的充分必要条件为 P(AB)=P(A)P(A \mid B) = P(A).

结论 1: 若事件 AABB 相互独立,则 AAB\overline{B}BBA\overline{A}A\overline{A}B\overline{B} 均独立。

证明:因为

P(AB)=P(A)P(AB)=P(A)P(A)P(B)=P(A)[1P(B)]=P(A)P(B)\begin{aligned} P(A \overline{B}) &= P(A) - P(AB)\\ &= P(A) - P(A) P(B)\\ &= P(A)[1 - P(B)]\\ &= P(A) P(\overline{B}) \end{aligned}

因此事件 AA 与事件 B\overline{B} 是相互独立的。同理,可证明其它结论。

结论 2:P(A)>0,P(B)>0P(A) > 0, P(B) > 0,则事件 A,BA, B 相互独立与 A,BA, B 互不相容不能同时成立。

结论 3: 必然事件 Ω\varOmega 和不可能事件 Φ\varPhi 与任何事件都是独立的。

三个事件的独立性: 对于同一试验的三个事件 A,B,CA, B, C,如果满足

P(AB)=P(A)P(B)P(BC)=P(B)P(C)P(AC)=P(A)P(C)P(AB) = P(A) P(B)\\ P(BC) = P(B) P(C)\\ P(AC) = P(A) P(C)

则称三个事件 A,B,CA, B, C两两相互独立的。此时如果

P(ABC)=P(A)P(B)P(C)P(ABC) = P(A) P(B) P(C)

则称三个事件 A,B,CA, B, C相互独立的。

n 个事件的独立性:A1,A2,,AnA_1, A_2, \cdots, A_n 是同一试验的 nn 个事件,如果对于任意正整数 kk 及这 nn 个事件中的任意 k (2kn)k \ (2 \leq k \leq n) 个事件 Ai1,Ai2,,AikA_{i_1}, A_{i_2}, \cdots, A_{i_k},都有等式

P(Ai1Ai2Aik)=P(Ai1)P(Ai2)+P(Aik)P(A_{i_1} A_{i_2} \cdots A_{i_k}) = P(A_{i_1}) P(A_{i_2}) \cdots + P(A_{i_k})

成立,则称 nn 个事件 A1,A2,,AnA_1, A_2, \cdots, A_n相互独立的。

n 个事件相互独立的含义是它们中的几个事件发生与否不影响另外一些事件发生的概率。在实际问题中,往往不是用定义去判断事件的独立性,而是根据问题的实际意义去确定。

伯努利(Bernoulli)概型

n 重伯努利试验概型:

  • 试验重复做 n 次;
  • 每次试验只有两个可能的结果:AAA\overline{A},且 P(A)=p (0<p<1)P(A) = p \ (0 < p < 1)
  • 每次试验的结果与其他次试验无关,称这 n 次试验是相互独立的。

在 n 重伯努利概型中,事件 A 发生 k 次的概率为

Pn(k)=Cnkpk(1p)nk,k=0,1,2,,n.P_n(k) = C_n^k p^k (1 - p)^{n - k}, \quad k = 0, 1, 2, \cdots, n.

该式又称为二项概率公式

概率很小的随机事件在一次试验中实际上几乎不发生,这一原理称为小概率事件的实际不可能原理,又称为实际推断原理。

例: 设每次试验的成功率为 p (0<p<1)p \ (0 < p < 1),将试验独立重复地进行,求直到第 nn 次才取得第 r (1<r<n)r \ (1 < r < n) 次成功的概率。

解:总共进行 nn 次试验,其中有 rr 次成功,且第 rr 次成功,有 nrn - r 次失败,所以令 k=nrk = n - r,即得所求概率为

Cr+k1r1pr(1p)k=Cr+(nr)1r1pr(1p)nr=Cn1r1pr(1p)nrC_{r + k - 1}^{r - 1} p^r (1 - p)^k = C_{r + (n - r) - 1}^{r - 1} p^r (1 - p)^{n - r} = C_{n - 1}^{r - 1} p^r (1 - p)^{n - r}

随机变量及其分布

随机变量的分布函数

定义: 设随机试验 EE 的样本空间为 Ω={ω}\varOmega = \{ \omega \}. 如果对每一个 ωΩ\omega \in \varOmega,都有一个实数 X(ω)X(\omega) 与之对应,则称 X=X(ω) (ωΩ)X = X(\omega)\ (\omega \in \varOmega)随机变量

随机变量 XX 是样本点 ω\omega 的函数,其定义域为 Ω\varOmega,值域为 RR.

随机变量随着试验结果的不同而变化,常用大写字母 X,Y,ZX, Y, Z 等表示。

由于随机变量的取值依赖于随机试验的结果,因此在试验之前我们只能知道它的所有可能取值,而不能预先知道它究竟取哪个值。因为试验的各个结果的出现都有一定的概率,所以随机变量取相应的值也有确定的概率。

用随机变量可以表示随机事件,{X=k},{a<X<b},{X>b}\{ X = k \}, \{ a < X < b \}, \{ X > b \} 等都是随机事件。

为了研究随机变量,我们引入分布函数的概念。

定义:XX 是一个随机变量,对于任意实数 xx,令

F(x)=P{Xx},xR,F(x) = P\{ X \leq x \}, \quad x \in \R,

F(x)F(x) 为随机变量 XX分布函数

随机变量 XX 的分布函数 F(x)F(x) 是定义在 (.+)(-\infty. +\infty) 上的函数,是随机事件 {Xx}\{ X \leq x \} 发生的概率。如果 XX 是数轴上随机点的坐标,则分布函数值 F(a)F(a) 表示 XX 落在区间 (,a](-\infty, a] 上的概率,如下图所示:

由于分布函数是普通的一元函数,因此通过它我们可以利用数学分析的方法来研究随机变量。

分布函数的基本性质:

性质 1: 对于任意实数 xx,有 0F(x)10 \leq F(x) \leq 1.

性质 2: 对于任意实数 x1,x2 (x1<x2)x_1, x_2\ (x_1 < x_2),有

P{x1<Xx2}=F(x2)F(x1).P\{ x_1 < X \leq x_2 \} = F(x_2) - F(x_1).

性质 3:

F()=limx(x)=0F(+)=limx+(x)=1F(-\infty) = \lim_{x \to -\infty}(x) = 0\\ F(+\infty) = \lim_{x \to +\infty}(x) = 1

性质 4: F(x)F(x) 处处右连续,即 F(x+)=F(x)F(x^+) = F(x).

性质 5: 对于任意实数 x1,x2 (x1<x2)x_1, x_2\ (x_1 < x_2),由于

F(x2)F(x1)=P{x1<Xx2}0,F(x_2) - F(x_1) = P\{ x_1 < X \leq x_2 \} \geq 0,

因此 F(x)F(x) 是一个单调不减函数。

由于 F(a)=P{Xa}F(a) = P\{ X \leq a \},则

P{X=a}=F(a)F(a0)P{X<a}=F(a0)P{X>a}=1F(a)P{Xa}=1F(a0)\begin{aligned} P\{ X = a \} &= F(a) - F(a - 0)\\ P\{ X < a \} &= F(a - 0)\\ P\{ X > a \} &= 1 - F(a)\\ P\{ X \geq a \} &= 1 - F(a - 0)\\ \end{aligned}

离散型随机变量及其概率分布

定义 1: 如果一个随机变量 XX 所有可能取到的不相同的值是有限个或可列无限多个,并且以确定的概率取这些不同的值,则称 XX离散型随机变量

定义 2: 设离散型随机变量 XX 所有可能取的值为 xk (k=1,2,)x_k\ (k = 1, 2, \cdots)XX 取各个可能值的概率,即事件 {X=xk}\{ X = x_k \} 的概率为

P{X=xk}=pk,k=1,2,P\{ X = x_k \} = p_k, \quad k = 1, 2, \cdots

并且 pkp_k 满足以下两个条件:

  1. pk0,k=1,2,p_k \geq 0, \quad k = 1, 2, \cdots
  2. k=1pk=1\sum_{k = 1}^\infty p_k = 1

则称上式为离散型随机变量 XX概率分布分布律

概率分布也可以用如下的表格来表示:

离散随机变量的分布函数:

F(x)=P{Xx}=xkxP{X=xk}=xkxpkF(x) = P\{ X \leq x \} = \sum_{x_k \leq x} P\{ X = x_k \} = \sum_{x_k \leq x} p_k

离散型随机变量 XX 的分布函数是一条阶梯形曲线,所有可能取值点为间断点,每个间断点处的跳跃度恰为取该点的概率值。

常用的离散型随机变量

离散型均匀分布:

若随机变量 XX 的分布律为

P{X=xk}=1n,k=1,2,,n.P\{ X = x_k \} = \frac{1}{n}, \quad k = 1, 2, \cdots, n.

iji \neq j 时xixjx_i \neq x_j,则称 XX 服从离散型均匀分布

(0-1) 分布:

若随机变量 XX 只取 0 与 1 两个值,其概率分布为

P{X=0}=1p,P{X=1}=p,P\{ X = 0 \} = 1 - p, \quad P\{ X = 1 \} = p,

其中 0<p<10 < p < 1,则称 XX 服从参数为 pp(0-1) 分布两点分布

分布律也可以写成

P{X=x}=px(1p)1x,x=0,1.P\{ X = x \} = p^x (1 - p)^{1 - x}, \quad x = 0, 1.

二项分布:

在 n 重伯努利试验中,以 XX 表示事件 AA 发生的次数,所有可能取值为 0,1,2,,n0, 1, 2, \cdots, n,设 P(A)=p (0<p<1)P(A) = p \ (0 < p < 1).
若随机变量 XX 的分布律为

P{X=k}=Cnkpk(1p)nk,k=0,1,2,,n,P\{ X = k \} = C_n^k p^k (1 - p)^{n - k}, \quad k = 0, 1, 2, \cdots, n,

其中 0<p<10 < p < 1,则称 XX 服从参数为 n,pn, p二项分布,记作 XB(n,p)X \sim B(n, p).

易得

P{X=k}=Cnkpk(1p)nk>0,k=0nP{X=k}=k=0nCnkpk(1p)nk=[p+(1p)]n=1.P\{ X = k \} = C_n^k p^k (1 - p)^{n - k} > 0,\\ \sum_{k = 0}^n P\{ X = k \} = \sum_{k = 0}^n C_n^k p^k (1 - p)^{n - k} = [p + (1 - p)]^n = 1.

特别地,当 n=1n = 1 时,二项分布 B(1,p)B(1, p) 就是 (0-1) 分布。

二项分布的最可能取值

k0={(n+1)p(n+1)p1,(n+1)p是整数时[(n+1)p](n+1)p不是整数时k_0 = \begin{cases} (n + 1)p \text{与} (n + 1)p - 1, & \text{当} (n + 1)p \text{是整数时} \\ [(n + 1)p] & \text{当} (n + 1)p \text{不是整数时} \end{cases}

对于任意 jXj \in X,若 P{X=k}P{X=k}P\{ X = k \} \geq P\{ X = k \},则称 kk 为最可能出现的次数。记 pk=P{X=k}=Cnkpk(1p)nkp_k = P\{ X = k \} = C_n^k p^k (1 - p)^{n - k},考虑

δ=pkpk1=(nk+1)pk(1p)=1+(n+1)pkk(1p)1\delta = \frac{p_k}{p_{k - 1}} = \frac{(n - k + 1)p}{k(1 - p)} = 1 + \frac{(n + 1)p - k}{k(1 - p)} \geq 1

k(n+1)pk \leq (n + 1)p.

几何分布:

假设试验只有两个可能的对立结果 AAA\overline{A},且 P(A)=p,P(A)=1pP(A) = p, P(\overline{A}) = 1 - p,将试验独立重复地进行下去,直到 AA 发生为止,以 XX 表示实验进行的次数,这样得到的分布结果就是几何分布。

若随机变量 XX 的分布律为

P{X=k}=(1p)k1p,k=1,2,,P\{ X = k \} = (1 - p)^{k - 1} p, \quad k = 1, 2, \cdots,

其中 0<p<10 < p < 1,则称随机变量 XX 服从几何分布,记作 XG(p)X \sim G(p).

易得

(1p)k1p>0,k=0n(1p)k1p=1.(1 - p)^{k - 1} p > 0,\\ \sum_{k = 0}^n (1 - p)^{k - 1} p = 1.

几何分布的无记忆性:XG(p)X \sim G(p)n,mn, m
为任意的两个自然数,则

P{X>n+mX>n}=P{X>m}P\{ X > n + m \mid X > n \} = P\{ X > m \}

泊松(Poisson)分布:

如果随机变量 XX 的所有可能取值为 0,1,2,0, 1, 2, \cdots,并且

P{x=k}=λkeλk!,k=0,1,2,,P\{ x = k \} = \frac{\lambda^k \mathrm{e}^{-\lambda}}{k!}, \quad k = 0, 1, 2, \cdots,

其中 λ>0\lambda > 0 是常数,则称随机变量 XX 服从参数为 λ\lambda泊松分布,记作 XP(λ)X \sim P(\lambda)π(λ)\pi(\lambda).

易得

λkeλk!>0,k=0nλkeλk!=eλeλ=1.\frac{\lambda^k \mathrm{e}^{-\lambda}}{k!} > 0,\\ \sum_{k = 0}^n \frac{\lambda^k \mathrm{e}^{-\lambda}}{k!} = \mathrm{e}^{-\lambda} \cdot \mathrm{e}^{\lambda} = 1.

泊松定理:λ>0\lambda > 0 是常数,nn 为任意正整数,npn=λn p_n = \lambda,则对任一固定的非负整数 kk,有

limnCnkpk(1pn)nk=λkeλk!\lim_{n \to \infty} C_n^k p^k (1 - p_n)^{n - k} = \frac{\lambda^k \mathrm{e}^{-\lambda}}{k!}

泊松定理表明,当 nn 很大而 pp 很小时,有下面的近似公式

Cnkpk(1pn)nkλkeλk!C_n^k p^k (1 - p_n)^{n - k} \approx \frac{\lambda^k \mathrm{e}^{-\lambda}}{k!}

连续型随机变量及其概率密度

定义: 对于随机变量 XX 的分布函数 F(x)F(x),如果存在非负函数 f(x)f(x),使得对任意的 xx,都有

F(x)=xf(t)dt,F(x) = \int_{-\infty}^x f(t) \mathrm{d}t,

则称随机变量 XX连续型随机变量,其中函数 f(x)f(x) 叫作 XX概率密度函数,简称为概率密度,记作 Xf(x)X \sim f(x).

由定义可知,连续型随机变量的分布函数处处连续。

概率密度的性质:

性质 1: f(x)0f(x) \geq 0.

性质 2: +f(x)dx=1\int_{-\infty}^{+\infty} f(x) \mathrm{d}x = 1.

性质 3: 对于任意实数 a,b (a<b)a, b \ (a < b),有

P{a<Xb}=F(b)F(a)=abf(x)dx.P\{ a < X \leq b \} = F(b) - F(a) = \int_{a}^{b} f(x) \mathrm{d}x.

性质 4: 如果 f(x)f(x) 在点 xx 处连续,则有

F(x)=f(x).F'(x) = f(x).

性质 5: 对任意实数 aaP{X=a}=0P\{ X = a \} = 0. 即连续型随机变量取任意指定实数的概率均为零,因此有

P{x1<X<x2}=P{x1Xx2}=x1x2f(x)dxP\{ x_1 < X < x_2 \} = P\{ x_1 \leq X \leq x_2 \} = \int_{x_1}^{x_2} f(x) \mathrm{d}x

几何意义: 概率密度曲线总是位于 x 轴上方,并且介于它和 x 轴之间的面积等于 1;随机变量落在区间 (a,b](a, b] 的概率 P{a<Xb}P\{ a < X \leq b \} 等于区间 (a,b](a, b] 上曲线 y=f(x)y = f(x) 之下的曲边梯形的面积。

常用的连续型随机变量

均匀分布:

如果连续型随机变量 XX 的概率密度为

f(x)={1ba,a<x<b,0,Other,f(x) = \begin{cases} \dfrac{1}{b - a}, & a < x < b, \\ 0, & Other, \end{cases}

则称 XX 在区间 (a,b)(a, b) 上服从均匀分布,记作 XU(a,b)X \sim U(a, b). XX 的分布函数为

F(x)={0,a<x<b,xaba,ax<b,1,xb.F(x) = \begin{cases} 0, & a < x < b, \\ \dfrac{x - a}{b - a}, & a \leq x < b, \\ 1, & x \geq b. \end{cases}

对于任意 x1,x2(a,b)x_1, x_2 \in (a, b),若 x1<x2x_1 < x_2,则有

P{x1<X<x2}=x1x2f(x)dx=x1x21badx=x2x1ba.P\{ x_1 < X < x_2 \} = \int_{x_1}^{x_2} f(x) \mathrm{d}x = \int_{x_1}^{x_2} \frac{1}{b - a} \mathrm{d}x = \frac{x_2 - x_1}{b - a}.

这说明 XX 取值落在 (a,b)(a, b) 内任一子区间 (x1,x2)(x_1, x_2) 内的概率,只依赖于子区间的长度 x2x1x_2 - x_1,与子区间位置无关。

指数分布:

如果连续型随机变量 XX 的概率密度为

f(x)={λeλx,x>0,0,x0,f(x) = \begin{cases} \lambda \mathrm{e}^{-\lambda x}, & x > 0, \\ 0, & x \leq 0, \end{cases}

其中 λ>0\lambda > 0 为常数,则称 XX 服从参数为 λ\lambda指数分布XX 的分布函数为

F(x)={1eλx,x>0,0,x0.F(x) = \begin{cases} 1 - \mathrm{e}^{-\lambda x}, & x > 0, \\ 0, & x \leq 0. \end{cases}

指数分布无记忆性,例如:设随机变量 XX 服从参数 θ=1\theta = 1 的指数分布,则 P{2<X<3X1}=P{2<X<3}=e2e3P\{ 2 < X < 3 \mid X \geq 1 \} = P\{ 2 < X < 3 \} = \mathrm{e}^{-2} - \mathrm{e}^{-3}.

正态分布

定义: 如果随机变量 XX 具有概率密度

f(x)=12πσe(xμ)22σ2,<x<+,f(x) = \frac{1}{\sqrt{2\pi}\sigma} \mathrm{e}^{-\frac{(x - \mu)^2}{2\sigma^2}}, \quad -\infty < x < +\infty,

其中 μ,σ (σ>0)\mu, \sigma \ (\sigma > 0) 为常数,则称 XX 服从参数为 μ,σ\mu, \sigma正态分布,记作 XN(μ,σ2)X \sim N(\mu, \sigma^2).

XX 的分布函数为

F(x)=12πσxe(tμ)22σ2dt,<x<+,F(x) = \frac{1}{\sqrt{2\pi}\sigma} \int_{-\infty}^x \mathrm{e}^{-\frac{(t - \mu)^2}{2\sigma^2}} \mathrm{d}t, \quad -\infty < x < +\infty,

XX 的概率密度及分布函数的图像分别如下图所示:

非负性: f(x)0f(x) \geq 0.

规范性: +12πσe(xμ)22σ2dx=1\int_{-\infty}^{+\infty} \dfrac{1}{\sqrt{2\pi}\sigma} \mathrm{e}^{-\frac{(x - \mu)^2}{2\sigma^2}} \mathrm{d}x = 1.

正态分布的密度曲线特点:

  1. 概率密度曲线 y=f(x)y = f(x) 关于直线 x=μx = \mu 对称,以 xx 轴为水平渐近线,并在 x=μx = \mu 处取得最大值 12πσ\dfrac{1}{\sqrt{2\pi}\sigma}
  2. 在横坐标 x=μ±σx = \mu \pm \sigma 处有拐点;
  3. F(μ)=12F(\mu) = \dfrac{1}{2}.

正态分布的参数:

  • 位置参数 μ\mu:如果固定 σ\sigma,改变 μ\mu 的值,则概率密度曲线沿着 xx 轴平移,但形状不变,如下图所示:

  • 形状参数 σ\sigma:如果固定 μ\mu,改变 σ\sigma 的值,则当 σ\sigma 越小时概率密度曲线在 x=μx = \mu 附近越陡峭,XX 落在 x=μx = \mu 附近的概率越大;当 σ\sigma 越大时概率密度曲线越平坦。如下图所示:

正态分布的线性函数仍然服从正态分布,即若 XN(μ,σ2)X \sim N(\mu, \sigma^2),则 Y=aX+bN(aμ+b,(aσ)2)Y= aX + b \sim N\left( a\mu + b, (a\sigma)^2 \right).

证明:XX 的概率密度函数为

fX(x)=12πσe(xμ)22σ2,<x<+f_X(x) = \frac{1}{\sqrt{2\pi} \sigma} \mathrm{e}^{-\frac{(x - \mu)^2}{2\sigma^2}}, \quad -\infty < x < +\infty

y=g(x)=ax+by = g(x) = ax + b,由此解得 x=h(y)=ybax = h(y) = \dfrac{y - b}{a},且 h(y)=1ah'(y) = \dfrac{1}{a}.

Y=aX+b (a0)Y = aX + b \ (a \neq 0) 的概率密度函数为

fY(y)=1afX(yba),<y<+f_Y(y) = \frac{1}{|a|} f_X\left( \frac{y - b}{a} \right), \quad -\infty < y < +\infty

fY(y)=1a2πσe(ybaμ)22σ2=1a2πσe(y(b+aμ))22(aσ)2,<x<+f_Y(y) = \frac{1}{|a| \sqrt{2\pi} \sigma} \mathrm{e}^{-\frac{(\frac{y - b}{a} - \mu)^2}{2\sigma^2}} = \frac{1}{|a| \sqrt{2\pi} \sigma} \mathrm{e}^{-\frac{(y - (b + a\mu))^2}{2(a\sigma)^2}}, \quad -\infty < x < +\infty

因此有

Y=aX+bN(aμ+b,(aσ)2)Y = aX + b \sim N\left( a\mu + b, (a\sigma)^2 \right)

标准正态分布:

定义:XN(μ,σ2)X \sim N(\mu, \sigma^2),如果 μ=0,σ=1\mu = 0, \sigma = 1,则称 XX 服从标准正态分布,记作 XN(0,1)X \sim N(0, 1).

服从标准正态分布的随机变量 XX 的概率密度及分布函数分别记作 φ(x)\varphi(x)Φ(x)\varPhi(x),即

φ(x)=12πex22,<x<+,Φ(x)=12πxet22dt,<x<+.\varphi(x) = \frac{1}{\sqrt{2\pi}} \mathrm{e}^{-\frac{x^2}{2}}, \quad -\infty < x < +\infty,\\ \varPhi(x) = \frac{1}{\sqrt{2\pi}}\int_{-\infty}^x \mathrm{e}^{-\frac{t^2}{2}} \mathrm{d}t, \quad -\infty < x < +\infty.

它们的图形分别如下所示:

标准正态分布的性质:

  1. φ(x)\varphi(x) 是偶函数,φ(x)=φ(x),φ(0)=12π\varphi(-x) = \varphi(x), \varphi(0) = \dfrac{1}{\sqrt{2\pi}}.
  2. Φ(0)=12,Φ(x)=1Φ(x)\varPhi(0) = \dfrac{1}{2}, \varPhi(-x) = 1 - \varPhi(x).
  3. P{xα}=2Φ(a)1P\{ |x| \leq \alpha \} = 2\varPhi(a) - 1P{x>α}=1P{xα}=22Φ(a)P\{ |x| > \alpha \} = 1 - P\{ |x| \leq \alpha \} = 2 - 2\varPhi(a).

正态分布的标准化:XN(μ,σ2)X \sim N(\mu, \sigma^2),分布函数为 F(x)F(x),则令

Y=XμσN(0,1),F(x)=Φ(xμσ)Y = \frac{X - \mu}{\sigma} \sim N(0, 1),\\ F(x) = \varPhi\left( \frac{x - \mu}{\sigma} \right)

所以有

P{a<X<b}=F(b)F(a)=Φ(bμσ)Φ(aμσ)=P(aμσYbμσ)\begin{aligned} P\{ a < X < b \} &= F(b) - F(a) = \varPhi\left( \frac{b - \mu}{\sigma} \right) - \varPhi\left( \frac{a - \mu}{\sigma} \right)\\ &= P\left( \frac{a - \mu}{\sigma} \leq Y \leq \frac{b - \mu}{\sigma} \right) \end{aligned}

设随机变量 XX 服从正态分布 N(μ,σ2)N(\mu, \sigma^2),则随着 σ2\sigma^2 的增大,概率 P{Xμ<σ}P\{ |X - \mu| < \sigma\} 保持不变。

证:令 Y=XμσY = \dfrac{X - \mu}{\sigma},因为 XN(μ,σ2)X \sim N(\mu, \sigma^2),故 YN(0,1)Y \sim N(0, 1),从而有 P{Xμ<σ}=P{Y<1}P\{ |X - \mu| < \sigma\} = P\{ |Y| < 1 \} 是一个固定值。

标准正态分布上的 α\alpha 分位点:

定义:XN(0,1)X \sim N(0, 1). 对于给定的 α (0<α<1)\alpha \ (0 < \alpha < 1),如果 uαu_{\alpha} 满足条件

P{Xuα}=12πex22dx=α,P\{ X \geq u_{\alpha} \} = \frac{1}{\sqrt{2\pi}} \mathrm{e}^{-\frac{x^2}{2}} \mathrm{d}x = \alpha,

则称点 uαu_{\alpha} 为标准正态分布的α\alpha 分位点

Φ(uα)=1α\varPhi(u_\alpha) = 1 - \alphauα=u1α-u_\alpha = u_{1 - \alpha}.

常用数据:u0.05=1.645,u0.025=1.96u_{0.05} = 1.645, u_{0.025} = 1.96.

随机变量函数的分布

XX 是随机变量,函数 y=g(x)y = g(x),如果 Y=g(x)Y = g(x)YY 是随机变量 XX 的函数。已知 XX 的分布,考虑 YY 的分布。

方法:将与 YY 有关的事件转化成与 XX 有关的事件。

离散型随机变量的函数的分布: 确定随机变量 Y=g(X)Y = g(X) 的所有可能取值及相应概率,遇到多个 XX 取值造成同一个 YY 值,则将对应概率加和。

连续型随机变量的函数的分布:XX 为连续型随机变量,其概率密度为 fX(x)f_X(x),又 Y=g(X)Y = g(X),且 YY 也是连续型随机变量,求 YY 的概率密度 fY(y)f_Y(y).

方法一(分布函数法):先借助 XX 分布函数表达 YY 的分布函数 FY(y)F_Y(y),再对 yy 求导,得 YY 的概率密度函数

fY(y)={dFY(y)dy,可导,0,不可导.f_Y(y) = \begin{cases} \dfrac{\mathrm{d}F_Y(y)}{\mathrm{d}y}, & \text{可导}, \\ 0, & \text{不可导}. \end{cases}

方法二(公式法):连续型随机变量 XX 的概率密度 fX(x) (<x<+)f_X(x) \ (-\infty < x < +\infty),函数 g(x)g(x) 是处处可导的单调函数,则随机变量 Y=g(X)Y = g(X) 的概率密度为

fY(y)={fX(h(y))h(y),α<y<β,0,else,f_Y(y) = \begin{cases} f_X(h(y))|h'(y)|, & \alpha < y < \beta, \\ 0, & else, \end{cases}

其中 α=min{g(),g(+)},β=max{g(),g(+)}\alpha = \min\{ g(-\infty), g(+\infty) \}, \beta = \max\{ g(-\infty), g(+\infty) \}h(y)h(y)g(x)g(x) 的反函数。

f(x)f(x) 在有限区间 [a,b][a, b] 外为 0 时,只需考虑 g(x)g(x) 在区间内是否单调可导。

例(对数正态分布):如果 Y=lnXY = \ln{X} 服从正态部分 N(μ,σ2)N(\mu, \sigma^2),则称随机变量 XX 为服从参数为 μ,σ2\mu, \sigma^2 的对数正态分布。求对数正态分布的概率密度。

解:由 X=eY,Y=lnXN(μ,σ2)X = \mathrm{e}^Y, Y = \ln{X} \sim N(\mu, \sigma^2)

x>0x > 0 时,

FX(x)=P{Xx}=P{eYx}=P{Ylnx}=FY(lnx)F_X(x) = P\{ X \leq x \} = P\{ e^Y \leq x \} = P\{ Y \leq \ln{x} \} = F_Y(\ln{x})

x0x \leq 0 时,FX(x)=0F_X(x) = 0.

综上所述,

fx(x)=FX(x)={1xfY(lnx),x>0,0,x0,={12πσxe(lnxμ)22σ2,x>0,0,x0.f_x(x) = F_X'(x) = \begin{cases} \dfrac{1}{x} f_Y(\ln x), & x > 0, \\ 0, & x \leq 0, \end{cases} = \begin{cases} \dfrac{1}{\sqrt{2\pi} \sigma x} \mathrm{e}^{-\frac{(\ln{x} - \mu)^2}{2\sigma^2}}, & x > 0, \\ 0, & x \leq 0. \end{cases}

二维随机变量及其分布

二维随机变量及其分布函数

定义 1: 设随机试验 EE 的样本空间为 Ω\varOmegaXXYY 是定义在 Ω\varOmega 上的两个随机变量,由它们构成的向量 (X,Y)(X, Y) 称为二维随机变量二维随机向量

定义 2:(X,Y)(X, Y) 是二维随机变量。对于任意实数 xxyy,记事件 {Xx}\{ X \leq x \}{Yy}\{ Y \leq y \} 的交事件为 {Xx,Yy}\{ X \leq x, Y \leq y \},称二元函数

F(x,y)=P{Xx,Yy},(x,y)R2F(x, y) = P\{ X \leq x, Y \leq y \}, \quad (x, y) \in \R^2

二维随机变量 (X,Y)(X, Y) 的分布函数,或称为随机变量 XXYY联合分布函数

几何意义: 如果我们将二维随机变量 (X,Y)(X, Y) 看作 xOyxOy 平面上随机点的坐标,则分布函数 F(x,y)F(x, y) 在点 (x,y)(x, y) 处的函数值就是随机点落在以 (x,y)(x, y) 为顶点且位于该点左下方的无界域内的概率,如下图所示

分布函数 F(x,y)F(x, y) 具有以下性质:

性质 1: 0F(x,y)10 \leq F(x, y) \leq 1,且

F(,)=limxy(x,y)=0F(+,+)=limx+y+(x,y)=1F(-\infty, -\infty) = \lim_{x \to -\infty \atop y \to -\infty} (x, y) = 0\\ F(+\infty, +\infty) = \lim_{x \to +\infty \atop y \to +\infty} (x, y) = 1

对于任意固定的 xx,有 F(x,)=limy(x,y)=0F(x, -\infty) = \lim\limits_{y \to -\infty} (x, y) = 0
对于任意固定的 yy,有 F(,y)=limx(x,y)=0F(-\infty, y) = \lim\limits_{x \to -\infty} (x, y) = 0.

性质 2: F(x,y)F(x, y) 对于每个变量是单调不减函数,即对于任意固定的 yy,当 x1<x2x_1 < x_2 时,有

F(x1,y)F(x2,y)F(x_1, y) \leq F(x_2, y)

对于任意固定的 xx,当 y1<y2y_1 < y_2 时,有

F(x,y1)F(x,y2)F(x, y_1) \leq F(x, y_2)

性质 3: F(x,y)F(x, y) 关于 xx 右连续,关于 yy 右连续,即有

F(x+,y)=F(x,y),F(x,y+)=F(x,y)F(x^+, y) = F(x, y), \quad F(x, y^+) = F(x, y)

性质 4: 对于任意的 x1<x2,y1<y2x_1 < x_2, y_1 < y_2,有

P{x1<Xx2,y1<Yy2}=F(x2,y2)F(x1,y2)F(x2,y1)+F(x1,y1).P\{ x_1 < X \leq x_2, y_1 < Y \leq y_2 \}\\ = F(x_2, y_2) - F(x_1, y_2) - F(x_2, y_1) + F(x_1, y_1).

该公式用于利用分布函数计算概率。

边缘分布

定义: 设二维随机变量 (X,Y)(X, Y) 的分布函数为 F(x,y)F(x, y),记随机变量 XX 的分布函数为 FX(x)F_X(x),随机变量 YY 的分布函数为 FY(y)F_Y(y),分别称为二维随机变量 (X,Y)(X, Y) 关于 XX 和关于 YY边缘分布函数

我们有

FX(x)=P{Xx}=P{Xx,Y<+}=F(x,+),xR,FY(x)=P{Yy}=P{X<+,Yy}=F(+,y),yR.F_X(x) = P\{ X \leq x \} = P\{ X \leq x, Y < +\infty \} = F(x, +\infty), \quad x \in \R,\\ F_Y(x) = P\{ Y \leq y \} = P\{ X < +\infty, Y \leq y \} = F(+\infty, y), \quad y \in \R.

二维随机变量 (X,Y)(X, Y) 的边缘分布函数完全由联合分布函数 F(x,y)F(x, y) 确定;但是,一般地,由边缘分布函数不能确定联合分布函数。

二维随机变量的独立性

定义: 设二维随机变量 (X,Y)(X, Y) 的分布函数及其关于 XX 和关于 YY 的边缘分布函数为 F(x,y)F(x, y)Fx(x)F_x(x)FY(y)F_Y(y),如果对于任意实数 x,yx, y,都有

F(x,y)=FX(x)FY(y)F(x, y) = F_X(x) F_Y(y)

则称随机变量 XXYY相互独立的。

XXYY 相互独立就是对 x,yR\forall x, y \in R,两事件 {Xx}\{ X \leq x \}{Yy}\{ Y \leq y \} 相互独立。

XXYY 相互独立,对任意连续函数 h,gh, g,则 h(X)h(X)g(Y)g(Y) 也相互独立。

二维离散型随机变量及其概率分布

定义 1: 若二维随机变量 (X,Y)(X, Y) 所有可能取得值是有限对或可列无限对,则称 (X,Y)(X, Y)二维离散型随机变量

(X,Y)(X, Y) 所有可能取的值为 (xi,yj) (i,j=1,2,)(x_i, y_j) \ (i, j = 1, 2, \cdots),记事件 {X=xi}\{ X = x_i \} 与事件 {Y=yj}\{ Y = y_j \} 的交事件为 {X=xi,Y=yj} (i,j=1,2,)\{ X = x_i, Y = y_j \} \ (i, j = 1, 2, \cdots).

定义 2:(X,Y)(X, Y) 所有可能取的值为 (xi,yj) (i,j=1,2,)(x_i, y_j) \ (i, j = 1, 2, \cdots),如果

P{X=xi,Y=yj}=pij,i,j=1,2,,P\{ X = x_i, Y = y_j \} = p_{ij}, \quad i, j = 1, 2, \cdots,

且满足

  1. pij0, (i,j=1,2,)p_{ij} \geq 0, \ (i, j = 1, 2, \cdots)
  2. i=1j=1pij=1\sum_{i = 1}^\infty \sum_{j = 1}^\infty p_{ij} = 1

则称上式为二维离散型随机变量 (X,Y)(X, Y)概率分布,或称为随机变量 XX 和随机变量 YY联合概率分布联合分布律

(X,Y)(X, Y) 的概率分布可以用如下的表格表示:

边缘概率分布: 设二维离散型随机变量 (X,Y)(X, Y) 的概率分布为

P{X=xi,Y=yj}=pij,i,j=1,2,P\{ X = x_i, Y = y_j \} = p_{ij}, \quad i, j = 1, 2, \cdots

随机变量 XXYY 的概率分布

P{X=xi}=pi,i=1,2,,P{Y=yj}=pj,j=1,2,,P\{ X = x_i \} = p_{i \cdot}, \quad i = 1, 2, \cdots,\\ P\{ Y = y_j \} = p_{\cdot j}, \quad j = 1, 2, \cdots,

分别称为 (X,Y)(X, Y) 关于 XX 和关于 YY 的边缘概率分布或边缘分布律

边缘分布律都是一维随机变量的分布;联合分布律可确定边缘分布律,反之不一定。

随机变量的独立性: 设二维离散型随机变量 (X,Y)(X, Y)XXYY 相互独立的充要条件是对任意的 i,ji, j,都有

P{X=xi,Y=yj}=P{X=xi}P{Y=yj},i,j=1,2,,P\{ X = x_i, Y = y_j \} = P\{ X = x_i \}P\{ Y = y_j \}, \quad i, j = 1, 2, \cdots,

pij=pipji,j=1,2,.p_{ij} = p_{i \cdot} p_{\cdot j} \quad i, j = 1, 2, \cdots.

二维连续型随机变量及其概率密度

定义: 设二维随机变量 (X,Y)(X, Y) 的分布函数为 F(x,y)F(x, y),如果存在非负的二元函数 f(x,y)f(x, y),对于任意实数 x,yx, y,有

F(x,y)=xyf(u,v)dudvF(x, y) = \int_{-\infty}^x \int_{-\infty}^y f(u, v) \mathrm{d}u\mathrm{d}v

则称 (X,Y)(X, Y)二维连续型随机变量f(x,y) ((x,y)R2)f(x, y) \ ((x, y) \in \R^2) 称为二维连续型随机变量 (X,Y)(X, Y)概率密度,或称为随机变量 XXYY联合概率密度

概率密度 f(x,y)f(x, y) 具有下列性质:

性质 1: f(x,y)0f(x, y) \geq 0.

性质 2: ++f(x,y)dxdy=1\int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty} f(x, y) \mathrm{d}x\mathrm{d}y = 1.

性质 3: 如果 f(x,y)f(x, y) 在点 (x,y)(x, y) 处连续,则

f(x,y)=2F(x,y)xy.f(x, y) = \frac{\partial^2{F(x, y)}}{\partial{x}\partial{y}}.

性质 4:GGxOyxOy 面上的一个区域,则

P{(X,Y)G}=Gf(x,y)dxdy.P\{ (X, Y) \in G \} = \iint_G f(x, y) \mathrm{d}x\mathrm{d}y.

几何意义: z=f(x,y)z = f(x, y) 表示空间 OxyzOxyz 中的一张曲面。性质 1 和性质 2 表明,曲面 z=f(x,y)z = f(x, y) 位于 xOyxOy 平面上方,介于它和 xOyxOy 平面之间的体积为 1. 性质 4 表明,随机点 (X,Y)(X, Y) 落在区域 GG 内的概率 P{(X,Y)G}P\{ (X, Y) \in G \} 等于以 GG 为底、以曲面 z=f(x,y)z = f(x, y) 为顶的曲顶柱体体积的数值。

边缘概率密度: 设二维连续型随机变量 (X,Y)(X, Y) 的概率密度为 f(x,y),(x,y)R2f(x, y), (x, y) \in \R^2,将一元函数

fX(x)=+f(x,y)dy,xR,fY(y)=+f(x,y)dx,yRf_X(x) = \int_{-\infty}^{+\infty} f(x, y) \mathrm{d}y, \quad x \in \R, \\ f_Y(y) = \int_{-\infty}^{+\infty} f(x, y) \mathrm{d}x, \quad y \in \R

分别称为二维随机变量 (X,Y)(X, Y) 关于 XX 和关于 YY边缘概率密度

随机变量的独立性: 设二维连续型随机变量 (X,Y)(X, Y)XXYY 相互独立的充要条件是对任意实数 x,yx, y,都有

f(x,y)=fX(x)fY(y),(<x,y<+).f(x, y) = f_X(x) f_Y(y), \quad (-\infty < x, y < +\infty).

二维均匀分布与正态分布

二维均匀分布:

DDxOyxOy 面上的有界区域,其面积为 AA. 如果二维随机变量 (X,Y)(X, Y) 具有概率密度

f(x,y)={1A,(x,y)D,0,else,f(x, y) = \begin{cases} \dfrac{1}{A}, & (x, y) \in D, \\ 0, & else, \end{cases}

则称 (X,Y)(X, Y) 在区域 DD 上服从均匀分布

例:设二维随机变量 (X,Y)(X, Y) 在区域 D={(x,y)0yx,0x1}D = \{ (x, y) \mid 0 \leq y \leq x, 0 \leq x \leq 1 \} 上服从均匀分布。

(1)求 (X,Y)(X, Y) 的概率密度。

解:区域 DD 的面积为 S=12S = \dfrac{1}{2},故 (X,Y)(X, Y) 的概率密度为

f(x,y)={2,(x,y)D,0,else.f(x, y) = \begin{cases} 2, & (x, y) \in D, \\ 0, & else. \end{cases}

(2)求 P{12<X<34,0<Y<34}P\left\{ \dfrac{1}{2} < X < \dfrac{3}{4}, 0 < Y < \dfrac{3}{4} \right\}.

P{12<X<34,0<Y<34}=梯形f(x,y)dxdy=1234dxox2dy=516\begin{aligned} P\left\{ \dfrac{1}{2} < X < \dfrac{3}{4}, 0 < Y < \dfrac{3}{4} \right\} &= \iint\limits_{\text{梯形}} f(x, y) \mathrm{d}x\mathrm{d}y \\ &= \int_{\frac{1}{2}}^{\frac{3}{4}} \mathrm{d}x \int_o^x 2 \mathrm{d}y \\ &= \frac{5}{16} \end{aligned}

二维正态分布:

设二维随机变量 (X,Y)(X, Y) 的概率密度为

f(x,y)=12πσ1σ21ρ2exp{12(1ρ2)[(xμ1)2σ122ρ(xμ1)(yμ2)σ1σ2+(yμ2)2σ22]},(x,y)R2\begin{aligned} f(x, y) = & \frac{1}{2\pi \sigma_1 \sigma_2 \sqrt{1 - \rho^2}} \cdot \exp\bigg\{ \frac{-1}{2(1 - \rho^2)} \bigg[ \frac{(x - \mu_1)^2}{\sigma_1^2} \\ & - 2\rho \frac{(x - \mu_1)(y - \mu_2)}{\sigma_1 \sigma_2} + \frac{(y - \mu_2)^2}{\sigma_2^2} \bigg] \bigg\}, \quad (x, y) \in \R^2 \end{aligned}

其中 μ1,μ2,σ1,σ2,ρ\mu_1, \mu_2, \sigma_1, \sigma_2, \rho 都是常数,且 σ1>0,σ2>0,1<ρ<1\sigma_1 > 0, \sigma_2 > 0, -1 < \rho < 1,则称 (X,Y)(X, Y) 服从参数为 μ1,μ2,σ1,σ2,ρ\mu_1, \mu_2, \sigma_1, \sigma_2, \rho二维正态分布,记作 (X,Y)N(μ1,μ2,σ12,σ22,ρ)(X, Y) \sim N(\mu_1, \mu_2, \sigma_1^2, \sigma_2^2, \rho).

二维正态分布的性质:

  1. (X,Y)N(μ1,μ2,σ12,σ22,ρ)(X, Y) \sim N(\mu_1, \mu_2, \sigma_1^2, \sigma_2^2, \rho),则 XN(μ1,σ12)X \sim N(\mu_1, \sigma_1^2)YN(μ2,σ22)Y \sim N(\mu_2, \sigma_2^2),即二维正态分布的边缘分布式一维正态分布。

只要给定二维正态分布的前 4 个参数,关于 XX 和关于 YY 的边缘分布也就确定了,和 ρ\rho 无关。

  1. XXYY 相互独立的充要条件是 ρ=0\rho = 0.

独立正态随机变量的非零线性组合仍服从正态分布,例如若 X,YX, Y 相互独立,XN(μ1,σ12),YN(μ2,σ22)X \sim N(\mu_1, \sigma_1^2), Y \sim N(\mu_2, \sigma_2^2),则

  • X+YN(μ1+μ2,σ12+σ22)X + Y \sim N(\mu_1 + \mu_2, \sigma_1^2 + \sigma_2^2)
  • aX+bYN(aμ1+bμ2,a2σ12+b2σ22)aX + bY \sim N(a\mu_1 + b\mu_2, a^2 \sigma_1^2 + b^2 \sigma_2^2)
  • XYN(μ1μ2,σ12+σ22)X - Y \sim N(\mu_1 - \mu_2, \sigma_1^2 + \sigma_2^2).

条件分布

离散型随机变量的条件分布:(X,Y)(X, Y) 是二维离散型随机变量,对于固定的 jj,如果 P{Y=yj}=pj>0P\{ Y = y_j \} = p_{\cdot j} > 0,则称

P{X=xiY=yj}=P{X=xi,Y=yj}P{Y=yj}=pijpj,i=1,2,P\{ X = x_i \mid Y = y_j \} = \frac{P\{X = x_i, Y = y_j\}}{P\{ Y = y_j \}} = \frac{p_{ij}}{p_{\cdot j}}, \quad i = 1, 2, \cdots

为在条件 Y=yjY = y_j 下随机变量 XX条件概率分布。对于固定的 ii,若 pi>0p_{i \cdot} > 0,则称

P{Y=yjX=xi}=P{X=xi,Y=yj}P{X=xi}=pijpi,j=1,2,P\{ Y = y_j \mid X = x_i \} = \frac{P\{X = x_i, Y = y_j\}}{P\{ X = x_i \}} = \frac{p_{ij}}{p_{i \cdot}}, \quad j = 1, 2, \cdots

为在条件 X=xiX = x_i 下随机变量 YY条件概率分布

连续型随机变量的条件分布:

Y=yY = y 条件下,XX条件概率密度

fXY(xy)=f(x,y)fY(y)f_{X \mid Y}(x \mid y) = \frac{f(x, y)}{f_Y(y)}

XX条件分布函数

FXY(xy)=xfXY(xy)dxF_{X \mid Y}(x \mid y) = \int_{-\infty}^x f_{X \mid Y}(x \mid y) \mathrm{d}x

X=xX = x 条件下,YY条件概率密度

fYX(yx)=f(x,y)fX(x)f_{Y \mid X}(y \mid x) = \frac{f(x, y)}{f_X(x)}

YY条件分布函数

FYX(yx)=yfYX(yx)dyF_{Y \mid X}(y \mid x) = \int_{-\infty}^y f_{Y \mid X}(y \mid x) \mathrm{d}y

二维随机变量的函数的分布

(X,Y)(X, Y) 为二维随机变量,g(x,y)g(x, y) 为二元函数,则一维随机变量 Z=g(X,Y)Z = g(X, Y) 是二维随机变量 (X,Y)(X, Y) 的函数。

二维离散型随机变量的函数的分布:

如果 XXYY 相互独立,则有

XP(λ1),YP(λ2)X+YP(λ1+λ2)XB(m,p),YB(n,p)X+YB(m+n,p)X \sim P(\lambda_1), Y \sim P(\lambda_2) \Rightarrow X + Y \sim P(\lambda_1 + \lambda_2)\\ X \sim B(m, p), Y \sim B(n, p) \Rightarrow X + Y \sim B(m + n, p)

二维连续型随机变量的函数的分布:

已知 (X,Y)(X, Y) 概率密度,求 Z=g(X,Y)Z = g(X, Y) 的概率密度。

1. 一般情形:

  1. 先求随机变量 ZZ 的分布函数

FZ(z)=P{Zz}=P{g(X,Y)z}=g(x,y)zf(x,y)dxdy\begin{aligned} F_Z(z) &= P\{ Z \leq z \} = P\{ g(X, Y) \leq z \} \\ &= \iint\limits_{g(x, y) \leq z} f(x, y) \mathrm{d}x\mathrm{d}y \end{aligned}

  1. FZ(z)F_Z(z) 求导,得 ZZ 的概率密度函数为

fZ(z)=FZ(z)f_Z(z) = F_Z'(z)

2. Z=X+YZ = X + Y 的概率密度:

fZ(z)=+f(zy,y)dyf_Z(z) = \int_{-\infty}^{+\infty} f(z - y, y) \mathrm{d}y

fZ(z)=+f(x,zx)dxf_Z(z) = \int_{-\infty}^{+\infty} f(x, z - x) \mathrm{d}x

X,YX, Y 相互独立,则 Z=X+YZ = X + Y 的密度函数为

fZ(z)=+fX(zy)fY(y)dyf_Z(z) = \int_{-\infty}^{+\infty} f_X(z - y) f_Y(y) \mathrm{d}y

fZ(z)=+fX(x)fY(zx)dyf_Z(z) = \int_{-\infty}^{+\infty} f_X(x) f_Y(z - x) \mathrm{d}y

上式称为卷积公式

3. XXYY 相互独立,M=maxX,YM = \max{X, Y} 的分布函数:

Fmax(z)=P{Mz}=P{Xz,Yz}=P{Xz}P{Yz}=FX(z)FY(z)\begin{aligned} F_{\max}(z) &= P\{ M \leq z \} = P\{ X \leq z, Y \leq z \} \\ &= P\{ X \leq z \} P\{ Y \leq z \} \\ &= F_X(z) F_Y(z) \end{aligned}

4. XXYY 相互独立,N=minX,YN = \min{X, Y} 的分布函数:

Fmin(z)=P{Nz}=1P{N>z}=1P{X>z,Y>z}=1P{X>z}P{Y>z}=1[1P{Xz}][1P{Yz}]=1[1FX(z)][1FY(z)]\begin{aligned} F_{\min}(z) &= P\{ N \leq z \} = 1 - P\{ N > z \} \\ &= 1 - P\{ X > z, Y > z \} \\ &= 1 - P\{ X > z \} P\{ Y > z \} \\ &= 1 - [1 - P\{ X \leq z \}] [1 - P\{ Y \leq z \}] \\ &= 1 - [1 - F_X(z)] [1 - F_Y(z)] \end{aligned}

推广:nn 个随机变量 X1,X2,,XnX_1, X_2, \cdots, X_n 相互独立,分布函数分别为 FX1(x1),FX2(x2),,FXn(xn)F_{X_1}(x_1), F_{X_2}(x_2), \cdots, F_{X_n}(x_n),则 M=max{X1,X2,,Xn}M = \max\{ X_1, X_2, \cdots, X_n \} 的分布函数为

Fmax(z)=FX1(z)FX2(z)FXn(z)F_{\max}(z) = F_{X_1}(z) F_{X_2}(z) \cdots F_{X_n}(z)

N=min{X1,X2,,Xn}N = \min\{ X_1, X_2, \cdots, X_n \} 的分布函数为

Fmin(z)=1[1FX1(z)][1FX2(z)][1FXn(z)]F_{\min}(z) = 1 - [1 - F_{X_1}(z)] [1 - F_{X_2}(z)] \cdots [1 - F_{X_n}(z)]

随机变量的数字特征

数学期望

数学期望是一个常数,它是一种加权平均,表示随机变量所有取值的平均值,具有重要的统计意义。

定义: 设离散型随机变量 XX 的概率分布为 P{X=xk}=pk,k=1,2,P\{ X = x_k \} = p_k, k = 1, 2, \cdots,如果无穷级数 k=1xkpk\sum_{k = 1}^\infty x_k p_k 绝对收敛,则称无穷级数 k=1xkpk\sum_{k = 1}^\infty x_k p_k 的和为离散型随机变量 XX数学期望均值,记作 E(X)E(X)EXEX,即

E(X)=k=1xkpk.E(X) = \sum_{k = 1}^\infty x_k p_k.

设连续型随机变量 XX 的概率密度为 f(x)f(x),如果反常积分 +xf(x)dx\int_{-\infty}^{+\infty} xf(x) \mathrm{d}x 绝对收敛,则称 +xf(x)dx\int_{-\infty}^{+\infty} xf(x) \mathrm{d}x 的值为连续型变量 XX 的数学期望或均值,即

E(X)=+xf(x)dx.E(X) = \int_{-\infty}^{+\infty} xf(x) \mathrm{d}x.

随机变量的函数的数学期望:

定理 1: 设随机变量 YY 是随机变量 XX 的函数:Y=g(X)Y = g(X),其中 gg 是一元连续函数。

  1. XX 是离散型随机变量,其概率分布为 P{X=xk}=pk,k=1,2,P\{ X = x_k \} = p_k, k = 1, 2, \cdots,如果无穷级数 k=1xkpk\sum_{k = 1}^\infty x_k p_k 绝对收敛,则随机变量 YY 的数学期望为

E(Y)=E[g(X)]=k=1g(xk)pkE(Y) = E[g(X)] = \sum_{k = 1}^\infty g(x_k) p_k

  1. XX 是连续型随机变量,其概率分布为 f(x)f(x),如果反常积分 +g(x)f(x)dx\int_{-\infty}^{+\infty} g(x)f(x) \mathrm{d}x 绝对收敛,则随机变量 YY 的数学期望为

E(Y)=E[g(X)]=+g(x)f(x)dxE(Y) = E[g(X)] = \int_{-\infty}^{+\infty} g(x)f(x) \mathrm{d}x

根据这一定理求随机变量 Y=g(X)Y = g(X) 的数学期望时,只需直到 XX 的分布,无需求出 YY 的分布。

定理 2:ZZ 是随机变量 XXYY 的函数:Z=g(X,Y)Z = g(X, Y),其中 gg 是二元连续函数。

  1. X,YX, Y 的联合分布律为 P{X=xi,Y=yj}=pij,i,j=1,2,P\{ X = x_i, Y = y_j \} = p_{ij}, i, j = 1, 2, \cdots,如果 i=1j=1g(xi,yj)pij\sum_{i = 1}^\infty\sum_{j = 1}^\infty g(x_i, y_j) p_{ij} 绝对收敛,则

E(Z)=E[g(X,Y)]=j=1i=1g(xi,yj)pijE(Z) = E[g(X, Y)] = \sum_{j = 1}^\infty \sum_{i = 1}^\infty g(x_i, y_j) p_{ij}

  1. X,YX, Y 的联合概率密度为 f(x,y)f(x, y),如果广义积分 ++g(x,y)f(x,y)dxdy\int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty} g(x, y)f(x, y) \mathrm{d}x\mathrm{d}y 绝对收敛,则

E(Z)=E[g(X,Y)]=++g(x,y)f(x,y)dxdyE(Z) = E[g(X, Y)] = \int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty} g(x, y)f(x, y) \mathrm{d}x\mathrm{d}y

数学期望的性质:

性质 1:cc 为常数,E(c)=cE(c) = c.

性质 2:cc 为常数,E(cX)=cE(X)E(cX) = cE(X).

性质 3: E(X±Y)=E(X)±E(Y)E(X \pm Y) = E(X) \pm E(Y).

性质 4:X,YX, Y 相互独立,则 E(XY)=E(X)E(Y)E(XY) = E(X) E(Y).

性质 5(柯西-施瓦茨不等式): 设随机变量 X,YX, Y,若 E(X2)E(X^2)E(Y2)E(Y^2) 都存在,则

[E(XY)]2E(X2)E(Y2).[E(XY)]^2 \leq E(X^2) E(Y^2).

方差

在许多问题中,除了要直到随机变量的均值外,还要研究随机变量与其均值之间的偏离程度,因此引入了方差的概念。

定义:XX 是一个随机变量,如果 E{[XE(X)]2}E\{ [X - E(X)]^2 \} 存在,则称之为随机变量 XX方差,记作 D(X)D(X)DXDX,即

D(X)=E{[XE(X)]2}.D(X) = E\{ [X - E(X)]^2 \}.

D(X)\sqrt{D(X)} 为随机变量 XX标准差均方差,记作 σ(X)\sigma(X),即

σ(X)=D(X).\sigma(X) = \sqrt{D(X)}.

随机变量 XX 的方差反映了 XX 与其数学期望 E(X)E(X) 的偏离程度。如果 XX 取值集中在 E(X)E(X) 附近,则 D(X)D(X) 较小;如果 XX 取值比较分散,则 D(X)D(X) 较大。

离散型随机变量:

D(X)=k=1[xkE(X)]2pkD(X) = \sum_{k = 1}^\infty [x_k - E(X)]^2 p_k

连续型随机变量:

D(X)=+[xE(X)]2f(x)dxD(X) = \int_{-\infty}^{+\infty} [x - E(X)]^2 f(x) \mathrm{d}x

方差的性质:

性质 1:cc 为常数,则 D(c)=0D(c) = 0.

性质 2: D(cX)=c2D(X)D(cX) = c^2 D(X)

性质 3: D(X+c)=D(X)D(X + c) = D(X)

性质 4:XXYY 相互独立,则 D(X±Y)=D(X)+D(Y)D(X \pm Y) = D(X) + D(Y).

性质 5: 随机变量 XX 的方差 D(X)=0D(X) = 0 的充分必要条件是 XX 以概率 1 取常数 C=E(X)C = E(X),即 P{X=C}=1P\{ X = C \} = 1.

常见分布的数学期望与方差:

分布 数学期望 方差
(0-1)分布 pp p(1p)p(1 - p)
二项分布 B(n,p)B(n, p) npnp np(1p)np(1 - p)
泊松分布 P(λ)P(\lambda) λ\lambda λ\lambda
几何分布 1/p1 / p (1p)/p2(1 - p) / p^2
均匀分布 U(a,b)U(a, b) (a+b)/2(a + b) / 2 (ba)2/12(b - a)^2 / 12
正态分布 N(μ,σ2)N(\mu, \sigma^2) μ\mu σ2\sigma^2
指数分布 E(λ)E(\lambda) 1/λ1 / \lambda 1/λ21 / \lambda^2

随机变量的标准化: 设随机变量 XX 具有数学期望 E(X)=μE(X) = \mu 及方差 D(X)=σ2>0D(X) = \sigma^2 > 0,则称

X=XμσX^* = \frac{X - \mu}{\sigma}

XX标准化随机变量

显然,E(X)=0E(X^*) = 0D(X)=1D(X^*) = 1.

XN(μ,σ2) (σ>0)X \sim N(\mu, \sigma^2) \ (\sigma > 0),则

X=XμσN(0,1).X^* = \frac{X - \mu}{\sigma} \sim N(0, 1).

协方差与相关系数

定义: 设随机变量 XXYY 的数学期望 E(X)E(X)E(Y)E(Y) 都存在,称 E{[XE(X)][YE(Y)]}E\{ [X - E(X)][Y - E(Y)] \}XXYY协方差,记作

Cov(X,Y)=E{[XE(X)][YE(Y)]}Cov(X, Y) = E\{ [X - E(X)][Y - E(Y)] \}

将上式展开,易得

Cov(X,Y)=E(XY)E(X)E(Y)Cov(X, Y) = E(XY) - E(X)E(Y)

这是常用的协方差计算公式。

Cov(X,X)=D(X)Cov(X, X) = D(X)

协方差的性质:

性质 1: Cov(X,Y)=Cov(Y,X)Cov(X, Y) = Cov(Y, X)

性质 2: Cov(aX,bY)=abCov(X,Y)Cov(aX, bY) = abCov(X, Y)

性质 3: Cov(X+Y,Z)=Cov(X,Z)+Cov(Y,Z)Cov(X + Y, Z) = Cov(X, Z) + Cov(Y, Z)

性质 4: D(X±Y)=D(X)+D(Y)±2Cov(X,Y)D(X \pm Y) = D(X) + D(Y) \pm 2Cov(X, Y)

定义: 设随机变量 XXYY 的方差都存在且不等于零,协方差 Cov(X,Y)Cov(X, Y) 存在,称

ρXY=Cov(X,Y)D(X)D(Y)\rho_{XY} = \dfrac{Cov(X, Y)}{\sqrt{D(X)} \sqrt{D(Y)}}

X,YX, Y相关系数。当 ρXY=0\rho_{XY} = 0 时,称 XXYY 不相关

相关系数的常用计算公式为

ρXY=E(XY)E(X)E(Y)D(X)D(Y).\rho_{XY} = \dfrac{E(XY) - E(X)E(Y)}{\sqrt{D(X)} \sqrt{D(Y)}}.

相关系数的性质:

性质 1: ρXY1|\rho_{XY}| \leq 1

性质 2: 如果 XXYY 相互独立,则 ρXY=0|\rho_{XY}| = 0.

性质 3: ρXY=1|\rho_{XY}| = 1 的充分必要条件是:存在常数 a,ba, b,使得

P{Y=a+bX}=1.P\{ Y = a + bX \} = 1.

相关系数的含义: 相关系数表示了随机变量 XXYY 的线性相关的程度:

  1. ρXY=1|\rho_{XY}| = 1 时,X,YX, Y 几乎就是线性关系;
  2. ρXY|\rho_{XY}| 较大时,X,YX, Y 的线性关系较为密切;
  3. ρXY|\rho_{XY}| 较小时,X,YX, Y 的线性关系较弱;
  4. ρXY=0\rho_{XY} = 0 时,X,YX, Y 之间没有线性关系(不相关)。

ρXY>0\rho_{XY} > 0,则 X,YX, Y 正相关;若 ρXY<0\rho_{XY} < 0,则 X,YX, Y 负相关。

对于随机变量 X,YX, Y,以下命题等价:

  1. Cov(X,Y)=0Cov(X, Y) = 0
  2. XXYY 不相关;
  3. E(XY)=E(X)E(Y)E(XY) = E(X) E(Y)
  4. D(X+Y)=D(X)+D(Y)D(X + Y) = D(X) + D(Y).

这里假设上述命题中出现的数字特征都是存在的。

二维正态分布的结论: 如果 (X,Y)N(μ1,μ2,σ12,σ22,ρ)(X, Y) \sim N(\mu_1, \mu_2, \sigma_1^2, \sigma_2^2, \rho),则

  1. ρ=ρXY\rho = \rho_{XY}

即二维正态随机变量 (X,Y)(X, Y) 的分布完全由 XXYY 的数学期望,方差以及 XXYY 的相关系数所确定。

  1. XXYY 相互独立     \iff XXYY 不相关;

  2. Cov(X,Y)=ρσ1σ2Cov(X, Y) = \rho \sigma_1 \sigma_2.

定义: 对于随机变量 X,YX, Y

  1. E(Xk)=μk (k=1,2,)E(X^k) = \mu_k \ (k = 1, 2, \cdots)XXkk 阶原点矩
  2. E{[XE(X)]k}E\{ [X - E(X)]^k \}XXkk 阶中心矩
  3. E(XkYl) (k,l=1,2,)E(X^k Y^l) \ (k, l = 1, 2, \cdots)XXYYk+lk + l 阶混合原点矩
  4. E{[XE(X)]k[YE(Y)]l}E\{ [X - E(X)]^k [Y - E(Y)]^l \}XXYYk+lk + l 阶混合中心矩

显然,随机变量 XX 的数学期望 E(X)E(X)XX 的一阶原点矩,方差 D(X)D(X)XX 的二阶中心矩,随机变量 XXYY 的协方差 Cov(X,Y)Cov(X, Y)XXYY 的二阶混合中心矩。

协方差矩阵:

对于二维随机变量 (X1,X2)(X_1, X_2),记

c11=E{(X1EX1)2}=D(X1)c12=E{(X1EX1)(X2EX2)}=Cov(X1,X2)c21=E{(X2EX2)(X1EX1)}=Cov(X2,X1)c22=E{(X2EX2)2}=Cov(X2,X1)\begin{aligned} c_{11} &= E\{ (X_1 - EX_1)^2 \} = D(X_1) \\ c_{12} &= E\{ (X_1 - EX_1)(X_2 - EX_2) \} = Cov(X_1, X_2) \\ c_{21} &= E\{ (X_2 - EX_2)(X_1 - EX_1) \} = Cov(X_2, X_1) \\ c_{22} &= E\{ (X_2 - EX_2)^2 \} = Cov(X_2, X_1) \\ \end{aligned}

如果它们都存在,写成矩阵的形式

C=[c11c12c21c22]\bm{C} = \begin{bmatrix} c_{11} & c_{12} \\ c_{21} & c_{22} \end{bmatrix}

称此矩阵为二维随机变量 (X1,X2)(X_1, X_2)协方差矩阵

nn 维随机变量 (X1,X2,,Xn)(X_1, X_2, \cdots, X_n) 关于 X1,X2,,XnX_1, X_2, \cdots, X_n 的二阶中心矩和二阶混合中心矩

cij=E{[XiE(Xi)][XjE(Xj)]},i,j=1,2,,nc_{ij} = E\{ [X_i - E(X_i)][X_j - E(X_j)] \}, \quad i, j = 1, 2, \cdots, n

都存在,则称矩阵

C=[c11c12c1nc21c22c2ncn1cn2cnn]\bm{C} = \begin{bmatrix} c_{11} & c_{12} & \cdots & c_{1n} \\ c_{21} & c_{22} & \cdots & c_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ c_{n1} & c_{n2} & \cdots & c_{nn} \\ \end{bmatrix}

nn 维随机变量 (X1,X2,,Xn)(X_1, X_2, \cdots, X_n)协方差矩阵。由于 cij=cji (ij,i,j=1,2,,n)c_{ij} = c_{ji} \ (i \neq j, i, j = 1, 2, \cdots, n),所以 C\bm{C} 是对称矩阵。

二维正态分布的矩阵表达式:(X,Y)N(μ1,μ2,σ12,σ22,ρ)(X, Y) \sim N(\mu_1, \mu_2, \sigma_1^2, \sigma_2^2, \rho),密度函数也表示为

f(x,y)=1(2π)2/2C1/2exp{12(Xμ)TC1(Xμ)}f(x, y) = \frac{1}{(2\pi)^{2 / 2}|\bm{C}|^{1/2}}\exp\left\{-\frac{1}{2} (X - \mu)^T \bm{C}^{-1}(X - \mu)\right\}

其中 C=[σ12ρσ1σ2ρσ1σ2σ22]\bm{C} = \begin{bmatrix} \sigma_1^2 & \rho \sigma_1 \sigma_2 \\ \rho \sigma_1 \sigma_2 & \sigma_2^2 \end{bmatrix} 为协方差矩阵,行列式为

C=σ12σ22(1ρ2)|\bm{C}| = \sigma_1^2 \sigma_2^2 (1 - \rho^2)

并且

C1=1C[σ22ρσ1σ2ρσ1σ2σ12],X=[xy],μ=[μ1μ2]\bm{C}^{-1} = \frac{1}{|\bm{C}|}\begin{bmatrix} \sigma_2^2 & -\rho \sigma_1 \sigma_2 \\ -\rho \sigma_1 \sigma_2 & \sigma_1^2 \end{bmatrix}, X = \begin{bmatrix} x \\ y \end{bmatrix}, \mu = \begin{bmatrix} \mu_1 \\ \mu_2 \end{bmatrix}

n 维正态随机变量的性质:

  1. nn 维随机变量 (X1,X2,,Xn)(X_1, X_2, \cdots, X_n) 服从 nn 维正态分布的充分必要条件是:X1,X2,,XnX_1, X_2, \cdots, X_n 的任意线性组合

    k1X1+k2X2++knXnk_1 X_1 + k_2 X_2 +\cdots + k_n X_n

都服从一维正态分布,其中 k1,k2,,knk_1, k_2, \cdots, k_n 是不全为零的常数。

  1. (X1,X2,,Xn)(X_1, X_2, \cdots, X_n) 服从 nn 维正态分布,如果 Y1,Y2,,YmY_1, Y_2, \cdots, Y_mXi (i=1,2,,n)X_i \ (i = 1, 2, \cdots, n) 的线性函数,则 Y1,Y2,,YmY_1, Y_2, \cdots, Y_m 也服从多维正态分布。

  2. 如果 (X1,X2,,Xn)(X_1, X_2, \cdots, X_n) 服从 nn 维正态分布,则 “随机变量 X1,X2,,XnX_1, X_2, \cdots, X_n 相互独立” 与 “X1,X2,,XnX_1, X_2, \cdots, X_n 两两不相关” 等价。

大数定律与中心极限定理

切比雪夫不等式

定理: 设随机变量 XX 具有数学期望 E(X)=μE(X) = \mu 和方差 D(X)=σ2D(X) = \sigma^2,则对于任意给定的整数 ε\varepsilon,有

P{XE(X)ε}D(X)ε2P\{ |X - E(X)| \geq \varepsilon \} \leq \frac{D(X)}{\varepsilon^2}

这一不等式称为切比雪夫(Chebyshev)不等式,它的等价形式是

P{XE(X)<ε}1D(X)ε2P\{ |X - E(X)| < \varepsilon \} \geq 1 - \frac{D(X)}{\varepsilon^2}

大数定律

定义:X1,X2,,Xn,X_1, X_2, \cdots, X_n, \cdots 是一个随机变量序列,aa 是一个常数。如果对于任意给定的正数 ε\varepsilon,有

limnP{Xna<ε}=1\lim_{n \to \infty} P\{ |X_n - a| < \varepsilon \} = 1

则称随机变量序列 X1,X2,,Xn,X_1, X_2, \cdots, X_n, \cdots 依概率收敛aa,记作 XnPaX_n \stackrel{P}{\to} a.

性质:XnPa,YnPbX_n \stackrel{P}\to a, Y_n \stackrel{P}\to b,函数 g(x,y)g(x, y) 在点 (a,b)(a, b) 连续,则 g(Xn,Yn)Pg(a,b)g(X_n, Y_n) \stackrel{P}\to g(a, b).

切比雪夫大数定律:X1,X2,,Xn,X_1, X_2, \cdots, X_n, \cdots相互独立的随机变量序列,具有数学期望 E(Xk)E(X_k) 及方差 D(Xk)D(X_k),且方差一致有上界,即存在正数 MM 使得 D(Xk)MD(X_k) \leq M. 则对任意给定的正数 ε\varepsilon,恒有

limnP{1nk=1nXk1nk=1nE(Xk)<ε}=1\lim_{n \to \infty} P\left\{ \left| \frac{1}{n} \sum_{k = 1}^n X_k - \frac{1}{n} \sum_{k = 1}^n E(X_k) \right| < \varepsilon \right\} = 1

nn 充分大时,1ni=1nXi\dfrac{1}{n} \sum_{i = 1}^n X_i 差不多不再是随机的了,取值接近于其数学期望的概率接近于 1.

伯努利大数定理:nAn_Ann 次独立重复试验中事件 AA 发生的次数,pp 是事件 AA 在每次试验中发生的概率,则对任给的正数 ε\varepsilon,有

limnP{nAnp<ε}=1\lim_{n \to \infty} P\left\{ \left| \frac{n_A}{n} - p \right| < \varepsilon \right\} = 1

等价形式:

limnP{nAnpε}=0\lim_{n \to \infty} P\left\{ \left| \frac{n_A}{n} - p \right| \geq \varepsilon \right\} = 0

伯努利定理表明,当 nn 充分大时,事件 AAnn 次独立重复试验中发生的频率 nAn\dfrac{n_A}{n} 依概率收敛于其发生的概率 pp.

辛钦(Khintchine)大数定理:X1,X2,,Xn,X_1, X_2, \cdots, X_n, \cdots相互独立的随机变量序列,服从同一分布,具有数学期望 E(Xk)=μ,k=1,2,E(X_k) = \mu, k = 1, 2, \cdots,则对任意给定的正数 ε\varepsilon,有

limnP{1nk=1nXkμ<ε}=1\lim_{n \to \infty} P\left\{ \left| \frac{1}{n} \sum_{k = 1}^n X_k - \mu \right| < \varepsilon \right\} = 1

中心极限定理