具有约束项的自校准卡尔曼实时滤波定轨方法

针对在传统卫星定轨方法中存在的轨道根数解算复杂 ,实时定轨精度低 ,以及系统误差校不准等问题 ,提出并详细阐述了具有约束项的自校准卡尔曼实时滤波定轨方法。该方法复杂度低、实时性强、计算量小 ,仿真结果表明该方法同时具有收敛速度快、定轨精度高的优点 ,具有很好的应用前景     

刚性方程的二阶Gear方法的注记

构造了一个求解刚性方程的二阶差分格式,它是隐式的,并采用简单迭代求解,对于稳定系统,迭代收敛性条件与步长无关,而对不稳定系统,迭代收敛性条件对步长却有所限制,利用这一限制性条件去克服解刚性方程的二阶Gear方法时出现的“危险性”。     

高超声速飞行器滑模控制参数整定方法设计

本文针对高超声速飞行器滑模控制人工试凑的参数整定方法较为繁琐、效率低的问题,改进了连续动作学习自动机（Continuous action reinforcement learning automata, CARLA）算法并将其应用于滑模控制参数整定问题中,改进后的算法通过对回报函数的设计有效克服了常规CARLA算法收敛速度慢、易受干扰、求解效率低的问题。该算法引入控制性能指标评价函数,在迭代中学习阶跃响应的经验数据,实现了控制参数自整定。仿真表明,对于阶跃响应问题,本文提出的算法能够在100次迭代中整定出一组高品质的控制参数,完成对给定指令的快速准确跟踪,与遗传算法、模拟退火算法相比,在求解速度上具有显著优势。由于该方法不依赖于模型,除了滑模控制器参数整定外对其他控制方法的控制参数整定问题也有一定适用性,具有推广应用价值。     

基于固定模糊度PPP的低轨卫星实时精密定轨

传统动力学定轨法受制于动力学模型精度，传统几何定轨法精度受限，只能达到亚米级，而基于精密单点定位（PPP）模式的几何定轨法一般采用浮点解，定轨精度及可靠性较基于双差模式的相对定位较差。为提高PPP模式低轨定轨的定位性能，利用中国区域内外的IGS测站计算出当前所有卫星的宽巷和窄巷相位小数偏差产品，对经过中国大陆区域上空的国产低轨卫星海洋二号（HY-2）和资源三号 (ZY-3) 卫星进行固定模糊度PPP的定轨解算，与事后精密轨道结果进行比较，分析其外符合精度。结果表明：仅利用约10min弧段的HY-2和ZY-3卫星数据，切向与径向的定轨精度可达2cm左右，法向为5cm左右，较浮点解定轨精度大幅提升。基于固定模糊度PPP的定轨方法能够满足厘米级的实时精密定轨。     

基于多模型最优融合的双星定位系统一体化精密定轨方法

针对卫星动力学模型复杂且不准,考虑到卫星定轨中待估参数在时间和空间上的相关性,提出了一种基于多模型最优融合的双星定位系统一体化参数建模的近地卫星精密定轨新方法.利用节点自由分布B样条描述卫星运动,实现了对卫星粗略动力学模型的抑制作用;同时结合双星观测模型,使该方法转化为关于求解卫星轨道样条表示参数和定轨系统误差的多模型融合的非线性优化问题;通过引入模型结构确定最优融合权值的选取准则,在最小二乘准则下,采用非线性最优化方法搜寻样条的最优节点分布,得到了待估参数的最优估计,完成了近地卫星的精密定轨.理论分析和仿真计算表明,该方法确实有效,不仅提高了卫星的定轨精度,而且使状态估计的结构更加稳定.      

“北斗”导航卫星定轨残差特征提取与分布检验

针对"北斗"导航卫星定轨残差中非建模系统误差问题,提出基于总体平均经验模式分解(EEMD)与Hilbert谱分析相结合的定轨残差特征提取方法,并对定轨残差正态分布特性进行卡方检验。首先,分析总体经验模式分解原理,提出滤波辅助的改进EEMD方法与Hilbert谱特征提取结合,建立定轨残差特征提取模型,并阐述卡方检验正态分布原理;然后对北斗导航卫星伪距与载波相位数据定轨残差特征进行分析;最后将提出的方法应用于多测站、多北斗GEO卫星的定轨残差分析中。结果表明,EEMD与Hilbert谱方法有效提取出了定轨残差中1天的轨道周期项,且EEMD处理后定轨残差的卡方统计量为5.5,其值小于卡方检验临界,可视为正态分布。该方法可为北斗定轨残差中非建模系统误差分离、提高定轨内符合精度提供有效技术手段。     

电推力器中复杂电场求解的超松弛并行算法

目前,电推力器中复杂电场的三维泊松方程求解仍存在速度偏慢、并行性欠佳等缺陷。通过分解三维求解区域,针对复杂电场提出一种超松弛并行迭代算法。将求解区域划分为多个子块,利用超松弛迭代格式构造出若干分组显式格式,分别给出不同迭代步数下的求解方程以进行并行计算。通过数值模拟,计算时间至少缩短一半以上,该P SOR方法较传统迭代格式在快速性方面取得较大进展,对电推进领域的数值仿真研究起到促进作用。     

径向速率估值的双屏电子篱笆数据关联方法

针对双屏电子篱笆系统在初轨确定前需要完成同一目标观测数据的数据关联问题,文章提出了一种计算形式简单且不依赖于循环迭代定轨计算径向速率的估值方法,并设计了相应的数据关联准则。仿真分析表明,这种数据关联准则与依赖定轨计算径向速率的关联方法相比,判别正确率相近,但计算时间不到后者的10%,算法的效率更高。     

基于目标特性的航天器太阳光压面积计算方法

太阳光压是影响深空探测航天器轨道确定与预报精度最主要的摄动力.针对实际任务需求，采用了一种基于目标特性的光压面积建模与计算方法，根据航天器形状、尺寸、表面材料以及材料光学特性等信息，实现了分析型光压模型的建立与求解，提高了计算效率和精度，可快速计算目标在光照方向上的光压面积、投影面积以及光压比例因子等参数.通过长方体光压面积理论值与仿真值的对比，验证了该方法的准确性和有效性.针对复杂结构探测器开展了光压面积计算，可为深空探测航天器精密定轨中的光压模型解算、定轨及预报提供参考.      

火焰筒头部流场研究之二——二元轴对称等温廻流区计算

本文将火焰筒头部的等温流场,用二元轴对称不可压定常流来处理并用数值方法解出初步结果。基本方法是用涡量流函数法,双方程粘度模型,采用交替方向隐含格式差分求解,在边界及迭代初值提出自己的看法,初步计算结果与[10]之实验定性一致。     

基于GPS的中、高轨道航天器定轨研究

利用GPS确定空间航天器轨道，尤其是GPS在中、高空间航天器定轨中的应用，针对中、高轨道中可见星少及GPS信号微弱等问题，提出了有效的微弱信号检测算法，提高了GPS接收机的接收能力，同时提出了利用伪加速度测量和轨道约束进行定轨解算的新方法，使得在HEO(高偏心轨道)中利用两颗可见星可以解算定位解，对GEO(静止地球轨道)只用一颗星就可求解，仿真结果表明该方法具有很好的效果。     

基于双星系统定轨的数学方法

通过对双星系统定轨原理的分析,结合卫星的动力学模型,采用联合多时刻观测系统的数据融合方法,得到了一种解算卫星轨道状态的解析方法,同时还给出了误差传递方式。该方法不需要标称轨道,计算方便。理论分析和仿真计算表明,方法确实可行,且定轨精度有所提高。     

基于Bayes统计模型的卫星融合定轨方法研究

根据目前天基导航系统现状,结合中国对低、中轨卫星精密定轨的要求,给出了天地基信息融合定轨的原理;结合卫星待估融合参数的先验信息,提出了基于Bayes统计模型的卫星精密定轨方法;在卫星观测的线性化融合模型中引入观测噪声,利用概率估计融合模型,根据Bayes理论进行卫星状态改进量的最大后验估计,并分析了Bayes估计方法的定轨精度;依据期望融合的待估改进量方差最小规则建立了相应的参数求解算法;最后以导航融合测控系统中测距和测速数据的融合定轨为例进行了仿真实验,表明该融合方法能够得到很好的定轨效果。     

基于机动力模型的GEO卫星恢复期间定轨

北斗卫星导航系统(BDS)中GEO卫星频繁的轨道机动对高精度、实时不间断的导 航服务需求提出了更高要求, 如何在短弧跟踪条件下提高GEO卫星轨道快速 恢复能力, 是提升导航系统服务精度的关键因素. 针对该问题, 本文提出了基 于机动力模型的动力学定轨方法, 尝试利用高精度的C波段转发式测距数据, 辅 以机动期间的遥测遥控信息建立机动力模型, 联合轨控前后的观测数据进行动 力学长弧定轨. 利用BDS中GEO卫星实测数据进行了定轨试验与分析, 结果表明, 恢复期间需要采用解算机动推力的定轨方法, 联合机动前、机动期间和机 动后4h数据定轨的轨道位置精度在20m量级, 径向精度优于2.5m. 该方 法克服了短弧跟踪条件下动力学法定轨和单点定位中的诸多问题, 提供了解决 GEO卫星机动后轨道快速恢复问题的技术方法.      

2016HO3探测任务星载光学观测量建模及定轨

以中国首次小行星探测任务为背景，根据星载相机获取的光学影像构造三种观测量，分别为小行星相对于航天器的高度角和方位角、赤经赤纬以及探测器与行星/小行星之间的夹角，分析了其在探测器定轨中的作用。仿真定轨结果表明，观测时长为100 h，探测器三轴位置误差小于50 km，满足工程上对巡航段的轨道精度要求，但x和y方向的位置和速度分量具有较强自相关性。此外还发现，使用单一观测数据类型比联合观测量的定轨精度低3～4个量级，第三类观测量相对于其他两类观测量在定轨精度方面具有显著作用，这表明在2016HO3探测中，利用太阳系大天体的位置信息有助于约束探测器轨道，提高探测器的定轨精度。     

马尔可夫更新方程的迭代算法研究

冷储备系统的动态结构变化使得系统可用度分析必须基于瞬态求解,针对频域两部件冷储备系统马尔可夫更新方程解算困难的问题,提出基于时域伏尔特拉方程的迭代算法,证明了该迭代方程的唯一性和收敛性,并给出迭代算法的误差分析.理论和实例分析表明,当冷储备系统部件的工作寿命和维修时间分布是有界函数时,利用迭代算法可以容易实现时域马尔可夫更新方程的求解.      

人工释能法求解线性代数方程组

本文建立求解线性代数方程组的人工释能法迭代格式,并给出该迭代格式收敛的充分必要条件和最佳释能时机,同时给出实现人工释能法的数值计算方法及其数值稳定性条件。     

火星探测器行星际转移段定轨精度评估

天问一号是中国第一次实现地火转移行星际飞行的探测器，在长达202天的行星际转移飞行期间，共经历了4次中途修正及1次深空机动控制，在2021年2月10日成功进行了近火制动，被火星捕获而进入环火轨道。本文对探测器行星际转移期间的动力学模型进行了分析，制定了转移飞行期间定轨积分中心转换原则：在探测器飞出地球影响球后，定轨积分中心需要由地心更换为日心；对不同版本行星星历表的使用进行了分析，确定了使用DE436行星历表进行计算对定轨影响最小。根据探测器行星际转移飞行的特点，制定了一种基于逐日迭代定轨策略的精度评估方法。基于实测数据分析，验证了该方法的有效性，火星探测器行星际转移期间定轨位置误差优于2 km，速度误差优于20 mm·s^–1 (1σ)。     

基于改进多目标遗传算法的舰尾紊流模拟方法

为提高舰尾紊流自由大气紊流分量仿真的可信度,提出了一种紊流数值模拟的新方法.首先,使用带有修正系数的Euler前向差分格式表示各个方向的紊流序列,同时结合智能算法的思想,把紊流相关性检验中的均方差误差和相关函数误差作为待优化目标函数,将修正系数的选择看成一个多目标优化问题,并采用改进的多目标遗传算法进行求解.最后,通过仿真算例验证了本文方法的正确性与合理性,计算结果表明该方法可以根据不同的采样步长灵活地生成所需紊流.尤其在小步长情况下,亦可得到很好符合理论值的紊流序列,可以满足虚拟飞行实验的要求.      

非定常N－S方程的线化及其差分求解

在微幅简谐振荡的假设下，对计算平面内非定常Ｎ－Ｓ方程进行了线比，并在ＭａｃＣｏｒｍａｒｋ于１９８０年提出的半隐格式的基础上，发展了一种适用于求解计算平面内线化Ｎ－Ｓ方程的差分格式．用该差分格式求解计算平面内的线化Ｎ－Ｓ方程时，可以避免求解五对角方程；且时间步长可以取得较大，计算效率高．