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.
Fundamental Research Funds for the Central UniversitiesBDY06Innovation Project of Science of Technology Commission of the Central Military Commission∗∗-H863-∗∗-XJ-001-∗∗∗-021. 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 ⋅H, respectively. The symbol ⊗ represents the Kronecker product. ℜ⋅, ℑ⋅, and det⋅ 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 N=2M+1 dual-polarization sensors which is placed along the y-axis, and its sensors position is −Md,…,0,…,Md, 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.
Symmetric uniform linear array with dual-polarization sensors.
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 as(1)E=exey=−cosαcosθsinαejβ,where α 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 as(2)umx=−∑k=1Ksktejτmkcosαk+nmx,(3)umy=∑k=1Ksktejτmkcosθksinαkejβk+nmy,where skt denotes the kth signal, k=1,…,K, τmk represents the phase shift related to the kth signal’s propagation time delay from the reference point to the mth sensor, αk, βk, and θk are the auxiliary polarization angle, polarization phase difference, and DOA of kth signal, and nmx and nmy symbolize the additive noise.
Consider that the sources are located in the NF, the time delay τmk can be approximated as [4](4)τmk=mωk+m2ϕk,where ωk=−2πdsinθk/λ, ϕk=πd2cos2θk/λrk, λ denotes the wavelength of signal, and rk represents the range parameter of the kth source.
When the source signals are noncircular signals, it can be obtained that skt=ejψks0,kt; herein, ψk:k=1,…,K is the noncircular phase of the kth signals and s0,kt is the zero-phased version of the source signal.
In a matrix form, (2) and (3) can be written as(5)ut=Ast+nt,where(6)ut=uxtuytTuxt=u−Mx⋯u0x⋯uMx,uyt=u−My⋯u0y⋯uMy,A=a1a2⋯aK,st=s1ts2t⋯sKtT,nt=nxtnytT.
In the above equations, ut is array output vector, A represents the array steering matrix, st symbolizes the source signal vector, and nt denotes the additive noise matrix.
With a total of L snapshots taken at the distinct instants tl:l=1,…,L, the problem is to determine the localization parameters from the array output data. Throughout this paper, the following hypotheses are assumed to hold:
The incoming signals are mutually independent, narrowband stationary noncircular
The sensor noise is additive zero-mean white Gaussian and independent of the source signals
In order to avoid the phase ambiguity, the intersensor spacing d should be within a quarter wavelength
The range parameter r lies in 0.62D3/λ,2D2/λ 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 xt, respectively, by the following equations:(7)ℜxt=xt+x∗t2,ℑxt=xt−x∗t2j.
According to the signal model in (5), we can construct the real-valued augmented data matrix as(8)Yt=ℜuxtℑuxtℜuytℑuyt=Arst+nrt,where Ar is the augmented steering matrix and nr is the augmented noise matrix. And the th row of Ar can be expressed as(9)arθk,rk=ackx,askx,acky,askyT.
The elements in arθk,rk is formulated in (10), where m=−M,…,0,…,M:(10)ackx=−cosαkcosmωk+m2ϕk+ψk1×2M+1,acky=−cosαksinmωk+m2ϕk+ψk1×2M+1,ackx=cosθksinαkcosmωk+m2ϕk+βk+ψk1×2M+1,acky=cosθksinαksinmωk+m2ϕk+βk+ψk1×2M+1.
Then, the real-valued augmented covariance matrix of Yt can be expressed as R=EYtYHt. By taking the eigen-decomposition (EVD) of R, we have(11)R=UsΛsUsH+NsΛnNsH,where Λs and Λn are the diagonal matrixes which contain K largest eigenvalues and other 4N−K smallest eigenvalues. Us is a 4N×K matrix spanning the signal subspace of R, and Un is a 4N×4N−K matrix spanning the noise subspace of R,
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:(12)P=aHθ,rUnUnHaθ,r.
It is obvious that K sets of parameters can be achieved by finding the K minimum value of Pθ,r. 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:(13)aθ,r=VθFθ,rg,where(14)Vθ=I2×2⊗V1−V2V2V1,Fθ,r=I2×2⊗f1−f2f2f1.Herein,(15)V1=cos−Mωkcos−M−1ωk⋯1⋯cosM−1ωkcosMωk2M+1×M+1,(16)V2=sin−Mωksin−M−1ωk⋯0⋯sinM−1ωksinMωk2M+1×M+1,(17)f1=cosM2ϕkcosM−12ϕk⋯1M+1×1T,(18)f2=sinM2ϕksinM−12ϕk⋯0M+1×1T,(19)g=−cosαkcosψk−cosαksinψksinαkcosθkcosψk+βksinαkcosθksinψk+βk.
If and only if θ,r=θk,rk,k=1,…,K, the following equation holds:(20)aHθ,rUnUnHaθ,r=0.
Equation (20) can be rewritten as(21)QHθ,rC1θQθ,r=0.where Qθ,r=Fθ,rg and C1θ=VHθUnUnHVθ.
It is noteworthy that C1 is only related to θ and independent with other parameters. Since C1θ=VHθUnUnHVθ and Vθ∈R42M+1×4M+1, Un∈R42M+1×42M+1−K, so if 42M+1−K≥4M+1i.e., K≤4M, C1θ is generally of full column rank. Thus, only if θ=θk,k=1,…,K, (21) becomes true because the rank of C1θ drops according to the RARE criterion. Therefore, the DOAs of sources can be estimated by the following spectrum function:(22)f1θ=1detC1θ.
After finding K peaks of f1θ 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 C2θ,r=FHθ,rVHθUnUnHVθFθ,r, then (20) can be expressed as(23)gHC2θ,rg=0.
Similarly, since g≠0, (23) can hold true if C2θ drops rank. Therefore, with obtained θ^k, the range of sources can be estimated by of the following 1-D spectrum searching function for K times:(24)f2r=1detC2θ^,r.
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(25)R^Y=1L∑t=1LYtYTt.
To summarize, the procedures of the proposed method are shown as follows:
Reconstruct the real-valued augmented data matrix Yt based on array output matrix and Euler equation.
Take eigen-decomposition operation of the covariance matrix RY, and obtain noise subspace matrix Un.
Compute rank reduction matrix C1, and estimate DOAs of signals by searching the K highest peaks of (22).
Compute rank reduction matrix C2, with obtained DOAs θ^k,k=1,…,K estimate the range parameter r^k by searching the highest peak of (24), repeat this step from k=1 to K. If the polarization parameters of signals need to be estimated, do step 5.
Decouple the vector g like (13) further and separate the noncircular parameter from polarization; the polarization can be obtained by another RARE estimator.
Insert estimated θ^k, r^k, α^k, and β^k into (12). By searching the highest peaks of 1/P, the noncircular parameter is estimated. Repeat this step from k=1 to K.
4. Discussion4.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 N=2M+1 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 2M sources. For the proposed method, the noncircularity has been utilized to construct an extended subspace; hence, it can estimate 4M 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 θ∈−90∘,90∘ and range parameter 0.62D3/λ,2D2/λ are denoted as Δθ and Δr. 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 N×N second-order covariance matrices, two EVDs, and twice spectral searches for DOA and range estimation. FR-RARE builds a 2N×2N 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 4N×4N second-order real-valued covariance matrix, perform EVD on this matrix, and execute 1-D spectral search (Table 1).
Computational complexity of different algorithms.
Method
Matrices
EVDs
Spectral search
GESPRIT
O2N2/L
8N3/3
KN2/ΔrRmax−Rmin+360N2/Δθ
FR-RARE
O4N2/L
32N3/3
4KN2/ΔrRmax−Rmin+720N2/Δθ
Proposed
O4N2/L
32N3/3
16KN2/ΔrRmax−Rmin+2880N2/Δθ
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 M=2 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 as(26)RMSE=1500∑n=1500y^k−yk2,where yk is the exact DOA θk or the range rk, and y^k denotes the estimation of yk.
In the first experiment, we suppose that five BPSK equi-power signals are located at −20∘,0.2λ,−20∘,−20∘, −10∘,0.3λ,−10∘,−10∘, 0∘,0.4λ,0∘,0∘, 10∘,0.5λ,10∘,10∘, and 20∘,0.6λ,20∘,20∘, respectively. The snapshot number and signal-to-noise ratio (SNR) are set as 500 and 20 dB. 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 4M sources while FR-RARE can only handle 2M sources. Since GESPRIT brings half aperture loss, it can handle M sources at most.
DOA spectrum of the proposed method.
DOA spectrum of the proposed method.
DOA spectrum of the proposed method, FR-RARE and GESPRIT.
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 −5∘,1.2λ,−5∘,−5∘ and 2∘,1.8λ,2∘,2∘. The snapshots number is fixed at 500 and the SNR varies from 0 dB to 40 dB in a stepsize of 5 dB. 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.
RMSEs of DOA estimates versus SNR.
RMSEs of range estimates versus SNR.
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 20 dB, 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.
RMSEs of DOA estimates versus snapshots.
RMSEs of range estimates versus snapshots.
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.
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-∗∗-XJ-001-∗∗-02).
KrimH.VibergM.Two decades of array signal processing research: the parametric approach1996134679410.1109/79.5268992-s2.0-0030193445SchmidtR.Multiple emitter location and signal parameter estimation198634327628010.1109/tap.1986.1143830RoyR.KailathT.ESPRIT-estimation of signal parameters via rotational invariance techniques198937798499510.1109/29.322762-s2.0-0024702515HuangY.-D.BarkatM.Near-field multiple source localization by passive sensor array199139796897510.1109/8.869172-s2.0-0026186340ChallaR.-N.ShamsunderS.High-order subspace-based algorithms for passive localization of near-field sources2Proceedings of the Conference Record of the Twenty-Ninth Asilomar Conference on Signals, Systems and ComputersOctober 1995Pacific Grove, CA, USA77778110.1109/acssc.1995.5408062-s2.0-0003736502LeeJ.-H.ChenY.YehC.-C.A covariance approximation method for near-field direction-finding using a uniform linear array19954351293129810.1109/78.3824212-s2.0-0029309388NohH.LeeC.A covariance approximation method for near-field coherent sources localization using uniform linear array201540118719510.1109/joe.2013.22498722-s2.0-84920947818GrosickiE.Abed-MeraimK.Yingbo HuaY.A weighted linear prediction method for near-field source localization200553103651366010.1109/tsp.2005.8551002-s2.0-27844457979ZhiW.ChiaM. Y.-W.Near-field source localization via symmetric subarrays200714640941210.1109/lsp.2006.8883902-s2.0-34249803220ChenX.GeJ.WangC.ZhangX.YuQ.Noncircular DOA estimation algorithm via propagator method and Euler transformationProceedings of the 2017 3rd IEEE International Conference on Computer and Communications (ICCC)December 2017Chengdu, China83584210.1109/compcomm.2017.83226612-s2.0-85049742494TanW.FengX.TanW.XieD.FanL.LiJ.Direction finding for non-circular sources based on weighted unitary nuclear norm minimizationProceedings of the 2017 IEEE 17th International Conference on Communication Technology (ICCT)October 2017Chengdu, China1214121810.1109/icct.2017.83598282-s2.0-85047729672XieJ.TaoH.RaoX.SuJ.Efficient method of passive localization for near-field noncircular sources2015141223122610.1109/lawp.2015.23995342-s2.0-84924769433XieJ.SuJ.RaoX.TaoH.Real-valued localisation algorithm for near-field non-circular sources201551171330133110.1049/el.2015.04542-s2.0-84939142431ChenH.WangW.LiuW.Joint DOA, range, and polarization estimation for rectilinear sources with a cold array2019851398140110.1109/lwc.2019.29195422-s2.0-85067984638ObeidatB. A.ZhangY.AminM. G.Range and DOA estimation of polarized near-field signals using fourth-order statistics2Proceedings of the 2004 IEEE International Conference on Acoustics, Speech, and Signal ProcessingMay 2004Montreal, Canada29710.1109/icassp.2004.1326203WuY.WangH.ZhangY.WangY.Multiple near-field source localisation with uniform circular array201349241509151010.1049/el.2013.20122-s2.0-84890035483TaoJ.-W.LiuL.LinZ.-Y.Joint DOA, range, and polarization estimation in the fresnel region20114742657267210.1109/taes.2011.60346572-s2.0-80054002202HeJ.AhmadM. O.SwamyM. N. S.Near-field localization of partially polarized sources with a cross-dipole array201349285787010.1109/taes.2013.64943852-s2.0-84876141843XuY.MaJ.LiuZ.Polarization sensitive PARAFAC beamforming for near-field/far-field signals using co-centered orthogonal loop and dipole pairsProceedings of the 2013 IEEE China Summit and International Conference on Signal and Information ProcessingJuly 201361662010.1109/chinasip.2013.66254152-s2.0-84889589995VolakisJ.20094thNew York, NY, USAMcGraw-HillFerréolA.BoyerE.LarzabalP.Low-cost algorithm for some bearing estimation methods in presence of separable nuisance parameters2004401596696710.1049/el:200405372-s2.0-3843135128