Tracking multiple targets with correlated measurements andmaneuvers

The problem of tracking N targets with correlation in both measurement and maneuver statistics is solved by transforming to a coordinate frame in which the N targets are decoupled. For the case of N identical targets, the decoupling is shown to coincide with a transformation to a set of nested center-of-mass coordinates. Absolute and differential tracking accuracies are compared with suboptimal results to show the improvement that is achieved by properly exploiting the correlation between targets     

Finding the best set of K paths through a trellis withapplication to multitarget tracking

A solution is presented to the problem of finding the best set of K completely unmerged paths through a trellis with M _i⩾K states at depth i in the trellis, i=0, 1, 2, . . ., N. Here, `best set' means that the sum of the metrics of all K paths in the set is minimized, and `completely unmerged' means that no two paths pass through a common state. The solution involves using the Viterbi algorithm on an expanded trellis. This result is then used to separate the tracks of K targets optimally in a simplified model of a multitarget radar system. The model includes measurement errors and false alarms, but it does not include the effects of missing detections or merged measurements     

Convergence bounds of an SMI/Gram-Schmidt canceler in colored noise

The performance of the sampled matrix inversion (SMI) adaptive algorithm in colored noise is investigated using the Gram-Schmidt (GS) canceler as an analysis tool. Lower and upper bounds of average convergence are derived, indicating that average convergence slows as the input time samples become correlated. When the input samples are uncorrelated, the fastest SMI algorithm convergence occurs. When the input samples are correlated then the convergence bounds depend on the number of channels N, the number of samples per channels K , and the eigenvalues associated with K×K correlation matrix of the samples in a given channel. This matrix is assumed identical for all channels     

Efficient algorithms for finding the K best paths througha trellis

A set of algorithms is presented for finding the best set of K mutually exclusive paths through a trellis of N nodes, with worst-case computation time bounded by N³log n for a fixed-precision computation. The algorithms are based on a transformation of the K-path trellis problem into an equivalent minimum-cost network flow (MCNF) problem. The approach allows the application of efficient MCNF algorithms, which can obtain optimal solutions orders of magnitude faster than the algorithm proposed by J.K. Wolf et al. (1989). The resulting algorithms extend the practicality of the trellis formulation (in terms of required computations) to multiobject tracking problems with much larger numbers of targets and false alarms. A response by Wolf et al. is included     

Binary integration of fluctuating targets

The detection performance of a binary integrator (M-out-of-N detector) against nonfluctuating, slowly fluctuating, and quickly fluctuating targets is given. Since the solution for the slowly fluctuating target is numerically intensive, a simpler approximate solution is developed. This approximation is very accurate and is valid even when the noise power varies from pulse to pulse within a single antenna scan     

Periodic ambiguity function of CW signals with perfect periodicautocorrelation

A periodic ambiguity function (PAF) is discussed which describes the response of a correlation receiver to a CW signal modulated by a periodic waveform, when the reference signal in the receiver is constructed from an integral number N, of periods T, of the transmitted signal. The PAF is a generalization of the periodic autocorrelation function, to the case of non-zero Doppler shift. It is shown that the PAF of N periods is obtained by multiplying the PAF of a single period by the universal function sin(Nπν T)/N sin(πνT), where ν is the Doppler shift, to phase-modulated signals which exhibit perfect periodic autocorrelation when there is no Doppler shift. The PAF of these signals exhibits universal cuts along the delay and Doppler axes. These cuts are functions only of t, N and the number M, the modulation bits in one period     

Adaptive canceler and pulse compressor interactions

Performance results for the sidelobe level of a compressed pulse that has been preprocessed through an adaptive canceler are obtained. The adaptive canceler is implemented using the sampled matrix inversion algorithm. Because of finite sampling, the quiescent compressed pulse sidelobe levels are degraded due to the preprocessing of the main channel input data stream (the uncompressed pulse) through an adaptive canceler. It is shown that if N is the number of input canceler channels (main and auxiliaries) and K is the number of independent samples per channel, then K/N can be significantly greater than one in order to retain sidelobes that are close to the original quiescent sidelobe level (with no adaptive canceler). Also it is shown that the maximum level of degradation is independent of whether pulse compression occurs before or after the adaptive canceler if the uncompressed pulse is completely contained within the K samples that are used to calculate the canceler weights. This same analysis can be used to predict the canceler noise power level that is induced by having the desired signal present in the canceler weight calculation     

Multiple target tracking using products of position measurements

The continued development of the symmetric measurement equation (SME) filter for track maintenance in multiple target tracking (MTT) is considered, focusing on the case in which the SMEs are generated by forming sums of products of the original position measurements. The SME filter is developed for the case of N targets whose motions consist of random perturbations about constant-velocity trajectories. It is assumed that measurements of x-coordinate positions are available, and that the number of measurements is equal to the number of targets. Various analytical properties of the SME filter are studied. It is shown that under a very weak condition, the estimation error equation is locally exponentially stable. The performance of the SME filter is investigated by comparing it with an optimal (minimum-variance) estimator and by generating a computer simulation in the six-target case     

CFAR data fusion center with inhomogeneous receivers

Detection systems with distributed sensors and data fusion are increasingly used by surveillance systems. A system formed by N inhomogeneous constant false alarm rate (CFAR) detectors (cell-averaging (CA) and ordered statistic (OS) CFAR detectors) is studied. A recursive formulation of an algorithm that permits a fixed level of false alarms in the data fusion center is presented, to set the optimum individual threshold levels in the CFAR receivers and the optimum `K out of N' decision rule in order to maximize the total probability of detection. The algorithm also considers receivers of different quality or with different communication channel qualities connecting them with the fusion center. This procedure has been applied to several hypothetical networks with distributed CA-CFAR and OS-CFAR receivers and for Rayleigh targets and interference, and it was seen that in general the fusion decision OR rule is not always the best     

Cascaded detector for multiple high-PRF pulse Doppler radars

A postdetection design methodology for a multiple high-pulse-repetition frequency (PRF) pulse Doppler radar has been developed. The postdetection processor consists of an M out of N detector where range and target ambiguities are resolved, followed by a square-law detector which enhances the minimum signal-to-noise (S/N) power-ratio per pulse burst performance. For given probabilities of false alarm and detection, formulas are derived from which the three thresholds associated with the cascaded detector can be found. Fundamental tradeoffs between the minimum S/N required, number of ghosts, and the number of operations (NOPs) that the cascaded detector must perform are identified. It is shown that the NOPs and the number of ghosts increase and the minimum S/N required decreases as the binary M out of N detector passes more detections to the square-law detector     

Optimal speckle reduction in polarimetric SAR imagery

Speckle is a major cause of degradation in synthetic aperture radar (SAR) imagery. With the availability of fully polarimetric SAR data, it is possible to use the three complex elements (HH, HV, VV) of the polarimetric scattering matrix to reduce speckle. The optimal method for combining the elements of the scattering matrix to minimize image speckle is derived, and the solution is shown to be a polarimetric whitening filter (PWF). A simulation of spatially correlated, K-distributed, fully polarimetric clutter is then used to compare the PWF with other, suboptimal speckle-reduction methods. Target detection performance of the PWF, span, and single-channel |HH|² detectors is compared with that of the optimal polarimetric detector (OPD). A novel, constant-false-alarm-rate (CFAR) detector (the adaptive PWF) is as a simple alternative to the OPD for detecting targets in clutter. This algorithm estimates the polarization covariance of the clutter, uses the covariance to construct the minimum-speckle image, and then tests for the presence of a target. An exact theoretical analysis of the adaptive PWF is presented; the algorithm is shown to have detection performance comparable with that of the OPD     

Tracking multitarget in cluttered environment

A method for multitarget tracking and initiating tracking in a cluttered environment is proposed. The algorithm uses a sliding window of length uT (T is the sampling time) to keep the measurement sequence at time k. Instead of solving a large problem, the entire set of targets and measurements is divided into several clusters so that a number of smaller problems are solved independently. When a set of measurements is received, a new set of data-association hypotheses is formed for all the measurements lying in the validation gates within each cluster from time K-u+1 to K. The probability of each track history is computed, and, choosing the largest of these histories, the target measurement is updated with an adaptive state estimator. A covariance-matching technique is used to improve the accuracy of the adaptive state estimator. In several examples, the algorithm successfully tracks targets over a wide range of conditions     

Detection performance in K-distributed and correlatedRayleigh clutters

Performance prediction for a detection system employing noncoherent integration is carried out for a chi-square family of fluctuating targets in K-distributed clutter plus noise. The detection performance for Swerling 11 targets in the K-distributed clutter plus noise is compared with that in exponentially correlated Rayleigh clutter. The results show that the performance prediction based on N pulses integrated in clutter plus noise using the K-distributed clutter model may be approximately equivalent to that using the exponentially correlated Rayleigh-distributed clutter model     

A parallel square-root algorithm for modified extended Kalmanfilter

A parallel square-root algorithm and its systolic array implementation are proposed for performing modified extended Kalman filtering (MEKF). The proposed parallel square-root algorithm is designed based on the singular value decomposition (SVD) and the Faddeev algorithm, and a very large scale integration (VLSI) systolic array architecture is developed for its implementation. Compared to other square root Kalman filtering algorithms, the proposed method is more numerically stable. The VLSI architecture described has good parallel and pipelining characteristics in applying to the MEKF and achieves higher efficiency. For n-dimensional state vector estimations, the proposed architecture consists of O(2n²) processing elements and uses O ((s+17)n) time-steps for a complete iteration at each instant, in contrast to the complexity of O((s+6) n³) time-steps for a sequential implementation, where s≈log n     

False alarm control using Doppler estimation

A technique which uses maximum-likelihood estimates (MLEs) of target Doppler and target amplitude is developed for rejecting clutter residues. Multiple estimates are made and consistency checks are applied to the estimates. Simulation results indicate that for large clutter-to-noise ratios (C/N⩾55 dB) the probability of false alarm from clutter residues is reduced from 1.0 to below 0.01     

Low-altitude target detection by coastline operated marine radar

Relevant to a Richian family of fluctuating targets with a composite background of sea-plus-land clutter, the performance prediction of a radar operating in near-coastal regions is elucidated by assuming noncoherent integration of the pulses. Considering the dominance of land clutter, a modified K-distributed statistic is indicated for the overall clutter envelope; and the corresponding probability of false alarm and probability of detection are deduced for fixed threshold detection (s) based on N pulses integrated in the presence of the sea-plus-land clutter and the noise. Even when the target offers a dominant scattered echo, the worst situations of the land clutter affecting the detection performance are indicated     

On combining multidimensional target location ellipsoids

N.M. Blachman (1989) presented a method for combining the two-dimensional target location ellipses using the equation of an ellipse. A method is presented for combining N-dimensional location ellipsoids by using the location covariance matrix. The general result is formulated in simple matrix form. The results presented contain Blachman's results as a special case. Examples for both two and three dimensions are presented     

Optimal serial distributed decision fusion

The problem of distributed detection involving N sensors is considered. The configuration of sensors is serial in the sense that the Jth sensor decides using the decision it receives along with its own observation. When each sensor uses the Neyman-Pearson test, the probability of detection is maximized for a given probability of false alarm, at the Nth stage. With two sensors, the serial scheme has a performance better than or equal to the parallel fusion scheme analyzed in the literature. Numerical examples illustrate the global optimization by the selection of operating thresholds at the sensors     

Distribution of general noncentral positive definite quadratic formin K-dimensions

The authors present a series solution using Hermite polynomials to the long-standing problem of computing the probability P that positive definite noncentral quadratic form d(x) of a Gaussian random vector x∈R satisfies d( x)⩽r² for any given r∈R. This problem has wide applications in radar, tracking, air traffic control, etc. The fast-converging series solution presented is very accurate and can be performed rapidly using the recursion relations for Hermite polynomials     

Performance of receivers with linear detectors

The false-alarm and detection probabilities of a receiver summing M independent outputs of a linear detector are calculated by numerical saddlepoint integration. The saddlepoint approximation is also considered. Both constant-amplitude and Rayleigh-fading signals are treated, and the relative efficiency of the quadratic and the linear detectors for these is calculated for a broad range of values of M . The numerical integration method is the more efficient, the smaller the false-alarm probability or the false-dismissal probability, that is, under just those conditions for which the terms in the Gram-Charlier series oscillate most violently and the series becomes least reliable. The simpler saddlepoint approximation yields values that in those same regions have been found close enough to the exact probabilities to be adequate for most engineering purposes. The larger the number M of samples, the more efficient methods are