CN116552812B - 一种电推进geo卫星自学习定轨方法 - Google Patents
一种电推进geo卫星自学习定轨方法 Download PDFInfo
- Publication number
- CN116552812B CN116552812B CN202310390472.1A CN202310390472A CN116552812B CN 116552812 B CN116552812 B CN 116552812B CN 202310390472 A CN202310390472 A CN 202310390472A CN 116552812 B CN116552812 B CN 116552812B
- Authority
- CN
- China
- Prior art keywords
- geo satellite
- electric propulsion
- orbit
- state variable
- satellite
- 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 58
- 230000004927 fusion Effects 0.000 claims abstract description 13
- 239000013598 vector Substances 0.000 claims description 66
- 230000001133 acceleration Effects 0.000 claims description 46
- 239000011159 matrix material Substances 0.000 claims description 42
- 238000005259 measurement Methods 0.000 claims description 20
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 230000007774 longterm Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000001914 filtration Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 241000127225 Enceliopsis nudicaulis Species 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/242—Orbits and trajectories
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/10—Artificial satellites; Systems of such satellites; Interplanetary vehicles
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
-
- 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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
- G01S19/47—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being an inertial measurement, e.g. tightly coupled inertial
-
- 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
Abstract
本发明公开了一种电推进GEO卫星自学习定轨方法,本发明对电推进GEO卫星通过加速度计对卫星的加速度进行测量,将卫星的控制加速度信息引入系统状态方程,减小轨道动力学模型误差,使轨道递推结果更精确。当有GNSS定轨信息时,本发明利用LSTM算法对GNSS导航的定轨误差进行估计,然后与UKF定轨结果进行二次融合,进一步提升卫星自主定轨精度。当GNSS不可用时,本发明在不额外增加星载敏感器的条件下,利用卫星上的太阳帆板输出功率信息,获取真实的地影进入时刻和离开时刻,从而构建地影信息。根据卫星飞行每圈的进出地影时间,构建卫星真经度角速度补偿量,从而减小轨道动力学模型误差。
Description
技术领域
本发明涉及卫星定轨领域,具体涉及一种电推进GEO卫星自学习定轨方法。
背景技术
传统的卫星定轨主要依赖地面测控,然而,电推进GEO卫星在入轨过程中需要进行长期持续轨道机动,这使得基于地面测控的定轨方式具有诸多局限性。其一,地面测控站资源有限,对于电推进GEO卫星的轨道转移全过程,难以做到实时全程监控;其二,地面测控对于长期持续机动卫星的定轨精度不够高,已经严重制约了电推进GEO卫星的最终入轨定点捕获精度;其三,一旦地面测控站受到干扰或破坏,卫星将难以正常的飞行和执行预定的任务。因此,卫星自主定轨技术是未来电推进GEO卫星定轨的必然发展方向。自主定轨技术是航天器自主运行的核心技术,是航天器实现轨道姿态自主控制、精确定点捕获和执行在轨服务等空间任务的前提。拥有自主定轨能力的电推进GEO卫星不需要通过全球布站来提供全过程定轨支持,从而可以大幅减轻地面测控的负担,提高定轨的实时性、准确性和系统的安全性。
卫星自主定轨是指卫星利用自身携带的测量设备来实时确定自己的位置和速度的一项技术。惯性导航系统依靠携带的IMU(惯性导航单元),通过测量卫星的加速度,并根据轨道动力学模型进行积分运算,获得卫星瞬时速度和瞬时位置数据。惯性导航具有完全自主的优势,但是随着运行时间增加,轨道动力学模型积分导致的系统误差将不断增大,因此需要借助其他导航信息修正。利用全球卫星导航系统,可通过GEO卫星携带的GNSS接收机接受GNSS信号(如北斗、GPS的信号等),从而解算获得卫星的位置和速度信息。利用GNSS的定轨能够提供高精度的导航信息,但是相比惯性导航,其信息输出频率较低,且受到导航卫星可见性限制,GNSS信号对高轨卫星难以全覆盖,在高动态条件下不能保证接收到足够数量的GNSS信号。
然而,当前的基于组合导航系统的自主定轨方法大多适用于不进行机动的卫星。若卫星进行长期持续轨道机动,则轨道动力学模型积分所得数据的误差累积,严重影响定轨精度。此外,当前导航滤波算法大多基于卡尔曼滤波理论设计,例如EKF(扩展卡尔曼滤波)、UKF(无迹卡尔曼滤波)等。基于卡尔曼滤波理论设计的导航滤波器依赖于先验已知的系统噪声方差和测量噪声方差。但在实际的卫星运行环境中,受地球非球形引力摄动、大气阻力摄动、太阳光压摄动等多种摄动、空间热环境变化、卫星机动导致的观测几何变化等影响,系统噪声和测量噪声的统计特性存在不确定性,这导致传统导航滤波算法的精度上限难以进一步提升。
发明内容
针对现有技术中的上述不足,本发明提供的一种电推进GEO卫星自学习定轨方法解决了现有电推进GEO卫星自主定轨精度低的问题。
为了达到上述发明目的,本发明采用的技术方案为:
提供一种电推进GEO卫星自学习定轨方法,其包括以下步骤:
S1、对电推进GEO卫星建立带动力改进春分点轨道根数模型;获取电推进GEO卫星的加速度信息;
S2、根据带动力改进春分点轨道根数模型获取电推进GEO卫星以位置和速度矢量为状态变量的系统状态方程;
S3、根据加速度信息和系统状态方程对位置和速度矢量进行时间更新,获取电推进GEO卫星的理论轨道信息;
S4、判断电推进GEO卫星是否接收到足够的GNSS信号,若是则进入步骤S5;否则进入步骤S8;
S5、基于GNSS信号,采用UKF算法进行第一次多源信息融合,获取电推进GEO卫星的第一状态变量估计值和第一状态变量协方差估计矩阵;
S6、将上一时刻的状态变量值和加速度信息输入LSTM网络,获取对应的误差预测值,并结合第一状态变量估计值和第一状态变量协方差估计矩阵进行第二次多源信息融合,获取电推进GEO卫星的第二状态变量估计值,并将其作为电推进GEO卫星的最终状态估计值;
S7、判断是否结束自学习定轨,若是则结束;否则返回步骤S3;
S8、根据星历和电推进GEO卫星上的太阳帆板输出功率信息,计算当前卫星飞行周期内的地影信息;
S9、通过地影信息分段补偿法修正带动力改进春分点轨道根数模型,得到修正后的系统状态方程;
S10、通过修正后的系统状态方程进行轨道递推,对位置和速度矢量进行时间更新,获取电推进GEO卫星的修正后理论轨道信息,并将其作为电推进GEO卫星的最终状态估计值;返回步骤S7。
进一步地,步骤S1中带动力改进春分点轨道根数模型为:
其中[p,f,g,h,k,L]T为带动力改进春分点轨道根数模型的状态变量,p为轨道半通径,f和g为偏心率矢量分量,h和k为轨道倾角矢量分量,L为真经度,[·]T表示矩阵的转置;μ为地球引力常数;w为常数,w=1+fcosL+gsinL;ft为电推进GEO卫星的切向加速度;fr为电推进GEO卫星的径向加速度;fn为电推进GEO卫星的法向加速度;s为常数,s2=1+h2+k2;m为电推进GEO卫星质量;T为电推进GEO卫星的推力大小;g0为海平面处重力加速度;Isp为电推进系统的比冲;字符上一点表示求导。
进一步地,步骤S2的具体方法为:
根据变量转换关系将改进春分点轨道根数模型的状态变量[p,f,g,h,k,L]T转换为位置、速度矢量[r,v]T,得到以位置和速度矢量为状态变量的系统状态方程。
进一步地,步骤S5的具体方法包括以下子步骤:
S5-1、根据公式:
获取状态初始值x0的期望值和状态变量协方差矩阵初始值P0,并将期望值/>作为初始状态变量的均值;其中(·)T表示矩阵的转置;E(·)表示求期望;
S5-2、根据公式:
构造2n+1个sigma点;其中χk为k时刻的sigma矢量;为k时刻状态变量的均值;Pk xx为k时刻状态变量的协方差矩阵;n为状态向量长度;λ为比例参数;/>表示Pk xx的Cholesky分解;
S5-3、根据公式:
计算每个sigma点对应的权重;其中Wi m为计算均值的权重系数;Wi c为计算协方差的权重系数;α和κ均为尺度参数;
S5-4、根据公式:
建立预测方程;其中为k+1时刻状态变量的sigma矢量;Φ(χi,k,ui,k)表示非线性系统状态转移方程;χi,k为k时刻状态变量的sigma矢量;ui,k表示k时刻的控制参数;/>为k+1时刻通过sigma矢量得到的状态变量的先验均值;/>为k+1时刻通过sigma矢量得到的状态变量的先验协方差矩阵;Q为系统噪声协方差矩阵;
S5-5、基于GNSS信号获取k+1时刻的状态变量的先验均值
S5-6、根据公式:
建立测量修正方程;其中H(·)为观测方程,Yi,k+1为k+1时刻对测量得到的状态变量进行转移得到的sigma矢量;为k+1时刻测量得到的状态变量的先验协方差矩阵;R为量测噪声协方差矩阵;/>为k+1时刻的交叉协方差阵;Kk+1为k+1时刻的卡尔曼增益系数;为k+1时刻更新后的状态变量值,即第一状态变量估计值;yk+1为k+1时刻通过GNSS测量得到的状态变量;/>为k+1时刻更新后的状态变量协方差矩阵;/>为Kk+1的转置。
进一步地,步骤S6的具体方法包括以下子步骤:
S6-1、将上一时刻的卫星位置、速度矢量和加速度信息输入至训练好的LSTM模型中,获取LSTM模型预测的状态变量值和LSTM模型预测的定轨误差估计向量;
S6-2、根据公式:
进行第二次多源信息融合,获取电推进GEO卫星的第二状态变量估计值并将其作为电推进GEO卫星的最终状态估计值;其中/>为电推进GEO卫星的第一状态变量估计值;PUKF(tk)为tk时刻下UKF算法中的状态变量的协方差矩阵;PLSTM(tk)为tk时刻下LSTM模型预测得到的定轨误差估计向量;/>为tk时刻下LSTM模型预测得到的状态变量值。
进一步地,步骤S8的具体方法包括以下子步骤:
S8-1、以地心为起点,以太阳光线照射方向为方向作射线,记为地影中线;
S8-2、根据电推进GEO卫星上的太阳帆板输出功率信息,获取真实的地影进入时刻ten和离开时刻tex;
S8-3、根据公式:
获取电推进GEO卫星到达地影中线的时刻tmid;
S8-4、根据公式:
获取地影中线真经度其中L0为初始状态的真经度;w*为地影中线转动的角速度。
进一步地,步骤S9的具体方法包括以下子步骤:
S9-1、根据公式:
获取电推进GEO卫星绕地球旋转的当前圈的真实真经度平均角速度和电推进GEO卫星绕地球旋转的当前圈的估计真经度角速度;其中/>为电推进GEO卫星绕地球旋转的上一圈数据对应的地影中线真经度;/>电推进GEO卫星绕地球旋转的上一圈数据对应的地影中线时刻;Les为通过带动力改进春分点轨道根数模型递推获得的对应时刻的卫星真经度估计值;
S9-2、根据公式:
获取电推进GEO卫星绕地球旋转的当前圈的真经度角速度补偿值wc;其中为电推进GEO卫星绕地球旋转的上一圈的真经度角速度补偿值;
S9-3、根据公式:
修正带动力改进春分点轨道根数模型,得到修正后的系统状态方程;其中p为轨道半通径,f和g为偏心率矢量分量,h和k为轨道倾角矢量分量,L为真经度,[·]T表示矩阵的转置;μ为地球引力常数;w为常数,w=1+fcosL+gsinL;ft为电推进GEO卫星的切向加速度;fr为电推进GEO卫星的径向加速度;fn为电推进GEO卫星的法向加速度;s为常数,s2=1+h2+k2;m为电推进GEO卫星的质量;T为电推进GEO卫星的推力大小;g0为海平面处重力加速度;Isp为电推进系统的比冲;字符上一点表示求导。
本发明的有益效果为:
1、本发明对电推进GEO卫星通过加速度计对卫星的加速度进行测量,将卫星的控制加速度信息引入系统状态方程,减小轨道动力学模型误差,使轨道递推结果更精确。
2、当有GNSS定轨信息时,本发明利用LSTM算法对GNSS导航的定轨误差进行估计,然后将UKF融合结果进行二次融合,进一步提升卫星自主定轨精度。
3、当GNSS不可用时,本发明在不额外增加星载敏感器的条件下,利用卫星上的太阳帆板输出功率信息,获取真实的地影进入时刻和离开时刻,从而构建地影信息。根据卫星飞行每圈的进出地影时间,构建卫星真经度角速度补偿量,从而减小轨道动力学模型误差。
附图说明
图1为本方法的流程示意图;
图2为地球阴影区域示意图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
如图1所示,该电推进GEO卫星自学习定轨方法包括以下步骤:
S1、对电推进GEO卫星建立带动力改进春分点轨道根数模型;获取电推进GEO卫星的加速度信息;
S2、根据带动力改进春分点轨道根数模型获取电推进GEO卫星以位置和速度矢量为状态变量的系统状态方程;
S3、根据加速度信息和系统状态方程对位置和速度矢量进行时间更新,获取电推进GEO卫星的理论轨道信息;
S4、判断电推进GEO卫星是否接收到足够的GNSS信号,若是则进入步骤S5;否则进入步骤S8;
S5、基于GNSS信号,采用UKF算法进行第一次多源信息融合,获取电推进GEO卫星的第一状态变量估计值和第一状态变量协方差估计矩阵;
S6、将上一时刻的状态变量值和加速度信息输入LSTM网络,获取对应的误差预测值,并结合第一状态变量估计值和第一状态变量协方差估计矩阵进行第二次多源信息融合,获取电推进GEO卫星的第二状态变量估计值,并将其作为电推进GEO卫星的最终状态估计值;
S7、判断是否结束自学习定轨,若是则结束;否则返回步骤S3;
S8、根据星历和电推进GEO卫星上的太阳帆板输出功率信息,计算当前卫星飞行周期内的地影信息;
S9、通过地影信息分段补偿法修正带动力改进春分点轨道根数模型,得到修正后的系统状态方程;
S10、通过修正后的系统状态方程进行轨道递推,对位置和速度矢量进行时间更新,获取电推进GEO卫星的修正后理论轨道信息,并将其作为电推进GEO卫星的最终状态估计值;返回步骤S7。
步骤S1中带动力改进春分点轨道根数模型为:
其中[p,f,g,h,k,L]T为带动力改进春分点轨道根数模型的状态变量,p为轨道半通径,f和g为偏心率矢量分量,h和k为轨道倾角矢量分量,L为真经度,[·]T表示矩阵的转置;μ为地球引力常数;w为常数,w=1+fcosL+gsinL;ft为电推进GEO卫星的切向加速度;fr为电推进GEO卫星的径向加速度;fn为电推进GEO卫星的法向加速度;s为常数,s2=1+h2+k2;m为电推进GEO卫星质量;T为电推进GEO卫星的推力大小;g0为海平面处重力加速度;Isp为电推进系统的比冲;字符上一点表示求导。
步骤S2的具体方法为:根据变量转换关系将改进春分点轨道根数模型的状态变量[p,f,g,h,k,L]T转换为位置、速度矢量[r,v]T,得到以位置和速度矢量为状态变量的系统状态方程。
步骤S5的具体方法包括以下子步骤:
S5-1、根据公式:
获取状态初始值x0的期望值和状态变量协方差矩阵初始值P0,并将期望值/>作为初始状态变量的均值;其中(·)T表示矩阵的转置;E(·)表示求期望;
S5-2、根据公式:
构造2n+1个sigma点;其中χk为k时刻的sigma矢量;为k时刻状态变量的均值;Pk xx为k时刻状态变量的协方差矩阵;n为状态向量长度; λ为比例参数;/>表示Pk xx的Cholesky分解;
S5-3、根据公式:
计算每个sigma点对应的权重;其中Wi m为计算均值的权重系数;Wi c为计算协方差的权重系数;α和κ均为尺度参数;β为常数;i表示sigma点的序号;α用于确定周围sigma点的分布情况,一般取一小正数(0.0001≤α≤1);κ通常为0或3-n;β是一个与状态的先验分布信息有关的参数,对于高斯分布β=2最优。
S5-4、根据公式:
建立预测方程;其中为k+1时刻状态变量的sigma矢量;Φ(χi,k,ui,k)表示非线性系统状态转移方程;χi,k为k时刻状态变量的sigma矢量;ui,k表示k时刻的控制参数;/>为k+1时刻通过sigma矢量得到的状态变量的先验均值;/>为k+1时刻通过sigma矢量得到的状态变量的先验协方差矩阵;Q为系统噪声协方差矩阵;
S5-5、基于GNSS信号获取k+1时刻的状态变量的先验均值
S5-6、根据公式:
建立测量修正方程;其中H(·)为观测方程,Yi,k+1为k+1时刻对测量得到的状态变量进行转移得到的sigma矢量;为k+1时刻测量得到的状态变量的先验协方差矩阵;R为量测噪声协方差矩阵;/>为k+1时刻的交叉协方差阵;Kk+1为k+1时刻的卡尔曼增益系数;为k+1时刻更新后的状态变量值,即第一状态变量估计值;yk+1为k+1时刻通过GNSS测量得到的状态变量;/>为k+1时刻更新后的状态变量协方差矩阵;/>为Kk+1的转置。
步骤S6的具体方法包括以下子步骤:
S6-1、将上一时刻的卫星位置、速度矢量和加速度信息输入至训练好的LSTM模型中,获取LSTM模型预测的状态变量值和LSTM模型预测的定轨误差估计向量;
S6-2、根据公式:
进行第二次多源信息融合,获取电推进GEO卫星的第二状态变量估计值并将其作为电推进GEO卫星的最终状态估计值;其中/>为电推进GEO卫星的第一状态变量估计值;PUKF(tk)为tk时刻下UKF算法中的状态变量的协方差矩阵;PLSTM(tk)为tk时刻下LSTM模型预测得到的定轨误差估计向量;/>为tk时刻下LSTM模型预测得到的状态变量值。
步骤S8的具体方法包括以下子步骤:
S8-1、以地心为起点,以太阳光线照射方向为方向作射线,记为地影中线;
S8-2、根据电推进GEO卫星上的太阳帆板输出功率信息,获取真实的地影进入时刻ten和离开时刻tex;
S8-3、根据公式:
获取电推进GEO卫星到达地影中线的时刻tmid;
S8-4、根据公式:
获取地影中线真经度其中L0为初始状态的真经度;w*为地影中线转动的角速度。
步骤S9的具体方法包括以下子步骤:
S9-1、根据公式:
获取电推进GEO卫星绕地球旋转的当前圈的真实真经度平均角速度和电推进GEO卫星绕地球旋转的当前圈的估计真经度角速度;其中/>为电推进GEO卫星绕地球旋转的上一圈数据对应的地影中线真经度;/>电推进GEO卫星绕地球旋转的上一圈数据对应的地影中线时刻;Les为通过带动力改进春分点轨道根数模型递推获得的对应时刻的卫星真经度估计值;
S9-2、根据公式:
获取电推进GEO卫星绕地球旋转的当前圈的真经度角速度补偿值wc;其中为电推进GEO卫星绕地球旋转的上一圈的真经度角速度补偿值;
S9-3、根据公式:
/>
修正带动力改进春分点轨道根数模型,得到修正后的系统状态方程;其中p为轨道半通径,f和g为偏心率矢量分量,h和k为轨道倾角矢量分量,L为真经度,[·]T表示矩阵的转置;μ为地球引力常数;w为常数,w=1+fcosL+gsinL;ft为电推进GEO卫星的切向加速度;fr为电推进GEO卫星的径向加速度;fn为电推进GEO卫星的法向加速度;s为常数,s2=1+h2+k2;m为电推进GEO卫星的质量;T为电推进GEO卫星的推力大小;g0为海平面处重力加速度;Isp为电推进系统的比冲;字符上一点表示求导。
在具体实施过程中,电推进GEO卫星三个方向上的加速度可表示如下:
fr=fre+frJ2
ft=fte+ftJ2
fn=fne+fnJ2
其中fe=[fre,fte,fne]T表示电推力产生的加速度矢量,fJ2=[frJ2,ftJ2,fnJ2]T表示J2摄动产生的加速度。
步骤S3中加速度信息的量测方程为:
ab=at+ac+as
其中ab为卫星本体系下测量的加速度矢量;at为真实的加速度矢量;ac为加速度计的常值漂移,可以通过在轨标定的方式获得较为准确的常漂估计值;as为测量噪声。
步骤S3中时间更新的方法可以采用龙格库塔45数值积分方法,也可以采用其他数值积分方法。
步骤S4中是否接收到足够的GNSS信号的判断标准为是否接收到以解算GNSS定位数据的GNSS信号。
LSTM网络在使用之前可以先进行训练,训练的方式为:构建样本数据集,包括卫星的位置矢量(3维)、速度矢量(3维)和加速度计测量得到的卫星测量加速度矢量(3维),采用样本数据集对LSTM网络进行训练,并使LSTM网络的输出状态变量值定轨误差估计向量。
如图2所示,地影信息是地球围绕太阳公转而产生的周期性的阴影区域,Len和Lex分别为卫星进入地影区域和离开地影区域的真赤经。本方法定义一条地影中线,是以地心为起点,太阳光线照射方向为方向的一条射线,Lmid为卫星到达地影中线的真赤经。由于电推进GEO卫星在阴影区域的时间短,轨道变化不明显,因此视为真经度L均匀变化。由于地球绕太阳公转,周期为1年,则可视为均速绕行转动。地影中线的真赤经可简化为
综上所述,本发明提出了一种结合LSTM算法的基于GNSS和惯性导航的组合导航方法,该方法利用UKF算法和LSTM进行两次多源信息融合,以减小电推进GEO卫星自主定轨误差。并针对GNSS信号弱或不可用的场景,提出了利用地影信息修正轨道积分模型的方法,实现在GNSS信号弱或GNSS卫星可见性差条件下的电推进GEO卫星自主定轨。
Claims (7)
1.一种电推进GEO卫星自学习定轨方法,其特征在于,包括以下步骤:
S1、对电推进GEO卫星建立带动力改进春分点轨道根数模型;获取电推进GEO卫星的加速度信息;
S2、根据带动力改进春分点轨道根数模型获取电推进GEO卫星以位置和速度矢量为状态变量的系统状态方程;
S3、根据加速度信息和系统状态方程对位置和速度矢量进行时间更新,获取电推进GEO卫星的理论轨道信息;
S4、判断电推进GEO卫星是否接收到足够的GNSS信号,若是则进入步骤S5;否则进入步骤S8;
S5、基于GNSS信号,采用UKF算法进行第一次多源信息融合,获取电推进GEO卫星的第一状态变量估计值和第一状态变量协方差估计矩阵;
S6、将上一时刻的状态变量值和加速度信息输入LSTM网络,获取对应的误差预测值,并结合第一状态变量估计值和第一状态变量协方差估计矩阵进行第二次多源信息融合,获取电推进GEO卫星的第二状态变量估计值,并将其作为电推进GEO卫星的最终状态估计值;
S7、判断是否结束自学习定轨,若是则结束;否则返回步骤S3;
S8、根据星历和电推进GEO卫星上的太阳帆板输出功率信息,计算当前卫星飞行周期内的地影信息;
S9、通过地影信息分段补偿法修正带动力改进春分点轨道根数模型,得到修正后的系统状态方程;
S10、通过修正后的系统状态方程进行轨道递推,对位置和速度矢量进行时间更新,获取电推进GEO卫星的修正后理论轨道信息,并将其作为电推进GEO卫星的最终状态估计值;返回步骤S7。
2.根据权利要求1所述的电推进GEO卫星自学习定轨方法,其特征在于,步骤S1中带动力改进春分点轨道根数模型为:
其中[p,f,g,h,k,L]T为带动力改进春分点轨道根数模型的状态变量,p为轨道半通径,f和g为偏心率矢量分量,h和k为轨道倾角矢量分量,L为真经度,[·]T表示矩阵的转置;μ为地球引力常数;w为常数,w=1+fcosL+gsinL;ft为电推进GEO卫星的切向加速度;fr为电推进GEO卫星的径向加速度;fn为电推进GEO卫星的法向加速度;s为常数,s2=1+h2+k2;m为电推进GEO卫星质量;T为电推进GEO卫星的推力大小;g0为海平面处重力加速度;Isp为电推进系统的比冲;字符上一点表示求导。
3.根据权利要求2所述的电推进GEO卫星自学习定轨方法,其特征在于,步骤S2的具体方法为:
根据变量转换关系将改进春分点轨道根数模型的状态变量[p,f,g,h,k,L]T转换为位置、速度矢量[r,v]T,得到以位置和速度矢量为状态变量的系统状态方程。
4.根据权利要求1所述的电推进GEO卫星自学习定轨方法,其特征在于,步骤S5的具体方法包括以下子步骤:
S5-1、根据公式:
获取状态初始值x0的期望值和状态变量协方差矩阵初始值P0,并将期望值/>作为初始状态变量的均值;其中(·)T表示矩阵的转置;E(·)表示求期望;
S5-2、根据公式:
构造2n+1个sigma点;其中χk为k时刻的sigma矢量;为k时刻状态变量的均值;Pk xx为k时刻状态变量的协方差矩阵;n为状态向量长度;λ为比例参数;/>表示Pk xx的Cholesky分解;
S5-3、根据公式:
计算每个sigma点对应的权重;其中Wi m为计算均值的权重系数;Wi c为计算协方差的权重系数;α和κ均为尺度参数;
S5-4、根据公式:
建立预测方程;其中为k+1时刻状态变量的sigma矢量;Φ(χi,k,ui,k)表示非线性系统状态转移方程;χi,k为k时刻状态变量的sigma矢量;ui,k表示k时刻的控制参数;/>为k+1时刻通过sigma矢量得到的状态变量的先验均值;/>为k+1时刻通过sigma矢量得到的状态变量的先验协方差矩阵;Q为系统噪声协方差矩阵;
S5-5、基于GNSS信号获取k+1时刻的状态变量的先验均值
S5-6、根据公式:
建立测量修正方程;其中H(·)为观测方程,Yi,k+1为k+1时刻对测量得到的状态变量进行转移得到的sigma矢量;为k+1时刻测量得到的状态变量的先验协方差矩阵;R为量测噪声协方差矩阵;/>为k+1时刻的交叉协方差阵;Kk+1为k+1时刻的卡尔曼增益系数;/>为k+1时刻更新后的状态变量值,即第一状态变量估计值;yk+1为k+1时刻通过GNSS测量得到的状态变量;/>为k+1时刻更新后的状态变量协方差矩阵;/>为Kk+1的转置。
5.根据权利要求1所述的电推进GEO卫星自学习定轨方法,其特征在于,步骤S6的具体方法包括以下子步骤:
S6-1、将上一时刻的卫星位置、速度矢量和加速度信息输入至训练好的LSTM模型中,获取LSTM模型预测的状态变量值和LSTM模型预测的定轨误差估计向量;
S6-2、根据公式:
进行第二次多源信息融合,获取电推进GEO卫星的第二状态变量估计值并将其作为电推进GEO卫星的最终状态估计值;其中/>为电推进GEO卫星的第一状态变量估计值;PUKF(tk)为tk时刻下UKF算法中的状态变量的协方差矩阵;PLSTM(tk)为tk时刻下LSTM模型预测得到的定轨误差估计向量;/>为tk时刻下LSTM模型预测得到的状态变量值。
6.根据权利要求1所述的电推进GEO卫星自学习定轨方法,其特征在于,步骤S8的具体方法包括以下子步骤:
S8-1、以地心为起点,以太阳光线照射方向为方向作射线,记为地影中线;
S8-2、根据电推进GEO卫星上的太阳帆板输出功率信息,获取真实的地影进入时刻ten和离开时刻tex;
S8-3、根据公式:
获取电推进GEO卫星到达地影中线的时刻tmid;
S8-4、根据公式:
获取地影中线真经度其中L0为初始状态的真经度;w*为地影中线转动的角速度。
7.根据权利要求6所述的电推进GEO卫星自学习定轨方法,其特征在于,步骤S9的具体方法包括以下子步骤:
S9-1、根据公式:
获取电推进GEO卫星绕地球旋转的当前圈的真实真经度平均角速度和电推进GEO卫星绕地球旋转的当前圈的估计真经度角速度;其中/>为电推进GEO卫星绕地球旋转的上一圈数据对应的地影中线真经度;/>电推进GEO卫星绕地球旋转的上一圈数据对应的地影中线时刻;Les为通过带动力改进春分点轨道根数模型递推获得的对应时刻的卫星真经度估计值;
S9-2、根据公式:
获取电推进GEO卫星绕地球旋转的当前圈的真经度角速度补偿值wc;其中为电推进GEO卫星绕地球旋转的上一圈的真经度角速度补偿值;
S9-3、根据公式:
修正带动力改进春分点轨道根数模型,得到修正后的系统状态方程;其中p为轨道半通径,f和g为偏心率矢量分量,h和k为轨道倾角矢量分量,L为真经度,[·]T表示矩阵的转置;μ为地球引力常数;w为常数,w=1+fcosL+gsinL;ft为电推进GEO卫星的切向加速度;fr为电推进GEO卫星的径向加速度;fn为电推进GEO卫星的法向加速度;s为常数,s2=1+h2+k2;m为电推进GEO卫星的质量;T为电推进GEO卫星的推力大小;g0为海平面处重力加速度;Isp为电推进系统的比冲;字符上一点表示求导。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310390472.1A CN116552812B (zh) | 2023-04-12 | 2023-04-12 | 一种电推进geo卫星自学习定轨方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310390472.1A CN116552812B (zh) | 2023-04-12 | 2023-04-12 | 一种电推进geo卫星自学习定轨方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116552812A CN116552812A (zh) | 2023-08-08 |
CN116552812B true CN116552812B (zh) | 2024-01-23 |
Family
ID=87495544
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310390472.1A Active CN116552812B (zh) | 2023-04-12 | 2023-04-12 | 一种电推进geo卫星自学习定轨方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116552812B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675861A (zh) * | 2013-11-18 | 2014-03-26 | 航天恒星科技有限公司 | 一种基于星载gnss多天线的卫星自主定轨方法 |
CN108120994A (zh) * | 2017-10-30 | 2018-06-05 | 千寻位置网络(浙江)有限公司 | 一种基于星载gnss的geo卫星实时定轨方法 |
CN108692729A (zh) * | 2018-05-04 | 2018-10-23 | 北京空间飞行器总体设计部 | 一种空间非合作目标相对导航协方差自适应修正滤波方法 |
CN109677636A (zh) * | 2018-12-11 | 2019-04-26 | 上海航天控制技术研究所 | 一种高精度双星编队对地观测轨道控制方法 |
CN110058287A (zh) * | 2019-05-16 | 2019-07-26 | 北京合众思壮科技股份有限公司 | 一种低轨卫星定轨方法、装置及系统 |
CN111942618A (zh) * | 2020-07-08 | 2020-11-17 | 北京控制工程研究所 | 一种适用于动中成像的基于gnss数据的轨道获取方法 |
CN115327587A (zh) * | 2022-07-29 | 2022-11-11 | 北京自动化控制设备研究所 | 基于gnss定位信息的低轨卫星轨道误差修正方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106291645B (zh) * | 2016-07-19 | 2018-08-21 | 东南大学 | 适于高维gnss/ins深耦合的容积卡尔曼滤波方法 |
-
2023
- 2023-04-12 CN CN202310390472.1A patent/CN116552812B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675861A (zh) * | 2013-11-18 | 2014-03-26 | 航天恒星科技有限公司 | 一种基于星载gnss多天线的卫星自主定轨方法 |
CN108120994A (zh) * | 2017-10-30 | 2018-06-05 | 千寻位置网络(浙江)有限公司 | 一种基于星载gnss的geo卫星实时定轨方法 |
CN108692729A (zh) * | 2018-05-04 | 2018-10-23 | 北京空间飞行器总体设计部 | 一种空间非合作目标相对导航协方差自适应修正滤波方法 |
CN109677636A (zh) * | 2018-12-11 | 2019-04-26 | 上海航天控制技术研究所 | 一种高精度双星编队对地观测轨道控制方法 |
CN110058287A (zh) * | 2019-05-16 | 2019-07-26 | 北京合众思壮科技股份有限公司 | 一种低轨卫星定轨方法、装置及系统 |
CN111942618A (zh) * | 2020-07-08 | 2020-11-17 | 北京控制工程研究所 | 一种适用于动中成像的基于gnss数据的轨道获取方法 |
CN115327587A (zh) * | 2022-07-29 | 2022-11-11 | 北京自动化控制设备研究所 | 基于gnss定位信息的低轨卫星轨道误差修正方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN116552812A (zh) | 2023-08-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111680354B (zh) | 近地回归轨道卫星星下点和摄影点轨迹自交点的计算方法 | |
CN107797130B (zh) | 低轨航天器多点多参数轨道上行数据计算方法 | |
CN102878995B (zh) | 一种静止轨道卫星自主导航方法 | |
CN110567454B (zh) | 一种复杂环境下sins/dvl紧组合导航方法 | |
CN112161632B (zh) | 一种基于相对位置矢量测量的卫星编队初始定位方法 | |
CN110243377B (zh) | 一种基于分层式结构的集群飞行器协同导航方法 | |
CN111552003B (zh) | 基于球卫星编队的小行星引力场全自主测量系统及方法 | |
CN111156986B (zh) | 一种基于抗差自适应ukf的光谱红移自主组合导航方法 | |
CN111947653A (zh) | 一种月面巡视探测器双模式惯性/视觉/天文导航方法 | |
CN111121766A (zh) | 一种基于星光矢量的天文与惯性组合导航方法 | |
CN112713922A (zh) | 一种多波束通讯卫星的可见性快速预报算法 | |
CN111207773B (zh) | 一种用于仿生偏振光导航的姿态无约束优化求解方法 | |
CN110554443A (zh) | 基于载波相位观测值和点加速度法确定地球重力场的方法 | |
CN116552812B (zh) | 一种电推进geo卫星自学习定轨方法 | |
CN113008229A (zh) | 一种基于低成本车载传感器的分布式自主组合导航方法 | |
CN110388942B (zh) | 一种基于角度和速度增量的车载姿态精对准系统 | |
CN114771873B (zh) | 一种超低轨卫星轨道自主精确维持方法 | |
CN115420280A (zh) | 一种航天器脉冲星/星光折射组合导航方法 | |
Kuroda et al. | Accurate localization in combination with planet observation and dead reckoning for lunar rover | |
CN116125503A (zh) | 一种高精度卫星轨道确定及预报算法 | |
CN115392540A (zh) | 一种用于月球轨道交会制导的快速预报方法 | |
CN106767824B (zh) | 一种计算双探测器在地外天体表面相对位置的方法 | |
CN112393835B (zh) | 一种基于扩展卡尔曼滤波的小卫星在轨推力标定方法 | |
CN115014321A (zh) | 一种基于自适应鲁棒滤波的仿生偏振多源融合定向方法 | |
CN111308570B (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 |