(57)【要約】【課題】誤警報の回数を抑制できる信号追尾装置を得る。【解決手段】入力される信号強度分布データに対して、方位毎にセンサアレイの基準方位からの相対方位の半減半角ビーム幅に基づくフィルタ長のフィルタ係数を用いて、不要ピークを除去する不要ピーク除去手段と、不要ピーク除去後の信号強度分布データを用いて、隣接する方位の信号強度を相対的に比較してピークを検出するピーク検出手段と、前回の周期までに検出されたラインとピーク検出手段によって検出されたピークとが方位方向の予測範囲に収まる場合に検出されたピークをラインに統合し、いずれのラインにも統合されないピークを新規のラインに設定する方位追尾手段と、ラインが時間的に連ねたピークのレベルを時間方向に平滑化した値が決められた検出閾値以上である場合にラインを信号成分と判定する検出判定手段とを有するものである。【選択図】図１
2020085525.tif 000002
【特許請求の範囲】
【請求項１】
センサアレイに到来する信号に対して一定の周期で整相処理および広帯域分析処理を行うことで得られる、方位空間上の信号強度を示す信号強度分布データを用いた信号追尾装置であって、
前記信号強度分布データに対して、方位毎に前記センサアレイの基準方位からの相対方位の半減半角ビーム幅に基づくフィルタ長のフィルタ係数を用いて、前記信号強度分布データから不要ピークを除去する不要ピーク除去手段と、
前記不要ピーク除去後の信号強度分布データを用いて、隣接する方位の信号強度を相対的に比較してピークを検出するピーク検出手段と、
前回の周期までに検出された前記ピークの時間的なつながりであるラインと前記ピーク検出手段によって検出されたピークとが方位方向の予測範囲に収まる場合に前記検出されたピークを前記ラインに統合し、前記検出されたピークのうち、いずれのラインにも統合されないピークを新規のラインに設定する方位追尾手段と、
前記ラインが時間的に連ねたピークのレベルを時間方向に平滑化した値が決められた検出閾値以上である場合に前記ラインを信号成分と判定する検出判定手段と、
を有する信号追尾装置。
【請求項２】
前記不要ピーク除去手段は、
前記方位をｍとし、時間ｔにおける前記信号強度分布データをＸｔ（ｍ）とし、前記信号強度分布データの長さをＭとし、前記方位ｍにおける前記半減半角ビーム幅をＤｒｅｌａｔｉｖｅ（ｍ）とし、調整係数をＣとし、前記フィルタ長を（２×ｓ（ｍ）＋１）とし、ｓ（ｍ）＝四捨五入（Ｄｒｅｌａｔｉｖｅ（ｍ）×Ｃ÷３６０×Ｍ）とし、前記フィルタ係数をＣ＿ｆｉｌｔｅｒとするとき、次式
【数１】
2020085525.tif 000004
に基づいて、前記不要ピーク除去後の信号強度分布データであるＸｆｉｌｔｅｒｔ（ｍ）を生成する、請求項１に記載の信号追尾装置。
【請求項３】
センサアレイに到来する信号に対して一定の周期で整相処理および広帯域分析処理を行うことで得られる、方位空間上の信号強度を示す信号強度分布データを用いた信号追尾方法であって、
前記信号強度分布データに対して、方位毎に前記センサアレイの基準方位からの相対方位の半減半角ビーム幅に基づくフィルタ長のフィルタ係数を用いて、前記信号強度分布データから不要ピークを除去し、
前記不要ピーク除去後の信号強度分布データを用いて、隣接する方位の信号強度を相対的に比較してピークを検出し、
前回の周期までに検出された前記ピークの時間的なつながりであるラインと検出されたピークとが方位方向の予測範囲に収まる場合に前記検出されたピークを前記ラインに統合し、前記検出されたピークのうち、いずれのラインにも統合されないピークを新規のラインに設定し、
前記ラインが時間的に連ねたピークのレベルを時間方向に平滑化した値が決められた検出閾値以上である場合に前記ラインを信号成分と判定する、
信号追尾方法。
【請求項４】
前記信号強度分布データから前記不要ピークを除去する際、
前記方位をｍとし、時間ｔにおける前記信号強度分布データをＸｔ（ｍ）とし、前記信号強度分布データの長さをＭとし、前記方位ｍにおける前記半減半角ビーム幅をＤｒｅｌａｔｉｖｅ（ｍ）とし、調整係数をＣとし、前記フィルタ長を（２×ｓ（ｍ）＋１）とし、ｓ（ｍ）＝四捨五入（Ｄｒｅｌａｔｉｖｅ（ｍ）×Ｃ÷３６０×Ｍ）とし、前記フィルタ係数をＣ＿ｆｉｌｔｅｒとするとき、次式
【数２】
2020085525.tif 000005
に基づいて、前記不要ピーク除去後の信号強度分布データであるＸｆｉｌｔｅｒｔ（ｍ）を生成する、請求項３に記載の信号追尾方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の受波器で構成されるセンサアレイに到来する信号を分析して信号成分を検出する信号追尾装置および信号追尾方法に関する。
【背景技術】
【０００２】
従来の信号追尾方法として、センサアレイに到来する信号であるアレイ入力データに対して整相処理および広帯域分析処理を行うことにより得られる方位空間上の信号強度分布データを用いて、到来する信号の方位を検出する方法が提案されている（例えば、非特許文献１参照）。
【０００３】
近年、信号成分に対する検出能力の向上を目的として、従来の広帯域分析処理の代わりに適応型広帯域分析処理が採用されるようになった。適応型広帯域分析処理による信号強度分布データを使用することで、従来の広帯域分析処理に比べて、より低いＳ／Ｎ比（信号対雑音比）の信号の検出が可能になったり、異なる音源から発生する信号同士の方位の分離能力が向上する。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Robert E. Zarnich, “A Fresh Look At ‘Broadband’ Passive Sonar Processing”, in Proceedings of 1999 Adaptive Sensor Array Processing Workshop(ASAP `99), MIT Lincoln Laboratory, Lexington, MA, pp. 99-104.
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかし、適応型広帯域分析処理による信号強度分布データは、信号成分だけでなく、雑音を含むピークが全周に渡り散在しやすくなる。そのため、雑音によるピークが信号成分と誤検出され、誤警報が多発してしまうおそれがある。
【課題を解決するための手段】
【０００６】
本発明に係る信号追尾装置は、センサアレイに到来する信号に対して一定の周期で整相処理および広帯域分析処理を行うことで得られる、方位空間上の信号強度を示す信号強度分布データを用いた信号追尾装置であって、前記信号強度分布データに対して、方位毎に前記センサアレイの基準方位からの相対方位の半減半角ビーム幅に基づくフィルタ長のフィルタ係数を用いて、前記信号強度分布データから不要ピークを除去する不要ピーク除去手段と、前記不要ピーク除去後の信号強度分布データを用いて、隣接する方位の信号強度を相対的に比較してピークを検出するピーク検出手段と、前回の周期までに検出された前記ピークの時間的なつながりであるラインと前記ピーク検出手段によって検出されたピークとが方位方向の予測範囲に収まる場合に前記検出されたピークを前記ラインに統合し、前記検出されたピークのうち、いずれのラインにも統合されないピークを新規のラインに設定する方位追尾手段と、前記ラインが時間的に連ねたピークのレベルを時間方向に平滑化した値が決められた検出閾値以上である場合に前記ラインを信号成分と判定する検出判定手段と、を有するものである。
【０００７】
本発明に係る信号追尾方法は、センサアレイに到来する信号に対して一定の周期で整相処理および広帯域分析処理を行うことで得られる、方位空間上の信号強度を示す信号強度分布データを用いた信号追尾方法であって、前記信号強度分布データに対して、方位毎に前記センサアレイの基準方位からの相対方位の半減半角ビーム幅に基づくフィルタ長のフィルタ係数を用いて、前記信号強度分布データから不要ピークを除去し、前記不要ピーク除去後の信号強度分布データを用いて、隣接する方位の信号強度を相対的に比較してピークを検出し、前回の周期までに検出された前記ピークの時間的なつながりであるラインと検出されたピークとが方位方向の予測範囲に収まる場合に前記検出されたピークを前記ラインに統合し、前記検出されたピークのうち、いずれのラインにも統合されないピークを新規のラインに設定し、前記ラインが時間的に連ねたピークのレベルを時間方向に平滑化した値が決められた検出閾値以上である場合に前記ラインを信号成分と判定するものである。
【発明の効果】
【０００８】
本発明によれば、信号強度分布データに対して、方位毎の半減半角ビーム幅に基づくフィルタ長のフィルタ係数を用いることで不要ピークが除去されるため、雑音によるピークの誤検出が抑制され、誤警報の回数を抑制することができる。
【図面の簡単な説明】
【０００９】
【図１】本発明の実施の形態１に係る信号追尾装置の一構成例を示すブロック図である。
【図２】図１に示した信号追尾装置のハードウェア構成の一例を示す図である。
【図３】アレイ入力データに対して整相処理および広帯域分析処理を行って信号強度分布データを得る方法を説明するための図である。
【図４】比較例１の広帯域分析処理による信号強度分布データの一例を示す図である。
【図５】比較例１の信号追尾装置の一構成例を示すブロック図である。
【図６】比較例１の信号追尾装置が実行する信号追尾方法の手順を示すフローチャートである。
【図７】比較例２の適応型広帯域分析処理による信号強度分布データの一例を示す図である。
【図８】図１に示した信号追尾装置が実行する信号追尾方法の手順を示すフローチャートである。
【図９】本発明の実施の形態１に係る信号追尾方法において、不要ピーク除去前の信号強度分布データＸｔ（ｍ）の一例を示す図である。
【図１０】図９に示した信号強度分布データから不要ピークを除去した後の信号強度分布データの一例を示す図である。
【発明を実施するための形態】
【００１０】
実施の形態１．
本実施の形態１の信号追尾装置および信号追尾方法は、アレイ入力データに対して整相処理および広帯域分析処理を行うことで得られる方位空間上の信号強度分布データを用いて、時間的に連続性のある広帯域信号成分を自動的に検出および追尾するものである。本実施の形態１の信号追尾装置および信号追尾方法は、例えば、潜水艦または水上艦等のソーナー装置に用いられる。
【００１１】
本実施の形態１の信号追尾装置の構成を説明する。図１は、本発明の実施の形態１に係る信号追尾装置の一構成例を示すブロック図である。図１に示すように、信号追尾装置１は、制御部１１および表示部１２を有する。表示部１２は、例えば、液晶ディスプレイ装置である。制御部１１は、データ入力手段２１と、不要ピーク除去手段２２と、ピーク検出手段２３と、方位追尾手段２４と、検出判定手段２５と、表示制御手段２６とを有する。
【００１２】
図２は、図１に示した信号追尾装置のハードウェア構成の一例を示す図である。信号追尾装置１は、例えば、コンピュータを含む情報処理装置である。制御部１１は、プログラムを記憶するメモリ５２と、プログラムにしたがって処理を実行するＣＰＵ（Central Processing Unit）５１とを有する。ＣＰＵ５１がプログラムを実行することで、図１に示した、データ入力手段２１、不要ピーク除去手段２２、ピーク検出手段２３、方位追尾手段２４、検出判定手段２５および表示制御手段２６が構成される。
【００１３】
なお、データ入力手段２１、不要ピーク除去手段２２、ピーク検出手段２３、方位追尾手段２４、検出判定手段２５および表示制御手段２６が備える機能のうち、一部または全部が専用回路で構成されてもよい。専用回路は、例えば、ＡＳＩＣ（Application Specific Integrated Circuit）である。信号処理の一部または全部を専用回路で構成することで、信号処理の高速化を図ることができる。
【００１４】
本実施の形態１の信号追尾装置１の構成を詳細に説明する前に、比較例１の信号追尾装置を説明する。比較例１の信号追尾装置は、アレイ入力データに対して広帯域分析処理を用いるものである。はじめに、信号追尾装置の処理対象となる信号強度分布データについて説明する。図３は、アレイ入力データに対して整相処理および広帯域分析処理を行って信号強度分布データを得る方法を説明するための図である。
【００１５】
整相方位数を全周でＢ個とすると、センサアレイを基準とした所定の相対方位に対して行われた整相処理の結果に対し、方位方向の補間処理が行われる。補間処理により、整相処理後のデータは等θ間隔のデータに変換される。ここで、補間後の全周の方位分割数をＭとすると、Ｂ＜Ｍの関係である。センサアレイが直線状の場合、整相処理結果による半減半角ビーム幅は、センサアレイからの相対方位に応じてその大きさは異なる。
【００１６】
その後、周波数ビン数をＮとすると、等θ間隔の各方位ｍ（ｍ＝１，２，・・・，Ｍ）のグリッドで周波数分析を行うことで２次元データＷ（ｍ，ｎ）（ｍ＝１，２，・・・，Ｍ、ｎ＝１，２，・・・，Ｎ）が生成される。さらに、１≦Ｎ１＜Ｎ２≦Ｎとすると、２次元データＷ（ｍ，ｎ）に対して各方位ｍに所望の周波数ビンＮ１〜Ｎ２の強度を全て加算することで、信号強度分布データＸｔ（ｍ）が生成される。そのため、信号成分が含まれる周波数ビンが限定的であるほど検出能力は劣化する。また、全周の方位分割数Ｍを大きくしても、図４に示すように、方位空間上に形成されるピークの分解能は整相方位の分解能程度に留まる。図４は、比較例１の広帯域分析処理による信号強度分布データの一例を示す図である。図４の横軸は方位ｍであり、縦軸は信号強度を示すレベルである。
【００１７】
比較例１の信号追尾装置の構成を説明する。図５は、比較例１の信号追尾装置の一構成例を示すブロック図である。図１に示した構成と同様な構成について同一の符号を付し、その詳細な説明を省略する。比較例１の信号追尾装置１００は、制御部１３および表示部１２を有する。制御部１３は、データ入力手段２１、ピーク検出手段２３、方位追尾手段２４、検出判定手段２５および表示制御手段２６を有する。また、制御部１３は、図に示さない、メモリおよびＣＰＵを有する。メモリが記憶するプログラムをＣＰＵが実行することで、データ入力手段２１、ピーク検出手段２３、方位追尾手段２４、検出判定手段２５および表示制御手段２６が構成される。
【００１８】
データ入力手段２１は、一定の周期のサンプリング時間ｔに信号強度分布データＸｔ（ｍ）が入力されると、信号強度分布データＸｔ（ｍ）をピーク検出手段２３に出力する。ただし、ｍ＝１，２，・・・，Ｍである。
【００１９】
ピーク検出手段２３は、信号強度分布データＸｔ（ｍ）に対して、隣接する方位の信号強度を相対的に比較してピークを検出する。具体的には、ピーク検出手段２３は、信号強度分布データＸｔ（ｍ）の各方位ｍに対し、次の式（１）を用いて、ピークか否かを判定する。ピーク検出手段２３は、式（１）が満たされる場合、Ｘｔ（ｍ）をピークと判定し、式（１）が満たされない場合、Ｘｔ（ｍ）をピークではないと判定する。
【００２０】
Ｘｔ（ｍ−１）＜Ｘｔ（ｍ）、かつＸｔ（ｍ）＞Ｘｔ（ｍ＋１）  ・・・（１）
なお、Ｘｔ（ｍ）が観測点を中心にして全周（３６０度）のデータである場合、Ｘｔ（１）とＸｔ（Ｍ）とが隣接する。そのため、ピーク検出手段２３は、Ｘｔ（１）についてピーク判定する場合、式（１）におけるＸｔ（０）をＸｔ（Ｍ）に置き換える。また、ピーク検出手段２３は、Ｘｔ（Ｍ）についてピーク判定する場合、式（１）におけるＸｔ（Ｍ＋１）をＸｔ（１）に置き換える。
【００２１】
方位追尾手段２４は、前回の周期（ｔ−１）までに検出されたピークの時間的なつながりであるラインがある場合、ラインと現在の周期のサンプリング時間ｔで検出されたピークとが同一であるか否かを、統合基準にしたがって判定する。統合基準は、ラインに関して予め算出した予測方位ｄｐｒｅｄを基準とし、所定の方位差ｄｓｅａｒｃｈの予測範囲内に新たに検出されたピークの方位ｄｐｅａｋが収まることである。統合基準を式（２）に示す。
（ｄｐｒｅｄ−ｄｓｅａｒｃｈ） ≦ ｄｐｅａｋ ≦ （ｄｐｒｅｄ＋ｄｓｅａｒｃｈ）  ・・・（２）
予測方位ｄｐｒｅｄの算出方法は、例えば、カルマンフィルタおよびαβトラッキング法などを用いることが考えられる。
【００２２】
方位追尾手段２４は、式（２）に示す統合基準を満たす場合、前回の周期（ｔ−１）から更新されるラインと現在の周期ｔで検出されたピークとが同一と判定し、検出されたピークをラインに統合させる。一方、方位追尾手段２４は、現在の周期ｔで検出されたピークのうち、どのラインにも統合されないピークを新規ラインに設定する。各ラインは、情報として、「未検出」、「検出」および「失探」のうち、いずれかのステータスと、信号レベル平滑値を持つ。新規ラインは、ステータスが「未検出」であり、信号レベル平滑値は０である。
【００２３】
検出判定手段２５は、方位追尾手段２４によって更新されたラインおよび方位追尾手段２４によって設定された新規ラインの各ラインについて、時間方向に信号レベルを平滑化する。具体的には、検出判定手段２５は、各ラインの現在の周期ｔにおける信号レベル平滑値Ｐｓｍｏｏｔｈ（ｔ）を、統合されたピークのレベルＸｔ（ｍ）、信号レベル平滑値の前回の値Ｐｓｍｏｏｔｈ（ｔ−１）および信号レベルの平滑化の所望の時定数αを用いて算出する。検出判定手段２５が用いる算出式を、式（３）に示す。
Ｐｓｍｏｏｔｈ（ｔ）＝α×Ｘｔ（ｍ）＋（１−α）×Ｐｓｍｏｏｔｈ（ｔ−１）・・・（３）
【００２４】
続いて、検出判定手段２５は、信号レベル平滑値Ｐｓｍｏｏｔｈ（ｔ）と判定基準となる検出閾値Ｖｔｈとを比較する。比較の結果、
Ｐｓｍｏｏｔｈ（ｔ） ≧ Ｖｔｈ    ・・・（４）
である場合、検出判定手段２５は、ラインのステータスを、次のように更新する。
・前回のステータスが「未検出」である場合、ステータスを「検出」に更新する。
・前回のステータスが「検出」である場合、ステータスを「検出」に維持する。
【００２５】
一方、比較の結果、
Ｐｓｍｏｏｔｈ（ｔ） ＜ Ｖｔｈ    ・・・（５）
である場合、検出判定手段２５は、ラインのステータスを、次のように更新する。
・前回のステータスが「未検出」である場合、ステータスを「未検出」に維持する。
・前回のステータスが「検出」である場合、ステータスを「失探」に更新する。
【００２６】
検出判定手段２５は、ステータスが「検出」を示すラインの情報のみを表示制御手段２６に出力する。ここで、検出判定手段２５は、ステータスが「失探」となったラインを、信号成分ではない、または信号成分ではなくなったと判断し、そのラインの情報を削除する。
【００２７】
表示制御手段２６は、検出判定手段２５によって「検出」と判定されたラインの方位とその検出時刻に基づいて、例えば、縦軸を時間とし、横軸を方位とした２次元の表示領域中の該当する座標にプロットする画像を表示部１２に表示させる。その画像が表示部１２によって表示させることで、オペレータが目標の存在およびその運動状況を把握することができる。
【００２８】
次に、比較例１の信号追尾装置１００の動作を説明する。図６は、比較例１の信号追尾装置が実行する信号追尾方法の手順を示すフローチャートである。
【００２９】
データ入力手段２１は、入力される信号強度分布データＸｔ（ｍ）をピーク検出手段２３に出力する（ステップＳ２０１）。ｍ＝１，２，・・・，Ｍである。ピーク検出手段２３は、信号強度分布データＸｔ（ｍ）を用いて、隣接する方位の信号強度を相対的に比較してピークを検出する（ステップＳ２０２）。ステップＳ２０３において、方位追尾手段２４は、前回の周期から維持されるラインとピーク検出手段２３によって検出されたピークとが方位方向の予測範囲に収まる場合に新規に検出されたピークをラインに統合する。一方、方位追尾手段２４は、ピーク検出手段２３によって検出されたピークのうち、いずれのラインにも統合されないピークを新規ラインに設定する。
【００３０】
ステップＳ２０４において、検出判定手段２５は、各ラインについて、ピークのレベルを時間方向に平滑化した信号レベル平滑値Ｐｓｍｏｏｔｈ（ｔ）が検出閾値Ｖｔｈ以上であるか否かを判定する。検出判定手段２５は、信号レベル平滑値Ｐｓｍｏｏｔｈ（ｔ）が検出閾値Ｖｔｈ以上である場合、判定対象のラインを信号成分と判定する。表示制御手段２６は、縦軸を時間とし、横軸を方位とした２次元の表示領域に、検出判定手段２５によって「検出」と判定されたラインの方位とその検出時刻とを座標としてプロットした画像を表示部１２に表示させる。信号追尾装置１００はステップＳ２０１〜Ｓ２０５の処理を一定の周期で繰り返す。
【００３１】
次に、比較例２として、適応型広帯域分析処理の場合を説明する。適応型広帯域分析処理は、従来の広帯域分析処理と同様に、２次元データＷ（ｍ，ｎ）（ｍ＝１，２，・・・，Ｍ、ｎ＝１，２，・・・，Ｎ）を入力とする。一方、適応型広帯域分析処理では、方位方向にピークを成す強度、すなわち任意のｍ、ｎに対し、式（６）を満たす２次元データＷ（ｍ，ｎ）が示す強度のみが方位ｍの加算対象となり、信号強度分布データＸｔ（ｍ）が生成される。
Ｗ（ｍ−１，ｎ）＜Ｗ（ｍ，ｎ）、かつＷ（ｍ，ｎ）＞Ｗ（ｍ＋１，ｎ）・・・（６）
【００３２】
そのため、従来の広帯域分析処理に比べて、異なる音源から発生する信号同士の方位の分離能力が向上する一方、信号強度分布データ上のピーク点が全周に渡り散在しやすくなる。全周の方位分割数Ｍが大きくなるにしたがって、この現象が顕著になる。図７は、比較例２の適応型広帯域分析処理による信号強度分布データの一例を示す図である。図７の横軸は方位ｍであり、縦軸は信号強度を示すレベルである。図７に示すグラフでは、図４に示したグラフに比べて、多くのピーク点が散在していることがわかる。
【００３３】
よって、適応型広帯域分析処理により生成された信号強度分布データＸｔ（ｍ）を入力データとして、比較例１の信号追尾装置１００で広帯域信号の検出および追尾を試みると、あらゆる方位においてピーク点が検出されることになる。検出されるピークには、信号成分だけでなく、雑音も多く含まれる可能性がある。その結果、雑音によるピークが信号成分と誤検出され、誤警報が多発する問題が生じる。そのため、適応型広帯域分析処理により生成されたデータに対し、目標信号の方位のみのピーク点を残し、それ以外の方位に散在する不要なピーク点を除去することが求められている。
【００３４】
本実施の形態１の信号追尾装置１は、広帯域分析処理により生成された信号強度分布データに対して、方位方向のローパスフィルタを実行することにより、方位全体における散在する雑音となるピーク点を除去するものである。これにより、誤警報ラインの発生を抑制し、信号成分の広帯域信号を自動的に検出および追尾できる。以下に、本実施の形態１の信号追尾装置１について詳しく説明する。
【００３５】
本実施の形態１の信号追尾装置１の構成を、図１を参照して説明する。ここでは、図５を参照して説明した構成と同様な構成についての詳細な説明を省略する。図１に示した信号追尾装置１は、図５に示した比較例１の信号追尾装置１００と比較すると、不要ピーク除去手段２２が追加された構成である。
【００３６】
不要ピーク除去手段２２は、信号強度分布データＸｔ（ｍ）に対して、方位ｍ毎にセンサアレイの基準方位からの相対方位であるアレイ相対方位の半減半角ビーム幅に基づくフィルタ長のフィルタ係数を用いて、入力されたデータから不要ピークを除去する。具体的には、不要ピーク除去手段２２は、信号強度分布データＸｔ（ｍ）の各方位ｍに対して、式（７）により、方位ｍにおけるアレイ相対方位の半減半角ビーム幅を考慮したフィルタ長によるローパスフィルタを実行することで、Ｘｆｉｌｔｅｒｔ（ｍ）を得る。ここで、ｍ＝１，２，・・・，Ｍである。Ｘｆｉｌｔｅｒｔ（ｍ）は、入力された信号強度分布データＸｔ（ｍ）から不要ピークを除去した後のデータである。式（７）に、Ｘｆｉｌｔｅｒｔ（ｍ）の算出式を示す。なお、Ｘｔ（ｍ）が観測点を中心にして全周（３６０度）のデータである場合、Ｘｔ（１）とＸｔ（Ｍ）とが隣接する。そのため、不要ピーク除去手段２２は、式（７）において、Ｘｔ（０）をＸｔ（Ｍ）に置き換え、Ｘｔ（Ｍ＋１）をＸｔ（１）に置き換える。
【００３７】
【数１】
2020085525.tif 000003
【００３８】
式（７）において、Ｃ＿ｆｉｌｔｅｒはフィルタ係数である。フィルタ係数Ｃ＿ｆｉｌｔｅｒは、フィルタ長が（２×ｓ（ｍ）＋１）である。フィルタ係数Ｃ＿ｆｉｌｔｅｒの形状は、例えば、ハニング窓およびＳｉｎｃ関数などであるが、用途に応じて選定することができる。ｓ（ｍ）は、次の式（８）で表される。
ｓ（ｍ）＝四捨五入（Ｄｒｅｌａｔｉｖｅ（ｍ）×Ｃ÷３６０×Ｍ）  ・・・（８）
【００３９】
式（８）において、Ｄｒｅｌａｔｉｖｅ（ｍ）は、方位ｍにおけるアレイ基準方位からの相対方位の半減半角ビーム幅である。Ｃは調整係数である。調整係数Ｃは、正の値であり、全ての方位ｍに共通の値である。Ｍは、全周の方位分割数であり、信号強度分布データＸｔ（ｍ）の長さに相当する。
【００４０】
次に、本実施の形態１の信号追尾装置１の動作を説明する。図８は、図１に示した信号追尾装置が実行する信号追尾方法の手順を示すフローチャートである。図８に示すステップＳ１０１およびＳ１０３〜Ｓ１０６の処理は図６を参照して説明したステップＳ２０１〜Ｓ２０５の処理と同様なため、ここでは、その詳細な説明を省略する。
【００４１】
ステップＳ１０１において、データ入力手段２１は、入力される信号強度分布データＸｔ（ｍ）を不要ピーク除去手段２２に出力する。ｍ＝１，２，・・・，Ｍである。ステップＳ１０２において、不要ピーク除去手段２２は、方位毎の信号強度分布データＸｔ（ｍ）に対して、式（７）にしたがって、フィルタ係数Ｃ＿ｆｉｌｔｅｒによるローパスフィルタを実行し、受け取ったデータから不要ピークを除去する。不要ピーク除去手段２２は、不要ピーク除去後の信号強度分布データＸｆｉｌｔｅｒｔ（ｍ）をピーク検出手段２３に出力する。
【００４２】
ここで、信号処理の一例として、図７に示した、適応型広帯域分析処理により生成された信号強度分布データＸｔ（ｍ）を、図８に示した信号追尾方法に適用した場合を説明する。図９は、本発明の実施の形態１に係る信号追尾方法において、不要ピーク除去前の信号強度分布データＸｔ（ｍ）の一例を示す図である。図９の横軸は方位ｍであり、縦軸は信号強度を示すレベルである。
【００４３】
図９に示すように、各方位ｍにおける相対方位に応じてフィルタ長が異なっている。不要ピーク除去手段２２は、式（７）にしたがって、方位ｍ毎に、近傍の方位の信号強度とフィルタ係数との乗算結果の和を算出することで、不要ピークを除去した信号強度分布データＸｆｉｌｔｅｒｔ（ｍ）を得る。
【００４４】
図１０は、図９に示した信号強度分布データから不要ピークを除去した後の信号強度分布データの一例を示す図である。図１０の横軸は方位ｍであり、縦軸は信号強度を示すレベルである。図１０を図９と比較すると、不要なピークが除去されていることがわかる。
【００４５】
本実施の形態１の信号追尾装置１は、広帯域分析処理により生成される信号強度分布データに対して、方位毎にアレイ相対方位の半減半角ビーム幅に基づくフィルタ長のフィルタ係数を用いて、不要ピークを除去する不要ピーク除去手段２２を有する。
【００４６】
本実施の形態１によれば、広帯域分析処理により生成された信号強度分布データに対して、不要ピーク除去手段２２が方位方向のローパスフィルタを実行することにより、方位全体における散在する雑音となるピーク点が除去される。入力された信号強度分布データ中に散在している不要なピーク点が除去され、本来の信号成分のピーク点のみを残すことができる。その結果、雑音によるピークの誤検出が抑制され、誤警報の回数を抑制することができる。
【００４７】
また、本実施の形態１によれば、適応型広帯域分析処理により生成された信号強度分布データが入力される場合でも、方位全体における散在する雑音となるピーク点を除去することができる。そのため、誤警報の発生が抑制されるだけでなく、従来の広帯域分析処理によるデータが入力される場合と比べて、異なる音源から発せられる信号同士の分離能力が向上した状態で広帯域信号成分をラインとして検出および追尾することができる。
【符号の説明】
【００４８】
１  信号追尾装置
１１  制御部
１２  表示部
１３  制御部
２１  データ入力手段
２２  不要ピーク除去手段
２３  ピーク検出手段
２４  方位追尾手段
２５  検出判定手段
２６  表示制御手段
５１  ＣＰＵ
５２  メモリ
１００  信号追尾装置
【図１】
2020085525.tif 000006
【図２】
2020085525.tif 000007
【図３】
2020085525.tif 000008
【図４】
2020085525.tif 000009
【図５】
2020085525.tif 000010
【図６】
2020085525.tif 000011
【図７】
2020085525.tif 000012
【図８】
2020085525.tif 000013
【図９】
2020085525.tif 000014
【図１０】
2020085525.tif 000015