CN105277958A - 一种采用三频数据进行相对定位解算的方法 - Google Patents
一种采用三频数据进行相对定位解算的方法 Download PDFInfo
- Publication number
- CN105277958A CN105277958A CN201410339029.2A CN201410339029A CN105277958A CN 105277958 A CN105277958 A CN 105277958A CN 201410339029 A CN201410339029 A CN 201410339029A CN 105277958 A CN105277958 A CN 105277958A
- Authority
- CN
- China
- Prior art keywords
- observation
- blur level
- wide item
- ultra
- frequency
- 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
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种采用三频数据进行相对定位解算的方法,所述方法包括以下步骤:获得站点已知信息、基准站三频观测数据和流动站三频观测数据并计算卫星位置;利用所述站点已知信息和所述基准站、流动站的三频观测数据进行组合构建超宽项组合观测值与宽项组合观测值后通过若干历元求解超宽项模糊度和宽项模糊度的值;通过固定了模糊度的超宽项组合观测值和宽项组合观测值与Geometric-free观测值进行组合求解出所有频率的模糊度;根据各固定模糊度的多频观测值组建求解弱电离层观测值模糊度。本发明提供的采用三频数据进行相对定位解算的方法,通过采用三频数据进行相对定位解算,模糊度固定速度快,相对定位基线更长。
Description
技术领域
本发明涉及卫星导航技术领域,尤其涉及一种采用三频数据进行相对定位解算的方法。
背景技术
在BDS(BeiDouNavigationSatelliteSystem,北斗卫星导航系统)投入使用之前,现有的导航系统采用的是基于双频数据进行结算,采用这种方法进行解算一般耗时较长,且模糊度固定结果的成功率容易受伪距精度影响,并且基线长度越长,固定的成功率越低,进而影响相对定位结算结果的准确性。
随着支持三频数据的BDS系统正式投入使用和GPS系统现代化增加L5波段,使用三频数据进行结算成为可能。
发明内容
鉴于目前卫星导航技术领域存在的上述不足,本发明提供一种采用三频数据进行相对定位解算的方法,模糊度固定速度快、相对定位基线更长。
为达到上述目的,本发明的实施例采用如下技术方案:
一种采用三频数据进行相对定位解算的方法,所述采用三频数据进行相对定位解算的方法包括以下步骤:
获得站点已知信息、基准站三频观测数据和流动站三频观测数据并计算卫星位置;
利用所述站点已知信息和所述基准站、流动站的三频观测数据进行组合构建超宽项组合观测值与宽项组合观测值后通过若干历元求解超宽项模糊度和宽项模糊度的值;
通过固定了模糊度的超宽项组合观测值和宽项组合观测值与Geometric-free观测值进行组合求解出所有频率的模糊度;
根据各固定模糊度的多频观测值组建求解弱电离层观测值模糊度。
依照本发明的一个方面,所述步骤利用三频观测数据进行组合构建超宽项组合观测值与宽项组合观测值后通过若干历元求解超宽项模糊度和宽项模糊度的值具体包括以下步骤:
假设三个频率相位观测值对应波长分别为λ1、λ2、λ5,频率分别为f1、f2、f5,构建双差观测方程如下:
其中:
Δ▽N(i,j,k)=i*Δ▽N1+j*Δ▽N2+k*Δ▽N5
Δ▽为双差算子,为载波相位观测值,δtr为对流层延迟误差,δion为L1相位观测值的电离层延迟误差,N为模糊度,ξ为噪声。
构建超宽项观测值和宽项观测值然后代入所述的双差观测方程中,计算得到对应的超宽项模糊度Δ▽N(0,1,-1)和宽项模糊度Δ▽N(1,-6,5)。
依照本发明的一个方面,所述步骤利用三频观测数据进行组合构建超宽项组合观测值与宽项组合观测值后通过若干历元求解超宽项模糊度和宽项模糊度的值实施之前需执行步骤:选取参考卫星。
依照本发明的一个方面,所述步骤利用三频观测数据进行组合构建超宽项组合观测值与宽项组合观测值后通过若干历元求解超宽项模糊度和宽项模糊度的值还包括以下步骤:
进行粗差探测;
进行周跳探测与修复;
删除短弧段历元观测值。
依照本发明的一个方面,所述步骤通过固定了模糊度的超宽项组合观测值和宽项组合观测值与Geometric-free观测值进行组合求解出所有频率的模糊度包括以下步骤:
根据下列方程:
结合超宽项模糊度Δ▽N(0,1,-1)的值和宽项模糊度Δ▽N(1,-6,5)的值,求出其中一频率的模糊度Δ▽N(1,0,0);
再求出另两频率的模糊度Δ▽N(0,1,0)和Δ▽N(0,0,1)。
依照本发明的一个方面,所述步骤根据各固定模糊度的多频观测值组建求解弱电离层观测值模糊度包括以下步骤:
组建L1、L2的弱电离层观测值的公式如下:
同理可得其他弱电离层观测值。
本发明实施的优点:本方法所述的采用三频数据进行相对定位解算的方法通过采集三频数据,然后利用三频观测数据进行组合构建超宽项组合观测值与宽项组合观测值后通过若干历元求解超宽项模糊度和宽项模糊度的值,然后通过固定了模糊度的超宽项组合观测值和宽项组合观测值与Geometric-free观测值进行组合求解出所有频率的模糊度,最后根据各固定模糊度的多频观测值组建求解弱电离层观测值模糊度,模糊度固定速度快,相对定位基线更长,可实现中长基线模糊度快速固定,最终实现相对定位解算。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明所述的一种采用三频数据进行相对定位解算的方法示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,一种采用三频数据进行相对定位解算的方法,所述采用三频数据进行相对定位解算的方法包括以下步骤:
步骤S1:获得站点已知信息、基准站三频观测数据和流动站三频观测数据并计算卫星位置;
所述步骤S1获得站点已知信息、基准站三频观测数据和流动站三频观测数据并计算卫星位置的具体实施方式可为:通过BDS(BeiDouNavigationSatelliteSystem,北斗卫星导航系统)获得站点已知信息、基准站三频观测数据和流动站三频观测数据,所述站点已知信息包括BDS广播星历,然后通过现有的数据计算卫星的位置。
步骤S2:利用所述站点已知信息和所述基准站、流动站的三频观测数据进行组合构建超宽项组合观测值与宽项组合观测值后通过若干历元求解超宽项模糊度和宽项模糊度的值;
在所述步骤S2利用所述站点已知信息和所述基准站、流动站的三频观测数据进行组合构建超宽项组合观测值与宽项组合观测值后通过若干历元求解超宽项模糊度和宽项模糊度的值执行前需选取参考卫星,以方便进行定位解算。
所述步骤S2利用所述站点已知信息和所述基准站、流动站的三频观测数据进行组合构建超宽项组合观测值与宽项组合观测值后通过若干历元求解超宽项模糊度和宽项模糊度的值的具体实施步骤如下:
假设三个频率相位观测值对应波长分别为λ1、λ2、λ5,频率分别为f1、f2、f5,构建双差观测方程如下:
其中:
Δ▽N(i,j,k)=i*Δ▽N1+j*Δ▽N2+k*Δ▽N5
Δ▽为双差算子,为载波相位观测值,δtr为对流层延迟误差,δion为L1相位观测值的电离层延迟误差,N为模糊度,ξ为噪声。
构建超宽项观测值和宽项观测值然后代入所述的双差观测方程中,计算得到对应的超宽项模糊度Δ▽N(0,1,-1)和宽项模糊度Δ▽N(1,-6,5)。
其中,还需执行以下步骤:进行粗差探测,进行周跳探测与修复和删除短弧段历元观测值。
步骤S3:通过固定了模糊度的超宽项组合观测值和宽项组合观测值与Geometric-free观测值进行组合求解出所有频率的模糊度;
所述步骤S3通过固定了模糊度的超宽项组合观测值和宽项组合观测值与Geometric-free观测值进行组合求解出所有频率的模糊度的具体实施步骤如下:
根据下列方程:
结合超宽项模糊度Δ▽N(0,1,-1)的值和宽项模糊度Δ▽N(1,-6,5)的值,求出其中一频率的模糊度Δ▽N(1,0,0);
再求出另两频率的模糊度Δ▽N(0,1,0)和Δ▽N(0,0,1)。
步骤S4:根据各固定模糊度的多频观测值组建求解弱电离层观测值模糊度;
所述步骤S4根据各固定模糊度的多频观测值组建求解弱电离层观测值模糊度的具体实施步骤可为:
组建L1、L2的弱电离层观测值的公式如下:
同理可得其他弱电离层观测值。
本发明实施的优点:本方法所述的采用三频数据进行相对定位解算的方法通过采集三频数据,然后利用三频观测数据进行组合构建超宽项组合观测值与宽项组合观测值后通过若干历元求解超宽项模糊度和宽项模糊度的值,然后通过固定了模糊度的超宽项组合观测值和宽项组合观测值与Geometric-free观测值进行组合求解出所有频率的模糊度,最后根据各固定模糊度的多频观测值组建求解弱电离层观测值模糊度,模糊度固定速度快,相对定位基线更长,可实现中长基线模糊度快速固定,最终实现相对定位解算。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本领域技术的技术人员在本发明公开的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (6)
1.一种采用三频数据进行相对定位解算的方法,其特征在于,所述采用三频数据进行相对定位解算的方法包括以下步骤:
获得站点已知信息、基准站三频观测数据和流动站三频观测数据并计算卫星位置;
利用所述站点已知信息和所述基准站、流动站的三频观测数据进行组合构建超宽项组合观测值与宽项组合观测值后通过若干历元求解超宽项模糊度和宽项模糊度的值;
通过固定了模糊度的超宽项组合观测值和宽项组合观测值与Geometric-free观测值进行组合求解出所有频率的模糊度;
根据各固定模糊度的多频观测值组建求解弱电离层观测值模糊度。
2.根据权利要求1所述的采用三频数据进行相对定位解算的方法,其特征在于,所述步骤利用三频观测数据进行组合构建超宽项组合观测值与宽项组合观测值后通过若干历元求解超宽项模糊度和宽项模糊度的值具体包括以下步骤:
假设三个频率相位观测值对应波长分别为λ1、λ2、λ5,频率分别为f1、f2、f5,构建双差观测方程如下:
其中:
Δ▽N(i,j,k)=i*Δ▽N1+j*Δ▽N2+k*Δ▽N5
Δ▽为双差算子,为载波相位观测值,δtr为对流层延迟误差,δion为L1相位观测值的电离层延迟误差,N为模糊度,ξ为噪声。
构建超宽项观测值和宽项观测值然后代入所述的双差观测方程中,计算得到对应的超宽项模糊度Δ▽N(0,1,-1)和宽项模糊度Δ▽N(1,-6,5)。
3.根据权利要求2所述的采用三频数据进行相对定位解算的方法,其特征在于,所述步骤利用三频观测数据进行组合构建超宽项组合观测值与宽项组合观测值后通过若干历元求解超宽项模糊度和宽项模糊度的值实施之前需执行步骤:选取参考卫星。
4.根据权利要求3所述的采用三频数据进行相对定位解算的方法,其特征在于,所述步骤利用三频观测数据进行组合构建超宽项组合观测值与宽项组合观测值后通过若干历元求解超宽项模糊度和宽项模糊度的值还包括以下步骤:
进行粗差探测;
进行周跳探测与修复;
删除短弧段历元观测值。
5.根据权利要求2至4之一所述的采用三频数据进行相对定位解算的方法,其特征在于,所述步骤通过固定了模糊度的超宽项组合观测值和宽项组合观测值与Geometric-free观测值进行组合求解出所有频率的模糊度包括以下步骤:
根据下列方程:
结合超宽项模糊度Δ▽N(0,1,-1)的值和宽项模糊度Δ▽N(1,-6,5)的值,求出其中一频率的模糊度Δ▽N(1,0,0);
再求出另两频率的模糊度Δ▽N(0,1,0)和Δ▽N(0,0,1)。
6.根据权利要求5所述的采用三频数据进行相对定位解算的方法,其特征在于,所述步骤根据各固定模糊度的多频观测值组建求解弱电离层观测值模糊度包括以下步骤:
组建L1、L2的弱电离层观测值的公式如下:
同理可得其他弱电离层观测值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410339029.2A CN105277958A (zh) | 2014-07-16 | 2014-07-16 | 一种采用三频数据进行相对定位解算的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410339029.2A CN105277958A (zh) | 2014-07-16 | 2014-07-16 | 一种采用三频数据进行相对定位解算的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105277958A true CN105277958A (zh) | 2016-01-27 |
Family
ID=55147299
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410339029.2A Pending CN105277958A (zh) | 2014-07-16 | 2014-07-16 | 一种采用三频数据进行相对定位解算的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105277958A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106932788A (zh) * | 2016-12-30 | 2017-07-07 | 广州市中海达测绘仪器有限公司 | 一种gnss参考站网三频解算方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006108227A1 (en) * | 2005-04-11 | 2006-10-19 | Yanming Feng | Improved phase ambiguity resolution using three gnss signals |
US20090262013A1 (en) * | 2003-10-28 | 2009-10-22 | Trimble Navigation Limited, A California Corporation | Ambiguity estimation of GNSS signals for three or more carriers |
CN103605145A (zh) * | 2013-12-04 | 2014-02-26 | 上海华测导航技术有限公司 | 基于gnss多频数据和cors实现网络实时动态定位的方法 |
CN103675874A (zh) * | 2013-12-20 | 2014-03-26 | 北京遥测技术研究所 | 一种北斗导航系统三频载波相位整周模糊度确定方法 |
CN103728643A (zh) * | 2014-01-20 | 2014-04-16 | 东南大学 | 附有宽巷约束的北斗三频网络rtk模糊度单历元固定方法 |
-
2014
- 2014-07-16 CN CN201410339029.2A patent/CN105277958A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090262013A1 (en) * | 2003-10-28 | 2009-10-22 | Trimble Navigation Limited, A California Corporation | Ambiguity estimation of GNSS signals for three or more carriers |
WO2006108227A1 (en) * | 2005-04-11 | 2006-10-19 | Yanming Feng | Improved phase ambiguity resolution using three gnss signals |
CN103605145A (zh) * | 2013-12-04 | 2014-02-26 | 上海华测导航技术有限公司 | 基于gnss多频数据和cors实现网络实时动态定位的方法 |
CN103675874A (zh) * | 2013-12-20 | 2014-03-26 | 北京遥测技术研究所 | 一种北斗导航系统三频载波相位整周模糊度确定方法 |
CN103728643A (zh) * | 2014-01-20 | 2014-04-16 | 东南大学 | 附有宽巷约束的北斗三频网络rtk模糊度单历元固定方法 |
Non-Patent Citations (5)
Title |
---|
余学祥: "《GPS测量与数据处理》", 31 October 2013, 中国矿业大学出版社 * |
李博峰等: ""中长基线三频GNSS模糊度的快速算法"", 《测绘学报》 * |
李博峰等: ""利用三频GNSS进行长距离实时精密导航"", 《武汉大学学报信息科学版》 * |
李博峰等: ""基于北斗卫星系统的长距离精密RTK定位"", 《中国卫星导航学术年会》 * |
李金龙: ""GNSS三频精密定位数据处理方法研究"", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106932788A (zh) * | 2016-12-30 | 2017-07-07 | 广州市中海达测绘仪器有限公司 | 一种gnss参考站网三频解算方法及系统 |
CN106932788B (zh) * | 2016-12-30 | 2019-10-25 | 广州市中海达测绘仪器有限公司 | 一种gnss参考站网三频解算方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhao et al. | Real-time detection and repair of cycle slips in triple-frequency GNSS measurements | |
CN110007320B (zh) | 网络rtk解算方法 | |
Chen et al. | A double-differenced cycle slip detection and repair method for GNSS CORS network | |
CN105738926A (zh) | 一种glonass系统接收机间相位频间偏差标定方法 | |
CN103605145B (zh) | 基于gnss多频数据和cors实现网络实时动态定位的方法 | |
DE102013206544A1 (de) | Fortschrittliche Positionierung mit globalen Navigationssatellitensystemen (GNSS) unter Verwendung von präziser Satelliteninformation | |
CN104102822A (zh) | 一种多频gnss观测值随机特性建模方法 | |
Paziewski et al. | Selected properties of GPS and Galileo-IOV receiver intersystem biases in multi-GNSS data processing | |
CN105158783A (zh) | 一种实时动态差分定位方法及其设备 | |
CN104483690A (zh) | 一种gnss三频精密单点定位模糊度固定方法 | |
CN105158782A (zh) | 一种bds和gps观测信息融合的宽巷模糊度解算方法 | |
CN106896386A (zh) | Glonass频间偏差精确估计方法 | |
CN111913201B (zh) | 一种gnss差分定位方法、装置及计算机可读存储介质 | |
CN102650694A (zh) | 基于北斗四频信号的中长基线模糊度解算方法 | |
DE102013003265A1 (de) | Verfahren und vorrichtung zur synchronisation und genauen zeitberechnung für gnss-empfänger | |
CN104483689A (zh) | 一种bds参考站三频载波相位整周模糊度确定方法 | |
CN104180781A (zh) | 一种单双频gps混合网变形监测数据处理方法 | |
CN114859387A (zh) | 静态站周跳探测的方法、装置、设备及计算机存储介质 | |
CN103760586B (zh) | 一种在gps姿态测量中快速探测与修复周跳的方法 | |
CN104007447A (zh) | 北斗、gps及glonass基线向量联合解算的方法及系统 | |
CN105277958A (zh) | 一种采用三频数据进行相对定位解算的方法 | |
CN105445771A (zh) | 一种单频rtk融合测试分析方法 | |
JP2000193733A (ja) | 測位装置及び測位方法 | |
EP3620825B1 (en) | Precise point positioning method and positioning apparatus and recording medium thereof | |
RU2592077C1 (ru) | Способ измерения взаимного положения объектов по сигналам глобальных навигационных спутниковых систем |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160127 |
|
RJ01 | Rejection of invention patent application after publication |