CN108345018A - 一种自适应卫星导航多模式相对定位方法 - Google Patents
一种自适应卫星导航多模式相对定位方法 Download PDFInfo
- Publication number
- CN108345018A CN108345018A CN201711469167.2A CN201711469167A CN108345018A CN 108345018 A CN108345018 A CN 108345018A CN 201711469167 A CN201711469167 A CN 201711469167A CN 108345018 A CN108345018 A CN 108345018A
- Authority
- CN
- China
- Prior art keywords
- difference
- relative positioning
- carrier phase
- pseudo range
- pseudorange
- 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.)
- Pending
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
Abstract
本发明提供了一种自适应卫星导航多模式相对定位方法,按照距离段将相对定位分为伪距差分、载波相位平滑码伪距差分和载波相位差分模式。根据三种定位模式自身的精度水平以及不同引导阶段的距离要求,将伪距差分定位模式作用距离划分为40km~20km,载波相位平滑码伪距差分定位模式作用距离划分为20~5km,载波相位差分定位模式作用距离划分为5~0km,在切换时对两种不同模式的定位结果进行平滑滤波处理,保证定位结果的稳定性。本发明不但减轻了算法运算压力,而且能够提升整个引导过程中相对定位结果的稳定性,具有较高的工程应用价值。
Description
技术领域
本发明属于卫星导航领域,是卫星导航相对定位系统中一种自适应多模式相对定位方法。
背景技术
卫星导航相对定位系统是一种基于卫星导航差分相对定位技术原理的系统,可用于精密着舰、飞船交会对接、空中自主加油等应用中,能够为这些应用提供用以引导的相对位置信息。
卫星导航相对定位系统按照技术分类可分为伪距差分、载波相位平滑伪距差分、载波相位差分技术。在三种相对定位技术中,伪距差分简单可靠,无需载波相位测量信息,不需求解整周模糊度,不受周跳的影响,精度稍低;载波相位平滑伪距差分不需求解整周模糊度,但是受周跳的影响,精度中等;载波相位差分需求解整周模糊度,运算复杂,其求解不稳定易受周跳的影响,但精度最高。
对于精密着舰、飞船交会对接、空中自主加油等应用,引导过程中的不同距离段对相对定位精度有着不同的要求,距离越近精度要求越高。因此,如果在引导全过程中选择伪距差分模式进行相对引导,由于其相对定位精度不高,在引导最后阶段不能够满足相对引导要求;如果选择载波相位差分模式进行相对引导,虽然相对定位精度最高,但增加了运算复杂性且全过程中定位求解不稳定。
目前,针对精密着舰、飞船交会对接、空中自主加油等应用所提出的相对定位精度指标仅为最后引导阶段的最高精度指标要求,对于从远距离引导到近距离的引导全过程中使用的定位模式没有统一的要求。因此,需要在满足相对定位精度要求的同时,根据距离远近合理选择相对定位模式。
发明内容
为了克服现有技术的不足,本发明提供一种自适应卫星导航多模式相对定位方法,根据引导距离的远近合理的选择相对定位方式,能够在满足相对定位精度要求的条件下,有效提升定位稳定性,减轻算法运算压力。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤一,进行伪距差分相对定位,将基准站的伪距观测信息发给移动站,然后建立双差伪距观测方程,通过最小二乘求得基准站与移动站的基线矢量;
步骤二,伪距差分与载波相位平滑伪距差分相对定位间切换,在模式切换前进行设定时长的滤波处理,在切换过程中以滤波处理后的定位结果作为相对定位结果值;
步骤三,进行载波平滑码伪距差分相对定位,依据前一时刻伪距双差平滑值,以及前一时刻和当前时刻的载波相位双差值,预估当前时刻伪距双差值;然后,由当前时刻伪距双差观测值和预估值,计算得到当前时刻伪距双差平滑值;最后,利用当前时刻伪距双差平滑值求解当前时刻动态站位置;
步骤四,载波平滑码伪距差分与载波相位差分相对定位间切换,在模式切换前进行设定时长的滤波处理,在切换过程中以滤波处理后的定位结果作为相对定位结果值;
步骤五,进行载波相位差分相对定位,基准站测量所有可见卫星的载波相位,并通过通信数据链实时将观测的载波相位发送给移动站;移动站同时测量可见卫星的载波相位,选定某一颗卫星作为参考卫星,并对同一颗卫星、相同观测历元的载波相位观测值进行站间星间求差,得到载波相位双差观测值;利用载波相位双差观测方程求得精密相对定位结果。
本发明的有益效果是:根据引导距离的远近选择差分定位方式,在为整个引导过程提供连续、满足精度要求的相对定位信息的同时,不但减轻了算法运算压力,而且能够提升整个引导过程中相对定位结果的稳定性,具有较高的工程应用价值。
附图说明
图1是自适应卫星导航多模式相对定位方法总体框图;
图2是载波相位平滑伪距差分相对定位方法流程图。
具体实施方式
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
本发明提供的自适应卫星导航多模式相对定位方法按照距离段将相对定位分为伪距差分、载波相位平滑码伪距差分和载波相位差分模式。根据三种定位模式自身的精度水平以及不同引导阶段的距离要求,将伪距差分定位模式作用距离划分为40km~20km,载波相位平滑码伪距差分定位模式作用距离划分为20~5km,载波相位差分定位模式作用距离划分为5~0km,总体框图如图1所示。
由于精密着舰、飞船交会对接、空中自主加油等应用的引导都是由远至近的过程,因此自适应卫星导航多模式相对定位方法在应用过程中首先进入伪距差分模式,其次进入载波相位平滑伪距差分模式,最后进入载波相位差分模式。同时由于三种相对定位模式定位精度差别较大,其中伪距差分相对定位精度为米量级,载波相位平滑伪距差分相对定位精度为亚米量级,载波相位差分相对定位精度为厘米量级。为保证在不同模式切换时相对定位结果不发生大的跳变,须在切换时对两种不同模式的定位结果进行平滑滤波处理,保证定位结果的稳定性。自适应卫星导航多模式相对定位方法共分为5个步骤,分别为步骤1、伪距差分相对定位;步骤2、伪距差分与载波平滑码伪距差分相对定位间切换;步骤3、载波平滑码伪距差分相对定位;步骤4、载波平滑码伪距差分与载波相位差分相对定位间切换;步骤5、载波相位差分相对定位。
本发明是一种根据引导距离的远近合理的选择相对定位方式的自适应卫星导航多模式相对定位方法,具体步骤如下。
步骤一:伪距差分相对定位
将移动基准站的伪距观测信息发给移动站,然后建立双差伪距观测方程,通过最小二乘可求得基准站与移动站的基线矢量。伪距双差定位时需要基准站的单点定位坐标以及卫星坐标,必须先进行伪距单点定位。伪距单点定位模型为:
其中上标j为卫星编号;下标1为观测站编号;为观测站1至j卫星的伪距测量值;为j卫星的位置矢量,已知;为测站1的位置矢量,未知;c为真空中的光速;δt1为观观测站1的钟差,未知;δtj为j卫星钟差,已知;为电离层延迟误差;为对流层延迟误差;为测量误差。
观测n(n≥4)颗卫星可列n个观测方程,采用高斯-牛顿迭代最小二乘方法求得观测站的位置。
然后根据观测站的位置及卫星坐标,可建立伪距双差观测方程。移动基准站1跟踪参考星的编号为k,同时移动站2也跟踪到这颗卫星,在40Km基线内,双差处理后可以认为电离层延迟误差、对流层延迟误差和星历误差消除或削弱得很充分,伪距双差观测方程如下:
其中,为观测站1号、2号的j卫星和i卫星双差伪距观测信息,为伪距测量误差,其他符号同上。
将移动基准站1的伪距单点定位结果作为已知值,移动站2伪距单点定位结果作为初始值线性式,有
式中,Yρ为双差伪距与计算值的差值,dX为基线改正矢量;r120为基线矢量初值;A为线性化的双差几何矩阵。
观测n颗卫星,可列(n-1)个伪距双差观测方程,根据最小二乘,可得基线矢量改正数为dX=(ATA)-1ATY,则伪距差分基线矢量为r12=r120+dX。
步骤二:伪距差分与载波平滑码伪距差分相对定位间切换
由于伪距差分模式与载波平滑码伪距差分模式定位精度差别较大,直接由伪距差分模式切换至载波平滑码伪距模式会造成相对定位结果较大的跳变。因此在模式切换前进行一段时间的滤波处理,在切换过程中以滤波处理后的定位结果作为相对定位结果值。根据模式作用距离划分当步骤一r12基线距离达到25km开始启动滤波处理,滤波处理计算公式如下:
式中,K(k)称为衰减记忆滤波增益因子;为本时刻平滑位置,为上一时刻的平滑位置;Pdρ(tk)为当前时刻伪距差分定位结果;δP(tk-1,tk)为前后时刻载波平滑码伪距差分相对定位结果Pdcρ(tk-1)、Pdcρ(tk)的变化值。k为平滑的历元个数,初值取0。
步骤三:载波平滑码伪距差分相对定位
当载波跟踪连续,且前后历元移动基准站位置变化不大时,载波相位测量能给出两个测量历元之间精度很高的伪距变化量。可以用此伪距变化量来修正当前历元的伪距观测量,得到更加准确的当前历元的伪距观测量。站间星间伪距双差、载波相位双差可以很好地消除或削弱星历误差、卫星钟差、接收机钟差的影响。由于载波相位的测量噪声比伪距测量噪声小,同时,多路径效应对载波相位的影响要小于对伪距的影响,因此,载波相位差分平滑伪距差分相对定位方法有效地消除了多径效应和测量噪声对伪距的影响,提高了导航精度。
载波相位差分平滑伪距差分相对定位算法如图2所示。依据前一时刻伪距双差平滑值,以及前一时刻和当前时刻载波相位双差值,预估当前时刻伪距双差值。然后,由当前时刻伪距双差观测值和预估值,计算得到当前时刻伪距双差平滑值。最后,利用当前时刻伪距双差平滑值求解当前时刻动态站位置。平滑伪距计算公式如下:
式中,K(k)称为衰减记忆滤波增益因子;为上一时刻的平滑双差伪距;为当前时刻双差伪距;δφ(tk-1,tk)为前后时刻载波相位双差的变化;k为平滑的历元个数,初值取0。
步骤四:载波平滑码伪距差分与载波相位差分相对定位间切换
载波平滑码伪距差分模式与载波相位差分模式定位精度差别较大,直接由载波平滑码伪距差分模式切换至载波相位差分模式会造成相对定位结果较大的跳边。因此在模式切换前进行一段时间的滤波处理,在切换过程中以滤波处理后的定位结果作为相对定位结果值。根据模式作用距离划分当步骤三基线距离达到10km开始启动滤波处理,滤波处理计算公式如下:
式中,K(k)称为衰减记忆滤波增益因子;为本时刻平滑位置,为上一时刻的平滑位置;Pdcρ(tk)为当前时刻载波平滑码伪距差分定位结果;δP(tk-1,tk)为前后时刻载波相位差分相对定位结果的变化值。k为平滑的历元个数,初值取0。
步骤五:载波相位差分相对定位
将移动基准站测量所有可见卫星的载波相位,并通过通信数据链实时将观测的载波相位发送给移动站。移动站同时测量可见卫星的载波相位,选定某一颗卫星作为参考卫星,并对同一颗卫星、相同观测历元的载波相位观测值进行站间星间求差,得到载波相位双差观测值。在实时正确得到整周模糊度的前提下,利用载波相位双差观测方程即可求得精密相对定位结果。载波相位双差观测方程如下:
其中,为移动基准站1,移动站2号的j卫星和i卫星双差载波相位观测信息,为载波噪声误差,其他符号同步骤一中定义。将移动基准站1的伪距单点定位结果作为已知值,移动站2伪距单点定位结果作为初始值线性式,有
式中,Yφ为双差载波相位与计算值的差值,dX为基线改正矢量;r120为基线矢量初值;A为线性化的双差几何矩阵。
观测n颗卫星,可列(n-1)个载波相位双差观测方程,根据最小二乘,可得基线矢量改正数为dX=(ATA)-1ATY,则载波相位差分基线矢量为r12=r120+dX。
Claims (1)
1.一种自适应卫星导航多模式相对定位方法,其特征在于包括下述步骤:
步骤一,进行伪距差分相对定位,将基准站的伪距观测信息发给移动站,然后建立双差伪距观测方程,通过最小二乘求得基准站与移动站的基线矢量;
步骤二,伪距差分与载波相位平滑伪距差分相对定位间切换,在模式切换前进行设定时长的滤波处理,在切换过程中以滤波处理后的定位结果作为相对定位结果值;
步骤三,进行载波平滑码伪距差分相对定位,依据前一时刻伪距双差平滑值,以及前一时刻和当前时刻的载波相位双差值,预估当前时刻伪距双差值;然后,由当前时刻伪距双差观测值和预估值,计算得到当前时刻伪距双差平滑值;最后,利用当前时刻伪距双差平滑值求解当前时刻动态站位置;
步骤四,载波平滑码伪距差分与载波相位差分相对定位间切换,在模式切换前进行设定时长的滤波处理,在切换过程中以滤波处理后的定位结果作为相对定位结果值;
步骤五,进行载波相位差分相对定位,基准站测量所有可见卫星的载波相位,并通过通信数据链实时将观测的载波相位发送给移动站;移动站同时测量可见卫星的载波相位,选定某一颗卫星作为参考卫星,并对同一颗卫星、相同观测历元的载波相位观测值进行站间星间求差,得到载波相位双差观测值;利用载波相位双差观测方程求得精密相对定位结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711469167.2A CN108345018A (zh) | 2017-12-29 | 2017-12-29 | 一种自适应卫星导航多模式相对定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711469167.2A CN108345018A (zh) | 2017-12-29 | 2017-12-29 | 一种自适应卫星导航多模式相对定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108345018A true CN108345018A (zh) | 2018-07-31 |
Family
ID=62963359
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711469167.2A Pending CN108345018A (zh) | 2017-12-29 | 2017-12-29 | 一种自适应卫星导航多模式相对定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108345018A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109541658A (zh) * | 2018-12-27 | 2019-03-29 | 中国电子科技集团公司第二十研究所 | 一种用于卫星导航相对定位的虚拟点融合观测量生成方法 |
CN109633616A (zh) * | 2018-12-11 | 2019-04-16 | 中国人民解放军63921部队 | 一种基于非相干测量的载波平滑伪码测距装置及方法 |
CN110687559A (zh) * | 2019-11-04 | 2020-01-14 | 中国电子科技集团公司第五十四研究所 | 适用于机载的gnss无缝高精度定位和完好性评估方法 |
CN111077550A (zh) * | 2019-12-26 | 2020-04-28 | 广东星舆科技有限公司 | 一种应用于智能终端rtd定位的粗差探测方法及系统 |
CN112987068A (zh) * | 2019-12-16 | 2021-06-18 | 中国电子科技集团公司第二十研究所 | 一种基于测距信息辅助的卫星导航高精度相对定位方法 |
CN113009537A (zh) * | 2021-02-18 | 2021-06-22 | 中国人民解放军国防科技大学 | 一种惯导辅助卫导相对定位单历元部分模糊度求解方法 |
CN113419263A (zh) * | 2021-05-25 | 2021-09-21 | 武汉导航与位置服务工业技术研究院有限责任公司 | 一种分离模糊度与位置解算的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104459737A (zh) * | 2014-12-02 | 2015-03-25 | 东南大学 | 一种基于信标实时差分的卫星定位方法 |
CN105242292A (zh) * | 2015-10-30 | 2016-01-13 | 中国电子科技集团公司第二十研究所 | 一种长基线的伪距差分定位方法 |
CN107229061A (zh) * | 2017-07-18 | 2017-10-03 | 武汉大学 | 一种基于低轨卫星的星地差分实时精密定位方法 |
-
2017
- 2017-12-29 CN CN201711469167.2A patent/CN108345018A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104459737A (zh) * | 2014-12-02 | 2015-03-25 | 东南大学 | 一种基于信标实时差分的卫星定位方法 |
CN105242292A (zh) * | 2015-10-30 | 2016-01-13 | 中国电子科技集团公司第二十研究所 | 一种长基线的伪距差分定位方法 |
CN107229061A (zh) * | 2017-07-18 | 2017-10-03 | 武汉大学 | 一种基于低轨卫星的星地差分实时精密定位方法 |
Non-Patent Citations (2)
Title |
---|
范利涛 等: "GPS差分相对定位应用于航天器自主交会对接的研究", 《飞行器测控学报》 * |
范利涛 等: "载波相位差分相对定位在航天器对接中的应用 ", 《中国惯性技术学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109633616A (zh) * | 2018-12-11 | 2019-04-16 | 中国人民解放军63921部队 | 一种基于非相干测量的载波平滑伪码测距装置及方法 |
CN109541658A (zh) * | 2018-12-27 | 2019-03-29 | 中国电子科技集团公司第二十研究所 | 一种用于卫星导航相对定位的虚拟点融合观测量生成方法 |
CN110687559A (zh) * | 2019-11-04 | 2020-01-14 | 中国电子科技集团公司第五十四研究所 | 适用于机载的gnss无缝高精度定位和完好性评估方法 |
CN110687559B (zh) * | 2019-11-04 | 2021-07-13 | 中国电子科技集团公司第五十四研究所 | 适用于机载的gnss无缝高精度定位和完好性评估方法 |
CN112987068A (zh) * | 2019-12-16 | 2021-06-18 | 中国电子科技集团公司第二十研究所 | 一种基于测距信息辅助的卫星导航高精度相对定位方法 |
CN112987068B (zh) * | 2019-12-16 | 2023-11-10 | 中国电子科技集团公司第二十研究所 | 一种基于测距信息辅助的卫星导航高精度相对定位方法 |
CN111077550A (zh) * | 2019-12-26 | 2020-04-28 | 广东星舆科技有限公司 | 一种应用于智能终端rtd定位的粗差探测方法及系统 |
CN113009537A (zh) * | 2021-02-18 | 2021-06-22 | 中国人民解放军国防科技大学 | 一种惯导辅助卫导相对定位单历元部分模糊度求解方法 |
CN113009537B (zh) * | 2021-02-18 | 2023-10-31 | 中国人民解放军国防科技大学 | 一种惯导辅助卫导相对定位单历元部分模糊度求解方法 |
CN113419263A (zh) * | 2021-05-25 | 2021-09-21 | 武汉导航与位置服务工业技术研究院有限责任公司 | 一种分离模糊度与位置解算的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108345018A (zh) | 一种自适应卫星导航多模式相对定位方法 | |
CN104536027B (zh) | 一种实时北斗精密相对定位的方法 | |
CN107607934B (zh) | 一种时差、频差、频差变化率联合估计方法 | |
CN101915928B (zh) | 双星时差/频差联合定位的方法及装置 | |
CN108363084A (zh) | 利用卫星定位的方法和装置、卫星导航接收机、存储介质 | |
CN102590840B (zh) | 一种卫星定位载波相位差分方法 | |
CN106054214B (zh) | 一种可削弱多路径效应的gnss数据处理方法 | |
CN104199280B (zh) | 一种基于差分gps的时间同步误差测量方法 | |
CN102269819B (zh) | 一种用于精密单点定位方法的估计方法 | |
CN106896386B (zh) | Glonass频间偏差精确估计方法 | |
CN101943749A (zh) | 基于星型结构的虚拟参考站网络rtk定位方法 | |
CN108196267B (zh) | 一种基于gnss cp技术的不间断时间传递方法 | |
CN105044741B (zh) | 一种伪距相位综合广域差分改正值的求解方法 | |
CN107193023B (zh) | 一种具有闭式解的高精度北斗卫星系统单点定位方法 | |
CN105242292A (zh) | 一种长基线的伪距差分定位方法 | |
CN106646564A (zh) | 一种基于低轨卫星增强导航方法 | |
CN104656108A (zh) | 一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法 | |
CN104730551B (zh) | 一种星地双基地差分干涉基线坐标以及形变量测量方法 | |
CN104503223B (zh) | Gnss 三频高精度卫星钟差估计与服务方法 | |
CN105700000A (zh) | 一种北斗导航接收机实时动态精密定位方法 | |
CN108427132A (zh) | 一种多频多系统网络rtk基准站模糊度的快速确定方法 | |
CN101614802A (zh) | 一种导航卫星姿态测量方法 | |
CN109507704A (zh) | 一种基于互模糊函数的双星定位频差估计方法 | |
Shao et al. | Study on differential GPS positioning methods | |
CN106597492B (zh) | 卫星导航接收机及其抗远近效应的方法和室内定位方法 |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180731 |