2024.1.30 GNSS 学习笔记

站星双差Kalman滤波伪距差分定位流程

1.

---

RTK定位技术（实时载波相位差分技术）原理-站间单差浮点解

1.RTK技术其实就是在RTD技术的基础上增加载波观测值的使用。由于伪距的噪声在分米量级，即使我们通过站间单差消除了绝大部分的误差影响，但受限于伪距的精度，我们也只能达到分米量级的定位水平。
但载波不同，载波的精度在毫米量级，所以如果其他误差完全消除的话，理论上定位精度可达到的毫米水平。毫米水平的定位精度，完全可以满足我们日常生活中绝大多数应用的需求。但实际上，站间单差后依然存在误差，定位精度在厘米量级。

2.载波观测方程是怎么建立的？

由于载波频率高，波⻓短，所以载波相位测量精度高。

但是方程中还有一项误差没有表示出来，即相位硬件延迟。前面我们了解到不管是卫星端还是接收机端，均存在因信号通道引入的时间延迟，我们称之为伪距硬件延迟。同样的对于载波，也存在相位硬件延迟，产生机理与伪距相同。
我们是否可以通过类似伪距码延迟处理的手段，来处理相位硬件延迟？
首先，对于卫星端相位硬件延迟，已经具备成熟的技术对其进行估计，终端通过使用相应的产品，可以实现模糊度固定，称之为 PPPAR技术。
在RTK技术中，卫星端的相位硬件延迟，可通过站间单差消除。但接收机端的相位硬件延迟，因为只关注不满足一周的部分，所以量级太小，无法估计。但可以通过星间单差消除其影响，但是何时进行星间单差，有一定讲究。

为什么我们一定要消除相位硬件延迟的影响？
因为只有消除影响，模糊度才具有整数特性，我们才可以通过模糊度固定，进一步提升定位结果的精度和可靠性。否则如果不固定模糊度，即模糊度是浮点解，定位精度一般在cm-dm量级，相对于伪距差分有一定的精度提升；但如果模糊度固定成功，那么定位精度就是cm量级。

---

站间单差浮点解方程

1.

2.再增加载波观测值，首先是状态量，每个卫星每个频点均要增加一个模糊度参数，因为有12个观测值，所以要增加12个模糊度参数，也要增加一倍的观测值数目。

3.

问题：

理论上，对于单个系统指定卫星指定频率的伪距和载波的接收机钟肯定不相等的，但我们只估计了伪距的接收机钟参数，同样的载波也是估计该接收机钟参数，这样设计是否合理？

回答：
对于同一系统同一频点的伪距和载波，肯定存在接收机钟的差异，但该偏差会被我们新增的每个观测值的模糊度吸收，所以这也是模糊度失去整数特性的原因之一。

4.

同样的，设计矩阵也更加庞大

---

RTK定位技术原理-站星双差浮点解

1.站星双差浮点解方程

2.kalman滤波以及双差模糊度参数转换

3.双差观测值但使用站间单差的模糊度参数

---

数据预处理包含哪几个方面？

1.观测值剔除

卫星信号微弱，极易受环境影响。比如，如果有树荫遮挡或者高楼的信号反射，卫星观测值的噪声会明显增大，也会引入粗差。考虑到卫星高度⻆较低的卫星容易受环境影响，所以我们一般设置一定的卫星高度截止⻆，用来剔除可能存在误差的卫星观测值。
我们一般设置卫星高度截止⻆为10°，但在一些复杂场景下可能会提高卫星高度截止⻆，以达到剔除粗差的目的。
同时，对于信噪比较低的卫星，其观测值也容易引入粗差，所以我们也会设置一定的信噪比阈值，对观测值进行剔除。针对测量性接收机，一般可设置38作为信噪比阈值。对于消费级模组或者芯片，一般其观测值信噪比均相对较低，可以适当将信噪比阈值降低。不同的接收机观测值特性和信噪比存在差异，所以算法需要不同的接收机设置合理的阈值。

2.周跳探测

产生原因：

周跳的实质就是模糊度发生变化，并不仅是产生了整周的变化，也会产生小数部分的变化。客观原因可能是卫星信号被障碍物遮挡，或者外界干扰导致卫星信号失锁。

个人经验，现在的周跳基本均无法修复，因为并不是整周跳变。所以对于出现周跳的模糊度，直接重置。即将模糊度设置一个新的数值，给一个较大的方差。

周跳探测的方法：

注：上面GF组合观测值的公式中，有地方写错了，应该是频率1的波长乘以模糊度值。

---

RTK定位技术原理-模糊度固定技术

1.双差浮点模糊度

其实上述内容中，已经讲了两种形式的模糊度估计方式--站间单差和站星双差。
如果我们的模糊度参数是站星双差形式，表明已消除卫星端和接收端的相位硬件延迟的影响，我们可以直接将双差模糊度进行搜索，尝试得到双差模糊度的固定解。
如果我们的模糊度参数是站间单差的形式，那么我们还要在模糊度参数层面组星间单差，进而消除接收机端相位硬件延迟的影响。再组星间单差，和之前的模糊度预测类似，填充一个对应的状态转移矩阵，但并不需要对方差协方差矩阵增加Q阵。

2.模糊度搜索

现在使用较多且最有效得模糊度固定方法，即LAMBDA方法（Leastsquare AMBiguity Decorrelation Adjustment，LAMBDA，最小二乘模糊度降相关平差法）。
该方法主要由两部分组成：
（1）为降低模糊度参数之间相关性而进行的多维整数变换；
（2）在转换后的空间内进行模糊度搜索，然后再将结果转换回模糊度空间中，进而求得模糊度整数解

2.1 模糊度参数降相关

由于模糊度之间具有相关性，一个模糊度参数的变化会影响其他模糊度的搜索，使得搜索算法的计算量巨大。如果能够设法降低模糊度参数之间的相关性，使得某一模糊度的变化对其他模糊度的取值的影响尽可能小，就能大大加快模糊度的搜索过程。

LAMBDA算法通过对模糊度参数及其方差-协方差阵进行整数高斯变换（也称z变换），将他们从原空间变换到新的空间中去，实现模糊度的降相关：

2.2 模糊度搜索

模糊度搜索与前述的降相关是LAMBDA算法的两个独立的部分，就算不做降相关操作，依然可以进行模糊度搜索，整数模糊度搜索条件：

参考文献：GNSS整周模糊度固定算法（LAMBDA算法）_lambda算法原理gps-CSDN博客

3.模糊度固定解

经过模糊度搜索后，我们可以得到双差模糊度的固定解，恢复了模糊度的整数特性。我们可以直接将固定后的模糊度值作为虚拟观测值，对参数进行约束，直接得到模糊度固定解。
所谓的作为虚拟观测值对参数进行约束，即设计矩阵相应的系数为1，残差为观测值减去预测值。同时给予该观测值比较大的权重，即给比较小的噪声（比如设置为1毫米）。

3.1 使用滤波更新固定解

3.2 直接更新固定解

---

RTK定位流程

1.RTK定位引擎

2.在线定位

3.离线定位（后处理定位）

4.单个历元的定位流程