基于位置观测的北斗广播星历速度精度分析

利用广播星历计算导航卫星的速度向量是GNSS高精度实时测速的必要条件.本文分析了仅以卫星位置向量为观测量的北斗广播星历的速度计算精度.从广播星历拟合过程出发,推导了北斗18参数模型的速度向量计算公式.基于北斗13颗在轨卫星一年的实际轨道数据,分析了全年广播星历计算卫星速度向量的精度.结果表明,利用18参数模型计算的速度误差最大在10-4m·-1量级;在相同拟合时段条件下,地球静止轨道（GEO）和倾斜地球同步轨道（IGSO）卫星的速度精度相当,高于中圆地球轨道（MEO）卫星.通过对位置残差序列分析,得出位置残差误差较小且变化趋势平稳是广播星历计算速度精度较高的原因.分析和计算结果验证了仅用位置观测量拟合北斗广播星历算法的有效性.      

基于状态转移矩阵的航天器多脉冲悬停方法

基于航天器相对运动的状态转移矩阵描述,研究了空间相对悬停的多脉冲控制方法,解决了工程实践中连续推力悬停轨道控制技术对航天器控制推进系统要求较高的难题。给出了两航天器在圆、椭圆和双曲线等圆锥曲线参考轨道上相对运动的状态转移矩阵描述。在此基础上,定性分析了椭圆参考轨道偏心率对悬停精度的影响,推导了航天器多脉冲悬停速度脉冲控制量的计算方法。数值仿真算例显示,该方法可有效实现一定悬停精度要求下的空间相对悬停控制,且随着一个轨道周期内脉冲数的增加,相对悬停的效果得到提升。     

时间与燃料约束的参数自主寻优变轨滑模控制

航天器相对运动轨控采用滑模控制具有较好的抗扰能力，但参数设置复杂。为贴近工程实际，引入燃料最优约束和寻优算法，提出一种综合考虑时间、燃耗以及误差的参数自主寻优滑模控制。首先，基于线性相对运动方程与指数趋近的滑模控制，建立相对运动滑模控制器模型，并由能量最优的轨迹规划器给出收敛约束时间，实现高效机动；然后，分析滑模控制器中可调参数与时间、误差的约束条件，制定了参数量级寻优规则；最后，通过惯性权值改进的粒子群算法，将误差允许范围内的最少燃料消耗作为寻优评价标准，输出最优量级与系数组合的控制参数，实现滑模的最优控制。仿真表明，使用粒子寻优器得到的参数组合，可使滑模偏差控制器在规定时间内通过最小燃料消耗令位置与速度误差稳定收敛，增加航天器在轨寿命。     

有限推力下的航天器绕飞轨道保持与控制

基于线性动力学模型,研究了有限推力下的航天器绕飞轨道保持与控制问题.首先,针对绕飞轨道控制问题推导了精确的动力学模型,给出了经过线性化的动力学方程.然后,在线性二次性最优控制理论的基础上,设计了一种针对绕飞轨道的反馈控制系统,同时为了实现有限推力,设计了一个有界的非线性环节对控制加速度进行限制.数值仿真结果证明了所设计的控制系统渐进稳定、收敛快速,并且具有良好的控制精度.最后,通过计算不同相对位置误差时轨道控制系统需用的发动机工作时间,估算了轨道维持与轨道控制的燃料消耗量.      

连续常值推力机动分析与应用

连续常值推力机动是空间飞行常用的轨道机动方式。其中,小推力适合于地球轨道航天器交会机动,而切向或周向推力以及较大的正径向推力可用于脱离地球引力场的逃逸飞行,执行星际交会使命。应用常推力作用下的质心运动方程,对机动推力的量值没有限制;在航天器交会应用中,对相对距离也无要求。这种方法可直接获得向径、轨道速度等参数随时间或极角(绕地心的转动角)的变化,便于分析轨道转移与逃逸运动,有助于飞行使命与运动轨迹的设计。特别是,若机动转移的初轨为圆轨道,在推力较小、飞行时间不长的情况下,应用量纲1形式的运动方程,可获得具有工程应用价值的近似解。     

太阳帆航天器编队维持和重构的方法研究

针对高面质比航天器可以利用太阳光压进行轨道控制的特点,本文提出一种太阳帆航天器编队构型维持和重构的方法.该方法通过控制主从航天器太阳帆姿态角和反射系数,调整主从航天器之间的光压差,产生抵消编队成员间相对运动受到摄动差或进行轨道机动时所需的连续小推力,从而实现编队构型的维持和重构.仿真结果表明,在主航天器太阳帆的姿态角和反射系数相对固定的条件下,对于太阳同步轨道上的高面质比太阳帆航天器编队,使用滑模控制方法,能够调整编队中从航天器太阳帆的姿态角和反射系数产生推力抵消摄动力影响,达到长期维持太阳帆航天器编队构型的目的;通过开环控制方法,能够调整编队中从航天器太阳帆的姿态角和反射系数产生连续小推力,在较长时间周期内实现编队重构.     

星载GNC技术现状及发展浅析

<正>星载制导、导航与控制(GNC)技术是指为了实现特定目标而对航天器进行的一系列的轨道规划、轨道机动和姿轨联合控制。星载GNC系统是卫星系统的重要组成部分,是进行任何航天任务必不可少的先决条件。其中,导航定义为测定航天器的位置和速度,并且计算出达到期望位置和速度的相关调整量;制导与控制定义为通过对航天器进行控制,以保持航天器满足所需的指向精度要求,同时,在必要时,还要进行跟踪导航计算。     

高精度在轨实时轨道机动决策

为保证在轨机动实时性和高精度的要求，提出了一种基于机器学习的在轨实时机动决策方法。通过优化算法离线获得摄动下的精确解，减去二体解得到速度增量差，将其投影到轨道坐标系获得速度增量摄动修正项，以此作为神经网络输出，设计网络参数并训练得到摄动修正网络、组合应用摄动修正网络和二体解实现高精度的在轨实时轨道机动决策。仿真结果表明:卫星按照该决策机动完成后的终端位置偏差与按照优化算法给出的决策机动完成后终端位置偏差精度一致，且前者决策耗时仅为后者决策耗时的0.01%左右。所提轨道机动决策方法兼顾了精度与实时性，适用于星上决策。      

静止轨道卫星东西位置保持控制参数的优化方法

为提高静止轨道卫星东西位置保持的轨道控制精度,提出一种优化轨道控制参数和推进剂消耗量估算的方法.这种方法考虑卫星轨道机动期间其姿态控制喷气对卫星轨道的摄动力,据此对位置保持参数中的轨控推力器点火时间进行补偿,并根据所有推力器的喷气时间估算推进剂消耗量,以修正轨控策略,实现更加精准的轨道控制效果.以在轨运行的静止轨道卫星的多次向西位置保持的变轨参数优化过程为实例,验证本方法的有效性和正确性,可应用于后续静止轨道卫星东西位置的实际操作中.     

基于经验加速度的低轨卫星轨道预报新方法

研究将定轨过程中的经验加速度应用于地球低轨卫星轨道预报的新方法. 利用GPS伪距观测数据和简化动力学最小二乘批处理方法对地球低轨卫星定 轨, 其中卫星位置、速度及大气阻力系数和辐射光压系数可以直接用于轨道预报. 作为简化动力学最重要特征的经验加速度呈现准周期、余弦曲线特点, 可通过 傅里叶级数拟合建模. 确定性动力学模型与补偿大气阻力模型误差的切向经验 加速度级数拟合模型组成增强型动力学模型用于提高轨道预报精度. 应用 GRACE-A星载GPS伪距观测数据和IGS超快星历定轨并进行轨道预报, 结果表明 轨道预报初值位置精度达到0.2m, 速度精度达到1×10-4m·s-1, 预报3天位置精度优于60m, 比只利用确定性动力学模型进行预报精度平 均提高2.3倍. 先定轨后预报的模式可用在星上自主精确导航系统中.      

基于线性协方差方法的交会对接误差分析

将线性协方差分析方法和蒙特卡罗仿真相结合,按交会任务和飞行特征把交会过程分为变轨飞行、自由飞行和中途速度修正三种特征段,研究了状态误差的传播规律和交会过程中各种误差对交会对接精度的影响。在变轨飞行段,分析了追踪航天器的姿态误差、控制系统性能状态估计误差,以及目标航天器轨道摄动对状态误差传播的影响。在自由飞行段,分析了追踪航天器估计状态误差的先验值和测轨误差对状态误差传播的影响。在中途速度修正段,分析了追踪航天器姿态误差和控制系统性能误差对状态误差传播的影响。仿真结果表明,误差分析方法设计合理,可以指导交会对接的轨道设计工作,能对已经设计好的交会策略进行误差分析和设计验证。     

空间激光干涉引力波探测器轨道修正方法

针对空间激光干涉引力波探测器轨道修正问题，提出一种基于虚拟编队构型设计的航天器轨道修正方法。空间激光干涉引力波探测器由3颗航天器组成等边三角形构型。由于入轨误差和摄动的影响，探测器的构型不稳定。假设名义轨道上运行着一颗理想航天器，实际轨道上的真实航天器与之组成虚拟编队，探测器的3颗真实航天器分别与对应的理想航天器组成3个虚拟编队。考虑探测器构型稳定性要求和摄动的影响，对虚拟编队的构型进行设计，进而求解航天器平均轨道要素修正量。求解得到的航天器平均轨道要素修正量小于偏差量，轨道修正通过四脉冲控制实现。数值仿真结果表明，该方法通过部分轨道修正满足了探测器的构型稳定性要求，具有减少燃料消耗、延长任务寿命的潜力。      

LSTM神经网络在太阳F10.7射电流量中期预报中的应用

F_10.7指数作为大气密度经验模型的重要输入参量，其预报精度直接影响航天器轨道预报精度.研究发现，太阳活动表现出长时间尺度上平均11年和中短时间尺度平均27天的周期性变化特征.依据这一观测事实，基于长短期记忆单元（Long Short-term Memory，LSTM）递归神经网络方法进行F_10.7指数未来27天的中期预报.利用一个连续长时段F_10.7数据作为训练数据，构建LSTM神经网络训练和预测模型，分别预测太阳活动高低年未来27天的F_10.7指数.结果表明，太阳活动高年的第27天F_10.7指数预报平均相对误差最优可达10%以内，低年最优可达2%以内.      

TLE预报精度改进及碰撞预警中的应用

双行根数（Two Line Elements,TLE）是目前唯一公开发布且编目最完备的地球轨道空间目标编目数据,由于其根数的特殊性,必须配合SGP4/SDP4 使用导致其预报精度有限. 通过对TLE+SGP4/SDP4预报误差源的分析,利用历史TLE数据生成准观测数据重新进行轨道拟合,针对不同轨道高度和面积质量比情况,确定相应拟合周期,得出了基于弹道系数先验信息有效确定面积质量比的经验方法,对多个独立时段拟合后轨道根数结合数值方法的预报误差进行统计,结果表明,拟合后的TLE根数预报精度与稳定性有很大程度提升. 利用实例对拟合TLE在碰撞预警中的应用进行分析,认为该方法对航天工程中碰撞预警置信度的提高具有重要意义.      

返回式航天器轨道周期变率确定及预报的数值方法

对返回式航天器进行变轨控制计算时要用到平均轨道周期变率,而数值法定轨通常求解的是大气阻力系数,无法直接得到平均轨道周期变率。文章通过建立航天器精密动力学模型、数值积分器、瞬时轨道根数到平均轨道根数的转换算法和平均周期序列多项式拟合算法,提出了一种基于数值法精密轨道确定和预报计算平均轨道周期变率的数值方法。     

基于GNSS测量的天宫二号质心确定

由于轨道机动燃料消耗，科学载荷加载、分离，以及伴飞小卫星在轨释放等原因引起天宫二号空间站质心(COM)发生位移，从而影响天宫二号的动力学质心定轨精度。针对这一问题，提出了基于全球导航卫星系统(GNSS)测量数据的简化动力学质心估计方法。燃料消耗是引起天宫二号质心发生位移的主要原因，质心在本体坐标系X轴方向位移最为显著。利用GNSS测量数据对天宫二号进行质心估计和精密定轨，在三轴对地稳定姿态下，本体坐标系X轴方向与轨道切向重合，定轨结果对本体坐标系X轴方向的质心位移并不敏感。但在连续偏航模式下，本体坐标系X轴在轨道法向上有较大分量，X轴方向的质心位移对基于GNSS测量计算的精密定轨结果有较大影响。定性和定量分析结果表明:偏航姿态模式下天宫二号本体坐标系X轴方向质心位移估计具有可行性。天宫二号实测数据计算结果表明:与未做质心估计的定轨结果进行对比，质心估计后表征轨道动力学建模误差的经验加速度补偿水平在轨道径向、切向和法向上分别降低62%、50%和65%；载波相位后验残差标准差降低0.04 cm；精密轨道与全球激光测距数据比较精度提高0.86 cm。所提方法可以应用于大型低轨航天器在轨质心估计。      

基于电推进卫星飞往Halo轨道的转移轨道研究

利用电推进及轨道力学的特性实现节能优化,将限制性三体问题中的稳定不变流形与小推力轨道优化相结合,研究全电推进卫星从地球停泊轨道飞向日地拉格朗日L2点Halo轨道的低消耗转移轨道.航天器的转移轨道分为逃逸段、拼接段与无动力滑行段.在逃逸段卫星沿速度方向加速脱离地球引力,拼接段采用Radau伪谱法进行优化,使航天器以最短时间到达目标Halo轨道的稳定不变流形上,随后航天器电推进系统关机,沿稳定不变流形无动力滑行至目标轨道.基于雅克比积分常数给出拼接段轨道初始猜测值,以先提高切向方向航天器能量避免了全程优化离散点过多难以求解的问题.仿真结果表明,该方法收敛速度较快,对平动点工程任务的初期轨道特性计算具有实际意义.      

基于CO2测量数据的大气辐射传输模型LBLRTM优化

根据TIMED/SABER 2002—2018年的CO₂观测数据，分析CO₂浓度的变化特征.依据变化特征给出了CO₂浓度随时间、高度、纬度变化的月平均拟合公式，利用非线性最小二乘拟合法，对不同高度和不同纬度的CO₂浓度数据分别进行拟合，生成相应的拟合参数.然后，将所有拟合参数汇总并生成拟合参数文件，结合拟合公式构建全球CO₂浓度经验计算模块，并将该模块应用到大气辐射传输模型LBLRTM中，对该模型进行优化.将优化前与优化后的LBLRTM模型模拟结果分别与TIMED/SABER观测数据进行比较发现，优化前的LBLRTM模型模拟结果与观测值的均方根误差为15.4%，而优化后的LBLRTM模型模拟结果与观测值的均方根误差由15.4%下降至8.91%.结果表明该优化方法可以提高LBLRTM模型在红外波段的辐射模拟精度.      

基于多项式近似的空间碎片群体轨道预报算法

快速准确地分析空间碎片群轨道演化行为对于其他在轨航天器碰撞规避至关重要。在各摄动力的作用下,空间碎片群演化运动呈现出复杂的非线性特征。空间碎片群体个体数量巨大,如果通过对空间碎片群中每个空间碎片进行轨道积分来分析群体预报的方法会导致计算量过大。针对该问题,提出一种基于多项式近似的轨道快速预报分析方法。该方法将空间碎片群分为少量的标称碎片和其他大量关联碎片。针对标称碎片的轨道预报采用数值积分求解保证预报精度;而针对其他大量的关联碎片轨道预报问题,采用多项式泰勒展开半解析方法求解,从而在保证预报精度的前提下有效减少空间碎片群轨道预报的计算量。为了验证方法的有效性,对不同空间碎片群进行了轨道预报仿真。仿真结果表明,当轨道预报精度设定在1m范围内时,多项式近似算法的计算量较蒙特卡洛方法计算效率提高了2.2~17.2倍,验证了所提出方法的有效性。