CN112083452B - 卫星导航接收机的环路跟踪系统及方法 - Google Patents
卫星导航接收机的环路跟踪系统及方法 Download PDFInfo
- Publication number
- CN112083452B CN112083452B CN202010888948.0A CN202010888948A CN112083452B CN 112083452 B CN112083452 B CN 112083452B CN 202010888948 A CN202010888948 A CN 202010888948A CN 112083452 B CN112083452 B CN 112083452B
- Authority
- CN
- China
- Prior art keywords
- spacecraft
- signal
- loop
- coherent integration
- code
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000010354 integration Effects 0.000 claims abstract description 116
- 230000001427 coherent effect Effects 0.000 claims abstract description 114
- 238000006243 chemical reaction Methods 0.000 claims abstract description 29
- 238000001914 filtration Methods 0.000 claims abstract description 21
- 238000005070 sampling Methods 0.000 claims abstract description 13
- 230000001133 acceleration Effects 0.000 claims description 48
- 230000036461 convulsion Effects 0.000 claims description 48
- 238000012937 correction Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 8
- 230000008569 process Effects 0.000 description 7
- 230000035945 sensitivity Effects 0.000 description 7
- 230000003044 adaptive effect Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/29—Acquisition or tracking or demodulation of signals transmitted by the system carrier including Doppler, related
-
- 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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/30—Acquisition or tracking or demodulation of signals transmitted by the system code related
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
本发明涉及导航技术领域,提供一种卫星导航接收机的环路跟踪系统及方法,包括:信号采样模块,获取导航卫星的下变频信号,并将下变频信号进行模数转换输出数字中频信号;相干积分模块,根据本地载波、本地伪码对数字中频信号的I路和Q路均进行相干积分;码环鉴别模块,根据相干积分的信号进行滤波鉴相并修正码数字振荡器的频率字;载波鉴别模块,获取航天器的动态参数和信噪比,根据动态参数、信噪比和相干积分的信号进行载波鉴相,并修正载波数字振荡器的频率字,使本地载波实时跟踪导航卫星。本发明能够接收到弱卫星信号,兼顾航天器的高动态,解决了强弱卫星信号共存造成的远近效应问题,在近远地轨道均能持续的卫星跟踪。
Description
技术领域
本发明涉及导航技术领域,尤其涉及一种卫星导航接收机的环路跟踪系统及方法。
背景技术
随着高轨卫星需求的逐年增加以及深空探测的发展,我国对导航定位技术提出了更高的要求,由于卫星导航具有定位精度高、无累积误差等优点,成为航天器导航的有效手段。但对于“嫦娥”等高轨飞行器,飞行轨道距离导航卫星远,卫星信号弱且容易受到地球的遮挡,同时高轨飞行器在飞行过程的速度快,因此为实现GNSS(Global NavigationSatellite System,全球卫星导航系统)定位以及对导航卫星的漏导航信号及旁瓣信号的接受,需要提高接收机的灵敏度。
但在航天器高动态条件下,接收到的卫星载波频率和相位有大幅度的变动,另外在高轨道条件下,不同导航卫星到接收机的距离差别大,造成强弱卫星信号共存,将会引起高轨接收机的远近效应问题。
发明内容
本发明实施例提供了一种卫星导航接收机的环路跟踪系统及方法。
本发明实施例的第一方面提供了一种卫星导航接收机的环路跟踪系统,包括:信号采样模块、相干积分模块、码环鉴别模块和载波鉴别模块;
所述信号采样模块用于获取导航卫星的下变频信号,并将所述下变频信号进行模数转换输出数字中频信号;
所述相干积分模块用于获取本地载波和本地伪码,并根据本地载波和本地伪码对所述数字中频信号的I/Q路均进行相干积分,相干积分后的信号输出给所述码环鉴别模块和所述载波鉴别模块;
所述码环鉴别模块用于对所述相干积分后的信号进行滤波鉴相得到码相位误差信号,所述码相位误差信号用于修正码数字振荡器的频率字;
所述载波鉴别模块用于获取航天器的动态参数和所述下变频信号的信噪比,根据所述动态参数、所述信噪比和所述相干积分后的信号进行载波鉴相得到载波相位误差信号,所述载波相位误差信号用于修正载波数字振荡器的频率字,使所述本地载波实时跟踪所述导航卫星。
进一步地,所述码环鉴别模块包括:码环鉴别器和码环滤波器;
所述码环鉴别器用于对所述相干积分后的信号进行码环鉴相;
所述码环滤波器用于对码环鉴相后的信号进行滤波得到码相位误差信号,所述码相位误差信号用于修正码数字振荡器的频率字。
进一步地,在修正码数字振荡器的频率字之后,所述相干积分模块还用于:
获取频率字修正后的所述码数字振荡器的第二伪码,根据所述本地载波和所述第二伪码对所述数字中频信号的Q路进行相干积分,相干积分后的信号输出给所述码环鉴别模块和所述载波鉴别模块。
进一步地,所述载波鉴别模块包括:载波环鉴别器、自适应控制模块、二阶环路滤波器和三阶环路滤波器;
所述载波环鉴别器用于对所述相干积分后的信号进行载波鉴相;
所述自适应控制模块用于获取航天器的动态参数、信噪比和所述载波环鉴别器的初始相干积分时间,根据所述动态参数、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽;
所述二阶环路滤波器用于在所述环路阶数为二阶时,根据所述环路带宽对载波鉴相后的信号进行滤波得到载波相位误差信号,所述载波相位误差信号用于修正数字振荡器的频率字;
所述三阶环路滤波器用于在所述环路阶数为三阶时,根据所述环路带宽对载波鉴相后的信号进行滤波得到载波相位误差信号,所述载波相位误差信号用于修正数字振荡器的频率字。
进一步地,所述航天器的动态参数包括:所述航天器的高度、所述航天器的加速度和所述航天器的加加速度;
所述自适应控制模块具体用于:
在所述航天器的高度低于第一高度阈值时,所述相干积分时间确定为第一时间,环路带宽确定为第一频率;
在所述航天器的高度大于第二高度阈值时,所述相干积分时间确定为第二时间,环路带宽确定为第二频率;
在所述航天器的高度在所述第一高度阈值和所述第二高度阈值之间时,根据所述航天器的加速度、所述航天器的加加速度、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽;
其中,所述第一时间小于所述第二时间,所述第一频率大于所述第二频率。
进一步地,所述自适应控制模块具体用于:
在所述航天器的高度在所述第一阈值和所述第二阈值之间时,且所述航天器的加速度小于第一加速度阈值,以及所述航天器的加加速度小于第一加加速度阈值时,所述环路阶数确定为二阶环路,并通过
确定二阶环路带宽BL2;其中,C/N0为所述信噪比,T为所述初始相干积分时间,
a为所述航天器的加速度;
在所述航天器的高度在所述第一阈值和所述第二阈值之间时,且所述航天器的加速度小于第一加速度阈值,以及所述航天器的加加速度大于第一加加速度阈值时,所述环路阶数确定为三阶环路,并通过
确定三阶环路带宽BL3;其中,C/N0为所述信噪比,T为所述初始相干积分时间,
b为所述航天器的加加速度。
本发明实施例的第二方面提供了一种卫星导航接收机的环路跟踪方法,包括:
通过信号采样模块获取导航卫星的下变频信号,将所述下变频信号进行模数转换输出数字中频信号;
获取本地载波和本地伪码,并根据本地载波和本地伪码对所述数字中频信号的I/Q路均进行相干积分;
对所述相干积分后的信号进行滤波鉴相得到码相位误差信号,所述码相位误差信号用于修正码数字振荡器的频率字;
获取航天器的动态参数和所述下变频信号的信噪比,根据所述动态参数、所述信噪比和所述相干积分后的信号进行载波鉴相得到载波相位误差信号,所述载波相位误差信号用于修正载波数字振荡器的频率字,使所述本地载波实时跟踪所述导航卫星。
进一步地,所述获取航天器的动态参数和所述下变频信号的信噪比,根据所述动态参数、所述信噪比和所述相干积分后的信号进行载波鉴相得到载波相位误差信号,所述载波相位误差信号用于修正载波数字振荡器的频率字,使所述本地载波实时跟踪所述导航卫星,包括:
获取航天器的动态参数、所述下变频信号的信噪比和所述载波环鉴别器的初始相干积分时间;
根据所述动态参数、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽;
根据所述环路阶数和所述环路带宽对载波鉴相后的信号进行滤波得到载波相位误差信号,所述载波相位误差信号用于修正数字振荡器的频率字,使所述本地载波实时跟踪所述导航卫星。
进一步地,所述航天器的动态参数包括:所述航天器的高度、所述航天器的加速度和所述航天器的加加速度;
所述根据所述动态参数、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽,包括:
在所述航天器的高度低于第一高度阈值时,所述相干积分时间确定为第一时间,环路带宽确定为第一频率;
在所述航天器的高度大于第二高度阈值时,所述相干积分时间确定为第二时间,环路带宽确定为第二频率;
在所述航天器的高度在所述第一高度阈值和所述第二高度阈值之间时,根据所述航天器的加速度、所述航天器的加加速度、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽;
其中,所述第一时间小于所述第二时间,所述第一频率大于所述第二频率。
进一步地,所述在所述航天器的高度在所述第一高度阈值和所述第二高度阈值之间时,根据所述航天器的加速度、所述航天器的加加速度、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽,包括:
在所述航天器的高度在所述第一阈值和所述第二阈值之间时,且所述航天器的加速度小于第一加速度阈值,以及所述航天器的加加速度小于第一加加速度阈值时,所述环路阶数确定为二阶环路,并通过
确定二阶环路带宽BL2;其中,C/N0为所述信噪比,T为所述初始相干积分时间,a为所述航天器的加速度;
在所述航天器的高度在所述第一阈值和所述第二阈值之间时,且所述航天器的加速度小于第一加速度阈值,以及所述航天器的加加速度大于第一加加速度阈值时,所述环路阶数确定为三阶环路,并通过
确定三阶环路带宽BL3;其中,C/N0为所述信噪比,T为所述初始相干积分时间,b为所述航天器的加加速度。
本发明实施例的卫星导航接收机的环路跟踪系统及方法与现有技术相比存在的有益效果是:
通过信号采样模块获取导航卫星的下变频信号,并将下变频信号进行模数转换输出数字中频信号,然后相干积分模块根据本地载波、本地伪码对数字中频信号的I路和Q路均进行相干积分,码环鉴别模块根据相干积分的信号进行滤波和鉴相并修正码数字振荡器的频率字,使得本系统能够接收到弱卫星信号,可避免强弱卫星信号共存造成的远近效应问题;最后载波鉴别模块获取航天器的动态参数和信噪比,根据动态参数、信噪比和相干积分的信号进行载波鉴相得到载波相位误差信号,所述载波相位误差信号用于修正载波数字振荡器的频率字,使本地载波实时跟踪导航卫星,兼顾了航天器的高动态,在近远地轨道均能持续的卫星跟踪。
附图说明
图1是本发明实施例提供的卫星导航接收机的环路跟踪系统的结构示意图;
图2是本发明实施例提供的另一种卫星导航接收机的环路跟踪系统的结构示意图;
图3是本发明实施例提供的卫星定位的示意图;
图4是本发明实施例提供的一种卫星导航接收机的环路跟踪方法的实现流程示意图;
图5是图4中的步骤S404的具体实现流程示意图。
具体实施方式
以下描述中,为了说明而不是为了限定,提出了诸如特定系统结构、技术之类的具体细节,以便透彻理解本发明实施例。然而,本领域的技术人员应当清楚,在没有这些具体细节的其它实施例中也可以实现本发明。
为了说明本发明所述的技术方案,下面通过具体实施例来进行说明。
参见图1,为本实施例提供的卫星导航接收机的环路跟踪系统的结构示意图,一种飞行高度在20000km以上的航天器用卫星导航接收机的环路跟踪系统,可以应用于“嫦娥”飞行器或空天飞机等深空探测器的卫星导航定位,可实现航天器在地球轨道及地月转移轨道上的高精度定位,适用于高轨航天器的北斗或GPS(Global Positioning System,全球定位系统)卫星导航接收机终端设备等。为了便于说明,仅示出了与本实施例相关的部分。
所述卫星导航接收机的环路跟踪系统主要包括:信号采样模块100、相干积分模块200、码环鉴别模块300和载波鉴别模块400。
信号采样模块100用于获取导航卫星的下变频信号,并将下变频信号进行模数转换输出数字中频信号。相干积分模块200用于获取本地载波和本地伪码,并根据本地载波和本地伪码对数字中频信号的I/Q路均进行相干积分,相干积分后的信号输出给码环鉴别模块300和载波鉴别模块400。码环鉴别模块300用于对相干积分后的信号进行滤波鉴相得到码相位误差信号,所述码相位误差信号用于修正码数字振荡器的频率字。载波鉴别模块400用于获取航天器的动态参数和下变频信号的信噪比,根据动态参数、信噪比和相干积分后的信号进行载波鉴相得到载波相位误差信号,所述载波相位误差信号用于修正载波数字振荡器的频率字,使本地载波实时跟踪导航卫星。
为实现卫星导航接收机的高灵敏度,需要采用长的相干积分时间并采用较小的环路带宽。但在航天器高动态条件下,接收到的卫星载波频率和相位有大幅度的变动,需要有较大的环路带宽能容忍高动态引起的载波频率和相位的波动。这两者是卫星导航接收机在跟踪环路设计时一个需要协调的矛盾。另外,在高轨道条件下不同导航卫星到接收机的距离差别大,造成强弱卫星信号共存,高轨高灵敏度接收机将引起远近效应。为调和卫星导航接收机的高灵敏度与高动态之间的矛盾,并解决强弱卫星信号并存时的远近效应问题,本发明对卫星跟踪环路带宽进行优化设计。
本实施例的高轨航天器的卫星导航接收机的环路跟踪系统,高轨航天器飞行过程动态大并且导航卫星信号弱,设计了载波鉴别模块400,即卫星跟踪环路,根据航天器的动态参数及卫星信号的信噪比调整载波环阶数及环路带宽参数,既保证能够接收到弱卫星信号又兼顾了航天器的高动态,并能够抑制强弱卫星信号共存造成的远近效应,既能保证卫星导航接收机的灵敏度又能适应大动态环境,适用于北斗、GPS等多种卫星导航系统的接收机终端,具有较强的实用性。
本实施例对信号采样模块100和相干积分模块200的具体结构不进行限定。参见如2,相干积分模块200可以包括多个调制器、两个积分器以及两个清零器等,积分器和清零器用于积分清零,还可以包括sin映射单元和cos映射单元等等。
进一步地,参见图2,本实施例的码环鉴别模块300可以包括:码环鉴别器和码环滤波器。码环鉴别器的输入端与相干积分模块200的输出端连接,码环鉴别器的输出端与码环滤波器的输入端连接,码环滤波器的输出端与码NCO(Numerically ControlledOscillator,数字振荡器)的输入端连接。
码环鉴别器用于对所述相干积分后的信号进行码环鉴相;码环滤波器用于对码环鉴相后的信号进行滤波得到码相位误差信号,所述码相位误差信号用于修正码数字振荡器的频率字。
可选的,在修正码数字振荡器的频率字之后,相干积分模块200还可以用于:
获取频率字修正后的所述码数字振荡器的第二伪码,根据所述本地载波和所述第二伪码对所述数字中频信号的Q路进行相干积分,相干积分后的信号输出给所述码环鉴别模块300和所述载波鉴别模块400,即对修正后的相干积分后的信号再进行码环鉴相和载波鉴相,使得本实施例的环路跟踪系统更加精准。
本实施例针对高轨航天器的飞行高度跨度大,近地过程卫星信号强,远离地球时信号异常弱,并且存在大动态的特殊应用环境问题,本实施例采用自适应调节卫星跟踪环路的积分时间、环路阶数与带宽的方法,实现了高轨航天器飞行全程稳定的卫星导航定位。
具体的,参见图2,本实施例的载波鉴别模块400包括:载波环鉴别器、自适应控制模块、二阶环路滤波器和三阶环路滤波器;载波环鉴别器的输入端与相干积分模块200的输出端连接,载波环鉴别器的输出端与自适应控制模块的输入端连接,自适应控制模块的输出端与二阶环路滤波器的输入端和三阶环路滤波器的输入端均连接,二阶环路滤波器的输出端和三阶环路滤波器的输出端均与载波NCO的输入端连接。
载波环鉴别器用于对所述相干积分后的信号进行载波鉴相;自适应控制模块用于获取航天器的动态参数、信噪比和所述载波环鉴别器的初始相干积分时间,根据所述动态参数、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽;二阶环路滤波器用于在所述环路阶数为二阶时,根据所述环路带宽对载波鉴相后的信号进行滤波得到载波相位误差信号,载波相位误差信号用于修正数字振荡器的频率字;三阶环路滤波器用于在所述环路阶数为三阶时,根据所述环路带宽对载波鉴相后的信号进行滤波得到载波相位误差信号,载波相位误差信号用于修正数字振荡器的频率字。
本实施例的卫星导航信号采用同相-正交环路实现对伪码和载波的跟踪。导航卫星的下变频信号经模数采样变为数字中频信号,数字中频信号、本地载波和本地伪码的I、Q两路相乘后经过相干积分,分别送入码环(码环鉴别模块300)以及载波环(载波鉴别模块400)进行相位鉴别。
码环部分,采用延迟锁定环路实现鉴相与跟踪。具体的,数字中频信号和本地载波I、Q两路相乘,然后经积分清除后送入码环鉴别器得到码相位误差信号,码相位误差信号经过码环滤波器后用于修正码NCO的频率字。其中,滤波后的码相位误差信号,得到的是所述本地伪码和卫星码的相位偏差,码NCO通过数字形式产生本地伪码的频率字。可选的,用码相位误差信号修正码NCO的频率字的方法,可以是直接加减操作,例如,已知GPS卫星信号共1023个码相位,当前码NCO的频率字的码相位为X,当得到的本地伪码和卫星码的相位偏差为1时,就修正码NCO的频率字为X+1。
载波环部分,相干积分后的两路信号即时累加,累加后的结果送入载波环鉴别器进行载波鉴相,自适应控制模块根据获取的高轨航天器的动态参数以及输入信号的信噪比,选择合适的载波环路阶数与环路滤波器参数(环路带宽),对载波鉴相后的信号进行滤波得到载波相位误差信,用于修正载波NCO的频率字,使本地载波实时跟踪卫星载波频率的动态变化。
可选的,本实施例的航天器的动态参数可以包括但不限定于以下参数:所述航天器的高度、所述航天器的速度、所述航天器的加速度和所述航天器的加加速度。
进一步地,自适应控制模块具体用于:
在所述航天器的高度低于第一高度阈值时,所述相干积分时间确定为第一时间,环路带宽确定为第一频率。
在所述航天器的高度大于第二高度阈值时,所述相干积分时间确定为第二时间,环路带宽确定为第二频率。
在所述航天器的高度在所述第一高度阈值和所述第二高度阈值之间时,根据所述航天器的加速度、所述航天器的加加速度、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽。其中,所述第一时间小于所述第二时间,所述第一频率大于所述第二频率。
示例性的,如图3,在航天器的飞行高度低于300km时,卫星信号强度强,相干积分时间采用1ms,环路带宽采用10Hz;航天器的飞行高度大于3000km时,卫星信号较弱,相干积分时间采用20ms,环路带宽采用1Hz;当飞行高度在300km~3000km之间时,通过自适应控制模块自动调节环路带宽,即通过信噪比、加速度和加加速度值确定环路阶数与环路带宽保证信号的稳定跟踪。
进一步地,本实施例的自适应控制模块具体用于:
在所述航天器的高度在所述第一阈值和所述第二阈值之间时,且所述航天器的加速度小于第一加速度阈值,以及所述航天器的加加速度小于第一加加速度阈值时,所述环路阶数确定为二阶环路,并通过
确定二阶环路带宽BL2;其中,C/N0为所述信噪比,T为所述初始相干积分时间,a为所述航天器的加速度。
示例性的,当航天器的飞行高度在300km~3000km之间时,且航天器的加速度小于4g,以及航天器的加加速度小于0.1g/s时采用二阶环路,并根据上述公式确定二阶环路滤波器的环路带宽。
在所述航天器的高度在所述第一阈值和所述第二阈值之间时,且所述航天器的加速度小于第一加速度阈值,以及所述航天器的加加速度大于第一加加速度阈值时,所述环路阶数确定为三阶环路,并通过
确定三阶环路带宽BL3;其中,C/N0为所述信噪比,T为所述初始相干积分时间,b为所述航天器的加加速度。
示例性的,当航天器的飞行高度在300km~3000km之间时,且航天器的加速度小于4g,以及航天器的加加速度大于0.1g/s时采用三阶环路,并根据上述公式确定三阶环路滤波器的环路带宽。
本实施例的自适应控制以稳定跟踪和环路测量误差最小为目标。在卫星导航接收机跟踪环路中,码环的鲁棒性和动态适应能力较强,所以本实施例对C/A码的跟踪仅对积分时间进行优化,保证跟踪灵敏度。而对于载波跟踪环路,则以减少环路噪声和保证环路稳定性为出发点,将载波环路带宽和阶数作为自适应控制的两个重要参数。
具体的,自适应控制模块的输入变量为信噪比、航天器的动态参数和初始的相干积分时间,最优环路带宽的选择以环路测量误差最小为目标,环路测量误差σ的计算公式如下:
σ=σtPLL+θe
其中,σtPLL为热噪声误差,θe为动态应力误差。
热噪声误差σtPLL的估算公式如下:
其中,BL为环路带宽,C/N0为输入信号的信噪比;T为相干积分时间。
动态应力误差θe的计算公式如下;
其中,ωn为本实施例的系统中环路特征频率,N为环路阶数,R为卫星导航接收机与卫星的距离。
对于二阶环,N=2,则根据动态应力误差θe的计算公式可得:
对于三阶环,N=3,则根据动态应力误差θe的计算公式可得:
为使环路测量误差σ取最小值,将σ=σtPLL+θe对带宽BL求导数,并使得
得到二阶环路的最优带宽,如下:
三阶环路的最优带宽,如下:
高轨航天器的卫星导航接收机根据信噪比、航天器的动态参数以及初始相干积分时间,根据上式可将环路带宽自适应地调节到最优值,从而使环路一直保持在最优跟踪状态。
上述实施例,根据高轨航天器的特点在不同的飞行高度采用不同的相干积分时间,能够提高卫星导航接收机的灵敏度,能够接收到漏导航及旁瓣信号;且飞行过程根据动态与信号强弱,自适应调整跟踪环路的参数,能够满足高轨航天器高动态的要求并能够解决强弱信号并存时的远近效应问题;同时,环路参数的选取结合工程实践经验,简单高效,既满足型号任务需要又有利于接收机的工程化实现,适用于北斗、GPS等多种卫星导航系统的接收机终端,具有较强的实用性。
对应于上文实施例所述的卫星导航接收机的环路跟踪系统,本实施例提供了一种卫星导航接收机的环路跟踪方法。参见图4,为本实施例提供的卫星导航接收机的环路跟踪方法的一个实施例的流程示意图,详述如下:
步骤S401,通过信号采样模块获取导航卫星的下变频信号,将所述下变频信号进行模数转换输出数字中频信号。
步骤S402,获取本地载波和本地伪码,并根据本地载波和本地伪码对所述数字中频信号的I/Q路均进行相干积分。
步骤S403,对所述相干积分后的信号进行滤波鉴相得到码相位误差信号,所述码相位误差信号用于修正码数字振荡器的频率字。
步骤S404,获取航天器的动态参数和所述下变频信号的信噪比,根据所述动态参数、所述信噪比和所述相干积分后的信号进行载波鉴相得到载波相位误差信号,所述载波相位误差信号用于修正载波数字振荡器的频率字,使所述本地载波实时跟踪所述导航卫星。
进一步地,在修正码数字振荡器的频率字之后,还包括:
获取频率字修正后的码数字振荡器的第二伪码,根据所述本地载波和所述第二伪码对所述数字中频信号的Q路进行相干积分,根据所述动态参数、所述信噪比和相干积分后的信号进行载波鉴相得到载波相位误差信号,所述载波相位误差信号用于修正载波数字振荡器的频率字,使所述本地载波实时跟踪所述导航卫星。
可选的,参见图5,本实施例的步骤S404的具体实现流程可以包括:
步骤S501,获取航天器的动态参数、所述下变频信号的信噪比和所述载波环鉴别器的初始相干积分时间。
步骤S502,根据所述动态参数、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽。
步骤S503,根据所述环路阶数和所述环路带宽对载波鉴相后的信号进行滤波得到载波相位误差信号,所述载波相位误差信号用于修正数字振荡器的频率字,使所述本地载波实时跟踪所述导航卫星。
进一步地,所述航天器的动态参数包括:所述航天器的高度、所述航天器的加速度和所述航天器的加加速度;
所述根据所述动态参数、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽,包括:
在所述航天器的高度低于第一高度阈值时,所述相干积分时间确定为第一时间,环路带宽确定为第一频率;
在所述航天器的高度大于第二高度阈值时,所述相干积分时间确定为第二时间,环路带宽确定为第二频率;
在所述航天器的高度在所述第一高度阈值和所述第二高度阈值之间时,根据所述航天器的加速度、所述航天器的加加速度、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽;
其中,所述第一时间小于所述第二时间,所述第一频率大于所述第二频率。
进一步地,所述在所述航天器的高度在所述第一高度阈值和所述第二高度阈值之间时,根据所述航天器的加速度、所述航天器的加加速度、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽,包括:
在所述航天器的高度在所述第一阈值和所述第二阈值之间时,且所述航天器的加速度小于第一加速度阈值,以及所述航天器的加加速度小于第一加加速度阈值时,所述环路阶数确定为二阶环路,并通过
确定二阶环路带宽BL2;其中,C/N0为所述信噪比,T为所述初始相干积分时间,a为所述航天器的加速度;
在所述航天器的高度在所述第一阈值和所述第二阈值之间时,且所述航天器的加速度小于第一加速度阈值,以及所述航天器的加加速度大于第一加加速度阈值时,所述环路阶数确定为三阶环路,并通过
确定三阶环路带宽BL3;其中,C/N0为所述信噪比,T为所述初始相干积分时间,b为所述航天器的加加速度。
上述卫星导航接收机的环路跟踪方法,根据本地载波、本地伪码对数字中频信号的I路和Q路均进行相干积分,根据相干积分的信号进行滤波和鉴相并修正码数字振荡器的频率字,使得本系统能够接收到弱卫星信号,可避免强弱卫星信号共存造成的远近效应问题;然后获取航天器的动态参数和信噪比,根据动态参数、信噪比和相干积分的信号进行载波鉴相得到载波相位误差信号,所述载波相位误差信号用于修正载波数字振荡器的频率字,使本地载波实时跟踪导航卫星,兼顾了航天器的高动态,在近远地轨道均能持续的卫星跟踪。
本领域技术人员可以理解,上述实施例中各步骤的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本发明实施例的实施过程构成任何限定。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围,均应包含在本发明的保护范围之内。
Claims (4)
1.一种卫星导航接收机的环路跟踪系统,其特征在于,包括:信号采样模块、相干积分模块、码环鉴别模块和载波鉴别模块;
所述信号采样模块用于获取导航卫星的下变频信号,并将所述下变频信号进行模数转换输出数字中频信号;
所述相干积分模块用于获取本地载波和本地伪码,并根据本地载波和本地伪码对所述数字中频信号的I/Q路均进行相干积分,相干积分后的信号输出给所述码环鉴别模块和所述载波鉴别模块;
所述码环鉴别模块用于对所述相干积分后的信号进行滤波鉴相得到码相位误差信号,所述码相位误差信号用于修正码数字振荡器的频率字;
所述载波鉴别模块用于获取航天器的动态参数和所述下变频信号的信噪比,根据所述动态参数、所述信噪比和所述相干积分后的信号进行载波鉴相得到载波相位误差信号,所述载波相位误差信号用于修正载波数字振荡器的频率字,使所述本地载波实时跟踪所述导航卫星;
其中,所述载波鉴别模块包括:载波环鉴别器、自适应控制模块、二阶环路滤波器和三阶环路滤波器;所述载波环鉴别器用于对所述相干积分后的信号进行载波鉴相;所述自适应控制模块用于获取航天器的动态参数、信噪比和所述载波环鉴别器的初始相干积分时间,根据所述动态参数、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽;所述二阶环路滤波器用于在所述环路阶数为二阶时,根据所述环路带宽对载波鉴相后的信号进行滤波得到载波相位误差信号,所述载波相位误差信号用于修正数字振荡器的频率字;所述三阶环路滤波器用于在所述环路阶数为三阶时,根据所述环路带宽对载波鉴相后的信号进行滤波得到载波相位误差信号,所述载波相位误差信号用于修正数字振荡器的频率字;
所述航天器的动态参数包括:所述航天器的高度、所述航天器的加速度和所述航天器的加加速度;所述自适应控制模块具体用于:在所述航天器的高度低于第一高度阈值时,所述相干积分时间确定为第一时间,环路带宽确定为第一频率;在所述航天器的高度大于第二高度阈值时,所述相干积分时间确定为第二时间,环路带宽确定为第二频率;在所述航天器的高度在所述第一高度阈值和所述第二高度阈值之间时,根据所述航天器的加速度、所述航天器的加加速度、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽;其中,所述第一时间小于所述第二时间,所述第一频率大于所述第二频率;
所述自适应控制模块具体用于:在所述航天器的高度在所述第一高度阈值和所述第二高度阈值之间时,且所述航天器的加速度小于第一加速度阈值,以及所述航天器的加加速度小于第一加加速度阈值时,所述环路阶数确定为二阶环路,并通过
确定二阶环路带宽BL2;其中,C/N0为所述信噪比,T为所述初始相干积分时间,a为所述航天器的加速度;在所述航天器的高度在所述第一高度阈值和所述第二高度阈值之间时,且所述航天器的加速度小于第一加速度阈值,以及所述航天器的加加速度大于第一加加速度阈值时,所述环路阶数确定为三阶环路,并通过
确定三阶环路带宽BL3;其中,C/N0为所述信噪比,T为所述初始相干积分时间,b为所述航天器的加加速度。
2.如权利要求1所述的卫星导航接收机的环路跟踪系统,其特征在于,所述码环鉴别模块包括:码环鉴别器和码环滤波器;
所述码环鉴别器用于对所述相干积分后的信号进行码环鉴相;
所述码环滤波器用于对码环鉴相后的信号进行滤波得到码相位误差信号,所述码相位误差信号用于修正码数字振荡器的频率字。
3.如权利要求1所述的卫星导航接收机的环路跟踪系统,其特征在于,在修正码数字振荡器的频率字之后,所述相干积分模块还用于:
获取频率字修正后的所述码数字振荡器的第二伪码,根据所述本地载波和所述第二伪码对所述数字中频信号的Q路进行相干积分,相干积分后的信号输出给所述码环鉴别模块和所述载波鉴别模块。
4.一种卫星导航接收机的环路跟踪方法,其特征在于,包括:
通过信号采样模块获取导航卫星的下变频信号,将所述下变频信号进行模数转换输出数字中频信号;
获取本地载波和本地伪码,并根据本地载波和本地伪码对所述数字中频信号的I/Q路均进行相干积分;
对所述相干积分后的信号进行滤波鉴相得到码相位误差信号,所述码相位误差信号用于修正码数字振荡器的频率字;
获取航天器的动态参数和所述下变频信号的信噪比,根据所述动态参数、所述信噪比和所述相干积分后的信号进行载波鉴相得到载波相位误差信号,所述载波相位误差信号用于修正载波数字振荡器的频率字,使所述本地载波实时跟踪所述导航卫星;
所述获取航天器的动态参数和所述下变频信号的信噪比,根据所述动态参数、所述信噪比和所述相干积分后的信号进行载波鉴相得到载波相位误差信号,所述载波相位误差信号用于修正载波数字振荡器的频率字,使所述本地载波实时跟踪所述导航卫星,包括:获取航天器的动态参数、所述下变频信号的信噪比和载波环鉴别器的初始相干积分时间;根据所述动态参数、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽;根据所述环路阶数和所述环路带宽对载波鉴相后的信号进行滤波得到载波相位误差信号,所述载波相位误差信号用于修正数字振荡器的频率字,使所述本地载波实时跟踪所述导航卫星;
所述航天器的动态参数包括:所述航天器的高度、所述航天器的加速度和所述航天器的加加速度;所述根据所述动态参数、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽,包括:在所述航天器的高度低于第一高度阈值时,所述相干积分时间确定为第一时间,环路带宽确定为第一频率;在所述航天器的高度大于第二高度阈值时,所述相干积分时间确定为第二时间,环路带宽确定为第二频率;在所述航天器的高度在所述第一高度阈值和所述第二高度阈值之间时,根据所述航天器的加速度、所述航天器的加加速度、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽;其中,所述第一时间小于所述第二时间,所述第一频率大于所述第二频率;
所述在所述航天器的高度在所述第一高度阈值和所述第二高度阈值之间时,根据所述航天器的加速度、所述航天器的加加速度、所述信噪比和所述初始相干积分时间确定环路阶数和环路带宽,包括:在所述航天器的高度在所述第一高度阈值和所述第二高度阈值之间时,且所述航天器的加速度小于第一加速度阈值,以及所述航天器的加加速度小于第一加加速度阈值时,所述环路阶数确定为二阶环路,并通过
确定二阶环路带宽BL2;其中,C/N0为所述信噪比,T为所述初始相干积分时间,a为所述航天器的加速度;在所述航天器的高度在所述第一高度阈值和所述第二高度阈值之间时,且所述航天器的加速度小于第一加速度阈值,以及所述航天器的加加速度大于第一加加速度阈值时,所述环路阶数确定为三阶环路,并通过
确定三阶环路带宽BL3;其中,C/N0为所述信噪比,T为所述初始相干积分时间,b为所述航天器的加加速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010888948.0A CN112083452B (zh) | 2020-08-28 | 2020-08-28 | 卫星导航接收机的环路跟踪系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010888948.0A CN112083452B (zh) | 2020-08-28 | 2020-08-28 | 卫星导航接收机的环路跟踪系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112083452A CN112083452A (zh) | 2020-12-15 |
CN112083452B true CN112083452B (zh) | 2024-03-01 |
Family
ID=73728160
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010888948.0A Active CN112083452B (zh) | 2020-08-28 | 2020-08-28 | 卫星导航接收机的环路跟踪系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112083452B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113093238B (zh) * | 2021-03-25 | 2024-01-30 | 中国人民解放军国防科技大学 | 载波数控振荡器及其优化方法、导航接收机 |
CN114063122B (zh) * | 2021-12-02 | 2022-08-02 | 武汉大学 | 电推进转移轨道航天器星载gnss在轨实时定轨方法 |
CN114296110A (zh) * | 2021-12-29 | 2022-04-08 | 航天恒星科技有限公司 | 多频点联合载波频率跟踪方法及装置 |
CN115712135A (zh) * | 2022-11-24 | 2023-02-24 | 中国人民解放军63921部队 | 一种gnss/ins矢量深组合的系统及方法 |
CN116170036A (zh) * | 2022-12-30 | 2023-05-26 | 西安空间无线电技术研究所 | 一种连续跳时扩频信号的自适应载波跟踪装置及方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102435999A (zh) * | 2011-10-26 | 2012-05-02 | 浙江理工大学 | Gps接收机基带模块和gps信号捕获及跟踪方法 |
CN104570016A (zh) * | 2014-12-30 | 2015-04-29 | 北京航天科工世纪卫星科技有限公司 | 一种高动态运动载体的北斗信号捕获跟踪接收方法 |
CN105629268A (zh) * | 2015-02-15 | 2016-06-01 | 航天恒星科技有限公司 | 一种环路跟踪方法及系统 |
CN107765270A (zh) * | 2016-08-17 | 2018-03-06 | 中国航空工业集团公司西安飞行自动控制研究所 | 一种基于卡尔曼滤波的卫星导航接收机跟踪环 |
CN111045040A (zh) * | 2019-12-09 | 2020-04-21 | 北京时代民芯科技有限公司 | 一种适用于动态弱信号的卫星导航信号跟踪系统及方法 |
CN111158022A (zh) * | 2019-12-27 | 2020-05-15 | 中国人民解放军军事科学院国防科技创新研究院 | 基于低轨卫星的接收机跟踪方法 |
CN111308520A (zh) * | 2020-04-15 | 2020-06-19 | 上海无线电设备研究所 | 高动态低信噪比下卫星信号的跟踪环路、方法及接收机 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7639181B2 (en) * | 2005-07-01 | 2009-12-29 | Sirf Technology Holdings, Inc. | Method and device for tracking weak global navigation satellite system (GNSS) signals |
US10031236B2 (en) * | 2015-10-06 | 2018-07-24 | Topcon Positioning Systems, Inc. | Navigation receiver with an adaptive system for tracking carrier phases received from a constellation of navigation satellites |
-
2020
- 2020-08-28 CN CN202010888948.0A patent/CN112083452B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102435999A (zh) * | 2011-10-26 | 2012-05-02 | 浙江理工大学 | Gps接收机基带模块和gps信号捕获及跟踪方法 |
CN104570016A (zh) * | 2014-12-30 | 2015-04-29 | 北京航天科工世纪卫星科技有限公司 | 一种高动态运动载体的北斗信号捕获跟踪接收方法 |
CN105629268A (zh) * | 2015-02-15 | 2016-06-01 | 航天恒星科技有限公司 | 一种环路跟踪方法及系统 |
CN107765270A (zh) * | 2016-08-17 | 2018-03-06 | 中国航空工业集团公司西安飞行自动控制研究所 | 一种基于卡尔曼滤波的卫星导航接收机跟踪环 |
CN111045040A (zh) * | 2019-12-09 | 2020-04-21 | 北京时代民芯科技有限公司 | 一种适用于动态弱信号的卫星导航信号跟踪系统及方法 |
CN111158022A (zh) * | 2019-12-27 | 2020-05-15 | 中国人民解放军军事科学院国防科技创新研究院 | 基于低轨卫星的接收机跟踪方法 |
CN111308520A (zh) * | 2020-04-15 | 2020-06-19 | 上海无线电设备研究所 | 高动态低信噪比下卫星信号的跟踪环路、方法及接收机 |
Non-Patent Citations (2)
Title |
---|
Direct position tracking loop based on linearised signal model for global navigation satellite system receivers;Jing Liu 等;《IET Radar, Sonar and Navigation》;20130801;第7卷(第7期);789-799 * |
基于环路相关积分观测的SINS/GPS深组合导航算法;谢非 等;《中国惯性技术学报》;20130815;第21卷(第04期);472-477 * |
Also Published As
Publication number | Publication date |
---|---|
CN112083452A (zh) | 2020-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112083452B (zh) | 卫星导航接收机的环路跟踪系统及方法 | |
US6240367B1 (en) | Full fusion positioning method for vehicle | |
US9513129B2 (en) | Low authority GPS aiding of navigation system for anti-spoofing | |
EP2064568B1 (en) | Highly integrated gps, galileo and inertial navigation system | |
US6311129B1 (en) | Positioning process and system thereof | |
US6449559B2 (en) | Fully-coupled positioning process and system thereof | |
CN111505679B (zh) | 一种基于星载gnss的leo初轨确定方法 | |
CN107656300B (zh) | 基于北斗/gps双模软件接收机的卫星/惯性超紧组合方法 | |
CN111190198B (zh) | Gbas测试设备的卫星接收机及其伪距生成方法与系统 | |
CN111896987A (zh) | 一种低轨导航增强下进行gnss/ins组合导航的方法和装置 | |
US8228229B2 (en) | Space navigation receiver operating with the assistance of navigation information | |
Soloviev et al. | Deeply Integrated GPS/Low‐Cost IMU for Low CNR Signal Processing: Concept Description and In‐Flight Demonstration | |
CN114646315A (zh) | 使用机会信号进行导航的系统、方法以及网络操作中心 | |
Capuano et al. | GNSS/INS/Star tracker integrated navigation system for Earth-Moon transfer orbit | |
US20210055425A1 (en) | Flexible device for synchronizing multi-antenna gnss measurements | |
TW448304B (en) | Fully-coupled positioning process and system | |
RU2187127C2 (ru) | Способ автономного снижения порогов обнаружения и отслеживания несущих сигналов, принимаемых на орбите | |
Ren et al. | A GEO satellite position and beam features estimation method based on beam edge positions | |
WO2022266204A1 (en) | Systems and methods for acquisition and tracking of unknown leo satellite signals | |
Buechler et al. | Integration of GPS and strapdown inertial subsystems into a single unit | |
Tu et al. | Research on satellite signal vector tracking based on prefilter under high dynamic conditions | |
CN106772473B (zh) | 一种实现导航跟踪的方法和装置 | |
EP0944842B1 (en) | High precision hardware carrier frequency and phase aiding in a gps receiver | |
JPH11500231A (ja) | 軌道で受信したスペクトル拡散コードの閾値の収集を自動的に減少し追跡する方法 | |
DiEsposti et al. | The benefits of integrating GPS, INS, and PCS |
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 |