知识点

协方差矩阵

定义

设 X=(X1,…,Xn)T\boldsymbol X=(X_1,\ldots,X_n)^{\mathsf T},每个分量的二阶矩都有限。定义均值向量与协方差矩阵

μ=E(X),Σ=E[(X−μ)(X−μ)T].\boldsymbol\mu=E(\boldsymbol X),\qquad \Sigma=E\bigl[(\boldsymbol X-\boldsymbol\mu)(\boldsymbol X-\boldsymbol\mu)^{\mathsf T}\bigr].

按分量写,μi=E(Xi)\mu_i=E(X_i),且

Σij=Cov⁡(Xi,Xj),i,j=1,…,n.\Sigma_{ij}=\operatorname{Cov}(X_i,X_j),\qquad i,j=1,\ldots,n.

所以 Σ\Sigma 是一个 n×nn\times n 矩阵:对角元 Σii=D(Xi)\Sigma_{ii}=D(X_i) 是各分量的方差,非对角元是两两协方差。协方差矩阵概括 nn 个变量的两两协方差,记法与多元正态分布 Nn(μ,Σ)N_n(\boldsymbol\mu,\Sigma) 的写法一致。

性质

对称与半正定

协方差矩阵满足

ΣT=Σ,\Sigma^{\mathsf T}=\Sigma,

且对任意 a∈Rn\boldsymbol a\in\mathbb R^n,

aTΣa≥0.\boldsymbol a^{\mathsf T}\Sigma\boldsymbol a\ge0.

也就是说,Σ\Sigma 是对称半正定矩阵。

证明

对称性来自

ΣT=E[(X−μ)(X−μ)T]T=E[(X−μ)(X−μ)T]=Σ.\Sigma^{\mathsf T} =E\bigl[(\boldsymbol X-\boldsymbol\mu)(\boldsymbol X-\boldsymbol\mu)^{\mathsf T}\bigr]^{\mathsf T} =E\bigl[(\boldsymbol X-\boldsymbol\mu)(\boldsymbol X-\boldsymbol\mu)^{\mathsf T}\bigr] =\Sigma.

对半正定性,注意 aT(X−μ)\boldsymbol a^{\mathsf T}(\boldsymbol X-\boldsymbol\mu) 是随机变量,于是

aTΣa=E[aT(X−μ)(X−μ)Ta]=E[(aT(X−μ))2]=D(aTX)≥0.\boldsymbol a^{\mathsf T}\Sigma\boldsymbol a =E\bigl[\boldsymbol a^{\mathsf T}(\boldsymbol X-\boldsymbol\mu)(\boldsymbol X-\boldsymbol\mu)^{\mathsf T}\boldsymbol a\bigr] =E\bigl[\bigl(\boldsymbol a^{\mathsf T}(\boldsymbol X-\boldsymbol\mu)\bigr)^2\bigr] =D(\boldsymbol a^{\mathsf T}\boldsymbol X)\ge0.

Σ\Sigma 只能保证半正定,不一定正定:只要存在非零向量 a\boldsymbol a 使 D(aTX)=0D(\boldsymbol a^{\mathsf T}\boldsymbol X)=0,aTΣa\boldsymbol a^{\mathsf T}\Sigma\boldsymbol a 就等于零,此时 Σ\Sigma 退化。

对角矩阵与不相关

Σ\Sigma 为对角矩阵,当且仅当各分量两两不相关。对角矩阵一般不能推出各分量相互独立;只有联合正态等额外条件成立时,不相关才等价于独立。

线性变换

若 Y=AX+b\boldsymbol Y=A\boldsymbol X+\boldsymbol b,其中 AA 是常数矩阵、b\boldsymbol b 是常数向量,则

E(Y)=Aμ+b,ΣY=AΣXAT.E(\boldsymbol Y)=A\boldsymbol\mu+\boldsymbol b,\qquad \Sigma_Y=A\Sigma_XA^{\mathsf T}.
证明

期望的线性性给出 E(Y)=AE(X)+b=Aμ+bE(\boldsymbol Y)=AE(\boldsymbol X)+\boldsymbol b=A\boldsymbol\mu+\boldsymbol b。又

Y−E(Y)=A(X−μ),\boldsymbol Y-E(\boldsymbol Y)=A(\boldsymbol X-\boldsymbol\mu),

所以

ΣY=E[A(X−μ)(X−μ)TAT]=AΣXAT.\Sigma_Y=E\bigl[A(\boldsymbol X-\boldsymbol\mu)(\boldsymbol X-\boldsymbol\mu)^{\mathsf T}A^{\mathsf T}\bigr] =A\Sigma_XA^{\mathsf T}.

取 A=aTA=\boldsymbol a^{\mathsf T}、b=0\boldsymbol b=0,得到线性组合的方差

D(aTX)=aTΣa.D(\boldsymbol a^{\mathsf T}\boldsymbol X)=\boldsymbol a^{\mathsf T}\Sigma\boldsymbol a.

这就是一维公式 D(∑iaiXi)=∑i∑jaiajCov⁡(Xi,Xj)D\left(\sum_i a_iX_i\right)=\sum_i\sum_j a_ia_j\operatorname{Cov}(X_i,X_j) 的矩阵写法。代入时检查 AA 与 Σ\Sigma 的维数以及转置的位置,AΣATA\Sigma A^{\mathsf T} 和 ATΣAA^{\mathsf T}\Sigma A 一般不相等。

例子

由协方差矩阵计算线性组合的方差

设

Σ=(4119),\Sigma=\begin{pmatrix}4&1\\1&9\end{pmatrix},

求 D(2X1−X2)D(2X_1-X_2)。取 a=(2,−1)T\boldsymbol a=(2,-1)^{\mathsf T},则

D(2X1−X2)=aTΣa=22×4+(−1)2×9+2×2×(−1)×1=16+9−4=21.D(2X_1-X_2)=\boldsymbol a^{\mathsf T}\Sigma\boldsymbol a =2^2\times4+(-1)^2\times9+2\times2\times(-1)\times1 =16+9-4=21.

对角项按系数的平方加权,交叉项要乘 2a1a22a_1a_2,这里的交叉项为负,所以结果小于 16+916+9。

交互探索

改变协方差和线性组合系数,计算相应的二次型。

协方差矩阵怎样决定线性组合的方差?

改变协方差与线性组合系数,计算对应的二次型。

正在准备演示…

题集

练习

  1. 设 Σ=(4119)\Sigma=\begin{pmatrix}4&1\\1&9\end{pmatrix},a=(1,2)T\boldsymbol a=(1,2)^{\mathsf T},则 aTΣa=\boldsymbol a^{\mathsf T}\Sigma\boldsymbol a= 。

    参考解答

    逐项展开:12×4+22×9+2×1×2×1=4+36+4=441^2\times4+2^2\times9+2\times1\times2\times1=4+36+4=44。

  2. 设 Σ\Sigma 是随机向量 X\boldsymbol X 的协方差矩阵。下列叙述正确的是

    选择答案
    参考解答

    二次型就是线性组合的方差,这一条对任意向量成立。Σ\Sigma 只保证半正定,可能退化,所以A错;对角只说明两两不相关,一般不推出独立,所以B错;线性变换的协方差矩阵是 AΣXATA\Sigma_XA^{\mathsf T},转置位置不对,所以D错。