# A Novel Sparse Array for Localization of Mixed Near-Field and Far-Field Sources

Authors: Yinsheng Wang; WeiJia Cui; Yuxi Du; Bin Ba; Fengtong Mei
Journal: International Journal of Antennas and Propagation
 (2021)

---

## Abstract

As we all know, nested array can obtain a larger array aperture and more degrees of freedom using fewer sensors. In this study, we not only designed an enhanced symmetric nested array (ESNA), which achieved more consecutive lags and more unique lags compared with a generalized nested array but also developed a special cumulant matrix, in the case of a given number of sensors, which can automatically generate the largest consecutive lags of the array. First, the direction-of-arrivals (DOAs) of mixed sources are estimated using the special cumulant matrix. Then, we can estimate the range of the near-field source in the mixed source using a one-dimensional spectral search through estimated DOAs, and in the mixed sources, the near-field and far-field sources are classified by bringing in the range parameter. The largest consecutive lags and composition method of ESNA are also given, under a given number of sensors.Our algorithm has moderate computation complexity, which provides a higher resolution and improves the parameters’ estimation accuracy. Numerical simulation results demonstrate that the proposed array showed an outstanding performance under estimation accuracy and resolution ability for both DOA and range estimation compared with existing arrays of the same physical array sensors.


---

## Body

## 1. Introduction

Mixed source localization is an important problem in the field of array signal processing such as radar, sonar, and communications [1–5]. Therefore, to solve this problem, many algorithms have been proposed, such as the ESPRIT algorithm [6], the MUSIC algorithm [7], and so on. However, most of them focused on far-field (FF) assumption to solve the far-field problem [8]. When a source is located in the Fresnel region, which belongs to the near-field (NF) source, and the wavefront is no longer a plane wave but a spherical wave [9], the wavefront is composed of DOA and range [10].Recently, many algorithms based on NF source localization have been developed to estimate DOA and the range [11–15]; however, the algorithms proposed above mainly focuses on NF or FF sources. In some other situations, for example, using microphone arrays to locate speakers and communications, mixed sources are received by an array [16]. Fortunately, many algorithms have been proposed to deal with this situation. Based on the classical second-order statistics algorithms, an oblique projection algorithm was proposed by He et al. [17] and Zuo et al. [18] to eliminate the subspace of FF sources from the covariance matrix and obtain the corresponding subspace of NF sources. Corresponding NF and FF sources could be well separated; then with the NF DOA estimates, corresponding ranges parameters could be obtained via 1-D peak searching. However, when there are only FF sources, this method fails. To solve this problem, a RARE-based localization algorithm was proposed by Hua et al. [19]. This method eliminated the range parameters in the signal subspace, then obtained the NF DOA, and finally substituting NF DOA into a 1-D peak searching to obtain the range parameters.It is generally known that high-order statistics (HOS) are getting more and more attention. HOS can not only increase the estimation accuracy but also extend an aperture of the given array. In addition, HOS is insensitive to Gaussian noise. A high-order MUSIC algorithm, proposed by Liang and Ding [20], DOA, and range are solved by constructing cumulant matrix, and it needed to construct two fourth-order cumulant matrices, which leads to high computation complexity. However, the methods mentioned above are all based on a uniform symmetric linear array to estimate the DOA and ranges of the mixed sources.Noteworthy, in the case of coexistence of the near- and far-field sources, existing works mentioned above can be used. However, due to the use of uniform symmetric linear array, when physical array sensors are certain, consecutive lag numbers will be limited [21]. The problem of using fewer physical sensors to estimate more sources is becoming more and more important [22]. To solve this problem, many sparse linear arrays, such as nested arrays [23] and coprime arrays [24], have been proposed. The sparse array can produce a larger array aperture and more consecutive lags compared with uniform linear array with the same array sensors.Some algorithms based on nested array for mixed sources localization have been developed. For example, using convolution neural networks to localization mixed sources [25], using exact spatial propagation geometry localization mixed sources [26], and components separation for localization and classification of mixed sources [27]. The methods mentioned above all used general nested arrays. Although the number of consecutive lags can be increased, they were not an optimal array.According to the characteristics of a sparse array, an enhanced symmetric nested array (ESNA) and a novel algorithm for localization of mixed sources are proposed in this paper. By exploiting subarray partition method of the ESNA, one cumulant matrix that contains only DOA information to estimate the DOA of mixed sources was constructed using the conventional MUSIC method. Then, with the estimated DOA, the range parameters can be obtained by a 1-D peak searching. Compared with the existing arrays, ESNA not only exhibits a larger number of consecutive lags but also provides a larger array aperture under the given same sensor number. Compared with methods of Wang et al. in [28] and Akbar et al. in [10], due to use the ESNA and subarray partition method, we can obtain more continuous lags and better performance under the condition of the same sensor number.In this paper, we designed the ESNA and a new algorithm to solve the mixed-field sources localization. Our main contributions are as follows:(1)
We proposed an enhanced symmetric nested array (ESNA). We also give the ESNA optimal array configuration parameters. Compared with uniform symmetric array, ESNA can provide higher spatial resolution and enhanced degrees of freedom with same sensors.(2)
Based on ESNA, we proposed a new algorithm, subarray partition, and compared with the existing similar algorithms, the complexity of this algorithm is greatly reduced.(3)
We verified the superiority of the proposed algorithm based on the ESNA in terms of mixed-field accuracy, resolution capacity, and many more DOFs and analyzed the range of ambiguity problem.The remaining part of this paper is organized as follows: Section2 introduces the signal model. Section 3 introduces the ESNA configuration and proposed algorithm. Section 4 provides the simulation results; the conclusions are given in Section 5.


## 2. Signal Model

ConsiderK narrow band and independent sources parameterized by (θk,rk), k=1,2,…,K impinge on the ESNA with 2M+2N+1 sensors, from Figure 1, which consists of two symmetric sparse arrays with intersensor spacing 2Nd contains 2M sensors; and a symmetric uniform linear array with intersensor spacing d contains 2N−1 sensors; and a left sensor, a right sensor with intersensor spacing 2N−1d on both sides of the array. The distance between two sparse symmetric linear subarrays and a uniform linear array is 2Nd. The sensors index can be expressed as −M−N,−M−N+1,…,M+N−1,M+N. We use set W to represent the position of the sensors:(1)W=p−L⋅d,…,p0⋅d,…,pL⋅d,where L=M+N. The position of the ith sensor is pi⋅d, i ∈−L,L. The value of pi is defined as(2)pi=i,i≤N−1,N−1+2Ni−N+1,N≤i≤M+N−1,1−N+2Ni+N−1,1−M−N≤i≤−N,N−2+2NM+1,i=M+N,2−N−2NM+1,i=−M−N.Figure 1
The ESNA geometry.The reference point of the phase is the center of the array, that is, 0 array element is the reference point. When theith sensor receives a signal, it can be expressed as [20](3)xit=∑k=1Ksktejτik+nit,where skt represents the kth narrowband source and nik represents the additive Gaussian noise. In addition, τik represents the propagation time delay of the kth source arriving the ith sensor relative to the 0th sensor. If kth source represents near-field one, τik can be written in the following form:(4)τik=iγk+i2ϕk,where γk and ϕk are given by(5)γk=−2πdλsinθk,ϕk=πd2λrkcos2θk,where θk and rk represent the azimuth DOA of kth and range of the kth near-field source. According to the definition, near-field source lies in the Fresnel zone 0.62D3/λ1/2,2D2/λ [14], where D=2pLd is the array aperture. If the kth source is a far-field one, then rk⟶∞, and τik can be written the following form:(6)τik=iγk,where ϕk≈0 due to beyond range of far-field source.Assume that theK1 sources are near-field sources and the rest of K sources are far-field sources, then the vector form of equation (3) can be expressed as(7)xt=ANsNt+AFsFt+nt,where(8)xt=x−Lt,…,x0t,…,xLtT,AN=aθ1,r1,…,aθK1,rK1,AF=aθK1+1,…,aθK,sNt=s1t,…,sK1tT,sFt=sK1+1t,…,sKtT,and(9)nt=n−Lt,…,n0t,…,nLtT.Note that the forms of steering vectoraθk,rk and aθk have the following form:(10)aθk,rk=ejp−Lγk+p−L2ϕk,…,1,…,ejpLγk+pL2ϕkT,aθk=ejp−Lγk,…,1,…,ejpLγkT.Some assumptions are required to hold in this paper:(1)
The source signals are narrowband, statistically mutually independent, and non-Gaussian processes with nonzero kurtosis(2)
The sensor noise is the additive white Gaussian process and independent of the source signals(3)
Using the ESNA array, in which the interelement spacingd of the uniform linear array is set to be λ/4


## 3. Proposed Algorithm

### 3.1. Analysis of ESNA in the Cumulant Domain

Cumulant has many advantages, it not only improves the estimation performance but also expands the array aperture [29–31]. In addition, cumulant can also provide an increased number of consecutive lags, which is insensitive to Gaussian noise.Based on the assumptions, we define the fourth-order cumulantcumxmt,xn∗t,xρ∗t,xqt as(11)cumxmt,xn∗t,xρ∗t,xqt=cum∑k=1Ksktejpmγk+pm2ϕk,∑k=1Ksktejpnγk+pn2ϕk∗,∑k=1Ksktejpργk+pρ2ϕk∗,∑k=1Ksktejpqγk+pq2ϕk=∑k=1Kc4,skejpm−pn−pρ−pqγk+pm2−pn2−pρ2−pq2ϕk×cumskt,sk∗t,sk∗t,skt=∑k=1Kc4,skejpm−pn−pρ−pqγk+pm2−pn2−pρ2−pq2ϕk,where c4,sk =  cumskt,sk∗t,sk∗t,skt is the kurtosis of the kth signal and m,n,ρ,q∈−L,L.In equation (11), we remove the ϕk term, while keeping the γk term, it is assumed that n=−m,q=−ρ. Thus, equation (11) becomes(12)cumxmt,x−m∗t,xρ∗t,x−ρt=∑k=1Kc4,skej2pm−pργk.Obviously, equation (12) only contains DOA parameters, which we can estimate DOA of all sources first by constructing a special cumulant matrix.Proposition 1.
LetL=lρ|−M−N≤ρ≤M+N represents the positions of the sensor in the ESNA array, and Lc=lc|lc=pm−pρ represents the lags in the ESNA. Based on the ESNA, we define the following facts to hold:(a)
There are4MN+8N−3 consecutive lags in the range −4N−2MN+2,4N+2MN−2Proof (a): see Appendix.(b)
WhenN=⌈L+2/2⌉ and M=⌊L+2/2⌋, the largest consecutive lags can be obtained, where ⌈⋅⌉ and ⌊⋅⌋ are the round-towards infinite and round-towards zeroProof: we definefN=4NM+8N−3, where N=L−M. f′N=4L−8N+8 can be obtained by deriving N. Notably, N is an integer, when N=⌈L+2/2⌉ or M=⌊L−2/2⌋, f can get the maximum value.As a comparison, four array configurations with2L+1 sensors were considered: (1) symmetric nested array (SNA) [28], (2) symmetric double nested array (SDNA) [23], (3) improved symmetric nested array (ISNA) [10], and (4) the proposed enhanced symmetric nested array (ESNA). In Table 1, in the case of the same number of sensors, we compare the number of consecutive lags and the values of M and N. Among all the nested arrays, the total sensors of the uniform linear array in the middle of a nested array are 2N−1, and the sparse array on both side of a uniform linear array is 2M+2. From Table 1, it is obvious that the proposed ESNA has the maximum number of consecutive lags when the number of array sensors is the same. As the results show, the ESNA can offer more consecutive lags for DOA estimation in mixed source scenario.In order to provide an intuitive impression of above-mentioned four array configurations, herein, we give an example of 13 sensors and mark their positions in Figure2. The red circles represent the lags of four virtual coarrays, and the blue crosses represent the holes for four virtual coarrays. The interval between positions is d.Table 1
Consecutive lags numbers for different arrays.
Physical numberSNASDNAISNAESNAMaximum number of consecutive lagsNMNMNMNMSNASDNAISNAESNA711223130715192192223323111212933112233424117294145133334434232375561153344535231477177174445545339578997Figure 2
Sensor locations.


### 3.2. DOA Estimation of All Sources

In this section, we construct a special cumulant, which can be used to eliminate the range parameter in the near-field, and then use the cumulant to estimate the DOA of the mixed field.In the proposed algorithm, we use the subarray partition method to combine the sparse array and the uniform linear array in the ESNA form, and the corresponding cumulant output and the output cumulant are formed by the largest virtual consecutive lags.Based on equation (11), first, let m=M+N, p∈1−N,…,N−1, and the 2N−1×1 cumulant vector c1 can be obtained, and its N−ρ th element is given by(13)c1N−ρ=cumxM+Nt,x−M−N∗t,xρ∗t,x−ρ=∑k=1Kc4,skej2N−2+2NM+1−ργk.Similarly, since we are using a symmetric enhanced nested array, we can easily obtain a symmetrical part. Letm=−M−N, ρ∈1−N,…,N−1, and we get another 2N−1×1 cumulant vector c¯1, and its N−ρth element is given by(14)c¯1N−ρ=cumx−M−Nt,xM+N∗t,x−ρ∗t,xρt=∑k=1Kc4,skej2−2N−2NM+1−ργk.Second, letm∈N,N+M−1, p∈1−N,…,N−1, we get another M2N−1×1 cumulant vector c2_m, and its 2Nm−N+N−ρ th element is given by(15)c2_m2Nm−N+N−ρ=cumxmt,x−m∗t,x−ρ∗t,xρt=∑k=1Kc4,skej2N−1+2Nm−N+1−ργk.Similarly, letm∈N,N+M−2, ρ=M+N and construct M−1×1 vector c¯2_m, and its m−N+1th element is given by(16)c¯2_mm−N+1=cumxM+Nt,x−M−N∗t,xm∗t,x−mt=∑k=1Kc4,skej22NM+N−m−1γk.Now, we construct the2MN−1×1 vector c2,(17)c2=c2_1,c¯2_M+N−2,…,c2_M+N−2,c¯2_1,c2_M+N−1.Similarly, since we are using a symmetric enhanced nested array, we can easily obtain a symmetrical part. Letm∈1−M−N,−N, ρ∈1−N,…,N−1, we get another M2N−1×1 cumulant vector c2_m∗, and its −2Nm+N+N−ρth element is given by(18)c2_m∗−2Nm+N+N−ρ=cumxmt,x−m∗t,xρ∗t,xρt=∑k=1Kc4,skej21−N+2Nm+N−1−ργk.Similarly, letm∈2−M−N,−N, ρ=−M−N and construct M−1×1 vector c¯2_m∗, and its m+N+2th element is given by(19)c¯2_m∗m+N+2=cumx−M−Nt,xM+N∗t,xm∗t,x−mt=∑k=1Kc4,skej21−2NM+N+mγk.Now, we construct the2MN−1×1 vector c¯2,(20)c¯2=c2_1∗,c¯2_M+N−2∗,…,c2_M+N−2∗,c¯2_1∗,c2_M+N−1∗.Third, letm∈1−N,N−1, ρ=0 and construct the 2N−1×1 vector c3, and its mth element is given by(21)c3m+N=cumxmt,x−m∗t,x0∗t,x0t=∑k=1Kc4,skej2mγk.Fourth, letm=N−1, ρ∈1−N,−1, we get another N−1×1 cumulant vector c4, and its N+ρth element is given by(22)c4N+ρ=cumxmt,x−m∗t,xρ∗t,xρt=∑k=1Kc4,skej2m−ργk.Similarly, since we are using a symmetric enhanced nested array, we can easily obtain a symmetrical part, letm=1−N, ρ∈1,N−1, we construct the N−1×1 vector c¯4, and its ρth element is given by(23)c¯4ρ=cumxmt,x−m∗t,xρ∗t,xρt=∑k=1Kc4,skej2m−ργk.Fifth, letm=M+N, ρ=M+N−1 and m=−M−N, ρ=1−M−N, and we construct the last two cumulants c5 and c¯5,(24)c5=cumxM+Nt,x−M−N∗t,xM+N−1∗t,x1−M−Nt=∑k=1Kc4,skej22N−1γk,and(25)c¯5=cumx−M−Nt,xM+N∗t,x1−M−N∗t,xM+N−1t=∑k=1Kc4,skej21−2Nγk.Finally, letm=M+N−1, ρ=−N, we construct cend, and it is given by(26)cend=cumxM+N−1t,x1−M−N∗t,x−N∗t,xNt=∑k=1Kc4,skej24N+2NM−2γk.Similarly, since we are using a symmetric enhanced nested array, we can easily obtain a symmetrical part. Letm=1−M−N, ρ=N, we construct c¯end, and it is given by(27)c¯end=cumx1−M−Nt,xM+N−1∗t,xN∗t,x−Nt=∑k=1Kc4,skej24N+2NM−2γk.Now, we combine the above vectors to form a24N+2MN−1−1×1 vector c, which is expressed as(28)c=c¯end,c¯1T,c¯2T,c¯5,c¯4T,c3T,c4T,c5,c2T,c1T,cendT.A Toeplitz matrixC is constructed through the vector c, and matrix C is 4N+2MN−1×4N+2MN−1, whose mth column of C consists of the 4N+2MN−1+1−mth to 24N+2MN−1+1−mth element of c. The Toeplitz matrix C can be expressed in a compact form as(29)C=AθCsAHθ,Cs=diagc4s1,…,c4sK,Aθ=aθ1,…,aθK,and(30)aθk=1,ej2γ1,…,ej24N+2MN−2γkT.Since the sources signals have nonzero fourth-order cumulant, the matrixC is full rank. Therefore, we estimate the common θ of mixed sources by performing the eigenvalue decomposition of C.Based on the EVD ofC, using the MUSIC algorithm [6] to estimate DOA and we can get the following formula:(31)θk=minaθHUnUnHaθ,k=1,2,…,K,where Un represents the noise subspace of C.


### 3.3. Range Estimation and Source Type Classification

Onceθ¯k,k=1,…,K is estimated, then we can estimate rk. When we get the convariance matrix R=ExtxHt, perform eigen-decomposition on it, and we can get(32)R=U¯sV¯sU¯sH+U¯nV¯nU¯nH,where U¯n spans the noise subspace, which contains 4N+2MN−1−K eigenvectors. By taking the estimated θ¯k into aθk,r, the rk can be obtained by the following:(33)rk=minaθk,rU¯nU¯nHaHθk,r,k=1,…,K.Obviously,K peaks can be obtained by searching in the 2-D MUSIC spectrum [11]. When the range rn lies in the Fresnel region rn∈0.62D3/λ1/2,2D2/λ, through spectrum search, we can get K1 near-field sources, when the range rn exceeds the Fresnel region rf=∞, K−K1 far-field sources we can get through the spectrum search [28].Finally,θ¯k and r¯k are paired one by one, and no operation is required.


### 3.4. Ambiguity Analysis of Range Estimation

In this part, the ambiguity of range needs us to discuss because the ranges of near-field are estimated by ESNA, where some intersensor spacings are more than a half wavelength. In [21, 28, 32], the use of a sparse array to estimate ambiguity problem of range of a given DOA has been discussed. Our analysis is similar to that of [28]. We analyze the ambiguity of the range under the ESNA.Like [28], according to equation (33), the steering vector aθk,rk contains the range parameter rk. If rk is unique, then there is no ambiguity in the steering vector aθk,rk.Similarly, suppose thataθk,rk is ambiguous, from equations (3), (4), and (6), we have(34)πd2pm2cos2θkλrk=πd2pm2cos2θkλr¯k+h2π,where rk denotes the true range and the r¯k denotes the false range, rk,r¯k∈0.62D3/λ1/2,2D2/λ. When the integer h≥1, the range must be ambiguous. Thus, equation (34) can be transformed into(35)pm2≥2λd2cos2θk⋅11/rk−1/r¯k.By minimizing the right-side part of equation (35), we can obtain the maximum range of pm2. Let cos2θk=1, rk=0.62D3/λ1/2, and r¯k=∞, then we have(36)pm2≥2λd20.62D3λ1/2.If we take the array aperture equal toλ and d equal to λ/4, we can get(37)pm≥4.45.which implies that when pm≥5 or pm≤−5, the mth element of the steering vector aθk,r is not unique. However, when D=λ and d=λ/4, pm is less than 5 in the ESNA. Thus, when −4≤pm≤4, the mth element of aθk,r is unique.Consequently, through the analysis of the above, we can obtain that the rangerk is unique through the DOA estimation θk from equation (33).


### 3.5. Computational Complexity

Regarding the computational complexity, the main computation load of the proposed algorithm includes the cumulant matrix and covariance matrices that need to be constructed. For2L+1 sensors and L=M+N. For M=N, the multiplications corresponding to the proposed method are O924M+2M2−1−1T+4M+12T, where T is the snapshot number and the multiplications corresponding to the high-order MUSIC mixed-order MUSIC [20] are O94M+12T+98M+12T and mixed-order MUSIC [21] are O924M+12+1T+4M+12, respectively. Hence, the complexity of the proposed algorithm is the lowest compared with the above two classical algorithms.


## 3.1. Analysis of ESNA in the Cumulant Domain

Cumulant has many advantages, it not only improves the estimation performance but also expands the array aperture [29–31]. In addition, cumulant can also provide an increased number of consecutive lags, which is insensitive to Gaussian noise.Based on the assumptions, we define the fourth-order cumulantcumxmt,xn∗t,xρ∗t,xqt as(11)cumxmt,xn∗t,xρ∗t,xqt=cum∑k=1Ksktejpmγk+pm2ϕk,∑k=1Ksktejpnγk+pn2ϕk∗,∑k=1Ksktejpργk+pρ2ϕk∗,∑k=1Ksktejpqγk+pq2ϕk=∑k=1Kc4,skejpm−pn−pρ−pqγk+pm2−pn2−pρ2−pq2ϕk×cumskt,sk∗t,sk∗t,skt=∑k=1Kc4,skejpm−pn−pρ−pqγk+pm2−pn2−pρ2−pq2ϕk,where c4,sk =  cumskt,sk∗t,sk∗t,skt is the kurtosis of the kth signal and m,n,ρ,q∈−L,L.In equation (11), we remove the ϕk term, while keeping the γk term, it is assumed that n=−m,q=−ρ. Thus, equation (11) becomes(12)cumxmt,x−m∗t,xρ∗t,x−ρt=∑k=1Kc4,skej2pm−pργk.Obviously, equation (12) only contains DOA parameters, which we can estimate DOA of all sources first by constructing a special cumulant matrix.Proposition 1.
LetL=lρ|−M−N≤ρ≤M+N represents the positions of the sensor in the ESNA array, and Lc=lc|lc=pm−pρ represents the lags in the ESNA. Based on the ESNA, we define the following facts to hold:(a)
There are4MN+8N−3 consecutive lags in the range −4N−2MN+2,4N+2MN−2Proof (a): see Appendix.(b)
WhenN=⌈L+2/2⌉ and M=⌊L+2/2⌋, the largest consecutive lags can be obtained, where ⌈⋅⌉ and ⌊⋅⌋ are the round-towards infinite and round-towards zeroProof: we definefN=4NM+8N−3, where N=L−M. f′N=4L−8N+8 can be obtained by deriving N. Notably, N is an integer, when N=⌈L+2/2⌉ or M=⌊L−2/2⌋, f can get the maximum value.As a comparison, four array configurations with2L+1 sensors were considered: (1) symmetric nested array (SNA) [28], (2) symmetric double nested array (SDNA) [23], (3) improved symmetric nested array (ISNA) [10], and (4) the proposed enhanced symmetric nested array (ESNA). In Table 1, in the case of the same number of sensors, we compare the number of consecutive lags and the values of M and N. Among all the nested arrays, the total sensors of the uniform linear array in the middle of a nested array are 2N−1, and the sparse array on both side of a uniform linear array is 2M+2. From Table 1, it is obvious that the proposed ESNA has the maximum number of consecutive lags when the number of array sensors is the same. As the results show, the ESNA can offer more consecutive lags for DOA estimation in mixed source scenario.In order to provide an intuitive impression of above-mentioned four array configurations, herein, we give an example of 13 sensors and mark their positions in Figure2. The red circles represent the lags of four virtual coarrays, and the blue crosses represent the holes for four virtual coarrays. The interval between positions is d.Table 1
Consecutive lags numbers for different arrays.
Physical numberSNASDNAISNAESNAMaximum number of consecutive lagsNMNMNMNMSNASDNAISNAESNA711223130715192192223323111212933112233424117294145133334434232375561153344535231477177174445545339578997Figure 2
Sensor locations.


## 3.2. DOA Estimation of All Sources

In this section, we construct a special cumulant, which can be used to eliminate the range parameter in the near-field, and then use the cumulant to estimate the DOA of the mixed field.In the proposed algorithm, we use the subarray partition method to combine the sparse array and the uniform linear array in the ESNA form, and the corresponding cumulant output and the output cumulant are formed by the largest virtual consecutive lags.Based on equation (11), first, let m=M+N, p∈1−N,…,N−1, and the 2N−1×1 cumulant vector c1 can be obtained, and its N−ρ th element is given by(13)c1N−ρ=cumxM+Nt,x−M−N∗t,xρ∗t,x−ρ=∑k=1Kc4,skej2N−2+2NM+1−ργk.Similarly, since we are using a symmetric enhanced nested array, we can easily obtain a symmetrical part. Letm=−M−N, ρ∈1−N,…,N−1, and we get another 2N−1×1 cumulant vector c¯1, and its N−ρth element is given by(14)c¯1N−ρ=cumx−M−Nt,xM+N∗t,x−ρ∗t,xρt=∑k=1Kc4,skej2−2N−2NM+1−ργk.Second, letm∈N,N+M−1, p∈1−N,…,N−1, we get another M2N−1×1 cumulant vector c2_m, and its 2Nm−N+N−ρ th element is given by(15)c2_m2Nm−N+N−ρ=cumxmt,x−m∗t,x−ρ∗t,xρt=∑k=1Kc4,skej2N−1+2Nm−N+1−ργk.Similarly, letm∈N,N+M−2, ρ=M+N and construct M−1×1 vector c¯2_m, and its m−N+1th element is given by(16)c¯2_mm−N+1=cumxM+Nt,x−M−N∗t,xm∗t,x−mt=∑k=1Kc4,skej22NM+N−m−1γk.Now, we construct the2MN−1×1 vector c2,(17)c2=c2_1,c¯2_M+N−2,…,c2_M+N−2,c¯2_1,c2_M+N−1.Similarly, since we are using a symmetric enhanced nested array, we can easily obtain a symmetrical part. Letm∈1−M−N,−N, ρ∈1−N,…,N−1, we get another M2N−1×1 cumulant vector c2_m∗, and its −2Nm+N+N−ρth element is given by(18)c2_m∗−2Nm+N+N−ρ=cumxmt,x−m∗t,xρ∗t,xρt=∑k=1Kc4,skej21−N+2Nm+N−1−ργk.Similarly, letm∈2−M−N,−N, ρ=−M−N and construct M−1×1 vector c¯2_m∗, and its m+N+2th element is given by(19)c¯2_m∗m+N+2=cumx−M−Nt,xM+N∗t,xm∗t,x−mt=∑k=1Kc4,skej21−2NM+N+mγk.Now, we construct the2MN−1×1 vector c¯2,(20)c¯2=c2_1∗,c¯2_M+N−2∗,…,c2_M+N−2∗,c¯2_1∗,c2_M+N−1∗.Third, letm∈1−N,N−1, ρ=0 and construct the 2N−1×1 vector c3, and its mth element is given by(21)c3m+N=cumxmt,x−m∗t,x0∗t,x0t=∑k=1Kc4,skej2mγk.Fourth, letm=N−1, ρ∈1−N,−1, we get another N−1×1 cumulant vector c4, and its N+ρth element is given by(22)c4N+ρ=cumxmt,x−m∗t,xρ∗t,xρt=∑k=1Kc4,skej2m−ργk.Similarly, since we are using a symmetric enhanced nested array, we can easily obtain a symmetrical part, letm=1−N, ρ∈1,N−1, we construct the N−1×1 vector c¯4, and its ρth element is given by(23)c¯4ρ=cumxmt,x−m∗t,xρ∗t,xρt=∑k=1Kc4,skej2m−ργk.Fifth, letm=M+N, ρ=M+N−1 and m=−M−N, ρ=1−M−N, and we construct the last two cumulants c5 and c¯5,(24)c5=cumxM+Nt,x−M−N∗t,xM+N−1∗t,x1−M−Nt=∑k=1Kc4,skej22N−1γk,and(25)c¯5=cumx−M−Nt,xM+N∗t,x1−M−N∗t,xM+N−1t=∑k=1Kc4,skej21−2Nγk.Finally, letm=M+N−1, ρ=−N, we construct cend, and it is given by(26)cend=cumxM+N−1t,x1−M−N∗t,x−N∗t,xNt=∑k=1Kc4,skej24N+2NM−2γk.Similarly, since we are using a symmetric enhanced nested array, we can easily obtain a symmetrical part. Letm=1−M−N, ρ=N, we construct c¯end, and it is given by(27)c¯end=cumx1−M−Nt,xM+N−1∗t,xN∗t,x−Nt=∑k=1Kc4,skej24N+2NM−2γk.Now, we combine the above vectors to form a24N+2MN−1−1×1 vector c, which is expressed as(28)c=c¯end,c¯1T,c¯2T,c¯5,c¯4T,c3T,c4T,c5,c2T,c1T,cendT.A Toeplitz matrixC is constructed through the vector c, and matrix C is 4N+2MN−1×4N+2MN−1, whose mth column of C consists of the 4N+2MN−1+1−mth to 24N+2MN−1+1−mth element of c. The Toeplitz matrix C can be expressed in a compact form as(29)C=AθCsAHθ,Cs=diagc4s1,…,c4sK,Aθ=aθ1,…,aθK,and(30)aθk=1,ej2γ1,…,ej24N+2MN−2γkT.Since the sources signals have nonzero fourth-order cumulant, the matrixC is full rank. Therefore, we estimate the common θ of mixed sources by performing the eigenvalue decomposition of C.Based on the EVD ofC, using the MUSIC algorithm [6] to estimate DOA and we can get the following formula:(31)θk=minaθHUnUnHaθ,k=1,2,…,K,where Un represents the noise subspace of C.


## 3.3. Range Estimation and Source Type Classification

Onceθ¯k,k=1,…,K is estimated, then we can estimate rk. When we get the convariance matrix R=ExtxHt, perform eigen-decomposition on it, and we can get(32)R=U¯sV¯sU¯sH+U¯nV¯nU¯nH,where U¯n spans the noise subspace, which contains 4N+2MN−1−K eigenvectors. By taking the estimated θ¯k into aθk,r, the rk can be obtained by the following:(33)rk=minaθk,rU¯nU¯nHaHθk,r,k=1,…,K.Obviously,K peaks can be obtained by searching in the 2-D MUSIC spectrum [11]. When the range rn lies in the Fresnel region rn∈0.62D3/λ1/2,2D2/λ, through spectrum search, we can get K1 near-field sources, when the range rn exceeds the Fresnel region rf=∞, K−K1 far-field sources we can get through the spectrum search [28].Finally,θ¯k and r¯k are paired one by one, and no operation is required.


## 3.4. Ambiguity Analysis of Range Estimation

In this part, the ambiguity of range needs us to discuss because the ranges of near-field are estimated by ESNA, where some intersensor spacings are more than a half wavelength. In [21, 28, 32], the use of a sparse array to estimate ambiguity problem of range of a given DOA has been discussed. Our analysis is similar to that of [28]. We analyze the ambiguity of the range under the ESNA.Like [28], according to equation (33), the steering vector aθk,rk contains the range parameter rk. If rk is unique, then there is no ambiguity in the steering vector aθk,rk.Similarly, suppose thataθk,rk is ambiguous, from equations (3), (4), and (6), we have(34)πd2pm2cos2θkλrk=πd2pm2cos2θkλr¯k+h2π,where rk denotes the true range and the r¯k denotes the false range, rk,r¯k∈0.62D3/λ1/2,2D2/λ. When the integer h≥1, the range must be ambiguous. Thus, equation (34) can be transformed into(35)pm2≥2λd2cos2θk⋅11/rk−1/r¯k.By minimizing the right-side part of equation (35), we can obtain the maximum range of pm2. Let cos2θk=1, rk=0.62D3/λ1/2, and r¯k=∞, then we have(36)pm2≥2λd20.62D3λ1/2.If we take the array aperture equal toλ and d equal to λ/4, we can get(37)pm≥4.45.which implies that when pm≥5 or pm≤−5, the mth element of the steering vector aθk,r is not unique. However, when D=λ and d=λ/4, pm is less than 5 in the ESNA. Thus, when −4≤pm≤4, the mth element of aθk,r is unique.Consequently, through the analysis of the above, we can obtain that the rangerk is unique through the DOA estimation θk from equation (33).


## 3.5. Computational Complexity

Regarding the computational complexity, the main computation load of the proposed algorithm includes the cumulant matrix and covariance matrices that need to be constructed. For2L+1 sensors and L=M+N. For M=N, the multiplications corresponding to the proposed method are O924M+2M2−1−1T+4M+12T, where T is the snapshot number and the multiplications corresponding to the high-order MUSIC mixed-order MUSIC [20] are O94M+12T+98M+12T and mixed-order MUSIC [21] are O924M+12+1T+4M+12, respectively. Hence, the complexity of the proposed algorithm is the lowest compared with the above two classical algorithms.


## 4. Simulation Results

In the last section, the performance of the proposed array is evaluated by simulation results, which are compared with the other three array geometries, SNA [28], SDNA [29], and ISNA [21], and we also compared with HOS MUSIC [20] and mixed-order MUSIC [21] algorithms to demonstrate the performance of the proposed algorithm. The total number of elements is Q=13. The sources signal is ejζt, where the phases ζt are uniformly distributed in 0,2π. However, the four arrays have different array apertures; therefore, they have different Fresnel regions. Here, taken the common Fresnel region for NF region 29λ<r<70λ to all arrays. The performance of the four arrays and proposed algorithm are evaluated through the root-mean-square error (RMSE) curve, and the number of Monte-Carlo trials is Ns=500.(38)RMSE=1NsK∑k=1K∑ns=1Nsα¯km−αk2,where αk represents the DOA θk or the range rk and α¯k is the estimate of αk for the mth trial. The resolution probability can be defined as if the estimated DOA of two sources satisfies θ¯k−θk<ε/2, where ε is the difference between the two angles, it is considered that it can successfully distinguish two angles.


### 4.1. Experiment 1

In the first simulation, we used two NF sources located at20°,30λ,50°,35λ and two FF sources located at −40°,∞,20°,∞, respectively. The SNR is fixed at 15 dB, and the number of snapshots is set at 2000. The spatial spectrum of mixed sources is drawn in Figure 3. We can observe from Figure 3(a) that a total of four obvious peaks are estimated, two of which are near-field sources and two are far-field sources. We can also see from Figure 3(b) that there aretwo peaks located in the NF region and another two peaks in the FF region. According to Figure 3, the proposed algorithm can effectively distinguish the near-field and far-field in the mixed field.Figure 3
Spatial spectrum of mixed sources. SNR = 15 dB andT = 2000 snapshots. (a) Spatial spectrum of the NF and FF estimators. (b) Spatial spectrum of range estimators.
(a)(b)


### 4.2. Experiment 2

In the second simulation, we compare the probability of arrays for two closely spaced sources against SNR in Figure4. In Figure 4(a), two pure NF sources are located at 40°,31λ,44°,32λ, and the angle interval ε of the two is 4°. In Figure 4(b), two pure FF sources are located at −20°,∞,−16°,∞, and the angle interval ε of the two is 4°. In Figure 4(c), two mixed sources and one NF source are located at 34°,31λ, another FF source is located at 30°,∞, and the angle interval ε of the two is also 4°. When T=2000 snapshots, and SNR ∈−6,24 dB. From the simulation result, compared with other arrays, the proposed ESNA can obtain better estimation accuracy under the same conditions, and it can successfully distinguish the close sources at lower SNR.Figure 4
Probability versus SNR for closely spaced sources. (a) Near-field. (b) Far-field. (c) Mixed sources.
(a)(b)(c)


### 4.3. Experiment 3

In the third simulation, we used two pure near-field sources located at40°,31λ and 50°,32λ, respectively:(1)
RMSE versus SNR: when snapshots = 2000, and SNR ∈−5,30 dB, from Figure 5, we can see the RMSE curves of DOA and range. As can be seen from Figures 5(a) and 5(b), in all SNRs, the ESNA can provide better estimation accuracy than the other arrays. As is expected, when the SNR increases, the RMSE curves of DOA and range show a downward trend.(2)
RMSE versus snapshots: when SNR = 15 dB, and snapshots ∈50,2000, from Figure 6, we can see the RMSE curves of DOA and range. From Figures 6(a) and 6(b), we can see that the proposed algorithm still has a better performance than the other algorithms for all available snapshots. This is also because the ESNA can obtain more consecutive lags compared with other arrays. As is expected, with the increase of snapshot, the RMSE curves show a downward trend.Figure 5
RMSE curves of DOA and range estimation versus SNR. (a) DOA estimations for near-field. (b) Range estimations for near-field.
(a)(b)Figure 6
RMSE curves of DOA and range estimation versus snapshots. (a) DOA estimations for near-field. (b) Range estimations for near-field.
(a)(b)


### 4.4. Experiment 4

In the fourth simulation, we used mixed sources. One near-field source is located at30°,30λ, and another far-field source is located at −10°,∞:(1)
RMSE versus SNR: when snapshots = 2000, and SNR ∈−5,30 dB, from Figure 7, we can see the RMSE curves of DOA and range about mixed sources (range estimations are not performed for far-field source). As can be seen from Figures 7(a) and 7(b), in all SNRs, the ESNA can obtain better estimation accuracy than the other arrays. As is expected, when the SNR increases, the RMSE curves of DOA and range show a downward trend.(2)
RMSE versus snapshots: when SNR = 15 dB, and snapshots ∈50,2000, from Figure 8, we can see the RMSE curves of DOA and range. As can be seen from Figures 8(a) and 8(b), in all snapshots, the ESNA can obtain better estimation accuracy than the other arrays. As is expected, with the increase of snapshot, the RMSE of DOA and range show a downward trend.Figure 7
RMSE curves of DOA and range estimation versus SNR. (a) DOA estimations for NF and FF. (b) Range estimations for NF.
(a)(b)Figure 8
RMSE curves of DOA and range estimation versus snapshots. (a) DOA estimations for NF and FF. (b) Range estimations for NF.
(a)(b)


### 4.5. Experiment 5

In the fifth simulation, we used to deal with two pure far-field sources. Two FF sources are located at40°,∞,50°,∞, respectively:(1)
RMSE versus SNR: when snapshots = 2000, and SNR ∈−5,30 dB, from Figure 9, we can see the RMSE curves of DOA and range. As can be seen from Figure 9(a), in all SNRs, the ESNA can obtain better estimation accuracy than the other arrays. This is because the ESNA can obtain more consecutive lags compared with other arrays. As is expected, when the SNR increases, the RMSE of DOA and range show a downward trend.(2)
RMSE versus snapshots: when SNR = 15 dB, and snapshots ∈50,2000, from Figure 9, we can see the RMSE curves of DOA and range. As can be seen from Figure 9(b), we can see that the proposed algorithm still has a better performance than the other algorithms for all available snapshots. This is also because the ESNA can obtain more consecutive lags compared with other arrays. As is expected, with the increase of snapshot, the RMSE of DOA and range show a downward trend.Figure 9
RMSE curves of DOA estimation for pure FF versus SNR and snapshots. (a) DOA estimations for pure FF. (b) DOA estimations for pure FF.
(a)(b)


### 4.6. Experiment 6

In the last simulation, we present simulations to prove that our proposed algorithm outperforms HOS MUSIC and mixed-order MUSIC algorithms. We used mixed sources. One near-field source is located at20°,31λ, and another far-field source is located at 30°,∞:(1)
RMSE versus SNR: when snapshots = 2000, and SNR ∈−5,30 dB, from Figure 10, we can see the RMSE curves of DOA (range estimation is not performed for mixed-field source). As can be seen from Figure 10(a), in all SNRs, the proposed algorithm can obtain better estimation accuracy than the other algorithms. As is expected, when the SNR increases, the RMSE curves of DOA show a downward trend.(2)
RMSE versus snapshots: when SNR = 15 dB, and snapshots ∈50,2000, from Figure 10, we can see the RMSE curves of DOA (range estimation is not performed for mixed-field source). As can be seen from Figure 10(b), in all snapshots, the proposed algorithms can obtain better estimation accuracy than the other algorithms. As is expected, with the increase of snapshot, the RMSE of DOA show a downward trend.(3)
Computational complexity versus array sensors: when array sensors varied from 13 to 25 in 4 steps, from Figure11(a), we can see that the computational complexity increases with the increase of the number of array sensors, but the computational complexity of the proposed algorithm is still the lowest.(4)
Estimate sources versus array sensors: when snapshots = 250 000, and SNR = 10 dB, fifteen sources uniformly distributed from−70° to 70° were considered. From Figure 11(b), we can see that our proposed array was suitable for underdetermined cases. We used 13 sensors to estimate 15 sources.(5)
DOF comparison versus other arrays: when the array sensors varied from 7 to 55 in 2 steps, from Figure11(c), we compare the number of DOF achieved from the four arrays. We can see that the number of DOF increases with array sensors, and proposed array provided a higher number of DOF than other arrays.Figure 10
RMSE curves of DOA estimation for pure FF versus SNR and snapshots. (a) DOA estimations for SNR. (b) DOA estimations for snapshots.
(a)(b)Figure 11
Computational complexity, underdetermined estimate and DOF comparison. (a) Computation complexity for array sensors. (b) Estimate sources versus array sensors. (c) DOF comparison versus different arrays.
(a)(b)(c)


## 4.1. Experiment 1

In the first simulation, we used two NF sources located at20°,30λ,50°,35λ and two FF sources located at −40°,∞,20°,∞, respectively. The SNR is fixed at 15 dB, and the number of snapshots is set at 2000. The spatial spectrum of mixed sources is drawn in Figure 3. We can observe from Figure 3(a) that a total of four obvious peaks are estimated, two of which are near-field sources and two are far-field sources. We can also see from Figure 3(b) that there aretwo peaks located in the NF region and another two peaks in the FF region. According to Figure 3, the proposed algorithm can effectively distinguish the near-field and far-field in the mixed field.Figure 3
Spatial spectrum of mixed sources. SNR = 15 dB andT = 2000 snapshots. (a) Spatial spectrum of the NF and FF estimators. (b) Spatial spectrum of range estimators.
(a)(b)


## 4.2. Experiment 2

In the second simulation, we compare the probability of arrays for two closely spaced sources against SNR in Figure4. In Figure 4(a), two pure NF sources are located at 40°,31λ,44°,32λ, and the angle interval ε of the two is 4°. In Figure 4(b), two pure FF sources are located at −20°,∞,−16°,∞, and the angle interval ε of the two is 4°. In Figure 4(c), two mixed sources and one NF source are located at 34°,31λ, another FF source is located at 30°,∞, and the angle interval ε of the two is also 4°. When T=2000 snapshots, and SNR ∈−6,24 dB. From the simulation result, compared with other arrays, the proposed ESNA can obtain better estimation accuracy under the same conditions, and it can successfully distinguish the close sources at lower SNR.Figure 4
Probability versus SNR for closely spaced sources. (a) Near-field. (b) Far-field. (c) Mixed sources.
(a)(b)(c)


## 4.3. Experiment 3

In the third simulation, we used two pure near-field sources located at40°,31λ and 50°,32λ, respectively:(1)
RMSE versus SNR: when snapshots = 2000, and SNR ∈−5,30 dB, from Figure 5, we can see the RMSE curves of DOA and range. As can be seen from Figures 5(a) and 5(b), in all SNRs, the ESNA can provide better estimation accuracy than the other arrays. As is expected, when the SNR increases, the RMSE curves of DOA and range show a downward trend.(2)
RMSE versus snapshots: when SNR = 15 dB, and snapshots ∈50,2000, from Figure 6, we can see the RMSE curves of DOA and range. From Figures 6(a) and 6(b), we can see that the proposed algorithm still has a better performance than the other algorithms for all available snapshots. This is also because the ESNA can obtain more consecutive lags compared with other arrays. As is expected, with the increase of snapshot, the RMSE curves show a downward trend.Figure 5
RMSE curves of DOA and range estimation versus SNR. (a) DOA estimations for near-field. (b) Range estimations for near-field.
(a)(b)Figure 6
RMSE curves of DOA and range estimation versus snapshots. (a) DOA estimations for near-field. (b) Range estimations for near-field.
(a)(b)


## 4.4. Experiment 4

In the fourth simulation, we used mixed sources. One near-field source is located at30°,30λ, and another far-field source is located at −10°,∞:(1)
RMSE versus SNR: when snapshots = 2000, and SNR ∈−5,30 dB, from Figure 7, we can see the RMSE curves of DOA and range about mixed sources (range estimations are not performed for far-field source). As can be seen from Figures 7(a) and 7(b), in all SNRs, the ESNA can obtain better estimation accuracy than the other arrays. As is expected, when the SNR increases, the RMSE curves of DOA and range show a downward trend.(2)
RMSE versus snapshots: when SNR = 15 dB, and snapshots ∈50,2000, from Figure 8, we can see the RMSE curves of DOA and range. As can be seen from Figures 8(a) and 8(b), in all snapshots, the ESNA can obtain better estimation accuracy than the other arrays. As is expected, with the increase of snapshot, the RMSE of DOA and range show a downward trend.Figure 7
RMSE curves of DOA and range estimation versus SNR. (a) DOA estimations for NF and FF. (b) Range estimations for NF.
(a)(b)Figure 8
RMSE curves of DOA and range estimation versus snapshots. (a) DOA estimations for NF and FF. (b) Range estimations for NF.
(a)(b)


## 4.5. Experiment 5

In the fifth simulation, we used to deal with two pure far-field sources. Two FF sources are located at40°,∞,50°,∞, respectively:(1)
RMSE versus SNR: when snapshots = 2000, and SNR ∈−5,30 dB, from Figure 9, we can see the RMSE curves of DOA and range. As can be seen from Figure 9(a), in all SNRs, the ESNA can obtain better estimation accuracy than the other arrays. This is because the ESNA can obtain more consecutive lags compared with other arrays. As is expected, when the SNR increases, the RMSE of DOA and range show a downward trend.(2)
RMSE versus snapshots: when SNR = 15 dB, and snapshots ∈50,2000, from Figure 9, we can see the RMSE curves of DOA and range. As can be seen from Figure 9(b), we can see that the proposed algorithm still has a better performance than the other algorithms for all available snapshots. This is also because the ESNA can obtain more consecutive lags compared with other arrays. As is expected, with the increase of snapshot, the RMSE of DOA and range show a downward trend.Figure 9
RMSE curves of DOA estimation for pure FF versus SNR and snapshots. (a) DOA estimations for pure FF. (b) DOA estimations for pure FF.
(a)(b)


## 4.6. Experiment 6

In the last simulation, we present simulations to prove that our proposed algorithm outperforms HOS MUSIC and mixed-order MUSIC algorithms. We used mixed sources. One near-field source is located at20°,31λ, and another far-field source is located at 30°,∞:(1)
RMSE versus SNR: when snapshots = 2000, and SNR ∈−5,30 dB, from Figure 10, we can see the RMSE curves of DOA (range estimation is not performed for mixed-field source). As can be seen from Figure 10(a), in all SNRs, the proposed algorithm can obtain better estimation accuracy than the other algorithms. As is expected, when the SNR increases, the RMSE curves of DOA show a downward trend.(2)
RMSE versus snapshots: when SNR = 15 dB, and snapshots ∈50,2000, from Figure 10, we can see the RMSE curves of DOA (range estimation is not performed for mixed-field source). As can be seen from Figure 10(b), in all snapshots, the proposed algorithms can obtain better estimation accuracy than the other algorithms. As is expected, with the increase of snapshot, the RMSE of DOA show a downward trend.(3)
Computational complexity versus array sensors: when array sensors varied from 13 to 25 in 4 steps, from Figure11(a), we can see that the computational complexity increases with the increase of the number of array sensors, but the computational complexity of the proposed algorithm is still the lowest.(4)
Estimate sources versus array sensors: when snapshots = 250 000, and SNR = 10 dB, fifteen sources uniformly distributed from−70° to 70° were considered. From Figure 11(b), we can see that our proposed array was suitable for underdetermined cases. We used 13 sensors to estimate 15 sources.(5)
DOF comparison versus other arrays: when the array sensors varied from 7 to 55 in 2 steps, from Figure11(c), we compare the number of DOF achieved from the four arrays. We can see that the number of DOF increases with array sensors, and proposed array provided a higher number of DOF than other arrays.Figure 10
RMSE curves of DOA estimation for pure FF versus SNR and snapshots. (a) DOA estimations for SNR. (b) DOA estimations for snapshots.
(a)(b)Figure 11
Computational complexity, underdetermined estimate and DOF comparison. (a) Computation complexity for array sensors. (b) Estimate sources versus array sensors. (c) DOF comparison versus different arrays.
(a)(b)(c)


## 5. Conclusion

In this paper, we proposed an enhanced symmetric nested array (ESNA) to solve localization problems about mixed near-field and far-field sources; we also proposed a novel subarray partition method to construct a special cumulant to eliminate the range parameter in the near-field compared with the existing cumulant. By using the ESNA, the effective number of consecutive lags and virtual aperture are increased greatly compared with other arrays. Using a special cumulant, the conventional MUSIC method is used to estimate the mixed sources, and the computation complexity is moderate. Numerical results show that the proposed array has higher resolution ability and provides the improved accuracy of parameter estimation in almost all SNRs levels, and snapshot numbers are compared with the existing arrays. Therefore, the proposed array evidently outperforms the existing arrays for scenarios such as pure NF, pure FF, and mixed sources.



---

*Source: 3960361-2021-12-20.xml*