CN101545967A - 卫星导航完好性参数的解算方法和监测*** - Google Patents
卫星导航完好性参数的解算方法和监测*** Download PDFInfo
- Publication number
- CN101545967A CN101545967A CN200910131376A CN200910131376A CN101545967A CN 101545967 A CN101545967 A CN 101545967A CN 200910131376 A CN200910131376 A CN 200910131376A CN 200910131376 A CN200910131376 A CN 200910131376A CN 101545967 A CN101545967 A CN 101545967A
- Authority
- CN
- China
- Prior art keywords
- satellite
- delta
- pseudorange
- monitoring station
- error
- 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
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种卫星导航完好性参数的解算方法和监测***。所述监测***包括一个以上的监测站,监测站得到气象数据和接收相应的GNSS卫星的信号,监测站根据所述GNSS卫星信号和所述气象数据分别得到其能够监测到的GNSS卫星的伪距,选择一个监测站作为主控站,其他监测站将得到的伪距发送给该作为主控站的监测站,该作为主控站的监测站根据其自身得到的伪距和接收到的其他监测站的伪距计算每颗卫星的伪距误差,再根据所述伪距误差计算每颗GNSS卫星的空间信号误差SISE。本发明提供了一种新的计算空间信号误差SISE的方法,能够提高该参数的计算精度和实时性。
Description
一、技术领域
本发明涉及一种卫星导航***信号完好性的解算方法和监测***,特别是一种可在区域范围内监测导航卫星信号的方法和***。
二、背景技术
卫星导航完好性监测技术是用在增强***中,对星座卫星的健康状况进行实时监测,在星座卫星发生故障时,能够将报警信息及时传递给特定用户的一种技术手段。尤其是生命安全相关的用户,如航空、航海,不仅希望得到精确的定位测量数据而且也希望及时获得定位数据的精度及其风险有多大。为满足这些用户的需求,完好性监测技术应运而生。
为了实现对卫星导航星座的完好性的监测,地面***需要计算一系列相关参数,其中空间信号误差SISE(Signal in space error)是一个承上启下的重要参数。有了此参数,才可以进一步实现完好性监测***的其它参数(如SISMA、IF、XPL等)的解算。
参数SISE的解算,涉及轨道估计、钟差估计、电离层误差消除、数学模型建立、解算数据验证等一系列问题,国内还没有机构对此参数全面分析和建模进行***研究,未见相关文献的深入介绍。国外有机构解算过类似参数,但没有公开说明,其所用技术方法对我国实施了保密政策,没有渠道获得技术引进。
对SISE的解算,解决了卫星导航完好性监测技术的一个基础性问题,对我国建设卫星导航增强***具有重大意义。
三、发明内容
本发明所要解决的技术问题是:高精度实时估计空间信号误差(SISE)。
针对现存的问题,我们的方法分离多种误差因素,以最小二乘为基础方法,实时计算SISE,使SISE实时反映空间信号的综合误差。
一个以上的监测站得到气象数据和接收相应的GNSS卫星的信号,监测站根据所述GNSS卫星信号和所述气象数据得到其能够监测到的GNSS卫星的伪距,选择一个监测站作为主控站,其他监测站将得到的伪距发送给该作为主控站的监测站,该作为主控站的监测站根据其自身得到的伪距和接收到的其他监测站的伪距计算每颗卫星的伪距误差,再根据所述伪距误差计算每颗GNSS卫星的空间信号误差SISE。
四、附图说明
图1为地面监测***结构图;
图2为监测站计算伪距的方法流程图;
图3为作为主控站的监测站计算空间信号误差SISE的方法流程图。
五、具体实施方式
监测站计算得到每颗卫星的伪距的处理步骤如下。
第一步:监测站接收相应的GNSS卫星(可多颗)的信号,解码得到相应的原始观测量,包括伪距,相位,星历。主控站和监测站的气象监测设备记录监测站附近的气象数据,包括大气压强,湿度和温度。
GNSS***中伪距和载波相位观测量的模型:
GNSS的L1频率上的伪距和相位分别是ρf1和φf1,其模型是,
ρf1=|rs-rr|+T+I1+R+c(δr-δs)+ε1 (1.1)
φf1=|rs-rr|+T-I1+R+c(δr-δs)+λ1N1+υ1 (1.2)
GNSS的L2频率上的伪距和相位分别是ρf2和φf2,其模型是,
ρf2=|rs-rr|+T+I2+R+c(δr-δs)+ε2 (1.3)
φf2=|rs-rr|+T-I2+R+c(δr-δs)+λ2N2+υ2 (1.4)
其中,rs是卫星真实位置矢量,rr是接收机真实位置矢量,矢量的起点为地心。T是对流层误差,I1和I2分别是L1和L2频率上的电离层误差,R是相对论误差。c为光速,δr是接收机钟差,δs是卫星钟差。ε1,ε2分别是伪距噪声。υ1,υ2分别是相位噪声。
第二步:监测站进行相位平滑伪距,得到平滑后的伪距。载波相位平滑公式:
其中,
将公式(1.1)-(1.4)代入上式,简化为,
第三步:监测站由经过平滑后的双频伪距组合成消电离层延迟伪距。由两个频率上平滑伪距组成消电离层伪距ρIFree:
第四步:监测站计算相对论改正误差。
相对论R的计算公式:
其中,G为物理常数,a是卫星轨道的长半轴,e是卫星轨道的偏心率,E为卫星轨道的偏近点角。
第五步:监测站计算对流层改正误差。
对流层改正采用Hopfield模型,计算公式:
其中,el为卫星仰角,Ts是温度,p是气压,hr是相对湿度,h是接收机高度。
第六步:监测站计算地球自转改正误差。
地球自转改正计算公式:
其中,ω为地球自转角速度;c为光速;xsi和xri分别代表卫星位置矢量和测站位置矢量的分量,i=1,2对应x,y分量。
第七步:监测站计算卫星钟差。
改正公式为:
δs=a0+a1(t-t0)+a2(t-t0)2
其中,a0,a1,a2分别是钟差,频差和频漂系数,t0为参考时刻,t为当前时刻。
第八步:监测站计算消除了相对论改正,对流层改正,电离层改正,地球自转改正,卫星钟差改正后的伪距。经过对流层、相对论改正、地球自转改正和卫星钟差改正之后的消电离层伪距为,
ρIFree=|rs-rr|+cδr-cΔδs+υ
Δδs为卫星钟差的残差。
作为主控站的监测站根据其自身的伪距和接收到的其他监测站的伪距计算伪距误差,并根据伪距误差计算每颗GNSS卫星的空间信号误差SISE的步骤如下:
第一步:选择一个监测站作为主控站,主控站接收其他监测站传输的伪距。
第二步:主控站由其他监测站伪距和其自身的伪距计算其他监测站接收机相对主控站接收机相对钟差。计算公式:
设监测站j和主控站M,均监测到同一导航卫星i,经过前述计算后,得到的消电离层伪距站间差分为
变形为,
假设同一时刻监测站j和基准站M共视Kj颗卫星,则接收机相对钟差的估计为,
第三步:将接收机相对钟差代入监测站的伪距,消去接收机相对钟差,同步接收机钟差。
第四步:主控站计算主控站钟的钟差,并将该钟差代入主控站和其他监测站伪距,消去主控站钟差。
主钟钟差计算方法如下:
设主控站同时观测到N颗(N>=4)导航卫星,
经过平滑伪距,组成消电离层组合,对流层改正,相对论改正,地球自转改正,卫星钟差改正后,伪距分别为ρ1,ρ2,...,ρN,形成如下的定位方程:
在地心地固坐标系下,几何矩阵
则最小二乘解为
将主控站钟差代入主控站伪距和其他监测站伪距,消去主控站钟差,方法如下:
第五步:根据所有在同一时刻监测到同一颗卫星的监测站伪距计算伪距误差方法如下。
主控站由主控站伪距误差和监测站伪距误差计算空间信号误差SISE,可以采用不同的数学模型,得到不同表达形式的SISE。第一种模型采用标量形式,公式为: 或
第二种模型采用矢量模型计算SISE,公式如下:
在卫星星固坐标系下,数学模型为
设
Claims (10)
1.一种卫星导航完好性参数的解算方法,其特征在于:
A.一个以上的监测站接收GNSS卫星信号,并获得其附近的气象数据,监测站根据所述GNSS卫星信号和所述气象数据得到其能够监测到的GNSS卫星的伪距,选择一个监测站作为主控站,其他监测站将其得到的GNSS卫星的伪距发送给该作为主控站的监测站;
B.该作为主控站的监测站根据其自身得到的伪距和接收到的其他监测站的伪距计算每颗GNSS卫星的伪距误差,再根据所述伪距误差计算每颗GNSS卫星的空间信号误差SISE。
4.如权利要求1所述的方法,其特征在于:所述作为主控站的监测站根据伪距误差计算每颗GNSS卫星的空间信号误差SISE的过程具体为:
B1:作为主控站的监测站根据其他监测站发来的伪距和其自身的伪距计算其他监测站接收机与该作为主控站的监测站的接收机的相对钟差;
B2:将所述接收机的相对钟差分别代入相应的其他监测站的伪距以消去相对钟差;
B3:作为主控站的监测站计算其自身接收机的钟差,并将该钟差代入其自身的伪距,消去钟差以计算作为主控站的监测站的伪距误差;
B4:将作为主控站的监测站的接收机钟差代入由B2得到的其他监测站伪距,消去作为主控站的监测站的接收机钟差,以计算其他监测站的伪距误差;
B5:根据所有在同一时刻监测到同一颗卫星的所有监测站伪距误差,计算空间信号误差SISE。
5.如权利要求4所述的方法,监测站根据所述GNSS卫星信号分别得到双频伪距、双频相位和星历,监测站获得的气象数据包括大气压强、湿度和温度,其特征在于:监测站计算得到每颗其所能监测到的GNSS卫星的伪距具体为:
A1:对每颗GNSS卫星的双频相位进行实时的周跳探测;
A2:用步骤A1得到的相位对对应卫星的双频伪距进行平滑;
A3:由经过平滑后的双频伪距组合成消电离层延迟伪距;
A4:计算相对论改正误差、对流层改正误差、地球自转改正误差和卫星钟差;
A5:应用步骤A4中的相对论改正,对流层改正,电离层改正,地球自转改正和卫星钟差修正步骤A3得到的伪距。
6.一种卫星导航完好性参数的监测***,其特征在于:
所述监测***包括一个以上的监测站,监测站得到气象数据和接收相应的GNSS卫星的信号,监测站根据所述GNSS卫星信号和所述气象数据得到其能够监测到的GNSS卫星的伪距,选择一个监测站作为主控站,其他监测站将得到的伪距发送给该作为主控站的监测站,该作为主控站的监测站根据其自身得到的伪距和接收到的其他监测站的伪距计算每颗卫星的伪距误差,再根据所述伪距误差计算每颗GNSS卫星的空间信号误差SISE。
7.如权利要求6所述的***,其特征在于:作为主控站的监测站根据伪距误差计算每颗GNSS卫星的空间信号误差SISE具体为:
作为主控站的监测站根据其他监测站发来的伪距和其自身的伪距计算其他监测站接收机与该作为主控站的监测站的接收机的相对钟差;
将所述接收机相对钟差分别代入相应的其他监测站的伪距,以消去相对钟差;
作为主控站的监测站计算其自身接收机的钟差,并将该钟差代入其自身的伪距,消去钟差以计算作为主控站的监测站的伪距误差;
将作为主控站的监测站的接收机钟差分别代入由上一步得到的其他监测站伪距,消去作为主控站的监测站的接收机钟差以计算监测站伪距误差;
根据所有在同一时刻监测到同一颗卫星的所有监测站伪距误差,计算空间信号误差SISE。
8.如权利要求7所述的***,监测站根据所述GNSS卫星信号得到双频伪距、双频相位和星历,监测站获得气象数据,该气象数据包括大气压强、湿度和温度,其特征在于:监测站计算得到每颗其所能监测到GNSS卫星的伪距具体为:
对每颗GNSS卫星的双频相位进行实时的周跳探测;
用上一步骤得到的相位对对应卫星的双频伪距进行平滑;
由经过平滑后的双频伪距组合成消电离层延迟的伪距;
计算相对论改正误差、对流层改正误差、地球自转改正误差和卫星钟差;
应用上一步骤中的相对论改正,对流层改正,电离层改正,地球自转改正和卫星钟差修正所述消电离层延迟的伪距。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101313765A CN101545967B (zh) | 2009-04-16 | 2009-04-16 | 卫星导航完好性参数的解算方法和监测*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101313765A CN101545967B (zh) | 2009-04-16 | 2009-04-16 | 卫星导航完好性参数的解算方法和监测*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101545967A true CN101545967A (zh) | 2009-09-30 |
CN101545967B CN101545967B (zh) | 2011-11-16 |
Family
ID=41193207
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009101313765A Expired - Fee Related CN101545967B (zh) | 2009-04-16 | 2009-04-16 | 卫星导航完好性参数的解算方法和监测*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101545967B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102004259A (zh) * | 2010-09-17 | 2011-04-06 | 浙江大学 | 高灵敏度环境下基于多普勒平滑伪距的卫星导航解算定位方法 |
CN102788985A (zh) * | 2012-08-31 | 2012-11-21 | 中国科学院国家授时中心 | 广域加权伪距差分对定位误差的改正方法 |
CN103941266A (zh) * | 2014-03-30 | 2014-07-23 | 桂林电子科技大学 | 一种基于互差的卫星完好性监测方法及*** |
CN104297557A (zh) * | 2014-10-08 | 2015-01-21 | 北京航空航天大学 | 一种适用于多飞行器自由飞行的联合导航自主完好性监测方法 |
CN104330805A (zh) * | 2014-10-24 | 2015-02-04 | 北京航空航天大学 | 卫星导航增强***中随机误差检测方法 |
CN104483678A (zh) * | 2014-12-04 | 2015-04-01 | 北京航空航天大学 | 一种空地协同的多星座卫星导航完好性多级监测方法 |
CN104330805B (zh) * | 2014-10-24 | 2017-01-04 | 北京航空航天大学 | 卫星导航增强***中随机误差检测方法 |
CN109100748A (zh) * | 2018-08-14 | 2018-12-28 | 西安空间无线电技术研究所 | 一种基于低轨星座的导航完好性监测***及方法 |
CN109901204A (zh) * | 2019-03-27 | 2019-06-18 | 北京航空航天大学 | 一种基于伪距误差分布模型的gbas完好性性能评估方法 |
CN110361755A (zh) * | 2019-03-12 | 2019-10-22 | 中国矿业大学 | 一种基于oedop因子的多卫星导航***监测站优化选取方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1637899A1 (en) * | 2004-09-20 | 2006-03-22 | EADS Astrium GmbH | Method and apparatus for providing integrity information for users of global navigation system |
DE102006046001A1 (de) * | 2006-09-27 | 2008-04-03 | Eads Astrium Gmbh | Verfahren zum Einstellen der Integritätsalarmschranke in einem Satellitennavigationssystem |
CN100523858C (zh) * | 2007-11-22 | 2009-08-05 | 北京航空航天大学 | 一种应用于卫星导航***的定位及监测方法 |
-
2009
- 2009-04-16 CN CN2009101313765A patent/CN101545967B/zh not_active Expired - Fee Related
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102004259A (zh) * | 2010-09-17 | 2011-04-06 | 浙江大学 | 高灵敏度环境下基于多普勒平滑伪距的卫星导航解算定位方法 |
CN102788985A (zh) * | 2012-08-31 | 2012-11-21 | 中国科学院国家授时中心 | 广域加权伪距差分对定位误差的改正方法 |
CN103941266A (zh) * | 2014-03-30 | 2014-07-23 | 桂林电子科技大学 | 一种基于互差的卫星完好性监测方法及*** |
CN104297557A (zh) * | 2014-10-08 | 2015-01-21 | 北京航空航天大学 | 一种适用于多飞行器自由飞行的联合导航自主完好性监测方法 |
CN104297557B (zh) * | 2014-10-08 | 2017-04-26 | 北京航空航天大学 | 一种适用于多飞行器自由飞行的联合导航自主完好性监测方法 |
CN104330805B (zh) * | 2014-10-24 | 2017-01-04 | 北京航空航天大学 | 卫星导航增强***中随机误差检测方法 |
CN104330805A (zh) * | 2014-10-24 | 2015-02-04 | 北京航空航天大学 | 卫星导航增强***中随机误差检测方法 |
CN104483678A (zh) * | 2014-12-04 | 2015-04-01 | 北京航空航天大学 | 一种空地协同的多星座卫星导航完好性多级监测方法 |
CN104483678B (zh) * | 2014-12-04 | 2017-03-01 | 北京航空航天大学 | 一种空地协同的多星座卫星导航完好性多级监测方法 |
CN109100748A (zh) * | 2018-08-14 | 2018-12-28 | 西安空间无线电技术研究所 | 一种基于低轨星座的导航完好性监测***及方法 |
CN110361755A (zh) * | 2019-03-12 | 2019-10-22 | 中国矿业大学 | 一种基于oedop因子的多卫星导航***监测站优化选取方法 |
CN110361755B (zh) * | 2019-03-12 | 2023-04-07 | 中国矿业大学 | 一种基于oedop因子的多卫星导航***监测站优化选取方法 |
CN109901204A (zh) * | 2019-03-27 | 2019-06-18 | 北京航空航天大学 | 一种基于伪距误差分布模型的gbas完好性性能评估方法 |
CN109901204B (zh) * | 2019-03-27 | 2020-12-04 | 北京航空航天大学 | 一种基于伪距误差分布模型的gbas完好性性能评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101545967B (zh) | 2011-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110031881B (zh) | 高精度星间激光测距辅助精密单点定位的方法 | |
CN109061696B (zh) | 一种确定导航卫星轨道和钟差的方法 | |
CN101545967A (zh) | 卫星导航完好性参数的解算方法和监测*** | |
CN104483678B (zh) | 一种空地协同的多星座卫星导航完好性多级监测方法 | |
CN106468774B (zh) | 一种应用于星基增强***的星历星钟改正参数及空间信号完好性参数方法 | |
CN101950025B (zh) | 用于局域增强***的数据质量监测方法 | |
CN104865587B (zh) | 基于geo星座粗定位的北斗接收机快速定位方法 | |
CN105044747B (zh) | 一种基于多星共视和滤波的时间同步装置及其方法 | |
CN104714244A (zh) | 一种基于抗差自适应Kalman滤波的多***动态PPP解算方法 | |
CN101419275B (zh) | 基于多接收机的局域机场监视方法和*** | |
WO2016034252A1 (en) | Ionospheric scintillation prediction | |
CN103760571A (zh) | 用于gps的基于影响因素特征的脆弱性监测***及方法 | |
CN102608633B (zh) | 一种卫星定位伪距差分方法 | |
CN110007326B (zh) | 一种用于星基增强***的双频测距误差参数生成方法 | |
CN109613579B (zh) | 一种基于最小二乘算法计算整周模糊度的方法和*** | |
CN103529482B (zh) | 一种高精度确定载体动态加速度的方法 | |
Wielgosz et al. | Troposphere modeling for precise GPS rapid static positioning in mountainous areas | |
US20130090858A1 (en) | System for measuring coseismic movements or vibrations of structures based on global navigation satellite systems-gnss and/or pseudolites | |
CN101893712B (zh) | 用于地球静止卫星精密定轨的选权拟合方法 | |
CN111596315A (zh) | 一种用于实时监测双频多星座星基增强***性能的方法 | |
CN105044741A (zh) | 一种伪距相位综合广域差分改正值的求解方法 | |
CN104316943A (zh) | 一种伪距离和多普勒组合差分定位***及方法 | |
CN113253303A (zh) | 一种用于实时监测单频星基增强***性能的方法 | |
Gunning et al. | Multi-GNSS constellation anomaly detection and performance monitoring | |
CN111505677A (zh) | 一种基于地面参考站观测的geo卫星轨道机动修复方法 |
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: 20111116 Termination date: 20180416 |
|
CF01 | Termination of patent right due to non-payment of annual fee |