一种适用于BOC调制的卫星导航授权信号估计方法

一种适用于boc调制的卫星导航授权信号估计方法
技术领域
1.本发明属于卫星导航信号扩频码盲估计技术领域，特别是涉及一种适用于boc调制的卫星导航授权信号估计方法。


背景技术：

2.随着卫星导航相应技术的迅猛发展，卫星导航系统已经广泛应用于社会的各行各业，随着gnss的广泛应用和数字信号处理技术的进步，基于先进集成电路的接收机不断接近导航信号的性能的性能极限，卫星导航传统民用公开频点信号的局限性也逐渐显现。二进制偏移载波（binary offset carrier, boc）扩频调制技术做为新一代gnss信号的典型新技术，使用了副载波的扩频调制方式，一方面在同一个中心频点约束下，最大程度实现了与旧导航信号的频谱分离；另一方面，更宽的频谱带宽提高了新信号的跟踪精度上限。以gps新信号为例，l1m信号由于采用了boc(10,5)调制，其相关峰在正负一个码片范围内会有多个峰值，而主峰相比于公开信号l1c/a信号的相关峰具备更加陡峭的边沿，从而测距偏差表现出更为灵敏的变化，在测距精度上具备更大的优势。
3.目前针对boc调制的授权码估计技术主要基于单边带滤波后的接收处理，在利用相应的大口径天线高增益条件下，通过上或下边带滤波得到类似传统bpsk调制的信号，从而完成相应码流观测量提取。该方法利用单边带滤波避免了boc调制的影响，但是损失了一半的有效信号能量，影响了授权码估计性能的潜在上限。


技术实现要素：

4.针对以上技术问题，本发明提供一种适用于boc调制的卫星导航授权信号估计方法。
5.本发明解决其技术问题采用的技术方案是：一种适用于boc调制的卫星导航授权信号估计方法，方法包括以下步骤：步骤s100：采集导航信号，对导航信号中同频段的公开信号进行捕获和跟踪，得到公开信号的码相位、载波相位和多普勒测量值；步骤s200：根据公开信号的载波相位和多普勒测量值对公开信号进行载波剥离，对导航信号中同频段的boc调制授权信号的双边带主瓣进行带通滤波；步骤s300：根据公开信号码相位测量值找到公开信号码相位起始位置点，截取带通滤波后的boc调制授权信号采样点，使授权信号的起始采样点与公开信号的码相位起始位置点对齐；步骤s400：剥离boc调制授权信号的副载波，对剥离副载波后的boc调制授权信号进行主瓣低通滤波得到基带iq复信号；步骤s500：将基带iq复信号加窗去掉码片边缘的采样点，对加窗内复信号的采样点求平均得到每个授权码的有效采样点，将所有有效采样点在复数平面上进行层次聚类分析得到两个簇，根据两个簇与预设的簇-授权码
±
1的对应关系，得到授权码流的估计结果。
6.优选地，步骤s200中根据公开信号的载波相位和多普勒测量值对公开信号进行载波剥离，包括：根据公开信号的载波相位和多普勒测量值对公开信号进行数字下变频至零频实现载波剥离。
7.优选地，数字下变频具体为：其中，表示采样信号，表示频段的中心频率，为公开信号的多普勒测量值，为载波相位测量值。
8.优选地，步骤s400中剥离boc调制授权信号的副载波具体为：其中，为剥离副载波后的信号，为对齐公开信号码相位起始位置点的授权信号，是符号函数，表示副载波频率。
9.优选地，步骤s400中的低通滤波实现方式包括快速傅里叶变换滤波器、有限长单位脉冲响应滤波器和无限长单位脉冲响应滤波器。
10.优选地，步骤s200中双边带主瓣进行带通滤波的实现方式为快速傅里叶变换滤波器。
11.优选地，步骤s500中将基带iq复信号加窗去掉码片边缘的采样点，对加窗内复信号的采样点求平均得到每个授权码的有效采样点，将所有有效采样点在复数平面上进行层次聚类分析得到两个簇，包括：取基带iq复信号码片正中间的一半宽度内采样点进行加窗，对加窗内复信号的采样点求平均得到得到每个授权码的有效采样点；使用ward法对所有有效采样点在复数平面上进行层次聚类分析得到两个簇。
12.优选地，步骤s520包括：步骤s521：将每个有效采样点当成一个初始簇；步骤s522：使用ward法计算任意两个簇的距离，并找到最近的两个簇；步骤s523：合并距离最近的两个簇，生成新的簇的集合；步骤s524：返回步骤s522直到最后簇的个数为2。
13.优选地，步骤s522具体为：其中，表示簇r和簇s的ward距离，表示欧几里德距离，
和分别表示簇r和簇s的中心点，和分别表示簇r和簇s的数据点数。
14.相对于现有技术，本发明考虑到boc调制对信号频谱的影响，采用信号对齐后boc解调和聚类分析的的方法，提出了一种适用于boc调制的卫星导航授权信号估计方法，通过对boc调制授权信号的双边带滤波和boc解调，充分利用信号能量，构建聚类分析结果与调制星座图的映射关系，估计出授权信号的码流信息，具有方法简单可靠，易于实现等优点。
附图说明
15.图1为本发明的一种适用于boc调制的卫星导航授权信号估计方法流程图；图2为本发明一实施例中提供的数字下变频后的频谱图；图3为本发明一实施例中授权码信号双边带主瓣带通滤波后的频谱图；图4为本发明一实施例中授权码信号剥离副载波后的频谱图；图5为本发明一实施例中低通滤波后的频谱图；图6为本发明一实施例中低通滤波后的调制星座图；图7为本发明一实施例中低通滤波后的加窗示意图；其中，（a）为加窗前波形图；（b）为加窗波形图；（c）为加窗后波形图；图8为本发明一实施例中加窗求平均后的调制星座图；图9为本发明一实施例中授权码正确估计后的信号相关峰。
具体实施方式
16.为了使本技术领域的人员更好地理解本发明的技术方案，下面结合附图对本发明作进一步的详细说明。
17.在一个实施例中，如图1所示，一种适用于boc调制的卫星导航授权信号估计方法，方法包括以下步骤：步骤s100：采集导航信号，对导航信号中同频段的公开信号进行捕获和跟踪，得到公开信号的码相位、载波相位和多普勒测量值。
18.具体地，抛物面高增益天线对准指定的gnss卫星后，信号采集器可实时采集高载噪比的gnss信号，该采样信号可以是数字中频信号或射频直采信号，量化位数通常在8bit以上，采样数字信号可以是实信号，也可以是复信号。在本例中，信号的采样率是1gsps，量化位数12bit，实信号射频直采，覆盖了整个l1频段。在采集导航信号后，对同频段公开信号捕获与跟踪，获得公开信号的码相位、载波相位和多普勒测量值。
19.步骤s200：根据公开信号的载波相位和多普勒测量值对公开信号进行载波剥离，对导航信号中同频段的boc调制授权信号的双边带主瓣进行带通滤波。
20.在一个实施例中，步骤s200中根据公开信号的载波相位和多普勒测量值对公开信号进行载波剥离，包括：根据公开信号的载波相位和多普勒测量值对公开信号进行数字下变频至零频实现载波剥离。
21.在一个实施例中，数字下变频具体为：
其中，表示采样信号，表示频段的中心频率，为公开信号的多普勒测量值，为载波相位测量值。
22.进一步地，数字下变频后的频谱图如图2所示。
23.在一个实施例中，步骤s200中双边带主瓣进行带通滤波的实现方式为快速傅里叶变换滤波器。
24.具体地，对授权码信号l1m的两个主瓣进行双边带通滤波，本例中选用1ms采样点数据，即1000000点快速傅里叶变换（fft）实现双边带中心频点分别为
±
10.23mhz，带宽都为10.23mhz的带通滤波，带通滤波后记为，其频谱图如图3所示。带通滤波后保留了双边带boc调制的授权信号，带外抑制大于40db，能够有效抑制副载波剥离引入的信号混叠。
25.步骤s300：根据公开信号的码相位测量值找到公开信号码相位起始位置点，截取带通滤波后的boc调制授权信号采样点，使授权信号的起始采样点与公开信号的码相位起始位置点对齐。
26.步骤s400：剥离boc调制授权信号的副载波，对剥离副载波后的boc调制授权信号进行主瓣低通滤波得到基带iq复信号。
27.在一个实施例中，步骤s400中剥离boc调制授权信号的副载波具体为：其中，为剥离副载波后的信号，为对齐公开信号码相位起始位置点的授权信号，是符号函数，表示副载波频率。
28.具体地，以l1m信号为例，。剥离载波后的信号频谱如图4所示，将双边带信号全部集中在基带上，信号功率有近3db的提升。
29.在一个实施例中，步骤s400中的低通滤波实现方式包括快速傅里叶变换滤波器、有限长单位脉冲响应滤波器和无限长单位脉冲响应滤波器。
30.具体地，为抑制剥离副载波导致的高频分量影响，采用主瓣低通滤波，本例中数字低通滤波的截止频率为5.115mhz，低通滤波器的实现可采用fft滤波器、fir滤波器、iir滤波器等多种实现方式，本例中选用1ms采样点数据，即1000000点fft实现通带范围为-5.115~ 5.115mhz的频率低通滤波，滤波后的信号频谱如图5所示，仅保留授权码主瓣带宽内信号，带外抑制大于40db，有效抑制了带外噪声对授权码估计的影响。
31.步骤s500：将基带iq复信号加窗去掉码片边缘的采样点，对加窗内复信号的采样点求平均得到每个授权码的有效采样点，将所有有效采样点在复数平面上进行层次聚类分析得到两个簇，根据两个簇与预设的簇-授权码
±
1的对应关系，得到授权码流的估计结果。
32.在一个实施例中，步骤s500中将基带iq复信号加窗去掉码片边缘的采样点，对加窗内复信号的采样点求平均得到每个授权码的有效采样点，将所有有效采样点在复数平面上进行层次聚类分析得到两个簇，包括：取基带iq复信号码片正中间的一半宽度内采样点进行加窗，对加窗内复信号的采样点求平均得到得到每个授权码的有效采样点；使用ward法对所有有效采样点在复数平面上进行层次聚类分析得到两个簇。
33.具体地，低通滤波后，信号调制星座图如图6所示，基带iq复信号对应为码片，码片在翻转点的变化趋于平缓，不利于码片正负的正确判决，为降低这种影响，采用码片中间部分采样点求平均的方法。如图7所示，从两个码片波形随时间变化，在两个码片交界处取值，容易导致码片值判断错误，因此可以通过加窗方式减少码片交界处的影响，仅取加窗内的数值代表码片内有效数值，图7中可以看出加窗后，使得码片内的有效值更加远离0值，有利于后续正确估计授权码值。
34.本例中取码片正中间一半码片宽度内采样点求平均。加窗求平均后的调制星座图如图8所示，对比图6加窗前的结果，显然图8更有利于信号的分类。
35.进一步地，将码片内加窗求平均结果的复信号的实部和虚部构成一个矢量，通过使用ward法计算簇之间的最短距离，创建凝聚层次树结构。ward距离计算使用递增的平方和，也就是说，由于连接两个簇，簇内总平方和的增加。簇内平方和定义为簇内所有对象与簇质心之间距离的平方和。
36.在一个实施例中，步骤s520包括：步骤s521：将每个有效采样点当成一个初始簇；步骤s522：使用ward法计算任意两个簇的距离，并找到最近的两个簇；步骤s523：合并距离最近的两个簇，生成新的簇的集合；步骤s524：返回步骤s522直到最后簇的个数为2。
37.在一个实施例中，步骤s522具体为：其中，表示簇r和簇s的ward距离，表示欧几里德距离，和分别表示簇r和簇s的中心点，和分别表示簇r和簇s的数据点数。
38.每个簇中包括多个授权码的有效采样点，根据两个簇与预设的簇-授权码
±
1的对应关系，得到授权码流的估计结果。得到授权码流后，通过对信号的跟踪处理可以得到授权信号的相关峰如图9所示，相关峰符合boc调制的特点，验证了该授权码估计方法的正确性。
39.相对于现有技术，本发明考虑到boc调制对信号频谱的影响，采用信号对齐后boc解调和聚类分析的的方法，提出了一种适用于boc调制的卫星导航授权信号估计方法，通过对高增益信号的双边带滤波和boc解调，充分利用信号能量，构建聚类分析结果与调制星座图的映射关系，估计出授权信号的码流信息，具有方法简单可靠，易于实现等优点。
40.以上对本发明所提供的一种适用于boc调制的卫星导航授权信号估计方法进行了详细介绍。本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的核心思想。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。