WO2021240927A1 - レーダ装置 - Google Patents

レーダ装置 Download PDF

Info

Publication number
WO2021240927A1
WO2021240927A1 PCT/JP2021/007215 JP2021007215W WO2021240927A1 WO 2021240927 A1 WO2021240927 A1 WO 2021240927A1 JP 2021007215 W JP2021007215 W JP 2021007215W WO 2021240927 A1 WO2021240927 A1 WO 2021240927A1
Authority
WO
WIPO (PCT)
Prior art keywords
received signal
radar device
target
processing unit
virtual
Prior art date
Application number
PCT/JP2021/007215
Other languages
English (en)
French (fr)
Inventor
幸徳 赤峰
美裕 中尾
堅一 嶋田
浩司 黒田
Original Assignee
日立Astemo株式会社
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 日立Astemo株式会社 filed Critical 日立Astemo株式会社
Priority to DE112021001201.0T priority Critical patent/DE112021001201T5/de
Publication of WO2021240927A1 publication Critical patent/WO2021240927A1/ja

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/93Radar or analogous systems specially adapted for specific applications for anti-collision purposes
    • G01S13/931Radar or analogous systems specially adapted for specific applications for anti-collision purposes of land vehicles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06Systems determining position data of a target
    • G01S13/42Simultaneous measurement of distance and other co-ordinates

Definitions

  • the present invention is a technique relating to a radar device.
  • millimeter-wave radar For autonomous driving of vehicles, high-resolution millimeter-wave radar is required in scenes such as free space detection at intersections and parking lots and vehicle detection at the end of traffic jams on highways.
  • the MUSIC method (Multiple Signal Classification) is known as a method for achieving high resolution with a limited aperture length.
  • the MUSIC method is a method that utilizes the orthogonality between the eigenvectors of the signal and the eigenvectors of noise. After calculating the correlation matrix of the array received signal, the number of wave sources is estimated by expanding the eigenvalues, and null detection caused by the orthogonality is performed while sweeping the angle. Perform the calculation.
  • Patent Document 1 is a method of virtually increasing the number of received antennas by the Katri-Lao product, and after generating a correlation matrix by array received signals, an extended signal is generated by non-overlapping elements of the matrix.
  • the MUSIC method has the disadvantage that the calculation load is heavy while high resolution can be obtained because null detection is performed by angle sweep.
  • DSP Digital Signal Processor
  • Patent Document 1 after the correlation matrix based on the array received signal is generated, the extended signal is generated by the non-overlapping elements of the matrix, but an error occurs in the direction detection of the target due to the correlation between the targets. In order to suppress this error, Patent Document 1 avoids the problem by repeatedly processing the generation of the extended signal and the direction estimation while narrowing the search range. However, even for one target, it is necessary to repeatedly execute the generation of the extended signal and the direction estimation as the loop processing, and there is a concern that the calculation load becomes very heavy depending on the convergence of the loop processing.
  • An object of the present invention is to provide a radar device that reduces a computational load and enables high-resolution processing.
  • a receiving antenna for receiving a received signal from a target, a processing unit for generating a virtual received signal when received by a virtual antenna based on the received signal, and the above-mentioned It has a target direction calculation unit that calculates the direction of the target based on the received signal and the virtual received signal, and the processing unit estimates the direction of the target from the received signal. It is a radar device that performs estimation processing and generates the virtual received signal based on the estimated orientation.
  • FIG. It is a block block diagram of the radar apparatus in Example 1.
  • FIG. It is a figure which imaged the received signal from a certain angle direction by a receiving array antenna and a virtual antenna. It is an image diagram when the amplitude / phase information in the received signal of the receiving array antenna is plotted by dividing into the real number / imaginary part, and the polynomial approximation is made. It is an angular directional spectrum obtained by processing the received signal of the receiving array antenna in the spatial direction. It is an image diagram when the amplitude / phase information in the received signal of the receiving array antenna is plotted by dividing into the real number / imaginary part, and the polynomial approximation is made.
  • FIG. 1 is a block configuration diagram of a radar device according to the first embodiment.
  • the radar device is composed of an array antenna / analog part and a digital part.
  • the array antenna / analog section outputs a transmission signal from the transmission antenna (101).
  • a synthesizer (103) is used as the transmission signal, and a chirp signal whose frequency is linearly changed with time is often used.
  • the transmitted transmission signal is reflected by the target target, and a part of the reflected wave is returned to the reception array antenna (102) and received.
  • the receiving array antenna (102) is composed of a plurality of antenna elements.
  • the signal received by the reception array antenna (102) is down-converted by the mixer (104), which is a frequency converter, to convert the millimeter wave reception signal into a low frequency signal.
  • the signal of the synthesizer (103) is used as the local signal of the mixer (104).
  • the time difference between the transmission signal and the reception signal that is, the frequency corresponding to the distance of the target target is output from the mixer (104).
  • the output signal of the mixer (104) is converted into a digital signal by the A / D converter (105) and transmitted to the digital unit.
  • the time / frequency FFT processing unit (106) processes the time / frequency FFT with respect to the received signal received by each antenna element constituting the reception array antenna (102).
  • FFT represents a Fast Fourier Transform.
  • the processing of each block constituting the digital unit may be executed by a CPU such as a microcomputer by reading a program recorded in a recording device and executing it as software, or by using dedicated hardware such as a circuit. You may do so.
  • the distance distribution can be grasped by performing the frequency FFT processing, and the velocity distribution can be grasped by performing the time FFT processing. be able to.
  • FIG. 2 is an image of a received signal from a certain angle direction by the receiving array antenna (102) and the virtual antenna in the first embodiment.
  • the time / frequency FFT processing unit (106) has complex information (201) x1, x2, x3, x4, regarding the received signal from the angular orientation (203) of the target. Generate x5.
  • the electromagnetic wave of the received signal is the wavelength (204).
  • a phase difference occurs in the signal received by each antenna element when the angular direction of the target is oblique from the front.
  • a virtual received signal is generated by the virtual antenna (205).
  • the virtual antenna (205) is arranged on the right side of the receiving array antenna (102) made of a real antenna element, but it may be arranged on the left side.
  • the complex information (201) is determined by the arrival direction of the electromagnetic wave, that is, the angular direction of the target (203), the wavelength of the electromagnetic wave (204), and the arrangement of the receiving array antenna (102).
  • the complex information (201) changes regularly among the receiving antenna elements constituting the receiving array.
  • the angular direction (203) of the target can be obtained by the regular change of the complex information (201).
  • the aperture length (202) is the distance from one end of the receiving array antenna (102) to the other.
  • In-vehicle radars are often installed inside vehicle bumpers or headlights, and the installation space is limited, making it difficult to increase the opening length (202).
  • the MUSIC method is well known as an algorithm for obtaining high resolution with a limited aperture length, but it is a problem in an automatic driving system that has a relatively heavy calculation load and needs to process multiple targets in a short time. Become.
  • the complex information of the virtual antenna (205) virtually extended outside the aperture length (202) is used by using regular changes in the complex information (201). It would be good if it could be estimated by a simple calculation as much as possible. Examples of simple calculations include a method of performing polynomial approximation on complex information (201) and estimating a virtual antenna (205), and a method using a learner.
  • the polynomial approximation processing unit (301) of the real number part and the polynomial approximation processing unit (302) of the imaginary number part separate the complex information (201) into the real number part and the imaginary number part, and then horizontally as shown in FIG. Plot is performed with the axis as the antenna position (303) and the vertical axis as the real value (304) or the imaginary value (305).
  • FIG. 3 is an image diagram when the amplitude / phase information in the received signal of the receiving array antenna is plotted by dividing it into real / imaginary parts and polynomial approximation is performed.
  • the vertical axis shows a real value or an imaginary value
  • the horizontal axis shows the antenna position.
  • complex information (201) at the antenna positions of the five antenna elements is plotted (306) in a circle.
  • the complex information of the virtual antenna position is estimated by drawing a polynomial approximation curve (307) to the outside of the aperture length (202) with respect to the plotted complex information.
  • the polynomial approximation is an example of a method of estimating the extended range from the plotted complex information (201), and as another estimation method, the extended range is calculated by using a learner trained by machine learning using data having a known orientation. It is also possible to estimate.
  • FIG. 4 is an angular directional spectrum obtained by processing the received signal of the receiving array antenna in the spatial direction. Further, FIG. 4 shows the complex information (201) using the target angle directional spectrum (401) calculated by spatial FFT processing and the virtual complex information extended by the polynomial approximation curve (307). It is a figure which shows the target angle direction spectrum (402) which processed the space FFT. The vertical axis shows the power level (dB) for the received electromagnetic wave, and the horizontal axis shows the angular direction (deg.) With the front as 0 degree.
  • FIG. 4 shows the simulation results when the target exists at the positions of +0.8 degrees and -0.8 degrees from the front, and the resolution is obtained by expanding the aperture length by 2.5 times with the polynomial approximation curve. Is improved, and it can be seen that the two targets can be detected separately.
  • the target detection near the front which is +/- 0.8 degrees, is taken as an example, but the deviation of the complex information (201) between adjacent antennas becomes large at an angle away from the front, and the polynomial. It becomes difficult to approximate and estimate the extended range by the learner.
  • FIG. 5 is an image diagram when the amplitude / phase information in the received signal of the receiving array antenna is plotted by dividing it into real / imaginary parts and polynomial approximation is performed.
  • the vertical axis shows a real value or an imaginary value
  • the horizontal axis shows the antenna position.
  • FIG. 5 is a diagram showing an example in the case where the target is present in the direction away from the front. Unlike the case of FIG. 3, since the deviation of the adjacent complex information (201) is large, the change between the plots (306) is large, and the polynomial approximation curve (307) does not fit the plot (306) well. Recognize.
  • the complex information (201) of the real antenna is subjected to spatial FFT processing by the first spatial FFT processing unit (107) in FIG. 1 before expansion by the approximate curve or the learner. ..
  • the first angle / direction estimation process the rough angle / direction of the target can be estimated in advance from the received signal from the target.
  • the spatial FFT process is executed as the first angle direction estimation process because the angle direction can be estimated by a lightweight process.
  • the angle direction may be estimated by digital beamforming instead of the first spatial FFT processing unit (107).
  • the first phase conversion processing unit (108) phase-converts the complex information (201) into the estimated direction.
  • This phase conversion is a simple phase rotation process, and the estimated azimuth becomes pseudo front in this conversion.
  • the deviation of the adjacent plot (306) shown in FIG. 5 is suppressed.
  • FIG. 6 is an angular directional spectrum obtained by processing the received signal of the receiving array antenna in the spatial direction. Further, FIG. 6 is a target obtained by subjecting the complex information obtained by combining the plots (306) of the real part and the imaginary part in FIG. 5 to the spatial FFT processing unit (107) in the first spatial FFT processing unit (107). It is an angular azimuth spectrum (601). The vertical axis of FIG. 6 indicates the power level (dB), and the horizontal axis indicates the angular direction (deg.) With the front as 0 degree.
  • the peak (602) of the target angle directional spectrum is about 3 degrees in this case, and the first phase conversion processing unit (108) in FIG. 1 performs a phase rotation process for converting the angle directional 3 degrees to 0 degrees. conduct.
  • the polynomial approximation processing unit (301) of the real number part and the polynomial approximation processing unit (302) of the imaginary number part The phase-transformed real part polynomial approximation process and the phase-converted imaginary number polynomial approximation process are performed, respectively.
  • the extrapolation / interpolation processing unit (109) and the extrapolation / interpolation processing unit (112) form the complex information of the real part at the virtual antenna position and the virtual antenna position based on the polynomial approximation curve (307), respectively. Calculate the complex information of the imaginary part in.
  • FIG. 7 is an image diagram when the amplitude / phase information in the received signal of the receiving array antenna is plotted by dividing it into real / imaginary parts and polynomial approximation is performed.
  • the vertical axis shows a real value or an imaginary value
  • the horizontal axis shows the antenna position.
  • FIG. 7 is a result of polynomial approximation by dividing the complex information (201) into a real part and an imaginary part after phase-converting the complex information (201) from the estimation result of the angle orientation of 3 degrees in FIG.
  • the plot (306) having a large change in FIG. 5 is suppressed, and the polynomial approximation curve (307) in FIG. 7 works well with the plot (306). You can see that they match.
  • the deviation of the complex information plotted (307) and plotted (306) is calculated as the accuracy in the virtual antenna.
  • the polynomial approximation curve (307) is extended to a range of 2.5 times the opening length (202) as in FIG.
  • extrapolation processing is performed to extend to the outside of the antenna arrangement in order to obtain high resolution, but depending on the antenna arrangement, complex information between the antennas is estimated as a measure to suppress the virtual image of the angular direction. It is also possible to easily perform the interpolation process.
  • This extrapolation and interpolation processing is executed by the extrapolation / interpolation processing unit (109) in FIG.
  • the complex information calculated by the extrapolation / interpolation processing is phase-converted to the original direction.
  • the second spatial FFT processing unit (111) performs spatial FFT processing based on the received signal and the virtual received signal, and executes the second angular direction estimation processing.
  • the second spatial FFT processing unit (111) outputs position information such as the angular direction of the target, the distance to the target, and the speed of the target as radar output.
  • Digital beamforming may be used instead of spatial FFT processing.
  • a second angular directional estimation process such as spatial FFT processing or digital beamforming is used to calculate the angular directional of a target with high resolution and accuracy.
  • FIG. 8 is an angular directional spectrum obtained by processing the received signal of the receiving array antenna in the spatial direction.
  • the vertical axis of FIG. 8 shows the power level (dB), and the horizontal axis shows the angular direction (deg.) With the front as 0 degree.
  • FIG. 8 is an example of the angle-direction spectrum after the spatial FFT processing when the target is present at the positions of the angle directions of 2.2 degrees and 3.8 degrees.
  • the target angle-direction spectrum (601) having one peak at the position of 3 degrees is the result of spatial FFT processing of the complex information of the real antenna, that is, the plot (306) of the real part and the imaginary part in FIG. ..
  • the target angle azimuth spectrum (802) having one peak near the angle azimuth 0 degrees is the result of spatial FFT processing of the virtual complex information expanded by the polynomial approximation curve (307) in FIG. It can be seen that it differs from the target position of.
  • the calculation load of the high resolution processing is reduced, the noise immunity is enhanced, and the effect of extending the detection distance can be obtained.
  • Example 2 will be described. The same contents as in the first embodiment will be omitted.
  • the window function processing unit (1201) may execute the window function processing in the processing before the second spatial FFT processing unit (111). That is, in the second phase conversion processing unit (110), the window function processing may be executed for the received signal and the virtual received signal based on the complex information calculated by the extrapolation / interpolation processing.
  • the window function processing unit (1201) performs the window function processing, the angular azimuth spectrum is blunted and the resolution is deteriorated.
  • the area of the virtual antenna (205) shown in FIG. 2 may be expanded, which is synonymous with expanding the expansion range by the polynomial approximation curve of FIGS. 3 and 7.
  • the range far from the reference plot (306) is estimated, there is a possibility that the estimation deviates greatly from the correct answer.
  • FIG. 9 is an image of preventing the divergence of the polynomial approximation by arranging zeros in a pseudo manner.
  • the vertical axis shows the real value
  • the horizontal axis shows the antenna position.
  • a pseudo plot (901) with a zero value is arranged outside the plot of complex information (306) as shown in FIG.
  • the credibility of the polynomial approximation curve (903) arises as a concern.
  • the complex information about the received signal is emphasized and the pseudo plot (901) is disregarded when the polynomial approximation is performed.
  • a polynomial approximation with varying weights may be performed between the plot of information (306) and the pseudo plot (901). In this case, the process is the same as obtaining the polynomial approximation curve (903) after overwriting the plot (306) of the complex information.
  • FIG. 10 is an image diagram when the amplitude / phase information in the received signal of the receiving array antenna is plotted by dividing it into real / imaginary parts and polynomial approximation is performed.
  • the vertical axis shows a real value or an imaginary value
  • the horizontal axis shows the antenna position. More specifically, FIG. 10 shows a polynomial approximation curve (903) in which the zero value is pseudo-plotted outside the complex information plot (306) and the weighted plot (306) is emphasized.
  • the weighting is 100: 1, and the extension range of the polynomial approximation curve is about 20 times the opening length.
  • the polynomial approximation curve (903) does not diverge even in a range far from the aperture length, and the fitting to the plot (306) is also good, and it can be inferred that the answer is close to the correct answer.
  • FIG. 11 is an angular directional spectrum obtained by processing the received signal of the receiving array antenna in the spatial direction.
  • the vertical axis of FIG. 1 indicates the power level (dB), and the horizontal axis indicates the angular direction (deg.) With the front as 0 degree.
  • the target angle azimuth spectrum (401) calculated by spatial FFT processing of the complex information (201) as in FIG. 4 and the polynomial approximation curve (903) described in FIG. 10 are used to open the aperture.
  • the target angle azimuth spectrum (1101) subjected to spatial FFT processing using virtual complex information expanded to about 20 times the length is shown.
  • FIG. 11 shows the simulation results when the target is present at the positions of +0.8 degrees and -0.8 degrees from the front, respectively, from the suppression of the side lobes and the steep spectrum of the target angle directional spectrum (1101). It can be seen that the characteristics that are resistant to noise and have very high resolution are obtained.
  • FIG. 12 is a block configuration diagram of the radar device according to the second embodiment.
  • FIG. 12 shows the overall configuration of the processing content described with reference to FIGS. 9, 10 and 11.
  • the functional blocks added to FIG. 1 described in the first embodiment include the zero padding and weighted processing unit (1202) and the window function processing unit (1201) described in FIG.
  • the first phase conversion processing unit (108) to the second space for each detected estimated angular direction.
  • the process by the FFT processing unit (111) may be processed by the iterative loop processing unit (1203) so as to perform the iterative loop processing according to the number of detections.
  • Example 2 it is possible to obtain characteristics that are resistant to noise and have very high resolution due to the suppression of the side lobes of the angular azimuth spectrum and the steep spectrum.

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)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本発明の目的は、計算負荷を低減し、高分解能処理を可能とするレーダ装置を提供することにある。 本発明のレーダ装置は、物標からの受信信号を受信する受信アンテナ(102)と、前記受信信号に基づいて、仮想的なアンテナ(205)により受信した場合の仮想的な受信信号を生成する処理部と、前記受信信号および前記仮想的な受信信号に基づいて、前記物標の方位を算出する物標方位算出部と、を有する。前記処理部は、前記受信信号から前記物標の方位を推定する角度方位推定処理を行い、推定した前記方位に基づいた前記仮想的な受信信号を生成する。

Description

レーダ装置
 本発明は、レーダ装置に関する技術である。
 車両の自動運転向けに、交差点や駐車場でのフリースペース検知や、高速道路の渋滞末尾の車両検知といったシーンで、高分解能なミリ波レーダが要求される。
 街中の交差点では、交差点の周囲を横断する歩行者や交差点を横切る車両、渋滞により交差点の向こう側で停車している車両など多くの物標を検知する必要がある。例として、60m幅の交差点の向こう側に停車している車両の車線を検知するためには、角度方位として約3度の分解能が必要である。また、高速道路の渋滞末尾の車両検知では、車両の制動距離により150m以上の遠方から渋滞末尾を検知する必要がある。このとき停車している車両の車線を検知するためには1度以下の分解能が必要となる。
 高分解能を実現するためには、受信アンテナの端から端までの距離である開口長を長くすることで実現可能であるが、開口長を長くするためにはアンテナ搭載基板のサイズを長くする必要がありコストが高くなることと、サイズにより設置できる場所が限定される欠点が生じる。
 限定された開口長で、高分解能を実現する方法としてはMUSIC法(Multiple Signal Classification)が知られている。MUSIC法は信号の固有ベクトルと雑音の固有ベクトルの直交性を利用した方法で、アレー受信信号の相関行列を算出した後、固有値展開により波源数推定を行い、角度をスイープさせながら直交性により生じるヌル検出演算を行う。
 特許文献1は、カトリ・ラオ積により仮想的に受信アンテナ数を増やす方式であり、アレー受信信号による相関行列を生成後、行列の非重複要素により拡張信号を生成する。
特開2018-48978
 MUSIC法は、角度スイープによるヌル検出のため高い分解能が得られる一方で計算負荷が重いことが欠点としてあげられる。汎用的なDSP(Digital Signal Processor)を用いた場合には、50ミリ秒のフレームレートに対して20ターゲット程度の検出に留まることが見込まれるため、街中の交差点といったターゲットが多いシーンでは、高分解能処理が追い付かない可能性がある。
 特許文献1では、アレー受信信号による相関行列を生成後、行列の非重複要素により拡張信号を生成するが、物標間の相関によって物標の方位検出に誤差が生じる。この誤差を抑圧するために、特許文献1では、拡張信号の生成と方位推定を、探索範囲を狭めながら繰り返し処理を行うことで課題を回避している。しかし、1つの物標に対してでも、拡張信号の生成と方位推定をループ処理として繰り返し実行する必要があり、ループ処理の収束によっては計算負荷が非常に重くなる懸念がある。
 本発明の目的は、計算負荷を低減し、高分解能処理を可能とするレーダ装置を提供することにある。
 本発明の好ましい一例としては、物標からの受信信号を受信する受信アンテナと、前記受信信号に基づいて、仮想的なアンテナにより受信した場合の仮想的な受信信号を生成する処理部と、前記受信信号および前記仮想的な受信信号に基づいて、前記物標の方位を算出する物標方位算出部と、を有し、前記処理部は、前記受信信号から物標の方位を推定する角度方位推定処理を行い、推定した前記方位に基づいた前記仮想的な受信信号を生成するレーダ装置である。
 本発明によれば、計算負荷を低減し、高分解能処理を可能とするレーダ装置を実現できる。
実施例1におけるレーダ装置のブロック構成図である。 受信アレーアンテナおよび仮想アンテナによる、ある角度方位からの受信信号をイメージした図である。 受信アレーアンテナの受信信号における振幅/位相情報を実数/虚数部にわけてプロットし、多項式近似したときのイメージ図である。 受信アレーアンテナの受信信号を空間方向に処理することで得られる角度方位スペクトルである。 受信アレーアンテナの受信信号における振幅/位相情報を実数/虚数部にわけてプロットし、多項式近似したときのイメージ図である。 受信アレーアンテナの受信信号を空間方向に処理することで得られる角度方位スペクトルである。 受信アレーアンテナの受信信号における振幅/位相情報を実数/虚数部にわけてプロットし、多項式近似したときのイメージ図である。 受信アレーアンテナの受信信号を空間方向に処理することで得られる角度方位スペクトルである。 疑似的にゼロを配置することで多項式近似の発散を防ぐことをイメージした図である。 受信アレーアンテナの受信信号における振幅/位相情報を実数/虚数部にわけてプロットし、多項式近似したときのイメージ図である。 受信アレーアンテナの受信信号を空間方向に処理することで得られる角度方位スペクトルである。 実施例2におけるレーダ装置のブロック構成図である。
 以下、図面を用いて実施例を説明する。
 図1は、実施例1におけるレーダ装置のブロック構成図である。レーダ装置は、アレーアンテナ/アナログ部とデジタル部で構成される。
 アレーアンテナ/アナログ部では、送信アンテナ(101)から送信信号を出力する。ここで、送信信号はシンセサイザ(103)を使い、周波数を時間で線形遷移させたチャープ信号がよく用いられる。送信された送信信号は、ターゲットである物標で反射し、その反射波の一部が受信アレーアンテナ(102)に戻って受信される。受信アレーアンテナ(102)は複数のアンテナ素子から構成される。
 受信アレーアンテナ(102)で受信した信号は、周波数変換器であるミキサ(104)により、ミリ波の受信信号を、低い周波数の信号に変換するダウンコンバートがなされる。このとき、ミキサ(104)のローカル信号として、シンセサイザ(103)の信号を用いる。このことで、送信信号と受信信号の時間差、すなわちターゲットである物標の距離に応じた周波数がミキサ(104)から出力される。ミキサ(104)の出力信号はA/Dコンバータ(105)でデジタル信号に変換されデジタル部に伝送される。
 デジタル部では、受信アレーアンテナ(102)を構成する各アンテナ素子で受信した受信信号に対して、時間/周波数FFT処理部(106)が時間/周波数FFTの処理を行う。ここで、FFTは高速フーリエ変換(Fast Fourier Transform)を表す。デジタル部を構成する各ブロックの処理は、マイコンなどのCPUが、記録装置に記録したプログラムを読み出してソフトウエアとして実行されるようにしてもよいし、回路などの専用のハードウエアで実行されるようにしてもよい。
 周波数FFT処理を行うことで距離分布、時間FFT処理を行うことで速度分布を把握することができ、時間/周波数の2次元FFT処理を行えば、距離と速度の2軸での分布を把握することができる。
 ターゲットとする距離と速度を決定すれば、所定の距離/速度での受信アレーアンテナ(102)を構成する各アンテナ素子の受信信号についての振幅と位相の複素情報(201)を把握することができる。
 図2は、実施例1における受信アレーアンテナ(102)および仮想アンテナによる、ある角度方位からの受信信号をイメージした図である。
受信アレーアンテナ(102)における各アンテナ素子には、物標の角度方位(203)からの受信信号について、時間/周波数FFT処理部(106)が複素情報(201)x1、x2、x3、x4、x5を生成する。受信信号の電磁波は波長(204)である。各アンテナ素子で受信する信号には、物標の角度方位が正面から斜めにある場合には位相差が生じる。本実施例では仮想アンテナ(205)により仮想的な受信信号を生成する。図2では実アンテナ素子からなる受信アレーアンテナ(102)の右側に仮想アンテナ(205)が配置されているが、左側に配置するようにしてもよい。
 複素情報(201)は、電磁波の到来方向すなわち物標の角度方位(203)と電磁波の波長(204)と受信アレーアンテナ(102)の配置により決まる。複素情報(201)は、受信アレーを構成する受信アンテナ素子の間で規則的に変化する。この複素情報(201)の規則的な変化により物標の角度方位(203)をもとめることができる。
 高分解能なレーダ装置を実現するためには、受信アレーアンテナ(102)の端から端の距離である開口長(202)を長くすることで実現が可能である。車載向けレーダでは、車両のバンパー内部やヘッドライト内部などに設置することが多く、設置スペースが限られており開口長(202)を長くとることが難しい。限られた開口長で高分解能を出すためのアルゴリズムとしてMUSIC法がよく知られているが、比較的計算負荷が重く、短時間で複数物標を処理する必要がある自動運転システムなどでは課題となる。
 計算負荷が軽い高分解能を実現するためには、複素情報(201)の規則的な変化を用いて、開口長(202)の外側に仮想的に拡張された仮想アンテナ(205)の複素情報を、なるべく簡単な計算により推定できればよい。簡単な計算の例としては、複素情報(201)について多項式近似を行い、仮想アンテナ(205)を推定する方法や学習器を用いた方法が挙げられる。
 具体的には、実数部の多項式近似処理部(301)と虚数部の多項式近似処理部(302)が複素情報(201)を実数部と虚数部に分離した上で、図3のように横軸をアンテナ位置(303)、縦軸を実数値(304)もしくは虚数値(305)としてプロットを実施する。
 図3は、受信アレーアンテナの受信信号における振幅/位相情報を実数/虚数部にわけてプロットし、多項式近似したときのイメージ図である。図3では、縦軸に実数値もしくは虚数値を示し、横軸にはアンテナ位置を示す。
 図3では、5つのアンテナ素子のアンテナ位置における複素情報(201)を円でプロット(306)している。このプロットした複素情報に対して多項式近似曲線(307)を開口長(202)の外側まで描くことで仮想的なアンテナ位置の複素情報を推定する。
 図3では、開口長(202)の約2.5倍の範囲まで拡張している。ここで多項式近似はプロットした複素情報(201)から拡張範囲を推定する方法の一例であり、別の推定方法として方位が既知のデータにより、機械学習により学習させた学習器を用いて拡張範囲を推定することも可能である。
 図4は、受信アレーアンテナの受信信号を空間方向に処理することで得られる角度方位スペクトルである。さらに説明すると、図4は、複素情報(201)を、空間FFT処理することで算出した物標角度方位スペクトル(401)と、多項式近似曲線(307)により拡張した仮想的な複素情報を用いて、空間FFTの処理をした物標角度方位スペクトル(402)を示す図である。縦軸は受信する電磁波について電力レベル(dB)であり、横軸は正面を0度とした角度方位(deg.)を示す。
 図4は、正面から+0.8度と-0.8度の位置に、それぞれ物標が存在した場合のシミュレーション結果であり、多項式近似曲線で開口長を2.5倍に拡張したことで分解能が改善され2つの物標が分離して検知できていることがわかる。
 図3および図4では、+/-0.8度という正面に近い物標検知を例としたが、正面から離れた角度になると隣接するアンテナ間の複素情報(201)の偏差が大きくなり多項式近似や学習器による拡張範囲の推定が困難になる。
 図5は、受信アレーアンテナの受信信号における振幅/位相情報を実数/虚数部にわけてプロットし、多項式近似したときのイメージ図である。図5では、縦軸に実数値もしくは虚数値を示し、横軸にはアンテナ位置を示す。また図5は、正面から離れた方位に物標が存在した場合の例を示す図である。図3の場合と異なり、隣接した複素情報(201)の偏差が大きいためプロット(306)の間の変化が大きくなり、多項式近似曲線(307)がうまくプロット(306)にフィットしていないことがわかる。
 この課題に対する解決策としては、近似曲線や学習器による拡張を行う前に、実アンテナの複素情報(201)を、図1における第1の空間FFT処理部(107)にて空間FFT処理を行う。この第1の角度方位推定の処理により、物標からの受信信号から、物標のおおまかな角度方位を事前に推定することができる。ここで、第1の角度方位推定処理として空間FFT処理を実行したのは、軽量な処理で角度方位が推定できるためである。第1の角度方位推定の処理として、第1の空間FFT処理部(107)の代わりに、デジタルビームフォーミングにより角度方位を推定してもよい。
 次に、第1の位相変換処理部(108)が、複素情報(201)を推定方位に位相変換する。この位相変換は、単純な位相回転処理で、この変換で推定方位が疑似的に正面となる。この位相変換により、図5で記した隣接するプロット(306)の偏差は抑圧されることになる。
 図6は、受信アレーアンテナの受信信号を空間方向に処理することで得られる角度方位スペクトルである。さらに説明すると、図6は、図5における実数部と虚数部のプロット(306)を合わせた複素情報を、第1の空間FFT処理部(107)にて空間FFT処理をして得られる物標角度方位スペクトル(601)である。図6の縦軸は電力レベル(dB)を示し、横軸は正面を0度とした角度方位(deg.)を示す。
 物標角度方位スペクトルのピーク(602)はこの場合は約3度であり、図1の第1の位相変換処理部(108)では、この角度方位3度を0度に変換する位相回転処理を行う。
 図1における第1の位相変換処理部(108)による位相変換の後は、先に述べたように、実数部の多項式近似処理部(301)と虚数部の多項式近似処理部(302)が、それぞれ、位相変換された実数部の多項式近似の処理と、位相変換された虚数の多項式近似の処理を行う。次に外挿/内挿処理部(109)と外挿/内挿処理部(112)が、それぞれ、多項式近似曲線(307)にもとづき仮想アンテナ位置での実数部の複素情報と、仮想アンテナ位置での虚数部の複素情報を算出する。
 図7は、受信アレーアンテナの受信信号における振幅/位相情報を実数/虚数部にわけてプロットし、多項式近似したときのイメージ図である。図7では、縦軸に実数値もしくは虚数値を示し、横軸にはアンテナ位置を示す。さらに説明すると、図7は、図6の角度方位3度の推定結果から複素情報(201)を0度に位相変換した後に実数部と虚数部にわけて多項式近似した結果である。
 第1の位相変換処理部(108)による位相変換をしたことで、図5では変化が大きかったプロット(306)が抑圧されて、図7の多項式近似曲線(307)はプロット(306)とうまく一致していることがわかる。
 ここで、実数部の多項式近似処理部(301)と虚数部の多項式近似処理部(302)において、多項式近似曲線(307)とプロット(306)した複素情報の偏差を、仮想アンテナにおける確度として算出し、算出した確度を利用することも可能である。
 この確度の利用例としては、確度が悪い場合には、仮想アンテナを利用せずに実アンテナのみで方位推定するように切替を行うようにできる。別の利用例としては、確度をセンサ間フュージョン処理といった後段の処理に与えることで誤検知を効率よく抑圧できる可能性がある。多項式近似曲線の代わりに、機械学習による推定結果と複素情報の偏差から、仮想アンテナにおける確度とすることも可能である。
 図7では、図3と同様に開口長(202)の2.5倍の範囲まで多項式近似曲線(307)を拡張している。図3や図7では、高分解能を得るためにアンテナ配置の外側へ拡張する外挿処理を行っているが、アンテナ配置によっては角度方位の虚像を抑圧する対策として、アンテナ間の複素情報を推定する内挿処理を行うことも容易に可能である。
 この外挿および内挿処理は、図1の外挿/内挿処理部(109)にて実行される。第2の位相変換処理部(110)において、外挿/内挿処理で算出した複素情報は、元の方位に位相変換される。その位相変換した後、第2の空間FFT処理部(111)が、受信信号および仮想的な受信信号に基づいて、空間FFT処理をし、第2の角度方位推定処理を実行する。第2の空間FFT処理部(111)は、物標(ターゲット)の角度方位、物標までの距離、物標の速度などの位置情報をレーダ出力として出力する。
 空間FFT処理に代わりデジタルビームフォーミングでもよい。空間FFT処理もしくはデジタルビームフォーミングといった第2の角度方位推定処理により、高分解能で正確な物標の角度方位を算出する。
 図8は、受信アレーアンテナの受信信号を空間方向に処理することで得られる角度方位スペクトルである。図8の縦軸は電力レベル(dB)を示し、横軸は正面を0度とした角度方位(deg.)を示す。さらに説明すると、図8は、角度方位2.2度と3.8度の位置に物標が存在する場合の空間FFT処理後の角度方位スペクトルの例である。角度方位3度の位置に1つのピークをもつ物標角度方位スペクトル(601)は、図5における、実アンテナの複素情報すなわち実数部と虚数部のプロット(306)を空間FFT処理した結果である。
 角度方位0度付近に1つのピークをもつ物標角度方位スペクトル(802)は、図5における、多項式近似曲線(307)により拡張した仮想的な複素情報を空間FFT処理した結果であり、大きく本来の物標位置から異なることがわかる。
 方位角度2.2度と3.8度の位置に2つのピークをもつ物標角度方位スペクトル(803)は、図7の多項式近似曲線(307)により拡張した仮想的な複素情報を空間FFT処理した結果であり、総じていえば図1に記した構成による処理結果だといえる。正面とは異なる角度方位で高分解能を実現している。
 実施例1によれば、高分解能処理の計算負荷を軽減するとともに、雑音耐性が強化され検知距離を延伸する効果が得られる。
 次に、実施例2について説明する。実施例1と同じ内容については、説明は省略する。
 図8の物標角度方位スペクトル(803)では、角度方位1度や5度付近にサイドローブがあり、検知閾値の設定によっては雑音環境下で誤検知となる可能性がある。サイドローブの抑圧には、第2の空間FFT処理部(111)の前の処理で窓関数処理部(1201)が窓関数処理を実行すればよい。つまり、第2の位相変換処理部(110)において、外挿/内挿処理で算出した複素情報に基づいて、受信信号と仮想的な受信信号について窓関数処理を実行すればよい。ただし、窓関数処理部(1201)が窓関数処理を行った場合には、角度方位スペクトルがなまり、分解能が劣化する。
 分解能の改善には、図2に示した仮想アンテナ(205)の領域を広げればよく、これは、図3や図7の多項式近似曲線による拡張範囲を広げることと同義である。この場合、基準となるプロット(306)から遠く離れた範囲を推定するため、正解から大きく外れた推定を行う恐れがある。
 図9は、疑似的にゼロを配置することで多項式近似の発散を防ぐことをイメージした図である。図9では、縦軸に実数値を示し、横軸にはアンテナ位置を示す。正解から大きく外れないための工夫として、図9のように、複素情報のプロット(306)の外側に、ゼロ値で疑似的なプロット(901)を配置する。このことで、プロット(306)から離れた多項式近似曲線(902)が極端に大きく外れることが回避でき、多項式近似曲線(903)のような現実的な曲線となる。
 ここで、疑似的なプロット(901)を置いたことによって、多項式近似曲線(903)の信憑性が懸念として生じる。疑似的なプロット(901)の影響を抑えるためには、多項式近似の際に複素情報のプロット(306)を重視して、疑似的なプロット(901)を軽視するように、受信信号についての複素情報のプロット(306)と疑似的なプロット(901)との間で重み付きを変更した多項式近似を実行すればよい。この場合、複素情報のプロット(306)を重ね書きした上で多項式近似曲線(903)を求めるのと同じ処理になる。
 図10は、受信アレーアンテナの受信信号における振幅/位相情報を実数/虚数部にわけてプロットし、多項式近似したときのイメージ図である。図10では、縦軸に実数値もしくは虚数値を示し、横軸にはアンテナ位置を示す。さらに説明すると、図10には、複素情報のプロット(306)の外側にゼロ値を疑似的にプロットし、重み付きでプロット(306)を重視した場合の多項式近似曲線(903)を示す。
 重み付きとしては100:1で行っており、多項式近似曲線の拡張範囲としては開口長の約20倍である。多項式近似曲線(903)は、開口長から離れた範囲でも発散することはなく、プロット(306)へのフィッティングも良好であり、正解に近いと推測できる。
 図11は、受信アレーアンテナの受信信号を空間方向に処理することで得られる角度方位スペクトルである。図1の縦軸は電力レベル(dB)を示し、横軸は正面を0度とした角度方位(deg.)を示す。さらに説明すると、図11には、図4と同じく複素情報(201)を空間FFT処理することで算出した物標角度方位スペクトル(401)と、図10で記した多項式近似曲線(903)により開口長の約20倍に拡張した仮想的な複素情報を用いて空間FFT処理した物標角度方位スペクトル(1101)を示す。
 図11は、物標が正面から+0.8度と-0.8度の位置にそれぞれ存在した場合のシミュレーション結果であり、物標角度方位スペクトル(1101)のサイドローブの抑圧と急峻なスペクトルより、雑音に強く非常に高分解能な特性が得られていることがわかる。
 図12は、実施例2におけるレーダ装置のブロック構成図である。図12には、図9と図10と図11を用いて説明した処理内容の全体構成を記載している。実施例1で説明した図1に対して追加される機能ブロックとしては、図9にて説明した、ゼロ埋めと重み付き処理部(1202)と窓関数処理部(1201)がある。
 第1の空間FFT処理部(107)にて事前に検出した推定角度方位が複数ある場合には、検出した各推定角度方位に対して第1の位相変換処理部(108)から第2の空間FFT処理部(111)による処理を、検出数に応じて繰り返しループ処理をするように繰り返しループ処理部(1203)が処理すればよい。
 実施例2によれば、角度方位スペクトルのサイドローブの抑圧と急峻なスペクトルより、雑音に強く非常に高分解能な特性を得ることができる。
101 送信アンテナ
102 受信アレーアンテナ
106 時間/周波数FFT処理部
107 第1の空間FFT処理部
108 第1の位相変換処理部
109、112 外挿/内挿処理部
110 第2の位相変換処理部
111 第2の空間FFT処理部
203 物標の角度方位
205 仮想アンテナ
301 実数部の多項式近似処理部
302 虚数部の多項式近似処理部
1201 窓関数処理部
1202 ゼロ埋めと重み付き処理部

Claims (8)

  1. 物標からの受信信号を受信する受信アンテナと、
    前記受信信号に基づいて、仮想的なアンテナにより受信した場合の仮想的な受信信号を生成する処理部と、
    前記受信信号および前記仮想的な受信信号に基づいて、前記物標の方位を算出する物標方位算出部と、を有し、
    前記処理部は、
    前記受信信号から前記物標の方位を推定する角度方位推定処理を行い、推定した前記方位に基づいた前記仮想的な受信信号を生成するレーダ装置。
  2. 請求項1に記載のレーダ装置において、
    角度方位推定処理は、
    空間FFT処理もしくはデジタルビームフォーミングであるレーダ装置。
  3. 請求項2に記載のレーダ装置において、
    前記処理部は、
    推定した前記方位に基づいて、前記受信信号の位相変換を行うレーダ装置。
  4. 請求項1に記載のレーダ装置において、
    前記処理部は、
    前記受信信号に基づいて多項式近似または機械学習を用いることで前記仮想的な受信信号を生成するレーダ装置。
  5. 請求項4に記載のレーダ装置において、
    前記処理部は、
    前記多項式近似または前記機械学習による推定結果と前記受信信号の偏差から、前記物標の方位推定の確度を算出するレーダ装置。
  6. 請求項4に記載のレーダ装置において、
    前記処理部は、
    前記受信信号の外側の範囲をゼロ埋めし、
    前記受信信号と前記ゼロ埋めした信号との間で重み付けを変更し、
    前記多項式近似または前記機械学習を用いることで前記仮想的な受信信号を生成するレーダ装置。
  7. 請求項4記載のレーダ装置において、
    前記受信信号と前記仮想的な受信信号について窓関数による処理を行うレーダ装置。
  8. 請求項1に記載のレーダ装置において、
    前記物標で反射される送信信号を出力する送信アンテナを有し、
    前記受信アンテナは、前記物標で反射した反射波を受信する受信アレーアンテナであるレーダ装置。
PCT/JP2021/007215 2020-05-27 2021-02-25 レーダ装置 WO2021240927A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
DE112021001201.0T DE112021001201T5 (de) 2020-05-27 2021-02-25 Radarvorrichtung

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2020091849A JP7382898B2 (ja) 2020-05-27 2020-05-27 レーダ装置
JP2020-091849 2020-05-27

Publications (1)

Publication Number Publication Date
WO2021240927A1 true WO2021240927A1 (ja) 2021-12-02

Family

ID=78744241

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2021/007215 WO2021240927A1 (ja) 2020-05-27 2021-02-25 レーダ装置

Country Status (3)

Country Link
JP (1) JP7382898B2 (ja)
DE (1) DE112021001201T5 (ja)
WO (1) WO2021240927A1 (ja)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001183437A (ja) * 1999-12-27 2001-07-06 Toshiba Corp 到来波の入射方位角測定装置
JP2010071865A (ja) * 2008-09-19 2010-04-02 Fujitsu Ten Ltd 信号処理装置、及びレーダ装置。
JP2017090229A (ja) * 2015-11-10 2017-05-25 富士通テン株式会社 到来方向推定装置、到来方向推定方法、到来方向推定プログラム
JP2019070558A (ja) * 2017-10-06 2019-05-09 日本無線株式会社 アレーアンテナおよびアレーアンテナの信号処理装置
US20190212411A1 (en) * 2018-01-08 2019-07-11 Hyundai Mobis Co., Ltd. Method and apparatus for estimating direction of arrival using generation of virtual received signals

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6843568B2 (ja) 2016-09-23 2021-03-17 株式会社デンソーテン レーダ装置および到来方向推定方法
JP7252052B2 (ja) 2019-05-13 2023-04-04 株式会社デンソーテン Mimoレーダ装置、および、信号処理方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001183437A (ja) * 1999-12-27 2001-07-06 Toshiba Corp 到来波の入射方位角測定装置
JP2010071865A (ja) * 2008-09-19 2010-04-02 Fujitsu Ten Ltd 信号処理装置、及びレーダ装置。
JP2017090229A (ja) * 2015-11-10 2017-05-25 富士通テン株式会社 到来方向推定装置、到来方向推定方法、到来方向推定プログラム
JP2019070558A (ja) * 2017-10-06 2019-05-09 日本無線株式会社 アレーアンテナおよびアレーアンテナの信号処理装置
US20190212411A1 (en) * 2018-01-08 2019-07-11 Hyundai Mobis Co., Ltd. Method and apparatus for estimating direction of arrival using generation of virtual received signals

Also Published As

Publication number Publication date
JP2021188943A (ja) 2021-12-13
DE112021001201T5 (de) 2023-01-12
JP7382898B2 (ja) 2023-11-17

Similar Documents

Publication Publication Date Title
JP6567832B2 (ja) レーダシステム、レーダ信号処理装置、車両走行制御装置および方法、ならびにコンピュータプログラム
CN109407093B (zh) 解决宽孔径雷达到达角模糊度的多普勒测量
JP6331195B2 (ja) レーダ装置
US9470782B2 (en) Method and apparatus for increasing angular resolution in an automotive radar system
US9041596B2 (en) Multibeam radar apparatus for vehicle, multibeam radar method and multibeam radar program
JP7117064B2 (ja) レーダ装置、レーダシステム
WO2017159735A1 (ja) 物標検出装置
WO2018034754A1 (en) Saliency based beam-forming for object detection
WO2013175558A1 (ja) レーダ装置、角度検証方法
JP2013200255A (ja) レーダ装置およびターゲット検出方法
JP5811931B2 (ja) 位相モノパルスレーダ装置
KR101953697B1 (ko) 균일 선형 배열 안테나 기반의 가상 수신신호 생성을 이용한 도래각 추정 방법 및 장치
Choi et al. High angular resolution estimation methods for vehicle FMCW radar
US11754671B2 (en) Incoming wave count estimation apparatus and incoming wave count incoming direction estimation apparatus
Li et al. Multipath ghost target identification for automotive MIMO radar
WO2021240927A1 (ja) レーダ装置
JP4232628B2 (ja) レーダ装置
KR102099388B1 (ko) 안테나 어레이 외삽을 이용한 레이더 수신신호의 도착방향 추정 방법 및 장치
JP7160561B2 (ja) 方位演算装置及び方位演算方法
Fischer et al. Evaluation of different super-resolution techniques for automotive applications
WO2021196165A1 (zh) 频率分析方法、装置及雷达
WO2020105451A1 (ja) レーダ装置、車両および物体位置検出方法
TWI808874B (zh) 用於交通工具的雷達系統及偵測方法
WO2019091783A1 (en) Detecting a parking row with a vehicle radar system
WO2024029100A1 (ja) レーダ装置

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: 21813835

Country of ref document: EP

Kind code of ref document: A1

122 Ep: pct application non-entry in european phase

Ref document number: 21813835

Country of ref document: EP

Kind code of ref document: A1