CN106507769B - 轨道机动受限的空间交会地面导引段轨道控制方法 - Google Patents
轨道机动受限的空间交会地面导引段轨道控制方法Info
- Publication number
- CN106507769B CN106507769B CN200910121326.9A CN200910121326A CN106507769B CN 106507769 B CN106507769 B CN 106507769B CN 200910121326 A CN200910121326 A CN 200910121326A CN 106507769 B CN106507769 B CN 106507769B
- Authority
- CN
- China
- Prior art keywords
- delta
- theta
- spacecraft
- orbit
- calculating
- 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 54
- 238000004088 simulation Methods 0.000 claims description 21
- 238000005259 measurement Methods 0.000 claims description 10
- 238000012546 transfer Methods 0.000 claims description 10
- 206010048669 Terminal state Diseases 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000001174 ascending effect Effects 0.000 claims description 2
- 230000004323 axial length Effects 0.000 claims 1
- 238000003032 molecular docking Methods 0.000 abstract description 7
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008570 general process Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Abstract
本发明公开了一种轨道机动受限的空间交会地面导引段轨道控制方法,目的是在轨道机动受限的情况下实现空间交会对接地面导引段高精度轨道控制。技术方案是先确定追踪航天器轨道机动方案为五次轨道机动方案,确定需要计算的轨道控制参数M1~M5;追踪航天器入轨后计算轨道控制参数M1,并上传到追踪航天器,由轨道机动系统执行第一次轨道机动;在第一、二、三次轨道机动完成后分别计算轨道控制参数M2、M3、M4,并上传到追踪航天器,由轨道机动系统执行第二、三、四次轨道机动;最后根据轨道控制参数M5(计算M4时得到)实施第五次轨道控制。本发明轨道机动次数少,对姿态控制系统要求低,在轨道机动受限的情况下提高了地面导引段轨道控制精度,且计算速度快,实时性好。
Description
技术领域
本发明涉及空间交会轨道机动控制方法,尤其是轨道机动受限的空间交会地面导引段轨道控制方法。
背景技术
空间交会对接过程一般可以分为地面导引、近程导引、平移靠拢和对接等阶段。在地面导引段,追踪航天器在地面测控的支持下进行数次轨道机动到达目标航天器后下方几十千米的瞄准点。地面导引段轨道控制方法是交会对接轨道设计与控制的关键问题。
当前地面导引段轨道控制方法主要有两类:一类是特殊点变轨轨道控制,另一类是综合变轨轨道控制。
特殊点变轨轨道控制方法主要在美国交会对接工程实践中得到应用和发展,由5次以上在远地点、近地点、升交点等轨道特殊点执行的机动完成,每次机动的方向被限制沿轨道切向或法向。该方法的缺点是轨道机动次数较多、终端控制精度有限。
综合变轨轨道控制方法主要在俄罗斯联盟/进步飞船交会对接工程实践中得到应用和发展,由4~5次不限定方向的轨道机动构成。该方法终端控制精度较高,但由于没有限制轨道机动方向,对航天器姿态机动要求较高、控制相对复杂。
如何在轨道机动受限(即轨道机动次数有限、轨道机动方向是特定方向)情况下实现高精度地面导引段轨道控制一直是困绕本领域技术人员的重要问题。本发明提出一个由不超过5次沿特定方向轨道机动组成的地面导引段轨道控制方法,目前尚无公开文献涉及此类方法。论述中提到的轨道仿真和坐标转换是本领域人员公知的知识,不进行具体阐述。
发明内容
本发明要解决的技术问题是在轨道机动受限的情况下,实现空间交会对接地面导引段高精度轨道控制。具体技术方案如下:
第一步,确定轨道机动方案
采用经典轨道根数形式E=(a,e,i,Ω,w,v)T表述航天器绝对状态,其中a为长半轴,e为偏心率,i为轨道倾角,Ω为升交点赤经,w为近拱点角,v为真近点角。
设追踪航天器入轨时刻为t0,圈次为N0,目标航天器和追踪航天器的初始状态分别为Etar(t0)、Echa(t0),地面导引段的终端时刻为tf,圈次为Nf,终端要求追踪航天器进入目标航天器后下方几十千米处的瞄准点(由用户给定),与目标航天器的相对位置速度一定。
目标航天器不进行轨道机动,依据航天器轨道力学基本原理和地面导引段轨道控制任务需求,经过大量分析与仿真计算,确定追踪航天器轨道机动方案为五次轨道机动方案:
第N1圈,在远地点(对应纬度幅角)沿切向机动,冲量大小为Δvt1;
第N2圈,在纬度幅角处沿法向机动,冲量大小Δvz2;
第N3圈,在近地点(对应纬度幅角)沿切向机动,冲量大小为Δvt3;
第N4圈,在远地点附近纬度幅角处沿切向机动,冲量大小为Δvt4;
第N5圈,在纬度幅角处沿切向机动,冲量大小为Δvt5。
其中,N1,…,N5是各次轨道机动的圈数,具体数值由实际飞行任务测控条件确定,满足N1<N2<N3<N4<N5、N5-N4=1、Nf-N5=1。其中第N5圈机动目的是修正前面四次轨道机动误差。
根据五次轨道机动方案,确定需要计算的轨道控制参数M为以Mi表示第i次变轨前需要计算的轨道控制参数,1≤i≤5,具体形式为
第二步,追踪航天器入轨后,根据实测轨道数据,采用第一迭代方法计算轨道控制参数M1,并上传到追踪航天器,由轨道机动系统执行第一次轨道机动。
步骤2.1采用第一迭代方法计算轨道控制参数M1
2.1.1根据实测的t01时刻追踪航天器、目标航天器初始状态Echa(t01)和Etar(t01),计算参考航天器初始状态Eref(t01),方法如下:
其中,μ为地心引力常数,为追踪航天器初始位置到瞄准点的角距离,为瞄准点纬度幅角,N01为t01时刻追踪航天器轨道飞行圈数。
令参考航天器轨道半径r0=aref(t01),角速度定义X=(δr,r0δθ,δz,δvr,δvt,δvz)T为航天器在参考航天器轨道坐标系中相对状态矢量,轨道坐标系定义如下:x沿径向,由地心指向航天器质心,z沿轨道面法向,y由物理学右手法则确定,δr、δθ、δz、δvr、δvt和δvz分别为径向位置、纬度幅角、法向位置、径向速度、切向速度和法向速度。
2.1.2将Echa(t01)、Etar(t01)和Eref(t01)从t01时刻分别轨道仿真到tf时刻获得Echa(tf)1、Etar(tf)1和Eref(tf)1,分别计算tf时刻瞄准点和追踪航天器在参考航天器轨道坐标系中的状态和
2.1.3根据计算轨道控制参数具体方法如下:
其中,Δθ1、Δθ3、Δθ4分别为N1、N3、N4圈变轨点到终端的角度距离
Δθ4通过求解如公式6所示的一元非线性方程确定
公式6中,
Δvt5目的是对前面四次机动的误差进行修正,在进行第一次轨道控制参数计算时,由于轨道控制尚未引入误差,因此Δvt5=0,
2.1.4根据轨道控制参数轨道仿真获得追踪航天器终端状态
2.1.5根据与Eref(tf)1,计算终端时刻参考航天器轨道坐标系中追踪航天器的状态计算脱靶量
2.1.6如果脱靶量小于迭代收敛标准δ(δ根据实际飞行任务控制精度要求确定),即或立,退出迭代;如果 k=k+1,转步骤2.1.3。
步骤2.2,实施第一次轨道控制。
根据步骤2.1获得的Δvt1,追踪航天器在N1圈远地点沿轨道切向进行第一次轨道机动。
第三步,第一次轨道机动完成后,根据实测轨道数据采用迭代方法计算轨道控制参数M2,并上传到追踪航天器,由轨道机动系统执行第二次轨道机动。
步骤3.1采用第二迭代方法计算精确轨道控制参数M2
3.1.1根据实测的t02时刻追踪航天器、目标航天器初始状态Echa(t02)和Etar(t02),计算参考航天器初始状态Eref(t02),方法如下:
其中,N02为t02时刻追踪航天器轨道飞行圈数。计算r0=aref(t02),
3.1.2将Echa(t02)、Etar(t02)和Eref(t02)从t02时刻分别轨道仿真到tf时刻获得Echa(tf)2、Etar(tf)2和Eref(tf)2,分别计算tf时刻瞄准点和追踪航天器在参考航天器轨道坐标系中的状态和
3.1.3根据计算轨道控制参数具体方法如下:
其中,Δθ3由公式5确定,Δθ5由公式9确定
Δθ4由求解公式10所示的一元方程获得
公式10中,
3.1.4根据轨道控制参数轨道仿真获得追踪航天器终端状态
3.1.5根据与Eref(tf)2,计算终端时刻参考航天器轨道坐标系中追踪航天器状态计算脱靶量
3.1.6如果成立,退出迭代;如果 k=k+1,转步骤3.1.3。
步骤3.2,实施第二次轨道控制。
根据步骤3.1获得的Δvz2,追踪航天器在N2圈纬度幅角处执行第二次轨道机动。
第四步,第二次轨道机动完成后,根据实测轨道数据采用迭代方法计算轨道控制参数M3,并上传到追踪航天器,由轨道机动系统执行第三次轨道机动。
步骤4.1采用第三迭代方法计算精确轨道控制参数M3
4.1.1根据实测的t03时刻追踪航天器、目标航天器初始状态Echa(t03)和Etar(t03),计算参考航天器初始状态Eref(t03),方法如下:
其中,N03为t03时刻追踪航天器轨道飞行圈数。计算r0=aref(t03),
4.1.2将Echa(t03)、Etar(t03)和Eref(t03)从t03时刻分别轨道仿真到tf时刻获得Echa(tf)3、Etar(tf)3和Eref(tf)3,分别计算tf时刻瞄准点和追踪航天器在参考航天器轨道坐标系中的状态和
4.1.3根据计算轨道控制参数具体方法如下:
其中,Δθ3、Δθ4、Δθ5分别由公式5、10和9确定。
4.1.4根据轨道控制参数轨道仿真获得追踪航天器终端状态
4.1.5根据和Eref(tf)3,计算终端时刻参考航天器轨道坐标系中追踪航天器的状态计算脱靶量
4.1.6如果成立,退出迭代;如果 k=k+1,转步骤4.1.3。
步骤4.2,实施第三次轨道控制。
根据步骤4.1获得的Δvt3,追踪航天器在N3圈近地点处执行第三次轨道机动。
第五步,第三次轨道机动完成后,根据实测轨道数据采用迭代方法计算轨道控制参数M4,并上传到追踪航天器,由轨道机动系统执行第四次轨道机动。
步骤5.1采用第四迭代方法计算精确轨道控制参数M4
5.1.1根据实测t04时刻追踪航天器、目标航天器初始状态Echa(t04)和Etar(t04),计算参考航天器初始状态Eref(t04),具体方法如下:
其中,N04为t04时刻追踪航天器轨道飞行圈数。计算r0=aref(t03),
5.1.2将Echa(t04)、Etar(t04)和Eref(t04)从t04时刻分别轨道仿真到tf时刻获得Echa(tf)4、Etar(tf)4和Eref(tf)4,分别计算tf时刻瞄准点和追踪航天器在参考航天器轨道坐标系中的状态和
5.1.3根据计算轨道控制参数具体方法如下:
其中,Δθ4、Δθ5由使用拟牛顿法求解公式15所示的二维非线性方程组确定
5.1.4根据轨道控制参数轨道仿真获得追踪航天器终端状态
5.1.5根据及Eref(tf)4,计算终端时刻参考航天器轨道坐标系中追踪航天器的状态并计算脱靶量
5.1.6如果成立,退出迭代;如果 k=k+1,转步骤5.1.3。
步骤5.2,实施第四次轨道控制。
根据步骤5.1获得的Δvt4,追踪航天器在N4圈纬度幅角处执行第四次轨道机动。
第六步,根据第五步计算获得的轨道控制参数M5实施第五次轨道控制。
第五次轨道控制参数不再重新计算,根据第五步计算获得的和Δvt5,追踪航天器在N5圈纬度幅角处执行第五次轨道机动。
采用本发明可以达到以下技术效果:
1.采用本发明可以在轨道机动受限的情况下提高地面导引段轨道控制精度。由于将第四和第五次轨道机动位置(纬度幅角)增加为轨道控制迭代变量,虽然机动方向受限,但仍然可以直接瞄准终端相对位置速度;由于在每次轨道机动之前基于实测轨道数据重新计算轨道控制参数,后面的机动可以修正前面机动的误差;最后两次机动距离地面导引终端不超过两圈,误差传播时间较短,从而有效提高了地面导引段轨道控制精度。本发明应用于我国载人航天二期工程神舟8号交会对接轨道方案初样设计中。
2.本发明采用的算法简单,计算速度快,实时性好,可以应用于实际轨道飞行控制。提出的基于参考轨道的迭代方法简单,收敛性得到了有效保证。在进行神舟8号初样方案数千次打靶试验中,本方法的收敛性在99.9%以上,在CPU 2.8GHz计算机上进行一次飞行控制模拟,总耗时不超过2分钟,效率高。
3.本发明轨道机动次数少,对姿态控制系统要求低,没有附加成本。相比美国航天飞机的调相变轨方案(6次以上轨道机动),本发明减少了1次以上轨道机动,本发明的轨道机动次数同俄罗斯联盟/进步号飞船相当,但本发明的轨道机动方向主要沿切向或法向方向,不同于俄罗斯的轨道机动方向是空间任意方向,可以有效减轻姿态控制系统负担。
总之,本发明具有精度高、执行简单和计算速度快等的优点,具有广阔的应用前景。
附图说明
图1是本发明总体流程图。
图2是本发明计算精确轨道控制参数的迭代方法流程图,其中下标i=1,2,3,4时分别对应第一、二、三、四迭代方法。
具体实施方式
图1是本发明总体流程图,具体流程是:
第一步,设计轨道机动方案。
第二步,追踪航天器入轨后,根据实测轨道数据采用第一迭代方法计算第一~五次机动的轨道控制参数,并上传到追踪航天器,由轨道机动系统执行第一次轨道机动。
第三步,第一次轨道机动完成后,根据实测轨道数据采用第二迭代方法计算第二~五次机动的轨道控制参数,并上传到追踪航天器,由轨道机动系统执行第二次轨道机动。
第四步,第二次轨道机动完成后,根据实测轨道数据采用第三迭代方法计算第三~五次机动的轨道控制参数,并上传到追踪航天器,由轨道机动系统执行第三次轨道机动。
第五步,第三次轨道机动完成后,根据实测轨道数据采用第四迭代方法计算第四~五次机动的轨道控制参数,并上传到追踪航天器,由轨道机动系统执行第四次轨道机动。
第六步,第四次轨道机动完成后,根据第五步计算获得的第五次轨道机动控制参数,追踪航天器轨道机动系统执行第五次轨道机动。
图2是本发明计算轨道控制参数的迭代方法流程图。第一、第二、第三、第四迭代方法的总体流程是一致的,只是初始时刻追踪航天器和目标航天器状态不同,轨道控制参数及计算方法不同。它们的流程都包括以下六步:
1)根据实测的t0i时刻追踪航天器、目标航天器初始状态Echa(t0i)和Etar(t0i),计算参考航天器初始状态Eref(t0i);
2)轨道仿真到tf时刻获得Echa(tf)i、Etar(tf)i和Eref(tf)i,计算tf时刻瞄准点和追踪航天器在参考航天器轨道坐标系中的状态和
3)根据计算轨道控制参数
4)根据轨道仿真获得追踪航天器终端状态
5)根据及Eref(tf)i,计算终端时刻参考航天器轨道坐标系中追踪航天器的状态计算脱靶量
6)判断脱靶量是否小于允许值,如果成立,退出迭代;如果 k=k+1,转步骤3)。
Claims (1)
1.一种轨道机动受限的空间交会地面导引段轨道控制方法,其特征在于包括以下步骤:
第一步,确定轨道机动方案:
采用经典轨道根数形式E=(a,e,i,Ω,w,v)T表述航天器绝对状态,其中a为半长轴长,e为偏心率,i为轨道倾角,Ω为升交点赤经,w为近拱点角,v为真近点角;设追踪航天器入轨时刻为t0,圈次为N0,目标航天器和追踪航天器的初始状态分别为Etar(t0)、Echa(t0),地面导引段的终端时刻为tf,圈次为Nf;目标航天器不进行轨道机动,依据航天器轨道力学基本原理和地面导引段轨道控制任务需求,确定追踪航天器轨道机动方案为五次轨道机动方案:
第N1圈,在远地点沿切向机动,冲量大小为Δvt1;
第N2圈,在纬度幅角处沿法向机动,冲量大小Δvz2;
第N3圈,在近地点沿切向机动,冲量大小为Δvt3;
第N4圈,在远地点附近纬度幅角处沿切向机动,冲量大小为Δvt4;
第N5圈,在纬度幅角处沿切向机动,冲量大小为Δvt5;
N1,…,N5是各次轨道机动的圈数,具体数值由实际飞行任务测控条件确定,满足N1<N2<N3<N4<N5、N5-N4=1、Nf-N5=1;
根据五次轨道机动方案,确定需要计算的轨道控制参数M为以Mi表示第i次变轨前需要计算的轨道控制参数,1≤i≤5,具体形式为
第二步,追踪航天器入轨后,根据实测轨道数据,采用第一迭代方法计算轨道控制参数M1,并上传到追踪航天器,由轨道机动系统执行第一次轨道机动:
步骤2.1采用第一迭代方法计算轨道控制参数M1:
2.1.1根据实测的t01时刻追踪航天器、目标航天器初始状态Echa(t01)和Etar(t01),计算参考航天器初始状态Eref(t01):
其中,μ为地心引力常数,为追踪航天器初始位置到瞄准点的角距离,为瞄准点纬度幅角,N01为t01时刻追踪航天器轨道飞行圈数;
2.1.2将Echa(t01)、Etar(t01)和Eref(t01)从t01时刻分别轨道仿真到tf时刻获得的Echa(tf)1、Etar(tf)1和Eref(tf)1,分别计算tf时刻瞄准点和追踪航天器在参考航天器轨道坐标系中的状态和
其中,X=(δr,r0δθ,δz,δvr,δvt,δvz)T为航天器在参考航天器轨道坐标系中相对状态矢量,轨道坐标系定义如下:x沿径向,由地心指向航天器质心,z沿轨道面法向,y由右手法则确定,δr、δθ、δz、δvr、δvt和δvz分别为径向位置、纬度幅角、法向位置、径向速度、切向速度和法向速度;航天器轨道半径r0=aref(t01),角速度
2.1.3根据计算轨道控制参数
其中,Δθ1、Δθ3、Δθ4分别为N1、N3、N4圈变轨点到终端的角度距离,
Δθ4通过求解如公式6所示的一元非线性方程确定
公式6中,
Δvt5=0,
2.1.4根据轨道控制参数轨道仿真获得追踪航天器终端状态
2.1.5根据与Eref(tf)1,计算终端时刻参考航天器轨道坐标系中追踪航天器的状态计算脱靶量
2.1.6如果脱靶量小于迭代收敛标准δ,δ根据实际飞行任务控制精度要求确定,即成立,退出迭代;如果 k=k+1,转步骤2.1.3;
步骤2.2实施第一次轨道控制:根据步骤2.1获得的Δvt1,追踪航天器在N1圈远地点沿轨道切向进行第一次轨道机动;
第三步,第一次轨道机动完成后,根据实测轨道数据采用第二迭代方法计算轨道控制参数M2,并上传到追踪航天器,由轨道机动系统执行第二次轨道机动:
步骤3.1采用第二迭代方法计算精确轨道控制参数M2:
3.1.1根据实测的t02时刻追踪航天器、目标航天器初始状态Echa(t02)和Etar(t02),计算参考航天器初始状态Eref(t02):
其中,N02为t02时刻追踪航天器轨道飞行圈数;计算r0=aref(t02),
3.1.2将Echa(t02)、Etar(t02)和Eref(t02)从t02时刻分别轨道仿真到tf时刻获得Echa(tf)2、Etar(tf)2和Eref(tf)2,分别计算tf时刻瞄准点和追踪航天器在参考航天器轨道坐标系中的状态和
3.1.3根据计算轨道控制参数
其中,Δθ3由式5确定,Δθ5由式9确定
Δθ4由求解公式10所示的一元方程获得
公式10中,
3.1.4根据轨道控制参数轨道仿真获得追踪航天器终端状态
3.1.5根据与Eref(tf)2计算终端时刻参考航天器轨道坐标系中追踪航天器状态计算脱靶量
3.1.6如果成立,退出迭代;如果 k=k+1,转步骤3.1.3;
步骤3.2,实施第二次轨道控制:根据步骤3.1获得的Δvz2,追踪航天器在N2圈纬度幅角处执行第二次轨道机动;
第四步,第二次轨道机动完成后,根据实测轨道数据采用第三迭代方法计算轨道控制参数M3,并上传到追踪航天器,由轨道机动系统执行第三次轨道机动:
步骤4.1采用第三迭代方法计算精确轨道控制参数M3:
4.1.1根据t03时刻实测的追踪航天器、目标航天器初始状态Echa(t03)和Etar(t03),计算参考航天器初始状态Eref(t03):
其中,N03为t03时刻追踪航天器轨道飞行圈数,计算r0=aref(t03),
4.1.2将Echa(t03)、Etar(t03)和Eref(t03)从t03时刻分别轨道仿真到tf时刻获得Echa(tf)3、Etar(tf)3和Eref(tf)3,分别计算tf时刻瞄准点和追踪航天器在参考航天器轨道坐标系中的状态和
4.1.3根据计算轨道控制参数
其中,Δθ3、Δθ4、Δθ5分别由式5、10和9确定;
4.1.4根据轨道控制参数轨道仿真获得追踪航天器终端状态
4.1.5根据和Eref(tf)3计算终端时刻参考航天器轨道坐标系中追踪航天器的状态计算脱靶量
4.1.6如果成立,退出迭代;如果 k=k+1,转步骤4.1.3;
步骤4.2,实施第三次轨道控制:根据步骤4.1获得的Δvt3,追踪航天器在N3圈近地点处执行第三次轨道机动;
第五步,第三次轨道机动完成后,根据实测轨道数据采用迭代方法计算轨道控制参数M4,并上传到追踪航天器,由轨道机动系统执行第四次轨道机动:
步骤5.1采用第四迭代方法计算精确轨道控制参数M4:
5.1.1根据实测的t04时刻追踪航天器、目标航天器初始状态Echa(t04)和Etar(t04),计算参考航天器初始状态Eref(t04):
其中,N04为t04时刻追踪航天器轨道飞行圈数;计算r0=aref(t03),
5.1.2将Echa(t04)、Etar(t04)和Eref(t04)从t04时刻分别轨道仿真到tf时刻获得Echa(tf)4、Etar(tf)4和Eref(tf)4,分别计算tf时刻瞄准点和追踪航天器在参考航天器轨道坐标系中的状态和
5.1.3根据计算轨道控制参数
其中,Δθ4、Δθ5由使用拟牛顿法求解公式15所示的二维非线性方程组确定,
5.1.4根据轨道控制参数轨道仿真获得追踪航天器终端状态
5.1.5根据及计算终端时刻参考航天器轨道坐标系中追踪航天器的状态并计算脱靶量
5.1.6如果成立,退出迭代;如果 k=k+1,转步骤5.1.3;
步骤5.2,实施第四次轨道控制:根据步骤5.1获得的Δvt4,追踪航天器在N4圈纬度幅角处执行第四次轨道机动;
第六步,根据第五步计算获得的轨道控制参数M5实施第五次轨道控制:根据第五步计算获得的和Δvt5,追踪航天器在N5圈纬度幅角处执行第五次轨道机动。
Publications (1)
Publication Number | Publication Date |
---|---|
CN106507769B true CN106507769B (zh) | 2013-01-23 |
Family
ID=
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108490966A (zh) * | 2018-01-31 | 2018-09-04 | 中国人民解放军国防科技大学 | 基于微分代数的静止轨道摄动相对轨迹高阶制导方法 |
CN109484674A (zh) * | 2018-10-12 | 2019-03-19 | 湖北航天技术研究院总体设计所 | 一种基于目标轨道参数的实时轨道机动控制方法 |
CN110765504A (zh) * | 2019-10-29 | 2020-02-07 | 北京空间技术研制试验中心 | 航天器环月轨道交会对接的轨道设计方法 |
CN111504330A (zh) * | 2020-05-07 | 2020-08-07 | 中国人民解放军63768部队 | 一种基于实测数据的低轨航天器准实时轨道机动检测方法 |
CN113148232A (zh) * | 2021-02-24 | 2021-07-23 | 中国人民解放军国防科技大学 | 机动大小固定的单脉冲轨道瞄准方法和装置 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108490966A (zh) * | 2018-01-31 | 2018-09-04 | 中国人民解放军国防科技大学 | 基于微分代数的静止轨道摄动相对轨迹高阶制导方法 |
CN108490966B (zh) * | 2018-01-31 | 2021-02-05 | 中国人民解放军国防科技大学 | 基于微分代数的静止轨道摄动相对轨迹高阶制导方法 |
CN109484674A (zh) * | 2018-10-12 | 2019-03-19 | 湖北航天技术研究院总体设计所 | 一种基于目标轨道参数的实时轨道机动控制方法 |
CN109484674B (zh) * | 2018-10-12 | 2020-12-25 | 湖北航天技术研究院总体设计所 | 一种基于目标轨道参数的实时轨道机动控制方法 |
CN110765504A (zh) * | 2019-10-29 | 2020-02-07 | 北京空间技术研制试验中心 | 航天器环月轨道交会对接的轨道设计方法 |
CN110765504B (zh) * | 2019-10-29 | 2022-01-18 | 北京空间技术研制试验中心 | 航天器环月轨道交会对接的轨道设计方法 |
CN111504330A (zh) * | 2020-05-07 | 2020-08-07 | 中国人民解放军63768部队 | 一种基于实测数据的低轨航天器准实时轨道机动检测方法 |
CN113148232A (zh) * | 2021-02-24 | 2021-07-23 | 中国人民解放军国防科技大学 | 机动大小固定的单脉冲轨道瞄准方法和装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104898642B (zh) | 一种用于航天器姿态控制算法的集成测试仿真系统 | |
CN107618678B (zh) | 卫星姿态角度偏差下的姿控信息联合估计方法 | |
CN103662091B (zh) | 一种基于相对导航的高精度安全着陆制导方法 | |
CN103955223B (zh) | 一种深空探测软着陆过程的姿轨耦合控制方法 | |
CN105865459B (zh) | 一种考虑视线角约束的小天体接近段制导方法 | |
CN105222780B (zh) | 一种基于Stirling插值多项式逼近的椭球集员滤波方法 | |
CN104960674A (zh) | 一种运动目标的指向跟踪控制方法 | |
CN104076819B (zh) | 一种圆参考轨道下卫星有界伴飞的边界控制方法 | |
CN101794154A (zh) | 一种编队卫星相对轨道与姿态解耦控制方法 | |
CN109839116B (zh) | 一种速度增量最小的交会接近方法 | |
CN102116628A (zh) | 一种着陆或附着深空天体探测器的高精度导航方法 | |
CN103019247B (zh) | 一种火星探测器无陀螺自主空间姿态机动控制方法 | |
CN104015191B (zh) | 基于基座卫星角速度的空间机械臂工具坐标下的运动补偿方法 | |
CN104281150A (zh) | 一种姿态机动的轨迹规划方法 | |
Yoshimitsu et al. | Hayabusa-final autonomous descent and landing based on target marker tracking | |
CN107300700B (zh) | 敏捷合成孔径雷达卫星聚束模式姿态机动需求计算方法 | |
CN103591956B (zh) | 一种基于可观测性分析的深空探测器自主导航方法 | |
CN105930305B (zh) | 一种三脉冲交会接近制导方法 | |
CN112198885B (zh) | 一种满足机动平台自主降落需求的无人机控制方法 | |
CN111605737A (zh) | 一种航天器三超控制多级协同规划与敏捷机动方法 | |
CN102878997A (zh) | 一种大偏心率轨道的星上快速高精度外推方法 | |
CN115343960B (zh) | 一种地月系统中航天器光照阴影规避控制方法 | |
CN105334731B (zh) | 一种远地点点火高精度解析轨道自主预报方法 | |
CN108082538B (zh) | 一种考虑始末约束的多体系统低能量捕获轨道方法 | |
CN106507769B (zh) | 轨道机动受限的空间交会地面导引段轨道控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
GR03 | Grant of secret patent right | ||
DC01 | Secret patent status has been lifted | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130123 Termination date: 20180508 |