Efficient Beamspace Eigen-Based Direction of Arrival Estimation schemes

The Multiple SIgnal Classification (MUSIC) algorithm developed in the late 70\u27s was the first vector subspace approach used to accurately determine the arrival angles of signal wavefronts impinging upon an array of sensors. As facilitated by the geometry associated with the common uniform linear array of sensors, a root-based formulation was developed to replace the computationally intensive spectral search process and was found to offer an enhanced resolution capability in the presence of two closely-spaced signals. Operation in beamspace, where sectors of space are individually probed via a pre-processor operating on the sensor data, was found to offer both a performance benefit and a reduced computationa1 complexi ty resulting from the reduced data dimension associated with beamspace processing. Little progress, however, has been made in the development of a computationally efficient Root-MUSIC algorithm in a beamspace setting. Two approaches of efficiently arriving at a Root-MUSIC formulation in beamspace are developed and analyzed in this Thesis. In the first approach, a structura1 constraint is placed on the beamforming vectors that can be exploited to yield a reduced order polynomial whose roots provide information on the signal arrival angles. The second approach is considerably more general, and hence, applicable to any vector subspace angle estimation algorithm. In this approach, classical multirate digital signal processing is applied to effectively reduce the dimension of the vectors that span the signal subspace, leading to an efficient beamspace Root-MUSIC (or ESPRIT) algorithm. An auxiliaay, yet important, observation is shown to allow a real-valued eigenanalysis of the beamspace sample covariance matrix to provide a computational savings as well as a performance benefit, particularly in the case of correlated signal scenes. A rigorous theoretical analysis, based upon derived large-sample statistics of the signal subspace eigenvectors, is included to provide insight into the operation of the two algorithmic methodologies employing the real-valued processing enhancement. Numerous simulations are presented to validate the theoretical angle bias and variance expressions as well as to assess the merit of the two beamspace approaches

Sensor array signal processing : two decades later

Caption title.Includes bibliographical references (p. 55-65).Supported by Army Research Office. DAAL03-92-G-115 Supported by the Air Force Office of Scientific Research. F49620-92-J-2002 Supported by the National Science Foundation. MIP-9015281 Supported by the ONR. N00014-91-J-1967 Supported by the AFOSR. F49620-93-1-0102Hamid Krim, Mats Viberg

Sensor Array Signal Processing via Eigenanalysis of Matrix Pencils Composed of Data Derived from Translationally Invariant Subarrays

An algorithm is developed for estimating characteristic parameters associated with a scene of radiating sources given the data derived from a pair of translationally invariant arrays, the X and Y arrays, which are displaced relative to one another. The algorithm is referred to as PR O—E SPRIT and is predicated on invoking two recent mathematical developments: (1) the SVD based solution to the Procrustes problem of optimally approximating an invariant subspace rotation and (2) the Total Least Squares method for perturbing each of the two estimates of a common subspace in a minimal fashion until the two perturbed spaces are the same. For uniform linear array scenarios, the use of forward-backward averaging (FBAVG) in conjunction with PR O—E S PR IT is shown to effect a substantial reduction in the computational burden, a significant improvement in performance, a simple scheme for estimating the number of sources and source decorrelation. These gains may be attributed to FBAVG’s judicious exploitation of the diagonal invariance operator relating the Direction of Arrival matrix of the Y array to that associated with the X array. Similar gains may be achieved in the case where the X and Y arrays are either not linear or not uniformly spaced through the use of pseudo-forward-backward averaging (PFBAVG). However, the use of PFBAVG does not effect source decorrelation and reduces the maximum number of resolvable sources by a factor of two. Simulation studies and the results of applying PR O—E S PR IT to real data demonstrate the excellent performance of the method

Wavefield modeling and signal processing for sensor arrays of arbitrary geometry

Sensor arrays and related signal processing methods are key technologies in many areas of engineering including wireless communication systems, radar and sonar as well as in biomedical applications. Sensor arrays are a collection of sensors that are placed at distinct locations in order to sense physical phenomena or synthesize wavefields. Spatial processing from the multichannel output of the sensor array is a typical task. Such processing is useful in areas including wireless communications, radar, surveillance and indoor positioning. In this dissertation, fundamental theory and practical methods of wavefield modeling for radio-frequency array processing applications are developed. Also, computationally-efficient high-resolution and optimal signal processing methods for sensor arrays of arbitrary geometry are proposed. Methods for taking into account array nonidealities are introduced as well. Numerical results illustrating the performance of the proposed methods are given using real-world antenna arrays. Wavefield modeling and manifold separation for vector-fields such as completely polarized electromagnetic wavefields and polarization sensitive arrays are proposed. Wavefield modeling is used for writing the array output in terms of two independent parts, namely the sampling matrix depending on the employed array including nonidealities and the coefficient vector depending on the wavefield. The superexponentially decaying property of the sampling matrix for polarization sensitive arrays is established. Two estimators of the sampling matrix from calibration measurements are proposed and their statistical properties are established. The array processing methods developed in this dissertation concentrate on polarimetric beamforming as well as on high-resolution and optimal azimuth, elevation and polarization parameter estimation. The proposed methods take into account array nonidealities such as mutual coupling, cross-polarization effects and mounting platform reflections. Computationally-efficient solutions based on polynomial rooting techniques and fast Fourier transform are achieved without restricting the proposed methods to regular array geometries. A novel expression for the Cramér-Rao bound in array processing that is tight for real-world arrays with nonidealities in the asymptotic regime is also proposed. A relationship between spherical harmonics and 2-D Fourier basis, called equivalence matrix, is established. A novel fast spherical harmonic transform is proposed, and a one-to-one mapping between spherical harmonic and 2-D Fourier spectra is found. Improvements to the minimum number of samples on the sphere that are needed in order to avoid aliasing are also proposed

Signal Subspace Processing in the Beam Space of a True Time Delay Beamformer Bank

A number of techniques for Radio Frequency (RF) source location for wide bandwidth signals have been described that utilize coherent signal subspace processing, but often suffer from limitations such as the requirement for preliminary source location estimation, the need to apply the technique iteratively, computational expense or others. This dissertation examines a method that performs subspace processing of the data from a bank of true time delay beamformers. The spatial diversity of the beamformer bank alleviates the need for a preliminary estimate while simultaneously reducing the dimensionality of subsequent signal subspace processing resulting in computational efficiency. The pointing direction of the true time delay beams is independent of frequency, which results in a mapping from element space to beam space that is wide bandwidth in nature. This dissertation reviews previous methods, introduces the present method, presents simulation results that demonstrate the assertions, discusses an analysis of performance in relation to the Cramer-Rao Lower Bound (CRLB) with various levels of noise in the system, and discusses computational efficiency. One limitation of the method is that in practice it may be appropriate for systems that can tolerate a limited field of view. The application of Electronic Intelligence is one such application. This application is discussed as one that is appropriate for a method exhibiting high resolution of very wide bandwidth closely spaced sources and often does not require a wide field of view. In relation to system applications, this dissertation also discusses practical employment of the novel method in terms of antenna elements, arrays, platforms, engagement geometries, and other parameters. The true time delay beam space method is shown through modeling and simulation to be capable of resolving closely spaced very wideband sources over a relevant field of view in a single algorithmic pass, requiring no course preliminary estimation, and exhibiting low computational expense superior to many previous wideband coherent integration techniques

Signal Subspace Processing in the Beam Space of a True Time Delay Beamformer Bank

A number of techniques for Radio Frequency (RF) source location for wide bandwidth signals have been described that utilize coherent signal subspace processing, but often suffer from limitations such as the requirement for preliminary source location estimation, the need to apply the technique iteratively, computational expense or others. This dissertation examines a method that performs subspace processing of the data from a bank of true time delay beamformers. The spatial diversity of the beamformer bank alleviates the need for a preliminary estimate while simultaneously reducing the dimensionality of subsequent signal subspace processing resulting in computational efficiency. The pointing direction of the true time delay beams is independent of frequency, which results in a mapping from element space to beam space that is wide bandwidth in nature. This dissertation reviews previous methods, introduces the present method, presents simulation results that demonstrate the assertions, discusses an analysis of performance in relation to the Cramer-Rao Lower Bound (CRLB) with various levels of noise in the system, and discusses computational efficiency. One limitation of the method is that in practice it may be appropriate for systems that can tolerate a limited field of view. The application of Electronic Intelligence is one such application. This application is discussed as one that is appropriate for a method exhibiting high resolution of very wide bandwidth closely spaced sources and often does not require a wide field of view. In relation to system applications, this dissertation also discusses practical employment of the novel method in terms of antenna elements, arrays, platforms, engagement geometries, and other parameters. The true time delay beam space method is shown through modeling and simulation to be capable of resolving closely spaced very wideband sources over a relevant field of view in a single algorithmic pass, requiring no course preliminary estimation, and exhibiting low computational expense superior to many previous wideband coherent integration techniques

3-D Beamspace ML Based Bearing Estimator Incorporating Frequency Diversity and Interference Cancellation

The problem of low-angle radar tracking utilizing an array of antennas is considered. In the low-angle environment, echoes return from a low flying target via a specular path as well as a direct path. The problem is compounded by the fact that the two signals arrive within a beamwidth of each other and are usually fully correlated, or coherent. In addition, the SNR at each antenna element is typically low and only a small number of data samples, or snapshots, is available for processing due to the rapid movement of the target. Theoretical studies indicates that the Maximum Likelihood (ML) method is the only reliable estimation procedure in this type of scenario. However, the classical ML estimator involves a multi-dimensional search over a multi-modal surface and is consequently computationally burdensome. In order to facilitate real time processing, we here propose the idea of beamspace domain processing in which the element space snapshot vectors are first operated on by a reduced Butler matrix composed of three orthogonal beamforming weight vectors facilitating a simple, closed-form Beamspace Domain ML (BDML) estimator for the direct and specular path angles. The computational simplicity of the method arises from the fact that the respective beams associated with the three columns of the reduced Butler matrix have all but three nulls in common. The performance of the BDML estimator is enhanced by incorporating the estimation of the complex reflection coefficient and the bisector angle, respectively, for the symmetric and nonsymmetric multipath cases. To minimize the probability of track breaking, the use of frequency diversity is incorporated. The concept of coherent signal subspace processing is invoked as a means for retaining the computational simplicity of single frequency operation. With proper selection of the auxiliary frequencies, it is shown that perfect focusing may be achieved without iterating. In order to combat the effects of strong interfering sources, a novel scheme is presented for adaptively forming the three beams which retains the feature of common nulls

Hybrid Beamforming via the Kronecker Decomposition for the Millimeter-Wave Massive MIMO Systems

Despite its promising performance gain, the realization of mmWave massive MIMO still faces several practical challenges. In particular, implementing massive MIMO in the digital domain requires hundreds of RF chains matching the number of antennas. Furthermore, designing these components to operate at the mmWave frequencies is challenging and costly. These motivated the recent development of hybrid-beamforming where MIMO processing is divided for separate implementation in the analog and digital domains, called the analog and digital beamforming, respectively. Analog beamforming using a phase array introduces uni-modulus constraints on the beamforming coefficients, rendering the conventional MIMO techniques unsuitable and call for new designs. In this paper, we present a systematic design framework for hybrid beamforming for multi-cell multiuser massive MIMO systems over mmWave channels characterized by sparse propagation paths. The framework relies on the decomposition of analog beamforming vectors and path observation vectors into Kronecker products of factors being uni-modulus vectors. Exploiting properties of Kronecker mixed products, different factors of the analog beamformer are designed for either nulling interference paths or coherently combining data paths. Furthermore, a channel estimation scheme is designed for enabling the proposed hybrid beamforming. The scheme estimates the AoA of data and interference paths by analog beam scanning and data-path gains by analog beam steering. The performance of the channel estimation scheme is analyzed. In particular, the AoA spectrum resulting from beam scanning, which displays the magnitude distribution of paths over the AoA range, is derived in closed-form. It is shown that the inter-cell interference level diminishes inversely with the array size, the square root of pilot sequence length and the spatial separation between paths.Comment: Submitted to IEEE JSAC Special Issue on Millimeter Wave Communications for Future Mobile Networks, minor revisio

High resolution beamforming of SIMRAD EM3000 bathymetric multibeam sonar data

Originally, beamforming meant forming the transmit beam for some purpose. Today, we use the word beamforming both for the transmit – and the receive beam. The beam is simply the mainlobe of the array’s radiation – or response pattern. Changing this pattern, either by weighting the array elements or by steering its direction, is what beamforming is all about. In most signal processing applications, there is a tradeoff between different needs. The tradeoff in bathymetric sonars is between cross-track coverage, accuracy and resolution. The goal of this thesis, is to improve the resolution of a bathymetric sonar, working only with the output of the array – after the tradeoff is done. To achieve this, the Fourier transform method for beamforming and some high-resolution, also called adaptive, methods have been studied. In this thesis, the performance of several high-resolution beamformers are studied: Minimum-Variance, Eigen-Vector, MUSIC, Minimum-Norm, root-MUSIC and ESPRIT. The conventional Fourier transform, used for beamforming in many array systems today, suffers from low resolution. It has low ability of resolving closely spaced sources. For bathymetric sonars, this may lead to incorrect mapping of the sea-floor. Different ways of estimating the spatial covariance matrix are investigated, and a method aimed to reduce structured noise in this matrix is developed. Some quantitative measures are developed and used as part of the evaluation of the high-resolution beamformers. The measures show that estimation of the spatial covariance matrix has a large impact on the performance of the beamformers, and that most methods perform better when the structured noise in this matrix is reduced. An alternative approach for a high-resolution beamformer (root-MUSIC) is also developed, and is found to perform better than the one given in literature. It is found that among the high-resolution beamformers studied here, MUSIC and ESPRIT perform best with regard to the quantitative measures used for evaluation. These results are supported by visual examination of different spectra for the beamformers. Finally, some interesting research areas based on the observations in the thesis are proposed, and a new way of obtaining bathymetric measures is suggested