基于最大似然估计的TDOA/FDOA无源定位偏差补偿算法

在到达时差/到达频差(TDOA/FDOA)无源定位系统中,定位问题的非线性使得定位的结果存在偏差,特别是在噪声较大或者接收站布站不合理的情况下,定位的偏差尤其显著。针对这一问题,提出了一种基于最大似然估计的偏差补偿算法。该方法分为3步:首先,利用最大似然估计器对目标的位置和速度进行求解;其次,通过利用目标定位的估计值和含噪的测量值,对目标的位置和速度偏差值进行理论分析和推导;最后,将最大似然估计解减去理论偏差值,得到经过偏差补偿的新的目标定位解。理论分析和实验仿真证明,在一定噪声的情况下,所推导的目标位置和速度的理论偏差值与实际偏差值相符,并且经过偏差补偿后的定位算法,在保持目标定位的均方根误差(RMSE)与原最大似然算法一致的情况下,目标的位置和速度偏差值远远小于原最大似然算法的偏差值,目标定位精度得到了有效的提高。     

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

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

Emitter localization using clustering-based bearing association

A closed-form emitter location estimator using time difference of arrival (TDOA) measurements is developed based on triangulation of hyperbolic asymptotes. The problem of associating the asymptotes with the emitter is solved by clustering the bearing angles of the linear asymptotes using a kernel density estimate. A closed-form estimate of the emitter location is obtained from triangulation of the clustered bearings using a weighted version of the pseudolinear estimator. By way of simulation examples, the proposed closed-form estimator is shown to outperform the computationally demanding and divergence-prone maximum likelihood (ML) estimator at moderate TDOA noise levels.     

一种新的基于TDOA与GROA的无源定位算法

针对基于到达时间差(TDOA)与到达增益比(GROA)的辐射源无源定位问题,提出了一种新的定位算法。首先通过引入一个中间变量,根据TDOA和GROA测量模型,构造一个约束加权最小二乘(CWLS)估计。由于这个CWLS问题是非凸的优化问题,现有的方法不能很好地求解。为此,提出了一种有效的方法可以求解到其全局最优解。最后,所提算法被推广到观测站存在自定位误差时的定位求解。计算机仿真结果验证了所提算法能够获得优于传统两步加权最小二乘法(2WLS)的定位性能,能够在更大的噪声条件下达到克拉美罗下界(CRLB)。     

DOA estimation by exploiting the amplitude modulation induced by antenna scanning

We propose a beamsplitting-like approach to estimate the directions of arrival (DOA) of multiple radar targets present in the mainlobe of a rotating antenna. The proposed method is based on the maximum likelihood (ML) technique and it avoids the need for a difference channel by exploiting knowledge of the antenna main beam pattern. Two scenarios are considered: multiple targets with unknown deterministic complex amplitudes and multiple targets with Gaussian distributed random complex amplitudes. The performance of the proposed estimator is investigated through Monte Carlo simulation and it is compared with the Cramer-Rao lower bound (CRLB).     

Estimation of singularities for intercept point forecasting

In this paper, a recursive approach (an algorithm for estimation of singularity (AES)) is proposed to forecast the intercept point of a target and the pursuing interceptor recognized as the estimated singularity of a nodal cubic curve fitted to the data. The data comprises direction of arrival (DOA) estimates of both target and interceptor obtained at regular intervals of time using the maximum likelihood (ML) DOA estimation method. The estimates of coefficients of the cubic polynomial are given by a recursive least squares solution. From these coefficients, closed-form solutions for angle of interception and intercept time are obtained which are the forecasted coordinates of the intercept point. Experimental results demonstrate performance of the proposed algorithm     

传感器位置误差条件下的约束总体最小二乘时差定位算法

 现代定位系统中,传感器往往被安放在运动平台上,其位置无法精确得知,存在估计误差,将严重影响对目标的定位精度。针对这一问题,提出基于约束总体最小二乘(CTLS)的到达时差(TDOA)定位算法。首先通过引入中间变量,将非线性TDOA定位方程转化为伪线性方程,再利用CTLS技术,全面考虑伪线性方程所有系数中的噪声。在此基础上推导了定位方程的目标函数,再根据牛顿迭代方法,进行数值迭代,快速得到精确解。采用一阶小噪声扰动分析方法,对该算法的理论性能进行了分析,证明了算法的无偏性和逼近克拉美-罗下限(CRLB)。仿真实验表明,该算法克服了现有总体最小二乘(TLS)算法不能达到CRLB、两步加权最小二乘(two-step WLS)算法在较高噪声时性能发散的缺陷,在较高噪声时定位精度仍然能达到CRLB。     

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.     

Radar power multiplier for acquisition of low observables using anESA radar

In this paper the acquisition of a low observable (LO) incoming tactical ballistic missile using the measurements from a surface based electronically scanned array (ESA) radar is presented. We present a batch maximum likelihood (ML) estimator to acquire the missile while it is exo-atmospheric. The proposed estimator, which combines ML estimation with the probabilistic data association (PDA) approach resulting in the ML-PDA algorithm to handle false alarms, also uses target features. The use of features facilitates target acquisition under low signal-to-noise ratio (SNR) conditions. Typically, ESA radars operate at 13-20 dB, whereas the new estimator is shown to be effective even at 4 dB SNR (in a resolution cell, at the end of the signal processing chain) for a Swerling III fluctuating target, which represents a significant counter-stealth capability. That is, this algorithm acts as an effective “power multiplier” for the radar by about an order of magnitude. An approximate Cramer-Rao lower bound (CRLB), quantifying the attainable estimation accuracies and shown to be met by the proposed estimator, is derived as well     

A direct position determination method with combined TDOA and FDOA based on particle filter

The localization of a stationary transmitter using moving receivers is considered. The original Direct Position Determination (DPD) methods, with combined Time Difference of Arrival (TDOA) and Frequency Difference of Arrival (FDOA), do not perform well under low Signal-to-Noise Ratio (SNR), and worse still, the computation cost is difficult to accept when the computational capabilities are limited. To get better positioning performance, we present a new DPD algorithm that proves to be more computationally efficient and more precise for weak signals than the conventional approach. The algorithm partitions the signal received with the same receiver into multiple non-overlapping short-time signal segments, and then uses the TDOA, the FDOA and the coherency among the short-time signals to locate the target. The fast maximum likelihood estimation, one iterative method based on particle filter, is designed to solve the problem of high computation load. A secondary but important result is a derivation of closed-form expressions of the Cramer-Rao Lower Bound (CRLB). The simulation results show that the algorithm proposed in this paper outperforms the traditional DPD algorithms with more accurate results and higher computational efficiency, and especially at low SNR, it is more close to the CRLB.     

一种针对恒模信号的运动单站直接定位算法

相比于常规的"测向+位置估计"两步定位模式,以Weiss等提出的目标直接位置确定(DPD)算法具有估计精度高、分辨能力强和无需数据关联等诸多优点。基于该类定位算法的基本理念,提出了一种利用单个运动天线阵列对恒模(即相位调制)信号的DPD算法。首先,依据最大似然(ML)准则以及恒模信号的恒包络特征,建立了相应的直接定位优化模型;接着,根据优化函数的代数特征提出了一种有效的多参量交替迭代算法,用以获得ML估计器的最优数值解;此外,推导了针对恒模信源的位置直接估计方差的克拉美罗界(CRB),从而为新算法的定位精度提供定量的理论下界。仿真实验表明:相比于已有的基于单个运动天线阵列的直接定位算法以及传统的两步定位算法,通过利用恒模信号的恒包络特征可以明显提高目标直接定位的估计精度。     

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.     

WGS-84模型下时差频差半定规划定位算法

利用多颗卫星的时差频差对辐射源进行位置和速度的测量,其本质意义上是一个含有噪声项的高度非线性方程组求解问题,针对地面目标而言,可以采用基于WGS-84地球模型作为目标位置和速度约束,更进一步的增加了定位系统的复杂性。提出了一种基于半定规划(SDP)的定位解算法,将非线性方程求解问题通过适当的松弛,转化为半定优化(SDO)的问题,借助于业界较为成熟的CVX等优化软件进行定位求解,并研究了该模型条件下的克拉美罗下界(CRLB)。仿真结果表明,该算法能够较好地逼近克拉美罗下界。     

非线性最优解在气动参数辨识中的应用

由于采用最大似然算法进行了扭动参数辨识通常存在着局极值问题，为了寻求非线性全局最优解，提高估计精度，介绍了以混沌优化算法作为搜索工具。与最大似然处算法相结合，求解非凸优化问题全局最优解的方法，算例表明，该方法是有效可行的，并且估计精度也较高。     

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

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

基于定位误差修正的运动目标TDOA/FDOA无源定位方法

针对时差(TDOA)、频差(FDOA)无源定位的两步加权最小二乘(TSWLS)方法定位均方根误差(RMSE)和定位偏差适应测量噪声能力差的问题,在分析了影响两步法定位性能的因素基础上提出一种基于一阶泰勒级数展开的定位误差修正方法。该方法的第1步和两步法相同;其第2步避免了两步法第2步中引入估计偏差的平方运算,利用一阶泰勒级数展开得到第1步定位误差的线性最小均方估计,修正第1步定位结果得到目标位置和速度的最终估计,从理论上证明了该方法可以达到定位的克拉美罗下限(CRLB)。计算机仿真对比了新方法和TSWLS方法、基于泰勒级数(TS)展开的迭代极大似然(ML)方法以及约束总体最小二乘(CTLS)方法的定位性能,新算法复杂度和两步法相当,且均方误差和定位偏差低于两步法、泰勒级数法和CTLS方法。     

A fast pulse phase estimation method for X-ray pulsar signals based on epoch folding

X-ray pulsar-based navigation (XPNAV) is an attractive method for autonomous deep-space navigation in the future. The pulse phase estimation is a key task in XPNAV and its accuracy directly determines the navigation accuracy. State-of-the-art pulse phase estimation techniques either suffer from poor estimation accuracy, or involve the maximization of generally non-convex object function, thus resulting in a large computational cost. In this paper, a fast pulse phase estimation method based on epoch folding is presented. The statistical properties of the observed profile obtained through epoch folding are developed. Based on this, we recognize the joint prob-ability distribution of the observed profile as the likelihood function and utilize a fast Fourier transform-based procedure to estimate the pulse phase. Computational complexity of the proposed estimator is analyzed as well. Experimental results show that the proposed estimator significantly outperforms the currently used cross-correlation (CC) and nonlinear least squares (NLS) estima-tors, while significantly reduces the computational complexity compared with NLS and maximum likelihood (ML) estimators.     

传感器参数误差下的运动目标TDOA/FDOA无源定位算法

在运动目标无源定位系统中，许多算法的前提是精确已知传感器的位置以及速度，但实际情况下可利用的传感器的参数均会存在一些噪声扰动。针对这一问题，提出一种改进的两步加权最小二乘（TSWLS）时差（TDOA）与频差（FDOA）定位算法。该算法是一种闭式算法并且分为2步。第1步与经典的两步加权最小二乘算法相同，第2步进一步研究了额外变量与目标参数之间的关系并且建立了新的矩阵方程。随后，利用加权最小二乘技术给出了最终解。理论分析证明了在测量噪声较小时该算法能够达到克拉美罗界（CRLB）。所提算法具有计算复杂度低，实时性高的优点；另外，经过适当的维度调整，该算法同样适用于对多非相交源进行定位求解。计算机仿真进一步证明了理论分析的正确性。     

Moving scanning emitter tracking by a single observer using time of interception: Observability analysis and algorithm

The target motion analysis (TMA) for a moving scanning emitter with known fixed scan rate by a single observer using the time of interception (TOI) measurements only is investigated in this paper.By transforming the TOI of multiple scan cycles into the direction difference of arrival (DDOA) model,the observability analysis for the TMA problem is performed.Some necessary conditions for uniquely identifying the scanning emitter trajectory are obtained.This paper also proposes a weighted instrumental variable (WIV) estimator for the scanning emitter TMA,which does not require any initial solution guess and is closed-form and computationally attractive.More importantly,simulations show that the proposed algorithm can provide estimation mean square error close to the Cramer-Rao lower bound (CRLB) at moderate noise levels with significantly lower estimation bias than the conventional pseudo-linear least square (PLS) estimator.     

无需辅助数据的分布式目标自适应检测器

在非高斯背景和没有辅助数据的条件下,研究了高分辨率雷达分布式目标的自适应检测问题.首先采用有序检测理论和协方差矩阵的迭代估计方法粗略估计散射点集合,进一步利用迭代估计方法获得协方差矩阵的近似最大似然估计,提出了无需辅助数据的自适应检测器(ADWSD).ADWSD在非高斯背景下具有近似恒虚警率特性,且检测性能远好于修正的...