Research Article | Open Access

# Computationally Efficient DOA Tracking Algorithm in Monostatic MIMO Radar with Automatic Association

**Academic Editor:**Frankie KitWing Chan

#### Abstract

We consider the problem of tracking the direction of arrivals (DOA) of multiple moving targets in monostatic multiple-input multiple-output (MIMO) radar. A low-complexity DOA tracking algorithm in monostatic MIMO radar is proposed. The proposed algorithm obtains DOA estimation via the difference between previous and current covariance matrix of the reduced-dimension transformation signal, and it reduces the computational complexity and realizes automatic association in DOA tracking. Error analysis and Cramér-Rao lower bound (CRLB) of DOA tracking are derived in the paper. The proposed algorithm not only can be regarded as an extension of array-signal-processing DOA tracking algorithm in (Zhang et al. (2008)), but also is an improved version of the DOA tracking algorithm in (Zhang et al. (2008)). Furthermore, the proposed algorithm has better DOA tracking performance than the DOA tracking algorithm in (Zhang et al. (2008)). The simulation results demonstrate effectiveness of the proposed algorithm. Our work provides the technical support for the practical application of MIMO radar.

#### 1. Introduction

Multiple-input multiple-output (MIMO) radar employs multiple antennas to simultaneously transmit diverse waveforms and utilizes multiple antennas to receive the reflected signals [1–5]. Direction of arrival (DOA) algorithms in MIMO radar have been recently investigated in [6–24], which contain estimation of signal parameters via rotational invariance technique (ESPRIT) algorithms [6–8], Capon algorithms [10, 11], multiple signal classification (MUSIC) algorithms [9, 12–15], parallel factor (PARAFAC) analysis algorithms [16–18], propagator method [19, 20], quaternion method [21, 22], compressive sampling methods [23, 24], and so on. However, the algorithms mentioned above are, generally, used in offline situation, and they are not applicable for tracking moving targets. The online algorithms can be used for real-time application of MIMO radar.

DOA tracking for array antenna has been investigated for a long time, which contains projection approximation and subspace tracking (PAST) algorithm [25], projection approximation and subspace tracking of deflation (PASTd) algorithm [26], Bi-iterative least-square method [27], Bi-iteration single value decomposition [28], and others [29, 30]. The DOA tracking algorithms in [25–30] have a high computational complexity. A low-complexity method was proposed in [31, 32] to track DOA of moving sources using the elements of the covariance matrix of the received signal in array signal processing, and the DOA tracking algorithm can implement automatic association, which is a key technique in DOA tracking [33].

DOA tracking for MIMO radar is to track the DOA of the moving targets. PARAFAC adaptive algorithms [17, 18] and PASTd [34] were used for DOA tracking for MIMO radar, but an extra data association is required. Kalman-PASTd was proposed in [35] for DOA tracking for monostatic MIMO Radar with automatic association. PARAFAC adaptive algorithm, PASTd algorithm, and Kalman-PASTd algorithm have high computational complexity. In the paper, we propose a computationally efficient DOA tracking algorithm in monostatic MIMO radar with automatic association.

In this paper, we reference the array-signal-processing DOA tracking idea in [31] to propose a DOA tracking algorithm which is suitable for MIMO radar. Using the reduced-dimension transformation for the received signal of MIMO radar, we obtain the covariance matrix of reduced-dimension transformation signal, and then we adopt an improved version for DOA tracking. The proposed algorithm can realize automatic association in DOA tracking. Error analysis and Cramér-Rao lower bound (CRLB) are also derived in this paper. Finally, the simulation results demonstrate the effectiveness and robustness of the proposed algorithm.

There are some differences between the work in [31] and the proposed algorithm. Reference [31] proposed an effective DOA tracking algorithm in array signal processing, while we address DOA tracking problem for MIMO radar in the paper. The DOA tracking algorithm in array signal processing in [31] requires the direction matrix of Vandermonde form. The direction matrix in MIMO radar is not a Vandermonde matrix, and the DOA tracking algorithm in [31] cannot be used directly for DOA tracking in MIMO radar. We employ the reduced-dimension transformation for the received signal to obtain the reduced-dimension direction matrix of Vandermonde form. Our work improves the DOA tracking algorithm in [31] to enhance the DOA tracking performance in MIMO radar, since it fully uses the Toeplitz matrix property to eliminate the noise. Therefore, the proposed algorithm not only can be regarded as an extension of the work in [31], but also is an improved algorithm. Simulation results show that the proposed algorithm has much better DOA tracking performance than the DOA tracking algorithm in [31].

The reminder of this paper is structured as follows. Section 2 develops the data model for monostatic MIMO radar. Section 3 establishes our DOA tracking algorithm based on the elements of the covariance matrix of the reduced-dimension transformation signal. In Section 4, error analysis and Cramer-Rao lower bound (CRLB) are derived. In Section 5, simulation results are presented to verify the effectiveness of the proposed algorithm, while the conclusions are made in Section 6.

*Notation*. , and denote transpose, conjugate-transpose, inverse, and pseudoinverse operations, respectively. stands for diagonal matrix whose diagonal is a vector is a identity matrix; , , and are the Kronecker product, Khatri-Rao product, and Hadamard product, respectively. denotes the expectation operator. denotes the estimated value of .

#### 2. Data Model

We consider a monostatic MIMO radar system equipped with both uniform linear arrays (ULAs) for the transmit/receive array, in which elements and elements are arranged with half-wavelength spacing between adjacent antennas, respectively. The array structure of monostatic MIMO radar is shown in Figure 1. We assume that there are uncorrelated targets. At time , the output of the matched filters at the receiver can be expressed as [9] where is the DOA of the th target at time . . , and with being Doppler frequency and being the amplitude. We assume that and are uncorrelated for the different targets. is the received additive white Gaussian noise with noise vector of zeros mean and covariance matrix an noise vector. stands for Kronecker product. The matrix is During the interval , is a constant and snapshots of sensor data are available for the signal processing.

#### 3. DOA Tracking Algorithm

The DOA tracking algorithm of array signal processing in [31] requires the direction matrix of Vandermonde form. The direction matrix in (2) is not a Vandermonde matrix, and the DOA tracking algorithm of array signal processing in [31] cannot be used directly for DOA tracking in MIMO radar. We employ the reduced-dimension transformation for the received signal to obtain the reduced-dimension direction matrix of Vandermonde form. In this paper, we reference the array-signal-processing DOA tracking idea in [31] to propose a low complexity DOA tracking algorithm which is suitable for MIMO radar.

##### 3.1. Reduced-Dimension Transformation

can be expressed by [9] where , and

Then we define as follows: Using the reduced-dimension transformation for the receiver signal , we obtain where

is the direction matrix of Vandermonde form. Since the reduced-dimension matrix is sparse, its transformation adds less computational load.

The covariance matrix of in (6) is . Consider where = and with .

##### 3.2. DOA Tracking

The direction matrix in the reduced-dimension signal in (6) is a Vandermonde matrix, and we use the improved version of the DOA tracking in [31] for DOA tracking of MIMO radar.

We assume that is slowly varying. The DOA of the targets at time is . Similarly, the DOA at time is , with being the DOA of the th target at time . We define where with .

We define as the covariance matrices of the signal at time . The covariance matrix is where is the direction matrix at time and . is the covariance matrix of the noise.

Then we can obtain [31, 36] is denoted by

We assume that the noise covariance matrix at time is approximately equal to that at time , and then we have

Using the Vandermonde characteristic of the matrices and , the noiseless in (13) can be denoted by [31, 32, 36]

where where is the element of the matrix and . in (15) can be expanded as . Considering that is very small,

Then in (15) can be denoted by

Considering that is very small, . And then in (17) can be rewritten as

According to (18), we can construct the following matrix:

Equation (19) can be rewritten with matrix form

From (14), we find that is a Toeplitz matrix, whose elements in a straight line paralleled to the principal diagonal are equal. We take the following processing to estimate and () to eliminate the noise. Consider where is the element of the matrix . In the DOA tracking algorithm in [31], and are obtained through and , respectively. The proposed algorithm fully uses the Toeplitz matrix property of to eliminate the noise and improves the parameter estimation performance.

Using the method of least square, we get where .

Through the above analysis, the angles at time of and are automatically associated, and the proposed algorithm avoids an extra data association.

For the finite samples, the covariance matrix in (8) is estimated by , and we also get , which is the estimate of the covariance matrix at time . We estimate and as follows:where . We use (22) to estimate , where , whose element is estimated via (23a) and (23b).

Till now, we show the major steps of DOA tracking algorithm for monostatic radar as follows.

*Step 1. *Use reduced-dimension matrix for the received signal and estimate the covariance matrix .

*Step 2. *Get covariance matrix at time , and we calculate .

*Step 3. *We compute and via (23a) and (23b) and construct .

*Step 4. *We estimate via (22), and the DOA at time is .

*Step 5. *Repeat Steps 1–4 to estimate DOA at next time.

*Remark 1. *In this paper, the number of targets in MIMO radar is assumed to be preknown. If we have no knowledge about it, we may use the existing source-number estimation technique in [37] to obtain an estimate of the number of targets.

*Remark 2. *The initial angles in DOA tracking are obtained by ESPRIT algorithm or other DOA algorithms. The initial DOA is , and we get angles .

*Remark 3. *We assume that the noise covariance matrices of adjacent time are approximately equal. The noise component in (12) will be eliminated no matter which type of noise. The proposed algorithm still works well for the case of the colored noise.

*Remark 4. *The proposed algorithm works well in condition of small value of . When becomes large, the proposed algorithm may fail to work.

##### 3.3. Complexity Analysis and Advantages of the Proposed Algorithm

Since the reduced-dimension matrix is sparse, its transformation adds less computational load. The proposed algorithm does not require eigenvalue decomposition of the covariance matrix and avoids an extra data association. For the proposed algorithm, the calculation of the covariance matrix requires , and the computation of needs . The major computational complexity of the proposed algorithm is . Table 1 and Figure 2 show the computational complexity comparison among the proposed algorithm and other DOA tracking algorithms. We find that the proposed algorithm has much lower computational load than PAST in [25] and PASTd in [26].

The advantages of the proposed algorithm can be presented as follows.(1)The proposed algorithm does not require eigenvalue decomposition of the covariance matrix and has lower complexity than the conventional DOA tracking algorithms including PAST and PASTd.(2)The proposed algorithm can implement automatic association of DOA for monostatic MIMO radar.(3)The proposed algorithm has much better DOA tracking performance than the DOA tracking algorithm in [31], which will be shown in Section 5.

#### 4. Error Analysis and CRLB

In this section, we derive the variance of DOA tracking and CRLB. We assume that the observation noise variances are almost the same at adjacent time. When computing , we use some approximate calculations such as and when is smaller. This leads to a little difference compared with the perfect value. Consider where are the high-order expansion terms.

According to (16), (17), (18), (24a), (24b), and (24c), we have

is estimation error of , and can be shown as follows:

According to (26) and (22), the variance of is denoted by

And we have

We define that is the estimation error of . . Then the variance of is

According to [38], we get the variance of where denotes the th row of and is the pseudoinverse of . According to (27a), (27b), and (27c), we get where is shown in (28).

According to [39], we derive the CRLB for DOA estimation in monostatic MIMO-radar at time . Consider where , , , and

. stands for Hadamard product.

Then we can define the average CRLB as follows: where is the total tracking time.

#### 5. Simulation Results

The Monte Carlo simulations are adopted to assess DOA tracking performance of the proposed algorithm. We suppose that there are three moving targets. We define root-mean square error (RMSE) as where is the estimate of DOA of the th Monte Carlo trial. is the times of Monte Carlo trial and . We note that , , and are the number of transmit antennas, receive antennas, and targets, respectively. During the tracking procedure, the targets are tracked over an interval of 100 s with = 1 s. During each 1 s interval, snapshots of sensor data are generated and used to estimate DOA.

Figure 3 depicts the DOA tracking result of the proposed algorithm for , , , , and dB. It is shown that the proposed algorithm can track the DOA in monostatic MIMO radar effectively.

Figures 4-5 show the DOA tracking result of the proposed algorithm for other moving trajectories of targets with , , , , and dB. From Figures 4-5 we find that the proposed algorithm can work well, which proves the robustness of the proposed algorithm.

Figure 6 shows the DOA tracking performance comparison with , , and , where we compare the proposed algorithm against the DOA tracking algorithm in [31], Kalman-PASTd algorithm in [35], and CRLB. It is shown in Figure 6 that the proposed algorithm has better DOA tracking performance than the DOA tracking algorithm in [31], since the proposed algorithm fully uses the Toeplitz matrix property of to eliminate the noise and improves the estimation performance. The proposed algorithm has less DOA tracking performance than Kalman-PASTd algorithm, which has a heavy computational complexity.

Figure 7 investigates the DOA tracking performance of the proposed algorithm with different values of . From Figure 7 we find that DOA tracking performance of the proposed algorithm improves with the decreasing of *. *When increases, the interference will enhance.

Figures 8-9 show DOA tracking performance with different . It is clear that the DOA tracking performance of the proposed algorithm is improved as the number of transmit/receive antennas increases. Multiple transmit/receive antennas improve DOA tracking performance because of diversity gain.

Figure 10 depicts the DOA tracking performance of the proposed algorithm with , , , , and dB, and the different DOA spacing between two consecutive observations is considered. It is shown in Figure 10 that the angle tracking performance of the proposed algorithm degrades with increasing of DOA spacing between two consecutive observations.

#### 6. Conclusions

In this paper we have presented the DOA tracking of multiple moving targets for monostatic MIMO radar. The proposed algorithm obtains DOA estimation of the target via the difference in previous and current covariance matrix of the reduced-dimension transformation signals, and the proposed algorithm reduces the computational complexity and realizes automatic association of DOA. Error analysis and CRLB of DOA tracking are derived in the paper. The simulation results demonstrate effectiveness and robustness of the proposed algorithm. Our research provides technical support for the practical application of MIMO radar.

#### 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 China NSF Grants (61371169, 61301108, and 61071164), Jiangsu Planned Projects for Postdoctoral Research Funds (1201039C), China Postdoctoral Science Foundation (2012M521099 and 2013M541661), Open Project of Key Laboratory of Modern Acoustic of Ministry of Education (Nanjing University), Aeronautical Science Foundation of China (20120152001), Qing Lan Project, Priority Academic Program Development of Jiangsu High Education Institutions, and Fundamental Research Funds for the Central Universities (NS2013024, kfjj130114, and kfjj130115).

#### References

- E. Fishler, A. Haimovich, R. Blum, D. Chizhik, L. Cimini, and R. Valenzuela, “MIMO radar: an idea whose time has come,” in
*Proceedings of the IEEE Radar Conference*, pp. 71–78, April 2004. View at: Google Scholar - X. Li, Z. Zhang, W. Mao, X. Wang, J. Lu, and W. Wang, “A study of frequency diversity MIMO radar beamforming,” in
*Proceedings of the IEEE 10th International Conference on Signal Processing (ICSP '10)*, pp. 352–356, Beijing, China, October 2010. View at: Publisher Site | Google Scholar - J. Li and P. Stoica, “MIMO radar-diversity means superiority,” in
*Proceedings of 14th Annual Workshop Adaptive Sensor Array Processing*, Lexington, Mass, USA, June 2006. View at: Google Scholar - J. Li and P. Stoica, “MIMO radar with colocated antennas,”
*IEEE Signal Processing Magazine*, vol. 24, no. 5, pp. 106–114, 2007. View at: Publisher Site | Google Scholar - A. M. Haimovich, R. S. Blum, and L. J. Cimini, “MIMO radar with widely separated antennas,”
*IEEE Signal Processing Magazine*, vol. 25, no. 1, pp. 116–129, 2008. View at: Publisher Site | Google Scholar - C. Jinli, G. Hong, and S. Weimin, “Angle estimation using ESPRIT without pairing in MIMO radar,”
*Electronics Letters*, vol. 44, no. 24, pp. 1422–1423, 2008. View at: Publisher Site | Google Scholar - J. Li and X. Zhang, “Performance analysis of ESPRIT for angle estimation in bistatic MIMO radar,” in
*Proceedings of the International Conference on Computer, Electrical, and Systems Sciences (CESSE '11)*, pp. 400–403, Wuhan, China, 2011. View at: Google Scholar - M. L. Bencheikh and Y. Wang, “Joint DOD-DOA estimation using combined ESPRIT-MUSIC approach in MIMO radar,”
*Electronics Letters*, vol. 46, no. 15, pp. 1081–1083, 2010. View at: Publisher Site | Google Scholar - M. L. Bencheikh, Y. Wang, and H. He, “Polynomial root finding technique for joint DOA DOD estimation in bistatic MIMO radar,”
*Signal Processing*, vol. 90, no. 9, pp. 2723–2730, 2010. View at: Publisher Site | Google Scholar - X. Zhang and D. Xu, “Angle estimation in MIMO radar using reduced-dimension Capon,”
*Electronics Letters*, vol. 46, no. 12, pp. 860–861, 2010. View at: Publisher Site | Google Scholar - X. Zhang, Y. Huang, C. Chen, J. Li, and D. Xu, “Reduced-complexity Capon for direction of arrival estimation in a monostatic multiple-input multiple-output radar,”
*IET Radar, Sonar and Navigation*, vol. 6, no. 8, pp. 796–801, 2012. View at: Publisher Site | Google Scholar - X. Gao, X. Zhang, G. Feng, Z. Wang, and D. Xu, “On the MUSIC-derived approaches of angle estimation for bistatic MIMO radar,” in
*Proceedings of the International Conference on Wireless Networks and Information Systems (WNIS '09)*, pp. 343–346, Shanghai, China, December 2009. View at: Publisher Site | Google Scholar - J. Li and X. Zhang, “Improved Joint DOD and DOA estimation for mimo array with velocity receive sensors,”
*IEEE Signal Processing Letters*, vol. 18, no. 12, pp. 717–720, 2011. View at: Publisher Site | Google Scholar - X. Zhang, L. Xu, and D. Xu, “Direction of departure (DOD) and direction of arrival (DOA) estimation in MIMO radar with reduced-dimension MUSIC,”
*IEEE Communications Letters*, vol. 14, no. 12, pp. 1161–1163, 2010. View at: Publisher Site | Google Scholar - F. Liu and J. Wang, “AD-MUSIC for jointly DOA and DOD estimation in bistatic MIMO radar system,” in
*Proceedings of the International Conference on Computer Design and Applications (ICCDA '10)*, pp. 455–458, Qinhuangdao, China, June 2010. View at: Publisher Site | Google Scholar - X. Zhang, Z. Xu, L. Xu, and D. Xu, “Trilinear decomposition-based transmit angle and receive angle estimation for multiple-input multiple-output radar,”
*IET Radar, Sonar & Navigation*, vol. 5, no. 6, pp. 626–631, 2011. View at: Publisher Site | Google Scholar - H. Wu and X. Zhang, “DOA tracking in monostatic MIMO radar using PARAFAC-RLST algorithm,” in
*Proceedings of the 3rd International Conference on Information and Systems Engineering (ICISE '11)*, pp. 958–961, Amsterdam, The Netherlands, July 2011. View at: Google Scholar - D. Nion and N. D. Sidiropoulos, “Adaptive algorithms to track the PARAFAC decomposition of a third-order tensor,”
*IEEE Transactions on Signal Processing*, vol. 57, no. 6, pp. 2299–2310, 2009. View at: Publisher Site | Google Scholar | MathSciNet - X. Zhang, H. Wu, J. Li, and D. Xu, “Computationally efficient DOD and DOA estimation for bistatic MIMO radar with propagator method,”
*International Journal of Electronics*, vol. 99, no. 9, pp. 1207–1221, 2012. View at: Publisher Site | Google Scholar - Z. D. Zheng and J. Y. Zhang, “Fast method for multi-target localisation in bistatic MIMO radar,”
*Electronics Letters*, vol. 47, no. 2, pp. 138–139, 2011. View at: Publisher Site | Google Scholar - J. Li and X. Zhang, “2D-angle estimation algorithm using quaternion theory in bistatic MIMO-radar,” in
*Proceedings of the 3rd International Conference on Information Science and Engineering (ICISE '11)*, pp. 752–755, 2011. View at: Google Scholar - J. Li, X. Zhang, and F. Wang, “Quaternion root-MUSIC algorithm for angle estimation in bistatic MIMO radar,”
*Journal of Electronics and Information Technology*, vol. 34, no. 2, pp. 300–304, 2012. View at: Publisher Site | Google Scholar - Y. Yu, A. P. Petropulu, and H. V. Poor, “MIMO radar using compressive sampling,”
*IEEE Journal on Selected Topics in Signal Processing*, vol. 4, no. 1, pp. 146–163, 2010. View at: Publisher Site | Google Scholar - Y. Liu, M. Y. Wu, and S. J. Wu, “Fast OMP algorithm for 2D angle estimation in MIMO radar,”
*Electronics Letters*, vol. 46, no. 6, pp. 444–445, 2010. View at: Publisher Site | Google Scholar - B. Yang, “Projection approximation subspace tracking,”
*IEEE Transactions on Signal Processing*, vol. 43, no. 1, pp. 95–107, 1995. View at: Publisher Site | Google Scholar - B. Yang, “Extension of the PASTd algorithm to both rank and subspace tracking,”
*IEEE Signal Processing Letters*, vol. 2, no. 9, pp. 179–182, 1995. View at: Publisher Site | Google Scholar - S. Ouyang and Y. Hua, “Bi-iterative least-square method for subspace tracking,”
*IEEE Transactions on Signal Processing*, vol. 53, no. 8, part 2, pp. 2984–2996, 2005. View at: Publisher Site | Google Scholar | MathSciNet - P. Strobach, “Bi-iteration SVD subspace tracking algorithms,”
*IEEE Transactions on Signal Processing*, vol. 45, no. 5, pp. 1222–1240, 1997. View at: Publisher Site | Google Scholar - S. B. Park, C. S. Ryu, and K. K. Lee, “Multiple target angle tracking algorithm using predicted angles,”
*IEEE Transactions on Aerospace and Electronic Systems*, vol. 30, no. 2, pp. 643–648, 1994. View at: Publisher Site | Google Scholar - C. K. Sword, M. Simaan, and E. W. Kamen, “Multiple target angle tracking using sensor array outputs,”
*IEEE Transactions on Aerospace and Electronic Systems*, vol. 26, no. 2, pp. 367–373, 1990. View at: Publisher Site | Google Scholar - H. Zhang, L. Zhang, S. Wu, and Y. Liu, “A new multiple targets angle tracking method,”
*Journal of Electronics and Information Technology*, vol. 29, no. 12, pp. 2840–2842, 2007. View at: Google Scholar - H. Zhang, S. Wu, L. Zhang, and Y. Liu, “Multiple targets angle tracking algorithm based on the elements of the covariance matrix,”
*Journal of Xidian University*, vol. 35, no. 5, pp. 785–792, 2008. View at: Google Scholar - C. R. Rao, L. Zhang, and L. C. Zhao, “Multitarget angle tracking an algorithm for data association,”
*IEEE Transactions on Signal Processing*, vol. 42, no. 2, pp. 459–462, 1994. View at: Publisher Site | Google Scholar - H. Wu and X. Zhang, “DOD and DOA tracking algorithm for bistatic MIMO radar using PASTD without additional angles pairing,” in
*Proceedings of the IEEE 5th International Conference on Advanced Computational Intelligence (ICACI '12)*, pp. 1132–1136, October 2012. View at: Publisher Site | Google Scholar - X. F. Zhang, J. F. Li, G. P. Feng, and H. L. Wu, “Kalman-PASTd based DOA tracking algorithm for monostatic MIMO radar,” in
*Proceedings of the International Conference on Information, Services and Management Engineering (ISME '11)*, pp. 220–224, 2011. View at: Google Scholar - W. Hailang,
*Multi-target angle estimation and tracking in monostatic MIMO radar system [M.S. thesis]*, Nanjing University of Aeronautics and Astronautics, 2013. - J. Xin, N. Zheng, and A. Sano, “Simple and efficient nonparametric method for estimating the number of signals without eigendecomposition,”
*IEEE Transactions on Signal Processing*, vol. 55, no. 4, pp. 1405–1420, 2007. View at: Publisher Site | Google Scholar | MathSciNet - B. D. Rao and K. V. S. Hari, “Performance analysis of ESPRIT and TAM in determining the direction of arrival of plane waves in noise,”
*IEEE Transactions on Speech and Signal Processing*, vol. 37, no. 12, pp. 1990–1995, 1989. View at: Publisher Site | Google Scholar - P. Stoica and A. Nehorai, “Performance study of conditional and unconditional direction-of-arrival estimation,”
*IEEE Transactions on Acoustics, Speech, and Signal Processing*, vol. 38, no. 10, pp. 1783–1795, 1990. View at: Publisher Site | Google Scholar

#### Copyright

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