M05 · 第 4 章 · 第二编 随机变量

期望、方差、协方差与条件分解

在矩存在的条件下定义期望、方差、协方差与相关系数,再把条件期望视为随机变量,推导塔式法则、全期望和全方差公式。

报告页面错误
预备知识随机变量与分布条件概率与独立性

本章目标

  1. 在离散、绝对连续和一般分布下写出期望,并先检查绝对可积条件。
  2. 使用 LOTUS、指示变量与期望线性性计算函数和计数的平均值。
  3. 推导方差、协方差、相关系数与线性组合方差,并保持量纲解释正确。
  4. 用协方差矩阵的半正定性检查多变量二阶计算,并区分零协方差与独立。
  5. 把条件期望识别为随机变量,使用塔式法则、全期望和全方差完成分层计算。
页面阅读位置0% · 仅保存在此浏览器
章节未开始
本册完成进度0/6 章 · 0%
本页目录

分布摘要保留什么信息

随机变量的分布决定所有 Borel 集的概率,信息完整却未必便于比较。期望把取值按概率权重压缩成一个中心,方差记录中心周围的平方波动,协方差记录两个中心化变量乘积的平均符号与尺度。条件期望进一步保留分层信息:它先在已知条件下求平均,再让条件本身随试验结果变化。

有限个矩不能唯一确定一般分布。两个分布可能具有相同均值和方差,却拥有不同的偏斜、尾部或点质量;零协方差也只排除二阶线性共同变化。使用摘要之前要先问三个问题:相关矩是否存在,单位是否一致,当前结论是否真正只依赖这些矩。

期望的存在条件与 LOTUS

有限期望

XX 的分布为 PX\mathbb P_X。若

E[X]=RxdPX(x)<,\mathbb E[|X|] =\int_{\mathbb R}|x|\,\mathrm d\mathbb P_X(x) <\infty,

则称 XX 可积,并定义

E[X]=RxdPX(x).\mathbb E[X] =\int_{\mathbb R}x\,\mathrm d\mathbb P_X(x).

离散情形为 E[X]=xxpX(x)\mathbb E[X]=\sum_x x\,p_X(x);绝对连续情形为 E[X]=xfX(x)dx\mathbb E[X]=\int_{-\infty}^{\infty}x f_X(x)\,\mathrm dx

绝对可积条件排除“正无穷与负无穷相减”。标准 Cauchy 分布关于原点对称,但 E[X]=\mathbb E[|X|]=\infty,所以普通意义下的期望不存在。对称截断积分的主值为零,不能替代有限期望。

Y=g(X)Y=g(X),函数期望由 LOTUS 直接计算:

E[g(X)]=Rg(x)dPX(x).\mathbb E[g(X)] =\int_{\mathbb R}g(x)\,\mathrm d\mathbb P_X(x).

离散时对 g(x)pX(x)g(x)p_X(x) 求和,绝对连续时对 g(x)fX(x)g(x)f_X(x) 积分。该公式不要求先求出 YY 的完整分布,但仍要检查 E[g(X)]<\mathbb E[|g(X)|]<\infty。变量变换中的 Jacobian 已经包含在 YY 的密度里;使用 XX 的分布直接算 LOTUS 时不再重复加入。

期望不必属于支持。公平六面骰子的期望为 3.53.5,一次结果却只能是整数。期望是概率加权中心;把样本平均与该中心连接起来,还需要后续的大数定律及其条件。

指示变量把计数拆成加法

事件 AA 的指示变量记为

1A(ω)={1,ωA,0,ωA.\mathbf 1_A(\omega)= \begin{cases} 1,&\omega\in A,\\ 0,&\omega\notin A. \end{cases}

它满足

E[1A]=P(A).\mathbb E[\mathbf 1_A]=\mathbb P(A).

NN 统计若干事件中发生的个数,则 N=i=1n1AiN=\sum_{i=1}^n\mathbf 1_{A_i}。这个分解把计数期望转化为事件概率之和;事件之间是否独立不影响期望线性性。

有限期望的线性性

X,YX,Y 可积,a,ba,b 为常数,则

E[aX+bY]=aE[X]+bE[Y].\mathbb E[aX+bY] =a\mathbb E[X]+b\mathbb E[Y].
证明

可积性给出 E[aX+bY]aEX+bEY<\mathbb E[|aX+bY|]\le |a|\mathbb E|X|+|b|\mathbb E|Y|<\infty。 相对于共同概率空间积分,并使用积分线性性,便得到结论。证明没有使用 X,YX,Y 独立。

线性性对有限和可反复使用。无穷级数与期望交换需要单调收敛、控制收敛或绝对可积等额外依据;有限和公式本身不提供交换极限的许可。

非负整数值随机变量还有一条实用的尾和公式:

E[X]=k=1P(Xk),\mathbb E[X] =\sum_{k=1}^{\infty}\mathbb P(X\ge k),

允许等式两侧同时为正无穷。逐点恒等式 X=k11{Xk}X=\sum_{k\ge1}\mathbf 1_{\{X\ge k\}} 把取值 mm 写成前 mm 个指示量之和;对非负部分和使用单调收敛即可交换求和与期望。几何等待时间、排队长度和寿命的期望常由尾概率求得,这条公式也说明“尾部很小”必须与无穷项累加一起判断。

B1,,BnB_1,\ldots,B_n 是成功指示变量且每个成功概率为 pp,则 N=iBiN=\sum_iB_i 满足 E[N]=np\mathbb E[N]=np。独立性决定 NN 是否服从二项分布以及方差是否等于 np(1p)np(1-p),但均值 npnp 只依赖各项边缘概率。

一张分布表的前两阶矩

由离散分布计算 E=1.7、二阶矩 3.7 与方差 0.81

设一天的需求量 XX0,1,2,30,1,2,3,概率分别为 0.1,0.3,0.4,0.20.1,0.3,0.4,0.2。概率非负且总和为一。直接求和得到

E[X]=0×0.1+1×0.3+2×0.4+3×0.2=1.7.\mathbb E[X] =0\times0.1+1\times0.3+2\times0.4+3\times0.2 =1.7.

使用 LOTUS 取 g(x)=x2g(x)=x^2

E[X2]=02×0.1+12×0.3+22×0.4+32×0.2=3.7.\mathbb E[X^2] =0^2\times0.1+1^2\times0.3+2^2\times0.4+3^2\times0.2 =3.7.

因此

Var(X)=E[X2](E[X])2=3.71.72=0.81.\operatorname{Var}(X) =\mathbb E[X^2]-(\mathbb E[X])^2 =3.7-1.7^2 =0.81.

回到定义式核对:

x(x1.7)2pX(x)=0.289+0.147+0.036+0.338=0.81.\sum_x(x-1.7)^2p_X(x) =0.289+0.147+0.036+0.338 =0.81.

两条路径一致,标准差为 σX=0.9\sigma_X=0.9

方差、标准差与高阶矩

矩、方差与标准差

E[Xk]<\mathbb E[|X|^k]<\infty,则第 kk 阶原点矩记为

mk=E[Xk].m_k=\mathbb E[X^k].

若二阶矩有限,令 μX=E[X]\mu_X=\mathbb E[X],方差与标准差为

Var(X)=E[(XμX)2],σX=Var(X).\operatorname{Var}(X) =\mathbb E[(X-\mu_X)^2], \qquad \sigma_X=\sqrt{\operatorname{Var}(X)}.

展开平方并使用期望线性性得到

Var(X)=E[X2]μX2.\operatorname{Var}(X) =\mathbb E[X^2]-\mu_X^2.

定义式立即说明方差非负;若方差为零,则 (XμX)2=0(X-\mu_X)^2=0 几乎处处,因此 X=μXX=\mu_X 几乎处处。平移和缩放满足

Var(aX+b)=a2Var(X).\operatorname{Var}(aX+b)=a^2\operatorname{Var}(X).

XX 以米为单位,期望和标准差以米为单位,方差以平方米为单位。方差计算式中的两个大数 E[X2]\mathbb E[X^2]μX2\mu_X^2 可能非常接近,浮点实现宜使用中心化或稳定的在线更新,不能把灾难性消减误认为真实的负方差。

平方函数是凸函数。只要两侧定义,Jensen 不等式给出

(E[X])2E[X2],(\mathbb E[X])^2\le\mathbb E[X^2],

这也是方差非负的另一条证明。更高阶矩可描述偏斜和尾部,但本章只在公式真正需要时引入,不用一串矩代替完整分布。

方差还能把“偏离均值的概率”转成通用上界。若 Var(X)<\operatorname{Var}(X)<\infty,则对每个 t>0t>0

P(XμXt)Var(X)t2.\mathbb P(|X-\mu_X|\ge t) \le \frac{\operatorname{Var}(X)}{t^2}.

这是 Chebyshev 不等式。证明把非负变量 (XμX)2(X-\mu_X)^2 与阈值 t2t^2 比较:

Var(X)=E[(XμX)2]t2P(XμXt).\operatorname{Var}(X) =\mathbb E[(X-\mu_X)^2] \ge t^2\mathbb P(|X-\mu_X|\ge t).

该界只使用前两阶矩,因而适用于许多不同形状的分布;代价是它通常不如已知具体分布时的精确尾概率尖锐。下一章会把同一不等式用于样本平均,证明有限方差版本的弱大数定律。

线性成本怎样改变中心与尺度

需求量映射为成本后的均值 285 与方差 2025

沿用上一例的需求量,设每日成本为

C=50X+200 元.C=50X+200\ \text{元}.

期望成本为

E[C]=50E[X]+200=50×1.7+200=285 元.\mathbb E[C] =50\mathbb E[X]+200 =50\times1.7+200 =285\ \text{元}.

固定成本改变中心,不贡献随机波动。方差为

Var(C)=502Var(X)=2500×0.81=2025 2.\operatorname{Var}(C) =50^2\operatorname{Var}(X) =2500\times0.81 =2025\ \text{元}^2.

标准差为

σC=50σX=45 元.\sigma_C=50\sigma_X=45\ \text{元}.

直接把成本的四个取值 200,250,300,350200,250,300,350 代入概率表,也会得到同一均值和方差。期望与标准差保留原量的单位,方差保留单位平方,三种解释不能互换。

协方差组织线性组合的波动

协方差

X,YX,Y 都有有限二阶矩,则

Cov(X,Y)=E[(XμX)(YμY)].\operatorname{Cov}(X,Y) =\mathbb E[(X-\mu_X)(Y-\mu_Y)].

展开后得到

Cov(X,Y)=E[XY]E[X]E[Y].\operatorname{Cov}(X,Y) =\mathbb E[XY]-\mathbb E[X]\mathbb E[Y].

有限二阶矩由 Cauchy–Schwarz 不等式保证乘积 XYXY 可积。协方差对两个变量对称且双线性,常数平移不改变它。对任意 a,ba,b

Var(aX+bY)=a2Var(X)+b2Var(Y)+2abCov(X,Y).\operatorname{Var}(aX+bY) =a^2\operatorname{Var}(X) +b^2\operatorname{Var}(Y) +2ab\operatorname{Cov}(X,Y).

正协方差的严格含义是中心化乘积的平均值为正;多数样本点是否落在同号象限,不能由协方差符号单独判定。协方差大小随单位缩放:把米换成厘米会把相应协方差乘一百,所以不同变量对之间不能脱离单位比较绝对数值。

对有限个随机变量,逐项展开得到

Var ⁣(i=1nXi)=i=1nVar(Xi)+21i<jnCov(Xi,Xj).\operatorname{Var}\!\left(\sum_{i=1}^nX_i\right) =\sum_{i=1}^n\operatorname{Var}(X_i) +2\sum_{1\le i<j\le n}\operatorname{Cov}(X_i,X_j).

独立且二阶矩有限时,联合分布因子化给出 E[XiXj]=E[Xi]E[Xj]\mathbb E[X_iX_j]=\mathbb E[X_i]\mathbb E[X_j],所有交叉协方差消失。反方向一般不成立,因而“方差相加”必须由独立性、两两零协方差或其他明确结构支撑,不能只凭变量名称不同。

交叉项把线性组合方差改为 37

已知

Var(X)=4,Var(Y)=9,Cov(X,Y)=3.\operatorname{Var}(X)=4,\qquad \operatorname{Var}(Y)=9,\qquad \operatorname{Cov}(X,Y)=-3.

Z=2XY+5Z=2X-Y+5。常数五不影响方差,故

Var(Z)=22Var(X)+(1)2Var(Y)+22(1)Cov(X,Y)=16+9+12=37.\begin{aligned} \operatorname{Var}(Z) &=2^2\operatorname{Var}(X) +(-1)^2\operatorname{Var}(Y) +2\cdot2\cdot(-1)\operatorname{Cov}(X,Y)\\ &=16+9+12\\ &=37. \end{aligned}

交叉项系数是 2ab=42ab=-4,再乘负协方差得到正十二。若误把线性组合当成独立和而删掉协方差项,答案会错成二十五。

相关系数有明确的定义域

0<σX<0<\sigma_X<\infty0<σY<0<\sigma_Y<\infty 时,相关系数定义为

ρX,Y=Cov(X,Y)σXσY.\rho_{X,Y} =\frac{\operatorname{Cov}(X,Y)}{\sigma_X\sigma_Y}.

Cauchy–Schwarz 不等式给出 1ρX,Y1-1\le\rho_{X,Y}\le1。相关系数消除了正单位缩放,概括线性关联强度;ρ=1|\rho|=1 表示两个中心化变量之间存在几乎处处的精确线性关系。

等号条件来自 Cauchy–Schwarz 的等号情形:存在非零常数 aa,使 YμY=a(XμX)Y-\mu_Y=a(X-\mu_X) 几乎处处。a>0a>0 时相关系数为一,a<0a<0 时为负一。若散点近似落在直线上但仍有非零残差,相关系数的绝对值只会接近一,不会精确等于一。

零方差时相关系数未定义。若 C5C\equiv5,则 Var(C)=0\operatorname{Var}(C)=0Cov(C,X)=0\operatorname{Cov}(C,X)=0,但

ρC,X=00σX\rho_{C,X} =\frac{0}{0\cdot\sigma_X}

没有数值。不能因为协方差等于零便填写相关系数零。相关系数也不给出因果方向;共同原因、选择机制和时间趋势都可能产生非零相关。

零协方差仍可保留确定依赖

平方关系产生零协方差但不产生独立

XX{1,0,1}\{-1,0,1\} 上均匀,定义 Y=X2Y=X^2。由于 YYXX 唯一确定,二者不独立。例如

P(Y=0X=0)=1P(Y=0)=13.\mathbb P(Y=0\mid X=0)=1 \ne \mathbb P(Y=0)=\frac13.

利用对称性,

E[X]=0,E[Y]=E[X2]=23.\mathbb E[X]=0,\qquad \mathbb E[Y]=\mathbb E[X^2]=\frac23.

又因 XY=X3XY=X^3

E[XY]=1+0+13=0,\mathbb E[XY] =\frac{-1+0+1}{3}=0,

所以

Cov(X,Y)=00×23=0.\operatorname{Cov}(X,Y) =0-0\times\frac23=0.

进一步核对 Var(X)=2/3\operatorname{Var}(X)=2/3Var(Y)=2/9\operatorname{Var}(Y)=2/9,两个方差都为正,相关系数确实等于零。零相关来自对称性抵消,平方依赖并未消失。

协方差矩阵检查多变量结果

随机向量的协方差矩阵

对具有有限二阶矩的随机向量

X=(X1,,Xd)T,μ=E[X],\mathbf X=(X_1,\ldots,X_d)^{\mathsf T}, \qquad \boldsymbol\mu=\mathbb E[\mathbf X],

协方差矩阵定义为

Σ=Cov(X)=E[(Xμ)(Xμ)T].\Sigma =\operatorname{Cov}(\mathbf X) =\mathbb E[(\mathbf X-\boldsymbol\mu) (\mathbf X-\boldsymbol\mu)^{\mathsf T}].

对角元是各分量方差,非对角元是两两协方差,因此 Σ\Sigma 对称。对任意确定向量 a\mathbf a

aTΣa=Var(aTX)0.\mathbf a^{\mathsf T}\Sigma\mathbf a =\operatorname{Var}(\mathbf a^{\mathsf T}\mathbf X) \ge0.

这证明协方差矩阵半正定。手算或程序若产生显著负特征值,至少有一处定义、数据处理或数值计算不一致。微小负值有时来自舍入,但必须结合容差和构造方式判断。

线性变换 Y=AX+b\mathbf Y=A\mathbf X+\mathbf b 满足

Cov(Y)=AΣAT.\operatorname{Cov}(\mathbf Y)=A\Sigma A^{\mathsf T}.

公式不要求正态分布。正态假设只在更强的结论中出现,例如均值与协方差共同确定多元正态分布。

二维协方差矩阵的半正定性还给出

Cov(X,Y)2Var(X)Var(Y).\operatorname{Cov}(X,Y)^2 \le \operatorname{Var}(X)\operatorname{Var}(Y).

这正是相关系数绝对值不超过一的矩阵版本:二阶主子式 Var(X)Var(Y)Cov(X,Y)2\operatorname{Var}(X)\operatorname{Var}(Y)-\operatorname{Cov}(X,Y)^2 不能为负。给定两个方差时,任意超出该范围的协方差都不可能来自合法联合分布。

条件期望本身是随机变量

关于随机变量的条件期望

XX 可积,即 E[X]<\mathbb E[|X|]<\infty。条件期望 E[XY]\mathbb E[X\mid Y] 是满足以下三项条件的随机变量 ZZZZ 关于 σ(Y)\sigma(Y) 可测,E[Z]<\mathbb E[|Z|]<\infty,并且对每个 Aσ(Y)A\in\sigma(Y) 都有

E[1AZ]=E[1AX].\mathbb E[\mathbf 1_A Z] =\mathbb E[\mathbf 1_A X].

这样的 ZZ 存在,并且在几乎处处相等的意义下唯一。在离散情形,若 pY(y)>0p_Y(y)>0,可取 Z=m(Y)Z=m(Y),其中

m(y)=E[XY=y]=xxpXY(xy).m(y)=\mathbb E[X\mid Y=y] =\sum_x x\,p_{X\mid Y}(x\mid y).

连续且条件密度存在时,把求和换成

m(y)=xfXY(xy)dx.m(y)=\int_{-\infty}^{\infty} x f_{X\mid Y}(x\mid y)\,\mathrm dx.

固定某个 yy 后,E[XY=y]\mathbb E[X\mid Y=y] 是数;试验前让 YY 保持随机,m(Y)m(Y) 便随条件层变化。它只能使用 YY 提供的信息,因此是 YY 的可测函数。事件条件 E[XB]\mathbb E[X\mid B]P(B)>0\mathbb P(B)>0 时也是一个数,不能与随机变量 E[XY]\mathbb E[X\mid Y] 混写。

定义中的积分等式表示:在任何只由 YY 决定的事件上,条件期望与原变量具有相同的加权平均。离散求和与条件密度积分都是这个平均保持性质的具体实现。

XX 有有限二阶矩,则对每个平方可积函数 h(Y)h(Y) 还有

E ⁣[(XE[XY])h(Y)]=0.\mathbb E\!\left[(X-\mathbb E[X\mid Y])h(Y)\right]=0.

证明只需先对 YY 条件化:h(Y)h(Y) 在已知 YY 后是常量,而 E[XE[XY]Y]=0\mathbb E[X-\mathbb E[X\mid Y]\mid Y]=0。因此条件期望留下的残差与任何只依赖 YY 的平方可积函数正交;这个性质解释了它为何是利用 YY 信息预测 XX 时的均方误差最优函数。

更具体地,记 M=E[XY]M=\mathbb E[X\mid Y]。对任意平方可积的 h(Y)h(Y),正交性消去交叉项并给出

E[(Xh(Y))2]=E[(XM)2]+E[(Mh(Y))2].\mathbb E[(X-h(Y))^2] =\mathbb E[(X-M)^2] +\mathbb E[(M-h(Y))^2].

因此 MM 在所有平方可积的 YY 可测预测量中最小化均方误差,并且极小解在几乎处处意义下唯一。这个最优性依赖平方损失和有限二阶矩;换成绝对损失等其他目标时,最优预测量可以不同。

条件期望保留线性性。若 X,ZX,Z 可积,a,ba,b 为常数,则

E[aX+bZY]=aE[XY]+bE[ZY].\mathbb E[aX+bZ\mid Y] =a\mathbb E[X\mid Y]+b\mathbb E[Z\mid Y].

G(Y)G(Y) 是有界的 YY 的函数,还可把已知因子提出:

E[G(Y)XY]=G(Y)E[XY].\mathbb E[G(Y)X\mid Y] =G(Y)\mathbb E[X\mid Y].

这些等式都在“几乎处处”意义下成立,因为条件期望本身只由平均保持性质确定到零概率集合上的差异。计算时不同版本在零概率层取值可能不同,却给出相同的积分和概率结论。

塔式法则合并分层平均

塔式法则与全期望公式

XX 可积,则

E[E[XY]]=E[X].\mathbb E[\mathbb E[X\mid Y]] =\mathbb E[X].

离散 YY 时,它写成

E[X]=yE[XY=y]pY(y).\mathbb E[X] =\sum_y \mathbb E[X\mid Y=y]\,p_Y(y).
证明

离散情形直接展开:

yE[XY=y]pY(y)=yxxpXY(xy)pY(y)=yxxpX,Y(x,y)=xxpX(x)=E[X].\begin{aligned} \sum_y\mathbb E[X\mid Y=y]p_Y(y) &=\sum_y\sum_x x\,p_{X\mid Y}(x\mid y)p_Y(y)\\ &=\sum_y\sum_x x\,p_{X,Y}(x,y)\\ &=\sum_x x\,p_X(x)\\ &=\mathbb E[X]. \end{aligned}

一般情形在条件期望的平均保持等式中取 A=ΩA=\Omega,便得到同一结论。

塔式法则先在每一层求条件均值,再按层概率加权。它与全概率公式结构相同:全概率分解事件概率,全期望分解数值平均。若分层变量进一步嵌套,还可逐层使用 E[E[XY,Z]Y]=E[XY]\mathbb E[\mathbb E[X\mid Y,Z]\mid Y]=\mathbb E[X\mid Y]

对有限或可数划分 B1,B2,B_1,B_2,\ldots,只要每个使用到的条件事件满足 P(Bi)>0\mathbb P(B_i)>0,全期望也写成

E[X]=iE[XBi]P(Bi).\mathbb E[X] =\sum_i\mathbb E[X\mid B_i]\mathbb P(B_i).

分层必须互斥并覆盖整个样本空间;重叠分组会重复计权,遗漏分组会丢失概率质量。零概率层不使用普通条件概率比值,应从划分中删除或改用更一般的条件分布表述。

全方差分解层内与层间波动

全方差公式

XX 有有限二阶矩,则

Var(X)=E[Var(XY)]+Var(E[XY]).\operatorname{Var}(X) = \mathbb E[\operatorname{Var}(X\mid Y)] + \operatorname{Var}(\mathbb E[X\mid Y]).
证明

M=E[XY]M=\mathbb E[X\mid Y]μ=E[X]\mu=\mathbb E[X]。分解

Xμ=(XM)+(Mμ).X-\mu=(X-M)+(M-\mu).

平方取期望得到两个平方项和一个交叉项。对交叉项先对 YY 条件化:

E[(XM)(Mμ)Y]=(Mμ)E[XMY]=0.\mathbb E[(X-M)(M-\mu)\mid Y] =(M-\mu)\mathbb E[X-M\mid Y] =0.

第一个平方项的期望是 E[Var(XY)]\mathbb E[\operatorname{Var}(X\mid Y)];第二个是 Var(M)\operatorname{Var}(M)。相加即得结论。

第一项是各层内部的平均波动,第二项是各层均值之间的波动。删掉任一项都会低估总方差,除非对应项确实为零。

XX 在每一层内都是常数,第一项为零,总波动全部来自层均值差异;若所有层均值相同,第二项为零,但层内仍可有随机波动。这个分解只陈述联合分布下的二阶恒等式。分层变量可能是标签、共同结果或选择机制,因果解释必须另由研究设计和机制假设提供。

分层总体的全方差复算

层内方差 6 与层间方差 24 合成总方差 30

设分层变量 YY 取零、一,概率分别为 0.6,0.40.6,0.4。条件矩为

Y=0Y=1E[XY]1020Var(XY)49\begin{array}{c|cc} &Y=0&Y=1\\\hline \mathbb E[X\mid Y]&10&20\\ \operatorname{Var}(X\mid Y)&4&9 \end{array}

全期望给出

E[X]=0.6×10+0.4×20=14.\mathbb E[X] =0.6\times10+0.4\times20 =14.

平均层内方差为

E[Var(XY)]=0.6×4+0.4×9=6.\mathbb E[\operatorname{Var}(X\mid Y)] =0.6\times4+0.4\times9 =6.

条件均值的方差为

0.6(1014)2+0.4(2014)2=9.6+14.4=24.0.6(10-14)^2+0.4(20-14)^2 =9.6+14.4 =24.

所以

Var(X)=6+24=30.\operatorname{Var}(X)=6+24=30.

再用二阶矩独立复核。条件二阶矩等于“条件方差加条件均值平方”,故

E[X2]=0.6(4+102)+0.4(9+202)=226.\mathbb E[X^2] =0.6(4+10^2)+0.4(9+20^2) =226.

于是 E[X2]E[X]2=226142=30\mathbb E[X^2]-\mathbb E[X]^2=226-14^2=30,两种计算完全一致。

矩摘要的三个边界

期望就是最可能出现的值

期望是概率加权中心,众数是概率质量或密度峰值所在位置。需求量例题的期望为 1.71.7,支持却只有整数。报告中心时要根据问题选择均值、中位数或众数。

相关为零便可按独立处理

零相关只给出协方差为零。平方关系例题保留确定依赖,联合概率不能因子化。只有在额外结构下,例如联合正态,零协方差才推出独立。

相同均值和方差仍对应不同分布

XX 以相同概率取 1,1-1,1;令 YY 以概率 1/8,3/4,1/81/8,3/4,1/82,0,2-2,0,2。两者均值都是零、方差都是一,但 P(X=0)=0\mathbb P(X=0)=0P(Y=0)=3/4\mathbb P(Y=0)=3/4。前两阶矩无法区分这两个分布。

重尾分布还会让某些矩不存在。软件对有限样本总能输出一个算术平均和平方离差,这些数字不保证对应总体具有有限期望或方差。先检查模型条件,再解释矩。

矩还会放大尾部。均值对极端值的权重与取值成正比,方差对偏离的权重按平方增长。若问题关心极端损失、分位数或超过阈值的概率,只报告均值和标准差通常不够;应同时回到 CDF、尾概率或条件分布。矩是为特定运算服务的摘要,不承担完整重建分布的职责。

练习:矩、条件化与二阶结构

练习 1:计算非对称离散分布的矩

随机变量 XX 以概率 1/4,1/2,1/41/4,1/2,1/4 取值 1,0,2-1,0,2。求期望、二阶矩、方差和标准差。

查看提示
先算 E[X] 与 E[X2]E[X^{2}],再使用 Var(X)=E[X2]E[X]2\operatorname{Var}(X)=E[X^{2}]-E[X]^{2} 复核定义式。
查看解答
E[X]=14+0+24=14,\mathbb E[X] =-\frac14+0+\frac{2}{4} =\frac14,
E[X2]=14+0+44=54.\mathbb E[X^2] =\frac14+0+\frac{4}{4} =\frac54.

因此

Var(X)=54(14)2=1916,σX=194.\operatorname{Var}(X) =\frac54-\left(\frac14\right)^2 =\frac{19}{16}, \qquad \sigma_X=\frac{\sqrt{19}}4.

三个概率和为一,方差为正,结果符合基本检查。

练习 2:不借助独立性求失效数期望

二十个元件中,每个元件在一天内失效的概率都是 0.030.03,但元件失效可能相关。求当天失效元件总数 NN 的期望。

查看提示
把总数写成二十个失效指示变量之和;题目没有要求计算方差。
查看解答

IiI_i 表示第 ii 个元件失效,则

N=i=120Ii.N=\sum_{i=1}^{20}I_i.

期望线性性不要求独立,所以

E[N]=i=120E[Ii]=i=120P(Ii=1)=20×0.03=0.6.\mathbb E[N] =\sum_{i=1}^{20}\mathbb E[I_i] =\sum_{i=1}^{20}\mathbb P(I_i=1) =20\times0.03 =0.6.

相关性会影响方差和联合失效概率,但不改变该期望。

练习 3:线性变换后的协方差

已知 Cov(X,Y)=2\operatorname{Cov}(X,Y)=2。求

Cov(3X1,42Y).\operatorname{Cov}(3X-1,\,4-2Y).

XX 以米为单位、YY 以秒为单位,写出结果单位。

查看提示
常数平移不影响协方差,缩放系数从两个变量分别提出。
查看解答

利用双线性与平移不变性,

Cov(3X1,42Y)=3(2)Cov(X,Y)=12.\operatorname{Cov}(3X-1,4-2Y) =3(-2)\operatorname{Cov}(X,Y) =-12.

结果单位为米乘秒。负号来自第二个变量乘以负二;常数 1-144 不贡献协方差。

练习 4:用两项公式计算分层波动

P(Y=0)=P(Y=1)=1/2\mathbb P(Y=0)=\mathbb P(Y=1)=1/2,并且

E[XY=0]=2,E[XY=1]=6,\mathbb E[X\mid Y=0]=2,\quad \mathbb E[X\mid Y=1]=6,
Var(XY=0)=1,Var(XY=1)=4.\operatorname{Var}(X\mid Y=0)=1,\quad \operatorname{Var}(X\mid Y=1)=4.

E[X]\mathbb E[X]Var(X)\operatorname{Var}(X)

查看提示
先算条件均值的加权平均,再分别计算平均条件方差和条件均值方差。
查看解答

全期望给出

E[X]=12×2+12×6=4.\mathbb E[X]=\frac12\times2+\frac12\times6=4.

平均条件方差为

12×1+12×4=52.\frac12\times1+\frac12\times4=\frac52.

条件均值以相同概率取二、六,其方差为

12(24)2+12(64)2=4.\frac12(2-4)^2+\frac12(6-4)^2=4.

因此

Var(X)=52+4=132=6.5.\operatorname{Var}(X)=\frac52+4=\frac{13}{2}=6.5.

两项均非负,总方差大于每个单独贡献。

练习 5:用二次型核对协方差矩阵

设随机向量 (X,Y)T(X,Y)^{\mathsf T} 的协方差矩阵为

Σ=[4229].\Sigma= \begin{bmatrix} 4&2\\ 2&9 \end{bmatrix}.

求相关系数与 Var(2XY)\operatorname{Var}(2X-Y),并核对 Σ\Sigma 半正定。

查看提示
先检查主对角元与行列式,再用线性组合方差公式计算 a=(2,1)a=(2,-1)
查看解答

标准差分别为二、三,所以

ρX,Y=22×3=13.\rho_{X,Y}=\frac{2}{2\times3}=\frac13.

线性组合方差为

Var(2XY)=4Var(X)+Var(Y)4Cov(X,Y)=16+98=17.\operatorname{Var}(2X-Y) =4\operatorname{Var}(X)+\operatorname{Var}(Y)-4\operatorname{Cov}(X,Y) =16+9-8 =17.

矩阵的首个主子式为四,行列式为 4×922=32>04\times9-2^2=32>0,故它正定,从而半正定。二次型在向量 (2,1)T(2,-1)^{\mathsf T} 上的值正是十七,与方差计算一致。

概念关系指向后续定理

  • 随机变量与分布 提供 PMF、PDF、联合与条件分布,决定期望中使用的概率权重。
  • 联合分布与条件分布 支撑乘积期望、协方差和条件矩。
  • 大数定律 说明样本平均在明确条件下趋近总体期望。
  • 中心极限定理 使用均值和方差标准化独立同分布和。
  • 线性回归 把协方差、线性组合和最小二乘联系起来。
  • 批归一化 在神经网络中使用批量矩;批量矩是随机统计量,固定总体参数则由总体分布定义。

从矩理论进入极限定理

课程 · 2013

MIT 6.041SC Probabilistic Systems Analysis and Applied Probability

John Tsitsiklis

课程包含讲授、习题、解答和考试,适合系统检验概率计算。

打开官方来源

MIT 6.041SC 的公开课程材料把期望、随机变量函数、条件化、条件期望和大数定律安排在连续课程链中,适合核对本章的存在条件、塔式法则与下一章的衔接。

书籍 · 2023

Introductory Statistics 2e

Barbara Illowsky, Susan Dean

提供大量分步例题和练习,适合核对分布、样本统计量、置信区间、假设检验和简单线性回归的基础计算与解释。

打开官方来源

OpenStax 开放教材第 3–4 章以概率表和离散分布计算期望与标准差,提供带实际单位的分步例题。它适合复算本章的需求量和成本例题;条件期望与全方差的严格推导仍以正文和概率课程材料为主。

两项资源侧重点不同:概率课程适合继续学习条件期望与分层恒等式,开放教材适合练习带单位的离散矩计算。均值 1.71.7、二阶矩 3.73.7、方差 0.810.81、成本方差 20252025、线性组合方差 3737 与全方差 3030 可作为重做本章例题时的结果自检。

下一章进入 大数定律

中心极限定理 。大数定律回答样本平均是否接近期望,中心极限定理描述标准化误差的分布形状;两者都要保留独立同分布、矩存在和标准化条件,不能由“样本量大”一句话代替。