王 中，陈永兴，樊书宏，岳 玲

(1.中国船舶集团有限公司 第七〇五研究所，陕西 西安 710077；2.水下信息与控制国防重点实验室，陕西 西安 710077)

0 引 言

水声通信是实现水下各类装备信息化发展的重要基础。近年来，各种水下无人航行器快速发展，其智能化和信息化对移动水声通信技术提出了日益迫切的需求。移动水声通信是水下移动节点执行任务、传输信息的重要途径，也是构建水下网络、提升水下信息化整体发展能力的必备技术[1]。但目前定点水声通信研究较多，移动水声通信技术研究相对而言尚存不足，特别是收发两端相对运动速度较快的高速移动水声通信技术研究甚少。

移动水声通信中收发两端的相对运动会导致通信信号出现时间压缩或扩展而产生多普勒效应，由于声波的传播速度低，因此移动水声通信中的多普勒效应影响相对严重。较大的多普勒频偏将引起严重的载波频率偏移，直接影响通信系统性能。通过采用帧同步信号作为同步码，为接收信息码信号准确开窗，并进行多普勒因子准确估计与补偿，是实现移动水声通信和保证高质量通信的前提。

针对移动水声通信中的帧同步和多普勒估计，国内外学者已经开展了一系列的研究，并提出不同的处理方法。Sharif 等[2]提出了在数据块前后添加LFM 信号，接收端根据相关峰值估计多普勒因子。李云飞[3]提出了利用Zadoff-Chu 序列，通过给训练序列分别添加循环前缀和后缀，进行多普勒因子估计。岳玲等[4]提出了利用HFM 信号和单频脉冲信号等估计多普勒频移。同时，还将宽带模糊度函数方法引入水声通信同步技术，以实现高精度多普勒因子估计。

总体来看，可以将当前这些多普勒估计方法分别归结为利用LFM 信号、m 序列、HFM 信号的3 类移动水声通信帧同步方法，这些方法均从不同角度给出了移动水声通信的帧同步解决途径。但对于特定的移动水声通信应用需求，还需要充分了解各种方法的具体适用环境和技术优势与特点等，这就需要对不同方法进行全面深入的性能分析与对比，为各方法优化选取和实际应用提供有效指导。

本文在分别研究LFM 信号、m 序列信号、HFM 信号3 种帧同步方法实现原理的基础上，针对AWGN 信道和实测多途衰落信道，通过计算机仿真分析，分别从信噪比、信道特性、收发相对多普勒速度、计算复杂度等方面，对这3 种帧同步方法在不同条件和应用环境的多普勒因子估计性能进行全面对比分析，为水下各类无人航行器移动水声通信帧同步方法的具体选择和综合应用提供参考。

1 移动水声通信的帧同步

对于发射和接收端存在相对运动速度的移动水声通信系统，在收发端多普勒效应的影响下，水声信号不仅在频域发生频率偏移，同时也表现为时域信号波形的压缩或扩展效果[5]。因此，在水声通信中，接收信号通常建模为时域波形上的伸缩效应，即

式中：s(t)为发射信号；r(t)为接收信号；α为多普勒因子，其定义为收发两端相对径向运动速度v与声速c之比，即

一般情况下，当通信双方相向移动时，v取正，反之取负。

水声信号经过不同的路径到达接收端，接收信号是所有路径传送信号的叠加。一般情况下，每条传播路径上的多普勒因子不同，但其主要部分都是由收发两端相对运动引起的，故假设所有传播路径上的多普勒因子均相等，且多普勒因子在一帧信号持续时间内保持不变。

因此，对于频率为f的发射信号，接收信号的频率改变为：

水声通信信号的时间展宽或压缩，将导致码元信号同步误差的累积。例如，对于水下相对运动速度30 kn 的收发系统，则会产生1% 的时间压缩或扩展量，其影响表现为在传输100 个码元信息后，同步误差累计为1 个码元宽度[5]，严重影响后续码元信息的解调和解码。因此，高速移动水声通信必须减缓严重多普勒效应的影响。

帧同步是移动水声通信中消除严重多普勒效应影响的重要方法。发射端在码元信号前后附加一定的帧同步信号作为同步码，接收端通过解析同步码，为接收信息码信号准确开窗，并对多普勒因子进行准确估计。因此，多普勒因子估计性能是帧同步方法的主要性能指标。在此基础上，利用多普勒补偿技术，可以消除多普勒效应对码元信号的影响，从而保证高质量的信息通信。

接收端根据帧同步信号的多普勒因子估计结果，可以通过重采样的方式进行多普勒补偿，也就是通过多普勒效应的逆过程实现对多普勒因子的补偿，即

2 帧同步的多普勒因子估计方法

2.1 基于LFM 信号的块估计方法

LFM 信号是指瞬时频率随着时间线性变化的一种信号，其时间函数可以表示为：

式中：A为信号幅度；f0为信号的中心频率；k=B/T为调频斜率；B为调频宽度；T为信号时间。

块估计方法是在传输数据前后各加上一段相同的多普勒不敏感信号，常用LFM 信号，通过比较接收信号与发射信号的帧长度变化估计多普勒因子[2,6]。

假设发射信号的帧长度为Ttp，接收信号的帧长度为Trp，则多普勒因子可以通过下式估计：

式中，Trp为通过接收端2 个相关器检测接收信号前后2 个相关峰值的间隔得到。图1 给出了块估计方法示意图。

图1 块估计方法Fig.1 Block estimation method

图中，多普勒估计的分辨率取决于LFM 信号的带宽B和发射数据帧长度Ttp，分辨率[7]近似为：

同时，宽带LFM 信号的多普勒容限[7]为：

根据式(7) 和式(8) 可知，发射信号带宽B越宽，数据帧Ttp越长，多普勒估计精度越高，多普勒容限越低。因此，对于收发端存在相对高速运动情况，接收信号相关峰将会出现分裂，导致多普勒估计性能下降。

2.2 基于m 序列的模糊度函数方法

m序列是一种具有良好自相关特性的伪随机序列，其循环移位自相关[8]定义为：

式中：s[k]为序列中的元素；*为复共轭；[k+n]N为循环移位。

宽带模糊度函数定义为：

式中，τ为时间延迟。模糊度函数描述了由多普勒因子和时间延迟共同导致的匹配滤波器失配情况。

基于模糊度函数估计多普勒因子的帧同步方法，其数据帧格式如图2 所示，前导信号使用对多普勒敏感的m 序列。

图2 模糊度函数方法的帧格式Fig.2 Frame format for the ambiguity function method

在实际应用中，由于接收机无法知道r(t)的到达时间，故需要进行时延 τ和多普勒因子 α的二维搜索，该过程如图3 所示。

图3 多普勒因子估计结构图Fig.3 Doppler factor estimation structure diagram

由图3 可知，在接收端，将本地保存的受不同多普勒效应影响的副本信号与接收信号中的前导信号作拷贝相关，相关器组输出的最大峰值便是多普勒估计结果，即

该方法可以实时处理数据，其多普勒估计精度主要依赖于前导信号对多普勒的敏感程度以及接收端相关器的个数。为了获得最佳的多普勒分辨率，前导信号的自模糊度函数应具有非常尖锐的主峰，因此常用伪随机信号。接收端相关器的个数主要由多普勒估计范围与精度决定，如果多普勒因子的搜索范围为±1%，精度要求为0.02%，则需要100 个相关器。实际应用中，该方法采用两步搜索方式，即先进行粗搜索，再进行细估计，虽然可以减少相关器的个数，但是计算开销仍然很大。

2.3 基于双HFM 信号的时延差估计方法

HFM 信号的瞬时频率随着时间呈对数变化，其上扫频和下扫频HFM 信号的时间函数[8]分别为：

以上扫频HFM 信号为例，它的瞬时频率为：

当存在多普勒效应时，该信号的瞬时频率改变为：

比较式(14)和式(15)，可以发现存在一个时延差τ，使得fdop(t-τ)=fup(t)，其中：

由此可见，通过解算时延τ，即可估计多普勒因子α，而时延τ可以利用副本相关处理方法进行精确解算。

通过时序连续发射或者同时叠加发射上扫频与下扫频HFM 信号作为前导信号，同步帧结构分别如图4所示。其中，图4(b)的帧结构中前导信号持续时间更短，可以降低发送端能耗，提高帧的数据利用率。

图4 发射帧结构Fig.4 The structure of transmitted frame

接收端设置2 个并行相关器对接收信号进行拷贝相关处理[9]，分别对接收信号中的HFM+和HFM-信号进行匹配滤波，实现过程如图5 所示。

图5 多普勒因子估计结构图Fig.5 The structure of Doppler factor estimation

图中，通过拷贝相关处理，2 个并行相关器分别输出上扫频和下扫频HFM 信号的相关峰偏移量 τ1和τ2，根据式(16)可以得到：

由式(17)和式(18)可得，2 个相关峰之间的差值Δτ可以表示为：

因此，通过式(19)可以解算出多普勒因子 α的估计值为：

由式(20)可以看出，多普勒因子 α的估计值与匹配滤波峰值时延差 Δτ，以及同步导引信号的持续时间T和带宽B有关。

3 性能仿真及对比分析

3.1 仿真条件

针对AWGN 信道和实测水声多径信道条件下高速移动节点水声通信的2 种典型应用场景，利用3 种多普勒因子估计方法进行计算机仿真。仿真时，发射信号的帧结构分别如图1 和图2 及图4(b)所示。LFM 信号和HFM 信号的参数设置均为起始频率fmin=18 kHz，终止频率fmax=22 kHz，信号持续时间T=50 ms，接收端采样频率为100 kHz。基于m 序列的模糊度函数方法采用DPSK 调制的511 点m 序列作为前导信号，载波频率为fc=20 k H z，信号带宽和持续时间与LFM，HFM 信号保持一致，采用两步搜索。

分别针对不同信噪比、信道特性、收发相对多普勒速度、计算复杂度等条件，进行仿真分析，每种情况下仿真次数均为1 000 次，统计分析各种方法对多普勒因子估计的均方误差性能。

3.2 不同信道特性条件下对比分析

实测水声多径信道参数如表1 所示，其中信道的多普勒扩展为1 Hz。

表1 水声信道仿真参数Tab.1 Underwater acoustic channel simulation parameters

在-10 dB 信噪比情况下，分别针对AWGN 信道和时变多途衰落信道，各种方法多普勒估计性能仿真结果如图6 所示。

图6 多普勒因子估计性能Fig.6 Doppler factor estimation performance

对比分析图6 可以发现，多途衰落信道条件下的多普勒估计性能明显低于AWGN 信道，这是由于多径效应的存在，增加了各路径信号之间的干扰，导致多普勒因子估计误差增大。从图6 均可以看出，m 序列模糊度函数算法的多普勒因子估计性能最好，双HFM 信号估计方法次之，LFM 信号块估计方法最差。

3.3 不同信噪比条件下对比分析

在实测水声信道中，为分析3 种方法在不同信噪比条件下的多普勒因子估计性能，将各方法分别在-10 dB 和0 dB 下进行仿真对比，结果如图7 所示。

图7 多普勒因子估计性能Fig.7 Doppler factor estimation performance

可以看出，当信噪比从-10 dB 增大到0 dB，m 序列模糊度函数方法的估计性能提升到10-5以下，双HFM 信号多普勒估计方法的估计性能提升到10-4以下。而LFM 信号块多普勒估计方法的估计性能变化不明显，这是由于在此仿真条件下，LFM 信号的多普勒容限约为vd=13.05 kn，即6.525 m/s，当相对运动速度超过容限值时，LFM 信号的相关峰将出现分裂，导致多普勒估计精度降低。

3.4 不同多普勒速度条件下对比分析

在实测水声信道中，为分析3 种方法在不同多普勒速度条件下的多普勒估计性能，将3 种多普勒估计方法分别在相对运动速度15 m/s 和30 m /s 条件下作对比，仿真结果如图8 所示。

图8 多普勒因子估计性能Fig.8 Doppler factor estimation performance

可知，给定多普勒速度，3 种多普勒估计方法的估计误差随信噪比增大而减小。由于仿真的多普勒速度远远超过LFM 信号的多普勒容限，即使增大信噪比，块估计方法也不能得到理想的估计性能。模糊度函数方法在多普勒速度15 m/s 和30 m/s 时均表现出良好的估计性能，但是付出了计算复杂度大的代价，工程实现起来较为困难。当多普勒速度分别为15 m/s和30 m/s 时，基于双HFM 信号的多普勒估计方法分别在-14 dB 和-12 dB 信噪比实现10-4量级的估计误差。

3.5 计算复杂度对比分析

假定各方法使用拷贝相关器的计算复杂度均相等。LFM 信号块多普勒估计方法和双HFM 信号多普勒估计方法均仅需要2 次拷贝相关，而对于m 序列模糊度函数方法，采用粗细两步搜索，即先以Δv1=2 m/s在速度区间[0,30]进行粗搜索得到v1，需要16 个拷贝相关器，再以Δv2=0.2 m/s 在区间[v1-Δv,v1+Δv]进行细搜索得到v2，需要21 个拷贝相关器，总共需要37 个拷贝相关器。表2 列出了3 种多普勒估计方法的计算复杂度比较。

表2 三种多普勒估计方法的计算复杂度比较Tab.2 Comparison of computational complexity of three Doppler estimation methods

3.6 仿真结果综合分析

根据以上仿真分析可以得出下面结论：基于LFM信号的块多普勒估计方法适用于静止节点或者移动速度较低，且对多普勒估计精度要求不高的水声通信设备；基于m 序列的模糊度函数方法计算复杂度高、运算量大，对移动水声通信接收机设备提出了十分严苛的要求；基于双HFM 信号多普勒估计方法的计算复杂度与块估计方法相同，且该方法在低信噪比和高多普勒速度条件下，仍能实现10-4量级的误差性能。

综上，双HFM 信号多普勒估计方法估计精度高、计算复杂度低、帧数据利用率高，可以在-10 dB 信噪比下实现对多普勒速度30 m/s 的多普勒因子精确估计，同时也可以实现复杂信道移动水声通信的帧同步。

4 结 语

对于移动水声通信系统，利用帧同步信号进行精确多普勒估计是不可或缺的。本文在分析多普勒效应对水声通信信号影响的基础上，分别研究LFM 信号、m 序列信号、HFM 信号3 种帧同步方法，给出了各种方法的具体实现原理，针对AWGN 信道和实测多途衰落信道，通过计算机仿真分析，分别从信噪比、信道特性、收发相对多普勒速度、计算复杂度等方面，全面对比分析了这3 种帧同步方法在不同条件和应用环境的多普勒因子估计性能，明确了各方法的优缺点和可适用的应用场景。仿真结果表明，3 种方法各有优势，均可适用不同的移动水声通信应用环境。但相对而言，HFM 信号同步方法多普勒估计精度高、数据传输速率快、计算复杂度低、稳健性强，可以实现复杂信道移动水声通信的帧同步。研究结果可为水下各类无人航行器移动水声通信帧同步方法的具体选择和综合应用提供参考。