CN103051572B - 一种lte信道估计方法和装置 - Google Patents

Abstract

本发明公开了一种LTE信道估计方法和装置。其中，所述方法包括：获取LTE信道的噪声功率信号功率以及导频点的发送信号矢量X、接收端的接收信号矢量Y；根据所述以及X、Y，应用最小均方误差MMSE准则<maths num="0001"></maths>计算得到信道频率响应估计值，其中，信道协方差矩阵R_hh为降低维数的简化矩阵，F是傅立叶变换矩阵。通过本发明技术方案，可以降低LTE信道估计的复杂度。

Description

一种LTE信道估计方法和装置

技术领域

本发明涉及通信处理技术领域，更具体地说，涉及一种LTE信道估计方法和装置。

背景技术

LTE(LongTermEvolution，长期演进技术)是3GPP(The3rdGenerationPartnershipProject，第三代合作伙伴计划)制定的移动通信系统标准。

LTE系统上行传输方案采用SC-FDMA(Signal-CarrierFrequencyDivisionMultipleAccess，带循环前缀的单载波频分多址)技术，而下行传输方案则采用OFDMA(OrthogonalFrequencyDivisionMultiplexingAccess，带循环前缀的正交频分多址)技术；此外，同时支持FDD(FrequencyDivisionDuplexing，频分双工)和TDD(TimeDivisionDuplexing，时分双工)技术。

无线通信系统的性能在很大程度上受到无线信道的影响，如：阴影衰落和频率选择性衰落等，使得发射机和接收机之间的传播路径非常复杂。无线信道不像有线信道固定并可预见，而是有很大的随机性，这就对接收机的设计提出了很大的挑战。为了能在接收端准确地恢复发射端的发送信号，通常需要采用各种措施来抵抗多径效应对传输信号的影响，信道估计技术的实现需要知道无线信道的信息，进而需要对信道的参数进行估计。因此，信道参数估计是实现无线通信系统的一项关键技术。能否获得较为详细且准确的信道信息，从而在接收端正确地解调出发射信号，是衡量一个无线通信系统性能的重要指标。

目前最常用的信道估计方法是LS(LeastSquare，最小平方)法、基于DFT(DiscreteFourierTransform，离散傅里叶变换)的信道估计和MMSE(MinimumMeanSquareError，最小均方误差)法。

然而，通过发明人的研究发现，现有的信道估计方法存在以下缺陷：

由于忽略了噪声的影响，使得LS法计算量小，但是信道估计值对噪声干扰的影响比较敏感，在信噪比较低的情况下，信道估计较差；

MMSE法性能优于LS法，但有很高的运算复杂度，这是由于，MMSE基于对导频的IFFT来计算信道时域冲激响应，由于实际的多径信道不可能每一多径时延都恰好落在接收机采样间隔的整数倍上，这必然导致多径在无限长时间上分裂，这种效应会导致估计的误差增大，同时由于需要做N点的FFT(FastFourierTransform，快速傅里叶变换)和IFFT(InverseFast-Fourier-Transformation，快速傅里叶反变换)，其复杂度也较高；

而基于DFT的方法则容易引起信息泄露。

发明内容

有鉴于此，本发明实施例提供一种LTE信道估计方法和装置，以降低LTE信道估计的复杂度。

本发明实施例提供一种LTE信道估计方法，包括：

获取LTE信道的噪声功率信号功率以及导频点的发送信号矢量X、接收端的接收信号矢量Y；

根据所述以及X、Y，应用最小均方误差MMSE准则 ${\hat{H}}_{MMSE}={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1}{X}^{H}Y,$ 计算得到信道频率响应估计值，其中，信道协方差矩阵R_hh为降低维数的简化矩阵，F是傅立叶变换矩阵；

所述信道协方差矩阵R_hh为简化后的对角矩阵，所述对角矩阵中的对角元素信道的自相关值由下式给出：

其中，D为衰减因子，L为信道的抽头数，对应于循环前缀的长度。

优选的，所述计算得到信道频率响应估计值，包括：

根据具体信噪比SNR，查找预先计算得到的相应的

$S={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1};$

将查找到的S值代入所述 ${\hat{H}}_{MMSE}={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1}{X}^{H}Y,$ 得到信道频率响应估计值。

优选的，所述计算得到信道频率响应估计值，包括：

将不同SINR分为M段区间，每段区间对应一个S_m，

$S_m={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+{\left(\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}\right)}_{m}I)}^{-1},$ 其中，m＝1,…,M，预置计算得到的相应的S_m；

根据具体SNR所属的区间，获取相应的S_m；

根据得到信道频率响应估计值。

一种LTE信道估计装置，包括：

参数获取模块，用于获取LTE信道的噪声功率信号功率以及导频点的发送信号矢量X、接收端的接收信号矢量Y；

计算模块，用于根据所述以及X、Y，应用最小均方误差MMSE准则 ${\hat{H}}_{MMSE}={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1}{X}^{H}Y,$ 计算得到信道频率响应估计值，其中，信道协方差矩阵R_hh为降低维数的简化矩阵，F是傅立叶变换矩阵；

所述信道协方差矩阵R_hh为简化后的对角矩阵，所述对角矩阵中的对角元素信道的自相关值由下式给出：

其中，D为衰减因子，L为信道的抽头数，对应于循环前缀的长度。

优选的，所述计算模块，包括：

第一预置单元，用于预置预先计算得到的对应具体SNR的

$S={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1};$

第一查找单元，用于根据具体SNR，查找预先计算得到的相应的

$S={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1};$

第一计算单元，用于将查找到的S值代入所述

${\hat{H}}_{MMSE}={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1}{X}^{H}Y,$ 得到信道频率响应估计值。

优选的，所述计算模块，包括：

第二预置单元，用于将不同SINR分为M段区间，每段区间对应一个S_m，其中，m＝1,…,M，预置计算得到的相应的S_m；

第二查找单元，用于根据具体SNR所属的区间，查找获取相应的S_m；

第二计算单元，用于根据得到信道频率响应估计值。

同现有技术相比，本发明提供的技术方案在采用MMSE准则进行信道估计的过程中，对MMSE准则所应用的信道协方差矩阵R_hh进行降低维数的简化处理，由于维数的降低，使得信道频率响应估计的运算量大大降低。

附图说明

为了更清楚地说明本发明实施例的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例提供的一种LTE信道估计方法流程示意图；

图2为本发明实施例提供的一种计算得到信道频率响应估计值的方式流程示意图；

图3为本发明实施例提供的另一种计算得到信道频率响应估计值的方式流程示意图；

图4为本发明实施例提供的一种LTE信道估计装置结构示意图；

图5为图4中一种计算模块结构示意图；

图6为图4中另一种计算模块结构示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明公开了一种LTE信道估计方法和装置，以降低LTE信道估计的复杂度。

首先对本发明提供的一种LTE信道估计方法进行说明，参照图1所示，所述方法包括：

步骤101、获取LTE信道的噪声功率信号功率以及导频点的发送信号矢量X、接收端的接收信号矢量Y；

步骤102、根据所述以及X、Y，应用MMSE(MinimalMeanSquaredError，最小均方误差)准则计算得到信道频率响应估计值，其中，信道协方差矩阵R_hh为降低维数的简化矩阵，F是傅立叶变换矩阵。

本发明实施例中，信道估计采用MMSE准则，MMSE准则的原理是使估计值和真实值之间的平方差期望最小，即：

$E\left\{\right|h-\hat{h}{|}^2\}=min---\left(1\right)$

当h具有高斯随机特性并与噪声向量n互相独立时，基于MMSE准则的CFR(ChannelFrequencyResponse，信道频率响应)估计值为：

${\hat{H}}_{MMSE}={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1}{X}^{H}Y---\left(2\right)$

其中，X是导频点的发送信号矢量，具体为一(NpxNp)的对角矩阵；Y是接收端的接收信号矢量，接收端FFT的输出，具体为一(Npx1)的矩阵；F是(NpxL)傅立叶变换矩阵；为噪声功率；为信号功率；I是(NpxNp)的单位矩阵；Np是导频符号的个数，即子载波的个数；L是信道的抽头数，对应于循环前缀的长度。

现有MMSE准则中，R_hh是(LxL)的信道协方差矩阵，而本发明实施例中，R_hh为降低维数的简化矩阵。

本发明提供的技术方案在采用MMSE准则进行信道估计的过程中，对MMSE准则所应用的信道协方差矩阵R_hh进行降低维数的简化处理，由于维数的降低，使得信道频率响应估计的运算量大大降低。

为了便于对本发明进一步的理解，下面结合本发明的具体实施方式对本发明进行详细描述。

X^HY生成列矢量Z，表示根据LS准则，利用Y/X得到的原始信道估计。对于导频符号复数相除等价于乘导频符号的复共轭。由此，等式(2)可以简化为:

(3)式中，S为(NpxNp)的矩阵，可以看做是对带有噪声的原始信道估计Z进行修正的平滑滤波器。

$S={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1}---\left(4\right)$

由(4)式可以看出，根据MMSE准则进行信道估计，需要获知信道协方差矩阵R_hh。

另外，前述傅立叶变换矩阵F可以用下式描述：

$F=\exp (j\&CenterDot;\frac{2\mathrm{\&pi;}}{N_{FFT}}\&CenterDot;k\&CenterDot;l^T)---\left(5\right)$

其中，k为频率索引矢量，k＝0…Np-1，l为时间索引矢量，l＝0…L-1。

时域信道冲击响应h可以表示为下式：

其中，α_n是幅度，τ_nT_s是第n条路径的延迟，T_s是采样间隔。通过对信道脉冲进行采样可以得到如下时域信道冲击响应的采样结果：

$h_k=\frac{1}{\sqrt{N}}\underset{n}{\&Sigma;}{\mathrm{\&alpha;}}_n{e}^{-j\frac{\mathrm{\&pi;}}{N}(k+(N-1){\mathrm{\&tau;}}_n)}\frac{sin\left({\mathrm{\&pi;\&tau;}}_n\right)}{sin\left(\frac{\mathrm{\&pi;}}{N}\right({\mathrm{\&tau;}}_n-k))}---\left(7\right)$

其中，k的取值为[1,N]。

需要说明的是，R_hh的元素值为h_m*h_n，R_hh的主要能量部分(多径信号的能量)在循环前缀之内，而低能量的部分则可以近似为零，因此，可以只保留信道前L个抽头作为有用项，即，对R_hh(m,n)＝0。

同时，简化R_hh为一个对角矩阵，其对角元素信道的自相关值按负指数函数分布，基本上符合典型的信道冲击响应功率延迟谱分布。

R_hh矩阵的元素由下式给出：

$r(l,l)=\frac{\exp (-D\&CenterDot;l)}{\underset{l=1}{\overset{L}{\&Sigma;}}\exp (-D\&CenterDot;l)},l=1,...,L---\left(8\right)$

衰减因子D的取值可以根据链路仿真结果进行确定，链路仿真可以考虑不同信道环境下各种延迟扩展情况。

MMSE准则算法中涉及到矩阵求逆的运算，这需要相当大的计算量，在实际应用中很难满足系统的实时性。为了解决这一问题，本发明实施例中，所述计算得到信道频率响应估计值的方式，如图2所示，可以包括以下步骤：

步骤201、根据具体SNR(Signal-to-NoiseRatio，信噪比)，查找预先计算得到的相应的 $S={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1}---\left(9\right);$

步骤202、将查找到的S值代入所述 ${\hat{H}}_{MMSE}={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1}{X}^{H}Y,$ 得到信道频率响应估计值。

该实施例中，可以针对一些典型的SNR值，可以预先计算出矩阵的逆，进而可以计算得到将计算得到的S值进行存储。存储时，可以建立S值与SNR值的映射关系。由此，对应某一SNR，可以通过查找映射关系，得到对应该SNR的S值，直接代入(1)式，即可得到信道频率响应估计值。通过该方式，可以进一步降低信道频率响应估计的运算复杂程度。

在本发明提供的另一个实施例中，所述计算得到信道频率响应估计值的方式，如图3所示，可以包括以下步骤：

步骤301、将不同SNR分为M段区间，每段区间对应一个S_m，

$S_m={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+{\left(\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}\right)}_{m}I)}^{-1}---\left(10\right),$

其中，m＝1,…,M，预置计算得到的相应的S_m；

步骤302、根据具体SNR所属的区间，获取相应的S_m；

步骤303、根据 ${\hat{H}}_{MMSE}=S_m{X}^{H}Y---\left(11\right),$ 得到信道频率响应估计值。

该实施例中，可以将不同SINR分为M段区间，每段区间对应一个S_m，其中，m＝1,…,M，具体地，可以根据实际应用场景进行SNR的划分，例如：SNR值的范围为0～20，可以每隔2划分出一个区间。预先计算出每个区间的S_m，并进行存储。对应具体的SNR，可以根据该SNR所属的区间，获取该区间对应的S_m。进而，根据H_MMSE＝S_mX^HY，直接代入对应的S_m值，即可得到信道频率响应估计值。通过该方式，也可进一步降低信道频率响应估计的运算复杂程度。

此外，还需要考虑Np的取值，当Np取值比较大时，所需要的运算量依旧很大。为此，可以一个PRB(PhysicalResourceBlock，物理资源块)，即12个子载波为单位进行信道估计，相当于将矩阵分成了维数较小的矩阵，可以降低运算量。

当计算获得单个PRB的信道估计之后，将所有的PRB排列起来就是完整的信道估计，由此，完整的信道估计为：

$\hat{H}={\&lsqb;{\hat{H}}_{MMSE,1}{\hat{H}}_{MMSE,2}...{\hat{H}}_{MMSE,K}\&rsqb;}^T---\left(12\right)$

其中，K为完整的频带所占用的PRB数。

相应上述方法实施例，本发明还提供了一种LTE信道估计装置，如图4所示，该装置具体可以包括：

参数获取模块401，用于获取LTE信道的噪声功率信号功率以及导频点的发送信号矢量X、接收端的接收信号矢量Y；

计算模块402，用于根据所述以及X、Y，应用最小均方误差MMSE准则 ${\hat{H}}_{MMSE}={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1}{X}^{H}Y,$ 计算得到信道频率响应估计值，其中，信道协方差矩阵R_hh为降低维数的简化矩阵，F是傅立叶变换矩阵。

本发明提供的技术方案在采用MMSE准则进行信道估计的过程中，对MMSE准则所应用的信道协方差矩阵R_hh进行降低维数的简化处理，由于维数的降低，使得信道频率响应估计的运算量大大降低。

R_hh的主要能量部分(多径信号的能量)在循环前缀之内，而低能量的部分则可以近似为零，因此，可以只保留信道前L个抽头作为有用项，即，对R_hh(m,n)＝0。由此，所述信道协方差矩阵R_hh为简化后的对角矩阵，所述对角矩阵中的对角元素信道的自相关值由下式给出：

其中，D为衰减因子，L为信道的抽头数，对应于循环前缀的长度。

在本发明的一个优选实施例中，如图5所示，所述计算模块402具体可以包括：

第一预置单元501，用于预置预先计算得到的对应具体SNR的

$S={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1};$

第一查找单元502，用于根据具体SNR，查找预先计算得到的相应的

$S={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1};$

第一计算单元503，用于将查找到的S值代入所述

${\hat{H}}_{MMSE}={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1}{X}^{H}Y,$ 得到信道频率响应估计值。

在本发明的另一个优选实施例中，如图6所示，所述计算模块402具体可以包括：

第二预置单元601，用于将不同SINR分为M段区间，每段区间对应一个其中，m＝1,…,M，预置计算得到的相应的S_m；

第二查找单元602，用于根据具体SNR所属的区间，查找获取相应的S_m；

第二计算单元603，用于根据得到信道频率响应估计值。

对于装置实施例而言，由于其基本相应于方法实施例，所以描述得比较简单，相关之处参见方法实施例的部分说明即可。以上所描述的设备实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-OnlyMemory，ROM)或随机存储记忆体(RandomAccessMemory，RAM)等。

对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明实施例的精神或范围的情况下，在其它实施例中实现。因此，本发明实施例将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (6)

1.一种LTE信道估计方法，其特征在于，包括：

获取LTE信道的噪声功率信号功率以及导频点的发送信号矢量X、接收端的接收信号矢量Y；

根据所述以及X、Y，应用最小均方误差MMSE准则 ${\hat{H}}_{MMSE}={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1}{X}^{H}Y,$ 计算得到信道频率响应估计值，其中，信道协方差矩阵R_hh为降低维数的简化矩阵，F是傅立叶变换矩阵；

所述信道协方差矩阵R_hh为简化后的对角矩阵，所述对角矩阵中的对角元素信道的自相关值由下式给出：

其中，D为衰减因子，L为信道的抽头数，对应于循环前缀的长度。

2.根据权利要求1所述的方法，其特征在于，所述计算得到信道频率响应估计值，包括：

根据具体信噪比SNR，查找预先计算得到的相应的 $S={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1};$

将查找到的S值代入所述 ${\hat{H}}_{MMSE}={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1}{X}^{H}Y,$ 得到信道频率响应估计值。

3.根据权利要求1所述的方法，其特征在于，所述计算得到信道频率响应估计值，包括：

将不同SINR分为M段区间，每段区间对应一个S_m，其中，m＝1,…,M，预置计算得到的相应的S_m；

根据具体SNR所属的区间，获取相应的S_m；

根据得到信道频率响应估计值。

4.一种LTE信道估计装置，其特征在于，包括：

参数获取模块，用于获取LTE信道的噪声功率信号功率以及导频点的发送信号矢量X、接收端的接收信号矢量Y；

计算模块，用于根据所述以及X、Y，应用最小均方误差MMSE准则计算得到信道频率响应估计值，其中，信道协方差矩阵R_hh为降低维数的简化矩阵，F是傅立叶变换矩阵；

所述信道协方差矩阵R_hh为简化后的对角矩阵，所述对角矩阵中的对角元素信道的自相关值由下式给出：

其中，D为衰减因子，L为信道的抽头数，对应于循环前缀的长度。

5.根据权利要求4所述的装置，其特征在于，所述计算模块，包括：

第一预置单元，用于预置预先计算得到的对应具体SNR的 $S={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1};$

第一查找单元，用于根据具体SNR，查找预先计算得到的相应的 $S={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1};$

第一计算单元，用于将查找到的S值代入所述 ${\hat{H}}_{MMSE}={\mathrm{FR}}_{hh}{F}^H{({\mathrm{FR}}_{hh}{F}^H+\frac{{{\widehat{\mathrm{\&sigma;}}}_n}^2}{{{\widehat{\mathrm{\&sigma;}}}_s}^2}I)}^{-1}{X}^{H}Y,$ 得到信道频率响应估计值。

6.根据权利要求5所述的装置，其特征在于，所述计算模块，包括：

第二预置单元，用于将不同SINR分为M段区间，每段区间对应一个S_m，其中，m＝1,…,M，预置计算得到的相应的S_m；

第二查找单元，用于根据具体SNR所属的区间，查找获取相应的S_m；

第二计算单元，用于根据得到信道频率响应估计值。