实高斯随机向量与复高斯随机向量

本文深入探讨了高斯随机向量,包括实高斯随机向量和复高斯随机向量的概念。实高斯随机向量在任何正交基上分布相同,线性变换后仍保持高斯分布特性。复高斯随机向量的统计特性由均值、协方差矩阵和伪协方差矩阵决定,循环对称的复高斯随机向量仅由协方差矩阵完全刻画。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

本文内容主要来自David Tse etc., Fundamentals of Wireless Communication,Cambridge University Press, 2004中的附录A。

1、实高斯随机向量

  若 w 1 , w 2 , … , w n w_1,w_2,\ldots,w_n w1,w2,,wn为i.i.d.标准高斯分布随机变量,则标准高斯分布向量 w = [ w 1 , w 2 , … , w n ] T {\bf w}=[w_1,w_2,\ldots,w_n]^{\rm T} w=[w1,w2,,wn]T的概率密度函数为
p ( w ) = 1 ( 2 π ) n exp ⁡ ( − ∣ ∣ w ∣ ∣ 2 2 ) , w ∈ R n . p({\bf w})=\frac{1}{(\sqrt{2\pi})^n}\exp(-\frac{||{\bf w}||^2}{2}),\qquad {\bf w}\in {\mathbb R}^n. p(w)=(2π )n1exp(2w2),wRn.这里 ∣ ∣ w ∣ ∣ = ∑ i = 1 n w i 2 ||{\bf w}||=\sqrt{\sum_{i=1}^{n}w_i^2} w=i=1nwi2 表示从原点到 w : = [ w 1 , w 2 , … , w n ] T {\bf w}:=[w_1,w_2,\ldots,w_n]^{\rm T} w:=[w1,w2,,wn]T的距离。不难发现,这里的概率密度函数只与幅度有关。
  进一步,由于正交变换 O {\bf O} O【即 O O T = O T O = I {\bf O}{\bf O}^{\rm T}={\bf O}^{\rm T}{\bf O}={\bf I} OOT=OTO=I】不改变信号幅度,可得如下结论

如果 w \bf w w为标准高斯随机向量,则 O W \bf OW OW也为标准高斯分布随机变量。

上述结果表明 w \bf w w在任何正态基上都具有相同分布。从几何上看, w \bf w w的分布不随旋转和反射改变,因此 w \bf w w不偏好任何特定方向。此外,由于矩阵 O \bf O O的行是标准正交的,标准高斯随机向量在正交方向上的投影是独立的。 ∣ ∣ w ∣ ∣ 2 ||{\bf w}||^2 w2等于 n n n个i.i.d.零均值高斯随机变量的平方和,称为具有 n n n个自由度的卡方分布,用 χ n 2 \chi_{n}^2 χn2表示。

  高斯分布随机向量可以看作标准高斯分布随机向量的线性变换,加上常数向量,即
x = A w + μ . \bf x=Aw+\bm \mu. x=Aw+μ.几个性质如下:
1)对于任意的 c ∈ R n {\bf c}\in {\mathbb R}^n cRn,都有随机变量
c T x ∼ N ( c T μ , c T A A T c ) . {\bf c}^{\rm T}{\bf x}\sim {\mathcal N}({\bf c}^{\rm T}{\bm \mu},{\bf c}^{\rm T}{\bf A}{\bf A}^{\rm T}{\bf c}). cTxN(cTμ,cTAATc).因此,高斯随机向量( x {\bf x} x)中元素的任意线性组合都为高斯随机变量( c T x {\bf c}^{\rm T}{\bf x} cTx)。更一般来说,高斯随机向量的任意线性组合也是高斯的。
2) 如果 A \bf A A可逆,则对于 x = A w + μ {\bf x}={\bf Aw+\bm \mu} x=Aw+μ,其中 w ∼ N ( 0 , 1 ) {\bf w}\sim \mathcal{N}({\bf 0},{\bf 1}) wN(0,1),有PDF为
p ( x ) = 1 ( 2 π ) n det ⁡ ( A A T ) exp ⁡ ( − 1 2 ( x − μ ) T ( A A T ) − 1 ( x − μ ) ) , x ∈ R n . p({\bf x})=\frac{1}{(\sqrt{2\pi})^n\sqrt{\det({\bf AA}^{\rm T})}}\exp {\large(} -\frac{1}{2}({\bf x}-{\bm \mu})^{\rm T}({\bf AA}^{\rm T})^{-1}({\bf x}-{\bm \mu}){\large)},\quad {\bf x}\in \mathbb{R}^n. p(x)=(2π )ndet(AAT) 1exp(21(xμ)T(AAT)1(xμ)),xRn.这里的 A A T {\bf AA}^{\rm T} AAT x \bf x x的协方差矩阵:
K : = E [ ( x − μ ) ( x − μ ) T ] = A A T {\bf K}:=\mathbb{E}{\large [} ({\bf x}-{\bm \mu}) ({\bf x}-{\bm \mu})^{\rm T} {\large ]}={\bf AA}^{\rm T} K:=E[(xμ)(xμ)T]=AAT对于可逆矩阵 A \bf A A,高斯随机向量可以完全由其均值向量 μ \bm \mu μ和协方差矩阵 K = A A T {\bf K}={\bf AA}^{\rm T} K=AAT来刻画,这里的 K \bf K K为对称且非负定矩阵。我们再来看几个推导:

  • 由于 O w \bf Ow Ow w \bf w w同分布,因此 A O w A\bf Ow AOw A w \bf Aw Aw同分布。
  • 如果 K \bf K K为对角阵,则高斯随机向量中的元素统计独立(不相关)。这样的随机向量也称为白高斯随机向量。当协方差矩阵 K \bf K K为单位阵时,则该高斯随机向量为标准高斯随机向量。
  • 现在考虑 A \bf A A不可逆的情况。此时, A w \bf Aw Aw将标准高斯随机向量 w \bf w w映射到维度小于 n n n的子空间上。这意味着 A w \bf Aw Aw中的一些成分可以表示为其它部分的线性组合。这样我们可以只关注 A w \bf Aw Aw中线性独立的元素,将其表示为具有更低维度的向量 x ~ \bf \tilde{x} x~,同时将 A w \bf Aw Aw中的其它元素表示为 x ~ \bf \tilde{x} x~的线性组合。这样,我们可以使得协方差矩阵 K \bf K K可逆。

通常来说,高斯随机向量可以由其均值向量 μ \bm \mu μ和协方差矩阵 K \bf K K来完全刻画,我们表示为 x ∼ N ( μ , K ) {\bf x}\sim \mathcal{N}({\bm \mu}, {\bf K}) xN(μ,K)

2、复高斯随机向量

  下面我们来考虑复随机向量 x = x R + x I {\bf x}={\bf x}_R+{\bf x}_I x=xR+xI。如果向量 [ x R , x I ] T [{\bf x}_R,{\bf x}_I]^{\rm T} [xR,xI]T为实高斯随机向量,则 x {\bf x} x为复高斯向量。分布完全由实向量 [ x R , x I ] T [{\bf x}_R, {\bf x}_I]^{\rm T} [xR,xI]T的均值和协方差矩阵来决定。可以证明,同样的信息包含在 x \bf x x的均值向量 μ \bm \mu μ、协方差矩阵 K \bf K K以及伪协方差矩阵 J \bf J J中,即
μ : = E [ x ] \bm \mu:={\mathbb E}[\bf x] μ:=E[x] K : = E [ ( x − μ ) ( x − μ ) ∗ ] {\bf K}:={\mathbb E} {\large [ } (\bf x-\bm \mu)(x-\mu)^{*} {\large]} K:=E[(xμ)(xμ)] J : = E [ ( x − μ ) ( x − μ ) T ] . \bf J:={\mathbb E} {\large [} (\bf x-\bm \mu)(x-\mu)^{\rm T} {\large]}. J:=E[(xμ)(xμ)T].注意通常来说,复随机向量 x \bf x x的协方差矩阵 K \bf K K是不能完全刻画其二阶统计特性的。事实上,由于 K \bf K K为Hermitian矩阵(复共轭对称矩阵),则它是由 n 2 n^2 n2个实参数来决定的。另外一方面, x \bf x x的全部二阶统计量是由 2 n × 2 n 2n \times 2n 2n×2n维的协方差矩阵 [ x R , x I ] T [{\bf x}_R,{\bf x}_I]^{\rm T} [xR,xI]T中的 n ( 2 n + 1 ) n(2n+1) n(2n+1)个实参数决定的。

循环对称特性:如果对于任意 θ \theta θ e j θ x e^{j\theta}\bf x ejθx都与 x \bf x x x \bf x x同分布,则称$\bf $具有循环对称性。

对于任意 θ \theta θ,循环对称复随机变量 x \bf x x都有
E [ x ] = E [ e j θ x ] = e j θ E [ x ] , {\mathbb E}[{\bf x}]={\mathbb E}[e^{j\theta}{\bf x}]=e^{j\theta}{\mathbb E}[{\bf x}], E[x]=E[ejθx]=ejθE[x],因此均值 μ = 0 \bm \mu=\bf 0 μ=0。进一步,对于任意 θ \theta θ,有
E [ x x T ] = E [ e j θ x ( e j θ x ) T ] = e j 2 θ E [ x x T ] , {\mathbb E}[{\bf xx}^{\rm T}]={\mathbb E}[e^{j\theta}{\bf x}(e^{j\theta}{\bf x})^{\rm T}]=e^{j2\theta}{\mathbb E}[{\bf xx}^{\rm T}], E[xxT]=E[ejθx(ejθx)T]=ej2θE[xxT],因此伪协方差矩阵 J \bf J J也为零。故协方差矩阵 K \bf K K可以完全刻画循环对称随机向量的一阶和二阶统计特性。并且如果随机向量为高斯的,则 K \bf K K可以刻画所有统计特性。协方差矩阵为 K \bf K K的循环对称高斯随机向量记为 C N ( 0 , K ) {\mathcal CN}(0,\bf K) CN(0,K)
  下面来看几个特例:

  • 复高斯随机变量 w = w R + j w I w=w_R+jw_I w=wR+jwI,这里的实部和虚部为.零均值高斯i.i.d实随机变量,则 w w w为循环对称,而其循环对称性其实就是实高斯随机向量 [ w R , w I ] T [w_R,w_I]^{\rm T} [wR,wI]T的旋转不变形。事实上,循环对称复高斯随机变量的实部和虚部一定是零均值i.i.d.的。其统计特性可以完全由方差 σ 2 : = E [ ∣ w ∣ 2 ] \sigma^2:={\mathbb E}[|w|^2] σ2:=E[w2]来决定,记为 w ∼ C N ( 0 , σ 2 ) w\sim {\mathcal CN}(0,\sigma^2) wCN(0,σ2)。注意如果不是循环对称的,则复高斯随机变量需要用五个参数来给定:实部均值、虚部均值、实部方差、虚部方差以及虚实部相关。 w w w的相位在 [ 0 , 2 π ] [0,2\pi] [0,2π]上均匀分布,并与幅度 r = ∣ ∣ w ∣ ∣ r=||w|| r=w独立,后者的PDF为
    P ( r ) = r σ 2 exp ⁡ { − r 2 2 σ 2 } , r ≥ 0 P(r)=\frac{r}{\sigma^2}\exp{\Large \{} \frac{-r^2}{2\sigma^2}{\Large \}},\quad r\ge 0 P(r)=σ2rexp{2σ2r2},r0满足瑞利分布。幅度的平方,即 r 2 = w R 2 + w I 2 r^2=w_R^2+w_I^2 r2=wR2+wI2,满足 χ x 2 \chi_x^2 χx2分布即指数分布。
  • n n n个i.i.d. C N ( 0 , 1 ) {\mathcal CN}(0,1) CN(0,1)随机变量构成标准循环对称高斯随机向量 w ∼ C N ( 0 , I ) {\bf w}\sim {\mathcal CN}(0,\bf{I}) wCN(0,I),其PDF为
    P ( w ) = 1 π n exp ⁡ ( − ∣ ∣ w ∣ ∣ 2 ) , w ∈ C n . P({\bf w})=\frac{1}{\pi^n}\exp (-||{\bf w}||^2),\quad {\bf w}\in {\mathbb C}^n. P(w)=πn1exp(w2),wCn.与实高斯随机向量 N ( 0 , I ) {\mathcal N}(0,{\bf I}) N(0,I)类似,我们有如下特性特性

U w \bf Uw Uw w \bf w w同分布,这里的 U \bf U U为酉矩阵。

与实数情况相似, w \bf w w的幅度满足 χ 2 n 2 \chi_{2n}^2 χ2n2分布。

  • 如果 w ∼ C N ( 0 , I ) {\bf w}\sim {\mathcal CN}(0,\bf I) wCN(0,I)并且 A \bf A A为复矩阵,则 x = A w \bf x=Aw x=Aw也是循环对称高斯的,且 x ∼ C N ( 0 , K ) {\bf x}\sim {\mathcal CN}(0,\bf K) xCN(0,K)。如果 A \bf A A可逆,则 x \bf x x的PDF为
    p ( x ) = 1 π n det ⁡ K exp ⁡ ( − x ∗ K − 1 x ) , x ∈ C n . p({\bf x})=\frac{1}{\pi^n \sqrt{\det {\bf K}}}\exp(-\bf x^*K^{-1}x),\quad x\in \mathbb C^n. p(x)=πndetK 1exp(xK1x),xCn.
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值