信道估计---LS、MMSE、LMMSE准则

本期目录

引言
基本假设
LS信道估计
- LS信道估计工程实现
MMSE信道估计
LMMSE信道估计
- LMMSE实现

引言

信道估计主要分为非盲信道估计和盲信道估计。顾名思义，非盲信道估计需要使用基站和接收机均已知的导频序列进行信道估计，并使用不同的时频域插值技术来估计导频之间或者符号之间的子载波上的信道响应。目前主要使用的非盲信道估计包括最小二乘（LS）信道估计、最小均方误差（MMSE）信道估计、基于DFT的信道估计以及基于判决反馈信道估计等；而盲信道估计不需要已经已知的导频序列，主要包括基于最大期望的信道估计、基于子空间的信道估计技术等。本文主要介绍非盲信道估计

训练符号可以用于信道估计，通常能够提供较好的性能。然而，除了发射数据符号外，还需要发射前导或导频信号,由此产生的负荷会降低传输效率。当可以获得训练符号时,最小二乘(LS)和最小均方误差(MMSE)技术被广泛应用于信道估计。

基本假设

假设所有子载波是正交的，即没有载频间干扰（ICI），那么可以将N个子载波的训练符号表示成矩阵形式:
X=[X[0]0...00X[1]...0............0...0X[N−1]]X= \begin{bmatrix} X[0] & 0 &... &0 \\ 0 & X[1] &... &0\\ ... &... &... &...\\ 0 &... &0 &X[N-1] \end{bmatrix} X=⎣⎢⎢⎡X[0]0...00X[1]...............000...X[N−1]⎦⎥⎥⎤
其中，X[k]X[k]X[k]表示第kkk个子载波上的导频信号，满足E{X[k]}=0E\{X[k]\}=0E{X[k]}=0 , Var{X[k]}=σ2Var\{X[k]\}=\sigma^2Var{X[k]}=σ2，k=0,1,2.…,N−1k=0,1,2.…,N-1k=0,1,2.…,N−1。因为假设所有的子载波都是正交的，所以X是一个对角矩阵。给定第kkk 个载波的信道增益H[k]H[k]H[k]，接收到的训练信号Y[k]Y[k]Y[k]能够表示为
Y=[Y[0]Y[1]...Y[N−1]]=[X[0]0...00X[1]...0............0...0X[N−1]][H[0]H[1]...H[N−1]]+[Z[0]Z[1]...Z[N−1]]Y= \begin{bmatrix} Y[0] \\ Y[1] \\ ...\\ Y[N-1] \\ \end{bmatrix} =\begin{bmatrix} X[0] & 0 &... &0 \\ 0 & X[1] &... &0\\ ... &... &... &...\\ 0 &... &0 &X[N-1] \end{bmatrix} \begin{bmatrix} H[0] \\ H[1] \\ ...\\ H[N-1] \\ \end{bmatrix} +\begin{bmatrix} Z[0] \\ Z[1] \\ ...\\ Z[N-1] \\ \end{bmatrix} Y=⎣⎢⎢⎡Y[0]Y[1]...Y[N−1]⎦⎥⎥⎤=⎣⎢⎢⎡X[0]0...00X[1]...............000...X[N−1]⎦⎥⎥⎤⎣⎢⎢⎡H[0]H[1]...H[N−1]⎦⎥⎥⎤+⎣⎢⎢⎡Z[0]Z[1]...Z[N−1]⎦⎥⎥⎤

其中，HHH为信道向量，H=[H[0,H[1],...,H[N−1]]TH=[H[0,H[1],...,H[N-1]]^TH=[H[0,H[1],...,H[N−1]]T；ZZZ为噪声向量Z=[Z[0,Z[1],...,Z[N−1]]TZ=[Z[0,Z[1],...,Z[N-1]]^TZ=[Z[0,Z[1],...,Z[N−1]]T，满足E{Z[k]}=0E\{Z[k]\}=0E{Z[k]}=0 , Var{Z[k]}=σz2Var\{Z[k]\}=\sigma_z^2Var{Z[k]}=σz2，k=0,1,2.…,N−1k=0,1,2.…,N-1k=0,1,2.…,N−1。在下文中用H^\hat{H}H^表示对信道HHH的估计。

LS信道估计

LS信道估计是根据最小二乘准则的信道估计方法。在无线系统中，接收信号可表示为：
　　Y=XH+ZY=XH+ZY=XH+Z
其中，X表示原始发射信号矢量(即导频信号)、H表示信道响应矢量、Z表示噪声矢量，Y表示接收信号矢量。我们的估计信道可以表示为:
Y=XH^Y=X\hat{H}Y=XH^

根据最小二乘准则，有如下目标函数:
　　J(H^)LS=∣∣Y−XH^∣∣2=(Y−XH^)H(Y−XH^)=YHY−YHXH^−H^HXHY+H^HXHXH^\begin{aligned} 　　J(\hat{H})_{LS}&= || Y-X\hat{H}||\raisebox{0.5em}{2}\\ 　　&=(Y-X\hat{H})\raisebox{0.5em}{H}(Y-X\hat{H})\\ 　　&=Y\raisebox{0.5em}{H}Y-Y\raisebox{0.5em}{H}X\hat{H}-\hat{H}\raisebox{0.5em}{H}X\raisebox{0.5em}{H}Y+\hat{H}\raisebox{0.5em}{H}X\raisebox{0.5em}{H}X\hat{H} 　　\end{aligned}　　J(H^)LS　　　　=∣∣Y−XH^∣∣2=(Y−XH^)H(Y−XH^)=YHY−YHXH^−H^HXHY+H^HXHXH^　　
　　
　　为了使得误差平方和最小，对上述目标函数求关于H^\hat{H}H^的一阶偏导数：
　　∂(J(H^))∂(H^)=∂(YHY−YHXH^−H^HXHY+H^HXHXH^)∂(H^)=∂(−YHXH^))∂(H^)+∂(−H^HXHY))∂(H^)+∂(H^HXHXH^))∂(H^)=−YHX−(XHY)H+(XHXH^)H+H^HXHX=−2YHX+2H^HXHX\begin{aligned} \frac{\partial(J(\hat{H}))}{\partial(\hat{H})}&=\frac{\partial(Y\raisebox{0.5em}{H}Y-Y\raisebox{0.5em}{H}X\hat{H}-\hat{H}\raisebox{0.5em}{H}X\raisebox{0.5em}{H}Y+\hat{H}\raisebox{0.5em}{H}X\raisebox{0.5em}{H}X\hat{H})}{\partial(\hat{H})}\\ &= \frac{\partial(-Y\raisebox{0.5em}{H}X\hat{H}))}{\partial(\hat{H})}+\frac{\partial(-\hat{H}\raisebox{0.5em}{H}X\raisebox{0.5em}{H}Y))}{\partial(\hat{H})}+\frac{\partial(\hat{H}\raisebox{0.5em}{H}X\raisebox{0.5em}{H}X\hat{H}))}{\partial(\hat{H})}\\ &=-Y\raisebox{0.5em}{H}X-(X\raisebox{0.5em}{H}Y)\raisebox{0.5em}{H}+(X\raisebox{0.5em}{H}X\hat{H})\raisebox{0.5em}{H}+\hat{H}\raisebox{0.5em}{H}X\raisebox{0.5em}{H}X\\ &=-2Y\raisebox{0.5em}{H}X+2\hat{H}\raisebox{0.5em}{H}X\raisebox{0.5em}{H}X 　　\end{aligned}∂(H^)∂(J(H^))=∂(H^)∂(YHY−YHXH^−H^HXHY+H^HXHXH^)=∂(H^)∂(−YHXH^))+∂(H^)∂(−H^HXHY))+∂(H^)∂(H^HXHXH^))=−YHX−(XHY)H+(XHXH^)H+H^HXHX=−2YHX+2H^HXHX　　
　　
　　令一阶偏导数为0，则有
　　2YHX+2H^HXHX=0H^HXHX=YHXXHXH^=XHYH^=(XHX)-1XHY\begin{aligned} 　　2Y\raisebox{0.5em}{H}X+2\hat{H}\raisebox{0.5em}{H}X\raisebox{0.5em}{H}X&=0\\ 　　\hat{H}\raisebox{0.5em}{H}X\raisebox{0.5em}{H}X&=Y\raisebox{0.5em}{H}X\\ 　　X\raisebox{0.5em}{H}X\hat{H}&=X\raisebox{0.5em}{H}Y\\ 　　\hat{H}&=(X\raisebox{0.5em}{H}X)\raisebox{0.5em}{-1}X\raisebox{0.5em}{H}Y 　　\end{aligned}　　2YHX+2H^HXHX　　H^HXHX　　XHXH^　　H^=0=YHX=XHY=(XHX)-1XHY　　
　　所以，当H^=(XHX)-1XHY\hat{H}=(X\raisebox{0.5em}{H}X)\raisebox{0.5em}{-1}X\raisebox{0.5em}{H}YH^=(XHX)-1XHY时，所得估计的误差平方和最小，此时得到的估计信道H^\hat{H}H^是LS信道估计的解(最小范数解或最佳逼近解，广义逆矩阵知识)
　　由上式可以得到目标函数的最小值，即LS信道估计的解为：
　　H^=(XHX)-1XHY=X-1Y(仅当X是满秩矩阵时候才可这么化简，否则XHX作为整体无法被分割)=X−1(XH+Z)=H+X−1Z\begin{aligned} 　　\hat{H}&=(X\raisebox{0.5em}{H}X)\raisebox{0.5em}{-1}X\raisebox{0.5em}{H}Y\\ 　　&=X\raisebox{0.5em}{-1}Y\\ 　　&(仅当X是满秩矩阵时候才可这么化简，否则X^HX作为整体无法被分割)\\ 　　&=X^{-1}(XH+Z)\\ 　　&=H+X^{-1}Z 　　\end{aligned}　　H^　　　　　　　　=(XHX)-1XHY=X-1Y(仅当X是满秩矩阵时候才可这么化简，否则XHX作为整体无法被分割)=X−1(XH+Z)=H+X−1Z　　

LS信道估计算法，实现比较简单，计算复杂度低，但是忽略了噪声的影响。LS信道估计的均方误差（MSE）为:
MSELS=E{(H−H^)H(H−H^)}=E{(H−X−1Y)H(H−X−1Y)}=E{(X−1Z)H(X−1Z)}=E{ZH(XXH)−1Z}=σz2σx2=1SNR\begin{aligned} MSE_{LS} &= E\{(H-\hat{H})^{H}(H-\hat{H})\}\\ &=E\{(H-X^{-1}Y)^H(H-X^{-1}Y)\}\\ &=E\{(X^{-1}Z)^H(X^{-1}Z)\}\\ &=E\{Z^H(XX^H)^{-1}Z\}\\ &= \frac{\sigma^2_z}{\sigma^2_x}\\ &=\frac{1}{SNR} \end{aligned}MSELS=E{(H−H^)H(H−H^)}=E{(H−X−1Y)H(H−X−1Y)}=E{(X−1Z)H(X−1Z)}=E{ZH(XXH)−1Z}=σx2σz2=SNR1

从上式可以看出LS估计信道的均方误差MSE与信噪比SNR成反比，意味着LS信道估计增强了噪声，尤其在信道处于深度衰落时，即低SNR的情况下更是如此，信道估计的精度会受到较大影响。虽然如此，但由于实现简单，此方法仍在实际中大规模使用。

LS信道估计工程实现

在实际工程应用中，由于矩阵求逆运算量很大，遇到大规模矩阵无法求解，因此可以采取其他求解方法。

方法一：
Y=XH^XHY=XHXH^H^=XHY∣∣X∣∣2\begin{aligned} Y&=X\hat{H}\\ X^HY&=X^HX\hat{H}\\ \hat{H}&=\frac{X^HY}{||X||^2} \end{aligned}YXHYH^=XH^=XHXH^=∣∣X∣∣2XHY

方法二：
令H^LS[k]\hat{H}_{LS}[k]H^LS[k]表示H^LS\hat{H}_{LS}H^LS中的元素，k=0,1,2,...,N−1（N表示含导频载波个数）k=0,1,2,...,N-1（N表示含导频载波个数）k=0,1,2,...,N−1（N表示含导频载波个数）。若是没有载波间干扰（ICI），就可以直接求出每个子载波上的LS信道估计：

H^LS[k]=Y[k]X[k],k=0,1,2,...,N−1\hat{H}_{LS}[k]=\frac{Y[k]}{X[k]}, k=0,1,2,...,N-1H^LS[k]=X[k]Y[k],k=0,1,2,...,N−1

MMSE信道估计

MMSE估计是在LS估计的基础上增加了加权矩阵W，改用最小均方误差准则进行优化。
考虑LS估计的最优解，即H^LS=X−1Y\hat{H}_{LS}=X^{-1}YH^LS=X−1Y。利用加权矩阵W，定义MMSE的估计为H^=H^MMSE=WH^LS\hat{H}=\hat{H}_{MMSE}=W\hat{H}_{LS}H^=H^MMSE=WH^LS

根据最小均方误差准则，有如下目标函数:
J(H^)MSE=E{∣∣e∣∣2}=E{∣∣H−H^∣∣2}=E{∣∣H−WH^LS∣∣2}\begin{aligned} J(\hat{H})_{MSE}&=E\{ ||e||^2\}\\ &=E\{||H-\hat{H}||^2\} \\ &=E\{||H-W\hat{H}_{LS}||^2\} \\ \end{aligned}J(H^)MSE=E{∣∣e∣∣2}=E{∣∣H−H^∣∣2}=E{∣∣H−WH^LS∣∣2}
由矩阵论知识可知，当e=H−H^=H−WH^LSe=H-\hat{H}=H-W\hat{H}_{LS}e=H−H^=H−WH^LS与H^LS\hat{H}_{LS}H^LS正交时，J(H^)MSEJ(\hat{H})_{MSE}J(H^)MSE可取得最小值。即满足
E{eH^LSH}=E{(H−H^)H^LSH}=E{(H−WH^LS)H^LSH}=E{HH^LSH}−WE{H^LSH^LSH}=RHH^LS−WRH^LSH^LS=0\begin{aligned} E\{e\hat{H}_{LS}^H\}&=E\{(H-\hat{H})\hat{H}_{LS}^H\}\\ &=E\{(H-W\hat{H}_{LS})\hat{H}_{LS}^H\}\\ &=E\{H\hat{H}_{LS}^H\}-WE\{\hat{H}_{LS}\hat{H}_{LS}^H\}\\ &=R_{H\hat{H}_{LS}}-WR_{\hat{H}_{LS}\hat{H}_{LS}}=0 \end{aligned}E{eH^LSH}=E{(H−H^)H^LSH}=E{(H−WH^LS)H^LSH}=E{HH^LSH}−WE{H^LSH^LSH}=RHH^LS−WRH^LSH^LS=0
即
W=RHH^LSRH^LSH^LS−1W=R_{H\hat{H}_{LS}}R_{\hat{H}_{LS}\hat{H}_{LS}}^{-1}W=RHH^LSRH^LSH^LS−1
其中RH^LSH^LSR_{\hat{H}_{LS}\hat{H}_{LS}}RH^LSH^LS为矩阵H^LS\hat{H}_{LS}H^LS的自相关矩阵，考虑到H^LS=X−1Y=H+X−1Z\hat{H}_{LS}=X^{-1}Y=H+X^{-1}ZH^LS=X−1Y=H+X−1Z，则有
RH^LSH^LS=E{H^LSH^LSH}=E{X−1Y(X−1Y)H}=E{(H+X−1Z)(H+X−1Z)H}=E{HHH+X−1ZHH+HZH(X−1)H+X−1ZZH(X−1)H}=E{HHH}+E{X−1ZZH(X−1)H}=RHH+σz2σx2I\begin{aligned} R_{\hat{H}_{LS}\hat{H}_{LS}}&=E\{\hat{H}_{LS}\hat{H}_{LS}^H\}\\ &=E\{X^{-1}Y(X^{-1}Y)^H\}\\ &=E\{(H+X^{-1}Z)(H+X^{-1}Z)^H\}\\ &=E\{HH^H+X^{-1}ZH^H+HZ^H(X^{-1})^H+X^{-1}ZZ^H(X^{-1})^H\}\\ &= E\{HH^H\}+E\{X^{-1}ZZ^H(X^{-1})^H\}\\ &= R_{HH}+\frac{\sigma_z^2}{\sigma_x^2}I \end{aligned}RH^LSH^LS=E{H^LSH^LSH}=E{X−1Y(X−1Y)H}=E{(H+X−1Z)(H+X−1Z)H}=E{HHH+X−1ZHH+HZH(X−1)H+X−1ZZH(X−1)H}=E{HHH}+E{X−1ZZH(X−1)H}=RHH+σx2σz2I
RHH^LSR_{H\hat{H}_{LS}}RHH^LS为矩阵HHH和H^LS\hat{H}_{LS}H^LS的互相关矩阵。由上式，我们最终可得：
H^=WH^LS=RHH^LSRH^LSH^LS−1H^LS=RHH^LS(RHH+σz2σx2I)−1H^LS\begin{aligned} \hat{H}&=W\hat{H}_{LS}=R_{H\hat{H}_{LS}}R_{\hat{H}_{LS}\hat{H}_{LS}}^{-1}\hat{H}_{LS}\\ &=R_{H\hat{H}_{LS}}(R_{HH}+\frac{\sigma_z^2}{\sigma_x^2}I)^{-1}\hat{H}_{LS} \end{aligned}H^=WH^LS=RHH^LSRH^LSH^LS−1H^LS=RHH^LS(RHH+σx2σz2I)−1H^LS

LMMSE信道估计

LMMSE信道估计是在MMSE信道估计的基础上做了一次线性平滑。考虑到MMSE需要计算(RHH+σz2σx2I)−1(R_{HH}+\frac{\sigma_z^2}{\sigma_x^2}I)^{-1}(RHH+σx2σz2I)−1，随着噪声变化和输入信号x的变化，该矩阵求逆的运算量很大而且需要不停地重新计算，极大地占据了计算资源，实时性很差。
因此LMMSE采用了一种线性最小均方误差的估算方法，用期望值的形式代替了σz2σx2\frac{\sigma_z^2}{\sigma_x^2}σx2σz2，可将MMSE的估计准则简化为：
H^(k)=RHH^LS(RHH+βSNR)−1H^LS(k)\begin{aligned} \hat{H}(k)=R_{H\hat{H}_{LS}}(R_{HH}+\frac{\beta}{SNR})^{-1}\hat{H}_{LS}(k) \end{aligned}H^(k)=RHH^LS(RHH+SNRβ)−1H^LS(k)
如此一来，因为RHHR_{HH}RHH和βSNR\frac{\beta}{SNR}SNRβ都是常数，在一次信息传输过程中保持不变，只需要计算一次其逆矩阵即可，大大减小了计算量。
β\betaβ为信道调制类型参数，且不同类型的调制信道其调制参数也不同。
若信道采用 16QAM 调制，则β\betaβ取 179\frac{17}{9}917。若采用 QPSK 调制，则β\betaβ取 1。

LMMSE信道估计的协方差矩阵σ2\sigma^2σ2可表示为
σ2=RHH^LS−[RHH^LS(RHH+βSNR)−1H^LS]\sigma^2=R_{H\hat{H}_{LS}}-[R_{H\hat{H}_{LS}}(R_{HH}+\frac{\beta}{SNR})^{-1}\hat{H}_{LS}] σ2=RHH^LS−[RHH^LS(RHH+SNRβ)−1H^LS]
LMMSE信道估计的均方误差可表示为：
MSELMMSE=1Ntr{σ2}MSE_{LMMSE}=\frac{1}{N}tr\{{\sigma^2}\} MSELMMSE=N1tr{σ2}
tr{}tr\{\}tr{}表示求迹运算

LMMSE实现

见最新文章