CN111102982B - 高轨目标的抵近方法 - Google Patents
高轨目标的抵近方法 Download PDFInfo
- Publication number
- CN111102982B CN111102982B CN201911268099.2A CN201911268099A CN111102982B CN 111102982 B CN111102982 B CN 111102982B CN 201911268099 A CN201911268099 A CN 201911268099A CN 111102982 B CN111102982 B CN 111102982B
- Authority
- CN
- China
- Prior art keywords
- target
- plane
- orbit
- approaching
- track
- 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
Images
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/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
-
- 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/20—Instruments for performing navigational calculations
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Navigation (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供了一种高轨目标的抵近方法,本发明通过分析不同高轨目标的运动特征,把目标分为面内目标和面外目标两大类情况;对于面内目标,在约定抵近距离范围后以抵近末端状态光照情况较好作为约束条件;对于面外目标,在约定抵近距离范围后以共同过赤道作为约束。针对上述两种情况,首先分析面内目标的抵近窗口。每个轨道周期内有一个最优成像窗口,通过计算变轨速度增量需求并优选变轨时刻,使得抵近到目标附近时恰好满足最优成像条件。再分析面外目标的抵近窗口,每个轨道周期内有两个最近距离窗口,通过计算变轨速度增量需求并优选变轨时刻,使得抵近到目标附近时恰好满足共同过赤道并有相对较优的成像条件。
Description
技术领域
本发明涉及民用航天空间在轨服务任务技术领域,具体地,涉及对高轨目标或碎片的抵近机动方法。
背景技术
高轨在轨服务任务包括对目标抵近观察、对目标实施维修。为了高效的抵近到碎片目标附近,服务星一般运行在漂飞轨道,目标运行在标称GEO轨道。需要到更近距离以获得高分辨图像,感知目标的状态信息,本技术可实现对目标近距离抵近甚至伴飞在目标附近,从而实施信息获取并进而对目标维修、加注等在轨服务任务。
申请号为201910472758.8的发明专利公开了高轨卫星SAR动目标检测方法,考虑到了在高轨模式下的两个主要误差,一是传统机载SAR的“走-停”假设不成立,二是由于地球自转和轨道轨迹带来的非线性变化。本发明将DPCA方法应用于高轨卫星SAR系统中,能够在抑制杂波的同时保留淹没于杂波区的慢速运动目标,并同时考虑到高轨模型下的非线性斜距变化,构造频域匹配函数对目标的距离走动和多普勒走动进行补偿,相比于传统的仅仅最多考虑到二阶项有更好的运动目标检测与成像效果。但是无法针对不同类型的目标实现不同的抵进策略。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种高轨目标的抵近方法。
根据本发明提供的一种高轨目标的抵近方法,包括如下步骤:
目标类型确定步骤:确定目标类型,是面内目标还是面外目标;
速度增量大小确认步骤:针对抵近高度,确定霍曼变轨两次实施的速度增量大小;
求解步骤:对面内目标或者对面外目标,建立满足约束的时间方程,找到可以进行变轨的时刻为解,并优选光照条件好的解;
验证步骤:对解进行验证。
优选地,目标类型确定步骤包括:确定目标为面内目标还是面外目标,通过如下公式计算我星和目标的轨道面夹角:
cos(di)=cos(i0).cos(i1)+sin(i0).sin(i1).cos(Ω0-Ω1)
其中,di为轨道面夹角,i0为目标轨道倾角,i1为我星轨道倾角,Ω0为目标升交点赤经,Ω1为我星升交点赤经;
进行计算后判断,当di≤0.1°时,目标类型为面内目标;当di>0.1°时,目标类型为面外目标。
优选地,所述速度增量大小确认步骤包括:根据抵近轨道的机动高度差计算所需的速度增量V:
V=sqrt(μ/a);
a=Re+h;
其中:Re为地球半径,取6378.137km;h为轨道高度,μ为地球常数398600.5。
优选地,所述求解步骤包括:
其中:△t为目标回到轨道面交点处所需时间;T表示一个轨道周期;ω1和ω2为轨道调整前后我星的漂移角速度;θ表示地心角;t1表示实施霍曼变轨的时刻,t1的可取范围为[0,θ/ω1]。
优选地,所述验证步骤包括:建立场景,进行求解验证。
与现有技术相比,本发明具有如下的有益效果:
1、本发明提供了一种高轨目标或碎片的抵近方法,通过对高轨抵近任务的需求进行分析,计算对不同类型的目标采取不同的抵近策略。
2、本发明首先按照轨道面夹角对目标进行分类,然后按照目标类型分别计算可选的变轨窗口。按照轨道高度差计算出变轨所需的速度增量后,分两次完成霍曼变轨。
3、本发明具有重要的应用价值,对高轨卫星目标观测和维修、在轨加注等任务提供了有价值的抵近策略和方法。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为不同地方时太阳和目标的关系示意图;
图2为面内目标抵近方法示意图;
图3为面外目标抵近方法示意图;
图4为抵近面内目标JCSAT3A示意图;
图5为抵近过程相对距离和JCSAT3A星下点地方时变化图;
图6为抵近过程JCSAT3A视线和太阳矢量夹角变化图;
图7为抵近面外目标DSP22示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
如图1至图7所示,根据本发明提供的一种高轨目标的抵近方法,包括如下步骤:
步骤A:确定目标类型,是面内目标还是面外目标;确定目标为面内目标还是面外目标,通过如下公式计算我和目标的轨道面夹角,
cos(di)=cos(i0).cos(i1)+sin(i0).sin(i1).cos(Ω0-Ω1)
其中,di为轨道面夹角,i0为目标轨道倾角,i1为我星轨道倾角,Ω0为目标升交点赤经,Ω1为我星升交点赤经。进行计算后判断,当di≤0.1°时,可认为是面内目标。当di>0.1°时,当作面外目标处理。
步骤B:针对抵近高度,确定霍曼变轨两次实施的速度增量大小;根据抵近轨道的机动高度差计算所需的速度增量。
a=Re+h;
V=sqrt(μ/a);
其中Re为地球半径,取6378.137km;h为轨道高度,μ为地球常数398600.5。
步骤C:对面内目标,建立满足约束的时间方程,找到可以进行变轨的时刻为解;对面外目标,建立满足约束的时间方程,找到可以进行变轨的时刻为解,并优选光照条件好的解;建立变轨机动总耗时与目标到约定地方时的时间的等时方程。不同地方时太阳与目标的关系如图1所示,对于轨道面内目标,选择合理的抵近目标时的地方时,以下漂为例,可选择21:00~03:00左右。如图2所示,卫星通过调整轨道高度,选择24点左右的地方时经过目标正下方(B点),从而完成对目标的优良光照条件的成像观测。卫星通过调整轨道高度,卫星选择24点左右的地方时经过目标正下方(B′点),从而完成对目标的优良光照条件的成像观测。
假定T0时卫星在A′点,目标当前地方时和期望地方时相差Δt。经过时间t1,卫星运行到B′点,这时候实施霍曼变轨,将卫星轨道高度变到距离目标轨道的约定高度差,两次霍曼变轨共耗时一天。再经过一段时间t2,在期望地方时卫星经过目标正下方。轨道调整前后卫星的漂移角速度分别为ω1和ω2,则卫星漂过地心角θ的时间应该等于目标n个轨道周期加上Δt,即有
ttarget=△t+nT
式中T表示一个轨道周期即24小时,建立上述两式的等式
随着n的取值不同,起始时间选择也不同,t1的可取范围为[0,θ/ω1],通过选取合适的n,决定变轨实施的时刻,从而完成抵近变轨过程。
求解面外目标变轨脉冲实施时刻计算方法时,建立变轨机动总耗时与目标到约定地方时的时间的等时方程。
如图3所示,卫星通过调整轨道高度,使得当目标经过B点(升/降交点)时,卫星同时经过目标正下方(C点),从而完成对目标的抵近。假定T0时卫星在A′点,目标在A点,目标回到轨道面交点处所需时间为Δt。经过时间t1,卫星运行到B′点,这时候实施霍曼变轨,将卫星轨道高度变到距离静止轨道约20公里左右,两次霍曼变轨共耗时一天。再经过一段时间t2,卫星和目标同时过交点。轨道调整前后卫星的漂移角速度分别为ω1和ω2,则卫星漂过地心角θ的时间应该等于目标从A点第n+1次(n=0,1,2…)经过B点的时间,即有
式中即24小时,建立上述两式的等式
随着n的取值不同,起始时间选择也不同,t1的可取范围为[0,θ/ω1],通过选取合适的n,决定变轨实施的时刻,从而完成抵近变轨过程。
步骤E:建立场景,对解进行验证。
下面给出具体实施例:根据本发明提供的高轨目标的抵近方法,包括如下步骤:
步骤A:建立我星和目标的相对运动场景,分别建立目标JCSAT-3A和目标DSP22;
步骤B:计算我星和目标的轨道面夹角关系,与JCSAT-3A夹角0.02°,为面内目标;与DSP22夹角5.8°,为面外目标;
步骤C:根据轨道高度差计算霍曼变轨所需的速度增量,达到最近距离约20多公里则高度变化量分别为75km和65km,对应速度增量分别为1.37m/s和1.19m/s;
步骤D:对目标,分别根据时间约束方程求解变轨时刻;
步骤E:分析抵近终端状态,对JCSAT-3A,最近距离27km,地方时为1:00,太阳光照角15°,观测光照条件良好;对DSP22最近距离22km,与目标同时过赤道;
表1为JCSAT3A卫星参数;表2为抵近JCSAT3A时变轨参数;表3为DSP22卫星参数;表4为抵近DSP22时变轨参数。
表1
时间 | 2019.2.20 03:00:00 | 倾角 | 0.021° |
定点位置 | 127.910° | 升交点 | 220.855° |
半长轴 | 42165.382km | 近地点幅角 | 356.466° |
偏心率 | 0.000171 | 真近点角 | 105.123° |
表2
面外目标参数、对其抵近变轨策略如下所示。
表3
表4
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (4)
1.一种高轨目标的抵近方法,其特征在于,包括如下步骤:
目标类型确定步骤:确定目标类型,是面内目标还是面外目标;
速度增量大小确认步骤:针对抵近高度,确定霍曼变轨两次实施的速度增量大小;
求解步骤:对面内目标或者对面外目标,建立满足约束的时间方程,找到可以进行变轨的时刻为解;
验证步骤:对解进行验证;
目标类型确定步骤包括:确定目标为面内目标还是面外目标,通过如下公式计算卫 星和目标的轨道面夹角:
cos(di)=cos(i0).cos(i1)+sin(i0).sin(i1).cos(Ω0-Ω1)
其中,di为轨道面夹角,i0为目标轨道倾角,i1为卫 星轨道倾角,Ω0为目标升交点赤经,Ω1为卫 星升交点赤经;
进行计算后判断,当di≤0.1°时,目标类型为面内目标;当di>0.1°时,目标类型为面外目标。
2.根据权利要求1所述的高轨目标的抵近方法,其特征在于,所述速度增量大小确认步骤包括:根据抵近轨道的机动高度差计算所需的速度增量V:
V=sqrt(μ/a);
a=Re+h;
其中:Re为地球半径,取6378.137km;h为轨道高度,μ为地球常数398600.5。
4.根据权利要求1所述的高轨目标的抵近方法,其特征在于,所述验证步骤包括:建立场景,进行求解验证。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911268099.2A CN111102982B (zh) | 2019-12-11 | 2019-12-11 | 高轨目标的抵近方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911268099.2A CN111102982B (zh) | 2019-12-11 | 2019-12-11 | 高轨目标的抵近方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111102982A CN111102982A (zh) | 2020-05-05 |
CN111102982B true CN111102982B (zh) | 2021-09-24 |
Family
ID=70422437
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911268099.2A Active CN111102982B (zh) | 2019-12-11 | 2019-12-11 | 高轨目标的抵近方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111102982B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113636106B (zh) * | 2021-09-15 | 2023-06-30 | 上海卫星工程研究所 | 连续小推力高轨目标变轨抵近方法及系统 |
CN114537714B (zh) * | 2022-02-28 | 2022-11-25 | 中国人民解放军63921部队 | 一种高轨卫星变轨调控方法及系统 |
CN117009606B (zh) * | 2023-09-22 | 2023-12-29 | 中科星图测控技术股份有限公司 | 一种考虑约束的观测星轨道机动进入点选择方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH02200600A (ja) * | 1988-12-09 | 1990-08-08 | General Electric Co <Ge> | 角度測定を使用した宇宙船の接近/分離 |
CN108279703A (zh) * | 2018-01-26 | 2018-07-13 | 河南工程学院 | 一种用于非合作机动目标拦截的轨道控制方法 |
CN109774974A (zh) * | 2019-01-31 | 2019-05-21 | 上海微小卫星工程中心 | 一种用于空间碎片接近的轨道设计方法 |
CN109839116A (zh) * | 2019-01-31 | 2019-06-04 | 上海微小卫星工程中心 | 一种速度增量最小的交会接近方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003084813A2 (en) * | 1999-03-11 | 2003-10-16 | Constellation Services International | Method of using dwell times in intermediate orbits to optimise orbital transfers and method and apparatus for satellite repair |
CN102538759B (zh) * | 2011-12-29 | 2015-01-07 | 中国科学院紫金山天文台 | 近地和中高轨空间目标实时全自动捕获方法 |
CN103728980B (zh) * | 2014-01-08 | 2016-08-31 | 哈尔滨工业大学 | 航天器相对轨道的控制方法 |
CN109592079A (zh) * | 2018-12-03 | 2019-04-09 | 哈尔滨工业大学 | 一种限定时间的航天器共面交会变轨策略确定方法 |
-
2019
- 2019-12-11 CN CN201911268099.2A patent/CN111102982B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH02200600A (ja) * | 1988-12-09 | 1990-08-08 | General Electric Co <Ge> | 角度測定を使用した宇宙船の接近/分離 |
CN108279703A (zh) * | 2018-01-26 | 2018-07-13 | 河南工程学院 | 一种用于非合作机动目标拦截的轨道控制方法 |
CN109774974A (zh) * | 2019-01-31 | 2019-05-21 | 上海微小卫星工程中心 | 一种用于空间碎片接近的轨道设计方法 |
CN109839116A (zh) * | 2019-01-31 | 2019-06-04 | 上海微小卫星工程中心 | 一种速度增量最小的交会接近方法 |
Non-Patent Citations (1)
Title |
---|
接近和跟踪非合作机动目标的非线性最优控制;高登巍等;《宇航学报》;20130630;第34卷(第6期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111102982A (zh) | 2020-05-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111102982B (zh) | 高轨目标的抵近方法 | |
CN106595674B (zh) | 基于星敏感器和星间链路的heo卫星编队飞行自主导航方法 | |
CN106595673B (zh) | 面对地球静止轨道目标操作的空间多机器人自主导航方法 | |
CN106643741B (zh) | 一种卫星相对小行星视觉自主导航方法 | |
CN106564622B (zh) | 一种基于视线测量的同步轨道卫星相对倾角远距离修正方法 | |
CN112591146B (zh) | 一种高轨目标分钟级快速遍历的观测方法及系统 | |
CN106767768A (zh) | 一种双星编队的自主导航方法 | |
CN102878995A (zh) | 一种静止轨道卫星自主导航方法 | |
CN109269504B (zh) | 一种具有末端约束的姿态机动路径规划方法 | |
CN106679653A (zh) | 一种基于星敏感器和星间链路的heo卫星群相对测量方法 | |
CN109269510A (zh) | 基于星敏感器和星间链路的heo卫星编队飞行自主导航方法 | |
CN103900577A (zh) | 一种面向编队飞行的相对导航测速及组合导航方法 | |
CN110053788A (zh) | 一种考虑复杂摄动的星座长期保持控制频次估计方法 | |
Chen et al. | Orbital evasive target tracking and sensor management | |
Baranov et al. | Ballistic aspects of large-size space debris flyby at low Earth near-circular orbits | |
CN113343442A (zh) | 一种求解固定时间有限燃料多脉冲转移轨道的方法及系统 | |
Ender et al. | Bistatic exploration using spaceborne and airborne SAR sensors: A close collaboration between FGAN, ZESS, and FOMAAS | |
CN109269508A (zh) | 一种卫星相对小行星视觉自主导航方法 | |
Cui et al. | Real-time navigation for Mars final approach using X-ray pulsars | |
CN106643742B (zh) | 一种卫星自主连续观测小行星的方法 | |
CN117649077A (zh) | 一种地球同步轨道目标成像卫星重点详查任务规划方法 | |
CN102607563A (zh) | 利用背景天文信息对于航天器进行相对导航的系统 | |
CN106855408A (zh) | 一种面向geo卫星在轨服务的空间多机器人相对观测方法 | |
CN106871890A (zh) | 一种用于卫星编队的星间相对观测方法 | |
Kominato et al. | Optical hybrid navigation and station keeping around Itokawa |
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 |