JP2010271337A - Radar device - Google Patents

Radar device Download PDF

Info

Publication number
JP2010271337A
JP2010271337A JP2010203174A JP2010203174A JP2010271337A JP 2010271337 A JP2010271337 A JP 2010271337A JP 2010203174 A JP2010203174 A JP 2010203174A JP 2010203174 A JP2010203174 A JP 2010203174A JP 2010271337 A JP2010271337 A JP 2010271337A
Authority
JP
Japan
Prior art keywords
incoming
incoming waves
signal
eigenvalue
wave
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
JP2010203174A
Other languages
Japanese (ja)
Other versions
JP4780240B2 (en
Inventor
Kazuma Natsume
一馬 夏目
Yasuyuki Miyake
康之 三宅
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.)
Denso Corp
Original Assignee
Denso 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 Denso Corp filed Critical Denso Corp
Priority to JP2010203174A priority Critical patent/JP4780240B2/en
Publication of JP2010271337A publication Critical patent/JP2010271337A/en
Application granted granted Critical
Publication of JP4780240B2 publication Critical patent/JP4780240B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To accurately estimate the number of incoming waves, in a radar system for estimating the number of the incoming waves by MUSIC method, Esprit method or the like. <P>SOLUTION: The distances and relative speeds (or directions) of targets are obtained based on beat signals, and stored as history information. Target positions in the incoming wave estimation are predicted from the stored history information. Prediction is made about whether two or more targets are aligned, from the predicted target positions, and the number of incoming waves is estimated based on the prediction result. Thus, based on the history information of the targets obtained in the last cycle, prediction is made about whether targets, which have not travelled in parallel in the last cycle, travel in parallel with each other in this cycle, or whether targets, which have travelled in parallel with each other in the last cycle, stop travelling in parallel in this cycle, and the number of incoming waves is revised. The number of incoming waves is thereby estimated more accurately. <P>COPYRIGHT: (C)2011,JPO&INPIT

Description

本発明は、周波数変調されたレーダ波を送受信することにより、目標物体と自分自身との位置関係を求めるレーダ装置に関し、例えば、車両などの移動体の衝突防止等に使用されるレーダ装置に適用して好適である。   The present invention relates to a radar apparatus that obtains a positional relationship between a target object and itself by transmitting and receiving a frequency-modulated radar wave, and is applied to, for example, a radar apparatus used for collision prevention of a moving body such as a vehicle. It is preferable.

従来、FMCW方式のレーダ装置が知られている(例えば、特許文献1参照)。このレーダ装置は、周波数変調されたレーダ波を送受信することにより、目標物体との相対距離や相対速度に関する情報を取り出すことで、例えば前方車両と自車両との関係を検出するようになっている。このようなレーダ装置では、周波数変調されたレーダ波を送受信したときに、自車両の前方にどれだけ他車両等の目標物体が存在し、その目標物体がどの方位に存在するかを求める必要がある。その方位を求める手段として、ビームフォーマー法やMUSIC法、Esprit法などが一般に知られている。この中で、MUSIC(Multiple Signal Classification)法やEsprit法を用いる場合、高い方位分解能が得られるという特徴があるが、到来波方向を演算する過程で到来波数を正確に推定する必要がある。   Conventionally, an FMCW radar device is known (see, for example, Patent Document 1). This radar apparatus detects the relationship between the preceding vehicle and the host vehicle, for example, by extracting information on the relative distance and relative speed with the target object by transmitting and receiving a frequency-modulated radar wave. . In such a radar apparatus, when a frequency-modulated radar wave is transmitted / received, it is necessary to determine how much a target object such as another vehicle exists in front of the host vehicle and in which direction the target object exists. is there. The beam former method, the MUSIC method, the Esprit method, etc. are generally known as means for obtaining the orientation. Among them, when using the MUSIC (Multiple Signal Classification) method or the Esprit method, there is a feature that a high azimuth resolution can be obtained, but it is necessary to accurately estimate the number of incoming waves in the process of calculating the direction of the incoming wave.

この到来波数の推定方法としてAIC(Akaike Information Criteria)やMDL(Minimum Description Length)などがある。これらの手法では、複数回のデータ収集を行い、その分散を評価することで到来波数の推定を行っている。   There are AIC (Akaike Information Criteria) and MDL (Minimum Description Length) as methods for estimating the number of incoming waves. In these methods, the number of incoming waves is estimated by collecting data a plurality of times and evaluating the dispersion.

また一方で、一定の閾値を設け固有値が雑音電力相当か信号電力相当かを区別する手法があり、この手法を閾値法と呼ぶ。例えば、特許文献2において、このMUSIC法の場合の閾値法による到来波数の推定手法が示されている。   On the other hand, there is a method of setting a certain threshold value and distinguishing whether the eigenvalue is equivalent to noise power or signal power, and this method is called a threshold method. For example, Patent Document 2 discloses a method for estimating the number of incoming waves by the threshold method in the case of the MUSIC method.

このMUSIC法の場合の閾値法による到来波数の推定について、図13を参照して説明する。   The estimation of the number of incoming waves by the threshold method in the case of the MUSIC method will be described with reference to FIG.

図13は、複数のアンテナが等間隔にアレイ状に並べられたアダプティブアンテナに対して、到来波が入力されるときの様子を示した図である。   FIG. 13 is a diagram illustrating a state where an incoming wave is input to an adaptive antenna in which a plurality of antennas are arranged in an array at equal intervals.

この図に示されるように、等間隔に並べられたK個のアンテナが車両に搭載されている場合に、各アンテナに車両遠方から複数の電波が到来波として入力された場合を想定する。この場合、各アンテナと到来波の発信源となる目標物体との距離が各アンテナの間隔より十分に長いため、各アンテナに到来波が同じ角度で入力されると想定できる。また、各アンテナに入力される到来波の位相は到来波の方位によって変化する。   As shown in this figure, it is assumed that when K antennas arranged at equal intervals are mounted on a vehicle, a plurality of radio waves are input to each antenna as incoming waves from a distance from the vehicle. In this case, since the distance between each antenna and the target object that is the transmission source of the incoming wave is sufficiently longer than the distance between the antennas, it can be assumed that the incoming wave is input to each antenna at the same angle. In addition, the phase of the incoming wave input to each antenna varies depending on the direction of the incoming wave.

このため、各アンテナに入力された到来波の受信データ(ベクトル)がxi(t)、(i=1〜k)で表されたとすると、各到来波を出力電圧Poutに変換したときにその出力電圧Poutが最小となる場合のウェイトWminは次式のように表される。   For this reason, if the received data (vector) of the incoming wave input to each antenna is represented by xi (t), (i = 1 to k), the output is obtained when each incoming wave is converted into the output voltage Pout. The weight Wmin when the voltage Pout is minimized is expressed by the following equation.

Figure 2010271337
ただし、Rxxは到来波の受信データxi(t)の自己相関行列を示している。また、Wはウェイト、WHはWの共役転換行列を示しており、WHW=1である。
Figure 2010271337
Here, Rxx represents an autocorrelation matrix of incoming wave reception data xi (t). W is a weight, W H is a conjugate conversion matrix of W, and W H W = 1.

ここで、Lagrangeの未定係数法により、固有値λを用いて次式が定義できる。   Here, the following equation can be defined using the eigenvalue λ by Lagrange's undetermined coefficient method.

(数2)
RxxW=λW
そして、この数式の両辺にWHを掛け合わせることで、次式が導出される。
(Equation 2)
RxxW = λW
Then, the following equation is derived by multiplying both sides of this equation by WH .

(数3)
HRxxW=λWHW=λ
このとき、λはウェイトWを掛けたときの出力電圧Poutの2倍に相当する。つまり、到来波の方向に指向性のヌル点(零点)が向けられ、到来波の方向の出力電圧Poutが最小値となるようなWが設定できた場合、すべての到来波からの信号を打ち消すことができるため、λは雑音に等しくなる。逆に、ヌル点が向けられない到来波がある場合には、ヌル点が向けられなかった到来波からの信号を打ち消すことができず、λは値を持つ。したがって、自己相関行列Rxxの固有値λを表すと、K本のアンテナそれぞれで得た受信データから求まる固有値λ1〜λKは次のような関係となる。
(Equation 3)
W H RxxW = λW H W = λ
At this time, λ corresponds to twice the output voltage Pout when the weight W is multiplied. That is, when a directional null point (zero point) is directed in the direction of the incoming wave and W can be set so that the output voltage Pout in the direction of the incoming wave becomes the minimum value, the signals from all the incoming waves are canceled. Can be equal to noise. Conversely, when there is an incoming wave to which the null point is not directed, the signal from the incoming wave to which the null point is not directed cannot be canceled, and λ has a value. Therefore, when the eigenvalue λ of the autocorrelation matrix Rxx is expressed, the eigenvalues λ1 to λK obtained from the reception data obtained by the K antennas have the following relationship.

Figure 2010271337
なお、σ2は、熱雑音電力に相当する。
Figure 2010271337
Note that σ 2 corresponds to thermal noise power.

このため、自己相関行列Rxxの固有値を求め、熱雑音電力σ2より大きい固有値の数から到来波数Lを推定することができる。つまり、数4に示されるλLとλL+1との間にスレッショルドを設定し、固有値λがスレッショルドλTHを超えるか否かにより、到来波数Lを推定することが可能となる。また、スレッショルドλTHを超える領域は、すべての到来波からの信号を打ち消すことができた雑音空間、λTHよりも下回る領域は、少なくとも1つの到来波からの信号を打ち消すことができなかった信号空間となる。したがって、すべての到来波からの信号を打ち消せた雑音空間は、到来波数分の信号を打ち消すことができたものであるため、これに基づいて到来波の方位の算出を行うことが可能となる。 For this reason, the eigenvalue of the autocorrelation matrix Rxx can be obtained, and the number of incoming waves L can be estimated from the number of eigenvalues greater than the thermal noise power σ 2 . That is, it is possible to estimate the number of incoming waves L by setting a threshold between λL and λL + 1 shown in Equation 4 and determining whether the eigenvalue λ exceeds the threshold λTH. Further, the region exceeding the threshold λTH is a noise space in which signals from all incoming waves can be canceled, and the region lower than λTH is a signal space in which signals from at least one incoming wave cannot be canceled. Become. Therefore, the noise space in which the signals from all the incoming waves can be canceled is the one in which the signals corresponding to the number of incoming waves can be canceled, and based on this, the direction of the incoming wave can be calculated. .

特開平11−030663号公報Japanese Patent Laid-Open No. 11-030663 特開2000−121716号公報JP 2000-121716 A

菊間信良、アダプティブアンテナ技術、オーム社、P.138−P.139(MUSIC)、P.145−P.150(Esprit)Nobuyoshi Kikuma, Adaptive Antenna Technology, Ohm, P.A. 138-P. 139 (MUSIC), P.I. 145-P. 150 (Esprit)

これらに示した手法はいずれの手法も、データ収集の回数(スナップショット)が多ければ多いほど、また受信信号のSN比(信号−ノイズ比)が大きければ大きいほど、到来波数をより正確に推定できる。つまり、データ数が少ない場合、正確な到来波数推定ができない。   In any of these methods, the more the number of data collections (snapshots) and the larger the received signal SN ratio (signal-noise ratio), the more accurately the number of incoming waves can be estimated. it can. That is, when the number of data is small, accurate estimation of the number of incoming waves cannot be performed.

しかしながら、車載用レーダ装置にMUSIC法やEsprit法を適用する場合のように、ターゲットとレーダの相対関係の動きが速く、また使用できるデバイスの演算能力が低い場合、同じ状況でのデータを数多く収集することは不可能である。また更に、周囲に複数の反射物体が存在する場合は、マルチパスの影響などにより安定したSN比を確保することも困難である。これらのような状況では、正確な到来波数推定を行うことは困難であり、ターゲットの誤検出や未検出の原因となるという問題があった。   However, when the MUSIC method or Esprit method is applied to an on-vehicle radar device, if the relative relationship between the target and the radar is fast and the usable computing power is low, a large amount of data in the same situation is collected. It is impossible to do. Furthermore, when there are a plurality of reflecting objects around, it is difficult to secure a stable SN ratio due to the influence of multipath. Under these circumstances, it is difficult to accurately estimate the number of incoming waves, and there is a problem in that it may cause a false detection or no detection of the target.

本発明は上記点に鑑みて、MUSIC法やEsprit法など、固有値分解を用いて到来波推定を行った上で、到来波方向を推定するレーダ装置において、精度良く到来波数の推定が行えるようにすることを目的とし、新たなる手法を提案するものである。   In view of the above points, the present invention enables estimation of the number of incoming waves with high accuracy in a radar apparatus that estimates the direction of the incoming wave after estimating the incoming wave using eigenvalue decomposition such as the MUSIC method and the Esprit method. The purpose is to propose a new method.

上記目的を達成するため、請求項1に記載の発明では、ビート信号に基づいて、ターゲットの距離および相対速度(もしくは方位)を求め、それらを履歴情報として記憶する履歴情報記憶手段を有し、到来波数推定手段は、履歴情報記憶手段に記憶された履歴情報から、今回の到来波推定の際におけるターゲットの位置を予測するターゲット位置予測手段(400)を有し、該ターゲット位置予測手段により予測されるターゲットの位置からターゲットの複数が並ぶか否かを予測し、その予測結果に基づいて到来波の数の推定を行うことを特徴としている。   In order to achieve the above object, the invention according to claim 1 has history information storage means for obtaining the distance and relative speed (or direction) of the target based on the beat signal and storing them as history information, The arrival wave number estimation means has target position prediction means (400) for predicting the position of the target in the current arrival wave estimation from the history information stored in the history information storage means. It is characterized by predicting whether or not a plurality of targets are arranged from the position of the target to be estimated, and estimating the number of incoming waves based on the prediction result.

このように、前回のサイクルの際に得られたターゲットの履歴情報に基づいて、前回のサイクルの際に並走していなかったターゲット同士が今回のサイクルの際に並走するか否か、もしくは、前回のサイクルの際に並走していたターゲット同士が今回のサイクルの際に並走しなくなるか否かを予測し、その予測に基づいて到来波数を修正するようにしている。これにより、より正確に到来波数の推定を行うことが可能となる。   Thus, based on the target history information obtained during the previous cycle, whether targets that were not running in parallel during the previous cycle run in parallel during the current cycle, or In addition, it is predicted whether or not the targets that have been running in parallel in the previous cycle will not run in parallel in the current cycle, and the number of incoming waves is corrected based on the prediction. As a result, the number of incoming waves can be estimated more accurately.

例えば、請求項2に示されるように、到来波数推定手段は、ターゲット位置予測手段により予測されるターゲットの位置から今回の到来波推定の際に並ぶと予測されるターゲットの数を求め、この数と前回の到来波推定の際に並んでいたターゲットの数との差が、前回の到来波推定の際に推定された到来波の数と今回の到来波推定の際に推定された到来波の数との差と一致していない場合には、前回の到来波推定の際に推定された到来波の数に対して、今回の到来波推定の際に並ぶと予測されるターゲットの数と前回の到来波推定の際に並んでいたターゲットの数との差を加算もしくは減算することで、到来波の数の推定を行う。   For example, as described in claim 2, the arrival wave number estimation means obtains the number of targets predicted to be arranged in the current arrival wave estimation from the target position predicted by the target position prediction means, and this number And the number of targets lined up at the time of the previous arrival wave estimation are the number of arrival waves estimated at the previous arrival wave estimation and the arrival wave estimated at the current arrival wave estimation. If the number of targets does not match the number of incoming signals, the number of incoming waves estimated during the previous incoming wave estimation is The number of incoming waves is estimated by adding or subtracting the difference from the number of targets lined up at the time of incoming wave estimation.

また、請求項3に示されるように、ターゲット位置予測手段により予測されるターゲットの位置からターゲットの複数が並ぶか否かを予測し、その予測結果に基づいてスレッショルド設定手段が設定するスレッショルド(λ)の値を調整するようにしても、請求項1または2と同様の効果を得ることができる。   According to a third aspect of the present invention, it is predicted whether or not a plurality of targets are arranged from the target positions predicted by the target position predicting unit, and the threshold (λ that is set by the threshold setting unit based on the prediction result) Even if the value of) is adjusted, the same effect as in the first or second aspect can be obtained.

なお、上記各手段の括弧内の符号は、後述する実施形態に記載の具体的手段との対応関係を示すものである。   In addition, the code | symbol in the bracket | parenthesis of each said means shows the correspondence with the specific means as described in embodiment mentioned later.

本発明の第1実施形態におけるレーダ装置のブロック構成を示す図である。It is a figure which shows the block configuration of the radar apparatus in 1st Embodiment of this invention. 図1に示すレーダ装置の信号処理部が実行するMUSICの処理のフローチャートである。It is a flowchart of the process of MUSIC which the signal processing part of the radar apparatus shown in FIG. 1 performs. 図1に示すレーダ装置の信号処理部が実行する到来波数推定処理のフローチャートである。It is a flowchart of the arrival wave number estimation process which the signal processing part of the radar apparatus shown in FIG. 1 performs. (a)は、ビート周波数[bin]と固有値λとの相関関係を示した図であり、(b)は、ビート周波数に応じた雑音とビート信号との関係を示した図である。(A) is the figure which showed the correlation of beat frequency [bin] and eigenvalue (lambda), (b) is the figure which showed the relationship between the noise according to a beat frequency, and a beat signal. 送信アンテナから出力されるミリ波帯の電波の出力方位エリアとそのエリア内において車両が並走できる最大数との関係を示した模式図である。It is the schematic diagram which showed the relationship between the output azimuth | direction area | region of the electromagnetic wave of the millimeter wave band output from a transmission antenna, and the maximum number which a vehicle can run in parallel in the area. 第2実施形態におけるレーダ装置の信号処理部が実行する到来波数推定処理のフローチャートを示したものである。10 is a flowchart of arrival wave number estimation processing executed by a signal processing unit of a radar apparatus according to a second embodiment. 前回のサイクルと今回のサイクルの自車両と他車両との位置関係を示した模式図である。It is the schematic diagram which showed the positional relationship of the own vehicle and other vehicle of the last cycle and this cycle. 第3実施形態におけるレーダ装置の信号処理部が実行する到来波数推定処理のフローチャートを示したものである。10 is a flowchart of arrival wave number estimation processing executed by a signal processing unit of a radar apparatus according to a third embodiment. 第4実施形態におけるレーダ装置の信号処理部が実行する到来波数推定処理のフローチャートを示したものである。10 is a flowchart of arrival wave number estimation processing executed by a signal processing unit of a radar apparatus according to a fourth embodiment. 第5実施形態における到来歯数推定処理のフローチャートである。It is a flowchart of the arrival tooth number estimation process in 5th Embodiment. 正しい到来波数の判定とその到来波の方向および受信パワーを求める処理のフローチャートである。It is a flowchart of the process which calculates | requires the determination of the correct number of incoming waves, the direction of the incoming waves, and reception power. 第6実施形態における正しい到来波数の判定とその到来波の方向および受信パワーを求める処理のフローチャートである。It is a flowchart of the process which calculates | requires the determination of the correct number of incoming waves, the direction of the incoming waves, and reception power in 6th Embodiment. 複数のアンテナが等間隔にアレイ状に並べられたアダプティブアンテナに対して、到来波が入力されるときの様子を示した図である。It is the figure which showed a mode when an incoming wave was input with respect to the adaptive antenna with which the several antenna was arranged in the array form at equal intervals.

以下、本発明の実施形態について図に基づいて説明する。なお、以下の各実施形態相互において、互いに同一もしくは均等である部分には、図中、同一符号を付してある。   Hereinafter, embodiments of the present invention will be described with reference to the drawings. In the following embodiments, the same or equivalent parts are denoted by the same reference numerals in the drawings.

(第1実施形態)
本発明の一実施形態を適用した車載用のレーダ装置をブロック構成を図1に示す。以下、この図を参照して本実施形態のレーダ装置の構成について説明する。
(First embodiment)
A block configuration of an in-vehicle radar device to which an embodiment of the present invention is applied is shown in FIG. Hereinafter, the configuration of the radar apparatus according to the present embodiment will be described with reference to FIG.

レーダ装置は、車両の進行方向、例えば前方に向けてミリ波帯の電波を出射すると共に、先行車両などの目標物体によって反射した電波を到来波として受信することで、その到来波の数および方位を求めると共に、目標物体までの距離および自車両に対する相対速度を求めるものである。   The radar device emits millimeter wave radio waves in the traveling direction of the vehicle, for example, forward, and receives the radio waves reflected by the target object such as the preceding vehicle as incoming waves. As well as the distance to the target object and the relative speed with respect to the host vehicle.

図1に示されるように、レーダ装置には、変調信号生成部1、電圧制御発振器2、分配器3、送信アンテナ4、受信部5、ミキサ6、A/Dコンバータ7、および、信号処理部8が備えられている。   As shown in FIG. 1, the radar apparatus includes a modulation signal generation unit 1, a voltage controlled oscillator 2, a distributor 3, a transmission antenna 4, a reception unit 5, a mixer 6, an A / D converter 7, and a signal processing unit. 8 is provided.

変調信号生成部1は、例えば、時間に対して直線的に周波数が増減するような変調を行わせる所定周期の三角波状の変調信号を生成するものである。この変調信号生成部1が生成する変調信号が電圧制御発信器2に出力されるようになっている。   For example, the modulation signal generation unit 1 generates a triangular wave-shaped modulation signal having a predetermined period for performing modulation such that the frequency linearly increases or decreases with respect to time. The modulation signal generated by the modulation signal generator 1 is output to the voltage control transmitter 2.

電圧制御発振器2は、送信信号となるミリ波帯の高周波信号を生成するものである。この電圧制御発振器2は、変調信号生成部1から出力された変調信号に基づいて高周波信号の周波数を変調信号に応じて制御し、分配器3に出力するようになっている。   The voltage controlled oscillator 2 generates a high frequency signal in the millimeter wave band as a transmission signal. The voltage-controlled oscillator 2 controls the frequency of the high-frequency signal according to the modulation signal based on the modulation signal output from the modulation signal generator 1 and outputs it to the distributor 3.

分配器3は、電圧制御発振器2から出力された送信信号を電力分配し、ローカル信号を生成するもので、この分配器3からのローカル信号がミキサ6に入力されるようになっている。   The distributor 3 distributes the power of the transmission signal output from the voltage controlled oscillator 2 and generates a local signal. The local signal from the distributor 3 is input to the mixer 6.

送信アンテナ4は、分配器3を経て、電圧制御発振器2から出力された送信信号で示されるミリ波帯の電波を車両の進行方向、例えば前方に向けて出力するものである。   The transmission antenna 4 outputs a millimeter wave band radio wave indicated by the transmission signal output from the voltage controlled oscillator 2 via the distributor 3 in the traveling direction of the vehicle, for example, forward.

受信部5は、k本のアンテナ5aが等間隔にアレイ状に並べられたアダプティブアンテナを構成すると共に、複数のアンテナ5aに入力された到来波の受信データ(受信信号)となるxi(t)、(i=1〜K)を得る受信部5bと、受信部5bが得た各受信データxi(t)それぞれを結合して高周波ミキサ6に伝えるためのスイッチ5cを備えた構成となっている。   The receiving unit 5 constitutes an adaptive antenna in which k antennas 5a are arranged in an array at equal intervals, and xi (t) that is received data (received signals) of incoming waves input to the plurality of antennas 5a. , (I = 1 to K), and a switch 5c for combining the received data xi (t) obtained by the receiving unit 5b and transmitting them to the high-frequency mixer 6. .

ミキサ6は、受信部5から伝えられる受信データxi(t)の結合と、分配器3から伝えられるローカル信号とを混合し、これらの信号の差の周波数成分であるビート信号を生成するものである。このとき生成されるビート信号の周波数がビート周波数と呼ばれるもので、送信信号の周波数が増加する時のビート周波数を上り変調時のビート周波数、送信信号の周波数が減少する時のビート周波数を下り変調時のビート周波数と呼び、FMCW方式による目標物体の距離および相対速度の演算に用いられる。   The mixer 6 mixes the combination of the received data xi (t) transmitted from the receiving unit 5 and the local signal transmitted from the distributor 3, and generates a beat signal that is a frequency component of the difference between these signals. is there. The frequency of the beat signal generated at this time is called the beat frequency. The beat frequency when the frequency of the transmission signal increases and the beat frequency when the frequency of the transmission signal decreases and the beat frequency when the frequency of the transmission signal decreases and downstream modulation. This is called the beat frequency of the hour, and is used to calculate the distance and relative speed of the target object by the FMCW method.

A/Dコンバータ7は、ミキサ6から出力されたアナログ値として示されるビート信号をデジタル値に変換するためのものである。   The A / D converter 7 is for converting the beat signal shown as an analog value output from the mixer 6 into a digital value.

信号処理部8は、CPU、ROM、RAM、I/Oなどを備えた周知のマイクロコンピュータによって構成され、ROMなどに記憶されたプログラムにしたがった処理を実行するものである。具体的には、信号処理部8は、A/Dコンバータ7によってデジタル値に変換されたビート信号に基づいてMUSICの処理に基づく到来波の数および方位の推定(算出)を行うと共に、FMCWによる距離や相対速度の算出を行うようになっている。なお、この信号処理部8で算出された到来波の数および方位、目標物体の距離や相対速度に関する情報は、例えばRAMに逐次記憶され、少なくとも次の算出サイクルまで履歴情報として残されるようになっている。なお、信号処理部8のうち、この履歴情報を記憶する部分が履歴情報算出手段に相当する。   The signal processing unit 8 is configured by a known microcomputer including a CPU, a ROM, a RAM, an I / O, and the like, and executes a process according to a program stored in the ROM or the like. Specifically, the signal processing unit 8 estimates (calculates) the number and directions of incoming waves based on the MUSIC process based on the beat signal converted into a digital value by the A / D converter 7 and uses FMCW. The distance and relative speed are calculated. Note that the information on the number and direction of incoming waves, the distance and relative velocity of the target object calculated by the signal processing unit 8 is sequentially stored in, for example, a RAM, and remains as history information at least until the next calculation cycle. ing. In the signal processing unit 8, the part that stores the history information corresponds to the history information calculation means.

次に、上記のように構成されるレーダ装置が実行する処理について、図2および図3に示される処理のフローチャートを参照して説明する。   Next, processing executed by the radar apparatus configured as described above will be described with reference to flowcharts of processing shown in FIGS.

図2は、MUSICの処理に基づく到来波の数や方位の推定を示したものであり、図3は、MUSICの処理中に実行される到来波数推定処理の詳細を示したものである。これら各図に示される処理はレーダ装置における信号処理部8で実行されるもので、レーダ装置の使用期間中、例えば、オートクルーズ制御を実行するためのスイッチ(図示せず)が押下されたときに実行される。   FIG. 2 shows the estimation of the number and direction of incoming waves based on the MUSIC process, and FIG. 3 shows the details of the incoming wave number estimation process executed during the MUSIC process. The processing shown in each of these figures is executed by the signal processing unit 8 in the radar apparatus. For example, when a switch (not shown) for executing auto-cruise control is pressed during the use period of the radar apparatus. To be executed.

まず、ステップ100では、受信データxi(t)の取得処理が実行される。この処理は、ビート信号から受信データxi(t)を抽出することによって行われる。   First, in step 100, a process for obtaining received data xi (t) is executed. This process is performed by extracting the received data xi (t) from the beat signal.

続く、ステップ110では、受信データxi(t)に基づいて自己相関行列Rxxが算出され、ステップ120では、ステップ110で算出された自己相関行列Rxxとなるときの固有値λ(λ1〜λK)および固有ベクトルEn、Esが算出される。なお、ここでいう固有ベクトルEnとは、それぞれ上述した雑音空間の固有値λを算出する上で掛けたウェイトWに相当するもので、固有ベクトルEsとは、信号空間の固有値λを算出する上で掛けたウェイトWに相当するものである。これら自己相関行列Rxxの算出手法、固有値λおよび固有ベクトルEn、Esの算出手法に関しては従来と同様であるため、ここでは説明を省略する。このように、信号処理部8における自己相関行列Rxxを算出する部分が自己相関行列算出手段に相当するもので、固有値を算出する部分が固有値算出手段に相当するものである。   In step 110, the autocorrelation matrix Rxx is calculated based on the received data xi (t). In step 120, the eigenvalue λ (λ1 to λK) and the eigenvector when the autocorrelation matrix Rxx calculated in step 110 is obtained. En and Es are calculated. Here, the eigenvector En here corresponds to the weight W multiplied in calculating the eigenvalue λ of the noise space described above, and the eigenvector Es is multiplied in calculating the eigenvalue λ of the signal space. This corresponds to the weight W. Since the calculation method of the autocorrelation matrix Rxx and the calculation method of the eigenvalue λ and the eigenvectors En and Es are the same as those in the related art, the description thereof is omitted here. As described above, the part for calculating the autocorrelation matrix Rxx in the signal processing unit 8 corresponds to the autocorrelation matrix calculation means, and the part for calculating the eigenvalue corresponds to the eigenvalue calculation means.

ステップ130では、到来波数推定処理が実行される。この到来波数推定処理では、図3のフローチャートに示される各処理が実行される。信号処理部8のうち、この到来波数推定処理を実行する部分が到来波数推定手段に相当するものである。   In step 130, an incoming wave number estimation process is executed. In this incoming wave number estimation process, each process shown in the flowchart of FIG. 3 is executed. The portion of the signal processing unit 8 that executes this arrival wave number estimation process corresponds to the arrival wave number estimation means.

到来波数推定処理が実行されると、まず、ステップ200において、ビート周波数ごとのスレッショルドλTH1〜λTH256の読み出しが行われる。ここでいうビート周波数ごとのスレッショルドλTH1〜λTH256とは、図4(a)で示されるビート周波数[bin]と固有値λとの相関関係から求められるものであり、例えば予め行った実験結果から求められる。例えば、ビート周波数に応じた雑音として、図4(b)に示されるようにFM−AM変換雑音や1/f雑音、熱雑音、路面からの反射などの不要信号が発生する。これらがビート周波数と相関関係があることから、各不要信号を足し合わせた値とビート信号との間に線を引くと、図4(a)に示されるスレッショルドλTH1〜λTH256が求まる。このスレッショルドλTH1〜λTH256は、例えば、信号処理部8のROMなどにマップとして記憶されている。したがって、ステップ200において、このマップとして記憶されているスレッショルドλTH1〜λTH256の読み出しが行われる。   When the arrival wave number estimation process is executed, first, in step 200, the thresholds λTH1 to λTH256 for each beat frequency are read. The thresholds λTH1 to λTH256 for each beat frequency here are obtained from the correlation between the beat frequency [bin] and the eigenvalue λ shown in FIG. 4A, and are obtained, for example, from the results of experiments performed in advance. . For example, as noise corresponding to the beat frequency, unnecessary signals such as FM-AM conversion noise, 1 / f noise, thermal noise, and reflection from the road surface are generated as shown in FIG. 4B. Since these are correlated with the beat frequency, the thresholds λTH1 to λTH256 shown in FIG. 4A are obtained by drawing a line between the sum of the unnecessary signals and the beat signal. The thresholds λTH1 to λTH256 are stored as maps in the ROM of the signal processing unit 8, for example. Accordingly, in step 200, the thresholds λTH1 to λTH256 stored as this map are read.

続くステップ210では、信号処理部8に入力されたビート信号のビート周波数に応じたスレッショルドλTHがマップに記憶されたビート周波数とスレッショルドλTH1〜λTH256との相関関係から読み取られ、そのときのビート周波数に応じたスレッショルドλTHとして設定される。以上のようなスレッショルドの設定を行う部分がスレッショルド設定手段に相当する。   In the following step 210, the threshold λTH corresponding to the beat frequency of the beat signal input to the signal processing unit 8 is read from the correlation between the beat frequency stored in the map and the thresholds λTH1 to λTH256, and the beat frequency at that time is read. The corresponding threshold λTH is set. The portion for setting the threshold as described above corresponds to the threshold setting means.

そして、ステップ220にて、上述したステップ120で求められた固有値λ(λ1〜λK)を入力し、ステップ230において、固有値λ(λ1〜λK)それぞれがスレッショルドλTHよりも大きいか否かが求められる。これにより、ステップ240およびステップ250に示されるように、各固有値λ(λ1〜λK)が信号空間の固有値λl〜λLと雑音空間の固有値λL+1〜λKとに分別される。なお、図3に示されるフローチャート中では、任意(x本目)の受信アンテナ5aの受信データをから求められた固有値λxがスレッショルドλthよりも大きいか否かとして示してあるが、これは、各固有値λ(λ1〜λK)を代表的に表したものである。   Then, in step 220, the eigenvalues λ (λ1 to λK) obtained in step 120 described above are input, and in step 230, it is determined whether or not each of the eigenvalues λ (λ1 to λK) is larger than the threshold λTH. . As a result, as shown in step 240 and step 250, each eigenvalue λ (λ1 to λK) is divided into eigenvalues λ1 to λL in the signal space and eigenvalues λL + 1 to λK in the noise space. In the flowchart shown in FIG. 3, the eigenvalue λx obtained from the received data of the arbitrary (xth) receiving antenna 5a is shown as whether or not the eigenvalue λx is larger than the threshold λth. λ (λ1 to λK) is representatively represented.

そして、ステップ260において、信号空間の固有値λ1〜λLの数から到来波数Lを推定する。これにより、雑音などの不要信号を考慮にいれた正確な到来波数Lの推定を行うことが可能となる。   In step 260, the number of incoming waves L is estimated from the number of eigenvalues λ1 to λL in the signal space. This makes it possible to accurately estimate the number of incoming waves L taking into account unnecessary signals such as noise.

このように到来波数の推定が行われると、ここで求められた到来波数に基づき、ステップ140において、MUSICスペクトルが算出される。このMUSICスペクトルの算出手法に関しても従来と同様である。   When the arrival wave number is estimated in this way, a MUSIC spectrum is calculated in step 140 based on the obtained arrival wave number. The method for calculating the MUSIC spectrum is the same as in the prior art.

続いて、ステップ150において、ステップ140で求められた本来の到来波のピーク波に基づいて到来波の方位が算出される。すなわち、到来波のピーク波の分布を見れば、到来波のピーク波に対応する受信アンテナの持つビームが指向する方位を求めることが可能となる。この方位が目標物体の存在する方位として算出される。   Subsequently, in step 150, the direction of the incoming wave is calculated based on the peak wave of the original incoming wave obtained in step 140. That is, by looking at the distribution of the peak wave of the arriving wave, it is possible to obtain the direction in which the beam of the receiving antenna corresponding to the peak wave of the arriving wave is directed. This direction is calculated as the direction in which the target object exists.

ステップ160では、受信パワーの算出が行われる。この受信パワーの情報を、FMCWにおける上り/下り両変調時のビート周波数のペアを探索する処理(ペアリング処理という名称で周知)の情報として用い、目標物体からの距離、相対速度を算出する。   In step 160, reception power is calculated. The received power information is used as information for a process of searching for a pair of beat frequencies in both up / down modulation in FMCW (known as a pairing process), and a distance from the target object and a relative speed are calculated.

以上説明したように、本実施形態では、到来波数推定処理の際に、雑音などの不要信号を考慮してスレッショルドλTHを設定し、このスレッショルドλTHに基づいて到来波の固有値λ1〜λKが信号空間の固有値λ1〜λLか雑音空間の固有値λL+1〜λKかを分別するようにしている。   As described above, in the present embodiment, in the arrival wave number estimation process, the threshold λTH is set in consideration of unnecessary signals such as noise, and the eigenvalues λ1 to λK of the incoming wave are determined based on the threshold λTH in the signal space. The eigenvalues λ1 to λL of the noise and the eigenvalues λL + 1 to λK of the noise space are discriminated.

これにより、雑音を考慮に入れた上で、信号空間の固有値λ1〜λLの数から到来波数Lを正確に推定することができる。したがって、正確に推定された到来波数Lに基づいて、オートクルーズ制御を実行することが可能となり、より適切なオートクルーズ制御の実行が可能となる。   As a result, the number of incoming waves L can be accurately estimated from the number of eigenvalues λ1 to λL in the signal space in consideration of noise. Therefore, it is possible to execute the auto cruise control based on the accurately estimated arrival wave number L, and it is possible to execute more appropriate auto cruise control.

(第2実施形態)
本発明の第2実施形態について説明する。本実施形態のレーダ装置は、第1実施形態に対して信号処理部8が実行する到来波数推定処理を変更したものであり、その他は第1実施形態と同様であるため、異なる部分についてのみ説明する。
(Second Embodiment)
A second embodiment of the present invention will be described. The radar apparatus of this embodiment is obtained by changing the arrival wave number estimation process executed by the signal processing unit 8 with respect to the first embodiment, and the others are the same as those of the first embodiment, and therefore only different parts will be described. To do.

本実施形態は、上記のように、雑音を考慮した到来波推定に加えて、送信アンテナ4から出力されるミリ波帯の電波の出力方位エリアに応じて推定される到来波数の最大値を制限するものである。これについて、図5を参照して説明する。   As described above, the present embodiment limits the maximum value of the number of incoming waves estimated according to the output direction area of the radio wave in the millimeter wave band output from the transmitting antenna 4 in addition to the incoming wave estimation considering noise. To do. This will be described with reference to FIG.

図5は、送信アンテナ4から出力されるミリ波帯の電波の出力方位エリアとそのエリア内において車両が並走できる最大数との関係を示した模式図である。   FIG. 5 is a schematic diagram showing a relationship between an output azimuth area of millimeter wave radio waves output from the transmission antenna 4 and a maximum number of vehicles that can run in parallel in the area.

この図に示されるように、出力方位エリアの横幅は自車両からの距離から一義的に決まるため、その出力方位エリア内で並走可能な車両(ターゲット)の数が必然的に決まってくる。例えば、出力方位エリアの角度が車両前方を中心とした左右10度程度であった場合、車両の横幅が1.7m程度、並走している車両間のピッチが3.1m程度であると想定すると、自車両から17〜18m程度先であれば3台までしか複数の車両が並走できない。したがって、出力方位エリアごとに到来波数の最大値Lmaxが決まり、それを超えるような到来波数が推定されるような場合には、それは適切であるとは言えない。このため、本実施形態では、出力方位エリアごとに到来波数の最大値Lmaxを制限する。   As shown in this figure, since the width of the output azimuth area is uniquely determined from the distance from the own vehicle, the number of vehicles (targets) that can run in parallel in the output azimuth area is inevitably determined. For example, when the angle of the output azimuth area is about 10 degrees on the left and right with respect to the front of the vehicle, it is assumed that the lateral width of the vehicle is about 1.7 m and the pitch between parallel running vehicles is about 3.1 m. Then, if it is about 17-18m away from the own vehicle, only a plurality of vehicles can run in parallel. Therefore, when the maximum value Lmax of the number of incoming waves is determined for each output azimuth area and the number of incoming waves exceeding the maximum value is estimated, it cannot be said that it is appropriate. For this reason, in this embodiment, the maximum value Lmax of the number of incoming waves is limited for each output azimuth area.

図6は、本実施形態のレーダ装置における信号処理部8が実行する到来波数推定処理のフローチャートを示したものである。   FIG. 6 shows a flowchart of the arrival wave number estimation process executed by the signal processing unit 8 in the radar apparatus of this embodiment.

まず、ステップ200〜250において、上記第1実施形態と同様の処理が実行される。そして、ステップ300において、出力方位エリアごとの到来波数の最大値Lmax(Lmax1、Lmax2・・・)の読み出しが行われる。この最大値Lmaxは、上述したような考察に基づいて自車両からの距離に応じて決められており、マップとして信号処理部8に記憶されている。このため、信号処理部8に記憶されたマップの読み出しが行われる。   First, in steps 200 to 250, processing similar to that in the first embodiment is executed. In step 300, the maximum value Lmax (Lmax1, Lmax2,...) Of the number of incoming waves for each output azimuth area is read out. The maximum value Lmax is determined according to the distance from the host vehicle based on the above-described consideration, and is stored in the signal processing unit 8 as a map. For this reason, the map stored in the signal processing unit 8 is read.

続く、ステップ310において、読み出したマップに基づいて、最大値Lmaxが設定される。具体的には、レーダ装置では、図2に示したように受信パワーに基づいてFMCW方式による距離および相対速度の演算が行われている。このため、ここでの演算結果をターゲットとなる目標物体のエリア情報(距離情報)として用い、そのエリア情報に応じた最大値Lmaxが設定される。以上のような、到来波の数の最大値を設定する部分が到来波最大値設定手段に相当する。   In step 310, the maximum value Lmax is set based on the read map. Specifically, in the radar apparatus, distance and relative speed are calculated by the FMCW method based on the received power as shown in FIG. For this reason, the calculation result here is used as area information (distance information) of the target object to be a target, and a maximum value Lmax corresponding to the area information is set. The part for setting the maximum value of the number of incoming waves as described above corresponds to the incoming wave maximum value setting means.

そして、ステップ320において、ステップ240のように分別された信号空間の固有値λ1〜λLから求められる到来波数Lが最大値Lmax以下になっているか否かが判定される。ここで肯定判定された場合には、到来波数Lが物理的に並走しうる車両の数以下になっており、矛盾していないものとして、ステップ260において、そのときの到来波数Lがそのまま到来波数Lとして出力される。   In step 320, it is determined whether or not the number of incoming waves L obtained from the eigenvalues λ1 to λL of the signal space sorted as in step 240 is equal to or less than the maximum value Lmax. If an affirmative determination is made here, the incoming wave number L is equal to or less than the number of vehicles that can physically run in parallel, and it is assumed that there is no contradiction. Output as wave number L.

逆に、ここで否定判定された場合には、到来波数Lが物理的に並走しうる車両の数以上になっており、矛盾しているものとして、ステップ330に進む。そして、このステップでLmaxが到来波数Lとして設定され、ステップ260において出力される。   On the other hand, if a negative determination is made here, the number of incoming waves L is equal to or greater than the number of vehicles that can physically run in parallel, and the process proceeds to step 330 as inconsistent. In this step, Lmax is set as the incoming wave number L and is output in step 260.

以上説明したように、本実施形態では、出力方位エリア内に物理的に並走しうる車両の数を到来波数Lの最大値Lmaxとして制限し、それを超えるような到来波数が推定されるような場合には、矛盾しているものとして、最大値Lmaxを到来波数Lとして設定するようにしている。これにより、より正確に到来波数Lの推定を行うことが可能となる。   As described above, in the present embodiment, the number of vehicles that can physically run in the output azimuth area is limited as the maximum value Lmax of the number of incoming waves L, and the number of incoming waves exceeding that is estimated. In such a case, the maximum value Lmax is set as the number of incoming waves L as contradictory. This makes it possible to estimate the arrival wave number L more accurately.

(第3実施形態)
本発明の第3実施形態について説明する。本実施形態のレーダ装置は、第1実施形態に対して信号処理部8が実行する到来波数推定処理を変更したものであり、その他は第1実施形態と同様であるため、異なる部分についてのみ説明する。
(Third embodiment)
A third embodiment of the present invention will be described. The radar apparatus of this embodiment is obtained by changing the arrival wave number estimation process executed by the signal processing unit 8 with respect to the first embodiment, and the others are the same as those of the first embodiment, and therefore only different parts will be described. To do.

本実施形態は、上記第2実施形態に加え、レーダ装置の履歴情報から次サイクルで並走するか否かを予測し、その予測結果に応じた到来波数推定を行うものである。これについて図7を参照して説明する。   In the present embodiment, in addition to the second embodiment, whether or not to run in parallel in the next cycle is predicted from the history information of the radar apparatus, and the number of incoming waves is estimated according to the prediction result. This will be described with reference to FIG.

図7は、複数の先行車両がある場合に、前回のサイクルにおいて、自車両から近い側に位置していた他車両が遠い側に位置していた他車両から離れていたものの、今回のサイクルにおいて、その近い側に位置していた他車両が遠い側に位置していた他車両と並走するような場合を示したものである。   FIG. 7 shows that when there are a plurality of preceding vehicles, in the previous cycle, the other vehicle located on the side closer to the own vehicle was separated from the other vehicle located on the far side. This shows a case where the other vehicle located on the near side runs in parallel with the other vehicle located on the far side.

前回のサイクルの時には、自車両と各他車両との距離および相対速度で決まるビート周波数が異なるため、それぞれの他車両からの到来波は違うエリアからのものとして検出され、それぞれのエリアでそれぞれの到来波数が推定されることになる。つまり、自車両から近い側の他車両が存在するエリアにおいて到来波数が1、自車両から遠い側の他車両が存在するエリアにおいて到来波数が1という形態で到来波数が推定される。   In the previous cycle, the beat frequency determined by the distance and relative speed between the host vehicle and each other vehicle is different, so the incoming waves from each other vehicle are detected as coming from different areas, The number of incoming waves will be estimated. That is, the arrival wave number is estimated in a form where the arrival wave number is 1 in the area where the other vehicle nearer to the own vehicle is present and the arrival wave number is 1 in the area where the other vehicle far from the own vehicle exists.

これに対し、今回のサイクルの時には、自車両と2台の他車両が両方との距離および相対速度で決まるビート周波数が一致し、2台の他車両からの到来波が同じエリアのものとして到来波数が推定されることになる。   On the other hand, at the time of this cycle, the beat frequencies determined by the distance and relative speed between the host vehicle and the two other vehicles are the same, and the incoming waves from the two other vehicles arrive in the same area. The wave number will be estimated.

したがって、前回のサイクルの時に求められた先行車両の距離および相対速度から、今回のサイクルの時には、先行車両が他の先行車両と並走するであろうと予測される場合には、それを到来波数推定に用いるようにする。   Therefore, if it is predicted from the distance and relative speed of the preceding vehicle obtained during the previous cycle that the preceding vehicle will run in parallel with other preceding vehicles during the current cycle, this is calculated as the number of incoming waves. Used for estimation.

図8は、本実施形態のレーダ装置における信号処理部8が実行する到来波数推定処理のフローチャートを示したものである。   FIG. 8 shows a flowchart of the arrival wave number estimation process executed by the signal processing unit 8 in the radar apparatus of this embodiment.

まず、ステップ200〜250およびステップ300〜330において、上記第2実施形態と同様の処理が実行される。そして、ステップ400において、レーダ装置で前回のサイクルの時に求められたターゲット(先行車両)の履歴情報から、今回のサイクルにおけるターゲットの位置予測が行われる。これにより、前回のサイクルで並走していなかった他車両同士が今回のサイクルで並走すると予測される場合には、前回の到来波数を1つ加算したものが予測された到来波数として設定されることになる。なお、このようなターゲット位置の予測を行う部分がターゲット位置予測手段に相当するものである。   First, in steps 200 to 250 and steps 300 to 330, processing similar to that in the second embodiment is executed. In step 400, the target position prediction in the current cycle is performed from the history information of the target (preceding vehicle) obtained at the previous cycle by the radar apparatus. As a result, when other vehicles that were not running in parallel in the previous cycle are predicted to run in parallel in the current cycle, the value obtained by adding the previous number of incoming waves is set as the predicted number of incoming waves. Will be. Note that such a target position prediction portion corresponds to the target position prediction means.

そして、ステップ410において、前回の到来波数に対して1つ加算した値が設定されているかが判定されると共に、ステップ320もしくはステップ330で設定された到来波数Lが前回の到来波数に対して1つ加算された値になっていないか否かが判定される。これにより、前回の到来波数に対して1つ加算した値が設定されており、かつ、到来波数Lが前回の到来波数に対して1つ加算された値になっていない場合には、ステップ320もしくはステップ330で設定された到来波数Lが誤っているものとして、ステップ420に進んで到来波数LがL+1に変えられる。また、前回の到来波数に対して1つ加算した値が設定されていない場合、もしくは、前回の到来波数に対して1つ加算した値が設定されていても、到来波数Lが前回の到来波数に対して1つ加算された値になっている場合には、ステップ320もしくはステップ330で設定された到来波数Lが正しいものとして、そのままステップ260に進む。そして、到来波数Lがステップ260において出力される。   In step 410, it is determined whether a value obtained by adding one to the previous arrival wave number is set, and the arrival wave number L set in step 320 or step 330 is 1 with respect to the previous arrival wave number. It is determined whether or not a value obtained by adding the two is not reached. As a result, when a value obtained by adding one to the previous number of incoming waves is set and the incoming wave number L is not a value obtained by adding one to the previous number of incoming waves, step 320 is performed. Alternatively, assuming that the incoming wave number L set in step 330 is incorrect, the routine proceeds to step 420 where the incoming wave number L is changed to L + 1. In addition, when the value obtained by adding one to the previous arrival wave number is not set, or even when the value obtained by adding one to the previous arrival wave number is set, the arrival wave number L is the previous arrival wave number. If the value is a value obtained by adding one to, the number of incoming waves L set in step 320 or step 330 is assumed to be correct, and the process proceeds to step 260 as it is. Then, the incoming wave number L is output at step 260.

以上説明したように、本実施形態では、前回のサイクルの時に得られたターゲットの履歴情報に基づいて、前回のサイクルの時に並走していなかったターゲット同士が今回のサイクルの時に並走するかを予測し、その予測に基づいて到来波数Lを修正するようにしている。これにより、より正確に到来波数Lの推定を行うことが可能となる。   As described above, in this embodiment, based on the target history information obtained at the previous cycle, whether the targets that were not running at the previous cycle run in parallel at the current cycle. And the number of incoming waves L is corrected based on the prediction. This makes it possible to estimate the arrival wave number L more accurately.

(第4実施形態)
本発明の第4実施形態について説明する。本実施形態のレーダ装置は、第3実施形態と同様に、前回のサイクルの時に得られたターゲットの履歴情報に基づいて、前回のサイクルの時に並走していなかったターゲット同士が今回のサイクルの時に並走するかを予測し、その予測に基づいて到来波数Lを求めるものであるが、ここではスレッショルドλTHを下げることで、その予測に基づく到来波数推定を行う。
(Fourth embodiment)
A fourth embodiment of the present invention will be described. As in the third embodiment, the radar apparatus of this embodiment is based on the target history information obtained at the previous cycle, and targets that were not running in parallel at the previous cycle are Whether to run in parallel is predicted, and the number of incoming waves L is obtained based on the prediction. Here, the number of incoming waves is estimated based on the prediction by lowering the threshold λTH.

図9は、本実施形態のレーダ装置における信号処理部8が実行する到来波数推定処理のフローチャートを示したものである。   FIG. 9 shows a flowchart of the arrival wave number estimation process executed by the signal processing unit 8 in the radar apparatus of the present embodiment.

まず、ステップ200〜210において、上記第1実施形態と同様の処理が実行される。そして、ステップ500において、レーダ装置で前回のサイクルの時に求められたターゲット(先行車両)の履歴情報から、今回のサイクルにおけるターゲットの位置予測が行われる。これにより、前回のサイクルで並走していなかった他車両同士が今回のサイクルで並走すると予測される場合には、前回の到来波数を1つ加算したものが予測された到来波数として設定されることになる。   First, in steps 200 to 210, processing similar to that in the first embodiment is executed. In step 500, the target position in the current cycle is predicted from the history information of the target (preceding vehicle) obtained at the previous cycle by the radar apparatus. As a result, when other vehicles that were not running in parallel in the previous cycle are predicted to run in parallel in the current cycle, the value obtained by adding the previous number of incoming waves is set as the predicted number of incoming waves. Will be.

そして、ステップ510において、前回の到来波数に対して1つ加算した値が設定されているかが判定される。このステップで肯定判定された場合には、ステップ520に進み、スレッショルドλTHをαだけ下げたのちステップ230に進み、否定判定された場合には、そのままステップ230に進む。   In step 510, it is determined whether or not a value obtained by adding one to the previous number of incoming waves is set. If an affirmative determination is made in this step, the process proceeds to step 520, the threshold λTH is decreased by α, and then the process proceeds to step 230. If a negative determination is made, the process proceeds to step 230 as it is.

ここでいうαは、信号空間の固有値λ1〜λLの数が1つ分増えると想定される程度の値に設定されるもので、例えば実験的に求められる。   Here, α is set to such a value that the number of eigenvalues λ1 to λL in the signal space is assumed to increase by one, and is obtained experimentally, for example.

この後は、ステップ230において、固有値λ(λ1〜λK)それぞれがスレッショルドλTHよりも大きいか否かが求められる。これ以降は、第2実施形態と同様の処理が実行される。   Thereafter, in step 230, it is determined whether each eigenvalue λ (λ1 to λK) is larger than the threshold λTH. Thereafter, the same processing as in the second embodiment is executed.

以上説明したように、本実施形態も、前回のサイクルの時に得られたターゲットの履歴情報に基づいて、前回のサイクルの時に並走していなかったターゲット同士が今回のサイクルの時に並走するかを予測し、その予測に基づいてスレッショルドλTHを修正することで、到来波数Lを修正するようにしている。これにより、より正確に到来波数Lの推定を行うことが可能となる。   As described above, according to the present embodiment, based on the target history information obtained at the previous cycle, whether the targets that were not running at the previous cycle run in parallel at the current cycle. The number of incoming waves L is corrected by correcting the threshold λTH based on the prediction. This makes it possible to estimate the arrival wave number L more accurately.

(第5実施形態)
本発明の第5実施形態について説明する。本実施形態のレーダ装置は、第3実施形態とほぼ同様であるが、到来波数を任意の初期値(L=L0、例えばL0=2、3)に固定しておき、その初期値に基づいて仮の到来波数Lを求めると共に、受信パワーに基づいて仮の到来波数Lから実際の到来波数Lを最終的に求めるものである。
(Fifth embodiment)
A fifth embodiment of the present invention will be described. The radar apparatus of this embodiment is almost the same as that of the third embodiment, but the number of incoming waves is fixed to an arbitrary initial value (L = L0, for example, L0 = 2, 3), and based on the initial value. The temporary arrival wave number L is obtained, and the actual arrival wave number L is finally obtained from the temporary arrival wave number L based on the received power.

図10は、本実施形態のレーダ装置における信号処理部8が実行する到来波数推定処理のフローチャートを示したものである。本実施形態の場合、この到来波数推定処理によって、まず仮の到来波数Lが求められることになる。なお、信号処理部8のうち、この処理を実行する部分が仮到来波数推定手段に相当する。   FIG. 10 shows a flowchart of the arrival wave number estimation process executed by the signal processing unit 8 in the radar apparatus of this embodiment. In the case of the present embodiment, a temporary incoming wave number L is first obtained by this incoming wave number estimation process. In addition, the part which performs this process among the signal processing parts 8 is equivalent to temporary arrival wave number estimation means.

ステップ220では、上記第3実施形態と同様に、固有値λ(λ1〜λK)が入力される。そして、ステップ600にて、到来波数Lの初期値が設定され、到来波数L=L0として到来波数Lの初期値が固定される。これに基づき、ステップ610では、上述した図3に示したステップ240、250と同様に、信号空間の固有値λ(λ1〜λL)と雑音空間の固有値λ(λL+1〜λK)とが求められる。   In step 220, the eigenvalue λ (λ1 to λK) is input as in the third embodiment. In step 600, the initial value of the incoming wave number L is set, and the initial value of the incoming wave number L is fixed as the incoming wave number L = L0. Based on this, in step 610, the eigenvalue λ (λ1 to λL) of the signal space and the eigenvalue λ (λL + 1 to λK) of the noise space are obtained in the same manner as in steps 240 and 250 shown in FIG.

この後、第3実施形態で示したステップ300〜330およびステップ400から420およびステップ260の処理が実行されることで、到来波数Lが推定される。そして、このようにして求められた値が仮の到来波数Lとされる。これにより、到来波数推定処理が終了となる。   Thereafter, the processing of steps 300 to 330 and steps 400 to 420 and step 260 shown in the third embodiment is executed, so that the incoming wave number L is estimated. And the value calculated | required in this way is made into the temporary number L of incoming waves. Thereby, the arrival wave number estimation process is completed.

そして、この到来波数推定処理が完了すると、図2のステップ140に示されるように、求められた仮の到来波数に基づき、従来の手法を用いてMUSICスペクトルが算出される。なお、信号処理部8のうち、この処理を実行する部分がMISICスペクトル算出手段に相当する。   When this arrival wave number estimation process is completed, a MUSIC spectrum is calculated using a conventional method based on the calculated temporary arrival wave number, as shown in step 140 of FIG. In addition, the part which performs this process among the signal processing parts 8 is equivalent to a MISIC spectrum calculation means.

図11は、本実施形態のレーダ装置における正しい到来波数Lの判定とその到来波の方向および受信パワーを求める処理のフローチャートである。この処理は、図2におけるステップ150およびステップ160の各処理に相当するものとして実行されるものである。なお、信号処理部8のうち、この処理を実行する部分が方位検出手段および受信パワー検出手段に相当する。   FIG. 11 is a flowchart of processing for determining the correct number of incoming waves L and obtaining the direction of the incoming waves and the received power in the radar apparatus of this embodiment. This process is executed as equivalent to each of the processes in step 150 and step 160 in FIG. Of the signal processing unit 8, the part that executes this processing corresponds to the azimuth detecting means and the received power detecting means.

まず、ステップ700では、ステップ140で求められたMUSICスペクトルから確認される到来波のピーク波の方位(以下、ピークサーチ方位という)D1〜DLと受信パワーP1〜PLが検出される。なお、ここでは、上述したステップ260において、仮の到来波数Lとされていることから、ピーク波の数もそれと同じLとされる。なお、ピークサーチ方位D1〜DLと受信パワーP1〜PLの算出手法に関しては、上述したステップ150、160と同様である。   First, in step 700, peak wave directions (hereinafter referred to as peak search directions) D1 to DL and received powers P1 to PL confirmed from the MUSIC spectrum obtained in step 140 are detected. Here, since the number of provisional incoming waves is set to L in step 260 described above, the number of peak waves is also set to the same L. The method for calculating the peak search directions D1 to DL and the received powers P1 to PL is the same as in steps 150 and 160 described above.

次に、ステップ710では、ビート周波数ごとの受信パワーのスレッショルドPTH1〜PTH256の読み出しが行われる。ここでいうビート周波数ごとの受信パワーのスレッショルドPTH1〜PTH256とは、ビート周波数[bin]と受信パワーとの相関関係から求められるものであり、例えば予め行った実験結果から求められる。   Next, in step 710, the received power thresholds PTH1 to PTH256 for each beat frequency are read. The thresholds PTH1 to PTH256 of the reception power for each beat frequency here are obtained from the correlation between the beat frequency [bin] and the reception power, and are obtained from, for example, a result of an experiment performed in advance.

続くステップ720では、信号処理部8に入力されたビート信号のビート周波数に応じた受信パワーのスレッショルドPTHがマップに記憶されたビート周波数と受信パワーのスレッショルドPTH1〜PTH256との相関関係から読み取られ、そのときのビート周波数に応じた受信パワーのスレッショルドPTHとして設定される。以上のようなスレッショルドの設定を行う部分がスレッショルド設定手段に相当する。   In the following step 720, the received power threshold PTH corresponding to the beat frequency of the beat signal input to the signal processing unit 8 is read from the correlation between the beat frequency stored in the map and the received power thresholds PTH1 to PTH256, It is set as a threshold PTH of received power corresponding to the beat frequency at that time. The portion for setting the threshold as described above corresponds to the threshold setting means.

そして、ステップ730において、ステップ700で求められた受信パワーP1〜PLがステップ720で設定されたビート周波数に応じた受信パワーのスレッショルドPTHを超えているか否かが判定され、受信パワーP1〜PLのうちスレッショルドPTHを超えていたもののみが実際に存在する到来波であるものとされる。そして、その実際の到来波の方位と受信パワーとの情報を用いて、目標物体の存在する方位および目標物体からの距離、相対速度が求められる。   In step 730, it is determined whether or not the received power P1 to PL obtained in step 700 exceeds the received power threshold PTH corresponding to the beat frequency set in step 720. Of these, only the wave exceeding the threshold PTH is regarded as an actually existing incoming wave. Then, using the information on the actual direction of the incoming wave and the received power, the direction in which the target object exists, the distance from the target object, and the relative velocity are obtained.

以上説明したように、本実施形態では、到来波数Lを任意の初期値に固定しておき、この初期値に基づいて仮の到来波数Lを求めたのち、仮の到来波の各方位の受信パワーがスレッショルドPTHを超えているか否か判定し、超えていないものを到来波から除くことで、実際の到来波数Lを求めるようにしている。これにより、実際の到来波数Lを正確に求めることが可能となり、上記第1実施形態と同様の効果を得ることが可能となる。   As described above, in the present embodiment, the incoming wave number L is fixed to an arbitrary initial value, and the temporary incoming wave number L is obtained based on the initial value, and then reception of each direction of the temporary incoming wave is received. It is determined whether or not the power exceeds the threshold PTH, and those that do not exceed the threshold are removed from the incoming wave, thereby obtaining the actual incoming wave number L. As a result, the actual number of incoming waves L can be obtained accurately, and the same effect as in the first embodiment can be obtained.

(第6実施形態)
本発明の第6実施形態について説明する。本実施形態のレーダ装置は、第5実施形態に対して、図11に示した処理に変えて図12に示す処理を実行するものである。その他に関しては第5実施形態と同様であるため、説明を省略する。
(Sixth embodiment)
A sixth embodiment of the present invention will be described. The radar apparatus according to the present embodiment executes the process shown in FIG. 12 instead of the process shown in FIG. 11 with respect to the fifth embodiment. Since other aspects are the same as those of the fifth embodiment, description thereof is omitted.

図12は、本実施形態のレーダ装置における正しい到来波数Lの判定とその到来波の方向および受信パワーを求める処理のフローチャートである。この処理は、図2におけるステップ150およびステップ160の各処理に相当するものとして実行されるものである。   FIG. 12 is a flowchart of processing for determining the correct number of incoming waves L and obtaining the direction of the incoming waves and the received power in the radar apparatus of this embodiment. This process is executed as equivalent to each of the processes in step 150 and step 160 in FIG.

まず、ステップ700では、上記ステップ140で求められたMUSICスペクトルからピークサーチ方位D1〜DLが検出される。   First, in step 700, peak search directions D1 to DL are detected from the MUSIC spectrum obtained in step 140.

次に、ステップ740では、各ピークサーチ方位D1〜DLでの受信パワーP1〜PLが算出される。具体的には、受信パワーP1〜PLは次式により求められ、Sの対角成分から抽出される(菊間信良、アダプティブアンテナ技術、オーム社、P.137−P.141参照)。   Next, in step 740, received powers P1 to PL at the peak search directions D1 to DL are calculated. Specifically, the received powers P1 to PL are obtained by the following equation and extracted from the diagonal components of S (see Nobuyoshi Kikuma, Adaptive Antenna Technology, Ohmsha, P.137-P.141).

(数5)
S=(AHA)-1H(Rxx−σ2I)A(AHA)-1
ただし、Sは信号(波源)相関行列、Aは方向行列、Rxxは到来波の受信データxi(t)の自己相関行列、Iは単位行列である。
(Equation 5)
S = (A H A) −1 A H (Rxx−σ 2 I) A (A H A) −1
Here, S is a signal (wave source) correlation matrix, A is a direction matrix, Rxx is an autocorrelation matrix of received data xi (t) of incoming waves, and I is a unit matrix.

また、ステップ710および720にて、第5実施形態と同様の処理が行われる。そして、ステップ730において、ステップ750で求められた受信パワーP1〜PLがステップ720で設定されたビート周波数に応じた受信パワーのスレッショルドPTHを超えているか否かが判定され、受信パワーP1〜PLのうちスレッショルドPTHを超えていたもののみが実際に存在する到来波であるものとされる。その実際の到来波の方位と受信パワーとの情報を用いて、目標物体の存在する方位および目標物体からの距離、相対速度が求められる。   In steps 710 and 720, the same processing as in the fifth embodiment is performed. In step 730, it is determined whether or not the received power P1 to PL obtained in step 750 exceeds the received power threshold PTH corresponding to the beat frequency set in step 720. Of these, only the wave exceeding the threshold PTH is regarded as an actually existing incoming wave. Using the information on the actual direction of the incoming wave and the received power, the direction in which the target object exists, the distance from the target object, and the relative velocity are obtained.

以上説明したように、本実施形態でも、第5実施形態と同様に、到来波数Lを任意の初期値に固定しておき、この初期値に基づいて仮の到来波数Lを求めたのち、仮の到来波の各方位の受信パワーがスレッショルドPTHを超えているか否か判定し、超えていないものを到来波から除くことで、実際の到来波数Lを求めるようにしている。これにより、実際の到来波数Lを正確に求めることが可能となり、上記第1実施形態と同様の効果を得ることが可能となる。   As described above, in this embodiment as well, as in the fifth embodiment, the incoming wave number L is fixed to an arbitrary initial value, and the temporary incoming wave number L is obtained based on this initial value. It is determined whether or not the received power in each direction of the incoming wave exceeds the threshold PTH, and those not exceeding are removed from the incoming wave to obtain the actual incoming wave number L. As a result, the actual number of incoming waves L can be obtained accurately, and the same effect as in the first embodiment can be obtained.

(他の実施形態)
上記第1〜第4実施形態では、車載用のレーダ装置を例に挙げて説明したが、他のレーダ装置に対しても本発明を適用することができる。ただし、高速で移動する車両に搭載され、かつ、高いデバイス性能が期待できない車載用のレーダ装置に本発明を適用すると特に好適である。
(Other embodiments)
In the first to fourth embodiments, the on-vehicle radar device has been described as an example, but the present invention can also be applied to other radar devices. However, it is particularly preferable to apply the present invention to an on-vehicle radar device that is mounted on a vehicle that moves at high speed and that cannot be expected to have high device performance.

第3、第4実施形態では、前回のサイクルの時に得られたターゲットの履歴情報に基づいて、前回のサイクルの時に並走していなかったターゲット同士が今回のサイクルの時に並走するかを予測し、その予測に基づいて到来波数Lを修正するようにしている。これに対し、前回のサイクルの時に得られたターゲットの履歴情報に基づいて、前回のサイクルの時に並走していたターゲット同士が今回のサイクルの時に並走しないことを予測し、その予測に基づいて到来波数Lを修正することも可能である。   In the third and fourth embodiments, based on the target history information obtained during the previous cycle, it is predicted whether the targets that were not running in parallel during the previous cycle will run in parallel during the current cycle. The incoming wave number L is corrected based on the prediction. On the other hand, based on the target history information obtained during the previous cycle, it is predicted that the targets that were running in parallel during the previous cycle will not run in parallel during the current cycle. It is also possible to correct the number of incoming waves L.

この場合、第3実施形態では到来波数LがL−1に修正され、第4実施形態ではスレッショルドλTHが上げられることで固有値λの数が減らされることになる。   In this case, in the third embodiment, the number of incoming waves L is corrected to L−1, and in the fourth embodiment, the number of eigenvalues λ is reduced by increasing the threshold λTH.

もちろん、ここでは並走するターゲットの数が1つ増減する場合についてのみ説明しているが、複数の場合にはその数分だけ増減することになる。   Of course, only the case where the number of parallel running targets increases or decreases is described here, but in the case of a plurality of targets, the number increases or decreases by that number.

つまり、今回のサイクルの際に到来波推定の際に並ぶと予測されるターゲットの数を求め、この数と前回のサイクルの際に並んでいたターゲットの数との差を求める。そして、前回のサイクルの際に推定された到来波の数と今回のサイクルの際に推定された到来波の数との差と一致していない場合には、前回のサイクルの際に推定された到来波の数に対して、今回のサイクルの際に並ぶと予測されるターゲットの数と前回のサイクルの際に並んでいたターゲットの数との差を加算もしくは減算する。これにより、並走するターゲットの数に応じた到来波の数の推定を行うことが可能となる。   That is, the number of targets predicted to be arranged in the arrival wave estimation in the current cycle is obtained, and the difference between this number and the number of targets arranged in the previous cycle is obtained. And if it does not match the difference between the number of incoming waves estimated during the previous cycle and the number of incoming waves estimated during the current cycle, it was estimated during the previous cycle. The difference between the number of targets expected to be arranged in the current cycle and the number of targets arranged in the previous cycle is added to or subtracted from the number of incoming waves. This makes it possible to estimate the number of incoming waves according to the number of parallel targets.

また、上記各実施形態では、MUSIC法を例に挙げて説明したが、固有値分解を用いて到来波数の推定を行った上で、到来波方向を推定するレーダ装置、例えばEsprit法、ユニタリMUSIC法、ユニタリEsprit法などが用いられるレーダ装置についても本発明を適用することが可能である。   Further, in each of the above embodiments, the MUSIC method has been described as an example. However, a radar device that estimates the direction of the incoming wave after estimating the number of incoming waves using eigenvalue decomposition, such as the Esprit method and the unitary MUSIC method The present invention can also be applied to a radar apparatus using the unitary Esprit method or the like.

さらに、上記実施形態では、FMCWにてビート信号から距離や相対速度を求める場合について例示したが、これに限るものではなく、2周波CWやパルス方式、スペクトラム拡散方式を採用することも可能である。そして、FMCW方式が適用される場合には、ビート周波数ごとのスレッショルドλTH1〜λTH256が設定されることになるが、2周波方式や多周波CW方式が適用される場合には、ビート信号の位相ごとに、パルス方式やスペクトラム拡散方式が適用される場合には、ビート信号の遅延時間、すなわち物標までの距離ごとにスレッショルドが設定されることになる。   Further, in the above embodiment, the case where the distance and relative speed are obtained from the beat signal by FMCW is exemplified, but the present invention is not limited to this, and it is also possible to adopt a two-frequency CW, a pulse method, and a spread spectrum method. . When the FMCW method is applied, the thresholds λTH1 to λTH256 for each beat frequency are set. When the two-frequency method or the multi-frequency CW method is applied, for each phase of the beat signal. In addition, when the pulse method or the spread spectrum method is applied, the threshold is set for each delay time of the beat signal, that is, for each distance to the target.

そして、上記実施形態では、ビート周波数ごとのスレッショルドλTH1〜λTH256、PTH1〜PTH256というように、256個のスレッショルドλTH、PTHを設定する場合について説明したが、スレッショルドλTH、PTHの数は任意の数nであり、ビート周波数ごとのスレッショルドλTH1〜λTHn、PTH1〜PTHnとすることができる。   In the above embodiment, the case where 256 thresholds λTH and PTH are set, such as the thresholds λTH1 to λTH256 and PTH1 to PTH256 for each beat frequency, has been described, but the number of thresholds λTH and PTH is an arbitrary number n. The thresholds λTH1 to λTHn and PTH1 to PTHn for each beat frequency can be set.

なお、各図中に示したステップは、各種処理を実行する手段に対応するものである。   The steps shown in each figure correspond to means for executing various processes.

1…変調信号生成部、2…電圧制御発振器、3…分配器、4…送信アンテナ、5…受信部、6…ミキサ、7…A/Dコンバータ、8…信号処理部。   DESCRIPTION OF SYMBOLS 1 ... Modulation signal production | generation part, 2 ... Voltage control oscillator, 3 ... Distributor, 4 ... Transmission antenna, 5 ... Reception part, 6 ... Mixer, 7 ... A / D converter, 8 ... Signal processing part.

Claims (3)

送信アンテナ(4)を通じて送信信号を送信すると共に、互いに異なる位置に配設された複数のアンテナ(5a)を介して複数の方向から到来する複数の到来波を受信し、各アンテナ毎の受信信号から、前記送信信号を反射したターゲットからの到来波の方向の推定を行うレーダ装置において、
前記送信信号と前記受信信号の差の周波数成分であるビート信号を生成するビート信号手段(6)と、
前記ビート信号に基づいて自己相関行列(Rxx)を算出する自己相関行列算出手段(110)と、
前記自己相関行列算出手段で算出された自己相関行列から固有値(λ)を算出する固有値算出手段(120)と、
前記固有値算出手段で算出された固有値(λ)に基づいて前記到来波の数の推定を行う到来波数推定手段(130)と、
前記ビート信号に基づいて、前記ターゲットの距離および相対速度を求め、それらを履歴情報として記憶する履歴情報記憶手段を有し、
前記到来波数推定手段は、前記履歴情報記憶手段に記憶された前記履歴情報から、今回の到来波推定の際における前記ターゲットの位置を予測するターゲット位置予測手段(400)を有し、該ターゲット位置予測手段により予測される前記ターゲットの位置から前記ターゲットの複数が並ぶか否かを予測し、その予測結果に基づいて前記到来波の数の推定を行うようになっていることを特徴とするレーダ装置。
A transmission signal is transmitted through the transmission antenna (4), and a plurality of incoming waves arriving from a plurality of directions are received via a plurality of antennas (5a) arranged at different positions, and a reception signal for each antenna is received. From the radar apparatus for estimating the direction of the incoming wave from the target reflecting the transmission signal,
Beat signal means (6) for generating a beat signal that is a frequency component of a difference between the transmission signal and the reception signal;
Autocorrelation matrix calculating means (110) for calculating an autocorrelation matrix (Rxx) based on the beat signal;
Eigenvalue calculating means (120) for calculating the eigenvalue (λ) from the autocorrelation matrix calculated by the autocorrelation matrix calculating means;
Arriving wave number estimating means (130) for estimating the number of incoming waves based on the eigenvalue (λ) calculated by the eigenvalue calculating means;
Based on the beat signal, the distance and relative speed of the target are obtained, and history information storage means for storing them as history information is provided.
The arrival wave number estimation means includes target position prediction means (400) for predicting the position of the target at the time of current arrival wave estimation from the history information stored in the history information storage means, and the target position A radar which predicts whether or not a plurality of targets are arranged from the target positions predicted by a prediction means, and estimates the number of incoming waves based on the prediction result apparatus.
前記到来波数推定手段は、前記ターゲット位置予測手段により予測される前記ターゲットの位置から今回の到来波推定の際に並ぶと予測される前記ターゲットの数を求め、この数と前回の到来波推定の際に並んでいた前記ターゲットの数との差が、前記前回の到来波推定の際に推定された到来波の数と前記今回の到来波推定の際に推定された到来波の数との差と一致していない場合には、前記前回の到来波推定の際に推定された到来波の数に対して、前記今回の到来波推定の際に並ぶと予測される前記ターゲットの数と前回の到来波推定の際に並んでいた前記ターゲットの数との差を加算もしくは減算することで、前記到来波の数の推定を行うようになっていることを特徴とする請求項1に記載のレーダ装置。   The arrival wave number estimation means obtains the number of targets predicted to be arranged in the current arrival wave estimation from the target position predicted by the target position prediction means, and this number and the previous arrival wave estimation The difference between the number of targets lined up at the time is the difference between the number of incoming waves estimated during the previous arrival wave estimation and the number of arrival waves estimated during the current arrival wave estimation. Is not equal to the number of incoming waves estimated during the previous arrival wave estimation, the number of targets predicted to be aligned during the current arrival wave estimation and the previous arrival wave estimation The radar according to claim 1, wherein the number of the incoming waves is estimated by adding or subtracting a difference from the number of the targets arranged in the case of the incoming wave estimation. apparatus. 送信アンテナ(4)を通じて送信信号を送信すると共に、互いに異なる位置に配設された複数のアンテナ(5a)を介して複数の方向から到来する複数の到来波を受信し、各アンテナ毎の受信信号から、前記送信信号を反射したターゲットからの到来波の方向の推定を行うレーダ装置において、
前記送信信号と前記受信信号の差の周波数成分であるビート信号を生成するビート信号生成手段(6)と、
前記ビート信号に基づいて自己相関行列(Rxx)を算出する自己相関行列算出手段(110)と、
前記自己相関行列算出手段で算出された自己相関行列から固有値(λ)を算出する固有値算出手段(120)と、
前記固有値算出手段で算出された固有値(λ)に基づいて前記到来波の数の推定を行う到来波数推定手段(130)と、
前記ビート信号に基づいて、前記ターゲットの距離および相対速度を求め、それらを履歴情報として記憶する履歴情報記憶手段とを備え、
前記到来波数推定手段は、前記固有値(λ)のスレッショルド(λTH)を設定するスレッショルド設定手段(200、210)を有し、該スレッショルド設定手段によって設定された前記スレッショルド(λTH)に基づいて前記固有値(λ)を信号空間の固有値(λ1〜λL)と雑音空間の固有値(λL+1〜λLK)に分別し、前記到来波の数の推定を行うと共に、前記履歴情報記憶手段に記憶された前記履歴情報から、今回の到来波推定の際における前記ターゲットの位置を予測するターゲット位置予測手段(400)を有し、該ターゲット位置予測手段により予測される前記ターゲットの位置から前記ターゲットの複数が並ぶか否かを予測し、その予測結果に基づいて前記スレッショルド設定手段が設定する前記スレッショルド(λ)の値を調整するようになっていることを特徴とするレーダ装置。
A transmission signal is transmitted through the transmission antenna (4), and a plurality of incoming waves arriving from a plurality of directions are received via a plurality of antennas (5a) arranged at different positions, and a reception signal for each antenna is received. From the radar apparatus for estimating the direction of the incoming wave from the target reflecting the transmission signal,
Beat signal generating means (6) for generating a beat signal which is a frequency component of a difference between the transmission signal and the reception signal;
Autocorrelation matrix calculating means (110) for calculating an autocorrelation matrix (Rxx) based on the beat signal;
Eigenvalue calculating means (120) for calculating the eigenvalue (λ) from the autocorrelation matrix calculated by the autocorrelation matrix calculating means;
Arriving wave number estimating means (130) for estimating the number of incoming waves based on the eigenvalue (λ) calculated by the eigenvalue calculating means;
Based on the beat signal, the distance and relative speed of the target are obtained, and history information storage means for storing them as history information,
The arrival wave number estimating means includes threshold setting means (200, 210) for setting a threshold (λTH) of the eigenvalue (λ), and the eigenvalue is based on the threshold (λTH) set by the threshold setting means. (Λ) is divided into eigenvalues (λ1 to λL) of the signal space and eigenvalues (λL + 1 to λLK) of the noise space to estimate the number of incoming waves, and the history information stored in the history information storage means To target position prediction means (400) for predicting the position of the target at the time of current arrival wave estimation, and whether or not a plurality of the targets are arranged from the target positions predicted by the target position prediction means. Of the threshold (λ) set by the threshold setting means based on the prediction result. Radar apparatus, characterized by being adapted to adjust.
JP2010203174A 2004-07-06 2010-09-10 Radar equipment Active JP4780240B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2010203174A JP4780240B2 (en) 2004-07-06 2010-09-10 Radar equipment

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2004199327 2004-07-06
JP2004199327 2004-07-06
JP2010203174A JP4780240B2 (en) 2004-07-06 2010-09-10 Radar equipment

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2005155627A Division JP4867200B2 (en) 2004-07-06 2005-05-27 Radar equipment

Publications (2)

Publication Number Publication Date
JP2010271337A true JP2010271337A (en) 2010-12-02
JP4780240B2 JP4780240B2 (en) 2011-09-28

Family

ID=42015556

Family Applications (3)

Application Number Title Priority Date Filing Date
JP2009270401A Active JP4780229B2 (en) 2004-07-06 2009-11-27 Radar equipment
JP2010203174A Active JP4780240B2 (en) 2004-07-06 2010-09-10 Radar equipment
JP2010203175A Pending JP2011007810A (en) 2004-07-06 2010-09-10 Radar system

Family Applications Before (1)

Application Number Title Priority Date Filing Date
JP2009270401A Active JP4780229B2 (en) 2004-07-06 2009-11-27 Radar equipment

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2010203175A Pending JP2011007810A (en) 2004-07-06 2010-09-10 Radar system

Country Status (1)

Country Link
JP (3) JP4780229B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015029339A1 (en) * 2013-08-29 2015-03-05 パナソニックIpマネジメント株式会社 Radar system and target detection method
US11125871B2 (en) 2018-02-20 2021-09-21 Denso Corporation Azimuth estimation device and method

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5628732B2 (en) * 2011-04-04 2014-11-19 富士通テン株式会社 Arithmetic apparatus for radar apparatus, radar apparatus, arithmetic method and program for radar apparatus
KR101231799B1 (en) 2011-12-02 2013-02-08 충남대학교산학협력단 System for calculating self position information
JP5978754B2 (en) * 2012-05-16 2016-08-24 株式会社デンソー Radar equipment
JP6997589B2 (en) * 2017-10-24 2022-01-17 三菱重工業株式会社 Direction estimation device, direction estimation method and program

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04348292A (en) * 1991-05-27 1992-12-03 Fujitsu Ten Ltd Measuring apparatus of distance between objects
JP2000121716A (en) * 1998-10-13 2000-04-28 Anritsu Corp Radio wave propagation estimating equipment
JP2001091639A (en) * 1999-09-27 2001-04-06 Toyota Central Res & Dev Lab Inc Fm-cw radar device
JP2001281326A (en) * 2000-03-29 2001-10-10 Toyota Central Res & Dev Lab Inc Radar signal processing circuit
JP4577742B2 (en) * 2000-11-14 2010-11-10 住友電気工業株式会社 Radio wave arrival direction estimation device
JP3784640B2 (en) * 2000-12-01 2006-06-14 三菱電機株式会社 Wave source detection apparatus and wave source detection method
JP3777576B2 (en) * 2002-08-27 2006-05-24 学校法人金沢工業大学 Position estimation method and apparatus
JP2004108942A (en) * 2002-09-18 2004-04-08 Toshiba Corp Determination method for number of estimable incoming directions and number of separable signal
JP4220821B2 (en) * 2003-04-03 2009-02-04 三菱電機株式会社 Sensor signal processing system

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015029339A1 (en) * 2013-08-29 2015-03-05 パナソニックIpマネジメント株式会社 Radar system and target detection method
US11125871B2 (en) 2018-02-20 2021-09-21 Denso Corporation Azimuth estimation device and method

Also Published As

Publication number Publication date
JP4780229B2 (en) 2011-09-28
JP4780240B2 (en) 2011-09-28
JP2010044095A (en) 2010-02-25
JP2011007810A (en) 2011-01-13

Similar Documents

Publication Publication Date Title
JP4867200B2 (en) Radar equipment
JP6313981B2 (en) Radar apparatus, vehicle control system, and signal processing method
JP6304257B2 (en) Radar apparatus, vehicle, and moving body speed detection method
US7912680B2 (en) Direction-of-arrival estimation apparatus
JP4780240B2 (en) Radar equipment
US9310470B2 (en) Radar apparatus and signal processing method
JP5078637B2 (en) Radar apparatus and target detection method
JP6520203B2 (en) Mounting angle error detection method and device, vehicle-mounted radar device
US9874634B1 (en) Radar apparatus
JP4829517B2 (en) Radar signal processing device
JP6413457B2 (en) Precipitation determination device
JP2001242242A (en) Millimeter-wave radar device with function for improving detecting performance
US9285467B2 (en) Radar apparatus, vehicle control system, and signal processing method
US9383440B2 (en) Radar apparatus and signal processing method
JP7014041B2 (en) Radar device
JP2003167047A (en) Miss-pairing determination and signal processing method for rm-cw radar
US9689975B2 (en) Radar apparatus
US20190293759A1 (en) Object detection apparatus
EP2834665A1 (en) Phase monopulse radar system and target detecting method
JP2019200083A (en) Radar device and signal processing method
US20210405177A1 (en) Object tracking device
JP5462452B2 (en) Signal processing apparatus and radar apparatus
JP7188894B2 (en) Radar device and signal processing method
JP2009128016A (en) Radar system, its control method, and radar control apparatus
US11892557B2 (en) Radar device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20100913

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

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20110620

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140715

Year of fee payment: 3

R151 Written notification of patent or utility model registration

Ref document number: 4780240

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140715

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250