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

MIMO MMSE 的公式的推导

2023-06-11 08:05 作者:乐吧的数学  | 我要投稿

(录制的视频:待上传)

这篇文章推导一下 MIMO MMSE 检测算法的数学公式。MMSE:Minimum Mean Squared Error.

Y%3DHX%2BW%20%20%20%5Ctag1

其中 H 是信道系数矩阵,是 MxN 的,X 是发送信号向量,Nx1的列向量,Y 是接收到的信号,是 Mx1的列向量,W 是加性高斯白噪声,也是 Mx1 的列向量。


我们的目标是找一个矩阵 G,用 G 作用在 Y 上来估计发送的 X,即:

%5Chat%20X%20%3D%20GY%20%20%5Ctag%202

其中%5Chat%20X 是对发送向量 X 的估计,由于矩阵乘法是一个线性操作,因此这个算法也更准确地被称为 Linear MMSE 检测算法。


我们要找一个合适的 G,使得下式最小:

E_%7BX%2CW%7D%7C%7CX-%5Chat%20X%7C%7C%5E2%20%5Ctag%203

这个式子的含义是,X-%5Chat%20X  是估计的向量与原始的向量之间的误差向量,再对这个向量取其模长。


而这套系统中,我们是假定信道系数矩阵 H 已知,在这些条件下,那么只有 发送向量 X  和 噪声W 是未知的(由于 Y = HX +W ,因此 Y 也是未知的,但是可以用 X 和 W 计算出来),我们把他们看成随机变量,所以,我们从这两个随机变量的角度,看统计意义下的误差最小,即名称中 Mean 这个单词的含义。


也就是说,我们不是针对某个特定的发送向量,或者某个特定的白噪声,来计算 G,使得误差最小(误差的模长最小),而是看所有发送的向量的各种可能,所有可能的白噪声,考虑这些所有的可能后来使得误差最小。这是与 Zero Forcing 算法不太一样的地方(待详细叙述)。


把 (2) 代入 (3) ,同时,用数学语言来表示找最小,我们得到如下这个公式:

%5Chat%20G%20%3D%20%5Cunderset%7BG%7D%7Bargmin%7D%20%5C%7BE_%7BX%2CW%7D%7C%7CX-GY%7C%7C%5E2%5C%7D%20%5Ctag%204

接下来就是纯数学推导了:

%5Cbegin%7Baligned%7D%0A%0A%7C%7CX-GY%7C%7C%5E2%20%26%3D%20(X-GY)%5EH%20(X-GY)%20%3D%20(X%5EH-Y%5EH%20G%5EH)(X-GY)%20%5C%5C%0A%0A%26%3DX%5EH%20X%20-%20X%5EH%20GY%20-%20Y%5EH%20G%5EH%20X%20%2B%20Y%5EH%20G%5EH%20GY%20%20%20%0A%0A%5Cend%7Baligned%7D%20%20%5Ctag%205

对公式 (5) 相对 G 求导,利用如下公式 (a,b 均为列向量,X 为矩阵)

%5Cfrac%7B%20%5Cpartial%7Ba%5EH%20X%20b%7D%20%7D%20%7B%20%5Cpartial%20X%7D%20%3D%20a%20b%5EH%20%20%5C%5C%0A%0A%5Cfrac%7B%20%5Cpartial%7Ba%5EH%20X%5EH%20b%7D%20%7D%20%7B%5Cpartial%20X%7D%20%3D%20b%20a%5EH%20%20%5C%5C%0A%0A%5C%5C%0A%0A%5Cfrac%7B%20%5Cpartial%7Ba%5EH%20X%20X%5EH%20b%7D%20%7D%20%7B%5Cpartial%20X%7D%20%3D%20(a%20b%5EH%20%2B%20ba%5EH)%20X%20%20%5C%5C%0A%0A%5Cfrac%7B%20%5Cpartial%7Ba%5EH%20X%5EH%20X%20b%7D%20%7D%20%7B%5Cpartial%20X%7D%20%3D%20X(a%20b%5EH%20%2B%20ba%5EH)

则:

%5Cfrac%7B%5Cpartial%20X%5EH%20GY%7D%7B%5Cpartial%20G%7D%20%3D%20X%20Y%5EH%20%20%5C%5C%0A%0A%5Cfrac%7B%5Cpartial%20Y%5EH%20G%5EH%20X%7D%7B%5Cpartial%20G%7D%20%3D%20XY%5EH%20%5C%5C%0A%0A%5C%5C%0A%0A%5Cfrac%7B%20%5Cpartial%20Y%5EH%20G%5EH%20GY%7D%7B%5Cpartial%20G%7D%20%3D%202G%20YY%5EH

那么

%5Cfrac%7B%5Cpartial%20E_%7BX%2CW%7D%7C%7CX-GY%7C%7C%5E2%7D%7B%5Cpartial%20G%7D%20%3D%20E_%7BX%2CW%7D(-2XY%5EH%20%2B%202%20G%20YY%5EH)%20%5Ctag%206

令公式(6)等于 0 ,即相当于求解公式 (5) 的极值问题,则:

G%20E_%7BX%2CW%7D(YY%5EH)%20%3D%20E_%7BX%2CW%7D(XY%5EH)%20%20%5Ctag%207

其中两个相关矩阵的计算如下:

%0A%5Cbegin%7Baligned%7D%0A%0AE_%7BX%2CW%7D(YY%5EH)%20%26%3D%20E_%7BX%2CW%7D(HX%2BW)(HX%2BW)%5EH%20%5C%5C%0A%0A%26%3D%20E_%7BX%2CW%7D%20(HX%2BW)(X%5EH%20H%5EH%2BW%5EH)%20%5C%5C%0A%0A%26%3D%20E_%7BX%2CW%7D%20(HXX%5EH%20H%5EH%2BHXW%5EH%20%2B%20WX%5EH%20H%5EH%20%2B%20W%20W%5EH)%20%5C%5C%0A%0A%26%3D%20E_%7BX%2CW%7D%20(HXX%5EH%20H%5EH)%2BE_%7BX%2CW%7D%20(HXW%5EH)%20%2B%20E_%7BX%2CW%7D%20(WX%5EH%20H%5EH)%20%2B%20E_%7BX%2CW%7D%20(W%20W%5EH)%20%5C%5C%0A%0A%26%3D%20%20H%20E_%7BX%2CW%7D(XX%5EH)%20H%5EH%2BH%20E_%7BX%2CW%7D%20(XW%5EH)%20%2B%20E_%7BX%2CW%7D%20(WX%5EH)%20H%5EH%20%2B%20E_%7BX%2CW%7D%20(W%20W%5EH)%20%5C%5C%0A%0A%26%3D%20H%20I%20H%5EH%20%2B%20H%20%5Ctimes0%2B%200%20%5Ctimes%20H%5EH%20%2B%20%5Csigma%5E2%20I%20%5C%5C%0A%0A%26%3D%20H%20H%5EH%20%2B%20%5Csigma%5E2%20I%0A%0A%5Cend%7Baligned%7D%20%20%5Ctag%208

以及

%0A%5Cbegin%7Baligned%7D%0A%0AE_%7BX%2CW%7D(XY%5EH)%20%26%3D%20E_%7BX%2CW%7DX(HX%2BW)%5EH%20%5C%5C%0A%0A%26%3D%20E_%7BX%2CW%7D%20X(X%5EH%20H%5EH%2BW%5EH)%20%5C%5C%0A%0A%26%3D%20E_%7BX%2CW%7D%20(XX%5EH%20H%5EH%2BXW%5EH)%20%5C%5C%0A%0A%26%3D%20E_%7BX%2CW%7D%20(XX%5EH%20H%5EH)%2BE_%7BX%2CW%7D%20(XW%5EH)%20%5C%5C%0A%0A%26%3D%20%20E_%7BX%2CW%7D(XX%5EH)%20H%5EH%2B%200%20%20%5C%5C%0A%0A%26%3D%20I%20H%5EH%20%20%5C%5C%0A%0A%26%3D%20H%5EH%20%0A%0A%5Cend%7Baligned%7D%20%20%5Ctag%209



把 (8) (9) 代入 (7) :

G(H%20H%5EH%20%2B%20%5Csigma%5E2%20I)%3DH%5EH%20%5Ctag%20%7B10%7D

最终得到:

G%3DH%5EH%20(H%20H%5EH%20%2B%20%5Csigma%5E2%20I)%5E%7B-1%7D%20%5Ctag%20%7B11%7D

需要注意的是,公式 (11) 也可以写成:

G%3D%20(H%5EH%20H%20%2B%20%5Csigma%5E2%20I)%5E%7B-1%7D%20H%5EH%20%5Ctag%20%7B12%7D

即:

H%5EH%20(H%20H%5EH%20%2B%20%5Csigma%5E2%20I)%5E%7B-1%7D%20%3D%20(H%5EH%20H%20%2B%20%5Csigma%5E2%20I)%5E%7B-1%7D%20H%5EH%20%20%5Ctag%7B13%7D


MIMO MMSE 的公式的推导的评论 (共 条)

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