# A joint Doppler frequency shift and DOA estimation algorithm based on sparse representations for colocated TDM-MIMO radar.

1. IntroductionIt is a novel technology and worth studying deeply that multiple-inputmultiple-output (MIMO) radar architecture is applied to autocruise and safety driving system [1, 2]. Using this technology we hope distance, velocity, and direction of arrival (DOA) of multisource signals reflected by other automobiles near your own car can be measured simultaneously. According to these estimated parameters the auto can take measures to deal with all kinds of situations automatically or warn the driver. MIMO radar architecture with colocated transmitter (TX) and receiver (RX) antennas has a larger virtual aperture [3]. The basic idea of virtual aperture is shown in Figure 1. Here, we consider a uniform linear array (ULA) with M = 4 transmitters, N = 4 receiver antennas, and the receiver array spacing d. According to MIMO theory we know that its virtual array consists of 16 virtual antennas in Figure 1. So it can achieve a more accurate DOA and improve multisource resolution in the same time. In recent years sparse representations and compressed sensing algorithm have been widely studied and some algorithms based on sparse representations have been applied to DOA estimation for MIMO radar [3-7]. DOA estimation algorithm based on sparse representations and reconstruction theory has two important advantages. One is that it can process coherent signal; the other is that it can complete DOA estimation in few snapshots.

Generally speaking, MIMO radars require more complex hardware [8] when they can be realized by code, frequency. MIMO system requires that the number of receiver channels should be equal to the number of receiver antennas. Moreover, a difference in gain and phase between the receiver channels deteriorates the accuracy significantly [8]. Meanwhile customers also demand an increased performance of automotive radar sensor at low sensor costs with the development of SiGe MMIC technology [2, 9]. The architecture of the time division modulated (TDM) for MIMO (TDM-MIMO) radars may be utilized for the purpose of economic and simple system design for autocruise and safety driving system. Some receiver or transmitter channels may be cut only by depending on switching different RX or TX antennas in time domain.

A MIMO radar has the ability to transmit different, noncoherent orthogonal signals from different positions using an array of multiple transmitter antennas [10]. If TX antennas could transmit signal at different times, only one TX antenna is permitted to transmit signal during each period. This MIMO radar system can be called as TDM-MIMO radar system. For example, 4 TX antennas may transmit signal at the schedule in Figure 2. Note that the signal phase related to other receiver channels can be varying even if we are facing one incident signal just because TX antennas transmit signal at different time instants in one cycle and the target moving conduces Doppler frequency shift. The steering vector-model of colocated TDM-MIMO radar is "time varying steering vector" model [8]. Hence, the conventional algorithms such as Capon and MUSIC [11-13] cannot process this kind of input signal model directly. It is also very difficult to precisely estimate DFS in advance in a short time for the practical application of TDM-MIMO radar.

In recent years a new highly integrated SiGe MIMO radar sensor has been developed [2, 14-16]. In these systems the frequency-modulated continuous wave (FMCW) will be radiated by transmitters. So DFS and high resolution DOA estimation technology research are very important requirements for the application in autocruise and safety driving system [16, 17].The specification of the automotive long range radar (LRR) sensor includes a highDOA estimation accuracy of [[theta].sub.err] [less than or equal to] 0.4[degrees], an angular resolution of [DELTA][theta] [less than or equal to] 3[degrees], and target ranges R up to 200m with relative target velocities V between -60 m/s and 20 m/s [18].

Facing the aforementioned applications and problems of parameters estimation for TDM-MIMO radar, some researchers proposed some methods to solve these problems. For time-invariant steering vector model, the conventional algorithms such as Capon beamformer [8], AM [18], and MUSIC [16, 18] were applied. For time varying steering vector model, the Cramer-Rao bound (CRB) of parameter estimation for a TDM-MIMO radar with colocated antennas was computed and an optimal TDM scheme was derived [1, 19]. ISI-SAGE algorithm was tested to estimate DFS and DOA jointly [20]. A new nonstationary version of Capon beam forming was developed for the time varying steering vector model [8].

In this paper, we offer a new insight into the joint DFS And DOA estimation based on sparse representations for time varying steering vector model of colocated TDM-MIMO radar. An ultracomplete redundancy dictionary for angle-frequency space is founded in order to complete sparse representations for the input signal. The SVD-SR algorithm which stands for joint estimation based on sparse representations using SVD decomposition with OMP algorithm and the improved M-FOCUSS algorithm which combines the classical M-FOCUSS algorithm with joint sparse recovery spectrum are applied to the new model to solve MMV problem. By using these algorithms, both the coherent signals and uncorrelated signals can be distinguished successfully and estimation accuracy can be further improved in the case of MMV with low SNR.

This paper is concerned with issues of signal model and signal processing for TDM-MIMO. The remainder of this paper is organized as follows. We present the signal model in Section 2 firstly. Then sparse representations' principle for TDM-MIMO model is discussed and angle-frequency dictionary is designed in Section 3. How to deal with joint DFS and DOA estimation based MMV model using SVD decomposition and improved M-FOCUSS algorithm is demonstrated in detail in Section 4. In Section 5 simulation results and some discussions are given. Finally, we conclude our work in Section 6.

2. Signal Model

We consider a colocated MIMO radar consisting of a linear receiver and transmitter array with isotropic antennas. The moving target is modeled as a point source and the transmitted signal is narrowband. The positions of the [N.sub.RX] RX antennas and [M.sub.TX] TX antennas are given in [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], respectively, in units of 2[pi]/[lambda], where [lambda] is the carrier wavelength. The positions of the TX antennas in the sequence in which they transmit are given in [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. [N.sub.Pulse] is the number of transmitted pulses in a cycle. The time instance at which the antennas transmit in a cycle is given [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. This means that the antenna with position [d.sup.Pulse.sub.i] transmits at time [t.sub.i]; see Figure 2 for an example. So the steering vector of the virtual array is given by

[a.bar](u) = exp (j x [[d.bar].sup.Virt]u) = exp (j x [[d.bar].sup.Pulse]u)[cross product] exp(j x [[d.bar].sup.RX]u). (1)

Here u = sin([THETA]) is the DOA of targets. [y.bar] = exp([x.bar]) is understood as an element-by-element operation [y.sub.i] = exp([x.sub.i]). [cross product] is the Kronecker tensor product. -stands for column vector.

We assume a number of L measurement cycles where one cycle consists of [N.sub.Pulse] pulses. See Figure 3 for an example. [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. So the complex baseband signal [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] of cycle l is given as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (2)

Here [??] is the entrywise Hadamard product and [N.sub.Virt] = [N.sub.Pulse] x [N.sub.RX] is the number of virtual antennas. s(l) [member of] C is the target echo signal. [N.bar](l) is additive complex Gaussian noise with zero mean, spatially and temporally uncorrelated with E([N.bar](l)[[N.bar].sup.H] (m)) = [[delta].sub.lxm][[sigma].sup.2]I ([sup.H] stands for Hermitian transformation). The factor 1/[square root of [N.sub.Pulse]] adjusts the signal strength of one cycle, consisting of [N.sub.Pulse] pulses, according to a constant transmitting energy. [[omega].sub.d] stands for Doppler frequency shift. Consider

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (3)

Here [[1.bar].sub.k] is a vector of length [N.sub.RX] with all elements equal to 1. With the definition of a new steering vector,

[b.bar](u, [[omega].sub.d]) := exp (-j x [[gamma].bar][[omega].sub.d]) [??] [a.bar] (u). (4)

the signal model can be written as

[y.bar](l) = [b.bar](u, [[omega].sub.d]) x 1 [square root of [N.sub.Pulse]]s (l)+ [N.bar](l), (5)

From formula (4), we will see that steering vector for input model of colocated TDM-MIMO radar is "time varying steering vector" model since there is [bar.t] item existing. In order to complete the parameters estimation wemake the following assumptions.

(1) The far field condition is satisfied. So DOA of one target echo signal is the same for all receiver antennas.

(2) The DOA and the moving velocity of the target will not change during one complete measurement. So both DOA and DFS are constants during L measurement cycles. In other words, the parameters' changes have smaller effect on the accuracy of parameters estimation and these changes can be ignored.

(3) We consider a MIMO radar system with colocated antennas; that is, the direction of departure (DOD) and DOA are the same.

(4) Time of switching different TX and RX antennas or channels is accurate; that is, the estimation error of [t.bar] vector is ignored.

3. Sparse Representations for TDM-MIMO Receiver Signal

The signal sparse representations mean that a signal may be represented in an ultracomplete redundancy dictionary. Then we find the best linear combination of atoms to represent original signal. Inspired by the literature [7], we

can establish an ultracomplete redundancy dictionary called angle-frequency space to represent original signal [y.bar](l) of formula (5).The following is the angle-frequency dictionary [PSI] we designed:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (6)

According to formula (6) we can rewrite the signal model of formula (5) to group L snapshots as

Y = [PSI] x Z + N, (7)

where Y = [[y.bar](1), [y.bar] (2), ..., [y.bar] (L)] is a matrix of size [N.sub.Virt] x L, Z = [[z.bar] (1), [z.bar] (2), ..., [z.bar] (L)] is a matrix of size PQ x L, and N is a noise matrix of size [N.sub.Virt] x L. Let [[[phi].bar].sub.1], [[[phi].bar].sub.2], ..., [[[phi].bar].sub.pq] be a set of joint angle-frequency spaces of DFS and DOA space. Here PQ is the total number of vectors of joint angle-frequency space matrix that the algorithm can scan. Assuming that K incident signals exist in the noise-free case, we can find K nonzero coefficients [z.sub.[theta]](t) = [s.sub.k](t) at positions corresponding to the atom [theta] = [[theta].sub.k], (k = 1, 2, ..., K, [[theta].sub.k] = [[[phi].bar].sup.q,p] = index(DOA([s.sub.k])), q = index(DFS(sk))) and zero coefficients at the remaining PQ - K; that is, the sparsity in the angle frequency space denotes that only a few atoms from the dictionary will be required to match the measurements. Here K [much less than to] PQ is needed to be satisfied.

The number of nonzero rows of Z matrix is referred to as the sparsity of the solution. There are many recovery algorithms to reconstruct original signal. Popular approaches include greedy sequential search techniques such as matching pursuit and methods based on minimizing diversity measures such as basis pursuit and FOCUSS [21]. Firstly we use orthogonal matching pursuit (OMP) algorithm [22] described as follows.

Input. The observation vector [y.bar] [member of] [C.sup.n] and a measurement matrix [PSI] = {[[[phi].bar].sub.i], i = 1, 2, ..., m} [member of] [C.sup.nxm].

Initialization. Index I = 0, residual [r.bar] = [y.bar], and sparse representation [??] = 0 [member of] [C.sup.m].

Iteration.

While (stopping criterion false)

i = arg max [absolute value of <[r.bar], [[[phi].bar].sub.j]>];

I = I [union] {i};

[r.bar] = [y.bar] - [[PSI](:;I)][[[PSI](:;I)].sup.[dagger]] [y.bar].

End while

[[??].bar] (I) = [[PSI][(:, I)].sup.[dagger]] y.

Output. Sparse representations [[??].bar] and the original signal [y.bar] = [PSI] [[??].bar].

Here [[*].sup.[dagger]] is pseudoinverse matrix transformation. Stopping criterion is that iteration times are equal to the known number K of incident sources (i.e., sparsity). In the TDMMIMO estimation model, we can calculate the DOA and [[omega].sub.d] according to I vector. From I vector we can determine which atom or atoms will be represented; then p and q may be ascertained. From p value DOA maybe estimated and [[omega].sub.d] may be estimated from q value successfully.

When TDM-MIMO estimation model has L cycles, that is, Y = [[y.bar](1), [y.bar] (2), ..., [y.bar] (L)], we can solve Z efficiently using joint sparse recovery algorithms [23]. The joint sparse recovery spectrum of the DOA and DFS can be calculated by

[P.sub.y]([theta]) = 1/L [L.summation over (l=1)] [parallel][[??].sub.[theta]](l)[[parallel].sup.2], [theta] = [[theta].sub.1], [[theta].sub.2], ..., [[theta].sub.PQ]. (8)

So if there are multiple snapshots of the received array signals, We will use joint recovery algorithms to estimate the spectrum as formula (8) [24, 25].

4. Signal Reconstructions Algorithm for MMV Problem

Only one snapshot (L = 1) is needed to estimate parameters in sparse representations theory, but in actual application more snapshots are sampled in order to improve SNR and estimation accuracy. In formula (8) computation burden can become larger and larger with the increase of snapshots and the estimation results can be influenced more easily by SNR of input signals. This problem is so called MMV problem. In order to reduce computational load and improve estimation precision, there are two classical kinds of methods to solve MMV problem. One is dependent on SVD algorithm such as [l.sub.1]-SVD [26, 27]; the other is dependent on M-FOCUSS approach. These two algorithms both take datum correlation of different snapshots into consideration.

4.1. Signal Reconstruction Based on SVD. We investigate SVD subspace decomposition method [28, 29] firstly. The following formula is SVD subspace decomposition for Y matrix:

Y = U[LAMBDA][V.sup.H] = [[U.sub.s] [U.sub.N]][LAMBDA][V.sup.H]. (9)

Size of matrices [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is signal subspace corresponding to the larger singular values (assuming r) in [LAMBDA] matrix. [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is noise subspace corresponding to the smaller singular values (assuming [N.sub.Virt] - r) in [LAMBDA] matrix. So we can define

[Y.sub.S] = [YVD.sub.r]. (10)

Here [D.sub.r] = [[[[LAMBDA].sub.rxr] [O.sub.rx(L-r)]].sup.H] is a matrix composed by the larger singular values. Also we define [Z.sub.S] = [ZVD.sub.r], [N.sub.S] = [NVD.sub.R]. So the formula (7) can be espressed as

[Y.sub.S] = [PSI] x [Z.sub.S] + [N.sub.S], (11)

Size of matrices [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] have the same sparsity level. By using SVD algorithm the dimension of observation matrix has been reduced from [N.sub.Virt] x L to [N.sub.Virt] x r. In general r [much less than to] L and r < [N.sub.Virt] is conditioned. Meanwhile we will achieve good performance of parameters estimation in lower SNR. [l.sub.1]-SVD algorithm is a successful DOA estimation method for sparse solution. But convex optimization problem such as a second-order cone (SOC) programming has to be utilized in [l.sub.1]-SVD algorithm, so its computational complexity will increase. Here, we will continue to use OMP greedy algorithm after SVD is done. We generalize the scheme of algorithm based on SVD decomposition with OMP.

Input. The observation matrix Y and a sparse representations matrix [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Consider the following:

(1) SVD decomposition for matrix Y;

(2) determining the number of the larger singular values (r);

(3) calculating [Y.sub.S] = [YVD.sub.r];

(4) using OMP algorithm for each column of matrix [Y.sub.S];

(5) searching the index atoms in atomic dictionary.

Output. DOA and DFS.

4.2. Signal Reconstruction Based onM-FOCUSS. M-FOCUSS algorithm is developed based on the focal underdetermined system solver (FOCUSS) algorithm [30] developed for the single measurement case and it has an ability to utilize multiple measurement vectors to accurately identify the sparsity structure [21]. According to formula (7), we can summarize the scheme as follows so as to compute Z sparse solution.

Input. The observation matrix Y and a measurement matrix [PSI] = {[[bar.[phi]].sub.i], i = 1, 2, ..., m} [member of] [C.sup.nxm]

Initialization. [X.sub.1] = [1.sup.mxL], [epsilon], p, iternum = 0.

Iteration.

While ([parallel][X.sub.k+1] - [X.sub.k][[parallel].sub.2]/[parallel][X.sub.k][[parallel].sub.2] > [epsilon], iternum < max iteration)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

End while.

Output. Sparse representations matrix [??] = [X.sub.k+1].

Here, [parallel]x[[parallel].sub.2] = [([SIGMA][absolute value of [[x.sub.i].sup.2]).sup.1/2] is Frobenius norm and [epsilon] is threshold to stop the whole algorithm. The parameter p lies in [0, 2]. According to experiences, the best selection for p lies in [0.8, 1] so that convergence rate is good and optimum solution can be achieved easily. The optimal solution is conditioned by the minimum Frobenius norm with weighted values, so matrix [W.sub.k] is the kth iterative weighted matrix.

By using M-FOCUSS algorithm we can achieve sparse representations matrix [??] and deal with MMV problem successfully. Furthermore, we can use formula (8) [P.sub.y]([theta]) = (1/L) [[SIGMA].sup.L.sub.l=1] [parallel][[??].sub.[theta]](l)[[parallel].sup.2], where [[??].sub.[theta]](l) is the row vector of matrix [??] to calculate the joint sparse recovery spectrum function. Then the positions of peaks of spectrum function are the DFS and DOA estimation values and number of peaks is the number of input signals with different DFSs and DOAs.

5. Simulation Experiments

We present some simulations to verify the theoretical results below.We consider a TDM-MIMOradar with 4 RX and 2 TX antennas, uniformly spaced with an antenna distance of [lambda]/2; that is, [[d.bar].sup.RX] = [[1 2 3 4].sup.T] and [[d.bar].sup.TX] = [[0 5].sup.T]. The antennas are laid as shown in Figure 4.TheTXantennas transmit signal in TDM schedule and transmit schedule is shown in Figure 5.

We choose 1 cycle with [N.sub.Pulse] = 4 pulses and set the transmitting time instants to [t.bar] = [0, 1, 2, 3] in units of ms. So [[d.bar].sup.Pulse] = [[d.sup.TX.sub.1], [d.sup.TX.sub.2], [d.sup.TX.sub.2], [d.sup.TX.sub.1]] will be achieved. The basic simulation parameters are set as follows. The carrier frequency is equal to 77GHz. DOA search range lies in [-60[degrees], 60[degrees]] with step size 0.1[degrees] and velocity search range lies in [0 m/s, 100m/s] with step size 1 m/s. Snapshots are equal to L = 200. We do two experiments: one is for measurement of the root mean square errors (RMSE) of velocity and DOA estimation in different SNRs; the other is for multisource resolution of the coherent and uncorrelated signals.

5.1. Simulation Experiments of Joint DFS and DOA Estimation RMSE. 500 independent Monte Carlo simulations are carried out for each SNR varying from 0 dB to 30 dB with 5 dB step. The true DOA of incident signal and velocity of targetmoving are selected randomly in each trial. But DOA and velocity parameters set in one trial are constant. Velocity estimation is shown in Figure 6 and DOA estimation is shown in Figure 7. As we know, the formula [[omega].sub.d] = 2V/[lambda], where V is moving velocity of the target. So velocity simulation results are given directly instead of DFS due to the consideration of more concerned velocity parameter in the practical application. JS-RS stands for algorithm of the joint sparse recovery spectrum of DFS and DOA according to formula (8). SVD-RS stands for algorithm of joint estimation based on sparse representations using SVD with OMP algorithm to reconstruction; see formula (11) and scheme. MFOCUSSSR stands for improved M-FOCUSS algorithm with joint sparse recovery spectrum. Root mean square error (RMSE) is defined by

RMSE = [(1/J [J.summation over (j=1)] ([[??].sub.j] - [[theta].sub.j]).sup.1/2], (12)

where J is the times of independent Monte Carlo trials, [[??].sub.j] stands for parameter estimation, and [[theta].sub.j] is the true DOA or velocity for the jth trial.

Following [1], we compute CRLB [J.sup.-1] for the unknown Parameter [[THETA].bar] = [[DOA, DFS].sup.T]. The CRLB is a lower bound for the covariance matrix of any unbiased estimator [[??].bar.]. We use CRLB formula for DOA estimation ([CRB.sub.DOA]) according to paper [1].Thefollowing formulas are the Cramer-Rao bounds of [[??].bar]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (13)

[CRB.sub.DOA]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (14)

Here [Var.sup.S](*) stands for sample variance and [Cov.sup.S](*, *) stands for sample covariance. [[sigma].sup.2.sub.s]/[[sigma].sup.2] denotes the overall SNR. If an optimal TDM scheme is designed [1], [CRB.sub.DOA] will be minimized and [absolute value of [Cov.sup.S] ([[d.bar].sup.Pulse], [t.bar])] will also be minimized to zero.

From the simulation results, we know that the error of angle measurements can reach below 0.2[degrees] when SNR > 5 dB is satisfied for all three algorithms. But the error of velocity is not ideal for the JS-SR method under the same conditions. Only when SNR > 20 dB is satisfied, the error of velocity will reach V = 2m/s below. But the SVDRS and MFOCUSS-RS algorithms can achieve significant performance improvement. DOA estimation results of these two algorithms are closer to the corresponding Cramer-Rao lower bound (CRLB) compared to JS-RS algorithm. So both MFOCUSS-RS and SVD-SR are superior to JS-SR algorithm.

5.2. Simulation Experiments of Multisignals Resolution. Here, we have utilized these three algorithms to test coherent and uncorrelated signals resolution. Experiments have the same input conditions that are SNR = 20 dB and the true DOA of three simultaneous incident signals is [-25[degrees] 0[degrees] 25[degrees]].

Correlated coefficient is defined by

[[rho].sub.ik] = E[[s.sub.i](t)[s.sup.H.sub.k](t)]/[square root of E[[[absolute value of [s.sub.i](t)].sup.2]] E[[[absolute value of [s.sub.i](t)].sup.2]]]. (15)

Here E[*] denotes statistical expectation operator. i, k denote two input signals. If [[rho].sub.ik] which lies in [0, 1] is equal to 1, two input signals will be called as coherent signals.

The signal resolution simulation results are presented in Figure 8 to Figure 13. In Figures 8, 9, 10, 11, 12, and 13, E_DOA stands for estimation DOA and R DOA stands for true DOA. All three algorithms can process uncorrelated signals. But there are some spurious peaks for resolution results of JS-SR approach and the magnitude of spurious peaks may be larger than the magnitude of true input signal in some cases like Figure 8. So it is very complex to seek the different peak values of true DOAs. The improved M-FOCUSS algorithm can work well. The magnitude of spurious peaks is smaller than the magnitude of true DOAs like in Figure 12. With the increase of SNR, the magnitude of spurious peaks can also become smaller and smaller. In the situation of coherent signals impinging, all three algorithms have an ability to complete multisource resolution. But the problem of spurious peaks still exists for JS-SR algorithm. Accuracy of DOA estimation for SVD-SR algorithm has deteriorated since the rank of matrix [LAMBDA] is reduced and then we have to select one vector corresponding to the largest singular value to calculate the sparse solution.

Above all, the improved M-FOCUSS algorithm can work more robust than SVD-SR and JS-SR algorithms. Multisource resolution is successful and DFS and DOA estimation values are very precise regardless of correlated or uncorrelated input signals when the improved M-FOCUSS algorithm is utilized to deal with TDM-MIMO model.

6. Conclusions

In this paper, we propose joint DFS and DOA estimation algorithm for colocated TDM-MIMO radar based signal sparse representations. Facing the problem of time varying steering vector, we establish input signal mode and depict the algorithms of sparse representations and signal reconstruction in detail. Three algorithms including joint sparse recovery spectrum, SVD decomposition with OMP algorithm, and M-FOCUSS are discussed to deal with MMV problem. The joint sparse recovery spectrum and M-FOCUSS methods are combined to improve M-FOCUSS performance in the case of MMV problem. Finally, the simulation results of multisource resolution and accuracy of estimations and some discussions are demonstrated. Using sparse representations based on the proposed angle-frequency space, the coherent signals can be distinguished and the Doppler frequency shift and DOA can be estimated jointly. Particularly, the joint estimation algorithm based on sparse representations using improved M-FOCUSS reconstruction can deal with MMV problem and achieve significant performance improvement of DOA and DFS estimation accuracy. So the proposed algorithm and model are effective and can be further applied for practical application.

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

Conflict of Interests

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

Acknowledgments

This research paper is sponsored by China Scholarship Fund (CSC), in part supported by the Fundamental Research Funds for the Central Universities (HEUCF1408) and by the National Science Foundation of China under Grant no. 61201410.

References

[1] K. Rambach and B. Yang, "Colocated MIMO radar: cramerRao bound and optimal time division multiplexing for DOA estimation of moving targets," in Proceedings of the 38th IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '13), pp. 4006-4010, Vancouver, Canada, May 2013.

[2] S. Lutz and T. Walter, "Lens based 77 GHz TDM MIMO radar sensor for angular estimation in multitarget environments," in Proceedings of the 10th European Radar Conference (EuRAD '13), pp. 212-215, October 2013.

[3] M. Rossi, A. M. Haimovich, and Y. C. Eldar, "Spatial compressive sensing for MIMO radar," IEEE Transactions on Signal Processing, vol. 62, no. 2, pp. 419-430, 2014.

[4] 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.

[5] J. Liu, C. Han, and Y. Hu, "A novel coherent column replacement method in compressed sensing for DOA estimation," in Proceedings of the 32nd Chinese Control Conference (CCC '13), pp. 4832-4837, IEEE, July 2013.

[6] J. Yin and T. Chen, "Direction-of-arrival estimation using a sparse representation of array covariance vectors," IEEE Transactions on Signal Processing, vol. 59, no. 9, pp. 4489-4493, 2011.

[7] J.-F. Gu, W.-P. Zhu, and M. N. S. Swamy, "Compressed sensing for DOA estimation with fewer receivers than sensors," in Proceedings of the IEEE International Symposium of Circuits and Systems (ISCAS '11), pp. 1752-1755, IEEE, May 2011.

[8] M.-S. Lee, V. Katkovnik, and Y.-H. Kim, "System modeling and signal processing for a switch antenna array radar," IEEE Transactions on Signal Processing, vol. 52, no. 6, pp. 1513-1523, 2004.

[9] S. Yin, G. Wang, and X. Yang, "Robust PLS approach for KPIrelated prediction and diagnosis against outliers and missing data," International Journal of Systems Science, vol. 45, no. 7, pp. 1375-1382, 2014.

[10] A. Zwanetski and H. Rohling, "Continuous wave MIMO radar based on time division multiplexing," in Proceedings of the 13th International RadarSymposium (IRS '12), pp. 119-121, IEEE, May 2012.

[11] L. Liu, Y. Jiang, L. Wan, and Z. Tian, "Beamforming of joint polarization-space matched filtering for conformal array," The Scientific World Journal, vol. 2013, Article ID 589675,10 pages, 2013.

[12] L.-T. Wan, L.-T. Liu, W.-J. Si, and Z.-X. Tian, "Joint estimation of 2D-DOA and frequency based on space-time matrix and conformal array," The Scientific World Journal, vol. 2013, Article ID 463828,10 pages, 2013.

[13] H. Feng, L. Liu, and B. Wen, "2D-DOA estimation for cylindrical array with mutual coupling," Mathematical Problems in Engineering, vol. 2014, Article ID 716978, 8 pages, 2014.

[14] S. Yin, G. Wang, and H. R. Karimi, "Data-driven design of robust fault detection system for wind turbines," Mechatronics, vol. 24, no. 4, pp. 298-306, 2014.

[15] M. Jahn, R. Feger, C. Wagner, Z. Tong, and A. Stelzer, "A four-channel 94-GHz sige-based digital beamforming FMCW radar," IEEE Transactions on Microwave Theory and Techniques, vol. 60, no. 3, pp. 861-869, 2012.

[16] P. Schmalenberg, J. S. Lee, and K. Shiozaki, "A SiGe-based 16channel phased array radar system at W-Band for automotive applications," in Proceedings of the 10th European Radar Conference (EuRAD '13), pp. 299-302, October 2013.

[17] S. Yin, X. Yang, and H. R. Karimi, "Data-driven adaptive observer for fault diagnosis," Mathematical Problems in Engineering, vol. 2012, Article ID 832836, 21 pages, 2012.

[18] P. Wenig, M. Schoor, O. Gunther, B. Yang, and R. Weigel, "System design of a 77 GHz automotive radar sensor with superresolution DOA estimation," in Proceedings of the IEEE International Symposium on Signals, Systems and Electronics (ISSSE '07), pp. 537-540, August 2007.

[19] K. Rambach and B. Yang, "Direction of Arrival estimation of two moving targets using a time division multiplexed colocated MIMO radar," in Proceedings of the IEEE Radar Conference, pp. 1118-1123, Cincinnati, Ohio, USA, May 2014.

[20] T. Pedersen, C. Pedersen, X. Yin et al., "Joint estimation of doppler frequency and directions in channel sounding using switched Tx and Rx arrays," in Proceedings of the IEEE Global Telecommunications Conference (GLOBECOM '04), vol. 4, pp. 2354-2360, December 2004.

[21] B. D. Rao, K. Engan, and S. Cotter, "Diversity measure minimization based method for computing sparse solutions to linear inverse problems with multiple measurement vectors," in Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '04), vol. 2, pp. II369II372, IEEE, May 2004.

[22] Y. Wang, T. Fu, M. Gao, and S. Ding, "DOA estimation by exploiting spatial and Doppler sparsity," in Proceedings of the IETInternational Radar Conference, pp. 1-4, Xi'an, China, April 2013.

[23] B. Li, Y. Zou, and Y. Zhu, "Direction estimation under compressive sensing framework: a review and experimental results," in Proceedings of the International Conference on Information and Automation (ICIA '11), pp. 63-68, June 2011.

[24] Y. He, K. Hueske, E. Coersmeier, and J. Gdtze, "Efficient computation of joint direction-of-arrival and frequency estimation," in Proceedings of the 8th IEEE International Symposium on Signal Processing and Information Technology (ISSPIT '08), pp. 144149, IEEE, December 2008.

[25] Y. Wang, G. Leus, and A. Pandharipande, "Direction estimation using compressive sampling array processing," in Proceedings of the IEEE/SP 15th Workshop on Statistical Signal Processing (SSP '09), pp. 626-629, September 2009.

[26] A. Gretsistas and D. M. Plumbley, "An alternating descent algorithm for the off-grid DOA estimation problem with sparsity constraints," in Proceedings of the 20th European Signal Processing Conference (EUSIPCO '12), pp. 874-878, August 2012.

[27] D. Malioutov, M. Cetin, and A. S. Willsky, "A sparse signal reconstruction perspective for source localization with sensor arrays," IEEE Transactions on Signal Processing, vol. 53, no. 8, pp. 3010-3022, 2005.

[28] S. Yin, X. Li, H. Gao, and O. Kaynak, "Data-based techniques focused on modern industry: an overview," IEEE Transactions on Industrial Electronics, 2014.

[29] S. Yin, S. X. Ding, X. Xie, and H. Luo, "A review on basic datadriven approaches for industrial process monitoring," IEEE Transactions on Industrial Electronics, vol. 61, no. 11, pp. 64186428, 2014.

[30] I. F. Gorodnitsky and B. D. Rao, "Sparse signal reconstruction from limited data using FOCUSS: a re-weighted minimum norm algorithm," IEEE Transactions on Signal Processing, vol. 45, no. 3, pp. 600-616, 1997.

Tao Chen, HuanxinWu, and Lutao Liu

College of Information and Communication Engineering, Harbin Engineering University, Harbin 150001, China

Correspondence should be addressed to Lutao Liu; liulutao@msn.com

Received 26 July 2014; Accepted 29 September 2014; Published 16 November 2014

Academic Editor: Zhiguang Feng

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article; direction of arrival; time division multiplexing and multiple-input multiple-output |
---|---|

Author: | Chen, Tao; Wu, Huanxin; Liu, Lutao |

Publication: | Journal of Applied Mathematics |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 5350 |

Previous Article: | Reliability assessment of active distribution system using Monte Carlo simulation method. |

Next Article: | Video object tracking in neural axons with fluorescence microscopy images. |

Topics: |