专利名称：Ofdm同步训练序列的生成方法和基于该训练序列的同步方法
技术领域：
本发明涉及通信技术领域，特别涉及一种OFDM(正交频分复用)同步训练序列的生成方法和基于该训练序列的同步方法。
背景技术：
目前，OFDM技术在越来越多的有线、无线通信领域得到应用，如DAB(数字音频广播)、DVB(数字视频广播)、IEEE802.11a、HIPERLAN/2、ADSL(非对称用户环路)等等。这主要由于OFDM技术具有许多优势有效对抗多径干扰和窄带干扰，频谱利用率高，数据传输速率高等。
OFDM系统对于同步偏差非常敏感。同步偏差主要分为符号定时偏差和频率偏差。符号定时偏差是指解调OFDM符号时，FFT窗口提前或是滞后。如果符号定时偏差过大，使定时的偏移量与最大时延扩展的长度之和小于循环前缀的长度，会产生符号间干扰，破坏OFDM符号的完整性，降低系统的性能。频率偏差fΔ是指接收端的解调载频与发送端的调制载频不一致，这主要由于晶振不稳定或是多普勒频移。频率偏差又分为子载波间隔的小数倍频率偏差fF和子载波间隔的整数倍频率偏差fI在下文中分别简称为小数频偏和整数频偏。其中，小数频偏会造成子载波间干扰(ICI)；整数频偏不会引起ICI，但会引起接收数据符号的循环移位，使得解调出来的信息符号的错误概率为50％。因此，有效地估计符号定时偏差和频率偏差，是OFDM系统能否正常工作的关键因素。
图1是一个通用OFDM基带系统的框图。图1中的发送端串并转换、IDFT(傅里叶反变换)、发送端并串转换、插入循环前缀，插入同步信息等模块，代表OFDM基带调制过程；图1中的分离同步信息，去除循环前缀，接收端串并转换，DFT(傅里叶变换)以及频域均衡，接收端并串转换等模块是与之对应的OFDM基带解调过程；图1中的数模转换，发送滤波处理，信道，接收滤波处理，模数转换等模块是模拟发射与信道环节；图1中的同步单元模块即是实现OFDM同步的部件。现结合图1对OFDM系统的调制和解调过程进行简要介绍
设基带采样时间间隔为Ts，OFDM系统的有效符号点数为N(一般情况下，N＝2β)，与之对应的有效符号周期为T＝NT；Nu是N个OFDM子载波中有效子载波的数目。在发送端，频域数据{ai，k|(k＝-N/2，...，0，1，...，N/2-1)}(ai，k是第i个符号，第k个子载波上加载的数据)被置入图1中的发送端串并转换模块；经过IDFT模块，在发送端并串转换模块的输出端得到时域数据{bi，l|(l＝0，1，2，...，N-1)}(bi，j是第i个符号，第l个采样点的数据)；为了对抗符号间干扰(ISI)，添加Ng点的循环前缀(图1中的插入循环前缀模块)，于是，每个OFDM符号含点数为Nsym(Nsym＝N+Ng)；再插入同步训练序列(图1中的插入同步信息模块)，然后整个信号通过信道到达接收端；在接收端，分离出同步训练序列(图1中的分离同步信息模块)送入图1中的同步单元模块；再删去循环前缀(图1中的去除循环前缀模块)，在图1中的接收端串并转换模块的输出端得到时域数据{ri，m|(m＝0，1，2，...，N-1)}(ri，m是第i个符号，第m个采样点的数据)；最后经过图1中的DFT和频域均衡模块，在图1中的接收端并串转换模块的输出端得到解调的频域数据{zi，n|(n＝-N/2，...，0，1，...，N/2-1)}(zi，n是第i个符号，第n个子载波上解调的数据)。
常见的OFDM同步方案有两种(1)基于特定频域序列和2等分结构的同步训练序列，该方法通过计算训练符号的半个符号延时的自相关，寻找相关峰值获得符号定时，求相关峰值处的相角估计小数频偏，再对同步训练符号进行小数频偏补偿后，作快速傅里叶变换(FFT)，然后与已知频域序列作循环移位相关，通过寻找相关峰来估计整数频偏。该方法所能估计的频偏范围，随着相关搜索范围的扩大而增大，但是，其计算复杂度也会不断上升。因此，由于计算复杂度的原因，该方法难以应用到实际系统中。参见文献Schmidl，T.M.等“Low-overhead，low-complexity [burst]synchronization for OFDM”，IEEE International Conference onCommunications，Volume 3，June 1996，Page(s)1301-1306 (“低数据开销、低复杂度的OFDM同步方法”IEEE国际通信技术会议)。
(2)基于L等分(L＝2v)的OFDM同步训练符号结构，通过计算训练符号的特定延时的自相关，寻找相关峰值获得符号定时，再求相关峰值处的相角来估计频偏。该方法虽然计算复杂度低，但其时间同步性能较差，且频偏估计范围较小，仅为[-L/2，L/2)。在实际中，L的取值一般不能太大。因为过大的L，会造成OFDM系统的峰均功率比过高，并产生较大的频偏估计误差，以及训练符号的频域数据稀疏性导致其抵抗频率选择性衰落的能力下降。参见文献Heiskala J，Terry JOFDM Wireless LANs-A Theoretical and Practical Guide.[M].Indianapolis USAPearson Education Inc，2002.70-73(《OFDM无线局域网——理论与实践的指导》)。
因此，方法(1)属于频域算法，需要对整个同步训练符号作FFT，其主要缺点是复杂度高；方法(2)属于时域算法，虽然不需要对整个同步训练符号作FFT，具有较低的复杂度，但其主要缺点是性能较差，这使上述两种方法在实际应用中受到限制。
发明内容
本发明的目的在于针对现有技术的不足，提供一种OFDM同步训练序列的生成方法和基于该训练序列的同步方法。本发明所述的OFDM同步训练序列的生成方法为由2个OFDM符号级联生成OFDM同步训练序列，所述2个OFDM符号分别是非正则OFDM符号与正则OFDM符号。本发明所述的同步方法基于根据所述OFDM同步训练序列的生成方法生成的OFDM同步训练序列，包括以下步骤粗定时同步，小数频偏估计，整数频偏估计，细定时同步。其中采用双符号联合判决算法获得较优的粗定时同步估计，又采用余数定理算法对频偏进行估计。相对于现有方案，本发明具有计算复杂度较低、频偏估计范围较大、符号定时同步及频率同步较为准确等优点。
本发明是通过以下技术方案实现的本发明所述的OFDM同步训练序列的生成方法，是由2个OFDM符号级联生成OFDM同步训练序列，而现有的一些同步训练序列结构往往需要更多符号(如无线局域网需要4个OFDM符号)，因此，用本发明的OFDM同步训练序列的生成方法所生成的训练序列具有数据开销低的优点。该序列的构成如下第1个符号这是一个非正则符号(irregular symbol)，记为Sa，Sa的有效符号部分的点数为Na(Na为整数，且含有奇素数因子q，Na≠N，N为除非正则符号外，包括正则符号在内的其他OFDM符号的有效符号点数)，其被等分为La份(La为整数，也含有奇素数因子q)，通常，每一份信号称为一个slot。每一个slot含Ma点(Ma为整数，且Ma＝Na/La)，Sa的循环前缀部分含点数为Nag(Nag＝Nsym-Na)，因此，Sa的总点数仍为Nsym。
第2个符号这是一个正则符号(regular symbol)，记为Sn，其有效部分的点数为N，其被等分为L份(L为整数，且L＝2v)，每一个slot含M点(M为整数，且M＝N/L)，L的一些典型值为4、8、16等。
以下对本发明的OFDM同步训练序列的生成方法进行进一步的描述设Sa、Sn对应的OFDM符号标号分别为1、2。
Sn是一个L等分的正则OFDM同步训练符号，这是一种相当常见的序列结构，经常用于符号定时同步算法和频偏估计算法中。例如，IEEE802.11a标准中取L＝4。其生成方法是首先按照式(1)插入频域数据a2,k≠0k(modL)=0a2,k=0k(modL)≠0(0<|k|Nu2)---(1)]]>然后进行N点傅里叶反变换产生对应的时域信号，之后在时域信号的前端加上循环前缀，即可生成Sn。在L等分的情况下，Sn的每一个slot所含的点数为M(M＝N/ L)。
符号Sa含有La个slot，设La含有奇素数因子q，且有，2α-1＜La≤2α，(α为一个正整数)，Na(mod La)＝0以及N(mod La)≠0，取Na为La的整数倍数中与N最接近的自然数，即Na=argminμ{|μ-N||μ(modLa)=0}---(2)]]>Sa的生成方法是将La和Nu-2α+1代替式(1)中的L和Nu，得a1,k≠0k(modLa)=0a1,k=0k(modLa)≠0(0<|k|≤Nμ2α+12)---(3)]]>按照式(3)插入频域数据，然后进行Na点傅里叶反变换后可得b1，lb1,l=1NΣk=0Na-1aI,kexp(j2πklNa)(l=0,1,2,...,Na-1)---(4)]]>
再添加Nag点循环前缀，就能得到Sa。需要说明的是，式(3)的Nu，比式(1)少2α+1个，这是为了补偿式(2)造成的频谱畸变。
本发明还提供基于上述同步训练序列的同步方法，具体如下步骤一粗定时同步，对Sa作延时为Ma点，窗长为(Nsym-Ma)点的自相关，得Pa(d)。对Sn作延时为M点，窗长为(Nsym- M)点的自相关，得Pn(d)。对Sn作延时为N/2点，窗长为(Nsym-N/2)点的自相关，得Pf(d)，再对Pa(d)、Pn(d)与Pf(d)求能量之和，除以相关窗内的信号总能量，得平均自相关P(d)，寻找P(d)的峰值，初步估计出OFDM符号的起始位置 这一算法称为“双符号联合判决算法”。该算法可以大大提高系统粗定时同步的性能，而且，该算法具有较低的计算复杂度。
步骤二小数频偏估计，即根据步骤一得到的Pf(d)，计算其在位置 处的相角，从而估计频偏的小数部分 和整数频偏的奇偶性η。
步骤三整数频偏估计，即根据步骤一得到的Pa(d)和Pn(d)，计算其在位置 处的相角，利用数论中的余数定理，估计整数频偏 该算法称为“余数定理算法”。将步骤二得到的 与 相加，得到频偏估计值 其估计范围可以达到[-Lx/2，Lx/2)(Lx是L和La的最小公倍数)。该算法同时具有计算复杂度低和频偏估计范围大的优点。
步骤四细定时同步，即根据步骤一得到的 抽取Sn的第2个slot，然后根据步骤三得到的 对其进行频偏补偿，再把其与本地slot样本作循环移位相关，寻找相关峰，从而估计出OFDM符号的精确起始位置 以下对本发明的同步方法进行进一步的描述(1)粗定时同步本发明提供的粗定时同步算法称为“双符号联合判决法”。
首先，计算Pa(d)、Pn(d)和Pf(d)。如下对Sa作延时为Ma点，窗长为(Nsym-Ma)点的自相关，得Pa(d)
Pa(d)=Σm=dd+Nsym-Mar1,m*r1,m+Ma---(5)]]>对Sn作延时为M点，窗长为(Nsym-M)点的自相关，得Pn(d)Pn(d)=Σm=dd+Nsym-Mr2,m*r2,m+M---(6)]]>对Sn作延时为N/2点，窗长为(Nsym-N/2)点的自相关，得Pf(d)Pf(d)=Σm=dd+Nsym-N/2r2,m*r2,m+N/2---(7)]]>接着对式(5)、式(6)与式(7)求能量之和，再除以相关窗内的信号总能量，得平均自相关P(d)，P(d)=|Pa(d)|2+|Pn(d)|2+|Pf(d)|2(Σm=dd+Nsym-Ma|r1,m+Ma|2)2+(Σm=dd+Nsym-M|r2,m+M|2)2+(Σm=dd+Nsym-N/2|r2,m+N/2|2)2---(8)]]>然后寻找P(d)的峰值，从而估计出OFDM符号的粗定时起始点 d^=argmaxd{P(d)}---(9)]]>式(8)中，|Pa(d)|2和|Pn(d)|2的相关峰不尖锐，在低信噪比时具有较强的抗噪性能；|Pf(d)|2的相关峰较尖锐，在高信噪比时可以对符号定时同步作出较为精确的估计。因此，本发明提供的“双符号联合判决法”使用多个相关峰合并，大大提高了系统粗定时同步的性能。而且，式(5)至式(8)的计算可以递归进行，对每个采样点的计算量等效为9次复乘和若干次复数加减。
需要指出的是，“双符号联合判决法”只适用于本发明设计的同步序列结构。如果采用另一个正则符号Sn代替Sa，假设其也被等分为L份，那么，“双符号联合判决法”会在d＝Nsym和d＝-Nsym(假定理想的符号起始位置d＝0)处分别出现了一个相关峰，这会严重干扰符号定时同步的判断。
(2)小数频偏估计本发明采用的小数频偏算法是一种传统的小数频偏估计算法。该方法可参见文献Schmidl，T.M.等“Low-overhead，low-complexity[burst]synchronization for OFDM”，IEEE International Conference onCommunications，Volume 3，June 1996，Page(s)1301-1306(“低数据开销、低复杂度的OFDM同少方法”，IEEE国际通信技术会议)。具体实现方法是把N点OFDM符号等分为2份，对其作自相关，再求相角，从而估计小数频偏。
由于本发明中的Sn是一个正则符号，L为偶数，故Sn也可以视作2等分的OFDM符号，时以利用式(7)中Pf(d)的相角来估计小数频偏ξ=1πangle{Pf(d^)}---(10)]]>ξ的范围是[-1，1)，而小数频偏的范围是[-0.5，0.5)，所以需要分类讨论，得到小数频偏估计值 和整数频偏估计值 的奇偶性η 其中，η＝1表示 为奇数；η＝0表示 为偶数。在本方案中， 和η是重要的辅助信息，该信息帮助整数频偏估计模块实现大范围的频偏估计。
(3)整数频偏估计在未来的OFDM系统中，载频将高达1GHz至100GHz，而子载波的间隔相对较小，只有几kHz至几十kHz，当接收端的晶振不是非常稳定时，就要求OFDM系统的频率同步模块能够纠正较大范围的频偏。以载频为10GHz，晶振频偏为20ppm，子载波间隔为5KHz的系统为例，其最大可能fΔ为±40。目前，纠正大范围频偏的技术主要有2种●基于特定频域序列的同步训练序列，该方法需要对同步训练符号进行小数频偏补偿，再作FFT，至少需要 次复数乘法，然后与已知频域序列作循环移位相关，通过寻找相关峰来估计整数频偏。该方法所能估计的频偏范围，随着相关搜索范围的扩大而增大，但是，其计算复杂度也会不断上升。每搜索一个频点，就需要约 次复乘。因此，由于计算复杂度的原因，该方法难以应用到实际系统中。参见文献Schmidl，T.M.等“Low-overhead，low-complexity[burst]synchronization for OFDM”，IEEE International Conference onCommunications，Volume 3，June 1996，Page(s)1301-1306(“低数据开销、低复杂度的OFDM同步方法”IEEE国际通信技术会议)。
●基于L等分(L＝2v)的OFDM 同步训练符号结构，通过计算训练符号的特定延时的自相关，再求相角来估计整数频偏。该方法的最大估计范围为[-L/2，L/2)，但是，L的实际取值一般不能太大。因为过大的L，会造成OFDM系统的峰均功率比过高，并产生较大的频偏估计误差，以及训练符号的频域数据稀疏性导致其抵抗频率选择性衰落的能力下降。参见文献Heiskala J，Terry JOFDM Wireless LANs-A Theoretical and Practical Guide.[M].Indianapolis USAPearson Education Inc，2002.70-73(《OFDM无线局域网——理论与实践的指导》)。
基于本发明的OFDM同步训练序列的生成方法生成的OFDM同步训练序列，设计一种频偏估计范围较大，且复杂度较低的整数频偏估计算法，称之为“余数定理算法”。该算法首先利用式(11)得到的 和η，求出式(5)中的Pa(d)和式(6)中的Pn(d)在位置 处的相角，从而对整数频偏作出2个估计(fa与fn)fa=argminλ{|λ-(La2πangle{Pa(d^)}-f^F)|}La(mod2)=1argminλ{|λ-(La2πangle{Pa(d^)}-f^F)||λ(mod2)=η}La(mod2)=0---(12)]]>fn=argminλ{|λ-(La2πangle{Pn(d^)}-f^F)||λ(mod2)=η}---(13)]]>fa(mod La)有La种可能的取值{0，1，...，La-1}；fn(mod L)有L种可能的取值{0，1，...，L-1}。由fa(mod La)和fn(mod L)，可以导出一个不定方程组f^I(modLa)=fa(modLa)f^I(modL)=fn(modL)---(14)]]>根据数论中的余数定理，由fa(mod La)和fn(mod L)的不同组合，可以解出 其有Lx种可能的整数解{-Lx/2，-Lx/2+1，...，Lx/2-1}(Lx是L和La的最小公倍数)。因此，本发明所能估计的频偏范围达到[-Lx/2，Lx/2)。为了方便硬件实现，可以使用查表法来解方程组(14)先对Lx个可能解 进行枚举，计算出其对应的
和
并把这两个值作为
的表格的入口地址，再将
填进入口地址对应的存储单元，这样就制成了方程组(14)的求解表。在估计整数频偏过程中，只需用fa(mod La)和fn(mod L)作为表格索引，查表求出
即可。
本发明提供的“余数定理算法”同时具有计算复杂度低和频偏估计范围大的优点，具体来讲，本算法的式(12)和式(13)只需要2次求相角的运算，然后可以用查表法解不定方程组(14)，就能完成对整数频偏的估计f^I=LUT(fa(modLa),fn(modL))---(15)]]>式(15)中的函数LUT代表查表函数，其计算复杂度可以忽略不计，然而，本算法的频偏估计范围却非常大。以La＝7，L＝8为例，由于Lx＝56，所以频偏估计范围为[-28，28)，再以La＝15，L＝8为例，由于Lx＝120，所以频偏估计范围达到[-60，60)，这样大的频偏估计范围可以满足未来OFDM系统的需要。
另外，本发明具有灵活性，选取不同的La和L，就能满足不同场合的整数频偏估计的要求。比如，取La＝3，L＝2，则Lx＝6，其频偏估计范围只有[-3，3)，可以用于小范围的频偏估计中。而且，方程组(14)所用的求解表是相当简单的，只需一个大小为(La×L)Byte的二维数组来存储Lx个整数解即可。以La＝7，L＝8为例，方程组(14)的求解表如下表所示
将式(11)得到的
与式(15)得到的
相加，作出完整的频偏估计
f^Δ=f^I+f^F---(16)]]>需要特别指出的是，由于Sa的非正则性，式(12)中的 是频偏估计的近似公式。实际上，利用 的相角进行频偏作估计的正确公式如式(17)所示ξa=N2π(Na/La)angle{Pa(d^)}---(17)]]>当式(2)确定的Na与N非常接近时ξa≈La2πangle{Pa(d^)}---(18)]]>显然，式(18)的估计范围是[-La/2，La/2)，满足余数定理的要求，而式(17)的估计范围是[-NLa/2Na，NLa/2Na)，与余数定理的内在要求不符，因此，本发明中，采用式(18)来代替式(17)，由此引入的误差因子为error(ξa)=|NaN-1|---(19)]]>由于error(ξa)是一个乘性因子，所以其在频偏较大时的影响会略微严重一些，其造成的最大频偏误差为 例如，当Na＝1022，La＝7，N＝1024，L＝8时，error(ξa)＝0.00195，La2error(ξa)=0.006825.]]>考虑到ξa只是用于整数频偏估计，因此，error(ξa)几乎不会对fa的估计产生影响。
(4)细定时同步为了达到较好的定时同步估计性能，本发明在采用“双符号联合判决法”取得较好的粗定时同步估计 后，再采用细同步算法对符号定时同步作出精细估计，具体方法为抽取正则符号Sn中的第2个slot，并用式(16)得到的 对其进行相位补偿，得x＝(x0，x1，...，xm，...，xM-1)，其中，xm=r2,d^+M+mexp(-j2πmf^Δ/N)(m=0,1,2,...,M-1)---(20)]]>
将x与本地slot样本作时域的循环移位相关，寻找相关峰，从而达到细定时同步的目的ϵ^=d^+argmaxu{R(u)=|Σm=0M-1b2,m+u(modM)*xm|2|u=-M2,-M2+1,...,M2-1}---(21)]]>只要
的估计误差在±M/2以内，式(21)可以对OFDM符号的起始时刻作出相当准确的估计，其主要计算量是M(M+1)次复数乘法。
本发明的创新之处在于设计了非正则OFDM符号Sa与正则OFDM符号Sn级联的同步训练序列的生成方法，并基于该训练序列，提出“双符号联合判决法”，不会出现2个相关峰，防止了对符号定时同步的判断的干扰，改善了粗定时同步的估计性能，同时又提出“余数定理算法”，将频偏估计范围扩大至[-Lx/2，Lx/2)，Lx是Sa的等分数La与Sn等分数L的乘积，一些常用的La与L的组合参数及其频偏估计范围和误差如下表所示
特别需要指出的是，本发明的计算复杂度非常低，如下表所示
而且，本发明的符号定时同步错误率和频偏估计错误率较低，因此，本发明在OFDM系统中具有很高的应用价值。
图1为OFDM基带调制解调框2为本发明同步训练序列的示意3为Na＝1022，La＝7时，符号Sa的生成示意4为N＝1024，L＝8时，符号Sn的生成示意5为本发明的同步方法框6为本发明的同步方法的具体实施流程7为本发明提供的“双符号联合判决算法”在Na＝1022，La＝7，N＝1024，L＝8，Nsym＝1088时，相关运算及合并的示意8为本发明提供的“双符号联合判决算法”在Na＝1024，La＝8，N＝1024，L＝8，Nsym＝1088时，相关运算及合并的示意9为本发明在0dB高斯白噪声情况下，某一次仿真时P(d)的数值图形图10为本发明在0dB高斯白噪声情况下，某一次仿真时angle{Pf(d)}/π的数值图形图11为本发明在0dB高斯白噪声情况下，某一次仿真时R(u)的数值图形图12为本发明与传统的正则符号训练序列在粗定时同步方面的性能比较图13为本发明与传统的正则符号训练序列在定时同步估计标准差方面的性能比较图14为本发明与传统的正则符号训练序列在整数频偏估计错误率方面的性能比较。
图15为本发明与频域估计算法在整数频偏估计错误率方面的性能比较图16为SNR＝2dB，采用本发明的整数频偏估计算法在估计范围内的错误率平坦性具体实施方式
下面给出一个具体的OFDM参数配置实施例，来阐述本发明的实现步骤。需要说明的是，下例中的参数并不影响本发明的一般性。
3GPP组织的文档TR 25.892 V6.0.0，“Feasibility Study for OrthogonalFrequency Division Multiplexing(OFDM)for UTRAN enhancement (Release 6)”(通用移动通信系统及陆基无线电接入的技术升级中采用OFDM方案的可靠性研究(第6版))，给出的一组OFDM参数，如下载频fc 3.5GHz系统带宽B 6.528MHz子载波数N 1024有效子载波数N 705有效带宽 4.495MHz子载波间隔Δf 6.375kHz循环扩展CP 64(9.803us)符号周期Ts 156.85+9.81＝166.66us保护载波数 319在上述参数条件下，取Na＝1024，La＝7，N＝1024，L＝8，Nsym＝1088，仿真时，设系统频偏fΔ＝19.3，采用的信道为8路瑞利衰落信道
本发明的实现步骤如下(1)生成OFDM同步训练符号，图2为由本发明的OFDM同步训练序列的生成方法所生成的同步训练序列的示意图。其具体生成过程如下依照式(2)得Na＝1022，根据式(3)在子载波序号为170、177、184、191、198、205、212、219、226、233、240、247、254、261、268、275、282、289、296、303、310、317、324、331、338、345、352、359、366、373、380、387、394、401、408、415、422、429、436、443、450、457、464、471、478、485、492、499、506、520、527、534、541、548、555、562、569、576、583、590、597、604、611、618、625、632、639、646、653、660、667、674、681、688、695、702、709、716、723、730、737、744、751、758、765、772、779、786、793、800、807、814、821、828、835、842、849、856的频率上加载随机生成的恒包络频域数据，在其他子载波上放置零数据，将这些频域数据通过图1的发送端串并转换、IDFT、发送端并串转换、插入循环前缀等模块，得到式(4)所对应的7等分的同步训练符号Sa。图3是Na＝1022，La＝7时，Sa的生成示意图。然后，依照式(1)在子载波序号为160、168、176、184、192、200、208、216、224、232、240、248、256、264、272、280、288、296、304、312、320、328、336、344、352、360、368、376、384、392、400、408、416、424、432、440、448、456、464、472、480、488、496、504、520、528、536、544、552、560、568、576、584、592、600、608、616、624、632、640、648、656、664、672、680、688、696、704、712、720、728、736、744、752、760、768、776、784、792、800、808、816、824、832、840、848、856、864的频率上加载随机生成的恒包络频域数据，其他子载波上放置零数据，将这样的频域数据通过图1的发送端串并转换、IDFT、发送端并串转换、插入循环前缀等模块，就能生成8等分的同步训练符号Sn。图4是N＝1024，L＝8时，Sn的生成示意图。
(2)运用“双符号联合判决法”进行粗定时同步，该过程位于图5的粗定时同步模块，其详细过程位于图6的粗定时同步模块，即根据式(5)至式(8)，估计出粗定时位置 图7是“双符号联合判决算法”在Na＝1022，La＝7，N＝1024，L＝8，Nsym＝1088时，相关运算及合并的示意图(假定理想的符号起始位置d＝0)。|Pa(d)|2和|Pn(d)|2的相关峰不尖锐，在低信噪比时具有较强的抗噪性能；|Pf(d)|2的相关峰较尖锐，在高信噪比时可以对符号定时同步作出较为精确的估计。因此，本发明提供的“双符号联合判决法”使用多个相关峰合并，大大提高了系统粗定时同步的性能。如果采用另一个正则符号Sn代替Sa，假设其也被等分为L＝8份，那么，“双符号联合判决法”的相关运算及合并的示意图如图8所示，此时，在d＝Nsym和d＝-Nsym处分别出现了一个相关峰，这会严重干扰符号定时同步的判断。图9是本发明在0dB高斯白噪声情况下，某一次仿真时P(d)的数值图形，可见其在正确的定时点附近出现了峰值。在图9中，d＝0是理想的符号起始点，此时根据式(8)得到的的粗同步估计值为d^=-5.]]>(3)小数频偏估计，该过程位于图5的小数频偏估计模块，其详细过程位于图6的小数频偏估计模块，即根据式(10)和式(11)，估计出小数频偏 和整数频偏估计值 的奇偶性η。图10是本发明在0dB高斯白噪声情况下，某一次仿真时angle{Pf(d)}/π的数值图形，由式(10)计算得ξ=angle{Pf(d^)}/π=angle{Pf(-5)}/π=-0.6772]]>再由式(11)得f^F=ξ-ξ/|ξ|=-0.6772-(-0.6772)/abs(-0.6772)=0.3228]]>η{(ξ-f^f)(mod2)}={(-0.6772-0.3228)(mod2)}=1]]>(4)运用“余数定理算法”进行整数频偏估计，该过程位于图5的整数频偏估计模块，其详细过程位于图6的整数频偏估计模块，即根据式(12)和式(13)，得到2个整数频偏fa和fn，代入式(14)，把fa(mod 7)和fn(mod 8)作为入口地址，查询不定方程的求解表，从而对整数频偏作出估计 再由式(16)得到完整的频偏估计值 本发明在0dB高斯白噪声情况下，某一次仿真时，d^=-5,]]>f^F=0.3228,]]>η＝1，由式(12)得fa＝5，由式(13)得fn＝3，查询不定方程的求解表，在纵坐标为5，横坐标为3处得到f^I=19,]]>再由式(16)得f^Δ=f^I+f^F=19.3228,]]>频偏仙计误差仅为19.3228-19.3＝0.0228。
(5)细定时同步，该过程位于图5的细定时同步模块，其详细过程位于图6的细定时同步模块，即根据 的信息，抽取正则符号Sn中的第2个slot，运用式(20)对其进行频偏补偿，再通过式(21)，对符号定时位置作出精细估计 本发明在0dB高斯白噪声情况下，某一次仿真时，d^=-5,]]>f^Δ=19.3228,]]>运用式(20)对接收到的Sn中第2个slot进行频偏补偿后，再计算式(21)中的R(u)，R(u)的数值图形如图11所示，其在u＝5出现峰值，所以，根据式(21)得ϵ^=d^+argmaxu{R(u)}=-5+5=0]]>可见，此时的定时估计误差为0。
以下考察本发明的性能表现图12是本发明与传统的正则符号训练序列在粗定时同步方面的性能比较。传统的正则符号训练序列的参数为Na＝1024，La＝64，N＝1024，L＝64，Nsym＝1088。定义粗定时同步错误率为概率Pr(|d^|≤M2),]]>M＝128。图12表明，本发明的粗定时同步性能优于传统方案超过4dB，这是本发明的“双符号联合判决算法”所产生的效果。
图13是本发明与传统的正则符号训练序列在定时同步估计标准差方面的性能比较。传统的正则符号训练序列的参数仍为Na＝1024，La＝64，N＝1024，L＝64，Nsym＝1088。其细同步也采用与本发明相同的细同步方法。图13表明，本发明的定时同步性能优于传统方案超过3dB，这主要是本发明的“双符号联合判决算法”使粗同步较为准确所产生的效果。
图14是本发明与传统的正则符号训练序列在整数频偏估计错误率方面的性能比较。传统的正则符号训练序列的参数仍为Na＝1024，La＝64，N＝1024，L＝64，Nsym＝1088，其小数频偏估计算法与本发明完全相同，其整数频偏估计算法也与本发明类似，因而计算复杂度也基本相同，不同之处是传统的正则符号训练序列无法采用“余数定理算法”，其通过式(12)和式(13)得到fa与fn之后，形成式(14)的方程组，然后直接将fa与fn作平均，得到整数频偏的估计值。图14表明，本发明的整数频偏估计性能优于传统方案超过6dB，这主要是本发明的“余数定理算法”使整数频偏估计较为准确。
图15是本发明与频域估计算法在整数频偏估计错误率方面的性能比较。频域估计算法的具体实现方法参见文献Schmidl，T.M.等“Low-overhead，low-complexity [burst] synchronization for OFDM”，IEEE InternationalConference on Communications，Volume 3，June 1996，Page(s)1301-1306(“低数据开销、低复杂度的OFDM同步方法”IEEE国际通信技术会议)。由于只考察频偏估计性能，所以，令两种方法的时间同步完全正确，且该频域估计算法中获得小数频偏估计值的方法以及频偏估计范围与本发明完全相同。图15表明，本发明的整数频偏估计性能几乎与频域估计算法相同，而本发明的计算复杂度远远低于该频域估计算法，在当前的参数条件下，本发明的计算复杂度大约只有法频域估计算法的5％。
图16是本发明在SNR＝2dB时，在频偏估计范围内的错误率平坦性仿真结果，该图表明，尽管非正则符号Sa引入了估计误差error(ξa)，但其对整数频偏估计几乎没有影响。
仿真结果表明，本发明具有计算复杂度较低，且估计错误率较低的优点，在OFDM系统中具有很高的应用价值。
权利要求
1.一种OFDM同步训练序列的生成方法，其特征在于由2个OFDM符号级联生成OFDM同步训练序列，所述2个OFDM符号分别是非正则OFDM符号Sa与正则OFDM符号Sn，所述的非正则符号Sa的有效符号部分的点数为Na，Na≠N，Na含有奇素数因子q，所述的非正则符号Sa被分为La等分，La含有奇素数因子q，所述每一等分含Ma点，其中Ma＝Na/La，且Na、La、Ma均为整数，Sa的循环前缀部分含点数为Nag，Nag＝Nsym-Na，Nsym是OFDM符号的总点数；所述正则符号Sn被分为L等分，L＝2v，所述每一个等分含M点，其中M＝N/L，N为除非正则符号外，包括正则符号在内的其他OFDM符号的有效符号点数，且N、L、M均为整数。
2.根据权利要求1所述的OFDM同步训练序列的生成方法，其特征在于所述非正则符号Sa，其生成方法是，选取Na为La的整数倍数中与N最接近的自然数Na=argminμ{|μ-N||μ(modLa)=0}]]>然后，在频域上每间隔La个子载波插入数据，进行Na点傅里叶反变换，再添加Nag点循环前缀。
3.一种基于根据权利要求1所述的方法生成的OFDM同步训练序列的同步方法，其特征在于，依次通过粗定时同步，小数频偏估计，整数频偏估计，细定时同步，来完成对OFDM符号的精确起始位置的估计，包括如下步骤步骤一对非正则符号Sa作延时为一个等分信号长度的自相关，得Pa(d)；对正则符号Sn作延时为一个等分信号长度的自相关，得Pn(d)；对正则符号Sn作延时为半个有效符号长度的自相关，得Pf(d)，再对Pa(d)、Pn(d)与Pf(d)求能量之和，除以相关窗内的信号总能量，得平均自相关P(d)，寻找P(d)的峰值，初步估计出OFDM符号的起始位置 步骤二根据步骤一得到的Pf(d)和 ，计算 的相角，从而估计频偏的小数部分 和整数频偏的奇偶性η；步骤三根据步骤一得到的Pa(d)，Pn(d)和 ，计算 和 的相角，再结合步骤二得到的 和和η，估计整数频偏 步骤四根据步骤一得到的 ，抽取正则符号Sn的第2个等分信号，然后对其进行频偏补偿，再把其与本地正则符号Sn的第2个等分信号样本作循环移位相关，寻找相关峰，从而估计出OFDM符号的精确起始位置
4.根据权利要求3所述的同步方法，其特征是，所述步骤一，对Sa作延时为Ma点，窗长为(Nsym-Ma)点的自相关，得Pa(d)；对Sn作延时为M点，窗长为(Nsym-M)点的自相关，得Pn(d)；对Sn作延时为LM/2点，窗长为(Nsym-LM/2)点的自相关，得Pf(d)。
5.根据权利要求3所述的同步方法，其特征是，所述步骤一，对Pa(d)、Pn(d)与Pf(d)求能量之和，再除以相关窗内的信号总能量，得平均自相关P(d)，P(d)=|Pa(d)|2+|Pn(d)|2+|Pf(d)|2(Σm=dd+Nsym-Ma|r1,m+Ma|2)2+(Σm=dd+Nsym-M|r2,m+M|2)2+(Σm=dd+Nsym-N/2|r2,m+N/2|2)2,]]>其中，ri，m表示第i个符号，第m个采样点的数据。
6.根据权利要求3所述的同步方法，其特征是，所述步骤三，估计整数频偏 时，先由 对整数频偏作出估计fa，fa=argminλ{|λ-(La2πangle{Pa(d^)}-f^F)|}La(mod2)=1argminλ{|λ-(La2πangle{Pa(d^)}-f^F)||λ(mod2)=η}La(mod2)=0]]>再由 对整数频偏作出估计fn，fn=argminλ{|λ-(La2πangle{Pn(d^)}-f^F)||λ(mod2)=η}]]>得到关于整数频偏的不定方程组f^1(modLa)=fa(modLa)f^I(modL)=fn(modL)]]>最后，运用数论中的余数定理，解这个不定方程组，从而估计出整数频偏
全文摘要
本发明提供一种OFDM同步训练序列的生成方法和基于该训练序列的同步方法，用于通信技术领域。OFDM同步训练序列的生成方法为由2个OFDM符号级联生成OFDM同步训练序列，所述2个OFDM符号分别是非正则OFDM符号S
文档编号H04L7/00GK101014029SQ20061003033
公开日2007年8月8日 申请日期2006年8月24日 优先权日2006年8月24日
发明者丁铭, 吴赟, 罗汉文, 佘锋, 张霆蔚, 张海滨 申请人:上海交通大学, 夏普株式会社