主权项 |
一种深空探测器与小行星交会成像的相对构像方法,其特征在于包括如下步骤:1)针对深空探测器在无姿态调整、无轨道维持的自由飞行状态下与小行星的近距离交会模式,通过以下六步建立深空探测器与小行星之间的空间相对运动关系:(1)在质心天球参考系下,基于深空探测器的轨道根数和小行星的轨道根数,按照以下公式分别计算交会成像期间深空探测器和小行星各自的引力加速度,以及深空探测器和小行星之间的相对加速度;所述引力加速度和相对加速度均由质心天球参考系中太阳引力产生;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>a</mi><mo>→</mo></mover><mi>s</mi></msub><mo>=</mo><mo>-</mo><mfrac><mrow><msub><mi>GM</mi><mrow><mi>S</mi><mi>u</mi><mi>n</mi></mrow></msub></mrow><mrow><mo>|</mo><msub><mover><mi>r</mi><mo>→</mo></mover><mi>s</mi></msub><msup><mo>|</mo><mn>3</mn></msup></mrow></mfrac><msub><mover><mi>r</mi><mo>→</mo></mover><mi>s</mi></msub></mrow>]]></math><img file="FDA0000994815950000011.GIF" wi="333" he="168" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>a</mi><mo>→</mo></mover><mi>a</mi></msub><mo>=</mo><mo>-</mo><mfrac><mrow><msub><mi>GM</mi><mrow><mi>S</mi><mi>u</mi><mi>n</mi></mrow></msub></mrow><mrow><mo>|</mo><msub><mover><mi>r</mi><mo>→</mo></mover><mi>a</mi></msub><msup><mo>|</mo><mn>3</mn></msup></mrow></mfrac><msub><mover><mi>r</mi><mo>→</mo></mover><mi>a</mi></msub></mrow>]]></math><img file="FDA0000994815950000012.GIF" wi="341" he="167" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>Δ</mi><mi>a</mi><mo>=</mo><msub><mover><mi>a</mi><mo>→</mo></mover><mi>s</mi></msub><mo>-</mo><msub><mover><mi>a</mi><mo>→</mo></mover><mi>a</mi></msub><mo>=</mo><mo>-</mo><mfrac><mrow><msub><mi>GM</mi><mrow><mi>S</mi><mi>u</mi><mi>n</mi></mrow></msub></mrow><mrow><mo>|</mo><msub><mover><mi>r</mi><mo>→</mo></mover><mi>s</mi></msub><msup><mo>|</mo><mn>3</mn></msup></mrow></mfrac><msub><mover><mi>r</mi><mo>→</mo></mover><mi>s</mi></msub><mo>+</mo><mfrac><mrow><msub><mi>GM</mi><mrow><mi>S</mi><mi>u</mi><mi>n</mi></mrow></msub></mrow><mrow><mo>|</mo><msub><mover><mi>r</mi><mo>→</mo></mover><mi>a</mi></msub><msup><mo>|</mo><mn>3</mn></msup></mrow></mfrac><msub><mover><mi>r</mi><mo>→</mo></mover><mi>a</mi></msub></mrow>]]></math><img file="FDA0000994815950000013.GIF" wi="820" he="174" /></maths>其中,a<sub>s</sub>为质心天球参考系中太阳对深空探测器产生的引力加速度,a<sub>a</sub>为质心天球参考系中同一时刻太阳对小行星产生的引力加速度,Δa为深空探测器和小行星间的相对加速度;GM<sub>Sun</sub>为太阳的引力常数;<img file="FDA0000994815950000014.GIF" wi="38" he="63" />和<img file="FDA0000994815950000015.GIF" wi="40" he="71" />分别为同一时刻探测器和小行星相对于太阳质心的位置矢量;(2)基于小行星已知质量,计算交会成像期间小行星和深空探测器之间的万有引力对深空探测器速度的影响量;具体步骤如下:建立相对引力坐标系O<sub>a</sub>‑xyz,其原点O<sub>a</sub>为小行星质心,x方向为相对运动方向,y方向为交会时刻卫星相对位置的反方向,z轴与x,y构成右手系;将交会时刻记为t<sub>0</sub>,根据坐标系定义,该时刻小行星坐标可记为(0,y<sub>0</sub>,0);在O<sub>a</sub>‑xyz坐标系中,小行星引力产生的速度变化表达如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mi>δ</mi><msub><mi>v</mi><mi>x</mi></msub><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mi>δ</mi><msub><mi>v</mi><mi>y</mi></msub><mo>=</mo><mfrac><mrow><mn>2</mn><msub><mi>μy</mi><mn>0</mn></msub><mover><mi>x</mi><mo>·</mo></mover><mi>Δ</mi><mi>t</mi></mrow><mrow><mover><mi>x</mi><mo>·</mo></mover><msubsup><mi>y</mi><mn>0</mn><mn>2</mn></msubsup><msqrt><mrow><msup><mrow><mo>(</mo><mover><mi>x</mi><mo>·</mo></mover><mi>Δ</mi><mi>t</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msubsup><mi>y</mi><mn>0</mn><mn>2</mn></msubsup></mrow></msqrt></mrow></mfrac></mtd></mtr><mtr><mtd><mrow><msub><mi>δv</mi><mi>z</mi></msub><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000994815950000016.GIF" wi="525" he="391" /></maths>其中,μ=GM<sub>A</sub>为小行星的引力常数,<img file="FDA0000994815950000017.GIF" wi="37" he="47" />为x坐标方向的速度,Δt为交会时刻到飞掠过后时刻的总时间;(3)如果步骤(1)中相对加速度小于0.01m/s<sup>2</sup>,并且步骤(2)中小行星和深空探测器之间的万有引力对深空探测器速度的影响量小于0.1%,则用匀速直线运动来描述深空探测器和小行星相对位置关系:小行星为静止质点,深空探测器相对其作匀速直线运动;或深空探测器为静止质点,小行星相对其作匀速直线运动;(4)在质心天球参考系下,根据以下公式计算交会成像期间小行星的空间姿态变化量Δθ;<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>Δ</mi><mi>θ</mi><mo>=</mo><mfrac><mi>t</mi><msub><mi>T</mi><mi>i</mi></msub></mfrac><mi>π</mi></mrow>]]></math><img file="FDA0000994815950000021.GIF" wi="222" he="127" /></maths>其中,T<sub>i</sub>为小行星绕其第i个旋转轴的旋转周期,t为交会成像时间;(5)在质心天球参考系下,根据交会成像期间深空探测器的姿态遥测数据和已知的星载相机安装矩阵,按照以下公式计算交会成像期间深空探测器相机的空间角度变化量;M(Δα<sub>c</sub>,Δβ<sub>c</sub>,Δχ<sub>c</sub>)=M<sub>bp</sub>·C(Δα)·C(Δβ)·C(Δχ)其中,M<sub>bp</sub>为相机相对本体系的安装矩阵,Δα,Δβ,Δχ分别为遥测得到的深空探测器交会成像期间俯仰、偏航和滚转姿态变化量,C(Δα)、C(Δβ)、C(Δχ)分别为由Δα,Δβ,Δχ生成的旋转矩阵,Δα<sub>c</sub>,Δβ<sub>c</sub>,Δχ<sub>c</sub>分别为成像期间相机的俯仰、偏航和滚转姿态变化量;(6)如果步骤(4)中小行星的空间姿态变化量小于0.1度,并且步骤(5)中相机的空间角度变化量小于0.1度,则用相对无旋转模式来描述交会成像期间深空探测器和小行星的相对姿态关系:小行星相对深空探测器的姿态角为常值,二者之间不存在相对旋转;2)通过以下三步建立深空探测器相机载荷所拍摄的小行星序列图象的构像关系:(1)对于纯平移运动成像,根据CCD相机中心投影规律,建立平移运动下成像投影关系,具体如下:平移运动物体上任何一对对应点的连线与平移方向平行,所有对应点的连线在空间中相交于无穷远点P<sub>∞</sub>,该无穷远点P<sub>∞</sub>表示平移方向;目标在不同位置对应点的成像点连线在图像中交于汇聚点e,该汇聚点e是平移方向P<sub>∞</sub>在图像平面上的投影;(2)结合深空探测器与小行星相对运动关系,建立交会成像期间相机序列图像上的小行星目标构像描述关系:按照步骤2)中第(1)步所述方法建立深空探测器对小行星的构像关系;所述汇聚点e=(u<sub>e</sub>,v<sub>e</sub>,1)′;其中:u<sub>e</sub>、v<sub>e</sub>分别为小行星序列图像中汇聚点的像面坐标;(3)建立小行星序列图象间基础矩阵F<sub>m,n</sub>的表达关系:F<sub>m,n</sub>=[e<sub>n,m</sub>]<sub>×</sub>F<sub>m,n</sub>=F=[e]<sub>×</sub><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mrow><mo>[</mo><msub><mi>e</mi><mrow><mi>n</mi><mo>,</mo><mi>m</mi></mrow></msub><mo>]</mo></mrow><mo>×</mo></msub><mo>=</mo><msub><mrow><mo>[</mo><mi>e</mi><mo>]</mo></mrow><mo>×</mo></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd><mtd><msub><mi>v</mi><mi>e</mi></msub></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><msub><mi>u</mi><mi>e</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><msub><mi>v</mi><mi>e</mi></msub></mrow></mtd><mtd><msub><mi>u</mi><mi>e</mi></msub></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000994815950000031.GIF" wi="662" he="244" /></maths>其中,[e<sub>n,m</sub>]<sub>×</sub>表示由汇聚点e构造的反对称矩阵。 |