发明内容
为解决上述技术问题,本申请实施方式提供了一种基于CP的频偏估计方法及其装置,以提高频偏估计的准确性。
本申请实施方式提供的基于CP的频偏估计方法包括:
在时域抽样提取一个符号的CP(循环前缀)信息序列以及位于该符号末尾处与该CP信息序列对应的第二信息序列;
对CP信息序列和第二信息序列分别补零,使补零后的信息序列长度为一个符号的长度;
将补零后的CP信息序列和第二信息序列变换到频域,得到CP信息序列频谱和第二信息序列频谱;
对CP信息序列频谱和第二信息序列频谱进行频谱扩展抑制,得到修正后的CP信息序列频谱和第二信息序列频谱;
根据修正后的CP信息序列频谱和第二信息序列频谱估计频偏。
优选地,所述对CP信息序列频谱和第二信息序列频谱进行频谱扩展抑制得到修正后的CP信息序列频谱和第二信息序列频谱具体包括:
去除CP信息序列频谱和第二信息序列频谱的带宽边缘处的子载波,得到修正后的CP信息序列频谱和第二信息序列频谱。
进一步优选地,所述去除CP信息序列频谱和第二信息序列频谱的带宽边缘处的子载波具体包括:
根据当前UE带宽大小和/或当前小区与其它小区的相邻方式确定待去除的带宽边缘处的子载波数;
根据所述子载波数去除CP信息序列频谱和第二信息序列频谱的带宽边缘处的子载波。
优选地,在对一个子帧内的N个符号进行频偏估计后,对估计得到的N个频偏值进行加权平均,将加权平均后得到的频偏作为最终的频偏,所述N大于1。
本申请实施方式还提供了一种基于CP的频偏估计装置。该装置包括:序列提取单元、序列补零单元、频域变换单元、频谱抑制单元和频偏估计单元,其中:
所述序列提取单元,用于在时域抽样提取一个符号的CP(循环前缀)信息序列以及位于该符号末尾处与该CP信息序列对应的第二信息序列;
所述序列补零单元,用于对CP信息序列和第二信息序列分别补零,使补零后的信息序列长度为一个符号的长度;
所述频域变换单元,用于将补零后的CP信息序列和第二信息序列变换到频域,得到CP信息序列频谱和第二信息序列频谱;
所述频谱抑制单元,用于对CP信息序列频谱和第二信息序列频谱进行频谱扩展抑制,得到修正后的CP信息序列频谱和第二信息序列频谱;
所述频偏估计单元,用于根据修正后的CP信息序列频谱和第二信息序列频谱估计频偏。
进一步优选地,所述频谱抑制单元包括载波去除单元,用于去除CP信息序列频谱和第二信息序列频谱的带宽边缘处的子载波,得到修正后的CP信息序列频谱和第二信息序列频谱。
进一步优选地,所述载波去除单元包括载波数确定子单元和载波去除子单元,其中:
所述载波数确定子单元,用于根据当前UE带宽大小和/或当前小区与其它小区的相邻方式确定待去除的带宽边缘处的子载波数;
所述载波去除子单元,用于根据所述子载波数去除CP信息序列频谱和第二信息序列频谱的带宽边缘处的子载波。
优选地,所述装置还包括加权平均单元,用于在对一个子帧内的N个符号进行频偏估计后,对估计得到的N个频偏值进行加权平均,将加权平均后得到的频偏作为最终的频偏,所述N大于1。
本申请实施方式在时域抽样提取CP信息序列和对应的第二信息序列,在对这些信息序列进行补零、变换到频域后,对频谱进行扩展抑制,并基于频谱抑制之后的频谱进行频偏估计,与现有技术相比,由于在频域进行了频谱扩展抑制,使带宽相邻的UE之间的干扰减小,从而有利于更为准确地进行频偏估计。而且,由于一个子帧的符号均带有CP,其间隔短,估计范围在7500Hz以内,经过对多个符号的加权估计后,使频偏估计结果更准确。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
参见图1,该图示出了本申请的基于CP的频偏估计方法的一个实施例流程,该流程包括:
步骤S11:时域抽样提取一个符号的CP信息序列以及位于该符号末尾处与该CP信息序列对应的第二信息序列;
如前所述,在LTE系统中,每个时域符号的起始位置均插入了CP信息,该CP信息来自于该符号尾部的一部分信息,即相当于将一个符号尾部的数据信息复制到符号前端。在不同情况下,插入符号中的CP信息的长度可能不同,比如,CP信息长度的典型值可以是160、144、512Ts。在本步骤中,基站侧对一个符号进行时域抽样,提取出一个符号的CP信息序列和该符号对应的尾部的信息序列(为便于区分,以下称为“第二信息序列”),提取的这两个信息序列的长度一般情况下应当相等,即如果提取的CP信息序列的长度为144个点,那么第二信息序列也应当为144点。
步骤S12:对CP信息序列和第二信息序列分别补零,使补零后的信息序列长度为一个符号的长度;
按照前述步骤提取出两个序列后对其分别进行补零操作,补零操作在频域上将产生“加窗效应”,即出现频谱扩展,频谱扩展的范围与窗函数的宽度大小密切相关。通常情况下,补零操作可通过与窗函数进行运算实现,这种运算可以直接在时域上进行,也可以在频域上实现,如果在时域上,则直接将CP信息序列和第二信息序列分别与窗函数相乘,如果在频域上,则将两个序列的频谱与窗函数频谱进行卷积。对于补零的位数,通常情况下,为便于后续的操作,需要将CP信息序列和第二信息序列通过补零操作将其长度延伸到一个符号的长度。比如,假设一个符号的抽样长度为2048Ts,那么针对144Ts的CP信息序列,应当补充1904个零位。
步骤S13:将补零后的CP信息序列和第二信息序列变换到频域,得到CP信息序列频谱和第二信息序列频谱;为便于后面叙述,这里将变换到频域的CP信息序列和第二信息序列频谱分别记为H0和H1。
步骤S14:对CP信息序列频谱和第二信息序列频谱进行频谱扩展抑制,得到修正后的CP信息序列频谱和第二信息序列频谱;
如前所述,由于经过补零以及由时域到频域的变换等操作,将出现“加窗效应”,频谱发生扩展,使当前UE受到相邻UE的干扰,无法实现准确的频偏估计。假设接收频谱H'(k),则由于相邻UE的干扰,该接收频谱可表示为:
H'(k)=Hr(k)+X(k)
式中:Hr(k)为当前UE的自身频谱,X(k)为干扰当前UE的另一个UE形成的干扰频谱。由该公式可见,通过对干扰频谱的抑制则可以还原本来的频谱。进行频谱扩展抑制后的频谱为修正后的频谱,修正后的CP信息序列频谱和第二信息序列频谱分别记为H'0和H'1。
步骤S15:根据修正后的CP信息序列频谱和第二信息序列频谱估计频偏。
利用修正后的H'0和H'1,根据如下公式计算当前符号频偏大小:
上述实施例在时域抽样提取CP信息序列和对应的第二信息序列,在对这些信息序列进行补零、变换到频域后,对频谱进行扩展抑制,并基于频谱抑制之后的频谱进行频偏估计,与现有技术相比,由于在频域进行了频谱扩展抑制,使带宽相邻的UE之间的干扰减小,从而有利于更为准确地进行频偏估计。此外,本实施例的方法与另一种直接在时域上进行频偏估计的方法相比,效果也更好。这种“时域频偏估计法”利用CP信息和对应符号末端信息在时域上进行相关运算来估计频偏,如果在多UE情况下,由于每个UE频偏不同,各自产生的相位旋转将在时域信号上叠加,使得在时域信号上进行相关运算无法区分开自身UE的频偏和其他UE的干扰,从而导致估计的频偏误差较大,而本申请的上述实施例是将待处理的信息序列变换到频域,补零后的序列经过傅里叶变换能保证频域信息只在当前UE带宽内产生,这就使得从频域中可以明显区分来自不同UE的频偏,从而有利于将干扰频谱给抑制掉,进一步对UE带宽内干扰进行抑制得到修正后的频谱,基于该修正后的频谱进行的频偏估计将更为准确。
在前述实施例的步骤S14中,提到需要对频域的CP序列频谱和第二信息序列频谱进行频谱扩展抑制。通常情况下,频谱扩展抑制的具体方式与补零操作的实现方式相关联。比如,如果在补零操作中采取窗函数进行运算实现补零操作,那么这里可以基于窗函数的特点进行频谱扩展抑制。参见图2,该图示出了一种窗函数的频谱幅值谱,从该图可以看出,窗函数频谱与功率谱一样,属于典型的U型谱,该谱的特点是频谱能量集中在靠窗边的两侧,而中间部分基本上均为零,也就是说,经过窗函数卷积后,频谱能量仍相对集中,扩展产生在频谱两侧附近,干扰频谱对当前UE自身频谱的影响将集中在带宽两侧,如果仅对这部分进行带宽抑制,即可达到对CP信息序列频谱和第二信息序列频谱进行修正的目的。当然,如果补零操作不是采取窗函数的形式,这里的频谱抑制可以根据需要选定需要抑制的频谱的范围,只要抑制掉部分频谱后不影响到正常的频偏估计即可。
通过对图2所示的窗函数功率谱的功率分布范围的统计,可以发现,窗函数功率谱的左右窗边各较少数量内的点包含了绝大部分的能量,比如,一种情况下,窗两边各11点可以包含总能量的89.71%,那么可以考虑对这11个点进行频谱抑制。假设当前UE的CP信息序列频谱和第二信息序列频谱边缘带宽受到来自相邻UE的干扰分别为:
X0((end-10)~end)
X1((end-10)~end)
式中:end表示当前UE的带宽边缘,由于X是来自相邻UE的信号展宽,这里令δ为相邻UE的频偏角,则X0和X10之间的关系可表示为:X0e2jπδ=X1。由于干扰频谱与自身频谱来自不同的UE,两者不相关,在统计条件下,存在如下的等式:
E[(H'0)HH'1]=E[(Hr 0)HHr 1]+E[(X0((end-10)~end))HX1((end-10)~end)]
同理,若当前UE的带宽两侧都具有相邻UE,则:
E[(H'0)HH'1]=E[(Hr 0)HHr 1]+E[(X0(1~11))HX1(1~11))]+
E[(X0((end-10)~end))HX1((end-10)~end)]
为了说明通过抑制频谱扩展可以实现对干扰频谱的消除,下面结合图3所示的数学模型来进行分析:图3所示的数学模型与前述两个公式具有对应关系,上述式子中右边第一项的模对应当前UE的频谱功率大小,即半径A,弧角为当前UE的频偏角2πε,第二(前一公式)或第三项(后一公式)的模对应当前UE的相邻UE的干扰频谱功率大小,即半径B,弧角为该相邻UE的频偏角2πδ,图中表示最大频偏干扰理论值的大小,其最大为(B/A)。如果对当前UE的带宽边缘(即当前UE受到相邻UE干扰的部分)进行抑制,且使则由图可见,δ的影响将忽略不计,从而完全消除了相邻UE的频谱干扰。
基于前述对频谱扩展进行抑制原理的介绍,本领域技术人员可推知其具体实现方式可以有多种,下面仅示例性的给出两种,本领域技术人员可在此基础之上进行各种变形或改进,从而得到更多的实现方式,这些实现方式均在本申请的保护范围之内。
示例性方式之一:子载波去除顺推法
该方法通过去除当前UE边缘部分的子载波来实现对频谱扩展的抑制,对于去除子载波的数量可通过如下的方法推理得到。这里假设基站接收到的两个UE信号平均到每个载波上的功率相等为1个单位,令P(k)表示窗函数功率谱,其中窗函数定义为:
如果将去除当前UE带宽边缘的d个子载波,则可得到去除带宽边缘的d个子载波后的真实信号(当前UE的频谱)能量A和干扰信号(当前UE的相邻UE的频谱)能量B分别为:
式中:N为当前UE的带宽大小,单位为子载波。基于则可以得到如下的公式:
上述不等式的左边为常数,右边只含有变量d,通过解算该不等式可得到变量d的范围,即d≤dmax,从而确定了将去除的子载波数。
示例性方式之二:子载波去除试验法
通过上述的示例性方法能够得到将要去除的子载波数量,但是,前述示例性方法中的不等式的解算过程相当复杂,将耗费较多的计算资源。为此,可以考虑通过试验的方法得到需要去除的子载波数,即将去除的子载波数分别依次代入公式中,分别计算出最大频偏干扰理论值,如果该频偏干扰理论值小于预设的频偏干扰理论值,则这样的子载波数去除将能达到要求。比如,假设预设的最大频偏干扰理论值如果限制在21Hz以内,当前UE的带宽为2~9PRB(一个PRB=12子载波),则将d=7、6、5、5、4、4、3、3(子载波)分别代入上述示例性方式的公式之中,计算得到的最大频偏干扰理论值分别为:
17.7002、17.1683、19.4499、15.0257、18.7424、15.8430、20.0211、17.6732
这些计算出来的最大频偏干扰理论值均小于预设的最大频偏干扰理论值,因此,去除这些数量的子载波数将能够满足预定要求。下述表格示例性地列明某种情况下的当前UE带宽与去除的子载波数之间的对应关系:
带宽(RB) |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10以上 |
去掉的子载波数 |
7 |
6 |
5 |
5 |
4 |
4 |
3 |
3 |
0 |
上述两种示例性方式均提及到去除当前UE带宽边缘部分的子载波,尽管可以通过不同的方式确定去除的子载波的数量,但是,在实际应用过程中,去除子载波的数量还需要考虑一些其他因素。比如,考虑去除子载波后剩余子载波对于频偏估计的影响程度,通常情况下,不是去除越多的子载波越好,如果去除过多子载波可能使得参与估计频偏的子载波的数量不足影响到频偏的估计精度。还比如,考虑某个小区内的当前UE受到其他相邻UE的干扰方式(或相邻方式),不同的干扰方式可能使去除子载波的数量不同。子载波的干扰方式(相邻方式)通常可以包括如下三种,上边缘相邻、下边缘相邻、两边缘相邻,由于基站拥有所有UE带宽的分配信息,这些干扰方式可由基站作出判断。一般而言,如果仅存在上边相邻或下边缘相邻,则去除其一侧的子载波即可,如果两边缘相邻,则需要去除两侧的子载波。
为了更进一步说明本申请的技术方案,下面给出一个更为详细的实施例。参见图4,该图示出了该实施例的流程,该流程包括:
步骤S40:初始化i=1,i为当前处理的子帧中的符号序数;
步骤S41:基站侧时域抽样提取第i个符号的144点CP信息和对应符号末尾处的144点信息:
步骤S42:在两个抽样序列前面补零至长度为一个符号的抽样长度NFFT,形成两个新序列h'0和h'1;
步骤S43:将补零后两新序列经傅里叶变换变换到频域,分别记为H′0和H′1;
步骤S44:判断当前小区内是否有其他干扰UE(下称UE2)与当前UE(下称UE1)带宽相邻,并判断相邻方式;
步骤S45:根据当前UE带宽大小计算需要去除的带宽边缘的子载波数;
步骤S46:根据干扰UE相邻方式和需要去除的子载波数,去除当前UE带宽边缘一定数量的子载波,得到修正后的H′0和H′1;
步骤S47:利用修正后的H′0和H′1,根据如下公式计算当前符号的频偏大小:
步骤S48:判断i是否达到子帧内的符号个数最大值,如果否,则i=i+1,返回步骤S41;如果是,则执行步骤S49;
步骤S49:将子帧内的全部符号的频偏值频偏值进行加权平均,将加权平均的结果作为当前UE的最终频偏值。
上述内容详细叙述了本申请的基于CP的频偏估计方法的实施例,与之对应地,本申请还提供了一种基于CP的频偏估计方法装置实施例。参见图5,该图示出了该基于CP的频偏估计装置实施例的组成结构。该装置包括:序列提取单元U51、序列补零单元U52、频域变换单元U53、频谱抑制单元U54和频偏估计单元U55,其中:
序列提取单元U51,用于在时域抽样提取一个符号的CP信息序列以及位于该符号末尾处与该CP信息序列对应的第二信息序列;
序列补零单元U52,用于对CP信息序列和第二信息序列分别补零,使补零后的信息序列长度为一个符号的长度;
频域变换单元U53,用于将补零后的CP信息序列和第二信息序列变换到频域,得到CP信息序列频谱和第二信息序列频谱;
频谱抑制单元U54,用于对CP信息序列频谱和第二信息序列频谱进行频谱扩展抑制,得到修正后的CP信息序列频谱和第二信息序列频谱;
频偏估计单元U55,用于根据修正后的CP信息序列频谱和第二信息序列频谱估计频偏。
上述装置实施例的工作过程是:序列提取单元U51在时域抽样提取一个符号的CP信息序列以及位于该符号末尾处与该CP信息序列对应的第二信息序列,然后由序列补零单元U52对CP信息序列和第二信息序列分别补零,使补零后的信息序列长度为一个符号的长度;再由频域变换单元U53将补零后的CP信息序列和第二信息序列变换到频域,得到CP信息序列频谱和第二信息序列频谱,频谱抑制单元U54对CP信息序列频谱和第二信息序列频谱进行频谱扩展抑制,得到修正后的CP信息序列频谱和第二信息序列频谱;最后由频偏估计单元U55根据修正后的CP信息序列频谱和第二信息序列频谱估计频偏。本装置实施例可以取得与前述方法实施例相同或类似的技术效果,为避免重复,这里不再赘言。
在前述装置实施例的基础之上,本领域技术人员可以基于各种实际需要对上述实施例进行改进或变形。比如,频谱抑制单元U54可以进一步包括载波去除单元U541,用于去除CP信息序列频谱和第二信息序列频谱的带宽边缘处的子载波,得到修正后的CP信息序列频谱和第二信息序列频谱,而载波去除单元U541可以进一步包括载波数确定子单元U5411和载波去除子单元U5412,其中:
载波数确定子单元U5411,用于根据当前UE带宽大小和/或当前小区与其它小区的相邻方式确定待去除的带宽边缘处的子载波数;
载波去除子单元U5412,用于根据所述子载波数去除CP信息序列频谱和第二信息序列频谱的带宽边缘处的子载波。
此外,上述装置实施例还可以包括加权平均单元U56,用于在对一个子帧内的N个符号进行频偏估计后,对估计得到的N个频偏值进行加权平均,将加权平均后得到的频偏作为最终的频偏,LTE中所述N取14或12。这样通过多次加权平均后,可得到更为准确的频偏估计值。
需要说明的是:为了叙述的简便,本说明书的上述实施例以及实施例的各种变形实现方式重点说明的都是与其他实施例或变形方式的不同之处,各个情形之间相同相似的部分互相参见即可。尤其,对于装置实施例的几个改进方式而言,由于其基本相似于方法实施例,所以描述得比较简单,相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例的各单元可以是或者也可以不是物理上分开的,既可以位于一个地方,或者也可以分布到多个网络环境下。在实际应用过程中,可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的,本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
以上所述仅是本申请的具体实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本申请原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本申请的保护范围。