PPP时间传递总结

PPP时间传递总结

对于精密单点定位，解算过程中需要较快地进行定位，参数就必须要具备快速收敛和收敛的稳定性。因为在模糊度的固定的过程中，由于各个历元观测值中的模糊度具有较强的相关性，在对条件模糊度按照序列顺序进行固定时，需要一个一个地固定模糊度的值，在固定了一个之后，便会通过其余电离层、对流层等参数验证固定值是否正确，否则取下一个可选的整数模糊度的值，重复上述的过程，所以一般PPP收敛时间为20-30分钟（而在时间传递的板块：当钟差参数的估计方差的平方根小于1ns时，则视作达到了收敛的状态）

PPP 时间传递是一种远距离时间传递方法，在事后精密星历和钟差产品的辅助下，PPP 时间传递可以达到亚纳秒量级的时间传递精度以及 $10^{-15}$/天的频率稳定度。

时间传递的结果一般采用修正的Allan偏差来对时间的传递结果的频率稳定度进行评估
$$\begin{array}{c}\begin{array}{rl}{\sigma }^2(\tau )=& \frac{1}{2{\tau }^2(N-2m)}\sum _{k=1}^{N-2m}(X_{k+2m}-{2X_{k+m}+X_k)}^2\end{array}\end{array}$$
其中$\tau$为历元的间隔，而$N$是总的历元的个数（一天之内取30s为间隔，则N=2880），$X_k$为第k个历元的接收机的钟差的解，m为表示倍数的参数，一般取做10。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-u9SRAyMS-1664801432494)(20220821.assets/image-20220821165008085.png)]

PPP方法计算本地时钟相对于CLK产品基准的钟差，得到单向时间的传递的结果，然后通过计算参与时间传递的两测站钟差结果的二次差，得到站间的时间传递结果。

单系统时间传递

采用无电离层的组合观测方程：
$$\begin{array}{c}\begin{array}{c}{P}_{\mathrm{IF}}=\rho +cd{t}_r-cd{t}^s+d_{r,\mathrm{IF}}-d_{\mathrm{IF}}^s+M_w\cdot {\mathrm{ZTD}}_w+{\epsilon }_{\mathrm{IF}}\\ {\lambda }_{\mathrm{IF}}{L}_{\mathrm{IF}}=\rho +cd{t}_r-cd{t}^s+{\lambda }_{\mathrm{IF}}\left(D_{r,\mathrm{IF}}-D_{\mathrm{IF}}^s\right)+{\lambda }_{\mathrm{IF}}{N}_{\mathrm{IF}}+M_w\cdot {\mathrm{ZTD}}_w+e_{\mathrm{IF}}\end{array}\end{array}$$
d代表伪距硬件延迟，D为载波相位硬件延迟，尾项为为建模的误差和噪声。利用卫星产品可以将卫星端的钟差部分消除，得到：
$$\begin{array}{c}\begin{array}{c}{\tilde{P}}_{\mathrm{IF}}=\rho +cd{t}_r+d_{r,\mathrm{IF}}+M_w\cdot {\mathrm{ZTD}}_w+{\epsilon }_{\mathrm{IF}}\\ {\lambda }_{\mathrm{IF}}{\tilde{L}}_{\mathrm{IF}}=\rho +cd{t}_r+{\lambda }_{\mathrm{IF}}\left(D_{r,\mathrm{IF}}-D_{\mathrm{IF}}^s\right)+d_{\mathrm{IF}}^s+{\lambda }_{\mathrm{IF}}{N}_{\mathrm{IF}}+M_w\cdot {\mathrm{ZTD}}_w+e_{\mathrm{IF}}\end{array}\end{array}$$
从而可以利用Kalman filter来估计所需要的参数$\mathbf{X}=[c\overline{d{t}_r},ZT{D}_w,{\lambda }_{IF}{\hat{\mathbf{N}}}_{IF}]$，这些参数作为上面参数的合并量。
$$\begin{array}{c}\{\begin{array}{l}c\bar{d}{t}_r=cd{t}_r+d_{r,\mathrm{IF}}\\ {\lambda }_{\mathrm{IF}}{\hat{N}}_{\mathrm{IF}}={\lambda }_{\mathrm{IF}}\left(N_{\mathrm{IF}}+D_{r,\mathrm{IF}}-D_{\mathrm{IF}}^s\right)+d_{\mathrm{IF}}^s-d_{r,\mathrm{IF}}\end{array}\end{array}$$

多模GNSS PPP

对于多系统而言，其无电离层组合的伪距和载波相位的方程组为：
{ P ~ I F G = ρ + c d t r + d r , I F G + M w G ⋅ Z T D w + ε I F G P ~ I F R k = ρ + c d t r + d r , I F R k + M w R k ⋅ Z T D w + ε I F R k P ~ I F C = ρ + c d t r + d r , F F C + M w C ⋅ Z T D w + ε I F C P ~ I F E = ρ + c d t r + d r , F F E + M w E ⋅ Z T D w + ε I F E \begin{equation} \left\{\begin{array}{l} \tilde{P}_{\mathrm{IF}}^{\mathrm{G}}=\rho+c d t_{r}+d_{r, \mathrm{IF}}^{\mathrm{G}}+M_{w}^{\mathrm{G}} \cdot \mathrm{ZTD}_{w}+\varepsilon_{\mathrm{IF}}^{\mathrm{G}} \\ \tilde{P}_{\mathrm{IF}}^{\mathrm{R}_{k}}=\rho+c d t_{r}+d_{r, \mathrm{IF}}^{\mathrm{R}_{k}}+M_{w}^{\mathrm{R}_{k}} \cdot \mathrm{ZTD}_{w}+\varepsilon_{\mathrm{IF}}^{\mathrm{R}_{k}} \\ \tilde{P}_{\mathrm{IF}}^{\mathrm{C}}=\rho+c d t_{r}+d_{r, \mathrm{FF}}^{\mathrm{C}}+M_{w}^{\mathrm{C}} \cdot \mathrm{ZTD}_{w}+\varepsilon_{\mathrm{IF}}^{\mathrm{C}} \\ \tilde{P}_{\mathrm{IF}}^{\mathrm{E}}=\rho+c d t_{r}+d_{r, \mathrm{FF}}^{\mathrm{E}}+M_{w}^{\mathrm{E}} \cdot \mathrm{ZTD}_{w}+\varepsilon_{\mathrm{IF}}^{\mathrm{E}} \end{array}\right. \end{equation} ⎩ ⎨ ⎧​P~IFG​=ρ+cdtr​+dr,IFG​+MwG​⋅ZTDw​+εIFG​P~IFRk​​=ρ+cdtr​+dr,IFRk​​+MwRk​​⋅ZTDw​+εIFRk​​