CN107356947A - 基于单频导航卫星数据确定卫星差分伪距偏差的方法 - Google Patents
基于单频导航卫星数据确定卫星差分伪距偏差的方法 Download PDFInfo
- Publication number
- CN107356947A CN107356947A CN201710400859.5A CN201710400859A CN107356947A CN 107356947 A CN107356947 A CN 107356947A CN 201710400859 A CN201710400859 A CN 201710400859A CN 107356947 A CN107356947 A CN 107356947A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- msub
- mtd
- mtr
- 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.)
- Granted
Links
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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
-
- 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
Abstract
一种基于单频导航卫星数据确定卫星差分伪距偏差的方法,涉及卫星导航应用中卫星差分伪距偏差确定与修正技术,该方法包括以下步骤:A、单频GNSS原始观测数据的采集以及精密卫星轨道、钟差产品的获取;B、单频GNSS非组合伪距和相位观测方程的构建;C、单频非组合精密单点定位函数模型的构建;D、单频非组合精密单点定位卡尔曼滤波解算;E、电离层延迟建模以及导航卫星差分伪距偏差的确定;利用单频接收机实现了导航卫星差分伪距偏差参数的确定,可将现有导航卫星差分伪距偏差估计方法的硬件成本降低90﹪以上,同时,该方法设计合理简单,提高了效率。本设计不仅成本低,而且效率高。
Description
技术领域
本发明涉及卫星导航应用中卫星差分伪距偏差确定与修正技术,尤其涉及一种基于单频导航卫星数据确定卫星差分伪距偏差的方法,主要适用于降低成本、提高效率。
背景技术
卫星差分伪距偏差(DCB)是全球卫星导航系统(GNSS)测距信号在卫星硬件通道内的时延偏差,该参数是严重影响GNSS电离层监测与建模精度的系统误差,同时,它还是多频GNSS观测数据综合处理过程中必须消除的误差之一。导航卫星DCB参数的确定方法有硬件标定法与软件估计法两种,导航卫星通常会直接采用硬件标定法对其DCB参数进行初始测定。但是DCB受到硬件性能、外界环境等多种因素的影响通常会发生变化,实际中必须采用软件法精确估计DCB参数,以监测和修正其对卫星导航应用的影响。软件估计法通常是基于实测的双频GNSS数据精确确定卫星DCB参数,即用双频“无几何影响”观测值提取的电离层延迟量对TEC进行全球或区域建模,在建模过程中将DCB作为时不变参数进行同步估计。基于此方法,国际GNSS服务(IGS)组织联合全球多家电离层分析中心利用全球分布的IGS跟踪站双频观测数据连续计算并定期发布了GNSS卫星的差分伪距偏差产品。但是,实施该方法需要布设多台测地型双频GNSS接收机,用以连续地采集导航卫星的双频伪距和相位观测值,硬件投入较大。
发明内容
本发明的目的是克服现有技术中存在的成本高、效率低的缺陷与问题,提供一种成本低、效率高的基于单频导航卫星数据确定卫星差分伪距偏差的方法。
为实现以上目的,本发明的技术解决方案是:一种基于单频导航卫星数据确定卫星差分伪距偏差的方法,该方法包括以下步骤:
A、单频GNSS原始观测数据的采集以及精密卫星轨道、钟差产品的获取;
B、单频GNSS非组合伪距和相位观测方程的构建
经由线性化,原始非组合单频GNSS伪距和相位观测方程可表示为:
式(1)中,s和r分别为卫星、接收机;i为历元号;j为频率号;和分别为伪距与相位观测量;包含了与频率无关项:站星距、对流层延迟、观测噪声及若干非模型化误差;dtr,i和分别为接收机钟差和卫星钟差;dr,i和分别为接收机和卫星伪距硬件延迟;为第i频率上电离层斜延迟;为模糊度参数;
C、单频非组合精密单点定位函数模型的构建
a、基于各导航系统首个频率j=1上的观测数据和精密卫星轨道、钟差产品构建单频非组合精密单点定位函数模型;
b、消秩亏处理
合并和可得卫星差分伪距偏差参数
式(2)中,为消电离层组合卫星伪距硬件延迟参数,为首个频率上的卫星伪距硬件延迟;
合并两类接收机端的未知参数dtr,i和dr,1,得形为的有偏接收机钟差;
定义和为基准,则可得形为的可估电离层斜延迟,其中,为首历元有偏接收机钟差;
c、至少联合前两个历元的观测数据进行滤波初始化,建立满秩的单频非组合精密单点定位函数模型:
单频非组合精密单点定位函数模型的设计矩阵为:
式(4)中,第一列对应位置参数与对流层参数,第二列对应接收机钟差项,第三列对应模糊度参数,第四、五列分别对应第一、二个历元电离层参数,从第二个历元开始逐历元滤波解算,且接收机钟差参数可估;
初始化之后,假定当前历元即第二个历元之后观测到m颗卫星,联合所有卫星可得2m首个频率上的伪距P与相位Φ观测方程,具体形式如下:
式(5)中,Yi为第i个历元观测值向量,形式如
式(5)中,A,i为第i个历元状态转移矩阵,为待估参数向量,εy为未被模型化的噪声;
D、单频非组合精密单点定位卡尔曼滤波解算
从第二个历元开始逐历元卡尔曼滤波解算,过程描述如下:
时间预报:
式(7)中,和Di,i-1分别为待估参数一步时间预报值及其预报协方差阵,Φ为状态转移矩阵;
新息向量:
式(8)中,和为新息向量及其协方差阵,Y为观测值向量,A为单频非组合精密单点定位函数模型的设计矩阵,为先验单位权方差,Q为观测值噪声的协方差阵;
增益矩阵K为:
状态向量更新:
式(10)中,和Di即为参数X在第i历元的卡尔曼滤波值及协方差阵;
E、电离层延迟的建模以及导航卫星差分伪距偏差的确定
将由步骤D估计得到的作为本步骤的输入信息,进行电离层延迟的建模,估计电离层延迟模型系数,同步实现卫星差分伪距偏差参数的确定。
步骤A中,所述单频GNSS原始观测数据的采集是指:利用单频GNSS接收机EVK–M8T作为硬件平台,采集导航卫星的单频观测数据,采样率为30s;所述精密卫星轨道、钟差产品的获取是指:通过国际GNSS服务组织IGS网站获取精密卫星轨道、钟差产品。
与现有技术相比,本发明的有益效果为:
本发明一种基于单频导航卫星数据确定卫星差分伪距偏差的方法,相对于传统的导航卫星差分伪距偏差确定方法,在保证结果精度的前提下显著降低了(降低约90﹪)硬件的投入成本;同时,该方法与搭建支持该技术的软件平台也较为简单,对现有支持传统方法的成熟软件作较小的修改(如植入非组合PPP模块)即可实现。因此,本发明不仅成本低,而且效率高。
附图说明
图1是本发明基于单频导航卫星数据确定卫星差分伪距偏差的方法的流程图。
图2是单频测站(CUAU站)数据利用本发明和邻近双频测站(CUCC站)利用传统载波相位平滑方法分别解算的GPS卫星DCB与CODE的月产品差异比较图。
具体实施方式
以下结合附图说明和具体实施方式对本发明作进一步详细的说明。
参见图1,一种基于单频导航卫星数据确定卫星差分伪距偏差的方法,该方法包括以下步骤:
A、单频GNSS原始观测数据的采集以及精密卫星轨道、钟差产品的获取;
B、单频GNSS非组合伪距和相位观测方程的构建
经由线性化,原始非组合单频GNSS伪距和相位观测方程可表示为:
式(1)中,s和r分别为卫星、接收机;i为历元号;j为频率号;和分别为伪距与相位观测量;包含了与频率无关项:站星距、对流层延迟、观测噪声及若干非模型化误差;dtr,i和分别为接收机钟差和卫星钟差;dr,i和分别为接收机和卫星伪距硬件延迟;为第i频率上电离层斜延迟;为模糊度参数;
C、单频非组合精密单点定位函数模型的构建
a、基于各导航系统首个频率j=1上的观测数据和精密卫星轨道、钟差产品构建单频非组合精密单点定位函数模型;
b、消秩亏处理
合并和可得卫星差分伪距偏差参数
式(2)中,为消电离层组合卫星伪距硬件延迟参数,为首个频率上的卫星伪距硬件延迟;
合并两类接收机端的未知参数dtr,i和dr,1,得形为的有偏接收机钟差;
定义和为基准,则可得形为的可估电离层斜延迟,其中,为首历元有偏接收机钟差;
c、至少联合前两个历元的观测数据进行滤波初始化,建立满秩的单频非组合精密单点定位函数模型:
单频非组合精密单点定位函数模型的设计矩阵为:
式(4)中,第一列对应位置参数与对流层参数,第二列对应接收机钟差项,第三列对应模糊度参数,第四、五列分别对应第一、二个历元电离层参数,从第二个历元开始逐历元滤波解算,且接收机钟差参数可估;
初始化之后,假定当前历元即第二个历元之后观测到m颗卫星,联合所有卫星可得2m首个频率上的伪距P与相位Φ观测方程,具体形式如下:
式(5)中,Yi为第i个历元观测值向量,形式如
式(5)中,A,i为第i个历元状态转移矩阵,为待估参数向量,εy为未被模型化的噪声;
D、单频非组合精密单点定位卡尔曼滤波解算
从第二个历元开始逐历元卡尔曼滤波解算,过程描述如下:
时间预报:
式(7)中,和Di,i-1分别为待估参数一步时间预报值及其预报协方差阵,Φ为状态转移矩阵;
新息向量:
式(8)中,和为新息向量及其协方差阵,Y为观测值向量,A为单频非组合精密单点定位函数模型的设计矩阵,为先验单位权方差,Q为观测值噪声的协方差阵;
增益矩阵K为:
状态向量更新:
式(10)中,和Di即为参数X在第i历元的卡尔曼滤波值及协方差阵;
E、电离层延迟建模以及导航卫星差分伪距偏差的确定
将由步骤D估计得到的作为本步骤的输入信息,进行电离层延迟的建模,估计电离层延迟模型系数,同步实现卫星差分伪距偏差参数的确定。
步骤A中,所述单频GNSS原始观测数据的采集是指:利用单频GNSS接收机EVK–M8T作为硬件平台,采集导航卫星的单频观测数据,采样率为30s;所述精密卫星轨道、钟差产品的获取是指:通过国际GNSS服务组织IGS网站获取精密卫星轨道、钟差产品。
本发明的原理说明如下:
本设计提供一种基于单频GNSS数据的卫星差分伪距偏差确定方法,主要服务于GNSS卫星差分伪距偏差确定,该方法对传统单频精密单点定位算法模型进行了改进,构建了非组合的单频精密单点定位算法模型,实现了导航卫星差分伪距偏差的确定。通过设计合理的消秩亏策略,构建了非组合单频精密单点定位(SF-PPP)函数模型,并基于SF-PPP和单频GNSS数据提取了站星方向电离层延迟量,通过后续的电离层建模进而实现了导航卫星的差分伪距偏差估计。通过实验验证,基于本设计的方法利用成本较低的单频GNSS接收机数据解算的卫星码偏差精度与传统方法采用高成本的双频GNSS接收机数据解算的结果相当,从而为导航卫星差分伪距偏差的确定提供了一种低成本、高效率的解决方案。
实施例:
参见图1,一种基于单频导航卫星数据确定卫星差分伪距偏差的方法,该方法包括以下步骤:
A、单频GNSS原始观测数据的采集以及精密卫星轨道、钟差产品的获取
目前市场售价仅为249美元的UBLOXEVK–M8T接收机,可同时捕获三个导航系统调制于各自首个频率上的测距信号,且最多跟踪卫星数为25颗,利用单频GNSS接收机EVK–M8T作为硬件平台,采集导航卫星的单频观测数据,采样率为30s;精密卫星轨道、钟差产品采用国际GNSS服务组织IGS网站提供的多GNSS精密卫星产品;
B、单频GNSS非组合伪距和相位观测方程的构建
经由线性化,原始非组合单频GNSS伪距和相位观测方程可表示为:
式(1)中,s和r分别为卫星、接收机;i为历元号;j为频率号;和分别为伪距与相位观测量;包含了与频率无关项:站星距、对流层延迟、观测噪声及若干非模型化误差;dtr,i和分别为接收机钟差和卫星钟差;dr,i和分别为接收机和卫星伪距硬件延迟;为第i频率上电离层斜延迟;为模糊度参数;
C、单频非组合精密单点定位(SF–PPP)函数模型的构建
a、基于各导航系统首个频率j=1上的观测数据和IGS提供的精密卫星轨道、钟差产品构建单频非组合精密单点定位函数模型;
b、消秩亏处理
由式(1)可知,原始非组合伪距和相位观测方程中部分未知参数间相互不可分离,导致观测方程秩亏,不可分离的参数有:首个频率上的卫星伪距硬件延迟接收机伪距硬件延迟dr,1,接收机钟差dtr,i,电离层斜延迟模糊度参数和经IGS精密卫星钟差产品改正后引入的消电离层组合卫星伪距硬件延迟参数其消秩亏策略如下:
已知,IGS提供的精密卫星钟差是基于消电离层组合观测值解算的,因此其精密钟差产品中含有消电离层组合卫星伪距硬件延迟参数具体形式可表示为:
合并和可得卫星差分伪距偏差参数
式(2)中,为消电离层组合卫星伪距硬件延迟参数,为首个频率上的卫星伪距硬件延迟;
合并两类接收机端的未知参数dtr,i和dr,1,得形为的有偏接收机钟差;
和之间存在秩亏,定义和为基准,则可得形为的可估电离层斜延迟,其中,为首历元有偏接收机钟差;
c、为增加模型冗余度,需至少联合前两个历元的观测数据进行滤波初始化,建立满秩的单频非组合精密单点定位函数模型:
单频非组合精密单点定位函数模型的设计矩阵为:
式(4)中,第一列对应位置参数与对流层参数,第二列对应接收机钟差项,第三列对应模糊度参数,第四、五列分别对应第一、二个历元电离层参数,从第二个历元开始逐历元滤波解算,且接收机钟差参数可估;
初始化之后,假定当前历元即第二个历元之后观测到m颗卫星,联合所有卫星可得2m首个频率上的伪距P与相位Φ观测方程,具体形式如下:
式(5)中,Yi为第i个历元观测值向量,形式如
式(5)中,A,i为第i个历元状态转移矩阵,为待估参数向量,εy为未被模型化的噪声;
D、单频非组合精密单点定位卡尔曼滤波解算
从第二个历元开始逐历元卡尔曼滤波解算,过程描述如下:
时间预报:
式(7)中,和Di,i-1分别为待估参数一步时间预报值及其预报协方差阵,Φ为状态转移矩阵;
新息向量:
式(8)中,和为新息向量及其协方差阵,Y为观测值向量,A为单频非组合精密单点定位函数模型的设计矩阵,为先验单位权方差,Q为观测值噪声的协方差阵;
增益矩阵K为:
状态向量更新:
式(10)中,和Di即为参数X在第i历元的卡尔曼滤波值及协方差阵;
E、电离层延迟建模以及导航卫星差分伪距偏差的确定
将由步骤D估计得到的作为本步骤的输入信息,进行电离层延迟的建模,估计电离层延迟模型系数,同步实现卫星差分伪距偏差参数的确定。
利用本设计基于单频GNSS观测数据进行GPS卫星差分伪距偏差估计,估计精度如图2所示,由图2可以看出,利用相对较低成本的单频GPS接收机(单频接收机价格约为双频接收机价格的10﹪)实现了对GPS卫星差分伪距偏差的估计,其估计精度与传统基于双频观测数据的载波相位平滑法(CCL)的估计精度基本相当。
Claims (2)
1.一种基于单频导航卫星数据确定卫星差分伪距偏差的方法,其特征在于,该方法包括以下步骤:
A、单频GNSS原始观测数据的采集以及精密卫星轨道、钟差产品的获取;
B、单频GNSS非组合伪距和相位观测方程的构建
经由线性化,原始非组合单频GNSS伪距和相位观测方程可表示为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>dt</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>dt</mi>
<mrow>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>d</mi>
<mrow>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>s</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&phi;</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>dt</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>dt</mi>
<mrow>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>z</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>s</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式(1)中,s和r分别为卫星、接收机;i为历元号;j为频率号;和分别为伪距与相位观测量;包含了与频率无关项:站星距、对流层延迟、观测噪声及若干非模型化误差;dtr,i和分别为接收机钟差和卫星钟差;dr,i和分别为接收机和卫星伪距硬件延迟;为第i频率上电离层斜延迟;为模糊度参数;
C、单频非组合精密单点定位函数模型的构建
a、基于各导航系统首个频率j=1上的观测数据和精密卫星轨道、钟差产品构建单频非组合精密单点定位函数模型;
b、消秩亏处理
合并和可得卫星差分伪距偏差参数
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>d</mi>
<mrow>
<mi>G</mi>
<mi>F</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>d</mi>
<mrow>
<mi>I</mi>
<mi>F</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>d</mi>
<mrow>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>f</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mrow>
<msubsup>
<mi>f</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>f</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<msubsup>
<mi>d</mi>
<mrow>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</msubsup>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>f</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mrow>
<msubsup>
<mi>f</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>f</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<msubsup>
<mi>d</mi>
<mrow>
<mo>,</mo>
<mn>2</mn>
</mrow>
<mi>s</mi>
</msubsup>
<mo>-</mo>
<mfrac>
<mrow>
<msubsup>
<mi>f</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>f</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>f</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>f</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<msubsup>
<mi>d</mi>
<mrow>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>f</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mrow>
<msubsup>
<mi>f</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>f</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msubsup>
<mi>d</mi>
<mrow>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</msubsup>
<mi>-</mi>
<msubsup>
<mi>d</mi>
<mrow>
<mo>,</mo>
<mi>2</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式(2)中,为消电离层组合卫星伪距硬件延迟参数,为首个频率上的卫星伪距硬件延迟;
合并两类接收机端的未知参数dtr,i和dr,1,得形为的有偏接收机钟差;
定义和为基准,则可得形为的可估电离层斜延迟,其中,为首历元有偏接收机钟差;
c、至少联合前两个历元的观测数据进行滤波初始化,建立满秩的单频非组合精密单点定位函数模型:
单频非组合精密单点定位函数模型的设计矩阵为:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式(4)中,第一列对应位置参数与对流层参数,第二列对应接收机钟差项,第三列对应模糊度参数,第四、五列分别对应第一、二个历元电离层参数,从第二个历元开始逐历元滤波解算,且接收机钟差参数可估;
初始化之后,假定当前历元即第二个历元之后观测到m颗卫星,联合所有卫星可得2m首个频率上的伪距P与相位Φ观测方程,具体形式如下:
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>A</mi>
<mrow>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mrow>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<mi>y</mi>
</msub>
<mo>,</mo>
<mi>&epsiv;</mi>
<mo>&CenterDot;</mo>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>,</mo>
<msub>
<mi>Q</mi>
<mi>y</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
式(5)中,Yi为第i个历元观测值向量,形式如
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>p</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>i</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mn>1</mn>
</msubsup>
<mo>,</mo>
<mn>...</mn>
<msubsup>
<mi>p</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>i</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</msubsup>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&Phi;</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>&phi;</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>i</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mn>1</mn>
</msubsup>
<mo>,</mo>
<mn>...</mn>
<msubsup>
<mi>&phi;</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>i</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</msubsup>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
2
式(5)中,A,i为第i个历元状态转移矩阵,为待估参数向量,εy为未被模型化的噪声;
D、单频非组合精密单点定位卡尔曼滤波解算
从第二个历元开始逐历元卡尔曼滤波解算,过程描述如下:
时间预报:
<mrow>
<mfenced open = "" close = "}">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&Phi;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&Phi;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>&Phi;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mi>n</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
式(7)中,和Di,i-1分别为待估参数一步时间预报值及其预报协方差阵,Φ为状态转移矩阵;
新息向量:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>v</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>Y</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<msub>
<mover>
<mi>v</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
</msub>
<mo>=</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>A</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式(8)中,和为新息向量及其协方差阵,Y为观测值向量,A为单频非组合精密单点定位函数模型的设计矩阵,为先验单位权方差,Q为观测值噪声的协方差阵;
增益矩阵K为:
<mrow>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&Phi;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>A</mi>
<mi>j</mi>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>D</mi>
<msub>
<mover>
<mi>v</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
状态向量更新:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mover>
<mi>v</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>-</mo>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>-</mo>
<msubsup>
<mi>A</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>K</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>&sigma;</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>K</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
式(10)中,和Di即为参数X在第i历元的卡尔曼滤波值及协方差阵;
E、电离层延迟建模以及导航卫星差分伪距偏差的确定
将由步骤D估计得到的作为本步骤的输入信息,进行电离层延迟的建模,估计电离层延迟模型系数,同步实现卫星差分伪距偏差参数的确定。
2.根据权利要求1所述的一种基于单频导航卫星数据确定卫星差分伪距偏差的方法,其特征在于:步骤A中,所述单频GNSS原始观测数据的采集是指:利用单频GNSS接收机EVK–M8T作为硬件平台,采集导航卫星的单频观测数据,采样率为30s;所述精密卫星轨道、钟差产品的获取是指:通过国际GNSS服务组织IGS网站获取精密卫星轨道、钟差产品。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710400859.5A CN107356947B (zh) | 2017-05-31 | 2017-05-31 | 基于单频导航卫星数据确定卫星差分伪距偏差的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710400859.5A CN107356947B (zh) | 2017-05-31 | 2017-05-31 | 基于单频导航卫星数据确定卫星差分伪距偏差的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107356947A true CN107356947A (zh) | 2017-11-17 |
CN107356947B CN107356947B (zh) | 2019-06-18 |
Family
ID=60271986
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710400859.5A Active CN107356947B (zh) | 2017-05-31 | 2017-05-31 | 基于单频导航卫星数据确定卫星差分伪距偏差的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107356947B (zh) |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108873029A (zh) * | 2018-06-07 | 2018-11-23 | 南京航空航天大学 | 一种实现导航接收机钟差建模的方法 |
CN109061696A (zh) * | 2018-09-28 | 2018-12-21 | 中国人民解放军61540部队 | 一种确定导航卫星轨道和钟差的方法 |
CN109683186A (zh) * | 2018-12-20 | 2019-04-26 | 中国科学院国家授时中心 | 一种消除多卫星导航系统载波相位时间传递天跳变的方法 |
CN109709579A (zh) * | 2019-02-21 | 2019-05-03 | 哈尔滨工程大学 | 一种基于用户测距误差实时估计的gnss卫星星历故障检测方法 |
CN110161547A (zh) * | 2019-06-26 | 2019-08-23 | 哈尔滨工程大学 | 一种自适应电离层估计模型的中长基线模糊度解算方法 |
CN110208835A (zh) * | 2019-05-21 | 2019-09-06 | 哈尔滨工程大学 | 一种基于消电离层组合的跨系统紧组合差分定位方法 |
CN110764122A (zh) * | 2019-11-22 | 2020-02-07 | 武汉云图智信科技有限公司 | 一种单频gps接收机的精密单点定位方法 |
CN111175789A (zh) * | 2020-01-16 | 2020-05-19 | 中国民用航空总局第二研究所 | 地基增强系统的电离层异常监测方法、装置以及系统 |
CN111323796A (zh) * | 2020-03-18 | 2020-06-23 | 中国科学院国家空间科学中心 | 一种gnss接收机高采样钟差解算方法 |
CN111947667A (zh) * | 2020-06-24 | 2020-11-17 | 火眼位置数智科技服务有限公司 | 一种基于运动学和动力学组合的低轨卫星实时高精度定轨方法 |
CN111965673A (zh) * | 2020-06-24 | 2020-11-20 | 中山大学 | 基于多gnss的单频精密单点定位算法的时间频率传递方法 |
WO2020233158A1 (zh) * | 2019-05-22 | 2020-11-26 | 东南大学 | 一种基于智能手机的高精度单点定位方法与装置 |
CN112014860A (zh) * | 2020-07-20 | 2020-12-01 | 中国科学院空天信息创新研究院 | 一种基于北斗ppp-rtk的低轨卫星时空基准建立方法 |
CN112433234A (zh) * | 2020-11-19 | 2021-03-02 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种适用于中低纬区域gnss接收机的电离层tec实时估计方法 |
CN113126128A (zh) * | 2021-04-20 | 2021-07-16 | 东方红卫星移动通信有限公司 | 基于dcb改正北斗卫星导航系统频点精密单点定位方法 |
CN113204042A (zh) * | 2021-05-21 | 2021-08-03 | 北京交通大学 | 一种基于精密单点定位的多星座联合列车定位方法 |
CN113552606A (zh) * | 2020-04-23 | 2021-10-26 | 千寻位置网络有限公司 | 比特模糊度的确定方法 |
CN113671534A (zh) * | 2020-05-15 | 2021-11-19 | 华为技术有限公司 | 一种定位补偿方法、车载单元、介质及系统 |
CN113703017A (zh) * | 2020-05-22 | 2021-11-26 | 中国人民解放军战略支援部队信息工程大学 | 一种卫星天线相位中心偏差计算方法及装置 |
US20220171076A1 (en) * | 2020-11-30 | 2022-06-02 | Qualcomm Incorporated | Global navigation satellite system precise positioning engine with estimated ionosphere |
CN115267843A (zh) * | 2022-06-14 | 2022-11-01 | 中国科学院精密测量科学与技术创新研究院 | 多频多模gnss高频精密卫星钟差的实时非差估计方法 |
CN115390104A (zh) * | 2022-10-26 | 2022-11-25 | 中国人民解放军国防科技大学 | 一种导航卫星时延偏差建模方法 |
CN115712130A (zh) * | 2022-11-18 | 2023-02-24 | 中国科学院空天信息创新研究院 | 一种伪距校正方法、装置、电子设备及存储介质 |
CN116299617A (zh) * | 2023-02-16 | 2023-06-23 | 北京航空航天大学 | 一种gnss星地原子钟实时动态模型的建立方法 |
WO2023123147A1 (zh) * | 2021-12-30 | 2023-07-06 | 华为技术有限公司 | 一种数据预处理方法、数据预处理装置和芯片 |
CN116699663A (zh) * | 2023-08-08 | 2023-09-05 | 中国西安卫星测控中心 | 一种用于gnss观测融合定位的系统间偏差参数确定方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103558609A (zh) * | 2013-11-08 | 2014-02-05 | 中国科学院测量与地球物理研究所 | 全球卫星导航系统对流层天顶延迟修正的方法 |
CN104101888A (zh) * | 2014-05-28 | 2014-10-15 | 中国科学院光电研究院 | 一种适合于中国区域的gnss电离层延迟精确建模方法 |
CN106569242A (zh) * | 2016-11-04 | 2017-04-19 | 中国科学院测量与地球物理研究所 | 固定参考星的gnss单差处理方法 |
-
2017
- 2017-05-31 CN CN201710400859.5A patent/CN107356947B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103558609A (zh) * | 2013-11-08 | 2014-02-05 | 中国科学院测量与地球物理研究所 | 全球卫星导航系统对流层天顶延迟修正的方法 |
CN104101888A (zh) * | 2014-05-28 | 2014-10-15 | 中国科学院光电研究院 | 一种适合于中国区域的gnss电离层延迟精确建模方法 |
CN106569242A (zh) * | 2016-11-04 | 2017-04-19 | 中国科学院测量与地球物理研究所 | 固定参考星的gnss单差处理方法 |
Non-Patent Citations (2)
Title |
---|
张宝成 等: ""GPS接收机仪器偏差的短期时变特征提取与建模"", 《地球物理学报》 * |
张宝成 等: ""利用非组合精密单点定位技术确定斜向电离层中电子含量和站星差分码偏差"", 《测绘学报》 * |
Cited By (42)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108873029A (zh) * | 2018-06-07 | 2018-11-23 | 南京航空航天大学 | 一种实现导航接收机钟差建模的方法 |
CN108873029B (zh) * | 2018-06-07 | 2022-04-22 | 南京航空航天大学 | 一种实现导航接收机钟差建模的方法 |
CN109061696A (zh) * | 2018-09-28 | 2018-12-21 | 中国人民解放军61540部队 | 一种确定导航卫星轨道和钟差的方法 |
CN109683186A (zh) * | 2018-12-20 | 2019-04-26 | 中国科学院国家授时中心 | 一种消除多卫星导航系统载波相位时间传递天跳变的方法 |
CN109709579A (zh) * | 2019-02-21 | 2019-05-03 | 哈尔滨工程大学 | 一种基于用户测距误差实时估计的gnss卫星星历故障检测方法 |
CN109709579B (zh) * | 2019-02-21 | 2023-02-14 | 哈尔滨工程大学 | 一种基于用户测距误差实时估计的gnss卫星星历故障检测方法 |
CN110208835A (zh) * | 2019-05-21 | 2019-09-06 | 哈尔滨工程大学 | 一种基于消电离层组合的跨系统紧组合差分定位方法 |
CN110208835B (zh) * | 2019-05-21 | 2023-05-05 | 哈尔滨工程大学 | 一种基于消电离层组合的跨系统紧组合差分定位方法 |
US11709281B2 (en) | 2019-05-22 | 2023-07-25 | Southeast University | High-precision point positioning method and device based on smartphone |
WO2020233158A1 (zh) * | 2019-05-22 | 2020-11-26 | 东南大学 | 一种基于智能手机的高精度单点定位方法与装置 |
US20220155465A1 (en) * | 2019-05-22 | 2022-05-19 | Southeast University | High-precision point positioning method and device based on smartphone |
CN110161547A (zh) * | 2019-06-26 | 2019-08-23 | 哈尔滨工程大学 | 一种自适应电离层估计模型的中长基线模糊度解算方法 |
CN110161547B (zh) * | 2019-06-26 | 2022-09-16 | 哈尔滨工程大学 | 一种自适应电离层估计模型的中长基线模糊度解算方法 |
CN110764122A (zh) * | 2019-11-22 | 2020-02-07 | 武汉云图智信科技有限公司 | 一种单频gps接收机的精密单点定位方法 |
CN111175789A (zh) * | 2020-01-16 | 2020-05-19 | 中国民用航空总局第二研究所 | 地基增强系统的电离层异常监测方法、装置以及系统 |
CN111175789B (zh) * | 2020-01-16 | 2022-03-04 | 中国民用航空总局第二研究所 | 地基增强系统的电离层异常监测方法、装置以及系统 |
CN111323796B (zh) * | 2020-03-18 | 2021-11-09 | 中国科学院国家空间科学中心 | 一种gnss接收机高采样钟差解算方法 |
CN111323796A (zh) * | 2020-03-18 | 2020-06-23 | 中国科学院国家空间科学中心 | 一种gnss接收机高采样钟差解算方法 |
CN113552606A (zh) * | 2020-04-23 | 2021-10-26 | 千寻位置网络有限公司 | 比特模糊度的确定方法 |
CN113671534A (zh) * | 2020-05-15 | 2021-11-19 | 华为技术有限公司 | 一种定位补偿方法、车载单元、介质及系统 |
CN113703017B (zh) * | 2020-05-22 | 2023-06-13 | 中国人民解放军战略支援部队信息工程大学 | 一种卫星天线相位中心偏差计算方法及装置 |
CN113703017A (zh) * | 2020-05-22 | 2021-11-26 | 中国人民解放军战略支援部队信息工程大学 | 一种卫星天线相位中心偏差计算方法及装置 |
CN111965673A (zh) * | 2020-06-24 | 2020-11-20 | 中山大学 | 基于多gnss的单频精密单点定位算法的时间频率传递方法 |
CN111947667A (zh) * | 2020-06-24 | 2020-11-17 | 火眼位置数智科技服务有限公司 | 一种基于运动学和动力学组合的低轨卫星实时高精度定轨方法 |
CN111965673B (zh) * | 2020-06-24 | 2023-06-20 | 中山大学 | 基于多gnss的单频精密单点定位算法的时间频率传递方法 |
CN112014860B (zh) * | 2020-07-20 | 2023-07-14 | 中国科学院空天信息创新研究院 | 一种基于北斗ppp-rtk的低轨卫星时空基准建立方法 |
CN112014860A (zh) * | 2020-07-20 | 2020-12-01 | 中国科学院空天信息创新研究院 | 一种基于北斗ppp-rtk的低轨卫星时空基准建立方法 |
CN112433234A (zh) * | 2020-11-19 | 2021-03-02 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种适用于中低纬区域gnss接收机的电离层tec实时估计方法 |
US11550065B2 (en) * | 2020-11-30 | 2023-01-10 | Qualcomm Incorporated | Global navigation satellite system precise positioning engine with estimated ionosphere |
US20220171076A1 (en) * | 2020-11-30 | 2022-06-02 | Qualcomm Incorporated | Global navigation satellite system precise positioning engine with estimated ionosphere |
CN113126128A (zh) * | 2021-04-20 | 2021-07-16 | 东方红卫星移动通信有限公司 | 基于dcb改正北斗卫星导航系统频点精密单点定位方法 |
CN113204042A (zh) * | 2021-05-21 | 2021-08-03 | 北京交通大学 | 一种基于精密单点定位的多星座联合列车定位方法 |
WO2023123147A1 (zh) * | 2021-12-30 | 2023-07-06 | 华为技术有限公司 | 一种数据预处理方法、数据预处理装置和芯片 |
CN115267843A (zh) * | 2022-06-14 | 2022-11-01 | 中国科学院精密测量科学与技术创新研究院 | 多频多模gnss高频精密卫星钟差的实时非差估计方法 |
CN115390104A (zh) * | 2022-10-26 | 2022-11-25 | 中国人民解放军国防科技大学 | 一种导航卫星时延偏差建模方法 |
CN115390104B (zh) * | 2022-10-26 | 2023-04-07 | 中国人民解放军国防科技大学 | 一种导航卫星时延偏差建模方法 |
CN115712130A (zh) * | 2022-11-18 | 2023-02-24 | 中国科学院空天信息创新研究院 | 一种伪距校正方法、装置、电子设备及存储介质 |
CN115712130B (zh) * | 2022-11-18 | 2023-11-07 | 中国科学院空天信息创新研究院 | 一种伪距校正方法、装置、电子设备及存储介质 |
CN116299617A (zh) * | 2023-02-16 | 2023-06-23 | 北京航空航天大学 | 一种gnss星地原子钟实时动态模型的建立方法 |
CN116299617B (zh) * | 2023-02-16 | 2023-10-31 | 北京航空航天大学 | 一种gnss星地原子钟实时动态模型的建立方法 |
CN116699663A (zh) * | 2023-08-08 | 2023-09-05 | 中国西安卫星测控中心 | 一种用于gnss观测融合定位的系统间偏差参数确定方法 |
CN116699663B (zh) * | 2023-08-08 | 2023-10-17 | 中国西安卫星测控中心 | 一种用于gnss观测融合定位的系统间偏差参数确定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107356947B (zh) | 2019-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107356947A (zh) | 基于单频导航卫星数据确定卫星差分伪距偏差的方法 | |
CN104714244B (zh) | 一种基于抗差自适应Kalman滤波的多系统动态PPP解算方法 | |
CN102288978B (zh) | 一种cors基站周跳探测与修复方法 | |
CN108226985B (zh) | 基于精密单点定位的列车组合导航方法 | |
CN105158780B (zh) | 一种基于多种导航卫星可互换的导航定位方法 | |
CN104102822B (zh) | 一种多频gnss观测值随机特性建模方法 | |
CN109061696A (zh) | 一种确定导航卫星轨道和钟差的方法 | |
CN104656108B (zh) | 一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法 | |
CN103630914B (zh) | 一种gnss基线解算参考卫星选择方法 | |
CN110045407A (zh) | 一种分布式伪卫星/gnss优化定位方法 | |
CN108196284B (zh) | 一种进行星间单差模糊度固定的gnss网数据处理方法 | |
CN108254773A (zh) | 一种多gnss的实时钟差解算方法 | |
CN105182388B (zh) | 一种快速收敛的精密单点定位方法 | |
CN105467415B (zh) | 一种基于差分气压高度约束的小型无人机rtk相对定位方法 | |
CN103728643B (zh) | 附有宽巷约束的北斗三频网络rtk模糊度单历元固定方法 | |
CN106970398A (zh) | 顾及卫星遮挡条件的卫星可见性分析及星历预报方法 | |
CN104597465A (zh) | 一种提高gps与glonass组合精密单点定位收敛速度的方法 | |
CN106569242A (zh) | 固定参考星的gnss单差处理方法 | |
CN107728171A (zh) | 基于粒子滤波的gnss相位系统间偏差实时追踪和精密估计方法 | |
CN102998681A (zh) | 一种卫星导航系统的高频钟差估计方法 | |
CN107861131A (zh) | 一种斜路径电离层延迟的获取方法及系统 | |
CN104865587A (zh) | 基于geo星座粗定位的北斗接收机快速定位方法 | |
CN104614741A (zh) | 一种不受glonass码频间偏差影响的实时精密卫星钟差估计方法 | |
CN104035113A (zh) | 基于伪距的多模gnss接收机的可靠定位方法 | |
Banville et al. | Defining the basis of an integer-levelling procedure for estimating slant total electron content |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |