(57)【要約】【課題】  追尾フィルタに入力する追尾信号を有効に活用し、また各センサに対する時刻毎のそれぞれの潮流の影響を考慮して目標状態量の推定を精度良く行う。【解決手段】  追尾信号として実際には存在しない遅延時刻（ｔ−１）から現在時刻ｔまでの追尾信号、すなわち周波数情報ｆ３（ｔ−１），ｆ３（ｔ），ｆ２（ｔ）を、予測追尾信号外挿処理６において目標の状態量から予測し、この予測した周波数情報を新たな周波数情報として追尾フィルタ３に供給することにより、現在時刻ｔにおける周波数情報を得る。また、浮遊するＧＰＳ搭載のセンサから時々刻々の各センサ位置を正確に入手し、センサ毎、時刻毎の潮流情報ｕ３（ｔ−１），ｕ３（ｔ），ｕ２（ｔ）を周波数情報に加えて追尾信号として追尾処理を行う。時間遅延に対して追尾信号を予測する際、目標の状態量に加え、得られた各センサの潮流情報を含めて、周波数情報を算出する。【選択図】      図２
2008008691.tif 000002
【特許請求の範囲】
【請求項１】
出力間に互いに遅延時間差が生じる複数のセンサによって、目標の運動により発せられる信号を受波して得られる観測追尾信号を追尾フィルタに入力して、目標状態量を出力するようにした追尾処理装置であって、
前記観測追尾信号のうち実際に存在する遅延時刻の信号を利用して、実際には存在しない前記遅延時刻の信号より新しい時刻の信号を予測する予測手段と、
この予測追尾信号を新たな追尾信号として前記追尾フィルタへ入力する入力手段とを含むことを特徴とする追尾処理装置。
【請求項２】
前記観測追尾信号は、前記センサの各出力を周波数分析することにより得られた周波数情報であり、前記センサの各々の時々刻々の位置から得られた潮流情報を前記追尾信号に付与するようにしたことを特徴とする請求項１記載の追尾処理装置。
【請求項３】
前記入力手段は、遅延時刻が同じ前記観測追尾信号及び前記予測追尾信号に対して前記潮流情報を付与して追尾信号セットとして生成して前記追尾フィルタに入力するよう構成されており、
前記予測手段は、前記追尾フィルタから出力された目標状態量を用いて、観測追尾信号として実際には存在しない時刻の追尾信号を予測して前記予測追尾信号として前記第二の手段へ渡すよう構成されており、
前記入力手段及び前記予測手段に対して、一番古い時刻の信号から現在時刻の信号まで順次繰り返し処理するよう制御する制御手段を、更に含むことを特徴とする請求項１または２記載の追尾処理装置。
【請求項４】
前記入力手段は、遅延時刻が同じ前記観測追尾信号及び前記予測追尾信号に対して前記潮流情報を付与して追尾信号セットとして生成して前記追尾フィルタに入力するよう構成されており、
前記予測手段は、前記追尾フィルタから出力された時刻における目標状態量を用いて、観測追尾信号として実際には存在しない時刻の追尾信号を予測して前記予測追尾信号として前記入力手段へ渡すよう構成されており、
前記入力手段及び前記予測手段に対して、一番古い時刻の信号から現在時刻の直前の時刻の信号まで順次繰り返し処理するよう制御する制御手段と、
前記直前の時刻の信号を入力としたときの前記追尾フィルタから出力された前記目標状態量に対して、時間遷移処理を行って現在時刻の状態目標量を算出する時間遷移処理手段とを、更に含むことを特徴とする請求項１または２記載の追尾処理装置。
【請求項５】
前記追尾信号セットが前記観測追尾信号のみからなる場合に、前記追尾フィルタの追尾パラメータを格納する格納手段を更に含み、前記追尾フィルタの前記現在時刻に続く追尾パラメータとして、前記格納手段に格納されたパラメータを用いることを特徴とする請求項３または４記載の追尾処理装置。
【請求項６】
出力間に互いに遅延時間差が生じる複数のセンサによって、目標の運動により発せられる信号を受波して得られる観測追尾信号を追尾フィルタに入力して、目標状態量を出力するようにした追尾処理方法であって、
前記観測追尾信号のうち実際に存在する遅延時刻の信号を利用して、実際には存在しない前記遅延時刻の信号より新しい時刻の信号を予測する予測ステップと、
この予測追尾信号を新たな追尾信号として前記追尾フィルタへ入力する入力ステップとを含むことを特徴とする追尾処理方法。
【請求項７】
前記観測追尾信号は、前記センサの各出力を周波数分析することにより得られた周波数情報であり、前記センサの各々の時々刻々の位置から得られた潮流情報を前記追尾信号に付与するようにしたことを特徴とする請求項６記載の追尾処理方法。
【請求項８】
前記入力ステップは、遅延時刻が同じ前記観測追尾信号及び前記予測追尾信号に対して前記潮流情報を付与して追尾信号セットとして生成して前記追尾フィルタに入力するステップを有し、
前記予測ステップは、前記追尾フィルタから出力された目標状態量を用いて、観測追尾信号として実際には存在しない時刻の追尾信号を予測して前記予測追尾信号として前記入力ステップへ渡すステップを有しており、
前記入力ステップ及び前記予測ステップに対して、一番古い時刻の信号から現在時刻の信号まで順次繰り返し処理するよう制御する制御ステップを、更に含むことを特徴とする請求項６または７記載の追尾処理方法。
【請求項９】
前記入力ステップは、遅延時刻が同じ前記観測追尾信号及び前記予測追尾信号に対して前記潮流情報を付与して追尾信号セットとして生成して前記追尾フィルタに入力するステップを有し、
前記予測ステップは、前記追尾フィルタから出力された時刻における目標状態量を用いて、観測追尾信号として実際には存在しない時刻の追尾信号を予測して前記予測追尾信号として前記入力ステップへ渡すステップを有し
前記入力ステップ及び前記予測ステップに対して、一番古い時刻の信号から現在時刻の直前の時刻の信号まで順次繰り返し処理するよう制御する制御ステップと、
前記直前の時刻の信号を入力としたときの前記追尾フィルタから出力された前記目標状態量に対して、時間遷移処理を行って現在時刻の状態目標量を算出する時間遷移処理ステップとを、更に含むことを特徴とする請求項６または７記載の追尾処理方法。
【請求項１０】
前記追尾信号セットが前記観測追尾信号のみからなる場合に、前記追尾フィルタの追尾パラメータを格納する格納手段を設けておき、前記追尾フィルタの前記現在時刻に続く追尾パラメータとして、前記格納手段に格納されたパラメータを用いることを特徴とする請求項８または９記載の追尾処理方法。
【請求項１１】
出力間に互いに遅延時間差が生じる複数のセンサによって、目標の運動により発せられる信号を受波して得られる観測追尾信号を追尾フィルタに入力して、目標状態量を出力するようにした追尾処理方法をコンピュータにより実行させるためのプログラムであって、
前記観測追尾信号のうち実際に存在する遅延時刻の信号を利用して、実際には存在しない前記遅延時刻の信号より新しい時刻の信号を予測する処理と、
この予測追尾信号を新たな追尾信号として前記追尾フィルタへ入力する処理とを含むことを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は追尾処理装置及びその方法並びにプログラムに関し、特に出力間に互いに遅延時間差が生じる複数のセンサによって、目標の運動により発せられる信号を受波して得られる観測追尾信号を追尾フィルタに入力して、目標状態量を出力するようにした追尾処理方式に関するものである。
【背景技術】
【０００２】
目標のリアルタイムな追尾において、目標の運動により発せられる信号の周波数情報を追尾信号とする場合、そのデータに時間遅延が発生する。周波数情報は従来からＦＦＴ（Fast Fourier Transform）処理により算出するのが一般的である。しかしながら、例えば、０．１Ｈｚの周波数分解能を得ようとする場合、必要なデータ時間長は１０秒となり、その１／２である５秒の時間遅延が見込まれる。つまり、時刻ｔに出力された周波数情報は時刻（ｔ−５秒）の情報にもかかわらず、時刻ｔにおける目標推定位置を出力していることになる。
【０００３】
これに関して、複数のセンサの各々から得られた信号が同じ周波数分解能で出力されれば信号の時刻は整合されているので、問題ない。しかし、周波数分解能が異なるものがある場合（例えば、０．１Ｈｚと０．２Ｈｚ）、周波数は５秒前と２．５秒前（０．２Ｈｚの時のデータ時間長は５秒）の信号になり、時間整合されないのである。よって、各センサの出力信号である追尾信号をそのまま使用すれば、追尾精度が悪くなる。
【０００４】
この問題を回避するために、従来、図９に示すように、各センサから得られた追尾信号のうち、時刻ｔに対し、時間遅延があれば一番古い時刻の追尾信号に合わせて追尾信号セットとし、追尾フィルタに入力することにより、目標の状態量（位置・速度等）を推定している。そして、時刻（ｔ−２）の状態量から時刻ｔの状態を算出するために、遅延時間分目標位置を遷移させることにより、現在時刻における目標の状態量を出力している。
【０００５】
図９を用いてより詳述する。図９は従来技術を説明するためのイメージ図である。図９を参照すると、複数のセンサ＃１〜＃３（本例では、３個とする）によるセンサ出力は信号処理部１０１へ入力されて各センサの観測追尾信号が生成されて出力される。各センサの追尾信号は、図示する如く、目標の運動により発せられる周波数情報ｆｎ（ｔ）と、各センサが受ける潮流情報ｕｎ（ｔ）とからなる。なお、ｎはセンサの番号を示しており、図９の例では、ｎ＝１〜３であり、ｔは現在時刻を示す。また、潮流情報ｕｎ（ｔ）は、１時刻前と現在取得したセンサ位置から１時刻分のセンサの移動量である。
【０００６】
図９の例では、各センサ＃１〜＃３の周波数分解能に差があり、時刻ｔでは、センサ＃１による観測追尾信号はｆ１（ｔ），ｕ１（ｔ）であり、センサ＃２による観測追尾信号はｆ２（ｔ−１），ｕ２（ｔ−１）であり、センサ＃３による観測追尾信号はｆ３（ｔ−２），ｕ３（ｔ−２）である。また、時刻ｔ−１及び時刻ｔ−２における各センサによる観測追尾信号はそれぞれ図示の如くなる。
【０００７】
このような状態において、時刻ｔに対して、一番古い時刻ｔ−２の追尾信号に合わせて、追尾信号セットとして、追尾フィルタ３へ入力するようにし、この追尾フィルタ３において、Ｆに示す如く、時刻ｔ−２の目標の状態量Ｘ（ｔ−２）を推定している。そして、時刻ｔ−２の状態量から時刻ｔの状態量Ｘ（ｔ）を算出するために、Ｇに示すように、遅延時間分目標位置を遷移させることにより、現在時刻における目標の状態量（位置及び速度など）を出力するようになっている。
【０００８】
また、海中等に浮遊するセンサを用いた場合、従来では、各センサ位置は固定として潮流の影響を考慮しないか、または全てのセンサに同一の潮流の影響があるとして追尾処理を行っている（図９において、ｕ１＝ｕ２＝ｕ３）。しかし、実際には、各センサは独立した潮流の影響が存在するために、各センサに対して、各時刻毎に、それぞれ異なる潮流の影響があって、追尾精度は悪くなる。
【０００９】
なお、関連技術としては、特許文献１〜５があげられる。
【特許文献１】特開２０００−２４１５２２号公報
【特許文献２】特開２００２−２６７７４５号公報
【特許文献３】特開２００４−０３７２６２号公報
【特許文献４】特開平８−３３４２９８号公報
【特許文献５】特開平１０−０６２５１１号公報
【発明の開示】
【発明が解決しようとする課題】
【００１０】
従来技術の第一の問題点は、各追尾信号間に時間遅延差がある場合、一番古い時刻の追尾信号セットを入力とするために、現在時刻に近い情報は使用されず、追尾精度が落ちる点である。その理由は、追尾信号の時刻整合を図るために、一番古い情報のみに頼ることとなり、精度の良い新しい情報が得られている場合でもその信号が有効利用できないからである。
【００１１】
また、第二の問題点は、水中目標の追尾において、全てのセンサ・時刻に同一の潮流の影響があるとして追尾処理を行う場合、追尾精度が落ちる点である。その理由は、実際には各センサに対して各時刻毎に、それぞれに独立した潮流の影響があるにもかかわらず、全てのセンサで同一時刻に同一の潮流の影響があるとして追尾処理を行うからである。
【００１２】
本発明の目的は、追尾フィルタに入力する追尾信号を有効に活用すること及び各センサに対する時刻毎のそれぞれの潮流の影響を考慮することにより、目標状態量（位置・速度等）の推定を精度良く行うことが可能な追尾処理装置及びその方法並びにプログラムを提供することである。
【課題を解決するための手段】
【００１３】
本発明による追尾処理装置は、出力間に互いに遅延時間差が生じる複数のセンサによって、目標の運動により発せられる信号を受波して得られる観測追尾信号を追尾フィルタに入力して、目標状態量を出力するようにした追尾処理装置であって、前記観測追尾信号のうち実際に存在する遅延時刻の信号を利用して、実際には存在しない前記遅延時刻の信号より新しい時刻の信号を予測する予測手段と、この予測追尾信号を新たな追尾信号として前記追尾フィルタへ入力する入力手段とを含むことを特徴とする。
【００１４】
本発明による追尾処理方法は、出力間に互いに遅延時間差が生じる複数のセンサによって、目標の運動により発せられる信号を受波して得られる観測追尾信号を追尾フィルタに入力して、目標状態量を出力するようにした追尾処理方法であって、前記観測追尾信号のうち実際に存在する遅延時刻の信号を利用して、実際には存在しない前記遅延時刻の信号より新しい時刻の信号を予測する予測ステップと、この予測追尾信号を新たな追尾信号として前記追尾フィルタへ入力する入力ステップとを含むことを特徴とする。
【００１５】
本発明によるプログラムは、出力間に互いに遅延時間差が生じる複数のセンサによって、目標の運動により発せられる信号を受波して得られる観測追尾信号を追尾フィルタに入力して、目標状態量を出力するようにした追尾処理方法をコンピュータにより実行させるためのプログラムであって、前記観測追尾信号のうち実際に存在する遅延時刻の信号を利用して、実際には存在しない前記遅延時刻の信号より新しい時刻の信号を予測する処理と、この予測追尾信号を新たな追尾信号として前記追尾フィルタへ入力する処理とを含むことを特徴とする。
【００１６】
本発明の作用を述べる。追尾信号として実際には存在しない遅延時刻（ｔ−１）から現在時刻ｔまでの追尾信号、すなわち周波数情報ｆ３（ｔ−１），ｆ３（ｔ），ｆ２（ｔ）を、予測追尾信号外挿処理において目標の状態量から予測し、この予測した周波数情報を新たな周波数情報として追尾フィルタに供給することにより、現在時刻ｔにおける周波数情報を得る。また、浮遊するＧＰＳ搭載のセンサから時々刻々の各センサ位置を正確に入手し、センサ毎、時刻毎の潮流情報ｕ３（ｔ−１），ｕ３（ｔ），ｕ２（ｔ）を周波数情報に加えて追尾信号として追尾処理を行う。時間遅延に対して追尾信号を予測する際、目標の状態量に加え、得られた各センサの潮流情報を含めて、周波数情報を算出する。
【発明の効果】
【００１７】
本発明による第１の効果は、従来利用されていなかった最新情報（最も古い時刻以降にある情報）を利用し、現在時刻での目標状態量の推定精度を向上できることである。その理由は、時間遅延が大きく、情報として古い追尾信号を利用して、遅延時刻から現在までの追尾信号を予測し、新たな追尾信号を生成して追尾フィルタに入力することにより、より正確な状態量を推定できるためである。
【００１８】
本発明による第２の効果は、従来、ブイ位置は固定又はブイ間の相対位置は固定として推定されていた追尾状態を、潮流状態も含めて推定し、現在時刻での目標状態量の推定精度を向上できることである。その理由は、センサ毎の時々刻々の潮流情報を付与して追尾処理を行うこと、及び潮流状態を考慮して遅延時間分の予測追尾信号を生成し、追尾フィルタに入力することにより、より正確な状態量を推定できるためである。
【発明を実施するための最良の形態】
【００１９】
以下に、図面を参照して本発明の実施例について説明する。図１は本発明の実施例の機能ブロック図である。図１を参照すると、本発明の実施例は、図示せぬセンサから得られる観測データを信号処理し追尾信号を生成すると同時に、ＧＰＳから得られる位置情報から潮流情報を算出する信号入力部１と、生成された各センサの追尾信号を周波数分解能等の違いにより発生する遅延時間に合わせて潮流情報を付与し、追尾信号セットを生成する時間整合処理部２と、入力された追尾信号から追尾結果を算出する追尾フィルタ３と、生成された追尾信号セットが観測データのみから得られたものかどうかを判断する全観測追尾信号判断部４とを含んでいる。
【００２０】
更に、本発明の実施例は、追尾フィルタ３から出力された追尾結果に対し、現在時刻から遅延しているかどうかを判断する現在時刻判断部５と、遅延している場合、実際には存在しない遅延時刻から１時刻進めた追尾信号を潮流情報を含めて予測し、予測した追尾信号を新たな追尾信号として入力する予測追尾信号外挿処理部６と、観測データのみから生成された追尾信号セットにより算出された追尾フィルタ出力（追尾パラメータ）を格納する追尾パラメータ格納部７と、追尾結果を出力する追尾結果出力部８と、これら各部を制御する制御部（ＣＰＵ）９と、この制御部の制御動作をプログラムとして格納したメモリ１０とを含んで構成されている。
【００２１】
図２は本発明の動作イメージ図であり、図１及び図２を参照すると、信号入力部１は、複数のセンサ＃１〜＃ｎ（ｎはセンサ番号）が受波した信号を、それぞれ周波数分析することにより、観測周波数情報ｆ１〜ｆｎを取得する。また、各センサが受ける潮流情報を算出しｕ１〜ｕｎを取得する。潮流情報の算出は、１時刻前と現在取得したセンサ位置から１時刻分のセンサの移動量を求める。もし、短い時間内で潮流情報に大きなばらつきがある場合、積分等を行い平滑化を行う。
【００２２】
時間整合処理部２は、取得した各センサの信号（以下、観測追尾信号と呼ぶ）に対して先ず音波伝播時間分の補正を行い、次に周波数分解能の違いにより発生する遅延時間に合わせて潮流情報を付与し、遅延時刻が揃っている信号群を１セットとして追尾フィルタ３に入力する。追尾フィルタ３は、例えば拡張カルマンフィルタや最小二乗フィルタを想定する。
【００２３】
追尾フィルタ３から出力された目標の推定状態量をＸとする。一般的に、Ｘは成分ｘ（目標位置のｘ軸成分）、ｙ（目標位置のｙ軸成分）、ｖｘ（目標速度のｘ軸成分）、ｖｙ（目標速度のｙ軸成分）、ｆ０（目標の運動により発された真の周波数）を含むものである。また、追尾パラメータとして目標状態量及び入力信号の分散情報も合わせて出力される。以下、この目標状態量Ｘと追尾パラメータを用いて説明することにする。
【００２４】
追尾フィルタ３の出力は全観測追尾信号判断部４に入力される。この判断部４は状態量Ｘが全て観測追尾信号のみからなる追尾信号セットにより得られたものかどうかを判断し、そうであれば、状態量Ｘと追尾パラメータとを一時的に追尾パラメータ格納部７に格納して現在時刻判断部５へ処理を渡す。そうでなければ、追尾パラメータは追尾パラメータ格納部７に格納することなく現在時刻判断部５へ処理移行する。
【００２５】
現在時刻判断部５は、追尾フィルタ出力が現在時刻であるかどうかを判断する。現在時刻ｔであれば、追尾フィルタ出力を追尾結果として出力し、次の追尾信号入力（ｔ＋１）を待つ。時刻（ｔ＋１）の追尾信号が入力された場合、同処理を繰り返す。一方、追尾フィルタ出力が現在時刻ｔではない場合、追尾フィルタ出力を予測追尾信号外挿処理部６へ入力する。
【００２６】
予測追尾信号外挿処理部６は、追尾フィルタ出力から予測追尾信号を算出する。図２の本発明の動作イメージ図に示すように、時刻ｔで得られた最古の追尾信号の周波数をｆｎ（ｔ−２）とすると、１時刻進めたｆｎ（ｔ−１）を以下の式（１）及び式（２）によって算出する。
【００２７】
【数１】
2008008691.tif 000003
【００２８】
【数２】
2008008691.tif 000004
【００２９】
但し、例えばセンサ＃１、センサ＃２のＳ／Ｎが悪く、これに対してセンサ＃３のＳ／Ｎが非常に良かった場合等全ての信号を用いて追尾信号を予測するよりも、センサ＃３の情報のみを用いて予測追尾信号を算出する方が精度が高いと判断される場合は、センサ３の情報のみから予測追尾信号を算出する。これは、例えば最小２乗法等を用いて算出する。ここでこの判断は、各センサが得ている信号の観測誤差分散及び目標の状態量に対する誤差分散等から判断する。
【００３０】
次に、算出されたｆｎ（ｔ−１）は新たに追尾信号として時間整合処理部２へ入力される。時間整合処理部２は時刻（ｔ−１）の追尾信号セットを生成し追尾フィルタ部３に入力する。これを現在時刻ｔまで繰り返す。追尾フィルタ出力が現在時刻ｔまで達したら、現在時刻ｔの追尾結果として出力し時刻（ｔ＋１）での観測追尾信号入力を待つ。
【００３１】
ここで、予測追尾信号外挿処理を行った場合、時刻（ｔ＋１）での追尾フィルタに使用する追尾パラメータは、繰り返しで得られた追尾結果に付随したパラメータではなく、追尾パラメータ格納部７に格納された追尾パラメータ（観測追尾信号セット：時刻（ｔ−２））を用いる。これにより、観測追尾信号のみで計算された信頼性の高い追尾パラメータを用いることになるため、もし予測追尾信号外挿処理によって得られた予測追尾信号の精度が悪い場合でも、この結果を引きつぐことがないので、従来の追尾性能を下回ることなく追尾が可能となる。
【００３２】
以下に、センサが＃１〜＃３の３個の場合について、図２を参照して説明する。信号入力部１は図示せぬ複数のセンサから得られた観測データ（観測追尾信号）から周波数情報ｆｎ（ｔ）及び潮流情報ｕｎ（ｔ）を算出する。各センサにはＧＰＳが搭載されているものとし、各時刻毎にセンサ位置を取得できるようになっているものとする。
【００３３】
追尾信号情報としては、周波数分析によって得られる周波数情報や、積分のかかった方位情報等を想定する。また潮流情報は１時刻前と現在取得したセンサ位置から１時刻分のセンサの移動量を求めることにより、算出される。もし、短い時間内で潮流情報に大きなバラツキがあれば、積分などを行って平滑化を行う。なお、センサ＃１は信号処理による時刻遅延なし、センサ＃２は（ｔ−１）の時刻遅延、センサ＃３は（ｔ−２）の時刻遅延を有するものとする。
【００３４】
時間整合処理部２は、時刻ｔに３つの観測追尾信号が入力されると、観測追尾信号のうち一番古いデータ時刻（ｔ−２）に合わせた追尾信号セットを図２に示す様に生成し、追尾フィルタ３へ供給する。追尾フィルタ３は、Ａに示すように、時刻（ｔ−２）における目標状態量Ｘと追尾パラメータとを出力する。
【００３５】
追尾フィルタ３として、拡張カルマンフィルタ（非線形フィルタ）とすると、その内部パラメータ（追尾パラメータ）としては、状態量を示す行列Ｘ（出力である位置情報ｘ，ｙ、速度情報ｖｘ，ｖｙ、目標の運動により発せられる真の周波数ｆ０で示される）、この状態量の誤差分散を示す行列、及び観測値の誤差分散を示す行列（入力信号分散量）がある。これらの値をフィルタ３が保持しており、通常は、観測値が入力される毎に、これらパラメータ値は更新されるようになっている。
【００３６】
全観測追尾信号判断部４は、この追尾フィルタ出力が観測追尾信号のみからなるセット（予測追尾信号外挿処理部６により予測された予測追尾信号を含まないセット）により出力されていると判断した場合、追尾パラメータ格納部７に追尾パラメータを格納する。現在時刻判断部５は、追尾フィルタ出力が現在時刻ｔのものではないことを判断し、予測追尾信号外挿処理部６に追尾フィルタ出力を供給する。この予測追尾信号外挿処理部６は、１時刻だけ進めたセンサの（ｔ−２）の予測追尾信号を算出し、時間整合処理部２へ供給する。
【００３７】
時間整合処理部２はデータ時刻（ｔ−１）に合わせた追尾信号セットを追尾フィルタ３へ供給する。追尾フィルタ３はＢに示す時刻（ｔ−１）における目標状態量Ｘと追尾パラメータを出力する。全観測追尾信号判断部４は、追尾フィルタ出力が予測追尾信号を含んだデータセットにより出力されているものと判断し、当該出力を現在時刻判断部５へ供給する。この判断部５は追尾フィルタ出力が現在時刻ｔのものではないと判断し、予測追尾信号外挿処理部６に追尾フィルタ出力を供給する。
【００３８】
予測追尾信号外挿処理部６は、更に１時刻だけ進めたセンサ＃２，＃３の時刻ｔの予測追尾信号を算出し、時間整合処理部２へ供給する。時間整合処理部２はデータ時刻ｔに合わせた追尾信号セットを追尾フィルタ３へ供給する。追尾フィルタ３はＣに示す時刻ｔにおける目標状態量Ｘと追尾パラメータを出力する。全観測追尾信号判断部４は、追尾フィルタ出力が予測追尾信号を含んだセットにより出力されているものと判断し、現在時刻判断部５へこの出力を供給する。この判断部５は、追尾フィルタ出力が現在時刻ｔのものであることを判断し、これを最終的な時刻ｔにおける追尾結果として、追尾結果出力部８を介して出力する。
【００３９】
時刻（ｔ＋１）に新たな観測追尾信号が入力されると、追尾パラメータ格納部７に格納されている追尾パラメータを用いて上述の処理を繰返し、時々刻々の追尾結果を出力するようになっている。
【００４０】
次に、本発明の実施例の動作の詳細を図３のフローチャートを参照しつつ説明する。時刻ｔに、信号入力部１は各センサから得た観測データに対して、指定されたモードで周波数分析を行い周波数情報を算出する。また、同時に各センサ位置を取得し、潮流情報を算出する。算出された結果は、時間整合処理部２に入力される（ステップＳ１）。
【００４１】
ここで、各観測値に対して異なる周波数分析モードが設定されていた場合、時間整合処理部２には各センサ間で時間遅延のある観測追尾信号が入力される。この場合、時間整合処理部２は観測追尾信号のうち一番古いデータ時刻（例としてｔ−２）に合わせた追尾信号セットを生成し、追尾フィルタ３に入力する（ステップＳ２）。
【００４２】
追尾フィルタ３は時刻（ｔ−２）における目標状態量Ｘと追尾パラメータを出力する（ステップＳ３）。ここで、判断部４はこの追尾フィルタ出力が観測追尾信号セットにより出力されているものと判断し（ステップＳ４）、追尾パラメータ格納部７に追尾パラメータを保存する（ステップＳ５）。次に、判断部５は追尾フィルタ出力が現在時刻ｔのものではないことを判断し、予測追尾信号外挿処理部６に追尾フィルタ出力を入力する（ステップＳ６）。予測追尾信号外挿処理部６は１時刻だけ進めたｔ−１の予測追尾信号（周波数情報）を算出する。算出された予測追尾信号は、時間整合処理部２へ入力される（ステップＳ７）。
【００４３】
時間整合処理部２は、算出された予測追尾信号と時刻ｔまでに存在する観測追尾信号により、新たに時刻（ｔ−１）の追尾信号セットを生成し、追尾フィルタ３に入力する。ここまでの処理を繰り返し、現在時刻ｔにおける追尾信号セットを追尾フィルタ部３に入力し、出力された結果が最終的な時刻ｔにおける追尾結果として出力される（ステップＳ８）。時刻（ｔ＋１）に新たな追尾信号（観測追尾信号）が入力されたら、追尾パラメータ格納部７に保存しておいた追尾パラメータを用いて、ここまでの処理を繰り返し、時々刻々の追尾結果を出力する。
【００４４】
ここで、本発明手法を用いてシミュレーションを行った結果を図４及び図５のグラフに示す。図４のグラフは従来手法と本手法との比較を行ったものである。シミュレーション方法としては、センサ３本を用いて目標を追尾する場面を想定し、各センサが得る追尾信号として、模擬信号（方位情報及び周波数情報）を作成し、この信号に対して追尾処理を行った。ここで、センサ＃１〜＃３の周波数情報はそれぞれ分解能が異なるものとした。また信号のＳ／Ｎは比較的良いものと仮定した。
【００４５】
追尾処理の手法としては、補正を行わない（時間補正しない）追尾方式、従来の追尾方式及び本発明手法による追尾方式を用いた。この３種類の追尾処理結果（位置情報）と正解位置とのズレ（＝距離誤差）を測定することにより、評価を行った。なお、追尾フィルタは拡張カルマンフィルタを用いた。この結果より、補正がないものよりも補正を行っている２手法の結果の方が誤差が小さいことが分かる。また、従来手法と本発明手法を比較した場合、本発明手法による追尾結果は従来手法に比べて全体的に誤差が小さいことが分かる。
【００４６】
特徴としては、従来手法では急に誤差が大きくなる箇所があるのに対し、本発明ではこの箇所が平滑化され良好な結果となっている。これは、一時的にセンサに観測誤差が大きい信号が入ってきた場合でも、遅延時間分予測追尾信号を作成し、従来方式よりも多く追尾フィルタを回すことによって目標状態量の誤差分散を小さくし、追尾誤差を小さくしたものと考えられる。
【００４７】
図５は、Ｓ／Ｎが悪い状況下でも本発明手法が有効であるかを評価したものである。Ｓ／Ｎが悪い状況の場合、追尾が出来る最低限（２本）のセンサ以外のセンサ信号（３本目以降）を用いて追尾すると、３本目以降のセンサの誤差が起因し、追尾結果が悪くなることが考えられる。特に、本発明では、新たに信号を予測して追尾を行うために、予測追尾信号の誤差が正解値から大きく外れてしまうと、発散してしまうことも考えられる。
【００４８】
そこで、Ｓ／Ｎの悪い状況下で、２本のセンサで追尾を行う場合と、３本のセンサにより本発明手法を用いて追尾処理を行う場合を比較した。ここで、３本目のセンサは２本のセンサよりも更にＳ／Ｎの悪い信号を与えてシミュレーションを行った。この結果より、Ｓ／Ｎが悪い３本目のセンサを用いた場合でも、２本のセンサで追尾するよりも誤差が小さいことがわかる。
【００４９】
これは、たとえ観測誤差分散が大きな値をとったとしても、追尾フィルタにより状態量は収束していくために、比較的良好な予測追尾信号を作成することができ、且つ遅延時間分従来よりも多く追尾フィルタを回すことにより、目標状態量の分散を小さすることが出来るためと考えられる。特に、観測誤差分散値が大きくても適切なカルマンゲインを算出できる拡張カルマンフィルタ等を追尾フィルタに用いる場合に、より効果的であるといえる。以上の結果より、本発明手法の有効性が確認できることになる。
【００５０】
図６は本発明の他の実施例の構成を示すブロック図であり、図１と同等部分は同一符号により示している。本例では、図１の構成に、時間遷移処理部１１を追加したものであり、他は図１と同一である。
【００５１】
図６を参照すると、時間整合処理部２に現在時刻ｔに得られた観測追尾信号のうち、現在時刻ｔのデータがない場合（図７のイメージ図に示す関係にある場合）、先の実施例と同様な処理を行い、最新のデータ時刻（ｔ−１）まで予測追尾信号を外挿し、現在時刻判断部５に出力する。この判断部５は、追尾フィルタ出力が最新時刻のものであるかどうかを判断し、最新時刻（ｔ−１）であれば、追尾フィルタ出力結果を時間遷移処理部１１に入力する。この時間遷移処理部１１は時刻（ｔ−１）から現在時刻ｔまで時間遷移処理を行うことにより、時刻ｔにおける追尾結果を出力する。時間遷移は式（１）を用いる。
【００５２】
次に、海中に浮遊するセンサ三本で受波があった場合の動作について、図８を用いて説明する。信号入力部１は、先の実施例と同様な処理を行うことにより、観測追尾信号及び潮流情報を算出する（ステップＳ１）。ここで、センサ＃１は信号処理により（ｔ−１）の時刻遅延、センサ＃２は（ｔ−２）の時刻遅延、センサ＃３は（ｔ−３）の時刻遅延をもつものとする。
【００５３】
いま、時刻ｔに、３本分の信号が観測追尾信号として時間整合処理部２に入力されると、時間整合処理部２は観測追尾信号のうち一番古いデータ時刻（ｔ−３）に合わせた追尾信号セットを生成し、追尾フィルタ３に入力する（ステップＳ２）。追尾フィルタ３は時刻（ｔ−３）における目標状態量Ｘ（図７のＥ）と追尾パラメータを出力する（ステップＳ３）。
【００５４】
全観測追尾信号判断部４は、この追尾フィルタ出力が観測追尾信号セットにより出力されているものと判断し、追尾パラメータ格納部７に追尾パラメータを保存する（ステップＳ４，Ｓ５）。次に、現在時刻判断部５は、追尾フィルタ出力が最新データ時刻（ｔ−１）のものではないと判断し（ステップＳ６）、予測追尾信号外挿処理部６に追尾フィルタ出力を入力する。予測追尾信号外挿処理部６は１時刻だけ進めたセンサ３の時刻（ｔ−２）の予測追尾信号を算出し、時間整合処理部２に出力する（ステップＳ７）。
【００５５】
次に、時間整合処理部２はデータ時刻（ｔ−２）に合わせた追尾信号セットを追尾フィルタ３に入力する。追尾フィルタ３は時刻（ｔ−２）における目標状態量Ｘ（図７のＡ）と追尾パラメータを出力する。
【００５６】
この後、データ時刻（ｔ−１）まで、先の実施例と同様の繰り返し処理を行う。最新データ時刻（ｔ−１）に達したら、判断部５は追尾フィルタ出力が最新データ時刻（ｔ−１）のものであることを判断し（図７のＢ）、時間遷移処理部１１に入力する。時間遷移処理部１１は時間遷移処理により、図７のＤに示すように、１刻分進めた時刻ｔの追尾結果を算出し（ステップＳ９）、これを最終的な現在時刻ｔにおける追尾結果として出力する（ステップＳ８）。
【００５７】
時刻（ｔ＋１）に新たな観測追尾信号が入力されたら、追尾パラメータ格納部７に保存しておいた追尾パラメータを用いてここまでの処理を繰り返し、時々刻々の追尾結果を出力する。なお、上記の各実施例では、センサは三個としているが、２個以上であれば多数存在しても可能である。
【図面の簡単な説明】
【００５８】
【図１】本発明の一実施例の機能ブロック図である。
【図２】本発明の一実施例の動作イメージ図である。
【図３】本発明の一実施例の動作を示すフローチャートである。
【図４】従来手法と本発明手法の距離誤差比較を示す図である。
【図５】Ｓ／Ｎが悪い状況における従来手法と本発明手法の距離誤差比較を示す図である。
【図６】本発明の他の実施例の機能ブロック図である。
【図７】本発明の他の実施例の動作イメージ図である。
【図８】本発明の他の実施例の動作を示すフローチャートである。
【図９】従来技術を説明するためのイメージ図である。
【符号の説明】
【００５９】
１  信号入力部
２  時間整合処理部
３  追尾フィルタ
４  全観測追尾信号判定部
５  現在時刻判定部
６  予測追尾信号外挿部
７  追尾パラメータ格納部
８  追尾結果出力部
９  制御部（ＣＰＵ）
１０  メモリ
１１  時間遷移処理部
【図１】
2008008691.tif 000005
【図２】
2008008691.tif 000006
【図３】
2008008691.tif 000007
【図４】
2008008691.tif 000008
【図５】
2008008691.tif 000009
【図６】
2008008691.tif 000010
【図７】
2008008691.tif 000011
【図８】
2008008691.tif 000012
【図９】
2008008691.tif 000013