CN104125188B - 一种基于Zadoff‑Chu序列的OFDM频率同步方法 - Google Patents
一种基于Zadoff‑Chu序列的OFDM频率同步方法 Download PDFInfo
- Publication number
- CN104125188B CN104125188B CN201410394369.5A CN201410394369A CN104125188B CN 104125188 B CN104125188 B CN 104125188B CN 201410394369 A CN201410394369 A CN 201410394369A CN 104125188 B CN104125188 B CN 104125188B
- Authority
- CN
- China
- Prior art keywords
- zadoff
- frequency
- chu sequence
- ofdm
- reception signal
- 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
Links
Landscapes
- Synchronisation In Digital Transmission Systems (AREA)
Abstract
本发明公开了一种基于Zadoff‑Chu序列的OFDM频率同步方法，属于无线通信技术领域，针对传统频率同步算法的精度受限于定时同步准确度且实现复杂的缺点，利用两段重复序列估计小数频偏，根据频偏对Zadoff‑Chu序列自相关的影响估计整数频偏，估计范围可达整个OFDM带宽，且所有操作均在时域完成，无需时频变换，与传统的频偏估计算法相比，运算复杂度较低，鲁棒性好。通过验证表明，该方法在高速移动多径环境下具有较好的性能。
Description
技术领域
本发明涉及到OFDM***频率同步技术，特别是涉及一种基于Zadoff-Chu序列的OFDM频率同步方法，属于通信技术领域。
背景技术
正交频分复用(Orthogonal Frequency Division Multiplexing,OFDM)***可提供更高的数据传输速率和频谱利用率，抗多径衰落能力强，是新一代宽带移动通信***的核心技术。载波同步是OFDM的关键技术，一般情况下，高斯信道下频率偏移误差要小于子载波间隔的4％，多径信道环境下频率偏移误差要小于子载波间隔的2％。因此，频偏估计算法的研究具有重要的现实意义。
现有的OFDM频率同步方法有：文献[1]提出了一种S&C算法，可以估计出整数和小数倍频偏，但当子载波数小于***IFFT点数时，载波同步性能会迅速下降；文献[2]提出的利用L个重复序列的最佳线性无偏估计(Best Linear Unbiased Estimation,BLUE)算法，估计范围可达L/2个子载波间隔，但为保证算法性能，每段重复序列的长度必须大于信道最大时延扩展，降低了传输效率，相同条件下其估计比S&C算法精确，但复杂度大大提高；文献[3]虽然得到了载波频率偏移的闭合表达式，但算法中指数因子p的设置受限于估计精度、整数倍频偏估计范围和***可容忍的峰均比(Peak-Average Power Ratio,PAPR)，需要通过仿真反复验证，不具有通用性；文献[4]使用恒包络零自相关(Constant EnvelopeZeroAuto-Correlation,CAZAC)序列进行频率同步，算法估计精度高，但其从频偏估计到频偏补偿的反馈时间较长，不适用于突发通信；文献[5]提出了一种Ren算法，通过对CAZAC序列加权进行频率同步，具有较大的频偏估计范围，但该算法在独立瑞利衰落信道下性能较差。
上述参考文献分别为：
[1]Schmidl T M,Cox D C.Robust frequency and timing synchronizationfor OFDM[J].IEEE Transactions on Communications,1997,45(12):1613-1621.
[2]Morelli M,Mengali U.An improved frequency offset estimator forOFDM applications[J].IEEE Communications Letters,1999,3(3):75-77.
[3]Laourine A,Stephenne A,Affes S.A new OFDM synchronization symbolfor carrier frequency offset estimation[J].IEEE Signal Processing Letters,2007,14(5):321-324.
[4]严春林,李少谦,唐友喜,等.利用CAZAC序列的OFDM频率同步方法[J].电子与信息学报,2006,28(1),139-142.
[5]Ren G,Chang Y,Zhang H,et al.Synchronization method based on a newconstant envelop preamble for OFDM systems[J].IEEE Transactions onBroadcasting,2005,51(1):139-143.
发明内容
针对上述缺陷，本发明提出了一种基于Zadoff-Chu序列的OFDM频率同步方法，该方法充分利用Zadoff-Chu序列良好的自相关特性，无需进行时频变换，在保证算法估计精度的同时，其频偏估计范围可达N个子载波间隔，可实现OFDM***精确频率同步。
为了达到上述目的，本发明所采用的具体技术方案如下：
一种基于Zadoff-Chu序列的OFDM频率同步方法，其关键在于按照以下步骤进行：
步骤S1：采用两组结构相同的训练符号TS1和TS2组成训练序列，其中TS1按照“Cn,Cn”构成，TS2按照“Cn *,Cn *”构成，其中Cn为本地Zadoff-Chu序列，周期为Np，且Np＝N/2，N为***子载波个数，Cn *为Cn的复共轭；
步骤S2：利用TS1的接收信号求自相关得到小数频偏并进行补偿；
步骤S3：利用TS1的接收信号循环前缀CP和TS1的最后Ng个样点估计残留频偏并进行补偿，Ng为循环前缀CP的长度；
步骤S4：利用TS1的接收信号定时度量函数M1(d)和TS2的接收信号定时度量函数M2(d)，并按照计算整数倍频偏并进行补偿，其中表示M1(d)中最大值所对应的位置，表示M2(d)中最大值所对应的位置。
作为一种优选方案，按照计算TS1的接收信号定时度量函数，按照计算TS2的接收信号定时度量函数，其中r(d+n)为第d+n个基带等效接收信号的采样值，C(n)为序列Cn中的第n个采样点的值，C*(n)为序列Cn*中的第n个采样点的值。
在实施过程中，所述本地Zadoff-Chu序列可以按照生成，其中j为虚数单位，r为正整数，且r与Np互质，在这里r设为1。
本发明的显著效果是：本发明提出了一种利用Zadoff-Chu序列的OFDM频率同步方法，该方法频率估计范围可达整个OFDM带宽，且所有操作均在时域完成，无需时频变换，与传统的频偏估计算法相比，运算复杂度较低，鲁棒性好。通过验证表明，该方法在存在大频偏的多径信道环境中具有良好的性能。
附图说明
图1为频偏对Zadoff-Chu序列相关性的影响曲线；
图2为本地训练符号结构示意图；
图3为SUI3信道下，各种算法的MSE结果对比分析图；
图4为ITU-VB信道下，各种算法的MSE结果对比分析图。
具体实施方式
下面结合附图和实施例对本发明的工作原理和显著效果作进一步说明：
在OFDM***的复基带等效模型中，发送端OFDM符号的复基带数据的采样为：
其中，N是载波总数，X(k)表示调制在第k个子载波上的数据。加入循环前缀CP后，一个OFDM符号可以表示为其中Ng是循环前缀CP的长度。
假设收发两端仅存在载波频偏，则接收信号可以表示为：
r(n)＝x(n)ej2πnε/N+w(n) (2)
其中ε表示被子载波间隔归一化后的频率偏移，w(n)是信道中的加性高斯白噪声。
如果采用Zadoff-Chu序列作为训练序列，周期为Np(Np为偶数)的Zadoff-Chu序列定义为：
上式j为虚数单位，r为正整数，r与Np互质。若Np＝N/2，r＝1，其自相关函数为：
上式中mod表示取模。由式(4)可知Zadoff-Chu序列的循环自相关具有旁瓣为0的特性。此外，Zadoff-Chu序列还有低峰均比，恒包络的特点，而且其FFT变换还是一个Zadoff-Chu序列，这些特性使得Zadoff-Chu序列广泛应用于OFDM***中。
实际通信***中存在载波频偏，且移动情况下会带来多普勒频移。由式(2)(4)可知，当存在载波频偏时，接收信号可以表示为：
由上式可知，对于生成式为(3)式的Zadoff-Chu序列而言，载波频率偏移会导致接收信号在时域循环移位。当ε为偶数时，上式忽略噪声的影响，接收信号相对于原始序列移位了ε/2个样点，图1给出了频偏对Zadoff-Chu序列相关峰值影响的示意图。
因此，本实施例中，按照ε＝εi+εf+εr将OFDM***中频率偏移分为三部分，其中，εi为整数倍频偏，εf为小数倍频偏，εr为残留频偏。其频率同步的具体过程如下：
步骤S1：采用两组结构相同的训练符号TS1和TS2组成训练序列，其中TS1按照“Cn,Cn”构成，TS2按照“Cn *,Cn *”构成，其中Cn为本地Zadoff-Chu序列，按照生成，其中j为虚数单位，r＝1，周期为Np，且Np＝N/2，N为***子载波个数，Cn *为Cn的复共轭，具体结构如图2所示。
步骤S2：利用TS1的接收信号求自相关得到小数频偏并进行补偿，通常按照公式(6)计算，这里的频偏估计范围为±1：
步骤S3：补偿小数频偏εf之后，可以利用TS1的接收信号循环前缀CP和TS1的最后Ng个样点估计残留频偏并进行补偿，本例中采用公式(7)计算；
步骤S4：先按照计算TS1的接收信号定时度量函数M1(d)，再按照计算TS2的接收信号定时度量函数M2(d)，其中，r(d+n)为第d+n个基带等效接收信号的采样值，C(n)为序列Cn中的第n个采样点的值，C*(n)为序列Cn *中的第n个采样点的值。
根据(5)式的分析可知，频率偏移会使M1(d)和M2(d)的峰值位置产生移位。由于TS1和TS2互为复共轭，所以M1(d)和M2(d)的峰值移位方向相反，因此可以按照计算整数倍频偏并进行补偿，其中表示M1(d)中最大值所对应的位置，表示M2(d)中最大值所对应的位置。
为验证本方法的有效性，本实施例选择SUI-3准静态多径信道和ITU-M.1225Vehicle Channel B信道进行仿真。比较了文献[5]所提算法与本方法在不同多普勒频移下均方误差(MSE)性能。仿真参数为：FFT点数N＝1024，子载波间隔Δf＝97.66kHz，***采样速率fs＝100MHz，循环前缀长度Ncp＝N4。两种信道环境下载波频偏均设置为34.2547个子载波间隔。信道参数见表1。
表1信道参数
图3和图4分别比较了本方法与文献[5]所提算法在SUI-3静态多径信道环境和ITU-VB移动信道环境中的检测性能曲线。从图中可以看到，本方法的检测性能明显优于Ren算法。静态和低速移动环境中，由于信道环境几乎不变，本方法的估计方差近似随信噪比增加而线性下降，没有地板效应。高速移动环境下，由于受到时变信道的影响，算法的估计方差较相同信噪比条件下低速移动环境中的性能稍差，且随信噪比增加估计方差的下降速度减缓，当信噪比大于一定程度时，估计方差不再随信噪比的增加而下降，但仍优于Ren算法性能。
Ren算法利用接收到的前导r(k)与加入不同频偏的本地前导s′(k)进行相关，将相关函数取得最大值时本地前导所加频偏作为整数频偏的估计值。当信道环境不断变化(存在多普勒频移)且各径能量相近时，Ren所提的定时同步算法会存在偏差，r(k)与s′(k)间的互相关也必然存在偏差，导致频偏估计性能较差。本方法的频偏估计是根据两个训练符号与本地符号相关峰值之间的相对距离得到，定时偏移对这两个峰值影响相同，因此对定时偏差不敏感，在两种信道环境中性能均较好。
最后需要说明的是，以上详细描述仅仅为本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。
Claims (3)
1.一种基于Zadoff-Chu序列的OFDM频率同步方法，其特征在于按照以下步骤进行：
步骤S1：采用两组结构相同的训练符号TS1和TS2组成训练序列，其中TS1按照“Cn,Cn”构成，TS2按照“Cn *,Cn *”构成，其中Cn为本地Zadoff-Chu序列，周期为Np，且Np＝N/2，N为***子载波个数，Cn *为Cn的复共轭；
步骤S2：利用TS1的接收信号求自相关得到小数频偏并进行补偿；
步骤S3：利用TS1的接收信号循环前缀CP和TS1的最后Ng个样点估计残留频偏并进行补偿，Ng为循环前缀CP的长度；
步骤S4：利用TS1的接收信号定时度量函数M1(d)和TS2的接收信号定时度量函数M2(d)，并按照计算整数倍频偏并进行补偿，其中表示M1(d)中最大值所对应的位置，表示M2(d)中最大值所对应的位置；
具体按照计算TS1的接收信号定时度量函数，按照计算TS2的接收信号定时度量函数，其中r(d+n)为第d+n个基带等效接收信号的采样值，C(n)为序列Cn中的第n个采样点的值，C*(n)为序列Cn *中的第n个采样点的值。
2.根据权利要求1所述的一种基于Zadoff-Chu序列的OFDM频率同步方法，其特征在于：所述本地Zadoff-Chu序列按照生成，其中j为虚数单位，r为正整数，且r与Np互质。
3.根据权利要求2所述的一种基于Zadoff-Chu序列的OFDM频率同步方法，其特征在于：本地Zadoff-Chu序列中r设为1。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410394369.5A CN104125188B (zh) | 2014-08-12 | 2014-08-12 | 一种基于Zadoff‑Chu序列的OFDM频率同步方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410394369.5A CN104125188B (zh) | 2014-08-12 | 2014-08-12 | 一种基于Zadoff‑Chu序列的OFDM频率同步方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104125188A CN104125188A (zh) | 2014-10-29 |
| CN104125188B true CN104125188B (zh) | 2017-03-22 |
Family
ID=51770455
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410394369.5A Active CN104125188B (zh) | 2014-08-12 | 2014-08-12 | 一种基于Zadoff‑Chu序列的OFDM频率同步方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104125188B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107431679B (zh) * | 2015-09-24 | 2020-09-25 | 华为技术有限公司 | 同步信号传输方法及装置 |
| CN107204945B (zh) * | 2016-03-18 | 2020-11-03 | 富士通株式会社 | 频差估计装置、信道间隔估计装置、方法和*** |
| CN106101043B (zh) * | 2016-05-31 | 2019-05-17 | 中国航天科技集团公司第九研究院第七七一研究所 | 宽带无线通信***中时频联合估计方法 |
| CN107835142B (zh) * | 2017-11-22 | 2020-07-07 | 中国地质大学（武汉） | 一种相干光ofdm通信***中定时同步方法、设备及存储设备 |
| CN109561042B (zh) * | 2018-12-17 | 2021-07-02 | 电子科技大学 | 一种ofdm***接收机的定时频率同步方法 |
| CN110881012B (zh) * | 2019-11-18 | 2022-05-31 | 展讯通信（上海）有限公司 | 一种ofdm调制信号的时间估计方法及装置 |
| CN111539323B (zh) * | 2020-04-22 | 2024-03-15 | 中国人民解放军战略支援部队信息工程大学 | 一种循环前缀线性调频信号的频率估计方法与装置 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102932307A (zh) * | 2012-10-26 | 2013-02-13 | 中国科学院上海微***与信息技术研究所 | 一种利用cazac序列的ofdm***时域同步方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080310539A1 (en) * | 2007-06-15 | 2008-12-18 | Kowalski John M | Systems and methods for generating an orthogonal signal from sequences that are not multiples of 2n |
-
2014
- 2014-08-12 CN CN201410394369.5A patent/CN104125188B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102932307A (zh) * | 2012-10-26 | 2013-02-13 | 中国科学院上海微***与信息技术研究所 | 一种利用cazac序列的ofdm***时域同步方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104125188A (zh) | 2014-10-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104125188B (zh) | 一种基于Zadoff‑Chu序列的OFDM频率同步方法 | |
| CN101588338B (zh) | 一种适用于分组传输的ofdm载波频偏估计方法 | |
| CN102694763B (zh) | 一种用于td-lte***整数倍频偏估计的方法 | |
| CN101083645B (zh) | 一种低复杂度ofdm快速同步的方法 | |
| CN105357160B (zh) | 发送参考信号的方法及装置、接收参考信号的方法及装置 | |
| CN101437005A (zh) | 通信同步过程中存在定时误差的整数倍频偏估计方法 | |
| CN102882670A (zh) | 一种基于cmmb信号的同步处理方法 | |
| CN102664859A (zh) | 多频带ofdm超宽带接收机同步与信道估计方案 | |
| CN102413091B (zh) | 一种基于对称cazac序列的ofdm***同步方法 | |
| CN103441778B (zh) | 一种基于cp平顶法的td-lte小区搜索同步方法 | |
| CN103379082A (zh) | Lte通信***中的时频同步方法和装置 | |
| CN104320367B (zh) | 一种对接收信号进行定时同步频偏估计和信道估计的方法 | |
| CN105516045A (zh) | 一种ofdm训练序列构造及同步方法 | |
| CN107086974A (zh) | 一种高动态环境下的ofdm同步方法 | |
| Zhou et al. | OFDMA initial ranging for IEEE 802.16 e based on time-domain and frequency-domain approaches | |
| CN101001235A (zh) | 一种弱能量并行pn序列时间同步与频率同步方法 | |
| CN103095613A (zh) | Sc-fde***中具有超低复杂度的整数倍频偏与信道参数联合估计算法 | |
| CN104270333A (zh) | 产生ofdm同步训练序列的方法及ofdm同步方法 | |
| Shao et al. | Robust timing and frequency synchronization based on constant amplitude zero autocorrelation sequence for OFDM systems | |
| CN101883067A (zh) | 一种适用于ofdm时间同步的训练序列构造及同步算法 | |
| CN104717168A (zh) | 正交频分复用超宽带***抗多径定时同步方案 | |
| CN100355255C (zh) | 基于统计平均的大搜索范围ofdm***载波同步方法 | |
| CN112910815B (zh) | 基于正负极性前导码的广义频分复用***时频同步方法 | |
| CN103139132A (zh) | 一种基于循环前缀和导频序列的联合同步方法 | |
| CN103095627A (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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Wu Yucheng Inventor after: Wang Yan Inventor after: Xia Yang Inventor after: Hu Zhen Inventor after: Li Weiqi Inventor after: Cao Meng Inventor after: He Xiang Inventor before: Wu Yucheng Inventor before: Hu Zhen Inventor before: Li Weiqi Inventor before: Cao Meng Inventor before: He Xiang |