Legendre-Gauss拟谱法求解最优控制问题

 提出一种新的求解基于常微分方程(ODE)和微分代数方程(DAE)的最优控制问题的数值方法。本方法基于直接配置法,通过Legendre-Gauss拟谱法同时离散化状态变量和控制变量把最优控制问题转化为一个非线性规划问题。与传统的直接转换法相比,本方法具有精度高、计算量小、结构简单的特点,而且可以求解最优控制“多相”问题。数值结果表明,本方法是一种通用的精度较高的最优控制直接数值求解法,可用于求解ODE/DAE最优控制问题。     

求解退化线性代数方程组的小参数法

文章介绍了求解退化线性代数方程组的小参数法,推导出简便的求解公式,给出了小参数的选择方法,并且证明了方法的稳定性和高精度,最后给出了一些数值结果.     

求解非线性代数方程组的一种建议方法

飞行动力学研究中常遇到的求解非线性代数方程组的问题，现提出一种求解实函数非线性代数方程组的组合迭代方法，以多步交替使用的方式将梯度下降法和拟牛顿代法结合起来，综合两者的优点，组成一种对大范围内偏离精确解的任意初值均能收敛，且有一定收敛速度的迭代法，通过算例，对三种方法进行了对比和分析，计算结果证明，该方法是优越的。     

用泰勒级数求解非线性代数和微分方程组

提出一种求解非线性代数方程和非线性常微分方程的新方法。所论方程组被转换成增量形式，而解被表示成泰勒级数。在代数方程（多项式形式）的情况下，此算法归结为求解一系列线性方程组，而且在每一增量步中系数矩阵是不变的；在常微分方程组的情况下，此算法归结为一组递归算式，并不要求解方程组。此方法固有的自动走步功能，可保证得到收敛的解，从而大大节省计算时间。     

带约束多体系统Lyapunov指数的数值计算方法

应用第1类Lagrange方程建立的带约束多体系统动力学方程为非线性微分-代数方程组.利用增广法,将其转化成了常微分方程组,并表示成矩阵形式.根据方程的特点,给出了方程的Jacobi矩阵的具体表达式,提高了计算效率.在此基础上,给出了Lyapunov指数的数值计算方法,并通过对具体的树形和非树形多体系统进行Lyapunov指数数值计算,结合相图和Poincare映射对系统的动力学特性进行了分析,表明了该方法的可行性和有效性.      

喷管温度场计算的加权余数有限元方法

喷管温度场计算是固体火箭发动机喷管设计中遇到的一个重要问题。利用差分法计算温度场已有很多文献报导,本文提供一个加权余数有限元法计算喷管的温度场,其突出的优点是能与热结构计算程序联用,组成热分析和热结构的统一程序,因为加权余数有限元与变分原理有限元在节点网格的划分,有限元几何参数的计算,变带宽的一维储存和代数方程组的求解都有相同之处。为了说明方法的精度和可靠性,本文给出了有限元计算与解析解的比较,有限元计算与地面试验结果的比较,这些比较表明:本文提供的方法是可靠的,可供