OSB方法进行模糊度固定

使用武大产品进行模糊度固定的准备工作

数据下载

需要ftp
WHU:ftp://igs.gnsswhu.cn/pub/gnss/products/mgex/2304/WUM0MGXFIN_20240650000_01D_05M_ORB.SP3
WHU:ftp://igs.gnsswhu.cn/pub/gnss/products/mgex/2304/WUM0MGXFIN_20240650000_01D_30S_CLK.CLK
WHU:ftp://igs.gnsswhu.cn/pub/gnss/products/mgex/2304/WUM0MGXFIN_20240650000_01D_30S_OSB.BIA
使用30S的bias文件应该比01D的精度高一些。

使用方法

1、首先通过MW组合产生宽巷模糊度（注意BDS系统建议在进行MW组合时候改正PCO），将OSB文件中对应频点的CODE_BIAS和PHASE_BIAS直接加在原始观测值上；

注意：使用SP3,ERP,OSB,CLK，此时不需要使用CLK头文件中WL-UPD !!!

2、再通过单差无电离层模糊度与经过WL-UPD改正后的单差宽巷模糊度度（此时经过固定），导出单差窄巷模糊度浮点解

3、后续通过LAMBDA方法将单差窄巷浮点模糊度固定，（注意此时单差窄巷模糊度不需要额外的NL-UPD）
不需要再加入额外的WL-UPD,NL-UPD偏差了。

程序修改

之前的程序在computeFixSDNW函数中计算时使用的是：

    NWLd = (
                (F1*(lamda1*oneSatlitData.L1 + oneSatlitData.C1L1C2L2[1]) - F2*(lamda2*oneSatlitData.L2 + oneSatlitData.C1L1C2L2[3]))/(F1-F2) -
                (F1*(oneSatlitData.C1 + oneSatlitData.C1L1C2L2[0] + oneSatlitData.P1_C1) + F2*(oneSatlitData.C2 + oneSatlitData.C1L1C2L2[2]))/(F1+F2)
            ) / lamdaMW;
    currEpoch[i].AmbWL = NWLd;

其中oneSatlitData.C1L1C2L2为OSB文件中的码和相位偏差。

修改为：

NWLd = (
                (F1*oneSatlitData.LL1 - F2*oneSatlitData.LL2)/(F1-F2) -
                (F1*(oneSatlitData.CC1 + oneSatlitData.P1_C1) + F2*oneSatlitData.CC2)/(F1+F2)
            ) / lamdaMW;
    currEpoch[i].AmbWL = NWLd;
    
tempSatlitData.LL1 = tempSatlitData.C1L1C2L2[1] + Lamta1*(tempSatlitData.L1 + tempSatlitData.L1Offset + tempSatlitData.SatL1Offset - tempSatlitData.AntWindup);
tempSatlitData.LL2 = tempSatlitData.C1L1C2L2[3] + Lamta2*(tempSatlitData.L2 + tempSatlitData.L2Offset + tempSatlitData.SatL2Offset - tempSatlitData.AntWindup);
tempSatlitData.CC1 = tempSatlitData.C1 + tempSatlitData.P1_C1 + tempSatlitData.C1L1C2L2[0] + Lamta1*tempSatlitData.L1Offset + Lamta1*tempSatlitData.SatL1Offset;
tempSatlitData.CC2 = tempSatlitData.C2 + tempSatlitData.C1L1C2L2[2] + Lamta2 *tempSatlitData.L2Offset + Lamta2*tempSatlitData.SatL2Offset;

LL1 、LL2 、CC1 、CC2为进行PCO改正后的，加入OSB的原始观测值

修改效果

宽巷由不能收敛固定（上）变为可以固定（下）