General direction-of-arrival estimation: a signal subspace approach

A high-resolution algorithm is presented for resolving multiple incoherent and coherent plane waves that are incident on an array of sensors. The incident sources can be a mixture of narrowband and broadband sources, and, the geometry of the array is unrestricted. The algorithm makes use of a fundamental property possessed by those eigenvectors of the array spectral density matrix that are associated with eigenvalues that are larger than the sensor noise level. Specifically, it is shown that these eigenvectors can each be represented as linear combinations of the steering vectors identifying the incident plane waves. This property is then used to solve the important special cases of incoherent sources incident on a general array and coherent sources incident on an equispaced linear array. Simulation results are presented to illustrate the high-resolution performance achieved with this approach relative to that obtained with MUSIC and spatial smoothed MUSIC in which the coherent-signal-subspace focusing method is used     

Comparison of signal estimation using calibrated and uncalibrated arrays

We consider the problem of separating and estimating the waveform of superimposed signals received by an array of sensors. If the array is well calibrated it is possible to first estimate the directions of arrival (DOA) of the signals and then use this information to separate the signals. When the array is not calibrated, but the array elements have the same unknown gain pattern, up to an unknown multiplicative factor and the phases of the elements are arbitrary and unknown, it is possible to estimate the array steering vectors and then use this information for signal estimation. We compare the quality of the estimated signals in the calibrated case with the quality of the estimated signals in the uncalibrated case, in terms of the output signal-to-interference ratio (SIRO) and output signal-to-noise ratio (SNRO).     

Efficient wideband signal parameter estimation using a radon-ambiguity transform slice

A novel efficient technique based on a single slice Radon-ambiguity transform (RAT) for time-delay and time-scale estimation is proposed. The proposed approach combines the narrowband cross-ambiguity function (NBCAF), the wideband cross-ambiguity function (WBCAF), and a single slice RAT to estimate multiple target parameters in noisy environments. The square modulus of Gaussian-enveloped linear frequency modulated (GLFM) signals has high-energy centrality in the ambiguity plane. Its peaks in the NBCAF fall along nearly straight lines whose slopes depend on the Doppler rates of the moving targets. These lines could be effectively detected by computing the entire Radon transform of the NBCAF for all possible angles; however, it is a computationally intensive procedure. It is shown that without calculating the entire RAT, it is possible to estimate target parameters using only a single slice of the RAT, i.e., using an appropriate projection of the NBCAF. It is demonstrated that the proposed method can successfully separate overlapping targets efficiently. The efficiency is achieved due to fast Fourier transform (FFT)-bascd processing, use of a single slice of RAT, and the use of only one-dimensional (1-D) searches.     

Noise subspace techniques in non-gaussian noise using cumulants

We consider noise subspace methods for narrowband direction-of-arrival or harmonic retrieval in colored linear non-gaussian noise of unknown covariance and unknown distribution. The non-gaussian noise covariance is estimated via higher order cumulants and combined with correlation information to solve a generalized eigenvalue problem. The estimated eigenvectors are used in a variety of noise subspace methods such as multiple signal classification (MUSIC), MVDR and eigenvector. The noise covariance estimates are obtained in the presence of the harmonic signals, obviating the need for noise-only training records. The covariance estimates may be obtained nonparametrically via cumulant projections, or parametrically using autoregressive moving average (ARMA) models. An information theoretic criterion using higher order cumulants is presented which may be used to simultaneously estimate the ARMA model order and parameters. Third- and fourth-order cumulants are employed for asymmetric and symmetric probability density function (pdf) cases, respectively. Simulation results show considerable improvement over conventional methods with no prewhitening. The effects of prewhitening are particularly evident in the dominant eigenvalues, as revealed by singular value decomposition (SVD) analysis     

Linear FM signal parameter estimation from discrete-timeobservations

The authors consider the problem of estimating the parameters of a complex linear FM signal from a finite number of noisy discrete-time observations. An estimation algorithm is proposed, and its asymptotic (large sample) performance is analyzed. The algorithm is computationally simple, consisting of two fast Fourier transforms (FFTs) accompanied by one-dimensional searches for maxima. The variance of the estimates is shown to be close to the Cramer-Rao lower bound when the signal-to-noise ratio is 0 dB and above. The authors applied the algorithm to the problem of estimating the kinematic parameters of an accelerating target by pulse-Doppler radar. A representative test case was used to exhibit the usefulness of the algorithm for this problem, and to verify the analytical results by Monte Carlo simulations     

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     

基于共轭循环平稳特性的二维波达方向估计方法

 提出了一种具有阵列孔径扩展能力的共轭循环平稳二维波达方向（DOA）估计方法。该方法通过引入伪采样,充分利用了阵元输出之间共轭循环互相关函数的时间维信息,扩展了阵列有效孔径,提高了算法DOA估计精度,具有较好的低信噪比适应能力。该方法避免了最优时延选择问题,无需谱峰搜索,且其二维角度参数能够自动配对,具有较低的计算量。计算机仿真结果验证了该方法的有效性。     

Resolution of signal sources via spectral moment estimation

In different practical situations it is desired to estimate the number of signal sources and their positions in space or in frequency domain. The first problem is known as the detection or the order estimation and the second one as the resolution. For the resolution problem techniques such as nonlinear least squares (NLSM), high-order Yule-Walker method (HOYW), multiple signal classification (MUSIC), Pisarenko harmonic retrieval method, min-norm method, estimation of signal parameters by rotational invariance technique (ESPRIT), were proposed (Marple, 1987 and Stoica and Moses, 1997). All these high-resolution methods are based on the analysis of the signal covariance matrix. But the covariance matrix is not the only choice to represent the signal spectrum. In different applications (weather radars, synthetic aperture radar (SAR) signal processing, ultrasound imaging in medicine, atmospheric turbulence measurements) the signal spectrum can be modeled through its algebraic moments. Recently a number of efficient nonparametric methods have been proposed to estimate the algebraic spectral moments (Monakov, 1999). The presented paper is an attempt to solve the direction of arrival (DOA) problem via estimation of the algebraic spectral moments. A method proposed in the article is comparable in its accuracy with the MUSIC method. At the same time its computational burden is much lower. The method permits to estimate the signal power of sources easily to complete the full spectral line analysis. Additionally the method shows good robustness in situations when signal sources have noticeable spatial extend     

Angle and polarization estimation in a coherent signal environment

It is shown how a uniform linear array of crossed dipoles may be used with the ESPRIT algorithm and spatial smoothing techniques to estimate the arrival directions and polarizations of incoming coherent plane waves. Some examples showing typical performance are presented. One method of smoothing can be used where it is necessary to estimate both the arrival angles and polarizations of signals. Two other methods can be used when only the arrival angles are of interest     

欺骗环境下GNSS信号估计与定位修正技术

欺骗信号以其极强的隐蔽性使卫星导航接收机难以察觉并迅速定位到错误位置,严重影响了卫星导航的安全性。现有抗欺骗技术需要其他导航系统辅助来修正受欺骗影响的定位解算,针对该问题,本文提出了一种GNSS欺骗信号参数估计与辨识方法,能够在欺骗干扰环境下估计并辨识出真实信号所对应的伪距,进而解算出接收机真实位置。该方法通过研究欺骗干扰下接收机相关值模型,在信号跟踪阶段建立真实与欺骗双信号状态模型与基于九路相关器输出的观测模型,利用扩展卡尔曼滤波（EKF）估计真实信号与欺骗信号的伪码延时与信号相关幅值,进而获得真实与欺骗伪距,在定位解算阶段利用改进观测量残差检测方法辨识出真实与欺骗伪距,最终使用真实伪距定位获得真实位置。仿真结果表明对相对码延时介于0.3~0.9 chip之间且欺骗/真实信号幅度比介于1~5之间的隐蔽欺骗攻击,所提方法的码延时估计误差约0.1 chip,可有效估计真实信号与欺骗信号参数,辨识出真实伪距,并使被欺骗的定位结果重新回到真实位置结果,改善GNSS接收机抗欺骗能力,提高卫星导航安全性。     

The modeling and estimation of asynchronous multirate multisensor dynamic systems

An asynchronous data fusion problem based on a kind of multirate multisensor dynamic system is studied. The system is observed by multirate sensors independently, with the state model known at the finest scale. Under the assumption that the sampling rates of sensors decrease successively by any positive integers, the discrete dynamic system models are established based on each single sensor and an asynchronous multirate multisensor state fusion estimation algorithm is presented. Theoretically, the estimate is proven to be unbiased and the optimal in the sense of linear minimum covariance, the fused estimate is better than the Kalman filtering results based on each single sensor, and the accuracy of the fused estimate will decrease if any of the sensors' information is neglected. The feasibility and effectiveness of the algorithm are shown through simulations.     

DOA and polarization estimation via signal reconstruction with linear polarization-sensitive arrays

This paper addresses the problem of direction-of-arrival(DOA) and polarization estimation with polarization sensitive arrays(PSA), which has been a hot topic in the area of array signal processing during the past two or three decades. The sparse Bayesian learning(SBL) technique is introduced to exploit the sparsity of the incident signals in space to solve this problem and a new method is proposed by reconstructing the signals from the array outputs first and then exploiting the reconstructed signals to realize parameter estimation. Only 1-D searching and numerical calculations are contained in the proposed method, which makes the proposed method computationally much efficient. Based on a linear array consisting of identically structured sensors, the proposed method can be used with slight modifications in PSA with different polarization structures. It also performs well in the presence of coherent signals or signals with different degrees of polarization. Simulation results are given to demonstrate the parameter estimation precision of the proposed method.     

子空间辨识方法的改进

针对现代航空发动机的自适应控制和故障监控领域建立物理意义明确的系统动态模型的需要，基于数字子空间状态空间系统辨识（N4SID）方法，提出并推导了指定状态为量的子空间辨识方法，并在某些双转子涡喷发动机气动热力学模型上进行了仿真实验，说明通过指定状态变量的子空间辨识方法可以得到较好的辨识结果。     

Distributed decision fusion using empirical estimation

The problem of optimal data fusion in multiple detection systems is studied in the case where training examples are available, but no a priori information is available about the probability distributions of errors committed by the individual detectors. Earlier solutions to this problem require some knowledge of the error distributions of the detectors, for example, either in a parametric form or in a closed analytical form. Here we show that, given a sufficiently large training sample, an optimal fusion rule can be implemented with an arbitrary level of confidence. We first consider the classical cases of Bayesian rule and Neyman-Pearson test for a system of independent detectors. Then we show a general result that any test function with a suitable Lipschitz property can be implemented with arbitrary precision, based on a training sample whose size is a function of the Lipschitz constant, number of parameters, and empirical measures. The general case subsumes the cases of nonindependent and correlated detectors.     

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     

Modified input estimation technique using pseudoresiduals

An adaptive tracking filter for maneuvering targets is proposed using modified input estimation technique. Pseudoresiduals are defined using measurements and the velocity estimate at the hypothesized maneuver onset time. With the pseudoresiduals and a new target model representing transitions of nominal accelerations, a new input estimation method for tracking a maneuvering target is derived. Since the proposed detection technique is more sensitive to maneuvers than previous work, the shorter window length can be employed to detect and compensate target maneuvers. Also shown is that the tracking performance of the proposed filter is similar to that of interacting multiple model method (IMM) with 3 models, while computational loads of our method are drastically reduced     

Maneuver estimation using measurements of orientation

An image-based algorithm which provides an estimate of the current radial acceleration of a target is studied. Since the data in a single image frame provide information only on the orientation of the target, a sequence of frames must be processed to detect maneuvers. The estimation problem is parameterized in terms of natural groupings of measurement errors, and the influence of these errors on estimator fidelity is studied     

DOA estimation using one-bit quantized measurements

The effects 1-bit quantization of the input samples has on the direction-of-arrival (DOA) estimation accuracy are considered. The signal model assumes a single stochastic Gaussian point source that is embedded in white Gaussian noise (WGN). The inherent limitations governed by the extreme clipping of the input data are analyzed using the Cramer-Rao bound (CRB) that is derived for a two-sensor array. In addition, several estimators for the I-bit estimation are discussed. Numerical and analytical analyses of the estimation error reveal weak dependency on signal-to-noise ratio (SNR) with singular behavior of the estimation error in certain DOA angles.     

An approach for parameter estimation of combined CPPM and LFM radar signal

In this paper, the problem of parameter estimation of the combined radar signal adopting chaotic pulse position modulation (CPPM) and linear frequency modulation (LFM), which can be widely used in electronic countermeasures, is addressed. An approach is proposed to estimate the initial frequency and chirp rate of the combined signal by exploiting the second-order cyclostationarity of the intra-pulse signal. In addition, under the condition of the equal pulse width, the pulse repetition interval (PRI) of the combined signal is predicted using the low-order Volterra adaptive filter. Simulations demonstrate that the proposed cyclic autocorrelation Hough transform (CHT) algorithm is theoretically tolerant to additive white Gaussian noise. When the value of signal noise to ratio (SNR) is less than 4 dB, it can still estimate the intra-pulse parameters well. When SNR = 3 dB, a good prediction of the PRI sequence can be achieved by the Volterra adaptive filter algorithm, even only 100 training samples.