Abstract:
An electronic scan type radar device which uses a high resolution performance processing to estimate directions of arrival of radio waves, wherein powers of arrival waves received for targets are accurately calculated, that is, a vehicle-mounted radar device utilizing electronic scan which uses a predetermined angle estimation system to estimate directions of arrival of reflected waves, comprising finding mode vectors for angles calculated from the receive signals of the antennas, decomposing a vector of the receive signals into directions of the mode vectors, and defining the lengths of the decomposed vectors the receive powers of the reflected waves arriving from the targets. Due to this method, even if there are targets, it is possible to accurately calculate the powers of the arrival waves, whereby pairing is accurately performed, the precision of detection of targets is improved, and erroneous operation of a vehicle-mounted radar device utilizing electronic scans is prevented.

Description:
CROSS-REFERENCE TO RELATED APPLICATION 
       [0001]    This application claims priority from, and incorporates by reference the entire disclosure of, Japanese Patent Application No. 2011-020058 filed on Feb. 1, 2011. 
       BACKGROUND OF THE INVENTION 
       [0002]    1. Field of the Invention 
         [0003]    The present invention relates to a radar device and a method of calculation of receive power in that radar device, more particularly relates to a radar device which receives radio waves sent from a vehicle and reflected back from a target at a plurality of receive antennas to detect a position of the target and a method of calculation of receive power in that radar device. 
         [0004]    2. Description of the Related Art 
         [0005]    In the past, there have been radar devices which continuously measure the distance and direction of preceding vehicles or targets in the front (obstacles) or of vehicles or other targets approaching from the rear relative to a vehicle being driven by a user and which prevent collisions or enable automatic driving. Such a radar device sent out transmit waves from an antenna installed in the vehicle, received reflected waves bounced back after hitting a target at another antenna, processed the receive signals, and estimated the direction of arrival the reflected waves to detect the target. As the method for estimating the direction of arrival of reflected waves, the DBF method, Capon method, LP method, minimum norm method, MUSIC method, ESPRIT method, and PRISM method are known. 
         [0006]    DBF: Digital Beam Forming 
         [0007]    LP: Linear Prediction 
         [0008]    MUSIC: Multiple Signal Classification 
         [0009]    ESPRIT: Estimation of Signal Parameters via Rotational Invariance Techniques 
         [0010]    PRISM: Propagator method based on an Improved Spatial-Smoothing Matrix 
         [0011]    In the method of estimation of the direction of arrival of reflected waves, the device sends out transmit waves consisting of radio waves modulated in frequency up and down, receives waves reflected from a target by electronic scan radar using a plurality of receive antennas, estimates the direction of arrival of the target at each peak from the peak frequency signals at the up time and down time of the frequency, and then uses the powers as the basis for pairing to calculate the distance and relative speed. The electronic scan radar described in Japanese Patent No. 4415040 found the vectors of the frequency spectrums which were detected for the antennas in the up section of two frequencies, found the combined value when making the vectors of the spectrums of the different frequencies the same angle, and found the distance etc. of a target from the combined value of the vectors. In high resolution processing, as representative systems comparing the spectrum values with the power, there are the DBF method and the Capon method. For this reason, when finding the power of a direction of arrival, the DBF method or the Capon method has mostly been used. 
         [0012]    Further, as the method of calculation of the power of a target, up until now the Capon method has mostly been used. The reason is that in the DBF method, other arrival waves are received at the side lobes, so sometimes the power of the reflected waves could not be accurately calculated. As opposed to this, in the Capon method, by turning the main lobe toward a certain direction and simultaneously minimizing the contribution from other directions to the output, the power of the reflected waves is found more accurately than the DBF method. That is, compared with the DBF method, in the Capon method, the spectrum values are output proportionally by the power of the receive waves, so the power of the reflected waves is accurately calculated. Note that, when estimating the direction of arrival of the reflected waves, usually the spatial average method is used for removing the cross correlation components by preprocessing of angle estimation. 
         [0013]    In this regard, when using the Capon method to calculate the power of the reflected waves, if there are a plurality of targets, reflected waves from the plurality of targets end up being received by the antennas and the cross correlation components of the targets are not removed at the time of preprocessing of angle estimation using the spatial average method, so there is the problem that the receive power cannot be accurately calculated. This becomes remarkable when there are a small number of receive antennas or when there are a small number of snapshots. 
         [0014]    Further, in an electronic scan radar device using a plurality of receive antennas, the receive power value of the arrival waves is used conditional on pairing of the targets obtained in the up section of the frequency and targets obtained in the down section of the frequency. Therefore, if calculation of the accurate receive power is not possible, the following problems arise. 
         [0015]    (1) The above condition of pairing includes the condition of “close power values of angles”. This is no longer satisfied and the possibility of mispairing occurring rises. This ends up causing the electronic scan radar device to malfunction. 
         [0016]    (2) In composite processing of high resolution processing and a system which requires precision of the power such as an amplitude monopulse, it is not possible to improve the precision of the amplitude monopulse. 
       SUMMARY OF THE INVENTION 
       [0017]    Therefore, the present invention has as its object to solve the problems in the related art and provide a radar device which utilizes electronic scan for estimating the direction of arrival of reflected waves, wherein when there are a plurality of targets, reflected waves from the plurality of targets end up being received by the antennas, and the cross correlation components of the targets cannot be removed at the time of preprocessing of angle estimation using the spatial average method, it is possible to accurately calculate the receive powers of the plurality of reflected waves, for example, it is possible to prevent the occurrence of mispairing, and provide a method of calculation of receive power in such a radar device. 
         [0018]    To achieve the above object, the radar device of the present invention provides a radar device which sends out radio waves, receives radio waves reflected from at least one target at a plurality of antennas, and has a signal processing unit process the receive signals by a predetermined angle estimation system, estimate directions of arrival of the reflected waves, and calculate receive powers of the reflected waves to detect the targets, which signal processing unit is provided with an angle detecting means for detecting angles based on the receive signals of the antennas, a mode vector calculating means for finding mode vectors for detected angles, a receive signal decomposing means for using the mode vectors as the basis to decompose a vector which combines the receive signals into vectors of the directions of the mode vectors, and a power value calculating means for calculating the magnitudes of the vectors of the mode vector directions of the decomposed receive signal as power values of the detected angles. 
         [0019]    Further, a method of calculation of receive power in a radar device of the present invention for achieving the above object is a method of calculation of receive power in a radar device comprising sending out radio waves, receiving radio waves which are reflected from at least one target at a plurality of antennas, processing the receive signals using a predetermined angle estimation system, estimating the directions of arrival of the reflected waves, and calculating receive powers of the reflected waves to detect the targets, characterized by detecting angles from receive signals of the antennas, finding mode vectors for the detected angles, using these mode vectors as the basis to decompose a vector of the receive signals into vectors of directions of the mode vectors, and defining the magnitudes of the vectors of the mode vector directions of the decomposed receive signals as the power values of the detected angles. 
         [0020]    According to the vehicle-mounted radar device of the present application and the method of calculation of receive power in such a radar device, even if there are a plurality of targets, it is possible to accurately calculate the receive powers of a plurality of antennas and, for example, possible to accurately perform pairing. As a result, there are the advantageous effects that the precision of detection of targets is improved and erroneous operation of a vehicle-mounted radar device utilizing electronic scans is prevented. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0021]    The present invention is illustrated by way of example, and not limitation, in the figures of the accompanying drawings in which like references indicate similar elements. Note that the following figures are not necessarily drawn to scale. 
           [0022]      FIG. 1  is a view of the configuration which shows the configuration of an FMCW type radar device according to an embodiment of the present invention. 
           [0023]      FIG. 2A  is a waveform diagram of a transmit signal which is emitted from a sending unit of the radar device which is shown in  FIG. 1 . 
           [0024]      FIG. 2B  is a waveform diagram which shows a change in frequency with respect to time of a reflected signal which is bounded back from a target and beat signals in a mixer. 
           [0025]      FIG. 2C  is a waveform diagram of a frequency spectrum which shows the results of FFT processing of an UP beat which is output from an individual receiving unit of one receive antenna which is shown in  FIG. 1 . 
           [0026]      FIG. 2D  is a waveform diagram of a frequency spectrum which shows the results of FFT processing of a DOWN beat which is output from an individual receiving unit of one receive antenna which is shown in  FIG. 1 . 
           [0027]      FIG. 3  is a graph which shows the magnitude of a spectrum with respect to angle in the case of receiving reflected waves which are bounded back from a target by a plurality of antennas which are arranged at equal intervals and of processing this by a conventional angle estimation system. 
           [0028]      FIG. 4  is a view which explains a state of pairing of peak information at an UP beat side and peak information of a DOWN beat side in a bearing computation unit of  FIG. 1 . 
           [0029]      FIG. 5  is a view showing phase differences and mode vectors of antennas in the case of receiving reflected waves which are bounced back from a target by a plurality of antennas which are arranged at equal intervals. 
           [0030]      FIG. 6  is a block diagram which shows one example of a receive power calculation unit which is provided at the inside of the bearing computation unit of  FIG. 1 . 
           [0031]      FIG. 7  is a view decomposing a receive signal in the mode vector directions based on the mode vectors for the angles found from the graph of  FIG. 3 . 
           [0032]      FIG. 8A  is a view of the configuration of an evaluation system in which a corner reflector is arranged as a target at a predetermined distance in a 0 [deg] direction of the radar device. 
           [0033]      FIG. 8B  is a distribution chart which shows the distribution of values of angle powers in power values which are found by the conventional Capon method from receive signals when sending radio waves 100 times from the evaluation system of  FIG. 8A . 
           [0034]      FIG. 8C  is a view of the configuration of an evaluation system in which a corner reflector is arranged as a target at a predetermined distance in a −15 [deg] direction of the radar device. 
           [0035]      FIG. 8D  is a distribution chart which shows the distribution of values of angle powers in power values which are found by the conventional Capon method from receive signals when sending radio waves 100 times from the evaluation system of  FIG. 8C . 
           [0036]      FIG. 9A  is a distribution chart which shows the distribution of values of angle powers in power values which are found by the method of the present invention from receive signals when sending radio waves 100 times from the evaluation system of  FIG. 8A . 
           [0037]      FIG. 9B  is a distribution chart which shows the distribution of values of angle powers in power values of the mode vector directions at positions which were calculated by the device of the present invention from receive signals when sending radio waves 100 times from the evaluation system of  FIG. 8C . 
           [0038]      FIG. 10A  is a view showing an evaluation system in which one target is set. 
           [0039]      FIG. 10B  is a view showing an evaluation system in which two targets are set. 
           [0040]      FIG. 10C  is a view showing an evaluation system in which three targets are set. 
           [0041]      FIG. 10D  is a view showing an evaluation system in which four targets are set. 
           [0042]      FIG. 10E  is a view showing an evaluation system in which five targets are set. 
           [0043]      FIG. 11  is a distribution chart which shows the distribution of values of angle powers in power values which are found by the present invention from receive signals when, as shown in  FIG. 10A  to  FIG. 10E , increasing the number of targets placed in front of the evaluation system one at a time from one to five and sending radio waves 100 times in the evaluation systems. 
       
    
    
     DESCRIPTION OF THE PREFERRED EMBODIMENTS 
       [0044]    Below, the attached drawings will be used to explain in detail the aspects of the present application based on specific embodiments. 
         [0045]      FIG. 1  shows the configuration of a radar device  100  of one embodiment of the present invention. The radar device  100  of this embodiment is comprised of a sending unit S, a receiving unit R, and a signal processing system P. The signal processing system P, while illustration of the detailed configuration is omitted, is configured provided with a microprocessor. There are also a Fourier transform unit  9 , a distance/relative speed processing unit  10 , and a sending/receiving control unit  20 . 
         [0046]    The sending unit S is provided with an oscillator  5  and a signal generation unit  25 . The signal generation unit  25  is controlled by a sending/receiving control unit  20  located at the signal processing system P. The signal generation unit  25  supplies a triangular wave-shaped modulated signal (triangular waves) as the transmit signal to the oscillator  5  where frequency modulation is performed. A transmit antenna  1  transmits radio waves (transmit waves) W. In this embodiment, the FMCW system is used. The oscillator  5  generates transmit waves which change by a certain repeating frequency due to the triangular waves of the signal generation unit  25 . Therefore, the transmit waves W are FMCW waves where the frequency rises and falls about an emitted frequency at the time of no modulation by the oscillator  5  by a predetermined repeated period. This transmit waves W are amplified in power by a not shown transmitter then are sent (emitted) from the transmit antenna  1  toward targets. 
         [0047]    The radar device  100  of this embodiment is mounted in a car. The transmit waves W are sent toward the front or the rear of the vehicle which mounts the radar device  100 . The transmit waves W which are sent from the transmit antenna  1  to the front are reflected at a not shown target, for example, a preceding vehicle or stationary object etc. The reflected waves RW return toward the vehicle and are received by the receiving unit R of the radar device  100 . 
         [0048]    The receiving unit R is comprised of an array antenna  3  which is provided with n number of receive antennas A 1  to An and individual receiving units R 1  to Rn which are connected to the same. The individual receiving units R 1  to Rn respectively have mixers M 1  to Mn and A/D converters (described as A/D in the figure) C 1  to Cn. The receive signals which are obtained from the reflected waves RW 1  to RWn which are received at the array antenna  3  are amplified by a not shown low noise amplifier, then sent to the mixers M 1  to Mn. The mixers M 1  to Mn receive as input the transmit signal from the oscillator  5  of the sending unit S. At the mixers M 1  to Mn, the transmit signal and the receive signals are respectively mixed whereby beat signals having as frequencies the differences of the frequency of the transmit signal and the frequencies of the receive signals are obtained. The beat signals from the mixers M 1  to Mn are converted at the A/D converters C 1  to Cn to digital receive signals X 1  to Xn, then are supplied to a fast Fourier transform device of the Fourier transform unit  9  where each of the digital receive signals X 1  to Xn is analyzed for frequency by a fast Fourier transform (FFT processing). 
         [0049]    In the radar device  100  of this embodiment, when the preceding targets move together, the frequencies of the reflected waves RW include Doppler frequency components which are proportional to the relative speeds of the targets and the vehicle. Further, in the present embodiment, FMCW is employed as the modulation scheme, so when this frequency trend is a linear chirp, the frequencies of the reflected waves RW include, in addition to Doppler components, frequency components which reflect the delay times which are added due to the transmit waves being propagated over the relative distances between the targets and the vehicle. As explained above, the transmit signal is a linear chirp signal, so the frequency of the transmitted waves W, as shown by the solid line in the waveform diagram of  FIG. 2A , repeats a period in which the frequency rises in a straight line (up section) and a period in which it descends (down section). Further, the reflected waves RW, as shown by the broken line in the waveform diagram of  FIG. 2A , compared with the transmit waves W, are simultaneously affected by both the Doppler frequency trend due to the relative speed and the time delay due to the relative distance, so the difference in frequencies between the transmit wave W and the reflected waves RW generally takes different values in the up section and the down section. 
         [0050]    That is, the frequency of the difference of frequencies between the transmit wave W and the reflected waves RW becomes fup in the up section and fdown in the down section. Therefore, at the mixers M 1  to Mn, a beat signal shown in the waveform diagram of  FIG. 2B  in which the Doppler frequency is superposed on the frequency based on the delay time is obtained. The beat signal in the up section is called the “UP beat”, while the beat signal in the down section is called the “DOWN beat”. Note that, in the case of  FIG. 2A  and  FIG. 2B , the frequency fdown of the DOWN beat is larger than the frequency fup of the UP beat. The relative speed of the direction of a smaller relative distance from the target (approach direction) is shown. 
         [0051]    The beat signals of the UP beat and the DOWN beat which are obtained at the mixers M 1  to Mn, as explained above, are converted at the A/D converters C 1  to Cn to the digital receive signals X 1  to Xn which are then supplied to the Fourier transform unit  9 . At the Fourier transform unit  9 , the UP beat frequency fup components and the DOWN beat frequency fdown components from the mixers M 1  to Mn are supplied to the fast Fourier transform device where a fast Fourier transform is applied for frequency analysis (FFT processing). Here, the results of the FFT processing of the receive antenna A 1  are shown in  FIG. 2C  and  FIG. 2D . The waveform diagram of  FIG. 2C  illustrates the frequency spectrum which is obtained from the UP beat frequency fup component, while the waveform diagram of  FIG. 2D  illustrates the frequency spectrum which is obtained from the DOWN beat frequency fdown component. 
         [0052]    As illustrated in  FIG. 2C , in the frequency spectrum of the FFT results of the UP beat of the antenna A 1 , there are peaks Pu 11 , Pu 12 , and Pu 13  at the UP frequencies fu 1 , fu 2 , and fu 3 . For the receive antennas A 2  to An as well, similar FFT results having the same peak frequencies are obtained. For example, at the antenna A 2 , FFT results are obtained with peaks Pu 21 , Pu 22 , and Pu 23  at the UP frequencies fu 1 , fu 2 , and fu 3 . Further, at the frequency spectrum of the FTT results of the DOWN beat of the antenna A 1 , there are peaks Pd 11  and Pd 12  at the DOWN frequencies fd 1  and fd 2 . For the receive antennas A 2  to An as well, similar FFT results having the same peak frequencies are obtained. For example, at the antenna A 2 , FFT results are obtained with peaks Pd 21  and Pd 22  at the DOWN frequencies fd 1  and fd 2 . 
         [0053]    That is, the receive antennas A 1  to An receive reflected radio waves RW from the same targets, so frequency spectrums of the same shapes having the same peak frequencies are obtained by FFT processing. However, the phases of the reflected waves differ depending on the receive antennas, so the phase information of the peaks of the same frequencies will differ at the different receive antennas. 
         [0054]    Returning to  FIG. 1 , the output of the Fourier transform unit  9  is supplied to a peak extraction unit  13  of the distance/relative speed processing unit  10 . At the peak extraction unit  13 , for each of the receive antennas A 1  to An, at the frequency spectrum obtained by the FFT processing, the peaks of a predetermined power or more are extracted at the UP beat and DOWN beat and the frequency, power, and phase information of the extracted peaks (below, called the “peak frequency information”) are extracted. The peak frequency information which was extracted at the peak extraction unit  13  is supplied to a bearing computation unit  15 . 
         [0055]    One peak in the frequency spectrum usually includes information of a plurality of targets, so it is necessary to separate the targets from the single peak and estimate the angles of the separated targets. For this reason, at the bearing computation unit  15 , the peak frequency information of the peaks having the same frequencies at the UP side and DOWN side at all receive antennas A 1  to An (for example, in the case of an UP beat, Pu 11 , Pu 21 , . . . Pun 1 , while in the case of a DOWN beat, Pd 11 , Pd 21 , . . . Pdn 1 ) are used as the basis to find the angle spectrums such as shown in  FIG. 3  by computation. As the method of finding the angle spectrums, the Capon method, DBF method, or other systems may be used. The solid line in  FIG. 3  illustrates the angle spectrum of the UP peak frequency fu 1  (Pu 11 , Pu 21 , . . . Pun 1 ), while the broken line illustrates the angle spectrum of the DOWN peak frequency fd 1  (Pd 11 , Pd 21 , . . . Pdn 1 ). 
         [0056]    At the bearing computation unit  15 , in the angle spectrums which are shown in  FIG. 3 , the peaks having powers of a threshold value or more, here, the peaks P 1  and P 2 , are judged to be targets and their angles and powers are extracted. As explained in further detail, the angle spectrum is found for each peak frequency in FFT processing. In the examples shown in  FIG. 2C  and  FIG. 2D , five angle spectrums at five frequencies fu 1 , fu 2 , fu 3 , fd 1 , and fd 2  are calculated.  FIG. 3  shows the angle spectrum which is found from the peaks of the UP peak frequency fu 1  and the angle spectrum which is found from the peaks of the DOWN peak frequency fd 1 . This shows that there are two targets P 1  (angle 0 [Deg]) and P 2  (angle about 3 [Deg]) at both the UP peak frequency fu 1  and at the DOWN peak frequency fd 1 . The results which are obtained at the bearing computation unit  15  become as illustrated in  FIG. 4 . 
         [0057]    The distance/relative speed computation unit  30  uses the data which is shown in  FIG. 4  as the basis to pair data having close angles and powers between the UP beat side target information and the DOWN beat side target information. In  FIG. 4 , the target U 1  of the angle θu 1  of the UP beat side frequency fu 1  and the target D 2  of the angle θd 2  of the DOWN beat side frequency fd 1  are shown paired. Detection of five targets is shown. The distances and relative speeds are calculated at the UP frequency and DOWN frequency obtained by the pairing. For the angles of the targets, the average values of the UP beat side and DOWN beat side angles are obtained. The distance and relative speed are found from the UP peak frequency fu 1  and the DOWN peak frequency fd 1 , while the angle is found by (θu 1 +θd 2 )/2. 
         [0058]    Here, as illustrated in  FIG. 5 , the receive antennas A 1  to An comprise the six antennas A 1  to A 6 . The operation of the bearing computation unit  15  in the case of assuming the arrival of only one radio wave at the antennas A 1  to A 6  will be explained. Note that, the distance between adjacent antennas is defined as “d”, the direction of arrival of the arrival wave with respect to a direction vertical to the line connecting the six antennas A 1  to A 6  is defined as “θ”, and the wavelength of arrival waves is defined as “λ”. In this case, the phase difference φ the between adjacent antennas becomes φ=(2π/λ)d sin(θ). Therefore, if the amplitude of the arrival waves at a certain point of time at the first antenna A 1  is A(t), the amplitude of the arrival waves at the second antenna A 2  at the same point of time becomes A(t)exp[j(2π/λ)d sin(θ)]. 
         [0059]    To facilitate the explanation, an array of ideal signals of the different antennas of the signals of the amplitude  1  which arrive from the direction θ is defined as a “mode vector a(θ)”. Further, if considering that the reference point of the equiphase plane at the timing t 1  is at the antenna A 1 , the phases at the antennas A 2  to A 6  with respect to the antenna A 1  at the same timings become as follows: 
         [0060]    antenna A 2 : exp[−j(2π/λ)d sin(θ)] 
         [0061]    antenna A 3 : exp[−j(2π/λ)2d sin(θ)] 
         [0062]    antenna A 4 : exp[−j(2π/λ)3d sin(θ)] 
         [0063]    antenna A 5 : exp[−j(2π/λ)4d sin(θ)] 
         [0064]    antenna A 6 : exp[−j(2π/λ)5d sin(θ)] 
         [0065]    Accordingly, the mode vector a(θ) at this time becomes a(θ)=(1, exp[−j(2π/λ)d sin(θ)], exp[−j(2π/λ)2d sin(θ)], exp[−j(2π/λ)3d sin(θ)], exp [−j(2π/λ))4d sin(θ)], and exp[−j(2π/λ)5d sin(θ)]) t . 
         [0066]    At the bearing computation unit  15 , the angle estimation system Capon method and spatial average method were used to calculate the powers of the directions of arrival of the arrival waves. However, as explained above, with the conventional method of calculation of power of the Capon method, the method was adopted of calculating the power after obtaining the spatial average, so if there were a plurality of targets, cross correlation components of the targets could not be removed at the time of the processing by the spatial average method and the receive powers of the arrival waves could not be accurately calculated. 
         [0067]    Therefore, in the present invention, the receive signals at the stage before the spatial average is taken are used as they are to calculate the receive powers. This device and method will be explained using  FIG. 6  and  FIG. 7 .  FIG. 6  shows an example of the configuration of a receive power calculation unit  50 P which is provided inside of the bearing computation unit  15  of  FIG. 1 . The receive power calculation unit  50 P to which the peak frequency information which is extracted by the peak extraction unit  13  (frequency, power, and phase information of peaks) is input includes an angle detection unit  51 , mode vector calculation unit  52 , receive signal decomposition unit  53 , and power value calculation unit  54 . 
         [0068]    The angle detection unit  51 , as explained using  FIG. 3  and  FIG. 4 , extracts the angles of the peaks of the receive signals at the receive antennas which are input from the peak extraction unit  13 . At this time, if there are a plurality of targets, a plurality of angles are detected from the receive signals. The mode vector calculation unit  52  finds the mode vectors for the detected angles. Further, the mode vectors are found as theoretical values or measured values for predetermined angles (for example, each 1 [Deg]) and are stored in a not shown memory. By reading out a mode vector corresponding to an angle which is detected by the angle detection unit  51 , the mode vector for the detected angle can be found. Further, the receive signal decomposition unit  53  uses the mode vectors as the basis to decompose a vector obtained by combining the peak powers after FFT processing of the receive signals used for the above angle detection into vectors of the directions of the mode vectors and directions of orthogonal components of the same. Note that, at this time, if there are a plurality of peaks after FFT processing of the receive signals, this is performed on the vector combining those peak powers. Further, such vector decomposition is performed for each of the UP section and DOWN section of the triangular waves. The power value calculation unit  54  outputs the magnitudes of the mode vector directions of the decomposed receive signals as the power values of the receive signals. 
         [0069]    If using a predetermined angle estimation system, for example, the conventional PRISM method, to find the magnitude of the spectrum with respect to the angles of receive signals, the result becomes as illustrated in  FIG. 3 . As explained above, the solid line in  FIG. 3  shows the angle spectrum of the UP peak frequency fu 1  (Pu 11 , Pu 21 , . . . Pun 1 ), while the broken line shows the angle spectrum of the DOWN peak frequency fd 1  (Pd 11 , Pd 21 , . . . Pdn 1 ). Further, the abscissa of  FIG. 3  shows the left-right deflection angle [deg] of the transmit wave when defining the front surface of the radar device as 0 [deg]. In this embodiment, corner reflectors (explained later) are assumed to be placed as targets at the 0 [deg] position and the 3 [deg] position. In the present invention, the magnitude of the spectrum showing the power values in  FIG. 3  is ignored. Only the angles corresponding to the peak parts P 1 , P 2 , and P 3  of the magnitude of the spectrum with the circle marks in  FIG. 3  are found by calculation. However, the peak P 3  is not a peak due to a target, but is used for calculation of the power. 
         [0070]    Next, the mode vectors corresponding to the found angles are found, the receive signals are deemed to be a vector of a complex number, the peak power after FFT processing which was used for angle calculation is decomposed into the mode vector directions, and the coefficients at that time are defined as the power values. The specific method of calculation of the receive power will be explained for the case where there are four receive antennas. This calculation method multiplies the receive signals of the basis transition matrix (standard basis is mode vector basis) to find the coefficients of the mode vector bases and finds the norms of the same. Here, the method of finding the orthogonal components will also be explained. 
         [0071]    First, the method defines θ 1  to θ 3  as the angles found by the angle estimation system and codes the mode vectors as follows: 
         [0000]        a (θ1)=( a 11,  a 12,  a 13,  a 14)
 
         [0000]        a (θ2)=( a 21,  a 22,  a 23,  a 24)
 
         [0000]        a (θ3)=( a 31,  a 32,  a 33,  a 34)
 
         [0072]    Here, it sets “A” as follows: 
         [0000]    
       
         
           
             A 
             = 
             
               
                 
                   
                     a 
                      
                     
                         
                     
                      
                     11 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     21 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     31 
                   
                 
               
               
                 
                   
                     a 
                      
                     
                         
                     
                      
                     12 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     22 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     32 
                   
                 
               
               
                 
                   
                     a 
                      
                     
                         
                     
                      
                     13 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     23 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     33 
                   
                 
               
               
                 
                   
                     a 
                      
                     
                         
                     
                      
                     14 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     24 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     34 
                   
                 
               
             
           
         
       
     
         [0073]    Next, it calculates the vectors P which are orthogonal to three mode vectors using the propagator method while defining A 1  and A 2  as follows: 
         [0000]    
       
         
           
             
               
                 A 
                  
                 
                     
                 
                  
                 1 
               
               = 
               
                 
                   
                     
                       a 
                        
                       
                           
                       
                        
                       11 
                     
                   
                   
                     
                       a 
                        
                       
                           
                       
                        
                       21 
                     
                   
                   
                     
                       a 
                        
                       
                           
                       
                        
                       31 
                     
                   
                 
                 
                   
                     
                       a 
                        
                       
                           
                       
                        
                       12 
                     
                   
                   
                     
                       a 
                        
                       
                           
                       
                        
                       22 
                     
                   
                   
                     
                       a 
                        
                       
                           
                       
                        
                       32 
                     
                   
                 
                 
                   
                     
                       a 
                        
                       
                           
                       
                        
                       13 
                     
                   
                   
                     
                       a 
                        
                       
                           
                       
                        
                       23 
                     
                   
                   
                     
                       a 
                        
                       
                           
                       
                        
                       33 
                     
                   
                 
               
             
             _ 
           
         
       
       
         
           
             
               A 
                
               
                   
               
                
               2 
             
             = 
             
               
                 
                   
                     a 
                      
                     
                         
                     
                      
                     14 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     24 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     34 
                   
                 
               
             
           
         
       
       
         
           
             
               
                 
                   P 
                   = 
                   
                     ( 
                     
                       
                         
                           
                             
                               ( 
                               
                                 A 
                                  
                                 
                                     
                                 
                                  
                                 2 
                                 * 
                                 A 
                                  
                                 
                                     
                                 
                                  
                                 
                                   1 
                                   
                                     - 
                                     1 
                                   
                                 
                               
                               ) 
                             
                             H 
                           
                         
                       
                       
                         
                           
                             - 
                             1 
                           
                         
                       
                     
                     ) 
                   
                 
               
             
             
               
                 
                   = 
                   
                     
                       ( 
                       
                         
                           p 
                            
                           
                               
                           
                            
                           1 
                         
                         , 
                         
                           p 
                            
                           
                               
                           
                            
                           2 
                         
                         , 
                         
                           p 
                            
                           
                               
                           
                            
                           3 
                         
                         , 
                         
                           p 
                            
                           
                               
                           
                            
                           4 
                         
                       
                       ) 
                     
                     t 
                   
                 
               
             
           
         
       
     
         [0074]    Here, it determines the matrix comprised of the mode vectors plus the orthogonal components as follows: 
         [0000]    
       
         
           
             B 
             = 
             
               
                 
                   
                     a 
                      
                     
                         
                     
                      
                     11 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     21 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     31 
                   
                 
                 
                   
                     P 
                      
                     
                         
                     
                      
                     1 
                   
                 
               
               
                 
                   
                     a 
                      
                     
                         
                     
                      
                     12 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     22 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     32 
                   
                 
                 
                   
                     P 
                      
                     
                         
                     
                      
                     2 
                   
                 
               
               
                 
                   
                     a 
                      
                     
                         
                     
                      
                     13 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     23 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     33 
                   
                 
                 
                   
                     P 
                      
                     
                         
                     
                      
                     3 
                   
                 
               
               
                 
                   
                     a 
                      
                     
                         
                     
                      
                     14 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     24 
                   
                 
                 
                   
                     a 
                      
                     
                         
                     
                      
                     34 
                   
                 
                 
                   
                     P 
                      
                     
                         
                     
                      
                     4 
                   
                 
               
             
           
         
       
     
         [0075]    Here, it defines the 1ch to 4ch receive signals as Y 1  to Y 4  and converts the receive signals to vectors as follows: 
         [0000]    
       
         
           
             Y 
             = 
             
               
                 
                   
                     Y 
                      
                     
                         
                     
                      
                     1 
                   
                 
               
               
                 
                   
                     Y 
                      
                     
                         
                     
                      
                     2 
                   
                 
               
               
                 
                   
                     Y 
                      
                     
                         
                     
                      
                     3 
                   
                 
               
               
                 
                   
                     Y 
                      
                     
                         
                     
                      
                     4 
                   
                 
               
             
           
         
       
     
         [0076]    Further, if using B and Y for the following calculations, the coefficients of the bases of the mode vectors and orthogonal components are found. 
         [0000]    
       
         
           
             
               
                 B 
                 
                   - 
                   1 
                 
               
                
               Y 
             
             = 
             
               
                 
                   
                     
                       bn 
                        
                       
                           
                       
                        
                       1 
                        
                       
                         ( 
                         
                           = 
                           
                             coefficient 
                              
                             
                                 
                             
                              
                             of 
                              
                             
                                 
                             
                              
                             θ 
                              
                             
                                 
                             
                              
                             1 
                              
                             
                                 
                             
                              
                             bearing 
                           
                         
                         ) 
                       
                     
                   
                 
                 
                   
                     
                       bn 
                        
                       
                           
                       
                        
                       2 
                        
                       
                         ( 
                         
                           = 
                           
                             coefficient 
                              
                             
                                 
                             
                              
                             of 
                              
                             
                                 
                             
                              
                             θ 
                              
                             
                                 
                             
                              
                             2 
                              
                             
                                 
                             
                              
                             bearing 
                           
                         
                         ) 
                       
                     
                   
                 
                 
                   
                     
                       bn 
                        
                       
                           
                       
                        
                       3 
                        
                       
                         ( 
                         
                           = 
                           
                             coefficient 
                              
                             
                                 
                             
                              
                             of 
                              
                             
                                 
                             
                              
                             θ 
                              
                             
                                 
                             
                              
                             3 
                              
                             
                                 
                             
                              
                             bearing 
                           
                         
                         ) 
                       
                     
                   
                 
                 
                   
                     
                       bn 
                        
                       
                           
                       
                        
                       4 
                        
                       
                         ( 
                         
                           = 
                           
                             coefficient 
                              
                             
                                 
                             
                              
                             of 
                              
                             
                                 
                             
                              
                             orthogonal 
                              
                             
                                 
                             
                              
                             components 
                           
                         
                         ) 
                       
                     
                   
                 
               
                
               
                 
                   
                     
                       
                         
                           
                             
                               
                                 
                                     
                                 
                               
                             
                             
                               
                                 
                                   ( 
                                   
                                     = 
                                     BN 
                                   
                                   ) 
                                 
                               
                             
                           
                         
                       
                       
                         
                           
                               
                           
                         
                       
                     
                   
                 
                 
                   
                     
                         
                     
                   
                 
               
             
           
         
       
     
         [0000]    p The norms of the bn 1  to bn 4  found here become the power values of the corresponding bearings of the arrival waves and the power values of the orthogonal components. 
         [0077]    On the other hand, the specific method of calculation in the case where there are no orthogonal components of the receive powers becomes as follows. However, this is the case where there are 4 ch receive antennas. 
         [0078]    When there are orthogonal components, the following formula stands assuming BN is the value sought. 
         [0000]    
       
      
       Y=B*BN  
      
     
         [0079]    Here, B is a regular matrix, so the following formula may be calculated. 
         [0000]    
       
      
       BN=B 
       −1 
       *Y  
      
     
         [0080]    In this regard, where no orthogonal components are included, the following formula stands: 
         [0000]    
       
      
       Y=A*BN  
      
     
         [0081]    “A” is not a square matrix, so cannot have an inverse matrix. For this reason, least square approximation is used to find the BN. This method first defines the following evaluation function: 
         [0000]        Q ( BN )=∥ Y−A*BN∥ (where, ∥·∥ is a vector norm)  
 
         [0082]    Here, it is sufficient to find BN so that Q(BN) becomes 0. Here, if developing Q(BN) further, it becomes as follows where “tr” is the trace of the matrix and “H” is the adjoint matrix. 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       Q 
                        
                       
                         ( 
                         BN 
                         ) 
                       
                     
                     = 
                       
                      
                     
                        
                       
                         Y 
                         - 
                         
                           A 
                           * 
                           BN 
                         
                       
                        
                     
                   
                 
               
               
                 
                   
                     = 
                       
                      
                     
                       tr 
                       ( 
                       
                         
                           
                             ( 
                             
                               Y 
                               - 
                               
                                 A 
                                 * 
                                 BN 
                               
                             
                             ) 
                           
                           H 
                         
                         * 
                         
                           ( 
                           
                             Y 
                             - 
                             
                               A 
                               * 
                               BN 
                             
                           
                           ) 
                         
                       
                     
                   
                 
               
               
                 
                   
                     = 
                       
                      
                     
                       
                         tr 
                          
                         
                           ( 
                           
                             
                               Y 
                               H 
                             
                             * 
                             Y 
                           
                           ) 
                         
                       
                       - 
                       
                         tr 
                          
                         
                           ( 
                           
                             
                               Y 
                               H 
                             
                             * 
                             A 
                             * 
                             BN 
                           
                           ) 
                         
                       
                       - 
                       
                         tr 
                          
                         
                           ( 
                           
                             
                               BN 
                               H 
                             
                             * 
                             
                               A 
                               H 
                             
                             * 
                             Y 
                           
                           ) 
                         
                       
                       + 
                     
                   
                 
               
               
                 
                   
                       
                      
                     
                       tr 
                        
                       
                         ( 
                         
                           
                             BN 
                             H 
                           
                           * 
                           
                             A 
                             H 
                           
                           * 
                           A 
                           * 
                           BN 
                         
                         ) 
                       
                     
                   
                 
               
             
               
           
         
       
     
         [0083]    Therefore, it is sufficient that the value of this differentiated by BN H  (−A H *Y+A H *A*BN) become 0. From the above, the following formula may be used to calculate BN. 
         [0000]        BN =( A   H   *A ) −1   *A   H   *Y    
         [0084]      FIG. 7  illustrates the present invention which finds the mode vectors MV 1  and MV 2  corresponding to the angles of the peak parts P 1  and P 2  found from the “angle-spectrum magnitude” characteristic shown in  FIG. 3 , deems the vector RS combining the peak powers of the receive signals of the antennas after FFT processing to be a vector of a complex number, decomposes the vector RS in the directions of the mode vectors MV 1  and MV 2 , and defines the magnitudes of the vectors at that time as the power values. From  FIG. 7 , it is learned that the power value corresponding to the peak part P 1  shown in  FIG. 3  is shown by the vector PV 1 , while the power value corresponding to the peak part P 2  is shown by the vector PV 2 . 
         [0085]    At this time, by calculating not only the mode vectors MV 1  and MV 2 , but also the vector RV in the direction orthogonal to these mode vectors MV 1  and MV 2  and the orthogonal component RE, it is possible to further raise the precision. Further, by viewing the magnitude of the orthogonal component RE, it is also possible to judge if arrival waves of more than the separation limit are being received. In addition, as an additional effect, it is possible to suppress the effects of noise included in the receive signal on the power. 
         [0086]    As explained above, in the past, receive signals were processed by FFT, then some sort of system was used to calculate the bearings of the received waves and the spectrum values were used as is as the power values. When precision of the power values was particularly required, the powers were recalculated from the spectrum values using the DBF method, Capon method, etc. As opposed to this, in the present invention, the bearings are calculated in the same way as the past to obtain the angles, but “vector decomposition of the receive signals” is employed for calculation of the powers. This is the biggest feature. Below, the effect obtained by such a method of calculation of the power of the present invention will be explained. 
         [0087]    Here, for the radar device  100 , two evaluation systems were prepared by setting corner reflectors CR as targets exactly predetermined distances away from the radar device  100  at a position of 0 [deg] such as shown in  FIG. 8A  and a position of −15 [deg] such as shown in  FIG. 8C . Further, data showing the values of the angle powers at power values of receive signals obtained by sending out radio waves 100 times from this radar device  100  and receiving the reflected waves at the antennas was calculated using the conventional Capon method. 
         [0088]      FIG. 8B  shows 100 sets of data for the case of using the evaluation system shown in  FIG. 8A , while  FIG. 8D  shows 100 sets of data for the case of using the evaluation system shown in  FIG. 8C . In  FIG. 8B  and  FIG. 8D , the ordinate shows the values of the angle powers [dB] at the power values of the receive signals, while the abscissa shows 100 receive signals as the “number of data”. The straight line parts shown by the two broken lines in the data shown in  FIG. 8B  and  FIG. 8D  show the true values of the angle powers. As will be understood from the data of  FIG. 8B  and  FIG. 8D , with the conventional system, the fluctuation in the angle powers is extremely great. 
         [0089]    On the other hand, the data showing the angle powers shown in  FIG. 9A  shows the angle powers at the power values which were measured by using the vector decomposition method of the present invention for 100 receive signals in the case of using the evaluation system shown in  FIG. 8A . Further, the data showing the angle powers shown in  FIG. 9B  shows the angle powers at the power values of the receive signals which were measured by using the vector decomposition method of the present invention for 100 receive signals in the case of using the evaluation system shown in  FIG. 8C . In  FIG. 9A  and  FIG. 9B  as well, the straight line parts shown by the two broken lines show the true values of the angle powers due to the reflected waves from a single target. It is learned that by using the method of calculation of power of the present invention, the angle powers at the power values of the receive signals can be detected considerably stably and that the detected values are values equivalent to the straight line parts shown by the two broken lines. 
         [0090]    Here, as shown in  FIG. 10A  to  FIG. 10E , evaluation systems are prepared in which one target ( FIG. 10A ), two targets ( FIG. 10B ), three targets ( FIG. 10C ), four targets ( FIG. 10D ), and five targets ( FIG. 10E ) are arranged. Further, each evaluation system is sent radio waves from the radar device  100  a total of 100 times.  FIG. 11  shows all together the angle powers at the power values of the orthogonal components obtained by using the method of calculation of power of the present invention for receive signals obtained by reception of reflected waves at antennas in the different evaluation systems. The horizontal line shown in bold in the figure is the threshold value R for judging the separable number and is a value similar to the threshold value of FFT. 
         [0091]    As will be understood from  FIG. 11 , when reflected waves from the separable number or less (in embodiment of  FIG. 11 , three or less) of targets are received as arrival waves at the antennas and the angles can be estimated accurately, the angle powers at the power values of the orthogonal components are smaller than the threshold value R. On the other hand, when reflected waves from the separable number or more (in embodiment of  FIG. 11 , four or more) of targets are received as arrival waves at the antennas and the angles cannot be estimated accurately, the angle powers at the power values of the orthogonal components are larger than the threshold value R. 
         [0092]    By viewing the angle powers at the power values of the orthogonal components, it is possible to judge if the separable number or more of arrival waves of received waves have arrived. Note that, when reflected waves from the separable number or more of targets have been received as arrival waves at the antennas, the angles cannot be accurately obtained even by the method of the present invention. Such a judgment can be used for judgment for changing to a separating means with a higher separable number of received waves. 
         [0093]    Although only some exemplary embodiments of this invention have been described in detail above, those skilled in the art will readily appreciate that many modifications are possible in the exemplary embodiments without materially departing from the novel teachings and advantages of this invention. Accordingly, all such modifications are intended to be included within the scope of this invention.