Abstract:
A radar system includes: a transmission antenna outputting transmission signals having multiple frequencies; multiple reception antennas receiving reflected waves of the transmission signals, reflected from an object; a mixer mixing the transmission signals with reception signals received by the reception antennas to generate beat signals; and a signal processing unit detecting Doppler frequency by analyzing frequencies of the beat signals, detecting phase information of the Doppler frequency for each of combinations of the reception antennas and the transmission signal frequencies, constructing a matrix having the pieces of phase information arranged in a predetermined order with respect to the reception antennas and the frequencies of the transmission signals, obtaining a correlation matrix from the matrix and its complex conjugate transposed matrix, and estimating at least one of a distance, direction and relative velocity of the object based on the correlation matrix.

Description:
INCORPORATION BY REFERENCE 
       [0001]    The disclosure of Japanese Patent Application No. 2009-156482 filed on Jul. 1, 2009 including the specification, drawings and abstract is incorporated herein by reference in its entirety. 
       BACKGROUND OF THE INVENTION 
       [0002]    1. Field of the Invention 
         [0003]    The invention relates to a CW radar system that uses a plurality of reception antennas, and a signal processing method for the radar system. 
         [0004]    2. Description of the Related Art 
         [0005]    In order to detect a distance to a stationary or moving object, a direction to the object and a moving velocity of the object, various radar systems have been developed. 
         [0006]    For example, Japanese Patent Application Publication No. 2008-145425 (JP-A-2008-145425) describes a radar system that outputs transmission signals having three or more different frequencies from an oscillator, receives signals reflected from a target, mixes the reception signals with the transmission signals by a mixer to generate beat signals, detects Doppler frequency signals from the beat signals through fast Fourier transform (FFT), or the like, and then obtains a distance to the target on the basis of complex signal components of the Doppler frequency signals of the respective transmission signals. 
         [0007]    In the CW radar system, reflected signals from the target in correspondence with the transmission signals having a plurality of frequencies are received by a plurality of reception antennas and are analyzed. Then, in the CW radar system, in order to obtain distance information to the target in high resolution, pieces of phase information obtained from the respective reception antennas (reception channels) are used to estimate a distance to the target. In this estimation, a correlation matrix is calculated using phase information for each reception channel. Therefore, as the number of reception channels increases, a computational load increases. 
         [0008]    In addition, when there is only one target, it is obvious that beat signals obtained through the respective reception channels are signals from the same target, so it is possible to accurately estimate a direction to the target from phase differences between the reception channels. However, when there are a plurality of targets having different relative velocities, beat signals of the number of targets are detected through each reception channel, and it is necessary to associate (pair) the beat signals among the reception channels. 
         [0009]    For example, when there are two targets having different relative velocities, two beat signals are detected through each of the two reception channels. Where the beat signals detected through a reception channel  1  are L 1  and L 2 , and the beat signals detected through a reception channel  2  are R 1  and R 2 , there are two combination patterns of the beat signals between the reception channels  1  and  2 , that is, ( 1 ) (L 1 , R 1 ) and (L 2 , R 2 ) or ( 2 ) (L 1 , R 2 ) and (L 2 , R 1 ). Here, if an erroneous combination is made, the directions to the targets are also erroneously estimated. In addition, when the number of targets increases, a processing load for associating beat signals increases. 
       SUMMARY OF THE INVENTION 
       [0010]    A first aspect of the invention provides a radar system. The radar system includes: a transmission antenna that outputs transmission signals having a plurality of frequencies as transmission waves; a plurality of reception antennas that receive reflected waves of the transmission signals, reflected from an object; a mixer that mixes the transmission signals with reception signals received by the reception antennas to generate beat signals of the reception signals received by the respective reception antennas for each of the transmission signals; and a signal processing unit that detects Doppler frequency by analyzing frequencies of the beat signals, that detects phase information of the Doppler frequency for each of combinations of the reception antennas and the frequencies of the transmission signals, that constructs a matrix, in which the pieces of phase information are arranged in a predetermined order with respect to the reception antennas and the frequencies of the transmission signals, that obtains a correlation matrix from the matrix and a complex conjugate transposed matrix of the matrix, and that estimates at least one of a distance to the object, a direction to the object and a relative velocity of the object on the basis of the correlation matrix. 
         [0011]    Here, the signal processing unit may estimate the at least one of the distance to the object, the direction to the object and the relative velocity of the object after the correlation matrix has been averaged by at least one of forward-backward averaging and spatial moving average. 
         [0012]    A second aspect of the invention provides a signal processing method for a radar system that includes a transmission antenna that outputs transmission signals having a plurality of frequencies as transmission waves and a plurality of reception antennas that receive reflected waves of the transmission signals, reflected from an object. The signal processing method includes: mixing the transmission signals with reception signals received by the reception antennas to generate beat signals of the reception signals received by the respective reception antennas for each of the transmission signals having the plurality of frequencies; detecting Doppler frequency by analyzing frequencies of the beat signals; detecting phase information of the Doppler frequency for each of combinations of the reception antennas and the frequencies of the transmission signals; constructing a matrix, in which the pieces of phase information are arranged in a predetermined order with respect to the reception antennas and the frequencies of the transmission signals; obtaining a correlation matrix from the matrix and a complex conjugate transposed matrix of the matrix; and estimating at least one of a distance to the object, a direction to the object and a relative velocity of the object on the basis of the correlation matrix. 
         [0013]    According to the aspects of the invention, it is possible to reduce processing load on the radar system. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0014]    The foregoing and further objects, features and advantages of the invention will become apparent from the following description of example embodiments with reference to the accompanying drawings, wherein like numerals are used to represent like elements and wherein: 
           [0015]      FIG. 1  is a view that shows the configuration of a radar system according to an embodiment of the invention; 
           [0016]      FIG. 2  is a graph that shows changes in frequency of transmission signals according to the embodiment of the invention; and 
           [0017]      FIG. 3  is a view that shows an example of analyzing frequencies of reception signals according to the embodiment of the invention. 
       
    
    
     DETAILED DESCRIPTION OF EMBODIMENTS 
     System Configuration 
       [0018]    As shown in  FIG. 1 , a radar system  100  according to an embodiment of the invention includes an oscillator  10 , a directional coupler  12 , a transmission antenna  14 , reception antennas  16 - k  (k is an integer larger than or equal to 2), a switch  18 , a mixer  20 , a band pass filter (BPF)  22 , an analog/digital converter (ADC)  24  and a signal processing unit  26 . 
         [0019]    The oscillator  10  generates and outputs transmission signals. The transmission signals are radiated from the transmission antenna  14  as transmission waves. The oscillator  10  is able to change the oscillatory frequency. In the present embodiment, the oscillator  10  generates and outputs N types (where N is 2 or above) of continuous waves respectively having a fundamental frequency f 0  to a frequency f 0 +(N−1)Δf at a predetermined frequency interval Δf. When N is 3, the oscillator  10  outputs transmission waves respectively having frequencies f 0 , f 0 +Δf and f 0 +2Δf. 
         [0020]    The directional coupler  12  demultiplexes the transmission signals output from the oscillator  10 , and outputs the demultiplexed transmission signals to both the transmission antenna  14  and the mixer  20 . The transmission antenna  14  outputs the transmission signals demultiplexed by the directional coupler  12  into space in a radiation pattern corresponding to the antenna characteristic. As shown in  FIG. 2 , transmission waves having frequencies of the fundamental frequency f 0  to the frequency f 0 +(N−1)Δf are sequentially and repeatedly transmitted from the transmission antenna  14  at a period of time T. 
         [0021]    The reception antennas  16 - k  each receive radio waves in accordance with the antenna characteristics from space. At least two or more reception antennas  16 - k  are provided (k is an integer larger than or equal to 2). In the present embodiment, K reception antennas  16 - 1  to  16 -K are provided. The reception antennas  16 - k  are spaced apart from each other. A reception signal received by each reception antenna  16 - k  includes components of reflected waves that a target  200  reflects the transmission signals radiated from the transmission antenna  14 . The frequencies of reflected waves shift from the frequencies of the transmission signals by a Doppler frequency in accordance with a relative velocity between the radar system  100  and the target  200 . Hereinafter, the reception antennas  16 - 1  to  16 -K may be expressed as reception channels ch 1  to chK. 
         [0022]    The switch  18  exclusively switches among reception signals received by the respective reception antennas  16 - 1  to  16 -K, and then outputs any one of the reception signals to the mixer  20 . By so doing, the reception signals received by the respective reception antennas  16 - 1  to  16 -K are sequentially output from the switch  18 . That is, transmission waves having frequencies of the fundamental frequency f 0  to the frequency f 0 +(N−1)Δf are sequentially radiated, signals containing components of reflected waves reflected by the target  200  are received by the reception antennas  16 - 1  to  16 -K, and then a reception signal received by one of the reception antennas  16 - 1  to  16 -K, selected by the switch  18 , is sequentially output to the mixer  20 . 
         [0023]    The mixer  20  mixes the transmission signal output from the directional coupler  12  with any one of the reception signals of the reception channels ch 1  to chK, output from the switch  18 , and outputs the mixed signal to the BPF  22 . The signal output from the mixer  20  contains a beat signal having a frequency corresponding to a difference between the frequency of the transmission signal and the frequency of the reception signal. That is, when there is a relative velocity between the target  200  and the radar system  100 , there occurs a frequency shift due to Doppler effect. This causes a difference in frequency between the transmission signal and the reception signal. A signal having a frequency corresponding to this difference is output as a beat signal. 
         [0024]    The BPF  22  removes an unnecessary signal, other than a component of a beat signal that indicates a frequency shift due to Doppler effect, from a signal generated by the mixer  20 , and then outputs the resultant signal to the ADC  24 . The ADC  24  converts the signal output from the BPF  22  from an analog signal into a digital signal and outputs the converted signal to the signal processing unit  26 . 
         [0025]    The signal processing unit  26  receives an output signal from the ADC  24 , and then estimates, for example, a distance from the radar system  100  to the target  200 , a direction from the radar system  100  to the target  200  and a relative velocity between the radar system  100  and the target  200  on the basis of the output signal. The signal processing unit  26  may be implemented by executing a program, which executes the following arithmetic processing, in a general computer provided with a CPU, a memory, an input/output device, and the like. Alternatively, the signal processing unit  26  may be formed of a logic circuit that executes the following arithmetic processing. 
         [0026]    Note that, in the present embodiment, a signal digitized by the ADC  24  is processed; instead, it is also applicable that the signal processing unit  26  is formed of an analog circuit and then an analog signal is directly processed. 
       Signal Processing 
       [0027]    Hereinafter, signal processing executed by the radar system  100  will be described. The following process will be executed by the signal processing unit  26 . Note that there may be a plurality of targets  200  and it is assumed that the location and velocity of each target  200  do not change throughout the observing time. 
         [0028]    The signal processing unit  26  obtains a frequency spectrum on the basis of a signal received from the ADC  24  through fast Fourier transform, or the like.  FIG. 3  shows an example in which, while the transmission signals are being transmitted, frequency spectra of beat signals generated by the mixer  20  for reception signals of the reception antennas  16 - k  (reception channels chk) that have received reflected waves from the targets  200  are obtained. Here, the transmission signals respectively having N (where N is 2 or above) types of frequencies of the fundamental frequency f 0  to the frequency f 0 +(N−1)Δf at the frequency interval Δf are transmitted. When there are a plurality of targets  200  having different velocities, respective reflected waves have different Doppler frequencies with respect to the radar system  100 , so signals of Doppler frequencies for respective velocities appear. In addition, for reflected waves of the targets  200  having no relative velocities with respect to the radar system  100 , the outputs of the mixer  20  are direct-current components and then the direct-current components are removed by the BPF  22 . 
         [0029]    In the example of  FIG. 3 , for respective transmission signals of the fundamental frequency f 0  to the frequency f 0 +(N−1)Δf, Doppler frequencies f 1  to f m  generated on the basis of the relative velocities between the targets  200  and the radar system  100  each have a peak. As shown in  FIG. 3 , the Doppler frequencies f 1  to f m , change in proportion to not only the relative velocities between the targets  200  and the radar system  100  but also the frequencies f 0  to f 0 +(N−1)Δf of the transmission signals. For example, in 76 GHz millimeter wave band, the Doppler frequency only changes by 1.3% even when the frequency changes by 1 GHz. Thus, differences in frequency between the transmission signals almost do not influence the Doppler frequencies f 1  to f m . 
         [0030]    The following analysis is applied to each of the thus obtained Doppler frequencies f 1  to f m , and then the distances, directions and relative velocities to the targets  200  corresponding to the respective Doppler frequencies f 1  to f m , are estimated. 
         [0031]    First, a complex signal component (phase information) of the spectrum of each Doppler frequency f j  (j is an integer ranging from 1 to m and specifies the Doppler frequency) is detected for each of combinations of the reception antennas  16 - 1  to  16 -K (reception channels ch 1  to chK) and the frequencies f 0  to f 0 +(N−1)Δf of the transmission signals. Then, the complex signal components (a pieces of phase information) of the spectra of the respective Doppler frequencies f 1  are arranged in predetermined orders with respect to the reception antennas  16 - 1  to  16 -K (reception channels ch 1  to chK) and the frequencies f 0  to f 0 +(N−1)Δf of the transmission signals to construct a matrix B j . 
         [0032]    The predetermined order with respect to the reception antennas  16 - 1  to  16 -K (reception channels ch 1  to chK) are desirably an order in which, for example, the switch  18  switches among the reception antennas  16 - 1  to  16 -K. More specifically, the predetermined order is desirably the order of the reception antenna  16 - 1 , the reception antenna  16 - 2 , . . . , the reception antenna  16 -K. In addition, the predetermined order with respect to the frequencies f 0  to f 0 +(N−1)Δf of the transmission signals is desirably an order in which, for example, the oscillator  10  generates the frequencies of the transmission signals. More specifically, the predetermined order is desirably the order of the frequency f 0 , the frequency f 0 +Δf, . . . , the frequency f 0 +(N−1)Δf. However, the predetermined order is not limited to the above; it is only necessary that the respective orders in each row and each column of the matrix B j  are kept unchanged. 
         [0033]    When the above predetermined orders are applied, as shown in the mathematical expression (1), an element b nk  of the matrix B j  is a complex signal component (phase information) of the Doppler frequency f j  in the frequency spectrum obtained by analyzing the reception signal received by the reception antenna  16 - k  (reception channel chk) while the transmission signal having the frequency f 0 +(n−1)Δf is being transmitted. That is, n is an integer ranging from 1 to N for specifying the frequency f 0 +(n−1)Δf of the transmission signal. In addition, k is an integer ranging from 1 to K for specifying the reception antenna  16 - k  (reception channel chk). 
         [0000]    
       
         
           
             
               
                 
                   
                     B 
                     j 
                   
                   = 
                   
                     ( 
                     
                       
                         
                           
                             b 
                             11 
                           
                         
                         
                           
                             b 
                             12 
                           
                         
                         
                           Λ 
                         
                         
                           
                             b 
                             
                               1 
                                
                               k 
                             
                           
                         
                         
                           Λ 
                         
                         
                           
                             b 
                             
                               1 
                                
                               K 
                             
                           
                         
                       
                       
                         
                           
                             b 
                             21 
                           
                         
                         
                           ⋱ 
                         
                         
                           Λ 
                         
                         
                           
                             b 
                             
                               2 
                                
                               k 
                             
                           
                         
                         
                           Λ 
                         
                         
                           
                             b 
                             
                               2 
                                
                               K 
                             
                           
                         
                       
                       
                         
                           ⋮ 
                         
                         
                           
                               
                           
                         
                         
                           ⋱ 
                         
                         
                           
                               
                           
                         
                         
                           
                               
                           
                         
                         
                           ⋮ 
                         
                       
                       
                         
                           
                             b 
                             
                               n 
                                
                               
                                   
                               
                                
                               1 
                             
                           
                         
                         
                           
                               
                           
                         
                         
                           
                               
                           
                         
                         
                           ⋱ 
                         
                         
                           
                               
                           
                         
                         
                           ⋮ 
                         
                       
                       
                         
                           ⋮ 
                         
                         
                           
                               
                           
                         
                         
                           
                               
                           
                         
                         
                           
                               
                           
                         
                         
                           ⋱ 
                         
                         
                           ⋮ 
                         
                       
                       
                         
                           
                             b 
                             
                               N 
                                
                               
                                   
                               
                                
                               1 
                             
                           
                         
                         
                           
                               
                           
                         
                         
                           Λ 
                         
                         
                           
                               
                           
                         
                         
                           
                               
                           
                         
                         
                           
                             b 
                             NK 
                           
                         
                       
                     
                     ) 
                   
                 
               
               
                 
                   ( 
                   1 
                   ) 
                 
               
             
           
         
       
     
         [0034]    For example, when N and K each are 3, the matrix B 1  corresponding to the Doppler frequency f 1  has three rows and three columns as shown in the mathematical expression (2). The element b 11  is a complex signal component (phase information) of the Doppler frequency f 1  in the frequency spectrum obtained by analyzing the reception signal received by the reception antenna  16 - 1  (reception channel ch 1 ) while the transmission signal having the frequency f 0  is being transmitted. In addition, the element b 12  is a complex signal component (phase information) of the Doppler frequency f 1  in the frequency spectrum obtained by analyzing the reception signal received by the reception antenna  16 - 2  (reception channel ch 2 ) while the transmission signal having the frequency f 0  is being transmitted. In addition, the element b 21  is a complex signal component (phase information) of the Doppler frequency f 1  in the frequency spectrum obtained by analyzing the reception signal received by the reception antenna  16 - 1  (reception channel ch 1 ) while the transmission signal having the frequency f 0 +Δf is being transmitted. The other elements are also similar to the above elements. 
         [0000]    
       
         
           
             
               
                 
                   
                     B 
                     j 
                   
                   = 
                   
                     ( 
                     
                       
                         
                           
                             b 
                             11 
                           
                         
                         
                           
                             b 
                             12 
                           
                         
                         
                           
                             b 
                             13 
                           
                         
                       
                       
                         
                           
                             b 
                             21 
                           
                         
                         
                           
                             b 
                             22 
                           
                         
                         
                           
                             b 
                             23 
                           
                         
                       
                       
                         
                           
                             b 
                             31 
                           
                         
                         
                           
                             b 
                             32 
                           
                         
                         
                           
                             b 
                             33 
                           
                         
                       
                     
                     ) 
                   
                 
               
               
                 
                   ( 
                   2 
                   ) 
                 
               
             
           
         
       
     
         [0035]    In the matrix B j , the element b nk  of the column vector, which corresponds to the reception antenna  16 - k  (reception channel chk), indicates a complex signal component (phase information) of the Doppler frequency f j  in each of the frequencies f 0  to f 0 +(N−1)Δf of the transmission signals. Thus, the phase differences between the elements b nk  of the column vector occur because of the frequencies f 0  to f 0 +(N−1)Δf of the transmission signals, and do not depend on the location of the reception antenna  16 - k . In addition, phase differences due to optical path differences between the reception antennas  16 - 1  to  16 -K and each target  200  depend on the locations of the reception antennas  16 - 1  to  16 -K. Thus, the phase differences between the elements b np  of the column vector with respect to a selected reception antenna  16 - p  (p is any one of integers ranging from 1 to K) is equal to the phase differences between the elements b nq  of the column vector with respect to another reception antenna  16 - q  (q is any one of integers ranging from 1 to K other than p). 
         [0036]    Where the phase differences between the elements of the column vector obtained from a selected reception antenna are denoted by a reference vector C j  and the phase differences due to optical path differences caused by the locations of the reception antennas are denoted by a vector D j  the matrix B j  may be expressed as C j ×D j  from the above described characteristic. 
         [0037]    Then, a correlation matrix Rxx j  for the matrix B j  may be expressed as the mathematical expression (3). Note that the matrix B j   H , the vector C j   H  and the vector D j   H  respectively denote complex conjugate transposed matrices (vectors) of the matrix B j , reference vector C j  and vector D j . 
         [0000]        Rxx   J   =B   J   ×B   J   H   =C   J   ×D   J   ×D   j   H   ×C   J   H   (3)
 
         [0038]    Here, D j ×D j   H  is a constant α j , so the mathematical expression (3) may be further transformed into the mathematical expression (4). 
         [0000]        Rxx   J   =B   J   ×B   J   H =α J   ×C   J   ×C   j   H   (4)
 
         [0039]    The mathematical expression (4) indicates that a mathematical expression for obtaining the correlation matrix Rxx j  is the same as a mathematical expression for obtaining a correlation matrix using the column vector of each reception antenna  16 - k  (reception channel chk). However, the correlation matrix Rxx j  contains complex signal components (phase information) of the Doppler frequencies f 1  obtained by all the reception antennas  16 - 1  to  16 -K (all the reception channels ch 1  to chK), so the S/N ratio of a signal spectrum obtained thereafter for the correlation matrix Rxx j  is higher than that of the correlation matrix obtained for each reception antenna  16 - k  (reception channel chk). 
         [0040]    The thus obtained correlation matrix Rxx j  is utilized to estimate information about each target  200 . A high-resolution estimation method, such as the MUSIC method, the ESPIRIT method and the Capon method, may be desirably employed. 
         [0041]    Hereinafter, a distance estimation method using the Capon method will be described as an example. In the Capon method, a mathematical expression for calculating a spectrum amplitude is expressed as the mathematical expression (5). Here, a(r) is a mode vector that depends on a distance r, for which a spectrum is obtained, and the frequencies f 0  to f 0 +(N−1)Δf of the transmission signals, and a(r) H  is a complex conjugate transposed matrix of a(r). However, the elements of a(r) are arranged in the order of the frequencies of the matrix B j . 
         [0000]    
       
         
           
             
               
                 
                   
                     Pw 
                      
                     
                       ( 
                       r 
                       ) 
                     
                   
                   = 
                   
                     
                       
                         1 
                         
                           
                             
                               a 
                               11 
                             
                              
                             
                               ( 
                               r 
                               ) 
                             
                           
                            
                           
                             Rxx 
                             j 
                             
                               - 
                               1 
                             
                           
                            
                           
                             a 
                              
                             
                               ( 
                               r 
                               ) 
                             
                           
                         
                       
                        
                       
                           
                       
                        
                       
                         a 
                          
                         
                           ( 
                           r 
                           ) 
                         
                       
                     
                     = 
                     
                       ( 
                       
                         
                           
                             
                                
                               
                                 j2π 
                                  
                                 
                                   
                                     2 
                                      
                                     r 
                                   
                                   c 
                                 
                                  
                                 
                                   f 
                                   0 
                                 
                               
                             
                           
                         
                         
                           
                             
                                
                               
                                 j2π 
                                  
                                 
                                   
                                     2 
                                      
                                     r 
                                   
                                   c 
                                 
                                  
                                 
                                   ( 
                                   
                                     
                                       f 
                                       0 
                                     
                                     + 
                                     
                                       Δ 
                                        
                                       
                                           
                                       
                                        
                                       f 
                                     
                                   
                                   ) 
                                 
                               
                             
                           
                         
                         
                           
                             M 
                           
                         
                         
                           
                             
                                
                               
                                 j2π 
                                  
                                 
                                   
                                     2 
                                      
                                     r 
                                   
                                   c 
                                 
                                  
                                 
                                   ( 
                                   
                                     
                                       f 
                                       0 
                                     
                                     + 
                                     
                                       
                                         ( 
                                         
                                           N 
                                           - 
                                           1 
                                         
                                         ) 
                                       
                                        
                                       Δ 
                                        
                                       
                                           
                                       
                                        
                                       f 
                                     
                                   
                                   ) 
                                 
                               
                             
                           
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   ( 
                   5 
                   ) 
                 
               
             
           
         
       
     
         [0042]    The mathematical expression (5) is used while changing the distance r at a selected distance interval to obtain power Pw(r), and then the distance r at which the power Pw(r) indicates a peak is estimated as the distance to the target  200 . 
         [0043]    The above process is carried out for each of the Doppler frequencies f 1  to f m  to thereby make it possible to estimate the distance and direction to the target  200 , and the relative velocity of the target  200 , which cause the peak of the spectrum to be formed for each of the Doppler frequencies f 1  to f m . 
       Alternative Embodiment 
       [0044]    When a correlation between the elements of the matrix B j  is high because the observing time is short, for example, the correlation matrix Rxx j  may be subjected to averaging. For example, averaging, such as forward-backward averaging and spatial moving average, may be applied to the correlation matrix Rxx j . These processes may be applied solely or in combination. 
         [0045]    A specific example of a method of calculating a forward-backward average for a correlation matrix Ru is shown by the mathematical expression (6). Note that r* denotes a complex conjugate of r. 
         [0000]    
       
         
           
             
               
                 
                   Ru 
                   = 
                   
                     
                       ( 
                       
                         
                           
                             
                               r 
                               11 
                             
                           
                           
                             
                               r 
                               12 
                             
                           
                           
                             
                               r 
                               13 
                             
                           
                         
                         
                           
                             
                               r 
                               21 
                             
                           
                           
                             
                               r 
                               22 
                             
                           
                           
                             
                               r 
                               23 
                             
                           
                         
                         
                           
                             
                               r 
                               31 
                             
                           
                           
                             
                               r 
                               32 
                             
                           
                           
                             
                               r 
                               33 
                             
                           
                         
                       
                       ) 
                     
                     ⇒ 
                     
                       Rus 
                        
                       
                         ( 
                         
                           
                             
                               
                                 
                                   r 
                                   11 
                                 
                                 + 
                                 
                                   r 
                                   33 
                                   * 
                                 
                               
                             
                             
                               
                                 
                                   r 
                                   12 
                                 
                                 + 
                                 
                                   r 
                                   32 
                                   * 
                                 
                               
                             
                             
                               
                                 
                                   r 
                                   13 
                                 
                                 + 
                                 
                                   r 
                                   31 
                                   * 
                                 
                               
                             
                           
                           
                             
                               
                                 
                                   r 
                                   21 
                                 
                                 + 
                                 
                                   r 
                                   23 
                                   * 
                                 
                               
                             
                             
                               
                                 
                                   r 
                                   22 
                                 
                                 + 
                                 
                                   r 
                                   22 
                                   * 
                                 
                               
                             
                             
                               
                                 
                                   r 
                                   23 
                                 
                                 + 
                                 
                                   r 
                                   23 
                                   * 
                                 
                               
                             
                           
                           
                             
                               
                                 
                                   r 
                                   31 
                                 
                                 + 
                                 
                                   r 
                                   13 
                                   * 
                                 
                               
                             
                             
                               
                                 
                                   r 
                                   32 
                                 
                                 + 
                                 
                                   r 
                                   12 
                                   * 
                                 
                               
                             
                             
                               
                                 
                                   r 
                                   33 
                                 
                                 + 
                                 
                                   r 
                                   11 
                                   * 
                                 
                               
                             
                           
                         
                         ) 
                       
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
           
         
       
     
         [0046]    In addition, in the moving average, a plurality of sub-arrays are defined along a diagonal line of the correlation matrix Rxx j , and then those components are averaged to calculate a new matrix. A specific example of the moving average for the correlation matrix Ru is shown by the mathematical expression (7). 
         [0000]    
       
         
           
             
               
                 
                   Ru 
                   = 
                   
                     
                       ( 
                       
                         
                           
                             
                               r 
                               11 
                             
                           
                           
                             
                               r 
                               12 
                             
                           
                           
                             
                               r 
                               13 
                             
                           
                         
                         
                           
                             
                               r 
                               21 
                             
                           
                           
                             
                               r 
                               22 
                             
                           
                           
                             
                               r 
                               23 
                             
                           
                         
                         
                           
                             
                               r 
                               31 
                             
                           
                           
                             
                               r 
                               32 
                             
                           
                           
                             
                               r 
                               33 
                             
                           
                         
                       
                       ) 
                     
                     ⇒ 
                     
                       Rus 
                        
                       
                         ( 
                         
                           
                             
                               
                                 
                                   r 
                                   11 
                                 
                                 + 
                                 
                                   r 
                                   22 
                                 
                               
                             
                             
                               
                                 
                                   r 
                                   12 
                                 
                                 + 
                                 
                                   r 
                                   23 
                                 
                               
                             
                           
                           
                             
                               
                                 
                                   r 
                                   21 
                                 
                                 + 
                                 
                                   r 
                                   32 
                                 
                               
                             
                             
                               
                                 
                                   r 
                                   22 
                                 
                                 + 
                                 
                                   r 
                                   33 
                                 
                               
                             
                           
                         
                         ) 
                       
                     
                   
                 
               
               
                 
                   ( 
                   7 
                   ) 
                 
               
             
           
         
       
     
         [0047]    Here, a sub-array  1  S 1  and a sub-array  2  S 2  are respectively defined as follows. 
         [0000]    
       
         
           
             
               S 
               1 
             
             = 
             
               
                 
                   ( 
                   
                     
                       
                         
                           r 
                           11 
                         
                       
                       
                         
                           r 
                           12 
                         
                       
                     
                     
                       
                         
                           r 
                           21 
                         
                       
                       
                         
                           r 
                           22 
                         
                       
                     
                   
                   ) 
                 
                  
                 
                     
                 
                  
                 
                   S 
                   2 
                 
               
               = 
               
                 ( 
                 
                   
                     
                       
                         r 
                         22 
                       
                     
                     
                       
                         r 
                         23 
                       
                     
                   
                   
                     
                       
                         r 
                         32 
                       
                     
                     
                       
                         r 
                         33 
                       
                     
                   
                 
                 ) 
               
             
           
         
       
     
         [0048]    The thus obtained new correlation matrix Rus is utilized to estimate information about each target  200 . A high-resolution estimation method, such as the MUSIC method, the ESPIRIT method and the Capon method, may be desirably employed for estimation.