非均匀噪声稀疏均匀圆阵的二维DOA估计

针对在机载雷达、通信等领域有着广泛应用的均匀圆阵(UCA),研究了非均匀噪声下稀疏均匀圆阵的二维波达方向(DOA)估计.首先采用改进的相位模式方法构造非均匀噪卢稀疏均匀圆阵的波柬空间似然函数;之后,在分析非均匀噪声稀疏均匀圆阵的波束空间似然函数特点的基础上,修改了Burg的逆迭代算法以适应稀疏均匀圆阵下非均匀噪声自相关...     

基于协方差矩阵高阶幂的二维DOA估计新算法

 针对稀疏分解方法进行均匀圆阵(UCA)的二维波达方向(DOA)估计运算复杂度大的问题,提出了一种基于协方差矩阵高阶幂稀疏分解的二维DOA估计新算法。该算法首先利用协方差矩阵高阶幂无需进行特征值分解和信源数估计的特性,构建了协方差矩阵高阶幂的稀疏分解向量;然后运用粒度分层思想,构造了粗区域估计和细方位估计的分层多粒度的快速分解模型,分层字典的长度大大减少,在保持估计精度的前提下,算法运算时间远小于现有的恒定冗余字典的稀疏分解方法,从而解决了基于稀疏分解的圆阵二维DOA估计问题。论文提出的算法与二维MUSIC算法相比,估计精度高,且能满足对相干信号的估计。仿真结果验证了算法的有效性和可行性。     

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.     

Adapitive waveform design for distributed OFDM MIMO radar system in multi-target scenario

In order to improve detection and estimation performance of distributed Orthogonal-Frequency-Division Multiplexing (OFDM) Multiple-Input Multiple-Output (MIMO) radar system in multi-target scene, we propose a novel approach of Adaptive Waveform Design (AWD) based on a constrained Multi-Objective Optimization (MOO). The sparse measurement model of this radar system is derived, and the method based on decomposed Dantzig selectors is applied for the sparse recovery according to the block structures of the sparse vector and the system matrix. An AWD approach is proposed, which optimizes two objective functions, namely minimizing the upper bound of the recovery error and maximizing the weakest-target return, by adjusting the complex weights of the emitting waveform amplitudes. Several numerical simulations are provided and their results show that the detection and estimation performance of the radar system is improved significantly when this MOO-based AWD approach is applied to the distributed OFDM MIMO radar system. Especially, we verify the effectiveness of our AWD approach when the available samples are reduced severally and the technique of compressed sensing is introduced.     

Direction of arrival tracking below the ambiguity threshold

We present an algorithm for direction-of-arrival (DOA) tracking that allows operation below the ambiguity threshold of the direction-finding (DF) system. Using multiple target tracking techniques, the algorithm turns the most likely DOAs of each measurement into multiple potential tracks and then selects the true track as that with the maximum cumulative likelihood. The improvement offered by the algorithm, namely the extension of the ambiguity-free signal-to-noise ratio (SNR) domain, is demonstrated in several simulated experiments using several array structures, including a sparse array and a uniform linear array     

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.     

一种稀疏阵列下的二维DOA估计方法

研究了稀疏阵列下二维波达方向（DOA）的估计问题,提出一种基于不动点迭代的空间谱估计（FPC-MUSIC）算法。首先建立基于矩阵填充的DOA估计信号模型,并验证该信号模型满足零空间性质（NSP）,其次通过不动点迭代算法将稀疏阵列信号恢复为完整信号,最后利用恢复信号估计二维DOA。该算法可在稀疏阵列下大幅度降低谱估计平均副瓣,在大幅度降低阵元数的同时具有较高的估计精度。计算机仿真表明：FPC-MUSIC算法可在稀疏阵列下准确估计二维DOA,验证了该算法的有效性和优越性。     

Multipath mitigation method based on Gaussian mixture model in RF relative measurement

Radio Frequency (RF) technology represents a high-precision relative navigation solution that has significant potential for application to earth-orbiting satellites. In precision applications, multipath errors dominate the total error because observables, which are used to estimate carrier-phase integer ambiguity, are not always subject to a Gaussian distribution when dual-frequency ambiguity estimation methods are used in the presence of multipath. As it has been shown that ranging observables obey a Gaussian mixture distribution, this study proposes improvements to the accuracy of estimation based on multipath mitigation founded on the Gaussian mixture model. To this end, such a model is created for integer ambiguity resolution in the presence of multipath, using which the theoretical error in dual-frequency ambiguity estimation is derived. Expectation Maximization (EM), which aids dual-frequency ambiguity estimation, is subsequently proposed to reduce the effect of multipath errors. Finally, two experimental scenarios are implemented to test the performance of the proposed method. The results show that EM-aided dual-frequency ambiguity estimation reduces the range error to approximately 20% in comparison with simple dual-frequency ambiguity estimation. Therefore the proposed technique is effective for multipath mitigation in RF relative measurement.     

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

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

Asymptotic Performance of Sparse Signal Detection Using Convex Programming Method

The detection of sparse signals against background noise is considered. Detecting signals of such kind is difficult since only a small portion of the signal carries information. Prior knowledge is usually assumed to ease detection. In this paper, we consider the general unknown and arbitrary sparse signal detection problem when no prior knowledge is available. Under a Neyman-Pearson hypothesis-testing framework, a new detection scheme is proposed by combining a generalized likelihood ratio test (GLRT)-like test statistic and convex programming methods which directly exploit sparsity in an underdetermined system of linear equations. We characterize large sample behavior of the proposed method by analyzing its asymptotic performance. Specifically, we give the condition for the Chernoff-consistent detection which shows that the proposed method is very sensitive to the 2 norm energy of the sparse signals. Both the false alarm rate and the miss rate tend to zero at vanishing signal-to-noise ratio (SNR), as long as the signal energy grows at least logarithmically with the problem dimension. Next we give a large deviation analysis to characterize the error exponent for the Neyman-Pearson detection. We derive the oracle error exponent assuming signal knowledge. Then we explicitly derive the error exponent of the proposed scheme and compare it with the oracle exponent. We complement our study with numerical experiments, showing that the proposed method performs in the vicinity of the likelihood ratio test (LRT) method in the finite sample scenario and the error probability degrades exponentially with the number of observations.     

An operational system implementation of the ESPIRIT DF algorithm

The practical implementation of the ESPRIT algorithm into direction-finding system architecture is considered. In particular, the problems of phase ambiguity resolution for extended array separations and multidimensional azimuth/elevation estimation are addressed. Generalized solutions for these problems are proposed, and a bearing quality indicator that provides a measure of credibility associated with each angle-of-arrival estimate is developed. Plots of computer simulated performance are presented in order to assess the probability of detection, of failure to alarm, and of false alarm as a function of azimuthal separation and period of signal observation     

Model-based multifrequency array signal processing for low-angletracking

Tracking low-altitude targets over the sea is problematic because of interference between the direct and reflected signal. Standard monopulse trackers can experience large errors because of multipath maximum likelihood estimation (MLE) has been used to more accurately estimate the target height in the presence of multipath MLE is a model-fitting technique where the model parameters are chosen to maximize the likelihood function. It is shown that the type of observation model has a large effect on performance. Tracking performance is compared using three different observation models employing varying amounts of a priori information. Results are presented for different array sizes: eight and 32-element arrays and two-element subarrays typical of phase monopulse. Performance is compared with that of standard techniques such as Fourier beamforming and phase monopulse     

Resolution and synthetic aperture characterization of sparse radar arrays

The concept of radar satellite constellations, or clusters, for synthetic aperture radar (SAR), moving target indicator (MTI), and other radar modes has been proposed and is currently under research. These constellations form an array that is sparsely populated and irregularly spaced; therefore, traditional matched filtering is inadequate for dealing with the constellation's radiation pattern. To aid in the design, analysis, and signal processing of radar satellite constellations and sparse arrays in general, the characterization of the resolution and ambiguity functions of such systems is investigated. We project the radar's received phase history versus five sensor parameters: time, frequency, and three-dimensional position, into a phase history in terms of two eigensensors that can be interpreted as the dimensions of a two-dimensional synthetic aperture. Then, the synthetic aperture expression is used to derive resolution and the ambiguity function. Simulations are presented to verify the theory.     

基于均匀圆阵的信号二维方向角高精度估计

基于均匀圆阵(UCA),提出一种信号二维方向角的高精度估计方法。首先,充分利用了信号的空、时域信息,构造时空旋转矩阵,从而达到多个信源分离的目的。其次,直接利用UCA的阵列流形,采用最小二乘法估计信号二维方向角。最后,给出解相位整周期性模糊的方法。计算机仿真表明了此方法具有估计精度高,对阵元的幅相误差鲁棒性强等特点。     

On cross-ambiguity properties of Welch-Costas arrays

We discuss cross-ambiguity properties of a specific family of Costas arrays called Welch-Costas (W-C) arrays. These properties are of interest in multiuser radar and sonar system, especially since Costas arrays are known to possess ideal auto-ambiguity functions. The theory of W-C arrays is reviewed. It is then proved that only pairs of W-C arrays can have at most two hits in their cross-ambiguity function (best possible case). The maximum number of hits in the cross-ambiguity functions of a family of W-C arrays is shown to be a function of the number of W-C arrays in the family. The upper bound on the number of hits in the cross ambiguity functions for a family of W-C arrays is also derived. Specific examples of how reducing the number of W-C arrays improves the cross-ambiguity properties are given for various types of prime numbers     

Small Aperture Angle Measurement for Active Echo Location Systems

The problem of bearing estimation for active systems is examined from the point of view of the generalized wideband ambiguity function (GAF). The maximum likelihood (ML) estimator is derived and its local and global properties are discussed. A structure is proposed which searches the three-dimensional ambiguity surface in two stages first, in range-Doppler, and then, in bearing with the goal of reducing search complexity when utilizing highly resolvent waveforms. Comparisons are made between the ML estimators and structures utilizing phase information which generate closed form estimators. The beneficial results of full bandwidth utilization are discussed in terms of both local and global properties of the GAF.     

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

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

Comparison of various periodograms for sinusoid detection andfrequency estimation

With the advent of the fast Fourier transform (FFT) algorithm, the periodogram and its variants such as the Bartlett's procedure and Welch method, have become very popular for spectral analysis. However, there has not been a thorough comparison of the detection and estimation performances of these methods. Different forms of the periodogram are studied here for single real tone detection and frequency estimation in the presence of white Gaussian noise. The threshold effect in frequency estimation, that is, when the estimation errors become several orders of magnitude greater than the Cramer-Rao lower bound (CRLB), is also investigated. It is shown that the standard periodogram gives the optimum detection performance for a pure tone while the Welch method is the best detector when there is phase instability in the sinusoid. As expected, since the conventional periodogram is a maximum likelihood estimator of frequency, it generally provides the minimum mean square frequency estimation errors     

Range, radial velocity, and acceleration MLE using radar LFM pulsetrain

An efficient implementation of the maximum likelihood estimator (MLE) is presented for the estimation of target range, radial velocity, and acceleration when the radar waveform consists of a wideband linear frequency modulated (LFM) pulse train. Analytic properties of the associated wideband ambiguity function are derived; in particular the ambiguity function, with acceleration set to zero, is derived in closed form. Convexity and symmetry properties of the ambiguity function over range, velocity, and acceleration are presented; these are useful for determining region and speed of convergence for recursive algorithms used to compute the MLE. In addition, the Cramer-Rao bound (CRB) is computed in closed form which shows that the velocity bound is decoupled from the corresponding bounds in range and acceleration. A fast MLE is then proposed which uses the Hough transform (HT) to initialize the MLE algorithm. Monte Carlo simulations show that the MLE attains the CRB for low to moderate signal-to-noise depending on the a priori estimates of range, velocity, and acceleration