CN101753176A - 一种干扰抑制合并方法及系统 - Google Patents

Abstract

本发明涉及一种干扰抑制合并方法，包括：分别计算在基站侧接收端ant数量的天线接收的物理资源块中导频符号的子载波上的干扰加噪声信号；根据所述干扰加噪声信号，计算得到所述物理资源块中的干扰加噪声协方差矩阵；计算所述物理资源块中子载波上的信道矩阵；根据所述干扰加噪声协方差矩阵和信道矩阵，计算所述物理资源块中子载波上权值合并矩阵；根据所述信道矩阵和所述上权值合并矩阵计算所述物理资源块中子载波上合并值调整系数；根据计算出的所述权值合并矩阵，计算物理资源块中非导频符号的子载波上的接收信号合并值；根据计算出的所述合并值调整系数调整物理资源块中非导频符号的接收信号合并值，得到调整结果。

Description

一种干扰抑制合并方法及系统

技术领域

本发明涉及通信领域，特别是涉及一种干扰抑制合并方法及系统。

背景技术

ICI(Inter cell interference，小区之间干扰)是蜂窝移动通信系统的一个固有问题，传统的解决办法是采用频率复用，复用系数只有特定的几个选择，如1、3、7等。复用系数＝1即相邻小区都使用相同的频率资源，这时在小区边缘的干扰很严重。较低的复用系数(3或7)可以有效地抑制ICI，但频谱效率将减低。

ICI抑制的技术包括：干扰随机化(ICI randomization)，干扰协调(ICIcoordination)和干扰消除(ICI cancellation)。

干扰随机化是指通过随机化干扰信号来进行干扰压制的一类技术。包括小区特定的加扰和小区特定的交织，也叫IDMA(Interlerved Division MutipleAccess，空分多址)。干扰随机化计算复杂度较小，但只是白化了干扰，并没有真正减少系统的干扰信号，因此带来的信噪比改善程度有限。研究结果表明，单独应用干扰随机化技术不能满足未来移动通信系统的信噪比要求。

干扰消除技术是利用对于干扰信号的估计，在进行信号检测时尽可能地消除干扰小区的信号，从而在接收端得到处理增益。干扰消除技术复杂度高、应用条件严格，而且只能消除一些强干扰源，而实际系统中干扰主要是由很多小的干扰叠加在一起产生的，因此这项技术在实际应用中研究效果有限。

干扰协调技术则是通过在小区间协调分配资源，尽量使相邻小区所用资源正交，从而使小区间的干扰减少。方法包括软频率复用和部分功率控制，前者是对资源调度的限制(即避免干扰小区使用可能造成干扰的资源块)，后者是对某个资源块的限制(如控制干扰小区在可能造成干扰的资源块内的发射功率)。

IRC(Interference Rejection Combining，干扰抑制合并)算法，这是一种估计出干扰(认为是有色噪声)和噪声的相关矩阵，从而对干扰起到一种抑制作用的分集合并技术。

IRC系统技术的背景描述：

将无线窄带无线系统建模，在模型中信道的影响可以用一个复数乘以传输信号，这个模型是对OFDM和OFDMA中每一个子载波适用的，设x是任一子载波上的传输信号，认为系统中有一个发射天线还有L个接收天线，接收的信号可以表示为：

r＝hx+n    (1)

h＝[h₀，h₁，...h_L-1]^T是用户信道接收向量，其中是接收天线1的复数信道增益，n＝[n₀，n₁，...n_L-1]^T是噪声加干扰向量，其中每一个n_l的方差是ρ_l ²。信号加噪声部分由下式给出：

$n_l={\mathrm{\υ}}_l+\underset{m=1}{\overset{M}{\mathrm{\Σ}}}{P}_m{h}_{l,m}---\left(2\right)$

其中υ_l元素是热噪声，具有方差σ_υ ²，对所有的天线都是一样的。M是干扰信号总数，P_m是第m个干扰的平均功率，h_l，m是第m个干扰源到1个接收天线间的信道增益，通用的合并消除干扰可以用下式表示：

y＝w^Tr＝w^Thx+w^Tn    (3)

信号比噪声加干扰比由下式给出：

$\mathrm{\γ}=\frac{E\left[w^T{\mathrm{hh}}^H{w}^{*}{\left|x\right|}^2\right]}{E\left[w^T{\mathrm{nn}}^H{w}^{*}\right]}=\frac{w^T{\mathrm{hh}}^H{w}^{*}{P}_0}{w^T{\mathrm{Rw}}^{*}}---\left(4\right)$

其中P₀是用户设计平均功率，R＝E[nn^H]是噪声加干扰协方差矩阵，例如我们有2个接收天线，协方差矩阵由下式给出：

$R=\left[\begin{array}{cc}{\mathrm{\ρ}}_1^2& {\mathrm{\ρ}}_{\mathrm{1,2}}\\ {\mathrm{\ρ}}_{\mathrm{1,2}}^{*}& {\mathrm{\ρ}}_2^2\end{array}\right]---\left(5\right)$

其中，

${\mathrm{\ρ}}_l^2={\mathrm{\σ}}_{\mathrm{\υ}}^2+\underset{m=1}{\overset{M}{\mathrm{\Σ}}}{P}_m{\left|h_{l,m}\right|}^2---\left(6\right)$

${\mathrm{\ρ}}_{\mathrm{1,2}}=\underset{m=1}{\overset{M}{\mathrm{\Σ}}}{P}_m{h}_{1,m}{h}_{2,m}^{*}---\left(7\right)$

如果有1个占主导地位的干扰的情况下，非对角线元素就更加有意义。我们假设系数h_m仅代表快衰落的情况，功率P_m代表已经结合了路径损耗和阴影衰落的效果。假设第1个干扰(m＝1)是最强的干扰源，我们可以定义主要干扰比率(DIR)为：

$\mathrm{DIR}=\frac{P_1}{{\mathrm{\σ}}_{\mathrm{\υ}}^2+{\mathrm{\Σ}}_{m=2}^M{P}_m}---\left(8\right)$

如果假设其他非主导地位干扰源可以看成是白噪声的一部分，可将式(6)和(7)重新写为：

${\mathrm{\ρ}}_l^2\≈{\mathrm{\σ}}_{\mathrm{\υ}}^2+P_1{\left|h_{l,1}\right|}^2---\left(9\right)$

${\mathrm{\ρ}}_{\mathrm{1,2}}\≈P_1{h}_{\mathrm{1,1}}{h}_{\mathrm{2,1}}^{*}---\left(10\right)$

由最大化γ(信号比干扰加噪声比率)准则最后可得到合并权值矩阵：

$w_{\mathrm{IRC}}^T=h^H{R}^{-1}---\left(11\right)$

对于IRC算法而言，最关键的技术是求解噪声加干扰的协方差矩阵(式(5))，根据IRC原理其协方差矩阵求解需要得到干扰功率和干扰的信道信息。但实际上，IRC系统中这些信息是无法求得的。目前近似的方法可以直接采用接收信号的自相关矩阵近似估计干扰与噪声的协方差矩阵，但是估计出的噪声加干扰协方差矩阵效果不佳。

发明内容

本发明的目的是提供一种干扰抑制合并方法及系统，解出噪声加干扰矩阵协方差矩阵，进而实现IRC方法对干扰进行抑制，即得到一套在上行信道中可行的IRC处理方法。

本发明提供一种干扰抑制合并方法，包括：

A、分别计算在基站侧接收端的ant个数天线接收的物理资源块中导频符号子载波上的干扰加噪声信号n_ant，SCn(N)；

B、根据所述干扰加噪声信号n_ant，SCn(N)，计算得到所述物理资源块的干扰加噪声协方差矩阵R(N)；

C、计算所述物理资源块中子载波上的信道矩阵H_SCn(N)；

D、根据所述干扰加噪声协方差矩阵R(N)和信道矩阵H_SCn(N)，计算所述物理资源块中子载波上权值合并矩阵w_SCn(N)；

E、根据所述信道矩阵H_SCn(N)和所述上权值合并矩阵w_SCn(N)计算所述物理资源块中子载波上合并值调整系数α_SCn(N)；

F、根据计算出的所述权值合并矩阵w_SCn(N)，计算物理资源块中非导频符号的子载波上的接收信号合并值c_SCn(N)；

G、根据计算出的所述合并值调整系数α_SCn(N)调整物理资源块中非导频符号的子载波上的接收信号合并值c_SCn(N)，得到调整结果c′_SCn(N)。

优选地，所述步骤A具体为：

计算一个物理资源块中12个子载波上的干扰加噪声信号n_ant，SCn(N)；

所述步骤A通过公式：

n_ant，SCn(N)＝r_ant，SCn(N)-h_ant，SCn(N)x_SCn(N)计算；

其中，r_ant，SCn(N)为载波上的接收信号；

h_ant，SCn(N)为信道系数；

x_SCn(N)为本地产生的导频序列；

所述ant取值0或1；

SCn取值1至12。

优选地，所述步骤B通过公式：

$R\left(N\right)=\frac{1}{12}\underset{\mathrm{SCn}=1}{\overset{12}{\mathrm{\Σ}}}\left[\begin{array}{c}{n}_{0,\mathrm{SCn}}\left(N\right)\\ n_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]{\left[\begin{array}{c}{n}_{0,\mathrm{SCn}}\left(N\right)\\ n_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]}^H$ 计算；

其中，n_ant，SCn(N)为ant个数天线，在同一物理资源块中同一频域子载波上的干扰加噪声信号；上标H表示共轭转置。

优选地，所述步骤C通过公式：

$H_{\mathrm{SCn}}\left(N\right)=\left[\begin{array}{c}{h}_{0,\mathrm{SCn}}\left(N\right)\\ h_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]$ 计算；

其中，H_SCn(N)为1个物理资源块中12个子载波上的信道矩阵。

优选地，所述步骤D通过公式：

w_SCn(N)＝(H_SCn(N)^HR(N)^-1)^T

                          计算；

＝[w_SCn，0(N)w_SCn，1(N)]

其中，上标T表示矩阵转置操作，上标-1表示矩阵求逆操作。

优选地，所述步骤E通过公式：

α_SCn(N)＝w_SCn(N)H_SCn(N)计算。

优选地，所述步骤F通过公式：

c_SCn(N)＝w_SCn，0(N)r_0，SCn(N)+w_SCn，1(N)r_1，SCn(N)计算。

优选地，所述步骤H通过公式：

c′_SCn(N)＝c_SCn(N)/α_SCn(N)计算。

本发明还提供一种干扰抑制合并系统，包括：

第一计算单元，用于基站侧接收端的ant个数天线分别计算物理资源块中导频符号子载波上的干扰加噪声信号n_ant，SCn(N)；

注：ant为基站侧接收端接收天线编号数；SCn为频域子载波的编号数；

第二计算单元，用于根据所述第一计算单元得到的ant个数的干扰加噪声信号n_ant，SCn(N)，计算得到所述物理资源块的干扰加噪声协方差矩阵R(N)；

第三计算单元，用于计算所述物理资源块中子载波上的信道矩阵H_SCn(N)；

第四计算单元，用于根据所述第二计算单元得到的干扰加噪声协方差矩阵R(N)和第三计算单元得到的信道矩阵H_SCn(N)，计算所述物理资源块中子载波上权值合并矩阵w_SCn(N)；

第五计算单元，用于根据所述第三计算单元得到的信道矩阵H_SCn(N)和所述第四计算单元得到的上权值合并矩阵w_SCn(N)，计算所述物理资源块中子载波上合并值调整系数α_SCn(N)；

第六计算单元，用于根据所述第四计算单元计算出的权值合并矩阵w_SCn(N)，计算非导频符号的物理资源块的子载波上的接收信号合并值c_SCn(N)；

第七计算单元，用于根据所述第五计算单元计算出的所述合并值调整系数α_SCn(N)调整所述第六计算单元非导频符号的物理资源块的子载波上的接收信号合并值c_SCn(N)，得到调整结果c′_SCn(N)。

优选地，所述第一计算单元，具体计算一个物理资源块中12个子载波上的干扰加噪声信号n_ant，SCn(N)；

所述第二计算单元的具体计算公式为：

$R\left(N\right)=\frac{1}{12}\underset{\mathrm{SCn}=1}{\overset{12}{\mathrm{\Σ}}}\left[\begin{array}{c}{n}_{0,\mathrm{SCn}}\left(N\right)\\ n_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]{\left[\begin{array}{c}{n}_{0,\mathrm{SCn}}\left(N\right)\\ n_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]}^H$

其中，n_ant，SCn(N)为ant个数天线，在同一物理资源块中同一频域子载波上的干扰加噪声信号；上标H表示共轭转置。

与现有技术相比，本发明具有以下优点：

本发明实施例干扰抑制合并方法，具体通过分别计算基站侧接收端的ant个数天线接收的物理资源块中非导频符号的子载波上的干扰加噪声信号n_ant，SCn(N)；再根据所述干扰加噪声信号n_ant，SCn(N)，计算所述物理资源块的干扰加噪声协方差矩阵R(N)；然后计算所述物理资源块中子载波上的信道矩阵H_SCn(N)；再根据所述信道矩阵H_SCn(N)，计算所述物理资源块中子载波上权值合并矩阵w_SCn(N)；然后根据所述信道矩阵H_SCn(N)和所述上权值合并矩阵w_SCn(N)计算所述物理资源块中子载波上合并值调整系数α_SCn(N)；根据计算出的所述权值合并矩阵w_SCn(N)，计算非导频符号的物理资源块的子载波上的接收信号合并值c_SCn(N)；最后根据计算出的所述合并值调整系数α_SCn(N)调整非导频符号的物理资源块的子载波上的接收信号合并值c_SCn(N)，得到调整结果c′_SCn(N)。所述调整结果c′_SCn(N)可以直接传递给整个上行链路的下一级使用。

本发明实施例所述方法可以通过上行信道中连续存放的导频序列，以物理资源块为单位，解出噪声加干扰矩阵协方差矩阵，进而实现IRC方法对干扰进行抑制，即得到一套在上行信道中可行的IRC处理方法。

附图说明

图1为本发明实施例所述干扰抑制合并方法流程图；

图2为上行方向小区间用户干扰示图；

图3为上行无线帧结构和时隙PRB结构图；

图4为本发明所述方法在一个时隙编号为N的PRB上的导频符号上IRC操作流程图；

图5为本发明所述方法在一个时隙编号为N的PRB上的非导频符号上IRC操作流程图；

图6为本发明实施例所述干扰抑制合并系统结构图。

具体实施方式

本发明的目的是提供一种干扰抑制合并方法，解出噪声加干扰矩阵协方差矩阵，进而实现IRC方法对干扰进行抑制，即得到一套在上行信道中可行的IRC处理方法。

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

为了便于本领域技术人员理解本发明的应用环境，下面结合图2具体说明。

参见图2，该图为上行方向小区间用户干扰示图。

UE1(第一无线终端)在小区1和小区2边缘向基站1发送信号。上行方向UE端为1根发送天线，基站端为2根接收天线，发送方式为1x2分集发送。与此同时，在小区1和小区2之间又有1个(或多个)UE处于基站边缘，对UE1向基站1发送的部分PRB(物理资源块)进行干扰。(例如图2中，UE2处于小区1和小区2边缘对UE1进行干扰)。本发明TD-LTE(时分-长期演进)上行部分中IRC算法正是应用于这一场景。

如图2所示系统模型，IRC算法主要作用是计算基站1上2根接收天线的权值，求出权值后分别乘以2根天线各自对应子载波上的接收信号，对应子载波信号相加后即可抑制干扰。具体过程参见图1所示。

参见图1，该图为本发明实施例所述干扰抑制合并方法流程图。

本发明实施例所述干扰抑制合并方法，其特征在于，所述方法包括：

S100、计算物理资源块中导频符号的子载波上的干扰加噪声信号n_ant，SCn(N)；

分别计算在基站侧接收端的ant个数天线接收的物理资源块中导频符号的子载波上的干扰加噪声信号n_ant，SCn(N)。

参见图3，该图为上行无线帧结构和时隙PRB结构图。

根据TD-LTE协议版本ts36.211-860。一个上行传输帧由10个上行子帧构成，1个上行子帧由2个时隙构成。1个时隙包括7个(或6)传输符号，其中包含1个导频符号。每一个传输符号的1个PRB中包含12个子载波。

在上行PUSCH信道中，导频符号中导频序列是连续存放的。即在导频符号上，每1个PRB上的子载波上导频序列都存在。

步骤S100具体可以将12个载波上的接收信号r_ant，SCn(N)分别对应减去信道系数h_ant，SCn(N)对应乘以本地产生的导频序列x_SCn(N)，其结果为噪声加干扰信号n_ant，SCn(N)。

对1个PRB的12个频域载波上，将信道系数h_ant，SCn(N)对应乘以本地产生的导频序列x_SCn(N)，其中基站端本地产生的导频序列数据x_SCn(N)，在天线0与天线1产生的值相同。

步骤S100的具体公式如下：n_ant，SCn(N)＝r_ant，SCn(N)-h_ant，SCn(N)x_SCn(N)

其中，信道系数h_ant，SCn(N)，其中ant为基站端接收天线编号数，ant具体为2根，取值为0，1。

SCn为频域子载波的编号数，取值1到12。

N为TD-LTE协议中上行部分physical resource block物理资源块的编号数，其取值由用户所占PRB数所定，取值可以为1到maxN，maxN最大值为100。

h_ant，SCn(N)以子载波为单位。

S200、根据所述干扰加噪声信号n_ant，SCn(N)，计算所述物理资源块的干扰加噪声协方差矩阵R(N)；

将ant个数(2根)天线，同1个PRB，同一的频域子载波上的n_ant，SCn(N)，求解信号加噪声的协方差矩阵，其中上标H表示共轭转置(赫密特转置)：

步骤S200的具体公式如下：

$R\left(N\right)=\frac{1}{12}\underset{\mathrm{SCn}=1}{\overset{12}{\mathrm{\Σ}}}\left[\begin{array}{c}{n}_{0,\mathrm{SCn}}\left(N\right)\\ n_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]{\left[\begin{array}{c}{n}_{0,\mathrm{SCn}}\left(N\right)\\ n_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]}^H$

S300、计算所述物理资源块中子载波上的信道矩阵H_SCn(N)；

步骤S300的具体公式为：

$H_{\mathrm{SCn}}\left(N\right)=\left[\begin{array}{c}{h}_{0,\mathrm{SCn}}\left(N\right)\\ h_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]$

S400、根据所述干扰加噪声协方差矩阵R(N)和信道矩阵H_SCn(N)，计算所述物理资源块中子载波上权值合并矩阵w_SCn(N)；

步骤S400的具体公式如下：

w_SCn(N)＝(H_SCn(N)^HR(N)^-1)^T

＝[w_SCn，0(N)w_SCn，1(N)]

其中，上标T表示矩阵转置操作，-1表示矩阵求逆操作。

S500、根据所述信道矩阵H_SCn(N)和所述上权值合并矩阵w_SCn(N)计算所述物理资源块中子载波上合并值调整系数α_SCn(N)；

步骤S500的具体公式如下：

α_SCn(N)＝w_SCn(N)H_SCn(N)

步骤S100-S500的过程可以参见图4所示，该图为本发明所述方法在一个时隙编号为N的PRB上的导频符号上IRC操作流程图。

S600、根据计算出的所述权值合并矩阵w_SCn(N)，计算物理资源块中非导频符号的子载波上的接收信号合并值c_SCn(N)；

步骤S600的具体公式如下：

c_SCn(N)＝w_SCn，0(N)r_0，SCn(N)+w_SCn，1(N)r_1，SCn(N)

S700、根据计算出的所述合并值调整系数α_SCn(N)调整物理资源块中非导频符号的子载波上的接收信号合并值c_SCn(N)，得到调整结果c′_SCn(N)。

步骤S700的具体公式如下：

c′_SCn(N)＝c_SCn(N)/α_SCn(N)

S600-S700的具体过程可以参见图5所示，该图为本发明所述方法在一个时隙编号为N的PRB上的非导频符号上IRC操作流程图。

调整结果c′_SCn(N)就是通过乘以权值抑制了ICI干扰后，并且调整了系数的最终值，可以直接传递给整个上行链路的下一级使用。

本发明实施例所述干扰抑制合并方法可以以1个PRB位单位操作的，利用1个时隙中导频符号连续存放的导频序列求解出非导频符号对应子载波位置的合并权值，再乘以对应非导频符号子载波位置2天线接收到的不同信号相加后抑制干扰。

本发明的目的是提供一种干扰抑制合并系统，解出噪声加干扰矩阵协方差矩阵，进而实现IRC方法对干扰进行抑制，即得到一套在上行信道中可行的IRC处理方法。

参见图6，该图为本发明实施例所述干扰抑制合并系统结构图。

本发明实施例所述干扰抑制合并系统，包括第一计算单元11、第二计算单元12、第三计算单元13、第四计算单元14、第五计算单元15、第六计算单元16，以及第七计算单元17。

第一计算单元11，用于分别计算基站侧接收端ant个数天线接收的物理资源块中子载波上的干扰加噪声信号n_ant，SCn(N)。

第一计算单元11，具体可以分别计算基站侧接收端ant个数天线中一个物理资源块中12个子载波上的干扰加噪声信号n_ant，SCn(N)。

第一计算单元11的具体计算公式为：

n_ant，SCn(N)＝r_ant，SCn(N)-h_ant，SCn(N)x_SCn(N)

其中，r_ant，SCn(N)为载波上的接收信号；ant为基站端接收天线编号数；SCn为频域子载波的编号数；

h_ant，SCn(N)为信道系数；

x_SCn(N)为本地产生的导频序列。

所述ant取值0或1，SCn取值1至12。

第二计算单元12，用于根据所述第一计算单元11得到的干扰加噪声信号n_ant，SCn(N)，计算所述物理资源块的干扰加噪声协方差矩阵R(N)；

第二计算单元12的具体计算公式为：

$R\left(N\right)=\frac{1}{12}\underset{\mathrm{SCn}=1}{\overset{12}{\mathrm{\Σ}}}\left[\begin{array}{c}{n}_{0,\mathrm{SCn}}\left(N\right)\\ n_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]{\left[\begin{array}{c}{n}_{0,\mathrm{SCn}}\left(N\right)\\ n_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]}^H$

其中，n_ant，SCn(N)为ant个数(2根)天线，同一PRB，同一频域子载波上的干扰加噪声信号；上标H表示共轭转置。

第三计算单元13，用于计算所述物理资源块中子载波上的信道矩阵H_SCn(N)。

第三计算单元13的具体计算公式为：

$H_{\mathrm{SCn}}\left(N\right)=\left[\begin{array}{c}{h}_{0,\mathrm{SCn}}\left(N\right)\\ h_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]$

其中，H_SCn(N)为1个PRB中12个子载波上的信道矩阵。

第四计算单元14，用于根据所述第二计算单元12得到的干扰加噪声协方差矩阵R(N)和第三计算单元得到的信道矩阵H_SCn(N)，计算所述物理资源块中子载波上权值合并矩阵w_SCn(N)；

第四计算单元的具体计算公式：

w_SCn(N)＝(H_SCn(N)^HR(N)^-1)^T

＝[w_SCn，0(N)w_SCn，1(N)]

其中，上标T表示矩阵转置操作，上标-1表示矩阵求逆操作。

第五计算单元15，用于根据所述第三计算单元13得到的信道矩阵H_SCn(N)和所述第四计算单元14得到的上权值合并矩阵w_SCn(N)，计算所述物理资源块中子载波上合并值调整系数α_SCn(N)；

第五计算单元15的具体计算公式为：α_SCn(N)＝w_SCn(N)H_SCn(N)。

第六计算单元16，用于根据所述第四计算单元14计算出的权值合并矩阵w_SCn(N)，计算物理资源块中非导频符号的子载波上的接收信号合并值c_SCn(N)；

第六计算单元16的具体计算公式为：

c_SCn(N)＝w_SCn，0(N)r_0，SCn(N)+w_SCn，1(N)r_1，SCn(N)。

第七计算单元17，用于根据所述第五计算单元15计算出的所述合并值调整系数α_SCn(N)，调整所述第六计算单元16物理资源块中非导频符号的子载波上的接收信号合并值c_SCn(N)，得到调整结果c′_SCn(N)。

第七计算单元17的具体计算公式为：c′_SCn(N)＝c_SCn(N)/α_SCn(N)。

本发明实施例所述干扰抑制合并系统，通过第七计算单元17得到调整结果c_SCn(N)就是通过乘以权值抑制了ICI干扰后，并且调整了系数的最终值，可以直接传递给整个上行链路的下一级使用。

以上对本发明所提供的所述干扰抑制合并方法及系统进行了详细介绍，本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (10)

1.一种干扰抑制合并方法，其特征在于，所述方法包括：

A、分别计算在基站侧接收端的ant个数天线接收的物理资源块中导频符号子载波上的干扰加噪声信号n_ant，SCn(N)；

B、根据所述干扰加噪声信号n_ant，SCn(N)，计算得到所述物理资源块的干扰加噪声协方差矩阵R(N)；

C、计算所述物理资源块中子载波上的信道矩阵H_SCn(N)；

D、根据所述干扰加噪声协方差矩阵R(N)和信道矩阵H_SCn(N)，计算所述物理资源块中子载波上权值合并矩阵w_SCn(N)；

E、根据所述信道矩阵H_SCn(N)和所述上权值合并矩阵w_SCn(N)计算所述物理资源块中子载波上合并值调整系数α_SCn(N)；

F、根据计算出的所述权值合并矩阵w_SCn(N)，计算物理资源块中非导频符号的子载波上的接收信号合并值c_SCn(N)；

G、根据计算出的所述合并值调整系数α_SCn(N)调整物理资源块中非导频符号的子载波上的接收信号合并值c_SCn(N)，得到调整结果c′_SCn(N)。

2.根据权利要求1所述的方法，其特征在于，所述步骤A具体为：

计算一个物理资源块中12个子载波上的干扰加噪声信号n_ant，SCn(N)；

所述步骤A通过公式：

n_ant，SCn(N)＝r_ant，SCn(N)-h_ant，SCn(N)x_SCn(N)计算；

其中，r_ant，SCn(N)为载波上的接收信号；

h_ant，SCn(N)为信道系数；

x_SCn(N)为本地产生的导频序列；

所述ant取值0或1；

SCn取值1至12。

3.根据权利要求2所述的方法，其特征在于，所述步骤B通过公式：

$R\left(N\right)=\frac{1}{12}\underset{\mathrm{SCn}=1}{\overset{12}{\mathrm{\Σ}}}\left[\begin{array}{c}{n}_{0,\mathrm{SCn}}\left(N\right)\\ n_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]{\left[\begin{array}{c}{n}_{0,\mathrm{SCn}}\left(N\right)\\ n_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]}^H$ 计算；

其中，n_ant，SCn(N)为ant个数天线，在同一物理资源块中同一频域子载波上的干扰加噪声信号；上标H表示共轭转置。

4.根据权利要求3所述的方法，其特征在于，所述步骤C通过公式：

$H_{\mathrm{SCn}}\left(N\right)=\left[\begin{array}{c}{h}_{0,\mathrm{SCn}}\left(N\right)\\ h_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]$ 计算；

其中，H_SCn(N)为1个物理资源块中12个子载波上的信道矩阵。

5.根据权利要求4所述的方法，其特征在于，所述步骤D通过公式：

w_SCn(N)＝(H_SCn(N)^HR(N)^-1)^T

                                计算；

      ＝[w_SCn，0(N)w_SCn，1(N)]

其中，上标T表示矩阵转置操作，上标-1表示矩阵求逆操作。

6.根据权利要求1所述的方法，其特征在于，所述步骤E通过公式：

α_SCn(N)＝w_SCn(N)H_SCn(N)计异。

7.根据权利要求1所述的方法，其特征在于，所述步骤F通过公式：

c_SCn(N)＝w_SCn，0(N)r_0，SCn(N)+w_SCn，1(N)r_1，SCn(N)计算。

8.根据权利要求1所述的方法，其特征在于，所述步骤H通过公式：

c′_SCn(N)＝c_SCn(N)/α_SCn(N)计算。

9.一种干扰抑制合并系统，其特征在于，所述系统包括：

第一计算单元，用于基站侧接收端的ant个数天线分别计算物理资源块中导频符号子载波上的干扰加噪声信号n_ant，SCn(N)；

注：ant为基站侧接收端接收天线编号数；SCn为频域子载波的编号数；

第二计算单元，用于根据所述第一计算单元得到的ant个数的干扰加噪声信号n_ant，SCn(N)，计算得到所述物理资源块的干扰加噪声协方差矩阵R(N)；

第三计算单元，用于计算所述物理资源块中子载波上的信道矩阵H_SCn(N)；

第四计算单元，用于根据所述第二计算单元得到的干扰加噪声协方差矩阵R(N)和第三计算单元得到的信道矩阵H_SCn(N)，计算所述物理资源块中子载波上权值合并矩阵w_SCn(N)；

第五计算单元，用于根据所述第三计算单元得到的信道矩阵H_SCn(N)和所述第四计算单元得到的上权值合并矩阵w_SCn(N)，计算所述物理资源块中子载波上合并值调整系数α_SCn(N)；

第六计算单元，用于根据所述第四计算单元计算出的权值合并矩阵w_SCn(N)，计算非导频符号的物理资源块的子载波上的接收信号合并值c_SCn(N)；

第七计算单元，用于根据所述第五计算单元计算出的所述合并值调整系数α_SCn(N)调整所述第六计算单元非导频符号的物理资源块的子载波上的接收信号合并值c_SCn(N)，得到调整结果c′_SCn(N)。

10.根据权利要求9所述的系统，其特征在于，所述第一计算单元，具体计算一个物理资源块中12个子载波上的干扰加噪声信号n_ant，SCn(N)；

所述第二计算单元的具体计算公式为：

$R\left(N\right)=\frac{1}{12}\underset{\mathrm{SCn}=1}{\overset{12}{\mathrm{\Σ}}}\left[\begin{array}{c}{n}_{0,\mathrm{SCn}}\left(N\right)\\ n_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]{\left[\begin{array}{c}{n}_{0,\mathrm{SCn}}\left(N\right)\\ n_{1,\mathrm{SCn}}\left(N\right)\end{array}\right]}^H$

其中，n_ant，SCn(N)为ant个数天线，在同一物理资源块中同一频域子载波上的干扰加噪声信号；上标H表示共轭转置。

Images