一种GNSS三频精密单点定位模糊度固定方法技术

本发明专利技术涉及一种GNSS三频精密单点定位模糊度固定方法，属于卫星导航定位技术领域。该方法基于IGS提供的精密卫星星历与钟差以及网络解算的相位延迟改正，利用三频伪距与载波相位观测值实现精密单点定位的解算方法；利用超宽巷相位伪距组合解算超宽巷整周模糊度，并结合解算的无电离层组合观测值模糊度，解算窄巷整周模糊，最终获得模糊度固定的高精度单点定位结果。本发明专利技术提供的一种GNSS三频精密单点定位模糊度固定方法，减少了待估参数的个数，提高了计算效率，同时提高了精密单点定位模糊度解算成功率，以及定位精度与可靠性。

【技术实现步骤摘要】
一种GNSS三频精密单点定位模糊度固定方法
本专利技术属于卫星导航定位
，涉及一种GNSS三频精密单点定位模糊度固定方法。
技术介绍
全球卫星导航系统(GlobalNavigationSatelliteSystem，GNSS)能够为地球表面和近地空间的广大用户提供全天时、全天候、高精度的定位、导航和授时服务，已广泛应用于国防、海陆空交通运输、测绘、移动通信、电力、电子金融、精细农业和减灾救灾等领域，是拓展人类活动和促进社会发展的重要空间基础设施。GNSS精密单点定位(PPP)技术是全球卫星导航系统提供精密位置服务的关键技术之一，它利用卫星精密星历及卫星钟差，以单台双频接收机采集的非差相位和伪距数据作为主要观测值来进行单点定位，可以直接得到分米级乃至厘米级的高精度的三维位置，极大的提高了作业效率，降低作业成本。因而在其出现之初，即受到科研与应用领域的广泛的关注，获得了较快发展。在过去10年间，PPP技术从双频发展到单频、从单系统发展到多系统组合、从事后发展到实时、从浮点解发展到固定解。其中精密单点定位的模糊度固定解可以获得厘米级的定位精度，是目前的卫星导航领域中的研究热点。但是基于目前在轨GPS卫星大多发送双频伪距与载波信号，使得精密单点定位的模糊度固定需要较长的初始化时间，从而限制了其应用范围。随着美国GPS的现代化、俄罗斯的GLONASS更新、欧盟的GALILEO发展、中国的BDS后来居上，未来的卫星导航系统都将在三个或三个以上频点播发卫星导航信号，因而多系统多频精密定位将成为未来GNSS高精度测量应用的发展方向。相对于双频PPP而言，三频PPP的数据处于理(周跳探测)将更加可靠，且可有效减少模糊度解算的初始化时间，不仅对实时高精度导航定位服务具有重要意义，而且对于北斗在高精度测量领域的产业化推广具有重要意义。
技术实现思路
有鉴于此，本专利技术的目的在于提供一种GNSS三频精密单点定位模糊度固定方法，该方法利用三频载波相位实现高精度单点定位。为达到上述目的，本专利技术提供如下技术方案：一种GNSS三频精密单点定位模糊度固定方法，该方法包括以下步骤：步骤一：基于第i历元原始的载波相位与伪距观测值，形成两个独立的无电离层组合观测值；步骤二：基于对流层模型、地球固体潮模型、卫星天线相位中心模型修正误差源的影响；步骤三：计算观测方程系数矩阵，并利用Kalman滤波技术估计接收机三维位置与无电离层组合观测值的载波相位模糊度；步骤四：基于原始的载波相位与伪距观测值，形成超宽巷伪距相位组合，并采用移动平滑技术，获得超宽巷整周模糊度及其模糊度解算的成功率；步骤五：基于步骤三解算的无电离层组合观测值模糊度及步骤四解算的超宽巷模糊度，通过以下公式，形成第一与第二频段的模糊度值及其方差协方差矩阵；其中，f1，f2，f3为观测值的频率，NIF1，NIF2为三频载波相位观测值形成的两个独立的无电离层组合观测值的模糊度，为超宽巷组合观测值的模糊度；步骤六：解算载波第一与第二频段的窄巷整周模糊度并检验，更新定位信息。进一步，所述步骤一通过以下公式形成两个独立的无电离层组合观测值，其中，为伪距组合观测值，为相位组合观测值，分别为伪距与相位组合观测值的噪声；为无电离层组合观测值的模糊度；为无电离层组合观测值的波长；ρ为接收机到卫星的距离；m为对流层延迟的投影函数；τ为天顶对流层延迟，tIFj为无电离层组合观测值的接收机钟差。进一步，所述步骤三Kalman滤波技术的Kalman滤波的过程为，其中，yk为无电离层伪距与载波相位组合观测值，Ak为观测值的系数矩阵，为k-1历元的状态向量，为预测的状态向量，为当前k历元的状态向量，Tk|k-1为状态转移矩阵。进一步，所述步骤四中计算超宽巷模糊度及其模糊度解算的成功率具体包括以下步骤：1)通过以下公式形成超宽巷伪距相位组合，其中，L(0，1，-1)为超宽巷伪距载波相位组合观测值，下标(0，1，-1)表示组合观测值的系数，λ(0，1，-1)为超宽巷组合观测值的模糊度，组合观测值的模糊度N(0，1，-1)=(N2-N3)；ε(0，1，-1)为组合观测值噪声；2)采用移动平滑技术，得到平滑的组合观测值，其中，k为平滑历元个数；为k历元平滑的超宽巷伪距相位组合观测值；为的方差；3)通过以下公式计算超宽巷模糊度，式中代表取整操作，S(0，1，-1)表示网络解算的超宽巷UPD值；4)通过以下公式计算宽巷模糊度解算成功率，其中，进一步，所述步骤五更新位置信息具体包括以下步骤：1)基于MLAMBDA，通过以下公式解算整周模糊度，其中，为GNSS网络解算的两个频率段UPD值；2)通过以下公式计算模糊度的成功率与Ratio值，其中，PBT为对应的Bootstrapping成功率；为模糊度方差矩阵经LDLT分解的条件方差；其中，r为Ratio值；为与实数模糊度最近的整周模糊度矢量，为与实数模糊度次最近整周模糊度矢量；3)如果通过整周模糊度正确检验，则通过以下公式更新位置信息，其中，b表示模糊度浮点解的接收机三维位置矢量，表示三维位置矢量与模糊度之间的方差协方差矩阵，Pbb为Kalman滤波获得的三维位置的方差协方差矩阵，为模糊度固定后三维位置的方差协方差矩阵，sIF为GNSS网络解算的无电离层组合观测值的UPD；本专利技术的有益效果在于：本专利技术所提供的一种GNSS三频精密单点定位模糊度固定方法，该方法基于IGS提供的精密卫星星历与钟差以及网络解算的相位延迟改正，利用三频伪距与载波相位观测值实现精密单点定位的解算方法；利用超宽巷相位伪距组合解算超宽巷整周模糊度，并结合两个独立无电离层组合观测值，解算宽巷与窄巷整周模糊，最终获得模糊度固定的高精度单点定位结果，本专利技术所提供的一种GNSS三频精密单点定位模糊度固定方法，减少了待估参数的个数，提高了计算效率，同时提高精密单点定位模糊度解算成功率，提高定位精度与可靠性。附图说明为了使本专利技术的目的、技术方案和优点更加清楚，下面将结合附图对本专利技术作进一步的详细描述，其中：图1为本专利技术所述方法的流程图。具体实施方式下面将结合附图，对本专利技术的优选实施例进行详细的描述。步骤1：利用原始的伪距与载波相位观测值形成两个无电离层组合观测值GNSS三个频率原始的伪距与载波相位观测值为：式中,i表示载波相位与伪距的频率数,i＝1,2,3,Pi为原始的伪距观测值,Li表示以米为单位的原始的载波相位观测值,ρ为接收机到卫星的距离(包含误差源，如相对论、地球固体潮、卫星相位中心、硬件延迟等),tr，i为不同频率接收机钟差，λi为载波波长,I1第一频段载波的电离层延迟,为电离层系数，τ为天顶对流层延迟,m为对流层延迟的投影函数，Ni为Li载波的模糊度值,∈i表示Pi的观测噪声，εi表示Li的观测噪声。表1GNSS载波的频率(Mhz)GPS/GALILEO/BDS的三个频率的信号如表1所示，利用三个频率的观测值可以形成两个独立的无电离层组合观测值，如下所示：式中,为伪距组合观测值，tIFj为无电离层组合观测值的接收机钟差。为相位组合观测值，分别为伪距与相位组合观测值的噪声与,为无电离层组合观测值的模糊度,为无电离层组合观测值的波长。步骤2和步骤3：基于对流层模型、地球固体潮模型、卫星天线相位中心模型修正本文档来自技高网...

【技术保护点】
一种GNSS三频精密单点定位模糊度固定方法，其特征在于：该方法包括以下步骤：步骤一：基于第i历元原始的载波相位与伪距观测值，形成两个独立的无电离层组合观测值；步骤二：基于对流层模型、地球固体潮模型、卫星天线相位中心模型修正误差源的影响；步骤三：计算观测方程系数矩阵，并利用Kalman滤波技术估计接收机三维位置与无电离层组合观测值的载波相位模糊度；步骤四：基于原始的载波相位与伪距观测值，形成超宽巷伪距相位组合，并采用移动平滑技术，获得超宽巷整周模糊度及其模糊度解算的成功率；步骤五：基于步骤三解算的无电离层组合观测值模糊度及步骤四解算的超宽巷模糊度，通过以下公式，形成第一与第二频段的模糊度值及其方差协方差矩阵；N1N2=f1f12-f22-f2f12-f220-f2f32-f22+f3f32-f22-1(NIF1NIF2-0-f3f32-f22N‾(0,1,-1)),]]>其中，f1，f2，f3为观测值的频率，NIF1，NIF2为三频载波相位观测值形成的两个独立的无电离层组合观测值的模糊度，为超宽巷组合观测值的模糊度；步骤六：解算载波第一与第二频段的窄巷整周模糊度并检验，更新定位信息。...

【技术特征摘要】
1.一种GNSS三频精密单点定位模糊度固定方法，其特征在于：该方法包括以下步骤：步骤一：基于第i历元原始的载波相位与伪距观测值，形成两个独立的无电离层组合观测值；步骤二：基于对流层模型、地球固体潮模型、卫星天线相位中心模型修正误差源的影响；步骤三：计算观测方程系数矩阵，并利用Kalman滤波技术估计接收机三维位置与无电离层组合观测值的载波相位模糊度；步骤四：基于原始的载波相位与伪距观测值，形成超宽巷伪距相位组合，并采用移动平滑技术，获得超宽巷整周模糊度及其模糊度解算的成功率；步骤五：基于步骤三解算的无电离层组合观测值模糊度及步骤四解算的超宽巷模糊度，通过以下公式，形成第一与第二频段的模糊度值及其方差协方差矩阵；其中，f1，f2，f3为观测值的频率，NIF1，NIF2为三频载波相位观测值形成的两个独立的无电离层组合观测值的模糊度，为超宽巷组合观测值的模糊度；步骤六：解算载波第一与第二频段的窄巷整周模糊度并检验，更新定位信息；所述步骤一通过以下公式形成两个独立的无电离层组合观测值，其中，为伪距组合观测值，为相位组合观测值，分别为伪距与相位组合观测值的噪声；为无电离层组合观测值的模糊度；为无电离层组合观测值的波长；ρ为接收机到卫星的距离；m为对流层延迟的投影函数；τ为天顶对流层延迟；tIFj为无电离层组合观测值的接收机钟差；所述步骤三Kalman滤波技术的Kalman滤波的过程为，其中，yk为无电离层伪距与载波相位组合观测值，Ak为观测值的系数矩阵，为k-1历元的状态向量，为预测的状态向量，为当前k历元的状态向量；Tk|k-1为状态转移矩阵；所述步骤四中计算超宽巷模糊度及其模糊度解算的成功率具体包括以下步骤：1)通过以下公式形成超宽巷伪距相位组合，其中，L(0,1,-1)为超宽巷伪距载波相位组合观测值，下标(0,1,-1)表示组合观测值的系数，λ(0,1,-1)为超宽巷组合观测值的模糊度，组合观测值的模糊度N(0,1,-1)＝(N2-N3)；ε(0,1,-1)为组合观测值噪声；2)采用移动平滑...

【专利技术属性】
技术研发人员：于兴旺，刘刚，卢艳娥，何盼，
申请(专利权)人：中国科学院重庆绿色智能技术研究院，
类型：发明
国别省市：重庆;85