西安测绘招聘信息（论文推荐）

《测绘学报》

构建与学术的桥梁 拉近与权威的距离

复制链接，关注《测绘学报》抖音！

【测绘学报的个人主页】长按复制此条消息，长按复制打开抖音查看TA的更多作品##7NsBSynuc88##[抖音口令]

本文内容来源于《测绘学报》2020年第8期，审图号GS（2020）4062号。

GNSS模糊度降相关性能的条件方差平稳度评价法

卢立果^1,2,刘万科³,鲁铁定²,马立烨⁴,吴汤婷²,杨元喜¹

1.西安测绘研究所地理信息工程国家重点实验室, 陕西 西安 710054;

2.东华理工大学测绘工程学院, 江西 南昌 330013;

3.武汉大学测绘学院, 湖北 武汉 430079;

4.武汉大学卫星导航定位技术研究中心, 湖北 武汉 430079

收稿日期：2019-10-12；修回日期：2020-05-26

基金项目：国家自然科学基金（41804020；41774031）；国家重点研发计划（2016YFB0501405）；江西省自然科学基金（20202BAB212010；20192BAB217011）

摘要：GNSS模糊度降相关通过整数变换优化条件方差的排列顺序，提高搜索效率。降相关和条件方差的关系及其评价是关键问题之一。针对这一问题，本文从理论上分析了排序后模糊度降相关与条件方差之间的数值关系，发现降相关性能与条件方差数值序列的平稳性有关，降相关性能越强，条件方差数值序列越平稳。基于这一理论关系，给出了“条件方差平稳度”定义，并将其作为评价降相关性能的指标。通过模拟和实测数据验证，并采用条件方差变化趋势图和搜索时间来定性和定量评价降相关性能，用以判定条件方差平稳度的合理性。试验结果表明，条件方差平稳度可以较精确直观地衡量模糊度的降相关性能。本文定义的指标揭示了模糊度降相关的本质。

关键词：GNSS 模糊度 降相关 条件方差 评价指标

引文格式：卢立果, 刘万科, 鲁铁定, 等. GNSS模糊度降相关性能的条件方差平稳度评价法. 测绘学报，2020，49(8)：955-964. DOI: 10.11947/j.AGCS.2020.20190417.

阅读全文：http://xb.sinomaps.com/article/2020/1001-1595/2020-8-955.htm

全文概述

载波相位整周模糊度的快速、准确解算是GNSS实时高精度定位的关键。在众多整周模糊度解算方法中，基于整数最小二乘搜索的估计方法成功率最高，而被广泛采用^[1-7]。为提高模糊度的搜索效率，文献[5]提出了基于整数变换的LAMBDA降相关算法。随后，许多学者也基于降相关思想开展了相关研究^[8-11]。通常认为，模糊度降相关准则是通过降低模糊度分量之间的相关性，实现搜索椭球的压缩，以便提高搜索效率。基于这一准则，降相关系数和条件数常被用作评价降相关性能的指标^[8-21]。近来，文献[22-24]通过理论和试验对比发现，条件数和降相关系数并不能准确反映模糊度的降相关性能，模糊度降相关的目的在于通过整数变换，优化条件方差的排序结构，使其尽可能按照降序或升序的方式排序。

尽管降相关性能与排序程度直接相关，但是不同降相关算法排序策略的差异性，致使无法直接判断不同算法性能的优劣^{[5, 8-10]}。为合理评价降相关算法的性能，文献[25-26]通过绘制条件方差的变化趋势图，来直接判读不同算法的性能。但这种方式仅限于单个模糊度方差的比较分析，对降相关性能差异较明显的算法管用；文献[27-30]直接采用搜索时间来衡量不同算法的降相关性能，但是需要经过模糊度搜索环节才能判定算法的优劣，尤其当模糊度精度较高时，不同算法间的搜索时间由于解算耗时量级较小^[29]，容易受计算机运行误差的影响，不利于算法间性能的合理比较。此外，搜索时间的长短不是模糊度固定的全部性能。

由于评价指标的选取对于正确评价降相关算法的性能至关重要，且便于理解降相关的本质，因此，需要定义合理的评价指标来衡量降相关性能。下文在阐述降相关基本理论基础上，重点从理论上分析排序与模糊度条件方差的数值关系，并定义相应指标，进而通过试验验证定义的指标的合理性。

1 降相关原理

假定^[5]

(1)

式中，a为模糊度整数候选向量。

通常人们采用整数变换来快速获取模糊度整数向量^{[5, 8-11]}。其核心思想是通过采用整数变换矩阵Z对

1.1 整数高斯消元

对L^TDL)分解

(2)

式中，L=[li,j]，i>j是单位下三角矩阵；D=[d_i]是对角阵，d_i为模糊度分量的条件方差。

当li,j> 0.5时，按照下式进行高斯消元

(3)

式中，μ=[li,j]_int表示对li,j取整。

1.2 条件方差排序

当_i,di+1)进行交换，其交换矩阵为

(4)

式中，

分别为i-1阶和n-i-1阶的单位阵。

采用Pi,i+1对L^TDL分解

(5)

式中

(6)

式中

(7)

2 条件方差平稳度评价指标

不同降相关算法由于采用的整数变换策略不同，导致算法性能各异，需要对其进行有效评价，而当前常用的条件数和降相关系数两个评价指标已不能全面地衡量降相关算法的性能^[22-24]。由于评价指标的选取对正确评价降相关算法的性能至关重要，为合理定义评价指标，下文首先分析条件方差的排序结构是影响降相关性能的关键，然后从理论上分析排序后模糊度降相关与条件方差之间的数值关系，最后在此基础上定义合理的评价指标。

不同于文献[22-24]采用的降相关的条件方差的排序分析法，本文将采用二维图形变换的方式对整数高斯消元和条件方差排序的整数变换过程进行分析，剖析二者之间的关系。

给定模糊度方差-协方差矩阵χ²=0.15。图 1为模糊度整数变换前后椭圆形状变换，表 1为搜索节点数的变化(先对模糊度协方差矩阵进行L^TDL分解，再将模糊度从后到前搜索，即先搜索N₂然后再进行N₁的搜索)。结合图表，可以看出，通过整数高斯消元法降低了分量间的相关性，实现对椭球的压缩，但不会改变搜索节点数；条件方差排序(可以看作矩阵的正交变换)并不改变椭球的大小，而是旋转了椭球方向，减少N₂的节点个数。因此，排序是决定搜索效率的关键因素。根据判别条件

图 1 整数变换前后椭圆形状变化

Fig. 1 Elliptic shape change before and after integer transformation

图选项

表 1 整数变换前后搜索节点数

Tab. 1 Number of search nodes before and after integer transformation

表选项

尽管降相关是通过整数变换优化条件方差的排列顺序，但是，并不能确保严格依据条件方差的降序排序(针对L^TDL分解)，即并非完全满足

(8)

具体交换关系见图 2。

图 2 相邻条件方差交换

Fig. 2 Diagram of adjacent conditional variance exchange

图选项

根据整数变换不改变矩阵的行列式这一性质，条件方差交换后满足等式

(9)

不等式(9)证明如下

(10)

由式(8)—式(9)可知，条件方差交换的实质是减少条件方差数值间的差异。降相关性能越强，条件方差交换程度越大，交换的条件方差数值间差异越小，获得的条件方差数值序列累积越小，进而条件方差的整体数值平稳性越好。于是，首先定义“条件方差平稳度”。

定义：条件方差平稳度定义为

(11)

式中，S即为“条件方差平稳度”，以此作为衡量模糊度向量排序的合理度。

性质：S满足不等式

(12)

当d₁=d₂=…=d_n时，不等式取等号，因此S的取值范围为

(13)

上述性质可根据

由式(9)—式(12)可以看到，S值取决于条件方差的平稳程度，条件方差数值序列波动越小，解算的S值越小，当条件方差相等时不等式取等号。因此，S值越小，条件方差平稳度越好，此时降相关质量越好，模糊度搜索效率越高；反之亦然。

3 试验数据分析

为验证定义的条件方差平稳度是否可以作为衡量模糊度降相关性能的有效评价指标，采用3种排序形式不同的降相关算法来验证降相关性能的差异性，分别为基于自然升序排序的ASCE算法^[10]、基于下三角Cholesky(LDL^T)分解的LIGT算法^[30]和基于最小列旋转的SEQR算法^[10]。为保证3种算法降相关性能比较的合理性，3种算法均采用对模糊度协方差矩阵进行LDL^T分解的升序排序算法，其中ASCE是对LDL^T分解的相邻条件方差交换算法，SEQR是将相邻条件方差交换条件拓展到整个模糊度向量。因此理论上3种算法的降相关性能为SEQR > LIGT > ASCE。需要说明的是，整数高斯消元步骤导致整数变换过程需要循环迭代，因此实际解算过程中3种算法并不是完全服从上述不等式关系。

本文采用条件方差变化趋势图^[25-26]和搜索耗时^[27-30]来定性和定量评价降相关性能，用以判定条件方差平稳度的合理性。其中，搜索算法采用当前最为流行的SE-VB搜索算法^{[27, 31-33]}。

本文的模糊度搜索耗时是在PC机上基于Matlab 2016a软件进行统计的，其软硬件配置为：Intel Core i7-7700 CPU，2.80 GHz主频，8 GB内存，Win10操作系统(64位)。

3.1 仿真试验

借鉴文献[27]的模拟方法随机模拟模糊度的方差-协方差阵

(14)

式中，L=[li,j]，i>j是单位下三角矩阵；D=[d_i]是对角阵。考虑到实际模糊度前3个条件方差远大于其余模糊度的条件方差^[5]，条件方差d_i

(15)

式中，rand表示产生在(0, 1)之间均匀分布的随机数。

模糊度浮点解

(16)

式中，mvnrnd (·)表示多元正态分布函数；I_n为n阶单位阵。

鉴于实际解算的模糊度精度和维数是随历元变化的，为从普遍意义下验证条件方差平稳度评价指标的合理性，本文设计了模糊度不同维数和精度两个角度的仿真试验。

3.1.1 不同维数下仿真试验

试验1分别采用20、30、40、50维共4组数据作仿真试验，由于载波相位观测量越多，模糊度维数越大，浮点模糊度精度越高^{[29, 34]}，因此4组d₁分别设置为500/100/50/10。为保证搜索耗时的差异性和稳定性，搜索候选向量个数设置为5，每组数据均搜索5次，并统计其平均搜索耗时。

表 2给出了原始模糊度(origin)、ASCE、LIGT和SEQR 3种降相关算法的模糊度条件方差平稳度和搜索耗时，图 3为4种模糊度的条件方差变化趋势图。由图表可见，降相关算法获得的条件方差平稳度越小，搜索效率越高，条件方差变化趋势越平稳，整体波动越小。其中，在20维处ASCE的条件方差平稳度小于LIGT、50维处LIGT的条件方差平稳度小于SEQR，此时搜索耗时与条件方差平稳度均呈一一对应关系。

表 2 采用不同降相关算法的模糊度条件方差平稳度和搜索耗时

Tab. 2 Ambiguity conditional variance stationarity and search time under different decorrelation algorithms

表选项

图 3 采用不同降相关算法的模糊度条件方差变化趋势

Fig. 3 The trend of ambiguity conditional variance under different decorrelation algorithms

图选项

3.1.2 不同精度下仿真试验

试验2采用4组维数在10~50维变化的不同精度的d₁分别设置为1/5/10/25，搜索候选向量个数设置为5，每一维模糊度数据均搜索5次。由图 4可见，在d₁为1/5/10时，3种算法的条件方差平稳度和搜索耗时一一对应；在d₁为25时，模拟的

图 4 不同精度下模糊度条件方差和搜索时间变化趋势

Fig. 4 The trend of ambiguity conditional variance and search time under different precision

图选项

3.2 实测试验

选取两组不同观测环境下的动态车载实测数据(试验基本情况见表 3)，并同时统计不同降相关算法的条件数和降相关系数，通过分析比较这两种评价指标，从不同维度对降相关性能评价。其中，通常认为降相关后获得的条件数数值越小、降相关系数数值越大，则降相关算法的降相关性能越强^[8-21]。双差模糊度的浮点解

表 3 试验基本情况

Tab. 3 Basic information of the experiment

表选项

图 5是模糊度维数和ADOP随历元变化的趋势图。由图 5可见，模糊度维数变化范围为25~46，ADOP数值维持在0.10以下。ADOP值取决于两个卫星系统的模糊度维数、站星几何关系和载波相位观测值精度等因素，ADOP并不严格对应模糊度维数。由于ADOP值越小，模糊度精度较高，搜索耗时数量级相对较小^[29]。基于此，本文为突出3种算法的时间差异性，模糊度搜索的候选向量个数均设置为20，每组数据搜索5次，并统计每一历元的平均搜索耗时。

图 5 不同历元下的模糊度维数与ADOP值

Fig. 5 Ambiguity dimensions and the value of ADOP under different epochs

图选项

图 6给出了ASCE、LIGT和SEQR 3种算法在不同历元下的解算结果。由图 6可见，在条件方差平稳度方面，SEQR具有最小的条件方差平稳度，LIGT存在部分历元的条件方差平稳度大于ASCE；在搜索耗时方面，考虑到计算机运行误差，LIGT存在个别历元略小于SEQR的现象，整体上3种算法的搜索耗时趋势基本和条件方差平稳度一致。对于条件数和降相关系数两个指标的评价效果，可以通过直接对比ASCE和SEQR两种算法来获得相关结论。根据图 6，ASCE相对于SEQR整体上具有相对较小的条件数，与文献[10]结论一致；ASCE相对于SEQR整体上具有较大的降相关系数。

图 6 不同历元下模糊度评价指标和搜索时间的关系

Fig. 6 The relationship between ambiguity evaluation index and search time under different epochs

图选项

图 7 降相关系数比值(ASCE/SEQR)的概率密度函数

Fig. 7 Probability density functions of the ratio of decorrelation coefficients for ASCE relative to SEQR

图选项

图 7中，P_better表示采用降相关系数统计的ASCE降相关性能优于SEQR的概率；P_worse表示ASCE劣于SEQR的概率。单纯从这两方面比较来看ASCE在降相关性能上优于SEQR。但是，ASCE的搜索效率明显低于SEQR，说明这两种指标在评价降相关质量方面具有一定的片面性，与文献[22-24]结论一致。因此，实际数据结果表明，相对于常用的条件数和降相关系数评价指标，“条件方差平稳度”更能够准确可靠衡量

4 结论

本文从理论上分析了排序过程对模糊度条件方差的影响，在此基础上，定义了评定降相关性能的“条件方差平稳度”指标。采用模拟和实测数据对其效果进行了验证，可以得到以下两点结论：

(1) 降相关通常采用整数变换，优化条件方差的排序结构，促使条件方差尽量按照升序或降序排序，为此需要对条件方差进行交换。由于两个条件方差交换后将会缩小二者数值间的差异，因此当完成整个降相关过程后，条件方差序列的整体数值差异相对较小，数值序列趋势较平稳。此时降相关性能增强，进而提高搜索效率。

(2) 由于整数变换前后矩阵行列式大小不变，因此变换后的条件方差序列数值之和小于交换前数值之和。根据这一性质，本文定义的条件方差平稳度评价指标更能可靠地反映模糊度条件方差的变化趋势，并能够较直观地衡量模糊度的降相关性能。当获得的平稳度数值越小且趋近于1时，条件方差的趋势变化越平稳，此时降相关性能越好。

理论与实际计算表明，以整数变换优化条件方差的排序结构实现的降相关，可有效缩小条件方差数值间差异，获得一组更平稳的数值序列，进而提高模糊度的搜索效率。

作者简介

第一作者简介：卢立果（1989-），男，博士，讲师，研究方向为GNSS模糊度解算。E-mail: lglu66@163.com

第二作者简介：（通讯作者）刘万科（1978-），男，博士，教授，研究方向为GNSS精密导航定位、GNSS/INS/视觉等多源信息协同精密定位、导航卫星自主定轨。E-mail: wkliu@sgg.whu.edu.cn

第三作者简介：鲁铁定（1974-），男，博士，教授，研究方向为误差理论与数据处理、GNSS数据处理。E-mail: tdlu@ecut.edu.cn

第四作者简介：马立烨（1993-），男，博士生，研究方向为GNSS精密定位。E-mail: liyema@whu.edu.cn

第五作者简介：吴汤婷（1990-），女，博士，讲师，研究方向为卫星大地测量学。E-mail: ttwu@ecut.edu.cn

第六作者简介：杨元喜(1956-)，男，博士，研究员，中国科学院院士，研究方向为动态大地测量数据与卫星导航数据处理。E-mail: yuanxi_yang@163.com

《测绘学报（英文版）》（JGGS）专刊征稿：LiDAR数据处理

论文推荐 | 胡敏章, 张胜军,金涛勇,文汉江,褚永海,姜卫平,李建成:新一代全球海底地形模型BAT_WHU2020

人物 | 无止境的科学追求——记大地测量学家陈俊勇院士

招聘启事丨事业编制！海南热带海洋学院2020年公开招聘工作人员计划表公布

权威 | 专业 | 学术 | 前沿

欢迎加入《测绘学报》作者QQ群： 751717395

进群请备注：姓名+单位+稿件编号