JP2023162528A - Radar device - Google Patents

Radar device Download PDF

Info

Publication number
JP2023162528A
JP2023162528A JP2022072896A JP2022072896A JP2023162528A JP 2023162528 A JP2023162528 A JP 2023162528A JP 2022072896 A JP2022072896 A JP 2022072896A JP 2022072896 A JP2022072896 A JP 2022072896A JP 2023162528 A JP2023162528 A JP 2023162528A
Authority
JP
Japan
Prior art keywords
calculation unit
eigenvalue
arrival
signal
wave number
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2022072896A
Other languages
Japanese (ja)
Other versions
JP7191262B1 (en
Inventor
祐翔 中島
Yuto Nakajima
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
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 Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP2022072896A priority Critical patent/JP7191262B1/en
Priority to PCT/JP2022/025052 priority patent/WO2023210030A1/en
Application granted granted Critical
Publication of JP7191262B1 publication Critical patent/JP7191262B1/en
Publication of JP2023162528A publication Critical patent/JP2023162528A/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00

Landscapes

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

Abstract

To obtain a radar device capable of excellently estimating the number of incoming waves even when intensity of a received signal is increased, and calculating an incoming angle with high accuracy.SOLUTION: A radar device 100 receives a reflection signal reflected by an object by a plurality of receiving antennas 6, discriminates an eigenvalue of a correlation matrix of a complex spectrum of a frequency-analyzed reflected signal into a noise eigenvalue and a signal eigenvalue by using two threshold values, estimates the number of incoming signals from the number of discriminated signal eigenvalues, and calculates an incoming angle based on the incoming signals and the estimated signal eigenvalue.SELECTED DRAWING: Figure 4

Description

本願は、レーダ装置に関するものである。 The present application relates to a radar device.

複数の受信機を用いて複数の受信信号の到来角度を求めるレーダ装置が知られている(例えば、特許文献1参照)。特許文献1においては、複数の受信機による複数の受信信号の周波数解析結果から雑音レベルを算出し、雑音レベルに基づいて受信機ごとに対象目標物体のピーク信号を抽出して共分散行列を生成し、共分散行列の各固有値を閾値により信号固有値と雑音固有値とに判別し、信号固有値の数から入射信号数(到来波数)を推定することが開示されている。 2. Description of the Related Art A radar device that uses a plurality of receivers to determine the angle of arrival of a plurality of received signals is known (see, for example, Patent Document 1). In Patent Document 1, a noise level is calculated from frequency analysis results of a plurality of received signals by a plurality of receivers, and a covariance matrix is generated by extracting a peak signal of a target object for each receiver based on the noise level. However, it is disclosed that each eigenvalue of a covariance matrix is distinguished into a signal eigenvalue and a noise eigenvalue by a threshold value, and the number of incident signals (number of arriving waves) is estimated from the number of signal eigenvalues.

特許第4351266号公報Patent No. 4351266

しかしながら、複数の受信機を備えたレーダ装置にあっては、受信機間の電磁的な結合により、あるいはレーダ装置を構成する部品により受信信号に振幅誤差、位相誤差を与える。この振幅誤差、位相誤差が固有値に影響を与え、雑音レベルと比べ受信信号強度が大きいときには、雑音レベルから求めた閾値を超過する固有値が増え、到来波数を誤推定する虞がある。 However, in a radar device equipped with a plurality of receivers, an amplitude error and a phase error are imparted to the received signal due to electromagnetic coupling between the receivers or by components constituting the radar device. This amplitude error and phase error affect the eigenvalues, and when the received signal strength is greater than the noise level, the number of eigenvalues exceeding the threshold determined from the noise level increases, leading to a risk of erroneously estimating the number of arriving waves.

そのため、受信信号強度が大きい場合であっても、受信信号の到来波数を誤推定しない技術が求められている。 Therefore, there is a need for a technique that does not erroneously estimate the number of arriving waves of a received signal even when the received signal strength is high.

本願は、上記のような課題を解決するための技術を開示するものであり、受信信号強度が小さい場合に良好な到来波数推定精度を維持しつつ、受信信号強度が大きくなった場合にも良好に到来波数を推定することができ、高精度に到来角度の算出が可能なレーダ装置を得ることを目的とする。 This application discloses a technique for solving the above-mentioned problems, which maintains good arriving wave number estimation accuracy when the received signal strength is low, and also achieves good accuracy when the received signal strength increases. The present invention aims to provide a radar device that can estimate the number of arriving waves and calculate the angle of arrival with high accuracy.

本願に開示されるレーダ装置は、
送信信号を送信する送信アンテナと、
単数又は複数の物体で反射された反射信号を受信する複数の受信アンテナと、
複数の前記受信アンテナで受信した前記反射信号を周波数解析し、前記物体のそれぞれについて、自装置との角度である到来角度を算出する制御部と、を備えたレーダ装置であって、
前記制御部は、
前記周波数解析により得られた複素スペクトラムの振幅値から雑音レベルを演算する雑音レベル演算部と、
前記周波数解析により得られた前記複素スペクトラムの振幅値及び前記雑音レベル演算部で演算された前記雑音レベルに基づいてピーク値を判定するピーク値判定部と、
前記ピーク値と判定された複素スペクトラムに対する相関行列を演算する相関行列演算部と、
前記相関行列の固有値を演算する固有値演算部と、
前記固有値を雑音に起因する雑音固有値と前記物体に対応する信号固有値とに判別するための閾値を演算する到来波数推定閾値演算部と、
前記到来波数推定閾値演算部で演算された閾値に基づいて、判別された前記信号固有値の数から到来信号の数を推定し、前記到来信号と推定された前記信号固有値に基づいて前記到来角度を演算する到来角度演算部と、を有し、
前記到来波数推定閾値演算部は、
前記雑音レベル演算部で演算された雑音レベルに基づいて第1到来波数推定閾値を演算するとともに、前記固有値演算部で演算された前記固有値に基づいて第2到来波数推定閾値を演算し、
前記到来角度演算部は、前記第1到来波数推定閾値及び前記第2到来波数推定閾値を用いて判別された前記信号固有値の数から前記到来信号の数を推定する、ように構成されている。
The radar device disclosed in this application includes:
a transmitting antenna for transmitting a transmitting signal;
a plurality of receiving antennas that receive reflected signals reflected by one or more objects;
A radar device comprising: a control unit that frequency-analyzes the reflected signals received by a plurality of the receiving antennas and calculates an angle of arrival, which is an angle with respect to the own device, for each of the objects,
The control unit includes:
a noise level calculation unit that calculates a noise level from the amplitude value of the complex spectrum obtained by the frequency analysis;
a peak value determination unit that determines a peak value based on the amplitude value of the complex spectrum obtained by the frequency analysis and the noise level calculated by the noise level calculation unit;
a correlation matrix calculation unit that calculates a correlation matrix for the complex spectrum determined to be the peak value;
an eigenvalue calculation unit that calculates an eigenvalue of the correlation matrix;
an arriving wave number estimation threshold calculation unit that calculates a threshold for discriminating the eigenvalue into a noise eigenvalue caused by noise and a signal eigenvalue corresponding to the object;
The number of arriving signals is estimated from the number of determined signal eigenvalues based on the threshold calculated by the arriving wave number estimation threshold calculation unit, and the arrival angle is determined based on the incoming signal and the estimated signal eigenvalue. an arrival angle calculation unit that calculates,
The arriving wave number estimation threshold calculation unit includes:
calculating a first arrival wave number estimation threshold based on the noise level calculated by the noise level calculation unit, and calculating a second arrival wave number estimation threshold based on the eigenvalue calculated by the eigenvalue calculation unit;
The arrival angle calculation unit is configured to estimate the number of arriving signals from the number of signal eigenvalues determined using the first arriving wave number estimation threshold and the second arrival wave number estimation threshold.

本願に開示されるレーダ装置によれば、2つの閾値を用いて信号固有値を判別するので、受信信号強度が小さい場合に良好な到来波数推定精度を維持しつつ、受信信号強度が大きくなった場合にも良好に到来波数を推定することができ、高精度に到来角度の算出が可能なレーダ装置を得ることが可能となる。 According to the radar device disclosed in the present application, since the signal eigenvalue is determined using two threshold values, good incoming wave number estimation accuracy is maintained when the received signal strength is small, and when the received signal strength becomes large. Therefore, it is possible to obtain a radar device that can estimate the number of arriving waves well and can calculate the angle of arrival with high accuracy.

実施の形態1に係るレーダ装置の概略構成図である。1 is a schematic configuration diagram of a radar device according to Embodiment 1. FIG. 実施の形態1に係るレーダ装置の信号処理部の機能ブロック図である。FIG. 3 is a functional block diagram of a signal processing unit of the radar device according to the first embodiment. 実施の形態1に係るレーダ装置の制御部のハードウエア構成図である。FIG. 3 is a hardware configuration diagram of a control unit of the radar device according to the first embodiment. 実施の形態1に係るレーダ装置における到来角度を算出する処理手順を示すフローチャートである図である。FIG. 3 is a flowchart illustrating a processing procedure for calculating an angle of arrival in the radar device according to the first embodiment. 実施の形態1に係るレーダ装置における複素スペクトラムの振幅演算結果を示した図である。FIG. 3 is a diagram showing amplitude calculation results of a complex spectrum in the radar device according to the first embodiment. 図5Aのうち、距離周波数と振幅との関係を示す図である。FIG. 5A is a diagram showing the relationship between distance frequency and amplitude in FIG. 5A. 図5Aのうち、相対速度周波数と振幅との関係を示す図である。FIG. 5A is a diagram showing the relationship between relative velocity frequency and amplitude in FIG. 5A. 実施の形態1に係るレーダ装置における相関行列演算時のサブアレイを示した概念図である。FIG. 2 is a conceptual diagram showing subarrays during correlation matrix calculation in the radar device according to the first embodiment. 実施の形態1に係るレーダ装置における固有値と到来角度との関係を示した図である。FIG. 3 is a diagram showing the relationship between eigenvalues and arrival angles in the radar device according to the first embodiment. 実施の形態1に係るレーダ装置における到来角度を算出する別の処理手順を示すフローチャートである図である。7 is a flowchart illustrating another processing procedure for calculating the angle of arrival in the radar device according to the first embodiment. FIG. 実施の形態1に係るレーダ装置における固有値及び到来波数推定閾値と受信信号強度との関係を示した図である。FIG. 3 is a diagram showing the relationship between the eigenvalue, the arrival wave number estimation threshold, and the received signal strength in the radar device according to the first embodiment.

以下、本願で開示されるレーダ装置の実施の形態について図を参照して説明する。なお、本実施の形態に係るレーダ装置は車両に搭載され、車両の周囲に存在する物体を検知することを一例として想定しているが、車両以外の装置(例えば、航空機などの移動体、監視装置等)に搭載されてもよい。各図中、同一符号は、同一または相当部分を示すものとする。 Embodiments of the radar device disclosed in this application will be described below with reference to the drawings. Note that the radar device according to this embodiment is mounted on a vehicle and is assumed to detect objects existing around the vehicle, but it is assumed that the radar device according to the present embodiment is installed in a vehicle and detects objects existing around the vehicle. device, etc.). In each figure, the same reference numerals indicate the same or corresponding parts.

実施の形態1.
以下に、実施の形態1に係るレーダ装置について図を用いて説明する。
図1は、実施の形態1に係るレーダ装置の概略構成図である。レーダ方式には、FM-CW(Frequency Modulated Continuous Wave)方式、FCM(Fast Chirp Modulation)方式、パルス・ドップラー方式などが挙げられるが、本実施の形態1はどの方式にも適用可能である。本実施の形態1では、一例としてFCM方式の周波数変調レーダ装置100(以下、単に「レーダ装置」と称する)について説明する。
Embodiment 1.
Below, a radar device according to Embodiment 1 will be explained using figures.
FIG. 1 is a schematic configuration diagram of a radar device according to a first embodiment. Examples of the radar system include an FM-CW (Frequency Modulated Continuous Wave) system, an FCM (Fast Chirp Modulation) system, and a pulse Doppler system, but the first embodiment is applicable to any of the systems. In the first embodiment, an FCM frequency modulation radar device 100 (hereinafter simply referred to as "radar device") will be described as an example.

図1において、レーダ装置100は、6チャネル分の受信信号として反射信号を受信する場合を例示しているが6チャネルに限るものではない。レーダ装置100は、信号処理部1、制御電圧発生器2、VCO(Voltage Controlled Oscillator;発振器)3、分配器4、送信アンテナ5、受信アンテナ6、ミキサ7、A/Dコンバータ8、及びFFT演算部(Fast Fourier Transform;高速フーリエ変換)9を備えている。ここで、A/Dコンバータ8、FFT演算部9及び信号処理部1で制御部10を構成する。なお、送信アンテナ5の数は複数としても良い。また、受信アンテナ6、ミキサ7、A/Dコンバータ8、及びおよびFFT演算部9はそれぞれ6チャネルに対応しているが、例えば受信アンテナ6は総称する場合に受信アンテナ6と表記し、個々を示す場合は受信アンテナ6_n(nは1から6の自然数)と記す。ミキサ7、A/Dコンバータ8、及びおよびFFT演算部9も同様の表記とする。 In FIG. 1, a case is illustrated in which the radar device 100 receives reflected signals as reception signals for six channels, but the number of channels is not limited to six. The radar device 100 includes a signal processing section 1, a control voltage generator 2, a VCO (Voltage Controlled Oscillator) 3, a distributor 4, a transmitting antenna 5, a receiving antenna 6, a mixer 7, an A/D converter 8, and an FFT calculation. (Fast Fourier Transform) 9. Here, the A/D converter 8, the FFT calculation section 9, and the signal processing section 1 constitute a control section 10. Note that the number of transmitting antennas 5 may be plural. In addition, the receiving antenna 6, mixer 7, A/D converter 8, and FFT calculation unit 9 each correspond to 6 channels, but for example, the receiving antenna 6 is referred to collectively as receiving antenna 6, and individually When shown, it is written as receiving antenna 6_n (n is a natural number from 1 to 6). The mixer 7, A/D converter 8, and FFT calculation section 9 are also expressed in the same manner.

次に、本実施の形態1に係るレーダ装置100の動作について説明する。まず始めに、信号処理部1から、制御電圧発生器2に対して変調開始命令が出力される。制御電圧発生器2は、変調開始命令に応じて、所望の電圧波形を発生し、VCO3に印加する。VCO3は、制御電圧に従って周波数変調を施した送信信号を出力する。出力された送信信号は分配器4により、送信アンテナ5および6個のミキサ7_1~7_6に分配される。 Next, the operation of radar device 100 according to the first embodiment will be explained. First, the signal processing section 1 outputs a modulation start command to the control voltage generator 2. Control voltage generator 2 generates a desired voltage waveform and applies it to VCO 3 in response to the modulation start command. The VCO 3 outputs a frequency-modulated transmission signal according to the control voltage. The output transmission signal is distributed by the distributor 4 to the transmission antenna 5 and six mixers 7_1 to 7_6.

送信アンテナ5は、目標物体(図示せず)に向けて送信信号を放射する。目標物体で反射された信号は、6個の受信アンテナ6_1~6_6により、6チャネル分(CH1~CH6とする)の受信信号として受信される。なお、複数の送信アンテナを使用する場合は、受信信号は送信アンテナ数と受信アンテナ数の積で表される数だけ取得することができる。この技術はMIMO(Multiple Input Multiple Output:マルチ入力マルチ出力)技術と呼ばれ、MIMO技術により形成される仮想受信アンテナも単に「受信アンテナ」と称する。すなわち、以下の実施の形態では、「受信アンテナ」は実配置された「受信アンテナ」に限らず「仮想受信アンテナ」を含んでいてもよい。 Transmission antenna 5 radiates a transmission signal toward a target object (not shown). The signals reflected by the target object are received by six receiving antennas 6_1 to 6_6 as received signals for six channels (referred to as CH1 to CH6). Note that when using a plurality of transmitting antennas, it is possible to obtain as many received signals as the product of the number of transmitting antennas and the number of receiving antennas. This technology is called MIMO (Multiple Input Multiple Output) technology, and the virtual receiving antenna formed by MIMO technology is also simply referred to as "receiving antenna." That is, in the following embodiments, the "receiving antenna" is not limited to the actually arranged "receiving antenna" but may include a "virtual receiving antenna".

6個の受信アンテナ6_1~6_6に対応してそれぞれ個別に設けられたミキサ7_1~7_6により、受信アンテナ6_1~6_6で受信されたそれぞれの受信信号と、分配器4により分配された送信信号とがミキシングされ、6チャネル分のビート信号が生成される。 Mixers 7_1 to 7_6 provided individually corresponding to the six receiving antennas 6_1 to 6_6 mix the respective received signals received by the receiving antennas 6_1 to 6_6 and the transmitted signals distributed by the distributor 4. The signals are mixed and beat signals for six channels are generated.

ミキシングされた各ビート信号はA/Dコンバータ8_1~8_6により、それぞれディジタルデータに変換される。 Each mixed beat signal is converted into digital data by A/D converters 8_1 to 8_6, respectively.

A/Dコンバータ8_1~8_6で変換されたディジタルデータは、FFT演算部9_1~9_6により、それぞれFFTを用いた周波数解析が施される。FFT演算部9_1~9_6により算出され、6チャネル分の周波数解析結果である複素スペクトラムは、信号処理部1に入力される。 The digital data converted by the A/D converters 8_1 to 8_6 is subjected to frequency analysis using FFT by FFT calculation units 9_1 to 9_6, respectively. The complex spectra calculated by the FFT calculation units 9_1 to 9_6 and which are the frequency analysis results for six channels are input to the signal processing unit 1.

図2は、信号処理部1の構成を示す機能ブロック図である。信号処理部1では、算出された複素スペクトラムに基づいて、目標物体の到来角度を演算する。本実施の形態では、各受信アンテナ6が有する振幅誤差、位相誤差が固有値に与える影響を考慮した到来波数推定閾値を設定することにより、到来波数を精度よく推定できる点を技術的特徴としている。そのため、図2に示すように、信号処理部1は、複素スペクトラムの振幅を演算する振幅演算部101、演算された複素スペクトラムの振幅における雑音レベルを演算する雑音レベル演算部102、演算された複素スペクトラムの振幅からピーク値を判定するピーク値判定部103、判定されたピーク値から目標物体の距離及び相対速度を求める距離及び相対速度演算部104、複数の受信アンテナからサブアレイを形成し、ピークと判定された複素スペクトラムの相関行列(共分散行列)を演算する相関行列演算部105、相関行列から固有値を演算する固有値演算部106、固有値を信号固有値と雑音固有値とに判別するための閾値を演算する到来波数推定閾値演算部107、到来波数推定閾値に基づいて信号固有値と雑音固有値とに判別する固有値判定部108、及び信号固有値に基づいて目標物体の到来波数を推定し到来角度を演算する到来角度演算部109を備える。この信号処理部1の各機能部101~109による具体的な処理については後述する。 FIG. 2 is a functional block diagram showing the configuration of the signal processing section 1. As shown in FIG. The signal processing unit 1 calculates the arrival angle of the target object based on the calculated complex spectrum. A technical feature of this embodiment is that the number of arriving waves can be estimated accurately by setting a threshold for estimating the number of arriving waves in consideration of the influence of the amplitude error and phase error of each receiving antenna 6 on the eigenvalue. Therefore, as shown in FIG. 2, the signal processing unit 1 includes an amplitude calculation unit 101 that calculates the amplitude of a complex spectrum, a noise level calculation unit 102 that calculates the noise level in the amplitude of the calculated complex spectrum, and a noise level calculation unit 102 that calculates the noise level in the amplitude of the calculated complex spectrum. A peak value determination unit 103 determines the peak value from the amplitude of the spectrum; a distance and relative velocity calculation unit 104 determines the distance and relative velocity of the target object from the determined peak value; a subarray is formed from a plurality of receiving antennas; A correlation matrix calculation unit 105 calculates a correlation matrix (covariance matrix) of the determined complex spectrum, an eigenvalue calculation unit 106 calculates eigenvalues from the correlation matrix, and a threshold value is calculated for discriminating eigenvalues into signal eigenvalues and noise eigenvalues. an incoming wave number estimation threshold calculation unit 107 for estimating the number of arriving waves; an eigenvalue determination unit 108 for discriminating between signal eigenvalues and noise eigenvalues based on the incoming wave number estimation threshold; An angle calculation unit 109 is provided. Specific processing by each of the functional units 101 to 109 of the signal processing unit 1 will be described later.

図3は、レーダ装置100の制御部10のハードウエアの構成例を示す図である。図10に示すように、制御部10は処理回路として、DSP(Digital Signal Processor)等の演算処理装置11(コンピュータ)、演算処理装置11とデータのやり取りをする記憶装置12、演算処理装置11にビート信号を入力するA/Dコンバータ8、演算処理装置11から制御電圧発生器2に周波数変調信号の指令値を出力するD/Aコンバータ14、及び通信回路13等を備えている。 FIG. 3 is a diagram showing an example of the hardware configuration of the control unit 10 of the radar device 100. As shown in FIG. 10, the control unit 10 serves as a processing circuit, and includes an arithmetic processing unit 11 (computer) such as a DSP (Digital Signal Processor), a storage device 12 that exchanges data with the arithmetic processing unit 11, and an arithmetic processing unit 11. It includes an A/D converter 8 that inputs a beat signal, a D/A converter 14 that outputs a command value of a frequency modulation signal from an arithmetic processing unit 11 to a control voltage generator 2, a communication circuit 13, and the like.

演算処理装置11として、CPU(Central Processing Unit)、ASIC(Application Specific Integrated Circuit)、IC(Integrated Circuit)、FPGA(Field Programmable Gate Array)、各種の論理回路、及び各種の信号処理回路等が備えられてもよい。また、演算処理装置11として、同じ種類のもの又は異なる種類のものが複数備えられ、各処理が分担して実行されてもよい。記憶装置12として、演算処理装置11からデータを読み出し及び書き込みが可能に構成されたRAM(Random Access Memory)、及び演算処理装置11からデータを読み出し可能に構成されたROM(Read Only Memory)等が備えられている。通信回路13は、例えば、レーダ装置100が車両に搭載されている場合、車両制御装置等の外部の制御装置(図示せず)と、通信線を介して接続され、CAN(Controller Area Network(登録商標))等の通信プロトコルに基づいて有線通信を行う。 The arithmetic processing unit 11 includes a CPU (Central Processing Unit), an ASIC (Application Specific Integrated Circuit), an IC (Integrated Circuit), an FPGA (Field Programmable Gate Array), various logic circuits, various signal processing circuits, etc. You can. Further, a plurality of arithmetic processing units 11 of the same type or different types may be provided, and each process may be shared and executed. As the storage device 12, a RAM (Random Access Memory) configured to be able to read and write data from the arithmetic processing unit 11, a ROM (Read Only Memory) configured to be able to read data from the arithmetic processing unit 11, etc. It is equipped. For example, when the radar device 100 is mounted on a vehicle, the communication circuit 13 is connected to an external control device (not shown) such as a vehicle control device via a communication line, and is connected to a CAN (Controller Area Network). Wired communication is performed based on communication protocols such as (Trademark)).

そして、制御部10が備えるFFT演算部9、送信信号生成部(図示せず)、信号処理部1の各機能部101~109は、演算処理装置11が、ROM等の記憶装置12に記憶されたソフトウエア(プログラム)を実行し、記憶装置12、A/Dコンバータ8、D/Aコンバータ14、及び通信回路13等の制御部10の各ハードウエアと協働することにより実現される。なお、信号処理部1が用いる設定データは、ソフトウエア(プログラム)の一部として、ROM等の記憶装置12に記憶されている。 The functional units 101 to 109 of the FFT calculation unit 9, the transmission signal generation unit (not shown), and the signal processing unit 1 included in the control unit 10 are stored in a storage device 12 such as a ROM. This is realized by executing the software (program) and cooperating with each hardware of the control unit 10 such as the storage device 12, the A/D converter 8, the D/A converter 14, and the communication circuit 13. Note that the setting data used by the signal processing section 1 is stored in a storage device 12 such as a ROM as part of software (program).

次に、信号処理部1による処理、すなわち到来波数を精度よく推定し、それに基づいて到来角度を算出する具体的な処理について、フローチャートに基づいて説明する。
図4は、本実施の形態1に係るレーダ装置100の信号処理部1における到来角度算出の処理を表すフローチャートである。
まず、ステップS1において、振幅演算部101は、FFT演算部9から入力された6チャネル分のビート信号の複素スペクトラムのそれぞれについて、振幅を求める。
Next, specific processing by the signal processing unit 1, that is, estimating the number of arriving waves with high accuracy and calculating an angle of arrival based on the estimation, will be described based on a flowchart.
FIG. 4 is a flowchart showing the process of calculating the angle of arrival in the signal processing unit 1 of the radar device 100 according to the first embodiment.
First, in step S1, the amplitude calculation unit 101 calculates the amplitude of each of the complex spectra of the six channels of beat signals input from the FFT calculation unit 9.

図5A~図5Cは、実施の形態1に係るレーダ装置100の振幅演算部101における複素スペクトラムの振幅演算結果を示した図である。図5Aは、距離周波数及び相対速度周波数に対応する振幅を示す図、図5Bは、図5Aのうち、距離周波数と振幅との関係を示す図、図5Cは、図5Aのうち、相対速度周波数と振幅との関係を示す図である。チャネルごとに、各複素スペクトラムの絶対値の二乗を計算することにより、図5A~図5Cに示したような、距離周波数、相対速度周波数に対応する振幅値を得ることができる。 5A to 5C are diagrams showing amplitude calculation results of a complex spectrum in amplitude calculation section 101 of radar device 100 according to the first embodiment. 5A is a diagram showing the amplitude corresponding to the distance frequency and the relative velocity frequency, FIG. 5B is a diagram showing the relationship between the distance frequency and the amplitude in FIG. 5A, and FIG. 5C is a diagram showing the relative velocity frequency in FIG. 5A. FIG. 3 is a diagram showing the relationship between and amplitude. By calculating the square of the absolute value of each complex spectrum for each channel, amplitude values corresponding to the distance frequency and relative velocity frequency as shown in FIGS. 5A to 5C can be obtained.

次にステップS2において、雑音レベル演算部102は、雑音振幅値の代表値として雑音レベルPを演算する。ただし、雑音振幅値の代表値は平均値算出、中央値算出、最頻値算出等によって得られるものとする。具体的な演算例としては、先のステップS1で得られた振幅値を全チャネルで加算し、距離周波数ごとに相対速度周波数軸上の振幅値の平均値Pを求める。Pに予め定めた値を乗じた値を仮閾値Tとし、仮閾値T以下の振幅値が雑音であるとみなし、改めて雑音のみの平均値P′を求める。上述の計算を繰り返し、収束した平均値を雑音レベルPとする。この場合、雑音レベルPは距離周波数ごとに算出される。 Next, in step S2, the noise level calculating section 102 calculates a noise level Pn as a representative value of the noise amplitude values. However, the representative value of the noise amplitude value shall be obtained by calculating the average value, median value, mode value, etc. As a specific calculation example, the amplitude values obtained in the previous step S1 are added for all channels, and the average value P a of the amplitude values on the relative velocity frequency axis is determined for each distance frequency. A value obtained by multiplying P a by a predetermined value is set as a temporary threshold T t , amplitude values below the temporary threshold T t are considered to be noise, and an average value P a ′ of only the noise is determined again. The above calculation is repeated and the converged average value is defined as the noise level Pn . In this case, the noise level P n is calculated for each distance frequency.

次にステップS3において、ピーク値判定部103は、先のステップS1で得られた振幅値の中からピーク信号を判定する。まず、ピーク信号判別用閾値Athを式(1)のように算出する。
th=α・P ・・・(1)
ここで、αは予め定められた誤警報確率によって決まる定数である。
thよりも大きな振幅値であり、かつ前後の周波数の振幅値よりも大きいものを目標物体に対応するピーク値と判定する。
Next, in step S3, the peak value determination unit 103 determines a peak signal from among the amplitude values obtained in the previous step S1. First, the peak signal discrimination threshold value A th is calculated as shown in equation (1).
A th = α・P n ...(1)
Here, α is a constant determined by a predetermined false alarm probability.
An amplitude value that is larger than A th and larger than the amplitude values of the preceding and following frequencies is determined to be the peak value corresponding to the target object.

次にステップS4において、距離及び相対速度演算部104は、ステップS3で得られたピーク値に対応する距離周波数及び相対速度周波数から、目標物体の距離及び相対速度を演算する。さらに、検出された目標物体の数をカウントし、記憶装置12に記憶しておく。 Next, in step S4, the distance and relative velocity calculation unit 104 calculates the distance and relative velocity of the target object from the distance frequency and relative velocity frequency corresponding to the peak value obtained in step S3. Furthermore, the number of detected target objects is counted and stored in the storage device 12.

次にステップS5において、相関行列演算部105は、ステップS4においてピーク値と判定された複素スペクトラムに対して相関行列を演算する。複数の受信アンテナ6から任意形状のサブアレイをM個抽出し、それらを平均化し、平均化された相関行列Rxxを次の式(2)により演算する。ここで任意形状とは、複数の受信アンテナ6から異なるアンテナの組でサブアレイを形成することを言う。また、サブアレイを形成する手法は空間平均を算出する方法であり、この手法に限らず、異なる時刻の複素スペクトラムを用いた時間平均であってもよい。例えば、同一物体から異なる時刻に受信された信号から周波数解析によって得られ、ピーク値と判定された複素スペクトラムに対応する相関行列をそれぞれ演算するとともにそれらを平均化して、平均化された相関行列を用いればよい。 Next, in step S5, the correlation matrix calculation unit 105 calculates a correlation matrix for the complex spectrum determined to have a peak value in step S4. M arbitrary-shaped subarrays are extracted from a plurality of receiving antennas 6, and the subarrays are averaged, and an averaged correlation matrix Rxx is calculated using the following equation (2). Here, the arbitrary shape means that a subarray is formed by a set of different antennas from the plurality of receiving antennas 6. Further, the method of forming the subarray is a method of calculating a spatial average, and is not limited to this method, but may be a time average using complex spectra at different times. For example, correlation matrices corresponding to complex spectra determined to be peak values obtained by frequency analysis from signals received from the same object at different times are computed, and the averaged correlation matrices are calculated by averaging them. Just use it.

Figure 2023162528000002

ここで、xは複素スペクトラムからなるベクトルで、式(3)のように表される。
Figure 2023162528000002

Here, x m is a vector consisting of a complex spectrum, and is expressed as in equation (3).

Figure 2023162528000003

Kはサブアレイ内の複素スペクトラムの総数であり、Hは複素共役転置、Tは転置を示す。
Figure 2023162528000003

K is the total number of complex spectra in the subarray, H indicates complex conjugate transposition, and T indicates transposition.

ここで、サブアレイ抽出について例を挙げて説明する。図6は、本実施の形態1に係るレーダ装置100における相関行列算出時のサブアレイ抽出を示す概念図である。図1で示した受信アンテナ6_1~6_6に対し、サブアレイ_1は、受信アンテナ6_1から受信アンテナ6_4で構成され、サブアレイ_2は、受信アンテナ6_2から受信アンテナ6_5で構成され、サブアレイ_3は受信アンテナ6_3から受信アンテナ6_6で構成される例である。従って、図5では、サブアレイの個数M=3、サブアレイ内の複素スペクトラムの総数K=4である(M、Kはそれぞれ自然数)。このサブアレイの形成には、上述したように受信アンテナ6に限らず仮想受信アンテナを含んでいてもよい。 Here, subarray extraction will be explained using an example. FIG. 6 is a conceptual diagram showing subarray extraction during correlation matrix calculation in radar device 100 according to the first embodiment. Regarding the receiving antennas 6_1 to 6_6 shown in FIG. 1, subarray_1 is composed of receiving antennas 6_1 to 6_4, subarray_2 is composed of receiving antennas 6_2 to 6_5, and subarray_3 is composed of receiving antennas 6_3 to 6_5. This is an example configured with a receiving antenna 6_6. Therefore, in FIG. 5, the number of subarrays M=3, and the total number of complex spectra in the subarrays K=4 (M and K are each natural numbers). Formation of this subarray may include not only the receiving antenna 6 but also a virtual receiving antenna as described above.

次にステップS6において、固有値演算部106は、相関行列Rxxについて固有値及び固有ベクトルの解析を行い、固有値λ~λを求める。ただし、固有値は式(4)のように大きさが降順になるように並び替えられている。
λ ≧ λ ≧・・・≧ λ ・・・(4)
Next, in step S6, the eigenvalue calculation unit 106 analyzes the eigenvalues and eigenvectors of the correlation matrix R xx to obtain eigenvalues λ 1 to λ K. However, the eigenvalues are sorted in descending order of magnitude as shown in equation (4).
λ 1 ≧ λ 2 ≧...≧ λ K ... (4)

雑音レベルに相当する固有値λは相関行列Rxxの生成方法に基づく係数cを用いて式(5)のように表すことができる。
λ =c・P ・・・(5)
The eigenvalue λ n corresponding to the noise level can be expressed as in equation (5) using a coefficient cm based on the method of generating the correlation matrix R xx .
λ n = cm・P n ...(5)

また、反射信号の信号強度が大きい場合における1波入射時の各固有値λ′(θ)は、アンテナの振幅誤差、位相誤差により式(6)のようにλ′(θ)の定数倍r(θ)として表すことができる。

Figure 2023162528000004
In addition, when the signal strength of the reflected signal is large, each eigenvalue λ k ′(θ) at the time of one wave incidence is a constant multiple of λ 1 ′(θ) due to the amplitude error and phase error of the antenna, as shown in equation (6). It can be expressed as r k (θ).
Figure 2023162528000004

ここで、θは到来角度を示す。図7は、固有値と到来角度との関係を示した図である。式(6)で示したr(θ)は到来角度によって異なる。 Here, θ indicates the angle of arrival. FIG. 7 is a diagram showing the relationship between eigenvalues and arrival angles. r k (θ) shown in equation (6) varies depending on the arrival angle.

次にステップS7において、到来波数推定閾値演算部107は、第1到来波数推定閾値λth1を演算する。具体的には、式(5)を用いて第1到来波数推定閾値λth1を式(7)のように演算する。
λth1 =λ +λ ・・・(7)
ここで、λ は式(8)で表される。
λ =c・P ・・・(8)
は相関行列Rxxの生成方法に基づく係数である。
すなわち、
λth1 =c・P +c・P =(c+c)P
であり、第1到来波数推定閾値λth1は、雑音レベルPに予め設定された値(c+c)を乗算した値である。
Next, in step S7, the arriving wave number estimation threshold calculation unit 107 calculates the first arriving wave number estimation threshold λ th1 . Specifically, the first arriving wave number estimation threshold λ th1 is calculated as shown in equation (7) using equation (5).
λ th1 = λ n + λ s ...(7)
Here, λ s is expressed by equation (8).
λ s =c s・P n ...(8)
cs is a coefficient based on the method of generating the correlation matrix Rxx .
That is,
λ th1 = cm・P n +c s・P n =(c m +c s )P n
The first incoming wave number estimation threshold λ th1 is a value obtained by multiplying the noise level P n by a preset value (c m +c s ).

式(7)において、相関行列Rxxに関連するλを加算することで、雑音レベルの時間ばらつき、チャネル間ばらつきを考慮してλth1が設定される。式(7)の第1到来波数推定閾値λth1を使用することで、固有値判定部108において、雑音レベルに基づいて、固有値を信号固有値と雑音固有値とに判別することができる。 In Equation (7), by adding λ s related to the correlation matrix R xx , λ th1 is set in consideration of the time variation of the noise level and the inter-channel variation. By using the first incoming wave number estimation threshold λ th1 in equation (7), the eigenvalue determination unit 108 can discriminate the eigenvalue into a signal eigenvalue and a noise eigenvalue based on the noise level.

次に、ステップS8において、到来波数推定閾値演算部107は、第2到来波数推定閾値λth2を演算する。具体的には、式(6)を用いて式(9)のように演算する。
λth2 =r′・λ ・・・(9)
ここで、r′は図7で示すとおり、ある到来角度範囲におけるrの最大の値とする。式(9)の第2到来波数推定閾値λth2は信号強度が大きい領域での到来波数推定閾値に相当するため、この領域において、第2到来波数推定閾値λth2を使用することで、アンテナの振幅誤差、位相誤差を考慮して、固有値を信号固有値と雑音固有値とに判別することができる。
Next, in step S8, the arriving wave number estimation threshold calculating section 107 calculates a second arriving wave number estimation threshold λ th2 . Specifically, equation (6) is used to calculate as shown in equation (9).
λ th2 = r 2 ′・λ 1 ...(9)
Here, r 2 ' is the maximum value of r 2 in a certain arrival angle range, as shown in FIG. Since the second arrival wave number estimation threshold λ th2 in equation (9) corresponds to the arrival wave number estimation threshold in a region where the signal strength is high, by using the second arrival wave number estimation threshold λ th2 in this region, the antenna The eigenvalues can be distinguished into signal eigenvalues and noise eigenvalues by considering the amplitude error and phase error.

次に、ステップS9において、固有値判定部108は、第1到来波数推定閾値λth1及び第2到来波数推定閾値λth2を用い、固有値を信号固有値と雑音固有値とに判別する。すなわち、それぞれの閾値より大きい固有値を信号固有値と判別し、その判別された信号固有値の数から到来波数、すなわち到来信号の数を推定する。 Next, in step S9, the eigenvalue determination unit 108 determines the eigenvalue into a signal eigenvalue and a noise eigenvalue using the first arrival wave number estimation threshold λ th1 and the second arrival wave number estimation threshold λ th2 . That is, the eigenvalues larger than the respective thresholds are determined to be signal eigenvalues, and the number of arriving waves, that is, the number of arriving signals is estimated from the number of determined signal eigenvalues.

次に、ステップS10において、固有値判定部108で判定された信号固有値の数から到来波数、すなわち到来信号の数を推定する。到来角度演算部109は、推定された到来信号に対応する固有値を用い、MUSIC(Multiple Signal Classification)あるいはESPRIT(Estimation Signal Parameters via Rotational Invariance Technique)等の高分解能到来方向推定法を実施することにより、到来角度を高精度に演算することが可能となる。従って、受信信号強度が極めて小さい場合、あるいは大きい場合において精度の高い到来波数推定結果を得ることができ、その到来波推定結果を用いて到来角度を高精度に算出することが可能となる。 Next, in step S10, the number of arriving waves, that is, the number of arriving signals, is estimated from the number of signal eigenvalues determined by the eigenvalue determining section 108. The arrival angle calculation unit 109 uses the eigenvalues corresponding to the estimated arrival signals to implement a high-resolution direction of arrival estimation method such as MUSIC (Multiple Signal Classification) or ESPRIT (Estimation Signal Parameters via Rotational Invariance Technique). It becomes possible to calculate the angle of arrival with high precision. Therefore, it is possible to obtain a highly accurate estimation result of the number of arriving waves when the received signal strength is extremely small or large, and it is possible to calculate the angle of arrival with high accuracy using the estimation result of the arriving wave.

また、到来波数推定閾値演算部107において、第2到来波数推定閾値λth2を演算する場合、到来角度θを予め見積もる(推定する)ことが可能な場合は、その角度に応じてr′を変更することでより高い到来波数推定精度を実現できる。 In addition, when calculating the second arrival wave number estimation threshold λ th2 in the arriving wave number estimation threshold calculation unit 107, if it is possible to estimate (estimate) the arrival angle θ in advance, r 2 ' is calculated according to the angle. By changing this, higher accuracy in estimating the number of arriving waves can be achieved.

到来角度を予め推定する方法としては、例えばビームフォーマー(Beam former)法のような到来波数情報を必要としない到来方向推定法を用いて得られる結果を暫定到来角度として用いる。あるいは、同物体の到来角度の前回値を暫定到来角度として用いてもよい。 As a method for estimating the angle of arrival in advance, a result obtained using a direction of arrival estimation method that does not require information on the number of arriving waves, such as a beam former method, is used as the provisional angle of arrival. Alternatively, the previous value of the arrival angle of the same object may be used as the provisional arrival angle.

′の変更例として、雑音レベルが極めて小さい場合における1波入射時を考える。到来角度θ=θのとき、第2到来波数推定閾値λth2(θ)は、式(9)から式(10)にように得られる。
λth2(θ)=r(θ)・λ ≦ r′・λ・・・(10)
従って、r′よりもr(θ)を使用したほうが閾値を小さく設定することができるため、到来角度θ付近に2波以上入射する場合の到来波数推定精度が高くなる。
As an example of changing r 2 ', consider the case where one wave is incident when the noise level is extremely small. When the arrival angle θ=θ 1 , the second arrival wave number estimation threshold λ th21 ) is obtained from equation (9) to equation (10).
λ th21 )=r 21 )・λ 1 ≦ r 2 ′・λ 1 (10)
Therefore, the threshold value can be set smaller when r 21 ) is used than r 2 ′, so the accuracy of estimating the number of arriving waves when two or more waves are incident near the arrival angle θ 1 becomes higher.

また、図8で示すように、ステップS8の後に、第3到来波数推定閾値λth3を演算するステップS8aを設けておくこともできる。すなわち、ステップS8aにおいて、到来波数推定閾値演算部107は、第3到来波数推定閾値λth3を式(11)のように演算する。 Further, as shown in FIG. 8, a step S8a for calculating a third arriving wave number estimation threshold λ th3 may be provided after step S8. That is, in step S8a, the arriving wave number estimation threshold calculation unit 107 calculates the third arriving wave number estimation threshold λ th3 as shown in equation (11).

λth3 =p・λth1 +q・λth2
=p・c・P +q・r′・λ ・・・(11)
ここで、p、qは任意定数、c=c+cである。式(11)の第3到来波数推定閾値λth3は、第1到来波数推定閾値λth1及び第2到来波数推定閾値λth2にそれぞれ重み付けをつけて加算した閾値である。この第3到来波数推定閾値λth3を使用することで、受信信号強度によらず、固有値を信号固有値と雑音固有値に判別することができる。
λ th3 =p・λ th1 +q・λ th2
= p・c・P n +q・r 2 ′・λ 1 ...(11)
Here, p and q are arbitrary constants, and c=c m +c s . The third arriving wave number estimation threshold λ th3 in Equation (11) is a threshold obtained by weighting and adding the first arriving wave number estimation threshold λ th1 and the second arriving wave number estimation threshold λ th2 , respectively. By using this third arriving wave number estimation threshold λ th3 , it is possible to distinguish the eigenvalue into a signal eigenvalue and a noise eigenvalue, regardless of the received signal strength.

例として、1波入射時の電波環境を考える。p、q=1としたときの受信信号強度に対する固有値及び各到来波数推定閾値の関係を図9に示す。このとき、式(11)は式(12)のように表される。
λth3 =c・P +r′・λ >λ ・・・(12)
As an example, consider the radio wave environment when one wave is incident. FIG. 9 shows the relationship between the eigenvalue and each arriving wave number estimation threshold with respect to the received signal strength when p, q=1. At this time, equation (11) is expressed as equation (12).
λ th3 =c・P n +r 2 ′・λ 12 ...(12)

図9は、雑音レベルが一定の例であり、第1到来波数推定閾値λth1は受信信号強度に依存せず、一定の値である。第2到来波数推定閾値λth2及は受信信号強度の増加とともに大きな値を取る。受信信号強度が小さい領域ではλth1>λで、λth3>λであり、1波入射時に到来波数は1と推定することが可能であるが、λth2<λで第2到来波数推定閾値λth2では、到来波数は1と推定されない。受信信号強度が大きい領域では、λth2>λで、λth3>λであり、1波入射時に到来波数は1と推定することが可能であるが、λth1<λで第1到来波数推定閾値λth1では、到来波数は1と推定されない。また、信号強度によっては、第1到来波数推定閾値λth1及び第2到来波数推定閾値λth2とも到来波数は1でなく、2と推定する領域が生じる。 FIG. 9 is an example in which the noise level is constant, and the first arriving wave number estimation threshold λ th1 is a constant value independent of the received signal strength. The second arriving wave number estimation threshold λ th2 takes a larger value as the received signal strength increases. In the region where the received signal strength is small, λ th1 > λ 2 and λ th3 > λ 2 , and it is possible to estimate the number of arriving waves as 1 when one wave is incident, but if λ th2 < λ 2 , the number of arriving waves is 1. At the estimated threshold value λ th2 , the number of arriving waves is not estimated to be 1. In a region where the received signal strength is high, λ th2 > λ 2 and λ th3 > λ 2 , and it is possible to estimate the number of arriving waves as 1 when one wave is incident, but when λ th1 < λ 2 , the number of arriving waves is 1. At the wave number estimation threshold λ th1 , the number of arriving waves is not estimated to be 1. Furthermore, depending on the signal strength, there may be a region in which the number of arriving waves is estimated to be 2 instead of 1 for both the first arriving wave number estimation threshold λ th1 and the second arriving wave number estimation threshold λ th2 .

そこで、第3到来波数推定閾値λth3を用いることで、受信信号強度によらず大きさがλ以下となる固有値は雑音固有値として判別することができ、図9で示すように1波入射時に到来波数は1と推定することが可能である(ステップS9)。 Therefore, by using the third arrival wave number estimation threshold λ th3 , eigenvalues whose magnitude is λ 2 or less regardless of the received signal strength can be determined as noise eigenvalues, and as shown in FIG. 9, when one wave is incident, The number of arriving waves can be estimated to be 1 (step S9).

しかしながら、実際には演算処理装置11の演算精度及びPの算出誤差等が内在するため、所望の精度となるようにp、qを調整する必要がある。 However, in reality, there is a calculation accuracy of the arithmetic processing unit 11 and a calculation error of Pn , so it is necessary to adjust p and q to obtain the desired accuracy.

図8のフローチャートにおいても、図4のフローチャートと同様に、ステップS10において、到来角度演算部109は、推定された到来信号の数及び、推定された到来信号に対応する信号固有値を用いてMUSICあるいはESPRIT等の高分解能到来方向推定法を実施することにより、到来角度を高精度に算出することができる。 Also in the flowchart of FIG. 8, similarly to the flowchart of FIG. By implementing a high resolution direction of arrival estimation method such as ESPRIT, the angle of arrival can be calculated with high accuracy.

また、図4及び図8のステップS5において、異なる形状すなわち、受信アンテナの組み合わせ方の異なるものからなるサブアレイを複数群用いてもよい。上述のサブアレイ_1~3を第1の群として、他の群のサブアレイに対し、同様に式(11)の到来波数推定閾値を演算し、推定された到来波数を比較した後、複数の群の結果を統合して最終到来波数を確定する。このように異なる形状のサブアレイを複数群用いることにより、サブアレイ形状に依存する到来波数推定精度の低下を防ぐことができる。 Further, in step S5 of FIGS. 4 and 8, a plurality of groups of subarrays having different shapes, that is, different combinations of receiving antennas may be used. With the above-mentioned subarrays_1 to 3 as the first group, the threshold for estimating the number of arriving waves in equation (11) is calculated in the same way for the subarrays in other groups, and after comparing the estimated numbers of arriving waves, The results are integrated to determine the final number of arriving waves. By using a plurality of groups of subarrays having different shapes in this way, it is possible to prevent the accuracy of estimating the number of arriving waves from decreasing depending on the shape of the subarray.

図4及び図8のフローチャートにおいて、上述したステップS5からステップS10の処理を、目標物体数分だけ繰り返し、全ての目標物体の距離、相対速度、到来角度を目標物体情報として算出し、外部装置(図示せず)に算出した目標物体情報を出力する。 In the flowcharts of FIGS. 4 and 8, the processes from step S5 to step S10 described above are repeated for the number of target objects, the distances, relative velocities, and angles of arrival of all target objects are calculated as target object information, and the external device ( (not shown), the calculated target object information is output.

以上のように、本実施の形態1に係るレーダ装置100によれば、受信信号強度が小さい場合に良好な到来波数推定精度を維持しつつ、受信信号強度が大きくなった場合にも良好に到来波数を推定することができる。すなわち、本実施の形態1に係るレーダ装置100は、物体から反射される反射波を複数の受信アンテナ6で受信し、受信した信号を周波数解析し、得られた複素スペクトラムに基づいて、雑音レベルを演算するとともにピーク値を判定する。そして、同一物体から異なる時刻に受信された信号から周波数解析によって得られた複素スペクトラムであって、ピーク値と判定された複素スペクトラムに対応する相関行列をそれぞれ演算するとともにそれらを平均化し、平均化された相関行列に対し、固有値及び固有値ベクトルの解析を行う。あるいは、複数の受信アンテナ6のうち異なる複数の受信アンテナ6の組で形成された複数のサブアレイに対し、サブアレイ毎にそのサブアレイに属する複素スペクトラムに対し相関行列を演算し、サブアレイ毎の相関行列を平均化することで得られるサブアレイ群の相関行列に対し、固有値及び固有値ベクトルの解析を行う。ここで、サブアレイを形成する受信アンテナは仮想受信アンテナであってもよい。本実施の形態1に係るレーダ装置100では、到来波数推定閾値演算部107において、信号強度が小さい領域に対し、雑音レベルを考慮して第1到来波数推定閾値λth1を設定する。さらに、到来波数推定閾値演算部107において、信号強度が大きい領域に対し、受信アンテナの振幅誤差、位相誤差を考慮して固有値に基いた第2到来波数推定閾値λth2を設定する。信号強度に応じて、第1到来波数推定閾値λth1及び第2到来波数推定閾値λth2のいずれかを用いて固有値から信号固有値を判別し、その信号固有値を用い、到来角度演算部109は、到来信号の数を推定し、到来信号に対応する信号固有値を用いて、高分解能到来方向推定法を実施することにより、到来角度を高精度に算出することが可能となる。 As described above, according to the radar apparatus 100 according to the first embodiment, good incoming wave number estimation accuracy is maintained when the received signal strength is small, and good incoming wave number estimation accuracy is maintained even when the received signal strength becomes large. The wave number can be estimated. That is, the radar device 100 according to the first embodiment receives reflected waves reflected from an object using a plurality of receiving antennas 6, performs frequency analysis on the received signals, and determines the noise level based on the obtained complex spectrum. is calculated and the peak value is determined. Then, correlation matrices corresponding to the complex spectra that are determined to be peak values, which are complex spectra obtained by frequency analysis from signals received from the same object at different times, are calculated and averaged. The eigenvalues and eigenvalue vectors of the correlation matrix are analyzed. Alternatively, for a plurality of subarrays formed by different sets of receiving antennas 6 among the plurality of receiving antennas 6, a correlation matrix is calculated for the complex spectrum belonging to the subarray for each subarray, and the correlation matrix for each subarray is calculated. Eigenvalues and eigenvalue vectors are analyzed for the correlation matrix of the subarray group obtained by averaging. Here, the receiving antennas forming the subarray may be virtual receiving antennas. In the radar device 100 according to the first embodiment, the arriving wave number estimation threshold calculation unit 107 sets the first arriving wave number estimation threshold λ th1 for a region where the signal strength is low, taking into consideration the noise level. Furthermore, the arriving wave number estimation threshold calculation unit 107 sets a second arriving wave number estimation threshold λ th2 based on the eigenvalue, taking into account the amplitude error and phase error of the receiving antenna, for a region where the signal strength is high. Depending on the signal strength, the signal eigenvalue is determined from the eigenvalues using either the first arrival wave number estimation threshold λ th1 or the second arrival wave number estimation threshold λ th2 , and using the signal eigenvalue, the arrival angle calculation unit 109 calculates the following: By estimating the number of arriving signals and performing a high-resolution direction-of-arrival estimation method using signal eigenvalues corresponding to the arriving signals, it becomes possible to calculate the angle of arrival with high precision.

また、本実施の形態1に係るレーダ装置100の到来波数推定閾値演算部107は、第1到来波数推定閾値λth1及び第2到来波数推定閾値λth2に対し、それぞれ重み付けして加算した第3到来波数推定閾値λth3を演算する。この第3到来波数推定閾値λth3を用いて、固有値から信号固有値を判別するようにしたので、信号強度によらず固有値から信号固有値を高精度に判別することが可能となる。その信号固有値を用いて、到来角度演算部109で、高分解能到来方向推定法を実施することにより、到来角度をさらに高精度を算出することが可能となる。 In addition, the arriving wave number estimation threshold calculation unit 107 of the radar device 100 according to the first embodiment also calculates a third value that is weighted and added to the first arriving wave number estimation threshold λ th1 and the second arriving wave number estimation threshold λ th2 . A threshold value λ th3 for estimating the number of arriving waves is calculated. Since the signal eigenvalue is determined from the eigenvalue using this third arrival wave number estimation threshold λ th3 , it becomes possible to determine the signal eigenvalue from the eigenvalue with high accuracy regardless of the signal strength. Using the signal eigenvalues, the angle of arrival calculating section 109 performs a high-resolution direction of arrival estimation method, thereby making it possible to calculate the angle of arrival with even higher accuracy.

雑音レベルに基づく閾値の設定のみでは、信号強度が大きい領域では、到来信号の数を誤推定し、その結果到来角度の算出の精度が低下する虞があったが、本実施の形態では、第2到来波数推定閾値λth2を設定するので、信号強度が増加しても到来信号の数を高精度に推定可能となる。 If only the threshold value is set based on the noise level, there is a risk that the number of arriving signals will be erroneously estimated in areas where the signal strength is large, resulting in a decrease in the accuracy of calculating the arrival angle. Since the threshold value λ th2 for estimating the number of arriving waves is set, the number of arriving signals can be estimated with high accuracy even if the signal strength increases.

本願は、例示的な実施の形態が記載されているが、実施の形態に記載された様々な特徴、態様、及び機能は特定の実施の形態の適用に限られるのではなく、単独で、または様々な組み合わせで実施の形態に適用可能である。
従って、例示されていない無数の変形例が、本願明細書に開示される技術の範囲内において想定される。例えば、少なくとも1つの構成要素を変形する場合、追加する場合または省略する場合が含まれるものとする。
Although this application describes exemplary embodiments, the various features, aspects, and functions described in the embodiments are not limited to the application of particular embodiments, and may be used alone or It is applicable to the embodiments in various combinations.
Accordingly, countless variations not illustrated are envisioned within the scope of the technology disclosed herein. For example, this includes cases in which at least one component is modified, added, or omitted.

以下、本開示の諸態様を付記としてまとめて記載する。 Hereinafter, various aspects of the present disclosure will be collectively described as supplementary notes.

(付記1)
送信信号を送信する送信アンテナと、
単数又は複数の物体で反射された反射信号を受信する複数の受信アンテナと、
複数の前記受信アンテナで受信した前記反射信号を周波数解析し、前記物体のそれぞれについて、自装置との角度である到来角度を算出する制御部と、を備えたレーダ装置であって、
前記制御部は、
前記周波数解析により得られた複素スペクトラムの振幅値から雑音レベルを演算する雑音レベル演算部と、
前記周波数解析により得られた前記複素スペクトラムの振幅値及び前記雑音レベル演算部で演算された前記雑音レベルに基づいてピーク値を判定するピーク値判定部と、
前記ピーク値と判定された複素スペクトラムに対する相関行列を演算する相関行列演算部と、
前記相関行列の固有値を演算する固有値演算部と、
前記固有値を雑音に起因する雑音固有値と前記物体に対応する信号固有値とに判別するための閾値を演算する到来波数推定閾値演算部と、
前記到来波数推定閾値演算部で演算された閾値に基づいて、判別された前記信号固有値の数から到来信号の数を推定し、前記到来信号と推定された前記信号固有値に基づいて前記到来角度を演算する到来角度演算部と、を有し、
前記到来波数推定閾値演算部は、
前記雑音レベル演算部で演算された雑音レベルに基づいて第1到来波数推定閾値を演算するとともに、前記固有値演算部で演算された前記固有値に基づいて第2到来波数推定閾値を演算し、
前記到来角度演算部は、前記第1到来波数推定閾値及び前記第2到来波数推定閾値を用いて判別された前記信号固有値の数から前記到来信号の数を推定する、レーダ装置。
(付記2)
前記到来波数推定閾値演算部は、前記雑音レベルに予め設定された第1の値を乗算することで前記第1到来波数推定閾値を演算する付記1に記載のレーダ装置。
(付記3)
前記到来波数推定閾値演算部は、前記固有値演算部で演算された前記固有値のうち最大の固有値に予め設定された第2の値を乗算することで前記第2到来波数推定閾値を演算する付記1または2に記載のレーダ装置。
(付記4)
前記反射信号の到来角度を予め推定し、前記第2の値は、予め推定された前記反射信号の到来角度に応じた値を用いる、付記3に記載のレーダ装置。
(付記5)
前記第1到来波数推定閾値は前記反射信号の信号強度の小さい領域で前記信号固有値の判別に用いられ、前記第2到来波数推定閾値は前記第1到来波数推定閾値が用いられるよりも前記反射信号の信号強度の大きい領域で前記信号固有値の判別に用いられる、付記1から4のいずれか1つの付記に記載のレーダ装置。
(付記6)
前記到来波数推定閾値演算部は、前記第1到来波数推定閾値及び前記第2到来波数推定閾値をそれぞれ重み付けして加算した第3到来波数推定閾値を演算し、
前記到来角度演算部は、前記第3到来波数推定閾値を用いて判別された前記信号固有値の数から到来信号数を推定し、前記到来角度を演算する、付記1から5のいずれか1つの付記に記載のレーダ装置。
(付記7)
前記相関行列演算部は、複数の前記受信アンテナのうち異なる複数の受信アンテナの組で形成された複数のサブアレイ毎に前記サブアレイに属し前記ピーク値と判定された複素スペクトラムに対する相関行列を演算し、得られた複数の前記サブアレイに対応する相関行列を平均化し、
前記固有値演算部は、平均化された前記相関行列の固有値を演算する、付記1から6のいずれか1つの付記に記載のレーダ装置。
(付記8)
複数の前記受信アンテナのうち異なる複数の受信アンテナの組で形成された複数のサブアレイを第1群アレイとし、前記第1群アレイとは複数の前記受信アンテナの組み合わせ方の異なる複数の受信アンテナの組で形成された複数のサブアレイを第2群アレイとし、
前記到来角度演算部は、
前記各群アレイに対し、それぞれ前記到来信号の数を推定し、推定された複数の前記到来信号の数をもとに最終の到来角度を決定する、付記7に記載のレーダ装置。
(付記9)
前記相関行列演算部は、同一の前記物体から反射され異なる時刻に受信した複数の前記反射信号を周波数解析することにより得られ、前記ピーク値と判定された複数の前記複素スペクトラムに対する相関行列をそれぞれ演算し、得られた前記相関行列を平均化し、
前記固有値演算部は、平均化された前記相関行列の固有値を演算する、付記1から6のいずれか1つの付記に記載のレーダ装置。
(Additional note 1)
a transmitting antenna for transmitting a transmitting signal;
a plurality of receiving antennas that receive reflected signals reflected by one or more objects;
A radar device comprising: a control unit that frequency-analyzes the reflected signals received by a plurality of the receiving antennas and calculates an angle of arrival, which is an angle with respect to the own device, for each of the objects,
The control unit includes:
a noise level calculation unit that calculates a noise level from the amplitude value of the complex spectrum obtained by the frequency analysis;
a peak value determination unit that determines a peak value based on the amplitude value of the complex spectrum obtained by the frequency analysis and the noise level calculated by the noise level calculation unit;
a correlation matrix calculation unit that calculates a correlation matrix for the complex spectrum determined to be the peak value;
an eigenvalue calculation unit that calculates an eigenvalue of the correlation matrix;
an arriving wave number estimation threshold calculation unit that calculates a threshold for discriminating the eigenvalue into a noise eigenvalue caused by noise and a signal eigenvalue corresponding to the object;
The number of arriving signals is estimated from the number of determined signal eigenvalues based on the threshold calculated by the arriving wave number estimation threshold calculation unit, and the arrival angle is determined based on the incoming signal and the estimated signal eigenvalue. an arrival angle calculation unit that calculates,
The arriving wave number estimation threshold calculation unit includes:
calculating a first arrival wave number estimation threshold based on the noise level calculated by the noise level calculation unit, and calculating a second arrival wave number estimation threshold based on the eigenvalue calculated by the eigenvalue calculation unit;
The arrival angle calculation unit is a radar device that estimates the number of arriving signals from the number of signal eigenvalues determined using the first arriving wave number estimation threshold and the second arrival wave number estimation threshold.
(Additional note 2)
The radar device according to supplementary note 1, wherein the arriving wave number estimation threshold calculation unit calculates the first arriving wave number estimation threshold by multiplying the noise level by a first value set in advance.
(Additional note 3)
Supplementary note 1, wherein the arriving wave number estimation threshold calculation unit calculates the second arrival wave number estimation threshold by multiplying the largest eigenvalue among the eigenvalues calculated by the eigenvalue calculation unit by a second value set in advance. or the radar device described in 2.
(Additional note 4)
The radar device according to supplementary note 3, wherein the angle of arrival of the reflected signal is estimated in advance, and the second value uses a value corresponding to the previously estimated angle of arrival of the reflected signal.
(Appendix 5)
The first arrival wave number estimation threshold is used to determine the signal eigenvalue in a region where the signal strength of the reflected signal is small, and the second arrival wave number estimation threshold is used to determine the reflected signal in a region where the signal intensity of the reflected signal is small. The radar device according to any one of appendices 1 to 4, which is used to determine the signal eigenvalue in a region where the signal strength is large.
(Appendix 6)
The arriving wave number estimation threshold calculation unit calculates a third arriving wave number estimation threshold obtained by weighting and adding the first arriving wave number estimation threshold and the second arriving wave number estimation threshold, respectively,
The arrival angle calculation unit estimates the number of arriving signals from the number of signal eigenvalues determined using the third arrival wave number estimation threshold, and calculates the arrival angle. The radar device described in .
(Appendix 7)
The correlation matrix calculation unit calculates a correlation matrix for a complex spectrum that belongs to the subarray and is determined to have the peak value, for each of a plurality of subarrays formed by different sets of reception antennas among the plurality of reception antennas, averaging the obtained correlation matrices corresponding to the plurality of subarrays;
The radar device according to any one of Supplementary Notes 1 to 6, wherein the eigenvalue calculation unit calculates the eigenvalue of the averaged correlation matrix.
(Appendix 8)
A plurality of sub-arrays formed by different sets of receiving antennas among the plurality of receiving antennas are referred to as a first group array, and the first group array is defined as a plurality of receiving antennas in which the plurality of receiving antennas are combined in different ways. A plurality of subarrays formed in a group are set as a second group array,
The arrival angle calculation unit is
The radar device according to supplementary note 7, wherein the number of the arriving signals is estimated for each group array, and the final angle of arrival is determined based on the estimated number of the plurality of arriving signals.
(Appendix 9)
The correlation matrix calculation unit calculates a correlation matrix for each of the plurality of complex spectra determined to have the peak value, which is obtained by frequency-analyzing the plurality of reflected signals reflected from the same object and received at different times. calculate and average the obtained correlation matrix,
The radar device according to any one of Supplementary Notes 1 to 6, wherein the eigenvalue calculation unit calculates the eigenvalue of the averaged correlation matrix.

1:信号処理部、 2:制御電圧発生器、 3:VCO、 4:分配器、 5:送信アンテナ、 6、6_1~6_6:受信アンテナ、 7、7_1~7_6:ミキサ、 8、8_1~8_6:A/Dコンバータ、 9、9_1~9_6:FFT演算部、 10:制御部、 11:演算処理装置、 12:記憶装置、 13:通信回路、 14:D/Aコンバータ、 100:レーダ装置、 101:振幅演算部、 102:雑音レベル演算部、 103:ピーク値判定部、 104:距離及び相対速度演算部、 105:相関行列演算部、 106:固有値演算部、 107:到来波数推定閾値演算部、 108:固有値判定部、 109:到来角度演算部。 1: Signal processing unit, 2: Control voltage generator, 3: VCO, 4: Distributor, 5: Transmitting antenna, 6, 6_1 to 6_6: Receiving antenna, 7, 7_1 to 7_6: Mixer, 8, 8_1 to 8_6: A/D converter, 9, 9_1 to 9_6: FFT calculation unit, 10: Control unit, 11: Processing unit, 12: Storage device, 13: Communication circuit, 14: D/A converter, 100: Radar device, 101: Amplitude calculation unit, 102: Noise level calculation unit, 103: Peak value determination unit, 104: Distance and relative velocity calculation unit, 105: Correlation matrix calculation unit, 106: Eigenvalue calculation unit, 107: Arrival wave number estimation threshold calculation unit, 108 : Eigenvalue determination unit, 109: Arrival angle calculation unit.

本願に開示されるレーダ装置は、
送信信号を送信する送信アンテナと、
単数又は複数の物体で反射された反射信号を受信する複数の受信アンテナと、
複数の前記受信アンテナで受信した前記反射信号を周波数解析し、前記物体のそれぞれについて、自装置との角度である到来角度を算出する制御部と、を備えたレーダ装置であって、
前記制御部は、
前記周波数解析により得られた複素スペクトラムの振幅値から雑音レベルを演算する雑音レベル演算部と、
前記周波数解析により得られた前記複素スペクトラムの振幅値及び前記雑音レベル演算部で演算された前記雑音レベルに基づいてピーク値を判定するピーク値判定部と、
前記ピーク値と判定された複素スペクトラムに対する相関行列を演算する相関行列演算部と、
前記相関行列の固有値を演算する固有値演算部と、
前記固有値を雑音に起因する雑音固有値と前記物体に対応する信号固有値とに判別するための閾値を演算する到来波数推定閾値演算部と、
前記到来波数推定閾値演算部で演算された閾値に基づいて、判別された前記信号固有値の数から到来信号の数を推定し、前記到来信号と推定された前記信号固有値に基づいて前記到来角度を演算する到来角度演算部と、を有し、
前記到来波数推定閾値演算部は、
前記雑音レベル演算部で演算された雑音レベルに基づいて第1到来波数推定閾値を演算するとともに、前記固有値演算部で演算された前記固有値に基づいて第2到来波数推定閾値を演算し、
前記到来角度演算部は、前記反射信号の信号強度の小さい領域では前記第1到来波数推定閾値に基づいて前記信号固有値を判別し、前記第1到来波数推定閾値が用いられるよりも前記反射信号の信号強度の大きい領域では前記第2到来波数推定閾値に基づいて前記信号固有値を判別し、判別された前記信号固有値の数から前記到来信号の数を推定する、ように構成されている。
The radar device disclosed in this application includes:
a transmitting antenna for transmitting a transmitting signal;
a plurality of receiving antennas that receive reflected signals reflected by one or more objects;
A radar device comprising: a control unit that frequency-analyzes the reflected signals received by a plurality of the receiving antennas and calculates an angle of arrival, which is an angle with respect to the own device, for each of the objects,
The control unit includes:
a noise level calculation unit that calculates a noise level from the amplitude value of the complex spectrum obtained by the frequency analysis;
a peak value determination unit that determines a peak value based on the amplitude value of the complex spectrum obtained by the frequency analysis and the noise level calculated by the noise level calculation unit;
a correlation matrix calculation unit that calculates a correlation matrix for the complex spectrum determined to be the peak value;
an eigenvalue calculation unit that calculates an eigenvalue of the correlation matrix;
an arriving wave number estimation threshold calculation unit that calculates a threshold for discriminating the eigenvalue into a noise eigenvalue caused by noise and a signal eigenvalue corresponding to the object;
The number of arriving signals is estimated from the number of determined signal eigenvalues based on the threshold calculated by the arriving wave number estimation threshold calculation unit, and the arrival angle is determined based on the incoming signal and the estimated signal eigenvalue. an arrival angle calculation unit that calculates,
The arriving wave number estimation threshold calculation unit includes:
calculating a first arrival wave number estimation threshold based on the noise level calculated by the noise level calculation unit, and calculating a second arrival wave number estimation threshold based on the eigenvalue calculated by the eigenvalue calculation unit;
The arrival angle calculation unit determines the signal eigenvalue based on the first arrival wave number estimation threshold in a region where the signal strength of the reflected signal is small, and the arrival angle calculation unit determines the signal eigenvalue based on the first arrival wave number estimation threshold, and In a region where the signal strength is high, the signal eigenvalues are determined based on the second arriving wave number estimation threshold, and the number of the arriving signals is estimated from the number of determined signal eigenvalues.

図3は、レーダ装置100の制御部10のハードウエアの構成例を示す図である。図に示すように、制御部10は処理回路として、DSP(Digital Signal Processor)等の演算処理装置11(コンピュータ)、演算処理装置11とデータのやり取りをする記憶装置12、演算処理装置11にビート信号を入力するA/Dコンバータ8、演算処理装置11から制御電圧発生器2に周波数変調信号の指令値を出力するD/Aコンバータ14、及び通信回路13等を備えている。 FIG. 3 is a diagram showing an example of the hardware configuration of the control unit 10 of the radar device 100. As shown in FIG. 3 , the control unit 10 serves as a processing circuit, and includes an arithmetic processing device 11 (computer) such as a DSP (Digital Signal Processor), a storage device 12 that exchanges data with the arithmetic processing device 11, and an arithmetic processing device 11. It includes an A/D converter 8 that inputs a beat signal, a D/A converter 14 that outputs a command value of a frequency modulation signal from an arithmetic processing unit 11 to a control voltage generator 2, a communication circuit 13, and the like.

Claims (9)

送信信号を送信する送信アンテナと、
単数又は複数の物体で反射された反射信号を受信する複数の受信アンテナと、
複数の前記受信アンテナで受信した前記反射信号を周波数解析し、前記物体のそれぞれについて、自装置との角度である到来角度を算出する制御部と、を備えたレーダ装置であって、
前記制御部は、
前記周波数解析により得られた複素スペクトラムの振幅値から雑音レベルを演算する雑音レベル演算部と、
前記周波数解析により得られた前記複素スペクトラムの振幅値及び前記雑音レベル演算部で演算された前記雑音レベルに基づいてピーク値を判定するピーク値判定部と、
前記ピーク値と判定された複素スペクトラムに対する相関行列を演算する相関行列演算部と、
前記相関行列の固有値を演算する固有値演算部と、
前記固有値を雑音に起因する雑音固有値と前記物体に対応する信号固有値とに判別するための閾値を演算する到来波数推定閾値演算部と、
前記到来波数推定閾値演算部で演算された閾値に基づいて、判別された前記信号固有値の数から到来信号の数を推定し、前記到来信号と推定された前記信号固有値に基づいて前記到来角度を演算する到来角度演算部と、を有し、
前記到来波数推定閾値演算部は、
前記雑音レベル演算部で演算された雑音レベルに基づいて第1到来波数推定閾値を演算するとともに、前記固有値演算部で演算された前記固有値に基づいて第2到来波数推定閾値を演算し、
前記到来角度演算部は、前記第1到来波数推定閾値及び前記第2到来波数推定閾値を用いて判別された前記信号固有値の数から前記到来信号の数を推定する、レーダ装置。
a transmitting antenna for transmitting a transmitting signal;
a plurality of receiving antennas that receive reflected signals reflected by one or more objects;
A radar device comprising: a control unit that frequency-analyzes the reflected signals received by a plurality of the receiving antennas and calculates an angle of arrival, which is an angle with respect to the own device, for each of the objects,
The control unit includes:
a noise level calculation unit that calculates a noise level from the amplitude value of the complex spectrum obtained by the frequency analysis;
a peak value determination unit that determines a peak value based on the amplitude value of the complex spectrum obtained by the frequency analysis and the noise level calculated by the noise level calculation unit;
a correlation matrix calculation unit that calculates a correlation matrix for the complex spectrum determined to be the peak value;
an eigenvalue calculation unit that calculates an eigenvalue of the correlation matrix;
an arriving wave number estimation threshold calculation unit that calculates a threshold for discriminating the eigenvalue into a noise eigenvalue caused by noise and a signal eigenvalue corresponding to the object;
The number of arriving signals is estimated from the number of determined signal eigenvalues based on the threshold calculated by the arriving wave number estimation threshold calculation unit, and the arrival angle is determined based on the incoming signal and the estimated signal eigenvalue. an arrival angle calculation unit that calculates,
The arriving wave number estimation threshold calculation unit includes:
calculating a first arrival wave number estimation threshold based on the noise level calculated by the noise level calculation unit, and calculating a second arrival wave number estimation threshold based on the eigenvalue calculated by the eigenvalue calculation unit;
The arrival angle calculation unit is a radar device that estimates the number of arriving signals from the number of signal eigenvalues determined using the first arriving wave number estimation threshold and the second arrival wave number estimation threshold.
前記到来波数推定閾値演算部は、前記雑音レベルに予め設定された第1の値を乗算することで前記第1到来波数推定閾値を演算する請求項1に記載のレーダ装置。 The radar device according to claim 1, wherein the arriving wave number estimation threshold calculation unit calculates the first arriving wave number estimation threshold by multiplying the noise level by a first value set in advance. 前記到来波数推定閾値演算部は、前記固有値演算部で演算された前記固有値のうち最大の固有値に予め設定された第2の値を乗算することで前記第2到来波数推定閾値を演算する請求項1に記載のレーダ装置。 2. The arriving wave number estimation threshold calculation unit calculates the second arrival wave number estimation threshold by multiplying the largest eigenvalue among the eigenvalues calculated by the eigenvalue calculation unit by a preset second value. 1. The radar device according to 1. 前記反射信号の到来角度を予め推定し、前記第2の値は、予め推定された前記反射信号の到来角度に応じた値を用いる、請求項3に記載のレーダ装置。 The radar device according to claim 3, wherein the angle of arrival of the reflected signal is estimated in advance, and the second value uses a value corresponding to the previously estimated angle of arrival of the reflected signal. 前記第1到来波数推定閾値は前記反射信号の信号強度の小さい領域で前記信号固有値の判別に用いられ、前記第2到来波数推定閾値は前記第1到来波数推定閾値が用いられるよりも前記反射信号の信号強度の大きい領域で前記信号固有値の判別に用いられる、請求項1から4のいずれか1項に記載のレーダ装置。 The first arrival wave number estimation threshold is used to determine the signal eigenvalue in a region where the signal strength of the reflected signal is small, and the second arrival wave number estimation threshold is used to determine the reflected signal in a region where the signal intensity of the reflected signal is small. The radar device according to any one of claims 1 to 4, wherein the radar device is used to determine the signal eigenvalue in a region where the signal intensity is large. 前記到来波数推定閾値演算部は、前記第1到来波数推定閾値及び前記第2到来波数推定閾値をそれぞれ重み付けして加算した第3到来波数推定閾値を演算し、
前記到来角度演算部は、前記第3到来波数推定閾値を用いて判別された前記信号固有値の数から到来信号数を推定し、前記到来角度を演算する、請求項1から3のいずれか1項に記載のレーダ装置。
The arriving wave number estimation threshold calculation unit calculates a third arriving wave number estimation threshold obtained by weighting and adding the first arriving wave number estimation threshold and the second arriving wave number estimation threshold, respectively,
4. The angle of arrival calculation unit estimates the number of arriving signals from the number of signal eigenvalues determined using the third arrival wave number estimation threshold, and calculates the angle of arrival. The radar device described in .
前記相関行列演算部は、複数の前記受信アンテナのうち異なる複数の受信アンテナの組で形成された複数のサブアレイ毎に前記サブアレイに属し前記ピーク値と判定された複素スペクトラムに対する相関行列を演算し、得られた複数の前記サブアレイに対応する相関行列を平均化し、
前記固有値演算部は、平均化された前記相関行列の固有値を演算する、請求項1から3のいずれか1項に記載のレーダ装置。
The correlation matrix calculation unit calculates a correlation matrix for a complex spectrum that belongs to the subarray and is determined to have the peak value, for each of a plurality of subarrays formed by different sets of reception antennas among the plurality of reception antennas, averaging the obtained correlation matrices corresponding to the plurality of subarrays;
The radar device according to any one of claims 1 to 3, wherein the eigenvalue calculation unit calculates an eigenvalue of the averaged correlation matrix.
複数の前記受信アンテナのうち異なる複数の受信アンテナの組で形成された複数のサブアレイを第1群アレイとし、前記第1群アレイとは複数の前記受信アンテナの組み合わせ方の異なる複数の受信アンテナの組で形成された複数のサブアレイを第2群アレイとし、
前記到来角度演算部は、
前記各群アレイに対し、それぞれ前記到来信号の数を推定し、推定された複数の前記到来信号の数をもとに最終の到来角度を決定する、請求項7に記載のレーダ装置。
A plurality of sub-arrays formed by different sets of receiving antennas among the plurality of receiving antennas are referred to as a first group array, and the first group array is defined as a plurality of receiving antennas in which the plurality of receiving antennas are combined in different ways. A plurality of subarrays formed in a group are set as a second group array,
The arrival angle calculation unit is
8. The radar device according to claim 7, wherein the number of the arriving signals is estimated for each of the array groups, and the final angle of arrival is determined based on the estimated number of the plurality of arriving signals.
前記相関行列演算部は、同一の前記物体から反射され異なる時刻に受信した複数の前記反射信号を周波数解析することにより得られ、前記ピーク値と判定された複数の前記複素スペクトラムに対する相関行列をそれぞれ演算し、得られた前記相関行列を平均化し、
前記固有値演算部は、平均化された前記相関行列の固有値を演算する、請求項1から3のいずれか1項に記載のレーダ装置。
The correlation matrix calculation unit calculates a correlation matrix for each of the plurality of complex spectra determined to have the peak value, which is obtained by frequency-analyzing the plurality of reflected signals reflected from the same object and received at different times. calculate and average the obtained correlation matrix,
The radar device according to any one of claims 1 to 3, wherein the eigenvalue calculation unit calculates an eigenvalue of the averaged correlation matrix.
JP2022072896A 2022-04-27 2022-04-27 radar equipment Active JP7191262B1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2022072896A JP7191262B1 (en) 2022-04-27 2022-04-27 radar equipment
PCT/JP2022/025052 WO2023210030A1 (en) 2022-04-27 2022-06-23 Radar device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2022072896A JP7191262B1 (en) 2022-04-27 2022-04-27 radar equipment

Publications (2)

Publication Number Publication Date
JP7191262B1 JP7191262B1 (en) 2022-12-16
JP2023162528A true JP2023162528A (en) 2023-11-09

Family

ID=84488902

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022072896A Active JP7191262B1 (en) 2022-04-27 2022-04-27 radar equipment

Country Status (2)

Country Link
JP (1) JP7191262B1 (en)
WO (1) WO2023210030A1 (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004309209A (en) * 2003-04-03 2004-11-04 Mitsubishi Electric Corp Sensor signal processing system
JP2008281424A (en) * 2007-05-10 2008-11-20 Mitsubishi Electric Corp Frequency modulation radar device
JP2013152239A (en) * 2013-03-29 2013-08-08 Honda Elesys Co Ltd Electronic scanning radar device, control method of the same and program
JP2016125981A (en) * 2015-01-08 2016-07-11 古野電気株式会社 Signal processing device, underwater sensing device and radar device
WO2019215790A1 (en) * 2018-05-07 2019-11-14 三菱電機株式会社 Incoming wave number estimation apparatus and incoming wave number incoming direction estimation apparatus
JP2020020697A (en) * 2018-08-01 2020-02-06 株式会社デンソーテン Arrival direction estimating device and arrival direction estimating method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004309209A (en) * 2003-04-03 2004-11-04 Mitsubishi Electric Corp Sensor signal processing system
JP2008281424A (en) * 2007-05-10 2008-11-20 Mitsubishi Electric Corp Frequency modulation radar device
JP2013152239A (en) * 2013-03-29 2013-08-08 Honda Elesys Co Ltd Electronic scanning radar device, control method of the same and program
JP2016125981A (en) * 2015-01-08 2016-07-11 古野電気株式会社 Signal processing device, underwater sensing device and radar device
WO2019215790A1 (en) * 2018-05-07 2019-11-14 三菱電機株式会社 Incoming wave number estimation apparatus and incoming wave number incoming direction estimation apparatus
JP2020020697A (en) * 2018-08-01 2020-02-06 株式会社デンソーテン Arrival direction estimating device and arrival direction estimating method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
BJORKLUND, SVANTE: "Estimating the Number of Signal Sources in Array Antenna Radar by an Empirical Threshold", 2019 IEEE RADAR CONFERENCE (RADARCONF), JPN6022033562, April 2019 (2019-04-01), US, ISSN: 0004850063 *

Also Published As

Publication number Publication date
WO2023210030A1 (en) 2023-11-02
JP7191262B1 (en) 2022-12-16

Similar Documents

Publication Publication Date Title
CN111373282B (en) Radar processing chain for FMCW radar system
CN108885254B (en) Object detection device
JP4679632B2 (en) Radar equipment
JP4368852B2 (en) Radar equipment
JP6489589B2 (en) Radar signal processing device
JP7108929B2 (en) Radar device and target determination method
JP2009025195A (en) Method of estimating number of incoming waves, and radar device
JP4972852B2 (en) Radar equipment
JP4926155B2 (en) Radar equipment
US20160377713A1 (en) Method for estimating reflected wave arrival direction, and program
US10712428B2 (en) Radar device and target detecting method
US11754671B2 (en) Incoming wave count estimation apparatus and incoming wave count incoming direction estimation apparatus
US8116169B2 (en) Active sonar system and active sonar method using noise reduction techniques and advanced signal processing techniques
WO2023210030A1 (en) Radar device
KR20190040637A (en) SYSTEM AND METHOD FOR ESTIMATING RADAR DoA
Chudnikov et al. Doa estimation in radar sensors with colocated antennas
Wu et al. An improved angle estimation algorithm for millimeter-wave radar
KR20190124488A (en) Method of signal subspace based DoA estimation for automotive radar system
WO2021196165A1 (en) Frequency analysis method, device and radar
JP7056212B2 (en) Direction estimation method and equipment
TWI808874B (en) Radar system for vehicle and detecting method
WO2022249552A1 (en) Information processing device and information processing method
US11876582B2 (en) Hybrid beam steering radar
JP2005121420A (en) Radar system
JP7294781B2 (en) Radar device and object discrimination method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220427

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220816

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220929

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20221108

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20221206

R151 Written notification of patent or utility model registration

Ref document number: 7191262

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151