Fact-checked by Grok 2 weeks ago

Direction of arrival

Direction of arrival (DOA) estimation is a core technique in array signal processing used to determine the angular directions from which propagating waves, such as electromagnetic, acoustic, or seismic signals, impinge on a sensor array, typically by analyzing phase differences or time delays across the sensors.^[1] This process enables the localization of signal sources and is fundamental to applications requiring spatial awareness of incoming signals.^[2] The origins of DOA estimation trace back to World War II-era beamforming methods, exemplified by the conventional Bartlett beamformer, which applies Fourier analysis to the array output for directional scanning but suffers from limited resolution for closely spaced sources.^[1] Significant progress occurred in the late 20th century with the advent of subspace-based algorithms in the 1970s and 1980s, leveraging eigenvalue decomposition of the signal covariance matrix to separate signal and noise subspaces.^[2] A landmark development was the Multiple Signal Classification (MUSIC) algorithm, introduced by Schmidt in 1986, which provides high-resolution DOA estimates by identifying peaks in a pseudospectrum derived from the noise subspace, achieving asymptotic unbiasedness and consistency even in noisy environments.^[3] Building on this, the Estimation of Signal Parameters via Rotational Invariance Techniques (ESPRIT) algorithm, proposed by Roy and Kailath in 1989, exploits the translational invariance of subarrays to compute DOAs directly from the signal subspace eigenvectors, offering reduced computational complexity compared to search-based methods like MUSIC without sacrificing accuracy.^[4] Other notable classical approaches include Capon's minimum variance distortionless response (MVDR) beamformer for improved sidelobe suppression and maximum likelihood estimators for optimal performance under Gaussian noise assumptions, though they often require intensive optimization.^[1] Contemporary advancements address limitations of traditional methods, such as sensitivity to coherent signals, low signal-to-noise ratios (SNR), and scenarios with more sources than sensors, through techniques like spatial smoothing for decorrelation, compressed sensing for sparse representations, and deep learning-based estimators that learn complex array responses from data.^[5] These evolutions have expanded DOA estimation's utility in radar and sonar for target tracking, wireless communications for smart antenna beamforming and user localization, acoustic systems for sound source separation, seismology for earthquake monitoring, and emerging fields like 5G/6G networks and autonomous vehicles.^[2] Despite these strides, challenges persist in real-time implementation, computational efficiency for large arrays, and robustness in multipath or non-uniform environments.^[5]

Overview

Definition and Principles

Direction of arrival (DOA) estimation is a fundamental technique in array signal processing that determines the angle at which a signal wavefront impinges on an array of sensors, typically measured relative to a reference axis such as the array's broadside direction.^[6]^[7] This process leverages the spatial differences in signal reception across multiple sensors to infer the incoming direction, enabling applications in radar, sonar, wireless communications, and acoustic localization.^[7] The core principles of DOA estimation rely on several key assumptions to model signal propagation accurately. Under the far-field assumption, sources are sufficiently distant from the array such that the incoming signals can be approximated as plane waves, where wavefronts are planar and parallel across the sensor array.^[6]^[7] Additionally, the narrowband approximation is employed, positing that the signal's bandwidth is much smaller than the reciprocal of the maximum propagation delay across the array, which simplifies the phase relationships and allows frequency-independent processing.^[6]^[7] A central concept in DOA estimation is the steering vector, which captures the array's response to a plane wave arriving from direction \theta. For a uniform linear array (ULA) with M elements spaced d apart and operating at wavelength \lambda, the steering vector is given by

\mathbf{a}(\theta) = \begin{bmatrix} 1 \\ e^{-j 2\pi d \sin\theta / \lambda} \\ \vdots \\ e^{-j 2\pi (M-1) d \sin\theta / \lambda} \end{bmatrix}.

^[7] This vector arises from the phase differences induced by the signal's path length variations to each sensor, where the phase shift at the m-th element is proportional to (m-1) d \sin\theta / \lambda, forming the basis for angle estimation through correlation or spectral analysis of the received signals.^[6]^[7]

Historical Development

The origins of direction of arrival (DOA) estimation trace back to the late 19th century with early experiments in radio direction finding (RDF). In 1888, Heinrich Hertz discovered the directional properties of radio waves using a simple loop antenna, demonstrating that signal strength varied with orientation, which laid the groundwork for basic RDF systems.^[8] Practical advancements followed in the early 20th century, particularly with single-antenna techniques. By 1907, Italian engineers Ettore Bellini and Alessandro Tosi developed the Bellini-Tosi direction finder, which combined two perpendicular loop antennas with a goniometer to determine signal bearings more accurately without mechanical rotation, marking a significant improvement in maritime and aviation navigation.^[9] These analog methods relied on amplitude comparisons and were foundational for locating transmitters in wireless telegraphy. Developments during World War II in the 1940s shifted focus toward array-based systems for radar applications, driven by military needs for precise targeting. The MIT Radiation Laboratory, established in 1940, pioneered advancements in phased array radars, enabling electronic beam steering to estimate arrival directions without physical movement of antennas.^[10] Classical beamforming techniques emerged during this era as a core method for DOA estimation in radar, where arrays formed directional beams to scan for signals, improving resolution over single-antenna RDF.^[10] These innovations, tested in systems like the SCR-584 gun-laying radar, established array processing as essential for high-accuracy applications in defense and early sonar. The digital era transformed DOA estimation in the 1970s and 1980s through subspace-based methods, leveraging computational power for super-resolution beyond classical limits. A pivotal milestone was Ralph O. Schmidt's introduction of the MUSIC algorithm in 1986, which exploited the eigenstructure of the signal covariance matrix to resolve closely spaced sources with high precision.^[3] Building on this, Richard Roy and Thomas Kailath proposed the ESPRIT algorithm in 1989, offering a computationally efficient alternative by using rotational invariance in uniform linear arrays to estimate DOAs without spectral searching.^[11] These parametric techniques revolutionized array signal processing, finding widespread adoption in radar, sonar, and communications. Since the early 2000s, DOA estimation has increasingly integrated machine learning, particularly for challenging environments in 5G and beyond-5G networks. Neural network-based estimators emerged around 2015–2020, using deep learning to handle non-stationary signals, multipath propagation, and low SNR conditions, often outperforming traditional methods in massive MIMO systems.^[12] For instance, convolutional neural networks trained on array snapshots have demonstrated robust DOA resolution in dynamic scenarios, supporting beamforming in millimeter-wave 5G deployments.^[13] This trend reflects a broader shift toward data-driven approaches, enhancing adaptability for emerging wireless standards.

Fundamentals

Signal Model

In array signal processing, the narrowband signal model describes the reception of plane waves from multiple distant sources at a sensor array of M elements. The received signal vector at time t is given by

\mathbf{x}(t) = \mathbf{A}(\boldsymbol{\theta}) \mathbf{s}(t) + \mathbf{n}(t),

where \mathbf{x}(t) \in \mathbb{C}^{M \times 1} is the observation vector, \mathbf{A}(\boldsymbol{\theta}) = [\mathbf{a}(\theta_1), \dots, \mathbf{a}(\theta_K)] \in \mathbb{C}^{M \times K} is the array manifold matrix with columns being the steering vectors \mathbf{a}(\theta_k) for the directions-of-arrival (DOAs) \boldsymbol{\theta} = [\theta_1, \dots, \theta_K]^T, \mathbf{s}(t) \in \mathbb{C}^{K \times 1} contains the complex source signals, and \mathbf{n}(t) \in \mathbb{C}^{M \times 1} is the additive noise vector.^[14]^[15] The steering vector \mathbf{a}(\theta_k) captures the phase shifts due to the signal's propagation delays across the array elements for the k-th source, serving as the fundamental building block of the manifold matrix \mathbf{A}(\boldsymbol{\theta}).^[14] Under the narrowband assumption—where the signal bandwidth is much smaller than the inverse of the maximum propagation delay across the array—the second-order statistics are characterized by the covariance matrix

\mathbf{R}_x = E[\mathbf{x}(t) \mathbf{x}^H(t)] = \mathbf{A}(\boldsymbol{\theta}) \mathbf{R}_s \mathbf{A}^H(\boldsymbol{\theta}) + \mathbf{R}_n,

with \mathbf{R}_s = E[\mathbf{s}(t) \mathbf{s}^H(t)] as the source covariance matrix and \mathbf{R}_n = E[\mathbf{n}(t) \mathbf{n}^H(t)] as the noise covariance. For uncorrelated sources, \mathbf{R}_s is diagonal with distinct powers on the diagonal; the noise is typically modeled as spatially white Gaussian with \mathbf{R}_n = \sigma^2 \mathbf{I}, where \sigma^2 is the noise variance and \mathbf{I} is the identity matrix. These formulations assume K < M known sources, zero-mean signals uncorrelated with noise, and far-field plane-wave incidence.^[14]^[15] When sources are coherent—such as in multipath propagation, where \mathbf{R}_s becomes rank-deficient—the standard model fails as the signal subspace dimension reduces below K. This issue is addressed via spatial smoothing preprocessing, which constructs L = M - J + 1 overlapping subarrays of size J from the original array, averages their covariance matrices to restore full rank, and applies subsequent estimation on the smoothed covariance \tilde{\mathbf{R}}_x = \frac{1}{L} \sum_{l=1}^L \mathbf{R}_{x_l}, where \mathbf{R}_{x_l} is the covariance of the l-th subarray.^[16]

Array Configurations

In direction of arrival (DOA) estimation, the geometric arrangement of sensors in an array significantly influences the accuracy, resolution, and field of coverage for locating incoming signals. Common configurations are designed to balance simplicity, computational efficiency, and robustness to ambiguities such as grating lobes, which arise from spatial undersampling. These geometries adapt the steering vector to the array's layout, enabling tailored signal processing for one- or two-dimensional angle estimation.^[14] The uniform linear array (ULA) consists of sensors equally spaced along a straight line, typically with inter-element spacing of half the signal wavelength to avoid grating lobes. This configuration offers simplicity in implementation and processing, making it a foundational choice for one-dimensional DOA estimation, particularly for azimuth angles in narrow fields of view. However, ULAs suffer from ambiguities in endfire directions (near the array axis), where signals from opposite ends produce identical phase patterns, limiting unambiguous coverage to less than 180 degrees. Such ambiguities can be resolved using techniques like root-MUSIC, which roots a polynomial derived from the array's response to identify unique directions.^[17] The uniform circular array (UCA) arranges sensors equidistantly on a circular aperture, providing omnidirectional 360-degree coverage ideal for azimuth-only estimation without directional bias. Unlike the ULA, the UCA avoids endfire ambiguities due to its symmetric geometry, enabling reliable DOA performance across the full azimuthal range. To leverage efficient linear processing algorithms, the UCA's covariance matrix is often transformed into that of a virtual ULA through phase mode excitation or forward-backward averaging, facilitating application of methods like MUSIC or ESPRIT.^[18] Other notable geometries include the uniform rectangular array (URA), which extends the ULA concept to a planar grid for two-dimensional DOA estimation of both azimuth and elevation angles, suitable for applications like radar or sonar requiring joint angular resolution. URAs exploit separability in the x-y plane for reduced-complexity processing, though they demand larger apertures and more sensors than linear arrays. Sparse arrays, such as nested or coprime configurations, intentionally irregularize sensor positions to enlarge the effective aperture while minimizing redundancy, thereby increasing the degrees of freedom for resolving more sources than physical sensors. These designs suppress grating lobes by ensuring non-uniform sampling that fills holes in the spatial spectrum, enhancing resolution at the cost of increased design complexity.^[19]^[20]^[21] Key trade-offs among these configurations involve snapshot efficiency, where ULAs excel by requiring fewer temporal samples for covariance estimation due to their structured manifold; degrees of freedom, maximized in sparse arrays to support super-resolution beyond sensor count; and calibration requirements, which are minimal for uniform arrays like ULAs and UCAs but essential for non-uniform or sparse geometries to compensate for gain-phase mismatches and mutual coupling. Non-uniform arrays often necessitate precise sensor position knowledge and error modeling to maintain estimation accuracy, particularly in coherent signal environments.^[14]^[21]

Classical Estimation Methods

Beamforming Techniques

Beamforming techniques provide a foundational spatial filtering method for direction of arrival (DOA) estimation, where the array's response is steered toward hypothesized directions to maximize output power from potential sources while attenuating interference from other angles. DOA estimates are obtained by scanning the beam pattern and detecting peaks in the resulting spatial spectrum, which represents the array's power output as a function of angle θ. This approach treats DOA estimation as a spectrum search problem, relying on the array geometry and signal statistics without assuming specific source models beyond narrowband assumptions. The Bartlett beamformer, a classic implementation of conventional beamforming, originates from the delay-and-sum principle, which compensates for differential time delays across array elements to coherently align and sum signals from a target direction. For a uniform linear array of M sensors spaced by d, the steering vector \mathbf{a}(\theta) encodes the phase shifts \exp(-j 2\pi (m-1) d \sin\theta / \lambda) for m = 1, \dots, M, where \lambda is the signal wavelength. The beamformer output for a received snapshot \mathbf{x}(t) is y(\theta, t) = \frac{1}{M} \mathbf{a}^H(\theta) \mathbf{x}(t), and the power spectrum is P(\theta) = \mathbf{a}^H(\theta) \mathbf{R}_x \mathbf{a}(\theta), with \mathbf{R}_x = \mathbb{E}[\mathbf{x}(t) \mathbf{x}^H(t)] the covariance matrix derived from the array signal model. Peaks in P(\theta) correspond to the DOAs of impinging sources. The conventional discrete Fourier transform (DFT) beamformer discretizes this process by applying a DFT across the array elements to generate the beam pattern over a grid of angles, enabling efficient angular scanning via fast Fourier transform implementations. Its resolution is constrained by the array aperture L \approx M d, yielding an angular separation limit of \Delta\theta \approx \lambda / L radians between resolvable sources, beyond which peaks merge due to the beam's mainlobe width. These techniques offer low computational cost, typically O(M^2) per direction in direct computation or O(M \log M) with FFT acceleration for full scans, facilitating real-time processing on modest hardware. However, they exhibit poor resolution for closely spaced sources and elevated sidelobe levels that can mask weak signals with interference from nearby directions. To address sidelobe issues in practical implementations, windowing functions are applied to the steering vector or sensor weights, such as the Hamming window w(m) = 0.54 - 0.46 \cos(2\pi m / (M-1)), which suppresses sidelobes by 40–50 dB at the expense of a slightly broadened mainlobe and reduced directivity.^[22]

Capon Method

The Capon method, also known as the minimum variance distortionless response (MVDR) beamformer, is an adaptive spectral estimation technique that enhances direction-of-arrival (DOA) estimation by minimizing the array output power while preserving a unity gain in the direction of interest.^[23] Introduced by Jack Capon in 1969, it leverages the data covariance matrix to suppress interferers and noise, providing higher resolution than conventional methods.^[23] This approach is particularly effective for scenarios with correlated sources or closely spaced signals, as it adaptively forms nulls toward interference directions.^[14] The core principle of the Capon beamformer involves solving an optimization problem that minimizes the variance of the beamformer output subject to a distortionless constraint in the look direction \theta. The optimal weight vector is given by

\mathbf{w}(\theta) = \frac{\mathbf{R}_x^{-1} \mathbf{a}(\theta)}{\mathbf{a}^H(\theta) \mathbf{R}_x^{-1} \mathbf{a}(\theta)},

where \mathbf{R}_x is the spatial covariance matrix of the received signals, \mathbf{a}(\theta) is the steering vector for direction \theta, and ^H denotes the Hermitian transpose.^[14] The resulting power spectrum estimate, which peaks at the true DOA(s), is

P(\theta) = \frac{1}{\mathbf{a}^H(\theta) \mathbf{R}_x^{-1} \mathbf{a}(\theta)}.

This formulation weights the array elements inversely to the signal covariance, effectively reducing contributions from noise and interferers while maintaining the desired signal integrity.^[23] Computing the Capon spectrum requires inverting the M \times M covariance matrix \mathbf{R}_x, where M is the number of array sensors, leading to a computational complexity of O(M^3) per estimate, plus O(M^2 I) for evaluating over I angular grid points.^[14] In practice, when \mathbf{R}_x is ill-conditioned due to limited snapshots or strong correlations, diagonal loading is applied by adding a small positive constant \epsilon to the diagonal: \mathbf{R}_{x, \text{loaded}} = \mathbf{R}_x + \epsilon \mathbf{I}, which stabilizes the inversion and improves robustness without significantly degrading performance.^[14] Compared to the Bartlett (conventional) beamformer, the Capon method offers superior sidelobe suppression and DOA resolution, especially for correlated sources, by adaptively adjusting weights based on data statistics rather than using fixed uniform weighting.^[24] This results in narrower main lobes and deeper nulls toward interferers, though at the cost of higher computational demands.^[14]

Parametric Estimation Methods

MUSIC Algorithm

The Multiple Signal Classification (MUSIC) algorithm represents a seminal high-resolution subspace-based approach to direction-of-arrival (DOA) estimation, exploiting the eigendecomposition of the array covariance matrix to distinguish signal and noise subspaces.^[3] Introduced by Schmidt in 1986, it achieves super-resolution by leveraging the asymptotic orthogonality between the steering vectors of incident signals and the noise subspace, enabling the localization of multiple sources with angular separations finer than those resolvable by conventional beamforming methods.^[3]^[14] The algorithm proceeds in several key steps. First, the sample covariance matrix \mathbf{R}_x is computed from the array snapshots as \mathbf{R}_x = \frac{1}{N} \sum_{n=1}^N \mathbf{x}(n) \mathbf{x}^H(n), where N denotes the number of snapshots and \mathbf{x}(n) the M \times 1 observation vector at an M-element array.^[3] Next, the eigendecomposition \mathbf{R}_x = \mathbf{U} \boldsymbol{\Lambda} \mathbf{U}^H is performed, partitioning the eigenvectors into the signal subspace \mathbf{U}_s (spanned by the K eigenvectors corresponding to the largest eigenvalues) and the noise subspace \mathbf{U}_n (spanned by the remaining M - K eigenvectors associated with near-equal noise eigenvalues).^[3] The MUSIC pseudospectrum is then formed as

P_{\text{MUSIC}}(\theta) = \frac{1}{\mathbf{a}^H(\theta) \mathbf{U}_n \mathbf{U}_n^H \mathbf{a}(\theta)},

where \mathbf{a}(\theta) is the M \times 1 steering vector for potential DOA \theta.^[3] Peaks in P_{\text{MUSIC}}(\theta) are identified by scanning over \theta, with the K highest peaks yielding the DOA estimates \hat{\theta}_k.^[3] These peaks emerge because the true source steering vectors \mathbf{a}(\theta_k) lie in the signal subspace and are thus orthogonal to \mathbf{U}_n, causing the denominator to approach zero at the correct DOAs and producing sharp maxima in the spectrum.^[3] The resolution stems directly from this subspace orthogonality, which persists asymptotically as the number of snapshots increases, allowing MUSIC to resolve closely spaced sources even in low signal-to-noise ratio conditions, provided K < M.^[3]^[14] MUSIC requires knowledge of the source count K, typically estimated via information-theoretic criteria such as the Akaike Information Criterion (AIC) or Minimum Description Length (MDL), which minimize expressions balancing model fit and complexity based on the eigenvalue distribution of \mathbf{R}_x.^[25] It excels with uncorrelated sources, yielding consistent and unbiased estimates under white noise assumptions, but performance can degrade for coherent or correlated signals without modifications like spatial smoothing.^[3]^[14] For uniform linear arrays (ULAs), the Root-MUSIC variant circumvents the fine angular grid search of standard MUSIC by formulating a polynomial from the noise subspace eigenvectors and finding its roots. Specifically, the roots z_k nearest the unit circle in the complex plane correspond to the DOAs via \theta_k = \sin^{-1} \left( -\frac{\lambda \angle z_k}{2\pi d} \right), where \angle z_k is the phase angle of the root, \lambda is the signal wavelength, and d is the inter-element spacing, offering computational efficiency for one-dimensional searches.^[26]

ESPRIT Algorithm

The ESPRIT (Estimation of Signal Parameters via Rotational Invariance Techniques) algorithm is a subspace-based method for direction-of-arrival (DOA) estimation that leverages the geometric structure of sensor arrays to achieve high-resolution performance without exhaustive spectral searches. Introduced as an efficient alternative to methods like MUSIC, ESPRIT exploits the translational invariance in the array configuration by dividing the sensors into two overlapping subarrays that are identical in structure but shifted by a fixed inter-element spacing d. This invariance induces a rotational relationship in the signal subspace, enabling direct computation of DOA parameters from the phase rotations between subarrays.^[4] The core principle relies on the eigendecomposition of the array covariance matrix \mathbf{R}_x = \mathbf{E}\{\mathbf{x}(t) \mathbf{x}^H(t)\}, which yields the signal subspace \mathbf{U}_s spanning the range of the steering matrix. For arrays with translational invariance, such as uniform linear arrays (ULAs), the signal subspace induces subarray subspaces \mathbf{U}_{s1} = \mathbf{J}_1 \mathbf{U}_s and \mathbf{U}_{s2} = \mathbf{J}_2 \mathbf{U}_s, where \mathbf{J}_1 and \mathbf{J}_2 are selection matrices that select the sensors corresponding to the first and second subarrays, respectively. These satisfy the approximate relation

\mathbf{U}_{s2} \approx \mathbf{U}_{s1} \mathbf{\Psi},

where \mathbf{\Psi} is a nonsingular diagonal matrix whose diagonal elements are the phase shifts \psi_k = e^{-j 2\pi d \sin \theta_k / \lambda} for each of the K sources, with \lambda denoting the signal wavelength and \theta_k the DOA of the k-th source. This equation captures how the subarray displacement introduces predictable phase rotations in the signal subspace, distinct from the noise subspace.^[4] DOA estimation proceeds by estimating \mathbf{\Psi} via a least-squares approach, yielding

\mathbf{\Psi} = (\mathbf{U}_{s1}^H \mathbf{U}_{s1})^{-1} \mathbf{U}_{s1}^H \mathbf{U}_{s2},

where the eigenvalues of \mathbf{\Psi} provide the \psi_k. The angles are then recovered using

\theta_k = \arcsin\left( -\frac{\lambda \angle \psi_k}{2\pi d} \right).

This closed-form solution avoids the one-dimensional search required in spectral methods, making ESPRIT computationally efficient with complexity on the order of O(M^3 + K^3), where M is the number of sensors and K the number of sources—dominated by the eigendecomposition step for typical scenarios where K \ll M. ESPRIT performs robustly at low signal-to-noise ratios and closely approaches the Cramér-Rao bound for uncorrelated sources under ideal conditions.^[4] A key advantage of ESPRIT is its applicability to ULAs, where the regular spacing naturally provides the required shift invariance, facilitating real-time implementation in resource-constrained systems like radar or wireless communications. However, the algorithm's reliance on translational invariance imposes limitations, rendering it less flexible for arbitrary or nonuniform array geometries, such as circular or sparse arrays, where the subarray selection may not preserve the exact rotational structure. Extensions like unitary ESPRIT address some numerical stability issues but retain this geometric constraint.^[4]

Advanced Topics

Wideband DOA Estimation

Wideband direction-of-arrival (DOA) estimation addresses the limitations of narrowband methods when signals span a significant bandwidth relative to the center frequency, where the steering vector varies with frequency, leading to smearing of the signal subspace and reduced estimation resolution, especially as the time-bandwidth product increases. This frequency dependence causes the narrowband far-field approximation to break down, distorting the orthogonality between signal and noise subspaces in classical algorithms like MUSIC.^[27] Incoherent methods mitigate these issues by decomposing the wideband signal into narrowband frequency bins via Fourier transform or filter banks, applying a narrowband DOA estimator (such as MUSIC) to each bin to obtain per-frequency estimates, and then averaging them across bins to yield a final DOA.^[28] To account for varying signal-to-noise ratios (SNRs) across frequencies, weighted averaging schemes can be employed, where weights are derived from the subspace eigenvalues or SNR estimates to emphasize reliable bins and suppress noisy ones, improving overall accuracy at moderate SNRs. The coherent signal subspace method (CSSM) offers a more integrated approach by transforming the frequency-dependent signal subspaces from all bins into a common reference frequency using focusing matrices, thereby preserving subspace coherence and enabling high-resolution subspace-based estimation on the aligned data. These focusing matrices, typically designed via Taylor series approximation of the steering vector phase or least-squares optimization, align the wideband manifold to the reference steering vectors, after which standard subspace methods like MUSIC are applied to the focused covariance matrix.^[27] CSSM performs robustly at low SNRs by exploiting the full bandwidth information without averaging losses, though it requires accurate initial DOA guesses for matrix design to avoid focusing errors. Advanced true wideband variants extend parametric methods like ESPRIT to broadband scenarios by employing focusing matrices to align signal subspaces across frequencies prior to rotational invariance exploitation, avoiding the need for per-bin processing or averaging.^[27] These approaches, such as focused wideband ESPRIT, construct a broadband signal model where the focusing transformation ensures consistent shift-invariance properties, yielding computationally efficient closed-form DOA estimates even for coherent sources when combined with spatial smoothing.^[29] Recent advancements incorporate deep learning techniques for wideband DOA estimation, leveraging neural networks to learn complex frequency-dependent array responses directly from data. For instance, ConvNeXt-based methods enhance accuracy for wideband signals by processing raw array data through convolutional layers, achieving superior performance in low-SNR and coherent scenarios as of 2024. Similarly, support vector regression frameworks have been adapted for efficient wideband DOA, showing improved resolution in multidisciplinary applications.^[30]^[31]

Robustness to Impairments

Real-world direction-of-arrival (DOA) estimation often encounters impairments such as mutual coupling, sensor gain and phase errors, and non-Gaussian noise, which degrade the accuracy of classical and parametric methods by distorting the array response or covariance structure. Robust techniques address these by incorporating models of imperfections into the estimation framework or using statistical methods that downweight anomalous data. These approaches enhance reliability in practical applications like radar and acoustics, where array imperfections are inevitable. Mutual coupling arises from electromagnetic interactions between closely spaced array elements, altering the effective steering vector and introducing biases in DOA estimates. This effect is typically modeled using the array's impedance matrix, which captures the voltage-current relationships across elements, allowing the coupling matrix to be incorporated into the signal model as a transformation on the ideal response.^[32] Mitigation strategies include calibration methods that estimate the coupling parameters from calibration signals or self-calibration using the received data alone, such as time-frequency distribution-based approaches that jointly estimate DOAs and coupling coefficients without prior knowledge.^[33] Physical decoupling networks, like parasitic elements or metamaterial-based isolators, can also reduce coupling at the hardware level by modifying the current distribution to minimize inter-element interactions. Sensor gain and phase errors represent manufacturing or environmental mismatches that perturb the array manifold, often modeled as stochastic multiplicative factors with Gaussian-distributed uncertainties to reflect variability. Robust estimators, such as weighted subspace fitting (WSF), mitigate these by formulating the estimation as a weighted least-squares problem over the signal subspace, where weights inversely proportional to error variances reduce sensitivity to perturbations. This approach outperforms standard subspace methods in mismatched scenarios, achieving near-optimal performance when errors are small, as demonstrated in analyses of uniform linear arrays under gain-phase uncertainties.^[34] Non-Gaussian noise, characterized by heavy-tailed distributions or sparse outliers, violates the Gaussian assumptions of conventional covariance estimators, leading to degraded subspace decomposition. Tyler's M-estimator provides a robust alternative by estimating the scatter matrix through an iterative fixed-point equation that scales the sample covariance to unit determinant, offering consistency for elliptical distributions without assuming a specific noise pdf. In DOA contexts, integrating Tyler's estimator into methods like MUSIC or ESPRIT enhances outlier resistance, as shown in formulations for complex elliptically symmetric noise where it maintains estimation accuracy under impulsive interference.^[35] Compressive sensing complements this by enabling sparse recovery of DOAs amid outliers, modeling noise as sparse corruptions and using optimization to separate signal and anomaly components. Compressive sensing approaches leverage the spatial sparsity of signal sources, assuming few active directions relative to the grid, to formulate DOA estimation as a sparse recovery problem resilient to impairments. For off-grid sources, where true DOAs fall between discrete grid points causing basis mismatch, \ell_1-norm minimization solves \min_{\mathbf{s}} \|\mathbf{s}\|_1 subject to \|\mathbf{y} - \mathbf{A}(\boldsymbol{\theta}) \mathbf{s}\|_2 \leq \epsilon, with \mathbf{A}(\boldsymbol{\theta}) as the dictionary of steering vectors; this promotes sparsity while the grid can be refined iteratively. Such methods, often combined with reweighted schemes, improve resolution and robustness in underdetermined scenarios with coherent sources or low snapshots.^[36] Emerging robust methods utilize deep learning to handle impairments more adaptively. For example, CRDCNN-LSTM frameworks jointly optimize channel attention and recurrent processing for DOA estimation under mutual coupling and noise, demonstrating enhanced robustness in 2025 simulations. Additionally, hyperdimensional computing-based estimators like HYPERDOA provide efficient, lightweight solutions resilient to sensor errors and non-Gaussian noise, suitable for real-time applications as of October 2025.^[37]^[38]

Applications

Radar and Sonar Systems

In radar systems, monopulse techniques provide precise angle tracking by estimating the direction of arrival (DOA) of target echoes through phase or amplitude comparisons across multiple antenna elements, enabling rapid and accurate target localization without mechanical scanning.^[39] This method has been integral to surveillance radars since the mid-20th century, offering high angular resolution for tracking fast-moving targets in real-time applications.^[40] Phased array radars, particularly active electronically scanned arrays (AESAs), leverage DOA estimation for electronic beam steering, allowing adaptive formation of beams toward detected targets to enhance signal-to-noise ratios and support multi-function operations like search and track.^[41] AESAs have been deployed in military platforms since the 1990s, exemplified by systems like the AN/APG-77 on the F-22 Raptor, where DOA-derived steering improves jamming resistance and simultaneous target engagement.^[42] In sonar systems, passive sonar employs hydrophone arrays to estimate the bearing of underwater targets, such as submarines, by applying beamforming to DOA measurements derived from time differences of arrival across array elements.^[43] This approach is critical for stealthy detection in noisy oceanic environments, where conventional beamforming scans the bearing space to identify acoustic signatures without emitting signals that could reveal the platform's position.^[44] For instance, towed array sonars on naval vessels use DOA estimation to achieve high-resolution bearing accuracy, often integrating adaptive beamforming to suppress ambient noise and focus on target directions.^[45] DOA measurements from radar and sonar are integrated into multi-target tracking frameworks, such as Kalman filters, to fuse angular data with range and velocity estimates for robust state prediction amid occlusions or maneuvers.^[46] In the AN/SPY-1 radar, a phased-array system on Aegis-equipped destroyers, DOA-based tracking supports simultaneous monitoring of hundreds of air and surface threats, incorporating Kalman filtering to refine trajectories and assign threats to interceptors.^[47] This integration enhances overall system performance by reducing false alarms and improving prediction accuracy in dynamic scenarios.^[48] Key challenges in these applications include clutter in radar, where ground or sea returns degrade DOA accuracy, necessitating adaptive filtering to distinguish targets from environmental echoes.^[49] In sonar, multipath propagation from surface or bottom reflections causes ambiguous DOA estimates, particularly in shallow waters, requiring robust algorithms like sparse Bayesian learning to resolve true bearings.^[50] Addressing these impairments is essential for maintaining operational reliability in contested environments.

Acoustic Source Localization

Microphone arrays enable direction of arrival (DOA) estimation in acoustic source localization by exploiting phase differences across sensors to determine the incoming sound direction, facilitating applications in audio processing such as speech enhancement and source tracking. Beamforming techniques, which steer spatial filters toward estimated DOA, are widely used to suppress noise and reverberation while amplifying the target signal. For instance, delay-and-sum beamforming constructs a weighted sum of microphone signals aligned by time delays derived from DOA, improving signal-to-noise ratio (SNR) through array gain in noisy environments.^[51] In teleconferencing systems, microphone array beamforming enhances speech clarity by focusing on the speaker's direction, reducing interference from background noise or other participants. The Microsoft Kinect, introduced in 2010, integrates a four-microphone linear array for real-time DOA estimation and beamforming, enabling robust speech capture in interactive scenarios like gaming and video calls, where it achieves localization accuracy within 5 degrees in quiet rooms.^[52] This approach has influenced subsequent devices by combining DOA with acoustic echo cancellation for hands-free communication. Reverberant environments pose challenges due to multipath propagation, which distorts DOA estimates by creating multiple signal arrivals. Subspace methods, such as the MUSIC algorithm adapted for acoustics, mitigate this by decomposing the signal covariance matrix into signal and noise subspaces, identifying DOA peaks in the noise subspace projection even under moderate reverberation (reverberation time up to 0.5 seconds). To further handle multipath, time-difference-of-arrival (TDOA) estimation integrates with generalized cross-correlation phase transform (GCC-PHAT), a robust method that normalizes the cross-spectrum to emphasize phase differences. These techniques often combine subspace DOA with TDOA for hybrid localization, improving resolution in non-ideal acoustics. Smart assistants like the Amazon Echo employ real-time DOA estimation using circular microphone arrays (e.g., seven microphones) to detect voice activity and direct beamforming toward the user, enhancing wake-word detection and command recognition in far-field settings up to 3 meters.^[53] Similarly, Google Home (now Nest) devices use multi-microphone configurations for DOA-based voice activity detection, steering beams to isolate speech from household noise and enabling multi-user interaction. These systems process wideband audio signals, adapting narrowband DOA methods to broadband speech spectra for low-latency performance. In biomedical applications, DOA principles are adapted to electroencephalography (EEG) and magnetoencephalography (MEG) arrays for brain source localization, treating sensors as virtual microphone arrays to estimate neural activity origins. Beamforming techniques, such as linearly constrained minimum variance (LCMV), spatially filter EEG/MEG data to localize dipolar sources with high precision, suppressing interference from unrelated brain regions.^[54] Subspace extensions like recursive MUSIC further enable multi-source localization in MEG, resolving correlated neural signals during tasks like auditory processing.^[55] These adaptations highlight DOA's versatility beyond acoustics, aiding in epilepsy mapping and cognitive neuroscience studies.

Other Applications

In wireless communications, DOA estimation supports smart antenna systems for beamforming, directing signals to specific users to enhance spectral efficiency and reduce interference. This is particularly vital in 5G and 6G networks, where massive multiple-input multiple-output (MIMO) configurations use DOA to enable user localization and adaptive beam management in dense environments.^[56] Seismology employs DOA estimation with sensor arrays to determine the direction of incoming seismic waves, aiding in earthquake epicenter localization and monitoring tectonic activity. Techniques like beamforming on tri-axial geophones provide bearing information for event characterization in real-time seismic networks.^[57] In autonomous vehicles, DOA estimation is applied in radar and acoustic systems for obstacle detection and localization. Automotive radars use high-resolution DOA methods to resolve closely spaced targets, supporting safe navigation and collision avoidance in dynamic traffic scenarios.^[58]

Performance Evaluation

Resolution Limits

The resolution of direction-of-arrival (DOA) estimation refers to the minimum angular separation between two sources that can be reliably distinguished by an array processing algorithm. In conventional beamforming techniques, such as the Bartlett or Capon beamformer, the fundamental limit is imposed by the Rayleigh criterion, which arises from the diffraction-like behavior of the array's spatial response. For a uniform linear array (ULA) with M elements spaced by d, the Rayleigh resolution is approximately \Delta \theta = 0.89 \frac{\lambda}{M d \cos \theta} radians, where \lambda is the signal wavelength and \theta is the angle of incidence from broadside; this limit scales inversely with the array aperture M d and broadens at off-broadside angles due to the \cos \theta projection factor.^[59] Subspace-based methods, such as MUSIC and ESPRIT, enable super-resolution by exploiting the eigenstructure of the array covariance matrix, allowing distinction of sources separated by angles much smaller than the Rayleigh limit—potentially on the order of \Delta \theta \approx \frac{1}{M \sqrt{\mathrm{SNR}}} in high signal-to-noise ratio (SNR) asymptotic regimes for large arrays.^[60] Several practical factors influence the achievable resolution beyond these theoretical expressions. The number of snapshots (temporal samples) used to estimate the covariance matrix is critical, as insufficient snapshots lead to poor eigenvalue separation and degraded resolution, particularly for closely spaced sources. Source correlation, such as from multipath propagation, reduces effective degrees of freedom in the signal subspace, causing resolution failure unless decorrelation techniques like spatial smoothing are applied. The array aperture directly determines baseline resolution, with larger apertures providing finer discrimination, though practical constraints like mutual coupling can counteract this benefit. Asymptotic analysis for large M reveals that resolution improves proportionally with array size, but performance saturates or degrades if snapshots or SNR are inadequate relative to M. Simulations demonstrate threshold effects where resolution abruptly fails at low SNR; for instance, in a ULA with M=8 elements and two uncorrelated sources separated by 5°, MUSIC successfully resolves them above 0 dB SNR but exhibits breakdown—manifesting as merged peaks in the spatial spectrum—below -5 dB SNR with 100 snapshots, highlighting the SNR-dependent threshold.^[61] Across methods, MUSIC achieves near-optimal resolution for uncorrelated sources by accurately identifying the noise subspace, outperforming beamforming by factors of 5–10 in angular separation under moderate SNR, though it requires more snapshots than root-MUSIC variants for correlated scenarios.^[60]

Cramer-Rao Bound

The Cramér-Rao lower bound (CRB) represents the theoretical minimum variance achievable by any unbiased estimator of the direction-of-arrival (DOA) parameters in array signal processing. It is derived from the inverse of the Fisher information matrix \mathbf{I}(\theta), which measures the sensitivity of the likelihood function to changes in the unknown DOA vector \theta. In the standard narrowband DOA model, the array observations are \mathbf{x}(t) = \mathbf{A}(\theta)\mathbf{s}(t) + \mathbf{n}(t) for t = 1, \dots, N snapshots, where \mathbf{A}(\theta) is the steering matrix, \mathbf{s}(t) are the source signals, and \mathbf{n}(t) \sim \mathcal{CN}(\mathbf{0}, \mathbf{R}_n) is additive Gaussian noise. The CRB states that for an unbiased estimator \hat{\theta}, \mathrm{var}(\hat{\theta}_k) \geq [\mathbf{I}^{-1}(\theta)]_{kk} for the k-th DOA component.^[62] For the general case where both the mean \boldsymbol{\mu} = \mathbf{A}\mathbf{s} and covariance \mathbf{R}_x depend on \theta, the Fisher information matrix under the complex Gaussian assumption is

\mathbf{I}(\theta) = 2 \Re\left\{ \frac{\partial \boldsymbol{\mu}^H}{\partial \theta} \mathbf{R}_n^{-1} \frac{\partial \boldsymbol{\mu}}{\partial \theta} \right\} + \mathrm{tr}\left( \mathbf{R}_n^{-1} \frac{\partial \mathbf{R}_x}{\partial \theta} \mathbf{R}_n^{-1} \frac{\partial \mathbf{R}_x}{\partial \theta} \right),

with the factor of N scaling for multiple snapshots. In the deterministic signal model (treating \mathbf{s}(t) as unknown constants), \mathbf{R}_x = \mathbf{R}_n and the trace term vanishes, reducing to the mean-contribution term scaled by N, while concentrating out the nuisance signal parameters yields \mathbf{I}(\theta) = 2N \Re \left\{ \left( \frac{\partial \mathbf{A}^H}{\partial \theta_k} \right) (\mathbf{I} - \mathbf{A}(\mathbf{A}^H \mathbf{A})^{-1} \mathbf{A}^H) \mathbf{R}_n^{-1} \frac{\partial \mathbf{A}}{\partial \theta_k} \right\} P_k for the k-th source power P_k under white noise \mathbf{R}_n = \sigma^2 \mathbf{I}. In the stochastic model ( \mathbf{s}(t) zero-mean Gaussian), the mean term is zero, and \mathbf{R}_x = \mathbf{A} \mathbf{R}_s \mathbf{A}^H + \mathbf{R}_n, leading to the trace term N \mathrm{tr} \left( \mathbf{R}_x^{-1} \frac{\partial \mathbf{R}_x}{\partial \theta_i} \mathbf{R}_x^{-1} \frac{\partial \mathbf{R}_x}{\partial \theta_j} \right). The resulting CRB covariance matrix is \mathrm{CRB}(\theta) = \mathbf{I}^{-1}(\theta).^[62] Maximum likelihood (ML) estimators for DOA are asymptotically unbiased and efficient, achieving the CRB as the number of snapshots N \to \infty, provided the Fisher information matrix is positive definite (requiring distinct DOAs and sufficient sensors). This efficiency holds under both deterministic and stochastic models, with the deterministic CRB typically tighter at low SNR or few snapshots, while the stochastic CRB better reflects random signal fluctuations. Extensions of the CRB account for practical complications. For correlated sources, where \mathbf{R}_s has off-diagonal structure or reduced rank (e.g., coherent signals), the Fisher information decreases due to signal subspace ambiguity, inflating the CRB; a closed-form adjustment involves the eigendecomposition of \mathbf{R}_s. When noise covariance \mathbf{R}_n is unknown, nuisance parameter estimation loosens the bound, often requiring structured assumptions like spatial whiteness. For a uniform linear array (ULA) with M sensors spaced d = \lambda/2, a closed-form CRB for a single uncorrelated source at angle \theta (from broadside) and white noise is

\mathrm{var}(\hat{\theta}) \geq \frac{6 \sigma^2}{N P \left( \frac{2\pi d}{\lambda} \right)^2 \cos^2 \theta \cdot M (M^2 - 1)},

highlighting the quadratic scaling with array aperture.^[62] The CRB interpretation reveals key performance drivers: variance inversely scales with signal-to-noise ratio (SNR = P / \sigma^2), snapshot count N, and source power P, while array geometry influences tightness via the steering vector derivatives—e.g., larger apertures or non-linear geometries (like UCAs) reduce the bound by enhancing angular sensitivity, though grating lobes can degrade it at endfire angles. At low SNR, the bound emphasizes the need for coherent processing over snapshots.^[62]

References

[1]
Direction of Arrival - an overview | ScienceDirect Topics
The direction of arrival of waves can be simply determined by the relative time delay of signals with similar shapes recorded by a group of adjacent sensors.<|control11|><|separator|>
[2]
An Overview of Direction-of-Arrival Estimation Methods Using ...
Apr 21, 2022 · The direction of arrival (DOA) is one of the critical topics in array signal processing. Initially, the DOA was estimated for wireless signals ...
[3]
https://ieeexplore.ieee.org/document/1143830
[4]
ESPRIT-estimation of signal parameters via rotational invariance ...
ESPRIT-estimation of signal parameters via rotational invariance techniques | IEEE Journals & Magazine | IEEE Xplore ... Date of Publication: 31 July 1989.Missing: pdf | Show results with:pdf
[5]
https://ieeexplore.ieee.org/document/10522654
[6]
[PDF] Direction of Arrival Estimation 1 Introduction
We have seen that there is a one-to-one relationship between the direction of a signal and the associated received steering vector.
[7]
https://arxiv.org/pdf/2508.11675.pdf
[8]
Radio Direction Finder - Engineering and Technology History Wiki
Nov 23, 2017 · In 1910, in Boulogne, France, Bellini and Tosi develop one of the earliest successful radio direction finding antennas.
[9]
[PDF] radio direction finding
the loop. In 1907, the Bellini-Tosi antenna sys- tem, named after its inventors, was developed. The principle of the Bellini-Tosi is used in most. Page 7. (+).
[10]
[PDF] The development of phased-array radar technology
Early phased-array research was centered on improving the national capability in phased-array radars. The Laboratory has developed several test-bed phased ...
[11]
[PDF] ESPRIT-estimation of signal parameters via rotational invariance ...
Such applica- tions include direction-of-arrival (DOA) estimation, system identifica- tion, and time series analysis. A novel approach to the general problem of ...
[12]
DOA Estimation in B5G/6G: Trends and Challenges - PMC - NIH
Jul 8, 2022 · This article discusses various open issues related to DOA estimation for B5G/6G communication networks.
[13]
Direction of Arrival Estimation: A Tutorial Survey of Classical ... - arXiv
Aug 8, 2025 · Direction of arrival (DOA) estimation is the process of determining the spatial directions from which signals impinge on an array of sensors ...
[14]
None
### Summary of DOA Estimation from arXiv:2508.11675
[15]
Array Signal Processing | SpringerLink
Johnson, D.H. and D.E. Dudgeon, Array Signal Processing: Concepts and Techniques, Prentice-Hall: Englewood Cliffs, NJ, 1993. MATH Google Scholar. Kanasewich ...
[16]
https://ieeexplore.ieee.org/document/1164649
[17]
Direction of Arrival Estimation Using Sparse Arrays - MathWorks
This example constructs several popular sparse array architectures and shows how they can be used to perform DOA estimation.Minimum Redundancy Array · Nested Array · Coprime ArraysMissing: review rectangular<|control11|><|separator|>
[18]
Transformation of a circular array to a linear ... - ScienceDirect.com
We address the transformation of a uniform circular array to a virtual uniform linear array. This phase mode transformation is widely used.
[19]
https://www.worldscientific.com/doi/10.1142/S021812661950049X
[20]
[PDF] Direction-of-Arrival Estimation for Uniform Rectangular Array
As one typical kind of planar arrays, URA could be deployed to estimate 2-D DOA, i.e., elevation and azimuth [1],. [2]. Several DOA estimation methods, e.g., ...
[21]
[PDF] Sparse Arrays and Sampling for Interference Mitigation and DOA ...
These arrays have many advantages, including reduced redundancy, larger physical and virtual apertures, and avoidance of grating lobes [15]–[17]. Direction of ...
[22]
A Comparison of Beamforming and Direction Finding Algorithms ...
Usually, a Hamming or an ultraspherical window (wj) is applied before beamforming (see Gurgel et al. 1999b; Helzel and Kniephoff 2010) to increase suppression ...
[23]
High-resolution frequency-wavenumber spectrum analysis
High-resolution frequency-wavenumber spectrum analysis | IEEE Journals ... Date of Publication: 31 August 1969. ISSN Information: Print ISSN: 0018-9219.
[24]
[PDF] Doubly Constrained Robust Capon Beamformer
The Capon beamformer has better resolution and much better interference rejection capability than the data-independent beamformer, provided that the array ...
[25]
https://ieeexplore.ieee.org/document/1164557
[26]
Focussing matrices for coherent signal-subspace processing
Aug 31, 1988 · Abstract: A class of focussing matrices is proposed for use in the coherent signal-subspace method (CSM) (H. Wang and M. Kaveh, ibid., vol.
[27]
An Overview of Different Wideband Direction of Arrival(DOA ...
The Direction of Arrival (DOA) estimation methods are useful in Sonar, Radar and Advanced Satellite and Cellular Communication systems. In this paper ...Missing: survey | Show results with:survey
[28]
TOPS: new DOA estimator for wideband signals - ResearchGate
Aug 6, 2025 · This paper introduces a new direction-of-arrival (DOA) estimation algorithm for wideband sources called test of orthogonality of projected subspaces (TOPS).Missing: seminal | Show results with:seminal
[29]
https://www.researchgate.net/publication/3319689_TOPS_new_DOA_estimator_for_wideband_signals
[30]
DOA Estimation and Self-Calibration under Unknown Mutual Coupling
Feb 25, 2019 · This paper proposes a new self-calibration method based on the time-frequency distributions (TFDs) in the presence of direction-dependent MC.Missing: mitigation seminal
[31]
Direction of arrival estimation for sparse arrays with gain-phase ...
In this work, we address the issue of DOA estimation for sparse arrays under the coexistence of gain-phase errors and nonuniform noise.
[32]
Robust and sparse M-estimation of DOA - ScienceDirect.com
A robust and sparse Direction of Arrival (DOA) estimator is derived for array data that follows a Complex Elliptically Symmetric (CES) distribution.
[33]
Grid-less T.V minimization for DOA estimation - ScienceDirect
In this paper we propose a method to alleviate this difficulty. This technique allows us to formulate the sparse recovery methods like L1-SVD or SPICE in a ...Missing: influential | Show results with:influential
[34]
Monopulse DOA estimation using adjacent matched filter samples
In modern surveillance radar systems, amplitude comparison monopulse techniques improve target angle estimation. By comparing the difference and sum of ...
[35]
Monopulse DOA Estimation for High-Speed Targets in Frequency ...
Apr 14, 2025 · The modified monopulse method facilitates rapid DOA estimation at low beam pointing angles, while the maximum likelihood estimation (MLE) method ...
[36]
The Development of Phased-Array Radar Technology - ResearchGate
The T/R module significantly enhances the EIRP of AESA radar, enables electronic steering and beam shaping/broadening, ensures linearity, optimizes the system's ...
[37]
[PDF] Active Electronically Scanned Arrays (AESA) Architectures
Jul 24, 2025 · AESAs are prevalent now in many different applications such as Defense, Communications, Automotive, etc. • Currently digital AESAs are being ...
[38]
Underwater Source Localization via Spectral Element Acoustic Field ...
DOA is determined using the time differences in the sound arrival among the spatially separated hydrophones of the array. This is the main difference from ...
[39]
Performance Analysis for Focused Beamformers in Passive ...
Apr 23, 2018 · ABSTRACT Focused beamformers are widely used in passive underwater acoustic localization. Many pseudo-peaks (spurious peaks not ...
[40]
Direction-of-Arrival Estimation of Far-Field Sources Under Near ...
Feb 23, 2021 · The paper proposes a near-field interference mitigation (NFIM) beamformer using subarrays and an azimuth-domain filter to improve DOA ...
[41]
[PDF] TARGET TRACKING: IMPLEMENTING THE KALMAN FILTER
In this study of linear filtering, the Kalman filter, a recursive linear filtering model, was used to estimate tracks. Various situations were examined, ...
[42]
AN/SPY-1 Radar - Missile Defense Advocacy Alliance
It operates in S-band and is a multi-function phased-array radar capable of search, automatic detection, transition to track, tracking of air and surface ...
[43]
[PDF] MULTI-TARGET, MULTI-MODEL TRACKING - Dr. Gregory L. Plett
It may be helpful to restart the Kalman filter upon reacquisition. Track crossings. □ Sometimes two targets will cross one another. □ When this happens, the ...
[44]
Overview of radar detection methods for low altitude targets in ...
In general, clutter and multipath scattering affect the background and target signals of the received radar sig- nal, respectively. The key challenge in clutter ...
[45]
DOA estimation based on multi-frequency joint sparse Bayesian ...
This paper develops a multi-frequency DOA estimation algorithm using sparse Bayesian learning, exploiting target sparsity, to improve accuracy and resolution.
[46]
https://depts.drew.edu/govschl/njgss2009/journal/TeamPapers/team5.pdf
[47]
https://missiledefenseadvocacy.org/defense-systems/anspy-1-radar/
[48]
Alexa? Techniques to Pull Clear Speech from Noise – Part II
Aug 25, 2021 · Direction of arrival detection (DOA) This algorithm compares slight delays in arrival time for a detected signal at different microphones.
[49]
A unified view on beamformers for M/EEG source reconstruction
Feb 1, 2022 · Beamforming is a popular method for functional source reconstruction using magnetoencephalography (MEG) and electroencephalography (EEG) data.
[50]
https://ieeexplore.ieee.org/iel7/5971804/9940202/09940163.pdf
[51]
Single-snapshot DOA estimation by using Compressed Sensing
Jul 30, 2014 · This paper deals with the problem of estimating the directions of arrival (DOA) of multiple source signals from a single observation vector of an array data.
[52]
https://www.research.unipd.it/retrieve/e14fb26f-e122-3de1-e053-1705fe0ac030/riccardo_levorato_tesi.pdf
[53]
[PDF] Performance analysis of an improved MUSIC DoA estimator - arXiv
Abstract—This paper adresses the statistical performance of subspace DoA estimation using a sensor array, in the asymptotic regime where the number of samples ...
[54]
https://www.sciencedirect.com/science/article/pii/S1053811921010612