基于ZC序列的周期前导序列生成方法
技术领域
本发明涉及一种基于ZC(Zadoff-Chu)序列的周期前导序列生成方法,属于无线通信技术领域。
背景技术
在无线通信***中前导序列的作用主要是利用前导序列的相关性完成帧检测、自动增益控制(AGC,Automatic Gain Control)、粗时频同步、细时频同步和信道估计等。同时,与单载波***相比,正交频分复用(OFDM,Orthogonal Frequency DivisionMultiplexing)***的峰均功率比(PAPR,Peak to Average Power Ratio)较高,会显著降低发射机功率放大器的效率。因此,前导序列的设计需要遵循具有良好相关特性和低PAPR两个原则。
长期演进(LTE,Long Term Evolution)通信***的主同步信号(PSS,PrimarySynchronization Signal)是由满足恒幅零自相关(CAZAC,Constant Amplitude ZeroAutocorrelation)特性的ZC序列组成。它具有较低的PAPR,同时也具有理想的周期自相关性以及良好的周期互相关性、自相关性和互相关性,适合于设计前导序列。
前导序列的周期特性可以用来进行帧检测,即将连续两个周期长度的序列进行自相关运算,当出现峰值时即可判断为检测到数据帧。由于数据帧到达时间是未知及不同的信噪比(SNR,Signal to Noise Ratio)下的自相关峰值浮动,难以利用峰值直接检测是否有数据帧到来。为了解决这个问题,很多研究成果提出在不同的SNR下设计不同的阈值,当在某个特定SNR下自相关值大于对应的阈值时,则判断为数据帧的到来。显然,这里存在的问题是假设帧检测时已经知道SNR,这与实际应用不符合。同时由于循环前缀(CP,CyclicPrefix)的存在,自相关峰值会出现平台现象,很多研究成果没有利用它进行帧检测,导致性能受到阈值的影响较大。
发明内容
发明目的:针对现有技术中存在的问题与不足,本发明提供一种基于ZC序列的周期前导序列生成方法。采用满足CAZAC特性的ZC序列生成前导序列的频域信号,并进行升采样技术以得到周期性时域序列。发射节点通过空间映射模块将它映射到多个天线端口,再进行资源映射并生成OFDM符号,最终形成周期前导序列进行发送。任意SNR条件下,接收节点的每根天线利用前导序列的周期性求解自相关值与同一个阈值比较,并结合自相关峰值多峰平台现象判断是否有数据帧到来,再将所有接收天线的帧检测结果进行集中处理得到最终帧检测结果。此发明可确保阈值的设置不受SNR影响,同时显著提高帧检测的可靠性。
为了更清楚的理解本发明内容,首先介绍一下本发明所涉及的ZC序列的相关知识:
ZC序列是非二进制单位振幅序列,长度为NZC的ZC序列zμ的第k个元素zμ(k)可以表示为
其中j是虚数单位,μ是ZC序列的根指数(μ与NZC互质且1≤μ≤NZC-1),q可以是任意整数。为了简单起见通常设置q=0。
ZC序列具有以下特性:
(1)ZC序列具有恒定振幅1,限制了PAPR和对其他用户产生的边界和时间平坦性干扰。
(2)ZC序列具有理想的周期自相关性以及良好的周期互相关性、自相关性和互相关性。
(3)ZC序列经过离散傅里叶变换(DFT,Discrete Fourier Transform)或逆离散傅里叶变换(IDFT,Inverse Discrete Fourier Transform)后仍然是ZC序列。
技术方案:一种基于ZC序列的周期前导序列生成方法,在发射端,包括以下步骤:
(1)生成前导序列的频域信号。前导序列包括NP个快速傅里叶变换(FFT,FastFourier Transform)点数为NFFT的相同OFDM符号,每个OFDM符号所对应的时域序列由T(T为偶数)个周期长度为的序列和长度为NCP的CP构成。令所使用的频域ZC序列zμ的长度NZC为素数,每个OFDM符号的直流子载波数为NDC(NDC为奇数)、数据子载波总数为NST(包括直流子载波,NST为奇数),他们之间的关系满足NST-TNZC<24。
为了避免直流子载波,将zμ的中间NDC个元素置为0;然后对改变后的zμ进行升T倍采样,即在每个元素后面***T-1个0;接着将升T倍采样结果补零至NST长度,具体操作是zμ(0)之前补个0,zμ(NZC-1)后面的T-1个0之后补个0,得到NST×1维向量它的第k个元素表示为
考虑到直流子载波,的中间NDC个元素不使用,因此将的元素进行重新编号,生成序列x,它的第k个元素表示为
(2)空间映射。将前导序列的频域信号空间映射到NTx根发射天线上,即:
s(:,k)=p(k)x(k),k=0,1,…,NST-NDC-1 (4)
其中是第k个元素x(k)的NTx×1维空间映射向量,δu表示第u根发射天线的映射系数,u=1,2,...,NTx,δ1=0;s(:,k)是空间映射输出的NTx×1维向量,它的第u个元素对应x(k)在第u根发射天线上的空间映射结果。
(3)资源映射。将空间映射输出结果资源映射到对应子载波上。前导序列的资源映射描述如下:
r(:,c(k))=s(:,k),k=0,1,…,NST-NDC-1 (4)
其中c(k)是向量的第k个元素,即子载波序号;r(:,c(k))是资源映射输出的第c(k)个子载波上的NTx×1维向量,它的第u个元素的值为s(u,k)。
(4)对每根发射天线上的资源映射结果做NFFT点逆快速傅里叶变换(IFFT,InverseFast Fourier Transform)得到周期时域序列,***CP后组成OFDM符号,并将该OFDM符号重复NP次发送。
在任意SNR下,接收节点的每根天线利用前导序列的周期性进行自相关运算且与同一个阈值比较,并结合NP个自相关峰值平台现象判断数据帧的到来,再将所有接收天线的帧检测结果进行集中处理得到最终帧检测结果。
每根接收天线的帧检测过程如下:
(1)生成信号采样时刻计数器n,其初始值为0;
(2)对第v根天线的接收采样信号yv(n)从n采样时刻起连续两个长度的序列进行自相关运算,v=1,2,...,NRx,NRx是接收天线数,得到在n采样时刻处的判决变量:
其中表示yv(n+m)的共轭。
(3)将M(n)与阈值λ进行比较,0<λ≤0.05是一个经验值且在所有接收天线的任意SNR下都是相同的。如果M(n)>λ,执行步骤(4);否则n=n+1并执行步骤(2);
(4)生成迭代计数器t,其初始值为0;
(5)如果t≤Niter,其中Niter是最大迭代次数,令t=t+1并执行步骤(6);否则提示未检测到数据帧。
(6)令为中间变量并进行赋值:
(7)令为均值变量并进行赋值:
(8)如果执行步骤(9);否则n=n+NCP并执行步骤(5);
(9)第v根接收天线检测到数据帧到来且帧起始位置nv=n-NCP。
得到NRx根接收天线的帧检测结果后进行集中处理。如果有多于根接收天线检测到数据帧到来,则最终帧检测结果是检测到数据帧且最终帧起始位置是这些接收天线的帧起始位置的平均值;否则未检测到数据帧,***等待下一次帧检测。
有益效果:与现有技术相比,本发明提供的基于ZC序列的周期前导序列生成方法充分利用ZC序列的低PAPR和良好相关特性并通过空间映射和资源映射等来设计周期前导序列,在自相关值与阈值进行比较并初步判断数据帧到来的情况下进一步利用自相关峰值多峰平台现象进行判断,从而确保阈值的设置不受SNR影响,同时有效提高帧检测的可靠性。
附图说明
图1是实施例中前导序列的结构;
图2是实施例中发射端的具体过程;
图3是实施例中第1根发射天线上周期时域序列的幅值波形图;
图4是实施例中第v根接收天线的帧检测过程。
具体实施方式
下面结合具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
基于ZC序列的周期前导序列生成方法,发射节点首先生成前导序列的频域信号,包括频域ZC序列的置零、升采样、补零;然后进行空间映射,将信号映射到不同的天线端口上;接着进行资源映射,将信号映射到对应的子载波上;最后进行IFFT变换生成周期时域序列,***CP后组成OFDM符号,并重复多次发送;对于任意SNR,接收节点的每根天线设置同一个阈值,利用前导序列的周期性进行自相关运算并与该阈值比较,再检测自相关峰值多峰平台现象进行帧检测,最后将所有接收天线的帧检测结果进行集中处理得到最终帧检测结果。
在本发明实例中,前导序列包括NP=2个FFT点数NFFT=1024的相同OFDM符号,每个OFDM符号所对应的时域序列由T=2个周期长度为512的序列和长度为NCP=80的CP构成,如图1所示。令所使用的频域ZC序列z1的长度NZC=317,根指数μ=1。每个OFDM符号的直流子载波数NDC=3、数据子载波总数NST=643。发射天线数NTx=8,接收天线数NRx=8。
发射端的具体过程如图2所示,描述如下:
(1)生成前导序列的频域信号。为了避免直流子载波,将z1的中间3个元素z1(157),z1(158),z1(159)置为0;然后对改变后的z1进行升2倍采样,即在每个元素后面***1个0;接着将升2倍采样结果补零至643长度,具体操作是z1(0)之前补5个0,z1(642)后面的1个0之后补4个0,得到643×1维向量它的第k个元素表示为
考虑到直流子载波,的中间3个元素不使用,因此将的元素进行重新编号,生成序列x,它的第k个元素表示为
(2)空间映射。将前导序列的频域信号空间映射到8根发射天线上,即:
s(:,k)=p(k)x(k),k=0,1,…,639 (11)
其中p(k)=[1,ej4πk/1024,ej7πk/1024,ej12πk/1024,ej16πk/1024,ej19πk/1024,ej24πk/1024,ej28 πk/1024]T是第k个元素x(k)的8×1维空间映射向量,即第1到8根发射天线的映射系数分别为δ1=0、δ2=2、δ3=3.5、δ4=6、δ5=8、δ6=9.5、δ7=12、δ8=14;s(:,k)是空间映射输出的8×1维向量,它的第u个元素对应x(k)在第u根发射天线上的空间映射结果,u=1,2,...,8。
(3)资源映射。将空间映射输出结果资源映射到对应子载波上。前导序列的资源映射描述如下:
r(:,c(k))=s(:,k),k=0,1,…,639 (12)
其中c(k)是向量c=[-321,-320,…,-3,-2,2,3,…,320,321]的第k个元素,即子载波序号;r(:,c(k))是资源映射输出的第c(k)个子载波上的8×1维向量,它的第u个元素的值为s(u,k)。
(4)对每根发射天线上的资源映射结果做1024点IFFT得到周期时域序列,其在第1根发射天线上幅值波形图如图3所示,***CP后组成OFDM符号,并将该OFDM符号重复2次发送。
在任意SNR下,接收节点的每根天线利用前导序列的周期性进行自相关运算且与同一个阈值比较,并结合2个自相关峰值平台现象判断数据帧的到来,再将8根接收天线的帧检测结果进行集中处理得到最终帧检测结果。
每根接收天线的帧检测过程如图4所示,描述如下:
(1)生成信号采样时刻计数器n,其初始值为0;
(2)对第v根天线的接收采样信号yv(n)从n采样时刻起连续两个512长度的序列进行自相关运算,v=1,2,...,8,得到在n采样时刻处的判决变量:
(3)将M(n)与阈值λ=0.03进行比较。如果M(n)>λ,执行步骤(4);否则n=n+1并执行步骤(2);
(4)生成迭代计数器t,其初始值为0;
(5)如果t≤100,令t=t+1并执行步骤(6);否则提示未检测到数据帧;
(6)令temp1,temp2为中间变量并进行赋值:
(7)令mean1,mean2为均值变量并进行赋值:
(8)如果mean1+mean2<temp1+temp2,执行步骤(9);否则n=n+80并执行步骤(5);
(9)第v根接收天线检测到数据帧到来且帧起始位置是nv=n-80。
得到8根接收天线的帧检测结果后进行集中处理。如果有多于4根接收天线检测到数据帧到来,则最终帧检测结果是检测到数据帧且最终帧起始位置是这些接收天线的帧起始位置的平均值;否则未检测到数据帧,***等待下一次帧检测。