レーダ装置
【課題】反射波の受信信号から反射元の種類を判別可能な技術を提供すること。
【解決手段】レーダ装置の信号処理部は、反射波のパワー（受信電力）Ｐと、このパワーＰに対応する反射波の受信時点から観測周期Ｔｓｗ＝Ｔｍ遡った時点で受信された反射波のパワーＰｂと、に基づき、パワーＰｂに対するパワーＰの変化量Ｙ＝ｌｏｇ（Ｐ）−ｌｏｇ（Ｐｂ）を算出する（Ｓ３１０）。一方、当該レーダ装置は、判別対象とする反射元の種類毎に、レーダ波が当該種類の反射元で反射して到来する反射波から算出される変化量Ｙの確率分布を記憶する。そして、算出された変化量Ｙと、上記反射元種類毎の確率分布とから、今回算出された変化量Ｙに対応するレーダ波の反射元の種類を判別する。即ち、反射元種類毎に、当該変化量Ｙが得られる確率Ｐｒを算出し（Ｓ３２０）、確率Ｐｒの最も高い種類を、今回のレーダ波の反射元であると判別する（Ｓ３３０）。
【解決手段】レーダ装置の信号処理部は、反射波のパワー（受信電力）Ｐと、このパワーＰに対応する反射波の受信時点から観測周期Ｔｓｗ＝Ｔｍ遡った時点で受信された反射波のパワーＰｂと、に基づき、パワーＰｂに対するパワーＰの変化量Ｙ＝ｌｏｇ（Ｐ）−ｌｏｇ（Ｐｂ）を算出する（Ｓ３１０）。一方、当該レーダ装置は、判別対象とする反射元の種類毎に、レーダ波が当該種類の反射元で反射して到来する反射波から算出される変化量Ｙの確率分布を記憶する。そして、算出された変化量Ｙと、上記反射元種類毎の確率分布とから、今回算出された変化量Ｙに対応するレーダ波の反射元の種類を判別する。即ち、反射元種類毎に、当該変化量Ｙが得られる確率Ｐｒを算出し（Ｓ３２０）、確率Ｐｒの最も高い種類を、今回のレーダ波の反射元であると判別する（Ｓ３３０）。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、レーダ装置に関する。
【背景技術】
【０００２】
従来、レーダ波を発射して、このレーダ波が前方物体に反射して戻ってくる反射波を受信し、送信波に対する受信波の遅延量から前方物体までの距離を検出し、送信波に対する受信波の周波数変化量から前方物体との相対速度を検出するレーダ装置が知られている（例えば特許文献１参照）。この種のレーダ装置は、例えば車両に搭載されて、車両制御に使用される。
【０００３】
この他、レーダ装置としては、物標の存在有無を受信信号の振幅と閾値との比較により判定するものが知られている。この手法では、クラッタにより物標の存在を誤判定する確率（誤警報確率）を一定に抑えるために、例えば、レーダ波の受信信号を対数増幅し、この対数増幅後の受信信号を時間的に前後する受信信号の対数増幅後の平均値で減算する（非特許文献１参照）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１０−０３２３１４号公報
【非特許文献】
【０００５】
【非特許文献１】関根松夫著、「レーダ信号処理技術」、社団法人 電子情報通信学会、１９９１年９月、ｐ．９６−１０３
【発明の概要】
【発明が解決しようとする課題】
【０００６】
ところで、従来知られるレーダ装置では、前方物体までの距離や前方物体の方位、前方物体との相対速度を検出することはできるものの、レーダ波を反射した前方物体の種類までを特定することはできない。従来技術によれば、受信波が物標からの反射により生じたものであるのか、それともクラッタであるのかを判別することができる程度である。
【０００７】
本発明は、こうした問題に鑑みなされたものであり、反射波の受信信号から、反射元の物体の種類を判別可能な技術を提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記目的を達成するためになされた本発明のレーダ装置は、レーダ波を発射し、当該レーダ波の反射波を受信する送受信手段と、送受信手段が受信した反射波の受信信号に基づき、反射波の信号強度を検出する検出手段と、を備えるレーダ装置であって、更に、変化量算出手段、確率分布記憶手段、及び、判別手段、を備えるものである。
【０００９】
このレーダ装置が備える変化量算出手段は、検出手段により検出された反射波の信号強度Ｐと、当該信号強度Ｐに対応する反射波の受信時点から所定の観測周期Ｔｓｗ遡った時点で送受信手段により受信された反射波の検出手段により検出された信号強度Ｐｂと、に基づき、信号強度Ｐｂに対する信号強度Ｐの変化量Ｙを算出する。変化量Ｙとしては、例えば信号強度Ｐ，Ｐｂの比Ｐ／Ｐｂを挙げることができ、信号強度Ｐ，Ｐｂを対数増幅する場合には、変化量Ｙとして、対数増幅後の信号強度ｌｏｇ（Ｐ）及びｌｏｇ（Ｐｂ）の差分ｌｏｇ（Ｐ）−ｌｏｇ（Ｐｂ）を算出することができる。
【００１０】
一方、確率分布記憶手段は、判別対象とする反射元の種類毎に、レーダ波が当該種類の反射元で反射して到来する反射波から算出される変化量Ｙの確率分布を記憶する。例えば、確率分布記憶手段は、上記確率分布を、累積分布関数や確率密度関数等の関数情報、又は、これら関数に代わるテーブル情報の形態で記憶する。
【００１１】
そして、判別手段は、確率分布記憶手段が記憶する上記種類毎の確率分布に基づき、変化量算出手段で算出された変化量Ｙに対応するレーダ波の反射元の種類を判別する。この他、このレーダ装置では、変化量Ｙからレーダ波の反射元の種類を判別可能とするため、判別対象とする各種類の反射元にて発生する揺らぎの特徴が変化量Ｙに反映されるような関係で、レーダ波の送信波長λ及び観測周期Ｔｓｗが定められている。
【００１２】
反射波の位相は、反射面の位置に応じて変化するため、反射面が揺らげば（振動すれば）、その揺らぎによる反射面の位置変化が反射波の位相変化に現れることになる。一方、反射元は、微小な複数の反射面の集合にて構成されるものと解釈することができ、その反射面の揺らぎ方は、反射元が金属物のような表面が硬いものなのか、人のような表面が柔らかいものなのか、各部位（例えば腕や脚）が変位するものなのか等によって異なる。そして、揺らぎ方が異なれば、変化量Ｙの確率分布に差異が生じる。本発明では、このような現象を利用し、上記変化量Ｙに基づいて、反射元の種類を判別する。
【００１３】
但し、観測周期Ｔｓｗの間に生じる各反射面から到来する反射波成分の位相変化量Δφがπ以上となると、各反射面の反射波成分が混合する上記送受信手段での反射波の受信信号においては、その振幅がレイリー分布に従うようになり、反射面の揺らぎを観測することが難しい。換言すれば、判別対象とする各種類の反射元にて発生する揺らぎの特徴が変化量Ｙに反映されなくなり、変化量Ｙに揺らぎ方による差異が表れないため、変化量Ｙに基づいて、反射元の種類を判別することが難しい。
【００１４】
そこで、本発明では、レーダ波の送信波長λ及び観測周期Ｔｓｗを、判別対象とする各種類の反射元にて発生する揺らぎの特徴が変化量Ｙに反映される関係となるように定めることで、揺らぎを、変化量Ｙにより観測することができるようにしている。
【００１５】
尚、観測周期Ｔｓｗ間の位相変化量Δφは、当然のごとく観測周期Ｔｓｗとレーダ波の送信波長λに依存する。更に、位相変化量Δφは、反射元の揺らぎの大きさ（反射面の変位幅）に依存する。よって、送信波長λ及び観測周期Ｔｓｗについては、判別対象の各種類の反射元での揺らぎの大きさを考慮して、判別対象とする各種類の反射元にて発生する揺らぎの特徴が変化量Ｙに反映される関係となるように定める。具体的には、判別対象とする複数種類の反射元のいずれにおいても、その反射元での各反射面から反射する反射波の位相変化量Δφがπを超えないように、送信波長λ及び観測周期Ｔｓｗを定めることで、上記関係を満足するように送信波長λ及び観測周期Ｔｓｗを定めることができる。
【００１６】
このように構成された本発明のレーダ装置によれば、反射波の受信信号から、反射元の物体の種類を判別することができる。従って、本発明のレーダ装置を用いれば、前方物体の種類を判別し、この判別結果に基づき、従来では実現することのできなかった処理を行うことができる。例えば、前方物体の種類に応じた処理の切り替えを行うことができる。従って、本発明によれば、従来よりも優れたレーダ装置を提供することができる。
【００１７】
尚、確率分布記憶手段は、上記種類毎の確率分布として少なくとも、反射元が金属物であるときの確率分布、及び反射元が人であるときの確率分布を記憶した構成にされるとよい。そして、判別手段は、これらの確率分布に基づき、上記反射元の種類として、少なくとも金属物及び人の夫々を判別可能な構成にされるとよい（請求項２）。
【００１８】
このようにレーダ装置を構成すれば、例えばレーダ装置による判別結果を車両制御に応用する場合に、前方物体が金属物（車両やガードレール等）であるのか人であるのかに応じて、車両制御の内容を切り替えることができ、危険回避等に大変役立つ。特に、従来技術では衝突回避が不可能となったときに、人を回避してガードレールに衝突するように車両制御を行うといったことは不可能であったが、本発明のレーダ装置を用いれば、そのような車両制御が実現可能となる。
【００１９】
また、反射元が人であるときの確率分布としては、人が静止しているときの確率分布及び人が移動しているときの確率分布の夫々を確率分布記憶手段に記憶させて、これらの確率分布に基づいて、静止している人及び移動している人を、上記反射元の種類として判別できるように上記判別手段を構成すると一層好ましい。このようにレーダ装置を構成すれば、反射元の種類についての判別結果を危険回避等に一層役立てることができる。
【００２０】
また、判別手段は、次のようにして反射元の種類を判別する構成にすることができる。即ち、判別手段は、判別対象とする反射元の種類毎に、確率分布記憶手段が記憶する当該種類の確率分布に基づき、上記変化量算出手段で算出された変化量Ｙの発生確率であって当該変化量Ｙに対応するレーダ波の反射元が仮に当該種類の反射元であるときの発生確率を導出し、導出した発生確率が最大である反射元の種類を、変化量算出手段で算出された変化量Ｙに対応する反射元の種類として判別する構成にすることができる（請求項３）。
【００２１】
この他、車載型のレーダ装置では、自車両周辺の道路環境を加味して反射元の種類を判別するように上記判別手段を構成すると一層好ましい。
即ち、車載型のレーダ装置には、自車両周辺の道路環境を推定する環境推定手段を設けて、判別手段は、判別対象とする反射元の種類毎に、確率分布記憶手段が記憶する当該種類の確率分布と、環境推定手段により推定された道路環境と、に基づき、変化量算出手段で算出された変化量Ｙに対応するレーダ波の反射元が当該種類の反射元である確度を算出し、反射元の種類毎の確度に基づき、変化量算出手段で算出された変化量Ｙに対応するレーダ波の反射元の種類を判別する構成にされると好ましい（請求項４）。
【００２２】
各種類の反射元が存在する可能性は、道路環境によって異なる。換言すれば、道路環境を指標にすれば、各種類の反射元が存在する可能性をある程度、事前予測することができる。例えば、車両が、街中を走行しているのか、郊外を走行しているのか、高速道路を走行しているのかによって、車両周辺に人が存在する可能性は大きく変化する。従って、このように道路環境を加味して、反射元の種類を判別するようにレーダ装置を構成すると、一層高精度に反射元の種類を判別することができる。
【００２３】
尚、環境推定手段による道路環境の推定結果を用いる場合、判別手段は、具体的に次のように構成することができる。即ち、判別手段は、判別対象とする反射元の種類毎に、確率分布記憶手段が記憶する当該種類の確率分布に基づき、変化量算出手段で算出された変化量Ｙの発生確率であってレーダ波の反射元が当該種類の反射元であると仮定した場合での発生確率を導出し、この導出した発生確率に環境推定手段により推定された道路環境に対応した補正を加えることで、変化量算出手段で算出された変化量Ｙに対応するレーダ波の反射元が当該種類の反射元である確度を算出する構成にすることができる（請求項５）。例えば、上記導出した発生確率に道路環境に対応した補正係数を掛けて、変化量Ｙに対応するレーダ波の反射元が上記種類の反射元である確度を算出するといった具合である。
【００２４】
また、環境推定手段は、各区域の道路環境を表すデータを記憶する道路環境データベースから直接又は間接的にデータ取得可能な構成であることを前提として、次のように構成することができる。即ち、環境推定手段は、道路環境データベースから自車両の現在位置に対応する区域の道路環境を表すデータを取得し、この取得データが表す道路環境を、自車両周辺の道路環境であると推定する構成にすることができる（請求項６）。
【００２５】
例えば、環境推定手段は、上記道路環境データベースとしての機能を有する道路地図データベースから、ＧＰＳ受信機等の位置検出器にて検出される自車両の現在位置に対応する区域の道路環境を表すデータ（例えば、街中、郊外、高速道路等の道路環境のカテゴリを表すデータ）を取得し、この取得データに基づき、自車両周辺の道路環境を推定する構成することができる。車両には、道路地図データベースやＧＰＳ受信機等の位置検出器が備えられている例も多いため、環境推定手段をこのように構成すれば、車内の機器を利用して、高精度に自車両周辺の道路環境を推定することができる。
【００２６】
但し、このような手法では、道路環境データベースをレーダ装置や車両内に設けたり、外部の道路環境データベースからデータを取得するための通信インタフェースをレーダ装置に設けたりする必要があるため、廉価に道路環境を加味した反射元種類の判別動作を実現するのが難しい。
【００２７】
そこで、環境推定手段は、自車両の走行速度に基づき、道路環境を推定する構成にされてもよい（請求項７）。交差点の多い環境において、車両は低速で走行することが多く、交差点の少ない環境においては、高速に走行することが多くなる。また、交差点のない高速道路等では、車両は更に高速に走行する。よって、車両の走行速度に基づけば、現在の道路環境が、どのような種類の道路環境であるのかを推定することができ、道路環境データベースを必要とせず、廉価に道路環境を加味した反射元種類の判別動作を実現することができる。
【００２８】
また、本発明は、ＦＭＣＷ（Frequency Modulated Continuous Wave）方式にてレーダ波の送受信を行うＦＭＣＷレーダ装置や、ＳＦＷ（Stepped Frequency Waveforms）方式にてレーダ波の送受信を行うＳＦＷレーダ装置に適用することができる（請求項８）。この種のレーダ装置によれば、レーダ波の反射元までの距離、反射元の方位、反射元の移動速度等を検出することができるので、本発明を適用して、更に反射元の種類を判別できるようにレーダ装置を構成すれば、その判別結果を、例えば当該レーダ装置を搭載する車両の走行制御等に役立てることができる。
【００２９】
また、本発明を、ＦＭＣＷ方式やＳＦＷ方式にて、レーダ波を送信アンテナから発射し、当該レーダ波の反射波をアレーアンテナで受信する送受信手段を有したＦＭＣＷレーダ装置やＳＦＷレーダ装置に適用する場合、当該レーダ装置には、アレーアンテナを構成する各アンテナ素子からの受信信号にレーダ波の送信信号を混合して、アンテナ素子毎のビート信号を生成するビート信号生成手段と、ビート信号生成手段により生成されたアンテナ素子毎のビート信号に基づき、レーダ波の反射元までの距離及び当該反射元の方位を検出する信号解析手段とを設け、更に上記検出手段、変化量算出手段、及び、判別手段を次のように構成することができる。
【００３０】
即ち、検出手段は、ビート信号生成手段により生成されたアンテナ素子毎のビート信号に基づき、距離及び方位の異なる上記反射元毎に、当該反射元から到来した反射波の信号強度を検出し、変化量算出手段は、上記反射元毎に、当該反射元から到来した反射波の上記検出手段により検出された信号強度Ｐと、当該信号強度Ｐに対応する反射波の受信時点から観測周期Ｔｓｗ遡った時点で当該反射元から到来し送受信手段により受信された反射波の上記検出手段により検出された信号強度Ｐｂと、に基づき、信号強度Ｐｂに対する信号強度Ｐの変化量Ｙを算出し、判別手段は、上記反射元毎に、当該反射元の種類を、上記変化量算出手段で算出された変化量Ｙと、確率分布記憶手段が記憶する上記種類毎の確率分布とに基づき判別する構成にすることができる（請求項９）。
【００３１】
このレーダ装置によれば、前方に存在する複数の物体の位置（距離及び方位）並びに種類を特定することができるので、レーダ装置周辺の状況を詳細に把握することができる。従って、この判別結果を例えば車両制御に利用すれば、前方に位置する物体の位置及び種類に応じて、適切な車両制御を実現でき、危険回避等に大変役立つ。
【図面の簡単な説明】
【００３２】
【図１】レーダ装置１の構成を表すブロック図である。
【図２】レーダ波の送受信号Ｓｓ，Ｓｒ及びビート信号ＢＴを示したグラフ（ａ）並びに、ビート信号ＢＴのパワースペクトルを示したグラフ（ｂ）である。
【図３】前方物体方位とＭＵＳＩＣスペクトルとの対応関係を示した説明図である。
【図４】確率分布テーブルの構成を表す図である。
【図５】信号処理部３０が実行するメインルーチンを表すフローチャートである。
【図６】信号処理部３０が実行する物標推定処理を表すフローチャートである。
【図７】第一変形例のレーダ装置２の構成を表す図である。
【図８】第一変形例の物標推定処理の一部を表すフローチャートである。
【図９】第二変形例の物標推定処理の一部を表すフローチャートである。
【発明を実施するための形態】
【００３３】
以下に本発明の実施例について、図面と共に説明する。
本実施例のレーダ装置１は、ＦＭＣＷ方式によりレーダ波を送受信する車載型のＦＭＣＷレーダ装置であり、図１に示すように、周波数が時間に対して直線的に漸次増減するミリ波帯の高周波信号を生成する発振器１１と、発振器１１が生成する高周波信号を増幅する増幅器１３と、増幅器１３の出力を送信信号Ｓｓとローカル信号Ｌとに電力分配する分配器１５と、送信信号Ｓｓに応じたレーダ波を発射する送信アンテナ１７と、を備える。尚、図２（ａ）は、レーダ波の送信信号Ｓｓ及び受信信号Ｓｒを示したグラフ（上段）及びビート信号ＢＴを示したグラフ（下段）である。図２（ａ）に示すように、本実施例のレーダ装置１は、レーダ波を、三角形状に周波数変調して発射する。
【００３４】
更に、レーダ装置１は、前方物体にて反射されたレーダ波（反射波）を受信するＫ個のアンテナ素子が一列に配置されたリニアアレーアンテナとしての受信アンテナ１９と、受信アンテナ１９を構成するアンテナ素子の一つを順次選択し、選択したアンテナ素子からの受信信号Ｓｒを後段に供給する受信スイッチ２１と、受信スイッチ２１から供給される受信信号Ｓｒを増幅する増幅器２３と、増幅器２３にて増幅された受信信号Ｓｒ及びローカル信号Ｌを混合して、ビート信号ＢＴを生成（図２（ａ）下段参照）するミキサ２５と、ミキサ２５が生成したビート信号ＢＴから不要な信号成分を除去するフィルタ２７と、フィルタ２７の出力をサンプリングし、ディジタルデータに変換するＡ／Ｄ変換器２９とを備える。以下では、Ｋ個のアンテナ素子の夫々を番号付けして、第ｋアンテナ素子（ｋ＝１，２，…，Ｋ）と表現する。
【００３５】
この他、レーダ装置１は、信号処理や装置内各部の制御を行うための構成として、マイクロコンピュータを中心に構成される信号処理部３０と、確率分布テーブル（詳細後述）を記憶する記憶部４０とを備える。
【００３６】
信号処理部３０は、発振器１１の起動／停止等を制御すると共に、マイクロコンピュータでのプログラム実行により、Ａ／Ｄ変換器２９から入力されるビート信号ＢＴのサンプリングデータ（ディジタルデータ）を用いた信号処理や、当該信号処理により得られる前方物体の位置・相対速度・方位・種類等の情報を、図示しない通信インタフェースを介して車内の電子制御装置（ＥＣＵ）に送信する処理を行う。例えば、レーダ装置１は、車間制御を行う電子制御装置である車間制御ＥＣＵと通信ケーブルを通じてシリアル通信可能に接続されて、車間制御ＥＣＵに対して上記情報を提供する。この他、レーダ装置１は、車内ＬＡＮに接続されて、同一の車内ＬＡＮに接続された電子制御装置に上記情報を提供する構成にされてもよい。
【００３７】
このレーダ装置１では、信号処理部３０からの指令に従って発振器１１が起動し、当該起動により発振器１１が生成した高周波信号は、増幅器１３にて増幅された後、分配器１５に入力され、分配器１５によって電力分配される。これにより、レーダ装置１では、送信信号Ｓｓ及びローカル信号Ｌが生成され、送信信号Ｓｓは、送信アンテナ１７を介し、レーダ波として発射される。
【００３８】
一方、送信アンテナ１７から発射され前方物体に反射して戻ってきたレーダ波（反射波）は、受信アンテナ１９を構成する各アンテナ素子にて受信され、各アンテナ素子からは、受信スイッチ２１に向けて、その受信信号Ｓｒが出力される。
【００３９】
また、受信スイッチ２１からは、受信スイッチ２１によって選択された第ｋアンテナ素子（ｋ＝１，…，Ｋ）の受信信号Ｓｒのみが増幅器２３に出力され、増幅器２３で増幅された受信信号は、ミキサ２５に供給される。
【００４０】
ミキサ２５では、受信信号Ｓｒに分配器１５からのローカル信号Ｌが混合されて、信号Ｓｒ，Ｓｓの差の周波数成分であるビート信号ＢＴが生成される。このビート信号ＢＴは、フィルタ２７にて不要な信号成分が除去された後、Ａ／Ｄ変換器２９を通じて、ディジタルデータとして信号処理部３０に取り込まれる。
【００４１】
但し、受信スイッチ２１は、レーダ波の変調周期Ｔｍの間に、全てのアンテナ素子ＡＮ＿１〜ＡＮ＿Ｋを所定回ずつ選択するように、切り替えられる。そして、Ａ／Ｄ変換器２９は、この切替タイミングに同期してデータサンプリングを行う。
【００４２】
信号処理部３０は、このようにＡ／Ｄ変換器２９にてサンプリングされるビート信号ＢＴの当該サンプリングデータ（ディジタルデータ）を、アンテナ素子毎のサンプリングデータに分離し、当該アンテナ素子毎のビート信号ＢＴのサンプリングデータを解析することにより、周知の手法で、反射波の発生元である前方物体までの距離Ｒ、自車両に対する前方物体の相対速度Ｖ、及び、自車両の進行方向を基準とした前方物体の方位θを推定し、更には、後述する本実施例に特有な手法で、前方物体の種類を判別する。
【００４３】
送信信号Ｓｓに基づくレーダ波を送信アンテナ１７が送信したことに起因して、受信アンテナ１９が、反射波を受信すると、受信信号Ｓｒは、図２上段に点線で示すように、レーダ波が前方物体との間を往復するのに要した時間、即ち、前方物体までの距離Ｒに応じた時間Ｔｒ分遅延し、前方物体との相対速度Ｖに応じた周波数ｆｄ分ドップラシフトする。そして、この時間Ｔｒ及び周波数ｆｄの情報は、図２（ａ）下段からも理解できるように、送信信号Ｓｓの周波数が上昇する上り区間及び送信信号Ｓｓの周波数が下降する下り区間の各区間におけるビート信号ＢＴの周波数の違いとして表れる。
【００４４】
このため、信号処理部３０は、１変調周期Ｔｍにおける上り区間及び下り区間（図２（ａ）上段参照）の夫々において、上記サンプリングデータ（ディジタルデータ）として得られるアンテナ素子毎のビート信号ＢＴを周波数変換（フーリエ変換）し、区間毎に、アンテナ素子毎のパワースペクトルを求める（図２（ｂ）参照）。
【００４５】
そして、各アンテナ素子のパワースペクトルに共通するピーク周波数ｆｐを、区間毎に検出し、上り区間のピーク周波数ｆｐ＝ｆｂ１と、下り区間のピーク周波数ｆｐ＝ｆｂ２とから、前方物体との距離Ｒ及び相対速度Ｖを、式（１）〜（４）に従って算出する。但し、ｃは電波伝搬速度，ｆｍは送信信号Ｓｓの変調周波数、Δｆは送信信号Ｓｓの周波数変動幅、ｆ０は送信信号Ｓｓの中心周波数である。
【００４６】
【数１】
【００４７】
一方、受信アンテナ１９の各アンテナ素子が受信する反射波には、到来方向に応じた位相差が生じる。信号処理部３０は、ビート信号ＢＴに含まれるこのような位相差の情報に基づいて、前方物体の方位θを推定する。方位推定の手法としては、ディジタルビームフォーミング（ＤＢＦ）法、ＣＡＰＯＮ法、ＥＳＰＲＩＴ法やＭＵＳＩＣ法など種々の方法が知られているが、例えば、ＭＵＳＩＣ法によって方位推定を行う場合、信号処理部３０は、上記検出されたピーク周波数ｆｐ毎に、各アンテナ素子のビート信号ＢＴの周波数変換値（フーリエ変換値）であって、当該ピーク周波数ｆｐでの周波数変換値を用いて自己相関行列を生成し、この自己相関行列の固有値・固有ベクトルに基づいてＭＵＳＩＣスペクトルを算出する。
【００４８】
このようにＭＵＳＩＣスペクトルを算出すると、例えば、自車両前方に先行車両が存在する場合には、図３に示すように、先行車両に対応する方位に鋭いピークが立つ。本実施例では、このようなＭＵＳＩＣスペクトルのピークを抽出し、各ピーク周波数ｆｐに対応するレーダ波の反射波成分の到来方位θ1，…，θMを推定する。
【００４９】
尚、上述したように、距離Ｒ及び相対速度Ｖの算出の際には、上り区間のピーク周波数ｆｂ１と下り区間のピーク周波数ｆｂ２とのペアを特定する必要があるが、前方物体が複数ある場合には、区間毎に複数のピーク周波数ｆｐが検出されることになる。このような区間毎に複数のピーク周波数ｆｐが検出される環境での上り区間のピーク周波数ｆｂ１と下り区間のピーク周波数ｆｂ２とのペアマッチングは、例えば、ピーク周波数ｆｐ毎の方位θ1，…，θMの推定結果を指標に、反射波の到来方位が対応するピーク周波数ｆｂ１，ｆｂ２の組をマッチングする。このような処理により、信号処理部３０は、前方物体の距離Ｒ及び相対速度Ｖ及び方位θを推定する。
【００５０】
この他、信号処理部３０は、各方位θ1，…，θMから到来する反射波のパワー（受信電力）Ｐ1，…，ＰMを所定の演算式に従って算出し、各パワーＰ1，…，ＰMの１変調周期Ｔｍ前のパワーＰｂ1，…，ＰｂMからの変化量Ｙを算出することで、当該変化量Ｙを指標に、対応する方位θから到来するレーダ波を反射した前方物体の種類（以下、「反射元種類」とも表現する。）を判別する。
【００５１】
反射波の位相は、反射面の位置に応じて変化するため、反射面が揺らげば（振動すれば）、その揺らぎによる反射面の位置変化が反射波の位相変化に現れることになる。一方、レーダ波を反射する前方物体は、微小な複数の反射面の集合にて構成されるものと解釈することができ、その反射面の揺らぎ方は、前方物体が金属物のような表面が硬いものなのか、人のような表面が柔らかいものなのか、各部位（例えば腕や脚）が変位するものなのか等によって異なる。そして、揺らぎ方が異なれば、変化量Ｙの確率分布に差異が生じる。本実施例では、このような現象を利用し、上記変化量Ｙに基づいて、反射元の種類を判別する。
【００５２】
ここで、一つの物体からの反射波について考えてみると、時刻ｔでの反射波は、次式により表現することができる。
【００５３】
【数２】
【００５４】
但し、上式は、物体の反射面が複数の独立に動く等ＲＣＳの反射面から構成されることを前提としたものであり、Ｃは、レーダ波の送信ゲイン、受信ゲイン、送信電力及びＲＣＳ値によって定まる値である。また、ここでは、時刻ｔにおけるｉ番目の反射面までの距離Ｒi（ｔ）を、Ｒi（ｔ）＝Ｒ（ｔ）＋ΔＲi（ｔ）で表現する。Ｒ（ｔ）は、複数反射面までの距離の平均値を表し、ΔＲi（ｔ）は、平均値Ｒ（ｔ）からの距離差分を表す。また、ｆは、レーダ波の送信周波数を表し、λは、レーダ波の送信波長を表す。尚、近似式は、ΔＲi（ｔ）が微小であることを利用して、分母のＲ（ｔ）＋ΔＲi（ｔ）を、Ｒ（ｔ）に近似したものである。
【００５５】
また、時刻ｔから、微小な観測周期Ｔｓｗ分進んだ時点での反射波は、次式により表現することができる。
【００５６】
【数３】
【００５７】
但し、上式で用いるｖ（ｔ）は、レーダ波を反射した前方物体の相対速度を表し、Δｖi（ｔ）は、ｉ番目の反射面での相対速度Ｖ（ｔ）に対する速度差分、即ち、揺らぎ量を表す。近似式では、時刻ｔ＋Ｔｓｗでの周波数ｆTswを、時刻ｔでの周波数ｆと同一値に近似している。
【００５８】
時刻ｔ及び時刻ｔ＋Ｔｓｗでの反射波は、上記のように近似できることから、例えば、変化量Ｙとして、レーダ波の受信信号のパワー比を算出すると、変化量Ｙは次式で表現される。
【００５９】
【数４】
【００６０】
上式では、パワー比を対数増幅したものを変化量Ｙとして表す。変化量Ｙとして、時刻ｔ＋Ｔｓｗでのレーダ波の受信信号のパワーと、時刻ｔでのレーダ波の受信信号のパワーとの比を算出すれば、変化量Ｙは、ΔＲi（ｔ）及びΔＶi（ｔ）によって変化する値として算出される。このため、変化量Ｙを、予め実験で得た反射元種類毎の変化量Ｙの確率分布と比較することで、反射元（前方物体）の種類を判別することができるのである。
【００６１】
但し、観測周期Ｔｓｗの間に生じる各反射面からの反射波成分の位相変化量Δφがπ以上となると、各反射面の反射波成分が混合する受信アンテナ１９の反射波の受信信号では、その振幅がレイリー分布に従うようになり、反射面の揺らぎを観測することが難しい。
【００６２】
そこで、本実施例では、レーダ波の送信波長λ及び観測周期Ｔｓｗを、判別対象とする各種類の前方物体のいずれにおいても、その前方物体での各反射面から反射する反射波の位相変化量Δφがπを超えないように、送信波長λ及び観測周期Ｔｓｗを定める。
【００６３】
ここで、送信波長λ及び観測周期Ｔｓｗの設定条件について詳述する。時刻ｔにおけるｉ番目の反射面からの反射波の位相をφi（ｔ）で表現し、時刻ｔ＋Ｔｓｗにおけるｉ番目の反射面からの反射波の位相をφi（ｔ＋Ｔｓｗ）で表現すると、観測周期Ｔｓｗの間でのｉ番目の反射面からの反射波の位相変化量Δφiは、Δφi＝φi（ｔ＋Ｔｓｗ）−φi（ｔ）で表現される。そして、観測周期Ｔｓｗでの各反射面からの反射波の位相変化量Δφiが独立に−πからπまでの一様分布に従うとき、これら反射波の混合波の振幅分布は、レイリー分布に従う。
【００６４】
換言すれば、各反射面からの反射波の位相Δφiがπを超えないように送信波長λ及び観測周期Ｔｓｗを定めれば、前方物体の揺らぎを、受信アンテナ１９による反射波の受信信号にて観測できることになる。即ち、送信波長λ及び観測周期Ｔｓｗについては、Δφiの最大値ｍａｘ｛Δφi｝がπ未満となるように定めればよい。
【００６５】
【数５】
【００６６】
一方、Δφiは、次式によって表すことができる。
【００６７】
【数６】
【００６８】
各反射面の揺らぎ量ΔＶi（ｔ）は、平均ゼロ、分散σv2の正規分布Ｎ（０，σv2）に従うと近似的にみなすことができ、分散σv2で揺らぐ前方物体の揺らぎを、３σvの範囲（９９．７％の確率で起こる範囲）で観測できるように、送信波長λ及び観測周期Ｔｓｗを定める場合には、次の関係式を満足すればよい。
【００６９】
【数７】
【００７０】
換言すると、送信波長λ及び観測周期Ｔｓｗは、次の関係式を満足すればよい。
【００７１】
【数８】
【００７２】
本実施例では、このような関係式に従って、判別対象の種類の前方物体の内、揺らぎが最も大きい種類の前方物体の分散σv2を基準に、送信波長λ及び観測周期Ｔｓｗを定める。尚、揺らぎの大きさ（分散σv2）については、実験により計測すればよい。
【００７３】
本実施例では、このような関係式に従って、送信波長λ及び観測周期Ｔｓｗを定めつつ、記憶部４０に、判別対象の反射元種類毎の変化量Ｙの確率分布を表すテーブルを用意する。これによって、変化量Ｙから反射元（前方物体）の種類を判別することができるようにする。
【００７４】
図４には、記憶部４０が記憶する確率分布テーブルの構成を示す。図４に示すように、本実施例の確率分布テーブルは、判別対象の反射元種類毎に、レーダ波の反射元が当該種類の反射元であったときに得られる変化量Ｙの確率分布が記述されてなるテーブルである。具体的に、この確率分布テーブルには、離散的な変化量Ｙの各値に対して、変化量Ｙがその値以下となる確率（累積確率）が、判別対象の反射元種類毎に記述されている。
【００７５】
このテーブルについては、例えば、判別対象の前方物体の種類毎に、当該種類の前方物体に対してレーダ波を発射して反射波の計測データを取得し、この反射波の計測データから得られた変化量Ｙの確率分布を、テーブルに記述することで作成することができる。また、別手法としては、反射波の計測データを取得する一方、シミュレーションにより、分散σv2を変化させて各分散σv2での変化量Ｙの確率分布を求め、反射波の計測データと最も近い分散σv2での変化量Ｙの確率分布を、テーブルに記述することで作成することができる。
【００７６】
ちなみに、本実施例では、確率分布テーブルに、判別対象の反射元種類毎の累積確率分布を記述したが、前方物体の種類を判別する際には、受信したレーダ波から得られた変化量Ｙの累積確率と、隣接する変化量Ｙに対して記述された累積確率との差分を採ることで、累積確率分布を微分し、変化量Ｙに対応する確率密度を求める。このため、確率分布テーブルには、累積確率分布に代えて確率密度分布を記述してもよい。この他、記憶部４０には、確率分布テーブルに代えて、累積分布関数及び確率密度関数の情報を記憶させることも可能である。
【００７７】
また、本実施例では、前方物体として、車両やガードレール、人等を判別して、車両制御に役立てることを考えている。このため、記憶部４０が記憶する確率分布テーブルには、反射元種類として、金属物及び人の夫々についての確率分布が記述されている。特に、本実施例の確率分布テーブルには、「静止している人」及び「移動している人」の夫々についての確率分布が記述されており、レーダ装置１は、前方物体として、「金属物」、「静止している人」、「移動している人」の夫々を判別可能に構成されている。このようなレーダ装置１によれば、前方に移動する人がいると判定された場合には、車両速度を緩める等の車両制御を行うことができて、危険回避等に大変役立つ。この他、本実施例のレーダ装置１によれば、ガードレールと人を見分けることができるので、例えば、衝突が回避不可能となったときに、人を回避してガードレールに衝突するといったような従来では不可能な車両制御を実現することが可能である。
【００７８】
続いては、以上に概略的に説明した信号処理部３０の処理動作を、図５及び図６を用いて詳細に説明する。信号処理部３０は、図５に示すメインルーチンを変調周期Ｔｍ毎に繰返し実行することにより、変調周期Ｔｍ毎に、一連の信号処理（Ｓ１１０〜Ｓ１５０）を実行して、前方物体までの距離Ｒ、前方物体との相対速度Ｖ及び前方物体の方位θを推定すると共に、前方物体の種類を判別する。
【００７９】
具体的には、まず１変調周期Ｔｍにおける上り区間を処理対象区間として、この処理対象区間（上り区間）に受信した反射波よって生成されたビート信号ＢＴのサンプリングデータに基づき、図６に示す物標推定処理を実行する（Ｓ１１０）。その後、処理対象区間を下り区間に変更し、この処理対象区間（下り区間）に受信した反射波によって生成されたビート信号ＢＴのサンプリングデータに基づき、図６に示す物標推定処理を実行する（Ｓ１２０）。
【００８０】
物標推定処理では、処理対象区間におけるアンテナ素子毎のビート信号ＢＴのサンプリングデータを取得し（Ｓ２１０）、その後、これらサンプリングデータに基づき、各アンテナ素子のビート信号ＢＴを周波数変換して、アンテナ素子毎のパワースペクトルを算出する（Ｓ２２０）。その後、各アンテナ素子共通のピーク周波数ｆｐを検出する。ここでは、アンテナ素子毎のパワースペクトルに共通するピークであって予め定められた閾値以上のピークを検出し、検出した閾値以上の各ピークの周波数を、ピーク周波数ｆｐとして検出する（Ｓ２３０）。
【００８１】
ピーク周波数ｆｐの検出後には、検出したピーク周波数ｆｐの一つを処理対象周波数に選択し（Ｓ２４０）、各アンテナ素子における処理対象周波数の周波数変換値（フーリエ変換値）ｘkを配列してなる受信ベクトルＸ＝（ｘ1，ｘ2，…ｘK）Tを生成して、この受信ベクトルＸについての自己相関行列Ｒｘｘ１＝Ｅ［ＸＸH］及びＲｘｘ２＝ＸＸHを算出する（Ｓ２５０）。ここで、Ｔは、ベクトル転置を示し、Ｈは、複素共役転置を示す。また、ｘkは、ｋ番目のアンテナ素子の受信信号から生成された処理対象区間のビート信号ＢＴの周波数変換値であって、処理対象周波数での周波数変換値を示す。従来知られる手法では、方位推定のために自己相関行列Ｒｘｘとして、ＸＸHの期待値Ｅ［ＸＸH］を算出するが、本実施例では、上述の変化量Ｙを算出する際に、各方位から到来する反射波の瞬時パワーを計算する必要があるので、自己相関行列Ｒｘｘ１として、ＸＸHの期待値Ｅ［ＸＸH］を算出する一方、パワー推定のために自己相関行列Ｒｘｘ２として、最新の受信ベクトルＸを用いたＸＸHの瞬時値を算出する。期待値Ｅ［ＸＸH］としては、ＸＸHの最新値を含む過去所定個のＸＸHの平均値を算出することができる。
【００８２】
その後、自己相関行列Ｒｘｘ１についての固有値λ1，…，λK（但し、λ1≧λ2≧…λk）及び固有ベクトルｅ1，…，ｅKを求め、熱雑音電力に対応する閾値λthより大きい固有値の数から到来波数Ｍを推定する（Ｓ２６０）。
【００８３】
そして、熱雑音電力以下となる（Ｋ−Ｍ）個の固有値λM+1，…，λKに対応した固有ベクトルｅM+1，…，ｅKからなる雑音固有ベクトルＥNと、方位θに対するアレーアンテナの複素応答、即ち、ステアリングベクトルａ（θ）とから、角度スペクトルとして、下式の評価関数ＰMU（θ）で表されるＭＵＳＩＣスペクトルを求める。
【００８４】
【数９】
【００８５】
そして、この評価関数ＰMU（θ）で表されるＭＵＳＩＣスペクトルから、到来波数Ｍ分のピークを抽出し、各ピークに対応する方位θ1，…，θMを、反射波の到来方位θ1，…，θM、換言すれば、レーダ波の反射元である前方物体（物標）の方位として推定する（Ｓ２７０）。
【００８６】
また、推定した各方位θ1，…，θMから到来する反射波のパワー（受信電力）を次のように推定して、これを一時記憶する（Ｓ２８０）。
具体的には、まずパワー（受信電力）を推定する対象であるＳ２７０で推定した反射波の各到来方位θ1，…，θMに対応するステアリングベクトルａ（θ1），…，ａ（θM）を用いて、次の方向行列Ａ（所謂アレー応答行列）を生成する。
【００８７】
【数１０】
【００８８】
そして、この方向行列Ａ及び瞬時値としての自己相関行列Ｒｘｘ２を用いて、次式で表される行列Ｓの対角成分を算出する。
【００８９】
【数１１】
【００９０】
行列Ｓの第ｍ対角成分は、方位θmの受信電力Ｐmに対応する（但し、ｍ＝１，…，Ｍ）。即ち、ここでは、行列Ｓの対角成分を算出することにより、Ｓ２７０で推定した各到来方位θ1，…，θMから到来した反射波のパワー（受信電力）を推定する。但し、Ｓ２８０では、Ｓ２７０で反射波の到来方位として推定したものの、上式により算出されたパワーが閾値以下の方位を、反射波の到来方位から外す処理を併せて行う。
【００９１】
また、この処理を終えると、信号処理部３０は、反射波の到来方位θ1，…，θMの一つを、処理対象方位に選択し（Ｓ２９０）、選択した処理対象方位から到来した反射波のパワーの前回推定値が存在するか否かを判断する（Ｓ３００）。
【００９２】
即ち、ここでは、１変調周期Ｔｍ前の処理対象区間における反射波の受信信号に基づき、今回の処理対象周波数と同一のピーク周波数ｆｐでの方位推定及びパワー推定が行われて、処理対象方位と同一方位から到来した反射波のパワー推定がなされているか否かを判断し、なされている場合には、パワーの前回推定値が存在すると判断し、なされていない場合には、パワーの前回推定値が存在しないと判断する。そして、パワーの前回推定値が存在しないと判断すると（Ｓ３００でＮｏ）、Ｓ３４０に移行し、パワーの前回推定値が存在すると判断すると（Ｓ３００でＹｅｓ）、Ｓ３１０に移行する。
【００９３】
Ｓ３１０に移行すると、信号処理部３０は、処理対象方位について今回推定したパワーの前回推定値からの変化量Ｙを算出する。ここで、今回推定したパワーの推定値をＰと表現し、前回推定したパワーの推定値をＰｂと表現すると、変化量Ｙについては、次式に従って算出することになる。尚、ここで対数増幅を行っているのは、パワーの推定値をデシベル単位で記憶保持しておくことで、差分により変化量Ｙを簡単に算出することができるためである。
【００９４】
【数１２】
【００９５】
ちなみに、本実施例では、このように、前回のパワー推定値として、１変調周期Ｔｍ前の値を用いるので、上述した観測周期Ｔｓｗは、変調周期Ｔｍに一致する。
また、この処理を終えると、信号処理部３０は、Ｓ３２０に移行し、判別対象とする反射元種類毎に、Ｓ３１０で算出された変化量Ｙが仮に当該種類の反射元から到来したレーダ波に基づき算出された変化量Ｙである場合に当該変化量Ｙが得られる確率Ｐｒ（以下、変化量Ｙの発生確率Ｐｒとも言う。）を、記憶部４０が記憶する確率分布テーブルに記された該当種類の確率分布の情報に基づいて導出する（Ｓ３２０）。尚、上述したように確率分布テーブルには、確率分布として各変化量Ｙでの累積確率が離散的に記されていることから、ここでは、Ｓ３１０で算出された変化量Ｙに対応する累積確率と隣接する変化量Ｙの累積確率との差分を算出することで、変化量Ｙに対応する累積確率の微分である確率密度を算出する。即ち、Ｓ３２０では、上記手法で、Ｓ３１０で算出された変化量Ｙの確率密度を算出することで、変化量Ｙの発生確率Ｐｒを反射元種類毎に算出する。
【００９６】
このようにして、判別対象とする反射元種類毎の発生確率Ｐｒを算出すると、信号処理部３０は、上記種類毎の発生確率Ｐｒの中から、最も値の高いものを抽出し、この最も発生確率Ｐｒの高い反射元の種類を、処理対象方位から到来した反射波の反射元の種類であると判別する（Ｓ３３０）。
【００９７】
その後、Ｓ２７０で反射波の到来方位として推定した全方位θ1，…，θM（但しＳ２８０で外されたものを除く。）に関して、当該方位から到来するレーダ波の反射元の種類を判別したか否かを判断し（Ｓ３４０）、判別していない場合には（Ｓ３４０でＮｏ）、Ｓ２９０に移行し、処理対象方位を切り替えて、Ｓ３００以降の処理を実行する。
【００９８】
一方、全方位に関して上記判別を行った場合には（Ｓ３４０でＹｅｓ）、全ピーク周波数ｆｐについて、これを処理対象周波数として選択し、Ｓ２５０以降の処理を実行したか否かを判断し、全ピーク周波数ｆｐに対して実行していなければ（Ｓ３５０でＮｏ）、Ｓ２４０に移行し、処理対象周波数を切り替えて、Ｓ２５０以降の処理を実行する。そして、全ピーク周波数ｆｐに対して実行すると、当該物標推定処理を終了する。
【００９９】
このようにして上り区間及び下り区間の物標推定処理を終了すると、信号処理部３０は、Ｓ１３０（図５参照）に移行し、上り区間でのピークと、下り区間でのピークとのペアマッチングを行い、上り区間でのピーク周波数ｆｐ＝ｆｂ１と下り区間でのピーク周波数ｆｐ＝ｆｂ２との組合せを特定する。そして特定した組合せに基づいて、上述した手法で前方物体までの距離Ｒ及び前方物体との相対速度Ｖを推定する（Ｓ１４０）。
【０１００】
そして、各ピーク周波数ｆｐに対して推定された反射波の到来方位θ1，…，θM及び各到来方位θ1，…，θMから到来するレーダ波の反射元種類の情報を、上記推定された距離Ｒ及び相対速度Ｖの情報と組み合わせて、自車両周囲に位置する各前方物体の距離Ｒ、相対速度Ｖ、方位θ及び種類の情報を、通信インタフェースを通じて外部の電子制御装置に出力する（Ｓ１５０）。信号処理部３０は、このような処理を変調周期Ｔｍ毎に繰返し実行する。
【０１０１】
以上、本実施例のレーダ装置１の構成について説明したが、このレーダ装置１によれば、従来のレーダ装置１とは異なり、前方物体の距離Ｒ、相対速度Ｖ及び方位θに加えて、前方物体の種類を推定することが可能である。特に本実施例によれば、前方物体の種類として、「金属物（車両やガードレール等）」、「静止している人」及び「移動している人」の夫々を判別することができる。
【０１０２】
従って、本実施例によれば、前方物体が金属物か人であるかによって車両制御のパターンを切り替えることができ、更には人が移動中であるか静止しているかによっても車両制御のパターンを切り替えることができて、上述したように危険回避等に大変役立つ。
【０１０３】
さて、上記実施例では、確率分布テーブルが示す確率分布とレーダ波の受信信号に基づく変化量Ｙとの比較により、反射元の種類を判別するようにしたが、各種類の反射元が存在する可能性については、道路環境によっても事前にある程度予測することができる。例えば、街中では、人が存在する可能性が高く、郊外では、街中よりも人が存在する可能性が低く、高速道路においては、車外に人が存在することは通常ありえない。
【０１０４】
そこで、変化量Ｙの発生確率Ｐｒと道路環境とから、変化量Ｙに対応するレーダ波の反射元が判別対象の各種類の反射元である確度を評価し、この確度に基づいて、反射元の種類を判別するように、レーダ装置を構成してもよい（第一変形例及び第二変形例）。
［第一変形例］
続いて、第一変形例について説明する。第一変形例のレーダ装置２は、記憶部４０に対応する記憶部４１が確率分布テーブルに加えて補正テーブルを記憶する点、信号処理部３０に対応する信号処理部３１が車内ＬＡＮを通じてナビゲーションＥＣＵ１００と通信可能に構成されている点、及び、信号処理部３１が実行する物標推定処理の一部が上記実施例と異なる点を除けば、基本的に上記実施例のレーダ装置１と同一構成にされている。従って、以下では、第一変形例のレーダ装置２の説明として、上記実施例のレーダ装置１と同一構成の説明を省略する。
【０１０５】
図７に示すように第一変形例のレーダ装置２が備える信号処理部３１は、車内ＬＡＮを通じてナビゲーションＥＣＵ１００と通信可能に接続され、ナビゲーションＥＣＵ１００から、自車両の現在位置に対応する道路環境のカテゴリ情報を取得可能な構成にされている。
【０１０６】
ナビゲーションＥＣＵ１００は、ＧＰＳ受信機等から構成される自車両の現在位置を検出可能な位置検出器１０１、及び、道路地図データベース１０３に接続されており、図７に示す道路環境情報送信処理を繰返し実行して、位置検出器１０１から自車両の現在位置情報を取得し（Ｓ４１０）、この取得情報が示す自車両の現在位置に対応する道路環境のカテゴリを、道路地図データベース１０３が記憶する情報に基づいて判定する（Ｓ４２０）。そして、判定した自車両の現在位置に対応する道路環境のカテゴリ情報をレーダ装置２に車内ＬＡＮを通じて送信する（Ｓ４３０）。
【０１０７】
尚、道路地図データベース１０３には、地図収録エリア内の各地域毎に、当該地域が「街中」及び「郊外」のいずれであるかを表すエリア情報が記憶されている。更に、高速道路に該当する道路に対しては「高速道路」であることを示すエリア情報が記憶されている。即ち、Ｓ４２０では、自車両の現在位置に対応する走行道路を特定して、走行道路が「高速道路」であるか否かを判断し、走行道路が「高速道路」である場合には、自車両の現在位置に対応する道路環境のカテゴリを「高速道路」であると判定し、走行道路が「高速道路」でない場合には、道路地図データベース１０３が記憶する自車両の現在位置に対応する地域のエリア情報が示すカテゴリ（「街中」又は「郊外」）を、自車両の現在位置に対応する道路環境のカテゴリであると判定する。
【０１０８】
この動作によって、ナビゲーションＥＣＵ１００からは、自車両の現在位置に対応する道路環境のカテゴリとして「街中」「郊外」「高速道路」のいずれかのカテゴリを示す情報がレーダ装置２の信号処理部３１に提供される。
【０１０９】
レーダ装置２の信号処理部３１は、このような道路環境のカテゴリ情報をナビゲーションＥＣＵ１００から逐次取得し、物標推定処理のＳ３２０（図６参照）の処理を終えると、Ｓ３３０の処理に代えて、Ｓ５１０〜Ｓ５４０（図８参照）の処理を実行する。
【０１１０】
具体的に、信号処理部３０は、ナビゲーションＥＣＵ１００から得られた最新の上記カテゴリ情報に基づき、現在位置に対応する道路環境のカテゴリを特定する（Ｓ５１０）。そして、記憶部４１が記憶する補正テーブルから、上記特定したカテゴリに対応する補正係数Ｗの組であって反射元種類毎の補正係数Ｗを読み出し（Ｓ５２０）、Ｓ３２０の処理で反射元種類毎に導出した変化量Ｙの発生確率（確率密度）Ｐｒの夫々に対し、対応する反射元種類の補正係数Ｗを乗算する。これによって、反射元種類毎に、変化量Ｙを算出する元となったレーダ波を反射してきた前方物体（反射元）が当該種類の反射元である確度Ｚ＝Ｗ・Ｐｒを算出する（Ｓ５３０）。
【０１１１】
尚、補正テーブルは、図７に示すように、「街中」「郊外」「高速道路」の各道路環境のカテゴリ毎に、補正係数Ｗの組として、反射元種類毎の補正係数Ｗが記述されてなるテーブルである。例えば、「街中」の補正係数Ｗの組としては、「静止している人」に対応する補正係数Ｗ１１＝１／３、「移動している人」に対応する補正係数Ｗ１２＝１／３、「金属物」に対応する補正係数Ｗ１３＝１／３が補正テーブルに記述され、「郊外」の補正係数の組としては、「静止している人」に対応する補正係数Ｗ２１＝０．２５、「移動している人」に対応する補正係数Ｗ２２＝０．２５、「金属物」に対応する補正係数Ｗ２３＝０．５が補正テーブルに記述される。この他、「高速道路」の補正係数の組としては、例えば「静止している人」に対応する補正係数Ｗ３１＝０、「移動している人」に対応する補正係数Ｗ３２＝０、「金属物」に対応する補正係数Ｗ３３＝１が補正テーブルに記述される。
【０１１２】
自車両の現在位置が「郊外」に該当する場合には、「静止している人」の確率分布に基づき導出した変化量Ｙの発生確率ＰｒにＷ２１＝０．２５を乗算して、反射元が「静止している人」である確度Ｚ＝０．２５・Ｐｒを算出するといった具合である。同様に、「移動している人」の確率分布に基づき導出した変化量Ｙの発生確率Ｐｒに、Ｗ２２＝０．２５を乗算して、反射元が「移動している人」である確度Ｚ＝０．２５・Ｐｒを算出し、「金属物」の確率分布に基づき導出した変化量Ｙの発生確率Ｐｒに、Ｗ２３＝０．５を乗算して、反射元が「金属物」である確度Ｚ＝０．５・Ｐｒを算出する。
【０１１３】
そして、算出した確度Ｚが最も大きい反射元種類を、上記処理対象方位から到来した反射波の反射元種類であると判別する（Ｓ５４０）。その後、信号処理部３１は、上記実施例で説明したＳ３４０以降の処理を実行する。
【０１１４】
以上、第一変形例のレーダ装置２について説明したが、この変形例によれば、道路環境を加味して、反射元の種類を判別するようにレーダ装置２を構成したので、道路環境による事前予測に基づいて、高精度に反射元の種類を判別することができる。
【０１１５】
また、第一変形例では、道路地図データベース１０３と位置検出器１０１とを用いて、自車両の道路環境を判定するようにしたが、自車両の走行速度に基づき、自車両の道路環境を判定するように、レーダ装置は、構成されてもよい。
［第二変形例］
続いて、第二変形例のレーダ装置３について説明する。但し、第二変形例のレーダ装置３は、信号処理部３１に対応する信号処理部３２がナビゲーションＥＣＵ１００を通じて上記道路環境のカテゴリ情報を取得する代わりに速度センサ１０５から自車両の速度情報を取得する構成にされている点、及び、信号処理部３２が実行する物標推定処理の一部が上記第一変形例と異なる点を除けば、基本的に第一変形例のレーダ装置２と同一構成にされている。従って、以下では、第二変形例のレーダ装置３の説明として、上述のレーダ装置１，２と同一構成の説明を省略する。
【０１１６】
第二変形例のレーダ装置３の信号処理部３２は、上述したように、道路環境のカテゴリ情報をナビゲーションＥＣＵ１００から取得する代わりに、速度センサ１０５から逐次自車両の速度情報を取得する。一方で、物標推定処理のＳ３２０の処理を終えると、Ｓ５１０〜Ｓ５４０の処理に代えて、Ｓ６１０〜Ｓ６４０（図９参照）の処理を実行する。
【０１１７】
具体的に、信号処理部３０は、Ｓ６１０に移行すると、速度センサ１０５から得られた速度情報に基づき、自車両周辺の道路環境を推定する（Ｓ６１０）。具体的には、現在から過去の所定時点までの期間（例えば現在から過去５分間）における車両の平均速度を算出し、この平均速度を指標に、自車両周囲の道路環境のカテゴリを判定する。
【０１１８】
平均速度が第一の境界値（例えば３０ｋｍ）未満であれば、自車両周辺の道路環境のカテゴリを「街中」と判定し、平均速度が第一の境界値（例えば３０ｋｍ）以上第二の境界値(例えば６０ｋｍ)未満であれば、自車両周辺の道路環境のカテゴリを「郊外」と判定し、平均速度が第二の境界値以上である場合には、自車両周辺の道路環境のカテゴリを「高速道路」と判定するといった具合である。
【０１１９】
このようにして道路環境を推定すると、信号処理部３２は、Ｓ５２０，Ｓ５３０での処理と同様に、記憶部４０が記憶する補正テーブルから、自車両周辺の道路環境のカテゴリに対応する反射元種類毎の補正係数Ｗを読み出し（Ｓ６２０）、反射元種類毎に確率分布テーブルから導出した変化量Ｙの発生確率（確率密度）Ｐｒの夫々に対し、上記読み出した反射種類毎の補正係数Ｗの内、対応する反射元種類の補正係数Ｗを乗算する。これによって、反射元種類毎に、変化量Ｙを算出する元となったレーダ波を反射してきた前方物体（反射元）が当該種類の反射元である確度Ｚ＝Ｗ・Ｐｒを算出する（Ｓ６３０）。
【０１２０】
そして、Ｓ５４０の処理と同様、算出した確度Ｚが最も大きい反射元種類を上記処理対象方位から到来した反射波の反射元種類であると判別し（Ｓ６４０）、その後、上記実施例で説明したＳ３４０以降の処理を実行する。
【０１２１】
以上、第二変形例のレーダ装置２について説明したが、第二変形例によれば、道路環境を加味して、反射元の種類を判別するようにレーダ装置３を構成したので、道路環境による事前予測に基づいて、高精度に反射元種類を判別することができる。特に、第二変形例によれば、車両速度に基づき道路環境を推定するので、道路地図データベース１０３や位置検出器１０１が搭載されていない車両においても、道路環境を加味した反射元種類の判別を行うことができ、安価に様々な車両に適用することができる。
【０１２２】
以上、本発明の実施例について説明したが、レーダ装置１，２，３が備える発振器１１、増幅器１３及び送信アンテナ１７、並びに、受信アンテナ１９、受信スイッチ２１及び増幅器２３は、本発明の送受信手段の一例に対応する。また、信号処理部３０，３１，３２が実行するＳ２８０の処理は、本発明の検出手段により実行される処理の一例に対応する。この他、信号処理部３０，３１，３２が実行するＳ３１０の処理は、本発明の変化量算出手段により実行される処理の一例に対応する。
【０１２３】
この他、記憶部４０における確率分布テーブルの記憶領域は、本発明の確率分布記憶手段の一例に対応し、信号処理部３０が実行するＳ３２０，３３０の処理、又は、信号処理部３１が実行するＳ３２０，Ｓ５２０〜Ｓ５４０の処理、又は、信号処理部３２が実行するＳ３２０，Ｓ６２０〜Ｓ６４０の処理は、本発明の判別手段にて実行される処理の一例に対応する。この他、信号処理部３１が実行するＳ５１０の処理又は信号処理部３２が実行するＳ６１０の処理は、本発明の環境推定手段にて実行される処理の一例に対応する。
【０１２４】
また、レーダ装置１，２，３が備えるミキサ２５は、本発明のビート信号生成手段に対応し、信号処理部３０，３１，３２が実行するＳ１１０〜Ｓ１４０及び物標推定処理におけるＳ２１０〜Ｓ２７０の処理は、本発明の信号解析手段により実行される処理に対応する。
【０１２５】
また、本発明は、上記実施例に限定されるものではなく、種々の態様を採ることができる。例えば、本発明は、ＦＭＣＷレーダ装置に限らず、その他の方式にてレーダ波の送受信を行うレーダ装置に適用することができる。具体的には、ＳＦＷ（Stepped Frequency Waveforms）方式にてレーダ波の送受信を行うＳＦＷレーダ装置に適用することができる。この他、本発明は、車載用のレーダ装置に限らず、その他の用途に用いられるレーダ装置にも適用することができる。また、上記実施例では、人や金属物を判別可能なようにレーダ装置１，２，３を構成したが、判別対象とする反射元の種類は、これに限定されない。また、上記実施例では、道路地図データベース１０３又は速度センサ１０５のいずれか一方を用いて道路環境を推定するようにレーダ装置２，３を構成したが、レーダ装置は、これら両者の情報を用いて道路環境を推定するように構成されてもよい。
【符号の説明】
【０１２６】
１，２，３…レーダ装置、１１…発振器、１３…増幅器、１５…分配器、１７…送信アンテナ、１９…受信アンテナ、２１…受信スイッチ、２３…増幅器、２５…ミキサ、２７…フィルタ、２９…Ａ／Ｄ変換器、３０，３１，３２…信号処理部、４０，４１…記憶部、１００…ナビゲーションＥＣＵ、１０１…位置検出器、１０３…道路地図データベース、１０５…速度センサ
【技術分野】
【０００１】
本発明は、レーダ装置に関する。
【背景技術】
【０００２】
従来、レーダ波を発射して、このレーダ波が前方物体に反射して戻ってくる反射波を受信し、送信波に対する受信波の遅延量から前方物体までの距離を検出し、送信波に対する受信波の周波数変化量から前方物体との相対速度を検出するレーダ装置が知られている（例えば特許文献１参照）。この種のレーダ装置は、例えば車両に搭載されて、車両制御に使用される。
【０００３】
この他、レーダ装置としては、物標の存在有無を受信信号の振幅と閾値との比較により判定するものが知られている。この手法では、クラッタにより物標の存在を誤判定する確率（誤警報確率）を一定に抑えるために、例えば、レーダ波の受信信号を対数増幅し、この対数増幅後の受信信号を時間的に前後する受信信号の対数増幅後の平均値で減算する（非特許文献１参照）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１０−０３２３１４号公報
【非特許文献】
【０００５】
【非特許文献１】関根松夫著、「レーダ信号処理技術」、社団法人 電子情報通信学会、１９９１年９月、ｐ．９６−１０３
【発明の概要】
【発明が解決しようとする課題】
【０００６】
ところで、従来知られるレーダ装置では、前方物体までの距離や前方物体の方位、前方物体との相対速度を検出することはできるものの、レーダ波を反射した前方物体の種類までを特定することはできない。従来技術によれば、受信波が物標からの反射により生じたものであるのか、それともクラッタであるのかを判別することができる程度である。
【０００７】
本発明は、こうした問題に鑑みなされたものであり、反射波の受信信号から、反射元の物体の種類を判別可能な技術を提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記目的を達成するためになされた本発明のレーダ装置は、レーダ波を発射し、当該レーダ波の反射波を受信する送受信手段と、送受信手段が受信した反射波の受信信号に基づき、反射波の信号強度を検出する検出手段と、を備えるレーダ装置であって、更に、変化量算出手段、確率分布記憶手段、及び、判別手段、を備えるものである。
【０００９】
このレーダ装置が備える変化量算出手段は、検出手段により検出された反射波の信号強度Ｐと、当該信号強度Ｐに対応する反射波の受信時点から所定の観測周期Ｔｓｗ遡った時点で送受信手段により受信された反射波の検出手段により検出された信号強度Ｐｂと、に基づき、信号強度Ｐｂに対する信号強度Ｐの変化量Ｙを算出する。変化量Ｙとしては、例えば信号強度Ｐ，Ｐｂの比Ｐ／Ｐｂを挙げることができ、信号強度Ｐ，Ｐｂを対数増幅する場合には、変化量Ｙとして、対数増幅後の信号強度ｌｏｇ（Ｐ）及びｌｏｇ（Ｐｂ）の差分ｌｏｇ（Ｐ）−ｌｏｇ（Ｐｂ）を算出することができる。
【００１０】
一方、確率分布記憶手段は、判別対象とする反射元の種類毎に、レーダ波が当該種類の反射元で反射して到来する反射波から算出される変化量Ｙの確率分布を記憶する。例えば、確率分布記憶手段は、上記確率分布を、累積分布関数や確率密度関数等の関数情報、又は、これら関数に代わるテーブル情報の形態で記憶する。
【００１１】
そして、判別手段は、確率分布記憶手段が記憶する上記種類毎の確率分布に基づき、変化量算出手段で算出された変化量Ｙに対応するレーダ波の反射元の種類を判別する。この他、このレーダ装置では、変化量Ｙからレーダ波の反射元の種類を判別可能とするため、判別対象とする各種類の反射元にて発生する揺らぎの特徴が変化量Ｙに反映されるような関係で、レーダ波の送信波長λ及び観測周期Ｔｓｗが定められている。
【００１２】
反射波の位相は、反射面の位置に応じて変化するため、反射面が揺らげば（振動すれば）、その揺らぎによる反射面の位置変化が反射波の位相変化に現れることになる。一方、反射元は、微小な複数の反射面の集合にて構成されるものと解釈することができ、その反射面の揺らぎ方は、反射元が金属物のような表面が硬いものなのか、人のような表面が柔らかいものなのか、各部位（例えば腕や脚）が変位するものなのか等によって異なる。そして、揺らぎ方が異なれば、変化量Ｙの確率分布に差異が生じる。本発明では、このような現象を利用し、上記変化量Ｙに基づいて、反射元の種類を判別する。
【００１３】
但し、観測周期Ｔｓｗの間に生じる各反射面から到来する反射波成分の位相変化量Δφがπ以上となると、各反射面の反射波成分が混合する上記送受信手段での反射波の受信信号においては、その振幅がレイリー分布に従うようになり、反射面の揺らぎを観測することが難しい。換言すれば、判別対象とする各種類の反射元にて発生する揺らぎの特徴が変化量Ｙに反映されなくなり、変化量Ｙに揺らぎ方による差異が表れないため、変化量Ｙに基づいて、反射元の種類を判別することが難しい。
【００１４】
そこで、本発明では、レーダ波の送信波長λ及び観測周期Ｔｓｗを、判別対象とする各種類の反射元にて発生する揺らぎの特徴が変化量Ｙに反映される関係となるように定めることで、揺らぎを、変化量Ｙにより観測することができるようにしている。
【００１５】
尚、観測周期Ｔｓｗ間の位相変化量Δφは、当然のごとく観測周期Ｔｓｗとレーダ波の送信波長λに依存する。更に、位相変化量Δφは、反射元の揺らぎの大きさ（反射面の変位幅）に依存する。よって、送信波長λ及び観測周期Ｔｓｗについては、判別対象の各種類の反射元での揺らぎの大きさを考慮して、判別対象とする各種類の反射元にて発生する揺らぎの特徴が変化量Ｙに反映される関係となるように定める。具体的には、判別対象とする複数種類の反射元のいずれにおいても、その反射元での各反射面から反射する反射波の位相変化量Δφがπを超えないように、送信波長λ及び観測周期Ｔｓｗを定めることで、上記関係を満足するように送信波長λ及び観測周期Ｔｓｗを定めることができる。
【００１６】
このように構成された本発明のレーダ装置によれば、反射波の受信信号から、反射元の物体の種類を判別することができる。従って、本発明のレーダ装置を用いれば、前方物体の種類を判別し、この判別結果に基づき、従来では実現することのできなかった処理を行うことができる。例えば、前方物体の種類に応じた処理の切り替えを行うことができる。従って、本発明によれば、従来よりも優れたレーダ装置を提供することができる。
【００１７】
尚、確率分布記憶手段は、上記種類毎の確率分布として少なくとも、反射元が金属物であるときの確率分布、及び反射元が人であるときの確率分布を記憶した構成にされるとよい。そして、判別手段は、これらの確率分布に基づき、上記反射元の種類として、少なくとも金属物及び人の夫々を判別可能な構成にされるとよい（請求項２）。
【００１８】
このようにレーダ装置を構成すれば、例えばレーダ装置による判別結果を車両制御に応用する場合に、前方物体が金属物（車両やガードレール等）であるのか人であるのかに応じて、車両制御の内容を切り替えることができ、危険回避等に大変役立つ。特に、従来技術では衝突回避が不可能となったときに、人を回避してガードレールに衝突するように車両制御を行うといったことは不可能であったが、本発明のレーダ装置を用いれば、そのような車両制御が実現可能となる。
【００１９】
また、反射元が人であるときの確率分布としては、人が静止しているときの確率分布及び人が移動しているときの確率分布の夫々を確率分布記憶手段に記憶させて、これらの確率分布に基づいて、静止している人及び移動している人を、上記反射元の種類として判別できるように上記判別手段を構成すると一層好ましい。このようにレーダ装置を構成すれば、反射元の種類についての判別結果を危険回避等に一層役立てることができる。
【００２０】
また、判別手段は、次のようにして反射元の種類を判別する構成にすることができる。即ち、判別手段は、判別対象とする反射元の種類毎に、確率分布記憶手段が記憶する当該種類の確率分布に基づき、上記変化量算出手段で算出された変化量Ｙの発生確率であって当該変化量Ｙに対応するレーダ波の反射元が仮に当該種類の反射元であるときの発生確率を導出し、導出した発生確率が最大である反射元の種類を、変化量算出手段で算出された変化量Ｙに対応する反射元の種類として判別する構成にすることができる（請求項３）。
【００２１】
この他、車載型のレーダ装置では、自車両周辺の道路環境を加味して反射元の種類を判別するように上記判別手段を構成すると一層好ましい。
即ち、車載型のレーダ装置には、自車両周辺の道路環境を推定する環境推定手段を設けて、判別手段は、判別対象とする反射元の種類毎に、確率分布記憶手段が記憶する当該種類の確率分布と、環境推定手段により推定された道路環境と、に基づき、変化量算出手段で算出された変化量Ｙに対応するレーダ波の反射元が当該種類の反射元である確度を算出し、反射元の種類毎の確度に基づき、変化量算出手段で算出された変化量Ｙに対応するレーダ波の反射元の種類を判別する構成にされると好ましい（請求項４）。
【００２２】
各種類の反射元が存在する可能性は、道路環境によって異なる。換言すれば、道路環境を指標にすれば、各種類の反射元が存在する可能性をある程度、事前予測することができる。例えば、車両が、街中を走行しているのか、郊外を走行しているのか、高速道路を走行しているのかによって、車両周辺に人が存在する可能性は大きく変化する。従って、このように道路環境を加味して、反射元の種類を判別するようにレーダ装置を構成すると、一層高精度に反射元の種類を判別することができる。
【００２３】
尚、環境推定手段による道路環境の推定結果を用いる場合、判別手段は、具体的に次のように構成することができる。即ち、判別手段は、判別対象とする反射元の種類毎に、確率分布記憶手段が記憶する当該種類の確率分布に基づき、変化量算出手段で算出された変化量Ｙの発生確率であってレーダ波の反射元が当該種類の反射元であると仮定した場合での発生確率を導出し、この導出した発生確率に環境推定手段により推定された道路環境に対応した補正を加えることで、変化量算出手段で算出された変化量Ｙに対応するレーダ波の反射元が当該種類の反射元である確度を算出する構成にすることができる（請求項５）。例えば、上記導出した発生確率に道路環境に対応した補正係数を掛けて、変化量Ｙに対応するレーダ波の反射元が上記種類の反射元である確度を算出するといった具合である。
【００２４】
また、環境推定手段は、各区域の道路環境を表すデータを記憶する道路環境データベースから直接又は間接的にデータ取得可能な構成であることを前提として、次のように構成することができる。即ち、環境推定手段は、道路環境データベースから自車両の現在位置に対応する区域の道路環境を表すデータを取得し、この取得データが表す道路環境を、自車両周辺の道路環境であると推定する構成にすることができる（請求項６）。
【００２５】
例えば、環境推定手段は、上記道路環境データベースとしての機能を有する道路地図データベースから、ＧＰＳ受信機等の位置検出器にて検出される自車両の現在位置に対応する区域の道路環境を表すデータ（例えば、街中、郊外、高速道路等の道路環境のカテゴリを表すデータ）を取得し、この取得データに基づき、自車両周辺の道路環境を推定する構成することができる。車両には、道路地図データベースやＧＰＳ受信機等の位置検出器が備えられている例も多いため、環境推定手段をこのように構成すれば、車内の機器を利用して、高精度に自車両周辺の道路環境を推定することができる。
【００２６】
但し、このような手法では、道路環境データベースをレーダ装置や車両内に設けたり、外部の道路環境データベースからデータを取得するための通信インタフェースをレーダ装置に設けたりする必要があるため、廉価に道路環境を加味した反射元種類の判別動作を実現するのが難しい。
【００２７】
そこで、環境推定手段は、自車両の走行速度に基づき、道路環境を推定する構成にされてもよい（請求項７）。交差点の多い環境において、車両は低速で走行することが多く、交差点の少ない環境においては、高速に走行することが多くなる。また、交差点のない高速道路等では、車両は更に高速に走行する。よって、車両の走行速度に基づけば、現在の道路環境が、どのような種類の道路環境であるのかを推定することができ、道路環境データベースを必要とせず、廉価に道路環境を加味した反射元種類の判別動作を実現することができる。
【００２８】
また、本発明は、ＦＭＣＷ（Frequency Modulated Continuous Wave）方式にてレーダ波の送受信を行うＦＭＣＷレーダ装置や、ＳＦＷ（Stepped Frequency Waveforms）方式にてレーダ波の送受信を行うＳＦＷレーダ装置に適用することができる（請求項８）。この種のレーダ装置によれば、レーダ波の反射元までの距離、反射元の方位、反射元の移動速度等を検出することができるので、本発明を適用して、更に反射元の種類を判別できるようにレーダ装置を構成すれば、その判別結果を、例えば当該レーダ装置を搭載する車両の走行制御等に役立てることができる。
【００２９】
また、本発明を、ＦＭＣＷ方式やＳＦＷ方式にて、レーダ波を送信アンテナから発射し、当該レーダ波の反射波をアレーアンテナで受信する送受信手段を有したＦＭＣＷレーダ装置やＳＦＷレーダ装置に適用する場合、当該レーダ装置には、アレーアンテナを構成する各アンテナ素子からの受信信号にレーダ波の送信信号を混合して、アンテナ素子毎のビート信号を生成するビート信号生成手段と、ビート信号生成手段により生成されたアンテナ素子毎のビート信号に基づき、レーダ波の反射元までの距離及び当該反射元の方位を検出する信号解析手段とを設け、更に上記検出手段、変化量算出手段、及び、判別手段を次のように構成することができる。
【００３０】
即ち、検出手段は、ビート信号生成手段により生成されたアンテナ素子毎のビート信号に基づき、距離及び方位の異なる上記反射元毎に、当該反射元から到来した反射波の信号強度を検出し、変化量算出手段は、上記反射元毎に、当該反射元から到来した反射波の上記検出手段により検出された信号強度Ｐと、当該信号強度Ｐに対応する反射波の受信時点から観測周期Ｔｓｗ遡った時点で当該反射元から到来し送受信手段により受信された反射波の上記検出手段により検出された信号強度Ｐｂと、に基づき、信号強度Ｐｂに対する信号強度Ｐの変化量Ｙを算出し、判別手段は、上記反射元毎に、当該反射元の種類を、上記変化量算出手段で算出された変化量Ｙと、確率分布記憶手段が記憶する上記種類毎の確率分布とに基づき判別する構成にすることができる（請求項９）。
【００３１】
このレーダ装置によれば、前方に存在する複数の物体の位置（距離及び方位）並びに種類を特定することができるので、レーダ装置周辺の状況を詳細に把握することができる。従って、この判別結果を例えば車両制御に利用すれば、前方に位置する物体の位置及び種類に応じて、適切な車両制御を実現でき、危険回避等に大変役立つ。
【図面の簡単な説明】
【００３２】
【図１】レーダ装置１の構成を表すブロック図である。
【図２】レーダ波の送受信号Ｓｓ，Ｓｒ及びビート信号ＢＴを示したグラフ（ａ）並びに、ビート信号ＢＴのパワースペクトルを示したグラフ（ｂ）である。
【図３】前方物体方位とＭＵＳＩＣスペクトルとの対応関係を示した説明図である。
【図４】確率分布テーブルの構成を表す図である。
【図５】信号処理部３０が実行するメインルーチンを表すフローチャートである。
【図６】信号処理部３０が実行する物標推定処理を表すフローチャートである。
【図７】第一変形例のレーダ装置２の構成を表す図である。
【図８】第一変形例の物標推定処理の一部を表すフローチャートである。
【図９】第二変形例の物標推定処理の一部を表すフローチャートである。
【発明を実施するための形態】
【００３３】
以下に本発明の実施例について、図面と共に説明する。
本実施例のレーダ装置１は、ＦＭＣＷ方式によりレーダ波を送受信する車載型のＦＭＣＷレーダ装置であり、図１に示すように、周波数が時間に対して直線的に漸次増減するミリ波帯の高周波信号を生成する発振器１１と、発振器１１が生成する高周波信号を増幅する増幅器１３と、増幅器１３の出力を送信信号Ｓｓとローカル信号Ｌとに電力分配する分配器１５と、送信信号Ｓｓに応じたレーダ波を発射する送信アンテナ１７と、を備える。尚、図２（ａ）は、レーダ波の送信信号Ｓｓ及び受信信号Ｓｒを示したグラフ（上段）及びビート信号ＢＴを示したグラフ（下段）である。図２（ａ）に示すように、本実施例のレーダ装置１は、レーダ波を、三角形状に周波数変調して発射する。
【００３４】
更に、レーダ装置１は、前方物体にて反射されたレーダ波（反射波）を受信するＫ個のアンテナ素子が一列に配置されたリニアアレーアンテナとしての受信アンテナ１９と、受信アンテナ１９を構成するアンテナ素子の一つを順次選択し、選択したアンテナ素子からの受信信号Ｓｒを後段に供給する受信スイッチ２１と、受信スイッチ２１から供給される受信信号Ｓｒを増幅する増幅器２３と、増幅器２３にて増幅された受信信号Ｓｒ及びローカル信号Ｌを混合して、ビート信号ＢＴを生成（図２（ａ）下段参照）するミキサ２５と、ミキサ２５が生成したビート信号ＢＴから不要な信号成分を除去するフィルタ２７と、フィルタ２７の出力をサンプリングし、ディジタルデータに変換するＡ／Ｄ変換器２９とを備える。以下では、Ｋ個のアンテナ素子の夫々を番号付けして、第ｋアンテナ素子（ｋ＝１，２，…，Ｋ）と表現する。
【００３５】
この他、レーダ装置１は、信号処理や装置内各部の制御を行うための構成として、マイクロコンピュータを中心に構成される信号処理部３０と、確率分布テーブル（詳細後述）を記憶する記憶部４０とを備える。
【００３６】
信号処理部３０は、発振器１１の起動／停止等を制御すると共に、マイクロコンピュータでのプログラム実行により、Ａ／Ｄ変換器２９から入力されるビート信号ＢＴのサンプリングデータ（ディジタルデータ）を用いた信号処理や、当該信号処理により得られる前方物体の位置・相対速度・方位・種類等の情報を、図示しない通信インタフェースを介して車内の電子制御装置（ＥＣＵ）に送信する処理を行う。例えば、レーダ装置１は、車間制御を行う電子制御装置である車間制御ＥＣＵと通信ケーブルを通じてシリアル通信可能に接続されて、車間制御ＥＣＵに対して上記情報を提供する。この他、レーダ装置１は、車内ＬＡＮに接続されて、同一の車内ＬＡＮに接続された電子制御装置に上記情報を提供する構成にされてもよい。
【００３７】
このレーダ装置１では、信号処理部３０からの指令に従って発振器１１が起動し、当該起動により発振器１１が生成した高周波信号は、増幅器１３にて増幅された後、分配器１５に入力され、分配器１５によって電力分配される。これにより、レーダ装置１では、送信信号Ｓｓ及びローカル信号Ｌが生成され、送信信号Ｓｓは、送信アンテナ１７を介し、レーダ波として発射される。
【００３８】
一方、送信アンテナ１７から発射され前方物体に反射して戻ってきたレーダ波（反射波）は、受信アンテナ１９を構成する各アンテナ素子にて受信され、各アンテナ素子からは、受信スイッチ２１に向けて、その受信信号Ｓｒが出力される。
【００３９】
また、受信スイッチ２１からは、受信スイッチ２１によって選択された第ｋアンテナ素子（ｋ＝１，…，Ｋ）の受信信号Ｓｒのみが増幅器２３に出力され、増幅器２３で増幅された受信信号は、ミキサ２５に供給される。
【００４０】
ミキサ２５では、受信信号Ｓｒに分配器１５からのローカル信号Ｌが混合されて、信号Ｓｒ，Ｓｓの差の周波数成分であるビート信号ＢＴが生成される。このビート信号ＢＴは、フィルタ２７にて不要な信号成分が除去された後、Ａ／Ｄ変換器２９を通じて、ディジタルデータとして信号処理部３０に取り込まれる。
【００４１】
但し、受信スイッチ２１は、レーダ波の変調周期Ｔｍの間に、全てのアンテナ素子ＡＮ＿１〜ＡＮ＿Ｋを所定回ずつ選択するように、切り替えられる。そして、Ａ／Ｄ変換器２９は、この切替タイミングに同期してデータサンプリングを行う。
【００４２】
信号処理部３０は、このようにＡ／Ｄ変換器２９にてサンプリングされるビート信号ＢＴの当該サンプリングデータ（ディジタルデータ）を、アンテナ素子毎のサンプリングデータに分離し、当該アンテナ素子毎のビート信号ＢＴのサンプリングデータを解析することにより、周知の手法で、反射波の発生元である前方物体までの距離Ｒ、自車両に対する前方物体の相対速度Ｖ、及び、自車両の進行方向を基準とした前方物体の方位θを推定し、更には、後述する本実施例に特有な手法で、前方物体の種類を判別する。
【００４３】
送信信号Ｓｓに基づくレーダ波を送信アンテナ１７が送信したことに起因して、受信アンテナ１９が、反射波を受信すると、受信信号Ｓｒは、図２上段に点線で示すように、レーダ波が前方物体との間を往復するのに要した時間、即ち、前方物体までの距離Ｒに応じた時間Ｔｒ分遅延し、前方物体との相対速度Ｖに応じた周波数ｆｄ分ドップラシフトする。そして、この時間Ｔｒ及び周波数ｆｄの情報は、図２（ａ）下段からも理解できるように、送信信号Ｓｓの周波数が上昇する上り区間及び送信信号Ｓｓの周波数が下降する下り区間の各区間におけるビート信号ＢＴの周波数の違いとして表れる。
【００４４】
このため、信号処理部３０は、１変調周期Ｔｍにおける上り区間及び下り区間（図２（ａ）上段参照）の夫々において、上記サンプリングデータ（ディジタルデータ）として得られるアンテナ素子毎のビート信号ＢＴを周波数変換（フーリエ変換）し、区間毎に、アンテナ素子毎のパワースペクトルを求める（図２（ｂ）参照）。
【００４５】
そして、各アンテナ素子のパワースペクトルに共通するピーク周波数ｆｐを、区間毎に検出し、上り区間のピーク周波数ｆｐ＝ｆｂ１と、下り区間のピーク周波数ｆｐ＝ｆｂ２とから、前方物体との距離Ｒ及び相対速度Ｖを、式（１）〜（４）に従って算出する。但し、ｃは電波伝搬速度，ｆｍは送信信号Ｓｓの変調周波数、Δｆは送信信号Ｓｓの周波数変動幅、ｆ０は送信信号Ｓｓの中心周波数である。
【００４６】
【数１】
【００４７】
一方、受信アンテナ１９の各アンテナ素子が受信する反射波には、到来方向に応じた位相差が生じる。信号処理部３０は、ビート信号ＢＴに含まれるこのような位相差の情報に基づいて、前方物体の方位θを推定する。方位推定の手法としては、ディジタルビームフォーミング（ＤＢＦ）法、ＣＡＰＯＮ法、ＥＳＰＲＩＴ法やＭＵＳＩＣ法など種々の方法が知られているが、例えば、ＭＵＳＩＣ法によって方位推定を行う場合、信号処理部３０は、上記検出されたピーク周波数ｆｐ毎に、各アンテナ素子のビート信号ＢＴの周波数変換値（フーリエ変換値）であって、当該ピーク周波数ｆｐでの周波数変換値を用いて自己相関行列を生成し、この自己相関行列の固有値・固有ベクトルに基づいてＭＵＳＩＣスペクトルを算出する。
【００４８】
このようにＭＵＳＩＣスペクトルを算出すると、例えば、自車両前方に先行車両が存在する場合には、図３に示すように、先行車両に対応する方位に鋭いピークが立つ。本実施例では、このようなＭＵＳＩＣスペクトルのピークを抽出し、各ピーク周波数ｆｐに対応するレーダ波の反射波成分の到来方位θ1，…，θMを推定する。
【００４９】
尚、上述したように、距離Ｒ及び相対速度Ｖの算出の際には、上り区間のピーク周波数ｆｂ１と下り区間のピーク周波数ｆｂ２とのペアを特定する必要があるが、前方物体が複数ある場合には、区間毎に複数のピーク周波数ｆｐが検出されることになる。このような区間毎に複数のピーク周波数ｆｐが検出される環境での上り区間のピーク周波数ｆｂ１と下り区間のピーク周波数ｆｂ２とのペアマッチングは、例えば、ピーク周波数ｆｐ毎の方位θ1，…，θMの推定結果を指標に、反射波の到来方位が対応するピーク周波数ｆｂ１，ｆｂ２の組をマッチングする。このような処理により、信号処理部３０は、前方物体の距離Ｒ及び相対速度Ｖ及び方位θを推定する。
【００５０】
この他、信号処理部３０は、各方位θ1，…，θMから到来する反射波のパワー（受信電力）Ｐ1，…，ＰMを所定の演算式に従って算出し、各パワーＰ1，…，ＰMの１変調周期Ｔｍ前のパワーＰｂ1，…，ＰｂMからの変化量Ｙを算出することで、当該変化量Ｙを指標に、対応する方位θから到来するレーダ波を反射した前方物体の種類（以下、「反射元種類」とも表現する。）を判別する。
【００５１】
反射波の位相は、反射面の位置に応じて変化するため、反射面が揺らげば（振動すれば）、その揺らぎによる反射面の位置変化が反射波の位相変化に現れることになる。一方、レーダ波を反射する前方物体は、微小な複数の反射面の集合にて構成されるものと解釈することができ、その反射面の揺らぎ方は、前方物体が金属物のような表面が硬いものなのか、人のような表面が柔らかいものなのか、各部位（例えば腕や脚）が変位するものなのか等によって異なる。そして、揺らぎ方が異なれば、変化量Ｙの確率分布に差異が生じる。本実施例では、このような現象を利用し、上記変化量Ｙに基づいて、反射元の種類を判別する。
【００５２】
ここで、一つの物体からの反射波について考えてみると、時刻ｔでの反射波は、次式により表現することができる。
【００５３】
【数２】
【００５４】
但し、上式は、物体の反射面が複数の独立に動く等ＲＣＳの反射面から構成されることを前提としたものであり、Ｃは、レーダ波の送信ゲイン、受信ゲイン、送信電力及びＲＣＳ値によって定まる値である。また、ここでは、時刻ｔにおけるｉ番目の反射面までの距離Ｒi（ｔ）を、Ｒi（ｔ）＝Ｒ（ｔ）＋ΔＲi（ｔ）で表現する。Ｒ（ｔ）は、複数反射面までの距離の平均値を表し、ΔＲi（ｔ）は、平均値Ｒ（ｔ）からの距離差分を表す。また、ｆは、レーダ波の送信周波数を表し、λは、レーダ波の送信波長を表す。尚、近似式は、ΔＲi（ｔ）が微小であることを利用して、分母のＲ（ｔ）＋ΔＲi（ｔ）を、Ｒ（ｔ）に近似したものである。
【００５５】
また、時刻ｔから、微小な観測周期Ｔｓｗ分進んだ時点での反射波は、次式により表現することができる。
【００５６】
【数３】
【００５７】
但し、上式で用いるｖ（ｔ）は、レーダ波を反射した前方物体の相対速度を表し、Δｖi（ｔ）は、ｉ番目の反射面での相対速度Ｖ（ｔ）に対する速度差分、即ち、揺らぎ量を表す。近似式では、時刻ｔ＋Ｔｓｗでの周波数ｆTswを、時刻ｔでの周波数ｆと同一値に近似している。
【００５８】
時刻ｔ及び時刻ｔ＋Ｔｓｗでの反射波は、上記のように近似できることから、例えば、変化量Ｙとして、レーダ波の受信信号のパワー比を算出すると、変化量Ｙは次式で表現される。
【００５９】
【数４】
【００６０】
上式では、パワー比を対数増幅したものを変化量Ｙとして表す。変化量Ｙとして、時刻ｔ＋Ｔｓｗでのレーダ波の受信信号のパワーと、時刻ｔでのレーダ波の受信信号のパワーとの比を算出すれば、変化量Ｙは、ΔＲi（ｔ）及びΔＶi（ｔ）によって変化する値として算出される。このため、変化量Ｙを、予め実験で得た反射元種類毎の変化量Ｙの確率分布と比較することで、反射元（前方物体）の種類を判別することができるのである。
【００６１】
但し、観測周期Ｔｓｗの間に生じる各反射面からの反射波成分の位相変化量Δφがπ以上となると、各反射面の反射波成分が混合する受信アンテナ１９の反射波の受信信号では、その振幅がレイリー分布に従うようになり、反射面の揺らぎを観測することが難しい。
【００６２】
そこで、本実施例では、レーダ波の送信波長λ及び観測周期Ｔｓｗを、判別対象とする各種類の前方物体のいずれにおいても、その前方物体での各反射面から反射する反射波の位相変化量Δφがπを超えないように、送信波長λ及び観測周期Ｔｓｗを定める。
【００６３】
ここで、送信波長λ及び観測周期Ｔｓｗの設定条件について詳述する。時刻ｔにおけるｉ番目の反射面からの反射波の位相をφi（ｔ）で表現し、時刻ｔ＋Ｔｓｗにおけるｉ番目の反射面からの反射波の位相をφi（ｔ＋Ｔｓｗ）で表現すると、観測周期Ｔｓｗの間でのｉ番目の反射面からの反射波の位相変化量Δφiは、Δφi＝φi（ｔ＋Ｔｓｗ）−φi（ｔ）で表現される。そして、観測周期Ｔｓｗでの各反射面からの反射波の位相変化量Δφiが独立に−πからπまでの一様分布に従うとき、これら反射波の混合波の振幅分布は、レイリー分布に従う。
【００６４】
換言すれば、各反射面からの反射波の位相Δφiがπを超えないように送信波長λ及び観測周期Ｔｓｗを定めれば、前方物体の揺らぎを、受信アンテナ１９による反射波の受信信号にて観測できることになる。即ち、送信波長λ及び観測周期Ｔｓｗについては、Δφiの最大値ｍａｘ｛Δφi｝がπ未満となるように定めればよい。
【００６５】
【数５】
【００６６】
一方、Δφiは、次式によって表すことができる。
【００６７】
【数６】
【００６８】
各反射面の揺らぎ量ΔＶi（ｔ）は、平均ゼロ、分散σv2の正規分布Ｎ（０，σv2）に従うと近似的にみなすことができ、分散σv2で揺らぐ前方物体の揺らぎを、３σvの範囲（９９．７％の確率で起こる範囲）で観測できるように、送信波長λ及び観測周期Ｔｓｗを定める場合には、次の関係式を満足すればよい。
【００６９】
【数７】
【００７０】
換言すると、送信波長λ及び観測周期Ｔｓｗは、次の関係式を満足すればよい。
【００７１】
【数８】
【００７２】
本実施例では、このような関係式に従って、判別対象の種類の前方物体の内、揺らぎが最も大きい種類の前方物体の分散σv2を基準に、送信波長λ及び観測周期Ｔｓｗを定める。尚、揺らぎの大きさ（分散σv2）については、実験により計測すればよい。
【００７３】
本実施例では、このような関係式に従って、送信波長λ及び観測周期Ｔｓｗを定めつつ、記憶部４０に、判別対象の反射元種類毎の変化量Ｙの確率分布を表すテーブルを用意する。これによって、変化量Ｙから反射元（前方物体）の種類を判別することができるようにする。
【００７４】
図４には、記憶部４０が記憶する確率分布テーブルの構成を示す。図４に示すように、本実施例の確率分布テーブルは、判別対象の反射元種類毎に、レーダ波の反射元が当該種類の反射元であったときに得られる変化量Ｙの確率分布が記述されてなるテーブルである。具体的に、この確率分布テーブルには、離散的な変化量Ｙの各値に対して、変化量Ｙがその値以下となる確率（累積確率）が、判別対象の反射元種類毎に記述されている。
【００７５】
このテーブルについては、例えば、判別対象の前方物体の種類毎に、当該種類の前方物体に対してレーダ波を発射して反射波の計測データを取得し、この反射波の計測データから得られた変化量Ｙの確率分布を、テーブルに記述することで作成することができる。また、別手法としては、反射波の計測データを取得する一方、シミュレーションにより、分散σv2を変化させて各分散σv2での変化量Ｙの確率分布を求め、反射波の計測データと最も近い分散σv2での変化量Ｙの確率分布を、テーブルに記述することで作成することができる。
【００７６】
ちなみに、本実施例では、確率分布テーブルに、判別対象の反射元種類毎の累積確率分布を記述したが、前方物体の種類を判別する際には、受信したレーダ波から得られた変化量Ｙの累積確率と、隣接する変化量Ｙに対して記述された累積確率との差分を採ることで、累積確率分布を微分し、変化量Ｙに対応する確率密度を求める。このため、確率分布テーブルには、累積確率分布に代えて確率密度分布を記述してもよい。この他、記憶部４０には、確率分布テーブルに代えて、累積分布関数及び確率密度関数の情報を記憶させることも可能である。
【００７７】
また、本実施例では、前方物体として、車両やガードレール、人等を判別して、車両制御に役立てることを考えている。このため、記憶部４０が記憶する確率分布テーブルには、反射元種類として、金属物及び人の夫々についての確率分布が記述されている。特に、本実施例の確率分布テーブルには、「静止している人」及び「移動している人」の夫々についての確率分布が記述されており、レーダ装置１は、前方物体として、「金属物」、「静止している人」、「移動している人」の夫々を判別可能に構成されている。このようなレーダ装置１によれば、前方に移動する人がいると判定された場合には、車両速度を緩める等の車両制御を行うことができて、危険回避等に大変役立つ。この他、本実施例のレーダ装置１によれば、ガードレールと人を見分けることができるので、例えば、衝突が回避不可能となったときに、人を回避してガードレールに衝突するといったような従来では不可能な車両制御を実現することが可能である。
【００７８】
続いては、以上に概略的に説明した信号処理部３０の処理動作を、図５及び図６を用いて詳細に説明する。信号処理部３０は、図５に示すメインルーチンを変調周期Ｔｍ毎に繰返し実行することにより、変調周期Ｔｍ毎に、一連の信号処理（Ｓ１１０〜Ｓ１５０）を実行して、前方物体までの距離Ｒ、前方物体との相対速度Ｖ及び前方物体の方位θを推定すると共に、前方物体の種類を判別する。
【００７９】
具体的には、まず１変調周期Ｔｍにおける上り区間を処理対象区間として、この処理対象区間（上り区間）に受信した反射波よって生成されたビート信号ＢＴのサンプリングデータに基づき、図６に示す物標推定処理を実行する（Ｓ１１０）。その後、処理対象区間を下り区間に変更し、この処理対象区間（下り区間）に受信した反射波によって生成されたビート信号ＢＴのサンプリングデータに基づき、図６に示す物標推定処理を実行する（Ｓ１２０）。
【００８０】
物標推定処理では、処理対象区間におけるアンテナ素子毎のビート信号ＢＴのサンプリングデータを取得し（Ｓ２１０）、その後、これらサンプリングデータに基づき、各アンテナ素子のビート信号ＢＴを周波数変換して、アンテナ素子毎のパワースペクトルを算出する（Ｓ２２０）。その後、各アンテナ素子共通のピーク周波数ｆｐを検出する。ここでは、アンテナ素子毎のパワースペクトルに共通するピークであって予め定められた閾値以上のピークを検出し、検出した閾値以上の各ピークの周波数を、ピーク周波数ｆｐとして検出する（Ｓ２３０）。
【００８１】
ピーク周波数ｆｐの検出後には、検出したピーク周波数ｆｐの一つを処理対象周波数に選択し（Ｓ２４０）、各アンテナ素子における処理対象周波数の周波数変換値（フーリエ変換値）ｘkを配列してなる受信ベクトルＸ＝（ｘ1，ｘ2，…ｘK）Tを生成して、この受信ベクトルＸについての自己相関行列Ｒｘｘ１＝Ｅ［ＸＸH］及びＲｘｘ２＝ＸＸHを算出する（Ｓ２５０）。ここで、Ｔは、ベクトル転置を示し、Ｈは、複素共役転置を示す。また、ｘkは、ｋ番目のアンテナ素子の受信信号から生成された処理対象区間のビート信号ＢＴの周波数変換値であって、処理対象周波数での周波数変換値を示す。従来知られる手法では、方位推定のために自己相関行列Ｒｘｘとして、ＸＸHの期待値Ｅ［ＸＸH］を算出するが、本実施例では、上述の変化量Ｙを算出する際に、各方位から到来する反射波の瞬時パワーを計算する必要があるので、自己相関行列Ｒｘｘ１として、ＸＸHの期待値Ｅ［ＸＸH］を算出する一方、パワー推定のために自己相関行列Ｒｘｘ２として、最新の受信ベクトルＸを用いたＸＸHの瞬時値を算出する。期待値Ｅ［ＸＸH］としては、ＸＸHの最新値を含む過去所定個のＸＸHの平均値を算出することができる。
【００８２】
その後、自己相関行列Ｒｘｘ１についての固有値λ1，…，λK（但し、λ1≧λ2≧…λk）及び固有ベクトルｅ1，…，ｅKを求め、熱雑音電力に対応する閾値λthより大きい固有値の数から到来波数Ｍを推定する（Ｓ２６０）。
【００８３】
そして、熱雑音電力以下となる（Ｋ−Ｍ）個の固有値λM+1，…，λKに対応した固有ベクトルｅM+1，…，ｅKからなる雑音固有ベクトルＥNと、方位θに対するアレーアンテナの複素応答、即ち、ステアリングベクトルａ（θ）とから、角度スペクトルとして、下式の評価関数ＰMU（θ）で表されるＭＵＳＩＣスペクトルを求める。
【００８４】
【数９】
【００８５】
そして、この評価関数ＰMU（θ）で表されるＭＵＳＩＣスペクトルから、到来波数Ｍ分のピークを抽出し、各ピークに対応する方位θ1，…，θMを、反射波の到来方位θ1，…，θM、換言すれば、レーダ波の反射元である前方物体（物標）の方位として推定する（Ｓ２７０）。
【００８６】
また、推定した各方位θ1，…，θMから到来する反射波のパワー（受信電力）を次のように推定して、これを一時記憶する（Ｓ２８０）。
具体的には、まずパワー（受信電力）を推定する対象であるＳ２７０で推定した反射波の各到来方位θ1，…，θMに対応するステアリングベクトルａ（θ1），…，ａ（θM）を用いて、次の方向行列Ａ（所謂アレー応答行列）を生成する。
【００８７】
【数１０】
【００８８】
そして、この方向行列Ａ及び瞬時値としての自己相関行列Ｒｘｘ２を用いて、次式で表される行列Ｓの対角成分を算出する。
【００８９】
【数１１】
【００９０】
行列Ｓの第ｍ対角成分は、方位θmの受信電力Ｐmに対応する（但し、ｍ＝１，…，Ｍ）。即ち、ここでは、行列Ｓの対角成分を算出することにより、Ｓ２７０で推定した各到来方位θ1，…，θMから到来した反射波のパワー（受信電力）を推定する。但し、Ｓ２８０では、Ｓ２７０で反射波の到来方位として推定したものの、上式により算出されたパワーが閾値以下の方位を、反射波の到来方位から外す処理を併せて行う。
【００９１】
また、この処理を終えると、信号処理部３０は、反射波の到来方位θ1，…，θMの一つを、処理対象方位に選択し（Ｓ２９０）、選択した処理対象方位から到来した反射波のパワーの前回推定値が存在するか否かを判断する（Ｓ３００）。
【００９２】
即ち、ここでは、１変調周期Ｔｍ前の処理対象区間における反射波の受信信号に基づき、今回の処理対象周波数と同一のピーク周波数ｆｐでの方位推定及びパワー推定が行われて、処理対象方位と同一方位から到来した反射波のパワー推定がなされているか否かを判断し、なされている場合には、パワーの前回推定値が存在すると判断し、なされていない場合には、パワーの前回推定値が存在しないと判断する。そして、パワーの前回推定値が存在しないと判断すると（Ｓ３００でＮｏ）、Ｓ３４０に移行し、パワーの前回推定値が存在すると判断すると（Ｓ３００でＹｅｓ）、Ｓ３１０に移行する。
【００９３】
Ｓ３１０に移行すると、信号処理部３０は、処理対象方位について今回推定したパワーの前回推定値からの変化量Ｙを算出する。ここで、今回推定したパワーの推定値をＰと表現し、前回推定したパワーの推定値をＰｂと表現すると、変化量Ｙについては、次式に従って算出することになる。尚、ここで対数増幅を行っているのは、パワーの推定値をデシベル単位で記憶保持しておくことで、差分により変化量Ｙを簡単に算出することができるためである。
【００９４】
【数１２】
【００９５】
ちなみに、本実施例では、このように、前回のパワー推定値として、１変調周期Ｔｍ前の値を用いるので、上述した観測周期Ｔｓｗは、変調周期Ｔｍに一致する。
また、この処理を終えると、信号処理部３０は、Ｓ３２０に移行し、判別対象とする反射元種類毎に、Ｓ３１０で算出された変化量Ｙが仮に当該種類の反射元から到来したレーダ波に基づき算出された変化量Ｙである場合に当該変化量Ｙが得られる確率Ｐｒ（以下、変化量Ｙの発生確率Ｐｒとも言う。）を、記憶部４０が記憶する確率分布テーブルに記された該当種類の確率分布の情報に基づいて導出する（Ｓ３２０）。尚、上述したように確率分布テーブルには、確率分布として各変化量Ｙでの累積確率が離散的に記されていることから、ここでは、Ｓ３１０で算出された変化量Ｙに対応する累積確率と隣接する変化量Ｙの累積確率との差分を算出することで、変化量Ｙに対応する累積確率の微分である確率密度を算出する。即ち、Ｓ３２０では、上記手法で、Ｓ３１０で算出された変化量Ｙの確率密度を算出することで、変化量Ｙの発生確率Ｐｒを反射元種類毎に算出する。
【００９６】
このようにして、判別対象とする反射元種類毎の発生確率Ｐｒを算出すると、信号処理部３０は、上記種類毎の発生確率Ｐｒの中から、最も値の高いものを抽出し、この最も発生確率Ｐｒの高い反射元の種類を、処理対象方位から到来した反射波の反射元の種類であると判別する（Ｓ３３０）。
【００９７】
その後、Ｓ２７０で反射波の到来方位として推定した全方位θ1，…，θM（但しＳ２８０で外されたものを除く。）に関して、当該方位から到来するレーダ波の反射元の種類を判別したか否かを判断し（Ｓ３４０）、判別していない場合には（Ｓ３４０でＮｏ）、Ｓ２９０に移行し、処理対象方位を切り替えて、Ｓ３００以降の処理を実行する。
【００９８】
一方、全方位に関して上記判別を行った場合には（Ｓ３４０でＹｅｓ）、全ピーク周波数ｆｐについて、これを処理対象周波数として選択し、Ｓ２５０以降の処理を実行したか否かを判断し、全ピーク周波数ｆｐに対して実行していなければ（Ｓ３５０でＮｏ）、Ｓ２４０に移行し、処理対象周波数を切り替えて、Ｓ２５０以降の処理を実行する。そして、全ピーク周波数ｆｐに対して実行すると、当該物標推定処理を終了する。
【００９９】
このようにして上り区間及び下り区間の物標推定処理を終了すると、信号処理部３０は、Ｓ１３０（図５参照）に移行し、上り区間でのピークと、下り区間でのピークとのペアマッチングを行い、上り区間でのピーク周波数ｆｐ＝ｆｂ１と下り区間でのピーク周波数ｆｐ＝ｆｂ２との組合せを特定する。そして特定した組合せに基づいて、上述した手法で前方物体までの距離Ｒ及び前方物体との相対速度Ｖを推定する（Ｓ１４０）。
【０１００】
そして、各ピーク周波数ｆｐに対して推定された反射波の到来方位θ1，…，θM及び各到来方位θ1，…，θMから到来するレーダ波の反射元種類の情報を、上記推定された距離Ｒ及び相対速度Ｖの情報と組み合わせて、自車両周囲に位置する各前方物体の距離Ｒ、相対速度Ｖ、方位θ及び種類の情報を、通信インタフェースを通じて外部の電子制御装置に出力する（Ｓ１５０）。信号処理部３０は、このような処理を変調周期Ｔｍ毎に繰返し実行する。
【０１０１】
以上、本実施例のレーダ装置１の構成について説明したが、このレーダ装置１によれば、従来のレーダ装置１とは異なり、前方物体の距離Ｒ、相対速度Ｖ及び方位θに加えて、前方物体の種類を推定することが可能である。特に本実施例によれば、前方物体の種類として、「金属物（車両やガードレール等）」、「静止している人」及び「移動している人」の夫々を判別することができる。
【０１０２】
従って、本実施例によれば、前方物体が金属物か人であるかによって車両制御のパターンを切り替えることができ、更には人が移動中であるか静止しているかによっても車両制御のパターンを切り替えることができて、上述したように危険回避等に大変役立つ。
【０１０３】
さて、上記実施例では、確率分布テーブルが示す確率分布とレーダ波の受信信号に基づく変化量Ｙとの比較により、反射元の種類を判別するようにしたが、各種類の反射元が存在する可能性については、道路環境によっても事前にある程度予測することができる。例えば、街中では、人が存在する可能性が高く、郊外では、街中よりも人が存在する可能性が低く、高速道路においては、車外に人が存在することは通常ありえない。
【０１０４】
そこで、変化量Ｙの発生確率Ｐｒと道路環境とから、変化量Ｙに対応するレーダ波の反射元が判別対象の各種類の反射元である確度を評価し、この確度に基づいて、反射元の種類を判別するように、レーダ装置を構成してもよい（第一変形例及び第二変形例）。
［第一変形例］
続いて、第一変形例について説明する。第一変形例のレーダ装置２は、記憶部４０に対応する記憶部４１が確率分布テーブルに加えて補正テーブルを記憶する点、信号処理部３０に対応する信号処理部３１が車内ＬＡＮを通じてナビゲーションＥＣＵ１００と通信可能に構成されている点、及び、信号処理部３１が実行する物標推定処理の一部が上記実施例と異なる点を除けば、基本的に上記実施例のレーダ装置１と同一構成にされている。従って、以下では、第一変形例のレーダ装置２の説明として、上記実施例のレーダ装置１と同一構成の説明を省略する。
【０１０５】
図７に示すように第一変形例のレーダ装置２が備える信号処理部３１は、車内ＬＡＮを通じてナビゲーションＥＣＵ１００と通信可能に接続され、ナビゲーションＥＣＵ１００から、自車両の現在位置に対応する道路環境のカテゴリ情報を取得可能な構成にされている。
【０１０６】
ナビゲーションＥＣＵ１００は、ＧＰＳ受信機等から構成される自車両の現在位置を検出可能な位置検出器１０１、及び、道路地図データベース１０３に接続されており、図７に示す道路環境情報送信処理を繰返し実行して、位置検出器１０１から自車両の現在位置情報を取得し（Ｓ４１０）、この取得情報が示す自車両の現在位置に対応する道路環境のカテゴリを、道路地図データベース１０３が記憶する情報に基づいて判定する（Ｓ４２０）。そして、判定した自車両の現在位置に対応する道路環境のカテゴリ情報をレーダ装置２に車内ＬＡＮを通じて送信する（Ｓ４３０）。
【０１０７】
尚、道路地図データベース１０３には、地図収録エリア内の各地域毎に、当該地域が「街中」及び「郊外」のいずれであるかを表すエリア情報が記憶されている。更に、高速道路に該当する道路に対しては「高速道路」であることを示すエリア情報が記憶されている。即ち、Ｓ４２０では、自車両の現在位置に対応する走行道路を特定して、走行道路が「高速道路」であるか否かを判断し、走行道路が「高速道路」である場合には、自車両の現在位置に対応する道路環境のカテゴリを「高速道路」であると判定し、走行道路が「高速道路」でない場合には、道路地図データベース１０３が記憶する自車両の現在位置に対応する地域のエリア情報が示すカテゴリ（「街中」又は「郊外」）を、自車両の現在位置に対応する道路環境のカテゴリであると判定する。
【０１０８】
この動作によって、ナビゲーションＥＣＵ１００からは、自車両の現在位置に対応する道路環境のカテゴリとして「街中」「郊外」「高速道路」のいずれかのカテゴリを示す情報がレーダ装置２の信号処理部３１に提供される。
【０１０９】
レーダ装置２の信号処理部３１は、このような道路環境のカテゴリ情報をナビゲーションＥＣＵ１００から逐次取得し、物標推定処理のＳ３２０（図６参照）の処理を終えると、Ｓ３３０の処理に代えて、Ｓ５１０〜Ｓ５４０（図８参照）の処理を実行する。
【０１１０】
具体的に、信号処理部３０は、ナビゲーションＥＣＵ１００から得られた最新の上記カテゴリ情報に基づき、現在位置に対応する道路環境のカテゴリを特定する（Ｓ５１０）。そして、記憶部４１が記憶する補正テーブルから、上記特定したカテゴリに対応する補正係数Ｗの組であって反射元種類毎の補正係数Ｗを読み出し（Ｓ５２０）、Ｓ３２０の処理で反射元種類毎に導出した変化量Ｙの発生確率（確率密度）Ｐｒの夫々に対し、対応する反射元種類の補正係数Ｗを乗算する。これによって、反射元種類毎に、変化量Ｙを算出する元となったレーダ波を反射してきた前方物体（反射元）が当該種類の反射元である確度Ｚ＝Ｗ・Ｐｒを算出する（Ｓ５３０）。
【０１１１】
尚、補正テーブルは、図７に示すように、「街中」「郊外」「高速道路」の各道路環境のカテゴリ毎に、補正係数Ｗの組として、反射元種類毎の補正係数Ｗが記述されてなるテーブルである。例えば、「街中」の補正係数Ｗの組としては、「静止している人」に対応する補正係数Ｗ１１＝１／３、「移動している人」に対応する補正係数Ｗ１２＝１／３、「金属物」に対応する補正係数Ｗ１３＝１／３が補正テーブルに記述され、「郊外」の補正係数の組としては、「静止している人」に対応する補正係数Ｗ２１＝０．２５、「移動している人」に対応する補正係数Ｗ２２＝０．２５、「金属物」に対応する補正係数Ｗ２３＝０．５が補正テーブルに記述される。この他、「高速道路」の補正係数の組としては、例えば「静止している人」に対応する補正係数Ｗ３１＝０、「移動している人」に対応する補正係数Ｗ３２＝０、「金属物」に対応する補正係数Ｗ３３＝１が補正テーブルに記述される。
【０１１２】
自車両の現在位置が「郊外」に該当する場合には、「静止している人」の確率分布に基づき導出した変化量Ｙの発生確率ＰｒにＷ２１＝０．２５を乗算して、反射元が「静止している人」である確度Ｚ＝０．２５・Ｐｒを算出するといった具合である。同様に、「移動している人」の確率分布に基づき導出した変化量Ｙの発生確率Ｐｒに、Ｗ２２＝０．２５を乗算して、反射元が「移動している人」である確度Ｚ＝０．２５・Ｐｒを算出し、「金属物」の確率分布に基づき導出した変化量Ｙの発生確率Ｐｒに、Ｗ２３＝０．５を乗算して、反射元が「金属物」である確度Ｚ＝０．５・Ｐｒを算出する。
【０１１３】
そして、算出した確度Ｚが最も大きい反射元種類を、上記処理対象方位から到来した反射波の反射元種類であると判別する（Ｓ５４０）。その後、信号処理部３１は、上記実施例で説明したＳ３４０以降の処理を実行する。
【０１１４】
以上、第一変形例のレーダ装置２について説明したが、この変形例によれば、道路環境を加味して、反射元の種類を判別するようにレーダ装置２を構成したので、道路環境による事前予測に基づいて、高精度に反射元の種類を判別することができる。
【０１１５】
また、第一変形例では、道路地図データベース１０３と位置検出器１０１とを用いて、自車両の道路環境を判定するようにしたが、自車両の走行速度に基づき、自車両の道路環境を判定するように、レーダ装置は、構成されてもよい。
［第二変形例］
続いて、第二変形例のレーダ装置３について説明する。但し、第二変形例のレーダ装置３は、信号処理部３１に対応する信号処理部３２がナビゲーションＥＣＵ１００を通じて上記道路環境のカテゴリ情報を取得する代わりに速度センサ１０５から自車両の速度情報を取得する構成にされている点、及び、信号処理部３２が実行する物標推定処理の一部が上記第一変形例と異なる点を除けば、基本的に第一変形例のレーダ装置２と同一構成にされている。従って、以下では、第二変形例のレーダ装置３の説明として、上述のレーダ装置１，２と同一構成の説明を省略する。
【０１１６】
第二変形例のレーダ装置３の信号処理部３２は、上述したように、道路環境のカテゴリ情報をナビゲーションＥＣＵ１００から取得する代わりに、速度センサ１０５から逐次自車両の速度情報を取得する。一方で、物標推定処理のＳ３２０の処理を終えると、Ｓ５１０〜Ｓ５４０の処理に代えて、Ｓ６１０〜Ｓ６４０（図９参照）の処理を実行する。
【０１１７】
具体的に、信号処理部３０は、Ｓ６１０に移行すると、速度センサ１０５から得られた速度情報に基づき、自車両周辺の道路環境を推定する（Ｓ６１０）。具体的には、現在から過去の所定時点までの期間（例えば現在から過去５分間）における車両の平均速度を算出し、この平均速度を指標に、自車両周囲の道路環境のカテゴリを判定する。
【０１１８】
平均速度が第一の境界値（例えば３０ｋｍ）未満であれば、自車両周辺の道路環境のカテゴリを「街中」と判定し、平均速度が第一の境界値（例えば３０ｋｍ）以上第二の境界値(例えば６０ｋｍ)未満であれば、自車両周辺の道路環境のカテゴリを「郊外」と判定し、平均速度が第二の境界値以上である場合には、自車両周辺の道路環境のカテゴリを「高速道路」と判定するといった具合である。
【０１１９】
このようにして道路環境を推定すると、信号処理部３２は、Ｓ５２０，Ｓ５３０での処理と同様に、記憶部４０が記憶する補正テーブルから、自車両周辺の道路環境のカテゴリに対応する反射元種類毎の補正係数Ｗを読み出し（Ｓ６２０）、反射元種類毎に確率分布テーブルから導出した変化量Ｙの発生確率（確率密度）Ｐｒの夫々に対し、上記読み出した反射種類毎の補正係数Ｗの内、対応する反射元種類の補正係数Ｗを乗算する。これによって、反射元種類毎に、変化量Ｙを算出する元となったレーダ波を反射してきた前方物体（反射元）が当該種類の反射元である確度Ｚ＝Ｗ・Ｐｒを算出する（Ｓ６３０）。
【０１２０】
そして、Ｓ５４０の処理と同様、算出した確度Ｚが最も大きい反射元種類を上記処理対象方位から到来した反射波の反射元種類であると判別し（Ｓ６４０）、その後、上記実施例で説明したＳ３４０以降の処理を実行する。
【０１２１】
以上、第二変形例のレーダ装置２について説明したが、第二変形例によれば、道路環境を加味して、反射元の種類を判別するようにレーダ装置３を構成したので、道路環境による事前予測に基づいて、高精度に反射元種類を判別することができる。特に、第二変形例によれば、車両速度に基づき道路環境を推定するので、道路地図データベース１０３や位置検出器１０１が搭載されていない車両においても、道路環境を加味した反射元種類の判別を行うことができ、安価に様々な車両に適用することができる。
【０１２２】
以上、本発明の実施例について説明したが、レーダ装置１，２，３が備える発振器１１、増幅器１３及び送信アンテナ１７、並びに、受信アンテナ１９、受信スイッチ２１及び増幅器２３は、本発明の送受信手段の一例に対応する。また、信号処理部３０，３１，３２が実行するＳ２８０の処理は、本発明の検出手段により実行される処理の一例に対応する。この他、信号処理部３０，３１，３２が実行するＳ３１０の処理は、本発明の変化量算出手段により実行される処理の一例に対応する。
【０１２３】
この他、記憶部４０における確率分布テーブルの記憶領域は、本発明の確率分布記憶手段の一例に対応し、信号処理部３０が実行するＳ３２０，３３０の処理、又は、信号処理部３１が実行するＳ３２０，Ｓ５２０〜Ｓ５４０の処理、又は、信号処理部３２が実行するＳ３２０，Ｓ６２０〜Ｓ６４０の処理は、本発明の判別手段にて実行される処理の一例に対応する。この他、信号処理部３１が実行するＳ５１０の処理又は信号処理部３２が実行するＳ６１０の処理は、本発明の環境推定手段にて実行される処理の一例に対応する。
【０１２４】
また、レーダ装置１，２，３が備えるミキサ２５は、本発明のビート信号生成手段に対応し、信号処理部３０，３１，３２が実行するＳ１１０〜Ｓ１４０及び物標推定処理におけるＳ２１０〜Ｓ２７０の処理は、本発明の信号解析手段により実行される処理に対応する。
【０１２５】
また、本発明は、上記実施例に限定されるものではなく、種々の態様を採ることができる。例えば、本発明は、ＦＭＣＷレーダ装置に限らず、その他の方式にてレーダ波の送受信を行うレーダ装置に適用することができる。具体的には、ＳＦＷ（Stepped Frequency Waveforms）方式にてレーダ波の送受信を行うＳＦＷレーダ装置に適用することができる。この他、本発明は、車載用のレーダ装置に限らず、その他の用途に用いられるレーダ装置にも適用することができる。また、上記実施例では、人や金属物を判別可能なようにレーダ装置１，２，３を構成したが、判別対象とする反射元の種類は、これに限定されない。また、上記実施例では、道路地図データベース１０３又は速度センサ１０５のいずれか一方を用いて道路環境を推定するようにレーダ装置２，３を構成したが、レーダ装置は、これら両者の情報を用いて道路環境を推定するように構成されてもよい。
【符号の説明】
【０１２６】
１，２，３…レーダ装置、１１…発振器、１３…増幅器、１５…分配器、１７…送信アンテナ、１９…受信アンテナ、２１…受信スイッチ、２３…増幅器、２５…ミキサ、２７…フィルタ、２９…Ａ／Ｄ変換器、３０，３１，３２…信号処理部、４０，４１…記憶部、１００…ナビゲーションＥＣＵ、１０１…位置検出器、１０３…道路地図データベース、１０５…速度センサ
【特許請求の範囲】
【請求項１】
レーダ波を発射し、当該レーダ波の反射波を受信する送受信手段と、
前記送受信手段が受信した前記反射波の受信信号に基づき、前記反射波の信号強度を検出する検出手段と、
前記検出手段により検出された前記反射波の信号強度Ｐと、前記信号強度Ｐに対応する前記反射波の受信時点から所定の観測周期Ｔｓｗ遡った時点で前記送受信手段により受信された前記反射波の前記検出手段により検出された信号強度Ｐｂと、に基づき、前記信号強度Ｐｂに対する前記信号強度Ｐの変化量Ｙを算出する変化量算出手段と、
判別対象とする反射元の種類毎に、前記レーダ波が当該種類の反射元で反射して到来する前記反射波から算出される前記変化量Ｙの確率分布を記憶する確率分布記憶手段と、
前記確率分布記憶手段が記憶する前記種類毎の確率分布に基づき、前記変化量算出手段で算出された前記変化量Ｙに対応する前記レーダ波の反射元の種類を判別する判別手段と、
を備え、前記レーダ波の送信波長λ及び前記観測周期Ｔｓｗは、前記判別対象とする各種類の前記反射元にて発生する揺らぎの特徴が前記変化量Ｙに反映される関係で定められていること
を特徴とするレーダ装置。
【請求項２】
前記確率分布記憶手段は、前記種類毎の確率分布として少なくとも、前記反射元が金属物であるときの前記確率分布、及び前記反射元が人であるときの前記確率分布を記憶し、
前記判別手段は、前記確率分布に基づき、前記反射元の種類として、少なくとも前記金属物及び人の夫々を判別可能な構成にされていること
を特徴とする請求項１記載のレーダ装置。
【請求項３】
前記判別手段は、前記判別対象とする反射元の種類毎に、前記確率分布記憶手段が記憶する当該種類の前記確率分布に基づき、前記変化量算出手段で算出された前記変化量Ｙの発生確率であって当該変化量Ｙに対応する前記レーダ波の反射元が仮に当該種類の反射元であるときの発生確率を導出し、前記導出した発生確率が最大である前記反射元の種類を、前記変化量算出手段で算出された前記変化量Ｙに対応する反射元の種類として判別すること
を特徴とする請求項１又は請求項２記載のレーダ装置。
【請求項４】
車載型のレーダ装置であって、
自車両周辺の道路環境を推定する環境推定手段
を備え、
前記判別手段は、
前記判別対象とする反射元の種類毎に、前記確率分布記憶手段が記憶する当該種類の前記確率分布と、前記環境推定手段により推定された前記道路環境と、に基づき、前記変化量算出手段で算出された前記変化量Ｙに対応する前記レーダ波の反射元が当該種類の反射元である確度を算出し、
前記反射元の種類毎の前記確度に基づき、前記変化量算出手段で算出された前記変化量Ｙに対応する前記レーダ波の反射元の種類を判別すること
を特徴とする請求項１又は請求項２記載のレーダ装置。
【請求項５】
前記判別手段は、前記判別対象とする反射元の種類毎に、前記確率分布記憶手段が記憶する当該種類の前記確率分布に基づき、前記変化量算出手段で算出された前記変化量Ｙの発生確率であって前記レーダ波の反射元が当該種類の反射元であると仮定した場合での発生確率を導出し、当該導出した発生確率に、前記環境推定手段により推定された前記道路環境に対応した補正を加えることで、前記変化量算出手段で算出された前記変化量Ｙに対応する前記レーダ波の反射元が当該種類の反射元である確度を算出すること
を特徴とする請求項４記載のレーダ装置。
【請求項６】
前記環境推定手段は、各区域の道路環境を表すデータを記憶する道路環境データベースから、自車両の現在位置に対応する区域の道路環境を表すデータを取得して、この取得データが表す道路環境を、自車両周辺の道路環境であると推定すること
を特徴とする請求項４又は請求項５記載のレーダ装置。
【請求項７】
前記環境推定手段は、自車両の走行速度に基づき、前記道路環境を推定すること
を特徴とする請求項４又は請求項５記載のレーダ装置。
【請求項８】
前記レーダ装置は、ＦＭＣＷ（Frequency Modulated Continuous Wave）方式又はＳＦＷ（Stepped Frequency Waveforms）方式にて前記レーダ波の送受信を行うＦＭＣＷレーダ装置又はＳＦＷレーダ装置であること
を特徴とする請求項１〜請求項７のいずれか一項記載のレーダ装置。
【請求項９】
前記送受信手段は、ＦＭＣＷ（Frequency Modulated Continuous Wave）方式又はＳＦＷ（Stepped Frequency Waveforms）方式にて、前記レーダ波を送信アンテナから発射し、当該レーダ波の反射波をアレーアンテナで受信する手段であり、
前記レーダ装置は、
前記アレーアンテナを構成する各アンテナ素子からの受信信号に前記レーダ波の送信信号を混合して、前記アンテナ素子毎のビート信号を生成するビート信号生成手段と、
前記ビート信号生成手段により生成された前記アンテナ素子毎のビート信号に基づき、前記レーダ波の反射元までの距離及び前記反射元の方位を検出する信号解析手段と、
備え、
前記検出手段は、前記ビート信号生成手段により生成された前記アンテナ素子毎のビート信号に基づき、前記距離及び方位の異なる反射元毎に、当該反射元から到来した前記反射波の信号強度を検出し、
前記変化量算出手段は、前記反射元毎に、当該反射元から到来した前記反射波の前記検出手段により検出された信号強度Ｐと、当該信号強度Ｐに対応する前記反射波の受信時点から前記観測周期Ｔｓｗ遡った時点で当該反射元から到来し前記送受信手段により受信された前記反射波の前記検出手段により検出された信号強度Ｐｂと、に基づき、前記信号強度Ｐｂに対する前記信号強度Ｐの変化量Ｙを算出し、
前記判別手段は、前記反射元毎に、当該反射元の種類を、前記変化量算出手段で算出された前記変化量Ｙと、前記確率分布記憶手段が記憶する前記種類毎の確率分布とに基づき判別すること
を特徴とする請求項１〜請求項８のいずれか一項記載のレーダ装置。
【請求項１】
レーダ波を発射し、当該レーダ波の反射波を受信する送受信手段と、
前記送受信手段が受信した前記反射波の受信信号に基づき、前記反射波の信号強度を検出する検出手段と、
前記検出手段により検出された前記反射波の信号強度Ｐと、前記信号強度Ｐに対応する前記反射波の受信時点から所定の観測周期Ｔｓｗ遡った時点で前記送受信手段により受信された前記反射波の前記検出手段により検出された信号強度Ｐｂと、に基づき、前記信号強度Ｐｂに対する前記信号強度Ｐの変化量Ｙを算出する変化量算出手段と、
判別対象とする反射元の種類毎に、前記レーダ波が当該種類の反射元で反射して到来する前記反射波から算出される前記変化量Ｙの確率分布を記憶する確率分布記憶手段と、
前記確率分布記憶手段が記憶する前記種類毎の確率分布に基づき、前記変化量算出手段で算出された前記変化量Ｙに対応する前記レーダ波の反射元の種類を判別する判別手段と、
を備え、前記レーダ波の送信波長λ及び前記観測周期Ｔｓｗは、前記判別対象とする各種類の前記反射元にて発生する揺らぎの特徴が前記変化量Ｙに反映される関係で定められていること
を特徴とするレーダ装置。
【請求項２】
前記確率分布記憶手段は、前記種類毎の確率分布として少なくとも、前記反射元が金属物であるときの前記確率分布、及び前記反射元が人であるときの前記確率分布を記憶し、
前記判別手段は、前記確率分布に基づき、前記反射元の種類として、少なくとも前記金属物及び人の夫々を判別可能な構成にされていること
を特徴とする請求項１記載のレーダ装置。
【請求項３】
前記判別手段は、前記判別対象とする反射元の種類毎に、前記確率分布記憶手段が記憶する当該種類の前記確率分布に基づき、前記変化量算出手段で算出された前記変化量Ｙの発生確率であって当該変化量Ｙに対応する前記レーダ波の反射元が仮に当該種類の反射元であるときの発生確率を導出し、前記導出した発生確率が最大である前記反射元の種類を、前記変化量算出手段で算出された前記変化量Ｙに対応する反射元の種類として判別すること
を特徴とする請求項１又は請求項２記載のレーダ装置。
【請求項４】
車載型のレーダ装置であって、
自車両周辺の道路環境を推定する環境推定手段
を備え、
前記判別手段は、
前記判別対象とする反射元の種類毎に、前記確率分布記憶手段が記憶する当該種類の前記確率分布と、前記環境推定手段により推定された前記道路環境と、に基づき、前記変化量算出手段で算出された前記変化量Ｙに対応する前記レーダ波の反射元が当該種類の反射元である確度を算出し、
前記反射元の種類毎の前記確度に基づき、前記変化量算出手段で算出された前記変化量Ｙに対応する前記レーダ波の反射元の種類を判別すること
を特徴とする請求項１又は請求項２記載のレーダ装置。
【請求項５】
前記判別手段は、前記判別対象とする反射元の種類毎に、前記確率分布記憶手段が記憶する当該種類の前記確率分布に基づき、前記変化量算出手段で算出された前記変化量Ｙの発生確率であって前記レーダ波の反射元が当該種類の反射元であると仮定した場合での発生確率を導出し、当該導出した発生確率に、前記環境推定手段により推定された前記道路環境に対応した補正を加えることで、前記変化量算出手段で算出された前記変化量Ｙに対応する前記レーダ波の反射元が当該種類の反射元である確度を算出すること
を特徴とする請求項４記載のレーダ装置。
【請求項６】
前記環境推定手段は、各区域の道路環境を表すデータを記憶する道路環境データベースから、自車両の現在位置に対応する区域の道路環境を表すデータを取得して、この取得データが表す道路環境を、自車両周辺の道路環境であると推定すること
を特徴とする請求項４又は請求項５記載のレーダ装置。
【請求項７】
前記環境推定手段は、自車両の走行速度に基づき、前記道路環境を推定すること
を特徴とする請求項４又は請求項５記載のレーダ装置。
【請求項８】
前記レーダ装置は、ＦＭＣＷ（Frequency Modulated Continuous Wave）方式又はＳＦＷ（Stepped Frequency Waveforms）方式にて前記レーダ波の送受信を行うＦＭＣＷレーダ装置又はＳＦＷレーダ装置であること
を特徴とする請求項１〜請求項７のいずれか一項記載のレーダ装置。
【請求項９】
前記送受信手段は、ＦＭＣＷ（Frequency Modulated Continuous Wave）方式又はＳＦＷ（Stepped Frequency Waveforms）方式にて、前記レーダ波を送信アンテナから発射し、当該レーダ波の反射波をアレーアンテナで受信する手段であり、
前記レーダ装置は、
前記アレーアンテナを構成する各アンテナ素子からの受信信号に前記レーダ波の送信信号を混合して、前記アンテナ素子毎のビート信号を生成するビート信号生成手段と、
前記ビート信号生成手段により生成された前記アンテナ素子毎のビート信号に基づき、前記レーダ波の反射元までの距離及び前記反射元の方位を検出する信号解析手段と、
備え、
前記検出手段は、前記ビート信号生成手段により生成された前記アンテナ素子毎のビート信号に基づき、前記距離及び方位の異なる反射元毎に、当該反射元から到来した前記反射波の信号強度を検出し、
前記変化量算出手段は、前記反射元毎に、当該反射元から到来した前記反射波の前記検出手段により検出された信号強度Ｐと、当該信号強度Ｐに対応する前記反射波の受信時点から前記観測周期Ｔｓｗ遡った時点で当該反射元から到来し前記送受信手段により受信された前記反射波の前記検出手段により検出された信号強度Ｐｂと、に基づき、前記信号強度Ｐｂに対する前記信号強度Ｐの変化量Ｙを算出し、
前記判別手段は、前記反射元毎に、当該反射元の種類を、前記変化量算出手段で算出された前記変化量Ｙと、前記確率分布記憶手段が記憶する前記種類毎の確率分布とに基づき判別すること
を特徴とする請求項１〜請求項８のいずれか一項記載のレーダ装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【公開番号】特開２０１２−１１２６５３（Ｐ２０１２−１１２６５３Ａ）
【公開日】平成２４年６月１４日（２０１２．６．１４）
【国際特許分類】
【出願番号】特願２０１０−２５９０７０（Ｐ２０１０−２５９０７０）
【出願日】平成２２年１１月１９日（２０１０．１１．１９）
【出願人】（０００００４２６０）株式会社デンソー (27,639)
【出願人】（５０２３２４０６６）株式会社デンソーアイティーラボラトリ (332)
【Ｆターム（参考）】
【公開日】平成２４年６月１４日（２０１２．６．１４）
【国際特許分類】
【出願日】平成２２年１１月１９日（２０１０．１１．１９）
【出願人】（０００００４２６０）株式会社デンソー (27,639)
【出願人】（５０２３２４０６６）株式会社デンソーアイティーラボラトリ (332)
【Ｆターム（参考）】
[ Back to top ]