异常值和未知观测噪声鲁棒的非线性滤波器

针对含异常观测值的非线性系统滤波问题,以Huber损失函数替代推导滤波器最大后验准则中观测误差的l₂范数,构造出了一种新的优化准则函数,从而给出了一种对异常值鲁棒的非线性后验线性化滤波器。分析表明：由于Huber损失函数兼具l₁和l₂范数的性质,从而使得由这个新准则推导出的滤波器,不仅具有l₂范数的低误差拟合性,也具备l₁范数对异常值的鲁棒性。而当观测噪声的分布未知时,通过引入箱线图法检测异常值,并对噪声统计分布的参数进行估计,进一步提出了对异常值和未知观测噪声分布鲁棒的非线性后验线性化滤波器。仿真实验验证了分析结果的有效性,并表明本文算法的性能优于现有文献报道的非线性滤波算法。     

基于H∞算法的飞机机翼结冰气动参数辨识

针对国内大型飞机结冰防护需求,开展针对大型结冰研究样机的H_∞算法参数辨识结冰探测研究。首先通过参数调节选取一组合适的H_∞算法参数,利用考虑测量噪声的结冰研究样机飞行仿真数据验证H_∞算法的辨识能力,由结果对比发现辨识算法能够跟踪飞机气动导数随结冰累积过程的变化趋势,辨识精度较高,其最大归一化平方根（RMS）误差仅为真值的11%;分析了H_∞算法对81种不同结冰累积过程的辨识能力,通过结果分析发现结冰累积时间较长且结冰速度较慢的情况辨识效果较差,结冰累积时间在100~300 s之间辨识精度较高;最后利用蒙特卡罗仿真分析了不同测量噪声大小对H_∞算法辨识精度和跟踪延时的影响,给出了3个纵向气动导数在随机误差影响下的辨识误差和跟踪延时的统计结果,发现在给定噪声标准差变化范围内,升力和俯仰力矩关于迎角的导数能够得到较为准确的辨识结果,二者的归一化平方根误差均值仅为各自真值的1.8%和4%,其预报延时均值最大仅为3 s和9.5 s。     

基于隶属度和LMK-ELM的航空电子部件诊断方法

为提高航空电子部件模块级故障诊断精度，提出一种新的面向"软聚类"的局部多核学习（LMKL）-超限学习机（ELM）离线诊断方法。通过引入模糊C均值聚类对样本空间进行模糊划分，挖掘聚类内部多样性的同时，实现了对过学习的抑制；将模糊划分产生的隶属度信息融入LMKL-ELM的优化过程，运用基于初始-对偶混合优化问题的三步优化策略克服了局部核权重二次非凸的问题，在l₁-范数与l₂-范数约束下分别给出了相应的更新方法。将所提方法应用于某型机前端接收机，结果表明：与4种流行的多核诊断方法相比，该方法可有效避免漏警、抑制虚警，在l₁-范数和l₂-范数约束下，其诊断精度比其他方法的平均值分别提升了4.09%和5.13%。     

靶场光电经纬仪跟踪精度评估技术

靶场鉴定光测设备跟踪精度以往采用静态拍星和动态误差增量合成法,然而这种方法与执行飞行试验任务时的真实工作条件差别较大,鉴定的内容和结果都有一定的局限性。针对光电经纬仪设备的实际跟踪状态,以GPS(全球定位系统)测量信息作为比较标准,采用散度、偏度、精度来综合评估设备跟踪精度,通过对各通道测量数据进行散度、偏度分析,建立散度统计模型检验方法和偏度统计处理算法,系统评估与分析光测设备跟踪质量,实现对光电经纬仪跟踪精度的综合评估,并利用实际数据进行仿真计算和分析,仿真分析结果验证了其方法的可行性。     

一种管道测绘系统螺旋误差精密补偿方法

针对石油/天然气管道测绘系统的测量精度需求,分析了测量数据中存在的螺旋误差对参数精度的影响,研究了测量过程中产生螺旋误差的机理,推导了螺旋误差的数学模型;在此基础上,提出了一种螺旋误差补偿方法,进而给出了模型参数的标定和补偿方案。实际的管道测量验证试验结果表明,本误差补偿方法,可以基本消除测量数据中的螺旋误差,有效提高管道参数的测量精度。     

复杂曲面铝反射镜磁流变抛光工艺优化

采用超精密车削加工的复杂曲面铝反射镜只能满足红外光学系统的应用需求,若要满足更高需求的应用场合,需要进一步提升反射镜面形精度。磁流变抛光能够进行确定性修形,在复杂曲面加工中具有独特优势,但是复杂曲面连续变化的面形特征,在磁流变抛光时会导致去除函数不稳定,影响误差收敛效率和加工精度。从高精度复杂曲面铝反射镜的应用需求出发,提出了复杂曲面局部区域磁流变抛光去除函数的动态建模方法,给出了驻留时间求解算法,以平均曲率变化最小为原则,设计了抛光路径优化算法,针对该算法计算速度慢的问题,提出了优化策略,并通过试验进行了验证,最终加工的复杂曲面铝反射镜的面形误差为0.216λ PV、0.033λ RMS （λ=632.8 nm）。     

经纬仪横轴与竖轴垂直度误差分析与检定

分析了经纬仪横轴与竖轴垂直度误差的产生原因及该误差对水平方向观测值的影响，给出了该误差的检定方法，对高低点法和高平点法检定精度进行了分析，并以实验进行了验证。     

CFD计算网格误差分析的一个算例

为了给一般工程问题的CFD仿真的网格划分提供参考,研究了光滑直圆管中粘性稳态等温流动在不同尺寸、不同纵横比网格情况下的CFD计算数值解的误差.采用FLUENT软件进行数值计算,应用CFD仿真验证和确认的基本方法,与精确解比较,计算了各种网格的数值计算结果的误差,得到了网格独立解的网格尺寸,给出了误差随网格尺寸变化的规律.     

基于广义Richardson外插方法的阻力预测精度分析

利用计算流体力学(CFD)方法预测阻力是飞行器气动设计中的关键环节。采用广义Richardson外插方法分别对数值方法预测二维简单构型的压差阻力、摩擦阻力和三维复杂构型的总阻力的精度进行了分析。在NACA 0012翼型无黏绕流和平板湍流边界层两个算例中验证了NSAWET程序和广义Richardson外插方法,分别得到数值算法预测压差阻力和摩擦阻力能达到的名义精度。进而模拟三维通用研究模型(CRM)翼身组合体绕流,得到的阻力名义精确值在DPW 5的统计误差带范围之内;综合DPW 5的计算结果来看,不同CFD解算器的结果之间存在一定差别,阻力预测精度总体上不符合二阶。可见,标准Richardson方法采用的二阶精度假设难以普遍适用,有必要采用广义Richardson外插方法得到名义精度。针对不合理的名义精度,采用Roache建议的方法加以限制。广义Richardson外插方法有助于提高误差分析的合理性,可以进一步降低网格对阻力预测的影响。     

针栓式喷注单元雾化角模型分析

为了实现不同径向孔形的针栓式喷注器雾化角的准确预测，从动量守恒方程出发建立了液膜撞击液膜和液膜撞击液束的雾化角理论修正模型。对于液膜撞击液膜的喷注单元，模型中通过理论推导引入了2个变形因子，将撞击的几何变形效应与雾化角关联；对于液膜撞击液束，通过引入阻塞率定义有效撞击动量比，同时将液束入口孔形的影响隐含考虑在变形因子中，最后根据高速摄影试验结果和数值仿真结果获得了对应的变形因子组合系数，使得新的雾化角模型适应性更广、准确性更高。结果表明：引入变形因子和阻塞率的理论模型预测值与试验及数值仿真结果吻合很好；对于液膜撞击液膜，变形因子基本维持在0.9~1.1，根据试验结果及仿真结果，变形因子推荐值为C₁=0.99和C₂=1.06；对于液膜撞击液束，变形因子推荐值为C₁=0.75和C₂=1.25。该模型根据实际出口的轴向动量和合成总动量计算雾化角，隐含考虑了撞击作用造成的影响，较根据撞击前入口的轴向动量和合成总动量计算雾化角的常用模型预测值准确度显著提高，为针栓式喷注器的理论研究和工程设计提供了重要参考。     

样本量为2的极小样本相容性检验方法

在航空航天领域由于成本、时间周期等原因进行疲劳寿命及可靠性评估时样本量通常极少（m=1或2）,利用相容性检验方法可对样本量进行扩充。常规的Wilcoxon秩和检验和K-S（Kolmogorov-Smirnov）检验适用于小样本情形,而极小样本相容性检验方面研究较少,且缺乏对方法合理性的详细说明和对不同方法检验功效优劣的比较。航空航天产品疲劳寿命多服从正态分布,因此本文主要以正态分布作为研究对象。利用Monte Carlo仿真发现从某一正态分布N（μ,σ²）中随机抽取两个样本x₁、x₂计算均值μ₁和标准差σ₁后构建新正态分布N（μ₁,σ₁²）,其±σ₁、±2σ₁和±3σ₁范围内的点落在原正态分布N（μ,σ²）±3σ范围内的概率依次为99.80%、98.13%和97.37%。在此基础上针对现场试验数据样本量为2的情况,本文提出利用3σ原则对先验信息数据进行相容性检验从而扩充样本量的方法。将该方法与两种文献方法对比后发现其误差率明显更低并呈现出检验性能随先验数据增加而不断提高的优势。     

边界层内WENO格式特征重构精度的理论分析

WENO格式因其高精度以及良好的激波捕捉能力而被广泛地用于可压缩湍流的数值模拟过程中。为了在有强间断的流场中仍然能保证数值稳定，通常将流动物理量投影到特征空间后再进行WENO重构。然而这种特征重构方法在边界层内的精度问题仍需要进一步分析与研究。通过可压缩边界层自相似解构造平板边界层内法向的流动参数，并在该方向上进行WENO重构，得到求解半点数值通量的WENO非线性权重。通过分析实际得到的权重与理想权重的偏差大小发现在边界层流动中特征重构的精度明显低于分量重构。为了更深入地分析特征重构精度降低的原因，通过对曲线坐标系下的流动方程进行理论分析，推导得到特征重构时的特征变量形式。根据特征变量表达式以及自相似解中权重误差分布发现引起特征重构精度降低的主要原因是特征投影过程在半点处出现了额外的极值点，同时这种因特征投影过程而产生极值的特性将不会随着左右特征矩阵的不同选取而发生改变。进一步地，基于理论分析的普适性，对于任何光滑流场，只要使用特征重构的投影过程，都会在半点处产生极值从而使得WENO系列格式产生精度降低的现象。     

基于特征模型的Halo轨道维持控制

平动点，尤其是共线平动点轨道在未来深空探测活动中具有重要的应用价值，但由于共线平动点轨道不稳定，运行在其上的航天器在无控情况下将很快偏离标称轨道，因此在实际任务中，轨道维持必不可少。针对地月系L₂点附近的Halo轨道维持问题，首先在圆型限制性三体模型下，利用Richardson三阶近似解析解、微分修正以及打靶法获得了用于维持控制的标称轨道；然后设计了基于特征模型理论的黄金分割控制器用于速度跟踪以及PD控制器用于位置跟踪；最后分别在圆型限制性三体模型和双圆限制性四体模型下进行了仿真分析。结果表明：在两种模型下，位置和速度的跟踪精度分别优于100 m和0.003 m/s，但双圆限制性四体模型下所需总的速度增量比圆型限制性三体模型下所需总的速度增量高一个数量级。     

高速流场中变刚度复合材料层合板颤振分析

变刚度复合材料层合板在高速流场中的颤振行为是设计中需要考虑的问题。本文研究了高速流场中的曲线纤维变刚度层合复合材料壁板非线性颤振响应,分析了边界条件和纤维方向对颤振特性的影响。利用von-Karman大变形应变-位移关系,采用气动力活塞理论,根据虚功原理和有限元法建立变刚度复合材料壁板颤振的气动弹性力学模型,采用Newmark法对壁板的颤振方程求解。给出了不同边界条件和纤维方向条件下层合复合材料壁板的颤振特性。计算结果表明：随着纤维在板中心处或在边界±a/2处与x方向夹角（T₀或T₁）的增大,颤振临界动压减小;相同动压下,随着T₀或T₁的增大,极限环振幅增大。研究表明采用曲线纤维进一步提高了复合材料层合板的可设计性,通过调整曲线纤维路径可以改变复合材料壁板的颤振特性。     

基于内圈滚道廓形的圆柱滚子轴承旋转精度

在轴承几何关系和运动关系分析基础上,建立了考虑内圈滚道廓形的圆柱滚子轴承径向跳动数学模型,分析了内圈滚道圆度误差幅值、谐波阶次、滚子个数以及径向游隙对轴承旋转精度的影响规律,并验证了模型的正确性.结果表明:内圈滚道的圆度误差对轴承旋转精度的影响较大,内圈滚道圆度误差幅值对轴承旋转精度影响较大，内圈滚道低阶谐波对轴承旋转精度的影响较高阶谐波显著；当内圈滚道谐波阶次为滚子个数的整数倍时，轴承旋转精度显著下降；减小内圈滚道圆度误差幅值能有效提高圆柱滚子轴承的旋转精度.      

复杂建筑物对近场扩散影响的数值与风洞模拟的比较分析

采用CFD数值模拟技术与风洞试验相结合的手段模拟了复杂建筑物群对周围流动与污染物扩散的影响。CFD技术采用k-ε（RNG）湍流模型模拟了建筑物群对流动与扩散的影响,风洞试验通过采用多种探测手段,分析建筑物群周围的流动和弥散,并将统计学方法应用于风洞试验结果,验证了CFD数值模拟结果的合理性。研究结果表明,CFD技术能较好地模拟建筑物群对污染物弥散的影响,并与风洞试验吻合,有、无建筑物B的情况下,吻合因子FAC2与FAC5均大于50％,归一化均方误差（NMSE）均小于4,部分偏差（｜FB｜）均小于0．3;建筑物对污染物弥散的影响非常复杂,建筑物尾流对污染物浓度分布影响较为显著,建筑物后方空腔区污染物的数值模拟结果略高于风洞试验结果。综合分析表明,风洞试验与数值模拟相结合的方法是研究这类微小尺度湍流扩散问题的有效手段。     

一种制导误差分离的新方法

针对制导误差分离模型中环境矩阵S存在严重病态性,从而影响分离结果精度问题,提出了一种基于动力系统求解的制导误差分离方法。该方法从分析线性迭代求解方法入手,将具有病态特性的线性方程组求解问题转化为对相应刚性动力系统的求解问题。这里给出了该方法收敛性及其他特性的证明。为了验证该方法效果,在遥外测视速度误差分别为0.01m/s、0.02m/s以及0.03 m/s的条件下,选用PB(Primary Bayesian,主成分贝叶斯)估计方法与其进行比较,数值结果表明,该方法可有效地降低环境矩阵病态性对误差分离结果的影响,且分离结果的稳健性和精度都优于PB估计方法得到的结果。