基于全伪距和部分伪距的定位方法技术

本申请公开了一种基于全伪距和部分伪距的定位方法。该定位方法包括：获得n颗卫星的全伪距观测量和星历，n≥4；获得m颗卫星的部分伪距观测量和星历，m≥1；构建联合状态向量，所述联合状态向量包括接收机位置、接收机时钟偏差、以及所述获得部分伪距观测量的m颗卫星所对应的伪距整数模糊度；以及根据全伪距观测量与接收机位置和接收机时钟偏差之间的关系，以及根据部分伪距观测量与接收机位置、接收机时钟偏差和与所述伪距整数模糊度之间的关系，以及所有卫星的星历，对联合状态向量进行迭代更新，获得联合状态向量的估计值。通过该定位方法，能够利用全伪距观测量和部分伪距观测量联合完成定位解算，从而提高了定位精度，缩短了定位时间。

【技术实现步骤摘要】

本申请涉及卫星导航
，具体地，涉及基于全伪距和部分伪距的定位方法。
技术介绍
卫星导航接收机启动后需要依次完成卫星信号捕获、码同步、比特同步、帧同步、导航电文解调等信号处理操作，在获取足够多颗卫星的全伪距观测量并具备有效星历时，可以解算得到接收机位置。星历已知时，首次定位时间主要受限于足够数量的全伪距观测量的获取，一般情况下接收机至少需要完成帧同步才能获得全伪距，完成码同步、比特同步只能获得部分伪距观测量。传统卫星导航接收机只能在获取不少于4颗卫星的全伪距观测量后才有可能完成定位解算。但是，如果4颗获得了全伪距观测量的卫星几何分布较差，则单纯依靠该4个全伪距观测量得到的定位结果误差较大或根本无法获取定位结果。如需改善定位结果，可使用其他卫星全伪距观测量参与定位解算，但必需等待对其他卫星信号的帧同步完成后才能获得其全伪距观测量参与定位解算，该过程中这些其他卫星的部分伪距观测量无法参与定位解算，使得改善定位精度的等待时间较长。尤其，北斗卫星系统中的GEO卫星和非GEO卫星具有不同的导航电文广播速率，非GEO卫星的导航电文广播速率为50bps，完成帧同步并获得全伪距观测量至少需要6s时间；而GEO卫星导航电文广播速率为500bps，是非GEO速率的10倍，因此完成帧同步获得全伪距观测量的时间理论上仅需要0.6s。北斗系统的GEO卫星的空间位置几何分布较差，仅通过GEO卫星的全伪距观测量无法实现精度较高的定位，而等待其他非GEO卫星全伪距观测量的获取则需要至少6s时间，定位时间较长。
技术实现思路
本申请的目的是提供一种基于全伪距和部分伪距的定位方法，能够利用全伪距观测量和部分伪距观测量联合完成定位解算，从而提高定位精度，缩短定位时间。根据本申请，公开了一种基于全伪距和部分伪距的定位方法，包括：获得n颗卫星的全伪距观测量和星历，n≥4；获得m颗卫星的部分伪距观测量和星历，m≥1；构建联合状态向量，所述联合状态向量包括接收机位置、接收机时钟偏差、以及所述获得部分伪距观测量的m颗卫星所对应的伪距整数模糊度；以及根据全伪距观测量与接收机位置和接收机时钟偏差之间的关系，以及根据部分伪距观测量与接收机位置、接收机时钟偏差和与所述伪距整数模糊度之间的关系，以及所有卫星的星历，对联合状态向量进行迭代更新，获得联合状态向量的估计值。附图说明图1示出了根据本申请的一种实施方式的定位方法。图2示出了全伪距观测量与部分伪距观测量和伪距整数模糊度之间的关系。具体实施方式下面参照附图对本申请公开的基于全伪距和部分伪距的定位方法进行详细说明。为简明起见，本申请各实施例的说明中，相同或类似的装置使用了相同或相似的附图标记。图1示出了根据本申请的一种实施方式的定位方法，根据该定位方法，可以基于全伪距观测量和部分伪距观测量进行联合定位方法。如图所示，在S101中，获得n颗卫星的全伪距观测量和星历，其中n≥4。在S102中，获得m颗卫星的部分伪距观测量和星历，其中m≥1。图2示出了全伪距观测量与部分伪距观测量和伪距整数模糊度之间的关系。接收机完成帧同步后获得全伪距。接收机完成码同步或者比特同步，则获得缺少伪距整数模糊度N个整周距离cT的部分伪距观测量。整周距离cT是码周期或者比特周期时间T内光传播距离。在S103中，构建联合状态向量。联合状态向量包括接收机位置、接收机时钟偏差、与获得部分伪距观测量的m颗卫星所对应的伪距整数模糊度。在S104中，根据全伪距观测量与接收机位置和接收机时钟偏差之间的关系，以及根据部分伪距观测量与接收机位置、接收机时钟偏差和与获得部分伪距观测量的卫星所对应的整数模糊度之间的关系，以及所有卫星星历，对联合状态向量进行迭代更新，获得联合状态向量的估计值。这样，也就获得了收机位置估计值和接收机时钟偏差估计值。例如，当获得全伪距观测量的n颗卫星是北斗GEO卫星，获得部分伪距观测量的m颗卫星是北斗非GEO卫星时，在具备有效星历并获取不少于4颗GEO卫星的全伪距观测量和至少1颗非GEO卫星的部分伪距观测量时即可得到较高精度的定位结果，并且理论上可以将原有需要6s以上时间完成的定位解算缩短至0.6s，具有10倍的提升效果。下面，结合GEO卫星和非GEO卫星的应用场景对根据本申请的定位方法进行详细说明。根据本申请的一种实施方式，对于接收到的GEO卫星信号，完成帧同步并获得n颗GEO卫星的全伪距观测量ρi，其中，i＝1,…,n，n≥4。对于接收到的非GEO卫星信号，仅完成码同步或者比特同步，获得m颗非GEO卫星的部分伪距观测量ρj，j＝1,…,m，m≥1。全伪距观测量ρi与接收机位置和接收机时钟偏差之间的关系可以表示为：ρi＝|Pi-xr|+b+εi，部分伪距观测量ρj与接收机位置、接收机时钟偏差和伪距整数模糊度之间的关系可以表示为：ρj＝|Pj-xr|+b-NjcT+εj，其中，Pi和Pj分别表示GEO卫星和非GEO卫星的卫星位置，xr表示接收机位置，b表示接收机时钟偏差(单位：米)，εi和εj分别表示与GEO卫星和非GEO卫星对应的其他误差，Nj表示非GEO卫星的伪距整数模糊度，cT表示整周距离。由于获得了全伪距观测量的GEO卫星几何分布较差，如果根据全伪距观测量解算获得定位结果，则定位结果的误差较大或根本无法获取定位结果。对于获得了部分伪距观测量的非GEO卫星信号，由于整数模糊度未知，无法直接进行伪距观测量解算获得定位。如果等待帧同步完成获得其全伪距观测量后参与定位解算，将使得改善定位精度的等待时间较长。根据本申请的一种实施方式，构建联合状态向量X：X=xrbN=xyzbN1...Nm,]]>联合状态向量X包括接收机位置xr(三维)、接收机时钟偏差b、非GEO卫星所对应的伪距整数模糊度N。构建联合状态向量X后，可以通过迭代更新的方式，获得联合状态向量的估计值。这样，也就可以获得接收机位置估计值和接收机时钟偏差估计值，从而实现了接收机位置的精确估计。联合状态向量X的初始值可以任意设置，例如，可以将联合状态向量X设置为零向量。可以根据全伪距观测量与接收机位置和接收机时钟偏差之间的关系，以及根据部分伪距观测量与接收机位置、接收机时钟偏差和伪距整数模糊度之间的关系，以及卫星星历，对联合状态向量进行迭代更新，获得联合状态向量的估计值。为了获得每次迭代更新中联合状态向量的迭代更新量δX，δX＝[δxδyδzδbδN1…δNm]T(上标T表示转置)，需要在每一次迭代运算中计算发射时刻的卫星位置，因此，需要对卫星发射时刻进行估计。根据一种实施方式，GEO卫星和非GEO卫星的发射时刻可以通过如下方式计算：tt,i＝tr,i-ρi/c-te,i，tt,j＝tr,j-(ρj+NjcT)/c-te,j，其中，tt,i表示第i颗GEO卫星的发射时刻，tt,j表示第j颗非GEO卫星的发射时刻，tr,i和tr,j分别表示与第i颗GEO卫星和第j颗非GEO卫星相对应的接收机接收时刻，te,i和te,j分别表示与第i颗GEO卫星和第j颗非GEO卫星相对应的卫星时钟误差修正量，可以通过广播卫星星历获得，c表示光速。由于非GEO卫星伪距观测量采用了部分伪距观测量ρj，因此，在计算非GEO卫星的发射时刻时，引入了伪距整数模糊度估计值Nj。随着联合状本文档来自技高网...

【技术保护点】
一种基于全伪距和部分伪距的定位方法，包括：获得n颗卫星的全伪距观测量和星历，n≥4；获得m颗卫星的部分伪距观测量和星历，m≥1；构建联合状态向量，所述联合状态向量包括接收机位置、接收机时钟偏差、以及所述获得部分伪距观测量的m颗卫星所对应的伪距整数模糊度；以及根据全伪距观测量与接收机位置和接收机时钟偏差之间的关系，以及根据部分伪距观测量与接收机位置、接收机时钟偏差和与所述伪距整数模糊度之间的关系，以及所有卫星的星历，对联合状态向量进行迭代更新，获得联合状态向量的估计值。

【技术特征摘要】
1.一种基于全伪距和部分伪距的定位方法，包括：获得n颗卫星的全伪距观测量和星历，n≥4；获得m颗卫星的部分伪距观测量和星历，m≥1；构建联合状态向量，所述联合状态向量包括接收机位置、接收机时钟偏差、以及所述获得部分伪距观测量的m颗卫星所对应的伪距整数模糊度；以及根据全伪距观测量与接收机位置和接收机时钟偏差之间的关系，以及根据部分伪距观测量与接收机位置、接收机时钟偏差和与所述伪距整数模糊度之间的关系，以及所有卫星的星历，对联合状态向量进行迭代更新，获得联合状态向量的估计值。2.如权利要求1所述的定位方法，其中，所述每一次迭代更新包括：根据所述全伪距观测量计算所述n颗卫星的发射时刻估计值，根据卫星星历和所述n颗卫星的发射时刻估计值获得所述n颗卫星的卫星位置估计值；根据所述部分伪距观测量和伪距整数模糊度估计值迭代估计所述m颗卫星的发射时刻估计值，根据卫星星历和所述m颗卫星的发射时刻估计值获得所述m颗卫星的卫星位置估计值；以及根据全伪距观测量与接收机位置和接收机时钟偏差之间的关系，部分伪距观测量与接收机位置、接收机时钟偏差和伪距整数模糊度之间的关系，以及所有卫星的卫星位置估计值，通过最小二乘法计算伪距向量的在当前迭代更新中的迭代更新量。3.如权利要求2所述的定位方法，其中，所构建的联合状态向量X为：X＝[xyzbN1…Nm]T；所述联合状态向量的迭代更新量δX通过以下方式计算：δX＝(HTH)-1HTδρ其中，矩阵H表示联合设计矩阵，其中...

【专利技术属性】
技术研发人员：赵思浩，崔晓伟，陆明泉，
申请(专利权)人：清华大学，
类型：发明
国别省市：北京;11