(57)【特許請求の範囲】
【請求項１】
複数の受信素子のそれぞれで受波された受信波から得られる受信信号に基づいて受信信号ベクトルを生成する信号ベクトル生成部と、
異なるタイミングで生成された複数の前記受信信号ベクトルの中から対象受信信号ベクトルを選択する対象信号ベクトル選択部と、
複数の前記対象受信信号ベクトルに基づいて、前記受信波が到来した方位を推定する物標方位推定部と、
を備えていることを特徴とする、信号処理装置。
【請求項２】
請求項１に記載の信号処理装置において、
前記対象信号ベクトル選択部は、前記受信信号ベクトルに含まれる複数の前記受信信号の強度のばらつき度合いに基づいて、前記対象受信信号ベクトルを選択することを特徴とする、信号処理装置。
【請求項３】
請求項２に記載の信号処理装置において、
前記ばらつき度合いは、標準偏差、最大値と最小値との差としての範囲、又は分散であることを特徴とする、信号処理装置。
【請求項４】
請求項２又は請求項３に記載の信号処理装置において、
前記対象信号ベクトル選択部は、前記ばらつき度合いに基づく値が閾値以下となる前記受信信号ベクトルを前記対象受信信号ベクトルとして選択することを特徴とする、信号処理装置。
【請求項５】
請求項１に記載の信号処理装置において、
前記対象信号ベクトル選択部は、各前記受信信号ベクトルの相関行列を生成する相関行列生成部を有し、該相関行列生成部で生成された複数の前記相関行列の中から対象相関行列を選択することにより、前記対象受信信号ベクトルを選択することを特徴とする、信号処理装置。
【請求項６】
請求項５に記載の信号処理装置において、
前記対象信号ベクトル選択部は、複数の前記相関行列のそれぞれの時間的変化に基づいて、前記対象相関行列を選択することを特徴とする、信号処理装置。
【請求項７】
請求項６に記載の信号処理装置において、
前記時間的変化は、各前記相関行列の微分値であることを特徴とする、信号処理装置。
【請求項８】
請求項７に記載の信号処理装置において、
前記対象信号ベクトル選択部は、前記微分値に基づく値が閾値以下となる前記相関行列を前記対象相関行列として選択することを特徴とする、信号処理装置。
【請求項９】
請求項４又は請求項８に記載の信号処理装置において、
前記閾値は、所定値として設定され、又は、ノイズレベルに応じて設定されることを特徴とする、信号処理装置。
【請求項１０】
請求項１から請求項９のいずれか１項に記載の信号処理装置において、
前記物標方位推定部は、ビームフォーミング法又は適応ビームフォーミング法を用いて前記受信波が到来した方位を推定することを特徴とする、信号処理装置。
【請求項１１】
請求項１から請求項１０のいずれか１項に記載の信号処理装置において、
前記複数の受信素子は、送信部から送信された送信波が物標で反射して帰来する反射波を、前記受信波として受波し、
前記信号ベクトル生成部は、前記反射波から得られる前記受信信号に基づいて前記受信信号ベクトルを生成することを特徴とする、信号処理装置。
【請求項１２】
請求項１から請求項１１のいずれか１項に記載の信号処理装置を備えた、水中探知装置。
【請求項１３】
請求項１から請求項１１のいずれか１項に記載の信号処理装置を備えた、レーダ装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の受信素子のそれぞれで受波された受信波から得られる受信信号を処理する信号処理装置、当該信号処理装置を備えた水中探知装置、及びレーダ装置に関する。
【背景技術】
【０００２】
従来より、複数の方位から到来する到来波を互いに判別するための様々な手法が知られている。例えば、非特許文献１には、相関行列に周波数平均を用いることにより、複数の方位から到来する到来波を互いに判別する手法が開示されている。また、非特許文献２には相関行列に空間平均を用いる手法が、非特許文献３には相関行列に時間平均を用いる手法が、それぞれ開示されている。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】H.Wang他、「Coherent Signal-Subspace Processing for the Detection and Estimation of Angles of Arrival of Multiple Wide-Band Sources」、IEEE TRANSACTIONS ON ACOUSTICS, SPEECH, AND SIGNAL PROCESSING、IEEE、1985年8月、vol.33、No.4、P823-P831
【非特許文献２】Harry L. Van Trees、「Optimum Array Processing Part IV of Detection, Estimation, and Modulation Theory」、WILEY-INTERSCIENCE、chapter6.12.4、P605-P607
【非特許文献３】D.R. Van Rheeden他、「A Temporal Smoothing Approach to Direction of Arrival Estimation of Coherent Signals in Fading Channels」、Wireless Communications and Networking Conference、1999年、vol.1、P286-P290
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところで、上述のように単に時間平均、空間平均、周波数平均を利用して物標の方位を推定した場合、十分な分解能を得ることができない。
【０００５】
本発明は、上記課題を解決するためのものであり、その目的は、分解能の劣化を抑制しつつ、複数の方位から到来する到来波を互いに判別することである。
【課題を解決するための手段】
【０００６】
（１）上記課題を解決するために、この発明のある局面に係る信号処理装置は、複数の受信素子のそれぞれで受波された受信波から得られる受信信号に基づいて受信信号ベクトルを生成する信号ベクトル生成部と、異なるタイミングで生成された複数の前記受信信号ベクトルの中から対象受信信号ベクトルを選択する対象信号ベクトル選択部と、複数の前記対象受信信号ベクトルに基づいて、前記受信波が到来した方位を推定する物標方位推定部と、を備えている。
【０００７】
（２）好ましくは、前記対象信号ベクトル選択部は、前記受信信号ベクトルに含まれる複数の前記受信信号の強度のばらつき度合いに基づいて、前記対象受信信号ベクトルを選択する。
【０００８】
（３）更に好ましくは、前記ばらつき度合いは、標準偏差、最大値と最小値との差としての範囲、又は分散である。
【０００９】
（４）好ましくは、前記対象信号ベクトル選択部は、前記ばらつき度合いに基づく値が閾値以下となる前記受信信号ベクトルを前記対象受信信号ベクトルとして選択する。
【００１０】
（５）好ましくは、前記対象信号ベクトル選択部は、各前記受信信号ベクトルの相関行列を生成する相関行列生成部を有し、該相関行列生成部で生成された複数の前記相関行列の中から対象相関行列を選択することにより、前記対象受信信号ベクトルを選択する。
【００１１】
（６）更に好ましくは、前記対象信号ベクトル選択部は、複数の前記相関行列のそれぞれの時間的変化に基づいて、前記対象相関行列を選択する。
【００１２】
（７）更に好ましくは、前記時間的変化は、各前記相関行列の微分値である。
【００１３】
（８）更に好ましくは、前記対象信号ベクトル選択部は、前記微分値に基づく値が閾値以下となる前記相関行列を前記対象相関行列として選択する。
【００１４】
（９）好ましくは、前記閾値は、所定値として設定され、又は、ノイズレベルに応じて設定される。
【００１５】
（１０）好ましくは、前記物標方位推定部は、ビームフォーミング法又は適応ビームフォーミング法を用いて前記受信波の方位を推定する。
【００１６】
（１１）好ましくは、前記複数の受信素子は、送信部から送信された送信波が物標で反射して帰来する反射波を、前記受信波として受波し、前記信号ベクトル生成部は、前記反射波から得られる前記受信信号に基づいて前記受信信号ベクトルを生成する。
【００１７】
（１２）上記課題を解決するために、この発明のある局面に係る水中探知装置は、上述したいずれかの信号処理装置を備えている。
【００１８】
（１３）上記課題を解決するために、この発明のある局面に係るレーダ装置は、上述したいずれかの信号処理装置を備えている。
【発明の効果】
【００１９】
本発明によれば、分解能の劣化を抑制しつつ、複数の方位から到来する到来波を互いに判別できる。
【図面の簡単な説明】
【００２０】
【図１】本発明の実施形態に係る水中探知装置の構成を示すブロック図である。
【図２】図１に示す信号処理部において到来波の方位推定が行われるまでに行われる処理を説明するための図である。
【図３】（Ａ）は、本発明の実施形態に係る水中探知装置によって算出されるエコー強度のシミュレーション結果であって、方位方向及び距離方向に対するエコー強度を示すグラフであり、（Ｂ）は、従来の水中探知装置によって算出されるエコー強度のシミュレーション結果であって、（Ａ）に対応させて示す図である。
【図４】各時刻での受信信号ベクトルごとに生成された相関行列を用いて算出された角度スペクトラムのグラフを説明するための図であって、（Ａ）は各超音波振動子によって受信される受信信号の振幅の時間的変化の一例を示す図、（Ｂ）から（Ｄ）は、各時刻のときの相関行列に基づいて算出された角度スペクトラムのグラフである。
【図５】変形例に係る水中探知装置の構成を示すブロック図である。
【図６】図５に示す信号処理部において到来波の方位推定が行われるまでに行われる処理を説明するための図である。
【図７】変形例に係る水中探知装置の構成を示すブロック図である。
【図８】図７に示す信号処理部において到来波の方位推定が行われるまでに行われる処理を説明するための図である。
【図９】本発明の実施形態に係るレーダ装置の構成を示すブロック図である。
【図１０】本発明の実施形態に係る信号処理装置の構成を示すブロック図である。
【発明を実施するための形態】
【００２１】
本発明の実施形態に係る信号処理部３０（信号処理装置）を備えた水中探知装置１について、図を参照して説明する。本発明の実施形態に係る水中探知装置１は、水中の物標（例えば魚群）を探知するために用いられる。水中探知装置１は、例えば漁船等の船舶に装備される。
【００２２】
［構成］
図１は、水中探知装置１の構成を示すブロック図である。水中探知装置１は、図１に示すように、トランスデューサ１０（送信部、受信部）と、送信機１１と、送受切替部１５と、受信機２０と、操作・表示装置４０とを備えている。
【００２３】
トランスデューサ１０は、超音波の送受信を行うため、例えば、船底に装備されて海中に露出されている。本実施形態におけるトランスデューサ１０は、送信素子アレイ及び受信素子アレイの双方として動作する送受信素子アレイとして設けられている。トランスデューサ１０は、複数の超音波振動子Ａｊ（送受信素子、受信素子）を有している。複数の超音波振動子Ａｊは、例えば一例として、互いに間隔をあけて直線状に配列されている。しかし、複数の超音波振動子Ａｊの配置については、これに限らず、例えば、２次元状又は３次元状に配置されていてもよい。
【００２４】
送受切替部１５は、送信時には、送信機１１からトランスデューサ１０に送信信号（パルス波）が送られる接続に切り替える。また、送受切替部１５は、受信時には、トランスデューサ１０によって超音波から変換された電気信号がトランスデューサ１０から受信機２０に送られる接続に切り替える。
【００２５】
送信機１１は、操作・表示装置４０において設定された条件に基づいてパルス波を生成する。また、送信機１１は、送受切替部１５を介して、該パルス波が所定時刻毎に各超音波振動子Ａｊから送信されるように、該パルス波をトランスデューサ１０に対して出力する。このパルス波により、各超音波振動子Ａｊが駆動され、トランスデューサ１０から所定時刻毎にパルス波（超音波）が放射される。
【００２６】
受信機２０は、増幅部２１と、Ａ／Ｄ変換部２２と、信号処理部３０とを有している。増幅部２１は、各超音波振動子Ａｊから送受切替部１５を介して送られてきた電気信号を増幅し、それらをＡ／Ｄ変換部２２に対して出力する。Ａ／Ｄ変換部２２は、増幅部２１から出力される受信信号をＩＱ検波等により複素デジタル信号に変換し、この複素デジタル信号に変換後の受信信号を信号処理部３０に出力する。
【００２７】
信号処理部３０は、Ａ／Ｄ変換部２２から出力される受信信号を処理して角度スペクトラムを生成するとともに、該角度スペクトラムに基づき到来波の到来方向および強度を推定する処理を行う。また、信号処理部３０は、角度スペクトラムに基づいて映像信号を生成する処理を行う。信号処理部３０は、生成した映像信号を操作・表示装置４０に出力する。信号処理部３０については、詳しくは後述する。
【００２８】
操作・表示装置４０は、表示部４１と、操作部４２とを有している。表示部４１は、映像信号に応じた映像を表示画面に表示する。操作部４２は、例えば種々の入力キーなどを有し、超音波の送受信や映像表示に必要な種々の設定や種々のパラメータなどを入力できるように構成されている。
【００２９】
［信号処理部］
信号処理部３０は、受信信号を処理して到来波の到来方向および強度を推定する処理を行う。信号処理部３０は、詳しくは後述するが、到来波の到来方向を推定する際に用いられる受信信号ベクトルのうち、所定の条件を満たす受信信号ベクトルを選択する。そして、信号処理部３０は、当該受信信号ベクトルを用いて到来波の到来方向を推定することにより、物標の方位を正確に導出することができる。以下、信号処理部３０について詳細に説明する。
【００３０】
信号処理部３０は、図１に示すように、信号ベクトル生成部３１と、相関行列生成部３２と、閾値設定部３３と、対象信号ベクトル選択部３４と、時間平均処理部３６と、到来波演算部３７（物標方位推定部）とを備えている。
【００３１】
図２は、信号処理部３０において到来波の方位推定が行われるまでに行われる処理を説明するための図である。図２では、トランスデューサ１０を構成する各超音波振動子Ａｊ（ｊ＝１，２，…）が各時刻ｔｉ（ｉ＝１，２，…）で信号を受信している状態を模式的に示している。
【００３２】
信号ベクトル生成部３１は、各超音波振動子Ａｊ（ｊ＝１，２，…）によって各時刻ｔｉ（ｉ＝１，２，…）に受信された受信信号ｘｊ（ｔｉ）に基づき、複数の超音波振動子によって時刻ｔｉ毎に得られる受信信号ベクトルＸｉ（＝[ｘ１（ｔｉ），ｘ２（ｔｉ），ｘ３（ｔｉ），…，ｘｊ（ｔｉ），…］Ｔ）を生成する。
【００３３】
相関行列生成部３２は、信号ベクトル生成部３１によって生成された受信信号ベクトルＸｉに基づき、時刻ｔｉ毎に相関行列Ｒｉ（＝［ＸｉＸｉＨ］）を生成する。
【００３４】
閾値設定部３３は、相関行列生成部３２で生成された相関行列Ｒｉに基づき、閾値ｄＲｍａｘを設定する。閾値ｄＲｍａｘは、到来波の到来方向を推定するための対象相関行列を選択するために用いられるものである。閾値ｄＲｍａｘは、各相関行列Ｒｉ（ｉ＝１，２，…）に対応して生成される。具体的には、閾値ｄＲｍａｘは、相関行列Ｒｉの微分行列のノルム||ｄＲｉ／ｄｔ||Ｆの標準偏差ＳＤｉに所定の係数αを乗算した値である（ｄＲｍａｘ＝α×ＳＤｉ（||ｄＲｉ／ｄｔ||Ｆ））。相関行列Ｒｉの微分行列の要素は相関行列Ｒｉの要素の微分となる。これにより、閾値設定部３３によって、ノイズの影響及び各エコーのパワーのばらつきが加味された適切な閾値が設定される。
【００３５】
対象信号ベクトル選択部３４は、対象相関行列選択部３５を有している。対象相関行列選択部３５は、相関行列生成部３２で生成された時刻ｔｉ毎の相関行列Ｒｉの中から、到来波の到来方向を推定するために用いられる対象相関行列を選択する。具体的には、対象相関行列選択部３５は、各相関行列Ｒｉのうち、該各相関行列Ｒｉの微分のノルムが閾値ｄＲｍａｘ以下の相関行列を、対象相関行列ＲＳＪＴｎ（ｎ＝１，２，…）として選択する。一方、各相関行列Ｒｉのうち、該各相関行列の微分のノルムが閾値ｄＲｍａｘを超える相関行列については、対象相関行列ＲＳＪＴｎから除外する。これにより、対象信号ベクトル選択部３４は、結果的に、複数の受信信号ベクトルの中から、到来方向の推定に用いられる対象受信信号ベクトルを選択していることになる。
【００３６】
時間平均処理部３６は、対象相関行列ＲＳＪＴｎ（ｎ＝１，２，…）を平均する処理（時間平均処理）を行う。これにより、時間平均処理部３６は、平均相関行列ＲＡＶＥを算出する。
【００３７】
到来波演算部３７は、時間平均処理部３６で算出された平均相関行列ＲＡＶＥに基づいて到来波の方向及び強度を算出する。到来波演算部３７は、平均相関行列ＲＡＶＥを用いてＭＵＳＩＣ法、Ｃａｐｏｎ法、サブスペース空間の分析を使用したビームフォーマ等によって、到来波の角度スペクトラムを導出する。到来波演算部３７は、このように算出した到来波の角度スペクトラムから映像信号を生成し、この映像信号を操作・表示装置４０に出力する。
【００３８】
図３（Ａ）は、本実施形態に係る水中探知装置１によって算出されるエコー強度のシミュレーション結果であって、方位方向及び距離方向に対するエコー強度を示すグラフである。また、図３（Ｂ）は、従来の水中探知装置によって算出されるエコー強度のシミュレーション結果であって、図３（Ａ）に対応させて示す図である。図３では、２０°方向及び−６０°方向から反射波が帰来している例を示している。すなわち、図３に示す例では、２０°方向及び−６０°方向に物標が存在している。また、図３では、エコーの強度を、ハッチングの濃度に対応させて示している。具体的には、エコー強度が強い領域ではハッチングの濃度が濃くなっており、エコー強度が弱い領域ではハッチングの高度が薄くなっているか、又は、ハッチングが省略されている。
【００３９】
図３（Ｂ）に示すように、従来の水中探知装置（対象相関行列の選択を行わず、生成した全ての相関行列に基づいてエコー強度を算出するように構成された水中探知装置）では、本実施形態に係る水中探知装置１と比べて、方位分解能が悪く（すなわち、物標に起因するエコー強度の方位方向の幅が広く）なっている。更に、従来の水中探知装置では、本来は存在しない物標の虚像（図３（Ｂ）の場合、７０°〜９０°付近に発生するエコー像）が生じる場合がある。これに対して、本実施形態に係る水中探知装置１によれば、受信信号に対して上述のような処理を行うことにより、図３（Ａ）に示すように、方位分解能を向上することができ、更に、虚像の発生を抑制することができる。従って、探知対象となる物標の方位をより正確に把握することができる。
【００４０】
［対象相関行列について］
上述のように、対象相関行列選択部３５では、各相関行列の微分のノルムが閾値ｄＲｍａｘ以下の相関行列を、対象相関行列ＲＳＪＴｎとして選択する。言い換えれば、対象相関行列選択部３５において選択されない相関行列は、上述したノルムが閾値を超える相関行列である。このような相関行列は、到来波の波面の立ち上がり部分（又は立ち下がり部分）における相関行列であり、例えば図２の場合では、ｔ２〜ｔ４の範囲内の各時刻における受信信号ベクトルに基づいて生成された相関行列Ｒ２〜Ｒ４である。このような相関行列が、到来方向を推定するために用いられる対象相関行列から除外される理由を、以下で説明する。
【００４１】
図４は、各時刻での受信信号ベクトルごとに生成された相関行列を用いて角度スペクトラムを算出した場合の角度スペクトラムのグラフを説明するための図であって、（Ａ）は各超音波振動子によって受信される受信信号の振幅の時間的変化の一例を示す図、（Ｂ）から（Ｄ）は、各時刻（ｔ２，ｔ４，ｔｉ）のときの相関行列に基づいて算出された角度スペクトラムのグラフである。
【００４２】
例えば、複数の超音波振動子Ａｊのうち１つの超音波振動子のみが物標からの反射波を受信している場合、その時刻ｔ２における相関行列によって算出される角度スペクトラムは、図４（Ｂ）に示すように、方位分解能が非常に低いグラフとなる。また、複数の超音波振動子のうちのいくつかが物標からの反射波を受信している場合、その時刻ｔ４における相関行列によって算出される角度スペクトラムは、図４（Ｃ）に示すようなグラフとなる。このグラフは、図４（Ｂ）のグラフに比べると方位分解能が高くなっているものの、方位分解能は十分とは言えない。これらに対して、複数の超音波振動子Ａｊの全てで反射波を受信している場合、その時刻ｔｉにおける相関行列によって算出される角度スペクトラムは、図４（Ｄ）に示すように、方位分解能が非常に高いグラフとなる。
【００４３】
すなわち、本実施形態の対象相関行列選択部３５は、方位分解能が低い角度スペクトラムを算出する原因となる相関行列（図２に示す例の場合、ｔ２〜ｔ４のときの相関行列）を対象相関行列から除外し、残りの相関行列（対象相関行列）に基づいて物標の方位を推定している。これにより、水中探知装置１の方位分解能を向上することができる。
【００４４】
［効果］
以上のように、本実施形態に係る水中探知装置１の信号処理部３０では、複数の受信信号ベクトルの中から、物標の方位を推定するために用いられる対象受信信号ベクトルが選択される。こうすると、物標の方位分解能が劣化する要因となる受信信号ベクトルを除外し、残りの受信信号ベクトルで物標の方位を推定することが可能になる。
【００４５】
従って、信号処理部３０では、分解能の劣化を抑制しつつ、複数の方位から到来する到来波を互いに判別できる。
【００４６】
また、信号処理部３０では、複数の相関行列の時間平均を行う際に生成される各時刻での相関行列から必要な相関行列が選択されるため、受信信号ベクトルの信号処理を行う一連の流れの中で対象相関行列を選択することができる。
【００４７】
また、信号処理部３０では、各相関行列の時間変化に基づいて、対象相関行列が選択される。これにより、分解能が劣化する要因となる相関行列を、対象相関行列から適切に除外することができる。
【００４８】
また、信号処理部３０では、各相関行列の微分値に基づいて、対象相関行列が選択される。これにより、各相関行列の時間変化を、容易に導出することができる。
【００４９】
また、信号処理部３０では、各相関行列の微分値に基づく値（具体的には、各相関行列の微分値のノルム）が閾値以下となる相関行列が、対象相関行列として選択される。言い換えると、上述した微分値に基づく値が閾値を超える相関行列（すなわち、方位分解能の劣化の要因となる相関行列）が対象相関行列から除外されるため、対象相関行列をより適切に選択することができる。
【００５０】
また、信号処理部３０では、上述した閾値がノイズレベルに応じて設定される。具体的には、信号処理部３０では、閾値ｄＲｍａｘは、相関行列Ｒｉの微分のノルムの標準偏差に基づいて設定される。これにより、ノイズレベルに応じた適切な閾値を設定することができる。
【００５１】
また、信号処理部３０では、適応ビームフォーミング法を用いて到来波の到来方向及びエコー強度の算出が行われるため、到来波の到来方向を正確に把握することができる。
【００５２】
また、水中探知装置１では、分解能の劣化を抑制しつつ、複数の方位から到来する到来波を互いに判別できる水中探知装置を提供できる。
【００５３】
［変形例］
以上、本発明の実施形態について説明したが、本発明はこれらに限定されるものではなく、本発明の趣旨を逸脱しない限りにおいて種々の変更が可能である。
【００５４】
（１）図５は、変形例に係る水中探知装置１ａの構成を示すブロック図である。本変形例に係る水中探知装置１ａは、上記実施形態に係る水中探知装置１の場合と異なり、相関行列を生成する前に、装置の方位分解能に悪影響を及ぼす受信信号ベクトルを予め取り除いた後、その残りの受信信号ベクトル（対象受信信号ベクトル）から相関行列（対象相関行列）を生成するように構成されている。
【００５５】
図６は、図５に示す信号処理部３０ａにおいて到来波の方位推定が行われるまでに行われる処理を説明するための図である。本変形例に係る水中探知装置１ａは、上記実施形態に係る水中探知装置１と比べて、信号処理部の構成が異なっている。具体的には、本変形例の信号処理部３０ａは、図５に示すように、信号ベクトル生成部３１、閾値設定部３８、対象信号ベクトル選択部３９、対象相関行列生成部４５、時間平均処理部３６、及び到来波演算部３７、を有している。これらのうち、信号ベクトル生成部３１、時間平均処理部３６、及び到来波演算部３７については、上記実施形態の場合と構成及び動作が同じであるため、その説明を省略する。
【００５６】
閾値設定部３８は、受信信号ベクトルＸｉ（＝[ｘ１（ｔｉ），ｘ２（ｔｉ），ｘ３（ｔｉ），…，ｘｊ（ｔｉ），…］Ｔ）に基づき、対象信号ベクトルを選択するために用いられる閾値ｄＸｍａｘを設定する。具体的には、閾値ｄＸｍａｘは、受信信号ベクトルＸｉに含まれる全ての要素（受信信号）の強度の標準偏差ＳＤｊの時間的な標準偏差ＳＤｉに所定の係数βを乗算した値である。本変形例では、閾値ｄＸｍａｘは、ｄＸｍａｘ＝β×ＳＤｉ（ＳＤｊ（ｘｊ（ｔｉ）））、と表すことができる。
【００５７】
対象信号ベクトル選択部３９は、時刻ｔｉ（ｉ＝１，２，…）毎の受信信号ベクトルＸｉの中から、到来波の到来方向を推定するために用いられる対象信号ベクトルを選択する。具体的には、対象信号ベクトル選択部３９は、各受信信号ベクトルＸｉのうち、該各受信信号ベクトルＸｉに含まれる全ての要素の強度の標準偏差が閾値ｄＸｍａｘ以下の受信信号ベクトルを、対象受信信号ベクトルＸＳＪＴｎ（ｎ＝１，２，…）として選択する。一方、各受信信号ベクトルＸｉのうち、該各受信信号ベクトルＸｉに含まれる全ての要素の強度の標準偏差が閾値ｄＸｍａｘを超える受信信号ベクトルについては、対象受信信号ベクトルＸＳＪＴｎから除外する。
【００５８】
対象相関行列生成部４５は、対象信号ベクトル選択部３９で選択された対象受信信号ベクトルＸＳＪＴｎのそれぞれにおいて、相関行列を生成する。この対象相関行列生成部４５で生成される相関行列は、到来波の到来方位を推定するために用いられる対象相関行列ＲＳＪＴｎである。
【００５９】
以上のように、本変形例に係る信号処理部３０ａは、上記実施形態に係る信号処理部３０の場合と比べて、対象相関行列ＲＳＪＴｎを生成する過程が異なるものの、上記実施形態の場合と同様の対象相関行列ＲＳＪＴｎを生成することができる。これにより、上記実施形態の場合と同様、分解能の劣化を抑制しつつ、複数の方位から到来する到来波を互いに判別できる。
【００６０】
また、信号処理部３０ａでは、各受信信号ベクトルに含まれる複数の受信信号の強度のばらつき度合いに基づいて、対象受信信号ベクトルが選択される。これにより、分解能が劣化する要因となる受信信号ベクトルを、対象受信信号ベクトルから適切に除外することができる。
【００６１】
また、信号処理部３０ａでは、各受信信号ベクトルに含まれる全ての要素の強度の標準偏差に基づいて、対象受信信号ベクトルが選択される。これにより、各受信信号ベクトルに含まれる各受信信号の時間的変化を、容易に導出することができる。
【００６２】
また、信号処理部３０ａでは、各受信信号ベクトルに含まれる全ての要素の強度の標準偏差に基づく値が閾値以下となる受信信号ベクトルが、対象受信信号ベクトルとして選択される。言い換えると、上述した標準偏差に基づく値が閾値を超える受信信号ベクトル（すなわち、方位分解能の劣化の要因となる受信信号ベクトル）が対象受信信号ベクトルから除外されるため、対象受信信号ベクトルをより適切に選択することができる。
【００６３】
なお、本変形例では、各受信信号ベクトルに含まれる全ての要素のばらつき度合いとして、標準偏差を用いたが、これに限らず、その他のばらつき度合い、例えば、範囲（最大値と最小値の差）、分散等を利用することもできる。
【００６４】
（２）上記実施形態では、対象相関行列ＲＳＪＴｎを選択するための閾値ｄＲｍａｘを、各相関行列Ｒｉの微分値に基づいて算出したが、これに限らず、各相関行列Ｒｉの全要素の値のばらつき度合い（例えば、標準偏差、範囲（最大値と最小値の差）、分散等）に基づいて設定してもよい。
【００６５】
図２及び図４を参照して説明すると、対象相関行列として選択したい相関行列は、例えば時刻ｔ１及び時刻ｔ５のそれぞれで生成される相関行列のように、受信信号ベクトルに含まれる各受信信号の値のばらつき度合いが小さい相関行列である。言い換えれば、対象相関行列から除外したい相関行列は、例えば時刻ｔ２〜ｔ４のそれぞれで生成される相関行列のように、受信信号ベクトルに含まれる各受信信号の値のばらつき度合いが大きい相関行列である。よって、本変形例のように、各受信信号の値のばらつき度合いが大きい相関行列を対象相関行列から除外することで、上記実施形態の場合と同様、適切に対象相関行列を選択することができる。
【００６６】
（３）上記実施形態では、複数の対象相関行列を時間平均して平均相関行列ＲＡＶＥを算出したが、これに限らず、空間平均を組み合わせて平均相関行列を算出してもよい。
【００６７】
（４）上記実施形態では、平均相関行列ＲＡＶＥを算出する際、各対象相関行列ＲＳＪＴｎに重み付けを行わずに平均処理を行ったが、これに限らず、各対象相関行列ＲＳＪＴｎに重み付けを行った後、平均処理を行ってもよい。
【００６８】
（５）図７は、変形例に係る水中探知装置１ｂの構成を示すブロック図である。上述した実施形態及び変形例における到来波の到来方向の推定手法の一例としては、トランスデューサ１０によって受信された信号に基づいて平均相関行列ＲＡＶＥを生成した後、この平均相関行列ＲＡＶＥの固有値及び固有ベクトルを算出し、これらの固有値及び固有ベクトルに基づいて到来波の到来方向を推定する手法が挙げられる。しかし、本変形例に係る水中探知装置１ｂでは、平均相関行列ＲＡＶＥを生成することなく、当該平均相関行列ＲＡＶＥの固有値及び固有ベクトルが算出される。具体的には、本変形例では、特異値分解を行うことにより、平均相関行列ＲＡＶＥの固有値及び固有ベクトルが算出される。
【００６９】
本変形例に係る水中探知装置１ｂの信号処理部３０ｂは、図７に示すように、信号ベクトル生成部３１、閾値設定部３８、対象信号ベクトル選択部３９、データ行列生成部４６、特異値分解部４７、及び到来波演算部３７ａ、を有している。これらのうち、信号ベクトル生成部３１、閾値設定部３８、及び対象信号ベクトル選択部３９については、図５の場合と同様であるため、その説明を省略する。
【００７０】
図８は、図７に示す信号処理部３０ｂにおいて到来波の方位推定が行われるまでに行われる処理を説明するための図である。データ行列生成部４６は、対象信号ベクトル選択部３９によって選択された対象受信信号ベクトルＸＳＪＴｎ（ｎ＝１，２，…，Ｎ）から、データ行列Ｘ（＝［ＸＳＪＴ１，ＸＳＪＴ２，…］）を生成する。
【００７１】
特異値分解部４７は、データ行列生成部４６で生成されたデータ行列Ｘを、次の式（１）のように分解（特異値分解）する。
【００７２】
［数１］
Ｘ＝ＵＳＶＨ    ・・・（１）
【００７３】
上式（１）において、Ｓは対角行列であり、その対角要素が特異値と呼ばれる。また、Ｕ及びＶは、列ベクトルが互いに直交するユニタリ行列であり、それぞれが特異ベクトルと呼ばれる。なお、式（１）におけるＵは左特異ベクトル、Ｖは右特異ベクトル、とも呼ばれる。
【００７４】
到来波演算部３７ａは、特異値分解部４７で算出された左特異ベクトルＵから、平均相関行列ＲＡＶＥの固有ベクトルを算出する。また、到来波演算部３７ａは、特異値分解部４７で算出された特異値から、平均相関行列ＲＡＶＥの固有値を算出する。
【００７５】
ここで、データ行列Ｘと平均相関行列ＲＡＶＥとの関係は、次の式（２）で与えられる。
【００７６】
【数２】
0006403327.tif 000002
【００７７】
式（２）において、ＸＸＨ＝Ｒ、とした場合、データ行列Ｘの特異値の２乗は、Ｒの固有値に等しい。また、データ行列Ｘの左特異ベクトルＵは、Ｒの固有ベクトルに等しい。到来波演算部３７ａは、これらの関係を用いて、平均相関行列ＲＡＶＥの固有値及び固有ベクトルを算出する。そして、到来波演算部３７ａは、上述のように算出された固有値及び固有ベクトルに基づいて、上記実施形態の場合と同様、到来波の到来方向、及び到来波のエコー強度の算出を行う。
【００７８】
以上のように、本変形例に係る水中探知装置１ｂのように特異値分解を利用しても、上記実施形態の場合と同様、分解能の劣化を抑制しつつ、複数の方位から到来する到来波を互いに判別できる。
【００７９】
（６）上記実施形態及び上記変形例では、適応ビームフォーミング法を用いて到来波の到来方向を推定する例を挙げて説明したが、これに限らず、ビームフォーミング法を用いて到来波の到来方向を推定してもよい。
【００８０】
（７）上記実施形態及び上記変形例では、探知装置の一例として水中探知装置を挙げて説明したが、これに限らず、レーダ装置に適用することもできる。
【００８１】
図９は、本発明の実施形態に係るレーダ装置１ｃの構成を示すブロック図である。レーダ装置１ｃは、海上の物標（例えば他船）を探知するために用いられる。レーダ装置１ｃは、例えば漁船等の船舶に装備される。
【００８２】
以下、レーダ装置１ｃについて、上記実施形態に係る水中探知装置１と異なる点について主に説明し、水中探知装置１と同様の構成の説明については、図面において同一の符号を付すことで又は同一の符号を引用して説明することで、省略する。
【００８３】
［構成］
図９に示すように、レーダ装置１ｃは、アンテナ１０ａ（送信部、受信部）と、送信機１１と、送受切替部１５と、受信機２０と、操作・表示装置４０とを備えている。
【００８４】
アンテナ１０ａは、電磁波の送受信を行うため、例えば船上に装備されている。本変形例におけるアンテナ１０ａは、送信素子アレイ及び受信素子アレイの双方として動作する。アンテナ１０ａは、送信素子及び受信素子の双方として動作するアンテナ素子Ｂｊ（送受信素子、受信素子）を複数有している。これらのアンテナ素子Ｂｊは、例えば、直線状に配列されている。しかし、複数のアンテナ素子Ｂｊの配置については、これに限らず、例えば、２次元状又は３次元状に配置されていてもよい。
【００８５】
そして、レーダ装置１ｃの送信機１１、送受切替部１５、受信機２０、操作・表示装置４０は、送信波及び受信波として電磁波を取り扱う点を除き、上記実施形態に係る水中探知装置１の場合と同様に動作する。
【００８６】
従って、レーダ装置１ｃによれば、上記実施形態に係る水中探知装置１の場合と同様、分解能の劣化を抑制しつつ、複数の方位から到来する到来波を互いに判別できる。
【００８７】
（８）図１０は、本発明の実施形態に係る信号処理装置３０ｄの構成を示すブロック図である。信号処理装置３０ｄは、信号ベクトル生成部３１と、対象信号ベクトル選択部３９と、物標方位推定部３７ｂと、を有している。
【００８８】
上記実施形態及び変形例では、相関行列の各要素の値、又は受信信号ベクトルの各受信信号の値に基づいて閾値を設定したが、これに限らない。具体的には、本実施形態では、閾値として、予め設定された所定値が用いられる。
【００８９】
また、上記実施形態及び変形例では、適応ビームフォーミング法又はビームフォーミング法を用いて到来波の到来方向を推定する例を挙げて説明したが、これに限らず、到来波の到来方向を推定できる方法であれば、どのような方法であってもよい。
【００９０】
（９）上述した実施形態及び変形例では、送信部として設けられたトランスデューサ又はアンテナから送信された送信波の反射波を受信し、当該反射波の方位を推定した。すなわち、上述した実施形態及び変形例では、本発明の適用例として、いわゆるアクティブ方式の探知装置を例に挙げて説明したが、これに限らない。具体的には、本発明は、送信素子を有さずに受信素子のみを有する、いわゆるパッシブ方式の探知装置に適用することもできる。
【００９１】
（１０）上記実施形態及び上記変形例では、探知装置の具体例として、水中探知装置及びレーダ装置を例に挙げて説明したが、これに限らず、例えば、超音波診断装置等の医療機器に適用することもできる。
【産業上の利用可能性】
【００９２】
本発明は、複数の受信素子のそれぞれで受波された受信波から得られる受信信号を処理する信号処理装置、当該信号処理装置を備えた水中探知装置、及びレーダ装置として広く適用することができるものである。
【符号の説明】
【００９３】
１，１ａ，１ｂ          水中探知装置（探知装置）
１ｃ                    レーダ装置（探知装置）
３０，３０ａ〜３０ｃ    信号処理部（信号処理装置）
３０ｄ                  信号処理装置
３１                    信号ベクトル生成部
３４，３９              対象信号ベクトル選択部
３７，３７ａ            到来波演算部（物標方位推定部）
３７ｂ                  物標方位推定部
Ａｊ                    超音波振動子（受信素子）
Ｂｊ                    アンテナ素子（受信素子）
【図１】
0006403327.tif 000003
【図２】
0006403327.tif 000004
【図３】
0006403327.tif 000005
【図４】
0006403327.tif 000006
【図５】
0006403327.tif 000007
【図６】
0006403327.tif 000008
【図７】
0006403327.tif 000009
【図８】
0006403327.tif 000010
【図９】
0006403327.tif 000011
【図１０】
0006403327.tif 000012