发明名称 兼顾多任务高度的太阳同步轨道倾角偏置方法
摘要 本发明公开了兼顾多任务高度的太阳同步轨道倾角偏置方法,根据子任务的工作高度H<sub>j</sub>,计算对应的考虑J<sub>4</sub>项设计满足太阳同步轨道特性的标称倾角<img file="DDA00002991098100011.GIF" wi="71" he="78" />降交点地方时初值使用总任务LTDN<sup>0</sup>作为输入,使用迭代方法,计算在一种倾角偏置量Δi下的第一次子任务降交点地方时LTDN<sub>1</sub>(t)和倾角i<sub>1</sub>(t),并将上一次子任务结束时的LTDN<sub>j</sub>(T<sub>j</sub>)和i<sub>j</sub>(T<sub>j</sub>)作为本次子任务的输入条件,按照j=1~N的排列,逐段完成所有子任务的降交点地方时漂移分析;最终以Δi为总任务的自变量,寿命期内降交点地方时漂移正负幅值差作为目标函数,当目标函数取值最小时,得到的倾角偏置量Δi为最优方案。本发明能够使得总任务期间的降交点地方时综合漂移变化最小。
申请公布号 CN103235870B 申请公布日期 2016.06.01
申请号 CN201310108729.6 申请日期 2013.03.29
申请人 航天东方红卫星有限公司 发明人 李志武;白照广;谭田;张燕;吕秋杰
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 中国航天科技专利中心 11009 代理人 褚鹏蛟
主权项 一种兼顾多任务高度的太阳同步轨道倾角偏置方法,所述多任务高度包括N次子任务的轨道高度H<sub>j</sub>,j=1~N;第j次子任务所占用的时间区间为T<sub>j‑1</sub>≤t≤T<sub>j</sub>,T<sub>0</sub>=0,其特征在于,步骤如下:(1)根据子任务的轨道高度H<sub>j</sub>计算该子任务的太阳同步轨道标称倾角<img file="FDA0000914939380000011.GIF" wi="70" he="79" />j=1~N;所述标称倾角<img file="FDA0000914939380000012.GIF" wi="46" he="79" />为在考虑J<sub>4</sub>项设计下满足太阳同步轨道特性的标称倾角;J<sub>4</sub>为地球引力四阶带谐项系数;(2)根据标称降交点地方时LTDN<sup>0</sup>设定太阳同步轨道倾角偏置量Δi的初始值和收敛方向;并设定倾角偏置量调整步长初值;所述倾角偏置量Δi为入轨倾角相对第一次子任务标称倾角<img file="FDA0000914939380000015.GIF" wi="40" he="64" />的太阳同步轨道倾角偏置量;(3)以倾角偏置量Δi和标称降交点地方时LTDN<sup>0</sup>作为输入,迭代计算第一次子任务期间倾角i<sub>1</sub>(t)和第一次子任务期间的降交点地方时LTDN<sub>1</sub>(t),t∈[0,T<sub>1</sub>];所述步骤(3)中迭代计算的结束条件为:t∈[0,T<sub>1</sub>]内第n+1次循环计算的所有时刻相对第n次循环计算的相同时刻均满足<img file="FDA0000914939380000013.GIF" wi="623" he="93" />(4)从第二次子任务开始,将上一次子任务结束时的降交点地方时LTDN<sub>j</sub>(T<sub>j</sub>)和倾角i<sub>j</sub>(T<sub>j</sub>)作为本次子任务的输入条件,迭代计算本次子任务期间的降交点地方时LTDN<sub>j+1</sub>(t)和倾角i<sub>j+1</sub>(t),T<sub>j</sub>≤t≤T<sub>j+1</sub>,直至完成所有子任务期间的降交点地方时和倾角计算;所述步骤(4)中迭代计算的结束条件为:T<sub>j</sub>≤t≤T<sub>j+1</sub>内第n+1次循环计算的所有时刻相对第n次循环计算的相同时刻均满足<maths num="0001"><math><![CDATA[<mrow><mo>|</mo><msubsup><mi>LTDN</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>LTDN</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo><mo>&lt;</mo><mn>1</mn><mi>s</mi><mo>;</mo></mrow>]]></math><img file="FDA0000914939380000014.GIF" wi="646" he="94" /></maths>(5)按时间顺序串联每一次子任务期间的降交点地方时LTDN<sub>j</sub>(t),T<sub>j‑1</sub>≤t≤T<sub>j</sub>,得出在当前倾角偏置量Δi下执行所有子任务期间降交点地方时LTDN(t),0≤t≤T<sub>N</sub>;将所有子任务期间降交点地方时LTDN(t)减去标称降交点地方时LTDN<sup>0</sup>获得所有子任务期间降交点地方时漂移;(6)判断所有子任务期间降交点地方时漂移的最大值的绝对值与最小值的绝对值之差的绝对值是否小于1s,如果小于1s,则此时的Δi为最佳倾角偏置量Δi<sub>f</sub>,转入步骤(7);如果不小于1s,则沿着所述收敛方向重新计算Δi,然后返回步骤(3);(7)以最佳倾角偏置量Δi<sub>f</sub>作为任务开始的倾角偏置量,入轨时的倾角为<maths num="0002"><math><![CDATA[<mrow><msubsup><mi>i</mi><mn>1</mn><mn>0</mn></msubsup><mo>+</mo><msub><mi>&Delta;i</mi><mi>f</mi></msub><mo>.</mo></mrow>]]></math><img file="FDA0000914939380000021.GIF" wi="187" he="85" /></maths>
地址 100094 北京市海淀区北京市5616信箱