CN102571679B - 一种用于突发ofdm系统中的符号粗同步的方法 - Google Patents
一种用于突发ofdm系统中的符号粗同步的方法 Download PDFInfo
- Publication number
- CN102571679B CN102571679B CN201210054559.3A CN201210054559A CN102571679B CN 102571679 B CN102571679 B CN 102571679B CN 201210054559 A CN201210054559 A CN 201210054559A CN 102571679 B CN102571679 B CN 102571679B
- Authority
- CN
- China
- Prior art keywords
- sequence
- delay
- calculation
- mean
- snr
- 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
Landscapes
- Synchronisation In Digital Transmission Systems (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明具体涉及突发OFDM系统的帧同步方法,本发明提出的方法,通过判断若干段接收到的数据序列的信噪比的平均值是否超过预设门限的办法,进而确定帧的起始位置。本发明提出的方法,很好的解决了SC算法的平台效应,并且不用计算较复杂的相关性,比起SC算法在计算量上也有优势。
Description
技术领域
本发明属于无线通信系统中的同步技术领域,具体涉及突发OFDM系统的符号粗同步方法。
背景技术
对于突发OFDM系统而言,接收端的首要工作是检测帧的起始位置,即:符号粗同步,又称帧同步。目前,研究帧同步方法的文献很多,大致可以分为两种,一是基于训练序列的帧同步,二是基于循环前缀的帧同步。由于基于循环前缀的帧同步在多径快变信道下的性能不理想,所以大部分的研究者都集中研究基于训练序列的帧同步。1997年12月,Timothy M.Schmidl和Donald C.Cox在IEEE Transactions on Communications第45卷第12期上发表的“RobustFrequency and Timing Synchronization for OFDM”中提出的基于训练序列的帧同步算法(以下简称SC算法)以及训练序列的时域结构被后来很多研究者借鉴,本发明也不例外。
SC算法中训练序列的大小是一个OFDM符号,其时域结构如附图2所示。其中,N是OFDM子载波个数,NCP是循环前缀的长度,训练序列的时域结构由两个完全相同的部分组成,图中两个B组成了训练序列,A表示循环前缀(cyclicprefix,简称CP),是由B的最后NCP个数据点组成。
SC算法利用训练序列前后两部分完全相同的特性来检测帧的起始位置。实际计算中,SC算法主要包括以下几个步骤:
步骤1:初始化接收序列的起始位置d为1,即d=1。
步骤2:从接收序列R中选取以d为起点,长度为N的一段序列S,即S=R(d:d+N-1)。
步骤3:计算序列S的前N/2点和后N/2点的相关值P(d)。P(d)的计算公式是
步骤4:计算序列S的后N/2点的能量和Q(d)。Q(d)的计算公式是
步骤5:计算时间判决变量MSC(d)。MSC(d)的计算公式是
步骤6:将MSC(d)和事先设定的门限值Th比较,若MSC(d)≥Th,则判断位置d帧的起始位置;若MSC(d)<Th,则更新d=d+1,返回步骤2。
经过分析发现,SC算法在实际应用中存在如下难点:
一是由于循环前缀CP的存在,时间判决变量MSC(d)会出现一个平台,且平台的长度与CP的长度成正比例关系,平台的存在会增加帧同步的判决困难。二是由于要计算序列S的前后两部分的相关值P(d)、及其与能量值Q(d)的比值,增加了实现过程中的复杂性。三是门限Th的设定通常随信噪比的变化而变化,因此,设定固定的门限值难以满足不同信道环境中不同信噪比条件下的性能要求。
发明内容
为了解决现有的方法存在的上述技术问题,本发明提出了一种用于突发的OFDM系统帧同步的方法,依次包括:
1)初始化:OFDM系统子载波数为N,循环前缀的长度为NCP,选取数据序列的起始位置d=1,判决门限Th,延迟处理点数Delay;
2)从接收到的数据序列R中选取以d为起点,长度为N的序列S,即S=R(d:d+N-1);
还依次包括如下步骤:
3)计算并保存序列S的信噪比
4)若d≥Delay,则进入5)步;否则更新d=d+1,返回2)步;
5)计算最近保存下来的Delay个snr(d)中最先保存的NCP个snr(d)的算术平均值Mean;
6)若Mean≥Th,则进入7)步;若Mean<Th,则更新d=d+1,返回2)步;
7)找出最近保存下来的Delay个snr(d)中的最大值所对应的序号dmax,判定dmax为符号粗同步位置。
为了解决上述技术问题,与上述第一方面相结合,本发明的第二方面在于提供一种用于突发的OFDM系统帧同步的方法,所述的判决门限Th的取值范围0.1≤Th≤33。
为了解决上述技术问题,与上述第一方面相结合,本发明的第三方面在于提供一种用于突发的OFDM系统帧同步的方法所述的延迟处理点数Delay≥2×NCP。
为了解决上述技术问题,与上述第一方面相结合,本发明的第四方面在于提供一种用于突发的OFDM系统帧同步的方法所述的计算序列S的信噪比snr(d)进一步包括如下步骤:
a)计算序列S的噪声功率
b)计算序列S的总功率
c)计算序列S中有用信号的功率s(d)=E(d)-n(d);
d)计算序列S的信噪比
附图3给出了现有的SC算法的时间判决变量MSC(d)与本发明技术方案的时间判决变量snr(d)的仿真结果的对比图,从图3(a)中可以看出,现有的SC算法的时间判决变量MSC(d)会产生平台效应,而从图3(b)可以看出,本发明技术方案中的时间判决变量snr(d)是序列S的有用信号的功率s(d)和噪声功率n(d)的比值,只有一个峰值,没有平台效应。所以本发明技术方案很好的解决了SC算法的平台效应。SC算法的门限值随信噪比的变化而变化,而本发明技术方案通过选择合理的恒定门限值Th,可以在不同的信噪比条件下工作。此外,本发明技术方案不用计算较复杂的相关性,比起SC算法在计算量上有优势。
附图说明
附图1为本发明技术方案的总体流程图。
附图2为现有的SC算法使用的训练序列的时域结构图;
附图3为本发明技术方案相对于现有技术的有益技术效果的仿真结果对比图。
具体实施方式
下面具体说明本发明技术方案的实施方式。
本实施例中,OFDM子载波个数N为128,有效子载波个数为116,循环前缀的长度NCP为N/4,符号间隔0.25μs。多普勒频率为900Hz,载波频率偏移位1KHz,OFDM符号采用BPSK调制,信道采用城市环境(TU)信道模型。训练序列后面跟着25个OFDM符号。信噪比设置为4-16dB。门限值设置为Th=0.9。
本实施例中用到的训练序列的时域结构与SC算法中提出的相同,如附图1所示,训练序列的具体内容按照如下步骤产生:
1:选取训练序列的频域值为:
F=[0 0-1 0-1 0-1 0-1 0-1 0-1 0-1 0 1 0 1 0-1 0 1 0 1 0 1 0 1 0 1 0-1 0-1 0-1 0-1 0 1 0 1 0 1 0 1 0-1 0-1 0 1 0 1 0 1 0-1 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0-1 0-1 0 1 0-10-1 0-1 0 1 0-1 0 1 0 1 0-1 0 1 0-1 0-1 0 1 0 1 0 1 0-1 0 1 0 1 0 1 0-1 0 1 0-1 0 10-1 0 1 0]。对F归一化处理,使得有效子载波上的平均功率为1,即将F乘以归一化因子
2:由训练序列的频域值得到时域值:对进行IFFT运算得到训练序列的时域序列,为了保证时域与频域能量一致,将得到的时域序列乘以归一化因子最终得到的序列记作T,即为附图1中两个B级联的部分。
T=[0.2500 -0.3490+0.0437i -0.6970-0.9694i -1.1761-0.4145i-0.3808-0.8921i 0.4204-1.0557i 0.9670+0.8320i 0.4718-0.0256i -0.7803-0.8536i 1.1242-0.3020i 0.3511+1.2303i-0.4240-0.9399i -0.1371-0.5475i -0.6635+0.0272i 0.0713-0.5307i 0.7629+1.0435i 1.0000-1.0000i 0.3955+0.3119i0.0371-0.3431i -1.0605-0.2520i -0.2164+0.6596i 0.4464-1.2384i 0.1938-0.4298i 0.5510+0.7979i 0.2803+0.1464i0.7323-0.8496i -0.8048-0.1210i -0.1678+0.0285i 0.7344-0.6850i -0.2035-1.3111i -0.1185-0.0746i -0.8603-0.6118i-1.2500 -0.8603+0.6118i -0.1185+0.0747i -0.2035+1.3111i0.7344+0.6850i -0.1678-0.0285i -0.8048+0.1210i 0.7324+0.8496i 0.2803-0.1464i 0.5510-0.7979i 0.1938+0.4298i0.4464+1.2384i -0.2164-0.6596i -1.0605+0.2520i 0.0371+0.3431i 0.3955-0.3119i 1.0000+1.0000i 0.7629-1.0435i0.0713+0.5307i -0.6635-0.0272i -0.1371+0.5475i -0.4240+0.9399i 0.3511-1.2303i 1.1242+0.3020i -0.7803+0.8536i0.4718+0.0256i 0.9670-0.8320i 0.4204+1.0557i -0.3808+0.8921i -1.1761+0.4145i -0.6970+0.9694i -0.3490-0.0437i0.2500 -0.3490+0.0437i -0.6970-0.9694i -1.1761-0.4145i-0.3808-0.8921i 0.4204-1.0557i 0.9670+0.8320i 0.4718-0.0256i -0.7803-0.8536i 1.1242-0.3020i 0.3511+1.2303i-0.4240-0.9399i -0.1371-0.5475i -0.6635+0.0272i 0.0713-0.5307i 0.7629+1.0435i 1.0000-1.0000i 0.3955+0.3119i0.0371-0.3431i -1.0605-0.2520i -0.2164+0.6596i 0.4464-1.2384i 0.1938-0.4298i 0.5510+0.7979i 0.2803+0.1464i0.7324-0.8496i -0.8048-0.1210i -0.1678+0.0285i 0.7344-0.6850i -0.2035-1.3111i -0.1185-0.0747i -0.8603-0.6118i-1.2500 -0.8603+0.6118i -0.1185+0.0747i -0.2035+1.3111i0.7344+0.6850i -0.1678-0.0285i -0.8048+0.1210i 0.7324+0.8496i 0.2803-0.1464i 0.5510-0.7979i 0.1938+0.4298i0.4464+1.2384i -0.2164-0.6596i -1.0605+0.2520i 0.0371+0.3431i 0.3955-0.3119i 1.0000+1.0000i 0.7629-1.0435i0.0713+0.5307i -0.6635-0.0272i -0.1371+0.5475i -0.4240+0.9399i 0.3511-1.2303i 1.1242+0.3020i -0.7803+0.8536i0.4718+0.0256i 0.9670-0.8320i 0.4204+1.0557i -0.3808+0.8921i -1.1761+0.4145i -0.6970+0.9694i -0.3490-0.0437i]。
下面开始针对上述训练序列进行符号粗同步运算:
步骤1:初始化d为1,即d=1,初始化延迟处理的点数Delay。
步骤2:从接收到的数据序列R中选取以d为起点,长度为N的序列S,即S=R(d:d+N-1)。
步骤3:计算序列S的噪声功率n(d),n(d)的计算公式为
步骤4:计算序列S的功率E(d),E(d)的计算公式为
步骤5:计算序列S中有用信号的功率s(d),s(d)的计算公式为s(d)=E(d)-n(d)。
步骤6:计算序列S中时间判决变量snr(d),snr(d)的计算公式为若d≥Delay,则进入步骤7;若d<Delay,则更新d=d+1,返回步骤2。
步骤7:计算最近保存下来的Delay个snr(d)中最先保存的NCP个snr(d)的平均值Mean,Mean的计算公式为若Mean≥Th,则进入步骤8;若Mean<Th,则更新d=d+1,返回步骤2。
步骤8:找出最近保存下来的Delay个snr(d)中的最大值所对应的序号dmax,判定dmax为符号粗同步位置。
实现过程中,需要存储的snr(d)的个数是Delay,即延迟处理点数。Delay的值越大,需要的存储空间越大;Delay的值越小,需要的存储空间越小,但是Delay必须满足Delay≥2×NCP。门限值Th的范围是0.1≤Th≤33。当Delay较小时,Th应该设置为较大的值,但都不应超过Th和Delay的范围。
Claims (1)
1.一种用于突发OFDM系统中的符号粗同步的方法,依次包括:
1)初始化:OFDM系统子载波数为N,循环前缀的长度为NCP,选取数据序列的起始位置d=1,数据序列的时域结构由两个完全相同的部分组成,判决门限Th,延迟处理点数Delay,所述延迟处理点数Delay满足Delay≥2×NCP,所述门限值Th的范围是0.1≤Th≤33;
2)从接收到的数据序列R中选取以d为起点,长度为N的序列S,即S=R(d:d+N-1);
其特征是还依次包括如下步骤:
3)计算并保存序列S的信噪比
4)若d≥Delay,则进入5)步;否则更新d=d+1,返回2)步;
5)计算最近保存下来的Delay个snr(d)中最先保存的NCP个snr(d)的算术平均值Mean;
6)若Mean≥Th,则进入7)步;若Mean<Th,则更新d=d+1,返回2)步;
7)找出最近保存下来的Delay个snr(d)中的最大值所对应的序号dmax,判定dmax为符号粗同步位置;
其中,所述3)步进一步包括:
a)计算序列S的噪声功率
b)计算序列S的总功率
c)计算序列S中有用信号的功率s(d)=E(d)-n(d);
d)计算序列S的信噪比
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210054559.3A CN102571679B (zh) | 2012-02-23 | 2012-02-23 | 一种用于突发ofdm系统中的符号粗同步的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210054559.3A CN102571679B (zh) | 2012-02-23 | 2012-02-23 | 一种用于突发ofdm系统中的符号粗同步的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102571679A CN102571679A (zh) | 2012-07-11 |
CN102571679B true CN102571679B (zh) | 2015-05-06 |
Family
ID=46416168
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210054559.3A Expired - Fee Related CN102571679B (zh) | 2012-02-23 | 2012-02-23 | 一种用于突发ofdm系统中的符号粗同步的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102571679B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105024966B (zh) * | 2015-08-12 | 2018-05-25 | 西安电子科技大学 | Ofdm系统的符号定时同步方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5732113A (en) * | 1996-06-20 | 1998-03-24 | Stanford University | Timing and frequency synchronization of OFDM signals |
CN101883069A (zh) * | 2009-05-09 | 2010-11-10 | 电子科技大学中山学院 | 一种训练序列构造及其ofdm时间同步方法 |
CN102185819A (zh) * | 2011-05-09 | 2011-09-14 | 中国科学院上海微系统与信息技术研究所 | 基于共轭对称序列的ofdm时间同步算法 |
-
2012
- 2012-02-23 CN CN201210054559.3A patent/CN102571679B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5732113A (en) * | 1996-06-20 | 1998-03-24 | Stanford University | Timing and frequency synchronization of OFDM signals |
CN101883069A (zh) * | 2009-05-09 | 2010-11-10 | 电子科技大学中山学院 | 一种训练序列构造及其ofdm时间同步方法 |
CN102185819A (zh) * | 2011-05-09 | 2011-09-14 | 中国科学院上海微系统与信息技术研究所 | 基于共轭对称序列的ofdm时间同步算法 |
Non-Patent Citations (5)
Title |
---|
An improved preamble-based SNR estimation algorithm for OFDM systems;Zivkovic, M等;《Personal Indoor and Mobile Radio Communications (PIMRC), 2010 IEEE 21st International Symposium on》;20100930;172 - 176 * |
Low-complexity time-domain SNR estimation for OFDM systems;Ijaz, A等;《Electronics Letters 》;20110929;第47卷(第20期);1154 - 1156 * |
OFDM基带系统关键技术及其FPGA设计与实现;李斌;《湖南大学硕士学位论文》;20110509;正文第4.3.1节-4.4.7节 * |
Robust Frequency and Timing Synchronization for OFDM;Timothy M;《IEEE TRANSACTIONS ON COMMUNICATIONS》;19971212;第57卷(第12期);1613-1621 * |
SNR estimation algorithm based on the preamble for OFDM systems in frequency selective channels;Guangliang等;《Communications, IEEE Transactions on 》;20090818;第57卷(第8期);2230-2234 * |
Also Published As
Publication number | Publication date |
---|---|
CN102571679A (zh) | 2012-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101079688B (zh) | 一种正交频分复用系统中实现同步的方法 | |
CN101018087B (zh) | 移动通信系统的下行同步方法及装置 | |
CN103023853B (zh) | 联合训练序列互相关信息的正交频分复用接收帧同步方法 | |
CN102882670A (zh) | 一种基于cmmb信号的同步处理方法 | |
CN104767706A (zh) | 一种mimo ofdm定时同步装置 | |
CN103259756A (zh) | 一种应用于ofdm系统的符号定时同步和载波同步方法 | |
CN102377712A (zh) | 辅同步信号检测方法和装置 | |
CN106685880A (zh) | 物理帧中前导符号的接收处理方法 | |
CN105991501B (zh) | 前导符号的接收方法及装置 | |
CN101651650A (zh) | 一种同步与频偏的联合估计方法及装置 | |
CN103581072B (zh) | 多载波系统的频偏估计方法及装置 | |
CN102026231B (zh) | 一种用于无线通信系统的随机接入检测方法 | |
CN102461102A (zh) | Lte通信网络中根据pss信号的信道估计过程及其接收器 | |
CN101325575B (zh) | 一种正交频分复用系统中帧同步的实现方法 | |
CN105516045A (zh) | 一种ofdm训练序列构造及同步方法 | |
CN100477655C (zh) | 一种估计ofdm整数倍频偏的方法 | |
CN104902560A (zh) | 一种下行定时同步方法及装置 | |
CN102143574B (zh) | 一种适合于IEEE802.16m前导字结构的定时同步方法 | |
CN102377726B (zh) | Ofdm系统的定时同步方法 | |
CN102571679B (zh) | 一种用于突发ofdm系统中的符号粗同步的方法 | |
CN101047686A (zh) | 一种用于WiMAX系统的接收符号同步方法 | |
CN106230758A (zh) | 一种lte‑a系统整数倍频偏估计方法 | |
CN104717168B (zh) | 正交频分复用超宽带系统抗多径定时同步方案 | |
CN102594762B (zh) | 一种用于突发的ofdm系统帧同步的方法 | |
CN101599933A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150506 Termination date: 20190223 |