发明名称 一种卫星对月定标的时间规划方法
摘要 一种卫星对月定标的时间规划方法,包括以下步骤:利用卫星工具软件初步确定卫星对月定标时间;根据公式确定对月定标具体开始时间、结束时间和相机曝光时间;根据星历获取月球位置参数和卫星轨道要素,由公式求得卫星定标期间各个时刻的补偿角速度,求取该角速度的最大和最小值;选取回到定标初始位置的不同方式,计算两种方式下的机动时间,最终确定卫星定标次数,完成卫星对月定标的时间规划。
申请公布号 CN105823485A 申请公布日期 2016.08.03
申请号 CN201610262686.0 申请日期 2016.04.25
申请人 航天东方红卫星有限公司 发明人 白照广;高涵;邸国栋;单悌磊
分类号 G01C21/24(2006.01)I;G01C25/00(2006.01)I 主分类号 G01C21/24(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 陈鹏
主权项 一种卫星对绝对月定标的时间规划方法,其特征在于步骤如下:(1)使用STK建立仿真场景模型,模型中包括月球,地球,太阳,卫星和星上传感器;设置卫星轨道参数、姿态机动能力参数、星上传感器参数和卫星初始姿态;(2)根据步骤(1)中的场景特征,初步确定对月定标时间窗口;利用STK设置月相角在[3°,7°]时卫星对月球全部可见时间段T1、T2……Tn,作为初步卫星对月定标时间窗口;(3)以步骤(2)中确定的定标时间窗口为基础,精确确定对月定标时间窗口;求取卫星进入地影区时间t1和进入阳照区时间t2,[t<sub>1</sub>,t<sub>2</sub>]即为精确对月定标时间窗口;(4)计算t时刻卫星需要补偿的角速度dθ,t∈[t1,t2],具体计算方法如下:求取卫星与月球的距离:L=S‑(R<sub>E</sub>+h)sin(ωt+β);其中S为地心到月心的距离,R<sub>E</sub>为地球平均半径,h为卫星轨道高度,ω为角速度,β是黄白交角;求取卫星定标期间在月球表面的地面采样间隔GSD:<maths num="0001"><math><![CDATA[<mrow><mi>G</mi><mi>S</mi><mi>D</mi><mo>=</mo><mfrac><mrow><msub><mi>&pi;R</mi><mi>M</mi></msub><msup><mn>10</mn><mrow><mo>-</mo><mn>3</mn></mrow></msup></mrow><mn>180</mn></mfrac><mo>&lsqb;</mo><mi>arcsin</mi><mfrac><mrow><mi>L</mi><mi> </mi><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mi>&alpha;</mi><mo>+</mo><mfrac><mi>p</mi><mrow><mn>2</mn><mi>f</mi></mrow></mfrac><mo>)</mo></mrow></mrow><msub><mi>R</mi><mi>M</mi></msub></mfrac><mo>-</mo><mi>arcsin</mi><mfrac><mrow><mi>L</mi><mi> </mi><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mi>&alpha;</mi><mo>-</mo><mfrac><mi>p</mi><mrow><mn>2</mn><mi>f</mi></mrow></mfrac><mo>)</mo></mrow></mrow><msub><mi>R</mi><mi>M</mi></msub></mfrac><mo>-</mo><mfrac><mi>p</mi><mi>f</mi></mfrac><mo>&rsqb;</mo><mo>;</mo></mrow>]]></math><img file="FDA0000973739530000011.GIF" wi="1438" he="269" /></maths>其中R<sub>M</sub>是月球平均半径;α为卫星侧摆角;p为像元尺寸,f为焦距;求取卫星的星下点地速:<maths num="0002"><math><![CDATA[<mrow><msub><mi>V</mi><mi>n</mi></msub><mo>=</mo><msub><mi>R</mi><mi>E</mi></msub><msqrt><mfrac><mi>&mu;</mi><msup><mrow><mo>(</mo><msub><mi>R</mi><mi>E</mi></msub><mo>+</mo><mi>h</mi><mo>)</mo></mrow><mn>3</mn></msup></mfrac></msqrt><mo>;</mo></mrow>]]></math><img file="FDA0000973739530000012.GIF" wi="405" he="150" /></maths>其中μ为地球引力常数,μ=3.986×10<sup>5</sup>km<sup>3</sup>/s<sup>2</sup>;求取理论积分时间:<maths num="0003"><math><![CDATA[<mrow><msub><mi>T</mi><mi>int</mi></msub><mo>=</mo><mfrac><mrow><msup><mi>GSD</mi><mo>&prime;</mo></msup></mrow><msub><mi>V</mi><mi>n</mi></msub></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000973739530000021.GIF" wi="253" he="134" /></maths>其中GSD′是相机的地面采样间隔;求取实际推扫速度和补偿速度:v=ωL;<maths num="0004"><math><![CDATA[<mrow><mi>V</mi><mo>=</mo><mfrac><mrow><mi>G</mi><mi>S</mi><mi>D</mi></mrow><msub><mi>T</mi><mi>int</mi></msub></mfrac></mrow>]]></math><img file="FDA0000973739530000022.GIF" wi="196" he="127" /></maths>计算获得补偿角速度大小:<maths num="0005"><math><![CDATA[<mrow><mi>d</mi><mi>&theta;</mi><mo>=</mo><mfrac><mrow><mi>V</mi><mo>&PlusMinus;</mo><mi>v</mi></mrow><mrow><mi>L</mi><mo>-</mo><msub><mi>R</mi><mi>M</mi></msub><mo>-</mo><msub><mi>R</mi><mi>E</mi></msub></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000973739530000023.GIF" wi="420" he="135" /></maths>(5)根据两种不同的成像方式,计算单次定标时间TT1,具体为:定义对月定标中成像方式分为两种,第一种为先进行滚转机动,再进行俯仰机动,推扫方向与卫星运动方向相同,第二种为先进行俯仰机动,再进行滚转机动,推扫方向与卫星运动方向相反;当选择第一种时,卫星补偿角速度大小<img file="FDA0000973739530000024.GIF" wi="406" he="135" />当选择第二种时,卫星补偿角速度大小<img file="FDA0000973739530000025.GIF" wi="387" he="127" />设ω<sub>min(补)</sub>为卫星最小补偿角速度,即第一种方式处于赤道时的补偿角速度,以及第二种方式处于两极时的补偿角速度;ω<sub>max(补)</sub>为卫星最大补偿角速度,即第一种方式处于两极时的补偿角速度,以及第二种方式处于赤道时的补偿角速度;设卫星最大俯仰机动角速度为ω<sub>Y_max</sub>,ω<sub>Y_max</sub>&gt;ω<sub>max(补)</sub>,俯仰最大姿态机动角加速度为a<sub>Y</sub>,t<sub>a</sub>为卫星加速时间,t<sub>s</sub>为卫星稳定时间,t<sub>d</sub>为卫星减速时间,则在精确对月定标时间窗口[t1,t2]内单次定标时间TT1的求法如下:<img file="FDA0000973739530000026.GIF" wi="555" he="127" /><img file="FDA0000973739530000027.GIF" wi="566" he="126" />在月面推扫的最长用时t<sub>scan_max</sub>为:<maths num="0006"><math><![CDATA[<mrow><msub><mi>t</mi><mrow><mi>s</mi><mi>c</mi><mi>a</mi><mi>n</mi><mo>_</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn>2</mn><msub><mi>R</mi><mi>M</mi></msub></mrow><mrow><mi>m</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000973739530000031.GIF" wi="358" he="127" /></maths>在月面推扫的最短用时t<sub>scan_min</sub>为:<maths num="0007"><math><![CDATA[<mrow><msub><mi>t</mi><mrow><mi>s</mi><mi>c</mi><mi>a</mi><mi>n</mi><mo>_</mo><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn>2</mn><msub><mi>R</mi><mi>M</mi></msub></mrow><mrow><mi>m</mi><mi>a</mi><mi>x</mi><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000973739530000032.GIF" wi="362" he="133" /></maths>当定标次数为理论最少定标次数时,最大单次定标时间TT1<sub>max</sub>:TT1<sub>max</sub>=t<sub>amax</sub>+t<sub>dmax</sub>+t<sub>scan_max</sub>+t<sub>s</sub>,<img file="FDA0000973739530000033.GIF" wi="469" he="124" />当定标次数为理论最多定标次数时,最小单次定标时间TT1<sub>min</sub>:TT1<sub>min</sub>=t<sub>amin</sub>+t<sub>dmin</sub>+t<sub>scan_min</sub>+t<sub>s</sub>,<img file="FDA0000973739530000034.GIF" wi="445" he="119" />(6)使卫星视轴回到定标初始位置,并计算卫星回到初始位置的最长机动时间TT2<sub>max</sub>和最短机动时间TT2<sub>min</sub>;(7)计算得到在一次定标过程中,卫星完成的定标次数。
地址 100094 北京市海淀区北京市5616信箱
您可能感兴趣的专利