一种TD‑AltBOC信号的时延标定方法与流程

本发明涉及卫星导航、时间频率、数字信号处理等领域，特别涉及卫星导航td-altboc信号的发射设备时延标定问题。

背景技术：

卫星导航系统基于到达时间(timeofarrival,toa)原理进行测距定位，即采用伪码测距的技术测量信号从已知位置的发射源至用户接收机所经历的时间，使用伪距测量的方法可以获得从发射源至接收机的组合时延。由于系统处理数据时，需要准确处理导航系统中的收发设备时延值，因此收发设备时延是卫星定位导航最主要的误差项之一。收发设备时延不仅会影响到卫星导航系统的定位、授时精度以及导航性能，同时也会影响卫星导航系统正常运行工作。

既有的卫星导航信号采用bpsk或qpsk调制方式，对于这两类信号已有成熟的时延标定方法，如魏海涛等人提出的“卫星导航设备时延精密标定方法与测试技术研究”以及钟兴旺等人提出的“bpsk调制器发射时延的一种高精度标定方法”等。该标定方法通过对时延位置处的相位反转点进行多次采样，进一步采用最小包络检测法计算其统计特性，从而获得设备时延的均值及稳定性，该方法已广泛应用在工程实践当中。

随着全球卫星导航系统的建设，一系列新体制导航信号投入应用，其中包括td-altboc等信号。由于信号体制发生变化，既有的bpsk信号时延标定方法无法直接应用于新一批卫星导航设备的时延标定中，这对卫星导航td-altboc等信号时延标定方法研究提出了急迫的需求，围绕此需求，本发明方法对相关内容展开了一系列探索。

在实际的卫星导航时延标定工作中，一般分为发射设备时延标定和接收设备时延标定。在目前的标定技术中，接收设备时延采用时延传递法获得，即组合时延扣除发射设备时延。因而，在时延标定过程中，发射设备时延的标定是最关键的。本发明方法主要研究td-altboc信号的发射设备时延标定方法。

技术实现要素：

本发明的技术解决问题是：针对td-altboc信号，提出一种有效的时延标定方法，该方法能够解决td-altboc信号的上、下边带各支路发射时延标定问题，同时该方法能够兼容bpsk信号以及qpsk信号的时延标定。

本发明的技术解决方案是：一种td-altboc信号的时延标定方法，步骤如下：

(1)搭建td-altboc信号时延标定系统；所述的td-altboc信号时延标定系统包括td-altboc信号发射设备、时频信号产生设备、高速信号采集装置、时延处理装置、触发信号线缆以及射频信号线缆；

(2)高速信号采集装置采集并存储预设长度的td-altboc信号；所述的td-altboc信号包含上边带信号和下边带信号；

(3)根据td-altboc信号4个支路的扩频码和时隙分配生成上、下边带信号各支路的等效扩频码；所述的4个支路的扩频码由上边带数据支路扩频码、上边带导频支路扩频码和下边带数据支路扩频码、下边带导频支路扩频码组成；

(4)将信号中的上边带信号和下边带信号分别进行下变频一一对应得到上边带基带信号和下边带基带信号；将上边带基带信号和下边带基带信号分别与上边带各支路等效扩频码和下边带各支路等效扩频码一一对应作滑动相关运算得到滑动相关运算结果，并搜索得到各支路最大相关峰值；

(5)在所有滑动相关运算结果中分别搜索各支路最大相关峰值所在位置，根据各支路最大相关峰值所在位置计算得到td-altboc信号的链路时延值；

(6)用链路时延值扣除触发信号线缆的时延值和射频信号线缆的时延值后，得到td-altboc信号的时延值。

其中，步骤(2)所述的高速信号采集装置以10mhz信号作为参考信号，以1pps信号作为触发信号，对足够长度的td-altboc信号以满足时延计算精度的采样率进行采集存储。

其中，步骤(3)所述的根据td-altboc信号4个支路的扩频码和时隙分配生成上、下边带信号各支路的等效扩频码的方法具体为，在原数据支路扩频码的每个扩频码片之后插入一个零值码片生成数据支路等效扩频码；在原导频支路扩频码的每个扩频码片之前插入一个零值码片生成导频支路等效扩频码。

其中，：步骤(4)具体包括以下步骤：

(401)分别递增上边带信号和下边带信号的下变频相位，并分别根据递增后的下变频相位一一对应对上边带信号和下边带信号进行下变频得到上边带基带信号和下边带基带信号；

(402)对上边带基带信号和下边带基带信号分别与各支路等效扩频码一一对应进行滑动相关运算，记录相关峰值大小，返回步骤(401)，直到下变频相位递增满2π；

(403)搜索得到记录的最大相关峰值。

其中，步骤(4)所述的滑动相关运算具体为：

其中，r(τ)为滑动相关运算结果，n为基带信号的点数，fs为高速信号采集装置的采样率，fc为码速率，t为扩频码周期，prn为扩频码，为向上取整，为向下取整，x(i)为基带信号。

其中，步骤(5)所述的链路时延值为：

其中，td_calc为链路时延值，pos_max为相关峰值所在位置序号，n为基带信号的点数，fs为高速信号采集装置的采样率；

其中，步骤(6)所述的用链路时延值扣除触发信号线缆的时延值和射频信号线缆的时延值后，得到td-altboc信号的时延值；具体为：td＝td_calc+tline1pps-tlinerf；

其中，td_calc为链路时延值，tline1pps为触发信号线缆，tlinerf为射频信号线缆时延值。

其中，重复执行步骤(2)至步骤(6)，获得多个时延值，对所有时延值求平均值，将该平均值作为td-altboc信号的最终时延值。

本发明相对背景技术的优点在于：

(1)本发明方法率先提出了一种针对卫星导航td-altboc信号的时延标定技术手段，为该领域填充了空白。

(2)本发明方法中的时延标定方法能够在采集一次信号的条件下，对各个支路的时延值进行全部标定。

附图说明

图1为本发明方法所涉及设备的连接图。

图2为本发明方法的实现流程图。

图3为本发明的td-altboc信号等效扩频码生成示意图。

图4为本发明方法所提出的相关峰搜索时延标定原理图。

图5为本发明的链路时延值与实际发射时延值之间的关系图。

具体实施方式

下面结合附图对本发明做进一步详细的说明。

如图1所示，本发明方法涉及的td-altboc信号时延标定系统包括td-altboc信号发射设备、时频信号产生设备、高速信号采集装置、时延处理装置、触发信号线缆以及射频信号线缆。

如图2所示，本发明方法的主要步骤如下：

(1)搭建td-altboc信号时延标定系统：按照图1所示连接td-altboc信号发射设备、时频信号产生设备、高速信号采集装置、时延处理装置、触发信号线缆以及射频信号线缆。

在工作状态下，时频信号产生设备为td-altboc信号发射设备和高速信号采集装置提供10mhz信号；td-altboc信号发射设备为高速信号采集装置提供1pps触发信号和td-altboc射频信号；时延处理装置控制高速信号采集装置采集td-altboc信号数据，并对数据进行分析和处理。在开始测试前，通过矢量网络分析仪或其它设备测试触发信号线缆以及射频信号线缆的时延值，分别记为tline1pps以及tlinerf。

(2)使用高速信号采集装置采集并存储td-altboc信号；所述的td-altboc信号包括上边带信号和下边带信号。

采集信号过程须满足以下两个方面要求：①采集信号起始时刻为1pps上升沿触发时刻；②高速信号采集装置的采样率应足够高，使得采样时间间隔大于标定精度要求。建议在实际使用的过程中，设置采集装置以1pps作为触发信号，采样率为10ghz，采集并存储一个伪码周期的td-altboc信号。在这样的设置条件下，既保证了时延标定精度，又减小了数据文件存储大小，减小了后续的计算量。

(3)根据td-altboc信号4个支路的扩频码和时隙分配生成上、下边带信号各支路的等效扩频码。

td-altboc信号上、下边带信号的扩频码都由数据支路扩频码和导频支路扩频码按时隙交替合成，其中奇数位置处为数据支路扩频码片，偶数位置处为导频支路扩频码片。据此规律，可按如图3所示，在原数据支路扩频码的每个扩频码片之后插入一个零值码片生成数据支路等效扩频码；同理，在原导频支路扩频码的每个扩频码片之前插入一个零值码片生成导频支路等效扩频码。

(4)将td-altboc信号中的上边带信号和下边带信号分别进行下变频一一对应得到上边带基带信号和下边带基带信号；将td-altboc上、下边带各支路等效扩频码以高速信号采集装置采样率作采样处理，采样后将上边带基带信号和下边带基带信号分别与上边带各支路等效扩频码和下边带各支路等效扩频码一一对应作滑动相关运算得到滑动相关运算结果，并搜索得到各支路最大相关峰值；

以高速信号采集装置采样率作采样处理，采用数字信号处理的方法对采集到的td-altboc信号进行下变频运算，下变频过程中对各信号初始相位进行搜索，以保证得到的相关峰值最大。

设td-altboc信号数据为rf(i),i＝1,2,...,n，td-altboc信号上边带信号中心频率为f1，下边带信号中心频率为f2，采样率为fs，则td-altboc上边带信号下变频得到的基带信号为：

td-altboc下边带信号下变频得到的基带信号为：

在公式(1)、(2)中的δθ1和δθ2分别是所采集td-altboc信号中td-altboc上边带信号初始相位以及td-altboc下边带信号初始相位。δθ1和δθ2的具体值通过循环执行步骤(4)、步骤(5)来确定，每次执行步骤(4)，使δθ1和δθ2的值递增，使得步骤(5)中相关峰值最高时的δθ1和δθ2为信号的初始相位。

步骤(4)具体包括以下步骤：

(401)分别递增上边带信号和下边带信号的下变频相位，并分别根据递增后的下变频相位一一对应对上边带信号和下边带信号进行下变频得到上边带基带信号和下边带基带信号；

(402)对上边带基带信号和下边带基带信号分别与各支路等效扩频码一一对应进行滑动相关运算，记录相关峰值大小，返回步骤(401)，直到下变频相位递增满2π；

(403)搜索得到记录的最大相关峰值。

相关运算过程如图4所示。图中的基带信号数据可用下式表示：

设高速信号采集装置采样率为fs，扩频码周期为t，扩频码码片序列为prn，码速率为fc，n为基带信号的点数，x(i)的长度为n，则图中所示的扩频码数据prn(i)的生成方法可用公式(4)表示：

其中，符号为向上取整，为向下取整。

滑动相关计算方法如公式(5)所示：

(5)根据滑动相关运算结果r(τ)，分别搜索各支路最大相关峰值所在位置，计算得到td-altboc信号上、下边带各个支路链路时延值。

设在滑动相关运算结果r(τ)中，最大值出现的位置序号为pos_max，则相关峰搜索计算得到td-altboc信号的链路时延值td_calc如公式(6)所示：

由此计算得到的发射设备的链路时延值td_calc单位为秒。

(6)用链路时延值扣除触发信号线缆的时延值和射频信号线缆的时延值后，得到td-altboc信号的时延值。

如图5所示，td-altboc信号实际时延td需要按照公式(7)扣除步骤(1)中获得的td-altboc信号发射设备与高速信号采集装置相连的1pps触发信号线缆tline1pps以及td-altboc射频信号线缆tlinerf的时延值获得：

td＝td_calc+tline1pps-tlinerf公式(7)

上述td-altboc信号时延标定方法只分析了在一个时刻采集到的td-altboc信号，所得到的时延值没有考虑时延的统计特性，在实际使用过程中，需多次重复计算求均值。

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