• 中国科学引文数据库(CSCD)
  • 中文科技期刊数据库
  • 中国核心期刊(遴选)数据库
  • 日本科学技术振兴机构数据库(JST)
  • 中国学术期刊(网络版)(CNKI)
  • 中国学术期刊综合评价数据库(CAJCED)
  • 中国超星期刊域出版平台

广义阻尼LAMBDA方法及其在长基线实时定位中的应用

曹士龙, 蔚保国, 伍蔡伦, 贾浩男

曹士龙, 蔚保国, 伍蔡伦, 贾浩男. 广义阻尼LAMBDA方法及其在长基线实时定位中的应用[J]. 全球定位系统, 2023, 48(4): 69-74. DOI: 10.12265/j.gnss.2023043
引用本文: 曹士龙, 蔚保国, 伍蔡伦, 贾浩男. 广义阻尼LAMBDA方法及其在长基线实时定位中的应用[J]. 全球定位系统, 2023, 48(4): 69-74. DOI: 10.12265/j.gnss.2023043
CAO Shilong, YU Baoguo, WU Cailun, JIA Haonan. Generalized damped LAMBDA method and the application in long baseline real-time positioning[J]. GNSS World of China, 2023, 48(4): 69-74. DOI: 10.12265/j.gnss.2023043
Citation: CAO Shilong, YU Baoguo, WU Cailun, JIA Haonan. Generalized damped LAMBDA method and the application in long baseline real-time positioning[J]. GNSS World of China, 2023, 48(4): 69-74. DOI: 10.12265/j.gnss.2023043

广义阻尼LAMBDA方法及其在长基线实时定位中的应用

基金项目: 石家庄市科技计划项目(226130417A)
详细信息
    作者简介:

    曹士龙: (1986—),男,博士,工程师,研究方向为GNSS精密定位技术及相关嵌入式软件开发

    蔚保国: (1966—),男,博士,研究员,研究方向为卫星导航系统体制、地面运控、终端芯片、测试评估以及位置服务系统

    伍蔡伦: (1981—),男,博士,研究员,研究方向为卫星导航与信号处理

    贾浩男: (1990—),男,硕士,工程师,研究方向为室内外无缝定位技术及智能导航硬件开发

    通信作者:

    曹士龙 E-mail: caoshilong2006@163.com

  • 中图分类号: P228.4

Generalized damped LAMBDA method and the application in long baseline real-time positioning

  • 摘要: 在长基线实时定位中,相邻观测值之间具有强相关性导致法方程出现病态,此时载波相位整周模糊度难以快速固定. 针对这一问题,在传统阻尼最小二乘模糊度降相关平差(least-squares ambiguity decorrelation adjustment, LAMBDA)方法基础上,提出了广义阻尼LAMBDA方法,并推导出载波相位模糊度的单历元解表达式. 该方法利用了坐标和载波相位模糊度的先验信息,在改善法方程病态性的同时也提高了浮点模糊度的解算精度,有助于整周模糊度的快速固定. 采用两条长基线(266 km和456 km)实测数据进行验证,改进后的方法相比于传统方法,载波相位模糊度固定率提高了28.8%;当放大因子c=0时,广义阻尼LAMBDA方法与传统阻尼LAMBDA方法等价. 选择澳大利亚连续运行参考站系统(continuous operation reference system,CORS)网中四条长基线实时数据流进行测试,实验结果表明:在基线长度小于1 000 km的多系统实时定位中可以获得水平优于0.02 m,高程优于0.04 m定位精度. 基于广义阻尼LAMBDA方法在长基线实时定位中具有一定的应用参考价值.
    Abstract: In long baseline real-time positioning, strong correlation between adjacent observations resulting in ill-conditioned normal equation which is adverse to fast integer solution of carrier phase ambiguity. In this contribution, a generalized damped least-squares ambiguity decorrelation adjustment (LAMBDA) method was be proposed based on traditional LAMBDA method and the solution expression for carrier phase ambiguity in single epoch also be derived. This method uses the prior information of coordinates and carrier phase ambiguity to improve the ill-condition of the normal equation and the solution accuracy of floating ambiguity, which is helpful for the fast integer ambiguity solution. Two long baselines (266 km and 456 km) measured data were used to verify. Compared with the traditional method, the ambiguity fixed rate of the improved method was increased by 28.8%. When the amplification factor c=0, the generalized damped LAMBDA method was equivalent to the traditional damped LAMBDA method. Four long baseline real-time data streams from the Australian continuous operation reference system (CORS) network were selected for testing. The experimental results show that multi-system real-time positioning accuracy was better than 2 cm in the horizontal direction and 4 cm in the elevation direction with a baseline length of less than 1 000 km. The generalized damped LAMBDA method has certain application reference value for long baseline real-time positioning.
  • GNSS实时相对定位技术能够为用户提供精确的三维坐标、速度和时间等信息,已被广泛应用于精密测量、精细农业、气候变化监测和地质灾害预警等诸多领域[1-3]. 整周模糊度解算是实现GNSS精密相对定位的关键. 如何提高定位精度和整周模糊度的收敛速度,国内外学者进行了深入而广泛的研究[4-6]. 在基于模糊度域搜索方法中,最小二乘模糊度降相关平差(least-squares ambiguity decorrelation adjustment, LAMBDA)法是公认搜索最快、理论最为严密的方法[7].

    在长基线实时定位中,短时间内的观测量之间存在相关性,会导致法方程病态,当引入对流层延迟参数后,进一步加剧了法方程的病态性. 另外,随着基线长度的增加,与空间相关的残余误差影响也逐渐增大,通常无法实现单历元模糊度固定,需要一定的收敛时间[8-10]. 如何缩短模糊度的收敛时间,提高模糊度固定的正确率,仍然值得进一步研究. 金星等[11]针对单频单历元组合载波相位差分技术定位过程中存在的秩亏及病态问题,提出将伪距观测值引入,采用经验分权法定权的一种新的模糊度降相关方法. 易重海等[12]对历元间坐标差虚拟观测误差方程进行重构,提出一种改进的历元间坐标差方法,削弱了法方程的病态性,提高了模糊度固定的稳定性和效率. 当观测数据质量不佳时,采用经典LAMBDA法计算的模糊度固定解仍然可能出现错误值. 为此一些学者提出了部分模糊度固定算法[13-15],只求解具有最大成功率的模糊度固定解,从而提高参数的估计精度. 高旺等[16]提出一种基于部分固定策略的多系统组合的长基线模糊度快速解算方法,实现了长距离基准站间模糊度快速固定. 刘根友等[17]在传统LAMBDA方法的基础上结合参数约束平差,提出了阻尼LAMBDA方法,改善了法方程的病态性,提高了模糊度解算效率. 本文借鉴该算法的思路,将部分模糊度固定与阻尼LAMBDA方法结合起来,提出了一种改进的LAMBDA方法,实现了阻尼LAMBDA方法和部分模糊度固定算法的统一. 在此之前,简要介绍一下阻尼LAMBDA方法.

    在最优化方法中,为了解决法方程的病态问题,通过适当加大矩阵主对角元素改善法方程的条件数,降低法方程的病态性. 阻尼LAMBDA方法将坐标先验的权阵作为阻尼因子添加到载波相位观测方程的法方程中,既改善了法方程的病态性,也实现了单历元固定载波相位模糊度.

    阻尼LAMBDA方法的观测方程为

    $$ \left\{ \begin{array}{*{20}{l}} {{\boldsymbol{V}}_{\boldsymbol{X}}} = {\boldsymbol{X}},& {{\boldsymbol{P}}_{\boldsymbol{X}}} \\ {{\boldsymbol{V}}_\varphi } = {\boldsymbol{AX}} + {\boldsymbol{BY}} - {L_\varphi },& {\boldsymbol{P}} \end{array} \right. . $$ (1)

    式中:$ {{\boldsymbol{V}}_{\boldsymbol{X}}} $为坐标先验值的误差项;$ {{\boldsymbol{V}}_\varphi } $为双差载波相位观测值的残差项;$ {\boldsymbol{X}} $为非模糊度参数向量,$ {\boldsymbol{Y}} $为模糊度参数向量,$ {\boldsymbol{A}} $$ {\boldsymbol{B}} $分别为相应的系数矩阵;$ {{\boldsymbol{L}}_\varphi } $为双差载波相位观测值的常数项;$ {{\boldsymbol{P}}_{\boldsymbol{X}}} $为坐标参数的先验权阵,也称为阻尼矩阵;$ {\boldsymbol{P}} $为载波相位观测值先验权阵.

    阻尼LAMBDA方法的法方程为

    $$ \left[ {\begin{array}{*{20}{c}} {{{\boldsymbol{A}}^{\rm{T}}}{\boldsymbol{PA}} + {{\boldsymbol{P_X}}}}&{{{\boldsymbol{A}}^{\rm{T}}}{\boldsymbol{PB}}} \\ {{{\boldsymbol{B}}^{\rm{T}}}{\boldsymbol{PA}}}&{{{\boldsymbol{B}}^{\rm{T}}}{\boldsymbol{PB}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\boldsymbol{X}} \\ {\boldsymbol{Y}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{{\boldsymbol{A}}^{\rm{T}}}{\boldsymbol{P}}{{\boldsymbol{L}}_\varphi }} \\ {{{\boldsymbol{B}}^{\rm{T}}}{\boldsymbol{P}}{{\boldsymbol{L}}_\varphi }} \end{array}} \right] $$ (2)

    由于添加了阻尼矩阵,该法方程变为满秩方程. 这里的先验坐标值一般通过伪距单点定位或者伪距差分定位计算得到,先验权阵$ {{\boldsymbol{P}}_{\boldsymbol{X}}} $的形式为

    $$ {{\boldsymbol{P}}_{\boldsymbol{X}}} = {\boldsymbol{Q}}_{{\boldsymbol{XX}}}^{ - 1} \approx {{\boldsymbol{A}}^{\rm{T}}}{{\boldsymbol{P}}_\rho }{\boldsymbol{A}} $$ (3)

    式中:A为观测方程的系数阵;$ {{\boldsymbol{P}}_\rho } $为伪距观测量权阵,与载波相位观测量的权阵$ {\boldsymbol{P}} $关系为

    $$ {{\boldsymbol{P}}_\rho } = \alpha {\boldsymbol{P}}. $$ (4)

    式中,$ \alpha $为伪距与载波相位观测量的权比值,需要根据客观的先验精度来确定.

    阻尼LAMBDA方法本质上是附有坐标约束条件的整数最小二乘法. 通过对坐标施加约束条件,改善法方程的病态性. 单频观测数据采用单历元阻尼LAMBDA方法定位时,一般只有在0.5 m精度以内的坐标约束才可能获得比较可靠的结果,而双频接收机可以适当放宽到1 m[17]. 文献[17]中详细对比分析了不同精度的初始坐标约束对模糊度解算的影响,本文不再赘述.

    而在实际解算过程中,随着参与解算的历元数增加,模糊度的浮点解精度会逐渐提高,最后收敛于整数值. 因此,为了充分利用已有坐标和模糊度的先验信息,本文提出了一种改进的阻尼LAMBDA方法,也称为广义阻尼LAMBDA方法.

    广义阻尼LAMBDA方法的观测方程为

    $$ \left\{ \begin{array}{*{20}{l}} {{\boldsymbol{V}}_{\boldsymbol{X}}} = {\boldsymbol{X}}, &{{\boldsymbol{P}}_{\boldsymbol{X}}} \\ {{\boldsymbol{V}}_{\boldsymbol{B}}} = {\boldsymbol{Y}}, &{{\boldsymbol{P}}_{\boldsymbol{B}}} \\ {{\boldsymbol{V}}_\varphi } = {\boldsymbol{A}}X + {\boldsymbol{B}}{\boldsymbol{Y}} - {{{L}}_\varphi },&{\boldsymbol{P}} \end{array} \right. $$ (5)

    平差准则为

    $$ {\boldsymbol{V}}_\varphi ^{\rm{T}} {\boldsymbol{P}}{{\boldsymbol{V}}_\varphi } + {\boldsymbol{V}}_{\boldsymbol{X}}^{\rm{T}}{{\boldsymbol{P}}_{\boldsymbol{X}}}{{\boldsymbol{V}}_{\boldsymbol{X}}} + {\boldsymbol{V}}_{\boldsymbol{B}}^{\rm{T}}{{\boldsymbol{P}}_{\boldsymbol{B}}}{{\boldsymbol{V}}_{\boldsymbol{B}}} = \min $$ (6)

    广义阻尼LAMBDA方法的法方程为

    $$ \left[ {\begin{array}{*{20}{c}} {{{\boldsymbol{A}}^{\rm{T}}}{\boldsymbol{PA}} + {{\boldsymbol{P}}_{\boldsymbol{X}}}}&{{{\boldsymbol{A}}^{\rm{T}}}{\boldsymbol{PB}}} \\ {{{\boldsymbol{B}}^{\rm{T}}}{\boldsymbol{PA}}}&{{{\boldsymbol{B}}^{\rm{T}}}{\boldsymbol{PB}} + {{\boldsymbol{P}}_{\boldsymbol{B}}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\boldsymbol{ X }}\\ {\boldsymbol{ Y}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{{\boldsymbol{A}}^{\rm{T}}}{\boldsymbol{P}}{{\boldsymbol{L}}_\varphi }} \\ {{{\boldsymbol{B}}^{\rm{T}}}{\boldsymbol{P}}{{\boldsymbol{L}}_\varphi }} \end{array}} \right]$$ (7)

    简记为

    $$ \left[ {\begin{array}{*{20}{c}} {{{\boldsymbol{N}}_{aa}} + {{\boldsymbol{P}}_{\boldsymbol{X}}}}&{{{\boldsymbol{N}}_{ab}}} \\ {{{\boldsymbol{N}}_{ba}}}&{{{\boldsymbol{N}}_{bb}} + {{\boldsymbol{P}}_{\boldsymbol{B}}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\boldsymbol{X}} \\ {\boldsymbol{Y}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{{\boldsymbol{U}}_a}} \\ {{{\boldsymbol{U}}_b}} \end{array}} \right] $$ (8)

    进一步简记为

    $$ {\boldsymbol{Q}}= \left[ {\begin{array}{*{20}{c}} {{{\boldsymbol{N}}_{aa}} + {{\boldsymbol{P}}_{\boldsymbol{X}}}}&{{{\boldsymbol{N}}_{ab}}} \\ {{{\boldsymbol{N}}_{ba}}}&{{{\boldsymbol{N}}_{bb}} + {{\boldsymbol{P}}_{\boldsymbol{B}}}} \end{array}} \right]=\left[ {\begin{array}{*{20}{c}} {\boldsymbol{Q}}_{\boldsymbol{\hat X \hat X}} & {\boldsymbol{Q}}_{\boldsymbol{\hat X \hat Y}} \\ {\boldsymbol{Q}}_{\boldsymbol{\hat Y \hat X}} & {\boldsymbol{Q}}_{\boldsymbol{\hat Y \hat Y}} \end{array}} \right] $$ (9)

    这里的$ {{\boldsymbol{P}}_{\boldsymbol{X}}} $为坐标参数的先验权矩阵,$ {{\boldsymbol{P}}_{\boldsymbol{B}}} $为模糊度参数的先验权矩阵. 对坐标和模糊度参数取较大权时,称为强约束;当取较小权时,称为弱约束. 因此在广义阻尼LAMBDA方法中,除了对坐标参数给予先验权阵外,对所有模糊度参数也给予先验权阵. 初始模糊度参数的先验方差矩阵具有如下形式:

    $$ {{\boldsymbol{P}}_{\boldsymbol{B}}} = {\text{diag}}\left(\begin{array}{*{20}{c}} {\displaystyle\frac{{\sigma _0^2}}{{\sigma _1^2}}}&{\displaystyle\frac{{\sigma _0^2}}{{\sigma _2^2}}}& \cdots &{\displaystyle\frac{{\sigma _0^2}}{{\sigma _m^2}}} \end{array}\right) $$ (10)

    式中:$ \sigma _0^2 $为载波相位单位权方差;$ \sigma _m^2 $为第m个模糊度的先验方差. 具体可采用式(11)计算:

    $$ \sigma _m^2 = \frac{{c{{{\text{(1}} - \left| \delta \right|)}^2}}}{{{{\sin }^2}(El)}} + 1 $$ (11)

    式中:$ \left| \delta \right| $为先验浮点模糊度参数的小数部分;El为卫星高度角;c为放大因子(通常取0或1). 初始模糊度先验方差计算公式不仅考虑了高度角的影响,而且考虑了浮点解靠近整数的程度. 特别指出,当c=0时,广义阻尼LAMBDA方法与传统阻尼LAMBDA方法等价;当部分模糊度参数的c=0时,广义阻尼LAMBDA方法等价于部分模糊度固定方法.

    广义阻尼LAMBDA方法不仅充分利用了先验坐标信息,而且充分利用了载波相位模糊度的先验信息,在改善法方程病态性的同时也提高了浮点模糊度的解算精度,有助于载波相位模糊度的快速固定.

    广义阻尼LAMBDA方法中浮点解$\hat {\boldsymbol{Y}}$的单历元解表达式为

    $$ \begin{aligned} \hat {\boldsymbol{Y}} =& {[({{\boldsymbol{N}}_{bb}} + {{\boldsymbol{P}}_{\boldsymbol{B}}}) - {\boldsymbol{N}}_{ab}^{\rm{T}}{({{\boldsymbol{N}}_{aa}} + {{\boldsymbol{P}}_{\boldsymbol{X}}})^{ - 1}}{{\boldsymbol{N}}_{ab}}]^{ - 1}}[{{\boldsymbol{U}}_b} \\&- {\boldsymbol{N}}_{ab}^{\rm{T}}{({{\boldsymbol{N}}_{aa}} + {{\boldsymbol{P}}_{\boldsymbol{X}}})^{ - 1}}{{\boldsymbol{U}}_a}] \end{aligned} $$ (12)

    浮点解$\hat {\boldsymbol{Y}}$的协因数阵为

    $$ {{\boldsymbol{Q}}_{\hat {\boldsymbol{Y}}\hat {\boldsymbol{Y}}}} = {[({{\boldsymbol{N}}_{bb}} + {{\boldsymbol{P}}_{\boldsymbol{B}}}) - {\boldsymbol{N}}_{ab}^{\rm{T}}{({{\boldsymbol{N}}_{aa}} + {{\boldsymbol{P}}_{\boldsymbol{X}}})^{ - 1}}{{\boldsymbol{N}}_{ab}}]^{ - 1}} $$ (13)

    浮点解$\hat {\boldsymbol{X}}$的表达式为

    $$ \hat {\boldsymbol{X }}= {({{\boldsymbol{N}}_{aa}} + {{\boldsymbol{P}}_{\boldsymbol{X}}})^{ - 1}}({{\boldsymbol{U}}_a} - {{\boldsymbol{N}}_{ab}}\hat {\boldsymbol{Y}}) $$ (14)

    模糊度固定后,按照下式更新浮点解$ \hat {\boldsymbol{X}} $,得到固定解$\breve{{\boldsymbol{X}}}$

    $$ \breve{{\boldsymbol{X}}} = \hat {\boldsymbol{X}} - {({{\boldsymbol{N}}_{aa}} + {{\boldsymbol{P}}_{\boldsymbol{X}}})^{ - 1}}{{\boldsymbol{N}}_{ab}}(\hat {\boldsymbol{Y}} - \breve{{\boldsymbol{Y}}} ) $$ (15)

    固定解$\breve{\boldsymbol{X}}$其协方差矩阵

    $$ \begin{aligned} {{\boldsymbol{Q}}_{\breve{{\boldsymbol{X}}} \breve{{\boldsymbol{X}}} }} =& {{\boldsymbol{Q}}_{\hat {\boldsymbol{X}}\hat {\boldsymbol{X}}}} - {{\boldsymbol{Q}}_{\hat {\boldsymbol{X}}\hat {\boldsymbol{Y}}}}{\boldsymbol{Q}}_{\hat {\boldsymbol{Y}}\hat {\boldsymbol{Y}}}^{ - 1}{{\boldsymbol{Q}}_{\hat {\boldsymbol{Y}}\hat {\boldsymbol{X}}}} \\ = &{({{\boldsymbol{N}}_{aa}} + {{\boldsymbol{P}}_{\boldsymbol{X}}})^{ - 1}} - {{\boldsymbol{Q}}_{\hat {\boldsymbol{X}}\hat {\boldsymbol{Y}}}}{({{\boldsymbol{N}}_{bb}} + {{\boldsymbol{P}}_{\boldsymbol{B}}})^{ - 1}}{{\boldsymbol{Q}}_{\hat {\boldsymbol{Y}}\hat {\boldsymbol{X}}}} \end{aligned} $$ (16)

    为了测试所提出的广义阻尼LAMBDA方法的性能,本节选取澳大利亚连续运行参考站系统(continuous operational reference system,CORS)网中的4个测站(STR1、PARK、 MOBS和 MCHL)的实测数据进行验证分析. 基线名称及长度如表1所示.

    表  1  基线名称及长度 km
    名称 测站 长度
    基线1 STR1-PARK 266
    基线2 STR1-MOBS 456
    基线3 MCHL-PARK 736
    基线4 MCHL-STR1 996
    下载: 导出CSV 
    | 显示表格

    本节选择2019年10月1日基线1和基线2的24 h实测数据进行事后仿动态处理. 对比测试方案如表2所示.

    表  2  广义阻尼LAMBDA测试方案
    方案 具体策略
    方案1 传统阻尼LAMBDA方法
    方案2 广义阻尼LAMBDA方法(全部c=0)
    方案3 广义阻尼LAMBDA方法(高度角<30,
    小数部分>0.25周时,c=0)
    下载: 导出CSV 
    | 显示表格

    以GPS和北斗卫星导航系统(BeiDou Navigation Satellita System,BDS)定位模式为例,分别采用上述三种方案求解整周模糊度,统计模糊度固定成功率. 对比结果如图1所示.

    图  1  不同方案的模糊度固定率统计对比

    图1可知,方案1的平均固定率为68.7%,方案2的平均固定率为68.9%,二者结果基本相同,这也验证了当c=0时广义阻尼LAMBDA方法与传统阻尼LAMBDA方法等价;方案3的平均固定率为88.6%,相比方案1,方案3的固定率提高了28.8%,说明广义阻尼LAMBDA方法可有效提高载波相位模糊度的固定率.

    为了测试广义阻尼LAMBDA方法在实时长基线定位中的定位精度,本节利用澳大利亚CORS网实时数据流进行测试. 测试时间为2021年3月29日四条长基线的实测数据. 其中状态空间表示(state space representation, SSR)产品采用法国空间研究中心(Centre National d’Etudes Spatiales, CNES)播发的实时流. 数据处理策略如表3所示.

    图2是统计的各测站的定位误差随时间变化的序列.

    将各基线在N、E、U三个方向误差的均方根(root mean square,RMS)进行统计对比,如图3所示.

    表  3  实时长基线解算策略
    项目处理策略
    定位模式GPS+Galileo+BDS
    观测频率GPS:L1+L2; Galileo:E1+E5b; BDS:B1+B2
    系统权比GPS∶Galileo∶BDS=3∶3∶2
    估计方法Kalman滤波
    高度角/采样率15°、 1 s
    待估参数3个坐标分量、双差模糊度(载波伪距权比:100∶1)、天顶对流层延迟(GMF投影函数+随机游走过程噪声1.0 e−8m2/s)
    电离层延迟消电离层组合
    卫星星历广播星历+ CNES的SSR产品(5 s)
    地球潮汐改正IERS conventions2010
    下载: 导出CSV 
    | 显示表格
    图  2  各个测站实时动态定位误差
    图  3  各个测站实时动态定位误差

    统计各个基线中的固定成功率和首次固定时间(判断标准为连续10个历元满足ratio值大于3),对比结果如图4所示.

    从统计结果看出,四条基线的水平定位误差平均值小于0.02 m,高程定位误差平均值小于0.04 m. 第一条基线的首次固定时间最短(10.5 min),随着基线长度增加,首次固定时间也逐渐增加. 其余三条基线的首次固定时间均超过15 min. 四条基线的平均首次固定时间为16.75 min,平均固定率为78.8%.

    图  4  首次固定时间(左)与固定率(右)对比

    本文针对长基线定位中由于法方程病态性导致模糊度收敛时间较长,模糊度固定率较低等问题,在传统阻尼LAMBDA方法基础上,结合部分模糊度固定算法,充分利用坐标和载波相位模糊度的先验信息,提出广义阻尼LAMBDA方法. 既实现了二者形式上的统一,又改善了方程的病态性,有助于载波相位模糊度的快速固定. 基于澳大利亚CORS站的长基线实测数据,将新方法与传统阻尼LAMBDA方法进行对比测试,结果表明,广义阻尼LAMBDA方法可有效提高模糊度的固定率. 采用四条长基线的实时数据流测试了新方法的定位精度. 结果表明,在基线长度小于1000 km的实时定位中,采用广义阻尼LAMBDA方法可获得水平方向优于0.02 m,高程方向优于0.04 m的定位精度,具有一定的应用参考价值.

  • 图  1   不同方案的模糊度固定率统计对比

    图  2   各个测站实时动态定位误差

    图  3   各个测站实时动态定位误差

    图  4   首次固定时间(左)与固定率(右)对比

    表  1   基线名称及长度 km

    名称 测站 长度
    基线1 STR1-PARK 266
    基线2 STR1-MOBS 456
    基线3 MCHL-PARK 736
    基线4 MCHL-STR1 996
    下载: 导出CSV

    表  2   广义阻尼LAMBDA测试方案

    方案 具体策略
    方案1 传统阻尼LAMBDA方法
    方案2 广义阻尼LAMBDA方法(全部c=0)
    方案3 广义阻尼LAMBDA方法(高度角<30,
    小数部分>0.25周时,c=0)
    下载: 导出CSV

    表  3   实时长基线解算策略

    项目处理策略
    定位模式GPS+Galileo+BDS
    观测频率GPS:L1+L2; Galileo:E1+E5b; BDS:B1+B2
    系统权比GPS∶Galileo∶BDS=3∶3∶2
    估计方法Kalman滤波
    高度角/采样率15°、 1 s
    待估参数3个坐标分量、双差模糊度(载波伪距权比:100∶1)、天顶对流层延迟(GMF投影函数+随机游走过程噪声1.0 e−8m2/s)
    电离层延迟消电离层组合
    卫星星历广播星历+ CNES的SSR产品(5 s)
    地球潮汐改正IERS conventions2010
    下载: 导出CSV
  • [1]

    ODOLINSKI R, TEUNISSEN P J G, ODIJK D. Combined BDS, Galileo, QZSS and GPS single-frequency RTK[J]. GPS solutions, 2015, 19(1): 151-163. DOI: 10.1007/s10291-014-0376-6

    [2]

    SHI C, ZHAO Q L, HU Z G, et al. Precise relative positioning using real tracking data from COMPASS GEO and IGSO satellites[J]. GPS solutions, 2013, 17(1): 103-119. DOI: 10.1007/s10291-012-0264-x

    [3] 王生亮, 刘根友, 高铭, 等. GPS RTK/UWB紧组合精密动态定位性能分析[J]. 全球定位系统, 2021, 46(2): 69-76.
    [4] 祝会忠, 李军, 蔚泽然, 等. 长距离GPS/BDS参考站网多频载波相位整周模糊度解算方法[J]. 测绘学报, 2020, 49(3): 300-311.
    [5]

    TEUNISSEN P J G, JONGE P D , TIBERIUS C C J M. The least-squares ambiguity decorrelation adjustment: its performance on short GPS baselines and short observation spans[J]. Journal of geodesy, 1997, 71(10): 589-602. DOI: 10.1007/s001900050127

    [6] 孙妍艳, 刘翠芝. 基于自适应差分进化算法的高维模糊度搜索[J]. 全球定位系统, 2018, 43(1): 36-42.
    [7]

    TEUNISSEN P J G. GPS carrier phase ambiguity fixing concepts[M]. GPS for Geodesy, Springer, Berlin, Heidelberg, 1998: 319-388.

    [8] 曹士龙, 刘根友, 王生亮, 等. GPS超长基线解算的误差特性与精度分析[J]. 武汉大学学报(信息科学版), 2023, 48(2): 260-267.
    [9] 张钰玺, 张小红, 刘全海, 等. 航空测量场景下的中长基线动态定位方法[J]. 测绘学报, 2019, 48(7): 871-878.
    [10] 于文浩, 张云, 韩彦岭, 等. 北斗三频中长基线差分定位性能研究[J]. 全球定位系统, 2018, 43(4): 102-109.
    [11] 金星, 王玲, 黄文德, 等. 一种单频单历元BDS/GPS组合整周模糊度解算方法[J]. 全球定位系统, 2018, 43(1): 75-80.
    [12] 易重海, 陈源军. 顾及历元间坐标差信息的GPS模糊度快速固定改进方法[J]. 武汉大学学报(信息科学版), 2019, 44(4): 489-494.
    [13]

    BRACK, A. Reliable GPS + BDS RTK positioning with partial ambiguity resolution[J]. GPS solutions, 2017, 21(1): 1083-1092. DOI: 10.1007/810291-016-0594-1

    [14] 舒宝, 刘晖, 张晋升, 等. 基于BDS/GPS组合定位的部分模糊度固定效果分析[J]. 武汉大学学报(信息科学版), 2017, 42(7): 989-994, 1001.
    [15]

    LI P, ZHANG X. Precise point positioning with partial ambiguity fixing[J]. Sensors, 2015, 15(6): 13627-13643. DOI: 10.3390/s150613627

    [16] 高旺, 高成发, 潘树国, 等. 基于部分固定策略的多系统长距离基准站间模糊度快速解算[J]. 武汉大学学报(信息科学版), 2017, 42(4): 558-562.
    [17] 刘根友, 朱耀仲, 韩保民. GPS单历元定位的阻尼LAMBDA算法[J]. 武汉大学学报(信息科学版), 2004, 29(3): 195-197.
  • 期刊类型引用(2)

    1. 孙胜浩,曹新亮,吴佳辉. 基于北斗的相对位移测量方法在滑坡监测中的设计与实现. 计算机测量与控制. 2024(12): 57-66 . 百度学术
    2. 阳仁贵,魏蒙恩. 一种提高模糊度解算成功率的新方法——优选拟合因子获取最大Ratio值. 全球定位系统. 2024(06): 1-7 . 本站查看

    其他类型引用(0)

图(4)  /  表(3)
计量
  • 文章访问数:  241
  • HTML全文浏览量:  83
  • PDF下载量:  30
  • 被引次数: 2
出版历程
  • 收稿日期:  2023-03-13
  • 网络出版日期:  2023-08-21
  • 刊出日期:  2023-09-14

目录

/

返回文章
返回
x 关闭 永久关闭