CN110132261B - 一种基于数值拟合的高精度星上轨道预报方法 - Google Patents

一种基于数值拟合的高精度星上轨道预报方法 Download PDF

Info

Publication number
CN110132261B
CN110132261B CN201811365986.7A CN201811365986A CN110132261B CN 110132261 B CN110132261 B CN 110132261B CN 201811365986 A CN201811365986 A CN 201811365986A CN 110132261 B CN110132261 B CN 110132261B
Authority
CN
China
Prior art keywords
orbit
forecasting
satellite
instantaneous
precision
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
CN201811365986.7A
Other languages
English (en)
Other versions
CN110132261A (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.)
China Xian Satellite Control Center
Original Assignee
China Xian Satellite Control Center
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 China Xian Satellite Control Center filed Critical China Xian Satellite Control Center
Priority to CN201811365986.7A priority Critical patent/CN110132261B/zh
Publication of CN110132261A publication Critical patent/CN110132261A/zh
Application granted granted Critical
Publication of CN110132261B publication Critical patent/CN110132261B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/02Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/04Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means
    • G01C21/08Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means involving use of the magnetic field of the earth

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • General Physics & Mathematics (AREA)
  • Astronomy & Astrophysics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于数值拟合的高精度星上轨道预报方法,首先获得包括最新太阳辐射流量和地磁指数在内的大气环境参数,卫星星体参数以及卫星最新的精密轨道根数,结合星体和空间环境参数,利用地面高精度数值法进行轨道预报,获得一组预报弧段内卫星瞬时轨道根数,然后转换为平均根数进行数值拟合,得到每个轨道根数的变化项系数,根据星上预报精度要求上注多组初轨进行分段外推,根据上注的初始轨道根数及变化项系数进行轨道预报,得到预报时刻卫星的平均根数,如果预报精度要求较高,注入多组初值进行分阶段预报,将预报时刻卫星的平均根数转换为瞬时根数,从而得到卫星的瞬时位置,本发明解决了现有技术中存在的星上轨道预报精度差的问题。

Description

一种基于数值拟合的高精度星上轨道预报方法
技术领域
本发明属于航天测量与控制技术领域,具体涉及一种基于数值拟合的高精度星上轨道预报方法。
背景技术
由于受星载计算机等条件的限制,星上轨道预报通常采用简化的分析法模型,预报精度较差,对于低轨卫星,加之大气密度的影响,预报精度更差。
发明内容
本发明的目的是提供一种基于数值拟合的高精度星上轨道预报方法,解决了现有技术中存在的星上轨道预报精度差的问题。
本发明所采用的技术方案是,一种基于数值拟合的高精度星上轨道预报方法,具体按照以下步骤实施:
步骤1、获得包括最新太阳辐射流量和地磁指数在内的大气环境参数,卫星星体参数以及卫星最新的精密轨道根数,结合星体和空间环境参数,利用地面高精度数值法进行轨道预报,获得一组预报弧段内卫星瞬时轨道根数;
步骤2、将步骤1获得的预报弧段内的一组瞬时轨道根数转换为平均根数;
步骤3、用步骤2得到的预报弧段内的一组平均根数进行数值拟合,得到每个轨道根数的变化项系数;
步骤4、根据星上预报精度要求上注多组初轨进行分段外推,根据上注的初始轨道根数及变化项系数用步骤3进行轨道预报,得到预报时刻卫星的平均根数,如果预报精度要求较高,注入多组初值进行分阶段预报;
步骤5、将预报时刻卫星的平均根数转换为瞬时根数,从而得到卫星的瞬时位置。
本发明的特点还在于,
步骤1具体如下:
在地心天球坐标系中,卫星运动方程为:
Figure RE-RE-GDA0002017045140000021
其中,
Figure RE-RE-GDA0002017045140000022
为体问题下地球对卫星的引力加速度:
Figure RE-RE-GDA0002017045140000023
Figure RE-RE-GDA0002017045140000024
为除二体问题的各种摄动源共同作用下的摄动加速度:
Figure RE-RE-GDA0002017045140000025
Figure RE-RE-GDA0002017045140000026
为卫星的地心位置矢量,
Figure RE-RE-GDA0002017045140000027
为卫星的地心速度矢量,
Figure RE-RE-GDA0002017045140000028
为卫星的地心加速度矢量,μ为地心引力常数;
根据公式(1)~公式(3)进行高精度数值法轨道预报,生成一组卫星的瞬时轨道根数。
步骤2中将瞬时轨道根数转换为平均根数公式如下,卫星t时刻平均根数
Figure RE-RE-GDA0002017045140000029
为:
Figure RE-RE-GDA0002017045140000031
其中,σ(t)为卫星t时刻瞬时根数,
σ=(a,e,i,Ω,ω,M) (5)
其中,a轨道半长径,Ω升交点经度,e轨道偏心率,ω近地点幅角,i轨道倾角,M为平近点角,
Figure RE-RE-GDA0002017045140000032
Figure RE-RE-GDA0002017045140000033
为轨道变化中的一价和二价短周期项,通过以上公式将一组瞬时轨道根数转换为平均根数。
步骤3具体如下:
具体拟合方法为:设拟合数据序列包含k个数据,其中轨道半长轴分别为a1,a2,...,ak,则平均半长轴为:
Figure RE-RE-GDA0002017045140000034
轨道角速度n为:
Figure RE-RE-GDA0002017045140000035
对任意时刻的每个轨道根数,ti=[t0,tf],存在
f(ti)=ΦT(ti)P (8)
其中,i=1,2,…n,
系数向量P满足的线性方程组为:
ΦP=F (9)
Figure RE-RE-GDA0002017045140000036
Figure RE-RE-GDA0002017045140000037
P=[Pk0 Pk1 Pk1s Pk1c Pk2s Pk2c Pk3s Pk3c]T (12)
对上式中系数矩阵P进行估计求解,得到卫星每个根数的变化项系数。
步骤4中上注初轨数量为3组/天,即外推8小时换一组初轨,具体预报方法如下:
Figure RE-RE-GDA0002017045140000041
根据上式利用注入的初值和变化项系数进行轨道预报。
步骤5具体如下:
Figure RE-RE-GDA0002017045140000042
最后根据上式将预报结果转换为瞬时根数,即完成高精度星上轨道预报。
本发明的有益效果是,一种基于数值拟合的高精度星上轨道预报方法,先采用高精度数值法进行轨道预报,并将得到的一组瞬时根数转换为平均根数,基于这组平均根数进行多项式拟合,给出每个轨道根数的一阶和二阶变化项系数,再根据星上轨道预报精度要求,分几组初轨和拟合系数进行上注,最后根据预报时间判断用哪组初轨进行分段预报,从而达到提高星上轨道预报精度的目的。
附图说明
图1是本发明一种基于数值拟合的高精度星上轨道预报方法中外推24 小时的误差分析情况图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明所采用的技术方案是,一种基于数值拟合的高精度星上轨道预报方法,具体按照以下步骤实施:
步骤1、获得包括最新太阳辐射流量和地磁指数在内的大气环境参数,卫星星体参数以及卫星最新的精密轨道根数,结合星体和空间环境参数,利用地面高精度数值法进行轨道预报,获得一组预报弧段内卫星瞬时轨道根数, 具体如下:
在地心天球坐标系中,卫星运动方程为:
Figure RE-RE-GDA0002017045140000051
其中,
Figure RE-RE-GDA0002017045140000052
为体问题下地球对卫星的引力加速度:
Figure RE-RE-GDA0002017045140000053
Figure RE-RE-GDA0002017045140000054
为除二体问题的各种摄动源共同作用下的摄动加速度:
Figure RE-RE-GDA0002017045140000055
Figure RE-RE-GDA0002017045140000056
为卫星的地心位置矢量,
Figure RE-RE-GDA0002017045140000057
为卫星的地心速度矢量,
Figure RE-RE-GDA0002017045140000058
为卫星的地心加速度矢量,μ为地心引力常数;
根据公式(1)~公式(3)进行高精度数值法轨道预报,生成一组卫星的瞬时轨道根数;
步骤2、将步骤1获得的预报弧段内的一组瞬时轨道根数转换为平均根数,公式如下,卫星t时刻平均根数
Figure RE-RE-GDA0002017045140000061
为:
Figure RE-RE-GDA0002017045140000062
其中,σ(t)为卫星t时刻瞬时根数,
σ=(a,e,i,Ω,ω,M) (5)
其中,a轨道半长径,Ω升交点经度,e轨道偏心率,ω近地点幅角,i轨道倾角,M为平近点角,
Figure RE-RE-GDA0002017045140000063
Figure RE-RE-GDA0002017045140000064
为轨道变化中的一价和二价短周期项,通过以上公式将一组瞬时轨道根数转换为平均根数;
步骤3、用步骤2得到的预报弧段内的一组平均根数进行数值拟合,得到每个轨道根数的变化项系数,具体如下:
具体拟合方法为:设拟合数据序列包含k个数据,其中轨道半长轴分别为a1,a2,...,ak,则平均半长轴为:
Figure RE-RE-GDA0002017045140000065
轨道角速度n为:
Figure RE-RE-GDA0002017045140000066
对任意时刻的每个轨道根数,ti=[t0,tf],存在
f(ti)=ΦT(ti)P (8)
其中,i=1,2,…n,
系数向量P满足的线性方程组为:
ΦP=F (9)
Figure RE-RE-GDA0002017045140000067
Figure RE-RE-GDA0002017045140000071
P=[Pk0 Pk1 Pk1s Pk1c Pk2s Pk2c Pk3s Pk3c]T (12)
对上式中系数矩阵P进行估计求解,得到卫星每个根数的变化项系数;
步骤4、根据星上预报精度要求上注多组初轨进行分段外推,根据上注的初始轨道根数及变化项系数用步骤3进行轨道预报,得到预报时刻卫星的平均根数,如果预报精度要求较高,注入多组初值进行分阶段预报,其中上注初轨数量为3组/天,即外推8小时换一组初轨,具体预报方法如下:
Figure RE-RE-GDA0002017045140000072
根据上式利用注入的初值和变化项系数进行轨道预报;
步骤5、将预报时刻卫星的平均根数转换为瞬时根数,从而得到卫星的瞬时位置,具体如下:
Figure RE-RE-GDA0002017045140000073
最后根据上式将预报结果转换为瞬时根数,即完成高精度星上轨道预报。
实施例
采取如表1所示的初始轨道数据:
表1:初始轨道数据
Figure RE-RE-GDA0002017045140000081
用高精度轨道预报软件预报两天的弹道,并转换为平根数,采用二次多项式函数拟合后得到每个轨道根数的变化项的系数,得到预报公式为:
Figure RE-RE-GDA0002017045140000082
选择弹道第一点为平根数注入需求点,相应的变率输出结果如表2所示:
表2:注入需求点与变率输出结果对应表
Figure RE-RE-GDA0002017045140000083
用拟合后的系数进行外推48小时,将外推结果与标准弹道比较,其位置误差最大为355m,主要是沿迹向,在径向和法线方向的误差较小,约为 20米;若外推24小时,其最大位置误差小于100m;若外推12小时约为20m;外推6小时约为10m。如图1所示,用简化的分析法模型进行预报24小时误差约为500m,由此可见,用该方法可大大提高了星上的轨道预报精度。
本发明一种基于数值拟合的高精度星上轨道预报方法,能够解决由于星上计算机条件限制采用简化分析法轨道预报精度较差的问题,使用地面高精度数值法轨道预报,并用预报弧段内轨道平均根数进行拟合轨道变化项,并按星上轨道预报的精度要求进行多组初轨预报的方法,达到提高星上轨道预报精度的目的。

Claims (3)

1.一种基于数值拟合的高精度星上轨道预报方法,其特征在于,具体按照以下步骤实施:
步骤1、获得包括最新太阳辐射流量和地磁指数在内的大气环境参数,卫星星体参数以及卫星最新的精密轨道根数,结合星体和空间环境参数,利用地面高精度数值法进行轨道预报,获得一组预报弧段内卫星瞬时轨道根数;
所述步骤1具体如下:
在地心天球坐标系中,卫星运动方程为:
Figure FDA0003932508400000011
其中,
Figure FDA0003932508400000012
为体问题下地球对卫星的引力加速度:
Figure FDA0003932508400000013
Figure FDA0003932508400000014
为除二体问题的各种摄动源共同作用下的摄动加速度:
Figure FDA0003932508400000015
Figure FDA0003932508400000016
为卫星的地心位置矢量,
Figure FDA0003932508400000017
为卫星的地心速度矢量,
Figure FDA0003932508400000018
为卫星的地心加速度矢量,μ为地心引力常数;
根据公式(1)~公式(3)进行高精度数值法轨道预报,生成一组卫星的瞬时轨道根数;
步骤2、将步骤1获得的预报弧段内的一组瞬时轨道根数转换为平均根数;
所述步骤2中将瞬时轨道根数转换为平均根数公式如下,卫星t时刻平均根数
Figure FDA0003932508400000021
为:
Figure FDA0003932508400000022
其中,σ(t)为卫星t时刻瞬时根数,
σ=(a,e,i,Ω,ω,M) (5)
其中,a轨道半长径,Ω升交点经度,e轨道偏心率,ω近地点幅角,i轨道倾角,M为平近点角,
Figure FDA0003932508400000023
Figure FDA0003932508400000024
为轨道变化中的一价和二价短周期项,通过以上公式将一组瞬时轨道根数转换为平均根数;
步骤3、用步骤2得到的预报弧段内的一组平均根数进行数值拟合,得到每个轨道根数的变化项系数;
所述步骤3具体如下:
具体拟合方法为:设拟合数据序列包含k个数据,其中轨道半长轴分别为a1,a2,...,ak,则平均半长轴为:
Figure FDA0003932508400000025
轨道角速度n为:
Figure FDA0003932508400000026
对任意时刻的每个轨道根数,ti=[t0,tf],存在
f(ti)=ΦT(ti)P (8)
其中,i=1,2,…n,
系数向量P满足的线性方程组为:
ΦP=F (9)
Figure FDA0003932508400000031
Figure FDA0003932508400000032
P=[Pk0 Pk1 Pk1s Pk1c Pk2s Pk2c Pk3s Pk3c]T (12)
对上式中系数矩阵P进行估计求解,得到卫星每个根数的变化项系数;
步骤4、根据星上预报精度要求上注多组初轨进行分段外推,根据上注的初始轨道根数及变化项系数用步骤3进行轨道预报,得到预报时刻卫星的平均根数,如果预报精度要求较高,注入多组初值进行分阶段预报;
步骤5、将预报时刻卫星的平均根数转换为瞬时根数,从而得到卫星的瞬时位置。
2.根据权利要求1所述的一种基于数值拟合的高精度星上轨道预报方法,其特征在于,所述步骤4中上注初轨数量为3组/天,即外推8小时换一组初轨,具体预报方法如下:
Figure FDA0003932508400000033
根据上式利用注入的初值和变化项系数进行轨道预报。
3.根据权利要求1所述的一种基于数值拟合的高精度星上轨道预报方法,其特征在于,所述步骤5具体如下:
Figure FDA0003932508400000041
最后根据上式将预报结果转换为瞬时根数,即完成高精度星上轨道预报。
CN201811365986.7A 2018-11-16 2018-11-16 一种基于数值拟合的高精度星上轨道预报方法 Active CN110132261B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811365986.7A CN110132261B (zh) 2018-11-16 2018-11-16 一种基于数值拟合的高精度星上轨道预报方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811365986.7A CN110132261B (zh) 2018-11-16 2018-11-16 一种基于数值拟合的高精度星上轨道预报方法

Publications (2)

Publication Number Publication Date
CN110132261A CN110132261A (zh) 2019-08-16
CN110132261B true CN110132261B (zh) 2023-03-14

Family

ID=67568181

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811365986.7A Active CN110132261B (zh) 2018-11-16 2018-11-16 一种基于数值拟合的高精度星上轨道预报方法

Country Status (1)

Country Link
CN (1) CN110132261B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110816896B (zh) * 2019-10-28 2021-06-11 中国空间技术研究院 一种卫星星上简易轨道外推方法
CN111060934B (zh) * 2019-12-25 2023-08-04 中国西安卫星测控中心 低轨航天器短弧定轨瞬时根数精度评估方法
CN111854765B (zh) * 2020-06-08 2022-04-26 中国人民解放军战略支援部队航天工程大学 一种中轨道导航卫星轨道长期预报方法
CN111953401B (zh) * 2020-07-28 2022-06-07 中国西安卫星测控中心 一种微小卫星自主请求式轨道服务系统
CN112949067B (zh) * 2021-03-05 2022-11-25 航天科工空间工程发展有限公司 一种星载空间目标轨道平滑方法
CN113297745B (zh) * 2021-05-28 2022-09-02 中国人民解放军63921部队 一种基于短弧拟合位置的双弧段轨道改进方法
CN113447957B (zh) * 2021-06-23 2022-11-25 北京航空航天大学 一种高精度实时卫星轨道瞬时半长轴的去振方法
CN113641949B (zh) * 2021-08-05 2023-03-28 中国西安卫星测控中心 一种地球同步转移段轨道根数高精度拟合方法
CN114383619B (zh) * 2021-12-07 2023-09-05 上海航天控制技术研究所 一种高精度轨道计算方法
CN114282166B (zh) * 2021-12-22 2024-04-19 中国西安卫星测控中心 基于环境函数矩阵的卫星经度计算方法、变轨判别方法
CN115402539B (zh) * 2022-11-02 2023-01-31 宁波天巡科技有限公司 一种基于轨道和空间环境数据的卫星变轨检测方法及系统
CN115828035B (zh) * 2023-02-16 2023-04-21 北京控制工程研究所 静止轨道卫星太阳光压干扰力矩的在轨拟合估计方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002249100A (ja) * 2001-02-21 2002-09-03 Mitsubishi Electric Corp 衛星追尾システム及び衛星追尾方法
CN102591343A (zh) * 2012-02-09 2012-07-18 航天东方红卫星有限公司 基于两行根数的卫星轨道维持控制方法
WO2012125293A2 (en) * 2011-03-11 2012-09-20 Sorce4 Llc. Offline ephemeris prediction
CN106959457A (zh) * 2017-04-25 2017-07-18 清华大学 一种用于卫星导航的glonass历书参数估计方法
CN107797130A (zh) * 2017-10-16 2018-03-13 中国西安卫星测控中心 低轨航天器多点多参数轨道上行数据计算方法
CN108761507A (zh) * 2018-05-21 2018-11-06 中国人民解放军战略支援部队信息工程大学 基于短弧定轨和预报的导航卫星轨道快速恢复方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002249100A (ja) * 2001-02-21 2002-09-03 Mitsubishi Electric Corp 衛星追尾システム及び衛星追尾方法
WO2012125293A2 (en) * 2011-03-11 2012-09-20 Sorce4 Llc. Offline ephemeris prediction
CN102591343A (zh) * 2012-02-09 2012-07-18 航天东方红卫星有限公司 基于两行根数的卫星轨道维持控制方法
CN106959457A (zh) * 2017-04-25 2017-07-18 清华大学 一种用于卫星导航的glonass历书参数估计方法
CN107797130A (zh) * 2017-10-16 2018-03-13 中国西安卫星测控中心 低轨航天器多点多参数轨道上行数据计算方法
CN108761507A (zh) * 2018-05-21 2018-11-06 中国人民解放军战略支援部队信息工程大学 基于短弧定轨和预报的导航卫星轨道快速恢复方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于GPS测量数据的卫星在轨轨道预报算法研究;刘燎等;《上海航天》;20170430(第02期);全文 *
基于动力学轨道拟合的LEO卫星轨道预报精度分析;王友存等;《全球定位系统》;20180831(第04期);全文 *

Also Published As

Publication number Publication date
CN110132261A (zh) 2019-08-16

Similar Documents

Publication Publication Date Title
CN110132261B (zh) 一种基于数值拟合的高精度星上轨道预报方法
Müller et al. High performance clocks and gravity field determination
CN107797130B (zh) 低轨航天器多点多参数轨道上行数据计算方法
CN107329146B (zh) 一种导航卫星低轨监测星座的优化设计方法
CN112257343B (zh) 一种高精度地面轨迹重复轨道优化方法及系统
CN102591343B (zh) 基于两行根数的卫星轨道维持控制方法
Canuto et al. Spacecraft dynamics and control: the embedded model control approach
CN100501331C (zh) 基于x射线脉冲星的导航卫星自主导航系统与方法
CN102495950B (zh) 一种适用于太阳同步轨道的倾角偏置量获取方法
CN110068845B (zh) 一种基于平根数理论确定卫星理论轨道的方法
CN110609972B (zh) 一种指定发射仰角的自由弹道构造方法
CN109323698B (zh) 一种空间目标陨落多模型跟踪引导方法
CN110816896B (zh) 一种卫星星上简易轨道外推方法
Verma et al. Electron density distribution and solar plasma correction of radio signals using MGS, MEX, and VEX spacecraft navigation data and its application to planetary ephemerides
CN102878997A (zh) 一种大偏心率轨道的星上快速高精度外推方法
CN104751012A (zh) 沿飞行弹道的扰动引力快速逼近方法
CN111238489B (zh) 一种低轨卫星大气阻力摄动建模和计算方法
CN104501804A (zh) 一种基于gps测量数据的卫星在轨轨道预报方法
CN112769466B (zh) 低轨卫星星座构型保持方法
CN112240941B (zh) 一种重力卫星星载加速度计的相对标校方法及系统
CN108875174A (zh) 一种基于多段打靶法的不变拟周期轨道确定方法
CN110059285B (zh) 考虑j2项影响的导弹自由段弹道偏差解析预报方法
Rutkowska et al. SLR technique used for description of the Earth elasticity
CN114771873B (zh) 一种超低轨卫星轨道自主精确维持方法
CN103196451A (zh) 一种基于Fisher信息矩阵的导航脉冲星选择方法

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