CN108366032B - 一种用于纠正大频偏的广义频分复用时频同步方法 - Google Patents

一种用于纠正大频偏的广义频分复用时频同步方法 Download PDF

Info

Publication number
CN108366032B
CN108366032B CN201810117264.3A CN201810117264A CN108366032B CN 108366032 B CN108366032 B CN 108366032B CN 201810117264 A CN201810117264 A CN 201810117264A CN 108366032 B CN108366032 B CN 108366032B
Authority
CN
China
Prior art keywords
frequency
sequence
sampling
value
correlation
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
CN201810117264.3A
Other languages
English (en)
Other versions
CN108366032A (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN201810117264.3A priority Critical patent/CN108366032B/zh
Publication of CN108366032A publication Critical patent/CN108366032A/zh
Application granted granted Critical
Publication of CN108366032B publication Critical patent/CN108366032B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/2657—Carrier synchronisation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • 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
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W56/00—Synchronisation arrangements
    • H04W56/0035—Synchronisation arrangements detecting errors in frequency or phase
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H04L2027/0024—Carrier regulation at the receiver end
    • H04L2027/0026—Correction of carrier offset

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Synchronisation In Digital Transmission Systems (AREA)

Abstract

本发明公开了一种用于纠正大频偏的广义频分复用时频同步方法，主要解决现有方法，受大频偏影响导致同步性能急剧下降的问题。具体步骤包括：(1)接收电信号；(2)对采样序列进行粗符号定时同步；(3)纠正采样序列的小数倍频偏；(4)选择路径候选定时时刻；(5)绘制二维时频度量曲面；(6)估计路径定时时刻(7)纠正采样序列的整数倍频偏；(8)估计第一径到达时刻。本发明存在大频偏时，载波频率同步和符号定时同步性能远远优于现有方法；本发明的频偏估计范围远大于现有广义频分复用GFDM系统同步方法。

Description

一种用于纠正大频偏的广义频分复用时频同步方法
技术领域
本发明属于通信技术领域，更进一步涉及无线通信技术领域中的一种用于纠正大频偏的广义频分复用GFDM(generalized frequency division multiplexing)时频同步方法。本发明可用于在广义频分复用GFDM系统中的载波频率同步和符号定时同步，提升系统在衰落信道下的同步性能。
背景技术
同步是信号均衡和解调的前提。同步误差对广义频分复用GFDM系统的误码率的影响远大于对正交频分复用OFDM(orthogonal frequency division multiplexing)系统的误码率的影响。
Vodafone GmbH提出的专利申请文献“GFDM radio transmission using apseudo circular preamble”(申请日：2014年12月12日，申请号：14/568570,公开号：US9236981B2)中公开了一种广义频分复用GFDM同步方法。该方法具体步骤为：第一步，在采用循环前缀和循环后缀的互相关得到尖锐的狄利克雷脉冲，狄利克雷脉冲对应的位置为符号定时同步点；第二步，获得定时同步点后，采用循环前缀和循环后缀的自相关获得小数倍频偏；第三步，去除循环前缀和循环后缀，利用解调的频域伪随机序列获得整数倍频偏。该方法存在的不足之处是，由于无线通信中信号传输受到信道衰落影响，接收端在计算循环前缀和循环后缀的互相关后，观察不到尖锐的狄利克雷脉冲，找不到正确的符号定时同步点。
Ivan S Gaspar等人在其发表的论文“A synchronization technique forgeneralized frequency division multiplexing”(Eurasip Journal on Advances inSignal Processing,2014,2014(1):67)中提出了一种适用于衰落信道的基于独立前导码的广义频分复用GFDM时频同步方法。该方法具体步骤为：第一步，发送端生成前后重复的前导码；第二步，接收端利用接收序列的自相关获得粗符号定时点和小数倍频偏，并纠正接收序列的小数倍频偏；第三步，计算纠正小数倍频偏后的接收序列和本地前导码的互相关；第四步，将未纠正小数倍频偏的接收序列的自相关与纠正小数倍频偏的接收序列和本地前导码的互相关相乘，获得最强路径定时点度量；第五步，找到最强路径定时点的位置，利用阈值准则搜索第一径定时点。该方法存在的不足之处是，第一，当采样序列中存在大频偏时，纠正小数倍频偏后的采样序列和本地前导码的互相关将会受到整数倍频偏的影响，导致第一径定时时刻不准确，符号定时同步性能急剧下降的；第二，Ivan S Gaspar的方法采用一个前导码序列仅能估计一个子载波带宽以内的频偏，导致频谱资源浪费。
发明内容
本发明的目的在于针对上述已有技术的不足，提出了一种用于纠正大频偏的广义频分复用时频同步方法。
实现本发明的思路是，将接收的广义频分复用GFDM模拟电信号处理后的采样序列，依次经粗符号定时同步、小数倍频偏纠正、整数倍频偏纠正，得到纠正大频偏后的无频偏采样序列，实现载波频率同步，从路径到达时刻向前搜索第一径到达时刻，实现符号定时同步。
本发明的具体步骤包括如下：
(1)接收电信号：
(1a)广义频分复用GFDM的接收机检测广义频分复用GFDM的发送机发送的模拟电信号；
(1b)对检测到的模拟电信号进行模数转换，得到实数字信号；
(1c)对实数字信号进行希尔伯特变换，得到复数信号；
(1d)对复数字信号进行数字下变频处理，得到采样序列；
(2)对采样序列进行粗符号定时同步：
(2a)利用自相关公式，计算采样序列中每一个采样点的自相关值，将所有的自相关值组成自相关序列；
(2b)利用能量值公式，计算采样序列中每一个采样点的能量值，将所有的能量值组成能量序列；
(2c)将自相关序列中的每一个自相关值依次作为截至自相关值，向前截取与循环前缀等长的子自相关序列，获得多个子自相关序列；其中，所述循环前缀的长度由广义频分复用GFDM系统参数决定；
(2d)将截至自相关值对应的采样点的序号作为子自相关序列的编号；
(2e)将能量序列中的每一个能量值依次作为截至能量值，向前截取与循环前缀等长的子能量序列，获得多个子能量序列；
(2f)将截至能量值对应的采样点的序号作为子能量序列的编号；
(2g)将具有相同编号的子自相关序列的自相关值和子能量序列的能量值进行相除操作，并将相除操作后的结果进行取绝对值操作，获得多个归一化子自相关序列；
(2h)对每一个归一化子自相关序列的子自相关值进行相加操作，得到对应采样点的粗符号定时度量值，将所有粗符号定时度量值组成粗符号定时度量序列；
(2i)找出粗符号定时度量序列中最大值对应的采样点，该采样点在采样序列中出现的时刻为粗符号定时同步时刻；
(3)纠正采样序列的小数倍频偏：
(3a)找出粗符号定时同步时刻对应的采样点的自相关值；
(3b)对粗符号定时同步时刻对应的采样点的自相关值进行取相位操作，获得该自相关值的相位，并将该自相关值的相位与圆周率进行相除操作，获得采样序列的小数倍频偏估计值；
(3c)利用小数倍频偏纠正公式，纠正采样序列的小数倍频偏，获得无小数倍频偏的采样序列；
(4)选择路径候选定时时刻：
(4a)对本地前导码序列进行共轭操作，得到共轭前导码序列；
(4b)将无小数倍频偏的采样序列中的每一个采样点依次作为起始点，向后截取与共轭前导码序列等长的子采样序列，将每一个子采样序列与共轭前导码序列进行相乘操作，获得多个子序列；
(4c)利用差分互相关公式，计算采样序列中每一个采样点的差分互相关值，将所有差分互相关值组成差分互相关序列；
(4d)对差分互相关序列进行取绝对值操作，将取绝对值操作后的差分互相关序列中的每一个差分互相关值进行平方操作，获得对应的路径候选定时度量值，将所有的路径候选定时度量值组成路径候选定时度量序列；
(4e)将路径候选定时度量序列，按照从大到小排列，找出前64个路径候选定时度量值对应的64个采样点，将64个采样点在采样序列中出现的时刻作为路径候选定时时刻；
(5)绘制二维时频度量曲面：
(5a)将64个路径候选定时时刻依次送入二维时频估计器；
(5b)二维时频估计器根据输入的路径候选定时时刻，找到该时刻的采样点，再找出该采样点对应的子序列；
(5c)对子序列进行快速傅里叶变换；
(5d)对快速傅里叶变换后的结果进行取绝对值操作，得到二维时频度量子序列；
(5e)判断64个路径候选定时时刻是否全部送入二维时频估计器，若是，则执行步骤(5f)，否则，执行步骤(5b)；
(5f)64个路径候选定时时刻全部送入二维时频估计器后，得到64个路径候选定时时刻对应的64个二维时频度量子序列，绘制由64个二维时频度量子序列构成的二维度时频度量曲面；
(6)估计路径定时时刻：
找出二维度时频度量曲面的最大值，将最大值所在的二维时频度量子序列对应的路径候选定时时刻，作为路径到达时刻；
(7)纠正采样序列的整数倍频偏：
(7a)找出二维度时频度量曲面的最大值对应的快速傅里叶变换的频率点值，将该频率点值作为采样序列的整数倍频偏估计值；
(7b)利用整数倍频偏纠正公式，纠正采样序列的整数倍频偏，得到纠正大频偏后的无频偏采样序列，实现载波频率同步；
(8)估计第一径到达时刻：
(8a)将无频偏采样序列中的每一个采样点依次作为起始点，向后截取与共轭前导码序列等长的无频偏子采样序列，获得多个无频偏子采样序列；
(8b)对每一个采样点对应的无频偏子采样序列与共轭前导码序列进行相乘操作，将相乘后的结果相加，得到一个互相关值；
(8c)将所有采样点对应的互相关值组成互相关序列；
(8d)将路径到达时刻对应的采样点所对应的互相值，作为截至互相关值；
(8e)在互相关序列中，从截至互相关值开始，向前截取与循环前缀等长的互相关子序列；
(8f)利用第一径定时估计阈值公式，计算第一径定时估计阈值；
(8g)对互相关子序列的每一个互相关值取绝对值后，依次与第一径定时估计阈值比较，找出互相关子序列中，第一个大于第一径定时估计阈值的互相关值，将该互相关值对应的采样点在采样序列中出现的时刻作为第一径到达时刻，实现符号定时同步。
与现有技术相比，本发明具有如下优点：
第一，由于本发明使用二维时频度量曲面寻找广义频分复用GFDM采样序列的整数倍频偏，并与纠正采样序列的小数倍频偏的结合，纠正了采样序列的大频偏，克服了现有技术Ivan S Gaspar的方法中由于采样序列受到残余整数倍频偏的影响，而导致符号定时同步性能急剧下降的问题，使得本发明具有更准确的符号定时同步的优点。
第二，由于本发明使用纠正采样序列的小数倍频偏，估计小于一个子载波范围的频偏，再使用二维时频度量曲面，估计整数倍频偏，使得本发明的频偏估计范围覆盖整个广义频分复用GFDM系统带宽，克服了现有技术Ivan S Gaspar的方法的频偏估计范围限制在一个子载波带宽内，而导致频谱资源浪费的问题，使得本发明的频偏估计范围远大于一个子载波带宽，节约频谱资源的优点。
附图说明
图1为本发明的流程图；
图2为本发明的仿真图。
具体实施方式
下面结合附图，对本发明的作进一步的详细描述。
参照图1，本发明的实现步骤作进一步的详细描述
步骤1，接收电信号。
广义频分复用GFDM的接收机检测广义频分复用GFDM的发送机发送的模拟电信号。
对检测到的模拟电信号进行模数转换，得到实数字信号。
对实数字信号进行希尔伯特变换，得到复数信号。
对复数字信号进行数字下变频处理，得到采样序列。
步骤2，对采样序列进行粗符号定时同步。
利用自相关公式，计算采样序列中每一个采样点的自相关值，将所有的自相关值组成自相关序列，该自相关序列中有一个与循环前缀等长的定时平台序列。
所述自相关公式为：
Figure GDA0002274081860000061
其中，Pd表示采样序列中第d个采样点的自相关值，N0表示计算每一个采样点的自相关值需要的采样点的总数，该总数的取值由系统参数决定的广义频分复用GFDM前导码序列的子载波数和时隙数确定，∑表示求和操作，k0表示自相关操作中采样点的序号，r(·)表示采样点，T表示共轭操作，m表示采样序列中采样点的序号，该值等于d的大小，*表示相乘操作，K表示由系统参数决定的广义频分复用GFDM的子载波数。
利用能量值公式，计算采样序列中每一个采样点的能量值，将所有的能量值组成能量序列。
所述能量值公式为：
Figure GDA0002274081860000062
其中，Rd表示采样序列中第d个采样点的能量值，N1表示计算每一个采样点的能量值时需要的采样点的总数，该总数的取值由系统参数决定的广义频分复用GFDM前导码序列的子载波数和时隙数确定，|·|表示绝对值操作，k1表示能量值操作中采样点的序号。
将自相关序列中的每一个自相关值依次作为截至自相关值，向前截取与循环前缀等长的子自相关序列，获得多个子自相关序列；其中，所述循环前缀的长度由广义频分复用GFDM系统参数决定。
将截至自相关值对应的采样点的序号作为子自相关序列的编号。
将能量序列中的每一个能量值依次作为截至能量值，向前截取与循环前缀等长的子能量序列，获得多个子能量序列。
将截至能量值对应的采样点的序号作为子能量序列的编号。
将具有相同编号的子自相关序列和子能量序列进行相除操作，并将相除操作后的结果进行取绝对值操作，获得多个归一化子自相关序列。
对每一个归一化子自相关序列进行相加操作，得到对应采样点的粗符号定时度量值，将所有粗符号定时度量值组成消除定时平台序列的粗符号定时度量序列。
找出粗符号定时度量序列中最大值对应的采样点，该采样点在采样序列中出现的时刻为粗符号定时同步时刻。
步骤3，纠正采样序列的小数倍频偏。
找出粗符号定时同步时刻对应的采样点的自相关值。
对粗符号定时同步时刻对应的采样点的自相关值进行取相位操作，获得该自相关值的相位，并将该自相关值的相位与圆周率进行相除操作，获得采样序列的小数倍频偏估计值。其中，所述小数倍频偏为经过广义频分复用GFDM子载波带宽归一化处理后的小数倍频偏。
利用小数倍频偏纠正公式，纠正采样序列的小数倍频偏，获得无小数倍频偏的采样序列。
所述小数倍频偏纠正公式为：
Figure GDA0002274081860000071
其中，rc(·)表示纠正小数倍频偏后的采样点，e表示自然底数，
Figure GDA0002274081860000072
表示小数倍频偏估计值。
步骤4，选择路径候选定时时刻。
对本地前导码序列进行共轭操作，得到共轭前导码序列。其中，所述本地前导码序列是具有两段重复结构的序列。
将无小数倍频偏的采样序列中的每一个采样点依次作为起始点，向后截取与共轭前导码序列等长的子采样序列，将每一个子采样序列与共轭前导码序列进行相乘操作，获得多个子序列。
利用差分互相关公式，计算采样序列中每一个采样点的差分互相关值，将所有差分互相关值组成差分互相关序列。
所述差分互相关公式为：
Figure GDA0002274081860000081
其中，Qd表示在第d个纠正小数倍频偏的采样点的差分互相关值，Ud(·)表示第d个纠正小数倍频偏的采样点对应的子序列的元素，k3表示子序列中元素的序号。
对差分互相关序列进行取绝对值操作，将取绝对值操作后的差分互相关序列中的每一个差分互相关值进行平方操作，获得对应的路径候选定时度量值，将所有的路径候选定时度量值组成路径候选定时度量序列。
将路径候选定时度量序列，按照从大到小排列，找出前64个路径候选定时度量值对应的64个采样点，将64个采样点在采样序列中出现的时刻作为路径候选定时时刻。
步骤5，绘制二维时频度量曲面。
第1步，将64个路径候选定时时刻依次送入二维时频估计器。
第2步，二维时频估计器根据输入的路径候选定时时刻，找到该时刻的采样点，再找出该采样点对应的子序列。
第3步，对子序列进行快速傅里叶变换。
第4步，对快速傅里叶变换后的结果进行取绝对值操作，得到二维时频度量子序列。
第5步，判断64个路径候选定时时刻是否全部送入二维时频估计器，若是，则执行本步骤的第6步，否则，执行本步骤的第2步。
第6步，64个路径候选定时时刻全部送入二维时频估计器后，得到64个路径候选定时时刻对应的64个二维时频度量子序列，绘制由64个二维时频度量子序列构成的二维度时频度量曲面。
步骤6，估计路径定时时刻。
找出二维度时频度量曲面的最大值，将最大值所在的二维时频度量子序列对应的路径候选定时时刻，作为路径到达时刻。
步骤7，纠正采样序列的整数倍频偏。
找出二维度时频度量曲面的最大值对应的快速傅里叶变换的频率点值，将该频率点值作为采样序列的整数倍频偏估计值。其中，所述整数倍频偏为经过广义频分复用GFDM子载波带宽归一化处理后的整数倍频偏。
利用整数倍频偏纠正公式，纠正采样序列的整数倍频偏，得到纠正大频偏后的无频偏采样序列，实现载波频率同步。
所述整数倍频偏纠正公式为：
Figure GDA0002274081860000091
其中，ri(·)表示无频偏采样序列的采样点，
Figure GDA0002274081860000092
表示整数倍频偏估计值。
步骤8，估计第一径到达时刻。
将无频偏采样序列中的每一个采样点依次作为起始点，向后截取与共轭前导码序列等长的无频偏子采样序列，获得多个无频偏子采样序列。
对每一个采样点对应的无频偏子采样序列与共轭前导码序列进行相乘操作，将相乘后的结果相加，得到一个互相关值。
将所有采样点对应的互相关值组成互相关序列。
将路径到达时刻对应的采样点所对应的互相值，作为截至互相关值。
在互相关序列中，从截至互相关值开始，向前截取与循环前缀等长的互相关子序列。
利用第一径定时估计阈值公式，计算第一径定时估计阈值。
所述第一径定时估计阈值公式为：
Figure GDA0002274081860000093
其中，TTh表示第一径定时估计阈值，ln(·)表示取自然对数操作，PFA表示错误预警概率，该值由系统参数决定的广义频分复用GFDM系统参数决定，ρ表示衰落信道多径参数，该衰落信道多径参数的取值范围在信道长度到循环前缀长度范围之间，Px(·)互相关子序列中的互相关值，
Figure GDA0002274081860000094
表示路径定时时刻对应的采样序列中的采样点的序号，k4表示互相关子序列中互相关值的序号。
对互相关子序列的每一个互相关值取绝对值后，依次与第一径定时估计阈值比较，找出互相关子序列中，第一个大于第一径定时估计阈值的互相关值，将该互相关值对应的采样点在采样序列中出现的时刻作为第一径到达时刻，实现符号定时同步。
本发明的效果可以通过下面的仿真得到进一步证明。
1.仿真条件：
本发明仿真实验采用MATLAB仿真软件实现，仿真实验1和仿真实验2的条件相同，仿真实验3的条件与仿真实验1和仿真实验2的条件不同；
仿真实验1和仿真实验2的条件是：广义频分复用GFDM子载波数为128，时隙数为2，循环前缀长度为32，错误检测概率为10-6，广义频分复用GFDM前导码序列的成型滤波器为矩形滤波器，频偏为2.2，信道环境为瑞利衰落信道，每一径的信道抽头为0.65，0，0，0，0.43，0，0，0，0.2，每一径抽头的瑞利随机变量服从均值为0，方差为
Figure GDA0002274081860000101
的瑞利分布。
仿真实验3的条件是：广义频分复用GFDM子载波数为128，时隙数为2，循环前缀长度为32，错误检测概率为10-6，广义频分复用GFDM前导码序列的成型滤波器为矩形滤波器，频偏范围为-10:0.5:10，无信道环境的影响。
二、仿真内容与结果分析：
下面结合附图2的仿真图，对本发明的效果作进一步的描述。
仿真实验1：
利用广义频分复用GFDM系统，对本发明的方法和现有的Ivan S Gaspar的方法的频偏估计均方误差性能进行仿真，仿真实验的结果如图2(a)所示。
图2(a)的横坐标表示广义频分复用GFDM系统的信噪比，单位dB，纵坐标表示频偏估计均方误差。图2(a)中以圆圈标示的曲线，表示采用本发明的方法得到的频偏估计均方误差与信噪比之间关系的曲线。图2(a)中以正方形标示的曲线，表示采用Ivan S Gaspar的方法得到的频偏估计均方误差与信噪比之间关系的曲线。
从图2(a)中可以看出：在信噪比为0dB时，采用本发明的方法得到的频偏估计均方误差接近10-3，而Ivan S Gaspar的方法得到的频偏估计均方误差大于100，而且，随着信噪比的增加，本发明的方法得到的频偏估计均方误差不断减小，而Ivan S Gaspar的方法得到的频偏估计均方误差几乎保持在100以上。
仿真实验2：
利用广义频分复用GFDM系统，对本发明的方法和现有的Ivan S Gaspar的方法的时偏估计均方误差性能进行仿真，仿真结果如图2(b)所示。
图2(b)的横坐标表示广义频分复用GFDM系统的信噪比，单位dB，纵坐标表示时偏估计均方误差。图2(b)中以圆圈标示的曲线，表示采用本发明的方法得到的时偏估计均方误差与信噪比之间关系的曲线。图2(b)中以正方形标示的曲线，表示采用Ivan S Gaspar的方法得到的时偏估计均方误差与信噪比之间关系的曲线。
图2(b)中可以看出：在信噪比为0dB时，采用本发明的方法得到的时偏估计均方误差小于100，而Ivan S Gaspar的方法得到的时偏估计均方误差大于102，而且，随着信噪比的增加，本发明的方法得到的时偏估计均方误差不断减小，而Ivan S Gaspar的方法得到的频偏估计均方误差几乎保持在102以上。
仿真实验3：
对本发明的方法和现有的Ivan S Gaspar的方法的频偏估计范围进行仿真，仿真结果如图2(c)所示。
图2(c)的横坐标表示实际频偏值，纵坐标表示频偏估计值。图2(c)中以三角形标示的曲线，表示采用本发明的方法的频偏估计值与实际频偏值关系的曲线。图2(c)中以正方形标示的曲线，表示Ivan S Gaspar的方法的频偏估计值与实际频偏值关系的曲线。
图2(c)中可以看出：当实际频偏值在-1:0.5:1范围内时，本发明的方法和Ivan SGaspar的方法均能得到正确的估计频偏值；当实际频偏值在-10:0.5:-1和1：0.5：10范围内时，本发明的方法能够得到正确的估计频偏值，而Ivan S Gaspar的方法的得到估计频偏值在-1:0.5:1范围内变化，无法得到正确的估计频偏值。
综上所述，采用本发明的一种用于纠正大频偏的广义频分复用GFDM时频同步方法，能很好的消除广义频分复用GFDM的整数倍频偏，消除了大频偏对广义频分复用GFDM同步性能的影响，使得第一径定时时刻估计更准确，符号定时同步性能更好，而且频偏估计范围远大于一个子载波带宽。

Claims (10)

1.一种用于纠正大频偏的广义频分复用时频同步方法，其特征在于，将接收的广义频分复用GFDM模拟电信号处理后的采样序列，依次经粗符号定时同步、小数倍频偏纠正、整数倍频偏纠正，得到纠正大频偏后的无频偏采样序列，实现载波频率同步，从路径到达时刻向前搜索第一径到达时刻，实现符号定时同步，该方法的具体步骤包括如下：
(1)接收电信号：
(1a)广义频分复用GFDM的接收机检测广义频分复用GFDM的发送机发送的模拟电信号；
(1b)对检测到的模拟电信号进行模数转换，得到实数字信号；
(1c)对实数字信号进行希尔伯特变换，得到复数信号；
(1d)对复数字信号进行数字下变频处理，得到采样序列；
(2)对采样序列进行粗符号定时同步：
(2a)利用自相关公式，计算采样序列中每一个采样点的自相关值，将所有的自相关值组成自相关序列；
(2b)利用能量值公式，计算采样序列中每一个采样点的能量值，将所有的能量值组成能量序列；
(2c)将自相关序列中的每一个自相关值依次作为截至自相关值，向前截取与循环前缀等长的子自相关序列，获得多个子自相关序列；其中，所述循环前缀的长度由广义频分复用GFDM系统参数决定；
(2d)将截至自相关值对应的采样点的序号作为子自相关序列的编号；
(2e)将能量序列中的每一个能量值依次作为截至能量值，向前截取与循环前缀等长的子能量序列，获得多个子能量序列；
(2f)将截至能量值对应的采样点的序号作为子能量序列的编号；
(2g)将具有相同编号的子自相关序列和子能量序列进行相除操作，并将相除操作后的结果进行取绝对值操作，获得多个归一化子自相关序列；
(2h)对每一个归一化子自相关序列进行相加操作，得到对应采样点的粗符号定时度量值，将所有粗符号定时度量值组成粗符号定时度量序列；
(2i)找出粗符号定时度量序列中最大值对应的采样点，该采样点在采样序列中出现的时刻为粗符号定时同步时刻；
(3)纠正采样序列的小数倍频偏：
(3a)找出粗符号定时同步时刻对应的采样点的自相关值；
(3b)对粗符号定时同步时刻对应的采样点的自相关值进行取相位操作，获得该自相关值的相位，并将该自相关值的相位与圆周率进行相除操作，获得采样序列的小数倍频偏估计值；
(3c)利用小数倍频偏纠正公式，纠正采样序列的小数倍频偏，获得无小数倍频偏的采样序列；
(4)选择路径候选定时时刻：
(4a)对本地前导码序列进行共轭操作，得到共轭前导码序列；
(4b)将无小数倍频偏的采样序列中的每一个采样点依次作为起始点，向后截取与共轭前导码序列等长的子采样序列，将每一个子采样序列与共轭前导码序列进行相乘操作，获得多个子序列；
(4c)利用差分互相关公式，计算采样序列中每一个采样点的差分互相关值，将所有差分互相关值组成差分互相关序列；
(4d)对差分互相关序列进行取绝对值操作，将取绝对值操作后的差分互相关序列中的每一个差分互相关值进行平方操作，获得对应的路径候选定时度量值，将所有的路径候选定时度量值组成路径候选定时度量序列；
(4e)将路径候选定时度量序列，按照从大到小排列，找出前64个路径候选定时度量值对应的64个采样点，将64个采样点在采样序列中出现的时刻作为路径候选定时时刻；
(5)绘制二维时频度量曲面：
(5a)将64个路径候选定时时刻依次送入二维时频估计器；
(5b)二维时频估计器根据输入的路径候选定时时刻，找到该时刻的采样点，再找出该采样点对应的子序列；
(5c)对子序列进行快速傅里叶变换；
(5d)对快速傅里叶变换后的结果进行取绝对值操作，得到二维时频度量子序列；
(5e)判断64个路径候选定时时刻是否全部送入二维时频估计器，若是，则执行步骤(5f)，否则，执行步骤(5b)；
(5f)64个路径候选定时时刻全部送入二维时频估计器后，得到64个路径候选定时时刻对应的64个二维时频度量子序列，绘制由64个二维时频度量子序列构成的二维度时频度量曲面；
(6)估计路径定时时刻：
找出二维度时频度量曲面的最大值，将最大值所在的二维时频度量子序列对应的路径候选定时时刻，作为路径到达时刻；
(7)纠正采样序列的整数倍频偏：
(7a)找出二维度时频度量曲面的最大值对应的快速傅里叶变换的频率点值，将该频率点值作为采样序列的整数倍频偏估计值；
(7b)利用整数倍频偏纠正公式，纠正无小数倍频偏的采样序列的整数倍频偏，得到纠正大频偏后的无频偏采样序列，实现载波频率同步；
(8)估计第一径到达时刻：
(8a)将无频偏采样序列中的每一个采样点依次作为起始点，向后截取与共轭前导码序列等长的无频偏子采样序列，获得多个无频偏子采样序列；
(8b)对每一个采样点对应的无频偏子采样序列与共轭前导码序列进行相乘操作，将相乘后的结果相加，得到一个互相关值；
(8c)将所有采样点对应的互相关值组成互相关序列；
(8d)将路径到达时刻对应的采样点所对应的互相值，作为截至互相关值；
(8e)在互相关序列中，从截至互相关值开始，向前截取与循环前缀等长的互相关子序列；
(8f)利用第一径定时估计阈值公式，计算第一径定时估计阈值；
(8g)对互相关子序列的每一个互相关值取绝对值后，依次与第一径定时估计阈值比较，找出互相关子序列中，第一个大于第一径定时估计阈值的互相关值，将该互相关值对应的采样点在采样序列中出现的时刻作为第一径到达时刻，实现符号定时同步。
2.根据权利要求1所述一种用于纠正大频偏的广义频分复用时频同步方法，其特征在于，步骤(2a)中所述自相关公式为：
Figure FDA0002274081850000041
其中，Pd表示采样序列中第d个采样点的自相关值，N0表示计算每一个采样点的自相关值需要的采样点的总数，该总数的取值由系统参数决定的广义频分复用GFDM前导码序列的子载波数和时隙数确定，Σ表示求和操作，k0表示自相关操作中采样点的序号，r(·)表示采样点，T表示共轭操作，m表示采样序列中采样点的序号，该值等于d的大小，*表示相乘操作，K表示由系统参数决定的广义频分复用GFDM的子载波数。
3.根据权利要求2所述一种用于纠正大频偏的广义频分复用时频同步方法，其特征在于，步骤(2b)中所述能量值公式为：
Figure FDA0002274081850000042
其中，Rd表示采样序列中第d个采样点的能量值，N1表示计算每一个采样点的能量值时需要的采样点的总数，该总数的取值由系统参数决定的广义频分复用GFDM前导码序列的子载波数和时隙数确定，|·|表示绝对值操作，k1表示能量值操作中采样点的序号。
4.根据权利要求1所述一种用于纠正大频偏的广义频分复用时频同步方法，其特征在于，步骤(3)中所述小数倍频偏为经过广义频分复用GFDM子载波带宽归一化处理后的小数倍频偏。
5.根据权利要求2所述一种用于纠正大频偏的广义频分复用时频同步方法，其特征在于，步骤(3c)中所述小数倍频偏纠正公式为：
Figure FDA0002274081850000043
其中，rc(·)表示纠正小数倍频偏后的采样点，e表示自然底数，
Figure FDA0002274081850000055
表示小数倍频偏估计值。
6.根据权利要求1所述一种用于纠正大频偏的广义频分复用时频同步方法，其特征在于，步骤(4a)中所述本地前导码序列是具有两段重复结构的序列。
7.根据权利要求1所述一种用于纠正大频偏的广义频分复用时频同步方法，其特征在于，步骤(4c)中所述差分互相关公式为：
Figure FDA0002274081850000051
其中，Qd表示在第d个纠正小数倍频偏的采样点的差分互相关值，Ud(·)表示第d个纠正小数倍频偏的采样点对应的子序列的元素，k3表示子序列中元素的序号。
8.根据权利要求1所述一种用于纠正大频偏的广义频分复用时频同步方法，其特征在于，步骤(7)中所述整数倍频偏为经过广义频分复用GFDM子载波带宽归一化处理后的整数倍频偏。
9.根据权利要求5所述一种用于纠正大频偏的广义频分复用时频同步方法，其特征在于，步骤(7b)中所述整数倍频偏纠正公式为：
Figure FDA0002274081850000052
其中，ri(·)表示无频偏采样序列的采样点，
Figure FDA0002274081850000053
表示整数倍频偏估计值。
10.根据权利要求1所述一种用于纠正大频偏的广义频分复用时频同步方法，其特征在于，步骤(8f)中所述第一径定时估计阈值公式为：
Figure FDA0002274081850000054
其中，TTh表示第一径定时估计阈值，ln(·)表示取自然对数操作，PFA表示错误预警概率，该值由系统参数决定的广义频分复用GFDM系统参数决定，ρ表示衰落信道多径参数，该衰落信道多径参数的取值范围在信道长度到循环前缀长度范围之间，Px(·)互相关子序列中的互相关值，
Figure FDA0002274081850000061
表示路径定时时刻对应的采样序列中的采样点的序号，k4表示互相关子序列中互相关值的序号。
CN201810117264.3A 2018-02-06 2018-02-06 一种用于纠正大频偏的广义频分复用时频同步方法 Active CN108366032B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810117264.3A CN108366032B (zh) 2018-02-06 2018-02-06 一种用于纠正大频偏的广义频分复用时频同步方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810117264.3A CN108366032B (zh) 2018-02-06 2018-02-06 一种用于纠正大频偏的广义频分复用时频同步方法

Publications (2)

Publication Number Publication Date
CN108366032A CN108366032A (zh) 2018-08-03
CN108366032B true CN108366032B (zh) 2020-04-07

Family

ID=63004670

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810117264.3A Active CN108366032B (zh) 2018-02-06 2018-02-06 一种用于纠正大频偏的广义频分复用时频同步方法

Country Status (1)

Country Link
CN (1) CN108366032B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112910815B (zh) * 2021-01-15 2022-03-04 西安电子科技大学 基于正负极性前导码的广义频分复用系统时频同步方法
CN112910817B (zh) * 2021-01-18 2022-01-11 西安电子科技大学 基于一维搜索和能量检测的ofdm时间同步方法
CN113438730B (zh) * 2021-06-08 2022-05-17 武汉大学 一种基于gfdm信号的无线定位方法
CN114070685B (zh) * 2021-11-01 2023-12-08 泰凌微电子(上海)股份有限公司 频偏估计方法和装置
CN114928520B (zh) * 2022-04-22 2023-12-12 南京邮电大学 一种用于补偿频偏的广义频分复用时频同步方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102694763A (zh) * 2012-05-31 2012-09-26 重庆邮电大学 一种用于td-lte系统整数倍频偏估计的方法
CN103220252A (zh) * 2013-04-10 2013-07-24 安徽华东光电技术研究所 编码正交频分复用的无线信号接收处理装置及其处理方法
CN103379082A (zh) * 2012-04-25 2013-10-30 马维尔国际有限公司 Lte通信系统中的时频同步方法和装置
CN103825859A (zh) * 2014-03-10 2014-05-28 江苏物联网研究发展中心 一种正交频分复用信号的同步捕获方法及接收端设备
CN106549896A (zh) * 2015-09-23 2017-03-29 中南大学 一种通用零自相关码的新型定时同步与频偏估计算法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9742611B2 (en) * 2013-10-23 2017-08-22 Silicon Laboratories Inc. Synchronizing orthogonal frequency division multiplexed (OFDM) symbols in a receiver

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103379082A (zh) * 2012-04-25 2013-10-30 马维尔国际有限公司 Lte通信系统中的时频同步方法和装置
CN102694763A (zh) * 2012-05-31 2012-09-26 重庆邮电大学 一种用于td-lte系统整数倍频偏估计的方法
CN103220252A (zh) * 2013-04-10 2013-07-24 安徽华东光电技术研究所 编码正交频分复用的无线信号接收处理装置及其处理方法
CN103825859A (zh) * 2014-03-10 2014-05-28 江苏物联网研究发展中心 一种正交频分复用信号的同步捕获方法及接收端设备
CN106549896A (zh) * 2015-09-23 2017-03-29 中南大学 一种通用零自相关码的新型定时同步与频偏估计算法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Effect of Timing and FrequErrors on GFDM Systemsency Synchronization;Jae-Hyung Choi 等;《2015 International Conference on Information and Communication Technology Convergence (ICTC)》;20151217;1322-1325 *
一种改进的GFDM时频同步算法;吴虹 等;《电讯技术》;20161228;第56卷(第12期);1322-1326 *

Also Published As

Publication number Publication date
CN108366032A (zh) 2018-08-03

Similar Documents

Publication Publication Date Title
CN108366032B (zh) 一种用于纠正大频偏的广义频分复用时频同步方法
CN110224968B (zh) 一种ofdm通信系统中的帧定时同步方法和装置
KR100947794B1 (ko) 미세 타이밍 획득
KR100808949B1 (ko) Ｌｃｉｃ-ｄｆｅ를 이용한 채널 추정 방법 및 그 장치
US7627059B2 (en) Method of robust timing detection and carrier frequency offset estimation for OFDM systems
US8300713B2 (en) Preamble sequence detection and integral carrier frequency offset estimation method for OFDM/OFDMA wireless communication system
KR100579531B1 (ko) Ｏｆｄｍ 수신기에 적용되는 심볼시간 동기 장치 및 그 방법
US20070047433A1 (en) Method and circuit for fine timing synchronization in the orthogonal frequency division multiplexing baseband receiver for IEEE 802.11a/g wireless LAN standard
CN101836411A (zh) 使用时分多路复用导频在广播ofdm系统中的同步
CN113518052B (zh) 一种正交频分复用通信的鲁棒频率偏移估计方法及装置
US20070019538A1 (en) Symbol Synchronization for OFDM Systems
CN104836770B (zh) 一种基于相关平均与加窗的定时估计方法
WO2007112681A1 (fr) Procédé d'estimation de décalage de fréquence multiple entier ofdm
CN113438730B (zh) 一种基于gfdm信号的无线定位方法
KR100664018B1 (ko) 직교주파수 분할다중화 수신기의 동기 검출장치
CN101836410A (zh) 使用时分多路复用导频在广播ofdm系统中的同步
CN101039300A (zh) 一种同步定时的方法
CN110943949B (zh) 一种基于sc-fde传输体制的宽带角跟踪方法
JP5657733B2 (ja) 無線通信システムにおけるタイミング情報の獲得
EP2245814B1 (en) Frame timing and carrier frequency recovery for frequency selective signals
KR20110135662A (ko) 다중경로 페이딩 채널에서 타이밍 오차 추정이 가능한 ｏｆｄｍ 수신기, 이를 포함하는 ｏｆｄｍ 시스템 및 이들의 타이밍 오차 추정방법
KR101100460B1 (ko) Ｏｆｄｍ 시스템에서의 심벌 및 프레임의 시간 동기 방법및 장치
KR20050003663A (ko) Ｏｆｄｍ 수신기의 심벌 타이밍 복원장치 및 방법
CN112910817B (zh) 基于一维搜索和能量检测的ofdm时间同步方法
CN112134823B (zh) 基于cir峰值的ofdm信号符号同步误差估计方法

Legal Events

Date Code Title Description
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