Research Article  Open Access
LowComplexityBased RDMUSIC with Extrapolation for Joint TOA and DOA at Automotive FMCW Radar Systems
Abstract
Lowcomplexitybased reduceddimension–multiplesignal classification (RDMUSIC) is proposed with extrapolation for joint time delay of arrivals (TOA) and direction of arrivals (DOA) at automotive frequencymodulated continuouswave (FMCW) radar systems. When a vehicle is driving on the road, the automotive FMCW radar can estimate the position of multiple other vehicles, because it can estimate multiple parameters, such as TOA and DOA. Over time, the requirement of the accuracy and resolution parameters of automotive FMCW radar is increasing. To accurately estimate the parameters of multiple vehicles, such as range and angle, it is difficult to use a lowresolution algorithm, such as the twodimensional fast Fourier transform. To improve parameter estimation performance, highresolution algorithms, such as the 2DMUSIC, are required. However, the conventional highresolution methods have a high complexity and, thus, are not applicable to a realtime radar system for a vehicle. Therefore, in this work, a lowcomplexity RDMUSIC with extrapolation algorithm is proposed to have a resolution similar to that of a highresolution algorithm to estimate the position of other vehicles. Compared with conventional low complexity high resolution, in experimental results, the proposed method had better performance.
1. Introduction
Frequencymodulated continuouswave (FMCW) radar systems have many advantages, including lower cost and complexity, over equivalent pulse radar systems [1–3]. For FMCW radar, spatial–temporal parameters, such as multipath time delays of arrivals (TOA) and directions of arrivals (DOA), have been widely studied [4]. These two parameters are useful to estimate the position of moving targets in FMCW radar systems. Especially, a beamforming technique based on phased arrays for vehicle radar sensor systems [5] has been used for smart cruise control, traffic monitoring, and collision avoidance [2, 6, 7]. The FMCW radar has the characteristic that it decreases with the bandwidth corresponding to the frequency of the dechirped received signal irrespective of the transmitted bandwidth. Therefore, the signalprocessing complexity is significantly lower than that of conventional ultrawideband radar. By means of a dechirping method of an FMCW radio frequency (RF) module, the received signals can be transformed into sinusoidal waveforms to acquire TOAs and DOAs information. We can define these sinusoidal signals as beat signals.
The resolution parameters of the FMCW radar are estimated through a variety of algorithms, from conventional fast Fourier transform (FFT) to multiple signal classification (MUSIC) algorithms, as vehicle requirements increase. A onedimensional parameter estimator cannot be used to estimate multiple parameters simultaneously. For this reason, it is necessary to consider a twodimensional parameter estimator. In the case of a twodimensional parameter estimator, conventional 2DFFT has performance degradation for automotive radars requiring highresolution parameters, such as range and angle. To improve the estimation performance of range and angle parameters jointly, conventional twodimensional highresolution algorithms such as twodimensional estimation of signal parameters via a rotational invariant technique (2DESPRIT) and 2DMUSIC were used to estimate the parameters of multiple targets. Many studies represent the estimation accuracy of the superresolution algorithm is more two or three times than the conventional FFT algorithm. However, due to the high complexity of the highresolution algorithms, realtime automotive radar has difficulties in applying conventional algorithms. For this reason, various lowcomplexitybased 2D highresolution algorithms have been proposed.
In [8], this paper proposed a lowcomplexity superresolution joint angle and delay estimation algorithm for rangeazimuth FMCW radar. The phase shifts in time and arrays are exploited by the temporal and spatialtemporal smoothing technique of the proposed method. The proposed method is designed to estimate the range and angle sequentially without singular value decomposition (SVD) and eigenvalue decomposition (EVD) to reduce the computational burden. In case of complexity, the SVD or EVD to obtain eigenvector and eigenvalue has many computational burdens. In order not to use this SVD or EVD, this paper uses the squareinverse operator to separate signal and noise eigenvalue. However, this method degrades parameter estimation performance. In [9], in order to obtain a highresolution range profile (HRRP), the highresolution algorithms are required. This paper also needs low complexity highresolution algorithm because of applying the realtime system. So, the RELAX algorithm is used in this paper. However, since this paper focuses on the onedimensional spectrum, it is not suitable for the twodimensional parameter estimator required in this paper. Another study [10] presented collocated multipleinput multipleoutput (MIMO) radar that utilizes a lowcomplexity ESPRITbased DOA estimator. To reduce complexity, the fact that the Kronecker product results of transmit and receive steering vectors can be transformed to other steering vectors was used. However, this result has problems in that the DOA estimator is only used for MIMO radar. In [11], a crosscorrelation function of a received signal that can find a noise subspace without singular value decomposition (SVD) or eigenvalue decomposition (EVD) was proposed. However, this method can only be used for the twocomponent uniform linear array (ULA) of an Lshaped array. Therefore, conventional lowcomplexity 2D highresolution algorithms still are difficult to apply to automotive radar systems, and there is the limitation of only operating in certain conditions, such as Lshaped ones. In another study [4], a lowcomplexity algorithm using the combination of DFT and MUSIC was found to have the disadvantage that it is difficult to apply to the trend of increasing range resolution in automotive environments. Therefore, a lowcomplexitybased dimension RDMUSIC is proposed with an extrapolation algorithm for joint TOAs and DOAs in automotive FMCW radar systems. In Section 2, the system model is explained. In Section 3, the structure of the conventional 2DMUSIC estimator is analyzed. In Section 4, the lowcomplexitybased dimension RDMUSIC with the extrapolation algorithm for automotive FMCW is proposed. In Section 5, a simulation conducted to assess the estimation results of the proposed algorithm is described. And, the CPU execution time measured in MATLAB is used to evaluate the complexity reduction of the proposed algorithm. In Section 6, the estimation performance of the proposed method assessed through real experiments is discussed. Section 7 concludes this work.
1.1. Notation
(.)^{T} and (.)^{H} denote transpose and conjugate transpose, respectively. Operator is the Kronecker product, and is the identity matrix.
2. System Model
This section addresses the system model of the FMCW radar systems. Figure 1 shows the TX signal of FMCW radar in terms of change of frequency according to time. As shown in Figure 1, TX frequency linearly increases as the symbol period of FMCW radar, . Let us denote the initial frequency as , and the frequency at a time is denoted by and expressed as [2]: where is the bandwidth of the system. The chirp rate according to time is denoted by , i.e., . By employing frequency according to time, the instantaneous phase for is expressed into integral term as follows:
Consequently, the FMCW radar signal is expressed as:
targets over for the th target for a ULA consisting of K elements are considered and distance between adjacent arrays is set to , where is the carrier frequency’s wavelength as shown in Figure 2. One can represent the received signal in a linear and timeinvariant environment at the th antenna element with the complex amplitude in previous work [2] such that where is the complex amplitude of the th target, is the TOA, defines the DOA, and means the th antenna position’s additive white Gaussian noise (AWGN).
(a) Top view of ULA
(b) Structure of ULA
In the receiver (RX) part by dechirping, the RX signal from targets can multiply with the conjugation of the TX signals such that where means the th antenna array’s output results of the dechirping method. The beat signal for the th chirp and the th array is obtained and expressed as the product of the TOA, Doppler, and DOA terms as follows:
From an analogtodigital (ADC) converter for Nyquist sampling, where the sampling frequency needed is , the digital received signal of can be represented such that where
3. Conventional 2DMUSIC Estimator
The MUSIC algorithm introduced in [12] is a popular superresolution parameter estimation algorithm. The MUSIC algorithm makes use of the signal and noise subspaces of the timeaverage covariance matrix of the received signal. The conventional 2DMUSIC estimation algorithm is based on a stacking received data vector such that where , and means the parameter of selection for . Using the stacking received data matrix, the covariance matrix can be estimated as:
Using covariance matrix , we should obtain the signal and noise eigenvector to establish the superresolution spectrum. Through EVD, one can obtain the signal subspace matrix and noise subspace matrix by where the signal subspace matrix is the eigenvectors corresponding to the largest eigenvalues of , and the noise subspace matrix is the eigenvectors corresponding to the smallest eigenvalues of . The largest eigenvalues of , …, correspond to the eigenvectors of . The other eigenvalues correspond to the eigenvectors of such that . The MUSIC algorithm uses the characteristics that the steering vectors and are spanned by the signal subspace and the steering vectors are orthogonal to the columns of noise eigenvectors such that where and are defined by
Using the noise subspace matrix , the 2DMUSIC TOA and DOA spectrum function can be established as where and are the steering vectors of TOA and DOA, respectively. Here, the largest peaks of are taken as the estimates of the TOA and DOA for the automotive targets. For 2DMUSIC, because a huge 2D search and EVD are needed, it is inefficient because of the high computational cost. In the following sections, a lowcomplexitybased RDMUSIC with extrapolation is proposed for joint TOAs and DOAs at automotive FMCW radar systems.
4. Proposed LowComplexityBased RDMUSIC with Extrapolation for Joint TOAs and DOAs
The proposed estimation algorithm is based on the stacking received data matrix of Eq. (8). The signal subspace matrix and the noise subspace matrix can separate using the EVD of . The noise subspace matrix is orthogonal to the actual TOA and DOA steering vector at target. The denominator of Eq. (11) can be rewritten as
It is possible to transform the two dimensions to one dimension through a dimension reduction method, which can separate the TOA and the DOA parameters, by using Eq. (15). The TOAs are first determined by the extrapolation results of the received signal, as in a previous study [12]. In the first antenna element, the indices of the received signal’s extrapolation estimate the TOAs for multiple targets. To achieve extrapolation results, the FFT method is needed as where is the steering vector of FFT for the TOAs, i.e., . In practice, however, the finite data length causes smearing and leakage in the FFT.
Therefore, it is necessary to process the next step for the extrapolation algorithm using autoregression (AR). For extrapolation, each chirp signal’s AR parameters can be estimated, and linear prediction with the estimated parameters can be achieved. To estimate AR parameters, the modelbased techniques depend on modeling received signal with in the first antenna as the output of a linear system of a rational system form as follows:
This mathematical model to represent the given received signal by the polezero linear model in (17) is called an autoregressivemoving average (ARMA) model. This model is categorized as the ARMA model, AR model, and movingaverage (MA) model by the values of and . Among these three linear models, the AR model with is the most widely used approach to model nonwhite random process, because the AR model results in very simple linear equations for the AR parameters [13].
Using the AR model in Figure 3, the received signal in the first antenna element can be modeled as
(a)
(b)
(c)
(d)
In this case, the power spectrum of a th order AR process is defined as follows: where and can be estimated from the data and the accuracy of spectrum estimation, will depend on how accurately the model parameters may be estimated, and is the modeling error, which is assumed to be a random noise. Since AR spectrum estimation requires that an allpole model be found for the process, a variety of techniques may be used to estimate the allpole parameter [14]. The typical method to estimate an AR parameter for estimating distance information is the covariance method [14]. To find the th order AR parameters , the covariance method requires a set of linear equations,
The autocorrelation sequence is defined as
For distance spectrum estimation using the AR model, the order of the AR process should be determined by the number of targets . When the order is smaller than the number of targets , the results of the distance spectrum will be smoothed and will have low resolution. When the order is larger than the number of targets , the results of the distance spectrum will have spurious peaks. Another method of determining the model order involves the use of the Akaike Information Criterion (AIC) and Minimum Description Length (MDL) [15, 16]. The optimal modeling order can be selected by changing the model order until the values of the MDL or the AIC are minimized.
For the range estimation, a 1DFFT [2] is performed on the extrapolated chirp signals with the chirp index obtained. The 1DFFT results of the extrapolated received signal in the first array are expressed as where , and denotes the number of extrapolated samples. After 1DFFT is established, the peaks among the distance spectrum are detected, and the peak index means the distance results of targets.
After the extrapolation in the first antenna element is accomplished among the multiple targets, the TOA index vector is obtained, , where denotes . Therefore, the estimated TOA direction vector for the th target can be written as:
For each target DOA, by substituting the aboveestimated vector in Eq. (15), a normal MUSIC estimator at the th target is obtained such that where After the MUSIC estimation of the th target and peak detection processing, through the DOA index matrix of the multiple targets, the DOA of th target is accomplished by where is the inverse operator of the sine function.
Thus, the proposed algorithm TOA and DOA information for vehicle FMCW radar is outlined in Figure 4. The major steps of the proposed algorithm are as follows.
Step 1. Construct the spatial–temporal covariance matrix .
Step 2. Through the EVD of , obtain the signal subspace matrix and the noise subspace matrix , respectively.
Step 3. Through Eq. (15), separate the steering vectors of TOAs and the DOAs, respectively, using reduced dimension.
Step 4. Search for through the extrapolation and peak detection of extrapolation results, so one can obtain the estimated TOA terms of multiple targets, i.e., for .
Step 5. Substitute the estimated TOA direction vector of Eq. (23) in Eq. (15).
Step 6. Search for through the MUSIC estimator of Eq. (24) and peak detection, one can find the estimated DOA terms of multiple targets, i.e., for .
5. Simulation Results and Complexity Analysis
A MATLAB simulation was performed to verify the estimation results of the proposed algorithm compared with the conventional 2DMUSIC. In this section, the proposed algorithm and 2DMUSIC are only analyzed, because other algorithms are not based spectrum. A snapshot of the time delay–angle map is represented in the case of one target and two targets. In this section, the parameters are set as shown in Table 1. It was assumed that a 24GHz RF module of an FMCW radar with one transmitting channel and two receiving channels is used to obtain the simulation results.

For the simulation, as presented in Figure 5, a single target was placed at and . As in Figure 5, the proposed algorithm estimates similar angle and time delay results compared with the conventional algorithm.
(a)
(b)
Figure 6 shows a simulation snapshot of two targets in an AWGN channel. Two different targets were placed at and . Because this result is similar to that in Figure 5, the proposed algorithm and the conventional algorithm obtain similar angles and time delay results correctly.
(a)
(b)
In root mean squared error (RMSE), the proposed algorithm is analyzed for the RMSE of the estimated TOA and DOA compared to other algorithms. We compare the RMSE of the ideal RELAX and ESPRIT algorithm with RMSE of the proposed estimator in TOA and DOA. In Figures 7(a) and 7(b), two targets were placed at and . Two different targets are located at a greater distance interval than the distance resolution. Because the distance interval is greater than distance resolution, the RMSE results of all algorithm have similar results. In Figures 7(c) and 7(d), two targets were placed at and . Two different targets are located at lower distance interval than the distance resolution. Since the distance interval is lower than distance resolution, the RMSE results of MUSIC, RELAX, and the proposed algorithm have similar results while that of FFT has low performance.
(a)
(b)
(c)
(d)
For complexity, the proposed algorithm, 2DMUSIC, and 2DFFT are analyzed. The computational complexity of the algorithms is consisted of the primary multiplication operation. The 2DMUSIC and 2DFFT use full search to detect target for the ranges and angles, whereas the proposed algorithm requires the memoryefficient search to estimate multiple parameters. In case of 2DFFT method, times point FFT for distance estimation and times point FFT for angle estimation are accomplished. The required number of multiplication operation of 2DFFT is described by (28). And, the algorithm part that does not require SVD is only the FFT. The 2DMUSIC algorithm’s complexity is composed of autocorrelation matrix, EVD, and spectrum generation by the orthogonality in (27). In case of the proposed algorithm, 1Dextrapolation and RDMUSIC for range and angle are performed in (26). Here, denotes the number of spectrum samples of the MUSIC algorithm. When the complexity of the proposed structure in (26) is compared with that of the 2DMUSIC in (27), the proposed structure has a complexity of the , while the 2DMUSIC is composed of .
The complexity of the proposed algorithm is verified by measuring the CPU execution time in MATLAB. Figure 8 shows the measured CPU time based on the number of time samples . With , the complexity of the proposed method is at least 50 times lower than that of the conventional 2DMUSIC algorithm. The proposed method dramatically reduces the complexity burden by applying a reduced dimension method compared with the conventional 2DMUSIC algorithm while providing similar performance to the 2D MUSIC algorithm. 2DFFT do not analyze because 2DFFT’s complexity is much lower than other algorithms.
6. Experiments
In various experiments, the estimation results of the proposed method were evaluated in an anechoic chamber located at Daegu Gyeongbuk Institute of Science and Technology (DGIST) in South Korea. Using one TX channel and two received (RX) channels, an FMCW RF module was made with a 24GHz carrier frequency. The transmitter involved an oscillator with 26 MHz, a frequency synthesizer, and a voltagecontrolled oscillator. From the frequency synthesizer, an FMCW TX signal was produced for the 200MHz bandwidth in the range of 24.05–24.25 GHz. The receiver consisted of two highpass filters, two lowpass filters, two lownoise amplifiers (LNAs), and two mixers. The noise figures of the receiver and the LNA were 8 and 14 dB, respectively. The gain was set to 2.5 dB. By the mixer, the RF signal was changed to an intermediatefrequency signal as the beat signal. The developed 24GHz FMCW RF module is shown in Figure 9 [2].
To avoid an undesired reflected signal, the experiment was conducted in an anechoic chamber. For the experiments, a single target was placed at and . The experimental results of the timeangle map were estimated, as given in Figure 10. As shown in Figure 10, the proposed method and the conventional method estimated similar angles and time delays correctly, but the proposed algorithm had a lower complexity load than 2DMUSIC.
(a)
(b)
7. Conclusions
A lowcomplexitybased RDMUSIC with extrapolation for joint TOAs and DOAs at automotive FMCW radar systems was proposed. The proposed method considerably reduces the complexity of using RDMUSIC to reduce the dimensions from two dimensions to one dimension for the automotive FMCW radar system. The proposed method solves the lowcomplexity problem by using the extrapolation method instead of 1DMUSIC in the range estimation.
The simulation results showed that 2DMUSIC and the proposed algorithm have similar estimation performance, while the complexity of the proposed method is at least 50 times lower than that of the conventional 2DMUSIC algorithm with . Experimental analyses showed that the proposed algorithm provided a similar performance compared with 2DMUSIC. Therefore, the proposed method is applicable to automotive radar because of its low complexity for parameter estimation.
Data Availability
Access to data is restricted: the data that support the findings of this study are not available on request from our institution. The data are not publicly available due to restrictions, e.g., their containing information to allow for the commercialization of research findings.
Conflicts of Interest
The authors declare no conflict of interest.
Authors’ Contributions
Sangdong Kim and Bongseok Kim are cofirst authors.
Acknowledgments
This work was supported by the DGIST R&D Program of the Ministry of Science, ICT and Future Planning, Korea (20IT02 and 1902HRHR03).
References
 M. Dudek, I. Nasr, G. Bozsik, M. Hamouda, D. Kissinger, and G. Fischer, “System analysis of a phasedarray radar applying adaptive beamcontrol for future automotive safety applications,” IEEE Transactions on Vehicular Technology, vol. 64, no. 1, pp. 34–47, 2015. View at: Publisher Site  Google Scholar
 S. Kim, D. Oh, and J. Lee, “Joint DFTESPRIT estimation for TOA and DOA in vehicle FMCW radars,” IEEE Antennas and Wireless Propagation Letters, vol. 14, pp. 1710–1713, 2015. View at: Publisher Site  Google Scholar
 S. Saponara and B. Neri, “Radar sensor signal acquisition and multidimensional FFT processing for surveillance applications in transport systems,” IEEE Transactions on Instrumentation and Measurement, vol. 66, no. 4, pp. 604–615, 2017. View at: Publisher Site  Google Scholar
 S. Kim, Y. Ju, and J. Lee, “A lowcomplexity joint TOAs and AOAs parameter estimator using dimension reduction for FMCW radar systems,” Elektronika ir Elektrotechnika, vol. 24, no. 4, 2018. View at: Publisher Site  Google Scholar
 C. Schroeder and H. Rohling, “Xband FMCW radar system with variable chirp duration,” in 2010 IEEE Radar Conference, pp. 1255–1259, Washington, DC, USA, 2010. View at: Publisher Site  Google Scholar
 A. Townley, P. Swirhun, D. Titz et al., “A 94GHz 4TX4RX phasedarray FMCW radar transceiver with antennainpackage,” IEEE Journal of SolidState Circuits, vol. 52, no. 5, pp. 1245–1259, 2017. View at: Publisher Site  Google Scholar
 S. Lee, Y. J. Yoon, J. E. Lee, and S. C. Kim, “Humanvehicle classification using featurebased SVM in 77GHz automotive FMCW radar,” IET Radar, Sonar & Navigation, vol. 11, no. 10, pp. 1589–1596, 2017. View at: Publisher Site  Google Scholar
 D. Oh and J. H. Lee, “Lowcomplexity rangeazimuth FMCW radar sensor using joint angle and delay estimation without SVD and EVD,” IEEE Sensors Journal, vol. 15, no. 9, pp. 4799–4811, 2015. View at: Publisher Site  Google Scholar
 S. J. Lee, S. J. Jeong, E. Yang, and K. T. Kim, “Target identification using bistatic highresolution range profiles,” IET Radar, Sonar & Navigation, vol. 11, no. 3, pp. 498–504, 2017. View at: Publisher Site  Google Scholar
 X. Zhang and D. Xu, “Lowcomplexity ESPRITbased DOA estimation for colocated MIMO radar using reduceddimension transformation,” Electronics Letters, vol. 47, no. 4, pp. 283284, 2011. View at: Publisher Site  Google Scholar
 N. Xi and L. Liping, “A computationally efficient subspace algorithm for 2D DOA estimation with Lshaped array,” IEEE Signal Processing Letters, vol. 21, no. 8, pp. 971–974, 2014. View at: Publisher Site  Google Scholar
 S. Kim and K. K. Lee, “LowComplexity Joint ExtrapolationMUSICBased 2D Parameter Estimator for Vital FMCW Radar,” IEEE Sensors Journal, vol. 19, no. 6, pp. 2205–2216, 2019. View at: Publisher Site  Google Scholar
 B. W. Choi, E. H. Bae, J. S. Kim, and K. K. Lee, “Improved prewhitening method for linear frequency modulation reverberation using dechirping transformation,” The Journal of the Acoustical Society of America, vol. 123, no. 3, pp. EL21–EL25, 2008. View at: Publisher Site  Google Scholar
 S. L. Marple, Digital Spectral Analysis with Applications, Prentice Hall, Englewood Cliffs, New Jersey, 1987.
 H. Akaike, “Fitting autoregressive models for prediction,” Annals of the Institute of Statistical Mathematics, vol. 21, no. 1, pp. 243–247, 1969. View at: Publisher Site  Google Scholar
 J. Rissanen, “Modeling by shortest data description,” Automatica, vol. 14, no. 5, pp. 465–471, 1978. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Sangdong Kim 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.