一种有效提高超大偏心率轨道初轨确定精度的新方法

针对某型号航天器的运行轨道偏心率大、轨控结束后可用于初轨确定的测轨数据少的特点,提出了一种能有效提高超大偏心率短弧段初轨确定精度的新方法。该新方法将轨道改进的思想融入初轨计算方法中,使之既适用于各种偏心率轨道的初轨计算,也适用于各种偏心率轨道的长弧段、多圈、多站数据的轨道改进;解决了传统初轨确定方法超大偏心率轨道短弧段初轨确定误差偏大的问题。通过任务实测数据的检验,新方法适用于各种偏心率轨道,并且超大偏心率轨道的初轨确定精度明显优于传统方法。     

编队小卫星星间基线估计算法研究

H ill方程假设卫星轨道是小偏心率和无摄动的,实际应用时难以获得高精度的星间基线估计。本文提出了一种基于主星位置展开的扩展Kalm an滤波算法,该算法采用含摄动的卫星运动模型,综合利用星间无线电和红外仪器的原始量测,实现了小卫星星间基线的融合估计,提高了估计精度。仿真实验针对J2摄动模型验证了算法的可行性和有效性。     

一种数值定轨方案

本文给出一种有别于习惯采用的人造卫星精密定轨方法，该方法仍以分析定轨法的轨道根数作为待估状态量，但却采用数值法计算精度要求较高的卫星受摄星历，而用分析法去计算条件方程中所需的状态转移矩阵。这样既保持了清晰的轨道图象，又使计算卫星受摄星历和状态转移矩阵得以简化，相应的软件形式简单。只要定轨弧段Ｓ不太长（Ｓ〈１０＾３），就能在一定精度条件下明显提高卫星定轨的计算效率。     

一种低轨卫星轨道预报算法

在基于EMBET（Error Model Best Estimation of Trajectory，误差模型最佳弹道估计）的雷达自校准处理过程中，需要基于给定的低轨卫星初始状态完成轨道预报，得到与雷达测量序列对应的卫星位置序列，进而通过迭代收敛过程得到误差模型参数估计。针对低轨卫星轨道迭代预报的需要，提出了一种边积分边插值的轨道预报方法。即将低轨卫星轨道积分计算过程的加速度计算结果直接应用到三弯矩插值，省略了三弯矩插值中节点矩的计算过程，同时将积分步长以及插值区间步长设置为1 s，减少了积分及插值中乘法和除法的计算次数，进一步提高了计算效率。经仿真计算表明，该方法在保证轨道预报精度的情况下，计算效率比传统低轨卫星轨道预报方法提高了约20倍，有效提高了雷达自校准处理的时效性。     

一种卫星精密星历的插值方法

在卫星各种应用中通常需要对卫星的星历进行插值,而且要求插值后星历精度仍保持与原轨道精度相当。本文采用在拉格朗日多项式插值方法基础上改进的Neville算法,对一颗低轨道地球卫星进行了仿真计算,分别对采用轨道根数和直角坐标速度形式的轨道进行插值。仿真结果表明选取合理的插值方法后,可以得到高精度的卫星星历,具有工程应用价值。     

大型航天器无控飞行再入时间短期预报的轨道摄动方法研究

针对寿命末期大型航天器临近再入前轨道衰降受到跨流域多尺度非平衡环境下气动力/力矩影响显著的问题,基于体坐标系下的气动力和力矩分量,建立了航天器轨道-姿态动力学摄动模型;结合航天器再入跨流域空气动力学统一计算理论与数值模拟手段,提出了大型航天器无控飞行轨道衰降过程气动特性一体化快速算法。基于地心惯性坐标系气动融合轨道动力学方程数值积分方法对航天器的轨道姿态摄动模型进行了外推计算,结合无迹卡尔曼滤波方法对外测轨道星历观测数据初值误差和模型误差进行了滤波。针对250～120 km高度范围的大型航天器短期轨道衰降变化进行仿真预测,对比了外测轨道星历数据,结果表明:本文提出的经卡尔曼滤波误差分析的轨道摄动模型方法对寿命末期大型航天器短期轨道预报具有较强的模拟能力与有效性,有助于推动近地轨道空间目标气动融合轨道数值预报领域发展。     

基于北斗混合星座的星历参数适用性分析

针对北斗全球卫星导航系统混合星座条件下不同星历参数表达方法的性能评估问题,采用了从信息占用量、轨道拟合精度、接收机首次定位时间3个方面统筹考虑的评价方法。针对北斗系统的混合星座构型,比较北斗广播星历16参数和新设计的广播星历18参数的参数定义及用户算法异同,分析了不同星历参数表达所带来的性能变化。仿真分析结果表明,两种星历参数表达虽然都能描述北斗混合星座中的地球同步轨道(GEO)、中间地球轨道(MEO)、倾斜地球同步卫星轨道(IGSO)的卫星轨道特性,但是它们的具体性能有明显的差异,星历18参数的轨道拟合精度比16参数的提高了2~3cm,但其信息表达却多占用了63bit,且接收机首次定位时间多消耗了1.26s。     

北斗GEO卫星实时求速算法研究

以北斗系统空间信号接口控制文件中规定的卫星位置计算标准算法为基础,研究了采用北斗广播星历实时计算北斗GEO(Geostationary Earth Orbit,地球静止轨道)卫星速度的算法。首先,通过将开普勒参数和摄动改正项对时间求导,给出了计算北斗GEO卫星瞬时速度的直接求速算法并进行了验证测试。接着,探讨了微分步长对北斗GEO卫星位置微分求速精度的影响。进一步,采集真实北斗GEO卫星广播星历数据,选择代表卫星,采取不同微分步长,分别按照位置微分与直接求速计算卫星速度,并进行了比对测试。验证测试结果证明了北斗GEO卫星直接求速算法的理论正确性,该算法可作为计算GEO卫星瞬时速度的标准算法;比对测试结果则表明适用于北斗GEO卫星位置微分求速的合理微分步长应为1s。     

大偏心率火星探测器轨道外推的近似方法

对于快自转天体,如地球和火星,由于其非球形引力位中田谐项的振动频率与卫星运动频率相当,会给构造大偏心率环绕型探测器轨道分析解带来无法克服的困难。本文以轨道偏心率e≈0.9的＂萤火虫一号＂（YH-1）火星环绕探测器为例,探讨在精度要求不高的情况下,给出近似分析解和计算效率与其相当的数值解计算方案,以满足在计算资源有限的制约下快速轨道外推的要求。     

利用同步测距测速资料的卫星初轨计算新方法

本文给出了利用同步测距测速资料计算卫星初轨的一种新方法,其创新之处在于: (1)改进了初值选择方式。不仅适用于任意偏心率轨道;而且精度较高,从而可显著地改善迭代收敛过程。 (2)变换了条件方程的形式,使测距数据p和测速数据p分别独立地出现在方程右端,从而可方便地消减系统误差,并可进行加权处理,有效地提高初轨计算的定轨精度。实际计算表明,该方法是可靠的。     

基于改进粒子群算法辨识Volterra级数的目标机动轨迹预测

目标机动轨迹预测是空战态势感知和目标威胁评估的重要前提。针对传统目标机动轨迹预测模型复杂度大、预测精度低等问题，结合目标机动轨迹时间序列的混沌特性，提出一种基于相空间重构理论和Volterra泛函级数的目标机动轨迹预测模型。该模型首先采用0-1检测法验证了目标机动轨迹时间序列具有混沌特性；其次，利用C-C法确定嵌入维数和时间延迟，对目标机动轨迹时间序列进行了相空间重构；然后，引入Volterra泛函级数预测模型，为了克服高阶Volterra核函数求解复杂的难题，提出一种混沌变异自适应粒子群算法，构建一种基于改进粒子群算法辨识的Volterra级数预测模型，并将其应用于目标机动轨迹预测；最后，将所提算法与卡尔曼滤波算法以及机器学习算法进行单步和多步预测对比，同时将改进粒子群算法与其他智能算法进行性能比较。仿真结果表明：所提预测模型具有良好的单步和多步预测性能，改进的粒子群算法具有参数辨识精度高、收敛速度快的优点。     

Autonomous navigation to quasi-periodic orbits near translunar Ubration points

Libration-point missions have been very useful and successful. Due to the unstable natures of most of these orbits, the long-time stationkeeping demands frequent maneuvers and precise orbit determinations. Earth-based tracking will have to undertake much more responsibilities with the increasing number of libration missions. An autonomous navigation system could offer a better way to decrease the need for Earth-based tracking. Nevertheless, when an autonomous navigation system is applied, there are three important factors affecting autonomous navigation accuracy, i.e., the accuracy of initial conditions, the accuracy of measurements, and the accuracy of onboard dynamics for propagation. This paper focuses on analyzing the influence from the third factor and finding an appropriate navigation dynamics, which can satisfy the requirement of estimation accuracy but not cause too much burden for onboard computation. When considering the restricted three-body model and the bicircular restricted four-body model as navigation dynamics, the astrin- gency is not shown during the simulations. Meanwhile, when considering the influences of the Sun＇s direct and indirect perturbations and the eccentricity of the Moon＇s orbit, a new navigation dynamic model with the standard ephemerides is proposed. The simulation shows the feasibility of the proposed model.     

基于QSS的自适应多步校正算法及其在柔性航天器动力学中的应用

量化状态系统（QSS）算法是一种基于状态变量离散化的数值积分方法,该方法与基于时间离散的传统方法显著不同,QSS通过计算状态变量每次跃迁所需的时间来推进下一步的积分。在求解非刚性常微分方程时,QSS算法比传统算法更具优势,但它不适合求解刚性问题。为此提出一种基于QSS的自适应多步校正算法（AMCQSS）,该算法以QSS为基础、结合隐式多步法思想,在计算过程中可以自适应选择二步法或三步法,以有效提高求解刚性问题的精度和效率。通过对柔性航天器动力学的仿真求解,验证了算法的可行性。将该算法与ODE23tb、ODE15s、ODE45以及QSS等算法进行对比,结果表明AMCQSS算法既能保证求解的效率及精度,又具有较好的收敛性和稳定性。     

非均匀网格湍流大涡模拟高精度有限体积解法

为准确预测不可压复杂湍流,提出了一种可用于大涡模拟均匀或非均匀网格上的高精度有限体积法。该方法空间离散采用有限体四阶紧致格式,时间推进采用四阶Runge-Kutta法,压力-速度耦合应用四阶紧致格式的动量插值。通过直接求解顶盖驱动方腔流动证实了该方法具有近四阶的空间精度;并在此基础上,采用动态Smagor-insky亚格子应力模式,成功地实现了充分发展槽道湍流和后台阶湍流流动的大涡模拟计算,所得结果与直接数值模拟结果吻合良好,且采用非均匀网格可在比均匀网格数少的离散系统上得到同样满意的结果。结果表明,该方法是实现高精度湍流大涡数值模拟的一个有效途径。     

MIMO随机振动试验控制的逆多步预测模型法

多输入多输出（MIMO）随机振动试验在产品可靠性验证中扮演重要角色。逆系统方法作为一种时域方法，直接由参考信号（即期望得到的响应信号）生成驱动信号，但容易产生不稳定的结果。本文以系统的有限差分模型为基础，推导了逆多步预测模型，并建立了该逆模型与矩阵幂次算法结合的MIMO随机振动试验方法。通过悬臂梁模型进行双输入双输出仿真试验，验证了逆多步预测模型进行MIMO随机振动试验驱动生成的可行性。在三轴振动台上进一步验证该方法进行MIMO随机振动试验控制的能力。     

靶场T型架光电经纬仪精度检测方法的研究与应用

根据T型架光电经纬仪结构特点，设计了偏心方法，利用动态精度靶标对T型架光电经纬仪进行脱靶量、静态、动态及跟踪随机均方差的检测，通过某型号光电经纬仪的实际检测验证了该偏心方法的正确性。     

Dynamical behaviors and relative trajectories of the spacecraft formation flying

To describe the relative motion of spacecraft formation flying, this paper presents a method based on relative orbital elements, which is suitable to elliptical orbit with arbitrary eccentricity. The long time formation flying conditions are theoretically derived taken into account the relationship between relative motion and relative orbital elements. These conditions include that both the orbital periods of all participating spacecrafts should be the same and other relative orbital elements should be small enough. The expected relative distance of the spacecrafts would determine the magnitudes of such relative orbital elements. Theoretical analysis and numerical simulation results show that the spacecrafts with sufficient small relative orbital elements can keep long time formation flying without any active control when the orbital perturbations are not considered. The results also show that Hill's equation is only suitable for describing a short time formation flying.     

导航星座自主导航技术--自主星历更新

自主星历更新是导航星座自主导航的关键技术之一,包括卫星轨道精密确定和轨道短时预报两个方面的内容。本文在系统地论述导航星座自主导航的信息处理流程的基础上,重点提出星座卫星轨道精密确定的自适应Kalman滤波算法。系统仿真结果表明：通过星间双向测量数据的滤波处理,能够获得高精度的卫星星历和较小的用户测距误差（URE）,满足高精度导航定位需求,初步证明导航星座自主导航信息处理流程及其星历更新算法的合理性和可行性。     

不规则地形条件下双向DMFT电波传播特性算法研究

 在传统中心差分离散混合傅里叶变换(DMFT)抛物方程模型基础上,针对传统算法计算动态阻抗边界条件时的数值振荡及忽略后向传播与散射造成计算误差增大的问题,提出基于前-后向混合差分方法的递归双向DMFT(TW-DMFT)模型,以一阶前-后向差分方程代替传统的中心差分方程拟合阻抗边界条件,以递归双向传播模型改进传统单向预测模型.仿真分析了双向DMFT模型计算不同边界及媒介条件下单刃峰和实际地形的电波传播分布特性,与传统算法及几何绕射理论(GTD)进行比较,新算法提高稳定性及精度的同时弥补了传统算法的不足,更具有普适性.