一种基于离散导频辅助的载波同步方法

摘要

本发明涉及一种基于离散导频辅助的载波同步方法。该方法要求输入信号具有非连续、等间隔分布的离散导频点的结构特点，方法实施步骤为：(1)从采集的数据符号中抽取导频符号；(2)根据抽取的导频符号获取载波频偏变化率，并根据所述载波频偏变化率对包含有所述导频符号的数据符号进行频偏变化率补偿；(3)根据频偏变化率补偿后的导频符号获取载波频偏，并根据所述载波频偏对包含有所述导频符号的数据符号进行频偏补偿；(4)将所述频偏补偿后的导频符号进行相干累加，得到载波初始相位；(5)通过载波初始相位对进行频偏补偿后的各数据符号进行相位校正，完成载波同步。

说明书

技术领域

本发明涉及猝发通信技术领域，特别涉及一种基于离散导频辅助的载波同步方法。

背景技术

在TDMA网络系统中，数据帧的长短决定了组网节点数目上限，进而决定着通信网络的容量。为保证网络具有较大容量并且兼顾数据帧的传输效率，帧长设计通常为10ms级别，并且在一些特殊应用领域如深远海通信，出于保密的需要，通信方式多选择猝发通信体制。由于数据帧持续时间较短，故在信噪比较低，多普勒动态较大场景下，传统的闭环跟踪方式受限于环路入锁时间，并不适用，因此基于导频辅助的开环跟踪方式是猝发通信的重要应用方向。

目前，在进行导频辅助的短时数据帧同步算法时，导频符号非连续、均匀分布在整个信号帧内，基于此通信格式，首先估计出导频符号的相位值，然后利用插值的方法，插出相邻导频符号之间的所有数据符号值，然而如果导频符号相位估计错误，则数据符号的相位也都校正错误，造成载波同步误差较大的问题。

发明内容

本发明的目的在于克服现有技术的不足，提供一种基于离散导频辅助的载波同步方法，利用导频符号进行频偏变化率估计，并对数据符号进行频偏变化率剥离，然后对频偏变化率剥离后的数据符号进行载波频偏估计，并对数据符号进行载波频偏剥离，最后进行初始相位校正，完成载波同步。可以解决现有载波同步误差较大，无法适用于大频率动态应用场景下的载波同步。

本发明的上述目的是通过如下技术方案予以实现：

一种基于离散导频辅助的载波同步方法，包括：

(1)从采集的数据符号中抽取导频符号；

(2)根据抽取的导频符号获取载波频偏变化率，并根据所述载波频偏变化率对包含有所述导频符号的数据符号进行频偏变化率补偿；

(3)根据频偏变化率补偿后的导频符号获取载波频偏，并根据所述载波频偏对包含有所述导频符号的数据符号进行频偏补偿；

(4)将所述频偏补偿后的导频符号进行相干累加，得到载波初始相位；

(5)通过载波初始相位对进行频偏补偿后的各数据符号进行相位校正，完成载波同步。

进一步地，所述数据符号的长度为N,其中，N＝K×(P+L)，K表示将长度为N的数据符号均匀化分为K段，L为每一段数据符号中包含的数据符号个数，P为每一段数据符号中导频符号的个数。

进一步地，所述步骤(2)中获取载波频偏变化率的方法为：

将每一段数据符号中的P个导频符号进行相干累加，得到K个累加结果p₁,p₂,…p_K；

将所述K个累加结果划分为两组，并将两组累加结果进行共轭相乘得到自相关函数序列其中，第一组累加结果为第二组累加结果为

对所述自相关函数序列进行FFTFFT运算，得到所述自相关函数序列对应的FFT谱线，其中，FFT运算的采样率其中R_s表示数据符号速率；

从所述FFT谱线中获取峰值谱线以及与所述峰值谱线前后相邻两根谱线分别对应的FFT变换值X_k，X_k-1，X_k+1，并且根据公式计算峰值位置的修正项δ，其中，Re[·]表示对[·]取实部；

根据公式计算修正后的载波频偏，其中，f_d为修正后的载波频偏，k_peak为修正后的峰值谱线频率，k_peak＝k+δ，k为修正前的峰值谱线频率；

根据公式计算载波频偏变化率，其中，f_d'为偏变化率。

进一步地，所述步骤(3)中获取载波频偏的方法为：

将进行频偏变化率补偿后的每一段数据符号中的P个导频符号进行相干累加，得到K个累加结果p₁,p₂,…p_K；

对累加结果p₁,p₂,…p_K进行FFT运算，得到所述自相关函数序列对应的FFT谱线；

从所述FFT谱线中获取峰值谱线以及与所述峰值谱线前后相邻两根谱线分别对应的FFT变换值X_k，X_k-1，X_k+1，并且根据公式计算峰值位置的修正项δ，其中，Re[·]表示对[·]取实部；

根据公式计算修正后的载波频偏，其中，f_d为修正后的载波频偏，k_peak为修正后的峰值谱线频率，k_peak＝k+δ，k为修正前的峰值谱线频率。

进一步地，所述步骤(4)中获取载波初始相位的方法为：

根据公式进行计算，其中，为初始相位，s_p,0为导频符号的初始值，S'_p为频偏补偿后的导频符号进行相干累加后的结果。

进一步地，所述步骤(5)中通过载波初始相位对进行频偏补偿后的各数据符号进行相位校正，完成载波同步的方法为：

将进行频偏补偿后的各数据符号与所述初始相位的倒数进行相乘运算。

本发明与现有技术相比具有如下有益效果：

(1)本发明提供的基于离散导频辅助的载波同步方法，在估计多普勒频率变化率和多普勒频偏时，均采用Jacobsen估计算法而非传统的FFT估计算法，即获取到FFT谱线后，对峰值谱线频率进行修正，并采用修正后的峰值谱线频率获取载波频偏变化率和载波频偏，从而可以保证该同步方法在低信噪比条件下，通过极小的FFT运算规模获得比传统FFT运算方案更高的估计精度，从而获得更理想的载波同步效果。

(2)本发明在频偏变化率估计、频偏估计以及最终的载波初始相位校正上，均采用相干处理方式，利用多个导频符号的相干累加提升判决信噪比，降低了错误判决概率，从而保证本发明提供的方法可以工作于低信噪比条件下。

(3)本发明提供的基于离散导频辅助的载波同步方法，首先利用导频符号进行多普勒频偏变化率估计，在多普勒频偏变化率补偿之后，再进行多普勒频偏估计。从而可以保证本发明提供的方法可以应用在更大的多普勒动态条件下，进而可以适用于大频率动态应用场景下的载波同步。

附图说明

图1为本发明的离散导频辅助的信号格式；

图2为本发明的信号处理流程图。

具体实施方式

下面结合附图2本发明的实现过程进行进一步详细阐述：

(1)从采集的数据符号中抽取导频符号。

对于本发明，要求输入的数据符号具有非连续、等间隔导频符号分布的特点，即每隔L个数据符号，有P个导频符号，信号结构可以如图1所示。则采集的数据符号符号长度N＝K×(P+L)，其中K表示将长度为N的数据符号均匀化分为K段，K的取值应满足累计导频符号信噪比的要求，通常需要满足其中，表示用dB表示的单个导频符号信噪比。并且预置输入信号已经完成粗同步，残余的多普勒频偏在导频图案频率估计的范围内。

(2)根据抽取的导频符号获取载波频偏变化率，并根据所述载波频偏变化率对包含有所述导频符号的数据符号进行频偏变化率补偿。

本发明中的频偏变化率估计方法采用延迟自相关FFT算法，算法原理描述如下：

设基带接收信号为：

其中，f_d为载波频偏，f′_d为载波频偏变化率，则信号自相关函数为：

其中，(·)^*表示(·)共轭函数。由此可见自相关函数实际上为频率为f′_dτ的单载波信号，因此频偏变化率估计的问题实际为频偏估计的问题。

本发明的频偏估计算法采用Jacobsen估计算法，具体为：

a、将每一段数据符号中的P个导频符号进行相干累加，得到K个累加结果p₁,p₂,…p_K。

b、将所述K个累加结果划分为两组，并将两组累加结果进行共轭相乘得到自相关函数序列其中，第一组累加结果为第二组累加结果为

c、对所述自相关函数序列进行FFT运算，得到所述自相关函数序列对应的FFT谱线。

其中，FFT在工程实现中可采用FPGA自带的IP核，进行FFT运算的点的个数为2的整数幂，如果不满足2的整数幂，可以通过补0的方式进行序列延长。此时，FFT运算的采样率其中R_s表示数据符号速率。

d、从所述FFT谱线中获取峰值谱线以及与所述峰值谱线前后相邻两根谱线分别对应的FFT变换值X_k，X_k-1，X_k+1，并且根据公式计算峰值位置的修正项δ，其中，Re[·]表示对[·]取实部。

e、根据公式计算修正后的载波频偏，其中，f_d为修正后的载波频偏，k_peak为修正后的峰值谱线频率，k_peak＝k+δ，k为修正前的峰值谱线频率。

f、根据公式计算载波频偏变化率，其中，f′_d为偏变化率。

对于本发明，此处进行频偏变化率估计，是为了解决大多普勒频偏变化率条件下，频偏无法直接估计的问题，同时，此处采用Jacobsen估计算法进行频率估计，进而计算出频偏变化率，是因为Jacobsen估计算法较传统的FFT算法相比，具有更小的FFT规模，并且在相同信噪比条件下，获得更高的估计精度。

(3)根据频偏变化率补偿后的导频符号获取载波频偏，并根据所述载波频偏对包含有所述导频符号的数据符号进行频偏补偿。

本发明中载波频偏算法采用Jacobsen估计算法，具体步骤为：

a、将进行频偏变化率补偿后的每一段数据符号中的P个导频符号进行相干累加，得到K个累加结果p₁,p₂,…p_K。

b、对累加结果p₁,p₂,…p_K进行FFT运算，得到所述自相关函数序列对应的FFT谱线。

c、从所述FFT谱线中获取峰值谱线以及与所述峰值谱线前后相邻两根谱线分别对应的FFT变换值X_k，X_k-1，X_k+1，并且根据公式计算峰值位置的修正项δ，其中，Re[·]表示对[·]取实部。

d、根据公式计算修正后的载波频偏，其中，f_d为修正后的载波频偏，k_peak为修正后的峰值谱线频率，k_peak＝k+δ，k为修正前的峰值谱线频率。

(4)将所述频偏补偿后的导频符号进行相干累加，得到载波初始相位。

具体地，根据公式进行计算，其中，为初始相位，s_p,0为导频符号的初始值，通常s_p,0＝1+j，S'_p为频偏补偿后的导频符号进行相干累加后的结果。其中，s′_p,j第j个导频符号，PK为存储的N个接收信号中，导频符号的总个数。

(5)通过载波初始相位对进行频偏补偿后的各数据符号进行相位校正，完成载波同步。

具体地，将进行频偏补偿后的各数据符号与所述初始相位的倒数进行相乘运算。然后进行后面的解调译码等操作。

实施例

假设接收信号具有图1所示结构，其中导频符号初始值为s₀＝1+j，业务数据符号采用QPSK调制方式，则接收信号可用下式描述：

其中，f_d为捕获之后残余的多普勒频偏，f′_d为多普勒频偏变化率，T_s为符号周期，s_n＝I_n+jQ_n为发送的第n个数据符号，n＝1,2,…,N。特别地，对于导频符号，s_n＝1+j。

按照步骤(1)提取数据符号中的所有导频符号，得到序列集合{s′_1,1,s'_2,1,…,s'_P,1,s₁'_,2,s'_2,2,…,s'_P,K}。

按照步骤(2)利用抽取的导频符号，通过Jacobsen估计算法计算出载波频偏变化率：

a、将抽取出来的导频符号按照P个为一组进行相干累加，得到K个累加结果p₁,p₂,…p_K，其中m＝1,2,…K。一般情况下，残余多普勒在连续P个导频符号内累计相位忽略不计，则其中，T′_s＝(P+L)T_s表示导频累加后样点序列p₁,p₂,…p_K的周期。为方便讨论，省略T′_s,于是m＝1,2,…K。

b、将K个样点划分为两组，第一组为第二组为：对应共轭相乘，得到自相关函数序列

c、对做N_FFT点的FFT运算。FFT采用FPGA自带的IP核，N_FFT为2的整数幂，如果不满足2的整数幂，可以通过补0的方式进行序列延长。

d、找出FFT的峰值谱线k其前后相邻两根谱线k-1,k+1所对应的FFT变换值X_k，X_k-1，X_k+1，然后按照下式得到峰值位置的修正项δ，δ可以为正数或者负数。

其中，Re[·]表示对[·]取实部。

e、修正后的峰值谱线k_peak及其对应的估计频偏分别为下式所示：

k_peak＝k+δ

f、根据f_d计算出频偏变化率的估计值为

此时，对式1的数据符号进行频偏变化率补偿，补偿后的接收符号序列为

其中，θ₀为补偿后的初始相位。

按照步骤(3)，利用经频偏变化率补偿后的导频符号，通过Jacobsen估计算法计算出载波频偏，具体地：

A、将经频偏变化率补偿后的导频符号按照P个为一组进行相干累加，得到K个累加结果p₁,p₂,…p_K，m＝1,2,…K。

B、对p₁,p₂,…p_K做N_FFT点的FFT运算。FFT在工程实现中可采用FPGA自带的IP核，N_FFT为2的整数幂，如果不满足2的整数幂，可以通过补0的方式进行序列延长。

C、获取FFT的峰值谱线k其前后相邻两根谱线k-1,k+1所对应的FFT变换值X_k，X_k-1，X_k+1，然后按照下式得到峰值位置的修正项δ，δ可以为正数或者负数。

其中，Re[·]表示对[·]取实部。

D、修正后的峰值谱线k_peak及其对应的估计频偏分别为下式所示。

k_peak＝k+δ

此时，对式2的数据符号进行频偏补偿，补偿后的接收符号序列为

按照步骤(4)，将频偏变化率和频偏都补偿过的导频符号抽取出来，进行相干累加得到然后跟导频符号的真实值进行比较，得到初始相位

按照步骤(5)，对式3数据符号进行初始相位校正，补偿后的接收符号序列为

s””(n)＝s_n>

s””(n)进行解调、译码等操作。

本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。