1. Introduction
In recent years, radar systems have been installed in automobiles to detect targets located in multiple directions. Typically, automotive radar systems use frequencies in the 24-GHz or 77–81-GHz band. Since such a high frequency band is used, the miniaturization of the radar antenna system has become possible. In the automotive radar system, the number of receiving antenna elements is gradually decreasing to reduce the manufacturing cost of the radar. Therefore, various techniques have been proposed for accurately estimating the direction of arrival (DOA) of the target with the limited antenna elements as much as possible.
Array interpolation is one of the methods for improving the DOA estimation accuracy using a limited number of antenna elements. Several studies on enhanced array interpolation methods have been conducted [
1,
2,
3,
4,
5]. In [
1], the Taylor series approximation was used to generate interpolated array elements in a uniform circular array and achieved improved DOA estimation performance; however, because the order of the series is limited to one less than the maximum number of array elements, the approximation performance is not guaranteed for automotive radar systems that use only a few (e.g., four or eight) array elements. In addition, the norm-constrained least squares method was used to find the interpolated microphone array in [
2], but the problem-solving process was heuristic because the proper norm constraint parameter was determined empirically. Recently, the effective interpolation method that changes the nonuniform co-prime array to the uniform linear array was proposed in [
3]. Furthermore, the enhanced DOA estimation method using the characteristics of the covariance matrix generated by the virtual array interpolation for the co-prime array has been proposed in [
4]. In [
5], the authors improved the angle estimation performance by linearly predicting the received signals and virtually extending the array.
Among the various array interpolation methods, an array interpolation method that moves array elements from an original location to a desired location using a transformation matrix is widely used [
6,
7,
8]. To this end, the linear least squares (LLS) method has been widely used to identify the proper transformation matrix [
6,
7,
8,
9,
10,
11,
12]. However, a transformation matrix obtained by means of the LLS method is not the best solution for interpolating array elements. When this transformation matrix is applied, interpolated array elements are generated by linear combinations of original array elements. In this case, the amplitudes of interpolated array elements can be different from those of original array elements. If amplitude differences exist among the array elements, the performance of DOA estimation algorithms is degraded [
13]. In addition, because the solution derived from the LLS method is obtained in the process of simultaneously minimizing differences in amplitudes and phases, the phase information of the interpolated array elements is not accurately formulated, which is a critical factor for DOA estimation.
Thus, in this paper, we propose a transformation matrix in the logarithmic domain for the array interpolation. We focus on minimizing the phase differences between the original and the interpolated array elements. First, we take logarithms for the array elements and extract the phase information from them. We then apply the LLS method to the logarithmic-domain matrices to find an appropriate transformation matrix. Finally, the interpolated array elements are generated by the new matrix, and the DOA estimation is conducted. Based on a comparison of interpolation errors of the proposed and conventional transformation methods, our array transformation method successfully interpolates newly-produced array elements with more elaborate phases. In addition, the proposed array transformation does not affect the amplitudes of the interpolated array elements; they are conserved even after the transformation. Moreover, for a given antenna array, the proposed transformation matrix for that array is calculated and stored (offline) in advance. Thus, we do not have to calculate the transformation matrix in real time.
We also extend the proposed array interpolation scheme to received signal interpolation. When we use the transformation matrices obtained by the LLS method and our proposed method, the powers of the interpolated received signals are not uniform over all array elements. In this case, the effect of the array interpolation and the performance of the DOA estimation are not fully ensured. Thus, to mitigate this problem, we also propose a calibration method for the interpolated received signal powers. The simulation results confirm that our proposed method performs DOA estimation better than the conventional array interpolation method. In addition, based on actual measurement data acquired using an automotive radar, our method shows improved angular resolution and estimation performance. The proposed method can be effectively used in a radar system using a small number of array elements.
The remainder of this paper is organized as follows. In
Section 2, we introduce the basic signal model for the array antenna, as well as the conventional array interpolation technique using the LLS method. Next, the proposed array interpolation method is described in
Section 3. In this section, we also propose a method of calibrating the interpolated received signals for more accurate DOA estimation. Simulation and measurement results are provided in
Section 4 and
Section 5, respectively. We conclude this paper in
Section 6.
4. Simulation Results
Many studies have been conducted on the location in which to interpolate array elements to improve the accuracy of DOA estimation algorithms. In [
6,
8], the authors located the interpolated array elements that minimized interpolation errors within given conditions. In addition, the array searching method proposed in [
12] revealed enhanced DOA estimation accuracy with the interpolated array. However, this method was deemed too heuristic and time consuming. In this paper, to verify the DOA estimation accuracy resulting from our proposed interpolation method, we transformed the original array elements to the minimum-redundancy linear arrays, while maintaining identical apertures. In general, minimum-redundancy linear arrays show the maximum resolution for a given number of array elements by minimizing the number of redundant spacings in the array [
15,
16]. Moreover, previous studies have reported that non-uniform linear arrays perform better at DOA estimation than do uniform linear arrays that have the same apertures [
17,
18]. Therefore, in our simulation, by transforming the original array to the non-uniform minimum-redundancy linear array, we analyzed the performance improvement in the DOA estimation.
In the simulation, we used four array elements (
) that are widely used in automotive long-range radar (LRR). The location of the original array elements was
. It is well known that the minimum-redundancy linear array location of four array elements is
[
15,
16]. Thus, using the array transformation matrices, we interpolated array elements in the location
. Here, we assumed that two targets were located at
and adopted the Bartlett method [
14] as the DOA estimation algorithm. In addition, the signal-to-noise ratio (SNR) at the array elements was set to 10 dB, and 1000 time samples were used to construct the correlation matrix used in the Bartlett algorithm. The FOV was given as
, which was equivalent to the FOV of the LRR. Since
and
were calculated and stored only once when the number of array elements and the FOV were given, the stored values can be used repeatedly without having to identify another
and
.
First, under these simulation conditions, we calculated two types of interpolation errors, which were given as:
The smaller the error values were calculated based on (24), the more accurate the array interpolation was conducted. For both transformation matrices,
and
, we calculated the interpolation errors by changing the size of the FOV. The result is given in
Figure 1. As the figure shows, the interpolation errors calculated from
were almost close to zero. In addition, for the FOV of the LRR (i.e., the size of the FOV being
), the errors are given as
and
. Therefore, judging from both types of interpolation errors, our proposed array transformation matrix
was more approximate to
than was
. In other words, the interpolated array elements were accurately generated when the proposed interpolation method was employed. For larger FOV sizes, the interpolation errors of the conventional method became larger because the interpolation matrix was calculated more accurately when the DOA range of the targets was tightly within the FOV.
Using these transformation matrices, we formulated the received signals and conducted the DOA estimation. As shown in
Figure 2, with the original received signals, the Bartlett method could not resolve the two targets, and the estimated DOA was
. In general, when we used four array elements with
spacing, the half-power beamwidth became
. Therefore, the difficulty in distinguishing those given DOAs was reasonable. Even with the interpolated received signals from
, two different DOAs were not estimated, and the estimated DOA was
, which was not the exact value. However, with the interpolated received signals from
, the Bartlett method showed enhanced angular resolution, and we can find two different DOAs such as
. Moreover, when using the interpolated received signal vector with the power calibration,
, the best estimation result was achieved, and the estimated DOA values were
, which were close to the actual DOA values.
For the statistical performance evaluation, we calculated the resolution probability
for the conventional Bartlett algorithm and the Bartlett with array interpolation methods. This probability is defined as:
where
indicates the number of times that two distinct DOAs were extracted from the received signals and
denotes the number of simulations. Since we conducted this simulation 1000 times under the same conditions,
became 1000. In addition, we calculated the root mean squared error (RMSE) defined as:
where
is the estimated value of
in the
simulation. When the number of the estimated targets was one, we used this as
. The results are shown in
Table 1. Considering the resolution probability and the RMSE, our proposed method performed better than the conventional Bartlett and the Bartlett with the transformation matrix
. In addition, while maintaining the simulation conditions, except the array SNR values, we calculated the resolution probability and the RMSE. As
Figure 3 and
Figure 4 show, our proposed method yielded good estimation results despite the different array SNR values. Moreover, after changing the number of time samples used to build the correlation matrix, a performance comparison among the interpolation methods was conducted, and the results of which are given in
Figure 5 and
Figure 6. Even though only a few time samples were used, our proposed array transformation showed improved estimation performance.
We also conducted a simulation for a case in which three targets existed in the FOV of the radar. The simulation was conducted while maintaining the same simulation conditions given in
Figure 1, except for the target information, and the result is shown in
Figure 7. Here, the targets were located at
. The conventional Bartlett and the Bartlett with the transformation matrix
each estimated only two DOAs:
and
, respectively. Thus, these methods failed to resolve the targets placed at
. However, when applying our proposed transformation matrix, we could identify the three different DOAs. Moreover, from the power calibrated interpolated received signal vector
, the DOAs were estimated as
, which were the most exact estimated values. We also compared the performance of the proposed method to that of the multiple signal classification (MUSIC) algorithm, which is known as a high-resolution DOA estimation algorithm [
19]. To apply the MUSIC algorithm, the number of targets must be estimated in advance using the Akaike information criterion or minimum description length [
20,
21]. If the number of targets is well estimated (
), the most exact performance occurs. However, if the number is not accurately estimated (e.g.,
or
), the estimation performance deteriorates considerably, and it cannot be used as shown in
Figure 7. In addition, since the MUSIC algorithm performs the eigenvalue decomposition and the multiplication of matrices spanned by the noise eigenvectors, additional computational complexity
occurs compared to the conventional beamforming algorithm (i.e., the Bartlett method) [
22,
23]. Moreover, the Bartlett method is more robust to noise variance than the MUSIC algorithm [
14]. Thus, for automotive radars, the Bartlett algorithm may be more appropriate for stably estimating the DOA of a target under noisy road environments.
Under the same simulation conditions, we also applied the total least squares estimation of signal parameters via rotational invariance techniques (TLS EPSRIT) [
24]. The TLS ESPRIT method is one of subspace-based DOA estimation algorithms like the MUSIC and is a nonparametric DOA estimation method. When we used the TLS ESPRIT algorithm, the DOAs were estimated as
, which shows a large difference from the actual values. Because the TLS ESPRIT performs the eigenvalue decomposition three times to estimate the DOA, it requires more computation than our method. In addition, when using a small number of antenna elements, proper DOA estimation performance is not guaranteed with the ESPRIT method. Moreover, we compared the performance of the beamspace MUSIC algorithm [
25] with that of our proposed method. We needed to find the approximate DOA of the target in the beamspace MUSIC method, which was an unnecessary process in our proposed method. After finding the approximate angle, a beamforming matrix was generated based on that angle. Generating the appropriate beamforming matrix is the most important point of the beamspace DOA estimation algorithms. For example, if the beamforming matrix was generated in the range of
–
, the DOA was estimated as
, which was close to the real value. However, if the beamforming matrix was formed between
and
, the DOA was estimated as
. In addition, when using beamspace DOA estimation algorithm, the algorithm had to be repeated as many times as the number of targets.
Furthermore, simulations were conducted not only for the four array elements, but also for three and five array elements. When the number of array elements was three, the original location of the array elements was given as
. This array was transformed to the minimum-redundancy array, and interpolated array elements were located at
[
15,
16]. In addition, we assumed that targets were located at
and that the FOV ranged from
–
. Since the half-power beamwidth for the given array was
, the array had a very low angular resolution, and the given DOAs were difficult to distinguish from the conventional Bartlett algorithm. In addition, for the five array elements, the location of the original array elements was given as
, and it was transformed to the location
[
15,
16]. For this case, the FOV was equal to that of the LRR, and targets were placed at
. These DOAs were hard to separate out using the conventional Bartlett because the half-power beamwidth for the given array was
. For both cases of three and five array elements, the resolution probability and the RMSE were as given in
Figure 8,
Figure 9,
Figure 10 and
Figure 11, respectively, by increasing the array SNR from 0 dB–10 dB. As shown in the figures, our method also performed better for cases in which the number of array elements was three and five.
5. Measurement Results
To verify the performance of our proposed method, we also conducted actual measurements on a testing ground of the Mando Corporation using its automotive LRR. In the measurement, a single-element transmit antenna and four-element receiving uniform linear array antenna () were used, and the spacing between adjacent elements was 1.8. In addition, the half-power beamwidth of the array antenna was , and the FOV of the LRR ranged from –. This antenna system was equipped with an automotive radar and transmitted a 76.5-GHz frequency-modulated continuous wave signal. The transmitted signal was reflected from the front targets, and then, the reflected signals were received by the array antenna.
Figure 12 shows the measurement environment. Two identical target vehicles were located at
and were 40 m away from a radar-equipped vehicle. In this measurement, we also used the Bartlett algorithm for the DOA estimation method and calculated the resolution probability and the RMSE for the original received signals and the interpolated received signals derived from the array interpolation methods. Under the same measurement environment, we recorded 600 radar scans. Thus,
in (25) and (26) became 600 in this case. The results are listed in
Table 2. Similar to the simulation results, based on both measures, the DOA estimation with the proposed transformation matrix
showed better angular resolution and estimation accuracy than that of the conventional Bartlett and Bartlett method with the transformation matrix
. Furthermore, the estimation with
showed the most improved resolution and estimation performance.
Using the same automotive radar, the measurements were conducted on the expressway, as shown in
Figure 13. In the experimental data, 100 radar scans were extracted when two targets were almost at the same distance and were close each other. For those cases, the DOA estimation methods were applied, and the results are listed in
Table 3. Although the performance of the proposed algorithm was slightly lower than in the environment of
Figure 12, the proposed algorithm showed improved angular resolution and lower RMSE over the other algorithms. It can be seen that the overall angular estimation performance degraded from the actual experimental measurement results rather than the simulation results. This was because the quality of the received signal was degraded due to the clutter caused by the surrounding road structures in an actual road environment [
26,
27]. In the environment shown in
Figure 13, because the radar signals reflected from the wall of the tunnel were received with those from the desired targets, the quality of the received signal was inevitably deteriorated.