本发明涉及水声通信领域，尤其涉及一种基于零相关带序列的水声通信系统多普勒扩展估计方法。
背景技术：
：带宽和频谱利用率是影响通信系统信息传输速率的两个关键因素，MIMO技术能够通过多根天线同时发射和接收信号而显著提高频谱利用率，同时，OFDM技术将频率选择性衰落信道划分为多个正交子信道，每一个子信道频率平坦衰落，以避免符号间干扰。因此，对于多径干扰严重的水声通信系统，MIMO-OFDM技术具有良好的应用前景。水声信道在物理上可以看成是具有不同时延、不同频移、不同起始角的许多传播路径的总和，其复杂性和多变性限制了水声通信的性能，特别是移动水声通信系统，由于声波在海水中的传播速度仅为1500m/s，远低于空气中电磁波的传播速度(3×108m/s)，收发端的移动导致通信信号在时间上的扩展或压缩远大于无线电通信，从而在接收端将严重恶化接收机的载频跟踪和相位符号同步，导致误码率增大。因此,对于移动水声通信系统,对多普勒效应进行估计和补偿是必不可少的，而进行多普勒补偿的前提是精确的估计出多普勒扩展因子。关于水声通信系统的多普勒估计与补偿，现已有比较丰富的研究成果。其中，使用线性调频(LFM,linearlyfrequencymodulated)信号作为前后同步信号，接收端利用已知的LFM序列与接收信号进行互相关运算，根据前后同步信号输出的峰值间隔与实际信号间隔的差值计算多普勒扩展因子的方法存在两个主要的缺点：接收机需要缓存所有的接收数据才能计算前后同步信号的峰值，不利于实时信号处理；由于接收信号与本地已知信号作互相关运算，无法消除由于收发端晶振频率差异等引起的信号频率偏移的影响。另一种使用具有较高多普勒容忍性的双曲线调频信号(HFM,hyperbolicfrequencymodulation)作为前后同步信号的方法，也不能进行实时的信号处理。使用双PN序列作为训练序列，接收端利用二维搜索算法估计多普勒扩展因子的方法，计算复杂度较高。使用带有循环前缀(CP,cyclicprefix)的正交频分复用(OFDM,orthogonalfrequencydivisionmultiplexing)符号作为训练序列的方法，在低信噪比时估计精度不高。本发明提出使用ZCZ序列作为训练序列，接收端使用多个并行相关器对接收信号进行延时自相关处理，根据最大输出结果的相关器的窗口长度进行多普勒扩展因子估计。ZCZ序列具有良好的自相关特性：主瓣尖锐，旁瓣为0，利于接收端峰值的检测；另外，ZCZ序列具有恒模特性，因而有平坦的频谱响应，适合OFDM系统。同时，该方法在信号实时处理和消除收发端固定频偏方面具有优势：重复的序列结构位于数据帧前端，不需缓存整帧数据即可进行同步和多普勒估计等；两段重复序列受到相同的收发端固定频偏影响，自相关运算时取共轭可消除固定频偏。技术实现要素：技术问题：为了克服水声通信系统中现有的多普勒估计方案估计精度较低的问题，本发明提供一种基于零相关带序列的水声通信系统多普勒扩展估计方法，充分利用ZCZ序列良好的自相关特性，在接收端进行实时的信号处理，同时大大提高了多普勒扩展因子的估计精度。技术方案：为实现上述目的，本发明采用的一种基于零相关带序列的水声通信系统多普勒扩展估计方法，使用具有恒模特性和良好的自相关特性的ZCZ序列作为MIMO-OFDM水声通信系统的训练序列，用以估计信号经过时变多径时延信道而引起的扩展或压缩。所述方案包括如下步骤：1)在发送数据帧前插入两段相同的ZCZ序列，每根发射天线上使用ZCZ序列集合中的不同序列作为训练序列；2)信号经过时变多径时延信道，造成时间上的扩展或压缩以及子载波的频率偏移；3)接收端采用多个并行相关器对接收信号进行延时自相关运算；4)根据最大输出结果的相关器的窗口长度进行多普勒扩展因子估计。所述步骤1)中，数据帧前插入的ZCZ序列是利用交织因子化方法生成的，令ZCZ序列集合表示为则Ψ(L,M,Z)的周期相关函数满足：Rci,cjP(τ)=Eci,cj,τ=0,i=j0,0<|τ|≤Z,i≠j0,0<|τ|≤Z,i=j]]>其中，上标P表示周期，L为ZCZ序列长度，M为ZCZ序列条数，Z为零相关带的长度，ci,cj表示序列元素。所述步骤2)中，时变多径时延信道是：冲激响应函数可以表示为：h(τ,t)=ΣpAp(t)δ(τ-τp(t))]]>其中，下标p表示多径数，Ap(t)是路径增益，τp(t)是路径时延，假定所有路径的多普勒扩展因子a相同，路径时延τp(t)，路径增益Ap(t)在一帧符号持续时间内保持不变，即记为τp,Ap。所述步骤2)中，时间上的扩展或压缩以及子载波的频率偏移是：接收端第j根接收天线上的接收信号可以表示为:y~j(t)=Σi=1Nt(Re{Σk=0K-1si[k]Hkej2π(kT+fc)(1+a)t}+w~(t)),t∈[tm,tm+2T(1+a)]]]>其中，Nt为发射天线数，K为子载波个数，s为发送信号，T为一个OFDM符号的持续时间，fc为载波频率，a为多普勒扩展因子，是加性高斯白噪声，tm为第m帧信号起始时间，信道传输函数Hk的定义为：Hk=ΣpApe-j2π(kT+fc)τp,tm≤t≤tm+2T(1+a)]]>信号在时间上的扩展或压缩，即持续时间从2T变为2T/(1+a)，同时，造成每一个子载波发生了的频率偏移。所述步骤3)中，接收端采用多个并行相关器对接收信号进行延时自相关运算的方法是：Cj(n)=Σk=0Kl-1yj(n+k)·(yj)*(n+k+Kl)]]>其中下标j代表第j根接收天线，Kl表示相关器的窗口长度，y表示接收信号，n表示第n个采样点。能量函数P为：Pj(n)=12Σk=0Kl-1[yj(n+k)·(yj)*(n+k)+yj(n+k+Kl)·(yj)*(n+k+Kl)]]]>所述步骤4)中，最大输出结果的判定方法是：Mn=Σj=1Nr|Cj(n)|2/Σj=1NrPj2(n)]]>Nr为接收天线数，选取所有相关器输出结果Mn中的最大值，记为Mmax。所述步骤4)中，进行多普勒扩展因子估计的方法是：取输出最大判决变量的窗口长度值可以计算出多普勒扩展因子为：K^=KF(1+a^)⇒a^=KF-K^K^]]>其中，KF为发送训练序列的实际长度，为多普勒扩展因子的估计值。相应的，令海水中声速为c,可得收发机的相对运动速度为：v^=ca^]]>有益效果：本发明提供的基于ZCZ序列的水声通信系统多普勒扩展估计方法，在接收端使用延时自相关算法进行多普勒扩展因子估计，不需要缓存整帧数据即可进行多普勒扩展因子估计，便于实时的信号处理；并且由于两段重复序列受到相同的收发端固定频偏影响，自相关运算时取共轭可消除收发端固定频偏带来的影响。本发明所提出的方法在时变多径时延信道下仍能比较精确的进行多普勒扩展因子估计，信道适应性强，稳定度高，估计精度远高于传统的多普勒扩展因子估计方法。附图说明图1：ZCZ序列的自相关性质示意图；图2：ZCZ序列的互相关性质示意图；图3：ZCZ序列的恒模特性示意图；图4为本发明中的发送信号的数据帧结构图；图5为本发明中的接收端多个并行相关器示意图；图6为多普勒扩展因子a＝0.005，多径数path＝3，估计误差随信噪比的变化而变化的仿真曲线，其中，估计误差是指多普勒扩展因子的估计值与实际值之间的偏差。图7为多普勒扩展因子a＝0.005，多径数path＝7，估计误差随信噪比的变化而变化的仿真曲线，图8为信噪比为0dB，多径数目为7的情况下，估计误差随运动速度的变化而变化的仿真曲线。具体实施方式下面结合附图对本发明作更进一步的说明。设MIMO-OFDM系统发射天线数为Nt，接收天线数为Nr，信号采用基于循环前缀(CP,cyclicprefix)的OFDM调制方式，以防止码间干扰，令B为信道带宽，K为子载波个数，则子载波间隔为Δf＝B/K，一个OFDM符号持续时间为T＝1/Δf＝K/B，每个OFDM符号的循环前缀时间长度为Tg。如图1、2，3所示，本发明所采用的训练序列为用交织因子化方法生成的ZCZ序列。令ZCZ序列集合表示为则Ψ(L,M,Z)的周期相关函数满足：Rci,cjP(τ)=Eci,cj,τ=0,i=j0,0<|τ|≤Z,i≠j0,0<|τ|≤Z,i=j]]>其中，上标P表示周期，L为ZCZ序列长度，M为ZCZ序列条数，Z为零相关带的长度，ci,cj表示序列元素。当τ∈[-Z,Z]时，ZCZ序列集合的同一序列具有良好的自相关性质，相关峰尖锐；而不同序列之间的相关结果接近于0，另外，ZCZ序列具有恒模特性，因而有平坦的频谱响应，作为OFDM系统的同步符号是非常合适的。如图4所示，训练序列包含两段相同的ZCZ序列，令s＝[s[0],s[1],…,s[K-1]]T表示一个OFDM符号的频域序列，则第i根发射天线上第m帧数据的基带发射信号可以写成：xi(t)=Σk=0K-1si[k]ej2πkTt,t∈[tm,tm+2T]]]>其中，tm为第m帧信号起始时间。相应的，经频率为fc的载波上变频，得到的带通信号为：x~i(t)=Re{Σk=0K-1si[k]ej2π(kT+fc)t},t∈[tm,tm+2T]]]>时变多径时延水声信道的冲激响应函数可以表示为：h(τ,t)=ΣpAp(t)δ(τ-τp(t))]]>其中，Ap(t)是路径增益，τp(t)是路径时延，我们假定：1)所有路径的多普勒扩展因子a相同，即：τp(t)≈τp-at2)路径时延τp(t)，路径增益Ap(t)在一帧符号持续时间内保持不变，即记为τp,Ap。信号经过上述信道,接收端第j根接收天线上的接收信号可以表示为:y~j(t)=Σi=1Nt(Re{Σk=0K-1si[k]Hkej2π(kT+fc)(1+a)t}+w~(t)),t∈[tm,tm+2T(1+a)]]]>其中，Nt为发射天线数，K为子载波个数，s为发送信号，T为一个OFDM符号的持续时间，fc为载波频率，a为多普勒扩展因子，是加性高斯白噪声，信道传输函数Hk的定义为：Hk=ΣpApe-j2π(kT+fc)τp,tm≤t≤tm+2T(1+a)]]>将接收信号下变频到基带信号为：yj(t)=ejΩtΣi=1Nt(Σk=0K-1si[k]Hkej2πkT(1+a)t+w(t)),t∈[tm,tm+2T(1+a)]]]>w(t)为基带高斯白噪声。Ω＝2πafc是载波频偏(CFO,carrierfrequencyoffset)。可以看出，多普勒效应对传输信号造成了两点影响：1)引起了信号在时间上的扩展或压缩，即持续时间从2T变为2T/(1+a)。2)多普勒效应造成每一个子载波发生了的频率偏移，而由于水声信道的带宽与载波频率相差不大，不同的子载波的频率偏移不同，所以不能将其作为窄带信号处理。因此在进行OFDM解调前必须要对信号进行多普勒频偏补偿，以防止子载波间干扰(ICI,inter-carrierinterference)。在接收端进行采样，则可以得到基带数字信号为：yj(n)=yj(t)|t=tm+nTs]]>其中KF为发送训练序列的实际长度，是采样间隔，是过采样率。如图5所示，本发明采用在接收端使用多个并行相关器对信号进行处理的方法，实现多普勒扩展因子的估计。具体来说，每个相关器的相关窗口取不同的长度，对接收信号进行延时自相关运算，设相关器的窗口长度为Kl，利用接收信号中两段重复序列作自相关，得到：Cj(n)=Σk=0Kl-1yj(n+k)·(yj)*(n+k+Kl)]]>其中j代表第j根接收天线，Kl表示相关器的窗口长度，y表示接收信号,n表示第n个采样点。能量函数P为：Pj(n)=12Σk=0Kl-1[yj(n+k)·(yj)*(n+k)+yj(n+k+Kl)·(yj)*(n+k+Kl)]]]>判决变量为：Mn=Σj=1Nr|Cj(n)|2/Σj=1NrPj2(n)]]>取输出最大判决变量的窗口长度值则认为此值与受多普勒影响的信号长度是最接近的，所以，可以计算出多普勒扩展因子为：K^=KF(1+a^)⇒a^=KF-K^K^]]>其中，KF为发送训练序列的实际长度，为多普勒扩展因子的估计值。相应的，可得收发机的相对运动速度为：v^=ca^]]>其中，c为声波在海水中的传播速度。接收机利用估计出的对接收信号进行重采样，即可消除多普勒效应对传输信号的影响。图6—图7给出了不同信道多径数目的情况下，估计误差随信噪比的变化而变化的仿真曲线。可见，在所有的实施例中本发明的性能都明显优于传统的LFM方法和OFDM方法，随着多径数目的增加，LFM方法和OFDM方法的估计误差都略有增加，而本发明所提出的ZCZ序列具有不同序列之间的相关结果接近于0的性质，可有效抑制多径干扰，估计误差几乎不受多径数目增加的影响，稳定性较强，估计精度高。图8是在信噪比为0dB，多径数目为7的条件下，估计误差随收发端移动速度变化的仿真图，v取正值表示收发端相互靠近，取负值表示收发端相互远离。由图可知，LFM方法在移动速度增大时估计误差明显上升；OFDM方法的估计性能在移动速度较大时，也有一定的下降；而ZCZ方法的估计精度一直保持在10-4左右，且估计误差受移动速度增加的影响不大，具有很强的稳定性。当前第1页1 2 3