WO2019163632A1 - 方位推定方法及び方位推定装置 - Google Patents

方位推定方法及び方位推定装置 Download PDF

Info

Publication number
WO2019163632A1
WO2019163632A1 PCT/JP2019/005315 JP2019005315W WO2019163632A1 WO 2019163632 A1 WO2019163632 A1 WO 2019163632A1 JP 2019005315 W JP2019005315 W JP 2019005315W WO 2019163632 A1 WO2019163632 A1 WO 2019163632A1
Authority
WO
WIPO (PCT)
Prior art keywords
matrix
unit
bin
peak
transmission
Prior art date
Application number
PCT/JP2019/005315
Other languages
English (en)
French (fr)
Inventor
尭之 北村
Original Assignee
株式会社デンソー
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 株式会社デンソー filed Critical 株式会社デンソー
Publication of WO2019163632A1 publication Critical patent/WO2019163632A1/ja
Priority to US16/996,622 priority Critical patent/US11125871B2/en

Links

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
    • 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/50—Systems of measurement based on relative movement of target
    • G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S13/583—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets
    • G01S13/584—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets adapted for simultaneous range and velocity measurements
    • 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/66—Radar-tracking systems; Analogous systems
    • G01S13/72—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
    • 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/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/93—Radar or analogous systems specially adapted for specific applications for anti-collision purposes
    • G01S13/931—Radar or analogous systems specially adapted for specific applications for anti-collision purposes of land vehicles
    • 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/10—Systems for measuring distance only using transmission of interrupted, pulse modulated waves
    • G01S13/18—Systems for measuring distance only using transmission of interrupted, pulse modulated waves wherein range gates are used
    • 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
    • G01S13/44—Monopulse radar, i.e. simultaneous lobing
    • G01S13/4418—Monopulse radar, i.e. simultaneous lobing with means for eliminating radar-dependent errors in angle measurements, e.g. multipath effects
    • 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
    • G01S13/44—Monopulse radar, i.e. simultaneous lobing
    • G01S13/4427—Monopulse radar, i.e. simultaneous lobing with means for eliminating the target-dependent errors in angle measurements, e.g. glint, scintillation effects
    • 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/50—Systems of measurement based on relative movement of target
    • G01S13/52—Discriminating between fixed and moving objects or between objects moving at different speeds
    • G01S13/522—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves
    • G01S13/524—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi
    • G01S13/53—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi performing filtering on a single spectral line and associated with one or more range gates with a phase detector or a frequency mixer to extract the Doppler information, e.g. pulse Doppler radar
    • 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/50—Systems of measurement based on relative movement of target
    • G01S13/52—Discriminating between fixed and moving objects or between objects moving at different speeds
    • G01S13/522—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves
    • G01S13/524—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi
    • G01S13/53—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi performing filtering on a single spectral line and associated with one or more range gates with a phase detector or a frequency mixer to extract the Doppler information, e.g. pulse Doppler radar
    • G01S13/532—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi performing filtering on a single spectral line and associated with one or more range gates with a phase detector or a frequency mixer to extract the Doppler information, e.g. pulse Doppler radar using a bank of range gates or a memory matrix
    • 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/50—Systems of measurement based on relative movement of target
    • G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S13/60—Velocity or trajectory determination systems; Sense-of-movement determination systems wherein the transmitter and receiver are mounted on the moving object, e.g. for determining ground speed, drift angle, ground track
    • 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/50—Systems of measurement based on relative movement of target
    • G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S13/64—Velocity measuring systems using range gates

Definitions

  • the present disclosure relates to a technique for executing a direction estimation calculation using a correlation matrix.
  • a high-resolution arrival direction estimation technique such as MUSIC or ESPRIT is known.
  • MUSIC is an abbreviation for Multiple Signal Classification
  • ESPRIT is an abbreviation for Estimate of Signal Parameters via Invariance Techniques.
  • these methods are referred to as high-resolution methods.
  • Patent Document 1 describes using a spatial averaging method to suppress the cross-correlation component.
  • antenna arrays arranged at equal intervals hereinafter referred to as equal intervals array
  • equal intervals array antenna arrays arranged at equal intervals
  • a correlation matrix is generated for each sub-array.
  • Average the correlation matrix In addition to this, a time averaging method is also known that averages a plurality of correlation matrices obtained by a plurality of measurements performed in a short time. In any case, it is necessary to secure the number of snapshots, which is the number of correlation matrices to be averaged.
  • the spatial averaging method has a problem that the receiving array antenna can only be used if it is an equally spaced array that can be divided into a plurality of small subarrays having the same shape.
  • FCM is an abbreviation for Fast-Chirp-Modulation. That is, in the FCM method, a chirp whose frequency changes with time is repeatedly transmitted, and the target distance is separated by executing a fast Fourier transform process for each chirp for each channel. Further, for each channel, the target frequency is separated by collecting the values of the same frequency bin of the distance spectrum from all the chirps and executing the FFT process. As described above, in the FCM method, it is necessary to individually perform FFT processing on all the chirps, and thus it is not possible to use a time averaging method that averages correlation matrices calculated respectively for a plurality of chirps.
  • One aspect of the present disclosure is to improve the accuracy of azimuth separation and the accuracy of azimuth estimation of a plurality of targets by providing a method for suppressing cross-correlation of a correlation matrix different from the conventional one.
  • the orientation estimation method includes an analysis step, an extraction step, a center generation step, a periphery generation step, an integration step, and an estimation step.
  • a frequency spectrum is generated for each transmission / reception channel by performing frequency analysis on signals obtained from a plurality of transmission / reception channels.
  • a peak bin that is a frequency bin that becomes a peak is extracted from the frequency spectrum generated in the analysis step.
  • a center matrix which is a correlation matrix, is generated using the same peak bin values collected from all transmission / reception channels for each peak bin extracted in the extraction step.
  • a peripheral matrix that is a correlation matrix is generated using the same peripheral bin values collected from all transmission / reception channels for each peripheral bin for each peak bin.
  • the peripheral bins are one or more frequency bins that exist within a preset range around the peak bin.
  • an integration matrix that is a correlation matrix is generated for each peak bin by weighted addition of the center matrix generated in the center generation step and one or more peripheral matrices generated in the periphery generation step. To do.
  • the direction estimation calculation is performed using the integration matrix generated in the integration step.
  • the orientation estimation apparatus includes an analysis unit, an extraction unit, a center generation unit, a peripheral generation unit, an integration unit, and an estimation unit.
  • Each unit of the azimuth estimation device performs the same processing as each step of the azimuth estimation method, that is, an analysis step, an extraction step, a center generation step, a periphery generation step, an integration step, and an estimation step.
  • the cross-correlation of the correlation matrix The component can be suppressed.
  • estimation accuracy can be improved by azimuth estimation calculation using a correlation matrix such as MUSIC.
  • An azimuth estimation apparatus 1 shown in FIG. 1 is a so-called millimeter wave radar apparatus that is mounted on a vehicle and detects various targets existing around the vehicle using millimeter waves.
  • the direction estimation apparatus 1 includes a signal processing unit 4.
  • the direction estimation apparatus 1 may further include an antenna unit 2 and a transmission / reception unit 3.
  • the antenna unit 2 includes a transmission antenna unit 21 and a reception antenna unit 22.
  • the transmission antenna unit 21 has one or more antennas used for transmission.
  • the reception antenna unit 22 has one or more antennas used for reception. However, at least one of the transmission antenna unit 21 and the reception antenna unit 22 has a plurality of antennas.
  • the antenna unit 2 is configured to have a plurality of transmission / reception channels by using a combination of an antenna belonging to the transmission antenna unit 21 and an antenna belonging to the reception antenna unit 22 as a transmission / reception channel.
  • L1 is an integer of 2 or more.
  • the arrangement interval of the antennas belonging to the transmission antenna unit 21 and the arrangement interval of the antennas belonging to the reception antenna unit 22 may both be equal intervals or unequal intervals.
  • the transmission / reception unit 3 includes an analog circuit that transmits a transmission signal modulated by the FCM method via the transmission antenna unit 21.
  • the transmission / reception unit 3 includes an analog circuit that generates a beat signal representing a frequency difference between the transmission signal and the reception signal for each of a plurality of transmission / reception channels based on the reception signal received via the reception antenna unit 22.
  • the transmission / reception unit 3 includes an AD converter that performs analog-digital conversion on the beat signal generated for each transmission / reception channel, and outputs AD conversion data that is an output of the AD converter to the signal processing unit 4.
  • the transmission / reception unit 3 sets a measurement cycle Tf as one frame, and chirps a frequency that changes linearly with respect to time during the first measurement period Tm of each frame. Only the set number of repetitions L2 is continuously transmitted.
  • the measurement period Tf may be 50 ms, for example.
  • the measurement period Tm may be 10 ms, for example.
  • the number of repetitions L2 may be, for example, 256 times.
  • the signal processing unit 4 includes a microcomputer having a CPU 41 and a semiconductor memory (hereinafter, memory 42) such as a RAM or a ROM.
  • the signal processing unit 4 includes a digital signal processor (hereinafter referred to as DSP) 43 that implements fast Fourier transform (hereinafter referred to as FFT) processing by hardware.
  • DSP 43 is a general-purpose hardware calculator. However, the DSP 43 is not necessarily a general-purpose one.
  • Each function other than the DSP 43 in the signal processing unit 4 is realized by the CPU 41 executing a program stored in a non-transitional physical recording medium.
  • the memory 42 corresponds to a non-transitional tangible recording medium that stores a program. Also, by executing this program, a method corresponding to the program is executed.
  • the signal processing unit 4 may include one microcomputer or a plurality of microcomputers.
  • the method of realizing functions other than the DSP 43 of the signal processing unit 4 is not limited to software, and part or all of the functions may be realized using one or a plurality of hardware.
  • the function when the function is realized by an electronic circuit that is hardware, the electronic circuit may be realized by a digital circuit, an analog circuit, or a combination thereof.
  • the function of the DSP 43 may be realized by software.
  • This process is repeatedly executed every processing cycle of the measurement cycle Tf.
  • the CPU 41 acquires AD conversion data obtained by AD-converting the beat signal from each of the L1 transmission / reception channels in S110.
  • the CPU 41 selects, as a target channel, one of the L1 transmission / reception channels that has not been subjected to the processing of S130 to S160 described later.
  • the CPU 41 uses the DSP 43 for the AD conversion data of the target channel, and performs FFT processing for each of the L2 chirps to generate L2 FFT spectra.
  • data to be processed is multiplied by a preset window function.
  • the FFT spectrum generated here is referred to as a distance spectrum in the following because the frequency is associated with the distance to the target.
  • a frequency bin in the distance spectrum is called a distance bin.
  • the number of distance bins is Nd.
  • the CPU 41 uses the DSP 43 to perform FFT processing on M data belonging to the same distance bin collected from the L2 distance spectrums generated in S130.
  • Nd FFT spectra are generated.
  • data to be processed is multiplied by a preset window function.
  • the FFT spectrum generated here is referred to as a speed spectrum below because the frequency is associated with the relative speed between the vehicle and the target.
  • a frequency bin in the velocity spectrum is called a velocity bin.
  • the number of velocity bins is Nv.
  • a two-dimensional FFT spectrum having Nd ⁇ Nv frequency bins having a distance and a relative velocity as two axes, represented by Nd velocity spectra generated for each distance bin will be described below. This is called a two-dimensional spectrum.
  • a bin of a two-dimensional spectrum is called a two-dimensional bin.
  • the CPU 41 determines whether or not the above-described processing of S130 to S140 has been executed for all of the L1 transmission / reception channels. If there is an unprocessed transmission / reception channel, the CPU 41 returns to S120, and if the processing has been executed for all transmission / reception channels, the CPU 41 proceeds to S160.
  • the CPU 41 extracts, as a peak bin, a two-dimensional bin in which a peak is detected in the L1 two-dimensional spectrum generated for each transmission / reception channel.
  • the number of detected peak bins is P.
  • the CPU 41 estimates a two-dimensional bin in which the peak bin is detected in the current processing cycle based on the peak bin detected in the previous processing cycle. Furthermore, if the distance on the two-dimensional spectrum between the detected peak bin and the estimated peak bin is within the threshold range, the CPU 41 performs correlation as a peak bin based on the same object. This process is a technique called tracking.
  • the CPU 41 correlates each channel of the transmission and reception channels based on L1 data belonging to the same peak bin collected from each of the L1 two-dimensional spectrums. Generate a correlation matrix representing the relationship. This correlation matrix is called a central matrix. The correlation is expressed by a coefficient representing the degree of similarity of signals obtained in each channel. The correlation is weaker as the absolute value of the coefficient approaches 0, and the correlation is stronger as the absolute value of the coefficient approaches 1. .
  • the CPU 41 extracts, for each of the P peak bins, a plurality of two-dimensional bins existing around the peak bin of interest as peripheral bins, and for each of the extracted peripheral bins, the central matrix described in S180 A correlation matrix is generated by a similar method. This correlation matrix is called a marginal matrix.
  • the peripheral bins may include a total of four bins, for example, two bins adjacent in the velocity bin direction of the peak bin and two bins adjacent in the distance bin direction as shown in the middle of FIG. Further, as shown in the lower part of FIG. 5, it may further include a total of eight bins including four bins located at the diagonal upper right, diagonal lower right, diagonal upper left, and diagonal lower left of the peak bin. . Not limited to this, the peripheral bins only need to include one or more frequency bins that exist within a preset range around the peak bin.
  • the CPU 41 In S200, the CPU 41 generates one averaged correlation matrix by weighting and adding the central matrix generated in S180 and the peripheral matrix generated in S190 for each peak bin. This correlation matrix is called an integrated matrix. As a result, P integration matrices corresponding to each of the P peak bins are generated. This integrated matrix is obtained by suppressing the value of the off-diagonal component of the integrated matrix representing the cross-correlation between channels as compared to the individual central matrix and the peripheral matrix. Further, the CPU 41 stores the calculated integration matrix in a preset storage area on the memory 42 in association with the two-dimensional bin that specifies the peak bin.
  • the CPU 41 further suppresses cross-correlation between channels using the time average method for each of the P integrated matrices. Specifically, as shown in FIG. 6, in the tracking process of S170, the peak bins for the past S cycles that are associated with the peak bin of interest are extracted, and the integrated matrix generated for each extracted peak bin is acquired. . S is an integer of 1 or more. Furthermore, a time-averaged integrated matrix is generated by weighted addition of the acquired integrated matrix and the integrated matrix calculated in S190.
  • the CPU 41 executes a direction estimation calculation such as MUSIC for each peak bin using the time-averaged integrated matrix generated in S210.
  • a direction estimation calculation such as MUSIC for each peak bin using the time-averaged integrated matrix generated in S210.
  • a plurality of directions may be calculated from one peak bin. Therefore, the number Q ( ⁇ P) obtained by summing the number of azimuths detected in the azimuth estimation calculation for each peak bin for all peak bins is the number of targets.
  • the CPU 41 In S230, the CPU 41 generates target information for each of the Q targets, and ends this process.
  • the target information is estimated by using a distance calculated from a distance bin corresponding to the peak bin of interest, a relative speed calculated from a velocity bin corresponding to the peak bin of interest, and a correlation matrix of the peak bin of interest. An azimuth angle is included.
  • This target information is distributed to each part of the vehicle via an in-vehicle LAN (not shown) and used for various driving support controls.
  • S130 and S140 correspond to an analysis step and an analysis unit.
  • S160 corresponds to a rounding step and an S extraction unit.
  • S180 corresponds to a center generation step and a center generation unit.
  • S190 corresponds to a peripheral generation step and a peripheral generation unit.
  • S200 corresponds to an integration step and an integration unit.
  • S220 corresponds to an estimation step and an estimation unit.
  • S200 also corresponds to a storage unit.
  • S170 corresponds to the tracking unit.
  • S210 corresponds to the time average part.
  • a multi-point reflection model is applied to the target.
  • a multipoint reflection model is effective because a plurality of reflections are returned from one target. That is, the reflection peak observed for one target is a combination of N reflection peaks from N reflection points.
  • reflected waves from N reflection points on the same target identified by i are incident on M antennas from the same direction, and a certain amount of phase rotation is performed between adjacent antennas. ⁇ i is generated.
  • the N reflection points are slightly different in distance, the N reflection peaks based on the N reflection points are centered on a certain distance bin determined for each target on the distance spectrum as shown in FIG. Stochastic variation.
  • Equation (1) is obtained by adding a vector in which window function filter responses H ( ⁇ ij) of signals from the j-th reflection point regarding the i-th target are sequentially arranged for the M antennas for the N reflection points. Means.
  • the received signal vector xi ( ⁇ (k)) at the angular frequency ⁇ (k) in the i-th target is expressed by equation (2).
  • ⁇ is the phase change amount for each angular frequency step ⁇ due to the filter response of the window function.
  • ⁇ is a constant amount.
  • Equation (3) represents a composite received signal x ( ⁇ (k)) obtained by synthesizing peak signals from two targets.
  • the cross-correlation term between the targets shown in the third and following terms on the right side of equation (4) can be positive or negative. Get closer.
  • the autocorrelation terms of each target shown in the first and second terms on the right side in the equation (4) always have a positive amplitude, and when added, the value increases unilaterally. Thereby, as a result, the cross-correlation between the incoming waves of the correlation matrix is reduced.
  • W (s) is the weight of the integration matrix generated before s cycles. Also in this case, the cross-correlation term is suppressed as in the case of averaging using the correlation matrix of the peripheral bins.
  • the azimuth estimation apparatus 1 secures the number of snapshots of the correlation matrix using the peripheral bins adjacent to the peak bin, and averages them to suppress the cross-correlation component of the correlation matrix. Therefore, the azimuth estimation apparatus 1 suppresses the cross-correlation component of the correlation matrix even when modulation that cannot apply the time average method using the result of FFT processing for each chirp, such as the FCM method, is used. it can. Further, according to the azimuth estimation apparatus 1, since the number of snapshots can be secured without using the spatial averaging method, the cross-correlation component of the correlation matrix can be suppressed even when an unequally spaced antenna array is used. As a result, according to the azimuth estimation apparatus 1, it is possible to improve the accuracy of azimuth separation by azimuth estimation calculation using a correlation matrix such as MUSIC and the accuracy of azimuth estimation.
  • FIG. 10 shows the result of calculating by simulation the detection rate of the target simulating two vehicles running in parallel while changing the S / N of the target using MUSIC as the direction estimation calculation.
  • the correlation matrix is generated with only peak bins
  • the cross-correlation component of the correlation matrix is suppressed using four peripheral bins
  • the cross-correlation component of the correlation matrix is suppressed using eight peripheral bins Shows the case.
  • FIG. 10 shows that the cross-correlation component can be further suppressed by increasing the number of peripheral pins to be used.
  • the range of peripheral bins to be used may be set in consideration of this because it may deteriorate due to the influence of another peak that does not overlap.
  • the azimuth estimation apparatus 1 extracts peak bins based on the same target detected during the past S cycle by tracking peak bins, and performs time averaging of the integration matrix calculated for the extracted peak bins. .
  • the cross-correlation component of the correlation matrix used for the distance estimation calculation is suppressed. That is, when the change in the azimuth angle of the target during the measurement period Tf is sufficiently small relative to the distance resolution in the azimuth estimation calculation, the cross-correlation component of the correlation matrix is further suppressed by further applying this method. Is done. As a result, it is possible to further improve the accuracy of orientation separation and the accuracy of orientation estimation by the orientation estimation calculation.
  • the present invention may be applied when generating a correlation matrix for peak bins extracted from a distance spectrum.
  • weighted addition is used as an operation for suppressing the cross-correlation component of the correlation matrix, but simple addition in which all weights are the same may be used.
  • a plurality of functions of one constituent element in the above embodiment may be realized by a plurality of constituent elements, or a single function of one constituent element may be realized by a plurality of constituent elements. . Further, a plurality of functions possessed by a plurality of constituent elements may be realized by one constituent element, or one function realized by a plurality of constituent elements may be realized by one constituent element. Moreover, you may abbreviate
  • the present disclosure can be implemented in various forms.

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

中心生成ステップ（Ｓ１８０）では、抽出ステップ（Ｓ１６０）にて抽出されたピークビン毎に、全ての送受信チャネルから集めた同一のピークビンの値を用いて、相関行列である中心行列を生成する。周辺生成ステップ（Ｓ１９０）では、ピークビンのそれぞれについて周辺ビン毎に、全ての送受信チャネルから集めた同一の周辺ビンの値を用いて、相関行列である周辺行列を生成する。統合ステップ（Ｓ２００）では、ピークビン毎に、中心生成ステップにて生成された中心行列と、周辺生成ステップにて生成された１つ以上の周辺行列とを重み付け加算して統合行列を生成する。推定ステップ（Ｓ２２０）では、統合ステップにて生成された統合行列を用いて方位推定演算を実行する。

Description

方位推定方法及び方位推定装置 関連出願の相互参照
　本国際出願は、２０１８年２月２０日に日本国特許庁に出願された日本国特許出願第２０１８－０２７９０７号に基づく優先権を主張するものであり、日本国特許出願第２０１８－０２７９０７号の全内容を本国際出願に参照により援用する。
　本開示は、相関行列を用いて方位推定演算を実行する技術に関する。
　レーダ装置において、複数のアンテナから得られる受信信号から相関行列を生成し、その相関行列を用いて、ターゲットの方位を演算する手法として、ＭＵＳＩＣ及びＥＳＰＲＩＴ等の高分解能な到来方向推定手法が知られている。ＭＵＳＩＣは、Multiple Signal Classificationの略であり、ＥＳＰＲＩＴは、Estimation of Signal Parameters via Invariance Techniquesの略である。以下、これらの手法を高分解能手法と言う。
　高分解能手法において、複数ターゲットの方位分離を行うには、相関行列における相互相関成分である、相関行列の非対角項の値を抑圧する必要がある。下記特許文献１には、相互相関成分を抑圧するために空間平均法を用いることが記載されている。空間平均法では、等間隔に配置されたアンテナアレー（以下、等間隔アレー）を、同一相似形状を持つ複数の小さなサブアレーに分割し、サブアレー毎に相関行列を生成して、生成された複数の相関行列を平均化する。これ以外にも、短時間の間に行われる複数回の計測で得られた複数の相関行列を平均化する時間平均法も知られている。いずれにしても、平均化の対象となる相関行列の数であるスナップショット数を確保する必要がある。
特開２０１７－９０２２９号公報
　しかしながら、開示者の詳細な検討の結果、従来技術には以下の課題が見出された。
　即ち、空間平均法では、受信アレーアンテナが、同一相似形状を持つ複数の小さなサブアレーに分割できる等間隔アレーでなければ用いることができないという課題があった。
　また、レーダの変調方式として、ＦＣＭ方式を用いた場合、時間平均法を用いることができないという課題もあった。ＦＣＭは、Fast-Chirp Modulationの略である。即ち、ＦＣＭ方式では、時間と共に周波数が変化するチャープを繰り返し送信し、各チャネルについて、チャープ毎に高速フーリエ変換処理を実行することによってターゲットの距離を分離する。更に、チャネル毎に、全てのチャープから距離スペクトラムの同一周波数ビンの値を集めてＦＦＴ処理を実行することによってターゲットの速度を分離する。このように、ＦＣＭ方式では、全てのチャープを個別にＦＦＴ処理する必要があるため、複数のチャープにてそれぞれ算出される相関行列を平均化する時間平均法を用いることができなかった。
　本開示の１つの局面は、従来とは異なる相関行列の相互相関を抑圧する手法を提供することによって、複数ターゲットの方位分離の精度と方位推定の精度とを向上させることにある。
　本開示の一態様による方位推定方法は、解析ステップと、抽出ステップと、中心生成ステップと、周辺生成ステップと、統合ステップと、推定ステップとを備える。
　解析ステップでは、複数の送受信チャネルから得られる信号を周波数解析することで、送受信チャネル毎に周波数スペクトラムを生成する。抽出ステップでは、解析ステップにて生成された周波数スペクトラムからピークとなる周波数ビンであるピークビンを抽出する。中心生成ステップでは、抽出ステップにて抽出されたピークビン毎に、全ての送受信チャネルから集めた同一のピークビンの値を用いて、相関行列である中心行列を生成する。周辺生成ステップでは、ピークビンのそれぞれについて周辺ビン毎に、全ての送受信チャネルから集めた同一の周辺ビンの値を用いて、相関行列である周辺行列を生成する。周辺ビンは、ピークビンを中心として予め設定された範囲内に存在する一つ以上の周波数ビンである。統合ステップでは、ピークビン毎に、中心生成ステップにて生成された中心行列と、周辺生成ステップにて生成された１つ以上の周辺行列とを重み付け加算することで、相関行列である統合行列を生成する。推定ステップでは、統合ステップにて生成された統合行列を用いて方位推定演算を実行する。
　本開示の別の態様による方位推定装置は、解析部と、抽出部と、中心生成部と、周辺生成部と、統合部と、推定部とを備える。方位推定装置の各部は、それぞれ、方位推定方法の各ステップ、即ち、解析ステップと、抽出ステップと、中心生成ステップと、周辺生成ステップと、統合ステップと、推定ステップと同様の処理を実行する。
　このような構成によれば、不等間隔に配置されたアンテナアレーが用いられる場合、又は、ＦＣＭ方式等のようにスナップショット数を１より大きくできない変調が用いられる場合でも、相関行列の相互相関成分を抑圧することができる。その結果、ＭＵＳＩＣなど、相関行列を利用した方位推定演算により推定精度を向上させることができる。
方位推定装置の構成を示すブロック図である。 送信信号の変調パターンを示す説明図である。 ターゲット検出処理のフローチャートである。 距離スペクトラムと速度スペクトラムとの関係を示す説明図である。 ２次元スペクトラム上のピークビンと周辺ビンとを例示する説明図である。 時間平均の対象となるピークビンを示す説明図である。 ２次元スペクトラム上のピークビンと一つのピークビンから求められるＭＵＳＩＣスペクトラムとを例示する説明図である。 ターゲット毎に反射波に基づくアンテナ間の位相回転が異なることを示す説明図である。 距離スペクトラム上で検出される個々の反射点に基づくピークと、原理の説明に使用するパラメータとを示す説明図である。 周辺ビンにより相関行列の相互相関成分を抑圧することによる効果を示すグラフである。
　以下、図面を参照しながら、本開示の実施形態を説明する。
　［１．構成］
　図１に示す方位推定装置１は、車両に搭載され、ミリ波を使用して車両の周囲に存在する様々なターゲットを検出する、いわゆるミリ波レーダ装置である。方位推定装置１は、信号処理部４を備える、方位推定装置１は、更に、アンテナ部２と、送受信部３と、を備えてもよい。
　アンテナ部２は、送信アンテナ部２１と、受信アンテナ部２２とを備える。送信アンテナ部２１は、送信に使用される一つ以上のアンテナを有する。受信アンテナ部２２は、受信に使用される一つ以上のアンテナを有する。但し、送信アンテナ部２１と受信アンテナ部２２とのうち、少なくとも一方は複数のアンテナを有する。
　つまり、送信アンテナ部２１に属するアンテナと、受信アンテナ部２２に属するアンテナとの組み合わせを送受信チャネルとして、アンテナ部２は、複数の送受信チャネルを有するように構成される。ここでは、送信アンテナ部２１に属するアンテナを１個、受信アンテナ部２２に属するアンテナをＬ個備え、Ｌ１個の送受信チャネルを有するものとする。Ｌ１は２以上の整数である。なお、送信アンテナ部２１に属するアンテナの配置間隔及び受信アンテナ部２２に属するアンテナの配置間隔は、いずれも、等間隔であってもよいし、不等間隔であってもよい。
　送受信部３は、ＦＣＭ方式で変調された送信信号を、送信アンテナ部２１を介して送信するアナログ回路を備える。送受信部３は、受信アンテナ部２２を介して受信される受信信号に基づき、複数の送受信チャネルのそれぞれについて、送信信号と受信信号との周波数差を表すビート信号を生成するアナログ回路を備える。更に、送受信部３は、送受信チャネル毎に生成されたビート信号を、それぞれアナログ－デジタル変換するＡＤ変換器を備え、ＡＤ変換器の出力であるＡＤ変換データを信号処理部４に出力する。
　送受信部３は、具体的には図２に示すように、測定周期Ｔｆを１フレームとして、各フレームの先頭の測定期間Ｔｍの間、時間に対して周波数が直線的に変化するチャープを、予め設定された繰返数Ｌ２だけ連続的に送信する。測定周期Ｔｆは、例えば５０ｍｓとしてもよい。測定期間Ｔｍは、例えば１０ｍｓとしてもよい。繰返数Ｌ２は、例えば２５６回としてもよい。
　信号処理部４は、ＣＰＵ４１と、例えば、ＲＡＭ又はＲＯＭ等の半導体メモリ（以下、メモリ４２）と、を有するマイクロコンピュータを備える。また、信号処理部４は、高速フーリエ変換（以下、ＦＦＴ）処理をハードウェアによって実現するデジタルシグナルプロセッサ（以下、ＤＳＰ）４３を備える。ＤＳＰ４３は、汎用のハードウェア演算器である。但し、ＤＳＰ４３は、必ずしも汎用のものである必要はない。
　信号処理部４におけるＤＳＰ４３以外の各機能は、ＣＰＵ４１が非遷移的実体的記録媒体に格納されたプログラムを実行することにより実現される。この例では、メモリ４２が、プログラムを格納した非遷移的実体的記録媒体に該当する。また、このプログラムが実行されることで、プログラムに対応する方法が実行される。なお、信号処理部４は、１つのマイクロコンピュータを備えてもよいし、複数のマイクロコンピュータを備えてもよい。
　信号処理部４のＤＳＰ４３以外の機能を実現する手法は、ソフトウェアに限るものではなく、その一部又は全部の機能は、一つあるいは複数のハードウェアを用いて実現されてもよい。例えば、上記機能がハードウェアである電子回路によって実現される場合、その電子回路は、デジタル回路、又はアナログ回路、あるいはこれらの組合せによって実現されてもよい。また、ＤＳＰ４３の機能が、ソフトウェアによって実現されてもよい。
　［２．処理］
　次に、信号処理部４のＣＰＵ４１が実行するターゲット検出処理について、図３のフローチャートを用いて説明する。
　本処理は、測定周期Ｔｆの処理サイクル毎に繰り返し実行される。
　本処理が開始されるとＣＰＵ４１は、Ｓ１１０にて、Ｌ１個の送受信チャネルのそれぞれから、ビート信号をＡＤ変換したＡＤ変換データを取得する。
　ＣＰＵ４１は、Ｓ１２０では、Ｌ１個の送受信チャネルのうち、後述するＳ１３０～Ｓ１６０の処理が実行されていない１つの送受信チャネルを、対象チャネルとして選択する。
　ＣＰＵ４１は、Ｓ１３０では、対象チャネルのＡＤ変換データについて、ＤＳＰ４３を使用し、Ｌ２個のチャープのそれぞれについてＦＦＴ処理を実行することにより、Ｌ２個のＦＦＴスペクトラムを生成する。なお、ＦＦＴ処理の際には、処理対象となるデータに予め設定された窓関数が乗じられる。ここで生成されるＦＦＴスペクトラムは、周波数がターゲットまでの距離に対応づけられるため、以下では、距離スペクトラムという。また、距離スペクトラムにおける周波数ビンを距離ビンという。ここでは距離ビンの数をＮｄ個とする。
　ＣＰＵ４１は、Ｓ１４０では、Ｓ１３０にて生成されたＬ２個の距離スペクトラムから集めた同一距離ビンに属するＭ個のデータについて、ＤＳＰ４３を使用してＦＦＴ処理を実行する。このＦＦＴ処理を、距離スペクトラムが有するＮｄ個の距離ビンのそれぞれについて実行することにより、Ｎｄ個のＦＦＴスペクトラムを生成する。なお、ＦＦＴ処理の際には、処理対象となるデータに予め設定された窓関数が乗じられる。ここで生成されるＦＦＴスペクトラムは、周波数が自車とターゲットとの相対速度に対応づけられるため、以下では速度スペクトラムという。また、速度スペクトラムにおける周波数ビンを速度ビンという。ここでは速度ビンの数をＮｖとする。
　図４に示すように、自車に対して一定の相対速度を有した同一ターゲットに基づくビート信号の位相は、時間の経過に伴って相対速度に応じた一定の割合で変化するため、チャープ毎に異なった値が検出される。その結果、同一距離ビンのデータを集めてＦＦＴ処理を実行することにより、相対速度に応じた周波数成分が抽出される。
　以下では、距離ビン毎に生成されたＮｄ個の速度スペクトラムによって表される、距離と相対速度とを２軸とするＮｄ×Ｎｖ個の周波数ビンを有した２次元的なＦＦＴスペクトラムを、以下では、２次元スペクトラムという。また２次元スペクトラムのビンを２次元ビンという。
　ＣＰＵ４１は、Ｓ１５０では、Ｌ１個の送受信チャネルの全てについて、上述のＳ１３０～Ｓ１４０の処理が実行済みであるか否かを判断する。ＣＰＵ４１は、未処理の送受信チャネルが存在する場合にはＳ１２０に戻り、全ての送受信チャネルについて処理が実行済みであればＳ１６０に移行する。
　ＣＰＵ４１は、Ｓ１６０では、送受信チャネル毎に生成されたＬ１個の２次元スペクトラムにおいて、ピークが検出される２次元ビンをピークビンとして抽出する。ここでは、検出されるピークビンの数をＰ個とする。
　ＣＰＵ４１は、Ｓ１７０では、前回までの処理サイクルで検出されたピークビンに基づいて、今回の処理サイクルでピークビンが検出される２次元ビンを推定する。更に、ＣＰＵ４１は、検出されたピークビンと推定されたピークビンとの２次元スペクトラム上での距離が閾値範囲内であれば、同一の物体に基づくピークビンとして対応づけを行う。この処理は、いわゆるトラッキングと呼ばれる技術である。
　ＣＰＵ４１は、Ｓ１８０では、Ｓ１６０で検出されたＰ個のピークビンのそれぞれについて、Ｌ１個の２次元スペクトラムのそれぞれから集めた同一ピークビンに属するＬ１個のデータに基づいて、送受信チャネルの各チャネル間の相関関係を表す相関行列を生成する。この相関行列を中心行列という。なお、相関関係は、各チャネルで得られる信号の類似度合いを表す係数で表現され、係数の絶対値が０に近づくほど相関が弱く、係数の絶対値が１に近づくほど相関が強いことを表す。
　ＣＰＵ４１は、Ｓ１９０では、Ｐ個のピークビンのそれぞれについて、着目するピークビンの周囲に存在する複数の２次元ビンを周辺ビンとして抽出し、抽出された周辺ビンのそれぞれについて、Ｓ１８０で説明した中心行列と同様の手法で、相関行列を生成する。この相関行列を周辺行列という。
　周辺ビンは、例えば、図５の中段に示すように、ピークビンの速度ビン方向に隣接する２つのビン及び距離ビン方向に隣接する２つのビンの合計４つのビンを含んでいてもよい。また、図５の下段に示すように、更に、ピークビンの斜め右上と、斜め右下と、斜め左上と、斜め左下とに位置する４つのビンを加えた合計８つのビンを含んでいてもよい。これに限らず、周辺ビンには、ピークビンを中心として予め設定された範囲内に存在する一つ以上の周波数ビンが含まれていればよい。
　ＣＰＵ４１は、Ｓ２００では、ピークビン毎に、Ｓ１８０で生成された中心行列と、Ｓ１９０で生成された周辺行列とを重み付け加算することで、平均化された一つの相関行列を生成する。この相関行列を統合行列という。これにより、Ｐ個のピークビンのそれぞれに対応したＰ個の統合行列が生成される。この統合行列は、個々の中心行列や周辺行列と比較して、チャネル間の相互相関を表す統合行列の非対角成分の値が抑圧されたものとなる。また、ＣＰＵ４１は、算出された統合行列を、ピークビンを特定する２次元ビンと対応づけて、メモリ４２上の予め設定された記憶領域に記憶させる。
　ＣＰＵ４１は、Ｓ２１０では、Ｐ個の統合行列のそれぞれについて、時間平均法を用いて、更に、チャネル間の相互相関を抑圧する。具体的には、図６に示すように、Ｓ１７０のトラッキング処理で、着目するピークビンに対応づけられる、過去Ｓサイクル分のピークビンを抽出し、抽出された各ピークビンについて生成された統合行列を取得する。Ｓは１以上の整数である。更に、取得した統合行列と、Ｓ１９０で算出された統合行列とを、重み付け加算することによって、時間平均された統合行列を生成する。
　ＣＰＵ４１は、Ｓ２２０では、Ｓ２１０にて生成された時間平均された統合行列を用いて、ピークビン毎に、ＭＵＳＩＣ等の方位推定演算を実行する。この場合、図７に示すように、一つのピークビンから複数の方位が算出される場合がある。従って、ピークビン毎に方位推定演算にて検出された方位の数を、全てのピークビンについて合計した数Ｑ（≧Ｐ）が、ターゲットの数となる。
　ＣＰＵ４１は、Ｓ２３０では、Ｑ個のターゲットのそれぞれについて、ターゲット情報を生成して、本処理を終了する。なお、ターゲット情報には、着目するピークビンに対応する距離ビンから算出された距離と、着目するピークビンに対応する速度ビンから算出された相対速度と、着目するピークビンの相関行列を用いて推定された方位角度とが含まれる。このターゲット情報は、図示しない車載ＬＡＮを介して、車両の各部に配信され、様々な運転支援制御に使用される。
　本処理において、Ｓ１３０とＳ１４０とが解析ステップと解析部とに相当する。Ｓ１６０が周出ステップとＳ抽出部とに相当する。Ｓ１８０が中心生成ステップと中心生成部とに相当する。Ｓ１９０が周辺生成ステップと周辺生成部とに相当する。Ｓ２００が統合ステップと統合部とに相当する。Ｓ２２０が推定ステップと推定部とに相当する。Ｓ２００は、記憶部にも相当する。Ｓ１７０が追跡部に相当する。Ｓ２１０が時間平均部に相当する。
　［３．原理］
　以下では、中心行列と周辺行列との重み付け加算によって、チャネル間の相互相関が抑圧される原理について説明する。ここでは、簡単のため、２次元スペクトラムの代わりに、１次元の距離スペクトラムの場合について説明する。
　等間隔に配置されたＭ個のアンテナで２個のターゲットを検出する場合を仮定する。２個のターゲットは、距離がほぼ等しく、各アンテナから見た方位は、それぞれ少しずつ異なるものとする。ターゲットには多点反射モデルを適用する。一般的に、ターゲットとして大きさのある車両等を想定した場合、１つのターゲットからは複数の反射が返ってくるため、多点反射モデルが有効である。つまり、一つのターゲットについて観測される反射ピークは、Ｎ個の反射点からのＮ点の反射ピークが合成されたものとなる。また、図８に示すように、ｉで識別される同一ターゲット上のＮ個の反射点からの反射波は、Ｍ個のアンテナに同一方位から入射され、隣接するアンテナ間で一定量の位相回転θiを生じさせる。但し、Ｎ個の反射点は微妙に距離が異なるため、Ｎ個の反射点に基づくＮ個の反射ピークは、距離スペクトラム上では、図９に示すように、ターゲット毎に決まるある距離ビンを中心に確率的にばらつく。
　ここで、ピークからδω離れた周波数での窓関数フィルタ応答をＨ（δω）とすると、ｉ番目ターゲットにおける角周波数ω（０）での受信信号ベクトルｘｉ（ω（０））は、（１）式で表される。（１）式は、i番目のターゲットに関するj番目の反射点からの信号の窓関数フィルタ応Ｈ（δωij）を、Ｍ個のアンテナについて順次並べたベクトルを、Ｎ個の反射点について足し合わせることを意味する。
Figure JPOXMLDOC01-appb-M000001
 　これを一般化して、i番目のターゲットにおける角周波数ω（ｋ）での受信信号ベクトルｘｉ（ω（ｋ））は（２）式で表される。
Figure JPOXMLDOC01-appb-M000002
 　ここで、ａjkはω（ｋ）でのω（０）からのフィルタ応答における振幅変化量であり、ａj0＝１となる。φは、窓関数のフィルタ応答による角周波数ステップΔω毎の位相変化量である。但し、窓関数は、直線位相特性を持つため、φは一定量となる。
　次に、（３）式は２つのターゲットからのピーク信号を合成したコンポジット受信信号ｘ（ω（ｋ））を表す。
Figure JPOXMLDOC01-appb-M000003
 　これを用いて、角周波数ω（ｋ）での相関行列Ｒ(k)ｘｘのｐ行ｑ列の成分を求めて整理すると（４）式が得られる。但し、ｐ≠ｑである。
Figure JPOXMLDOC01-appb-M000004
 　これを用いて、ピーク周辺の合計Ｋ個のビンの相関行列を足し合わせると、最終的な相関行列Ｒｘｘは、（５）式で表される。
Figure JPOXMLDOC01-appb-M000005
 　（５）式の演算を行ったときに、（４）式における右辺第３項以下に示すターゲット間の相互相関項は、振幅がプラスにもマイナスにもなり得るため、加算すると値がゼロに近づく。一方、（４）式における右辺第１項と第２項とに示す各ターゲットの自己相関項は、振幅が常にプラスとなり、加算すると一方的に値が増大する。これにより、結果的に相関行列の到来波同士の相互相関が低減される。以上のことから、ターゲットが多点反射モデルで表される場合には、ピークビンの相関行列である中心行列と、周辺ビンの相関行列である周辺行列とを重み付け加算することによって、チャネル間の相互相関が抑圧されることがわかる。
　（６）式は、Ｓ２１０で説明した統合行列を時間平均する処理を表す。Ｗ（ｓ）は、ｓサイクル前に生成された統合行列の重みである。この場合も、周辺ビンの相関行列を用いて平均化する場合と同様に、相互相関項が抑圧される。
Figure JPOXMLDOC01-appb-M000006
 　［４．効果］
　以上詳述した第１実施形態によれば、以下の効果を奏する。
　（４ａ）方位推定装置１では、ピークビンに隣接する周辺ビンを用いて、相関行列のスナップショット数を確保し、これを平均化することで、相関行列の相互相関成分を抑圧する。従って、方位推定装置１によれば、ＦＣＭ方式等のようにチャープ毎のＦＦＴ処理の結果を利用した時間平均法を適用することができない変調が用いられる場合でも、相関行列の相互相関成分を抑圧できる。また、方位推定装置１によれば、空間平均法を用いることなくスナップショット数を確保できるため、不等間隔のアンテナアレーが用いられる場合でも、相関行列の相互相関成分を抑圧できる。その結果、方位推定装置１によれば、ＭＵＳＩＣなど、相関行列を利用した方位推定演算による方位分離の精度と、方位推定の精度とを向上させることができる。
　ここで、図１０は、方位推定演算としてＭＵＳＩＣを使用し、ターゲットのＳ／Ｎを変化させながら、並走する２台の車両を模したターゲットの検出率をシミュレーションによって算出した結果を示す。ここでは、ピークビンのみで相関行列を生成した場合と、４つの周辺ビンを利用して相関行列の相互相関成分を抑圧した場合と、８つの周辺ビンを利用して相関行列の相互相関成分を抑圧した場合とを示す。図１０からは、利用する周辺ピンの数を増やすことで、相互相関成分をより抑圧できることがわかる。但し、利用する周辺ビンの範囲を広げ過ぎると、重なっていない別のピークの影響を受けて劣化する場合があるため、これを考慮して、利用する周辺ビンの範囲を設定してもよい。
　（４ｂ）方位推定装置１では、ピークビンのトラッキングを行うことで、過去Ｓサイクルの間に検出される同一ターゲットに基づくピークビンを抽出し、抽出されたピークビンについて算出された統合行列の時間平均を行う。これより、距離推定演算に使用する相関行列の相互相関成分が抑圧される。つまり、測定周期Ｔｆの間におけるターゲットの方位角度の変化が、方位推定演算における距離分解能に対して十分に小さい場合には、この手法を更に適用することによって、相関行列の相互相関成分が更に抑圧される。その結果、方位推定演算による方位分離の精度と方位推定の精度とをより向上させることができる。
　［５．他の実施形態］
　以上、本開示の実施形態について説明したが、本開示は上述の実施形態に限定されることなく、種々変形して実施することができる。
　（５ａ）上記実施形態では、２次元スペクトラムから抽出されるピークビンについて相関行列を生成する場合について説明したが、本開示は、これに限定されるものではない。例えば、距離スペクトラムから抽出されるピークビンについて相関行列を生成する場合に適用してもよい。
　（５ｂ）上記実施形態では、変調方式としてＦＣＭ方式を用いる場合において、２次元ビンが距離と速度とに対応した２次元スペクトラムであるレンジドップラマップを生成する場合について説明したが、本開示は、これに限定されるものではない。レンジドップラマップを得られるのであれば、どのような変調方式が用いられてもよい。
　（５ｃ）上記実施形態では、相関行列の相互相関成分を抑圧する演算として重み付け加算を用いるが、全ての重みを同一とした単純加算を用いてもよい。
　（５ｄ）上記実施形態では、Ｓ２１０にて統合行列を時間平均する処理を実行するが、この処理は、省略されてもよい。
　（５ｅ）上記実施形態における１つの構成要素が有する複数の機能を、複数の構成要素によって実現したり、１つの構成要素が有する１つの機能を、複数の構成要素によって実現したりしてもよい。また、複数の構成要素が有する複数の機能を、１つの構成要素によって実現したり、複数の構成要素によって実現される１つの機能を、１つの構成要素によって実現したりしてもよい。また、上記実施形態の構成の一部を省略してもよい。また、上記実施形態の構成の少なくとも一部を、他の上記実施形態の構成に対して付加又は置換してもよい。
　（５ｆ）上述した方位推定方法及び方位推定装置の他、種々の形態で本開示を実現することもできる。例えば、当該方位推定装置を構成要素とするシステム、当該方位推定装置としてコンピュータを機能させる又は方位推定方法を実現するためのプログラム、このプログラムを記録した半導体メモリ等の非遷移的実態的記録媒体などとして実現されてもよい。

Claims (6)

  1. 　複数の送受信チャネルから得られる信号を周波数解析することで、前記送受信チャネル毎に周波数スペクトラムを生成する解析ステップ（Ｓ１３０、Ｓ１４０）と、
    　前記解析ステップにて生成された周波数スペクトラムからピークとなる周波数ビンであるピークビンを抽出する抽出ステップ（Ｓ１６０）と、
    　前記抽出ステップにて抽出されたピークビン毎に、全ての送受信チャネルから集めた同一の前記ピークビンの値を用いて、相関行列である中心行列を生成する中心生成ステップ（Ｓ１８０）と、
    　前記ピークビンを中心として予め設定された範囲内に存在する一つ以上の周波数ビンを周辺ビンとして、前記ピークビンのそれぞれについて前記周辺ビン毎に、全ての送受信チャネルから集めた同一の前記周辺ビンの値を用いて、相関行列である周辺行列を生成する周辺生成ステップ（Ｓ１９０）と、
    　前記ピークビン毎に、前記中心生成ステップにて生成された中心行列と、前記周辺生成ステップにて生成された１つ以上の周辺行列とを重み付け加算することで、相関行列である統合行列を生成する統合ステップ（Ｓ２００）と、
    　前記統合ステップにて生成された統合行列を用いて方位推定演算を実行する推定ステップ（Ｓ２２０）と、
    　を備える方位推定方法。
  2. 　複数の送受信チャネルから得られる信号を周波数解析することで、前記送受信チャネル毎に周波数スペクトラムを生成するように構成された解析部（Ｓ１３０、Ｓ１４０）と、
    　前記解析部にて生成された周波数スペクトラムからピークとなる周波数ビンであるピークビンを抽出するように構成された抽出部（Ｓ１６０）と、
    　前記抽出部にて抽出されたピークビン毎に、全ての送受信チャネルから集めた同一の前記ピークビンの値を用いて、相関行列である中心行列を生成するように構成された中心生成部（Ｓ１８０）と、
    　前記ピークビンを中心として予め設定された範囲内に存在する一つ以上の周波数ビンを周辺ビンとして、前記ピークビンのそれぞれについて前記周辺ビン毎に、全ての送受信チャネルから集めた同一の前記周辺ビンの値を用いて、相関行列である周辺行列を生成するように構成された周辺生成部（Ｓ１９０）と、
    　前記中心生成部が生成する中心行列と、前記周辺生成部が生成する１つ以上の周辺行列とを重み付け加算することで、相関行列である統合行列を生成するように構成された統合部（Ｓ２００）と、
    　前記統合部で生成された統合行列を用いた方位推定演算を実行するように構成された推定部（Ｓ２２０）と、
    　を備える方位推定装置。
  3. 　請求項２に記載の方位推定装置であって、
    　前記統合部にて算出された統合行列を、複数の処理サイクルに渡って、予め設定された記憶領域に記憶させる記憶部（Ｓ２００）と、
    　前記複数の処理サイクルに渡って、同一の物体に基づくと推定される前記ピークビンを互いに対応づける追跡部（Ｓ１７０）と、
    　前記追跡部により対応づけられた前記ピークビンについて前記記憶部に記憶された統合行列を、重み付け加算する時間平均部（Ｓ２１０）と、
    　を更に備え、
    　前記推定部は、前記時間平均部により重み付け加算された統合行列を用いて、前記方位推定演算を実行するように構成された、
    　方位推定装置。
  4. 　請求項２又は請求項３に記載の方位推定装置であって、
    　複数のアンテナを有するアンテナ部（２）と、
    　前記アンテナ部が有する前記複数のアンテナのうち、送信に用いられるアンテナを送信アンテナ、受信に用いられるアンテナを受信アンテナとして、前記送信アンテナと前記受信アンテナとの組み合わせである複数の送受信チャネルを介して、予め設定された変調方式で変調された信号を送受信し、送信信号と受信信号との周波数差を表すビート信号を生成する送受信部（３）と、
    　を更に備え、
    　前記解析部は、前記送受信部にて生成されたビート信号を周波数解析するように構成された
    　方位推定装置。
  5. 　請求項４に記載の方位推定装置であって、
    　前記送受信部は、ＦＣＭ方式で変調し、
    　前記解析部は、チャープ毎に周波数解析をした結果である距離スペクトラムと、前記距離スペクトラムを、該距離スペクトラムの周波数ビン毎に、複数のチャープに渡って周波数解析した結果である速度スペクトラムとにより表される２次元スペクトラムを生成する
    　方位推定装置。
  6. 　請求項４又は請求項５に記載の方位推定装置であって、
    　前記アンテナ部に属する前記複数のアンテナが不等間隔に配置された
    　方位推定装置。
PCT/JP2019/005315 2018-02-20 2019-02-14 方位推定方法及び方位推定装置 WO2019163632A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US16/996,622 US11125871B2 (en) 2018-02-20 2020-08-18 Azimuth estimation device and method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2018-027907 2018-02-20
JP2018027907A JP7056212B2 (ja) 2018-02-20 2018-02-20 方位推定方法および装置

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US16/996,622 Continuation US11125871B2 (en) 2018-02-20 2020-08-18 Azimuth estimation device and method

Publications (1)

Publication Number Publication Date
WO2019163632A1 true WO2019163632A1 (ja) 2019-08-29

Family

ID=67688427

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2019/005315 WO2019163632A1 (ja) 2018-02-20 2019-02-14 方位推定方法及び方位推定装置

Country Status (3)

Country Link
US (1) US11125871B2 (ja)
JP (1) JP7056212B2 (ja)
WO (1) WO2019163632A1 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3839561A1 (en) * 2019-12-18 2021-06-23 Imec VZW Radar ranging
WO2023021586A1 (ja) * 2021-08-18 2023-02-23 三菱電機株式会社 Ｍｉｍｏレーダ信号処理装置及びその受信信号処理装置、並びに着目受信信号ベクトルの伝搬モード判別方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006145251A (ja) * 2004-11-16 2006-06-08 Denso Corp 電波到来方向推定装置
US20070080855A1 (en) * 2003-09-03 2007-04-12 Gerlach Karl R Radar processor system and method
JP2010271115A (ja) * 2009-05-20 2010-12-02 Toshiba Corp レーダ装置
WO2014199609A1 (ja) * 2013-06-13 2014-12-18 パナソニック株式会社 レーダ装置
JP2017090229A (ja) * 2015-11-10 2017-05-25 富士通テン株式会社 到来方向推定装置、到来方向推定方法、到来方向推定プログラム

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003240842A (ja) * 2002-02-14 2003-08-27 Murata Mfg Co Ltd レーダ
US7714782B2 (en) * 2004-01-13 2010-05-11 Dennis Willard Davis Phase arrays exploiting geometry phase and methods of creating such arrays
JP4780229B2 (ja) 2004-07-06 2011-09-28 株式会社デンソー レーダ装置
JP4867200B2 (ja) 2004-07-06 2012-02-01 株式会社デンソー レーダ装置
JP4498269B2 (ja) * 2005-11-30 2010-07-07 株式会社デンソーアイティーラボラトリ レーダ信号処理装置
JP4946101B2 (ja) * 2006-03-08 2012-06-06 株式会社デンソー 方位検出装置
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
JP2009180538A (ja) * 2008-01-29 2009-08-13 Mitsubishi Electric Corp レーダ装置
DE102009024918B4 (de) * 2008-06-17 2022-04-28 Denso Corporation Zielobjekterkennungsvorrichtung zur Erlangung von Informationen betreffend Zielobjekte basierend auf einer Korrelationsmatrix, die aus Signalwerten erhalten wird, welche reflektierten elektromagnetischen Wellen entsprechen
DE102008038365A1 (de) * 2008-07-02 2010-01-07 Adc Automotive Distance Control Systems Gmbh Fahrzeug-Radarsystem und Verfahren zur Bestimmung einer Position zumindest eines Objekts relativ zu einem Fahrzeug
JP6092596B2 (ja) * 2012-11-28 2017-03-08 富士通テン株式会社 レーダ装置、および、信号処理方法
CN114185042A (zh) * 2015-09-17 2022-03-15 松下电器产业株式会社 雷达装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070080855A1 (en) * 2003-09-03 2007-04-12 Gerlach Karl R Radar processor system and method
JP2006145251A (ja) * 2004-11-16 2006-06-08 Denso Corp 電波到来方向推定装置
JP2010271115A (ja) * 2009-05-20 2010-12-02 Toshiba Corp レーダ装置
WO2014199609A1 (ja) * 2013-06-13 2014-12-18 パナソニック株式会社 レーダ装置
JP2017090229A (ja) * 2015-11-10 2017-05-25 富士通テン株式会社 到来方向推定装置、到来方向推定方法、到来方向推定プログラム

Also Published As

Publication number Publication date
JP7056212B2 (ja) 2022-04-19
US20200379100A1 (en) 2020-12-03
JP2019144077A (ja) 2019-08-29
US11125871B2 (en) 2021-09-21

Similar Documents

Publication Publication Date Title
US11874395B2 (en) Radar processing chain for frequency-modulated continuous wave radar systems
EP3471210B1 (en) Radar apparatus
JP6677408B2 (ja) 到来方向推定装置、到来方向推定方法、到来方向推定プログラム
US7847733B2 (en) Direction-of-arrival estimating device and program
CN111736131B (zh) 一种剔除一比特信号谐波虚假目标的方法及相关组件
JP2009162688A (ja) 電子走査型レーダ装置、受信波方向推定方法及び受信波方向推定プログラム
Northardt et al. Spatial compressive sensing for direction-of-arrival estimation with bias mitigation via expected likelihood
JP2012132846A (ja) 電子走査型レーダ装置、受信波方向推定方法及び受信波方向推定プログラム
JP2006258529A (ja) 電波到来方向推定装置及び方法
JP6556399B2 (ja) レーダ装置
CN106486769B (zh) 用于线性相控阵天线的空间插值方法和设备
JP5114217B2 (ja) 電子走査型レーダ装置、受信波方向推定方法及び受信波方向推定プログラム
JP2013152239A (ja) 電子走査型レーダ装置、その制御方法及びプログラム
JP4757629B2 (ja) 到来方位推定装置
WO2019163632A1 (ja) 方位推定方法及び方位推定装置
JP2014228291A (ja) 無線検出装置及び無線検出方法
JP5553980B2 (ja) 電波方向探知装置およびビーム形成装置
US11754671B2 (en) Incoming wave count estimation apparatus and incoming wave count incoming direction estimation apparatus
KR102099388B1 (ko) 안테나 어레이 외삽을 이용한 레이더 수신신호의 도착방향 추정 방법 및 장치
Wu et al. Switched-element direction finding
CN114265058A (zh) Mimo雷达目标测角方法、装置、电子设备及存储介质
JP7191262B1 (ja) レーダ装置
CN109752688B (zh) 一种针对传感器阵列系统的临近信源角度差值计算方法
Coventry Fine-time resolution passive RF source imaging
Chang et al. DOA and DOD estimation using orthogonal projection approach for bistatic MIMO radars

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 19757248

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 19757248

Country of ref document: EP

Kind code of ref document: A1