(57)【要約】【課題】信号パワーと雑音パワーとが同等になっても、信号を継続的に検出できる信号分析装置を得る。【解決手段】受波信号と基底関数セットとを用いて、基底関数セットからの基底関数の選択と選択された基底関数を受波信号から差し引く演算処理とを繰り返すことで、スナップショット毎に受波信号を表す基底関数と基底関数の振幅との組み合わせを推定する信号分析装置であって、スナップショット毎に選択された基底関数を含む基底関数リストを記憶する保持手段と、未選択の基底関数について分析対象信号の存在の可能性を示す存在可能性値を算出する可能性値算出手段と、算出された存在可能性値に基づいて未選択の基底関数を基底関数リストに追加するか否かを判定する追加判定手段と、追加された基底関数を含む基底関数リストを用いて追加された基底関数の振幅を計算する再推定手段とを有するものである。【選択図】図１
2020085656.tif 000002
【特許請求の範囲】
【請求項１】
受波信号と基底関数セットとを用いて、前記基底関数セットからの基底関数の選択と選択された基底関数を前記受波信号から差し引く演算処理とを繰り返すことで、観測時間単位であるスナップショット毎に前記受波信号を表す基底関数と該基底関数の振幅との組み合わせを推定する信号分析装置であって、
前記スナップショット毎に選択された基底関数を含む基底関数リストを記憶する保持手段と、
未選択の基底関数について分析対象信号の存在の可能性を示す存在可能性値を算出する可能性値算出手段と、
算出された前記存在可能性値に基づいて前記未選択の基底関数を前記基底関数リストに追加するか否かを判定する追加判定手段と、
追加された基底関数を含む基底関数リストを用いて前記追加された基底関数の振幅を計算する再推定手段と、
を有する信号分析装置。
【請求項２】
前記保持手段は、過去に選択された基底関数の振幅の統計量を記憶し、
前記再推定手段によって算出された振幅と前記統計量とに基づいて、前記追加された基底関数を前記基底関数リストから除外するか否かを判定する除外判定手段を、さらに有する、請求項１に記載の信号分析装置。
【請求項３】
前記基底関数セットは正弦関数および余弦関数のペアで構成され、
前記受波信号は狭帯域信号に属する、請求項１または２に記載の信号分析装置。
【請求項４】
前記受波信号は、複数のセンサを含むセンサアレイで観測される信号である、請求項１または２に記載の信号分析装置。
【請求項５】
前記基底関数はステアリングベクトルであり、
前記センサアレイが直線等間隔アレイであり、
前記スナップショット毎に前記受波信号を表す前記ステアリングベクトルに基づいて前記振幅を推定する推定手段をさらに有し、
前記推定手段は、前記受波信号の到来方位を推定する、請求項４に記載の信号分析装置。
【請求項６】
前記可能性値算出手段は、
過去に選択された基底関数の重み付け平均値から、前記未選択の基底関数について前記存在可能性値を算出する、請求項１〜５のいずれか１項に記載の信号分析装置。
【請求項７】
受波信号と基底関数セットとを用いて、前記基底関数セットからの基底関数の選択と選択された基底関数を前記受波信号から差し引く演算処理とを繰り返すことで、観測時間単位であるスナップショット毎に前記受波信号を表す基底関数と該基底関数の振幅との組み合わせを推定する信号分析方法であって、
前記スナップショット毎に選択された基底関数を含む基底関数リストを記憶し、
未選択の基底関数について分析対象信号の存在の可能性を示す存在可能性値を算出し、
算出された前記存在可能性値に基づいて前記未選択の基底関数を前記基底関数リストに追加するか否かを判定し、
追加された基底関数を含む基底関数リストを用いて前記追加された基底関数の振幅を計算する、
信号分析方法。
【請求項８】
過去に選択された基底関数の振幅の統計量を記憶し、
前記追加された基底関数の振幅と前記統計量とに基づいて、前記追加された基底関数を前記基底関数リストから除外するか否かを判定する、請求項７に記載の信号分析方法。
【請求項９】
前記基底関数セットは正弦関数および余弦関数のペアで構成され、
前記受波信号は狭帯域信号に属する、請求項７または８に記載の信号分析方法。
【請求項１０】
前記受波信号は複数のセンサを含むセンサアレイで観測される信号である、請求項７または８に記載の信号分析方法。
【請求項１１】
前記基底関数はステアリングベクトルであり、
前記センサアレイが直線等間隔アレイであり、
前記スナップショット毎に前記受波信号を表す前記ステアリングベクトルに基づいて前記振幅を推定する際、前記受波信号の到来方位を推定する、請求項１０に記載の信号分析方法。
【請求項１２】
過去に選択された基底関数の重み付け平均値から、前記未選択の基底関数について前記存在可能性値を算出する、請求項７〜１１のいずれか１項に記載の信号分析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、信号を追尾する信号分析装置および信号分析方法に関する。
【背景技術】
【０００２】
従来、狭帯域信号処理の解析手法として、ＦＦＴ（Fast Fourier Transform）が広く用いられている。一方、低ＳＮＲ（Signal-to-Noise Ratio）環境でも信号検出が可能な手法として、ＦＯＳ（Fast Orthogonal Search）が提案されている（例えば、非特許文献１参照）。
【０００３】
ＦＯＳはデータ長に依らず多数の基底関数セットを用意でき、そのセットから基底関数の選択、残差の計算、基底関数の選択を繰り返して周波数解析を行う。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Abdalla Osman, Aboelamgd Nourledin, Naser El-Sheimy, Jim Theriault and Scott Campbell, “Improved target detection and bearing estimation utilizing fast orthogonal search for real-time spectral analysis”, Measurement Science and Technology, 2009
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、非特許文献１に開示された方法では、信号パワーと雑音パワーとが同等になるにつれて、信号も雑音も選択される確率は同等となるため、信号の継続的な検出は難しくなってくる。
【課題を解決するための手段】
【０００６】
本発明に係る信号分析装置は、受波信号と基底関数セットとを用いて、前記基底関数セットからの基底関数の選択と選択された基底関数を前記受波信号から差し引く演算処理とを繰り返すことで、観測時間単位であるスナップショット毎に前記受波信号を表す基底関数と該基底関数の振幅との組み合わせを推定する信号分析装置であって、前記スナップショット毎に選択された基底関数を含む基底関数リストを記憶する保持手段と、未選択の基底関数について分析対象信号の存在の可能性を示す存在可能性値を算出する可能性値算出手段と、算出された前記存在可能性値に基づいて前記未選択の基底関数を前記基底関数リストに追加するか否かを判定する追加判定手段と、追加された基底関数を含む基底関数リストを用いて前記追加された基底関数の振幅を計算する再推定手段と、を有するものである。
【０００７】
本発明に係る信号分析方法は、受波信号と基底関数セットとを用いて、前記基底関数セットからの基底関数の選択と選択された基底関数を前記受波信号から差し引く演算処理とを繰り返すことで、観測時間単位であるスナップショット毎に前記受波信号を表す基底関数と該基底関数の振幅との組み合わせを推定する信号分析方法であって、前記スナップショット毎に選択された基底関数を含む基底関数リストを記憶し、未選択の基底関数について分析対象信号の存在の可能性を示す存在可能性値を算出し、算出された前記存在可能性値に基づいて前記未選択の基底関数を前記基底関数リストに追加するか否かを判定し、追加された基底関数を含む基底関数リストを用いて前記追加された基底関数の振幅を計算するものである。
【発明の効果】
【０００８】
本発明によれば、信号の存在する可能性が高い基底関数を基底関数リストに追加して振幅を再推定しているので、追尾対象の信号を継続的に検出することができる。
【図面の簡単な説明】
【０００９】
【図１】本発明の実施の形態１に係る信号分析装置の一構成例を示す図である。
【図２】図１に示した信号分析装置のハードウェア構成の一例を示すブロック図である。
【図３】ＦＯＳ処理におけるスナップショットを説明するための図である。
【図４】比較例の信号分析装置の一構成例を示すブロック図である。
【図５】図４に示した信号分析装置が実行する信号分析方法の手順を示すフロー図である。
【図６】ＦＯＳ処理におけるスナップショット毎に変動する雑音パワーの一例を示す図である。
【図７】比較例において、時間経過に伴って振幅が減衰する信号に対するＦＯＳの結果を示す図である。
【図８】図１に示した信号分析装置が実行する信号分析方法の手順を示すフロー図である。
【図９】本実施の形態１において、時間経過に伴って振幅が減衰する信号に対するＦＯＳの結果を示す図である。
【図１０】本発明の実施の形態２に係る信号分析装置の一構成例を示す図である。
【図１１】図１０に示した信号分析装置が実行する信号分析方法の手順を示すフロー図である。
【図１２】本発明の実施の形態３に係る信号分析装置の一構成例を示す図である。
【図１３】図１２に示した信号分析装置が実行する信号分析方法の手順を示すフロー図である。
【図１４】本実施の形態３において、Ｓ個のターゲットからの音を等間隔で直線に配置したＲ個のセンサで受波するモデルを示す図である。
【発明を実施するための形態】
【００１０】
実施の形態１．
本実施の形態１の信号分析装置について説明する。図１は、本発明の実施の形態１に係る信号分析装置の一構成例を示す図である。図１に示すように、信号分析装置１は、記憶部１１および制御部１２を有する。記憶部１１は保持手段１３を有する。制御部１２は、ＦＯＳ処理部２０と、可能性値算出手段１４と、追加判定手段１５と、再推定手段１６とを有する。ＦＯＳ処理部２０は、基底関数選択手段２１、推定手段２２、残差計算手段２３および終了判定手段２４を有する。
【００１１】
図２は、図１に示した信号分析装置のハードウェア構成の一例を示すブロック図である。信号分析装置１は、例えば、コンピュータを含む情報処理装置である。記憶部１１は、制御部１２が実行する演算処理の結果を記憶する。記憶部１１は、例えば、ＨＤＤ（Hard Disk Drive）装置である。制御部１２は、プログラムを記憶するメモリ５２と、プログラムにしたがって処理を実行するＣＰＵ（Central Processing Unit）５１とを有する。ＣＰＵ５１がプログラムを実行することで、図１に示した、基底関数選択手段２１、推定手段２２、残差計算手段２３、終了判定手段２４、可能性値算出手段１４、追加判定手段１５および再推定手段１６が構成される。
【００１２】
なお、基底関数選択手段２１、推定手段２２、残差計算手段２３、終了判定手段２４、可能性値算出手段１４、追加判定手段１５および再推定手段１６が備える機能のうち、一部または全部が専用回路で構成されてもよい。専用回路は、例えば、ＡＳＩＣ（Application Specific Integrated Circuit）である。信号処理の一部または全部を専用回路で構成することで、信号処理の高速化を図ることができる。
【００１３】
本実施の形態１の信号分析装置１の構成を詳細に説明する前に、比較例の信号分析装置を説明する。比較例の信号分析装置は、スペクトログラム等の時間、周波数および信号成分の強さで表される三次元グラフ上にＦＯＳの計算結果を表示する。
【００１４】
はじめに、比較例の信号分析装置が行う信号分析処理に用いられるＦＯＳについて説明する。ＦＯＳはセンサの受波信号
【００１５】
【数１】
2020085656.tif 000003
を、事前に用意したＭ個の基底関数からなる基底関数セット
【００１６】
【数２】
2020085656.tif 000004
の中から、受波信号ｘとの誤差を小さくするＬ個の基底関数
【００１７】
【数３】
2020085656.tif 000005
を当てはめて、狭帯域信号の捜索を行う。なお、文章中においては、基底関数Ｐのベクトルを「Ｐ＾」で表す。例えば、任意のｌの基底関数ＰをＰ＾ｌと表記し、Ｌ個の基底関数ＰのリストをＰ＾Ｌと表記する。
【００１８】
時系列で変化する振幅を示す受波信号に対する観測時間単位をスナップショットと称する。ここで、Ｎは、１スナップショットの処理で用いる受波信号のサンプル数である。図３は、ＦＯＳ処理におけるスナップショットを説明するための図である。図３に示すように、１スナップショットは、０〜１００％未満のオーバーラップ処理で設定される。受波信号に対するＬ個の基底関数の当てはめを、式（１）で表す。
【００１９】
【数４】
2020085656.tif 000006
式（１）において、
【００２０】
【数５】
2020085656.tif 000007
は、ｘが式（１）の第１項で構成されると仮定したときの真値との誤差を表す。式（１）におけるａｌは基底関数Ｐ＾ｌに対する重みである。
【００２１】
基底関数の一例として、解析対象の信号が狭帯域信号であるとすると、正弦関数と余弦関数とのペアとなる。このとき、Ｍ個の、正弦関数および余弦関数のペアで構成される基底関数セットＰは、次の式（２）で表される。
【００２２】
【数６】
2020085656.tif 000008
式（２）において、ωｍ＝２πｆｍは角周波数であり、Ｍ、ｆｍは任意に設定することができる値である。式（２）の基底関数を用いる場合、式（１）は次の式（３）で表せる。
【００２３】
【数７】
2020085656.tif 000009
ＦＯＳは、基底関数セットＰおよび受波信号ｘから、式（３）に示すεを最小とする基底関数
【００２４】
【数８】
2020085656.tif 000010
の組み合わせとその基底関数の個数Ｌとを推定する手法である。
【００２５】
図４は、比較例の信号分析装置の一構成例を示すブロック図である。比較例の信号分析装置１００は、保持手段１３を含む記憶部１０１と、ＦＯＳ処理部２０を含む制御部１０２とを有する。制御部１０２は、図に示さない、メモリおよびＣＰＵを有する。メモリが記憶するプログラムをＣＰＵが実行することで、基底関数選択手段２１、推定手段２２、残差計算手段２３および終了判定手段２４が構成される。記憶部１０１は、例えば、ＨＤＤ装置である。
【００２６】
保持手段１３は、予め用意された複数の基底関数を含む基底関数セットを記憶する。また、保持手段１３は、受波信号に対する捜索回数が１から（ｔ−１）までの信号分析処理でスナップショット毎に選択された基底関数を格納した基底関数リストを記憶する。保持手段１３は、残差計算手段２３が算出する残差信号ｘｅ（ｔ）を記憶する。
【００２７】
基底関数選択手段２１は、捜索回数１回目では受波信号ｘを入力とし、捜索回数２回目以降では残差信号ｘｅ（ｔ−１）を入力とし、捜索回数ｔにおける基底関数の選択処理を行う。基底関数選択手段２１は、選択した基底関数を含む基底関数リストＰ＾ｔと、基底関数に対する重みｇｔおよび係数行列αｔとを出力する。
【００２８】
推定手段２２は、基底関数リストＰ＾ｔと重みｇｔおよび係数行列αｔとを入力とし、選択された基底関数の振幅Ａｔを計算する。推定手段２２は、基底関数リストＰ＾ｔと振幅Ａｔとを出力する。残差計算手段２３は、基底関数リストＰ＾ｔと振幅Ａｔとを入力とし、受波信号ｘとの残差計算を行い、残差信号ｘｅ（ｔ）を出力する。
【００２９】
終了判定手段２４は、残差信号ｘｅ（ｔ）を入力とし、終了条件が満たされるか否かを判定する。終了条件は、例えば、１つ前のスナップショットで推定された振幅が予め決められた閾値以下になることである。また、終了条件の別の例として、残差信号の残差パワーが予め決められた閾値以下になるという条件が挙げられる。終了判定の結果、終了条件が満たされない場合、終了判定手段２４は、基底関数選択手段２１に残差信号ｘｅ（ｔ−１）を出力する。終了判定の結果、終了条件が満たされる場合、終了判定手段２４は、捜索終了時の回数ＬにおけるＰ＾Ｌと振幅ＡＬとを出力する。
【００３０】
図５は、図４に示した信号分析装置が実行する信号分析方法の手順を示すフロー図である。図５に示す手順の概要を説明する。
【００３１】
信号分析装置１００は、基底関数の選択（ステップＳ５０１）、振幅の計算（ステップＳ５０２）および残差信号の計算（ステップＳ５０３）で狭帯域信号から追尾対象の信号を捜索する。信号分析装置１００は、ステップＳ５０４において、検出すべき基底関数の成分がないと判定するまで、捜索回数を順次更新し、ステップＳ５０１〜Ｓ５０３の処理を繰り返す。捜索が終了すると、信号分析装置１００は、分析結果としてスペクトルを出力する。信号分析装置１００は、スペクトル出力の後、図３に示したように、スナップショットｋを（ｋ＋１）に更新し、スナップショット（ｋ＋１）の受波信号に対して、ステップＳ５０１〜Ｓ５０４の処理を繰り返す。ｋは任意の正の整数である。その後、信号分析装置１００は、スナップショット毎にスペクトルをスペクトログラム等の三次元グラフに表示する（ステップＳ５０５）。
【００３２】
以下に、図５に示した手順を詳しく説明する。なお、本実施の形態１において、スペクトルは、各周波数ｆｌに対する振幅
【００３３】
【数９】
2020085656.tif 000011
のことである。
【００３４】
［基底関数の選択（ステップＳ５０１）］
基底関数選択手段２１は、捜索回数（ｔ−１）においてステップＳ５０３で算出された残差信号ｘｅ（ｔ−１）（ｔ＝０でｘｅ（ｔ−１）＝ｘ）から最もよく当てはまる基底関数を基底関数セットＰの中から選択する。
ここで、捜索回数（ｔ−１）までに選択した基底関数を含む基底関数リスト
【００３５】
【数１０】
2020085656.tif 000012
は、次のように表される。
【００３６】
【数１１】
2020085656.tif 000013
また、捜索回数（ｔ−１）まで、基底関数Ｐうち、未選択の基底関数を含む未選択基底関数リストＵは、以下の通りである。
【００３７】
【数１２】
2020085656.tif 000014
【００３８】
未選択基底関数リストＵにおいて、Ｑは用意したＭ個の基底関数の内の未選択の個数で、Ｑ＝Ｍ−（ｔ−１）である。未選択基底関数リストＵの中からｑ個目の基底関数を捜索する場合、基底関数選択手段２１は、次式の行列Ｚｑを作成する。
【００３９】
【数１３】
2020085656.tif 000015
ただし、
【００４０】
【数１４】
2020085656.tif 000016
【００４１】
である。ここで、基底関数セットＰ内の基底関数は基底関数同士で相関をもつ場合がある。相関を持つとスペクトル漏れが生じ、漏れ込んだ成分に微弱な信号が埋もれ、見逃す場合がある。そこで、Ｚｑ内の基底関数同士を直交化させた行列をＷｑとすると、基底関数選択手段２１は、以下のように基底関数同士の直交化を図る。
【００４２】
【数１５】
2020085656.tif 000017
この行列Ｗｑにおいて、
【００４３】
【数１６】
2020085656.tif 000018
である。式（３）の直交化後の式を、式（４）に示す。
【００４４】
【数１７】
2020085656.tif 000019
【００４５】
式（４）において、ｇｑ，ｌは直交化したｗｑ，ｌに与える重みを表す。ｅは直交化後の誤差εを表す。係数αｑ，ｌ，ｒを定義すると、高速化のために、基底関数選択手段２１は、ｗｑ，ｌを直接求めることなく、ｇｑ，ｌとαｑ，ｌ，ｒを、以下のように再帰的に求める。
【００４６】
【数１８】
2020085656.tif 000020
【００４７】
【数１９】
2020085656.tif 000021
【００４８】
【数２０】
2020085656.tif 000022
【００４９】
【数２１】
2020085656.tif 000023
【００５０】
ここで、式（５）〜式（８）において｛ｌ＝０，１，・・・，２ｔ＋１｝である。
基底関数選択手段２１は、未選択基底関数リストＵの中から残差信号との誤差パワーを最小とする基底関数
【００５１】
【数２２】
2020085656.tif 000024
を選択する。誤差パワーｅ２は次式で算出する。
【００５２】
【数２３】
2020085656.tif 000025
【００５３】
式（９）におけるｘｅ（ｔ−１）は残差信号であり、ステップＳ５０３で詳細に説明する。基底関数選択手段２１は、式（９）によって選択した基底関数
【００５４】
【数２４】
2020085656.tif 000026
を用いて、次式のように基底関数リストＰ＾ｔを更新する。
【００５５】
【数２５】
2020085656.tif 000027
このとき、基底関数リストＰ＾ｔ内の各基底関数に対応した重み
【００５６】
【数２６】
2020085656.tif 000028
および係数
【００５７】
【数２７】
2020085656.tif 000029
はステップＳ５０２の処理で用いられる。そのため、基底関数選択手段２１は、
【００５８】
【数２８】
2020085656.tif 000030
を重みベクトルとして
【００５９】
【数２９】
2020085656.tif 000031
を保持手段１３に保存し、
【００６０】
【数３０】
2020085656.tif 000032
を次式の下三角行列αｔとして保持手段１３に保存する。
【００６１】
【数３１】
2020085656.tif 000033
【００６２】
［振幅の計算（ステップＳ５０２）］
ＦＯＳは、ステップＳ５０１で選択された基底関数リスト
【００６３】
【数３２】
2020085656.tif 000034
に対応する周波数
【００６４】
【数３３】
2020085656.tif 000035
のみに推定された振幅
【００６５】
【数３４】
2020085656.tif 000036
をもち、ステップＳ５０１で未選択の基底関数に対応する周波数の振幅値はゼロである。
推定手段２２は、選択された各周波数に対する振幅
【００６６】
【数３５】
2020085656.tif 000037
を、直交化後の重みｇｔと係数の行列αｔとを用いて、次の式（１１）のように再帰的に求める。
【００６７】
【数３６】
2020085656.tif 000038
式（１１）において、
【００６８】
【数３７】
2020085656.tif 000039
である。式（１１）で求められたａｌは振幅ベクトルＡｔとして保持手段１３に、以下のように保存される。
【００６９】
【数３８】
2020085656.tif 000040
【００７０】
［残差信号の計算（ステップＳ５０３）］
残差計算手段２３は、ステップＳ５０１で選択された基底関数リストＰ＾ｔとステップＳ５０２で推定された振幅Ａｔとを用いて、受波信号ｘから、選択された基底関数リストの成分を除くことで、残差信号ｘｅ（ｔ）を計算する。ｘｅ（ｔ）の計算は、次式で行われる。
【００７１】
【数３９】
2020085656.tif 000041
【００７２】
［終了判定（ステップＳ５０４）］
ＦＯＳは、受波信号と残差信号との誤差パワーを最小とする基底関数を選択する演算処理なので、捜索終盤は雑音のみしかない可能性が高い。捜索回数を必要以上に多くしてしまうと、無駄な計算と誤警報とが増加する要因となる。そのため、終了条件を設け、捜索の途中で打ち切ることが妥当である。終了判定手段２４は、終了条件が満たされるか否かを判定する（ステップＳ５０４）。
【００７３】
ステップＳ５０４で、終了条件が満たされる場合、終了判定手段２４は、ステップＳ５０１で選択された基底関数リストＰ＾ｔに対応する周波数とステップＳ５０２で推定された振幅Ａｔとをスペクトログラム等の三次元グラフ上に表示する（ステップＳ５０５）。このときの捜索回数ｔが選択された基底関数の総数Ｌ＝ｔとなる。一方、ステップＳ５０４において、終了条件が満たされない場合、終了判定手段２４は次の捜索を基底関数選択手段２１に指示する（ステップＳ５０１）。ステップＳ５０３で算出された残差信号ｘｅ（ｔ）は次のステップＳ５０１〜Ｓ５０３の処理で用いられる。
【００７４】
ＦＯＳは、ＦＦＴと比較して、同じサンプル数で周波数分解能を詳細に設定できる。さらに、ＦＯＳは、窓関数を必要としないことからＳＮＲの劣化を防げる。そのため、ＦＯＳは、ＦＦＴに比べて低ＳＮＲ環境での信号検出能力の向上が期待できる。その結果、狭帯域信号検出の改善が期待できる。
【００７５】
しかし、低ＳＮＲ環境におけるＦＯＳは、信号成分が残っていたとしても雑音が選択される確率が高まるため、継続的な信号検出は難しい。
【００７６】
この問題の原因は、図５に示したステップＳ５０４における終了判定にあると考えられる。終了条件は、１つ前のスナップショットで推定された振幅値に対する閾値判定および残差パワーに対する閾値判定等が挙げられる。これらの閾値判定によって、残差信号に信号が含まれていない状況での早期の捜索終了が期待できる。しかし、定常性の雑音であっても、スナップショット毎では局所的な雑音パワーは変動する。図６は、ＦＯＳ処理におけるスナップショット毎に変動する雑音パワーの一例を示す図である。図６に示すように、雑音が信号パワーを上回るようなスナップショットが存在し、そのようなスナップショットでは信号を捜索する前に上述した閾値に達してしまうことが考えられる。
【００７７】
図７は、比較例において、時間経過に伴って振幅が減衰する信号に対するＦＯＳの結果を示す図である。比較例の信号分析方法では、信号の振幅が小さくなるに伴って、信号が存在する周波数を選択する前に捜索を終了してしまうことが多発し、図７に示すように時間方向で継続的に信号を検出することが難しくなる。
【００７８】
ＦＯＳの終了条件を緩和することで、信号が存在する周波数を選択することにつながるが、比較例の信号分析方法の説明で述べたように捜索回数を単純に増やすことは、誤検出の増加につながる。本実施の形態１は、捜索の終了条件を緩和せずに継続的に分析対象信号が存在する周波数を選択することで、従来のＦＯＳにおける課題を解決するものである。
【００７９】
本実施の形態１の信号分析装置１の構成を、図１を参照して詳しく説明する。なお、比較例の信号分析装置１００と同様な構成については同一の符号を付し、その詳細な説明を省略する。本実施の形態１の信号分析装置１は、図１に示すように、図４および図５を参照して説明した信号分析装置１００に、可能性値算出手段１４、追加判定手段１５および再推定手段１６が追加された構成である。
【００８０】
可能性値算出手段１４は、捜索終了した基底関数リストＰ＾Ｌと全ての基底関数を含む基底関数セットＰとを入力とし、基底関数リストＰ＾Ｌに対する周波数ビンの検出回数のカウントを行い、基底関数リストＰ＾Ｌと重み付き検出回数ＤＴｍｅａｎとを出力する。
【００８１】
追加判定手段１５は、基底関数リストＰ＾Ｌと重み付き検出回数ＤＴｍｅａｎとを入力とし、重み付き検出回数ＤＴｍｅａｎに基づいて未選択の基底関数を新規に追加するか否かを判定する。追加判定手段１５は、未選択の基底関数を追加する場合、基底関数リストＰ＾Ｌを更新し、更新された基底関数リストＰ＾Ｌを出力する。
【００８２】
再推定手段１６は、更新された基底関数リストＰ＾Ｌを入力とし、基底関数リストＰ＾Ｌに基づいて振幅ＡＬを推定する。そして、再推定手段１６は、スペクトログラム等の三次元グラフ上の推定結果を表示する。
【００８３】
次に、本実施の形態１の信号分析装置１の動作を説明する。図８は、図１に示した信号分析装置が実行する信号分析方法の手順を示すフロー図である。図８に示すステップＳ１０１〜Ｓ１０４およびＳ１０８の処理は図５を参照して説明したステップＳ５０１〜Ｓ５０５と同様なため、その詳細な説明を省略する。
【００８４】
［検出回数のカウント（ステップＳ１０５）］
可能性値算出手段１４は、事前に用意されたＭ個の基底関数のセット
【００８５】
【数４０】
2020085656.tif 000042
に対して、ステップＳ１０４の終了判定で得られた基底関数リスト
【００８６】
【数４１】
2020085656.tif 000043
を用いて、どの基底関数が検出されたか、検出回数をカウントする。具体的には、
【００８７】
【数４２】
2020085656.tif 000044
の関係であるため、
【００８８】
【数４３】
2020085656.tif 000045
に一致する
【００８９】
【数４４】
2020085656.tif 000046
の周波数ｆｍに対し、可能性値算出手段１４は、次式のように検出回数バッファＤＴへ検出回数をカウントする。
【００９０】
【数４５】
2020085656.tif 000047
ここで、ｋはスナップショット番号である。可能性値算出手段１４は、α個前までのスナップショットのＤＴを用いて、ｋ番目のスナップショットで検出されていない基底関数の（ｋ−１）から（ｋ−α）番目までの検出回数を次式のように見積もる。
【００９１】
【数４６】
2020085656.tif 000048
式（１４）において、｜ＤＴ（ｆｍ，ｋ）−１｜は現スナップショットのＦＯＳで選択済みの基底関数を除くための演算である。β（ｋ）は重み係数であり、過去α個の平均的な検出回数で、可能性値算出手段１４は、分析対象となる狭帯域信号の存在を判定する。この方法に限らず、直近の検出回数の重みの比率を上げるなどして、直近の検出回数を重視するなど自由に設定できる。重み付き検出回数ＤＴｍｅａｎは、分析対象信号の存在の可能性を示す存在可能性値に相当する。
【００９２】
［基底関数リストへ新規追加（ステップＳ１０６）］
ステップＳ１０５において可能性値算出手段１４が式（１４）で作成した重み付き検出回数
【００９３】
【数４７】
2020085656.tif 000049
に基づいて、追加判定手段１５は、新しく基底関数を追加して基底関数リストを更新する。具体的には、追加判定手段１５は、新規追加する基底関数として、ＤＴｍｅａｎを用いて重み付き検出回数が多い周波数ｆｍに対応する基底関数
【００９４】
【数４８】
2020085656.tif 000050
を選択し、次式のように基底関数リストＰ＾Ｌを更新する。
【００９５】
【数４９】
2020085656.tif 000051
ここで、θＤＴは閾値である。追加判定手段１５は、閾値θＤＴよりも重み付き検出回数の多い周波数ｆｍに対応する基底関数は分析対象信号が存在する可能性が高いと判断し、その基底関数をＰ＾Ｌに追加する。
【００９６】
［振幅の再計算（ステップＳ１０７）］
再推定手段１６は、ステップＳ１０６において基底関数リストに新規で追加された基底関数の振幅を、ステップＳ１０６で更新された基底関数リストＰ＾Ｌと式（５）〜式（８）および式（１１）〜式（１２）とを用いて計算する。再推定手段１６は、更新された基底関数リストＰ＾Ｌの基底関数に対応する周波数と基底関数リストＰ＾Ｌの基底関数に対応する振幅Ａｔとをスペクトログラム等の三次元グラフ上に表示する。
【００９７】
本実施の形態１の信号分析装置１は、スナップショット毎に受波信号を表す基底関数と基底関数の振幅との組み合わせを推定するＦＯＳ処理部２０と、保持手段１３と、可能性値算出手段１４と、追加判定手段１５と、再推定手段１６とを有する。保持手段１３は、スナップショット毎に選択された基底関数を含む基底関数リストを記憶する。可能性値算出手段１４は、未選択の基底関数について分析対象信号の存在の可能性を示す存在可能性値を算出する。追加判定手段１５は、算出された存在可能性値に基づいて未選択の基底関数を基底関数リストに追加するか否かを判定する。再推定手段１６は、追加された基底関数を含む基底関数リストを用いて、追加された基底関数の振幅を計算する。
【００９８】
本実施の形態１では、信号の存在する周波数を継続的に検出するために、信号の有無の判定にＦＯＳの周波数ビン毎の検出回数を利用している。無相関雑音を仮定した場合、雑音を信号として誤検出した場合の検出周波数はスナップショット毎で異なる。そのため、雑音に対する検出周波数は時間的に不安定になる。一方、信号が存在する周波数は継続的に検出することが困難であっても、検出周波数は時間的に安定する。そのため、周波数ビン毎に検出回数をカウントすることで、信号の存在する周波数ビンの検出回数が増加し、それ以外の雑音のみの周波数ビンの検出回数は、信号の存在する周波数ビンよりは増加しないと考えられる。本実施の形態１の信号分析装置１は、過去の複数のスナップショットの捜索に基づいた検出回数を利用することで、信号の存在する可能性が高い周波数ビンを優先的に捜索することができるため、信号検出の時間的な継続性が向上する。
【００９９】
図９は、本実施の形態１において、時間経過に伴って振幅が減衰する信号に対するＦＯＳの結果を示す図である。図９を図７と比較すると、図７に比べて、信号をより継続的に検出できていることがわかる。本実施の形態１では、捜索の終了条件を緩和せずに、信号の存在する可能性が高い基底関数を基底関数リストに追加して振幅を再推定しているので、追尾対象の信号を継続的に検出することができる。
【０１００】
実施の形態２．
本実施の形態２は、実施の形態１で説明した信号分析方法において、追加した基底関数が信号分析対象でない場合、追加した基底関数を基底関数リストから除外するものである。
【０１０１】
本実施の形態２の信号分析装置の構成を説明する。本実施の形態２では、実施の形態１で説明した構成と同一の構成に同一の符号を付し、その詳細な説明を省略する。図１０は、本発明の実施の形態２に係る信号分析装置の一構成例を示す図である。
【０１０２】
図１０に示すように、信号分析装置１ａは、記憶部１１および制御部１２を有する。制御部１２は、ＦＯＳ処理部２０と、可能性値算出手段１４と、追加判定手段１５と、再推定手段１６と、除外判定手段３１とを有する。保持手段１３は、過去に選択された基底関数の振幅の統計量が格納される振幅値バッファを記憶する。
【０１０３】
除外判定手段３１は、基底関数リストＰ＾Ｌと振幅ＡＬおよび振幅値バッファとを入力とし、追加された基底関数の除外判定を行い、判定結果を反映させた基底関数リストＰ＾Ｌと振幅ＡＬとを出力する。
【０１０４】
次に、本実施の形態２の信号分析装置１ａの動作を説明する。図１１は、図１０に示した信号分析装置が実行する信号分析方法の手順を示すフロー図である。図１１に示すステップＳ２０１〜Ｓ２０４およびＳ２１０の処理は図５を参照して説明したステップＳ５０１〜Ｓ５０５と同様なため、その詳細な説明を省略する。また、図１１に示すステップＳ２０５〜Ｓ２０７の処理は図８を参照して説明したステップＳ１０５〜１０７と同様なため、その詳細な説明を省略する。
【０１０５】
ステップＳ２０４の後、終了判定手段２４は、基底関数リストＰ＾Ｌと振幅ＡＬとを出力した後、振幅値バッファ
【０１０６】
【数５０】
2020085656.tif 000052
への振幅値の保存を、次式のように行う。
【０１０７】
【数５１】
2020085656.tif 000053
【０１０８】
［追加した基底関数の除外判定（ステップＳ２０８）］
除外判定手段３１は、ステップＳ２０６で更新された基底関数リストＰ＾Ｌと、ステップＳ２０７で算出された基底関数リストＰ＾Ｌに対する振幅ＡＬと、振幅値バッファＡｍｐとを用いて、追加された基底関数の除去判定を行う。具体的には、除外判定手段３１は、追加した基底関数
【０１０９】
【数５２】
2020085656.tif 000054
に対する過去の振幅の統計量として平均および分散を、振幅値バッファＡｍｐから次式のように計算する。
【０１１０】
【数５３】
2020085656.tif 000055
このとき、
【０１１１】
【数５４】
2020085656.tif 000056
の平均と分散は非零成分のみで算出され、αは非零のみの過去のスナップショット数とする。除外判定手段３１は、算出した平均と分散に基づいて次式のような
【０１１２】
【数５５】
2020085656.tif 000057
に対する除外判定を行う。
【０１１３】
【数５６】
2020085656.tif 000058
ここで、Ｖは重みである。２つの関係式のうち、上段の関係式が満たされる場合、除外判定手段３１は、追加された基底関数を基底関数リストＰ＾Ｌから除外する（ステップＳ２０９）。一方、下段の関係式が満たされる場合、除外判定手段３１は、追加された基底関数を基底関数リストＰ＾Ｌに維持する。除外判定手段３１は、基底関数リストＰ＾Ｌの基底関数に対応する周波数と基底関数リストＰ＾Ｌの基底関数に対する振幅ＡＬとをスペクトログラム等の三次元グラフ上に表示する（ステップＳ２１０）。
【０１１４】
本実施の形態２の信号分析装置１ａは、再推定手段１６によって算出された振幅と過去の基底関数の振幅の統計量とに基づいて、追加された基底関数を基底関数リストから除外するか否かを判定する除外判定手段３１を有するものである。
【０１１５】
実施の形態１では、追加判定で条件を満たす基底関数が追加されるが、追加された基底関数の信号が雑音であるかについての判別は行われていない。そのため、現在のスナップショットにおいて、追加された基底関数の信号が雑音である可能性もある。そこで、狭帯域信号は定常的な信号であるため、本実施の形態２の信号分析装置１ａでは、過去に検出された基底関数の振幅値から、追加された基底関数の統計量を算出し、その統計量に基づいて、追加された基底関数の信号が雑音であるかを判定している。信号分析対象となった信号が雑音である場合、算出した信号の統計量にしたがわないため、本実施の形態２では、その信号が雑音か否かを判別できる。そのため、実施の形態１と同様な効果が得られるだけでなく、追加された基底関数の信号が雑音である場合の誤警報を抑制することができる。その結果、誤警報の増加を回避できる。
【０１１６】
実施の形態３．
本実施の形態３は、整相処理による信号到来方位推定に対して、実施の形態１を適用したものである。
【０１１７】
本実施の形態３の信号分析装置の構成を説明する。本実施の形態３では、実施の形態１および２で説明した構成と同一の構成に同一の符号を付し、その詳細な説明を省略する。図１２は、本発明の実施の形態３に係る信号分析装置の一構成例を示す図である。
図１２に示すように、信号分析装置１ｂは、記憶部１１および制御部１２を有する。記憶部１１は保持手段１３を有する。制御部１２は、ＦＯＳ処理部２０と、可能性値算出手段１４と、追加判定手段１５と、再推定手段１６とを有する。信号分析装置１ｂには、高速フーリエ変換部（ＦＦＴ部）４０を介して、複数のセンサを有するセンサアレイ５０が接続されている。センサアレイ５０の各センサで観測される信号は、ＦＦＴ部４０によってフーリエ変換された後、制御部１２に入力される。センサアレイ５０の形状は、例えば、直線等間隔アレイである。
【０１１８】
センサアレイ５０が直線等間隔アレイである場合、基底関数は、一例として、ステアリングベクトル（ＳＶ）が考えられる。ステアリングベクトルについては、信号分析装置１ｂの動作の説明で述べる。保持手段１３は、複数のステアリングベクトルを含むＳＶセットを記憶する。
【０１１９】
基底関数選択手段２１は、捜索回数１回目では各センサの受波信号Ｘ（ｆ）を入力とし、捜索回数２回目以降では残差信号Ｘｅ（ｆ，ｔ−１）を入力として、捜索回数ｔにＳＶセットからステアリングベクトルを選択する。基底関数選択手段２１は、選択したステアリングベクトルを含むＳＶリストＰ＾ｔと各ステアリングベクトルに対する重みｇｔおよび係数行列αｔとを出力する。
【０１２０】
推定手段２２は、ＳＶリストＰ＾ｔと重みｇｔおよび係数行列αｔとを入力として、ＳＶリストＰ＾ｔと振幅Ａｔとを出力する。本実施の形態３では、推定手段２２は、入力されるＳＶリストＰ＾ｔを用いて、信号の到来方位を含む整相出力Ａｔを算出してもよい。整相出力の推定については後述する。残差計算手段２３は、ＳＶリストＰ＾ｔと振幅Ａｔとを入力として、残差信号Ｘｅ（ｆ，ｔ）を算出する。残差計算手段２３は、ＳＶリストＰ＾ｔと整相出力Ａｔとを入力として、残差信号Ｘｅ（ｆ，ｔ）を算出してもよい。
【０１２１】
終了判定手段２４は、残差信号Ｘｅ（ｆ，ｔ）を入力とし、終了条件が満たされるか否かを判定する。終了判定手段２４は、終了条件が満たされない場合、残差信号Ｘｅ（ｆ，ｔ−１）を基底関数選択手段２１に出力し、終了条件が満たされる場合、捜索終了時におけるＳＶリストＰ＾Ｌと振幅ＡＬを出力する。終了判定手段２４は、振幅ＡＬの代わりに、信号を推定方位へ整相した結果を出力してもよい。
【０１２２】
可能性値算出手段１４は、捜索終了したＳＶリストＰ＾Ｌと全てのステアリングベクトルを含むＳＶセットＰとを入力とし、ＳＶリストＰ＾Ｌに対する周波数ビンの検出回数のカウントを行い、ＳＶリストＰ＾Ｌと重み付き検出回数ＤＴｍｅａｎとを出力する。
【０１２３】
追加判定手段１５は、ＳＶリストＰ＾Ｌと重み付き検出回数ＤＴｍｅａｎとを入力とし、重み付き検出回数ＤＴｍｅａｎに基づいて未選択のステアリングベクトルを新規に追加するか否かを判定する。追加判定手段１５は、未選択のステアリングベクトルを追加する場合、ＳＶリストＰ＾Ｌを更新し、更新されたＳＶリストＰ＾Ｌを出力する。
【０１２４】
再推定手段１６は、更新されたＳＶリストＰ＾Ｌを入力とし、時間、方位および整相出力の三次元グラフ上に、ＳＶリストＰ＾Ｌに基づいて振幅ＡＬを推定し、推定結果を表示する。再推定手段１６は、更新されたＳＶリストＰ＾ｔを用いて、整相出力ＡＬを算出してもよい。
【０１２５】
次に、本実施の形態３の信号分析装置１ｂの動作を説明する。図１３は、図１２に示した信号分析装置が実行する信号分析方法の手順を示すフロー図である。ここでは、推定手段２２および再推定手段１６は、ＡｔおよびＡＬについて、振幅の代わりに、信号を推定方位へ整相した結果を出力する場合で説明する。
【０１２６】
ステップＳ３０１において、基底関数選択手段２１は、捜索回数１回目では各センサの受波信号Ｘ（ｆ）を入力とし、捜索回数ｔにＳＶセットからステアリングベクトルを選択する。基底関数選択手段２１は、選択したステアリングベクトルを含むＳＶリストＰ＾ｔと各ステアリングベクトルに対する重みｇｔおよび係数行列αｔとを出力する。なお、捜索回数２回目以降では、ステップＳ３０１において、基底関数選択手段２１は、残差信号Ｘｅ（ｆ，ｔ−１）を入力とする。
【０１２７】
ステップＳ３０２において、推定手段２２は、ＳＶリストＰ＾ｔと重みｇｔおよび係数行列αｔとを入力として、ＳＶリストＰ＾ｔと整相出力Ａｔとを出力する。ステップＳ３０３において、残差計算手段２３は、ＳＶリストＰ＾ｔと整相出力Ａｔとを入力として、残差信号Ｘｅ（ｆ，ｔ）を算出する。ステップＳ３０４において、終了判定手段２４は、残差信号Ｘｅ（ｆ，ｔ）が入力されると、終了条件が満たされるか否かを判定する。
【０１２８】
ステップＳ３０４の判定の結果、終了条件が満たされない場合、終了判定手段２４は、残差信号Ｘｅ（ｆ，ｔ−１）を基底関数選択手段２１に出力する（ステップＳ３０１）。一方、ステップＳ３０４の判定の結果、終了条件が満たされる場合、終了判定手段２４は、捜索終了時におけるＳＶリストＰ＾Ｌと整相出力ＡＬを出力する。
【０１２９】
ステップＳ３０５において、可能性値算出手段１４は、捜索終了したＳＶリストＰ＾Ｌと全てのステアリングベクトルを含むＳＶセットＰとを入力とし、ＳＶリストＰ＾Ｌに対する周波数ビンの検出回数のカウントを行う。そして、可能性値算出手段１４は、ＳＶリストＰ＾Ｌと重み付き検出回数ＤＴｍｅａｎとを出力する。
【０１３０】
ステップＳ３０６において、追加判定手段１５は、ＳＶリストＰ＾Ｌと重み付き検出回数ＤＴｍｅａｎとを入力とし、重み付き検出回数ＤＴｍｅａｎに基づいて未選択のステアリングベクトルを新規に追加するか否かを判定する。追加判定手段１５は、未選択のステアリングベクトルを追加する場合、ＳＶリストＰ＾Ｌを更新し、更新されたＳＶリストＰ＾Ｌを出力する。
【０１３１】
ステップＳ３０７において、再推定手段１６は、更新されたＳＶリストＰ＾Ｌを入力とし、ＳＶリストＰ＾Ｌに基づいて整相出力ＡＬを推定する。再推定手段１６は、時間、方位および整相出力の三次元グラフ上に推定結果のスペクトルを出力する（ステップＳ３０８）。ステップＳ３０８の後、信号分析装置１ｂは、次のスナップショットの受波信号に対して、ステップＳ３０１〜Ｓ３０８の処理を繰り返す。
【０１３２】
図１４は、本実施の形態３において、Ｓ個のターゲットからの音を等間隔で直線に配置したＲ個のセンサで受波するモデルを示す図である。図１４に示すモデルの場合を考える。このとき、ｒ番目のセンサのフーリエ変換後のある周波数ｆｎにおける受波信号Ｘｒ（ｆｎ）、｛ｒ＝１，２，・・・，Ｒ｝は、以下となる。
【０１３３】
【数５７】
2020085656.tif 000059
ここで、ａｒｓはｒ番目のセンサで受波したｓ番目のターゲットの振幅、ｄはセンサ間隔、θｒｓはｒ番目のセンサからみたｓ番目のターゲットの到来方向である。φｓはｓ番目のターゲットの初期位相、ｃは音速、ｎｏｉｓｅはｒ番目のセンサで受波した雑音である。同様に得られたＲ個のセンサ出力
【０１３４】
【数５８】
2020085656.tif 000060
を用いて、θｐ方向へ整相した出力Ｙ（θｐ，ｆｎ）は次式で表される。
【０１３５】
【数５９】
2020085656.tif 000061
ここで、ｂ（θｐ，ｆｎ）は各センサの受波信号をθｐ方向へ整相するステアリングベクトルであり、
【０１３６】
【数６０】
2020085656.tif 000062
と表される。出力Ｙ（θｐ，ｆｎ）は、θｐがセンサアレイからみたｓ番目のターゲットの到来方向θｓに一致した場合に大きくなるので、Ｙ（θｐ，ｆｎ）のピーク方向が推定方位となり、ピークの数が推定ターゲット数となる。Ｘ（ｆｎ）は，ステアリングベクトルを用いて以下のように表すことができる。
【０１３７】
【数６１】
2020085656.tif 000063
ここで、
【０１３８】
【数６２】
2020085656.tif 000064
であり、
【０１３９】
【数６３】
2020085656.tif 000065
である。式（１）と式（１８）を比較すると、ステアリングベクトルｂが基底関数に対応していることがわかる。また、周波数解析における１スナップショットの受波信号
【０１４０】
【数６４】
2020085656.tif 000066
は、推定手段２２および再推定手段１６によって行われる信号到来方位推定では、各センサの受波信号
【０１４１】
【数６５】
2020085656.tif 000067
に対応していることがわかる。以上より、基底関数セットＰは、次式に示すような整相方向θｍを任意に設定可能なステアリングベクトルのセットＢに置き換えられる。
【０１４２】
【数６６】
2020085656.tif 000068
【０１４３】
このようにして、ＦＯＳを信号の到来方位推定に適用することができる。実施の形態１で図８を参照して説明したステップＳ１０１〜Ｓ１０８と図１３に示したステップＳ３０１〜Ｓ３０８の動作原理は同じである。本実施の形態３において、実施の形態１と同様に、実施の形態２を適用することもできる。
【０１４４】
従来の整相方向の分解能は（π／Ｍ）でありセンサ数に依存する。そのため、詳細な方位推定を行う場合は、多くのセンサを必要とする。センサ数以上に分析対象を分割すると、ステアリングベクトル同士の直交関係を維持できず、整相方位以外に対して漏れこみが発生する。整相方位に隣接した方位から成分が漏れ込んでくるため、整相方位に関する成分と漏れ込んできた成分との区別ができず分解能は向上しない。これに対して、本実施の形態３の信号分析装置１ｂでは、ＦＯＳ処理でステアリングベクトルを残差計算と直交化するので、センサ数を増やすことなく、方位推定の精度向上が期待できる。
【０１４５】
なお、実施の形態１〜３では、ＦＯＳを用いた信号分析方法について説明したが、ＦＯＳの代わりにＯＭＰ（Orthogonal Matching Pursuit）を用いてもよい。ＯＭＰはＦＯＳに近似した方法である。ＯＭＰについて簡単に説明する。はじめに、ＯＭＰは、残差信号（捜索１回目は受波信号）との相関が最も高い基底関数を基底関数セットから選択する。続いて、ＯＭＰは、選択後に既に選択済みである基底関数リストと直交化を行い、直交化した基底関数リストと受波信号とから残差信号を計算する。基底関数リストを選択した結果を、上述した実施の形態１〜３に、ＦＯＳの場合と同様に適用することができる。ループの終了判定はＦＯＳと同様に残差信号から行う。
【０１４６】
また、実施の形態１〜３では、周波数解析のために基底関数の一例として正弦関数および余弦関数を用いて説明したが、正弦関数および余弦関数に限らずウェーブレット関数等の他の基底関数を用いてもよい。また、実施の形態１〜３では、説明のために、振幅、信号到来方位および整相出力を推定する手段として、推定手段２２および再推定手段１６を説明したが、推定手段２２が再推定手段１６の機能を備えていてもよい。
【符号の説明】
【０１４７】
１、１ａ、１ｂ  信号分析装置
１１  記憶部
１２  制御部
１３  保持手段
１４  可能性値算出手段
１５  追加判定手段
１６  再推定手段
２０  ＦＯＳ処理部
２１  基底関数選択手段
２２  推定手段
２３  残差計算手段
２４  終了判定手段
３１  除外判定手段
４０  ＦＦＴ部
５０  センサアレイ
５１  ＣＰＵ
５２  メモリ
１００  信号分析装置
１０１  記憶部
１０２  制御部。
【図１】
2020085656.tif 000069
【図２】
2020085656.tif 000070
【図３】
2020085656.tif 000071
【図４】
2020085656.tif 000072
【図５】
2020085656.tif 000073
【図６】
2020085656.tif 000074
【図７】
2020085656.tif 000075
【図８】
2020085656.tif 000076
【図９】
2020085656.tif 000077
【図１０】
2020085656.tif 000078
【図１１】
2020085656.tif 000079
【図１２】
2020085656.tif 000080
【図１３】
2020085656.tif 000081
【図１４】
2020085656.tif 000082