(57)【特許請求の範囲】
【請求項１】
レーダ装置（１００）を用いて少なくとも１つのレーダ標的（２００）の横方向速度成分を算出するための方法において、
同一に変調された送信信号を、所定数の送信素子を有する送信装置（１，１０）によって、所定の測定時間（Ｔｍｅａｓ）中に前記レーダ装置（１００）の所定の検出範囲内に周期的に送信するステップと、
前記レーダ標的（２００）で反射された少なくとも１つのレーダ受信信号を、所定数の受信素子を有する受信装置（２，２０）によって受信するステップと、
受信されたレーダ受信信号を、該受信されたレーダ受信信号のアナログデジタル変換器および評価装置（３０）に伝送するステップと、
送信素子および受信素子から成るそれぞれの組み合わせのためのデジタル測定値のそれぞれ１つの速度距離スペクトルを生成するための二次元フーリエ変換を実施するステップと、
前記速度距離スペクトルの振幅スペクトル内の所定の最高値を基に前記レーダ標的（２００）の少なくとも１つの標的反射を検出するステップと、
前記レーダ装置（１００）からの前記レーダ標的（２００）の距離、および前記レーダ装置（１００）に対して相対的な半径方向速度成分を、前記速度距離スペクトルから算出するステップと、
アンテナ（１０，２０）に対して相対的な前記レーダ標的（２００）の少なくとも１つの角度を決定するステップと、
前記横方向速度成分が算出されるべき前記レーダ標的（２００）を選択するステップと、
このような形式で選択された前記レーダ標的（２００）の標的反射の逆フーリエ変換を実施するステップと、
前記変換された測定値から前記レーダ標的（２００）の横方向速度成分を算出するステップと、
有する、少なくとも１つのレーダ標的の横方向速度成分を算出するための方法。
【請求項２】
前記同一に変調された送信信号が、ランプ状に周波数変調されたレーダ信号または非線形のランプ信号または周期的なパルス信号またはＯＦＤＭ信号である、請求項１に記載の方法。
【請求項３】
前記少なくとも１つのレーダ標的（２００）の前記横方向速度成分を、方位角方向および／または仰角方向で算出する、請求項１または２に記載の方法。
【請求項４】
物理的な所与性のモデリングのために、以下の数理モデルを使用する、
【数１】
0006821051.tif 000007
この数理モデルは次のパラメータを有している、
ｒ…        送受信アンテナとレーダ標的との間の距離
ｔ…        時間
Δｒバー…  測定時間に亘る距離変化の平均値
ｖｑ…      横方向速度（横方向の速度）
ｒ(ｔ)…    測定時間に亘る標的距離（ランプからランプ）
φ…        方位角および／または仰角
Δφ(ｔ)…  時間の経過に伴う方位角および／または仰角の変化
ｖｅ…      横方向速度から得られた半径方向の速度成分
ａｅ…      横方向速度から得られた半径方向の加速度成分
Ｔｍｅａｓ…  測定時間
ｂ…        補助パラメータ
請求項１〜３のいずれか一項に記載の方法。
【請求項５】
前記横方向速度成分の推定を、次の数学的関係を用いて品質関数（Ｑ）の最大化によって算出する、
【数２】
0006821051.tif 000008
この数学的関係は次のパラメータを有している、
Ｑ…      方位角および仰角における横方向速度成分の品質関数
ｖｑａ…    方位角における横方向速度成分
ｖｑｅ…    仰角における横方向速度成分
ｋ…      波数
ａｍ…    仮想の受信チャンネル（送信チャンネルと受信チャンネルとから成る組み合わせ）の数に相当する書き込みを有するベクトル
Ａｃ…    校正マトリックス
ｍ…      時間ｔｍの指数
ｒｍ…    時点ｔｍにおけるすべての受信チャンネルとの標的距離のベクトル
ωｄ…    推定されたドップラー角振動数
ｒバー…  測定区間内の平均的な標的間隔
ＴＲＲ…  ランプ信号間の時間的間隔
Δβｍ…  時点ｔｍにおける測定信号の位相位置の変化
ＸＴｍ…    実際に横方向速度が推定されるべき物理的標的の、できるだけ反射（ピーク）の信号成分だけを含有する、時点ｔｍにおける予め処理され選択された逆フーリエ変換測定信号を有する転置されたベクトル
請求項４に記載の方法。
【請求項６】
前記レーダ標的（２００）を算出するために、前記横方向速度を有する検出範囲を、前記検出範囲内の前記品質関数が凸状でありしかも前記検出範囲が完全にスキャニングされないように選択する、請求項５に記載の方法。
【請求項７】
標的反射信号の選択および再構成を、角度スペクトルを介して実施し、標的反射位置を方位角および／または仰角を介して互いに分離する、請求項１〜６のいずれか一項に記載の方法。
【請求項８】
前記方法を用途に応じて、選択された前記レーダ標的（２００）のために実施する、請求項１〜７のいずれか一項に記載の方法。
【請求項９】
トラッキングに基づく方法のために使用できるよう構成されていることを特徴とする、請求項１〜８のいずれか一項に記載の方法。
【請求項１０】
落ちてくる対象物を警告し、かつ／またはこの対象物との衝突を回避し、かつ／またはその結果発生する少なくとも事故の深刻さを最小にするよう構成されていることを特徴とする、請求項１〜９のいずれか一項に記載の方法。
【請求項１１】
レーダ装置（１００）を用いて少なくとも１つのレーダ標的（２００）の横方向速度成分を算出するための装置において、
送信装置（１，１０）を有していて、該送信装置（１，１０）が、同一に変調されたレーダ送信信号を所定の測定時間（Ｔｍｅａｓ）中に前記レーダ装置（１００）の所定の検出範囲内に周期的に送信することを実施するように構成された所定数の送信素子を備えており、
受信装置（２，２０）を有していて、該受信装置（２，２０）が、前記レーダ標的（２００）で反射された少なくとも１つのレーダ受信信号の受信を実施するように構成された所定数の受信素子を備えており、
処理装置（３０）を有していて、該処理装置（３０）は、受信されたレーダ受信信号のアナログデジタル変換を実施し、デジタル測定値の速度距離スペクトルを生成するための二次元フーリエ変換を実施し、前記速度距離スペクトルの振幅スペクトル内の所定の最高値に基づいて前記レーダ標的（２００）の少なくとも１つの標的反射の検出を実施し、アンテナ（１０，２０）からの前記レーダ標的（２００）の距離および前記アンテナ（１０，２０）に対して相対的な半径方向速度成分を前記速度距離スペクトルから算出することを実施し、前記アンテナ（１０，２０）に対して相対的な前記レーダ標的（２００）の少なくとも１つの角度の決定を実施するように、構成されており、
前記横方向速度成分が推定されるべき前記レーダ標的（２００）を除くすべての標的反射および妨害が抑制され、
このような形式で選択された前記レーダ標的（２００）の標的反射のために、ドップラー方向で逆フーリエ変換が実施され、変換された測定値から前記レーダ標的（２００）の横方向速度成分を算出する、
少なくとも１つのレーダ標的（２００）の横方向速度成分を算出するための装置。
【請求項１２】
落ちてくる対象物を警告し、かつ／または落ちてくる対象物との衝突を回避し、かつ／またはその結果発生する、少なくとも事故の深刻さを最小にするために構成されていることを特徴とする、請求項１１に記載の装置。
【請求項１３】
請求項１〜１０のいずれか一項に記載の方法を実行するためのプログラム。
【請求項１４】
請求項１〜１０のいずれか一項に記載の方法を実行するためのプログラムを記録した記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、レーダ標的の横方向相対速度成分を算出するための方法に関する。本発明はさらに、レーダ標的の横方向相対速度成分を算出するための装置に関する。
【背景技術】
【０００２】
従来技術では、例えば非特許文献１によるレーダシステムが公知である。図１は、公知のレーダ装置１００の原理図を示す。送信器１内で生成された変調されたレーダ信号が送信アンテナ１０を介して放射される。放射された電磁信号は次いで、場合によっては検知フィールド内に存在するレーダ標的２００（例えば、自動車、人間、柱、ガードレール、様々な物質間の移行部その他）で反射され、遅延時間Ｔ後に受信アンテナ２０を介して再び受信され、受信器２内で評価装置３０によってさらに処理される。
【０００３】
新型のレーダシステムは、送信信号として、例えば非特許文献２により公知であるいわゆる高速チャープ変調若しくはラピッドチャープ変調を使用する。この場合、図２に示されているように、例えば２０ｍｓ長さの測定時間Ｔｍｅａｓ＝ＭｘＴｒｒの間に、全部でＭ個の短いＦＭＣＷランプ（英語：frequency modulated continuous wave）が例えば１０μｓ〜１００μｓの時間Ｔｍｏｄで送信される。この場合、個別のランプＴｒｒの時間的な間隔は、同じ大きさであり、この場合、この時間的な間隔は、ランプ時間よりもやや大きいかまたは小さくてもよい。ランプは、時間的に非等間隔で配置されていない（図示せず）。例えば速度または距離の曖昧さを解像するために使用される、互いに入り組んだランプを有するレーダシステムは、例えば特許文献１、特許文献２または特許文献３により公知である。
【０００４】
結果としてまず、それぞれ個別のランプのための（変調）周波数を介して受信信号が得られる。次いで、それぞれのランプに所属する、通常はデジタル化された、測定された受信信号が逆フーリエ変換を介して時間レンジ内で変換される。しかしながら大抵はこのために、適切な窓掛けを有するデジタルフーリエ変換（ＤＦＴ）若しくは高速フーリエ変換（ＦＦＴ）が使用され、このようにして変換された領域が「うなり周波数範囲」と呼ばれる。さらに設けられたステップが、ランプからランプへのドップラー周波数範囲内のフーリエ変換である。このためにフーリエ変換が、個別ランプの対応する値に沿って（ドップラー方向若しくは速度方向で）実施される。２つの変換は、その連続を入れ替えるか、または二次元フーリエ変換とみなされてもよい。
【０００５】
続いて、典型的な別のステップ、例えば二次元スペクトルのピークの位置に基づいた、相対速度の半径方向成分および間隔の検出、推定および様々な誤差補正が行われる。この場合、ピークは標的反射に相当し、物理的標的（例えば自動車、人間、柱その他）は複数の標的反射を有していてよい。
【０００６】
このレーダシステムは、複数の送信および／または受信チャンネルを有するアンテナ（例えばパッチアンテナの個別のパッチによって実現される）、つまりアンテナアレイを有していて、追加的に標的反射の角度推定およびひいては空間内の３Ｄ標的位置の決定が実施され得る。このために、第１の信号処理ステップが、検出の前に、送信および受信チャンネルから成るそれぞれの組み合わせのために別個に実施される。次いで、個別のチャンネルの集合されたスペクトルを基に、角度推定が実施される。推定された標的（反射）パラメータは、様々な用途、例えばアダプティブクルーズコントロール（ＡＣＣ）、死角検出、自動非常ブレーキ機能その他のための、例えば次いで行われるトラッキング、クラスタリング、標的分類若しくはデータ融合のために使用される。
【０００７】
非特許文献３により、レーダ標的の横方向速度を推定するための方法が公知である。この場合、個別の剛体に割り当てられる複数の標的反射が組み合わせられることによって、横方向速度が決定される。これは可能である。何故ならば、相対速度の測定された半径方向成分の分布は空間内の反射の個別位置に依存するからである。このための前提は、相応の物理的標的のために三次元的に分布された多数の反射位置も測定され、互いに密接な関係にあるものとしてクラスタリングされ得る、ということである。しかしながら、クラスタリングにおいては、誤った割当ても発生し得る。
【０００８】
レーダ信号処理における公知のやり方は、レーダセンサと標的の個別反射位置との間の相対速度の半径方向成分若しくは長手方向成分を決定するだけである。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】独国特許出願公開第１０２０１２２２０８７９号明細書
【特許文献２】国際公開第２０１５／１９７２２９号
【特許文献３】国際公開第２０１５／１９７２２２号
【非特許文献】
【００１０】
【非特許文献１】「M. Skolnik,“Radar Handbook”, 3rd edition,2008（エム スコルニク、“レーダハンドブック”、第３版、２００８年）」
【非特許文献２】「Steffen Lutz, Daniel Ellenrieder, Thomas Walter, Robert Weigel:“On fast chirp Modulations and Compressed Sensing for Automotive Radar Applications”, International Radar Symposium, 2014（ステファン ルッツ、ダニエル エレンリーダー、トーマス ワルター、ロベルト ヴァイゲル、“自動車レーダアプリケーションのための高速チャープ変調および圧縮センシング”、インターナショナルレーダシンポジウム、２０１４）」
【非特許文献３】「D. Kellner et al. “Instantaneous Lateral Velocity Estimation of a Vehicle using Doppler Radar”, 16th International Conference on Information Fusion, Istanbul, Turkey, 2013（ディー ケルナーその他“ドップラーレーダを使用した自動車の瞬時横方向速度推定”、情報融合に関する第１６回国際会議、イスタンブール、トルコ、２０１３）」
【発明の概要】
【発明が解決しようとする課題】
【００１１】
本発明の課題は、レーダ標的の横方向速度を算出するための改善された方法を提供することである。
【課題を解決するための手段】
【００１２】
この課題は、第１の態様によれば、レーダ装置を用いて少なくとも１つのレーダ標的の横方向速度成分を算出するための方法によって解決され、この方法は、
同一に変調された送信信号を、所定数の送信素子を有する送信装置によって、所定の測定時間中にレーダ装置の所定の検出範囲内に周期的に送信するステップと、
レーダ標的で反射された少なくとも１つのレーダ受信信号を、所定数の受信素子を有する受信装置によって受信するステップと、
受信されたレーダ受信信号を、受信されたレーダ受信信号のアナログデジタル変換および評価装置に伝送するステップと、
送信素子および受信素子から成るそれぞれの組み合わせのためのデジタル測定値のそれぞれ１つの速度距離スペクトルを生成するための二次元フーリエ変換を実施するステップと、
速度距離スペクトルの振幅スペクトル内の所定の最高値を基にレーダ標的の少なくとも１つの標的反射を検出するステップと、
レーダ装置からのレーダ標的の距離、およびレーダ装置に対して相対的な半径方向速度成分を、速度距離スペクトルから算出するステップと、
アンテナに対して相対的なレーダ標的の少なくとも１つの角度を決定するステップと、
横方向速度成分が算出されるべきレーダ標的を選択するステップと、
このような形式で選択されたレーダ標的の標的反射の逆フーリエ変換を実施するステップと、
変換された測定値からレーダ標的の横方向速度成分を算出するステップと、
を有している。
【００１３】
好適にはこのような形式で、一度だけの測定によりレーダ標的の横方向の速度成分を測定することができる。好適には、これによって高価なトラッキング法を必要としないので、一度だけの測定によって完全な速度ベクトルが検出される。
【００１４】
この課題は、第２の態様によれば、レーダ装置によって少なくとも１つのレーダ標的の横方向速度成分を算出するための装置によって解決され、この装置は、
送信装置を有していて、この送信装置が、同一に変調されたレーダ送信信号を所定の測定時間中にレーダ装置の所定の検出範囲内に周期的に送信することを実施するように構成された所定数の送信素子を備えており、
受信装置を有していて、この受信装置が、レーダ標的で反射された少なくとも１つのレーダ受信信号の受信を実施するように構成された所定数の受信素子を備えており、
処理装置を有していて、この処理装置は、受信されたレーダ受信信号のアナログデジタル変換を実施し、デジタル測定値の速度距離スペクトルを生成するための二次元フーリエ変換を実施し、速度距離スペクトルの振幅スペクトル内の所定の最高値を基にレーダ標的の少なくとも１つの標的反射の検出を実施し、アンテナからのレーダ標的の距離およびアンテナに対して相対的な半径方向速度成分を、速度距離スペクトルから算出することを実施し、アンテナに対して相対的なレーダ標的の少なくとも１つの角度の決定を実施するように、構成されており、
この場合、横方向速度成分が推定されるべきレーダ標的を除くすべての標的反射および妨害が抑制され、
この場合、このような形式で選択されたレーダ標的の標的反射のために、逆フーリエ変換が実施され、変換された測定値からレーダ標的の横方向速度成分が算出される。
【００１５】
この方法の好適な発展形態は、従属請求項の対象である。
【００１６】
この方法の好適な発展形態は、同一に変調された送信信号が、ランプ状に周波数変調されたレーダ信号または非線形のランプ信号または周期的なパルス信号またはＯＦＤＭ信号であることを意図している。このような形式で、好適な様々なレーダ信号がこの方法を実施するために使用され得る。
【００１７】
この方法の別の好適な発展形態は、少なくとも１つのレーダ標的の横方向の速度成分が、方位角方向および／または仰角方向で算出されることが意図されている。これによって好適には、レーダ装置が有する検出特性が促進される。
【００１８】
この方法の別の好適な発展形態は、物理的な所与性のモデリングのために、以下の数理モデルが使用されることを意図している：
【００１９】
【数１】
0006821051.tif 000002
【００２０】
この数理モデルは次のパラメータを有している：
Δｒバー…  全測定時間に亘る距離変化の平均値
ｖｑ…      横方向速度（横方向の速度成分）
ｒ(ｔ)…    測定時間に亘る距離（ランプからランプ）
Δφ(ｔ)…  時間の経過に伴う方位角および仰角の変化
ｖｅ…      横方向速度から得られた半径方向の速度成分
ａｅ…      横方向速度から得られた半径方向の加速度成分
ｒ０…      時点ｔ＝０ｓにおける最小の標的距離。
これによって、横方向の速度成分の算出が、簡単な数理モデルに基づいて実施され得る。
【００２１】
この方法の別の好適な発展形態によれば、横方向速度成分の推定が、次の数学的関係を用いて品質関数の最大化によって算出されるようになっており：
【００２２】
【数２】
0006821051.tif 000003
【００２３】
この数学的関係は次のパラメータを有している：
Ｑ…    方位角および仰角における横方向速度成分の品質関数
ｖｑａ…  方位角における横方向速度成分
ｖｑｅ…  仰角における横方向速度成分
ｋ…    波数
ａｍ…    仮想の受信チャンネル（送信チャンネルと受信チャンネルとから成る組み合わせ）の数に相当する書き込みを有するベクトル
Ａｃ…    校正マトリックス
ｍ…      時間の指数
ｒｍ…    時点ｔｍにおけるすべての受信チャンネルに対する標的距離のベクトル
ωｄ…    推定されたドップラー角振動数
ｒバー…  測定区間内の平均的な標的間隔
ＴＲＲ…  ランプ信号間の時間間隔
Δβｍ…  時点ｔｍにおける測定信号の位相位置の変化
ＸＴｍ…    実際に横方向速度が推定されるべき物理的標的の、できるだけ反射（ピーク）の信号成分だけを含有する、時点ｔｍにおける予め処理され選択された逆フーリエ変換測定信号を有する転置されたベクトル。
【００２４】
このような形式で、簡単に実施することができる相関関係が提供され、この相関関係によって、レーダ標的２００の求められた相対速度成分および距離が測定信号の位相位置に換算される。
【００２５】
この方法の好適な発展形態によれば、標的反射信号の選択および再構成が、角度スペクトルを介して実施され、この場合、標的反射位置が方位角および／または仰角を介して互いに分離されるようになっている。このような形式で、距離速度空間内で複数の標的反射が同じ距離／速度位置に配置されている場合、レーダ標的の識別が実施される。
【００２６】
この方法の別の好適な発展形態によれば、レーダ標的を算出するために、横方向速度を有する検出範囲が、この検出範囲内の品質関数が凸状でありしかも検出範囲が完全にスキャニングされないように選択されるようになっている。このような形式で、完全な検出範囲を算出する必要はないので、局所的な最大「停滞」に留まることはない。この方法の効果的な実施は、このような形式で促進され、このことは、“Ｎｅｗｔｏｎ−Ｉｔｅｒａｔｉｏｎ（ニュートン反復法）”のような例えば公知の勾配法によって実施され得る。
【００２７】
この方法の別の好適な発展形態によれば、この方法は、用途に応じて選択されたレーダ標的のために実施されるようになっている。これによって、この方法は好適には、自動車のレーダシステムの多様な用途のために、例えばＡＣＣシステム、自動非常ブレーキ機能、死角検出その他のために使用され得る。
【００２８】
この方法の別の好適な発展形態によれば、この方法は、トラッキングに基づく方法のために使用されるようになっている。このような形式で、公知のトラッキングに基づく方法はより精確に実施されることができ、これによってレーダ標的の予測が改善されている。
【００２９】
本発明をその他の特徴および利点と共に複数の図面を用いて以下に詳しく説明する。この場合、明細書および図面に開示されたすべての特徴は、特許請求の範囲の引用とは無関係に、本発明の対象を形成する。図面は、特に本発明にとって必須の原理を説明するために役立つ。
【００３０】
開示された方法の特徴は、開示された対応する装置の特徴から同様に得られ、その逆でもある。このことは特に、方法に関する特徴、技術的な利点および実施例は、同様の形式で、装置に関する対応する特徴、技術的な利点および実施例から得られ、その逆でもある、ということである。
【図面の簡単な説明】
【００３１】
【図１】レーダシステムの原理的なブロック回路図である。
【図２】この方法に使用された例としてのランプ信号の変調周波数の経時変化である。
【図３】レーダ標的の横方向速度成分を算出するための信号処理の原理的なフローチャートである。
【図４】距離速度空間内の振幅スペクトルの複数のレーダ標的である。
【図５】レーダ標的の横方向速度の物理的関係の図である。
【図６】レーダ標的の変化するベクトルの図である。
【図７】レーダ標的の変化するパラメータの複数の経時変化である。
【図８】レーダ標的の横方向速度成分を算出するための別の信号処理の原理的なフローチャートである。
【図９】レーダ標的の窓掛け若しくは選択の図である。
【発明を実施するための形態】
【００３２】
本発明の課題は、特にレーダ標的の速度の横方向成分若しくはレーダ標的の速度の完全な速度ベクトルを、前もってのクラスタリング若しくはトラッキングなしに個別測定によって算出することである。
【００３３】
このために、図１のそれ自体公知のレーダ装置１００が使用され、この場合、送信装置１，１０から測定時間若しくは測定長さＴｍｅａｓで同一に変調された送信信号が検波範囲に送信される。この場合、同一に変調された送信信号は、ランプ状の周波数変調されたレーダ信号または非線形のランプ信号または周期的なパルス信号またはＯＦＤＭ信号であってよい。
【００３４】
横方向速度成分は、測定時間Ｔｍｅａｓ内の測定値に基づいて決定され、例えば後続のトラキング、クラスタリング、標的分類その他に追加的な情報として提供される。トラッキング中の測定毎の場合によっては誤った反射位置の割り当て、および反射位置の時間微分は、このような形式で好適には避けられる。
【００３５】
本発明によれば、横方向速度成分の追加的に得られた情報が、自動車の運転者支援システム、例えば自動非常ブレーキ機能の改善された機能性のために使用され得る。この場合、好適には、標的が例えば自車の車線外へ若しくは車線内に移動するかどうかが検知され得る。死角検出および／またはＡＣＣのためにも、追跡しようとする標的のための算出された横方向速度成分が使用されてよい。
【００３６】
任意に、提案された方法によって、落ちてくる対象物（例えば山の落石、橋からの落盤等）のための警告機能若しくは衝突回避機能を実現することも可能である。推定された横方向速度は、ここでは追加的にトラッキングのために使用されてよく、この場合、特に仰角方向で十分に広い検出フィールドを有する適切なアンテナ設計が必要となる。
【００３７】
提案された方法のために必要な信号処理の基本的なフローチャートは図３に示されている。ステップ３００でアナログからデジタルに変換された測定値を基に、ステップ３１０で実行された二次元フーリエ変換（２Ｄ−ＤＦＴ）に従って、それぞれのチャンネル毎に、つまりアンテナ（例えばパッチアレイ）１０，２０の送信素子および受信素子から成るそれぞれの組み合わせ毎に、うなり周波数ドップラー中の、若しくは距離速度空間内のうなり周波数（図示せず）内のドップラー成分の補正後の二次元スペクトルが、原則的に振幅スペクトルとして図４に示されているように得られる。図４は、単に例として全部で８つの標的反射を示しており、これらの標的反射は、すべてが同じ半径方向速度で異なる距離を保って距離速度空間内に配置されている。
【００３８】
標的反射のために、複数のチャンネルの個別スペクトルから算出された振幅スペクトル内の、ステップ３２０で検出されたピーク値（ピーク）が得られ、これらのピーク値の位置がステップ３３０で推定される。それぞれ検出しようとする反射位置に属する（例えばＭＩＭＯ−アンテナシステムの）複数のチャンネルの測定値に基づいて、次にステップ３４０で、当該の標的反射のための方位角方向および／または仰角方向での角度推定が実施される。
【００３９】
従って、ステップ３００乃至３４０は、それ自体公知の信号処理の標準ステップを有している。
【００４０】
本発明によれば、単数または複数の選択された検出された標的反射のために、追加的にさらに、横方向速度成分が推定されるか若しくは概ね算出される。このために、ステップ３５０で個別スペクトル内のそれぞれの標的反射若しくはピークに属していない値がゼロに設定され（図９に詳しく説明されているように“窓掛け”）、これによって結局、標的反射を表す測定値の選択が実施される。選択的に、このために複素値のピーク振幅および位置の精確な推定も実施され得る。次いで、この再構成は、ステップ３６０で、フーリエ変換とは逆のドップラー方向で行われるので、予め処理された複素信号値ｘｍの連続が得られ、これは、それぞれ時間的に相次いで連続するランプに対応する。
【００４１】
このような形式で得られかつ予め処理された測定値ｘｍを基に、ステップ３７０で、方位角方向および／または仰角方向における横方向速度成分が推定される。このために、測定時間Ｔｍｅａｓ＝ＭｘＴｒｒ（つまり、ランプｍからランプｍ＋１まで）に亘って変化する距離が、以下の簡単な数理モデルによって表される。
【００４２】
【数３】
0006821051.tif 000004
【００４３】
この場合、平均的な距離変化、半径方向の速度および加速度は、以下のように規定されている。
【００４４】
【数４】
0006821051.tif 000005
【００４５】
この数理モデルは次のパラメータを有している：
ｂ…        見やすさを改善するための補助パラメータ
ｒ０…      測定時間の中央での（時点ｔ＝０での）レーダ標的の最小距離
Δｒ(ｔ)…  レーダ標的の時間的な距離変化
ｖｑ…      方位角方向または仰角方向における横方向速度
ｖｅ…      横方向速度から得られた、角度に依存する相対速度の半径方向成分
ａｅ…      横方向速度から得られた、角度に依存する半径方向の加速度成分
Δφ…      方位角方向または仰角方向の若しくは両方向から成る一次結合の角度変化。
【００４６】
以下に図５、図６および図７に示されているように、横方向速度は、時間的な距離変化を介しても、また角度変化を介しても測定信号に作用する。
【００４７】
図５は、アンテナ１０，２０と、時点ｔ＝０におけるレーダ標的２００の間隔ｒ０と、時間ｔ後に横方向速度成分ｖｑに基づいて進んだ走行距離ｖｑ×ｔと、それから得られた、半径方向の速度成分ｖｅに基づく走行距離ｖｅ×ｔとを示す。
【００４８】
図６は、アンテナ１０，２０と、それから得られた、アンテナ１０，２０に対するレーダ標的２００の間隔ｒｍを有するランプ信号の中央時点ｔｍとを示す。
【００４９】
図７は、単に例として４つのグラフを示しており、上から下に向かって、横方向速度成分ｖｑに基づく半径方向の速度成分ｖｅの経時変化、それから得られた間隔変化Δｒ、それから得られた角度変化Δφ、それから得られたレーダ標的２００の加速度ａｅの変化を示す。
【００５０】
横方向速度成分ｖｑの推定のために、それぞれのランプｍの時点でｖｑに依存する標的位置が算出され、それから得られた個別の受信チャンネル内の位相変化Δβｍが求められ、この場合、ＭＩＭＯシステムにおいて、ここでは実際の送受信チャンネルの組み合わせから得られる仮想の受信チャンネルを意味する。
【００５１】
各時間ステップｍのためにベクトルａｍ（ｖｑ，ａ；ｖｑ，ｅ）が得られ、このベクトルのエレメントは受信チャンネルを表す。すべてのＭベクトルの集合が１つのマトリックスをもたらし、このマトリックスの複数の列が複数の受信チャンネルを表し、このマトリックスの複数の行が複数のランプ時点を表す。マトリックスはそれぞれ、方位角方向ｖｑ，ａおよび仰角方向ｖｑ，ｅにおける横方向速度の調査されるべき組み合わせのために算出され、予め処理された測定値ベクトルｘｍが乗算される。次いで、横方向速度の推定が、ｖｑ，ａおよび／またはｖｑ，ｅを介して品質関数Ｑ（ｖｑ，ａ；ｖｑ，ｅ）の最大化によって以下の数理モデルに従って行われる。
【００５２】
【数５】
0006821051.tif 000006
【００５３】
この数理モデルは次のパラメータを有している：
Ｑ…    方位角および仰角における横方向速度成分の品質関数
ｖｑａ…  方位角における横方向速度成分
ｖｑｅ…  仰角における横方向速度成分
ｋ…      波数
ａｍ…    仮想の受信チャンネル（送信チャンネルおよび受信チャンネルの組み合わせ）の数に相当する書き込みを有するベクトル
Ａｃ…    校正マトリックス
ｍ…      ランプ指数若しくはランプからランプへの時間の指数
ｒｍ…    時点ｔｍにおけるすべての受信チャンネルとの標的距離のベクトル
ωｄ…    推定されたドップラー角振動数
ｒバー…  測定区間内の平均的な標的間隔
ＴＲＲ…    ランプ信号間の時間的間隔
Δφｍ…    時点ｔｍにおける方位角および／または仰角の変化
ＸＴｍ…    実際に横方向速度が推定されるべき物理的標的の、できるだけ反射（ピーク）の信号成分だけを含有する、時点ｔｍにおける予め処理され選択された逆フーリエ変換測定信号を有する転置されたベクトル。
【００５４】
図８は、この方法の原理的なフローチャートを示し、この原理的なフローチャートは、図３のフローチャートに概ね相当するが、図８ではステップ３１０で、複数のチャンネル（若しくは仮想の受信チャンネル）のために二次元フーリエ変換が実施される点で異なっている。ステップ３３１でさらに、角度推定のための測定値の選択が示されている。このオプション的な追加ステップは、角度推定３４０のためのコストを低減するために、つまり後続の使用のために重要な標的に割り当てられている、距離速度スペクトル内のピーク位置のためにだけ角度推定を実施するために特に有効である。
【００５５】
仮想の受信アンテナの位置が互いに概ね等距離で配置されていれば、これによってサイドローブが例えばＦＦＴ窓掛け（例えば“Kaiser-Fenster「カイザー窓」”、“Dolph-Chebyshev-Fenster「ドルフ・チェビシェフ窓」”）によって抑制され、それによって角度スペクトル内の高いダイナミックスが促進され、ひいては横方向速度３５０の推定のためにレーダ標的の選択が、追加的に角度スペクトル内で簡単に可能となる。
【００５６】
これが行われない場合、角度を介してのみ切り離し可能な複数のレーダ標的における横方向速度成分ｖｑの推定は不都合な影響を受ける。何故ならば、これらのレーダ標的は同じ距離速度セル内に位置しているからである。
【００５７】
図３および図８に図示された信号処理ステップは、レーダ装置１００の所属の制御プログラムにより、コンピュータ装置（例えばマイクロコントローラ、ＤＳＰその他）のソフトウエアで実行され、またハードウエア（例えばＦＰＧＡ，ＡＳＩＣその他）でも実現され得る。
【００５８】
図９は、実際に横方向速度成分ｖｑが算出されるべきｋ番目の標的反射Ｒの、点として示された複数の反射位置を有する距離速度スペクトル内のレーダ標的２００の選択の原理を示す。推定された距離ｒｋにおけるベクトルｘｍの値だけが使用される。ｖｋだけピーク幅の外に位置するベクトルｘｍの値は、窓掛け若しくはゼロ設定によって抑制される。その結果、それにより、本当に実際の横方向速度が算出されるべき標的反射Ｒだけが残る。ベクトルｘｍのエレメントは、速度方向ｖにおけるスペクトルの指数ｍに対する仮想の角度スペクトルの値であるか、若しくはｘｍ＝ＩＦＦＴ（ｘｍ）だけ逆フーリエ変換後の時間的なランプ指数ｍに対する仮想の角度スペクトルの値である。窓掛けは、一次元若しくは二次元の角度スペクトル内でアナログに行われる。提案された方法は、このような形式で窓掛けされたそれぞれの標的反射Ｒのためにシーケンシャルに実施される。ステップ３６０で逆フーリエ変換によってドップラー方向でベクトルｘ＝ＩＦＦＴ（ｘ）が形成され、これは、次にステップ３７０で行われる横方向速度の推定のための入力値として用いられる。
【００５９】
本発明による方法の好適な適用は、落ちてくる対象物（例えば橋から投げ落とされた対象物または山の落石）を警告し、かつ／または前記対象物との衝突を回避し、かつ／または少なくともその結果発生する事故の深刻さを最小にするシステムを実現する点にある。
【００６０】
当業者は、本発明の核心から逸脱することなしに、この方法および装置のここに開示されていない実施例も実現可能である。
【符号の説明】
【００６１】
１，１０  送信装置
２，２０  受信装置
１０，２０  アンテナ
３０  評価装置、処理装置
１００  レーダ装置
２００  レーダ標的
３００，３１０，３２０，３３０，３３１，３４０，３５０，３６０，３７０  ステップ
【図１】
0006821051.tif 000009
【図２】
0006821051.tif 000010
【図３】
0006821051.tif 000011
【図４】
0006821051.tif 000012
【図５】
0006821051.tif 000013
【図６】
0006821051.tif 000014
【図７】
0006821051.tif 000015
【図８】
0006821051.tif 000016
【図９】
0006821051.tif 000017