Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Periodic Nonlinear Error Analysis and Compensation of a Single-Excited Petal-Shaped Capacitive Encoder to Achieve High-Accuracy Measurement
Next Article in Special Issue
An Approach for the Customized High-Dimensional Segmentation of Remote Sensing Hyperspectral Images
Previous Article in Journal
Ultraprecision Real-Time Displacements Calculation Algorithm for the Grating Interferometer System
Previous Article in Special Issue
Eliminating the Effect of Image Border with Image Periodic Decomposition for Phase Correlation Based Remote Sensing Image Registration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Logarithmic-Domain Array Interpolation for Improved Direction of Arrival Estimation in Automotive Radars

1
Machine Learning Lab, AI & SW Research Center, Samsung Advanced Institute of Technology (SAIT), Gyeonggi-do 16678, Korea
2
Department of Electrical and Computer Engineering and Institute of New Media & Communications (INMC), Seoul National University (SNU), Seoul 08826, Korea
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(10), 2410; https://doi.org/10.3390/s19102410
Submission received: 12 April 2019 / Revised: 18 May 2019 / Accepted: 22 May 2019 / Published: 27 May 2019
(This article belongs to the Special Issue Computational Intelligence in Remote Sensing)

Abstract

:
In automotive radar systems, a limited number of antenna elements are used to estimate the angle of the target. Therefore, array interpolation techniques can be used for direction of arrival (DOA) estimation to achieve high angular resolution. In general, to generate interpolated array elements from original array elements, the method of linear least squares (LLS) is used. When the LLS method is used, the amplitudes of the interpolated array elements may not be equivalent to those of the original array elements. In addition, through the transformation matrix obtained from the LLS method, the phases of the interpolated array elements are not precisely generated. Therefore, we propose an array transformation matrix that generates accurate phases for interpolated array elements to improve DOA estimation performance, while maintaining constant amplitudes of the array elements. Moreover, to enhance the effect of our interpolation method, a power calibration method for interpolated received signals is also proposed. Through the simulation, we confirm that the array interpolation accuracy and DOA estimation performance of the proposed method are improved compared to those of the conventional method. Moreover, the performance and effectiveness of our proposed method are also verified using data obtained from the commercial radar system. Because the proposed method exhibits better performance when applied to actual measurement data, it can be utilized in commercial automotive radar systems.

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.

2. Fundamentals of Array Interpolation

2.1. Signal Model for the Array Antenna

We assume that signals, coming from K directions of θ 1 , θ 2 , , θ K , are incident on N linearly placed array elements. The location of the array elements is [ d 1 , d 2 , , d N ] and θ k ( k = 1 , 2 , , K ) is defined from the boresight direction of the array antenna. Assuming far-field narrow-band signal sources, the received signal vector of the array at time t can be expressed as:
X ( t ) = A × S ( t ) + N ( t ) = [ x 1 ( t ) , x 2 ( t ) , , x N ( t ) ] T
where [ · ] T denotes the vector transpose operator and A = [ a ( θ 1 ) , a ( θ 2 ) , , a ( θ K ) ] is the steering matrix composed of the steering vectors a ( θ k ) given by:
a ( θ k ) = [ e j 2 π λ d 1 sin θ k , e j 2 π λ d 2 sin θ k , , e j 2 π λ d N sin θ k ] T .
λ denotes the wavelength corresponding to the carrier frequency. S ( t ) = [ s 1 ( t ) , s 2 ( t ) , , s K ( t ) ] T is the incident signal vector, where s k ( t ) ( k = 1 , 2 , , K ) is the complex amplitude of the incident signal from the kth signal source at time t. These amplitudes are assumed to be zero-mean complex Gaussians and uncorrelated with each other. In addition, the noise vector N ( t ) = [ n 1 ( t ) , n 2 ( t ) , , n N ( t ) ] T is assumed to be the zero-mean complex Gaussian vector, and its components are also uncorrelated with each other. The samples from S ( t ) and N ( t ) are also assumed to be uncorrelated with each other. If there is a correlation between the signals received from the radar system, the correlation can be lowered through signal processing techniques such as spatial smoothing [14].

2.2. Conventional Array Interpolation Method

In this section, we briefly introduce the conventional array interpolation technique using the transformation matrix derived from the LLS method. For the field of view (FOV) of an automotive radar, which is expressed as:
Θ = { θ p θ p = θ L + ( p 1 ) × θ R θ L P 1 , p = 1 , 2 , , P } ,
we find a suitable matrix that transforms the original array elements into the interpolated array elements ( θ L and θ R are angles that indicate the left and the right boundaries of the FOV). In this range, the steering matrix of the original array elements can be given as:
A ( Θ ) = [ a ( θ L ) , a ( θ L + Δ θ ) , , a ( θ L + ( P 2 ) Δ θ ) , a ( θ R ) ]
where Δ θ = θ R θ L P 1 is the angle step size. Then, if we want to interpolate array elements in the location [ g 1 , g 2 , , g M ] , the steering matrix of the interpolated array elements is determined as:
B ( Θ ) = [ b ( θ L ) , b ( θ L + Δ θ ) , , b ( θ L + ( P 2 ) Δ θ ) , b ( θ R ) ]
where b ( θ p ) = [ e j 2 π λ g 1 sin θ p , e j 2 π λ g 2 sin θ p , , e j 2 π λ g M sin θ p ] T ( p = 1 , 2 , , P ) . Assuming that a matrix T transforms the original steering matrix to the interpolated steering matrix, it can be expressed as:
B ( Θ ) = T × A ( Θ ) .
To find the proper transformation matrix T , the least squares method is used as:
T * = arg min T ( B ( Θ ) T × A ( Θ ) F )
where · F denotes the Frobenius matrix norm. Then, based on the method of LLS, the transformation matrix can be determined as:
T * = B ( Θ ) A ( Θ ) H ( A ( Θ ) A ( Θ ) H ) 1
where [ · ] H denotes the Hermitian matrix operator. Finally, from (6) and (8), the estimate of B ( Θ ) is given as:
B ^ T * ( Θ ) = T * × A ( Θ ) = B ( Θ ) A ( Θ ) H ( A ( Θ ) A ( Θ ) H ) 1 × A ( Θ ) .
In addition, the ( m , p ) th element of the matrix B ^ T * ( Θ ) can be expressed as:
B ^ ( m , p ) T * ( Θ ) = n = 1 N T * ( m , n ) × A ( n , p ) ( Θ ) = n = 1 N T * ( m , n ) × e j 2 π λ d n sin θ p , m { 1 , 2 , , M } , p { 1 , 2 , , P }
where T * ( m , n ) and A ( n , p ) ( Θ ) denote the ( m , n ) th and the ( n , p ) th elements of the matrices T * and A ( Θ ) , respectively.
Through this transformation matrix, received signals for the interpolated array elements can also be generated, which are defined as:
Y ^ ( t ) = T * × X ( t ) = [ y ^ 1 ( t ) , y ^ 2 ( t ) , , y ^ M ( t ) ] .
By utilizing these interpolated received signals, the authors of [6,7,8,9,10,11,12] conducted improved DOA estimations.

3. Logarithmic-Domain Array Interpolation

3.1. Proposed Array Interpolation Method

When we use the conventional transformation matrix for array interpolation, a major problem occurs. Based on the transformation matrix obtained from the LLS method, interpolated array elements are generated by linear combinations of original array elements. In this case, the amplitudes of the interpolated array elements may not be equivalent to those of the original array elements. In other words, B ^ ( m , p ) T * ( Θ ) does not always become unity. When the amplitudes of each array element are not uniform over the entire array, DOA estimation performance is degraded [13]. In addition, based on the solution derived from the LLS method, the phases of the interpolated array elements are not precisely generated. For DOA estimation, the phase information of the interpolated array elements is critical. Therefore, in this section, we propose a more effective array interpolation method that minimizes the phase differences between the original and the interpolated array elements while maintaining the equivalent amplitudes of the array elements.
For the elements of the steering matrices on both sides of (6), we take logarithms such as:
L O G ( A ( Θ ) ) = log ( A ( 1 , 1 ) ( Θ ) ) log ( A ( 1 , P ) ( Θ ) ) log ( A ( 2 , 1 ) ( Θ ) ) log ( A ( 2 , P ) ( Θ ) ) log ( A ( N , 1 ) ( Θ ) ) log ( A ( N , P ) ( Θ ) ) , L O G ( B ( Θ ) ) = log ( B ( 1 , 1 ) ( Θ ) ) log ( B ( 1 , P ) ( Θ ) ) log ( B ( 2 , 1 ) ( Θ ) ) log ( B ( 2 , P ) ( Θ ) ) log ( B ( M , 1 ) ( Θ ) ) log ( B ( M , P ) ( Θ ) ) ,
where L O G ( · ) denotes the operator that takes logarithms for each element in the matrix and B ( m , p ) ( Θ ) indicates the ( m , p ) th element of the matrix B ( Θ ) . All elements in matrices L O G ( A ( Θ ) ) and L O G ( B ( Θ ) ) have pure imaginary values. Then, in the logarithmic domain, we find a proper transformation matrix V that transforms L O G ( A ( Θ ) ) to L O G ( B ( Θ ) ) , which is expressed as:
L O G ( B ( Θ ) ) = V × L O G ( A ( Θ ) ) .
As in the original domain, the appropriate matrix V can be found using the LLS method, and the solution is given as:
V * = L O G ( B ( Θ ) ) L O G ( A ( Θ ) ) H × L O G ( A ( Θ ) ) L O G ( A ( Θ ) ) H 1 .
This matrix V * effectively transforms the phases of the original array elements into those of the interpolated array elements. However, since the matrix is defined in the logarithmic domain, it cannot be directly applied to the original array elements as in (9). In other words, this transformation matrix cannot be expressed with a linear operator. Instead, it can be written with the original array elements as:
B ^ ( m , p ) V * ( Θ ) = n = 1 N A ( n , p ) ( Θ ) V * ( m , n ) , m 1 , 2 , , M , p 1 , 2 , , P
where B ^ ( m , p ) V * ( Θ ) is the ( m , p ) th element of the newly-interpolated steering matrix B ^ V * ( Θ ) and V * ( m , n ) indicates the ( m , n ) th element of the matrix V * .
The conventional transformation matrix in (8) formulates the interpolated array elements with linear combinations of the original array elements. However, this new transformation matrix generates only the phase information of the interpolated array elements using combinations of the phases of the original array elements. In other words, based on our transformation, it conserves the amplitudes of the original array elements in the interpolated array elements because:
B ^ ( m , p ) V * ( Θ ) = n = 1 N A ( n , p ) ( Θ ) V * ( m , n ) = e n = 1 N j 2 π λ d n sin θ p V * ( m , n ) = 1 , m 1 , 2 , , M , p 1 , 2 , , P .
Therefore, the proposed array transformation affects only the phases of the interpolated array elements and generates more accurate phases for the interpolated array elements. The transformation matrix T * does not preserve the amplitudes of the original array elements because B ^ ( m , p ) T * ( Θ ) is not always unity. Thus, the interpolation accuracy derived from the new transform matrix V * is higher than that from the conventional matrix T * .

3.2. Enhanced Received Signal Interpolation

Similar to the received signal interpolation in (11), received signals of the interpolated array elements with the transformation matrix V * are expressed as:
z ^ m ( t ) = n = 1 N x n ( t ) V * ( m , n ) Z ^ ( t ) , m { 1 , 2 , , M } .
Using Z ^ ( t ) = [ z ^ 1 ( t ) , z ^ 2 ( t ) , , z ^ M ( t ) ] , we conduct the DOA estimation and can achieve improved performance compared to the estimation using Y ^ ( t ) .
For a much better DOA estimation, we also consider the power of the received signals. When we use the interpolated received signal vectors, Y ^ ( t ) and Z ^ ( t ) , power differences exist among the interpolated received signals. In other words,
y ^ m ( t ) 2 = y ^ m ( t ) 2 and z ^ m ( t ) 2 = z ^ m ( t ) 2 ( for m m , m { 1 , 2 , , M } )
does not always hold, because:
n = 1 N T * ( m , n ) 2 n = 1 N T * ( m , n ) 2 , T * ( m , n ) T * ( m , n ) ¯ T * ( m , n ) T * ( m , n ) ¯ and n = 1 N V * ( m , n ) n = 1 N V * ( m , n ) ( for n n , n { 1 , 2 , , N } , for m m , m { 1 , 2 , , M } )
(see Appendix A) where · ¯ denotes the complex conjugate of a complex number. This power imbalance can cause performance degradation in the DOA estimation [13]. In our proposed method, the amplitudes of B ^ ( m , p ) V * ( Θ ) are equivalent for all array elements. However, it is not directly related to the powers of the interpolated received signals, and power differences exist among the interpolated received signals. Therefore, to mitigate this problem, we propose an effective compensation method to formulate the received signals of each interpolated array element such that they have similar power levels while maintaining the effect of our proposed phase interpolation method. In other words, the compensated received signal is given as:
w ^ m ( t ) = n N ( m ) x n ( t ) 1 N ( m ) × exp j n N ( m ) V * ( m , n ) x n ( t ) W ^ ( t )
where:
N ( m ) = { n * n * = arg n V * ( m , n ) 0 , n = 1 , 2 , , N } , m 1 , 2 , , M ,
and N ( m ) denotes the cardinality of the set N ( m ) . When comparing w ^ m ( t ) with z ^ m ( t ) , the interpolated phase of w ^ m ( t ) is the same as that of z ^ m ( t ) . Therefore, the phase interpolation effect from the transformation matrix V * is maintained. In addition, when we use this compensated interpolated received signal, the following equation is always established as:
w ^ m ( t ) 2 = w ^ m ( t ) 2 ( for m m , m { 1 , 2 , , M } ) ,
because:
w ^ m ( t ) 2 = n N ( m ) x n ( t ) 2 N ( m ) = x 1 ( t ) × × x N ( t ) N ( m ) 2 N ( m ) = x 1 ( t ) 2 × × x N ( t ) 2 N ( m ) 1 N ( m ) γ , m 1 , 2 , , M .
In other words, the powers of the interpolated received signals are nearly equivalent among the array elements. Thus, if we use the received signal vector W ^ ( t ) for the DOA estimation, we can achieve more enhanced performance than when using Y ^ ( t ) and Z ^ ( t ) .

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 ( N = 4 ) that are widely used in automotive long-range radar (LRR). The location of the original array elements was [ d 1 , d 2 , d 3 , d 4 ] = [ 0 , 2 λ , 4 λ , 6 λ ] . It is well known that the minimum-redundancy linear array location of four array elements is [ 0 , 1 , 4 , 6 ] [15,16]. Thus, using the array transformation matrices, we interpolated array elements in the location [ g 1 , g 2 , g 3 , g 4 ] = [ 0 , 1 λ , 4 λ , 6 λ ] . Here, we assumed that two targets were located at [ θ 1 , θ 2 ] = [ 3.5 , 2.5 ] 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 Θ = { θ p θ p = 10 + ( p 1 ) × 0.1 , p = 1 , 2 , , 201 } , which was equivalent to the FOV of the LRR. Since T * and V * 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 T * and V * .
First, under these simulation conditions, we calculated two types of interpolation errors, which were given as:
E U = B ( Θ ) B ^ U ( Θ ) F 2 , E p h a s e U = B ( Θ ) B ^ U ( Θ ) F 2 , U { T * , V * } .
The smaller the error values were calculated based on (24), the more accurate the array interpolation was conducted. For both transformation matrices, T * and V * , 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 V * were almost close to zero. In addition, for the FOV of the LRR (i.e., the size of the FOV being 20 ), the errors are given as [ E T * , E V * ] = [ 1.240 , 4.719 × 10 28 ] and [ E p h a s e T * , E p h a s e V * ] = [ 1.004 , 4.700 × 10 28 ] . Therefore, judging from both types of interpolation errors, our proposed array transformation matrix B ^ V * ( Θ ) was more approximate to B ( Θ ) than was B ^ T * ( Θ ) . 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 0.1 . In general, when we used four array elements with 2 λ spacing, the half-power beamwidth became 6.5 . Therefore, the difficulty in distinguishing those given DOAs was reasonable. Even with the interpolated received signals from T * , two different DOAs were not estimated, and the estimated DOA was 1.2 , which was not the exact value. However, with the interpolated received signals from V * , the Bartlett method showed enhanced angular resolution, and we can find two different DOAs such as [ 2.8 , 2.0 ] . Moreover, when using the interpolated received signal vector with the power calibration, W ^ ( t ) , the best estimation result was achieved, and the estimated DOA values were [ 3.1 , 2.2 ] , which were close to the actual DOA values.
For the statistical performance evaluation, we calculated the resolution probability P r for the conventional Bartlett algorithm and the Bartlett with array interpolation methods. This probability is defined as:
P r = N r N t × 100 ( % )
where N r indicates the number of times that two distinct DOAs were extracted from the received signals and N t denotes the number of simulations. Since we conducted this simulation 1000 times under the same conditions, N t became 1000. In addition, we calculated the root mean squared error (RMSE) defined as:
RMSE = k = 1 K q = 1 N t θ k θ ^ k ( q ) 2 N t ( )
where θ ^ k ( q ) is the estimated value of θ k ( k = 1 , 2 ) in the q t h ( q = 1 , 2 , N t ) simulation. When the number of the estimated targets was one, we used this as θ ^ k ( q ) . 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 T * . 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 [ θ 1 , θ 2 , θ 3 ] = [ 8 , 1.5 , 7.5 ] . The conventional Bartlett and the Bartlett with the transformation matrix T * each estimated only two DOAs: [ 8.4 , 3.2 ] and [ 8.6 , 2.7 ] , respectively. Thus, these methods failed to resolve the targets placed at [ θ 2 , θ 3 ] = [ 1.5 , 7.5 ] . However, when applying our proposed transformation matrix, we could identify the three different DOAs. Moreover, from the power calibrated interpolated received signal vector W ^ ( t ) , the DOAs were estimated as [ 8.6 , 1.8 , 8.8 ] , 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 ( K ^ = 3 ), the most exact performance occurs. However, if the number is not accurately estimated (e.g., K ^ = 1 or K ^ = 2 ), 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 O N 3 + N 2 × ( 2 N 2 L 1 ) 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 [ 26.5 , 13.8 , 12.7 ] , 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 1 3 , the DOA was estimated as 1.8 , which was close to the real value. However, if the beamforming matrix was formed between 1 and 5 , the DOA was estimated as 2.3 . 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 [ d 1 , d 2 , d 3 ] = [ 0 , 1.5 λ , 3 λ ] . This array was transformed to the minimum-redundancy array, and interpolated array elements were located at [ g 1 , g 2 , g 3 ] = [ 0 , 1 λ , 3 λ ] [15,16]. In addition, we assumed that targets were located at [ θ 1 , θ 2 ] = [ 4 , 6.5 ] and that the FOV ranged from 15 15 . Since the half-power beamwidth for the given array was 12 , 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 [ d 1 , d 2 , d 3 , d 4 , d 5 ] = [ 0 , 2.25 λ , 4.5 λ , 6.75 λ , 9 λ ] , and it was transformed to the location [ g 1 , g 2 , g 3 , g 4 , g 5 ] = [ 0 , 1 λ , 4 λ , 7 λ , 9 λ ] [15,16]. For this case, the FOV was equal to that of the LRR, and targets were placed at [ θ 1 , θ 2 ] = [ 1 , 3 ] . These DOAs were hard to separate out using the conventional Bartlett because the half-power beamwidth for the given array was 4.5 . 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 ( N = 4 ) were used, and the spacing between adjacent elements was 1.8 λ . In addition, the half-power beamwidth of the array antenna was 7 , and the FOV of the LRR ranged from 10 10 . 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 [ θ 1 , θ 2 ] = [ 1.7 , 4.6 ] 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, N t 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 V * showed better angular resolution and estimation accuracy than that of the conventional Bartlett and Bartlett method with the transformation matrix T * . Furthermore, the estimation with W ^ ( t ) 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.

6. Conclusions

In this paper, we proposed a logarithmic-domain transformation matrix used for array interpolation to improve the accuracy of DOA estimation. Our transformation matrix was obtained by minimizing the differences between the phases of the original array elements and the interpolated array elements. Our proposed method identified a more accurate transformation matrix with less phase distortion, and the amplitudes of the array elements were maintained after the transformation. In addition, to improve the accuracy of the DOA estimation algorithm, we proposed a method for adjusting the powers of the interpolated received signals to a similar level. Finally, from the simulation and the measurement results, we verified that our new method showed much better angular resolution and estimation accuracy than did the DOA estimation using the conventional transformation matrix derived from the LLS method. The proposed method can be effectively applied to radar systems using a small number of antenna elements, such as automotive radar systems.

Author Contributions

Conceptualization, S.L. and S.-C.K.; methodology, S.L.; software, S.L.; validation, S.L.; formal analysis, S.L.; investigation, S.L.; resources, S.-C.K.; data curation, S.L.; writing, original draft preparation, S.L.; writing, review and editing, S.-C.K.; visualization, S.L.; supervision, S.-C.K.; project administration, S.-C.K.; funding acquisition, S.-C.K.

Funding

This work was supported by the Technology Innovation Program (or the Industrial Strategic Technology Development Program, No. 10080086, Development of 77/79 GHz Dual Band Radar) funded By the Ministry of Trade, Industry & Energy (MOTIE, Korea).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DOADirection of arrival
LLSLinear least squares
FOVField of view
LRRLong-range radar
SNRSignal-to-noise ratio
RMSERoot mean squared error
MUSICMultiple signal classification
TLS ESPRITTotal least squares estimation of signal parameters via rotational invariance techniques

Appendix A

First, y ^ m ( t ) and y ^ m ( t ) in (11) can be expressed as:
y ^ m ( t ) = n = 1 N T * ( m , n ) × x n ( t ) , y ^ m ( t ) = n = 1 N T * ( m , n ) × x n ( t ) .
If we assume that received signals of the original array elements have nearly equal power such as:
x n ( t ) 2 γ 0 ( n = 1 , 2 , , N ) ,
then the powers of each signal are given as:
y ^ m ( t ) 2 = n = 1 N T * ( m , n ) × x n ( t ) 2 = γ × n = 1 N T * ( m , n ) 2 + 2 n = 1 , n = 1 , n n N T * ( m , n ) T * ( m , n ) ¯ x n ( t ) x n ( t ) ¯ , y ^ m ( t ) 2 = n = 1 N T * ( m , n ) × x n ( t ) 2 = γ × n = 1 N T * ( m , n ) 2 + 2 n = 1 , n = 1 , n n N T * ( m , n ) T * ( m , n ) ¯ x n ( t ) x n ( t ) ¯
where { · } is the operator that takes the real part of a complex number. Since:
n = 1 N T * ( m , n ) 2 n = 1 N T * ( m , n ) 2 , T * ( m , n ) T * ( m , n ) ¯ T * ( m , n ) T * ( m , n ) ¯ ( for n n , n { 1 , 2 , , N } , for m m , m { 1 , 2 , , M } ) ,
the powers of the interpolated received signals from the transformation matrix T * are not equivalent in all array elements.
In addition, the interpolated received signals in (17) can be expressed as:
z ^ m ( t ) = n = 1 N x n ( t ) V * ( m , n ) , z ^ m ( t ) = n = 1 N x n ( t ) V * ( m , n ) .
Therefore, the powers of each term are given as:
z ^ m ( t ) 2 = n = 1 N x n ( t ) V * ( m , n ) 2 = x 1 ( t ) V * ( m , 1 ) 2 × × x N ( t ) V * ( m , N ) 2 = n = 1 N x n ( t ) V * ( m , n ) 2 , z ^ m ( t ) 2 = n = 1 N x n ( t ) V * ( m , n ) 2 = x 1 ( t ) V * ( m , 1 ) 2 × × x N ( t ) V * ( m , N ) 2 = n = 1 N x n ( t ) V * ( m , n ) 2 .
Since V * ( m , n ) R and V * ( m , n ) R , the above equation is redefined as:
z ^ m ( t ) 2 = n = 1 N x n ( t ) 2 V * ( m , n ) = n = 1 N γ V * ( m , n ) = γ n = 1 N V * ( m , n ) , z ^ m ( t ) 2 = n = 1 N x n ( t ) 2 V * ( m , n ) = n = 1 N γ V * ( m , n ) = γ n = 1 N V * ( m , n ) .
However, because:
n = 1 N V * ( m , n ) n = 1 N V * ( m , n ) ( for m m , m { 1 , 2 , , M } ) ,
the interpolated received signals have different powers for each element.

References

  1. Cao, M.-Y.; Huang, L.; Xie, W.-X.; So, H.C. Interpolation array technique for direction finding via Taylor series fitting. In Proceedings of the IEEE China Summit and International Conference on Signal and Information Processing (ChinaSIP), Chengdu, China, 12–15 July 2015. [Google Scholar]
  2. Doblinger, G. Optimized design of interpolated array and sparse array wideband beamformers. In Proceedings of the 16th European Signal Processing Conference, Lausanne, Switzerland, 25–29 August 2008. [Google Scholar]
  3. Zhou, C.; Gu, Y.; Shi, Z.; Zhang, Y.D. Off-grid direction-of-arrival estimation using coprime array interpolation. IEEE Signal Process. Lett. 2018, 11, 1710–1714. [Google Scholar] [CrossRef]
  4. Zhou, C.; Gu, Y.; Fan, X.; Shi, Z.; Mao, G.; Zhang, Y.D. Direction-of-arrival estimation for coprime array via virtual array interpolation. IEEE Trans. Signal Process. 2018, 11, 5956–5971. [Google Scholar] [CrossRef]
  5. Sim, H.; Lee, S.; Kang, S.; Kim, S.-C. Enhanced DOA estimation using linearly predicted array expansion for automotive radar systems. IEEE Access 2019, 4, 47714–47727. [Google Scholar] [CrossRef]
  6. Friedlander, B. Direction finding using an interpolated array. In Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP), Munich, Germany, 21–24 April 1990. [Google Scholar]
  7. Weiss, A.J.; Gavish, M. Direction finding using ESPRIT with interpolated arrays. IEEE Trans. Signal Process. 1991, 6, 1473–1478. [Google Scholar] [CrossRef]
  8. Friedlander, B. The root-MUSIC algorithm for direction finding with interpolated arrays. ELSEVIER Signal Process. 1993, 1, 15–29. [Google Scholar] [CrossRef]
  9. Friedlander, B.; Weiss, A.J. Direction finding for wide-band signals using an interpolated array. IEEE Trans. Signal Process. 1993, 4, 1618–1634. [Google Scholar] [CrossRef]
  10. Tuncer, T.E.; Yasar, T.K.; Friedlander, B. Direction of arrival estimation for nonuniform linear arrays by using array interpolation. Radio Sci. 2008, 8, 1–11. [Google Scholar] [CrossRef]
  11. Kim, K.; Sarkar, T.K.; Palma, M.S. Adaptive processing using a single snapshot for a nonuniformly spaced array in the presence of mutual coupling and near-field scatterers. IEEE Trans. Antennas Propag. 2002, 5, 582–590. [Google Scholar]
  12. Kang, S.; Lee, S.; Lee, J.-E.; Kim, S.-C. Improving the performance of DOA Estimation using virtual antenna in automotive radar. IEICE Trans. Commun. 2017, 5, 771–778. [Google Scholar] [CrossRef]
  13. Lee, S.; Yoon, Y.-J.; Kang, S.; Lee, J.-E.; Kim, S.-C. Enhanced performance of MUSIC algorithm using spatial interpolation in automotive FMCW radar systems. IEICE Trans. Commun. 2018, 1, 163–175. [Google Scholar] [CrossRef]
  14. Krim, H.; Viberg, M. Two decades of array signal processing research: The parametric approach. IEEE Signal Process. Mag. 1996, 7, 67–94. [Google Scholar] [CrossRef]
  15. Moffet, A.T. Minimum-redundancy linear arrays. IEEE Trans. Antennas Propag. 1968, 3, 172–175. [Google Scholar] [CrossRef]
  16. Vertatschitsch, E.; Haykin, S. Nonredundant arrays. Proc. IEEE 1986, 1, 217. [Google Scholar] [CrossRef]
  17. Kassis, C.E.; Picheral, J.; Mokbel, C. Advantages of nonuniform arrays using root-MUSIC. ELSEVIER Signal Process. 2010, 2, 689–695. [Google Scholar] [CrossRef]
  18. Mhamdi, A.; Samet, A. Direction of arrival estimation for nonuniform linear antenna. In Proceedings of the IEEE International Conference on Communications, Computing and Control Applications (CCCA), Hammamet, Tunisia, 3–5 March 2011. [Google Scholar]
  19. Schmidt, R.O. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 3, 276–280. [Google Scholar] [CrossRef]
  20. Akaike, H. A new look at the statistical model identification. IEEE Trans. Automatic Control 1974, 12, 716–723. [Google Scholar] [CrossRef]
  21. Wax, M.; Kailath, T. Detection of signals by information theoretic criteria. IEEE Trans. Acoust. Speech Signal Process. 1985, 4, 387–392. [Google Scholar] [CrossRef]
  22. Pan, V.Y.; Chen, Z.Q. The complexity of the matrix eigenproblem. In Proceedings of the Proceedings of the 31st annual ACM Symposium on Theory Of Computing (STOC), Marina del Rey, CA, USA, 5–7 May 1999. [Google Scholar]
  23. Rubsamen, M.; Gershman, A.B. Direction-of-arrival estimation for nonuniform sensor arrays: From manifold separation to Fourier domain MUSIC methods. IEEE Trans. Signal Process. 2009, 2, 588–599. [Google Scholar] [CrossRef]
  24. Ottersten, B.; Viberg, M.; Kailath, T. Performance analysis of the total Least squares ESPRIT algorithm. IEEE Trans. Signal Process. 1991, 5, 1122–1135. [Google Scholar] [CrossRef]
  25. Yeom, D.-J.; Park, S.-H.; Kim, J.-R.; Lee, M.-J. Performance analysis of beamspace MUSIC with beamforming angle. In Proceedings of the 8th IEEE International Conference on Signal Processing and Communication Systems (ICSPCS), Gold Coast, Australia, 15–17 December 2014. [Google Scholar]
  26. Diewald, F.; Klappstein, J.; Sarholz, F.; Dickmann, J.; Dietmayer, K. Radar-interference-based bridge identification for collision avoidance systems. In Proceedings of the IEEE Intelligent Vehicles Symposium (IV), Baden-Baden, Germany, 5–9 June 2011. [Google Scholar]
  27. Lee, S.; Lee, B.-H.; Lee, J.-E.; Kim, S.-C. Statistical characteristic-based road structure recognition in automotive FMCW radar systems. IEEE Trans. Intell. Transp. Syst. (Early Access) 2018, 9, 1–12. [Google Scholar] [CrossRef]
Figure 1. Two types of interpolation errors from T * and V * .
Figure 1. Two types of interpolation errors from T * and V * .
Sensors 19 02410 g001
Figure 2. Normalized Bartlett pseudospectrums for two adjacent targets located at [ 3.5 , 2.5 ] .
Figure 2. Normalized Bartlett pseudospectrums for two adjacent targets located at [ 3.5 , 2.5 ] .
Sensors 19 02410 g002
Figure 3. Resolution probabilities versus SNR ( N = 4 ).
Figure 3. Resolution probabilities versus SNR ( N = 4 ).
Sensors 19 02410 g003
Figure 4. Root mean squared errors versus SNR ( N = 4 ).
Figure 4. Root mean squared errors versus SNR ( N = 4 ).
Sensors 19 02410 g004
Figure 5. Resolution probabilities versus the number of time samples (SNR is 10 dB).
Figure 5. Resolution probabilities versus the number of time samples (SNR is 10 dB).
Sensors 19 02410 g005
Figure 6. Root mean squared errors versus the number of time samples (SNR is 10 dB).
Figure 6. Root mean squared errors versus the number of time samples (SNR is 10 dB).
Sensors 19 02410 g006
Figure 7. Normalized pseudospectrums for three targets located at [ 8 , 1.5 , 7.5 ] .
Figure 7. Normalized pseudospectrums for three targets located at [ 8 , 1.5 , 7.5 ] .
Sensors 19 02410 g007
Figure 8. Resolution probabilities versus SNR ( N = 3 ).
Figure 8. Resolution probabilities versus SNR ( N = 3 ).
Sensors 19 02410 g008
Figure 9. Root mean squared errors versus SNR ( N = 3 ).
Figure 9. Root mean squared errors versus SNR ( N = 3 ).
Sensors 19 02410 g009
Figure 10. Resolution probabilities versus SNR ( N = 5 ).
Figure 10. Resolution probabilities versus SNR ( N = 5 ).
Sensors 19 02410 g010
Figure 11. Root mean squared errors versus SNR ( N = 5 ).
Figure 11. Root mean squared errors versus SNR ( N = 5 ).
Sensors 19 02410 g011
Figure 12. Measurement environment of two target vehicles located at [ 1.7 , 4.6 ].
Figure 12. Measurement environment of two target vehicles located at [ 1.7 , 4.6 ].
Sensors 19 02410 g012
Figure 13. Measurement environment of two target vehicles (expressway).
Figure 13. Measurement environment of two target vehicles (expressway).
Sensors 19 02410 g013
Table 1. Resolution probabilities and root mean squared errors for two adjacent targets located at [ 3.5 , 2.5 ] .
Table 1. Resolution probabilities and root mean squared errors for two adjacent targets located at [ 3.5 , 2.5 ] .
DOA Estimation Method P r (%)RMSE ( )
Conventional Bartlett04.28
Bartlett with T * and Y ^ ( t ) 74.92.67
Bartlett with V * and Z ^ ( t ) 99.40.64
Bartlett with V * and W ^ ( t ) 99.90.45
Table 2. Resolution probabilities and root mean squared errors for two target vehicles located at [ 1.7 , 4.6 ] .
Table 2. Resolution probabilities and root mean squared errors for two target vehicles located at [ 1.7 , 4.6 ] .
DOA Estimation Method P r (%)RMSE ( )
Conventional Bartlett05.72
Bartlett with T * and Y ^ ( t ) 0.676.20
Bartlett with V * and Z ^ ( t ) 46.74.65
Bartlett with V * and W ^ ( t ) 68.03.90
Table 3. Resolution probabilities and root mean squared errors for two target vehicles (expressway).
Table 3. Resolution probabilities and root mean squared errors for two target vehicles (expressway).
DOA Estimation Method P r (%)RMSE ( )
Conventional Bartlett07.83
Bartlett with T * and Y ^ ( t ) 1.256.57
Bartlett with V * and Z ^ ( t ) 43.34.98
Bartlett with V * and W ^ ( t ) 65.14.21

Share and Cite

MDPI and ACS Style

Lee, S.; Kim, S.-C. Logarithmic-Domain Array Interpolation for Improved Direction of Arrival Estimation in Automotive Radars. Sensors 2019, 19, 2410. https://doi.org/10.3390/s19102410

AMA Style

Lee S, Kim S-C. Logarithmic-Domain Array Interpolation for Improved Direction of Arrival Estimation in Automotive Radars. Sensors. 2019; 19(10):2410. https://doi.org/10.3390/s19102410

Chicago/Turabian Style

Lee, Seongwook, and Seong-Cheol Kim. 2019. "Logarithmic-Domain Array Interpolation for Improved Direction of Arrival Estimation in Automotive Radars" Sensors 19, no. 10: 2410. https://doi.org/10.3390/s19102410

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop