Direct Position Determination of Coherently Distributed Noncircular Sources

The existing localization algorithms for coherently distributed (CD) sources fail to achieve the best localization performance because of position information loss and error accumulation. To solve this problem, this study proposed a novel direct position determination (DPD) algorithm that profits from the characteristics of noncircular signals. Based on the parameterization assumption of CD sources, the DPD localization model is initially constructed and an extended subspace data fusion-based DPD algorithm is subsequently proposed by decomposing the extended covariance matrices, which are constructed by combining the characteristics of noncircular sources. The algorithm achieves a low complexity and a high efficiency by avoiding the calculation of the intermediate variables. Specifically, the closed-form expression of the Cramer-Rao lower bound for CD noncircular sources is also derived. Simulation results show that, compared with existing DPD algorithms, the proposed algorithm improves the multitarget localization capability and achieves high-accuracy results.


Introduction
Wireless target location technology is of considerable value in emergency rescue, safety management, and navigation planning.Numerous studies on this subject have indicated that this field has attracted the attention of scholars in the academe and industry [1][2][3].In most existing localization algorithms, the target is modeled as a point source, which has the advantages of simple calculation and high precision [4][5][6][7].However, in location applications, such as cellular mobile communication, low elevation radar, and underwater positioning, the target is usually crowded with surrounding scatterers that reflect the signal toward the receiving array.In such case, the target is not perceived as a point but rather as a spatial coherently distributed (CD) source, and the received signals are composed of multiple paths with angular expansion [8][9][10].If the point source model is still used for this condition, then the algorithm will rarely provide a high fidelity location result.The point source model is ideal, whereas the CD source model is closer to the actual scene.Therefore, analyzing the localization method for CD sources is more practically useful than that for point sources.
Obtaining the CD source position in a complex environment is difficult.Thus, the parameterization assumption combined with spatial distribution characteristics is often used to facilitate theoretical analysis.Jantti first proposed the concept of spatially distributed source and used the angle density function to reflect the angle expansion characteristics of received signals [11,12].Based on this concept, the Vec-MUSIC [13], DSPE [14], and ESPRIT [15] algorithms have been designed to enable parameter estimation of CD sources.As the special characteristics of noncircular signals, such as the common communication signals like binary phase shift keying (BPSK), amplitude modulation (AM), and preprocessed minimum shift keying (MASK), can improve estimation performance, the study of noncircular signals is of considerable practical value.Thus, Yang proposed a method for CD noncircular signals based on sparse representation, which exhibits a good performance under the condition of small snapshots with low SNR [16].Wan put forward a DSPE algorithm based on the combination of circular and noncircular signals in a large MIMO system, which improves the parameter estimation precision [17].To reduce computation cost, Yang combined the characteristics of noncircular signals 2 Wireless Communications and Mobile Computing with the ESPRIT algorithm to avoid complex parameter search and obtain high-accuracy results [18].Compared with the algorithms of circular signals, the algorithms based on the characteristics of noncircular signals enhance the parameter estimation accuracy of CD sources significantly.
After achieving parameter estimation of CD sources, target localization based on certain methods, such as Chan algorithm, Taylor series expansion, and Kalman filter, should be accomplished.The above-mentioned conventional localization methods employ two-step processing, whereby the measurement parameters (e.g., AOA and TOA) are initially extracted, and the source positions are subsequently estimated [19].The two-step method has several shortcomings.In the first step, the measurement parameters are acquired separately and independently from each receiver, disregarding the intrinsic constraint that the measurements are consistent with the same target position, thereby leading to loss of location information.In the second step, parameter estimation accuracy is greatly influenced by complex environment, and it is not easy to get high-accuracy estimation results, which inevitably involves processing errors that affect positioning accuracy.Therefore, the two-step method is suboptimal, failing to achieve the best performance [19].Compared with the two-step localization algorithm, direct position determination (DPD) algorithms proposed in recent years do not need to calculate the position step-by-step.Improved positioning performance can be obtained using the signal directly to estimate the position of the target [20][21][22].Developing the study of DPD algorithms is more valuable than that of the two-step algorithm.
The DPD algorithm has been extensively investigated in the past decade because of its excellent performance.Weiss first expounded the principle of DPD algorithms and proposed the maximum likelihood (ML) estimation-based method [23,24].To reduce the computational load generated by ML estimation, subspace data fusion-(SDF-) [25] and minimum variance distortionless response-based DPD algorithms [26] are proposed to decrease the dimension of the parameters to be estimated in order to mitigate the computation cost.As the DPD algorithm directly uses array outputs to localize targets, certain signal properties can be exploited to improve localization performance.Thus, the DPD algorithm based on the constant modulus, cyclostationary, and orthogonal frequency division multiplexing properties of the received signals achieve good localization performance [27][28][29].Bar proposed the CML-DPD algorithm for CD sources based on the generalized covariance matching method to improve the localization accuracy of CD sources [30].Although the algorithm outperforms conventional two-step methods in terms of the localization accuracy of CD sources, it is not practically useful because of its huge computational complexity.Moreover, localization performance still has room for improvement because signal characteristics are unused.
In summary, although the two-step localization method for CD sources improves the parameter estimation accuracy using the noncircular signal features, it cannot achieve optimal performance because of its limitations.In addition, the existing DPD algorithm designed for CD sources do not utilize signal characteristics and the computation cost is high, such that localization performance should be improved.Therefore, this study presents a novel DPD algorithm that profits from the characteristics of noncircular signals to achieve improved positioning accuracy of CD sources.An extended SDF method is proposed to reduce the computational load to improve the practical value of the algorithm.
The contributions of this study can be summarized as follows.
(1) The proposed DPD model for multiple CD sources fully employs the characteristics of noncircular signals, thereby extending the virtual array aperture of the receivers to improve the localization accuracy and number of detectable sources.To reduce computation cost, an extended SDF method is developed, which requires only low-dimensional optimization, resulting in low complexity and high efficiency.
(2) We derive the closed-form expression of the CRLB on the position estimation variance for noncircular signals to present a benchmark of the highest localization accuracy that the DPD method can reach.In the simulation part, the localization performance of the proposed algorithm can be proven to be close to the CRLB, which shows the effectiveness of our method.
This paper is organized as follows: Section 2 lists the notations used throughout the study.Section 3 describes the signal model for CD sources and formulates the problem.Section 4 derives the extended SDF-DPD estimator based on the characteristics of noncircular signals and analyzes the computation cost.Section 5 derives the closed-form expression of the CRLB on the position estimation variance.Section 6 presents several numerical simulations to validate the accuracy and reliability of the theoretical analysis.Finally, Section 7 concludes the study.

Notations
The notations used throughout this study are listed in Table 1.

System Model and Problem Formulation
In this study, the transmitter is assumed to be a far-field narrowband source, and we only consider the localization of CD sources based on the DOA information. receivers and  transmitters, which are all stationary, are employed.
All receivers equipped with a uniform linear antenna array are synchronized.The array is composed of  sensors with the first sensor located at u  = ( , ,  , ),  = 1, 2, . . ., , and the adjacent sensors are spaced at   .The transmitters are located at o  = ( , ,  , ),  = 1, 2, . . ., .The signals of different transmitters are independent of each other.At the th receiver, the received signal can be expressed as [14] where a  ( , (o  )) is the steering vector of the th transmitter,   ( , (o  ) −  , (o  ), ) is the angular density signal, and w  () is the statistically independent complex Gaussian noise. , (o  ) is the central direction of arrival as shown in For the CD sources, the angular density signal can be expressed as where   ( , (o  ) −  , (o  )) is the deterministic distributed function with the symmetric center of  , (o  ) satisfying the following condition: We let b , (o  ) be the modified steering vector which can be expressed as Then the observed signal can be sampled as where   [] is the th samples of   , r  [] is the th samples of r  , and  is the number of samples.We let T denotes the position vector of all transmitters.The received signal r  is then expressed as Based on the received signal r  according to (7), the problem addressed here is to estimate the source positions directly without explicitly computing the DOAs.The subsequent section will describe the novel DPD algorithm in detail.According to (6), when

DPD Algorithm of CD Noncircular Sources
) can be obtained, which is the same as that of the point source model.When   ( , (o  ) −  , (o  )) is a complex function, obtaining the analytical expression of b , (o  ) is difficult.Thus, the proposed algorithm is analyzed using the Gaussian angle of arrival (GAA) channel model as an example, similar to that in [30], to discuss the localization performance conveniently.The process employed in this study can be used as reference for the analysis of the performance of other distributed functions.
Based on the GAA channel model, the distribution of the angle perturbations generated by the scatterers is Gaussian with the mean of  , (o  ) and the standard deviation of  ,, , such that the th element of b , (o  ) can be expressed as where [b , (o  )]  denotes the th element of b , (o  ). , is given by In general, the distribution of the scatterers around the transmitter is assumed to be circular symmetric, and the variance of the distribution radius is set as   .Thus  ,, can be approximated as [30] where  , is the distance between the th transmitter and the th receiver.We then have the simplified form of (8) as follows [12]: where

Wireless Communications and Mobile Computing
Through the previously presented analysis, the analytical expression of b , (o  ) is obtained.The covariance matrix of r  is expressed in the following form: where  2 , is the noise power, I  is the  ×  identity matrix, and R  = E[xx H ] is the covariance matrix of signals.
The position information of the transmitters is included in R  .In [30], the CML-DPD method is proposed; however, the signal waveform should be known in advance and the computation cost is high when dealing with multitarget localization.Thus, we propose a novel SDF-DPD algorithm based on the characteristics of noncircular signals to achieve improved positioning accuracy of CD sources.
We briefly describe the statistical properties of noncircular signals.
x  is a noncircular signal with zero mean, then its variance and unconjugated variance satisfy the following condition: where 0 ≤  ≤ 1 is the noncircularity rate and   is the noncircularity phase.Given that the common signals, such as BPSK, MASK, and AM, are the strict-sense noncircular signals with the maximum noncircularity rate ( = 1), we only analyze the localization of strict-sense noncircular signals in this study.By substituting  = 1 into (15), the unconjugated covariance matrix of x can be given as where The diagonal elements are related to the noncircularity phases.The unconjugated covariance matrix of r  can then be expressed as We extend the observed signal as r = [r T  , r H  ] T to apply the unconjugated property of noncircular signals.Then the covariance matrix of r is structured as We let We set where g  = [ 1 ,  2 , . . .,   ] and c  = [ 1 ,  2 , . . .,   ].Thus (18) can be rewritten as From the previously presented analysis, the dimension of the virtual steering vector b NC , is twice as more than that of the original steering vector b , .By exploiting the characteristics of noncircular signals, the array aperture can be expanded, which will improve localization accuracy.Given that the virtual steering vectors of each source are independent, the subspace method can be used to reduce the dimension of the parameters to be estimated, which will significantly reduce the computation cost.The subspace decomposition of R can be given by where   is the eigenvalue and   is the eigenvector.Let U , be the signal subspace and U , be the noise subspace.Given that B NC, is orthogonal to U , , (b NC , ) H U , = 0 can be obtained.We can then derive the following cost function: The target position can be obtained by determining the minimum of (õ, c  , g  ).In reducing the number of estimated parameters further, we determine that g  is unrelated to the target position, such that it can be removed [14].Based on the characteristics of noncircular signals, the following expression can be obtained: where D ,1 and D ,2 have the same dimension.Then ( 22) can be rewritten as where Thus, the position parameters of the CD sources can be obtained using the following equation: In summary, compared with the CML-DPD algorithm proposed in [30], the proposed algorithm in this study reduces the dimension of the estimated parameters of multiple CD sources by the SDF method, which can significantly decrease computational load.Moreover, by applying the characteristics of noncircular signals, the signal subspace is enhanced and the noise subspace is extended, which will improve localization accuracy and the number of detectable CD sources.The simulation experiments will further verify the performance of the algorithm.Given that this algorithm combines the noncircular signal features and SDF method, it is called the NCSDF-DPD algorithm.To present the algorithm better, the NCSDF-DPD algorithm based on the principle derived in this section can be implemented as in Algorithm 1.

Computational Complexity Analysis.
Table 2 lists the main formulas and the amount of calculation in the NCSDF-DPD and CML-DPD algorithms [30] to illustrate the difference of the computation cost between the two algorithms.
The exhaustive search method is used in the two algorithms for reasonable comparison and analysis because of the large gap in the computational load of different search methods.In Table 2, the number of grids divided on each parameter is , and the parameter dimension that needs to be estimated for one source is .
The comparison results shown in Table 2 indicate that the CML-DPD algorithm requires a larger amount of calculation than the NCSDF-DPD algorithm because the dimension of the estimated parameters of the CML-DPD algorithm is larger than that of the NCSDF-DPD algorithm.However, the computational load of the NCSDF-DPD algorithm is significantly decreased, as expected.Moreover, through several fast calculation methods, such as genetic algorithm, particle filtering method, and expectation maximization algorithm, the computation cost will be further reduced to improve the localization performance.

CRLB for the Position of CD Noncircular Sources
The CRLB is a valid benchmark used to evaluate the best localization performance that a method can reach.This section derives the closed-form expression of the CRLB on the position estimation variance for CD noncircular sources.
The parameters that should be estimated can be expressed as where z  is the vector of signal power, and k  is the vector of noise power.The Fisher information matrix (FIM) of unknown vector y is given by where y  and y  denote the th and th elements of y, respectively, and J , signifies the (, )th element of J.According to (20),  R /o  is the partial derivative with respect to o  , which can be obtained using the following equation: We let [B NC, ] , ,  = 1, 2, . . ., ,  = 1, 2, . . .,  denote the (, )th element of B NC, , which can be expressed as ) . ( We then derive the following expression: where Wireless Communications and Mobile Computing 7 The  sources are independent of each other; then there is R  = diag{z  }.Thus we obtain the following expression: where diag{[0, 0, . . ., 1, . . ., 0, 0]} is the th element on the diagonal of R  /  which is 1, and the other elements are 0.
R / 2 , is the partial derivative with respect to  2 , , which can be obtained as follows: In conclusion, the partial derivative with respect to any element of y can be obtained, such that the FIM can be derived according to (29).We then determine the CRLB which is the inverse of J, as follows: The CRLB denotes the best performance that any unbiased estimator can achieve.When the parameter estimation model is determined, the CRLB can effectively verify the effectiveness of one algorithm by comparing it with the CRLB.In this study, the numerical simulation results of the CRLB will be obtained and compared with all algorithms to assess the localization performance of each algorithm.

Simulation Results
The simulation experiments will be performed in detail in this section to verify the localization performance of the proposed algorithm in Section 4. To analyze the performance differences between the existing algorithms, the experiments will include the following three algorithms: (1) The NCSDF-DPD algorithm proposed in this study (2) The DPD algorithm for CD source without using the characteristics of noncircular signals proposed in [30]: but as the ML method is too complex to get the simulation result for multisources, we replace it with SDF method for easy simulation, which has almost the same performance with the ML method (the method is denoted as CSDF-DPD) (3) The DPD algorithm for point source without using the characteristics of noncircular signals proposed in [25] (represented as SDF-DPD).
Moreover, the CRLBs for the NCSDF-DPD algorithm derived in Section 5 and for the CSDF-DPD algorithm obtained in [30] are presented to make our results as comprehensive as possible.
We assume that a 2D planar geometry has four observation stations located at The evaluation of the inverse cost functions, that is, 1/ f(õ, c  ), is conducted under the condition that SNR = 15 dB as shown in Figure 1.The 3D plots and their top views are presented to show the results clearly, wherein the and axes are the position coordinates.Figure 1(a) shows that the position estimations of two sources are seriously mixed up with a large fuzzy region because the SDF-DPD algorithm utilizes the point source model to locate CD sources.Thus obtaining the exact source number and the accurate position estimation is impossible.The results show that the point source model is invalid and cannot be used to solve the localization of CD sources.Figure 1(b) indicates that a sharp peak near the real position is obtained by the CSDF-DPD algorithm, which considers the distribution characteristics of the CD sources to improve localization performance.However, the algorithm cannot distinguish the two sources effectively because they are close to each other, resulting in only one location result with a large deviation.Figure 1(c) shows that, compared with the two other algorithms, the NCSDF-DPD algorithm can clearly distinguish two sources and yield high-accuracy location results by applying the characteristics of noncircular signals.The simulation results indicate that the distribution characteristics of CD sources should be considered in the localization algorithm; otherwise, the positioning accuracy will be poor.By contrast, the localization performance of the proposed algorithm is improved, as expected, which proves the effectiveness of the method that uses noncircular signal features to enhance the positioning accuracy of CD sources.
Given that the dimension of the virtual steering vector is two times more than that of the original steering vector, the array aperture is expanded by exploiting the characteristics of noncircular signals, which improves localization performance compared with that of the CSDF-DPD algorithm under the same array scale.We assumed that the three CD where   is the positioning error.The CDFs of the three algorithms are shown in Figure 3.The findings indicate that the positioning performance of the SDF-DPD algorithm worsens with the enlargement of the variance of the distribution radius, which again shows that the distribution of the source adversely affects the performance of conventional point source localization algorithms.The positioning performance of the CSDF-DPD algorithm is improved after considering the distribution characteristics of the sources but is still reduced with the enlargement of the variance of the distribution radius.Compared with the two other algorithms, the NCSDF-DPD algorithm is insensitive to the variance of the distribution radius after using the noncircular signal features to extend the subspace.Thus, with the increase in variance, the localization performance of the NCSDF-DPD algorithm does not decline and is better than that of the SDF-DPD and CSDF-DPD algorithms under any distribution radius, which further illustrates the effectiveness of the algorithm.The root-mean-square error (RMSE) and CRLB are regarded as functions of SNR across 100 Monte Carlo experiments, as shown in Figure 4, to illustrate the localization performance of each algorithm further.The RMSEs of the SDF-DPD and CSDF-DPD algorithms are large at low SNR, and their localization performances are nearly the same because noise is the main factor affecting positioning accuracy.However, at high SNR, the RMSE of the CSDF-DPD algorithm can almost reach the CRLB, whereas the SDF-DPD algorithm performs poorly.This performance may be attributed to the distribution characteristics of the sources at this time being the main factors that affect positioning accuracy.Compared with the CSDF-DPD algorithm, although the performance of the NCSDF-DPD algorithm is slightly reduced using the subspace method to decrease the computation cost, its CRLB and RMSE are lower, particularly when the SNR is low.The simulation results again demonstrate the effectiveness of the proposed algorithm because of the use of the characteristics of noncircular signals.
Owing to the expansion of the virtual steering vector, the number of samples required should be decreased.Under the condition that SNR = 5 dB, the influence of the number of samples on the positioning performance is shown in Figure 5.The results indicate that the performance of the SDF-DPD algorithm is not significantly improved with the increase in the number of samples, whereas the RMSE of the NCSDF-DPD algorithm decreases obviously, such that its performance is better than that of the SDF-DPD and CSDF-DPD algorithms.Moreover, the NCSDF-DPD algorithm under the condition of small samples can outperform the CSDF-DPD algorithm under the condition of large samples, which indicate that the observation time can be shortened to increase the localization timeliness.Table 2 lists the main formulas and the amount of calculation in the proposed method and the DPD method in [30].Now we give the running-times to further illustrate the computation cost of each method.The simulation experiment environment is as follows: Intel Core i7 CPU 2.4 GHz and MATLAB2015.The running-times are regarded as functions of the grids number divided on each parameter.Simulation results are shown in Figure 6 (for single source).We can conclude from the results that, with the increase of grids number, the running-time of the DPD method in [30] increases sharply.However, that of the proposed method in this study increases slowly.Moreover, the more the grids are, the more effective the proposed method is than the DPD method in [30].So the proposed method is more practical.

Conclusions
The localization of CD sources is a difficult problem in complex environment.The algorithm proposed in this study is closer to real applications in terms of localization scenarios and signal models.In order to further improve the positioning performance, this paper presents a DPD algorithm with combined the characteristics of noncircular sources and subspace data fusion method.The closed-form expression of CRLB on the position estimation variance for noncircular sources is also derived.In the simulation part, the localization performance of the proposed algorithm can be proven to be better than the existing algorithms, with higher positioning accuracy and lower computation cost.The proposed algorithm exhibits good performance when    The grids number divided on each parameter dealing with multitarget localization under the condition of small samples, which is of considerable value in practical applications.
sources are located at o 1 = (−1200, 0) m, o 2 = (500, 0) m, and o 3 = (1000, 500) m to verify this performance.The variance of distribution radius is c  = [200, 210, 220], and the noncircularity phases are g  = [/3, /5, /4] rad.The evaluation of the inverse cost functions of the CSDF-DPD and NCSDF-DPD algorithms is obtained under the condition that SNR = 15 dB, as shown in Figure 2. The results show that the CSDF-DPD algorithm can obtain only two position estimations of three CD sources, with low localization accuracy, whereas the NCSDF-DPD algorithm can derive all the position estimations of three CD sources, with high localization accuracy.The results can be attributed to the CSDF-DPD algorithm only able to locate two sources when the number of array sensors is three, whereas the

Figure 1 :
Figure 1: Evaluation of the inverse cost functions of the three algorithms in the scenario where four receivers with three sensors intercept two CD sources.(a) SDF-DPD algorithm, (b) CSDF-DPD algorithm, and (c) NCSDF-DPD algorithm.

Figure 2 :
Figure 2: Evaluation of the inverse cost functions of CSDF-DPD and NCSDF-DPD algorithms in the scenario wherein four receivers with three sensors intercept three CD sources.(a) CSDF-DPD algorithm and (b) NCSDF-DPD algorithm.

Figure 3 :
Figure 3: The CDF of three algorithms with different distribution radius: (a) source 1; (b) source 2.

Figure 4 :
Figure 4: Square root of CRLBs and estimated RMSEs of three algorithms as a function of SNR: (a) source 1 and (b) source 2.

Figure 5 :
Figure 5: Estimated RMSEs of three algorithms as a function of : (a) source 1 and (b) source 2.