CN101487883A - 一种适用于多模卫星导航系统的电离层格网校正方法 - Google Patents
一种适用于多模卫星导航系统的电离层格网校正方法 Download PDFInfo
- Publication number
- CN101487883A CN101487883A CNA2009100793642A CN200910079364A CN101487883A CN 101487883 A CN101487883 A CN 101487883A CN A2009100793642 A CNA2009100793642 A CN A2009100793642A CN 200910079364 A CN200910079364 A CN 200910079364A CN 101487883 A CN101487883 A CN 101487883A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- ionospheric
- user
- msubsup
- 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
- 238000000034 method Methods 0.000 title claims abstract description 28
- 239000005433 ionosphere Substances 0.000 claims abstract description 34
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 230000035515 penetration Effects 0.000 claims description 56
- 238000010606 normalization Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims 1
- 230000008030 elimination Effects 0.000 abstract 1
- 238000003379 elimination reaction Methods 0.000 abstract 1
- 238000005303 weighing Methods 0.000 abstract 1
- 230000009977 dual effect Effects 0.000 description 7
- 238000012544 monitoring process Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- NUHSROFQTUXZQQ-UHFFFAOYSA-N isopentenyl diphosphate Chemical compound CC(=C)CCO[P@](O)(=O)OP(O)(O)=O NUHSROFQTUXZQQ-UHFFFAOYSA-N 0.000 description 2
- 230000003416 augmentation Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种应用于多模卫星导航系统的电离层格网校正方法,根据来自不同卫星导航系统的数据,将电离层时延改正值及其方差归一化到指定的频率上,根据归一化的电离层穿透点的时延值及其方差分别进行加权计算,确定指定频率下格网点的电离层天顶时延改正值;利用格网点电离层天顶时延改正值进行内插,确定用户在归一化频率下所接收的卫星测距信号天顶电离层时延改正值,最终用户确定接收不同频率导航测距信号所产生的电离层时延改正值,将用户接收到的导航测距信号产生的电离层时延改正值减去用户电离层时延值,得到经消除后的用户电离层时延值。本发明减小电离层时延改正误差,提高电离层格网校正精度,应用该用户伪距值进行定位计算时可以大大提高用户定位精度。
Description
技术领域
本发明属于卫星导航领域,具体地说,是指一种适用于多模卫星导航系统的电离层格网校正方法。
背景技术
卫星导航定位系统的主要误差源之一是电波传播的电离层延迟。广大C/A码接收机用户使用卫星广播的Kloubuchar模型参数法只能修正50%左右的电离层误差。
美国联邦航空局(FAA)提出的广域增强系统(WAAS),采用格网校正法,能够为单频用户提供实时、大范围精确的电离层修正。这种校正法将区域内电离层参考面上(通常为350~400km)经度线和纬度线按5度的间隔分割成格网。根据分布于WAAS系统区域内的各监测站的实时精确测定的电离层时延值,用一定算法获得格网点的天顶电离层时延值,然后通过静止卫星在区域内广播。用户根据自己接收的信号和格网点的数据,对测距信号进行电离层时延修正。
在这种单卫星导航系统中,一些格网点因周围缺乏足够的穿透点分布而不可用,使得电离层格网改正方法的服务范围受限。
发明内容
本发明为解决单个卫星导航系统存在时,电离层格网校正方法因穿透点数目不足而带来的局限性,提出了一种适用于多模卫星导航系统的电离层格网校正方法。
本发明一种应用于多模卫星导航系统的电离层格网校正方法,通过五个步骤来实现:步骤一:首先根据来自不同卫星导航系统的数据,将电离层时延改正值及其方差归一化到用户任意指定的频率f上;设第m个导航系统穿透点电离层时延改正值为和方差为相应频率为fm,将和归一化到指定频率f上,归一化后的电离层穿透点时延改正值为 方差为
其中, 为穿透点到格网点的球面距离,rE为地球半径,hI为电离层电子密度最大处的平均高度(通常为350~400km),与φk、与λk为穿透点和格网点的纬度和经度,M为M为计算中所用的导航系统的穿透点数目之和,和是由Klobuchar模型估算的格网点k和穿透点mj的电离层天顶时延改正值;Wp,k与是第mj个穿透点到第k个格网点的权,令
步骤三:利用格网点电离层天顶时延改正值进行内插,确定用户在归一化频率f下所接收的卫星测距信号天顶电离层时延改正值
其中,Wk为加权函数;为经纬度格网点的天顶电离层时延;K为包围穿透点的可用格点数为4;但当4个格点中的某个不可用时,如剩余的3个点包围了用户穿透点,则用这3个点计算;否则,用户穿透点的天顶时延无法得到。
步骤四:用户确定不同导航系统测距信号所产生的电离层时延值Iu′;
步骤五:确定用户伪距值。
用户接收到的测距信号伪距值减去用户电离层时延值Iu′,得到用户伪距值。
通过本发明的一种适用于多模卫星导航系统的电离层格网校正方法得到的用户伪距值,与现有技术相比,可有效减小电离层时延改正误差,提高电离层格网校正精度,应用该用户伪距值进行定位计算时可以大大提高用户定位精度。
本发明的优点在于:
(1)综合运用多个导航系统的数据,提高电离层格网校正法的服务时间和服务范围;
(2)减小电离层时延改正误差;
(3)提高电离层格网校正精度;
(4)提高用户定位精度。
附图说明
图1是本发明在频率归一化到指定频率后电离层格网校正方法的流程图;
图2是本发明中国区域单时刻可由IPP内插出的IGP分布;
图3是本发明IGP(40°N,115°E)电离层垂直改正值及其误差;
图4是本发明IGP(40°N,115°E)电离层垂直改正值误差百分比;
图5是本发明电离层穿透点示意图;
图6是本发明电离层四点内插法示意图;
图7是本发明电离层三点内插法示意图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明是一种应用于多模卫星导航系统的电离层格网校正方法,如图1所示,通过五个步骤来实现:
根据经纬度格网的分布,选择加权方法,拟和满足条件的格网点的电离层天顶时延改正值为:
为用户电离层穿透点的天顶电离层时延,Wk为加权函数;为经纬度格网点的天顶电离层时延。K为包围穿透点的可用格点数,通常为4。但当4个格点中的某个不可用时,如剩余的3个点包围了用户穿透点,则用这3个点进行计算。否则,用户穿透点的天顶时延不处理。
步骤四:用户确定不同导航系统测距信号所产生的电离层时延值Iu′;
步骤五:确定用户伪距值。用户接收到的测距信号伪距值减去用户电离层时延值Iu′,得到用户伪距值。
通过本发明的一种适用于多模卫星导航系统的电离层格网校正方法得到的用户伪距值,与现有技术相比,减小电离层时延改正误差,提高电离层格网校正精度,应用该用户伪距值进行定位计算时可以大大提高用户定位精度。
实施例:
已知广域监测站接收来自GPS和北斗(BD)卫星系统的信号,中国区域内经纬度格网分布为5°×5°,GPS卫星和广域监测站之间频率f1上电离层延迟改正值为方差为北斗卫星和广域监测站之间频率f2上的电离层延迟改正值为方差为计算中国区域内用户的电离层延迟值。
步骤一:确定观测的卫星导航系统,将电离层时延改正值及其方差归一化到指定的频率上。
归一化频率f=f1下的电离层时延改正值
方差为
步骤二:根据归一化的电离层穿透点的时延值及其方差分别进行加权计算,确定指定频率下格网点的电离层天顶时延改正值。
计算5°×5°经纬度格网点的电离层天顶时延改正值;
dij为穿透点i到格网点j的球面距离。
以单个观测时刻GPS时刻84600s为例,如图2所示,圆圈表示单由GPS系统的穿透点计算出的中国区域可用格网点,叉号表示融合GPS&BD双系统的穿透点计算出的中国区域可用格网点,可见,融合GPS与BD双系统使得中国边缘地区获得了更多的可用格网点,扩大了电离层格网校正方法的服务范围。
如图3所示,黑色曲线表示采用单GPS系统,采用各个采样时刻的穿透点内插的电离层延迟改正值与采用更新时刻的格网点的内插延迟值之差,灰色曲线表示采用BD&GPS双系统,采用各个采样时刻的穿透点内插的电离层延迟改正值与采用更新时刻的格网点的内插延迟值之差,图中两组数据统计如表1所示:
表1 GPS系统与BD&GPS双系统电离层垂直改正值及其误差对比
GPS系统 | GPS&BD系统 | |
最大值 | 1.234 | 0.7763 |
最小值 | -1.612 | -0.6958 |
均值 | -0.0007428 | 0.0007237 |
标准差 | 0.2951 | 0.1995 |
可以看出:采用双系统比采用单系统得到的格网点内插延迟值更平稳地贴近真值,从而有效提高电离层延迟改正精度。
如图4所示,黑色和灰色曲线分别表示采用GPS系统和BD&GPS双系统采样点时刻穿透点加权计算格网点的电离层垂直延迟值偏离电离层模拟的真实值IIGP的误差百分比,即:
图4中两组数据统计如表2所示:
表2 GPS系统与BD&GPS双系统电离层垂直改正值误差百分比对比
GPS系统 | BD&GPS双系统 | |
最大值 | 0.7025 | 0.3477 |
均值 | 0.07841 | 0.05532 |
标准差 | 0.06867 | 0.04682 |
由表2和图4可以看出,采用双系统比采用单系统得到的格网点电离层改正误差百分比大大减小,使其值更贴近真值,从而可以减小格网点电离层垂直改正误差,进而减小用户的电离层延迟改正误差,提高定位精度。
步骤三:利用格网点电离层天顶时延改正值进行内插,确定用户在归一化频率下所接收的卫星测距信号电离层时延改正值。
1)根据用户的经纬度坐标来计算导航卫星测距信号的电离层穿透点的足下点的地理纬度φpp和经度λpp,如图5所示:
φpp=λu+ψppcosAZ(deg)
其中,ψpp为地心角,φu和λu是用户的纬度和经度,E1和AZ分别是监测站至卫星的仰角和方位角。
Wj(x,y)=x2y2(9-6x-6y+4xy)
上述内插法中Wj分两种情况:
情况一:K=4,电离层四点内插法如图6所示:
x=Δλpp/(λ2-λ1)、y=Δφpp/(φ2-φ1)、Δλpp=λPP-λ1、Δφpp=φPP-φ1,λ2与λ1为用于内插的电离层格网点经度,φ1与φ2为用于内插的电离层格网点纬度。
情况二:K=3,当只有三个格点可用且穿透点在三个格点组成的三角形内时,使用三角形内插法计算穿透点的电离层延迟。
内插法如图7所示,加权函数为:
x=Δλpp/(λ2-λ1),y=Δφpp/(φ2-φ1),λ2,λ1为用于内插的电离层格网点经度,φ1,φ2为用于内插的电离层格网点纬度,Δλpp和Δφpp分别为穿透点到用于内插的三角形直角顶点的经度间隔和纬度间隔。若K≤2,则用户穿透点的天顶时延值不可得到。
步骤四:确定用户接收不同频率导航系统的导航测距信号所产生的电离层时延改正值。用户穿透点的电离层天顶时延改正值乘以倾斜因子 得到用户在f1下的电离层时延改正值Iu。若用户接收频率为fBD1,则相应的用户电离层时延值为
步骤五:用户接收到的导航测距信号产生的电离层时延改正值减去用户电离层时延值Iu′,得到经消除后的用户电离层时延值。
Claims (2)
1、一种应用于多模卫星导航系统的电离层格网校正方法,其特征在于通过以下步骤来实现:
步骤一:首先根据来自不同卫星导航系统的数据,将电离层时延改正值及其方差归一化到用户任意指定的频率f上;设第m个导航系统穿透点电离层时延改正值为和方差为相应频率为fm,将和归一化到指定频率f上,归一化后的电离层穿透点时延改正值为 方差为
根据经纬度格网的分布,选择加权方法,拟和满足条件的格网点的电离层天顶时延改正值为:
其中, 为穿透点到格网点的球面距离,rE为地球半径,hI为电离层电子密度最大处的平均高度,与φk、与λk为穿透点和格网点的纬度和经度,M为计算中所用的所有导航系统的穿透点数目之和,和是由Klobuchar模型估算的格网点k和穿透点mj的电离层天顶时延改正值;Wpk与是第mj个穿透点到第k个格网点的权,令
步骤四:用户确定不同导航系统测距信号所产生的电离层时延值
将用户穿透点的天顶电离层时延改正值乘以倾斜因子Fu,得到用户在频率f下的电离层时延值Iu,用户接收频率fu的导航系统测距信号所产生的电离层时延值为
步骤五:确定用户伪距值;
用户接收到的测距信号伪距值减去用户电离层时延值得到用户伪距值。
2、根据权利要求1所述一种应用于多模卫星导航系统的电离层格网校正方法,其特征在于:步骤三中所述当4个格点中的某个不可用时,如剩余的3个点包围了用户穿透点,则用这3个点计算;否则,用户穿透点的天顶时延不处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100793642A CN101487883B (zh) | 2009-03-09 | 2009-03-09 | 一种适用于多模卫星导航系统的电离层格网校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100793642A CN101487883B (zh) | 2009-03-09 | 2009-03-09 | 一种适用于多模卫星导航系统的电离层格网校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101487883A true CN101487883A (zh) | 2009-07-22 |
CN101487883B CN101487883B (zh) | 2012-11-14 |
Family
ID=40890850
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100793642A Expired - Fee Related CN101487883B (zh) | 2009-03-09 | 2009-03-09 | 一种适用于多模卫星导航系统的电离层格网校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101487883B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101893714A (zh) * | 2010-07-09 | 2010-11-24 | 中国科学院测量与地球物理研究所 | 全球卫星导航系统广播电离层时延修正方法 |
CN102928850A (zh) * | 2012-11-22 | 2013-02-13 | 桂林电子科技大学 | 一种广域电离层误差改正新方法 |
CN103076615A (zh) * | 2011-12-27 | 2013-05-01 | 东莞市泰斗微电子科技有限公司 | 一种北斗导航系统电离层延迟估计方法及相应接收装置 |
CN103353600A (zh) * | 2013-08-05 | 2013-10-16 | 东莞市泰斗微电子科技有限公司 | 一种用于卫星导航系统的信息处理方法、信息处理模块及终端 |
CN103443651A (zh) * | 2011-03-23 | 2013-12-11 | 古野电气株式会社 | 定位装置、全球导航卫星系统接收机、信息终端设备、定位方法以及程序 |
CN103792546A (zh) * | 2012-10-31 | 2014-05-14 | 中国科学院光电研究院 | 一种增量式电离层折射误差修正方法 |
CN105911572A (zh) * | 2016-04-11 | 2016-08-31 | 中国人民解放军国防科学技术大学 | 北斗接收机单频电离层改正选择方法 |
CN106154268A (zh) * | 2016-09-06 | 2016-11-23 | 合肥工业大学 | 一种基于子带分割的低波段isar电离层校正方法 |
CN106405589A (zh) * | 2016-06-24 | 2017-02-15 | 西安科技大学 | 一种全球电离层格网模型的确定方法及装置 |
CN108254762A (zh) * | 2016-12-28 | 2018-07-06 | 千寻位置网络有限公司 | 伪距差分定位方法及系统 |
CN108491479A (zh) * | 2018-03-12 | 2018-09-04 | 东南大学 | 基于地图采集的道路路段经纬度与桩号匹配方法 |
CN111149017A (zh) * | 2017-09-26 | 2020-05-12 | 罗伯特·博世有限公司 | 用于提供用于卫星导航的校正数据的方法和装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6356232B1 (en) * | 1999-12-17 | 2002-03-12 | University Corporation For Atmospheric Research | High resolution ionospheric technique for regional area high-accuracy global positioning system applications |
CN1182408C (zh) * | 2002-03-18 | 2004-12-29 | 中国科学院武汉物理与数学研究所 | 一种用于卫星导航系统中电离层折射时延改正的方法 |
FR2849209B1 (fr) * | 2002-12-19 | 2007-04-06 | Agence Spatiale Europeenne | Procede et systeme de navigation en temps reel a l'aide de signaux radioelectriques a trois porteuses emis par des satellites et de corrections ionospheriques |
CN100437142C (zh) * | 2007-07-12 | 2008-11-26 | 北京航空航天大学 | 基于地基增强系统的电离层延迟误差分离方法及地基增强系统 |
-
2009
- 2009-03-09 CN CN2009100793642A patent/CN101487883B/zh not_active Expired - Fee Related
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101893714A (zh) * | 2010-07-09 | 2010-11-24 | 中国科学院测量与地球物理研究所 | 全球卫星导航系统广播电离层时延修正方法 |
CN101893714B (zh) * | 2010-07-09 | 2013-01-23 | 中国科学院测量与地球物理研究所 | 全球卫星导航系统广播电离层时延修正方法 |
CN103443651B (zh) * | 2011-03-23 | 2016-08-10 | 古野电气株式会社 | 定位装置、全球导航卫星系统接收机、信息终端设备、定位方法 |
US9488735B2 (en) | 2011-03-23 | 2016-11-08 | Furuno Electric Co., Ltd. | Positioning apparatus, GNSS receiver, information terminal device, positioning method and program |
CN103443651A (zh) * | 2011-03-23 | 2013-12-11 | 古野电气株式会社 | 定位装置、全球导航卫星系统接收机、信息终端设备、定位方法以及程序 |
CN103076615B (zh) * | 2011-12-27 | 2015-09-23 | 泰斗微电子科技有限公司 | 一种北斗导航系统电离层延迟估计方法及相应接收装置 |
CN103076615A (zh) * | 2011-12-27 | 2013-05-01 | 东莞市泰斗微电子科技有限公司 | 一种北斗导航系统电离层延迟估计方法及相应接收装置 |
CN103792546A (zh) * | 2012-10-31 | 2014-05-14 | 中国科学院光电研究院 | 一种增量式电离层折射误差修正方法 |
CN103792546B (zh) * | 2012-10-31 | 2016-03-09 | 中国科学院光电研究院 | 一种增量式电离层折射误差修正方法 |
CN102928850B (zh) * | 2012-11-22 | 2014-04-09 | 桂林电子科技大学 | 一种广域电离层误差改正新方法 |
CN102928850A (zh) * | 2012-11-22 | 2013-02-13 | 桂林电子科技大学 | 一种广域电离层误差改正新方法 |
CN103353600A (zh) * | 2013-08-05 | 2013-10-16 | 东莞市泰斗微电子科技有限公司 | 一种用于卫星导航系统的信息处理方法、信息处理模块及终端 |
CN105911572A (zh) * | 2016-04-11 | 2016-08-31 | 中国人民解放军国防科学技术大学 | 北斗接收机单频电离层改正选择方法 |
CN105911572B (zh) * | 2016-04-11 | 2018-06-05 | 中国人民解放军国防科学技术大学 | 北斗接收机单频电离层改正选择方法 |
CN106405589A (zh) * | 2016-06-24 | 2017-02-15 | 西安科技大学 | 一种全球电离层格网模型的确定方法及装置 |
CN106154268A (zh) * | 2016-09-06 | 2016-11-23 | 合肥工业大学 | 一种基于子带分割的低波段isar电离层校正方法 |
CN108254762A (zh) * | 2016-12-28 | 2018-07-06 | 千寻位置网络有限公司 | 伪距差分定位方法及系统 |
CN108254762B (zh) * | 2016-12-28 | 2021-07-27 | 千寻位置网络有限公司 | 伪距差分定位方法及系统 |
CN111149017A (zh) * | 2017-09-26 | 2020-05-12 | 罗伯特·博世有限公司 | 用于提供用于卫星导航的校正数据的方法和装置 |
CN108491479A (zh) * | 2018-03-12 | 2018-09-04 | 东南大学 | 基于地图采集的道路路段经纬度与桩号匹配方法 |
CN108491479B (zh) * | 2018-03-12 | 2021-09-03 | 东南大学 | 基于地图采集的道路路段经纬度与桩号匹配方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101487883B (zh) | 2012-11-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101487883B (zh) | 一种适用于多模卫星导航系统的电离层格网校正方法 | |
EP2746811B1 (en) | Methods for generating accuracy information on an ionosphere model for satellite navigation applications | |
CN107861131B (zh) | 一种斜路径电离层延迟的获取方法及系统 | |
CN104101888B (zh) | 一种适合于中国区域的gnss电离层延迟精确建模方法 | |
JP6290539B2 (ja) | 電離層の電子数を推定する適応的方法 | |
JP5122066B2 (ja) | 積率限界を利用してデータ処理システムの動作のインテグリティを評価する方法及び装置 | |
KR101761782B1 (ko) | 보호 수준 결정 시스템 및 그 방법 | |
US11460583B2 (en) | Method and apparatus for providing correction data for satellite navigation | |
CN101419275B (zh) | 基于多接收机的局域机场监视方法和系统 | |
CN108919316B (zh) | 一种基于局部球对称假设的单站多系统硬件延迟估计方法 | |
EP3124998B1 (en) | Positioning device | |
CN103076615B (zh) | 一种北斗导航系统电离层延迟估计方法及相应接收装置 | |
CN109975749A (zh) | 一种校正源存在条件下的短波单站直接定位方法 | |
CN101140322A (zh) | 用于布局卫星导航地基区域增强系统的装置及方法 | |
CN110146904B (zh) | 一种适用于区域电离层tec的精确建模方法 | |
CN115047499B (zh) | 星载gnss-r土壤温度与湿度的反演方法、系统 | |
EP3285091B1 (en) | Information processing device | |
CN104391279A (zh) | 基于电离层传播特性的相径扰动抑制方法 | |
Perez et al. | Status of NeQuick G after the solar maximum of cycle 24 | |
CN101419274B (zh) | 电离层延迟误差的获取方法及系统 | |
Park et al. | Performance of ionospheric maps in support of long baseline GNSS kinematic positioning at low latitudes | |
Hoque et al. | A new approach for mitigating ionospheric mapping function errors | |
Hoque et al. | An ionosphere broadcast model for next generation GNSS | |
CN114155347B (zh) | 基于菲波那契网格的星载无源定位误差分布图投影方法 | |
CN115166779A (zh) | 一种北斗导航卫星信号地基监测方法及系统 |
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: 20121114 Termination date: 20160309 |
|
CF01 | Termination of patent right due to non-payment of annual fee |