Joint High-Resolution Range and DOA Estimation via MUSIC Method Based on Virtual Two-Dimensional Spatial Smoothing for OFDM Radar
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.
Recently, with the development of hardware technique for complicated waveform generator, it is possible to configure the radar waveform flexibly . 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 : 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  firstly considers the subcarrier phase modulation into multicarrier structure of OFDM radar; the spectrum utilization ratio can be improved. In , 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 . The theory of Doppler signal processing for the OFDM radar signal is proposed in . However, Doppler processing on receiver based on correlation function usually encounters the troublesome Doppler ambiguity. Tigrek et al.  and Garmatyuk et al.  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 , and the echoed signal model is established in vector expression. However, it does not take angle information of the target into consideration. Reference  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  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 , 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 , 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
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 . 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.
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.
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
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
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
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 . 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 . 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].
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.
|OFDM:||Orthogonal frequency-division multiplexing|
|DOA:||Direction of arrival|
|MUSIC:||Multiple signal classification|
|DOD:||Direction of departure|
|ULA:||Uniform linear array|
|IDFT:||Inverse discrete Fourier transform|
|FPGA:||Field programmable gate Array|
|GPU:||Graphics processing units.|
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 [email protected] and [email protected].
Conflicts of Interest
The authors declare that they have no competing interests.
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.
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.
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.View at: Publisher Site | Google Scholar
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.View at: Google Scholar
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.View at: Google Scholar