考虑如下的MIMO信道模型
Y m × 1 = H m × n X n × 1 + N m × 1 , N ∼ C N ( 0 , σ 2 I m ) , E [ X ] = 0 n , E [ X 2 ] = I n 。 Y_{m\times1}=H_{m\times n}X_{n\times1}+N_{m\times1},\\ N\sim CN(0,\sigma^2I_m),\\ E[X]=0_n,E[X^2]=I_n。 Ym×1=Hm×nXn×1+Nm×1,N∼CN(0,σ2Im),E[X]=0n,E[X2]=In。
如果做MMSE(最小均方误差)检测,检测算法不一定是线性的,也即不一定能用矩阵表示算法,需要写成函数形式
Y M M S E = f ( Y ) = E [ X ∣ Y = y ] = ∫ x P X ( x ) C N ( y − x ∣ 0 , σ 2 I m ) ∫ P X ( x ) C N ( y − x ∣ 0 , σ 2 I m ) d x d x , Y_{MMSE}=f(Y)=E[X|Y=y] \\ =\int x\frac{P_X(x)CN(y-x|0,\sigma^2I_m)}{\int P_X(x)CN(y-x|0,\sigma^2I_m)dx}dx, YMMSE=f(Y)=E[X∣Y=y]=∫x∫PX(x)CN(y−x∣0,σ2Im)dxPX(x)CN(y−x∣0,σ2Im)dx,
也就是条件期望可以达到MMSE。通常情况下(发送为高斯信源除外),条件期望计算很困难,因此MMSE检测很难实现。
LMMSE形式1
如果对他进行Linear MMSE(LMMSE)检测,那么
D = H H ( H H H + σ 2 I m ) − 1 , X ^ = D Y = H H ( H H H + σ 2 I m ) − 1 Y , D=H^H(HH^H+\sigma^2I_m)^{-1},\\ \hat X=DY=H^H(HH^H+\sigma^2I_m)^{-1}Y, D=HH(HHH+σ2Im)−1,X^=DY=HH(HHH+σ2Im)−1Y, 其中 D D D是LMMSE检测矩阵。具体推导如下,
e = X ^ − X = D Y − X M S E = E [ e H e ] = t r { E [ e e H ] } = t r { E [ ( D Y − X ) ( D Y − X ) H ] } = t r { E [ D Y Y H D H − X Y H D H − D Y X H + X X H ] } = t r { D R Y Y D H } − t r { R X Y D H } − t r { D R Y X } + t r { R X X } = T 1 − T 2 − T 3 + T 4 e={\hat X}-X=DY-X\\ MSE=E[e^He]=tr\{E[ee^H]\}\\ =tr\{E[(DY-X)(DY-X)^H]\}\\ =tr\{E[DYY^HD^H-XY^HD^H-DYX^H+XX^H]\}\\ =tr\{DR_{YY}D^H\}-tr\{R_{XY}D^H\}-tr\{DR_{YX}\}+tr\{R_{XX}\}\\ =T_1-T_2-T_3+T_4 e=X^−X=DY−XMSE=E[eHe]=tr{
E[eeH]}=tr{
E[(DY−X)(DY−X)H]}=tr{
E[DYYHDH−XYHDH−DYXH+XXH]}=tr{
DRYYDH}−tr{
RXYDH}−tr{
DRYX}+tr{
RXX}=T1−T2−T3+T4
根据《Complex valued matrix derivatives》中,表4.3的求导公式,可以得到
0 = ∂ M S E ∂ D = ∂ T 1 ∂ D − ∂ T 2 ∂ D − ∂ T 3 ∂ D + ∂ T 4 ∂ D , ∂ T 1 ∂ D = D ∗ R Y Y T ∂ T 2 ∂ D = 0 ∂ T 3 ∂ D = R Y X T ∂ T 4 ∂ D = 0 0=\frac{\partial MSE}{\partial D}=\frac{\partial T_1}{\partial D}-\frac{\partial T_2}{\partial D}-\frac{\partial T_3}{\partial D}+\frac{\partial T_4}{\partial D},\\ \frac{\partial T_1}{\partial D}=D^*R_{YY}^T\\ \frac{\partial T_2}{\partial D}=0\\ \frac{\partial T_3}{\partial D}=R_{YX}^T\\ \frac{\partial T_4}{\partial D}=0 0=∂D∂MSE=∂D∂T1−∂D∂T