哈工大深空探测轨道设计作业_地球至火星轨道设计

1、目 录1.1研究现状及分析21.1.1 发射窗口41.1.2火星探测轨道设计51.1.3火星探测轨道优化71.2轨道基础知识91.2.1时间系统91.2.2坐标系统101.2.3星历数据111.2.4B平面111.2.5Lambert问题121.3火星探测直接转移轨道的初步设计131.3.1日心轨道设计及发射窗口的搜索131.3.2地心段参数的确定151.3.3火心段参数的确定191.4 基于B平面参数的精确轨道设计201.4.1 问题描述201.4.2 制导方法211.4.3 轨道精确设计求解221.4仿真分析231.4.1初步轨道参数设计结果241.4.2 精确轨道参数设计结果261.5结

2、论27I深空探测轨道设计课程读书报告地球火星转移轨道设计轨道设计是火星探测任务的基础，在设计出精确轨道前，一般都忽略次要因素，以二体模型为基础设计一条简单的轨道来满足任务的要求。本章采用普适变量方法求解Lambert问题，并给出基于pork-chop图以及优化算法两种方法对发射窗口进行搜索，基于此窗口对转移轨道进行初步设计和精确设计。1.1 研究现状及分析近十年来火星探测已成为科学家们开展空间研究的主流趋势之一，火星是太阳系内与地球最接近的一颗行星，它们有很多共同特征。自从水被证实在其上存在后，有存在生命的可能是人类目前对火星感兴趣的主要原因之一，此推动了科学研究，在之后每一个合适的发射窗口，

3、都有新型的行星际探测器飞往火星，并携带科学设备用来研究火星的大气与表面，以及发现一些新奇的现象。在过去的50年里，仅美国在火星探测研究的经费已超过了100亿美金，而在不远的将来他们计划开展大量的火星科学探测活动。目前，包括俄罗斯航天局在内的世界各大航天机构正在考虑发射载人探测器到火星上的可能性，而确定这样的计划后使得火星探测基础理论研究、技术支持和工程实验迅猛发展，此时我国开展火星探测是及时的，在自主研发的基础上，借鉴外国经验，发展我国自己的火星探测技术，开拓空间资源和领域，促使太空经济蓬勃发展。截至到2013年，人类从“火星1号”开始共发射41颗火星探测器，其中按任务类型可分为：飞越型、环绕

4、型和着陆型三类。下表给出了部分抵达过火星的探测器数据（发射与到达时火星和地球的黄经差分别为和、转移时间和转移角度）。表 1 部分抵达火星的探测器数据探测器任务类型火星1号飞越46.34-53.01220161.61水手4号飞越58.74-46.6228178.22火星2号轨道/着陆器47.95-40.77192148.56水手9号轨道器42.83-34.60168130.86火星4号轨道器41.46-52.63204148.37火星7号着陆器32.64-65.1212143.75海盗1号轨道/着陆器50.59-103.95335226.16火卫一2号轨道/着陆器36.64-56.0820114

5、1.96环球勘测者轨道器73.51-68.77308234.90火星探路者着陆/火星车61.03-36.86212172.15希望号轨道器155.01-43.361990-奥德赛号轨道器49.46-42.98200154.25快车轨道/着陆器46.56-48.52206154.49勇气号火星车42.89-53.16208151.9机遇号火星车29.95-62.88201135.27侦查轨道器轨道器36.37-60.55210146.43凤凰号着陆器63.12-73.01295217.7上表中转移角基本上都在180度附近，类似于霍曼转移轨道的结果，反映出火星轨道设计优先考虑能量问题，这为我国自行

6、开展火星探测任务提供了参考信息。总结过去50多年的火星探测任务，人类借助多个探测器对火星进行了观测与探测，揭示了一个与地球相近而又有众多不同的新世界。针对之前所取得的成果，世界各个航天大国纷纷提出自己的火星探测计划，其中以美国NASA火星科学实验室和欧空局的天外火星最具有代表性，而其主要任务为：1）开展有针对性的火星侦测任务，寻找火星上的生命迹象和热液源，及火星大气高层探测，即空间生物研究实验室和火星大气高层探测卫星；2）火星样本返回计划，采用相对简单的方法在着陆点就地采集土壤和大气样品，通过返回式飞船带回地球进行详细分析，最理想情况下这一任务将于2016年欧空局开展；3）在火星安置寿命较长的

7、着陆器形成网络，进而开展地震学、地质化学和气象学研究，评估火星上微粒的特性及大气参数的动态变化，乃至于研究航天员到达火星的生物危险等级；4）开展深度钻孔任务，即配备钻孔设备对火星土壤深度挖掘，寻找可能存在的水资源并评估其特性，以期望在未来的载人登火任务中可以充分利用火星上的资源；5）开启载人登陆火星任务，主要分三个阶段，初始阶段主要验证就地资源的使用可行性和高空拦截技术（70°锥角）及在火星大气机动时用仪器测量环境参数，中间阶段主要解决精确着陆和火星表面土壤的辐射防护特性分析，及验证精简模式的载人火星着陆系统，最终阶段需要解决首次载人任务的着陆点勘测和关键任务的全比例实验。总的来说，

8、在载人探测之前的无人探测任务还有许多待解决的问题，火星上水和甲烷的发现加速了世界各国科学家们火星移民和开发火星的设想，根据过去五十年来对火星环境的信息已为我国提供了宝贵的数据支持，开展火星探测研究可以提升我国科学和经济发展，更加可以提高我国在社会的地位，所以在国家“十一五”计划中早早地规划了以月球探测为基础，继而开展火星探测为主线的深空探测任务。1.1.1 发射窗口大多数火星转移轨道均是采用霍曼过渡方式，而在地球上每隔26个月才会出现一个较为合适的发射窗口，此时地球与火星相对比较近，发射能量较少，一旦错过这个合适的发射窗口，地球与火星在日心坐标系下的位置发生了变化，飞行路线也发生变化，导致能量

9、不是最优，甚至不能临近火星，故发射窗口的选取是火星探测任务中重要的一环，主要是受到诸多限制，而这些限制条件与发射成本和工程复杂性有关，主要包括飞行时间、发射能量与任务总能量，甚至包括运载火箭的级别，而在合适的年份探测火星，火星探测器对运载火箭是有一些特定的要求，无论运载能力是否足够，一些运载火箭仍然不能发射火星探测器。选取最优发射窗口就显得格外重要，工程上应用最为广泛的是等高线图法，即绘制pork-chop图，其能够直观地描述出在既定的时间段内发射窗口的变化情况，为满足约束条件的最优发射窗口提供较为精确的初值，这种方法最早出现于1983年Sergeyevsky 给出了1991年至2005年金星

10、探测的最优发射窗口选取中，鉴于此优点，国内外大部分学者一般采用该方法设计行星探测发射窗口的初值，但针对较长时间段的搜索，这种穷举搜索法的计算量极大，给探测设计任务带来了诸多不便。针对此缺点，国内哈尔滨工业大学的乔栋博士提出了一种基于遗传算法的最优发射窗口搜索， 其通过对星历与Guess问题的解算，将上述问题简化为仅含有两个变量的寻优问题，该方法具有全局搜索的能力，并且收敛快，作者给出了基于此方法的小行星的最优发射窗口的搜索时间仅为传统方法的4.19%，极大地提高计算效率，但往往需要多次迭代才能够得到精确的发射窗口，之后，作者又提出了遗传算法与序列二次规划方法相结合来搜索发射窗口， 其能够快速得

11、到精确的发射窗口。而在小推力火星探测发射窗口搜索方面，国外的Petropoulos 提出了一种基于形状的方法对发射窗口快速搜索，其能在二维轨道面快速地搜索到一个发射窗口，但含有过多的冗余计算，随后国内的一些学者对此法进行合理选取搜索参数和加入一些约束要求，完整的给出了算法流程，并搜索了火星借力探测小行星的发射窗口，可仍然不是最优的发射窗口；此后，西北工业大学的岳晓奎教授结合了DE算法与间接法各自优点，使用这两种算法搜索到的发射窗口是一片区域，需根据实际情况综合考虑时间与燃料才能确定最优的发射窗口。我国作为火星探测任务的后起之秀，自主火星探测工程难度大，但其国际影响和战略意义甚大，而火星探测新一

12、轮浪潮中，我国要切实抓住宝贵的发射窗口，以确保在2015年实施自主火星能成功发射。1.1.2火星探测轨道设计对于火星探测轨道设计技术，按照能量获取方法分类，可大致分为三种方法:直接转移、小推力变轨和借力飞行的方法。直接转移方法是指，探测器在短时间内由大推力冲量的方法获得瞬时加速度，改变探测器的速度，通过一次大推力加速过程直接完成探测器飞往目标天体的转移。对于直接转移方法，可以分为初步设计和精确设计两个过程。对于初步设计，是以圆锥曲线拼接法(Patched Conic Method)为基本原理，将探测轨道划分为几个分段过程，每个分段过程可以近似成二体问题处理，然后通过接口，将每一段拼接在一起。B

13、reakwell和Perko证明了针对行星际轨道设计问题，应用圆锥曲线拼接法可以得到合理结果。对于二体问题，Sergeyevsky和Yin提出了针对不同发射时间和到达时间，绘制发射能量和到达能量以求得发射机会的方法Pork-Chop图法，来搜寻可能的发射机会，Pork-Chop图法是搜索发射机会的经典方法。能量等高线图一般以发射和到达时间为坐标轴绘制发射或到达能量的等高线图，图上的每一点都可以由相对应的兰伯特(Lambert )问题求得。由于初步设计，采用的是二体问题处理的方法，没有考虑其他星体的引力、太阳风等因素，所以尽管圆锥曲线拼接法对于深空轨道初始设计与任务验证来说可以提供足够的精度，在

14、实际工程探测中，依然需要进行基于精确模型的计算以求得符合精度要求的转移轨道。Angelo Miele和T. Wang曾经采用连续梯度修补C Sequential Gradient Restoration)非线性规划算法研究火星探测任务的轨道特性。此种方法采取最优控制原理求解，方法的非线性较强。此外，大部分精确动力学模型算法采取的均为轨道搜索算法。包括不需要偏导数信息的变步长折回爬山法、可变容差多面体算法和采用偏导数信息的牛顿微分校正算法、最速下降梯度法等方法。小推力方法采用的是高比冲推进，通过长时间的小推力加速，螺旋运动，到达目标星体。由于小推力方法动力学模型长时间存在小推力项，采用的转移方案

15、和设计方法将会与常规的方案存在差异。小推力转移的优化属于函数空间的最优化控制问题，目前解决此问题主要存在两种方法。一是基于Pontryagin极小值原理的间接方法，该方法主要是通过变分方法求解必要条件，然后再通过数值方法求解边值问题。目前提出的间接方法包括梯度方法、拟线性化(Quasilinearization)方法和有限差分方法(Finite Difference Method。现阶段应用更多的是直接方法，主要是通过离散化处理，将连续问题变为有限参数优化问题，通过迭代的方法寻求最优解。现阶段主要应用的直接法包括与间接法相结合的方法、直接打靶法和遗传算法等。借力飞行是探测器在飞往目标星体的过程

16、中，接近其他星体，并利用该星体的引力改变自身轨道。通过借力飞行的方法，可以使探测器改变到理想的轨道，尤其是在飞往距地球较远的星体时，需要较大的能量，通过引力辅助变轨可以获得速度增量，以减少发射能量和飞行时探测器所消耗的能量。早在二十世纪五十年代B attin就己经开始研究往返星际航行中通过借力飞行的方法节省探测器能量消耗。J.K.Mille:运用蒂塞朗准则(Tisserand's Criterion)针对不存在动力近拱点的借力飞行进行了分析，找到了发射星体和借力星体、借力星体和目标星体之间成对的可行轨道。此外，Longuski提出了一种设计多天体交会借力飞行轨道的方法，通过给定的初始发

17、射时间段和目标星体，运用自动寻找C3的匹配(即飞出借力星体和飞入借力星体之间的能量匹配)找出满足所需条件的发射机会。同时Williams和Longuski还证明了此种方法可以找到一些新的、效率更高的轨道设计方案。在纯借力飞行之外，带有轨道机动的借力飞行轨道设计同样做了大量研究。Gobetz研究了借力飞行时附加单个机动的轨道转移方法，通过研究得到:在借力飞行期间采取机动(包括单一冲量情况)可增加任务的应变性;单一冲量情况下，在近拱点处附加机动可以优化转移点;在借力星体影响球内，双曲线转移可被简化成二维问题处理。此外，Wlaton, Marchal和Culpye同样对借力飞行期间的轨道转移进行了研

18、究，讨论了包含多个冲量的情况。在附加深空机动的研究方面，R. E. Diehl和M. R. Myers证明了包含深空机动的借力飞行方案可以减少能量，提高效率，并结合实例进行了轨道方案设计。随后，Moonish和Longuski又提出了一种包括深空机动的借力飞行发射窗口搜索方法。这种方法能够找出所有附加深空机动的转移方案和发射窗口，通过此方法设计的深空机动点可以使两个星体之间飞行的能量最小。在借力飞行轨道优化设计方面，Carl G. Sauer提出了一种基于主矢量原理的无机动借力飞行优化方法，针对多天体交会借力飞行问题，LouisA .D'-Amari.和Dennis V. Byrnes

19、等提出了带有约束的优化方法，能量匹配和飞越高度的约束通过罚函数体现，转化成无约束问题。J. Schoemnaekers提出了针对多天体交会的线性优化方法，并结合罗塞塔任务进行了计算。基于最优控制理论Dario Pastrone,LorenzoCasalino和Guido Colasurdo提出了搜索探测火星轨道发射窗口的优化方法。在国内方面，近几年各研究单位也纷纷开展火星探测的研究。清华大学、上海航天控制技术研究所、上海航天技术研究院、装备指挥技术学院等单位开展了基于B平面的精确动力学模型火星轨道设计，提出了精确动力学模型设计方案和轨道修正方法。哈尔滨工业大学、北京理工大学进行了小推力转移轨道

20、方面的研究，分别提出了地球火星最省小推力优化方案和借力飞行小推力方案。哈尔滨工业大学还对摄动作用对火星探测器的影响做了研究，讨论了多种不同摄动作用对于火星探测器轨道设计的不同影响。中国运载火箭技术研究院也展开了火星探测轨道的研究工作，进行了基于霍曼轨道转移的火星探测直接转移设计。北京航天航空大学与国防科技大学进行了载人火星方案的研究，提出了返回式火星探测的轨道方案。南京大学对火星探测器轨道变化特征做出了研究分析了不同摄动产生的影响。此外，装备指挥技术学院还对火星探测发射时机进行了研究，分析了逃逸速度渐近线和地球赤道面之间的夹角(DLA)对发射窗口的影响。1.1.3火星探测轨道优化探测器轨道设计

21、是基于发射任务、能量限制和测控范围等条件，并以轨道动力学为基础理论进行轨道设计，确定发射窗口和轨道参数（即标称轨道），然后经过轨道优化处理，从而得到设计轨道的过程。其中初步轨道设计是简化了探测器动力学模型的情况下得到的轨道参数，一般用来分析其特性和优缺点；精确轨道设计是由初步轨道设计的初值，并采用精确的探测器动力学模型，使用数值分析的方法求解轨道的过程。其是一个两点边值问题，即已知初始条件（转移轨道的初始速度和地球停泊轨道参数），选择合理的设计参数，优化出满足终端约束的若干条轨道。主要考虑终端约束条件和求解方法，其中终端参数用B平面参数表示时，搜索算法具有良好的收敛性，而求解主要包括基于目标函

22、数的优化方法和基于偏导数矩阵的微分修正。前一种方法在于如何选取到一种快速收敛的全局优化算法，目前应用在轨道优化的算法主要包括拟牛顿法、遗传算法等；后一种方法多次出现在NASA的设计任务中，而在国内最早是北京空间飞行器总体的杨维廉研究员在极月轨道设计中首次采用微分修正的方法，详细地给出了计算过程与步骤，并成功地应用在探月转移轨道中；自此之后，该方法引起了大量学者的关注，其中国防科技大学的高玉东提出了分层搜索的方法，讲述了搜索过程，将地月轨道设计分成瞄准搜索、到达搜索和精化搜索三个部分，能够快捷有效地显示出任意时刻探测器的状态，这种新思路可以应用在火星探测器轨道设计中。对于火星探测精确轨道设计，考

23、虑到探测距离远，采用基于偏导数矩阵的微分修正计算时间较长，虽然能够计算出精确轨道，但给导航和误差分析带来了严重的困难，而且收敛性不够好，一旦出现故障将带来巨大的损失，对于这种大型工程项目，因其技术复杂、投入风险大等原因，在任务初始阶段一般需要数字仿真，美国AGI公司开发的STK软件已经应用微分修正法来精确确定轨道，并且能够支持整个探测任务周期的全过程，包括需求、设计、制造、测试、发射、运行和应用环节。当火星探测器进入日心轨道后，由于存在各种误差（时间误差、导航误差、修正执行误差等）必将远离标称轨道，为使探测器准确到达火星目标点，必须进行中途修正，主要研究内容是修正时机的选取、误差分析和修正方法

24、三个部分。首先由于不同的假设条件与优化目标，修正时机的选取产生了众多优化理论，其中Battin提出的方差比率法忽视了能量最优问题，Pfeiffer提出的最小误差理论是从动态规划的角度去优化最小均方差，却忽视了推进剂的约束，而BreakWell提出的间距比理论是在满足终端约束条件下寻求一系列修正时机的能量最优理论。在实际工程上，各种测量手段所获取的测量数据必定存在着误差，而这样的数据经过定轨算法处理后，必定使日心轨道参数存在误差，如若误差小到一定程度，对到达火星终端参数影响不大，即在工程上可忽略，则可近似认为是探测器的真实轨道，如果误差不可忽略，那么研究误差产生的影响尤为必要。上世纪60年代Ki

25、zner提出了用B平面参数描述终端参数，并发现其与探测器轨道参数存在线性关系。处理含随机误差的线性系统特性的统计问题，在众多理论中，采用协方差分析描述法体现出高精度和省时的明显优势。之后大量学者用此法进行误差分析，这些方法都是在理论上进行分析，对实际工程具有一定的指导作用。鉴于B平面参数的优势，国外将其应用在中途修正中，这加速了深空探测轨道设计的发展，Carter详细介绍了采用该理论所涉及到的约束条件与相关技术，并给出了设计过程。国内对中途修正分析理论研究不足10年，都是从探月工程实施开始的，周文艳给出了月球的中途修正的数学模型，并研究了发射初始误差和修正时机的选取对修正速度脉冲的影响分析，合

26、理的设计了两次修正脉冲的时机和大小；针对平动点卫星探测任务，李明涛研究了Halo轨道的中途修正问题，分析了首次与末次修正时刻对中途修正的影响；在深空探测方面张晓文率先研究了自主中途修正，给出了一种相对简单的基于脉冲控制的自主修正方法，并以美国“凤凰号”火星探测器进行仿真验证；上海航天控制技术研究所的周杰给出了火星探测器到达火星影响球概率的计算方法，以能量最优为目的选取适宜的入轨瞄准点，在考虑各种误差的情况下设计了两套4次中途修正的方案，对今后我国自主开展火星任务具有参考意义。1.2 轨道基础知识1.2.1 时间系统时间是描述运动和事件的关键独立变量，其包含了计量的起点和步长两方面。对于探测器的

27、轨道描述和计算有着及其重要的意义。针对不同的探测器和任务，选用的时间系统的是不同的。本文的研究中，主要关心的问题是探测器在太阳系内的运动轨迹以及其与地球、目标行星的相对关系。主要涉及到的时间系统有以下几种:1)质心坐标时:在质心参考架中用来计算太阳系各行星及其卫星位置的独立时间变量。2)地球力学时:在地心参考架中的动力学时，当探测器从地球表面发射时，可将地球力学时作为探测器运动方程的时间变量。3)协调世界时:各国的民用时间标准，是观测资料所用的基本时间系统，也是用于深空探测器与地球通信的时间系统。本文使用了(NASA)美国航空航天局的JPL(喷气推进实验室)给出的行星精密历表，使用的时间为地球

28、动力学时(TDT)。1.2.2 坐标系统火星探测器的运行轨道主要由地心逃逸轨道、日心转移轨道和火星遭遇轨道，如图 1所示，主要涉及的坐标系有地心惯性坐标系、日心惯性坐标系和渐行线坐标系。图 1 火星探测器轨迹示意图地心惯性坐标系：也称赤道惯性坐标系，坐标原点在地球质心；轴沿着地球赤道面与地球黄道面的交线，即指向春分点；轴指向北极；轴与另外两轴构成右手坐标系。地心大地坐标系：该坐标系假设地球椭圆中心和短轴分别与地球质心和自转轴重合，经度为过目标点的地球椭球面和本初子午面之间夹角，纬度为经过目标点的地球椭球法线和地球椭球赤道面之间的夹角，高度为目标所在点至地球椭球赤道面的法向距离，如下图 2所示。

29、日心惯性坐标系：坐标原点在太阳质心；轴沿着太阳赤道面与太阳黄道面的交线；轴指向黄北极，轴与另外两轴构成右手坐标系。渐行线坐标系：坐标原点位于探测器质心；轴指向速度方向；轴指向轨道的动量矩方向；轴与其他两轴构成右手坐标系。则在渐行线坐标系下，探测器的速度分量为。J2000火心惯性坐标系：J2000火星惯性坐标系原点为火星中心，基本平面是历元J2000对应的火星平赤道，X轴指向J2000对应的火星平赤道与J2000对应地球平赤道的升交点，Z轴垂直于基本平面指向火星北极，Y轴与Z轴和X轴构成右手直角坐标系。图 2地球惯性坐标系与大地坐标系1.2.3 星历数据本文使用了美国喷气推进实验室(JPL)的行

30、星精密历表来获得各行星在给定时刻的位置速度。为了计算方便，调用JPL星历的时间统一为相对J2000.0历元的简约儒略口。输出的行星位置对应的坐标系为J2000口心赤道坐标系。1.2.4 B平面B平面是20世纪60年代初由Kizner.W发现的，主要思想基于目标天体的B平面上参数与探测器飞行轨道状态参数之间存在很好的线性关系。通常是以火星探测器的渐近线方向（即速度无穷远方向）为法线，并且过火心的假想平面，如图 3所示。火星探测任务中轨道的目标参数通常采用B平面坐标系中的B平面参数，其B平面坐标系的原点选在火星中心，通过火星中心并垂直于双曲线无穷远速度的平面称为B平面。记探测器进入轨道渐近线方向的

31、单位矢量为轴，取某参考方向的单位矢量为，其理论上方向是任意的，但一般选为火星赤道的法线方向，和的叉乘是轴，轴与轴和轴构成右手坐标系，即图 3 B平面示意图1.2.5 Lambert问题Lambert问题：航天器的初始与终了位置矢量和两者之间的机动时间t为已知的，从而确定始末速度矢量。这个问题的几何描述如图 4所示，它最早由拉格朗日和高斯从几何方面提出的，Lambert问题可以确定一系列轨道制导律和控制策略，从而这个基本问题吸引了大量的学者研究。图 4 Lambert问题求解Lambert问题有诸多方法，例如传统的高斯方法、p迭代法、级数法、普适变量法，甚至基于进化-模拟退火求解。上述方法各有优

32、缺点，传统的高斯方法仅针对小于90度的转移有效，而p迭代法需对圆锥曲线进行讨论，级数法对测量时间间隔有较为苛刻的要求，普适变量法可以适用于所有圆锥曲线轨道，但转移角为180度时发生奇异，而基于进化-模拟退火求解Lambert问题能够很好的解决上述问题，但计算量较大，过程较复杂。本文是采用一种无奇异的普适变量方法来求解Lambert问题。1.3 火星探测直接转移轨道的初步设计进行火星探测时，探测器从地球加速逃逸出发后，经过口心段转移到达火星，经过制动捕获过程到工作轨道后进行探测任务，途中不经过其它天体的甩摆。直接转移的探测器运行周期短，工作相对简单，适用于我国的首次火星探测。直接转移方法是指，探

33、测器在短时间内由大推力冲量的方法获得瞬时加速度，改变探测器的速度，通过一次大推力加速过程直接完成探测器飞往目标天体的转移。对于直接转移方法，可以分为初步设计和精确设计两个过程。对于初步设计，是以圆锥曲线拼接法(Patched Conic Method)为基本原理，将探测轨道划分为几个分段过程，每个分段过程可以近似成二体问题处理，然后通过接口，将每一段拼接在一起。它是在二体模型假设下，通过求解Lambert问题确定发射窗口，利用圆锥曲线拼接法确定地心段、日心段、火星段的初始轨道参数。1.3.1 日心轨道设计及发射窗口的搜索发射机会的搜索是火星探测任务设计及其规划的关键问题。项目成本与搜索时间是攻

34、关人员参考的重要指标。而在设计过程中，会受到诸多限制，例如发射能量、发射时间段、任务过程中的总速度增量等，这些因素不仅仅与上述指标相关，还导致系统的复杂程度和运载火箭的级别。故如何搜索到最优的发射机会是火星探测任务的首要问题。对于发射窗口的搜索，它是通过限定发射日期的区间以及飞行时间的区间，通过求解日心转移段的 Lambert 问题，确定性能指标最优的发射日期和飞行时间。根据优化求解方法的不同，主要可以分为两种方法，一种是基于枚举的思想，通过绘制时间与能量的pork-chop图，确定最优的发射时间和飞行时间的方法，而另一种是基于优化方法，通过建立飞行时间、转移能量的不等式约束以及目标函数，将问

35、题转化为优化问题进行求解。 基于pork-chop图的发射窗口搜索等高线图法是火星探测工程实施中搜索发射窗口的经典方法之一，通过给出既定时间段内所有的发射和到达时间情况，获得初始和终了位置，进而解算Lambert问题，从而绘制“猪排”图，观察得到最优发射机会。易知其直观性好，但计算量庞大，不便于进行大规模的发射机会搜索。基于等高线图的最优发射机会搜索算法的主要步骤可以总结如下：（1） 根据任务的需要确定出发射机会搜索的目标函数（性能指标）、发射时间的区间以及飞行时间的区间； （2） 选取一组出发时刻t0和达到时刻tf，根据行星历表计算地球的位置RE(t0)、速度VE(t0)和火星

36、的位置RM(tf)、速度VM(tf)；（3） 利用转移时间（tf - t0）、RE(t0)、以及RM(tf)，通过求解Lambert问题，可以得到探测器在始末位置处的速度矢量V1(t0), V2(tf)；（4） 确定发射机会的目标函数，并绘制出等高线图；（5） 根据等高线图，找到目标函数取值较小的区域，确定出性能指标指标最优的发射时间；以上步骤中涉及到的目标函数（性能指标）通常是指双曲线超速v，发射能量C3，速度增量v，它们的具体定义为： 双曲线在达到地球引力影响球边缘时速度有剩余，这个双曲线剩余速度v通常称为双曲线超速，计算公式为 式中v1是飞行器的速度矢量，vLE是发射时刻地球绕太阳公转的

37、速度矢量。 发射能量是影响任务初始设计的关键参数，在飞行器质量一定的情况下，发射能量越大，所需运载火箭的运载能力越强，它的大小是发射时双曲线超速的平方 当假设飞行器从停泊轨道开始转移，这在停泊轨道上施加的速度增量v1为 类似地，可以定义飞行器到达目标天体的双曲线超速以及到达目标停泊轨道，交会时的速度增量 式中v2是飞行器到达目标行星时的速度矢量，vAT是发射时刻目标行星绕太阳公转的速度矢量，rTp为飞行器在目标天体俘获是近心点的高度。在整个飞行过程中，所需的总的速度增量vtotal为 基于优化算法的发射窗口搜索基于pork-chop图的发射窗口搜索方法的直观性好，通过等高线图可

38、以清楚地看到给定时间段内发射机会的变化情况，然而，由于该方法的本质是一种穷举算法，该方法的确定是计算量庞大，不便于进行大规模的发射机会搜索。针对上述缺点，有些学者提出了基于优化算法的发射机会搜索方法，其选择总速度增量或发射能量等为目标函数，将搜索问题简化成搜索发射日期和飞行时间两个变量。这种方法在计算时间上有较大的进步，并解决了等高线图法的缺点，然而需多次迭代才能精确求解最优发射机会。在搜索发射窗口时，主要考虑的约束：（1）发射日期的范围；（2）转移时间；（3）发射后所需的速度增量；（4）所需总的速度增量；而优化目标也主要有：（1）飞行器逃逸地球时的双曲线超速v；（2）发射能量C3；（3）速度

39、增量v等；通过对搜索问题进行简化，可以将问题简化为寻求发射日期tL和飞行时间t，使得优化目标达到最优的问题，这里若定义T=tL,tT，则待优化的目标函数可以表示为 而约束条件可以简化为 1.3.2 地心段参数的确定 对地心段参数进行初步设计时，假设在地球影响球内探测器只受到地球的引力作用，从而按照二体轨道特性进行轨道参数的计算，逃逸速度增量假设为脉冲。假设探测器的发射过程为:探测器从地面发射后，首先进入圆形的停泊轨道;然后在经历一段时间的滑行后，在特定的时间，特定的位置，经过末级火箭加速后进入双曲线轨道进行逃逸，如图 5所示。为了尽量避免不必要的燃料消耗，假定停泊轨道和双曲线轨道是共面的，且双

40、曲线的近地点半径与圆轨道半径相同，速度相切。因此，地心段参数初步设计的内容主要包括两个方面:停泊轨道参数的特性分析以及停泊轨道参数选择。图 5地心逃逸双曲线轨道面内参数的几何关系设停泊轨道半径(即双曲线近地点半径)为今1，根据双曲线轨道参数的相关公式容易算得双曲线轨道半长轴、偏心率、动量矩以及近地点速度: 其中a1为轨道半长轴，e1为偏心率，为地心引力常量，v1为地心双曲线剩余速度大小，hl为双曲线轨道的动量矩大小，vp1为双曲线近地点速度大小。 而根据剩余速度矢量v1可以计算出转移双曲线渐近线对应的赤经1（RLA: Right ascension）和赤纬1（DLA: Declination）

41、 进一步可以确定B平面中的渐近线的单位方向矢量 而B平面的其他两个坐标系可以表示为 式中，B平面的夹角可以由渐近线的赤纬和轨道倾角确定 双曲线的单位角动量矢量，可以表示为 当速度趋于无穷时，真近点角的正弦和余弦可以确定为 进一步可以确定转移双曲线轨道的近地点矢量 从而可以确定航天器在双曲线转移轨道的近地点处的位置矢量为，而速度的单位矢量为 而航天器在近地点的速度矢量可以表示为。根据近地点的位置矢量和速度矢量，就可以计算出双曲线轨道的轨道六根数。考虑到行星际距离比地球影响球大很多，故地心段双曲线轨道与地球影响球的交点位置对轨道几乎没有影响，而v1矢量的大小和方向完全决定了探测器在脱离地球作用范围

42、之外的运动，因此控制v1矢量大小和方向的精度非常重要。实际中存在无数条满足剩余速度大小方向约束的逃逸双曲线存在，这些双曲线位于将图 5中双曲线轨道绕过地心的v1矢量(即渐近线的方向)旋转一周组成的曲面上。并且所有双曲线的近地点组成的圆形轨迹称为轨道发射圆，因此停泊轨道必然经过OL点，等探测器运行到轨道发射圆的位置，沿着速度方向施加脉冲，进入双曲线轨道，如图 6所示。图 6轨道发射圆示意图设v1矢量在地心惯性下可以表示为: 则地心到OL点的矢量的赤经为（+1），赤纬为（-1），表示为： 由以上分析可知，停泊轨道的轨道倾角i必须满足如下关系式 根据轨道倾角i，发射场的纬度L以及发射方位角AZ三者之

43、间满足以下关系式 由以上关系分析可知，轨道倾角还应该满足如下关系式 假设在北半球某发射场的发射场的纬度为L，可以在发射方位角满足范围内任意发射。发射方位角的范围为AZ1<AZ<AZ2。图 7发射方位角和发射窗口示意图考虑地球自转的影响，发射场从经度为1，的位置开始，到（+1）的位置，在这半天的时间内，发射窗口如图 7所示。发射场从点1随地球自转到点2对应的时段就是发射窗口。其中，在点1对应的发射方位角为最小方位角，点2对应的发射方位角为最大方位角。在发射窗口内，每个时刻对应一条相应倾角的发射轨道。一般情况下，在这些窗口发射的探测器，可以在停泊轨道上停留的最小时间少于停泊轨道周期的一

44、半。而发射场从经度为（+1）的位置转到经度为1，位置的过程中同样存在类似的窗口，区别是在这些窗口发射的探测器在停泊轨道上停留的最小时间大于停泊轨道周期的一半。对于选定的发射场，实际工程中为利用地球自转能量，测控站点分布和安全考虑，发射的轨道倾角要满足一定的要求。深空任务本身对于地心段发射的轨道倾角没有特定要求，而只需考虑测控等因素，故而一般选取该发射场能发射的最小轨道倾角的顺行轨道。根据轨道倾角与v1的赤纬的约束，如果v1的赤纬小于发射场所在的纬度，则轨道倾角可选择为发射场可发射的最小倾角，如果v1的赤纬大于发射场所在的纬度，则需要发射的最小轨道倾角等于v1的赤纬。所以，在初步设计中，轨道倾角

45、取值为: 选定轨道倾角后，有两组值可以满足v1矢量的要求，设这两组双曲线轨道升交点赤经，近地点幅角分别为，则其满足如下公式: 其中 第一组参数表示双曲线的近地点速度的Z分量小于0，属于下降段入轨。发射场的位置在v1矢量沿着N轴顺时针旋转的角度小于180o，探测器可以在发射后小于半圈的停泊入轨。而第二组参数表示双曲线近地点的速度Z分量大于0，属于上升段。发射场的位置在v1矢量沿着N轴顺时针旋转的角度大于1800，探测器需要在发射后大于半圈的停泊后入轨。1.3.3 火心段参数的确定对火心段参数进行初步设计时，可以认为在火星影响球内探测器只受到火星的中心引力的作用，从而按照二体轨道特性进行轨道参数的

46、计算，捕获速度增量假设为脉冲。火心段双曲线轨道的特性与地心段双曲线的性质完全类似。在地心段双曲线轨道的剩余速度约束是设计的目标，而在火星段双曲线剩余速度的约束是初始条件，其大小和方向是固定的。火星段轨道参数的选择主要内容是确定双曲线的近火点半径和轨道倾角。假设进入双曲线的轨道和探测器的工作轨道是共面的。为了实现火星探测器对于火星的全面覆盖，火星探测的轨道选取为极地轨道，倾角约为i2=90o。双曲线剩余速度的赤纬为，故而轨道倾角的取值是满足式。假设双曲线的近火点高度为rp2。在进行轨道设计时，通常将目标的双曲线的参数用B平面参数来描述。探测器进入火星影响球时，可以根据火心双曲线进入速度v2与目标

47、火心轨道的倾角和近心距rp2计算标称B平面参数BT和BR。 其中为火星引力常量，为B矢量与T矢量的夹角，如图 8所示，i2为目标火心轨道的轨道倾角，速度v2为火心双曲线进入速度大小，2为该矢量的火星赤纬。图 8火心进入双曲线轨道与B平面示意图1.4 基于B平面参数的精确轨道设计在实际工程中，各种测量手段所获取的数据是存在误差（初始状态误差），带误差的数据经过处理后，所获得的轨道必然存在偏差。所以火星探测器的轨道误差分析极具有意义，在火星探测器发射后，可以全面地评估各种误差对任务的影响。1.4.1 问题描述探测器在各个天体(本文主要考虑地球、火星和太阳)的引力和众多摄动力的影响下，其动力学模型如

48、下其中是其他摄动力加速度之和。在探测任务中，火星探测器从地球逃逸轨道到与火星相遇为止，从燃料消耗的角度来看，火星探测中途修正不是把轨道修正成标称轨道，而是在有误差的位置上施加一个合适的控制参数增量，从而使火星探测器沿着一条新的轨道机动来满足对终了状态的要求。火星探测器的终端参数一般选择为目标轨道的倾角、近心距和B平面参数。选择好终端参数，记初始时刻为（选变量为位置和速度）。火星探测器抵达目标区间的终端参数记做，则火星探测器初始状态和终了状态之间可以用某个函数来表示，即微分修正问题就是采用一定的制导方法使探测器抵达目标区间的实际状态与期望状态的误差小于规定值1.4.2 制导方法我们将实际轨道在标

49、称轨道附近进行泰勒展开后只保留线性项，得其中，是被控制量，敏感矩阵。求解式(3.39)的方法主要考虑和的维数，其分别为和1）此种情况较容易，只需对敏感矩阵求逆便可求得2）此种情况控制量的数量多，在式的约束下，求得幅值最小的修正量即可。其性能指标为其中是Lagrange乘子。可利用变分法求解式，即 3）此种情况控制量较少，可采用最小二乘法求解该问题。即求式的最小值其解是，但此解不能保证小于规定值。1.4.3 轨道精确设计求解轨道精确求解是解两点边值问题，从前面的内容可知，即约束自由的多圈Lambert问题的时间自由下的最优双脉冲解，通过上一章的计算大致知道时间的具体范围，而且考虑能量问题，圈数只

50、可能为0，故本文不采用解析的方法求得，而是采用微分修正的方法，并给出一种数值求解式中的偏导数矩阵的方法。下面是以B平面参数为终端参数并采用微分修正求解精确轨道设计的步骤：基于B平面参数的火星探测器的精确轨道设计方法，其具体过程如下：步骤一、通过遗传算法确定的火星探测器的轨道参数计算初值，主要包括地球逃逸轨道参数和日心转移轨道参数，选取地球逃逸轨道的轨道倾角、轨道半径、轨道初速、升焦点赤经和近地点幅角为控制参数。步骤二、根据控制参数为初值在精确动力学模型下进行轨道数值积分运算，求得终端参数值，终端参数选取B平面参数（、）和到达时间。步骤三、计算得到的参数值与标准参数进行比较，获得参数偏差量，从而

51、求得新的控制参数。步骤四、利用新的控制参数重新对动力学模型进行轨道积分运算，得到新的终端参数值偏差，这个偏差量逐渐减小。步骤五、重复上述过程，直到终端参数满足精度要求。其程序流程图如图 9所示：图 9基于微分修正的精确轨道设计流程图设火星探测器在近地点和近火点的速度与位置构成的状态矩阵分别为，若出发时刻和飞行时间确定，则、及存在的关系为其中、均为多元函数组，由多元函数法则，考虑一阶项，得即其中、均为一阶偏导数构成的矩阵，具体表达形式如附录一所示则相对的偏导数矩阵1.4 仿真分析仿真任务：在未来3年内寻找发射机会，设计地球火星转移轨道，以总能量最小为指标。假设地球停泊轨道高度200km、轨道倾角为的圆轨道，目标轨道为高度500km、倾角为的火星绕飞轨道。轨道动力学模型只考虑太阳、地球、火星的引力作用及地球J2项摄动。分别给出初步轨道设计参数和精确轨道设计参数。1.4.1初步轨道参数设计结果地球和火星的会合周期为约为780天，为了寻找未来三年内的发射机会，设