CN114488245A - 一种系统间共频信号差分的gnss/ins组合定位方法 - Google Patents

一种系统间共频信号差分的gnss/ins组合定位方法 Download PDF

Info

Publication number
CN114488245A
CN114488245A CN202210084870.6A CN202210084870A CN114488245A CN 114488245 A CN114488245 A CN 114488245A CN 202210084870 A CN202210084870 A CN 202210084870A CN 114488245 A CN114488245 A CN 114488245A
Authority
CN
China
Prior art keywords
gnss
ins
parameters
disb
satellite
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210084870.6A
Other languages
English (en)
Inventor
李博峰
陈广鄂
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tongji University
Original Assignee
Tongji University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tongji University filed Critical Tongji University
Priority to CN202210084870.6A priority Critical patent/CN114488245A/zh
Publication of CN114488245A publication Critical patent/CN114488245A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/48—Determining position by combining or switching between position solutions derived from the satellite radio beacon positioning system and position solutions derived from a further system
    • G01S19/49—Determining position by combining or switching between position solutions derived from the satellite radio beacon positioning system and position solutions derived from a further system whereby the further system is an inertial position system, e.g. loosely-coupled
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明涉及一种系统间共频信号差分的GNSS/INS组合定位方法，该发明结合GNSS技术兼容与互操作的发展趋势，利用不同卫星系统间的共频信号，基于不同系统相同频率间的差分系统间偏差(Differential Inter‑System Biases,DISB)在时域上的优良稳定性，提出了一种新的GNSS/INS组合定位方法。首先通过超短基线静态实验，得到DISB参数的先验信息；利用该先验DISB信息，基于系统间共频信号差分的GNSS紧组合RTK实现GNSS/INS组合定位。该方法与传统GNSS/INS组合定位算法相比，能够充分挖掘多系统GNSS的信息，显著增加可用信息的个数，从而有效提升GNSS/INS组合系统在复杂环境中的定位性能。

Description

一种系统间共频信号差分的GNSS/INS组合定位方法
技术领域
本发明涉及GNSS/INS组合定位数据处理领域，尤其涉及一种系统间共频信号差分的GNSS/INS组合定位方法。
背景技术
高精度定位技术带动了一系列产业的迅速发展，是无人驾驶、智慧交通、智慧农业等美好愿景能够顺利实现的前提。基于全球卫星导航系统(Global Navigation SatelliteSystem,GNSS)定位技术和惯性导航技术(Inerital Navigation System,INS)的GNSS/INS组合技术毫无疑问是目前室外场景中使用最为广泛的高精度定位手段。虽然该技术投入实际工程使用已将近40年，但仍存在着一些急需解决的问题。特别是对于民用GNSS/INS组合系统而言，由于其严重依赖于GNSS信号，需要GNSS提供足够的信息用于抑制INS误差的累积，当GNSS信号受到干扰或遮挡时，INS误差将快速累积，这将导致组合系统在此类场景中也无法保持长时间高精度的定位定姿。学者们提出了许多办法来改善此情形下GNSS/INS组合系统的定位性能，比如利用车辆载体信息约束、加入其它类型传感器辅助等等。这些方法本质上都是为了引入新的观测信息，代替GNSS观测值用于抑制INS误差的发散。引入新的信息必须考虑新增设备的成本以及信息融合的问题，这将导致实现GNSS/INS组合系统的难度进一步提升。
一方面，目前已有中国的BDS，美国的GPS，欧盟的GALILEO和俄罗斯的GLONASS四大卫星导航系统可为全球用户提供PNT服务，并且日本的QZSS，印度的NAVIC两个区域导航系统可为亚太用户提供PNT服务。为了增加与其他卫星系统的兼容与互操作性，目前各GNSS系统均有意播发了部分与其他卫星系统频率重叠的信号，考虑到相同频率的卫星信号可以共用射频前端从而减少接收机的硬件成本，以及相同频率信号之间优良的兼容与互操作性，可以预见，支持多系统重叠频率将是未来GNSS接收机的发展趋势。
另一方面，目前商用GNSS/INS组合定位算法通常是基于(RealTime kinematic,RTK)+INS的组合，在该类RTK算法中，每个卫星系统会单独使用一颗参考卫星，该类RTK也被称为GNSS系统间松组合RTK(GNSS between-system Loose combination RTK,LRTK)。与之相对应的，还有一类RTK模式被称为GNSS系统间紧组合RTK(GNSS between-system Tightcombination RTK,TRTK)，在该RTK模式中，所有的卫星系统将共用一颗参考卫星。当然，由于不同卫星系统的观测值对应的波长、硬件延迟以及初始相位偏差可能存在差异，导致TRTK比LRTK实现起来更为麻烦。基于重叠频率的TRTK仅需要考虑不同卫星系统间差分硬件延迟(Differential Inter-system Bias,DISB)的影响，DISB在时域上具有非常好的稳定性，这使得基于重叠频率的TRTK比LRTK能够更加充分的利用GNSS观测信息。
发明内容
本发明的目的在于，针对目前GNSS/INS组合定位算法的痛点问题，提供一种系统间共频信号差分的GNSS/INS组合定位方法，用于提高GNSS/INS组合系统在复杂环境中的定位性能。
为了达到上述目的，本发明技术方案如下：
一种系统间共频信号差分的GNSS/INS组合定位方法，包括如下步骤：
步骤S1：将待测试设备与基准站分别放置在超短基线的两端，静止并同步采集10-60分钟的共频数据；
步骤S2：单历元估计DISB，根据DISB时间序列得到DISB的先验值；
步骤S3：接收并解码组合导航模块产生的IMU观测值、GNSS观测值，以及基准站播发的GNSS改正数；
步骤S4：利用事先校正得到的DISB信息，根据组合导航模块中内嵌的组合定位算法，进行GNSS/INS松组合或紧组合解算，待估参数滤波解算，模糊度参数固定与传递。
所述步骤S1包括：
步骤S11：将待测试设备与基准站布设为超短基线；
步骤S12：接收并解码待测试设备10-60分钟的GNSS共频观测值；
步骤S13：接收并解码基准站10-60分钟的GNSS共频观测值。
所述步骤S2包括：
步骤S21：预处理GNSS观测数据；
数据预处理包括但不限于异常数据的监测与剔除、周跳的探测与修复、监测站与参考站的单点定位、卫星截止高度角设置、数据截止信噪比设置、异步差分最大时差设置、参考卫星选取策略设置、按照高度角定权；
采用的高度角定权公式为：
Figure BDA0003487208020000031
其中，θi表示第i颗卫星的高度角，σ0是验前单位权中误差，σ是第i颗卫星观测值的验前标准差；
步骤S22：构建双差观测方程；
首先，已知GNSS原始观测值写作
Figure BDA0003487208020000032
Figure BDA0003487208020000033
其中，上下标s,K,j,r分别代表卫星、卫星所属系统、频率、以及接收机；φ,p分别表示载波相位与伪距观测值，单位为米；ρ,τ,ι分别表示卫地距、对流层延迟与电离层延迟；
Figure BDA0003487208020000034
其中
Figure BDA0003487208020000035
代表频率j；dtr,
Figure BDA0003487208020000036
则分别表示接收机钟差和卫星钟差；
Figure BDA0003487208020000037
表示接收机端相位硬件延迟与卫星端相位硬件延迟；
Figure BDA0003487208020000038
表示接收机端伪距硬件延迟与卫星端伪距硬件延迟；λj表示波长；
Figure BDA0003487208020000039
表示模糊度；
Figure BDA00034872080200000310
和
Figure BDA00034872080200000311
分别表示相位与伪距的观测噪声；参数的上下标反映了该参数受哪些因素的影响；
如果此时有一属于系统Q的卫星v，并且卫星v也播发频率j上的信号，那么卫星s与v之间的双差观测值则为：
Figure BDA00034872080200000312
Figure BDA00034872080200000313
为了推导的方便，这里以双系统单频短基线的情况为例进行推导，即忽略电离层与对流层延迟的影响，并且分别假设观测到了系统Q,K中的n,m颗卫星；那么单历元TRTK观测方程写为：
Figure BDA00034872080200000314
其中，
Figure BDA00034872080200000315
Figure BDA00034872080200000316
具有类似的结构，
Figure BDA00034872080200000317
为待估的坐标改正值；
步骤S23：估计参数浮点解；
对于式(4)，利用
Figure BDA0003487208020000041
进行等价变化，得到TRTK与LRTK之间的关系
Figure BDA0003487208020000042
其中
Figure BDA0003487208020000043
即为LRTK单历元的观测方程；TRTK相比于LRTK在引入两个新的观测值的同时，也引入了三个待估参数；由于DISBs参数在时域上具有非常好的稳定性，并且模糊度参数在没有周跳发生的情况下，在时域上也保持常数特性，因此，随着观测时间的累积，TRTK将比LRTK拥有越来越多的多余观测信息；但是由于TRTK中
Figure BDA0003487208020000044
与
Figure BDA0003487208020000045
成线性相关的关系，这两个参数将无法同时正确估计出来；为了能够正确计算出相位DISB的值，将
Figure BDA0003487208020000046
与
Figure BDA0003487208020000047
进行参数合并；用
Figure BDA0003487208020000048
来代替该合并之后的参数，则式(5)写为
Figure BDA0003487208020000049
最后利用最小二乘计算出每个历元中所有参数的浮点值；
步骤S24：获取DISB参数；
DISB参数按照定义指的是
Figure BDA0003487208020000051
和
Figure BDA0003487208020000052
但是由于在参数估计时
Figure BDA0003487208020000053
和
Figure BDA0003487208020000054
无法分离，使得
Figure BDA0003487208020000055
无法估计出来，因此通常的做法是将
Figure BDA0003487208020000056
的小数部分
Figure BDA0003487208020000057
看为等效的相位DISB；首先计算出每个历元中
Figure BDA0003487208020000058
与
Figure BDA0003487208020000059
的值，为了减少观测值噪声对DISB的影响，通过长时间取平均来获得最终的DISB。
所述步骤S3包括：
步骤S31：接收并解码IMU观测值，提供给S4；
步骤S32：接收并解码组合系统中GNSS观测值，提供给S4；
步骤S33：接收并解码基准站GNSS改正数，提供给S4。
所述步骤S4包括：
步骤S41：预处理IMU数据，依次进行异常数据的监测与剔除、单位转换、坐标转换；提供给S42；
同时，预处理GNSS观测数据，依次进行异常数据的监测与剔除、周跳的探测与修复、监测站与参考站的单点定位、卫星截止高度角设置、数据截止信噪比设置、异步差分最大时差设置、参考卫星选取策略设置、按照高度角定权；提供给S42；
步骤S42：对INS进行初始化：通过GNSS SPP解算出组合系统的大致位置、速度与姿态，设定滤波系统的初始方差协方差矩阵与过程噪声矩阵，进入S43；
步骤S43：进行INS机械编排，进入S44；
步骤S44：结合固定模糊度参数以构建GNSS/INS松/紧组合方程，包括：
根据空间向量的关系，得到：
Figure BDA00034872080200000510
其中下上标b,e分别代表IMU自身坐标系和地心地固坐标系，
Figure BDA00034872080200000511
代表e系下GNSS天线相位中心处的坐标，
Figure BDA00034872080200000512
代表e系下IMU质心处的坐标，
Figure BDA00034872080200000513
代表从b系到e系的坐标变换矩阵，lb代表在IMU坐标系下GNSS天线相位中心与IMU质心之间的空间距离，为常数；式(8)描述了理想情况下GNSS解算的位置与INS解算的位置之间的关系，由于实际中解算过程不可避免的将引入误差，因此对上式(8)进行误差扰动分析得到GNSS解算的位置误差与INS解算的位置误差的关系
Figure BDA00034872080200000514
其中δ表示误差，φ表示姿态的误差；那么基于TRTK的GNSS/INS紧组合方程写为
Figure BDA0003487208020000061
基于TRTK的GNSS/INS松组合方程写为
Figure BDA0003487208020000062
其中(～)代表观测值，
Figure BDA0003487208020000063
代表速度误差，bω,bf分别代表陀螺仪与加速度计的零偏，sω,sf分别代表陀螺仪与加速度计的比例因子；
进入步骤S45；
步骤S45:待估参数滤波解算，输出最终的位置、速度、姿态结算结果，同时把结算结果提供给S46，用于进一步循环迭代：
根据参数的性质分别建立不同的随机游走模型，利用kalman滤波对所有待估参数进行解算；
进入步骤S46；
步骤S46：模糊度参数固定与传递：
采用LAMBDA方法对模糊度参数进行固定，待模糊度正确固定后，得到固定模糊度；同时对固定的模糊度参数进行记录，用于后续历元，执行S44，直到所有数据解算完成。
综上所述，本发明系统间共频信号差分的GNSS/INS组合定位方法，能够挖掘多系统GNSS的信息，在不额外增强成本的情况下，为GNSS/INS组合系统创造新的观测信息，从而有效提高组合系统在复杂环境下的定位性能。
具体的，与现有技术相比，本发明具有以下优点：
现有的GNSS/INS组合定位算法，严重依赖于GNSS信息，当GNSS信号受到干扰或遮挡时，将无法抑制INS误差的累积，从而导致组合系统定位性能的下降。现有的GNSS/INS组合定位改善技术需要引入新的外部传感器，一方面增强了组合系统的硬件成本，另一方面也使得融合算法的实现难度进一步提高。本发明借助了DISB在时域上优越的稳定性，充分挖掘了多系统GNSS的信息，显著增加了组合系统可用信息的个数，能够有效的提升目前GNSS/INS组合系统在复杂环境中的定位性能。
附图说明
图1系统间共频信号差分的GNSS/INS组合定位方法流程示意图；
具体实施方式
下面将结合示意图对本发明的具体实施方式进行更详细的描述。根据下列描述和权利要求书，本发明的优点和特征将更清楚。
参考图1，本发明一优选实施例中，一种系统间共频信号差分的GNSS/INS组合定位方法包括：
步骤S1：将待测试设备(可用与待测试设备型号相同的设备替代)与基准站分别放置在超短基线的两端，静止并同步采集10-60分钟的共频数据；
具体的，步骤S1包括：
步骤S11：将待测试设备与基准站布设为超短基线；
步骤S12：接收并解码待测试设备10-60分钟的GNSS共频观测值；
步骤S13：接收并解码基准站10-60分钟的GNSS共频观测值；
步骤S2：单历元估计DISB，根据DISB时间序列得到DISB的先验值；具体的，步骤S2包括：
步骤S21：预处理GNSS观测数据；
数据预处理包括但不限于异常数据的监测与剔除、周跳的探测与修复、监测站与参考站的单点定位、卫星截止高度角设置、数据截止信噪比设置、异步差分最大时差设置、参考卫星选取策略设置、按照高度角定权。
优选的，采用的高度角定权公式为：
Figure BDA0003487208020000081
其中，θi表示第i颗卫星的高度角，σ0是验前单位权中误差，σ是第i颗卫星观测值的验前标准差。
步骤S22：构建双差观测方程；
首先，已知GNSS原始观测值可以写作
Figure BDA0003487208020000082
Figure BDA0003487208020000083
其中，上下标s,K,j,r分别代表卫星、卫星所属系统、频率、以及接收机；φ,p分别表示载波相位与伪距观测值，单位为米；ρ,τ,ι分别表示卫地距、对流层延迟与电离层延迟；
Figure BDA0003487208020000084
其中
Figure BDA0003487208020000085
代表频率j；dtr,
Figure BDA0003487208020000086
则分别表示接收机钟差和卫星钟差；
Figure BDA0003487208020000087
表示接收机端相位硬件延迟与卫星端相位硬件延迟；
Figure BDA0003487208020000088
表示接收机端伪距硬件延迟与卫星端伪距硬件延迟；λj表示波长；
Figure BDA0003487208020000089
表示模糊度；
Figure BDA00034872080200000810
和
Figure BDA00034872080200000811
分别表示相位与伪距的观测噪声。需要强调的是，参数的上下标反映了该参数受哪些因素的影响。
如果此时有一属于系统Q的卫星v，并且卫星v也播发频率j上的信号，那么卫星s与v之间的双差观测值则为：
Figure BDA00034872080200000812
Figure BDA00034872080200000813
为了推导的方便，我们这里以双系统单频短基线的情况为例进行推导，即忽略电离层与对流层延迟的影响，并且分别假设观测到了系统Q,K中的n,m颗卫星。那么单历元TRTK观测方程可以写为：
Figure BDA00034872080200000814
其中，
Figure BDA00034872080200000815
Figure BDA00034872080200000816
具有类似的结构，
Figure BDA00034872080200000817
为待估的坐标改正值。
步骤S23：估计参数浮点解；
对于上式(4)，可以利用
Figure BDA0003487208020000091
进行等价变化，得到TRTK与LRTK之间的关系
Figure BDA0003487208020000092
其中
Figure BDA0003487208020000093
即为LRTK单历元的观测方程。可以看到，TRTK相比于LRTK在引入两个新的观测值的同时，也引入了三个待估参数。貌似TRTK拥有更少的多余观测数，但是幸运的是，由于DISBs参数在时域上具有非常好的稳定性，并且模糊度参数在没有周跳发生的情况下，在时域上也保持常数特性。因此，随着观测时间的累积，TRTK将比LRTK拥有越来越多的多余观测信息。但是由于TRTK中
Figure BDA0003487208020000094
与
Figure BDA0003487208020000095
成线性相关的关系，因此这两个参数将无法同时正确估计出来。为了能够正确计算出相位DISB的值，我们不妨将
Figure BDA0003487208020000096
与
Figure BDA0003487208020000097
进行参数合并。用
Figure BDA0003487208020000098
来代替该合并之后的参数，则式(5)可以写为
Figure BDA0003487208020000099
最后利用最小二乘即可计算出每个历元中所有参数的浮点值。
步骤S24：获取DISB参数；
DISB参数按照定义指的是
Figure BDA0003487208020000101
和
Figure BDA0003487208020000102
但是由于在参数估计时
Figure BDA0003487208020000103
和
Figure BDA0003487208020000104
无法分离，使得
Figure BDA0003487208020000105
无法估计出来，因此通常的做法是将
Figure BDA0003487208020000106
的小数部分
Figure BDA0003487208020000107
看为等效的相位DISB。首先计算出每个历元中
Figure BDA0003487208020000108
与
Figure BDA0003487208020000109
的值，为了减少观测值噪声对DISB的影响，可以通过长时间取平均来获得最终的DISB。
步骤S3：接收并解码组合导航模块产生的IMU观测值、GNSS观测值，以及基准站播发的GNSS改正数；具体的，步骤S3包括：
步骤S31：接收并解码IMU观测值，提供给S4。
步骤S32：接收并解码组合系统中GNSS观测值，提供给S4。
步骤S33：接收并解码基准站GNSS改正数，提供给S4。
步骤S4：利用事先校正得到的DISB信息，根据组合导航模块中内嵌的组合定位算法，进行GNSS/INS松组合或紧组合解算，进入S45；具体的，步骤S4包括：
对GNSS、IMU数据进行数据预处理、INS初始化、INS机械编排、构建GNSS/INS松/紧组合方程、待估参数滤波解算、模糊度参数固定与传递。
步骤S41：预处理IMU数据，依次进行(包括但不限于)异常数据的监测与剔除、单位转换、坐标转换；提供给S42；
同时，预处理GNSS观测数据，依次进行(包括但不限于)异常数据的监测与剔除、周跳的探测与修复、监测站与参考站的单点定位、卫星截止高度角设置、数据截止信噪比设置、异步差分最大时差设置、参考卫星选取策略设置、按照高度角定权；提供给S42。
步骤S42：对INS进行初始化，进入S43：
通过GNSS SPP解算出组合系统的大致位置、速度与姿态，设定滤波系统的初始方差协方差矩阵与过程噪声矩阵。
步骤S43：进行INS机械编排，进入S44。该技术为GNSS/INS组合导航技术中的基础内容，可在现有文献及开源程序中找到，此处不再赘述。
步骤S44：结合S46固定模糊度参数以构建GNSS/INS松/紧组合方程，包括：
根据空间向量的关系，可以得到：
Figure BDA00034872080200001010
其中下上标b,e分别代表IMU自身坐标系和地心地固坐标系。
Figure BDA0003487208020000111
代表e系下GNSS天线相位中心处的坐标，
Figure BDA0003487208020000112
代表e系下IMU质心处的坐标，
Figure BDA0003487208020000113
代表从b系到e系的坐标变换矩阵，lb代表在IMU坐标系下GNSS天线相位中心与IMU质心之间的空间距离，可以看为常数。上式描述了理想情况下GNSS解算的位置与INS解算的位置之间的关系，由于实际中解算过程不可避免的将引入误差，因此对上式(8)进行误差扰动分析即可得到GNSS解算的位置误差与INS解算的位置误差的关系
Figure BDA0003487208020000114
其中δ表示误差，φ表示姿态的误差。那么基于TRTK的GNSS/INS紧组合方程可以写为
Figure BDA0003487208020000115
类似的，基于TRTK的GNSS/INS松组合方程可以写为
Figure BDA0003487208020000116
其中(～)代表观测值，
Figure BDA0003487208020000117
代表速度误差，bω,bf分别代表陀螺仪与加速度计的零偏，sω,sf分别代表陀螺仪与加速度计的比例因子。
步骤S45:待估参数滤波解算，输出最终的位置、速度、姿态结算结果，同时把结算结果提供给S46，用于进一步循环迭代：
根据参数的性质分别建立不同的随机游走模型，利用kalman滤波对所有待估参数进行解算。
步骤S46：模糊度参数固定与传递：
采用LAMBDA方法对模糊度参数进行固定，待模糊度正确固定后，得到固定模糊度。同时对固定的模糊度参数进行记录，用于后续历元，直到所有数据解算完成。
综上所述，本发明系统间共频信号差分的GNSS/INS组合定位方法，通过深入挖掘多系统GNSS的信息，为GNSS/INS组合系统创造新的观测信息，能够有效提高组合系统在复杂环境下的定位性能。
具体的，与现有技术相比，本发明具有以下优点：
现有的GNSS/INS组合定位算法，严重依赖于GNSS信息，当GNSS信号受到干扰或遮挡时，将无法抑制INS误差的累积，从而导致组合系统导航性能的下降。现有的GNSS/INS组合定位改善技术需要引入新的外部传感器，一方面增加了组合系统的硬件成本，另一方面也使得融合算法的实现难度进一步提高。本发明充分挖掘多系统GNSS的信息，利用DISB在时域上优越的稳定性，显著增加组合系统可用信息的个数，能够有效提升GNSS/INS组合系统在复杂环境中的定位性能。
上述仅为本发明的优选实施例而已，并不对本发明起到任何限制作用。任何所属技术领域的技术人员，在不脱离本发明的技术方案的范围内，对本发明揭露的技术方案和技术内容做任何形式的等同替换或修改等变动，均属未脱离本发明的技术方案的内容，仍属于本发明的保护范围之内。

Claims (5)

1.一种系统间共频信号差分的GNSS/INS组合定位方法，其特征在于，包括如下步骤：
步骤S1：将待测试设备与基准站分别放置在超短基线的两端，静止并同步采集10-60分钟的共频数据；
步骤S2：单历元估计DISB，根据DISB时间序列得到DISB的先验值；
步骤S3：接收并解码组合导航模块产生的IMU观测值、GNSS观测值，以及基准站播发的GNSS改正数；
步骤S4：利用事先校正得到的DISB信息，根据组合导航模块中内嵌的组合定位算法，进行GNSS/INS松组合或紧组合解算，待估参数滤波解算，模糊度参数固定与传递。
2.如权利要求1所述的系统间共频信号差分的GNSS/INS组合定位方法，其特征在于，所述步骤S1包括：
步骤S11：将待测试设备与基准站布设为超短基线；
步骤S12：接收并解码待测试设备10-60分钟的GNSS共频观测值；
步骤S13：接收并解码基准站10-60分钟的GNSS共频观测值。
3.如权利要求1所述的系统间共频信号差分的GNSS/INS组合定位方法，其特征在于，所述步骤S2包括：
步骤S21：预处理GNSS观测数据；
数据预处理包括但不限于异常数据的监测与剔除、周跳的探测与修复、监测站与参考站的单点定位、卫星截止高度角设置、数据截止信噪比设置、异步差分最大时差设置、参考卫星选取策略设置、按照高度角定权；
采用的高度角定权公式为：
Figure FDA0003487208010000011
其中，θi表示第i颗卫星的高度角，σ0是验前单位权中误差，σ是第i颗卫星观测值的验前标准差；
步骤S22：构建双差观测方程；
首先，已知GNSS原始观测值写作
Figure FDA0003487208010000012
Figure FDA0003487208010000013
其中，上下标s,K,j,r分别代表卫星、卫星所属系统、频率、以及接收机；φ,p分别表示载波相位与伪距观测值，单位为米；ρ,τ,ι分别表示卫地距、对流层延迟与电离层延迟；
Figure FDA0003487208010000021
其中
Figure FDA0003487208010000022
代表频率j；dtr,
Figure FDA0003487208010000023
则分别表示接收机钟差和卫星钟差；
Figure FDA0003487208010000024
表示接收机端相位硬件延迟与卫星端相位硬件延迟；
Figure FDA0003487208010000025
表示接收机端伪距硬件延迟与卫星端伪距硬件延迟；λj表示波长；
Figure FDA0003487208010000026
表示模糊度；
Figure FDA0003487208010000027
和
Figure FDA0003487208010000028
分别表示相位与伪距的观测噪声；参数的上下标反映了该参数受哪些因素的影响；
如果此时有一属于系统Q的卫星v，并且卫星v也播发频率j上的信号，那么卫星s与v之间的双差观测值则为：
Figure FDA0003487208010000029
Figure FDA00034872080100000210
为了推导的方便，这里以双系统单频短基线的情况为例进行推导，即忽略电离层与对流层延迟的影响，并且分别假设观测到了系统Q,K中的n,m颗卫星；那么单历元TRTK观测方程写为：
Figure FDA00034872080100000211
其中，
Figure FDA00034872080100000212
Figure FDA00034872080100000213
Figure FDA00034872080100000214
具有类似的结构，
Figure FDA00034872080100000215
为待估的坐标改正值；
步骤S23：估计参数浮点解；
对于式(4)，利用
Figure FDA00034872080100000216
进行等价变化，得到TRTK与LRTK之间的关系
Figure FDA0003487208010000031
其中
Figure FDA0003487208010000032
即为LRTK单历元的观测方程；TRTK相比于LRTK在引入两个新的观测值的同时，也引入了三个待估参数；由于DISBs参数在时域上具有非常好的稳定性，并且模糊度参数在没有周跳发生的情况下，在时域上也保持常数特性，因此，随着观测时间的累积，TRTK将比LRTK拥有越来越多的多余观测信息；但是由于TRTK中
Figure FDA0003487208010000033
与
Figure FDA0003487208010000034
成线性相关的关系，这两个参数将无法同时正确估计出来；为了能够正确计算出相位DISB的值，将
Figure FDA0003487208010000035
与
Figure FDA0003487208010000036
进行参数合并；用
Figure FDA0003487208010000037
来代替该合并之后的参数，则式(5)写为
Figure FDA0003487208010000038
最后利用最小二乘计算出每个历元中所有参数的浮点值；
步骤S24：获取DISB参数；
DISB参数按照定义指的是
Figure FDA0003487208010000039
和
Figure FDA00034872080100000310
但是由于在参数估计时
Figure FDA00034872080100000311
和
Figure FDA00034872080100000312
无法分离，使得
Figure FDA00034872080100000313
无法估计出来，因此通常的做法是将
Figure FDA00034872080100000314
的小数部分
Figure FDA00034872080100000315
看为等效的相位DISB；首先计算出每个历元中
Figure FDA00034872080100000316
与
Figure FDA00034872080100000317
的值，为了减少观测值噪声对DISB的影响，通过长时间取平均来获得最终的DISB。
4.如权利要求1所述的系统间共频信号差分的GNSS/INS组合定位方法，其特征在于，所述步骤S3包括：
步骤S31：接收并解码IMU观测值，提供给S4；
步骤S32：接收并解码组合系统中GNSS观测值，提供给S4；
步骤S33：接收并解码基准站GNSS改正数，提供给S4。
5.如权利要求1所述的系统间共频信号差分的GNSS/INS组合定位方法，其特征在于，所述步骤S4包括：
步骤S41：预处理IMU数据，依次进行异常数据的监测与剔除、单位转换、坐标转换；提供给S42；
同时，预处理GNSS观测数据，依次进行异常数据的监测与剔除、周跳的探测与修复、监测站与参考站的单点定位、卫星截止高度角设置、数据截止信噪比设置、异步差分最大时差设置、参考卫星选取策略设置、按照高度角定权；提供给S42；
步骤S42：对INS进行初始化：通过GNSS SPP解算出组合系统的大致位置、速度与姿态，设定滤波系统的初始方差协方差矩阵与过程噪声矩阵，进入S43；
步骤S43：进行INS机械编排，进入S44；
步骤S44：结合固定模糊度参数以构建GNSS/INS松/紧组合方程，包括：
根据空间向量的关系，得到：
Figure FDA0003487208010000041
其中下上标b,e分别代表IMU自身坐标系和地心地固坐标系，
Figure FDA0003487208010000042
代表e系下GNSS天线相位中心处的坐标，
Figure FDA0003487208010000043
代表e系下IMU质心处的坐标，
Figure FDA0003487208010000044
代表从b系到e系的坐标变换矩阵，lb代表在IMU坐标系下GNSS天线相位中心与IMU质心之间的空间距离，为常数；式(8)描述了理想情况下GNSS解算的位置与INS解算的位置之间的关系，由于实际中解算过程不可避免的将引入误差，因此对上式(8)进行误差扰动分析得到GNSS解算的位置误差与INS解算的位置误差的关系
Figure FDA0003487208010000045
其中δ表示误差，φ表示姿态的误差；那么基于TRTK的GNSS/INS紧组合方程写为
Figure FDA0003487208010000051
基于TRTK的GNSS/INS松组合方程写为
Figure FDA0003487208010000052
其中(～)代表观测值，
Figure FDA0003487208010000053
代表速度误差，bω,bf分别代表陀螺仪与加速度计的零偏，sω,sf分别代表陀螺仪与加速度计的比例因子；
进入步骤S45；
步骤S45:待估参数滤波解算，输出最终的位置、速度、姿态结算结果，同时把结算结果提供给S46，用于进一步循环迭代：
根据参数的性质分别建立不同的随机游走模型，利用kalman滤波对所有待估参数进行解算；
进入步骤S46；
步骤S46：模糊度参数固定与传递：
采用LAMBDA方法对模糊度参数进行固定，待模糊度正确固定后，得到固定模糊度；同时对固定的模糊度参数进行记录，于后续历元，执行S44，直到所有数据解算完成。
CN202210084870.6A 2022-01-25 2022-01-25 一种系统间共频信号差分的gnss/ins组合定位方法 Pending CN114488245A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210084870.6A CN114488245A (zh) 2022-01-25 2022-01-25 一种系统间共频信号差分的gnss/ins组合定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210084870.6A CN114488245A (zh) 2022-01-25 2022-01-25 一种系统间共频信号差分的gnss/ins组合定位方法

Publications (1)

Publication Number Publication Date
CN114488245A true CN114488245A (zh) 2022-05-13

Family

ID=81474761

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210084870.6A Pending CN114488245A (zh) 2022-01-25 2022-01-25 一种系统间共频信号差分的gnss/ins组合定位方法

Country Status (1)

Country Link
CN (1) CN114488245A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115166799A (zh) * 2022-07-11 2022-10-11 中国科学院精密测量科学与技术创新研究院 顾及硬件延迟时变特征的gnss精密单点定位方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115166799A (zh) * 2022-07-11 2022-10-11 中国科学院精密测量科学与技术创新研究院 顾及硬件延迟时变特征的gnss精密单点定位方法
CN115166799B (zh) * 2022-07-11 2023-07-21 中国科学院精密测量科学与技术创新研究院 顾及硬件延迟时变特征的gnss精密单点定位方法

Similar Documents

Publication Publication Date Title
RU2749667C1 (ru) Способ и система быстрого и точного позиционирования
US10281587B2 (en) Navigation satellite system positioning involving the generation of correction information
Leandro et al. RTX positioning: The next generation of cm-accurate real-time GNSS positioning
US9201147B2 (en) GNSS signal processing with regional augmentation network
US10012738B2 (en) Positioning method and positioning apparatus using satellite positioning system
US9671501B2 (en) Global navigation satellite systems (GNSS) positioning using precise satellite data
CN111751853B (zh) 一种gnss双频载波相位整周模糊度解算方法
Li et al. Review of PPP–RTK: Achievements, challenges, and opportunities
WO2019203679A1 (en) High performance positioning system.
EP3901667A1 (en) Gnss receiver adapted to produce, use, and communicate software-generated satellite signal data
CN110824505B (zh) Gnss卫星接收机的偏差估计方法及系统、定位方法及终端
Liu et al. Single-frequency BDS/GPS RTK with low-cost U-blox receivers
CN108535749B (zh) 基于cors的定位增强方法及系统、定位系统
Cao et al. Performance evaluation of integrated GPS/GIOVE precise point positioning
CN116243591A (zh) 融合UTC(k)和北斗广播星历的亚纳秒级授时方法
US10534087B1 (en) Differential vector phase locked loop GPS reception method
Liu et al. Assessing partial ambiguity resolution and WZTD-constraint multi-frequency RTK in an urban environment using new BDS signals
Roesler et al. Tightly coupled processing of precise point positioning (PPP) and INS data
Tolman et al. Absolute precise kinematic positioning with GPS and GLONASS
CN114488245A (zh) 一种系统间共频信号差分的gnss/ins组合定位方法
Liu et al. Real-time multi-constellation precise point positioning with integer ambiguity resolution
Kazmierski Performance of absolute real-time multi-GNSS kinematic positioning
Shin et al. Inertially aided precise point positioning
CN113156477B (zh) 一种基于Android智能手机的高精度RTK定位方法
CN103389502B (zh) 基于多个地面基站高精度确定载体加速度的方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination