JP2019215275A - Target detection device, and target detection method - Google Patents

Target detection device, and target detection method Download PDF

Info

Publication number
JP2019215275A
JP2019215275A JP2018113047A JP2018113047A JP2019215275A JP 2019215275 A JP2019215275 A JP 2019215275A JP 2018113047 A JP2018113047 A JP 2018113047A JP 2018113047 A JP2018113047 A JP 2018113047A JP 2019215275 A JP2019215275 A JP 2019215275A
Authority
JP
Japan
Prior art keywords
target
unit
data
instantaneous data
vehicle
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
JP2018113047A
Other languages
Japanese (ja)
Other versions
JP7164327B2 (en
Inventor
聖也 藤津
Seiya Fujitsu
聖也 藤津
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 Ten Ltd
Original Assignee
Denso Ten Ltd
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 Ten Ltd filed Critical Denso Ten Ltd
Priority to JP2018113047A priority Critical patent/JP7164327B2/en
Priority to US16/437,044 priority patent/US11237260B2/en
Publication of JP2019215275A publication Critical patent/JP2019215275A/en
Application granted granted Critical
Publication of JP7164327B2 publication Critical patent/JP7164327B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

To provide a target detection device and target detection method that can upgrade responsiveness with respect to a target.SOLUTION: A target detection device 1 comprises: a generation unit 32a; a calculation unit; and an estimation unit 323b. The generation unit is configured to generate instant data including a relative speed toward an own vehicle of a target as to each of a plurality of reflection points where a transmitted radio wave is reflected upon a target. The calculation unit is configured to calculate a ground speed on the basis of the relative speed of the instant data generated by the generation unit and a speed of the own vehicle. The estimation unit is configured to make origins in a ground vector toward the own vehicle corresponding to each of a plurality of instant data uniform on the basis of the ground speed calculated by the calculation unit, and estimate a movement direction of the target from an intersection of a perpendicular line with respect to the ground vector.SELECTED DRAWING: Figure 2

Description

本発明は、物標検出装置および物標検出方法に関する。   The present invention relates to a target detection device and a target detection method.

従来、例えば車両の周囲へ電波を送信し、送信した電波が物標で反射した反射波に基づいて物標の自車両への向きの相対速度等を含む瞬時データを検出する物標検出装置が知られている(例えば、特許文献1参照)。   Conventionally, for example, a target detection device that transmits a radio wave around a vehicle and detects instantaneous data including a relative speed of a direction of the target toward the own vehicle based on a reflected wave of the transmitted radio wave reflected by the target is known. It is known (for example, see Patent Document 1).

特開2000−193745号公報JP 2000-193745 A

しかしながら、上記した物標検出装置において検出される相対速度は、自車両への向きの速度成分のみであり、物標の実際の移動向きの速度成分を得るためには、初回検知から数スキャン分の瞬時データが必要となる。従って、従来の物標検出装置では、数スキャン後に物標の実際の移動向きが分かるため、物標に対する応答性を上げられないという課題があった。   However, the relative speed detected by the target detection device described above is only the speed component in the direction to the own vehicle, and in order to obtain the speed component in the actual moving direction of the target, it takes several scans from the initial detection. Instantaneous data is required. Therefore, the conventional target detection device has a problem that the responsiveness to the target cannot be improved because the actual moving direction of the target is known after several scans.

本発明は、上記に鑑みてなされたものであって、物標に対する応答性を向上させることができる物標検出装置および物標検出方法を提供することを目的とする。   The present invention has been made in view of the above, and an object of the present invention is to provide a target detection device and a target detection method capable of improving responsiveness to a target.

上述した課題を解決し、目的を達成するために、本発明に係る物標検出装置は、生成部と、算出部と、推定部とを備える。前記生成部は、送信した電波が物標で反射した複数の反射点それぞれについて、前記物標の自車両への向きの相対速度を含む瞬時データを生成する。前記算出部は、前記生成部によって生成された前記瞬時データの前記相対速度と前記自車両の速度とに基づいて、対地速度を算出する。前記推定部は、前記算出部によって算出された前記対地速度に基づき、複数の前記瞬時データそれぞれに対応する前記自車両への向きの対地ベクトルにおける起点を揃え、当該対地ベクトルに対する垂線の交点から前記物標の移動向きを推定する。   In order to solve the above-described problem and achieve the object, a target detection device according to the present invention includes a generation unit, a calculation unit, and an estimation unit. The generation unit generates, for each of a plurality of reflection points at which the transmitted radio wave is reflected by the target, instantaneous data including a relative speed of a direction of the target toward the own vehicle. The calculation unit calculates a ground speed based on the relative speed of the instantaneous data generated by the generation unit and the speed of the host vehicle. The estimating unit, based on the ground speed calculated by the calculating unit, aligns the starting point in the ground vector of the direction to the own vehicle corresponding to each of the plurality of instantaneous data, from the intersection of the perpendicular to the ground vector. Estimate the moving direction of the target.

本発明によれば、物標に対する応答性を向上させることができる。   According to the present invention, responsiveness to a target can be improved.

図1Aは、実施形態に係る物標検出方法の概要を示す図である。FIG. 1A is a diagram illustrating an outline of a target detection method according to the embodiment. 図1Bは、実施形態に係る物標検出方法の概要を示す図である。FIG. 1B is a diagram illustrating an outline of a target detection method according to the embodiment. 図1Cは、実施形態に係る物標検出方法の概要を示す図である。FIG. 1C is a diagram illustrating an outline of the target detection method according to the embodiment. 図2は、実施形態に係るレーダ装置の構成を示すブロック図である。FIG. 2 is a block diagram illustrating a configuration of the radar device according to the embodiment. 図3は、信号処理部の前段処理から信号処理部におけるピーク抽出処理までの処理説明図である。FIG. 3 is an explanatory diagram of processing from the pre-processing of the signal processing unit to the peak extraction processing in the signal processing unit. 図4Aは、角度推定処理の処理説明図である。FIG. 4A is an explanatory diagram of the angle estimation process. 図4Bは、ペアリング処理の処理説明図(その1)である。FIG. 4B is an explanatory diagram (part 1) of the pairing process. 図4Cは、ペアリング処理の処理説明図(その2)である。FIG. 4C is an explanatory diagram (part 2) of the pairing process. 図5は、対地速度に基づく予測処理を示す図である。FIG. 5 is a diagram illustrating a prediction process based on a ground speed. 図6は、推定部の機能ブロック図である。FIG. 6 is a functional block diagram of the estimation unit. 図7は、バッファ更新部の更新処理を示す図である。FIG. 7 is a diagram illustrating an update process of the buffer update unit. 図8は、対地速度算出部による対地速度の算出処理を示す図である。FIG. 8 is a diagram illustrating a process of calculating a ground speed by a ground speed calculation unit. 図9Aは、交点が3つの場合における移動ベクトルの作成方法を示す図である。FIG. 9A is a diagram showing a method of creating a movement vector when there are three intersections. 図9Bは、交点が6つの場合における移動ベクトルの作成方法を示す図である。FIG. 9B is a diagram showing a method of creating a movement vector when there are six intersections. 図10は、移動ベクトルに基づく重み付け処理を示す図である。FIG. 10 is a diagram illustrating a weighting process based on a movement vector. 図11は、実施形態に係る物標検出装置が実行する処理の処理手順を示すフローチャートである。FIG. 11 is a flowchart illustrating a processing procedure of a process executed by the target detection device according to the embodiment. 図12は、実施形態に係る推定部が実行する推定処理の処理手順を示すフローチャートである。FIG. 12 is a flowchart illustrating a processing procedure of an estimation process performed by the estimation unit according to the embodiment. 図13は、移動ベクトルの検証処理を示す図である。FIG. 13 is a diagram illustrating the verification processing of the movement vector.

以下、添付図面を参照して、本願の開示する物標検出装置および物標検出方法の実施形態を詳細に説明する。なお、この実施形態により本発明が限定されるものではない。なお、以下では、物標検出装置がFM−CW(Frequency Modulated Continuous Wave)方式である場合を例に挙げて説明するが、物標検出装置は、例えばFCM(Fast-Chirp Modulation)方式といった他の方式であってもよい。   Hereinafter, embodiments of a target detection device and a target detection method disclosed in the present application will be described in detail with reference to the accompanying drawings. The present invention is not limited by the embodiment. In the following, a case where the target detecting device is an FM-CW (Frequency Modulated Continuous Wave) method will be described as an example. However, the target detecting device may use another method such as an FCM (Fast-Chirp Modulation) method. The system may be used.

まず、図1A〜図1Cを用いて、実施形態に係る物標検出方法の概要について説明する。図1A〜図1Cは、実施形態に係る物標検出方法の概要を示す図である。図1Aでは、実施形態に係る物標検出装置1を搭載した自車両MCと、所定の移動ベクトルVで移動する他車両LCとを示している。移動ベクトルVとは、他車両LCの実際の移動向きおよびかかる移動向きへの対地速度を含むベクトルである。対地速度とは、自車両MCおよび他車両LCの相対速度から自車両MCの速度成分を除いた速度である。また、自車両MCは、所定の移動ベクトルMV(以下、自車ベクトルMV)で移動することとする。   First, an outline of a target detection method according to the embodiment will be described with reference to FIGS. 1A to 1C. 1A to 1C are diagrams illustrating an outline of a target detection method according to the embodiment. FIG. 1A shows a host vehicle MC on which the target detection device 1 according to the embodiment is mounted, and another vehicle LC that moves with a predetermined movement vector V. The movement vector V is a vector including the actual movement direction of the other vehicle LC and the ground speed in the movement direction. The ground speed is a speed obtained by removing the speed component of the own vehicle MC from the relative speed of the own vehicle MC and the other vehicle LC. In addition, the own vehicle MC moves with a predetermined movement vector MV (hereinafter, own vehicle vector MV).

図1Aに示すように、物標検出装置1は、例えば自車両MCのフロントグリル内等に搭載され、自車両MCの進行方向に存在する物標(例えば、他車両LC等)を検出する。なお、物標検出装置1の搭載箇所は、例えばフロントガラスやリアグリル、左右の側部(例えば、左右のドアミラー)等他の箇所に搭載されてもよい。   As shown in FIG. 1A, the target detection device 1 is mounted in, for example, a front grill of the own vehicle MC, and detects a target (for example, another vehicle LC or the like) existing in the traveling direction of the own vehicle MC. Note that the target detection device 1 may be mounted on another location such as a windshield, a rear grill, and left and right side portions (for example, left and right door mirrors).

また、図1Aに示すように、物標検出装置1は、自車両MCの周囲に送信した電波が他車両LCで反射した複数の反射点それぞれについて、瞬時データ100を生成する。瞬時データ100には、例えば、自車両MCへの向きの相対速度(ベクトルRV)等が含まれる。   Further, as shown in FIG. 1A, the target detection device 1 generates instantaneous data 100 for each of a plurality of reflection points at which radio waves transmitted around the own vehicle MC are reflected by other vehicles LC. The instantaneous data 100 includes, for example, the relative speed (vector RV) of the direction to the host vehicle MC.

ここで、従来の物標検出装置について説明する。従来の物標検出装置では、1スキャンで検出される瞬時データは、自車両への向きの相対速度のみ(例えば、図1Aに示すベクトルRV)であり、上記した移動ベクトルを得ることができない。つまり、1スキャンでは、物標の実際の移動向きを得ることができず、物標の実際の移動向きを得るためには、初回検知から数スキャン分の瞬時データの変化から求めなければならない。従って、従来の物標検出装置では、数スキャン後に物標の実際の移動向きが分かるため、物標に対する応答性を上げられないという課題があった。   Here, a conventional target detection device will be described. In the conventional target detection device, the instantaneous data detected in one scan is only the relative speed in the direction toward the host vehicle (for example, the vector RV shown in FIG. 1A), and the above-described movement vector cannot be obtained. That is, in one scan, the actual moving direction of the target cannot be obtained, and in order to obtain the actual moving direction of the target, it is necessary to obtain the actual moving direction of the target from a change in instantaneous data for several scans from the first detection. Therefore, the conventional target detection device has a problem that the responsiveness to the target cannot be improved because the actual moving direction of the target is known after several scans.

そこで、実施形態に係る物標検出方法では、1スキャンで物標の実際の移動向きを検出可能とした。具体的には、図1Bに示すように、実施形態に係る物標検出方法では、まず、生成した複数の瞬時データ100のベクトルRV(以下、相対ベクトルRV)における相対速度から自車両MCの速度成分を除いた自車両MCへの向きの対地速度を示すベクトルGV(以下、対地ベクトルGV)を算出する。具体的には、対地ベクトルGVは、移動ベクトルVが自車両MCと瞬時データ100とを通る直線に射影されたベクトルである。また、自車両MCの速度成分は、自車ベクトルMVを上記した直線に射影したベクトルであるが、かかる点については図8で後述する。   Therefore, in the target detecting method according to the embodiment, the actual moving direction of the target can be detected in one scan. Specifically, as shown in FIG. 1B, in the target detection method according to the embodiment, first, the speed of the host vehicle MC is calculated from the relative speed in the vector RV (hereinafter, relative vector RV) of the plurality of generated instantaneous data 100. A vector GV (hereinafter, ground vector GV) indicating the ground speed in the direction toward the own vehicle MC excluding the components is calculated. Specifically, the ground vector GV is a vector obtained by projecting the movement vector V onto a straight line passing through the own vehicle MC and the instantaneous data 100. The speed component of the host vehicle MC is a vector obtained by projecting the host vehicle vector MV onto the above-described straight line. This point will be described later with reference to FIG.

つづいて、図1Cに示すように、実施形態に係る物標検出方法では、複数の瞬時データ100における対地速度に基づき、複数の瞬時データ100それぞれに対応する自車両MCへの向きの対地ベクトルGVにおける起点を揃え、対地ベクトルGVの終点から伸ばした垂線VLの交点CPから物標である他車両LCの移動向きを推定する。   Subsequently, as shown in FIG. 1C, in the target detection method according to the embodiment, based on the ground speed in the plurality of instantaneous data 100, the ground vector GV of the direction to the own vehicle MC corresponding to each of the plurality of instantaneous data 100. Are aligned, and the moving direction of the other vehicle LC as the target is estimated from the intersection CP of the perpendicular VL extending from the end point of the ground vector GV.

図1Cに示す例では、2つの瞬時データ100に対応する対地ベクトルGVの起点を揃えた2次元平面を示している。なお、図1Cにおいて、縦軸は、自車両MCに対して左右方向への対地速度(VY)を示し、横軸は、自車両MCに対して前後方向への対地速度(VX)を示す。   The example illustrated in FIG. 1C illustrates a two-dimensional plane in which the starting points of the ground vectors GV corresponding to the two instantaneous data 100 are aligned. In FIG. 1C, the vertical axis indicates the ground speed (VY) in the left-right direction with respect to the host vehicle MC, and the horizontal axis indicates the ground speed (VX) in the front-rear direction with respect to the host vehicle MC.

図1Cに示すように、実施形態に係る物標検出方法では、まず、起点が揃った2つの対地ベクトルGVの終点を通る垂線VLを引き、2つの垂線VLの交点CPを求める。そして、実施形態に係る物標検出方法では、交点CPが1つであった場合、対地ベクトルGVの起点を移動ベクトルVの起点とし、交点CPを移動ベクトルVの終点とする。すなわち、対地ベクトルGVの起点から交点CPへの向きが物標の実際の移動向きを示し、対地ベクトルGVの起点と交点CPとの間の距離が物標の移動向きへの対地速度を示す。   As shown in FIG. 1C, in the target detection method according to the embodiment, first, a perpendicular VL passing through the end points of two ground vectors GV whose starting points are aligned is drawn, and an intersection CP of the two perpendiculars VL is obtained. In the target detection method according to the embodiment, when there is one intersection CP, the starting point of the ground vector GV is set as the starting point of the movement vector V, and the intersection CP is set as the end point of the movement vector V. That is, the direction from the starting point of the ground vector GV to the intersection CP indicates the actual moving direction of the target, and the distance between the starting point of the ground vector GV and the intersection CP indicates the ground speed in the moving direction of the target.

このように、実施形態に係る物標検出方法によれば、複数の瞬時データ100に基づいて1スキャンで物標の移動向きを検出できるため、数スキャン必要とした従来に比べて、物標に対する応答性を向上させることができる。   As described above, according to the target detection method according to the embodiment, the moving direction of the target can be detected in one scan based on the plurality of instantaneous data 100. Responsiveness can be improved.

さらに、実施形態に係る物標検出方法では、瞬時データ100の対地速度に基づいて移動ベクトルVを推定することで、相対速度と比べて、自車両MCの変化(速度や位置等)を考慮する必要がないため、移動ベクトルVの推定精度をさらに高めることができる。   Further, in the target detection method according to the embodiment, by estimating the movement vector V based on the ground speed of the instantaneous data 100, a change (speed, position, and the like) of the own vehicle MC is considered as compared with the relative speed. Since there is no need, the estimation accuracy of the movement vector V can be further improved.

なお、実施形態に係る物標検出方法では、物標の移動向きのみが必要であれば、移動向きへの対地速度である対地ベクトルGVの起点と交点CPとの間の距離を必ずしも求める必要はない。つまり、実施形態に係る物標検出方法では、少なくとも物標の移動向きを推定する。   In the target detection method according to the embodiment, if only the moving direction of the target is necessary, it is not always necessary to obtain the distance between the starting point of the ground vector GV, which is the ground speed in the moving direction, and the intersection CP. Absent. That is, in the target detection method according to the embodiment, at least the moving direction of the target is estimated.

また、実施形態に係る物標検出方法では、交点CPが3つ以上の場合に、かかる交点CPにより形成される三角形の内心を求めることで、物標の移動向きを推定するが、かかる点については後述する。   In the target detection method according to the embodiment, when the number of intersections CP is three or more, the moving direction of the target is estimated by obtaining the inner center of a triangle formed by the intersections CP. Will be described later.

次に、図2を参照して、実施形態に係る物標検出装置1の構成について詳細に説明する。図2は、実施形態に係る物標検出装置1の構成を示すブロック図である。なお、図2では、本実施形態の特徴を説明するために必要な構成要素のみを機能ブロックで表しており、一般的な構成要素についての記載を省略している。   Next, the configuration of the target detection device 1 according to the embodiment will be described in detail with reference to FIG. FIG. 2 is a block diagram illustrating a configuration of the target detection device 1 according to the embodiment. In FIG. 2, only the components necessary for describing the features of the present embodiment are represented by functional blocks, and descriptions of general components are omitted.

換言すれば、図2に図示される各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。例えば、各機能ブロックの分散・統合の具体的形態は図示のものに限られず、その全部または一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的または物理的に分散・統合して構成することが可能である。   In other words, each component illustrated in FIG. 2 is a functional concept and does not necessarily need to be physically configured as illustrated. For example, the specific form of distribution / integration of each functional block is not limited to the illustrated one, and all or a part of the functional blocks may be functionally or physically distributed in arbitrary units according to various loads and usage conditions. -It is possible to integrate and configure.

図2に示すように、物標検出装置1は、送信部10と、受信部20と、処理部30とを備える。物標検出装置1は、自車両MCの挙動を制御する車両制御装置2に接続される。   As illustrated in FIG. 2, the target detection device 1 includes a transmission unit 10, a reception unit 20, and a processing unit 30. The target detection device 1 is connected to a vehicle control device 2 that controls the behavior of the host vehicle MC.

かかる車両制御装置2は、物標検出装置1による物標の検出結果に基づいて、PCS(Pre-crash Safety System)やAEB(Advanced Emergency Braking System)などの車両制御を行う。   The vehicle control device 2 performs vehicle control such as a PCS (Pre-crash Safety System) and an AEB (Advanced Emergency Braking System) based on the detection result of the target by the target detection device 1.

送信部10は、信号生成部11と、発振器12と、送信アンテナ13とを備える。信号生成部11は、後述する送受信制御部31の制御により、三角波で周波数変調されたミリ波を送信するための変調信号を生成する。発振器12は、かかる信号生成部11によって生成された変調信号に基づいて送信信号を生成し、送信アンテナ13へ出力する。なお、図2に示すように、発振器12によって生成された送信信号は、後述するミキサ22に対しても分配される。   The transmission unit 10 includes a signal generation unit 11, an oscillator 12, and a transmission antenna 13. The signal generation unit 11 generates a modulation signal for transmitting a millimeter wave frequency-modulated by a triangular wave under the control of a transmission / reception control unit 31 described later. The oscillator 12 generates a transmission signal based on the modulation signal generated by the signal generation unit 11, and outputs the transmission signal to the transmission antenna 13. Note that, as shown in FIG. 2, the transmission signal generated by the oscillator 12 is also distributed to a mixer 22, which will be described later.

送信アンテナ13は、発振器12からの送信信号を送信波へ変換し、かかる送信波を自車両MCの外部へ出力する。送信アンテナ13が出力する送信波は、三角波で周波数変調された連続波である。送信アンテナ13から自車両MCの外部、たとえば前方へ送信された送信波は、他車両LC等の物標で反射されて反射波となる。   The transmission antenna 13 converts a transmission signal from the oscillator 12 into a transmission wave, and outputs the transmission wave to the outside of the vehicle MC. The transmission wave output by the transmission antenna 13 is a continuous wave frequency-modulated by a triangular wave. A transmission wave transmitted from the transmission antenna 13 to the outside of the own vehicle MC, for example, forward, is reflected by a target such as another vehicle LC and becomes a reflected wave.

受信部20は、アレーアンテナを形成する複数の受信アンテナ21と、複数のミキサ22と、複数のA/D変換部23とを備える。ミキサ22およびA/D変換部23は、受信アンテナ21ごとに設けられる。   The receiving unit 20 includes a plurality of receiving antennas 21 forming an array antenna, a plurality of mixers 22, and a plurality of A / D converters 23. The mixer 22 and the A / D converter 23 are provided for each reception antenna 21.

各受信アンテナ21は、物標からの反射波を受信波として受信し、かかる受信波を受信信号へ変換してミキサ22へ出力する。なお、図2に示す受信アンテナ21の数は4つであるが、3つ以下または5つ以上であってもよい。   Each receiving antenna 21 receives a reflected wave from the target as a received wave, converts the received wave into a received signal, and outputs the signal to the mixer 22. Although the number of receiving antennas 21 shown in FIG. 2 is four, it may be three or less or five or more.

受信アンテナ21から出力された受信信号は、図示略の増幅器(たとえば、ローノイズアンプ)で増幅された後にミキサ22へ入力される。ミキサ22は、分配された送信信号と、受信アンテナ21から入力される受信信号との一部をミキシングし不要な信号成分を除去してビート信号を生成し、A/D変換部23へ出力する。   The reception signal output from the reception antenna 21 is input to the mixer 22 after being amplified by an amplifier (not shown) (not shown). The mixer 22 mixes a part of the distributed transmission signal and a part of the reception signal input from the reception antenna 21, removes unnecessary signal components, generates a beat signal, and outputs the beat signal to the A / D conversion unit 23. .

ビート信号は、送信信号の周波数(以下、「送信周波数」と記載する)と受信信号の周波数(以下、「受信周波数」と記載する)との差となるビート周波数を有する。ミキサ22で生成されたビート信号は、図示しない同期部によって受信アンテナ同士でタイミングを合わせた上でA/D変換部23でデジタル信号に変換された後に、処理部30へ出力される。   The beat signal has a beat frequency that is the difference between the frequency of the transmission signal (hereinafter, referred to as “transmission frequency”) and the frequency of the reception signal (hereinafter, referred to as “reception frequency”). The beat signal generated by the mixer 22 is converted into a digital signal by the A / D converter 23 after the timing is adjusted between the receiving antennas by a synchronization unit (not shown), and then output to the processing unit 30.

処理部30は、送受信制御部31と、信号処理部32と、記憶部33とを備える。信号処理部32は、生成部32aと、フィルタ処理部32bとを備える。   The processing unit 30 includes a transmission / reception control unit 31, a signal processing unit 32, and a storage unit 33. The signal processing unit 32 includes a generation unit 32a and a filter processing unit 32b.

記憶部33は、履歴データ33aを記憶する。履歴データ33aは、信号処理部32が実行する一連の信号処理における物標データの履歴や、瞬時データ100の履歴を含む情報である。また、履歴データ33aには、瞬時データ100の自車両MCへの向きの対地速度や自車両MCへの向きに相当する角度等といった情報も含まれるが、かかる点については後述する。   The storage unit 33 stores the history data 33a. The history data 33a is information including the history of the target data and the history of the instantaneous data 100 in a series of signal processing executed by the signal processing unit 32. The history data 33a also includes information such as the ground speed of the direction to the own vehicle MC of the instantaneous data 100 and an angle corresponding to the direction to the own vehicle MC. Such a point will be described later.

処理部30は、たとえば、CPU(Central Processing Unit)、記憶部33に対応するROM(Read Only Memory)やRAM(Random Access Memory)、レジスタ、その他の入出力ポートなどを含むマイクロコンピュータであり、物標検出装置1全体を制御する。   The processing unit 30 is, for example, a microcomputer including a CPU (Central Processing Unit), a ROM (Read Only Memory) and a RAM (Random Access Memory) corresponding to the storage unit 33, a register, and other input / output ports. The entire target detection device 1 is controlled.

かかるマイクロコンピュータのCPUがROMに記憶されたプログラムを読み出して実行することによって、送受信制御部31および信号処理部32として機能する。なお、送受信制御部31および信号処理部32は全部をASIC(Application Specific Integrated Circuit)やFPGA(Field Programmable Gate Array)などのハードウェアで構成することもできる。   When the CPU of the microcomputer reads out and executes the program stored in the ROM, the microcomputer functions as the transmission / reception control unit 31 and the signal processing unit 32. Note that the transmission / reception control unit 31 and the signal processing unit 32 may be entirely configured by hardware such as an ASIC (Application Specific Integrated Circuit) or an FPGA (Field Programmable Gate Array).

送受信制御部31は、信号生成部11を含む送信部10、および、受信部20を制御する。信号処理部32は、一連の信号処理を周期的に実行する。つづいて信号処理部32の各構成要素について説明する。   The transmission / reception control unit 31 controls the transmission unit 10 including the signal generation unit 11 and the reception unit 20. The signal processing unit 32 periodically executes a series of signal processing. Subsequently, each component of the signal processing unit 32 will be described.

生成部32aは、瞬時データ100を生成する。具体的には、生成部32aは、周波数解析処理と、ピーク抽出処理と、瞬時データ生成処理とを行うことで、瞬時データ100を生成する。   The generation unit 32a generates the instantaneous data 100. Specifically, the generation unit 32a generates the instantaneous data 100 by performing a frequency analysis process, a peak extraction process, and an instantaneous data generation process.

周波数解析処理では、各A/D変換部23から入力されるビート信号に対して高速フーリエ変換(FFT:Fast Fourier Transform)処理(以下、「FFT処理」と記載する)を行う。かかるFFT処理の結果は、ビート信号の周波数スペクトルであり、ビート信号の周波数ごと(周波数分解能に応じた周波数間隔で設定された周波数ビンごと)のパワー値(信号レベル)である。   In the frequency analysis processing, a fast Fourier transform (FFT: Fast Fourier Transform) processing (hereinafter, referred to as “FFT processing”) is performed on the beat signal input from each A / D conversion unit 23. The result of the FFT processing is a frequency spectrum of the beat signal, and a power value (signal level) for each frequency of the beat signal (for each frequency bin set at a frequency interval corresponding to the frequency resolution).

ピーク抽出処理では、周波数解析処理によるFFT処理の結果においてピークとなるピーク周波数を抽出する。なお、ピーク抽出処理では、後述するビート信号の「UP区間」および「DN区間」のそれぞれについてピーク周波数を抽出する。   In the peak extraction processing, a peak frequency that becomes a peak in the result of the FFT processing by the frequency analysis processing is extracted. In the peak extraction processing, a peak frequency is extracted for each of an “UP section” and a “DN section” of a beat signal described later.

瞬時データ生成処理では、ピーク抽出処理において抽出されたピーク周波数のそれぞれに対応する反射波の到来角度とそのパワー値を算出する角度推定処理を実行する。なお、角度推定処理の実行時点で、到来角度は、物標が存在すると推定される角度であることから、以下では「推定角度」と記載する場合がある。   In the instantaneous data generation process, an angle estimation process for calculating an arrival angle and a power value of a reflected wave corresponding to each of the peak frequencies extracted in the peak extraction process is executed. At the time of execution of the angle estimation process, the arrival angle is an angle at which the target is estimated to be present, and thus may be referred to as an “estimated angle” below.

また、瞬時データ生成処理では、算出した推定角度とパワー値との算出結果に基づいて「UP区間」および「DN区間」それぞれのピーク周波数の正しい組み合わせを判定するペアリング処理を実行する。   In the instantaneous data generation process, a pairing process of determining a correct combination of the peak frequencies of the “UP section” and the “DN section” based on the calculation result of the calculated estimated angle and the power value is executed.

また、瞬時データ生成処理では、判定した組み合わせ結果から各物標の自車両MCに対する距離および自車両MCへの向きの相対速度を算出する。また、瞬時データ処理では、算出した各物標の推定角度、距離および相対速度を、最新周期(最新スキャン)分の瞬時データ100としてフィルタ処理部32bへ出力するとともに、記憶部33の履歴データ33aとして記憶する。   In the instantaneous data generation processing, the distance of each target to the own vehicle MC and the relative speed of the direction to the own vehicle MC are calculated from the determined combination result. In the instantaneous data processing, the calculated estimated angle, distance, and relative speed of each target are output to the filter processing unit 32b as the instantaneous data 100 for the latest cycle (latest scan), and the history data 33a in the storage unit 33 is output. To be stored.

説明を分かりやすくするために、信号処理部32の前段処理から信号処理部32におけるここまでの処理の流れを図3〜図4Cに示す。図3は、信号処理部32の前段処理から生成部32aにおけるピーク抽出処理までの処理説明図である。   In order to make the description easy to understand, FIGS. 3 to 4C show the flow of processing from the pre-stage processing of the signal processing unit 32 to the processing in the signal processing unit 32. FIG. 3 is an explanatory diagram of a process from a pre-stage process of the signal processing unit 32 to a peak extraction process in the generation unit 32a.

また、図4Aは、角度推定処理の処理説明図である。また、図4Bおよび図4Cは、ペアリング処理の処理説明図(その1)および(その2)である。なお、図3は、2つの太い下向きの白色矢印で3つの領域に区切られている。以下では、かかる各領域を順に、上段、中段、下段と記載する。   FIG. 4A is an explanatory diagram of the angle estimation process. FIGS. 4B and 4C are explanatory diagrams (part 1) and (part 2) of the pairing process. FIG. 3 is divided into three regions by two thick downward white arrows. Hereinafter, each of these areas is described as an upper row, a middle row, and a lower row.

図3の上段に示すように、送信信号fs(t)は、送信アンテナ13から送信波として送出された後、物標において反射されて反射波として到来し、受信アンテナ21において受信信号fr(t)として受信される。   As shown in the upper part of FIG. 3, the transmission signal fs (t) is transmitted from the transmission antenna 13 as a transmission wave, is reflected by a target, arrives as a reflected wave, and is received by the reception antenna 21 at the reception signal fr (t). ).

このとき、図3の上段に示すように、受信信号fr(t)は、自車両MCと物標との距離に応じて、送信信号fs(t)に対して時間差Tだけ遅延している。この時間差Tと、自車両MCおよび物標の相対速度に基づくドップラー効果とにより、ビート信号は、周波数が上昇する「UP区間」の周波数fupと、周波数が下降する「DN区間」の周波数fdnとが繰り返される信号として得られる(図3の中段参照)。   At this time, as shown in the upper part of FIG. 3, the reception signal fr (t) is delayed from the transmission signal fs (t) by a time difference T according to the distance between the host vehicle MC and the target. Due to the time difference T and the Doppler effect based on the relative speeds of the host vehicle MC and the target, the beat signal has a frequency fup in an “UP section” where the frequency increases and a frequency fdn in a “DN section” where the frequency decreases. Are obtained as a repeated signal (see the middle part of FIG. 3).

図3の下段には、かかるビート信号を周波数解析処理においてFFT処理した結果を、「UP区間」側および「DN区間」側のそれぞれについて模式的に示している。   The lower part of FIG. 3 schematically shows the result of the FFT processing of the beat signal in the frequency analysis processing for each of the “UP section” side and the “DN section” side.

図3の下段に示すように、FFT処理後には、「UP区間」側および「DN区間」側のそれぞれの周波数領域における波形が得られる。ピーク抽出処理では、かかる波形においてピークとなるピーク周波数を抽出する。   As shown in the lower part of FIG. 3, after the FFT processing, waveforms in the respective frequency regions on the “UP section” side and the “DN section” side are obtained. In the peak extraction processing, a peak frequency that becomes a peak in such a waveform is extracted.

たとえば、図3の下段に示した例の場合、ピーク抽出閾値が用いられ、「UP区間」側においては、ピークPu1〜Pu3がそれぞれピークとして判定され、ピーク周波数fu1〜fu3がそれぞれ抽出される。   For example, in the case of the example shown in the lower part of FIG. 3, a peak extraction threshold is used. On the “UP section” side, peaks Pu1 to Pu3 are determined as peaks, respectively, and peak frequencies fu1 to fu3 are extracted.

また、「DN区間」側においては、同じくピーク抽出閾値により、ピークPd1〜Pd3がそれぞれピークとして判定され、ピーク周波数fd1〜fd3がそれぞれ抽出される。   On the “DN section” side, the peaks Pd1 to Pd3 are determined as peaks by the same peak extraction threshold, and the peak frequencies fd1 to fd3 are respectively extracted.

ここで、ピーク抽出処理により抽出した各ピーク周波数の周波数成分には、複数の物標からの反射波が混成している場合がある。そこで、瞬時データ生成処理では、各ピーク周波数のそれぞれについて方位演算する角度推定処理を行い、ピーク周波数ごとに対応する物標の存在を解析する。   Here, the frequency component of each peak frequency extracted by the peak extraction processing may include a mixture of reflected waves from a plurality of targets. Therefore, in the instantaneous data generation process, an angle estimation process for calculating the azimuth for each of the peak frequencies is performed, and the presence of a target corresponding to each peak frequency is analyzed.

なお、瞬時データ生成処理における方位演算は、たとえばESPRIT(Estimation of Signal Parameters via Rotational Invariance Techniques)などの公知の到来方向推定手法を用いて行うことができる。   The azimuth calculation in the instantaneous data generation processing can be performed using a known arrival direction estimation technique such as ESPRIT (Estimation of Signal Parameters via Rotational Invariance Techniques).

図4Aは、瞬時データ生成処理の方位演算結果を模式的に示すものである。瞬時データ生成処理では、かかる方位演算結果の各ピークPu1〜Pu3から、これらピークPu1〜Pu3にそれぞれ対応する各物標(各反射点)の推定角度を算出する。また、各ピークPu1〜Pu3の大きさがパワー値となる。瞬時データ生成処理では、図4Bに示すように、かかる角度推定処理を「UP区間」側および「DN区間」側のそれぞれについて行う。   FIG. 4A schematically shows the azimuth calculation result of the instantaneous data generation process. In the instantaneous data generation processing, from the peaks Pu1 to Pu3 of the azimuth calculation result, the estimated angles of the targets (each reflection point) respectively corresponding to these peaks Pu1 to Pu3 are calculated. The magnitude of each of the peaks Pu1 to Pu3 is a power value. In the instantaneous data generation process, as shown in FIG. 4B, the angle estimation process is performed on each of the “UP section” side and the “DN section” side.

そして、瞬時データ生成処理では、方位演算結果において、推定角度およびパワー値の近い各ピークを組み合わせるペアリング処理を行う。また、その組み合わせ結果から、瞬時データ生成処理では、各ピークの組み合わせに対応する各物標(各反射点)の距離および自車両MCへの向きの相対速度を算出する。   Then, in the instantaneous data generation process, a pairing process of combining peaks having close estimated angles and power values in the direction calculation result is performed. Further, from the combination result, in the instantaneous data generation processing, the distance of each target (each reflection point) corresponding to each combination of peaks and the relative speed of the direction to the own vehicle MC are calculated.

距離は、「距離∝(fup+fdn)」の関係に基づいて算出することができる。相対速度は、「速度∝(fup−fdn)」の関係に基づいて算出することができる。その結果、図4Cに示すように、自車両MCに対する、各反射点RPの推定角度、距離および相対速度の瞬時データ100を示すペアリング処理結果が得られる。   The distance can be calculated based on the relationship of “distance∝ (fup + fdn)”. The relative speed can be calculated based on the relationship of “speed∝ (fup−fdn)”. As a result, as shown in FIG. 4C, a pairing processing result indicating the instantaneous data 100 of the estimated angle, the distance, and the relative speed of each reflection point RP with respect to the host vehicle MC is obtained.

図2に戻って、フィルタ処理部32bについて説明する。図2に示すように、フィルタ処理部32bは、予測部321bと、割り当て部322bと、推定部323bと、重み付け部324bと、リサンプリング部325bと、物標データ生成部326bとを備える。   Returning to FIG. 2, the filter processing unit 32b will be described. As illustrated in FIG. 2, the filter processing unit 32b includes a prediction unit 321b, an assignment unit 322b, an estimation unit 323b, a weighting unit 324b, a resampling unit 325b, and a target data generation unit 326b.

フィルタ処理部32bは、生成部32aによって生成された瞬時データ100に対して所定数の粒子データを割り当てるパーティクルフィルタを施すことによって、瞬時データ100に対応する物標データを生成する。   The filter processing unit 32b generates target data corresponding to the instantaneous data 100 by applying a particle filter that allocates a predetermined number of particle data to the instantaneous data 100 generated by the generation unit 32a.

予測部321bは、パーティクルフィルタにおけるサンプル点(粒子データ)の予測処理を行う。具体的には、予測部321bは、最新の周期を時間tとし、時間tにおける粒子データの分布状態Xとした場合、前回の周期の時間t−1の分布状態Xt−1に基づく確率密度関数に基づいてN個の粒子データを配置(サンプリング)する。つまり、予測部321bは、予測処理において、時間t−1の粒子データから時間tにおいて瞬時データ100が現れそうな領域に粒子データを分布させる。 The prediction unit 321b performs a process of predicting sample points (particle data) in the particle filter. Probability Specifically, the prediction unit 321b includes the latest period time t, when the distribution X t of the particle data at time t, based on the distribution state X t-1 time t-1 of the previous cycle N pieces of particle data are arranged (sampled) based on the density function. That is, in the prediction process, the prediction unit 321b distributes the particle data from the particle data at the time t−1 to an area where the instantaneous data 100 is likely to appear at the time t.

また、予測部321bは、前回の物標データに基づいて今回の物標データに対応する予測データを生成する。具体的には、予測部321bは、前回の物標データの移動向きおよびかかる移動向きへの相対速度に基づいて予測データを生成する。   Further, the prediction unit 321b generates prediction data corresponding to the current target data based on the previous target data. Specifically, the prediction unit 321b generates prediction data based on the previous moving direction of the target data and the relative speed in the moving direction.

また、予測部321bは、新規の物標に対応する瞬時データ100については、前回の周期の粒子データが存在しないため、所定の分布状態の粒子データを分布させる。   Further, the prediction unit 321b distributes the particle data in a predetermined distribution state for the instantaneous data 100 corresponding to the new target because there is no particle data in the previous cycle.

なお、予測部321bは、移動向きへの相対速度に基づいて予測データを生成する場合に限定されるものではなく、例えば、移動向きへの対地速度に基づいて予測データを生成してもよい。ここで、対地速度を用いた予測処理について、図5を用いて説明する。   Note that the prediction unit 321b is not limited to the case where the prediction data is generated based on the relative speed in the moving direction. For example, the prediction unit 321b may generate the prediction data based on the ground speed in the moving direction. Here, the prediction process using the ground speed will be described with reference to FIG.

図5は、対地速度に基づく予測処理を示す図である。なお、対地速度の算出方法については後述する。図5では、時刻tにおける物標検出装置1(t)が右前方へ旋回して時刻t+1における物標検出装置1(t+1)へ移動する場合において、時刻tの物標データ50(t)から時刻t+1の予測データ60(t+1)を予測する予測処理について説明する。つまり、予測処理とは、物標データ自体の予測処理と、自車両MCの移動にともなう座標系の変換を行う処理である。以下では、物標データ自体の予測処理を先に行い、次いで座標系変換を行う手順を例に挙げて説明する。なお、座標系変換を先に行い、次いで予測処理を行っても問題はない。   FIG. 5 is a diagram illustrating a prediction process based on a ground speed. The method of calculating the ground speed will be described later. In FIG. 5, when the target detection device 1 (t) at time t turns rightward and moves to the target detection device 1 (t + 1) at time t + 1, the target data 50 (t) at time t A prediction process for predicting the prediction data 60 (t + 1) at the time t + 1 will be described. In other words, the prediction process is a process of predicting the target data itself and a process of performing a coordinate system conversion accompanying the movement of the vehicle MC. In the following, a description will be given by taking as an example a procedure of first performing prediction processing of the target data itself and then performing coordinate system conversion. Note that there is no problem even if the coordinate system conversion is performed first, and then the prediction processing is performed.

図5に示すように、予測部321bは、まず、物標データ50(t)を移動ベクトルV(対地速度)の分だけ移動させて座標系変換前の一時データ50(t+1)を生成する。つまり、前回測定時である時刻tの自車両MCの位置、方向に基づく座標系で表した場合の物標の移動先を一時データ50(t+1)として予測する。   As shown in FIG. 5, the prediction unit 321b first moves the target data 50 (t) by the movement vector V (ground speed) to generate temporary data 50 (t + 1) before the coordinate system conversion. In other words, the destination of the target in the coordinate system based on the position and direction of the vehicle MC at the time t at the previous measurement is predicted as the temporary data 50 (t + 1).

時刻tの自車両MCの状態に基づいて、座標系変換後の原点と基準方向となる時刻t+1における物標検出装置1(t+1)の位置と方向とを予測する。より具体的には、予測部321bは、自車両MCの走行速度および旋回半径に基づいて時刻t+1における物標検出装置1(t+1)の位置を予測する。さらに具体的には、予測部321bは、時刻tから時刻t+1へのXY平面における位置変化(ΔxおよびΔy)と、方向の変化、つまり回転角度θとを予測する。つまり、この時点において、自車両MCの並進移動量と、回転移動量とが求まる。   Based on the state of the host vehicle MC at time t, the position and direction of the target detection device 1 (t + 1) at time t + 1 which is the reference direction after the coordinate system conversion is predicted. More specifically, the prediction unit 321b predicts the position of the target detection device 1 (t + 1) at time t + 1 based on the traveling speed and the turning radius of the vehicle MC. More specifically, the prediction unit 321b predicts a position change (Δx and Δy) on the XY plane from time t to time t + 1 and a change in direction, that is, a rotation angle θ. That is, at this time, the translational movement amount and the rotational movement amount of the vehicle MC are obtained.

そして、予測部321bは、かかる位置変化(ΔxおよびΔy)および回転角度θに基づいて予測データ60(t+1)の予測位置を決定する。つまり、一時データ50(t+1)に対して、上記の並進移動量と回転移動量に基づいた座標系変換を行う。具体的には、予測部321bは、まず、一時データ50(t+1)に対して並進移動変換を行う。具体的には上記した位置変化(ΔyおよびΔx)を基に移動させる。次に、回転移動変換を行う。具体的には、並進移動変換後の一時データ50(t+1)に対して回転角度θだけ回転させる。並進移動変換を先に行うと、回転移動変換における回転中心、つまり原点と、実際の回転中心、つまり物標検出装置1(t+1)の位置とが一致して都合がよい。以上の操作はつまり、予測部321bは、時刻tから時刻t+1までの物標検出装置1の変化量に基づき、時刻t+1における物標検出装置1の位置ならびに方向を原点および基準方向した場合の予測データ60(t+1)の位置を生成する。これにより、物標検出装置1(t),1(t+1)をXY平面の原点に揃えた場合における物標データ50(t)から予測データ60(t+1)への位置変化を予測できる。なお、上記の例では並進移動変換と回転移動変換とを別々に適用したが、この限りではない。例えば並進移動変換や回転移動変換は、その操作を複合したものを一つの線形変換行列として記述することもできる。また回転移動量は回転角度θを変数として記載したが、この限りでもない。方向ベクトルや複素数、四元数などを用いて記述する手法をとってもなんら問題はない。また、座標系変換手法は種々あり、要求に応じて種々選択することができる。   Then, the prediction unit 321b determines the predicted position of the prediction data 60 (t + 1) based on the position change (Δx and Δy) and the rotation angle θ. That is, the coordinate system conversion is performed on the temporary data 50 (t + 1) based on the translation amount and the rotation amount. Specifically, the prediction unit 321b first performs translational conversion on the temporary data 50 (t + 1). Specifically, the movement is performed based on the position change (Δy and Δx) described above. Next, rotational movement conversion is performed. Specifically, the temporary data 50 (t + 1) after the translation conversion is rotated by the rotation angle θ. If the translational conversion is performed first, the rotation center in the rotational movement conversion, that is, the origin, and the actual rotation center, that is, the position of the target detection device 1 (t + 1) are convenient. That is, the prediction unit 321b performs the prediction when the position and the direction of the target detection device 1 at the time t + 1 are set to the origin and the reference direction based on the amount of change of the target detection device 1 from the time t to the time t + 1. The position of the data 60 (t + 1) is generated. This makes it possible to predict a position change from the target data 50 (t) to the prediction data 60 (t + 1) when the target detection devices 1 (t) and 1 (t + 1) are aligned with the origin on the XY plane. In the above example, the translation conversion and the rotation conversion are separately applied, but this is not a limitation. For example, in the case of the translational transformation and the rotational transformation, a combination of the operations can be described as one linear transformation matrix. The rotation amount is described using the rotation angle θ as a variable, but is not limited thereto. There is no problem in using a method of describing using a direction vector, a complex number, a quaternion, or the like. In addition, there are various coordinate system conversion methods, and various methods can be selected as required.

図2に戻って、割り当て部322bについて説明する。割り当て部322bは、最新の周期における瞬時データ100を、予測部321bの予測結果である最新の粒子データへ割り当てる処理を行う。具体的には、割り当て部322bは、予測部321bによって生成された予測データに対応する所定の割り当て範囲内に存在する瞬時データ100に対して、当該予測データに対応する粒子データを割り当てる。   Returning to FIG. 2, the assignment unit 322b will be described. The allocating unit 322b performs a process of allocating the instantaneous data 100 in the latest cycle to the latest particle data, which is the prediction result of the prediction unit 321b. Specifically, the allocating unit 322b allocates the particle data corresponding to the prediction data to the instantaneous data 100 existing within a predetermined allocation range corresponding to the prediction data generated by the prediction unit 321b.

なお、割り当て部322bは、いずれの物標データの割り当て範囲内にも存在しない瞬時データ100があった場合には、かかる瞬時データ100を新規の物標として扱う。   In addition, when there is the instantaneous data 100 that does not exist in the allocation range of any target data, the allocation unit 322b treats the instantaneous data 100 as a new target.

推定部323bは、生成部32aによって生成された複数の瞬時データ100に基づいて物標の移動ベクトルVを推定する。具体的には、推定部323bは、瞬時データ100における自車両MCへの向きの対地速度に基づき、複数の瞬時データ100それぞれに対応する自車両MCへの向きの対地速度を示す対地ベクトルGVにおける起点を揃え、対地ベクトルGVに対する垂線VLの交点CPから物標の移動ベクトルVを推定する。   The estimation unit 323b estimates the movement vector V of the target based on the plurality of instantaneous data 100 generated by the generation unit 32a. Specifically, the estimation unit 323b calculates the ground vector GV indicating the ground speed in the direction toward the own vehicle MC corresponding to each of the plurality of instantaneous data 100 based on the ground speed in the direction toward the own vehicle MC in the instantaneous data 100. The starting points are aligned, and the movement vector V of the target is estimated from the intersection CP of the perpendicular VL to the ground vector GV.

ここで、図6〜図10を用いて、推定部323bの処理内容について具体的には説明する。図6は、推定部323bの機能ブロック図である。図6に示すように、推定部323bは、瞬時データ選択部323baと、バッファ更新部323bbと、対地速度算出部323bc(算出部の一例)と、バッファリング部323bdと、垂線作成部323beと、交点作成部323bfと、内心作成部323bgと、移動ベクトル作成部323bhとを備える。   Here, the processing content of the estimation unit 323b will be specifically described with reference to FIGS. FIG. 6 is a functional block diagram of the estimation unit 323b. As illustrated in FIG. 6, the estimation unit 323b includes an instantaneous data selection unit 323ba, a buffer update unit 323bb, a ground speed calculation unit 323bc (an example of a calculation unit), a buffering unit 323bd, and a perpendicular line creation unit 323be. An intersection creation unit 323bf, an inner center creation unit 323bg, and a movement vector creation unit 323bh are provided.

瞬時データ選択部323baは、移動ベクトルVの推定処理に用いる瞬時データ100を選択する。具体的には、瞬時データ選択部323baは、瞬時データ100に含まれる自車両MCへの向きの相対速度(または、対地速度)や、自車両MCへの向きに相当する角度および距離に基づいて移動ベクトルVの推定処理では不要な瞬時データ100を除外する処理を行う。   The instantaneous data selection unit 323ba selects the instantaneous data 100 used for the process of estimating the movement vector V. Specifically, the instantaneous data selection unit 323ba determines the relative speed (or ground speed) of the direction to the own vehicle MC included in the instantaneous data 100, and the angle and the distance corresponding to the direction to the own vehicle MC. In the process of estimating the movement vector V, a process of excluding unnecessary instantaneous data 100 is performed.

例えば、瞬時データ選択部323baは、複数の瞬時データ100のうち、相対ベクトルRV(または、対地ベクトルGV)の自車両MCへの向きが類似する瞬時データ100を移動ベクトルVの推定処理から除外する。具体的には、瞬時データ選択部323baは、2つの瞬時データ100の角度差(例えば、略ゼロ)が所定の閾値未満の場合、2つの瞬時データ100のいずれか一方を除外する。   For example, the instantaneous data selection unit 323ba excludes the instantaneous data 100 in which the direction of the relative vector RV (or the ground vector GV) to the own vehicle MC is similar from the plurality of instantaneous data 100 from the estimation process of the movement vector V. . Specifically, the instantaneous data selection unit 323ba excludes one of the two instantaneous data 100 when the angle difference (for example, substantially zero) between the two instantaneous data 100 is less than a predetermined threshold.

かかる場合、瞬時データ選択部323baは、例えば、2つの瞬時データ100のうち、自車両MCまでの距離が近い方の瞬時データ100や、角度のパワー値が大きい方の瞬時データ100を残すようにする。これは、角度が略同じで、相対速度(または、対地速度)が異なる2つの瞬時データ100は、同一の物標からは得られないことに起因している。また、自車両MCへの向きに相当する角度が略同じで、相対速度(または、対地速度)が同じ2つの瞬時データ100は、物標において同じような位置の反射点に由来していると考えられ、移動ベクトルVの推定処理においては、いずれか一方のみを用いれば情報量として足りるため除外する。   In such a case, the instantaneous data selection unit 323ba, for example, leaves the instantaneous data 100 having a shorter distance to the host vehicle MC and the instantaneous data 100 having a larger angle power value among the two instantaneous data 100. I do. This is because two instantaneous data 100 having substantially the same angle but different relative speeds (or ground speeds) cannot be obtained from the same target. Further, it is assumed that two instantaneous data 100 having substantially the same angle corresponding to the direction to the host vehicle MC and the same relative speed (or ground speed) are derived from reflection points at similar positions on the target. It is conceivable that in the process of estimating the movement vector V, if only one of them is used, the amount of information is sufficient, so that it is excluded.

このように、角度が略同じ複数の瞬時データ100のうち、少なくとも1つを残して、残りを除外することで、移動ベクトルVの推定精度を落とすことなく、推定部323bの推定処理における処理量が嵩むことを防止できる。   As described above, by removing at least one of the plurality of instantaneous data 100 having substantially the same angle and excluding the rest, the processing amount in the estimation process of the estimation unit 323b can be performed without lowering the estimation accuracy of the movement vector V. Can be prevented from increasing.

また、瞬時データ選択部323baは、複数の瞬時データ100のうち、自車両MCへの向きの相対速度が類似、かつ、自車両MCまでの距離が類似の瞬時データ100を移動ベクトルVの推定処理から除外する。   Further, the instantaneous data selection unit 323ba estimates the instantaneous data 100 having a similar relative speed in the direction to the own vehicle MC and a similar distance to the own vehicle MC among the plurality of instantaneous data 100, by estimating the movement vector V. Exclude from

具体的には、瞬時データ選択部323baは、2つの瞬時データ100の相対速度差(例えば、略ゼロ)が所定の閾値未満、かつ、距離差(例えば、略ゼロ)が所定の閾値未満の場合、2つの瞬時データ100のいずれか一方を除外する。かかる場合、瞬時データ選択部323baは、角度のパワー値が大きい方の瞬時データ100や、角度が物標データ近い方の瞬時データ100を残すようにする。これは、相対速度が略同じ、かつ、距離が略同じ2つの瞬時データ100は原理上存在せず、異なる物標の瞬時データ100または角度割れである可能性が高いためである。   Specifically, the instantaneous data selecting unit 323ba determines that the relative speed difference (for example, approximately zero) between the two instantaneous data 100 is less than a predetermined threshold value and the distance difference (for example, approximately zero) is less than a predetermined threshold value. One of the two instantaneous data 100 is excluded. In such a case, the instantaneous data selection unit 323ba leaves the instantaneous data 100 having the larger angle power value or the instantaneous data 100 having the angle closer to the target data. This is because the two instantaneous data 100 having substantially the same relative speed and the same distance do not exist in principle, and there is a high possibility that the instantaneous data 100 or the angle crack of a different target is present.

これにより、他の物標の瞬時データ100を用いることによる移動ベクトルVの推定精度の低下を防ぐことができるとともに、推定部323bの推定処理における処理量が嵩むことを防止できる。   Thus, it is possible to prevent a decrease in the estimation accuracy of the movement vector V due to the use of the instantaneous data 100 of another target, and to prevent an increase in the amount of processing in the estimation processing of the estimation unit 323b.

次に、バッファ更新部323bbは、記憶部33に記憶された履歴データ33aの更新処理を行う。バッファ更新部323bbによる更新処理は、例えば、瞬時データ選択部323baによって選択された新たな瞬時データ100を履歴データ33aに加える前に行われるが、測定周期に1回行われればどのタイミングでもいい。処理内容とバッファ内容との関連に応じて適宜実施タイミングは設定するとよい。具体的には、バッファ更新部323bbは、新たな瞬時データ100を履歴データ33aとして記憶する場合に、自車両MCの回転角度の変化量に基づいて履歴データ33aにおける過去の瞬時データ100の自車両MCへの向きに相当する角度を補正する。ここで、図7を用いて、バッファ更新部323bbの更新処理についてより具体的に説明する。   Next, the buffer update unit 323bb updates the history data 33a stored in the storage unit 33. The update process by the buffer update unit 323bb is performed, for example, before adding the new instantaneous data 100 selected by the instantaneous data selection unit 323ba to the history data 33a, but may be performed at any timing as long as it is performed once in the measurement cycle. The execution timing may be appropriately set according to the relationship between the processing content and the buffer content. Specifically, when the new instantaneous data 100 is stored as the history data 33a, the buffer update unit 323bb uses the own vehicle of the past instantaneous data 100 in the history data 33a based on the change amount of the rotation angle of the own vehicle MC. The angle corresponding to the direction to the MC is corrected. Here, the update process of the buffer update unit 323bb will be described more specifically with reference to FIG.

図7は、バッファ更新部323bbの更新処理を示す図である。図7には、更新処理の前後における履歴データ33aを示している。具体的には、履歴データ33aには、「ID」、「相対速度」、「対地速度」、「角度」および「保存経過カウント」といった項目が含まれる。   FIG. 7 is a diagram illustrating an update process of the buffer update unit 323bb. FIG. 7 shows the history data 33a before and after the update processing. Specifically, the history data 33a includes items such as “ID”, “relative speed”, “ground speed”, “angle”, and “storage progress count”.

「ID」は、各瞬時データ100を識別する識別情報である。「相対速度」は、瞬時データ100における自車両MCへの向きの相対速度を示す。「対地速度」は、後述の対地速度算出部323bcによって算出される自車両MCへの向きの対地速度を示す。「角度」は、上記の瞬時データ生成処理の角度推定処理で推定した角度であり、自車両MCへの向きに相当する方向である。具体的には、「角度」は、自車両MCの正面方向をゼロとした場合の左右方向への角度を示す。「保存経過カウント」は、最新周期のスキャンから経過した周期を示す。例えば、保存経過カウントが「1」は、最新周期の1つ前の周期、すなわち時刻t−1で得られた瞬時データ100であることを示す。   “ID” is identification information for identifying each instantaneous data 100. “Relative speed” indicates the relative speed of the direction to the host vehicle MC in the instantaneous data 100. The “ground speed” indicates the ground speed in the direction toward the host vehicle MC calculated by the ground speed calculation unit 323bc described below. “Angle” is the angle estimated in the angle estimation processing of the instantaneous data generation processing, and is a direction corresponding to the direction toward the host vehicle MC. Specifically, the “angle” indicates an angle in the left-right direction when the front direction of the vehicle MC is set to zero. The “storage progress count” indicates a cycle that has elapsed since the scan of the latest cycle. For example, the storage progress count “1” indicates that the current cycle is one cycle before the latest cycle, that is, the instantaneous data 100 obtained at time t−1.

図7に示すように、バッファ更新部323bbは、自車両MCの回転角度の変化量に基づいて履歴データ33aの「角度」および「保存経過カウント」を更新する。自車両MCの回転角度の変化量とは、自車両MCの旋回に伴う水平方向への回転角度である。具体的には、回転角度の変化量(θとする)は、自車両MCの走行速度をVself、自車両MCの旋回半径をcrvR、時間の変化量をΔtとした場合、以下の式で算出される。すなわち、θ=Vself×Δt/crvRにより算出される。また、旋回の方向、つまり右向きか左向きか、つまりθの正負は例えばステアリング角などの情報を用いて決定する。   As illustrated in FIG. 7, the buffer update unit 323bb updates the “angle” and the “save progress count” of the history data 33a based on the amount of change in the rotation angle of the vehicle MC. The amount of change in the rotation angle of the host vehicle MC is a rotation angle in the horizontal direction accompanying the turning of the host vehicle MC. Specifically, the change amount of the rotation angle (referred to as θ) is calculated by the following equation, where Vself is the traveling speed of the own vehicle MC, crvR is the turning radius of the own vehicle MC, and Δt is the change amount of time. Is done. That is, it is calculated by θ = Vself × Δt / crvR. The direction of turning, that is, rightward or leftward, that is, positive or negative of θ is determined using information such as a steering angle.

図7に示す例では、自車両MCの回転角度の変化量(θ)が、−3度(例えば、右前方へ旋回)であったとする。かかる場合、バッファ更新部323bbは、履歴データ33aに含まれるすべての過去の瞬時データ100における「角度」を「−3度」減算する。変化量で言えば「3度」加算される。また、バッファ更新部323bbは、履歴データ33aに含まれるすべての過去の瞬時データ100における「保存経過カウント」を「1」増やす。これら過去の瞬時データ100は、後段の移動ベクトルVの推定処理において、新たな瞬時データ100が所定数未満の場合に用いられる。   In the example illustrated in FIG. 7, it is assumed that the change amount (θ) of the rotation angle of the own vehicle MC is −3 degrees (for example, turning right forward). In such a case, the buffer update unit 323bb subtracts "-3 degrees" from the "angle" in all the past instantaneous data 100 included in the history data 33a. In terms of the amount of change, "3 times" is added. Further, the buffer updating unit 323bb increases the “save progress count” in all the past instantaneous data 100 included in the history data 33a by “1”. These past instantaneous data 100 are used when the number of new instantaneous data 100 is less than a predetermined number in the subsequent process of estimating the movement vector V.

つまり、バッファ更新部323bbは、過去の瞬時データ100に対して、予測部321bで行ったような自車両MCの旋回に伴う座標系変換を適用するものである。つまり、記憶している過去の瞬時データ100も、現在の瞬時データ100と同様に、現在の自車両MCの座標系における方向を有することになる。そのため、後段の移動ベクトルVの推定処理において、過去の瞬時データ100を、最新の瞬時データ100の「角度」と同じ基準に揃えることができる。したがって、移動ベクトルVの推定処理において、過去の瞬時データ100を、現在の瞬時データ100と同様に使用でき、結果、推定精度を向上できる。   That is, the buffer update unit 323bb applies the coordinate system conversion accompanying the turning of the own vehicle MC to the past instantaneous data 100 as performed by the prediction unit 321b. That is, the stored past instantaneous data 100 also has the current direction in the coordinate system of the host vehicle MC, similarly to the present instantaneous data 100. Therefore, in the estimation processing of the movement vector V at the subsequent stage, the past instantaneous data 100 can be aligned with the same reference as the “angle” of the latest instantaneous data 100. Therefore, in the estimation process of the movement vector V, the past instantaneous data 100 can be used in the same manner as the current instantaneous data 100, and as a result, the estimation accuracy can be improved.

なお、図7では、瞬時データ100の「角度」を加算する場合について説明したが、自車両MCの回転角度の変化量が正の値(例えば、左前方へ旋回)であった場合、瞬時データ100の「角度」にかかる値が減算される。なお、上記では履歴データ33aの角度を都度補正する手法を例示したがこの限りではない。例えば各測定タイミングnにおける自車両MCの回転角度の変化量θnのみを記録しておき、履歴データ33aを使用する時にθnを減算しつづけていくという手法もある。この場合、時刻t−kの履歴データ33aに対する補正された角度θはθ=θn+θn−1+……+θn−kとなる。θnは回転角度の変化量ではなく積算値、すなわち各測定タイミングnにおける自車角度そのものでも構わない。その場合はθ=θn−θn−kとなる。また、速度と旋回半径からではなく、電子コンパスなどを用いて直接測定した自車の方向を用いてもよい。また角度として例示したが、方向ベクトルとしてもよいし、複素数や四元数を用いた演算を用いてもよい。回転移動の座標系変換には種々の方法があり、目的に応じて適宜選択できる。   FIG. 7 illustrates the case where the “angle” of the instantaneous data 100 is added. However, when the change amount of the rotation angle of the host vehicle MC is a positive value (for example, turning left forward), the instantaneous data is added. The value for the “angle” of 100 is subtracted. In the above, the method of correcting the angle of the history data 33a each time has been exemplified, but the present invention is not limited thereto. For example, there is a method in which only the change amount θn of the rotation angle of the host vehicle MC at each measurement timing n is recorded, and θn is continuously subtracted when using the history data 33a. In this case, the corrected angle θ for the history data 33a at the time tk is θ = θn + θn−1 +... + Θn−k. θn may be an integrated value, not the change amount of the rotation angle, that is, the own vehicle angle itself at each measurement timing n. In that case, θ = θn−θn−k. Instead of the speed and the turning radius, the direction of the own vehicle directly measured using an electronic compass or the like may be used. Although the angle has been exemplified, the direction vector may be used, or an operation using a complex number or a quaternion may be used. There are various methods for converting the coordinate system of the rotational movement, and the method can be appropriately selected according to the purpose.

また、図7に示すように、バッファ更新部323bbは、「保存経過カウント」が所定数(図7では「3」)以上となった瞬時データ100については、履歴データ33aから削除する。これは、情報として比較的古くなるなどによって、信頼性が低下した瞬時データ100を除外し、後段の移動ベクトルVの推定処理に使用しないための操作であり、これにより、移動ベクトルVの推定精度を向上できる。なお、上記では経過カウントを超過したら削除、としたが方法としてはこの限りではない。逆に寿命として記録し、更新処理ではカウントを減算し、0になったら除外する手法もある。また測定時刻のカウントのみを記録し、履歴データ33aの更新はせず、現在の測定カウントから所定のカウント数以上過去の履歴データ33aを除外するなどの手法もある。種々の方法があり、目的に応じて適宜選択できる。   Also, as shown in FIG. 7, the buffer update unit 323bb deletes the instantaneous data 100 for which the “storage progress count” has become a predetermined number (“3” in FIG. 7) or more from the history data 33a. This is an operation for excluding the instantaneous data 100 whose reliability has been reduced due to the fact that the information becomes relatively old, and not using the instantaneous data 100 for estimation processing of the movement vector V at the subsequent stage. Can be improved. In the above description, deletion is performed when the elapsed count is exceeded, but the method is not limited to this. Conversely, there is a method of recording the life as a lifetime, subtracting the count in the update process, and excluding the count when the count reaches 0. There is also a method of recording only the count of the measurement time, not updating the history data 33a, and excluding the history data 33a past a predetermined count or more from the current measurement count. There are various methods, which can be appropriately selected according to the purpose.

次に、対地速度算出部323bcは、瞬時データ選択部323baによって選択された瞬時データ100の自車両MCへの向きの相対速度から自車両MCの速度成分を除いた対地速度を算出する。つまり、瞬時データ100の相対速度には、自車両MCの移動に伴う見かけ上の速度成分と、他車両の実際の移動に伴う速度成分が含まれており、対地速度算出部323bcでは、後者の実際の他車両の移動に伴う速度成分を抽出する。すなわち、対地速度とは、静止座標系、つまり、地面などを基準とした座標系において、他車両の実際の移動速度を自車両方向に投影した速度成分ともいえる。ここで、図8を用いて、対地速度の算出処理について説明する。   Next, the ground speed calculation unit 323bc calculates the ground speed by removing the speed component of the own vehicle MC from the relative speed of the direction to the own vehicle MC of the instant data 100 selected by the instant data selection unit 323ba. That is, the relative speed of the instantaneous data 100 includes an apparent speed component associated with the movement of the own vehicle MC and a speed component associated with the actual movement of the other vehicle, and the ground speed calculation unit 323bc calculates the latter The speed component accompanying the actual movement of another vehicle is extracted. That is, the ground speed can be said to be a speed component obtained by projecting the actual moving speed of another vehicle in the direction of the own vehicle in a stationary coordinate system, that is, a coordinate system based on the ground or the like. Here, the ground speed calculation process will be described with reference to FIG.

図8は、対地速度算出部323bcによる対地速度の算出処理を示す図である。図8では、自車両MCの実際の移動向きおよびかかる移動向きへの実際の走行速度を示す自車ベクトルMVを示している。図8に示すように、対地速度算出部323bcは、まず、自車ベクトルMVに基づいて瞬時データ100への向きの自車ベクトルMaVを算出する。具体的には、自車ベクトルMaVは、自車ベクトルMVが自車両MCと瞬時データ100とを通る直線に射影されたベクトルである。この自車ベクトルMaVが、相対ベクトルRVにおける自車両MCの速度成分となる。   FIG. 8 is a diagram illustrating a ground speed calculation process performed by the ground speed calculation unit 323bc. FIG. 8 shows an own moving vector MV indicating the actual moving direction of the own vehicle MC and the actual traveling speed in the moving direction. As shown in FIG. 8, the ground speed calculation unit 323bc first calculates the own vehicle vector MaV in the direction to the instantaneous data 100 based on the own vehicle vector MV. Specifically, the own vehicle vector MaV is a vector obtained by projecting the own vehicle vector MV on a straight line passing through the own vehicle MC and the instantaneous data 100. The host vehicle vector MaV is the speed component of the host vehicle MC in the relative vector RV.

つまり、対地速度算出部323bcは、瞬時データ100の相対ベクトルRVと自車ベクトルMaVとを加算することで、自車両MCへの向きの対地速度を示す対地ベクトルGVを算出する。   That is, the ground speed calculation unit 323bc calculates the ground vector GV indicating the ground speed in the direction toward the own vehicle MC by adding the relative vector RV of the instantaneous data 100 and the own vehicle vector MaV.

次に、バッファリング部323bdは、瞬時データ選択部323baによって選択された新たな瞬時データ100を記憶部33の履歴データ33aに記憶する。なお、履歴データ33aにおける新たな瞬時データ100の「保存経過カウント」は、例えば「0」のように適宜初期値に設定される。   Next, the buffering unit 323bd stores the new instantaneous data 100 selected by the instantaneous data selection unit 323ba in the history data 33a of the storage unit 33. The “storage progress count” of the new instantaneous data 100 in the history data 33a is appropriately set to an initial value such as “0”.

次に、垂線作成部323beは、対地速度算出部323bcによって算出された瞬時データ100の自車両MCへの向きの対地速度を示す対地ベクトルGVに対する垂線VLを作成する。   Next, the perpendicular creating unit 323be creates a perpendicular VL to the ground vector GV indicating the ground speed in the direction toward the host vehicle MC of the instantaneous data 100 calculated by the ground speed calculating unit 323bc.

具体的には、垂線作成部323beは、複数の瞬時データ100それぞれに対応する対地ベクトルGVの起点を揃えた場合に、対地ベクトルGVの終点を通る垂線VLを作成する。例えば、垂線作成部323beは、自車両MCの左右方向への対地速度成分(VY軸)と、自車両MCに対して前後方向への対地速度成分(VX軸)とで表される平面の原点に対地ベクトルGVの起点を揃えたとする(図1C参照)。   Specifically, when the starting points of the ground vector GV corresponding to each of the plurality of instantaneous data 100 are aligned, the perpendicular creating unit 323be creates a perpendicular VL passing through the end point of the ground vector GV. For example, the perpendicular creating unit 323be calculates the origin of a plane represented by a ground speed component of the host vehicle MC in the left-right direction (VY axis) and a ground speed component of the host vehicle MC in the front-rear direction (VX axis). It is assumed that the starting points of the ground vector GV are aligned (see FIG. 1C).

かかる場合、垂線VLの直線式は、VY=a×VX+bで表される。具体的には、対地ベクトルGVのVX軸を基準とする角度αと、対地ベクトルGVの大きさを対地速度Taとした場合、係数aは、a=tan(α+π/2)で表わされ、係数bは、b=Ta/sin(α)で表される。なお、α=0または±2πの場合、係数aが発散するため、以降の処理においては、Vx=Vとして例外処理を行う。   In such a case, the linear equation of the perpendicular VL is represented by VY = a × VX + b. Specifically, assuming that the angle α of the ground vector GV with respect to the VX axis and the ground vector GV is the ground speed Ta, the coefficient a is represented by a = tan (α + π / 2), The coefficient b is represented by b = Ta / sin (α). Note that when α = 0 or ± 2π, the coefficient a diverges. In the subsequent processing, exception processing is performed with Vx = V.

なお、垂線作成部323beは、最新の瞬時データ100の数が不足している場合、履歴データ33aの過去の瞬時データ100を用いて垂線作成処理を行う。つまり、垂線作成部323beは、最新の瞬時データ100の数が所定数未満の場合、記憶部33に記憶された履歴データ33a(過去の瞬時データ100の一例)を移動ベクトルVの推定処理に用いる。   When the number of the latest instantaneous data 100 is insufficient, the perpendicular creating unit 323be performs the perpendicular creating process using the past instantaneous data 100 of the history data 33a. That is, when the number of the latest instantaneous data 100 is less than the predetermined number, the perpendicular line creating unit 323be uses the history data 33a (an example of the past instantaneous data 100) stored in the storage unit 33 for estimating the movement vector V. .

これにより、瞬時データ100の数が不足することで移動ベクトルVの推定精度が低下することを抑えることができる。   Accordingly, it is possible to suppress a decrease in the estimation accuracy of the movement vector V due to a shortage of the instantaneous data 100.

次に、交点作成部323bfは、垂線作成部323beによって作成された垂線VLの交点CPを作成する。具体的には、交点作成部323bfは、交点CPを求める2つの垂線VLをn1およびn2とし、n1の係数aをa[n1]、係数bをb[n1]と表した場合、交点CPの座標VXは、VX=(b[n2]−b[n1])/(a[n1]−a[n2])により算出され、座標VYは、VY=a×VX+bに上記のVXを代入することで算出される。   Next, the intersection creating unit 323bf creates an intersection CP of the perpendicular VL created by the perpendicular creating unit 323be. Specifically, the intersection creating unit 323bf sets the two perpendiculars VL for finding the intersection CP to n1 and n2, and when the coefficient a of n1 is represented by a [n1] and the coefficient b is represented by b [n1], The coordinates VX are calculated by VX = (b [n2] −b [n1]) / (a [n1] −a [n2]), and the coordinates VY are obtained by substituting the above VX into VY = a × VX + b. Is calculated.

次に、内心作成部323bgは、交点作成部323bfによって作成された交点CPが3つ以上であった場合、かかる3つの交点CPにより形成される三角形の内心IP(図9A参照)を作成する。なお、内心作成部323bgによる内心IPの作成方法については、図9Aで後述する。   Next, when there are three or more intersection points CP created by the intersection point creation unit 323bf, the inner center creation unit 323bg creates the triangle inner center IP (see FIG. 9A) formed by the three intersection points CP. Note that a method of creating the inner core IP by the inner core creating unit 323bg will be described later with reference to FIG. 9A.

次に、移動ベクトル作成部323bhは、交点作成部323bfによって作成された交点CPに基づいて物標の移動ベクトルVを作成する。例えば、移動ベクトル作成部323bhは、交点作成部323bfによって作成された交点CPが1つであった場合、ベクトルRVの起点(VY軸およびVX軸の原点)を移動ベクトルVの起点とし、交点CPを移動ベクトルVの終点として推定する。   Next, the movement vector creation unit 323bh creates the movement vector V of the target based on the intersection CP created by the intersection creation unit 323bf. For example, when there is one intersection CP created by the intersection creation unit 323bf, the movement vector creation unit 323bh sets the starting point of the vector RV (the origin of the VY axis and the VX axis) as the starting point of the movement vector V, and Is estimated as the end point of the movement vector V.

つまり、推定部323bは、移動ベクトル作成部323bhによって作成された移動ベクトルVに基づいて、対地ベクトルGVにおける起点から交点CPへの向きが物標の実際の移動向きであると推定し、対地ベクトルGVにおける起点および交点CPの距離が物標の移動向きへの対地速度であると推定する。   That is, based on the movement vector V created by the movement vector creation unit 323bh, the estimation unit 323b estimates that the direction from the starting point to the intersection CP in the ground vector GV is the actual movement direction of the target, and It is estimated that the distance between the starting point and the intersection CP in the GV is the ground speed in the moving direction of the target.

また、移動ベクトル作成部323bhは、交点作成部323bfによって作成された交点CPが3つ以上であった場合、内心作成部323bgによって作成された内心IPに基づいて移動ベクトルVを作成する。ここで、図9Aおよび図9Bを用いて、交点CPが3つ以上の場合における移動ベクトルVの作成方法について説明する。   In addition, when the number of intersection points CP created by the intersection point creation unit 323bf is three or more, the movement vector creation unit 323bh creates the movement vector V based on the inner center IP created by the inner center creation unit 323bg. Here, a method of creating the movement vector V when there are three or more intersection points CP will be described with reference to FIGS. 9A and 9B.

図9Aは、交点CPが3つの場合における移動ベクトルVの作成方法を示す図である。図9Bは、交点CPが6つの場合における移動ベクトルVの作成方法を示す図である。まず、図9Aを用いて、交点CPが3つの場合における移動ベクトルVの作成方法について説明する。   FIG. 9A is a diagram illustrating a method of creating a movement vector V when there are three intersection points CP. FIG. 9B is a diagram illustrating a method of creating the movement vector V when there are six intersection points CP. First, with reference to FIG. 9A, a method of creating a movement vector V when there are three intersection points CP will be described.

まず、内心作成部323bgは、3つの交点CP1,CP2,CP3の内心IPを作成する。具体的には、内心IPは、3つの交点CP1,CP2,CP3に形成される三角形の内接円の中心である。   First, the inner center creating unit 323bg creates inner centers IP of the three intersections CP1, CP2, and CP3. Specifically, the inner center IP is the center of a triangular inscribed circle formed at the three intersections CP1, CP2, and CP3.

より具体的には、内心作成部323bgは、上記したVY軸およびVX軸の平面における原点を起点として、各交点CP1,CP2,CP3を終点とする算出用ベクトルCV1,CV2,CV3を作成する。   More specifically, the inner center creating unit 323bg creates calculation vectors CV1, CV2, and CV3 starting from the origins on the planes of the VY axis and the VX axis and ending at the intersections CP1, CP2, and CP3.

そして、内心作成部323bgは、算出用ベクトルCV1,CV2,CV3と、各交点CP1,CP2,CP3の対辺に基づいて内心IPを作成する。具体的には、内心作成部323bgは、交点CP1,CP2,CP3の対辺をCP1a,CP2a,CP3aと表した場合、原点から内心IPへのベクトル(図9Aでは、移動ベクトルV)は、V=(CP1a×CV1+CP2a×CV2+CP3a×CV3)/(CP1a+CP2a+CP3a)で算出される。   Then, the inner center creating unit 323bg creates an inner center IP based on the calculation vectors CV1, CV2, and CV3 and the opposite sides of the intersections CP1, CP2, and CP3. Specifically, when the opposite sides of the intersections CP1, CP2, and CP3 are represented as CP1a, CP2a, and CP3a, the inner-center creating unit 323bg obtains a vector (movement vector V in FIG. 9A) from the origin to the inner center IP by V = It is calculated by (CP1a × CV1 + CP2a × CV2 + CP3a × CV3) / (CP1a + CP2a + CP3a).

なお、内心作成部323bgは、三角形を形成する3つの交点CPのバラつきが所定値以上の場合には、内心IPの算出処理から除外する。具体的には、内心作成部323bgは、3つの交点CPの不偏分散が所定値以上の場合には、内心IPの算出処理から除外する。   Note that, when the variation of the three intersection points CP forming the triangle is equal to or more than a predetermined value, the inner center creating unit 323bg excludes the variation from the inner center IP calculation process. Specifically, when the unbiased variance of the three intersection points CP is equal to or more than a predetermined value, the inner center creating unit 323bg excludes the unbiased variance from the inner center IP calculation process.

そして、移動ベクトル作成部323bhは、内心IPが1つの場合、原点から内心IPへのベクトルを移動ベクトルVとして作成する。このように、内心IPを算出して移動ベクトルVを作成することで、交点CPが複数の場合であっても、移動ベクトルVの推定精度が低下することを防止できる。   Then, when there is one inner core IP, the movement vector generator 323bh creates a vector from the origin to the inner core IP as the movement vector V. As described above, by calculating the inner center IP and creating the movement vector V, even if there are a plurality of intersections CP, it is possible to prevent the estimation accuracy of the movement vector V from being lowered.

なお、移動ベクトル作成部323bhは、算出用ベクトルCV1,CV2,CV3のうち、いずれか2つの算出用ベクトルCV1,CV2,CV3の誤差が所定値以上であることが明らかである場合等には、残りの1つの算出用ベクトルCV1,CV2,CV3を移動ベクトルVとしてもよい。   In addition, when it is clear that the error of any two calculation vectors CV1, CV2, and CV3 among the calculation vectors CV1, CV2, and CV3 is equal to or more than a predetermined value, the movement vector creation unit 323bh performs, for example, The remaining one calculation vector CV1, CV2, CV3 may be used as the movement vector V.

次に、図9Bを用いて、交点CPが6つの場合における移動ベクトルVの作成方法について説明する。図9Bに示すように、交点CP3〜CP8が6つの場合、かかる交点CP3〜CP8によって形成される三角形が複数となり、従って内心IP1〜IP4も複数となる。なお、図9Bに示す例では、3つの交点CP(例えば、交点CP3,CP5,CP8)により三角形が形成されない組み合わせもあり、かかる交点CPの組み合わせについては内心IPの作成から除外する。   Next, with reference to FIG. 9B, a method of creating the movement vector V in the case where the number of intersection points CP is six will be described. As shown in FIG. 9B, when there are six intersections CP3 to CP8, there are a plurality of triangles formed by the intersections CP3 to CP8, and thus there are also a plurality of inner centers IP1 to IP4. In the example shown in FIG. 9B, there is a combination in which a triangle is not formed by three intersections CP (for example, intersections CP3, CP5, and CP8), and such a combination of intersections CP is excluded from the creation of the inner center IP.

図9Bに示すように、移動ベクトル作成部323bhは、例えば、内心IP1〜IP4の平均値となる平均点IPAを算出し、原点から平均点IPAへのベクトルを移動ベクトルVとして算出する。   As illustrated in FIG. 9B, the movement vector creation unit 323bh calculates, for example, an average point IPA that is an average value of the inner centers IP1 to IP4, and calculates a vector from the origin to the average point IPA as the movement vector V.

このように、内心IPが複数作成された場合に、平均点IPAによる移動ベクトルVを作成することにより、内心IPのバラつきによる移動ベクトルVの推定誤差を最小限に抑えることができる。   As described above, when a plurality of inner core IPs are created, by creating the movement vector V based on the average point IPA, it is possible to minimize the estimation error of the movement vector V due to the variation of the inner core IP.

なお、移動ベクトル作成部323bhは、複数の内心IPのバラつきが所定値以上の場合には、移動ベクトルVの作成処理から除外する。具体的には、移動ベクトル作成部323bhは、複数の内心IPのうち、不偏分散が所定値以上の内心IPについては、移動ベクトルVの作成処理から除外する。   Note that, when the variation of the plurality of inner cores IP is equal to or more than a predetermined value, the movement vector creation unit 323bh excludes the variation from the creation processing of the movement vector V. Specifically, the movement vector creation unit 323bh excludes, from the plurality of inside IPs, the inside IP whose unbiased variance is equal to or more than a predetermined value from the creation processing of the movement vector V.

図2に戻って重み付け部324bについて説明する。重み付け部324bは、割り当て関係にある今回の瞬時データ100および今回の粒子データについて今回の粒子データそれぞれに対して重みを付ける。   Returning to FIG. 2, the weighting unit 324b will be described. The weighting unit 324b weights each of the current instantaneous data 100 and the current particle data in the assignment relationship with respect to the current particle data.

例えば、重み付け部324bは、今回の粒子データのうち、今回の瞬時データ100に近い粒子データの重みを大きくし、今回の瞬時データ100から遠い粒子データの重みを小さくする。なお、ここでいう「近い」および「遠い」は、マハラノビス距離が「近い」および「遠い」ことを指す。   For example, the weighting unit 324b increases the weight of the particle data close to the current instant data 100 among the current particle data, and decreases the weight of the particle data far from the current instant data 100. Here, “close” and “far” mean that the Mahalanobis distance is “close” and “far”.

また、重み付け部324bは、推定部323bによって推定された移動ベクトルVに基づいて重み付けを行う。ここで、かかる点について、図10を用いて説明する。図10は、移動ベクトルVに基づく重み付け処理を示す図である。   The weighting unit 324b performs weighting based on the movement vector V estimated by the estimation unit 323b. Here, such a point will be described with reference to FIG. FIG. 10 is a diagram illustrating a weighting process based on the movement vector V.

図10では、3つの粒子データに基づくベクトルPV1〜PV3の起点と、移動ベクトルVの起点とをVY軸およびVX軸平面の原点に揃えた場合を示している。なお、粒子データに基づくベクトルPV1〜PV3は、前回の粒子データから今回の粒子データ(予測部321bにより予測された粒子データ)へのベクトルである。   FIG. 10 shows a case where the starting points of the vectors PV1 to PV3 based on the three particle data and the starting point of the movement vector V are aligned with the origins of the VY-axis and VX-axis planes. The vectors PV1 to PV3 based on the particle data are vectors from the previous particle data to the current particle data (particle data predicted by the prediction unit 321b).

例えば、重み付け部324bは、ベクトルPV1〜PV3と移動ベクトルVとの類似度が高いほど重みを大きくする。例えば、重み付け部324bは、ベクトルPV1〜PV3および移動ベクトルVのVY軸およびVX軸の対地速度成分が近いほど、類似度を高くする。   For example, the weighting unit 324b increases the weight as the similarity between the vectors PV1 to PV3 and the movement vector V increases. For example, the weighting unit 324b increases the similarity as the ground speed components on the VY axis and the VX axis of the vectors PV1 to PV3 and the movement vector V are closer.

つまり、図10に示す例では、重み付け部324bは、ベクトルPV1,PV2に由来する粒子データの重みを大きくし、ベクトルPV3に由来する粒子データの重みを小さくする。これにより、後段のリサンプリング部325bにおけるリサンプリングの精度を向上させることができ、算出される物標データの精度を向上させることができる。   That is, in the example shown in FIG. 10, the weighting unit 324b increases the weight of the particle data derived from the vectors PV1 and PV2 and decreases the weight of the particle data derived from the vector PV3. Thereby, the accuracy of resampling in the resampling unit 325b at the subsequent stage can be improved, and the accuracy of the calculated target data can be improved.

次に、リサンプリング部325bは、今回の粒子データそれぞれの重みに基づいて粒子データを再配置(リサンプリング)する。具体的には、リサンプリング部325bは、重みが小さい粒子データを瞬時データ100の近くへ移動させる。   Next, the resampling unit 325b rearranges (resamples) the particle data based on the weights of the current particle data. Specifically, the resampling unit 325b moves the particle data having a small weight closer to the instantaneous data 100.

物標データ生成部326bは、リサンプリング部325bによって再配置された今回の粒子データに基づいて確率密度関数を再計算し、再計算された確率密度関数の重心に基づいて物標データを生成する。なお、物標データ生成部326bは、確率密度関数の重心に基づいて物標データを生成したが、例えば、確率密度関数の平均に基づいて物標データを生成してもよい。   The target data generation unit 326b recalculates the probability density function based on the current particle data rearranged by the resampling unit 325b, and generates target data based on the center of gravity of the recalculated probability density function. . Although the target data generation unit 326b generates the target data based on the barycenter of the probability density function, the target data generation unit 326b may generate the target data based on, for example, an average of the probability density function.

また、物標データ生成部326bは、粒子データが割り当てられなかった瞬時データ100を新規の物標として扱い、そのまま物標データとして出力する。すなわち、物標データ生成部326bは、新規の物標の場合、瞬時データ100=物標データとして出力する。   Further, the target data generation unit 326b treats the instantaneous data 100 to which the particle data has not been assigned as a new target, and outputs the data as it is as target data. That is, in the case of a new target, the target data generation unit 326b outputs the instantaneous data 100 = target data.

次に、図11を用いて、実施形態に係る物標検出装置1が実行する処理の処理手順について説明する。図11は、実施形態に係る物標検出装置1が実行する処理の処理手順を示すフローチャートである。   Next, a processing procedure of processing executed by the target detection device 1 according to the embodiment will be described with reference to FIG. FIG. 11 is a flowchart illustrating a processing procedure of processing performed by the target detection device 1 according to the embodiment.

図11に示すように、まず、生成部32aは、送信した電波が物標で反射した複数の反射点それぞれについて、物標の自車両MCへの向きの相対速度を含む瞬時データ100を生成する(ステップS101)。   As illustrated in FIG. 11, first, the generation unit 32a generates instantaneous data 100 including a relative speed of a direction of the target toward the own vehicle MC at each of a plurality of reflection points where the transmitted radio wave is reflected by the target. (Step S101).

つづいて、フィルタ処理部32bの予測部321bは、前回の粒子データに基づいて今回の粒子データを予測する予測処理を行う(ステップS102)。なお、フィルタ処理部32bは、前回の粒子データが存在しない場合、前回の瞬時データ100の相対速度に基づいて今回の粒子データを予測する。   Subsequently, the prediction unit 321b of the filter processing unit 32b performs a prediction process of predicting the current particle data based on the previous particle data (Step S102). When there is no previous particle data, the filter processing unit 32b predicts the current particle data based on the relative speed of the previous instantaneous data 100.

つづいて、割り当て部322bは、今回の粒子データに今回の瞬時データ100を割り当てる(ステップS103)。つづいて、推定部323bは、複数の今回の瞬時データ100に基づいて物標の移動ベクトルVを推定する(ステップS104)。なお、推定部323bによる推定処理の処理手順については図12で後述する。   Subsequently, the allocating unit 322b allocates the current instantaneous data 100 to the current particle data (Step S103). Subsequently, the estimating unit 323b estimates the movement vector V of the target based on the plurality of current instant data 100 (step S104). The processing procedure of the estimation processing by the estimation unit 323b will be described later with reference to FIG.

つづいて、割り当て部322bは、今回の粒子データが割り当てられなかった瞬時データ100の有無により新規の物標の有無を判定する(ステップS105)。割り当て部322bは、瞬時データ100が新規の物標であった場合(ステップS105,Yes)、新規の物標に対応する瞬時データ100に対して所定の粒子データ(例えば、初期状態の粒子データ)を設定する(ステップS106)。   Subsequently, the assigning unit 322b determines the presence or absence of a new target based on the presence or absence of the instantaneous data 100 to which the current particle data has not been assigned (step S105). When the instantaneous data 100 is a new target (Step S105, Yes), the allocating unit 322b assigns predetermined particle data (for example, initial state particle data) to the instantaneous data 100 corresponding to the new target. Is set (step S106).

つづいて、重み付け部324bは、瞬時データ100が新規の物標でなかった場合(ステップS105,No)、推定部323bによって推定された移動ベクトルVに基づいて今回の粒子データそれぞれに重み付けを行う(ステップS107)。   Subsequently, when the instantaneous data 100 is not a new target (No at Step S105), the weighting unit 324b weights each of the current particle data based on the movement vector V estimated by the estimation unit 323b ( Step S107).

つづいて、リサンプリング部325bは、重み付け部324bによる重み付けに基づいて今回の粒子データのリサンプリングを行う(ステップS108)。つづいて、物標データ生成部326bは、リサンプリングされた今回の粒子データの確率密度関数を更新し、かかる確率密度関数に基づいて物標データを生成し(ステップS109)、処理を終了する。   Subsequently, the resampling unit 325b performs resampling of the current particle data based on the weight by the weighting unit 324b (Step S108). Subsequently, the target data generation unit 326b updates the probability density function of the resampled particle data of this time, generates target data based on the probability density function (step S109), and ends the process.

次に、図12を用いて、実施形態に係る推定部323bが実行する推定処理の処理手順について説明する。図12は、実施形態に係る推定部323bが実行する推定処理の処理手順を示すフローチャートである。   Next, a processing procedure of an estimation process executed by the estimation unit 323b according to the embodiment will be described with reference to FIG. FIG. 12 is a flowchart illustrating a processing procedure of an estimation process performed by the estimation unit 323b according to the embodiment.

図12に示すように、まず、瞬時データ選択部323baは、複数の瞬時データ100のうち、移動ベクトルVの推定処理に用いる瞬時データ100を選択する選択処理を行う(ステップS201)。   As shown in FIG. 12, first, the instantaneous data selection unit 323ba performs a selection process of selecting the instantaneous data 100 used for the estimation process of the movement vector V from the plurality of instantaneous data 100 (step S201).

つづいて、バッファ更新部323bbは、自車両MCの回転角度の変化量に基づいて記憶部33の履歴データ33aにおける過去の瞬時データ100の自車両MCへの向きに相当する角度を補正することで、履歴データ33aの更新処理を行う(ステップS202)。   Subsequently, the buffer update unit 323bb corrects the angle corresponding to the direction of the past instantaneous data 100 to the own vehicle MC in the history data 33a of the storage unit 33 based on the change amount of the rotation angle of the own vehicle MC. The update process of the history data 33a is performed (step S202).

つづいて、対地速度算出部323bcは、瞬時データ100の自車両MCへの向きの相対速度から自車両MCの速度成分を除いた対地速度を算出する(ステップS203)。つづいて、バッファリング部323bdは、瞬時データ選択部323baによって選択された新たな瞬時データ100を履歴データ33aとして記憶するバッファリング処理を行う(ステップS204)。   Subsequently, the ground speed calculation unit 323bc calculates the ground speed by removing the speed component of the own vehicle MC from the relative speed of the direction to the own vehicle MC in the instantaneous data 100 (step S203). Subsequently, the buffering unit 323bd performs a buffering process of storing the new instantaneous data 100 selected by the instantaneous data selection unit 323ba as the history data 33a (Step S204).

つづいて、垂線作成部323beは、対地速度算出部323bcによって算出された瞬時データ100における自車両MCへの向きの対地速度に基づき、複数の瞬時データ100それぞれに対応する対地ベクトルGVにおける起点を揃え、対地ベクトルGVに対する垂線VLを作成する垂線作成処理を行う(ステップS205)。   Subsequently, the perpendicular creating unit 323be aligns the starting points in the ground vector GV corresponding to each of the plurality of instantaneous data 100 based on the ground speed in the direction to the host vehicle MC in the instantaneous data 100 calculated by the ground speed calculating unit 323bc. Then, a perpendicular line creating process for creating a perpendicular line VL with respect to the ground vector GV is performed (step S205).

つづいて、交点作成部323bfは、垂線作成部323beによって作成された垂線VLの交点CPを作成する交点作成処理を行う(ステップS206)。移動ベクトル作成部323bhは、交点作成部323bfによって作成された交点CPが1つであるか否かを判定する(ステップS207)。   Subsequently, the intersection creating unit 323bf performs an intersection creating process of creating an intersection CP of the perpendicular VL created by the perpendicular creating unit 323be (Step S206). The movement vector creation unit 323bh determines whether or not there is one intersection CP created by the intersection creation unit 323bf (Step S207).

移動ベクトル作成部323bhは、交点作成部323bfによって作成された交点CPが1つであった場合(ステップS207,Yes)、移動ベクトル作成処理を行い(ステップS208)、処理を終了する。なお、移動ベクトル作成部323bhは、かかる場合、対地ベクトルGVの起点から交点CPへのベクトルを移動ベクトルVとして作成する。   When there is one intersection CP created by the intersection creation unit 323bf (step S207, Yes), the movement vector creation unit 323bh performs a movement vector creation process (step S208), and ends the process. In this case, the movement vector creation unit 323bh creates a vector from the starting point of the ground vector GV to the intersection CP as the movement vector V.

一方、ステップS207において、内心作成部323bgは、交点CPが1つでない場合(ステップS207,No)、すなわち、交点CPが3つ以上の場合、3つの交点CPにより形成される三角形の内心IPを作成する内心作成処理を行い(ステップS209)、処理をステップS208へ移行する。なお、かかる場合、移動ベクトル作成部323bhは、ステップS208において、内心作成部323bgによって作成された内心IPに基づいて移動ベクトル作成処理を行う。   On the other hand, in step S207, if the number of intersections CP is not one (step S207, No), that is, if there are three or more intersections CP, the inner center creation unit 323bg determines the inner center IP of the triangle formed by the three intersections CP. An inner core creation process to be created is performed (step S209), and the process proceeds to step S208. In such a case, in step S208, the movement vector creation unit 323bh performs a movement vector creation process based on the inner core IP created by the inner core creation unit 323bg.

上述してきたように、実施形態に係る物標検出装置1は、生成部32aと、対地速度算出部323bc(算出部の一例)と、推定部323bとを備える。生成部32aは、送信した電波が物標で反射した複数の反射点それぞれについて、物標の自車両MCへの向きの相対速度を含む瞬時データ100を生成する。対地速度算出部323bcは、生成部32aによって生成された瞬時データ100の自車両MCへの向きの相対速度と自車両MCの速度とに基づいて、対地速度を算出する。推定部323bは、対地速度算出部323bcaによって算出された対地速度に基づき、複数の瞬時データ100それぞれに対応する自車両MCへの向きの対地ベクトルGVにおける起点を揃え、対地ベクトルGVに対する垂線VLの交点CPから物標の移動向き(移動ベクトルVの一要素)を推定する。これにより、物標に対する応答性を向上させることができる。   As described above, the target detection device 1 according to the embodiment includes the generation unit 32a, the ground speed calculation unit 323bc (an example of a calculation unit), and the estimation unit 323b. The generation unit 32a generates instantaneous data 100 including the relative speed of the direction of the target toward the own vehicle MC for each of the plurality of reflection points where the transmitted radio wave is reflected on the target. The ground speed calculating unit 323bc calculates the ground speed based on the relative speed of the direction to the own vehicle MC of the instantaneous data 100 generated by the generating unit 32a and the speed of the own vehicle MC. Based on the ground speed calculated by the ground speed calculation unit 323bca, the estimation unit 323b aligns the starting points in the ground vector GV of the direction to the own vehicle MC corresponding to each of the plurality of instantaneous data 100, and sets the vertical line VL to the ground vector GV. The movement direction of the target (one element of the movement vector V) is estimated from the intersection CP. Thereby, the response to the target can be improved.

上述した実施形態では、物標検出装置1は車両に設けられることとしたが、無論、車両以外の移動体、たとえば船舶や航空機などに設けられてもよい。   In the above-described embodiment, the target detection device 1 is provided in the vehicle. However, it is needless to say that the target detection device 1 may be provided in a moving object other than the vehicle, for example, a ship or an aircraft.

また、上述した実施形態では、物標検出装置1の用いる到来方向推定手法の例にESPRITを挙げたが、これに限られるものではない。たとえばDBF(Digital Beam Forming)や、PRISM(Propagator method based on an Improved Spatial-smoothing Matrix)、MUSIC(Multiple Signal Classification)などを用いてもよい。   Further, in the above-described embodiment, ESPRIT has been described as an example of the direction-of-arrival estimation method used by the target detection device 1, but the present invention is not limited to this. For example, DBF (Digital Beam Forming), PRISM (Propagator method based on an Improved Spatial-smoothing Matrix), MUSIC (Multiple Signal Classification), or the like may be used.

また、実施形態に係る物標検出装置1は、移動ベクトル作成部323bhによって作成された移動ベクトルVの確度を粒子データに基づいて検証してもよい。かかる点について、図13を用いて説明する。   The target detection device 1 according to the embodiment may verify the accuracy of the movement vector V created by the movement vector creation unit 323bh based on the particle data. This will be described with reference to FIG.

図13は、移動ベクトルVの検証処理を示す図である。図13では、移動ベクトルVと、粒子データに基づくベクトルPVとを示している。   FIG. 13 is a diagram illustrating the verification processing of the movement vector V. FIG. 13 shows a movement vector V and a vector PV based on the particle data.

図13に示すように、推定部323bは、作成した移動ベクトルVと粒子データのベクトルPVとの類似度により移動ベクトルVの信頼度を示す確度を算出する。具体的には、推定部323bは、移動ベクトルVおよびベクトルPVのVY軸およびVX軸の対地速度成分の類似度を確度として算出する。   As illustrated in FIG. 13, the estimation unit 323b calculates the accuracy indicating the reliability of the movement vector V based on the similarity between the created movement vector V and the particle data vector PV. Specifically, the estimating unit 323b calculates the similarity between the ground speed components on the VY axis and the VX axis of the movement vector V and the vector PV as accuracy.

例えば、推定部323bは、粒子データのベクトルPVと移動ベクトルVとの類似度が所定値以上である場合に、かかる移動ベクトルVを最終的な推定結果の確定値とする。一方、推定部323bは、粒子データのベクトルPVと移動ベクトルVとの類似度が所定値未満である場合に、かかる移動ベクトルVの使用を禁止する。これにより、より推定精度の高い移動ベクトルVを使用することができる。   For example, when the similarity between the particle data vector PV and the movement vector V is equal to or greater than a predetermined value, the estimation unit 323b sets the movement vector V as a final value of the final estimation result. On the other hand, when the similarity between the particle data vector PV and the movement vector V is less than a predetermined value, the estimation unit 323b prohibits the use of the movement vector V. As a result, it is possible to use the movement vector V with higher estimation accuracy.

なお、図13では、移動ベクトルVおよびベクトルPVのVY軸およびVX軸の対地速度成分の類似度を確度として算出したが、例えば、移動ベクトルVおよびベクトルPVの実際の移動向きの類似度を確度として算出してもよい。   In FIG. 13, the similarity between the ground speed components on the VY axis and the VX axis of the movement vector V and the vector PV is calculated as the accuracy. May be calculated as

なお、上述した実施形態では、フィルタ処理部32bは、パーティクルフィルタを用いた場合を示したが、用いる時系列フィルタは、パーティクルフィルタに限定されるものではなく、例えば、カルマンフィルタや、αβフィルタ等の時系列フィルタであってもよい。   In the above-described embodiment, the filter processing unit 32b shows the case where the particle filter is used. However, the time series filter to be used is not limited to the particle filter. For example, a Kalman filter, an αβ filter, or the like is used. It may be a time series filter.

さらなる効果や変形例は、当業者によって容易に導き出すことができる。このため、本発明のより広範な態様は、以上のように表しかつ記述した特定の詳細および代表的な実施形態に限定されるものではない。したがって、添付の特許請求の範囲およびその均等物によって定義される総括的な発明の概念の精神または範囲から逸脱することなく、様々な変更が可能である。   Further effects and modifications can be easily derived by those skilled in the art. Thus, the broader aspects of the present invention are not limited to the specific details and representative embodiments shown and described above. Accordingly, various modifications can be made without departing from the spirit or scope of the general inventive concept as defined by the appended claims and their equivalents.

1 物標検出装置
32 信号処理部
32a 生成部
32b フィルタ処理部
33 記憶部
33a 履歴データ
100 瞬時データ
321b 予測部
322b 割り当て部
323b 推定部
323ba 瞬時データ選択部
323bb バッファ更新部
323bc 対地速度算出部
323bd バッファリング部
323be 垂線作成部
323bf 交点作成部
323bg 内心作成部
323bh 移動ベクトル作成部
324b 重み付け部
325b リサンプリング部
326b 物標データ生成部
LC 他車両
MC 自車両
Reference Signs List 1 target detection device 32 signal processing unit 32a generation unit 32b filter processing unit 33 storage unit 33a history data 100 instantaneous data 321b prediction unit 322b allocation unit 323b estimation unit 323ba instantaneous data selection unit 323bb buffer update unit 323bc buffer for ground speed calculation 323bbd Ring part 323be Perpendicular line creation part 323bf Intersection creation part 323bg Inner center creation part 323bh Movement vector creation part 324b Weighting part 325b Resampling part 326b Target data generation part LC Other vehicle MC Own vehicle

Claims (9)

送信した電波が物標で反射した複数の反射点それぞれについて、前記物標の自車両への向きの相対速度を含む瞬時データを生成する生成部と、
前記生成部によって生成された前記瞬時データの前記相対速度と前記自車両の速度とに基づいて、対地速度を算出する算出部と、
前記算出部によって算出された前記対地速度に基づき、複数の前記瞬時データそれぞれに対応する前記自車両への向きの対地ベクトルにおける起点を揃え、当該対地ベクトルに対する垂線の交点から前記物標の移動向きを推定する推定部と
を備えることを特徴とする物標検出装置。
For each of the plurality of reflection points of the transmitted radio wave reflected by the target, a generation unit that generates instantaneous data including a relative speed of the direction of the target toward the own vehicle,
A calculating unit that calculates a ground speed based on the relative speed of the instantaneous data generated by the generating unit and the speed of the host vehicle;
Based on the ground speed calculated by the calculation unit, the starting points in the ground vector of the direction to the host vehicle corresponding to each of the plurality of instantaneous data are aligned, and the moving direction of the target from the intersection of the perpendicular to the ground vector. A target detecting device, comprising: an estimating unit for estimating.
前記推定部は、
前記交点が1つであった場合、前記対地ベクトルにおける前記起点から当該交点への向きが前記移動向きであると推定すること
を特徴とする請求項1に記載の物標検出装置。
The estimation unit includes
The target detecting device according to claim 1, wherein when the number of the intersections is one, the direction from the starting point to the intersection in the ground vector is estimated to be the moving direction.
前記推定部は、
前記交点が3つ以上であった場合、該交点により形成される三角形の内心に基づいて前記移動向きを推定すること
を特徴とする請求項1または2に記載の物標検出装置。
The estimation unit includes
3. The target detection device according to claim 1, wherein when the number of the intersections is three or more, the movement direction is estimated based on a center of a triangle formed by the intersections. 4.
前記推定部は、
前記複数の瞬時データのうち、前記対地ベクトルの前記自車両への向きが類似する前記瞬時データを前記移動向きの推定処理から除外すること
を特徴とする請求項1〜3のいずれか1つに記載の物標検出装置。
The estimation unit includes
The method according to any one of claims 1 to 3, wherein, among the plurality of instantaneous data, the instantaneous data in which the direction of the ground vector toward the own vehicle is similar is excluded from the movement direction estimation processing. The target detection device described in the above.
前記瞬時データには、
前記自車両から当該瞬時データの位置までの距離が含まれ、
前記推定部は、
前記生成部によって生成された前記複数の瞬時データのうち、前記相対速度が類似、かつ、前記距離が類似の前記瞬時データを前記移動向きの推定処理から除外すること
を特徴とする請求項1〜4のいずれか1つに記載の物標検出装置。
The instantaneous data includes
The distance from the own vehicle to the position of the instantaneous data is included,
The estimation unit includes
The instantaneous data having a similar relative speed and a similar distance among the plurality of instantaneous data generated by the generation unit are excluded from the estimation process of the moving direction. 4. The target detection device according to any one of 4.
前記生成部によって生成された前記瞬時データを履歴データとして記憶する記憶部をさらに備え、
前記推定部は、
前記生成部によって生成された最新の前記瞬時データの数が所定数未満であった場合、前記履歴データとして記憶された過去の前記瞬時データを前記移動向きの推定処理に用いること
を特徴とする請求項1〜5のいずれか1つに記載の物標検出装置。
A storage unit that stores the instantaneous data generated by the generation unit as history data,
The estimation unit includes
When the number of the latest instantaneous data generated by the generation unit is less than a predetermined number, the past instantaneous data stored as the history data is used for the estimating process of the moving direction. Item 7. The target detection device according to any one of Items 1 to 5.
前記推定部は、
前記自車両の方向の変化に基づき、前記履歴データにおける過去の前記瞬時データにおける前記自車両に対する方向を補正すること
を特徴とする請求項6に記載の物標検出装置。
The estimation unit includes
The target detection device according to claim 6, wherein the direction of the own vehicle in the past instantaneous data in the history data is corrected based on a change in the direction of the own vehicle.
前記生成部によって生成された前記瞬時データに対して複数の粒子データを割り当てるパーティクルフィルタを施すことによって当該瞬時データに対応する物標データを生成するフィルタ処理部をさらに備え、
前記推定部は、
前記複数の粒子データに基づくベクトルの向きと、前記移動向きとの類似度が所定値以上である場合に、当該移動向きを推定結果の確定値とすること
を特徴とする請求項1〜7のいずれか1つに記載の物標検出装置。
A filter processing unit that generates target data corresponding to the instantaneous data by applying a particle filter that assigns a plurality of particle data to the instantaneous data generated by the generation unit,
The estimation unit includes
When the degree of similarity between the direction of the vector based on the plurality of particle data and the movement direction is equal to or greater than a predetermined value, the movement direction is used as the final value of the estimation result. The target detection device according to any one of the above.
送信した電波が物標で反射した複数の反射点それぞれについて、前記物標の自車両への向きの相対速度を含む瞬時データを生成する生成工程と、
前記生成工程によって生成された前記瞬時データの前記相対速度と前記自車両の速度とに基づいて、対地速度を算出する算出工程と、
前記算出工程によって算出された前記対地速度に基づき、複数の前記瞬時データそれぞれに対応する前記自車両への向きの対地ベクトルにおける起点を揃え、当該対地ベクトルに対する垂線の交点から前記物標の移動向きを推定する推定工程と
を含むことを特徴とする物標検出方法。
For each of the plurality of reflection points of the transmitted radio wave reflected by the target, a generation step of generating instantaneous data including a relative speed of the direction of the target toward the own vehicle,
A calculating step of calculating a ground speed based on the relative speed of the instantaneous data generated by the generating step and the speed of the host vehicle;
Based on the ground speed calculated in the calculating step, the starting points in the ground vector of the direction to the host vehicle corresponding to each of the plurality of instantaneous data are aligned, and the moving direction of the target from the intersection of the perpendicular to the ground vector. An estimation step of estimating the target.
JP2018113047A 2018-06-13 2018-06-13 Target detection device Active JP7164327B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2018113047A JP7164327B2 (en) 2018-06-13 2018-06-13 Target detection device
US16/437,044 US11237260B2 (en) 2018-06-13 2019-06-11 Processing apparatus for estimating a movement direction of a target

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018113047A JP7164327B2 (en) 2018-06-13 2018-06-13 Target detection device

Publications (2)

Publication Number Publication Date
JP2019215275A true JP2019215275A (en) 2019-12-19
JP7164327B2 JP7164327B2 (en) 2022-11-01

Family

ID=68918602

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018113047A Active JP7164327B2 (en) 2018-06-13 2018-06-13 Target detection device

Country Status (1)

Country Link
JP (1) JP7164327B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022113472A1 (en) * 2020-11-30 2022-06-02 日立Astemo株式会社 Vehicle control device and vehicle control method
WO2023026920A1 (en) * 2021-08-26 2023-03-02 パナソニックIpマネジメント株式会社 Sensing device, processing device and data processing method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007333615A (en) * 2006-06-16 2007-12-27 Alpine Electronics Inc Radar measurement system
JP2008286582A (en) * 2007-05-16 2008-11-27 Fujitsu Ten Ltd Radar signal processing device and method
JP2013217837A (en) * 2012-04-11 2013-10-24 Furuno Electric Co Ltd Target motion estimating device, target motion estimating method, and radar device
WO2017207727A1 (en) * 2016-06-03 2017-12-07 Autoliv Development Ab Enhanced object detection and motion state estimation for a vehicle environment detection system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007333615A (en) * 2006-06-16 2007-12-27 Alpine Electronics Inc Radar measurement system
JP2008286582A (en) * 2007-05-16 2008-11-27 Fujitsu Ten Ltd Radar signal processing device and method
JP2013217837A (en) * 2012-04-11 2013-10-24 Furuno Electric Co Ltd Target motion estimating device, target motion estimating method, and radar device
WO2017207727A1 (en) * 2016-06-03 2017-12-07 Autoliv Development Ab Enhanced object detection and motion state estimation for a vehicle environment detection system
JP2019526781A (en) * 2016-06-03 2019-09-19 ヴィオニア スウェーデン エービー Improved object detection and motion state estimation for vehicle environment detection systems

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022113472A1 (en) * 2020-11-30 2022-06-02 日立Astemo株式会社 Vehicle control device and vehicle control method
JP7474352B2 (en) 2020-11-30 2024-04-24 日立Astemo株式会社 Vehicle control device and vehicle control method
WO2023026920A1 (en) * 2021-08-26 2023-03-02 パナソニックIpマネジメント株式会社 Sensing device, processing device and data processing method

Also Published As

Publication number Publication date
JP7164327B2 (en) 2022-11-01

Similar Documents

Publication Publication Date Title
US9869761B2 (en) Radar apparatus
JP6924046B2 (en) Radar device and elevation estimation method
JP6439728B2 (en) Mounted angle learning device
JP6993136B2 (en) Radar device and target detection method
JP2018115931A (en) Radar device and method for detecting target
JP7164327B2 (en) Target detection device
JP6924066B2 (en) Radar device and target detection method
JP7060441B2 (en) Radar device and target detection method
JP7244220B2 (en) Radar device and target data output method
JP6873009B2 (en) Radar device and target detection method
US11237260B2 (en) Processing apparatus for estimating a movement direction of a target
JP2020020607A (en) Radar device and axial deviation detection method
US11614533B2 (en) Radar signal processing device, radar device, and radar signal processing method
JP7111506B2 (en) Target detection device
Nadjiasngar et al. A performance model for target tracking with a radar network
JP6980979B2 (en) Radar device and target detection method
JP7181010B2 (en) Radar device and target detection method
WO2021075331A1 (en) Object tracking device
JP6157565B2 (en) Object detection device
JP6587784B2 (en) Radar signal processing apparatus and radar signal processing method
JP6929183B2 (en) Radar device and target detection method
JP7168353B2 (en) Radar device and target data allocation method
JP2019200079A (en) Target detector and target detection method
JP2019039750A (en) Rader system and method for detecting target
JP2018151327A (en) Radar device and method of combining bearings

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210429

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20211013

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220216

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20220217

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220322

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220523

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20221020

R150 Certificate of patent or registration of utility model

Ref document number: 7164327

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150