一种卫星导航信号各分量间载波相位关系检测方法
【技术领域】
[0001] 本发明涉及一种卫星导航信号各分量间载波相位关系检测方法，属于卫星导航技 术领域。
【背景技术】
[0002] 在导航卫星信号中，为满足不同服务需求，每个频点上加载了多个信号分量。例如 最早的GPS L1上为C/A码和P码信号，为QPSK调制复合方式，只有两路信号分量，而现代化后 的GPS L1上有C/A、P、M、L1C四路信号，Galileo系统、BDS等很多频点上信号分量都多于三 路，并使用了 InterpleX、CASM、POCET等较为复杂的调制复合方式。
[0003] 目前，可以通过以下方法来进行导航信号载波相位关系的检测：
[0004] 方法1，矢量信号分析仪；
[0005] 方法2,接收机跟踪；
[0006] 然而，上述方法存在如下问题：
[0007] (1)方法1使用通用仪器，这种方法仅能对QPSK、8PSK等标准数字调制信号的载波 相位关系的正交性等给出结果。而调制复合方式复杂的卫星导航无法给出结果。
[0008] (2)方法2是较常用的导航信号载波相位检测方法，但是闭环跟踪接收机仅能给出 正交、翻转等信息，无法给出复杂的诸如AltBOC这种含有异频分量的卫星导航信号载波相 位关系。
【发明内容】
[0009] 本发明的技术解决问题是：克服现有技术的不足，提供了一种卫星导航信号各分 量间载波相位关系检测的方法，可以准确判明多分量卫星导航信号的同频和异频分量载波 相位相对关系。
[0010]本发明的技术解决方案是：
[0011] -种卫星导航信号各分量间载波相位关系检测方法，步骤如下：
[0012] (1)采集导航信号：与导航信号生成源同源采集导航信号并进行存储，如不能同 源，则采集后剥离多普勒影响；
[0013] (2)选取最易获得起始时间信息的信号分量作为引导信号，通过频域搜索的方式 初步确定该信号载波相位和码相位，之后再进行时域搜索，获取引导信号最终的码相位和 载波相位；
[0014] (3)移动数据读取指针至码相位起始处，开环解调数据，获取该频点所有信号的共 同的起始时间信息；
[0015] (4)利用所述起始时间信息生成所有信号分量伪码；
[0016] (5)根据步骤(4)中生成的所有信号分量伪码，计算信号最终的载波相位记录值；
[0017] (6)各路信号分量进行开环解调，剥离各级编码，得到帧头或电文，判明电文翻转 状态；
[0018] (7)载波相位校正，利用步骤(6)中的电文翻转状态校正步骤(5)中的载波相位统 计值，从而得到各分量间载波相位关系。
[0019] 步骤(1)中采集导航信号时，采集时间大于帧长，采样ADC位数大于等于8bit，采样 速率的范围通过如下方式确定：
[0020] (2 · 1)导航信号频率为L波段的XI: 1575 · 42MHz 土BW1，X2:1191 · 795MHz 土BW2，X3: 1268 · 52MHz 土 BW3，其中，BW1、BW2和BW3均为单边带宽，BW1 >20 · 46MHz，BW2>35MHz，BW3> 20.46MHz;
[0021] (2.2)根据(2.1)中选定的导航信号频率和带宽，先依带通采样定理确定每个频点 不混叠采样频率范围，再选取重叠部分；
[0022] (2.3)根据下表确定导航信号的频谱和相位翻转状态以及采样速率范围：
[0023]表1相位翻转与采样频率选择
[0026] 其中频谱和相位翻转中的1表示未翻转、-1表示翻转。
[0027] 所述步骤(2)中进行频域搜索时，根据步骤(1)中相位翻转状态设置本地载波频率 符号，当该采样速率下信号频谱和相位翻转为-1时，本地信号中心频率符号取反。
[0028] 所述步骤(2)中进行时域搜索具体为：
[0029] (4.1)设定粗略初相为Φ 0，初始搜索范围为Φ ±0.1，搜索步长为0.0 lrad;
[0030] (4.2)将移动和采样的本地码与剥离载波的信号进行相关，求取最大值处的载波 相位Φ 1;
[0031] (4.3)以巾1为初相，搜索范围为巾±0.01，搜索步长为0.001瓜(1;
[0032] (4.4)迭代执行步骤(4.2)~(4.3)，直至载波相位搜索时码相关最大值大于等于 两个，结束迭代，此时求出的最大值处的载波相位即为引导信号最终的载波相位。
[0033] 所述步骤(3)开环解调数据，通过以下方法实现：
[0034] (5.1)偏移数据文件起点使数据从码相位起点开始读取，以一个数据bit或者子码 bit的时间作为固定的积分时间，根据所述固定的积分时间进行数据读取并利用载波相位 进行相关运算，提取相关峰的符号；
[0035] (5.2)在相关峰符号基础上搜索子码相位，剥离子码并进行子帧同步，获取该频点 所有信号的共同的起始时间信息。
[0036]所述步骤(5)通过以下方法实现：
[0037] 1)当信号分量伪码相位存在超过1个样点的以上偏差时，将计算载波相位的位置 选在一个样点处；
[0038] 2)本地生成完整的恒包络基带信号，以恒包络基带信号作为整体的一个伪码与信 号进行相关，求取该合路信号伪码相位，以此相位位置作为计算载波相位的时刻；
[0039] 3)根据步骤2)中得到的计算载波相位的时刻作为起始时间信息，通过频域搜索的 方式初步确定信号载波相位和码相位；
[0040] 4)再进行时域搜索，获取信号最终的码相位和载波相位，从而得到信号最终的载 波相位记录值。
[0041] 步骤(7)所述的载波相位校正，利用步骤(6)电文翻转状态结果校正步骤(5)中的 载波相位统计值，以下表为准进行校正：
[0042] 表2校正表
[0044] 若载波相位计算值大于2*π，则减去2*π后为最终的载波相位计算值。
[0045] 本发明与现有技术相比的有益效果是：
[0046] (1)
[0047] 本发明考虑了因采集、变频造成的相位翻转对载波相位的影响，规定了导航信号 该种性能检测时必需的采样频率的取值范围；通过明确频谱内容和相位翻转情况，可以准 确给出载波相位关系符号。
[0048] (2)本发明可以检测含长码信号分量的导航信号，包括因二级编码引入的长时间 周期码信号分量。通过解调引导信号，可以通过较短时间的采集数据得到电文翻转状态。
[0049] (3)本发明可以对含异频信号分量的导航信号各分量间载波相位关系进行检测， 包括A1 tBOC信号以及对频间信号时延有约束的多频段信号。
[0050] (4)本发明对载波相位评估精度较高，可达固定积分时间和采样频率下的极限。
【附图说明】
[0051 ]图1为本发明卫星导航信号各分量间载波相位关系检测方法流程图；
[0052]图2为码相位和载波相位粗搜索即频域搜索信号流图；
[0053]图3为精细码相位搜索示意图；
[0054]图4为信号分量相关峰不在一处的情况不意图；
[0055] 图5为GPS L1信号分量载波相位关系仿真示意图；
【具体实施方式】
[0056] 本发明提供了一种卫星导航信号各分量间载波相位关系检测的方法，可以准确判 明多分量卫星导航信号的同频和异频分量载波相位相对关系。当代BDS、GPS、Galileo系统、 GL0NASS、QZSS卫星导航系统及导航信号模拟源生成的导航信号都具有一个特点：一个码周 期对应整数个载波周期。这也是这种检测方法应用的基本条件。
[0057]如图1所示，本发明提供的载波相位关系检测方法包括以下步骤：
[0058] (1)导航信号采集、存储。
[0059]与导航信号生成源同源采集导航信号并进行存储，如不能同源，需要采集后剥离 多普勒影响;采集时间大于帧长;采样ADC位数大于等于8bit，采样速率由信号中心频率、带 宽决定，采样速率的范围通过如下三个小步骤实现：
[0060] ①导航信号频率和带宽选定。GPS、Galileo系统、GLONASS、QZSS、RNSS主要为L波段 的1575.4210^±8￥1(以乂1表示，8￥1为其单边带宽）、1191.7951他±8￥2(以乂2表示，8￥2为其 单边带宽）、1268.52MHz 土BW3(以X3表示，BW3为其单边带宽），每个频点信号带宽都在80MHz 以内。考虑到调制方式并留有一定余量，BW1 >20.46MHz，BW2>35MHz，BW3>20.46MHz。
[0061] ②不混叠采样频率选定。根据①中选定的导航信号频率和带宽，先依带通采样定 理确定每个频点不混叠采样频率范围，再选取重叠部分。
[0062] ③翻转状态确定。采样频率选取范围确定后，对于位于第1、3、5、7···等奇次奈奎斯 特区域的信号，信号上下边带频谱内容和载波相位不会发生翻转。对于位于第2、4、6、8···等 偶次奈奎斯特区域的信号，信号上下边带频谱内容和载波相位会发生翻转。对得到采样频 率选取范围、对应频谱和相位翻转状态如下表所示。其中频谱和相位翻转的1表示未翻转、-1表示翻转，表中同时给出了采样频率选择的典型值。
[0063]表1相位翻转与采样频率选择
[0065] 实际中根据硬件条件选取尽可能高的采样频率，且采样频率不能为码速率、副载 波速率的整数倍。
[0066] (2)引导信号码相位、载波相位搜索。选取最易获得起始时间信息的信号分量作为 引导信号，一般可以选择短周期民码信号，这些信号周期较短，容易捕获和解调电文，如GPS L1C/A码信号、BDS B11信号。使用频域粗略搜索该信号载波相位和码相位，之后在时域精细 搜索，获取精确的码相位和载波相位。
[0067] 首先使用频域方法粗搜索码相位和载波相位，如图2所示。该步码相位搜索精度到 达样点级。此步需要根据步骤（1)中翻转状态设置本地载波频率符号。当（1)中该采样频率 下信号频谱和相位翻转为-1时，本地信号中心频率符号取反。
[0068] 粗搜索后进行细搜索，如图