二阶中立型方程周期解的存在性

利用Ｆｏｕｒｉｅｒ级数理论研究二阶常第数线性中立型方程的周期解，得到了在｜ｃ｜＝１的条件下，二阶中立型方程周期解存在的充要条件；在｜ｃ｜＝１条件下，方程周期解存在的充分条件，并给出了具有特殊时滞的方程周期解存在的充分条件。     

二阶多项式自治系统代数曲线解的判定法

为了以解析的方法和代数的方法寻求二阶多项式自治系统的代数曲线解，并为进一步研究非线性系统的可积性和得域上多项系统的定理论，证明了关于复域上二元多项式的一个整除定理，并以ＶａｎＤｅｒＰｏｌ程及一个二阶二次自治系统为例，说明这一整除定量可以作为二阶多项自治系统是否存在代数曲线解的主要判据。     

环面Fuchs方程解的性质及其可积性

研究环面Ｔ＾２上只有一个正则奇点的Ｆｕｃｈｓ方程，得到了参数ｒ＝６时，方程有一个椭圆函数解，其任何解皆为半纯函数，以及方程的单值群为可解群的结果。在此基础上，将Ｒｉｅｍａｎｎ球面上Ｆｕｃｈｓ方程的可积性概念推广到环面上，并得到一系列环面Ｆｕｃｈｓ方程都是可积的结果。     

一种基于微分代数的任意阶空间目标轨道传播方法及其分析

微分代数方法可以在不改变当前算法计算过程的基础上给出函数对自变量任意阶导数的精确值。本文给出了一种基于微分代数的任意阶空间目标轨道传播方法。本方法首先将初始微分代数代入轨道传播方程，然后用获得的高阶导数构造新的高阶微分代数。用新的高阶微分代数迭代前述过程可求解空间目标状态对时间的任意阶导数。最后，将任意阶导数代入泰勒展开公式求解空间目标轨道单步传播。本方法要求轨道传播方程采用的摄动力模型在轨道传播积分区间上是解析的。本文通过仿真分析验证了所提方法的有效性。     

铁木辛柯梁固有振动频率的边界元解法

根据傅里叶变换推导具有两个广义位移的铁木辛柯梁固有振动的基本解.利用加权残量方法,从控制微分方程出发建立边界积分方程,进而根据边界条件得到频率方程,采用代数特征值方法和影响系数方法求解频率,并分析了两种方法的特点.以杆为例证明了对于一维均匀结构,对不同的边界条件利用边界元方法(BEM,Boundary Element Method)都可以得到精确频率.将铁木辛柯梁的BEM结果与有限元结果和精确解进行了比较.     

六杆机构四位置运动生成的解域综合理论与方法

针对Stephenson-Ⅱ型六杆机构四位置运动生成问题,给出基于布氏曲线的广义布氏曲线概念及其方程,同时给出解域综合方法。该方法首先求解出满足要求的解曲线,而后利用雅可比矩阵行列式值符号及机构运动连续性条件剔除解域中存在运动缺陷的机构,最后根据实际需求附加约束条件形成可行解域图。给出的综合示例验证了解域综合方法的可行性和有效性。该方法简单且便于编程实现,同时对其他多杆机构的综合也有参考借鉴意义。     

基于C-W制导的瞄准点修正技术

在分析C-W冲量制导误差源的基础上,指出主要摄动为C-W方程的线性化误差,相对J2项摄动误差和发动机的冲量假设误差.在求解得到主要摄动作用下相对运动方程解析解的基础上,提出了瞄准点修正技术,在计算量增加不大的前提下,不仅提高了C-W制导精度,而且实际制导曲线趋近于最优交会轨迹,燃料消耗接近最优解.     

可渗秀各向同性多孔材料动静压轴承性能研究

利用Ｒｅｙｎｏｌｄｓ和Ｌａｐｌａｃｅ方程以及Ｄａｒｃｙ定律，导出了可渗透各向同性动静压窄轴承的计算方程。计算结果显示出这种轴承具有承载力高和姿态角小的特点，并指出设计参数之间存在优化组合。     

自动控制下,飞机运动过渡过程曲线简捷计算法

为了验证自动控制下,飞机过渡过程曲线是否符合调节质量的要求,通常需联合解算飞机及自动器的运动方程。纵向运动中飞机及自动器的特性方程通常为六阶方程。侧向运动中如果考虑航向及横向自动器方程及飞机侧向运动方程,飞机及自动器的特性方程将达九阶的高次方程。参考中曾提出应用“带后时间”的等效作用,以简化飞机的纵向及侧向运动的计算工作。但进行计算时,要求适当估计“滞后时间”的大小。本文提出应用近似等式及自动器方程中的时间常数?与反馈时间常数?o.c即可在降低特性方程方次的条件下算出飞机运动的过渡过程曲线。     

可渗透各向同性多孔材料动静压轴承性能研究

利用Ｒｅｙｎｏｌｄｓ和Ｌａｐｌａｃｅ方程以及Ｄａｒｃｙ定律，导出了可渗透各向同性材料动静压窄轴承的计算方程．计算结果显示出这种轴承具有承载力高和姿态角小的特点，并指出设计参数之间存在优化组合．     

广义正交多项式在时变系统跟踪问题中的应用

采用广义正交多项式的展开技术,利用其积分、乘积运算阵,对于时变线性系统最优控制中的跟踪问题,直接从最优控制的性能泛函入手,将最优控制问题转化为代数极值问题,从而避免了求解非线性Riccati方程,得到了一个较为简便的算法,并结合实例说明了方法的可行性.      

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

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

基于Riccati方程解的高增益未知输入观测器设计

针对具有未知输入和测量噪声的一类Lipschitz非线性系统,研究了状态估计、噪声估计及未知输入重构问题.通过将输出噪声看作扩展状态,把原系统转化为描述系统.针对描述系统,首先基于Riccati方程的解,提出了一种高增益观测器设计方法,实现对系统状态的估计和测量噪声的重构;之后,设计二阶高增益滑模观测器精确估计输出的微分,并利用状态和输出微分的估计,提出了一种未知输入的重构方法.在一Riccati方程有解的前提下,所提出的未知输入和测量噪声的重构,均适用于强时变信号.最后,对一个实际模型仿真,验证所提出方法的有效性.      

Optimal guidance based on receding horizon control for low-thrust transfer to libration point orbits

This paper addresses the design and computation of a guidance law for a transfer mission from an orbit near the Earth to a halo orbit around the libration point L₂ in the Sun–Earth system. The guidance law, which is designed based on receding horizon control and compensates for launch velocity errors that are introduced by inaccuracies of the launch vehicle, is solved using the generating function method. During the design of the closed-loop guidance law, the entire transfer mission, which is considered a nonlinear optimal control problem, is evaluated to obtain a nominal reference trajectory. Using the launch velocity errors and the uncertainty of the model, a spacecraft controlled by the proposed guidance law tracks the reference trajectory. Furthermore, the original Riccati differential equation in the receding horizon control algorithm is replaced by an equivalent convenient form of the Riccati differential equation that is based on the generating function. The high-efficiency solution of the equivalent equation avoids the online direct integration of the original Riccati differential equation, which significantly increases the computational efficiency for the receding horizon control problem. Numerical simulations using a nonlinear bicircular four-body model demonstrate the capabilities of the proposed receding horizon guidance law for the transfer mission. In addition, the generating function method improves the computational efficiency by at least one order of magnitude over the backward sweep method in solving the receding horizon control problem.     

Solitary,explosive and periodic solutions for electron acoustic solitary waves with non-thermal hot ions

A theoretical investigation has been made for electron acoustic waves propagating in a system of unmagnetized collisionless plasma consists of a cold electron fluid and ions with two different temperatures in which the hot ions obey the non-thermal distribution. The reductive perturbation method has been employed to derive the Korteweg–de Vries equation for small but finite amplitude electrostatic waves. It is found that the presence of the energetic population of non-thermal hot ions δ, initial normalized equilibrium density of low temperature ions μ and the ratio of temperatures of low temperature ions to high temperature ions β do not only significantly modify the basic properties of solitary structure, but also change the polarity of the solitary profiles. At the critical hot ions density, the KdV equation is not appropriate for describing the system. Hence, a new set of stretched coordinates is considered to derive the modified KdV equation. In the vicinity of the critical hot ions density, neither KdV nor modified KdV equation is appropriate for describing the electron acoustic waves. Therefore, a further modified KdV equation is derived. An algebraic method with computerized symbolic computation, which greatly exceeds the applicability of the existing tanh, extended tanh methods in obtaining a series of exact solutions of the various KdV-type equations, is used here. Numerical studies have been reveals different solutions e.g., bell-shaped solitary pulses, singular solitary “blowup” solutions, Jacobi elliptic doubly periodic wave, Weierstrass elliptic doubly periodic type solutions, in addition to explosive pulses. The results of the present investigation may be applicable to some plasma environments, such as Earth’s magnetotail region.     

随机多变量连续系统的辨识

本文讨论了由随机向量微分方程描述的多变量连续系统的辨识问题。虽然假定微分方程中的噪声项是白噪声,但是用方脉冲函数积分微分方程后,所得到的代数方程中的噪声项却是有色噪声。为了得到无偏的参数估计,本文提出一种预滤波方法。仿真例子表明这种方法是很有效的。     

Multi-objective optimal formation reconfiguration with scalarization and stochastic boundaries

A new method is proposed for multi-objective optimal control of satellite formation reconfiguration. First, necessary optimality conditions of single-objective fuel-optimal reconfiguration are studied. Problem of initial guess of the costate and control discontinuity are addressed by using a new developed stochastic based method. Terminal conditions of Two Point Boundary Value Problem (TPBVP) of the optimal control equations are relaxed using augmented Gaussian. Variances or bandwidths of the terminal conditions are set to decision variables to be optimized. Responses of varied initial costate are modeled and propagated using Riccati equation. Using quadratic convolution of relaxed terminal conditions, objective function of the fuel-optimal problem is reformulated into a quadratic equation, and a gradient based optimizer is used to search the optimal initial guess. In the second part of the paper, optimal control of format flying reconfiguration with two competitive objectives:fuel and time of flight (ToF) is studied. With Pascoletti-Serafini scalarization, the Multi-objective Optimal Control Problem (MOCP) is scalarized into a set of weighted Single-objective Optimal Control Problems (SOCP). Hamiltonian and switch function for the weighted SOCP are developed and constructed. New developed stochastic based approach is then used to search the optimal Pareto solutions. Verification and performances of proposed algorithms are demonstrated through numerical simulations of single-objective and multi-objective optimal control of low-thrust elliptical orbit formation reconfiguration. Potential applications of the algorithms to some other preliminary space mission designs are also discussed.     

低副机构的一种平面五级组运动分析

采用复数与矢量结合的方法,对低副机构中一种8杆12副的平面五级组运动分析问题进行了求解.首先根据不同的环路列出矢量方程,然后将方程转化成复数形式,经过代数消元后得到了输入输出方程.对次数的判定后证明方程的解为42个,即该平面五级组可以有42种构型.最后对一个实例进行了计算,给出了4组实根值.推导过程中的符号运算采用计算机代数系统"Mathematica".     

一种高超声速飞行器的非线性再入姿态控制方法

针对高超声速飞行器的再入非线性动力学模型,利用SDRE（state dependent Riccati equation）设计姿态控制器。基于奇异摄动理论,把姿态动力学分解成姿态角和姿态角速度跟踪内、外环回路,同时把非线性动力学伪线性化。每个跟踪回路用SDRE获得控制律,考虑到SDRE局部渐近稳定的特点,可以保证系统闭环稳定。最后设计高超声速飞行器飞行控制系统,并在高超声速条件下进行仿真,验证了该方案的有效性。