GPS整周模糊度的求解方法解析
- 格式:doc
- 大小:1.25 MB
- 文档页数:18
实时GPS姿态测量中整周模糊度的快速解算方法彭晓刚;吕志平;王新山【期刊名称】《海洋测绘》【年(卷),期】2011(31)2【摘要】姿态测量核心问题是整周模糊度解算.提出了一种适合实时姿态测量的模糊度解算方法,利用单差平滑伪距进行解算,与传统的模糊度解算方法相比具有许多优点.通过仿真实验验证了该方法的有效性和实用性.%The most important problem of the GPS real-time attitude determination is integer ambiguity resolution. This article discusses a new algorithm for integer ambiguity resolution which can be used in real-time attitude determination that uses single-difference smoothed pseudo range measurements. Over the conventional algorithms, it has several advantages. And we investigate the validity and the performance of the new algorithm through simulations.【总页数】4页(P34-37)【作者】彭晓刚;吕志平;王新山【作者单位】解放军信息工程大学测绘学院,河南郑州450052;91404部队,河北秦皇岛066000;解放军信息工程大学测绘学院,河南郑州450052;91404部队,河北秦皇岛066000【正文语种】中文【中图分类】P228.43【相关文献】1.GPS姿态测量的载波相位整周模糊度快速解算 [J], 郑庆晖;张育林2.GPS载体姿态测量中的整周模糊度的快速求解 [J], 肖迪;王美玲3.单频GPS动态定位中整周模糊度的一种快速解算方法 [J], 刘宁;熊永良;冯威;徐韶光4.GPS姿态测量中的整周模糊度的实时解算方法 [J], 彭晓刚;杨靖宝;吝宁5.GPS接收机RTK定位中整周模糊度的快速解算方法 [J], 冯仲科;韩熙春因版权原因,仅展示原文概要,查看原文内容请购买。
GPS整周模糊度GPS整周模糊度的计算与确定引⾔精密型GPS信号接收机⼀般都具有伪距和载波相位两种基本观测量。
相对于伪噪声码观测量⽽⾔,GPS载波相位观测量能提供⾮常精确的相对定位。
但由于GPS载波相位测量存在整周模糊数较难解算的问题,致使它在快速定位及导航中的应⽤受到了限制。
因此,快速⽽准确地求解GPS载波相位测量的整周模糊度就成了它在快速定位及导航中应⽤的关键问题。
整周模糊度求解的理论及其实⽤研究是近⼀、⼆⼗年的研究热点和难点。
许多学者提出了⼀些解算⽅法,其中双频P码伪距法、整周模糊度函数法、最⼩⼆乘搜索法和整周模糊度协⽅差法应⽤较⼴泛。
整周模糊度的确定是GPS载波相位测量中的关键问题,其原因如下:精确地、不⾜⼀周的相位与修复周跳后的正确整周记数只有在与正确的整周模糊度配合使⽤才有意义。
整周模糊度参数⼀旦出现问题,就将导致⼤量的卫地距出现系统性的粗差,从⽽严重影响定位的精度和可靠性,正确确定整周模糊度N是获得⾼精度定位结果的必要条件。
在⼤量对精确确定整周模糊度的计算研究中不断推出了新的计算算法。
⼏种整周模糊度的确定⽅法:(⼀)快速求解整周模糊度伪距双差⽅程经过线性化之后如下[2],(1)其中,ρ表⽰实际观测值与计算值之差,A表⽰系数阵,δx表⽰坐标增量,v表⽰模型误差和测量噪声,N(·)表⽰正态分布,QDΨ表⽰伪距测量的协⽅差阵。
由式(1),根据最⼩⼆乘原理可得(2)对于载波相位,其双差模型线性化之后可得[3](3)其中,l表⽰实际观测值与计算值之差,λ表⽰L1载波波长,N表⽰载波相位双差模糊度,w 表⽰模型误差和测量噪声,QDφ表⽰载波相位测量的协⽅差阵。
由式(2)、(3),可得整周模糊度的浮点解N^。
(4)由式(4)根据协因数传播定律,此时整周模糊度N^的协⽅差阵QN^为(5)其中表⽰坐标增量的协⽅差阵;表⽰后验⽅差系数;表⽰残差;n表⽰卫星数;u= rank(A)表⽰系数阵A的秩。
第3卷第3期2006年6月CHIN ESE J OU RNAL OF EN GIN EERIN G GEOP H YSICSVol 13,No 13J un 1,2006文章编号:1672—7940(2006)03—0225—05GPS 整周模糊度解算的L AMBDA 法及程序实现林 丹,郭 敏,江 华,蒋旭惠(中国地质大学工程学院,武汉 430074)作者简介:林丹(1982—),男,湖北武汉人,中国地质大学(武汉)在读研究生,主要研究方向:GPS ,GIS 技术的工程应用。
E 2mail :lind56001@摘 要:对目前GPS 模糊度解算方法中搜索效率和成功率较好的L AMBDA 的理论进行深入探讨的同时,结合实例对L AMBDA 的解算流程进行了研究。
利用C #编程语言实现了卫星位置、导航位置、基线向量的浮点解和固定解,以及整周模糊度等未知量的可视化输出,为GPS 其它研究工作提供了研究基础。
关键词:GPS ;整周模糊度解算;L AMBDA ;整数最小二乘中图分类号:P228文献标识码:A收稿日期:2006—03—23GPS LAMB DA METH OD AN D ITS PR OGRAM REALIZATIONL IN Dan ,GUO Min ,J IAN G Hua ,J IAN G Xu 2hui(Facult y of Engineering ,China Universit y of Geosciences ,W uhan 430074,China )Abstract :In view of t he shortage of detailed programmed realization of met hods of GPS ambi 2guity resolution ,t his paper talks about L AMBDA met hod for it s high efficiency and success ,and st udies t he resolution p rogram of LAMBDA met hod wit h an example.Then C #p ro 2gram language is used to realize t he floating point resolution and t he fixed resolution of t he following factors such as t he satellite position ,navigation position and t he baseline vector.In addition ,t he outp ut visualization of t he integer ambiguity resolution is achieved.K ey w ords :GPS ;integer ambiguity resolution ;LAMBDA ;integer least square1 引 言过去的二十多年中,国内外许多学者对整周模糊度解算的理论进行了研究,提出了许多解算整周模糊度的方法,如模糊度函数法、最小二乘法和最小二乘模糊度去相关法。
整周模糊度的确定确定整周未知数,是基于载波相位测量进行相对定位,必需解决的另一个关键问题。
精确和快速地求解整周未知数,对于确保相对定位的高精度,提高作业效率,开拓高精度动态定位新方法,都是极其重要的。
确定整周未知数的方法许多,若按解算所需时间的长短区分,可分为经典静态相对定位法和快速解算模糊度(整周未知数)法,而快速解算模糊度法又包括交换天线法,P码双频法、滤波法,搜寻法和模糊函数法等等;若按确定整周未知数时gps接收机的运动状态区分,又可分为静态法和动态法。
上述各种快速解算法皆属于静态法的范畴。
所谓动态法,就是GPS接收机在运动状态中完成求解整周未知数,它是实施高精度实时动态定位的基础。
一、经典静态相对定位法确定整周未知数这种方法是将作为待定的未知参数,在基线平差中与其它未知参数(如δXi、δYi、δZi等)一并求解的方法。
一般是由载波相位观测值组成双差分观测方程式,并进行方程式线性化,得到双差分误差方程式,则该方程式中包含有待定测站三个坐标改正数δXi、δYi、δZi和整周未知数的线性组合这四个未知数[此处]。
只要在已知测站和待定测站上同步观测不少于4颗卫星,则可平差解出整周未知数。
用这种方法一般需观测较长时间(几非常钟至几小时),但解算的精度最高,常用于静态相对定位中,尤其是用于长距离相对定位中。
在平差计算中,依据对的取值方式不同,可分为“整数解”(固定解)和“实数解”(浮动解)两种。
整数解是利用应当是整数的特性[也应为整数],将解得的▽▽N(t0)值进行凑整(凑成最接近的整数),然后将凑整后的作为已知量再代入双差分误差方程,重新平差,解算待定测站坐标改正数。
这种方法,只有当观测误差和外界误差对观测值影响较小,解得的比较接近整数的状况下才有效,此时,它可以提高解算结果的精度。
整数解常用于四、五十公里以下的基线的相对定位。
实数解当联测基线较长时,某些外界误差(如大气折射误差、卫星星历误差等)对基线两端点观测值的影响差别较大(即相关性不强),这时,在两测站间求差分时,就不能较好地消退或减弱其影响,它们在基线平差解算中将被汲取进待定测站坐标改正数和整周未知数中,这样解算出来的整周未知数一般偏离整数值较远,且其精度较低,误差可能大于半周,这时,我们不再考虑的整数特性,而取其实际解算值―实数解。
GPS 整周模糊度解算方法探讨一、为什么要解算GPS 整周模糊度?整周模糊度的确定是载波相位测量中的关键问题,这是因为:(1)精确的、不足1周的相位观测值()φr F 和修复周跳后的正确的整周计数()φInt 只有与正确的整周模糊度配合使用才有意义。
模糊度参数一旦出错,就将导致大量的卫地距出现系统性的粗差,从而严重损害定位的精度和可靠性。
正确确定整周模糊度N 是获得高精度定位结果的必要条件。
(2)在一般精度的GPS 定位中,定位所需的时间实际上就是正确确定整周模糊度所需要的时间。
快速确定整周模糊度对提高GPS 定位的作业效率具有极其重要的作用;对开拓GPS 定位技术的应用领域,将其推广应用到低等级控制测量和一般的工程测量等领域也具有极其重要的作用。
二、GPS 整周模糊度解算方法1、LAMBDA 法1993年荷兰Delft 大学的Teunissen 教授提出了最小二乘模糊度降相关平差法,简称LAMBDA 法。
该方法可缩小搜索范围,加快搜索过程,是目前快速静态定位中最成功的一种模糊度搜索方法。
LAMBDA 法的基本原理: (1)整数变换在LAMBDA 法中,并不直接对整数模糊度参数N 进行搜索,而是先对初始解中的实数模糊度参数⎪⎭⎫⎝⎛=∧∧∧∧n N N N N ,......,,21及其协因数阵∧N Q 进行整数变换:∧∧⋅=N Z z TZ Q Z Q NT z⋅⋅=∧∧式中Z 为整数变换矩阵。
整数变换具有以下特点:当N 为整数时,变换后的参数z 也为整数;反之,当z 为整数时,经逆变换后所得的()z Z N T⋅=-1也为整数。
整数变换并不是唯一的。
我们希望整数变换后所得到的新参数⎪⎭⎫⎝⎛=∧∧∧∧n z z z z ,......,,21之间的相关性能显著减小,其协因数阵∧z Q 中的非对角线元素5.0≤,模糊度参数的方差也能大幅度减小。
注意,整数变换指的是具有上述特性的一种数学变换方法,但并非只能对整数进行变换。
GPS精密定位周跳检测与修复(Cycle slip detection and repair)完整的载波相位是由初始整周模糊度N、计数器记录的整周数INT和接收机基频信号与收到到卫星信号的小于一周部分相位差Δφ。
Δφ能以极高的精度测定,但这只有在N和INT都正确无误地确定情况下才有意义。
卫星在观测中失锁后,造成接收机载波整周计数INT误差,这种现象称为周跳。
当重新捕获卫星后,周跳给计数器造成的偏差即为中断期间丢失的整周数,小周跳可以通过检测方法发现后并加以修复,大的周跳或较长时间的失锁,周跳不易修复,需要重新固定整周模糊度。
周跳的探测及修复对于用载波相位精密定位至关重要,成功的修复才能获得高精度的结果。
周跳产生的原因:1.卫星信号暂时阻断;2.仪器线路暂时故障;3.外界环境的突变干扰,如电离层、动态变化。
检测周跳的主要方法:1.屏幕扫描法观测值中出现周跳后。
相位观测值的变化率就不再连续。
凡曲线出现不规则的突然变化时,就意味着在相应的相位观测值中出现了整周跳变。
早期进行GPS相位测量的数据处理时,就是靠作业人员坐在计算机屏幕前依次对每个站、每个时段、每个卫星的相位观测值的变化率的图像进行逐段检查来探测周跳,然后再加以修复。
这种方法比较直观,在早期曾广泛使用。
但由于工作繁琐枯燥乏味,而且需反复进行,所以这种手工编辑方法目前正逐步被淘汰,而很少使用了。
2.高次差或多项式拟合法由于卫星和接收机间的距离在不断变化,因而载波相位测量的观测值INT+Δφ也随时间在不断变化。
但这种变化应是有规律的、平滑的。
周跳将破坏这种规律性。
根据这一特性就能将一些大的周跳寻找出来(尤其是对采样率较高的数据)。
一般来说,一个测站S对同一卫星J的相位观测量,对不同历元间相位观测值取至4至5次差之后,距离变化对整周数的影响已可忽略,这时的差值主要是由于振荡器的随机误差而引起的,因而应具有随机的特性见下表。
但是,如果在观测过程中产生了周跳现象,那么便破坏了上述相位观测量的正常变化规率,从而使其高次差的随机特性也受到破坏。
胡国辉孟浩袁信摘要:对Cholesky分解整周模糊度的求解进行了改进,在求解整周模糊度的过程中,首先采用LAMBDA法对整周模糊度进行整数线性变换再作Cholesky分解,然后利用最优剪枝法(best cut)对整周模糊度进行搜索,实验结果表明该方法具有快速搜索整周模糊度的能力,可以满足采用GPS载波相位测量确定姿态以及GPS载波相位测量与INS组合的实时性。
关键词:导航整周模糊度载波相位Cholesky分解中图分类号:V241.5FAST CARRIER PHASE AMBIGUITY RESOLUTION FORDIFFERENCE GPSHu Guohui1, Meng Hao2, Yuan Xin11(Department of Automatic Control, Naijing University of Aeronautics & Astronautics,Nanjing,210016)2(Department of Automatic Control, Harbin EngineeringUniversity,Harbin,150001)Abstract The paper presents a new development method for Cholesky ambiguity search method. The method makes use of an ambiguity reparametrization, Cholesky decomposition and best cut. Experiment results show that the method can achieve fast search ability, and satisfy real time attitude determination and GPS/INS integration with GPS carrier phase measurement.Key words navigation, ambiguity, carrier phase, Cholesky factorization单纯采用Cholesky分解整周模糊度的求解[1]往往搜索次数较多,采用LAMBDA法[2]对整周模糊度进行整数线性变换再作Cholesky分解,使变换后的整周模糊度方差更小,有效的提高了搜索速度,实验结果表明该算法能快速确定整周模糊度,能满足采用载波相位的姿态确定以及与惯导组合着陆的实时性要求。