Research Article  Open Access
Efficient Localization Algorithm for NearField Noncircular Sources via DualPolarization Sensor Array
Abstract
Noncircular signals are widely used in the area of radar, sonar, and wireless communication array systems, which can offer more accurate estimates and detect more sources. In this paper, the noncircular signals are employed to improve source localization accuracy and identifiability. Firstly, an extended realvalued covariance matrix is constructed to transform complexvalued computation into realvalued computation. Based on the property of noncircular signals and symmetric uniform linear array (SULA) which consist of dualpolarization sensors, the array steering vectors can be separated into the source position parameters and the nuisance parameter. Therefore, the rank reduction (RARE) estimators are adopted to estimate the source localization parameters in sequence. By utilizing polarization information of sources and realvalued computation, the maximum number of resolvable sources, estimation accuracy, and resolution can be improved. Numerical simulations demonstrate that the proposed method outperforms the existing methods in both resolution and estimation accuracy.
1. Introduction
Passive source localization is a key problem in array signal processing for applications such as radar, sonar, microphone arrays, and communication [1]. In recent years, it has received more concern and has developed lots of methods to deal with this issue. Among them, the most typical algorithms are multiple signal classification (MUSIC) [2], estimating signal parameter via rotational invariance techniques (ESPRIT) [3], and their derivatives. Nevertheless, these methods are under the assumption that the sources are in the farfield (FF), which means the wavefronts are plane waves and therefore only directionofarrival (DOA) parameters are required to be estimated.
However, when the radiating sources are located in the nearfield (NF) of the array, whose wavefronts are spherical waves, both DOA and range parameters should be estimated to localize the sources. Thus, the traditional FF algorithms are no longer suitable for NF sources. Fortunately, many advanced algorithms have been presented under NF assumption. Huang and Barkat proposed a twodimensional (2D) MUSIC method in the anglerange domain to achieve NF superresolution localization, but the 2D joint search brought high computational complexity [4]. Challa and Shamsunder took the lead in introducing highorder cumulant into NF source parameter estimation problem. By constructing multiple cumulant matrixes, they proposed ESPRITLike to estimate DOA and range parameters of NF sources [5]. However, the ESPRITLike algorithm had many drawbacks, such as high computational complexity, parameter pairing, and array aperture loss. Lee et al. proposed a covariance approximation method (CA) [6]. The method reconstructed the elements of the NF covariance matrix, so that the NF source was converted into a virtual FF source (the DOA information was the same), and the traditional FF direction finding methods could apply to the DOA estimation of FF sources, avoiding multidimensional search. But the CA algorithm would produce an image source in the case of coherent sources. Noh and Lee analyzed the phenomenon and proposed a method to suppress the image source effectively [7]. Grosicki et al. proposed the weighted linear prediction (WLP) algorithm to obtain the DOA and range estimation [8]. Utilizing the symmetric linear array, Zhi and Chia proposed the classical generalized ESPRIT algorithm [9].
In recent years, many directionfinding methods which employ the noncircular signals received more concern, such as Binary Phase Shift Keying (BPSK), Pulse Amplitude Modulation (PAM), and Amplitude Shift Keying (ASK) signals. By taking use of the noncircularity of the signal, the array can benefit from the extended virtual aperture, which means that the resolution capability and the estimation accuracy can be improved. Chen et al. considered the DOA estimation of noncircular signal for uniform linear array via the propagation method and Euler transformation [10]. Tan et al. proposed a weighted unitary nuclear norm minimization approach for DOA estimation in the strictly noncircular sources case [11]. Xie et al. proposed a realvalued localization algorithm for noncircular signals using the uniform linear array [12]. Furthermore, Xie et al. proposed another nearfield localization method for noncircular sources via generalized ESPRIT [13]. Chen et al. proposed a novel localization method for NF rectilinear or strictly noncircular sources with a symmetric uniform linear array of cocentered orthogonal loop and dipole (COLD) antennas [14].
However, most of the abovementioned algorithms use the scalar sensors array, which cannot exploit the polarization information embedded in the electromagnetic waves. An array of vectors sensors can detect signals by utilizing the polarization diversity. For this reason, electromagnetic vector sensors array signal processing has attracted much attention in recent years. Obeidat et al. proposed polarization ESPRITLike algorithm by using polarization sensitive array [15]. However, it suffers from half aperture loss. To avoid the aperture loss problem, Wu et al. developed a least squaresvirtual ESPRIT algorithm (LSVESPA) [16]. But it involves extra parameter pairing procedure. Based on the symmetric sparse linear array with dualpolarization sensors, Tao et al. proposed the Fresnelregion rank reduction (FRRARA) algorithm [17] that enhanced array aperture and only required secondorder statistics. He et al. presented a NF localization of partially polarization sources with a crossdipole array [18].
In this paper, we construct an augmented covariance matrix which consists of the real part and imaginary part of array outputs data. Then, based on the noncircularity of signals and the property of symmetric uniform linear array (SULA), the array steering vector could be decoupled as the product of three realvalued matrixes including DOA, range, and other nuisance parameters, respectively. Consequently, a rank reduction (RARE) based localization method is derived, which translates multidimensional spectral search into multiple onedimensional (1D) spectral searches. The proposed method has the following advantages: (1) As a result of the exploitation of noncircularity, more sources can be resolved. (2) It is efficient since it avoids exhaustive complexvalued computation and multidimensional search. (3) The estimation accuracy and resolution are improved effectively by utilizing the noncircularity and polarization diversity.
The rest of this paper is organized as follows. In Section 2, the data model for NF noncircular signals which received by dualpolarization sensors SULA is formulated. The proposed localization algorithm for NF noncircular sources is developed in Section 3. In Section 4, we discuss the performance of the proposed algorithm and some newly developed algorithms. Then, numerical simulations are presented in Section 5. Conclusion is drawn in Section 6.
Notations: The transpose, conjugate, and conjugate transpose are denoted by, , and , respectively. The symbol represents the Kronecker product. , , and symbolize the real part operator, the imaginary part operator, and the determinant of a matrix.
2. Signal Model
We suppose that K independent NF noncircular signals impinge upon a SULA as shown in Figure 1. The array is composed of dualpolarization sensors which is placed along the yaxis, and its sensors position is , where d is the interelement spacing. The dualpolarization sensors used in this paper is crossdipole. This localization algorithm is achievable by other polarization sensors, such as cocentered orthogonal loop and dipole pair [19]. But the crossdipole is very small and easier to design, so it is more common in practice.
Note that the two polarization components of the crossdipole point to the xaxis and yaxis directions, respectively. Assuming that all sources are located in the yz plane, then the two direction components of electronic field can be described aswhere α denotes the auxiliary polarization angle, β represents the polarization phase difference, and θ signifies the DOA of the signal.
Let the array center (sensor 0) be the reference point; the output signal components in xpolarization and ypolarization received by the mth sensor at time t can be modeled aswhere denotes the kth signal, , represents the phase shift related to the kth signal’s propagation time delay from the reference point to the mth sensor, , , and are the auxiliary polarization angle, polarization phase difference, and DOA of kth signal, and and symbolize the additive noise.
Consider that the sources are located in the NF, the time delay can be approximated as [4]where , , λ denotes the wavelength of signal, and represents the range parameter of the kth source.
When the source signals are noncircular signals, it can be obtained that ; herein, is the noncircular phase of the kth signals and is the zerophased version of the source signal.
In a matrix form, (2) and (3) can be written aswhere
In the above equations, is array output vector, represents the array steering matrix, symbolizes the source signal vector, and denotes the additive noise matrix.
With a total of L snapshots taken at the distinct instants , the problem is to determine the localization parameters from the array output data. Throughout this paper, the following hypotheses are assumed to hold:(1)The incoming signals are mutually independent, narrowband stationary noncircular(2)The sensor noise is additive zeromean white Gaussian and independent of the source signals(3)In order to avoid the phase ambiguity, the intersensor spacing d should be within a quarter wavelength(4)The range parameter r lies in which means the signal sources lies in the NF [20], where D is the aperture of array
3. Proposed Algorithm
The traditional subspace estimation algorithm requires multidimensional spectral peak search that exhausts high computational burden. And to the best of our knowledge, there has been very limited work utilizing polarization diversity and the noncircularity of signals simultaneously. In order to overcome these shortages and improve estimation performance, we construct a realvalued augmented output matrix based on the Euler equation. Then, the steering vector is factorized with respect to the localization parameters and nuisance parameters. Based on the RARE criterion [21] and 1D search, the localization parameters can be estimated. Consequently, the multidimensional optimization problem could be accomplished by realvalued computation and 1D spectral searches.
3.1. RealValued Augmented Covariance Matrix
In order to transform the complexvalued data into realvalued domain, we achieve the real part and imaginary part of , respectively, by the following equations:
According to the signal model in (5), we can construct the realvalued augmented data matrix aswhere is the augmented steering matrix and is the augmented noise matrix. And the th row of can be expressed as
The elements in is formulated in (10), where :
Then, the realvalued augmented covariance matrix of can be expressed as . By taking the eigendecomposition (EVD) of , we havewhere and are the diagonal matrixes which contain K largest eigenvalues and other smallest eigenvalues. is a matrix spanning the signal subspace of , and is a matrix spanning the noise subspace of ,
3.2. Joint DOA, Range, and Polarization Estimation
According to the principle of MUSIC, the parameters can be estimated by multidimensional searching of the following spectrum function:
It is obvious that K sets of parameters can be achieved by finding the K minimum value of . However, the estimator in (12) is computationally intensive since it requires 2D spectral search. To avoid the highdimensional computation, the augmented steering vector can be factorized, leading to a simple 1D operation.
Due to the symmetric property of array steering vector, the augmented steering vector can be decoupled as the following formation:whereHerein,
If and only if , the following equation holds:
Equation (20) can be rewritten aswhere and .
It is noteworthy that is only related to θ and independent with other parameters. Since and , , so if , is generally of full column rank. Thus, only if , (21) becomes true because the rank of drops according to the RARE criterion. Therefore, the DOAs of sources can be estimated by the following spectrum function:
After finding K peaks of through a 1D search, the DOAs of sources are obtained. Therefore, the range parameters of sources can be achieved via another RARE estimator.
Define , then (20) can be expressed as
Similarly, since , (23) can hold true if drops rank. Therefore, with obtained , the range of sources can be estimated by of the following 1D spectrum searching function for K times:
3.3. Implementation of the Algorithm
Note that the exact covariance matrix and subspaces are utilized in the previous sections, but the theoretical covariance matrix is unavailable due to the limited snapshots. In practice, it can be estimated as
To summarize, the procedures of the proposed method are shown as follows:(1)Reconstruct the realvalued augmented data matrix based on array output matrix and Euler equation.(2)Take eigendecomposition operation of the covariance matrix , and obtain noise subspace matrix .(3)Compute rank reduction matrix , and estimate DOAs of signals by searching the K highest peaks of (22).(4)Compute rank reduction matrix , with obtained DOAs estimate the range parameter by searching the highest peak of (24), repeat this step from to K. If the polarization parameters of signals need to be estimated, do step 5.(5)Decouple the vector like (13) further and separate the noncircular parameter from polarization; the polarization can be obtained by another RARE estimator.(6)Insert estimated , , , and into (12). By searching the highest peaks of , the noncircular parameter is estimated. Repeat this step from to K.
4. Discussion
4.1. Maximum Number of Resolvable Sources
In this part, we discuss the maximum numbers of resolvable sources of GESPRIT [9], FRRARE [17], and the proposed method, respectively. To facilitate the analysis, a SULA is assumed to have elements. Since GESPRIT brings half aperture loss, it can handle M sources at most. Because the subspacebased algorithm needs at least one eigenvector to span noise space, FRRARE can resolve up to sources. For the proposed method, the noncircularity has been utilized to construct an extended subspace; hence, it can estimate sources at most, which is doubled, compared with FRRARE.
4.2. Computational Complexity
In this part, only the major computation complexity is considered, such as construction of statistical matrices, eigenvalue decomposition (EVD), and spectral search. The search stepsizes for the angle parameter and range parameter are denoted as and . We assume that the number of an array is N and the number of snapshots equals L. The GESPRIT algorithm requires the construction of two secondorder covariance matrices, two EVDs, and twice spectral searches for DOA and range estimation. FRRARE builds a secondorder covariance matrices, performs the EVD, and twice spectral searches for DOA and range estimation. For the proposed one, it needs to establish a secondorder realvalued covariance matrix, perform EVD on this matrix, and execute 1D spectral search (Table 1).

5. Numerical Simulation
In this section, numerical simulations are conducted to validate the performance of the proposed algorithm. Without loss of generality, we consider a uniform linear symmetric array composed of 5 dualpolarization sensors with the interelement spacing being a quarterwavelength. The NF signal sources impinged upon the array are equipowered, statistically independent BPSK signal (code rate = 0.1/Ts, noncircularity = 1). Moreover, the estimation performance is measured by the root meansquare error (RMSE) of independent 500 Monte Carlo trials. The RMSE is defined aswhere is the exact DOA or the range , and denotes the estimation of .
In the first experiment, we suppose that five BPSK equipower signals are located at , , , , and , respectively. The snapshot number and signaltonoise ratio (SNR) are set as 500 and . The DOA and range spectrums of proposed method are shown in Figures 2 and 3. And one can observe that all the sources have been resolved effectively. The DOA spectrums of proposed method, GESPRIT [9], and FRRARE [17] are plotted in Figure 4. By employing realvalued computations, the proposed method can resolve up to sources while FRRARE can only handle sources. Since GESPRIT brings half aperture loss, it can handle M sources at most.
In the second experiment, we investigate RMSEs of the proposed method, GESPRIT and FRRARE with the variation of SNR. We consider two BPSK signals are located at and . The snapshots number is fixed at 500 and the SNR varies from to in a stepsize of . The RMSEs versus SNR are illustrated in Figures 5 and 6. It is obvious that the estimation performance of two algorithms improves as the SNR increases. Moreover, as a result of higher degrees of freedom (DOFs), the estimation accuracy of the proposed method is higher than GESPRIT and FRRARE. Furthermore, there is a significant improvement of algorithm performance when the two sources are close to each other, which means the resolution of the proposed method is higher.
In the third experiment, we study the estimation performance with the variation of snapshots. The parameter settings of two signals are the same as that of the second experiment. The SNR is set to , and the number of snapshots varies from 1 to 1000. Figures 7 and 8 leads to a similar conclusion as in the second experiment that the RMSEs decrease with the increasing number of snapshots. This is because that larger samples will provide a better estimate of the covariance matrix for stationary data.
6. Conclusion
In this paper, a novel localization algorithm for the nearfield noncircular signals is presented by employing the realvalued computation and 1D search. The proposed method utilizes the polarization information and noncircularity, which improves the estimation performance significantly. Compared with some existing works, the proposed method has achieved more resolvable signals and improved estimation accuracy and resolution. The simulation results demonstrate the efficiency and effectiveness of the proposed method for the localization of noncircular sources in nearfield.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was supported by the Fundamental Research Funds for the Universities (No. BDY06) and Innovation Project of Science of Technology Commission of the Central Military Commission (No. ∗∗H863∗∗XJ001∗∗02).
References
 H. Krim and M. Viberg, “Two decades of array signal processing research: the parametric approach,” IEEE Signal Processing Magazine, vol. 13, no. 4, pp. 67–94, 1996. View at: Publisher Site  Google Scholar
 R. Schmidt, “Multiple emitter location and signal parameter estimation,” IEEE Transactions on Antennas and Propagation, vol. 34, no. 3, pp. 276–280, 1986. View at: Publisher Site  Google Scholar
 R. Roy and T. Kailath, “ESPRITestimation of signal parameters via rotational invariance techniques,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 37, no. 7, pp. 984–995, 1989. View at: Publisher Site  Google Scholar
 Y.D. Huang and M. Barkat, “Nearfield multiple source localization by passive sensor array,” IEEE Transactions on Antennas and Propagation, vol. 39, no. 7, pp. 968–975, 1991. View at: Publisher Site  Google Scholar
 R.N. Challa and S. Shamsunder, “Highorder subspacebased algorithms for passive localization of nearfield sources,” in Proceedings of the Conference Record of the TwentyNinth Asilomar Conference on Signals, Systems and Computers, vol. 2, pp. 777–781, Pacific Grove, CA, USA, October 1995. View at: Publisher Site  Google Scholar
 J.H. Lee, Y. Chen, and C.C. Yeh, “A covariance approximation method for nearfield directionfinding using a uniform linear array,” IEEE Transactions on Signal Processing, vol. 43, no. 5, pp. 1293–1298, 1995. View at: Publisher Site  Google Scholar
 H. Noh and C. Lee, “A covariance approximation method for nearfield coherent sources localization using uniform linear array,” IEEE Journal of Oceanic Engineering, vol. 40, no. 1, pp. 187–195, 2015. View at: Publisher Site  Google Scholar
 E. Grosicki, K. AbedMeraim, and Y. Yingbo Hua, “A weighted linear prediction method for nearfield source localization,” IEEE Transactions on Signal Processing, vol. 53, no. 10, pp. 3651–3660, 2005. View at: Publisher Site  Google Scholar
 W. Zhi and M. Y.W. Chia, “Nearfield source localization via symmetric subarrays,” IEEE Signal Processing Letters, vol. 14, no. 6, pp. 409–412, 2007. View at: Publisher Site  Google Scholar
 X. Chen, J. Ge, C. Wang, X. Zhang, and Q. Yu, “Noncircular DOA estimation algorithm via propagator method and Euler transformation,” in Proceedings of the 2017 3rd IEEE International Conference on Computer and Communications (ICCC), pp. 835–842, Chengdu, China, December 2017. View at: Publisher Site  Google Scholar
 W. Tan, X. Feng, W. Tan, D. Xie, L. Fan, and J. Li, “Direction finding for noncircular sources based on weighted unitary nuclear norm minimization,” in Proceedings of the 2017 IEEE 17th International Conference on Communication Technology (ICCT), pp. 1214–1218, Chengdu, China, October 2017. View at: Publisher Site  Google Scholar
 J. Xie, H. Tao, X. Rao, and J. Su, “Efficient method of passive localization for nearfield noncircular sources,” IEEE Antennas and Wireless Propagation Letters, vol. 14, pp. 1223–1226, 2015. View at: Publisher Site  Google Scholar
 J. Xie, J. Su, X. Rao, and H. Tao, “Realvalued localisation algorithm for nearfield noncircular sources,” Electronics Letters, vol. 51, no. 17, pp. 13301331, 2015. View at: Publisher Site  Google Scholar
 H. Chen, W. Wang, and W. Liu, “Joint DOA, range, and polarization estimation for rectilinear sources with a cold array,” IEEE Wireless Communications Letters, vol. 8, no. 5, pp. 1398–1401, 2019. View at: Publisher Site  Google Scholar
 B. A. Obeidat, Y. Zhang, and M. G. Amin, “Range and DOA estimation of polarized nearfield signals using fourthorder statistics,” in Proceedings of the 2004 IEEE International Conference on Acoustics, Speech, and Signal Processing, vol. 2, pp. 2–97, Montreal, Canada, May 2004. View at: Publisher Site  Google Scholar
 Y. Wu, H. Wang, Y. Zhang, and Y. Wang, “Multiple nearfield source localisation with uniform circular array,” Electronics Letters, vol. 49, no. 24, pp. 15091510, 2013. View at: Publisher Site  Google Scholar
 J.W. Tao, L. Liu, and Z.Y. Lin, “Joint DOA, range, and polarization estimation in the fresnel region,” IEEE Transactions on Aerospace and Electronic Systems, vol. 47, no. 4, pp. 2657–2672, 2011. View at: Publisher Site  Google Scholar
 J. He, M. O. Ahmad, and M. N. S. Swamy, “Nearfield localization of partially polarized sources with a crossdipole array,” IEEE Transactions on Aerospace and Electronic Systems, vol. 49, no. 2, pp. 857–870, 2013. View at: Publisher Site  Google Scholar
 Y. Xu, J. Ma, and Z. Liu, “Polarization sensitive PARAFAC beamforming for nearfield/farfield signals using cocentered orthogonal loop and dipole pairs,” in Proceedings of the 2013 IEEE China Summit and International Conference on Signal and Information Processing, pp. 616–620, July 2013. View at: Publisher Site  Google Scholar
 J. Volakis, Antenna Engineering Handbook, McGrawHill, New York, NY, USA, 4th edition, 2009.
 A. Ferréol, E. Boyer, and P. Larzabal, “Lowcost algorithm for some bearing estimation methods in presence of separable nuisance parameters,” Electronics Letters, vol. 40, no. 15, pp. 966967, 2004. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Jiaqi Song and Haihong Tao. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.