折线振子天线的时域有限差分法分析

研究了在手机和掌上电脑等移动通信终端上应用的平面折线振子天线．用三维变换方向隐式时域有限差分法（ADI—FDTD）进行仿真。利用各向异性参数得到改进的ADI—FDTD方法,并给出计算公式便于提高计算效率和精度,对新算法的稳定性、色散关系及参数的选取进行了研究。然后用这种方法研究了平面折线振子天线的反射系数、电压驻波比和带宽等特性．得到一些有参考价值的结论。     

固体中脉冲超声波传播的有限差分模拟

脉冲超声波传播过程的模拟研究对实际超声检测的信号分析具有重要的借鉴意义。对弹性固体中的声场方程应用时域有限差分法,进行了二维声场特性研究。通过完全匹配层方法对截断边界的吸收处理,解决了无限自由空间的计算问题。模拟计算表明该方法可用于各向异性介质的分析。分析了液固耦合情况下的超声波传播特性,取得了与光弹法物理试验相一致的结果。通过与平底孔试块实际缺陷回波信号的对比,验证了脉冲超声波有限差分法模拟计算的有效性和正确性。     

有效修正ADI-FDTD算法中数值色散的方法

针对变换方向隐式(ADI)时域有限差分(FDTD)算法提出了一种改进方法．通过添加各向异性参数来修正相速度误差，从而减少数值色散，提高计算精度。给出了改进后的ADI—FDTD计算公式，并对改进算法的稳定性、色散关系以及各向异性参数的选取进行了研究，数值结果证明了所给方法的有效性。     

聚束SAR扩展Chirp Scaling成像算法.

在合成孔径雷达（SAR Synthetic Aperture Radar)的成像算法中，Chrip Scaling成像算法具有计算效率高的优点，因此得到了较为广泛的应用。详细研究了子孔径扩展Chirp Scaling算法在高分辩率聚束模式SAR中的应用，包括子孔径划分和方位向处理问题，针对A.Moreira等1996年所提算法在处理聚束SAR数据时所产生的问题，给出了经过改进的适合于大斜视角处理的整个计算过程的完整表达式。在给出点目标仿真的同时，利用E-SAR实际数据对述方法进行了验证，结果对具体的兼容条带和聚束两种工作模式的SAR处理机设计具有一定的参考价值。     

一种改进的LMS算法及其在自适应消噪中的应用

介绍了一种改进的LMS自适应滤波算法及其在自适应噪声对消中的应用。该算法利用误差信号的相关值去调节步长,解决了算法收敛时间与稳态误差间的矛盾,为实际应用提供了更大的灵活性。它同时实现了均方误差小和收敛速度快,并且降低了LMS算法对噪声的敏感性。文中最后给出了仿真结果,仿真结果与理论分析是一致的。     

多层快速多极子算法的改进措施

 为精确求解散射问题,采用混合场积分方程（CFIE）、多层快速多极子算法（MLFMA）和共轭梯度算法（CG）的收敛技术。基于传统多层快速多级子算法,详细研究了二维拉格朗日插值节点数对计算精度的影响,并改进了插值方法,在不同的层采用不同的插值节点数;提出了在不同的层采用不同的精度控制来计算多级子模式数;分析了稀疏矩阵的对称性对内存使用的影响以及磁场积分方程对迭代初始值的选择。数值计算结果表明以上改进可较大幅度地提高计算精度和计算效率,同时降低内存使用,可满足复杂目标电磁散射计算要求。     

基于树形马氏链模型的可靠性分析方法

复杂系统的极限状态函数非线性程度较高，在进行可靠性分析时，易导致失效概率的计算误差大、效率低，针对上述问题，提出了树形马氏链（TMC）算法和基于该算法的可靠性分析方法。树形马氏链是对原始马尔可夫链的改进，其状态转移过程更加灵活，具备局部多链并行和自适应探索失效域边界的特性。树形马氏链通过多候选状态点扩大对失效域信息的收集，生成能充分反映失效分布特征的样本，对该样本进行自适应核密度估计得到近似最优的重要抽样分布密度函数，从而提高计算结果的准确度。文末的数值算例和工程算例验证了算法性能，计算结果表明算法对设计点、抽样起点的位置不敏感，处理强非线性及复杂串联系统问题时，能在少样本量下得到相对高准确度的计算结果，且在样本量改变时，计算结果相对稳定可靠；工程算例给出了所提方法在实际问题下的效率，体现了所提方法的工程应用价值。     

基于元胞自动机的飞机低空突防航路动态优化

提出了一种新型的航路动态优化算法——多层元胞自动机算法，描述了应用元胞自动机搜索飞机动态优化航路的策略。利用元胞自动机的计算并行性，使该算法大大减少了计算时间，提高了动态优化的实时性。最后用计算机进行了仿真，取得了较好的仿真结果。     

基于RBFNN自适应混合学习算法的航天测控系统任务可靠性分配

为解决执行航天测控任务的各设备存在复杂的时空关联、可视与信息关联等动态约束关系,使得航天测控系统任务可靠性分配建模和分析极其困难,同时模型求解效率低的问题,提出了自适应混合学习算法的径向基神经网络建模方法.算法通过训练样本相关性矩阵的主成分分析确定网络隐含层初始节点数;在此基础上,利用梯度信息衰减因子改进了迭代过程中网络参数的梯度信息计算方式,避免了学习过程早熟的不足,且加快了迭代收敛速度.最后,通过采集航天测控系统输入-输出数据,将自适应混合学习算法应用于参数训练,并给出了具体实现步骤.通过算例仿真,表明算法在解决航天测控系统任务可靠性分配问题时具有较高泛化能力和分配结果稳定等优点.      

基于进化策略的飞行控制系统多目标优化

针对传统飞机飞行控制律参数单目标优化设计不能同时满足多控制指标要求,且还有与飞行品质要求缺乏相关性和物理意义不明确等缺点,提出了一种基于改进多目标进化算法的飞行控制律多目标优化设计方法。该算法引入了决策者的偏好信息,用以指导算法的搜索过程,使算法在决策者感兴趣的区域进行搜索,不但缩小了算法的搜索空间,提高了算法的效率,而且一次运算只求得偏好区域内若干个折衷解,避免了决策者要在众多非劣解中作出困难的选择。在算法中,采用了新的偏好信息给定方法,并提出了一种模拟"投票选举"的偏好信息处理方法,该方法直观简便并易于实现。最后,使用改进的多目标进化算法对某型飞机纵向控制律设计进行数值仿真,结果显示该算法可有效提高控制律优化调参效率,结果满足期望的飞行品质要求。     

An SOR Implicit Time-accurate Scheme in Calculating Unsteady Flows

A new time-accurate marching scheme for unsteady flow calculations is proposed in the present work. This method is the combination of classical Successive Over-Relaxation (SOR) iteration method and Jacobian matrix diagonally dominant splitting method of LUSGS. One advantage of this algorithm is the second-order accuracy because of no factorization error. Another advantage is the low computational cost because the Jacobian matrices and fluxes are only calculated once in each physical time step. And, the SOR algorithm has better convergence property than Gauss-Seidel. To investigate its accuracy and convergency, several unsteady flow computa- tional tests are carried out by using the proposed SOR algorithm. Roe’s FDS scheme is used to discritize the inviscid flux terms. Un- steady computational results of SOR are compared with the experiment results and those of Gauss-Seidel. Results reveal that the numerical results agree well with the experimental data and the second-order accuracy can be obtained as the Gauss-Seidel for unsteady flow computations. The impact of SOR factor is investigated for unsteady computations by using different SOR factors in this algorithm to simulate each computational test. Different numbers of inner iterations are needed to converge to the same criterion for different SOR factors and optimal choice of SOR factor can improve the computational efficiency greatly.     

基于高阶耗散紧致格式的GMRES方法收敛特性研究

计算效率较低是当前限制高阶精度计算方法应用的重要因素。为了提高高阶精度混合型耗散紧致格式（HDCS）的计算效率,发展了适合多块对接网格的广义最小残值（GMRES）方法,并利用GMRES方法开展了HDCS格式的加速收敛研究。首先研究了GMRES的预处理方法、CFL数和内层迭代步数对HDCS数值模拟收敛特性的影响,计算结果显示：点松弛方法是一种高效的预处理方法;CFL数对计算收敛速度影响较大;GMRES方法存在最优的内层迭代步数。利用GMRES方法完成了NACA 0012翼型绕流、NLR 7301翼型绕流和DLR-F4翼身组合体绕流的数值模拟,并与其他隐式时间推进方法进行了对比,GMRES方法计算更加稳定,并且计算效率相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法可以提高5倍以上。研究结果表明,本文发展的GMRES方法在多块对接网格中具有良好的计算稳定性,计算结果的残差可以收敛到更低的量级,并且可以较大幅度地提高高阶精度数值模拟的计算效率。     

一种变步长LMS自适应滤波算法研究

为了解决最小均方(LMS)算法的稳定性以及收敛速度和稳态误差之间的矛盾,提出了一种复数型自适应变步长的LMS算法,它的权系数的调整取决于误差曲面在新权值点上的梯度。分析了新算法的收敛特性以及参数选择对算法性能的影响。算法具有较快的收敛速度、鲁棒稳定性且运算小易于实现的特点。计算机仿真的结果与理论分析相一致,证实了算法的收敛性能优于标准的LMS算法并且具有较好的实用性。     

叶轮机械全三维气动杂交问题的数值求解方法

本文介绍了一种求解叶轮机械全三维气动杂交问题新的数值方法。本方法使用了非正交曲线座标系和非正交速度分量的势函数方程, 并构造了守恒形式的近似分解因式格式, 通过不断修正叶片通道宽度, 使叶面上速度达到预先要求的速度分布。根据本方法使用 FORTRANTT编写了计算程序, 并对单转子压气机叶片进行了修正计算, 其结果十分令人满意。      

大范围收敛的摄动Lambert问题新型解法:拟线性化-局部变分迭代法

研究实时、高效、稳定性强的高性能空间轨道计算方法对于中国未来航天工程具有重大应用价值。针对强非线性系统的多维两点边值问题,提出了一种拟线性化-局部变分迭代法（QL-LVIM）,通过拟线性化（QL）思想,将非线性两点边值问题转化为一系列具有一定迭代格式,并且成对出现的初值问题,进而通过局部变分迭代法（LVIM）对其进行求解。利用拟线性化的大范围收敛特性和局部变分迭代法的快收敛、高精度特性,该方法能够在较大的时间和空间尺度下快速精确获得摄动Lambert问题的初速度和转移轨道,其收敛域远大于传统的牛顿打靶法,为航天器轨道转移提供了一种简便高效、稳定性强的新型计算方法。在不同轨道情形下,与几类参考方法对比,结果表明本方法能够在计算效率方面实现大幅提升,并且能够在大范围内实现快速收敛。方法的有效性在地-月系三体问题中得到了进一步验证。     

考虑性能及成本的固体火箭发动机多学科设计优化

为综合考虑固体火箭发动机的燃烧室、药柱、内弹道、喷管及成本等学科影响,梳理学科间耦合关系,并建立了以总冲最大、成本最小为优化目标的固体火箭发动机多学科设计优化（MDO）模型。为降低MDO问题的计算成本,提出一种基于Kriging代理模型的多目标自适应优化方法（KRG-MAOM）。优化过程中,分别对目标与约束构建Kriging模型,并采用多目标优化算法求解,在伪Pareto解中综合考虑支配关系与分布特性选取新增样本点,引导优化快速收敛。算例结果表明,KRG-MAOM算法在全局收敛性与优化效率方面具有显著优势。最后,采用KRG-MAOM算法求解该MDO问题,得到可行的Pareto解集方案,与初始方案相比,同性能情况下成本节省约3.36%;同成本情况下性能提升约10.93%,从而验证MDO模型合理性与KRG-MAOM算法有效性。     

基于凸优化的再入轨迹三维剖面规划方法

可重复使用飞行器一般采用大升阻比气动外形，再入轨迹三维剖面规划方法可充分发挥这类飞行器固有的机动能力。计算量大是制约三维剖面规划应用的难题，为提高计算效率，提出了一种基于凸优化的再入轨迹三维剖面规划方法。首先，分析运动方程特性，利用定义新的控制变量、约束松弛、连续线性化等技术，将原始非凸的三维剖面规划问题转化为一个凸优化问题。其次，将指令反解步骤嵌入至序列凸化算法中，通过迭代求解凸优化子问题，获得原问题的可行解。数值仿真结果表明所提方法具有较高的求解精度和确定的收敛性质，飞行器的机动能力得到充分发挥；与伪谱法的结果对比表明凸优化方法在轨迹规划问题上具有更高的求解效率。     

动态网格快速生成及组合体跨音速非定常气动力计算

采用代数方法、局部变形技术及快速动态弹性变形技术生成三维翼/体组合的贴体动态网格。用时间精确近似因式分解差分方法求解守恒型非定常全位势方程。考虑附面层影响，通过粘/位迭代得到跨音速翼/体的非定常气动力。所得结果与实验数据相比相当吻合。     

A 3D hybrid grid generation technique and a multigrid/parallel algorithm based on anisotropic agglomeration approach

A hybrid grid generation technique and a multigrid/parallel algorithm are presented in this paper for turbulence flow simulations over three-dimensional (3D) complex geometries. The hybrid grid generation technique is based on an agglomeration method of anisotropic tetrahedrons. Firstly, the complex computational domain is covered by pure tetrahedral grids, in which anisotropic tetrahedrons are adopted to discrete the boundary layer and isotropic tetrahedrons in the outer field. Then, the anisotropic tetrahedrons in the boundary layer are agglomerated to generate prismatic grids. The agglomeration method can improve the grid quality in boundary layer and reduce the grid quantity to enhance the numerical accuracy and efficiency. In order to accelerate the convergence history, a multigrid/parallel algorithm is developed also based on anisotropic agglomeration approach. The numerical results demonstrate the excellent accelerating capability of this multigrid method.     

Real-time trajectory planning for UCAV air-to-surface attack using inverse dynamics optimization method and receding horizon control

This paper presents a computationally efficient real-time trajectory planning framework for typical unmanned combat aerial vehicle (UCAV) performing autonomous air-to-surface (A/S) attack. It combines the benefits of inverse dynamics optimization method and receding horizon optimal control technique. Firstly, the ground attack trajectory planning problem is mathematically formulated as a receding horizon optimal control problem (RHC-OCP). In particular, an approximate elliptic launch acceptable region (LAR) model is proposed to model the critical weapon delivery constraints. Secondly, a planning algorithm based on inverse dynamics optimization, which has high computational efficiency and good convergence properties, is developed to solve the RHCOCP in real-time. Thirdly, in order to improve robustness and adaptivity in a dynamic and uncer- tain environment, a two-degree-of-freedom (2-DOF) receding horizon control architecture is introduced and a regular real-time update strategy is proposed as well, and the real-time feedback can be achieved and the not-converged situations can be handled. Finally, numerical simulations demon- strate the efficiency of this framework, and the results also show that the presented technique is well suited for real-time implementation in dynamic and uncertain environment.