CN106873009B - 利用非组合ppp辅助的长基线模糊度解算方法 - Google Patents
利用非组合ppp辅助的长基线模糊度解算方法 Download PDFInfo
- Publication number
- CN106873009B CN106873009B CN201710023120.7A CN201710023120A CN106873009B CN 106873009 B CN106873009 B CN 106873009B CN 201710023120 A CN201710023120 A CN 201710023120A CN 106873009 B CN106873009 B CN 106873009B
- Authority
- CN
- China
- Prior art keywords
- satellite
- ambiguity
- delay
- station
- double
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 239000005433 ionosphere Substances 0.000 claims abstract description 28
- 238000005516 engineering process Methods 0.000 claims abstract description 14
- 239000005436 troposphere Substances 0.000 claims abstract description 14
- 238000000926 separation method Methods 0.000 claims abstract description 11
- 239000012298 atmosphere Substances 0.000 claims abstract description 7
- 230000001934 delay Effects 0.000 claims abstract description 7
- 239000011541 reaction mixture Substances 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 11
- 239000004576 sand Substances 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000002407 reforming Methods 0.000 claims description 8
- 230000007812 deficiency Effects 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000003379 elimination reaction Methods 0.000 claims description 5
- NUHSROFQTUXZQQ-UHFFFAOYSA-N isopentenyl diphosphate Chemical compound CC(=C)CCO[P@](O)(=O)OP(O)(O)=O NUHSROFQTUXZQQ-UHFFFAOYSA-N 0.000 claims description 5
- 230000036962 time dependent Effects 0.000 claims description 5
- 238000010521 absorption reaction Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 4
- 238000002242 deionisation method Methods 0.000 claims description 4
- 230000011218 segmentation Effects 0.000 claims description 4
- 238000013461 design Methods 0.000 abstract description 9
- 238000011156 evaluation Methods 0.000 abstract description 4
- 239000010410 layer Substances 0.000 description 29
- 230000000694 effects Effects 0.000 description 7
- 238000012545 processing Methods 0.000 description 5
- CLOMYZFHNHFSIQ-UHFFFAOYSA-N clonixin Chemical compound CC1=C(Cl)C=CC=C1NC1=NC=CC=C1C(O)=O CLOMYZFHNHFSIQ-UHFFFAOYSA-N 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 2
- 240000004808 Saccharomyces cerevisiae Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000012299 nitrogen atmosphere Substances 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000979 retarding effect Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
一种利用非组合PPP辅助的长基线模糊度解算方法,该方法包括以下步骤:利用非组合PPP技术提取各站的天顶对流层延迟和站星斜向电离层总延迟;电离层总电子含量与卫星仪器偏差的分离,得到绝对斜向电离层延迟;利用对流层投影函数将天顶对流层延迟转换到站星方向的斜向对流层延迟;将所有共视卫星的电离层延迟和对流层延迟形成双差;在长基线双差定位模型中,将双差对流层延迟和双差电离层延迟作为约束观测,联合相位和伪距观测模型,进行双差模糊度快速解算。本设计克服了长基线大气残余误差对模糊度解算的影响,不仅缩短了模糊度解算时间,而且提高了固定成功率。
Description
技术领域
本发明属于“测绘科学与技术”学科中的“大地测量”技术领域,尤其涉及利用非组合PPP辅助的长基线模糊度解算方法,主要适用于缩短模糊度解算的时间、提高模糊度固定的成功率。
背景技术
GNSS长基线模糊度快速成功解算已成为相对定位、时间传输、大气遥感、地球动力学等大范围应用方面的瓶颈。为了获取高精度的基线分量等信息,一般采用双差定位技术,该技术的原理是利用至少两台GNSS接收机同步观测多颗卫星,通过站间差分、星间差分等,消除或减弱与卫星、接收机、大气等有关的误差,简化定位模型,获取高精度的基线分量。随着基线长度的增加,大气误差的相关性减弱,双差模型中残余大气误差对模糊度解算的影响不可忽略,使模糊度解算非常困难。
长基线模糊度快速和正确解算研究一直是GNSS定位领域的热点,其中BERNESE大学推出的BERNESE软件中长基线解算方法相对比较成熟,在国内外许多领域得到广泛应用。该软件可分别采用Quasi-Ionosphere-Free(QIF)和Wide Lane/Narrow Lane(WN)等技术实现长基线模糊度解算。概括而言,QIF技术通过附加电离层约束改善模糊度浮点解的精度,以期进一步提高模糊度固定成功率;另一方面,基于WN技术的长基线模糊度解算包括两步:第一步先利用直接归整法解算宽巷模糊度,第二步利用模糊度搜索法固定窄巷模糊度,根据基频模糊度与宽、窄巷模糊度的关系解算基频模糊度。
QIF技术的可靠性和成功率取决于电离层约束的准确性,该约束可被进一步细分为函数模型约束和随机模型约束两类。函数模型约束假定电离层集中在距地面350km的薄层上,其单层模型精度受重构水平、模型误差等影响,目前模型改正只能达到50%~60%,精度只能达到几十厘米;随机模型约束要求事先获取精度较高的电离层模型改正信息,这在实际中较为困难。WN技术通过形成组合观测值克服电离层延迟的影响,但代价是牺牲了大量的观测信息,且宽巷模糊度固定的成功率特别依赖高精度的伪距观测值。不仅如此,QIF和WN技术对对流层误差的处理是估计站间天顶对流层延迟,长基线的两个站,低高度角卫星不可能共视,形成双差时这些卫星的观测进行剔除,导致模糊度解算时间超长,一般达8-24个小时。
发明内容
本发明的目的是克服现有技术中存在的长基线模糊度解算时间长且固定成功率低的缺陷与问题,提供一种模糊度解算时间短、模糊度固定成功率高的利用非组合PPP辅助的长基线模糊度解算方法。
为实现以上目的,本发明的技术解决方案是:一种利用非组合PPP辅助的长基线模糊度解算方法,该方法包括以下步骤:
A、利用非组合PPP技术提取各站的天顶对流层延迟和斜向电离层总延迟
a、构建如下非组合PPP观测模型:
式(1)中,分别表示伪距和相位观测值,为由卫星精密星历和测站近似坐标计算的站星近似距离,S,r,j分别表示卫星、测站及频率,xr为接收机坐标改正数,为接收机和卫星之间单位矢量,τr为测站天顶对流层延迟,为投影函数,dtr,dtS分别表示接收机和卫星钟差,为测站与卫星之间的斜向电离层延迟,为系数阵,br,j,分别为接收机和卫星码偏差,为整数模糊度,δr,j、分别为接收机和卫星的相位延迟,为吸收了δr,j、的非整数模糊度,λj为波长;
b、通过参数重整消除伪距观测方程中参数之间的相关性
将式(1)中等价表示为:
其中,Br=br,2-br,1和分别为接收机和卫星的仪器偏差,br,IF和分别为接收机和卫星的消电离层组合伪距偏差,则br,IF、和(Br-BS)分别被dtr、dtS和吸收,形成如下三类新参数:
c、通过参数重整消除相位观测方程中参数之间的相关性
将式(1)中相位观测方程中的dtr、dtS和参数分别转化为和的形式,则相位观测方程中的等价表示为:
其中,为重新定义的模糊度参数;
d、各测站大气延迟误差参数估计
上述参数重整后,将卫星钟差dtS和卫星轨道采用IGS发布的卫星最终精密产品代入,并采用CODE发布的卫星码偏差以及通过接收机位置得出的则各站非组合PPP的满秩观测方程为:
其中,
经过参数消秩亏处理后,利用最小二乘法估计出测站r的天顶对流层延迟参数τr和各站星的斜向电离层延迟参数
B、电离层总电子含量与卫星仪器偏差的分离
分离斜向电离层总电子含量和卫星仪器偏差的观测方程为:
式(5)中,z′为穿刺点处的天顶距,t为观测时刻,[tk≤t≤tk+Δt],y=uIPP-uγ,λ、为地理经纬度,u为地磁纬度,下标IPP和r分别表示穿刺点和测站,参数a00、a10、a01为时间相关的分段函数,即在时间间隔[tk,tk+Δt]内为常数,每隔Δt估计一组参数,Δt取5分钟,f1、f2为观测频率;
利用式(5)求出及分段参数a00、a10、a01,进而得到测站与卫星之间的斜向电离层延迟
C、大气延迟约束辅助长基线模糊度固定
a、长基线双差观测方程为:
式(6)中,分别为双差伪距和双差相位观测值,为双差站星距,分别为双差对流层延迟误差和双差电离层延迟误差,为双差模糊度参数,j表示频率;
b、构建大气延迟约束模型
先利用对流层投影函数将τr转换到站星斜向对流层延迟再将两个测站所有共视卫星的斜向电离层延迟和对流层延迟误差形成双差,构建如下大气约束模型:
式(7)中,s1,s2表示卫星号,r1,r2表示测站号;
c、双差模糊度固定
先将式(7)中的双差对流层延迟和双差电离层延迟代入式(6)进行改正,再利用最小二乘法估计得到各卫星对的双差模糊度实数解和及其方差–协方差阵和然后定义模糊度搜索方法为:
式(8)中,为返回实数最近的整数值,为频率1的2im+1个模糊度整数备选值,为频率2的2km+1个模糊度整数备选值,im和km取值3~5;
d、双差模糊度固定正确性检验
对于搜索范围内的(2im+1)×(2km+1)个备选模糊度对先根据如下式(9)求解检验值d3(i);
其中,i=0…(2im+1)×(2km+1);
再对所有备选模糊度对的检验值d3(i)进行排序,然后将最小检验值与给定的阀值ε进行比较,ε=0.1周,当d3,min<ε,则其对应的备选模糊度对即为最终固定的正确模糊度整数解。
步骤B中,所述电离层总电子含量与卫星仪器偏差的分离先采用电离层薄层模型结合相应的投影函数实施分离,具体模型为:
其中,R为地球半径,H为电离层集中的薄层距离地球表面的高度,H=350km,z为接收机的天顶距,vTEC为斜向电离层总电子含量投影到穿刺点处的垂直电离层总电子含量;
再利用二次多项式函数描述单站电离层延迟在单天内的变化,具体模型为:
vTEC=a00(t)+a10(t)·x+a01(t)·y (11)
然后通过式(2)、(10)、(11)得到分离斜向电离层总电子含量和卫星仪器偏差的观测方程。
与现有技术相比,本发明的有益效果为:
本发明一种利用非组合PPP辅助的长基线模糊度解算方法中利用非组合PPP技术分别提取各站的大气延迟误差,可以利用每个站所有卫星的观测信息,获得的大气延迟误差比较精确,将两个差分站的精确大气延迟误差进行双差,形成的双差大气延迟误差作为长基线双差定位模型中的残余大气延迟误差估值,改正定位模型中的残余大气延迟误差,提高了模糊度参数估值的准确性,即缩小模糊度搜索的范围,提高搜索的成功率。因此,本发明不仅缩短了模糊度解算时间,而且提高了固定成功率。
附图说明
图1是本发明利用非组合PPP辅助的长基线模糊度解算方法的流程图。
图2是本发明的实施例中的以LHAS为基准站形成长基线的各GNSS测站分布图。
图3是本发明的实施例中的四个测站的天顶对流层延迟结果图。
图4是本发明的实施例中的四个测站的电离层斜向延迟结果图。
图5是本发明的实施例中的LHAS和URUM测站的CPIO及其改正效果图。
具体实施方式
以下结合附图说明和具体实施方式对本发明作进一步详细的说明。
参见图1,一种利用非组合PPP辅助的长基线模糊度解算方法,该方法包括以下步骤:
A、利用非组合PPP技术提取各站的天顶对流层延迟和斜向电离层总延迟
a、构建如下非组合PPP观测模型:
式(1)中,分别表示伪距和相位观测值,为由卫星精密星历和测站近似坐标计算的站星近似距离,S,r,j分别表示卫星、测站及频率,xr为接收机坐标改正数,为接收机和卫星之间单位矢量,τr为测站天顶对流层延迟,为投影函数,dtr,dtS分别表示接收机和卫星钟差,为测站与卫星之间的斜向电离层延迟,为系数阵,br,j,分别为接收机和卫星码偏差,为整数模糊度,δr,j、分别为接收机和卫星的相位延迟,为吸收了δr,j、的非整数模糊度,λj为波长;
b、通过参数重整消除伪距观测方程中参数之间的相关性
将式(1)中等价表示为:
其中,Br=br,2-br,1和分别为接收机和卫星的仪器偏差,br,IF和分别为接收机和卫星的消电离层组合伪距偏差,则br,IF、和(Br-BS)分别被dtr、dtS和吸收,形成如下三类新参数:
c、通过参数重整消除相位观测方程中参数之间的相关性
将式(1)中相位观测方程中的dtr、dtS和参数分别转化为和的形式,则相位观测方程中的等价表示为:
其中,为重新定义的模糊度参数;
d、各测站大气延迟误差参数估计
上述参数重整后,将卫星钟差dtS和卫星轨道采用IGS发布的卫星最终精密产品代入,并采用CODE发布的卫星码偏差以及通过接收机位置得出的则各站非组合PPP的满秩观测方程为:
其中,
经过参数消秩亏处理后,利用最小二乘法估计出测站r的天顶对流层延迟参数τr和各站星的斜向电离层延迟参数
B、电离层总电子含量与卫星仪器偏差的分离
分离斜向电离层总电子含量和卫星仪器偏差的观测方程为:
式(5)中,z′为穿刺点处的天顶距,t为观测时刻,[tk≤t≤tk+Δt],y=uIPP-uγ,λ、为地理经纬度,u为地磁纬度,下标IPP和r分别表示穿刺点和测站,参数a00、a10、a01为时间相关的分段函数,即在时间间隔[tk,tk+Δt]内为常数,每隔Δt估计一组参数,Δt取5分钟,f1、f2为观测频率;
利用式(5)求出及分段参数a00、a10、a01,进而得到测站与卫星之间的斜向电离层延迟
C、大气延迟约束辅助长基线模糊度固定
a、长基线双差观测方程为:
式(6)中,分别为双差伪距和双差相位观测值,为双差站星距,分别为双差对流层延迟误差和双差电离层延迟误差,为双差模糊度参数,j表示频率;
b、构建大气延迟约束模型
先利用对流层投影函数将τr转换到站星斜向对流层延迟再将两个测站所有共视卫星的斜向电离层延迟和对流层延迟误差形成双差,构建如下大气约束模型:
式(7)中,s1,s2表示卫星号,r1,r2表示测站号;
c、双差模糊度固定
先将式(7)中的双差对流层延迟和双差电离层延迟代入式(6)进行改正,再利用最小二乘法估计得到各卫星对的双差模糊度实数解和及其方差–协方差阵和然后定义模糊度搜索方法为:
式(8)中,为返回实数最近的整数值,为频率1的2im+1个模糊度整数备选值,为频率2的2km+1个模糊度整数备选值,im和km取值3~5;
d、双差模糊度固定正确性检验
对于搜索范围内的(2im+1)×(2km+1)个备选模糊度对先根据如下式(9)求解检验值d3(i);
其中,i=0…(2im+1)×(2km+1);
再对所有备选模糊度对的检验值d3(i)进行排序,然后将最小检验值与给定的阀值ε进行比较,ε=0.1周,当d3,min<ε,则其对应的备选模糊度对即为最终固定的正确模糊度整数解。
步骤B中,所述电离层总电子含量与卫星仪器偏差的分离先采用电离层薄层模型结合相应的投影函数实施分离,具体模型为:
其中,R为地球半径,H为电离层集中的薄层距离地球表面的高度,H=350km,z为接收机的天顶距,vTEC为斜向电离层总电子含量投影到穿刺点处的垂直电离层总电子含量;
再利用二次多项式函数描述单站电离层延迟在单天内的变化,具体模型为:
vTEC=a00(t)+a10(t)·x+a01(t)·y (11)
然后通过式(2)、(10)、(11)得到分离斜向电离层总电子含量和卫星仪器偏差的观测方程。
本发明的原理说明如下:
本设计提出先利用非组合PPP技术对各站进行精确的大气延迟误差提取,再对提取的各站大气误差形成双差,以辅助长基线模糊度快速固定的方法,主要适用于解决超长基线(大于1000km)模糊度解算中,由于站间大气误差的相关性减弱,双差后残余大气误差比较大,使得模糊度固定比较慢甚至不能正确固定的问题,从而获取高精度可靠的长基线信息。
式(1)中,对于双频GNSS观测,如果观测m颗卫星,可以组成4m个观测方程,需要解算的全部参数个数为7m+7,导致设计矩阵列秩亏,若不加以处理,则参数解不唯一。由于部分参数线性相关,在参数估计过程中是不可分的,本设计提出通过参数重整的分步消秩亏法方案,以确保参数的可估性。式(1)的伪距观测方程中,dtr-dtS和均与之间存在线性关系,则可以等价表示为:
对于相位观测方程,为了和伪距观测方程中公共未知参数保持相同的形式,则将dtr、dtS和参数分别转化为和的形式。
经过参数消秩亏处理,测站r双频接收机观测到m>4颗卫星,独立待估参数降为3m+2个,可以组成4m个类似式(4)的观测方程,利用最小二乘法可以估计出测站r的天顶对流层延迟参数τr和各卫星的斜向电离层延迟参数
由式(2)可知,步骤A中估计的电离层延迟误差参数包括绝对电离层延迟参数和吸收的仪器偏差(Br-BS),辅助超长基线模糊度解算时,必须将和仪器偏差(Br-BS)进行精确分离。一般采用电离层薄层模型结合相应的投影函数实施分离,电离层薄层模型是假设电离层集中在距离地球表面一定高度H的薄层上,同时将斜向电离层总电子含量投影到穿刺点(站星视线与薄层的交点)处的垂直电离层总电子含量vTEC,具体模型为:
式中,R为地球半径,H为薄层高度,这里选350km,z和z′分别为接收机和穿刺点处的天顶距。由于GNSS观测信号覆盖范围有限,利用二次多项式函数可有效描述单站电离层延迟在单天内的变化,模型为:
vTEC=a00(t)+a10(t)·x+a01(t)·y
其中,t为观测时刻,y=uIPP-uR,λ、为地理经纬度,u表示地磁纬度,下标IPP和r分别表示穿刺点和测站,未知参数a00、a10、a01为时间相关的分段函数,即在时间间隔[tk,tk+Δt]内为常数,每隔Δt估计一组参数,这里Δt取5分钟。
步骤C中,如果忽略表示接收机、卫星和历元的符号,则超长基线双差观测方程可以简写为:
如果两个测站双频接收机同步观测卫星为n颗,则可以形成4n-4个类似式(6)的观测方程。当基线比较短时,可以认为两个站的大气特性一致,和可以忽略,待估参数只有和基线分量参数;对于长基线,和不能忽略,使得实数解偏差比较大,给模糊度搜索带来困难。构建大气约束模型时,对于n颗同步观测卫星,可以形成2n–2个大气约束方程。双差模糊度固定时,将式(7)中的双差对流层延迟和双差电离层延迟代入式(6)进行改正,式(6)中仅剩余双差模糊度参数和基线分量参数。由于卫星和接收机的偏差参数通过双差完全消除,双差模糊度具有整数特性,则可定义式(8)所示的模糊度搜索方法。
本设计利用非组合PPP技术提取各站的大气延迟能利用全部可视卫星的观测信息,精度更高、可靠性更高,将站间精确的双差大气延迟作为虚拟观测,克服长基线残余大气误差过大使模糊度收敛时间过长(一般要几个小时)的瓶颈,加速了长基线模糊度的收敛,提高基线分量的解算精度。
实施例:
数据描述及处理方案:为了更好展示和分析本设计的效果,IGS站不仅有精确的坐标信息,而且还提供精确的大气信息,因此选择我国范围内6个IGS跟踪站2003年DOY96天(DOY:Day Of Year,年积日)的GPS数据进行处理,地理分布见图2,其中选择LHAS作为参考站,所形成的5条基线长度在1200km~3000km。各站利用非组合PPP技术提取大气延迟参数时采样的观测数据及参数设置见表一。
表一
表一中卫星截止高度角设置为5度而不是设置成以往常用的10度是确保各站提取大气延迟误差时充分利用低高度角卫星的观测值,以及ZTD参数与PPP确定的高程方向位置信息可分离。
利用非组合PPP技术提取的LHAS、URUM、TWTF和KUNM 4个IGS站的天顶对流层延迟结果见图3,图3中横轴表示UT时,纵轴为ZTD大小,单位为米,其中,灰线表示利用非组合PPP提取的ZTD(天顶对流层延迟),时间分辨率是30s,黑色三角形是IGS发布的ZTD估值,时间间隔是2h,精度达mm级,可以作为参考值。图3表明,利用非组合PPP提取的各站ZTD估值和IGS发布的ZTD参考值的变化趋势相同,且在同步时刻差异不超过1cm,各站的天顶对流层延迟一天内变化在dm量级,但是相同时刻各站的差异在dm–m量级,在长基线模糊度解算模型中,残余的对流层延迟误差是不可忽略的。
利用非组合PPP技术提取的LHAS、WUHN、TWTF和URUM 4个站的站星电离层斜向延迟见图4,图4中横轴表示地方时,纵轴表示提取的绝对电离层斜向延迟,其中,一种颜色代表一颗卫星的电离层延迟变化。图4表明,由于各站相距比较远(上千km),受地磁活动强度的差异,各站在相同时刻电离层延迟差异比较大:URUM和LHAS站在14:00前后各历元电离层延迟差异最大可达15-20m,即使对于一天内电离层活动比较平静的8:00附近,电离层延迟差异也可以达到2-3m,因此双差后残余电离层延迟对模糊度固定的影响也不可忽略。
参见图1,一种利用非组合PPP辅助的长基线模糊度解算方法,该方法包括以下步骤:
A、利用非组合PPP技术提取各站的天顶对流层延迟ZTD和斜向电离层总延迟PID
a、构建如下非组合PPP观测模型:
式(1)中,分别表示伪距和相位观测值,为由卫星精密星历和测站近似坐标计算的站星近似距离,S,r,j分别表示卫星、测站及频率,xr为接收机坐标改正数,为接收机和卫星之间单位矢量,τr为测站天顶对流层延迟,为投影函数,dtr,dtS分别表示接收机和卫星钟差,为测站与卫星之间的斜向电离层延迟,为系数阵,br,j,分别为接收机和卫星码偏差,为整数模糊度,δr,j、分别为接收机和卫星的相位延迟,为吸收了δr,j、的非整数模糊度,λj为波长;
b、通过参数重整消除伪距观测方程中参数之间的相关性
将式(1)中等价表示为:
其中,Br=br,2-br,1和分别为接收机和卫星的仪器偏差,br,IF和分别为接收机和卫星的消电离层组合伪距偏差,则br,IF、和(Br-BS)分别被dtr、dtS和吸收,形成如下三类新参数:
c、通过参数重整消除相位观测方程中参数之间的相关性
将式(1)中相位观测方程中的dtr、dtS和参数分别转化为和的形式,则相位观测方程中的等价表示为:
其中,为重新定义的模糊度参数;
d、各测站大气延迟误差参数估计
上述参数重整后,将卫星钟差dtS和卫星轨道采用IGS发布的卫星最终精密产品代入,并采用CODE发布的卫星码偏差以及通过接收机位置得出的则各站非组合PPP的满秩观测方程为:
其中,
经过参数消秩亏处理后,利用最小二乘法估计出测站r的天顶对流层延迟参数τr和各站星的斜向电离层延迟参数
B、电离层总电子含量与卫星仪器偏差的分离
先采用电离层薄层模型结合相应的投影函数实施分离,具体模型为:
其中,R为地球半径,H为电离层集中的薄层距离地球表面的高度,H=350km,z为接收机的天顶距,vTEC为斜向电离层总电子含量投影到穿刺点处的垂直电离层总电子含量;
再利用二次多项式函数描述单站电离层延迟在单天内的变化,具体模型为:
vTEC=a00(t)+a10(t)·x+a01(t)·y (11)
然后通过式(2)、(10)、(11)得到如下分离斜向电离层总电子含量和卫星仪器偏差的观测方程:
式(5)中,z′为穿刺点处的天顶距,t为观测时刻,[tk≤t≤tk+Δt],y=uIPP-uγ,λ、为地理经纬度,u为地磁纬度,下标IPP和r分别表示穿刺点和测站,参数a00、a10、a01为时间相关的分段函数,即在时间间隔[tk,tk+Δt]内为常数,每隔Δt估计一组参数,Δt取5分钟,f1、f2为观测频率;
利用式(5)求出及分段参数a00、a10、a01,进而得到测站与卫星之间的斜向电离层延迟
C、大气延迟约束辅助长基线模糊度固定
a、长基线双差观测方程为:
式(6)中,分别为双差伪距和双差相位观测值,为双差站星距,分别为双差对流层延迟误差和双差电离层延迟误差,为双差模糊度参数,j表示频率;
b、构建大气延迟约束模型
根据步骤A、B解算的两个测站天顶对流层延迟τr和斜向电离层延迟信息,先利用对流层投影函数将τr转换到站星斜向对流层延迟再将两个测站所有共视卫星的斜向电离层延迟和对流层延迟误差形成双差,构建如下大气约束模型:
式(7)中,s1,s2表示卫星号,r1,r2表示测站号;
c、双差模糊度固定
先将式(7)中的双差对流层延迟和双差电离层延迟代入式(6)进行改正,再利用最小二乘法估计得到各卫星对的双差模糊度实数解和及其方差–协方差阵和然后定义模糊度搜索方法为:
式(8)中,为返回实数最近的整数值,为频率1的2im+1个模糊度整数备选值,为频率2的2km+1个模糊度整数备选值,im和km取值3~5;
d、双差模糊度固定正确性检验
对于搜索范围内的(2im+1)×(2km+1)个备选模糊度对先根据如下式(9)求解检验值d3(i);
其中,i=0…(2im+1)×(2km+1);
再对所有备选模糊度对的检验值d3(i)进行排序,然后将最小检验值与给定的阀值ε进行比较,ε=0.1周,当d3,min<ε,则其对应的备选模糊度对即为最终固定的正确模糊度整数解。
对各站提取的大气延迟ZTD和CPIO结果进行精度评估。对各站提取的天顶对流层延迟的精度评估可以以IGS发布的结果作为参考,计算各站的外符合精度RMS,其结果见下述表二中第三列,其中SHAO站的IGS参考值缺失,相应的RMS无法评估,表二中结果表明,LHAS站的ZTD估值的外符合精度最高,达0.73cm,其余各站的ZTD精度一般均优于1.2cm。但是,对于各站的电离层延迟误差的精度评估缺乏外部电离层信息,本设计采用电离层修正后计算的L4组合观测量的稳定性作为电离层延迟误差的内符合精度STD,由于L4组合消除了与频率无关的项,如站星间的距离、钟差等,仅保留与频率有关的项,即电离层延迟误差I1和模糊度参数N1、N2,这种组合也称为Geotry-Free(GF)组合。L4组合观测模型为:
由上式可知,由于模糊度参数是时不变参数,当电离层延迟足够小时,L4应该是一个常数,表现在图中应该是平稳曲线,因此利用L4变化幅度可以作为电离层延迟改正效果的评价标准。利用原始相位观测值计算的L4组合观测量和经过前面提取的各站电离层延迟改正后的相位观测计算的L4组合观测量的结果见图5,这里只列出LHAS和URUM两个站的结果,所有站的电离层延迟误差的STD值见表二中第四列。图5中上面两幅图表示两个站原始相位观测值计算的L4组合观测量,下面两幅图是经过电离层延迟改正后的相位观测形成的L4组合观测量,由图5可知,由于LHAS和URUM两个站的地理位置比较特殊,一个在高原,一个在高纬地区,用原始观测计算L4组合观测量的电离层延迟误差比较大,而经过电离层延迟改正后的L4组合观测量在各卫星连续观测弧段表现出平稳的变化趋势,结合表二可知,其STD量级均优于8cm。
利用非组合PPP提取的大气延迟辅助长基线模糊度解算的效果。为了说明本设计方法固定的长基线模糊度的效果,将其与BERNESE中QIF模糊度解算的成功率进行对比,结果见表二中第5–6列。结果表明,本设计处理模糊度固定的成功率比BERNESE有所提高,和大气延迟提取精度联合分析我们可以看到各站模糊度固定的成功率与其电离层延迟精度相关,电离层延迟精度越高,模糊度固定的成功率越高。
表二 各站的基线解算结果汇总
Claims (2)
1.一种利用非组合PPP辅助的长基线模糊度解算方法,其特征在于,该方法包括以下步骤:
A、利用非组合PPP技术提取各站的天顶对流层延迟和斜向电离层总延迟
a、构建如下非组合PPP观测模型:
式(1)中,分别表示伪距和相位观测值,为由卫星精密星历和测站近似坐标计算的站星近似距离,S,r,j分别表示卫星、测站及频率,xr为接收机坐标改正数,为接收机和卫星之间单位矢量,τr为测站天顶对流层延迟,为投影函数,dtr,dtS分别表示接收机和卫星钟差,为测站与卫星之间的斜向电离层延迟,为系数阵,br,j,分别为接收机和卫星码偏差,为整数模糊度,δr,j、分别为接收机和卫星的相位延迟,为吸收了δr,j、的非整数模糊度,λj为波长;
b、通过参数重整消除伪距观测方程中参数之间的相关性
将式(1)中等价表示为:
其中,Br=br,2-br,1和分别为接收机和卫星的仪器偏差,br,IF和分别为接收机和卫星的消电离层组合伪距偏差,则br,IF、和(Br-BS)分别被dtr、dtS和吸收,形成如下三类新参数:
c、通过参数重整消除相位观测方程中参数之间的相关性
将式(1)中相位观测方程中的dtr、dtS和参数分别转化为和的形式,则相位观测方程中的等价表示为:
其中,为重新定义的模糊度参数;
d、各测站大气延迟误差参数估计
上述参数重整后,将卫星钟差dtS和卫星轨道采用IGS发布的卫星最终精密产品代入,并采用CODE发布的卫星码偏差以及通过接收机位置得出的则各站非组合PPP的满秩观测方程为:
其中,
经过参数消秩亏处理后,利用最小二乘法估计出测站r的天顶对流层延迟参数τr和各站星的斜向电离层延迟参数
B、电离层总电子含量与卫星仪器偏差的分离
分离斜向电离层总电子含量和卫星仪器偏差的观测方程为:
式(5)中,z′为穿刺点处的天顶距,α=(f2 2-f1 2)/f2 2,t为观测时刻,[tk≤t≤tk+Δt],y=uIPP-uγ,λ、为地理经纬度,u为地磁纬度,下标IPP和r分别表示穿刺点和测站,参数a00、a10、a01为时间相关的分段函数,即在时间间隔[tk,tk+Δt]内为常数,每隔Δt估计一组参数,Δt取5分钟,f1、f2为观测频率;
利用式(5)求出及分段参数a00、a10、a01,进而得到测站与卫星之间的斜向电离层延迟
C、大气延迟约束辅助长基线模糊度固定
a、长基线双差观测方程为:
式(6)中,分别为双差伪距和双差相位观测值,为双差站星距,分别为双差对流层延迟误差和双差电离层延迟误差,为双差模糊度参数,j表示频率;
b、构建大气延迟约束模型
先利用对流层投影函数将τr转换到站星斜向对流层延迟再将两个测站所有共视卫星的斜向电离层延迟和对流层延迟误差形成双差,构建如下大气约束模型:
式(7)中,s1,s2表示卫星号,r1,r2表示测站号;
c、双差模糊度固定
先将式(7)中的双差对流层延迟和双差电离层延迟代入式(6)进行改正,再利用最小二乘法估计得到各卫星对的双差模糊度实数解和及其方差–协方差阵和然后定义模糊度搜索方法为:
式(8)中,为返回实数最近的整数值,为频率1的2im+1个模糊度整数备选值,为频率2的2km+1个模糊度整数备选值,im和km取值3~5;
d、双差模糊度固定正确性检验
对于搜索范围内的(2im+1)×(2km+1)个备选模糊度对先根据如下式(9)求解检验值d3(i);
其中,i=0…(2im+1)×(2km+1);
再对所有备选模糊度对的检验值d3(i)进行排序,然后将最小检验值与给定的阀值ε进行比较,ε=0.1周,当d3,min<ε,则其对应的备选模糊度对即为最终固定的正确模糊度整数解。
2.根据权利要求1所述的一种利用非组合PPP辅助的长基线模糊度解算方法,其特征在于:步骤B中,所述电离层总电子含量与卫星仪器偏差的分离先采用电离层薄层模型结合相应的投影函数实施分离,具体模型为:
其中,R为地球半径,H为电离层集中的薄层距离地球表面的高度,H=350km,z为接收机的天顶距,vTEC为斜向电离层总电子含量投影到穿刺点处的垂直电离层总电子含量;
再利用二次多项式函数描述单站电离层延迟在单天内的变化,具体模型为:
vTEC=a00(t)+a10(t)·x+a01(t)·y (11)
然后通过式(2)、(10)、(11)得到分离斜向电离层总电子含量和卫星仪器偏差的观测方程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710023120.7A CN106873009B (zh) | 2017-01-12 | 2017-01-12 | 利用非组合ppp辅助的长基线模糊度解算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710023120.7A CN106873009B (zh) | 2017-01-12 | 2017-01-12 | 利用非组合ppp辅助的长基线模糊度解算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106873009A CN106873009A (zh) | 2017-06-20 |
CN106873009B true CN106873009B (zh) | 2019-08-02 |
Family
ID=59158424
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710023120.7A Active CN106873009B (zh) | 2017-01-12 | 2017-01-12 | 利用非组合ppp辅助的长基线模糊度解算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106873009B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108415049B (zh) * | 2018-04-19 | 2022-05-06 | 千寻位置网络有限公司 | 提高网络rtk双差宽巷模糊度固定正确率的方法 |
CN111158031B (zh) * | 2018-11-07 | 2022-09-06 | 千寻位置网络有限公司 | 电离层延迟的提取方法及装置、用户及服务端、定位系统 |
CN109521453B (zh) * | 2018-12-29 | 2022-01-25 | 广东电网有限责任公司 | 一种cors网长距离基线模糊度快速固定方法及装置 |
CN110596641B (zh) * | 2019-09-29 | 2023-04-18 | 中国科学院测量与地球物理研究所 | 混合los/nlos场景下基于toa定位模式的l1范数方法 |
CN111025354A (zh) * | 2019-12-20 | 2020-04-17 | 东南大学 | 基于单差分电离层加权模型的中长基线rtk定位方法 |
CN111538044B (zh) * | 2020-04-10 | 2023-03-07 | 苏州市高新北斗导航平台有限公司 | 低成本rtk接收机定位精度测试方法 |
CN112382844B (zh) * | 2020-11-13 | 2023-04-25 | 成都戎星科技有限公司 | 一种低轨卫星通信系统的天线伺服电机控制方法和系统 |
CN112649826A (zh) * | 2021-01-12 | 2021-04-13 | 中国石油大学(华东) | 长距离rtk与ppp-rtk整周模糊度解算方法 |
CN115144878B (zh) * | 2022-07-13 | 2023-05-30 | 武汉大学 | 基于ppp的短距离大高差nrtk对流层延迟改正方法 |
CN116609799B (zh) * | 2023-07-20 | 2023-10-20 | 武汉大学 | 一种厘米级斜向电离层延迟产品的生成方法及装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7961147B1 (en) * | 2008-07-25 | 2011-06-14 | Rockwell Collins, Inc. | Long baseline phase interferometer ambiguity resolution using frequency differences |
CN102650694B (zh) * | 2011-02-25 | 2014-03-12 | 中国人民解放军61081部队 | 基于北斗四频信号的中长基线模糊度解算方法 |
CN104749597B (zh) * | 2015-04-13 | 2017-05-03 | 中国科学院光电研究院 | 一种针对北斗系统的中长基线下模糊度解算方法 |
US11187808B2 (en) * | 2015-06-19 | 2021-11-30 | Profound Positioning Inc. | Methods and systems for performing global navigation satellite system (GNSS) orbit and clock augmentation and position determination |
CN105467406B (zh) * | 2015-12-18 | 2017-10-24 | 中国科学院测量与地球物理研究所 | 多测站广播星历确定用户测距精度上界值的方法 |
-
2017
- 2017-01-12 CN CN201710023120.7A patent/CN106873009B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106873009A (zh) | 2017-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106873009B (zh) | 利用非组合ppp辅助的长基线模糊度解算方法 | |
CN109581452B (zh) | 一种gnss参考站载波相位整周模糊度解算方法 | |
CN109828288A (zh) | 一种基于区域cors的实时电离层建模与监测方法 | |
CN105607089B (zh) | 一种基于伪距辅助载波相位的北斗共视时间传递的方法 | |
CN106842268B (zh) | 双gnss接收机载波相位双差整周模糊度浮点解向量估计方法 | |
CN110275186A (zh) | Leo卫星增强的gnss电离层归一化与融合建模方法 | |
CN105629263A (zh) | 一种对流层大气延迟误差估计改正方法和改正系统 | |
CN110109158A (zh) | 基于gps、glonass和bds多系统的事后超快速rtk定位算法 | |
CN114019584B (zh) | 一种大高差地区高精度cors网vrs解算方法 | |
CN116299615B (zh) | 一种实现单北斗实时ppp模糊固定的相位偏差估计方法 | |
CN111551971B (zh) | 一种支持异频gnss信号伪距差分定位的方法 | |
CN115079236B (zh) | 通过低轨增强缩短广域非差非组合ppp-rtk定位收敛时间的方法 | |
CN114779301B (zh) | 一种基于广播星历的卫星导航实时精密单点定位方法 | |
CN115902968A (zh) | 基于北斗三号geo播发增强信息的ppp终端定位方法 | |
Soycan et al. | Precise point positioning versus traditional solution for GNSS networks | |
CN116203598A (zh) | 基于地基与星基增强融合的电离层建模方法、装置和介质 | |
CN115933356A (zh) | 一种虚拟原子钟的高精度时间同步系统和方法 | |
CN105182367A (zh) | 一种提取bds卫星电离层穿刺点电子浓度的新方法 | |
CN104502943B (zh) | 一种室内伪卫星差分相对定位方法 | |
CN108508470A (zh) | 面向全球无缝的瞬时分米级导航定位方法 | |
CN100371731C (zh) | Gps和伪卫星组合定位方法 | |
CN115980317B (zh) | 基于修正相位的地基gnss-r数据土壤水分估算方法 | |
CN104991265A (zh) | 一种北斗卫星导航系统用户统一性定位方法 | |
CN106814376B (zh) | 一种快速精确厘米级单点定位方法 | |
CN109738912B (zh) | 一种基于gnss卫星信号实现固定点授时的方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |