International Journal of Antennas and Propagation

Volume 2014 (2014), Article ID 684591, 5 pages

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

## The PARAFAC-MUSIC Algorithm for DOA Estimation with Doppler Frequency in a MIMO Radar System

^{1}School of Electronic and Information Engineering, Beihang University, Beijing 100191, China^{2}College of Information Science & Technology, Beijing University of Chemical Technology, Beijing 100029, China

Received 26 February 2014; Accepted 10 June 2014; Published 23 June 2014

Academic Editor: Wei Xu

Copyright © 2014 Nan Wang 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

The PARAFAC-MUSIC algorithm is proposed to estimate the direction-of-arrival (DOA) of the targets with Doppler frequency in a monostatic MIMO radar system in this paper. To estimate the Doppler frequency, the PARAFAC (parallel factor) algorithm is firstly utilized in the proposed algorithm, and after the compensation of Doppler frequency, MUSIC (multiple signal classification) algorithm is applied to estimate the DOA. By these two steps, the DOA of moving targets can be estimated successfully. Simulation results show that the proposed PARAFAC-MUSIC algorithm has a higher accuracy than the PARAFAC algorithm and the MUSIC algorithm in DOA estimation.

#### 1. Introduction

Recently, there has been a growing interest in multiple-input multiple-output radar which utilizes multiple antennas at both transmitter and receiver. According to the array configuration, MIMO radar can be classified into two main types, the first of which is called distributed MIMO radar. It is composed of widely separated transmit antennas which transmit the linearly independent signals. Based on the sufficient distribution space and linearly independent signals, the distributed MIMO radar is able to obtain rich scattering properties of the targets and to mitigate radar cross-section (RCS) fluctuations. The second type is called collocated MIMO radar, in which the transmit elements are collocated. By transmitting independent waveforms and capitalizing on the MIMO spatial signature, the collocated MIMO can estimate the parameters of interest via coherent processing.

On the other hand, an accurate estimation of signal DOA has made a great sense in radar system of military and commercial application. Many algorithms were proposed on DOA estimation in MIMO radar system as it provides waveform diversity and spatial distribution of flexibility. In [1], the Capon algorithm, proposed in 1969, is a main lobe self-adaptive algorithm which is stable but poor at low SNR. ESPRIT (estimation signal parameter via rotational invariance techniques) algorithm in [2] exploits the invariance property for DOA estimation in MIMO radar system. In [3], MUSIC algorithm, one of the most popular algorithms for DOA estimation, utilizes subspace analysis and has a good performance.

Although these algorithms introduced above are useful for DOA estimation, they are poor when the Doppler frequency is taken into consideration. As the target always moves, DOA estimation of moving targets is extremely essential. In [4, 5], PARAFAC algorithm exploits the iteration of TALS (trilinear alternating least square) to estimate the frequency and DOA. PARAFAC algorithm has a better performance than Capon, ESPRIT, and MUSIC when considering Doppler frequency. PARAFAC-MUSIC algorithm, proposed in this paper, utilizes the PARAFAC algorithm to estimate the Doppler frequency and then exploits MUSIC algorithm to estimate the DOA after the Doppler effect is eliminated. Simulation shows that PARAFAC-MUSIC algorithm is also able to solve the problem generated from Doppler frequency and may have a better performance than PARAFAC algorithm.

#### 2. Signal Model

Assume that there is a monostatic MIMO radar system with collocated transmit antennas and collocated receive antennas, both of which are uniform linear arrays (ULA), shown in Figure 1. Since the transmit array and the receive array are assumed to be close to each other, the DOD (direction of departure) and the DOA are approximately equal, denoted by . and are the interelement spacing at the transmitter and receiver. The transmit antennas transmit orthogonal waveforms simultaneously while the receive antennas receive the signals reflected from the targets. It is supposed that there are targets in the far field, each of which can be considered as a point source. And the RCS of each target is assumed as a constant during a pulse period.

The received signal at time before match filters can be expressed as where is the receive steering vector and is the transmit steering vector. denotes the wavelength. Consider , where is the reflection coefficient depending on the RCS of the target , is Doppler frequency of the target , and is the transmit signal of each transmit antenna. is a complex Gaussian white noise vector with zero mean and covariance matrix . represents transpose operator.

As the transmit waveforms are orthogonal from each other, the output of the match filters at time is where , and “” represents the Khatri-Rao product.

So the received signal before the match filters can be expressed by matrix as where and denotes the number of snapshot. is complex Gaussian white noise matrix.

The output at match filters can be expressed by matrix as where is a matrix and is a matrix. Consider , in which is sampling time. is complex Gaussian white noise matrix.

#### 3. PARAFAC-MUSIC Algorithm

PARAFAC-MUSIC algorithm combines the advantages of PARAFAC algorithm and MUSIC algorithm. It exploits PARAFAC algorithm to estimate the Doppler frequency, then eliminates the effect in received signal generated from Doppler frequency, and finally utilizes MUSIC algorithm to estimate the DOA.

As known from formula (5), the output at match filters possesses the character of three-way model. Therefore, it can be expressed as set composed of sections , each of which is a matrix. can be expressed as where represents the diagonal matrix composed of the elements obtained from the th row of the matrix.

In the same way, can be expressed as a set composed of sections or a set composed of sections . and can be, respectively, expressed as So , , and can be, respectively, expressed as

When the three-way array is obtained, it is common to utilize TALS [5] to estimate , , and . TALS is a popular method in data detection of three-way array model which updates one estimated matrix in each step after obtaining the initial estimated matrixes. Least square (LS) algorithm always works in estimated matrix update that takes the estimated matrix to be updated as a variable and other estimated matrixes as constants in each step. When all the estimated matrixes are updated, it will carry on the next iteration until convergence. Specific steps of TALS are stated as follows.

*Step 1. *Construct the three-way array , , based on the output at match filters.

*Step 2. *Initialize all the estimated matrixes , , and .

*Step 3. *Substitute , into formula (9) to obtain shown as formula (10). , represents Frobenius norm and the pseudo inverse operation. represents the iterations

*Step 4. *Substitute , into formula (11) to obtain shown as formula (12). Consider

*Step 5. *Substitute , into formula (13) to obtain shown as formula (14). Then calculate . If ( is error threshold), repeat Step 3 to Step 5; otherwise, go to Step 6. Consider

*Step 6. *, , and are obtained after iteration. Then can be calculated by the following formula:
where presents the element in and and are the indexes of row and column. represents phase obtain operation.

When the Doppler frequency has been estimated, it can be used to eliminate the Doppler effect to received signal . The signal after preprocessing can be expressed as where , . , represents the compensatory phase of target and can be expressed as So the output of match filters can be expressed as where represents Hermit operation and represents vector obtain operation.

The covariance of can be expressed as

Then eigendecomposition is carried on for shown in formula (20). Consider where and denotes eigenvalue of .

Since the first eigenvalues are composed of variance of signal and Gaussian white noise while the remaining eigenvalues are only composed of variance of Gaussian white noise, the first eigenvalues are larger than the rest at high SNR. The first eigenvalues are defined as “signal eigenvalues” and the rest are defined as “noise eigenvalues.” So the eigenvector matrix can be classified into two parts in the following formula: where is composed of the signal eigenvectors and is composed of the noise eigenvectors.

So the spectrum estimation formula of MUSIC is expressed as where . “” represents the Kronecker product.

Then the estimated DOA can be obtained according to the position of the spectral peak.

#### 4. Results and Discussion

In order to state that PARAFAC-MUSIC algorithm has a better performance in DOA estimation with Doppler frequency, simulation is carried on with MATLAB software to compare PARAFAC-MUSIC algorithm with PARAFAC algorithm and MUSIC algorithm.

Assume that there are two targets in the far field and a MIMO radar system with 8 collocated transmit antennas and 8 collocated receive antennas. The frequency of the transmit signals is 4 GHz. The array structure is the same as Figure 1 with space between adjacent transmit elements and space between adjacent receive elements. The targets are at and , respectively, relative to the MIMO radar system with the same velocity and the same RCS . The sampling frequency is 100 KHz. The number of snapshots is . Additionally, results shown below are obtained from 500 Monte Carlo experiments at each SNR.

Figures 2 and 3 show the two targets’ mean errors of the estimated DOAs obtained from PARAFAC algorithm, MUSIC algorithm, and proposed algorithm. It is apparent to find that the mean error of target 1 obtained from MUSIC algorithm fluctuates near and converges on gradually with the increasing SNR, and the mean error of target 2 obtained from MUSIC algorithm fluctuates near and converges on gradually with the increasing SNR. The results show that the MUSIC algorithm cannot provide an accurate estimation because of the Doppler frequency. Nevertheless, as far as the PARAFAC algorithm and the proposed algorithm are concerned, the mean error of two targets fluctuates near and converges on gradually with the increasing SNR. But the mean error of the proposed algorithm in both figures is closer to than the one of PARAFAC algorithm which means that the proposed algorithm has a better performance than the PARAFAC algorithm.

Figures 4 and 5 illustrate the two targets’ root mean square errors (RMSE) of the estimated DOAs obtained from PARAFAC algorithm, MUSIC algorithm, and proposed algorithm. RMSE of the estimated DOAs in both figures decreases gradually with the increasing SNR. Though the RMSE of the estimated DOAs in both figures obtained from the MUSIC algorithm, the PARAFAC algorithm and the proposed algorithm tend to with the increasing SNR; the RMSE of the proposed algorithm is smaller than the other two algorithms which means that the proposed algorithm has a smaller fluctuation in DOA estimation.

Simulation results state that MUSIC algorithm suffers from Doppler frequency seriously and cannot provide an accurate DOA estimation to moving target. Besides, the PARAFAC-MUSIC is able to estimate the DOA with a higher accuracy than the PARAFAC algorithm to moving targets.

#### 5. Conclusions

In this paper, some algorithms in DOA estimation are reviewed first. Then a new algorithm named PARAFAC-MUSIC is proposed to estimate the DOA with Doppler frequency in a MIMO radar system. The proposed algorithm is able to estimate DOAs of multiple targets. Finally simulation results confirm that the MUSIC algorithm is not suitable for DOA estimation with Doppler frequency while PARAFAC-MUSIC algorithm is very popular for DOA estimation with Doppler frequency and has a higher accuracy than the PARAFAC algorithm and the MUSIC algorithm.

#### Conflict of Interests

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

#### Acknowledgment

The authors would like to thank the Pre-research Foundation for the support of their research (no. 9140C800202120C80279).

#### References

- J. Capon, “High-resolution frequency-wavenumber spectrum analysis,”
*Proceedings of the IEEE*, vol. 57, no. 8, pp. 1408–1418, 1969. View at Google Scholar · View at Scopus - C. Duofang, C. Baixiao, and Q. Guodong, “Angle estimation using ESPRIT in MIMO radar,”
*Electronics Letters*, vol. 44, no. 12, pp. 770–771, 2008. View at Publisher · View at Google Scholar · View at Scopus - A. Zahernia, M. J. Dehghani, and R. Javidan, “MUSIC algorithm for DOA estimation using MIMO arrays,” in
*Proceedings of the 6th International Conference on Telecommunication Systems, Services, and Applications (TSSA '11)*, pp. 149–153, October 2011. View at Publisher · View at Google Scholar · View at Scopus - D. Nion and N. D. Sidiropoulos, “A PARAFAC-based technique for detection and localization of multiple targets in a MIMO radar system,” in
*Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '09)*, pp. 2077–2080, April 2009. View at Publisher · View at Google Scholar · View at Scopus - Z. Jianyun, Z. Zhidong, and L. Xiaobo, “Joint DOD, DOA and doppler frequency estimation for bistatic MIMO radar system,”
*Journal of Electronics & Information Technology*, vol. 32, no. 8, pp. 1009–5896, 2010. View at Google Scholar