CN114383619B - 一种高精度轨道计算方法 - Google Patents
一种高精度轨道计算方法 Download PDFInfo
- Publication number
- CN114383619B CN114383619B CN202111485208.3A CN202111485208A CN114383619B CN 114383619 B CN114383619 B CN 114383619B CN 202111485208 A CN202111485208 A CN 202111485208A CN 114383619 B CN114383619 B CN 114383619B
- Authority
- CN
- China
- Prior art keywords
- orbit
- satellite
- precision
- vector matrix
- coefficient vector
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
- G01C21/16—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
- G01C21/165—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Remote Sensing (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Automation & Control Theory (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Astronomy & Astrophysics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种高精度轨道计算方法,包括:根据不同的轨道参数受摄动力的影响不同,以及星上轨道预报精度需求,选取不同的基函数,建立基于所述基函数的星上轨道预报模型;根据地面测控系统进行高精度轨道预报,获取所述基函数的系数向量矩阵的估计值;根据所述系数向量矩阵和所述星上轨道预报模型,通过递推得到惯性位置和惯性速度;上注系数向量矩阵和初始时刻的系统钟,本发明提高了卫星星上预报的精度,大幅度简化了卫星星上基于轨道动力学的轨道预报模型的复杂程度,同时避免了轨道参数地面频繁上注和上注数据过多的问题,降低了星载计算机和地面测控的压力,具有较强的工程可实现性。
Description
技术领域
本发明涉及卫星星上轨道递推算法技术领域,具体涉及一种高精度轨道计算方法。
背景技术
地球静止轨道卫星、大椭圆轨道等高轨道卫星,从入轨阶段,到进入标称轨道后,通常都由给出卫星初始的位置和速度,然后由星上高精度轨道动力学模型进行递推,实现星上轨道预报。
但是上述方法在卫星轨道精度要求越来越高时,卫星星上轨道预报模型变得越来越复杂;同时也要求地面测控对上注的卫星初始位置和速度的频度越来越快,增加了地面测控的任务量。因而迫切需要寻求一种高精度的轨道预报算法,在保证预报精度的前提下,实现预报模型简化同时减小地面上注频度。
发明内容
本发明的目的是为了提供一种高精度轨道计算方法。此方法旨在解决传统方法中卫星轨道精度要求越来越高时,卫星星上轨道预报模型变得越来越复杂、要求地面测控对上注的卫星初始位置和速度的频度越来越快,增加了地面测控的任务量的问题。
为达到上述目的,本发明提供了一种高精度轨道计算方法,包括:
步骤S1:根据不同的轨道参数受摄动力的影响不同,以及星上轨道预报精度需求,选取不同的基函数,建立基于所述基函数的星上轨道预报模型;
步骤S2:根据地面测控系统进行高精度轨道预报,获取所述基函数的系数向量矩阵的估计值;
步骤S3:根据所述系数向量矩阵和所述星上轨道预报模型,通过递推得到惯性位置和惯性速度;
步骤S4:上注系数向量矩阵和初始时刻的系统钟,根据所述星上轨道预报模型,以实现在轨高精度计算卫星轨道。
优选的,在所述步骤S1中,所述摄动力包括:地球非球形摄动、日月三体引力摄动和太阳光压摄动。
优选的,所述步骤S2具体包括:
步骤S2.1:由所述地面测控系统利用地面资源对卫星轨道进行高精度预报,得到在一定时间内卫星的轨道六根数;
步骤S2.2:假定一采样频率,并取若干个星历数据,根据所述卫星的轨道六根数,计算所述系数向量矩阵的估计值。
优选的,所述初始时刻为第一个星历数据的系统钟时刻。
优选的,将所述步骤S3中得到的惯性位置和惯性速度与精确星历表进行对比,若三轴位置和速度误差小于或等于设定的阈值,则系数向量矩阵计算正确;否则,则缩短拟合总时长。
优选的,所述步骤S1中,对于不同的轨道,选取不同的基函数,以地球静止轨道为例,建立的所述基函数的星上轨道预报模型的递推表达式为:
其中,a为半长轴,λ为平经度,ex为卫星轨道在x轴的偏心率,ey为卫星轨道在y轴的偏心率,ix为卫星轨道与x轴的倾角,iy为卫星轨道与y轴的倾角;
下角标0表示常数项系数;
下角标ic表示第i阶的cos函数系数;
下角标is表示第i阶的sin函数系数;
t表示相对参考时刻的时间;
l表示平赤经;
Ck表示参数k的系数,k为参数a,λ,ex,ey,ix和iy中的任意一个,其中,k分别表示为a,λ,ex,ey,ix和iy。
优选的,所述λ和l的表达式分别为:
λ=ω+Ω+M-θ
l=ω+Ω+M
其中,ω为轨道近地点辐角;
Ω为轨道升交点赤经;
M为轨道平近点角;
θ为格林尼治时角。
优选的,所述步骤S2具体包括:
步骤S2.1:根据地面测控系统进行高精度轨道预报,并根据地面动力学,获取卫星的精密星历参数;
步骤S2.2:根据所述卫星的精密星历参数,获取所述基函数的系数向量矩阵的估计值。
优选的,所述卫星的精密星历参数为fk(tj),则对星历表中的任意点,tj∈[t0,tf],存在
f(tj)=ΦT(tj)P,j=1,2,…,n
其中:P为系数向量矩阵,Φ(tj)为基函数;
tj为第j个采样时间点;
t0为精密星历参数对应的初始时间点;
tf为精密星历参数对应的结束时间点,
所述系数向量矩阵P满足的线性方程组为:
根据超定方程组的最小二乘法对系数矩阵P进行估计,得到:
其中:Φ(t)=(Φ0(t)…Φn(t))T。
优选的,所述基函数基于卫星空间摄动力对轨道根数的影响规律和所述星上轨道预报精度需求,选取不同的阶数。
与现有技术相比,本发明具有以下有益效果:
本发明充分利用地面高精度轨道预报的资源,结合卫星轨道参数自身在摄动影响下的变化规律,进行星上轨道预报,提高了卫星星上预报的精度,大幅度简化了卫星星上基于轨道动力学的轨道预报模型的复杂程度,同时避免轨道参数地面频繁上注和上注数据过多的问题,从而降低了星载计算机和地面测控的压力,具有较强的工程可实现性。
附图说明
为了更清楚地说明本发明的技术方案,下面将对描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一个实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图:
图1为本发明一实施例提供的一种高精度轨道计算方法的流程示意图。
具体实施方式
以下结合附图1和具体实施方式对本发明提出的一种高精度轨道计算方法作进一步详细说明。根据下面说明,本发明的优点和特征将更清楚。需要说明的是,附图采用非常简化的形式且均使用非精准的比例,仅用以方便、明晰地辅助说明本发明实施方式的目的。为了使本发明的目的、特征和优点能够更加明显易懂,请参阅附图。须知,本说明书所附图式所绘示的结构、比例、大小等,均仅用以配合说明书所揭示的内容,以供熟悉此技术的人士了解与阅读,并非用以限定本发明实施的限定条件,故不具技术上的实质意义,任何结构的修饰、比例关系的改变或大小的调整,在不影响本发明所能产生的功效及所能达成的目的下,均应仍落在本发明所揭示的技术内容能涵盖的范围内。
鉴于现有技术中卫星轨道精度要求越来越高时,卫星星上轨道预报模型变得越来越复杂、要求地面测控对上注的卫星初始位置和速度的频度越来越快,增加了地面测控的任务量的问题,为了实现对高轨道卫星在转移段以及入轨后对卫星轨道进行高精度预报,简化卫星星上轨道递推模型,降低地面上注频度,本实施例提供了一种高精度轨道计算方法,包括以下步骤:
步骤S1:建立基于基函数的星上高精度轨道预报模型:
针对不同的轨道,选取不同的轨道根数表示方式。根据卫星所处的空间环境,对于高轨卫星考虑地球非球形摄动、日月三体引力摄动、太阳光压摄动三个主要摄动力的影响,建立卫星轨道根数受摄动力的影响规律模型。(以地球静止轨道为例说明星上高精度轨道预报模型的建立过程)
步骤S2:针对地球静止轨道卫星,为了避免由于小偏心率和小倾角引起的歧义,选取以下的轨道根数表述方式:
选取半长轴a,平经度λ和平赤经l,则
λ=ω+Ω+M-θ (1)
l=ω+Ω+M (2)
其中,ω为轨道近地点辐角;
Ω为轨道升交点赤经;
M为轨道平近点角;
θ为格林尼治时角。
摄动力引起的平经度的变化率Δλ与平赤经变化量Δl一致,即Δλ=Δl。
卫星轨道偏心率矢量e为:
卫星轨道倾角矢量i为:
建立如下式(6)形式的星上轨道递推模型:
其中:下角标0表示常数项系数;
下角标ic表示第i阶的cos函数系数;
下角标is表示第i阶的sin函数系数;
t表示相对参考时刻的时间;
Ck表示参数k的系数,k为参数a,λ,ex,ey,ix和iy中的任意一个,其中,k分别表示为a,λ,ex,ey,ix和iy。
结合地面测控系统进行高精度轨道预报,由地面测控系统经过高精度定轨,并由地面动力学获取卫星的精密星历参数fk(ti)。
利用地面高精度轨道预报,获取基函数系数矩阵P。
以半长轴a为例,根据星上递推模型采用的基函数为:Φ(t)=(1 t cos nt sinnt cos 2nt sin 2nt)T,系数向量矩阵对星历表中的数据采样点,tj∈[t0,tf],存在
f(tj)=ΦT(tj)P,j=1,2,…,n (7)
tj为第j个采样时间点;
t0为精密星历参数对应的初始时间点;
tf为精密星历参数对应的结束时间点。
因此有,系数向量矩阵P满足的线性方程组为:
利用超定方程组的最小二乘法对系数向量矩阵P进行估计,得到:
其中:Φ(t)=(Φ0(t)…Φn(t))T为基函数,该基函数基于卫星空间摄动对轨道根数的影响规律,并根据卫星星上轨道预报精度要求,选取合适的阶数。
重复上述步骤S2,获取其他轨道根数的基函数系数向量矩阵的估计值。
步骤S3:根据计算的系数向量矩阵估计值和星上轨道预报模型,递推得到惯性位置、惯性速度,并与精确星历表进行比对,若三轴位置和速度误差小于设定的阈值,则系数向量矩阵计算正确,可以使用;若超出阈值,则需缩短拟合总时长。
步骤S4:上注系数向量矩阵和初始时刻的系统钟,利用星上轨道预报模型,在轨可实现高精度计算卫星轨道。
综上所述,本实施例充分利用地面高精度轨道预报的资源,结合卫星轨道参数自身在摄动影响下的变化规律,进行星上轨道预报,提高了卫星星上预报的精度,大幅度简化了卫星星上基于轨道动力学的轨道预报模型的复杂程度,同时避免轨道参数地面频繁上注和上注数据过多的问题,从而降低了星载计算机和地面测控的压力,具有较强的工程可实现性。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
应当注意的是,在本文的实施方式中所揭露的装置和方法,也可以通过其他的方式实现。以上所描述的装置实施方式仅仅是示意性的,例如,附图中的流程图和框图显示了根据本文的多个实施方式的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现方式中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用于执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (9)
1.一种高精度轨道计算方法,其特征在于,包括:
步骤S1:根据不同的轨道参数受摄动力的影响不同,以及星上轨道预报精度需求,选取不同的基函数,建立基于所述基函数的星上轨道预报模型,
对于不同的轨道,选取不同的基函数,以地球静止轨道为例,建立的所述基函数的星上轨道预报模型的递推表达式为:
其中,a为半长轴,λ为平经度,ex为卫星轨道在x轴的偏心率,ey为卫星轨道在y轴的偏心率,ix为卫星轨道与x轴的倾角,iy为卫星轨道与y轴的倾角;
下角标0表示常数项系数;
下角标ic表示第i阶的cos函数系数;
下角标is表示第i阶的sin函数系数;
t表示相对参考时刻的时间;
l表示平赤经;
Ck表示参数k的系数,k为参数a,λ,ex,ey,ix和iy中的任意一个参数,其中,k分别表示为a,λ,ex,ey,ix和iy;
步骤S2:根据地面测控系统进行高精度轨道预报,获取所述基函数的系数向量矩阵的估计值;
步骤S3:根据所述系数向量矩阵的估计值和所述星上轨道预报模型,通过递推得到惯性位置和惯性速度;
步骤S4:上注系数向量矩阵和初始时刻的系统钟,根据所述星上轨道预报模型,以实现在轨高精度计算卫星轨道。
2.如权利要求1所述的高精度轨道计算方法,其特征在于,在所述步骤S1中,所述摄动力包括:地球非球形摄动、日月三体引力摄动和太阳光压摄动。
3.如权利要求2所述的高精度轨道计算方法,其特征在于,所述步骤S2具体包括:
步骤S2.1:由所述地面测控系统利用地面资源对卫星轨道进行高精度预报,得到在一定时间内卫星的轨道六根数;
步骤S2.2:假定一采样频率,并取若干个星历数据,根据所述卫星的轨道六根数,计算所述系数向量矩阵的估计值。
4.如权利要求3所述的高精度轨道计算方法,其特征在于,所述初始时刻为第一个星历数据的系统钟时刻。
5.如权利要求4所述的高精度轨道计算方法,其特征在于,将所述步骤S3中得到的惯性位置和惯性速度与精确星历表进行对比,若惯性位置和惯性速度误差小于或等于设定的阈值,则所述系数向量矩阵的估计值计算正确;否则,则缩短拟合总时长。
6.如权利要求5所述的高精度轨道计算方法,其特征在于,所述平经度λ与平赤经l的表达式分别为:
λ=ω+Ω+M-θ
l=ω+Ω+M
其中,ω为轨道近地点辐角;
Ω为轨道升交点赤经;
M为轨道平近点角;
θ为格林尼治时角。
7.如权利要求6所述的高精度轨道计算方法,其特征在于,所述步骤S2具体包括:
步骤S2.1:根据地面测控系统进行高精度轨道预报,并根据地面动力学,获取卫星的精密星历参数;
步骤S2.2:根据所述卫星的精密星历参数,获取所述基函数的系数向量矩阵的估计值。
8.如权利要求7所述的高精度轨道计算方法,其特征在于,所述卫星的精密星历参数为fk(tj),则对星历表中的任意点,tj∈[t0,tf],存在
f(tj)=ΦT(tj)P,j=1,2,…,n
其中:P为系数向量矩阵,Φ(tj)为基函数;
tj为第j个采样时间点;
t0为精密星历参数对应的初始时间点;
tf为精密星历参数对应的结束时间点,
所述系数向量矩阵P满足的线性方程组为:
根据超定方程组的最小二乘法对系数矩阵P进行估计,得到:
其中:Φ(t)=(Φ0(t)…Φn(t))T。
9.权利要求8所述的高精度轨道计算方法,其特征在于,所述基函数基于卫星空间摄动力对轨道根数的影响规律和所述星上轨道预报精度需求,选取不同的阶数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111485208.3A CN114383619B (zh) | 2021-12-07 | 2021-12-07 | 一种高精度轨道计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111485208.3A CN114383619B (zh) | 2021-12-07 | 2021-12-07 | 一种高精度轨道计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114383619A CN114383619A (zh) | 2022-04-22 |
CN114383619B true CN114383619B (zh) | 2023-09-05 |
Family
ID=81196425
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111485208.3A Active CN114383619B (zh) | 2021-12-07 | 2021-12-07 | 一种高精度轨道计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114383619B (zh) |
Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102878997A (zh) * | 2012-10-24 | 2013-01-16 | 北京控制工程研究所 | 一种大偏心率轨道的星上快速高精度外推方法 |
CN103017760A (zh) * | 2011-09-27 | 2013-04-03 | 上海航天控制工程研究所 | 一种大椭圆轨道火星探测器自主对火定向方法 |
US8538682B1 (en) * | 2009-04-24 | 2013-09-17 | Qualcomm Incorporated | Systems and methods for satellite navigation using locally generated ephemeris data |
CN104048664A (zh) * | 2014-07-01 | 2014-09-17 | 南京航空航天大学 | 一种导航卫星星座自主定轨的方法 |
CN104298647A (zh) * | 2014-09-30 | 2015-01-21 | 北京航空航天大学 | 基于低轨道地球卫星的地影时刻预报的星上确定方法 |
CN104501804A (zh) * | 2014-12-17 | 2015-04-08 | 深圳航天东方红海特卫星有限公司 | 一种基于gps测量数据的卫星在轨轨道预报方法 |
CN106124170A (zh) * | 2016-08-26 | 2016-11-16 | 上海航天控制技术研究所 | 一种基于高精度姿态信息的相机光轴指向计算方法 |
CN107451656A (zh) * | 2017-08-08 | 2017-12-08 | 上海航天控制技术研究所 | 一种火星探测器器上自主轨道计算方法 |
CN109059937A (zh) * | 2018-08-21 | 2018-12-21 | 深圳市天智运控科技有限公司 | 一种星上自主卫星轨道预报方法 |
CN109774974A (zh) * | 2019-01-31 | 2019-05-21 | 上海微小卫星工程中心 | 一种用于空间碎片接近的轨道设计方法 |
CN109990779A (zh) * | 2019-04-30 | 2019-07-09 | 桂林电子科技大学 | 一种惯性导航系统及方法 |
CN110132261A (zh) * | 2018-11-16 | 2019-08-16 | 中国西安卫星测控中心 | 一种基于数值拟合的高精度星上轨道预报方法 |
CN110595485A (zh) * | 2019-08-08 | 2019-12-20 | 上海航天电子有限公司 | 基于两行根数的低轨卫星长期轨道预报方法 |
CN111102981A (zh) * | 2019-12-09 | 2020-05-05 | 北京空间机电研究所 | 一种基于ukf的高精度卫星相对导航方法 |
CN111947667A (zh) * | 2020-06-24 | 2020-11-17 | 火眼位置数智科技服务有限公司 | 一种基于运动学和动力学组合的低轨卫星实时高精度定轨方法 |
CN112257343A (zh) * | 2020-10-22 | 2021-01-22 | 上海卫星工程研究所 | 一种高精度地面轨迹重复轨道优化方法及系统 |
CN112713922A (zh) * | 2019-10-25 | 2021-04-27 | 北京航空航天大学 | 一种多波束通讯卫星的可见性快速预报算法 |
CN113447043A (zh) * | 2021-05-21 | 2021-09-28 | 北京控制工程研究所 | 一种基于gnss的卫星天文导航系统误差自主标定方法及系统 |
CN113624243A (zh) * | 2021-08-16 | 2021-11-09 | 中国科学院上海天文台 | 一种近地轨道卫星的星上实时轨道预报方法 |
CN113641949A (zh) * | 2021-08-05 | 2021-11-12 | 中国西安卫星测控中心 | 一种地球同步转移段轨道根数高精度拟合方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9075140B2 (en) * | 2009-09-23 | 2015-07-07 | Purdue Research Foundation | GNSS ephemeris with graceful degradation and measurement fusion |
-
2021
- 2021-12-07 CN CN202111485208.3A patent/CN114383619B/zh active Active
Patent Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8538682B1 (en) * | 2009-04-24 | 2013-09-17 | Qualcomm Incorporated | Systems and methods for satellite navigation using locally generated ephemeris data |
CN103017760A (zh) * | 2011-09-27 | 2013-04-03 | 上海航天控制工程研究所 | 一种大椭圆轨道火星探测器自主对火定向方法 |
CN102878997A (zh) * | 2012-10-24 | 2013-01-16 | 北京控制工程研究所 | 一种大偏心率轨道的星上快速高精度外推方法 |
CN104048664A (zh) * | 2014-07-01 | 2014-09-17 | 南京航空航天大学 | 一种导航卫星星座自主定轨的方法 |
CN104298647A (zh) * | 2014-09-30 | 2015-01-21 | 北京航空航天大学 | 基于低轨道地球卫星的地影时刻预报的星上确定方法 |
CN104501804A (zh) * | 2014-12-17 | 2015-04-08 | 深圳航天东方红海特卫星有限公司 | 一种基于gps测量数据的卫星在轨轨道预报方法 |
CN106124170A (zh) * | 2016-08-26 | 2016-11-16 | 上海航天控制技术研究所 | 一种基于高精度姿态信息的相机光轴指向计算方法 |
CN107451656A (zh) * | 2017-08-08 | 2017-12-08 | 上海航天控制技术研究所 | 一种火星探测器器上自主轨道计算方法 |
CN109059937A (zh) * | 2018-08-21 | 2018-12-21 | 深圳市天智运控科技有限公司 | 一种星上自主卫星轨道预报方法 |
CN110132261A (zh) * | 2018-11-16 | 2019-08-16 | 中国西安卫星测控中心 | 一种基于数值拟合的高精度星上轨道预报方法 |
CN109774974A (zh) * | 2019-01-31 | 2019-05-21 | 上海微小卫星工程中心 | 一种用于空间碎片接近的轨道设计方法 |
CN109990779A (zh) * | 2019-04-30 | 2019-07-09 | 桂林电子科技大学 | 一种惯性导航系统及方法 |
CN110595485A (zh) * | 2019-08-08 | 2019-12-20 | 上海航天电子有限公司 | 基于两行根数的低轨卫星长期轨道预报方法 |
CN112713922A (zh) * | 2019-10-25 | 2021-04-27 | 北京航空航天大学 | 一种多波束通讯卫星的可见性快速预报算法 |
CN111102981A (zh) * | 2019-12-09 | 2020-05-05 | 北京空间机电研究所 | 一种基于ukf的高精度卫星相对导航方法 |
CN111947667A (zh) * | 2020-06-24 | 2020-11-17 | 火眼位置数智科技服务有限公司 | 一种基于运动学和动力学组合的低轨卫星实时高精度定轨方法 |
CN112257343A (zh) * | 2020-10-22 | 2021-01-22 | 上海卫星工程研究所 | 一种高精度地面轨迹重复轨道优化方法及系统 |
CN113447043A (zh) * | 2021-05-21 | 2021-09-28 | 北京控制工程研究所 | 一种基于gnss的卫星天文导航系统误差自主标定方法及系统 |
CN113641949A (zh) * | 2021-08-05 | 2021-11-12 | 中国西安卫星测控中心 | 一种地球同步转移段轨道根数高精度拟合方法 |
CN113624243A (zh) * | 2021-08-16 | 2021-11-09 | 中国科学院上海天文台 | 一种近地轨道卫星的星上实时轨道预报方法 |
Non-Patent Citations (1)
Title |
---|
Equilibrium configurations of a four-body tethered system;CORRA A A,GMEZ G;《Journal of Guidance,Control,and Dynamics》;第29卷(第6期);1430-1435 * |
Also Published As
Publication number | Publication date |
---|---|
CN114383619A (zh) | 2022-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Giannitrapani et al. | Comparison of EKF and UKF for spacecraft localization via angle measurements | |
US6859727B2 (en) | Attitude change kalman filter measurement apparatus and method | |
US5909381A (en) | System of on board prediction of trajectories for autonomous navigation of GPS satellites | |
Gerlach et al. | CHAMP gravity field recovery using the energy balance approach | |
Abdelrahman et al. | Sigma-point Kalman filtering for spacecraft attitude and rate estimation using magnetometer measurements | |
Raol et al. | On the orbit determination problem | |
Montenbruck | An epoch state filter for use with analytical orbit models of low earth satellites | |
CN103148849A (zh) | 基于地月卫星联合测距和紫外敏感器的组合导航方法 | |
Shou | Orbit propagation and determination of low earth orbit satellites | |
CN110285815A (zh) | 一种可在轨全程应用的微纳卫星多源信息姿态确定方法 | |
CN112325886B (zh) | 一种基于重力梯度仪和陀螺仪组合的航天器自主定姿系统 | |
Noureldin et al. | Basic navigational mathematics, reference frames and the earth’s geometry | |
CN114383619B (zh) | 一种高精度轨道计算方法 | |
CN112304309B (zh) | 一种基于心动阵列的高超飞行器组合导航信息解算方法 | |
CN114802818A (zh) | 晨昏轨道卫星及其对日姿态计算方法、导引方法 | |
CN115127591A (zh) | 一种基于统计置信距离量测自举的机载dpos传递对准方法 | |
CN112327333A (zh) | 一种卫星位置计算方法及装置 | |
Reddy et al. | Transfer orbit stage gyrocompass alignment algorithm twist and sway environment for Mar's Observer mission on commercial Titan | |
Wu et al. | Strapdown INS/GPS integrated navigation using geometric algebra | |
Amert et al. | Hardware Demonstration and Improvements of the Stellar Positioning System | |
US20230152100A1 (en) | Navigation Methods Using Virtual Celestial Objects | |
CN111982128B (zh) | 一种陀螺安装矩阵的在轨标定方法及星载软件 | |
Vepa | Joint position localisation of spacecraft and debris for autonomous navigation applications using angle measurements only | |
Cha et al. | Coarse alignment of lunar exploration rover using accelerometer and sun sensor | |
Vijayanandh et al. | Prediction of spacecraft of position by Particle Filter based on GPS/INS integrated system |
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 |