US20090040097A1 - Radar apparatus - Google Patents
Radar apparatus Download PDFInfo
- Publication number
- US20090040097A1 US20090040097A1 US11/497,033 US49703306A US2009040097A1 US 20090040097 A1 US20090040097 A1 US 20090040097A1 US 49703306 A US49703306 A US 49703306A US 2009040097 A1 US2009040097 A1 US 2009040097A1
- Authority
- US
- United States
- Prior art keywords
- previous
- correlation matrices
- addition
- radar apparatus
- correlation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011159 matrix material Substances 0.000 claims abstract description 163
- 238000000605 extraction Methods 0.000 claims abstract description 44
- 238000004364 calculation method Methods 0.000 claims abstract description 29
- 238000001514 detection method Methods 0.000 claims abstract description 8
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims abstract description 7
- 239000000284 extract Substances 0.000 claims abstract description 3
- 238000000034 method Methods 0.000 description 97
- 238000001228 spectrum Methods 0.000 description 36
- 230000007423 decrease Effects 0.000 description 7
- 230000035559 beat frequency Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000010363 phase shift Effects 0.000 description 2
- 238000013144 data compression Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/74—Multi-channel systems specially adapted for direction-finding, i.e. having a single antenna system capable of giving simultaneous indications of the directions of different signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
- G01S13/32—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
- G01S13/34—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
- G01S13/32—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
- G01S13/34—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal
- G01S13/345—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal using triangular modulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/42—Simultaneous measurement of distance and other co-ordinates
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/46—Indirect determination of position data
- G01S13/48—Indirect determination of position data using multiple beams at emission or reception
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
Definitions
- the present invention relates to a radar apparatus.
- a frequency modulated continuous wave (FMCW) radar apparatus is a radar apparatus that continuously transmits a frequency-modulated radar signal to the object and detects the distance or a relative velocity of the object based on the transmitted radar signal reflected from the object.
- FMCW frequency modulated continuous wave
- a transmitting means for transmitting a radar signal mechanically turns and scans the transmitted radar signal reflected from the object.
- the transmitting means is fixed and an antenna having multiple elements arranged in an array receives the transmitted radar signal.
- the received radar signal is digitally processed and the direction of the object is detected based on the digital signal.
- an angular spectrum is generated based on the received radar signal on each of the elements and a peak of the angular spectrum is detected.
- the direction of the object is estimated based on the peak of the angular spectrum.
- a beamformer algorithm as the DBF algorithm, the angular spectrum is generated such that amplitudes of the received signal at a given time are connected as shown in FIGS. 10A and 10B .
- a multiple signal classification (MUSIC) algorithm is known as a high-resolution direction of arrival (DOA) estimation algorithm.
- DOA direction of arrival
- correlation matrices are calculated, eigenvalue expansions are performed on each of the correlation matrices, the angular spectrum is calculated from eigenvectors of the correlation matrices, and the direction of the object is calculated based on the angular spectrum.
- a FMCW radar apparatus disclosed in U.S. Pat. No. 6,121,917 corresponding to JP-A-H11-133142 detects the direction of the object by using the beamformer algorithm.
- a fast Fourier transform (FFT) is applied to the received wave signal to obtain a peak frequency of a distance power spectrum.
- the beamformer algorithm is applied to only the peak frequency component of the received signal so that the amount of calculation required to detect the direction of the object is reduced.
- the radar apparatus using the beamformer algorithm needs to be increased in size to obtain high resolution.
- the high-resolution DOA estimation algorithm such as the MUSIC algorithm achieves the high resolution without an increase in the number of the elements.
- the resolution may be increased by reducing noise with summation of the received signal with respect to time. The summation is performed such that a present angular spectrum calculated in a present process and a previous angular spectrum calculated in a previous process are summed up.
- the DOA estimation algorithm needs to be applied to all the frequency components of the received signal on each process to calculate the angular spectrum, and the calculated angular spectrum needs to be stored in a memory.
- the calculation of the angular spectrum requires an eigenvalue expansion that requires a lot of calculation. Therefore, when the angular spectrum is calculated on each frequency component, the amount of calculation is significantly increased. In the DOA estimation algorithm, the high resolution results in a significant increase in the amount of calculation.
- a radar apparatus uses a DOA estimation algorithm performed at a predetermined time interval to detect an object.
- the radar apparatus includes a transmitting means that transmits a wave signal to the object, a receiving means that receives the transmitted wave signal reflected from the object, a beat signal generation means that generates a plurality of beat signals based on the received wave signal, a correlation matrix generation means that calculates a plurality of correlation matrices based on the beat signals with respect to each frequency component of the beat signals, an addition means that calculates a plurality of addition correlation matrices, a storing means that stores at least a portion of a plurality of previous correlation matrices calculated by the correlation matrix generation means at a previous time, a detection means that detects at least one frequency component of the beat signals, the frequency component satisfying a predetermined condition, an extraction means that extracts an extraction matrix from the addition correlation matrices, the extraction matrix corresponding to a frequency component closest to the frequency component detected by the detection means, and a direction calculation means that calculates a
- the receiving means includes a plurality of elements arranged in an array and each of the beat signals generated by the beat signal generation means corresponds to each of the elements.
- the addition means calculates the addition correlation matrices by adding the correlation matrices calculated by the correlation matrix generation means to the previous correlation matrices stored by the storing means.
- correlation matrices are generated from beat signals. Then, angular spectrums are calculated from each eigenvector of each of the correlation matrices.
- the storing means stores the angular spectrums as previous information and the angular spectrums are used in a next process to reduce noise.
- a position of the object i.e., a frequency component indicating a presence of the object
- an angular spectrum corresponding to the position is calculated.
- Correlation matrices calculated in a present process and previous correlation matrices calculated in a previous process are summed up to reduce noise.
- the storing means stores the correlation matrices as the previous information, not the angular spectrums.
- the amount of calculation required to calculate the correlation matrices is about one-thirteenth the amount of calculation required to calculate MUSIC spectrums. Because the radar apparatus uses the correlation matrices instead of the angular spectrums to reduce the noise, a significant increase in the amount of calculation can be prevented. Thus, the radar apparatus achieves a high resolution without the significant increase in the amount of calculation.
- FIG. 1 is a block diagram showing a radar apparatus according to a first embodiment of the present invention
- FIG. 2 is an arrangement of elements of a receiving antenna of the radar apparatus of FIG. 1 ;
- FIGS. 3A and 3B are diagrams showing a principal of beat signals generated in the radar apparatus of FIG. 1 ;
- FIG. 4A is a diagram showing a reflected wave signal received by the elements of the receiving antenna of the radar apparatus of FIG. 1
- FIG. 4B is FFT beat signals generated by applying a FFT to beat signals of FIG. 3B
- FIG. 4C is a diagram showing a sum beat signal into which the FFT beat signals of FIG. 4B are summed;
- FIG. 5 is a flow chart illustrating a process performed by a microcomputer of the radar apparatus of FIG. 1 ;
- FIG. 6 is a table showing the amount of calculation performed by the microcomputer of the radar apparatus of FIG. 1 ;
- FIG. 7 is a flow chart illustrating a process performed by a microcomputer of a radar apparatus according to a second embodiment of the present invention.
- FIG. 8 is a flow chart illustrating a process performed by a microcomputer of a radar apparatus according to a third embodiment of the present invention.
- FIG. 9 is a flow chart illustrating a process performed by a microcomputer of a radar apparatus according to a fourth embodiment of the present invention.
- FIG. 10A is a diagram showing an angular spectrum generated when reflected wave signals arrive at elements from a front direction in a beam former method
- FIG. 10B is a diagram showing the angular spectrum generated when the reflected wave signals arrive at the elements from an oblique direction in the beam former method.
- the radar apparatus 100 includes a transmitting antenna 11 , a receiving antenna 12 having N elements E 1 -EN arranged in an array, where N is a positive integer, a high-frequency switch 13 , a mixer 14 , an oscillator 15 , a digital-to-analog (D/A) converter 16 , an analog-to-digital (A/D) converter 17 , a microcomputer 18 , a switch controller 19 , and a timer 20 .
- D/A digital-to-analog
- A/D analog-to-digital converter
- the microcomputer 18 outputs a digital transmit signal to the D/A converter 16 .
- the D/A converter 16 converts the digital transmit signal into an analog transmit signal and outputs the analog transmit signal to the oscillator 15 .
- the oscillator 15 outputs the analog transmit signal to the transmitting antenna 11 at a predetermined frequency.
- the transmitting antenna 11 converts the analog transmit signal into a radar signal and transmits the radar signal to an object.
- the receiving antenna 12 receives the transmitted radar signal reflected from the object.
- the received radar signal is converted into an analog receive signal on each of the elements E 1 -EN.
- the high-frequency switch 13 in turn sends the analog receive signal to the mixer 14 .
- the analog receive signal is mixed with the analog transmit signal to generate a beat signal on each of the elements E 1 -EN.
- the beat signal is input to the A/D converter 17 and converted into a digital receive signal.
- the digital signal is input to the microcomputer 18 .
- the microcomputer 18 controls the high-frequency switch 13 through the switch controller 19 and controls a sampling rate of the A/D converter 17 though the timer 20 .
- the microcomputer 18 has a memory 21 .
- the microcomputer 18 performs a distance calculation process for calculating the distance between the radar apparatus 100 and the object and a direction calculation process for calculating the direction between the radar apparatus 100 and the object.
- the elements E 1 -EN of the receiving antenna 12 are spaced from each other by a predetermined spacing S.
- the time delay corresponds to the distance between the radar apparatus 100 and the object
- the frequency shift corresponds to a relative velocity between the radar apparatus 100 and the object.
- a phase shift of the received radar signal with respect to the transmitted radar signal increases with the distance between the radar apparatus 100 and the object.
- the distance and relative velocity between the radar apparatus 100 and the object can be detected based on the phase shift. Therefore, a beat signal is calculated as a frequency difference between the transmitted radar signal and the received radar signal.
- the beat signal has a first beat frequency Bu in an increase area where the frequency of the transmitted radar signal increases, and has a second beat frequency Bd in a decrease area where the frequency of the transmitted radar signal decreases.
- the beat signal includes a first beat signal having the first beat frequency Bu and a second beat signal having the second beat frequency Bd.
- the beat signal is generated on each of the elements E 1 -EN.
- 2N beat signals are generated. Specifically, N first beat signals having the first beat frequency Bu and N second beat signals having the second beat frequency Bd are generated.
- V ⁇ C /(4 ⁇ F 0) ⁇ ( BuH ⁇ BdH ) (2)
- BuH is the first beat signal generated based on the received radar signal that is received by the element EH, where H is a positive integer less than or equal to N (i.e., 1 ⁇ H ⁇ N).
- BdH is the second beat signal generated based on the received radar signal that is received by the element EH.
- the microcomputer 18 performs the direction calculation process at a predetermined time interval Ts.
- a fast Fourier transform is applied to each of the N beat signals Bu 1 -BuN to generate N FFT beat signals Bfu 1 -BfuN.
- the fast Fourier transform is applied to the N beat signals Bd 1 -BdN to generate N FFT beat signals Bfd 1 -BfdN.
- the FFT beat signals Bfd 1 -BfdN are processed in the same way as the FFT beat signals Bfu 1 -BfuN.
- a correlation matrix group RG of correlation matrices are generated by using each of the FFT beat signals.
- the correlation matrix group RG includes M correlation matrices R(F 1 )-R(FM), where M is a positive integer greater than 1.
- a correlation matrix R(FI) corresponds to a frequency FI, where I is a positive integer less than or equal to M (i.e., 1 ⁇ I ⁇ M).
- each of the correlation matrices R(F 1 )-R(FM) is an N ⁇ N matrix.
- an addition correlation matrix group UG having addition correlation matrices U(F 1 )-U(FM) is generated such that the correlation matrix group RG generated in a present process is added to a previous correlation matrix group RoG that is the correlation matrix group RG generated in a previous process (i.e., Ts earlier) and stored in the memory 21 .
- the addition correlation matrix group UG has less noise than the correlation matrix group RG.
- the use of the previous correlation matrix group RoG as previous information reduces the noise.
- a MUSIC algorithm used in the radar apparatus 100 is described below.
- the MUSIC algorithm allows the radar apparatus 100 to achieve the high resolution without an increase in the amount of calculation.
- the FFT beat signals Bfu 1 -BfuN shown in FIG. 4B are summed up to generate a sum beat signal Bfu 0 shown in FIG. 4C .
- the sum beat signal Bfu 0 has less noise than each of the FFT beat signals Bfu 1 -BfuN.
- the FFT beat signals Bfd 1 -BfdN are summed up to generate a sum beat signal Bfd 0 having less noise than each of the FFT beat signals Bfd 1 -BfdN.
- each of the sum beat signals Bfu 0 , Bfd 0 has peak strength.
- the sum beat signal Bfu 0 has the peak strength at a frequency FP, where P is a positive integer less than or equal to M.
- the peak frequency FP is detected and an extraction matrix C(FP), which is an addition correlation matrix U(FP) corresponding to the peak frequency FP, is extracted from the addition correlation matrices U(F 1 )-U(FM) of the addition correlation matrix group UG.
- a process 500 including the distance calculating process and the direction calculating process is described.
- the beat signals Bu 1 -BuN are only discussed below, the beat signals Bd 1 -BdN are processed in the same way as the beat signals Bu 1 -BuN.
- the microcomputer 18 performs the process 500 as an interrupt process at the predetermined time interval Ts.
- the process 500 starts with step S 501 , where the microcomputer 18 obtains the beat signals Bu 1 -BdN.
- step S 502 the FFT is applied to each of the beat signals Bu 1 -BdN to generate the FFT beat signals Bfu 1 -BfuN.
- step S 503 where the FFT beat signals Bfu 1 -BfuN are summed into the sum beat signal Bfu 0 .
- step S 504 where the peak frequency FP of the sum beat signal Bfu 0 is detected.
- step S 505 where the correlation matrix group RG having the correlation matrices R(F 1 )-R(FM) is calculated from the FFT beat signals Bu 1 -BdN generated in step S 502 .
- step S 506 each of the correlation matrices R(F 1 )-R(FM) of the correlation matrix group RG is multiplied by a weighting factor (1 ⁇ K) and each of correlation matrices Ro(F 1 )-Ro(FM) of the previous correlation matrix group RoG is multiplied by a weighting factor K, where K is a fixed value between 0.0 and 1.0.
- the previous correlation matrix group RoG is the correlation matrix group RG that is generated in a previous loop (i.e., Ts earlier) of the process 500 .
- the correlation matrix group RG multiplied by the weighting factor (1 ⁇ K) and the previous correlation matrix group RoG multiplied by the weighting factor K are added together to produce the addition correlation matrix group UG having addition correlation matrices U(F 1 )-U(FM). Therefore, the addition correlation matrix U(FI), i.e., each of the addition correlation matrices U(F 1 )-(FM) of the addition correlation matrix group UG is given by:
- step S 507 the extraction matrix C(FP) is extracted from the addition correlation matrix group UG.
- the extraction matrix C(FP) is the addition correlation matrices U(FP) corresponding to the peak frequency FP detected in step S 504 ,
- step S 508 an eigenvalue expansion of the extraction matrix C(FP) is performed.
- step S 509 a MUSIC spectrum is calculated based on an eigenvector of the extraction matrix C(FP).
- step S 510 the direction of the object in the increase area is calculated based on the MUSIC spectrum. Because the beat signals Bd 1 -BdN are processed in the same way as the beat signals Bu 1 -BuN, the direction of the object in the decrease area is also calculated.
- step S 511 the correlation matrix group RG generated in step S 505 is stored in the memory 21 as the previous correlation matrix group RoG that is used in a next loop of the process 500 .
- step S 512 pair matching of the object is performed based on the strength of the sum beat signal Bfu 0 , the direction of the object in the increase area, the strength of the sum beat signal Bfd 0 , and the direction of the object in the decrease area.
- the distance and the relative velocity between the radar apparatus 100 and the object are detected.
- step S 512 After step S 512 is finished, the process 500 returns to step S 501 .
- the correlation matrix group RG generated in a present loop of the process 500 is added to the previous correlation matrix group RoG that is generated in the previous loop of the process 500 and stored in the memory 21 .
- the addition correlation matrix group UG can have less noise than the correlation matrix group RG.
- the peak frequency FP which indicates the presence of the object, is detected and the extraction matrix C(FP) corresponding to the peak frequency FP is extracted from the additional correlation matrix group UG.
- the MUSIC spectrum is generated by using the extraction matrix C(FP). Specifically, the eigenvalue expansion is performed on only the extraction matrix C(FP) to generate the MUSIC spectrum. Therefore, the amount of calculation executed by the microcomputer 18 is very small, as compared to when the eigenvalue expansion is performed on each of the addition correlation matrices U(F 1 )-U(FM).
- the previous information is stored in the memory 21 in the form of the correlation matrix for the following reason.
- FIG. 6 is a table showing the amount of calculation for converting the FFT beat signals Buf 0 -BufN into each form per frequency (i.e., one of the frequencies F 1 -FM).
- R ⁇ R represents an addition/subtraction of real numbers and R ⁇ R represents a multiplication of real numbers.
- 500 addition/subtractions and 500 multiplications are required to calculate one of the correlation matrices R(F 1 )-R(FM) from the FFT beat signals Bfu 1 -BfuN.
- 3500 addition/subtractions and 3500 multiplications are required to calculate the eigenvectors of one of the correlation matrices R(F 1 )-R(FM) from the FFT beat signals Bfu 1 -BfuN.
- 3000 addition/subtractions and 3000 multiplications are required to calculate the eigenvectors from one of the correlation matrices R(F 1 )-R(FM).
- 500 ⁇ M addition/subtractions and 500 ⁇ M multiplications are performed to calculate the correlation matrices R(F 1 )-R(FM) from the FFT beat signals Bfu 1 -BfuN.
- the number M is 10
- 50000 addition/subtractions and 50000 multiplications are performed to calculate the 10 correlation matrices R(F 1 )-R(F 10 ) from the FFT beat signals Bfu 1 -BfuN.
- the memory 21 stores the eigenvectors of each of the correlation matrices R(F 1 )-R(FM), 3000 ⁇ M addition/subtractions and 3000 ⁇ M multiplications are further performed to calculate the eigenvectors.
- the eigenvector of the correlation matrix R(FP) corresponding to the peak frequency FP is only used to generate the extraction matrix C(FP). Therefore, when the peak frequency in the previous loop is equal to that in the present loop, 3000 ⁇ (M ⁇ 1) addition/subtractions and 3000 ⁇ (M ⁇ 1) multiplications are wasted. Likewise, when the peak frequency in the previous loop is not equal to that in the present loop, 3000 ⁇ (M ⁇ 2) addition/subtractions and 3000 ⁇ (M ⁇ 2) multiplications are wasted.
- the previous information should be stored in the memory 21 in the form of the correlation matrices R(F 1 )-R(FM).
- the radar apparatus 100 achieves the high resolution without the increase in the amount of calculation.
- the radar apparatus 100 can detect two or more objects.
- the sum beat signal Bfu 0 has two peak strengths, i.e., two peak frequencies.
- two extraction matrices C(FP) one of which corresponds to one peak frequency and the other of which corresponds to the other peak frequency, are extracted from the additional correlation matrix group UG.
- the MUSIC spectrums are calculated based on each of the two extraction matrices C(FP) so that each direction of the two objects can be detected.
- the sum beat signal Bfu 0 has only one peak strength.
- one extraction matrix C(FP) corresponding to the peak frequency is extracted from the additional correlation matrix group UG.
- the MUSIC spectrum is calculated based on the extraction matrix C(FP). Because the MUSIC spectrum contains signals indicating each direction of the two objects, i.e., the MUSIC spectrum has two peaks, each direction of the two objects can be detected.
- the microcomputer 18 performs a process 700 shown in FIG. 7 instead of the process 500 shown in FIG. 5 .
- the process 700 includes steps S 705 -S 711 instead of steps S 505 -S 511 of the process 500 .
- the beat signals Bu 1 -BuN are only discussed below, the beat signals Bd 1 -BdN are processed in the same way as the beat signals Bu 1 -BuN.
- step S 705 the correlation matrix group RG having the correlation matrices R(F 1 )-R(FM) is calculated from the FFT beat signals Bfu 1 -BfuN generated in step S 502 .
- each of the correlation matrices R(F 1 )-R(FM) of the correlation matrix group RG is multiplied by the weighting factor (1 ⁇ K) and each of previous addition correlation matrices Uo(F 1 )-Uo(FM) of a previous correlation matrix group UoG is multiplied by a weighting factor K.
- the previous correlation matrix group RoG is the addition correlation matrix group UG generated in a previous loop (i.e., Ts earlier) of the process 700 .
- the correlation matrix group RG multiplied by the weighting factor (1 ⁇ K) and the previous addition correlation matrix group UoG multiplied by the weighting factor K are added together to produce the addition correlation matrix group UG having the addition correlation matrices U(F 1 )-U(FM). Therefore, the addition correlation matrix U(FI), i.e., each of the addition correlation matrices U(F 1 )-(FM) of the addition correlation matrix group UG is given by:
- step S 707 the extraction matrix C(FP) is extracted from the addition correlation matrix group UG.
- the extraction matrix C(FP) is the addition correlation matrices U(FP) corresponding to the peak frequency FP detected in step S 504 .
- step S 708 the eigenvalue expansion of the extraction matrix C(FP) is performed.
- step S 709 the MUSIC spectrum is calculated based on the eigenvector of the extraction matrix C(FP).
- step S 710 the direction of the object in the increase area is calculated based on the MUSIC spectrum. Because the beat signals Bd 1 -BdN are processed in the same way as the beat signals Bu 1 -BuN, the direction of the object in the decrease area is also calculated.
- step S 711 the addition correlation matrix group UG generated in step S 706 is stored in the memory 21 as the previous addition correlation matrix group UoG that is used in a next loop of the process 700 .
- step S 512 the process 700 proceeds to step S 512 .
- the previous additional correlation matrix UoG generated in the previous loop is used to generate the addition correlation matrix group UG.
- the addition correlation matrix UG can be generated based on two or more previous correlation matrix groups so that the addition correlation matrix group UG of the second embodiment can have less noise than that of the first embodiment.
- the microcomputer 18 performs a process 800 shown in FIG. 8 instead of the process 500 shown in FIG. 5 .
- the process 800 includes steps S 805 -S 811 instead of steps S 505 -S 511 of the process 500 .
- the beat signals Bu 1 -BuN are only discussed below, the beat signals Bd 1 -BdN are processed in the same way as the beat signals Bu 1 -BuN.
- step S 805 the correlation matrix group RG having the correlation matrices R(F 1 )-R(FM) is calculated from the FFT beat signals Bfu 1 -BfuN generated in step S 502 .
- step S 806 a correlation matrix R(FP) corresponding to the peak frequency FP detected in step S 504 is extracted from the correlation matrix group RG. Further, a previous correlation matrix Ro(FP) corresponding to the peak frequency FP is extracted from a previous correlation matrix group RoG having correlation matrices Ro(F 1 )-Ro(FM).
- the previous correlation matrix group RoG is the correlation matrix group RG that is generated in a previous loop (i.e., Ts earlier) of the process 800 and stored in the memory 21 .
- step S 807 the correlation matrix R(FP) is multiplied by the weighting factor (1 ⁇ K) and the previous correlation matrix Ro(FP) is multiplied by the weighting factor K. Then, the correlation matrix R(FP) multiplied by the weighting factor (1 ⁇ K) and the previous correlation matrix Ro(FP) multiplied by the weighting factor K are added together to produce the addition correlation matrix U(FP). Therefore, the addition correlation matrix U(FP) is given by:
- step S 808 the eigenvalue expansion of the addition correlation matrix U(FP) is performed.
- step S 809 the MUSIC spectrum is calculated based on the eigenvector of the addition correlation matrix U(FP).
- step S 810 the direction of the object in the increase area is calculated based on the MUSIC spectrum. Because the beat signals Bd 1 -BdN are processed in the same way as the beat signals Bu 1 -BuN, the direction of the object in the decrease area is also calculated.
- step S 811 the correlation matrix group RG generated in step S 805 is stored in the memory 21 as the previous correlation matrix group RoG that is used in a next loop of the process 800 .
- step S 512 the process 800 proceeds to step S 512 .
- the addition correlation matrix group UG having the addition correlation matrices U(F 1 )-U(FM) is generated such that the correlation matrix group RG is added to the previous correlation matrix group RoG.
- the sum beat signal Bfu 0 has only one peak frequency. Therefore, although each of the addition correlation matrices U(F 1 )-U(FM) is calculated, the addition correlation matrix U(FP) corresponding to the peak frequency FP is only used. In other words, the calculation of the addition correlation matrices U(F 1 )-U(FM) except for the addition correlation matrix U(FP) may result in waste.
- the addition correlation matrix U(FP) is generated such that the correlation matrix R(FP) is added to the previous correlation matrix Ro(FP). Thus, the wasted calculation can be avoided.
- the memory 21 needs to store the previous correlation matrix group RoG and the addition correlation matrix group UG at the same time.
- the memory 21 needs to store the previous correlation matrix group RoG and the addition correlation matrix U(FP) at the same time. Therefore, the memory 21 can have a small amount of storage capacity in the process 800 , as compared to in the process 500 .
- the memory 21 stores each of the M correlation matrices R(F 1 )-R(FM) for the next loop.
- the memory 21 stores M/2 correlation matrices R(F 1 ), R(F 3 ), R(F 5 ) • • • .
- the M correlation matrices R(F 1 )-R(FM) are thinned out to the M/2 correlation matrices R(F 1 ), R(F 3 ), R(F 5 ) • • • .
- the M correlation matrices R(F 1 )-R(FM) are alternately stored in the memory 21 such that the memory 21 stores the M/2 correlation matrices R(F 1 ), R(F 3 ), R(F 5 ) • • • • .
- the memory 21 stores the correlation matrices R(F 1 ), R(F 3 ) as the previous correlation matrices Ro(F 1 ), Ro(F 3 ).
- the correlation matrix R(F 2 ) is not stored in the memory 21 in the Sth loop.
- the peak frequency FP is F 2 in the (S+1)th loop
- the previous correlation matrix Ro(F 2 ) is generated such that a weighted average of the previous correlation matrix Ro(F 1 ) is added to a weighted average of the previous correlation matrix Ro(F 3 ).
- the microcomputer 18 performs a process 900 shown in FIG. 9 instead of the process 500 shown in FIG. 5 .
- the process 900 includes steps S 905 -S 912 instead of steps S 505 -S 511 of the process 500 .
- the beat signals Bu 1 -BuN are only discussed below, the beat signals Bd 1 -BdN are processed in the same way as the beat signals Bu 1 -BuN.
- step S 905 the correlation matrix group RG having the correlation matrices R(F 1 )-R(FM) is calculated from the FFT beat signals Bfu 1 -BfuN generated in step S 502 .
- step S 906 it is determined whether the previous correlation matrix Ro(FP) corresponding to the peak frequency FP is stored in the memory 21 .
- step S 908 If the previous correlation matrix Ro(FP) is stored in the memory 21 , the process 900 proceeds to step S 908 directly.
- step S 908 the previous correlation matrix Ro(FP) is generated such that a weighted average of a previous correlation matrix Ro(FP ⁇ 1) is added to a weighted average of a previous correlation matrix Ro(FP+1).
- the correlation matrix R(FP) corresponding to the peak frequency FP is extracted from the correlation matrix group RG.
- the correlation matrix R(FP) is multiplied by the weighting factor (1 ⁇ K) and the previous correlation matrix Ro(FP) is multiplied by the weighting factor K.
- the correlation matrix R(FP) multiplied by the weighting factor (1 ⁇ K) and the previous correlation matrix Ro(FP) multiplied by the weighting factor K are added together to produce the addition correlation matrix U(FP). Therefore, the addition correlation matrix U(FP) is given by:
- step S 909 the eigenvalue expansion of the addition correlation matrix U(FP) is performed.
- step S 910 the MUSIC spectrum is calculated based on the eigenvector of the addition correlation matrix U(FP).
- step S 911 the direction of the object in the increase area is calculated based on the MUSIC spectrum. Because the beat signals Bd 1 -BdN are processed in the same way as the beat signals Bu 1 -BuN, the direction of the object in the decrease area is also calculated.
- step S 912 the correlation matrix group RG is thinned out and stored in the memory 21 as the previous correlation matrix group RoG that is used in a next loop of the process 900 .
- step S 512 the process 900 proceeds to step S 512 .
- the correlation matrix group RG is thinned out and stored in the memory 21 as the previous correlation matrix group RoG.
- the memory 21 can have a small amount of storage capacity in the process 900 , as compared to in the process 800 according to the third embodiment. Even when the previous correlation matrix Ro(FP) is not stored in the memory 21 , the previous correlation matrix Ro(FP) is estimated from previous correlation matrices Ro(FP ⁇ 1), Ro(FP+1).
- the radar apparatus 100 may use the received radar signal received by some of the elements E 1 -EN of the receiving antenna 12 , not each of the elements E 1 -EN. In such an approach, the amount of calculation can be reduced.
- the correlation matrix group RG or the addition correlation matrix group UG may be stored in the memory 21 after being compressed by a data compression algorithm.
- the memory 21 may store the previous information in a form of data containing elements of the correlation matrix group RG or the addition correlation matrix group UG.
- the weighting factor K may be a variable. For example, when the received radar signal has considerable instantaneous noise, the weighing factor K may be increased. In such an approach, an influence of the noise can be reduced.
- the algorithm used in the present invention can be applied to various types of the DOA estimation algorithms such as unitary-MUSIC algorithm, ESPRIT algorithm, unitary-ESPRIT algorithm, Capon algorithm, and Beam Former algorithm.
- unitary-MUSIC or the unitary-ESPRIT is used, only the real part of the matrix is stored in the memory 21 . Therefore, the amount of calculation can be significantly reduced and the memory 21 can have a very small amount of storage capacity.
- the algorithm used in the present invention also can be applied to a spatial smoothing algorithm.
- the previous information (i.e., the previous correlation matrix group RoG, or the previous addition correlation matrix group UG) may be generated in two or more previous loop.
- the previous information may be generated in two previous loop (i.e., generated 2Ts earlier).
- the transmitting antenna 11 instead of the receiving antenna 12 may have the elements arranged in the array to generate the beat signals.
Abstract
Description
- This application is based on and incorporates herein by reference Japanese Patent Application No. 2005-224636 filed on Aug. 2, 2005.
- The present invention relates to a radar apparatus.
- There are various types of radar apparatus for detecting a distance and a direction of an object with respect to the radar apparatus. A frequency modulated continuous wave (FMCW) radar apparatus is a radar apparatus that continuously transmits a frequency-modulated radar signal to the object and detects the distance or a relative velocity of the object based on the transmitted radar signal reflected from the object.
- In one method for detecting the direction of the object, a transmitting means for transmitting a radar signal mechanically turns and scans the transmitted radar signal reflected from the object. In another method that uses a digital beam forming (DBF) algorithm, the transmitting means is fixed and an antenna having multiple elements arranged in an array receives the transmitted radar signal. The received radar signal is digitally processed and the direction of the object is detected based on the digital signal. Specifically, in the DBF algorithm, an angular spectrum is generated based on the received radar signal on each of the elements and a peak of the angular spectrum is detected. The direction of the object is estimated based on the peak of the angular spectrum.
- In a beamformer algorithm as the DBF algorithm, the angular spectrum is generated such that amplitudes of the received signal at a given time are connected as shown in
FIGS. 10A and 10B . A multiple signal classification (MUSIC) algorithm is known as a high-resolution direction of arrival (DOA) estimation algorithm. - In the DOA estimation algorithm, correlation matrices are calculated, eigenvalue expansions are performed on each of the correlation matrices, the angular spectrum is calculated from eigenvectors of the correlation matrices, and the direction of the object is calculated based on the angular spectrum.
- A FMCW radar apparatus disclosed in U.S. Pat. No. 6,121,917 corresponding to JP-A-H11-133142 detects the direction of the object by using the beamformer algorithm. In the FMCW radar apparatus, a fast Fourier transform (FFT) is applied to the received wave signal to obtain a peak frequency of a distance power spectrum. Then, the beamformer algorithm is applied to only the peak frequency component of the received signal so that the amount of calculation required to detect the direction of the object is reduced.
- However, when the beamformer algorithm is used in the radar apparatus, resolution of the radar apparatus depends on the number of elements arranged in the array. Therefore, the radar apparatus using the beamformer algorithm needs to be increased in size to obtain high resolution.
- The high-resolution DOA estimation algorithm such as the MUSIC algorithm achieves the high resolution without an increase in the number of the elements. In the DOA estimation algorithm, the resolution may be increased by reducing noise with summation of the received signal with respect to time. The summation is performed such that a present angular spectrum calculated in a present process and a previous angular spectrum calculated in a previous process are summed up.
- However, when the object moves, a frequency corresponding to the distance changes between in the previous process and in the present process. Therefore, the DOA estimation algorithm needs to be applied to all the frequency components of the received signal on each process to calculate the angular spectrum, and the calculated angular spectrum needs to be stored in a memory. The calculation of the angular spectrum requires an eigenvalue expansion that requires a lot of calculation. Therefore, when the angular spectrum is calculated on each frequency component, the amount of calculation is significantly increased. In the DOA estimation algorithm, the high resolution results in a significant increase in the amount of calculation.
- In view of the above-described problem, it is an object of the present invention to provide a radar apparatus that achieves a high resolution without a significant increase in the amount of calculation.
- A radar apparatus uses a DOA estimation algorithm performed at a predetermined time interval to detect an object. The radar apparatus includes a transmitting means that transmits a wave signal to the object, a receiving means that receives the transmitted wave signal reflected from the object, a beat signal generation means that generates a plurality of beat signals based on the received wave signal, a correlation matrix generation means that calculates a plurality of correlation matrices based on the beat signals with respect to each frequency component of the beat signals, an addition means that calculates a plurality of addition correlation matrices, a storing means that stores at least a portion of a plurality of previous correlation matrices calculated by the correlation matrix generation means at a previous time, a detection means that detects at least one frequency component of the beat signals, the frequency component satisfying a predetermined condition, an extraction means that extracts an extraction matrix from the addition correlation matrices, the extraction matrix corresponding to a frequency component closest to the frequency component detected by the detection means, and a direction calculation means that calculates a direction of the object with respect to the radar apparatus based on the extraction matrix. The receiving means includes a plurality of elements arranged in an array and each of the beat signals generated by the beat signal generation means corresponds to each of the elements. The addition means calculates the addition correlation matrices by adding the correlation matrices calculated by the correlation matrix generation means to the previous correlation matrices stored by the storing means.
- In a conventional DOA estimation algorithm, correlation matrices are generated from beat signals. Then, angular spectrums are calculated from each eigenvector of each of the correlation matrices. The storing means stores the angular spectrums as previous information and the angular spectrums are used in a next process to reduce noise.
- In the DOA estimation algorithm used in the radar apparatus, a position of the object (i.e., a frequency component indicating a presence of the object) is estimated from the beat signals and then an angular spectrum corresponding to the position is calculated. Correlation matrices calculated in a present process and previous correlation matrices calculated in a previous process are summed up to reduce noise. In the radar apparatus, therefore, the storing means stores the correlation matrices as the previous information, not the angular spectrums.
- For example, in a MUSIC algorithm, the amount of calculation required to calculate the correlation matrices is about one-thirteenth the amount of calculation required to calculate MUSIC spectrums. Because the radar apparatus uses the correlation matrices instead of the angular spectrums to reduce the noise, a significant increase in the amount of calculation can be prevented. Thus, the radar apparatus achieves a high resolution without the significant increase in the amount of calculation.
- The above and other objectives, features and advantages of the present invention will become more apparent from the following detailed description made with reference to the accompanying drawings. In the drawings:
-
FIG. 1 is a block diagram showing a radar apparatus according to a first embodiment of the present invention; -
FIG. 2 is an arrangement of elements of a receiving antenna of the radar apparatus ofFIG. 1 ; -
FIGS. 3A and 3B are diagrams showing a principal of beat signals generated in the radar apparatus ofFIG. 1 ; -
FIG. 4A is a diagram showing a reflected wave signal received by the elements of the receiving antenna of the radar apparatus ofFIG. 1 ,FIG. 4B is FFT beat signals generated by applying a FFT to beat signals ofFIG. 3B , andFIG. 4C is a diagram showing a sum beat signal into which the FFT beat signals ofFIG. 4B are summed; -
FIG. 5 is a flow chart illustrating a process performed by a microcomputer of the radar apparatus ofFIG. 1 ; -
FIG. 6 is a table showing the amount of calculation performed by the microcomputer of the radar apparatus ofFIG. 1 ; -
FIG. 7 is a flow chart illustrating a process performed by a microcomputer of a radar apparatus according to a second embodiment of the present invention; -
FIG. 8 is a flow chart illustrating a process performed by a microcomputer of a radar apparatus according to a third embodiment of the present invention; -
FIG. 9 is a flow chart illustrating a process performed by a microcomputer of a radar apparatus according to a fourth embodiment of the present invention; and -
FIG. 10A is a diagram showing an angular spectrum generated when reflected wave signals arrive at elements from a front direction in a beam former method, andFIG. 10B is a diagram showing the angular spectrum generated when the reflected wave signals arrive at the elements from an oblique direction in the beam former method. - A
radar apparatus 100 according to a first embodiment of the present invention will now be described withFIGS. 1 to 5 . Theradar apparatus 100 includes a transmittingantenna 11, a receivingantenna 12 having N elements E1-EN arranged in an array, where N is a positive integer, a high-frequency switch 13, amixer 14, anoscillator 15, a digital-to-analog (D/A)converter 16, an analog-to-digital (A/D)converter 17, amicrocomputer 18, aswitch controller 19, and atimer 20. - The
microcomputer 18 outputs a digital transmit signal to the D/A converter 16. The D/A converter 16 converts the digital transmit signal into an analog transmit signal and outputs the analog transmit signal to theoscillator 15. Theoscillator 15 outputs the analog transmit signal to the transmittingantenna 11 at a predetermined frequency. The transmittingantenna 11 converts the analog transmit signal into a radar signal and transmits the radar signal to an object. - The receiving
antenna 12 receives the transmitted radar signal reflected from the object. The received radar signal is converted into an analog receive signal on each of the elements E1-EN. The high-frequency switch 13 in turn sends the analog receive signal to themixer 14. In themixer 14, the analog receive signal is mixed with the analog transmit signal to generate a beat signal on each of the elements E1-EN. The beat signal is input to the A/D converter 17 and converted into a digital receive signal. The digital signal is input to themicrocomputer 18. - The
microcomputer 18 controls the high-frequency switch 13 through theswitch controller 19 and controls a sampling rate of the A/D converter 17 though thetimer 20. Themicrocomputer 18 has amemory 21. - The
microcomputer 18 performs a distance calculation process for calculating the distance between theradar apparatus 100 and the object and a direction calculation process for calculating the direction between theradar apparatus 100 and the object. - Referring to
FIGS. 2 and 3 , the distance calculation process is described. - As shown in
FIG. 2 , the elements E1-EN of the receivingantenna 12 are spaced from each other by a predetermined spacing S. - There arises a time delay and a frequency shift between the transmitted radar signal, which is frequently modulated by the
oscillator 15, and the received radar signal. The time delay corresponds to the distance between theradar apparatus 100 and the object, and the frequency shift corresponds to a relative velocity between theradar apparatus 100 and the object. A phase shift of the received radar signal with respect to the transmitted radar signal increases with the distance between theradar apparatus 100 and the object. The distance and relative velocity between theradar apparatus 100 and the object can be detected based on the phase shift. Therefore, a beat signal is calculated as a frequency difference between the transmitted radar signal and the received radar signal. - As shown in
FIGS. 3A and 3B , the beat signal has a first beat frequency Bu in an increase area where the frequency of the transmitted radar signal increases, and has a second beat frequency Bd in a decrease area where the frequency of the transmitted radar signal decreases. Thus, the beat signal includes a first beat signal having the first beat frequency Bu and a second beat signal having the second beat frequency Bd. - When the receiving
antenna 12 has the elements E1-EN, the beat signal is generated on each of the elements E1-EN. In the whole receivingantenna 12, therefore, 2N beat signals are generated. Specifically, N first beat signals having the first beat frequency Bu and N second beat signals having the second beat frequency Bd are generated. - The distance and the relative velocity between the
radar apparatus 100 and the object are given by the following equations: -
D={C·T/(4·ΔF)}·(BuH+BdH) (1) -
V={C/(4·F0)}·(BuH−BdH) (2) - In the above equations (1) (2), D represents the distance, V represents the relative velocity, C represents the speed of light, ΔF represents a frequency range of the transmitted radar signal, and F0 represents the center of the ΔF. BuH is the first beat signal generated based on the received radar signal that is received by the element EH, where H is a positive integer less than or equal to N (i.e., 1≦H≦N). BdH is the second beat signal generated based on the received radar signal that is received by the element EH.
- Referring to
FIGS. 4A-4C , the direction calculation process is described. Themicrocomputer 18 performs the direction calculation process at a predetermined time interval Ts. - First, a fast Fourier transform (FFT) is applied to each of the N beat signals Bu1-BuN to generate N FFT beat signals Bfu1-BfuN. Likewise, the fast Fourier transform (FFT) is applied to the N beat signals Bd1-BdN to generate N FFT beat signals Bfd1-BfdN. Although only the FFT beat signals Bfu1-BfuN are illustrated in
FIG. 4B , the FFT beat signals Bfd1-BfdN are processed in the same way as the FFT beat signals Bfu1-BfuN. - Next, a correlation matrix group RG of correlation matrices are generated by using each of the FFT beat signals. For example, in the case of
FIG. 4B , the correlation matrix group RG includes M correlation matrices R(F1)-R(FM), where M is a positive integer greater than 1. A correlation matrix R(FI) corresponds to a frequency FI, where I is a positive integer less than or equal to M (i.e., 1≦I≦M). When the receivingantenna 12 has the elements E1-EN, each of the correlation matrices R(F1)-R(FM) is an N×N matrix. - As described later, an addition correlation matrix group UG having addition correlation matrices U(F1)-U(FM) is generated such that the correlation matrix group RG generated in a present process is added to a previous correlation matrix group RoG that is the correlation matrix group RG generated in a previous process (i.e., Ts earlier) and stored in the
memory 21. The addition correlation matrix group UG has less noise than the correlation matrix group RG. The use of the previous correlation matrix group RoG as previous information reduces the noise. - A MUSIC algorithm used in the
radar apparatus 100 is described below. The MUSIC algorithm allows theradar apparatus 100 to achieve the high resolution without an increase in the amount of calculation. - The FFT beat signals Bfu1-BfuN shown in
FIG. 4B are summed up to generate a sum beat signal Bfu0 shown inFIG. 4C . As can be seen fromFIG. 4C , the sum beat signal Bfu0 has less noise than each of the FFT beat signals Bfu1-BfuN. Likewise, the FFT beat signals Bfd1-BfdN are summed up to generate a sum beat signal Bfd0 having less noise than each of the FFT beat signals Bfd1-BfdN. - When the received radar signal contains a reflected wave from the object, each of the sum beat signals Bfu0, Bfd0 has peak strength. For example, in
FIG. 4C , the sum beat signal Bfu0 has the peak strength at a frequency FP, where P is a positive integer less than or equal to M. The peak frequency FP is detected and an extraction matrix C(FP), which is an addition correlation matrix U(FP) corresponding to the peak frequency FP, is extracted from the addition correlation matrices U(F1)-U(FM) of the addition correlation matrix group UG. - Referring to
FIG. 5 , aprocess 500 including the distance calculating process and the direction calculating process is described. Although the beat signals Bu1-BuN are only discussed below, the beat signals Bd1-BdN are processed in the same way as the beat signals Bu1-BuN. Themicrocomputer 18 performs theprocess 500 as an interrupt process at the predetermined time interval Ts. - The
process 500 starts with step S501, where themicrocomputer 18 obtains the beat signals Bu1-BdN. - Then, the
process 500 proceeds to step S502, where the FFT is applied to each of the beat signals Bu1-BdN to generate the FFT beat signals Bfu1-BfuN. - Then, the
process 500 proceeds to step S503, where the FFT beat signals Bfu1-BfuN are summed into the sum beat signal Bfu0. - Then, the
process 500 proceeds to step S504, where the peak frequency FP of the sum beat signal Bfu0 is detected. - Then, the
process 500 proceeds to step S505, where the correlation matrix group RG having the correlation matrices R(F1)-R(FM) is calculated from the FFT beat signals Bu1-BdN generated in step S502. - Then, the
process 500 proceeds to step S506, where each of the correlation matrices R(F1)-R(FM) of the correlation matrix group RG is multiplied by a weighting factor (1−K) and each of correlation matrices Ro(F1)-Ro(FM) of the previous correlation matrix group RoG is multiplied by a weighting factor K, where K is a fixed value between 0.0 and 1.0. As described above, the previous correlation matrix group RoG is the correlation matrix group RG that is generated in a previous loop (i.e., Ts earlier) of theprocess 500. Then, the correlation matrix group RG multiplied by the weighting factor (1−K) and the previous correlation matrix group RoG multiplied by the weighting factor K are added together to produce the addition correlation matrix group UG having addition correlation matrices U(F1)-U(FM). Therefore, the addition correlation matrix U(FI), i.e., each of the addition correlation matrices U(F1)-(FM) of the addition correlation matrix group UG is given by: -
U(FI)=R(FI)·(1−K)+Ro(FI)·K (3) - Then, the
process 500 proceeds to step S507, where the extraction matrix C(FP) is extracted from the addition correlation matrix group UG. The extraction matrix C(FP) is the addition correlation matrices U(FP) corresponding to the peak frequency FP detected in step S504, - Then, the
process 500 proceeds to step S508, where an eigenvalue expansion of the extraction matrix C(FP) is performed. - Then, the
process 500 proceeds to step S509, where a MUSIC spectrum is calculated based on an eigenvector of the extraction matrix C(FP). - Then, the
process 500 proceeds to step S510, where the direction of the object in the increase area is calculated based on the MUSIC spectrum. Because the beat signals Bd1-BdN are processed in the same way as the beat signals Bu1-BuN, the direction of the object in the decrease area is also calculated. - Then, the
process 500 proceeds to step S511, where the correlation matrix group RG generated in step S505 is stored in thememory 21 as the previous correlation matrix group RoG that is used in a next loop of theprocess 500. - Then, the
process 500 proceeds to step S512, where pair matching of the object is performed based on the strength of the sum beat signal Bfu0, the direction of the object in the increase area, the strength of the sum beat signal Bfd0, and the direction of the object in the decrease area. Thus, the distance and the relative velocity between theradar apparatus 100 and the object are detected. - After step S512 is finished, the
process 500 returns to step S501. - The correlation matrix group RG generated in a present loop of the
process 500 is added to the previous correlation matrix group RoG that is generated in the previous loop of theprocess 500 and stored in thememory 21. In such an approach, the addition correlation matrix group UG can have less noise than the correlation matrix group RG. - The peak frequency FP, which indicates the presence of the object, is detected and the extraction matrix C(FP) corresponding to the peak frequency FP is extracted from the additional correlation matrix group UG. The MUSIC spectrum is generated by using the extraction matrix C(FP). Specifically, the eigenvalue expansion is performed on only the extraction matrix C(FP) to generate the MUSIC spectrum. Therefore, the amount of calculation executed by the
microcomputer 18 is very small, as compared to when the eigenvalue expansion is performed on each of the addition correlation matrices U(F1)-U(FM). - The previous information is stored in the
memory 21 in the form of the correlation matrix for the following reason. -
FIG. 6 is a table showing the amount of calculation for converting the FFT beat signals Buf0-BufN into each form per frequency (i.e., one of the frequencies F1-FM). In the table, R±R represents an addition/subtraction of real numbers and R×R represents a multiplication of real numbers. - As show in the table, 500 addition/subtractions and 500 multiplications are required to calculate one of the correlation matrices R(F1)-R(FM) from the FFT beat signals Bfu1-BfuN. Likewise, 3500 addition/subtractions and 3500 multiplications are required to calculate the eigenvectors of one of the correlation matrices R(F1)-R(FM) from the FFT beat signals Bfu1-BfuN. In other words, 3000 addition/subtractions and 3000 multiplications are required to calculate the eigenvectors from one of the correlation matrices R(F1)-R(FM).
- In the case of
FIG. 4 , 500×M addition/subtractions and 500×M multiplications are performed to calculate the correlation matrices R(F1)-R(FM) from the FFT beat signals Bfu1-BfuN. For example, when the number M is 10, 50000 addition/subtractions and 50000 multiplications are performed to calculate the 10 correlation matrices R(F1)-R(F10) from the FFT beat signals Bfu1-BfuN. - If the
memory 21 stores the eigenvectors of each of the correlation matrices R(F1)-R(FM), 3000×M addition/subtractions and 3000×M multiplications are further performed to calculate the eigenvectors. However, the eigenvector of the correlation matrix R(FP) corresponding to the peak frequency FP is only used to generate the extraction matrix C(FP). Therefore, when the peak frequency in the previous loop is equal to that in the present loop, 3000×(M−1) addition/subtractions and 3000×(M−1) multiplications are wasted. Likewise, when the peak frequency in the previous loop is not equal to that in the present loop, 3000×(M−2) addition/subtractions and 3000×(M−2) multiplications are wasted. - In view of the amount of calculation, therefore, it is appropriate that the previous information should be stored in the
memory 21 in the form of the correlation matrices R(F1)-R(FM). - Thus, the
radar apparatus 100 achieves the high resolution without the increase in the amount of calculation. - Although the case where one object is detected is discussed in the first embodiment, the
radar apparatus 100 can detect two or more objects. - For example, when the number of the objects is two and the distance between one object and the
radar apparatus 100 is not equal to that between the other object and theradar apparatus 100, the sum beat signal Bfu0 has two peak strengths, i.e., two peak frequencies. In this case, two extraction matrices C(FP), one of which corresponds to one peak frequency and the other of which corresponds to the other peak frequency, are extracted from the additional correlation matrix group UG. The MUSIC spectrums are calculated based on each of the two extraction matrices C(FP) so that each direction of the two objects can be detected. - In contrast, when the distance between one object and the
radar apparatus 100 is equal to that between the other object and theradar apparatus 100, the sum beat signal Bfu0 has only one peak strength. In this case, one extraction matrix C(FP) corresponding to the peak frequency is extracted from the additional correlation matrix group UG. The MUSIC spectrum is calculated based on the extraction matrix C(FP). Because the MUSIC spectrum contains signals indicating each direction of the two objects, i.e., the MUSIC spectrum has two peaks, each direction of the two objects can be detected. - Referring to
FIGS. 5 and 7 , a second embodiment of the present invention is described. In the second embodiment, themicrocomputer 18 performs aprocess 700 shown inFIG. 7 instead of theprocess 500 shown inFIG. 5 . As shown inFIG. 7 , theprocess 700 includes steps S705-S711 instead of steps S505-S511 of theprocess 500. Although the beat signals Bu1-BuN are only discussed below, the beat signals Bd1-BdN are processed in the same way as the beat signals Bu1-BuN. - After steps S501-S504 are finished, the
process 700 proceeds to step S705, where the correlation matrix group RG having the correlation matrices R(F1)-R(FM) is calculated from the FFT beat signals Bfu1-BfuN generated in step S502. - Then, the
process 700 proceeds to step S706, where each of the correlation matrices R(F1)-R(FM) of the correlation matrix group RG is multiplied by the weighting factor (1−K) and each of previous addition correlation matrices Uo(F1)-Uo(FM) of a previous correlation matrix group UoG is multiplied by a weighting factor K. The previous correlation matrix group RoG is the addition correlation matrix group UG generated in a previous loop (i.e., Ts earlier) of theprocess 700. Then, the correlation matrix group RG multiplied by the weighting factor (1−K) and the previous addition correlation matrix group UoG multiplied by the weighting factor K are added together to produce the addition correlation matrix group UG having the addition correlation matrices U(F1)-U(FM). Therefore, the addition correlation matrix U(FI), i.e., each of the addition correlation matrices U(F1)-(FM) of the addition correlation matrix group UG is given by: -
U(FI)=R(FI)·(1−K)+Uo(FI)·K (4) - Then, the
process 700 proceeds to step S707, where the extraction matrix C(FP) is extracted from the addition correlation matrix group UG. The extraction matrix C(FP) is the addition correlation matrices U(FP) corresponding to the peak frequency FP detected in step S504. - Then, the
process 700 proceeds to step S708, where the eigenvalue expansion of the extraction matrix C(FP) is performed. - Then, the
process 700 proceeds to step S709, where the MUSIC spectrum is calculated based on the eigenvector of the extraction matrix C(FP). - Then, the
process 700 proceeds to step S710, where the direction of the object in the increase area is calculated based on the MUSIC spectrum. Because the beat signals Bd1-BdN are processed in the same way as the beat signals Bu1-BuN, the direction of the object in the decrease area is also calculated. - Then, the
process 700 proceeds to step S711, where the addition correlation matrix group UG generated in step S706 is stored in thememory 21 as the previous addition correlation matrix group UoG that is used in a next loop of theprocess 700. - Then, the
process 700 proceeds to step S512. - In the
process 700, thus, the previous additional correlation matrix UoG generated in the previous loop is used to generate the addition correlation matrix group UG. In such an approach, the addition correlation matrix UG can be generated based on two or more previous correlation matrix groups so that the addition correlation matrix group UG of the second embodiment can have less noise than that of the first embodiment. - Referring to
FIGS. 5 and 8 , a third embodiment of the present invention is described. In the third embodiment, themicrocomputer 18 performs aprocess 800 shown inFIG. 8 instead of theprocess 500 shown inFIG. 5 . As shown inFIG. 8 , theprocess 800 includes steps S805-S811 instead of steps S505-S511 of theprocess 500. Although the beat signals Bu1-BuN are only discussed below, the beat signals Bd1-BdN are processed in the same way as the beat signals Bu1-BuN. - After steps S501-S504 are finished, the
process 800 proceeds to step S805, where the correlation matrix group RG having the correlation matrices R(F1)-R(FM) is calculated from the FFT beat signals Bfu1-BfuN generated in step S502. - Then, the
process 800 proceeds to step S806, where a correlation matrix R(FP) corresponding to the peak frequency FP detected in step S504 is extracted from the correlation matrix group RG. Further, a previous correlation matrix Ro(FP) corresponding to the peak frequency FP is extracted from a previous correlation matrix group RoG having correlation matrices Ro(F1)-Ro(FM). The previous correlation matrix group RoG is the correlation matrix group RG that is generated in a previous loop (i.e., Ts earlier) of theprocess 800 and stored in thememory 21. - Then, the
process 800 proceeds to step S807, where the correlation matrix R(FP) is multiplied by the weighting factor (1−K) and the previous correlation matrix Ro(FP) is multiplied by the weighting factor K. Then, the correlation matrix R(FP) multiplied by the weighting factor (1−K) and the previous correlation matrix Ro(FP) multiplied by the weighting factor K are added together to produce the addition correlation matrix U(FP). Therefore, the addition correlation matrix U(FP) is given by: -
U(FP)=R(FP)·(1−K)+Ro(FP)·K (5) - Then, the
process 800 proceeds to step S808, where the eigenvalue expansion of the addition correlation matrix U(FP) is performed. - Then, the
process 800 proceeds to step S809, where the MUSIC spectrum is calculated based on the eigenvector of the addition correlation matrix U(FP). - Then, the
process 800 proceeds to step S810, where the direction of the object in the increase area is calculated based on the MUSIC spectrum. Because the beat signals Bd1-BdN are processed in the same way as the beat signals Bu1-BuN, the direction of the object in the decrease area is also calculated. - Then, the
process 800 proceeds to step S811, where the correlation matrix group RG generated in step S805 is stored in thememory 21 as the previous correlation matrix group RoG that is used in a next loop of theprocess 800. - Then, the
process 800 proceeds to step S512. - In the
process 500 according to the first embodiment, the addition correlation matrix group UG having the addition correlation matrices U(F1)-U(FM) is generated such that the correlation matrix group RG is added to the previous correlation matrix group RoG. When the number of the objects is one, the sum beat signal Bfu0 has only one peak frequency. Therefore, although each of the addition correlation matrices U(F1)-U(FM) is calculated, the addition correlation matrix U(FP) corresponding to the peak frequency FP is only used. In other words, the calculation of the addition correlation matrices U(F1)-U(FM) except for the addition correlation matrix U(FP) may result in waste. In contrast, in theprocess 800 according to the third embodiment, the addition correlation matrix U(FP) is generated such that the correlation matrix R(FP) is added to the previous correlation matrix Ro(FP). Thus, the wasted calculation can be avoided. - In the
process 500, thememory 21 needs to store the previous correlation matrix group RoG and the addition correlation matrix group UG at the same time. In contrast, in theprocess 800, thememory 21 needs to store the previous correlation matrix group RoG and the addition correlation matrix U(FP) at the same time. Therefore, thememory 21 can have a small amount of storage capacity in theprocess 800, as compared to in theprocess 500. - A fourth embodiment of the present invention is described. In the third embodiment, the
memory 21 stores each of the M correlation matrices R(F1)-R(FM) for the next loop. In contrast, in the fourth embodiment, thememory 21 stores M/2 correlation matrices R(F1), R(F3), R(F5) • • • . Thus, the M correlation matrices R(F1)-R(FM) are thinned out to the M/2 correlation matrices R(F1), R(F3), R(F5) • • • . In other words, the M correlation matrices R(F1)-R(FM) are alternately stored in thememory 21 such that thememory 21 stores the M/2 correlation matrices R(F1), R(F3), R(F5) • • • . - For example, when the correlation matrices R(F1)-R(F3) are generated in the Sth loop, where S is a positive integer, the
memory 21 stores the correlation matrices R(F1), R(F3) as the previous correlation matrices Ro(F1), Ro(F3). In other words, the correlation matrix R(F2) is not stored in thememory 21 in the Sth loop. In this case, if the peak frequency FP is F2 in the (S+1)th loop, the previous correlation matrix Ro(F2) is generated such that a weighted average of the previous correlation matrix Ro(F1) is added to a weighted average of the previous correlation matrix Ro(F3). - In the fourth embodiment, the
microcomputer 18 performs aprocess 900 shown inFIG. 9 instead of theprocess 500 shown inFIG. 5 . As shown inFIG. 9 , theprocess 900 includes steps S905-S912 instead of steps S505-S511 of theprocess 500. Although the beat signals Bu1-BuN are only discussed below, the beat signals Bd1-BdN are processed in the same way as the beat signals Bu1-BuN. - After steps S501-S504 are finished, the
process 900 proceeds to step S905, where the correlation matrix group RG having the correlation matrices R(F1)-R(FM) is calculated from the FFT beat signals Bfu1-BfuN generated in step S502. - Then, the
process 900 proceeds to step S906, where it is determined whether the previous correlation matrix Ro(FP) corresponding to the peak frequency FP is stored in thememory 21. - If the previous correlation matrix Ro(FP) is stored in the
memory 21, theprocess 900 proceeds to step S908 directly. - If the previous correlation matrix Ro(FP) is not stored in the
memory 21, theprocess 900 proceeds to step S908 through step S907, where the previous correlation matrix Ro(FP) is generated such that a weighted average of a previous correlation matrix Ro(FP−1) is added to a weighted average of a previous correlation matrix Ro(FP+1). - At step S908, the correlation matrix R(FP) corresponding to the peak frequency FP is extracted from the correlation matrix group RG. The correlation matrix R(FP) is multiplied by the weighting factor (1−K) and the previous correlation matrix Ro(FP) is multiplied by the weighting factor K. Then, the correlation matrix R(FP) multiplied by the weighting factor (1−K) and the previous correlation matrix Ro(FP) multiplied by the weighting factor K are added together to produce the addition correlation matrix U(FP). Therefore, the addition correlation matrix U(FP) is given by:
-
U(FP)=R(FP)·(1−K)+Ro(FP)·K (6) - Then, the
process 900 proceeds to step S909, where the eigenvalue expansion of the addition correlation matrix U(FP) is performed. - Then, the
process 900 proceeds to step S910, where the MUSIC spectrum is calculated based on the eigenvector of the addition correlation matrix U(FP). - Then, the
process 900 proceeds to step S911, where the direction of the object in the increase area is calculated based on the MUSIC spectrum. Because the beat signals Bd1-BdN are processed in the same way as the beat signals Bu1-BuN, the direction of the object in the decrease area is also calculated. - Then, the
process 900 proceeds to step S912, where the correlation matrix group RG is thinned out and stored in thememory 21 as the previous correlation matrix group RoG that is used in a next loop of theprocess 900. - Then, the
process 900 proceeds to step S512. - As described above, in the
process 900, the correlation matrix group RG is thinned out and stored in thememory 21 as the previous correlation matrix group RoG. Thus, thememory 21 can have a small amount of storage capacity in theprocess 900, as compared to in theprocess 800 according to the third embodiment. Even when the previous correlation matrix Ro(FP) is not stored in thememory 21, the previous correlation matrix Ro(FP) is estimated from previous correlation matrices Ro(FP−1), Ro(FP+1). - (Modifications)
- The embodiments described above may be modified in various ways. For example, the
radar apparatus 100 may use the received radar signal received by some of the elements E1-EN of the receivingantenna 12, not each of the elements E1-EN. In such an approach, the amount of calculation can be reduced. - The correlation matrix group RG or the addition correlation matrix group UG may be stored in the
memory 21 after being compressed by a data compression algorithm. Thus, thememory 21 may store the previous information in a form of data containing elements of the correlation matrix group RG or the addition correlation matrix group UG. - The weighting factor K may be a variable. For example, when the received radar signal has considerable instantaneous noise, the weighing factor K may be increased. In such an approach, an influence of the noise can be reduced.
- The algorithm used in the present invention can be applied to various types of the DOA estimation algorithms such as unitary-MUSIC algorithm, ESPRIT algorithm, unitary-ESPRIT algorithm, Capon algorithm, and Beam Former algorithm. In particular, when the unitary-MUSIC or the unitary-ESPRIT is used, only the real part of the matrix is stored in the
memory 21. Therefore, the amount of calculation can be significantly reduced and thememory 21 can have a very small amount of storage capacity. The algorithm used in the present invention also can be applied to a spatial smoothing algorithm. - The previous information (i.e., the previous correlation matrix group RoG, or the previous addition correlation matrix group UG) may be generated in two or more previous loop. For example, the previous information may be generated in two previous loop (i.e., generated 2Ts earlier).
- The transmitting
antenna 11 instead of the receivingantenna 12 may have the elements arranged in the array to generate the beat signals.
Claims (14)
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005224636A JP4857644B2 (en) | 2005-08-02 | 2005-08-02 | Radar equipment |
JP2005-224636 | 2005-08-02 |
Publications (2)
Publication Number | Publication Date |
---|---|
US20090040097A1 true US20090040097A1 (en) | 2009-02-12 |
US7495605B1 US7495605B1 (en) | 2009-02-24 |
Family
ID=37670209
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/497,033 Active 2027-03-15 US7495605B1 (en) | 2005-08-02 | 2006-07-31 | Radar apparatus |
Country Status (3)
Country | Link |
---|---|
US (1) | US7495605B1 (en) |
JP (1) | JP4857644B2 (en) |
DE (1) | DE102006036099B4 (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090021422A1 (en) * | 2007-07-20 | 2009-01-22 | Denso Corporation | Device and method for estimating the number of arrival signals |
US20090309784A1 (en) * | 2008-06-17 | 2009-12-17 | Denso Corporation | Target object detection apparatus for acquiring information concerning target objects based on correlation matrix derived from signal values corresponding to reflected electromagnetic waves |
US20100134343A1 (en) * | 2008-12-01 | 2010-06-03 | Mitsubishi Electric Corporation | Radar device |
FR2950149A1 (en) * | 2009-09-17 | 2011-03-18 | Mbda France | METHOD AND SYSTEM FOR AVOIDING AN INTERCEPTION MACHINE BY AN AIR MOBILE |
US20110102242A1 (en) * | 2009-05-25 | 2011-05-05 | Kabushiki Kaisha Toshiba | Radar apparatus |
US20110156947A1 (en) * | 2009-12-25 | 2011-06-30 | Honda Elesys Co., Ltd. | Electronic scanning radar apparatus, receiving wave direction estimating method, and computer-readable storage media storing receiving wave direction estimation program |
US20110193740A1 (en) * | 2010-02-08 | 2011-08-11 | Honda Elesys Co., Ltd. | Electronic scanning radar apparatus, receiving wave direction estimating method, and computer-readable storage media storing receiving wave direction estimation program |
US20110309971A1 (en) * | 2010-05-19 | 2011-12-22 | Honda Elesys Co., Ltd. | Electronic scanning radar apparatus, receiving wave direction estimating method, and computer-readable storage media storing receiving wave direction estimation program |
US20120038506A1 (en) * | 2010-07-01 | 2012-02-16 | Honda Elesys Co., Ltd | Electronic scanning radar apparatus, received wave direction estimating method, and received wave direction estimation program |
US20120249360A1 (en) * | 2011-02-04 | 2012-10-04 | Honda Elesys Co., Ltd. | Electronic scanning radar apparatus, received wave direction estimating method, and received wave direction estimating program |
US10067227B2 (en) * | 2014-10-06 | 2018-09-04 | Nidec Corporation | Neural network-based radar system |
US10359511B2 (en) * | 2014-12-29 | 2019-07-23 | Sony Corporation | Surveillance apparatus having a radar sensor |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4656124B2 (en) * | 2007-11-09 | 2011-03-23 | 株式会社デンソー | Direction detection device |
US8446312B2 (en) | 2007-12-25 | 2013-05-21 | Honda Elesys Co., Ltd. | Electronic scanning type radar device, estimation method of direction of reception wave, and program estimating direction of reception wave |
JP5102165B2 (en) * | 2008-09-22 | 2012-12-19 | 株式会社デンソー | Radar equipment |
JP4926155B2 (en) * | 2008-11-04 | 2012-05-09 | 三菱電機株式会社 | Radar equipment |
JP5717950B2 (en) * | 2009-02-05 | 2015-05-13 | アズビル株式会社 | Model function processing apparatus and method |
JP2010237087A (en) * | 2009-03-31 | 2010-10-21 | Hitachi Automotive Systems Ltd | Radar apparatus and method for measuring radio wave arrival direction using the same |
JP5494567B2 (en) * | 2011-05-17 | 2014-05-14 | 株式会社デンソー | Radar apparatus, inspection system, and inspection method |
JP2013152239A (en) * | 2013-03-29 | 2013-08-08 | Honda Elesys Co Ltd | Electronic scanning radar device, control method of the same and program |
JP6278961B2 (en) * | 2013-06-13 | 2018-02-14 | パナソニック株式会社 | Radar equipment |
JP5925264B2 (en) * | 2014-09-10 | 2016-05-25 | 三菱電機株式会社 | Radar equipment |
JP6887571B2 (en) * | 2018-09-07 | 2021-06-16 | 三菱電機株式会社 | Signal processing circuits, radar devices, signal processing methods and signal processing programs |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6121917A (en) * | 1997-10-31 | 2000-09-19 | Toyota Jidosha Kabushiki Kaisha | FM-CW radar |
US20010015698A1 (en) * | 2000-02-07 | 2001-08-23 | Toyota Jidosha Kabushiki Kaisha | Radar apparatus |
US20030052813A1 (en) * | 2001-09-17 | 2003-03-20 | Kazuma Natsume | Radar designed to minimize error in detecting target |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000147083A (en) * | 1998-11-06 | 2000-05-26 | Ntt Mobil Communication Network Inc | Method and apparatus for measuring arrival angle |
JP3489114B2 (en) * | 2000-02-07 | 2004-01-19 | トヨタ自動車株式会社 | FM-CW radar device |
DE10105240A1 (en) | 2000-02-07 | 2001-08-16 | Toyota Motor Co Ltd | Radar apparatus for cars, has pair of bearing detectors which carry out signal processing, in different manners |
JP2001337148A (en) * | 2000-05-26 | 2001-12-07 | Mitsubishi Electric Corp | Extrapolation device for electromagnetic wave arrival direction |
JP3663623B2 (en) * | 2000-06-29 | 2005-06-22 | トヨタ自動車株式会社 | Radar equipment |
JP2002107440A (en) * | 2000-09-29 | 2002-04-10 | Rikogaku Shinkokai | Method and device for detecting direction |
JP4577742B2 (en) * | 2000-11-14 | 2010-11-10 | 住友電気工業株式会社 | Radio wave arrival direction estimation device |
JP3821688B2 (en) * | 2001-10-30 | 2006-09-13 | 三菱電機株式会社 | Radar equipment |
JP2003318792A (en) * | 2002-04-26 | 2003-11-07 | Matsushita Electric Ind Co Ltd | Apparatus for estimating arrival wave direction, and receiver provided with the same |
JP2006145251A (en) * | 2004-11-16 | 2006-06-08 | Denso Corp | Radio wave arrival direction prediction system |
JP4829517B2 (en) * | 2005-03-31 | 2011-12-07 | 株式会社デンソーアイティーラボラトリ | Radar signal processing device |
-
2005
- 2005-08-02 JP JP2005224636A patent/JP4857644B2/en active Active
-
2006
- 2006-07-31 US US11/497,033 patent/US7495605B1/en active Active
- 2006-08-02 DE DE102006036099.0A patent/DE102006036099B4/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6121917A (en) * | 1997-10-31 | 2000-09-19 | Toyota Jidosha Kabushiki Kaisha | FM-CW radar |
US20010015698A1 (en) * | 2000-02-07 | 2001-08-23 | Toyota Jidosha Kabushiki Kaisha | Radar apparatus |
US20030052813A1 (en) * | 2001-09-17 | 2003-03-20 | Kazuma Natsume | Radar designed to minimize error in detecting target |
US6646589B2 (en) * | 2001-09-17 | 2003-11-11 | Denso Corporation | Radar designed to minimize error in detecting target |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090021422A1 (en) * | 2007-07-20 | 2009-01-22 | Denso Corporation | Device and method for estimating the number of arrival signals |
US20090309784A1 (en) * | 2008-06-17 | 2009-12-17 | Denso Corporation | Target object detection apparatus for acquiring information concerning target objects based on correlation matrix derived from signal values corresponding to reflected electromagnetic waves |
US7724181B2 (en) * | 2008-06-17 | 2010-05-25 | Denso Corporation | Target object detection apparatus for acquiring information concerning target objects based on correlation matrix derived from signal values corresponding to reflected electromagnetic waves |
US20100134343A1 (en) * | 2008-12-01 | 2010-06-03 | Mitsubishi Electric Corporation | Radar device |
US8102309B2 (en) * | 2008-12-01 | 2012-01-24 | Mitsubishi Electric Corporation | Radar device |
US20110102242A1 (en) * | 2009-05-25 | 2011-05-05 | Kabushiki Kaisha Toshiba | Radar apparatus |
FR2950149A1 (en) * | 2009-09-17 | 2011-03-18 | Mbda France | METHOD AND SYSTEM FOR AVOIDING AN INTERCEPTION MACHINE BY AN AIR MOBILE |
WO2011033190A1 (en) * | 2009-09-17 | 2011-03-24 | Mbda France | Method and system for avoiding an intercepting vehicle by an airborne moving body |
EP2302411A1 (en) * | 2009-09-17 | 2011-03-30 | MBDA France | Method and system for an airborne mobile to avoid an intercepting vehicle |
US8718921B2 (en) | 2009-09-17 | 2014-05-06 | Mbda France | Method and system for avoiding an intercepting vehicle by an airborne moving body |
US8441395B2 (en) * | 2009-12-25 | 2013-05-14 | Honda Elesys Co., Ltd. | Electronic scanning radar apparatus, receiving wave direction estimating method, and computer-readable storage media storing receiving wave direction estimation program |
US20110156947A1 (en) * | 2009-12-25 | 2011-06-30 | Honda Elesys Co., Ltd. | Electronic scanning radar apparatus, receiving wave direction estimating method, and computer-readable storage media storing receiving wave direction estimation program |
US8466830B2 (en) * | 2010-02-08 | 2013-06-18 | Honda Elesys Co., Ltd. | Electronic scanning radar apparatus, receiving wave direction estimating method, and computer-readable storage media storing receiving wave direction estimation program |
US20110193740A1 (en) * | 2010-02-08 | 2011-08-11 | Honda Elesys Co., Ltd. | Electronic scanning radar apparatus, receiving wave direction estimating method, and computer-readable storage media storing receiving wave direction estimation program |
US20110309971A1 (en) * | 2010-05-19 | 2011-12-22 | Honda Elesys Co., Ltd. | Electronic scanning radar apparatus, receiving wave direction estimating method, and computer-readable storage media storing receiving wave direction estimation program |
US8552907B2 (en) * | 2010-05-19 | 2013-10-08 | Honda Elesys Co., Ltd. | Electronic scanning radar apparatus, receiving wave direction estimating method, and computer-readable storage media storing receiving wave direction estimation program |
US20120038506A1 (en) * | 2010-07-01 | 2012-02-16 | Honda Elesys Co., Ltd | Electronic scanning radar apparatus, received wave direction estimating method, and received wave direction estimation program |
US8581777B2 (en) * | 2010-07-01 | 2013-11-12 | Honda Elesys Co., Ltd. | Electronic scanning radar apparatus, received wave direction estimating method, and received wave direction estimation program |
US20120249360A1 (en) * | 2011-02-04 | 2012-10-04 | Honda Elesys Co., Ltd. | Electronic scanning radar apparatus, received wave direction estimating method, and received wave direction estimating program |
US8907840B2 (en) * | 2011-02-04 | 2014-12-09 | Honda Elesys Co., Ltd. | Electronic scanning radar apparatus, received wave direction estimating method, and received wave direction estimating program |
US10067227B2 (en) * | 2014-10-06 | 2018-09-04 | Nidec Corporation | Neural network-based radar system |
US10359511B2 (en) * | 2014-12-29 | 2019-07-23 | Sony Corporation | Surveillance apparatus having a radar sensor |
Also Published As
Publication number | Publication date |
---|---|
JP2007040806A (en) | 2007-02-15 |
DE102006036099B4 (en) | 2023-03-02 |
JP4857644B2 (en) | 2012-01-18 |
US7495605B1 (en) | 2009-02-24 |
DE102006036099A1 (en) | 2007-02-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7495605B1 (en) | Radar apparatus | |
JP5114187B2 (en) | Electronic scanning radar apparatus, received wave direction estimation method, and received wave direction estimation program | |
US8730096B2 (en) | Electronic scanning type radar device and method for estimating direction of reception wave | |
US7474252B2 (en) | Radar device | |
JP5062225B2 (en) | Target detection device | |
US7567201B2 (en) | Vehicle-installation direction detection apparatus enabling accurate detection of target body directions irrespective of vehicle speed | |
US7688255B2 (en) | Electronic scanning radar apparatus | |
US11275172B2 (en) | Target detection device | |
JP2009162688A (en) | Electronic scanning radar device, and received wave direction estimation method and program | |
WO2006085352A1 (en) | Target detecting device | |
JP3821688B2 (en) | Radar equipment | |
JP5114217B2 (en) | Electronic scanning radar apparatus, received wave direction estimation method, and received wave direction estimation program | |
JP2013152239A (en) | Electronic scanning radar device, control method of the same and program | |
JP5779370B2 (en) | Electronic scanning radar apparatus, received wave direction estimation method, and program | |
JP2006145251A (en) | Radio wave arrival direction prediction system | |
JP3505441B2 (en) | Peak frequency calculation method in FFT signal processing | |
JP7056212B2 (en) | Direction estimation method and equipment | |
KR102099388B1 (en) | Method of estimating direction of arrival of radar signal based on antenna array extrapolation and apparatus for the same | |
JP2001091639A (en) | Fm-cw radar device | |
JP2017211241A (en) | Incoming wave angle estimation method and incoming wave angle estimation device | |
JPH1031065A (en) | Fm-cw radar | |
WO2022249552A1 (en) | Information processing device and information processing method | |
JP6688977B2 (en) | Radar equipment | |
JP4564407B2 (en) | Antenna device | |
WO2021210480A1 (en) | Communication device and communication method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: DENSO CORPORATION, JAPAN Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:SAKAMOTO, MAI;NATSUME, KAZUMA;MIYAKE, YASUYUKI;REEL/FRAME:018115/0362;SIGNING DATES FROM 20060717 TO 20060718 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
FEPP | Fee payment procedure |
Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
FEPP | Fee payment procedure |
Free format text: PAYER NUMBER DE-ASSIGNED (ORIGINAL EVENT CODE: RMPN); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
FPAY | Fee payment |
Year of fee payment: 8 |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 12TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1553); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY Year of fee payment: 12 |