CN109159923B - 一种航天器追踪目标丢失段的相对可达区域计算方法 - Google Patents
一种航天器追踪目标丢失段的相对可达区域计算方法 Download PDFInfo
- Publication number
- CN109159923B CN109159923B CN201810837425.6A CN201810837425A CN109159923B CN 109159923 B CN109159923 B CN 109159923B CN 201810837425 A CN201810837425 A CN 201810837425A CN 109159923 B CN109159923 B CN 109159923B
- Authority
- CN
- China
- Prior art keywords
- spacecraft
- target spacecraft
- equation
- target
- follows
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 12
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Images
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
-
- 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)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Remote Sensing (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Combustion & Propulsion (AREA)
- Operations Research (AREA)
- Aviation & Aerospace Engineering (AREA)
- Algebra (AREA)
- Chemical & Material Sciences (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种航天器追踪目标丢失段的相对可达区域计算方法,包括以下步骤:1)建立地球惯性坐标系EXYZ及目标航天器轨道坐标系Txyz,并构建航天器相对运动轨道模型,求解航天器相对运动轨道模型的解;2)根据航天器相对运动轨道模型,在任一时刻加入目标航天器施加的机动力,以构造新的航天器相对运动轨道模型;3)求解新的航天器相对运动轨道模型的解析解,然后将求解得到的解析解作为目标航天器丢失段的相对可达区域的方程,该方法能够计算得到航天器追踪目标丢失段的相对可达区域。
Description
技术领域
本发明属于空间航天器相对运动技术领域,涉及一种航天器追踪目标丢失段的相对可达区域计算方法。
背景技术
近几年,随着在轨服务任务的需求,对目标航天器的轨道预测变得越来越重要。追踪航天器时时掌握目标航天器在轨信息,是进一步对目标进行任务的前提。但当追踪航天器对目标进行观测过程中,由于太阳光、地球遮蔽等因素,在一段时间内是无法观测到目标的,如果在该时间段内,目标航天器实施机动时,追踪航天器将无法成功预测目标航天器下一时段的位置信息,这将导致目标的丢失。因此利用相对可达区域的概念,对目标丢失段进行相对可达区域的计算,对目标在该时间段可能到达的区域进行描述,有利于追踪航天器对目标位置信息的把握。
对可达区域的研究是对航天器机动能力的解析表征,“Reachable Domain forInterception at Hyperbolic Speeds”首次将提出了可达区域的概念,用于洲际弹道导弹的拦截问题,研究了施加单脉冲后具有双曲线速度的拦截器的可达区域。“ReachableDomain for Spacecraft with a Single Impulse”将可达区域用于航天器的脉冲机动问题,研究了航天器单脉冲后的可达区域问题,考虑了三种情况:1)施加机动的位置固定,脉冲方向任意;2)脉冲方向固定,施加机动的位置任意;3)施加机动的位置和脉冲方向均为任意。但他们都没有研究相对可达区域的概念,这样无法将可达区域的概念运用到相对运动中,另外,对于丢失段的相对可达区域问题,也很少有文献进行研究。
发明内容
本发明的目的在于克服上述现有技术的缺点,提供了一种航天器追踪目标丢失段的相对可达区域计算方法,该方法能够计算得到航天器追踪目标丢失段的相对可达区域。
为达到上述目的,本发明所述的航天器追踪目标丢失段的相对可达区域计算方法包括以下步骤:
1)建立地球惯性坐标系EXYZ及目标航天器轨道坐标系Txyz,并构建航天器相对运动轨道模型,求解航天器相对运动轨道模型的解;
2)根据航天器相对运动轨道模型,在任一时刻加入目标航天器施加的机动力,以构造新的航天器相对运动轨道模型;
3)求解新的航天器相对运动轨道模型的解析解,然后将求解得到的解析解作为目标航天器丢失段的相对可达区域的方程。
步骤1)的具体操作为:
建立地球惯性坐标系EXYZ及目标航天器轨道坐标系Txyz,其中,地球惯性坐标系EXYZ的原点在地心E,地球惯性坐标系EXYZ的X轴指向春分点,地球惯性坐标系EXYZ的Z轴垂直于地球赤道平面向上,地球惯性坐标系EXYZ的Y轴满足右手法则;
目标航天器轨道坐标系Txyz的原点在目标航天器的质心T,目标航天器轨道坐标系Txyz的z轴沿目标航天器的矢径方向并指向地心,目标航天器轨道坐标系Txyz的y轴垂直于轨道平面,目标航天器轨道坐标系Txyz的x轴沿目标航天器的速度方向,且遵循右手法则;
当追踪航天器与目标航天器沿近圆轨道运行时,其相对运动方程能够用C-W方程表示,则在作用力为零的情况下的自由相对运动方程,即:
常值力下C-W方程存在漂移项的解为:
其中,
目标航天器在t0时刻在太阳阴影区内施加机动力,其中,机动方式为连续推力或脉冲推力。
当目标航天器在t0时刻在太阳阴影区内施加连续推力时,设目标航天器施加的连续推力为F,则将连续推力F分解到轨道坐标系下为:
单位质量下的连续推力在轨道坐标系下的分量为:
设目标航天器只施加x轴方向的连续推力,则初始条件为:
根据式(7),目标航天器只施加x轴方向连续推力的相对可达区域方程为:
设目标航天器只施加y轴方向的连续推力,则初始条件为:
根据式(9),目标航天器只施加y轴方向连续推力的相对可达区域方程为:
设目标航天器只施加z轴方向的连续推力,则初始条件为:
根据式(11),目标航天器只施加z轴方向连续推力的相对可达区域方程为:
设目标航天器施加连续推力的方向任意,则初始条件为:
根据式(13),目标航天器施加任意方向连续推力的相对可达区域方程为:
当目标航天器在t0时刻在太阳阴影区内施加脉冲推力时,设目标航天器施加的脉冲推力为ΔV,将脉冲推力ΔV分解到轨道坐标系下为:
设目标航天器只施加x轴方向的脉冲推力,则初始条件为:
根据式(16),目标航天器只施加x轴方向脉冲推力的相对可达区域方程为:
设目标航天器只施加y轴方向的脉冲推力,则初始条件为:
根据式(18),目标航天器只施加y方向脉冲推力的相对可达区域方程为:
设目标航天器只施加z轴方向的脉冲推力,则初始条件为:
根据式(20),目标航天器只施加z方向脉冲推力的相对可达区域方程为:
设目标航天器施加脉冲推力的方向任意,则初始条件为:
根据式(22),目标航天器施加任意方向脉冲推力的相对可达区域方程为:
本发明具有以下有益效果:
本发明所述的航天器追踪目标丢失段的相对可达区域计算方法在具体操作时,根据航天器相对运动轨迹模型在任一时刻加入目标航天器施加的机动力,从而构建出新的航天器相对运动轨道模型,再求解该新的航天器相对运动轨道模型,以得到目标航天器丢失段的相对可达区域的方程,从而确定出目标航天器丢失段的相对可达区域,在实际应用时,追踪航天器即使在没有观测到目标的情况下,同样能够充分把握目标航天器的位置信息,实现追踪航天器对目标航天器可能出现的相对位置区域的掌握。
附图说明
图1为本发明中目标航天器丢失段的示意图;
图2为本发明中构建的地球惯性坐标系EXYZ示意图;
图3为本发明中构建的目标航天器轨道坐标系Txyz示意图。
具体实施方式
下面结合附图对本发明做进一步详细描述:
如图1所示,由于地球围绕太阳转动的缘故,在地球轨道运行的航天器,受到的太阳光照总会有一段时间被地球遮挡,这种情况下,当目标航天器进入阴影区后,利用相机作为传感器的追踪航天器就会失去对目标航天器的观测,因此在该丢失段目标可能实施的机动,利用相对可达区域的概念对目标航天器与追踪航天器之间相对运动进行刻画,使得追踪航天器对目标航天器的位置信息有全面的把握。
本发明所述的航天器追踪目标丢失段的相对可达区域计算方法包括以下步骤:
1)建立地球惯性坐标系EXYZ及目标航天器轨道坐标系Txyz,并构建航天器相对运动轨道模型,求解航天器相对运动轨道模型的解;
2)根据航天器相对运动轨道模型,在任一时刻加入目标航天器施加的机动力,以构造新的航天器相对运动轨道模型;
3)求解新的航天器相对运动轨道模型的解析解,然后将求解得到的解析解作为目标航天器丢失段的相对可达区域的方程。
步骤1)的具体操作为:
建立地球惯性坐标系EXYZ及目标航天器轨道坐标系Txyz,其中,地球惯性坐标系EXYZ的原点在地心E,地球惯性坐标系EXYZ的X轴指向春分点,地球惯性坐标系EXYZ的Z轴垂直于地球赤道平面向上,地球惯性坐标系EXYZ的Y轴满足右手法则;
目标航天器轨道坐标系Txyz的原点在目标航天器的质心T,目标航天器轨道坐标系Txyz的z轴沿目标航天器的矢径方向并指向地心,目标航天器轨道坐标系Txyz的y轴垂直于轨道平面,目标航天器轨道坐标系Txyz的x轴沿目标航天器的速度方向,且遵循右手法则;
当追踪航天器与目标航天器沿近圆轨道运行时,其相对运动方程能够用C-W方程表示,则在作用力为零的情况下的自由相对运动方程,即:
常值力下C-W方程存在漂移项的解为:
其中,
目标航天器在t0时刻在太阳阴影区内施加机动力,其中,机动方式为连续推力或脉冲推力。
当目标航天器在t0时刻在太阳阴影区内施加连续推力时,设目标航天器施加的连续推力为F,则将连续推力F分解到轨道坐标系下为:
单位质量下的连续推力在轨道坐标系下的分量为:
设目标航天器只施加x轴方向的连续推力,则初始条件为:
根据式(7),目标航天器只施加x轴方向连续推力的相对可达区域方程为:
设目标航天器只施加y轴方向的连续推力,则初始条件为:
根据式(9),目标航天器只施加y轴方向连续推力的相对可达区域方程为:
设目标航天器只施加z轴方向的连续推力,则初始条件为:
根据式(11),目标航天器只施加z轴方向连续推力的相对可达区域方程为:
设目标航天器施加连续推力的方向任意,则初始条件为:
根据式(13),目标航天器施加任意方向连续推力的相对可达区域方程为:
当目标航天器在t0时刻在太阳阴影区内施加脉冲推力时,设目标航天器施加的脉冲推力为ΔV,将脉冲推力ΔV分解到轨道坐标系下为:
设目标航天器只施加x轴方向的脉冲推力,则初始条件为:
根据式(16),目标航天器只施加x轴方向脉冲推力的相对可达区域方程为:
设目标航天器只施加y轴方向的脉冲推力,则初始条件为:
根据式(18),目标航天器只施加y方向脉冲推力的相对可达区域方程为:
设目标航天器只施加z轴方向的脉冲推力,则初始条件为:
根据式(20),目标航天器只施加z方向脉冲推力的相对可达区域方程为:
设目标航天器施加脉冲推力的方向任意,则初始条件为:
根据式(22),目标航天器施加任意方向脉冲推力的相对可达区域方程为:
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。
Claims (3)
1.一种航天器追踪目标丢失段的相对可达区域计算方法,其特征在于,包括以下步骤:
1)建立地球惯性坐标系EXYZ及目标航天器轨道坐标系Txyz,并构建航天器相对运动轨道模型,求解航天器相对运动轨道模型的解;
2)根据航天器相对运动轨道模型,在任一时刻加入目标航天器施加的机动力,以构造新的航天器相对运动轨道模型;
3)求解新的航天器相对运动轨道模型的解析解,然后将求解得到的解析解作为目标航天器丢失段的相对可达区域的方程;
步骤1)的具体操作为:
建立地球惯性坐标系EXYZ及目标航天器轨道坐标系Txyz,其中,地球惯性坐标系EXYZ的原点在地心E,地球惯性坐标系EXYZ的X轴指向春分点,地球惯性坐标系EXYZ的Z轴垂直于地球赤道平面向上,地球惯性坐标系EXYZ的Y轴满足右手法则;
目标航天器轨道坐标系Txyz的原点在目标航天器的质心T,目标航天器轨道坐标系Txyz的z轴沿目标航天器的矢径方向并指向地心,目标航天器轨道坐标系Txyz的y轴垂直于轨道平面,目标航天器轨道坐标系Txyz的x轴沿目标航天器的速度方向,且遵循右手法则;
当追踪航天器与目标航天器沿近圆轨道运行时,其相对运动方程能够用C-W方程表示,得在作用力为零的情况下的自由相对运动方程,即:
常值力下C-W方程存在漂移项的解为:
其中,
mc表示追踪航天器的质量,目标航天器在t0时刻在太阳阴影区内施加机动力,其中,机动方式为连续推力或脉冲推力;
当目标航天器在t0时刻在太阳阴影区内施加连续推力时,设目标航天器施加的连续推力为F,则将连续推力F分解到轨道坐标系下为:
单位质量下的连续推力在轨道坐标系下的分量为:
3.根据权利要求2所述的航天器追踪目标丢失段的相对可达区域计算方法,其特征在于,当目标航天器在t0时刻在太阳阴影区内施加脉冲推力时,设目标航天器施加的脉冲推力为△V,将脉冲推力△V分解到轨道坐标系下为:
设目标航天器只施加x轴方向的脉冲推力,则初始条件为:
根据式(16),目标航天器只施加x轴方向脉冲推力的相对可达区域方程为:
设目标航天器只施加y轴方向的脉冲推力,则初始条件为:
根据式(18),目标航天器只施加y方向脉冲推力的相对可达区域方程为:
设目标航天器只施加z轴方向的脉冲推力,则初始条件为:
根据式(20),目标航天器只施加z方向脉冲推力的相对可达区域方程为:
设目标航天器施加脉冲推力的方向任意,则初始条件为:
根据式(22),目标航天器施加任意方向脉冲推力的相对可达区域方程为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810837425.6A CN109159923B (zh) | 2018-07-26 | 2018-07-26 | 一种航天器追踪目标丢失段的相对可达区域计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810837425.6A CN109159923B (zh) | 2018-07-26 | 2018-07-26 | 一种航天器追踪目标丢失段的相对可达区域计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109159923A CN109159923A (zh) | 2019-01-08 |
CN109159923B true CN109159923B (zh) | 2022-05-13 |
Family
ID=64898327
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810837425.6A Expired - Fee Related CN109159923B (zh) | 2018-07-26 | 2018-07-26 | 一种航天器追踪目标丢失段的相对可达区域计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109159923B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111324137B (zh) * | 2020-03-05 | 2022-08-05 | 中国工程物理研究院总体工程研究所 | 一种基于常值连续推力的区域悬停轨道控制方法 |
CN113401364B (zh) * | 2021-06-07 | 2022-06-14 | 中国人民解放军国防科技大学 | 基于日食效应的空间攻防方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0628899A1 (en) * | 1990-04-16 | 1994-12-14 | Space Systems / Loral Inc. | Method and apparatus for inclined orbit attitude control for momentum bias spacecraft |
EP0992429A2 (de) * | 1998-10-06 | 2000-04-12 | SEL Verteidigungssysteme GmbH | Satellitengestütztes Verteidigungssystem und Verfahren zur satellitengestützten Verteidigung |
CN103455707A (zh) * | 2013-07-22 | 2013-12-18 | 西北工业大学 | 基于凸优化技术的有限推力航天器自主交会轨迹规划方法 |
CN103950555A (zh) * | 2014-04-22 | 2014-07-30 | 北京控制工程研究所 | 超近距离的高精度相对位置保持控制方法 |
CN107529376B (zh) * | 2013-08-01 | 2015-12-30 | 上海新跃仪表厂 | 多模融合的微小卫星非合作目标相对导航的方法 |
CN106970530A (zh) * | 2017-04-28 | 2017-07-21 | 西北工业大学 | 空间非合作目标自主视线交会的无模型预设性能控制方法 |
CN107168345A (zh) * | 2017-04-07 | 2017-09-15 | 西北工业大学 | 一种针对空间翻滚目标上任意点的交会方法 |
-
2018
- 2018-07-26 CN CN201810837425.6A patent/CN109159923B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0628899A1 (en) * | 1990-04-16 | 1994-12-14 | Space Systems / Loral Inc. | Method and apparatus for inclined orbit attitude control for momentum bias spacecraft |
EP0992429A2 (de) * | 1998-10-06 | 2000-04-12 | SEL Verteidigungssysteme GmbH | Satellitengestütztes Verteidigungssystem und Verfahren zur satellitengestützten Verteidigung |
CN103455707A (zh) * | 2013-07-22 | 2013-12-18 | 西北工业大学 | 基于凸优化技术的有限推力航天器自主交会轨迹规划方法 |
CN107529376B (zh) * | 2013-08-01 | 2015-12-30 | 上海新跃仪表厂 | 多模融合的微小卫星非合作目标相对导航的方法 |
CN103950555A (zh) * | 2014-04-22 | 2014-07-30 | 北京控制工程研究所 | 超近距离的高精度相对位置保持控制方法 |
CN107168345A (zh) * | 2017-04-07 | 2017-09-15 | 西北工业大学 | 一种针对空间翻滚目标上任意点的交会方法 |
CN106970530A (zh) * | 2017-04-28 | 2017-07-21 | 西北工业大学 | 空间非合作目标自主视线交会的无模型预设性能控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109159923A (zh) | 2019-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Felicetti et al. | Design of robotic manipulators for orbit removal of spent launchers’ stages | |
Prakash et al. | Mars Science Laboratory entry, descent, and landing system overview | |
CN107270933B (zh) | 一种基于多星协同的空间碎片运动状态联合确定方法 | |
Dong et al. | Satellite proximate interception vector guidance based on differential games | |
Liang et al. | Angular velocity tracking for satellite rendezvous and docking | |
CN104252574A (zh) | 一种基于空间系绳抓捕系统的非合作目标质量辨识方法 | |
CN112607065B (zh) | 一种基于电推进系统的高精度相位控制方法 | |
CN109159923B (zh) | 一种航天器追踪目标丢失段的相对可达区域计算方法 | |
CN109190158B (zh) | 一种考虑非合作目标禁飞区约束的最优轨道设计方法 | |
CN109269504B (zh) | 一种具有末端约束的姿态机动路径规划方法 | |
CN112198885B (zh) | 一种满足机动平台自主降落需求的无人机控制方法 | |
CN106931967A (zh) | 一种助推‑滑翔式临近空间飞行器的捷联惯性导航方法 | |
Xu et al. | Modeling and planning of a space robot for capturing tumbling target by approaching the dynamic closest point | |
Zhao et al. | Image-based control for rendezvous and synchronization with a tumbling space debris | |
Wang et al. | Formation flying for electric sails in displaced orbits. Part I: geometrical analysis | |
CN110567462B (zh) | 一种近似自旋非合作航天器三轴转动惯量比的辨识方法 | |
CN103863578A (zh) | 火星着陆器喷气推力器和控制力矩陀螺复合控制系统 | |
CN115535304B (zh) | 一种针对多个编队卫星周期性重访的轨道设计与控制方法 | |
Tao et al. | Research on the evasive strategy of missile based on the theory of differential game | |
CN110466803A (zh) | 基于等倾角姿态控制的自旋稳定卫星姿态预测方法 | |
CN111596678B (zh) | 一种航天器逼近空间翻滚目标的相对位置自适应控制方法 | |
JP2015168315A (ja) | 誘導装置、誘導装置搭載宇宙機 | |
Zubov et al. | Analytical synthesis of a modal controller by output vector for attitude control of a descent module during its descent in the Earth’s atmosphere | |
Low | Nearly circular transfer trajectories for descending satellites | |
CN112518136A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220513 |
|
CF01 | Termination of patent right due to non-payment of annual fee |