CN107193020B - 一种基于熵权法的bds卫星轨道位置插值方法 - Google Patents
一种基于熵权法的bds卫星轨道位置插值方法 Download PDFInfo
- Publication number
- CN107193020B CN107193020B CN201710568613.9A CN201710568613A CN107193020B CN 107193020 B CN107193020 B CN 107193020B CN 201710568613 A CN201710568613 A CN 201710568613A CN 107193020 B CN107193020 B CN 107193020B
- Authority
- CN
- China
- Prior art keywords
- weighted
- interpolation
- chebyshev
- time
- satellite orbit
- 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
- 238000000034 method Methods 0.000 title claims abstract description 46
- 239000011159 matrix material Substances 0.000 claims abstract description 20
- 230000005484 gravity Effects 0.000 claims description 21
- 239000000126 substance Substances 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 description 9
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- DMBHHRLKUKUOEG-UHFFFAOYSA-N diphenylamine Chemical compound C=1C=CC=CC=1NC1=CC=CC=C1 DMBHHRLKUKUOEG-UHFFFAOYSA-N 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
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/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
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (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
本发明提出一种基于熵权法的BDS卫星轨道位置插值方法,获取BDS卫星轨道系统精密星历中的t时刻的卫星轨道位置的拉格朗日插值和t时刻的卫星轨道位置的切比雪夫插值,将卫星轨道位置的拉格朗日插值和对应时刻的卫星轨道位置的切比雪夫插值组成待加权项目数据矩阵,计算待加权项目数据矩阵中各个方向的各待加权项的拉格朗日加权比重和切比雪夫加权比重,计算各个方向上待加权值的熵值,确定各个方向上待加权值的熵权,最终确定各个方向基于熵权法的轨道位置插值,本发明提高了单一差值方法的插值精度,解决了使用单一插值方法可能导致的误差较大情况。
Description
技术领域
本发明属于卫星轨道插值的技术领域,具体涉及一种基于熵权法的BDS卫星轨道位置插值方法。
背景技术
北斗导航系统(beidou navigation satellite system,BDS)是中国自主研制的全球卫星定位系统,是继美国全球定位系统(Global Positioning System,GPS)和俄罗斯的格洛纳斯系统(Global NavigationSatellite System,GLONASS)之后第三个成熟的卫星导航系统。北斗导航系统是全天候,全天时提供卫星导航信息的区域导航系统,此系统由定位于赤道上空的两颗地球同步卫星,地面中心站,用户终端三部分组成。该系统能提供四大功能:快速定位、实时导航、简短通信、精密授时。
在北斗卫星定位中,卫星位置的计算是一个关键环节。在北斗卫星数据后处理中,德国地学研究中心(Geo Forschungs Zentrum,GFZ)提供的精密星历数据时间间隔过长,是不同时间段等间隔的卫星位置坐标,而在实际的应用中,需要10s,5s甚至更密集的时间点卫星位置坐标。因此根据GFZ提供的精密星历数据选择正确的插值方法解算任意观测时刻的卫星位置具有重大的意义。本发明使用了拉格朗日插值法、切比雪夫插值法计算不同轨道卫星任一时刻的卫星位置,再利用熵权法完成两种插值方法的组合得出最终插值结果。
发明内容
针对现有技术的不足,本发明提出一种基于熵权法的BDS卫星轨道位置插值方法。
一种基于熵权法的BDS卫星轨道位置插值方法,包括以下步骤:
步骤1:获取BDS卫星轨道系统精密星历中的等间隔的已知历元时刻ti和tj,根据已知历元时刻ti和tj构造任意插值时刻的历元基函数li(t),其中,t为任意插值时刻,i,j为0到n的整数,n为待插值阶数;
步骤2:根据BDS卫星轨道系统精密星历中的已知等间隔历元位置(Xi,Yi,Zi)和任意插值时刻的历元基函数li(t),确定在t时刻的卫星轨道位置的拉格朗日插值[Xl(t),Yl(t),Zl(t)],其中,Xl(t)=(xl1,xl2,···xln),Yl(t)=(yl1,yl2,···yln),Zl(t)=(zl1,zl2,···zln);xli为X方向拉格朗日插值的第i个待加权值,yli为Y方向拉格朗日插值的第i个待加权值,zli为Z方向拉格朗日插值的第i个待加权值;
步骤3:获取观测时段[t0,t0+Δt],对观测时段[t0,t0+Δt]进行归一化处理,得到归一化后的观测时段τ∈[-1,1],其中,t0为观测历元首时刻,Δt为观测时长;
步骤4:构造卫星轨道位置插值的切比雪夫误差方程,根据归一化后的观测时段τ计算切比雪夫误差方程的系数矩阵A;
步骤5:根据卫星轨道位置插值的切比雪夫误差方程,以及切比雪夫误差方程的系数矩阵A,得到等间隔历元时刻多项式系数C,从而确定在t时刻的卫星轨道位置的切比雪夫插值[Xq(t),Yq(t),Zq(t)],其中,Xq(t)=(xq1,xq2,···xqn),Yq(t)=(yq1,yq2,···yqn),Zq(t)=(zq1,zq2,···zqn),xqi为X方向切比雪夫插值的第i个待加权值,yqi为Y方向切比雪夫插值的第i个待加权值,zqi为Z方向切比雪夫插值的第i个待加权值;
步骤7:计算待加权项目数据矩阵中各个方向的各待加权项的拉格朗日加权比重lli=[lxli,lyli,lzli]和切比雪夫加权比重lqi=[lxqi,lyqi,lzqi],其中,lxli为X方向第i个待加权值的拉格朗日加权比重,lyli为Y方向第i个待加权值的拉格朗日加权比重,lzli为Z方向第i个待加权值的拉格朗日加权比重,lxqi为X方向第i个待加权值的切比雪夫加权比重,lyqi为Y方向第i个待加权值的切比雪夫加权比重,lzqi为Z方向第i个待加权值的切比雪夫加权比重;
步骤8:根据待加权项目数据矩阵中各个方向的各待加权项的拉格朗日加权比重lli=[lxli,lyli,lzli]和切比雪夫加权比重lqi=[lxqi,lyqi,lzqi],计算各个方向上待加权值的熵值si=[sxi,syi,szi];
所述各个方向上待加权值的熵值sxi,syi,szi的计算公式如下所示:
sxi=-k(lxlilnlxli+lxqilnlxqi);
syi=-k(lylilnlyli+lyqilnlyqi);
szi=-k(lzlilnlzli+lzqilnlzqi);
步骤9:根据各个方向上待加权值的熵值si=[sxi,syi,szi],计算各个方向上待加权值的熵权ei=[exi,eyi,ezi];
所述各个方向上待加权值的熵权exi,eyi,ezi的计算公式如下所示:
步骤10:根据各个方向上待加权值的熵权ei=[exi,eyi,ezi],以及在t时刻的卫星轨道位置的拉格朗日插值[Xl(t),Yl(t),Zl(t)]和切比雪夫插值[Xq(t),Yq(t),Zq(t)],确定各个方向基于熵权法的轨道位置插值(wxi,wyi,wzi),其中,wxi为X方向上的基于熵权法的轨道位置插值,wyi为Y方向上的基于熵权法的轨道位置插值,wzi为Z方向上的基于熵权法的轨道位置插值。
所述各个方向基于熵权法的轨道位置插值(wxi,wyi,wzi)的计算公式如下所示:
wxi=exixqi+(1-exi)xli;
wyi=eyiyqi+(1-eyi)yli;
wzi=ezizqi+(1-ezi)zli。
所述构造的任意插值时刻的历元基函数li(t)如下所示:
所述t时刻的卫星轨道位置的拉格朗日插值[Xl(t),Yl(t),Zl(t)]的计算公式如下所示:
所述构造的卫星轨道位置插值的切比雪夫误差方程V=[Vx,Vy,Vz]的计算公式如下所示:
Vx=ACx-X;
Vy=ACy-Y;
Vz=ACz-Z;
本发明的有益效果:
本发明提出一种基于熵权法的BDS卫星轨道位置插值方法,本发明使用切比雪夫和拉格朗日插值法完成BDS轨道插值,再利用熵权法求得改进插值算法的轨道插值,本发明提高了单一差值方法的插值精度,解决了使用单一插值方法可能导致的误差较大情况。
附图说明
图1为本发明具体实施方式中基于熵权法的BDS卫星轨道位置插值方法的流程图;
图2为本发明具体实施方式中得到的不同插值情况误差比较图;
其中,(a)为实验一得到的不同历元的插值误差值;
(b)为实验二得到的不同历元的插值误差值;
(c)为实验三得到的不同历元的插值误差值;
(d)为实验四得到的不同历元的插值误差值;
(e)为实验五得到的不同历元的插值误差值;
(f)为实验六得到的不同历元的插值误差值。
具体实施方式
下面结合附图对本发明具体实施方式加以详细的说明。
一种基于熵权法的BDS卫星轨道位置插值方法,如图1所示,包括以下步骤:
步骤1:获取BDS卫星轨道系统精密星历中的等间隔的已知历元时刻ti和tj,根据已知历元时刻ti和tj构造任意插值时刻的历元基函数li(t)。
本实施方式中,构造的任意插值时刻的历元基函数li(t)如式(1)所示:
其中,t为任意插值时刻,i,j为0到n的整数,n为待插值阶数。
步骤2:根据BDS卫星轨道系统精密星历中的已知等间隔历元位置(Xi,Yi,Zi)和任意插值时刻的历元基函数li(t),确定在t时刻的卫星轨道位置的拉格朗日插值[Xl(t),Yl(t),Zl(t)]。
本实施方式中,t时刻的卫星轨道位置的拉格朗日插值[Xl(t),Yl(t),Zl(t)]的计算公式如式(2)、(3)和(4)所示:
其中,Xl(t)=(xl1,xl2,···xln),Yl(t)=(yl1,yl2,···yln),Zl(t)=(zl1,zl2,···zln);xli为X方向拉格朗日插值的第i个待加权值,yli为Y方向拉格朗日插值的第i个待加权值,zli为Z方向拉格朗日插值的第i个待加权值。
步骤3:获取观测时段[t0,t0+Δt],对观测时段[t0,t0+Δt]进行归一化处理,得到归一化后的观测时段τ∈[-1,1]。
本实施方式中,归一化后的观测时段τ的计算公式如式(5)所示:
其中,t0为观测历元首时刻,Δt为观测时长。
步骤4:构造卫星轨道位置插值的切比雪夫误差方程,根据归一化后的观测时段τ计算切比雪夫误差方程的系数矩阵A。
本实施方式中,构造的卫星轨道位置插值的切比雪夫误差方程V=[Vx,Vy,Vz]的计算公式如式(6)、(7)和(8)所示:
Vx=ACx-X (6)
Vy=ACy-Y (7)
Vz=ACz-Z (8)
根据公式(6)可以推出其展开公式(9)如下所示:
同理,可得公式(7)和公式(8)的展开公式。
切比雪夫误差方程的系数矩阵A如式(10)所示:
步骤5:根据卫星轨道位置插值的切比雪夫误差方程,以及切比雪夫误差方程的系数矩阵A,得到等间隔历元时刻多项式系数C,从而确定在t时刻的卫星轨道位置的切比雪夫插值[Xq(t),Yq(t),Zq(t)]。
本实施方式中,将公式(6)采用最小二乘概念整理为公式(11)所示:
Vx TPVx=min (11)
其中,P=I为权值。
可以得到等间隔历元时刻多项式系数Cx如式(12)所示:
Cx=(ATA)-1(ATX) (12)
同理,可以得到公式等间隔历元时刻多项式系数Cy和Cz。
本实施方式中,根据公式等间隔历元时刻多项式系数Cx、Cy和Cz确定在t时刻的卫星轨道位置的切比雪夫插值Xq(t),Yq(t)和Zq(t)如式(13)、(14)和(15)所示:
Xq(t)=(xq1,xq2,···xqn) (13)
Yq(t)=(yq1,yq2,···yqn) (14)
Zq(t)=(zq1,zq2,···zqn) (15)
其中,X方向切比雪夫插值的第i个待加权值xqi如式(16)所示:
Y方向切比雪夫插值的第i个待加权值yqi如式(17)所示:
Z方向切比雪夫插值的第i个待加权值zqi如式(18)所示:
步骤6:将卫星轨道位置的拉格朗日插值和对应时刻的卫星轨道位置的切比雪夫插值组成待加权项目数据矩阵
如式(19)、(20)和(21)所示:
步骤7:计算待加权项目数据矩阵中各个方向的各待加权项的拉格朗日加权比重lli=[lxli,lyli,lzli]和切比雪夫加权比重lqi=[lxqi,lyqi,lzqi]。
本实施方式中,X方向第i个待加权值的拉格朗日加权比重lxli的计算公式如式(22)所示:
同理,可得Y方向第i个待加权值的拉格朗日加权比重lyli,Z方向第i个待加权值的拉格朗日加权比重lzli。
X方向第i个待加权值的切比雪夫加权比重lxqi的计算公式如式(23)所示:
同理,可得Y方向第i个待加权值的切比雪夫加权比重lyqi,Z方向第i个待加权值的切比雪夫加权比重lzqi。
步骤8:根据待加权项目数据矩阵中各个方向的各待加权项的拉格朗日加权比重lli=[lxli,lyli,lzli]和切比雪夫加权比重lqi=[lxqi,lyqi,lzqi],计算各个方向上待加权值的熵值si=[sxi,syi,szi]。
本实施方式中,各个方向上待加权值的熵值sxi,syi,szi的计算公式如式(24)、(25)和(26)所示:
sxi=-k(lxlilnlxli+lxqilnlxqi) (24)
syi=-k(lylilnlyli+lyqilnlyqi) (25)
szi=-k(lzlilnlzli+lzqilnlzqi) (26)
步骤9:根据各个方向上待加权值的熵值si=[sxi,syi,szi],计算各个方向上待加权值的熵权ei=[exi,eyi,ezi]。
本实施方式中,各个方向上待加权值的熵权exi,eyi,ezi的计算公式如下如式(27)、(28)和(29)所示:
步骤10:根据各个方向上待加权值的熵权ei=[exi,eyi,ezi],以及在t时刻的卫星轨道位置的拉格朗日插值[Xl(t),Yl(t),Zl(t)]和切比雪夫插值[Xq(t),Yq(t),Zq(t)],确定各个方向基于熵权法的轨道位置插值(wxi,wyi,wzi),其中,wxi为X方向上的基于熵权法的轨道位置插值,wyi为Y方向上的基于熵权法的轨道位置插值,wzi为Z方向上的基于熵权法的轨道位置插值。
本实施方式中,各个方向基于熵权法的轨道位置插值(wxi,wyi,wzi)的计算公式如式(30)、(31)和(32)所示:
wxi=exixqi+(1-exi)xli; (30)
wyi=eyiyqi+(1-eyi)yli; (31)
wzi=ezizqi+(1-ezi)zli。 (32)
本发明实施方式中,采用了六组实验得到的不同历元的插值误差值,如图2所示。不同实验的误差值平均值如表1所示:
表1不同实验的误差值平均值
Claims (4)
1.一种基于熵权法的BDS卫星轨道位置插值方法,其特征在于,包括以下步骤:
步骤1:获取BDS卫星轨道系统精密星历中的等间隔的已知历元时刻ti和tj,根据已知历元时刻ti和tj构造任意插值时刻的历元基函数li(t),其中,t为任意插值时刻,i,j为0到n的整数,n为待插值阶数;
步骤2:根据BDS卫星轨道系统精密星历中的已知等间隔历元位置(Xi,Yi,Zi)和任意插值时刻的历元基函数li(t),确定在t时刻的卫星轨道位置的拉格朗日插值[Xl(t),Yl(t),Zl(t)],其中,Xl(t)=(xl1,xl2,…xln),Yl(t)=(yl1,yl2,…yln),Zl(t)=(zl1,zl2,…zln);xli为X方向拉格朗日插值的第i个待加权值,yli为Y方向拉格朗日插值的第i个待加权值,zli为Z方向拉格朗日插值的第i个待加权值;
步骤3:获取观测时段[t0,t0+Δt],对观测时段[t0,t0+Δt]进行归一化处理,得到归一化后的观测时段τ∈[-1,1],其中,t0为观测历元首时刻,Δt为观测时长;
步骤4:构造卫星轨道位置插值的切比雪夫误差方程,根据归一化后的观测时段τ计算切比雪夫误差方程的系数矩阵A;
步骤5:根据卫星轨道位置插值的切比雪夫误差方程,以及切比雪夫误差方程的系数矩阵A,得到等间隔历元时刻多项式系数C,从而确定在t时刻的卫星轨道位置的切比雪夫插值[Xq(t),Yq(t),Zq(t)],其中,Xq(t)=(xq1,xq2,…xqn),Yq(t)=(yq1,yq2,…yqn),Zq(t)=(zq1,zq2,…zqn),xqi为X方向切比雪夫插值的第i个待加权值,yqi为Y方向切比雪夫插值的第i个待加权值,zqi为Z方向切比雪夫插值的第i个待加权值;
步骤7:计算待加权项目数据矩阵中各个方向的各待加权项的拉格朗日加权比重lli=[lxli,lyli,lzli]和切比雪夫加权比重lqi=[lxqi,lyqi,lzqi],其中,lxli为X方向第i个待加权值的拉格朗日加权比重,lyli为Y方向第i个待加权值的拉格朗日加权比重,lzli为Z方向第i个待加权值的拉格朗日加权比重,lxqi为X方向第i个待加权值的切比雪夫加权比重,lyqi为Y方向第i个待加权值的切比雪夫加权比重,lzqi为Z方向第i个待加权值的切比雪夫加权比重;
步骤8:根据待加权项目数据矩阵中各个方向的各待加权项的拉格朗日加权比重lli=[lxli,lyli,lzli]和切比雪夫加权比重lqi=[lxqi,lyqi,lzqi],计算各个方向上待加权值的熵值si=[sxi,syi,szi];
所述各个方向上待加权值的熵值sxi,syi,szi的计算公式如下所示:
sxi=-k(lxlilnlxli+lxqilnlxqi);
syi=-k(lylilnlyli+lyqilnlyqi);
szi=-k(lzlilnlzli+lzqilnlzqi);
步骤9:根据各个方向上待加权值的熵值si=[sxi,syi,szi],计算各个方向上待加权值的熵权ei=[exi,eyi,ezi];
所述各个方向上待加权值的熵权exi,eyi,ezi的计算公式如下所示:
步骤10:根据各个方向上待加权值的熵权ei=[exi,eyi,ezi],以及在t时刻的卫星轨道位置的拉格朗日插值[Xl(t),Yl(t),Zl(t)]和切比雪夫插值[Xq(t),Yq(t),Zq(t)],确定各个方向基于熵权法的轨道位置插值(wxi,wyi,wzi),其中,wxi为X方向上的基于熵权法的轨道位置插值,wyi为Y方向上的基于熵权法的轨道位置插值,wzi为Z方向上的基于熵权法的轨道位置插值;
所述各个方向基于熵权法的轨道位置插值(wxi,wyi,wzi)的计算公式如下所示:
wxi=exixqi+(1-exi)xli;
wyi=eyiyqi+(1-eyi)yli;
wzi=ezizqi+(1-ezi)zli。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710568613.9A CN107193020B (zh) | 2017-07-13 | 2017-07-13 | 一种基于熵权法的bds卫星轨道位置插值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710568613.9A CN107193020B (zh) | 2017-07-13 | 2017-07-13 | 一种基于熵权法的bds卫星轨道位置插值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107193020A CN107193020A (zh) | 2017-09-22 |
CN107193020B true CN107193020B (zh) | 2020-04-14 |
Family
ID=59883545
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710568613.9A Active CN107193020B (zh) | 2017-07-13 | 2017-07-13 | 一种基于熵权法的bds卫星轨道位置插值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107193020B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111209523B (zh) * | 2020-01-06 | 2020-12-29 | 中国科学院紫金山天文台 | 适用于大偏心率轨道密集星历精密计算的快速处理方法 |
CN112346086B (zh) * | 2020-10-23 | 2024-03-15 | 西安空间无线电技术研究所 | 一种基于临近空间浮空平台的高效快速掩星预报方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101403790A (zh) * | 2008-11-13 | 2009-04-08 | 浙江师范大学 | 单频gps接收机的精密单点定位方法 |
CN101968542A (zh) * | 2010-09-29 | 2011-02-09 | 中国科学院国家天文台 | 一种利用地球站对月球探测器进行跟踪的方法 |
KR20110047485A (ko) * | 2009-10-30 | 2011-05-09 | 한국해양연구원 | 아이지에스 위성궤도 정보로부터 보간 데이터를 산출하는 장치 및 방법 |
CN102968552A (zh) * | 2012-10-26 | 2013-03-13 | 郑州威科姆科技股份有限公司 | 一种卫星轨道数据预估与修正方法 |
CN106772488A (zh) * | 2015-11-20 | 2017-05-31 | 河南理工大学 | 一种航空重力测量gps后处理方法 |
-
2017
- 2017-07-13 CN CN201710568613.9A patent/CN107193020B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101403790A (zh) * | 2008-11-13 | 2009-04-08 | 浙江师范大学 | 单频gps接收机的精密单点定位方法 |
KR20110047485A (ko) * | 2009-10-30 | 2011-05-09 | 한국해양연구원 | 아이지에스 위성궤도 정보로부터 보간 데이터를 산출하는 장치 및 방법 |
CN101968542A (zh) * | 2010-09-29 | 2011-02-09 | 中国科学院国家天文台 | 一种利用地球站对月球探测器进行跟踪的方法 |
CN102968552A (zh) * | 2012-10-26 | 2013-03-13 | 郑州威科姆科技股份有限公司 | 一种卫星轨道数据预估与修正方法 |
CN106772488A (zh) * | 2015-11-20 | 2017-05-31 | 河南理工大学 | 一种航空重力测量gps后处理方法 |
Non-Patent Citations (3)
Title |
---|
A new interpolation method based on satellite physical character in using IGS precise ephemeris;Liu Weiping,et al;《Geodesy and Geodynamics》;20141231;p29-33 * |
北斗三类卫星精密星历内插方法分析比较;汪威等;《全球定位系统》;20160430;第60-65页 * |
基于阶次组合的 GPS精密星历插值研究;王晓明等;《大地测量与地球动力学》;20110831;第103-106页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107193020A (zh) | 2017-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107797126B (zh) | 基于星型网络的bds/gps广播式网络rtk算法 | |
EP4206743A1 (en) | Method for vehicle positioning, related apparatus, device, and storage medium | |
CN106468774B (zh) | 一种应用于星基增强系统的星历星钟改正参数及空间信号完好性参数方法 | |
CN106125110B (zh) | 基于分区改正的星基增强系统定位精度提高方法 | |
Zhou et al. | Positioning accuracy assessment for the 4GEO/5IGSO/2MEO constellation of COMPASS | |
CN103823222B (zh) | 通过扩展sps轨道信息进行定位的方法和装置 | |
Jiao et al. | Improving BDS-2 and BDS-3 joint precise point positioning with time delay bias estimation | |
CN108363084A (zh) | 利用卫星定位的方法和装置、卫星导航接收机、存储介质 | |
CN110261879B (zh) | 广域地基增强位置服务的格网化虚拟参考站方法 | |
Li et al. | Review of PPP–RTK: Achievements, challenges, and opportunities | |
CN110376621A (zh) | 一种基于北斗三号B2b信号的卫星定位方法及装置 | |
CN111856534B (zh) | 智能终端的双模gnss载波精密单点定位方法及系统 | |
CN107656294B (zh) | 一种基于选星模板的多卫星导航系统选星方法 | |
CN102305935B (zh) | 多卫星导航星基增强系统提高定位精度的方法和系统 | |
CN107193020B (zh) | 一种基于熵权法的bds卫星轨道位置插值方法 | |
CN112285752B (zh) | 一种定位精度高的单点定位的方法及装置 | |
CN107861131A (zh) | 一种斜路径电离层延迟的获取方法及系统 | |
Zhang et al. | Impact of tracking station distribution structure on BeiDou satellite orbit determination | |
Zhang et al. | Impact of BDS-3 experimental satellites to BDS-2: Service area, precise products, precise positioning | |
CN114791614A (zh) | 钟差预报方法和装置 | |
CN114966760A (zh) | 一种电离层加权的非差非组合ppp-rtk技术实现方法 | |
CN102707296B (zh) | 一种用于单星座卫星导航系统的选星方法 | |
Xu et al. | A new ionosphere-free ambiguity resolution method for long-range baseline with GNSS triple-frequency signals | |
CN106772501B (zh) | 一种具有时间维持特性的伪距差分方法 | |
CN114721012A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20170922 Assignee: Fuxin Tianhui Geographic Information Technology Co.,Ltd. Assignor: LIAONING TECHNICAL University Contract record no.: X2023980048587 Denomination of invention: A BDS satellite orbit position interpolation method based on entropy weight method Granted publication date: 20200414 License type: Common License Record date: 20231127 |