Abstract:
This invention is present an iterative method for joint antenna array calibration and direction of arrival estimation using millimeter-wave (mm-Wave) radar. The calibration compensates for array coupling, phase, and gain errors and does not require any training data. This method is well suited for applications where multiple antenna elements are packaged in a chip and where offline calibration is either expensive or is not possible. This invention is also effective when the array coupling is a function of direction of arriving waves from the object. It is also applicable to any two-dimensional array shape. Real experiment results demonstrate the viability of the algorithm using real data collected from a four-element array.

Description:
CLAIM OF PRIORITY 
       [0001]    This application claims priority under 35 U.S.C 119(e)(1) to U.S. Provisional Application No. 62/233,055 filed Sep. 25, 2015 
     
    
     TECHNICAL FIELD OF THE INVENTION 
       [0002]    The technical field of this invention is antenna array calibration. 
       BACKGROUND OF THE INVENTION 
       [0003]    This invention estimates the coupling between radar antenna array elements as well as phase and gain of each element. At the same time, the invention estimates the direction of arrival of objects in radar&#39;s field of view. The invention does not require any use of training data for antenna array calibration, and it is applicable to any array shape. 
       SUMMARY OF THE INVENTION 
       [0004]    An iterative method is shown for joint antenna array calibration and direction of arrival estimation on a millimeter wave radar. The calibration compensates for array coupling, phase and gain errors without the requirement for training data. The method shown is also effective when the array coupling is a function of the direction of the arriving waves, and is applicable to two dimensional array shapes. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0005]    These and other aspects of this invention are illustrated in the drawings, in which: 
           [0006]      FIG. 1  shows a block diagram of the invention; 
           [0007]      FIG. 2  illustrates an exemplary L-element antenna array; 
           [0008]      FIG. 3  is a block diagram of one example implementation; and 
           [0009]      FIG. 4  illustrates the estimated and true elevation and azimuth angles in degrees versus azimuth. 
       
    
    
     DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
       [0010]    Recent years have witnessed widespread use of millimeter wave (mm-Wave) radars for advanced driver assistance system (ADAS) applications. Compared with other sensing modalities such as a camera, radar has the ability to perform equally well during different times of the day and can be deployed out of sight behind the car bumper or the doors. In many ADAS applications such as parking, cruise control, and braking, the radar is primarily used to find the three-dimensional location of objects around the vehicle. This includes range, azimuth angle, and elevation angle. The range is computed from the round trip delay of the transmitted signal and the two-dimensional (2D) angle is estimated by using the data collected by an antenna array employing a beamforming-based or an eigen-decomposition based high-resolution frequency estimation method. It is well known that these directions of arrival (DOA) estimates are highly biased if the coupling between antenna array elements is not corrected and compensated for. Additionally, the phase and gain mismatch among antenna elements also adversely affects the estimation. The impact of these non-idealities is pronounced when the antennas are placed very close to each other, which generally is the case in automotive radars. 
         [0011]    Calibration for antenna coupling has been widely studied in the past and many methods have been proposed. Generally, training data is collected by L-element antenna array from radar field of view and a L×L calibration matrix is estimated; this matrix is then applied to signal received by the array. While this methodology works well for many cases, it is not very well suited for automotive radar applications for the following reasons. First, the antenna coupling for a 2D array changes with DOA from objects. It is, therefore, not possible to estimate a calibration matrix that would be applicable for all directions in radar field of view. Secondly, in automotive radars where multiple radars are placed on the car and are thus produced in high volume, it is desired that the calibration is done online without any need for training data. 
         [0012]    This invention is a joint array calibration and 2D angle estimation method of multiple objects around the vehicle. The method does not require any training data and needs minimal supervision. Whereas in the past, the joint estimation problem was solved for specific array shapes, the problem formulation and optimization proposed in this application can be applied to any array design and shape. We will present experimental results using data collected from a 77 GHz radar with a four-element antenna array to show efficacy of the proposed method. A simplified block diagram of such an implementation is shown in  FIG. 3 , where  301  is a processor,  302  is the radio frequency front end, and  303  is the antenna array. 
         [0013]    Though the analysis below applies to any 2D array shape. we consider a rectangular, L-element antenna array shown in  FIG. 2 , where L−1 elements  201  through  202  are in one direction, and the L th  element  203  is in the orthogonal direction. The array response at time n is given by 
         [0000]        x ( n )= BΓAs ( n ),  (1)
 
         [0000]    where x(n)=[x 1 (n), x 2 (n), . . . , xL(n)] T , B is a L×L coupling matrix, Γ={α 1 e −jωψ     1   , α 1 e −jωψ     2   , . . . , α 1 e −jωψ     L   , } is L×L matrix with antenna gains and phase values as diagonal elements. A=[a(θ 1 , φ 1 ), a(θ 2 , φ 2 ), . . . , a(θ K , φ K )] is a L×K matrix of steering vectors with a(θ l , φ l ) given by 
         [0000]    
       
         
           
             
               
                 [ 
                 
                   
                     
                       
                         1 
                         , 
                         
                            
                           
                             
                               - 
                               j 
                             
                              
                             
                               
                                 2 
                                  
                                 π 
                               
                               λ 
                             
                              
                             
                               d 
                               x 
                             
                              
                             
                               sin 
                                
                               
                                 ( 
                                 
                                   θ 
                                   l 
                                 
                                 ) 
                               
                             
                              
                             
                               cos 
                                
                               
                                 ( 
                                 
                                   φ 
                                   l 
                                 
                                 ) 
                               
                             
                           
                         
                         , 
                         
                            
                           
                             
                               - 
                               j 
                             
                              
                             
                               
                                 2 
                                  
                                 
                                   ( 
                                   
                                     L 
                                     - 
                                     1 
                                   
                                   ) 
                                 
                                  
                                 π 
                               
                               λ 
                             
                              
                             
                               d 
                               x 
                             
                              
                             
                               sin 
                                
                               
                                 ( 
                                 
                                   θ 
                                   l 
                                 
                                 ) 
                               
                             
                              
                             
                               cos 
                                
                               
                                 ( 
                                 
                                   φ 
                                   l 
                                 
                                 ) 
                               
                             
                           
                         
                         , 
                         … 
                       
                     
                   
                   
                     
                       
                          
                         
                           - 
                           
                             j 
                              
                             
                               ( 
                               
                                 
                                   
                                     
                                       2 
                                        
                                       π 
                                     
                                     λ 
                                   
                                    
                                   
                                     d 
                                     x 
                                   
                                    
                                   
                                     sin 
                                      
                                     
                                       ( 
                                       
                                         θ 
                                         l 
                                       
                                       ) 
                                     
                                   
                                    
                                   
                                     cos 
                                      
                                     
                                       ( 
                                       
                                         φ 
                                         l 
                                       
                                       ) 
                                     
                                   
                                 
                                 + 
                                 
                                   
                                     
                                       2 
                                        
                                       π 
                                     
                                     λ 
                                   
                                    
                                   
                                     d 
                                     y 
                                   
                                    
                                   
                                     sin 
                                      
                                     
                                       ( 
                                       
                                         φ 
                                         l 
                                       
                                       ) 
                                     
                                   
                                 
                               
                               ) 
                             
                           
                         
                       
                     
                   
                 
                 ] 
               
               T 
             
             , 
           
         
       
     
         [0000]    s(n)=[s 1 (n), s 2 (n), . . . , s K (n)] T  is complex signal amplitude, K is the number of objects, θ is azimuth direction of arrival, φ is elevation direction of arrival, and (•) T  denotes conjugate transpose. The parameters α l  and ψ l  are gain and delay associated with lth sensor. Biased DOA estimates are obtained if the effects of B and Γ are not compensated for in the received signal. 
         [0014]    Let λ i  and u i , i=1, 2, . . . , L be the eigenvalues and eigenvectors of the sample covariance matrix R x =E{xx H }. Collecting the set of eigenvectors belonging to noise subspace in matrix U=[u K+1 , u K+2 , . . . , u L ], the unknowns Γ, B, θ, and φ are obtained by minimizing the cost function: 
         [0000]    
       
         
           
             
               
                 
                   
                     J 
                     = 
                     
                       
                         ∑ 
                         
                           k 
                           = 
                           1 
                         
                         K 
                       
                        
                       
                           
                       
                        
                       
                         
                            
                           
                             
                               U 
                               H 
                             
                              
                             B 
                              
                             
                                 
                             
                              
                             Γ 
                              
                             
                                 
                             
                              
                             
                               a 
                                
                               
                                 ( 
                                 
                                   
                                     θ 
                                     k 
                                   
                                   , 
                                   
                                     φ 
                                     k 
                                   
                                 
                                 ) 
                               
                             
                           
                            
                         
                         2 
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   2 
                   ) 
                 
               
             
           
         
       
     
         [0000]               where denotes Frobenius norm. 
         [0015]    In the proposed method, the matrices Γ and B, and the angles (θ, φ) are estimated using the following iterative method as shown in  FIG. 1 : 
         [0016]    1) Initialization ( 101 ); i=0; Set B i  and Γ i  to initial values. 
         [0017]    2) Estimate data covariance matrix ( 102 ) 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       R 
                       ^ 
                     
                     x 
                   
                   = 
                   
                     
                       1 
                       N 
                     
                      
                     
                       
                         ∑ 
                         
                           n 
                           = 
                           1 
                         
                         N 
                       
                        
                       
                           
                       
                        
                       
                         
                           x 
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                          
                         
                           
                             x 
                              
                             
                               ( 
                               n 
                               ) 
                             
                           
                           H 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
           
         
       
     
         [0018]    3) Eigen-decompose {circumflex over (R)} x  ( 103 ), find U and search for K peaks ( 104 ) in the 2D spectrum defined by 
         [0000]        P   i (θ, φ)=∥ U   H   B   i Γ i   a (θ, φ)∥ −2   (4)
 
         [0000]    The peaks of equation (4) correspond to the DOA estimates {(θ k , φ k )} k=1   K . 
         [0019]    4) Use the estimated DOA to form the matrix A. ( 105 ) 
         [0020]    5) Under the constraint γ H w=1, where w=[1, 0, 0, . . . , 0] T , estimate L×1 vector γ ( 106 ) 
         [0000]      γ= Z   −1   w /( w   T   Z   −1   w ),  (5)
 
         [0000]    where the matrix Z is given by = 
         [0000]    
       
         
           
             
               
                 
                   Z 
                   = 
                   
                     
                       ∑ 
                       
                         k 
                         = 
                         1 
                       
                       K 
                     
                      
                     
                         
                     
                      
                     
                       
                         Q 
                         k 
                         H 
                       
                        
                       
                         B 
                         
                           i 
                           H 
                         
                       
                        
                       
                         UU 
                         H 
                       
                        
                       
                         B 
                         i 
                       
                        
                       
                         Q 
                         k 
                       
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
           
         
       
     
         [0000]    and the diagonal matrix Qk is formed using Qk=diag{a(θ k , φ k )}. 
         [0021]    6) Update the estimate of Γ ( 107 ) using the diagonal elements of γ as Γ i     +1   =diag{γ}. 
         [0022]    7) Finally, the cost function in equation (2) is minimized ( 108 ) in the least squares sense to solve for B under the constraint B 11 =1. This optimization is carried out as follows: 
         [0023]    Compute KL×L 2  matrix M=(A T B iT )          U T , where           defines the Kronecker product. 
         [0024]    Extract M 1 =M(:, 1) and M 2 =M(:, 2:end); i.e, M 1  contains only the first column of M and M 2  is M, except for its first column. 
         [0025]    Compute (L 2 −1)×1 vector b=−M 2   # M 1 , where (•) #  denotes pseudo-inverse. 
         [0026]    Compute L 2 ×1 vector {tilde over (b)}=[1 b T ] T . 
         [0027]    Re-arrange {tilde over (b)} in rows of L to form updated B i     +1   . 
         [0028]    The iterative calibration and angle estimation is continued until the cost function at (i+1)th iteration is smaller than what it was at the ith iteration by a pre-set threshold. 
         [0029]      FIG. 4  illustrates the case using a four-element array (L=4) with three elements in one direction and the fourth element in orthogonal direction with a carrier frequency of 77 GHz. The inter-element spacing in either direction was 2 mm. An object was placed at an elevation angle of −20 degrees and moved along azimuth direction from −40 degrees to 40 degrees in increments of 5 degrees. The matrices B and Γ were initialized with identity matrices. At each location of the object, the joint iterative algorithm was used to estimate the DOA&#39;s. It took a maximum of 5 iterations for the algorithm to converge to an estimate.  FIG. 4  illustrates the azimuth and elevation angles plotted against the azimuth angles. 
         [0030]    A joint calibration and angle estimation algorithm is presented. The method is especially suited for automotive applications where multiple sensors are installed around the vehicle and online calibration and angle estimation is highly desired. The method does not impose any constraint on the array shape and structure and takes only a few iterations to converge.