CN105337909B - 一种频偏估计的方法和装置 - Google Patents

一种频偏估计的方法和装置 Download PDF

Info

Publication number
CN105337909B
CN105337909B CN201410387070.7A CN201410387070A CN105337909B CN 105337909 B CN105337909 B CN 105337909B CN 201410387070 A CN201410387070 A CN 201410387070A CN 105337909 B CN105337909 B CN 105337909B
Authority
CN
China
Prior art keywords
frequency offset
current subframe
estimation value
offset estimation
value
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.)
Active
Application number
CN201410387070.7A
Other languages
English (en)
Other versions
CN105337909A (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.)
ZTE Corp
Original Assignee
ZTE Corp
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 ZTE Corp filed Critical ZTE Corp
Priority to CN201410387070.7A priority Critical patent/CN105337909B/zh
Priority to EP14899441.1A priority patent/EP3179682A4/en
Priority to PCT/CN2014/092359 priority patent/WO2016019657A1/zh
Priority to US15/329,065 priority patent/US9900184B2/en
Publication of CN105337909A publication Critical patent/CN105337909A/zh
Application granted granted Critical
Publication of CN105337909B publication Critical patent/CN105337909B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/0224Channel estimation using sounding signals
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2657Carrier synchronisation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2668Details of algorithms
    • H04L27/2669Details of algorithms characterised by the domain of operation
    • H04L27/2671Time domain
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2668Details of algorithms
    • H04L27/2673Details of algorithms characterised by synchronisation parameters
    • H04L27/2675Pilot or known symbols
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2668Details of algorithms
    • H04L27/2681Details of algorithms characterised by constraints
    • H04L27/2682Precision
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2668Details of algorithms
    • H04L27/2681Details of algorithms characterised by constraints
    • H04L27/2686Range of frequencies or delays tested
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2689Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
    • H04L27/2695Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation with channel estimation, e.g. determination of delay spread, derivative or peak tracking

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Circuits Of Receivers In General (AREA)

Abstract

本发明实施例公开了一种频偏估计的方法和装置,该方法包括:通过接收的探测参考信号SRS或Preamble码进行初始信道估计,获取当前子帧的第一频偏估计值;确定所述第一频偏估计值所处的频偏区间,并根据预设的频偏区间与初始频偏值之间的对应关系,获取所述第一频偏估计值对应的初始频偏值;基于接收的当前子帧包含的导频序列确定相邻导频位置的相位差,并在所述第一频偏估计值对应的初始频偏值限定的范围内,根据所述相位差获取双导频频偏估计值;根据所述第一频偏估计值对应的初始频偏值和所述双导频频偏估计值获取当前子帧的频偏估计值。

Description

一种频偏估计的方法和装置
技术领域
本发明涉及通信信号处理技术,尤其涉及一种频偏估计的方法和装置。
背景技术
正交频分复用(OFDM,Orthogonal Frequency Division Multiplexing)技术是一种子载波相互正交的多载波传输技术,这种技术要求子载波之间相互正交,并且对载波频率偏移非常敏感;所以,载波频率偏移问题成为当前OFDM技术的核心问题之一。
目前对于解决频偏较小的问题,已有很多精度高并且实现简单的方法;但是对于大频偏的问题,解决方案的精度都相对比较低,即使有些方法提高了频偏估计的精度,但同时也增加了实现的复杂度。
发明内容
为解决上述技术问题,本发明实施例期望提供一种频偏估计的方法和装置,针对大频偏的问题,既保证频偏估计精度,而且还降低了对于大频偏估计的复杂度。
本发明的技术方案是这样实现的:
第一方面,本发明实施例提供了一种频偏估计的方法,该方法可以包括:
通过接收的探测参考信号(SRS,Sounding Reference Signal)或Preamble码进行初始信道估计,获取当前子帧的第一频偏估计值;
确定所述第一频偏估计值所处的频偏区间,并根据预设的频偏区间与初始频偏值之间的对应关系,获取所述第一频偏估计值对应的初始频偏值;
基于接收的当前子帧包含的导频序列确定相邻导频位置的相位差,并在所述第一频偏估计值对应的初始频偏值限定的范围内,根据所述相位差获取双导频频偏估计值;
根据所述第一频偏估计值对应的初始频偏值和所述双导频频偏估计值获取当前子帧的频偏估计值。
进一步地,所述通过接收的SRS或Preamble码进行初始信道估计,获取当前子帧的第一频偏估计值,包括:
通过接收当前子帧的SRS或者当前子帧的Preamble码来获取信道估计值Hu和接收数据Yu
根据所述信道估计值Hu和所述接收数据Yu,并结合本地导频序列Pu和下式获取接收天线kaRx的当前子帧的时域采样值wu(n,kaRx):
Figure GDA0002255942960000021
其中,n表示采样点序号,u表示所述当前子帧对应的用户设备UE,k表示所述UE对应的子载波索引,kaRx表示接收天线索引,
Figure GDA0002255942960000023
表示所述UE所占用的资源块(RB,Resource Block)的大小,*表示共轭运算,IDFT表示傅里叶逆变换;
根据所述接收天线kaRx的当前子帧的时域采样值wu(n,kaRx)的相位信息获取接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx)。
进一步地,所述根据所述接收天线kaRx的当前子帧的时域采样值wu(n,kaRx)的相位信息获取接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx),包括:
将所述当前子帧的时域采样值wu(n,kaRx)的前半部分和后半部分的相位值分别进行平均后,获取前、后两个部分之间的相位差
Figure GDA0002255942960000024
根据下式获取所述接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx):
Figure GDA0002255942960000025
进一步地,所述方法还包括:
对所述接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx)进行平滑处理,得到平滑后的当前子帧的第一频偏估计值f1 (u)(kaRx)。
进一步地,所述基于接收的当前子帧包含的导频序列确定相邻导频位置的相位差,并在所述第一频偏估计值对应的初始频偏值限定的范围内,根据所述相位差获取双导频频偏估计值,包括:
根据所述接收天线kaRx的当前子帧包含的两个导频序列所得到的两个导频信道估计值
Figure GDA0002255942960000031
Figure GDA0002255942960000032
并按照下式获取接收天线kaRx的两个导频位置的相位差
Figure GDA0002255942960000034
其中,*表示共轭运算符,angle()表示求角度运算符;
通过所述两个导频序列的时间间隔Δt以及所述两个导频位置的相位差
Figure GDA0002255942960000035
按照下式获取接收天线kaRx的双导频频偏估计结果Δf2 (u)(kaRx):
Figure GDA0002255942960000036
在所述初始频偏值
Figure GDA0002255942960000037
限定的范围内,根据双导频频偏估计结果Δf2 (u)(kaRx)以及下式获取所述接收天线kaRx的双导频频偏估计值ΔΔf2(kaRx):
Figure GDA0002255942960000038
其中,T为所述双导频频偏估计值的最大范围值;l为预设整数,用于将所述双导频频偏估计值限定在预定的范围内。
进一步地,在根据所述第一频偏估计值对应的初始频偏值和所述双导频频偏估计值获取当前子帧的频偏估计值之后,所述方法还包括:
对当前子帧的频偏估计值fu(kaRx)通过滤波的方式进行平滑处理,得到平滑后的当前子帧的频偏估计值f(u)(kaRx)。
第二方面,本发明实施例提供了一种频偏估计的装置,所述装置包括:第一获取单元、第二获取单元、第三获取单元和第四获取单元,其中,
所述第一获取单元,用于通过接收的SRS或Preamble码进行初始信道估计,获取当前子帧的第一频偏估计值;
所述第二获取单元,用于确定所述第一频偏估计值所处的频偏区间,并根据预设的频偏区间与初始频偏值之间的对应关系,获取所述第一频偏估计值对应的初始频偏值;
所述第三获取单元,用于基于接收的所述当前子帧包含的导频序列确定相邻导频位置的相位差,并在所述第一频偏估计值对应的初始频偏值限定的范围内,根据所述相位差获取双导频频偏估计值;
所述第四获取单元,用于根据所述第一频偏估计值对应的初始频偏值和所述双导频频偏估计值获取当前子帧的频偏估计值。
进一步地,所述第一获取单元,用于:
通过接收当前子帧的SRS或者当前子帧的Preamble码来获取信道估计值Hu和接收数据Yu
以及,根据所述信道估计值Hu和所述接收数据Yu,并结合本地导频序列Pu和下式获取接收天线kaRx的当前子帧的时域采样值wu(n,kaRx):
其中,n表示采样点序号,u表示所述当前子帧对应的用户设备UE,k表示所述UE对应的子载波索引,kaRx表示接收天线索引,
Figure GDA0002255942960000043
表示所述UE所占用的RB的大小,*表示共轭运算,IDFT表示傅里叶逆变换;
以及,根据所述接收天线kaRx的当前子帧的时域采样值wu(n,kaRx)的相位信息获取接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx)。
进一步地,所述第一获取单元,用于:
将所述当前子帧的时域采样值wu(n,kaRx)的前半部分和后半部分的相位值分别进行平均后,获取前、后两个部分之间的相位差
Figure GDA0002255942960000044
根据下式获取所述接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx):
Figure GDA0002255942960000045
进一步地,所述第一获取单元,还用于对所述接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx)进行平滑处理,得到平滑后的当前子帧的第一频偏估计值f1 (u)(kaRx)。
进一步地,所述第三获取单元,用于:
根据所述接收天线kaRx的当前子帧包含的两个导频序列所得到的两个导频信道估计值
Figure GDA0002255942960000051
Figure GDA0002255942960000052
并按照下式获取接收天线kaRx的两个导频位置的相位差
Figure GDA0002255942960000053
Figure GDA0002255942960000054
其中,*表示共轭运算符,angle()表示求角度运算符;
以及,通过所述两个导频序列的时间间隔Δt以及所述两个导频位置的相位差
Figure GDA0002255942960000055
按照下式获取接收天线kaRx的双导频频偏估计结果Δf2 (u)(kaRx):
Figure GDA0002255942960000056
以及,在所述初始频偏值
Figure GDA0002255942960000057
限定的范围内,根据双导频频偏估计结果Δf2 (u)(kaRx)以及下式获取所述接收天线kaRx的双导频频偏估计值ΔΔf2(kaRx):
Figure GDA0002255942960000058
其中,T为所述双导频频偏估计值的最大范围值;l为预设整数,用于将所述双导频频偏估计值限定在预定的范围内。
进一步地,所述第四获取单元,还用于对当前子帧的频偏估计值fu(kaRx)通过滤波的方式进行平滑处理,得到平滑后的当前子帧的频偏估计值f(u)(kaRx)。
本发明实施例提供了一种频偏估计的方法和装置,通过根据UE发送的SRS或Preamble码所确定的第一频偏估计值得到初始频偏值,并且在初始频偏值的范围内结合双导频频偏估计值获取当前子帧的频偏估计值,既保证频偏估计精度,而且还降低了对于大频偏估计的复杂度。
附图说明
图1为本发明实施例提供的一种频偏估计的方法流程示意图;
图2为本发明实施例提供的一种频偏估计的装置结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。
参见图1,其示出了本发明实施例提供的一种频偏估计的方法流程,需要说明的是,为了能够清楚的对本发明实施例的技术方案进行说明,本实施例中以一个UE为例进行说明,可以理解的,本领域技术人员可以无需创造性劳动将一个UE为例的实施例应用于多个UE的实施例中。该方法可以包括:
S101:通过接收用户设备(UE,User Equipment)发送的SRS或Preamble码进行初始信道估计,获取所述UE当前子帧的第一频偏估计值;
示例性地,步骤S101具体可以包括:
首先,通过接收由UE发送的当前子帧的SRS或者当前子帧的Preamble码来获取所述UE的信道估计值Hu和接收数据Yu;具体的获取过程为本领域技术人员的常规实现手段,在此不再赘述;
然后,根据所述UE的信道估计值Hu和接收数据Yu,并结合本地导频序列Pu和式1获取所述UE在接收天线kaRx的当前子帧的时域采样值wu(n,kaRx):
Figure GDA0002255942960000061
其中,n表示采样点序号,u表示所述UE,k表示所述UE对应的子载波索引,
Figure GDA0002255942960000062
kaRx表示接收天线索引,
Figure GDA0002255942960000063
表示所述UE所占用的资源块(RB,ResourceBlock)的大小,*表示共轭运算,IDFT表示傅里叶逆变换;
接着,根据所述UE在接收天线kaRx的当前子帧的时域采样值wu(n,kaRx)的相位信息获取所述UE在接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx);详细地,在本实施例中,可以将所述当前子帧的时域采样值wu(n,kaRx)的前半部分和后半部分的相位值分别进行平均后,获取前、后两个部分之间的相位差
Figure GDA0002255942960000071
也可以对所述当前子帧的时域采样值wu(n,kaRx)中相邻采样点之间的相位差进行平均来得到随后根据式2获取所述UE在接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx):
Figure GDA0002255942960000073
优选地,可以对式2得到的f0 u(kaRx)进行平滑处理,得到平滑后的当前子帧的第一频偏估计值f0 (u)(kaRx);具体地,可以将当前子帧以前的子帧对应的第一频偏估计值与f0 u(kaRx)进行加权求和;可以理解的,时间上越接近当前子帧则权值越大;时间上越远离当前子帧则权值越小。
在本实施例中,可以根据当前子帧的前一个子帧的对应的第一频偏估计值
Figure GDA0002255942960000074
与f0 u(kaRx)按照式3进行加权求和来对f0 u(kaRx)进行平滑处理,得到平滑后的当前子帧的第一频偏估计值f0 (u)(kaRx):
Figure GDA0002255942960000075
其中,p1为滤波因子,且0≤p1≤1。
可以理解的,f1 (u)(kaRx)仅较f0 u(kaRx)更加的平滑,但是f0 u(kaRx)也同样可以用于后续的处理,而将f1 (u)(kaRx)用于后续的处理仅是本实施例的一个更为优选的方案,所以,后续处理过程中所使用的第一频偏估计值既可以选择f0 u(kaRx),也可以选择f1 (u)(kaRx),本发明实施例对此不作具体的限定。
S102:确定第一频偏估计值所处的频偏区间范围,并根据预设的频偏区间范围与初始频偏值之间的对应关系,获取第一频偏估计值对应的初始频偏值;
示例性地,预设的所述频偏区间范围与初始频偏值之间的对应关系可以是一张具体的频偏范围区间的上、下限与初始频偏值的映射表,如表1所示,
表1
区间编号 上限:Hz 下限:Hz 初始频偏值:Hz
0 T1 T2 f0
1 T3 T4 f1
2 T5 T6 f2
3 T7 T8 f3
4 T9 T10 f4
5 T11 T12 f5
由表1可知,可以将所述第一频偏估计值与表1中的各频偏范围区间的上、下限进行比较,确定第一频偏估计值所处的频偏范围区间,然后还可以相应地获取到所述第一频偏估计值对应的初始频偏值。在本实施例中,与所述UE在接收天线kaRx的第一频偏估计值所对应的初始频偏值可以用
Figure GDA0002255942960000081
表示。
S103:基于接收的当前子帧包含的导频序列确定相邻导频位置的相位差,并在所述第一频偏估计值对应的初始频偏值限定的范围内,根据所述相位差获取双导频频偏估计值。
具体地,在本实施例中,S103具体可以包括:
首先,根据所述UE在接收天线kaRx的当前子帧包含的两个导频序列所得到的两个导频信道估计值
Figure GDA0002255942960000082
Figure GDA0002255942960000083
按照式3获取所述UE在接收天线kaRx的两个导频位置的相位差
其中,*表示共轭运算符,angle()表示求角度运算符。
然后,通过两个导频序列的时间间隔Δt以及所述两个导频位置的相位差
Figure GDA0002255942960000086
按照式4获取所述UE在接收天线kaRx的双导频频偏估计结果Δf1 (u)(kaRx):
Figure GDA0002255942960000087
接着,在初始频偏值
Figure GDA0002255942960000088
限定的范围内,根据双导频频偏估计结果Δf2 (u)(kaRx)以及式5获取所述UE在接收天线kaRx的双导频频偏估计值ΔΔf2(kaRx):
Figure GDA0002255942960000091
其中,T为双导频频偏估计值的最大范围值;l为预设整数,用于将双导频频偏估计值限定在预定的范围内。
S104:根据所述第一频偏估计值对应的初始频偏值和所述双导频频偏估计值得到当前子帧的频偏估计值。
示例性地,在本实施例中,可以通过将所述初始频偏值
Figure GDA0002255942960000092
与所述双导频频偏估计值ΔΔf2(kaRx)相加获取所述UE在接收天线kaRx的当前子帧的频偏估计值fu(kaRx)。
优选地,在获取到所述UE在接收天线kaRx的当前子帧的频偏估计值fu(kaRx)之后,可以对fu(kaRx)通过滤波的方式进行平滑处理,得到平滑后的当前子帧的频偏估计值f(u)(kaRx);
具体地,本实施例中,可以通过所述UE在接收天线kaRx的前一子帧的频偏估计值
Figure GDA0002255942960000093
与fu(kaRx)按照式6进行加权求和的方式来得到平滑后的当前子帧的频偏估计值f(u)(kaRx):
其中,p2是滤波因子,且0≤p2≤1。
上述方案为本发明实施例提供的一种频偏估计的方法,通过根据UE发送的SRS或Preamble码所确定的第一频偏估计值得到初始频偏值,并且在初始频偏值的范围内结合双导频频偏估计值获取当前子帧的频偏估计值,既保证频偏估计精度,而且还降低了对于大频偏估计的复杂度。
基于图1所示实施例相同的技术构思,参见图2,其示出了本发明实施例提供的一种频偏估计的装置20,该装置20可以包括:第一获取单元201、第二获取单元202、第三获取单元203和第四获取单元204,其中,
所述第一获取单元201,用于通过接收UE发送的SRS或Preamble码进行初始信道估计,获取所述UE当前子帧的第一频偏估计值;
所述第二获取单元202,用于确定所述第一频偏估计值所处的频率范围,并根据所述频率范围以及预设的所述频率范围与初始频偏值之间的对应关系,获取所述第一频偏估计值对应的初始频偏值;
所述第三获取单元203,用于基于接收的所述UE当前子帧包含的两个导频序列计算两个导频位置的相位差,并在所述初始频偏值限定的范围内,根据所述相位差获取双导频频偏估计值;
所述第四获取单元204,用于根据所述初始频偏值和所述双导频频偏估计值获取当前子帧的频偏估计值。
示例性地,所述第一获取单元201,用于:
通过接收由UE发送的当前子帧的SRS或者当前子帧的Preamble码来获取所述UE的信道估计值Hu和接收数据Yu
以及,根据所述UE的信道估计值Hu和接收数据Yu,并结合本地导频序列Pu和下式获取所述UE在接收天线kaRx的当前子帧的时域采样值wu(n,kaRx):
Figure GDA0002255942960000101
其中,n表示采样点序号,u表示所述当前子帧对应的用户设备UE,k表示所述UE对应的子载波索引,
Figure GDA0002255942960000102
kaRx表示接收天线索引,
Figure GDA0002255942960000103
表示所述UE所占用的RB的大小,*表示共轭运算,IDFT表示傅里叶逆变换;
以及,根据所述UE在接收天线kaRx的当前子帧的时域采样值wu(n,kaRx)的相位信息获取所述UE接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx)。
具体地,所述第一获取单元201,用于:
将所述当前子帧的时域采样值wu(n,kaRx)的前半部分和后半部分的相位值分别进行平均后,获取前、后两个部分之间的相位差
Figure GDA0002255942960000104
根据下式获取所述UE在接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx):
Figure GDA0002255942960000111
优选地,所述第一获取单元201,还用于对所述UE在接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx)进行平滑处理,得到平滑后的当前子帧的第一频偏估计值f1 (u)(kaRx);具体地,所述第一获取单元201可以将当前子帧以前的子帧对应的第一频偏估计值与f0 u(kaRx)进行加权求和;可以理解的,时间上越接近当前子帧则权值越大;时间上越远离当前子帧则权值越小。
在本实施例中,所述第一获取单元201可以根据当前子帧的前一个子帧的对应的第一频偏估计值
Figure GDA0002255942960000112
与f0 u(kaRx)按照下式进行加权求和来对f0 u(kaRx)进行平滑处理,得到平滑后的当前子帧的第一频偏估计值f1 (u)(kaRx):
Figure GDA0002255942960000113
其中,p1为滤波因子,且0≤p1≤1。
可以理解的,f1 (u)(kaRx)仅较f0 u(kaRx)更加的平滑,但是f0 u(kaRx)也同样可以用于后续模块进行处理,而将f1 (u)(kaRx)用于后续模块的处理仅是本实施例的一个更为优选的方案,所以,后续模块处理过程中所使用的第一频偏估计值既可以选择f0 u(kaRx),也可以选择f1 (u)(kaRx),本发明实施例对此不作具体的限定。
示例性地,预设的所述频率范围与初始频偏值之间的对应关系可以是一张具体的频率范围区间的上、下限与初始频偏值的映射表,如表1所示,
具体地,由表1可知,第二获取单元202可以用于将所述第一频偏估计值与表1中的各频率范围区间的上、下限进行比较,确定第一频偏估计值所处的频率范围区间,然后还可以相应地获取到所述第一频偏估计值对应的初始频偏值。在本实施例中,与所述UE在接收天线kaRx的第一频偏估计值所对应的初始频偏值可以用
Figure GDA0002255942960000114
表示。
示例性地,所述第三获取单元203,用于:
根据所述UE在接收天线kaRx的当前子帧包含的两个导频序列所得到的两个导频信道估计值
Figure GDA0002255942960000121
Figure GDA0002255942960000122
按照下式获取所述UE在接收天线kaRx的两个导频位置的相位差
Figure GDA0002255942960000123
其中,*表示共轭运算符,angle()表示求角度运算符;
以及,通过所述两个导频序列的时间间隔Δt以及所述两个导频位置的相位差
Figure GDA0002255942960000125
按照下式获取所述UE在接收天线kaRx的双导频频偏估计结果Δf2 (u)(kaRx):
以及,在所述初始频偏值限定的范围内,根据双导频频偏估计结果Δf2 (u)(kaRx)以及下式获取所述UE在接收天线kaRx的双导频频偏估计值ΔΔf2(kaRx):
Figure GDA0002255942960000127
其中,T为所述双导频频偏估计值的最大范围值;l为预设整数,用于将所述双导频频偏估计值限定在预定的范围内。
具体地,在本实施例中,第四获取单元204可以通过将所述初始频偏值
Figure GDA0002255942960000128
与所述双导频频偏估计值ΔΔf2(kaRx)相加获取所述UE接收天线kaRx的当前子帧的频偏估计值fu(kaRx)。
优选地,在获取到所述UE接收天线kaRx的当前子帧的频偏估计值fu(kaRx)之后,所述第四获取单元204,还用于对当前子帧的频偏估计值fu(kaRx)通过滤波的方式进行平滑处理,得到平滑后的当前子帧的频偏估计值f(u)(kaRx)。
具体地,本实施例中,第四获取单元204可以通过所述UE接收天线kaRx的前一子帧的频偏估计值
Figure GDA0002255942960000129
与fu(kaRx)按照下式进行加权求和的方式来得到平滑后的当前子帧的频偏估计值f(u)(kaRx):
Figure GDA00022559429600001210
其中,p2是滤波因子,且0≤p2≤1。
本实施例提供了一种频偏估计的装置20,通过根据UE发送的SRS或Preamble码所确定的第一频偏估计值得到初始频偏值,并且在初始频偏值的范围内结合双导频频偏估计值获取当前子帧的频偏估计值,既保证频偏估计精度,而且还降低了对于大频偏估计的复杂度。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用硬件实施例、软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述,仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。

Claims (8)

1.一种频偏估计的方法,其特征在于,所述方法包括:
通过接收当前子帧的探测参考信号(SRS)或者当前子帧的Preamble码来获取信道估计值Hu和接收数据Yu
根据所述信道估计值Hu和所述接收数据Yu,并结合本地导频序列Pu和下式获取接收天线kaRx的当前子帧的时域采样值wu(n,kaRx):
其中,n表示采样点序号,u表示所述当前子帧对应的用户设备(UE),k表示所述UE对应的子载波索引,
Figure FDA0002255942950000011
kaRx表示接收天线索引,
Figure FDA0002255942950000012
表示所述UE所占用的资源块RB的大小,*表示共轭运算,IDFT表示傅里叶逆变换;
根据所述接收天线kaRx的当前子帧的时域采样值wu(n,kaRx)的相位信息获取接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx);
确定所述第一频偏估计值所处的频偏区间,并根据预设的频偏区间与初始频偏值之间的对应关系,获取所述第一频偏估计值对应的初始频偏值;
基于接收的当前子帧包含的导频序列确定相邻导频位置的相位差,并在所述第一频偏估计值对应的初始频偏值限定的范围内,根据所述相位差获取双导频频偏估计值;
根据所述第一频偏估计值对应的初始频偏值和所述双导频频偏估计值获取当前子帧的频偏估计值。
2.根据权利要求1所述的方法,其特征在于,所述根据所述接收天线kaRx的当前子帧的时域采样值wu(n,kaRx)的相位信息获取接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx),包括:
将所述当前子帧的时域采样值wu(n,kaRx)的前半部分和后半部分的相位值分别进行平均后,获取前、后两个部分之间的相位差
根据下式获取所述接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx):
Figure FDA0002255942950000022
3.根据权利要求1所述的方法,其特征在于,所述方法还包括:
对所述接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx)进行平滑处理,得到平滑后的当前子帧的第一频偏估计值f1 (u)(kaRx)。
4.根据权利要求1所述的方法,其特征在于,在根据所述第一频偏估计值对应的初始频偏值和所述双导频频偏估计值获取当前子帧的频偏估计值之后,所述方法还包括:
对当前子帧的频偏估计值fu(kaRx)通过滤波的方式进行平滑处理,得到平滑后的当前子帧的频偏估计值f(u)(kaRx)。
5.一种频偏估计的装置,其特征在于,所述装置包括:第一获取单元、第二获取单元、第三获取单元和第四获取单元,其中,
所述第一获取单元,用于通过接收当前子帧的探测参考信号(SRS)或者当前子帧的Preamble码来获取信道估计值Hu和接收数据Yu
以及,根据所述信道估计值Hu和所述接收数据Yu,并结合本地导频序列Pu和下式获取接收天线kaRx的当前子帧的时域采样值wu(n,kaRx):
Figure FDA0002255942950000025
其中,n表示采样点序号,u表示所述当前子帧对应的用户设备(UE),k表示所述UE对应的子载波索引,
Figure FDA0002255942950000023
kaRx表示接收天线索引,
Figure FDA0002255942950000024
表示所述UE所占用的资源块RB的大小,*表示共轭运算,IDFT表示傅里叶逆变换;
以及,根据所述接收天线kaRx的当前子帧的时域采样值wu(n,kaRx)的相位信息获取接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx);
所述第二获取单元,用于确定所述第一频偏估计值所处的频偏区间,并根据预设的频偏区间与初始频偏值之间的对应关系,获取所述第一频偏估计值对应的初始频偏值;
所述第三获取单元,用于基于接收的所述当前子帧包含的导频序列确定相邻导频位置的相位差,并在所述第一频偏估计值对应的初始频偏值限定的范围内,根据所述相位差获取双导频频偏估计值;
所述第四获取单元,用于根据所述第一频偏估计值对应的初始频偏值和所述双导频频偏估计值获取当前子帧的频偏估计值。
6.根据权利要求5所述的装置,其特征在于,所述第一获取单元,用于:
将所述当前子帧的时域采样值wu(n,kaRx)的前半部分和后半部分的相位值分别进行平均后,获取前、后两个部分之间的相位差
Figure FDA0002255942950000031
根据下式获取所述接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx):
7.根据权利要求5所述的装置,其特征在于,所述第一获取单元,还用于对所述接收天线kaRx的当前子帧的第一频偏估计值f0 u(kaRx)进行平滑处理,得到平滑后的当前子帧的第一频偏估计值f1 (u)(kaRx)。
8.根据权利要求5所述的装置,其特征在于,所述第四获取单元,还用于对当前子帧的频偏估计值fu(kaRx)通过滤波的方式进行平滑处理,得到平滑后的当前子帧的频偏估计值f(u)(kaRx)。
CN201410387070.7A 2014-08-07 2014-08-07 一种频偏估计的方法和装置 Active CN105337909B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201410387070.7A CN105337909B (zh) 2014-08-07 2014-08-07 一种频偏估计的方法和装置
EP14899441.1A EP3179682A4 (en) 2014-08-07 2014-11-27 Method, device and storage medium for frequency offset estimate
PCT/CN2014/092359 WO2016019657A1 (zh) 2014-08-07 2014-11-27 一种频偏估计的方法和装置、存储介质
US15/329,065 US9900184B2 (en) 2014-08-07 2014-11-27 Method, device and storage medium for frequency offset estimate

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410387070.7A CN105337909B (zh) 2014-08-07 2014-08-07 一种频偏估计的方法和装置

Publications (2)

Publication Number Publication Date
CN105337909A CN105337909A (zh) 2016-02-17
CN105337909B true CN105337909B (zh) 2020-02-14

Family

ID=55263077

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410387070.7A Active CN105337909B (zh) 2014-08-07 2014-08-07 一种频偏估计的方法和装置

Country Status (4)

Country Link
US (1) US9900184B2 (zh)
EP (1) EP3179682A4 (zh)
CN (1) CN105337909B (zh)
WO (1) WO2016019657A1 (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10439663B2 (en) * 2016-04-06 2019-10-08 Qualcomm Incorporated Methods and apparatus for phase noise estimation in data symbols for millimeter wave communications
CN107135182B (zh) * 2017-04-10 2020-06-30 上海顺久电子科技有限公司 调频信号的频偏计算方法及装置
CN109495414B (zh) * 2017-09-12 2020-07-03 电信科学技术研究院 一种频偏估计方法、装置、设备及计算机可读存储介质
CN109756435B (zh) * 2017-11-07 2021-10-15 北京中科晶上科技股份有限公司 一种对信号的频偏估计方法
CN110311874B (zh) * 2019-07-18 2022-02-15 锐迪科微电子(上海)有限公司 载波频偏确定方法及装置
CN114301743B (zh) * 2020-09-23 2023-10-27 紫光展锐(重庆)科技有限公司 频偏估计方法及装置、存储介质、终端
CN114500184B (zh) * 2020-10-23 2023-08-01 大唐移动通信设备有限公司 一种信道估计方法和装置及设备
CN114615122B (zh) * 2022-01-30 2023-11-03 北京邮电大学 通信信号的频偏确定方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101001231A (zh) * 2006-01-12 2007-07-18 上海原动力通信科技有限公司 宽带时分双工蜂窝系统的频偏校正方法及小区初搜方法
CN101309108A (zh) * 2008-07-10 2008-11-19 北京天碁科技有限公司 一种专用载波时分多媒体广播组播系统及通信方法
CN101547062A (zh) * 2009-03-03 2009-09-30 华为技术有限公司 频偏纠正方法和设备
CN102315908A (zh) * 2011-09-26 2012-01-11 华为技术有限公司 一种频偏估计方法和基站
CN102571647A (zh) * 2010-12-27 2012-07-11 中兴通讯股份有限公司 一种频偏估计方法及装置

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7027429B2 (en) * 2001-06-19 2006-04-11 Flarion Technologies, Inc. Method and apparatus for time and frequency synchronization of OFDM communication systems
US7257102B2 (en) * 2002-04-02 2007-08-14 Broadcom Corporation Carrier frequency offset estimation from preamble symbols
CN102347918B (zh) * 2010-07-30 2015-08-12 中兴通讯股份有限公司 精细频偏估计方法及装置
FR3010269B1 (fr) * 2013-09-04 2015-10-02 Commissariat Energie Atomique Recepteur fbmc a compensation d'offset de frequence porteuse
EP3197114B1 (en) * 2014-08-06 2022-01-19 Huawei Technologies Co., Ltd. Method and apparatus for transmitting uplink information in multi-user multiple-input multiple-output system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101001231A (zh) * 2006-01-12 2007-07-18 上海原动力通信科技有限公司 宽带时分双工蜂窝系统的频偏校正方法及小区初搜方法
CN101309108A (zh) * 2008-07-10 2008-11-19 北京天碁科技有限公司 一种专用载波时分多媒体广播组播系统及通信方法
CN101547062A (zh) * 2009-03-03 2009-09-30 华为技术有限公司 频偏纠正方法和设备
CN102571647A (zh) * 2010-12-27 2012-07-11 中兴通讯股份有限公司 一种频偏估计方法及装置
CN102315908A (zh) * 2011-09-26 2012-01-11 华为技术有限公司 一种频偏估计方法和基站

Also Published As

Publication number Publication date
US9900184B2 (en) 2018-02-20
US20170214552A1 (en) 2017-07-27
EP3179682A4 (en) 2017-08-16
EP3179682A1 (en) 2017-06-14
CN105337909A (zh) 2016-02-17
WO2016019657A1 (zh) 2016-02-11

Similar Documents

Publication Publication Date Title
CN105337909B (zh) 一种频偏估计的方法和装置
EP3174258A1 (en) Method and apparatus for generating and processing frequency division multiple waveform signal
US20180152276A1 (en) Method and Apparatus for Implementing Data Transmission
CN101989968A (zh) 频偏估计方法及装置
CN111988246B (zh) 一种广播信道解调参考信号检测方法、装置、设备和介质
CN105991484B (zh) 一种频偏估计的方法和装置
EP3179681B1 (en) Inter-carrier interference cancellation method, apparatus and computer storage media
CN110392003B (zh) 一种信号接收方法及装置
CN108075807B (zh) 一种小区标识的检测方法及装置
CN103475605B (zh) 一种3gpplte-a下行系统中基于用户专用参考信号的信道估计方法
CN101330488A (zh) 同步检测
TWI749036B (zh) 頻域載波間干擾補償方法與裝置、非暫時性電腦可讀媒體、製造處理器的方法及建構積體電路的方法
KR102093266B1 (ko) 신호 감지 방법 및 장치
WO2012171407A1 (zh) 一种确定时间同步位置的方法及设备
CN103856419B (zh) 对至eNB的上行信道进行信道估计的方法与设备
CN107018100B (zh) 一种时域去噪方法及装置
CN115412214B (zh) Rach信号的检测方法和装置
CN103248591A (zh) 基于频谱重心的粗频偏估计方法
CN111884958B (zh) 定时误差补偿方法、装置、设备和计算机可读存储介质
CN106921450B (zh) 一种信噪比估计方法及装置
CN104753834B (zh) 一种信道估计方法和装置
CN114301745A (zh) 载波频偏和采样频偏的确定方法及装置
CN112152950A (zh) Ofdm系统中的基于稀疏离散导频的信道估计方法和装置
CN106793144B (zh) 一种物理随机接入信道prach资源的分配方法和装置
US7876844B2 (en) Apparatus and method for determining positions of pilot sub-carriers in an OFDM symbol

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant