欢迎光临散文网 会员登陆 & 注册

MIMO 信道容量的计算(一)

2022-12-23 23:21 作者:乐吧的数学  | 我要投稿

这篇文章,我们想分析一下 MIMO 信道的信道容量,有个前提:是基于非频率选择性衰落信道来分析的,即频率平稳衰落的信道。

(录制的视频:https://www.bilibili.com/video/BV1K24y1S77P/


Y%20%3D%20%5Csqrt%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20HS%20%2B%20W%20%20%5Cquad%20------%20%5Cquad%20(1)


信道的容量,就是分析互信息的最大值,用数学公式表示为:

C%20%3D%20%5Cunderset%7Bf(S)%7D%7B%20max%7D%20I(S%3BY)%20%20%20%20%20%5Cquad%20------%20%5Cquad%20(2)


其中 f(S) 是向量 S 的 联合概率分布函数。上面的公式的含义,就是在所有可能的概率分布 f(S) 中,找到使得互信息最大的那种概率分布。

下面来推导一下互信息 I(S;Y) :

I(S%3BY)%20%3D%20H(Y)%20-%20H(Y%7CS)%20%20%20%20%5Cquad%20------%20%5Cquad%20(3)


先来分析公式 (3) 的后半部分 H(Y|S)
(注: S 是离散的, Y 是连续的)

%5Cbegin%7Baligned%7D%0AH(Y%7CS)%0A%26%3D%20%5Csum_s%20p(S%3Ds)%20H(Y%7CS%3Ds)%20%20%5C%5C%0A%26%3D%20%5Csum_s%20p(S%3Ds)%20%5Cint_y%20%20p(y%7CS%3Ds)%20log%5Cfrac%7B1%7D%7Bp(y%7CS%3Ds)%7D%20dy%0A%5Cend%7Baligned%7D%20%20%20%20%5Cquad%20------%20%5Cquad%20(4)


利用公式 (1) 有

y%20%3D%20%5Csqrt%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20Hs%20%2B%20w%20%20%5Cquad%20------%20%5Cquad%20(5)


则公式 (4) 继续推导为:

%5Cbegin%7Baligned%7D%0AH(Y%7CS)%0A%26%3D%20%5Csum_s%20p(S%3Ds)%20%5Cint_w%20%20p(y%20%3D%20%5Csqrt%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20Hs%20%2B%20w)%20log%5Cfrac%7B1%7D%7Bp(y%20%3D%20%5Csqrt%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20Hs%20%2B%20w)%7D%20dw%20%20%5C%5C%0A%26%3D%20%5Csum_s%20p(S%3Ds)%20%5Cint_w%20%20p(w%20%3D%20y%20-%20%5Csqrt%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20Hs%20)%20log%5Cfrac%7B1%7D%7Bp(w%3Dy%20-%20%5Csqrt%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20Hs%20)%7D%0Adw%0A%5Cend%7Baligned%7D%20%5C%5C%20%20%20%5Cquad%20------%20%5Cquad%20(6)


因为噪声 w 的概率分布,与 s 的取值无关,所以,公式 (6) 后面的积分就与 s 的取值无关,则求和与求积分就可以独立开来:

%5Cbegin%7Baligned%7D%0AH(Y%7CS)%0A%26%3D%20%5B%5Csum_s%20p(S%3Ds)%5D%20%5B%20%5Cint_w%20%20p(w%20%3D%20y%20-%20%5Csqrt%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20Hs%20)%20log%5Cfrac%7B1%7D%7Bp(w%3Dy%20-%20%5Csqrt%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20Hs%20)%7D%0Adw%20%5D%20%5C%5C%0A%26%3D%20%5Cint_w%20%20p(w%20%3D%20y%20-%20%5Csqrt%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20Hs%20)%20log%5Cfrac%7B1%7D%7Bp(w%3Dy%20-%20%5Csqrt%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20Hs%20)%7D%0Adw%20%5C%5C%0A%26%3D%20%5Cint_w%20%20p(w%20)%20log%5Cfrac%7B1%7D%7Bp(w)%7D%20dw%20%20%5C%5C%0A%26%3D%20%20H(W)%20%20%5Cquad%20------%20%5Cquad%20(7)%0A%5Cend%7Baligned%7D


把公式 (7) 代入公式 (3) :

I(S%3BY)%20%3D%20H(Y)%20-%20H(W)%20%20%20%5Cquad%20------%20%5Cquad%20(8)


现在就相当于找 H(Y) 的最大值。

这里给出两个不加证明(还不知道咋证明 ^_^)的定理:

定理 1: 当给定 Y 的自相关矩阵 R_%7BYY%7D,那么 当 Y 满足 Zero-Mean Circularly Symmetric Complex Gaussian 分布时,Y的 differential entroy( 即连续熵) H(Y) 取最大值。

推论:根据公式 (1), Y 是 Zero-Mean Circularly Symmetric Complex Gaussian 分布, 则 S 也需要是 Zero-Mean Circularly Symmetric Complex Gaussian 分布。

定理 2:
(1)当 Y 是Zero-Mean Circularly Symmetric Complex Gaussian 分布时:

H(Y)%20%3D%20log_2(%7C%5Cpi%20e%20R_%7BYY%7D%7C)%20%20%20%5Cquad%20%5Ctext%7Bbps%2FHz%7D%20%20%20%5Cquad%20------%20%5Cquad%20(9)


(2)当 W 是Zero-Mean Circularly Symmetric Complex Gaussian 分布时,因为各个分量的白噪声之间相互独立,则:

H(W)%20%3D%20log_2(%7C%5Cpi%20e%20%20N_0%20I_%7BN_R%7D%7C)%20%20%20%5Cquad%20%5Ctext%7Bbps%2FHz%7D%20%20%20%5Cquad%20------%20%5Cquad%20(10)


把 (9) (10) 代入公式 (8) :

%5Cbegin%7Baligned%7D%0AI(S%3BY)%0A%26%3D%20H(Y)%20-%20H(W)%20%20%5C%5C%0A%26%3D%20%20log_2(%7C%5Cpi%20e%20R_%7BYY%7D%7C)%20%20-%20log_2(%7C%5Cpi%20e%20%20N_0%20I_%7BN_R%7D%7C)%5C%5C%0A%26%3D%20log_2(%7C%5Cfrac%7BR_%7BYY%7D%7D%7BN_0%7D%7C)%0A%5Cend%7Baligned%7D%20%20%5Cquad%20------%20%5Cquad%20(11)


注意: 本文中的 || 都是表示取方阵的行列式.

推导一下 R_%7BYY%7D%20


%5Cbegin%7Baligned%7D%0AR_%7BYY%7D%20%0A%26%3D%20E%5B(Y-EY)(Y-EY)%5EH%5D%20%20%5C%5C%0A%26%3DE(YY%5EH)%20%20%5C%5C%0A%26%3D%20E%5B(%5Csqrt%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20HS%20%2B%20W)%20(%5Csqrt%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20HS%20%2B%20W)%5EH%5D%20%5C%5C%0A%26%3D%20E%5B(%5Csqrt%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20HS%20%2B%20W)%20(%5Csqrt%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20S%5EH%20H%5EH%20%2B%20W%5EH)%5D%20%20%5C%5C%0A%26%3D%20E%5B%5Cfrac%7BE_s%7D%7BN_T%7D%20HSS%5EH%20H%5EH%20%2B%20%5Csqrt%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20HS%20W%5EH%20%2B%20W%20%5Csqrt%7B%5Cfrac%7BE_s%7D%7BN_T%7D%7D%20S%5EH%20H%5EH%20%2B%20WW%5EH%5D%20%20%5C%5C%0A%26%3D%20E(%5Cfrac%7BE_s%7D%7BN_T%7D%20HSS%5EH%20H%5EH)%20%2B%20E(%20WW%5EH)%20%5C%5C%0A%26%3D%20%5Cfrac%7BE_s%7D%7BN_T%7D%20HE(SS%5EH)%20H%5EH%20%2B%20N_0%20I_%7BN_R%7D%20%5C%5C%0A%26%3D%20%5Cfrac%7BE_s%7D%7BN_T%7D%20HR_%7BSS%7D%20H%5EH%20%2B%20N_0%20I_%7BN_R%7D%20%5Cquad%20------%20%5Cquad%20(12)%20%20%0A%5Cend%7Baligned%7D


把公式 (12) 代入公式 (11):

%5Cbegin%7Baligned%7D%0AI(S%3BY)%0A%26%3D%20H(Y)%20-%20H(W)%20%20%5C%5C%0A%26%3D%20log_2(%7C%5Cfrac%7BR_%7BYY%7D%7D%7BN_0%7D%7C)%20%20%5C%5C%0A%26%20%3D%20log_2(%7C%5Cfrac%7B%20%5Cfrac%7BE_s%7D%7BN_T%7D%20HR_%7BSS%7D%20H%5EH%20%2B%20N_0%20I_%7BN_R%7D%20%7D%7BN_0%7D%7C)%20%20%5C%5C%0A%26%3D%20%20log_2(%7C%5Cfrac%7BE_s%7D%7BN_T%20N_0%7D%20HR_%7BSS%7D%20H%5EH%20%2B%20%20I_%7BN_R%7D%20%7C)%0A%5Cend%7Baligned%7D%20%20%5Cquad%20------%20%5Cquad%20(13)


从公式 (13) 可以看出,这个互信息的最大取值位置,只与 R_%7BSS%7D 有关。
这篇文章里面都假定 H 不是随机变量,对于接收方 H 是已知的,是确定的。
把公式(13) 代入公式(2)

C%20%3D%20%5Cunderset%7Bf(S)%7D%7B%20max%7D%5B%20log_2(%7C%5Cfrac%7BE_s%7D%7BN_T%20N_0%7D%20HR_%7BSS%7D%20H%5EH%20%2B%20%20I_%7BN_R%7D%20%7C)%20%20%20%5D%20%20%20%20%20%20%5Cquad%20------%20%5Cquad%20(14)


在实际应用中 f(S)  这个分布,其实是要满足一个总能量一定的约束条件,我们假定能量都归一化了。总能量是 R_%7BSS%7D 的迹:

Tr(R_%7BSS%7D)%20%3D%20N_T


那么公式 (14) 就变成

C%20%3D%20%5Cunderset%7BTr(R_%7BSS%7D)%20%3D%20N_T%7D%7B%20max%7D%20log_2(%7C%5Cfrac%7BE_s%7D%7BN_T%20N_0%7D%20HR_%7BSS%7D%20H%5EH%20%2B%20%20I_%7BN_R%7D%20%7C)%20%20%20%20%20%20%5Cquad%20------%20%5Cquad%20(15)





参考书:
Introduction to Space-Time Wireless Communications, Arogyaswami Paulraj,Cambridge University Press 2003

MIMO 信道容量的计算(一)的评论 (共 条)

分享到微博请遵守国家法律