(57)【要約】【課題】異なる方位から到来する複数の到来波の中から出力ピークが低い到来波を検出する。【解決手段】処理対象信号Ｘｎ（ｔ）に対してビームフォーミング処理が行われ、その処理結果である角度スペクトラム信号Ｐｎ（θ）に基づいて到来波データＰｎ，θｎ，φｎが導出される。これらの到来波データＰｎ，θｎ，φｎに基づいて算出された到来波信号成分Ｙｎ（ｔ）で処理対象信号Ｘｎ（ｔ）が補正され、補正された処理対象信号に基づいて到来波データが導出される。【選択図】図３
2014174093.tif 000002
【特許請求の範囲】
【請求項１】
入力される受信データを処理対象信号として保持可能な処理対象信号保持部と、
前記処理対象信号保持部で保持されている処理対象信号に対してビームフォーマ法によるビームフォーミング処理を行うビームフォーミング部と、
前記ビームフォーミング部から出力される角度スペクトラム信号に基づいて、該角度スペクトラム信号における最大ピーク値の強度であるピーク強度、該最大ピーク値を有する到来波が到来する方位である到来方位、及び該到来方位に基づく到来波の位相推定値、のそれぞれを到来波データとして導出する導出部と、
前記導出部から出力される前記到来波データに基づく到来波信号成分を算出する到来波信号成分算出部と、
前記処理対象信号保持部で保持されている処理対象信号に含まれる前記到来波信号成分を低減するように、該処理対象信号を前記到来波信号成分で補正する補正部と、
前記処理対象信号保持部で保持されている処理対象信号を、前記補正部で補正された処理対象信号で置換し、置換後の該処理対象信号を新たな処理対象信号として前記処理対象信号保持部に保持させる置換部と、
を備えている、到来波推定装置。
【請求項２】
請求項１に記載の到来波推定装置であって、
前記補正部は、前記処理対象信号保持部で保持されている処理対象信号から前記到来波信号成分を減算することにより該処理対象信号を補正する、到来波推定装置。
【請求項３】
請求項１又は請求項２に記載の到来波推定装置であって、
前記導出部で導出された前記到来波データを保持するとともに、保持している前記到来波データを表示部へ出力可能な到来波データ保持部を更に備えている、到来波推定装置。
【請求項４】
請求項３に記載の到来波推定装置であって、
前記到来波データ保持部は、前記導出部で前記到来波データが導出された導出回数が所定回数未満の場合、保持している前記到来波データのうち最新の到来波データであるピーク強度、到来方位、及び位相推定値を前記到来波信号成分算出部へ出力する一方、前記導出回数が所定回数に達した場合、保持している前記到来波データを前記表示部へ出力する、到来波推定装置。
【請求項５】
請求項３に記載の到来波推定装置であって、
前記到来波データ保持部は、保持している前記ピーク強度のうち強度が最も低いピーク強度が所定値を超えている場合、該ピーク強度と、該ピーク強度に対応する前記到来方位及び前記位相推定値とを、前記到来波信号成分算出部へ出力する一方、保持している前記ピーク強度のうち強度が最も低いピーク強度が前記所定値以下になった場合、保持している前記到来波データを前記表示部へ出力する、到来波推定装置。
【請求項６】
請求項１から請求項５のいずれか１項に記載の到来波推定装置と、
複数のアンテナ素子を有し、電波を受信するアレイアンテナと、
前記アレイアンテナから出力されたエコー信号をデジタル信号に変換して前記到来波推定装置に出力する送受信装置と
を備えている、レーダ装置。
【請求項７】
請求項１から請求項５のいずれか１項に記載の到来波推定装置と、
複数のセンサを有し、超音波を送受信するアレイセンサと、
前記アレイセンサから出力されたエコー信号をデジタル信号に変換して前記到来波推定装置に出力する送受信装置と
を備えている、ソナー装置。
【請求項８】
入力される受信データを処理対象信号として保持するステップと、
保持している処理対象信号に対してビームフォーマ法によるビームフォーミング処理を行うステップと、
前記ビームフォーミング処理を行うステップで生成された角度スペクトラム信号に基づいて、該角度スペクトラム信号における最大ピーク値の強度であるピーク強度、該最大ピーク値を有する到来波が到来する方位である到来方位、及び該到来方位に基づく到来波の位相推定値、のそれぞれを到来波データとして導出するステップと、
前記到来波データに基づく到来波信号成分を算出するステップと、
保持している処理対象信号に含まれる前記到来波信号成分を低減するように、保持している処理対象信号を前記到来波信号成分で補正するステップと、
保持している処理対象信号を、前記補正するステップで補正された処理対象信号で置換し、置換後の該処理対象信号を新たな処理対象信号として保持するステップと
を含む、到来波推定方法。
【請求項９】
入力される受信データを処理対象信号として保持するステップと、
保持している処理対象信号に対してビームフォーマ法によるビームフォーミング処理を行うステップと、
前記ビームフォーミング処理を行うステップで生成された角度スペクトラム信号に基づいて、該角度スペクトラム信号における最大ピーク値の強度であるピーク強度、該最大ピーク値を有する到来波が到来する方位である到来方位、及び該到来方位に基づく到来波の位相推定値、のそれぞれを到来波データとして導出するステップと、
前記到来波データに基づく到来波信号成分を算出するステップと、
保持している処理対象信号に含まれる前記到来波信号成分を低減するように、保持している処理対象信号を前記到来波信号成分で補正するステップと、
保持している処理対象信号を、前記補正するステップで補正された処理対象信号で置換し、置換後の該処理対象信号を新たな処理対象信号として保持するステップと
をコンピュータに実行させるための到来波推定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、到来波の到来方向を推定する到来波推定装置、到来波推定装置を備えるレーダ装置、到来波推定装置を備えるソナー装置、到来波の到来方向を推定する到来波推定方法、及び到来波の到来方向を推定する到来波推定プログラムに関する。
【背景技術】
【０００２】
従来より、到来波の到来方法を推定する方法として、例えば非特許文献１に示すようなビームフォーマ法が知られている。ビームフォーマ法では、方位毎にその方位にメインビームが向くような所定のビームパターンを持つ受信ビームを形成することで、方位毎のエコーの強さを推定することができる。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】菊間信良著、「アレーアンテナによる適応信号処理」、科学技術出版、２００４年、ｐ．１７８−１８１
【発明の概要】
【発明が解決しようとする課題】
【０００４】
上述のようなビームフォーマ法を用いて、異なる方向から到来する複数の到来波の到来方向を推定する場合を考える。この場合、メインローブが、ピーク強度の低い到来波の方位に向いているときに、サイドローブが、ピーク強度の高い到来波の方向に向くことが考えらえる。そうなると、ピーク強度の低い到来波の出力ピークがサイドローブに埋もれてしまうおそれがある。
【０００５】
本発明は、上記課題を解決するためのものであり、その目的は、異なる方位から到来する複数の到来波の中から出力ピークが低い到来波を検出することである。
【課題を解決するための手段】
【０００６】
（１）上記課題を解決するため、本発明のある局面に係る到来波推定装置は、入力される受信データを処理対象信号として保持可能な処理対象信号保持部と、前記処理対象信号保持部で保持されている処理対象信号に対してビームフォーマ法によるビームフォーミング処理を行うビームフォーミング部と、前記ビームフォーミング部から出力される角度スペクトラム信号に基づいて、該角度スペクトラム信号における最大ピーク値の強度であるピーク強度、該最大ピーク値を有する到来波が到来する方位である到来方位、及び該到来方位に基づく到来波の位相推定値、のそれぞれを到来波データとして導出する導出部と、前記導出部から出力される前記到来波データに基づく到来波信号成分を算出する到来波信号成分算出部と、前記処理対象信号保持部で保持されている処理対象信号に含まれる前記到来波信号成分を低減するように、該処理対象信号を前記到来波信号成分で補正する補正部と、前記処理対象信号保持部で保持されている処理対象信号を、前記補正部で補正された処理対象信号で置換し、置換後の該処理対象信号を新たな処理対象信号として前記処理対象信号保持部に保持させる置換部と、を備えている。
【０００７】
（２）好ましくは、前記補正部は、前記処理対象信号保持部で保持されている処理対象信号から前記到来波信号成分を減算することにより該処理対象信号を補正する。
【０００８】
（３）好ましくは、前記到来波推定装置は、前記導出部で導出された前記到来波データを保持するとともに、保持している前記到来波データを表示部へ出力可能な到来波データ保持部を更に備えている。
【０００９】
（４）更に好ましくは、前記到来波データ保持部は、前記導出部で前記到来波データが導出された導出回数が所定回数未満の場合、保持している前記到来波データのうち最新の到来波データであるピーク強度、到来方位、及び位相推定値を、前記到来波信号成分算出部へ出力する一方、前記導出回数が所定回数に達した場合、保持している前記到来波データを前記表示部へ出力する。
【００１０】
（５）更に好ましくは、前記到来波データ保持部は、保持している前記ピーク強度のうち強度が最も低いピーク強度が所定値を超えている場合、該ピーク強度と、該ピーク強度に対応する前記到来方位及び前記位相推定値とを、前記到来波信号成分算出部へ出力する一方、保持している前記ピーク強度のうち強度が最も低いピーク強度が前記所定値以下になった場合、保持している前記到来波データを前記表示部へ出力する。
【００１１】
（６）また、上記課題を解決するため、本発明のある局面に係るレーダ装置は、上述したいずれかの到来波推定装置と、複数のアンテナ素子を有し、電波を送受信するアレイアンテナと、前記アレイアンテナから出力されたエコー信号をデジタル信号に変換して前記到来波推定装置に出力する送受信装置とを備えている。
【００１２】
（７）また、上記課題を解決するため、本発明のある局面に係るソナー装置は、上述したいずれかの到来波推定装置と、複数のセンサを有し、超音波を送受信するアレイセンサと、前記アレイセンサから出力されたエコー信号をデジタル信号に変換して前記到来波推定装置に出力する送受信装置とを備えている。
【００１３】
（８）また、上記課題を解決するため、本発明のある局面に係る到来波推定方法は、入力される受信データを処理対象信号として保持するステップと、保持している処理対象信号に対してビームフォーマ法によるビームフォーミング処理を行うステップと、前記ビームフォーミング処理を行うステップで生成された角度スペクトラム信号に基づいて、該角度スペクトラム信号における最大ピーク値の強度であるピーク強度、該最大ピーク値を有する到来波が到来する方位である到来方位、及び該到来方位に基づく到来波の位相推定値、のそれぞれを到来波データとして導出するステップと、前記到来波データに基づく到来波信号成分を算出するステップと、保持している処理対象信号に含まれる前記到来波信号成分を低減するように、保持している処理対象信号を前記到来波信号成分で補正するステップと、保持している処理対象信号を、前記補正するステップで補正された処理対象信号で置換し、置換後の該処理対象信号を新たな処理対象信号として保持するステップとを含む。
【００１４】
（９）また、上記課題を解決するため、本発明のある局面に係る到来波推定プログラムは、コンピュータに、入力される受信データを処理対象信号として保持するステップと、保持している処理対象信号に対してビームフォーマ法によるビームフォーミング処理を行うステップと、前記ビームフォーミング処理を行うステップで生成された角度スペクトラム信号に基づいて、該角度スペクトラム信号における最大ピーク値の強度であるピーク強度、該最大ピーク値を有する到来波が到来する方位である到来方位、及び該到来方位に基づく到来波の位相推定値、のそれぞれを到来波データとして導出するステップと、前記到来波データに基づく到来波信号成分を算出するステップと、保持している処理対象信号に含まれる前記到来波信号成分を低減するように、保持している処理対象信号を前記到来波信号成分で補正するステップと、保持している処理対象信号を、前記補正するステップで補正された処理対象信号で置換し、置換後の該処理対象信号を新たな処理対象信号として保持するステップとを実行させるためのプログラムである。
【発明の効果】
【００１５】
本発明によれば、異なる方位から到来する複数の到来波の中から出力ピークが低い到来波を検出できる。
【図面の簡単な説明】
【００１６】
【図１】本発明の第１実施形態に係るレーダ装置の構成を示すブロック図である。
【図２】図１に示すレーダ装置のアレイアンテナの詳細を示す概念図である。
【図３】図１に示すレーダ装置の到来波推定装置の構成を示すブロック図である。
【図４】図１に示すレーダ装置の動作を示すフローチャートである。
【図５】物標の存在状態の一例を示す図である。
【図６】図５に示す物標についての受信データに対してビームフォーミング処理が行われて算出された角度スペクトラム信号Ｐ１（θ）の一例を示す図である。
【図７】減算処理が１回行われた後に算出された角度スペクトラム信号Ｐ２（θ）の一例を示す図である。
【図８】減算処理が２回行われた後に算出された角度スペクトラム信号Ｐ３（θ）の一例を示す図である。
【図９】表示部の出力結果の一例を示す概略図である。
【図１０】変形例に係るレーダ装置の構成を示すブロック図である。
【図１１】図１０に示すレーダ装置の到来波推定装置の構成を示すブロック図である。
【図１２】図１０に示すレーダ装置の動作を示すフローチャートである。
【図１３】変形例に係るソナー装置の構成を示すブロック図である。
【発明を実施するための形態】
【００１７】
［第１実施形態］
本発明に係るレーダ装置の第１実施形態について図面を参照しつつ説明する。図１は、第１実施形態に係るレーダ装置１の構成を示すブロック図である。本実施形態のレーダ装置１は、例えば、漁船などの船舶に装備され、海上の物標（他船等）の位置を把握するために用いられる。
【００１８】
［レーダ装置の概要］
図１に示すように、本実施形態に係るレーダ装置１は、アンテナ２と、送受信装置５と、到来波推定装置１０と、表示装置９とを備えている。なお、到来波推定装置１０は、送受信装置５内に設けられていてもよい。
【００１９】
アンテナ２は、送信用のアンテナ３と、受信用のアレイアンテナ４とを有している。アレイアンテナ４は、例えば、Ｋ個のアンテナ素子４１〜４Ｋが等間隔且つ直線状に配置されている。
【００２０】
図２に、受信ビームを形成するアレイアンテナ４の概念図を示す。図２においては、所望波が角度θでアンテナ素子４１〜４Ｋに入射している状態が示されている。時刻ｔにおけるアンテナ素子４１〜４Ｋの複素受信信号をまとめて入力ベクトルＸ１（ｔ）で表すことができ、Ｘ１（ｔ）＝［ｘ１１（ｔ）、ｘ１２（ｔ）、・・・、ｘ１Ｋ（ｔ）］Ｔである。なお、Ｔは転置を表している。
【００２１】
送受信装置５は、送信部６と、受信部７と、局部発振器８とを有している。送信部６は、マイクロ波を発振する電子素子として、例えば、Ｄ／Ａコンバータ、周波数変換器、及び電力増幅器を有している。送信部６は、レーダ送信信号を生成して送信用アンテナ３に送出する。送信部６からレーダ送信信号が出力されるタイミングは、送信部６から到来波推定装置１０に出力される。これにより、到来波推定装置１０は、送信用アンテナ３から出力された電磁波の反射波が受信用のアレイアンテナ４で受信されるまでの時間を検知し、レーダ受信信号に係る電磁波の伝搬距離を算出することができる。
【００２２】
受信部７は、アレイアンテナ４からレーダ受信信号を取り込み、増幅した後に、中間周波数に変換してデジタル信号に変換し、到来波推定装置１０へ出力する。受信部７は、アレイアンテナ４の各アンテナ素子４１〜４Ｋのレーダ受信信号を処理するために、Ｋ個のアンテナ素子４１〜４Ｋのそれぞれに対応するＫ個の受信機７１〜７Ｋを有している。各受信機７１〜７Ｋは、ＬＮＡ（ローノイズアンプ）７ａ、ミキサ７ｂ、ＡＭＰ（リニアアンプ）７ｃ、及びＡ／Ｄコンバータ７ｄを有している。
【００２３】
ＬＮＡ７ａで増幅されたレーダ受信信号は、ミキサ７ｂにおいて局部発振器８から出力されるローカル信号とミキシングされ、中間周波数に変換される。ＡＭＰ７ｃは、中間周波数に変換された受信信号を増幅する。Ａ／Ｄコンバータ７ｄは、ＡＭＰ７ｃが出力した信号をＩＱ検波等により複素デジタル信号に変換し、この複素デジタル信号に変換後の受信データを到来波推定装置１０に対して出力する。
【００２４】
到来波推定装置１０は、送受信装置５から出力される受信データを処理して到来波のピーク強度、到来方向、及び位相を推定する処理を行い、その処理結果を表示装置９に出力する。到来波推定装置１０の構成及び動作については、後述する。
【００２５】
表示装置９は、表示部９ａ及び操作部９ｂを有している。表示部９ａは、到来波推定装置１０から出力された到来波に関するデータを表示するように構成されている。操作部９ｂには、種々の入力キー等の入力手段が設けられており、電磁波の送受信、信号処理、又は映像表示に必要な種々の設定又は種々のパラメータなどを入力できるように構成されている。
【００２６】
［到来波推定装置の構成］
図３は、本実施形態に係る到来波推定装置１０の構成を示すブロック図である。到来波推定装置１０は、処理対象信号生成部１１と、ビームフォーミング部１５と、導出部２０と、到来波データ保持部２５と、到来波信号成分算出部３０とを備えている。この到来波推定装置１０は、例えば図示しないＣＰＵ、ＦＰＧＡ、及びメモリなどのデバイスで構成される。例えば、ＣＰＵがメモリからプログラムを読み出して実行することにより、処理対象信号生成部１１と、ビームフォーミング部１５と、導出部２０と、到来波データ保持部２５と、到来波信号成分算出部３０との機能を構成することができる。
【００２７】
処理対象信号生成部１１は、到来波に関するデータを導出するための処理対象となる信号を生成する。処理対象信号生成部１１は、処理対象信号保持部１２と、減算部１３と、置換部１４とを有している。
【００２８】
処理対象信号保持部１２は、到来波推定装置１０での処理対象となる処理対象信号Ｘｎ（ｔ）＝［ｘｎ１（ｔ）、ｘｎ２（ｔ）、・・・、ｘｎＫ（ｔ）］Ｔ（但し、ｎ＝１，２，…，Ｎ）を保持する。具体的には、処理対象信号保持部１２は、送受信装置５から出力される受信データＸ１（ｔ）、又は、減算部１３による減算処理後の信号であるＸ２（ｔ），Ｘ３（ｔ），…，ＸＮ（ｔ）のいずれか１つを、処理対象信号Ｘｎ（ｔ）として保持する。なお、Ｘ２（ｔ）は、受信データＸ１（ｔ）について減算処理が１回行われた処理対象信号であり、Ｘ３（ｔ）は、Ｘ１（ｔ）について減算処理が２回行われた処理対象信号であり、ＸＮ（ｔ）は、Ｘ１（ｔ）について減算処理がＮ−１回行われた処理対象信号を示す。処理対象信号保持部１２は、保持している処理対象信号を、ビームフォーミング部１５へ出力する。
【００２９】
減算部１３は、処理対象信号保持部１２で保持されている処理対象信号Ｘｎ（ｔ）を、到来波信号成分算出部３０で算出された到来波信号成分Ｙｎ（ｔ）（詳しくは後述する）で補正する補正部として設けられている。本実施形態では、減算部１３は、処理対象信号保持部１２で保持されている処理対象信号Ｘｎ（ｔ）を到来波信号成分Ｙｎ（ｔ）で減算する減算処理を行うことにより、処理対象信号を補正する。
【００３０】
置換部１４は、処理対象信号保持部１２で保持されている処理対象信号Ｘｎ（ｔ）を、減算部１３で減算処理された処理対象信号Ｘｎ＋１（ｔ）で置換する。これにより、この置換後の処理対象信号Ｘｎ＋１（ｔ）は、新たな処理対象信号Ｘｎ＋１（ｔ）として、処理対象信号保持部１２に保持される。これにより、処理対象信号保持部１２には、その段階で減算処理が行われた回数が最も多い処理対象信号が保持される。
【００３１】
ビームフォーミング部１５は、相関行列算出部１６と、角度スペクトラム算出部１７とを有している。ビームフォーミング部１５は、処理対象信号生成部１１から順次、出力される処理対象信号Ｘｎ（ｔ）（Ｘ１（ｔ），Ｘ２（ｔ），…，ＸＮ（ｔ））に対して、ビームフォーマ法によるビームフォーミング処理を行う。
【００３２】
具体的には、相関行列算出部１６は、処理対象信号生成部１１から順次、出力される処理対象信号Ｘｎ（ｔ）に基づき、時刻ｔにおける各アンテナ素子４１〜４Ｋ間の受信信号の相関値を示す相関行列Ｒｎ（ｔ）（Ｒ１（ｔ），Ｒ２（ｔ），…，ＲＮ（ｔ））を算出する。相関行列Ｒｎ（ｔ）は、次の（１）式で与えられる。
【００３３】
［数１］
Ｒｎ（ｔ）＝Ｅ［Ｘｎ（ｔ）ＸｎＨ（ｔ）］    ・・・（１）
【００３４】
ここで、Ｅ［・］は、期待値（アンサンブル平均）を求める操作を表し、Ｈは複素共役転置を表す。アンサンブル平均は、時間平均で代用される。つまり、時刻ｔにおける相関行列Ｒｎ（ｔ）を求めるにあたり、その時刻ｔの前後のサンプリング時刻における入力ベクトルも用いる。一様時間平均の場合、相関行列Ｒｎ（ｔ）は次の式（２）で与えられる。
【００３５】
【数２】
2014174093.tif 000003
ここで、時間平均の範囲を［ｔ−Ｔｍｉｎ，ｔ＋Ｔｍａｘ］とした。
【００３６】
相関行列Ｒｎ（ｔ）の時間平均は、一様ではなく適当な重みを付けた時間平均でもよい。また、時間平均は行わなくてもよい。なお、本実施形態では、相関行列Ｒｎ（ｔ）の時間平均は行わない。
【００３７】
角度スペクトラム算出部１７は、相関行列算出部１６で算出された相関行列Ｒｎ（ｔ）に基づいて、角度スペクトラム信号Ｐｎ（θ）を算出する。具体的には、角度スペクトラム算出部１７は、式（３）に基づいて算出される出力電力を、式（４）で示す式で正規化することにより、角度スペクトラム信号Ｐｎ（θ）を算出する。角度スペクトラム算出部１７は、算出した角度スペクトラム信号Ｐｎ（θ）（Ｐ１（θ），Ｐ２（θ），…，ＰＮ（θ））を順次、導出部２０へ出力する。
【００３８】
【数３】
2014174093.tif 000004
【００３９】
【数４】
2014174093.tif 000005
【００４０】
ここで、式（３）及び式（４）におけるｗ（θ）は、Ｋ個のアンテナ素子で得られる複素信号ベクトルＸｎ（ｔ）＝［ｘ１（ｔ），…，ｘＫ（ｔ）］Ｔ（Ｔは転置の意味）に乗ずる複素ウェイトであり、式（５）で表される。
【数５】
2014174093.tif 000006
但し、λは搬送波（マイクロ波）の波長、ｊは虚数単位、ｄｋ（ｋ＝１，２，・・・，Ｋ）は基準点から第ｋ番目の超音波振動子２１ｋの位置までの距離を表す。
【００４１】
導出部２０は、検出部２１と、位相算出部２２とを有している。導出部２０は、ピーク強度Ｐｎ、到来方位θｎ、及び位相推定値φｎのそれぞれを、到来波データとして導出する。導出部２０は、角度スペクトラム算出部１７から順次、出力される角度スペクトラム信号Ｐｎ（θ）に基づき、各角度スペクトラム信号Ｐｎ（θ）に対応する到来波データを導出する。
【００４２】
具体的には、検出部２１は、角度スペクトラム算出部１７で算出された角度スペクトラム信号Ｐｎ（θ）においてθを変化させて、Ｐｎ（θ）のピーク値（極大値）のうち最も信号レベルの値が大きい最大ピーク値を検出する。そして、その最大ピーク値をピーク強度Ｐｎとして検出し、そのときの方位を到来方位θｎとして検出する。
【００４３】
位相算出部２２は、検出部２１で検出された到来方位θｎから、該到来方位θｎから到来する到来波の位相推定値φｎを、式（６）に基づいて算出する。
【００４４】
［数６］
φｎ＝arg［ｗ（θｎ）HXｎ（ｔ）］    ・・・（６）
【００４５】
すなわち、導出部２０は、各角度スペクトラム信号Ｐｎ（θ）（Ｐ１（θ），Ｐ２（θ），…，ＰＮ（θ））に含まれる最大ピーク値を有する到来波に関する到来波データＰｎ，θｎ，φｎ（Ｐ１，θ１，φ１、Ｐ２，θ２，φ２、…、ＰＮ，θＮ，φＮ）を導出する。導出部２０は、導出した到来波データＰｎ，θｎ，φｎを順次、到来波データ保持部２５に出力する。
【００４６】
到来波データ保持部２５は、規定回数記憶部２６と、出力部２７とを有する。到来波データ保持部２５は、導出部２０から順次、出力される到来波データＰｎ，θｎ，φｎを蓄積する。
【００４７】
規定回数記憶部２６には、ユーザによって規定される規定回数Ｎが記憶される。ユーザは、例えば、表示装置９の操作部９ｂを操作することにより、上記規定回数Ｎを入力することができる。
【００４８】
出力部２７は、到来波データを、導出部２０で到来波データＰｎ，θｎ，φｎが導出された導出回数に応じて、到来波信号成分算出部３０及び表示装置９のいずれか一方へ出力する。具体的には、出力部２７は、上記導出回数が上記規定回数Ｎ未満の場合、到来波データ保持部２５で保持されている到来波データのうち最新の到来波データであるピーク強度Ｐｎ、到来方位θｎ、及び位相推定値φｎを、到来波信号成分算出部３０へ出力する。一方、出力部２７は、上記導出回数が規定回数Ｎに達した場合、到来波データ保持部２５で保持される全ての到来波データＰ１，θ１，φ１、Ｐ２，θ２，φ２、…、ＰＮ，θＮ，φＮを表示装置９へ出力する。
【００４９】
なお、導出回数とは、その段階において到来波データＰｎ，θｎ，φｎが導出された回数である。例えば、処理対象信号Ｘ１（ｔ）に基づいて到来波データＰ１，θ１，φ１が導出され、その後、処理対象信号Ｘ２（ｔ）に基づいて到来波データＰ２，θ２，φ２が導出された段階では、導出回数は２回となる。
【００５０】
到来波信号成分算出部３０は、到来波データ保持部２５から順次、出力される到来波データＰｎ，θｎ，φｎから、到来波信号成分Ｙｎ（ｔ）（Ｙ１（ｔ），Ｙ２（ｔ），…，ＹＮ（ｔ））を算出する。具体的には、到来波信号成分算出部３０は、式（７）に基づいて到来波信号成分Ｙｎ（ｔ）を算出する。到来波信号成分算出部３０は、算出した到来波信号成分を順次、処理対象信号生成部１１へ出力する。
【００５１】
［数７］
Ｙｎ（ｔ）＝Ａｎ・ｗ（θｎ）Ｔ    ・・・（７）
【００５２】
ここで、式（７）におけるＡｎは、ピーク強度Ｐｎと位相推定値φｎとで表される複素振幅であり、式（８）で与えられる。
【００５３】
【数８】
2014174093.tif 000007
【００５４】
［レーダ装置の動作］
図４は、本実施形態に係るレーダ装置１の動作を示すフローチャートである。以下では、本実施形態に係るレーダ装置１の動作について、具体例を用いて説明する。
【００５５】
図５は、物標の存在状態の一例を示す図である。本実施形態では、図５に示すような物標に起因する到来波が到来する状況について説明する。なお、図５において、０度の方向から到来する０ｄＢの到来波をＡ、２０度の方向から到来する−１０ｄＢの到来波をＢ、−２０度の方向から到来する−２０ｄＢの到来波をＣ、とする。
【００５６】
まず、ステップＳ１では、処理対象信号Ｘｎ（ｔ）に乗ずる複素ウェイトｗ（θ）が生成される。この複素ウェイトｗ（θ）は、モードベクトルとも呼ばれ、アレイアンテナ４の形状等の条件により求めることができる。一旦モードベクトルが生成されると、このモードベクトルによって一連の信号処理が行われる。なお、アレイアンテナ４についての形状等の条件は機器設計時に仕様として与えられるものであるので、モードベクトルはあらかじめ算出できる。よって、あらかじめ算出しておいたモードベクトルの値を数値データとしてメモリ等に保持しておき、到来波推定装置１０における信号処理実行時に適宜読みだすようにしてもよい。
【００５７】
次に、ステップＳ２では、送信用のアンテナ３から送信ビームが自船の周囲に送信され、所定時間経過毎に、入力ベクトルＸ１（ｔ１）、Ｘ１（ｔ２）、Ｘ１（ｔ３）、・・・、Ｘ１（ｔｍ）が生成される。そして、受信部７は、これらの入力ベクトルＸ１（ｔ）を処理して複素デジタル信号に変換し、受信データとして処理対象信号生成部１１に出力する。この受信データは、処理対象信号Ｘ１（ｔ）として、処理対象信号保持部１２に保持される（ステップＳ３）。処理対象信号生成部１１は、この処理対象信号Ｘ１（ｔ）をビームフォーミング部１５に出力する。
【００５８】
次に、ステップＳ４では、相関行列算出部１６が、上述した式（１）及び式（２）に基づいて相関行列Ｒ１（ｔ）を算出する。そして、角度スペクトラム算出部１７が角度スペクトラム信号Ｐ１（θ）を算出する。具体的には、角度スペクトラム算出部１７は、−９０度から９０度までの各方位の出力電力（例えば、Ｐ（−９０°），Ｐ（−８９°），…，Ｐ（９０°））を、式（３）及び式（４）に基づいて算出する。ビームフォーミング部１５は、このように算出された角度スペクトラム信号Ｐ１（θ）を、導出部２０へ出力する。
【００５９】
図６に、角度スペクトラム算出部１７によって算出された角度スペクトラム信号Ｐ１（θ）の一例を示す。本実施形態では、図６に示すように、信号レベルが低い到来波（到来波Ｂ及び到来波Ｃ）はサイドローブに埋もれた状態となる。
【００６０】
次に、ステップＳ５では、導出部２０が、角度スペクトラム信号Ｐ１（θ）に基づいて、該角度スペクトラム信号Ｐ１（θ）に含まれる最大到来波（到来波Ａ）に関するデータである到来波データを導出する。この到来波データには、角度スペクトラム信号Ｐ１（θ）における最大ピーク値の強度であるピーク強度Ｐ１、該最大ピーク値を有する到来波が到来する方位である到来方位θ１、及び該到来方位に基づく到来波の位相推定値φ１が含まれる。
【００６１】
具体的には、まず、検出部２１が、角度スペクトラム信号Ｐ１（θ）の最大ピーク値の到来方位θ１、及び該到来方位θ１におけるピーク強度Ｐ１を検出する。次に、位相算出部２２が、式（６）に基づいて位相推定値φ１を算出する。導出部２０は、これらの到来方位θ１、ピーク強度Ｐ１、及び位相推定値φ１を、到来波データ保持部２５に出力する。到来波データ保持部２５は、これらの到来波データＰ１，θ１，φ１を保持する。この到来波データＰ１，θ１，φ１は、アレイアンテナ４で受信した受信データ中に含まれる、ピーク強度が最も高い到来波Ａに関するデータである。
【００６２】
次に、ステップＳ６では、到来波データが導出された導出回数が、ユーザによって予め設定される規定回数Ｎ以上であるか否かが判定される。本実施形態では、ユーザが、到来波の数が３つであると予め予測できており、上記規定回数Ｎを３回に設定しているとする。そして、本実施形態のここまでの段階では、到来波データの導出回数は１回である。よって、出力部２７は、導出回数が規定回数未満であるため（ステップＳ６のＮｏ）、ピーク強度Ｐ１、到来方位θ１、及び位相推定値φ１を、到来波信号成分算出部３０へ出力する。
【００６３】
次に、ステップＳ７では、到来波信号成分算出部３０が、Ｐ１、θ１、及びφ１に基づいて到来波信号成分Ｙ１（ｔ）を算出する。具体的には、到来波信号成分算出部３０は、式（７）及び式（８）に基づいて、到来波信号成分Ｙ１（ｔ）を算出し、この到来波信号成分Ｙ１（ｔ）を処理対象信号生成部１１へ出力する。この到来波信号成分Ｙ１（ｔ）は、本実施形態におけるアレイアンテナ４が実際に受信している３つの到来波のうち到来波Ａのみを受信していると仮定した場合の、アレイアンテナ４の受信データを推定したものである。
【００６４】
次に、ステップＳ８では、処理対象信号生成部１１の減算部１３が、処理対象信号保持部１２が保持している処理対象信号Ｘ１（ｔ）から、到来波信号成分算出部３０から出力された到来波信号成分Ｙ１（ｔ）を減算する。このように減算処理された処理対象信号Ｘ１（ｔ）−Ｙ１（ｔ）は、置換部１４によって、処理対象信号保持部１２が保持している処理対象信号に置換される。これにより、処理対象信号保持部１２は、減算処理された処理対象信号Ｘ１（ｔ）−Ｙ１（ｔ）を、新たな処理対象信号Ｘ２（ｔ）として保持する（ステップＳ９）。この処理対象信号Ｘ２（ｔ）は、本実施形態におけるアレイアンテナ４が実際に受信している３つの到来波のうち到来波Ａを除く２つの到来波（到来波Ｂ及び到来波Ｃ）を受信していると仮定した場合の、アレイアンテナ４の受信データを推定したものである。
【００６５】
上記処理対象信号Ｘ２（ｔ）について、再びステップＳ４からステップＳ６の処理が施される。具体的には、ステップＳ４で、ビームフォーミング部１５が処理対象信号Ｘ２（ｔ）に対してビームフォーミング処理を行って、相関行列Ｒ２（ｔ）及び角度スペクトラム信号Ｐ２（θ）を算出する。この角度スペクトラム信号Ｐ２（θ）は、アレイアンテナ４が到来波Ｂ及び到来波Ｃを受信していると仮定した場合の角度スペクトラム信号を推定したものである。
【００６６】
図７に、角度スペクトラム算出部１７によって算出された角度スペクトラム信号Ｐ２（θ）の一例を示す。図７に示すように、到来波Ａに起因する到来波信号成分Ｙ１（ｔ）が減算された角度スペクトラム信号Ｐ２（θ）では、２番目に信号レベルが高い到来波Ｂを確認できる。しかし、３つの到来波の中で最も信号レベルが低い到来波Ｃについては、当該到来波Ｃに起因するピークなのか、サイドローブなのか、判別がつきにくくなっている。
【００６７】
次に、ステップＳ５では、導出部２０が、角度スペクトラム信号Ｐ２（θ）に基づいて、該角度スペクトラム信号Ｐ２（θ）に含まれる最大到来波（到来波Ｂ）に関するデータである到来波データＰ２，θ２，φ２を導出する。この到来波データＰ２，θ２，φ２は、到来波データ保持部２５に出力され、到来波データ保持部２５は、既に保持している到来波データＰ１，θ１，φ１とともに、上記Ｐ２，θ２，φ２を保持する。
【００６８】
次に、ステップＳ６で、到来波データの導出回数が規定回数Ｎ以上であるか否かが判定される。本実施形態のここまでの段階では、到来波データの導出回数は２回である。よって出力部２７は、保持している到来波データのうち最新の到来波データであるピーク強度、到来方位、及び位相推定値を、到来波信号成分算出部３０へ出力する。具体的には、出力部２７は、到来波Ｂに関する到来波データＰ２，θ２，φ２を、到来波信号成分算出部３０へ出力する。
【００６９】
上記Ｐ２，θ２，φ２について、再びステップＳ７からステップＳ９の処理が施される。具体的には、ステップＳ７で、到来波信号成分算出部３０が、Ｐ２，θ２，φ２に基づいて到来波信号成分Ｙ２（ｔ）を算出する。この到来波信号成分Ｙ２（ｔ）は、アレイアンテナ４が到来波Ｂのみを受信していると仮定した場合、アレイアンテナ４の受信データを推定したものである。
【００７０】
次に、ステップＳ８で、減算部１３が減算処理を行って新たな処理対象信号Ｘ３（ｔ）（＝Ｘ２（ｔ）−Ｙ２（ｔ））を生成する。処理対象信号保持部１２は、この新たな処理対象信号Ｘ３（ｔ）を保持する（ステップＳ９）。この処理対象信号Ｘ３（ｔ）は、アレイアンテナ４が到来波Ｃのみを受信していると仮定した場合の、アレイアンテナ４の受信データを推定したものである。
【００７１】
上記処理対象信号Ｘ３（ｔ）について、再びステップＳ４からステップＳ６の処理が施される。具体的には、ステップＳ４で、ビームフォーミング部１５が処理対象信号Ｘ３（ｔ）に対してビームフォーミング処理を行って、相関行列Ｒ３（ｔ）及び角度スペクトラム信号Ｐ３（θ）を算出する。この角度スペクトラム信号Ｐ３（θ）は、アレイアンテナ４が到来波Ｃを受信していると仮定した場合の角度スペクトラム信号を推定したものである。
【００７２】
図８に、角度スペクトラム算出部１７によって算出された角度スペクトラム信号Ｐ３（θ）の一例を示す。図８に示すように、到来波Ａに起因する到来波信号成分Ｙ１（ｔ）及び到来波Ｂに起因する到来波信号成分Ｙ２（ｔ）が減算された角度スペクトラム信号Ｐ３（θ）では、最も信号レベルが低い到来波Ｃを確認できる。
【００７３】
次に、ステップＳ５では、導出部２０が、角度スペクトラム信号Ｐ３（θ）に基づいて、該角度スペクトラム信号Ｐ３（θ）に含まれる最大到来波（到来波Ｃ）に関するデータである到来波データＰ３，θ３，φ３を導出する。この到来波データＰ３，θ３，φ３は、到来波データ保持部２５に出力され、到来波データ保持部２５は、既に保持している到来波データＰ１，θ１，φ１及びＰ２，θ２，φ２とともに、上記Ｐ３，θ３，φ３を保持する。
【００７４】
そして、ステップＳ６で、到来波データの導出回数が規定回数Ｎ以上であるか否かが判定される。本実施形態のこの段階では、到来波データの導出回数は３回となっている。よって、出力部２７は、保持している全ての到来波データＰ１〜Ｐ３，θ１〜θ３及びφ１〜φ３を表示装置９へ出力する（ステップＳ１０）。
【００７５】
図９に、表示装置９の表示部９ａの一例を示す。図９に示すように、表示部９ａには、３つの到来波に関するデータが、非常に幅の狭いパルス状の波形として表示される。このように表示される波形により、３つの到来波Ａ，Ｂ及びＣのそれぞれに関する到来波データを確認することができる。なお、位相推定値φ１〜φ３については、例えば図９における各波形付近に表示される。また、到来波データの表示方法は、図９に示すものに限らず、例えば、表示部９ａに、到来波データの値を表示してもよい。
【００７６】
また、上述した説明においては、到来波データの導出回数を３回と設定しているが、この限りでなく、ユーザが想定する到来波の数に応じて、ユーザが適宜、設定することができる。具体的には、到来波の数が１０程度であるとユーザが想定した場合、その数よりもやや多い１２程度の数を、規定回数Ｎとして設定してもよい。この場合、本実施形態に係る到来波推定装置１０は、ステップＳ６において導出回数が規定回数であるＮに到達するまでステップＳ４からステップＳ９の処理を繰り返す。そして、導出回数が１２回に到達すると、到来波データ保持部２５は、それまでに導出された全ての到来波データＰ１，θ１，φ１、Ｐ２，θ２，φ２、…、Ｐ１２，θ１２，φ１２、を表示装置９の表示部９ａへ出力する。そして、表示部９ａに、これらの到来波データに関する情報が表示される。これにより、ユーザは、強度の異なる複数の到来波のそれぞれに関する情報を知ることができる。
【００７７】
また、ユーザは、到来波の数を想定できない場合であっても、Ｎとして比較的大きい数を設定すればよい。この場合、ユーザは、ピーク強度がある程度のレベルに達しているものを到来波データと推定し、ピーク強度が極端に低いものをノイズと推定することにより、物標に起因する到来波に関する到来波データを推定できる。
【００７８】
［実施形態の効果］
以上のように、本実施形態に係る到来波推定装置１０では、処理対象信号Ｘｎ（ｔ）に含まれる前記到来波信号成分Ｙｎ（ｔ）を低減するように、処理対象信号Ｘｎ（ｔ）を、該処理対象信号Ｘｎ（ｔ）に含まれる最大到来波の寄与成分である到来波信号成分Ｙｎ（ｔ）で補正している。こうすると、処理対象信号Ｘｎ（ｔ）において信号レベルが最も高い到来波の影響を低減できる。これにより、信号レベルの低い到来波が、信号レベルの高い到来波のサイドローブに埋もれてしまうのを抑制できるため、信号レベルの低い到来波のピークを確認できる。
【００７９】
従って、本実施形態に係る到来波推定装置１０では、異なる方位から到来する複数の到来波の中から出力ピークが低い到来波を検出することができる。
【００８０】
また、到来波推定装置１０では、処理対象信号Ｘｎ（ｔ）から到来波信号成分Ｙｎ（ｔ）を減算することにより、処理対象信号Ｘｎ（ｔ）に含まれる最大到来波に起因する成分を除去している。これにより、処理対象信号Ｘｎ（ｔ）に含まれる最大到来波を適切に取り除くことができる。
【００８１】
また、到来波推定装置１０は、到来波データを保持するとともに該到来波データを表示部へ出力可能な到来波データ保持部２５を備えている。これにより、ユーザは、受信データＸ１（ｔ）に含まれる複数の到来波に関する情報である到来波データを確認できる。
【００８２】
また、到来波推定装置１０によれば、到来波データが導出された導出回数が規定回数Ｎ以上となるまで、到来波データの導出が繰り返される。これにより、ユーザは、到来波の数を推測した上で規定回数Ｎを適切な回数に設定することで、所望の到来波に関する到来波データを確認できる。また、到来波推定装置１０によれば、ユーザが到来波の数を推定せずとも、Ｎとして比較的大きい数を設定すればよい。この場合、ユーザは、ピーク強度がある程度のレベルに達しているものを到来波データと推定し、ピーク強度が極端に低いものをノイズと推定することができる。こうすると、到来波の数を推測する手間が省けるため、複数の物標に起因する到来波データを容易に推測できる。
【００８３】
また、本実施形態に係るレーダ装置１によれば、異なる方位から到来する複数の到来波の中から出力ピークが低い到来波を検出することができる。
【００８４】
以上、本発明の実施形態について説明したが、本発明はこれらに限定されるものではなく、本発明の趣旨を逸脱しない限りにおいて種々の変更が可能である。
【００８５】
（１）図１０は、変形例に係るレーダ装置１ａの構成を示すブロック図であり、図１１は、変形例に係るレーダ装置１ａの到来波推定装置４０の構成を示すブロック図である。図１０及び図１１に示すように、本変形例に係る到来波推定装置４０は、上記実施形態に係る到来波推定装置１０と比べると、到来波データ保持部の構成が異なる。以下、上記実施形態と異なる点について主に説明し、上記実施形態と同様の構成の説明については、図面において同一の符号を付すことで又は同一の符号を引用して説明することで、省略する。
【００８６】
本変形例に係る到来波推定装置４０の到来波データ保持部３５は、閾値記憶部３６と、出力部３７とを有している。閾値記憶部３６には、ユーザによって規定される閾値Ｐｔｈｒが記憶される。ユーザは、例えば、表示装置９の操作部９ｂを操作することにより、上記閾値Ｐｔｈｒを入力できる。
【００８７】
出力部３７は、到来波データ保持部３５が保持しているピーク強度の値に応じて、到来波データを到来波信号成分算出部３０及び表示装置９のいずれか一方へ出力する。具体的には、出力部２７は、保持しているピーク強度Ｐ１、Ｐ２、…、Ｐｑのうち強度が最も低いピーク強度Ｐｑが閾値Ｐｔｈｒを超えている場合、Ｐｑ、θｑ及びφｑを、到来波信号成分算出部３０へ出力する。一方、保持しているピーク強度Ｐ１、Ｐ２、…、Ｐｑのうち強度が最も低いピーク強度ＰＱが閾値Ｐｔｈｒ以下になった場合、保持している全ての到来波データＰ１，θ１，φ１、Ｐ２，θ２，φ２、…、ＰＱ，θＱ，φＱを表示部へ出力する。
【００８８】
図１２は、本変形例に係るレーダ装置の動作を示すフローチャートである。以下では、上記実施形態の場合と同様、図５に示すような物標に起因する到来波が到来する状況について説明する。
【００８９】
本変形例に係るレーダ装置は、ステップＳ１からステップＳ５までは、上記実施形態に係るレーダ装置１と同じように動作する。具体的には、レーダ装置は、処理対象信号Ｘｑ（ｔ）に乗ずるモードベクトルを生成し（ステップＳ１）、送信用のアンテナ３から送信ビームを自船の周囲に送信する（ステップＳ２）。次に、アレイアンテナ４で受信された受信信号が受信部７で処理されて複素デジタル信号に変換され、受信データとして処理対象信号生成部１１に出力される。到来波推定装置４０では、この受信データが処理対象信号Ｘ１（ｔ）として保持される（ステップＳ３）。そして、この処理対象信号Ｘ１（ｔ）に基づいて角度スペクトラム信号Ｐ１（θ）が算出され（ステップＳ４）、この角度スペクトラム信号Ｐ１（θ）に基づいて到来波データＰ１，θ１，φ１が導出される（ステップＳ５）。この到来波データＰ１，θ１，φ１は、到来波データ保持部３５に出力され、到来波データ保持部３５は、到来波データＰ１，θ１，φ１を保持する。
【００９０】
次に、本変形例に係る到来波推定装置４０では、保持しているピーク強度のうち強度が最も低いピーク強度が閾値Ｐｔｈｒ以下であるか否かが判定される。本変形例では、ユーザが閾値Ｐｔｈｒを−１８ｄＢと設定しているとする。本変形例のこの時点では、保持しているピーク強度はＰ１のみであり、Ｐ１は０ｄＢである。よって、Ｐ１はＰｔｈｒを上回っているため（ステップＳ１１のＮｏ）、出力部２７は、ピーク強度Ｐ１、到来方位θ１、及び位相推定値φ１を、到来波信号成分算出部３０へ出力する。
【００９１】
次に、到来波推定装置４０は、ステップＳ７からステップＳ９、ステップＳ４、及びステップＳ５について、上記実施形態の場合と同様に動作する。具体的には、到来波信号成分Ｙ１（ｔ）が算出され（ステップＳ７）、Ｘ１（ｔ）について減算処理が行われ（ステップＳ８）、新たな処理対象信号Ｘ２（ｔ）が生成されて処理対象信号生成部１１に保持される（ステップＳ９）。そして、この新たな処理対象信号Ｘ２（ｔ）に基づいて角度スペクトラム信号Ｐ２（θ）が算出され（ステップＳ４）、この角度スペクトラム信号Ｐ２（θ）に基づいて到来波データＰ２，θ２，φ２が導出される（ステップＳ５）。この到来波データＰ２，θ２，φ２は、到来波データ保持部３５に出力され、到来波データ保持部３５は、既に保持しているＰ１，θ１，φ１とともに、上記Ｐ２，θ２，φ２を保持する。
【００９２】
次に、再びステップＳ１１で、保持しているピーク強度のうち強度が最も低いピーク強度が閾値Ｐｔｈｒ以下であるか否かが判定される。本変形例のこの時点では、保持しているピーク強度のうち強度が最も低いピーク強度はＰ２であり、Ｐ２は−１０ｄＢである。よって、Ｐ２はＰｔｈｒを上回っているため（ステップＳ１１のＮｏ）、出力部２７は、ピーク強度Ｐ２、到来方位θ２、及び位相推定値φ２を、到来波信号成分算出部３０へ出力する。
【００９３】
次に、到来波推定装置４０は、ステップＳ７からステップＳ９、ステップＳ４、及びステップＳ５について、上記実施形態の場合と同様に動作する。具体的には、到来波信号成分Ｙ２（ｔ）が算出され（ステップＳ７）、Ｘ２（ｔ）について減算処理が行われ（ステップＳ８）、新たな処理対象信号Ｘ３（ｔ）が生成されて処理対象信号生成部１１に保持される（ステップＳ９）。そして、この新たな処理対象信号Ｘ３（ｔ）に基づいて角度スペクトラム信号Ｐ３（θ）が算出され（ステップＳ４）、この角度スペクトラム信号Ｐ３（θ）に基づいて到来波データＰ３，θ３，φ３が導出される（ステップＳ５）。この到来波データＰ３，θ３，φ３は、到来波データ保持部３５に出力され、到来波データ保持部３５は、既に保持しているＰ１，θ１，φ１、Ｐ２，θ２，φ２とともに、上記Ｐ３，θ３，φ３を保持する。
【００９４】
次に、再びステップＳ１１で、保持しているピーク強度のうち強度が最も低いピーク強度が閾値Ｐｔｈｒ以下であるか否かが判定される。本変形例のこの時点では、保持しているピーク強度のうち強度が最も低いピーク強度はＰ３であり、Ｐ３は−２０ｄＢである。よって、Ｐ３はＰｔｈｒ以下になるため（ステップＳ１１のＹｅｓ）、出力部２７は、保持している全ての到来波データＰ１〜Ｐ３，θ１〜θ３及びφ１〜φ３を表示装置９へ出力する（ステップＳ１０）。これにより、ユーザは、強度の異なる複数の到来波のそれぞれに関する情報を知ることができる。
【００９５】
以上のように、本変形例に係る到来波推定装置４０でも、上記実施形態の場合と同様、異なる方位から到来する複数の到来波の中から出力ピークが低い到来波を検出することができる。
【００９６】
また、到来波推定装置４０によれば、ピーク強度が閾値Ｐｔｈｒ以下になるまで、到来波データの導出が繰り返される。これにより、ユーザが閾値Ｐｔｈｒを適切なレベルに設定することで、ユーザは、所望の出力レベル以上の到来波データを確認できる。
【００９７】
（２）上記実施形態及び上記変形例では、到来波推定装置をレーダ装置に適用しているが、この限りでなく、ソナー装置に適用することもできる。
【００９８】
図１３は、変形例に係るソナー装置１ｂの  構成を示すブロック図である。以下、上記実施形態と異なる点について主に説明し、上記実施形態と同様の構成の説明については、図面において同一の符号を付すことで又は同一の符号を引用して説明することで、省略する。また、本変形例に係るソナー装置１ｂの到来波推定装置４０の構成及び動作についても、上記実施形態に係るレーダ装置１の構成及び動作と概ね同様であるため、その説明を省略する。
【００９９】
［ソナー装置の概要］
図１３に示すように、本実施形態に係るソナー装置１ｂは、トランスデューサ５０（アレイセンサ）と、送受信装置５１と、到来波推定装置１０と、表示装置９とを備えている。なお、到来波推定装置１０は、送受信装置５１内に設けられていてもよい。
【０１００】
トランスデューサ５０は、例えば、Ｍ個の超音波振動子が直線状に配置されているリニアアレイ型トランスデューサである。超音波振動子は、例えば、キャリア周波数の波長の２分の１の間隔で、等間隔に配置されている。各超音波振動子は、電気信号を超音波に変換して水中に超音波を送信するとともに、受信した超音波を電気信号（エコー信号）に変換する。トランスデューサ５０は、超音波の送受信を行うため、例えば、船底に装備されて海中に露出されている。
【０１０１】
送受信装置５１は、送受切替部５２と、送信部５３と、受信部５４とを備えている。送受切替部５２は、送信時には、送信部５３からトランスデューサ５０に送信信号が送られる接続に切り替える。また、送受切替部５２は、受信時には、トランスデューサ５０によって超音波から変換された電気信号がトランスデューサ５０から受信部５４に送られる接続に切り替える。
【０１０２】
送信部５３は、表示装置９において設定された条件に基づいて送信信号を生成する。また、送信部５３は、送受切替部５２を介して、生成した送信信号をトランスデューサ５０に対して出力する。この送信信号により、各超音波振動子が駆動され、トランスデューサ５０から各方位に超音波が放射される。
【０１０３】
受信部５４は、増幅部５４ａと、Ａ／Ｄ変換部５４ｂとを有している。増幅部５４ａは、Ｍ個の超音波振動子のそれぞれから送受切替部５２を介して送られてきた電気信号を増幅してＭ個の受信信号をＡ／Ｄ変換部５４ｂに対して出力する。
【０１０４】
Ａ／Ｄ変換部５４ｂは、増幅部５４ａから出力される受信信号をＩＱ検波等により複素デジタル信号に変換し、この複素デジタル信号に変換後の受信データを到来波推定装置１０に対して出力する。
【０１０５】
そして、本変形例に係る到来波推定装置１０は、上記実施形態の場合と同様の動作をすることで、複数の到来波のそれぞれに対応する到来波データＰ１，θ１，φ１、Ｐ２，θ２，φ２、…、Ｐｎ，θｎ，φｎを導出する。
【０１０６】
以上のように、本変形例に係るソナー装置１ｂでも、異なる方位から到来する複数の到来波の中から出力ピークが低い到来波を検出することができる。
【０１０７】
（３）また、到来波推定装置１０，４０は、上述したレーダ装置１，１ａ、及びソナー装置１ｂに限らず、超音波診断装置に適用することもできる。
【０１０８】
（４）また、上記実施形態及び変形例では、処理対象信号Ｘｎ（ｔ）から到来波信号成分Ｙｎ（ｔ）を減算することにより新たな処理対象信号を生成しているが、この限りでない。例えば、到来波信号成分Ｙｎ（ｔ）に所定の係数を乗じたものを、処理対象信号Ｘｎ（ｔ）から減算してもよい。更には、上記実施形態及び変形例では、処理対象信号Ｘｎ（ｔ）を到来波信号成分Ｙｎ（ｔ）で補正する補正部を、減算部１３で構成しているが、この限りでなく、補正部は、処理対象信号Ｘｎ（ｔ）に含まれる到来波信号成分Ｙｎ（ｔ）を低減できれば、どのような構成であってもよい。
【産業上の利用可能性】
【０１０９】
本発明は、アレイアンテナを有するレーダ装置、ソナー装置、又は超音波診断装置として広く適用することができるものである。
【符号の説明】
【０１１０】
１，１ａ        レーダ装置
１ｂ            ソナー装置
１２            処理対象信号保持部
１３            減算部（補正部）
１４            置換部
１５            ビームフォーミング部
２０            導出部
２５            到来波データ保持部
３０            到来波信号成分算出部
【図１】
2014174093.tif 000008
【図２】
2014174093.tif 000009
【図３】
2014174093.tif 000010
【図４】
2014174093.tif 000011
【図５】
2014174093.tif 000012
【図６】
2014174093.tif 000013
【図７】
2014174093.tif 000014
【図８】
2014174093.tif 000015
【図９】
2014174093.tif 000016
【図１０】
2014174093.tif 000017
【図１１】
2014174093.tif 000018
【図１２】
2014174093.tif 000019
【図１３】
2014174093.tif 000020