CN104309817A - 基于多台并址接收机的北斗导航卫星区域定轨方法 - Google Patents
基于多台并址接收机的北斗导航卫星区域定轨方法 Download PDFInfo
- Publication number
- CN104309817A CN104309817A CN201410535800.3A CN201410535800A CN104309817A CN 104309817 A CN104309817 A CN 104309817A CN 201410535800 A CN201410535800 A CN 201410535800A CN 104309817 A CN104309817 A CN 104309817A
- Authority
- CN
- China
- Prior art keywords
- orbit determination
- carrier phase
- baseline
- satellite
- beidou navigation
- 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.)
- Granted
Links
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供了一种基于多台并址接收机的北斗导航卫星区域定轨方法,针对中国区域网跟踪站,每个站点选取多台接收机数据,根据共视观测量最多原则生成中长基线和超短基线;基于上述中长基线和超短基线的单差数据,生成双差数据进行初步的卫星精密定轨,得到模糊度浮点解和估计参数;跟踪站点坐标,对于不同长度的基线采用不同的方法解算载波相位模糊度;利用生成的载波相位模糊度重新进行参数估计解算卫星轨道。本发明提出利用并址放置的多台接收机进行北斗导航卫星区域定轨,可以提高载波相位模糊度解算成功率,提高载波相位观测值在北斗导航卫星区域定轨的价值,最终提高北斗导航卫星定轨精度。
Description
技术领域
本发明涉及一种北斗导航卫星区域定轨方法,涉及卫星精密定轨领域,能够显著提高区域布站条件下的北斗导航卫星定轨精度。
背景技术
鉴于卫星导航***在军事、经济、政治领域的巨大作用,我国在积极发展自己的卫星导航***。1994年,北斗一号卫星导航***正式启动。该***包括3颗卫星,均为地球静止轨道卫星。与美国GPS相比,北斗一号覆盖范围小,定位精度低、无法在高速移动平台上使用。为了突破用户数量限制、提高***生存能力、扩大***覆盖范围、实现与GPS、Galileo等国外卫星导航***的兼容互操作等,中国考虑建设北斗卫星导航定位***(COMPASS***)。2004年,中国正式启动了具有全球导航能力的北斗卫星导航***的建设。并于2007年发射一颗中地球轨道卫星,并进行了大量试验。2009年起,后续卫星持续发射,并在2011年开始对中国和周边地区提供测试服务。2012年12月27日,北斗***空间信号接口控制文件正式版正式公布,北斗导航卫星***正式对亚太地区提供无源定位、导航、授时服务。截至目前,北斗卫星导航***已经发射14颗卫星,其中5颗GEO卫星、5颗IGSO卫星、4颗MEO卫星。
与GPS等其他卫星导航***不同,北斗卫星导航***发播B1(1561.098MHZ)、B2(1207.14MHZ)、B3(1268.52MHZ)三个频点数据。每个频点信号由I、Q支路的“测距码+导航电文”正交调制在载波上构成。其中I支路信号为民用信号,Q支路为军用信号。民用用户实际可以接受三个频点的载波相位观测值和伪距观测值。其中,伪距观测值精度为0.3m,载波相位观测值精度为毫米级。因此,相比伪距观测值,载波相位观测值对于北斗卫星导航***的高精度应用具有重要价值。但是由于载波相位观测值中存在整周模糊度项,因此整周模糊度的解算成为制约载波相位观测值应用的主要因素。
目前,IGS(International GNSS Service)的跟踪站已经实现全球分布,此外卫星轨道动力学模型也不断完善,IGS提供的产品的精度也不断提高。其中,IGS事后最终轨道精度可以达到5cm,事后最终钟差精度为0.1~0.2ns,快速精密星历的精度为5cm,快速实时预报星历的精度为10cm。与IGS不同,限于我国国情,完全可控的COMPASS***跟踪站只能在中国境内布设。相比跟踪站全球均匀分布,国内跟踪站对北斗卫星的可观测弧段较少,阻碍了跟踪站对卫星的全弧段跟踪。此外,虽然国内跟踪站均匀分布,但是由于我国地域广泛,站间距离相距较远,不利于双差模糊度解算,降低了载波相位模糊度的解算成功率,进而影响高精度载波相位观测值的应用价值。
发明内容
为了克服现有技术的不足,本发明提供一种基于多台并址接收机的北斗导航卫星区域定轨方法,在区域网跟踪站上安放了多套北斗接收机。该方案针对中国区域跟踪站,每个站点至少选取2台接收机数据。轨道确定过程中,由于站点并址放置了多套接收机,采用双差模式进行数据解算时可以形成超短基线,大幅度提高模糊度固定率了,提高了载波相位观测值的应用价值,进而提高卫星精密定轨精度。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤1:针对中国区域网跟踪站,每个站点选取至少两台接收机数据;根据共视观测量最多原则,生成中长基线和超短基线;
步骤2:基于上述中长基线和超短基线的单差数据,生成双差数据进行初步的卫星精密定轨,得到模糊度浮点解和估计参数,估计参数包括6个轨道根数、9个光压系数、地球自转参数、对流层参数以及跟踪站站坐标;
步骤3:采用步骤2生成的模糊度浮点解和对流层参数,跟踪站点坐标,对于不同长度的基线采用不同的方法解算载波相位模糊度,其中,对于基线长度小于等于10km的基线采用Sigma策略,对大于10km且小于2000km的基线采用QIF策略,对大于2000km的基线采用宽巷技术;
步骤4:利用步骤3生成的载波相位模糊度,重新进行参数估计解算卫星轨道。
本发明的有益效果是:采用一种区域网跟踪站布设方案进行北斗卫星精密定轨,利用新的区域网进行精密定轨时可以形成超短基线,解决了区域跟踪网模糊度固定率低的问题,能够明显的提高双差模糊度固定率,从而提高定轨精度。在原有区域网布站基础下,采用本发明可以明显提高卫星定轨精度。
附图说明
图1是区域网跟踪站布站示意图;
图2是北斗卫星精密定轨数据处理流程图。
具体实施方式
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
本发明包括以下步骤:
步骤1:针对中国区域网跟踪站,每个站点选取至少两台接收机数据。根据共视观测量最多原则,除了生成多条中长基线外,还可以形成多条超短基线。
步骤2:基于上述中长基线和超短基线的单差数据,生成双差数据进行初步的卫星精密定轨。主要目的是进行模糊度浮点解算。其他估计参数有:6个轨道根数、9个光压系数(其中D、Y常数项两个方向加1e-12大小约束)、地球自转参数(强约束)、对流层参数,以及跟踪站站坐标。D、Y、X为太阳光压模型中定义的三个单位矢量,其中D为太阳指向卫星方向的单位矢量,Y为太阳帆板的旋转轴方向,垂直于与D以及地心指向卫星方向的矢量决定的平面。
步骤3:采用步骤2生成的模糊度浮点解、对流层参数,跟踪站站坐标,对于不同长度的基线采用不同的方法进行模糊度固定。其中,对于基线长度小于等于10km的为超短基线采用Sigma策略,对大于10km且小于2000km的基线采用QIF策略,对大于2000km的基线采用宽巷技术解算模糊度。
步骤4:利用步骤3生成的载波相位模糊度,重新进行参数估计解算卫星轨道。参数估计时引入步骤3生成的载波相位模糊度,估计参数只有:6个轨道根数、9个光压系数(其中D、Y、X三个方向加1e-12大小约束)、地球自转参数(加强约束)、对流层参数。
本发明实施例的技术特征在于:突破了北斗区域网跟踪站布设条件对卫星定轨的局限性。在新的区域跟踪站布站方案下进行卫星精密定轨过程中,可以生成多条超短基线。利用超短基线可以明显提高模糊度固定率,进而可以提高北斗卫星定轨精度。
该方法的实施基于以下条件:1.一个站点并址安放两台北斗接收机;2.卫星精密定轨时采用双差数据模式,利用新的区域跟踪站布设方案,可以生成多条超短基线,这样可以提高模糊度固定率。
本发明包括以下步骤:
步骤1.选取区域网跟踪站接收机数据,生成双差观测文件。
非差观测方程:
其中:为站点K、L与卫星i的非差载波相位观测值;分别为站点K、L至卫星i的几何距离;δkδl为站点K、L的接收机钟差;δi为卫星i的卫星钟差;为卫星i至站点K、L传播路径上的电离层时延;为卫星i至站点K、L传播路径上的对流层时延;为F频点站点K、L接收卫星i的载波相位观测值的整周模糊度;λF为F频点的波长;c为光速;
由非差观测值,进行站间、星间单差生成非差观测值。观测方程如下:
其中:为F频点站点K、L与卫星i、j之间的双差载波相位观测值;为双差电离层延迟;为双差对流层延迟;为双差几何距离项;为双差载波相位模糊度;
由公式(1)(2)(3)可知,由非差观测值生成双差观测值后,卫星钟误差δi、δj、接收机钟误差项均消除。双差观测方程中只有电离层延迟对流层延迟
参看图1,现有北斗区域跟踪网已经广泛分布于中国境内,但是由于中国地域广阔,采用双差模式进行数据处理时,基线长度一般比较大。由于对流层延迟、电离层延迟空间相关度随基线长度的增加而减弱,因此双差观测值中电离层延迟、对流层延迟的影响仍然比较大,很大程度上影响了模糊度的解算。
针对中国区域网布站条件,每个站点选取多台北斗接收机数据进行数据处理。根据共视观测量最多的原则,生成的基线中包含多个超短基线。由于对流层延迟主要与传播路径上的大气压、温度、湿度相关;电离层延迟与信号传播路径上的电子含量有关,因此超短基线形成的双差观测值中对流层延迟、电离层延迟基本消除,有利于载波相位模糊度的解算。
步骤2,基于上述双差文件,数据处理进行模糊度浮点解算。
根据成熟的GPS模糊度解算理论,首先在卫星轨道确定过程中,将模糊度同卫星轨道根数、光压参数、对流层以及其他参数一并解算(按照加权最小二乘法进行解)。其中模糊度是浮点解。
步骤3.模糊度固定
将步骤2进行模糊度浮点解生成的对流层解算结果、卫星轨道等参数作为已知条件,针对不同的基线采用不同的方法进行模糊度固定。其中,由于跟踪网中各站点并址安放了两天接收机。因此,步骤一生成的单差文件中除了中长基线外(1000km~2000km),还包括有多条基线长度小于10km的超短基线。解算模糊度时,根据基线长度的不同,分别选取不同的模糊度解算方法。其中,对于超短基线采用成熟的Sigma方法,对于中长基线采用QIF方法。模糊度解算后保存在观测文件中。
步骤4.模糊度固定解情况下的,北斗卫星精密定轨。
利用步骤3生成的模糊度固定解,进行最终的卫星精密定轨数据解算。其中,模糊度做为已知值引入,其他待估参数有:6个轨道根数、9个光压系数(其中D、Y两个方向加1e-12大小约束)、地球自转参数(加强约束)、对流层参数。
Claims (1)
1.一种基于多台并址接收机的北斗导航卫星区域定轨方法,其特征在于包括下述步骤:
步骤1:针对中国区域网跟踪站,每个站点选取至少两台接收机数据;根据共视观测量最多原则,生成中长基线和超短基线;
步骤2:基于上述中长基线和超短基线的单差数据,生成双差数据进行初步的卫星精密定轨,得到模糊度浮点解和估计参数,估计参数包括6个轨道根数、9个光压系数、地球自转参数、对流层参数以及跟踪站站坐标;
步骤3:采用步骤2生成的模糊度浮点解和对流层参数,跟踪站点坐标,对于不同长度的基线采用不同的方法解算载波相位模糊度,其中,对于基线长度小于等于10km的基线采用Sigma策略,对大于10km且小于2000km的基线采用QIF策略,对大于2000km的基线采用宽巷技术;
步骤4:利用步骤3生成的载波相位模糊度,重新进行参数估计解算卫星轨道。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410535800.3A CN104309817B (zh) | 2014-10-11 | 2014-10-11 | 基于多台并址接收机的北斗导航卫星区域定轨方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410535800.3A CN104309817B (zh) | 2014-10-11 | 2014-10-11 | 基于多台并址接收机的北斗导航卫星区域定轨方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104309817A true CN104309817A (zh) | 2015-01-28 |
CN104309817B CN104309817B (zh) | 2016-08-24 |
Family
ID=52365181
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410535800.3A Expired - Fee Related CN104309817B (zh) | 2014-10-11 | 2014-10-11 | 基于多台并址接收机的北斗导航卫星区域定轨方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104309817B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107678048A (zh) * | 2017-10-10 | 2018-02-09 | 唐山学院 | 卫星钟的超快钟差预报模型的预先选择方法 |
CN110361755A (zh) * | 2019-03-12 | 2019-10-22 | 中国矿业大学 | 一种基于oedop因子的多卫星导航***监测站优化选取方法 |
CN112782710A (zh) * | 2020-12-29 | 2021-05-11 | 中国人民解放军63921部队 | 基于干涉测量体制的低轨空间目标跟踪方法及*** |
CN113581501A (zh) * | 2021-08-27 | 2021-11-02 | 重庆两江卫星移动通信有限公司 | 一种适用于组网低轨卫星联合定轨的***及方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5543804A (en) * | 1994-09-13 | 1996-08-06 | Litton Systems, Inc. | Navagation apparatus with improved attitude determination |
KR20030042543A (ko) * | 2001-11-23 | 2003-06-02 | 한국전자통신연구원 | 위성의 정밀 궤도 결정 시스템 및 그 방법 |
CN101435863A (zh) * | 2008-12-25 | 2009-05-20 | 武汉大学 | 一种导航卫星实时精密定轨的方法 |
CN103363994A (zh) * | 2012-04-09 | 2013-10-23 | 陈刘成 | 仅基于无线电载波相位观测的卫星精密定轨技术 |
CN103675861A (zh) * | 2013-11-18 | 2014-03-26 | 航天恒星科技有限公司 | 一种基于星载gnss多天线的卫星自主定轨方法 |
CN104048664A (zh) * | 2014-07-01 | 2014-09-17 | 南京航空航天大学 | 一种导航卫星星座自主定轨的方法 |
-
2014
- 2014-10-11 CN CN201410535800.3A patent/CN104309817B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5543804A (en) * | 1994-09-13 | 1996-08-06 | Litton Systems, Inc. | Navagation apparatus with improved attitude determination |
KR20030042543A (ko) * | 2001-11-23 | 2003-06-02 | 한국전자통신연구원 | 위성의 정밀 궤도 결정 시스템 및 그 방법 |
CN101435863A (zh) * | 2008-12-25 | 2009-05-20 | 武汉大学 | 一种导航卫星实时精密定轨的方法 |
CN103363994A (zh) * | 2012-04-09 | 2013-10-23 | 陈刘成 | 仅基于无线电载波相位观测的卫星精密定轨技术 |
CN103675861A (zh) * | 2013-11-18 | 2014-03-26 | 航天恒星科技有限公司 | 一种基于星载gnss多天线的卫星自主定轨方法 |
CN104048664A (zh) * | 2014-07-01 | 2014-09-17 | 南京航空航天大学 | 一种导航卫星星座自主定轨的方法 |
Non-Patent Citations (2)
Title |
---|
***,郝金明,李建文等: "一种北斗卫星精密定轨方法", 《测绘科学技术学报》 * |
施闯,赵齐乐,李敏等: "北斗卫星导航***的精密定轨与定位研究", 《中国科学:地球科学》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107678048A (zh) * | 2017-10-10 | 2018-02-09 | 唐山学院 | 卫星钟的超快钟差预报模型的预先选择方法 |
CN107678048B (zh) * | 2017-10-10 | 2021-01-15 | 唐山学院 | 卫星钟的超快钟差预报模型的预先选择方法 |
CN110361755A (zh) * | 2019-03-12 | 2019-10-22 | 中国矿业大学 | 一种基于oedop因子的多卫星导航***监测站优化选取方法 |
CN110361755B (zh) * | 2019-03-12 | 2023-04-07 | 中国矿业大学 | 一种基于oedop因子的多卫星导航***监测站优化选取方法 |
CN112782710A (zh) * | 2020-12-29 | 2021-05-11 | 中国人民解放军63921部队 | 基于干涉测量体制的低轨空间目标跟踪方法及*** |
CN112782710B (zh) * | 2020-12-29 | 2022-10-14 | 中国人民解放军63921部队 | 基于干涉测量体制的低轨空间目标跟踪方法及*** |
CN113581501A (zh) * | 2021-08-27 | 2021-11-02 | 重庆两江卫星移动通信有限公司 | 一种适用于组网低轨卫星联合定轨的***及方法 |
CN113581501B (zh) * | 2021-08-27 | 2023-02-28 | 重庆两江卫星移动通信有限公司 | 一种适用于组网低轨卫星联合定轨的***及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104309817B (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109782313A (zh) | 地面精密单点定位数据处理方法 | |
CN106556851B (zh) | 一种船载gnss辅助北斗导航卫星定轨方法 | |
CN108548479B (zh) | 基于gnss的桥梁实时监测快速初始化方法 | |
CN110208836B (zh) | 基于卡尔曼滤波的gnss高适应性周跳探测与修复方法 | |
Angrisano et al. | Performance assessment of aided global navigation satellite system for land navigation | |
CN102749637A (zh) | 一种车载gps精确定位的实现方法 | |
US9612340B1 (en) | Systems, methods, devices and subassemblies for creating and delivering crowd-sourced GNSS models | |
CN104309817B (zh) | 基于多台并址接收机的北斗导航卫星区域定轨方法 | |
Seepersad | Reduction of initial convergence period in GPS PPP data processing | |
Ning et al. | Single-frequency precise point positioning enhanced with multi-GNSS observations and global ionosphere maps | |
Tolman et al. | Absolute precise kinematic positioning with GPS and GLONASS | |
Jiang et al. | Precise orbit determination of Haiyang-2D using onboard BDS-3 B1C/B2a observations with ambiguity resolution | |
CN116299596B (zh) | 顾及测站基线长度和对流层约束的海上精密单点定位方法 | |
CN105510946A (zh) | 一种bds卫星载波相位整周模糊度快速解算方法 | |
AU2015227415A1 (en) | Method and system for dynamic-to-dynamic Precise Relative Positioning using Global Navigation Satellite Systems | |
CN105182387A (zh) | 一种北斗三频单历元动动模糊度确定方法 | |
CN105158781A (zh) | 一种地球gnss卫星实现外层空间用户定位的方法 | |
Ferrão | Positioning with combined gps and glonass observations | |
CN115639582B (zh) | 一种基于共视法授时的geo卫星轨道机动期间轨道测定方法 | |
Shin et al. | Inertially aided precise point positioning | |
CN105353385A (zh) | 基于北斗三频点的araim标称偏置估算方法及装置 | |
CN109116396A (zh) | 一种多天线gnss差分定位方法 | |
Kolb et al. | A new method to model the ionosphere across local area networks | |
CN107272024B (zh) | 一种飞行器在高动态环境下的基准站布设方法 | |
Elmunim et al. | Overview of the GPS |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160824 Termination date: 20181011 |
|
CF01 | Termination of patent right due to non-payment of annual fee |