International Journal of Antennas and Propagation
Volume 2018 (2018), Article ID 6012426, 9 pages
https://doi.org/10.1155/2018/6012426

Research Article

Joint High-Resolution Range and DOA Estimation via MUSIC Method Based on Virtual Two-Dimensional Spatial Smoothing for OFDM Radar

Rui Zhang,¹ Ying-Hui Quan,¹ Sheng-Qi Zhu,¹ Lei Yang,² Ya-chao Li,¹ and Meng-Dao Xing¹

¹National Laboratory of Radar Signal Processing, Xidian University, Shaanxi 710071, China
²Tianjin Key Laboratory for Advanced Signal Processing, Civil Aviation University of China, Tianjin 300300, China

Correspondence should be addressed to Ying-Hui Quan; quanyinghui@126.com

Received 29 June 2018; Accepted 24 September 2018; Published 21 November 2018

Academic Editor: Stefano Selleri

Copyright © 2018 Rui Zhang 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.

For the purpose of target parameter estimation of the orthogonal frequency-division multiplexing (OFDM) radar, a high-resolution method of joint estimation on range and direction of arrival (DOA) based on OFDM array radar is proposed in this paper. It begins with the design and analysis of an echo model of OFDM array radar. Since there is no coupling between range and angle parameter estimation for a narrow-band signal, a method which exploits the data of one snapshot to estimate the range and angle of the target by means of multiple signal classification (MUSIC) based on virtual two-dimensional spatial smoothing in range and angle dimensions is devised. The proposed method is capable of joint estimating the range and DOA of the target in a high resolution under a single snapshot circumstance. Simulation experiments demonstrate the validity of the proposal.

1. Introduction

Recently, with the development of hardware technique for complicated waveform generator, it is possible to configure the radar waveform flexibly [1]. Multicarrier signal has been widely applied to radar [2, 3], because it can provide more information about targets. Orthogonal frequency-division multiplexing (OFDM) radar modulates traditional single-carrier transmitting signal to subcarrier signals one by one, and then, it transmits the superposed signal, which greatly improves the performance of radar system and attracts thousands of scholars’ wide attention.

OFDM radar has the following advantages [4]: flexible waveform design, using orthogonal subcarrier signal as subchannel which improves the utilization ratio of the spectrum, enhancing the ability of resisting the pulse noise and channel fading, with large time-bandwidth product, and easy to be processed digitally. Reference [5] firstly considers the subcarrier phase modulation into multicarrier structure of OFDM radar; the spectrum utilization ratio can be improved. In [6], the concept of OFDM applied to the radar system is discussed in detail. Meanwhile, the framework of the OFDM radar system is proposed and the simulation analysis of each module is carried out. It has been proved that OFDM radar has strong abilities of antinoise and antijamming in [6]. The theory of Doppler signal processing for the OFDM radar signal is proposed in [7]. However, Doppler processing on receiver based on correlation function usually encounters the troublesome Doppler ambiguity. Tigrek et al. [8] and Garmatyuk et al. [9] make a feasible study about the influence of OFDM signal on radar imaging and establishment of the communication system network, which lays the foundation of the integration of radar and communication. Using OFDM radar signal to eliminate Doppler ambiguity is presented in [10], and the echoed signal model is established in vector expression. However, it does not take angle information of the target into consideration. Reference [11] utilizes the property that different carrier waveforms of OFDM radar have different responses to the target signal to improve the frequency diversity and solve the problem of multipath reflection. Reference [12] proposes an algorithm of joint estimation on direction of departure (DOD) and direction of angle (DOA) based on multiple-input multiple-output (MIMO) radar. The complexity of the algorithm is decreased by the dimension reduction, and it is suitable to estimate DOD and DOA under the condition of a single snapshot and multiple targets. A DOA estimation method based on subspace is put forward in [13], providing a fundamental for time-varying DOA estimation. References [14, 15] transform the DOA model based on sparse representation of covariance matrix into sparse representation model of a single vector, and the effect of imprecise estimation on noise power is mitigated by the dimension reduction. In [16], a method which determines the range and angle of target is presented, by sending double pulses based on traditional frequency diverse array radar.

In this paper, aiming at target parameter estimation in the OFDM array system, the characteristic of received signal is analyzed and a high resolution method of joint estimation on the range and DOA of the target in OFDM array radar is proposed. Under the narrow-band signal condition, since there is no coupling between the range and angle of the target in the echoed signal of the OFDM radar, the target information is included in the range and angle steering vector, respectively, which can be utilized to jointly estimate the range and angle of signal source. In single snapshot case, a joint high resolution estimation method for the range and the angle of signal source can be obtained by means of the multiple signal classification (MUSIC) algorithm and virtual two-dimensional spatial smoothing technique. In order to evaluate the performance of the proposed method in this paper, the root-mean-square errors (RMSE) of the estimates of the range and angle parameters are examined. Numerical simulation experiments are designed to prove the effectiveness of the proposed method.

2. Echoed Signal Model of OFDM Radar

Taking linear array as an example, it is assumed that the number of array elements is , the distance between two elements is , the propagation speed of electromagnetic wave is , and signal bandwidth is . Considering phased array radar and MIMO (multiple-input and multiple-output) radar system, there are far-field narrow-band signals in space. If the longest time which signal spends on traversing array aperture is much shorter than reciprocal of signal bandwidth, ; the signal can be called “narrow-band signal.”

Under narrow-band signal condition, the time delay between different array elements can be expressed as a formula about phase difference of signal. The response of different array elements to the same signal is to have different phase shifts.

Consider a uniform linear array (ULA) including elements, where each of them transmits OFDM waveform. There are subcarrier channels in the OFDM waveform, and the transmitting signal can be expressed as

where is amplitude, is the initial phase of the th subcarrier, denotes the frequency interval of two adjacent subcarriers, and is the central frequency. Furthermore, in order to ensure the orthogonality between subcarriers, needs to meet with , where is the radar pulse duration. Figure 1 shows the allocation of transmitting antennas and transmitting signals on each element. All the antennas transmit the same signal at the same time. And the signal of each subcarrier channel can be expressed as

Figure 1: Allocation of transmitting antennas.

Under narrow-band circumstance and for arbitrary target signal in space, the downconverted echoes of the OFDM array radar can be written as

where is the target distance, is the target radial velocity, is the speed of light, is the element interval of ULA, denotes DOA of the single point target, and stands for the number of array elements. Note that is sampled discretely, i.e., , where denotes the sampling index. And is steering vector,

where is wave length.

The signal model in the discretization form can be written as

where denotes a data matrix of discreted echoed signal, and

Assuming , and is with the same dimension as the inverse discrete Fourier transform (IDFT) matrix, where stands for the matrix transpose. The useful information about range and angle is contained in , , and . In order to obtain information about range and angle, the signal phase needs to be compensated. Firstly, should be compensated by Fourier transform matrix ; then, next step is phase compensation [10]. The signal phase can be compensated by , where gives the diagonal matrix. After the compensation, the expression of the echoed signal can be written according to (3) as

When the product of matrix and meets (15), it will cause cyclic shift in each row of . In order to compensate for the cyclic shift, Doppler compensation should be completed before the phase compensation.

where is an integer. Each column of the echoed signal in (5) is sequentially processed with Doppler compensation, compensation, and phase compensation. Then, each column of the data that has been compensated is rearranged as a data matrix. Therefore, the space-domain steering vector of OFDM radar echo signal includes direction steering vector in (14) and distance steering vector , that is,

where denotes Kronecker product and can be expressed as

From (17), it can be noted that there is no coupling between the distance and angle parameter estimation for OFDM array radar echo signal under the condition of narrow-band signal, and therefore, the joint estimation on the range and DOA can be carried out.

3. A Method of Joint Estimation on Range and Angle in High Resolution by Virtual Two-Dimensional Spatial Smoothing

It is well known that the Fourier-based beamforming method applies the Fourier transform to the echoed signal to estimate DOA. However, because Fourier transform is limited by Rayleigh limitation, the performance of the beamforming is with high side lobes and limited resolution. To this end, the high-resolution method based on the second-order statistical properties can overcome this limitation. The MUSIC (multiple signal classification) can make an asymptotic and unbiased angle estimation of the target in theory. Under a certain SNR (signal-noise ratio) threshold, the performance of MUSIC is closed to ML (maximum likelihood) estimation. It utilizes the orthogonality between signal subspace and noise subspace to construct spatial spectral function and searches array manifold vector which is orthogonal to noise subspace for the DOA estimation.

Under the condition of narrow-band signal, the distance and angle parameter estimation of target echoed signal in OFDM array radar have no coupling relationship. Therefore, a high-resolution estimation of range and angle can be realized by exploiting the data obtained from a single snapshot. Since the traditional MUSIC algorithm needs many snapshots to estimate covariance matrix, a novel method that uses a virtual two-dimensional spatial smoothing method to sample a single snapshot data is proposed in this paper. It is capable of estimating the covariance matrix by using sampled data which has been smoothed in two-dimensional space.

3.1. The Sampling Method of Virtual Two-Dimensional Spatial Smoothing

Consider the compensated data of one single snapshot. After match filtering and signal separation processing, the data is rearranged into a matrix, and then smooth the matrix spatially by column, and finally divide it into subarrays, and the dimension of a subarry is , where denotes the row number of subarrays and is the column number of subarrays. Then, rearrange all the subarrays into a matrix. Figure 2 shows the virtual two-dimensional spatial smoothing method.

Figure 2: Virtual two-dimensional spatial smoothing method.

The first subarray is

and the steering vector is .

The second subarray is

where is the rotation factor.

The subarray is

The subarray

where is the rotation factor.

The subarray is

where , and is the combined rotational factor,

Therefore, it is equivalent to translate each subarray spatially through above steps. The rest rotation factors can be combined into signal envelope . The data of each subarray can be regarded as sampling data which has been sampled times, that is,

It is reasonable to estimate covariance matrix by which has been two-dimensional spatial smoothed, and it is suitable to realize the joint estimation of range and angle in high resolution.

3.2. Procedures of Estimating Range and Angle by a MUSIC Algorithm

In this section, we can estimate the range and angle of the target by applying the MUSIC algorithm into , which has been smoothed in two-dimensional space in Section 3.1. The detailed processing procedures of the proposed method are given as follows:

(1)The covariance matrix can be estimated from , and

(2)Perform eigenvalue decomposition to covariance matrix (3)Sort eigenvalues and determine the number of signal sources , then extract eigenvectors corresponding to small eigenvalues to form the noise subspace (4)Project two-dimensional search vector onto noise subspace , and calculate spectral peak

(5)The position of the spectral peak is angle and range of the target

3.3. Performance Analysis of Joint Estimation on Range and Angle in High Resolution

In order to test the performance of joint estimation on range and angle in high resolution, this section proposes Cramer-Rao bound (CRB) of parameter estimation on range and angle in OFDM radar.

Consider that estimating distance and angle parameters of a single point source in the condition of Gaussian white noise. Assume that the value of signal-to-noise ratio is , and Fisher information matrix corresponding to angle and distance can be written as

then

where and . Similarly, there are

then .

4. Results and Discussion

The basic parameters of OFDM array radar and target parameters information used in the numerical simulation are shown in Table 1.

Table 1: Simulation parameters of joint estimation on range and angle in high resolution.

Parameter
Value

Center frequency
1 GHz

Channel number ()
10

Sub-carrier number ()
10

Array element distance ()

Frequency interval
100 kHz

Monte Carlo times ()
100

Pulse duration time

4.1. Performance Comparison of Range and Angle Estimation with Multiple Targets

Aiming to verify the effectiveness of range and angle estimation method in OFDM radar, Figure 3 shows the results of target parameter estimation, respectively, using two-dimensional (range-angle) common beamforming and two-dimensional MUSIC algorithm. In simulation experiments, assume that there are three noncoherent point targets, and the angle is, respectively, 14°, 41°, and −18°. The distance is, respectively, 30,550 m, 30,200 m, and 29,700 m, and of three point targets is 10 dB. Other parameters are shown in Table 1. Note that the experiment is carried out in the case of Gaussian white noise.

(a) Frequency spectrum of normalized common beamforming

(b) Frequency spectrum of normalized MUSIC

Figure 3: Spectrum comparison of two methods after normalized.

From Figure 3, since the common beamforming is limited by the aperture of the physical antenna, the resultant beam width is wide and the ability of target resolution is low. Because the target data is sampled for only once, the noise has a great influence on side lobe and target parameter estimation is effected badly by SNR. On the contrary, the MUSIC algorithm utilizes virtual two-dimensional smoothing technique. Although it loses certain array aperture and distance “aperture,” it has a good performance of noise statistics. According to the results of power spectrum estimation, the side lobe is very low and SNR effect on target parameter estimation is deducted. Additionally, the MUSIC algorithm can break the physical limitation of antenna aperture, and its resolution is higher than the common beamforming method. Therefore, MUSIC has a better ability of resolution in the estimation of range and angle.

4.2. Comparison about Spectrum Width of a Single Point Target with a Single Snapshot

For quantitative description of target parameter estimation accuracy, Figure 4 shows the comparing results of using the two-dimensional common beamforming method and virtual two-dimensional smoothing spectrum estimation in a high-resolution method, respectively, to estimate the range and angle in OFDM radar. Simulation parameters are shown in Table 1.

(a) Variation of spectral width with angle

(b) Variation of spectral width with distance

Figure 4: Spectrum comparison of two methods after normalized.

According to simulation experiments, whether in angle domain or in range domain, the spectrum peaks of the MUSIC algorithm are sharper than that of the common beamforming method. The spectrum peaks of the MUSIC algorithm almost have no side lobes. However, the side lobes are high by using the common beamforming method. The results not only demonstrate the effectiveness of analysis about echoed signal properties of OFDM radar but also show that in a single snapshot case, the signal space that is estimated by the sampled data, which is obtained from a virtual two-dimensional smoothing method, is accurate. It is successful to break the limitation that the performance of the MUSIC algorithm depends on a large number of snapshots. Although a virtual two-dimensional smoothing method loses a part of apertures, the resolution ability of target parameters can still reach 5 times more than that of the two-dimensional common beamforming method, as shown in Figure 4.

4.3. The Performance Analysis of Joint Estimation on Range and Angle in High Resolution of a Single Point Target with a Single Snapshot

Monte Carlo experiments are carried out to estimate angle and distance . It is obvious that RMSE of angle is , and RMSE of distance is , where is the number of Monte Carlo experiments, and are, respectively, the estimated angle and range value of th Monte Carlo experiment, and are, respectively, the true angle and range of the target, and RMSE functions of and vary with SNR.

Exploiting the simulation parameters in Table 1 to simulate how the RMSEs, which are estimated by using a high-resolution method of range and angle, varies with SNR. Then, compare the above simulation performance with the results when using the ML method to estimate RMSE and the outcomes that how Cramer-Rao bound varies with SNR, respectively. Figure 5 shows the comparison results.

(a) Performance comparison of angle estimation

(b) Performance comparison of range estimation

Figure 5: Performance comparison of range and angle joint estimation.

In a Gaussian white noise condition, the ML method is the most optimal estimation algorithm. It is found from Figure 5 that the performance of MUSIC is close to the ML method. That is to say, although the MUSIC algorithm is suboptimal in a Gaussian white noise case, its estimation accuracy loss can be negligible. Also, the MUSIC algorithm has a high resolution for the range and angle estimates. It is seen that both the performances of two algorithms are awful when is below −5 dB. The accuracy of two algorithms is close to Cramer-Rao bound when SNR is from −5 dB to 15 dB. Because with a low , the estimation performance of two methods is influenced greatly by , that is to say, the application of two methods needs a certain threshold; the two algorithms are not close to Cramer-Rao bound when is above 15 dB for the reason that there is inherent error caused by angle and distance increment. Therefore, under the condition of a certain value, as long as the step size in search is small enough, the MUSIC algorithm can make an accurate estimation of distance and angle.

4.4. Analysis of Speed and Computational Complexity for MUSIC

The first step of MUSIC is to compute covariance matrix . According to Section 3.1, the dimension of is so that the computational complexity of this step is . The second step is to perform eigenvalue decomposition. The computational complexity of eigenvalue decomposition is [17]. The third step is to sort the eigenvalues which have been obtained from the second step, and the average computational complexity is when utilizing bubble sort [18]. The fourth step is to calculate spectral peaks and search. According to formula (27) above, the computational complexity of every spectral peak searching is . Aiming to two-dimensional searching of and to meet the demand of real-time processing when it is put into practice, the proper computing choices include FPGA (field programmable gate array) and GPU (graphics processing units) for their parallel computing ability. The real-time processing speed depends on the processor type and the range of spectral peak searching [19, 20].

5. Conclusion

In this paper, aiming at the problem of radar target parameter estimation in the OFDM array system, the property of received signal has been analyzed and a method about high-resolution joint estimation on range and angle in OFDM radar is proposed. Under narrow-band signal condition, there is no coupling between the range and angle of a target in the echoed signal. The information about the range and angle is included, respectively, in its steering vector, which can be used to jointly estimate range and angle of the target. Under the single snapshot condition, combining the MUSIC algorithm with a virtual two-dimensional spatial smoothing method can jointly estimate range and angle in high resolution. In the background of the Gaussian white noise, the accuracy of the proposed method is close to that of the common beamforming method. Although the high-resolution method which utilizes a virtual two-dimensional smoothing method loses a part of aperture under the condition of a single snapshot, it still has high resolution in distance and angle. In the future, further investigation on the joint high-resolution estimation of OFDM radar based on sparse recovery will be performed.

Abbreviations

OFDM: | Orthogonal frequency-division multiplexing
DOA: | Direction of arrival
MUSIC: | Multiple signal classification
DOD: | Direction of departure
RMSE: | Root-mean-square errors
ULA: | Uniform linear array
IDFT: | Inverse discrete Fourier transform
SNR: | Signal-noise ratio
ML: | Maximum likelihood
CRB: | Cramer-Rao bound
FPGA: | Field programmable gate Array
GPU: | Graphics processing units.

Data Availability

We have presented some persuasive simulation results as figures based on MATLAB. Requests for access to these detailed data which used to support the findings of this study should be made to Ying-Hui Quan with email addresses yhquan@mail.xidian.edu.cn and quanyinghui@126.com.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

Rui Zhang, Ying-Hui Quan, and Sheng-Qi Zhu designed the algorithm scheme. Rui Zhang and Ying-Hui Quan performed the experiments and analyzed the experiment results. Lei Yang, Ya-chao Li, and Meng-Dao Xing contributed to the manuscript drafting and critical revision. All authors read and approved the final manuscript.

Acknowledgments

This work is partially supported by the Fund for Foreign Scholars in University Research and Teaching Programs (the 111 Project) (no. B18039). This work was supported by the National Natural Science Foundation of China under Grant nos. 61772397, 61303035, and 61471283 and Natural Science Foundation of Tianjin, China under Grant nos. 16JCYBJC41200 and 20162898.

References

R. Calderbank, S. Howard, and B. Moran, “Waveform diversity in radar signal processing,” IEEE Signal Processing Magazine, vol. 26, no. 1, pp. 32–41, 2009.

N. Levanon, “Multifrequency complementary phase-coded radar signal,” IEE Proceedings - Radar, Sonar and Navigation, vol. 147, no. 6, pp. 276–284, 2000.

G. Lellouch, P. Tran, R. Pribic, and P. van Genderen, “OFDM waveforms for frequency agility and opportunities for Doppler processing in radar,” in 2008 IEEE Radar Conference, pp. 1–6, Rome, Italy, May 2008.

Y. Liu, G. Liao, J. Xu, Z. Yang, and Y. Zhang, “Adaptive OFDM integrated radar and communications waveform design based on information theory,” IEEE Communications Letters, vol. 21, no. 10, pp. 2174–2177, 2017.

B. J. Donnet and I. D. Longstaff, “Combining MIMO radar with OFDM communications,” in 2006 European Radar Conference, pp. 37–40, Manchester, UK, September 2006.

N. N. S. S. R. K. Prasad, V. Shameem, U. B. Desai, and S. N. Merchant, “Improvement in target detection performance of pulse coded Doppler radar based on multicarrier modulation with fast Fourier transform (FFT),” IEE Proceedings - Radar, Sonar and Navigation, vol. 151, no. 1, pp. 11–17, 2004.

R. F. Tigrek, W. J. A. D. Heij, and P. V. Genderen, “Multi-carrier radar waveform schemes for range and Doppler processing,” in 2009 IEEE Radar Conference, pp. 1–5, Pasadena, CA, USA, May 2009.

R. F. Tigrek, W. J. A. D. Heij, and P. V. Genderen, “Solving Doppler ambiguity by Doppler sensitive pulse compression using multi-carrier waveform,” in 2008 European Radar Conference, pp. 72–75, Amsterdam, Netherlands, October 2008.

D. Garmatyuk, J. Schuerger, Y. T. Morton, K. Binns, M. Durbin, and J. Kimani, “Feasibility study of a multi-carrier dual-use imaging radar and communication system,” in 2007 European Radar Conference, pp. 194–197, Munich, Germany, October 2007.

R. F. Tigrek, W. J. A. De Heij, and P. Van Genderen, “OFDM signals as the radar waveform to solve Doppler ambiguity,” IEEE Transactions on Aerospace and Electronic Systems, vol. 48, no. 1, pp. 130–143, 2012.

S. Sen and A. Nehorai, “Adaptive OFDM radar for target detection in multipath scenarios,” IEEE Transactions on Signal Processing, vol. 59, no. 1, pp. 78–90, 2011.

Y. Jia, X. Zhong, Y. Guo, and W. Huo, “DOA and DOD estimation based on bistatic MIMO radar with co-prime array,” in 2017 IEEE Radar Conference (RadarConf), Seattle, WA, USA, May 2017.

C. Kusuma and G. Radha, “Robust direction of arrival estimation using subspace methods,” International Journal of Computer Applications, vol. 159, no. 1, pp. 24–29, 2017.

Z. Q. He, Q. H. Liu, L. N. Jin, and S. Ouyang, “Low complexity method for DOA estimation using array covariance matrix sparse representation,” Electronics Letters, vol. 49, no. 3, pp. 228–230, 2013.

T. Qian, J. Z. Xiang, and W. Cui, “Fast covariance matrix sparse representation for DOA estimation based on dynamic dictionary,” in 2016 IEEE 13th International Conference on Signal Processing (ICSP), pp. 138–143, Chengdu, China, November 2016.

W.-Q. Wang and H. Shao, “Range-angle localization of targets by a double-pulse frequency diverse array radar,” IEEE Journal of Selected Topics in Signal Processing, vol. 8, no. 1, pp. 106–114, 2014.

Y. Quan, Y. Li, X. Gao, and M. Xing, “FPGA implementation of real-time compressive sensing with partial Fourier dictionary,” International Journal of Antennas and Propagation, vol. 2016, no. 2, Article ID 1671687, 12 pages, 2016.

O. Astrachan, “Bubble sort: an archaeological algorithmic analysis,” in SIGCSE '03 Proceedings of the 34th SIGCSE technical symposium on Computer science education, pp. 1–5, ACM, New York, NY USA, February 2003.

Z. Zou, W. Hongyuan, and Y. Guowen, “An improved MUSIC algorithm implemented with high-speed parallel optimization for FPGA,” in 2006 7th International Symposium on Antennas, Propagation & EM Theory, pp. 1–4, Guilin, China, October 2006.

M. Kim, K. Ichige, and H. Arai, “Real-time smart antenna system incorporating FPGA-based fast DOA estimator,” in IEEE 60th Vehicular Technology Conference, 2004. VTC2004-Fall. 2004, pp. 160–164, Los Angeles, CA, 2004.