CN102123127B - 一种时域信道估计滤波器的系数计算方法和装置 - Google Patents
一种时域信道估计滤波器的系数计算方法和装置 Download PDFInfo
- Publication number
- CN102123127B CN102123127B CN2010106242966A CN201010624296A CN102123127B CN 102123127 B CN102123127 B CN 102123127B CN 2010106242966 A CN2010106242966 A CN 2010106242966A CN 201010624296 A CN201010624296 A CN 201010624296A CN 102123127 B CN102123127 B CN 102123127B
- Authority
- CN
- China
- Prior art keywords
- symbol
- information
- position information
- channel estimation
- vector
- 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.)
- Expired - Fee Related
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 85
- 239000013598 vector Substances 0.000 claims abstract description 73
- 239000011159 matrix material Substances 0.000 claims abstract description 70
- 238000000034 method Methods 0.000 claims abstract description 20
- 238000005314 correlation function Methods 0.000 claims description 12
- 238000006386 neutralization reaction Methods 0.000 claims description 4
- 101100460704 Aspergillus sp. (strain MF297-2) notI gene Proteins 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 238000001914 filtration Methods 0.000 description 13
- 238000010586 diagram Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 6
- 238000013461 design Methods 0.000 description 5
- 238000007781 pre-processing Methods 0.000 description 4
- 125000004122 cyclic group Chemical group 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2655—Synchronisation arrangements
- H04L27/2662—Symbol synchronisation
- H04L27/2663—Coarse synchronisation, e.g. by correlation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/0212—Channel estimation of impulse response
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/0224—Channel estimation using sounding signals
- H04L25/0228—Channel estimation using sounding signals with direct estimation from sounding signals
- H04L25/023—Channel estimation using sounding signals with direct estimation from sounding signals with extension to other symbols
- H04L25/0232—Channel estimation using sounding signals with direct estimation from sounding signals with extension to other symbols by interpolation between sounding signals
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
- H04L25/0242—Channel estimation channel estimation algorithms using matrix methods
- H04L25/0244—Channel estimation channel estimation algorithms using matrix methods with inversion
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
- H04L25/0256—Channel estimation using minimum mean square error criteria
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Mobile Radio Communication Systems (AREA)
- Synchronisation In Digital Transmission Systems (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供一种时域信道估计滤波器的系数计算方法和装置,所述方法包括:获取X1,X2…XP所对应的携带RS信息的符号的位置信息和权重信息;获取Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息;根据所述位置信息和权重信息,计算向量I的自相关矩阵,其中,I=[X1,X2…XP,Y1,Y2,…YQ];根据所述位置信息和权重信息,计算向量O与向量I之间的互相关矩阵,其中,O=[Y0],Y0为当前进行信道估计的符号的时域的理想值;计算所述互相关矩阵与所述自相关矩阵的逆矩阵之积,得到当前进行信道估计的符号对应的系数向量。本发明能够计算出较优的滤波器系数。
Description
技术领域
本发明属于无线通信技术领域,特别涉及一种应用于长期演进(LTE)***中的时域信道估计滤波器的系数计算方法和装置。
背景技术
在LTE无线通信***中,正交频分复用(OFDM)是其核心的物理层技术。它把数据流调制在多个正交的子载波上,正交子载波之间频谱可重叠,这样就大大提高了频谱利用率。
参照图1,在正交频分复用***中,传输的信息可以通过资源栅格来描述,资源栅格将整个时域和频域的资源以资源元素(RE,Resource Element)的形式表示出来,即图1中的一个最小方格,它对应于频域上的一个子载波和时域上的一个符号(symbol)的时间长度,所有需要传输的信息都是通过资源元素来承载。多个资源元素构成一个资源块,具体来说,在正常循环前缀(CP)的情况下,12(子载波数)×7(符号数)构成一个资源块;扩展循环前缀的情况下,12(子载波数)×6(符号数)构成一个资源块。
为了正确解调数据,需要知道每个RE位置上的信道情况。在LTE***中,每个RE位置上的信道估计是通过如下方式得到的:在每个资源块的特定的资源元素位置***已知的参考信号(RS,Reference Signal)以计算该RE位置上的信道估计,然后利用插值的方法来得到其它所有资源元素位置上的信道估计。
信道估计算法通常分两步来进行:先在频域方向上进行频域信道估计,即先对图1中具有RS的一列符号在频域上对RS信息的初始信道估计进行滤波插值,得到这一列每个符号的频域信道估计,然后利用该结果在时域方向上进行滤波插值来完成时域信道估计,即最终的信道估计。为使后文描述简单,用表示携带RS信息的符号的频域信道估计,包括通过频域插值得到的信道估计。
LTE***分为TDD和FDD两种模式。对于FDD模式来说,UE可以接收到连续的下行子帧;对于TDD来说,其接收到的下行子帧不是连续的,即RS的分布不是连续的,如图2所示,图中,用D表示下行子帧,用U表示上行子帧。
实际上,即使是FDD模式可以接收到连续的下行子帧,但是在正常CP模式下,其在一个子帧中RS分布是不均匀的,携带RS信息的符号间隔有3和4两种,如图3所示。
频域滤波插值和时域滤波插值的关键是选择合适的滤波器结构和滤波器系数。而LTE***中RS分布的不连续和不均匀的特性,给时域信道估计滤波器的选择和系数设计带来了挑战。
发明内容
本发明的目的是提供一种时域信道估计滤波器的系数计算方法和装置,以克服LTE中RS分布的不连续性和不均匀性带来的问题。
为实现上述目的,本发明提供技术方案如下:
一种时域信道估计滤波器的系数计算方法,应用于LTE***中,包括:
获取X1,X2…XP所对应的携带RS信息的符号的位置信息和权重信息,X1,X2…XP为与当前输入符号最接近的P个频域信道估计值,P为大于0的整数;
获取Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息,Y1,Y2,…YQ为与当前进行信道估计的符号最接近的Q个时域信道估计值,Q为大于0的整数;
根据所述位置信息和权重信息,计算向量I的自相关矩阵,其中,I=[X1,X2…XP,Y1,Y2,…YQ];
根据所述位置信息和权重信息,计算向量O与向量I之间的互相关矩阵,其中,O=[Y0],Y0为当前进行信道估计的符号的时域信道的理想值;
计算所述互相关矩阵与所述自相关矩阵的逆矩阵之积,得到当前进行信道估计的符号对应的系数向量。
一种时域信道估计滤波器的系数计算装置,应用于LTE***中,包括:
第一获取模块,用于获取X1,X2…XP所对应的携带RS信息的符号的位置信息和权重信息,X1,X2…XP为与当前输入符号最接近的P个频域信道估计值,P为大于0的整数;
第二获取模块,用于获取Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息,Y1,Y2,…YQ为与当前进行信道估计的符号最接近的Q个时域信道估计值,Q为大于0的整数;
自相关计算模块,用于根据所述位置信息和权重信息,计算向量I的自相关矩阵,其中,I=[X1,X2…XP,Y1,Y2,…YQ];
互相关计算模块,用于根据所述位置信息和权重信息,计算向量O与向量I之间的互相关矩阵,其中,O=[Y0],Y0为当前进行信道估计的符号的时域信道的理想值;
滤波系数计算模块,用于计算所述互相关矩阵与所述自相关矩阵的逆矩阵之积,得到当前进行信道估计的符号对应的系数向量。
附图说明
图1为正常循环前缀且单天线发送情况下,参考信号在资源块中的分布情况示意图;
图2为TDD模式,正常CP模式,上下行转换配置为1#的子帧分布示意图;
图3为FDD模式,正常CP模式,携带RS信息的符号在子帧中的分布示意图;
图4为本发明实施例中采用的时域信道估计滤波器的结构示意图;
图5为本发明实施例的时域信道估计滤波器的系数计算方法流程图;
图6为本发明实施例的时域信道估计滤波器的系数计算装置结构图;
图7为本发明实施例中LTE TDD帧结构的时域滤波场景划分示意图;
图8为利用图6所示的装置进行滤波器系数计算的一个详细流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图及具体实施例对本发明进行详细描述。
图4为本发明实施例中采用的时域信道估计滤波器的结构示意图,本发明是针对图4所示的时域信道估计滤波器结构,根据可用的资源、信道信息和***信息来设计时域信道估计滤波器的系数。
参照图4,该插值滤波器使用P个寄存器保存与当前输入符号(序号为n)最接近的P个携带RS信息的符号的频域信道估计值X1~XP(包括通过频域插值得到的信道估计值),并使用Q个寄存器保存与当前进行信道估计的符号(序号为n-L)最接近的Q个时域信道估计的历史值Y(n-L-1)~Y(n-L-Q)(即利用该滤波器进行插值滤波得到的Q个时域信道估计值,后文也可以用Y1,Y2,…YQ表示,Y1即为Y(n-L-1),Y2即为Y(n-L-2)…),根据保存的值对当前进行信道估计的符号进行时域信道估计的插值滤波。这一过程用公式表示为:
当Q≠0时,这是一个无限冲激响应(IIR)滤波器结构;当Q=0时,这是一个有限冲激响应(FIR)滤波器结构。这里的L表示滤波器的输入和输出之间有多少个符号的时延,通常,时延越多,可用于时域信道估计的频域信道估计越多,越有利于时域信道估计,但在另一方面它又增加了***时延,这是不利的一面,实际应用中可以根据具体需求来对这两方面进行折中。
与一般滤波器结构不同的是,这里的滤波器结构有一个更新模块。这是由于这里滤波器的输入是携带RS信息的符号的频域信道估计然而不是每个符号都携带RS信息的。参照图1,以正常CP模式并且单天线端口来说,一个子帧中的14个符号中只有4个符号是携带RS信息的,也即只有符号序号为1、5、8和12的4个输入是有效的(序号计数从1开始)。更新模块可有多种实现方式,这里介绍模式1和模式2两种方式。
对于模式1,更新模块就是把与当前输入符号最接近的携带RS信号的符号的频域信道估计值更新到P个寄存器当中。用LRS表示包含RS信息的符号的序号,其操作过程可用伪码表示如下:
其中,X(n)表示当前输入符号的频域信道估计值。
可见,按照模式1的更新为:每输入一个新的符号时,判断该符号是否为携带RS信息的符号,若是,从i=P开始,将Xi更新为Xi-1,直到i=1时,将X1更新为X(n);否则,不进行更新,即,X1~XP中的每一个都保持不变。
这样,与当前输入符号最接近的P个携带RS信息的符号(可以包括当前输入符号)的频域信道估计就被存放在X1~XP中了。
对于模式2,其操作过程可用伪码表示如下:
其中,X(n)表示当前输入符号的频域信道估计值。
可见,按照模式2的更新为:每输入一个新的符号时,从i=P开始,将Xi更新为Xi-1,直到i=1时,判断该符号是否为携带RS信息的符号,若是,将将X1更新为X(n);否则,不对X1进行更新。
这样,与当前输入符号最接近的P个携带RS信息的符号(可以包括当前输入符号)的频域信道估计就被存放在X1~XP中了,这里X1~XP的频域信道估计信息可能是有重复的(X1与X2可能存储的是同一个值)。
针对上述时域信道估计滤波器结构,本发明实施例提供如下的滤波器系数计算方法。
参照图5,本发明实施例的时域信道估计滤波器的系数计算方法,应用于LTE***中,主要包括如下步骤:
步骤501:获取X1,X2…XP所对应的携带RS信息的符号的位置信息和权重信息,X1,X2…XP为与当前输入符号最接近的P个频域信道估计值,P为大于0的整数;
步骤502:获取Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息,Y1,Y2,…YQ为与当前进行信道估计的符号最接近的Q个时域信道估计值,Q为大于0的整数;
步骤503:根据所述位置信息和权重信息,计算向量I的自相关矩阵,其中,I=[X1,X2…XP,Y1,Y2,…YQ];
步骤504:根据所述位置信息和权重信息,计算向量O与向量I之间的互相关矩阵,其中,O=[Y0],Y0为当前进行信道估计的符号的时域信道的理想值;
步骤505:计算所述互相关矩阵与所述自相关矩阵的逆矩阵之积,得到当前进行信道估计的符号对应的系数向量。
本发明实施例对步骤503和步骤504的执行顺序不做限制,即,可以先执行步骤504,再执行步骤503,另外,两者还可以并行执行。
上述方法各步骤的具体实现可参见后文对本发明实施例的时域信道估计滤波器的系数计算装置的描述。
图6为本发明实施例的时域信道估计滤波器的系数计算装置结构图,参照图6,所述系数计算装置主要包括:第一获取模块10、第二获取模块20、自相关计算模块30、互相关计算模块40、滤波系数计算模块50和滤波系数选择模块60。以下对上述各模块的工作原理进行详细说明。
第一获取模块10
第一获取模块10是根据时域信道估计滤波器所采用的更新模式,来记录当前X1~XP所对应的携带RS信息的符号的位置信息和权重信息。如前所述,X1,X2…XP为与当前输入符号(序号为n)最接近的P个频域信道估计值(这里,所谓最接近,可以包括当前输入符号自身),P为大于0的整数。将X1~XP对应的位置信息和权重信息分别表示为和位置信息和权重信息需要被初始化为无效值,即,
可见,按照模式1的更新为:每输入一个新的符号时,判断该符号是否为携带RS信息的符号,若是,从i=P开始,将更新为将更新为直到i=1时,将更新为n,将更新为1,n为该符号的序号;否则,不对所述位置信息和权重信息进行更新,即,和中的每一个都保持不变。
可见,按照模式2的更新为:每输入一个新的符号时,从i=P开始,将更新为将更新为直到i=1时,判断该符号是否为携带RS信息的符号,若是,将更新为n,将更新为1,n为该符号的序号;否则,不对和进行更新。这里X1~XP对应的位置信息和权重信息有可能是重复的,例如,X1与X2可能对应相同的位置信息和权重信息。
自相关计算模块30
自相关计算模块30是计算所有输入变量之间的相关性,并用矩阵的形式表示出来。输入变量可用向量I=[X1,X2…XP,Y1,Y2,…YQ]表示,如前所述,Y1,Y2,…YQ为与当前进行信道估计的符号最接近的Q个时域信道估计值,即,利用所述时域信道估计滤波器进行插值滤波得到的Q个历史值,Q为大于0的整数。
具体地,自相关计算模块30是根据X1,X2…XP所对应的携带RS信息的符号的位置信息和权重信息,以及,Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息,按照如下公式来计算向量I的自相关矩阵RII:
RII中各元素的计算过程可用如下公式表示:
其中,E[·]为求数学期望运算,为Xm对应的位置信息,为Xm对应的权重信息,m=1,…P,为Yk对应的位置信息,为Yk对应的权重信息,k=1,…Q,为的第j个元素,为的第j个元素,length(·)为求向量的元素个数,R(τi,j)为根据信道信息和***信息确定的相关函数,τi,j为第i个符号与第j个符号之间的时间距离。
相关函数R(τi,j)的计算除了需要和以外,还需要信道信息和***信息。信道信息包括多普勒频移、信噪比、所使用的信道衰落模型等,***信息包括CP模式和符号长度等。例如,滤波器的系数设计所使用的信道信息和***信息如下表所示:
那么,所述相关函数可表示为,
R(τi,j)=J0(2πfmaxτi,j)
J0(·)表示第一类0阶贝塞尔函数,τi,j=Tsym|i-j|,Tsym为一个符号的时间长度。
针对其他信道信息和***信息,本领域技术人员可以相应地选用其他类型的相关函数。
需要说明的是,在自相关矩阵的计算过程中,所使用的Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息为:当前进行信道估计的符号对应的Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息,所述位置信息和权重信息为对上一个符号进行系数后更新得到,并且,初始化时,所述位置信息和权重信息为空集。
在得到当前进行信道估计的符号对应的系数向量之后,本发明实施例还对Y1,Y2,…YQ所涉及的携带RS信息的符号位置信息和权重信息进行更新,以供下一个符号的系数计算使用,具体请参见后文对第二获取模块20的描述。
也就是说,本发明实施例对滤波系数的计算是一个迭代过程,本次使用的Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息,是上一次系数计算后更新的结果。
另外,本次使用的X1,X2…XP所对应的携带RS信息的符号的位置信息和权重信息,是在本次系数计算之前,根据输入的新的符号(即当前符号n)进行更新得到。
互相关计算模块40
互相关计算模块40是计算所有输入变量和输出变量Y0(当前进行信道估计的符号的时域信道的理想值)之间的相关性,并用矩阵的形式表示出来。输入变量可用向量I=[X1,X2…XP,Y1,Y2,…YQ]表示,如前所述,Y1,Y2,…YQ为与当前进行信道估计的符号最接近的Q个时域信道估计值,即,利用所述时域信道估计滤波器进行插值滤波得到的Q个历史值,Q为大于0的整数。
具体地,互相关计算模块40是根据X1,X2…XP所对应的携带RS信息的符号的位置信息和权重信息,以及,Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息,按照如下公式来计算向量O与向量I之间的互相关矩阵ROI,其中,O=[Y0]:
其中,E[·]为求数学期望运算,为Xp对应的位置信息,为Xp对应的权重信息,p=1,…P,为Yq对应的位置信息,为Yq对应的权重信息,q=1,…Q,为的第i个元素,为的第i个元素,length(·)为求向量的元素个数,R(τi,j)为根据信道信息和***信息确定的相关函数,τi,j为第i个符号与第j个符号之间的时间距离,n-L为当前进行信道估计的符号的序号,n为当前输入符号的序号,L表示滤波器的输入与输出之间有多少个符号的时延。
这里,相关函数R(τi,j)的选择请参见前文。
需要说明的是,在互相关矩阵的计算过程中,所使用的Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息为:当前进行信道估计的符号对应的Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息,所述位置信息和权重信息为对上一个符号进行系数计算后更新得到,并且,初始化时,所述位置信息和权重信息为空集。
另外,本次使用的X1,X2…XP所对应的携带RS信息的符号的位置信息和权重信息,是在本次系数计算之前,根据输入的新的符号(即当前符号n)进行更新得到。
还需要说明的是,这里的Y0代表当前进行信道估计的符号的时域信道的理想值,即,Y(n-L),其实,在互相关矩阵的计算中,是不需要知道Y(n-L)的实际值的,使用的只是其对应的位置信息,即,当前进行信道估计的符号的序号:n-L。
滤波系数计算模块50
滤波系数计算模块50是根据预定准则,利用自相关矩阵RII和互相关矩阵ROI来为当前进行信道估计的符号(第n-L个符号)的时域信道估计计算出一个较优的系数,原理为:计算所述互相关矩阵与所述自相关矩阵的逆矩阵之积,得到当前进行信道估计的符号对应的系数向量。
作为实施例,这里给出的滤波系数计算方案包括矩阵预处理、系数计算和系数补偿。
由于矩阵RII在某些情况下不满秩,需要删除特定的行和列,即通过缩维的方式来确保矩阵满秩。
矩阵预处理根据X1,X2…XP和Y1,Y2,…YQ中包含的RS信息来判断RII是否满秩,如果不满秩,需要删除特定的行和列,即通过缩维的方式来确保矩阵满秩。用R′II表示更新后的矩阵,用和表示被删掉的行和列所对应的变量,即被删除的行和列为和所对应的行和列,那么和在ROI中所对应的元素也应被删除,得到新的矩阵R′OI。
例如,X1,X2…XP和Y1,Y2,…YQ共包含M个RS的信息:
当M<P+Q时,这时需要删除RII的P+Q-M行和列以保证矩阵满秩,ROI中所对应的元素也应被删除,优先删除位置和权重均为0的Xi所对应的行和列,以及位置和权重均为空的Yj,被删除变量所对应的系数被置为0;
当M≥P+Q时,如果X1,X2…XP中存在多个变量对应同一个RS信息符号的情况时(例如,第一获取模块10使用模式2更新),同一个RS信息符号所对应的多个变量中,只保留一个变量,其他变量所对应的行和列在RII中删除,ROI中所对应的元素也应被删除,被删除变量所对应的系数被置为0。同样地,如果Y1,Y2,…YQ中存在多个变量包含相同RS信息的情况时(位置向量和权重向量都完全相同),只保留一个变量,其他变量所对应的行和列在RII中删除,ROI中所对应的元素也应被删除,被删除变量所对应的系数被置为0。
系数计算是根据最小均方误差(MMSE)准则,通过R′II和R′OI计算时域信道估计滤波器的系数。
矩阵预处理过程中没有变量被删除时,此时的系数计算可表示为:
系数补偿是对所计算的系数乘以一个因子β,使其归一化。例如,如果使用幅度归一化,β可按如下公式计算:
第二获取模块20
第二获取模块20用于获取Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息,Y1,Y2,…YQ为与当前进行信道估计的符号最接近的Q个时域信道估计值,Q为大于0的整数。
所述获取Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息为:获取当前进行信道估计的符号对应的Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息,所述位置信息和权重信息为对上一个符号进行系数计算后更新得到,并且,初始化时,所述位置信息和权重信息为空集。在得到当前进行信道估计的符号对应的系数向量之后,第二获取模块20还用于:对Y1,Y2,…YQ所涉及的携带RS信息的符号位置信息和权重信息进行更新,以供下一个符号的系数计算使用。
第二获取模块20是根据当前进行信道估计的符号(第n-L个符号)的滤波系数[b1,…bP,a1,…aQ],X1,X2…XP和Y1,Y2,…YQ中所包含的的位置信息、权重信息,更新当前Y(n-L-1)~Y(n-L-Q)中所包含的的位置信息的位置信息和权重信息。由于Y(n-L-q),1≤q≤Q可以表示成所有序号小于或等于n-q的RS的频域信道估计的一个合并,即,
hf(i)表示第i个符号的频域方向上的信道估计。第二获取模块20就是记录位置信息{i|i≤n-q,i∈LRS}以及权重信息{αi|i≤n-q,i∈LRS}。
第二获取模块20对Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息进行更新,可由如下5步完成:
步骤1:定义临时变量D(temp)和W(temp)
步骤2:
步骤3:
For q=2:Q
End
步骤4:
步骤5:
For q=Q:-1:2
End
上述过程用文字描述如下:
令
从q=2开始,对W(temp)进行如下的更新,直到q=Q:
从p=1开始,对W(temp)和D(temp)进行如下的更新,直到p=P:
其中,为Xp对应的位置信息,p=1,…P,为Yq对应的位置信息,为Yq对应的权重信息,q=1,…Q,[b1,…bP,a1,…aQ]为当前进行信道估计的符号对应的系数向量,length(·)为求向量的元素个数。
滤波系数选择模块60(可选模块)
LTE中RS的分布存在不连续和不均匀的特性,如果将同一组系数应用于所有的符号将不能达到较好的插值效果,因此,如果对属于不同场景的符号使用不同的系数则能达到较好的插值效果。采用滤波系数选择模块60后,本发明实施例中的时域信道估计滤波器就不等同于常规意义上的只有同一组系数的时域信道估计滤波器。
对于一个给定的帧结构、天线端口和CP模式,时域信道估计滤波器会对应特定的一系列滤波场景。这里所谓的滤波场景与时间距离向量一一对应,时间距离向量是指:用于当前符号信道估计的X1,X2…XP和Y1,Y2,…YQ中有效元素(非空非零元素)所对应的符号与当前进行信道估计的符号之间的距离向量T=[T1,T2…TP′,TP′+1,…TP′+Q′],0≤P′≤P,0≤Q′≤Q,对于每一个可能的向量T称为一个场景。滤波系数选择模块60就是根据需要选择特定场景的滤波器系数。
这里的场景主要分为两类,起始场景和稳定场景。起始场景是指在滤波器的初始阶段,X1,X2…XP和Y1,Y2,…YQ中有效元素从0变化到P或Q的一系列场景。稳定场景是指当n大于一定值的时候,系数计算过程进入到一个稳定状态,这时有一部分滤波场景会周期性地往复出现,与此相对应的滤波器的系数在时域上呈现出一定的周期往复特性,将这一部分场景成为稳定场景。初始场景的滤波系数需要在滤波器最开始的时候在滤波系数选择模块60中被记录;稳定场景的滤波系数需要在滤波器达到稳定状态的时候在滤波系数选择模块60中被记录。
以TDD帧结构,正常CP模式,上下行子帧配置为1#,特殊子帧配置为0#(DwPTS有3个符号)为例,那么可用图7所示的6组场景覆盖所有的情况:场景1:起始时隙;场景2:特殊时隙的前一个时隙;场景3:特殊子帧中的DwPTS;场景4:紧邻上行子帧的时隙;场景5:连续时隙;场景6:结束时隙。图中未标注数字的为上行时隙、特殊子帧中的UpPTS或者特殊子帧中的保护间隔。
对于第三组场景,需要存3(P+Q)个系数,3为DwPTS中的符号数;对于其他场景组需要存储7(P+Q)个系数,7为一个下行时隙的符号数。
实际上,滤波器从初始场景到稳定场景的过程中,还存在一个过渡阶段,这一阶段的符号所对应的时间距离向量均包含在稳定场景所对应的时间距离向量之中,但是对于给定的时间距离向量,它所对应的系数却与稳定场景对应的滤波系数存在一定差别,这一部分系数可以根据需要选择存储或不存储。
图8为利用图6所示的装置进行滤波器系数计算的一个详细流程图,包括如下步骤:
步骤801:初始化处理;
这里N值的设定应当满足两个条件:首先应当保证所有的滤波场景都能遍历到,该项决定于所应用的帧结构;其次是系数计算能够进入到稳定状态,该项决定于所设定的多普勒频移。
步骤803:判断n-L是否为期望的时域信道估计位置,若是,进入步骤804,否则,进入步骤810;
当符号n-L不是一个期望的时域信道估计位置时,例如当前位置是TDD模式的上行子帧,自相关计算模块30、互相关计算模块40、滤波系数计算模块50和第二获取模块20将被跳过而不进行任何计算。
步骤804:由自相关计算模块30计算向量I的自相关矩阵RII;
步骤805:由互相关计算模块40计算向量O与向量I之间的互相关矩阵ROI;
本发明实施例对步骤804和步骤805的执行顺序不做限制,即,可以先执行步骤805,再执行步骤804,另外,两者还可以并行执行。
步骤806:由滤波系数计算模块50计算滤波系数b1~bP、a1~aQ
步骤807:判断是否启用滤波系数选择模块60,若是,进入步骤808,否则,进入步骤810;
步骤808:由滤波系数选择模块60判断当前符号所属的场景,并将滤波系数与该场景进行对应;
步骤810:令n=n+1;
步骤811:判断n是否大于N,若是,结束,否则,返回步骤802。
综上所述,本发明实施例利用迭代运算的方式,记录当前可用资源X1,X2…XP和Y1,Y2,…YQ中所包含的的位置信息及权重信息,根据位置信息、权重信息、信道信息、***信息计算可用资源之间的相关性以及可用资源与当前符号的相关性,进而计算出较优的滤波器系数。如此,能够有效克服LTE中RS分布的不连续性和不均匀性带来的问题,适用于LTE***中任何帧结构、任何天线端口的时域信道估计滤波器的系数设计。
最后应当说明的是,以上实施例仅用以说明本发明的技术方案而非限制,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神范围,其均应涵盖在本发明的权利要求范围当中。
Claims (20)
1.一种时域信道估计滤波器的系数计算方法,应用于LTE***中,其特征在于,包括:
获取X1,X2…XP所对应的携带参考信号RS信息的符号的位置信息和权重信息,X1,X2…XP为与当前输入符号最接近的P个频域信道估计值,P为大于0的整数;
获取Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息,Y1,Y2,…YQ为与当前进行信道估计的符号最接近的Q个时域信道估计值,Q为大于0的整数;
根据所述位置信息和权 重信息,计算向量I的自相关矩阵,其中,I=[X1,X2…XP,Y1,Y2,…YQ];
根据所述位置信息和权重信息,计算向量O与向量I之间的互相关矩阵,其中,O=[Y0],Y0为当前进行信道估计的符号的时域信道的理想值;
计算所述互相关矩阵与所述自相关矩阵的逆矩阵之积,得到当前进行信道估计的符号对应的系数向量。
6.如权利要求4或5所述的系数计算方法,其特征在于,所述相关函数为:
R(τi,j)=J0(2πfmaxτi,j)
其中,J0(·)表示第一类0阶贝塞尔函数,fmax为最大多普勒频移,τi,j=Tsym|i-j|,Tsym为一个符号的时间长度。
7.如权利要求1所述的系数计算方法,其特征在于,当所述自相关矩阵不满秩时,所述系数计算方法还包括:
通过缩维的方式将所述自相关矩阵更新为满秩矩阵,并对所述互相关矩阵进行相应更新;
此时,所述计算所述互相关矩阵与所述自相关矩阵的逆矩阵之积为,计算更新后的互相关矩阵与更新后的自相关矩阵的逆矩阵之积。
8.如权利要求7所述的系数计算方法,其特征在于,还包括:
对所述系数向量进行归一化处理。
9.如权利要求1所述的系数计算方法,其特征在于:
所述获取Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息为:获取当前进行信道估计的符号对应的Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息,所述位置信息和权重信息为对上一个符号进行系数计算后更新得到,并且,初始化时,所述位置信息和权重信息为空集;
在得到当前进行信道估计的符号对应的系数向量之后,还包括:对Y1,Y2,…YQ所涉及的携带RS信息的符号位置信息和权重信息进行更新,以供下一个符号的系数计算使用。
10.如权利要求9所述的系数计算方法,其特征在于,所述对Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息进行更新,包括:
从q=2开始,对W(temp)进行如下的更新,直到q=Q:
从p=1开始,对W(temp)和D(temp)进行如下的更新,直到p=P:
11.如权利要求1所述的系数计算方法,其特征在于,还包括:
计算I中有效元素所对应的符号与当前进行信道估计的符号之间的距离向量;
根据所述距离向量确定当前进行信道估计的符号所属的场景;
系数计算进入稳定状态后,将计算得到的系数向量记录为相应场景的滤波器系数。
12.一种时域信道估计滤波器的系数计算装置,应用于LTE***中,其特征在于,包括:
第一获取模块,用于获取X1,X2…XP所对应的携带参考信号RS信息的符号的位置信息和权重信息,X1,X2…XP为与当前输入符号最接近的P个频域信道估计值,P为大于0的整数;
第二获取模块,用于获取Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息,Y1,Y2,…YQ为与当前进行信道估计的符号最接近的Q个时域信道估计值,Q为大于0的整数;
自相关计算模块,用于根据所述位置信息和权重信息,计算向量I的自相关矩阵,其中,I=[X1,X2…XP,Y1,Y2,…YQ];
互相关计算模块,用于根据所述位置信息和权重信息,计算向量O与向量 I之间的互相关矩阵,其中,O=[Y0],Y0为当前进行信道估计的符号的时域信道的理想值;
滤波系数计算模块,用于计算所述互相关矩阵与所述自相关矩阵的逆矩阵之积,得到当前进行信道估计的符号对应的系数向量。
17.如权利要求15或16所述的系数计算装置,其特征在于,所述相关函数为:
R(τi,j)=J0(2πfmaxτi,j)
其中,J0(·)表示第一类0阶贝塞尔函数,fmax为最大多普勒频移,τi,j=Tsym|i-j|,Tsym为一个符号的时间长度。
18.如权利要求12所述的系数计算装置,其特征在于,所述第二获取模块进一步用于:
获取当前进行信道估计的符号对应的Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息,所述位置信息和权重信息为对上一个符号进行系数计算后更新得到,并且,初始化时,所述位置信息和权重信息为空集;
在得到当前进行信道估计的符号对应的系数向量之后,对Y1,Y2,…YQ所涉及的携带RS信息的符号位置信息和权重信息进行更新,以供下一个符号的系数计算使用。
19.如权利要求18所述的系数计算装置,其特征在于,所述第二获取模块进一步用于,按照如下方式对Y1,Y2,…YQ所涉及的携带RS信息的符号的位置信息和权重信息进行更新:
从q=2开始,对W(temp)进行如下的更新,直到q=Q:
从p=1开始,对W(temp)和D(temp)进行如下的更新,直到p=P:
20.如权利要求12所述的系数计算装置,其特征在于,还包括滤波系数选择模块,用于:
计算I中有效元素所对应的符号与当前进行信道估计的符号之间的距离向量;
根据所述距离向量确定当前进行信道估计的符号所属的场景;
系数计算进入稳定状态后,将计算得到的系数向量记录为相应场景的滤波器系数。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010106242966A CN102123127B (zh) | 2010-12-31 | 2010-12-31 | 一种时域信道估计滤波器的系数计算方法和装置 |
EP11854492.3A EP2661041A4 (en) | 2010-12-31 | 2011-12-30 | Coefficient calculation method and device for time-domain channel estimation filter |
PCT/CN2011/085002 WO2012089156A1 (zh) | 2010-12-31 | 2011-12-30 | 一种时域信道估计滤波器的系数计算方法和装置 |
US13/977,375 US8908813B2 (en) | 2010-12-31 | 2011-12-30 | Coefficient calculation method and device for time-domain channel estimation filter |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010106242966A CN102123127B (zh) | 2010-12-31 | 2010-12-31 | 一种时域信道估计滤波器的系数计算方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102123127A CN102123127A (zh) | 2011-07-13 |
CN102123127B true CN102123127B (zh) | 2013-03-13 |
Family
ID=44251583
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010106242966A Expired - Fee Related CN102123127B (zh) | 2010-12-31 | 2010-12-31 | 一种时域信道估计滤波器的系数计算方法和装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US8908813B2 (zh) |
EP (1) | EP2661041A4 (zh) |
CN (1) | CN102123127B (zh) |
WO (1) | WO2012089156A1 (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102123127B (zh) * | 2010-12-31 | 2013-03-13 | 意法·爱立信半导体(北京)有限公司 | 一种时域信道估计滤波器的系数计算方法和装置 |
CN102892152B (zh) * | 2011-07-19 | 2017-11-03 | 中兴通讯股份有限公司 | 一种自适应传输方法及*** |
US9414372B2 (en) * | 2012-03-16 | 2016-08-09 | Qualcomm Incorporated | Digital filter control for filter tracking speedup |
KR102266595B1 (ko) * | 2014-05-02 | 2021-06-18 | 삼성전자주식회사 | 무선 통신 시스템에서 채널 추정 방법 및 장치 |
US9596102B2 (en) * | 2014-09-16 | 2017-03-14 | Samsung Electronics Co., Ltd. | Computing system with channel estimation mechanism and method of operation thereof |
JP6197916B1 (ja) * | 2016-06-03 | 2017-09-20 | Nttエレクトロニクス株式会社 | 補償係数算出方法 |
CN108075992B (zh) * | 2016-11-09 | 2020-05-15 | 电信科学技术研究院 | 一种信道估计方法及装置 |
US10305555B2 (en) | 2017-10-20 | 2019-05-28 | Micron Technology, Inc. | Autocorrelation and memory allocation for wireless communication |
WO2019172412A1 (ja) * | 2018-03-09 | 2019-09-12 | 日本電気株式会社 | 無線装置及びチャネル予測方法 |
CN110691052B (zh) * | 2018-07-06 | 2021-11-23 | 扬智科技股份有限公司 | 符号同步方法与信号接收电路 |
CN111935746B (zh) * | 2020-08-14 | 2024-01-09 | Oppo广东移动通信有限公司 | 获取通信参数的方法、装置、终端及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101534266A (zh) * | 2009-04-14 | 2009-09-16 | 北京天碁科技有限公司 | 一种正交频分复用***的信道估计方法和装置 |
CN101753498A (zh) * | 2008-12-05 | 2010-06-23 | 中兴通讯股份有限公司 | 正交频分复用信道估计结果的滤波方法与装置 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2264921B1 (en) * | 2008-04-10 | 2020-08-12 | Panasonic Intellectual Property Management Co., Ltd. | Receiving apparatus, receiving method, integrated circuit, digital television receiver, and program |
US8761274B2 (en) * | 2009-02-04 | 2014-06-24 | Acorn Technologies, Inc. | Least squares channel identification for OFDM systems |
CN101795246B (zh) * | 2010-01-07 | 2012-12-05 | 北京天碁科技有限公司 | 信道估计方法和信道估计装置 |
CN101815042B (zh) * | 2010-04-13 | 2012-12-05 | 新邮通信设备有限公司 | 一种正交频分复用***信道估计方法和装置 |
US8520781B2 (en) * | 2010-07-29 | 2013-08-27 | Blackberry Limited | System and method for channel estimation |
CN102123127B (zh) * | 2010-12-31 | 2013-03-13 | 意法·爱立信半导体(北京)有限公司 | 一种时域信道估计滤波器的系数计算方法和装置 |
-
2010
- 2010-12-31 CN CN2010106242966A patent/CN102123127B/zh not_active Expired - Fee Related
-
2011
- 2011-12-30 WO PCT/CN2011/085002 patent/WO2012089156A1/zh active Application Filing
- 2011-12-30 US US13/977,375 patent/US8908813B2/en not_active Expired - Fee Related
- 2011-12-30 EP EP11854492.3A patent/EP2661041A4/en not_active Withdrawn
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101753498A (zh) * | 2008-12-05 | 2010-06-23 | 中兴通讯股份有限公司 | 正交频分复用信道估计结果的滤波方法与装置 |
CN101534266A (zh) * | 2009-04-14 | 2009-09-16 | 北京天碁科技有限公司 | 一种正交频分复用***的信道估计方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
EP2661041A4 (en) | 2017-06-14 |
CN102123127A (zh) | 2011-07-13 |
EP2661041A1 (en) | 2013-11-06 |
US8908813B2 (en) | 2014-12-09 |
WO2012089156A1 (zh) | 2012-07-05 |
US20140056393A1 (en) | 2014-02-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102123127B (zh) | 一种时域信道估计滤波器的系数计算方法和装置 | |
CN108234364B (zh) | 一种lte-a***中基于小区参考信号的信道估计方法 | |
JP5152325B2 (ja) | 直交周波数分割多重におけるチャネル推定の方法及びそのシステム | |
US8130852B2 (en) | Method for estimating channel in radio communication system and device therefor | |
CN100499610C (zh) | 一种基于正交序列设计的低复杂度信道估计方法 | |
US20070098090A1 (en) | Systems, methods, and apparatus for mitigation of nonlinear distortion | |
CN110113286B (zh) | 一种基于正交匹配追踪的低复杂度水声信道估计算法 | |
CN103685088B (zh) | 稀疏信道的导频优化方法、装置和信道估计方法 | |
CN101827057B (zh) | 用于正交频分复用通信***的信道估计方法及信道估计器 | |
WO2010118588A1 (zh) | 一种正交频分复用***的信道估计方法和装置 | |
US7218693B2 (en) | Method and apparatus for channel estimation in a packet based transmission system having reduced complexity | |
JPWO2007020943A1 (ja) | Ofdm通信方法 | |
WO2011054279A1 (zh) | 基于多相分解的多天线信道估计方法 | |
CN101494468A (zh) | 多小区联合信道估计方法和装置 | |
Liu et al. | Iterative compensated MMSE channel estimation in LTE systems | |
JP5723371B2 (ja) | Ofdm通信システムにおいて通信路を推定するためのプロセス、およびそのための受信機 | |
CN102164110A (zh) | 频域均衡方法和*** | |
EP2536082A1 (en) | Iterative channel estimation method with compensated MMSE | |
CN102790746A (zh) | Ofdm***的信道估计方法 | |
JP4213747B2 (ja) | 複雑度を低減したスライディングウィンドウベースの等化器 | |
JP6110641B2 (ja) | 送信元と受信機との間のチャネルインパルス応答を推定する方法及びデバイス | |
CN113055318B (zh) | 一种信道估计方法 | |
CN113099427B (zh) | 一种基于稀疏度自适应的带内部署NB-IoT信号恢复方法 | |
KR101853184B1 (ko) | 하나 이상의 수신된 무선 신호를 처리하는 장치 및 방법 | |
CN104753834B (zh) | 一种信道估计方法和装置 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130313 Termination date: 20171231 |
|
CF01 | Termination of patent right due to non-payment of annual fee |