103-RTKLIB中对流层延迟模型

1、Troposphere Model

该模型实际上就是Saastamoinen模型，实现函数为tropmodel，求得的是信号传播路径上的电离层延迟。
值得注意的是，该函数中干延迟与湿延迟是分开求取的，然后相加，这样做的目的是可以方便的直接求取天顶干延迟，只需要将卫星高度角设置为90度，相对湿度设置为0.在第三个模型精细的对流层模型中就使用了该函数求天顶干延迟。
模型实现较为简单，只是最后面干延迟求取的公式有点疑问。对照《GPS测量与数据处理》第三版P118的公式，应该乘以(1.0-0.00266*cos(2.0*pos[0])-0.00028*hgt/1E3)部分而不是除以。

2、SBAS troposphere model

在程序中的实现函数为sbstropcorr。在一些材料中，该模型也被称为EGNOS模型，比如“三种对流层延迟改正模型精度评估，曲伟菁等，天文学报，2008”这篇文献。对照该文献，对流层的求取过程也较为容易理解。该模型最后使用了一个简化的投影函数转换为斜对流层延迟，m=1.001/sqrt(0.002001+sinel*sinel);，其实也较为容易理解,如图：

上图来自“丁晓光. 对流层延迟改正在GPS数据处理中的应用与研究[D].长安大学,2009.”。

容易理解sinel=cosz，所以程序中的投影函数就可以简单理解为图中的F，至于更深入的理解，做不到了。
该函数在ppp状态参数的赋值中，将高度角设置为90度，求得天顶方向的对流层延迟，赋值到状态向量x。

3、Precise troposphere model

该模型的实现函数为trop_model_prec。该函数在ppp中估计天顶对流层和东向北向梯度时会用到。
函数的实现较为简单，
调用Saastamoinen计算天顶对流层干延迟，调用NMF模型计算干湿投影函数，根据方位角高度角计算湿延迟及东、北梯度系数dtdx，根据干延迟投影函数、新计算的湿延迟投影函数及干延迟计算信号传播路径上的对流层延迟dtrp。

将对流层三个估值的系数简要说一下。
公式E.5.5可以这样简化：

$$m=m+m\ast cotEl\ast G_N\ast cosAz+m\ast cotEl\ast G_E\ast sinAz$$

等式左边就是顾及东向北向梯度情况下的对流层湿延迟投影函数（系数），右边第一项为NMF模型求得的对流层湿延迟投影函数（系数），第二项和第三项分别为北向和东向电离层梯度系数，后面会将这三个值添加到量测矩阵。

总结

在ppp中，对流层的三个参数作为状态进行估计，那么关于对流层的计算主要包括：
1、EGNOS模型（模型2）计算天顶对流层总延迟，并将这个值赋值给状态向量，北向东向梯度在状态向量中赋值为常数，在udtrop_ppp函数中进行；
2、在相位伪距残差计算函数ppp_res中，使用精细对流层模型计算上述三个参数的系数dtdx（用于确定H阵）及对流层延迟dtrp（用于残差计算）。