CN109391572A - 一种基于相位增量的载波频偏估计方法 - Google Patents

Abstract

本发明公开了一种基于相位增量的载波频偏估计方法，属于无线通信技术领域。本发明分为两步：第一步，通过设置有效延迟长度可得到能够估计较大载波频偏的自相关算法，再利用复信号指数化近似的原理可获得简化的L&R算法；第二步，借鉴自相关估计思想可设计出一种具有高精度的简化互相关算法。本发明可利用有限的导频开销估计出接近于50％符号速率的大频偏，且具有较低的复杂度和良好的解调性能，可用于包括航天测控、深空通信和卫星通信的导频受限短突发通信。

Description

一种基于相位增量的载波频偏估计方法

技术领域

本发明属于无线通信技术领域，特别是指一种基于相位增量的载波频偏估计方法，可用于包括航天测控、深空通信和卫星通信的导频受限短突发通信。

背景技术

近年来，短突发通信已经广泛应用于航天测控、卫星遥感以及无人机数据链等前沿领域，同时还将应用到第五代移动通信中。在这些通信领域中，有限的频谱资源会限制导频开销，同时通信双方的高速移动会产生多普勒效应，从而产生较大的多普勒扩展即载波频偏，进而造成接收端无法实现相干解调，以致通信质量严重下降。

为了矫正较大的载波频偏，传统的频偏估计算法可以分为数据辅助、非数据辅助和判决引导三大类。其中，非数据辅助和判决引导估计算法的信噪比门限和复杂度都高于数据辅助估计算法。因此，在短突发通信中，普遍采用数据辅助估计算法。这类估计算法一般利用已知的数据符号对载波参数进行估计。这些数据符号通常被称作导频符号。它们可以分插到数据帧的不同位置而形成不同的数据帧结构。

从采用的具体数据辅助估计算法而言，它可以分成频域数据辅助估计算法和时域数据辅助估计算法两小类。前者一般是基于搜索周期图的峰值来锁定载波频偏的，而后者一般是基于求解相关运算的相位增量来获得载波频偏的。在相同信噪比和导频开销下，前者的估计范围要比后者的大，而后者的估计精度要比前者的高。为了弥补这两种算法各自的缺陷，孙锦华等在“低信噪比下时频联合的载波估计算法”(西安交通大学学报,2015,49(2):62-68)一文中提出了一种时频联合的载波估计JTDFDCY算法。其基本原理是同时使用频域的旋转平均周期图算法和时域的互相关算法，获得了较大的估计范围和较高的估计精度，但同时也带来了较高的复杂度和时频域多参数优化配置的问题。另一方面，时域数据辅助估计算法可以进一步分成自相关算法和互相关算法。其中，自相关算法一般利用单个导频块对载波频偏进行估计。M.Luise et al.在“Carrier frequency recovery in all-digital modems for burst-mode transmissions”(IEEE Trans.Commun.,1995,43(234):1169-1178)一文中和U.Mengali et al.在“Data-aided frequency estimation forburst digital transmission”(IEEE Trans.Commun.,1997,45(1):23-25)一文中分别提出了利用单个导频块进行自相关运算来估计载波频偏的L&R算法和M&M算法。而互相关算法通常利用两个或多个不相交的导频块对载波频偏进行估计。在相同的信噪比和导频开销下，自相关算法和互相关算法具有以下的特点：前者的估计范围更大且信噪比门限及复杂度更低，后者的估计精度更高。因此，在导频受限的情况下，两者可以结合起来运用。

发明内容

有鉴于此，本发明提出一种基于相位增量的载波频偏估计方法，本方法能够在较少的导频开销下获得较大的估计范围、较低的复杂度和良好的解调性能，更加适应导频受限的短突发通信。

为了实现上述目的，本发明采用的技术方案是：

一种基于相位增量的载波频偏估计方法，其包括以下步骤：

(1)设定数据帧的格式，所述数据帧包括m-1个数据块和一个位于数据块之后的导频块，所述导频块的长度为L，每个数据块的长度为M_i，其中索引i从1到m-1，所述导频块用于存储导频信息；

(2)对数据帧进行加调得到调制信号s(k)；

(3)将调制信号s(k)通过空间信道传到接收端，在接收端得到导频信号r_p(k)；

(4)在采样时刻k对导频信号r_p(k)进行采样，并利用已知的导频信息对导频信号r_p(k)进行预处理，得到导频信号中第q个导频块的去调制信号z_q1(k)；

(5)利用去调制信号z_q1(k)及其延迟长度为l的延迟信号z_q1(k+1)，通过自相关运算，得到自相关值R(l)：

其中，z_q1(k)^*为去调制信号z_q1(k)的共轭；

(6)利用自相关值R(l)，通过自相关频偏估计，得到第一个频偏估计值

(7)利用第一个频偏估计值通过补偿去调制信号z_q1(k)，得到校正的去调制信号z_q2(k)：

其中，exp表示自然数e的指数，j为虚数单位，T_s为信道速率的符号周期，k为采样时刻；

(8)利用校正的去调制信号z_q2(k)，通过复信号指数化近似操作，得到简化的去调制信号z_q3(k)：

z_q3(k)＝|z_q2(k)|exp(-jarg{z_q2(k)})≈exp(-jarg{z_q2(k)})，

其中，arg{}为求辐角运算；

(9)根据去调制信号z_q3(k)和可变延迟长度为α的延迟信号z_q3(k+α)以及平滑噪声系数N，得到平滑自相关值R；

(10)利用平滑自相关值R，通过自相关频偏估计，得到第二个频偏估计值

(11)利用第一个频偏估计值和第二个频偏估计值通过补偿去调制信号z_q1(k)，得到校正的去调制信号z_q4(k)：

(12)利用校正的去调制信号z_q4(k)及其可变导频间隔为D_i的延迟信号z_q4(k+D_i)，通过互相关运算，得到互相关值R₁(D_i)：

其中，z_q4(k)^*为去调制信号z_q4(k)的共轭，可变导频间隔索引i＝1,2,...,m-1；

(13)利用互相关值R₁(D_i)，通过互相关频偏估计，得到第三个频偏估计值

(14)对得到的第一个频偏估计值第二个频偏估计值和第三个频偏估计值进行求和，得到最终的频偏估计值

具体的，所述步骤(4)的具体方式为：

(4a)在采样时刻k，按照数据帧的格式，对已接收到的导频信号r_p(k)进行遍历，得到对应于n个导频块的采样时刻集合k_p：

(4b)遍历采样时刻集合k_p，得到对应于第q个导频块的导频信号r_pq(k)，q＝1,2,...,n；

(4c)在接收端对已知的导频信息进行调制，得到导频信息的调制信号s_p(k)；

(4d)利用导频信号r_pq(k)以及导频信息的调制信号s_p(k)，通过共轭相乘运算，得到对应于第q个导频块的第一时刻去调制信号z_q1(k)：

z_q1(k)＝r_pq(k)*s_p(k)^*，

其中，s_p(k)^*为调制信号s_p(k)的共轭。

具体的，所述步骤(9)的具体方式为：

(9a)利用去调制信号z_q3(k)及其可变延迟长度为α的延迟信号z_q3(k+α)，通过自相关运算，得到自相关值R(α)：

其中，z_q3(k)^*为去调制信号z_q3(k)的共轭；

(9b)利用自相关值R(α)，通过对α求和，得到平滑自相关值R：

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

1.本发明由于采用了复信号指数化近似和自相关估计的思想，大大降低了对载波频偏估计的计算复杂度。

2.本发明由于先使用联合自相关估计再使用互相关估计，不仅在有限的导频开销下可以估计出较大的载波频偏，还可以获得很高的估计精度，提高了在短突发通信中的可行性。

总之，本发明可利用有限的导频开销估计出接近于50％符号速率的大频偏，且具有较低的复杂度和良好的解调性能，可用于包括航天测控、深空通信和卫星通信的导频受限短突发通信。

附图说明

图1是本发明的方法流程图；

图2是本发明中使用的突发帧格式；

图3是本发明与两种经典算法在不同信噪比下的频偏估计均方根误差性能仿真图；

图4是本发明在不同信噪比下的误比特性能仿真图。

具体实施方式

下面结合附图和具体实施方式对本发明进行详细描述。

一种基于相位增量的载波频偏估计方法，其包括以下步骤：

(1)设定数据帧的格式，所述数据帧包括m-1个数据块和一个位于数据块之后的导频块，所述导频块的长度为L，每个数据块的长度为M_i，其中索引i从1到m-1，所述导频块用于存储导频信息；

(2)对数据帧进行加调得到调制信号s(k)；

(3)将调制信号s(k)通过空间信道传到接收端，在接收端得到导频信号r_p(k)；

(4)在采样时刻k对导频信号r_p(k)进行采样，并利用已知的导频信息对导频信号r_p(k)进行预处理，得到导频信号中第q个导频块的去调制信号z_q1(k)；

(5)利用去调制信号z_q1(k)及其延迟长度为l的延迟信号z_q1(k+1)，通过自相关运算，得到自相关值R(l)：

其中，z_q1(k)^*为去调制信号z_q1(k)的共轭；

(6)利用自相关值R(l)，通过自相关频偏估计，得到第一个频偏估计值

(7)利用第一个频偏估计值通过补偿去调制信号z_q1(k)，得到校正的去调制信号z_q2(k)：

其中，exp表示自然数e的指数，j为虚数单位，T_s为信道速率的符号周期，k为采样时刻；

(8)利用校正的去调制信号z_q2(k)，通过复信号指数化近似操作，得到简化的去调制信号z_q3(k)：

z_q3(k)＝|z_q2(k)|exp(-jarg{z_q2(k)})≈exp(-jarg{z_q2(k)})，

其中，arg{}为求辐角运算；

(9)根据去调制信号z_q3(k)和可变延迟长度为α的延迟信号z_q3(k+α)以及平滑噪声系数N，得到平滑自相关值R；

(10)利用平滑自相关值R，通过自相关频偏估计，得到第二个频偏估计值

(11)利用第一个频偏估计值和第二个频偏估计值通过补偿去调制信号z_q1(k)，得到校正的去调制信号z_q4(k)：

(12)利用校正的去调制信号z_q4(k)及其可变导频间隔为D_i的延迟信号z_q4(k+D_i)，通过互相关运算，得到互相关值R₁(D_i)：

其中，z_q4(k)^*为去调制信号z_q4(k)的共轭，可变导频间隔索引i＝1,2,...,m-1；

(13)利用互相关值R₁(D_i)，通过互相关频偏估计，得到第三个频偏估计值

(14)对得到的第一个频偏估计值第二个频偏估计值和第三个频偏估计值进行求和，得到最终的频偏估计值

具体的，所述步骤(4)的具体方式为：

(4a)在采样时刻k，按照数据帧的格式，对已接收到的导频信号r_p(k)进行遍历，得到对应于n个导频块的采样时刻集合k_p：

(4b)遍历采样时刻集合k_p，得到对应于第q个导频块的导频信号r_pq(k)，q＝1,2,...,n；

(4c)在接收端对已知的导频信息进行调制，得到导频信息的调制信号s_p(k)；

(4d)利用导频信号r_pq(k)以及导频信息的调制信号s_p(k)，通过共轭相乘运算，得到对应于第q个导频块的第一时刻去调制信号z_q1(k)：

z_q1(k)＝r_pq(k)*s_p(k)^*，

其中，s_p(k)^*为调制信号s_p(k)的共轭。

具体的，所述步骤(9)的具体方式为：

(9a)利用去调制信号z_q3(k)及其可变延迟长度为α的延迟信号z_q3(k+α)，通过自相关运算，得到自相关值R(α)：

其中，z_q3(k)^*为去调制信号z_q3(k)的共轭；

(9b)利用自相关值R(α)，通过对α求和，得到平滑自相关值R：

本发明的技术思路是：在接收端，任取一段导频信号进行去调制操作，得到去调制信号，通过对这些去调制信号进行延迟长度为1的自相关运算，再利用复信号指数化近似操作得到简化的L&R算法，最后借鉴自相关估计思想获得简化的互相关算法，实现以较低的复杂度对较大的频偏进行高精度的估计。

本发明方法解决了导频受限短突发通信中存在经典载波频偏估计方法无法兼顾低复杂度和高精度的问题，其技术方案概括来说包括如下步骤：

1.在接收端，提取任意相邻的两段导频信号进行去调制操作，得到两段去调制信号；

2.对前一段去调制信号及其延迟长度为1的延迟信号进行自相关运算，利用该自相关运算的结果得到第一个频偏估计值；

3.利用第一个频偏估计值对前一段去调制信号进行补偿和复信号指数化近似操作，得到简化的去调制信号；

4.对该简化的去调制信号及其可变延迟长度为的延迟信号进行自相关运算并求和，利用该求和的结果得到第二个频偏估计值；

5.利用第一个频偏估计值和第二个频偏估计值对两段去调制信号进行补偿；

6.借鉴自相关估计的思想，对补偿后的两段去调制信号进行互相关运算，利用该互相关运算的结果得到第三个频偏估计值；

7.对这三个频偏估计值进行求和，得到最终的载波频偏估计值。

具体来说，参照图1，一种基于相位增量的载波频偏估计方法，其具体步骤如下：

步骤1，设置突发帧格式。

设定数据帧的格式，所述数据帧包括m-1个数据块和一个位于数据块之后的导频块，所述导频块的长度为L，每个数据块的长度为M_i，其中索引i从1到m-1，所述导频块用于存储导频信息；

步骤2，对数据帧进行加调得到信号调制信号s(k)；

步骤3，将调制信号s(k)通过空间传到接收端，在接收端得到导频信号r_p(k)；

步骤4，进行去调制操作。

利用接收的导频信号r_p(k)与其对应的调制信号s_p(k)，通过共轭相乘运算，得到去调制信号z(k)：

其中，为按照突发帧格式B的复用结构设置采样时刻k得到的对应于m个导频块的采样时刻集合，s_p(k)^*为能量归一化的调制信号s_p(k)的共轭，f_d为载波频偏，θ为相偏，T_s表示符号周期，为噪声项，n(k)～CN(0,N₀)表示均值为0、实部和虚部方差均为的圆对称复高斯随机变量，j是虚数单位，r_p(k)为接收信号，其形式如下：r_p(k)＝s_p(k)exp[j(2πf_dT_sk+θ)]+n(k)。

步骤5，进行自相关运算。

5a)将去调制信号z₁(k)延迟1个长度，得到延迟信号z₁(k+1)；

5b)对去调制信号z₁(k)及其延迟信号z₁(k+1)进行共轭相乘并求和，得到自相关值R(l)：

其中，为归一化因子，z₁(k)^*为去调制信号z₁(k)的共轭，Ψ(1)为噪声累加项，其表达形式如下：

5c)对自相关值R(l)进行取幅角运算，得到第一个频偏估计值

步骤6，进行频偏补偿和复信号指数化近似。

6a)利用得到的第一个频偏估计值通过补偿去调制信号z₁(k)，得到校正的去调制信号z₂(k)：

其中，

6b)利用校正的去调制信号z₂(k)，通过复信号指数化近似操作，得到简化去调制信号z₃(k)：

步骤7，进行自相关运算并求和。

7a)将去调制信号z₃(k)延迟可变的α个长度，得到延迟信号z₃(k+α)；

7b)对去调制信号z₃(k)及其延迟信号z₃(k+α)进行共轭相乘并求和，得到自相关值R(α)：

其中，为归一化因子，z₃(k)^*为去调制信号z₃(k)的共轭；

7c)对自相关值R(α)进行求和，得到平滑自相关值R：

其中，N为平滑噪声系数；

7d)对平滑自相关值R进行取幅角运算，得到第二个频偏估计值

步骤8，进行频偏补偿。

8a)利用得到的第一个频偏估计值和第二个频偏估计值通过补偿去调制信号z(k)，得到校正的去调制信号z₄(k)：

其中，采样时刻k∈k_p；

步骤9，进行基于自相关估计思想的互相关运算。

9a)将校正的去调制信号z₄(k)延迟可变的D_i个长度，得到延迟信号z₄(k+D_i)；

9b)对去调制信号z₄(k)及其延迟信号z₄(k+D_i)进行共轭相乘并求和，得到互相关值R₁(D_i)：

其中，Im为取虚部操作，z₄(k)^*为去调制信号z₄(k)的共轭；

9c)对互相关值R₁(D_i)进行取幅角运算，得到第三个频偏估计值

步骤10，进行频偏求和。

10a)对得到的第一个频偏估计值第二个频偏估计值和第三个频偏估计值进行求和，得到最终的频偏估计值

其中，为简化互相关算法利用第i个导频间隔D_i而获得的频偏估计值。

本发明的效果可通过如下仿真进一步说明：

1.仿真条件

调制方式为二进制相移键控BPSK，归一化频偏f_dT_s＝0.45，导频块数m＝2，数据块数m-1＝1，导频长度L＝20，数据长度M＝90。

2.仿真内容

仿真1：对2个导频块和1个数据块进行二进制相移键控BPSK调制，再经过加性高斯白噪声AWGN信道加噪处理，当归一化频偏f_dT_s＝0.45时，在不同信噪比E_s/N₀下，进行M&M算法估计、JTDFDCY算法估计和本文算法估计并通过蒙特卡罗仿真统计频偏估计均方根误差，仿真结果如图3所示。

图3中以星形标记的曲线表示在不同信噪比E_s/N₀下，进行M&M算法估计后的频偏估计均方根误差。

图3中以三角形标记的曲线表示在不同信噪比E_s/N₀下，进行JTDFDCY算法估计后的频偏估计均方根误差。

图3中以方形标记的曲线表示在不同信噪比E_s/N₀下，进行本文算法估计后的频偏估计均方根误差。

由图3可以看出，本文算法的估计性能与JTDFDCY算法的估计性能相当，且在较高的信噪比下都比M&M算法的估计性能高很多，即超过了1个数量级。

仿真2：对2个导频块和2个数据块进行二进制相移键控BPSK调制，再经过加性高斯白噪声AWGN信道加噪处理，当归一化频偏f_dT_s＝0.45时，在不同信噪比E_s/N₀下，进行本文算法估计并通过蒙特卡罗仿真统计错误比特数，仿真结果如图4所示。

图4中以方形标记的曲线表示在不同信噪比E_s/N₀下，进行本文算法估计后的误比特性能。

图4中虚线表示在不同信噪比E_s/N₀下，没有频偏的理想情况下的误比特性能。

由图4可以看出，在仅用40个导频符号同时携带590个数据符号的情况下，此时导频开销本文算法仍取得了良好的解调性能。

总之，本发明提出了一种适用于导频受限短突发通信的两步频偏估计算法。其实现分为两步：第一步，通过设置有效延迟长度可得到能够估计较大载波频偏的自相关算法，再利用复信号指数化近似的原理可获得简化的L&R算法；第二步，借鉴自相关估计思想可设计出一种具有高精度的简化互相关算法。最后，仿真结果表明，在仅用40个导频符号(6.4％的导频开销)的条件下，本文算法仍获得了较大的估计范围、较低的复杂度和良好的解调性能。

Claims (3)

1.一种基于相位增量的载波频偏估计方法，其特征在于，包括以下步骤：

(1)设定数据帧的格式，所述数据帧包括m-1个数据块和一个位于数据块之后的导频块，所述导频块的长度为L，每个数据块的长度为M_i，其中索引i从1到m-1，所述导频块用于存储导频信息；

(2)对数据帧进行加调得到调制信号s(k)；

(3)将调制信号s(k)通过空间信道传到接收端，在接收端得到导频信号r_p(k)；

(4)在采样时刻k对导频信号r_p(k)进行采样，并利用已知的导频信息对导频信号r_p(k)进行预处理，得到导频信号中第q个导频块的去调制信号z_q1(k)；

(5)利用去调制信号z_q1(k)及其延迟长度为l的延迟信号z_q1(k+1)，通过自相关运算，得到自相关值R(l)：

其中，z_q1(k)^*为去调制信号z_q1(k)的共轭；

(6)利用自相关值R(l），通过自相关频偏估计，得到第一个频偏估计值

(7)利用第一个频偏估计值通过补偿去调制信号z_q1(k)，得到校正的去调制信号z_q2(k)：

其中，exp表示自然数e的指数，j为虚数单位，T_s为信道速率的符号周期，k为采样时刻；

(8)利用校正的去调制信号z_q2(k)，通过复信号指数化近似操作，得到简化的去调制信号z_q3(k)：

z_q3(k)＝|z_q2(k)|exp(-jarg{z_q2(k)})≈exp(-jarg{z_q2(k)})，

其中，arg{}为求辐角运算；

(9)根据去调制信号z_q3(k)和可变延迟长度为α的延迟信号z_q3(k+α)以及平滑噪声系数N，得到平滑自相关值R；

(10)利用平滑自相关值R，通过自相关频偏估计，得到第二个频偏估计值

(11)利用第一个频偏估计值和第二个频偏估计值通过补偿去调制信号z_q1(k)，得到校正的去调制信号z_q4(k)：

(12)利用校正的去调制信号z_q4(k)及其可变导频间隔为D_i的延迟信号z_q4(k+D_i)，通过互相关运算，得到互相关值R₁(D_i)：

其中，z_q4(k)^*为去调制信号z_q4(k)的共轭，可变导频间隔索引i＝1,2,...,m-1；

(13)利用互相关值R₁(D_i)，通过互相关频偏估计，得到第三个频偏估计值

(14)对得到的第一个频偏估计值第二个频偏估计值和第三个频偏估计值进行求和，得到最终的频偏估计值

2.根据权利要求1所述的基于相位增量的载波频偏估计方法，其特征在于，所述步骤(4)的具体方式为：

(4a)在采样时刻k，按照数据帧的格式，对已接收到的导频信号r_p(k)进行遍历，得到对应于n个导频块的采样时刻集合k_p：

(4b)遍历采样时刻集合k_p，得到对应于第q个导频块的导频信号r_pq(k)，q＝1,2,...,n；

(4c)在接收端对已知的导频信息进行调制，得到导频信息的调制信号s_p(k)；

(4d)利用导频信号r_pq(k)以及导频信息的调制信号s_p(k)，通过共轭相乘运算，得到对应于第q个导频块的第一时刻去调制信号z_q1(k)：

z_q1(k)＝r_pq(k)*s_p(k)^*，

其中，s_p(k)^*为调制信号s_p(k)的共轭。

3.根据权利要求1所述的基于相位增量的载波频偏估计方法，其特征在于，所述步骤(9)的具体方式为：

(9a)利用去调制信号z_q3(k)及其可变延迟长度为α的延迟信号z_q3(k+α)，通过自相关运算，得到自相关值R(α)：

其中，z_q3(k)^*为去调制信号z_q3(k)的共轭；

(9b)利用自相关值R(α)，通过对α求和，得到平滑自相关值R：