基于自适应步长选择的周期格型线搜索估计

稀疏、含噪观测条件下周期点过程的周期估计是一个经典的信号处理问题。针对该问题,提出了一种格型线搜索(LLS)算法,该算法通过数值方式搜索似然函数的最大值,但其性能取决于人为预先选取的搜索步长。推导了一个步长计算公式,并利用该公式改进了LLS算法。改进的LLS算法能够自适应选择搜索步长,其达到的克拉美-罗界(CRLB)的信噪比(SNR)门限与最大似然估计(MLE)算法一致,但计算复杂度比后者低一个多的数量级。性能分析与仿真实验表明,所提算法比已有算法能更好地实现估计精度与复杂度的折中。     

利用重要性采样的时差-频差联合估计算法

针对无源定位中参考信号真实值未知的时差(TDOA)-频差(FDOA)联合估计问题,构建了一种新的时差-频差最大似然(ML)估计模型,并采用重要性采样(IS)方法求解似然函数极大值,得到时差-频差联合估计。算法通过生成时差-频差样本,并统计样本加权均值得到估计值,克服了传统互模糊函数(CAF)算法只能得到时域和频域采样间隔整数倍估计值的问题,且不存在期望最大化(EM)等迭代算法的初值依赖和收敛问题。推导了时差-频差联合估计的克拉美罗下界(CRLB),并通过仿真实验表明,算法的计算复杂度适中,估计精度优于CAF算法和EM算法,在不同信噪比条件下估计误差接近CRLB。     

Pursuer identification and time-to-go estimation using passive measurements from an evader

We present an algorithm for identifying the parameters of a proportional navigation guidance missile (pursuer) pursuing an airborne target (evader) using angle-only measurements from the latter. This is done for the purpose of classifying the missile so that appropriate counter-measures can be taken. Mathematical models are constructed for a pursuer with a changing velocity, i.e., a direction change and a speed change. Assuming the pursuer is launched from the ground with fixed thrust, its motion can be described by a four-dimensional parameter vector consisting of its proportional navigation constant and three parameters related to thrusting. Consequently, the problem can be solved as a parameter estimation problem, rather than state estimation and we provide an estimator based on maximum likelihood (ML) to solve it. The parameter estimates obtained can be mapped into the time-to-go until intercept estimation results are presented for different scenarios together with the Cramer-Rao lower bound (CRLB), which quantifies the best achievable estimation accuracy. The accuracy of the time-to-go estimate is also obtained. Simulation results demonstrate that the proposed estimator is efficient by meeting the CRLB.     

X射线脉冲星导航硬件脉冲轮廓累积研究

为了在X射线脉冲星地面实验系统仿真源模拟产生X射线的基础上,能够快速稳定地得到脉冲轮廓,采用硬件历元叠加的方法获得脉冲轮廓。研究了用硬件实现历元叠加及其数据整合的算法,该算法首先在MATLAB现场可编程逻辑阵列(FPGA)中实现,再通过MATLAB硬件描述语言(HDL)代码生成模块把算法转换成HDL,经编译后获得配置硬件的Bit文件,最终在开发板FPGA上实现数据处理的硬件模块。一段时间内的光子到达时间数据通过MATLAB算法得到的脉冲轮廓数据与通过硬件模块处理后得到的数据结果存在误差,在单个时间窗口内误差最大值为2个光子数,误差平均值占光子数统计平均值的0.084%;两组统计的脉冲轮廓数据中不同数据占总数据个数的9.481%,这样的误差不影响后端模拟导航模块的导航。利用硬件实现的历元叠加及其数据整合模块具有处理速度快、设备紧凑、功耗低的特点,为航天器利用X射线脉冲星导航提供了一种可行的硬件数据处理技术上的支持。     

基于Starlink机会信号/INS的组合导航方法

针对全球卫星导航系统(GNSS)因频点单一、落地功率低、易受电磁干扰以及存在覆盖较差区域等潜在的被拒止或被干扰导致的导航系统性能降低甚至失效的问题,提出了一种基于星链(Starlink)机会信号融合惯性导航系统(INS)的飞行器动态组合导航方法。首先分析了星链信号体制,建立了基于星链星座卫星下行机会信号的瞬时多普勒定位观测模型,设计了一种基于频率细分的快速最大似然多普勒频率估计方法,然后建立了基于扩展卡尔曼滤波(EKF)的Starlink机会信号/INS的组合导航模型,并对该导航方法进行了实验及分析。结果表明,该方法可为飞行器提供长航时、连续、高精度的导航。动态飞行情况下,该方法可实现平均优于25 m的三维定位精度和平均优于0.1 m/s的速度估计精度,比相同观测时间下的惯导精度提高了1～2个数量级,显著提高了飞行器的导航精度,可为战略导航提供方法和技术支撑。     

基于改进ICPF的多分量LFM信号参数估计

针对线性调频信号参数估计中搜索匹配过程带来的计算量问题,提出一种基于改进ICPF的多分量线性调频信号参数估计方法。给出了一种基于非均匀快速傅立叶变换的ICPF改进算法,利用改进的ICPF算法估计线性调频信号的调频率,利用解线频调方法估计中心频率和幅度。方法避免了参数估计过程中的搜索匹配过程,计算复杂度低,在低信噪比条件下具有较好的估计性能。仿真结果验证了方法的有效性。     

Maximum Likelihood DOA Estimation in Unknown Colored Noise Fields

Direction-of-arrival (DOA) estimation in unknown noise environments is an important but challenging problem. Several methods based on maximum likelihood (ML) criteria and parameterization of signals or noise covariances have been established. Generally, to obtain the exact ML (EML) solutions, the DOAs must be jointly estimated along with other noise or signal parameters by optimizing a complicated nonlinear function over a high-dimensional problem space. Although the computation complexity can be reduced via derivation of suboptimal approximate ML (AML) functions using large sample assumption or least square criteria, nevertheless the AML estimators still require multi-dimensional search and the accuracy is lost to some extent. A particle swarm optimization (PSO) based solution is proposed here to compute the EML functions and explore the potential superior performances. A key characteristic of PSO is that the algorithm itself is highly robust yet remarkably simple to implement, while processing similar capabilities as other evolutionary algorithms such as the genetic algorithm (GA). Simulation results confirm the advantage of paring PSO with EML, and the PSO-EML estimator is shown to significantly outperform AML-based techniques in various scenarios at less computational costs.     

Multiple radar targets estimation by exploiting induced amplitude modulation

This work deals with the problem of estimating complex amplitudes, Doppler frequencies, and directions of arrival (DOA) of multiple targets present in the same range-azimuth resolution cell of a surveillance radar. The maximum likelihood (ML) and the asymptotic (large sample size) ML (AML) estimators are derived. To reduce the computational complexity of the maximization of the nonlinear two-dimensional criterion function of the AML estimator, we propose a computationally efficient algorithm based on the RELAXation method. It allows decoupling the problem of jointly estimating the parameters of the signal components into a sequence of simpler problems, where the parameters of each component are separately and iteratively estimated. The proposed method overcomes the resolution limitation of the classical monopulse technique and resolves multiple targets exhibiting an arbitrarily small difference in azimuth as long as their Doppler frequencies differ at least by the inverse of the number of integrated pulses, provided that enough signal-to-noise ratio (SNR) per pulse is available. The performance of the proposed AML-RELAX estimator is numerically investigated through Monte Carlo simulation and Cramer-Rao lower bound (CRLB) calculation.     

Small Aperture Angle Measurement for Active Echo Location Systems

The problem of bearing estimation for active systems is examined from the point of view of the generalized wideband ambiguity function (GAF). The maximum likelihood (ML) estimator is derived and its local and global properties are discussed. A structure is proposed which searches the three-dimensional ambiguity surface in two stages first, in range-Doppler, and then, in bearing with the goal of reducing search complexity when utilizing highly resolvent waveforms. Comparisons are made between the ML estimators and structures utilizing phase information which generate closed form estimators. The beneficial results of full bandwidth utilization are discussed in terms of both local and global properties of the GAF.     

基于DCT-CS的脉冲星周期超快速估计方法

受快速傅里叶变换（fast Fourier transform，FFT）的影响，基于FFT和压缩感知（compressive sensing, CS）的脉冲星周期快速估计算法的计算量大。为进一步减小计算量并提高计算精度，利用离散余弦变换（discrete cosine transform，DCT)取代FFT，提出了一种基于DCT-CS的脉冲星周期超快速估计算法。在该方法中，利用DCT提取脉冲星信号的低频部分构建低频DCT矩阵；构建畸变轮廓字典并获取累积轮廓；提出了利用最大值超分辨率稀疏恢复估计脉冲星周期的方法。仿真结果表明，DCT-CS的脉冲星周期估计精度达到了3.82×10-12 s，计算时间达到了9.31 ms。与FFT-CS相比，周期估计精度提高了约16%，计算时间缩短了约37.5%，实现了实时高精度的脉冲星周期估计。     

Fast approximate maximum likelihood period estimation from incomplete timing data

 To estimate the period of a periodic point process from noisy and incomplete observations, the classical periodogram algorithm is modified. The original periodogram algorithm yields an estimate by performing grid search of the peak of a spectrum, which is equivalent to the periodogram of the periodic point process, thus its performance is found to be sensitive to the chosen grid spacing. This paper derives a novel grid spacing formula, after finding a lower bound of the width of the spectral mainlobe. By employing this formula, the proposed new estimator can determine an appropriate grid spacing adaptively, and is able to yield approximate maximum likelihood estimate (MLE) with a computational complexity of O(n²). Experimental results prove that the proposed estimator can achieve better trade-off between statistical accuracy and complexity, as compared to existing methods. Simulations also show that the derived grid spacing formula is also applicable to other estimators that operate similarly by grid search.     

载波相位平滑伪距在GPS/SINS紧组合导航系统中的应用

针对GPS/SINS紧组合导航系统中伪距噪声大从而引起组合导航系统精度低的问题,提出了将载波相位平滑伪距方法引入到组合导航系统中,利用具有较高精度的载波相位对低精度伪距进行平滑滤波,在建立观测方程的同时,为了减小Kalman滤波器计算量,选用最佳4颗卫星的伪距与伪距率作为观测值,并提出了一种四面体选星法,该方法不需要进行矩阵求逆运算,减小了运算量。利用光纤陀螺捷联惯导系统与GPS接收机搭建了实验验证系统,通过车载实验对所提出的方法进行了验证,实验结果表明,经过载波相位平滑后的伪距噪声得到了降低,从而能够进一步提高GPS/SINS组合导航系统的定位精度,其位置误差均方根值相比无载波相位平滑减小近40%。     

CFAR detection and estimation for STAP radar

The algorithm presented here provides both a constant false-alarm rate (CFAR) detection and a maximum likelihood (ML) Doppler-bearing estimator of a target in a background of unknown Gaussian noise. A target is detected, and its parameters estimated within each range gate by evaluating a statistical test for each Doppler-angle cell and by selecting the cell with maximum output and finally comparing it with a threshold. Its CFAR performance is analyzed by the use of the sample matrix inversion (SMI) method and is evaluated in the cases of a fully adaptive space-time adaptive processing (STAP) and two partially adaptive STAPs. The performances of these criteria show that the probability of detection is a function only of the sample size K used to estimate the covariance matrix and a generalized signal-to-noise ratio. The choice of the number K is a tradeoff between performance and computational complexity. The performance curves demonstrate that the finer the resolution is, the poorer the detection capability. That means that one can trade off the accuracy of ML estimation with the performance of the CFAR detection criterion     

宽带信号近似最大似然方位估计快速算法

 针对短采样宽带信号近似最大似然方位估计（AML）计算量大的问题,将马尔可夫蒙特卡罗(MCMC)方法与近似最大似然方位估计相结合,提出一种基于完美抽样的近似最大似然方位估计快速算法（PAML）。该算法将AML算法的空间谱函数作为信号的概率分布函数,并利用完美抽样方法从该概率分布函数中抽样。与AML和遗传算法的对比实验研究表明,两目标情况下PAML算法在中低信噪比条件下的估计性能与AML和遗传算法性能相当,而计算量分别是二者的1/24和1/3。随着目标个数的增加,PAML算法的计算量优势将更加明显。     

Range, radial velocity, and acceleration MLE using radar LFM pulsetrain

An efficient implementation of the maximum likelihood estimator (MLE) is presented for the estimation of target range, radial velocity, and acceleration when the radar waveform consists of a wideband linear frequency modulated (LFM) pulse train. Analytic properties of the associated wideband ambiguity function are derived; in particular the ambiguity function, with acceleration set to zero, is derived in closed form. Convexity and symmetry properties of the ambiguity function over range, velocity, and acceleration are presented; these are useful for determining region and speed of convergence for recursive algorithms used to compute the MLE. In addition, the Cramer-Rao bound (CRB) is computed in closed form which shows that the velocity bound is decoupled from the corresponding bounds in range and acceleration. A fast MLE is then proposed which uses the Hough transform (HT) to initialize the MLE algorithm. Monte Carlo simulations show that the MLE attains the CRB for low to moderate signal-to-noise depending on the a priori estimates of range, velocity, and acceleration     

修正容积卡尔曼滤波数据域直接定位算法

为解决数据域直接定位（DPD）算法面临的计算压力,提高算法效率,提出一种基于修正的容积卡尔曼滤波（MCKF）的DPD算法。首先,融合各观测信号波达方向信息,利用子空间数据融合方法建立一种基于间接观测量的DPD滤波模型;然后,根据模型特点设计MCKF算法进行求解,解决间接观测量带来的噪声累积问题;最后,对算法计算量进行分析和对比,说明计算效率的提升。仿真结果表明,所提算法与基于最大似然遍历搜索和遗传算法的DPD算法相比,在相同的估计性能下,计算量下降明显,时效性显著提升,增加了算法实用价值。     

Doppler ambiguity resolution using multiple PRF

An algorithm for velocity ambiguity resolution in coherent pulsed Doppler radar using multiple pulse repetition frequencies (PRF) is presented. It relies on the choice of particular values for the PRFs. The folded frequency of the target signal is obtained by averaging the folded frequency estimates for each PRF, and a quasi maximum likelihood criterion is maximized for ambiguity order estimation. The fast implementation of this nonambiguous estimation procedure is based on the fast Fourier transform (FFT), The proposed waveform allows full exploitation of any (even) number of PRFs, which appears to be important for estimation improvement. The effects of the waveform parameters and the folded frequency estimation variance on the performance of the ambiguity order estimation procedure are evaluated theoretically and through computer simulations. Mean square error (MSE) curves are given to assess the Doppler frequency estimation accuracy. Finally, the new method is compared with a classical technique and the implementation of the algorithm in a clutter environment is addressed.     

基于MCIS技术的相位差变化率单站无源定位

针对最大似然估计（ML）方法求解测相位差变化率单站无源定位问题计算量大、定位慢的问题,本文提出一种利用蒙特卡洛重要性抽样技术（MCIS）高精度、低复杂度的估计方法。根据Pincus定理推导出ML问题的近似全局解,利用重要性抽样（IS）技术构建符合高斯分布概率密度（PDF）的重要性函数,作为样本选取的依据,通过逆变换采样获得样本集,统计样本均值直接得到辐射源位置估计结果。MCIS方法简单易实现且运算量低,能够克服传统ML估计多维网格搜索耗时较长的缺陷,而且对目标位置初始估计误差有较低的敏感性。实验结果表明,MCIS算法在相同测量噪声水平下,定位精度优于EKF、NLS算法,有效减小了初始化估计误差对算法定位精度的影响,也进一步讨论分析了算法参数和不同观测条件对定位性能的影响。     

基于鲁棒扩展卡尔曼滤波的双目视觉测量相对导航技术研究

针对非合作目标相对导航问题,为提高相对导航的精度和可靠性,采用双目视觉测量方法实现追踪器与目标器之间相对状态的测量,并基于鲁棒扩展卡尔曼滤波方法,提出了一种鲁棒相对导航滤波方法。仿真结果表明,该方法对相对导航系统模型中的不确定性具有良好的鲁棒性,且滤波精度较高。     

基于协方差矩阵高阶幂的二维DOA估计新算法

 针对稀疏分解方法进行均匀圆阵(UCA)的二维波达方向(DOA)估计运算复杂度大的问题,提出了一种基于协方差矩阵高阶幂稀疏分解的二维DOA估计新算法。该算法首先利用协方差矩阵高阶幂无需进行特征值分解和信源数估计的特性,构建了协方差矩阵高阶幂的稀疏分解向量;然后运用粒度分层思想,构造了粗区域估计和细方位估计的分层多粒度的快速分解模型,分层字典的长度大大减少,在保持估计精度的前提下,算法运算时间远小于现有的恒定冗余字典的稀疏分解方法,从而解决了基于稀疏分解的圆阵二维DOA估计问题。论文提出的算法与二维MUSIC算法相比,估计精度高,且能满足对相干信号的估计。仿真结果验证了算法的有效性和可行性。