Precision large scale air traffic surveillance using IMM/assignmentestimators

We present the development and implementation of a multisensor-multitarget tracking algorithm for large scale air traffic surveillance based on interacting multiple model (IMM) state estimation combined with a 2-dimensional assignment for data association. The algorithm can be used to track a large number of targets from measurements obtained with a large number of radars. The use of the algorithm is illustrated on measurements obtained from 5 FAA radars, which are asynchronous, heterogeneous, and geographically distributed over a large area. Both secondary radar data (beacon returns from cooperative targets) as well as primary radar data (skin returns from noncooperative targets) are used. The target IDs from the beacon returns are not used in the data association. The surveillance region includes about 800 targets that exhibit different types of motion. The performance of an IMM estimator with linear motion models is compared with that of the Kalman filter (KF). A number of performance measures that can be used on real data without knowledge of the ground truth are presented for this purpose. It is shown that the IMM estimator performs better than the KF. The advantage of fusing multisensor data is quantified. It is also shown that the computational requirements in the multisensor case are lower than in single sensor case, Finally, an IMM estimator with a nonlinear motion model (coordinated turn) is shown to further improve the performance during the maneuvering periods over the IMM with linear models     

Two-step optimal estimator for three dimensional target tracking

This study presents an adaptation of a novel estimation methodology to the general nonlinear three-dimensional problem of tracking a maneuvering target. The two-step optimal estimator (TSE) suggests an attractive alternative to the standard extended Kalman filter (EKF). A superior performance is accomplished by dividing the estimation problem into two steps: a linear first step and a nonlinear second step. The target tracking performance of the TSE is shown to be better than an EKF implemented in either inertial or modified spherical coordinates. In the passive case, where bearing/elevation angles only are measured, the TSE yields excellent range and target acceleration estimates. In the active case, where range measurement is available as well, a homing missile employing closed-loop optimal guidance based on the TSE state estimates obtains smaller miss distances than with either versions of the EKF.     

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.     

A novel multistage estimation of signal parameters

A multistage estimation scheme is presented for estimating the parameters of a received carrier signal possibly phase-modulated by unknown data and experiencing very high Doppler, Doppler rate, etc. Such a situation arises, for example, in the case of the Global Positioning Systems (GPS). In the proposed scheme, the first-stage estimator operates as a coarse estimator of the frequency and its derivatives, resulting in higher RMS estimation errors but with a relatively small probability of the frequency estimation error exceeding one-half of the sampling frequency (an event termed cycle slip). The second stage of the estimator operates on the error signal available from the first stage, refining the overall estimates, and in the process also reduces the number of cycle slips. The first-stage algorithm is a modified least-squares algorithm operating on the differential signal model and referred to as differential least squares (DLS). The second-stage algorithm is an extended Kalman filter, which yields the estimate of the phase as well as refining the frequency estimate. A major advantage of the proposed algorithm is a reduction in the threshold for the received carrier power-to-noise power spectral density ratio (CNR) as compared with the threshold achievable by either of the algorithms alone     

EM-ML algorithm for track initialization using possibly noninformative data

Initializing and maintaining a track for a low observable (LO) (low SNR, low target detection probability and high false alarm rate) target can be very challenging because of the low information content of measurements. In addition, in some scenarios, target-originated measurements might not be present in many consecutive scans because of mispointing, target maneuvers, or erroneous preprocessing. That is, one might have a set of noninformative scans that could result in poor track initialization and maintenance. In this paper an algorithm based on the expectation-maximization (EM) algorithm combined with maximum likelihood (ML) estimation is presented for tracking slowly maneuvering targets in heavy clutter and possibly noninformative scans. The adaptive sliding-window EM-ML approach, which operates in batch mode, tries to reject or weight down noninformative scans using the Q-function in the M-step of the EM algorithm. It is shown that target features in the form of, for example, amplitude information (AI), can also be used to improve the estimates. In addition, performance bounds based on the supplemented EM (SEM) technique are also presented. The effectiveness of new algorithm is first demonstrated on a 78-frame long wave infrared (LWIR) data sequence consisting of an Fl Mirage fighter jet in heavy clutter. Previously, this scenario has been used as a benchmark for evaluating the performance of other track initialization algorithms. The new EM-ML estimator confirms the track by frame 20 while the ML-PDA (maximum likelihood estimator combined with probabilistic data association) algorithm, the IMM-MHT (interacting multiple model estimator combined with multiple hypothesis tracking) and the EVIM-PDA estimator previously required 28, 38, and 39 frames, respectively. The benefits of the new algorithm in terms of accuracy, early detection, and computational load are illustrated using simulated scenarios as well.     

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 Velocity Estimation Algorithm of Moving Targets using Single Antenna SAR

A new algorithm is proposed for velocity estimation of moving targets in single antenna synthetic aperture radar (SAR). Based on the fact that different velocity vectors cause different geometrical figures of the two-dimensional (2-D) signature in the range-Doppler (RD) domain, this algorithm estimates the azimuth and range velocities by a 2-D search such that the range cell migration correction (RCMC) and the second range compression (SRC) are correctly performed. It is shown that, using the proposed algorithm, the Doppler ambiguity problem can be avoided and satisfactory accurate velocity estimation can be obtained in high signal-to-clutter ratio (SCR) scenarios.     

一种新颖的超高分辨率SAR频率域反投影成像算法

超高分辨率条件下,机载合成孔径雷达（SAR）发射信号带宽大,合成孔径时间比较长,对成像处理算法的精度和效率要求较高。现有近似频率域处理和时间域滤波反投影（FBP）算法聚焦SAR数据时均存在诸多问题。基于微局部分析方法,提出了一种新颖的频率域滤波反投影（FD-FBP）成像处理方案。首先,利用Keystone变换简化了数据距离多普勒（RD）域徙动表达式。然后,在RD域进行反投影操作,对参考位置处反投影数据进行移位、相位补偿和FFT等操作即可以得到图像,从而在保证算法精确性的前提下有效降低了运算效率,实现了频率域方法的高效率和时间域方法的精确性特点的结合。最后,点目标仿真和实测数据处理以及与FBP等算法的对比验证了该方法的有效性。     

Angle of arrival estimation of coherent signals using an arraydoublet in motion

Shown here is how the estimation of signal parameters via relational invariance techniques (ESPRIT) algorithm may be used with a single pair of antennas in motion to estimate angles of arrival (AOA) for coherent signals. The approach exploits the Doppler frequency shifts caused by the doubler in motion. With this estimator, the number of signals that can be handled is not limited by the size of the array, as in the usual ESPRIT application, but by an adjustable parameter. A theoretical performance analysis of the estimator and typical examples showing the use of this estimator are given     

Low observable target motion analysis using amplitude information

In conventional passive and active sonar system, target amplitude information (AI) at the output of the signal processor is used only to declare detections and provide measurements. We show that the AI can be used in passive sonar system, with or without frequency measurements, in the estimation process itself to enhance the performance in the presence of clutter where the target-originated measurements cannot be identified with certainty, i.e., for “low observable” or “dim” (low signal-to-noise ratio (SNR)) targets. A probabilistic data association (PDA) based maximum likelihood (ML) estimator for target motion analysis (TMA) that uses amplitude information is derived. A track formation algorithm and the Cramer-Rao lower bound (CRLB) in the presence of false measurements, which is met by the estimator even under low SNR conditions, are also given. The CRLB is met by the proposed estimator even at 6 dB in a cell (which corresponds to 0 dB for 1 Hz bandwidth in the case of a 0.25 Hz frequency cell) whereas the estimator without AI works only down to 9 dB. Results demonstrate improved accuracy and superior global convergence when compared with the estimator without AI. The same methodology can be used for bistatic radar     

Satellite lithium-ion battery remaining useful life estimation with an iterative updated RVM fused with the KF algorithm

Lithium-ion batteries have become the third-generation space batteries and are widely utilized in a series of spacecraft. Remaining Useful Life (RUL) estimation is essential to a spacecraft as the battery is a critical part and determines the lifetime and reliability. The Relevance Vector Machine (RVM) is a data-driven algorithm used to estimate a battery’s RUL due to its sparse feature and uncertainty management capability. Especially, some of the regressive cases indicate that the RVM can obtain a better short-term prediction performance rather than long-term prediction. As a nonlinear kernel learning algorithm, the coefficient matrix and relevance vectors are fixed once the RVM training is conducted. Moreover, the RVM can be simply influenced by the noise with the training data. Thus, this work proposes an iterative updated approach to improve the long-term prediction performance for a battery’s RUL prediction. Firstly, when a new estimator is output by the RVM, the Kalman filter is applied to optimize this estimator with a physical degradation model. Then, this optimized estimator is added into the training set as an on-line sample, the RVM model is re-trained, and the coefficient matrix and relevance vectors can be dynamically adjusted to make next iterative prediction. Experimental results with a commercial battery test data set and a satellite battery data set both indicate that the proposed method can achieve a better performance for RUL estimation.     

基于先验误差模型的机载高分宽幅DBF-SAR自聚焦算法

高分宽幅（HRWS）数字波束形成（DBF）合成孔径雷达（SAR）利用多通道空间采样代替部分时域采样,可以有效缓解SAR成像时高分辨率与宽测绘带间的矛盾,具有重要的军用和民用价值。现有常规DBF-SAR成像算法都假设雷达传感器相对位置精确已知,实际应用中受传感器位置测量误差影响,由位置不精确导致的相位误差会严重影响DBF-SAR高精度成像能力。在极坐标格式算法（PFA）框架下,推导了DBF-SAR成像处理后,残留相位误差的解析模型,分析了该误差对成像质量的影响。依据推导的先验相位误差解析结构模型,提出了一种基于图像对比度最优化准则的自聚焦算法。新算法通过引入先验相位结构信息,极大降低了待估参数的空间维数,可以同时改善自聚焦算法的参数估计精度和计算效率。数据处理结果验证了理论分析的正确性和所提算法的有效性。     

基于改进容积卡尔曼滤波的奇异避免姿态估计

 利用矢量进行卫星姿态估计可以归结为非线性滤波问题。为了提高卫星姿态估计的精度,利用龙贝格-马尔塔(LM)迭代算法改进了容积卡尔曼滤波(CKF)。继而,提出改进容积卡尔曼滤波与四元数结合的容积四元数估计器(CQE),有效地避免了卫星大角度机动出现的奇异现象。进一步,给出了一种与影子修正罗德里格参数切换的容积修正罗德里格参数估计器(CME)。仿真对比表明,初始误差较大时容积修正罗德里格参数估计器具有更好的收敛速度和鲁棒性。     

Super resolution feature extraction of moving targets

High range resolution (HRR) moving target indicator (MTI) is becoming increasingly important for many military and civilian applications such as the detection and classification of moving targets in strong clutter background. We consider the problem of extracting the HRR features of moving targets with very closely spaced scatterers in the presence of strong stationary clutter, where the range migration and Doppler frequency are taken into account. A relaxation-based algorithm, which is robust and computationally simple, is proposed to deal with the above problem. Numerical results have shown that the proposed algorithm exhibits super resolution and excellent estimation performance     

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

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

基于自适应差分进化算法的变循环发动机模型求解方法研究

为了降低传统迭代算法在求解变循环发动机非线性模型时对初值的依赖性,将模型的求解问题转换为求最小值的优化问题,引入差分进化算法进行模型的求解,并提出一种自适应差分进化算法(ADE)。ADE借助轮盘赌选择法,利用种群的进化经验可以自适应的选择最适合当前种群的差分策略与算法控制参数。针对变循环发动机四个典型工作点的模型求解问题,研究了标准差分进化算法(SDE)的控制参数对其性能的影响,获取了SDE在求解四个典型工作点时的最优控制参数组合,对比分析了ADE与SDE的性能差异,最后研究了种群规模对ADE性能的影响。结果表明：SDE在求解发动机模型时具有较好的鲁棒性,在求解不同工作点时算法的最优控制参数并不完全相同;相比于使用最优控制参数的SDE,ADE可以在不影响算法鲁棒性的情况下提升效率50%以上;减少ADE的种群规模会在提升算法效率的同时破坏鲁棒性。     

Optimal sensor selection strategy for discrete-time stateestimators

A novel sensor selection strategy is introduced, which can be implemented on-line in time-varying discrete-time system. We consider a case in which several measurement subsystem are available, each of which may be used to drive a state estimation algorithm. However, due to practical implementation constraints (such as the ability of the on-board computer to process the acquired data), only one of these subsystems can actually by utilized at a measurement update. An algorithm is needed, by which the optimal measurement subsystem to be used is selected at each sensor selection epoch. The approach described is based on using the square root V-Lambda information filter as the underlying state estimation algorithm. This algorithm continuously provides its user with the spectral factors of the estimation error covariance matrix, which are used in this work as the basis for an on-line decision procedure by which the optimal measurement strategy is derived. At each sensor selection epoch, a measurement subsystem is selected, which contributes the largest amount of information along the principal state space direction associated with the largest current estimation error. A numerical example is presented, which demonstrates the performance of the new algorithm. The state estimation problem is solved for a third-order time-varying system equipped with three measurement subsystem, only one of which can be used at a measurement update. It is shown that the optimal measurement strategy algorithm enhances the estimator by substantially reducing the maximal estimation error     

High-resolution SAR imaging with angular diversity

We propose to use the APES (amplitude and phase estimation) approach for the spectral estimation of gapped data and synthetic aperture radar (SAR) imaging with angular diversity. A relaxation-based algorithm, referred to as GAPES (Gapped-data APES), is proposed, which includes estimating the spectrum via APES and filling in the gaps via a least squares (LS) fitting. For SAR imaging with angular diversity data fusion, we perform one-dimensional (1-D) windowed fast Fourier transforms (FFTs) in range, use the GAPES algorithm to interpolate the gaps in the aperture for each range, apply 1-D inverse FFTs (IFFTs) and dewindow in range, and finally apply the two-dimensional (2-D) APES algorithm to the interpolated matrix to obtain the 2-D SAR image. Numerical results are presented to demonstrate the effectiveness of the proposed algorithm     

一种改进的PGA最大似然相位估计法

 讨论了合成孔径雷达( Synthet ic Aper ture Radar ) SAR) 成像处理过程中一种不需要直接求解矩阵特征值及其特征向量, 但仍然基于最大似然估计的应用于相位梯度自聚焦算法( Phase Gr adient Algorithm)PGA) 中的相位估计方法, 且主要从计算量及数据处理精度等方面与传统的特征向量方法相比较。结果表明,这种基于迭代方法计算特征向量相位的算法在SAR 成像自聚焦过程中具有实用意义。     

On joint track initiation and parameter estimation under measurement origin uncertainty

The problem of joint detection and estimation for track initiation under measurement origin uncertainty is studied. The two well-known approaches, namely the maximum likelihood estimator with probabilistic data association (ML-PDA) and the multiple hypotheses tracking (MHT) via multiframe assignment, are characterized as special cases of the generalized likelihood ratio test (GLRT) and their performance limits indicated. A new detection scheme based on the optimal gating is proposed and the associated parameter estimation scheme modified within the ML-PDA framework. A simplified example shows the effectiveness of the new algorithm in detection performance under heavy clutter. Extension of the results to state estimation with measurement origin uncertainty is also discussed with emphasis on joint detection and recursive state estimation.