CN102932290A - Lte***干扰抑制接收方法及装置 - Google Patents

Abstract

本发明公开了一种LTE***干扰抑制接收方法及装置。基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并根据信道响应H计算所述各个资源块的平均协方差矩阵R（b）；依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）；确定所述接收信号Y中待补偿信号所在资源块b，依据所述待补偿信号所在资源块对应的干扰抑制补偿参数W（b）和所述信道响应H进行信号补偿，得到补偿后的接收信号Y’。基于上述方法及装置，能够获得更准确的平均协方差矩阵R（b），保证了干扰抑制接收的准确性以及整个通信网络的性能。

Description

LTE***干扰抑制接收方法及装置

技术领域

本发明涉及无线通信技术领域，更具体的说，是涉及一种LTE（Long TermEvolution system，长期演进***）***干扰抑制接收方法及装置。

背景技术

3GPP LTE项目是3G的演进，3GPP LTE是3G与4G技术之间的一个过渡，是当前被广泛认可的无线通信技术。

LTE***采用正交频分复用（Orthogonal Frequency Division Multiple，OFDM）的接入技术而构建的通信网络中，将用户的信息承载在相互正交的不同的载波上，能够有效对抗同小区的用户之间的干扰。但是同频组网下相邻小区干扰以及其他***干扰会带来很强的同频干扰。同频干扰会影响小区边缘用户的服务质量。

现有技术中，为了对抗同频干扰，对于配备了多根接收天线的终端，可以采用干扰抑制合并算法（Interference Rejection Combining，IRC）。

IRC是一种发射机或接收机信号处理的方法，可以提高LTE***中的信号接收质量。在基站的下行链路，通过IRC对接收天线接收到的信号进行干扰抑制合并后，再输出到均衡器进行信道解码。

以两接收天线为例，接收到的信号可以表示为：

$y=\left[\begin{array}{c}{y}_0\\ y_1\end{array}\right]=H*x+m---\left(1\right)$

在无干扰的情况下，m为高斯白噪声；在有干扰的情况下，m为有色噪声（即干扰加高斯白噪声），其协方差矩阵为R；x为发送信号，当发射天线为两个时，x的维数为2×1，当发射天线为单个时，x的维数为1×1；H为信道响应矩阵，当发射天线为两个时，H的维数为2×2，当发射天线为单个时，H的维数为2×1。

但是，现有的均衡器均是以白噪声为前提条件进行信道解码的。现有技术中利用IRC进行干扰抑制接收估计而得的协方差矩阵R并不准确，可能造成对接收的信号进行IRC处理后，不仅不能消除干扰，反而比处理前更糟。因此现有技术中，利用IRC进行干扰抑制接收，准确性低，严重时会影响整个通信网络的性能。

发明内容

有鉴于此，本发明提供了一种LTE***干扰抑制接收方法及装置，以克服现有技术中利用IRC进行干扰抑制接收，准确性低，严重时会影响整个通信网络的性能的问题。

为实现上述目的，本发明提供如下技术方案：

一种LTE***干扰抑制接收方法，包括：

基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并计算所述接收信号Y所在各个资源块的平均协方差矩阵R（b）；

依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）；

确定所述接收信号Y中待补偿信号所在资源块，依据所述待补偿信号所在资源块对应的干扰抑制补偿参数W（b）和所述信道响应H进行信号补偿，得到补偿后的接收信号Y’（b）。

优选地，在得到补偿后的接收信号之后，还包括：

选取所述待补偿信号所在资源块对应的等效高斯白噪声的方差NoiseVar（b），对所述补偿后的接收信号Y'（b）进行均衡处理。

优选地，基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并计算所述接收信号Y所在各个资源块的平均协方差矩阵R（b）的过程包括：

基于本地导频信号对接收信号Y进行信道估计，得到信道响应H；

依据所述信道响应H计算各个资源块的协方差矩阵；

判断各个所述资源块的协方差矩阵与N个相邻资源块的协方差矩阵之间的距离是否都小于预设距离门限；

如果是，则增多N，直到所述资源块的协方差矩阵与增多后的N个相邻资源块的协方差矩阵之间的距离不都小于预设距离门限，确定当前资源块与其增多后的N个相邻资源块构成的平均资源块集合；

如果否，则确定由当前资源块与N个相邻资源块构成的平均资源块集合；

根据所述平均资源块集合计算各个资源块平均协方差矩阵R（b）。

优选地，在所述判断各个所述资源块的协方差矩阵与N个相邻资源块之间的距离是否小于预设距离门限之前，所述计算所述估计信道响应H中各个资源块的协方差矩阵之后，还包括：

对各个所述资源块的协方差矩阵在时域上进行平滑处理。

优选地，基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并计算所述接收信号Y所在各个资源块的平均协方差矩阵R（b）的过程包括：

基于本地导频信号对接收信号Y进行信道估计，得到信道响应H；

依据所述信道响应H计算各个资源块的协方差矩阵；

根据预设平均资源块集合计算各个资源块平均协方差矩阵R（b）。

优选地，依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）的过程包括：将所述平均协方差矩阵R（b）转换为 $\mathrm{Rb}=\left[\begin{array}{cc}\mathrm{Rb}\left(\mathrm{1,1}\right)& \mathrm{Rb}\left(\mathrm{1,2}\right)\\ \mathrm{Rb}\left(\mathrm{2,1}\right)& \mathrm{Rb}\left(\mathrm{2,2}\right)\end{array}\right],$

其中，Rb(1,1)表示接收天线0的干扰噪声方差，Rb(1,2)表示接收天线0和接收天线1的干扰噪声协方差，Rb(2,1)表示接收天线1和接收天线0的干扰噪声协方差，Rb(2,2)表示接收天线1的干扰噪声方差；

当Rb（1,1）大于或等于Rb（2,2）时，依据 $W\left(b\right)=\left(\begin{array}{cc}\mathrm{\γ}& -\mathrm{\ρ}\\ 0& \sqrt{1-{\left|\mathrm{\ρ}\right|}^2}\end{array}\right),$ 获取干扰抑制补偿参数W（b）；

其中，γ、ρ为常量， $\mathrm{\γ}=\sqrt{\frac{\mathrm{Rb}\left(\mathrm{2,2}\right)}{\mathrm{Rb}\left(\mathrm{1,1}\right)}},$ $\mathrm{\ρ}=\frac{\mathrm{\γ}\·\mathrm{Rb}\left(\mathrm{1,2}\right)}{\mathrm{Rb}\left(\mathrm{2,2}\right)}.$

优选地，依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）的过程包括：将所述平均协方差矩阵R（b）转换为 $\mathrm{Rb}=\left[\begin{array}{cc}\mathrm{Rb}\left(\mathrm{1,1}\right)& \mathrm{Rb}\left(\mathrm{1,2}\right)\\ \mathrm{Rb}\left(\mathrm{2,1}\right)& \mathrm{Rb}\left(\mathrm{2,2}\right)\end{array}\right],$

其中，Rb(1,1)表示接收天线0的干扰噪声方差，Rb(1,2)表示接收天线0和接收天线1的干扰噪声协方差，Rb(2,1)表示接收天线1和接收天线0的干扰噪声协方差，Rb(2,2)表示接收天线1的干扰噪声方差；

当Rb（1,1）小于Rb（2,2）时，令R'=P·R·P，其中， $P=\left[\begin{array}{cc}0& 1\\ 1& 0\end{array}\right];$

依据 $V=\left(\begin{array}{cc}\mathrm{\γ}& -\mathrm{\ρ}\\ 0& \sqrt{1-{\left|\mathrm{\ρ}\right|}^2}\end{array}\right),$ 获取干扰抑制补偿参数W(b)＝P·V·P；

其中，γ、ρ为常量， $\mathrm{\γ}=\sqrt{\frac{\mathrm{Rb}\left(\mathrm{2,2}\right)}{\mathrm{Rb}\left(\mathrm{1,1}\right)}},$ $\mathrm{\ρ}=\frac{\mathrm{\γ}\·\mathrm{Rb}\left(\mathrm{1,2}\right)}{\mathrm{Rb}\left(\mathrm{2,2}\right)}.$

优选地，依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的等效高斯白噪声的方差NoiseVar（b）的过程包括：

依据NoiseVar(b)＝Rb(2,2)·(1-|ρ|²)得到对应各个所述资源块的等效高斯白噪声的方差NoiseVar（b）。

优选地，依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的等效高斯白噪声的方差NoiseVar（b）包括：

依据NoiseVar(b)＝Rb'(2,2)·(1-|ρ|²)得到对应各个所述资源块的等效高斯白噪声的方差NoiseVar（b）。

一种LTE***干扰抑制接收装置，包括：

协方差计算单元，用于基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并计算所述接收信号Y所在各个资源块的平均协方差矩阵R（b）；

干扰抑制参数计算单元，用于依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）；

预补偿单元，用于确定所述接收信号Y中待补偿信号所在资源块，依据所述待补偿信号所在资源块对应的干扰抑制补偿参数W（b）和所述信道响应H进行信号补偿，得到补偿后的接收信号Y’（b）。

优选地，上述装置还包括：

均衡单元，用于选取所述待补偿信号所在资源块对应的等效高斯白噪声的方差NoiseVar（b），对所述补偿后的接收信号Y'（b）进行均衡处理。

优选地，所述协方差计算单元包括：

信道估计单元，用于基于本地导频信号对接收信号Y进行信道估计，得到信道响应H；

每资源块计算单元，用于依据所述信道响应H计算各个资源块的协方差矩阵；

频域参数估计单元，用于判断各个所述资源块的协方差矩阵与N个相邻资源块的协方差矩阵之间的距离是否都小于预设距离门限；

频域平均单元，用于根据所述平均资源块集合计算各个资源块平均协方差矩阵R（b）。

优选地，还包括：

时域平滑单元，用于对各个所述资源块的协方差矩阵在时域上进行平滑处理。

优选地，还包括：

选择开关K，用于选择频域参数估计单元确定可平均资源块集合或选择预设平均资源块集合。

经由上述的技术方案可知，与现有技术相比，本发明公开了一种LTE***干扰抑制接收方法及装置。基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并根据信道响应H计算所述各个资源块的平均协方差矩阵R（b）；依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）；确定所述接收信号Y中待补偿信号所在资源块b，依据所述待补偿信号所在资源块对应的干扰抑制补偿参数W（b）和所述信道响应H进行信号补偿，得到补偿后的接收信号Y’。基于上述方法及装置，能够获得更准确的平均协方差矩阵R（b），保证了干扰抑制接收的准确性以及整个通信网络的性能。

附图说明

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

图1为本发明实施例一公开的一种LTE***干扰抑制接收方法流程图；

图2为本发明实施例二公开的另一种LTE***干扰抑制接收方法流程图；

图3为基站具备4个公共天线端口的下行公共参考信号在二维网格的资源位置的一个示意图；

图4为基站具备4个公共天线端口的下行公共参考信号在二维网格的资源位置的另一个示意图；

图5为本发明实施例二公开的另一种LTE***干扰抑制接收方法中计算平均协方差矩阵R的具体流程图；

图6为本发明实施例三公开的一种LTE***干扰抑制接收装置具体结构示意图；

图7为本发明实施例四公开的另一种LTE***干扰抑制接收装置具体结构示意图；

图8为本发明实施例四公开的另一种LTE***干扰抑制接收装置中平均协方差计算单元的具体结构示意图。

具体实施方式

为了引用和清楚起见，下文中使用的技术名词的说明、简写或缩写总结如下：

LTE：Long Term Evolution system，长期演进***；

IRC：Interference Rejection Combining，干扰抑制合并算法；

OFDM：Orthogonal Frequency Division Multiplexing，正交频分复用。

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

由背景技术可知，现有技术中利用IRC进行干扰抑制接收估计而得的协方差矩阵R并不准确，可能造成对接收的信号进行IRC处理后，不仅不能消除干扰，反而比处理前更糟。因此现有技术中，利用IRC进行干扰抑制接收，准确性低，严重时会影响整个通信网络的性能。

本发明公开了一种LTE***干扰抑制接收方法及装置。基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并根据信道响应H计算所述各个资源块的平均协方差矩阵R（b）；依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）；确定所述接收信号Y中待补偿信号所在资源块b，依据所述待补偿信号所在资源块对应的干扰抑制补偿参数W（b）和所述信道响应H进行信号补偿，得到补偿后的接收信号Y’。基于上述方法及装置，能够获得更准确的平均协方差矩阵R（b），保证了干扰抑制接收的准确性以及整个通信网络的性能。

关于该LTE***干扰抑制接收方法的具体流程以及该装置的具体结构将通过以下实施例进行详细描述。

实施例一

请参阅附图1，为本发明实施例一公开的一种LTE***干扰抑制接收方法流程图，主要包括以下步骤：

步骤S101：基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并计算所述接收信号Y所在各个资源块的平均协方差矩阵R（b），其中b为资源块索引；

LTE***采用OFDM技术，传输的资源位置由时间、频率二维资源网格指示，在网格中每个方格表示1个OFDM符号上的1个子载波。

LTE***中定义一个正常子帧的长度为1ms（毫秒），每个正常子帧由2个0.5ms的时隙构成，二维网格则由频域上12个子载波，时域上1个时隙构成并将其定义为资源块。接收信号Y分布在多个资源块中。

在执行步骤S101时，首先，基于本地导频信号，通过使用一些众所周知的信道估计方法，如最小均方误差估计、最小二乘、最大似然或者Kalman滤波，对接收信号Y进行信道估计，输出信道响应H，然后根据信道响应H计算各个资源块的平均协方差矩阵R（b）。

需要说明的是，基于干扰的频率选择性，在时域上，较短时间内干扰具有连续性的特点，也就是说相邻资源块的干扰情况相近。在计算各个资源块的平均协方差矩阵R（b）的过程中，会将与被计算的资源块相邻的、具有相近干扰的资源块进行组合，然后利用这个组合中的所有资源块计算该被计算资源块的平均协方差矩阵R（b），而将不相邻的、干扰差异特别大的资源块排除在外，得到更加准确的各个资源块的平均协方差矩阵R（b）。

步骤S102：依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）；

基于步骤S101，在执行步骤S102的过程中，根据步骤S101中计算得到的资源块b的平均协方差矩阵R（b），通过比较资源块b的平均协方差矩阵R（b）中的特定元素，而采用不同的方法计算得到对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）。

需要说明的是，由于步骤S101中对各个资源块的平均协方差矩阵R（b）的计算更加准确，再加上本发明在算法上的改进，能够更加快捷、准确获得对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）。计算而得的各个所述资源块的干扰抑制补偿参数W（b）的取值范围均在[0,1]之间，可以保证数据不会产生溢出。

步骤S103：确定所述接收信号Y中待补偿信号所在资源块，依据所述待补偿信号所在资源块对应的干扰抑制补偿参数W（b）和所述信道响应H进行信号补偿，得到补偿后的接收信号Y’（b）。

接收信号Y中存在有若干个待补偿信号，不同的待补偿信号存在于不同的资源块中。针对每个接收信号Y，首先需要确定其中的待补偿信号，以及该待补偿信号所在的资源块，然后获取根据步骤S101与步骤S102中计算而得的相应资源块的干扰抑制补偿参数W（b）和信道响应H（b）进行信号补偿，得到补偿后的接收信号Y’（b）。

通过上述实施例一中公开的LTE***干扰抑制接收方法，通过根据信道估计H计算各个资源块的平均协方差矩阵R（b），作为当前被计算的资源块的平均协方差矩阵R（b），能够得到更加准确的的各个资源块的平均协方差矩阵R（b），为整个干扰抑制接收的过程奠定基础，保证了干扰抑制接收的准确性以及整个通信网络的性能。

实施例二

请参阅附图2，为本发明实施例二公开的一种LTE***干扰抑制接收方法流程图，主要包括以下步骤：

步骤S101：基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并计算所述接收信号Y所在各个资源块的平均协方差矩阵R（b）；

LTE***采用OFDM技术，传输的资源位置由时间、频率二维资源网格指示，在网格中每个方格表示1个OFDM符号上的1个子载波。

LTE***中定义一个正常子帧的长度为1ms（毫秒），每个正常子帧由2个0.5ms的时隙构成，二维网格则由频域上12个子载波，时域上1个时隙构成并将其定义为资源块。接收信号Y分布在多个资源块中。

如图3给出的基站具备4个公共天线端口的下行公共参考信号在二维网格的资源位置。其中，R0、R1、R2、R3的分布分别对应表示发送天线端口0、1、2、3发送的参考信号的资源位置，黑色粗框中的区域为资源块0，其余资源块可以由此类推，以省略号表示，每个小方格称为资源单元，由坐标（k，l）唯一确定，其中，k表示频率索引，l表示时域符号索引，同时由坐标（k，l）可以计算出该资源单元所在的资源块。图3示出的是正常循环前缀，图4示出的是扩展循环前缀。

在终端侧具备接收天线端口0和1。我们令第j（0≤j<Nr，Nr为接收天线端口数）根接收天线在参考信号的资源位置（k，l）收到的来自第i（0≤i<Nt，Nt为接收天线端口数）根发送天线端口的参考信号为r_ji(k,l),第j根接收天线在非参考信号资源位置（k，l）收到的数据信号为y_j(k.l)，令Y＝{y_j(k,l)|j,k,l}表示接收信号。

在执行步骤S101时，首先，基于本地导频信号，通过使用一些众所周知的信道估计方法，如最小均方误差估计、最小二乘、最大似然或者Kalman滤波，对接收信号Y分布在资源块中的数据信号y_j(k.l)进行信道估计，输出估计信道响应

并令表示信道响应集合，然后计算各个资源块的平均协方差矩阵R（b）。

请参阅附图5，该过程具体由以下几步构成：

步骤S201：基于本地导频信号对接收信号Y进行信道估计，得到信道响应H；

本发明假设接收端具有两根接收天线，因此发送天线i(0,1,2,3)到接收天线j(0,1)的第k子载波第l个符号上的本地导频信号可以表述为

r_ji(k,l)=h_ji(k,l)+m_ji(k,l)

基于本地导频信号，通过使用一些众所周知的信道估计方法，如最小均方误差估计、最小二乘、最大似然或者Kalman滤波，对接收信号Y分布在资源块中的数据信号y_j(k.l)分别进行信道估计，输出估计信道响应

并令 $H=\{{\tilde{h}}_{\mathrm{ji}}(k,l)i,j,k.l\}$ 表示信道响应。

步骤S202：依据所述信道响应H计算各个资源块的协方差矩阵；

根据公式估计每个资源块的协方差矩阵，具体公式为

$\hat{R}\left(b\right)=\frac{1}{N_{\mathrm{tot}}}\underset{i\&Element;S_t}{\mathrm{\&Sigma;}}\underset{k\&Element;S_{k,l}}{\mathrm{\&Sigma;}}\left[\begin{array}{c}{r}_{0i}(k,l)-{\tilde{h}}_{0i}(k,l)\\ r_{1i}(k,l)-{\tilde{h}}_{1i}(k,l)\end{array}\right]{\left[\begin{array}{c}{r}_{0i}(k,l)-{\tilde{h}}_{0i}(k,l)\\ r_{1i}(k,l)-\widetilde{h_{1i}}(k,l)\end{array}\right]}^H$

其中，b表示资源块索引，S_kl为该资源块内的参与计算的本地导频信号集合，S_t为该资源块内的参与计算的发送天线集合，N^tot为所有参与计算的本地导频信号数量。需要说明的是，参与上述公式计算的本地导频信号集合S_k，l可以为整个资源块内所包含的全部，也可以根据预先设置的计算规则选取为其中的一部分。同样，参与上述公式计算的发送天线集合S_t可以为整个资源块内所包含的全部，也可以根据预先设置的计算规则选取为其中的一部分。

步骤S203：对各个所述资源块的协方差矩阵在时域上进行平滑处理；

根据公式对计算出来的每个资源快的协方差矩阵在时域上进行平滑，计算方式为：

其中，α(0<α≤1)为平滑因子，可根据LTE***的实际情况进行配置，既可以是预先设置好的，也可以是变化的。

需要说明的是，该步骤是为了进一步保证最后计算而得的各个资源块的平均协方差矩阵R（b）的准确性，对于应用要求不高的情况下，可以选择不采用该步骤对各个资源块的协方差矩阵进行平滑处理。

步骤S204：确定各个资源块的平均资源块集合；

该步骤可由两种方法实现，具体如下：

方法一：判断各个所述资源块的协方差矩阵与N个相邻资源块的协方差矩阵之间的距离是否都小于预设距离门限；

如果是，则增多N，直到所述资源块的协方差矩阵R与增多后的N个相邻资源块的协方差矩阵R之间的距离不都小于预设距离门限，确定当前资源块与其增多后的N个相邻资源块构成的平均资源块集合；

如果否，则确定由当前资源块与N个相邻资源块构成的平均资源块集合；

具体计算过程为：计算其协方差矩阵与相邻N个资源块的协方差矩阵的距离

其中‖·‖为Frobenius范数。判断上述距离是否小于预设距离门限F，若均满足，则扩大N，直至不能满足，并输出该资源块可平均资源块集合S_b。需要说明的是，上述距离门限F可根据LTE***的实际情况进行配置。

方法二：对每个资源块b，可根据LTE***的实际情况进行配置，确定该资源块可平均资源块集合S_b

需要说明的是，方法一和方法二都能获得该资源块可平均资源块集合Sb。只是方法一更准确，方法二则相对简单，具体可根据应用需求选择使用方法一或方法二。

步骤S205：根据所述平均资源块集合计算各个资源块的平均协方差矩阵R（b）。

根据公式计算每个资源块的平均协方差矩阵R（b），计算方式为：

其中N_b为平均资源块集合S_b中元素的个数。

进一步需要说明的是，基于干扰的频率选择性，在时域上，较短时间内干扰具有连续性的特点，也就是说相邻资源块的干扰情况相近。在计算各个资源块的平均协方差矩阵R（b）的过程中，会将与被计算的资源块相邻的、具有相近干扰的资源块进行组合，然后利用这个组合中的所有资源块计算该被计算资源块的平均协方差矩阵R（b），而将不相邻的、干扰差异特别大的资源块排除在外，得到更加准确的各个资源块的平均协方差矩阵R（b）。

步骤S102：依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）。

基于步骤S101，在执行步骤S102的过程中，根据步骤S101中计算得到的资源块b（b为资源块索引）的平均协方差矩阵R(b)，通过比较资源块b的平均协方差矩阵R(b)中的特定元素，而采用不同的方法计算得到对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）。

具体计算过程如下：

令Rb=R(b)，则所述平均协方差矩阵Rb可表示为 $\mathrm{Rb}=\left[\begin{array}{cc}\mathrm{Rb}\left(\mathrm{1,1}\right)& \mathrm{Rb}\left(\mathrm{1,2}\right)\\ \mathrm{Rb}\left(\mathrm{2,1}\right)& \mathrm{Rb}\left(\mathrm{2,2}\right)\end{array}\right],$

其中，Rb(1,1)表示接收天线0的干扰噪声方差，Rb(1,2)表示接收天线0和接收天线1的干扰噪声协方差，Rb(2,1)表示接收天线1和接收天线0的干扰噪声协方差，Rb(2,2)表示接收天线1的干扰噪声方差；

情况一：当Rb（1,1）大于或等于Rb（2,2）时，依据 $W\left(b\right)=\left(\begin{array}{cc}\mathrm{\&gamma;}& -\mathrm{\&rho;}\\ 0& \sqrt{1-{\left|\mathrm{\&rho;}\right|}^2}\end{array}\right),$ 获取干扰抑制补偿参数W(b)；

其中，γ、ρ为常量， $\mathrm{\&gamma;}=\sqrt{\frac{\mathrm{Rb}\left(\mathrm{2,2}\right)}{\mathrm{Rb}\left(\mathrm{1,1}\right)}},$ $\mathrm{\&rho;}=\frac{\mathrm{\&gamma;}\&CenterDot;\mathrm{Rb}\left(\mathrm{1,2}\right)}{\mathrm{Rb}\left(\mathrm{2,2}\right)}.$

依据NoiseVar(b)＝Rb(2,2)·(1-|ρ|²)得到对应各个所述资源块的等效高斯白噪声的方差NoiseVar(b)。

情况二：当Rb（1,1）小于Rb（2,2）时，令Rb'＝P·Rb·P，其中， $P=\left[\begin{array}{cc}0& 1\\ 1& 0\end{array}\right];$

依据 $V=\left(\begin{array}{cc}\mathrm{\&gamma;}& -\mathrm{\&rho;}\\ 0& \sqrt{1-{\left|\mathrm{\&rho;}\right|}^2}\end{array}\right),$ 获取干扰抑制补偿参数W(b)＝P·V·P；

其中，γ、ρ为常量， $\mathrm{\&gamma;}=\sqrt{\frac{\mathrm{Rb}\left(\mathrm{2,2}\right)}{\mathrm{Rb}\left(\mathrm{1,1}\right)}},$ $\mathrm{\&rho;}=\frac{\mathrm{\&gamma;}\&CenterDot;\mathrm{Rb}\left(\mathrm{1,2}\right)}{\mathrm{Rb}\left(\mathrm{2,2}\right)}.$

依据NoiseVar(b)＝Rb'(2,2)·(1-|ρ|²)得到对应各个所述资源块的等效高斯白噪声的方差NoiseVar(b)。

需要说明的是，由于步骤S101中对各个资源块的平均协方差矩阵R(b)的计算更加准确，再加上本发明在算法上的改进，能够更加快捷、准确获得对应各个所述资源块的干扰抑制补偿参数W(b)和等效高斯白噪声的方差NoiseVar(b)。计算而得的各个所述资源块的干扰抑制补偿参数W（b）的取值范围均在[0,1]之间，可以保证数据不会产生溢出。

步骤S103：确定所述接收信号Y中待补偿信号所在资源块b，依据所述待补偿信号所在资源块对应的干扰抑制补偿参数W（b）和所述信道响应H进行信号补偿，得到补偿后的接收信号Y’（b）。

接收信号Y中存在有若干个待补偿信号，不同的待补偿信号存在于不同的资源块中。针对每个接收信号Y，首先需要确定其中的待补偿信号y_j(k.l)，以及该待补偿信号y_j(k.l)所在的资源块b，然后获取根据步骤S101与步骤S102中计算而得的相应资源块b的干扰抑制补偿参数W(b)和信道响应H进行信号补偿，得到补偿后的待补偿信号y_j(k.l)’，进而得到补偿后的接收信号Y’（b）。

具体计算过程如下：

根据待补偿信号y_j(k.l)所在资源单元(k.l)确定所述待补偿信号所在的资源块b，并选取该资源块的补偿参数W(b)，令

$y(k,l)=\left[\begin{array}{c}{y}_0(k,l)\\ .\\ .\\ .\\ y_{(N_r-1)}(k,l)\end{array}\right]$

实施以下补偿：

$y{(k,l)}^{\&prime;}=\left[\begin{array}{c}{y}_0{(k,l)}^{\&prime;}\\ .\\ .\\ .\\ y_{(N_r-1)}{(k,l)}^{\&prime;}\end{array}\right]=W\left(b\right)\&CenterDot;y(k,l),$

得到补偿后的信道响应

和补偿后的待补偿信号y(k,l)′，待补偿信号补偿完后，与不需要补偿的信号共同组合得到该资源块中经过补偿后的接收信号Y'（b）。

步骤S104：选取所述待补偿信号所在资源块对应的等效高斯白噪声的方差NoiseVar（b），对所述补偿后的接收信号Y'（b）进行均衡处理。

通过前面各个步骤的方法对接收信号Y进行处理以后，接收信号Y'（b）中的干扰噪声已经具备白噪声性质，可以适用于所有的均衡方法，如最小均方误差均衡或最大似然均衡。采用上述均衡方法，基于所述待补偿信号所在资源块对应的等效高斯白噪声的方差NoiseVar（b）对接收信号进行均衡的方式与现有技术相同，这里不再赘述。

需要说明的是，对于所有资源块中的接收信号都是按照以上步骤进行处理，最终得到抑制干扰接收后的接收信号Y'。

通过上述实施例二中公开的LTE***干扰抑制接收方法，通过计算各个资源块的平均协方差矩阵R（b），作为当前被计算的资源块的平均协方差矩阵R（b），能够得到更加准确的的各个资源块的平均协方差矩阵R（b），为整个干扰抑制接收的过程奠定基础；同时，选取更简便的算法计算干扰抑制补偿参数W（b），且得出的干扰抑制补偿参数W（b）取值范围在[0,1]之间，可保证数据不会溢出，进一步保证了干扰抑制接收的准确性以及整个通信网络的性能。

上述本发明公开的实施例中详细描述了方法，对于本发明的方法可采用多种形式的装置实现，因此本发明还公开了一种装置，下面给出具体的实施例进行详细说明。

实施例三

请参阅附图6，为本发明实施例三公开的一种LTE***干扰抑制接收装置具体结构示意图，该装置由协方差计算单元10、干扰抑制参数计算单元11以及预补偿单元12组成，具体链接关系如下：

协方差计算单元10，用于基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并计算所述接收信号Y所在各个资源块的平均协方差矩阵R（b）其中b为资源块索引；

LTE***采用OFDM技术，传输的资源位置由时间、频率二维资源网格指示，在网格中每个方格表示1个OFDM符号上的1个子载波。

LTE***中定义一个正常子帧的长度为1ms（毫秒），每个正常子帧由2个0.5ms的时隙构成，二维网格则由频域上12个子载波，时域上1个时隙构成并将其定义为资源块，接收信号Y分布在多个资源块中。

该协方差计算单元10具体功能如下：

首先，基于本地导频信号，通过使用一些众所周知的信道估计方法，如最小均方误差估计、最小二乘、最大似然或者Kalman滤波，对接收信号Y进行信道估计，输出信道响应H，然后根据信道响应H计算各个资源块的平均协方差矩阵R（b）。

需要说明的是，基于干扰的频率选择性，在时域上，较短时间内干扰具有连续性的特点，也就是说相邻资源块的干扰情况相近。在计算各个资源块的平均协方差矩阵R（b）的过程中，会将与被计算的资源块相邻的、具有相近干扰的资源块进行组合，然后利用这个组合中的所有资源块计算该被计算资源块的平均协方差矩阵R（b），而将不相邻的、干扰差异特别大的资源块排除在外，得到更加准确的各个资源块的平均协方差矩阵R（b）。

干扰抑制参数计算单元11，用于依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）以及等效高斯白噪声的方差NoiseVar（b）；

该干扰抑制参数计算单元11的具体功能如下：

根据协方差计算单元10计算得到的资源块b的平均协方差矩阵R（b），通过比较各个资源块的平均协方差矩阵R（b）中的特定元素，而采用不同的方法计算得到对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）。

需要说明的是，由于协方差计算单元10中对各个资源块的平均协方差矩阵R（b）的计算更加准确，再加上该干扰抑制参数计算单元11采用的算法有所改进，能够更加快捷、准确获得对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）。计算而得的各个所述资源块的干扰抑制补偿参数W的取值范围均在[0,1]之间，可以保证数据不会产生溢出。

预补偿单元12，用于确定所述接收信号Y中待补偿信号所在资源块，依据所述待补偿信号所在资源块对应的干扰抑制补偿参数W（b）和所述信道响应H进行信号补偿，得到补偿后的接收信号Y’（b）。

接收信号Y中存在有若干个待补偿信号，不同的待补偿信号存在于不同的资源块中。针对每个接收信号Y，预补偿单元12首先需要确定其中的待补偿信号，以及该待补偿信号所在的资源块，然后获取根据协方差计算单元10与干扰抑制参数计算单元11中计算而得的相应资源块的干扰抑制补偿参数W（b）和信道响应H进行信号补偿，得到补偿后的接收信号Y’（b）。

通过上述实施例三中公开的LTE***干扰抑制接收装置，通过协方差计算单元计算各个资源块的平均协方差矩阵R（b），作为当前被计算的资源块的平均协方差矩阵R（b），能够得到更加准确的的各个资源块的平均协方差矩阵R（b），为整个干扰抑制接收的过程奠定基础，保证了干扰抑制接收的准确性以及整个通信网络的性能。

实施例四

请参阅附图3，为本发明实施例四公开的一种LTE***干扰抑制接收装置具体结构图，该装置包括：协方差计算单元10、干扰抑制参数计算单元11、预补偿单元12以及均衡单元13，具体连接关系如下：

协方差计算单元10，用于基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并计算所述接收信号Y所在各个资源块的平均协方差矩阵R（b）；

LTE***采用OFDM技术，传输的资源位置由时间、频率二维资源网格指示，在网格中每个方格表示1个OFDM符号上的1个子载波。

LTE***中定义一个正常子帧的长度为1ms（毫秒），每个正常子帧由2个0.5ms的时隙构成，二维网格则由频域上12个子载波，时域上1个时隙构成并将其定义为资源块，接收信号Y分布在多个资源块中。

该协方差计算单元10具体功能如下：

首先，基于本地导频信号，通过使用一些众所周知的信道估计方法，如最小均方误差估计、最小二乘、最大似然或者Kalman滤波，对接收信号Y分布在资源块中的数据信号y_j(k.l)进行信道估计，估计信道响应

并令

表示信道响应集合，然后计算各个资源块的平均协方差矩阵R（b）。

该功能又可分几步完成，请参阅附图8，为本发明实施例四中协方差计算单元10的具体结构，具体如下：

信道估计单元101，用于基于本地导频信号对接收信号Y进行信道估计，得到信道响应H；

本发明假设接收端具有两根接收天线，因此发送天线i(0,1,2,3)到接收天线j(0,1)的第k子载波第1个符号上的本地导频信号可以表述为

r_ji(k)=h_ji(k)+m_ji(k)

基于本地导频信号，通过使用一些众所周知的信道估计方法，如最小均方误差估计、最小二乘、最大似然或者Kalman滤波，对接收信号Y分布在资源块中的数据信号y_j(k.l)分别进行信道估计，输出估计信道响应并令

表示信道响应。

每资源块计算单元102，依据所述信道响应H计算各个资源块的协方差矩阵；

根据公式估计每个资源块的协方差矩阵，具体公式为

$\hat{R}\left(b\right)=\frac{1}{N_{\mathrm{tot}}}\underset{i\&Element;S_t}{\mathrm{\&Sigma;}}\underset{k\&Element;S_{k,l}}{\mathrm{\&Sigma;}}\left[\begin{array}{c}{r}_{0i}(k,l)-{\tilde{h}}_{0i}(k,l)\\ r_{1i}(k,l)-{\tilde{h}}_{1i}(k,l)\end{array}\right]{\left[\begin{array}{c}{r}_{0i}(k,l)-{\tilde{h}}_{0i}(k,l)\\ r_{1i}(k,l)-\widetilde{h_{1i}}(k,l)\end{array}\right]}^H$

其中，b表示资源块索引，S_k为该资源块内的参与计算的本地导频信号集合，S_t为该资源块内的参与计算的发送天线集合，N_tot为所有参与计算的本地导频信号数量。需要说明的是，参与上述公式计算的本地导频信号集合S_k可以为整个资源块内所包含的全部，也可以根据预先设置的计算规则选取为其中的一部分。同样，参与上述公式计算的发送天线集合S_t可以为整个资源块内所包含的全部，也可以根据预先设置的计算规则选取为其中的一部分。

时域平滑单元103，用于对各个所述资源块的协方差矩阵在时域上进行平滑处理；

根据公式对计算出来的每个资源块的协方差矩阵在时域上进行平滑，计算方式为：

其中，α(0<α≤1)为平滑因子，可根据LTE***的实际情况进行配置，既可以是预先设置好的，也可以是变化的。

需要说明的是，该单元是为了进一步保证最后计算而得的各个资源块的平均协方差矩阵R（b）的准确性，对于应用要求不高的情况下，可以选择不采用该单元对各个资源块的协方差矩阵进行平滑处理。

频域参数估计单元104，用于确定各个资源块的平均资源块集合；

判断各个所述资源块的协方差矩阵与N个相邻资源块的协方差矩阵之间的距离是否都小于预设距离门限；

如果是，则增多N，直到所述资源块的协方差矩阵R与增多后的N个相邻资源块的协方差矩阵R之间的距离不都小于预设距离门限，确定当前资源块与其增多后的N个相邻资源块构成的平均资源块集合；

如果否，则确定由当前资源块与N个相邻资源块构成的平均资源块集合；

具体计算过程为：计算其协方差矩阵与相邻N个资源块的协方差矩阵的距离

其中‖·‖为Frobenius范数。判断上述距离是否小于预设距离门限F，若均满足，则扩大N，直至不能满足，并输出该资源块可平均资源块集合S_b。需要说明的是，上述距离门限F可根据LTE***的实际情况进行配置。

需要说明的是，除了采用频域参数估计单元104获得可平均资源块集合S_b，还可以对每个资源块b，可根据LTE***的实际情况进行配置，确定该资源块可平均资源块集合S_b。在这可以通过选择开关K选择使用频域参数估计单元104获得可平均资源块集合S_b或者选择对每个资源块b，可根据LTE***的实际情况进行配置，确定该资源块可平均资源块集合S_b。

频域平均单元105，用于根据所述平均资源块集合计算各个资源块的平均协方差矩阵R（b）。

根据公式计算每个资源块的平均协方差矩阵R（b），计算方式为：

其中N_b为平均资源块集合S_b中元素的个数。

需要说明的是，基于干扰的频率选择性，在时域上，较短时间内干扰具有连续性的特点，也就是说相邻资源块的干扰情况相近。在计算各个资源块的平均协方差矩阵R（b）的过程中，会将与被计算的资源块相邻的、具有相近干扰的资源块进行组合，然后利用这个组合中的所有资源块计算该被计算资源块的平均协方差矩阵R（b），而将不相邻的、干扰差异特别大的资源块排除在外，得到更加准确的各个资源块的平均协方差矩阵R（b）。

干扰抑制参数计算单元11，用于依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）以及等效高斯白噪声的方差NoiseVar（b）；

该干扰抑制参数计算单元11的具体功能如下：

根据协方差计算单元10计算得到的各个资源块b的平均协方差矩阵R（b），通过比较各个资源块b的平均协方差矩阵R（b）中的特定元素，而采用不同的方法计算得到对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）。具体计算过程如下：

令Rb=R(b)，则所述平均协方差矩阵Rb可表示为 $\mathrm{Rb}=\left[\begin{array}{cc}\mathrm{Rb}\left(\mathrm{1,1}\right)& \mathrm{Rb}\left(\mathrm{1,2}\right)\\ \mathrm{Rb}\left(\mathrm{2,1}\right)& \mathrm{Rb}\left(\mathrm{2,2}\right)\end{array}\right],$

其中，Rb(1,1)表示接收天线0的干扰噪声方差，Rb(1,2)表示接收天线0和接收天线1的干扰噪声协方差，Rb(2,1)表示接收天线1和接收天线0的干扰噪声协方差，Rb(2,2)表示接收天线1的干扰噪声方差；

情况一：当Rb（1,1）大于或等于Rb（2,2）时，依据 $W\left(b\right)=\left(\begin{array}{cc}\mathrm{\&gamma;}& -\mathrm{\&rho;}\\ 0& \sqrt{1-{\left|\mathrm{\&rho;}\right|}^2}\end{array}\right),$ 获取干扰抑制补偿参数W(b)；

其中，γ、ρ为常量， $\mathrm{\&gamma;}=\sqrt{\frac{\mathrm{Rb}\left(\mathrm{2,2}\right)}{\mathrm{Rb}\left(\mathrm{1,1}\right)}},$ $\mathrm{\&rho;}=\frac{\mathrm{\&gamma;}\&CenterDot;\mathrm{Rb}\left(\mathrm{1,2}\right)}{\mathrm{Rb}\left(\mathrm{2,2}\right)}.$

依据NoiseVar(b)＝Rb(2,2)·(1-|ρ|²)得到对应各个所述资源块的等效高斯白噪声的方差NoiseVar(b)。

情况二：当Rb（1,1）小于Rb（2,2）时，令Rb'＝P·Rb·P，其中， $P=\left[\begin{array}{cc}0& 1\\ 1& 0\end{array}\right];$

依据 $V=\left(\begin{array}{cc}\mathrm{\&gamma;}& -\mathrm{\&rho;}\\ 0& \sqrt{1-{\left|\mathrm{\&rho;}\right|}^2}\end{array}\right),$ 获取干扰抑制补偿参数W(b)＝P·V·P；

其中，γ、ρ为常量， $\mathrm{\&gamma;}=\sqrt{\frac{\mathrm{Rb}\left(\mathrm{2,2}\right)}{\mathrm{Rb}\left(\mathrm{1,1}\right)}},$ $\mathrm{\&rho;}=\frac{\mathrm{\&gamma;}\&CenterDot;\mathrm{Rb}\left(\mathrm{1,2}\right)}{\mathrm{Rb}\left(\mathrm{2,2}\right)}.$

依据NoiseVar(b)＝Rb'(2,2)·(1-|ρ|²)得到对应各个所述资源块的等效高斯白噪声的方差NoiseVar(b)。

需要说明的是，由于协方差计算单元10中对各个资源块的平均协方差矩阵R（b）的计算更加准确，再加上该干扰抑制参数计算单元11采用的算法有所改进，能够更加快捷、准确获得对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）。计算而得的各个所述资源块的干扰抑制补偿参数W（b）的取值范围均在[0,1]之间，可以保证数据不会产生溢出。

预补偿单元12，用于确定所述接收信号Y中待补偿信号所在资源块b，依据所述待补偿信号所在资源块对应的干扰抑制补偿参数W（b）和所述信道响应H进行信号补偿，得到补偿后的接收信号Y’（b）。

接收信号Y中存在有若干个待补偿信号，不同的待补偿信号存在于不同的资源块中。针对每个接收信号Y，预补偿单元12首先需要确定其中的待补偿信号y_j(k.l)，以及该待补偿信号y_j(k.l)所在的资源块b，然后获取根据协方差计算单元10与干扰抑制参数计算单元11中计算而得的相应资源块的干扰抑制补偿参数W（b）和信道响应H（b）进行信号补偿，得到补偿后的待补偿信号y_j(k.l)’，进而得到补偿后的接收信号Y’（b）。

具体计算过程如下：

根据待补偿信号y_j(k.l)所在资源单元(k.l)确定所述待补偿信号所在的资源块b，并选取该资源块的补偿参数W(b)，令

$y(k,l)=\left[\begin{array}{c}{y}_0(k,l)\\ .\\ .\\ .\\ y_{(N_r-1)}(k,l)\end{array}\right]$

实施以下补偿：

$y{(k,l)}^{\&prime;}=\left[\begin{array}{c}{y}_0{(k,l)}^{\&prime;}\\ .\\ .\\ .\\ y_{(N_r-1)}{(k,l)}^{\&prime;}\end{array}\right]=W\left(b\right)\&CenterDot;y(k,l),$

得到补偿后的信道响应

和补偿后的待补偿信号y(k,l)′，待补偿信号补偿完后，与不需要补偿的信号共同组合得到该资源块中经过补偿后的接收信号Y'（b）。

均衡单元13，用于选取所述待补偿信号所在资源块对应的等效高斯白噪声的方差NoiseVar（b），对所述补偿后的接收信号Y'（b）进行均衡处理。

通过前面各个单元对接收信号Y进行处理以后，接收信号Y'（b）中的干扰噪声已经具备白噪声性质，可以适用于所有的均衡装置，如最小均方误差均衡器或最大似然均衡器。采用上述均衡装置，基于所述待补偿信号所在资源块对应的等效高斯白噪声的方差NoiseVar（b）对接收信号进行均衡的方式与现有技术相同，这里不再赘述。

需要说明的是，对于所有资源块中的接收信号都是按照以上步骤进行处理，最终得到抑制干扰接收后的接收信号Y'。

通过上述实施例四中公开的LTE***干扰抑制接收装置，通过协方差计算单元10各个资源块的平均协方差矩阵R（b），作为当前被计算的资源块的平均协方差矩阵R（b），能够得到更加准确的的各个资源块的平均协方差矩阵R（b），为整个干扰抑制接收的过程奠定基础；同时，干扰抑制参数计算单元11选取更简便的算法计算干扰抑制补偿参数W（b），且得出的干扰抑制补偿参数W（b）取值范围在[0,1]之间，可保证数据不会溢出，进一步保证了干扰抑制接收的准确性以及整个通信网络的性能。

综上所述：

本发明公开了一种LTE***干扰抑制接收方法及装置。基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并根据信道响应H计算所述各个资源块的平均协方差矩阵R（b）；依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）；确定所述接收信号Y中待补偿信号所在资源块b，依据所述待补偿信号所在资源块对应的干扰抑制补偿参数W（b）和所述信道响应H进行信号补偿，得到补偿后的接收信号Y’。基于上述方法及装置，能够获得更准确的平均协方差矩阵R（b），保证了干扰抑制接收的准确性以及整个通信网络的性能。

本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。

结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储器（RAM）、内存、只读存储器（ROM）、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。

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

Claims (14)

1.一种LTE（Long Term Evolution system，长期演进***）***干扰抑制接收方法，其特征在于，包括：

基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并计算所述接收信号Y所在各个资源块的平均协方差矩阵R（b）；

依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）；

确定所述接收信号Y中待补偿信号所在资源块，依据所述待补偿信号所在资源块对应的干扰抑制补偿参数W（b）和所述信道响应H进行信号补偿，得到补偿后的接收信号Y’（b）。

2.根据权利要求1所述的方法，其特征在于，在得到补偿后的接收信号之后，还包括：

选取所述待补偿信号所在资源块对应的等效高斯白噪声的方差NoiseVar（b），对所述补偿后的接收信号Y'（b）进行均衡处理。

3.根据权利要求1或2所述的方法，其特征在于，基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并计算所述接收信号Y所在各个资源块的平均协方差矩阵R（b）的过程包括：

基于本地导频信号对接收信号Y进行信道估计，得到信道响应H；

依据所述信道响应H计算各个资源块的协方差矩阵；

判断各个所述资源块的协方差矩阵与N个相邻资源块的协方差矩阵之间的距离是否都小于预设距离门限；

如果是，则增多N，直到所述资源块的协方差矩阵与增多后的N个相邻资源块的协方差矩阵之间的距离不都小于预设距离门限，确定当前资源块与其增多后的N个相邻资源块构成的平均资源块集合；

如果否，则确定由当前资源块与N个相邻资源块构成的平均资源块集合；

根据所述平均资源块集合计算各个资源块平均协方差矩阵R（b）。

4.根据权利要求3所述的方法，其特征在于，在所述判断各个所述资源块的协方差矩阵与N个相邻资源块之间的距离是否小于预设距离门限之前，所述计算所述估计信道响应H中各个资源块的协方差矩阵之后，还包括：

对各个所述资源块的协方差矩阵在时域上进行平滑处理。

5.根据权利要求1或2所述的方法，其特征在于，基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并计算所述接收信号Y所在各个资源块的平均协方差矩阵R（b）的过程包括：

基于本地导频信号对接收信号Y进行信道估计，得到信道响应H；

依据所述信道响应H计算各个资源块的协方差矩阵；

根据预设平均资源块集合计算各个资源块平均协方差矩阵R（b）。

6.根据权利要求1所述的方法，其特征在于，依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）的过程包括：

将所述平均协方差矩阵R（b）转换为 $\mathrm{Rb}=\left[\begin{array}{cc}\mathrm{Rb}\left(\mathrm{1,1}\right)& \mathrm{Rb}\left(\mathrm{1,2}\right)\\ \mathrm{Rb}\left(\mathrm{2,1}\right)& \mathrm{Rb}\left(\mathrm{2,2}\right)\end{array}\right],$

其中，Rb(1,1)表示接收天线0的干扰噪声方差，Rb(1,2)表示接收天线0和接收天线1的干扰噪声协方差，Rb(2,1)表示接收天线1和接收天线0的干扰噪声协方差，Rb(2,2)表示接收天线1的干扰噪声方差；

当Rb（1,1）大于或等于Rb（2,2）时，依据 $W\left(b\right)=\left(\begin{array}{cc}\mathrm{\&gamma;}& -\mathrm{\&rho;}\\ 0& \sqrt{1-{\left|\mathrm{\&rho;}\right|}^2}\end{array}\right),$ 获取干扰抑制补偿参数W（b）；

其中，γ、ρ为常量， $\mathrm{\&gamma;}=\sqrt{\frac{\mathrm{Rb}\left(\mathrm{2,2}\right)}{\mathrm{Rb}\left(\mathrm{1,1}\right)}},$ $\mathrm{\&rho;}=\frac{\mathrm{\&gamma;}\&CenterDot;\mathrm{Rb}\left(\mathrm{1,2}\right)}{\mathrm{Rb}\left(\mathrm{2,2}\right)}.$

7.根据权利要求1所述的方法，其特征在于，依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）的过程包括：

将所述平均协方差矩阵R（b）转换为 $\mathrm{Rb}=\left[\begin{array}{cc}\mathrm{Rb}\left(\mathrm{1,1}\right)& \mathrm{Rb}\left(\mathrm{1,2}\right)\\ \mathrm{Rb}\left(\mathrm{2,1}\right)& \mathrm{Rb}\left(\mathrm{2,2}\right)\end{array}\right],$

其中，Rb(1,1)表示接收天线0的干扰噪声方差，Rb(1,2)表示接收天线0和接收天线1的干扰噪声协方差，Rb(2,1)表示接收天线1和接收天线0的干扰噪声协方差，Rb(2,2)表示接收天线1的干扰噪声方差；

当Rb（1,1）小于Rb（2,2）时，令R'=P·R·P，其中， $P=\left[\begin{array}{cc}0& 1\\ 1& 0\end{array}\right];$

依据 $V=\left(\begin{array}{cc}\mathrm{\&gamma;}& -\mathrm{\&rho;}\\ 0& \sqrt{1-{\left|\mathrm{\&rho;}\right|}^2}\end{array}\right),$ 获取干扰抑制补偿参数W(b)＝P·V·P；

其中，γ、ρ为常量， $\mathrm{\&gamma;}=\sqrt{\frac{\mathrm{Rb}\left(\mathrm{2,2}\right)}{\mathrm{Rb}\left(\mathrm{1,1}\right)}},$ $\mathrm{\&rho;}=\frac{\mathrm{\&gamma;}\&CenterDot;\mathrm{Rb}\left(\mathrm{1,2}\right)}{\mathrm{Rb}\left(\mathrm{2,2}\right)}.$

8.根据权利要求6所述的方法，其特征在于，依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的等效高斯白噪声的方差NoiseVar（b）的过程包括：

依据NoiseVar(b)＝Rb(2,2)·(1-|ρ|²)得到对应各个所述资源块的等效高斯白噪声的方差NoiseVar（b）。

9.根据权利要求7所述的方法，其特征在于，依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的等效高斯白噪声的方差NoiseVar（b）包括：

依据NoiseVar(b)＝Rb'(2,2)·(1-|ρ|²)得到对应各个所述资源块的等效高斯白噪声的方差NoiseVar（b）。

10.一种LTE***干扰抑制接收装置，其特征在于，包括：

协方差计算单元，用于基于本地导频信号对接收信号Y进行信道估计，得到信道响应H，并计算所述接收信号Y所在各个资源块的平均协方差矩阵R（b）；

干扰抑制参数计算单元，用于依据各个所述资源块的所述平均协方差矩阵R（b）进行计算，得到对应各个所述资源块的干扰抑制补偿参数W（b）和等效高斯白噪声的方差NoiseVar（b）；

预补偿单元，用于确定所述接收信号Y中待补偿信号所在资源块，依据所述待补偿信号所在资源块对应的干扰抑制补偿参数W（b）和所述信道响应H进行信号补偿，得到补偿后的接收信号Y’（b）。

11.根据权利要求10所述的装置，其特征在于，还包括：

均衡单元，用于选取所述待补偿信号所在资源块对应的等效高斯白噪声的方差NoiseVar（b），对所述补偿后的接收信号Y'（b）进行均衡处理。

12.根据权利要求10或11所述的装置，其特征在于，所述协方差计算单元包括：

信道估计单元，用于基于本地导频信号对接收信号Y进行信道估计，得到信道响应H；

每资源块计算单元，用于依据所述信道响应H计算各个资源块的协方差矩阵；

频域参数估计单元，用于判断各个所述资源块的协方差矩阵与N个相邻资源块的协方差矩阵之间的距离是否都小于预设距离门限；

频域平均单元，用于根据所述平均资源块集合计算各个资源块平均协方差矩阵R（b）。

13.根据权利要求12所述的装置，其特征在于，还包括：

时域平滑单元，用于对各个所述资源块的协方差矩阵在时域上进行平滑处理。

14.根据权利要求12所述的装置，其特征在于，还包括：

选择开关K，用于选择频域参数估计单元确定可平均资源块集合或选择预设平均资源块集合。

Images