基于ZadoffChu序列的OFDM精确定时同步算法

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

㊀㊀文章编号:1009-2552(2019)01-0001-05㊀㊀DOI:10 13274/j cnki hdzj 2019 01 001
基于Zadoff ̄Chu序列的OFDM精确定时同步算法
周㊀平1ꎬ周思远1ꎬ刘思含2ꎬ吴玉成2
(1 扬州万方电子技术有限公司ꎬ江苏扬州225006ꎻ2 重庆大学微电子与通信工程学院ꎬ重庆400044)
摘㊀要:针对OFDM系统在多径衰落环境下通过能量大小定时而产生定时误差的问题ꎬ文中提出了一种基于Zadoff ̄Chu序列的第一径检测算法ꎬ利用序列相关性得到信道冲激响应估计值ꎬ此外分析了残留频偏对定时同步的影响ꎬ并通过改进自适应门限抑制残留频偏来判定准确定时位置ꎮ算法在时域完成ꎬ实现简单且运算量小ꎮ通过仿真对比其他算法可以看出ꎬ文中所提算法在高速移动多径环境中第一径检测概率高㊁均方误差小ꎬ特别是低信噪比下仍能准确同步ꎮ关键词:OFDMꎻ定时同步ꎻZadoff ̄Chu序列ꎻ多径衰落信道中图分类号:TN929 5㊀㊀文献标识码:A
AccuratetimingsynchronizationalgorithmbasedonZadoff ̄Chu
sequenceforOFDMsystem
ZHOUPing1ꎬZHOUSi ̄yuan1ꎬLIUSi ̄han2ꎬWUYu ̄cheng2
(1.YangzhouWanfangElectronTechnolCo.ꎬLtd.ꎬYangzhou225006ꎬJiangsuProvinceꎬChinaꎻ
2.SchoolofMicroelectronicsandCommunicationEngineeringꎬChongqingUniversityꎬChongqing400044ꎬChina)
Abstract:ThispaperproposedafirstpathdetectionalgorithmbasedonZadoff ̄Chusequenceinorderto
solvetimingerrorcausedbyOFDMsysteminmultipathchannel.ThecorrelationpropertyofZadoff ̄Chusequenceisutilizedtoestimatethechannelimpulseresponse.Inadditionꎬtheinfluenceofresidual
frequencyoffsetontimingsynchronizationisanalyzedꎬandtheaccuratetimingpositionisdeterminedbyimprovingtheadaptivethresholdandsuppressingtheresidualfrequencyoffset.Thealgorithmiscompletedinthetimedomainꎬwhichissimpletoimplementandsmallincomputation.Analysisandsimulationresultsshowthattheproposedalgorithmhashighprobabilityofdetectionandlowmeansquarederrorinhigh ̄speedmovingmultipathenvironmentꎬandatlowsignal ̄to ̄noiseratio.
Keywords:OFDMꎻtimingsynchronizationꎻZadoff ̄Chusequenceꎻmultipathfadingchannel
收稿日期:2018-10-09
基金项目:江苏省重大产业创新专项(BA2018114)
作者简介:周平(1951-)ꎬ男ꎬ高级工程师ꎬ研究方向为智能通信
技术ꎮ
0㊀引言
在新一代宽带移动通信中ꎬ正交频分复用(Or ̄
thogonalFrequencyDivisionMultiplexingꎬOFDM)系统以更高的数据传输速率和频谱利用率ꎬ较强的抗多径能力成为其核心技术[1]ꎮOFDM系统对精确定时同步要求严格ꎬ定时偏差会引入ISIꎬ因此同步技术的准确程度直接影响系统的性能[2-3]ꎮ
在多径衰落信道中ꎬ每条路径在不同时刻抵达
接收端ꎬ且能量最强径未必就是第一径ꎬ因此基于能量检测的同步算法[4-6]存在定时偏差ꎮ目前ꎬ很多学者对OFDM系统同步问题做了研究ꎬ其中最为经
典的是文献[4]所提方法ꎬ设计了具有两个相同序列的训练符号用于同步ꎬ但是该算法的定时度量函数在准确同步位置处没有明显峰值ꎬ其左右样点均有较高峰值导致定时精度不高ꎮ文献[5]在文献[4]算法基础上修改训练符号结构ꎬ将其变为相同的4段序列且前两段与后两段正负号相反ꎬ该方法有较大副峰ꎮ文献[6]则把结构设计为前后序列共轭对称ꎬ准确定时位置峰值更为明显但其左右仍

有副峰ꎮ文献[7]利用具有良好相关性的CAZAC序列进行符号定时和频偏估计ꎬ但其定时度量函数仍有明显旁瓣ꎬ致使低信噪比下检测概率不高ꎮ文献[8]中算法利用信道冲激响应估计定时位置ꎬ然而其定时估计存在峰值平台ꎮ文献[9]考虑了第一径不是能量最强径的情况ꎬ利用滤波器平滑定时函数后采用门限判决确定第一径ꎬ然而滤波器参数和门限值均需多次验证ꎬ适应性不强ꎮ文献[10]所提算法定时度量函数自相关峰值尖锐且旁瓣干扰小ꎬ但未考虑到达第一径问题ꎬ系统在多径衰落信道中定时精度不高ꎮ文献[11]未考虑频偏对定时的影响且门限不能自适应ꎬ因此在信噪比较低时检测性能不高ꎮ文献[12]所提算法采用了自适应门限检测第一径ꎬ但忽略了频偏对序列相关性的影响ꎮ本文主要研究基于Zadoff ̄Chu序列自相关性的定时同步算法ꎬ提出了一种精确定位多径衰落信道中第一径的同步算法ꎬ利用Zadoff ̄Chu序列自相关峰值尖锐的特性进行信道估计ꎬ设计一种可解决残留频偏带来虚警概率增加问题的自适应门限ꎬ提高了判决准确性ꎮ该方法解决了高速移动宽带通信系统定时不准确问题ꎬ且无需进行时频变换ꎬ方便实现ꎮ
1㊀OFDM系统模型
OFDM系统发送端对信息比特进行信道编码和交织ꎬ经串并转换后调制到各个子载波上ꎬ为减少符号间干扰(Inter ̄SymbolInterferenceꎬISI)ꎬ应添加循环前缀(CyclicPrefixꎬCP)ꎮ若考虑循环前缀ꎬ频域信号经过快速逆傅里叶变换(InverseFastFourierTransformationꎬIFFT)得到时域表达式为:
x(n)=1NðN-1k=0X(k)ej2πk(n-Ncp)/Nꎬ
0ɤnɤN+Ncp-1(1)式中ꎬX(k)为调制在第k个子载波上的数据ꎬN是系统子载波总数ꎬNcp是CP长度ꎮ
假设系统无载波频偏且收发两端采样时钟均已同步ꎬ经过多径信道并存在噪声影响ꎬ同时考虑到系统存在定时偏差dꎬ采样后的接收信号可表示为:r(n)=ðL-1l=0h(l)x(n-d-l)+w(n)(2)其中ꎬL为多径分量个数ꎻh(l)为信道冲激响应ꎻw(n)为加性高斯白噪声ꎮ
2㊀OFDM系统定时同步算法
2.1㊀训练符号结构
本文采用Zadoff ̄Chu序列构建训练序列ꎬ式(3)给出其序列值Zn为:
Zn=expjπrn2
Np
()(3)
式中ꎬNp(Np为偶数)为序列周期ꎬr为正整数ꎬr与Np互质ꎮ若Np=N/2ꎬr=1ꎬ其自相关函数为:Rτ=ðNp-1n=0ZnZ∗(n+τ)modNp=ðNp-τ-1n=0ZnZ∗(n+τ)+ðNp-1n=Np-τZnZ∗(n+τ)=Np
ꎬτ=0
0ꎬτʂ0
{(4)
上式中mod表示取模ꎮ由式(4)可以看出ꎬZadoff ̄Chu序列具有循环自相关[13]峰值尖锐㊁旁瓣为0的特点ꎬ在多径衰落信道下ꎬ利用该序列的优良特性可以改善同步效果ꎮ同时ꎬ峰均比低㊁模值恒定等特点都使其在同步算法中得到广泛应用ꎬ极大地增强系统定时准确性ꎮ
本文所采用的训练符号结构如图1所示ꎮ训练符号中将Zn重复两段ꎬ并截取Zn尾部序列作为CPꎬ
在定时同步过程中可利用前后序列相同的特性实现符号定时ꎮ
图1㊀训练符号结构
2.2㊀定时同步算法
定时度量函数的选取会影响系统定时准确程度ꎬ在准确定时点附近ꎬ定时度量函数应尽可能尖锐ꎮ根据图1训练符号的特点ꎬ利用训练序列前后重复的结构特点及其同本地已知序列的相关性ꎬ设计定时度量函数为:
M(d)=1NpðNp-1n=0r(d+n)Z∗(n)(5)式中ꎬr(n)是经过采样的接收信号ꎬZ∗(n)是本地训练序列的复共轭ꎬ由相关性可知在|M(d)|峰值处所对应的d就是定时位置ꎬ即:
d^=argmaxd{|M(d)|}(6)在实际的多径通信环境下ꎬ接收信号的能量最强径不一定是到达第一径ꎬ然而d
^求得的是能量最强径对应的定时位置ꎮ由图2可知ꎬ若能量最强径是到达第二径ꎬ定时位置出现在到达第一径的数据部分ꎬ在进行FFT变换时则会产生符号间干扰[14]ꎮ图2给出定时位置在第二径时引起ISI干扰的示意图ꎮ分析可知ꎬ为保证系统性能ꎬ定时位置必须在第一径上ꎮ
为确保定时位置在到达第一径上ꎬ下面给出第

图2㊀定时误差引起ISI的示意图
一径检测方法:由能量大小确定的定时位置d^
得到搜索区间[d^-Lꎬd^]ꎬ如果当前位置d^
不是第一径ꎬ那么在其前面L点内一定可以搜索到第一径ꎬ将每个样点处的定时度量函数M(d)都同所设定的门限相比ꎬ如果M(d)高于门限ꎬ则认为检测到了第一径ꎮ门限值的选取对定时位置判定至关重要ꎬ这就要求定时度量函数M(d)能正确反映信道冲激响应特征ꎮ
忽略发射机与接收机间的传输时延ꎬ将式(5)中定义的定时度量函数展开可得:
M(d)
=1
NpðNp-1
n=0
r(d+n)Z∗(n)=1NpðNp-1
n=0
ðL-1
l=0
h(l)Z(d+n-l)Z∗(n)+

NpðNp-1
n=0
w(d+
n)Z∗
(n)=
ðL-1
l=0
h(l)δ(d-l)+wᶄ(d)=h(d)+
wᶄ(d)
(7)
其中ꎬδ( )为单位冲激函数ꎮ由中心极限定理可知ꎬ当Np充分大时ꎬwᶄ(d)~CN(0ꎬσ2/Np)ꎮ根据上述M(d)和h(d)的关系得到定时度量
函数M(d)的条件概率密度函数为:p(M(d)|h(d))=

2π/Npσ
exp-Np2σ2[M(d)-h(d)]2
{}
(8)继而计算其最大似然估计方程为:
∂lnp(M|h)
∂h=Np
[M(d)-h(d)]σ2
h(d)=h^
ML(d)
=0(9)
因此ꎬM(d)是信道冲激响应h(d)的最大似然估计ꎬ即:
h^
ML(d)=M(d)
(10)
SUI ̄3信道下信道冲激响应估计h^
ML(d)的均
方误差(MeanSquaredErrorꎬMSE)曲线如图3所示ꎮ随着信噪比(SignaltoNoiseRatioꎬSNR)的提升ꎬMSE值大幅度减小ꎬ当SNR高于3dB时ꎬh^
ML
(d)已经趋近于0ꎮ
图3㊀信道冲激响应估计MSE
2.3㊀自适应门限设置
计算出信道冲激响应的估计值h^
ML后ꎬ根据设定的门限值判断是否为到达第一径ꎮ传统做法中使用的是固定门限ꎬ然而在实际通信环境中ꎬ传输信号会受到衰落及噪声的影响而产生动态变化ꎬ造成系统定时过程中虚警概率高或漏检概率高ꎮ该检测到达第一径的方法在本质上等效于峰值检测ꎬ利用匹配滤波器得到相关峰值ꎮ参考文献[15]中的自适应门限检测方法ꎬ本文将自适应门限设置为:
VT=
(P
-1/Jfa
-1)ðJ
d=1
ud(t)
(11)
其中ꎬPfa为系统虚警概率ꎬJ表示噪声功率估计中相互独立的统计判决量ud(t)的个数ꎬud(t)=r2dI(t)+r2dQ(t)=|M(d)|2
ꎬrdI(t)和rdQ(t)分别是
I㊁Q两路的接收信号ꎮ
在高速移动通信环境中系统存在多普勒频移ꎬ
定时位置d^
受频偏影响发生偏移而产生定时误差ꎬ因此在计算定时度量函数M(d)前必须对接收信号进行频偏补偿ꎮ在SNR较低情况下ꎬ噪声对定时同步影响很大ꎬ因此采用基于噪声功率的自适应门限ꎬ在虚警概率不变时ꎬ可以使检测概率提升至最高ꎮ可是在高SNR情况下ꎬ噪声对定时的影响可以忽略不计ꎬ此时影响定时精度的主要因素是系统中的残留频偏[16]ꎮ
此处不考虑信道衰落和噪声的影响ꎬ假设系统
存在残留频偏Δεꎬ那么在接收端通过匹配滤波器得到的相关函数为:
RΔε(τ)=ðNp-1
n=0
ejπn2
/Np
ej2πΔεmaxn/N
e-jπ(n+τ)
2/N

=e
-jπτ2/Np
1-ejπΔεmax
1-ejπ(Δεmax-2τ)/Np
=Np1-ejπΔεmax
1-ejπΔεmax/Npꎬ㊀㊀㊀τ=0e
-jπτ2/Np1-ejπΔεmax
1-ejπ(Δεmax-2τ)/Np
ꎬτʂ0ìîí
ïïïï(12)

式中ꎬτ是采样位置相对于定时同步点偏移的样点数ꎬΔεmax是将子载波归一化后的系统最大残留频偏ꎮ
对比式(4)与式(12)的相关函数ꎬ可以看出残
留频偏的存在导致自相关峰值降低且互相关性遭到破坏ꎬ严重影响了Zadoff ̄Chu序列的循环相关特性ꎮ当τʂ0时ꎬ若系统存在残留频偏ꎬ序列相关性遭到破坏ꎬ非准确同步位置附近会产生多个副峰ꎬ若副峰值大小超过门限值ꎬ则会引起系统虚警概率增加ꎮ为改善上述问题ꎬ本文通过在所求自适应门限VT上增加一固定偏移量ΔV得到新门限VᶄT
VᶄT=VT+ΔV
(13)
ΔV的取值与相关峰值有关ꎬ式(12)中ꎬ越靠近
主峰的峰值幅度越大ꎮ在τʂ0情况下ꎬτ=1时对应的RΔε(τ)取最大值ꎬ故设置偏移量ΔV=1
NpRΔε
(1)ꎬ对其能量最大径幅度进行归一化为:ΔV=1Npe
-jπ/Np1-ejπΔεmax
1-ejπ(Δεmax-2)/Np
M(d^)(14)
结合2.2节可知ꎬd^
是能量最强径对应的定时位置ꎬ当能量最强径不是到达第一径时ꎬ在搜索区间[d^-Lꎬd^
-1]内至少存在一点ꎬ其定时度量值高于门限值ꎬ即M(d)>VᶄTꎮ如果最大多径时延扩展L未知ꎬ就设定L=Ncpꎬ原因是一般情况下CP长度不小于最大时延扩展ꎬ在两者长度相等时系统性能达到最优ꎮ此时ꎬ到达第一径定时位置估计为:
d^ᶄ=argmin

{
M(d)
>Vᶄ

d^
-Ncp+1ɤdɤd

}
(15)
3㊀算法性能仿真
为验证本文算法有效性ꎬ利用MATLAB搭建
OFDM系统ꎮ参数设置如下:系统FFT点数为
1024ꎬ循环前缀长度Ncp=128ꎬ虚警概率Pfa=0 0001ꎬ每个SNR仿真10000次ꎬ采样时钟fs=
80MHzꎮ多径衰落信道模型采用SUI信道模型ꎬ多普勒频移设置较大ꎬ到达路径数为3ꎬ满足高速移动环境要求ꎮ为验证本文算法检测第一径的性能ꎬ设定第二径为能量最强径ꎬ具体参数如表1所示ꎮ
表1㊀多径信道参数
SUI ̄3信道
SUI ̄5信道
多径信道增益(dB)[-5㊀0㊀-10][-5㊀0㊀-10]
多径时延(μs)[0㊀0.5㊀1]
[0㊀5㊀10]
最大多普勒(Hz)
5001000K因子
[1㊀0㊀0]
[0㊀0㊀
0]
图4㊀SUI ̄3
信道下检测概率
图5㊀SUI ̄5信道下检测概率
㊀㊀图4和图5对比了本文算法与Kishore&Reddy
算法[8]和Yang算法[12]的第一径检测概率曲线ꎮ在SUI ̄3和SUI ̄5信道中ꎬ本文所提算法的检测概率都明显高于其他两种算法ꎬ不论高低信噪比情况下均有极大优势ꎮ在低信噪比情况下也能保持
90%以上的检测概率ꎬ表明基于噪声功率的自适应门限适用于当前条件ꎮ高信噪比情况下ꎬ本文算法在SUI ̄3信道下的检测概率更是接近100%ꎬ此时影响检测性能的主要因素就是残留频偏ꎬ而本文在门限设置过程中考虑了残留频偏问题ꎬ降低了其对系统性能的影响ꎮYang算法之所以在衰落信道中检测概率低ꎬ是因为其没有考虑频偏会恶化训练符号相关性ꎮ
对比图4和图5检测性能可知ꎬ本文算法在SUI ̄5信道下的检测概率较SUI3信道有所降低ꎬ原因是SUI ̄5信道环境更为恶劣ꎬ多普勒频移更为显著ꎬ此时系统残留频偏也相应变大ꎬ检测性能受其影响而下降ꎮ
均方误差(MeanSquareErrorꎬMSE)是衡量系统差异的重要参数ꎬ图6和图7对比了两种SUI信道下本文算法和另外两种算法定时位置的MSE值ꎮ本文算法和Yang算法在低信噪比下仍有较好的MSE性能ꎬ然而Kishore&Reddy算法在求信道冲激响应的估计值时忽略了噪声功率的影响ꎬ导致该算

图6㊀SUI ̄3信道下
MSE
图7㊀SUI ̄5信道下MSE
法在低信噪比条件下MSE值非常大ꎮ由于定时度量函数的副峰值在准确定时位置附近ꎬ所以Yang算法的MSE值非常低ꎮ本文算法考虑了噪声功率和残留频偏的影响ꎬ改进了基于噪声功率的自适应门限ꎬ使其在保证了高检测概率的同时ꎬ还有极低的MSE值ꎬ在高速移动环境下可实现精确定时ꎮ
4㊀结束语
在OFDM系统中ꎬ符号定时偏移会导致系统性
能下降ꎮ现有定时算法多以能量判定定时位置ꎬ然而到达第一径不一定是能量最强径ꎬ为解决以上情况引起系统定时误差的问题ꎬ本文提出一种基于Zadoff ̄Chu序列的定时同步算法ꎮ本文算法利用Zadoff ̄Chu序列良好的循环自相关特性定义定时度量函数ꎬ通过理论推导得到信道冲激响应估计与定时度量函数的关系ꎬ采用改进了的自适应门限技术实现第一径检测ꎬ有效降低了残留频偏对系统检测性能的影响ꎮ仿真结果表明ꎬ本文所提算法容易实现ꎬ在多径衰落环境下可保证较高检测性能㊁均方误差小ꎬ满足高速移动通信要求ꎮ参考文献:
[1]CuiTꎬZhouQꎬLiJꎬetal.Theanalysisandresearchofadaptive
transmissionalgorithmsinwirelesscommunicationsystem[C].IEEEInternationalConferenceonCommunicationProblem ̄Solving.IEEEꎬ2016:615-618.
[2]WangWJꎬTianYꎬAhnCJ.Blindsymboltimingsynchronization
schemeforreal ̄timeindustrialwirelesscontrolnetworkusinghigh ̄speedpreamblelessOFDMsystemsovermultipathfadingchannels[J].DigitalSignalProcessingꎬ2017ꎬ62:30-37.
[3]PandaMꎬMishraS.Performanceanalysisofconcomitantsynchroni ̄
zationalgorithmsforOFDMandnon ̄contiguousOFDMsystems
[C].InternationalConferenceonManandMachineInterfacing.IEEEꎬ2018.
[4]SchmidlTMꎬCoxDC.Robustfrequencyand ̄timingsynchroniza ̄
tionforOFDM[J].IEEETransactionsonCommunicationsꎬ1997ꎬ
45(12):1613-1621.
[5]MinnHꎬZengMꎬBhargavaVK.Ontimingoffsetestimationfor
OFDMsystems[J].IEEECommunicationsLettersꎬ2000ꎬ4(7):
242-244.
[6]ParkBꎬCheonHꎬKangCꎬetal.Anoveltimingestimationmethod
forOFDMsystems[J].IEEECommunicationsLettersꎬ2003ꎬ7(5):
239-241.
[7]FangYꎬZhangZꎬLiuG.ANovelSynchronizationAlgorithmBased
onCAZACSequenceforOFDMSystems[C].InternationalConfer ̄enceonWirelessCommunicationsꎬNETWORKINGandMobileComputing.IEEEꎬ2012:1-4.
[8]KishoreCNꎬReddyV.Atechniquefordominantpathdelayestima ̄
tioninanOFDMsystemanditsapplicationtoframesynchronization
inOFDMmodeofWMAN[J].EURASIPJournalonWirelessCom ̄municationsandNetworkingꎬ2006(2):1-8.
[9]张洁ꎬ蔡鹏ꎬ张平.一种适用于多径衰落信道的OFDM定时同步
算法[J].北京邮电大学学报ꎬ2005ꎬ28(2):105-108.
[10]刘彬ꎬ罗志年ꎬ彭疆.改进CAZAC序列的OFDM时频同步算
法[J].计算机工程与应用ꎬ2017ꎬ53(14):76-79ꎬ98.[11]DasAꎬMohantyBꎬSahuB.ANovelCAZACSequenceBasedTim ̄
ingSynchronizationSchemeforOFDMSystem[J].International
JournalofAdvancesinAppliedSciencesꎬ2018ꎬ7(1):66-72.
[12]YangFꎬZhangX.RobustTime ̄DomainFineSymbolSynchroniza ̄
tionforOFDM ̄BasedPacketTransmissionUsingCAZACPreamble[C].2013IEEEMilitaryCommunicationsConference(MILCOM
2013).SanDiegoꎬ2013:436-440.
[13]BlumensteinJꎬBobulaM.CoarseTimeSynchronizationUtilizing
SymmetricPropertiesofZadoff ̄ChuSequences[J].IEEECommu ̄nicationsLettersꎬ2018(99):1.
[14]GulMMUꎬMaXꎬLeeS.TimingandFrequencySynchronization
forOFDMDownlinkTransmissionsUsingZadoff ̄ChuSequences
[J].IEEETransactionsonWirelessCommunicationsꎬ2015ꎬ14(3):1716-1729.
[15]吴玉成ꎬ陈婷婷.在滑动窗口中判决的自适应门限检测方
法[J].系统仿真学报ꎬ2008ꎬ20(10):2770-2773.
[16]HuangHꎬWangWGJꎬHeJ.Phasenoiseandfrequencyoffset
compensationinhighfrequencyMIMO ̄OFDMsystem[C].IEEEInternationalConferenceonCommunications.IEEEꎬ2015:1280-
1285.
责任编辑:杨静
5。

相关文档
最新文档