Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2019 / Article

Research Article | Open Access

Volume 2019 |Article ID 5759427 | 8 pages |

Efficient Localization Algorithm for Near-Field Noncircular Sources via Dual-Polarization Sensor Array

Academic Editor: Fazal M. Mahomed
Received29 Jul 2019
Accepted27 Nov 2019
Published12 Dec 2019


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 real-valued covariance matrix is constructed to transform complex-valued computation into real-valued computation. Based on the property of noncircular signals and symmetric uniform linear array (SULA) which consist of dual-polarization 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 real-valued 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 far-field (FF), which means the wavefronts are plane waves and therefore only direction-of-arrival (DOA) parameters are required to be estimated.

However, when the radiating sources are located in the near-field (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 two-dimensional (2-D) MUSIC method in the angle-range domain to achieve NF super-resolution localization, but the 2-D joint search brought high computational complexity [4]. Challa and Shamsunder took the lead in introducing high-order cumulant into NF source parameter estimation problem. By constructing multiple cumulant matrixes, they proposed ESPRIT-Like to estimate DOA and range parameters of NF sources [5]. However, the ESPRIT-Like 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 direction-finding 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 real-valued localization algorithm for noncircular signals using the uniform linear array [12]. Furthermore, Xie et al. proposed another near-field 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 ESPRIT-Like 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 squares-virtual ESPRIT algorithm (LS-VESPA) [16]. But it involves extra parameter pairing procedure. Based on the symmetric sparse linear array with dual-polarization sensors, Tao et al. proposed the Fresnel-region rank reduction (FR-RARA) algorithm [17] that enhanced array aperture and only required second-order statistics. He et al. presented a NF localization of partially polarization sources with a cross-dipole 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 real-valued 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 one-dimensional (1-D) spectral searches. The proposed method has the following advantages: (1) As a result of the exploitation of non-circularity, more sources can be resolved. (2) It is efficient since it avoids exhaustive complex-valued 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 dual-polarization 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 dual-polarization sensors which is placed along the y-axis, and its sensors position is , where d is the interelement spacing. The dual-polarization sensors used in this paper is cross-dipole. This localization algorithm is achievable by other polarization sensors, such as cocentered orthogonal loop and dipole pair [19]. But the cross-dipole is very small and easier to design, so it is more common in practice.

Note that the two polarization components of the cross-dipole point to the x-axis and y-axis directions, respectively. Assuming that all sources are located in the y-z 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 x-polarization and y-polarization 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 zero-phased 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 zero-mean 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 real-valued 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 1-D search, the localization parameters can be estimated. Consequently, the multidimensional optimization problem could be accomplished by real-valued computation and 1-D spectral searches.

3.1. Real-Valued Augmented Covariance Matrix

In order to transform the complex-valued data into real-valued 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 real-valued 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 real-valued augmented covariance matrix of can be expressed as . By taking the eigen-decomposition (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 2-D spectral search. To avoid the high-dimensional computation, the augmented steering vector can be factorized, leading to a simple 1-D 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 1-D 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 1-D 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 real-valued augmented data matrix based on array output matrix and Euler equation.(2)Take eigen-decomposition 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], FR-RARE [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 subspace-based algorithm needs at least one eigenvector to span noise space, FR-RARE 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 FR-RARE.

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 second-order covariance matrices, two EVDs, and twice spectral searches for DOA and range estimation. FR-RARE builds a second-order covariance matrices, performs the EVD, and twice spectral searches for DOA and range estimation. For the proposed one, it needs to establish a second-order real-valued covariance matrix, perform EVD on this matrix, and execute 1-D spectral search (Table 1).

MethodMatricesEVDsSpectral search


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 dual-polarization sensors with the interelement spacing being a quarter-wavelength. The NF signal sources impinged upon the array are equipowered, statistically independent BPSK signal (code rate = 0.1/Ts, non-circularity = 1). Moreover, the estimation performance is measured by the root mean-square 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 equi-power signals are located at , , , , and , respectively. The snapshot number and signal-to-noise 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 FR-RARE [17] are plotted in Figure 4. By employing real-valued computations, the proposed method can resolve up to sources while FR-RARE 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 FR-RARE 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 FR-RARE. 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 near-field noncircular signals is presented by employing the real-valued computation and 1-D 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 near-field.

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.


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-∗∗-XJ-001-∗∗-02).


  1. 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
  2. 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
  3. R. Roy and T. Kailath, “ESPRIT-estimation 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
  4. Y.-D. Huang and M. Barkat, “Near-field 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
  5. R.-N. Challa and S. Shamsunder, “High-order subspace-based algorithms for passive localization of near-field sources,” in Proceedings of the Conference Record of the Twenty-Ninth Asilomar Conference on Signals, Systems and Computers, vol. 2, pp. 777–781, Pacific Grove, CA, USA, October 1995. View at: Publisher Site | Google Scholar
  6. J.-H. Lee, Y. Chen, and C.-C. Yeh, “A covariance approximation method for near-field direction-finding using a uniform linear array,” IEEE Transactions on Signal Processing, vol. 43, no. 5, pp. 1293–1298, 1995. View at: Publisher Site | Google Scholar
  7. H. Noh and C. Lee, “A covariance approximation method for near-field 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
  8. E. Grosicki, K. Abed-Meraim, and Y. Yingbo Hua, “A weighted linear prediction method for near-field source localization,” IEEE Transactions on Signal Processing, vol. 53, no. 10, pp. 3651–3660, 2005. View at: Publisher Site | Google Scholar
  9. W. Zhi and M. Y.-W. Chia, “Near-field source localization via symmetric subarrays,” IEEE Signal Processing Letters, vol. 14, no. 6, pp. 409–412, 2007. View at: Publisher Site | Google Scholar
  10. 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
  11. W. Tan, X. Feng, W. Tan, D. Xie, L. Fan, and J. Li, “Direction finding for non-circular 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
  12. J. Xie, H. Tao, X. Rao, and J. Su, “Efficient method of passive localization for near-field noncircular sources,” IEEE Antennas and Wireless Propagation Letters, vol. 14, pp. 1223–1226, 2015. View at: Publisher Site | Google Scholar
  13. J. Xie, J. Su, X. Rao, and H. Tao, “Real-valued localisation algorithm for near-field non-circular sources,” Electronics Letters, vol. 51, no. 17, pp. 1330-1331, 2015. View at: Publisher Site | Google Scholar
  14. 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
  15. B. A. Obeidat, Y. Zhang, and M. G. Amin, “Range and DOA estimation of polarized near-field signals using fourth-order 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
  16. Y. Wu, H. Wang, Y. Zhang, and Y. Wang, “Multiple near-field source localisation with uniform circular array,” Electronics Letters, vol. 49, no. 24, pp. 1509-1510, 2013. View at: Publisher Site | Google Scholar
  17. 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
  18. J. He, M. O. Ahmad, and M. N. S. Swamy, “Near-field localization of partially polarized sources with a cross-dipole array,” IEEE Transactions on Aerospace and Electronic Systems, vol. 49, no. 2, pp. 857–870, 2013. View at: Publisher Site | Google Scholar
  19. Y. Xu, J. Ma, and Z. Liu, “Polarization sensitive PARAFAC beamforming for near-field/far-field signals using co-centered 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
  20. J. Volakis, Antenna Engineering Handbook, McGraw-Hill, New York, NY, USA, 4th edition, 2009.
  21. A. Ferréol, E. Boyer, and P. Larzabal, “Low-cost algorithm for some bearing estimation methods in presence of separable nuisance parameters,” Electronics Letters, vol. 40, no. 15, pp. 966-967, 2004. View at: Publisher Site | Google Scholar

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.

195 Views | 59 Downloads | 0 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder
 Sign up for content alertsSign up