Open AccessArticle

Efficient Two-Dimensional Direction Finding Algorithm for Rectilinear Sources Under Unknown Mutual Coupling

Electronics and Information School, Northwestern Polytechnical University, Xi’an 710072, China

Author to whom correspondence should be addressed.

Sensors 2020, 20(7), 1914; https://doi.org/10.3390/s20071914

Submission received: 18 January 2020 / Revised: 19 March 2020 / Accepted: 20 March 2020 / Published: 30 March 2020

(This article belongs to the Special Issue Recent Advances in Sensor Array and Multichannel Signal Processing with Its Applications to IoT Security)

Download

Browse Figures

Figure 1
The structure of a uniform rectangular array (URA) with M elements. "> Figure 2
The spatial spectra of (a) Two-dimensional multiple signal classification (2D-MUSIC), (b) the auxiliary sensor-based method (AUX), (c) the 2D rank reduction-based method (2D-RARE) and (d) the first stage and (e) the second stage of the proposed method. "> Figure 2 Cont.
The spatial spectra of (a) Two-dimensional multiple signal classification (2D-MUSIC), (b) the auxiliary sensor-based method (AUX), (c) the 2D rank reduction-based method (2D-RARE) and (d) the first stage and (e) the second stage of the proposed method. "> Figure 3
RMSEs of Root mean square errors (RMSEs) of direction-of-arrival (DOA) estimates versus signal-to-noise ratio (SNR). The number of snapshots is 1000. "> Figure 4
RMSEs of DOA estimates versus the number of snapshots. The SNR is 10 dB. "> Figure 5
Computational complexity of the four methods versus number of antennas. ">

Versions Notes

Abstract

Digital communication signals in wireless systems may possess noncircularity, which can be used to enhance the degrees of freedom for direction-of-arrival (DOA) estimation in sensor array signal processing. On the other hand, the electromagnetic characteristics between sensors in uniform rectangular arrays (URAs), such as mutual coupling, may significantly deteriorate the estimation performance. To deal with this problem, a robust real-valued estimator for rectilinear sources was developed to alleviate unknown mutual coupling in URAs. An augmented covariance matrix was built up by extracting the real and imaginary parts of observations containing the circularity and noncircularity of signals. Then, the actual steering vector considering mutual coupling was reparameterized to make the rank reduction (RARE) property available. To reduce the computational complexity of two-dimensional (2D) spectral search, we individually estimated y-axis and x-axis direction-cosines in two stages following the principle of RARE. Finally, azimuth and elevation angle estimates were determined from the corresponding direction-cosines respectively. Compared with existing solutions, the proposed method is more computationally efficient, involving real-valued operations and decoupled 2D spectral searches into twice those of one-dimensional searches. Simulation results verified that the proposed method provides satisfactory estimation performance that is robust to unknown mutual coupling and close to the counterparts based on 2D spectral searches, but at the cost of much fewer calculations.

Keywords:

DOA estimation; mutual coupling; rectilinear sources; RARE.

1. Introduction

Uniform rectangular array (URA), due to its versatility in providing both azimuthal and elevational coverage simultaneously, has attracted much attention in many applications, such as radar, sonar, and wireless communication [1]. Two-dimensional (2D) direction finding with URA is thus of practical use in many situations. Based on the subspace principle, the classical 1D super resolution algorithms such as multiple signal classification (MUSIC) [2] and compressive sensing [3] can be readily generalized to URA. Moreover, to further reduce the computational complexity or increase the accuracy of the direction finding results, some effective schemes have been proposed [4,5,6,7,8]. In [4], the authors proposed a 2D unitary ESPRIT algorithm for URA, which can be realized in both element space and beamspace. Also, its main procedure is formulated in terms of real-valued operation and is therefore computationally efficient. In [5], a dual-size spatial invariance array structure was proposed. Based on the ESPRIT theory, it utilizes unambiguous but low-accuracy estimates to disambiguate high-accuracy but periodically ambiguous direction cosine estimates. In [6], a 2D matrix pencil method is developed to convert the complex data matrix to a real one, which can reduce the complexity of the computation significantly. In [7], to relieve the computational burden involved in 2D spectral searches, the authors decoupled the 2D MUSIC problem into two successive real-valued one-dimensional (1D) DOA estimation problems. In [8], a 2D unitary ESPRIT-like DOA estimation algorithm was developed for coherent signals based on a URA. It can provide better decorrelation performance and closed-form solutions due to the properties of the 2D ESPRIT-like method.

Recently, in modern wireless communication systems, rectilinear or strictly noncircular signals such as amplitude modulated (AM), binary phase shift keying (BPSK), and pulse-amplitude modulated (PAM) signals are commonly applied. Notice that both conjugated and unconjugated covariance matrices of rectilinear signals include the second-order statistical information [9]. Hence, the rectilinearity can be utilized to promote the estimation accuracy and increase the maximum resolvable sources. Accordingly, various algorithms have been developed to fully exploit the rectilinearity for DOA estimation performance enhancement. In [10], Abeida and Delmas proposed a tailored MUSIC algorithm for noncircular sources called NCMUSIC, which combines the statistics of conjugated and unconjugated covariance matrices and thus outperforms the standard MUSIC algorithm. In [11], the sparse representation technique was introduced for representing conjugated and unconjugated covariance jointly with a sparsity constraint, which surpasses the performance of NCMUSIC. However, it is computationally demanding since it requires complicate convex optimization procedures. In [12], the difference and sum coarray was utilized for the DOA estimation of noncircular quasi-stationary sources, which can extend the degrees of freedom (DOF) of the array effectively. In [13,14], the scenario of coexisting circular and noncircular sources is considered, and various algorithms have been further developed to distinguish the mixed sources and estimate the DOA parameters. Based on an URA, the authors proposed a root-MUSIC method to obtain the 2D DOA parameters and classification of the mixed sources [15]. In [16], the authors proposed a localization algorithm for rectilinear sources based on the rank reduction (RARE) principle, which can cope with the mixed near-field and far-field sources.

Generally, the above mentioned DOA estimation algorithms are based on the assumption that the array manifold is ideal without any perturbation. However, in practical direction finding systems, systematic errors such as mutual coupling, gain and phase errors, and sensor position errors are inevitable [17]. These model mismatches would substantially degrade the DOA estimation performance of these algorithms. Mutual coupling effect, as an electromagnetic interaction between the array elements, causes the current distribution of each sensor to depend on the excitations from both their own and other adjacent element. It deviates the signal subspace of the observed array covariance matrix from the ideal array steering vector, thus causing a decrease in direction finding performance. Therefore, the mutual coupling effect cannot be ignored in the development of these algorithms.

Fortunately, many effective schemes have been proposed by various researchers to compensate for the mutual coupling effect or to obtain the DOA estimates blindly. One class of compensation methods are based on calibration sources. Using a set of predefined sources in known locations, the mutual coupling coefficients can be calculated by the least square method [18]. However, the calibration sources are often inaccessible in real-time direction finding applications. Moreover, the mutual coupling characteristics of the array may vary as the environmental conditions change, and therefore the calibration source-based methods are only effective under certain circumstances. Consequently, another class of auto-calibration algorithms is advanced to compensate for the mutual coupling effect without any calibration sources [19,20,21,22,23,24,25]. In [19], an eigenstructure-based algorithm was developed for the DOA estimation and mutual coupling auto-calibration at the same time. However, it requires multidimensional alternating minimization which is computationally prohibitive. To avoid exhaustive multidimensional searches over the unknown DOA and mutual coupling parameters, the authors proposed a modified MUSIC algorithm [20]. By taking the sensors at either side of the array as auxiliary ones, the authors found that the MUSIC algorithm is resilient against unknown mutual coupling effects. However, the array aperture loss problem of these methods is severe particularly under the strong mutual coupling scenario, because only the observations from middle subarray array elements can be utilized for DOA estimation. Aiming at resolving this problem, Liao et al. utilized the banded symmetric Toeplitz property of the mutual coupling matrix and transformed this matrix into gain and phase perturbations of the array [21], where the perturbations of the middle subarray elements are identical. Consequently, the observation data from all the array elements can be fully utilized for DOA estimation via simple 1D RARE method. Furthermore, a fourth order cumulant (FOC) method was also proposed [22], which performs better especially for strong mutual coupling effects. In [23,24,25], the authors utilized sparse reconstruction methods to improve the DOA estimation performance under unknown mutual coupling. However, these methods are computationally demanding due to the l₁-norm minimization problem. In [26], the authors presented a method for the localization of mixed sources using the noncircular information of the impinging signals under unknown mutual coupling. In [27], an ESPRIT-based algorithm is proposed for the ODA estimation of strictly noncircular sources with unknown mutual coupling, which has large virtual array aperture and low computational load. According to the latest work by B. Friedlander [28], the mutual coupling matrix is a direction dependent matrix, and not a constant matrix as is commonly assumed, which may impact some of the work presented in the DOA estimation literature. However, few of the works on DOA estimation have taken into account the direction-dependent mutual coupling effect. In [29], a decoupled 2D DOA estimation algorithm is proposed for uniform circular arrays. Yet, it assumes that the mutual coupling effect is only elevation-dependent. Also, it is an offline technique, for which is costly to build a calibration setup. In [30], a unified transformation method is proposed for the DOA estimation under direction-dependent mutual coupling. This method can be used in many array structures, but it requires convex optimization and iteration processes, which is computationally inefficient. In [31], the authors proposed a RARE-based low complexity DOA estimation method to solve the direction-dependent mutual coupling problem. Moreover, coprime sparse arrays are proposed to reduce the mutual coupling effect between the elements [32,33,34,35].

For the 2D array calibration problem, several effective methods have been developed [36,37,38,39,40]. In [36], the authors proposed a blind calibration scheme for uniform circular array based on azimuthal symmetric structures in the mutual coupling matrix. However, it only provides the azimuth angle estimates of the radiating sources, while the elevation angles are assumed to be zeros. In [37], a robust method was proposed to calibrate the mutual coupling and channel gain/phase inconsistency of uniform circular array. Yet, it requires calibration sources and 2D spectral searches. In [38], the auxiliary sensor-based algorithm is generalized for URA and uniform hexagon arrays, where the elements at the boundary of the array are set as the auxiliary ones. In order to mitigate the array aperture loss problem in the auxiliary sensor-based method, Wu et al. proposes a RARE-based auto-calibration method for URA [40], which can fully utilize the array observations. However, all these planar array calibration methods require exhaustive 2D spectral searches for the DOA estimation, which is computationally demanding.

To the extent of our knowledge, few works have been devoted to the 2D direction finding algorithms for rectilinear sources in the presence of unknown mutual coupling effects. Despite the various advantages that were manifested in the approaches cited above, there are some fundamentally inherent problems to be overcome in view of the previous analyses: (1) all of the above-mentioned 2D DOA estimation works only use the covariance matrix for direction finding and mutual coupling calibration, the rectilinearity of the sources are not fully utilized. (2) These works require 2D spectral searches for the azimuth and elevation angle estimations, which is computationally demanding. (3) Generally, subspace-based direction finding algorithms require complex-valued operations for the calculation of covariance matrix, eigenvalue decomposition (EVD), and spectral searching. Many works have been devoted to converting the complex-valued covariance into a real-valued one. The computational complexity of EVD and spectral searches for this kind of transformation can be decreased by a factor of four. Yet, the transformation itself and the computation of the covariance still require additional complex-valued operations.

Therefore, aiming at addressing these problems, we developed a computationally efficient 2D direction finding algorithm for rectilinear signals under unknown mutual coupling. Based on the rectilinearity of the signal, the received data is transformed into the real-valued domain. Then the structure of the virtual steering vector is utilized to decouple the DOA parameter from other nuisance parameters. According to the RARE principle, the y-axis direction-cosine estimates are first obtained via a 1D spectral search. Consequently, the x-axis direction-cosines are then generated through another RARE estimator, which also only requires a 1D spectral search with the y-axis direction-cosine estimates.

The main contributions of this paper are as follows:

(1): The proposed DOA estimation method avoids the two-dimensional spectral search and effectively reduces the computational complexity;
(2): Only real-valued operations are needed, and its computation load is reduced to 1/4 of the complex value operations;
(3): The estimation accuracy is improved by the efficient use of the source rectilinearity.

The rest of this paper is organized as follows. In Section 2, we present the signal model for the rectilinear sources of a URA under unknown mutual coupling. In Section 3, derivation of the proposed real-valued two stage 2D DOA estimation algorithm is detailed. Then we discuss the computational complexity and the maximum number of resolvable sources in Section 4. In Section 5, several numerical simulations are carried out to validate the performance of the proposed algorithm. In Section 6, the paper ends with conclusions. Finally, in the Appendix A, we have derived the Cramér-Rao Bound (CRB) for the problem.

Throughout the paper, we make use of the following notations as listed in Table 1.

2. Signal Model

As shown in Figure 1, consider a URA composed of M = M_x × M_y omnidirectional elements, with the inter element spacing along the x-axis and the y-axis being d_x and d_y, respectively. Suppose that there are K far-field narrowband sources radiating onto the URA from directions

(θ_{1}, ϕ_{1})

, …,

(θ_{K}, ϕ_{K})

, where

θ_{k}

and

ϕ_{k}

are the azimuth and elevation angles of the k-th source.

Under the ideal scenario, the received data of each element only depends on the incident wavefield. Accordingly, the ideal array received data vector without mutual coupling can be modeled as

x (t) = \sum_{k = 1}^{K} a (θ_{k}, ϕ_{k}) s_{k} (t) + n (t) = A s (t) + n (t)

(1)

where

x (t) = {[x_{1} (t), \dots, x_{M_{x}} (t), x_{M_{x} + 1} (t), \dots, x_{2 M_{x}} (t), \dots, x_{M_{y} M_{x}} (t)]}^{T}

s (t) = {[s_{1} (t), s_{2} (t), \dots, s_{K} (t)]}^{T}

is the source signal vector,

n (t) = {[n_{1} (t), \dots, n_{M_{x}} (t), n_{M_{x} + 1} (t), \dots, n_{2 M_{x}} (t), \dots, n_{M_{y} M_{x}} (t)]}^{T}

denotes the array noise vector, and

A = [a (θ_{1}, ϕ_{1}), \dots, a (θ_{K}, ϕ_{K})]

represents the ideal steering matrix, with its steering vectors being

a (θ_{k}, ϕ_{k}) = a_{y} (ψ_{y, k}) \otimes a_{x} (ψ_{x, k}), k = 1, \dots, K

(2)

where

a_{x} (ψ_{x, k}) = {[1, e^{j ψ_{x, k}}, \dots, e^{j (M_{x} - 1) ψ_{x, k}}]}^{T}

(3)

a_{y} (ψ_{y, k}) = {[1, e^{j ψ_{y, k}}, \dots, e^{j (M_{y} - 1) ψ_{y, k}}]}^{T}

(4)

Herein,

ψ_{x, k}

and

ψ_{y, k}

are the electrical angles of the kth source with respect to the x-axis and the y-axis, with

ψ_{x, k} = 2 π d_{x} \cos θ_{k} \sin ϕ_{k} / λ

and

ψ_{y, k} = 2 π d_{y} \sin θ_{k} \sin ϕ_{k} / λ

where λ is the carrier wavelength.

However, in practical direction finding systems, a mismatch may be introduced by the electromagnetic interaction between the array elements. Therefore, taking into account the mutual coupling effect between the array elements, we can modify the received data model as

x (t) = C A s (t) + n (t)

(5)

where C is the

M_{x} M_{y} \times M_{x} M_{y}

mutual coupling matrix, with its elements being mutual coupling coefficients (MCCs). As illustrated in [19], the mutual coupling effect between the elements is inversely proportional to their distance. Thus, the MCCs between the neighboring sensors are approximately equal. When two elements are far apart, the mutual coupling effect is negligible, and the coefficients approach zero. As in [30], the general mutual coupling model for URA is developed by the following block matrix

C = [\begin{matrix} C_{1}^{} & C_{2}^{} & \dots & C_{P}^{} \\ C_{2}^{} & C_{1}^{} & C_{2}^{} & ⋱ \\ ⋮ & ⋱ & ⋱ & ⋱ \\ C_{P}^{} & C_{2}^{} & C_{1}^{} & C_{2}^{} & C_{P}^{} \\ ⋱ & ⋱ & ⋱ & ⋱ & ⋮ \\ ⋱ & C_{2}^{} & C_{1}^{} & C_{2}^{} \\ C_{P}^{} & \dots & C_{2}^{} & C_{1}^{} \end{matrix}]

(6)

and the

M_{x} \times M_{x}

submatrix C_p (p = 1,…,P) is expressed as [27]

C_{p} = [\begin{matrix} c_{p - 1, 0} & c_{p - 1, 1} & \dots & c_{p - 1, P - 1} \\ c_{p - 1, 1} & c_{p - 1, 0} & c_{p - 1, 1} & ⋱ \\ ⋮ & c_{p - 1, 1} & c_{p - 1, 0} & ⋱ & c_{p - 1, P - 1} \\ c_{p - 1, P - 1} & ⋱ & ⋱ & c_{p - 1, 1} & ⋮ \\ ⋱ & c_{p - 1, 1} & c_{p - 1, 0} & c_{p - 1, 1} \\ c_{p - 1, P - 1} & \dots & c_{p - 1, P - 1} & c_{p - 1, 0} \end{matrix}]

(7)

The complex mutual coupling coefficient

c_{m, n} = ρ_{m, n} e^{j ξ_{m, n}}

symbolizes the electromagnetic interaction between a reference sensor defined as the origin of a Cartesian coordinate system and the sensor located at (±md_x, ±nd_y) under the corresponding coordinate [30]. Note that

c_{0, 0}

is especially defined as 1. Therefore, the mutual coupling matrix C can be further expressed as

C = \sum_{p = 1}^{P} J_{p} \otimes C_{p}

(8)

and

J_{p}

is a

M_{y} \times M_{y}

selection matrix, whose elements are

{[J_{p}]}_{m, n} = {\begin{array}{l} 1, | m - n | = p - 1 \\ 0, otherwise \end{array}

(9)

Moreover, for rectilinear sources, the signal vector is of the following form [10]

s (t) = Φ s_{0}^{} (t)

(10)

where

s_{0}^{} (t) = {[s_{0, 1}^{} (t), \dots, s_{0, K}^{} (t)]}^{T}

, with

s_{0, k} (t)

being the real-valued zero-phase version of the k-th source signal,

Φ = diag {e_{}^{j φ_{1}}, \dots, e_{}^{j φ_{K}}}

contains the noncircular phase shifts of the sources [9].

Owing to the above analysis, the received data vector of the planar array can be formulated as

x (t) = C A Φ s_{0} (t) + n (t)

(11)

With the received data vector x(t), the main problem to be addressed in this paper is to obtain the 2D DOA estimates of rectilinear signals under unknown mutual coupling effects. Unlike the traditional methods which require multidimensional spectral searches and omit the rectilinearity of the signals, a novel algorithm is investigated in the next section for reducing the computational complexity and promoting the estimation performance.

3. Proposed Real-Valued 2D DOA Estimation Algorithm

In order to generate accurate 2D DOA estimates with low computational complexity, we creatively utilize the rectilinearity of signals and factorize the steering vector with respect to the DOA and other nuisance parameters. Our technique proceeds through the following stages: (1) based on the rectilinearity of the signal, the received data is transformed into the real-valued domain; (2) the structure of the virtual steering vector is utilized to decouple the DOA parameter from other nuisance parameters; (3) according to the RARE principle, the y-axis direction-cosine estimates can be obtained via a 1D spectral search; (4) with the y-axis direction-cosine estimates, the x-axis direction-cosine can be found through another RARE estimator, which also only requires a 1D spectral search. The following subsections give the details of the proposed algorithm.

3.1. Extended Real-Valued Signal Model

It is well known that by transforming complex-valued data into the real-valued domain, the computational cost can be effectively reduced by a factor of four. Therefore, it inspired us to look for a more efficient algorithm to reduce the complex-valued computational loads.

Based on the sources’ rectilinearity, the transformation can be realized by extracting the real and imaginary parts of x(t), respectively.

x_{R} (t) = Re [x (t)], x_{I} (t) = Im [x (t)]

(12)

From Equation (11),

x_{R} (t)

and

x_{I} (t)

can be further expressed as

x_{R} (t) = A_{R} s_{0} (t) + n_{R} (t)

(13)

x_{I} (t) = A_{I} s_{0} (t) + n_{I} (t)

(14)

where

A_{R} = [v_{R} (θ_{1}, ϕ_{1}, φ_{1}, C), \dots, v_{R} (θ_{K}, ϕ_{K}, φ_{K}, C)]

and

A_{I} = [v_{I} (θ_{1}, ϕ_{1}, φ_{1}, C), \dots, v_{I} (θ_{K}, ϕ_{K}, φ_{K}, C)]

are the real-valued manifold matrices of

x_{R} (t)

and

x_{I} (t)

n_{R} (t)

and

n_{I} (t)

symbolize the real and imaginary components of sensor noise vector

n (t)

, respectively. Accordingly, the real-valued virtual steering vectors in

A_{R}

is of the following form

v_{R} (θ_{k}, ϕ_{k}, φ_{k}, C) = [C a (θ_{k}, ϕ_{k}) \cdot e^{j φ_{k}} + C^{*} a^{*} (θ_{k}, ϕ_{k}) \cdot e^{- j φ_{k}}] / 2

(15)

Based on Equations (2) and (8), we can obtain the following equation from the Kronecker product property

(A \otimes B) (C \otimes D) = (A C) \otimes (B D)

[41]

C a (θ_{k}, ϕ_{k}) = (\sum_{p = 1}^{P} J_{p} \otimes C_{p}) (a_{y} (ψ_{y, k}) \otimes a_{x} (ψ_{x, k})) = \sum_{p = 1}^{P} (J_{p} a_{y} (ψ_{y, k})) \otimes (C_{p} a_{x} (ψ_{x, k}))

(16)

Therefore, the real-valued virtual steering vector

v_{R} (θ_{k}, ϕ_{k}, φ_{k}, C)

can be further written as

\begin{matrix} v_{R} (θ_{k}, ϕ_{k}, φ_{k}, C) = \sum_{p = 1}^{P} (J_{p} a_{y, c} (ψ_{y, k})) \otimes (Γ_{c, p} (ψ_{x, k}) a_{x, c} (ψ_{x, k})) - \sum_{p = 1}^{P} (J_{p} a_{y, s} (ψ_{y, k})) \otimes (Γ_{c, p} (ψ_{x, k}) a_{x, s} (ψ_{x, k})) \\ - \sum_{p = 1}^{P} (J_{p} a_{y, c} (ψ_{y, k})) \otimes (Γ_{s, p} (ψ_{x, k}) a_{x, s} (ψ_{x, k})) - \sum_{p = 1}^{P} (J_{p} a_{y, s} (ψ_{y, k})) \otimes (Γ_{s, p} (ψ_{x, k}) a_{x, c} (ψ_{x, k})) \end{matrix}

(17)

where

a_{x, c} (ψ_{x, k}) = {[1, \cos (ψ_{x, k}), \dots, \cos ((M_{x} - 1) ψ_{x, k})]}^{T}

(18)

a_{x, s} (ψ_{x, k}) = {[0, \sin (ψ_{x, k}), \dots, \sin ((M_{x} - 1) ψ_{x, k})]}^{T}

(19)

a_{y, c} (ψ_{y, k}) = {[1, \cos (ψ_{y, k}), \dots, \cos ((M_{y} - 1) ψ_{y, k})]}^{T}

(20)

a_{y, s} (ψ_{y, k}) = {[0, \sin (ψ_{y, k}), \dots, \sin ((M_{y} - 1) ψ_{y, k})]}^{T}

(21)

Γ_{c, p} (ψ_{x, k}) = h_{c}^{(p)} \cdot diag [α_{c, 1}^{(p)}, \dots, α_{c, P - 1}^{(p)}, 1, \dots, 1, β_{c, 1}^{(p)}, \dots, β_{c, P - 1}^{(p)}]

(22)

Γ_{s, p} (ψ_{x, k}) = h_{s}^{(p)} \cdot diag [α_{s, 1}^{(p)}, \dots, α_{s, P - 1}^{(p)}, 1, \dots, 1, β_{s, 1}^{(p)}, \dots, β_{s, P - 1}^{(p)}]

(23)

h_{c}^{(p)} = \sum_{i = 1 - P}^{P - 1} ρ_{p - 1, | i |} \cos (ξ_{p - 1, | i |} + φ_{k} + i \cdot 2 π λ^{- 1} ψ_{x, k})

(24)

h_{s}^{(p)} = \sum_{i = 1 - P}^{P - 1} ρ_{p - 1, | i |} \sin (ξ_{p - 1, | i |} + φ_{k} + i \cdot 2 π λ^{- 1} ψ_{x, k})

(25)

Note that

Γ_{c, p} (ψ_{x, k})

and

Γ_{s, p} (ψ_{x, k})

are

M_{x} \times M_{x}

diagonal matrices comprising

M_{x} - 2 P + 2

ones between

α_{c, P - 1}^{(p)}

and

β_{c, 1}^{(p)}

(also between

α_{s, P - 1}^{(p)}

and

β_{s, 1}^{(p)}

h_{c}^{(p)}

and

h_{s}^{(p)}

are two scalar parameters determined by the x-axis direction-cosine, rectilinearity phases, and MCCs. For

p = 1, 2, \dots, P - 1

, we have

α_{c, q}^{(p)} = \frac{h_{c}^{(p)} - \sum_{i = q}^{P - 1} ρ_{p - 1, i} \cos (ξ_{p - 1, i} + φ_{k} - i \cdot 2 π λ^{- 1} ψ_{x, k})}{h_{c}^{(p)}}, q = 1, \dots, P - 1

(26)

α_{s, q}^{(p)} = \frac{h_{s}^{(p)} - \sum_{i = q}^{P - 1} ρ_{p - 1, i} \sin (ξ_{p - 1, i} + φ_{k} - i \cdot 2 π λ^{- 1} ψ_{x, k})}{h_{s}^{(p)}}, q = 1, \dots, P - 1

(27)

β_{c, q}^{(p)} = \frac{h_{c}^{(p)} - \sum_{i = P - q}^{P - 1} ρ_{p - 1, i} \cos (ξ_{p - 1, i} + φ_{k} + i \cdot 2 π λ^{- 1} ψ_{x, k})}{h_{c}^{(p)}}, q = 1, \dots, P - 1

(28)

β_{s, q}^{(p)} = \frac{h_{s}^{(p)} - \sum_{i = P - q}^{P - 1} ρ_{p - 1, i} \sin (ξ_{p - 1, i} + φ_{k} + i \cdot 2 π λ^{- 1} ψ_{x, k})}{h_{s}^{(p)}}, q = 1, \dots, P - 1

(29)

Similarly, the real-valued virtual steering vectors of

A_{I}

can be defined as

\begin{matrix} v_{I} (θ_{k}, ϕ_{k}, φ_{k}, C) = \sum_{p = 1}^{P} (J_{p} a_{y, c} (ψ_{y, k})) \otimes (Γ_{s, p} (ψ_{x, k}) a_{x, c} (ψ_{x, k})) - \sum_{p = 1}^{P} (J_{p} a_{y, s} (ψ_{y, k})) \otimes (Γ_{s, p} (ψ_{x, k}) a_{x, s} (ψ_{x, k})) \\ + \sum_{p = 1}^{P} (J_{p} a_{y, c} (ψ_{y, k})) \otimes (Γ_{c, p} (ψ_{x, k}) a_{x, s} (ψ_{x, k})) + \sum_{p = 1}^{P} (J_{p} a_{y, s} (ψ_{y, k})) \otimes (Γ_{c, p} (ψ_{x, k}) a_{x, c} (ψ_{x, k})) \end{matrix}

(30)

Based on the above real-valued signal model, we can concatenate the real part and the imaginary part to generate an augmented data vector in the real-valued domain.

y (t) = [\begin{matrix} x_{R}^{} (t) \\ x_{I}^{} (t) \end{matrix}] = [\begin{matrix} A_{R}^{} \\ A_{I}^{} \end{matrix}] s_{0}^{} (t) + [\begin{matrix} n_{R}^{} (t) \\ n_{I}^{} (t) \end{matrix}] = A_{E}^{} s_{0}^{} (t) + n_{E}^{} (t)

(31)

where

A_{E}^{} = [a_{E}^{} (θ_{1}^{}, ϕ_{1}^{}, φ_{1}^{}, C), \dots, a_{E}^{} (θ_{K}^{}, ϕ_{K}^{}, φ_{K}^{}, C)]

is the augmented manifold matrix, and the corresponding real-valued steering vector of the k-th source is defined as

a_{E}^{} (θ_{k}^{}, ϕ_{k}^{}, φ_{k}^{}, C) = [\begin{matrix} v_{R}^{} (θ_{k}^{}, ϕ_{k}^{}, φ_{k}^{}, C) \\ v_{I}^{} (θ_{k}^{}, ϕ_{k}^{}, φ_{k}^{}, C) \end{matrix}]

(32)

Therefore, the covariance matrix of the augmented real-valued data matrix y(t) is of the following form

R_{y} = E {y (t) y^{T} (t)} = A_{E} R_{S_{0}} A_{E}^{T} + σ_{n}^{2} I_{2 M}

(33)

where

R_{S_{0}}^{} = E {s_{0}^{} (t) s_{0}^{T} (t)}

I_{2 M}

denotes a 2M × 2M identity matrix. In practical applications,

R_{y}

can be approximated from the snapshots generated at L distinct time instants

{\hat{R}}_{y} = \frac{1}{L} \sum_{t = 1}^{L} {y (t) y^{T} (t)}

(34)

Through the eigenvalue decomposition of

R_{y}^{}

, we can obtain that

R_{y}^{} = E_{s}^{} Λ_{s}^{} E_{s}^{H} + σ_{n}^{2} E_{n}^{} E_{n}^{H}

(35)

Herein,

Λ_{s}^{}

is a diagonal matrix comprised of the largest K eigenvalues of

R_{y}^{}

E_{s}^{}

is the signal subspace matrix which contains the eigenvectors corresponding to the largest K eigenvalues.

E_{n}^{}

is the noise subspace matrix, which is composed of the remaining 2M-K eigenvectors of

R_{y}^{}

Based on the subspace principle that the noise subspace is orthogonal to the signal subspace, we can construct the following MUSIC spectrum function for estimating the azimuth and elevation angles of the sources.

P (θ, ϕ, φ, C) = a_{E}^{T} (θ, ϕ, φ, C) E_{n} E_{n}^{T} a_{E} (θ, ϕ, φ, C)

(36)

Unfortunately, this multidimensional spectrum function presented above requires a (P + 2)-dimension joint spectral search for 2D DOAs, MCCs, and noncircular phases, which is computationally prohibitive. Therefore, to decrease the complexity of the algorithm, we first factorize the augmented steering vector

a_{E}^{} (θ_{k}^{}, ϕ_{k}^{}, φ_{k}^{}, C)

with respect to the y-axis direction-cosine and present a novel real-valued RARE estimator in the next subsection, which only requires a simple 1D spectral search procedure.

3.2. First Stage Rank Reduction DOA Estimation

In this stage, the real-valued steering vector

a_{E}^{}

in Equation (32) is first factorized with respect to the y-axis direction-cosine and other nuisance parameters. From Equations (17) and (30), the mutual coupling effect in the two virtual steering vectors can be deemed as angle-dependent gain and phase perturbations. However, it is clear that the central elements on the main diagonal of both

Γ_{c, p} (ψ_{x})

and

Γ_{s, p} (ψ_{x})

comprise a vector of all ones, which implies that the middle subarray can be handled as an unperturbed one. Therefore, we can obtain the following equations

Γ_{c, p} (ψ_{x}) a_{x, c} (ψ_{x}) = T_{x, c} (ψ_{x}) α_{x, c}^{(p)}

(37)

Γ_{c, p} (ψ_{x}) a_{x, s} (ψ_{x}) = T_{x, s} (ψ_{x}) α_{x, c}^{(p)}

(38)

Γ_{s, p} (ψ_{x}) a_{x, s} (ψ_{x}) = T_{x, s} (ψ_{x}) α_{x, s}^{(p)}

(39)

Γ_{s, p} (ψ_{x}) a_{x, c} (ψ_{x}) = T_{x, c} (ψ_{x}) α_{x, s}^{(p)}

(40)

where

T_{x, c} (ψ_{x})

and

T_{x, s} (ψ_{x})

are

M_{x} \times (2 P - 1)

real matrices, and

α_{x, c}^{(p)}

and

α_{x, s}^{(p)}

are

(2 P - 1) \times 1

real vectors, which can be expressed as

T_{x, c} (ψ_{x}) = [\begin{matrix} 1 \\ \cos (ψ_{x}) & 0 \\ ⋱ \\ \cos ((P - 1) ψ_{x}) \\ ⋮ \\ \cos ((M_{x} - P) ψ_{x}) \\ 0 & ⋱ \\ \cos ((M_{x} - 1) ψ_{x}) \end{matrix}]

(41)

T_{x, s} (ψ_{x}) = [\begin{matrix} 0 \\ \sin (ψ_{x}) & 0 \\ ⋱ \\ \sin ((P - 1) ψ_{x}) \\ ⋮ \\ \sin ((M_{x} - P) ψ_{x}) \\ 0 & ⋱ \\ \sin ((M_{x} - 1) ψ_{x}) \end{matrix}]

(42)

α_{x, c}^{(p)} = h_{c}^{(p)} {[α_{c, 1}^{(p)}, \dots, α_{c, P - 1}^{(p)}, 1, β_{c, 1}^{(p)}, \dots, β_{c, P - 1}^{(p)}]}^{T}

(43)

α_{x, s}^{(p)} = h_{s}^{(p)} {[α_{s, 1}^{(p)}, \dots, α_{s, P - 1}^{(p)}, 1, β_{s, 1}^{(p)}, \dots, β_{s, P - 1}^{(p)}]}^{T}

(44)

Similarly, for the y-axis direction-cosine related steering vector, we have that

J_{p} a_{y, c} (ψ_{y}) = T_{y, c} (ψ_{y}) α_{y, c}^{(p)} - T_{y, s} (ψ_{y}) α_{y, s}^{(p)}

(45)

J_{p} a_{y, s} (ψ_{y}) = T_{y, s} (ψ_{y}) α_{y, c}^{(p)} + T_{y, c} (ψ_{y}) α_{y, s}^{(p)}

(46)

where

T_{y, c} (ψ_{y})

and

T_{y, s} (ψ_{y})

are

M_{y} \times (2 P - 1)

real matrices, and

α_{y, c}^{(p)}

and

α_{y, s}^{(p)}

are

(2 P - 1) \times 1

real vectors, which can be expressed as

T_{y, c} (ψ_{y}) = [\begin{matrix} 1 \\ \cos (ψ_{y}) & 0 \\ ⋱ \\ \cos ((P - 1) ψ_{y}) \\ ⋮ \\ \cos ((M_{y} - P) ψ_{y}) \\ 0 & ⋱ \\ \cos ((M_{y} - 1) ψ_{y}) \end{matrix}]

(47)

T_{y, s} (ψ_{y}) = [\begin{matrix} 0 \\ \sin (ψ_{y}) & 0 \\ ⋱ \\ \sin ((P - 1) ψ_{y}) \\ ⋮ \\ \sin ((M_{y} - P) ψ_{y}) \\ 0 & ⋱ \\ \sin ((M_{y} - 1) ψ_{y}) \end{matrix}]

(48)

For simplicity, we also denote

T_{y, c} (ψ_{y})

T_{y, c}

in the following derivation. Based on the relationships derived above, we can factorize the virtual steering vectors

v_{R}

and

v_{I}

in Equations (17) and (30) as

\begin{array}{l} v_{R} (θ, ϕ, φ, C) = \\ (T_{y, c} α_{y, c}) \otimes (T_{x, c} α_{x, c}) - (T_{y, c} α_{y, c}) \otimes (T_{x, s} α_{x, s}) - (T_{y, s} α_{y, s}) \otimes (T_{x, c} α_{x, c}) + (T_{y, s} α_{y, s}) \otimes (T_{x, s} α_{x, s}) \\ - (T_{y, s} α_{y, c}) \otimes (T_{x, s} α_{x, c}) - (T_{y, s} α_{y, c}) \otimes (T_{x, c} α_{x, s}) - (T_{y, c} α_{y, s}) \otimes (T_{x, s} α_{x, c}) - (T_{y, c} α_{y, s}) \otimes (T_{x, c} α_{x, s}) \end{array}

(49)

\begin{array}{l} v_{I} (θ, ϕ, φ, C) = \\ (T_{y, c} α_{y, c}) \otimes (T_{x, c} α_{x, s}) + (T_{y, c} α_{y, c}) \otimes (T_{x, s} α_{x, c}) - (T_{y, s} α_{y, s}) \otimes (T_{x, c} α_{x, s}) - (T_{y, s} α_{y, s}) \otimes (T_{x, s} α_{x, c}) \\ - (T_{y, s} α_{y, c}) \otimes (T_{x, s} α_{x, s}) + (T_{y, s} α_{y, c}) \otimes (T_{x, c} α_{x, c}) - (T_{y, c} α_{y, s}) \otimes (T_{x, s} α_{x, s}) + (T_{y, c} α_{y, s}) \otimes (T_{x, c} α_{x, c}) \end{array}

(50)

where

α_{x, c}

α_{x, s}

α_{y, c}

, and

α_{y, s}

are defined as

α_{x, c} = \sum_{p = 1}^{P} α_{x, c}^{(p)}, α_{x, s} = \sum_{p = 1}^{P} α_{x, s}^{(p)}

(51)

α_{y, c} = \sum_{p = 1}^{P} α_{y, c}^{(p)}, α_{y, s} = \sum_{p = 1}^{P} α_{y, s}^{(p)}

(52)

Therefore, utilizing the Kronecker product property

(A B) \otimes (C D) = (A \otimes C) (B \otimes D)

, we can write the virtual steering vector

a_{E}^{}

as the following expression

\begin{array}{l} a_{E} (ψ_{x}, ψ_{y}, φ, C) & = [\begin{matrix} v_{R} \\ v_{I} \end{matrix}] \\ = [\begin{matrix} (T_{y, c} \otimes I_{x}) (I_{y} \otimes T_{x, c}) - (T_{y, s} \otimes I_{x}) (I_{y} \otimes T_{x, s}) & - (T_{y, c} \otimes I_{x}) (I_{y} \otimes T_{x, s}) - (T_{y, s} \otimes I_{x}) (I_{y} \otimes T_{x, c}) \\ (T_{y, c} \otimes I_{x}) (I_{y} \otimes T_{x, s}) + (T_{y, s} \otimes I_{x}) (I_{y} \otimes T_{x, c}) & (T_{y, c} \otimes I_{x}) (I_{y} \otimes T_{x, c}) - (T_{y, s} \otimes I_{x}) (I_{y} \otimes T_{x, s}) \end{matrix}] [\begin{matrix} α_{y, c} \otimes α_{x, c} - α_{y, s} \otimes α_{x, c} \\ α_{y, c} \otimes α_{x, s} - α_{y, s} \otimes α_{x, c} \end{matrix}] \end{array}

(53)

where I_x is a M_x × M_x identity matrix and I_y is a (2P − 1) × (2P − 1) identity matrix. By factorizing the first matrix, we can reorganize

a_{E}^{}

in the following form

\begin{array}{l} a_{E} (ψ_{x}, ψ_{y}, φ, C) & = [\begin{matrix} v_{R} \\ v_{I} \end{matrix}] \\ = \underset{T_{1} (ψ_{y})}{\underset{︸}{[\begin{matrix} T_{y, c} \otimes I_{x} & - T_{y, s} \otimes I_{x} \\ T_{y, s} \otimes I_{x} & T_{y, c} \otimes I_{x} \end{matrix}]}} \underset{α_{1} (ψ_{x}, ψ_{y}, φ, C)}{\underset{︸}{[\begin{matrix} I_{y} \otimes T_{x, c} & - I_{y} \otimes T_{x, s} \\ I_{y} \otimes T_{x, s} & I_{y} \otimes T_{x, c} \end{matrix}] [\begin{matrix} α_{y, c} \otimes α_{x, c} - α_{y, s} \otimes α_{x, c} \\ α_{y, c} \otimes α_{x, s} - α_{y, s} \otimes α_{x, c} \end{matrix}]}} \end{array}

(54)

where

T_{1} (ψ_{y})

is a

2 M_{x} M_{y} \times 2 M_{x} (2 P - 1)

matrix which only depends on the y-axis direction-cosine, and

α_{1} (ψ_{x}, ψ_{y}, φ, C)

is a

2 M_{x} (2 P - 1) \times 1

vector which depends on the x-axis and y-axis direction-cosines, MCCs, and noncircular phases.

Substituting the above equation into Equation (36), we have

α_{1}^{T} (ψ_{x}, ψ_{y}, φ, C) Q_{1} (ψ_{y}) α_{1} (ψ_{x}, ψ_{y}, φ, C) = 0

(55)

where

Q_{1} (ψ_{y})

is defined as

Q_{1} (ψ_{y}) = T_{1}^{T} (ψ_{y}) E_{n} E_{n}^{T} T_{1} (ψ_{y})

(56)

It is obvious that

Q_{1} (ψ_{y})

is a

2 M_{x} (2 P - 1) \times 2 M_{x} (2 P - 1)

real matrix,

T_{1} (ψ_{y})

is a

2 M_{x} M_{y} \times 2 M_{x} (2 P - 1)

real matrix, and

E_{n}

is the

2 M_{x} M_{y} \times (2 M_{x} M_{y} - K)

real-valued noise subspace matrix. Accordingly, if

2 M_{x} M_{y} - K \geq 2 M_{x} (2 P - 1)

Q_{1} (ψ_{y})

will generally be of full column rank. Note that

α_{1}

is presumed to be a non-zero vector, thus Equation (55) holds only when

Q_{1} (ψ_{y})

rank drops, which implies that

Q_{1} (ψ_{y})

becomes a rank-deficient matrix when

ψ_{y} = ψ_{y, k}, k = 1, \dots, K

. Therefore, based on the principle of RARE, the y-axis direction-cosines can be estimated from the highest peaks of the following 1D spectral search function:

P_{1} (ψ_{y}) = \frac{1}{\det [T_{1}^{T} (ψ_{y}) E_{n} E_{n}^{T} T_{1} (ψ_{y})]}

(57)

Since the y-axis direction-cosine

ψ_{y}

has been decoupled with other nuisance parameters in this stage, we can simplify the previous multidimensional optimization problem in Equation (36) into a 1D spectral search problem, which is computationally more efficient. Additionally, the operations involved in this stage are all real-valued, which means the computational load can be reduced by a factor of four compared with its complex-valued counterparts.

Our next objective is to estimate the values of x-axis direction-cosines and correctly pair them with the corresponding y-axis direction-cosines. Therefore, in the following subsection, we will present another RARE function which can generate the estimation of x-axis direction-cosines

ψ_{x, k}, k = 1, \dots, K

without any additional matching procedures.

3.3. Second Stage Rank Reduction DOA Estimation

In the second stage, the real-valued steering vector

a_{E}^{}

in Equation (32) is factorized in another form with respect to the x-axis direction-cosine and other nuisance parameters. Since the y-axis direction-cosine

ψ_{y}

has already been obtained in the previous stage, we can calculate the values of

\sum_{p = 1}^{P} (J_{p} a_{y, c} (ψ_{y, k}))

and

\sum_{p = 1}^{P} (J_{p} a_{y, s} (ψ_{y, k}))

. Also we denote that

β_{y, c} = \sum_{p = 1}^{P} (J_{p} a_{y, c} (ψ_{y, k})), β_{y, s} = \sum_{p = 1}^{P} (J_{p} a_{y, s} (ψ_{y, k}))

(58)

Therefore, the virtual steering vectors

v_{R}

and

v_{I}

in Equations (17) and (30) can be also factorized as

v_{R} (θ_{k}, ϕ_{k}, φ_{k}, C) = β_{y, c} \otimes [T_{x, c} α_{x, c}^{} - T_{x, s} α_{x, s}^{}] - β_{y, s} \otimes [T_{x, s} α_{x, c}^{} + T_{x, c} α_{x, s}^{}]

(59)

v_{I} (θ_{k}, ϕ_{k}, φ_{k}, C) = β_{y, c} \otimes [T_{x, c} α_{x, s}^{} + T_{x, s} α_{x, c}^{}] - β_{y, s} \otimes [T_{x, s} α_{x, s}^{} - T_{x, c} α_{x, c}^{}]

(60)

Utilizing the Kronecker product property

β \otimes (T \cdot α) = (β \otimes T) \cdot α

, we have the following expression

v_{R} (θ_{k}, ϕ_{k}, φ_{k}, C) = (β_{y, c} \otimes T_{x, c}) α_{x, c}^{} - (β_{y, c} \otimes T_{x, s}) α_{x, s}^{} - (β_{y, s} \otimes T_{x, s}) α_{x, c}^{} - (β_{y, s} \otimes T_{x, c}) α_{x, s}^{}

(61)

v_{I} (θ_{k}, ϕ_{k}, φ_{k}, C) = (β_{y, c} \otimes T_{x, c}) α_{x, s}^{} + (β_{y, c} \otimes T_{x, s}) α_{x, c}^{} - (β_{y, s} \otimes T_{x, s}) α_{x, s}^{} + (β_{y, s} \otimes T_{x, c}) α_{x, c}^{}

(62)

where

I_{β}

is a M_y × M_y identity matrix. Consequently, the extended real-valued steering vector

a_{E}

can be written in the following form

a_{E} (ψ_{x}, ψ_{y}, φ, C) = [\begin{matrix} v_{R} \\ v_{I} \end{matrix}] = \underset{T_{2} (ψ_{x})}{\underset{︸}{[\begin{matrix} β_{y, c} \otimes T_{x, c} - β_{y, s} \otimes T_{x, s} & - β_{y, c} \otimes T_{x, s} - β_{y, s} \otimes T_{x, c} \\ β_{y, c} \otimes T_{x, s} + β_{y, s} \otimes T_{x, c} & β_{y, c} \otimes T_{x, c} - β_{y, s} \otimes T_{x, s} \end{matrix}]}} \underset{α_{2} (ψ_{x}, φ, C)}{\underset{︸}{[\begin{matrix} α_{x, c} \\ α_{x, s} \end{matrix}]}}

(63)

where

T_{2} (ψ_{x})

is a

2 M_{x} M_{y} \times 2 (2 P - 1)

matrix which only depends on the x-axis direction-cosine, and

α_{2} (ψ_{x}, φ, C)

is a

2 (2 P - 1)

vector which depends on the x-axis direction-cosines, MCCs, and noncircular phases.

Again, according to the RARE principle, we substitute the above equation into Equation (36)

α_{2}^{T} (ψ_{x}, φ, C) Q_{2} (ψ_{x}) α_{2} (ψ_{x}, φ, C) = 0

(64)

where

Q_{2} (ψ_{x})

is defined as

Q_{2} (ψ_{x}) = T_{2}^{T} (ψ_{x}) E_{n} E_{n}^{T} T_{2} (ψ_{x}) .

(65)

Note that

Q_{2} (ψ_{x}) \in ℝ^{2 (2 P - 1) \times 2 (2 P - 1)}

and

T_{2} (ψ_{x}) \in ℝ^{2 M_{x} M_{y} \times 2 (2 P - 1)}

. If

2 M_{x} M_{y} - K \geq 2 (2 P - 1)

, then

Q_{2} (ψ_{x})

is of full column rank in general. As both

α_{x, c}

and

α_{x, s}

are supposed to be non-zero vectors,

α_{2}

must be non-zero as well. One can deduce that (64) holds only if the rank of

Q_{2} (ψ_{x})

drops, that is, rank deficiency occurs in

Q_{2} (ψ_{x})

provided that

ψ_{x} = ψ_{x, k}, k = 1, \dots, K

. By the principle of RARE, the x-axis direction-cosines can be obtained by searching for the highest peaks of the following spatial spectrum function:

P_{2} (ψ_{x}) = \frac{1}{\det [T_{2}^{T} (ψ_{x}) E_{n} E_{n}^{T} T_{2} (ψ_{x})]} .

(66)

Since the x-axis direction-cosine

ψ_{x}

has been decoupled with other nuisance parameters in the second stage, the multidimensional optimization problem can be reduced to a computationally efficient 1D spectral searching problem. Moreover, the computations involved are all real-valued, which means the computational load can be reduced by a factor of four compared with its complex-valued counterparts. Also, due to the orthogonality between the signal subspace and noise subspace, the 2D DOA parameters are automatically paired without any additional matching procedures.

Therefore, the 2D DOA estimates can be obtained from the following equations:

{\hat{θ}}_{k} = arc \cot ((\frac{λ {\hat{ψ}}_{x}}{2 π d_{x}}) / (\frac{λ {\hat{ψ}}_{y}}{2 π d_{y}}))

(67)

{\hat{ϕ}}_{k} = \arcsin (\sqrt{{(\frac{λ {\hat{ψ}}_{x}}{2 π d_{x}})}^{2} + {(\frac{λ {\hat{ψ}}_{y}}{2 π d_{y}})}^{2}})

(68)

3.4. Procedures of the Proposed Algorithm

The steps of the proposed algorithm can be summarized as follows.

(1) Collect the real and imaginary parts of x(t) based on Equation (12);

(2) construct the augmented data vector via Equation (31);

(3) calculate the covariance matrix of y(t) by Equation (34);

(4) perform the eigenvalue decomposition of

{\hat{R}}_{y}^{}

, and extract the noise subspace matrix

{\hat{E}}_{n}^{}

;

(5) formulate the matrix

T_{1} (ψ_{y})

for each y-axis direction-cosine searching grid according to Equation (54);

(6) find the estimation of y-axis direction-cosines

ψ_{y, k}, k = 1, \dots, K

through the 1D spectral search in Equation (57);

(7) compute the values of

β_{y, c}

β_{y, s}

and

T_{2} (ψ_{x})

for each x-axis direction-cosines estimates based on Equations (58) and (63);

(8) obtain the estimates of x-axis direction-cosines

ψ_{x, k}, k = 1, \dots, K

according to Equation (66);

(9) calculate the DOA estimates of the azimuth and elevation angles according to Equations (67) and (68).

4. Discussion

In this section, we compare the performance of the proposed algorithm with respect to two state-of-art 2D DOA estimation algorithms considering mutual coupling effects: the auxiliary sensor-based method (AUX) in [39] and the 2D rank reduction-based method (2D-RARE) in [40].

4.1. Computational Complexity

In this subsection, the computational complexity of the proposed algorithm is evaluated with respect to AUX and 2D-RARE. Note that the computational complexity is investigated in terms of the major multiplications in constructing the statistical matrices, EVDs, and spectral searches. Suppose that an M_x × M_y URA is utilized and the grid size for spectral searching is

Δ θ

. From [29], it can be observed that AUX constructs a (M_x − 2P + 2)(M_y − 2P + 2) × (M_x − 2P + 2)(M_y − 2P + 2) covariance matrix, takes the EVD of this matrix, and performs a 2D spectral search for DOA estimation, which require a computational complexity of

O (4 {(M_{x} - 2 P + 2)}^{2} {(M_{y} - 2 P + 2)}^{2} L)

O (4 {(M_{x} - 2 P + 2)}^{3} {(M_{y} - 2 P + 2)}^{3} / 3)

, and

O (4 {(M_{x} - 2 P + 2)}^{2} {(M_{y} - 2 P + 2)}^{2} (360 ° / Δ θ) (90 ° / Δ θ))

, respectively. 2D-RARE constructs a M_xM_y × M_xM_y covariance matrix, performs the EVD, calculates the rank reduction testing matrix, and performs a 2D RARE spectral search, which require a computational complexity of

O (4 {(M_{x} M_{y})}^{2} L)

O (4 {(M_{x} M_{y})}^{3} / 3)

, and

O ((360 ° / Δ θ) (90 ° / Δ θ)) (M_{x} M_{y} (M_{x} M_{y} - K) {(2 P - 1)}^{2} + (M_{x} M_{y} - K) {(2 P - 1)}^{4} + {(M_{x} M_{y})}^{3}))

, respectively [30]. For the proposed algorithm, the establishment of the 2M_xM_y × 2M_xM_y real-valued statistical matrix takes

O (4 {(M_{x} M_{y})}^{2} L)

operations, the EVD of this matrix requires

O (8 {(M_{x} M_{y})}^{3} / 3)

multiplications, the first stage of the 1D RARE spectral search takes

O ((360 ° / Δ θ) {(2 M_{x} (2 P - 1))}^{3})

, and the second stage of the 1D RARE spectral search has a complexity of

O (K \cdot (90 ° / Δ θ) {(4 P - 2)}^{3})

. It is worth noting that the calculations involved in the proposed algorithm are in the real-valued domain. Thus, with the same matrix dimensions, the computational load can be effectively reduced by a factor of four.

4.2. Maximum Number of Resolvable Sources

In this subsection, we compare the maximum number of resolvable sources for the three algorithms. Note that AUX only utilizes the (M_x − 2P + 2)(M_y − 2P + 2) sensors in the middle of the URA. As a result of the subspace identification technique, at least one eigenvector is required for spanning the noise subspace. Therefore, AUX is capable of resolving no more than (M_x − 2P + 2)(M_y − 2P + 2) − 1 sources. Similarly, 2D-RARE requires that (2P − 1)² ≤ M_xM_y-K. Thus, it can also resolve M_xM_y − (2P − 1)² sources at most. For the proposed method, it requires

2 M_{x} M_{y} - K \geq 2 M_{x} (2 P - 1)

and

2 M_{x} M_{y} - K \geq 2 (2 P - 1)

in the two consequent RARE estimators, which means the maximum number of resolvable sources for the proposed method is

2 M_{x} M_{y} - 2 M_{x} (2 P - 1)

. Therefore, it can be concluded that the alleviation of computational loads is at the cost of DOF reduction, since the signal subspace has been decoupled.

4.3. Dependence of Mutual Coupling on Directions

In this subsection, we discuss the dependence of antenna mutual coupling effects on the angles of signals. According to the latest work by B. Friedlander [28], the mutual coupling matrix is a direction dependent matrix, and not a constant matrix as is commonly assumed. In this work, the mutual coupling matrix of the URA is assumed to be direction-independent with a block banded Toeplitz structure, which has been simplified compared with the analytic manifold [28]. Therefore, the DOA estimation performance may deteriorate due to the mismatch between the assumed model and the physical characteristics of the array antennas.

5. Simulation Results

In this section, we present several numerical simulations to analyze the DOA estimation performance of the proposed algorithm with respect to 2D-MUSIC, AUX, and 2D-RARE. In the following simulations, a 10 × 10-elements URA is employed with the inter-element spacing being d_x = d_y =

λ

/2. The impinging sources are BPSK modulated with equal power and uncorrelated to each other, while the noise is amenable to the spatially white complex Gaussian distribution. The mutual coupling length is assumed to be

P = 2

, and MCCs are

c_{0, 1} = c_{1, 0} = 0.7527 + j 0.4854

c_{1, 1} = 0.2825 + j 0.2801

. The estimation performance is assessed in terms of spatial spectrum, root mean-square error (RMSE), and computational complexity. The RMSE is measured by 500 independent Monte Carlo runs:

R M S E = \sqrt{\frac{1}{1000 K} \sum_{n = 1}^{500} \sum_{k = 1}^{K} ({(θ_{k} - {\hat{θ}}_{n, k})}^{2} + {(ϕ_{k} - {\hat{ϕ}}_{n, k})}^{2})}

(69)

where

θ_{k}

is the DOA of the k-th source, and

{\hat{θ}}_{n, k}

stands for the estimate of in the nth trial.

In the first case, we investigate the DOA spectra of the three algorithms. We consider that three incident rectilinear signals radiate from (

θ_{1} = - 20 °, ϕ_{1} = 40 °

), (

θ_{2} = 40 °, ϕ_{2} = 30 °

), and (

θ_{3} = 80 °, ϕ_{3} = 50 °

), respectively. The corresponding x-axis direction-cosines are (0.604, 0.383, 0.133) in radians, and y-axis direction-cosines are (–0.220, 0.321, 0.754) in radians, respectively. The signal-to-noise ratio (SNR) is 10 dB, and the number of snapshots is fixed at 1000. In Figure 2d, the proposed algorithm generates three sharp peaks corresponding to the y-axis direction-cosines of the three BPSK sources in the first DOA estimation stage. In Figure 2e, the three spectra produce three sharp peaks at the desired locations of the x-axis direction-cosines, respectively. Note that the DOA estimates of in the two stages are automatically paired without any postprocessing. From Figure 2a, it is can be seen that 2D-MUSIC fails to distinguish all three sources, and its two resolvable spectra peaks are biased due to the mismatch caused by the mutual coupling effects. Although AUX and 2D-RARE can generate three peaks corresponding to the sources in the 2D DOA spectrum plain, the former sacrifices eight sensors (i.e., the auxiliary sensors) to make the remaining part of the URA coupling-free while the latter requires exhaustive 2D spectral searches to obtain the joint azimuth and elevation DOA estimates.

In the second experiment, we studied the RMSE of the DOA estimates as a function of the SNR. The number of snapshots remained at 1000. Other simulation settings were the same as the first experiment. Figure 3 depicts the RMSE of the proposed method while the RMSEs of the other three solutions, 2D-MUSIC, AUX, 2D-RARE, and the corresponding CRB for noncircular source DOA estimation in the presence of mutual coupling is also calculated for comparison. It is clear that 2D-RARE performs the best among the four algorithms, overall AUX takes second place, our method provides better accuracies than it at 10 dB and above but is inferior to AUX at low SNRs, and 2D-MUSIC has the worst performance due to the low robustness to the unknown mutual coupling. Besides, the performance of the proposed method, 2D-RARE, and AUX bears down on the CRB asymptotically at high SNRs, while the RMSE of 2D-MUSIC is stabilized at approximately

4^{\circ}

when the SNR is larger than 5dB.

In the third experiment, we examine the variation of RMSE with the increase of the number of snapshots from 10 to 1000, where the SNR is set as 10dB. Figure 4 demonstrates that RMSEs of the proposed method, 2D-RARE, and AUX decline steadily as the snapshot size raises. This is because that the covariance and elliptic covariance matrices can be estimated more accurately from a larger number of observations. Additionally, the performance of our developed approach is almost the same as 2D-RARE and strictly superior to AUX, asymptotically approaching the corresponding CRB, while the RMSE of 2D-MUSIC saturates at

4^{\circ}

as similar to the second scenario, which is the worst performance.

In the last experiment, the computational burdens of the three algorithms versus the variation of the number of sensors were compared. The searching step size for DOA estimation was set as 0.1°, and the number of snapshots was fixed at 200. Suppose the URA is square, i.e.,

M_{x} = M_{y}

, and the number of elements ranges from 10 to 40. It can be observed from Figure 5 that our method is the most efficient compared with other two as it is only reliant upon a 1D spectral search, followed by AUX, and 2D-RARE requires the highest complexity.

6. Conclusions

In this paper, a 2D real-valued DOA estimator that is robust to unknown mutual coupling was developed for rectilinear sources. Taking advantage of the noncircularity in signals and the structured mutual coupling matrix revealed as block Toeplitz, a RARE estimator in real-value field was proposed to extract the angular information from the coupling coefficients. Compared with state-of-the-art approaches, our solution has satisfactory estimation performance with more efficient computations. A series of numerical simulations showed that the proposed method has its own advantages on estimation accuracy and computational complexity over previous algorithms with a favorable robustness to unknown mutual coupling. To apply the proposed algorithm to further applications, we may extend the proposed algorithm to considering the physical characteristics of the array antennas, such as direction-dependent mutual coupling and unstructured mutual coupling matrix in further studies.

Author Contributions

All authors contributed extensively to this paper. Conceptualization, J.X. and Y.W.; methodology, J.X. and Q.W.; software, X.Y. and Q.W.; writing—original draft preparation, Y.W., and X.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China under Grant No.61971355, 61971357, 61901390, 61771404, 61871459, and in part by the Natural Science Basic Research Plan in Shaanxi Province of China under Grant 2019JM-350.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In the Appendix, we derive the unconditional or stochastic CRB of DOA estimation for noncircular sources under unknown mutual coupling effects. Since the signals are rectilinear (also called strictly second-order noncircular),

s_{0} (t)

is zero-mean real-valued Gaussian distributed with covariance

E (s_{0} (t) s_{0}^{T} (t)) = R_{s}

. Furthermore, the model can be parameterized by the following real-valued unknown parameter vector

Θ = [θ^{T}, ϕ^{T}, φ^{T}, c_{R}^{T}, c_{I}^{T}, ρ^{T}, σ_{n}^{2}]

, where

θ = {[θ_{1}, \dots, θ_{K}]}^{T}

denotes the azimuth angle vector,

ϕ = {[ϕ_{1}, \dots, ϕ_{K}]}^{T}

represents the elevation angle vector,

φ = {[φ_{1}, \dots, φ_{K}]}^{T}

contains the noncircular phases,

c_{R}

and

c_{I}

are the real and imaginary parts of the mutual coupling coefficients vector, and

ρ = {[ρ_{1}, \dots, ρ_{K}]}^{T}

is the K × 1 vector made from

{[R_{s}]}_{i, i}

, 1≤ i ≤K. Based on the results in [42], the Fisher information matrix (FIM) corresponding to the parameter vector

Θ

can be written as:

F_{m, n} = L \cdot Tr (\frac{\partial R_{\bar{x}}}{\partial Θ_{m}} R_{\bar{x}}^{- 1} \frac{\partial R_{\bar{x}}}{\partial Θ_{n}} R_{\bar{x}}^{- 1})

(A1)

where

R_{\bar{x}}

is the extended covariance matrix corresponding to the augmented signal

\bar{x} (t) = {[x_{}^{T} (t), x_{}^{H} (t)]}^{T}

, which is formed as:

R_{\bar{x}} = E {\bar{x} (t) {\bar{x}}^{H} (t)} = \bar{A} R_{s} {\bar{A}}^{H} + σ_{n}^{2} I

(A2)

where

\bar{A} = [\begin{matrix} C A Φ \\ C^{*} A^{*} Φ^{*} \end{matrix}] = [{\bar{a}}_{1}, \dots, {\bar{a}}_{K}]

, with

{\bar{a}}_{k} = [\begin{matrix} e^{j φ_{k}} C a (θ_{k}, ϕ_{k}) \\ e^{- j φ_{k}} C^{*} a^{*} (θ_{k}, ϕ_{k}) \end{matrix}]

The FIM for the unknown parameters and their cross terms can be readily derived following the steps in [43]. Since F is a symmetric matrix, we only need to find the upper triangular parts of F, and the block matrices in F can be obtained as:

F_{θ θ} = 2 L \cdot Re {(R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{θ}) ⊙ {(R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{θ})}^{T} + (R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} \bar{A} R_{s}) ⊙ {({\dot{\bar{A}}}_{θ}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{θ})}^{T}}

(A3)

F_{ϕ ϕ} = 2 L \cdot Re {(R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{ϕ}) ⊙ {(R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{ϕ})}^{T} + (R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} \bar{A} R_{s}) ⊙ {({\dot{\bar{A}}}_{ϕ}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{ϕ})}^{T}}

(A4)

F_{φ φ} = 2 L \cdot Re {(R_{S} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{φ}^{}) ⊙ {(R_{S} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{φ}^{})}^{T} + (R_{S}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} \bar{A} R_{S}) ⊙ {({\dot{\bar{A}}}_{φ}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{φ})}^{T}}

(A5)

F_{c_{R, m} c_{R, n}} = 2 L \cdot Re {Tr {R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{R, m}}^{} R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{R, n}}^{} R_{s}^{} {\bar{A}}_{}^{H}} + Tr {R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{R, m}}^{} R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} \bar{A} R_{s}^{} {\dot{\bar{A}}}_{c_{R, n}}^{H}}}

(A6)

F_{c_{I, m} c_{I, n}} = 2 L \cdot Re {Tr {R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{I, m}}^{} R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{I, n}}^{} R_{s}^{} {\bar{A}}_{}^{H}} + Tr {R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{I, m}}^{} R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} \bar{A} R_{s}^{} {\dot{\bar{A}}}_{c_{I, n}}^{H}}}

(A7)

F_{c_{R, m} c_{I, n}} = 2 L \cdot Re {Tr {R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{R, m}}^{} R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{I, n}}^{} R_{s}^{} {\bar{A}}_{}^{H}} + Tr {R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{R, m}}^{} R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} \bar{A} R_{s}^{} {\dot{\bar{A}}}_{c_{I, n}}^{H}}}

(A8)

F_{ρ ρ} = L \cdot ({\bar{A}}^{H} R_{\bar{x}}^{- 1} \bar{A}) ⊙ {({\bar{A}}^{H} R_{\bar{x}}^{- 1} \bar{A})}^{T}

(A9)

F_{θ φ} = 2 L \cdot Re {(R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{θ}) ⊙ {(R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{φ})}^{T} + (R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} \bar{A} R_{s}) ⊙ {({\dot{\bar{A}}}_{φ}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{θ})}^{T}}

(A10)

F_{ϕ φ} = 2 L \cdot Re {(R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{ϕ}) ⊙ {(R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{φ})}^{T} + (R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} \bar{A} R_{s}) ⊙ {({\dot{\bar{A}}}_{φ}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{ϕ})}^{T}}

(A11)

F_{θ c_{R, n}} = 2 L \cdot Re {diag {R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{R, n}}^{} R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{θ}^{}} + diag {R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} \bar{A} R_{s}^{} {\dot{\bar{A}}}_{c_{R, n}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{θ}^{}}}

(A12)

F_{ϕ c_{R, n}} = 2 L \cdot Re {diag {R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{R, n}}^{} R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{θ}^{}} + diag {R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} \bar{A} R_{s}^{} {\dot{\bar{A}}}_{c_{R, n}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{θ}^{}}}

(A13)

F_{φ c_{R, n}} = 2 L \cdot Re {diag {R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{R, n}}^{} R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{φ}^{}} + diag {R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} \bar{A} R_{s}^{} {\dot{\bar{A}}}_{c_{R, n}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{φ}^{}}}

(A14)

F_{θ c_{I, n}} = 2 L \cdot Re {diag {R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{I, n}}^{} R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{θ}^{}} + diag {R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} \bar{A} R_{s}^{} {\dot{\bar{A}}}_{c_{I, n}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{θ}^{}}}

(A15)

F_{ϕ c_{I, n}} = 2 L \cdot Re {diag {R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{I, n}}^{} R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{ϕ}^{}} + diag {R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} \bar{A} R_{s}^{} {\dot{\bar{A}}}_{c_{I, n}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{ϕ}^{}}}

(A16)

F_{φ c_{I, n}} = 2 L \cdot Re {diag {R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{I, n}}^{} R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{φ}^{}} + diag {R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} \bar{A} R_{s}^{} {\dot{\bar{A}}}_{c_{I, n}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{φ}^{}}}

(A17)

F_{θ ρ} = 2 L \cdot Re {(R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} \bar{A}) ⊙ {({\bar{A}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{θ})}^{T}}

(A18)

F_{ϕ ρ} = 2 L \cdot Re {(R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} \bar{A}) ⊙ {({\bar{A}}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{ϕ})}^{T}}

(A19)

F_{φ ρ} = 2 L \cdot Re {(R_{s}^{} {\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} \bar{A}) ⊙ {({\bar{A}}_{}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{φ}^{})}^{T}}

(A20)

F_{c_{R, m}, ρ_{n}} = 2 L \cdot Re {{\bar{a}}_{n}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{R, m}} R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} {\bar{a}}_{n}}

(A21)

F_{c_{I, m}, ρ_{n}} = 2 L \cdot Re {{\bar{a}}_{n}^{H} R_{\bar{x}}^{- 1} {\dot{\bar{A}}}_{c_{I, m}} R_{s} {\bar{A}}^{H} R_{\bar{x}}^{- 1} {\bar{a}}_{n}}

(A22)

where the partial derivative of A with respect to the parameters

θ, ϕ, φ, c_{R}^{}, c_{I}^{}

in the above equations can be calculated as:

{\dot{\bar{A}}}_{θ} = [\frac{\partial {\bar{a}}_{1}}{\partial θ_{1}}, \dots, \frac{\partial {\bar{a}}_{K}}{\partial θ_{K}}], {\dot{\bar{A}}}_{ϕ} = [\frac{\partial {\bar{a}}_{1}}{\partial ϕ_{1}}, \dots, \frac{\partial {\bar{a}}_{K}}{\partial ϕ_{K}}]

(A23)

{\dot{\bar{A}}}_{φ} = [\frac{\partial {\bar{a}}_{1}}{\partial φ_{1}}, \dots, \frac{\partial {\bar{a}}_{K}}{\partial φ_{K}}]

(A24)

{\dot{\bar{A}}}_{c_{R,}_{m}} = [\begin{matrix} \frac{\partial C}{\partial c_{R,}_{m}} A Φ \\ \frac{\partial C^{*}}{\partial c_{R,}_{m}} A^{*} Φ^{*} \end{matrix}], {\dot{\bar{A}}}_{c_{I,}_{m}} = [\begin{matrix} \frac{\partial C}{\partial c_{I,}_{m}} A Φ \\ \frac{\partial C^{*}}{\partial c_{I,}_{m}} A^{*} Φ^{*} \end{matrix}]

(A25)

Since the noise power is not coupled with other parameters, the cross terms with respect to

σ_{n}^{2}

are all zeros. Consequently, the FIM can be formulated as the following form:

F = [\begin{array}{l} F_{θ θ} & F_{θ ϕ} & F_{θ φ} & F_{θ c_{r}} & F_{θ c_{I}} & F_{θ ρ} & 0 \\ F_{ϕ θ} & F_{ϕ ϕ} & F_{ϕ φ} & F_{ϕ c_{r}} & F_{ϕ c_{I}} & F_{ϕ ρ} & 0 \\ F_{φ θ} & F_{φ ϕ} & F_{φ φ} & F_{φ c_{r}} & F_{φ c_{I}} & F_{φ ρ} & 0 \\ F_{c_{r} θ} & F_{c_{r} ϕ} & F_{c_{r} φ} & F_{c_{r} c_{r}} & F_{c_{r} c_{I}} & F_{c_{r} ρ} & 0 \\ F_{c_{I} θ} & F_{c_{I} ϕ} & F_{c_{I} φ} & F_{c_{I} c_{r}} & F_{c_{I} c_{I}} & F_{c_{I} ρ} & 0 \\ F_{ρ θ} & F_{ρ ϕ} & F_{ρ φ} & F_{ρ c_{r}} & F_{ρ c_{I}} & F_{ρ ρ} & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & F_{σ_{n}^{2} σ_{n}^{2}} \end{array}]

(A26)

Subsequently, the CRB for DOA parameters is defined as:

C R B_{D O A} = \sqrt{\frac{1}{2 K} \sum_{k = 1}^{2 K} {[F^{- 1}]}_{k k}}

(A27)

References

Ioannides, P.; Balanis, C.A. Uniform circular and rectangular arrays for adaptive beamforming applications. IEEE Antennas Wireless Propag. Lett. 2005, 4, 351–354. [Google Scholar] [CrossRef]
Zhou, C.; Shi, Z.; Gu, Y.; Shen, X. DECOM: DOA Estimation with Combined MUSIC for Coprime Array. In Proceedings of the International Conference on Wireless Communications and Signal Processing (WCSP), Hangzhou, China, 24–26 October 2013; pp. 1–5. [Google Scholar]
Zhou, C.; Gu, Y.; Zhang, Y.D.; Shi, Z.; Jin, T.; Wu, X. Compressive sensing-based coprime array direction-of-arrival estimation. IET Commun. 2017, 11, 1719–1724. [Google Scholar] [CrossRef]
Zoltowski, M.D.; Haardt, M.; Mathews, C.P. Closed-form 2-D angle estimation with rectangular arrays in element space or beamspace via unitary ESPRIT. IEEE Trans. Signal Process. 1996, 44, 316–328. [Google Scholar] [CrossRef] [Green Version]
Wong, K.T.; Zoltowski, M.D. Direction-finding with sparse rectangular dual-size spatial invariance array. IEEE Trans. Aerosp. Electron. Syst. 1998, 34, 1320–1336. [Google Scholar] [CrossRef]
Yilmazer, N.; Sarkar, T.K. 2-D unitary matrix pencil method for efficient direction of arrival estimation. Digital Signal Process. 2006, 16, 767–781. [Google Scholar] [CrossRef]
Zhang, W.; Liu, W.; Wang, J. Computationally efficient 2-D DOA estimation for uniform rectangular arrays. Multidimens. Syst. Signal Process. 2014, 25, 847–857. [Google Scholar] [CrossRef]
Ren, S.; Ma, X.; Yan, S. 2-D unitary ESPRIT-like direction-of-arrival DOA estimation for coherent signals with a uniform rectangular array. Sensors 2013, 13, 4272–4288. [Google Scholar] [CrossRef] [Green Version]
Schreier, P.J.; Scharf, L.L. Statistical Signal Processing of Complex-Valued Data: The Theory of Improper and Noncircular Signals; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
Abeida, H.; Delmas, J.P. MUSIC-like estimation of direction of arrival for noncircular sources. IEEE Trans. Signal Process. 2006, 54, 2678–2690. [Google Scholar] [CrossRef]
Liu, Z.M.; Huang, Z.T.; Zhou, Y.Y. Direction-of-arrival estimation of noncircular signals via sparse representation. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 2690–2698. [Google Scholar] [CrossRef]
Chen, Z.; Ding, Y.; Ren, S. A novel noncircular MUSIC algorithm based on the concept of the difference and sum coarray. Sensors 2018, 18, 344. [Google Scholar] [CrossRef] [Green Version]
Gao, F.; Nallanathan, A.; Wang, Y. Improved MUSIC under the coexistence of both circular and noncircular sources. IEEE Trans. Signal Process. 2008, 56, 3033–3038. [Google Scholar] [CrossRef]
Chen, H.; Hou, C.; Liu, W. Efficient two-dimensional direction-of-arrival estimation for a mixture of circular and noncircular sources. Sensors 2016, 16, 2527–2536. [Google Scholar] [CrossRef] [Green Version]
Wang, Q.; Zhu, X.; Chen, H. Computationally efficient direction finding for a mixture of circular and strictly noncircular sources with uniform rectangular arrays. Sensors 2017, 17, 1269. [Google Scholar] [CrossRef] [PubMed]
Chen, H.; Zhu, W.P.; Liu, W. RARE-based localization for mixed near-field and far-field rectilinear sources. Digital Signal Process. 2019, 85, 54–61. [Google Scholar] [CrossRef]
Makarov, S.N. Antenna and EM Modeling with MATLAB; Wiley-Interscience: Hoboken, NJ, USA, 2002. [Google Scholar]
Ng, B.C.; Samson, C.M. Sensor-array calibration using a maximum-likelihood approach. IEEE Trans. Antennas Propag. 1996, 44, 827–835. [Google Scholar]
Friedlander, B.; Weiss, A.J. Direction finding in the presence of mutual coupling. IEEE Trans. Antennas Propag. 1991, 39, 273–284. [Google Scholar] [CrossRef]
Ye, Z.; Dai, J.; Xu, X. DOA estimation for uniform linear array with mutual coupling. IEEE Trans. Aerosp. Electron. Syst. 2009, 45, 280–288. [Google Scholar]
Liao, B.; Zhang, Z.G.; Chan, S.C. DOA estimation and tracking of ULAs with mutual coupling. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 891–905. [Google Scholar] [CrossRef] [Green Version]
Liao, B.; Chan, S.C. A cumulant-based approach for direction finding in the presence of mutual coupling. Signal Process. 2014, 104, 197–202. [Google Scholar] [CrossRef]
Wang, Y.; Trinkle, M.; Ng, B.W.H. DOA estimation under unknown mutual coupling and multipath with improved effective array aperture. Sensors 2015, 15, 30856–30869. [Google Scholar] [CrossRef]
Wang, Y.; Wang, L.; Xie, J.; Trinkle, M.; Ng, B.W.H. DOA estimation under mutual coupling of uniform linear arrays using sparse reconstruction. IEEE Wireless Commun. Lett. 2019, 8, 1004–1007. [Google Scholar] [CrossRef]
Wang, Y.; Yang, X.; Xie, J.; Wang, L.; Ng, B.W.H. Sparsity-inducing DOA estimation of coherent signals under the coexistence of mutual coupling and nonuniform noise. IEEE Access. 2019, 7, 40271–40278. [Google Scholar] [CrossRef]
Chen, H.; Liu, W.; Zhu, W.P.; Swamy, M.N.S.; Wang, Q. Mixed rectilinear sources localization under unknown mutual coupling. J. Frankl. Inst. 2019, 356, 2372–2394. [Google Scholar]
Wang, Q.; Wang, X.; Chen, H. DOA Estimation Algorithm for Strictly Noncircular Sources with Unknown Mutual Coupling. IEEE Comm. Lett. 2019, 23, 2215–2218. [Google Scholar] [CrossRef]
Friedlander, B. Antenna array manifolds for high-resolution direction finding. IEEE Trans. Signal Process. 2017, 66, 923–932. [Google Scholar] [CrossRef]
Wang, B.H.; Hui, H.T.; Leong, M.S. Decoupled 2D direction of arrival estimation using compact uniform circular arrays in the presence of elevation-dependent mutual coupling. IEEE Trans. Antennas Propag. 2009, 58, 747–755. [Google Scholar] [CrossRef]
Elbir, A.M. Direction finding in the presence of direction-dependent mutual coupling. IEEE Antennas Wirel. Propag. Lett. 2017, 16, 1541–1544. [Google Scholar] [CrossRef]
Ge, Q.; Zhang, Y.; Wang, Y. A Low Complexity Algorithm for Direction of Arrival Estimation with Direction-Dependent Mutual Coupling. IEEE Comm. Lett. 2019, 24, 90–94. [Google Scholar] [CrossRef]
Shi, Z.; Zhou, C.; Gu, Y.; Goodman, N.A.; Qu, F. Source estimation using coprime array: A sparse reconstruction perspective. IEEE Sens. J. 2017, 17, 573–586. [Google Scholar] [CrossRef]
Zhou, C.; Gu, Y.; Fan, X.; Shi, Z.; Mao, G.; Zhang, Y.D. Direction-of-arrival estimation for coprime array via virtual array interpolation. IEEE Trans. Signal Process. 2018, 66, 5956–5971. [Google Scholar] [CrossRef]
Zhou, C.; Gu, Y.; He, S.; Shi, Z. A robust and efficient algorithm for coprime Array adaptive beamforming. IEEE Trans. Veh. Technol. 2018, 67, 1099–1112. [Google Scholar] [CrossRef]
Zhou, C.; Gu, Y.; Shi, Z.; Zhang, Y.D. Off-grid direction-of-arrival estimation using coprime array interpolation. IEEE Signal Process. Lett. 2018, 25, 1710–1714. [Google Scholar] [CrossRef]
Lin, M.; Yang, L. Blind calibration and DOA estimation with uniform circular arrays in the presence of mutual coupling. IEEE Antennas Wirel. Propag. 2006, 5, 315–318. [Google Scholar] [CrossRef]
Liu, S.; Yang, L.; Yang, S. Robust joint calibration of mutual coupling and channel gain/phase inconsistency for uniform circular array. IEEE Antennas Wireless Propag. Lett. 2015, 15, 1191–1195. [Google Scholar] [CrossRef]
Ye, Z.; Liu, C. 2-D DOA estimation in the presence of mutual coupling. IEEE Trans. Antennas Propag. 2008, 56, 3150–3158. [Google Scholar] [CrossRef]
Liu, C.; Ye, Z.; Zhang, Y. Autocalibration algorithm for mutual coupling of planar array. Signal Process. 2010, 90, 784–794. [Google Scholar] [CrossRef]
Wu, H.; Hou, C.; Chen, H. Direction finding and mutual coupling estimation for uniform rectangular arrays. Signal Process. 2015, 117, 61–68. [Google Scholar] [CrossRef]
Steeb, W.H.; Hardy, Y. Matrix Calculus and Kronecker Product: A Practical Approach to Linear and Multilinear Algebra, 2nd ed.; World Scientific Publishing Company: Singapore, 2011. [Google Scholar]
Stoica, P.; Larsson, E.G.; Gershman, A.B. The stochastic CRB for array processing: A textbook derivation. IEEE Signal Process. Lett. 2001, 8, 148–150. [Google Scholar] [CrossRef]
Xie, J.; Wang, L.; Wang, Y. Efficient Real-Valued Rank Reduction Algorithm for DOA Estimation of Noncircular Sources Under Mutual Coupling. IEEE Access 2018, 6, 64450–64460. [Google Scholar] [CrossRef]

Figure 1. The structure of a uniform rectangular array (URA) with M elements.

Figure 2. The spatial spectra of (a) Two-dimensional multiple signal classification (2D-MUSIC), (b) the auxiliary sensor-based method (AUX), (c) the 2D rank reduction-based method (2D-RARE) and (d) the first stage and (e) the second stage of the proposed method.

Figure 3. RMSEs of Root mean square errors (RMSEs) of direction-of-arrival (DOA) estimates versus signal-to-noise ratio (SNR). The number of snapshots is 1000.

Figure 4. RMSEs of DOA estimates versus the number of snapshots. The SNR is 10 dB.

Figure 5. Computational complexity of the four methods versus number of antennas.

Table 1. List of notations.

Notations	Explanation
${(\cdot)}^{*}$	conjugate operator
${(\cdot)}^{T}$	transpose operator
${(\cdot)}^{H}$	conjugate transpose operator
Re{·}	real part of the operand
Im{·}	imaginary part of the operand
$E {\cdot}$	expectation operator
vec(·)	vectorization operator
e	Hadamard product operator
$\otimes$	Kronecker product operator
Tr{·}	trace of the operand
det(·)	determinant of the operand
I_n	identity matrix of order n

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Xie, J.; Wang, Q.; Wang, Y.; Yang, X. Efficient Two-Dimensional Direction Finding Algorithm for Rectilinear Sources Under Unknown Mutual Coupling. Sensors 2020, 20, 1914. https://doi.org/10.3390/s20071914

AMA Style

Xie J, Wang Q, Wang Y, Yang X. Efficient Two-Dimensional Direction Finding Algorithm for Rectilinear Sources Under Unknown Mutual Coupling. Sensors. 2020; 20(7):1914. https://doi.org/10.3390/s20071914

Chicago/Turabian Style

Xie, Jian, Qiuping Wang, Yuexian Wang, and Xin Yang. 2020. "Efficient Two-Dimensional Direction Finding Algorithm for Rectilinear Sources Under Unknown Mutual Coupling" Sensors 20, no. 7: 1914. https://doi.org/10.3390/s20071914

APA Style

Xie, J., Wang, Q., Wang, Y., & Yang, X. (2020). Efficient Two-Dimensional Direction Finding Algorithm for Rectilinear Sources Under Unknown Mutual Coupling. Sensors, 20(7), 1914. https://doi.org/10.3390/s20071914

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Efficient Two-Dimensional Direction Finding Algorithm for Rectilinear Sources Under Unknown Mutual Coupling

Abstract

1. Introduction

2. Signal Model

3. Proposed Real-Valued 2D DOA Estimation Algorithm

3.1. Extended Real-Valued Signal Model

3.2. First Stage Rank Reduction DOA Estimation

3.3. Second Stage Rank Reduction DOA Estimation

3.4. Procedures of the Proposed Algorithm

4. Discussion

4.1. Computational Complexity

4.2. Maximum Number of Resolvable Sources

4.3. Dependence of Mutual Coupling on Directions

5. Simulation Results

6. Conclusions

Author Contributions

Funding

Conflicts of Interest

Appendix A

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI