Passive localization by a single moving observer using TOA only with unknown SRI: Observability analysis and algorithm

Passive localization by a single moving observer using Time of Arrival (TOA) only with an unknown Signal Repetition Interval (SRI) is investigated in this paper. Observability analysis is performed first. The observability condition for uniquely determining the emitter position and SRI is derived. The conditional Cramer-Rao Lower Bound (CRLB) is also analyzed. It is found that the ambiguity of the SRI integer of the first TOA does not affect the theoretical estimation precision of the emitter position and SRI. A Reference-Fixed Differential TOA (RFDTOA)-based Iterative Maximum Likelihood Estimator (IMLE) is proposed, which only needs O(M) computational operations. Theoretical analysis and simulation results show that the Mean Square Error (MSE) of the proposed algorithm could attain the CRLB with moderate Gaussian measurement noise.     

Performance of a Square Law Pseudonoise Ranging Time-of-Arrival Estimator

The performance of a square law time-of-arrival (TOA) estimator that has been proposed for use in ASTRO-DABS, part of a possible satellite-based fourth generation air traffic control system is considered. The transmitted message consists of a pulse amplitude modulated (PAM) ranging sequence that, due to transmitter characteristics, is corrupted by an unknown frequency offset. The optimum TOA estimator, for the case of no frequency uncertainty, is first presented, together with a lower bound on the variance of the estimate generated. This is followed by the consideration of a suboptimum TOA estimator for which a high signal-to-noise ratio (SNR) performance analysis is carried out; here, the effects of frequency uncertainty are included. Next, the zero-crossing properties of the derivative of the (suboptimum) estimation statistic are presented and the results used to derive an upper bound to the TOA estimate variance that is valid for all SNR values. This latter result is significant because it displays the system threshold effect and complements performance lower bounds that may be derived via other methods. In addition, the method presented here may be applied to other optimum and suboptimum systems where a discrete set of parameters is to be estimated.     

Hard-constrained versus soft-constrained parameter estimation

The paper aims at contrasting two different ways of incorporating a priori information in parameter estimation, i.e., hard-constrained and soft-constrained estimation. Hard-constrained estimation can be interpreted, in the Bayesian framework, as maximum a posteriori probability (MAP) estimation with uniform prior distribution over the constraining set, and amounts to a constrained least-squares (LS) optimization. Novel analytical results on the statistics of the hard-constrained estimator are presented for a linear regression model subject to lower and upper bounds on a single parameter. This analysis allows to quantify the mean squared error (MSE) reduction implied by constraints and to see how this depends on the size of the constraining set compared with the confidence regions of the unconstrained estimator. Contrastingly, soft-constrained estimation can be regarded as MAP estimation with Gaussian prior distribution and amounts to a less computationally demanding unconstrained LS optimization with a cost suitably modified by the mean and covariance of the Gaussian distribution. Results on the design of the prior covariance of the soft-constrained estimator for optimal MSE performance are also given. Finally, a practical case-study concerning a line fitting estimation problem is presented in order to validate the theoretical results derived in the paper as well as to compare the performance of the hard-constrained and soft-constrained approaches under different settings     

传感器位置误差条件下仅用到达频率差的无源定位性能分析

传感器自身的位置误差对辐射源的无源定位精度可能有较大影响,详细分析并推导了在存在传感器位置误差条件下,仅使用到达频率差(FDOA)参数对辐射源定位时,定位精度的克拉美罗下限(CRLB)和均方误差(MSE).在传感器位置误差条件下,提出仅用FDOA来同时估计辐射源和传感器位置的泰勒级数方法,并证明该定位方法理论上的MSE...     

一种频域插值变化迭代频率估计方法

针对现有频率估计算法存在的复杂度高、频率估计能力弱、估计结果均方差大等缺点,在固定迭代AM（Aboutanios—Mulgrew）无偏频率估计算法基础上,提出一种频域插值变化迭代频率估计算法,推导了不同迭代参数实现无偏估计的充分条件,证明了有偏估计时本算法的收敛性和偏离度,通过设置不同迭代参数,可以实现无偏或有偏估计。仿真分析表明：当具有较高信噪比时,在整个频率估计范围内,该方法均方误差接近CRLB（Cramer-RaoLowerBound,克拉美一罗下限）;当FFT（FastFourierTransform,快速傅里叶变换）粗估计残余频率接近0．5时,该方法的均方误差优于CRLB,为CRLB的96％。     

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.     

Optimal below Threshold Delay Estimation for Radio Signals

Time delay estimation of a random signal observed at twospatially separated sites is a problem of considerable practicalinterest. Radio-frequency (RF) signals pose a difficulty that does notarise while treating sonar or other low frequency signals: the signalprocessing can be performed only on low frequency versions of theobserved signals. We propose an estimation scheme, based oncomplex envelopes, the mean square error (MSE) of which achievesthe performance predicted by Barankin, Bellini-Tartara, and Weiss-Weinsteinlower bounds. The proposed scheme was tested byextensive computer simulations, and the results are in closeagreement with theory predictions. Moreover, a real-time hardwareimplementation of the algorithm was developed and testedsuccessfully using real RF signals.     

Finite data performance of MUSIC and minimum norm methods

In the direction of arrival (DOA) estimation problem, we encounter both finite data and insufficient knowledge of array characterization. It is therefore important to study how subspace-based methods perform in such conditions. We analyze the finite data performance of the multiple signal classification (MUSIC) and minimum norm (min. norm) methods in the presence of sensor gain and phase errors, and derive expressions for the mean square error (MSE) in the DOA estimates. These expressions are first derived assuming an arbitrary array and then simplified for the special case of an uniform linear array with isotropic sensors. When they are further simplified for the case of finite data only and sensor errors only, they reduce to the recent results given previously (1989, 1991). Computer simulations are used to verify the closeness between the predicted and simulated values of the MSE     

A novel pulse TOA analysis technique for radar identification

A new processing technique is presented for identifying individual pulse trains in a buffer load of pulse time-of-arrival (TOA) data from radar warning receiver sets currently in use for electronic warfare (EW). The method is based on the numerical characteristics of the matrix of the differences of the TOA data points. In particular, for a uniform pulse repetition interval (PRI) emitter, this matrix is a symmetrical harmonic Toeplitz matrix the inverse of which has banded tridiagonal form. This form leads to an estimator for the PRI in terms of the trace of the inverted matrix     

航空通信中基于PN码的信道估计算法研究

航空移动通信系统（空空信道）中,大的多径时延和多普勒频移会造成非常严重的符号间干扰,高速宽带下尤为突出,使得传统的频域信道估计算法不得不通过提高运算长度和复杂度来实现性能,而航空通信对实时性的要求使其难以工程实现。本文针对空空信道时延长而衰落小的条件,提出一种基于PN码的时域信道估计算法,连续发送N段m序列作为训练序列...     

Interlaced Kalman filtering of 3D angular motion based on Euler'snonlinear equations

A novel Kalman filtering technique is presented that reduces the mean-square-error (MSE) between three-dimensional (3D) actual angular velocity values and estimated ones by an order of magnitude (when compared with the MSE resulting from direct measurements) even under extremely low signal-to-noise ratio conditions. The filtering problem is nonlinear in nature because the dynamics of 3D angular motion are described by Euler's equations. This nonlinear set of differential equations state that the angular acceleration in one axis is proportional to the torque applied to that axis, and to the products of angular velocity components in the other two axes of rotation. Instead of using extended Kalman filtering techniques to solve this complex problem, the authors developed a new approach where the nonlinear Euler's model is decomposed into two pseudolinear models (primary and secondary). The first model describes the time progression of the state vector containing the linear terms, while the other characterizes the propagation of the state vector containing the nonlinearities. This makes it possible to run two interlaced discrete-linear Kalman filters simultaneously. One filter estimates the values of the state vector containing the linear terms, while the other estimates the values of the state vector containing the nonlinear terms in the system. These estimates are then recombined, solving the nonlinear estimation process without linearizing the system. Thus, the new approach takes advantage of the simplicity, computational efficiency and higher convergence speed of the linear Kalman filter form and it overcomes many of the drawbacks typical of conventional extended Kalman filtering techniques. The high performance and effectiveness of this method is demonstrated through a computer simulation case study     

Principles and Simulation of JTIDS Relative Navigation

The time-synchronous operation and high accuracy time-of arrival (TOA) measurement capability of Joint Tactical Information Distribution System (JTIDS) terminals makes possible a high performance relative nagivaton (RELNAV) function through addition of only software in the terminal's computer program. The principles of operation, the basic observation equations, and the system architecture for both absolute (geographic) and relative navigation are described. Sequential passive ranging by means of the TOA measurements, in conjunction with appropriate source selection logic and a recursive (e. g., Kalman) filter mechanization are employed to determine the user's position, velocity, and time bias. The filter algorithms and error sources, the software functional flow, and some simulation results are presented.     

大气湍流退化图像多帧盲复原

针对图像盲复原时PSF（点扩散函数）尺寸参数对复原效果影响较大,提出一种点扩散函数尺寸自适应估计算法,该算法以均方误差为代价函数,通过遍历搜索点扩散函数有限支持域,克服了点扩散函数尺寸参数选择难的问题。仿真结果表明,估计的点扩散函数尺寸接近真实尺寸,对于大气湍流变化快的情形,改进算法取得了较好的图像复原效果。     

Algebraic two-satellite TOA/FOA position solution on an ellipsoidal Earth

A direct algebraic solution is given for estimates of the ambiguous locations of an RF emitter on the surface of an ellipsoidal Earth given time-of-arrival (TOA) and frequency-of-arrival (FOA) measurements at two satellites.     

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.     

One-step solution for the multistep out-of-sequence-measurement problem in tracking

In multisensor target tracking systems measurements from the same target can arrive out of sequence. Such "out-of-sequence" measurement (OOSM) arrivals can occur even in the absence of scan/frame communication time delays. The resulting problem - how to update the current state estimate with an "older" measurement - is a nonstandard estimation problem. It was solved first (suboptimally, then optimally) for the case where the OOSM lies between the two last measurements, i.e, its lag is less than a sampling interval - the 1-step-lag case. The real world has, however, OOSMs with arbitrary lag. Subsequently, the suboptimal algorithm was extended to the case of an arbitrary (multistep) lag, but the resulting algorithm required a significant amount of storage. The present work shows how the 1-step-lag algorithms can be generalized to handle an arbitrary (multistep) lag while preserving their main feature of solving the update problem without iterating. This leads only to a very small (a few percent) degradation of MSE performance. The incorporation of an OOSM into the data association process is also discussed. A realistic example with two GMTI radars is presented. The consistency of the proposed algorithm is also evaluated and it is found that its calculated covariances are reliable.     

Novel Nonlinear Filtering & Prediction Method for Maneuvering Target Tracking

A new nonlinear filtering and prediction (NFP) algorithm with input es?imation is proposed for maneuvering target tracking. In the proposed method, the acceleration level is determined by a decision process, where a least squares (LS) estimator plays a major role in detecting target maneuvering within a sliding window. We first illustrate that the optimal solution to minimize the mean squared error (MSE) must consider a trade-off between the bias and error variance. For the application of target tracking, we then derive the MSE of target positions in a closed form by using orthogonal space decompositions. Then we discuss the NFP estimator, and evaluate how well the approach potentially works in the case of a set of given system parameters. Comparing with the traditional unbiased minimum variance filter (UMVF), Kalman filter, and interactive multiple model (IMM) algorithms, numerical results show that the newly proposed NFP method performs comparable or better in all scenarios with significantly less computational requirements.     

利用不机动单站DOA、TOA的测量实现被动定位与跟踪

 可以利用多个被动探测器或单个机动的探测器实现对运动辐射源的被动定位与跟踪。本文介绍的是一种利用非机动的单个探测器所获得的DOATOA测量值,对运动辐射源实现被动定位与跟踪的新方法。在二维平面内对运动辐射源使用这种方法,通过计算机仿真显示可以获得收敛而精确的跟踪数据。     

Reducing geometric dilution of precision using ridge regression

A signal processing technique is proposed for improving position-fix navigation system accuracy performance when the geometry of the navigation landmarks (e.g. sensors) are nearly collinear. In the navigation literature, the accuracy degradation associated with a nearly collinear measure geometry is termed the geometric dilution of precision (GDOP). Its presence causes not only the variance of the position estimates to be highly inflated but also any bias terms which may be present in the model. Since a nearly collinear predictor matrix is mathematically equivalent to GDOP, it is proposed to use the ridge regression technique in a navigation signal processor. A position-fix algorithm based on ridge regression reduces the bias and variance inflation caused by GDOP and the overall mean-squared position error as well. Ridge regression contains the GDOP-sensitive least-mean-square (LMS) estimator as a special case. Even with a matched model, GDOP can inflate the mean-square error (MSE) of the ordinary least-squares estimator, whereas the ridge regression technique chooses a suitable biased estimator that will reduce the MSE, which is the main goal. The ridge concept is extended to include GDOP-amplified bias errors. A simple range/range navigation system is analyzed to illustrate the underlying principles of ridge regression     

The Performance of Adaptive Array in a Dynamic Environment

The performance of adaptive arrays for jamming cancellation isdegraded due to relative angular movements between the arrayelements and the signal sources in space. We consider small angularvariations about the steady state of a linear (power inversion) arraycaused by either 1) a discrete angular jump, or 2) a continuousangular movement with constant angular velocity. Performancedegradation is assessed in terms of the resulting increase in the totaloutput power relative to its steady-state (minimum) value. Theresults reported are independent of the algorithm used to minimizethe output mean square error (MSE) and are supported bycomputer simulations.