CN107323689A - 在轨卫星对地勘察的轨道机动方法 - Google Patents

在轨卫星对地勘察的轨道机动方法 Download PDF

Info

Publication number
CN107323689A
CN107323689A CN201710260635.9A CN201710260635A CN107323689A CN 107323689 A CN107323689 A CN 107323689A CN 201710260635 A CN201710260635 A CN 201710260635A CN 107323689 A CN107323689 A CN 107323689A
Authority
CN
China
Prior art keywords
satellite
orbit
ground
substar
load
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.)
Granted
Application number
CN201710260635.9A
Other languages
English (en)
Other versions
CN107323689B (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.)
63789 Army Of Pla
Original Assignee
63789 Army Of Pla
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 63789 Army Of Pla filed Critical 63789 Army Of Pla
Priority to CN201710260635.9A priority Critical patent/CN107323689B/zh
Publication of CN107323689A publication Critical patent/CN107323689A/zh
Application granted granted Critical
Publication of CN107323689B publication Critical patent/CN107323689B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/242Orbits and trajectories

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Navigation (AREA)

Abstract

本发明提供了一种在轨卫星对地勘察的轨道机动方法,利用卫星地固系星历计算出卫星载荷波束在地面目标方向的覆盖能力,通过比较地面目标与卫星经过地面目标的同纬度圈的经度、轨道机动时刻计算得到地面目标与卫星载荷地面投影相遇的时间差,由相遇时间差及运行圈次计算出卫星满足勘察要求所要具备的轨道周期,进而得到卫星轨道机动半长轴所需调整量,筛选出在卫星控制窗口期最小控制量及对应的地面目标覆盖时间,此控制量即为约束条件下最节省燃料的轨道机动控制量。本发明实现了短时内以最节省燃料的方式进行轨道机动在特定时间段内对地面目标的勘察。

Description

在轨卫星对地勘察的轨道机动方法
技术领域
本发明属于在轨航天器测控管理领域,涉及一种实现在轨卫星对地面目标勘察的轨道机动计算方法。
背景技术
对地勘察卫星一般为回归轨道,设计的回归轨道可使有效载荷对地面固定区域进行反复勘察。通过多颗卫星组网以多种载荷可对地面目标进行短时勘察。由于卫星组网前勘察的区域不同,组网后对新目标勘察会出现某些时段勘察覆盖间隙过大。为了缩短对新目标的最大重访间隔,需要短时间内完成对组网内某些卫星的轨道机动,使得这些卫星在指定的时间段内覆盖勘察目标。由于需要短时间内完成对地面目标覆盖勘察,另外卫星携带燃料有限,在约束条件下节省燃料也是轨道机动的优选方式。现有技术只是针对卫星标准设计轨道的固定勘察目标以地球摄动力长时间进行轨道机动,不具备对特殊需求的其他目标进行轨道机动,同时也不能短时完成在指定时段内对新目标的覆盖勘察,在多星临时组网中无法大幅缩短最大勘察覆盖间隙。
发明内容
为了克服现有技术的不足,本发明提供一种在轨卫星对地勘察的轨道机动方法,通过计算卫星载荷与地面目标同纬度经度差,设计了卫星载荷对地投影与地面特定目标在指定的短时段内相遇,筛选给出符合约束条件最节省燃料的轨道机动控制量。
本发明解决其技术问题所采用的技术方案包括以下步骤:
已知勘查卫星在地固系的位置星历为S(xl,yl,zl),N为正整数;星载载荷在最大侧摆角φ下波束中心在地面的投影E,地球质心点o,地面目标R的地固系位置坐标为(xr,yr,zr),大地坐标为(lonr,latr,hr),lon表示大地经度,lat表示大地纬度,h表示大地高;设卫星载荷波束为圆锥型,锥角即载荷波束宽度为2η;当卫星在滚动方向侧摆时,卫星载荷波束与卫星地心矢量的夹角在(φ-η,φ+η)之间,则载荷波束在地球的投影为长半轴a1、短半轴b的半椭圆和长半轴a2、短半轴b的半椭圆组成的椭圆图形,a1>a2,椭圆图形的中心即星载载荷波束中心在地面的投影E;
根据轨道控制要求卫星联合勘查要求勘查目标的时间从控制时刻Tc开始后Td时间段内完成;由S(xl,yl,zl)通过拉格朗日插值计算出(Tc,Tc+Td)时段内每圈次卫星在处的卫星位置S(xsj,ysj,zsj)及星下点G(xgj,ygj,zgj),latΔ=0.05, j=1,…,10;
通过S(xsj,ysj,zsj)计算卫星与地心之间的距离
长半轴a1顶点与星下点G之间的弧长
长半轴a2顶点与星下点G之间的弧长
投影E与星下点G之间的弧长
短轴b顶点与星下点G之间的弧长
长半轴a1顶点与星下点G之间的弦长
长半轴a2顶点与星下点G之间的弦长
投影E与星下点G之间的弦长其中,Rej为卫星位于位置 S(xsj,ysj,zsj)时刻Tsj的星下点G与地心o之间的距离,
其中,ae为地球赤道半径,f为地球扁率;
已知载荷波束中心点E的地固系坐标为(xej,yej,zej)、星下点G及地面目标点R之间的夹角为∠GER,
如果∠GERj在投影椭圆的第二、三象限,则星载载荷在地面目标方向的覆盖能力如果∠GERj在投影椭圆的第一、四象限,则星载载荷在地面目标方向的覆盖能力
根据(FEk-1-FEk)(FEk-FEk+1)<0计算得到FEk,即为卫星此圈经过交点周期latr的最短弧长,FEk在地面的覆盖经度FElonk=FEk/loneq,其中,勘察点latr纬度圈的每度对应的弧长大地子午圈偏心率
卫星每圈次经过latr的交点周期为Tsm,m=1,2,…,n,n为Td时段内卫星运行的圈次,
其中,轨道周期GM为地心引力常数;J2=1082.636×10-6; a为轨道半长轴;e为卫星轨道偏心率;i为卫星轨道倾角;u*=ω+fs,ω为卫星近地点幅角,fs为卫星真近点角;
与轨道机动时刻Tc的时间差ΔTscm=Tsm-Tc
若卫星在Tsm时刻勘查目标,从轨道机动时刻Tc开始,卫星轨道调整运行的圈次
卫星载荷波束与地面目标的经度差Δlonm=lonsm-lonr-FElonm,其中,lonsm为Tsm时刻星下点G的大地坐标经度;
卫星轨道周期增量其中,ωe为地球自转角速度;
则在Tsm时刻卫星覆盖目标,需要的半长轴卫星轨道机动半长轴控制量Δam=a-am';
计算出(Tc,Tc+Td)时段内卫星轨道机动控制量Δam集合,从中筛选出最小控制量Δa,则此控制量是在勘查要求的启控时刻Tc及Td时间段内完成勘查所需的最小控制量;Δa对应的Tsm时刻即为轨道调整后卫星载荷勘查覆盖地面目标的时刻。
本发明的有益效果是:通过给出轨道机动方向及轨道半长轴控制量实现了对现有技术无法观测的地面特定目标的覆盖观测;依据卫星载荷投影每圈次与地面目标的相遇覆盖,从同纬度控制量中筛选出最小控制量,计算的轨道机动控制量在勘查时间段内最节省燃料。本发明实现了在指定的短时间段内,通过轨道机动对地面目标的勘察。
附图说明
图1卫星载荷侧摆覆盖区域示意图。
图2是卫星载荷投影示意图。
具体实施方式
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
本发明利用卫星地固系星历计算出卫星载荷波束在地面目标方向的覆盖能力,通过比较地面目标与卫星经过地面目标的同纬度圈的经度、轨道机动时刻计算得到地面目标与卫星载荷地面投影相遇的时间差,由相遇时间差及运行圈次计算出卫星满足勘察要求所要具备的轨道周期,进而得到卫星轨道机动半长轴所需调整量,筛选出在卫星控制窗口期最小控制量及对应的地面目标覆盖时间,此控制量即为约束条件下最节省燃料的轨道机动控制量,从而实现了短时内以最节省燃料的方式进行轨道机动在特定时间段内对地面目标的勘察。
1.计算每圈最短弧长
已知勘查卫星在地固系的位置星历为(其中:N为正整数;坐标原点为地球质心o,x的方向为地球质心指向格林尼治子午圈,z指向北极的国际习用原点,x,y,z成右手系)和星载载荷在最大侧摆角φ下波束中心在地面的投影E,地球质心点o,地面目标的R的地固系位置坐标为(xr,yr,zr),则大地坐标为(lonr,latr,hr)(其中:lon为大地经度,lat为大地纬度,h为大地高)。设卫星载荷波束为简单圆锥型,锥角即载荷波束宽度为2η。当卫星在滚动方向侧摆时,卫星载荷波束与卫星与地心矢量的夹角在(φ-η,φ+η)之间,则载荷波束在地球的投影为长半轴a1、短半轴b和长半轴a2(a1>a2)、短半轴b组成的椭圆图形,椭圆中心即星载载荷波束中心在地面的投影E。
根据轨道控制要求若卫星轨道控制时刻为Tc,卫星联合勘查要求勘查目标的时间从控制时刻Tc开始后Td时间段内完成。由(其中:N为正整数)通过拉格朗日插值计算出(Tc,Tc+Td)时段内每圈次卫星在(其中latΔ=0.05, j为正整数,j=1,…,10)处的卫星位置S(xsj,ysj,zsj)及星下点G(xgj,ygj,zgj)。
通过S(xsj,ysj,zsj)可得卫星与地心之间的距离so:
星载载荷在地面投影椭圆长半轴顶点a1与星下点G之间的弧长为:
星载载荷在地面投影椭圆长半轴顶点a2与星下点G之间的弧长为:
星载载荷波束中心在地面的投影E与星下点G之间的弧长为:
星载载荷在地面投影椭圆短轴顶点b与星下点G之间的弧长为:
计算星载载荷在地面的投影椭圆参数。星载载荷在地面投影椭圆长半轴顶点a1与星下点G之间的弦长为:
星载载荷在地面投影椭圆长半轴顶点a2与星下点G之间的弦长为:
星载载荷波束中心在地面的投影E与星下点G之间的弦长为:
其中,Rej为Tsj时刻卫星星下点G与地心o之间的距离:
其中,ae为地球赤道半径,ae=6378137.0m
f为地球扁率,
已知载荷波束中心点E的地固系坐标为(xej,yej,zej)(其中:坐标原点为地球质心o, x的方向为地球质心指向格林尼治子午圈,z指向北极的国际习用原点,x,y,z成右手系)、星下点G及地面目标点R之间的夹角∠GER:
星载载荷在地面目标方向的覆盖能力FE:
如果∠GERj在投影椭圆的第二、三象限,则
如果∠GERj在投影椭圆的第一、四象限,则
求得FEj(其中:j为正整数,j=1,…,10)序列,则由下式计算得到即为卫星此圈经过交点周期latr的最短弧长。
(FEk-1-FEk)(FEk-FEk+1)<0(k=2,3,…,9)
FEk在地面的覆盖经度为:
FElonk=FEk/loneq
其中,loneq为勘察点latr纬度圈的每度对应的弧长;
eE为大地子午圈偏心率,
通过计算可以得到(Tc,Tc+Td)时段内卫星每圈次卫星经过latr载荷的覆盖能力序列
2.轨道调整量计算与筛选
卫星每圈次经过latr的交点周期Tsm(m=1,2,…,n)(n为Td时段内卫星运行的圈次,),与轨道周期T有以下关系:
其中,GM为地心引力常数,GM=398600.5×109m3/s2
J2=1082.636×10-6
a为轨道半长轴;
e为卫星轨道偏心率;
i为卫星轨道倾角;
u*=ω+fs,ω为卫星近地点幅角,fs为卫星真近点角。
与轨道机动时刻Tc的时间差ΔTsc为:
ΔTscm=Tsm-Tc
若卫星在Tsm时刻勘查目标,从轨道机动时刻Tc开始,卫星轨道调整运行的圈次Δn:
卫星载荷波束与地面目标的经度差Δlon:
Δlonm=lonsm-lonr-FElonm
其中,lonsm为Tsm时刻星下点G的大地坐标经度。
卫星轨道周期增量为ΔT:
其中,ωe为地球自转角速度,ωe=7292115.1467×10-11rad/s
则在Tsm时刻卫星覆盖目标,需要的半长轴a':
其中,GM为地心引力常数,GM=398600.5×109m3/s2
卫星轨道机动半长轴控制量为Δam
Δam=a-am'
其中,a为卫星轨道半长轴。
计算出(Tc,Tc+Td)时段内卫星轨道机动控制量Δam集合,则从Δam中筛选出最小控制量Δa,则此控制量是在勘查要求的启控时刻Tc及Td时间段内完成勘查所需的最小控制量。
Δa=min{Δam}
Δa对应的Tsm时刻即为轨道调整后卫星载荷勘查覆盖地面目标的时刻。
已知A、B星的最大侧摆角为35°,轨道参数如表1所示,下面以A、B两颗卫星临时组网给出本发明的具体实施示例。
表1
地面目标R的大地坐标为(0.0,120.0,0)。
则两星对目标R的最大不可见时间为:2017-03-02 18:04:8.013到2017-03-0418:28:41.847,最大间隔约为2904分钟。
对A星轨道进行机动,指定轨道机动时刻为:2017-03-01 21:00:0.0。
计算最短弧长(单位:千米):
6618.28184
8770.1592
8648.75912
6435.16606
3919.13527
1352.6118
1224.08141
3791.54478
6312.32926
8561.32743
8844.86976
6740.90815
4238.78626
1675.55997
1900.87193
3470.88105
计算半长轴修正量:
1)计算的轨道机动半长轴调整序列为(单位:千米):
598.2143786
-1246.451076
326.992727
-1394.639875
86.77948588
-1527.685794
-71.80289151
1417.500139
-267.6155261
1134.866448
-444.5540333
881.6288482
-605.2777892
653.1737185
-751.9574975
445.8576701
2)轨道机动调整量为(单位:千米):
最节省燃料的轨道机动半长轴调整量为:-71.80289151。控制时刻2017-03-0121:00:0.0A星半长轴为6863.228985,调整后的半长轴为6791.42609349。
经过轨道调整后,则两星对目标R的最大不可见时间为:2017-03-02 18:04:8.013到2017-03-03 07:04:43.833,最大间隔约为780分钟。

Claims (1)

1.一种在轨卫星对地勘察的轨道机动方法,其特征在于包括下述步骤:
已知勘查卫星在地固系的位置星历为S(xl,yl,zl),N为正整数;星载载荷在最大侧摆角φ下波束中心在地面的投影E,地球质心点o,地面目标R的地固系位置坐标为(xr,yr,zr),大地坐标为(lonr,latr,hr),lon表示大地经度,lat表示大地纬度,h表示大地高;设卫星载荷波束为圆锥型,锥角即载荷波束宽度为2η;当卫星在滚动方向侧摆时,卫星载荷波束与卫星地心矢量的夹角在(φ-η,φ+η)之间,则载荷波束在地球的投影为长半轴a1、短半轴b的半椭圆和长半轴a2、短半轴b的半椭圆组成的椭圆图形,a1>a2,椭圆图形的中心即星载载荷波束中心在地面的投影E;
根据轨道控制要求卫星联合勘查要求勘查目标的时间从控制时刻Tc开始后Td时间段内完成;由S(xl,yl,zl)通过拉格朗日插值计算出(Tc,Tc+Td)时段内每圈次卫星在处的卫星位置S(xsj,ysj,zsj)及星下点G(xgj,ygj,zgj),latΔ=0.05,j=1,…,10;
通过S(xsj,ysj,zsj)计算卫星与地心之间的距离
长半轴a1顶点与星下点G之间的弧长
长半轴a2顶点与星下点G之间的弧长
投影E与星下点G之间的弧长
短轴b顶点与星下点G之间的弧长
长半轴a1顶点与星下点G之间的弦长
长半轴a2顶点与星下点G之间的弦长
投影E与星下点G之间的弦长其中,Rej为卫星位于位置 S(xsj,ysj,zsj)时刻Tsj的星下点G与地心o之间的距离,
其中,ae为地球赤道半径,f为地球扁率;
已知载荷波束中心点E的地固系坐标为(xej,yej,zej)、星下点G及地面目标点R之间的夹角为∠GER,
如果∠GERj在投影椭圆的第二、三象限,则星载载荷在地面目标方向的覆盖能力如果∠GERj在投影椭圆的第一、四象限,则星载载荷在地面目标方向的覆盖能力
根据(FEk-1-FEk)(FEk-FEk+1)<0计算得到FEk,即为卫星此圈经过交点周期latr的最短弧长,FEk在地面的覆盖经度FElonk=FEk/loneq,其中,勘察点latr纬度圈的每度对应的弧长大地子午圈偏心率
卫星每圈次经过latr的交点周期为Tsm,m=1,2,…,n,n为Td时段内卫星运行的圈次,
其中,轨道周期GM为地心引力常数;J2=1082.636×10-6;a为轨道半长轴;e为卫星轨道偏心率;i为卫星轨道倾角;u*=ω+fs,ω为卫星近地点幅角,fs为卫星真近点角;
与轨道机动时刻Tc的时间差ΔTscm=Tsm-Tc
若卫星在Tsm时刻勘查目标,从轨道机动时刻Tc开始,卫星轨道调整运行的圈次
卫星载荷波束与地面目标的经度差Δlonm=lonsm-lonr-FElonm,其中,lonsm为Tsm时刻星下点G的大地坐标经度;
卫星轨道周期增量其中,ωe为地球自转角速度;
则在Tsm时刻卫星覆盖目标,需要的半长轴卫星轨道机动半长轴控制量Δam=a-am';
计算出(Tc,Tc+Td)时段内卫星轨道机动控制量Δam集合,从中筛选出最小控制量Δa,则此控制量是在勘查要求的启控时刻Tc及Td时间段内完成勘查所需的最小控制量;Δa对应的Tsm时刻即为轨道调整后卫星载荷勘查覆盖地面目标的时刻。
CN201710260635.9A 2017-04-20 2017-04-20 在轨卫星对地勘察的轨道机动方法 Expired - Fee Related CN107323689B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710260635.9A CN107323689B (zh) 2017-04-20 2017-04-20 在轨卫星对地勘察的轨道机动方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710260635.9A CN107323689B (zh) 2017-04-20 2017-04-20 在轨卫星对地勘察的轨道机动方法

Publications (2)

Publication Number Publication Date
CN107323689A true CN107323689A (zh) 2017-11-07
CN107323689B CN107323689B (zh) 2019-11-05

Family

ID=60193473

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710260635.9A Expired - Fee Related CN107323689B (zh) 2017-04-20 2017-04-20 在轨卫星对地勘察的轨道机动方法

Country Status (1)

Country Link
CN (1) CN107323689B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108536990A (zh) * 2018-04-26 2018-09-14 上海微小卫星工程中心 一种计算重访卫星载荷入射角随轨道漂移量变化的方法
CN109684688A (zh) * 2018-12-11 2019-04-26 上海航天控制技术研究所 一种基准载荷与目标载荷交叉定标轨道设计方法
CN109697325A (zh) * 2018-12-27 2019-04-30 中国人民解放军战略支援部队信息工程大学 一种实时卫星对地覆盖时间窗口的确定方法及装置
CN111949922A (zh) * 2020-07-31 2020-11-17 上海卫星工程研究所 适用于对地探测任务多时间窗的星上快速计算的方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1110862B1 (en) * 1999-12-23 2004-03-10 The Boeing Company Combined roll-yaw spacecraft steering method for low earth orbit target trajectory compensation
CN102323571A (zh) * 2011-05-31 2012-01-18 北京航空航天大学 综合总体参数的星载双天线sar干涉定标器布放方法
CN102439393A (zh) * 2009-05-15 2012-05-02 密歇根宇航公司 范围成像激光雷达
CN104348539A (zh) * 2014-07-23 2015-02-11 中国人民解放军海军工程大学 卫星通信点波束覆盖区域计算方法
CN106134458B (zh) * 2013-07-31 2015-05-13 中国人民解放军63956部队 用于目标探测/伪装防护的飞艇载全极化sar系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1110862B1 (en) * 1999-12-23 2004-03-10 The Boeing Company Combined roll-yaw spacecraft steering method for low earth orbit target trajectory compensation
CN102439393A (zh) * 2009-05-15 2012-05-02 密歇根宇航公司 范围成像激光雷达
CN102323571A (zh) * 2011-05-31 2012-01-18 北京航空航天大学 综合总体参数的星载双天线sar干涉定标器布放方法
CN106134458B (zh) * 2013-07-31 2015-05-13 中国人民解放军63956部队 用于目标探测/伪装防护的飞艇载全极化sar系统
CN104348539A (zh) * 2014-07-23 2015-02-11 中国人民解放军海军工程大学 卫星通信点波束覆盖区域计算方法

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108536990A (zh) * 2018-04-26 2018-09-14 上海微小卫星工程中心 一种计算重访卫星载荷入射角随轨道漂移量变化的方法
CN108536990B (zh) * 2018-04-26 2022-07-12 上海微小卫星工程中心 一种计算重访卫星载荷入射角随轨道漂移量变化的方法
CN109684688A (zh) * 2018-12-11 2019-04-26 上海航天控制技术研究所 一种基准载荷与目标载荷交叉定标轨道设计方法
CN109684688B (zh) * 2018-12-11 2023-05-30 上海航天控制技术研究所 一种基准载荷与目标载荷交叉定标轨道设计方法
CN109697325A (zh) * 2018-12-27 2019-04-30 中国人民解放军战略支援部队信息工程大学 一种实时卫星对地覆盖时间窗口的确定方法及装置
CN109697325B (zh) * 2018-12-27 2022-10-25 中国人民解放军战略支援部队信息工程大学 一种实时卫星对地覆盖时间窗口的确定方法及装置
CN111949922A (zh) * 2020-07-31 2020-11-17 上海卫星工程研究所 适用于对地探测任务多时间窗的星上快速计算的方法及系统
CN111949922B (zh) * 2020-07-31 2023-11-28 上海卫星工程研究所 适用于对地探测任务多时间窗的星上快速计算的方法及系统

Also Published As

Publication number Publication date
CN107323689B (zh) 2019-11-05

Similar Documents

Publication Publication Date Title
CN107323689B (zh) 在轨卫星对地勘察的轨道机动方法
US5957409A (en) Ellipitical satellite system which emulates the characteristics of geosychronous satellites
CN106209205B (zh) 一种重点区域按需覆盖的全球通信星座设计方法
US9309010B2 (en) Methods and apparatus for controlling a plurality of satellites using node-synchronous eccentricity control
CN103675832B (zh) 一种快速重访离散目标的低轨卫星轨道设计方法
US6868316B1 (en) Satellite constellation system
CN113632090B (zh) 全球碳盘点卫星的轨道设计系统
US9533774B1 (en) Highly inclined elliptical orbit launch and orbit acquisition techniques
CN112623277B (zh) 一种对空间异面圆轨道目标快速抵达的变轨方法
WO2017147478A1 (en) Stationkeeping techniques for spacecraft in inclined geosynchronous orbit
CN107505948A (zh) 一种用于敏捷卫星机动中沿曲线条带成像的姿态调整方法
CN108613655B (zh) 一种敏捷卫星机动中沿斜条带成像的姿态调整方法
CN112498746B (zh) 一种自主规划卫星沿经线推扫时间及姿态的方法
CN112130590B (zh) 一种基于瞬时惯性系下速度补偿的星载天线对地指向确定方法
US20010051521A1 (en) Elliptical satellite system which emulates the characteristics of geosynchronous satellites
US6611683B1 (en) Elliptical satellite system which emulates the characteristics of geosynchronous satellites
Yang et al. Problem A of 9th China trajectory optimization competition: Results found at NUAA
WO2015160416A9 (en) Communication satellite system
Campagnola et al. Design of lunar gravity assist for the BepiColombo mission to Mercury
CN113014311A (zh) 一种覆盖局部区域的低轨共地面轨迹星座设计方法
JP2527895B2 (ja) 衛星管制方法
CN115936324A (zh) 一种地球同步轨道目标成像卫星巡视任务规划方法
Yahya et al. Delta-V and ground metric investigations for responsive satellite formation flying
Ivashkin Low energy trajectories for the Moon-to-Earth space flight
WO2016113586A1 (en) System of satellites

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20191105

Termination date: 20210420