实高斯随机向量与复高斯随机向量
本文内容主要来自David Tse etc., Fundamentals of Wireless Communication,Cambridge University Press, 2004中的附录A。文章目录1、实高斯随机向量2、复高斯随机向量1、实高斯随机向量 若w1,w2,…,wnw_1,w_2,\ldots,w_nw1,w2,…,wn为i.i.d.标准高斯分布随机变量,则标准...
本文内容主要来自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(−2∣∣w∣∣2),w∈Rn.这里 ∣ ∣ 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 ∣∣w∣∣2等于 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 c∈Rn,都有随机变量
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}). cTx∼N(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}) w∼N(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−μ)),x∈Rn.这里的 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}) x∼N(μ,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[∣w∣2]来决定,记为 w ∼ C N ( 0 , σ 2 ) w\sim {\mathcal CN}(0,\sigma^2) w∼CN(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σ2−r2},r≥0满足瑞利分布。幅度的平方,即 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}) w∼CN(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(−∣∣w∣∣2),w∈Cn.与实高斯随机向量 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) w∼CN(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) x∼CN(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)=πndetK1exp(−x∗K−1x),x∈Cn.
更多推荐



所有评论(0)