(57)【特許請求の範囲】
【請求項１】
受信信号を周波数分析し、該周波数分析により分析された周波数ごとに方位ベクトルを算出し、該算出された周波数ごとの方位ベクトルから目標信号の到来方位を算出する目標信号検出方法であって、
予め設定された下限周波数と上限周波数の間の前記方位ベクトルに対して、予め設定され、周波数帯域幅より狭い幅の周波数シフト幅ずつ前記下限周波数と上限周波数との間をシフトさせながら、各シフト位置において予め設定された前記周波数帯域幅内に含まれる前記方位ベクトルを加算して加算ベクトルを算出し、前記各シフト位置で算出された加算ベクトルの中から最大のベクトル値を示す加算ベクトルを選択し、該選択した加算ベクトルの方位及びベクトル値を目標の有意情報として出力することを特徴とする目標信号検出方法。
【請求項２】
前記下限周波数と上限周波数、および前記周波数帯域幅は、対象とする目標信号入力が想定される下限周波数と上限周波数、および対象とする目標信号からの周波数が集中する周波数帯域幅に設定されることを特徴とする請求項１に記載の目標信号検出方法。
【請求項３】
前記周波数帯域幅内に含まれる前記方位ベクトルの加算処理は、前記周波数ごとの方位ベクトルの大きさを全て１に正規化して加算することを特徴とする請求項１または２に記載の目標信号検出方法。
【請求項４】
前記出力された目標の有意情報を、時間軸と方位軸の平面上に輝度レベルによって表示することを特徴とする請求項１〜３のいずれか１項に記載の目標信号検出方法。
【請求項５】
前記受信信号は、水中音響信号であることを特徴とする請求項１〜４のいずれか１項に記載の目標信号検出方法。
【請求項６】
下限周波数と上限周波数、周波数帯域幅、および、前記周波数帯域幅より狭い幅の周波数シフト幅をそれぞれ設定する目標信号処理周波数領域設定手段と、
前記下限周波数と上限周波数の間を前記周波数シフト幅ずつシフトさせながら、各シフト位置において、前記周波数帯域幅内に含まれる受信信号について周波数ごとの方位ベクトルを算出して加算するベクトル加算手段と、
前記各シフト位置で求められた加算ベクトルの中から最大のベクトル値を示す加算ベクトルを選択する最大加算ベクトル選択手段と、
該選択された加算ベクトルの方位及び加算ベクトル値を目標の有意情報として出力する目標信号出力手段と、
を備えていることを特徴とする目標信号検出装置。
【請求項７】
前記ベクトル加算手段は、前記方位計算前処理部から出力される前記周波数ごとの方位ベクトルの大きさを全て１に正規化して加算処理を行う手段から成ることを特徴とする請求項６に記載の目標信号検出装置。
【請求項８】
前記目標信号出力手段から出力された目標の有意情報を入力して、時間軸と方位軸からなる平面上に輝度レベルに変換した表示を行う表示部を備えていることを特徴とする請求項６または７に記載の目標信号検出装置。
【請求項９】
前記受信信号は、水中音響信号であることを特徴とする請求項６〜８のいずれか１項に記載の目標信号検出装置。
【請求項１０】
周波数ごとの方位ベクトルに周波数分析された受信信号から目標信号の到来方位を算出する処理をコンピュータに実行させるためのプログラムであって、
該コンピュータに、予め設定された下限周波数と上限周波数の間の前記方位ベクトルに対して、予め設定され、周波数帯域幅より狭い幅の周波数シフト幅ずつ前記下限周波数と上限周波数との間をシフトさせながら、各シフト位置において予め設定された前記周波数帯域幅内に含まれる前記方位ベクトルを加算して加算ベクトルを算出する処理と、前記各シフト位置で算出された加算ベクトルの中から最大のベクトル値を示す加算ベクトルを選択する処理と、該選択した加算ベクトルの方位及びベクトル値を目標の有意情報として出力する処理を実行させることを特徴とするプログラム。
【請求項１１】
前記加算ベクトルを算出する処理は、前記周波数ごとの方位ベクトルの大きさを全て１に正規化して加算する処理であることを特徴とする請求項１０に記載のプログラム。
【請求項１２】
前記出力された目標の有意情報を、時間軸と方位軸の平面上に、輝度レベルに変換して表示する処理を含むことを特徴とする請求項１０または１１に記載のプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、目標信号検出方法、装置およびプログラムに関し、特に、受信した水中音響信号を解析して表示する技術に関する。
【背景技術】
【０００２】
水中に存在する音源から発生する音響信号を受信して信号処理することにより、音響信号の到来方位を探知し、該到来方位情報を輝度レベルに変換して、方位軸と時間軸の平面上に表示する技術が特許文献１等で提案されている。
【０００３】
図４は、本発明に関連する受信音響信号の処理の例を示すブロック図であり、図５は、音響信号の到来方位算出方法の一例を示すブロック図である。
【０００４】
図４（ａ）に示す目標信号表示処理装置は、複数の指向性ビームにより受信した音響信号のレベルを計測するレベル計算部４０１と、各ビームの音響信号のレベルや位相情報により、方位を計測する方位計算処理部４０２と、レベル計算部４０１から出力されたレベル計測結果と方位計算処理部４０２から出力された方位計測結果により、方位軸と時間軸の平面上に目標信号を輝度で表示する表示処理部４０３により、構成されている。
【０００５】
また、図４（ｂ）に示す別の音響信号表示処理装置は、複数の指向性ビームにより受信した音響信号を周波数分析する周波数分析部４１０と、周波数ごとに音響信号のレベルを計測するレベル計算部４１１と、各ビームの音響信号のレベルや位相情報により、周波数ごとに方位を計算する方位計算処理部４１２と、レベル計算部４１１から出力された周波数ごとのレベル計測結果と方位計算処理部４１２から出力された周波数ごとの方位計測結果により、レベルの高いものを抽出し、方位軸と時間軸の平面上に目標信号を輝度で表示する表示処理部４１３により、構成されている。
【０００６】
音響信号の到来方位の算出は、図５に示すように、全方位に同一感度である（以降、無指向性と呼ぶ）ＯＭＮＩ信号と南北方向に８の字の指向性を形成するＮＳ信号及び東西方向に８の字の指向性を形成するＥＷ信号のレベル及び位相関係によってその方位を算出することができる。図５において、ＮＳ信号及びＥＷ信号の破線による表示は、ＯＭＮＩ信号と比較した場合に位相が逆になっていることを示している。
【０００７】
例えば図５において、１２０度方向から音響信号が到来した場合には、ＮＳ信号のレベルは、ＯＡ［０．５］であり、ＥＷ信号のレベルは、ＯＢ［０．８６６］である。なお、［  ］内の数字は、ＯＭＮＩ信号のレベルを１としたときのＮＳ信号及びＥＷ信号のレベルの値を示している。また、この場合、位相関係は、ＯＭＮＩ信号を基準としＮＳ信号は逆位相、ＥＷ信号は同位相となっている。
【０００８】
このことから、位相関係により象限が決定され、次の式により方位が算出される。なお、象限は、平面を南北方向の軸及び東西方向の軸で４つに分割したもので、０度から９０度にあたる平面を象限１、９０度から１８０度にあたる平面を象限２、１８０度から２７０度にあたる平面を象限３、２７０度から３６０度にあたる平面を象限４とする。
【０００９】
θ＝｜Ｔａｎ−１（ＯＡ／ＯＢ）｜＋９０（象限数−１）
象限１：ＮＳ信号（同位相）、ＥＷ信号（同位相）
象限２：ＮＳ信号（逆位相）、ＥＷ信号（同位相）
象限３：ＮＳ信号（逆位相）、ＥＷ信号（逆位相）
象限４：ＮＳ信号（同位相）、ＥＷ信号（逆位相）
【００１０】
上記式に、それぞれの値を代入すると、θ≒１２０の結果が得られる。したがって、ＯＭＮＩ信号と南北方向に８の字の指向性を形成するＮＳ信号及び東西方向に８の字の指向性を形成するＥＷ信号のレベルと位相を計測することにより音響信号の到来方位を求めることができる。
【００１１】
また図４（ｂ）に示す、周波数分析を行い周波数ごとに方位を算出する方法の場合、狭帯域分析により広帯域にわたって分布する雑音のレベルを低下させているので、狭帯域の受信信号のときにはＳ／Ｎが改善されるが、受信信号が広帯域信号の場合には信号レベルも各帯域に分散して低下してしまうためにＳ／Ｎ改善の効果が得られない。このような広帯域性の水中音響信号に対して、Ｓ／Ｎを改善して検出を容易にする技術が特許文献２〜３等で提案されている。
【００１２】
例えば特許文献２では、指向性ソノブイの受波信号から周波数毎の信号方位値とパワースペクトル値を生成し、これらの信号方位とパワースペクトル値により周波数毎の信号ベクトルを合成し、これらの信号ベクトルを全周波数で加算して全周波数の和信号ベクトルを求め、次にこの和信号ベクトルのスレッショルドレベルを越える信号を検出する構成とし、信号成分は、全周波数で信号成分方位が等しいから全周波数のベクトル加算によりレベルの総和が得られるとともに信号成分方位が保存され、雑音成分は、各周波数で雑音成分方位がランダムであるため、全周波数のベクトル加算によりレベルが０となることを利用して、Ｓ／Ｎの改善を図っている。
【００１３】
また、特許文献３では、音響信号の周波数セル毎の信号エネルギーを算出するスペクトル分析手段と、該スペクトル分析手段にて算出された信号エネルギーに基づいて前記音響信号の周波数帯域を抽出し、抽出された周波数帯域における信号エネルギーのみを出力する周波数選択手段と、該周波数選択手段から出力された周波数セルの信号エネルギーをスカラー量として合成し、スカラーエネルギーとして出力する全エネルギー算出手段と、前記周波数選択手段から出力された周波数セルの信号エネルギーに基づいて、前記音響信号の到来方位を周波数セル毎に算出する方位算出手段と、前記方位算出手段にて算出された周波数セル毎の信号到来方位によって構成されるベクトルを合成し、ベクトルエネルギーとして出力するベクトル合成手段と、該ベクトル合成手段から出力されたベクトルエネルギーの大きさと前記全エネルギー算出手段から出力されたスカラーエネルギーの大きさとに基づいて、音響信号が発せられる音源の有無を判断する探知判定手段を備えることにより、同様の効果を得ている。
【００１４】
【特許文献１】特開２０００−９８０１６号公報
【特許文献２】特開平８−１５２４６５号公報
【特許文献３】特開平１１−５２０３８号公報
【発明の開示】
【発明が解決しようとする課題】
【００１５】
水中の音響信号を解析する目標信号表示処理装置の中に、一定の周波数帯域幅の特性をもつ信号、特に短時間信号や過渡信号を検出して表示する目標信号表示処理装置があるが、短時間信号や一過性信号の場合、一般的に信号の特性として一定幅の周波数帯域内にエネルギーが分散しており、そのため、上記図４に示す目標信号表示処理装置を用いて、短時間信号や一過性信号の検出処理を行った場合には、方位誤差が大きくなるとともに、表示輝度範囲が小さくなって信号の方位が強調できなくなるという問題がある。
【００１６】
すなわち、図４（ａ）に示す受信信号を直接用いて方位を算出する方法では、一般に受信周波数帯域幅が短時間信号や一過性信号の帯域幅よりも大きいために、雑音成分しかない帯域も受信してしまい、その結果、ＳＮ比が低下し方位誤差（偏差）が大きくなる。またＳＮ比の低下により十分なレベルの輝度情報が得られないために、輝度変化が小さくなって信号の方位を強調できなくなる。
【００１７】
また、図４（ｂ）に示す周波数分析を行い周波数ごとに方位を算出する方法では、信号エネルギーが各周波数に分散され、かつ、レベルの大きい周波数のみの方位（１サンプル）を採用して信号検出を行っているので、一定の周波数帯域幅の特性をもつ信号のエネルギーを十分に利用できず、結果として方位誤差（偏差）が大きくなる。また、周波数分析するために信号の溜め込みを行っているために時間的な平滑化が行われる結果、レベル変動が小さくなり輝度変化が小さくなってしまう。小さなレベル変動にて表示輝度を大きくする方法はあるが、その場合、雑音による変動も強調表示されてしまうために、目標からの信号と雑音との区別がつけ難くなる。
【００１８】
また上記特許文献２〜３に記載の目標信号表示処理装置の場合、基本的に、周波数分析された全帯域を対象にしたベクトル加算を行う構成であるため、連続性のある（常時発生している）広帯域信号を検出する場合には、その効果が発揮されるが、一定の周波数帯域幅の特性をもつ信号、特に短時間信号や過渡信号を検出する場合であっても、信号が存在しない帯域（雑音成分のみの帯域）も含めて全て加算された信号として出力するので、結果としてＳ／Ｎの低下を招くことになる。
【００１９】
本発明の目的は、上述した課題を解決することが可能な、目標信号検出方法、装置およびプログラムを提供することにある。
【課題を解決するための手段】
【００２０】
本発明の目標信号検出方法は、受信信号を周波数分析し、該周波数分析により分析された周波数ごとに方位ベクトルを算出し、該算出された周波数ごとの方位ベクトルから目標信号の到来方位を算出する目標信号検出方法であって、予め設定された下限周波数と上限周波数の間の前記方位ベクトルに対して、予め設定され、周波数帯域幅より狭い幅の周波数シフト幅ずつ前記下限周波数と上限周波数との間をシフトさせながら、各シフト位置において予め設定された前記周波数帯域幅内に含まれる前記方位ベクトルを加算して加算ベクトルを算出し、前記各シフト位置で算出された加算ベクトルの中から最大のベクトル値を示す加算ベクトルを選択し、該選択した加算ベクトルの方位及びベクトル値を目標の有意情報として出力することを特徴とする。
【００２１】
本発明の目標信号検出装置は、下限周波数と上限周波数、周波数帯域幅、および、前記周波数帯域幅より狭い幅の周波数シフト幅をそれぞれ設定する目標信号処理周波数領域設定手段と、前記下限周波数と上限周波数の間を前記周波数シフト幅ずつシフトさせながら、各シフト位置において、前記周波数帯域幅内に含まれる受信信号について周波数ごとの方位ベクトルを算出して加算するベクトル加算手段と、前記各シフト位置で求められた加算ベクトルの中から最大のベクトル値を示す加算ベクトルを選択する最大加算ベクトル選択手段と、該選択された加算ベクトルの方位及び加算ベクトル値を目標の有意情報として出力する目標信号出力手段と、を備えていることを特徴とする。
【００２２】
本発明のプログラムは、周波数ごとの方位ベクトルに周波数分析された受信信号から目標信号の到来方位を算出する処理をコンピュータに実行させるためのプログラムであって、該コンピュータに、予め設定された下限周波数と上限周波数の間の前記方位ベクトルに対して、予め設定され、周波数帯域幅より狭い幅の周波数シフト幅ずつ前記下限周波数と上限周波数との間をシフトさせながら、各シフト位置において予め設定された前記周波数帯域幅内に含まれる前記方位ベクトルを加算して加算ベクトルを算出する処理と、前記各シフト位置で算出された加算ベクトルの中から最大のベクトル値を示す加算ベクトルを選択する処理と、該選択した加算ベクトルの方位及びベクトル値を目標の有意情報として出力する処理を実行させることを特徴とする。
【発明の効果】
【００２３】
第１の効果は、目標信号が短時間信号や一過性信号の場合の方位誤差が小さくなることである。その理由は、方位算出に関し、周波数分析の分解能単位で得られる方位算出結果のベクトル加算及びベクトル値（最大）方位の抽出により、目標信号の集中している方位が有意情報として得られるためである。
【００２４】
第２の効果は、信号が強調され表示されることである。その理由は、方位算出に関し、周波数分析の分解能単位で得られる方位算出結果の信号帯域のベクトル加算及び方位集中度の算出結果を輝度として表示することにより、雑音に関しては、周波数ごとの方位算出結果がランダムな結果となるためにベクトル加算した場合にベクトルは小さく（ゼロに近づく）なり、信号に関しては、特定の方位に集中するためにベクトルは大きくなるためである。
【発明を実施するための最良の形態】
【００２５】
図１は、本発明の実施形態を示す目標信号表示処理装置の機能ブロック図である。
【００２６】
本実施形態の目標信号表示処理装置１００は、周波数分析部１０１、方位計算前処理部１０２、方位計算後処理部１０３、及び表示処理部１０４にて構成される。周波数分析部１０１及び方位計算前処理部１０２は、図４（ｂ）における周波数分析部４１０及び方位計算処理部４１２と同様の回路である。
【００２７】
周波数分析部１０１は、受信した音響信号を入力し、周波数分析を行い方位計算前処理部１０２に周波数分析結果を出力する。方位計算前処理部１０２は、入力した周波数分析結果から、図５で示した各ビームのレベル及び位相関係から周波数ごとに方位計算を行い方位計算後処理部１０３に周波数ごとの方位計算結果を出力する。
【００２８】
方位計算後処理部１０３は、入力した周波数ごとの方位算出結果から、対象とする目標信号の周波数帯域幅ごとの方位とその集中度を計算して求め、その中から集中度が最大となる方位及びその集中度を対象信号の有意情報として抽出し、表示処理部１０４に抽出した方位と集中度を出力する。表示処理部１０４は、入力した方位と集中度をもとに方位軸と時間軸の平面上に集中度を輝度情報に変換して表示させる。
【００２９】
この目標信号表示処理装置１００は、例えば図２に示すように、データバス１１０、受信音響信号をデジタルデータに変換してデータバスへ接続するＡ／Ｄ変換器１１１、ＣＰＵ１１２、ＲＯＭおよびＲＡＭを有する記憶部１１３、表示部２００及び入力部３００とデータバス間を接続する入出力インタフェース部１１４を含む構成として実現することができる。
【００３０】
この場合、デジタルデータに変換された受信音響信号は記憶部１１３のＲＡＭに一旦格納される。記憶部１１３のＲＯＭには、本実施形態の目標信号表示処理を実行するためのプログラムが格納されており、ＣＰＵ１１２は、このＲＯＭに格納されたプログラムを読み込み、このプログラムに従って、ＲＡＭに格納されているデジタル受信音響信号に対して本実施形態の目標信号表示処理を実行する。したがって、ＣＰＵ１１２および記憶部１１３は、本実施形態の目標信号表示処理プログラムを実行するためのコンピュータとしての機能を有している。
【００３１】
図３は、本実施形態の方位計算後処理部１０３の処理概要を示す図である。以下、図１〜図３を参照して本実施形態の動作について説明する。
【００３２】
周波数分析部１０１は、受信した音響信号を入力し、周波数分析を行い方位計算前処理部１０２に周波数分析結果を出力する。ＯＭＮＩ信号、ＮＳ信号及びＥＷ信号から方位を算出する際には３つの信号それぞれを周波数分析して結果を出力することになる。また、全周に複数の指向性ビームを形成し方位を算出するようなシステムである場合には、全てのビームについて周波数分析を行う。
【００３３】
方位計算前処理部１０２は、周波数分析部１０１から入力された周波数分析結果を元にして、各ビームのレベル及び位相関係から周波数ごとに方位計算を行って、方位計算後処理部１０３に周波数ごとの方位計算結果（周波数ごとの方位ベクトル）を出力する。全周に複数の指向性ビームを形成し方位を算出する場合には、レベルのみで方位を算出する方法を用いる。簡単に説明すると、レベルの高いビームを抽出しその両側の信号レベルの関係から二次曲線近似等によりピークレベルをもとめることで周波数ごとの方位ベクトルを算出する方法等を採用することができる。
【００３４】
方位計算後処理部１０３は、方位計算前処理部１０２から入力した周波数ごとの方位ベクトルを用いて、一定の周波数帯域幅ごとの方位とその集中度を計算し、集中度が最大となる方位及びその集中度を対象信号の有意情報として抽出し、表示処理部１０４に、抽出した方位と集中度を出力する。方位計算後処理部１０３の動作を、図３を用いて具体的に説明する。
【００３５】
受信した信号から対象信号（過渡信号）にあわせた周波数帯域のみを処理させるために対象信号（過渡信号）の入力が想定される下限周波数と上限周波数を設定する（例えば下限周波数を２００Ｈｚ、上限周波数を２０００Ｈｚなど）。また、対象信号（過渡信号）の周波数帯域幅を設定する（例えば周波数帯域幅２００Ｈｚなど）。さらに上下限周波数範囲内の帯域周波数間をシフトさせながら方位ベクトル加算処理を行うために周波数シフト幅を設定する（例えば周波数シフト幅１００Ｈｚなど）（Ｓ１）。
【００３６】
この周波数シフト幅は、大きすぎる場合、対象信号を適切に捕らえられない可能性があり、また小さすぎる場合には、処理負荷が大きくなるので、処理負荷と対象信号捕捉効果のバランスにより決定する必要がある。これらの設定値は、操作者が入力部３００から設定することが可能であり、設定データは記憶部１１３のＲＡＭ等に格納される。ＣＰＵ１１２は本実施形態の処理プログラムを実行する際にこの設定値をＲＡＭから読み出して参照する。
【００３７】
パラメータ設定完了後、方位計算前処理部１０２から入力した周波数ごとの方位ベクトルの中から、２００Ｈｚ〜４００Ｈｚの範囲（最初の周波数帯域幅２００Ｈｚ）に含まれる周波数ごとの方位ベクトルの加算処理を実施して、加算ベクトルの方向（信号方位）と加算ベクトルの大きさ（方位集中度）を求める（Ｓ２）。方位が、特定方位に集中している場合は、方向は平均的な方位を示し、ベクトルの大きさは方位の集中度を示すことになる。一方、方位が雑音などで一様に分布している場合には、加算ベクトルは中心部に近づき、ベクトルの大きさ（方位集中度）が小さくなる。
【００３８】
なお、本実施例の方位ベクトル加算処理では、周波数ごとに持つ方位計算結果（方位ベクトル）について、ベクトルの大きさを全て「１」に正規化して、方位方向に加算していく処理を行っている。これは、加算される周波数帯域内の特定の周波数に大きな信号（狭帯域信号）が存在した場合、受信対象信号（過渡信号）の方位がこの狭帯域信号の方位に依存してしまう可能性を防ぐためであるが、そのような影響が少ない場合には、周波数ごとに持つ方位計算結果（方位ベクトル）をそのまま加算してもよい。
【００３９】
続いて、設定パラメータに基づき、周波数を１００Ｈｚシフトして（Ｓ３）、方位計算前処理部１０２から入力した周波数ごとの方位ベクトルの中から、３００Ｈｚ〜５００Ｈｚの範囲（次の周波数帯域幅２００Ｈｚ）に含まれる周波数ごとの方位ベクトルの加算処理を実施して、加算ベクトルの方向（信号方位）と加算ベクトルの大きさ（方位集中度）を求める（Ｓ２）。
【００４０】
以下同様にして、方位計算前処理部１０２から入力した周波数ごとの方位ベクトルの中から、周波数を１００Ｈｚシフトしながら周波数帯域幅２００Ｈｚの範囲に含まれる周波数ごとの方位ベクトルの加算処理を、１８００Ｈｚ〜２０００Ｈｚの範囲に含まれる周波数ごとの方位ベクトルの加算処理が終了するまで、合計１６回の方位ベクトル加算処理を繰り返す（Ｓ２〜Ｓ４）。
【００４１】
設定された上限周波数（２０００Ｈｚ）までの方位ベクトル加算処理が終了した（Ｓ４，ＹＥＳ）後、求められた１６個の加算ベクトル値の中で、ベクトルの大きさが最大となる加算ベクトルの値及び加算ベクトルの方向を有意情報として選択するとともに、最大加算ベクトル値が求められた周波数帯域を抽出する（Ｓ５）。
【００４２】
有意情報として選択された加算ベクトルの方向および加算ベクトルの値は、目標の方位および集中度として表示処理部１０４に出力される（Ｓ６）。表示処理部１０４は、入力した方位と集中度をもとに方位軸と時間軸の平面上に集中度を輝度として表示するためのデータを生成して表示部２００へ出力する。表示部２００は、表示処理部１０４からのデータを入力して、方位軸と時間軸の平面上に前記集中度を輝度レベルに変換して表示する。
【００４３】
なお、本実施例において（  ）内に示された設定パラメータは一例であってこれに限定されるものではなく、下限周波数と上限周波数、周波数帯域幅、および周波数シフト幅は、受信対象信号（過渡信号）の予想される周波数範囲、周波数帯域特性、および処理能力等に合わせて、適宜変更して設定することができる。
【００４４】
本実施形態によれば、予め受信対象信号（過渡信号）の周波数帯域特性および周波数範囲を予想して方位ベクトル加算処理を行っているので、対象信号検出効率を高めることができ、かつ方位ベクトル加算処理によりＳ／Ｎも改善されるので、表示輝度範囲を大きくすることができ信号の方位を強調して表示することができる。
【００４５】
なお、上記実施形態は、受信信号として水中に存在する音源から発生する音響信号をその対象として適用したが、本発明は、空中に存在する音源から発生する音響信号、あるいはレーダ信号等の電波信号に対しても適用することができる。
【図面の簡単な説明】
【００４６】
【図１】本発明の実施形態を示す目標信号表示処理装置の機能ブロック図である。
【図２】本実施形態の目標信号表示処理装置をソフトウェアにより実現する場合の構成を示す概略ブロック図である。
【図３】本実施形態の方位計算後処理部の処理概要を示す図である。
【図４】本発明に関連する構成を示す図である。
【図５】方位算出方法の一例の概念を示す図である。
【符号の説明】
【００４７】
１００  目標信号表示処理装置
１０１、４１０  周波数分析部
１０２  方位計算前処理部
１０３  方位計算後処理部
１０４，４０３，４１３  表示処理部
４０１，４１１  レベル計算部
４０２，４１２  方位計算処理部
【図１】
0004900229.tif 000002
【図２】
0004900229.tif 000003
【図３】
0004900229.tif 000004
【図４】
0004900229.tif 000005
【図５】
0004900229.tif 000006