International Journal of Antennas and Propagation

Volume 2014 (2014), Article ID 849039, 7 pages

http://dx.doi.org/10.1155/2014/849039

Research Article

Joint 2D Direction-of-Arrival and Range Estimation for Nonstationary Sources

Jian Chen, Hui Zhao, Xiaoying Sun, and Guohong Liu

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

Received 10 March 2014; Accepted 1 July 2014; Published 17 July 2014

Academic Editor: Michelangelo Villano

Copyright © 2014 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.

Passive localization of nonstationary sources in the spherical coordinates (azimuth, elevation, and range) is considered, and a parallel factor analysis based method is addressed for the near-field parameter estimation problem. In this scheme, a parallel factor analysis model is firstly constructed by computing five time-frequency distribution matrices of the properly chosen observation data. In addition, the uniqueness of the constructed model is proved, and both the two-dimensional (2D) direction-of-arrival (DOA) and range can be jointly obtained via trilinear alternating least squares regression (TALS). The investigated algorithm is well suitable for near-field nonstationary source localization and does not require parameter-pairing or multidimensional search. Several simulation examples confirm the effectiveness of the proposed algorithm.

1. Introduction

Bearing estimation has been a strong interest in radar and sonar as well as communication. In the last three decades, various high-resolution algorithms for direction finding of multiple narrowband sources assume that the propagating waves are considered to be plane waves at the sensor array. However, when the sources are located in the Fresnel region [1] of the array aperture, the wavefronts emitted from these sources are spherical rather than planar at each sensor position and characterized by both the DOA and the range parameters; thus, the existing DOA estimation schemes, such as MUSIC and ESPRIT [2], would fail in estimating near-field localization parameters.

By applying the Fresnel approximation to the near-field sources localization, the two-dimensional (2D) MUSIC method, the high-order ESPRIT method, and the path-following method were, respectively, proposed in [3–7] in order to cope with the problem of estimating azimuth and range. In recent years, several methods that can obtain azimuth, elevation, and range have been developed. For instance, in [8], an original 3D higher order statistics based localization algorithm has been presented. By translating the 1D uniform linear array of near-field into a virtual rectangular array of virtual far-field, Challa and Shamsunder [9] proposed a unitary-ESPRIT method. Aberd Meraim and Hua [10] used only the second-order statistics and proposed a higher resolution 3D near-field source localization; however, a parameter-pairing process leading to the poor performance in lower signal-to-noise ratio had to be taken into account. References [11, 12] also proposed fourth-order cumulant based algorithms to estimate 2D DOA and range, but they all coincided with [10]. All abovementioned algorithms rely on the assumption that impinging source signals are stationary; when nonstationary FM signals exist, they will show an unsatisfactory performance.

While quadratic time-frequency distribution [13, 14] has been sought out and properly investigated into sensor and spatial signal processing, and its evaluation of the observation data across the sensor array yields spatial time-frequency distribution (STFDs), the main advantage of STFDs is that it can well handle signals of nonstationary waveforms that are highly localized in the time-frequency domain and effectively improve the robustness of localization methods by spreading the noise power into the whole time-frequency domain. The STFDs based algorithm to locate near-field nonstationary sources has been presented in [15] and showed a satisfactory parameters estimation accuracy; however, it required 2D search and only estimated 1D DOA and range.

In this paper, by exploiting favorable characteristics of a uniform cross array, we present a joint 2D DOA and range estimation algorithm. We first compute five time-frequency matrices to construct a parallel factor (PARAFAC) analysis model. Then, we obtain three-dimensional (3D) near-field parameters via trilinear alternating least squares regression (TALS). Compared with the other methods, the main contribution for the proposed method can be summarized as follows: we obtain 3D near-field sources parameters (elevations, azimuths, and ranges) of nonstationary signals rather than stationary waves; we creatively incorporate STFDs with parallel factor analysis to well avoid both parameter pairing and multidimensional search.

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 elements with interelement spacing . Let the array center be the phase point; the signals received by the th and the th can be, respectively, expressed as

with being additive noise, and

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

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

The th source signal with phase can be modeled as

where and is the amplitude of the th source signal.

2.2. Assumption of Signal Model

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

(1)The source signal is narrowband, independent, and nonstationary process.(2)The additive noise is spatially white Gaussian with zero-mean and independent from the source signals.(3)For unique estimation, we require , , and as well as .

2.3. Parallel Factor Analysis

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

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

where denotes the element of matrix , and similarly for the others. Equation (5) 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 [16]). 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 [16]). Consider a three-dimensional tensor as defined in (5), 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. Computation of the Spatial Time-Frequency Distribution Matrices

The discrete form of Cohen’s class of time-frequency distribution of a signal can be expressed as

where is the time-frequency kernel and the superscript denotes complex conjugate. Replacing by the data snapshot and , we obtain

Substituting (1) into (8), can be extended to the following form:

Under the assumptions and , it is obvious that

where indicates the STFDs of source .

Using a rectangular window of old length , the pseudo Wigner-Ville distribution (PWVD) of is given by

Assume that the third-order derivative of the phase can be negligible over the rectangular window length , and , ; then we obtain the approximated expression as

We construct matrix with the th element being given by

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

And their matrices form becomes

where , , , , , and .

3.2. Construction of the Parallel Factor Analysis Model

Considering the situation of limited samples, we build a parallel factor analysis model that uses the spatial time-frequency distribution as

Letting , the Khutrr-Rao product [16] for (16) shows

where

denoting a row vector consisting of diagonal matrix .

Similarly, (17) also yields

3.3. Estimation of 2D Direction-of-Arrival and Range

As it stands, and are both Vander-monde 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 get each pair as follows:

Finally, the sources parameters can be estimated as

4. Computer Simulation Results

In this section, we explicit several simulation results to evaluate the performance of proposed method. For all examples, a symmetrical cross array with a number of 17 elements and interelements spacing of , where is the wavelength of the narrowband source signals. The noise used in this section is zero-mean, Gaussian distributed, and temporally white, and the root mean square error (RMSE) is defined as (23). All the following presented results are obtained by averaging the results of 200 independent Monte Carlo simulations. Consider

where denotes the number of independent Monte Carlo simulations.

In the first example, we examine the performance of the elevation, azimuth, and range estimations accuracy versus the SNR. The snapshot number is set at 512. Two linear frequency-modulated signals arrival at the sensor array with start and end frequencies and , their 3D near-field parameters locate at and . For the comparison, the fourth-order cumulant based method [9] is also displayed, in which the source signals are non-Gaussian and stationary process. When the SNR varies from 0 dB to 25 dB, the RMSEs of the 2D direction-of-arrival and range estimations using the proposed method and the fourth-order cumulant based method are shown in Figure 2. From Figure 2, we can see that the proposed method outperforms the fourth-order cumulant based method in elevation and azimuth as well as range estimation for all available SNRs. In addition, the RMSE of range estimations for the first source that is closer to the array is less than the second one. This phenomenon is in well agreement with the theoretical analysis that the sources closer to sensor array would hold a smaller standard deviation than the ones far away from the array.

(a) Elevation

(b) Azimuth

(c) Range

Figure 2: The RMSEs of the elevation, azimuth, and range estimation using the proposed method and the fourth-order cumulant based method versus SNRs.

In the second example, the proposed method is used to deal with the situation that two near-field FM signals are impinging on the sensor array shown in Figure 1. The elevation, azimuth, and range of the near-field sources are located at and . Moreover, the snapshot number and SNR are set at 512 and 10 dB. Table 1 establishes the mean and variance of the elevation, azimuth, and range estimations using the proposed method. From Table 1, we can see that the proposed method shows a satisfactory performance in localizing the 3D near-field nonstationary sources.

Table 1: The mean and variance of the proposed method for the second example.

True
Mean
Variance

Source 1
Elevation (°)
35
34.1036
0.0854

Azimuth (°)
40
38.8652
0.0514

Range
1/6
0.1956

Source 2
Elevation (°)
20
19.8274
0.0963

Azimuth (°)
60
59.3610
0.1827

Range
0.4
0.4449

In the last example, we consider the situation when far-field and near-field nonstationary sources are incoming on the sensor array mentioned above, and they are located at and , respectively. The snapshot number is fixed at 512. When the SNR varies from 0 dB to 25 dB, the RMSEs of the 2D direction-of-arrival and range estimations using the proposed method are shown in Figure 3. In addition, the results of the fourth-order cumulant based method are also displayed in the same figure for comparison. From Figure 3, it can be seen that the proposed method still performs a satisfactory estimation accuracy for the case that far-field and near-field sources exist simultaneously. However, the fourth-order cumulant based method shows a poor performance in the same situation.

(a) Elevation

(b) Azimuth

(c) Range

Figure 3: The RMSEs of the elevation, azimuth, and range estimation using the proposed method and the fourth-order cumulant based method versus SNRs.

5. Conclusion

We have developed a spatial time-frequency distribution based algorithm for 3D near-field nonstationary source localization problems. Additionally, with parallel factor analysis technique, there is no parameter pairing or multidimensional searching. Finally, the computer simulation results indicate that using spatial time-frequency distribution and parallel factor together significantly solves the problem of the joint estimation of elevation, azimuth, and range of nonstationary signals. However, the spatial time-frequency averaging methods may lead to the additional computation load.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

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: A Review Journal, vol. 19, no. 4, pp. 596–611, 2009.

H. Krim and M. Viberg, “Two decades of array signal processing research,” IEEE Signal Processing Magazine, vol. 13, no. 4, pp. 67–94, 1996.

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.

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. J. Zhi and M. Y. Chia, “Near-field source localization via symmetric subarrays,” in Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '07), vol. 2, pp. II-1121–II-1124, Honolulu, Hawaii, USA, April 2007.

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

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. N. Challa and S. Shamsunder, “High-order subspace-based algorithms for passive localization of near-field sources,” in Proceedings of the 29th Asilomar Conference on Signals, Systems and Computers, pp. 771–781, 1995.

R. N. Challa and S. Shamsunder, “3-D spherical localization of multiple non-Gaussian sources using cumulants,” in Proceedings of the 8th IEEE Signal Processing Workshop on Statistical Signal and Array Processing (SSAP '96), pp. 101–104, June 1996.

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

C. M. Lee, K. S. Yoon, J. H. Lee, and K. K. Lee, “Efficient algorithm for localising 3-D narrowband multiple sources,” IEE Proceedings: Radar, Sonar and Navigation, vol. 148, no. 1, pp. 23–26, 2001.

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.

M. G. Amin and Y. Zhang, “Direction finding based on spatial time-frequency distribution matrices,” Digital Signal Processing: A Review Journal, vol. 10, no. 4, pp. 325–359, 2000.

Y. M. Zhang, W. F. Ma, and M. G. Amin, “Subspace analysis of spatial time-frequency distribution matrices,” IEEE Transactions on Signal Processing, vol. 49, no. 4, pp. 747–759, 2001.

L. A. Cirillo, A. M. Zoubir, and M. G. Amin, “Estimation of near-field parameters using spatial time-frequency distributions,” in Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '07), vol. 3, pp. III1141–III1144, Honolulu, Hawaii, USA, April 2007.

R. A. Harshman, “Foundation of the PARAFAC procedure: models and conditions for an explanatory multimodal factor analysis,” UCLA Working Papers in Phonetics, vol. 16, pp. 1–84, 1970.