JP6821051B2 - レーダ標的の横方向相対速度成分を算出するための方法および装置 - Google Patents

レーダ標的の横方向相対速度成分を算出するための方法および装置 Download PDF

Info

Publication number
JP6821051B2
JP6821051B2 JP2019548725A JP2019548725A JP6821051B2 JP 6821051 B2 JP6821051 B2 JP 6821051B2 JP 2019548725 A JP2019548725 A JP 2019548725A JP 2019548725 A JP2019548725 A JP 2019548725A JP 6821051 B2 JP6821051 B2 JP 6821051B2
Authority
JP
Japan
Prior art keywords
radar
target
velocity component
signal
radar target
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.)
Active
Application number
JP2019548725A
Other languages
English (en)
Other versions
JP2020509390A (ja
Inventor
ブローシェ,トーマス
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.)
Robert Bosch GmbH
Original Assignee
Robert Bosch GmbH
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 Robert Bosch GmbH filed Critical Robert Bosch GmbH
Publication of JP2020509390A publication Critical patent/JP2020509390A/ja
Application granted granted Critical
Publication of JP6821051B2 publication Critical patent/JP6821051B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/66Radar-tracking systems; Analogous systems
    • G01S13/72Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
    • G01S13/723Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar by using numerical data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S13/589Velocity or trajectory determination systems; Sense-of-movement determination systems measuring the velocity vector
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06Systems determining position data of a target
    • G01S13/42Simultaneous measurement of distance and other co-ordinates
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/285Receivers
    • G01S7/295Means for transforming co-ordinates or for evaluating data, e.g. using computers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/35Details of non-pulse systems
    • G01S7/352Receivers
    • G01S7/354Extracting wanted echo-signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/35Details of non-pulse systems
    • G01S7/352Receivers
    • G01S7/356Receivers involving particularities of FFT processing

Landscapes

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

Description

本発明は、レーダ標的の横方向相対速度成分を算出するための方法に関する。本発明はさらに、レーダ標的の横方向相対速度成分を算出するための装置に関する。
従来技術では、例えば非特許文献1によるレーダシステムが公知である。図1は、公知のレーダ装置100の原理図を示す。送信器1内で生成された変調されたレーダ信号が送信アンテナ10を介して放射される。放射された電磁信号は次いで、場合によっては検知フィールド内に存在するレーダ標的200(例えば、自動車、人間、柱、ガードレール、様々な物質間の移行部その他)で反射され、遅延時間T後に受信アンテナ20を介して再び受信され、受信器2内で評価装置30によってさらに処理される。
新型のレーダシステムは、送信信号として、例えば非特許文献2により公知であるいわゆる高速チャープ変調若しくはラピッドチャープ変調を使用する。この場合、図2に示されているように、例えば20ms長さの測定時間Tmeas=MxTrrの間に、全部でM個の短いFMCWランプ(英語:frequency modulated continuous wave)が例えば10μs〜100μsの時間Tmodで送信される。この場合、個別のランプTrrの時間的な間隔は、同じ大きさであり、この場合、この時間的な間隔は、ランプ時間よりもやや大きいかまたは小さくてもよい。ランプは、時間的に非等間隔で配置されていない(図示せず)。例えば速度または距離の曖昧さを解像するために使用される、互いに入り組んだランプを有するレーダシステムは、例えば特許文献1、特許文献2または特許文献3により公知である。
結果としてまず、それぞれ個別のランプのための(変調)周波数を介して受信信号が得られる。次いで、それぞれのランプに所属する、通常はデジタル化された、測定された受信信号が逆フーリエ変換を介して時間レンジ内で変換される。しかしながら大抵はこのために、適切な窓掛けを有するデジタルフーリエ変換(DFT)若しくは高速フーリエ変換(FFT)が使用され、このようにして変換された領域が「うなり周波数範囲」と呼ばれる。さらに設けられたステップが、ランプからランプへのドップラー周波数範囲内のフーリエ変換である。このためにフーリエ変換が、個別ランプの対応する値に沿って(ドップラー方向若しくは速度方向で)実施される。2つの変換は、その連続を入れ替えるか、または二次元フーリエ変換とみなされてもよい。
続いて、典型的な別のステップ、例えば二次元スペクトルのピークの位置に基づいた、相対速度の半径方向成分および間隔の検出、推定および様々な誤差補正が行われる。この場合、ピークは標的反射に相当し、物理的標的(例えば自動車、人間、柱その他)は複数の標的反射を有していてよい。
このレーダシステムは、複数の送信および/または受信チャンネルを有するアンテナ(例えばパッチアンテナの個別のパッチによって実現される)、つまりアンテナアレイを有していて、追加的に標的反射の角度推定およびひいては空間内の3D標的位置の決定が実施され得る。このために、第1の信号処理ステップが、検出の前に、送信および受信チャンネルから成るそれぞれの組み合わせのために別個に実施される。次いで、個別のチャンネルの集合されたスペクトルを基に、角度推定が実施される。推定された標的(反射)パラメータは、様々な用途、例えばアダプティブクルーズコントロール(ACC)、死角検出、自動非常ブレーキ機能その他のための、例えば次いで行われるトラッキング、クラスタリング、標的分類若しくはデータ融合のために使用される。
非特許文献3により、レーダ標的の横方向速度を推定するための方法が公知である。この場合、個別の剛体に割り当てられる複数の標的反射が組み合わせられることによって、横方向速度が決定される。これは可能である。何故ならば、相対速度の測定された半径方向成分の分布は空間内の反射の個別位置に依存するからである。このための前提は、相応の物理的標的のために三次元的に分布された多数の反射位置も測定され、互いに密接な関係にあるものとしてクラスタリングされ得る、ということである。しかしながら、クラスタリングにおいては、誤った割当ても発生し得る。
レーダ信号処理における公知のやり方は、レーダセンサと標的の個別反射位置との間の相対速度の半径方向成分若しくは長手方向成分を決定するだけである。
独国特許出願公開第102012220879号明細書 国際公開第2015/197229号 国際公開第2015/197222号
本発明の課題は、レーダ標的の横方向速度を算出するための改善された方法を提供することである。
この課題は、第1の態様によれば、レーダ装置を用いて少なくとも1つのレーダ標的の横方向速度成分を算出するための方法によって解決され、この方法は、
同一に変調された送信信号を、所定数の送信素子を有する送信装置によって、所定の測定時間中にレーダ装置の所定の検出範囲内に周期的に送信するステップと、
レーダ標的で反射された少なくとも1つのレーダ受信信号を、所定数の受信素子を有する受信装置によって受信するステップと、
受信されたレーダ受信信号を、受信されたレーダ受信信号のアナログデジタル変換および評価装置に伝送するステップと、
送信素子および受信素子から成るそれぞれの組み合わせのためのデジタル測定値のそれぞれ1つの速度距離スペクトルを生成するための二次元フーリエ変換を実施するステップと、
速度距離スペクトルの振幅スペクトル内の所定の最高値を基にレーダ標的の少なくとも1つの標的反射を検出するステップと、
レーダ装置からのレーダ標的の距離、およびレーダ装置に対して相対的な半径方向速度成分を、速度距離スペクトルから算出するステップと、
アンテナに対して相対的なレーダ標的の少なくとも1つの角度を決定するステップと、
横方向速度成分が算出されるべきレーダ標的を選択するステップと、
このような形式で選択されたレーダ標的の標的反射の逆フーリエ変換を実施するステップと、
変換された測定値からレーダ標的の横方向速度成分を算出するステップと、
を有している。
好適にはこのような形式で、一度だけの測定によりレーダ標的の横方向の速度成分を測定することができる。好適には、これによって高価なトラッキング法を必要としないので、一度だけの測定によって完全な速度ベクトルが検出される。
この課題は、第2の態様によれば、レーダ装置によって少なくとも1つのレーダ標的の横方向速度成分を算出するための装置によって解決され、この装置は、
送信装置を有していて、この送信装置が、同一に変調されたレーダ送信信号を所定の測定時間中にレーダ装置の所定の検出範囲内に周期的に送信することを実施するように構成された所定数の送信素子を備えており、
受信装置を有していて、この受信装置が、レーダ標的で反射された少なくとも1つのレーダ受信信号の受信を実施するように構成された所定数の受信素子を備えており、
処理装置を有していて、この処理装置は、受信されたレーダ受信信号のアナログデジタル変換を実施し、デジタル測定値の速度距離スペクトルを生成するための二次元フーリエ変換を実施し、速度距離スペクトルの振幅スペクトル内の所定の最高値を基にレーダ標的の少なくとも1つの標的反射の検出を実施し、アンテナからのレーダ標的の距離およびアンテナに対して相対的な半径方向速度成分を、速度距離スペクトルから算出することを実施し、アンテナに対して相対的なレーダ標的の少なくとも1つの角度の決定を実施するように、構成されており、
この場合、横方向速度成分が推定されるべきレーダ標的を除くすべての標的反射および妨害が抑制され、
この場合、このような形式で選択されたレーダ標的の標的反射のために、逆フーリエ変換が実施され、変換された測定値からレーダ標的の横方向速度成分が算出される。
この方法の好適な発展形態は、従属請求項の対象である。
この方法の好適な発展形態は、同一に変調された送信信号が、ランプ状に周波数変調されたレーダ信号または非線形のランプ信号または周期的なパルス信号またはOFDM信号であることを意図している。このような形式で、好適な様々なレーダ信号がこの方法を実施するために使用され得る。
この方法の別の好適な発展形態は、少なくとも1つのレーダ標的の横方向の速度成分が、方位角方向および/または仰角方向で算出されることが意図されている。これによって好適には、レーダ装置が有する検出特性が促進される。
この方法の別の好適な発展形態は、物理的な所与性のモデリングのために、以下の数理モデルが使用されることを意図している:
Figure 0006821051
この数理モデルは次のパラメータを有している:
Δrバー… 全測定時間に亘る距離変化の平均値
… 横方向速度(横方向の速度成分)
r(t)… 測定時間に亘る距離(ランプからランプ)
Δφ(t)… 時間の経過に伴う方位角および仰角の変化
… 横方向速度から得られた半径方向の速度成分
… 横方向速度から得られた半径方向の加速度成分
… 時点t=0sにおける最小の標的距離。
これによって、横方向の速度成分の算出が、簡単な数理モデルに基づいて実施され得る。
この方法の別の好適な発展形態によれば、横方向速度成分の推定が、次の数学的関係を用いて品質関数の最大化によって算出されるようになっており:
Figure 0006821051
この数学的関係は次のパラメータを有している:
Q… 方位角および仰角における横方向速度成分の品質関数
qa… 方位角における横方向速度成分
qe… 仰角における横方向速度成分
k… 波数
… 仮想の受信チャンネル(送信チャンネルと受信チャンネルとから成る組み合わせ)の数に相当する書き込みを有するベクトル
… 校正マトリックス
m… 時間の指数
… 時点tにおけるすべての受信チャンネルに対する標的距離のベクトル
ω… 推定されたドップラー角振動数
rバー… 測定区間内の平均的な標的間隔
RR… ランプ信号間の時間間隔
Δβ… 時点tにおける測定信号の位相位置の変化
実際に横方向速度が推定されるべき物理的標的の、できるだけ反射(ピーク)の信号成分だけを含有する、時点tにおける予め処理され選択された逆フーリエ変換測定信号を有する転置されたベクトル。
このような形式で、簡単に実施することができる相関関係が提供され、この相関関係によって、レーダ標的200の求められた相対速度成分および距離が測定信号の位相位置に換算される。
この方法の好適な発展形態によれば、標的反射信号の選択および再構成が、角度スペクトルを介して実施され、この場合、標的反射位置が方位角および/または仰角を介して互いに分離されるようになっている。このような形式で、距離速度空間内で複数の標的反射が同じ距離/速度位置に配置されている場合、レーダ標的の識別が実施される。
この方法の別の好適な発展形態によれば、レーダ標的を算出するために、横方向速度を有する検出範囲が、この検出範囲内の品質関数が凸状でありしかも検出範囲が完全にスキャニングされないように選択されるようになっている。このような形式で、完全な検出範囲を算出する必要はないので、局所的な最大「停滞」に留まることはない。この方法の効果的な実施は、このような形式で促進され、このことは、“Newton−Iteration(ニュートン反復法)”のような例えば公知の勾配法によって実施され得る。
この方法の別の好適な発展形態によれば、この方法は、用途に応じて選択されたレーダ標的のために実施されるようになっている。これによって、この方法は好適には、自動車のレーダシステムの多様な用途のために、例えばACCシステム、自動非常ブレーキ機能、死角検出その他のために使用され得る。
この方法の別の好適な発展形態によれば、この方法は、トラッキングに基づく方法のために使用されるようになっている。このような形式で、公知のトラッキングに基づく方法はより精確に実施されることができ、これによってレーダ標的の予測が改善されている。
本発明をその他の特徴および利点と共に複数の図面を用いて以下に詳しく説明する。この場合、明細書および図面に開示されたすべての特徴は、特許請求の範囲の引用とは無関係に、本発明の対象を形成する。図面は、特に本発明にとって必須の原理を説明するために役立つ。
開示された方法の特徴は、開示された対応する装置の特徴から同様に得られ、その逆でもある。このことは特に、方法に関する特徴、技術的な利点および実施例は、同様の形式で、装置に関する対応する特徴、技術的な利点および実施例から得られ、その逆でもある、ということである。
レーダシステムの原理的なブロック回路図である。 この方法に使用された例としてのランプ信号の変調周波数の経時変化である。 レーダ標的の横方向速度成分を算出するための信号処理の原理的なフローチャートである。 距離速度空間内の振幅スペクトルの複数のレーダ標的である。 レーダ標的の横方向速度の物理的関係の図である。 レーダ標的の変化するベクトルの図である。 レーダ標的の変化するパラメータの複数の経時変化である。 レーダ標的の横方向速度成分を算出するための別の信号処理の原理的なフローチャートである。 レーダ標的の窓掛け若しくは選択の図である。
本発明の課題は、特にレーダ標的の速度の横方向成分若しくはレーダ標的の速度の完全な速度ベクトルを、前もってのクラスタリング若しくはトラッキングなしに個別測定によって算出することである。
このために、図1のそれ自体公知のレーダ装置100が使用され、この場合、送信装置1,10から測定時間若しくは測定長さTmeasで同一に変調された送信信号が検波範囲に送信される。この場合、同一に変調された送信信号は、ランプ状の周波数変調されたレーダ信号または非線形のランプ信号または周期的なパルス信号またはOFDM信号であってよい。
横方向速度成分は、測定時間Tmeas内の測定値に基づいて決定され、例えば後続のトラキング、クラスタリング、標的分類その他に追加的な情報として提供される。トラッキング中の測定毎の場合によっては誤った反射位置の割り当て、および反射位置の時間微分は、このような形式で好適には避けられる。
本発明によれば、横方向速度成分の追加的に得られた情報が、自動車の運転者支援システム、例えば自動非常ブレーキ機能の改善された機能性のために使用され得る。この場合、好適には、標的が例えば自車の車線外へ若しくは車線内に移動するかどうかが検知され得る。死角検出および/またはACCのためにも、追跡しようとする標的のための算出された横方向速度成分が使用されてよい。
任意に、提案された方法によって、落ちてくる対象物(例えば山の落石、橋からの落盤等)のための警告機能若しくは衝突回避機能を実現することも可能である。推定された横方向速度は、ここでは追加的にトラッキングのために使用されてよく、この場合、特に仰角方向で十分に広い検出フィールドを有する適切なアンテナ設計が必要となる。
提案された方法のために必要な信号処理の基本的なフローチャートは図3に示されている。ステップ300でアナログからデジタルに変換された測定値を基に、ステップ310で実行された二次元フーリエ変換(2D−DFT)に従って、それぞれのチャンネル毎に、つまりアンテナ(例えばパッチアレイ)10,20の送信素子および受信素子から成るそれぞれの組み合わせ毎に、うなり周波数ドップラー中の、若しくは距離速度空間内のうなり周波数(図示せず)内のドップラー成分の補正後の二次元スペクトルが、原則的に振幅スペクトルとして図4に示されているように得られる。図4は、単に例として全部で8つの標的反射を示しており、これらの標的反射は、すべてが同じ半径方向速度で異なる距離を保って距離速度空間内に配置されている。
標的反射のために、複数のチャンネルの個別スペクトルから算出された振幅スペクトル内の、ステップ320で検出されたピーク値(ピーク)が得られ、これらのピーク値の位置がステップ330で推定される。それぞれ検出しようとする反射位置に属する(例えばMIMO−アンテナシステムの)複数のチャンネルの測定値に基づいて、次にステップ340で、当該の標的反射のための方位角方向および/または仰角方向での角度推定が実施される。
従って、ステップ300乃至340は、それ自体公知の信号処理の標準ステップを有している。
本発明によれば、単数または複数の選択された検出された標的反射のために、追加的にさらに、横方向速度成分が推定されるか若しくは概ね算出される。このために、ステップ350で個別スペクトル内のそれぞれの標的反射若しくはピークに属していない値がゼロに設定され(図9に詳しく説明されているように“窓掛け”)、これによって結局、標的反射を表す測定値の選択が実施される。選択的に、このために複素値のピーク振幅および位置の精確な推定も実施され得る。次いで、この再構成は、ステップ360で、フーリエ変換とは逆のドップラー方向で行われるので、予め処理された複素信号値xの連続が得られ、これは、それぞれ時間的に相次いで連続するランプに対応する。
このような形式で得られかつ予め処理された測定値xを基に、ステップ370で、方位角方向および/または仰角方向における横方向速度成分が推定される。このために、測定時間Tmeas=MxTrr(つまり、ランプmからランプm+1まで)に亘って変化する距離が、以下の簡単な数理モデルによって表される。
Figure 0006821051
この場合、平均的な距離変化、半径方向の速度および加速度は、以下のように規定されている。
Figure 0006821051
この数理モデルは次のパラメータを有している:
b… 見やすさを改善するための補助パラメータ
… 測定時間の中央での(時点t=0での)レーダ標的の最小距離
Δr(t)… レーダ標的の時間的な距離変化
… 方位角方向または仰角方向における横方向速度
… 横方向速度から得られた、角度に依存する相対速度の半径方向成分
… 横方向速度から得られた、角度に依存する半径方向の加速度成分
Δφ… 方位角方向または仰角方向の若しくは両方向から成る一次結合の角度変化。
以下に図5、図6および図7に示されているように、横方向速度は、時間的な距離変化を介しても、また角度変化を介しても測定信号に作用する。
図5は、アンテナ10,20と、時点t=0におけるレーダ標的200の間隔rと、時間t後に横方向速度成分vに基づいて進んだ走行距離v×tと、それから得られた、半径方向の速度成分vに基づく走行距離v×tとを示す。
図6は、アンテナ10,20と、それから得られた、アンテナ10,20に対するレーダ標的200の間隔rを有するランプ信号の中央時点tとを示す。
図7は、単に例として4つのグラフを示しており、上から下に向かって、横方向速度成分vに基づく半径方向の速度成分vの経時変化、それから得られた間隔変化Δr、それから得られた角度変化Δφ、それから得られたレーダ標的200の加速度aの変化を示す。
横方向速度成分vの推定のために、それぞれのランプmの時点でvに依存する標的位置が算出され、それから得られた個別の受信チャンネル内の位相変化Δβmが求められ、この場合、MIMOシステムにおいて、ここでは実際の送受信チャンネルの組み合わせから得られる仮想の受信チャンネルを意味する。
各時間ステップmのためにベクトルa(vq,a;vq,e)が得られ、このベクトルのエレメントは受信チャンネルを表す。すべてのMベクトルの集合が1つのマトリックスをもたらし、このマトリックスの複数の列が複数の受信チャンネルを表し、このマトリックスの複数の行が複数のランプ時点を表す。マトリックスはそれぞれ、方位角方向vq,aおよび仰角方向vq,eにおける横方向速度の調査されるべき組み合わせのために算出され、予め処理された測定値ベクトルxが乗算される。次いで、横方向速度の推定が、vq,aおよび/またはvq,eを介して品質関数Q(vq,a;vq,e)の最大化によって以下の数理モデルに従って行われる。
Figure 0006821051
この数理モデルは次のパラメータを有している:
Q… 方位角および仰角における横方向速度成分の品質関数
qa… 方位角における横方向速度成分
qe… 仰角における横方向速度成分
k… 波数
… 仮想の受信チャンネル(送信チャンネルおよび受信チャンネルの組み合わせ)の数に相当する書き込みを有するベクトル
… 校正マトリックス
m… ランプ指数若しくはランプからランプへの時間の指数
… 時点tにおけるすべての受信チャンネルとの標的距離のベクトル
ω… 推定されたドップラー角振動数
rバー… 測定区間内の平均的な標的間隔
RR… ランプ信号間の時間的間隔
Δφ… 時点tにおける方位角および/または仰角の変化
実際に横方向速度が推定されるべき物理的標的の、できるだけ反射(ピーク)の信号成分だけを含有する、時点tにおける予め処理され選択された逆フーリエ変換測定信号を有する転置されたベクトル。
図8は、この方法の原理的なフローチャートを示し、この原理的なフローチャートは、図3のフローチャートに概ね相当するが、図8ではステップ310で、複数のチャンネル(若しくは仮想の受信チャンネル)のために二次元フーリエ変換が実施される点で異なっている。ステップ331でさらに、角度推定のための測定値の選択が示されている。このオプション的な追加ステップは、角度推定340のためのコストを低減するために、つまり後続の使用のために重要な標的に割り当てられている、距離速度スペクトル内のピーク位置のためにだけ角度推定を実施するために特に有効である。
仮想の受信アンテナの位置が互いに概ね等距離で配置されていれば、これによってサイドローブが例えばFFT窓掛け(例えば“Kaiser-Fenster「カイザー窓」”、“Dolph-Chebyshev-Fenster「ドルフ・チェビシェフ窓」”)によって抑制され、それによって角度スペクトル内の高いダイナミックスが促進され、ひいては横方向速度350の推定のためにレーダ標的の選択が、追加的に角度スペクトル内で簡単に可能となる。
これが行われない場合、角度を介してのみ切り離し可能な複数のレーダ標的における横方向速度成分vの推定は不都合な影響を受ける。何故ならば、これらのレーダ標的は同じ距離速度セル内に位置しているからである。
図3および図8に図示された信号処理ステップは、レーダ装置100の所属の制御プログラムにより、コンピュータ装置(例えばマイクロコントローラ、DSPその他)のソフトウエアで実行され、またハードウエア(例えばFPGA,ASICその他)でも実現され得る。
図9は、実際に横方向速度成分vが算出されるべきk番目の標的反射Rの、点として示された複数の反射位置を有する距離速度スペクトル内のレーダ標的200の選択の原理を示す。推定された距離rにおけるベクトルxの値だけが使用される。vだけピーク幅の外に位置するベクトルxの値は、窓掛け若しくはゼロ設定によって抑制される。その結果、それにより、本当に実際の横方向速度が算出されるべき標的反射Rだけが残る。ベクトルxのエレメントは、速度方向vにおけるスペクトルの指数mに対する仮想の角度スペクトルの値であるか、若しくはx=IFFT(x)だけ逆フーリエ変換後の時間的なランプ指数mに対する仮想の角度スペクトルの値である。窓掛けは、一次元若しくは二次元の角度スペクトル内でアナログに行われる。提案された方法は、このような形式で窓掛けされたそれぞれの標的反射Rのためにシーケンシャルに実施される。ステップ360で逆フーリエ変換によってドップラー方向でベクトルx=IFFT(x)が形成され、これは、次にステップ370で行われる横方向速度の推定のための入力値として用いられる。
本発明による方法の好適な適用は、落ちてくる対象物(例えば橋から投げ落とされた対象物または山の落石)を警告し、かつ/または前記対象物との衝突を回避し、かつ/または少なくともその結果発生する事故の深刻さを最小にするシステムを実現する点にある。
当業者は、本発明の核心から逸脱することなしに、この方法および装置のここに開示されていない実施例も実現可能である。
1,10 送信装置
2,20 受信装置
10,20 アンテナ
30 評価装置、処理装置
100 レーダ装置
200 レーダ標的
300,310,320,330,331,340,350,360,370 ステップ

Claims (14)

  1. レーダ装置(100)を用いて少なくとも1つのレーダ標的(200)の横方向速度成分を算出するための方法において、
    同一に変調された送信信号を、所定数の送信素子を有する送信装置(1,10)によって、所定の測定時間(Tmeas)中に前記レーダ装置(100)の所定の検出範囲内に周期的に送信するステップと、
    前記レーダ標的(200)で反射された少なくとも1つのレーダ受信信号を、所定数の受信素子を有する受信装置(2,20)によって受信するステップと、
    受信されたレーダ受信信号を、該受信されたレーダ受信信号のアナログデジタル変換器および評価装置(30)に伝送するステップと、
    送信素子および受信素子から成るそれぞれの組み合わせのためのデジタル測定値のそれぞれ1つの速度距離スペクトルを生成するための二次元フーリエ変換を実施するステップと、
    前記速度距離スペクトルの振幅スペクトル内の所定の最高値を基に前記レーダ標的(200)の少なくとも1つの標的反射を検出するステップと、
    前記レーダ装置(100)からの前記レーダ標的(200)の距離、および前記レーダ装置(100)に対して相対的な半径方向速度成分を、前記速度距離スペクトルから算出するステップと、
    アンテナ(10,20)に対して相対的な前記レーダ標的(200)の少なくとも1つの角度を決定するステップと、
    前記横方向速度成分が算出されるべき前記レーダ標的(200)を選択するステップと、
    このような形式で選択された前記レーダ標的(200)の標的反射の逆フーリエ変換を実施するステップと、
    前記変換された測定値から前記レーダ標的(200)の横方向速度成分を算出するステップと、
    有する、少なくとも1つのレーダ標的の横方向速度成分を算出するための方法。
  2. 前記同一に変調された送信信号が、ランプ状に周波数変調されたレーダ信号または非線形のランプ信号または周期的なパルス信号またはOFDM信号である、請求項1記載の方法。
  3. 前記少なくとも1つのレーダ標的(200)の前記横方向速度成分を、方位角方向および/または仰角方向で算出する、請求項1または2記載の方法。
  4. 物理的な所与性のモデリングのために、以下の数理モデルを使用する、
    Figure 0006821051
    この数理モデルは次のパラメータを有している
    r… 送受信アンテナとレーダ標的との間の距離
    t… 時間
    Δrバー… 測定時間に亘る距離変化の平均値
    … 横方向速度(横方向の速度)
    (t)… 測定時間に亘る標的距離(ランプからランプ)
    φ… 方位角および/または仰角
    Δφ(t)… 時間の経過に伴う方位角および/または仰角の変化
    … 横方向速度から得られた半径方向の速度成分
    … 横方向速度から得られた半径方向の加速度成分
    meas… 測定時間
    b… 補助パラメー
    請求項1〜3のいずれか一項に記載の方法。
  5. 前記横方向速度成分の推定を、次の数学的関係を用いて品質関数(Q)の最大化によって算出する、
    Figure 0006821051
    この数学的関係は次のパラメータを有している
    Q… 方位角および仰角における横方向速度成分の品質関数
    qa… 方位角における横方向速度成分
    qe… 仰角における横方向速度成分
    k… 波数
    … 仮想の受信チャンネル(送信チャンネルと受信チャンネルとから成る組み合わせ)の数に相当する書き込みを有するベクトル
    … 校正マトリックス
    m… 時間tの指数
    … 時点tにおけるすべての受信チャンネルとの標的距離のベクトル
    ω… 推定されたドップラー角振動数
    rバー… 測定区間内の平均的な標的間隔
    RR… ランプ信号間の時間的間隔
    Δβ… 時点tにおける測定信号の位相位置の変化
    実際に横方向速度が推定されるべき物理的標的の、できるだけ反射(ピーク)の信号成分だけを含有する、時点tにおける予め処理され選択された逆フーリエ変換測定信号を有する転置されたベクト
    請求項4記載の方法。
  6. 前記レーダ標的(200)を算出するために、前記横方向速度を有する検出範囲を、前記検出範囲内の前記品質関数が凸状でありしかも前記検出範囲が完全にスキャニングされないように選択する、請求項5記載の方法。
  7. 標的反射信号の選択および再構成を、角度スペクトルを介して実施し、標的反射位置を方位角および/または仰角を介して互いに分離する、請求項1〜6のいずれか一項に記載の方法。
  8. 前記方法を用途に応じて、選択された前記レーダ標的(200)のために実施する、請求項1〜7のいずれか一項に記載の方法。
  9. ラッキングに基づく方法のために使用できよう構成されていることを特徴とする、請求項1〜8のいずれか一項に記載の方法。
  10. 落ちてくる対象物を警告し、かつ/またはこの対象物との衝突を回避し、かつ/またはその結果発生する少なくとも事故の深刻さを最小にするよう構成されていることを特徴とする、請求項1〜9のいずれか一項に記載の方法。
  11. レーダ装置(100)を用いて少なくとも1つのレーダ標的(200)の横方向速度成分を算出するための装置において、
    送信装置(1,10)を有していて、該送信装置(1,10)が、同一に変調されたレーダ送信信号を所定の測定時間(Tmeas)中に前記レーダ装置(100)の所定の検出範囲内に周期的に送信することを実施するように構成された所定数の送信素子を備えており、
    受信装置(2,20)を有していて、該受信装置(2,20)が、前記レーダ標的(200)で反射された少なくとも1つのレーダ受信信号の受信を実施するように構成された所定数の受信素子を備えており、
    処理装置(30)を有していて、該処理装置(30)は、受信されたレーダ受信信号のアナログデジタル変換を実施し、デジタル測定値の速度距離スペクトルを生成するための二次元フーリエ変換を実施し、前記速度距離スペクトルの振幅スペクトル内の所定の最高値に基づいて前記レーダ標的(200)の少なくとも1つの標的反射の検出を実施し、アンテナ(10,20)からの前記レーダ標的(200)の距離および前記アンテナ(10,20)に対して相対的な半径方向速度成分を前記速度距離スペクトルから算出することを実施し、前記アンテナ(10,20)に対して相対的な前記レーダ標的(200)の少なくとも1つの角度の決定を実施するように、構成されており、
    前記横方向速度成分が推定されるべき前記レーダ標的(200)を除くすべての標的反射および妨害が抑制され、
    このような形式で選択された前記レーダ標的(200)の標的反射のために、ドップラー方向で逆フーリエ変換が実施され、変換された測定値から前記レーダ標的(200)の横方向速度成分を算出する、
    少なくとも1つのレーダ標的(200)の横方向速度成分を算出するための装置。
  12. ちてくる対象物を警告し、かつ/または落ちてくる対象物との衝突を回避し、かつ/またはその結果発生する、少なくとも事故の深刻さを最小にするために構成されていることを特徴とする、請求項11記載の装置。
  13. 請求項1〜10のいずれか一項に記載の方法を実行するためのプログラム。
  14. 請求項1〜10のいずれか一項に記載の方法を実行するためのプログラムを記録した記録媒体。
JP2019548725A 2017-03-17 2018-01-25 レーダ標的の横方向相対速度成分を算出するための方法および装置 Active JP6821051B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
DE102017204495.0 2017-03-17
DE102017204495.0A DE102017204495A1 (de) 2017-03-17 2017-03-17 Verfahren und Vorrichtung zum Ermitteln von transversalen Relativgeschwindigkeitskomponenten von Radarzielen
PCT/EP2018/051796 WO2018166683A1 (de) 2017-03-17 2018-01-25 Verfahren und vorrichtung zum ermitteln von transversalen relativgeschwindigkeitskomponenten von radarzielen

Publications (2)

Publication Number Publication Date
JP2020509390A JP2020509390A (ja) 2020-03-26
JP6821051B2 true JP6821051B2 (ja) 2021-01-27

Family

ID=61027743

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019548725A Active JP6821051B2 (ja) 2017-03-17 2018-01-25 レーダ標的の横方向相対速度成分を算出するための方法および装置

Country Status (7)

Country Link
US (1) US11249180B2 (ja)
EP (1) EP3596488B1 (ja)
JP (1) JP6821051B2 (ja)
KR (1) KR102445130B1 (ja)
CN (1) CN110431437B (ja)
DE (1) DE102017204495A1 (ja)
WO (1) WO2018166683A1 (ja)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102018200755A1 (de) * 2018-01-18 2019-07-18 Robert Bosch Gmbh Verfahren und Vorrichtung zum Plausibilisieren einer Querbewegung
US11175382B2 (en) * 2019-06-14 2021-11-16 GM Global Technology Operations LLC Elevation angle estimation in horizontal antenna array with doppler and velocity measurements
CN113033586B (zh) * 2019-12-24 2024-04-16 大富科技(安徽)股份有限公司 目标识别方法及设备
CN113030890B (zh) * 2019-12-24 2023-11-21 大富科技(安徽)股份有限公司 基于车载雷达的目标识别方法及设备
CN111458703A (zh) * 2020-01-13 2020-07-28 武汉大学 一种测量多目标横向速度的方法及***
US11137488B1 (en) * 2020-03-10 2021-10-05 Nokia Technologies Oy Radar excitation signals for wireless communications system
CN111580060B (zh) 2020-04-21 2022-12-13 北京航空航天大学 目标姿态识别的方法、装置和电子设备
KR20210136631A (ko) * 2020-05-08 2021-11-17 주식회사 만도모빌리티솔루션즈 차량용 레이더의 수직 장착 오정렬 감지 장치, 방법 및 그를 포함하는 레이더 장치
WO2022020995A1 (zh) * 2020-07-27 2022-02-03 华为技术有限公司 一种信号处理方法、装置以及存储介质

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10100413A1 (de) * 2001-01-08 2002-07-11 Bosch Gmbh Robert Verfahren und Vorrichtung zur Schätzung von Bewegungsparametern von Zielen
DE102012220879A1 (de) 2012-11-15 2014-05-15 Robert Bosch Gmbh Rapid-Chirps-FMCW-Radar
DE102014212280A1 (de) 2014-06-26 2015-12-31 Robert Bosch Gmbh Radarmessverfahren
DE102013210256A1 (de) * 2013-06-03 2014-12-04 Robert Bosch Gmbh Interferenzunterdrückung bei einem fmcw-radar
US9594159B2 (en) * 2013-07-15 2017-03-14 Texas Instruments Incorporated 2-D object detection in radar applications
JP6364986B2 (ja) * 2014-06-13 2018-08-01 株式会社デンソー レーダ装置
DE102014212390A1 (de) 2014-06-27 2015-12-31 Robert Bosch Gmbh Verfahren zur Objektortung mit einem FMCW-Radar
CN106443671A (zh) * 2016-08-30 2017-02-22 西安电子科技大学 基于调频连续波的sar雷达动目标检测与成像方法

Also Published As

Publication number Publication date
KR102445130B1 (ko) 2022-09-21
EP3596488A1 (de) 2020-01-22
US11249180B2 (en) 2022-02-15
CN110431437B (zh) 2023-09-22
JP2020509390A (ja) 2020-03-26
DE102017204495A1 (de) 2018-09-20
US20210132212A1 (en) 2021-05-06
EP3596488B1 (de) 2022-03-09
WO2018166683A1 (de) 2018-09-20
KR20190125452A (ko) 2019-11-06
CN110431437A (zh) 2019-11-08

Similar Documents

Publication Publication Date Title
JP6821051B2 (ja) レーダ標的の横方向相対速度成分を算出するための方法および装置
US10234541B2 (en) FMCW radar device
US11280882B2 (en) Method and radar device for ascertaining radial relative acceleration of at least one target
US9897685B2 (en) Method for detecting interference in a received signal of a radar sensor, driver assistance device and motor vehicle
JP7173735B2 (ja) レーダ装置及び信号処理方法
US10473760B2 (en) Radar device and vertical axis-misalignment detecting method
CN111630406A (zh) 使用机动车辆雷达***的不同传输器的相位编码fmcw斜坡进行改善目标分离的方法
JP7019059B2 (ja) レーダセンサによる点目標の横速度またはデカルト速度の推定
JP6316390B1 (ja) 物体検知装置および物体検知方法
EP4350389A1 (en) Method for estimating target fuzzy speed
US11150338B2 (en) Azimuth calculation device and azimuth calculation method
US11971472B2 (en) Method for determining an ego-velocity estimated value and an angle estimated value of targets
EP3480624A1 (en) Detection of parking row orientation
US11808884B2 (en) Adaptive high point neighborhood range
US20230228840A1 (en) Device and method with radar signal processing
EP4336214A1 (en) Method and apparatus with radar signal processing
EP4206730A1 (en) Method and apparatus with radar signal processing
JP7274271B2 (ja) レーダ装置および物標判別方法
JP2024501466A (ja) 自動車レーダーにおける高速ターゲット検出のための動き補償

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190906

RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20200615

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200924

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20201009

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20201111

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210105

R150 Certificate of patent or registration of utility model

Ref document number: 6821051

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250