一种速度增量最小的交会接近方法与流程

本发明涉及航天
技术领域：
，尤其涉及一种机动卫星用最小速度增量交会接近目标卫星的方法。
背景技术：
：本发明中的交会接近是指通过轨道控制使机动卫星与目标卫星位置相同但对速度无约束。无论是同轨道的交会接近还是异面轨道的交会接近，都可以归结成lambert问题，也就是空间两点相对引力的位置矢量确定，当给定时间或给定合适半长轴，即能确定过这两点的转移轨道。lambert问题也常见于导弹/卫星拦截、星际航行、碎片接近等航天活动，是航天领域的研究热点之一。对于lambert问题，可以依据有无时间约束确定两种研究方向。前者如张鹏宇博士的《空间拦截最优轨道设计》，就是在给定时间的情况下研究速度增量最小的转移轨道，但这种分析未从机动卫星和目标卫星的轨道特性出发，速度增量也只能是某种时间约束条件下的最小，换一种时间约束又会得到另一种结果，对于一般卫星来说，携带的燃料有限，时间约束反而没有燃料约束重要，因此首要考虑的约束是燃料消耗最小，也就是后一种研究方向，即无时间约束的情况下能量最优问题。在齐映红博士的《小卫星追踪拦截制导问题研究》，初步探讨了这种情况，但只是针对机动卫星和目标卫星是共轨圆轨道的情况下，得出霍曼转移轨道是速度增量最小的转移轨道，对于机动卫星和目标卫星是非圆轨道以及轨道异面等情况，则无法适用。另外，对于交会接近必须考虑摄动、有限推力等因素，若不考虑，会使机动卫星接近目标卫星的精度难以达标。比如以模型差异为例，对于倾角45度、800公里的圆轨道高度，twobody模型和hpop模型在外推一天情况下，两者误差就可达到1000公里量级。技术实现要素：有鉴于上述技术问题，本发明提供了一种交会接近方法，包括机动卫星和目标卫星轨道异面时，根据所述机动卫星和所述目标卫星的参数，得到停泊轨道和目标轨道的位置和速度；求解所述停泊轨道和所述目标轨道的交线矢量和相关参数；计算所述停泊轨道和所述目标轨道的交线地心距差并确定所述机动卫星和所述目标卫星的交会点；根据摄动情况下的单脉冲轨道设计，确定最优轨控点；通过闭合轨道控制调整相位，使所述机动卫星和所述目标卫星最终在所述交会点交会接近。进一步地，所述相关参数包括停泊轨道纬度幅角、停泊轨道交线幅角、停泊轨道交线幅角差、目标轨道纬度幅角、目标轨道交线幅角、目标轨道交线幅角差。进一步地，所述交线地心距差是停泊轨道交点径向距与目标轨道交点径向距的差异。进一步地，所述交会点是所述目标轨道上对应于较小的交线地心距差的交线点。进一步地，在所述单脉冲轨道设计中，通过平面投影法克服所述停泊轨道和所述目标轨道的法向偏差。进一步地，在所述单脉冲轨道设计中，所述交会点采用相应的虚拟交会点，所述虚拟交会点是事先计算出二体模型与精确模型外推后的偏差矢量，然后将所述偏差矢量叠加到当前交会点后形成的交会点。进一步地，在所述单脉冲轨道设计中，从所述交会点对应时刻回推一轨，通过插值得到360个可能轨控点，每两个所述可能轨控点间隔1度，求解所述交会点和与各所述可能轨控点之间的最小速度增量，从中选择最小值对应的轨控点作为所述最优轨控点。进一步地，所述闭合轨道控制至少分解成两次控制。进一步地，所述闭合轨道控制中的最后一次控制时按照固定时长转移进行速度增量迭代求解，直至所述机动卫星和所述目标卫星满足交会时刻和距离的设定要求。本发明还提供一种交会接近方法，包括机动卫星和目标卫星轨道同面时，根据所述机动卫星和所述目标卫星的参数，得到停泊轨道和目标轨道的位置和速度；求解所述停泊轨道和所述目标轨道的纬度幅角；将所述目标轨道均分成若干个点作为可能交会点；从每个所述可能交会点对应时刻回推一轨，通过插值得到若干个可能轨控点；求解所述可能交会点和与其相应的所述可能轨控点之间的最小速度增量，从中选择最小值对应的轨控点和交会点分别作为最优轨控点和最优交会点；通过闭合轨道控制调整相位，使所述机动卫星和所述目标卫星最终在所述最优交会点交会接近。在本发明的一个实施例中，将所述目标轨道均分成360个点作为所述可能交会点；每个所述可能交会点对应时刻回推一轨，通过插值得到360个可能轨控点。进一步地，所述闭合轨道控制至少分解成两次控制。进一步地，所述闭合轨道控制中的最后一次控制时按照固定时长转移进行速度增量迭代求解，直至所述机动卫星和所述目标卫星满足交会时刻和距离的设定要求。无论机动卫星与目标卫星轨道是共面还是异面，轨道为近圆轨道还是椭圆轨道，本发明的交会接近方法都可以通过对机动卫星和目标卫星的轨道特性分析，得到速度增量最小的最优转移轨道，并在考虑推进能力的情况下兼顾时间约束，提出可行的控制策略实现最终的交会接近。本发明充分考虑了各种摄动力、有限推力及最小轨控间隔等实际实施因素，具有工程可实现性，适用于一般卫星甚至微纳卫星的交会接近任务。以下将结合附图对本发明的构思、具体结构及产生的技术效果作进一步说明，以充分地了解本发明的目的、特征和效果。附图说明图1是本发明的一个实施例中闭合轨道式的多次控制实施流程图；图2中本发明的一个实施例的仿真实验中机动卫星与目标卫星交会接近三维图；图3是本发明的一个实施例的仿真实验中机动卫星与目标卫星距离的变化曲线图。具体实施方式在本发明的实施方式的描述中，需要理解的是，术语“上”、“下”、“前”、“后”、“左”、“右”、“垂直”、“水平”、“顶”、“底”、“内”、“外”、“顺时针”、“逆时针”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对发明的限制。附图为原理图或者概念图，各部分厚度与宽度之间的关系，以及各部分之间的比例关系等等，与其实际值并非完全一致。实施例一本实施例的交会接近方法首先依据测定轨得到机动卫星和目标卫星的位置速度以及大气阻力、迎风面质比、光压系数、光压面质比，通过stk或matlab等工具精确外推，得到停泊轨道位置速度(rp,vp)，目标轨道位置速度(rt,vt)，下标p为停泊轨道，t为目标轨道，当机动卫星与目标卫星轨道异面时，整个交会接近的计算步骤如下：(一)求解交线矢量和相关参数机动卫星轨道面比角动量：hp＝rp×vp(1)设惯性坐标系下的坐标系的z轴矢量为：k＝[001]'停泊轨道节线方程：np＝k×hp(2)停泊轨道节线与停泊轨道位置的叉乘：rnp＝np×rp(3)停泊轨道节线到当前位置的夹角为纬度幅角，具体为：其中，np＝|np|,rp＝|rp|，||符号为对矢量求模，下同。目标卫星轨道面比角动量：ht＝rt×vt(5)机动卫星轨道面与目标卫星的轨道面相交的线为交线矢量，具体如下：b＝hp×ht(6)停泊轨道节线与交线的叉乘：bnp＝np×b(7)停泊轨道节线到交线矢量的夹角为停泊轨道交线幅角，具体为：其中，np＝|np|,b＝|b|停泊轨道交线幅角差为：δup＝up-upb(9)参照停泊轨道计算过程,计算目标轨道交线幅角差δut，并以δut计算圈数，设初始圈数为nt＝0，当δut每次覆盖一次完整的0～360度的范围时，当前圈次加1，即nt＝nt+1。(二)计算交线地心距差和确定交会点选定可能的交会点圈次为第n～m圈，初始圈次n需能满足最小轨控间隔要求，具体见第四部分内容。交线矢量b与轨道面的交点为交线点，对于一个轨道面，只有一个交线矢量及两个交线点，此时对应的交线幅角差为0,π，交线矢量会随着摄动发生变化。设当前轨为第k圈时，k的初值为n圈，通过拉格朗日插值法得到第k圈目标轨道与交线矢量b的交线点：时，目标轨道交线点为时，目标轨道交线点为当目标轨道交线幅角差为0,π对应的时刻，当前停泊轨道的交线幅角差不一定也正好0,π，因此对应目标轨道的第k圈，往回取出停泊轨道数据，插值得到与目标轨道交线幅角差0,π最接近时刻的停泊轨道与交线矢量b的交线点：时，停泊轨道交线点为时，停泊轨道交线点为交线矢量在停泊轨道和目标轨道各有交点，停泊轨道的交点径向距与目标轨道交点径向距的差异为交线地心距差，交线地心距差为：交会点为机动卫星与目标卫星相交的点，通过大量仿真数据知道，对于异面轨道交会接近来说，速度增量最小的转移轨道交会点出现在目标轨道的交线点上，并且是在交线地心距差较小的那个交线点上，设交会点为b，确定原则如下：若则定为第n圈对应交会时刻，交会点交线幅角差δub＝0，目标轨道交会点位置：停泊轨道在交会点方向对应的交线点为：否则，为第n圈对应交会时刻，交会点交线幅角差δub＝π，目标轨道交会点位置：停泊轨道在交会点方向对应的交线点为：在此需要说明的是，未轨控时，停泊轨道在交会点方向对应的交线点时刻与目标轨道交会点时刻一般情况下不在同一时刻。(三)摄动情况下单脉冲轨道设计对于机动卫星和目标卫星的交会接近来说误差自然越小越好，而误差传递又与时间有关，因此希望轨控点与交会点的时间越短越好，由此在单脉冲轨道设计时采用当轨转移。当停泊轨道与目标轨道为异面轨道时，目标轨道按上一步确定第k圈交会点，由于是当轨转移，因此停泊轨道轨控点与目标轨道第k圈交会点的时间少于1轨时间，所以停泊轨道以此点时刻往回推1轨插值得到360个可能轨控点，每两个可能轨控点间隔1度，设此类可能轨控点的下标为pa，插值后得到的停泊轨道数据记为：δupa,rpa,vpa，而交会点的轨道数据记为δub,rtb，为克服轨道偏转带来的法向偏差，用平面投影法得到停泊轨道轨控点δupa,和目标轨道交会点δub,然后对可能轨控点和交会点求解最小速度增量，求解的时候考虑摄动带来的径向偏差，最终从上述仿真数据中统计出速度增量最小的轨控点，需提醒一下的是这部分的计算只是从可能轨控点力找出最优轨控点，然后以脉冲的方式从最优轨控点到交会点上建立转移轨道，但在无额外的手段时，机动卫星与目标卫星即使已经有最优转移轨道，也并不能同时到交会点上。具体计算过程如下：1、平面投影法后的位置和速度由于摄动的关系，轨道面会发生偏转，即在交线点以外的轨道面上的点与交线点的密切根数的倾角或升交点赤经会不一样，若是直接由当前的根数与交会点求解转移轨道，会由于法向的偏差带来额外的速度增量。将交线点为交会点，无论机动点在停泊轨道面那一点上，本质上还是同轨道面转移，因此为避免由于轨道面偏转带来的额外的速度增量，在此采用平面投影法，即将目标轨道和停泊轨道都投影到倾角为0，升交点为0的虚拟平面上，并用交线幅角差计算当前对应的角度，然后算出平面投影后的轨道位置和速度，具体步骤如下：首先计算停泊轨道可能轨控点的轨道根数：径向距离：rpa＝|rpa|(11)比角动量：hpa＝rpa×vpa(12)倾角：设惯性坐标系下的坐标系的z轴矢量为：k＝[001]'停泊轨道节线方程：npa＝k×hpa(14)升交点赤经：偏心率矢量：其中，μ为重力常数，μ＝398600.4418km3/sec2。偏心率：epa＝|epa|(17)近地点幅角：径向速度大小：垂直径向沿着飞行方向的矢量为：飞行方向速度大小：vpa⊥＝vpa·j(21)飞行路径角：真近点角：半长轴：平面投影后的可能轨控点的位置矢量为：其中，投影后的可能轨控点的位置矢量的模注：将平面投影后的轨道参数用上标-标识，下同。平面投影后的可能轨控点的速度矢量为：2、虚拟交会点求解可能轨控点与交会点间的最小速度增量时，虽然已经用平面投影法克服法向偏差，但可能轨控点用二体模型外推到交会点方向和精确模型外推到交会点方向时两者还是会有径向偏差，这也是由于摄动带来的。而lagrange转移时间方程只针对二体情况建立两点转移方程，要克服摄动带来的径向偏差，通过建立虚拟交会点后与可能轨控点求解两点最小速度增量，所谓虚拟的交会点就是事先计算出二体模型与精确模型外推后的偏差矢量，然后将此偏差叠加到当前交会点后形成的交会点，具体计算如下：可能轨控点当前的交线幅角差为δupa，交会点对应的交线幅角差为δub，则两者角度差为δf0＝δub-δupa，从可能轨控点外推到交会点方向时的径向距离为：可能轨控点a用二体模型外推到交会点方向b时的位置矢量为：其中，可能轨控点a在考虑各种摄动情况下外推到交会点方向时的位置其中，两者差为：目标轨道交会点的位置为：其中，虚拟交会点为：虚拟交会点的模为3、转移角转移轨道的法向单位矢量：停泊轨道的法向单位矢量：停泊轨道与转移轨道的夹角为：依据速度增量小原则求解转移角δf以及轨道面的夹角余弦cosε，具体如下：4、求解速度增量最小的转移轨道在求解最小速度增量最小的转移轨道时，由于要用到lagrange转移时间方程，以及控后速度、速度增量、速度增量与转移半长轴等知识点，在此先阐述这里理论。■lagrange转移时间方程输入参数：两点弦长c：周长之半s：最小能量半长轴am：■lagrange当轨转移时间方程及α、β参数：lagrange转移时间方程为：ntr(tb-ta)＝(α-β)-(sinα-sinβ)(40)式中，ntr为转移轨道角速度，atr为转移轨道半长轴。记α0和β0为余弦主值，具体表示如下：其中，0≤α0≤π，0≤β0≤βm＜π，α和β依据象限判定如下：δf≤π且为第1类轨道转移，α＝α0，β＝β0；δf≤π且为第2类轨道转移，α＝2π-α0，β＝β0；δf＞π且为第1类轨道转移，α＝α0，β＝-β0；δf＞π且为第2类轨道转移，α＝α0，β＝-β0。■当转移半长轴已知，控后速度、飞行路径角求解如下：当转移角δf≠π：a点控后速度vtra为：转移轨道的半通径ptr为：etrsinftra为：转移轨道偏心率为：最终飞行路径角γtra：当转移角δf＝π：控后速度同(43)，转移轨道的半通径ptretrsinftra为：最终飞行路径角γtra：■速度增量δva求解：■速度增量δva与转移半长轴atr的导数求解：当δf≠π时，当δf＝π时，在和及转移角δf确定后，可求解出所有可能轨控点和交会点之间的最小速度增量，然后对最小速度增量统计，统计出速度增量最小的轨控点和对应的轨道，每一个可能轨控点与交会点的求解最小速度增量的步骤如下：1)计算lagrange转移时间方程的输入参数：s、c、am2)确定转移半长轴区间值，已知可能转移轨道的半长轴范围在[0.5ap1.5ap]，若0.5ap≤am，半长轴区间下限值为ax＝am，反之半长轴计算区间值下限值为ax＝0.5ap，若1.5ap≤am，则当前区间值无解，反之半长轴上限值为as＝1.5ap。3)当确定转移轨道半长轴区间[ax,as]，按照式(41)～(53)求解和然后根据和的符号特性进行求解轨控点和目标点的最小速度增量转移轨道。·当第1类轨道转移第2类轨道转移的则最小速度增量发生在atr＝ax时。·当第1类轨道转移第2类轨道转移的最小速度增量发生在第2类轨道转移：若第2类轨道转移的速度增量还是保持随半长轴的增加继续减小的特性，此时最小速度增量发生在atr＝as时，若最小速度增量发生在用对分法求解。·当第1类轨道转移第2类轨道转移的则最小速度增量在第1类轨道转移：若第1类轨道转移的速度增量还是保持随半长轴的增加继续减小的特性，此时最小速度增量在atr＝as时，若最小速度增量发生在用对分法求解。对分法的求解过程为：i.已知axi、asi，令azi＝(axi+asi)/2，计算其中i初值从0开始迭代，ax0＝ax,as0＝as；ii.若或|asi-axi|＜eps，停止计算，取atr＝azi，求取最小速度增量值iii.若取ax(i+1)＝axi,as(i+1)＝azi；若取ax(i+1)＝azi,as(i+1)＝asi，转i继续求解。按照上述过程，360个可能轨控点与交会点都能求解出最小速度增量，然后对360个最小速度增量统计，得出速度增量最小的值，记为min(δvmin)，对应的轨控点为最优轨控点，对应的转移为最优转移。(四)闭合轨道式的多次控制实施方案单脉冲轨道设计是在停泊轨道和目标轨道构形确定的情况，无时间的约束下求解出的速度增量最小的最优转移轨道和最优轨控点，通过这样的设计方式，只能保证此转移轨道能以最小的速度增量从最优轨控点转移到目标交会点上，但由于相位的关系，机动卫星与目标卫星一般情况下并不能在交会时刻交会，为此这里通过采用闭合轨道控制策略实现交会：将单脉冲设计时计算出的最小速度增量min(δvmin)在考虑推进能力的情况下分解成两次或更多次的控制，依据单脉冲设计时确定的轨控时刻计算出对应的停泊轨道的平跟数：平平近点角，第一次控制依据控制时长在上述的平平近点角对应时刻处对称控制，当有大于两次以上的控制时，则记录第一次控制结束时的平平近点角，依据第一次控制时长反推出控制中心点对应的平平近点角，然后将其作为下一次控制点，每次都按照上述步骤都依次类推计算对应下次的平平近点角。多次控制中的最后一次控制的中心点与单脉冲设计时的控制时刻相同，轨控时间间隔需满足最小轨控间隔等约束，由于增加的中间过渡轨道近似整轨，因此称为闭合轨道，通过闭合轨道调整相位，使机动卫星与目标卫星最终在交会点交会接近。参照伴星飞控经验，近地轨道控制一次轨控至少包含2轨测定轨数据，1轨计算复核时间，1轨注入变轨参数包，最后1轨执行轨控动作，因此这里暂定最小轨控间隔为5轨，若整个控制只有两次控制，则从开始下达命令到交会接近完成最小间隔最小值为10轨，则交会时刻的圈次k应满足k≥10的条件，当目标轨道在第k圈时，计算未轨控时停泊轨道对应值圈数，记为nz。在速度增量最小值和最优轨控点及交会点确定后，以两次控制为例，每次的速度增量求解如下：1)若最小轨控间隔为5圈，则闭合轨道圈数n1初始值为5，停泊轨道圈数为n0＝nz-n1，按照下述方法初步求解两次控制的速度增量，且求解出两次控制的速度增量总和不应超过min(δvmin)，为保证最后一次的控制精度，第一次控制的速度增量应大于一定的比例，因此δv1应与min(δvmin)同向同时，还应满足这里暂定bl＝0.5。当求解出的速度增量满足如上条件，则继续2)循环，否则n1＝n1+1，n0＝nt-n1，并判断n＜5是否成立，若成立，则k＝k+1，转到(三)继续求解，否则继续1)循环。当目标卫星到交会时刻的圈数为k，初始平平近点角为交会时刻的平平近点角有关平根数的参数用上标＝标识，忽略由于摄动引起的平近地点等偏转引起的时间差异，则整个交会时间为当最小速度增量最小值min(δvmin)未分解时，即单脉冲轨道方式控制时，相应的时间传递如下式：其中，是机动卫星的初始平平近点角；是机动卫星轨控位置的平平近点角；是停泊轨道的平半长轴；n0+n1是停泊轨道跑过的圈数；是最小速度增量最小值min(δvmin)未分解时到对应交会时刻的机动卫星在转移轨道的相位；是转移轨道轨控点平平近点角；是目标轨道交会时刻的平平近点角；是目标轨道初始平平近点角；是目标轨道的平半长轴；注：式(54)及式(55)的分析忽略由于摄动引起的平近地点等偏转引起的时间差异。最小速度增量最小值min(δvmin)分解成两次控制后，目标卫星和机动卫星正好到交会点上，对应的时间传递如下式：其中，是闭合轨道的平半长轴；是转移轨道交会时刻的平平近点角；n1是闭合轨道跑过的圈数；将式(55)减式(54)可得：整理可得：其中，为保证求解出的δv1与min(δvmin)同向，的符号应与min(δvmin)的切向分量符号相同，因此依据和min(δvmin)切向分量方向，进行如下处理：若min(δvmin)切向分量>0，则时不做处理，否则若min(δvmin)切向分量<0，则时不做处理，否则提取第一次轨控控前平半长轴平偏心率平真近点角则地心距为：第一次轨控控前半通径为：第一次轨控控前卫星速度矢量：控后卫星速度平方为：第一次轨控主要功能是调相，因此按切向控制，由此可得第一次控制的速度增量矢量(lvlh坐标系)：求出第一次轨控速度增量矢量后，第二次轨控的速度增量矢量为：δv2＝min(δvmin)-δv1(63)2)第1)步求解的两次控制速度增量是粗值，将其代入模型精确外推后，得到过渡轨道，将过渡轨道的平轨道参数用标识。从精确外推数据里取出交会时刻转移轨道的相位此时是小量，不用再像1)进行符号判断处理，计算然后取出过渡轨道半长轴转移轨道半长轴依据下式进行计算：其中，是上一次迭代后到对应交会时刻的机动卫星在转移轨道的相位。然后更新第一次速度增量：将更新后的δv1代入模型继续外推，δv2不变，取出第二次控后的数据到交会时刻的数据，依照平面投影法计算出以交线为起点的坐标系下的位置和数据，然后求解目标星在交会时刻时与机动卫星的距离，其值小于1m时，停止迭代，否则继续按照单脉冲轨道设计的方法，先求出虚拟交会点，然后求出第二次轨控结束时刻的位置和速度，依据此两点求解最小速度增量虚拟交会点和第二次轨控结束时刻的点依然采用平面投影后的位置和速度，最后更新第二次速度增量：将更新后的速度增量δv2代入模型继续精确外推，继续2)循环。整个流程见图1，上述整个求解是将单次控制分解成两次控制进行计算的，由于推力器执行时间有限的关系，有时候需要通过两次以上的控制实现，无论增加多少次控制，闭合控制主要功能是调相，使交会时刻的相位与预期相位相符，即mtrx＝mtrb。各个控制点都按照前述的平平近点角迭代方案确定。(五)摄动情况下固定时长的速度增量求解按照闭合轨道式的多次控制实施方案执行轨控后，由于测定轨误差、姿轨控控制误差会影响控制效果，因此最后一次控制时需按照固定时长转移进行速度增量迭代求解，其过程如下：1)由测定轨得到的精准轨道数据；2)轨控时刻初始值为上一步策略计算出的轨控时刻，初始速度增量为0，依据外推模型，计算出交会时刻以及轨控时刻到交会时刻的时间间隔，即转移时间。3)在平面投影后计算相同转移时间时内，由于停泊轨道摄动带来的位置矢量偏差，需注意的是，这一点与闭合轨道式的多次控制实施方案的计算方案不一样，在闭合轨道方案里，是计算相同转移角带来的位置矢量偏差。4)将位置矢量差叠加到交会位置上得到虚拟交会目标点，当转移两点及转移时间确定时，按照二分法求出转移轨道及对应的速度增量，将计算出的速度增量叠加到上一步迭代计算出的速度增量上。5)将更新后的速度增量代入模型，计算交会时刻和转移时间，当轨道机动卫星在交会时刻与目标卫星的距离，当此距离小于1m时，停止迭代，否则继续3)过程。以下是基于本实施例方法的一个仿真实验及其结果。目的：机动卫星与目标卫星交会到0km。(一)卫星参数机动卫星参数：1)卫星初始质量：100kg2)迎风面积：0.3m23)辐射面积：0.3m24)推进参数：5)推力：20n单次最长工作时长：1200s单次提供的最大速度增量：约240m/s初始燃料：31.5kg比冲：2156ns/kg6)最小轨控间隔：7h机动卫星的轨道参数是降交点为12点时的太阳同步轨道，到6oct201702:17:11.000utcg的轨道参数具体见表1。目标卫星参数：质量：1000kg迎风面积：20m2辐射面积：20m2机动卫星和目标卫星的初始参数见下表所示，实际实施时由地面测定轨提供。表1初始轨道参数(瞬根数)参数机动卫星目标卫星差值半长轴(km)7087.2807490.482403.202偏心率0.0004580.0100130.009555倾角(°)98.21463.394-34.820升交点赤经(°)14.701315.061300.360近地点幅角(°)359.726173.629-186.097纬度幅角(°)359.905359.781-0.124真近点角(°)0.179186.152185.973平近点角(°)0.178186.276186.98(二)仿真工具与模型整个仿真工具及参数设置、仿真工况如下：1)仿真工具stk10：astrogator模型matlab2010b2)模型参数引力模型采用wgs84_egm96，参数设置：degree:70，order:70。大气模型采用msise1990，参数按照stk默认设置。光压模型采用sphericalsrp，参数按照stk默认设置。第三体引力考虑太阳和月亮。机动卫星与目标卫星的大气阻力系数cd为2.2。机动卫星与目标卫星的辐射系数cr为1。(三)标称仿真结果按照摄动情况下单脉冲轨道设计，速度增量最小的控制点在δupa＝179°，交会点在δub＝0°，速度增量最小值矢量为min(δvmin)＝[-4.389×10-40.0980]km/s。通过matlab和stk的联合仿真，按照闭合轨道式多次轨道控制策略，可以得到：机动卫星经过2次地面控制，耗时0.9465天，在7oct201701:00:09.873utcg与目标卫星交会，总共需速度增量98.2290m/s，交会图见图2，其中标号，1表示目标轨道、2表示停泊轨道、3表示闭合轨道、4表示交会轨道、5表示交会点，具体控制参数和仿真结果如下：表2控制参数可以看出，机动卫星在第一次接近目标卫星时完成位置交会，交会精度0.0002111km。从表1，目标卫星半长轴比机动卫星半长轴高403km，目标卫星倾角比机动卫星倾角小-34.820°，目标卫星升交点赤经比机动卫星升交点赤经多300.360°，通过换算得到机动卫星与目标卫星轨道面夹角为67.462°，从上述可以看到，机动卫星与目标卫星初始轨道高度差异大，两者轨道面夹角也大，但通过本实施例的方法，只需98.2290m/s就可实现使机动卫星与目标卫星交会接近，这也再次印证了本发明适用于一般卫星甚至微纳卫星的交会接近任务。实施例二实施例一是机动卫星与目标卫星异轨面的计算步骤，本实施针对机动卫星与目标卫星同轨道面，其计算步骤与实施例一基本一致，但在单脉冲轨道设计寻找速度增量最小不一样，具体如下：异轨道面时只需将交线地心距差较小的那个交线点作为交会点即可，但同轨道时无这个规律，由于同轨道面无交线矢量，因此目标卫星轨数的统计直接以纬度幅角进行统计，即当纬度幅角在0～360度覆盖一次时，圈数增加一圈，当目标卫星将第k轨作为可能交会圈次时，将第k轨目标轨道均分360个点作为可能交会点，停泊轨道以目标轨道第k轨后的最后一点为结束时刻往回取出一轨数据，也均分成360个点作为可能轨控点，为了克服法向摄动的偏差，同样用平面投影法计算投影后的位置速度，但对应的角度不是交线幅角差而是纬度幅角。然后从目标轨道依次取出1点，从停泊轨道中也取出一点，用虚拟交会点克服径向偏差，然后求解虚拟交会点和可能轨控点这两点的最小速度增量，最终可得到360×360个最小速度增量值，然后从中统计出速度增量最小值min(δvmin)，由此对应的可能轨控点为最优轨控点，可能交会点为最优交会点。以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本
技术领域：
中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。当前第1页12