CN107544082B - 北斗igso/meo卫星伪距码偏差一步建模法 - Google Patents
北斗igso/meo卫星伪距码偏差一步建模法 Download PDFInfo
- Publication number
- CN107544082B CN107544082B CN201610494536.2A CN201610494536A CN107544082B CN 107544082 B CN107544082 B CN 107544082B CN 201610494536 A CN201610494536 A CN 201610494536A CN 107544082 B CN107544082 B CN 107544082B
- Authority
- CN
- China
- Prior art keywords
- epoch
- satellite
- frequency
- normal equation
- specified
- 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.)
- Active
Links
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种北斗IGSO/MEO卫星伪距码偏差一步建模法,本发明通过对原始相位观测值进行周跳探测,确定当前历元各频率相位观测值是否设置为参考历元,基于MP组合观测值,忽略接收机端伪距码偏差,将各历元各测站MP组合观测值依卫星高度角定权与相同弧段的参考历元做差,准确消去未知的整周模糊度参数,从而将单颗卫星各频率的伪距码偏差加入相同法方程综合求解三阶多项式改正模型系数,最大程度地精确改正各颗卫星各频率高度角相关的伪距码偏差,提高北斗卫星的码相组合观测值的稳定性,对缩短模糊度收敛时间有重要意义。
Description
技术领域
本发明属于全球导航系统领域的定位定轨模糊度固定技术,特别涉及一种北斗IGSO/MEO卫星伪距码偏差一步建模法。
背景技术
模糊度收敛时间是目前卫星导航定位定轨技术的研究热点和难点问题。北斗系统IGSO/MEO卫星的伪距码观测值存在系统性偏差,见图1,影响载波相位模糊度解算的收敛速度和可靠性。Andre′Hauschild、Oliver Montenbruck等学者提出北斗IGSO/MEO卫星的码-相组合观测值存在与高度角相关的系统性误差[1-2]。
此外,Lambert Wanninger学者通过分析消除了电离层延迟、对流层延迟、距离相关的误差的MP(Multi-Path)组合观测值,采用第一步分弧段建模和第二步对相同类型卫星所有弧段增加MP组合观测值之和为零的约束,分别建立了IGSO/MEO卫星与高度角及频率相关的两套伪距码偏差离散型改正模型[3]。Lou YiDong等学者按照相同的建模方法提出三阶多项式的连续型改正模型[4]。由于每个弧段建模时对未知的模糊度参数采取的是时序平均值作为真值消去的策略,因此模糊度消除误差会被吸收到建模结果中。其次,不同测站不同弧段能追踪到的可视卫星不同,对同类型卫星所有参与建模的数据整体加上一个所有观测值之和为零的约束,也会影响北斗卫星的伪距码偏差建模精度,最后,每颗卫星的唯一性也会导致相同类型不同卫星的伪距码偏差存在差异。
文中涉及如下参考文献:
[1]Andre′Hauschild,Oliver Montenbruck,Jean-Marie Sleewaegen,LennardHuisman,Peter J.G.Teunissen,Characterization of Compass M-1 signals.GPS Solut(2012)16:117–126,DOI 10.1007/s10291-011-0210-3.
[2]Oliver Montenbruck,Andre′Hauschild,Peter Steigenberger,UrsHugentobler,Peter Teunissen,Shinichi Nakamura,Initial assessment of theCOMPASS/BeiDou-2 regional navigation satellite system.GPS Solut(2013)17:211–222,DOI 10.1007/s10291-012-0272-x.
[3]Lambert Wanninger,Susanne Beer,BeiDou satellite-induced codepseudorange variations:diagnosis and therapy.GPS Solut(2015)19:639–648,DOI10.1007/s10291-014-0423-3.
[4]Lou Y,Gong X,Gu S,et al.Assessment of code bias variations of BDStriple-frequency signals and their impacts on ambiguity resolution for longbaselines[J].GpsSolutions,2016:1-10.
发明内容
针对北斗系统长距离双差模糊度和精密单点定位非差模糊度难以快速有效固定的问题,本发明提出了一种可快速有效固定模糊度的北斗IGSO/MEO卫星伪距码偏差一步建模法。
本发明思路为:
有效顾及了模糊度消除策略的正确性,及数据处理过程中各历元依据高度角定权、各弧段依据历元数定权、不同测站依据先验值定权的灵活性,并且充分利用法方程综合解算模型系数,便于各弧段有效衔接,合理规避了所有卫星观测值之和为零的约束条件下存在的忽略不同弧段可视卫星存在差异的客观问题。
本发明技术方案如下:
一种北斗IGSO/MEO卫星伪距码偏差一步建模法,包括步骤:
步骤1,对指定单颗北斗IGSO/MEO卫星指定测站,根据指定测站的观测数据进行法方程叠加,本步骤进一步包括:
1.1利用IGS精密卫星轨道、精密卫星钟差产品和观测数据,采用精密单点定位法解算指定测站坐标;
1.2对观测数据逐历元进行周跳探测,将第一历元或相位观测值发生周跳的当前历元记为参考历元t0,计算t0对应的卫星高度角Elv(t0)及各频率MP组合观测值MPi(t0),将MPi(t0)表示为Elv(t0)的三阶多项式,i表示频率;然后,将t0的下一历元记为当前历元t1,执行子步骤1.3;
1.3计算当前历元t1对应的卫星高度角Elv(t1)及各频率MP组合观测值MPi(t1),将MPi(t1)表示为Elv(t1)的三阶多项式;
1.4将MPi(t0)和MPi(t1)的三阶多项式按权重做差得:
其中,a1、a2、a3为三阶多项式系数;和分别表示t1和t0观测数据的权重;
历元观测数据的权重根据历元对应的卫星高度角确定,若历元t对应的卫星高度角Elv(t)不大于30°,该历元观测数据的权重为4sin2(Elv(t));否则,权重为1;
1.5构建矩阵X=[a1a2 a3]T、X为待估系数矩阵,计算当前历元t1频率i的法方程要素和将和按频率分别对应累加至前一历元更新的法方程要素,得当前历元t1更新的法方程要素;
1.6对当前参考历元t0后的观测数据逐历元重复执行子步骤1.3~1.5,直至探测到相位观测值发生周跳的历元,将以当前参考历元t0为起点、该发生周跳历元的前一历元为终点间的历元弧段记为当前历元弧段,该发生周跳历元的前一历元更新的法方程要素即当前历元弧段的法方程要素,保存当前历元弧段法方程要素和历元数;然后,以该发生周跳的历元重置参考历元t0,采用子步骤1.3~1.5方法获得下一历元弧段法方程要素;
1.7将各历元弧段的法方程要素叠加得指定测站对应的法方程要素和Ni和Li表示指定测站频率i对应的法方程要素,Ni,n和Li,n表示频率i对应的第n个历元弧段的法方程要素,k表示历元弧段数,pi,n为频率i对应的第n个历元弧段的权值,表示频率i对应的第n个历元弧段的历元数;
步骤2,对指定单颗北斗IGSO/MEO卫星变换指定测站,采用步骤1方法获得变换后指定测站各频率对应的法方程要素;
步骤3,根据需求以及指定测站所处环境或指定测站历史观测数据质量,自由设定不同指定测站观测数据的权重,对指定单颗卫星各指定测站对应的法方程要素按频率分别进行加权叠加,获得各指定单颗卫星的法方程要素;
步骤4,在各频率下,分别根据指定单颗北斗IGSO/MEO卫星的法方程要素构建法方程,整体求解法方程得单颗北斗IGSO/MEO卫星各频率对应的待估系数矩阵X=[a1 a2 a3]T,即三阶多项式改正模型系数。
上述周跳探测采用TurboEdit法进行。
步骤3中,各指定测站观测数据的权重采用等权策略设定。
上述方法还包括对步骤4所获得的三阶多项式改正模型系数进行循环迭代剔除粗差,具体为:
在各频率下分别进行:
5.1根据各指定北斗IGSO/MEO卫星的三阶多项式改正模型系数,计算三阶多项式改正模型系数的均方差误差RMS,根据RMS设置阈值;
5.2逐指定北斗IGSO/MEO卫星逐历元,利用步骤4估计的三阶多项式改正模型系数以及当前历元和当前参考历元的卫星高度角,根据MPi组合观测值的三阶多项式估算当前历元和当前参考历元的组合观测值之差,记为
5.3将利用原始观测值计算的当前历元和当前参考历元的实际MPi组合观测值做差,所得差值记为ΔMPi;
5.4将与ΔMPi做差,若差值绝对值小于阈值,则当前历元观测数据为有效数据;否则,剔除当前历元观测数据;
5.5根据保留的历元观测数据,重新执行步骤1~3,求解待估系数矩阵X=[a1 a2a3]T。
上述阈值为3~10倍的RMS。
与现有技术相比,本发明具有以下优点和有益效果:
(1)本发明利用相同弧段内各历元与参考历元做差,准确消除未知的整周模糊度参数,提高模型精度。
(2)本发明将所有建模数据加入特定的法方程进行整体解算,求取三阶多项式改正模型系数,对不同弧段依据历元数定权,不同测站利用历史观测数据质量的经验值定权,进一步提高建模精度。
(3)本发明依据一步建模法建立了各颗IGSO/MEO卫星的高度角及频率相关的三阶多项式改正模型,提高各颗卫星码-相观测值的可靠性,有益于减少模糊度收敛时间。
附图说明
图1为IGSO/MEO卫星伪距系统性偏差示意图;
图2为不同IGSO卫星及所有IGSO卫星统一的伪距码偏差随高度角变化的曲线图,其中,图(a)、(b)、(c)分别对应频率B1、B2、B3;
图3为不同MEO卫星及所有MEO卫星统一的伪距码偏差随高度角变化的曲线图,其中,图(a)、(b)、(c)分别对应频率B1、B2、B3;
图4为原始伪距码观测值经三阶多项式改正模型改正前后的MP组合观测值及高度角的时间序列,其中,图(a)为改正后的MP组合观测值及高度角的时间序列,图(b)为改正前的MP组合观测值及高度角的时间序列,历元采样间隔为30s;
图5为本发明的具体流程图。
具体实施方式
下面将结合具体实施方式进一步说明本发明,见图5,本发明包括建模和粗差剔除两部分。
建模部分的具体步骤如下:
步骤1,对指定单颗北斗IGSO/MEO卫星指定测站,根据指定测站的观测数据进行法方程叠加。
本步骤进一步包括:
步骤1.1,利用IGS精密卫星轨道、精密卫星钟差产品和指定测站的观测数据,采用精密单点定位(Precise Point positioning,简称PPP)法解算指定测站坐标。
步骤1.2,对指定测站的观测数据进行周跳探测,若当前历元为第一历元或当前历元的相位观测值发生了周跳,将当前历元设置为参考历元t0,计算参考历元t0对应的卫星高度角Elv(t0)及各频率MP组合观测值MPi(t0),MPi(t0)表示参考历元t0频率i的MP组合观测值;若当前历元不为第一历元且当前历元的相位观测值未发生周跳,计算当前历元对应的卫星高度角Elv(t1)及各频率MP组合观测值MPi(t1),t1表示当前历元,i表示频率。本实施例中,采用已有的TurboEdit法进行周跳探测。
记i(i=1,2,3)分别表示北斗系统的三个频率B1、B2、B3,B1、B2、B3分别表示1561.098MHz、1207.14MHz、1268.52MHz。j、k为与i相应的另外两个频率,即当i=1时,j=2,k=3;当i=2时,j=1,k=3;当i=3时,j=1,k=2。λi、λj、λk分别为频率i、j、k对应的波长,分别为频率i、j、k对应的相位观测值;Pi表示频率i对应的伪距码观测值,bi,r、mMP、δMP分别表示频率i对应的卫星端码偏差、接收机端码偏差、MP组合观测值的多路径、MP组合观测值的观测噪声。
当前历元t1频率i的MP组合观测值MPi(t1)计算公式如下:
式(1)~(2)中:
Pi表示当前历元t1频率i对应的伪距码观测值;
mijk为比例系数,由式(2)计算获得。
实际应用中为计算方便,通常令j=i或k=i,即MP组合观测值中一个相位观测值与伪距码观测值频率相同。
步骤1.3,将子步骤1.2所得MP组合观测值表示为卫星高度角的三阶多项式,见式(3)~(4),三阶多项式系数记为am,am即三阶多项式中m次项的系数,m=0,1,2,3。
参考历元t0的MP组合观测值MPi(t0)所表示的三阶多项式如下:
MPi(t0)=a0+a1·Elv(t0)+a2·Elv(t0)2+a3·Elv(t0)3 (3)
当前历元t1的MP组合观测值MPi(t1)所表示的三阶多项式如下:
MPi(t1)=a0+a1·Elv(t1)+a2·Elv(t1)2+a3·Elv(t1)3 (4)
步骤1.4,根据各历元的卫星高度角确定各历元观测数据的权重,当历元t的卫星高度角Elv(t)≤30°时,历元t观测数据的权重为4sin2(Elv(t));否则,历元t观测数据的权重为1;t表示任意历元。
记参考历元t0的MP组合观测值MPi(t0)依据卫星高度角确定的权重为记当前历元t1的MP组合观测值MPi(t1)依据卫星高度角确定的权重为则将MPi(t0)和MPi(t1)的三阶多项式按照卫星高度角定权后做差,得式(5),如下:
步骤1.5,以矩阵形式表示公式(5)并构建法方程。
构建矩阵构建系数矩阵X=[a1 a2 a3]T,构建矩阵则式(5)等效于:
因待估系数矩阵X在不同历元是公共的,因此本发明采取可进行方程整体求解的法方程叠加原理进行系数矩阵估计,将式(6)两侧均乘以矩阵的转置矩阵,有:
记和即法方程要素,i表示北斗系统的三个不同频率。
以下一历元为当前历元,对下一历元重复执行子步骤1.2~1.5,并将所得历元的法方程要素累加到前一历元更新的法方程要素,直至探测到相位观测值发生周跳的当前历元,此时,以相位观测值发生周跳的当前历元重置参考历元t0,然后,基于重置的参考历元重新执行子步骤1.2~15。
重置参考历元时,保存当前历元弧段最后更新的法方程要素和历元数,然后对下一历元弧段按照上述方法计算法方程要素和历元数。
步骤1.6,重置参考历元时,则表示当前历元弧段结束,保存该当前历元弧段各历元法方程要素叠加的结果,即矩阵Ni,n和Li,n,将Ni,n和Li,n记为当前历元弧段的法方程要素;同时,保存当前历元弧段历元数numi,n,n=1,2,3......k,n代表历元弧段编号,k表示历元弧段数量。
对当前历元弧段所有历元按照步骤1.2~1.6处理至最后历元,逐历元按照相位观测值的周跳情况更新法方程叠加要素。
步骤1.7采用子步骤1.2~1.5计算各历元弧段的法方程要素及历元数,求取各历元弧段权值pi,n为频率i对应的第n个历元弧段权值。
步骤1.8,叠加各历元弧段的法方程要素,见公式(8)~(9):
其中,Ni和Li表示指定测站频率i对应的法方程要素,Ni,n和Li,n表示频率i对应的第n个历元弧段法方程要素,k表示历元弧段数,pi,n为频率i对应的第n个历元弧段权值,numi,n表示频率i对应的第n个历元弧段历元数。
步骤2,对指定单颗北斗IGSO/MEO卫星变换指定测站,采用步骤1方法获得变换后指定测站对应的法方程要素。
步骤3,根据需求以及指定测站所处环境或指定测站历史观测数据质量自由设定不同指定测站观测数据的权重。本发明测站均为IGS测站,因此选取等权策略。
根据各指定测站观测数据的权重,对指定单颗卫星的各指定测站对应的法方程要素进行加权叠加,获得各指定单颗卫星对应的法方程要素。
步骤4,根据指定单颗北斗IGSO/MEO卫星对应的法方程要素,逐卫星构建与频率对应的法方程,整体求解法方程得各频率对应的待估系数矩阵X=[a1 a2 a3]T,即三阶多项式改正模型系数,从而获得各单颗北斗IGSO/MEO卫星各频率对应的三阶多项式改正模型。
本步骤进一步包括:
步骤4.1,各指定单颗北斗IGSO/MEO卫星各频率均对应一组三阶多项式改正模型系数。
步骤4.2,各指定单颗北斗IGSO/MEO卫星各频率按照步骤1~2处理。
步骤4.2,依据已有的综合解方法,整体求解各组三阶多项式改正模型系数。
由上述步骤有NiXi=Li,Xi表示频率i对应的一组三阶多项式改正模型系数,两边同时左乘N的逆矩阵,即可求解出三阶多项式改正模型系数。求解公式如下:
Xi prn=Ni -1prnLi prn (10)
式(10)中:
i表示北斗系统的频率;
prn表示各颗北斗IGSO/MEO卫星,即prn依次为C06、C07、C08、C09、C10、C11、C12、C13、C14;
Ni prn和Li prn为北斗IGSO/MEO卫星prn频率i对应的法方程要素;
Xi prn为北斗IGSO/MEO卫星prn频率i对应的三阶多项式改正模型系数。
粗差剔除部分的具体步骤如下:
步骤5,采用已有的循环迭代剔除粗差法,确定最终改正模型。
在各频率下分别进行如下步骤:
步骤5.1,采用步骤1~4的方法整体估计各指定北斗IGSO/MEO卫星与卫星高度角相关的伪距码偏差三阶多项式改正模型系数,计算三阶多项式改正模型系数的均方差误差RMS。
步骤5.2,设置10倍RMS的限制,即阈值,针对各卫星逐历元迭代剔除部分粗差数据,具体方法如下:
首先,逐卫星逐历元,利用估计的三阶多项式改正模型系数以及当前历元和当前参考历元的卫星高度角,根据公式(3)和(4)估计当前历元和当前参考历元的组合观测值之差,记为接着,将利用原始观测值按照公式(1)计算的当前历元和当前参考历元的实际MPi组合观测值做差,所得差值记为ΔMPi;最后,将与ΔMPi做差,若该差值绝对值小于阈值,则视该当前历元观测数据为有效数据;反之,剔除该当前历元观测数据,即该当前历元观测数据不参与重新估计三阶多项式改正模型系数的估计过程。
步骤5.3,由步骤5.2逐卫星逐历元处理后,得到新的法方程要素和按照公式(10)解算法方程,重新估计待估系数,即更新各卫星各频率的三阶多项式改正模型系数。
实施例中,依次设置10、7、5、3倍RMS为阈值,并重复子步骤5.1~5.3,更新伪距码偏差三阶多项式改正模型系数,并最终确定以3倍RMS为阈值的粗差剔除数据策略的解算结果为最终改正模型系数。
作为优选,所述的IGSO/MEO卫星端伪距码偏差改正模型采用三阶多项式函数,为精密定位用户提供对伪距毫米级的改正量。
本发明基于多路径观测值的特性,通过动态设定参考历元,利用与相同弧段内的参考历元做历元间差分,准确地消除整周模糊度参数。同时,采取各历元依据高度角定权、不同弧段依据历元个数定权、不同测站依据先验值定权的定权策略,灵活合理地确定数据处理过程中的权重。此外,在建模求解改正模型系数时还利用了法方程综合解算,便于各颗卫星在不同测站不同弧段的有效衔接,成功规避了所有卫星观测值之和为零的约束条件中存在的忽略不同弧段可视卫星存在差异的客观问题,从而精确地得到各颗卫星伪距码偏差改正模型,保证其在模糊度固定中的实用性。
本文中所描述的具体实施例中的三阶多项式改正模型仅仅是对本发明思想做举例说明。本发明所述技术领域的技术人员可以对所描述的具体实施案例做各种各样的修改或补充或采用其他的函数模型替代,但并不会偏离本发明的思想或者超越所附权利要求书所定义的范围。
Claims (4)
1.一种北斗IGSO/MEO卫星伪距码偏差一步建模法,其特征是,包括:
步骤1,对指定单颗北斗IGSO/MEO卫星指定测站,根据指定测站的观测数据进行法方程叠加,本步骤进一步包括:
1.1利用IGS精密卫星轨道、精密卫星钟差产品和观测数据,采用精密单点定位法解算指定测站坐标;
1.2对观测数据逐历元进行周跳探测,将第一历元或相位观测值发生周跳的当前历元记为参考历元t0,计算t0对应的卫星高度角Elv(t0)及各频率MP组合观测值MPi(t0),将MPi(t0)表示为Elv(t0)的三阶多项式,i表示频率;然后,将t0的下一历元记为当前历元t1,执行子步骤1.3;
1.3计算当前历元t1对应的卫星高度角Elv(t1)及各频率MP组合观测值MPi(t1),将MPi(t1)表示为Elv(t1)的三阶多项式;
1.4将MPi(t0)和MPi(t1)的三阶多项式按权重做差得:
其中,a1、a2、a3为三阶多项式系数;和分别表示t1和t0观测数据的权重;
历元观测数据的权重根据历元对应的卫星高度角确定,若历元t对应的卫星高度角Elv(t)不大于30°,该历元观测数据的权重为4sin2(Elv(t));否则,权重为1;
1.5构建矩阵X=[a1 a2a3]T、X为待估系数矩阵,计算当前历元t1频率i的法方程要素将和按频率分别对应累加至前一历元更新的法方程要素,得当前历元t1更新的法方程要素;
1.6对当前参考历元t0后的观测数据逐历元重复执行子步骤1.3~1.5,直至探测到相位观测值发生周跳的历元,将以当前参考历元t0为起点、该发生周跳历元的前一历元为终点间的历元弧段记为当前历元弧段,该发生周跳历元的前一历元更新的法方程要素即当前历元弧段的法方程要素,保存当前历元弧段法方程要素和历元数;然后,以该发生周跳的历元重置参考历元t0,采用子步骤1.3~1.5方法获得下一历元弧段法方程要素;
1.7将各历元弧段的法方程要素叠加得指定测站对应的法方程要素和Ni和Li表示指定测站频率i对应的法方程要素,Ni,n和Li,n表示频率i对应的第n个历元弧段的法方程要素,k表示历元弧段数,pi,n为频率i对应的第n个历元弧段的权值,numi,n表示频率i对应的第n个历元弧段的历元数;
步骤2,对指定单颗北斗IGSO/MEO卫星变换指定测站,采用步骤1方法获得变换后指定测站各频率对应的法方程要素;
步骤3,根据需求以及指定测站所处环境或指定测站历史观测数据质量,自由设定不同指定测站观测数据的权重,对指定单颗卫星各指定测站对应的法方程要素按频率分别进行加权叠加,获得各指定单颗卫星的法方程要素;
步骤4,在各频率下,分别根据指定单颗北斗IGSO/MEO卫星的法方程要素构建法方程,整体求解法方程得单颗北斗IGSO/MEO卫星各频率对应的待估系数矩阵X=[a1 a2 a3]T,即三阶多项式改正模型系数。
2.如权利要求1所述的北斗IGSO/MEO卫星伪距码偏差一步建模法,其特征是:
所述的周跳探测采用TurboEdit法进行。
3.如权利要求1所述的北斗IGSO/MEO卫星伪距码偏差一步建模法,其特征是:
步骤3中,各指定测站观测数据的权重采用等权策略设定。
4.如权利要求1所述的北斗IGSO/MEO卫星伪距码偏差一步建模法,其特征是:
还包括对步骤4所获得的三阶多项式改正模型系数进行循环迭代剔除粗差,具体为:
在各频率下分别进行:
5.1根据各指定北斗IGSO/MEO卫星的三阶多项式改正模型系数,计算三阶多项式改正模型系数的均方差误差RMS,根据RMS设置阈值;所述的阈值为3~10倍的RMS;
5.2逐指定北斗IGSO/MEO卫星逐历元,利用步骤4估计的三阶多项式改正模型系数以及当前历元和当前参考历元的卫星高度角,根据MPi组合观测值的三阶多项式估算当前历元和当前参考历元的组合观测值之差,记为
5.3将利用原始观测值计算的当前历元和当前参考历元的实际MPi组合观测值做差,所得差值记为ΔMPi;
5.4将与ΔMPi做差,若差值绝对值小于阈值,则当前历元观测数据为有效数据;否则,剔除当前历元观测数据;
5.5根据保留的历元观测数据,重新执行步骤1~3,求解待估系数矩阵X=[a1 a2 a3]T。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610494536.2A CN107544082B (zh) | 2016-06-29 | 2016-06-29 | 北斗igso/meo卫星伪距码偏差一步建模法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610494536.2A CN107544082B (zh) | 2016-06-29 | 2016-06-29 | 北斗igso/meo卫星伪距码偏差一步建模法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107544082A CN107544082A (zh) | 2018-01-05 |
CN107544082B true CN107544082B (zh) | 2019-07-09 |
Family
ID=60962936
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610494536.2A Active CN107544082B (zh) | 2016-06-29 | 2016-06-29 | 北斗igso/meo卫星伪距码偏差一步建模法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107544082B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108710144B (zh) * | 2018-04-02 | 2021-02-02 | 武汉大学 | 零/短基线glonass相位通道间偏差校准方法及系统 |
CN110907960B (zh) * | 2018-09-17 | 2022-06-28 | 千寻位置网络有限公司 | 基于K-Means动态聚类分析的周跳探测方法及装置 |
CN111505687B (zh) * | 2020-04-17 | 2021-12-21 | 中国科学院国家授时中心 | 一种基于gps卫星导航系统的原始观测值粗差剔除方法 |
CN113552606B (zh) * | 2020-04-23 | 2022-07-26 | 千寻位置网络有限公司 | 比特模糊度的确定方法 |
CN111766616A (zh) * | 2020-06-15 | 2020-10-13 | 中国人民解放军61081部队 | 一种北斗二号时间传递星端多径误差改正方法 |
CN113933868B (zh) * | 2021-12-16 | 2022-03-22 | 中南大学 | 一种北斗二号meo卫星频率间卫星钟偏差的建模方法 |
CN114355421B (zh) * | 2021-12-22 | 2023-08-18 | 杭州电子科技大学 | 一种基于北斗卫星l4和cmc组合观测值洪水探测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217698A (zh) * | 2013-04-16 | 2013-07-24 | 武汉大学 | 一种基于北斗卫星导航系统三频信号模糊度的确定方法 |
CN103323867A (zh) * | 2013-07-05 | 2013-09-25 | 武汉大学 | 一种北斗二代geo卫星精密定轨系统偏差改进方法 |
CN104297773A (zh) * | 2014-02-27 | 2015-01-21 | 北京航天时代光电科技有限公司 | 一种高精度北斗三频sins深组合导航系统 |
CN104714244A (zh) * | 2015-03-31 | 2015-06-17 | 东南大学 | 一种基于抗差自适应Kalman滤波的多系统动态PPP解算方法 |
CN104749594A (zh) * | 2015-04-10 | 2015-07-01 | 武汉大学 | 一种gps双频非差周跳探测与修复方法及装置 |
CN105700000A (zh) * | 2016-04-13 | 2016-06-22 | 武汉大学 | 一种北斗导航接收机实时动态精密定位方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8456353B2 (en) * | 2011-01-14 | 2013-06-04 | Deere & Company | Method and system for determining clock corrections |
-
2016
- 2016-06-29 CN CN201610494536.2A patent/CN107544082B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217698A (zh) * | 2013-04-16 | 2013-07-24 | 武汉大学 | 一种基于北斗卫星导航系统三频信号模糊度的确定方法 |
CN103323867A (zh) * | 2013-07-05 | 2013-09-25 | 武汉大学 | 一种北斗二代geo卫星精密定轨系统偏差改进方法 |
CN104297773A (zh) * | 2014-02-27 | 2015-01-21 | 北京航天时代光电科技有限公司 | 一种高精度北斗三频sins深组合导航系统 |
CN104714244A (zh) * | 2015-03-31 | 2015-06-17 | 东南大学 | 一种基于抗差自适应Kalman滤波的多系统动态PPP解算方法 |
CN104749594A (zh) * | 2015-04-10 | 2015-07-01 | 武汉大学 | 一种gps双频非差周跳探测与修复方法及装置 |
CN105700000A (zh) * | 2016-04-13 | 2016-06-22 | 武汉大学 | 一种北斗导航接收机实时动态精密定位方法 |
Non-Patent Citations (2)
Title |
---|
北斗多频伪距单点定位精度分析;唐卫明,惠孟堂,崔健慧,鹿德凯;《测绘科学》;20160331;第41卷(第3期);全文 |
基于北斗CORS系统的民用导航差分定位系统实现与测试分析;潘国富,鲍志雄;《卫星导航定位与北斗系统应用2014——壮大北斗产业 创新位置服务》;20140930;全文 |
Also Published As
Publication number | Publication date |
---|---|
CN107544082A (zh) | 2018-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107544082B (zh) | 北斗igso/meo卫星伪距码偏差一步建模法 | |
CN111045034B (zh) | 基于广播星历的gnss多系统实时精密时间传递方法及系统 | |
CN107728171B (zh) | 基于粒子滤波的gnss相位系统间偏差实时追踪和精密估计方法 | |
CN107728180B (zh) | 一种基于多维粒子滤波偏差估计的gnss精密定位方法 | |
CN111381262B (zh) | 一种北斗三号精密单点定位参数优化方法和装置 | |
CN111965673B (zh) | 基于多gnss的单频精密单点定位算法的时间频率传递方法 | |
CN105629279B (zh) | 一种网络基准站间的宽巷模糊度固定方法 | |
CN104215977B (zh) | 一种基于卫星导航系统的精度评估方法及系统 | |
CN101770033A (zh) | 连续运行参考站系统站间整周模糊度网络固定方法 | |
CN102608633B (zh) | 一种卫星定位伪距差分方法 | |
CN107966722B (zh) | 一种gnss钟差解算方法 | |
RU2010102324A (ru) | Способ определения координат мобильного приемника спутниковой радионавигационной системы (срнс) | |
CN104730551B (zh) | 一种星地双基地差分干涉基线坐标以及形变量测量方法 | |
CN112285745B (zh) | 基于北斗三号卫星导航系统的三频模糊度固定方法及系统 | |
CN105044741A (zh) | 一种伪距相位综合广域差分改正值的求解方法 | |
CN104503223A (zh) | Gnss 三频高精度卫星钟差估计与服务方法 | |
CN111060928A (zh) | 一种基于双向测距的空间航天器时间同步实现方法 | |
CN103698790A (zh) | 北斗与gps双系统宽巷载波相位混频星间差分组合方法 | |
CN103675858B (zh) | 北斗系统b1与gps系统l1载波相位混频差分方法 | |
CN106292265A (zh) | 一种基于导航卫星的多地时间同步方法 | |
CN111044972A (zh) | 基于gnss精密时间同步的飞行器时差定位方法及系统 | |
CN105045086B (zh) | 一种改进的高精度rdss授时方法 | |
CN115906496A (zh) | 一种卫星钟差估计方法 | |
CN106291637B (zh) | 基于全伪距和部分伪距的定位方法 | |
CN105510936B (zh) | 星载gnss联合定轨方法及装置 |
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 |