一种基于GNSS的单频高精度定位方法技术

一种基于GNSS的单频高精度定位方法，其步骤如下：一：准备工作；二：初始历元位置解算；三：构建载波相位历元差分观测模型，求解相对位置变化量；四：进行模糊度信息与位置信息传递；五：进行模糊度信息与方差信息调整；六：构建模糊度与位置信息虚拟观测方程；七：GRAPHIC组合观测方程与虚拟观测方程共同进行最小二乘估计，解算最终定位结果。通过上述步骤，本发明专利技术使用载波相位历元差分辅助的GRAPHIC组合运动学定位算法，提高了运动学导航定位精度；采用精密星历能对低轨卫星进行高精度定轨，采用广播星历则能实现实时导航定位；该方法，不受导航定位过程中换星影响，可靠度高，成本低廉，数据处理简单，利于普及与应用。

A single frequency high precision positioning method based on GNSS

A high precision positioning method based on single frequency GNSS, which comprises the following steps: A: preparation; two: initial epoch position solution; three: Construction of epoch difference carrier phase observation model, solving the relative position change; four: the ambiguity of information and position information transfer; five: the ambiguity of information and variance information adjustment; six: the construction of fuzzy degree and position information of the virtual observation equation; seven: GRAPHIC combination observation equation and observation equation of virtual joint least squares estimation, calculation of the final results. Through the above steps, the invention uses carrier phase epoch difference GRAPHIC positioning algorithm combined kinematic auxiliary, improve the accuracy of kinematic navigation and positioning; the precise ephemeris for high precision orbit determination for low earth orbit satellites, using broadcast ephemeris can realize real-time navigation and positioning; this method is not affected by the change of satellite navigation and positioning in the process of impact high reliability, low cost, simple data processing, conducive to the popularization and application of.

【技术实现步骤摘要】

本专利技术提供一种基于GNSS(GlobalNavigationSatelliteSystem)的高精度单频定位方法，它涉及一种对飞行器或地面，海面运动载体利用GNSS进行绝对位置解算的方法，属于导航

技术介绍
全球定位系统(GPS，GlobalPositioningSystem)是美国建立的一种卫星导航系统,能够为用户提供实时的三维导航服务。伴随着GPS的广泛应用，许多科研机构开始注意其单个接收机在低轨卫星精密定轨，以及地面和海面，静止或运动载体的高精度定位方面的应用。目前，上述方面的研究可以分为两个方向。第一，采用双频GPS的精密定位；该方法充分利用GPS播发的多频信息，通过各种不同的线性组合构建出无电离层观测量，从而消除电离层延迟的影响，以提高导航精度。采用双频GPS进行的运动学后处理低轨卫星定轨精度在厘米级。但双频方法对接收机要求较高，数据处理复杂，计算量大。第二，采用单频GPS的精密定位；单频GPS定位只利用GPS卫星播发的一个频率上的信息，其只包含伪距信息和载波信息，采用电离层模型或GRAPHIC(GroupandPhaseIonosphericCalibration)消电离层。相较于双频GPS精密定位，单频方法的精度要低一些，其定位精度一般在分米级，但由于单频接收机造价低廉，数据量小，处理简单，使得单个单频接收机的高精度定位方法研究仍具有巨大的实用意义。如今，各国的卫星导航系统都在如火如荼的建设当中，除了已经成熟完善的美国的GPS系统，我国的北斗卫星导航系统(BeiDouNavigationSatelliteSystem)，俄罗斯的格洛纳斯(GLONASS，GLOBALNAVIGATIONSATELLITESYSTEM)以及欧洲的伽利略卫星导航系统(Galileosatellitenavigationsystem)都在建设中。可以预想未来众多卫星导航系统的建立完善会使得基于GNSS的精密导航定位服务研究拥有更广阔的平台，而本专利所提出的一种基于GNSS的单频高精度定位方法便是其中之一。在进行精密定位的过程中，电离层延迟误差是最大的误差来源之一。在双频GPS中利用电离层延迟一阶项与载波频率的关系，通过两个载波上信息的线性组合来消除电离层延迟的影响。在单频GPS中利用电离层延迟一阶项在伪距观测量与载波观测量中大小相等符号相反的性质，构建GRAPHIC(GroupandPhaseIonosphericCalibration)观测量消除电离层延迟的影响。而由于伪距观测量的噪声为半米左右，载波观测量的噪声为毫米级，导致了构建的GRAPHIC组合观测量的噪声约为伪距噪声的一半，可达分米级。大噪声的观测量严重限制了定位结果的精度。在不同的情况与目的下，对所得到的伪距观测量与载波相位观测量有许多处理方法。目前在导航定位方面主要有星间单差消除接收机卫星钟差、站间单差消除GPS卫星钟差以及历元单差消除载波相位模糊度，三种单差又可以进行组合从而得到更加复杂的双差观测量以及三差观测量。综上，基于GNSS的精密定位具有巨大的研究与应用价值，其中的采用单接收机单频GPS测量量的定位算法，其设备简单，间隔低廉，数据处理容易，具有更大的普及与实用价值。但由于其定位精度相比双频GPS较差，从而限制了其在高精度情景下的应用。本专利技术提出了一种基于GNSS的单频高精度定位方法，提高了单频GPS定位精度，从而实现了低轨卫星的高精度定轨、地面、水面以及空中的静止或运动载体高精度定位与导航。
技术实现思路
(一)专利技术目的：本专利技术采用单接收机的单频GPS接收机数据，以运动学GRAPHIC组合定位方法为基础，创新的在每个历元定位过程中引入该历元时刻位置信息虚拟观测量，而当前历元位置信息虚拟观测量是使用上一历元的绝对位置加上载波相位历元差分所求解的相对位置变化量得到。由于载波相位观测量的观测精度很高，其历元差分求解的相对位置变化量误差在厘米级，在传统的GRAPHIC组合运动学定轨中加入这些信息，增加了观测方程，而没有引入大的误差，可以有效的抑制传统GRAPHIC组合定位结果的噪声，因此对于定位精度有很大的提高。(二)技术方案本专利技术一种基于GNSS的单频高精度定位方法，其步骤如下：步骤一：准备工作首先给出伪距观测方程：P＝ρ+c(δtr-δts)+δion+δrel+εp(1)式(1)中，P伪距观测量，以米(m)为单位；ρ为GPS卫星到接收机的真实几何距离；c为真空中的光速；δtr为接收机钟差；δts为GPS卫星钟差；δion为电离层延迟；δrel为相对论效应引起的误差；εp为伪距观测噪声，其符合方差为的白噪声；给出载波相位观测方程：式(2)中，Φ为载波相位观测量(以米为单位)，为载波相位观测量(以周为单位)；λ为载波波长；ρ为GPS卫星到接收机的真实几何距离；c为真空中光速；δtr为接收机钟差；δts为GPS卫星钟差；N为模糊度；δion为电离层延迟；δrel为相对论效应引起的误差；为载波相位观测噪声，其符合方差为的白噪声；为了消除电离层误差，利用伪距观测方程(1)和载波相位观测方程(2)中电离层延迟一阶项，大小相等，符号相反的性质，得到以下消除电离层误差项的GRAPHIC组合观测方程为：上式中δts可以通过外部星历得到，为已知项；δrel可以通过建模得到并进行修正。合并已知项的观测方程为：Gc＝ρ+cδtr+λNG+εG(4)其中Gc为误差修正后的GRAPHIC组合观测量，NG为GRAPHIC组合模糊度，表达式如下：下文在不引起歧义的情况下，统一将NG简写为N，并称为模糊度；εG为GRAPHIC组合观测噪声，表达式如下：其噪声方差如下：式(7)是GRAPHIC组合观测量的噪声方差，其中D(·)为求方差算子(下文出现，不再说明)，可见其是由伪距噪声方差与载波噪声方差组合得到；伪距的测量噪声可达米级，而载波相位的测量噪声一般为毫米级；由上式可以看到GRAPHIC组合的噪声水平处于伪距噪声和载波相位噪声之间，为分米级，使用这种组合进行传统的GRAPHIC组合运动学定轨、定位和导航，结果的误差在分米级，难以满足高精度定位要求；给出载波相位历元差分观测方程：上式(8)中，k表示当前历元，k-1表示上一历元；其他符号与上文相同。其中，GPS卫星钟差项δts可通过外部星历得到；相对论引起的误差项δrel可通过建模得到；而当历元时间较短时，相邻历元电离层变化不大，残留项(δion,k-δion,k-1)可以忽略；因此得到整理后的载波相位历元差分观测方程：上式中ΔΦc为GPS接收机钟差，相对论效应修正后的载波相位历元差分观测量；其观测噪声方差为：从(10)中可以看出，载波相位历元差分观测量的噪声水平与载波相位观测量的水平相当；以上，给出了载波相位历元差分测量模型，该模型用于相邻历元的相对位置求解，其结果将传递到下一历元用于构建位置信息虚拟观测量；使用最小二乘估计对当前时刻接收机的位置和钟差信息进行估计；以下给出一种基于GNSS的单频高精度定位方法具体步骤；步骤二：初始历元位置解算由于初始历元是导航定位过程的第一个历元，其没有来自上一历元的先验信息，因此对于初始历元要进行单独处理，初始化它的位置信息与模糊度信息；构建GRAPHIC组合观测本文档来自技高网...

【技术保护点】
一种基于GNSS的单频高精度定位方法，其特征在于：其步骤如下：步骤一：准备工作首先给出伪距观测方程：P＝ρ+c(δtr‑δts)+δion+δrel+εp·········(1)式(1)中，P伪距观测量，以米(m)为单位；ρ为GPS卫星到接收机的真实几何距离；c为真空中的光速；δtr为接收机钟差；δts为GPS卫星钟差；δion为电离层延迟；δrel为相对论效应引起的误差；εp为伪距观测噪声，其符合方差为的白噪声；给出载波相位观测方程：式(2)中，Φ为载波相位观测量，以米为单位，为载波相位观测量，以周为单位；λ为载波波长；ρ为GPS卫星到接收机的真实几何距离；c为真空中光速；δtr为接收机钟差；δts为GPS卫星钟差；N为模糊度；δion为电离层延迟；δrel为相对论效应引起的误差；为载波相位观测噪声，其符合方差为的白噪声；为了消除电离层误差，利用伪距观测方程(1)和载波相位观测方程(2)中电离层延迟一阶项，大小相等，符号相反的性质，得到以下消除电离层误差项的GRAPHIC组合观测方程为：上式中δts通过外部星历得到，为已知项；δrel通过建模得到并进行修正，合并已知项的观测方程为：Gc＝ρ+cδtr+λNG+εG···········(4)其中Gc为误差修正后的GRAPHIC组合观测量，NG为GRAPHIC组合模糊度，表达式如下：下文在不引起歧义的情况下，统一将NG简写为N，并称为模糊度；εG为GRAPHIC组合观测噪声，表达式如下：其噪声方差如下：式(7)是GRAPHIC组合观测量的噪声方差，其中D(·)为求方差算子是由伪距噪声方差与载波噪声方差组合得到；伪距的测量噪声能达米级，而载波相位的测量噪声一般为毫米级；由上式能看到GRAPHIC组合的噪声水平处于伪距噪声和载波相位噪声之间，为分米级，使用这种组合进行传统的GRAPHIC组合运动学定轨、定位和导航，结果的误差在分米级，难以满足高精度定位要求；给出载波相位历元差分观测方程：上式(8)中，k表示当前历元，k‑1表示上一历元；其他符号与上文相同；其中，GPS卫星钟差项δts通过外部星历得到；相对论引起的误差项δrel通过建模得到；而当历元时间较短时，相邻历元电离层变化不大，残留项(δion,k‑δion,k‑1)忽略；因此得到整理后的载波相位历元差分观测方程：上式中ΔΦc为GPS接收机钟差，相对论效应修正后的载波相位历元差分观测量；其观测噪声方差为：从(10)中看出，载波相位历元差分观测量的噪声水平与载波相位观测量的水平相当；以上，给出了载波相位历元差分测量模型，该模型用于相邻历元的相对位置求解，其结果将传递到下一历元用于构建位置信息虚拟观测量；使用最小二乘估计对当前时刻接收机的位置和钟差信息进行估计；以下给出一种基于GNSS的单频高精度定位方法具体步骤；步骤二：初始历元位置解算由于初始历元是导航定位过程的第一个历元，其没有来自上一历元的先验信息，因此对于初始历元要进行单独处理，初始化它的位置信息与模糊度信息；构建GRAPHIC组合观测量计算值Z；由公式(4)，得知GRAPHIC组合观测方程是一个关于状态变量x＝(xr,yr,zr,δtr,N1…Np)T的非线性方程，其中模糊度下标表示当前历元共观测到的GPS卫星数目；GRAPHIC组合观测量计算值Z，计算如下：其中的关系式如下：其中，上式符号“～”表示估计值，为接收机的三个位置估计分量，(xs,ys,zs)T为GPS卫星的三个位置分量；采用卡尔曼滤波器一步估计状态变量x＝(xr,yr,zr,δtr,N1…Np)T；因为GRAPHIC组合观测方程是关于状态变量x＝(xr,yr,zr,δtr,N1…Np)T的非线性方程，因此要对其线性化，线性化过程如下：将GRAPHIC组合观测方程(4)，写成如下形式：Gc＝h(x)+εG············(13)在参考点处线性化：其中便是上面所求的GRAPHIC组合观测量的计算值Z；Δx为状态变量x的实际值与估计值之差；H模型观测量对状态量x的偏导数雅克比矩阵，有：认为每个GRAPHIC组合观测量是相互独立的，观测噪声阵有如下形式：状态变量初始值的先验协方差阵P‑如下：式(17)，是一个分块对角阵；其中为位置信息初始先验方差阵；Pt‑为接收机钟差初始先验方差阵；为模糊度方差阵，取值如下：其中，取值一般为伪距单点定位误差的一半，现取30米；Pt‑＝0··············(19)式(19)，表示初始历元钟差δtr的初始带入值是准确的，在后面估计出的各个历元接收机钟差项全都是相对于第一历元时刻的相对钟差，后边不再叙述；式(20)，模糊度方差值为一个极大值，这是因为初始历元并没有模糊度信息的先验值，所以如此处理表示方差的初始带入值是完全不可信的；给状态变量的估计值赋初值,其中位置状态...

【技术特征摘要】
1.一种基于GNSS的单频高精度定位方法，其特征在于：其步骤如下：步骤一：准备工作首先给出伪距观测方程：P＝ρ+c(δtr-δts)+δion+δrel+εp·········(1)式(1)中，P伪距观测量，以米(m)为单位；ρ为GPS卫星到接收机的真实几何距离；c为真空中的光速；δtr为接收机钟差；δts为GPS卫星钟差；δion为电离层延迟；δrel为相对论效应引起的误差；εp为伪距观测噪声，其符合方差为的白噪声；给出载波相位观测方程：式(2)中，Φ为载波相位观测量，以米为单位，为载波相位观测量，以周为单位；λ为载波波长；ρ为GPS卫星到接收机的真实几何距离；c为真空中光速；δtr为接收机钟差；δts为GPS卫星钟差；N为模糊度；δion为电离层延迟；δrel为相对论效应引起的误差；为载波相位观测噪声，其符合方差为的白噪声；为了消除电离层误差，利用伪距观测方程(1)和载波相位观测方程(2)中电离层延迟一阶项，大小相等，符号相反的性质，得到以下消除电离层误差项的GRAPHIC组合观测方程为：上式中δts通过外部星历得到，为已知项；δrel通过建模得到并进行修正，合并已知项的观测方程为：Gc＝ρ+cδtr+λNG+εG···········(4)其中Gc为误差修正后的GRAPHIC组合观测量，NG为GRAPHIC组合模糊度，表达式如下：下文在不引起歧义的情况下，统一将NG简写为N，并称为模糊度；εG为GRAPHIC组合观测噪声，表达式如下：其噪声方差如下：式(7)是GRAPHIC组合观测量的噪声方差，其中D(·)为求方差算子是由伪距噪声方差与载波噪声方差组合得到；伪距的测量噪声能达米级，而载波相位的测量噪声一般为毫米级；由上式能看到GRAPHIC组合的噪声水平处于伪距噪声和载波相位噪声之间，为分米级，使用这种组合进行传统的GRAPHIC组合运动学定轨、定位和导航，结果的误差在分米级，难以满足高精度定位要求；给出载波相位历元差分观测方程：上式(8)中，k表示当前历元，k-1表示上一历元；其他符号与上文相同；其中，GPS卫星钟差项δts通过外部星历得到；相对论引起的误差项δrel通过建模得到；而当历元时间较短时，相邻历元电离层变化不大，残留项(δion,k-δion,k-1)忽略；因此得到整理后的载波相位历元差分观测方程：上式中ΔΦc为GPS接收机钟差，相对论效应修正后的载波相位历元差分观测量；其观测噪声方差为：从(10)中看出，载波相位历元差分观测量的噪声水平与载波相位观测量的水平相当；以上，给出了载波相位历元差分测量模型，该模型用于相邻历元的相对位置求解，其结果将传递到下一历元用于构建位置信息虚拟观测量；使用最小二乘估计对当前时刻接收机的位置和钟差信息进行估计；以下给出一种基于GNSS的单频高精度定位方法具体步骤；步骤二：初始历元位置解算由于初始历元是导航定位过程的第一个历元，其没有来自上一历元的先验信息，因此对于初始历元要进行单独处理，初始化它的位置信息与模糊度信息；构建GRAPHIC组合观测量计算值Z；由公式(4)，得知GRAPHIC组合观测方程是一个关于状态变量x＝(xr,yr,zr,δtr,N1…Np)T的非线性方程，其中模糊度下标表示当前历元共观测到的GPS卫星数目；GRAPHIC组合观测量计算值Z，计算如下：其中的关系式如下：其中，上式符号“～”表示估计值，为接收机的三个位置估计分量，(xs,ys,zs)T为GPS卫星的三个位置分量；采用卡尔曼滤波器一步估计状态变量x＝(xr,yr,zr,δtr,N1…Np)T；因为GRAPHIC组合观测方程是关于状态变量x＝(xr,yr,zr,δtr,N1…Np)T的非线性方程，因此要对其线性化，线性化过程如下：将GRAPHIC组合观测方程(4)，写成如下形式：Gc＝h(x)+εG············(13)在参考点处线性化：其中便是上面所求的GRAPHIC组合观测量的计算值Z；Δx为状态变量x的实际值与估计值之差；H模型观测量对状态量x的偏导数雅克比矩阵，有：认为每个GRAPHIC组合观测量是相互独立的，观测噪声阵有如下形式：状态变量初始值的先验协方差阵P-如下：式(17)，是一个分块对角阵；其中为位置信息初始先验方差阵；Pt-为接收机钟差初始先验方差阵；为模糊度方差阵，取值如下：其中，取值一般为伪距单点定位误差的一半，现取30米；Pt-＝0··············(19)式(19)，表示初始历元钟差δtr的初始带入值是准确的，在后面估计出的各个历元接收机钟差项全都是相对于第一历元时刻的相对钟差，后边不再叙述；式(20)，模糊度方差值为一个极大值，这是因为初始历元并没有模糊度信息的先验值，所以如此处理表示方差的初始带入值是完全不可信的；给状态变量的估计值赋初值,其中位置状态变量的估计值赋值如下：其中，(xspp,yspp,zspp)T表示由伪距单点定位得出的接收机位置，单点定位在此不再叙述；钟差状态变量的赋值如下：上式与式(19)联立，其意义表示初始历元钟差为0，并且是准确的；以后历元的钟差估计值皆表示相对于初始历元的相对钟差；模糊度状态变量估计值的赋值如下：卡尔曼滤波，滤波系数如下：K＝P-HT(HP-H...

【专利技术属性】
技术研发人员：陈培，张键，孙秀聪，魏华波，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11