CN113987407B - 一种地球静止卫星非球形引力切向漂移加速度计算方法 - Google Patents

一种地球静止卫星非球形引力切向漂移加速度计算方法 Download PDF

Info

Publication number
CN113987407B
CN113987407B CN202111260536.3A CN202111260536A CN113987407B CN 113987407 B CN113987407 B CN 113987407B CN 202111260536 A CN202111260536 A CN 202111260536A CN 113987407 B CN113987407 B CN 113987407B
Authority
CN
China
Prior art keywords
gravitational
geostationary satellite
drift
formula
drift acceleration
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
CN202111260536.3A
Other languages
English (en)
Other versions
CN113987407A (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 CN202111260536.3A priority Critical patent/CN113987407B/zh
Publication of CN113987407A publication Critical patent/CN113987407A/zh
Application granted granted Critical
Publication of CN113987407B publication Critical patent/CN113987407B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开的一种地球静止卫星非球形引力切向漂移加速度计算方法，其特征在于，包括以下步骤：建立地球非球形引力半长轴摄动方程，得到半长轴变化率；建立漂移加速度计算方程；将半长轴变化率代入漂移加速度计算方程，进而通过代入地球静止卫星当前平经度解析漂移加速度。本发明一种地球静止卫星非球形引力切向漂移加速度计算方法，不同于目前工程中使用的查表插值法，直接由地球非球形引力对静止卫星的半长轴摄动方程推导出解析计算公式，不依赖高精度漂移加速度表单文件，也避免了插值带来的误差，具有计算速度快、精度高的特点。

Description

一种地球静止卫星非球形引力切向漂移加速度计算方法
技术领域
本发明属于航天导航控制规划计算方法技术领域，具体涉及一种地球静止卫星非球形引力切向漂移加速度计算方法。
背景技术
由于地球非球形引力的作用，地球静止卫星始终受到一个切向漂移加速度的作用，因此半长轴会以恒定的速率增大或者减小，漂移加速度的大小因卫星当前的经度位置而异。半长轴变化导致轨道周期改变地球静止卫星不再“静止”，逐渐偏离定点经度，且随着半长轴变化量的累积，偏离速度(即漂移率的绝对值)会逐渐增大。地球静止卫星东西位保就是利用过调控制使卫星在尽可能长的一段时间里都位于一定经度范围内。
在地球静止轨道高度附近，卫星相对定点经度的漂移速度和卫星半长轴近似为线性关系，卫星半长轴随时间线性变化，因此卫星经度位置关于时间呈抛物线变化。所以在工程上，根据定点经度和经度保持范围设计抛物线保持环来进行东西位保计算。由公式推导易知，抛物线公式的二次项系数即为漂移加速度。因此定点经度切向加速度的计算精度直接决定东西位保控制计算的精度。目前工程上常用的方法是，在经度步长为1°的“经度——加速度”序列上进行插值计算，计算依赖序列文件，且精度较低。
发明内容
本发明的目的在于提供一种地球静止卫星非球形引力切向漂移加速度计算方法，具有计算速度快、精度高的特点。
本发明所采用的技术方案是：一种地球静止卫星非球形引力切向漂移加速度计算方法，包括以下步骤：
步骤1、建立地球非球形引力半长轴摄动方程，得到半长轴变化率；
步骤2、建立漂移加速度计算方程；
步骤3、将步骤1中得到的半长轴变化率代入步骤2中建立的漂移加速度计算方程，进而通过代入地球静止卫星当前平经度解析漂移加速度。
本发明的特点还在于，
步骤1中建立的地球非球形引力半长轴摄动方程为：
Figure BDA0003325257970000021
式(1)中，J22、J31、J33、J42、J44、l22、l31、l33、l42、l44为地球引力场4阶摄动系数，re为地球半径，a为轨道半长轴，l＝ω+M+Ω-Θ，l为平经度，ω为近地点幅角，M为平近点角，Ω为升交点赤经，Θ为格林尼治恒星时角，μ为地球引力常数。
步骤2中具体包括：
步骤2.1、设定漂移率公式：
Figure BDA0003325257970000022
式(2)中，ωe为地球自转角速度，ac为地球静止卫星标称半长轴；
步骤2.2、式(2)对时间t求导，得到漂移加速度计算公式：
Figure BDA0003325257970000031
式(3)中轨道半长轴a设定为标称值，即a＝ac。
步骤3具体包括将式(1)代入式(3)得到漂移加速度计算方程为：
Figure BDA0003325257970000032
式(4)中，代入地球静止卫星当前平经度l，即可解析计算出漂移加速度Dacc。
本发明的有益效果是：本发明一种地球静止卫星非球形引力切向漂移加速度计算方法，不同于目前工程中使用的查表插值法，直接由地球非球形引力对静止卫星的半长轴摄动方程推导出解析计算公式，不依赖高精度漂移加速度表单文件，也避免了插值带来的误差，具有计算速度快、精度高的特点。
附图说明
图1是本发明一种地球静止卫星非球形引力切向漂移加速度计算方法的流程图。
具体实施方式
下面结合附图以及具体实施方式对本发明进行详细说明。
本发明提供了一种地球静止卫星非球形引力切向漂移加速度计算方法，如图1所示，包括以下步骤：
1、建立地球非球形引力半长轴摄动方程
Figure BDA0003325257970000041
其中，J22、J31、J33、J42、J44、l22、l31、l33、l42、l44为地球引力场4阶摄动系数，re为地球半径，a为轨道半长轴，l＝ω+M+Ω-Θ，l为平经度，ω为近地点幅角，M为平近点角，Ω为升交点赤经，Θ为格林尼治恒星时角，μ为地球引力常数。
2、建立漂移加速度计算方程
漂移率公式
Figure BDA0003325257970000042
式(2)中，ωe为地球自转角速度，ac为地球静止卫星标称半长轴；
式(2)对时间t求导，得到漂移加速度计算公式
Figure BDA0003325257970000043
3、代入经度值解析计算
式(3)中，仅l一项变量，其他元素都是常量(半长轴设定为标称值，即a＝ac)，代入地球静止卫星当前平经度，即计算得到其漂移加速度。
以计算115°E的漂移加速度Dacc为例说明本方法的具体应用方式。
将式(1)带入式(3)得到，
Figure BDA0003325257970000051
带入变量l＝115°和常数l22＝-14.929°、l31＝6.968°、l33＝20.994°、l42＝31.024°、l44＝30.280°、J22＝1.811528e-6、J31＝2.2091169e-6、J33＝0.2213602e-6、J42＝1.6600442e-7、J44＝7.63937899e-9、re＝6378.1363、ac＝42165.700、μ＝398600.4415。
计算得到Dacc(115°E)＝-4.63093e-15，而通过查表插值法计算得到的115°E的漂移加速度为-4.6307e-15，两者差别极小，满足控制规划的经度要求。

Claims (3)

1.一种地球静止卫星非球形引力切向漂移加速度计算方法，其特征在于，包括以下步骤：
步骤1、建立地球非球形引力半长轴摄动方程，得到半长轴变化率；建立的地球非球形引力半长轴摄动方程为：
Figure FDA0003644590420000011
式(1)中，J22、J31、J33、J42、J44、l22、l31、l33、l42、l44为地球引力场4阶摄动系数，re为地球半径，a为轨道半长轴，l＝ω+M+Ω-Θ，l为平经度，ω为近地点幅角，M为平近点角，Ω为升交点赤经，Θ为格林尼治恒星时角，μ为地球引力常数；
步骤2、建立漂移加速度计算方程；具体包括：
步骤2.1、设定漂移率公式：
Figure FDA0003644590420000012
式(2)中，ωe为地球自转角速度，ac为地球静止卫星标称半长轴；
步骤2.2、式(2)对时间t求导，得到漂移加速度计算公式：
Figure FDA0003644590420000013
步骤3、将步骤1中得到的半长轴变化率代入步骤2中建立的漂移加速度计算方程，进而通过代入地球静止卫星当前平经度解析漂移加速度。
2.如权利要求1所述的一种地球静止卫星非球形引力切向漂移加速度计算方法，其特征在于，所述式(3)中轨道半长轴a设定为标称值，即a＝ac。
3.如权利要求2所述的一种地球静止卫星非球形引力切向漂移加速度计算方法，其特征在于，所述步骤3具体包括将式(1)代入式(3)得到漂移加速度计算方程为：
Figure FDA0003644590420000021
式(4)中，代入地球静止卫星当前平经度l，解析计算出漂移加速度Dacc。
CN202111260536.3A 2021-10-28 2021-10-28 一种地球静止卫星非球形引力切向漂移加速度计算方法 Active CN113987407B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111260536.3A CN113987407B (zh) 2021-10-28 2021-10-28 一种地球静止卫星非球形引力切向漂移加速度计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111260536.3A CN113987407B (zh) 2021-10-28 2021-10-28 一种地球静止卫星非球形引力切向漂移加速度计算方法

Publications (2)

Publication Number Publication Date
CN113987407A CN113987407A (zh) 2022-01-28
CN113987407B true CN113987407B (zh) 2022-07-08

Family

ID=79743130

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111260536.3A Active CN113987407B (zh) 2021-10-28 2021-10-28 一种地球静止卫星非球形引力切向漂移加速度计算方法

Country Status (1)

Country Link
CN (1) CN113987407B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001153682A (ja) * 1999-11-26 2001-06-08 Fujitsu Ltd 軌道計算装置および軌道計算方法
CN102935898A (zh) * 2012-11-30 2013-02-20 清华大学 空间纯引力轨道万有引力摄动在轨飞行验证方法
CN107179084A (zh) * 2017-06-27 2017-09-19 上海航天控制技术研究所 Gnss兼容机伪距与加表组合导航及漂移估计方法
CN108614575A (zh) * 2018-06-20 2018-10-02 北京电子工程总体研究所 一种卫星静止轨道定点位置调整方法
CN109255096A (zh) * 2018-07-25 2019-01-22 西北工业大学 一种基于微分代数的地球同步卫星轨道不确定演化方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7786929B2 (en) * 2007-11-26 2010-08-31 Mediatek Inc. Method and device for predicting GNSS satellite trajectory extension data used in mobile apparatus

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001153682A (ja) * 1999-11-26 2001-06-08 Fujitsu Ltd 軌道計算装置および軌道計算方法
CN102935898A (zh) * 2012-11-30 2013-02-20 清华大学 空间纯引力轨道万有引力摄动在轨飞行验证方法
CN107179084A (zh) * 2017-06-27 2017-09-19 上海航天控制技术研究所 Gnss兼容机伪距与加表组合导航及漂移估计方法
CN108614575A (zh) * 2018-06-20 2018-10-02 北京电子工程总体研究所 一种卫星静止轨道定点位置调整方法
CN109255096A (zh) * 2018-07-25 2019-01-22 西北工业大学 一种基于微分代数的地球同步卫星轨道不确定演化方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Orbital perturbation differential equations with non-linear corrections for champ-like satellite;Yu Jin-Hai et.al;《Chinese Journal of Geophysics》;20170531;全文 *
地球静止卫星电推进轨道保持策略优化;蒯政中 等;《中国空间科学技术》;20180630;全文 *

Also Published As

Publication number Publication date
CN113987407A (zh) 2022-01-28

Similar Documents

Publication Publication Date Title
CN110487301B (zh) 一种雷达辅助机载捷联惯性导航系统初始对准方法
CN110068845B (zh) 一种基于平根数理论确定卫星理论轨道的方法
CN109708663B (zh) 基于空天飞机sins辅助的星敏感器在线标定方法
CN108731674B (zh) 一种基于单轴旋转调制的惯性天文组合导航系统及计算方法
CN108875175A (zh) 一种高阶非中心引力场下的不变相对轨道初值确定方法
CN103925930B (zh) 一种重力仪双轴陀螺稳定平台航向误差效应的补偿方法
CN111268177B (zh) 一种静止轨道卫星分散式闭环自主位置保持控制方法
CN110727003A (zh) 一种北斗卫星导航系统的伪距仿真方法
CN112713922A (zh) 一种多波束通讯卫星的可见性快速预报算法
CN114684389B (zh) 考虑再入约束的月地转移窗口及精确转移轨道确定方法
CN103955005B (zh) 一种火箭橇轨道重力实时测量方法
CN113987407B (zh) 一种地球静止卫星非球形引力切向漂移加速度计算方法
CN112393835A (zh) 一种基于扩展卡尔曼滤波的小卫星在轨推力标定方法
CN115993777A (zh) 基于轨道摄动模型反演的径切联控解耦迭代标定方法
CN115542363A (zh) 一种适用于垂直下视航空吊舱的姿态测量方法
CN115510378A (zh) 地球静止轨道标称半长轴摄动计算方法
CN113218390B (zh) 基于姿态和恒星高度角融合的旋转惯性天文组合导航方法
CN114771873A (zh) 一种超低轨卫星轨道自主精确维持方法
Mikrin et al. Circumlunar spacecraft navigation using the measurements from global navigation satellite systems glonass, gps, galileo and beidou
CN103900569A (zh) 微惯导与dgps和电子罗盘组合导航姿态测量方法
CN104063582B (zh) 一种面外面内分步实施的绕飞构型建立方法
CN111547274A (zh) 一种航天器高精度自主目标预报方法
Braun et al. Performance analysis of an IMU-augmented GNSS tracking system on board the MAIUS-1 sounding rocket
CN111090830A (zh) 一种高轨非合作目标在轨光压辨识方法
CN114547549B (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