论文初稿晚稿.docVIP

参赛密码 （由组委会填写） 全 第九届“华为杯” 全国研究生数学建模竞赛 学 校 华北电力大学（保定） 参赛队号 队员姓名 高成彬 李少岩 刘贺晨 参赛密码 （由组委会填写） 第九届“华为杯”全国研究生数学建模竞赛 题 目 基于遗传算法的无源探测 空间飞行器主动段轨道估计与误差分析 摘 要： 关键词：遗传算法 、参数估计、PID、误差分析、轨道预测 一、问题重述 有些国家会发射特殊目的的空间飞行器，如弹道式导弹、侦察卫星等。对他国发射具有敌意的空间飞行器实施监控并作出快速反应，对于维护国家安全具有重要的战略意义。卫星居高临下，是当今探测空间飞行器发射与轨道参数的重要平台。装置于卫星上无源探测器被动地接收目标辐射，其目的是为了推断空间飞行器的轨道参数 在无源探测中，多使用红外光学探测器，单个红外光学探测器不具备测距能力，但借助多颗（含两颗）观测卫星的同步观测能够进行逐点定位，再结合空间飞行器的运动模型，可以进行轨道参数估计。在单星观测条件下，利用空间飞行器轨道的特殊性，结合较强的模型约束也可得到一定精度轨道参数估计。由于受大气影响，垂直上升段的火箭尾焰不易观测，程序拐弯段的运动方程又较为复杂，所以本题重点关注重力斜飞段的后程段，本题所附仿真数据也集中于此段。 本题以中低轨近圆轨道卫星为观测星座对假想的空间飞行器进行仿真观测，生成仿真观测数据，要求利用仿真观测数据，对假想空间飞行器的轨道参数进行估计。 我们将研究下列问题： 观测卫星在任意时刻的位置计算是估计的前提，根据satinfo.txt和观测卫星的简化运动方程（2），计算09号观测卫星在50.0s、100.0s、150.0s、200.0s、250.0s五个时刻的三维位置。结果保留6位有效数字。 给定的仿真数据下，06号和09号观测卫星对0号空间飞行器形成了立体交叠观测，我们将结合立体几何知识按照逐点交汇定位的思路，给出0号空间飞行器在公式（1）框架下的轨道估计，并选取适当的和的表示模型。按照从50.0s到170.0s间隔10.0s进行采样，计算并列表给出0号空间飞行器在各个采样点的位置和速度，并给出估计残差。结果保留6位有效数字。同时绘制0号空间飞行器的三个位置t-x、t-y、t-z和三个速度t-vx、t-vy、t-vz曲线示意图。 在同时考虑系统误差的条件下，进一步研究下列问题： 若06和09号两颗观测卫星均有可能带有一定的系统误差，对系统误差进行正确的估计能够有效提高精度。我们将考虑利用上述的逐点交汇方法能否同时对系统误差进行估计？若不能，我们将寻找其他的思路能够同时估计系统误差与轨道。并给出与第二问要求相同的结果，以及两颗观测卫星的系统误差估计结果，共六个数值，分别是两颗卫星的。 我们对只有09号观测卫星单星观测的01号空间飞行器进行轨道估计并进一步考虑在同时有多颗观测卫星观测多个空间飞行器的情况下能否联合进行系统误差估计？ 二、模型假设 本题假设为形式简单的严格单调递减的非负函数。的方向与飞行器的速度相反方向接近或相同，其大小较为稳定。 不考虑飞行器上各点的速度差异； 不考虑飞行器质心因燃料消耗而产生的位移； 作用于飞行器上的所有外力都通过质心，也就是说不计导弹在飞行过程中绕其质心的转动。 本题假设观测坐标系与基础坐标系之间的转化误差很小，可以忽略不计； 我们仅考虑单级火箭模型，对于多级火箭的模型不予考虑，以简化模型。 三、符号说明 ：基础坐标系 ：观测坐标系 ：空间飞行器在基础坐标系下的位置矢量 ：Vr(t)是燃料相对于火箭尾部喷口的喷射速度 ：飞行器瞬时质量 ：地球引力常数（本题中地球引力常数取） ：飞行器所受的外力加速度之和 ：火箭产生的推力加速度 ：观测坐标系下的两个无量纲比值 由于本文符号较多，其他符号在论文中第一次出现时，再进行说明。 四、问题分析： 1、问题1： 文件中已经给定基础坐标系下的位置和速度初值，可以利用常微分方程数值解方法计算空间飞行器的运动轨迹。09号观测卫星的运动方程满足题中形式较为简单的微分方程。此外，把09号观测卫星在基础坐标系中的位置和速度作为微分方程的初始条件，通过编程即可得到任意时刻的09号卫星的位置的数值。 2、问题2： 已知06号和09号观测卫星对0号空间飞行器的进行了立体交叠观测，形如meadata_i_j.txt的数据文件中存放着编号为i的卫星对编号为j的飞行器的仿真观测数据信息