应用GPS载波确定载体姿态

本文简述了利用GPS载波相位技术进行载体姿态确定的原理，比较了目前已有的多种整周模糊度求解算法，选用最小二乘法进行了载体定姿实践，减少了备选整周模糊度的组合数，并采用多种约束信息来剔除不正确的模糊度组合，结合实验，指出了载体姿态变化和整周模糊度搜索速度的关系，提出了一种模糊度确定的辅助方法，并给出了应用算例。     

A closed-form solution for moving source localization using LBI changing rate of phase difference only

Due to the deficiencies in the conventional multiple-receiver localization systems based on direction of arrival(DOA) such as system complexity of interferometer or array and amplitude/phase unbalance between multiple receiving channels and constraint on antenna configuration,a new radiated source localization method using the changing rate of phase difference(CRPD)measured by a long baseline interferometer(LBI) only is studied. To solve the strictly nonlinear problem, a two-stage closed-form solution is proposed. In the first stage, the DOA and its changing rate are estimated from the CRPD of each observer by the pseudolinear least square(PLS) method,and then in the second stage, the source position and velocity are found by another PLS minimization. The bias of the algorithm caused by the correlation between the measurement matrix and the noise in the second stage is analyzed. To reduce this bias, an instrumental variable(IV) method is derived. A weighted IV estimator is given in order to reduce the estimation variance. The proposed method does not need any initial guess and the computation is small. The Cramer–Rao lower bound(CRLB) and mean square error(MSE) are also analyzed. Simulation results show that the proposed method can be close to the CRLB with moderate Gaussian measurement noise.     

基于试验的电动舵机模型辨识

目前电动舵机模型的建立基本上是基于理论的数学推导,所建立的舵机模型与实际舵机模型存在一定的误差。针对此问题,提出一种基于试验实测数据建立舵机模型的方法,搭建自动数据采集试验平台,采用组合窗法对实测数据进行频域转换,得到舵机的频率响应函数并采用非线性最小二乘法辨识舵机模型。对辨识得到的数学仿真结果与实际舵机响应比较,在频域,整个频段内幅频特性最大绝对误差1.9 dB;在时域,仿真模型跟踪实际舵机角度的最大误差1.5°,延迟时间10 ms。结果表明,基于试验数据建立的舵机仿真模型能很好地模拟实际舵机。     

脉冲随机微分系统的均方指数稳定性分析

研究了脉冲随机时滞微分泛函方程的均方指数稳定性问题。利用Lyapunov-Razumikhin型方法及随机分析的一些技巧,建立了一类脉冲随机泛函微分方程的均方指数稳定性定理。     

总体最小二乘算法在自适应信道均衡中的应用研究

分析了最小二乘解产生偏差的原因,推导出了最小二乘解偏差项的数学表达式。提出了一种基于瑞利商的变步长总体最小二乘求解方法。仿真试验表明,总体最小二乘算法的性能优于加权递归最小二乘算法,是一种适用于输入信号端含有噪声的自适应信道均衡算法。     

Eliminating negative VTEC in global ionosphere maps using inequality-constrained least squares

Currently, ground-based Global Navigation Satellite System (GNSS) stations of the International GNSS Service (IGS) are distributed unevenly around the world. Most of them are located on the mainland, while only a small part of them are scattered on some islands in the oceans. As a consequence, many unreasonable zero values (in fact negative values) appear in Vertical Total Electron Content (VTEC) of European Space Agency (ESA) and Center for Orbit Determination in Europe (CODE) IONEX products, especially in 2008 and 2009 when the solar activities were rather quiet. To improve this situation, we directly implement non-negative physical constraints of ionosphere for global ionosphere maps (GIM) with spherical harmonic functions. Mathematically, we propose an inequality-constrained least squares method by imposing non-negative inequality constraints in the areas where negative VTEC values may occur to reconstruct GIM models. We then apply the new method to process the IGS data in 2008. The results have shown that the new algorithm efficiently eliminates the unwanted behavior of negative VTEC values, which could otherwise often be seen in the current CODE and ESA GIM products in both middle and high latitude areas of the Southern Hemisphere (45°S∼90°S) and the Northern Hemisphere (50°N∼90°N). About 64% of GPS receivers’ DCBs have been significantly improved. Finally, we compare the GIM results between with and without the inequality constraints, which has clearly shown that the GIM result with inequality constraints is significantly better than that without the inequality constraints. The inequality-constrained GIM result is also highly consistent with the final IGS products in terms of root mean squared (RMS) and mean VTEC.     

贴壁方柱湍流场TR—PIV实验研究

采用时变粒子图像速度场测试技术（TR—PIV），对低速循环水槽中贴壁二维方柱绕流湍流场进行了细致的测量。通过对实验得到的30000个连续瞬态速度场进行分析，得到了时均速度场和流线图谱，以及流向、法向速度分量的脉动强度场和涡量场。此外，对脉动速度时序信号进行谱分析，得到了流场中低频大尺度相干结构脱落频率。通过对瞬态速度场的分析，揭示了贴壁方柱绕流的旋涡脱落及其发展过程。     

无网格法耦合RNG k-ε湍流模型在亚、跨声速翼型黏性绕流中的数值模拟

基于移动最小二乘无网格方法,耦合RNG(Re-Normalisation Group)k-ε湍流模型求解雷诺平均Navier-Stokes方程。采用AUSM(Advection Upstream Splitting Method)+-up迎风格式求解数值通量,应用在高度各向异性点云结构中取得良好结果的点云重构技术结合移动最小二乘法拟合空间导数,并用三阶SSP(Strong Stability Preserving)型Runge-Kutta显式时间推进格式求解离散后的控制方程。在此基础之上,实现了对NACA0012、RAE2822翼型亚、跨声速黏性绕流的数值模拟,给出了翼型表面压力系数分布曲线、不同位置处的平均速度剖面、马赫数等值线等计算结果,并与实验值及相关文献数值模拟结果进行比较,结果吻合较好。表明所发展的结合点云重构技术的无网格方法耦合RNGk-ε湍流模型能够成功模拟翼型亚、跨声速黏性绕流,验证了所提算法的有效性,并拓展了无网格方法求解湍流流动的途径。     

CamShift在加油锥套识别跟踪中的应用

针对插头锥套式自主空中加油过程中加油锥套的识别跟踪问题,提出了一种基于Cam Shift的视觉识别跟踪算法。选取加油锥套跟踪视觉特征,将RGB空间图像转换到HSV色彩空间,并采用自适应阈值算法分割图像。运用最小二乘椭圆拟合算法解算锥套在图像中的中心点坐标;根据目标像素位置变换速度,对Mean Shift算法的目标搜索窗宽进行自适应修正;定义基于巴氏系数的阈值函数,处理锥套被遮挡的情况;将VC与Open CV相结合,设计了加油锥套视觉跟踪软件,对锥套位置进行实时解算和显示。试验结果表明,该算法可达到毫米级精度,满足自主空中加油视觉导航的技术要求。     

基于粒子群算法的极短弧定轨

针对经典的初轨计算方法在极短弧定轨中不适用的情况,建立了一种基于粒子群算法的极短弧(TooShort-Arc,TSA)定轨的计算方法。该方法将问题转化为两个三变量的分层优化问题,采用(a,e,M)作为优选变量,在保持问题维数较低的同时,实现了计算结果和观测资料的解耦。由于实测资料处理中的野值剔除方法不适用于粒子群算法,所以,采用稳健估计法,通过在适值函数中使用最小中值二乘准则,实现了稳健的极短弧计算方法。同时,应用MATLAB计算软件,选用缺省参数实现该算法,以进行数据验证。基于实测数据的数值验证表明,方法对于近圆轨道目标30s以下的弧段仍可以获得有效的结果,10s弧段误差仅为16km。此精度满足后续处理的需要,且方法稳健,具有很高的崩溃点。