International Journal of Antennas and Propagation

Volume 2013 (2013), Article ID 657653, 7 pages

http://dx.doi.org/10.1155/2013/657653

Research Article

Passive Localization of 3D Near-Field Cyclostationary Sources Using Parallel Factor Analysis

Jian Chen, Guohong Liu, and Xiaoying Sun

College of Communication and Engineering, Jilin University, Changchun 130022, China

Received 6 February 2013; Accepted 2 June 2013

Academic Editor: Krzysztof Kulpa

Copyright © 2013 Jian Chen et al. 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.

Abstract.

By exploiting favorable characteristics of a uniform cross-array, a passive localization algorithm of narrowband cyclostationary sources in the spherical coordinates (azimuth, elevation, and range) is proposed. Firstly, we construct a parallel factor (PARAFAC) analysis model by computing the third-order cyclic moment matrices of the properly chosen sensor outputs. Then, we analyze the uniqueness of the constructed model and obtain three-dimensional (3D) near-field parameters via trilinear alternating least squares regression (TALS). The investigated algorithm is well suitable for the localization of the near-field cyclostationary sources. In addition, it avoids the multidimensional search and pairing parameters. Results of computer simulations are carried out to confirm the satisfactory performance of the proposed method.

1. Introduction

There has been considerable interest in bearing estimation for radar, sonar, communication, and electronic surveillance [1]. Various high-resolution algorithms, such as MUSIC [2] and ESPRIT [3], have been proposed to obtain the direction-of-arrival estimation of the far-field sources. In addition, when the sources are localized at the Fresnel region [4] of the array aperture, both the azimuth and range should be estimated. Recently, a significant amount attention has been paid to this issue and several near-field sources localization algorithms [5–7] are also available. However, all these methods as mentioned above only address the 2D problem of estimating azimuth and range and rely on the assumption of the stationary sources.

In recent years, several 3D near-field sources localization methods have been developed to obtain azimuth, elevation, and range. Meraim and Hua [8] proposed a second-order based method to cope with this issue. By translating the 1D uniform linear array of near field into a virtual rectangular array of virtual far field, Challa and Shamsunder [9] developed a fourth-order cumulants based Unitary-ESPRIT method. Moreover, [10, 11] also introduced efficient localization methods for 3-D near-field non-Gaussian stationary sources. It is obviously seen that all these methods still require that the incoming signals should be stationary ones; in addition, a multidimensional search or pairing parameters is also failed to avoid.

Cyclostationarity, which is a statistical property processed by most man-made communication signals, is related to the underlying periodicity arising from cyclic frequency of based rates. Besson et al. [12] introduced an original far-field approximation [13] based method to deal with the bearings estimation of near-field cyclostationary sources. Due to the fact that the effect caused by the mismatch between the actual spherical wavefront phase vector and assumed planar wavefront vector was alleviated, the proposed technique showed a satisfactory performance for the near-field sources far away from the sensor array. The estimations for elevation and range, however, have not been well considered.

In this paper, we consider the problem of jointly estimating elevation, azimuth, and range of the near-field cyclostationary sources; what is more, a two-stage passive localization method has been proposed. In the first stage, several third-order cyclic moment matrices of a cross-array observations data are computed, and a parallel factor analysis model in the cyclic statistic domain is constructed. In the second stage, the uniqueness of the constructed model is proved; in addition, the 3-D localizing parameters for the near-field sources have also been obtained through TALS. The algorithm developed in this paper would be well suitable for near-field cyclostationary sources, and it does not require multidimensional search or pairing parameters; in addition, it also can effectively alleviates the array aperture loss.

The rest of this paper is organized as follows. Section 2 introduces the signal model of near-field localization based on cross-array. Section 3 develops a joint estimation algorithm of three parameters in near-field. Section 4 shows simulation results. Section 5 presents the conclusion of the whole paper.

2. Near-Field Signal Model Based Cross-Array

2.1. Near-Field Signal Model

We consider a near-field scenario of uncorrelated narrowband signals impinging on a cross-array signed with the and axes (Figure 1), which consists of sensors with element spacing .

Figure 1: Sensor-source configuration for the near-field problem.

Let the array center be the phase reference point, and the signals received by the th and the th can be, respectively, expressed as

where donates the sensor additive noise, and

where , , and indicate elevation, azimuth, and range of th signal, respectively, and is wavelength of source signal.

The th source signal with amplitude can be modeled as

where is the center frequency.

2.2. Assumption of Signal Model

The main problem addressed in this paper is to jointly estimate the sets of parameters , and then the following assumptions are assumed to hold:

(1)the envelope is non-Gaussian stationary random process with zero mean and nonzero skewness;(2)the sensor noise is the additive stationary one and independent from the source signals;(3)the sensor array is a uniform linear array with element spacing ; in addition, the source number is not more than sensor number .

2.3. Parallel Factor Analysis

We need to introduce the following notation that will be used in the sequel.

Definition 1 (see [14]). Let stand for the () element of a three-dimensional tensor , if

where denotes the element of matrix and similarly for the others. Equation (4) indicates as a sum of triple products, which is variably known as the trilinear model, trilinear decomposition, triple product decomposition, canonical decomposition, or parallel factor (PARAFAC) analysis.

Definition 2 (see [14]). For a matrix , if all columns of are linearly independent, but there exists a collection of linearly dependent columns of , then it has Kruskal-rank (-rank) .

Theorem 3 (see [15]). Consider a three dimensional tensor as defined in (4), and represents the common dimension, if

then , , and are unique up to permutation and (complex) scaling of columns.

3. PARAFAC Based 3D Near-Field Sources Localization

3.1. Compute the Third-Order Cyclic Moments Matrices

Let denote the third-order cyclic moment with cyclic frequency defined as

where the superscript “*” denotes the complex conjugate operation; in addition, can be given by

where , , and denote mean, variance, and the third-order moment of ; in addition, , , and represent mean, variance, and the third-order cyclic moment of , respectively.

With assumption being considered, we further obtain

where is impulse function. Letting and substituting (8) into (6) yield

Based on (9), we reconstruct the spatial third-order cyclic moment matrix , in which the th element can be expressed as follows:

In a matrix form, (10) can be written as

where the superscript “” represents conjugate transpose operation, denotes the third-order moment matrix of source signals, and

On the other hand, following the same process described above, we can easily obtain

And the elements of the other three matrices can be expressed as

Finally, we have

with

3.2. Build the Parallel Factor Analysis Model

Considering the situation of limited samples, we build a parallel factor analysis model that uses the third-order cyclic moments as

with , and the Khatri-Rao product [15] for (17) shows

where

with denoting a row vector consisting of diagonal matrix .

Similarly, (18) also yields

3.3. Solve the Parallel Factor Analysis Model, and Estimate 3-D Parameters

As it stands, and are both Vandermonde matrices, and then they have Kruskal-rank (-rank) . On the other hand, the -rank of will be . When the condition that the number of signals being holds, then , , and are unique up to permutation and scaling of columns. With trilinear alternating least squares regression, we obtain that

Then using these estimates, we can associate with each pair as follows:

Finally, the sources parameters can be estimated as

4. Computer Simulation Results

In this section, we explicit some simulation results to evaluate the performance of proposed method. For all examples, a uniform linear array with number of 15 sensors and element spacing 0.25 is displayed, where is the wavelength of the narrowband source signals. Two near-field or far-field equal power cyclostationary signals are impinging on the array with center frequency 0.25, and their envelopes can be modeled as exponentially distributed. To be compared, we simultaneously execute the method of [12], which is suitable for the 3-D near-field non-Gaussian stationary sources localization. The presented results are evaluated by the estimated root mean square error (RMSE) from the averaged results of 200 independent Monte-Carlo simulations.

In the first example, Two near-field sources are located at (35°, 46°, and 1/6) and (20°, 60°, and 2/5), respectively. The additive noise is time and spatially white Gaussian with zero mean and unit variance. In addition, the snapshot number is set equal to 1024. When SNR varies from 0 dB to 25 dB, the RMSE of the elevation, azimuth, and range estimations for two near-field cyclostationary sources can be shown in Figures 2, 3, and 4. For the comparison, the performance in the same situation for the fourth-order cumulants based method has also been displayed. From these three figures, it is obvious that the proposed method performs better in elevation, azimuth, and range estimation than the FFA based algorithm for all SNRs. What is more, the RMSE of range estimations for the first source that is closer to the array is less than the second one. This phenomenon is well agreement with the theoretical analysis that the sources closer to sensor array would hold a smaller standard deviation than the one far from the array [7].

Figure 2: The RMSE of elevation estimations for two near-field sources using the proposed method and the cumulant based method versus SNR.

Figure 3: RMSE of azimuth estimations for two near-field sources using the proposed method and the cumulant based method versus SNR.

Figure 4: The RMSE of range estimations for two near-field sources using the proposed method and the cumulant based method versus SNR.

In the second example, the simulation condition is similar to the first example, except that the SNR is set at 10 dB, and the snapshot number is varied from 400 to 2000. The RMSE of elevation, azimuth, and range estimations of two near-field cyclostationary sources obtained from the proposed method as well as the cumulant based method can be displayed in Figures 5, 6, and 7. From these three figures, we can see that the proposed method still shows a more satisfactory accuracy than the cumulant based method in all available snapshots. In addition, the range estimations for the first sources are much better than the second one, while the elevation and azimuth estimations for both near-field sources are similar for the proposed method.

Figure 5: The RMSE of elevation estimations for two near-field sources using the proposed method and the cumulant based method versus snapshot number.

Figure 6: The RMSE of azimuth estimations for two near-field sources using the proposed method and the cumulant based method versus snapshot number.

Figure 7: The RMSE of range estimations for two near-field sources using the proposed method and the cumulant based method versus snapshot number.

In the third example, the additive noise is time and spatially white Gaussian random process and the colored noise from a second-order AR model with parameter (−1.8,0.9), respectively. In addition, the other simulation condition is similar to the first example, except that the SNR is set at 10 dB, and snapshot number is equal to 1024. The mean and variance of the elevation, azimuth, and range estimations for two near-field cyclostationary sources can be shown in Table 1. From this table, we can easily see that the effectiveness of the proposed method is not little affected by the difference of the sensor noise.

Table 1: RMSE of 3D parameters for near-field localization in the third example.

Example
Source
True
Mean
Variance
Source
True
Mean
Variance

White noise

35°
34.95°
0.1923

20°
20.03°
0.1647

40°
39.99°
0.0419

60°
59.99°
0.1000

λ/6
0.1669λ
3.5744 × 10⁻⁵

2λ/5
0.3984λ
8.2107 × 10⁻⁴

Colored noise

35°
35.02°
0.2146

20°
19.98°
0.0618

40°
40.03°
0.0795

60°
60.02°
0.1165

λ/6
0.1662λ
2.8455 × 10⁻⁵

2λ/5
0.4036λ
8.4173 × 10⁻⁴

5. Conclusion

This paper considers the problem of the passive localization of 3-D near-field cyclostationary sources and proposes an efficient third-order cyclic moment based algorithm. The construction of the parallel factor analysis model in the cyclic domain avoids the multidimensional search and pairing parameters. Moreover, the utilization of trilinear alternating least squares regression deals with the joint estimation of elevation, azimuth, and range. From the simulation results mentioned above, the proposed method outperforms the cumulant based method in locating near-field cyclostationary sources; in addition, the stationary noise has a little influence on the estimated accuracy of the proposed method.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (61171137) and Specialized Research Fund for the Doctoral Program of Higher Education (20090061120042).

References

J. Liang, X. Zeng, B. Ji, J. Zhang, and F. Zhao, “A computationally efficient algorithm for joint range-DOA-frequency estimation of near-field sources,” Digital Signal Processing, vol. 19, no. 4, pp. 596–611, 2009.

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.

A. J. Weiss and B. Friedlander, “Range and bearing estimation using polynomial rooting,” IEEE Journal of Oceanic Engineering, vol. 18, no. 2, pp. 130–137, 1993.

W. Zhi and M. Y.-W. Chia, “Near-field source localization via symmetric subarrays,” in Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '07), pp. II1121–II1124, Honolulu, Hawaii, USA, April 2007.

D. Starer and A. Nehorai, “Path-following algorithm for passive localization of near-field sources,” in Proceedings of the 5th ASSP Workshop on Spectrum Estimation and Modeling, pp. 322–326, 1990.

D. Starer and A. Nehorai, “Passive localization on near-field sources by path following,” IEEE Transactions on Signal Processing, vol. 42, no. 3, pp. 677–680, 1994.

R. A. Harshman, “Foundation of the PARAFAC procedure: model and conditions for an explanatory mutil-mode factor analysis,” UCLA Working Papers in Phonetics, vol. 22, no. 7, pp. 111–117, 1972.

K. Aberd Meraim and Y. Hua, “3-D near field source localization using second order statistics,” in Proceedings of 31st Alisomar Conference on Signals, Systems and Computers, vol. 2, pp. 1307–1310, 1997.

R. N. Challa and S. Shamsunder, “Passive near-field localization of multiple non-Gaussian sources in 3-D using cumulants,” Signal Processing, vol. 65, no. 1, pp. 39–53, 1998.

K. Deng and Q. Yin, “Closed form parameters estimation for 3-D near field sources,” in Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '06), pp. IV1133–IV1136, Toulouse, France, May 2006.

C. M. Lee, K. S. Yoo, and K. K. Lee, “Efficient algorithm for localizing 3-D narrowband multiple sources,” IEE Proceedings of Radar, Sonar and Navigation, vol. 38, no. 1, pp. 23–26, 2001.

O. Besson, P. Stoica, and A. B. Gershman, “Simple and accurate direction of arrival estimator in the case of imperfect spatial coherence,” IEEE Transactions on Signal Processing, vol. 49, no. 4, pp. 730–737, 2001.

A. L. Swindlehurst and T. Kailath, “A performance analysis of subspace-based methods in the presence of model errors—I: the MUSIC algorithm,” IEEE Transactions on Signal Processing, vol. 40, no. 7, pp. 1758–1774, 1992.

N. D. Sidiropoulos, R. Bro, and G. B. Giannakis, “Parallel factor analysis in sensor array processing,” IEEE Transactions on Signal Processing, vol. 48, no. 8, pp. 2377–2388, 2000.

J.-L. Liang, B.-J. Ji, F. Zhao, and J.-Y. Zhang, “Near-field source localization algorithm using parallel factor analysis,” Acta Electronica Sinica, vol. 35, no. 10, pp. 1909–1915, 2007.