CN105182388B - 一种快速收敛的精密单点定位方法 - Google Patents
一种快速收敛的精密单点定位方法 Download PDFInfo
- Publication number
- CN105182388B CN105182388B CN201510660882.9A CN201510660882A CN105182388B CN 105182388 B CN105182388 B CN 105182388B CN 201510660882 A CN201510660882 A CN 201510660882A CN 105182388 B CN105182388 B CN 105182388B
- Authority
- CN
- China
- Prior art keywords
- msubsup
- mtd
- msub
- mrow
- mover
- 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 111
- 238000001914 filtration Methods 0.000 claims abstract description 19
- 239000005433 ionosphere Substances 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 22
- 239000013598 vector Substances 0.000 claims description 18
- 230000001934 delay Effects 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 6
- 239000005436 troposphere Substances 0.000 claims description 6
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000004804 winding Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 abstract description 6
- 238000012545 processing Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 5
- 230000003595 spectral effect Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000002242 deionisation method Methods 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000005295 random walk Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- 229910052791 calcium Inorganic materials 0.000 description 1
- 239000011575 calcium Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 239000011541 reaction mixture Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000007704 transition Effects 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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
- G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种快速收敛的精密单点定位方法,首先对全球导航卫星的观测数据进行预处理,再联合原始载波相位观测值的观测方程、码相位半合组合观测值的观测方程和几何无关组合观测值的观测方程构建改进的精密单点定位的观测方程,采用卡尔曼滤波对改进的精密单点定位的观测方程中的参数进行估计,得到接收机位置信息和对流层延迟、电离层延迟信息,实现GNSS快速单点精密定位。与现有技术相比:本发明一种快速收敛的精密单点定位方法能更好地降低观测噪声、轨道误差对滤波收敛速度与定位精度影响,并利用对流层延迟、电离层延迟历元间缓慢变化的特性作为约束条件,提高了精密单点定位解算的性能。
Description
技术领域
本发明涉及全球导航卫星定位技术领域,尤其涉及一种快速收敛的精密单点定位方法。
背景技术
精密单点定位(Precise Point Positioning,简称PPP)是基于状态空间域改正信息的高精度定位方法,基本实现方法为:用户采用一台GNSS(Global NavigationSatellite System)接收机的双频码伪距和载波相位观测值,通过国际GNSS服务组织(International GNSS Service,简称IGS)提供的GNSS精密轨道和精密钟差产品内插出相应观测时刻的卫星轨道和钟差信息,同时应用误差模型修正天线相位中心、地球自转、固体潮汐等误差源的影响,进行单站绝对定位。目前,位置解收敛到静态厘米级、动态分米级的精度通常需要较长的时间,如静态模式下一般需要20~30分钟才能收敛到厘米级。收敛速度慢是限制PPP技术更加广泛应用的主要原因。
自从精密单点定位技术被提出以来,许多专家学者提出了多种定位方法,如Kouba在期刊“GPS Solutions”2001年第5卷第2期中提出了标准消电离层组合方法(Un-Difference ionosphere-free combined method,简称UD方法),Gao Y在期刊“Navigation”2002年第49卷第2期中提出了UofC方法(University of Calgary method,简称UofC方法)。这两种方法均采用了消电离层组合,虽然降低了电离层延迟对定位的影响,但放大了观测噪声和多路径效应,不利于位置解的快速收敛;张宝成在“测绘学报”2010年第39卷第5期中提出了非差非组合方法(Un-Combined method,简称UC方法),该方法采用参数估计形式将电离层延迟模拟成随机游走,但采用了原始码伪距观测值,观测噪声仍然较大,且无法降低轨道误差对定位的影响。另外,Ge M在期刊“Journal of Geodesy”2008年82卷第7期中、张小红在期刊“武汉大学学报·信息科学版”2010年第35卷第6期中、Shi J在期刊“GPS Solutions”2014年第66卷第1期中均实现了基于相位偏差改正的精密单点定位固定解,虽在一定程度上提高了PPP收敛速度和定位精度,但参数的解算仍需基于上述三种基本方法来实现。
发明内容
本发明所要解决的技术问题是提供一种快速收敛的精密单点定位方法,有效降低观测噪声、轨道误差对滤波收敛速度与定位精度影响,并利用对流层延迟、电离层延迟历元间缓慢变化的特性作为约束条件,提高了精密单点定位解算的性能。
为了解决上述技术问题,本发明采用如下技术方案:
一种快速收敛的精密单点定位方法,包括如下步骤:
(S1)、对全球导航卫星的观测数据进行预处理:
采用M-W组合法和电离层残差法对原始载波相位观测值中存在的周跳进行探测,并修复;
对天线相位中心偏差、地球自转、相对论效应、天线相位缠绕误差源的影响进行修正;
对精密卫星轨道和卫星钟差数据进行拉格朗日多项式内插,获取观测时刻的卫星轨道和卫星钟差;
(S2)、顾及电离层延迟对不同频率观测值的影响关系,建立卫星的原始载波相位观测值和的观测方程,如公式(5)、公式(6)所示:
其中,s表示卫星,和表示载波相位观测值,ρ表示卫星与接收机之间的几何距离,c表示真空中光的速度,表示接收机钟差,表示卫星钟差,分别表示对流层干、湿迟延映射函数;δzpd、δzpw分别表示对流层天顶方向干延迟和湿延迟,为载波相位观测值中站星视线方向的电离层延迟,分别为原始载波相位观测值和的模糊度,f1、f2为载波相位的频率,λ1、λ2表示载波相位的波长,表示载波相位观测值的观测噪声;
(S3)、根据原始码伪距观测值和原始载波相位观测值构建码相位半合组合观测值的观测方程,如公式(7)所示:
(S4)、根据原始载波相位观测值和构建几何无关组合观测值的观测方程,如公式(8)所示:
(S5)、联合原始载波相位观测值和的观测方程、码相位半合组合观测值的观测方程和几何无关组合观测值的观测方程构建改进的精密单点定位的观测方程,如公式(9)所示,
(S6)、采用卡尔曼滤波对改进的精密单点定位的观测方程中的参数进行估计,得到接收机位置信息和对流层延迟、电离层延迟信息。
作为上述方案的进一步优化,对步骤(S5)的改进的精密单点定位的观测方程进行线性化处理,得到公式(10):
y(k)=A·X(k)+εy,εy~N(0,Ωy) (10)
其中,y(k)为观测向量,其维数为4j,j为历元k同步观测到的卫星数量;X(k)为待估参数向量,维数为(5+3j),A为设计矩阵,εy为测量噪声向量,Ωy为测量噪声向量εy的协方差阵;y(k)、X(k)和A表达式分别为
其中,X、Y、Z为接收机三维位置,表示克罗内克积,e4为各元素均为1的4维列向量,Ij为j维单位矩阵,B、C和D子矩阵的含义分别为:
其中,l、m、n分别为接收机到卫星的方向余弦。
作为上述方案的进一步优化,步骤(S5)得到的改进的精密单点定位的观测方程第k历元对应的协方差阵Ωy(k),如公式(11)所示:
其中, 表示克罗内克积,Ij为j维单位矩阵,分别为卫星原始载波相位观测值Φ1、Φ2,原始码伪距观测值P2的测量噪声的方差。
与已有技术相比,本发明的一种快速收敛的精密单点定位方法的有益效果体现在:
1、本发明的一种快速收敛的精密单点定位方法有效地降低观测噪声、轨道误差对滤波收敛速度与定位精度影响,提高了精密单点定位解算的性能。
2、本发明的一种快速收敛的精密单点定位方法,不仅降低了观测噪声和轨道误差对定位的影响,还可以利用大气延迟平稳变化的特性作为约束提高滤波收敛速度,其构建的改进的精密单点定位的观测方程,具有更好的模型结构,采用卡尔曼滤波对改进的精密单点定位的观测方程进行参数估计,能克服现有解算方法的技术缺陷。
3、本发明的一种快速收敛的精密单点定位方法,能快速地实现精密单点定位厘米级位置解,可提高精密单点定位收敛速度以及定位精度。
附图说明
图1是本发明的一种快速收敛的精密单点定位方法的流程图。
图2为本发明的应用实施例采用的IGS观测站的分布图。
图3为使用本发明方法和UC方法对公共滤波收敛数据处理后的收敛时间统计结果图。
图4为使用本发明方法和UC方法对公共滤波收敛数据处理后的平面位置偏差统计结果图。
图5为使用本发明方法和UofC方法对公共滤波收敛数据处理后的收敛时间统计结果图。
图6为使用本发明方法和UofC方法对公共滤波收敛数据处理后的平面位置偏差统计结果图。
图7为使用本发明方法和UD方法对公共滤波收敛数据处理后的收敛时间统计结果图。
图8为使用本发明方法和UD方法对公共滤波收敛数据处理后的平面位置偏差统计结果图。
具体实施方式
下面对本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
一种快速收敛的精密单点定位方法,包括如下步骤:参见图1,图1是本发明的一种快速收敛的精密单点定位方法的流程图。
(S1)、对全球导航卫星的观测数据进行预处理:
采用M-W组合法和电离层残差法对原始载波相位观测值中存在的周跳进行探测,并修复;
对天线相位中心偏差、地球自转、相对论效应、天线相位缠绕误差源的影响进行修正;
对精密卫星轨道和卫星钟差数据进行拉格朗日多项式内插,获取观测时刻的卫星轨道和卫星钟差。
(S2)、建立卫星的原始载波相位观测值和的观测方程:
在GNSS应用中,在卫星s的原始码伪距和载波相位观测值的观测方程表达式如下:
其中,s表示卫星(s=1,…,j),i表示载波相位和码伪距的频率数,表示码伪距观测值,表示以米为单位的载波相位观测值,ρ表示卫星与接收机之间的几何距离,c表示真空中光的速度,δtr表示接收机钟差,δts表示卫星钟差,分别表示对流层干、湿迟延映射函数,δzpd、δzpw分别表示对流层天顶方向干延迟和湿延迟,表示站星视线方向的电离层延迟,dr,i表示接收机端的硬件码延迟,表示卫星端的硬件码延迟,表示载波相位的模糊度,λi表示载波相位的波长,表示接收机端的相位延迟,表示卫星端的相位延迟,分别表示码伪距和载波相位观测值的观测噪声;
将式(1)和式(2)简化为
其中,分别表示重新定义后的接收机钟差、卫星钟差和载波相位模糊度,表达式分别为
顾及电离层延迟对不同频率观测值的影响关系为则得到原始载波相位观测值和的观测方程,表达式为
式(5)或式(6)中,和表示载波相位观测值,f1、f2为载波相位的频率,λ1、λ2表示载波相位的波长,表示载波相位观测值的观测噪声;I1为载波相位观测值中站星视线方向的电离层延迟,分别为原始载波相位观测值和的模糊度;
(S3)、根据原始码伪距观测值和原始载波相位观测值构建码相位半合组合观测值的观测方程,如公式(7)所示:
(S4)、根据原始载波相位观测值和构建几何无关组合观测值的观测方程,如公式(8)所示:
(S5)、联合原始载波相位观测值和的观测方程、码相位半合组合观测值的观测方程和几何无关组合观测值的观测方程构建改进的精密单点定位的观测方程,如公式(9)所示:
式(9)线性化后的矩阵形式为:
y(k)=A·X(k)+εy,εy~N(0,Ωy) (10)
式(10)中,y(k)为观测向量,其维数为4j,j为历元k同步观测到的卫星数量;X(k)为待估参数向量,维数为(5+3j),A为设计矩阵,εy为测量噪声向量,Ωy为测量噪声向量εy的协方差阵;y(k)、X(k)和A表达式分别为
其中,X、Y、Z为接收机三维位置,表示克罗内克积,e4为各元素均为1的4维列向量,Ij为j维单位矩阵,B、C和D子矩阵的含义分别为:
其中,l、m、n分别为接收机到卫星的方向余弦。
观测值协方差阵衡量了观测值的精度及相关性。由于卫星观测值受到噪声及多路径效应影响主要发生在低高度角卫星,为了不使高仰角卫星的观测值降低权重,精密单点定位中一般采用高度角定权。不考虑码和载波之间的相关性,即则步骤(S5)得到的改进的精密单点定位的观测方程第k历元对应的协方差阵Ωy(k),如公式(11)所示:
式(11)中, 表示克罗内克积,Ij为j维单位矩阵,分别为卫星原始载波相位观测值Φ1、Φ2,原始码伪距观测值P2的测量噪声的方差。
(S6)、采用卡尔曼滤波对改进的精密单点定位的观测方程中的参数进行估计,得到接收机位置信息和对流层延迟、电离层延迟信息,实现GNSS快速单点精密定位。
在PPP数据处理中,通常采用卡尔曼(Kalman)滤波对参数进行估计,状态方程表达为
Xk+1=Φk+1,kXk+wk,wk~N(0,Qk) (12)
式(12)中,Xk表示状态向量,Φk+1,k为状态转移矩阵,wk为动态噪声,Qk为动态噪声wk的协方差阵。
动态噪声协方差阵Qk反映了各待估参数的时变特性。在本发明方法中,主要参数包括接收机三维位置、接收机钟差、对流层天顶湿延迟、视线方向电离层延迟和载波相位模糊度,其中接收机钟差、对流层天顶湿延迟、视线方向电离层延迟模型化为随机游走,载波相位模糊度在参数估计时当作常数处理。于是,各待估参数动态噪声如下
式(13)中,qp为接收机三维位置谱密度,在静态和动态定位中分别为单位阵和零矩阵,为接收机钟差谱密度;qtrop为对流层天顶湿延迟的谱密度,qion为视线方向电离层延迟的谱密度,大气误差谱密度的选取与其变化率有关,Δt为相邻历元的时间间隔。
应用实施例
测试站点选取全球范围分布的144个IGS连续运行观测站,测站分布如图2,图2为本发明的应用实施例采用的IGS观测站的分布图。选取2013年6月17日00时至24时的观测数据,采样间隔为30s,间隔为15min的精密星历和间隔为5min的钟差产品,卫星截至高度角设置为5°。对144个IGS站24小时的观测数据以每0.5小时进行截取,得6912组数据。采用IGS发布的各个观测站坐标作为参考值,当水平位置偏差从某一历元开始均小于10cm时,认为滤波收敛,实现厘米级定位,而从第一个历元到该历元的时间间隔为收敛时间。为分析本发明提供PPP的定位方法的效果,处理步骤为:
步骤(1):对应用实施例中的6912组数据分别使用本发明的方法、UC方法、UofC方法、UD方法进行计算,对平面位置收敛到厘米级精度的百分比和平均收敛时间进行分析;得到本发明的方法、UC方法、UofC方法和UD方法的收敛百分比和收敛时间统计信息,参见表1。
从表1中,可以得到:本发明方法具有最高的滤波收敛比例和最少的收敛时间,分别为70.2%和12.9min,较UC方法、UofC方法和UD方法收敛比例提高了19.2、2.8和49.3个百分点,平均收敛时间分别缩短了18.4%、7.2%和29.5%。可见,本发明方法能够较高比例的收敛到厘米级,且收敛时间较短。
表1收敛百分比和收敛时间统计信息
步骤(2):选取本发明方法与UC方法公共滤波收敛部分,使用本发明方法与UC方法对收敛时间和平面位置偏差方面进行分析,得到收敛时间统计结果和平面位置偏差统计结果,参见图3和图4,图3为使用本发明方法和UC方法对公共滤波收敛数据处理后的收敛时间统计结果图;图4为使用本发明方法和UC方法对公共滤波收敛数据处理后的平面位置偏差统计结果图。从本发明方法与UC方法共同收敛部分(3007组数据)的收敛时间和平面位置偏差统计结果可以得到:本发明方法和UC方法的平均收敛时间分别为12.3min和13.6min,较UC方法缩短了9.6%,且分别有50%的收敛时间小于11.0min和15.5min;本发明方法和UC方法估计得到的接收机平面位置的平均偏差分别为3.7cm和4.3cm,且分别有75.4%和65.0%的平面位置偏差位于(0,5)cm内。
步骤(3):选取本发明方法与UofC方法公共滤波收敛部分,使用本发明方法与UofC方法对收敛时间和平面位置偏差方面进行分析,得到收敛时间统计结果和平面位置偏差统计结果,参见图5和图6,图5为使用本发明方法和UofC方法对公共滤波收敛数据处理后的收敛时间统计结果图;图6为使用本发明方法和UofC方法对公共滤波收敛数据处理后的平面位置偏差统计结果图。从本发明方法与UofC方法共同收敛部分(4090组数据)的收敛时间和平面位置偏差统计结果可以得到:本发明方法和UofC方法的平均收敛时间分别为12.1min和15.7min,较UofC方法缩短了22.9%,且分别有50%的收敛时间小于11.5min和13.0min;本发明方法和UofC方法估计得到的接收机平面位置的平均偏差分别为3.7cm和3.9cm,且分别有73.9%和71.5%的平面位置偏差位于(0,5)cm内。
步骤(4):选取本发明方法与UD方法公共滤波收敛部分,使用本发明方法与UD方法对收敛时间和平面位置偏差方面进行分析,得到收敛时间统计结果和平面位置偏差统计结果,参见图7和图8,图7为使用本发明方法和UD方法对公共滤波收敛数据处理后的收敛时间统计结果图;图8为使用本发明方法和UD方法对公共滤波收敛数据处理后的平面位置偏差统计结果图。从本发明方法与UD方法共同收敛部分(1180组数据)的收敛时间和平面位置偏差统计结果可以得到:本发明方法和UD方法的平均收敛时间分别为11.9min和18.4min,较UD方法缩短了35.3%,且分别有50%的收敛时间小于11.0min和19.0min;本发明方法和UD方法估计得到的接收机平面位置的平均偏差分别为3.7cm和5.2cm,且分别有74.7%和48.2%的平面位置偏差位于(0,5)cm内。
通过本应用实施例的数据处理,得到:使用本发明的一种快速收敛的精密单点定位方法能够快速、准确、高比例的估计接收机位置信息,相对于目前存在的UC方法、UofC方法和UD方法在定位精度和收敛速度方面均有明显提高。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同技术特征的含义和范围内的所有变化囊括在本发明内。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (2)
1.一种快速收敛的精密单点定位方法,其特征在于,包括如下步骤:
(S1)、对全球导航卫星的观测数据进行预处理:
采用M-W组合法和电离层残差法对原始载波相位观测值中存在的周跳进行探测,并修复;
对天线相位中心偏差、地球自转、相对论效应、天线相位缠绕误差源的影响进行修正;
对精密卫星轨道和卫星钟差数据进行拉格朗日多项式内插,获取观测时刻的卫星轨道和卫星钟差;
(S2)、顾及电离层延迟对不同频率观测值的影响关系,建立卫星的原始载波相位观测值和的观测方程,如公式(5)、公式(6)所示:
<mrow>
<msubsup>
<mi>&Phi;</mi>
<mn>1</mn>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<mi>&rho;</mi>
<mo>+</mo>
<mi>c</mi>
<mo>&CenterDot;</mo>
<mi>&delta;</mi>
<msub>
<mover>
<mi>t</mi>
<mo>&OverBar;</mo>
</mover>
<mi>r</mi>
</msub>
<mo>-</mo>
<mi>c</mi>
<mo>&CenterDot;</mo>
<mi>&delta;</mi>
<msup>
<mover>
<mi>t</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
</msup>
<mo>+</mo>
<msubsup>
<mi>M</mi>
<mi>d</mi>
<mi>s</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>p</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msubsup>
<mi>M</mi>
<mi>w</mi>
<mi>s</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>p</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>I</mi>
<mn>1</mn>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mover>
<mi>N</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<msubsup>
<mi>&Phi;</mi>
<mn>1</mn>
<mi>s</mi>
</msubsup>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>&Phi;</mi>
<mn>2</mn>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<mi>&rho;</mi>
<mo>+</mo>
<mi>c</mi>
<mo>&CenterDot;</mo>
<mi>&delta;</mi>
<msub>
<mover>
<mi>t</mi>
<mo>&OverBar;</mo>
</mover>
<mi>r</mi>
</msub>
<mo>-</mo>
<mi>c</mi>
<mo>&CenterDot;</mo>
<mi>&delta;</mi>
<msup>
<mover>
<mi>t</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
</msup>
<mo>+</mo>
<msubsup>
<mi>M</mi>
<mi>d</mi>
<mi>s</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>p</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msubsup>
<mi>M</mi>
<mi>w</mi>
<mi>s</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>p</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>f</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>/</mo>
<msubsup>
<mi>f</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>I</mi>
<mn>1</mn>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mover>
<mi>N</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<msubsup>
<mi>&Phi;</mi>
<mn>2</mn>
<mi>s</mi>
</msubsup>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,s表示卫星,和表示载波相位观测值,ρ表示卫星与接收机之间的几何距离,c表示真空中光的速度,表示接收机钟差,表示卫星钟差,分别表示对流层干、湿迟延映射函数;δzpd、δzpw分别表示对流层天顶方向干延迟和湿延迟,为载波相位观测值中站星视线方向的电离层延迟,分别为原始载波相位观测值和的模糊度,f1、f2为载波相位的频率,λ1、λ2表示载波相位的波长,表示载波相位观测值的观测噪声;
(S3)、根据原始码伪距观测值和原始载波相位观测值构建码相位半合组合观测值的观测方程,如公式(7)所示:
<mrow>
<msubsup>
<mi>L</mi>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&Phi;</mi>
<mn>2</mn>
</msub>
</mrow>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mn>2</mn>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&Phi;</mi>
<mn>2</mn>
<mi>s</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mo>=</mo>
<mi>&rho;</mi>
<mo>+</mo>
<mi>c</mi>
<mo>&CenterDot;</mo>
<mi>&delta;</mi>
<msub>
<mover>
<mi>t</mi>
<mo>&OverBar;</mo>
</mover>
<mi>r</mi>
</msub>
<mo>-</mo>
<mi>c</mi>
<mo>&CenterDot;</mo>
<mi>&delta;</mi>
<msup>
<mover>
<mi>t</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
</msup>
<mo>+</mo>
<msubsup>
<mi>M</mi>
<mi>d</mi>
<mi>s</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>p</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msubsup>
<mi>M</mi>
<mi>w</mi>
<mi>s</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>p</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mover>
<mi>N</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
<mi>s</mi>
</msubsup>
<mo>/</mo>
<mn>2</mn>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<msubsup>
<mi>L</mi>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&Phi;</mi>
<mn>2</mn>
</msub>
</mrow>
<mi>s</mi>
</msubsup>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
(S4)、根据原始载波相位观测值和构建几何无关组合观测值的观测方程,如公式(8)所示:
<mrow>
<msubsup>
<mi>&Phi;</mi>
<mrow>
<mi>G</mi>
<mi>F</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&Phi;</mi>
<mn>1</mn>
<mi>s</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&Phi;</mi>
<mn>2</mn>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>f</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>/</mo>
<msubsup>
<mi>f</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msubsup>
<mi>I</mi>
<mn>1</mn>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mover>
<mi>N</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
<mi>s</mi>
</msubsup>
<mo>-</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mover>
<mi>N</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<msubsup>
<mi>&Phi;</mi>
<mrow>
<mi>G</mi>
<mi>F</mi>
</mrow>
<mi>s</mi>
</msubsup>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
(S5)、联合原始载波相位观测值和的观测方程、码相位半合组合观测值的观测方程和几何无关组合观测值的观测方程构建改进的精密单点定位的观测方程,如公式(9)所示,
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&Phi;</mi>
<mn>1</mn>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<mi>&rho;</mi>
<mo>+</mo>
<mi>c</mi>
<mo>&CenterDot;</mo>
<mi>&delta;</mi>
<msub>
<mover>
<mi>t</mi>
<mo>&OverBar;</mo>
</mover>
<mi>r</mi>
</msub>
<mo>-</mo>
<mi>c</mi>
<mo>&CenterDot;</mo>
<mi>&delta;</mi>
<msup>
<mover>
<mi>t</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
</msup>
<mo>+</mo>
<msubsup>
<mi>M</mi>
<mi>d</mi>
<mi>s</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>p</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msubsup>
<mi>M</mi>
<mi>w</mi>
<mi>s</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>p</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>I</mi>
<mn>1</mn>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mover>
<mi>N</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<msubsup>
<mi>&Phi;</mi>
<mn>2</mn>
<mi>s</mi>
</msubsup>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&Phi;</mi>
<mn>2</mn>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<mi>&rho;</mi>
<mo>+</mo>
<mi>c</mi>
<mo>&CenterDot;</mo>
<mi>&delta;</mi>
<msub>
<mover>
<mi>t</mi>
<mo>&OverBar;</mo>
</mover>
<mi>r</mi>
</msub>
<mo>-</mo>
<mi>c</mi>
<mo>&CenterDot;</mo>
<mi>&delta;</mi>
<msup>
<mover>
<mi>t</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
</msup>
<mo>+</mo>
<msubsup>
<mi>M</mi>
<mi>d</mi>
<mi>s</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>p</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msubsup>
<mi>M</mi>
<mi>w</mi>
<mi>s</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>p</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>f</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>/</mo>
<msubsup>
<mi>f</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>I</mi>
<mn>1</mn>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mover>
<mi>N</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<msubsup>
<mi>&Phi;</mi>
<mn>2</mn>
<mi>s</mi>
</msubsup>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>L</mi>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&Phi;</mi>
<mn>2</mn>
</msub>
</mrow>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mn>2</mn>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&Phi;</mi>
<mn>2</mn>
<mi>s</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mo>=</mo>
<mi>&rho;</mi>
<mo>+</mo>
<mi>c</mi>
<mo>&CenterDot;</mo>
<mi>&delta;</mi>
<msub>
<mover>
<mi>t</mi>
<mo>&OverBar;</mo>
</mover>
<mi>r</mi>
</msub>
<mo>-</mo>
<mi>c</mi>
<mo>&CenterDot;</mo>
<mi>&delta;</mi>
<msup>
<mover>
<mi>t</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
</msup>
<mo>+</mo>
<msubsup>
<mi>M</mi>
<mi>d</mi>
<mi>s</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>p</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msubsup>
<mi>M</mi>
<mi>w</mi>
<mi>s</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>p</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mover>
<mi>N</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
<mi>s</mi>
</msubsup>
<mo>/</mo>
<mn>2</mn>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<msubsup>
<mi>L</mi>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&Phi;</mi>
<mn>2</mn>
</msub>
</mrow>
<mi>s</mi>
</msubsup>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&Phi;</mi>
<mrow>
<mi>G</mi>
<mi>F</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&Phi;</mi>
<mn>1</mn>
<mi>s</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&Phi;</mi>
<mn>2</mn>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>f</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>/</mo>
<msubsup>
<mi>f</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msubsup>
<mi>I</mi>
<mn>1</mn>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mover>
<mi>N</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
<mi>s</mi>
</msubsup>
<mo>-</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mover>
<mi>N</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<msubsup>
<mi>&Phi;</mi>
<mrow>
<mi>G</mi>
<mi>F</mi>
</mrow>
<mi>s</mi>
</msubsup>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
(S6)、采用卡尔曼滤波对改进的精密单点定位的观测方程中的参数进行估计,得到接收机位置信息和对流层延迟、电离层延迟信息;
对步骤(S5)的改进的精密单点定位的观测方程进行线性化处理,得到公式(10):
y(k)=A·X(k)+εy,εy~N(0,Ωy) (10)
其中,y(k)为观测向量,其维数为4j,j为历元k同步观测到的卫星数量;X(k)为待估参数向量,维数为(5+3j),A为设计矩阵,εy为测量噪声向量,Ωy为测量噪声向量εy的协方差阵;y(k)、X(k)和A表达式分别为
<mrow>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>&Phi;</mi>
<mn>1</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msubsup>
<mi>&Phi;</mi>
<mn>1</mn>
<mi>j</mi>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>&Phi;</mi>
<mn>2</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msubsup>
<mi>&Phi;</mi>
<mn>2</mn>
<mi>j</mi>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>L</mi>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&Phi;</mi>
<mn>2</mn>
</msub>
</mrow>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msubsup>
<mi>L</mi>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&Phi;</mi>
<mn>2</mn>
</msub>
</mrow>
<mi>j</mi>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>&Phi;</mi>
<mrow>
<mi>G</mi>
<mi>F</mi>
</mrow>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msubsup>
<mi>&Phi;</mi>
<mrow>
<mi>G</mi>
<mi>F</mi>
</mrow>
<mi>j</mi>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
</mrow>
<mrow>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
<mtd>
<mi>Y</mi>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
<mtd>
<mrow>
<mi>c</mi>
<mo>&CenterDot;</mo>
<mi>&delta;</mi>
<mover>
<msub>
<mi>t</mi>
<mi>r</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
<mtd>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>z</mi>
<mi>p</mi>
<mi>w</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msubsup>
<mi>I</mi>
<mn>1</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msubsup>
<mi>I</mi>
<mn>1</mn>
<mi>j</mi>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mover>
<mi>N</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msubsup>
<mover>
<mi>N</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
<mi>j</mi>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mover>
<mi>N</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msubsup>
<mover>
<mi>N</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
<mi>j</mi>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
</mrow>
<mrow>
<mi>A</mi>
<mo>=</mo>
<mo>&lsqb;</mo>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>e</mi>
<mn>4</mn>
</msub>
<mo>&CircleTimes;</mo>
<mi>B</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>C</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>I</mi>
<mi>j</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>D</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>I</mi>
<mi>j</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>&rsqb;</mo>
</mrow>
其中,X、Y、Z为接收机三维位置,表示克罗内克积,e4为各元素均为1的4维列向量,Ij为j维单位矩阵,B、C和D子矩阵的含义分别为:
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<mi>l</mi>
<mn>1</mn>
</msup>
</mtd>
<mtd>
<msup>
<mi>m</mi>
<mn>1</mn>
</msup>
</mtd>
<mtd>
<msup>
<mi>n</mi>
<mn>1</mn>
</msup>
</mtd>
<mtd>
<msubsup>
<mi>M</mi>
<mi>w</mi>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>l</mi>
<mi>j</mi>
</msup>
</mtd>
<mtd>
<msup>
<mi>m</mi>
<mi>j</mi>
</msup>
</mtd>
<mtd>
<msup>
<mi>n</mi>
<mi>j</mi>
</msup>
</mtd>
<mtd>
<msubsup>
<mi>M</mi>
<mi>w</mi>
<mi>j</mi>
</msubsup>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>C</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mo>-</mo>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>-</mo>
<msubsup>
<mi>f</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>/</mo>
<msubsup>
<mi>f</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>f</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>/</mo>
<msubsup>
<mi>f</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>D</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中,l、m、n分别为接收机到卫星的方向余弦。
2.根据权利要求1所述的一种快速收敛的精密单点定位方法,其特征在于:步骤(S5)得到的改进的精密单点定位的观测方程第k历元对应的协方差阵Ωy(k),如公式(11)所示:
<mrow>
<msub>
<mi>&Omega;</mi>
<mi>y</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>F</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>I</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
其中, 表示克罗内克积,Ij为j维单位矩阵,分别为卫星原始载波相位观测值Φ1、Φ2,原始码伪距观测值P2的测量噪声的方差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510660882.9A CN105182388B (zh) | 2015-10-10 | 2015-10-10 | 一种快速收敛的精密单点定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510660882.9A CN105182388B (zh) | 2015-10-10 | 2015-10-10 | 一种快速收敛的精密单点定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105182388A CN105182388A (zh) | 2015-12-23 |
CN105182388B true CN105182388B (zh) | 2017-08-25 |
Family
ID=54904584
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510660882.9A Active CN105182388B (zh) | 2015-10-10 | 2015-10-10 | 一种快速收敛的精密单点定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105182388B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108508470A (zh) * | 2018-05-16 | 2018-09-07 | 武汉大学 | 面向全球无缝的瞬时分米级导航定位方法 |
CN109655849A (zh) * | 2018-12-24 | 2019-04-19 | 中国电子科技集团公司第二十研究所 | 一种ppp快速定位收敛的方法 |
CN109633690B (zh) * | 2018-12-25 | 2020-11-17 | 中国电子科技集团公司第二十研究所 | 一种电离层梯度参数的确定方法、装置及系统 |
CN109861780B (zh) * | 2019-03-13 | 2020-06-23 | 浙江工业大学 | 一种wams时延信号补偿装置 |
CN111007556B (zh) * | 2019-12-10 | 2021-08-17 | 武汉大学 | 一种顾及方向约束信息的gps/bds单点测速方法 |
CN111273327B (zh) * | 2020-03-20 | 2022-04-08 | 中国人民解放军61081部队 | 一种基于组合和非组合混合观测模型的精密单点定位方法 |
CN114518586B (zh) * | 2021-03-17 | 2024-04-30 | 山东科技大学 | 一种基于球谐展开的gnss精密单点定位方法 |
CN113791434B (zh) * | 2021-09-14 | 2023-08-01 | 中国科学院国家授时中心 | 一种gnss对流层湿延迟估计方法、系统、设备及介质 |
CN114966760B (zh) * | 2022-06-15 | 2023-08-22 | 中国科学院精密测量科学与技术创新研究院 | 一种电离层加权的非差非组合ppp-rtk技术实现方法 |
CN115166799B (zh) * | 2022-07-11 | 2023-07-21 | 中国科学院精密测量科学与技术创新研究院 | 顾及硬件延迟时变特征的gnss精密单点定位方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1875291A (zh) * | 2003-10-28 | 2006-12-06 | 天宝导航有限公司 | 三个或更多载波的gnss信号的模糊度估计 |
CN102353969A (zh) * | 2011-09-02 | 2012-02-15 | 东南大学 | 精密单点定位技术中相位偏差的估计方法 |
-
2015
- 2015-10-10 CN CN201510660882.9A patent/CN105182388B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1875291A (zh) * | 2003-10-28 | 2006-12-06 | 天宝导航有限公司 | 三个或更多载波的gnss信号的模糊度估计 |
CN102353969A (zh) * | 2011-09-02 | 2012-02-15 | 东南大学 | 精密单点定位技术中相位偏差的估计方法 |
Non-Patent Citations (3)
Title |
---|
Precise Point Positioning Using IGS Orbit and Clock Products;Jan Kouba等;《GPS Solutions》;20011031;第5卷(第2期);全文 * |
基于GPS双频原始观测值的精密单点定位算法及应用;张宝成等;《测绘学报》;20101031;第39卷(第5期);第478-483页 * |
基于多导航卫星系统的精密单点定位模型与方法研究;徐宗秋;《中国博士学位论文全文数据库 基础科学辑》;20150215;正文第22-28、52-62页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105182388A (zh) | 2015-12-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105182388B (zh) | 一种快速收敛的精密单点定位方法 | |
CN111190203B (zh) | 基于非差观测模型和消秩亏的ppp-rtk定位方法 | |
EP3130943B1 (en) | Navigation satellite system positioning involving the generation of tropospheric correction information | |
EP1550241B1 (en) | Position estimation using a network of global-positioning receivers | |
CN105807300B (zh) | 一种用北斗双频接收机进行高精度动态单点定位的方法 | |
CN103760572B (zh) | 一种基于区域cors的单频ppp电离层加权方法 | |
CN104102822B (zh) | 一种多频gnss观测值随机特性建模方法 | |
CN105891860B (zh) | 一种基于误差分离模式的gnss区域伪距差分增强定位方法 | |
CN105044741B (zh) | 一种伪距相位综合广域差分改正值的求解方法 | |
WO2008008099A2 (en) | Gnss signal processing methods and apparatus | |
CN102096084A (zh) | 基于星间组合差分的精密单点定位方法 | |
CN110161547B (zh) | 一种自适应电离层估计模型的中长基线模糊度解算方法 | |
CN107193023B (zh) | 一种具有闭式解的高精度北斗卫星系统单点定位方法 | |
Liu et al. | Comparison of convergence time and positioning accuracy among BDS, GPS and BDS/GPS precise point positioning with ambiguity resolution | |
CN114966760B (zh) | 一种电离层加权的非差非组合ppp-rtk技术实现方法 | |
CN112526564A (zh) | 一种精密单点定位重新收敛方法 | |
Parvazi et al. | Role of the RLS-VCE-estimated stochastic model for improvement of accuracy and convergence time in multi-GNSS precise point positioning | |
Choy | High accuracy precise point positioning using a single frequency GPS receiver | |
Tolman et al. | Absolute precise kinematic positioning with GPS and GLONASS | |
Zhu et al. | Comparative analysis of four different single-frequency PPP models on positioning performance and atmosphere delay retrieval | |
Liu et al. | Combining GPS+ GLONASS observations to improve the fixing percentage and precision of long baselines with limited data | |
Basile et al. | Analysis on the potential performance of GPS and Galileo Precise Point Positioning using simulated Real-Time Products | |
Martin | GNSS precise point positioning: The enhancement with GLONASS | |
CN104502943B (zh) | 一种室内伪卫星差分相对定位方法 | |
Li et al. | A grid-based ionospheric weighted method for PPP-RTK with diverse network scales and ionospheric activity levels |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |