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     

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     

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     

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     

Optimum gain control for A/D conversion using digitized I/Q data inquadrature sampling

Calculation of optimum gain for minimum distortion due to A/D (analog-to-digital) conversion requires the estimation of the input signal strength. To use a common AGC (automatic gain control) for both the I/Q (in-phase and quadrature) signals, it is efficient to estimate the input signal strength using the quantized A/D output from both channels. Assuming a Gaussian input, the relationship between σ of the input of the A/D converter and E(|x|+|y|) and E(max(|x|,|y|)+1/2 min (|x|,|y|)) for t quantized I/Q output x and y is derived. Numerical results obtained using the derived expression and the statistical data obtained through simulation show excellent agreement. It is concluded that, because of its simplicity, the cubic equation obtained by fitting the numerical results should be useful     

Nonparametric rank detectors on quantized radar video signals

Uniform randomization of ties is required for defining distribution-free ranks of independent and identically distributed quantized samples. Formulas of rank probabilities are given and applied to radar detection under quantized video samples. For some detectors, and assuming Gaussian noise, the asymptotic loss L_∞(q) is calculated versus the normalized quantization step q, and the loss L(q) is estimated by Monte Carlo simulations. Both of these resulted in monotonic functions of q (0<q<1.1) that are independent of the other parameters. Furthermore, L(q)≈L_∞(q )⩽0.45 dB, as q<0.8. The quantization step q is normalized with respect to the noise standard deviation     

Detectability of signals in Laplace noise

Closed-form expressions for the false-alarm and detection probabilities attained by the optimum and the linear detectors of a positive signal in n independent samples of noise having a bilateral exponential or Laplace distribution require lengthy computation when n is large, and those for the optimum detector suffer from round-off error because their terms alternate in sign. It is shown how the method of saddlepoint integration can be conveniently applied to compute these probabilities, and numerical comparisons of the accuracies of the methods are presented. The relative efficiency of the two detectors is calculated as a function of n and found to approach the asymptotic value of 2 very slowly     

Systolic array implementation of second-order Kalman filter

The approach for design of parallel implementation of the second-order Kalman filter using systolic array processors is proposed. The original time complexity of more than 27n³ for sequential processing can be reduced to only 21n+10r-4 in the proposed parallel implementation. The improvement is gained from the utilization of interarray pipelining, concurrent processing and substitution forwarding techniques. The approach is especially suitable for real-time application of the Kalman filter in solving second-order linear system     

CCRS C/X-airborne synthetic aperture radar: An Rand D tool for the ERS-1 time frame

The airborne synthetic-aperture radar (SAR) system developed for the Canada Centre for Remote Sensing (CCRS) is described. It consists of two radars, at C-band and X-band. Each radar incorporates the following features: dual-channel receivers and dual-polarized antennas; a high quality, 7-look, real-time processor; a sensitivity time control for range-dependent gain control; a motion-compensation system for antenna steering in azimuth and elevation; and baseband I and Q signal phase rotation. The system also uses a high-power transmitter with a low-power back-up. The SAR maps to either side of the aircraft, at high or low resolution, at incidence angles which in high resolution span 0° to 80°. Radar operating parameters, data products, key specifications and the motion compensation scheme used are presented. Properties of the real-time imagery are discussed and examples of C-band SAR data in the three operating modes are given     

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     

Mismatched-filter design for periodic binary phased signals

A mismatched filter design is introduced whereby the receiver filter coefficients are optimized with the property that all sidelobes of the cross-correlation function (CCF) are zero. Some results are given for binary sequences and the associated mismatched filters. For each length N>2 there exist at least one sequence and an associated mismatched filter with an ideal impulselike CCF, which means that there are no sidelobes at all in the receiver output signal     

Identification of rational transfer function from frequencyresponse sample

A method for identifying a transfer function, H(z)=A(z)/B(z), from its frequency response values is presented. Identifying the transfer function involves determining the unknown degrees and coefficients of the polynomials A(z) and B( z), given the frequency response samples. The method for finding the parameters of the transfer function involves solving linear simultaneous equations only. An important aspect of the method is the decoupled manner in which the polynomials A(z) and B(z) are determined. The author presents two slightly different derivations of the linear equations involved, one based on the properties of divided differences and the other using Vandermonde matrices or, equivalently, Lagrange interpolation. A matrix synthesized from the given frequency response samples is shown to have a rank equal to the number of poles in the system     

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     

Estimation of radar detection and false alarm probability

The accuracy with which detection and false alarm probabilities can be estimated with a limited amount of measured radar data is addressed. A simple simulation method for estimating the statistical performance of a radar detection system is presented. Confidence limits and a rule of thumb for accuracy for the estimated probabilities are presented along with procedures for calculating them. It is concluded that the minimum value of N used in a detection radar signal simulation should be 10/P_FA when the simple simulation method is used, where P_FA is the probability of false alarm, and that a value closer to 100/P _FA is preferable     

Detection probabilities for correlated Rayleigh fading signals

Both the method of saddlepoint integration and its associated saddlepoint approximation are applied to calculating the probability of detecting correlated Rayleigh-fading signals in Gaussian noise by means of a detector that integrates M samples of the output of a quadratic rectifier. The quadrature components of the signal samples are modeled as an autoregressive moving-average process, and specific results are exhibited for a first-order Markov process. By these methods the fluctuation loss can be computed for much larger values of M and for larger values of the detection probability than previously. Values calculated by the saddlepoint approximation prove to be close enough to the exact values to be useful over a broad range of signal parameters     

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     

Continuous-time ECV and ECA track filters with colored measurementnoise

Analytical solutions are presented for two continuous-time track filters with correlation in both the measurement errors and the target maneuvers. The measurement error model is a first-order Markov process, characterized by a correlation time (γ) and a position error variance (r). The target motion is described by one of two dynamic models: the exponentially correlated velocity (ECV) model, and the exponentially correlated acceleration (ECA) model. For both models, the steady-state velocity variance exhibits resonant behavior. The analytical solutions include as special cases several previously published results     

The autocorrelation of m-sequences over nonprime finitefields

Nonbinary m-sequences (maximal length sequences) for spread-spectrum communication systems that have a two-level autocorrelation are presented. The autocorrelation function of an m -sequence over the Galois field of q elements GF(q), where q=p^k, for p a prime and k an integer greater than 1, is developed and shown to be bilevel when the elements of GF(q) are expressed as elements of a vector space over the pth roots of unity     

CW alternatives to the coherent pulse train-signals and processors

Continuous wave (CW) signals phase modulated by a periodic waveform, and their corresponding receivers, are discussed. The combined response in delay and Doppler is almost identical to the (ideal) response of the coherent pulse train. The receivers are matched to an integral number N of modulation periods of the transmitted signal. CW implies a duty cycle of 100%. However, the signal duration need not be longer than N+2 periods. The CW signals have the advantage that their peak power is equal to the average power. Their disadvantages are more complicated receiver processing and the need for two antennas     

A divide and conquer approach to least-squares estimation

The problem of estimating parameters &thetas; which determine the mean μ(&thetas;) of a Gaussian-distributed observation X is considered. It is noted that the maximum-likelihood (ML) estimate, in this case the least-squares estimate, has desirable statistical properties but can be difficult to compute when μ(&thetas;) is a nonlinear function of &thetas;. An estimate formed by combining ML estimates based on subsections of the data vector X is proposed as a computationally inexpensive alternative. The main result is that this alternative estimate, termed here the divide-and-conquer (DAC) estimate, has ML performance in the small-error region when X is appropriately subdivided. As an example application, an inexpensive range-difference-based position estimator is derived and shown by means of Monte-Carlo simulation to have small-error-region mean-square error equal to the Cramer-Rao lower bound