CN104243377A - 一种干扰抑制方法和装置 - Google Patents

Abstract

本发明公开了一种干扰抑制方法和装置，用于解决由于无法精确估计PMI，导致不能使用SLIC或者ML接收机消除干扰的问题。方法包括：获取信道参数，包括服务小区对应的第一信道矩阵和干扰小区对应的第二信道矩阵；对接收信号进行等效变换，得到其中，y为接收信号，x₀为服务小区发送的第一发射信号，x_I为干扰小区发送的第二发射信号，n为噪声信号，H₀为第一信道矩阵与服务小区对应的预编码矩阵的乘积，为第二信道矩阵，P_I为干扰小区对应的预编码矩阵；将P_Ix_I的乘积作为等效发射信号，对等效变换后的接收信号进行解调处理，以在不获取所述干扰小区对应的预编码矩阵的系数的前提下，消除所述等效发射信号的干扰。

Description

一种干扰抑制方法和装置

技术领域

本发明涉及通信技术领域，特别涉及一种干扰抑制方法和装置。

背景技术

用户对个人通信的需求不断增长，推动着无线传输技术不断向前演进。2004年底，第三代合作伙伴计划(3rd Generation Partnership Project，3GPP)开始了通用移动通信系统(Universal Mobile Telecommunications System，UMTS)技术的长期演进(Long Term Evolution，LTE)项目。在LTE实际场景中，用户设备(User Equipment，UE)可能受到邻小区的干扰，或受到本小区内配对用户的干扰，这些干扰将严重降低UE解调数据信道的性能。基站高密度和异构是LTE网络结构的演进方向，这显然会带来更加严重的小区间和多用户间的同频干扰。

在LTE版本Release 11阶段定义了抑制干扰的标准接收机，即如干扰抑制合并(Interference Rejection Combining，IRC)接收机。然而，IRC抑制小区间干扰的能力有限，且完全不能抑制用户间干扰。因此，在LTE Release 12阶段定义了增强的干扰抑制接收机，即符号级干扰消除(Symbol Level InterferenceCancellation，SLIC)接收机和最大似然(Maximum Likelihood，ML)接收机。这两类接收机具有更好的抑制用户间和小区间干扰的能力。但是，这两类接收机都需要已知干扰信号的某些参数，如调制阶数、传输模式和预编码矩阵指示(Pre-coder Matrix Indicator，PMI)等，而这些参数需要UE自主估计得到。

LTE多输入多输出(Multiple-Input Multiple-Output，MIMO)-正交频分复用(Orthogonal Frequency Division Multiplex，OFDM)系统中，假设UE侧接收天线个数为N，用x、y和H分别表示频域某个资源单元(Resource Element,RE)的发射信号、接收信号和等效频域信道矩阵(包含PMI所指示的预编码矩阵)，那么信道模型可以表示为：

y＝Hx+n  (1)

其中，假设n是N维加性高斯白噪声(均值为0，协方差矩阵为σ2I)；x为L维的列向量(L表示发射信号总的层数(即发射信号的秩)，包括服务小区和所有待消除的干扰小区)，即x＝[x₁,x₂,…,x_L]^T，信道矩阵H＝[h₁,h₂,…,h_L](h_k为N维列向量，表示第k个发射符号对应的等效频域信道)。

为了区分多个小区的信号和信道，可以把x和H拆成两部分，即：

y＝H₀x₀+H_Ix_I+n  (2)

其中，0下标表示服务小区，I下标表示干扰小区，不失一般性，假设前B层发射符号来源于服务小区，从B+1层到第L层发射符号来源于其它干扰小区(可以是一个干扰小区，也可以是多个干扰小区)或同小区配对的其他UE，即H₀＝[h₁,h₂,…,h_B]，H_I＝[h_B+1,h_B+2,…,h_L]，x₀＝[x₁,x₂,…,x_B]^T和x_I＝[x_B+1,x_B+2,…,x_L]^T。

通常采用能量算法或者广义最大似然算法，估计干扰的PMI参数。

1、能量算法：

$P_{\mathrm{opt}}=\arg \underset{P_i}{\min }\mathrm{\Σ}{\left|\right|R_{\mathrm{yy}}-H_0{H}_0^H-{\tilde{H}}_I{P}_I{P}_I^H{\tilde{H}}_I^H-{\mathrm{\σ}}^{2}I\left|\right|}^2---\left(3a\right)$

R_yy＝E(yy^H)  (3b)

其中，P_I是干扰小区对应的预编码矩阵，σ是噪声功率，(3a)中的求和在所有RE样点上进行，(3b)中数学期望可以用样点的代数平均获得。

2、广义最大似然算法：

$P_{\mathrm{opt}}=\arg \underset{P_i}{\min }\mathrm{\Σ}\{\ln \det (H_0{H}_0^H+{\tilde{H}}_I{P}_I{P}_I^H{\tilde{H}}_I^H+{\mathrm{\σ}}^{2}I)+y^H{(H_0{H}_0^H+{\tilde{H}}_I{P}_I{P}_I^H{\tilde{H}}_I^H+{\mathrm{\σ}}^{2}I)}^{-1}y\};$

对于传输模式(Transmission Mode，TM)3和TM4，当干扰信号使用2端口且层数为2，或者使用4端口且层数为4时，无论干扰小区对应的PMI取哪个值，都是单位矩阵。因此，无论是能量方法还是广义最大似然方法都是无法精确估计干扰小区对应的PMI的。由于上述场景下，无法精确估计干扰小区对应的PMI，因此，不能使用SLIC或者ML接收机消除干扰，从而限制了UE接收机抑制干扰的能力。

发明内容

本发明实施例提供了一种干扰抑制方法和装置，用于解决在干扰信号使用2端口且层数为2，或者使用4端口且层数为4场景下，由于无法精确估计干扰小区对应的PMI，而导致不能使用SLIC或者ML接收机消除干扰的问题。

第一方面，一种干扰抑制方法，该方法包括：

获取信道参数，所述信道参数包括服务小区对应的第一信道矩阵和干扰小区对应的第二信道矩阵；

对接收信号进行等效变换，得到其中，y为接收信号，x₀为所述服务小区发送的第一发射信号，x_I为所述干扰小区发送的第二发射信号，n为噪声信号，H₀为所述第一信道矩阵与所述服务小区对应的预编码矩阵的乘积，为所述第二信道矩阵，P_I为所述干扰小区对应的预编码矩阵；

将P_Ix_I的乘积作为等效发射信号，对等效变换后的接收信号进行解调处理，以在不获取所述干扰小区对应的预编码矩阵的系数的前提下，消除所述等效发射信号的干扰。

结合第一方面，在第一种可能的实现方式中，对等效变换后的接收信号进行解调处理之前，还包括：根据所述第二发射信号采用的调制方式，确定出所述等效发射信号对应的星座图；

根据以下步骤确定解调处理中使用的解调方式：

根据所述等效发射信号对应的星座图中的星座点的个数，确定解调处理中使用的解调方式；

若确定出的星座点的个数大于设定阈值，确定采用第一解调方式对所述发射信号进行解调处理，其中，所述第一解调方式包括最小均方差MMSE解调；或者，若确定出的星座点的个数小于或等于所述阈值，确定采用第二解调方式对所述发射信号进行解调处理，其中，所述第二解调方式包括符号级干扰消除SLIC解调和最大似然ML解调。

结合第一方面、或第一方面的第一种可能的实现方式，在第二种可能的实现方式中，若采用SLIC解调，对等效变换后的接收信号进行解调处理，包括：

根据如下公式，分别确定每次解调中所述等效发射信号的每层信号的均值和方差；

根据得到的均值和方差确定出每次解调中所述等效发射信号的每层信号的对数似然比LLR值；

$E\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{x}_j\mathrm{Pr}\left(x_j\right);$

$\mathrm{Var}\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{|x_j-E\left(x_j\right)|}^2\mathrm{Pr}\left(x_j\right);$

其中，E(x_j)表示所述等效发射信号的第j层信号的均值，Var(x_j)表示所述等效发射信号的第j层信号的方差，Pr(x_j)表示E(x_j)的先验概率，{x_j}表示所述等效发射信号的第j层信号在所述等效发射信号对应的星座图中所有可能的取值点，表示所述等效发射信号的第j层信号的加权向量，表示所述等效发射信号的第j层信号的等效噪声，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

结合第一方面、或第一方面的第一种可能的实现方式，在第三种可能的实现方式中，若采用SLIC解调，对等效变换后的接收信号进行解调处理，包括：

根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组，以在计算每组星座点的LLR、均值和方差时能够复用所述第二发射信号对应的星座图中的星座点的LLR、均值和方差；

针对所述SLIC解调中的每一次解调，根据上一次解调的LLR值，分别计算出每组星座点的均值和方差；将计算出的每组星座点的均值进行合并得到该次解调的均值，将计算出的每组星座点的方差进行合并得到该次解调的方差；根据该次解调的均值和方差，计算该次解调的LLR值。

结合第一方面的第三种可能的实现方式，在第四种可能的实现方式中，按照以下公式，将计算出的每组星座点的均值进行合并得到所述等效发射信号的每层信号在该次SLIC解调的均值：

$E\left(x_j\right)={\mathrm{\Σ}}_{m=1}^{M}E\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right);$

其中，E(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的均值，E(Ω_m)和Pr(Ω_m)分别表示第m组星座点的均值和第m组星座点的均值对应的先验概率，M表示将所述等效发射信号对应的星座图中的星座点划分为M组，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

结合第一方面的第三种可能的实现方式，在第五种可能的实现方式中，按照以下公式，将计算出的每组星座点的方差进行合并得到所述等效发射信号中的每层信号在该次SLIC解调的方差：

$\mathrm{Var}\left(x_j\right)={\mathrm{\Σ}}_{m=1}^M\mathrm{Var}\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)+{\mathrm{\Σ}}_{m=1}^M{E}^2\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)-E\left(x_j\right);$

其中，Var(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的方差，Var(Ω_m)表示第m组星座点的方差，Pr(Ω_m)表示第m组星座点的均值对应的先验概率，M表示将所述等效发射信号对应的星座图中的星座点划分为M组，j＝B+1,……,L，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

结合第一方面、或第一方面的第一种可能的实现方式，在第六种可能的实现方式中，若采用ML解调，对等效变换后的接收信号进行解调处理，包括：

根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组；

根据所述第二发射信号进行ML解调时所采用的节点选择方法，从每组星座点中，分别选择出设定数量的节点；

从已选择出的节点中，分别为所述等效发射信号的每层信号中的每个比特信号选择最终节点；以及

根据选择出的最终节点，计算所述等效发射信号的每层信号中的每个比特信号的LLR值。

结合第一方面的第六种可能的实现方式，在第七种可能的实现方式中，从已选择出的节点中，为所述等效发射信号的每层信号中的每个比特信号选择最终节点，包括：

从已选择出的节点中，选择出满足设定条件的节点，并将选择出的节点确定为所述等效发射信号的每层信号中的每个比特信号选择最终节点；

其中，设定条件为 $\left\{\begin{array}{c}\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\\ \underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\end{array}\right.,$ {x}表示所述等效发射信号在其对应的星座图中所有可能的取值点，b_j,i表示第j层第i个比特信号，σ²表示噪声功率，表示所述等效发射信号的均值对应的先验概率，y表示接收数据，H表示信道矩阵，x表示所述第二发射信号，j＝B+1,……,L，L表示所述发射信号的总层数，B表示所述第一发射信号的总层数。

结合第一方面的第七种可能的实现方式，在第八种可能的实现方式中，根据选择出的最终节点，按照以下公式，计算所述等效发射信号的每层信号中的每个比特信号的LLR值：

${\mathrm{\λ}}_{j,i}=\frac{1}{{\mathrm{\σ}}^2}\{\underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))-\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\};$

其中，λ_j,i表示所述等效发射信号的第j层的第i个比特信号的LLR值。

结合第一方面的第三种可能的实现方式、或第一方面的第六种可能的实现方式，在第九种可能的实现方式中，根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组，包括：

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过旋转相同角度得到星座点划分为一组；或者

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过压缩相同压缩量得到星座点划分为一组；或者

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过旋转相同角度并压缩相同压缩量得到星座点划分为一组。

第二方面，一种干扰抑制装置，该装置包括：

获取模块，用于获取信道参数，所述信道参数包括服务小区对应的第一信道矩阵和干扰小区对应的第二信道矩阵；

等效变换模块，用于对接收信号进行等效变换，得到其中，y为接收信号，x₀为所述服务小区发送的第一发射信号，x_I为所述干扰小区发送的第二发射信号，n为噪声信号，H₀为所述第一信道矩阵与所述服务小区对应的预编码矩阵的乘积，为所述第二信道矩阵，P_I为所述干扰小区对应的预编码矩阵；

解调模块，用于将P_Ix_I的乘积作为等效发射信号，对等效变换后的接收信号进行解调处理，以在不获取所述干扰小区对应的预编码矩阵的系数的前提下，消除所述等效发射信号的干扰。

结合第二方面，在第一种可能的实现方式中，所述解调模块对等效变换后的接收信号进行解调处理之前，还用于：根据所述第二发射信号采用的调制方式，确定出所述等效发射信号对应的星座图；

所述解调模块根据以下步骤确定解调处理中使用的解调方式：

根据所述等效发射信号对应的星座图中的星座点的个数，确定解调处理中使用的解调方式；

若确定出的星座点的个数大于设定阈值，确定采用第一解调方式对所述发射信号进行解调处理，其中，所述第一解调方式包括最小均方差MMSE解调；或者，若确定出的星座点的个数小于或等于所述阈值，确定采用第二解调方式对所述发射信号进行解调处理，其中，所述第二解调方式包括符号级干扰消除SLIC解调和最大似然ML解调。

结合第二方面、或第二方面的第一种可能的实现方式，在第二种可能的实现方式中，若采用SLIC解调，所述解调模块对等效变换后的接收信号进行解调处理，包括：

根据如下公式，分别确定每次解调中所述等效发射信号的每层信号的均值和方差；以及根据得到的均值和方差确定出每次解调中所述等效发射信号的每层信号的对数似然比LLR值；

$E\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{x}_j\mathrm{Pr}\left(x_j\right);$

$\mathrm{Var}\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{|x_j-E\left(x_j\right)|}^2\mathrm{Pr}\left(x_j\right);$

其中，E(x_j)表示所述等效发射信号的第j层信号的均值，Var(x_j)表示所述等效发射信号的第j层信号的方差，Pr(x_j)表示E(x_j)的先验概率，{x_j}表示所述等效发射信号的第j层信号在所述等效发射信号对应的星座图中所有可能的取值点，表示所述等效发射信号的第j层信号的加权向量，表示所述等效发射信号的第j层信号的等效噪声，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

结合第二方面、或第二方面的第一种可能的实现方式，在第三种可能的实现方式中，若采用SLIC解调，所述解调模块对等效变换后的接收信号进行解调处理，包括：

根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组，以在计算每组星座点的LLR、均值和方差时能够复用所述第二发射信号对应的星座图中的星座点的LLR、均值和方差；

针对所述SLIC解调中的每一次解调，根据上一次解调的LLR值，分别计算出每组星座点的均值和方差；将计算出的每组星座点的均值进行合并得到该次解调的均值，将计算出的每组星座点的方差进行合并得到该次解调的方差；根据该次解调的均值和方差，计算该次解调的LLR值。

结合第二方面的第三种可能的实现方式，在第四种可能的实现方式中，所述解调模块按照以下公式，将计算出的每组星座点的均值进行合并得到所述等效发射信号的每层信号在该次SLIC解调的均值：

$E\left(x_j\right)={\mathrm{\Σ}}_{m=1}^{M}E\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right);$

其中，E(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的均值，E(Ω_m)和Pr(Ω_m)分别表示第m组星座点的均值和第m组星座点的均值对应的先验概率，M表示将所述等效发射信号对应的星座图中的星座点划分为M组，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

结合第二方面的第三种可能的实现方式，在第五种可能的实现方式中，所述解调模块按照以下公式，将计算出的每组星座点的方差进行合并得到所述等效发射信号中的每层信号在该次SLIC解调的方差：

$\mathrm{Var}\left(x_j\right)={\mathrm{\Σ}}_{m=1}^M\mathrm{Var}\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)+{\mathrm{\Σ}}_{m=1}^M{E}^2\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)-E\left(x_j\right);$

其中，Var(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的方差，Var(Ω_m)表示第m组星座点的方差，Pr(Ω_m)表示第m组星座点的均值对应的先验概率，M表示将所述等效发射信号对应的星座图中的星座点划分为M组，j＝B+1,……,L，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

结合第二方面、或第二方面的第一种可能的实现方式，在第六种可能的实现方式中，若采用ML解调，所述解调模块对等效变换后的接收信号进行解调处理，包括：

根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组；根据所述第二发射信号进行ML解调时所采用的节点选择方法，从每组星座点中，分别选择出设定数量的节点；从已选择出的节点中，分别为所述等效发射信号的每层信号中的每个比特信号选择最终节点；以及根据选择出的最终节点，计算所述等效发射信号的每层信号中的每个比特信号的LLR值。

结合第二方面的第六种可能的实现方式，在第七种可能的实现方式中，所述解调模块从已选择出的节点中，为所述等效发射信号的每层信号中的每个比特信号选择最终节点，包括：

从已选择出的节点中，选择出满足设定条件的节点，并将选择出的节点确定为所述等效发射信号的每层信号中的每个比特信号选择最终节点；

其中，设定条件为 $\left\{\begin{array}{c}\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\\ \underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\end{array}\right.,$ {x}表示所述等效发射信号在其对应的星座图中所有可能的取值点，b_j,i表示第j层第i个比特信号，σ²表示噪声功率，表示所述等效发射信号的均值对应的先验概率，y表示接收数据，H表示信道矩阵，x表示所述第二发射信号，j＝B+1,……,L，L表示所述发射信号的总层数，B表示所述第一发射信号的总层数。

结合第二方面的第七种可能的实现方式，在第八种可能的实现方式中，所述解调模块根据选择出的最终节点，按照以下公式，计算所述等效发射信号的每层信号中的每个比特信号的LLR值：

${\mathrm{\λ}}_{j,i}=\frac{1}{{\mathrm{\σ}}^2}\{\underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))-\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\};$

其中，λ_j,i表示所述等效发射信号的第j层的第i个比特信号的LLR值。

结合第二方面的第三种可能的实现方式、或第二方面的第六种可能的实现方式，在第九种可能的实现方式中，所述解调模块根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组，包括：

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过旋转相同角度得到星座点划分为一组；或者

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过压缩相同压缩量得到星座点划分为一组；或者

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过旋转相同角度并压缩相同压缩量得到星座点划分为一组。

第三方面，一种通信设备，该通信设备包括：

处理器，用于获取信道参数，所述信道参数包括服务小区对应的第一信道矩阵和干扰小区对应的第二信道矩阵；

接收机，用于对接收信号进行等效变换，得到其中，y为接收信号，x₀为所述服务小区发送的第一发射信号，x_I为所述干扰小区发送的第二发射信号，n为噪声信号，H₀为所述第一信道矩阵与所述服务小区对应的预编码矩阵的乘积，为所述第二信道矩阵，P_I为所述干扰小区对应的预编码矩阵；将P_Ix_I的乘积作为等效发射信号，对等效变换后的接收信号进行解调处理，以在不获取所述干扰小区对应的预编码矩阵的系数的前提下，消除所述等效发射信号的干扰。

结合第三方面，在第一种可能的实现方式中，所述接收机对等效变换后的接收信号进行解调处理之前，还用于：根据所述第二发射信号采用的调制方式，确定出所述等效发射信号对应的星座图；

所述接收机根据以下步骤确定解调处理中使用的解调方式：

根据所述等效发射信号对应的星座图中的星座点的个数，确定解调处理中使用的解调方式；

若确定出的星座点的个数大于设定阈值，确定采用第一解调方式对所述发射信号进行解调处理，其中，所述第一解调方式包括最小均方差MMSE解调；或者，若确定出的星座点的个数小于或等于所述阈值，确定采用第二解调方式对所述发射信号进行解调处理，其中，所述第二解调方式包括符号级干扰消除SLIC解调和最大似然ML解调。

结合第三方面、或第三方面的第一种可能的实现方式，在第二种可能的实现方式中，若采用SLIC解调，所述接收机对等效变换后的接收信号进行解调处理，包括：

根据如下公式，分别确定每次解调中所述等效发射信号的每层信号的均值和方差；以及根据得到的均值和方差确定出每次解调中所述等效发射信号的每层信号的对数似然比LLR值；

$E\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{x}_j\mathrm{Pr}\left(x_j\right);$

$\mathrm{Var}\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{|x_j-E\left(x_j\right)|}^2\mathrm{Pr}\left(x_j\right);$

其中，E(x_j)表示所述等效发射信号的第j层信号的均值，Var(x_j)表示所述等效发射信号的第j层信号的方差，Pr(x_j)表示E(x_j)的先验概率，{x_j}表示所述等效发射信号的第j层信号在所述等效发射信号对应的星座图中所有可能的取值点，表示所述等效发射信号的第j层信号的加权向量，表示所述等效发射信号的第j层信号的等效噪声，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

结合第三方面、或第三方面的第一种可能的实现方式，在第三种可能的实现方式中，若采用SLIC解调，所述接收机对等效变换后的接收信号进行解调处理，包括：

根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组，以在计算每组星座点的LLR、均值和方差时能够复用所述第二发射信号对应的星座图中的星座点的LLR、均值和方差；

针对所述SLIC解调中的每一次解调，根据上一次解调的LLR值，分别计算出每组星座点的均值和方差；将计算出的每组星座点的均值进行合并得到该次解调的均值，将计算出的每组星座点的方差进行合并得到该次解调的方差；根据该次解调的均值和方差，计算该次解调的LLR值。

结合第三方面的第三种可能的实现方式，在第四种可能的实现方式中，所述接收机按照以下公式，将计算出的每组星座点的均值进行合并得到所述等效发射信号的每层信号在该次SLIC解调的均值：

$E\left(x_j\right)={\mathrm{\Σ}}_{m=1}^{M}E\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right);$

其中，E(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的均值，E(Ω_m)和Pr(Ω_m)分别表示第m组星座点的均值和第m组星座点的均值对应的先验概率，M表示将所述等效发射信号对应的星座图中的星座点划分为M组，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

结合第三方面的第三种可能的实现方式，在第五种可能的实现方式中，所述接收机按照以下公式，将计算出的每组星座点的方差进行合并得到所述等效发射信号中的每层信号在该次SLIC解调的方差：

$\mathrm{Var}\left(x_j\right)={\mathrm{\Σ}}_{m=1}^M\mathrm{Var}\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)+{\mathrm{\Σ}}_{m=1}^M{E}^2\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)-E\left(x_j\right);$

其中，Var(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的方差，Var(Ω_m)表示第m组星座点的方差，Pr(Ω_m)表示第m组星座点的均值对应的先验概率，M表示将所述等效发射信号对应的星座图中的星座点划分为M组，j＝B+1,……,L，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

结合第三方面、或第三方面的第一种可能的实现方式，在第六种可能的实现方式中，若采用ML解调，所述接收机对等效变换后的接收信号进行解调处理，包括：

根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组；根据所述第二发射信号进行ML解调时所采用的节点选择方法，从每组星座点中，分别选择出设定数量的节点；从已选择出的节点中，分别为所述等效发射信号的每层信号中的每个比特信号选择最终节点；以及根据选择出的最终节点，计算所述等效发射信号的每层信号中的每个比特信号的LLR值。

结合第三方面的第六种可能的实现方式，在第七种可能的实现方式中，所述接收机从已选择出的节点中，为所述等效发射信号的每层信号中的每个比特信号选择最终节点，包括：

从已选择出的节点中，选择出满足设定条件的节点，并将选择出的节点确定为所述等效发射信号的每层信号中的每个比特信号选择最终节点；

其中，设定条件为 $\left\{\begin{array}{c}\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\\ \underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\end{array}\right.,$ {x}表示所述等效发射信号在其对应的星座图中所有可能的取值点，b_j,i表示第j层第i个比特信号，σ²表示噪声功率，表示所述等效发射信号的均值对应的先验概率，y表示接收数据，H表示信道矩阵，x表示所述第二发射信号，j＝B+1,……,L，L表示所述发射信号的总层数，B表示所述第一发射信号的总层数。

结合第三方面的第七种可能的实现方式，在第八种可能的实现方式中，所述接收机根据选择出的最终节点，按照以下公式，计算所述等效发射信号的每层信号中的每个比特信号的LLR值：

${\mathrm{\λ}}_{j,i}=\frac{1}{{\mathrm{\σ}}^2}\{\underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))-\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\};$

其中，λ_j,i表示所述等效发射信号的第j层的第i个比特信号的LLR值。

结合第三方面的第三种可能的实现方式、或第三方面的第六种可能的实现方式，在第九种可能的实现方式中，所述接收机根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组，包括：

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过旋转相同角度得到星座点划分为一组；或者

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过压缩相同压缩量得到星座点划分为一组；或者

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过旋转相同角度并压缩相同压缩量得到星座点划分为一组。

本发明实施例提供的方法、装置和通信设备中，在对接收信号进行解调处理时，将干扰小区对应的预编码矩阵与干扰小区发送的第二发射信号乘积作为等效发射信号，并直接消除该等效发射信号的干扰。由于在第二发射信号里考虑了干扰小区对应的预编码矩阵，这样就不需要在信道估计里考虑干扰小区对应的预编码矩阵了，这样，在进行信道估计时就不用估计干扰小区对应的预编码矩阵的PMI的正确取值，从而提高了接收机抑制干扰的能力，提高了终端的吞吐量性能。

附图说明

图1为本发明实施例提供的一种干扰抑制方法的示意图；

图2A为本发明实施例提供的第二发射信号对应的星座图；

图2B为本发明实施例提供的等效发射信号对应的星座图；

图3为本发明实施例提供的SLIC解调过程的示意图；

图4为本发明实施例提供的ML解调过程的示意图；

图5为本发明实施例提供的一种干扰抑制装置的示意图；

图6为本发明实施例提供的一种通信设备的示意图。

具体实施方式

本发明通过在发射信号中考虑PMI所指示的预编码矩阵，而不在信道估计中考虑PMI所指示的预编码矩阵，由于信道估计时不用估计PMI的正确取值，从而提高了接收机抑制干扰的能力，提高了终端的吞吐量性能。

下面结合说明书附图对本发明实施例作进一步详细描述。应当理解，此处所描述的实施例仅用于说明和解释本发明，并不用于限定本发明。

如图1所示，本发明实施例提供的一种干扰抑制方法，该方法包括：

S11、获取信道参数，所述信道参数包括服务小区对应的第一信道矩阵和干扰小区对应的第二信道矩阵。

本步骤中，根据信道估计，可直接获取到服务小区对应的第一信道矩阵和干扰小区对应的第二信道矩阵第一信道矩阵和第二信道矩阵均为不包含预编码矩阵的信道矩阵。

S12、对接收信号进行等效变换，得到其中，y为接收信号，x₀为所述服务小区发送的第一发射信号，x_I为所述干扰小区发送的第二发射信号，n为噪声信号，H₀为所述第一信道矩阵与所述服务小区对应的预编码矩阵P₀的乘积，为所述第二信道矩阵，P_I为所述干扰小区对应的预编码矩阵。H₀

本步骤中，服务小区对应的预编码矩阵P₀为已知量，而干扰小区对应的预编码矩阵P_I是无法正确估计的量。

S13、将P_Ix_I的乘积作为等效发射信号，对等效变换后的接收信号进行解调处理，以在不获取所述干扰小区对应的预编码矩阵的系数的前提下，消除所述等效发射信号的干扰。

本发明实施例中，在进行解调处理时，将干扰小区对应的预编码矩阵与干扰小区发送的第二发射信号乘积作为等效发射信号，并直接消除该等效发射信号的干扰。由于在第二发射信号里考虑了干扰小区对应的预编码矩阵，这样就不需要在信道估计里考虑干扰小区对应的预编码矩阵了，这样，在进行信道估计时就不用估计干扰小区对应的预编码矩阵的PMI的正确取值，从而提高了接收机抑制干扰的能力，提高了终端的吞吐量性能。

本发明实施例中，各种信号(包括发射信号、第一发射信号、第二发射信号、等效发射信号、噪声信号等)都是以向量的形式表示的。

下面对本发明实施例中所使用的信道模型进行说明。对于背景技术中的公式(2)，可将公式(2)改写为如下公式：

$y=H_0{x}_0+{\tilde{H}}_I{P}_I{x}_I+n=H_0{x}_0+{\tilde{H}}_I{\hat{x}}_I+n---\left(4a\right)$

其中，为信道估计得到的干扰小区对应的第二信道矩阵，即不包含干扰小区对应的预编码矩阵的信道矩阵，为干扰的等效发射信号，具体的：

${\hat{x}}_I=P_I{x}_I[{\hat{x}}_{B+1},...,{\hat{x}}_L]---\left(4b\right)$

显然，等效发射信号的第j层(L≥j≥B+1)的表达式如下：

${\hat{x}}_j={\mathrm{\Σ}}_{i=B+1}^L{\mathrm{\α}}_i{x}_i---\left(4c\right)$

其中，公式(4c)中的α_i为干扰小区对应的预编码矩阵P_I中的向量，可见，等效发射信号的第j层是干扰小区原始发送的发射信号(即第二发射信号)中所有层(即(B+1)～L层)信号的线性组合。

从公式(4a)可以看出，接收机可以把等效符号看成是干扰小区发送的发射信号，例如，依据公式(4a)进行SLIC和ML检测，这样，在进行解调的过程中就没有干扰小区对应的PMI的影响了。

在实施中，S13中对等效变换后的接收信号进行解调处理时，可采用SLIC、ML、或最小均方误差(Minimum Mean Square Error，MMSE)等方式进行解调。具体采用哪种解调方式，可以在标准协议中规定、也可以通过出厂设置规定所使用的解调方式、也可以通过网络侧配置所使用的解调方式，还可以通过额外的信令指示所使用的解调方式，在此不做限定。

当第二发射信号的调制阶数较高(如64QAM)或者第二发射信号的层数较大(如层数为4)时，对应的星座图中的星座点个数就会很多，由此带来的额外复杂度开销也变得不能忽略了。为了进一步降低实现复杂度，较佳地，根据以下步骤确定S13中解调处理过程中所使用的解调方式：

根据等效发射信号对应的星座图中的星座点的个数，确定解调处理中使用的解调方式；

若确定出的星座点的个数大于设定阈值，确定采用第一解调方式对所述发射信号进行解调处理，其中，所述第一解调方式包括MMSE解调；或者，若确定出的星座点的个数小于或等于设定阈值，确定采用第二解调方式对所述发射信号进行解调处理，其中，所述第二解调方式包括SLIC解调和ML解调。

具体的，MMSE解调实际上就是一次迭代的SLIC解调，即通过以下公式计算LLR值：

${\mathrm{\λ}}_{j,i}=\underset{\left\{x_j\right\},b_{j,i}=1}{\max }\left\{\frac{{-|{\tilde{x}}_j-{\widetilde{\mathrm{\φ}}}_j{x}_j|}^2}{\mathrm{Var}\left({\widetilde{\mathrm{\ξ}}}_j\right)}\right\}-\underset{\left\{x_j\right\},b_{j,i}=0}{\max }\left\{\frac{-{|{\tilde{x}}_j-{\widetilde{\mathrm{\φ}}}_j{x}_j|}^2}{\mathrm{Var}\left({\widetilde{\mathrm{\ξ}}}_j\right)}\right\}---\left(5\right)$

其中， ${\hat{x}}_j={\tilde{w}}_j^H(y-\underset{k\≠j}{\mathrm{\Σ}}{h}_{k}E\left(x_k\right));$

${\widetilde{\mathrm{\φ}}}_j={\tilde{w}}_j^H{h}_j,{\tilde{w}}_j^H=h_j^H{({\mathrm{HVH}}^H+{\mathrm{\σ}}^{2}I)}^{-1},V=\mathrm{diag}(\mathrm{Var}\left(x_1\right),...,\mathrm{Var}\left(x_j\right),...,\mathrm{Var}\left(x_L\right));$

$\mathrm{Var}\left({\widetilde{\mathrm{\ξ}}}_j\right)={\widetilde{\mathrm{\φ}}}_j-\mathrm{Var}\left(x_j\right){\widetilde{\mathrm{\φ}}}_j^2;$

上面表达式中的均值和方差分别取值0和1，即E(x_k)＝0，Var(x_j)＝1。

由于在2端口层为2或4端口层为4等无法估计干扰小区对应的PMI的正确值的场景，干扰小区对应的预编码矩阵的自相关矩阵是单位矩阵，则：

${\tilde{w}}_j^H=h_j^H{({\mathrm{HVH}}^H+{\mathrm{\σ}}^{2}I)}^{-1}=h_j^H{(H_0{H_0}^H+{\tilde{H}}_I{{\tilde{H}}_I}^H+{\mathrm{\σ}}^{2}I)}^{-1};$

可见，MMSE解调不需要干扰信号的PMI参数。

在实施中，确定解调处理中使用的解调方式时，所使用的阈值为经验值，可以通过仿真等方式，确定出最优的阈值。

在实施中，根据等效发射信号对应的星座图中的星座点的个数，确定解调处理中使用的解调方式之前，该方法还包括：

根据第二发射信号采用的调制方式，确定出等效发射信号对应的星座图。

具体的，由于本发明实施例中使用等效发射信号替代了干扰小区原始发送的发射信号(即第二发射信号)，等效发射信号中考虑了干扰小区对应的预编码矩阵的影响，因此，等效发射信号对应的星座图相比于第二发射信号对应的星座图有所不同，以QPSK方式为例，对于2端口层为2的场景，第二发射信号对应的星座图如图2A所示，等效发射信号对应的星座图如图2B所示，从图中可以看出，由于干扰小区对应的预编码矩阵的影响，使得干扰小区发送的第二发射信号的星座图发生了改变。虽然不同传输模式使用不同的PMI参数，但只要发射信号使用的端口数和发射信号的层数、以及发射信号使用的调制解调方式确定后，干扰的等效发射信号对应的星座图是不会随着PMI的取值发生变化的。

基于上述任一实施例，作为第一种实现方式，若S13中采用SLIC解调，则对等效变换后的接收信号进行解调处理，包括：

根据如下公式，分别确定每次解调中等效发射信号的每层信号的均值和方差；以及

根据得到的均值和方差确定出每次解调中等效发射信号的每层信号的对数似然比LLR值；

$E\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{x}_j\mathrm{Pr}\left(x_j\right);$

$\mathrm{Var}\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{|x_j-E\left(x_j\right)|}^2\mathrm{Pr}\left(x_j\right);$

其中，E(x_j)表示等效发射信号的第j层信号的均值，Var(x_j)表示等效发射信号的第j层信号的方差，Pr(x_j)表示E(x_j)的先验概率，{x_j}表示等效发射信号的第j层信号在该等效发射信号对应的星座图中所有可能的取值点表示等效发射信号的第j层信号的加权向量，表示等效发射信号的第j层信号的等效噪声，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示第一发射信号的总层数。

基于上述任一实施例，作为第二种实现方式，若S13中采用SLIC解调，则对等效变换后的接收信号进行解调处理，包括：

根据等效发射信号对应的星座图中的星座点的位置，将等效发射信号对应的星座图中的星座点划分为M组，以在计算每组星座点的LLR、均值和方差时能够复用第二发射信号对应的星座图中的星座点的LLR、均值和方差；

针对SLIC解调中的每一次解调，根据上一次解调的LLR值，分别计算出每组星座点的均值和方差；将计算出的每组星座点的均值进行合并得到该次解调的均值，将计算出的每组星座点的方差进行合并得到该次解调的方差；根据该次解调的均值和方差，计算该次解调的LLR值。

具体处理过程如图3所示，对于第i次SLIC解调，先分别计算每组星座点的LLR、均值和方差、其中，每组星座点的LLR、均值和方差都可以复用第二发射信号的星座点的LLR、均值和方差的计算结果)，再将所有组星座点的均值和方差分别合并得到最终的均值和方差，其中，第i+1次SLIC解调在处理时，需要调用第i次SLIC解调处理得到的最终的均值和方差。

在该实现方式下，根据等效发射信号对应的星座图中的星座点的位置，将等效发射信号对应的星座图中的星座点划分为M组之前，该方法还包括：

根据第二发射信号采用的调制方式，确定出等效发射信号对应的星座图。

在该实现方式下，按照以下公式，将计算出的每组星座点的均值进行合并得到等效发射信号的每层信号在该次SLIC解调的均值：

$E\left(x_j\right)={\mathrm{\Σ}}_{m=1}^{M}E\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right);$

其中，E(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的均值，E(Ω_m)和Pr(Ω_m)分别表示第m组星座点的均值和第m组星座点的均值对应的先验概率，M表示将等效发射信号对应的星座图中的星座点划分为M组，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示第一发射信号的总层数。

在该实现方式下，按照以下公式，将计算出的每组星座点的方差进行合并得到等效发射信号中的每层信号在该次SLIC解调的方差：

$\mathrm{Var}\left(x_j\right)={\mathrm{\Σ}}_{m=1}^M\mathrm{Var}\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)+{\mathrm{\Σ}}_{m=1}^M{E}^2\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)-E\left(x_j\right);$

其中，Var(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的方差，Var(Ω_m)表示第m组星座点的方差，Pr(Ω_m)表示第m组星座点的均值对应的先验概率，M表示将等效发射信号对应的星座图中的星座点划分为M组，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示第一发射信号的总层数。

在该实现方式下，根据等效发射信号对应的星座图中的星座点的位置，将等效发射信号对应的星座图中的星座点划分为M组，包括但不限于以下方式中的一种：

等效发射信号对应的星座图中的星座点中，将第二发射信号对应的星座图中的星座点通过旋转相同角度得到星座点划分为一组；

等效发射信号对应的星座图中的星座点中，将第二发射信号对应的星座图中的星座点通过压缩相同压缩量得到星座点划分为一组；

等效发射信号对应的星座图中的星座点中，将第二发射信号对应的星座图中的星座点通过旋转相同角度并压缩相同压缩量得到星座点划分为一组。

需要说明的是，等效发射信号对应的星座图中的星座点中，将位于星座图的原点的星座点划分为一组，其均值和方差均为0。

举例说明，以图2B所示的星座点形状为例，说明第二种实现方式中计算均值和方差的流程。将图2B所示的星座图中的星座点可以分成三组，即“●”表示的星座点(用Ω_b表示该组星座点)、“o”表示的星座点(用Ω_r表示该组星座点)和“△”表示的星座点(用Ω_g表示该组星座点)。

Ω_b就是图2A所示的第二发射信号对应的星座图中的星座点，该组星座点的LLR、均值和方差分别为第二发射信号对应的星座图中的星座点的LLR、均值和方差；Ω_r中的各星座点可以看成是图2A所示的第二发射信号对应的星座图中的各星座点旋转45度再放缩得到的，因此，在复用第二发射信号对应的星座图中的星座点的LLR、均值和方差计算该组星座点的LLR、均值和方差时，先将该组星座点顺时针旋转45度并放大倍，即得到第二发射信号对应的星座图中的星座点，此时，该组星座点的LLR值即为第二发射信号对应的星座图中的星座点的LLR值，再根据该LLR值计算该组星座点的均值和方差，最后将得到的均值反向旋转45度并压缩并将得到的方差压缩为原来的1/2，从而得到Ω_r中的各星座点的均值和方差；Ω_g中的星座点均值和方差都为0。

根据公式(4c)，得到上述三组星座点的先验概率分别为：Pr(Ω_b)＝1/4，Pr(Ω_r)＝1/2，Pr(Ω_g)＝1/4，则：

$E\left(x_j\right)={\mathrm{\Σ}}_{m=1}^{M}E\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)=\frac{1}{4}E\left({\mathrm{\Ω}}_b\right)+\frac{1}{2}E\left({\mathrm{\Ω}}_r\right)+\frac{1}{4}E\left({\mathrm{\Ω}}_g\right),$ 由于E(Ω_g)＝0，因此， $E\left(x_j\right)=\frac{1}{4}E\left({\mathrm{\Ω}}_b\right)+\frac{1}{2}E\left({\mathrm{\Ω}}_r\right);$

$\mathrm{Var}\left(x_j\right)=\frac{1}{4}\mathrm{Var}\left({\mathrm{\Ω}}_b\right)+\frac{1}{2}\mathrm{Var}\left({\mathrm{\Ω}}_r\right)+\frac{1}{4}{E}^2\left({\mathrm{\Ω}}_b\right)+\frac{1}{2}{E}^2\left({\mathrm{\Ω}}_r\right)-E\left(x_j\right).$

该实现方式下，由于在计算每组星座点的均值和方差时都可以复用第二发射信号对应的星座图中星座点的LLR、均值和方差，因此实现复杂度较低。

基于上述任一实施例，作为第三种实现方式，若S13中采用ML解调，则对等效变换后的接收信号进行解调处理，包括：

从等效发射信号对应的星座图中的星座点中，选择出满足设定条件的节点，并将选择出的节点确定为该等效发射信号的每层信号中的每个比特信号的节点；以及

根据选择出的节点，计算等效发射信号的每层信号中的每个比特信号的LLR值；

其中，设定条件为 $\left\{\begin{array}{c}\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\\ \underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\end{array}\right.,$ {x}表示等效发射信号在其对应的星座图中所有可能的取值点，b_j,i表示第j层第i个比特信号，σ²表示噪声功率，表示等效发射信号的均值对应的先验概率，y表示接收信号，H表示信道矩阵，x表示第二发射信号，j＝B+1,……,L，L表示所述发射信号的总层数，B表示第一发射信号的总层数。

该实现方式下，根据选择出的节点，按照以下公式，计算等效发射信号的每层信号中的每个比特信号的LLR值：

${\mathrm{\λ}}_{j,i}=\frac{1}{{\mathrm{\σ}}^2}\{\underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))-\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\};$

其中，λ_j,i表示等效发射信号的第j层的第i个比特信号的LLR值。

基于上述任一实施例，作为第四种实现方式，若S13中采用ML解调，则对等效变换后的接收信号进行解调处理，包括：

根据等效发射信号对应的星座图中的星座点的位置，将等效发射信号对应的星座图中的星座点划分为M组；

根据第二发射信号进行ML解调时所采用的节点选择方法，从每组星座点中，分别选择出设定数量的节点；

从已选择出的节点中，为等效发射信号的每层信号中的每个比特信号选择最终节点；以及

根据选择出的最终节点，计算该等效发射信号的每层信号中的每个比特信号的LLR值。

具体处理过程如图4所示，首先，做预处理，如QR分解；然后，对PMI已知层做节点选择(即前B层x₁,x₂,…,x_B)；再对PMI未知层做节点选择(即)，即针对每组星座点分别进行节点选择(在进行节点选择时，每组星座点均可以复用原ML接收机中QPSK、16QAM或者64QAM的节点选择方法)；然后，进行最终节点选择，最终节点选择需要考虑偏移项

该实现方式下，根据等效发射信号对应的星座图中的星座点的位置，将等效发射信号对应的星座图中的星座点划分为M组之前，该方法还包括：

根据第二发射信号采用的调制方式，确定出等效发射信号对应的星座图。

该实现方式下，根据等效发射信号对应的星座图中的星座点的位置，将等效发射信号对应的星座图中的星座点划分为M组，包括但不限于以下方式中的一种：

等效发射信号对应的星座图中的星座点中，将第二发射信号对应的星座图中的星座点通过旋转相同角度得到星座点划分为一组；

等效发射信号对应的星座图中的星座点中，将第二发射信号对应的星座图中的星座点通过压缩相同压缩量得到星座点划分为一组；

等效发射信号对应的星座图中的星座点中，将第二发射信号对应的星座图中的星座点通过旋转相同角度并压缩相同压缩量得到星座点划分为一组。

需要说明的是，等效发射信号对应的星座图中的星座点中，将位于星座图的原点的星座点划分为一组，其均值和方差均为0。

该实现方式下，从已选择出的节点中，为等效发射信号的每层信号中的每个比特信号选择最终节点，包括：

从已选择出的节点中，选择出满足设定条件的节点，并将选择出的节点确定为等效发射信号的每层信号中的每个比特信号选择最终节点；

其中，设定条件为 $\left\{\begin{array}{c}\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\\ \underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\end{array}\right.,$ {x}表示等效发射信号在其对应的星座图中所有可能的取值点，b_j,i表示第j层第i个比特信号，σ²表示噪声功率，表示等效发射信号的均值对应的先验概率，y表示接收信号，H表示信道矩阵，x表示第二发射信号，，j＝B+1,……,L，L表示接收到的发射信号的总层数，B表示第一发射信号的总层数。

该实现方式下，根据选择出的最终节点，按照以下公式，计算等效发射信号的每层信号中的每个比特信号的LLR值：

${\mathrm{\λ}}_{j,i}=\frac{1}{{\mathrm{\σ}}^2}\{\underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))-\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\};$

其中，λ_j,i表示等效发射信号的第j层的第i个比特信号的LLR值。

举例说明，仍以图2B所示的星座点形状为例，说明第四种实现方式中ML解调的流程。将图2B所示的星座图中的星座点可以分成三组，即“●”表示的星座点(用Ω_b表示该组星座点)、“o”表示的星座点(用Ω_r表示该组星座点)和“△”表示的星座点(用Ω_g表示该组星座点)。

Ω_b就是图2A所示的第二发射信号对应的星座图中的星座点，该组星座点可以复用原始ML解调方式中的节点选择方法(即选择使得||y-Hx||²较小的几个节点，通常可以用查表的方法实现)，假设最终选择出两个节点；Ω_r中的各星座点可以看成是图2A所示的第二发射信号对应的星座图中的各星座点旋转45度再放缩得到的，因此，该组星座点在复用原始ML解调方式中的节点选择方法进行节点选择时，先将该组星座点顺时针旋转45度并放大倍，即得到第二发射信号对应的星座图中的星座点，此时，再复用原始ML解调方式中的节点选择方法进行节点选择，假设最终选择出两个节点；Ω_g中只有一个点，不需要复用原始ML解调方式中的节点选择方法进行节点选择。分组选择之后一共选择出5个节点，最后，从这5个节点中选择出等效发射信号的第j层的每个比特对应的节点，并计算该比特信号的LLR值。

该实现方式下，由于在组星座点在进行节点选择时，均可以复用原始ML解调方式中的节点选择方法，因此实现复杂度较低。

上述方法处理流程可以用软件程序实现，该软件程序可以存储在存储介质中，当存储的软件程序被调用时，执行上述方法步骤。

基于同一发明构思，本发明实施例中还提供了一种干扰抑制装置，由于该装置解决问题的原理与上述干扰抑制方法相似，因此该装置的实施可以参见方法的实施，重复之处不再赘述。

本发明实施例提供的一种干扰抑制装置，如图5所示，该装置包括：

获取模块51，用于获取信道参数，所述信道参数包括服务小区对应的第一信道矩阵和干扰小区对应的第二信道矩阵；

等效变换模块52，对接收信号进行等效变换，得到其中，y为接收信号，x₀为所述服务小区发送的第一发射信号，x_I为所述干扰小区发送的第二发射信号，n为噪声信号，H₀为所述第一信道矩阵与所述服务小区对应的预编码矩阵的乘积，为所述第二信道矩阵，P_I为所述干扰小区对应的预编码矩阵；

解调模块53，用于将P_Ix_I的乘积作为等效发射信号，对等效变换后的接收信号进行解调处理，以在不获取所述干扰小区对应的预编码矩阵的系数的前提下，消除所述等效发射信号的干扰。

本发明实施例提供的装置，可以位于终端设备中，也可以位于网络设备中。

本发明实施例中，各种信号(包括发射信号、第一发射信号、第二发射信号、等效发射信号、噪声信号等)都是以向量的形式表示的。

在实施中，解调模块53对等效变换后的接收信号进行解调处理时，可采用SLIC、ML、或最小均方误差(Minimum Mean Square Error，MMSE)等方式进行解调。具体采用哪种解调方式，可以在标准协议中规定、也可以通过出厂设置规定所使用的解调方式、也可以通过网络侧配置所使用的解调方式，还可以通过额外的信令指示所使用的解调方式，在此不做限定。

当第二发射信号的调制阶数较高(如64QAM)或者第二发射信号的层数较大(如层数为4)时，对应的星座图中的星座点个数就会很多，由此带来的额外复杂度开销也变得不能忽略了。为了进一步降低实现复杂度，较佳地，解调模块53对等效变换后的接收信号进行解调处理之前，还用于：根据所述第二发射信号采用的调制方式，确定出所述等效发射信号对应的星座图；

解调模块53根据以下步骤确定解调处理中使用的解调方式：

根据所述等效发射信号对应的星座图中的星座点的个数，确定解调处理中使用的解调方式；

若确定出的星座点的个数大于设定阈值，确定采用第一解调方式对所述发射信号进行解调处理，其中，所述第一解调方式包括最小均方差MMSE解调；或者，若确定出的星座点的个数小于或等于所述阈值，确定采用第二解调方式对所述发射信号进行解调处理，其中，所述第二解调方式包括符号级干扰消除SLIC解调和最大似然ML解调。

需要说明的是，虽然不同传输模式使用不同的PMI参数，但只要发射信号使用的端口数和发射信号的层数、以及发射信号使用的调制解调方式确定后，干扰的等效发射信号对应的星座图是不会随着PMI的取值发生变化的。

在实施中，确定解调处理中使用的解调方式时，所使用的阈值为经验值，可以通过仿真等方式，确定出最优的阈值。

基于上述任一实施例，作为第一种实现方式，若采用SLIC解调，解调模块53对等效变换后的接收信号进行解调处理，包括：

根据如下公式，分别确定每次解调中所述等效发射信号的每层信号的均值和方差；以及根据得到的均值和方差确定出每次解调中所述等效发射信号的每层信号的对数似然比LLR值；

$E\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{x}_j\mathrm{Pr}\left(x_j\right);$

$\mathrm{Var}\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{|x_j-E\left(x_j\right)|}^2\mathrm{Pr}\left(x_j\right);$

其中，E(x_j)表示所述等效发射信号的第j层信号的均值，Var(x_j)表示所述等效发射信号的第j层信号的方差，Pr(x_j)表示E(x_j)的先验概率，{x_j}表示所述等效发射信号的第j层信号在所述等效发射信号对应的星座图中所有可能的取值点，表示所述等效发射信号的第j层信号的加权向量，表示所述等效发射信号的第j层信号的等效噪声，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

基于上述任一实施例，作为第二种实现方式，若采用SLIC解调，解调模块53对等效变换后的接收信号进行解调处理，包括：

根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组，以在计算每组星座点的LLR、均值和方差时能够复用所述第二发射信号对应的星座图中的星座点的LLR、均值和方差；

针对所述SLIC解调中的每一次解调，根据上一次解调的LLR值，分别计算出每组星座点的均值和方差；将计算出的每组星座点的均值进行合并得到该次解调的均值，将计算出的每组星座点的方差进行合并得到该次解调的方差；根据该次解调的均值和方差，计算该次解调的LLR值。

在该实现方式下，解调模块53根据等效发射信号对应的星座图中的星座点的位置，将等效发射信号对应的星座图中的星座点划分为M组之前，该方法还包括：

根据所述第二发射信号采用的调制方式，确定出所述等效发射信号对应的星座图。

在该实现方式下，解调模块53按照以下公式，将计算出的每组星座点的均值进行合并得到所述等效发射信号的每层信号在该次SLIC解调的均值：

$E\left(x_j\right)={\mathrm{\Σ}}_{m=1}^{M}E\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right);$

其中，E(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的均值，E(Ω_m)和Pr(Ω_m)分别表示第m组星座点的均值和第m组星座点的均值对应的先验概率，M表示将所述等效发射信号对应的星座图中的星座点划分为M组，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

在该实现方式下，解调模块53按照以下公式，将计算出的每组星座点的方差进行合并得到所述等效发射信号中的每层信号在该次SLIC解调的方差：

$\mathrm{Var}\left(x_j\right)={\mathrm{\Σ}}_{m=1}^M\mathrm{Var}\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)+{\mathrm{\Σ}}_{m=1}^M{E}^2\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)-E\left(x_j\right);$

其中，Var(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的方差，Var(Ω_m)表示第m组星座点的方差，Pr(Ω_m)表示第m组星座点的均值对应的先验概率，M表示将所述等效发射信号对应的星座图中的星座点划分为M组，j＝B+1,……,L，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

在该实现方式下，解调模块53根据等效发射信号对应的星座图中的星座点的位置，将等效发射信号对应的星座图中的星座点划分为M组，包括但不限于以下方式中的一种：

等效发射信号对应的星座图中的星座点中，将第二发射信号对应的星座图中的星座点通过旋转相同角度得到星座点划分为一组；

等效发射信号对应的星座图中的星座点中，将第二发射信号对应的星座图中的星座点通过压缩相同压缩量得到星座点划分为一组；

等效发射信号对应的星座图中的星座点中，将第二发射信号对应的星座图中的星座点通过旋转相同角度并压缩相同压缩量得到星座点划分为一组。

需要说明的是，等效发射信号对应的星座图中的星座点中，将位于星座图的原点的星座点划分为一组，其均值和方差均为0。

该实现方式下，由于在计算每组星座点的均值和方差时都可以复用第二发射信号对应的星座图中星座点的LLR、均值和方差，因此实现复杂度较低。

基于上述任一实施例，作为第三种实现方式，若采用ML解调，解调模块53对等效变换后的接收信号进行解调处理，包括：

从等效发射信号对应的星座图中的星座点中，选择出满足设定条件的节点，并将选择出的节点确定为该等效发射信号的每层信号中的每个比特信号的节点；以及

根据选择出的节点，计算等效发射信号的每层信号中的每个比特信号的LLR值；

其中，设定条件为 $\left\{\begin{array}{c}\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\\ \underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\end{array}\right.,$ {x}表示等效发射信号在其对应的星座图中所有可能的取值点，b_j,i表示第j层第i个比特信号，σ²表示噪声功率，表示等效发射信号的均值对应的先验概率，y表示接收信号，H表示信道矩阵，x表示第二发射信号，j＝B+1,……,L，L表示所述发射信号的总层数，B表示第一发射信号的总层数。

该实现方式下，解调模块53根据选择出的节点，按照以下公式，计算等效发射信号的每层信号中的每个比特信号的LLR值：

${\mathrm{\λ}}_{j,i}=\frac{1}{{\mathrm{\σ}}^2}\{\underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))-\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\};$

其中，λ_j,i表示等效发射信号的第j层的第i个比特信号的LLR值。

基于上述任一实施例，作为第四种实现方式，若采用ML解调，解调模块53对等效变换后的接收信号进行解调处理，包括：

根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组；根据所述第二发射信号进行ML解调时所采用的节点选择方法，从每组星座点中，分别选择出设定数量的节点；从已选择出的节点中，分别为所述等效发射信号的每层信号中的每个比特信号选择最终节点；以及根据选择出的最终节点，计算所述等效发射信号的每层信号中的每个比特信号的LLR值。

该实现方式下，解调模块53从已选择出的节点中，为所述等效发射信号的每层信号中的每个比特信号选择最终节点，包括：

从已选择出的节点中，选择出满足设定条件的节点，并将选择出的节点确定为所述等效发射信号的每层信号中的每个比特信号选择最终节点；

其中，设定条件为 $\left\{\begin{array}{c}\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\\ \underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\end{array}\right.,$ {x}表示所述等效发射信号在其对应的星座图中所有可能的取值点，b_j,i表示第j层第i个比特信号，σ²表示噪声功率，表示所述等效发射信号的均值对应的先验概率，y表示接收数据，H表示信道矩阵，x表示所述第二发射信号，j＝B+1,……,L，L表示所述发射信号的总层数，B表示所述第一发射信号的总层数。

该实现方式下，解调模块53根据选择出的最终节点，按照以下公式，计算所述等效发射信号的每层信号中的每个比特信号的LLR值：

${\mathrm{\λ}}_{j,i}=\frac{1}{{\mathrm{\σ}}^2}\{\underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))-\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\};$

其中，λ_j,i表示所述等效发射信号的第j层的第i个比特信号的LLR值。

该实现方式下，解调模块53根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组，包括：

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过旋转相同角度得到星座点划分为一组；或者

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过压缩相同压缩量得到星座点划分为一组；或者

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过旋转相同角度并压缩相同压缩量得到星座点划分为一组。

基于同一发明构思，本发明实施例提供了一种通信设备，如图6所示，该通信设备包括：

处理器62，用于获取模块，用于获取信道参数，所述信道参数包括服务小区对应的第一信道矩阵和干扰小区对应的第二信道矩阵；

接收机61，用于等效变换模块，用于对接收信号进行等效变换，得到其中，y为接收信号，x₀为所述服务小区发送的第一发射信号，x_I为所述干扰小区发送的第二发射信号，n为噪声信号，H₀为所述第一信道矩阵与所述服务小区对应的预编码矩阵的乘积，为所述第二信道矩阵，P_I为所述干扰小区对应的预编码矩阵；将P_Ix_I的乘积作为等效发射信号，对等效变换后的接收信号进行解调处理，以在不获取所述干扰小区对应的预编码矩阵的系数的前提下，消除所述等效发射信号的干扰；以及将解调处理后的信号传输给处理器62进行处理。

接收机61与处理器62之间通过总线连接。

本发明实施例提供的通信设备，可以是终端设备，也可以是网络设备。

本发明实施例中，各种信号(包括发射信号、第一发射信号、第二发射信号、等效发射信号、噪声信号等)都是以向量的形式表示的。

在实施中，接收机61中对等效变换后的接收信号进行解调处理时，可采用SLIC、ML、或MMSE等方式进行解调。具体采用哪种解调方式，可以在标准协议中规定、也可以通过出厂设置规定所使用的解调方式、也可以通过网络侧配置所使用的解调方式，还可以通过额外的信令指示所使用的解调方式，在此不做限定。

当第二发射信号的调制阶数较高(如64QAM)或者第二发射信号的层数较大(如层数为4)时，对应的星座图中的星座点个数就会很多，由此带来的额外复杂度开销也变得不能忽略了。为了进一步降低实现复杂度，较佳地，接收机61对等效变换后的接收信号进行解调处理之前，还用于：根据所述第二发射信号采用的调制方式，确定出所述等效发射信号对应的星座图；

接收机61根据以下步骤确定解调处理中使用的解调方式：

根据所述等效发射信号对应的星座图中的星座点的个数，确定解调处理中使用的解调方式；

若确定出的星座点的个数大于设定阈值，确定采用第一解调方式对所述发射信号进行解调处理，其中，所述第一解调方式包括最小均方差MMSE解调；或者，若确定出的星座点的个数小于或等于所述阈值，确定采用第二解调方式对所述发射信号进行解调处理，其中，所述第二解调方式包括符号级干扰消除SLIC解调和最大似然ML解调。

需要说明的是，虽然不同传输模式使用不同的PMI参数，但只要发射信号使用的端口数和发射信号的层数、以及发射信号使用的调制解调方式确定后，干扰的等效发射信号对应的星座图是不会随着PMI的取值发生变化的。

在实施中，确定解调处理中使用的解调方式时，所使用的阈值为经验值，可以通过仿真等方式，确定出最优的阈值。

基于上述任一实施例，作为第一种实现方式，若采用SLIC解调，接收机61对等效变换后的接收信号进行解调处理，包括：

根据如下公式，分别确定每次解调中所述等效发射信号的每层信号的均值和方差；以及根据得到的均值和方差确定出每次解调中所述等效发射信号的每层信号的对数似然比LLR值；

$E\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{x}_j\mathrm{Pr}\left(x_j\right);$

$\mathrm{Var}\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{|x_j-E\left(x_j\right)|}^2\mathrm{Pr}\left(x_j\right);$

其中，E(x_j)表示所述等效发射信号的第j层信号的均值，Var(x_j)表示所述等效发射信号的第j层信号的方差，Pr(x_j)表示E(x_j)的先验概率，{x_j}表示所述等效发射信号的第j层信号在所述等效发射信号对应的星座图中所有可能的取值点，表示所述等效发射信号的第j层信号的加权向量，表示所述等效发射信号的第j层信号的等效噪声，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

基于上述任一实施例，作为第二种实现方式，若采用SLIC解调，接收机61对等效变换后的接收信号进行解调处理，包括：

根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组，以在计算每组星座点的LLR、均值和方差时能够复用所述第二发射信号对应的星座图中的星座点的LLR、均值和方差；

针对所述SLIC解调中的每一次解调，根据上一次解调的LLR值，分别计算出每组星座点的均值和方差；将计算出的每组星座点的均值进行合并得到该次解调的均值，将计算出的每组星座点的方差进行合并得到该次解调的方差；根据该次解调的均值和方差，计算该次解调的LLR值。

在该实现方式下，接收机61根据等效发射信号对应的星座图中的星座点的位置，将等效发射信号对应的星座图中的星座点划分为M组之前，该方法还包括：

根据所述第二发射信号采用的调制方式，确定出所述等效发射信号对应的星座图。

在该实现方式下，接收机61按照以下公式，将计算出的每组星座点的均值进行合并得到所述等效发射信号的每层信号在该次SLIC解调的均值：

$E\left(x_j\right)={\mathrm{\Σ}}_{m=1}^{M}E\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right);$

其中，E(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的均值，E(Ω_m)和Pr(Ω_m)分别表示第m组星座点的均值和第m组星座点的均值对应的先验概率，M表示将所述等效发射信号对应的星座图中的星座点划分为M组，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

在该实现方式下，接收机61按照以下公式，将计算出的每组星座点的方差进行合并得到所述等效发射信号中的每层信号在该次SLIC解调的方差：

$\mathrm{Var}\left(x_j\right)={\mathrm{\Σ}}_{m=1}^M\mathrm{Var}\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)+{\mathrm{\Σ}}_{m=1}^M{E}^2\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)-E\left(x_j\right);$

其中，Var(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的方差，Var(Ω_m)表示第m组星座点的方差，Pr(Ω_m)表示第m组星座点的均值对应的先验概率，M表示将所述等效发射信号对应的星座图中的星座点划分为M组，j＝B+1,……,L，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

在该实现方式下，接收机61根据等效发射信号对应的星座图中的星座点的位置，将等效发射信号对应的星座图中的星座点划分为M组，包括但不限于以下方式中的一种：

等效发射信号对应的星座图中的星座点中，将第二发射信号对应的星座图中的星座点通过旋转相同角度得到星座点划分为一组；

等效发射信号对应的星座图中的星座点中，将第二发射信号对应的星座图中的星座点通过压缩相同压缩量得到星座点划分为一组；

等效发射信号对应的星座图中的星座点中，将第二发射信号对应的星座图中的星座点通过旋转相同角度并压缩相同压缩量得到星座点划分为一组。

需要说明的是，等效发射信号对应的星座图中的星座点中，将位于星座图的原点的星座点划分为一组，其均值和方差均为0。

该实现方式下，由于在计算每组星座点的均值和方差时都可以复用第二发射信号对应的星座图中星座点的LLR、均值和方差，因此实现复杂度较低。

基于上述任一实施例，作为第三种实现方式，若采用ML解调，接收机61对等效变换后的接收信号进行解调处理，包括：

从等效发射信号对应的星座图中的星座点中，选择出满足设定条件的节点，并将选择出的节点确定为该等效发射信号的每层信号中的每个比特信号的节点；以及

根据选择出的节点，计算等效发射信号的每层信号中的每个比特信号的LLR值；

其中，设定条件为 $\left\{\begin{array}{c}\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\\ \underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\end{array}\right.,$ {x}表示等效发射信号在其对应的星座图中所有可能的取值点，b_j,i表示第j层第i个比特信号，σ²表示噪声功率，表示等效发射信号的均值对应的先验概率，y表示接收信号，H表示信道矩阵，x表示第二发射信号，j＝B+1,……,L，L表示所述发射信号的总层数，B表示第一发射信号的总层数。

该实现方式下，接收机61根据选择出的节点，按照以下公式，计算等效发射信号的每层信号中的每个比特信号的LLR值：

${\mathrm{\λ}}_{j,i}=\frac{1}{{\mathrm{\σ}}^2}\{\underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))-\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\};$

其中，λ_j,i表示等效发射信号的第j层的第i个比特信号的LLR值。

基于上述任一实施例，作为第四种实现方式，若采用ML解调，接收机61对等效变换后的接收信号进行解调处理，包括：

根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组；根据所述第二发射信号进行ML解调时所采用的节点选择方法，从每组星座点中，分别选择出设定数量的节点；从已选择出的节点中，分别为所述等效发射信号的每层信号中的每个比特信号选择最终节点；以及根据选择出的最终节点，计算所述等效发射信号的每层信号中的每个比特信号的LLR值。

该实现方式下，接收机61从已选择出的节点中，为所述等效发射信号的每层信号中的每个比特信号选择最终节点，包括：

从已选择出的节点中，选择出满足设定条件的节点，并将选择出的节点确定为所述等效发射信号的每层信号中的每个比特信号选择最终节点；

其中，设定条件为 $\left\{\begin{array}{c}\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\\ \underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\end{array}\right.,$ {x}表示所述等效发射信号在其对应的星座图中所有可能的取值点，b_j,i表示第j层第i个比特信号，σ²表示噪声功率，表示所述等效发射信号的均值对应的先验概率，y表示接收数据，H表示信道矩阵，x表示所述第二发射信号，j＝B+1,……,L，L表示所述发射信号的总层数，B表示所述第一发射信号的总层数。

该实现方式下，接收机61根据选择出的最终节点，按照以下公式，计算所述等效发射信号的每层信号中的每个比特信号的LLR值：

${\mathrm{\λ}}_{j,i}=\frac{1}{{\mathrm{\σ}}^2}\{\underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))-\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\};$

其中，λ_j,i表示所述等效发射信号的第j层的第i个比特信号的LLR值。

该实现方式下，接收机61根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组，包括：

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过旋转相同角度得到星座点划分为一组；或者

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过压缩相同压缩量得到星座点划分为一组；或者

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过旋转相同角度并压缩相同压缩量得到星座点划分为一组。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。

显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (20)

1.一种干扰抑制方法，其特征在于，该方法包括：

获取信道参数，所述信道参数包括服务小区对应的第一信道矩阵和干扰小区对应的第二信道矩阵；

对接收信号进行等效变换，得到其中，y为接收信号，x₀为所述服务小区发送的第一发射信号，x_I为所述干扰小区发送的第二发射信号，n为噪声信号，H₀为所述第一信道矩阵与所述服务小区对应的预编码矩阵的乘积，为所述第二信道矩阵，P_I为所述干扰小区对应的预编码矩阵；

将P_Ix_I的乘积作为等效发射信号，对等效变换后的接收信号进行解调处理，以在不获取所述干扰小区对应的预编码矩阵的系数的前提下，消除所述等效发射信号的干扰。

2.如权利要求1所述的方法，其特征在于，对等效变换后的接收信号进行解调处理之前，还包括：根据所述第二发射信号采用的调制方式，确定出所述等效发射信号对应的星座图；

根据以下步骤确定解调处理中使用的解调方式：

根据所述等效发射信号对应的星座图中的星座点的个数，确定解调处理中使用的解调方式；

若确定出的星座点的个数大于设定阈值，确定采用第一解调方式对所述发射信号进行解调处理，其中，所述第一解调方式包括最小均方差MMSE解调；或者，若确定出的星座点的个数小于或等于所述阈值，确定采用第二解调方式对所述发射信号进行解调处理，其中，所述第二解调方式包括符号级干扰消除SLIC解调和最大似然ML解调。

3.如权利要求1或2所述的方法，其特征在于，若采用SLIC解调，对等效变换后的接收信号进行解调处理，包括：

根据如下公式，分别确定每次解调中所述等效发射信号的每层信号的均值和方差；

根据得到的均值和方差确定出每次解调中所述等效发射信号的每层信号的对数似然比LLR值；

$E\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{x}_j\mathrm{Pr}\left(x_j\right);$

$\mathrm{Var}\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{|x_j-E\left(x_j\right)|}^2\mathrm{Pr}\left(x_j\right);$

其中，E(x_j)表示所述等效发射信号的第j层信号的均值，Var(x_j)表示所述等效发射信号的第j层信号的方差，Pr(x_j)表示E(x_j)的先验概率，{x_j}表示所述等效发射信号的第j层信号在所述等效发射信号对应的星座图中所有可能的取值点，表示所述等效发射信号的第j层信号的加权向量，表示所述等效发射信号的第j层信号的等效噪声，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

4.如权利要求1或2所述的方法，其特征在于，若采用SLIC解调，对等效变换后的接收信号进行解调处理，包括：

根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组，以在计算每组星座点的LLR、均值和方差时能够复用所述第二发射信号对应的星座图中的星座点的LLR、均值和方差；

针对所述SLIC解调中的每一次解调，根据上一次解调的LLR值，分别计算出每组星座点的均值和方差；将计算出的每组星座点的均值进行合并得到该次解调的均值，将计算出的每组星座点的方差进行合并得到该次解调的方差；根据该次解调的均值和方差，计算该次解调的LLR值。

5.如权利要求4所述的方法，其特征在于，按照以下公式，将计算出的每组星座点的均值进行合并得到所述等效发射信号的每层信号在该次SLIC解调的均值：

$E\left(x_j\right)={\mathrm{\Σ}}_{m=1}^{M}E\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right);$

其中，E(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的均值，E(Ω_m)和Pr(Ω_m)分别表示第m组星座点的均值和第m组星座点的均值对应的先验概率，M表示将所述等效发射信号对应的星座图中的星座点划分为M组，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

6.如权利要求4所述的方法，其特征在于，按照以下公式，将计算出的每组星座点的方差进行合并得到所述等效发射信号中的每层信号在该次SLIC解调的方差：

$\mathrm{Var}\left(x_j\right)={\mathrm{\Σ}}_{m=1}^M\mathrm{Var}\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)+{\mathrm{\Σ}}_{m=1}^M{E}^2\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)-E\left(x_j\right);$

其中，Var(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的方差，Var(Ω_m)表示第m组星座点的方差，Pr(Ω_m)表示第m组星座点的均值对应的先验概率，M表示将所述等效发射信号对应的星座图中的星座点划分为M组，j＝B+1,……,L，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

7.如权利要求1或2所述的方法，其特征在于，若采用ML解调，对等效变换后的接收信号进行解调处理，包括：

根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组；

根据所述第二发射信号进行ML解调时所采用的节点选择方法，从每组星座点中，分别选择出设定数量的节点；

从已选择出的节点中，分别为所述等效发射信号的每层信号中的每个比特信号选择最终节点；以及

根据选择出的最终节点，计算所述等效发射信号的每层信号中的每个比特信号的LLR值。

8.如权利要求7所述的方法，其特征在于，从已选择出的节点中，为所述等效发射信号的每层信号中的每个比特信号选择最终节点，包括：

从已选择出的节点中，选择出满足设定条件的节点，并将选择出的节点确定为所述等效发射信号的每层信号中的每个比特信号选择最终节点；

其中，设定条件为 $\left\{\begin{array}{c}\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\\ \underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\end{array}\right.,$ {x}表示所述等效发射信号在其对应的星座图中所有可能的取值点，b_j,i表示第j层第i个比特信号，σ²表示噪声功率，表示所述等效发射信号的均值对应的先验概率，y表示接收数据，H表示信道矩阵，x表示所述第二发射信号，j＝B+1,……,L，L表示所述发射信号的总层数，B表示所述第一发射信号的总层数。

9.如权利要求8所述的方法，其特征在于，根据选择出的最终节点，按照以下公式，计算所述等效发射信号的每层信号中的每个比特信号的LLR值：

${\mathrm{\λ}}_{j,i}=\frac{1}{{\mathrm{\σ}}^2}\{\underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))-\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\};$

其中，λ_j,i表示所述等效发射信号的第j层的第i个比特信号的LLR值。

10.如权利要求4或7所述的方法，其特征在于，根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组，包括：

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过旋转相同角度得到星座点划分为一组；或者

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过压缩相同压缩量得到星座点划分为一组；或者

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过旋转相同角度并压缩相同压缩量得到星座点划分为一组。

11.一种干扰抑制装置，其特征在于，该装置包括：

获取模块，用于获取信道参数，所述信道参数包括服务小区对应的第一信道矩阵和干扰小区对应的第二信道矩阵；

等效变换模块，用于对接收信号进行等效变换，得到其中，y为接收信号，x₀为所述服务小区发送的第一发射信号，x_I为所述干扰小区发送的第二发射信号，n为噪声信号，H₀为所述第一信道矩阵与所述服务小区对应的预编码矩阵的乘积，为所述第二信道矩阵，P_I为所述干扰小区对应的预编码矩阵；

解调模块，用于将P_Ix_I的乘积作为等效发射信号，对等效变换后的接收信号进行解调处理，以在不获取所述干扰小区对应的预编码矩阵的系数的前提下，消除所述等效发射信号的干扰。

12.如权利要求11所述的装置，其特征在于，所述解调模块对等效变换后的接收信号进行解调处理之前，还用于：根据所述第二发射信号采用的调制方式，确定出所述等效发射信号对应的星座图；

所述解调模块根据以下步骤确定解调处理中使用的解调方式：

根据所述等效发射信号对应的星座图中的星座点的个数，确定解调处理中使用的解调方式；

若确定出的星座点的个数大于设定阈值，确定采用第一解调方式对所述发射信号进行解调处理，其中，所述第一解调方式包括最小均方差MMSE解调；或者，若确定出的星座点的个数小于或等于所述阈值，确定采用第二解调方式对所述发射信号进行解调处理，其中，所述第二解调方式包括符号级干扰消除SLIC解调和最大似然ML解调。

13.如权利要求11或12所述的装置，其特征在于，若采用SLIC解调，所述解调模块对等效变换后的接收信号进行解调处理，包括：

根据如下公式，分别确定每次解调中所述等效发射信号的每层信号的均值和方差；以及根据得到的均值和方差确定出每次解调中所述等效发射信号的每层信号的对数似然比LLR值；

$E\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{x}_j\mathrm{Pr}\left(x_j\right);$

$\mathrm{Var}\left(x_j\right)=\underset{\left\{x_j\right\}}{\mathrm{\Σ}}{|x_j-E\left(x_j\right)|}^2\mathrm{Pr}\left(x_j\right);$

其中，E(x_j)表示所述等效发射信号的第j层信号的均值，Var(x_j)表示所述等效发射信号的第j层信号的方差，Pr(x_j)表示E(x_j)的先验概率，{x_j}表示所述等效发射信号的第j层信号在所述等效发射信号对应的星座图中所有可能的取值点，表示所述等效发射信号的第j层信号的加权向量，表示所述等效发射信号的第j层信号的等效噪声，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

14.如权利要求11或12所述的装置，其特征在于，若采用SLIC解调，所述解调模块对等效变换后的接收信号进行解调处理，包括：

根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组，以在计算每组星座点的LLR、均值和方差时能够复用所述第二发射信号对应的星座图中的星座点的LLR、均值和方差；

针对所述SLIC解调中的每一次解调，根据上一次解调的LLR值，分别计算出每组星座点的均值和方差；将计算出的每组星座点的均值进行合并得到该次解调的均值，将计算出的每组星座点的方差进行合并得到该次解调的方差；根据该次解调的均值和方差，计算该次解调的LLR值。

15.如权利要求14所述的装置，其特征在于，所述解调模块按照以下公式，将计算出的每组星座点的均值进行合并得到所述等效发射信号的每层信号在该次SLIC解调的均值：

$E\left(x_j\right)={\mathrm{\Σ}}_{m=1}^{M}E\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right);$

其中，E(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的均值，E(Ω_m)和Pr(Ω_m)分别表示第m组星座点的均值和第m组星座点的均值对应的先验概率，M表示将所述等效发射信号对应的星座图中的星座点划分为M组，j∈{B+1,……,L}，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

16.如权利要求14所述的装置，其特征在于，所述解调模块按照以下公式，将计算出的每组星座点的方差进行合并得到所述等效发射信号中的每层信号在该次SLIC解调的方差：

$\mathrm{Var}\left(x_j\right)={\mathrm{\Σ}}_{m=1}^M\mathrm{Var}\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)+{\mathrm{\Σ}}_{m=1}^M{E}^2\left({\mathrm{\Ω}}_m\right)\mathrm{Pr}\left({\mathrm{\Ω}}_m\right)-E\left(x_j\right);$

其中，Var(x_j)表示该次SLIC解调中所述等效发射信号的第j层信号的方差，Var(Ω_m)表示第m组星座点的方差，Pr(Ω_m)表示第m组星座点的均值对应的先验概率，M表示将所述等效发射信号对应的星座图中的星座点划分为M组，j＝B+1,……,L，L表示接收到的发射信号的总层数，B表示所述第一发射信号的总层数。

17.如权利要求11或12所述的装置，其特征在于，若采用ML解调，所述解调模块对等效变换后的接收信号进行解调处理，包括：

根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组；根据所述第二发射信号进行ML解调时所采用的节点选择方法，从每组星座点中，分别选择出设定数量的节点；从已选择出的节点中，分别为所述等效发射信号的每层信号中的每个比特信号选择最终节点；以及根据选择出的最终节点，计算所述等效发射信号的每层信号中的每个比特信号的LLR值。

18.如权利要求17所述的装置，其特征在于，所述解调模块从已选择出的节点中，为所述等效发射信号的每层信号中的每个比特信号选择最终节点，包括：

从已选择出的节点中，选择出满足设定条件的节点，并将选择出的节点确定为所述等效发射信号的每层信号中的每个比特信号选择最终节点；

其中，设定条件为 $\left\{\begin{array}{c}\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\\ \underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\end{array}\right.,$ {x}表示所述等效发射信号在其对应的星座图中所有可能的取值点，b_j,i表示第j层第i个比特信号，σ²表示噪声功率，表示所述等效发射信号的均值对应的先验概率，y表示接收数据，H表示信道矩阵，x表示所述第二发射信号，j＝B+1,……,L，L表示所述发射信号的总层数，B表示所述第一发射信号的总层数。

19.如权利要求18所述的装置，其特征在于，所述解调模块根据选择出的最终节点，按照以下公式，计算所述等效发射信号的每层信号中的每个比特信号的LLR值：

${\mathrm{\λ}}_{j,i}=\frac{1}{{\mathrm{\σ}}^2}\{\underset{\left\{x\right\}:b_{j,i}=1}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))-\underset{\left\{x\right\}:b_{j,i}=0}{\min }({\left|\right|y-\mathrm{Hx}\left|\right|}^2-{\mathrm{\σ}}^2\ln \mathrm{Pr}\left({\hat{x}}_I\right))\};$

其中，λ_j,i表示所述等效发射信号的第j层的第i个比特信号的LLR值。

20.如权利要求14或17所述的装置，其特征在于，所述解调模块根据所述等效发射信号对应的星座图中的星座点的位置，将所述等效发射信号对应的星座图中的星座点划分为M组，包括：

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过旋转相同角度得到星座点划分为一组；或者

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过压缩相同压缩量得到星座点划分为一组；或者

所述等效发射信号对应的星座图中的星座点中，将所述第二发射信号对应的星座图中的星座点通过旋转相同角度并压缩相同压缩量得到星座点划分为一组。