CN103384226A - 一种频域均衡检测的方法和设备 - Google Patents

一种频域均衡检测的方法和设备 Download PDF

Info

Publication number
CN103384226A
CN103384226A CN2012101354273A CN201210135427A CN103384226A CN 103384226 A CN103384226 A CN 103384226A CN 2012101354273 A CN2012101354273 A CN 2012101354273A CN 201210135427 A CN201210135427 A CN 201210135427A CN 103384226 A CN103384226 A CN 103384226A
Authority
CN
China
Prior art keywords
slot
pilot
carrier
sinr
time domain
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012101354273A
Other languages
English (en)
Other versions
CN103384226B (zh
Inventor
任斌
李琼
戴晓明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Academy of Telecommunications Technology CATT
Datang Mobile Communications Equipment Co Ltd
Original Assignee
China Academy of Telecommunications Technology CATT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Academy of Telecommunications Technology CATT filed Critical China Academy of Telecommunications Technology CATT
Priority to CN201210135427.3A priority Critical patent/CN103384226B/zh
Publication of CN103384226A publication Critical patent/CN103384226A/zh
Application granted granted Critical
Publication of CN103384226B publication Critical patent/CN103384226B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Mobile Radio Communication Systems (AREA)
  • Noise Elimination (AREA)

Abstract

本发明公开了一种频域均衡检测的方法和设备,该方法包括:计算导频子载波的时域信干噪比SINR以及导频子载波的噪声干扰协方差矩阵Re_pilot;利用所述导频子载波的噪声干扰协方差矩阵Re_pilot计算噪干比NIR;利用所述时域信干噪比SINR以及所述噪干比NIR确定采用干扰抑制合并IRC检测方式或者最大比值合并MRC检测方式进行频域均衡检测。本发明实施例中,通过自适应的选择MRC检测方式或者IRC检测方式来进行频域均衡检测,使得当干扰在接收信号中起主导作用时,可以将上述干扰和噪声区分处理,能够充分保证检测性能;当干扰在接收信号中不起主导作用时,不会由于信道估计误差的引入导致计算精度下降,能够充分保证检测性能。

Description

一种频域均衡检测的方法和设备
技术领域
本发明涉及通信技术领域,尤其涉及一种频域均衡检测的方法和设备。
背景技术
(1)LTE(Long Term Evolution,长期演进)系统的上行数据信号和导频信号的时频域结构。
子帧n中PUSCH(Physical Uplink Shared Channel,物理上行共享信道)的时频域结构为:时域上占用一个子帧(1ms),频域上占用N个PRB(PhysicalResource Block,物理资源块),每个PRB在频域上占用12个SC(Subcarrier,子载波)。每个子帧内包含两列PUSCH导频,在常规CP(Cyclic Prefix,循环前缀)下,PUSCH的导频信号位于第四列SC-FDMA(SingleCarrier-Frequency Division Multiple Access,单载波频分多址)符号;在扩展CP下,导频信号位于第三列SC-FDMA符号。
PUSCH跳频具有子帧内跳频和子帧间跳频两种模式,且有以下组合:组合一:子帧间不跳频、子帧内不跳频;组合二:子帧间不跳频、子帧内跳频;组合三:子帧间跳频、子帧内不跳频;组合四:子帧间跳频、子帧内跳频。
图1和图2分别给出了常规CP条件下、组合一和组合二对应的PUSCH数据符号和导频符号的位置结构。其中,反斜线标注的两列SC-FDMA符号为导频符号,其它12列SC-FDMA符号为数据符号。
(2)MRC(Maximum Ratio Combining,最大比值合并)检测方式。
eNodeB(演进型节点B,即基站)进行上行信道估计,得到目标小区i内的目标用户(i,1)的信道响应估计值
Figure BDA00001596586000011
频域均衡检测过程中,在最小二乘(LS)-MRC准则下,发送符号的估计值为:
s ^ i = w H r = [ H ^ i , 1 H H ^ i , 1 ] - 1 H ^ i , 1 H r ;
频域均衡检测过程中,在最小均方误差(MMSE)-MRC准则下,发送符号的估计值为:
s ^ i = w H r = [ H ^ i , 1 H H ^ i , 1 + δ n 2 I N T ] - 1 H ^ i , 1 H r .
(3)IRC(Interference Rejection Combining,干扰抑制合并)检测方式。
eNodeB进行上行信道估计,得到信道响应估计值
Figure BDA00001596586000023
并进一步计算得到噪声和干扰的协方差矩阵
Figure BDA00001596586000024
频域均衡检测过程中,在最小二乘(LS)-IRC准则下,发送符号的估计值为:
s ^ i = w H r = [ H ^ i , 1 H R ^ e - 1 H ^ i , 1 ] - 1 H ^ i , 1 H R ^ e - 1 r ;
频域均衡检测过程中,在最小均方误差(MMSE)-IRC准则下,发送符号的估计值为:
s ^ i = w H r = [ H ^ i , 1 H R ^ e - 1 H ^ i , 1 + I N T ] - 1 H ^ i , 1 H R ^ e - 1 r .
(4)LTE系统上行存在小区间干扰条件下的接收信号模型。
假设eNodeB最多支持K个小区,目标小区为i,干扰小区为j(0<=j<=K-1),来自目标小区i的UE(User Equipment,用户设备)信号为期望信号,来自干扰小区j的UE信号为干扰信号;则在目标小区i内,调度了用户ID为(i,1)的UE,该UE的发送信号为
Figure BDA00001596586000027
信道响应为
假定小区半径较小,eNodeB接收到的多小区UE的信号基本同步,且CP(Cyclic Prefix,循环前缀)的长度远大于无线信道的时延扩展,则每个子载波近似经历平坦衰落,且等效的频域数学模型可以表达为:
r = H cell _ id 0 s cell _ id 0 + H cell _ id 1 s cell _ id 1 + . . . + H cell _ id ( K - 1 ) s cell _ id ( K - 1 ) + N = Hs + N ; 其中,
H=[Hcell_id0,Hcell_id1,...,Hcell_id(K-1)],Hcell_id(k)表示小区ID为k的UE的信道响应,k=0,1...K-1,且K表示一个eNodeB内部支持的小区数。
此外,s=[scell_id0,scell_id1,...,scell_id(K-1)]T,scell_id(k)表示小区ID为k的UE的发送信号,k=0,1...K-1,scell_id(k)∈Ω,Ω表示星座点的集合,(·)T表示·的转置;N表示噪声,服从均值为0、方差为δ2的复高斯分布。
进一步的,定义 H i , 1 s i , 1 = H 1 cell _ idi s 1 cell _ id i , i=0,1,...K-1, I = &Sigma; j = 0 , j &NotEqual; i K - 1 H j s j , 则接收信号表达式为:
r = H i , 1 s i , 1 + &Sigma; j = 0 , j &NotEqual; i K - 1 H j s j + N = H i , 1 s i , 1 + I + N .
在实现本发明的过程中,发明人发现现有技术中至少存在以下问题:
现有技术中只能固定的选择采用MRC检测方式或者IRC检测方式进行检测,通过对接收信号表达式的分析可知,当干扰在接收信号起主导作用时,MRC检测方式仅仅将上述干扰当作噪声来处理,势必引起检测性能急剧恶化;当干扰在接收信号中不起主导作用时,IRC检测方式由于信道估计误差的引入,将导致计算精度下降,并会导致检测性能下降。
发明内容
本发明实施例提供一种频域均衡检测的方法和设备,以自适应的选择MRC检测方式或者IRC检测方式来进行频域均衡检测,保证检测性能。
为了达到上述目的,本发明实施例提供一种频域均衡检测的方法,包括:
计算导频子载波的时域信干噪比SINR以及导频子载波的噪声干扰协方差矩阵Re_pilot;
利用所述导频子载波的噪声干扰协方差矩阵Re_pilot计算噪干比NIR;
利用所述时域信干噪比SINR以及所述噪干比NIR确定采用干扰抑制合并IRC检测方式或者最大比值合并MRC检测方式进行频域均衡检测。
本发明实施例提供一种频域均衡检测的设备,包括:
第一计算模块,用于计算导频子载波的时域信干噪比SINR;
第二计算模块,用于计算导频子载波的噪声干扰协方差矩阵Re_pilot;
第三计算模块,用于利用所述导频子载波的噪声干扰协方差矩阵Re_pilot计算噪干比NIR;
确定模块,用于利用所述时域信干噪比SINR以及所述噪干比NIR确定采用干扰抑制合并IRC检测方式或者最大比值合并MRC检测方式进行频域均衡检测。
与现有技术相比,本发明实施例至少具有以下优点:本发明实施例中,通过自适应的选择MRC检测方式或者IRC检测方式来进行频域均衡检测,使得当干扰在接收信号中起主导作用时,可以将上述干扰和噪声区分处理,能够充分保证检测性能;当干扰在接收信号中不起主导作用时,不会由于信道估计误差的引入导致计算精度下降,能够充分保证检测性能。
附图说明
为了更清楚地说明本发明的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为现有技术中子帧间不跳频、子帧内不跳频时数据符号和导频符号位置图;
图2为现有技术中子帧间不跳频、子帧内跳频时数据符号和导频符号位置图;
图3是本发明实施例一提供的一种频域均衡检测的方法流程示意图;
图4是本发明实施例二提供的一种频域均衡检测的设备结构示意图。
具体实施方式
下面将结合本发明中的附图,对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一
在LTE系统的下行采用了OFDMA(Orthogonal Frequency DivisionMultiple Access,正交频分复用多址接入),上行采用SC-FDMA技术,小区内各用户相互正交,不存在干扰,但各小区之间使用了相同频带,小区间的干扰成为主要干扰,在很多实际场景下,LTE上行是干扰受限系统,且在某些PRB带宽可能存在窄带干扰;在干扰受限系统中(系统中存在较强的用户间干扰或者频率复用引起的同道干扰),当干扰在接收信号起主导作用时,MRC检测方式将上述干扰当作噪声来处理,此时IRC检测方式的性能优于MRC检测方式;当干扰在接收信号中不起主导作用时,IRC检测方式会导致检测性能下降,此时MRC检测方式的性能优于IRC检测方式。
针对上述发现,本发明实施例一提供一种频域均衡检测的方法,该方法可以应用在LTE系统中,通过自适应的选择MRC检测方式或者IRC检测方式来进行频域均衡检测,当干扰在接收信号中起主导作用时,IRC检测方式的性能优于MRC检测方式,因此通过IRC检测方式进行频域均衡检测,可以将干扰和噪声区分处理,充分保证检测性能;当干扰在接收信号中不起主导作用时,MRC检测方式的性能优于IRC检测方式,因此通过MRC检测方式进行频域均衡检测,从而不会导致计算精度下降,充分保证检测性能。
如图3所示,该频域均衡检测的方法包括以下步骤:
步骤301,计算导频子载波的时域信干噪比SINR以及导频子载波的噪声干扰协方差矩阵Re_pilot。
步骤302,利用计算得到的导频子载波的噪声干扰协方差矩阵Re_pilot计算噪干比NIR。
步骤303,利用时域信干噪比SINR以及噪干比NIR确定采用IRC检测方式或者MRC检测方式进行频域均衡检测。
本发明实施例中,假设LTE系统上行目标小区i内有一个目标用户(i,1),并且可以定义以下变量:i表示目标小区编号,slot表示时隙编号,symbol表示SC-FDMA符号编号,sc表示子载波编号,NR表示上行接收天线个数;以及,表示目标用户(i,1)在第slot个时隙内,第sc个导频子载波上的导频信号,且其为标量;
Figure BDA00001596586000062
表示第slot个时隙内,第sc个导频子载波上的接收信号,其为NR×1维的列向量;
Figure BDA00001596586000063
表示目标小区i内,第slot个时隙内,第sc个导频子载波上的目标用户(i,1)的信道响应估计值,且其为NR×1维的列向量;
Figure BDA00001596586000064
表示第symbol个SC-FDMA符号上,第sc个导频子载波上的接收信号,且其为NR×1维的列向量。
本发明实施例中,计算导频子载波的时域信干噪比SINR的过程,具体包括:根据导频子载波的接收信号
Figure BDA00001596586000065
以及分配给目标用户的导频信号
Figure BDA00001596586000066
进行基于IDFT(Inverse Discrete Fourier Transform,离散傅里叶逆变换)和噪声干扰抑制的信道估计,以获得时域信干噪比SINR。
进一步的,上述根据导频子载波的接收信号
Figure BDA00001596586000067
以及分配给目标用户的导频信号
Figure BDA00001596586000068
进行基于IDFT和噪声干扰抑制的信道估计,获得时域信干噪比SINR,具体包括如下步骤:
步骤1、通过如下公式对
Figure BDA00001596586000069
以及
Figure BDA000015965860000610
进行基于最小二乘LS准则的信道估计(即将导频子载波的频域接收信号
Figure BDA000015965860000611
Figure BDA000015965860000612
做LS准则的信道估计),得到初步信道响应
Figure BDA000015965860000613
H ^ slot , sc i , 1 ( V 1 ) = r slot , sc i , pilot s slot , sc i , 1 ;
步骤2、通过如下公式对初步信道响应进行IDFT变换,得到第一时域信道响应
Figure BDA00001596586000073
Figure BDA00001596586000074
步骤3、通过如下公式计算第一信号功率P1(即此时的总信号功率,其单位是瓦特):
P 1 = sum { | h ^ slot , sc i , 1 ( V 1 ) | 2 } ;
步骤4、通过如下公式对第一时域信道响应
Figure BDA00001596586000076
进行时域取窗抑噪处理,得到抑噪后的第二时域信道响应
h ^ slot , sc i , 1 ( V 2 ) = F { h ^ slot , sc i , 1 ( V 1 ) } ;
其中,F{□}表示时域取窗抑制噪声的函数;
步骤5、通过如下公式计算第二信号功率P2(即抑噪处理后的总信号功率,其单位是瓦特):
P 2 = sum { | h ^ slot , sc i , 1 ( V 2 ) | 2 } ;
步骤6、通过如下公式计算时域信干噪比SINR(其单位是dB):
SINR time = 10 &times; log 10 ( P 2 P 1 - P 2 ) .
本发明实施例中,计算导频子载波的噪声干扰协方差矩阵Re_pilot的过程,包括:根据导频子载波的接收信号
Figure BDA000015965860000711
以及分配给目标用户的导频信号
Figure BDA000015965860000712
进行基于IDFT和噪声干扰抑制的信道估计,获得信道响应估计值
Figure BDA000015965860000713
利用信道响应估计值
Figure BDA000015965860000714
计算导频子载波的噪声干扰协方差矩阵Re_pilot。
进一步的,上述根据导频子载波的接收信号
Figure BDA00001596586000081
以及分配给目标用户的导频信号
Figure BDA00001596586000082
进行基于IDFT和噪声干扰抑制的信道估计,获得信道响应估计值
Figure BDA00001596586000083
具体包括如下步骤:
步骤1、通过如下公式对
Figure BDA00001596586000084
以及
Figure BDA00001596586000085
进行基于最小二乘LS准则的信道估计(即将导频子载波的频域接收信号
Figure BDA00001596586000086
Figure BDA00001596586000087
做LS准则的信道估计),得到初步信道响应
Figure BDA00001596586000088
H ^ slot , sc i , 1 ( V 1 ) = r slot , sc i , pilot s slot , sc i , 1 ;
步骤2、通过如下公式对初步信道响应
Figure BDA000015965860000810
进行IDFT变换,得到第一时域信道响应
Figure BDA000015965860000811
h ^ slot , sc i , 1 ( V 1 ) = IDFT { H ^ slot , sc i , 1 ( V 1 ) } ;
步骤3、通过如下公式对第一时域信道响应进行时域取窗抑噪处理,得到抑噪后的第二时域信道响应
Figure BDA000015965860000814
h ^ slot , sc i , 1 ( V 2 ) = F { h ^ slot , sc i , 1 ( V 1 ) } ;
F{□}表示时域取窗抑制噪声的函数;
步骤4、将
Figure BDA000015965860000816
做DFT(Discrete Fourier Transform,离散傅立叶变换)变换到频域,得到处理后的信道响应估计值
Figure BDA000015965860000817
此外,上述利用信道响应估计值计算导频子载波的噪声干扰协方差矩阵Re_pilot的过程,具体包括:对每个slot内的每个导频子载波,按照如下公式计算每个导频子载波的噪声干扰协方差矩阵
Figure BDA000015965860000819
R ^ e slot , sc i , pilot = E [ | | r slot , sc i , pilot - H ^ slot , sc i , 1 s slot , sc i , 1 | | 2 2 ] ;
以及,按照如下公式计算一个PRB内所有导频子载波对应的平均后的导频子载波的噪声干扰协方差矩阵
R ^ e n RB i , pilot = 1 N sc RB &Sigma; sc = N sc RB n RB N sc RB ( n RB + 1 ) - 1 ( R ^ e slot , sc i , pilot ) ;
需要注意的是,
Figure BDA00001596586000093
表示向量A乘以A的共轭转置,即:
Figure BDA00001596586000094
E[A]表示向量A/矩阵A的数学期望(或者平均值);表示一个PRB内包含的子载波个数;nRB表示当前PRB的编号。
以一个PRB内包含的子载波个数为12为例,则以一个PRB的12个子载波为基本单位计算得到导频子载波的噪声干扰协方差矩阵
Figure BDA00001596586000096
具体的,对每个slot内的每个导频子载波,计算每个导频子载波的噪声干扰协方差矩阵
Figure BDA00001596586000097
之后以12个子载波为基本单位对导频子载波的
Figure BDA00001596586000098
进行平均,得到平均后的
Figure BDA00001596586000099
Figure BDA000015965860000910
表示一个PRB内包含的子载波个数。
进一步的,还可以有
Figure BDA000015965860000911
的简化计算方法,由LTE导频信号特点可知:导频序列si模值为1,即||si||=1,那么:
R e = E [ | | r - H i s i | | 2 2 ] = E [ | | r - H i s i s i | | 2 2 ] = E [ | | r S i - H i | | 2 2 ]
在上述公式中,
Figure BDA000015965860000913
表示向量A乘以A的共轭转置,即:
Figure BDA000015965860000914
E[A]表示向量A/矩阵A的数学期望(或者平均值),r表示频域接收信号;即Re可以通过用户i的初始信道估计值(接收信号除以导频序列)与用户i经过抑制噪声处理后的信道估计值的差异来得到。
本发明实施例中,利用导频子载波的噪声干扰协方差矩阵Re_pilot计算噪干比NIR的过程,具体包括:按照如下公式计算噪干比NIR(即噪声功率和干扰功率的比值,简称“噪干比”):
NIR = mean ( diag ( R e ) ) mean ( abs ( R e - diag ( diag ( R e ) ) ) ) - 1 = &beta; &alpha; - 1 ;
在上述公式中,mean(A)表示向量A平均值,diag(A)表示提取矩阵A的对角元素,abs(A)表示取向量A的绝对值;以及,
β=mean(diag(Re)),且其是Re矩阵对角线元素的平均值,可以看作是干扰和噪声的平均功率;以及,
α=mean(abs(Re-diag(diag(Re)))),且其是Re矩阵非对角线元素的模值的平均值,近似为干扰平均功率;以及,
Re由导频子载波的噪声干扰协方差矩阵Re_pilot(如
Figure BDA00001596586000101
)确定。
本发明实施例中,利用时域信干噪比SINR以及噪干比NIR确定采用IRC检测方式或者MRC检测方式进行频域均衡检测的过程,具体包括:
比较时域信干噪比SINR与预设门限值SINR_th1和SINR_th2,其中,SINR_th2小于SINR_th1;
当时域信干噪比SINR大于SINR_th1时,则认为信道质量较好,确定需要在每个子载波上采用MRC检测方式进行频域均衡检测;
当时域信干噪比SINR小于SINR_th2时,则认为信道质量较差,确定需要在每个子载波上采用IRC检测方式进行频域均衡检测;
当时域信干噪比SINR位于SINR_th2与SINR_th1之间时,基于每个PRB,比较噪干比NIR与预设门限值NIR_th;
当噪干比NIR大于NIR_th时,则认为噪声功率小于干扰功率,Re仍保持不变,确定需要在每个子载波上采用IRC检测方式进行频域均衡检测;
当噪干比NIR不大于NIR_th时,则认为噪声功率大于等于干扰功率,确定需要在每个子载波上采用MRC检测方式进行频域均衡检测。
需要注意的是,SINR_th1和SINR_th2的确定方式,具体包括:通过MCS(Modulation and Coding Style,调制编码方式)等级确定SINR_th1和SINR_th2;或者,通过不同的调制阶数确定SINR_th1和SINR_th2。
通过MCS等级确定SINR_th1和SINR_th2时,SINR_th1和SINR_th2与MCS等级挂钩,对于不同的MCS等级通过仿真,设定相对合理的SINR_th1和SINR_th2;通过不同的调制阶数确定SINR_th1和SINR_th2时,与不同的调制阶数相对应,针对QPSK(Quadrature Phase Shift Keying,四相相移键控信号)、16QAM(Quadrature Amplitude Modulation,正交幅度调制)和64QAM等,分别给出一个门限SINR_th1和SINR_th2的取值范围。
本发明实施例中,采用IRC检测方式或者MRC检测方式进行频域均衡检测的过程进一步包括:利用数据子载波的接收信号
Figure BDA00001596586000111
数据子载波的信道估计响应矩阵
Figure BDA00001596586000112
以及更新后的数据子载波的噪声干扰协方差矩阵
Figure BDA00001596586000113
对每个数据子载波进行频域均衡检测。
例如,在确定采用MRC检测方式进行频域均衡检测后,利用数据子载波的接收信号
Figure BDA00001596586000114
信道估计响应矩阵
Figure BDA00001596586000115
以及更新后的数据子载波的噪声干扰协方差矩阵
Figure BDA00001596586000116
等信息,进行最小二乘(LS)-MRC准则下的频域均衡检测,或者进行最小均方误差(MMSE)-MRC准则下的频域均衡检测。
进一步的,基于是否存在PUSCH子帧内跳频,则数据子载波的信道估计响应矩阵
Figure BDA00001596586000117
的确定方式,具体为:当不存在PUSCH子帧内跳频时,对进行时隙间插值运算,得到每个SC-FDMA符号上数据子载波的信道估计响应矩阵
Figure BDA00001596586000119
或者,当存在PUSCH子帧内跳频时,则对于一个slot内的所有SC-FDMA符号,在相同频带上的数据子载波的等效信道估计响应矩阵
Figure BDA000015965860001110
等于该slot内导频子载波的
Figure BDA000015965860001111
此外,更新后的数据子载波的噪声干扰协方差矩阵的确定方式具体为:按照如下公式确定数据子载波的噪声干扰协方差矩阵
R ^ e n RB i , Data ( new ) = A 1 &times; R ^ e n RB i , Data + A 2 &times; I ;
需要注意的是,当采用IRC检测方式进行频域均衡检测时,则A1=1.0,A2=0.0;当采用MRC检测方式进行频域均衡检测时,则A1=0.0,
Figure BDA00001596586000121
且其是数据子载波的噪声干扰协方差矩阵
Figure BDA00001596586000122
矩阵对角线元素的平均值,是干扰和噪声的平均功率;此外,I表示单位矩阵。
进一步的,基于是否存在PUSCH子帧内跳频,则数据子载波的噪声干扰协方差矩阵
Figure BDA00001596586000123
的确定方式,具体为:当不存在PUSCH子帧内跳频时,需要对导频子载波的噪声干扰协方差矩阵Re_pilot进行时隙间插值运算,以得到每个SC-FDMA符号上数据子载波的噪声干扰协方差矩阵
Figure BDA00001596586000124
或者,当存在PUSCH子帧内跳频时,则对于一个slot内的SC-FDMA符号,需要在相同频带上的数据子载波的噪声干扰协方差矩阵
Figure BDA00001596586000125
等于该slot内导频子载波的噪声干扰协方差矩阵Re_pilot。
本发明实施例的一种优选实施方式中,可以根据数据子载波的接收信号
Figure BDA00001596586000126
信道估计响应矩阵
Figure BDA00001596586000127
以及更新后的数据子载波的噪声干扰协方差矩阵
Figure BDA00001596586000128
对每个数据子载波进行LS准则下的频域均衡检测;且在LS准则下,得到上行发送的PUSCH数据符号的估计值为:
s ^ symbol , sc i , data = [ ( H ^ symbol , sc i , data ) H ( R ^ e n RB i , Data ( new ) ) - 1 H ^ symbol , sc i , data ] - 1 ( H ^ symbol , sc i , data ) H ( R ^ e n RB i , Data ( new ) ) - 1 r symbol , sc i , data
本发明实施例中,在频域均衡检测过程结束之后,还可以对频域均衡检测后的信号进行解调、解扰和译码等操作,得到信源比特。
实施例二
基于与上述方法同样的发明构思,本发明实施例中还提供了一种频域均衡检测的设备,如图4所示,该设备包括:
第一计算模块11,用于计算导频子载波的时域信干噪比SINR;
第二计算模块12,用于计算导频子载波的噪声干扰协方差矩阵Re_pilot;
第三计算模块13,用于利用所述导频子载波的噪声干扰协方差矩阵Re_pilot计算噪干比NIR;
确定模块14,用于利用所述时域信干噪比SINR以及所述噪干比NIR确定采用干扰抑制合并IRC检测方式或者最大比值合并MRC检测方式进行频域均衡检测。
所述第一计算模块11,具体用于根据导频子载波的接收信号
Figure BDA00001596586000131
以及分配给目标用户的导频信号
Figure BDA00001596586000132
进行基于离散傅里叶逆变换IDFT和噪声干扰抑制的信道估计,获得所述时域信干噪比SINR;
其中,(i,1)为上行目标小区i内的目标用户,且i表示目标小区编号,slot表示时隙编号,sc表示子载波编号;
Figure BDA00001596586000133
表示第slot个时隙内,第sc个导频子载波上的接收信号,其为NR×1维的列向量,且NR表示上行接收天线个数;
Figure BDA00001596586000134
表示目标用户(i,1)在第slot个时隙内,第sc个导频子载波上的导频信号,且其为标量。
所述第一计算模块11,进一步用于通过如下公式对以及
Figure BDA00001596586000136
进行基于最小二乘LS准则的信道估计,得到初步信道响应
Figure BDA00001596586000137
H ^ slot , sc i , 1 ( V 1 ) = r slot , sc i , pilot s slot , sc i , 1 ;
通过如下公式对所述初步信道响应
Figure BDA00001596586000139
进行IDFT变换,得到第一时域信道响应
Figure BDA000015965860001310
h ^ slot , sc i , 1 ( V 1 ) = IDFT { H ^ slot , sc i , 1 ( V 1 ) } ;
通过如下公式计算第一信号功率P1
P 1 = sum { | h ^ slot , sc i , 1 ( V 1 ) | 2 } ;
通过如下公式对所述第一时域信道响应
Figure BDA000015965860001313
进行时域取窗抑噪处理,得到抑噪后的第二时域信道响应
Figure BDA000015965860001314
Figure BDA00001596586000141
F{□}表示时域取窗抑制噪声的函数;
通过如下公式计算第二信号功率P2
P 2 = sum { | h ^ slot , sc i , 1 ( V 2 ) | 2 } ;
通过如下公式计算所述时域信干噪比SINR:
SINR time = 10 &times; log 10 ( P 2 P 1 - P 2 ) .
所述第二计算模块12,具体用于根据导频子载波的接收信号
Figure BDA00001596586000144
以及分配给目标用户的导频信号
Figure BDA00001596586000145
进行基于IDFT和噪声干扰抑制的信道估计,获得信道响应估计值
Figure BDA00001596586000146
利用所述信道响应估计值
Figure BDA00001596586000147
计算所述导频子载波的噪声干扰协方差矩阵Re_pilot;
其中,(i,1)为上行目标小区i内的目标用户,且i表示目标小区编号,slot表示时隙编号,sc表示子载波编号;
表示第slot个时隙内,第sc个导频子载波上的接收信号,其为NR×1维的列向量,且NR表示上行接收天线个数;
Figure BDA00001596586000149
表示目标用户(i,1)在第slot个时隙内,第sc个导频子载波上的导频信号,且其为标量;
表示目标小区i内,第slot个时隙内,第sc个导频子载波上的目标用户(i,1)的信道响应估计值,且其为NR×1维的列向量。
所述第二计算模块12,进一步用于通过如下公式对
Figure BDA000015965860001411
以及
Figure BDA000015965860001412
进行基于最小二乘LS准则的信道估计,得到初步信道响应
Figure BDA000015965860001413
H ^ slot , sc i , 1 ( V 1 ) = r slot , sc i , pilot s slot , sc i , 1 ;
通过如下公式对所述初步信道响应
Figure BDA00001596586000151
进行IDFT变换,得到第一时域信道响应
Figure BDA00001596586000152
h ^ slot , sc i , 1 ( V 1 ) = IDFT { H ^ slot , sc i , 1 ( V 1 ) } ;
通过如下公式对所述第一时域信道响应进行时域取窗抑噪处理,得到抑噪后的第二时域信道响应
Figure BDA00001596586000155
Figure BDA00001596586000156
F{□}表示时域取窗抑制噪声的函数;
做离散傅立叶变换DFT变换到频域,得到处理后的信道响应估计值
Figure BDA00001596586000158
所述第二计算模块12,进一步用于对每个slot内的每个导频子载波,按照如下公式计算每个导频子载波的噪声干扰协方差矩阵
R ^ e slot , sc i , pilot = E [ | | r slot , sc i , pilot - H ^ slot , sc i , 1 s slot , sc i , 1 | | 2 2 ] ;
以及,按照如下公式计算一个物理资源块PRB内所有导频子载波对应的平均后的导频子载波的噪声干扰协方差矩阵
R ^ e n RB i , pilot = 1 N sc RB &Sigma; sc = N sc RB n RB N sc RB ( n RB + 1 ) - 1 ( R ^ e slot , sc i , pilot ) ;
其中,表示向量A乘以A的共轭转置,且
Figure BDA000015965860001514
E[A]表示向量A/矩阵A的数学期望或者平均值;
Figure BDA000015965860001515
表示一个PRB内包含的子载波个数;
nRB表示当前PRB的编号。
所述第三计算模块13,具体用于按照如下公式计算噪干比NIR:
NIR = mean ( diag ( R e ) ) mean ( abs ( R e - diag ( diag ( R e ) ) ) ) - 1 = &beta; &alpha; - 1 ;
其中,mean表示平均值,diag表示提取对角元素,abs表示取绝对值;
β=mean(diag(Re)),且其是Re矩阵对角线元素的平均值,是干扰和噪声的平均功率;
α=mean(abs(Re-diag(diag(Re)))),且其是Re矩阵非对角线元素的模值的平均值,近似为干扰平均功率;
所述Re由导频子载波的噪声干扰协方差矩阵Re_pilot确定。
所述确定模块14,具体用于比较所述时域信干噪比SINR与预设门限值SINR_th1和SINR_th2,且所述SINR_th2小于所述SINR_th1;
当所述时域信干噪比SINR大于所述SINR_th1时,则确定采用MRC检测方式进行频域均衡检测;
当所述时域信干噪比SINR小于所述SINR_th2时,则确定采用IRC检测方式进行频域均衡检测;
当所述时域信干噪比SINR位于所述SINR_th2与所述SINR_th1之间时,比较所述噪干比NIR与预设门限值NIR_th;
当所述噪干比NIR大于所述NIR_th时,则确定采用IRC检测方式进行频域均衡检测;
当所述噪干比NIR不大于所述NIR_th时,则确定采用MRC检测方式进行频域均衡检测。
所述确定模块14,进一步用于通过调制编码方式MCS等级确定所述SINR_th1和SINR_th2;或者,通过不同的调制阶数确定所述SINR_th1和SINR_th2。
本发明实施例中,该设备还包括:
检测模块15,用于利用数据子载波的接收信号
Figure BDA00001596586000161
数据子载波的信道估计响应矩阵
Figure BDA00001596586000162
以及更新后的数据子载波的噪声干扰协方差矩阵
Figure BDA00001596586000163
对每个数据子载波进行频域均衡检测;
其中,i表示目标小区编号,sc表示子载波编号,symbol表示单载波频分多址SC-FDMA符号编号;nRB表示当前PRB的编号;
表示第symbol个SC-FDMA符号上,第sc个导频子载波上的接收信号,且其为NR×1维的列向量。
所述确定模块14,还用于确定所述数据子载波的信道估计响应矩阵
Figure BDA00001596586000172
且所述数据子载波的信道估计响应矩阵
Figure BDA00001596586000173
的确定方式,具体为:
当不存在物理上行共享信道PUSCH子帧内跳频时,对
Figure BDA00001596586000174
进行时隙间插值运算,得到每个SC-FDMA符号上数据子载波的信道估计响应矩阵
Figure BDA00001596586000175
或者,
当存在PUSCH子帧内跳频时,则对于一个slot内的所有SC-FDMA符号,在相同频带上的数据子载波的信道估计响应矩阵
Figure BDA00001596586000176
等于该slot内导频子载波的
Figure BDA00001596586000177
其中,
Figure BDA00001596586000178
表示目标小区i内,第slot个时隙内,第sc个导频子载波上的目标用户(i,1)的信道响应估计值,且其为NR×1维的列向量。
所述确定模块14,还用于确定所述更新后的数据子载波的噪声干扰协方差矩阵
Figure BDA00001596586000179
且所述更新后的数据子载波的噪声干扰协方差矩阵的确定方式,具体为:
按照如下公式确定所述数据子载波的噪声干扰协方差矩阵
Figure BDA000015965860001711
R ^ e n RB i , Data ( new ) = A 1 &times; R ^ e n RB i , Data + A 2 &times; I ;
其中,当采用IRC检测方式进行频域均衡检测时,则A1=1.0,A2=0.0;
当采用MRC检测方式进行频域均衡检测时,则A1=0.0,
Figure BDA000015965860001713
且其是数据子载波的噪声干扰协方差矩阵
Figure BDA000015965860001714
矩阵对角线元素的平均值;I表示单位矩阵。
所述确定模块14,还用于确定所述数据子载波的噪声干扰协方差矩阵
Figure BDA000015965860001715
且所述数据子载波的噪声干扰协方差矩阵
Figure BDA000015965860001716
的确定方式,具体为:
当不存在PUSCH子帧内跳频时,对所述导频子载波的噪声干扰协方差矩阵Re_pilot进行时隙间插值运算,得到每个SC-FDMA符号上数据子载波的噪声干扰协方差矩阵
Figure BDA00001596586000181
或者,
当存在PUSCH子帧内跳频时,则对于一个slot内的所有SC-FDMA符号,在相同频带上的数据子载波的噪声干扰协方差矩阵
Figure BDA00001596586000182
等于该slot内导频子载波的噪声干扰协方差矩阵Re_pilot。
其中,本发明装置的各个模块可以集成于一体,也可以分离部署。上述模块可以合并为一个模块,也可以进一步拆分成多个子模块。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到本发明可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。
本领域技术人员可以理解附图只是一个优选实施例的示意图,附图中的模块或流程并不一定是实施本发明所必须的。
本领域技术人员可以理解实施例中的装置中的模块可以按照实施例描述进行分布于实施例的装置中,也可以进行相应变化位于不同于本实施例的一个或多个装置中。上述实施例的模块可以合并为一个模块,也可以进一步拆分成多个子模块。
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
以上公开的仅为本发明的几个具体实施例,但是,本发明并非局限于此,任何本领域的技术人员能思之的变化都应落入本发明的保护范围。

Claims (26)

1.一种频域均衡检测的方法,其特征在于,包括:
计算导频子载波的时域信干噪比SINR以及导频子载波的噪声干扰协方差矩阵Re_pilot;
利用所述导频子载波的噪声干扰协方差矩阵Re_pilot计算噪干比NIR;
利用所述时域信干噪比SINR以及所述噪干比NIR确定采用干扰抑制合并IRC检测方式或者最大比值合并MRC检测方式进行频域均衡检测。
2.如权利要求1所述的方法,其特征在于,所述计算导频子载波的时域信干噪比SINR,包括:
根据导频子载波的接收信号
Figure FDA00001596585900011
以及分配给目标用户的导频信号进行基于离散傅里叶逆变换IDFT和噪声干扰抑制的信道估计,获得所述时域信干噪比SINR;
其中,(i,1)为上行目标小区i内的目标用户,且i表示目标小区编号,slot表示时隙编号,sc表示子载波编号;
Figure FDA00001596585900013
表示第slot个时隙内,第sc个导频子载波上的接收信号,其为NR×1维的列向量,且NR表示上行接收天线个数;
Figure FDA00001596585900014
表示目标用户(i,1)在第slot个时隙内,第sc个导频子载波上的导频信号,且其为标量。
3.如权利要求2所述的方法,其特征在于,根据导频子载波的接收信号
Figure FDA00001596585900015
以及分配给目标用户的导频信号
Figure FDA00001596585900016
进行基于离散傅里叶逆变换IDFT和噪声干扰抑制的信道估计,获得所述时域信干噪比SINR,具体包括:
通过如下公式对
Figure FDA00001596585900017
以及
Figure FDA00001596585900018
进行基于最小二乘LS准则的信道估计,得到初步信道响应
Figure FDA00001596585900019
H ^ slot , sc i , 1 ( V 1 ) = r slot , sc i , pilot s slot , sc i , 1 ;
通过如下公式对所述初步信道响应
Figure FDA00001596585900021
进行IDFT变换,得到第一时域信道响应
Figure FDA00001596585900022
h ^ slot , sc i , 1 ( V 1 ) = IDFT { H ^ slot , sc i , 1 ( V 1 ) } ;
通过如下公式计算第一信号功率P1
P 1 = sum { | h ^ slot , sc i , 1 ( V 1 ) | 2 } ;
通过如下公式对所述第一时域信道响应
Figure FDA00001596585900025
进行时域取窗抑噪处理,得到抑噪后的第二时域信道响应
Figure FDA00001596585900026
h ^ slot , sc i , 1 ( V 2 ) = F { h ^ slot , sc i , 1 ( V 1 ) } , F{□}表示时域取窗抑制噪声的函数;
通过如下公式计算第二信号功率P2
P 2 = sum { | h ^ slot , sc i , 1 ( V 2 ) | 2 } ;
通过如下公式计算所述时域信干噪比SINR:
SINR time = 10 &times; log 10 ( P 2 P 1 - P 2 ) .
4.如权利要求1所述的方法,其特征在于,所述计算导频子载波的噪声干扰协方差矩阵Re_pilot,包括:
根据导频子载波的接收信号
Figure FDA000015965859000210
以及分配给目标用户的导频信号
Figure FDA000015965859000211
进行基于IDFT和噪声干扰抑制的信道估计,获得信道响应估计值
Figure FDA000015965859000212
利用所述信道响应估计值
Figure FDA000015965859000213
计算所述导频子载波的噪声干扰协方差矩阵Re_pilot;
其中,(i,1)为上行目标小区i内的目标用户,且i表示目标小区编号,slot表示时隙编号,sc表示子载波编号;
Figure FDA000015965859000214
表示第slot个时隙内,第sc个导频子载波上的接收信号,其为NR×1维的列向量,且NR表示上行接收天线个数;
Figure FDA00001596585900031
表示目标用户(i,1)在第slot个时隙内,第sc个导频子载波上的导频信号,且其为标量;
表示目标小区i内,第slot个时隙内,第sc个导频子载波上的目标用户(i,1)的信道响应估计值,且其为NR×1维的列向量。
5.如权利要求4所述的方法,其特征在于,根据导频子载波的接收信号以及分配给目标用户的导频信号进行基于IDFT和噪声干扰抑制的信道估计,获得信道响应估计值具体包括:
通过如下公式对
Figure FDA00001596585900036
以及
Figure FDA00001596585900037
进行基于最小二乘LS准则的信道估计,得到初步信道响应
Figure FDA00001596585900038
H ^ slot , sc i , 1 ( V 1 ) = r slot , sc i , pilot s slot , sc i , 1 ;
通过如下公式对所述初步信道响应
Figure FDA000015965859000310
进行IDFT变换,得到第一时域信道响应
Figure FDA000015965859000311
h ^ slot , sc i , 1 ( V 1 ) = IDFT { H ^ slot , sc i , 1 ( V 1 ) } ;
通过如下公式对所述第一时域信道响应
Figure FDA000015965859000313
进行时域取窗抑噪处理,得到抑噪后的第二时域信道响应
Figure FDA000015965859000314
Figure FDA000015965859000315
F{□}表示时域取窗抑制噪声的函数;
Figure FDA000015965859000316
做离散傅立叶变换DFT变换到频域,得到处理后的信道响应估计值
Figure FDA000015965859000317
6.如权利要求4所述的方法,其特征在于,利用所述信道响应估计值
Figure FDA000015965859000318
计算所述导频子载波的噪声干扰协方差矩阵Re_pilot,具体包括:
对每个slot内的每个导频子载波,按照如下公式计算每个导频子载波的噪声干扰协方差矩阵
Figure FDA00001596585900041
R ^ e slot , sc i , pilot = E [ | | r slot , sc i , pilot - H ^ slot , sc i , 1 s slot , sc i , 1 | | 2 2 ] ;
以及,按照如下公式计算一个物理资源块PRB内所有导频子载波对应的平均后的导频子载波的噪声干扰协方差矩阵
Figure FDA00001596585900043
R ^ e n RB i , pilot = 1 N sc RB &Sigma; sc = N sc RB n RB N sc RB ( n RB + 1 ) - 1 ( R ^ e slot , sc i , pilot ) ;
其中,表示向量A乘以A的共轭转置,且
Figure FDA00001596585900046
E[A]表示向量A/矩阵A的数学期望或者平均值;
Figure FDA00001596585900047
表示一个PRB内包含的子载波个数;
nRB表示当前PRB的编号。
7.如权利要求1所述的方法,其特征在于,利用所述导频子载波的噪声干扰协方差矩阵Re_pilot计算噪干比NIR,包括:
按照如下公式计算噪干比NIR:
NIR = mean ( diag ( R e ) ) mean ( abs ( R e - diag ( diag ( R e ) ) ) ) - 1 = &beta; &alpha; - 1 ;
其中,mean表示平均值,diag表示提取对角元素,abs表示取绝对值;
β=mean(diag(Re)),且其是Re矩阵对角线元素的平均值,是干扰和噪声的平均功率;
α=mean(abs(Re-diag(diag(Re)))),且其是Re矩阵非对角线元素的模值的平均值,近似为干扰平均功率;
所述Re由导频子载波的噪声干扰协方差矩阵Re_pilot确定。
8.如权利要求1所述的方法,其特征在于,利用所述时域信干噪比SINR以及所述噪干比NIR确定采用干扰抑制合并IRC检测方式或者最大比值合并MRC检测方式进行频域均衡检测,包括:
比较所述时域信干噪比SINR与预设门限值SINR_th1和SINR_th2,且所述SINR_th2小于所述SINR_th1;
当所述时域信干噪比SINR大于所述SINR_th1时,则确定采用MRC检测方式进行频域均衡检测;
当所述时域信干噪比SINR小于所述SINR_th2时,则确定采用IRC检测方式进行频域均衡检测;
当所述时域信干噪比SINR位于所述SINR_th2与所述SINR_th1之间时,比较所述噪干比NIR与预设门限值NIR_th;
当所述噪干比NIR大于所述NIR_th时,则确定采用IRC检测方式进行频域均衡检测;
当所述噪干比NIR不大于所述NIR_th时,则确定采用MRC检测方式进行频域均衡检测。
9.如权利要求8所述的方法,其特征在于,所述SINR_th1和SINR_th2的确定方式,具体包括:
通过调制编码方式MCS等级确定所述SINR_th1和SINR_th2;或者,通过不同的调制阶数确定所述SINR_th1和SINR_th2。
10.如权利要求1所述的方法,其特征在于,采用IRC检测方式或者MRC检测方式进行频域均衡检测的过程,进一步包括:
利用数据子载波的接收信号
Figure FDA00001596585900051
数据子载波的信道估计响应矩阵
Figure FDA00001596585900052
以及更新后的数据子载波的噪声干扰协方差矩阵
Figure FDA00001596585900053
对每个数据子载波进行频域均衡检测;
其中,i表示目标小区编号,sc表示子载波编号,symbol表示单载波频分多址SC-FDMA符号编号;nRB表示当前PRB的编号;
Figure FDA00001596585900054
表示第symbol个SC-FDMA符号上,第sc个导频子载波上的接收信号,且其为NR×1维的列向量。
11.如权利要求10所述的方法,其特征在于,所述数据子载波的信道估计响应矩阵
Figure FDA00001596585900061
的确定方式,具体为:
当不存在物理上行共享信道PUSCH子帧内跳频时,对
Figure FDA00001596585900062
进行时隙间插值运算,得到每个SC-FDMA符号上数据子载波的信道估计响应矩阵
Figure FDA00001596585900063
或者,
当存在PUSCH子帧内跳频时,则对于一个slot内的所有SC-FDMA符号,在相同频带上的数据子载波的信道估计响应矩阵
Figure FDA00001596585900064
等于该slot内导频子载波的
Figure FDA00001596585900065
其中,表示目标小区i内,第slot个时隙内,第sc个导频子载波上的目标用户(i,1)的信道响应估计值,且其为NR×1维的列向量。
12.如权利要求10所述的方法,其特征在于,所述更新后的数据子载波的噪声干扰协方差矩阵
Figure FDA00001596585900067
的确定方式,具体为:
按照如下公式确定所述数据子载波的噪声干扰协方差矩阵
Figure FDA00001596585900068
R ^ e n RB i , Data ( new ) = A 1 &times; R ^ e n RB i , Data + A 2 &times; I ;
其中,当采用IRC检测方式进行频域均衡检测时,则A1=1.0,A2=0.0;
当采用MRC检测方式进行频域均衡检测时,则A1=0.0,
Figure FDA000015965859000610
且其是数据子载波的噪声干扰协方差矩阵
Figure FDA000015965859000611
矩阵对角线元素的平均值;I表示单位矩阵。
13.如权利要求12所述的方法,其特征在于,所述数据子载波的噪声干扰协方差矩阵
Figure FDA000015965859000612
的确定方式,具体为:
当不存在PUSCH子帧内跳频时,对所述导频子载波的噪声干扰协方差矩阵Re_pilot进行时隙间插值运算,得到每个SC-FDMA符号上数据子载波的噪声干扰协方差矩阵
Figure FDA000015965859000613
或者,
当存在PUSCH子帧内跳频时,则对于一个slot内的所有SC-FDMA符号,在相同频带上的数据子载波的噪声干扰协方差矩阵
Figure FDA00001596585900071
等于该slot内导频子载波的噪声干扰协方差矩阵Re_pilot。
14.一种频域均衡检测的设备,其特征在于,包括:
第一计算模块,用于计算导频子载波的时域信干噪比SINR;
第二计算模块,用于计算导频子载波的噪声干扰协方差矩阵Re_pilot;
第三计算模块,用于利用所述导频子载波的噪声干扰协方差矩阵Re_pilot计算噪干比NIR;
确定模块,用于利用所述时域信干噪比SINR以及所述噪干比NIR确定采用干扰抑制合并IRC检测方式或者最大比值合并MRC检测方式进行频域均衡检测。
15.如权利要求14所述的设备,其特征在于,
所述第一计算模块,具体用于根据导频子载波的接收信号
Figure FDA00001596585900072
以及分配给目标用户的导频信号
Figure FDA00001596585900073
进行基于离散傅里叶逆变换IDFT和噪声干扰抑制的信道估计,获得所述时域信干噪比SINR;
其中,(i,1)为上行目标小区i内的目标用户,且i表示目标小区编号,slot表示时隙编号,sc表示子载波编号;
Figure FDA00001596585900074
表示第slot个时隙内,第sc个导频子载波上的接收信号,其为NR×1维的列向量,且NR表示上行接收天线个数;
Figure FDA00001596585900075
表示目标用户(i,1)在第slot个时隙内,第sc个导频子载波上的导频信号,且其为标量。
16.如权利要求15所述的设备,其特征在于,
所述第一计算模块,进一步用于通过如下公式对
Figure FDA00001596585900076
以及
Figure FDA00001596585900077
进行基于最小二乘LS准则的信道估计,得到初步信道响应
Figure FDA00001596585900078
H ^ slot , sc i , 1 ( V 1 ) = r slot , sc i , pilot s slot , sc i , 1 ;
通过如下公式对所述初步信道响应
Figure FDA00001596585900082
进行IDFT变换,得到第一时域信道响应
Figure FDA00001596585900083
h ^ slot , sc i , 1 ( V 1 ) = IDFT { H ^ slot , sc i , 1 ( V 1 ) } ;
通过如下公式计算第一信号功率P1
P 1 = sum { | h ^ slot , sc i , 1 ( V 1 ) | 2 } ;
通过如下公式对所述第一时域信道响应
Figure FDA00001596585900086
进行时域取窗抑噪处理,得到抑噪后的第二时域信道响应
Figure FDA00001596585900087
Figure FDA00001596585900088
F{□}表示时域取窗抑制噪声的函数;
通过如下公式计算第二信号功率P2
P 2 = sum { | h ^ slot , sc i , 1 ( V 2 ) | 2 } ;
通过如下公式计算所述时域信干噪比SINR:
SINR time = 10 &times; log 10 ( P 2 P 1 - P 2 ) .
17.如权利要求14所述的设备,其特征在于,
所述第二计算模块,具体用于根据导频子载波的接收信号
Figure FDA000015965859000811
以及分配给目标用户的导频信号
Figure FDA000015965859000812
进行基于IDFT和噪声干扰抑制的信道估计,获得信道响应估计值
Figure FDA000015965859000813
利用所述信道响应估计值
Figure FDA000015965859000814
计算所述导频子载波的噪声干扰协方差矩阵Re_pilot;
其中,(i,1)为上行目标小区i内的目标用户,且i表示目标小区编号,slot表示时隙编号,sc表示子载波编号;
Figure FDA00001596585900091
表示第slot个时隙内,第sc个导频子载波上的接收信号,其为NR×1维的列向量,且NR表示上行接收天线个数;
Figure FDA00001596585900092
表示目标用户(i,1)在第slot个时隙内,第sc个导频子载波上的导频信号,且其为标量;
Figure FDA00001596585900093
表示目标小区i内,第slot个时隙内,第sc个导频子载波上的目标用户(i,1)的信道响应估计值,且其为NR×1维的列向量。
18.如权利要求17所述的设备,其特征在于,
所述第二计算模块,进一步用于通过如下公式对以及
Figure FDA00001596585900095
进行基于最小二乘LS准则的信道估计,得到初步信道响应
Figure FDA00001596585900096
H ^ slot , sc i , 1 ( V 1 ) = r slot , sc i , pilot s slot , sc i , 1 ;
通过如下公式对所述初步信道响应
Figure FDA00001596585900098
进行IDFT变换,得到第一时域信道响应
Figure FDA00001596585900099
h ^ slot , sc i , 1 ( V 1 ) = IDFT { H ^ slot , sc i , 1 ( V 1 ) } ;
通过如下公式对所述第一时域信道响应
Figure FDA000015965859000911
进行时域取窗抑噪处理,得到抑噪后的第二时域信道响应
Figure FDA000015965859000912
Figure FDA000015965859000913
F{□}表示时域取窗抑制噪声的函数;
Figure FDA000015965859000914
做离散傅立叶变换DFT变换到频域,得到处理后的信道响应估计值
Figure FDA000015965859000915
19.如权利要求17所述的设备,其特征在于,
所述第二计算模块,进一步用于对每个slot内的每个导频子载波,按照
如下公式计算每个导频子载波的噪声干扰协方差矩阵
Figure FDA00001596585900101
R ^ e slot , sc i , pilot = E [ | | r slot , sc i , pilot - H ^ slot , sc i , 1 s slot , sc i , 1 | | 2 2 ] ;
以及,按照如下公式计算一个物理资源块PRB内所有导频子载波对应的平均后的导频子载波的噪声干扰协方差矩阵
Figure FDA00001596585900103
R ^ e n RB i , pilot = 1 N sc RB &Sigma; sc = N sc RB n RB N sc RB ( n RB + 1 ) - 1 ( R ^ e slot , sc i , pilot ) ;
其中,
Figure FDA00001596585900105
表示向量A乘以A的共轭转置,且
Figure FDA00001596585900106
E[A]表示向量A/矩阵A的数学期望或者平均值;
表示一个PRB内包含的子载波个数;
nRB表示当前PRB的编号。
20.如权利要求14所述的设备,其特征在于,
所述第三计算模块,具体用于按照如下公式计算噪干比NIR:
NIR = mean ( diag ( R e ) ) mean ( abs ( R e - diag ( diag ( R e ) ) ) ) - 1 = &beta; &alpha; - 1 ;
其中,mean表示平均值,diag表示提取对角元素,abs表示取绝对值;
β=mean(diag(Re)),且其是Re矩阵对角线元素的平均值,是干扰和噪声的平均功率;
α=mean(abs(Re-diag(diag(Re)))),且其是Re矩阵非对角线元素的模值的平均值,近似为干扰平均功率;
所述Re由导频子载波的噪声干扰协方差矩阵Re_pilot确定。
21.如权利要求14所述的设备,其特征在于,
所述确定模块,具体用于比较所述时域信干噪比SINR与预设门限值SINR_th1和SINR_th2,且所述SINR_th2小于所述SINR_th1;
当所述时域信干噪比SINR大于所述SINR_th1时,则确定采用MRC检测方式进行频域均衡检测;
当所述时域信干噪比SINR小于所述SINR_th2时,则确定采用IRC检测方式进行频域均衡检测;
当所述时域信干噪比SINR位于所述SINR_th2与所述SINR_th1之间时,比较所述噪干比NIR与预设门限值NIR_th;
当所述噪干比NIR大于所述NIR_th时,则确定采用IRC检测方式进行频域均衡检测;
当所述噪干比NIR不大于所述NIR_th时,则确定采用MRC检测方式进行频域均衡检测。
22.如权利要求21所述的设备,其特征在于,
所述确定模块,进一步用于通过调制编码方式MCS等级确定所述SINR_th1和SINR_th2;或者,通过不同的调制阶数确定所述SINR_th1和SINR_th2。
23.如权利要求14所述的设备,其特征在于,还包括:
检测模块,用于利用数据子载波的接收信号
Figure FDA00001596585900111
数据子载波的信道估计响应矩阵
Figure FDA00001596585900112
以及更新后的数据子载波的噪声干扰协方差矩阵对每个数据子载波进行频域均衡检测;
其中,i表示目标小区编号,sc表示子载波编号,symbol表示单载波频分多址SC-FDMA符号编号;nRB表示当前PRB的编号;
Figure FDA00001596585900114
表示第symbol个SC-FDMA符号上,第sc个导频子载波上的接收信号,且其为NR×1维的列向量。
24.如权利要求23所述的设备,其特征在于,
所述确定模块,还用于确定所述数据子载波的信道估计响应矩阵且所述数据子载波的信道估计响应矩阵
Figure FDA00001596585900116
的确定方式,具体为:
当不存在物理上行共享信道PUSCH子帧内跳频时,对
Figure FDA00001596585900117
进行时隙间插值运算,得到每个SC-FDMA符号上数据子载波的信道估计响应矩阵或者,
当存在PUSCH子帧内跳频时,则对于一个slot内的所有SC-FDMA符号,在相同频带上的数据子载波的信道估计响应矩阵等于该slot内导频子载波的
其中,
Figure FDA00001596585900124
表示目标小区i内,第slot个时隙内,第sc个导频子载波上的目标用户(i,1)的信道响应估计值,且其为NR×1维的列向量。
25.如权利要求23所述的设备,其特征在于,
所述确定模块,还用于确定所述更新后的数据子载波的噪声干扰协方差矩阵
Figure FDA00001596585900125
且所述更新后的数据子载波的噪声干扰协方差矩阵
Figure FDA00001596585900126
的确定方式,具体为:
按照如下公式确定所述数据子载波的噪声干扰协方差矩阵
R ^ e n RB i , Data ( new ) = A 1 &times; R ^ e n RB i , Data + A 2 &times; I ;
其中,当采用IRC检测方式进行频域均衡检测时,则A1=1.0,A2=0.0;
当采用MRC检测方式进行频域均衡检测时,则A1=0.0,
Figure FDA00001596585900129
且其是数据子载波的噪声干扰协方差矩阵
Figure FDA000015965859001210
矩阵对角线元素的平均值;I表示单位矩阵。
26.如权利要求25所述的设备,其特征在于,
所述确定模块,还用于确定所述数据子载波的噪声干扰协方差矩阵
Figure FDA000015965859001211
且所述数据子载波的噪声干扰协方差矩阵
Figure FDA000015965859001212
的确定方式,具体为:
当不存在PUSCH子帧内跳频时,对所述导频子载波的噪声干扰协方差矩阵Re_pilot进行时隙间插值运算,得到每个SC-FDMA符号上数据子载波的噪声干扰协方差矩阵
Figure FDA000015965859001213
或者,
当存在PUSCH子帧内跳频时,则对于一个slot内的所有SC-FDMA符号,在相同频带上的数据子载波的噪声干扰协方差矩阵
Figure FDA00001596585900131
等于该slot内导频子载波的噪声干扰协方差矩阵Re_pilot。
CN201210135427.3A 2012-05-02 2012-05-02 一种频域均衡检测的方法和设备 Active CN103384226B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210135427.3A CN103384226B (zh) 2012-05-02 2012-05-02 一种频域均衡检测的方法和设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210135427.3A CN103384226B (zh) 2012-05-02 2012-05-02 一种频域均衡检测的方法和设备

Publications (2)

Publication Number Publication Date
CN103384226A true CN103384226A (zh) 2013-11-06
CN103384226B CN103384226B (zh) 2016-06-08

Family

ID=49491924

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210135427.3A Active CN103384226B (zh) 2012-05-02 2012-05-02 一种频域均衡检测的方法和设备

Country Status (1)

Country Link
CN (1) CN103384226B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104539363A (zh) * 2015-01-15 2015-04-22 东南大学 一种多色led可见光通信自适应传输方法
CN105450562A (zh) * 2014-09-23 2016-03-30 飞思卡尔半导体公司 用于联合均衡的均衡器
CN105577296A (zh) * 2014-10-17 2016-05-11 联芯科技有限公司 一种解调信干噪比测量方法及装置
CN105827559A (zh) * 2015-01-08 2016-08-03 中兴通讯股份有限公司 天线合并算法的选择方法及装置
CN106034091A (zh) * 2015-03-13 2016-10-19 电信科学技术研究院 一种检测方法及装置
CN108156107A (zh) * 2016-12-06 2018-06-12 普天信息技术有限公司 线性均衡算法的选择方法及装置
CN111465051A (zh) * 2020-03-25 2020-07-28 紫光展锐(重庆)科技有限公司 Sinr的估计方法、系统、电子设备和存储介质
CN116418656A (zh) * 2023-06-06 2023-07-11 东方空间技术(北京)有限公司 一种数据通信方法及设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080231500A1 (en) * 2007-02-23 2008-09-25 Nokia Corporation Receiving method and receiver
CN101667845A (zh) * 2008-09-05 2010-03-10 中兴通讯股份有限公司 多路信号自适应合并方法及装置
CN102215193A (zh) * 2011-05-26 2011-10-12 中兴通讯股份有限公司 频域均衡方法及装置
CN102324960A (zh) * 2011-05-13 2012-01-18 中兴通讯股份有限公司 一种干扰抑制合并方法和接收机

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080231500A1 (en) * 2007-02-23 2008-09-25 Nokia Corporation Receiving method and receiver
CN101667845A (zh) * 2008-09-05 2010-03-10 中兴通讯股份有限公司 多路信号自适应合并方法及装置
CN102324960A (zh) * 2011-05-13 2012-01-18 中兴通讯股份有限公司 一种干扰抑制合并方法和接收机
CN102215193A (zh) * 2011-05-26 2011-10-12 中兴通讯股份有限公司 频域均衡方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
TECHNICAL SPECIFICATION GROUP GSM/EDGE RADIO ACCESS NETWORK: "Feasibility study for evolved GSM/EDGE Radio Access Network", 《3GPP TR 45.912 V7.0.0》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105450562A (zh) * 2014-09-23 2016-03-30 飞思卡尔半导体公司 用于联合均衡的均衡器
CN105577296B (zh) * 2014-10-17 2018-06-22 联芯科技有限公司 一种解调信干噪比测量方法及装置
CN105577296A (zh) * 2014-10-17 2016-05-11 联芯科技有限公司 一种解调信干噪比测量方法及装置
CN105827559A (zh) * 2015-01-08 2016-08-03 中兴通讯股份有限公司 天线合并算法的选择方法及装置
CN104539363B (zh) * 2015-01-15 2017-04-19 东南大学 一种多色led可见光通信自适应传输方法
CN104539363A (zh) * 2015-01-15 2015-04-22 东南大学 一种多色led可见光通信自适应传输方法
CN106034091A (zh) * 2015-03-13 2016-10-19 电信科学技术研究院 一种检测方法及装置
CN108156107A (zh) * 2016-12-06 2018-06-12 普天信息技术有限公司 线性均衡算法的选择方法及装置
CN108156107B (zh) * 2016-12-06 2021-03-05 普天信息技术有限公司 线性均衡算法的选择方法及装置
CN111465051A (zh) * 2020-03-25 2020-07-28 紫光展锐(重庆)科技有限公司 Sinr的估计方法、系统、电子设备和存储介质
CN111465051B (zh) * 2020-03-25 2022-05-24 紫光展锐(重庆)科技有限公司 Sinr的估计方法、系统、电子设备和存储介质
CN116418656A (zh) * 2023-06-06 2023-07-11 东方空间技术(北京)有限公司 一种数据通信方法及设备
CN116418656B (zh) * 2023-06-06 2023-08-18 东方空间技术(北京)有限公司 一种数据通信方法及设备

Also Published As

Publication number Publication date
CN103384226B (zh) 2016-06-08

Similar Documents

Publication Publication Date Title
CN103384226B (zh) 一种频域均衡检测的方法和设备
CN103379078B (zh) 一种频域均衡检测的方法和设备
KR100996535B1 (ko) 다중입출력 안테나 직교주파수 분할 다중화 시스템의 셀간 간섭 완화를 위한 공간 공분산 행렬 추정 기법을 이용한 셀간 간섭 완화 방법 및 이를 이용한 셀간 간섭 완화 수신 장치
US9363817B2 (en) Mutual information based signal to interference plus noise ratio estimator for radio link monitoring
CN103155502A (zh) 干扰信号参数估计方法和装置
US20100260251A1 (en) Radio communication system, reception device, and reception method
US9036722B2 (en) Robust channel estimation and interpolation
CN101989968A (zh) 频偏估计方法及装置
EP2163057B1 (en) Memory-saving method for generating soft bit values from an ofdm signal
US10148463B2 (en) Method and device for estimating frequency offset
US20100266065A1 (en) Radio communication system, reception device, and reception method
US10129063B2 (en) Method and device for detecting signal of LTE uplink system in interference condition
CN106161322B (zh) 一种信号处理方法、上行资源分配方法及其装置
CN102480444A (zh) 一种宽带同频干扰噪声估计和干扰抑制的方法及相应系统
EP3993332A1 (en) Method and apparatus for determining doppler shift
CN102025426A (zh) 正交频分复用系统中估计载波干扰噪声比的方法及装置
CN103916229B (zh) 一种干扰抑制方法及装置
CN102111363B (zh) 一种信道估计的方法和装置
CN102932290A (zh) Lte系统干扰抑制接收方法及装置
CN103997473A (zh) 一种信号干扰的滤波方法及相关装置
US20120057543A1 (en) Wireless communication system, reception apparatus, transmission apparatus, communication method of wireless communication system, control program, and autonomous distributed network
CN102546483A (zh) 一种宽带同频干扰噪声估计和干扰抑制的方法及相应系统
CN102594737A (zh) 一种邻区干扰检测方法及系统
US20170223636A1 (en) Communications system and communications method
US10009076B2 (en) Method and apparatus for obtaining downlink data in a massive MIMO system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee after: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

Address before: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee before: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210608

Address after: 100085 1st floor, building 1, yard 5, Shangdi East Road, Haidian District, Beijing

Patentee after: DATANG MOBILE COMMUNICATIONS EQUIPMENT Co.,Ltd.

Address before: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee before: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY