CN114527499B - 用于卫星导航系统b1宽带复合信号的无模糊定位方法及其装置 - Google Patents

用于卫星导航系统b1宽带复合信号的无模糊定位方法及其装置 Download PDF

Info

Publication number
CN114527499B
CN114527499B CN202210156560.0A CN202210156560A CN114527499B CN 114527499 B CN114527499 B CN 114527499B CN 202210156560 A CN202210156560 A CN 202210156560A CN 114527499 B CN114527499 B CN 114527499B
Authority
CN
China
Prior art keywords
signal
propagation delay
code
ambiguity
determining
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210156560.0A
Other languages
English (en)
Other versions
CN114527499A (zh
Inventor
姚铮
漆耘含
陆明泉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN202210156560.0A priority Critical patent/CN114527499B/zh
Publication of CN114527499A publication Critical patent/CN114527499A/zh
Priority to PCT/CN2023/077393 priority patent/WO2023155920A1/zh
Application granted granted Critical
Publication of CN114527499B publication Critical patent/CN114527499B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

提供了一种用于卫星导航系统B1宽带复合信号的无模糊定位方法及其装置，能够通过对B1宽带复合信号中B1C信号和B1I信号进行跟踪，确定多颗卫星的观测量，并通过所述观测量确定各颗卫星到接收机的伪距值和码伪距值，利用所述伪距值和所述码伪距值求解B1宽带复合信号的传播延迟的模糊度浮点解，并通过模糊度浮点解得到模糊度整数解，以修正被估计错误的B1宽带复合信号的传播延迟模糊度，从而得到无模糊的伪距值，实现对接收机所处位置的无模糊且高精度的定位。

Description

用于卫星导航系统B1宽带复合信号的无模糊定位方法及其 装置
技术领域
本申请实施方式涉及卫星导航定位领域。具体地，提供了一种用于卫星导航系统B1宽带复合信号的无模糊定位方法及其装置。
背景技术
目前，基于位置信息的服务产业的发展对全球导航卫星系统(GNSS)定位的精度要求越来越高。为了满足不断出现的新应用对高精度定位的需求，新一代的卫星导航信号被设计，以提高测距精度的潜在能力上限。
新一代GNSS卫星导航信号普遍引入了副载波调制。与传统GNSS中的二进制相移键控(BPSK)信号相比，副载波调制的卫星导航信号，比如使用了方波副载波的BOC调制信号、使用了复用副载波的MBOC调制信号、使用了双边带副载波的AltBOC调制信号、非对称恒包络ACE-BOC调制信号以及单边带SCBOC调制信号等，具有更好的频域分离特性和更宽的Gabor带宽，从而拥有更高精度的测距潜力。在通过信号升级提升了测距精度的潜在能力上限的同时，还需要相应的先进的接收处理算法将升级后的信号中蕴含的潜在高性能实际发挥出来。一方面，接收技术需要能够高精度地跟踪信号，另一方面，由于副载波调制信号多峰的ACF为码跟踪引入了模糊度威胁，导致定位结果的不可靠，因此接收技术还需要能消除模糊威胁，实现可靠的跟踪和定位。
第三代北斗导航卫星系统(BDS3)B1频点通过CEMIC技术在BDS3卫星上播发的宽带复合信号由公开服务信号B1C和B1I组成，其中，为了与第二代北斗导航卫星系统(BDS2)B1I信号兼容，BDS3B1I信号采用了SCBOC(14,2)调制的方式。因为SCBOC(14,2)调制信号属于高阶BOC信号，所以BDS3 B1宽带复合信号具有极高精度的测距潜力。
目前，在对信号跟踪的过程中，卫星导航信号需要与本地载波相乘得到基带信号，基带信号再与本地码相乘并积分从而得到相关值，采用了SCBOC(14,2)调制的B1I信号的副载波相位与载波相位都在cos函数中，是严重耦合的，传统的跟踪技术难以在只对SCBOC(14,2)调制信号跟踪的情况下同时准确估计副载波相位和载波相位。现有的交叉辅助跟踪算法利用了BDS3 B1C和BDS3 B1I两个信号互相辅助跟踪的方式解开了SCBOC(14,2)副载波相位和载波相位严重耦合的问题，但该技术没有解决模糊跟踪问题，其测距和定位结果是不可靠的。
此外，SCBOC(14,2)信号的自相关函数是多峰的，这会为信号的跟踪带来严重的模糊跟踪问题。目前，解决模糊跟踪问题的技术主要分为三大类：第一类是一维无模糊跟踪技术，比如Bump Jump、BPSK-like、PUDLL，该类技术对低阶BOC非常有效，但当BOC信号阶数升高时，这类技术可能会失效；第二类是二维无模糊跟踪技术，比如针对简单BOC信号通用的DET、DPE、以及针对具有复杂结构的AltBOC的跟踪算法，与一维无模糊跟踪技术相比，二维无模糊跟踪技术更具有一般性并且能提供更可靠的跟踪结果，但当二维跟踪的信号属于高阶BOC信号、低载噪比或受到严重多径影响时，该技术仍然会发生模糊跟踪问题；第三类是基于参数估计的技术，该类方法是在二维无模糊跟踪技术的基础上，进一步修正出现模糊问题的副载波观测量，使得测距和定位结果更加可靠，但该方法求解的副载波模糊度误差较大，无法保证修复所有出现模糊问题的副载波观测量。
目前尚没有一种处理技术能够保证SCBOC(14,2)的无模糊测距和定位。
发明内容
本申请的实施方式提供了可至少部分解决相关技术中存在的上述问题、或其他问题的用于卫星导航系统B1宽带复合信号的无模糊定位方法及其装置。
本申请提供了一种用于卫星导航系统B1宽带复合信号的无模糊定位方法及其装置，能够通过对B1宽带复合信号中B1C信号和B1I信号进行跟踪，确定多颗卫星的观测量，并通过所述观测量确定各颗卫星到接收机的伪距值和码伪距值，利用所述伪距值和所述码伪距值求解B1宽带复合信号的传播延迟的模糊度浮点解，并通过模糊度浮点解得到模糊度整数解，以修正被估计错误的B1宽带复合信号的传播延迟模糊度，从而得到无模糊的伪距值，实现对接收机所处位置的无模糊且高精度的定位。
本申请一方面提供了一种用于卫星导航系统B1宽带复合信号的无模糊定位方法，包括：通过接收机接收多颗卫星的B1宽带复合信号，其中，所接收的B1宽带复合信号包括B1I信号和B1C信号；获取所述多颗卫星的所述B1I信号和所述B1C信号的观测量；基于所述观测量，确定各颗所述卫星到所述接收机的伪距值和码伪距值；基于所述伪距值和所述码伪距值，确定所述B1宽带复合信号的传播延迟的模糊度浮点解；通过所述模糊度整数解更正所述B1宽带复合信号的传播延迟的模糊度，获取无模糊的B1宽带复合信号的传播延迟，并基于所述无模糊的B1宽带复合信号的传播延迟，对所述接收机所处位置进行无模糊的定位。
在本申请一实施方式中，所述观测量包括：所述B1C信号的载波相位、B1I信号的码传播延迟、以及B1I信号的副载波传播延迟。
在本申请一实施方式中，通过合并所述码传播延迟和所述副载波传播延迟，确定所述B1I信号的传播延迟，基于所述B1I信号的传播延迟确定所述伪距值，以及基于所述码传播延迟确定所述码伪距值。
在本申请一实施方式中，获取所述B1I信号和所述B1C信号的观测量的步骤包括：基于所述B1C信号、所述B1I信号的码传播延迟、以及载波跟踪环路，确定所述B1C信号的载波频率和所述B1C信号的载波相位；基于所述B1C信号的载波频率和所述B1C信号的载波相位，确定所述B1C信号的载波相位；以及基于所述B1I信号的码跟踪环路确定所述B1I信号的码传播延迟。
在本申请一实施方式中，获取所述B1I信号和所述B1C信号的观测量的步骤包括：基于所述B1I信号、所述B1C信号的载波频率和相位、以及所述B1I信号的码跟踪环路，确定所述B1I信号的扩频码频率和相位；基于所述B1I信号的扩频码频率和相位，确定所述B1I信号的码传播延迟；基于所述B1I信号、所述B1C信号的载波频率和相位、以及副载波跟踪环路，确定所述B1I信号的副载波频率和相位；以及基于所述B1I信号的副载波频率和相位，确定所述B1I信号的副载波传播延迟。
在本申请一实施方式中，基于所述B1I信号和所述B1C信号的观测量，确定所述多颗卫星中各颗卫星到所述接收机的伪距值和码伪距值的步骤包括：通过所述B1C信号的载波相位和所述B1I信号的码传播延迟确定所述码伪距值；以及通过所述B1I信号的副载波传播延迟和所述B1I信号的码传播延迟确定所述伪距值。
在本申请一实施方式中，基于所述B1I信号和所述B1C信号的观测量，确定所述多颗卫星中各颗卫星到所述接收机的伪距值和码伪距值的步骤包括：基于所述B1C信号的载波相位，对所述码传播延迟进行平滑；基于所述B1I信号的传播延迟和平滑后的所述码传播延迟，确定所述多颗卫星中各颗卫星到所述接收机的所述伪距值和所述码伪距值。
在本申请一实施方式中，基于所述模糊度浮点解确定所述模糊度整数解的步骤包括：通过LAMBDA算法，基于所述模糊度浮点解确定所述模糊度整数解。
本申请另一方面提供了一种用于卫星导航系统B1宽带复合信号的无模糊定位装置，包括：存储器，存储有计算机可执行的指令；处理器，执行所述指令以执行如上所述的任一方法。
本申请另一方面提供了一种存储介质，存储有计算机可执行的指令，所述指令被一个或多个处理器执行时实现如上所述的任一方法。
附图说明
为了更全面地理解本发明及其优点，现在结合附图参考以下描述，其中相同的附图标记表示相同的部分：
图1示出了根据本申请实施例的用于卫星导航系统B1宽带复合信号的无模糊定位方法的流程图；
图2示出了根据本申请实施例的获取B1C信号的观测量的流程图；
图3示出了根据本申请实施例的获取B1I信号的观测量的流程图；
图4示出了根据本申请实施例的基于B1I信号和B1C信号的观测量确定卫星到接收机的伪距值和码伪距值的流程图；
图5示出了根据本申请实施例的用于卫星导航系统B1宽带复合信号的无模糊定位装置的示意图；
图6示出了根据本申请实施例的跟踪模块和模糊度解算模块的示意图。
具体实施方式
为了更好地理解本申请，将参考附图对本申请的各个方面做出更详细的说明。应理解，这些详细说明只是对本申请的示例性实施方式的描述，而非以任何方式限制本申请的范围。在说明书全文中，相同的附图标号指代相同的元件。表述“和/或”包括相关联的所列项目中的一个或多个的任何和全部组合。
还应理解的是，诸如“包括”、“包括有”、“具有”、“包含”和/或“包含有”等表述在本说明书中是开放性而非封闭性的表述，其表示存在所陈述的特征、元件和/或部件，但不排除一个或多个其它特征、元件、部件和/或它们的组合的存在。此外，当诸如“...中的至少一个”的表述出现在所列特征的列表之后时，其修饰整列特征，而非仅仅修饰列表中的单独元件。此外，当描述本申请的实施方式时，使用“可”表示“本申请的一个或多个实施方式”。并且，用语“示例性的”旨在指代示例或举例说明。
除非另外限定，否则本文中使用的所有措辞(包括工程术语和科技术语)均具有与本申请所属领域普通技术人员的通常理解相同的含义。还应理解的是，除非本申请中有明确的说明，否则在常用词典中定义的词语应被解释为具有与它们在相关技术的上下文中的含义一致的含义，而不应以理想化或过于形式化的意义解释。
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。另外，除非明确限定或与上下文相矛盾，否则本申请所记载的方法中包含的具体步骤不必限于所记载的顺序，而可以任意顺序执行或并行地执行。下面将参考附图并结合实施例来详细说明本申请。
本申请一方面提供了一种用于卫星导航系统B1宽带复合信号的无模糊定位方法。图1示出了根据本申请实施例的用于卫星导航系统B1宽带复合信号的无模糊定位方法的流程图。如图1所示，所述方法1000包括：
S110，通过接收机接收多颗卫星的B1宽带复合信号，其中，B1宽带复合信号包括：B1I信号和B1C信号；
S120，获取所述多颗卫星的所述B1I信号和所述B1C信号的观测量；
S130，基于所述B1I信号和所述B1C信号的观测量，确定所述多颗卫星中各颗卫星到所述接收机的伪距值和码伪距值；
S140，基于所述伪距值和所述码伪距值，确定所述B1宽带复合信号的传播延迟的模糊度整数解；
S150，通过所述模糊度整数解更正所述B1宽带复合信号的传播延迟的模糊度，获取无模糊的B1宽带复合信号的传播延迟，并基于所述无模糊的B1宽带复合信号的传播延迟，对所述接收机所处位置进行无模糊的定位。
下面将进一步对上述各个步骤进一步描述。
在步骤S110中，通过接收机接收多颗卫星的B1宽带复合信号，其中，B1宽带复合信号包括：B1I信号和B1C信号。在一个实施方式中，接收机适用于卫星导航信号，如GPS接收机等。
在步骤是S120中，获取步骤S110接收的B1I信号和B1C信号的观测量。在一个实施方式中，所述B1I信号和所述B1C信号的观测量可以包括：B1C信号的载波相位、B1I信号的码传播延迟、以及B1I信号的副载波传播延迟；通过B1C信号的载波相位和B1I信号的码传播延迟确定码伪距值，以及通过B1I信号的副载波传播延迟和码传播延迟确定伪距值。
图2示出了根据本申请实施例的获取B1C信号的观测量的流程图。对于同一颗卫星来说，BDS3 B1C信号和BDS3 B1I信号拥有相同的载波频率、多普勒频率，两者的载波相位严格相等，也即当确定任一信号的载波信息时，另一信号的载波信息也可以被确定。由于BDS3B1宽带复合信号中B1C信号是实信号，没有SCBOC(14,2)载波相位与副载波相位的耦合情况，所以B1C信号可以通过载波跟踪环路(PLL)获取其准确的载波信息，从而为BDS3 B1I的跟踪提供帮助。
在一个实施方式中，如图2所示，步骤S120可以包括：S1201，基于B1C信号、B1I信号的码传播延迟、以及载波跟踪环路，确定B1C信号的载波频率和相位；S1202，基于B1I信号的码跟踪环路，确定B1I信号的码传播延迟。
图3示出了根据本申请实施例的获取B1I信号的观测量的流程图。在一个实施方式中，如图3所示，步骤S120还可以包括：S1201’，基于B1I信号、B1C信号的载波频率和相位、以及B1I信号的码跟踪环路，确定B1I信号的扩频码频率和相位；S1202’，基于B1I信号的扩频码频率和相位，确定B1I信号的码传播延迟；S1203’，基于B1I信号、B1C信号的载波频率和相位、以及副载波跟踪环路，确定B1I信号的副载波频率和相位；以及S1204’，基于B1I信号的副载波频率和相位，确定B1I信号的副载波传播延迟。
在步骤S130中，基于所述B1I信号和所述B1C信号的观测量，确定各颗所述卫星到所述接收机的伪距值和码伪距值。在一个实施方式中，通过B1C信号的载波相位和B1I信号的码传播延迟确定码伪距值；以及通过B1I信号的副载波传播延迟和B1I信号的码传播延迟确定伪距值。
图4示出了根据本申请实施例的基于所述B1I信号和所述B1C信号的观测量确定卫星到接收机的伪距值和码伪距值的流程图。在一个实施方式中，如图4所示，步骤S130可以包括：S1301，通过合并码传播延迟和副载波传播延迟，确定B1I信号的传播延迟；S1302，基于B1C信号的载波相位，对码传播延迟进行平滑；S1303，基于B1I信号的传播延迟和平滑后的码传播延迟，确定各颗卫星到接收机的伪距值和码伪距值。
在步骤S140中，基于所述伪距值和所述码伪距值，确定所述B1宽带复合信号的传播延迟的模糊度整数解。在一个实施方式中，步骤S140可以包括：基于伪距值和码伪距值，确定B1宽带复合信号的传播延迟的模糊度浮点解；以及基于该模糊度浮点解，确定其模糊度整数解。
在步骤S150中，通过所述模糊度整数解更正所述B1宽带复合信号的传播延迟的模糊度，获取无模糊的B1宽带复合信号的传播延迟，以对所述接收机所处位置进行无模糊的定位。在一个实施方式中，可以通过LAMBDA算法，基于模糊度浮点解确定模糊度整数解。
具体地，将B1宽带复合信号的传播延迟的模糊度浮点解送入LAMBDA算法，利用LAMBDA算法搜索得到所述模糊度整数解，当由码相位估计的副载波模糊度无误时，求解出的传播延迟模糊度为零，而当由码相位估计的副载波模糊度错误时，求解出的传播延迟模糊度为非零整数，据此更正估计错误的传播延迟模糊度。
本申请另一方面提供了一种用于卫星导航系统B1宽带复合信号的无模糊定位装置。图5示出了根据本申请实施例的用于卫星导航系统B1宽带复合信号的无模糊定位装置的框图。如图5所示，该装置500包括：
接收机510，用于接收多颗卫星的B1宽带复合信号，其中，B1宽带复合信号包括：B1I信号和B1C信号；
跟踪模块520，用于获取所述多颗卫星的所述B1I信号和所述B1C信号的观测量；以及
模糊度解算模块530，其配置为：基于所述B1I信号和所述B1C信号的观测量，确定各颗所述卫星到所述接收机的伪距值和码伪距值；基于所述伪距值和所述码伪距值，确定所述B1宽带复合信号的传播延迟的模糊度整数解；以及通过所述模糊度整数解更正所述B1宽带复合信号的传播延迟的模糊度，获取无模糊的B1宽带复合信号的传播延迟，以对所述接收机所处位置进行无模糊的定位。
对于接收机510，用于接收多颗卫星的B1宽带复合信号，其中，B1宽带复合信号包括：B1I信号和B1C信号。在一个实施方式中，接收机适用于卫星导航信号，如GPS接收机等。
对于跟踪模块520，用于获取接收机510接收的B1I信号和B1C信号的观测量。在一个实施方式中，B1I信号和B1C信号的观测量可以包括：B1C信号的载波相位、B1I信号的码传播延迟、以及B1I信号的副载波传播延迟；并通过B1C信号的载波相位和B1I信号的码传播延迟确定码伪距值，以及通过B1I信号的副载波传播延迟和码传播延迟确定伪距值。
在一个实施方式中，跟踪模块520可以配置为：基于B1C信号、B1I信号的码传播延迟、以及载波跟踪环路，确定B1C信号的载波频率和相位；以及基于B1I信号的码跟踪环路，确定B1I信号的码传播延迟。
图6示出了根据本申请实施例的跟踪模块和模糊度解算模块的示意图。其中，跟踪模块520用于跟踪B1C信号的载波相位、B1I信号的码传播延迟、以及B1I信号的副载波传播延迟。
在一个实施方式中，如图6所示，当跟踪第i颗卫星信号时，对于BDS3 B1C信号来说，接收机接收的信号r(t)首先与载波NCO(数控振荡器)521产生的载波相乘，得到剥离掉载波的第i颗卫星的BDS3 B1C的基带信号gB1C,i(t)，其中是载波跟踪环路523估计出的第i颗卫星信号的B1C信号的载波频率和相位。对gB1C,i(t)取实部和虚部得到基带IQ信号IB1C,i,QB1C,i。
基带IQ信号IB1C,i,QB1C,i与B1C码-副载波NCO 522生成的第i颗卫星的本地即时B1C基带信号相乘并积分，确定信号IPB1C,i,QPB1C,i。其中，B1C码-副载波NCO 522的相位是由BDS3 B1I跟踪提供的。
基带IQ信号IPB1C,i,QPB1C,i经过载波跟踪环路523鉴相及滤波后,得到载波频率和相位的误差，以更正载波NCO 523以及B1C信号的载波相位当B1C信号的载波相位确定后，B1I信号的载波相位也被确定，由此只需要估计B1I信号的副载波相位，而无需同时估计其载波相位，有效地解决了调制信号的副载波相位与载波相位耦合的技术问题。
在一个实施方式中，跟踪模块520还可以用于，基于B1I信号、B1C信号的载波频率和相位、以及B1I信号的码跟踪环路，确定B1I信号的扩频码频率和相位；基于B1I信号的扩频码频率和相位，确定B1I信号的码传播延迟；基于B1I信号、B1C信号的载波频率和相位、以及副载波跟踪环路，确定B1I信号的副载波频率和相位；以及基于B1I信号的副载波频率和相位，确定B1I信号的副载波传播延迟。
在一个实施方式中，如图6所示，接收信号r(t)与载波NCO 521产生的载波相乘，得到基带信号gB1I,i(t)，该信号通过与B1I副载波NCO 524生成的B1I信号的副载波相乘剥离掉B1I信号的副载波后，分别与B1I码NCO 525生成的本地B1I提前码、即时码以及延迟码相乘并积分，其中信号IPB1I,i,QPB1I,i进入副载波跟踪环路527鉴相并滤波，以更新B1I副载波NCO624，并确定B1I信号的副载波传播延迟信号IEB1I,i,QEB1I,i,ILB1I,i,QLB1I,i进入码跟踪环路526鉴相并滤波，以更新B1I码NCO 525，确定B1I信号的码传播延迟/>并基于/>和/>获取B1I信号的传播延迟。
由于B1C信号与B1I信号具有相同的传播延迟，因此，跟踪B1I信号得到的B1I信号的传播延迟可以直接传递给B1C信号，使B1C码-副载波NCO 522可以准确地生成B1C信号的本地即时码，缓解了SCBOC(14,2)调制信号严重的模糊跟踪的技术问题。
对于模糊度解算模块530，在一个实施方式中，模糊度解算模块530可以配置为：通过合并码传播延迟和副载波传播延迟，确定B1I信号的传播延迟，合并方式为：其中round()函数返回对/>进行四舍五入的结果，B1I信号的传播延迟模糊度为：/>当码传播延迟/>为副载波传播延迟/>提供正确的副载波传播延迟模糊度时，B1I信号的传播延迟模糊度为0，当码传播延迟/>不能为副载波传播延迟/>提供正确的副载波传播延迟模糊度时，B1I信号的传播延迟模糊度不为0；基于B1C信号的载波相位，对码传播延迟进行平滑；基于B1I信号的传播延迟和平滑后的码传播延迟，确定多颗卫星中各颗卫星到接收机的伪距值和码伪距值。
在一个实施方式中，模糊度解算模块530可以配置为：基于伪距值和码伪距值，确定B1宽带复合信号的传播延迟的模糊度浮点解；以及基于该模糊度浮点解，确定其模糊度整数解。
在一个实施方式中，经由码传播延迟修正的副载波传播延迟模糊度具有整数性，同时参与位置解算的码伪距观测值是低精度的，其定位结果有较大误差，会扩散到副载波模糊度求解值中。如图6所示的模糊度解算模块530可以用于修正B1I信号的传播延迟的模糊度。
如图6所示，模糊度解算模块530通过组合模块531，对B1I信号的码传播延迟和其副载波传播延迟/>进行合并组合，生成B1I信号的传播延迟/>
通过伪距计算模块532，基于码跟踪环路526获取的B1I信号的码传播延迟建立如下无模糊的伪距方程：
其中，为根据/>计算得到的第i颗卫星的码伪距值，Xs,i为第i颗卫星的卫星坐标，Xu为接收机坐标，C为光速，δT为接收机钟差，Ii为第i颗卫星的电离层延迟，Ti为第i颗卫星的对流层延迟，/>为根据B1I信号的码传播延迟计算的第i颗卫星的码伪距噪声。
通过伪距计算模块532，基于B1I信号的传播延迟建立如下的有模糊的伪距方程：
其中，为根据/>计算得到的第i颗卫星的伪距值，/>为第i颗卫星的由B1I信号的传播延迟计算的伪距噪声，/>为一个BDS3B1I信号的副载波码片对应的距离，ΔNi为第i颗卫星的B1宽带复合信号的传播延迟的模糊度。将上述两个伪距方程一阶展开，可以得到：
其中，di为第i颗卫星与接收机的几何距离，ei为从接收机指向第i颗卫星的单位向量，(δx,δy,δz)T为接收机坐标的误差。
根据本申请一实施方式，为了提高副载波模糊度求解的正确性，模糊度解算模块530通过载波平滑模块533，采用载波平滑伪距的方式提高码伪距观测量的精度。载波平滑码伪距的具体方式为：
其中，M代表平滑窗长度，k代表当前时刻，k-1代表上一时刻，为由载波NCO 521估计出的第i颗卫星的载波相位。因此，模糊度解算模块530基于各颗卫星的经过载波平滑后的码伪距值/>和由B1I信号的传播延迟计算的伪距值/>获取B1宽带复合信号的传播延迟的模糊度ΔNi。
根据本申请一实施方式，对具有m颗卫星的系统，得到其伪距方程组如下：
从上式来看，方程待估计参数一共有(m+4)个，为了保证该方程是秩满的，卫星个数m需要大于等于4，与伪距直接定位需要的最少卫星数一致，这表明只要卫星数能达到定位的要求，模糊度解算模块530就能据此解算出B1宽带复合信号的传播延迟的模糊度。其定义如下：
由此，伪距方程组可以整理为：
Y＝Hα+v 8)
定义v的协方差矩阵为Q，求解上式等价于求解：
采用加权最小二乘求解，得到：
α＝(HTQH)-1HTQY 10)
因为接收机位置未知，所以加权最小二乘算法需要经过迭代可收敛，迭代公式为：
βk＝βk-1+δβk 11)
其中，βk＝(x,y,z,δT)T为第k次迭代的接收机坐标及钟差。当||δβ||<Th时，该算法收敛，得到B1宽带复合信号的传播延迟的模糊度浮点解其中，Th为设定阈值，用户可以根据实际情况进行设定，例如Th可以设置为1E-4。
根据本申请的一种实施方式，将该模糊度浮点解及其协方差矩阵输入LAMBDA模块535，解算出模糊度整数解/>
根据LAMBDA算法解算出的B1宽带复合信号的传播延迟的模糊度整数解，通过确定伪距值模块536得到无模糊的伪距值其可以表示为：
需要指出的是，虽然在解算模糊度浮点解的时候会得到接收机的定位结果，但因为没有消除副载波模糊度估计不正确的问题，所以该定位结果是有偏的，只用于辅助求解B1宽带复合信号的传播延迟的模糊度，而不作为最终的接收机定位结果。最终将无模糊的伪距值输入定位模块537，从而实现对接收机的无模糊定位。
基于同一历元中各颗卫星信号的无模糊的B1I信号的码传播延迟以及有模糊的B1I信号的传播延迟/>通过伪距计算模块532建立基于B1I信号的码传播延迟和其传播延迟的伪距方程；并利用B1C信号的载波相位，通过对B1I信号的码传播延迟进行平滑，以提高副载波模糊度求解的正确性；然后将各颗卫星的B1宽带复合信号的传播延迟的模糊度ΔNi作为待估计参数，与接收机位置以及钟差一同求解，得到传播延迟模糊度浮点解；通过LAMBDA算法搜索得到传播延迟模糊度整数解，当由B1I信号的码传播延迟估计的副载波模糊度无误时，求解出的B1宽带复合信号的传播延迟的模糊度为零，而当由B1I信号的码传播延迟估计的B1I信号的副载波的模糊度错误时，求解出的B1宽带复合信号的传播延迟模糊度为非零整数，据此更正估计错误的传播延迟模糊度，对接收机进行无模糊的定位。
本申请另一方面提供了一种用于卫星导航系统B1宽带复合信号的无模糊定位装置，所述装置包括通过总线或互连件连接的处理器、存储器。处理器可以表示其中包括单个处理器内核或多个处理器内核的单个处理器或多个处理器。处理器可以表示一个或多个通用处理器，诸如，微处理器、中央处理单元(CPU)等。更具体地，处理器可以是复杂指令集计算(CISC)微处理器、精简指令集计算(RISC)微处理器、超长指令字(VLIW)微处理器、或实施其它指令集的处理器、或实施指令集组合的处理器。处理器还可以是一个或多个专用处理器，诸如，专用集成电路(ASIC)、蜂窝或基带处理器、现场可编程门阵列(FPGA)、数字信号处理器(DSP)、网络处理器、图形处理器、通信处理器、加密处理器、协处理器、嵌入式处理器、或者能够处理指令的任何其它类型的逻辑。
处理器(其可以是低功率多核处理器套接口，诸如超低电压处理器)可以充当用于与所述系统的各种部件通信的主处理单元和中央集线器。这种处理器可以实施为片上系统(SoC)。处理器被配置成执行用于执行本申请所讨论的方法的指令。
处理器可以与存储器通信，存储器在一个实施方式中可以经由多个存储器装置实施以提供给定量的系统存储。存储器可以包括一个或多个易失性存储(或存储器)装置，诸如，随机存取存储器(RAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、静态RAM(SRAM)或者其它类型的存储装置。存储器可以存储包括由处理器或任何其它装置执行的指令序列的信息。例如，各种操作系统、装置驱动程序、固件(例如，输入输出基本系统或BIOS)和/或应用的可执行代码和/或数据可以加载到存储器中并由处理器执行。操作系统可以是任何类型的操作系统，例如，机器人操作系统(ROS)、来自公司的/>操作系统、来自苹果公司的Mac/>来自/>公司的/>LINUX、UNIX，或者其它实时或嵌入式操作系统。
为了提供对诸如数据、应用、一个或多个操作系统等信息的永久性存储，大容量存储设备也可以联接到处理器。在各种实施方式中，为了实现更薄且更轻的系统设计并且改进系统响应性，这种大容量存储设备可以经由固态装置(SSD)来实施。然而，在其它实施方式中，大容量存储设备可以主要使用硬盘驱动器(HDD)来实施，其中较小量的SSD存储设备充当SSD高速缓存以在断电事件期间实现上下文状态以及其它此类信息的非易失性存储，从而使得在系统活动重新启动时能够实现快速通电。另外，闪存装置可以例如经由串行外围接口(SPI)联接到处理器。这种闪存装置可以提供系统软件的非易失性存储，所述系统软件包括所述系统的BIOS以及其它固件。
本申请另一方面提供了一种计算机可读介质，计算机可读存储介质也可以用来永久性地存储以上描述的方法。虽然，计算机可读存储介质在示例性实施方式中被示为单个介质，但是术语“计算机可读存储介质”应当被认为包括存储所述一个或多个指令集的单个介质或多个介质(例如，集中式或分布式数据库和/或相关联的高速缓存和服务器)。术语“计算机可读存储介质”还应当被认为包括能够存储或编码指令集的任何介质，所述指令集用于由机器执行并且使得所述机器执行本公开的任何一种或多种方法。因此，术语“计算机可读存储介质”应当被认为包括但不限于固态存储器以及光学介质和磁性介质，或者任何其它非暂时性机器可读介质。
以上描述仅为本申请的较佳实施方式以及对所运用技术原理的说明。本领域技术人员应当理解，本申请中所涉及的保护范围，并不限于上述技术特征的特定组合而成的技术方案，同时也应涵盖在不脱离技术构思的情况下，由上述技术特征或其等同特征进行任意组合而形成的其它技术方案。例如上述特征与本申请中公开的(但不限于)具有类似功能的技术特征进行互相替换而形成的技术方案。

Claims (14)

1.一种用于卫星导航系统B1宽带复合信号的无模糊定位方法，包括：
接收多颗卫星的B1宽带复合信号，其中，所述B1宽带复合信号包括B1I信号以及与所述B1I信号具有相同载波相位的B1C信号；
获取所述B1I信号和所述B1C信号的观测量，其中所述观测量包括所述B1I信号的码传播延迟、所述B1I信号的副载波传播延迟以及所述B1C信号的载波相位；
通过所述B1I信号的码传播延迟和所述B1I信号的副载波传播延迟确定伪距值；
通过所述B1I信号的码传播延迟和所述B1C信号的载波相位确定码伪距值；
基于所述伪距值和所述码伪距值，确定所述B1宽带复合信号的传播延迟的模糊度整数解；以及
通过所述模糊度整数解更正所述B1宽带复合信号的传播延迟的模糊度，获取无模糊的B1宽带复合信号的传播延迟，以对接收机所处位置进行无模糊的定位。
2. 如权利要求1所述的方法，其中，基于所述伪距值和所述码伪距值，确定所述B1宽带复合信号的传播延迟的模糊度整数解的步骤包括：
基于所述伪距值和所述码伪距值，确定所述B1宽带复合信号的传播延迟的模糊度浮点解；以及
基于所述模糊度浮点解确定所述模糊度整数解。
3.如权利要求2所述的方法，其中，获取所述B1I信号和所述B1C信号的观测量的步骤包括：
基于所述B1C信号、所述B1I信号的码传播延迟、以及载波跟踪环路，确定所述B1C信号的载波频率和所述B1C信号的载波相位；以及
基于所述B1I信号的码跟踪环路，确定所述B1I信号的码传播延迟。
4.如权利要求3所述的方法，其中，获取所述B1I信号和所述B1C信号的观测量的步骤进一步包括：
基于所述B1I信号、所述B1C信号的载波频率和相位、以及所述B1I信号的码跟踪环路，确定所述B1I信号的扩频码频率和相位；
基于所述B1I信号的扩频码频率和相位，确定所述B1I信号的码传播延迟;
基于所述B1I信号、所述B1C信号的载波频率和相位、以及副载波跟踪环路，确定所述B1I信号的副载波频率和相位；以及
基于所述B1I信号的副载波频率和相位，确定所述B1I信号的副载波传播延迟。
5.如权利要求4所述的方法，其中，基于所述观测量，确定各颗所述卫星到所述接收机的伪距值和码伪距值的步骤包括：
通过合并所述B1I信号的码传播延迟和所述B1I信号的副载波传播延迟，确定所述B1I信号的传播延迟；
基于所述B1C信号的载波相位，对所述B1I信号的码传播延迟进行平滑；以及
基于所述B1I信号的传播延迟和平滑后的所述B1I信号的码传播延迟，确定各颗所述卫星到所述接收机的所述伪距值和所述码伪距值。
6.如权利要求2所述的方法，其中，基于所述模糊度浮点解确定所述模糊度整数解的步骤包括：
通过LAMBDA算法，基于所述模糊度浮点解确定所述模糊度整数解。
7.一种用于卫星导航系统B1宽带复合信号的无模糊定位装置，包括：
接收机，用于接收多颗卫星的B1宽带复合信号，其中，所述B1宽带复合信号包括B1I信号以及与所述B1I信号具有相同载波相位的B1C信号；
跟踪模块，用于获取所述B1I信号和所述B1C信号的观测量，其中所述观测量包括所述B1I信号的码传播延迟、所述B1I信号的副载波传播延迟以及所述B1C信号的载波相位；
模糊度解算模块，其配置为：
通过所述B1I信号的码传播延迟和所述B1I信号的副载波传播延迟确定伪距值；
通过所述B1I信号的码传播延迟和所述B1C信号的载波相位确定码伪距值；
基于所述伪距值和所述码伪距值，确定所述B1宽带复合信号的传播延迟的模糊度整数解；以及
通过所述模糊度整数解更正所述B1宽带复合信号的传播延迟的模糊度，获取无模糊的B1宽带复合信号的传播延迟，以对所述接收机所处位置进行无模糊的定位。
8. 如权利要求7所述的装置，其中，所述模糊度解算模块被配置为：
基于所述伪距值和所述码伪距值，确定所述B1宽带复合信号的传播延迟的模糊度浮点解；以及
基于所述模糊度浮点解确定所述模糊度整数解。
9.如权利要求8所述的装置，其中，
所述跟踪模块被配置为：
基于所述B1C信号、所述B1I信号的码传播延迟、以及载波跟踪环路，确定所述B1C信号的载波频率和所述B1C信号的载波相位；以及
基于所述B1I信号的码跟踪环路，确定所述B1I信号的码传播延迟。
10.如权利要求9所述的装置，其中，所述跟踪模块被进一步配置为：
基于所述B1I信号、所述B1C信号的载波频率和相位、以及所述B1I信号的码跟踪环路，确定所述B1I信号的扩频码频率和相位；
基于所述B1I信号的扩频码频率和相位，确定所述B1I信号的码传播延迟；
基于所述B1I信号、所述B1C信号的载波频率和相位、以及副载波跟踪环路，确定所述B1I信号的副载波频率和相位；以及
基于所述B1I信号的副载波频率和相位，确定所述B1I信号的副载波传播延迟。
11.如权利要求10所述的装置，其中，所述模糊度解算模块被配置为：
通过合并所述B1I信号的码传播延迟和所述B1I信号的副载波传播延迟，确定所述B1I信号的传播延迟；
基于所述B1C信号的载波相位，对所述B1I信号的码传播延迟进行平滑；
基于所述B1I信号的传播延迟和平滑后的所述B1I信号的码传播延迟，确定各颗所述卫星到所述接收机的所述伪距值和所述码伪距值。
12.如权利要求8所述的装置，其中，所述模糊度解算模块被配置为：
通过LAMBDA算法，基于所述模糊度浮点解确定所述模糊度的所述整数解。
13.一种用于卫星导航系统B1宽带复合信号的无模糊定位装置，包括：
存储器，存储有计算机可执行的指令；
处理器，执行所述指令以执行如权利要求1-6中任一项所述的方法。
14.一种存储介质，存储有计算机可执行的指令，所述指令被一个或多个处理器执行时实现如权利要求1-6中任一项所述的方法。
CN202210156560.0A 2022-02-21 2022-02-21 用于卫星导航系统b1宽带复合信号的无模糊定位方法及其装置 Active CN114527499B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202210156560.0A CN114527499B (zh) 2022-02-21 2022-02-21 用于卫星导航系统b1宽带复合信号的无模糊定位方法及其装置
PCT/CN2023/077393 WO2023155920A1 (zh) 2022-02-21 2023-02-21 用于卫星导航系统b1宽带复合信号的无模糊定位方法及其装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210156560.0A CN114527499B (zh) 2022-02-21 2022-02-21 用于卫星导航系统b1宽带复合信号的无模糊定位方法及其装置

Publications (2)

Publication Number Publication Date
CN114527499A CN114527499A (zh) 2022-05-24
CN114527499B true CN114527499B (zh) 2023-08-01

Family

ID=81624711

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210156560.0A Active CN114527499B (zh) 2022-02-21 2022-02-21 用于卫星导航系统b1宽带复合信号的无模糊定位方法及其装置

Country Status (2)

Country Link
CN (1) CN114527499B (zh)
WO (1) WO2023155920A1 (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114527499B (zh) * 2022-02-21 2023-08-01 清华大学 用于卫星导航系统b1宽带复合信号的无模糊定位方法及其装置
CN114675310B (zh) * 2022-05-30 2022-09-30 长沙金维信息技术有限公司 载波半周修复方法及其rtk整周模糊度固定方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2413160B1 (en) * 2010-07-28 2016-11-23 Airbus DS GmbH Method for resolving the binary offset carrier (BOC) signals tracking ambiguity in a double estimator receiver
EP3109673B1 (en) * 2015-06-24 2020-01-08 Centre National d'Etudes Spatiales Gnss receiver with an improved capability to resolve sub-carrier tracking ambiguities
EP3112905B1 (en) * 2015-07-03 2018-02-21 Centre National d'Etudes Spatiales Gnss receiver calculating a non-ambiguous discriminator to resolve subcarrier tracking ambiguities
EP3141931B1 (en) * 2015-09-14 2019-12-04 Airbus Defence and Space GmbH Tracking of signals with at least one subcarrier
CN106291637B (zh) * 2016-08-05 2018-12-11 清华大学 基于全伪距和部分伪距的定位方法
CN108008422B (zh) * 2016-11-02 2020-06-02 清华大学 伪卫星跳时信号捕获装置和方法
CN112578422B (zh) * 2019-09-27 2024-09-17 清华大学 复合导航信号接收方法及接收机
CN113703021B (zh) * 2021-07-29 2023-09-29 西安空间无线电技术研究所 一种基于码伪距的秒级实时高精度定位方法与系统
CN114063124B (zh) * 2021-11-08 2022-09-06 北京航空航天大学合肥创新研究院(北京航空航天大学合肥研究生院) 一种北斗b1c信号快速捕获方法及存储介质
CN114527499B (zh) * 2022-02-21 2023-08-01 清华大学 用于卫星导航系统b1宽带复合信号的无模糊定位方法及其装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于双环估计的CBOC调制信号无模糊跟踪方法;曲博等;华中科技大学学报（自然科学版）;第58-62页 *

Also Published As

Publication number Publication date
CN114527499A (zh) 2022-05-24
WO2023155920A1 (zh) 2023-08-24

Similar Documents

Publication Publication Date Title
CN114527499B (zh) 用于卫星导航系统b1宽带复合信号的无模糊定位方法及其装置
US10705223B2 (en) Low power asynchronous GPS baseband processor
Wendel et al. A robust technique for unambiguous BOC tracking
JPH0659013A (ja) Ｇｐｓ受信機の信号捕捉方法
Diessongo et al. Precise position determination using a Galileo E5 single-frequency receiver
US8391339B2 (en) Correlator sum method for spread spectrum signal receivers
Tamazin et al. Robust fine acquisition algorithm for GPS receiver with limited resources
Tabatabaei et al. Robust adaptive joint tracking of GNSS signal code phases in urban canyons
CN109991630B (zh) 用于组合信号的系统和方法
CN109564293B (zh) 用于全球导航卫星系统(gnss)信号跟踪的装置和方法
Sixin et al. Doppler frequency-code phase division multiple access technique for LEO navigation signals
Li et al. Rapid acquisition assisted by navigation data for inter-satellite links of navigation constellation
Yao et al. Spreading code phase measurement technique for snapshot GNSS receiver
CN105204050B (zh) 一种惯性辅助的多通道混合型矢量跟踪方法
Jovanovic et al. Two‐Step Galileo E1 CBOC Tracking Algorithm: When Reliability and Robustness Are Keys!
Yoo et al. A coarse-time positioning method for improved availability
Nie et al. A fast method for the acquisition of weak long-code signal
CN115774275A (zh) 导航卫星信号的双频多分量联合跟踪方法、设备和介质
Deng et al. An unambiguous tracking technique for sine-BOC (kn, n) modulated GNSS signals
US20210255333A1 (en) Ultralow power gnss receiver
CN105699997B (zh) 一种使用glonass单频信号进行差分定位的方法
Hu et al. An efficient method for GPS multipath mitigation using the Teager-Kaiser-operator-based MEDLL
US11747487B2 (en) GNSS receiver clock frequency drift detection
Chu et al. Mixed QZSS/GPS/Galileo RTK positioning with a trade-off between ionosphere-weighting and ionosphere-fixed models: An ambiguity resolution method for Japan
Kou et al. Processing GPS L2C signals under ionospheric scintillations

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant