CN107193020B - 一种基于熵权法的bds卫星轨道位置插值方法 - Google Patents

一种基于熵权法的bds卫星轨道位置插值方法 Download PDF

Info

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
Application number
CN201710568613.9A
Other languages
English (en)
Other versions
CN107193020A (zh
Inventor
王建敏
黄佳鹏
李亚博
董宏祥
谢栋平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Liaoning Technical University
Original Assignee
Liaoning Technical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Liaoning Technical University filed Critical Liaoning Technical University
Priority to CN201710568613.9A priority Critical patent/CN107193020B/zh
Publication of CN107193020A publication Critical patent/CN107193020A/zh
Application granted granted Critical
Publication of CN107193020B publication Critical patent/CN107193020B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/35Constructional details or hardware or software details of the signal processing chain
    • G01S19/37Hardware 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卫星轨道位置插值方法
技术领域
本发明属于卫星轨道插值的技术领域,具体涉及一种基于熵权法的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个待加权值;
步骤6:将卫星轨道位置的拉格朗日插值和对应时刻的卫星轨道位置的切比雪夫插值组成待加权项目数据矩阵
Figure GDA0002286884540000021
步骤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);
其中,
Figure GDA0002286884540000031
步骤9:根据各个方向上待加权值的熵值si=[sxi,syi,szi],计算各个方向上待加权值的熵权ei=[exi,eyi,ezi];
所述各个方向上待加权值的熵权exi,eyi,ezi的计算公式如下所示:
Figure GDA0002286884540000032
Figure GDA0002286884540000033
Figure GDA0002286884540000034
步骤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)如下所示:
Figure GDA0002286884540000035
所述t时刻的卫星轨道位置的拉格朗日插值[Xl(t),Yl(t),Zl(t)]的计算公式如下所示:
Figure GDA0002286884540000041
Figure GDA0002286884540000042
Figure GDA0002286884540000043
所述构造的卫星轨道位置插值的切比雪夫误差方程V=[Vx,Vy,Vz]的计算公式如下所示:
Vx=ACx-X;
Vy=ACy-Y;
Vz=ACz-Z;
其中,C=[Cx,Cy,Cz]为等间隔历元时刻多项式系数,
Figure GDA0002286884540000044
Figure GDA0002286884540000045
为切比雪夫误差方程的系数矩阵,
Figure GDA0002286884540000046
Figure GDA0002286884540000047
其中
Figure GDA0002286884540000048
τ∈[-1,1],n≥2。
本发明的有益效果:
本发明提出一种基于熵权法的BDS卫星轨道位置插值方法,本发明使用切比雪夫和拉格朗日插值法完成BDS轨道插值,再利用熵权法求得改进插值算法的轨道插值,本发明提高了单一差值方法的插值精度,解决了使用单一插值方法可能导致的误差较大情况。
附图说明
图1为本发明具体实施方式中基于熵权法的BDS卫星轨道位置插值方法的流程图;
图2为本发明具体实施方式中得到的不同插值情况误差比较图;
其中,(a)为实验一得到的不同历元的插值误差值;
(b)为实验二得到的不同历元的插值误差值;
(c)为实验三得到的不同历元的插值误差值;
(d)为实验四得到的不同历元的插值误差值;
(e)为实验五得到的不同历元的插值误差值;
(f)为实验六得到的不同历元的插值误差值。
具体实施方式
下面结合附图对本发明具体实施方式加以详细的说明。
一种基于熵权法的BDS卫星轨道位置插值方法,如图1所示,包括以下步骤:
步骤1:获取BDS卫星轨道系统精密星历中的等间隔的已知历元时刻ti和tj,根据已知历元时刻ti和tj构造任意插值时刻的历元基函数li(t)。
本实施方式中,构造的任意插值时刻的历元基函数li(t)如式(1)所示:
Figure GDA0002286884540000051
其中,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)所示:
Figure GDA0002286884540000052
Figure GDA0002286884540000053
Figure GDA0002286884540000054
其中,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)所示:
Figure GDA0002286884540000061
其中,t0为观测历元首时刻,Δt为观测时长。
步骤4:构造卫星轨道位置插值的切比雪夫误差方程,根据归一化后的观测时段τ计算切比雪夫误差方程的系数矩阵A。
本实施方式中,构造的卫星轨道位置插值的切比雪夫误差方程V=[Vx,Vy,Vz]的计算公式如式(6)、(7)和(8)所示:
Vx=ACx-X (6)
Vy=ACy-Y (7)
Vz=ACz-Z (8)
其中,C=[Cx,Cy,Cz]为等间隔历元时刻多项式系数,
Figure GDA0002286884540000062
根据公式(6)可以推出其展开公式(9)如下所示:
Figure GDA0002286884540000063
其中,
Figure GDA0002286884540000071
同理,可得公式(7)和公式(8)的展开公式。
切比雪夫误差方程的系数矩阵A如式(10)所示:
Figure GDA0002286884540000072
其中,
Figure GDA0002286884540000073
Figure GDA0002286884540000074
其中
Figure GDA0002286884540000075
τ∈[-1,1],n≥2。
步骤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)所示:
Figure GDA0002286884540000081
Y方向切比雪夫插值的第i个待加权值yqi如式(17)所示:
Figure GDA0002286884540000082
Z方向切比雪夫插值的第i个待加权值zqi如式(18)所示:
Figure GDA0002286884540000083
步骤6:将卫星轨道位置的拉格朗日插值和对应时刻的卫星轨道位置的切比雪夫插值组成待加权项目数据矩阵
如式(19)、(20)和(21)所示:
Figure GDA0002286884540000084
Figure GDA0002286884540000085
Figure GDA0002286884540000086
步骤7:计算待加权项目数据矩阵中各个方向的各待加权项的拉格朗日加权比重lli=[lxli,lyli,lzli]和切比雪夫加权比重lqi=[lxqi,lyqi,lzqi]。
本实施方式中,X方向第i个待加权值的拉格朗日加权比重lxli的计算公式如式(22)所示:
Figure GDA0002286884540000091
同理,可得Y方向第i个待加权值的拉格朗日加权比重lyli,Z方向第i个待加权值的拉格朗日加权比重lzli
X方向第i个待加权值的切比雪夫加权比重lxqi的计算公式如式(23)所示:
Figure GDA0002286884540000092
同理,可得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)
其中,
Figure GDA0002286884540000093
步骤9:根据各个方向上待加权值的熵值si=[sxi,syi,szi],计算各个方向上待加权值的熵权ei=[exi,eyi,ezi]。
本实施方式中,各个方向上待加权值的熵权exi,eyi,ezi的计算公式如下如式(27)、(28)和(29)所示:
Figure GDA0002286884540000101
Figure GDA0002286884540000102
Figure GDA0002286884540000103
步骤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不同实验的误差值平均值
Figure GDA0002286884540000104
Figure GDA0002286884540000111

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个待加权值;
步骤6:将卫星轨道位置的拉格朗日插值和对应时刻的卫星轨道位置的切比雪夫插值组成待加权项目数据矩阵
Figure FDA0002286884530000011
Figure FDA0002286884530000012
步骤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);
其中,
Figure FDA0002286884530000021
步骤9:根据各个方向上待加权值的熵值si=[sxi,syi,szi],计算各个方向上待加权值的熵权ei=[exi,eyi,ezi];
所述各个方向上待加权值的熵权exi,eyi,ezi的计算公式如下所示:
Figure FDA0002286884530000022
Figure FDA0002286884530000023
Figure FDA0002286884530000024
步骤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
2.根据权利要求1所述的基于熵权法的BDS卫星轨道位置插值方法,其特征在于,所述构造的任意插值时刻的历元基函数li(t)如下所示:
Figure FDA0002286884530000031
3.根据权利要求1所述的基于熵权法的BDS卫星轨道位置插值方法,其特征在于,所述t时刻的卫星轨道位置的拉格朗日插值[Xl(t),Yl(t),Zl(t)]的计算公式如下所示:
Figure FDA0002286884530000032
Figure FDA0002286884530000033
Figure FDA0002286884530000034
4.根据权利要求1所述的基于熵权法的BDS卫星轨道位置插值方法,其特征在于,所述构造的卫星轨道位置插值的切比雪夫误差方程V=[Vx,Vy,Vz]的计算公式如下所示:
Vx=ACx-X;
Vy=ACy-Y;
Vz=ACz-Z;
其中,C=[Cx,Cy,Cz]为等间隔历元时刻多项式系数,
Figure FDA0002286884530000035
Figure FDA0002286884530000041
为切比雪夫误差方程的系数矩阵,
Figure FDA0002286884530000042
Figure FDA0002286884530000043
其中
Figure FDA0002286884530000044
τ∈[-1,1],n≥2。
CN201710568613.9A 2017-07-13 2017-07-13 一种基于熵权法的bds卫星轨道位置插值方法 Active CN107193020B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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后处理方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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