一种基于三维转一维的载波相位整周模糊度搜索的定位方法
【技术领域】
[0001] 本发明涉及卫星导航定位技术领域，尤其涉及一种基于三维转一维的载波相位整 周模糊度搜索的定位方法。
【背景技术】
[0002] 精度是卫星导航应用所追求的重要指标之一，对载波相位观测量的合理利用是取 得分米级甚至厘米级高精度定位的关键。对于整周模糊度求解总体上可分为四大类：双 频伪距法、模糊度函数法、最小二乘搜索法和模糊度协方差阵方法，分别从测量域、位置域、 模糊度域出发搜索整周模糊度。基于位置域搜索的模糊度函数法被认为是最早的0TF方法 (Hatch等，1994)。最早由Counselman等人于1981年提出，其后由Mader和Remondi将之 用于动态定位（Mader, 1992,Remondi, 1991) 〇
[0003] 传统的基于位置域的整周模糊度搜索方法是基于模糊度点位三维空间搜索，基于 三维位置域整周模糊度搜索的优点是对载波相位的整周不敏感，缺点是计算量大搜索速度 慢，在三维几何空间进行搜索，假设位置误差为dR，贝lj整个模糊度搜索空间是dRXdRXdR， 尽管后续算法有所改进但是因为计算量大原因近10年没有得到更进一步的突破发展。
【发明内容】
[0004] 有鉴于此，本发明提供了一种基于三维转一维的载波相位整周模糊度搜索的定位 方法，可以从很大程度上减少模糊度搜索计算量，实现整周模糊度位置域搜索方法的快速 确定。该方法可用于包括GPS、GLONASS、Compass、Galileo等在内多星座GNSS兼容系统或 单系统，实现整周模糊度的快速确定。
[0005] -种基于三维转一维的载波相位整周模糊度搜索的定位方法，具体过程如下：
[0006] 步骤1、在可视范围内，取高度角最大的一颗卫星为基准星k，在用户u、基准星k与 非基准星i组成的平面中，坐标系选取Z'轴垂直于该平面，垂直于基准星k与非基准星i 连线方向为X'轴，建立右手螺旋坐标系0' -X'Y'Z';然后将卫星坐标和用户初始概 略坐标从原始坐标系转换到坐标系〇' -X'Y'Z'下；
[0007] 步骤2、在坐标系0'-X'Y'Z'下，将双差线性化定位观测方程近似为：
[0008]
[0009] 其中，VAW是双差载波相位观测值，A是载波相位波长，Sy' u指用户位置在坐 标系-X'Y'Z'下Y'轴坐标的误差改正数，是整周模糊度，是用户到基 准星k和非基准星i的单位方向余弦矢量差在Y'轴的表示；VApf表示利用卫星坐标和用 户坐标计算的用户及基准站到卫星k和卫星i的双差距离，Noise表示误差；
[0010] 步骤3、在坐标系O' -X'Y'Z'下基于用户Y'轴位置域或基于模糊度域进行 整周模糊度一维空间搜索，根据步骤2得到的双差线性化定位观测方程进行对应整周模糊 度和用户Y'轴位置的求解，并根据整周模糊度和用户Y'轴位置计算用户到卫星的双差 距离，再获得所述双差距离的计算值与测量值的差值，其中最小差值对应的用户Y'轴位置 即为最终结果；
[0011] 步骤4、获得步骤3的最终结果在坐标系(V- Z'中对应的用户三维位 置，将该用户三维位置转换到步骤1所述的原始坐标系中，完成定位过程。
[0012] 较佳的，所述波长A为GNSS载波波长，或者是其超宽巷和宽巷组合波长。
[0013] 较佳的，所述步骤3中基于用户Y'轴位置域进行整周模糊度一维空间搜索的具 体方法为：
[0014] 在坐标系O' -X'Y'Z'下，以整周模糊度初值▽AN。为中心，确定整周模糊 度搜索范围为：▽ANQ_dN到▽AN〇+dN，其中(^ =打〇〇1'((11?/人），；|^1〇〇1'(.)表示取整，(11? 为用户初始定位误差；将整周模糊度搜索范围内的每个整周模糊度代入到步骤2的双差线 性化定位观测方程（2.3)中，分别得到对应的误差改正数Sy'u，根据误差改正数Sy'u 计算用户Y'轴位置。
[0015] 较佳的，所述步骤3中基于模糊度域进行整周模糊度一维空间搜索的具体方法 为：
[0016] 在坐标系O' -X'Y'Z'下，以用户初始点位为中心，以设定的搜索步长在一维 空间y'uQ-dR到y'uQ+dR范围内对用户Y'轴位置进行分别取值，并代入到步骤2的双差 观测方程公式（2.3)中，得到对应的整周模糊度▽AN;其中，y' u。表示用户初始点位在 Y'轴的坐标，dR为用户初始定位误差。
[0017] 较佳的，在所述步骤4中，选取多颗非基准星，先分别采用步骤1到步骤3的方法， 获得各颗非基准星对应的用户Y'轴位置；然后再利用步骤4的方法将各个用户三维位置 转换到步骤1所述的原始坐标系中，选取这些位置的集中区域，确定集中区域中各位置点 的用户三维位置并求均值，则该均值即为定位结果。
[0018] 较佳的，在所述步骤4中，选取多颗非基准星，先分别采用步骤1到步骤3的方法， 获得各颗非基准星对应的用户Y'轴位置；然后再利用步骤4的方法将各个用户三维位置 转换到步骤1所述的原始的三维坐标系中；选取这些位置的集中区域，确定集中区域中各 位置点的用户三维位置及其对应整周模糊度值；
[0019] 利用集中区域内各位置点对应的整周模糊度值，再结合载波相位测量值，根据载 波相位双差观测方程公式，利用传统最小二乘进行定位计算。
[0020] 较佳的，所述集中区域内两点之间距离小于2倍的载波波长。
[0021] 较佳的，所述设定的搜索步长为0. 1入。
[0022] 本发明具有如下有益效果：
[0023] (1)本发明的定位方法克服了现有基于三维位置域整周模糊度搜索技术的不足， 可以从很大程度上减少模糊度搜索计算量，实现整周模糊度位置域搜索方法的快速确定。 同时该方法在整周模糊度确定的同时会排除观测量误差较大的卫星，完成高精度定位，有 效的保障定位精度和可靠性。该方法可用于包括GPS、GLONASS、Compass、Galileo等在内 多星座GNSS兼容系统或单系统，实现整周模糊度的快速确定和高精度定位。
[0024] (2)本发明的定位方法通过选择多颗非基准星分别计算用户位置，并排除集中区 域之外的用户位置，对集中区域之内的用户位置进行求平均作为最后的定位结果，使得定 位更加精准可靠。
【附图说明】
[0025] 图1为本发明的三维转一维位置域搜索及定位流程图；
[0026] 图2为本发明的三维转一维位置降维处理分析结果。
【具体实施方式】
[0027] 下面结合附图并举实施例，对本发明进行详细描述。
[0028] 卫星导航载波相位高精度定位通常采用站星双差方式，其线性化的载波相位双差 观测方程公式一般可以表示为公式（2. 1):
[0029]
[0030] 其中，u对应用户接收机，c对应基准站，k对应基准星，i对应非基准星， Sxu，Syu，5%指用户三维位置坐标的误差改正数，％禮是双差载波相位观测值，是 整周模糊度，A是载波相位波长（这里波长可以是GNSS载波波长，或者是其超宽巷、宽巷 组合载波波长），V'M/是利用卫星坐标和用户坐标（x，y，z)计算的用户到卫星的双差距离， ▽A 是双差电离层时延，▽AST是双差对流层时延，是用户到卫星k和 i的单位方向余弦矢量差在三维坐标系中的表示，该单位方向余弦矢量差值越大，对整周模 糊度影响越大；Noise表不误差。
[0031] 公式（2.1)是基于卫星导航系统所定义的坐标系，这里定义为〇->XYZ，在此 坐标系下，传统的基于位置域的整周模糊度搜索方法都是基于三维位置空间进行整周模 糊度搜索，基于三维位置域整周模糊度搜索的优点是对载波相位的整周不敏感，缺点是计 算量大搜索速度慢，在三维几何空间进行搜索，假设用户概略位置误差为dR，搜索步长为 入/1〇,则每一维搜索空间对应收索次数N=floor(10 ?dR/A)，整个模糊度搜索三维空间 是dRXdRXdR，对应搜索次数为N3,计算量非常大，尽管后续算法有所改进但是因为计算量 大原因近10年没有得到更进一步的突破发展。此外，基于传统三维位置域载波相位双差观 测方程的最小二乘定位计算方法在某颗卫星观测量误差较大的情况下无法自动检测出来， 进而会导致某一时间点的定位误差较大。
[0032] 基于此，本发明的核心思想是：克服现有基于三维位置域整周模糊度搜索技术的 不足，提出一种三维转一维的整周模糊度降维搜索及定位方法，其原理是通过建立一种新 的坐标系〇' - >X'Y'Z'，使新坐标下双差定位观测方程公式（2. 1)对应的X'，Z'方 向的单位方向余弦矢量差值远小于Y'方向，因此通过坐标转换后，三维位置对整周模糊度 的影响主要集中在Y'轴方向。故假设用户概略位置误差是dR，三维转一维位置最大误差 仍是初始位置定位精度误差，根据三维转一维的线性关系确定域模糊度搜索范围为dR，收 索次数N=floor(10 ?(!!?/A)，即搜索次数从三维位置域搜索的N3缩减为三维转一维位置 域的N，因此能够从更大程度上减少了计算量。该方法一方面能从很大程度上减少整周模糊 度搜索计算量，可实现整周模糊度的快速确定，另一方面该方法在整周模糊度确定的同时 会自动检测排除观测量误差较大的卫星，能够更有效的保障定位精度和可靠性。本方明可 用于包括GPS、GLONASS、Compass、Galileo等在内多星座GNSS兼容系统或单系统，实现整 周模糊度的快速确定和高精度定位。
[0033] 如图1所示，本发明的基于三维转一维的载波相位整周模糊度搜索的定位方法主 要包括如下步骤：
[0034] 通过非几何构型方法初步估算整周模糊度初值
[0035] 利用卫星位置和用户初始概略位置初步计算用户及基准站到卫星k和卫星i的双 差距离▽AR。，通过非几何构型方法利用直接取整的方法初步估算可视卫星对应的双差整 周模糊度初值：
[0036] ,.、 (2.2)
[0037] 其中，▽AR。是用卫星位置（Xk，yk，Zk)，（Xi，yi，Zi)和用户初始概略位置 (xuD，yu(]，zj初步计算的用户及基准站到卫星k和卫星i的双差距离，A是载波波长（这 里，该波长可以是GNSS载波波长，或者是其超宽巷组合或宽巷组合载波波长），▽A(}>是双 差载波相位观测值，▽AN。是所计算的双差整周模糊度初值。
[0038] 步骤1、坐标转换到新坐标系建立新的双差观测方程
[0039] 在可视范围内，一般取高度角最大的一颗可视卫星为基准星k，在用户u、基准星 k与非基准星i组成的平面中，坐标系选取Z'轴垂直于该平面，垂直于两颗卫星连线方向 为r轴，建立右手螺旋坐标系(V-Y'Z、将卫星坐标 用户初始概略坐标（xu。，yu。，zu。）从原坐标系0 -XYZ转到新坐标系0' -X'Y'Z'下(x' k，y' k，z' k)，（x' i，y' i，z' i)(x' u0，y' u0，z'u0)〇
[0040] 经过坐标转换，通过仿