Next Article in Journal
Waste Material Classification: A Short-Wave Infrared Discrete-Light-Source Approach Based on Light-Emitting Diodes
Previous Article in Journal
Validity of the CALERA Research Sensor to Assess Body Core Temperature during Maximum Exercise in Patients with Heart Failure
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

A Novel Sparse Array Configuration for Direction of Arrival Estimation with Increased Uniform Degrees of Freedom and Reduced Mutual Coupling

College of Information, Mechanical, and Electrical Engineering, Shanghai Normal University, Shanghai 200234, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(3), 808; https://doi.org/10.3390/s24030808
Submission received: 26 December 2023 / Revised: 23 January 2024 / Accepted: 23 January 2024 / Published: 26 January 2024
(This article belongs to the Section Sensor Networks)

Abstract

:
Sparse arrays are widely employed in array signal processing due to their obvious advantages in array element distribution and uniform degrees of freedom (uDOFs). In this paper, a generalized augmented multi-subarray nested array (GAMSNA-I) and its variant, GAMSNA-II are proposed, with the objective of increasing uDOFs and reducing mutual coupling. Based on two subarrays of the prototype nested array (NA), GAMSNA-I is constructed by reconfiguring the dense uniform linear array (ULA) and forward-shifting the sparse ULA. GAMSNA-II is obtained by sparsifying the dense part of GAMSNA-I, ensuring constant uDOFs while further reducing mutual coupling. Subsequently, the closed-form expression for the uDOFs of GAMSNA-I with an arbitrary number of sensors is derived, and the proof is provided that the uDOFs of GAMSNA-II remain unchanged relative to that of GAMSNA-I. Compared to some existing array configurations, both GAMSNA-I and GAMSNA-II exhibit improved uDOFs, with GAMSNA-II achieving lower mutual coupling. Simulation results show the superior performance of the proposed GAMSNA-I and GAMSNA-II.

1. Introduction

Direction of arrival (DOA) estimation is a critical signal-processing technique extensively applied in fields such as radar, sonar, and wireless communication [1,2,3,4,5] to determine the relative direction of signal sources. Uniform linear arrays (ULAs) [6,7,8,9] are frequently used sensor array types for DOA estimation, as they have inter-element spacing smaller than half-wavelength, which helps avoid angular ambiguities. However, the uniform degrees of freedom (uDOFs) and array aperture of ULAs are limited, impacting the ability to resolve the number of sources and the accuracy of DOA estimation significantly. Therefore, in order to ensure estimation performance, large-sized ULAs are often employed. Unfortunately, this method not only increases system complexity and hardware costs but also results in significant array mutual coupling effects.
Alternatively, the emergence of sparse arrays [10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25] provides a direction for solving the problem of ULAs. Via second-order differencing transformation, sparse arrays can obtain virtual ULAs with functionality similar to physical ULAs [10]. Specifically, a sparse array with M physical sensors can provide M 2 virtual sensors, theoretically allowing the identification of sources up to M 2 1 in number. Additionally, the large inter-element spacing in sparse arrays significantly mitigates array mutual coupling effects.
Typical sparse arrays include the minimum redundant array (MRA) [11], the nested array (NA) [12], and the coprime array (CPA) [13]. MRA features an optimal array geometry, as it offers the maximum contiguous lags for a given number of sensors. However, it lacks closed-form expressions for sensor positions and uDOFs, which necessitates exhaustive searches over all possible spatial combinations of available sensors. Motivated by this, NA and CPA were subsequently proposed. Both NA and CPA can be seen as combinations of two ULAs expanded in different ways. Specifically, one ULA remains fixed, while the inter-element spacing of another ULA is expanded to the number of elements in the first ULA plus one, resulting in the formation of NA. CPA is composed of two sparse ULAs, one consisting of M sensors with an inter-element spacing of N units, and the other consisting of N elements with an inter-element spacing of M units, where M and N are coprime. The sparser nature of CPA results in smaller mutual coupling effects; however, it exhibits a significant number of holes in its difference coarray (DCA), which prevents the full utilization of its virtual elements. In contrast, NA has a completely contiguous lags in its DCA, but the significant mutual coupling effect limits its application to a certain extent. Subsequent research predominantly concentrated on how to maintain the advantages of these two arrays while addressing their corresponding shortcomings. For instance, coprime array with displaced subarrays (CADiS) [14], enhanced CACIS [15], thinned coprime array (TCPA) [16], padded coprime array (PCPA) [17], and k times extended coprime array (kECPA) [18] have partially filled the holes in the DCA of CPA through various methods, thereby increasing its uDOFs. However, these strategies still do not achieve a hole-free DCA. Additionally, the generalized nested array (GNA) [19], enhanced nested array (ENA) [20], improved nested array (INA) [21], Iizuka NA [22], super nested array (SuperNA) [23,24], augmented nested array (ANA) [25], and nested array with three sub-ULAs (NA-TS) [26] reconfigured the array geometry of NA to reduce mutual coupling effects and enhance their uDOFs. However, ENA simply divides a dense ULA into two smaller dense ULAs, so its mutual coupling effect is hardly reduced. Meanwhile, INA and the Iizuka NA retain the dense subarray of the prototype NA and suffer from severe mutual coupling effects. Although SuperNA and GNA significantly reduce mutual coupling, they do not increase its uDOFs. In addition, the DCAs of GNA and ANA both have holes. NA-TS shows some improvement in both metrics, but the improvement is relatively modest.
In this paper, a generalized augmented multi-subarray nested array (GAMSNA-I) and its variant, GAMSNA-II are proposed. GAMSNA-I is obtained by decomposing and reconfiguring the dense ULA of NA into multiple parts, resulting in an increase in uDOFs and a reduction in mutual coupling effects. Considering the existence of a dense ULA within GAMSNA-I, further decomposition and reconfiguration of this part results in GAMSNA-II, which maintains the same uDOFs while significantly reducing mutual coupling effects. Additionally, closed-form expressions for the array configurations and uDOFs of both GAMSNA-I and GAMSNA-II are derived and analyzed. Simulation results are provided to demonstrate their superiority in DOA estimation performance.
The rest of this paper is organized as follows: Section 2 introduces the signal model for sparse arrays, relevant concepts of the DCA, and the DOA estimation process. The process of array configuration for the proposed GAMSNA-I and GAMSNA-II is presented in Section 3. The superior performance of these arrays is demonstrated through simulation results in Section 4. Section 5 provides the conclusion of the paper.

2. Preliminaries

2.1. Signal Model

Consider a sparse array with M sensors, where the set of sensor positions is S = { l i d | l i L , i = 1 , 2 , , M } . Here, d is the unit spacing, usually expressed as λ / 2 , where λ is the wavelength of the signal. L is a set of integers arranged in ascending order, with the first element serving as a reference point and assigned a value of 0, i.e., l 1 = 0 . Assume that K is far-field uncorrelated narrowband source signals imping on the array, with the DOAs of Θ = { θ k | k = 1 , 2 , , K } , where θ k ( π / 2 , π / 2 ) . Then, at time t, the received signal vector of the array is modeled as
x ( t ) = As ( t ) + n ( t ) = k K a ( θ k ) s k ( t ) + n ( t )
where A = [ a ( θ 1 ) , a ( θ 2 ) , , a ( θ K ) ] represents the array manifold matrix with the column a ( θ k ) = [ 1 , e j 2 π l 2 d sin θ k / λ , , e j 2 π l M d sin θ k / λ ] T . s ( t ) = [ s 1 ( t ) , s 2 ( t ) , , s K ( t ) ] T is the source vector, and the power of the k-th source is p k 2 . n ( t ) denotes zero-mean complex additive Gaussian white noise with variance of σ 2 I M , where I M represents the M × M -dimensional identity matrix.

2.2. Difference Coarray

The essence of DCA lies in the utilization of second-order statistics, namely the covariance matrix of the received signal [27]. According to (1), the covariance matrix of x ( t ) is calculated as
R xx = E [ x x H ] = A R ss A H + σ 2 I M 1 N n = 1 N x x H
where R ss = d i a g ( [ p 1 2 , p 2 2 , , p K 2 ] ) is the signal covariance matrix, and N represents the number of snapshots.
By vectorizing (2), it can be obtained that
Z = v e c ( R xx ) = B G + σ 2 W M
where B = A * A can be seen as an array manifold matrix, G = d i a g ( R ss ) = [ p 1 2 , p 2 2 , , p K 2 ] T and ⊙ is the Khatri-Rao product. Meanwhile, W M = v e c ( I M ) .
Comparing (1) and (3), it can be revealed that (3) is a received signal vector of the array under a single snapshot condition. In this case, B denotes the array manifold matrix, which is composed of the steering vectors corresponding to the DCA of L . The DCA of L is given by
D = l i l j i , j = 1 , 2 , , M
Definition 1. 
(Degrees of freedom, DOFs): Given a sparse array L , the DOFs is the cardinality of its DCA D .
Definition 2. 
(Uniform degrees of freedom, uDOFs): The central continuous part of D is named as U , the uDOFs is defined as the cardinality of U .
Definition 3. 
(Virtual element weight function): The number of repetitions of the element with a value of u in D is represented as the weight function w ( u ) for that element.
Therefore, w ( u ) can be written as
w ( u ) = { l i , l j L | l i l j = u , u D }

2.3. DOA Estimation

By analyzing (3), it can be seen that Z not only includes all the received data corresponding to the virtual ULA but also contains redundant elements. Furthermore, according to Definition 2, there exists a possibility that U is not equal to D ; in other words, the DCA may not be fully continuous lags. Therefore, to obtain the covariance data corresponding to the virtual ULA, based on Definitions 2 and 3, Z is subjected to filtering, sorting, and redundancy removal. The processing result of Z is denoted as Z ¯ , and it is evident that its dimension is Q × 1 , where Q stands for the uDOFs.
To construct the covariance matrix of the received signals corresponding to the virtual ULA, the spatial smoothing technique [28] is employed to perform the Toeplitz transformation on Z ¯ , and the result is as follows:
R v = Z ¯ Q + 1 2 Z ¯ Q + 1 2 + 1 Z ¯ Q Z ¯ Q + 1 2 1 Z ¯ Q + 1 2 Z ¯ Q 1 Z ¯ 1 Z ¯ 2 Z ¯ Q + 1 2
Subsequently, some covariance matrix-based algorithms can be directly applied to R v for DOA estimation, such as the MUSIC algorithm.

2.4. Mutual Coupling

The scenario corresponding to (1) does not consider the mutual coupling effects between physical sensors. However, in practice, the mutual coupling effects between sensors are one of the important factors affecting estimation results and cannot be directly ignored. In the context of the presence of mutual coupling effects, (1) is rewritten as
x ( t ) = CAs ( t ) + n ( t )
where C is the mutual coupling matrix with a dimension of M × M . According to [29], in the case of a linear array, the expression for C is simplified without the need to consider the influence of various factors such as humidity, operating frequency, adjacent objects, etc. Therefore, C can be approximated as a B-banded symmetric matrix as
C l i , l j = c l i l j , i f l i l j B 0 , e l s e w h e r e
where l i , l j L , and c i , i = 0 , 1 , , B are the coupling coefficients, satisfying c 0 = 1 > | c 1 | > | c 2 | > > | c B | and c i / c j = j / i , i , j = 0 , 1 , , B .
To objectively evaluate the mutual coupling effects, the coupling leakage is introduced, which is represented as
L = | | C d i a g ( C ) | | F | | C | | F
where | | · | | F is the Frobenius norm of a matrix. L reflects the energy proportion of the off-diagonal components of C , and a smaller value indicates weaker mutual coupling effects between elements in the sensor array.

3. Proposed Array Configurations

In this section, the structural compositions of the proposed GAMSNA-I and GAMSNA-II are primarily introduced, along with their closed-form expressions of sensor positions and achievable uDOFs. By reconfiguring two subarray elements of the prototype NA, GAMSNA-I is proposed, increasing uDOFs while simultaneously reducing mutual coupling influences. Furthermore, to further mitigate the impact of mutual coupling effects while maintaining uDOFs, improvements are made to GAMSNA-I to obtain GAMSNA-II.

3.1. GAMSNA-I

The normalized set of element positions of the prototype NA, formed by two subarrays belonging to a dense ULA and a sparse ULA respectively, is denoted as
L N A = { m 1 , 0 m 1 M 1 1 } { m 2 ( M 1 + 1 ) + M 1 , 0 m 2 M 2 1 }
where M 1 and M 2 represent the number of sensors in two subarrays, respectively, which satisfies the optimal configuration of NA, as shown in Table 1.
The mutual coupling between the elements of NA primarily arises from the dense ULA portion. Therefore, the main idea of GAMSNA-I is to reconfigure the structure of the dense ULA and make appropriate adjustments to the sparse ULA, achieving an increase in uDOFs while reducing mutual coupling.
The structure of GAMSNA-I consists of four parts, represented as
L G A M S N A I = L 1 L 2 L 3 L 4
where the closed-form expressions of L 1 , L 2 , L 3 , and L 4 are
L 1 = { 0 , 1 , M 1 } , M 1 4 L 2 = { 3 + ( M 1 + 1 ) | = 0 , 1 , , M 2 1 } L 3 = { M 2 ( M 1 + 1 ) M 1 + 3 } L 4 = { ( M 1 + 1 ) ( M 2 + 1 ) M 1 + 2 η | η = 1 , 2 , , M 1 4 }
Note in (12) that M 1 4 . Combining this with Table 1, it can be deduced that the minimum number of sensors in GAMSNA-I is 8, i.e., M 8 . Additionally, when M 1 = 4 , i.e., M = 8 or M = 9 , L 4 = 0 , the array aperture of GAMSNA-I is equal to that of NA, and both arrays have a DCA without holes, resulting in equal uDOFs for both. This implies that the advantage of GAMSNA-I in this case is the relatively less mutual coupling. When M 1 = 5 , i.e., M = 10 or M = 11 , L 4 = 1 , although the array aperture of GAMSNA-I is larger than that of NA, there are holes in the DCA of GAMSNA-I, and the uDOFs of GAMSNA-I are less than those of NA. In this case, GAMSNA-I sacrifices two virtual array elements to reduce mutual coupling.
When M 12 , GAMSNA-I demonstrates advantages both in terms of uDOFs and mutual coupling. The array aperture of GAMSNA-I is larger than that of NA, and its DCA is hole-free, as summarized in Proposition 1.
Proposition 1. 
The DCA of the GAMSNA-I as defined in (12) is a hole-free ULA, i.e., D = U = [ Q g , Q g ] with Q g = ( M 2 + 1 ) M 2 + 2 .
Proof of Proposition 1. 
Refer to Appendix A. □
From Proposition 1, Corollary 1 can be obtained.
Corollary 1. 
For the GAMSNA-I with M sensors, the closed-form expression of uDOFs is u D O F s = M 2 2 + M + 5 in the case of M being even, and u D O F s = M 2 1 2 + M + 6 in the case of M being odd.
Proof of Corollary 1. 
According to (12) and Proposition 1, it follows that
u D O F s = 2 Q g + 1 = 2 ( M 2 + 1 ) M 2 + 5
In the case of M being even, M 2 = M 2 = M 2 , so u D O F s = M 2 2 + M + 5 . In the case of M being odd, M 2 = M 1 2 , M 2 = M + 1 2 , so u D O F s = M 2 1 2 + M + 6 . Therefore, the proof is completed. □
Figure 1a illustrates the structural composition of GAMSNA-I with M = 15 , specifically L 1 = { 0 , 1 , 7 } , L 2 = { 3 , 11 , 19 , 27 , 35 , 43 , 51 , 59 } , L 3 = { 60 } , L 4 = { 64 , 65 , 66 } . Figure 1b shows the positive result of DCA for GAMSNA-I, which is hole-free.
Corollary 1 indicates that the improved uDOFs are greater than those of NA. Moreover, according to [30], it is known that w ( 1 ) , w ( 2 ) and w ( 3 ) dominate the mutual coupling. Based on the analysis of (12), for M 12 , the first three weight functions of GAMSNA-I are
w ( 1 ) = M 1 3 w ( 2 ) = M 1 5 w ( 3 ) = 2 , M 1 = 6 M 1 6 , M 1 > 6
Compared to the first three weight functions of NA, which are w ( 1 ) = M 1 , w ( 2 ) = M 1 1 and w ( 3 ) = M 1 2 , the results of (14) demonstrate that the mutual coupling of GAMSNA-I is less. Nevertheless, there is a linear increase in the correlation between the weight functions of GAMSNA-I and the number of sensors. As the number of sensors increases, the impact of mutual coupling becomes more significant, imposing substantial limitations on the accuracy improvement of DOA estimation. Therefore, there is a need to improve GAMSNA-I to reduce its mutual coupling.

3.2. GAMSNA-II

According to (12), it can be noted that L 4 is a dense ULA with strong mutual coupling. Therefore, on the basis of GAMSNA-I, by adjusting some sensor positions of L 4 to reduce mutual coupling, GAMSNA-II is obtained.
The composition of GAMSNA-II is as follows:
L 1 = { 0 , 1 , M 1 } , M 1 4 L 2 = { 3 + ( M 1 + 1 ) | = 0 , 1 , , M 2 1 } L 3 = { M 2 ( M 1 + 1 ) M 1 + 3 } L 4 = { ( M 1 + 1 ) ( M 2 + 1 ) M 1 + 2 η | η [ 1 , M 1 4 ] , η is odd . } L 5 = { ( M 1 + 1 ) ( M 2 + 1 ) + 3 v 5 | v 5 [ 1 , M 1 4 ] , v 5 is even . }
Note that GAMSNA-II is only feasible with | L 4 | 3 , i.e., M 14 . Figure 2 shows the array structure of GAMSNA-II. It is worth noting that despite the presence of holes in the DCA of GAMSNA-II, the uDOFs corresponding to the central ULA portion of the DCA remain consistent with those of GAMSNA-I.
Proposition 2. 
The continuous part of the DCA for GAMSNA-II defined in (15) is equal to that of GAMSNA-I, i.e., U G A M S N A I I = U G A M S N A I = [ Q g , Q g ] with Q g = ( M 2 + 1 ) M 2 + 2 .
Proof of Proposition 2. 
Refer to Appendix B. □
It can be inferred from (15) that the first three weight functions of GAMSNA-II are
w ( 1 ) = 2 w ( 2 ) = M 1 5 w ( 3 ) = 2 , M 1 = 6 1 , M 1 > 6
Compared with GAMSNA-I, the mutual coupling of GAMSNA-II is less, and as the number of sensors increases, their differences become increasingly apparent. In addition, the expressions of uDOFs and sensor number for different array configurations are compared in Table 2, where M i represents the number of sensors in subarray L i of the corresponding sensor array.

4. Numerical Simulations

In this section, numerical simulations are performed to assess the performance of the proposed GAMSNA-I and GAMSNA-II from three key aspects: uDOFs, coupling leakage, and DOA estimation accuracy.

4.1. uDOFs

In the first numerical simulation, we compare the uDOFs of CPA, NA, SuperNA, ENA, NA-TS, GAMSNA-I and GAMSNA-II. The numerical relationship between the uDOFs of different arrays and the sensor quantity M is shown in Figure 3, where M = [14:1:30]. It can be observed that due to the presence of a large number of holes, the uDOFs curve for CPA is the lowest and considerably lower than the other curves. The uDOFs curves of SuperNA and NA overlap, and they are lower compared to the uDOFs curves of other NA-type arrays, indicating that SuperNA does not have an advantage over NA in terms of uDOFs and performs worse than other NA-type arrays. The structural uniqueness of ENA makes its uDOFs dependent on the parity of M. Specifically, when M is even, its uDOFs is greater than that of NA, while when M is odd, its uDOFs is equal to that of NA. The uDOFs curve for NA-TS is higher than the uDOF curves for the first three NA-type arrays, indicating its superior uDOFs performance. Finally, the uDOFs curves for the proposed GAMSNA-I and GAMSNA-II overlap with each other and are higher than the other curves, suggesting that these two array configurations exhibit the same excellent performance in terms of uDOFs. These results validate the effectiveness of the uDOFs property summarized for different array configurations in Table 2.

4.2. Coupling Leakage

In the second numerical simulation, we compare the weight functions and coupling leakage of CPA, NA, SuperNA, ENA, NA-TS, MISC, GENAMS, GAMSNA-I, and GAMSNA-II. In the interest of fairness, consider that these arrays uniformly consist of 15 sensors, with the specific sensor positions listed in Table 3.
According to Definition 3, the results of the virtual element weight functions for seven kinds of arrays are shown in Figure 4. To minimize the impact of mutual coupling effects as much as possible, weight functions with larger numerical values should be located far away from u = 0 . By observing Figure 4a, it can be noted that the weight function values near u = 0 are all 1, indicating that the mutual coupling effect of CPA is relatively small. In Figure 4b, a large number of weight functions with large numerical values are concentrated near u = 0 , reflecting the significant mutual coupling in NA. Subsequent Figure 4c–i disperse weight functions with large numerical values to varying degrees, away from u = 0 , thus confirming that their corresponding arrays exhibit reduced mutual coupling compared to NA. By comparison, it can be observed that MISC and the proposed GAMSNA-II have fewer weight functions with large numerical values near u = 0 , and therefore, it can be concluded that they have the least mutual coupling among all array configurations being compared. To provide a more intuitive representation, Table 4 presents a numerical comparison of the first three weight functions, and coupling leakage. Note that the values of the coupling coefficients are c 1 = 0.3 e ( j π / 3 ) , c i = c 1 e [ j ( i 1 ) π / 8 ] / i , 2 i B , where B = 100 [31].

4.3. DOA Estimation Performance

The DOA estimation performances of various array configurations are compared in this subsection using the MUSIC algorithm, where the peak search step size is set to 0.1°. The array configurations are the same as those given in Table 3.

4.3.1. MUSIC Spectrum

This experiment involves comparing the MUSIC spectrum of CPA, NA, SuperNA, ENA, NA-TS, MISC, GENAMS, GAMSNA-I, and GAMSNA-II under the scenario where there are 13 sources with DOAs of Θ = [ 30 : 5 : 30 ] , a signal-to-noise ratio (SNR) of 0 dB, and a snapshot number of 500.
The results are presented in Figure 5, where the red dashed lines and the black solid lines, respectively, indicate the true DOAs and the estimated spectra. In Figure 5a, the number of peaks in the spectral curve is inconsistent with the number of sources, and the actual values and estimated values hardly match, indicating poor performance of CPA in this experiment. Furthermore, when observing Figure 5b,d, it can be noticed that although the estimated DOA values are near the true values, there are perturbations in the peaks. This is due to significant mutual coupling of NA and ENA, which will lead to significant errors in the final DOA estimation results. In Figure 5e, although there are no peak disturbances, the peaks are not sharp enough. This is because the larger uDOFs of NA-TS balance the effects of mutual coupling, but the degree of balance is relatively low. Based on the observations from Figure 5c,f–i, it is evident that the DOA estimation results for SuperNA, MISC, GENAMS, GAMSNA-I, and GAMSNA-II are notably accurate. This is consistent with their own high uDOFs and low mutual coupling, as detailed in Table 4. Among these, due to the fact that the mutual coupling of GAMSNA-I is higher than SuperNA, although it has greater uDOFs, the sharpness of its spectral curve peaks is still inferior to SuperNA. The significantly high uDOFs of GENAMS balance the influence of mutual coupling well, making its MUSIC spectrum sharper than GAMSNA-I, but still inferior to SuperNA. This result suggests that, when the difference in uDOFs is within a certain range, lower mutual coupling can appropriately compensate for the disadvantages of uDOFs. The DOA estimation results of MISC are superior to that of GENAMS, which also confirms this conclusion. Additionally, when comparing Figure 5f,i, it can be observed that the peak stability of the spectral curves for MISC and GAMSNA-II is essentially consistent. This suggests that both exhibit similar low mutual coupling effects, ultimately resulting in similar DOA estimation outcomes.

4.3.2. Root Mean Square Error Performance

The next simulation compares the root mean square error (RMSE) performance of CPA, NA, SuperNA, ENA, NA-TS, GAMSNA-I, and GAMSNA-II versus the input SNR and the number of snapshots. The fixed parameter settings are N = 500 snapshots, S N R = 0 dB, and K = 13 sources with θ k U ( 90 , 90 ) , where U ( · ) means uniform distribution. RMSE is expressed as follows:
R M S E = 1 C K ( c = 1 C k = 1 K ( θ k θ ^ k , c ) 2 ,
where C represents the number of Monte Carlo trials, and it is set to 500 in this subsection.
RMSE performance versus the input SNR is shown in Figure 6, where S N R = [ 10 : 5 : 20 ] dB. We observe that throughout the entire SNR range, the proposed GAMSNA-I and GAMSNA-II exhibit better performance than other arrays. Among these, due to the smaller mutual coupling effects, the advantage of GAMSNA-II is more pronounced, as evidenced by its lower RMSE curve. The results indicate that the performance of these two arrays is superior to other arrays, with GAMSNA-II demonstrating the best performance.
Figure 7 illustrates RMSE performance versus the number of snapshots, with the set of snapshot numbers as [ 100 : 100 : 1000 ] . By observation, it can be noticed that when the number of snapshots is less than 300, the differences in RMSE performance among the arrays are relatively pronounced, and the RMSE curves of GAMSNA-I and GAMSNA-II are lower than the others. As the number of snapshots increases—except for GAMSNA-II, which exhibits a relatively significant improvement in RMSE performance—the RMSE performance of the other arrays remains relatively constant with minimal differences between them. Throughout the experiment, GAMSNA-II consistently outperforms GAMSNA-I in RMSE performance. The above-mentioned phenomenon indicates that the proposed GAMSNA-I and GAMSNA-II exhibit superior DOA estimation performance, with GAMSNA-II demonstrating lower mutual coupling, thus offering a significant advantage in accurately estimating DOA.

5. Conclusions

In this paper, a generalized augmented multi-subarray nested array (GAMSNA-I) with a large number of uDOFs and small mutual coupling is proposed. Additionally, its variant, GAMSNA-II, is developed to further reduce mutual coupling. The array configuration and uDOFs of GAMSNA-I are described by closed-form expressions constrained by the number of sensors. Considering the existence of a dense ULA subarray in GAMSNA-I, GAMSNA-II is developed with a sparser array configuration, while its uDOFs is proven to remain unaltered. Theoretical analysis and numerical simulation results demonstrate the effectiveness and superiority of the proposed GAMSNA-I and GAMSNA-II in terms of uDOFs, coupling leakage, spectral estimation, and DOA estimation accuracy.

Author Contributions

Conceptualization, S.W. and G.Z.; methodology, S.W.; software, G.Z.; validation, S.W., G.Z. and Y.S.; formal analysis, S.W.; investigation, S.W.; resources, S.W. and Y.S.; data curation, S.W.; writing—original draft preparation, G.Z.; writing—review and editing, S.W. and Y.S.; visualization, S.W., G.Z. and Y.S.; supervision, S.W., G.Z. and Y.S.; project administration, Y.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of Shanghai under grant No. 19ZR1437600, and the Foundation of Shanghai Key Laboratory of Navigation and Location Based Services under grant No. C-6105-23-021.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Proof of Proposition 1

Since D is symmetric about position zero, for simplifying the proof steps, only the positive part of D is proven to be continuous. Defining symbol D i , j = L i L j , i j , according to (12), the DCAs between the four subarrays are
D 11 = { 0 , 1 , M 1 1 , M 1 } D 21 = { 3 + ( M 1 + 1 ) , 2 + ( M 1 + 1 ) , ( 3 M 1 ) + ( M 1 + 1 ) | = [ 0 , M 2 1 ] } D 31 = { ( M 2 1 ) ( M 1 + 1 ) + 4 , ( M 2 1 ) ( M 1 + 1 ) + 3 , ( M 2 2 ) ( M 1 + 1 ) + 5 } D 41 = { M 2 ( M 1 + 1 ) + 3 η , M 2 ( M 1 + 1 ) + 2 η , ( M 2 1 ) ( M 1 + 1 ) + 4 η | η = [ 1 , M 1 4 ] } D 22 = { ( M 1 + 1 ) | = [ 0 , M 2 1 ] } D 32 = { ( M 2 ) ( M 1 + 1 ) M 1 | = [ 0 , M 2 1 ] } D 42 = { ( M 2 ) ( M 1 + 1 ) η | = [ 0 , M 2 1 ] , η = [ 1 , M 1 4 ] } D 33 = { 0 } D 43 = { M 1 η | η = [ 1 , M 1 4 ] } D 44 = { [ 0 , M 1 5 ] }
Since L 2 is a set of arithmetic sequences with a common difference of M 1 + 1 , and the last element of L 1 , which is M 1 , falls between the first and second elements of L 2 , it is only necessary to ensure continuity from 0 to M 1 + 4 to guarantee continuity from 0 to 3 + ( M 2 1 ) ( M 1 + 1 ) .
When = 0 , the combination of D 11 and D 21 contains a continuous range [ 0 , 3 ] . D 21 has ‘4’ (when = 1 ). Considering = M 2 1 , D 42 has a continuous subset [ 5 , M 1 ] . D 22 contains ‘ M 1 + 1 ’. Let = M 2 2 ; the element of D 32 is ‘ M 1 + 2 ’. In addition, D 21 includes ‘ M 1 + 3 ’ and ‘ M 1 + 4 ’ (when = 1 ). In summary, the continuous range [ 0 , M 1 + 4 ] can be obtained through D 11 , D 21 , D 22 , D 32 and D 42 , that is, the continuous set [ 0 , 3 + ( M 2 1 ) ( M 1 + 1 ) ] can be obtained. Furthermore, D 31 contains ‘ 4 + ( M 2 1 ) ( M 1 + 1 ) ’. In D 42 , let = 0 obtain the continuous part [ 5 + ( M 2 1 ) ( M 1 + 1 ) , M 2 ( M 1 + 1 ) 1 ] . The continuous range of [ M 2 ( M 1 + 1 ) , M 2 ( M 1 + 1 ) + 2 ] can be provided by D 41 . Thus, the continuous range [ 0 , M 2 ( M 1 + 1 ) + 2 ] exists, which completes the proof.

Appendix B. Proof of Proposition 2

The L 4 and L 5 of GAMSNA-II originate from the reallocation of the L 4 of GAMSNA-I. Specifically, L 4 is the odd-numbered part of η in L 4 , while L 5 is the even-numbered part of η in L 4 plus M 1 + 1 . Therefore, it is only necessary to prove that L 5 has the same effect as the even-numbered part of η in L 4 . For simplicity, only the expressions for D 42 , D 43 and D 52 are given.
D 42 = { ( M 2 ) ( M 1 + 1 ) η | = [ 0 , M 2 1 ] , η = [ 1 , M 1 4 ] , η is odd } D 43 = { M 1 η | η = [ 1 , M 1 4 ] , η is odd } D 52 = { [ M 2 ( 1 ) ] ( M 1 + 1 ) η | = [ 0 , M 2 1 ] , η = [ 1 , M 1 4 ] , η is even }
It is easy to know that 1 = [ 1 , M 2 2 ] . So, when η is even, the elements that D 42 has but D 52 does not have are ( M 1 + 1 ) η . It is worth noting that D 43 contains all the elements of { ( M 1 + 1 ) η | η i s e v e n } . Therefore, when considering the contribution to ensure that continuities of DCA, D 42 , D 52 , and D 43 are equivalent to D 42 . It is worth noting that D 43 is a subset of D 43 , which was not involved in Proof of Proposition 1. Therefore, the use of elements of D 43 in this proof does not have an impact on the continuity of DCA. In conclusion, the continuous part of the DCA of GAMSNA-II is completely equal to that of GAMSNA-I.

References

  1. Chen, Y.; Yang, M.; Li, J.; Zhang, X. A Nested–Nested Sparse Array Specially for Monostatic Colocated MIMO Radar with Increased Degree of Freedom. Sensors 2023, 23, 9230. [Google Scholar] [CrossRef]
  2. Liu, W.; Haardt, M.; Greco, M.S.; Mecklenbräuker, C.F.; Willett, P. Twenty-Five Years of Sensor Array and Multichannel Signal Processing: A review of progress to date and potential research directions. IEEE Signal Process. Mag. 2023, 40, 80–91. [Google Scholar] [CrossRef]
  3. Li, J.; Li, P.; Li, P.; Tang, L.; Zhang, X.l.; Wu, Q. Self-position awareness based on cascade direct localization over multiple source data. IEEE Trans. Intell. Transp. Syst. 2022, 25, 796–804. [Google Scholar] [CrossRef]
  4. Zhang, Y.; Zhou, H.; Hu, G. Further Results on the Degrees of Freedom for Coprime MIMO Systems. IEEE Commun. Lett. 2023, 27, 3335–3339. [Google Scholar] [CrossRef]
  5. Shi, J.; Hu, G.; Zong, B.; Chen, M. DOA estimation using multipath echo power for MIMO radar in low-grazing angle. IEEE Sensors J. 2016, 16, 6087–6094. [Google Scholar] [CrossRef]
  6. Fang, Y.; Zhu, S.; Gao, Y.; Lan, L.; Zeng, C.; Liu, Z. Direction-of-arrival estimation of coherent signals for uniform linear antenna arrays with mutual coupling in unknown nonuniform noise. IEEE Trans. Veh. Technol. 2021, 71, 1656–1668. [Google Scholar] [CrossRef]
  7. 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 Wirel. Commun. Lett. 2019, 8, 1004–1007. [Google Scholar] [CrossRef]
  8. Li, B.; Zhang, X.; Li, J.; Ma, P. DOA estimation of non-circular source for large uniform linear array with a single snapshot: Extended DFT method. IEEE Commun. Lett. 2021, 25, 3843–3847. [Google Scholar] [CrossRef]
  9. Wu, R.; Wang, M.; Zhang, Z. Computationally efficient DOA and carrier estimation for coherent signal using single snapshot and its time-delay replications. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 2469–2480. [Google Scholar] [CrossRef]
  10. Patwari, A. Sparse Linear Antenna Arrays: A Review. In Antenna Systems; IntechOpen: London, UK, 2022. [Google Scholar]
  11. Moffet, A. Minimum-redundancy linear arrays. IEEE Trans. Antennas Propag. 1968, 16, 142–149. [Google Scholar] [CrossRef]
  12. 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]
  13. Vaidyanathan, P.P.; Pal, P. Sparse sensing with co-prime samplers and arrays. IEEE Trans. Signal Process. 2010, 59, 573–586. [Google Scholar] [CrossRef]
  14. Qin, S.; Zhang, Y.D.; Amin, M.G. Generalized coprime array configurations for direction-of-arrival estimation. IEEE Trans. Signal Process. 2015, 63, 1377–1390. [Google Scholar] [CrossRef]
  15. Zhang, Y.; Hu, G.; Zhang, F.; Zhou, H. Enhanced CACIS configuration for direction of arrival estimation. Electron. Lett. 2022, 58, 737–739. [Google Scholar] [CrossRef]
  16. Raza, A.; Liu, W.; Shen, Q. Thinned coprime array for second-order difference co-array generation with reduced mutual coupling. IEEE Trans. Signal Process. 2019, 67, 2052–2065. [Google Scholar] [CrossRef]
  17. Zheng, W.; Zhang, X.; Wang, Y.; Shen, J. Padded coprime arrays for improved DOA estimation: Exploiting hole representation and filling strategies. IEEE Trans. Signal Process. 2020, 68, 4597–4611. [Google Scholar] [CrossRef]
  18. Patra, R.K.; Dhar, A.S. A Novel k-times Extended Coprime Array for DOA Estimation With Increased Degrees of Freedom. IEEE Signal Process. Lett. 2022, 29, 1402–1406. [Google Scholar] [CrossRef]
  19. Shi, J.; Hu, G.; Zhang, X.; Zhou, H. Generalized nested array: Optimization for degrees of freedom and mutual coupling. IEEE Commun. Lett. 2018, 22, 1208–1211. [Google Scholar] [CrossRef]
  20. Zhao, P.; Hu, G.; Qu, Z.; Wang, L. Enhanced nested array configuration with hole-free co-array and increasing degrees of freedom for DOA estimation. IEEE Commun. Lett. 2019, 23, 2224–2228. [Google Scholar] [CrossRef]
  21. Yang, M.; Sun, L.; Yuan, X.; Chen, B. Improved nested array with hole-free DCA and more degrees of freedom. Electron. Lett. 2016, 52, 2068–2070. [Google Scholar] [CrossRef]
  22. Iizuka, Y.; Ichige, K. Extension of nested array for large aperture and high degree of freedom. IEICE Commun. Express 2017, 652, 381–386. [Google Scholar] [CrossRef]
  23. 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]
  24. 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]
  25. Liu, J.; Zhang, Y.D.; 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]
  26. Zhang, Y.; Hu, G.; Zhou, H.; Bai, J.; Zhan, C.; Guo, S. Hole-Free Nested Array with Three Sub-ULAs for Direction of Arrival Estimation. Sensors 2023, 23, 5214. [Google Scholar] [CrossRef]
  27. Shi, W.; Vorobyov, S.A.; Li, Y. ULA fitting for sparse array design. IEEE Trans. Signal Process. 2021, 69, 6431–6447. [Google Scholar] [CrossRef]
  28. 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]
  29. 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]
  30. Zheng, Z.; Wang, W.Q.; Kong, Y.; Zhang, Y.D. MISC array: A new sparse array design achieving increased degrees of freedom and reduced mutual coupling effect. IEEE Trans. Signal Process. 2019, 67, 1728–1741. [Google Scholar] [CrossRef]
  31. Zhou, C.; Gu, Y.; He, S.; Shi, Z. A robust and efficient algorithm for coprime array adaptive beamforming. IEEE Trans. Veh. Technol. 2017, 67, 1099–1112. [Google Scholar] [CrossRef]
Figure 1. An example of GAMSNA-I with M = 15 , and the positive result of its DCA. (a) Sensor positions of GAMSNA-I. (b) Positive DCA of GAMSNA-I.
Figure 1. An example of GAMSNA-I with M = 15 , and the positive result of its DCA. (a) Sensor positions of GAMSNA-I. (b) Positive DCA of GAMSNA-I.
Sensors 24 00808 g001
Figure 2. An example of GAMSNA-II with M = 15 , and the positive result of its DCA. (a) Sensor positions of GAMSNA-II. (b) Positive DCA of GAMSNA-II.
Figure 2. An example of GAMSNA-II with M = 15 , and the positive result of its DCA. (a) Sensor positions of GAMSNA-II. (b) Positive DCA of GAMSNA-II.
Sensors 24 00808 g002
Figure 3. uDOFs versus M for different arrays.
Figure 3. uDOFs versus M for different arrays.
Sensors 24 00808 g003
Figure 4. The virtual element weight functions of nine kinds of 15-element sparse arrays. (a) CPA. (b) NA. (c) Super NA. (d) ENA. (e) NA-TS. (f) MISC. (g) GENAMS. (h) GAMSNA-I. (i) GAMSNA-II.
Figure 4. The virtual element weight functions of nine kinds of 15-element sparse arrays. (a) CPA. (b) NA. (c) Super NA. (d) ENA. (e) NA-TS. (f) MISC. (g) GENAMS. (h) GAMSNA-I. (i) GAMSNA-II.
Sensors 24 00808 g004
Figure 5. MUSIC spectrum of nine kinds of 15-element sparse arrays. (a) CPA. (b) NA. (c) Super NA. (d) ENA. (e) NA-TS. (f) MISC. (g) GENAMS. (h) GAMSNA-I. (i) GAMSNA-II.
Figure 5. MUSIC spectrum of nine kinds of 15-element sparse arrays. (a) CPA. (b) NA. (c) Super NA. (d) ENA. (e) NA-TS. (f) MISC. (g) GENAMS. (h) GAMSNA-I. (i) GAMSNA-II.
Sensors 24 00808 g005
Figure 6. RMSE performance versus input SNR.
Figure 6. RMSE performance versus input SNR.
Sensors 24 00808 g006
Figure 7. RMSE performance versus number of snapshots.
Figure 7. RMSE performance versus number of snapshots.
Sensors 24 00808 g007
Table 1. Optimal sensor quantity configuration for two subarrays of NA.
Table 1. Optimal sensor quantity configuration for two subarrays of NA.
The Total Number of SensorsSubarray 1Subarray 2
M M 1 M 2
Odd ( M 1 ) / 2 ( M + 1 ) / 2
Even M / 2 M / 2
Table 2. The relationship between the number of sensors and the uDOFs for seven types of sparse arrays.
Table 2. The relationship between the number of sensors and the uDOFs for seven types of sparse arrays.
Array ConfigurationSensor NumberuDOFs
CPA M + N 2 ( M + N ) 1
NA M 1 + M 2 2 M 2 ( M 1 + 1 ) 1
SuperNA M 1 + M 2 2 M 2 ( M 1 + 1 ) 1
ENA M 1 + M 2 2 M 1 ( M 2 + 1 ) + 1
NA-TSeven: M 1 = ( M 2 M 3 ) / 2 ,   M 2 = M / 2 , 1 M 3 < M / 2 M 2 / 2 + 2 M 2 M 3 1
odd: M 1 = ( M 2 M 3 1 ) / 2 ,   M 2 = ( M + 1 ) / 2 ,   1 M 3 < ( M 1 ) / 2 M 2 / 2 + 2 M 2 M 3 3 / 2
GAMSNA-Ieven: M 1 = 3 ,   M 2 = M / 2 ,   M 3 = 1 , M 4 = M / 2 4 ,   M 12 M 2 / 2 + M + 5
odd: M 1 = 3 ,   M 2 = ( M + 1 ) / 2 ,   M 3 = 1 , M 4 = ( M + 1 ) / 2 5 ,   M 12 ( M 2 1 ) / 2 + M + 6
GAMSNA-IIeven: M 1 = 3 ,   M 2 = M / 2 ,   M 3 = 1 , M 4 = M / 4 2 ,   M 5 = M / 2 4 M 4 ,   M 14 M 2 / 2 + M + 5
odd: M 1 = 3 ,   M 2 = ( M + 1 ) / 2 ,   M 3 = 1 , M 4 = ( M + 1 ) / 4 5 / 2 ,   M 5 = ( M + 1 ) / 2 5 M 4 ,   M 14 ( M 2 1 ) / 2 + M + 6
Table 3. Sensors locations for nine sparse arrays.
Table 3. Sensors locations for nine sparse arrays.
Array ConfigurationSensor Locations
CPA{0, 7, 8, 14, 16, 21, 24, 28, 32, 35, 40, 42, 48, 49}
NA{0, 1, 2, 3, 4, 5, 6, 7, 15, 23, 31, 39, 47, 55, 63}
SuperNA{0, 2, 4, 6, 9, 11, 13, 15, 23, 31, 39, 47, 55, 62, 63}
ENA{0, 1, 2, 3, 4, 5, 6, 14, 21, 28, 35, 42, 49, 56, 63}
NA-TS{0, 1, 2, 3, 11, 19, 27, 35, 43, 51, 59, 63, 64, 65, 66}
MISC{0, 1, 6, 14, 22, 30, 38, 46, 54, 62, 64, 66, 69, 71, 73}
GENAMS{0, 1, 2, 5, 10, 15, 26, 37, 48, 59, 65, 71, 77, 78, 79}
GAMSNA-I{0, 1, 3, 7, 11, 19, 27, 35, 43, 51, 59, 60, 64, 65, 66}
GAMSNA-II{0, 1, 3, 7, 11, 19, 27, 35, 43, 51, 59, 60, 64, 66, 73}
Table 4. A summary of the uDOFs, virtual element weight functions, and coupling leakage for nine types of sparse arrays with 15 sensors.
Table 4. A summary of the uDOFs, virtual element weight functions, and coupling leakage for nine types of sparse arrays with 15 sensors.
Array ConfigurationCPANASuperNAENANA-TS
uDOFs29127127127133
w(1)27166
w(2)26654
w(3)25142
L0.20540.32360.19700.30100.2894
Array ConfigurationMISCGENAMSGAMSNA-IGAMSNA-II
uDOFs147159133133
w(1)1442
w(2)4222
w(3)1111
L0.17720.23680.24340.1928
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wei, S.; Zhu, G.; Su, Y. A Novel Sparse Array Configuration for Direction of Arrival Estimation with Increased Uniform Degrees of Freedom and Reduced Mutual Coupling. Sensors 2024, 24, 808. https://doi.org/10.3390/s24030808

AMA Style

Wei S, Zhu G, Su Y. A Novel Sparse Array Configuration for Direction of Arrival Estimation with Increased Uniform Degrees of Freedom and Reduced Mutual Coupling. Sensors. 2024; 24(3):808. https://doi.org/10.3390/s24030808

Chicago/Turabian Style

Wei, Shuang, Gencun Zhu, and Ying Su. 2024. "A Novel Sparse Array Configuration for Direction of Arrival Estimation with Increased Uniform Degrees of Freedom and Reduced Mutual Coupling" Sensors 24, no. 3: 808. https://doi.org/10.3390/s24030808

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