KdV方程的数值解及数值模拟

文章在已构造的一种分片光滑的Lagrange型插值多项式空间中利用Galerkin法对非线性KdV方程进行求解。分析了数值格式的线性稳定性并进行了数值模拟,结果表明格式精确有效。     

基于精密时频预测的定位域和测速域RAIM可用性增强算法

能提供精确定位解的观测值集合提供的测速解并不一定是可信的,首先将传统的定位域RAIM算法应用到测速域的完好性监测中,提出将频差预测值用于辅助增强测速域RAIM可用性,使RAIM技术能更好的满足各类用户对导航服务完好性监测的需求.通过对多种预测模型的比较,提出将三阶多项式模型和递推阻尼最小二乘法应用到接收机钟差和频差的在线精密预测算法中,使钟差预测精度提高6倍.克服了传统的钟差辅助增强RAIM可用性算法中对晶振漂移率要求较高的限制,使钟差和频差辅助增强RAIM可用性的方法能更广泛地应用到各类接收机中.     

椭圆轨道航天器相对导航的多项式插值滤波算法

针对运行于椭圆轨道的航天器的相对自主导航问题,建立了惯性坐标系下的近距离相对运动方程,方程中没有圆轨道的假设,相对Clohessy-Wiltshire(C-W)方程具有更广的应用范围。将二阶多项式插值(Divided Difference 2,DD2)滤波应用于相对导航算法中,和扩展卡尔曼滤波(Extended Kalman Filter,EKF)相比,前者不需要对状态方程组进行微分处理,且对非线性函数统计特性近似到二阶程度。数值仿真比较了DD2滤波和EKF两种算法的性能,验证了利用惯性坐标系下的近距离相对运动方程实现相对导航的可行性,并分析了主动航天器的轨道和姿态确定误差对相对导航精度的影响。     

空间太阳望远镜上摆镜装卡变形的研究

通过控制一个光学摆镜的运动来补偿空间太阳望远镜抖动,摆镜在装卡力作用下变形的求解至关重要。分别运用解析法和接触非线性有限元法进行计算,解析法通过引入弹性地基理论得到由Zernike多项式表示的摆镜装卡变形微分方程,运用Moore\|Penrose广义逆矩阵法求解所有可能变形,再根据最小势能原理得到外力作用下的实际变形,同时直接得到光学设计者惯用的Seidel像差系数;有限元法通过接触非线性来模拟螺栓装卡的作用力。两种算法对同一算例分别进行计算,比较结果可知两者仅相差10.7%,表明解析法正确有效,有限元法模型正确。     

基于直接配点法的滑翔轨迹快速优化设计

介绍了基于五次Gauss-Lobatto多项式的直接配点法在再入飞行器三维轨迹最优化问题中的应用。首先给出了再入飞行器轨迹优化问题模型,其中运动方程为三自由度模型,性能指标选为到达指定地点飞行时间最短,控制变量则为无量纲升力系数和倾侧角。再入飞行过程中受到加热率、过载和动压约束,终端状态受到目标位置约束。然后,应用直接配点法将最优控制问题离散化为非线性规划问题,将动态优化问题转化为静态参数最优化问题。选取各节点和配点上的状态量和控制量作为优化参数。最后应用SNOPT软件包对参数最优化问题进行求解。仿真结果表明直接配点法对于再入飞行器轨迹初始参数取值不敏感,且求解过程具有一定的实时性。     

IIR数字滤波器的序列约束Chebyshev设计方法

 进行无限冲击响应(IIR)数字滤波器的无约束Chebyshev设计时,得到的滤波器幅值响应通常在过渡带存在明显过冲现象,且通带边缘的群延迟误差较大。本文提出一种序列约束Chebyshev(SCC)方法,在通带和阻带最大频率响应(FR)误差小于某给定值ρ的条件下,使过渡带频率响应误差(相对某个具有单调幅值特性的过渡响应函数)的Chebyshev模最小,并通过收缩和对分找出ρ的最小值。该方法不但可以使通带和阻带最大频率响应误差尽可能小,还能减小过渡带幅值响应的过冲以及带边群延迟误差。仿真实例表明了该方法的有效性。     

低速流场全方位校准五孔探针测量数据处理方法

论述了低速不可压流场全方位校准五孔探针测量数据处理的计算机程序,并对计算结果进行了对比分析。研究表明,采用速度迭代法构筑二阶拟合多项式对校准和测量数据进行处理,角度偏差不大于1°,速度偏差不大于1％,完全可以满足工程上定量测量三维流场的要求,具有较高的工程实用价值。     

计及随机噪声的频域多输入多输出模态参数识别

提出了一种计及随机噪声的多输入多输出频域模态参数识别方法。该方法基于连续时间域的频率响应函数有理分式模型,利用Forsythe正交多项式改善求解性态,属频域的极大似然类估计方法。在最小二乘估计的基础上,利用"先验"的随机噪声信息进行优化迭代的极大似然估计,得到更优的模态参数识别结果。根据Cramer Rao下界不等式,在几乎不增加计算量的情况下获得所识别结果的不确定性信息,提高了模态参数识别的可靠性。采用GARTEUR飞机模型作为仿真算例,对本算法进行了仿真验证。     

随机响应面法在结构随机响应计算中的应用

传统的响应面方法以一般多项式逼近结构的随机响应,但这种方式并不能保证收敛性.以随机多项式为基础的随机响应面方法,可以弥补这一不足.两个数值例子和一个发动机轮盘实例来验证这一方法的有效性.结果表明:随着随机多项式次数的增加,随机响应面所得到的响应概率密度曲线愈加趋近于由蒙特卡罗方法所得到的概率密度曲线.      

基于径向基函数改进的Delaunay图映射动网格方法

为了提高大变形下动网格生成的效率和质量,基于径向基函数插值技术对原始Delaunay图映射动网格方法进行了改进。首先通过带多项式的径向基函数插值方法计算出背景网格远场边界点的位移;然后将背景网格位移插值到计算网格;最后利用衰减函数将计算网格远场位移衰减为零。基于矩形旋转网格变形实例比较了改进方法与原始Delaunay图映射动网格方法之间的差别,并研究了径向基基点数目和衰减函数类型对变形网格质量的影响。矩形旋转网格变形实例说明该方法能够精确恢复出Delaunay背景网格的旋转特性。NACA 0012翼型、NLR 7301两段翼和M6机翼网格变形算例进一步证明,通过添加控制点,该方法能够不重构背景网格实现大变形下高质量动网格的生成。