JP2010271337A - Radar device - Google Patents
Radar device Download PDFInfo
- 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
Links
Images
Abstract
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,
この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.
ここで、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)
WHRxxW=λ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.
このため、自己相関行列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
これらに示した手法はいずれの手法も、データ収集の回数(スナップショット)が多ければ多いほど、また受信信号の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
このように、前回のサイクルの際に得られたターゲットの履歴情報に基づいて、前回のサイクルの際に並走していなかったターゲット同士が今回のサイクルの際に並走するか否か、もしくは、前回のサイクルの際に並走していたターゲット同士が今回のサイクルの際に並走しなくなるか否かを予測し、その予測に基づいて到来波数を修正するようにしている。これにより、より正確に到来波数の推定を行うことが可能となる。 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
また、請求項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.
以下、本発明の実施形態について図に基づいて説明する。なお、以下の各実施形態相互において、互いに同一もしくは均等である部分には、図中、同一符号を付してある。 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
変調信号生成部1は、例えば、時間に対して直線的に周波数が増減するような変調を行わせる所定周期の三角波状の変調信号を生成するものである。この変調信号生成部1が生成する変調信号が電圧制御発信器2に出力されるようになっている。
For example, the modulation
電圧制御発振器2は、送信信号となるミリ波帯の高周波信号を生成するものである。この電圧制御発振器2は、変調信号生成部1から出力された変調信号に基づいて高周波信号の周波数を変調信号に応じて制御し、分配器3に出力するようになっている。
The voltage controlled
分配器3は、電圧制御発振器2から出力された送信信号を電力分配し、ローカル信号を生成するもので、この分配器3からのローカル信号がミキサ6に入力されるようになっている。
The
送信アンテナ4は、分配器3を経て、電圧制御発振器2から出力された送信信号で示されるミリ波帯の電波を車両の進行方向、例えば前方に向けて出力するものである。
The
受信部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
ミキサ6は、受信部5から伝えられる受信データxi(t)の結合と、分配器3から伝えられるローカル信号とを混合し、これらの信号の差の周波数成分であるビート信号を生成するものである。このとき生成されるビート信号の周波数がビート周波数と呼ばれるもので、送信信号の周波数が増加する時のビート周波数を上り変調時のビート周波数、送信信号の周波数が減少する時のビート周波数を下り変調時のビート周波数と呼び、FMCW方式による目標物体の距離および相対速度の演算に用いられる。
The
A/Dコンバータ7は、ミキサ6から出力されたアナログ値として示されるビート信号をデジタル値に変換するためのものである。
The A / D converter 7 is for converting the beat signal shown as an analog value output from the
信号処理部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
続く、ステップ110では、受信データxi(t)に基づいて自己相関行列Rxxが算出され、ステップ120では、ステップ110で算出された自己相関行列Rxxとなるときの固有値λ(λ1〜λK)および固有ベクトルEn、Esが算出される。なお、ここでいう固有ベクトルEnとは、それぞれ上述した雑音空間の固有値λを算出する上で掛けたウェイトWに相当するもので、固有ベクトルEsとは、信号空間の固有値λを算出する上で掛けたウェイトWに相当するものである。これら自己相関行列Rxxの算出手法、固有値λおよび固有ベクトルEn、Esの算出手法に関しては従来と同様であるため、ここでは説明を省略する。このように、信号処理部8における自己相関行列Rxxを算出する部分が自己相関行列算出手段に相当するもので、固有値を算出する部分が固有値算出手段に相当するものである。
In
ステップ130では、到来波数推定処理が実行される。この到来波数推定処理では、図3のフローチャートに示される各処理が実行される。信号処理部8のうち、この到来波数推定処理を実行する部分が到来波数推定手段に相当するものである。
In
到来波数推定処理が実行されると、まず、ステップ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
続くステップ210では、信号処理部8に入力されたビート信号のビート周波数に応じたスレッショルドλTHがマップに記憶されたビート周波数とスレッショルドλTH1〜λTH256との相関関係から読み取られ、そのときのビート周波数に応じたスレッショルドλTHとして設定される。以上のようなスレッショルドの設定を行う部分がスレッショルド設定手段に相当する。
In the
そして、ステップ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
そして、ステップ260において、信号空間の固有値λ1〜λLの数から到来波数Lを推定する。これにより、雑音などの不要信号を考慮にいれた正確な到来波数Lの推定を行うことが可能となる。
In
このように到来波数の推定が行われると、ここで求められた到来波数に基づき、ステップ140において、MUSICスペクトルが算出される。このMUSICスペクトルの算出手法に関しても従来と同様である。
When the arrival wave number is estimated in this way, a MUSIC spectrum is calculated in
続いて、ステップ150において、ステップ140で求められた本来の到来波のピーク波に基づいて到来波の方位が算出される。すなわち、到来波のピーク波の分布を見れば、到来波のピーク波に対応する受信アンテナの持つビームが指向する方位を求めることが可能となる。この方位が目標物体の存在する方位として算出される。
Subsequently, in
ステップ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
図5は、送信アンテナ4から出力されるミリ波帯の電波の出力方位エリアとそのエリア内において車両が並走できる最大数との関係を示した模式図である。
FIG. 5 is a schematic diagram showing a relationship between an output azimuth area of millimeter wave radio waves output from the
この図に示されるように、出力方位エリアの横幅は自車両からの距離から一義的に決まるため、その出力方位エリア内で並走可能な車両(ターゲット)の数が必然的に決まってくる。例えば、出力方位エリアの角度が車両前方を中心とした左右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
続く、ステップ310において、読み出したマップに基づいて、最大値Lmaxが設定される。具体的には、レーダ装置では、図2に示したように受信パワーに基づいてFMCW方式による距離および相対速度の演算が行われている。このため、ここでの演算結果をターゲットとなる目標物体のエリア情報(距離情報)として用い、そのエリア情報に応じた最大値Lmaxが設定される。以上のような、到来波の数の最大値を設定する部分が到来波最大値設定手段に相当する。
In
そして、ステップ320において、ステップ240のように分別された信号空間の固有値λ1〜λLから求められる到来波数Lが最大値Lmax以下になっているか否かが判定される。ここで肯定判定された場合には、到来波数Lが物理的に並走しうる車両の数以下になっており、矛盾していないものとして、ステップ260において、そのときの到来波数Lがそのまま到来波数Lとして出力される。
In
逆に、ここで否定判定された場合には、到来波数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
以上説明したように、本実施形態では、出力方位エリア内に物理的に並走しうる車両の数を到来波数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
そして、ステップ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
以上説明したように、本実施形態では、前回のサイクルの時に得られたターゲットの履歴情報に基づいて、前回のサイクルの時に並走していなかったターゲット同士が今回のサイクルの時に並走するかを予測し、その予測に基づいて到来波数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
そして、ステップ510において、前回の到来波数に対して1つ加算した値が設定されているかが判定される。このステップで肯定判定された場合には、ステップ520に進み、スレッショルドλTHをαだけ下げたのちステップ230に進み、否定判定された場合には、そのままステップ230に進む。
In
ここでいうαは、信号空間の固有値λ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
以上説明したように、本実施形態も、前回のサイクルの時に得られたターゲットの履歴情報に基づいて、前回のサイクルの時に並走していなかったターゲット同士が今回のサイクルの時に並走するかを予測し、その予測に基づいてスレッショルドλ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
この後、第3実施形態で示したステップ300〜330およびステップ400から420およびステップ260の処理が実行されることで、到来波数Lが推定される。そして、このようにして求められた値が仮の到来波数Lとされる。これにより、到来波数推定処理が終了となる。
Thereafter, the processing of
そして、この到来波数推定処理が完了すると、図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
図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
まず、ステップ700では、ステップ140で求められたMUSICスペクトルから確認される到来波のピーク波の方位(以下、ピークサーチ方位という)D1〜DLと受信パワーP1〜PLが検出される。なお、ここでは、上述したステップ260において、仮の到来波数Lとされていることから、ピーク波の数もそれと同じLとされる。なお、ピークサーチ方位D1〜DLと受信パワーP1〜PLの算出手法に関しては、上述したステップ150、160と同様である。
First, in
次に、ステップ710では、ビート周波数ごとの受信パワーのスレッショルドPTH1〜PTH256の読み出しが行われる。ここでいうビート周波数ごとの受信パワーのスレッショルドPTH1〜PTH256とは、ビート周波数[bin]と受信パワーとの相関関係から求められるものであり、例えば予め行った実験結果から求められる。
Next, in
続くステップ720では、信号処理部8に入力されたビート信号のビート周波数に応じた受信パワーのスレッショルドPTHがマップに記憶されたビート周波数と受信パワーのスレッショルドPTH1〜PTH256との相関関係から読み取られ、そのときのビート周波数に応じた受信パワーのスレッショルドPTHとして設定される。以上のようなスレッショルドの設定を行う部分がスレッショルド設定手段に相当する。
In the
そして、ステップ730において、ステップ700で求められた受信パワーP1〜PLがステップ720で設定されたビート周波数に応じた受信パワーのスレッショルドPTHを超えているか否かが判定され、受信パワーP1〜PLのうちスレッショルドPTHを超えていたもののみが実際に存在する到来波であるものとされる。そして、その実際の到来波の方位と受信パワーとの情報を用いて、目標物体の存在する方位および目標物体からの距離、相対速度が求められる。
In
以上説明したように、本実施形態では、到来波数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
まず、ステップ700では、上記ステップ140で求められたMUSICスペクトルからピークサーチ方位D1〜DLが検出される。
First, in
次に、ステップ740では、各ピークサーチ方位D1〜DLでの受信パワーP1〜PLが算出される。具体的には、受信パワーP1〜PLは次式により求められ、Sの対角成分から抽出される(菊間信良、アダプティブアンテナ技術、オーム社、P.137−P.141参照)。
Next, in
(数5)
S=(AHA)-1AH(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
以上説明したように、本実施形態でも、第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
Claims (3)
前記送信信号と前記受信信号の差の周波数成分であるビート信号を生成するビート信号手段(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.
前記送信信号と前記受信信号の差の周波数成分であるビート信号を生成するビート信号生成手段(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.
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)
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)
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)
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 |
-
2009
- 2009-11-27 JP JP2009270401A patent/JP4780229B2/en active Active
-
2010
- 2010-09-10 JP JP2010203174A patent/JP4780240B2/en active Active
- 2010-09-10 JP JP2010203175A patent/JP2011007810A/en active Pending
Cited By (2)
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 |