求解含多个右端向量非对称线性方程组的自适应块拟最小残量IOM(q)算法

许多实际应用问题需要求解含多个右端向量的大型非对称线性方程组 ,通常是把原来方程组分成单独几个含一个右端向量的方程组 ,再用某种迭代法分别单个求解 ,而更加经济有效的方法是应用能同时产生几个迭代向量的块迭代法来直接求解。本文在 IOM(q)算法的基础上 ,提出一种求解此类方程组的块拟最小残量 IOM(q)算法 ,讨论了如何收缩掉已收敛的部分方程组以及如何从产生的块 Krylov序列中删除线性相关或几乎线性相关向量的自适应技术。数值试验表明 ,此新的自适应块算法比块 GMRES算法及其他相关算法具有更好的收敛行为、更少的计算量和 CPU计算时间 ,是求解此类方程组的一种更加经济有效的算法。     

用Newton子迭代方法计算前飞旋翼粘性绕流

通过求解Navier-Stokes方程数值模拟了直升机旋翼前飞非定常流场，为了同时保证计算的时间精确性和计算效率，时间推进格式采用了双时间推进方法，在该方法中，子迭代过程由十分高效的LU-SSOR方法完成，且使子迭代过程成为Newton子迭代，空间上应用中心平均的有限体积法进行离散，为了模拟前飞桨时间的相对运动，网格布局采用了运动嵌套网格方法，应用本文方法对-悬停流场进行了数值计算，计算结果与实验吻合较好，尽管缺乏实验数据的验证，对-有升力前飞状态的数值模拟结果是可信的。     

基于Krylov子空间法的柔性航天器降阶研究

带大型柔性附件的复杂航天器,其弹性变形通常用有限单元法描述.由于有限元节点坐标数目庞大将会给动力学方程求解和控制策略的实现带来巨大负担,因此对柔性太阳电池阵的动力学模型降阶进行研究.利用单向递推组集推导带太阳电池阵的航天器的刚柔耦合动力学方程,分别采用模态截断法、模态价值分析法和Krylov子空间法对太阳电池阵进行模型降阶,比较不同降阶模型和全有限元模型的动力学仿真结果.仿真结果表明,无论是否考虑大范围运动的影响,采用Krylov方法只需要较低的自由度就可以得到和采用有限元方法完全一致的结果.说明Krylov方法能够有效地降低航天器柔性附件的自由度,提高动力学仿真的效率,便于驱动控制的实现.     

求解大型非对称线性方程组的拟最小残量IOM（q）算法

不完全正交化算法（IOM（q））由于存储量和计算量小，常用来求解大非对称线性方程组。而此方法收敛过程常出现不规划振荡现象，从而影响了收敛速度。本文将拟残量最小的化性质加到IMO（q）算法中，提出拟最小残量不完全正交化算法（QMRIOM（q），这样收敛曲线光滑无振荡，从而大大加快其收敛速度，而且保留其存储量和计算量小的性质。     

基于PSF相关的星点定位算法

为提高非高斯分布星点定位的精度,提出了一种新的PSFC (point spread function correlation)星点定位算法,具有很好的抗噪性能和精度水平,易于工程实现.该算法利用互相关中的定义,通过确定系统所测定的PSF与星像灰度值之间的最大相似度来定位星点位置,PSF的测定是在作互相关计算之前全视场范围内可按照整像素和亚像素两种方式完成的.仿真实验结果表明：在带有噪声的星像图条件下,本文PSF相关法的定位精度最高可达到0.034像素,比质心法的精度提高1个数量级,比曲面拟合法的精度提高约20%.     

Singular formalism and admissible control of spacecraft with rotating flexible solar array

This paper is concerned with the attitude control of a three-axis-stabilized spacecraft which consists of a central rigid body and a flexible sun-tracking solar array driven by a solar array drive assembly. Based on the linearization of the dynamics of the spacecraft and the modal identi- ties about the flexible and rigid coupling matrices, the spacecraft attitude dynamics is reduced to a formally singular system with periodically varying parameters, which is quite different from a space- craft with fixed appendages. In the framework of the singular control theory, the regularity and impulse-freeness of the singular system is analyzed and then admissible attitude controllers are designed by Lyapunov＇s method. To improve the robustness against system uncertainties, an H∞ optimal control is designed by optimizing the H∞ norm of the system transfer function matrix. Comparative numerical experiments are performed to verify the theoretical results.     

遥测伪周期时间序列子序列异常检测算法

针对现有异常检测算法用于伪周期时间序列异常序列检测时易造成误差累积,导致序列周期与特征值上显著差异的不足,文章以卫星遥测伪周期时序数据为对象,综合两种常规分段方法的优势,提出了最大周期窗宽内基于极值的模式子序列分段算法。在此基础上,给出了一种基于均序列动态生成模型的子序列异常检测方法（AnomalySubsequenceDetectionmethodbasedonOptimizedSequenceModel,ASD_OSM）,并采用2次四分位距准则（DoubleQuantilerangescriterion,2Q准则）设置距离检测门限阈值,将超出阈值的序列判定为异常序列。某航天器传感器遥测子序列异常检测试验结果表明,提出的检测方法能够有效减少漏判,提高卫星遥测伪周期数据异常序列检测的准确性。     

基于H∞最优化与层次结构动态逆的非线性飞行控制器设计

针对扰动环境下的导弹飞行控制问题提出一种基于H∞最优化与层次结构动态逆(HSDI)的鲁棒非线性控制方法.设计层次结构动态逆控制器时将层次结构优化为3层,并按照时标分离的思想调整了各层次的状态变量,使得控制器更为简洁.在此基础上,将层次结构动态逆控制器与导弹6自由度(DOF)非线性模型联合作为H..控制的广义被控对象,通...     

求解大规模矩阵问题的Krylov子空间方法

求解大规模矩阵问题包括线性方程组和特征值问题等是计算数学和科学工程计算中的重大课题，最近几年，其研究工作取得了许多重大进展。文中给出大型线性方程组和特征值问题Krylov子空间方法若干进展的一个概述，其中包括作者对这些问题的研究成果。涉及的专题包括求解大型线性方程组的共轭梯度法、SYMMLQ算法、MINRES算法、GMRES算法、Lanczos双正交化算法、QMR算法以及这些算法的块格式；求解大对称特征值问题的Lanczos算法和块Lqnczos算法；求解大型非对称特征问题的Lanczos算法、Arnodi算法以及这些算法的推广。讨论求解大规模矩阵问题的加速技术和预处理技术。了一些有待进一步研究的问题。     

非对称线性方程组的广义拟最小向后误差算法

正交投影方法已经广泛应用于求解线性方程组.人们很少注意到斜投影方法,事实上斜投影方法更适合于解大型非对称线性方程组.本文的目的是考虑一标准来判断是否一个给定的近似值是合适的,并给出一个算法来计算线性方程组Ax=b的解使得向后误差算法满足某个优化条件.