Joint Direction-of-Departure and Direction-of-Arrival Estimation in a UWB MIMO Radar Detecting Targets with Fluctuating Radar Cross Sections

This paper presents a joint direction-of-departure (DOD) and direction-of-arrival (DOA) estimation in a multiple-input multipleoutput (MIMO) radar utilizing ultra wideband (UWB) signals in detecting targets with fluctuating radar cross sections (RCS). The UWBMIMO radar utilized a combination of two-way MUSIC and majority decision based on angle histograms of estimated DODs and DOAs at each frequency of the UWB signal.The proposed angle estimation scheme was demonstrated to be effective in detecting targets with fluctuating RCS, compared to conventional spectra averaging method used in subband angle estimations. It was found that a wider bandwidth resulted in improved estimation performance. Numerical simulations along with experimental evaluations in a radio anechoic chamber are presented.


Introduction
The introduction of multiple-input multiple-output (MIMO) radar enables numerous improvements on the conventional single input single output radar system.The MIMO radar is typically defined as a radar system utilizing multiple transmitting and receiving antennas that are either widely distributed or colocated [1][2][3][4].The former exploits independent path of transmit-receive pairs to increase the probability of detection [1,2].On the other hand, the latter uses orthogonality of transmitting signals to obtain an increased degree of freedom and parameter identifiability [3,4], which are very useful in target localization applications.
Target localization in radar has been intensively studied in literatures since the early years of radar.In general, radar systems estimate the target position by means of trilateration or triangulation.Trilateration can be implemented by using a minimum of two stations; however, the localization resolution is limited by the signal bandwidth, and usage of multiple stations is required to avoid ambiguities (ghost targets).On the contrary, triangulation is based on the angles of targets observed from the radar stations, and hence it does not suffer from the bandwidth constraint.In MIMO radar, it is possible to jointly estimate the direction-of-departure (DOD) and direction-of-arrival (DOA) by implementing array processing at both of the transmitting and receiving arrays, as depicted in Figure 1.This makes it suitable for triangulation-based localization.
Numerous works on DOD and DOA estimation have been reported [5][6][7][8].Nevertheless, these studies were mainly based on narrowband signal assumption.Target localization using the narrowband signal, however, was unstable due to fluctuation of target's radar cross section (RCS).In order to counter this problem, MIMO radars using multiple subcarriers and orthogonal frequency division multiplexing (OFDM) waveforms were proposed [9,10].Those schemes, however, only discussed one-dimensional angle estimation of a single target and implementation of limited number of subbands.Utilization of wider signal bandwidth might be useful in localizing target with significantly small or severely fluctuating RCS.One of the most promising technologies with wideband capability is ultra wideband (UWB) systems.The usage of UWB signal for angle estimation in MIMO radar, however, has not been given much attention, limiting the study to medical imaging and through-the-wall radars [11,12].Application of conventional angle estimation methods (e.g., Capon, multiple signal classification (MUSIC), and estimation of signal parameters via rotational invariance (ESPRIT)) to a UWB signal is a challenge since those methods inherently assume narrowband signals.
The present authors proposed a joint DOD and DOA estimation in a UWB MIMO radar using the combination of a two-way MUSIC and angle histograms [13].The basic idea in the proposed scheme was to treat the UWB signal as a summation of sinusoidal waves swept over the frequency band, and angle estimation was done at each of the frequency by means of two-way MUSIC [14].The estimation results were then combined using majority decisions formulated using angle histograms [15].However, the report only investigated estimation performance when detecting targets with a constant RCS, which was unlikely in actual environments.
In this paper, we evaluate the performance of the proposed algorithm in detecting targets with fluctuating RCS.Here, it will be shown that combining the estimation at different frequencies through majority decisions will overcome the problem of poor estimation when detecting fluctuating targets.A comparison between the proposed technique and the conventional spectrum averaging is also presented.Numerical simulations and experimental results are presented.
It is important to mention here that the RCS fluctuation problem has been continuously studied in the radar community.Until recently, special attention has been given to the subject in the case of MIMO radar, since the usage of MIMO configuration offers further degrees of freedom in the forms of spatial, frequency, and also waveform diversity.For example, the works in [16][17][18] employed spatial and waveform diversity in MIMO radars to increase the probability of detection and direction finding performance when detecting fluctuating target.The present study employs different approach from those reports since we focused on the utilization of frequency diversity.Studies regarding MIMO radars utilizing frequency diversity in detecting fluctuating RCS were studied, for instance, in [19,20], where several subbands with substantially wide frequency spacing were used for angle estimation.This study differs from those works from a point of view that we employ the diversity among a large number of subbands throughout the frequency bandwidth of a UWB system to enhance target localization performance and specifically demonstrated the application to MIMO radar.
The remainder of this paper is organized as follows.The next section discusses the proposed algorithm.Section 3 describes the numerical simulations, and Section 4 explained the results of experimental evaluations.Finally, the concluding remarks are presented in Section 5.

Proposed Scheme
2.1.Proposed Joint DOD and DOA Estimation Scheme.Consider a MIMO radar with M transmitting and N receiving elements, illuminating L uncorrelated targets located at the far field of transmit and receive arrays.At the transmit side, M orthogonal UWB signals are emitted, and each consists of multiple sinusoidal waves swept over the UWB bandwidth.Here, we define the complex transmit signal waveform vector by s() = [s 1 (), . . ., s  ()], where each term contains K frequency components.The orthogonality between the transmit signals can be achieved through time division scheme, where each transmitting antenna emits the UWB signal in separate time slots.Figure 2 illustrates an example of the transmitting signal which uses 3.1 to 10.6 GHz of sinusoidal waves in 1 MHz intervals.The receiving signal of the kth frequency component can be expressed by   where ⊗ denotes the Kronecker product, a t and a r are the transmit and receive steering vectors, respectively,  and  are the corresponding transmit and receive angles, s (k) (t) represents the kth frequency component of the transmit signal vector s(t), and n(t)is the total additive white Gaussian noise.The receiving signal covariance matrix of the kth frequency component is given by where E[⋅] is the ensemble average and [⋅] H represents the conjugate transpose operation.Here, singular value decomposition (SVD) of the covariance matrix gives where V (k) is a diagonal matrix whose diagonal elements contain the signal and noise eigenvalues for the kth frequency and E (k) is the corresponding eigenvectors of the signal and noise components.The two-dimensional spatial MUSIC spectrum at the kth frequency component can be constructed using the function where e N e H N is the noise eigenvectors obtained from the eigendecomposition of the receive signal covariance matrix in (2).Here, we have the L largest peaks which correspond to the DOD and DOA of the targets at each frequency component.The wideband DOD and wideband DOA are decided by taking the majority of the estimated angles among all the estimates at each frequency.This is denoted by "majority decision" in this study.The majority decision is formulated by initially combining DOD and DOA estimates at the K frequencies into a vector and arranging them in the form of angle histograms [13], as depicted in Figure 3.The angle histogram can be viewed as a function of angle i from −90 ∘ to 90 ∘ at intervals of, for example, 0.5 ∘ .The normalized number of occurrence of the peak angle is then given by where p (i) is the number of occurrences of the angle i and z is the normalized coefficient given by The majority decision is obtained by searching the peak of the histogram.As a benchmark, the performance of the proposed scheme will be compared with the conventional spectrum averaging method [21] used in existing subband processing scheme for angle estimation.The spectrum averaging method adopts an approach which takes the average of the estimated MUSIC spectra of all K frequency components: Then the wideband DOD and DOA are estimated from the L largest peaks of  MU ave .We will demonstrate in Section 3 that the majority decision technique performs better than the spectrum averaging method.Figure 4 depicts the block diagrams of the proposed and conventional schemes.

MIMO Radar Array Configuration Used in this Study.
The main advantage of a MIMO radar system is that the degrees of freedom can be enhanced by using the concept of virtual array [22,23].When orthogonal signals were transmitted from different antennas, the back scatter returns of each orthogonal signal will carry independent phase information that forms a new virtual array steering vectors at the receiver.

Transmitting antennas
Receiving antennas The virtual array can be characterized by convolution of the transmitting and receiving antenna positions [22].Given a MIMO radar using M transmitting and N receiving antennas, it is possible to form a full virtual array with nonoverlapping MN elements, by optimizing the antenna positions.Larger numbers of M and N contribute to construction of longer virtual array, which means further enhancement in angle estimation performance.Thus, from a signal processing point of view, it is generally important to use a larger virtual array.However, in practical applications, it is often beneficial to limit the number of antennas for the sake of cost and space.In this study, we demonstrate the proposed algorithm using a 4 × 4 MIMO array and discuss the performance within that limitation.
A full MN virtual array can be constructed, for example, by using transmitting antennas with spacing of Nd and receiving antennas with  spacing [23], where  is the distance equal to the half wavelength.
However, as shown in Figure 5, utilizing a two-way MUSIC using this array configuration resulted in spurious peaks in the MUSIC spectrum, since the receiving array's spacing is much larger than the distance of half wavelength.This problem can be overcome by limiting the scan range, but in the cost of narrower coverage area.
In this study, we employed a nonuniform array configuration as shown in Figure 6 to reduce the spurious peaks.This array configuration was chosen due to the minimum number of redundant elements that could be obtained in the virtual MIMO array.Refer to the Appendix for further explanation.The basic antenna spacing of  = 15 mm (corresponds to half wavelength of 10 GHz) was used.Figure 7 depicts the resulting MUSIC spectrum using the nonuniform array, where we could observe that the spurious peaks were eliminated, and the resulting spectrum peak sharpness is comparable to a full MN array with limited range scan.The authors have previously reported in [13] that this array configuration outperforms the conventional uniform linear array.

Complexity Analysis.
This subsection presents the analysis of the computational complexity of the proposed scheme.The computational burden of a conventional 2D-MUSIC has   been reported in literatures, such as in [5].For the sake of clarity, we broke down the complexity analysis of the proposed scheme in terms of SVD operation and searching algorithm computational costs, expressed using the O notation [24].The conventional spectrum averaging method [21] was also analyzed for comparison.
Considering the dimensions of the covariance matrix ×V and the total K frequency components used, applying SVD operation on the receiving covariance matrix generally costs ({ 2 V + V 3 }).Here, from the receiving covariance matrix formulation,  and V are given by M 2 and N 2 , respectively.This is the same for both the proposed and the spectrum averaging schemes.In terms of peak search operation, the proposed scheme performs two-dimensional peak search on the MUSIC spectrum, which costs ( 2 ), where i is the number of angle bins during search operations.In addition, the majority decision routine costs (2{ + }) + (2{}), where the first term corresponds to the histogram formulation of K total estimates, and the second term represents the one-dimensional peak search to identify the angles of L targets.On the other hand, the computational burden of the spectrum averaging method is () + ( 2 ), corresponding to the spectrum averaging and two-dimensional peak search of the averaged spectrum, respectively.
The computational complexity against MN antennas is plotted in Figure 8.It was shown that the proposed scheme marked larger computational burden than the spectral averaging method, particularly when the number of antennas is small.However, the difference of complexity is reduced with increasing number of antennas and converges when MN approaches 81.This is because the cost of SVD operation becomes dominant with larger dimension of receiving covariance matrix, thus resulting in similar order of total computational complexity in both methods.

Simulation Results
The proposed algorithm was simulated according to the parameters listed in Table 1.As mentioned in the previous sections, the UWB signal considered in the proposed scheme contains multiple sinusoidal ranging from 3.1 to 10.6 GHz at 1 MHz intervals.Here, selection of the start and stop frequencies within the frequency range determines the total bandwidth of the signal used.The signal bandwidth is defined by where   and   are the highest and lowest frequency components, respectively.The simulations will be conducted while varying the signal bandwidth from 10 to 1000 MHz, at different center frequencies,   , for example, from 3.6, 5, and 8.4 to 10 GHz.Two different scenarios were simulated, where the MIMO radar was detecting either a fixed point target or a target with fluctuating RCS.The fixed point target was modeled with a constant RCS, normalized by the free space propagation loss coefficient in each frequency.The targets with fluctuating RCS were modeled by Weibull distribution, since it was shown in literatures that measured RCS of complex targets such as automobiles and small cars follows Weibull distribution [25].The measurement data also showed that at a given observation angle, the RCS against frequency (7 GHz bandwidth) also follows similar distribution.The targets considered in the simulation were modeled to have Weibull distributed RCS throughout the UWB frequency range, whose cumulative distribution function is given by where a and b are the shape and scale parameters, respectively.The values of a and b were selected so as to equal the RCS medians of the Weibull and fixed point targets.Figure 9 plots the RCS versus frequency (expressed in signal power), and Figure 10 depicts the corresponding cumulative distribution.We demonstrate the performance of the proposed scheme when detecting 2 targets located at ( 1 ,  1 ) = (10 ∘ , −28 ∘ ) and ( 2 ,  2 ) = (−10 ∘ , −10 ∘ ).The signal to noise ratio (SNR) was 15 dB, and the number of snapshots was 50.
Figure 11 shows the examples of angle histograms at   = 3.6 GHz, when the MIMO radar was detecting the fixed point targets.It can be observed that the angle histograms were more or less similar when using either 1000 MHz or 50 MHz of signal bandwidths.Both cases marked sharp histogram peaks, which corresponds to the respective targets.This deduced that, when detecting fixed point targets, the bandwidth has minimal effects on the angle estimation performance.The estimation using a 50 MHz signal bandwidth, however, yielded an error of 1 ∘ .On the other hand, the angle estimation of the Weibull targets marked outstanding outliers in the angle histograms, particularly when using a 50 MHz signal bandwidth, as depicted in Figure 12(a).These outliers were attributable to spurious peaks in the MUSIC spectrum due to small RCS values at certain frequencies.However, taking the majority decisions from the histograms, the targets were successfully estimated within 2 ∘ of estimation error.It was also shown in Figure 12(b) that, by taking a larger bandwidth, for example, 1 GHz, the peaks of the outliers were suppressed compared to using a 50 MHz signal bandwidth.The performance of the proposed scheme was evaluated in terms of estimation error, defined by where W is the number of independent trials and  est and  est are the estimated DOD and DOA, respectively.The estimation error against signal bandwidth using   = 3.6 GHz was plotted in Figure 13, comparing the performance when detecting the fixed point and Weibull targets with 50 independent trials and SNR = 15 dB.Comparison with a conventional spectra averaging method [21] was also presented.It can be observed that the proposed scheme yielded improved estimation performance when taking larger signal bandwidth.This was true when detecting both types International Journal of Antennas and Propagation  of targets; however, the improvement was more significant in the case of Weibull targets.On the other hand, the estimation error increased with larger bandwidth when conventional spectra averaging method was applied.This is attributed to large number of outliers which significantly influenced the averaging result.This result demonstrated the effectiveness of the proposed method in detecting targets with severely fluctuating RCS.
The performance of the proposed scheme in terms of root mean square error (RMSE) against SNR is plotted in Figure 14.Performance of the spectrum averaging method and another conventional method by de Silva and Seow [14] was also presented for benchmarking purposes.Simulation was carried out assuming that all methods use a 4 × 4 MIMO nonuniform array detecting two Weibull targets.The performance of the proposed and spectrum averaging schemes was shown when they are using a 1000 MHz bandwidth with   = 3.6 GHz, while Silva's method was presented at several frequencies since it is a narrowband-based angle estimation.The performance bound calculated from the Cramér-Rao bound (CRB) when the DOD and DOA of a stationary target are estimated individually [26,27] was also plotted for reference.
We could observe that the proposed scheme marked the best performance, which was the nearest to the CRB, and produced no estimation error when the SNR exceeds 14 dB.It was also shown that the performance of Silva method depended on the frequency for the targets with fluctuating RCS.
The impact of a number of antennas on the performance of the proposed scheme was plotted in Figure 15.The RMSE was computed in 50 independent trials against MN, when using a signal bandwidth of 1000 MHz and SNR = 8 dB.Specific cases of M =  = 3, 4, 5, and 6 were simulated.It is shown in the figure that increasing number of antennas yielded improvement in RMSE performance, even in the low SNR condition.

Experimental Results
Experiments were conducted to verify the results of the numerical simulations.The measurements were done in a radio anechoic chamber, using a measurement setup illustrated in Figure 16.The setup consists of a vector network analyzer (VNA) and GPIB-controlled scanners which were used to virtually construct the MIMO arrays at both transmitting and receiving sides.The VNA was used to generate UWB signals from 3.1 to 10.6 GHz in 1 MHz intervals, which were used as the transmitting signal.The VNA was adopted due to its good calibration functions and time gating capability to cancel out direct coupling effects between transmitting and receiving antennas.Wideband horn antennas with average gain of 12.5 dBi were used as both transmitting and receiving antennas.The SNR of the system is defined as the ratio of the receiving signal average power to the thermal noise.
The setup was used to localize two targets positioned at ( 1 ,  1 ) = (10 ∘ , −28 ∘ ) and ( 2 ,  2 ) = (−10 ∘ , −10 ∘ ), similar to the condition in previous simulations.Two types of targets were used, which were conductive spheres and complex shaped targets, both fabricated using polystyrenes and aluminum foil.The conductive spheres have 20 cm of diameters which yield RCS of approximately −10 dBm 2 .The complex targets were constructed so that they yield a fluctuating RCS in the frequency domain.An example of the complex target is shown in Figure 17. Figure 18 shows the measured frequencydomain data of the complex targets, where  is the monostatic angle of observation.From the figure, we could observe that the receive signal fluctuates severely against frequency, compared to that of sphere targets.The complex targets   The estimation errors from measurement campaign were plotted in Figure 20.In the figure, the estimation errors were shown using signal bandwidths from 10 to 1000 MHz and   of 3.6, 5.0 and 8.4 GHz.As demonstrated in the numerical simulations, the measurement results also indicated a decreasing trend in estimation errors with increasing signal bandwidth, especially in the case of complex targets.The estimation errors of complex targets were worse than the sphere targets, where maximum of 7 ∘ of estimation errors was observed.Both cases marked poorer performance compared to the simulations since the experimental measurement included plane wave modeling errors and other maneuvering factors.Furthermore, no significant difference in estimation performance was observed when choosing different center frequencies.
A series of experiments was conducted to evaluate the localization performance of the proposed scheme.A single complex target was positioned in several locations in the radio anechoic chamber.The positions of the target are summarized in Table 2.The target was localized using the standard triangulation method, adopting the proposed scheme.In addition, localization using DOA and time of where Δ and Δ are the ranging errors along the  and  axes.The results of the measurements were plotted in Figure 21.In the case of triangulation, less than 1 m of error was obtained when using signal bandwidth of 1000 MHz.Using a 50 MHz signal resulted in lower accuracy.The DOA-TDOA method also produced good results when using a 1000 MHz signal and marked poor accuracy when using a 50 MHz bandwidth, due to low range resolution of TDOA.Both localization methods showed slightly lower accuracy when the target is further away from the radar due to decrease in SNR.

Concluding Remarks
The performance of a joint DOD and DOA estimation in a UWB MIMO radar detecting fluctuating targets was evaluated through numerical simulations and experimental evaluations.From the investigation, it was found that in detecting targets with fluctuating RCS (in this case Weibull distributed RCS against frequency), it is essential to use large signal bandwidth to reduce the estimation error using the proposed algorithm.When taking wider signal bandwidth, the usage of majority decisions from the angle histograms resulted in good estimation performance compared to the conventional spectrum averaging method.We concluded based on the work that the proposed scheme was a suitable candidate to implement joint angle estimation in MIMO radar using ultra wideband signal.utilizes multiple subcarriers, it is suitable to be extended to an OFDM-based radar system.

Appendix MIMO Radar Virtual Array Using Nonuniform Array Configuration
The MIMO antennas could be arranged in such a way that they produce a filled virtual array; however, the filled virtual array tends to consist of less number of unique virtual elements due to redundant elements.A few examples are depicted in Figure 22.Here, we illustrate the physical transmitting and receiving antennas position using the notation {1111}, where each entry corresponds to the number of antennas at the particular location on the /2 grid, and * is the convolutional operator.The results of the convolution of the transmitting and receiving antenna positions are shown on the right side of each figure.As depicted in Figure 22(a), the usage of uniform linear array on both transmitting and receiving arrays resulted in virtual arrays with significant number of redundant elements.

Figure 3 :
Figure 3: Formulation of angle histograms from two-way MUSIC spectra at each frequency component.
spectrum (dB) D O D (d eg .)D OA (d eg .)

Figure 7 :
Figure 7: Example of a MUSIC spectrum using the nonuniform array configuration (  = 5.0 GHz).

Figure 9 :Figure 10 :
Figure 9: Simulated RCS of a simulated Weibull target in comparison with a fixed point target considering propagation loss in frequency domain.

Figure 11 :
Figure 11: Angle histograms of estimated DODs and DOAs in detecting fixed point targets using a signal bandwidth of (a) 50 and (b) 1000 MHz.

Figure 12 :
Figure 12: Angle histograms of estimated DODs and DOAs in detecting Weibull targets using a signal bandwidth of (a) 50 and (b) 1000 MHz.

Figure 17 :
Figure 17: Example of a fabricated complex target.

Figure 18 :Figure 19 :
Figure 18: Measured frequency-domain data of the fabricated targets.

)Figure 20 :
Figure 20: Estimation errors against bandwidth obtained from measurements: (a) spherical and (b) complex targets.

Table 2 :
Target positions in a radio anechoic chamber.