Next Article in Journal
Preparation of “Ion-Imprinting” Difunctional Magnetic Fluorescent Nanohybrid and Its Application to Detect Cadmium Ions
Next Article in Special Issue
An Array Switching Strategy for Direction of Arrival Estimation with Coprime Linear Array in the Presence of Mutual Coupling
Previous Article in Journal
Affiliated Fusion Conditional Random Field for Urban UAV Image Semantic Segmentation
Previous Article in Special Issue
Cumulant-Based DOA Estimation of Noncircular Signals against Unknown Mutual Coupling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

DOA Estimation Using Fourth-Order Cumulants in Nested Arrays with Structured Imperfections

1
National Key Laboratory of Science and Technology on UAV, Northwestern Polytechnical University, Xi’an 710072, China
2
School of Electronics and Information, Northwestern Polytechnical University, Xi’an 710072, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(4), 994; https://doi.org/10.3390/s20040994
Submission received: 24 December 2019 / Revised: 8 February 2020 / Accepted: 10 February 2020 / Published: 12 February 2020

Abstract

:
Recently developed super nested array families have drawn much attention owing to their merits on keeping the benefits of the standard nested arrays while further mitigating coupling in dense subarray portions. In this communication, a new mutual coupling model for nested arrays is constructed. Analyzing the structure of the newly formed mutual coupling matrix, a transformation of the distorted steering vector to separate angular information from the mutual coupling coefficients is revealed. By this property, direction of arrival (DOA) estimates can be determined via a grid search for the minimum of a determinant function of DOA, which is induced by the rank reduction property. We also extend the robust DOA estimation method to accommodate the unknown mutual coupling and gain-phase mismatches in the nested array. Compared with the schemes of super nested array families on reducing the mutual coupling effects, the solutions presented in this paper has two advantages: (a) It is applicable to the standard nested arrays without rearranging the configuration to increase the inter-element spacing, alleviating the cross talk in dense uniform linear arrays (ULAs) as well as gain-phase errors in sparse ULA parts; (b) Perturbations in nested arrays are estimated in colored noise, which is significant but rarely discussed before. Simulations results corroborate the superiority of the proposed methods using fourth-order cumulants.

1. Introduction

Sensor array signal processing is a significant research area owing to its wide applications to radar, sonar, navigation, wireless communications, et al. [1,2,3,4]. Virtues of array signal processing, such as direction of arrival (DOA) estimation and beamforming, lie in spatial diversity and capability to determine the angular information of electromagnetic waves and enhance the reception of signals of interest while mitigating the interferences [5,6,7,8,9]. The basic but widely employed array geometry is the uniform linear array (ULA), where the inter-element spacing is usually less than or equal to half a wavelength to avoid spatial aliasing. However, the performance of ULAs is subject to a limited degrees of freedom (DOFs) that is linearly dependent on to the number of sensors. For instance, an M-elements ULA can offer up to M 1 DOFs, i.e., only M 1 signals can be resolved at most [10,11]. To increase the DOFs, a simple way is to introduce extra sensors, which may result in to difficulties in practical implementations. Besides, ULAs with a large number of sensors may suffer from prominent mutual coupling effects between physical sensors.
Sparse linear arrays (SLAs), a promising countermeasure, can tackle the above issues to a large extent. The critical operation on the DOFs enhancement is to vectorize the covariance matrix of the observations, rendering O ( M 2 ) virtual sensors with M physical elements in the difference coarrays, which means one can probably resolve more DOAs than the physical sensors provided that there is no correlation between signals [12]. Compared with the ULA, this is a big leap in raising DOFs. In addition, the electromagnetic cross talk between sensors is alleviates since the inter-element spacing becomes larger in SLAs [13].
Recently, a new type of sparse array, so-called the nested array [14], is put forward in signal processing society and causes tremendous repercussions because the positions of sensors can be arranged in a regular design, the synthetic apertures becomes larger, and the resultant enhanced DOFs have a closed form expression with respect to the given number of sensors [15,16,17,18,19,20,21,22]. The nested array is constructed by cascading ULAs that differ in systematically designed inter-element spacing. In particular, the difference coarray formed from a two-level nested array has consecutive virtual sensors without any holes, while the hole-free property in the difference coarray cannot be guaranteed by higher level nested arrays that possess higher DOFs. The DOFs of nested arrays can be further increased by making the outer ULA much sparser, i.e., enlarging its aperture, as well as putting up an auxiliary sensor [21]. In [23,24], the one-dimensional nested array is extended to the two-dimensional case.
However, the mutual coupling effects of the nested array cannot be ignored at least as physical sensors in the inner ULA are deployed relatively close [25,26,27,28,29,30,31,32,33,34,35,36,37,38]. In [39], the super nested array in the context of the second-order statistics is designed to significantly mitigate the cross talk between sensors while preserving all advantages of the standard nested arrays, by increase the inter-element spacing of the inner ULA to maintain the coarray but alleviate the adverse electromagnetic effects. In [40], a high order super nested array is introduced to further reduce the mutual coupling between sensors. In [41], by splitting the inner ULA into two or four sections and rearranging them at the two sides of the remaining outer ULA, the so-called augmented nested array (ANA) is developed. Compared with the super nested array, the ANA can further reduce the unknown mutual coupling but at a cost of complicated design to guarantee the difference coarray hole-free. In [42], a new sparse array configuration, namely MISC array, is optimized to obtain maximum inter-element spacing, where a reasonable distribution of three sparse ULAs and two additional sensors is necessary. The MISC array is less susceptive to mutual coupling effects while achieving higher DOFs.
It should be noted that: (a) the aforementioned nested array configurations cannot fundamentally avoid the mutual coupling effects, i.e., all the DOA estimation methods performing on such arrays still work in the presence of perturbations, and (b) the existing nested array processing using the second-order statistics is vulnerable to colored noise. To obtain more robust estimates, we would deal with these two issues from another perspective—consider the array imperfections and jointly resolve DOAs and the parameterized perturbations—rather than reduce the errors by a new geometry design. This idea is reasonable in the sense that the super nested arrays and its successors occupying large areas are not applicable to some circumstances, such as airborne platforms, with a limited space for devices. In this correspondence, a new cumulant-based DOA estimation method for nested arrays with perturbations is proposed. The mutual coupling of the nested array is first analysed and parameterized in terms of the relationship of coupling range and inter-element spacing, giving a new structure of the mutual coupling matrix (MCM). Then we provide a proof that the distorted steering vector can be factorized into a new matrix, including angular information only, multiplied by the coupling coefficient vector. Leveraging this property, we develop a MUSIC-like estimator to determine the DOA estimates based on the rank reduction (RARE) technique in the context fourth-order cumulants (FOC). In addition, the proposed method can be extended to adapt to a more challenging scenario where mutual coupling and gain-phase errors coexist. Analytical specifications show that our solutions are more robust to array imperfections as compared to the standard nested array processing.
The remainder of this paper is organized as follows. In Section 2, an nested array model with unknown mutual coupling is introduced. In Section 3, the MUSIC-like estimators for nested arrays with perturbations using fourth-order cumulants are developed. Section 4 provides numerical examples for demonstrating the validity and efficiency of our proposed algorithms. Finally, some concluding remarks are given in Section 5.

2. Problem Formulation

2.1. Signal Model

Consider N narrowband non-Gaussian signals impinging on a linear M-elements nested array, composed of two different ULAs concatenating together. Assume that the inter-element spacing of the inner ULA with M 1 elements is d I and the inter-element spacing of the outer ULA with M 2 elements is d O satisfying d O = ( M 1 + 1 ) d I . Without loss of generality, we set d I = λ 2 with λ being the wavelength. To be specific, the sensors of a 2-level nested array are deployed at the positions that can be formulated as a union of the position set of inner ULA, S I , and the position set of outer ULA, S O , which can be expressed as
P = S I S O = { m d I | m = 0 , , M 1 1 } { ( n ( M 1 + 1 ) + M 1 ) d I | n = 0 , , M 2 1 }
If there exits mutual coupling between the array elements, the corresponding M × 1 array observation vector is then given by
x ( t ) = i = 1 N C a ( θ i ) s i ( t ) + n ( t ) = C A s ( t ) + n ( t )
where a ( θ i ) = 1 , e j 2 π p 2 λ sin θ i , , e j 2 π p M λ sin θ i T C M is the steering vector, A = a ( θ 1 ) , , a ( θ N ) is the array manifold, C denotes the MCM, and s ( t ) = s 1 ( t ) , , s N ( t ) T . We assume that the incident signals { s i ( t ) } i = 1 N are independent of each other, the noise vector n ( t ) is amenable to the Gaussian distribution, and it is independent of the signals. Besides, we assume that A is unambiguous, i.e., the steering vectors { a ( θ i ) } i = 1 N are linearly independent for any set of distinct { θ i } i = 1 N .
As described in [27,28,29,30,31,32,33,34,35,36,37,38,43], it may reasonably consider that there are remarkable negative correlation between the coupling strength and the inter-element spacing, and the cross talk effects can be neglected if the distance between two sensors is larger than several times the length of the minimum inter-element spacing d I . By assuming that the minimum inter-element spacing of two sensors without mutual coupling to be Q d I , the generalized MCM can be approximated by
C i j = c | p i p j | , | p i p j | Q d I 0 , | p i p j | > Q d I .
For the case of 2-level nested arrays, if Q < M 1 + 1 , the MCM can be specifically expressed as
C = 1 c 1 c P 1 0 c 1 1 c 1 c 2 c 1 1 c 1 c P 1 c 1 1 c P 2 c 1 1 0 1 = blkdiag Toeplitz 1 , c 1 , , c P 1 , 0 1 × ( M 1 + 1 P ) , I ( M 2 1 )
where Toeplitz r denotes a symmetric Toeplitz matrix constructed by the vector r , and 0 < | c 1 | , | c 2 | , , | c P 1 | < c 0 = 1 are the mutual coupling coefficients.

2.2. Parameters Setting

To make the signal model (simulations) much easier to follow and replicate, it should be emphasized herein that throughout this communication, the following parameters are variables in the performance analysis and assessment: signal-to-noise ratio (SNR), the number of snapshots L, and the number of Monte Carlo runs; the fixed parameters adopted in the numerical calculus, deductions, and simulations include: the number of total sensors M, inter-element spacing of the inner ULA d I , the number of sensors the inner ULA M 1 , the number of sensors the outer ULA M 2 , the number of signals N, DOAs θ i i = 1 N , the coupling range P, and the mutual coupling coefficients c i i = 1 P 1 .

3. Proposed FOC-Based DOA Estimator without Mutual Coupling Compensation

3.1. FOC Matrix Construction

Considering the received signals are assumed to be non-Gaussian, one can establish the array FOC matrix between the received data blocks as
B = cum x ( t ) , x H ( t ) , x * ( t ) , x T ( t )
The entry of B in the k 1 1 M + k 2 -th row and the l 1 1 M + l 2 -th column is defined as
B k 1 1 M + k 2 , l 1 1 M + l 2 = cum { x k 1 ( t ) , x l 1 * ( t ) , x k 2 * ( t ) , x l 2 ( t ) } = E [ x k 1 ( t ) x k 2 * ( t ) x l 1 * ( t ) x l 2 ( t ) ] E [ x k 1 ( t ) x k 2 * ( t ) ] E [ x l 1 * ( t ) x l 2 ( t ) ] E [ x k 1 ( t ) x l 1 * ( t ) ] E [ x k 2 * ( t ) x l 2 ( t ) ] E [ x k 1 ( t ) x l 2 ( t ) ] E [ x k 2 * ( t ) x l 1 * ( t ) ]
Substituting (2) into (5), and exploiting a series of properties of FOC in [44], one can further get
B = cum C A s ( t ) , C A s ( t ) H , C A s ( t ) * , C A s ( t ) T = cum p = 1 N C a ( θ p ) s p ( t ) , m = 1 N C a ( θ m ) H s m * ( t ) , q = 1 N C a * ( θ q ) s q * ( t ) , n = 1 N C a ( θ n ) T s n ( t ) = p = 1 N m = 1 N q = 1 N n = 1 N C a ( θ p ) C a ( θ q ) * C a ( θ m ) C a ( θ m ) * H cum s p ( t ) , s m * ( t ) , s q * ( t ) , s n ( t ) = i = 1 N C a ( θ i ) C a ( θ i ) * C a ( θ i ) C a ( θ i ) * H cum s i ( t ) , s i * ( t ) , s i * ( t ) , s i ( t ) = i = 1 N γ i C a ( θ i ) C a ( θ i ) * C a ( θ i ) C a ( θ i ) * H = C A C A * C s C A C A * H
where γ i cum s i ( t ) , s i * ( t ) , s i * ( t ) , s i ( t ) , and C s diag γ 1 , , γ N R N × N .
To resolve the DOA estimates from the FOC built above, some approaches [45,46,47,48] have been proposed. However, these methods are developed for ULAs in the presence of unknown mutual coupling, and cannot be directed applied to the case of nested arrays as the mutual coupling effects of the two array geometries are distinct from each other. By this reason, a new estimator using FOC should be investigated for combating the unknown mutual coupling.

3.2. Robust DOA Estimation Against Unknown Mutual Coupling

By examining the special structure of the MCM of the nested arrays, one finds that the actual steering vector can be factorized in another way by
C a ( θ ) = T ( θ ) c
where c = 1 , c 1 , , c P 1 C P , and T ( θ ) = T 1 + T 2 T 3 C M × P where T 1 , T 2 C M 1 + 1 × P , and T 3 C M 2 1 × P are given by
T 1 p , q = a 1 p + q 1 , p + q M 1 + 2 0 , otherwise
T 2 p , q = a 1 p q + 1 , p q 2 0 , otherwise
T 3 = a 2 , 0 M 2 1 × P 1
with a 1 = F 1 a ( θ ) , F 1 = I M 1 + 1 , 0 ( M 1 + 1 ) × M 2 1 , a 2 = F 2 a ( θ ) , and F 2 = 0 ( M 2 1 ) × M 1 + 1 , I M 2 1 .
Proof of (8).
One can first rewrite the left-hand side of (8) as
C a ( θ ) = C 1 0 0 I M 2 1 a 1 a 2 = C 1 a 1 a 2
where C 1 = Toeplitz 1 , c 1 , , c P 1 , 0 1 × ( M 1 + 1 P ) C M 1 + 1 × M 1 + 1 . By Lemma 3 in [25], C 1 a 1 = T ¯ ( θ ) c where T ¯ ( θ ) = T 1 + T 2 . Considering the first entry of c is 1, one can make an identical transformation between a 2 and c, that is, a 2 = a 2 , 0 M 2 1 × P 1 c = T 3 c . As a result, one has C a ( θ ) = T 1 + T 2 c T 3 c = T 1 + T 2 T 3 c = T ( θ ) c . This completes the proof of (8). □
Now performing the singular value decomposition (SVD) on B , one has
B = U Σ V H
where Σ = diag λ 1 , , λ M 2 consists of M 2 singular values satisfying λ 1 λ N > λ N + 1 = = λ M 2 = 0 . The columns of U s U ( : , 1 : N ) are the singular vectors corresponding to the N largest eigenvalues, while the columns of U n U : , N + 1 : M 2 are the singular vectors corresponding to the M 2 singular values which are all zeros. The signal subspace is spanned by the columns of U s , whereas the noise subspace is spanned by the columns of U n . One can construct the following function
p θ , c = C a ( θ ) C a ( θ ) * H U n 2 2 = T ( θ ) c T ( θ ) c * H U n 2 2 = c c * H T ( θ ) T * ( θ ) H U n U n H T ( θ ) T * ( θ ) c c * = g H c M ( θ ) g c
where g c = c c * C P 2 , and M ( θ ) = T ( θ ) T * ( θ ) H U n U n H T ( θ ) T * ( θ ) C P 2 × P 2 . Next one can estimate DOAs from the determinant M ( θ ) . It can be found that the size of T ( θ ) T * ( θ ) H U n is P 2 × ( M 2 N ) and if P 2 M 2 N , the matrix T ( θ ) T * ( θ ) H U n is normally of full row rank and M ( θ ) is of full rank. However, when θ is equal to any DOA of interest, i.e., θ = θ i , i = 1 , , N , Equation (14) is zero according to the orthogonality between the signal and noise subspaces. Since g c 0 , (14) is available only if the matrix M ( θ ) is rank deficient, i.e., the determinant of M ( θ ) is zero. Hence, one can determine DOA estimates by searching the minimum of det M ( θ ) that is close to zero as
θ ^ = arg min θ det M ( θ ) .

3.3. Mutual Coupling Coefficient Estimation

With the DOA estimates, the orthogonality between the signal and noise subspaces can be exploited to estimate the mutual coupling coefficients. To this end, one can stack a series of equations as
U n H C a ( θ ^ 1 ) C a ( θ ^ 1 ) * U n H C a ( θ ^ N ) C a ( θ ^ N ) * = 0 M 2 N × 1 0 M 2 N × 1 .
Substituting (8) into (16), one gets
U n H T ( θ ^ 1 ) T * ( θ ^ 1 ) U n H T ( θ ^ N ) T * ( θ ^ N ) H c c * = h , H ¯ 1 c ¯ = 0 N M 2 N × 1 .
where h is the first column of H while H ¯ contains the remaining columns of H , and c ¯ = 0 P 2 1 × 1 , I P 2 1 c c * C P 2 1 .
One can carry out the estimation of c ¯ in the least squares sense, that is
c ¯ = H ¯ + h .
Then the mutual coupling coefficients c 1 = c 1 , c 2 , c P 1 can be extracted from c ¯ by
c 1 = F ¯ 1 c ¯ *
c 1 = F ¯ 2 c ¯
where F ¯ 1 = I P 1 , 0 P 1 × P 2 P and F ¯ 2 = e P , e 2 P , , e P ( P 1 ) T with e i R P 2 1 , i = P , 2 P , , P ( P 1 ) , is a column vector with 1 at the i-th entry and 0 elsewhere. Finally, an improved estimate of c 1 can be obtained by averaging the above results as
c 1 = 1 2 F ¯ 1 c ¯ * + F ¯ 2 c ¯ .

3.4. Extension to Partly Calibrated Nested Array with Unknown Mutual Coupling

In this section, the nested array is partly calibrated, i.e., the inner ULA suffers from the mutual coupling effect while the outer ULA has unknown gain and phase errors.
In practice, one may encounter an even more challenging scenario where the mutual coupling and gain-phase mismatch coexist in the nested array. However, it is difficult or even not able to get the issue solved as there are too many unknowns to be handled by the limited DOFs. To make DOA estimation feasible, one has to relieve the nuisance imperfections on the array to some extend, by calibrating a portion of sensors well. This class of arrays, referred to as partly calibrated nested arrays, is common in practice as the aperture of the inner ULA is typically much smaller than that of the whole array and, hence, the inner ULA can be assumed to be well calibrated for gain-phase error while the calibration of the outer ULA may be poor due to unknown channel imbalances between sensors that are located far away from each other. Through above analysis, the MCM is rectified as
C ¯ = blkdiag Toeplitz 1 , c 1 , , c P 1 , 0 1 × ( M 1 + 1 P ) , Υ ( M 2 1 ) × ( M 2 1 )
where Υ = diag β with β = β 1 , β 2 , , β M 2 1 T C M 2 1 containing the element gains β i i = 1 M 2 1 and phases arg β i i = 1 M 2 1 . Following a similar idea, one can derive
C ¯ a ( θ ) = T ¯ ( θ ) α
where α = c T , β T T and T ¯ ( θ ) = blkdiag T 1 + T 2 , T ¯ 3 C M × P + M 2 1 with T ¯ 3 = diag e j 2 π p M 1 + 2 λ sin θ , , e j 2 π p M λ sin θ . To prove (23), one needs to prove Υ a 2 = T ¯ 3 β that can be easily verified by the fact that diag b d = diag d b holds for arbitrary vectors b and d with the same size. Then, by the same principle behind (13)–(15), one can determine the DOA estimates by
θ ^ = arg min θ det M ¯ ( θ ) ,
where M ¯ = T ¯ ( θ ) T ¯ * ( θ ) H U n U n H T ¯ ( θ ) T ¯ * ( θ ) C P + M 2 1 2 × P + M 2 1 2 .
Following the same procedures (16)–(21) and denoting α ¯ = 0 P + M 2 1 2 1 × 1 , I P + M 2 1 2 1 α α * C P + M 2 1 2 1 , one can get the estimate of α ¯ by
α ¯ = H ˜ + h ¯ .
where h ¯ is the first column of H ˘ that is defined as
H ˘ = U n H T ¯ ( θ ^ 1 ) T ¯ * ( θ ^ 1 ) U n H T ¯ ( θ ^ N ) T ¯ * ( θ ^ N ) ,
and H ˜ includes the rest columns of H ˘ . The final mutual coupling coefficient and gain-phase error estimates are determined by α ¯ as
c 1 = 1 2 F ˜ 1 α ¯ * + F ˜ 2 α ¯
β = 1 2 F ˜ 3 α ¯ * + F ˜ 4 α ¯
where F ˜ 1 = I P 1 , 0 P 1 × P + M 2 1 2 P , F ¯ 2 = e ¯ P + M 2 1 , e ¯ 2 P + M 2 1 , , e ¯ ( P 1 ) ( P + M 2 1 ) T , F ˜ 3 = 0 M 2 1 × P 1 , I M 2 1 , 0 M 2 1 × P + M 2 1 2 P M 2 + 1 , and F ¯ 4 = e ¯ P P + M 2 1 , e ¯ P + 1 P + M 2 1 , , e ¯ ( P + M 2 2 ) ( P + M 2 1 ) T with e ¯ i R ( P + M 2 1 ) 2 1 , i = P + M 2 1 , 2 ( P + M 2 1 ) , , ( P + M 2 2 ) ( P + M 2 1 ) , is a column vector with 1 at the i-th entry and 0 elsewhere.
The major steps of the proposed algorithm are summarised in Table 1.
Remark 1.
The methods presented in this paper are valid on condition that the MCM is isolated from any angular information of the incident signals and are not applicable to arbitrary geometries of antenna arrays. However, they still work for some commonly used antenna arrays like dipole, monopole, and slot arrays [49]. On the other hand, the MCM could be direction variant for directional antenna arrays. In such a case, the proposed methods preform well as long as the MCM possesses a similar structure to (4) or (22), even if the mutual coupling coefficients c 1 are direction-dependent. This is because the transformation of product of the nominal steering vector and the direction-independent MCM still holds for the direction-dependent MCM by (8) or (23). To be specific, the array observation vector, of size M × 1 , can be expressed as x ( t ) = C 1 a ( θ 1 ) , C 2 a ( θ 2 , , C N a ( θ N ) ) s ( t ) + n ( t ) where C i is reliant upon θ i , different from each other in general. It should be noted that if C i i = 1 N have the same banded symmetric Toeplitz structure as shown in (4), then one still has C i a ( θ i ) = T ( θ i ) c i . As a result, one can readily deduce that d e t M ( θ ) = 0 holds for θ = θ i , i = 1 , 2 , , N , and the DOA estimates can be obtained by searching the minimum of d e t M ( θ ) . It should be noted that the structure of MCM of directional antenna arrays may not stay the same for all angles in practice, causing an arbitrariness on the composition of M ( θ ) and, hence our solution is inapplicable to arbitrary direction-dependent mutual coupling via searching the spatial spectrum of a unified function d e t M ( θ ) . Only the direction-independent model is discussed herein while its direction-dependent counterpart will be investigated in the future.
Remark 2.
As with the previous work [49], the coupling range P in (4) is assumed to be known a priori, otherwise it has to be determined properly. If P is selected to be smaller than the true value, the remaining mutual coupling may still deteriorate the estimation performance, whereas if the chosen P is larger than the true one, the accuracy and resolution of estimates are confined due to the loss of array aperture. A pragmatic approach to identify P is to measure the minimal distance between two antennas where the coupling effect is negligible. In practice, the coupling range depends on various factors of antenna design, such as array geometry, ground structure, material characteristic, etc, and the effect of mutual coupling is weak enough to be neglected in many applications when the inter-element spacing of an antenna array is larger than half a wavelength, which means that P is relatively small.
Remark 3.
It should be noted that one may get some pseudo estimates in addition to the true DOAs, satisfying (15) or (24). Without loss of generality, we assume that there are N 0 N DOA estimates, pseudo and true, one obtains by a spectral search of (15) or (24). Substituting the DOA estimates into (21) or (27), on can determine N 0 mutual coupling coefficients vectors and reconstruct the corresponding MCM C ^ i i = 1 N 0 . Defining a spatial spectrum equation as follows
f = C ^ i a ( θ ^ i ) C ^ i a ( θ ^ i ) * H U n 2 2 , i = 1 , 2 , , N 0 ,
one can sift the estimates of true DOAs and mutual coupling coefficients from the pseudo ones by identifying whether f = 0 or not as the orthogonality between the signal subspace C ^ i a ( θ ^ i ) C ^ i a ( θ ^ i ) * and the noise subspace U n does not hold for the false DOAs and mutual coupling coefficients. In practice, only finite snapshots are available, so the true DOA and mutual coupling coefficient estimates correspond to the minima of (29) that returns large values for the rest N 0 N pseudo estimates. By this means, one can tell the true estimates from the false ones.
Remark 4.
It is known that T ( θ ) T * ( θ ) H U n C P 2 × ( M 2 N ) and T ¯ ( θ ) T ¯ * ( θ ) H U n C ( P + M 2 1 ) 2 × ( M 2 N ) , so the proposed method works on condition that P 2 M 2 N for the case of unknown mutual coupling only and ( P + M 2 1 ) 2 M 2 N for the scenario where unknown mutual coupling and gain-phase error occur on dense ULAs and sparse ULA parts, respectively. This provides upper bounds on the number of resolvable signals for these two circumstances: one is N M 2 P 2 , and the other is N M 2 ( P + M 2 1 ) 2 , which means that the developed algorithm can handle more signals in the former than in the latter, or equivalently, more DOFs and effective aperture are available when mutual coupling takes place only. Therefore, the proposed Method 1 has better performance than Method 2 in that intuitively, the more the unknowns are, the worse the estimation arises, while a reverse phenomenon can be observed since Method 1 does not take the gain-phase error into consideration, causing a considerable model mismatch.

4. Simulation Results and Discussion

In this section, various numerical experiments under different conditions are carried out to examine the performance of the proposed algorithms. Simulations are conducted for 2-level nested array with four elements in each level. For simplicity, we assume that all signals have an identical power. Similar to the settings in [50,51], the signals are modeled as s ( t ) = F ( t ) r ( t ) , where F ( t ) = diag f 1 ( t ) , , f K ( t ) , r ( t ) = r 1 ( t ) , , r K ( t ) T with f i ( t ) and r i ( t ) being zero-mean Gaussian processes with unit-variance and σ s 2 -variance, respectively. The noise is assumed to be a first order spatial autoregressive process, and the ( a , b ) -th entry of the noise covariance matrix is given by R ( a , b ) = σ n 2 0.8 | p a p b λ | [52,53]. The SNR is defined as SNR = 10 log 10 ( σ s 2 / σ n 2 ) . The accuracy of the DOA estimate, the statistical performance of the algorithms, is measured from 800 Monte Carlo runs in terms of the root mean square error (RMSE) which is defined as
RMSE θ = 1 800 N n = 1 800 i = 1 N ( θ ^ i ( n ) θ i ) 2
RMSE c = 1 800 c 1 2 2 n = 1 800 c ^ 1 ( n ) c 1 2 2 × 100
where θ ^ i ( n ) and c ^ 1 ( n ) are the estimates of θ i and c 1 , respectively, for the n-th trial, and N is the number of signals.
In the first scenario, we consider that two independent sources from 17 , 43 impinge on the nested array. The resultant RMSEs of DOA estimates are illustrated in Figure 1. It can be seen that the proposed Method 1 approaches the FOC-based estimator with known coupling asymptotically and significantly outperforms the Method 2 with the increase of SNR. The main reason is that Method 1 deals with the outer ULA as calibrated and, hence more DOFs are obtained for DOA estimation, whereas Method 2 considers a even worse case, where the outer ULA has gain-phase errors, and the corresponding DOFs are occupied by the nuisance parameters. Besides, the estimation error of Method 2 decreases clearly as the number of snapshots increases and stabilize at approximately 0.35 with the observation size larger than 500 snapshots, while both Method 1 and 4-MUSIC provide much better accuracy, saturating at RMSEs of 0.05 and 0.03 , respectively, than Method 1 through all snapshot sizes when SNR is fixed at 5 dB.
Table 2 and Table 3 list the RMSEs of the mutual coupling coefficient estimates under different SNRs and snapshot sizes. Similar to the performance of DOA estimation in the first scenario, the above tabular results show that the proposed Method 1 is able to offer satisfactory estimation of mutual coupling, especially for moderate SNRs and the numbers of snapshots. A plausible explanation is that compared with Method 2, using the calibrated subarray has implicit superiorities, such as a lower Cramer-Rao lower bound for mutual coupling estimation, since its unknowns are less than those with the partly calibrated array. As a result, Method 1 can achieve better performance on mutual coupling calibration.

5. Conclusions

In this paper, an FOC-based DOA estimation method for nested arrays is proposed to simultaneously improve the robustness of direction finding of non-Gaussian signals to array imperfections and mitigate the spatially colored noise. Under the assumptions that the mutual coupling effects are direction-independent and its range is known a priori, the mutual coupling effects of nested arrays is analytically specified from the aspect of the relationship between coupling strength and inter-element spacing, and a new factorization of the distorted steering vector into a matrix, containing bearing information only multiplied by the coupling coefficient vector. Combining this property with RARE technique, a MUSIC-like estimator is developed to obtain the DOA estimates in FOC. Additionally, the proposed estimator can be extended to handle a much harsher issue that mutual coupling and gain-phase mismatches separately occur on dense and sparse sunbarrays, respectively. Analytical specifications show that our solutions are more robust to array imperfections as compared to the standard nested array processing. Compared with the schemes of super nested array families on reducing the mutual coupling effects, our solutions enjoy two merits: (a) It is applicable to the standard nested arrays without rearranging the configuration to increase the inter-element spacing, isolating the cross talk in dense subarrays as well as gain-phase errors in sparse ULA parts; (b) Perturbations in nested arrays are estimated in colored noise, which has not been fully addressed. Simulations results corroborate that the proposed methods are advantageous to self array calibration as well as immunity to colored noise.

Author Contributions

All authors contributed extensively to the study presented in this manuscript. B.W. designed the main idea, methods and experiments, interpreted the results and wrote the paper. J.Z. carried out the experiments, edited the manuscript and provided many valuable suggestions to this study. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grant 61472324, and in part by the Seed Foundation of Innovation and Creation for Graduate Students in Northwestern Polytechnical University under Grant ZZ2019023.

Acknowledgments

The authors would like to thank Chow Yii Pui from the University of Adelaide for his many constructive suggestions and comments that helped to improve the quality of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Krim, H.; Viberg, M. Two decades of array signal processing research: The parametric approach. IEEE Signal Process. Mag. 1996, 13, 67–94. [Google Scholar] [CrossRef]
  2. Van Trees, H.L. Detection Estimation and Modulation Theory, Part IV, Optimum Array Processing; John Wiley and Sons: New York, NY, USA, 2002. [Google Scholar]
  3. Yan, J.; Pu, W.; Zhou, S.; Liu, H.; Bao, Z. Collaborative detection and power allocation framework for target tracking in multiple radar system. Inform. Fusion 2020, 55, 173–183. [Google Scholar] [CrossRef]
  4. Shi, C.; Wang, F.; Sellathurai, M.; Zhou, J.; Salous, S. Low probability of intercept-based optimal power allocation scheme for an integrated multistatic radar and communication system. IEEE Syst. J. 2019. [Google Scholar] [CrossRef] [Green Version]
  5. Shi, Z.; Zhou, C.; Gu, Y.; Goodman, N.A.; Qu, F. Source estimation using coprime array: A sparse reconstruction perspective. IEEE Sensors J. 2017, 17, 573–586. [Google Scholar] [CrossRef]
  6. 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]
  7. 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]
  8. 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]
  9. Wang, Y.; Trinkle, M.; Ng, B.W.H. Efficient DOA estimation of noncircular signals in the presence of multipath propagation. Signal Process. 2018, 149, 14–26. [Google Scholar] [CrossRef]
  10. Schmidt, R.O. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef] [Green Version]
  11. Roy, R.; Kailath, T. ESPRIT-Estimation of signal parameters via rotational invariance techniques. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 984–995. [Google Scholar] [CrossRef] [Green Version]
  12. Hoctor, R.T.; Kassam, S.A. The unifying role of the coarray in aperture synthesis for coherent and incoherent imaging. Proc. IEEE 1990, 78, 735–752. [Google Scholar] [CrossRef]
  13. BouDaher, E.; Ahmad, F.; Amin, M.G.; Hoorfar, A. Mutual coupling effect and compensation in non-uniform arrays for direction-of-arrival estimation. Digit. Signal Process. 2017, 61, 3–14. [Google Scholar] [CrossRef] [Green Version]
  14. Pal, P.; Vaidyanathan, P.P. Nested arrays: A novel approach to array processing with enhanced degrees of freedom. IEEE Trans. Signal Process. 2010, 58, 4167–4181. [Google Scholar] [CrossRef] [Green Version]
  15. Pal, P.; Vaidyanathan, P.P. Multiple level nested array: An efficient geometry for 2qth order cumulant based array processing. IEEE Trans. Signal Process. 2012, 60, 1253–1269. [Google Scholar] [CrossRef]
  16. Han, K.; Nehorai, A. Wideband Gaussian source processing using a linear nested array. IEEE Signal Process. Lett. 2013, 20, 1110–1113. [Google Scholar]
  17. Han, K.; Nehorai, A. Improved source number detection and direction estimation with nested arrays and ULAs using jackknifing. IEEE Trans. Signal Process. 2013, 61, 6118–6128. [Google Scholar] [CrossRef]
  18. Han, K.; Nehorai, A. Nested array processing for distributed sources. IEEE Signal Process. Lett. 2014, 21, 1111–1114. [Google Scholar]
  19. Han, K.; Yang, P.; Nehorai, A. Calibrating nested sensor arrays with model errors. IEEE Trans. Antennas Propag. 2015, 63, 4739–4748. [Google Scholar] [CrossRef]
  20. Liu, C.L.; Vaidyanathan, P.P. Remarks on the spatial smoothing step in coarray MUSIC. IEEE Signal Process. Lett. 2015, 22, 1438–1442. [Google Scholar] [CrossRef]
  21. Yang, J.; Liao, G.; Li, J. An efficient off-grid DOA estimation approach for nested array signal processing by using sparse Bayesian learning strategies. IEEE Signal Process. Lett. 2016, 128, 110–122. [Google Scholar] [CrossRef]
  22. Wang, Y.; Hashemi-Sakhtsari, A.; Trinkle, M.; Ng, B.W.H. Sparsity-aware DOA estimation of quasi-stationary signals using nested arrays. Signal Process. 2018, 144, 87–98. [Google Scholar] [CrossRef]
  23. Pal, P.; Vaidyanathan, P.P. Nested arrays in two dimensions, Part I: Geometrical considerations. IEEE Trans. Signal Process. 2012, 60, 4694–4705. [Google Scholar] [CrossRef]
  24. Pal, P.; Vaidyanathan, P.P. Nested arrays in two dimensions, Part II: Application in two dimensional array processing. IEEE Trans. Signal Process. 2012, 60, 4706–4718. [Google Scholar] [CrossRef]
  25. Friedlander, B.; Weiss, A.J. Direction finding in the presence of mutual coupling. IEEE Trans. Antennas Propag. 1991, 39, 273–284. [Google Scholar] [CrossRef]
  26. Sellone, F.; Serra, A. A novel online mutual coupling compensation algorithm for uniform and linear arrays. IEEE Trans. Antennas Propag. 2007, 55, 560–573. [Google Scholar] [CrossRef]
  27. Ye, Z.; Liu, C. On the resiliency of MUSIC direction finding against antenna sensor coupling. IEEE Trans. Antennas Propag. 2008, 56, 371–380. [Google Scholar] [CrossRef]
  28. Xu, X.; Ye, Z.F.; Zhang, Y.F. DOA estimation for mixed signals in the presence of mutual coupling. IEEE Trans. Signal Process. 2009, 57, 3523–3532. [Google Scholar] [CrossRef]
  29. Dai, J.; Xu, W.; Zhao, D. Real-valued DOA estimation for uniform linear array with unknown mutual coupling. Signal Process. 2012, 92, 2056–2065. [Google Scholar] [CrossRef]
  30. Dai, J.; Zhao, D.; Ji, X. A sparse representation method for DOA estimation with unknown mutual coupling. IEEE Antennas Wireless Propag. Lett. 2012, 11, 1210–1213. [Google Scholar] [CrossRef]
  31. Dai, J.; Xu, X.; Zhao, D. Direction-of-arrival estimation via real-valued sparse representation. IEEE Antennas Wireless Propag. Lett. 2013, 12, 376–379. [Google Scholar] [CrossRef]
  32. Liu, Z.; Zhou, Y. A unified framework and sparse Bayesian perspective for direction-of-arrival estimation in the presence of array imperfections. IEEE Trans. Signal Process. 2013, 61, 3786–3798. [Google Scholar] [CrossRef]
  33. 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] [PubMed]
  34. Liu, H.; Zhao, L.; Li, Y.; Jing, X.; Truong, T.K. A sparse-based approach for DOA estimation and array calibration in uniform linear array. IEEE Sensors J. 2016, 16, 6018–6027. [Google Scholar] [CrossRef]
  35. Wang, Q.; Dou, T.; Chen, H.; Yan, W.; Liu, W. Effective block sparse representation algorithm for DOA estimation with unknown mutual coupling. IEEE Commun. Lett. 2017, 21, 2622–2625. [Google Scholar] [CrossRef]
  36. Hawes, M.; Mihaylova, L.; Septier, F.; Godsill, S. Bayesian compressive sensing approaches for direction of arrival estimation with mutual coupling effects. IEEE Trans. Antennas Propag. 2017, 65, 1357–1368. [Google Scholar] [CrossRef] [Green Version]
  37. Wang, X.; Meng, D.; Huang, M.; Wan, L. Reweighted regularized sparse recovery for DOA estimation with unknown mutual coupling. IEEE Commun. Lett. 2019, 23, 290–293. [Google Scholar] [CrossRef]
  38. 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]
  39. Liu, C.L.; Vaidyanathan, P.P. Super nested arrays: Linear sparse arrays with reduced mutual coupling— Part I: Fundamentals. IEEE Trans. Signal Process. 2016, 64, 3997–4012. [Google Scholar] [CrossRef]
  40. Liu, C.L.; Vaidyanathan, P.P. Super nested arrays: Linear sparse arrays with reduced mutual coupling— Part II: High-order extensions. IEEE Trans. Signal Process. 2016, 64, 4203–4217. [Google Scholar] [CrossRef]
  41. Liu, J.; Zhang, Y.; Lu, Y.; Ren, S.; Cao, S. Augmented nested arrays with enhanced DOF and reduced mutual coupling. IEEE Trans Signal Process. 2017, 65, 5549–5563. [Google Scholar] [CrossRef]
  42. Zheng, Z.; Wang, W.; Kong, Y.; Zhang, Y.D. MISC Array: A new sparse array design achieving increased degrees of freedom and reduced mutual coupling dffect. IEEE Trans. Signal Process. 2019, 67, 1728–1741. [Google Scholar] [CrossRef]
  43. 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]
  44. Dǒgan, M.C.; Mendel, J.M. Applications of cumulants to array processing—Part I: Aperture extension and array calibration. IEEE Trans. Signal Process. 1995, 43, 1200–1216. [Google Scholar] [CrossRef] [Green Version]
  45. Liu, C.; Ye, Z.; Zhang, Y. DOA estimation based on fourth-order cumulants with unknown mutual coupling. Signal Process. 2009, 89, 1839–1843. [Google Scholar] [CrossRef]
  46. 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]
  47. Liao, B.; Chan, S.C. A cumulant-based method for direction finding in uniform linear arrays with mutual coupling. IEEE Antennas Wireless Propag. Lett. 2014, 13, 1717–1720. [Google Scholar] [CrossRef]
  48. Ahmed, A.; Tufail, M. Genetic algorithm-based improved DOA estimation using fourth-order cumulants. Int. J. Electron. 2017, 104, 747–760. [Google Scholar] [CrossRef]
  49. Craeye, C.; González-Ovejero, D. A review on array mutual coupling analysis. Radio Sci. 2011, 46, 1–25. [Google Scholar] [CrossRef] [Green Version]
  50. Ye, Z.; Zhang, Y. DOA estimation for non-Gaussian signals using fourth-order cumulants. IET Microw. Antennas Propag. 2009, 3, 1069–1078. [Google Scholar] [CrossRef]
  51. Wang, Y.; Trinkle, M.; Brian, B.W.H. Two-stage DOA estimation of independent and coherent signals in spatially coloured noise. Signal Process. 2016, 128, 350–359. [Google Scholar] [CrossRef]
  52. Wong, K.M.; Reilly, J.P.; Wu, Q.; Qiao, S. Estimation of the directions of arrival of signals in unknown correlated noise, part I: The MAP approach and its implementation. IEEE Trans. Signal Process. 1992, 40, 2007–2017. [Google Scholar] [CrossRef]
  53. Viberg, M.; Stoica, P.; Ottersten, B. Maximum likelihood array processing in spatially correlated noise fields using parameterized signals. IEEE Trans. Signal Process. 1997, 45, 996–1004. [Google Scholar] [CrossRef]
Figure 1. RMSE of the DOA estimates of 2 independent signals versus (a) SNR when N = 500 ; (b) N when SNR = 5 dB.
Figure 1. RMSE of the DOA estimates of 2 independent signals versus (a) SNR when N = 500 ; (b) N when SNR = 5 dB.
Sensors 20 00994 g001
Table 1. Summary of the proposed method against the perturbations of nested arrays.
Table 1. Summary of the proposed method against the perturbations of nested arrays.
Step 1.Compute the cumulant matrix B from observations x ( t ) according to (5).
Step 2.Extract the noise subspace U n by performing SVD on B.
Step 3.Construct T θ (or T ¯ θ ) composed by T 1 in (9), T 2 in (22),
and T 3 in (11) (or T ¯ 3 in (23)).
Step 4.Determine the DOA estimates θ ^ i i = 1 N by searching for N
minima of d e t M ( θ ) defined in (14) (or M ¯ ( θ ) in (24)).
Step 5.Reconstruct H with the DOA estimates θ ^ i in (17) (or H ˘ in (26)).
Step 6.Estimate c ¯ in (18) (or α ¯ in (25)).
Step 7.Obtain estimates of the mutual coupling coefficients c ^ 1 in (21)
(or c ^ 1 in (27) and estimates of gain-phase errors β ^ in (28)).
Table 2. Root mean square error (RMSE) of c 1 versus SNR when L = 500 .
Table 2. Root mean square error (RMSE) of c 1 versus SNR when L = 500 .
Method 1Method 2
SNRRMSE c RMSE c
5 dB8.534%35.94%
7 dB5.621%23.413%
9 dB3.7062%13.881%
11 dB2.5199%8.2154%
13 dB1.7612%4.9845%
15 dB1.2811%3.2019%
Table 3. RMSE of c 1 versus the number of snapshots when SNR = 5 dB.
Table 3. RMSE of c 1 versus the number of snapshots when SNR = 5 dB.
Method 1Method 2
The NumberRMSE c RMSE c
of Snapshots
10008.1907%35.118%
12008.0853%34.504%
14008.0204%34.884%
16008.0126%34.465%
18007.9647%34.305%
20007.9368%34.043%

Share and Cite

MDPI and ACS Style

Wang, B.; Zheng, J. DOA Estimation Using Fourth-Order Cumulants in Nested Arrays with Structured Imperfections. Sensors 2020, 20, 994. https://doi.org/10.3390/s20040994

AMA Style

Wang B, Zheng J. DOA Estimation Using Fourth-Order Cumulants in Nested Arrays with Structured Imperfections. Sensors. 2020; 20(4):994. https://doi.org/10.3390/s20040994

Chicago/Turabian Style

Wang, Baoping, and Junhao Zheng. 2020. "DOA Estimation Using Fourth-Order Cumulants in Nested Arrays with Structured Imperfections" Sensors 20, no. 4: 994. https://doi.org/10.3390/s20040994

APA Style

Wang, B., & Zheng, J. (2020). DOA Estimation Using Fourth-Order Cumulants in Nested Arrays with Structured Imperfections. Sensors, 20(4), 994. https://doi.org/10.3390/s20040994

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

Article Metrics

Back to TopTop