一种基于卡方检验的鲁棒自适应GNSS高精度定位方法技术

本发明专利技术涉及信号的测量、传输、定位技术领域，尤其涉及一种基于卡方检验的鲁棒自适应GNSS高精度定位方法。首先，在已知GNSS接收机k

【技术实现步骤摘要】
一种基于卡方检验的鲁棒自适应GNSS高精度定位方法


[0001]本专利技术涉及信号的测量、传输、定位
，尤其涉及一种基于卡方检验的鲁棒自适应GNSS高精度定位方法。

技术介绍

[0002]全球导航卫星系统(Global Navigation Satellite System，GNSS)精确单点定位(Precise Point Positioning，PPP)是一种基于GNSS卫星信号载波的高精度定位方法，具有覆盖范围广、长期精度高、成本低且不受基线长度约束的优势。GNSS PPP一般通过伪距和载波相位测量以及来自卫星的精确卫星时钟和轨道产品实现单点高精度定位。只需要一个接收器就可以提供静态毫米级和动态分米级定位。尽管PPP技术具有精度高、结构简单等优点，但由于受到多种因素，特别是异常测量值和模型精度影响，获得可靠和准确的定位结果仍然是一个难题。
[0003]目前，解决以上问题的方法主要是基于期望状态、方差分量比和预测残差的自适应因子生成方法，第一种方法一般利用状态模型的预测状态与估计状态之间的差值确定自适应因子，需要更多的数据点比状态特征的数量，才能得到一个准确的状态估计。第二种方法是基于Helmert的方差分量技术，通过投影状态方差与观测状态方差之比形成自适应因子，这类通过观测计算出预期状态的方法，也需要进行冗余观测，计算成本较高。第三种方法使用投影残差生成自适应因子，这样可以节省时间和精力，消除了收集重复数据点的需要，但是定位精度仍然可能受到观测噪声、建模精度和小周跳的影响，观测值的方差由随机模型决定，不精确的随机模型会降低卡尔曼滤波器的性能，甚至会导致滤波器发散。

技术实现思路

[0004]为克服传统的鲁棒卡尔曼滤波(RKF)所需观测值的方差由随机模型确定,不正确的随机模型会降低卡尔曼滤波器的性能的技术缺陷，本专利技术提供了一种基于卡方检验的鲁棒自适应GNSS高精度定位方法。
[0005]本专利技术提供了一种基于卡方检验的鲁棒自适应GNSS高精度定位方法，包括如下步骤：步骤一、在已知GNSS接收机k
‑
1时刻状态向量的前提下，经卡尔曼滤波的时间更新过程后得到k历元时刻的先验状态估计量和先验估计误差协方差矩阵，其子步骤为：GNSS接收机双频伪距和载波相位测量值组合形成组合测量值：（1）
（2）其中，是伪距P1和P2组合形成的电离层无关组合伪距测量值，是载波相位和组合形成的电离层无关组合载波测量值，和是GNSS播发的两个不同频率的频点，、分别是和对应的具体值，是卫星和接收机之间的几何距离，是光在真空中的传播速度，是接收机时钟误差，是对流层延迟，是电离层无关组合的整周模糊度，是卫星振荡器的初始相位，和分别是卫星和接收机误差，是多径效应，是测量中的综合噪声；将公式（1）和公式（2）中m颗卫星电离层无关组合伪距测量值和电离层无关组合载波测量值联合构建为测量向量：（3）其中，表示GNSS接收到的可用卫星数目；再将待估计的接收机位置、速度、接收机钟差、对流层延迟误差和整周模糊度组成状态向量：（4）其中，和=是接收器的位置和速度，是接收机钟差，是天顶对流层总延时误差，为整数模糊度；卡尔曼滤波的动态方程和测量方程为：（5）（6）其中，和分别表示k时刻的状态向量和测量向量，和分别为动态模型和测量模型；和分别为过程噪声和测量噪声；将公式（4）中的状态向量代入公式（5）经过时间更新得到先验估计量和对应的协方差矩阵分别为：（7）（8）其中，为先验状态估计值，为先验估计误差协方差矩阵；是转移矩阵，是过程噪声对应的协方差矩阵；步骤二、GNSS接收机在k时刻接收到m颗可用GNSS卫星的原始测量值并组合形成测量向量，再根据卫星定位解算方法求得GNSS定位结果；
步骤三、根据m颗可用卫星的测量向量解算得到的GNSS定位结果与公式（7）得到的先验状态估计量通过公式（9）计算得到卡方检验量：（9）步骤四、再根据公式（9）中的卡方检验量确定鲁棒自适应因子：（10）其中，表示鲁棒自适应因子，和是在自由度中预设的卡方检验临界值，是第j个GNSS测量值种类的个数；步骤五、根据(10)中的鲁棒自适应因子确定测量噪声的协方差矩阵为，当时，GNSS定位解算结果舍弃不用；当以及时，（11）其中为预设的GNSS观测环境理想情况下的测量噪声协方差矩阵；步骤六、根据公式（11）得到的计算卡尔曼增益矩阵：（12）步骤七、根据公式（12）得到的卡尔曼增益矩阵，计算时刻的后验状态估计量和对应的协方差矩阵：(13)(14)，步骤八，根据步骤七得到的后验状态估计量和对应的协方差矩阵，输出GNSS可靠高精度定位信息。
[0006]本专利技术提供的技术方案与现有技术相比具有如下优点：该方法可用于在城市地区运动模式下的GNSS精确单点定位（PPP），实现鲁棒自适应卡尔曼滤波和异常观测值的降权，有效提高GNSS PPP定位精度并缩短定位收敛时间。
附图说明
[0007]此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本专利技术的实施例，并与说明书一起用于解释本专利技术的原理。
[0008]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0009]图1为本专利技术所述一种基于卡方检验的鲁棒自适应GNSS高精度定位方法的流程图
具体实施方式
[0010]为了能够更清楚地理解本专利技术的上述目的、特征和优点，下面将对本专利技术的方案进行进一步描述。需要说明的是，在不冲突的情况下，本专利技术的实施例及实施例中的特征可以相互组合。
[0011]在下面的描述中阐述了很多具体细节以便于充分理解本专利技术，但本专利技术还可以采用其他不同于在此描述的方式来实施；显然，说明书中的实施例只是本专利技术的一部分实施例，而不是全部的实施例。
[0012]下面结合附图对本专利技术的具体实施例进行详细说明。
[0013]在一个实施例中，如图1所示，公开了一种基于卡方检验的鲁棒自适应GNSS高精度定位方法，包括如下步骤：步骤一、在已知GNSS接收机k
‑
1时刻状态向量的前提下，经卡尔曼滤波的时间更新过程后得到k历元时刻的先验状态估计量和先验估计误差协方差矩阵，其子步骤为：GNSS接收机双频伪距和载波相位测量值组合形成组合测量值：（1）（2）其中，是伪距P1和P2组合形成的电离层无关组合伪距测量值，是载波相位和组合形成的电离层无关组合载波测量值，和是GNSS播发的两个不同频率的频点，、分别是和对应的具体值，是卫星和接收机之间的几何距离，是光在真空中的传播速度，是接收机时钟误差，是对流层延迟，是电离层无关组合的整周模糊度，是卫星振荡器的初始相位，和分别是卫星和接收机误差，是多径效应，是测量中的综合噪声；将公式（1）和公式（2）中m颗卫星电离层无关组合伪距测量值和电离层无关组合载波测量值联合构建为测量向量：（3）其中，表示GNSS接收到的可用卫星数目；再将待估计的接收机本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于卡方检验的鲁棒自适应GNSS高精度定位方法，其特征在于，包括如下步骤：步骤一、在已知GNSS接收机k
‑
1时刻状态向量的前提下，经卡尔曼滤波的时间更新过程后得到k历元时刻的先验状态估计量和先验估计误差协方差矩阵，其子步骤为：GNSS接收机双频伪距和载波相位测量值组合形成组合测量值：（1）（2）其中，是伪距P1和P2组合形成的电离层无关组合伪距测量值，是载波相位和组合形成的电离层无关组合载波测量值，和是GNSS播发的两个不同频率的频点，、分别是和对应的具体值，是卫星和接收机之间的几何距离，是光在真空中的传播速度，是接收机时钟误差，是对流层延迟，是电离层无关组合的整周模糊度，是卫星振荡器的初始相位，和分别是卫星和接收机误差，是多径效应，是测量中的综合噪声；将公式（1）和公式（2）中m颗卫星电离层无关组合伪距测量值和电离层无关组合载波测量值联合构建为测量向量：（3）其中，表示GNSS接收到的可用卫星数目；再将待估计的接收机位置、速度、接收机钟差、对流层延迟误差和整周模糊度组成状态向量：（4）其中，和=是接收器的位置和速度，是接收机钟差，是天顶对流层总延时误差，为整数模糊度；卡尔曼滤波的动态方程和测量方程为：（5）（6）其中，和分别表示k时刻的状态向量和测量向量，和分别为动态模型和测量模型；和分别为过程噪声和测量噪声；将公式（4）中的状态向量代入公式（5）...

【专利技术属性】
技术研发人员：程俊兵，武杰，高云飞，
申请(专利权)人：太原理工大学，
类型：发明
国别省市：