CN107193020A - 一种基于熵权法的bds卫星轨道位置插值方法 - Google Patents
一种基于熵权法的bds卫星轨道位置插值方法 Download PDFInfo
- Publication number
- CN107193020A CN107193020A CN201710568613.9A CN201710568613A CN107193020A CN 107193020 A CN107193020 A CN 107193020A CN 201710568613 A CN201710568613 A CN 201710568613A CN 107193020 A CN107193020 A CN 107193020A
- Authority
- CN
- China
- Prior art keywords
- mrow
- interpolation
- msub
- chebyshev
- orbital position
- 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
- 238000000034 method Methods 0.000 title claims abstract description 17
- 239000011159 matrix material Substances 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000010276 construction Methods 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 6
- 238000013459 approach Methods 0.000 abstract description 2
- 238000002474 experimental method Methods 0.000 description 9
- DMBHHRLKUKUOEG-UHFFFAOYSA-N diphenylamine Chemical compound C=1C=CC=CC=1NC1=CC=CC=C1 DMBHHRLKUKUOEG-UHFFFAOYSA-N 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
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
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和ti，根据已知历元时刻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：将卫星轨道位置的拉格朗日插值和对应时刻的卫星轨道位置的切比雪夫插值组成待加权项目数据矩阵
步骤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(lxli ln lxli+lxqi ln lxqi).
syi＝-k(lyli ln lyli+lyqi ln lyqi).
szi＝-k(lzli ln lzli+lzqi ln lzqi).
其中，
步骤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＝eiXqi+(1-ei)Xli；
wyi＝exiYqi+(1-exi)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；
其中，C＝[Cx，Cy，Cz]为等间隔历元时刻多项式系数，
为切比雪夫误差方程的系数矩阵， |τ≤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)所示：
其中，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)
其中，C＝[Cx，Cy，Cz]为等间隔历元时刻多项式系数，
根据公式(6)可以推出其展开公式(9)如下所示：
其中，
同理，可得公式(7)和公式(8)的展开公式。
切比雪夫误差方程的系数矩阵A如式(10)所示：
其中， |τ≤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)
Xq(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(lxli ln lxli+lxqi ln lxqi) (24)
syi＝-k(lyli ln lyli+lyqi ln lyqi) (25)
szi＝-k(lzli ln lzli+lzqi ln lzqi) (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＝eiXqi+(1-ei)Xli (30)
wyi＝exiYqi+(1-exi)Yli (31)
wzi＝eziZqi+(1-ezi)Zli (32)
本发明实施方式中，采用了六组实验得到的不同历元的插值误差值，如图2所示。不同实验的误差值平均值如表1所示：
表1不同实验的误差值平均值
Claims (4)
1.一种基于熵权法的BDS卫星轨道位置插值方法，其特征在于，包括以下步骤：
步骤1：获取BDS卫星轨道系统精密星历中的等间隔的已知历元时刻ti和ti，根据已知历元时刻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：将卫星轨道位置的拉格朗日插值和对应时刻的卫星轨道位置的切比雪夫插值组成待加权项目数据矩阵
步骤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(lxli ln lxli+lxqi ln lxqi)；
syi＝-k(lyli ln lyli+lyqi ln lyqi)；
szi＝-k(lzli ln lzli+lzqi ln lzqi)；
其中，
步骤9：根据各个方向上待加权值的熵值si＝[sxi，syi，szi]，计算各个方向上待加权值的熵权ei＝[exi，eyi，ezi]；
所述各个方向上待加权值的熵权exi，eyi，ezi的计算公式如下所示：
<mrow>
<msub>
<mi>e</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>|</mo>
<mfrac>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>|</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>e</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>|</mo>
<mfrac>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>|</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>e</mi>
<mrow>
<mi>z</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>|</mo>
<mfrac>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>z</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>z</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>|</mo>
</mrow>
<mo>;</mo>
</mrow>
步骤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＝eiXqi+(1-ei)Xli；
wyi＝exiYqi+(1-exi)Yli；
wzi＝eziZqi+(1-ezi)Zli。
2.根据权利要求1所述的基于熵权法的BDS卫星轨道位置插值方法，其特征在于，所述构造的任意插值时刻的历元基函数li(t)如下所示：
<mrow>
<msub>
<mi>l</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<munder>
<mi>&Pi;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>t</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
3.根据权利要求1所述的基于熵权法的BDS卫星轨道位置插值方法，其特征在于，所述t时刻的卫星轨道位置的拉格朗日插值[Xl(t)，Yl(t)，Zl(t)]的计算公式如下所示：
<mrow>
<mi>X</mi>
<mi>l</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<msub>
<mi>l</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<mi>Y</mi>
<mi>l</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>Y</mi>
<mi>i</mi>
</msub>
<msub>
<mi>l</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<mi>Z</mi>
<mi>l</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>Z</mi>
<mi>i</mi>
</msub>
<msub>
<mi>l</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
4.根据权利要求1所述的基于熵权法的BDS卫星轨道位置插值方法，其特征在于，所述构造的卫星轨道位置插值的切比雪夫误差方程V＝[Vx，Vy，Vz]的计算公式如下所示：
Vx＝ACx-X；
Vy＝ACy-Y；
Vz＝ACz-Z；
其中，C＝[Cx，Cy，Cz]为等间隔历元时刻多项式系数， 为切比雪夫误差方程的系数矩阵， |τ≤1|，n≥2。
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 true CN107193020A (zh) | 2017-09-22 |
| CN107193020B 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) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111209523A (zh) * | 2020-01-06 | 2020-05-29 | 中国科学院紫金山天文台 | 适用于大偏心率轨道密集星历精密计算的快速处理方法 |
| CN112346086A (zh) * | 2020-10-23 | 2021-02-09 | 西安空间无线电技术研究所 | 一种基于临近空间浮空平台的高效快速掩星预报方法 |
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 |
|---|
| LIU WEIPING,ET AL: "A new interpolation method based on satellite physical character in using IGS precise ephemeris", 《GEODESY AND GEODYNAMICS》 * |
| 汪威等: "北斗三类卫星精密星历内插方法分析比较", 《全球定位系统》 * |
| 王晓明等: "基于阶次组合的 GPS精密星历插值研究", 《大地测量与地球动力学》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111209523A (zh) * | 2020-01-06 | 2020-05-29 | 中国科学院紫金山天文台 | 适用于大偏心率轨道密集星历精密计算的快速处理方法 |
| CN111209523B (zh) * | 2020-01-06 | 2020-12-29 | 中国科学院紫金山天文台 | 适用于大偏心率轨道密集星历精密计算的快速处理方法 |
| US11319094B2 (en) * | 2020-01-06 | 2022-05-03 | Purple Mountain Observatory, Chinese Academy Of Sciences | Method for accurately and efficiently calculating dense ephemeris of high-eccentricity orbit |
| CN112346086A (zh) * | 2020-10-23 | 2021-02-09 | 西安空间无线电技术研究所 | 一种基于临近空间浮空平台的高效快速掩星预报方法 |
| CN112346086B (zh) * | 2020-10-23 | 2024-03-15 | 西安空间无线电技术研究所 | 一种基于临近空间浮空平台的高效快速掩星预报方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107193020B (zh) | 2020-04-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105629263B (zh) | 一种对流层大气延迟误差估计改正方法和改正系统 | |
| CN103424116B (zh) | 一种适应轨道机动的地球同步卫星精密定轨方法 | |
| CN106250616B (zh) | 基于动态可靠性估计的导航卫星星座可用性评估方法 | |
| CN103760572A (zh) | 一种基于区域cors的单频ppp电离层加权方法 | |
| CN106468774A (zh) | 一种应用于星基增强系统的星历星钟改正参数及空间信号完好性参数方法 | |
| CN104865587B (zh) | 基于geo星座粗定位的北斗接收机快速定位方法 | |
| CN106168672A (zh) | 一种gnss多模单频rtk周跳探测方法及装置 | |
| CN103777518B (zh) | 基于改进的最小二乘法的卫星导航授时系统及方法 | |
| CN106125110A (zh) | 基于分区改正的星基增强系统定位精度提高方法 | |
| CN104483689B (zh) | 一种bds参考站三频载波相位整周模糊度确定方法 | |
| CN104090280B (zh) | 一种基于区域cors的电离层延迟修正预报方法 | |
| CN102591343A (zh) | 基于两行根数的卫星轨道维持控制方法 | |
| CN103499822A (zh) | 一种基于最优gdop和牛顿恒等式的双星座快速选星方法 | |
| CN105445765B (zh) | 一种北斗导航卫星系统混合星座的伪距测量方法 | |
| CN107193020A (zh) | 一种基于熵权法的bds卫星轨道位置插值方法 | |
| CN107656294B (zh) | 一种基于选星模板的多卫星导航系统选星方法 | |
| CN106093967A (zh) | 一种伪距相位综合的电离层延迟求解方法 | |
| CN101226062B (zh) | 一种星上实时计算环月轨道的方法 | |
| CN105372685A (zh) | 一种区域高精度位置增强系统及方法 | |
| CN106292265A (zh) | 一种基于导航卫星的多地时间同步方法 | |
| CN103698790A (zh) | 北斗与gps双系统宽巷载波相位混频星间差分组合方法 | |
| CN107071893A (zh) | 一种蜂窝网络rtk定位方法及系统 | |
| CN108415039B (zh) | iGMAS多分析中心多卫星系统精密轨道产品综合方法 | |
| EP2817660B1 (en) | Low power asynchronous gps baseband processor | |
| CN105629278A (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 | ||
| 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 |