结构输出响应概率密度估计中分数矩求解方法

鉴于概率不确定性背景下基于分数矩极大熵准则的结构可靠性分析方法具有较大的效率与精度优势,综合研究并给出了可以用于极大熵准则中约束条件输出响应分数矩求解的3种分数矩求解方法,包括降维积分(DRI)方法、稀疏网格积分(SGI)方法和无迹变换(UT)方法。阐述了分数矩求解原理及过程,给出了方法的计算效率,并分析了方法的适用性。3种分数矩求解方法在确保计算精度的同时可以很大程度减少结构输入-输出模型的调用次数,大幅提高统计分析效率。通过与Monte Carlo仿真分析法对比,验证了3种分数矩求解方法的正确性与高效性。      

基于MapReduce的CME参数识别模型并行计算技术

日冕物质抛射（Coronal Mass Ejection,CME）参数识别模型是太阳风预报过程的重要组成部分.在空间环境预报业务中,为提高太阳风预报的准确率,需要提高CME参数识别的精度.模型以计算任务串行的方式运行,运算效率低导致模型运算时间长,不能满足这种需求.CME参数识别模型的物理运算过程相互不独立,其在单节点上的运行方式不能满足并行化要求.基于MapReduce的并行计算框架,改进了CME参数识别模型的计算流程,提出CDMR（CME detection under MapReduce）方法,实现了CME参数识别模型的并行计算,并对比分析CME参数识别模型在串行计算和MapReduce并行计算下的运行时间,提高了模型的识别精度和计算效率.      

基于光电编码器的电机变周期速度观测器设计

针对光电编码器在伺服系统应用中速度估计易受系统振动干扰和存在量化误差的问题,设计了一种基于变采样周期观测器的速度估计算法.通过观测器的积分运算抑制振动干扰;通过采样周期和光电编码器脉冲的同步,使量化误差最小;为了减小观测器离散域极点变化对算法性能的影响,保证变周期观测器在较大速度范围保持较高精度,给出了一种观测器增益系数的计算方法.实验表明该算法可以在保证动态性能的同时,有效抑制振动和量化误差引起的估计误差,提高速度估计精度.      

基于CUDA的超声二维声场EFIT仿真

随着图形处理器（GPU）的快速发展,基于计算设备统一构架（CUDA）可以方便地将并行计算技术应用于超声声场数值仿真计算,极大地提升计算效率。阐述了弹性动力学有限积分算法（EFIT）的原理,在采用CPU实现带吸收边界的钢材料二维点源激励声场仿真的基础上,基于GPU实现了仿真模型的并行计算,介绍了GPU程序的设计流程和参数优化方法,包括纹理内存使用、吸收边界优化和数据传输优化。对比了相同条件下CPU和GPU仿真计算的耗时和平均计算效率,定量分析了GPU对于EFIT模型效率的提升。比对结果表明,EFIT具有良好的并行计算条件,采用并行计算方法能够有效提升模型计算速度,对于复杂声场仿真应用具有广阔的应用前景。      

不确定结构区间特征值上下界的并行解法

当工程结构参数包含不确定因素时,结构的固有频率也将是不确定的.这就需要讨论不确定性振动问题中广义区间特征值的求解方法.在Deif标准区间特征值求解定理的基础上,通过区间分析,将特征值的上下界分解成2个广义特征值问题进行求解.基于此求解方法的并行性分析,给出并行求解算法,克服了求解区间问题计算量大的缺点,使传统串行机或者串行算法难以解决的区间特征值问题得以较好的解决.      

一种基于DCT的SAR图像自聚焦算法

研究了一种基于离散余弦变换(DCT, Discrete Cosine Transform)的SAR(Synthetic Aperture Radar)图像自聚焦算法.该方法从复图像域出发,通过在距离压缩相位历史域引入相位误差模型,改变图像的聚焦程度直至一维像的DCT序列在高序数区域的权值达到最大,从而完成误差校正.同相位梯度自聚焦算法相比,该方法无需在图像域分离出强点目标,因此特别适用于无任何明显特征的图像.由于DCT存在快速算法,使得该自聚焦算法计算量较少,更易实现.实测数据及仿真数据的成像结果证明了此方法的可行性.      

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

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

基于显式制导的月地返回轨道中途修正研究

月地返回轨道存在各种摄动误差,终端约束复杂,有必要对其进行中途修正研究。显式制导法通过二体轨道与精确轨道之间的差别进行多次迭代求解给定时刻所需的修正速度。文章利用显式制导法,采用月球段及地球段分段进行中途修正的策略,给出了基于分段落点预报显式制导的月地返回轨道中途修正方案。该方案无需计算雅克比矩阵,算法简单、计算快速、实用性强,能满足再入点参数要求。算例仿真与蒙特卡洛仿真验证了该方案的适用性。     

高可靠性航空电子设备热分析中的有限体积法

阐述了热对系统可靠性的影响,对比了数值法求解温度场分布的优缺点.针对航空电子设备热分析中同时要求满足高的分析精度和快的计算速度的问题,提出用有限体积法求解温度场分布的计算方法.根据典型航空电子设备的结构特点,建立了温度场的数学模型,给出了边界条件.在此基础上,用有限体积法对待求解方程离散化,给出了离散化过程中网格划分及边界的处理方法.结合实例进行了温度场求解,计算结果与热测量结果进行了对比,得到了较高的分析精度,验证了方法的有效性.最后对计算误差进行了分析,并提出了进一步热设计的合理化建议.      

基于时域有限差分的微放电仿真算法

为克服传统微放电阈值预测方法建模粗糙、精度低的缺点,提高阀值预测精度和效率,研究了基于时域有限差分的精确微放电阀值预测算法。基于时域有限差分算法和粒子追踪算法,通过时空网格自洽互耦实现复杂结构微波器件微放电阈值准确计算。其中,共形和并行算法是提高微放电仿真精度和效率的关键。文章基于空间蛙跳策略实现了基于时域有限差分的微放电仿真算法,利用算法分析了典型微波器件的微放电阈值,仿真与实验结果吻合良好,误差小于0.3dB;同时,并行效率最高可达83%,验证了算法的准确性和高效性。     

按区域惩罚划分的并行多目标遗传算法

解决多学科设计优化问题的多目标遗传算法通常面临着大计算量的挑战,提出了一种新型的并行化算法来提高其效率.全局个体均匀的分布在各个进程,首先从所有的进程中获取全局范围的Pareto最优解极值,并发送给每个进程,再由这些极值来构造各个进程自己的惩罚函数.通过惩罚函数给个体添加约束来划分各个进程的收敛区域,同时采取优化措施保证每个进程加速收敛并且收敛区域没有重叠和遗漏,这样每个进程只需收敛到特定的一段Pareto最优解,降低了计算量;同时由于进程间交换的数据量小,保证了效率的提高.通过与串行算法(NSGA2)和其他的并行化算法比较,显示了该算法的有效性和先进性.     

多尺度子域基函数矩量法解的稳定性与迭代法

通过对几个具体实例的计算，证实了在电磁场边值问题的矩量法求解过程中，选用不同的基函数其解的数值稳定性会有很大的差别。从而有可能选择一种数值稳定性比较好的基函数来进行矩量法求解。另一方面，提出了选用多尺度子域基函数，并通过多次迭代来提高矩量法解的精度的算法。实例计算表明，该算法不仅可以有效地降低系数矩阵的条件数，而且还降低了对计算机内存容量的要求，是一种有效的、稳定的算法     

2-SGCMGs与磁力矩器的对地姿态混合控制方法

针对故障后仅剩两单框架控制力矩陀螺(SGCMG)可工作的对地定向三轴稳定卫星姿态控制问题进行了研究,提出了2-SGCMGs系统与磁力矩组合的混合控制策略及方法,以克服两SGCMG欠驱动控制的鲁棒性问题.首先,给出2-SGCMGs零动量方式的标称框架角构型选择计算过程.然后,结合标称框架角构型,构造了一种不同于沿传统体轴...     

单目SLAM直线匹配增强平面发现方法

针对微小型机器人及无人机系统日益迫切的轻量化视觉导航需求,提出了一种多维几何特征单目视觉三维环境建模方法。单一点特征单目SLAM制图方法地图描述效率相对较低,噪声容忍性能需要进一步提高。将线和面特征引入单目SLAM的三维地图构建过程,提高系统三维空间建模的搜索速度和稳定性。利用快速直线搜索算法,并基于二维直线匹配生成三维空间直线。现有基于三维空间特征点生成最小采样集的J-Linkage算法需要的倾向向量维数较高,完成单目SLAM常见场景三维平面聚类所需的计算量大。通过点线特征结合以及直线增强的J-Linkage算法可以提高特征平面聚类速度和稳定性,减少系统三维空间表达的冗余信息。     

基于非线性观测器的控制力矩陀螺操纵律设计

为改善操纵性能,将单框架控制力矩陀螺(SGCMG)操纵律设计问题转化为非线性系统的状态观测问题,并基于状态观测理论,导出了一种新型的SGCMG操纵律.通过适当的参数选择,SGCMG操纵律可使操纵误差在理论上渐近收敛至零;不用计算Jacobi阵的伪逆,而代之以Jacobi阵的转置,从而避免了由Jacobi阵求伪逆带来的一系列问题.同时,该操纵律形式简单,计算量小,易于实现.对应用在航天器上的某SGCMG系统的仿真结果表明,上述操纵律是可行的.      

模拟三维重力波传播过程的并行数值模式

从三维大气运动的基本控制方程出发, 基于MPI消息传递接口和区域分解的思想, 建立了模拟三维可压缩大气中重力波传播过程的并行数值模式. 在对垂直方向进行区域分解的基础上, 针对跳点网格的特点, 以垂直速度主格点线为子区域的下边界, 状态变量主格点线为子区域的上边界进行区域分解. 利用MPI消息传递接口来传递计算各子区域中变量时所要用到的相邻子区域上相应变量, 从而顺利地完成对整个计算域的并行计算. 根据线性重力波理论, 通过模拟小振幅重力波的传播过程对所建模式进行了验证. 结果表明, 模式可以很好地模拟小振幅重力波在三维空间中的传播过程, 模拟的重力波振幅随着传播高度的增加以指数形式增长, 对能量传播路径以及有效扰动位能和扰动动能的模拟结果均与理论预测的结果吻合很好. 另外, 随着进程个数的增加, 完成相同的计算量所用的计算时间也显著减少. 这表明本文建立的并行数值模式不但能够很好的模拟重力波的传播过程, 而且能够有效地节约计算时间.