"信号源が異なる複数の入射信号が混信している受信信号を白色化する白色化手段と、上記白色化手段により白色化された受信信号におけるピーク間の時間差を変数とする上記受信信号の時間差相互相関行列を生成し、所定の範囲内の任意の時間差を上記時間差相互相関行列に順次代入して、各時間差の相互相関値を算出する相互相関値算出手段と、上記相互相関値算出手段により算出された各時間差の相互相関値の中で、隣接している時間差の相互相関値より大きく、かつ、所定の閾値より大きい相互相関値を検出して、その検出した相互相関値に対応する時間差を特定するピーク相関値検出手段と、上記ピーク相関値検出手段により特定された複数の時間差の中で、同じ周期で現れる時間差を同一のグループに分類し、各グループに属している時間差の周期を推定する周期推定手段と、上記入射信号の周期が上記周期推定手段により推定された周期と一致すると仮定して、上記白色化手段により白色化された受信信号に混信している複数の入射信号の信号分離を行う信号分離手段とを備えた信号分離装置。信号分離手段は、独立指標として時間差相互相関値を用いる独立成分分析であるTDを実施するものであり、周期推定手段により推定された周期における受信信号の時間差相互相関値を最小化することで、上記受信信号に混信している複数の入射信号の信号分離を行うことを特徴とする請求項1記載の信号分離装置。周期推定手段は、あるグループに属している時間差の周期が、他のグループに属している時間差の周期の整数倍である場合、どちらか一方のグループに属している時間差の周期だけを推定結果として信号分離手段に出力することを特徴とする請求項1または請求項2記載の信号分離装置。周期推定手段は、ピーク相関値検出手段により特定された複数の時間差の中から任意の時間差を選択し、上記時間差を推定済みの周期で除算した場合の余りが所定値より小さければ、上記時間差を上記周期で現れる時間差のグループに分類し、上記余りが所定値より大きければ、その選択した時間差を上記周期以外の周期で現れる時間差のグループに分類することを特徴とする請求項1または請求項2記載の信号分離装置。周期推定手段は、ピーク相関値検出手段により特定された複数の時間差をPRI変換することで、上記複数の時間差の中で、同じ周期で現れる時間差を特定して、同じ周期で現れる時間差を同一のグループに分類することを特徴とする請求項1または請求項2記載の信号分離装置。白色化手段が、信号源が異なる複数の入射信号が混信している受信信号を白色化する白色化処理ステップと、相互相関値算出手段が、上記白色化処理ステップで白色化された受信信号におけるピーク間の時間差を変数とする上記受信信号の時間差相互相関行列を生成し、所定の範囲内の任意の時間差を上記時間差相互相関行列に順次代入して、各時間差の相互相関値を算出する相互相関値算出処理ステップと、ピーク相関値検出手段が、上記相互相関値算出処理ステップで算出された各時間差の相互相関値の中で、隣接している時間差の相互相関値より大きく、かつ、所定の閾値より大きい相互相関値を検出して、その検出した相互相関値に対応する時間差を特定するピーク相関値検出処理ステップと、周期推定手段が、上記ピーク相関値検出処理ステップで特定された複数の時間差の中で、同じ周期で現れる時間差を同一のグループに分類し、各グループに属している時間差の周期を推定する周期推定処理ステップと、信号分離手段が、上記入射信号の周期が上記周期推定処理ステップで推定された周期と一致すると仮定して、上記白色化処理ステップで白色化された受信信号に混信している複数の入射信号の信号分離を行う信号分離処理ステップとを備えた信号分離方法。従来の信号分離装置は以上のように構成されているので、ランダムに複数の時間差τを選択する方法では、最適な時間差τが必ずしも選択されない。入射信号の時間差自己相関値が最大となる0以外の時間差τを選択する方法では、入射信号の周期を事前に知っている必要があり、事前に入射信号の周期を知らなければ、時間差τを選択することができないなどの課題があった。この発明は上記のような課題を解決するためになされたもので、事前に入射信号の周期を知っていなくても、適正な時間差τを自動的に選択して、受信信号に混信している複数の信号信号源が異なる複数の入射信号を信号分離することができる信号分離装置及び信号分離方法を得ることを目的とする。この発明に係る信号分離装置は、信号源が異なる複数の入射信号が混信している受信信号を白色化する白色化手段と、白色化手段により白色化された受信信号におけるピーク間の時間差を変数とする上記受信信号の時間差相互相関行列を生成し、所定の範囲内の任意の時間差を上記時間差相互相関行列に順次代入して、各時間差の相互相関値を算出する相互相関値算出手段と、相互相関値算出手段により算出された各時間差の相互相関値の中で、隣接している時間差の相互相関値より大きく、かつ、所定の閾値より大きい相互相関値を検出して、その検出した相互相関値に対応する時間差を特定するピーク相関値検出手段と、ピーク相関値検出手段により特定された複数の時間差の中で、同じ周期で現れる時間差を同一のグループに分類し、各グループに属している時間差の周期を推定する周期推定手段とを設け、信号分離手段が、入射信号の周期が周期推定手段により推定された周期と一致すると仮定して、白色化手段により白色化された受信信号に混信している複数の入射信号の信号分離を行うようにしたものである。この発明によれば、信号源が異なる複数の入射信号が混信している受信信号を白色化する白色化手段と、白色化手段により白色化された受信信号におけるピーク間の時間差を変数とする上記受信信号の時間差相互相関行列を生成し、所定の範囲内の任意の時間差を上記時間差相互相関行列に順次代入して、各時間差の相互相関値を算出する相互相関値算出手段と、相互相関値算出手段により算出された各時間差の相互相関値の中で、隣接している時間差の相互相関値より大きく、かつ、所定の閾値より大きい相互相関値を検出して、その検出した相互相関値に対応する時間差を特定するピーク相関値検出手段と、ピーク相関値検出手段により特定された複数の時間差の中で、同じ周期で現れる時間差を同一のグループに分類し、各グループに属している時間差の周期を推定する周期推定手段とを設け、信号分離手段が、入射信号の周期が周期推定手段により推定された周期と一致すると仮定して、白色化手段により白色化された受信信号に混信している複数の入射信号の信号分離を行うように構成したので、事前に入射信号の周期を知っていなくても、適正な時間差を自動的に選択して、受信信号に混信している複数の入射信号を信号分離することができる効果がある。実施の形態1.図1はこの発明の実施の形態1による信号分離装置を示す構成図である。図1において、白色化部1は複数のアンテナにより受信された受信信号x信号源が異なる複数の入射信号sが混信している受信信号xを白色化し、白色化後の受信信号である白色化受信信号zを時間差相互相関計算部2及び信号分離部5に出力する処理を実施する。なお、白色化部1は白色化手段を構成している。時間差相互相関計算部2は白色化部1から出力された白色化受信信号zにおけるピーク間の時間差τを変数とする白色化受信信号zの時間差相互相関行列Rzτを生成し、所定の範囲内の任意の時間差τを時間差相互相関行列Rzτに順次代入して、各時間差τの相互相関値時間差相互相関の二乗和評価関数値Jτを算出する処理を実施する。なお、時間差相互相関計算部2は相互相関値算出手段を構成している。ピーク列抽出部3は時間差相互相関計算部2により算出された各時間差τの二乗和評価関数値Jτの中で、隣接している時間差τ-1,τ1の二乗和評価関数値Jτ-1,Jτ1より大きく、かつ、所定の閾値thresholdより大きい二乗和評価関数値Jτを検出して、その二乗和評価関数値Jτに対応する時間差τを特定する処理を実施する。なお、ピーク列抽出部3はピーク相関値検出手段を構成している。周期推定部4はピーク列抽出部3により特定された複数の時間差τの中で、同じ周期で現れる時間差τを同一のグループに分類し、各グループに属している時間差τの周期τpを推定する処理を実施する。p=1,2,,Pであり、Pは分類されたグループ数である。なお、周期推定部4は周期推定手段を構成している。信号分離部5は受信信号xに混信している複数の入射信号sの周期が周期推定部4により推定された周期τpと一致すると仮定して、白色化部1から出力された白色化受信信号zに混信している複数の入射信号sの信号分離を行う。即ち、信号分離部5は独立指標として時間差相互相関値を用いる独立成分分析であるTDを実施するものであり、周期推定部4により推定された周期τpにおける白色化受信信号zの時間差相互相関値を最小化することで、その白色化受信信号zに混信している複数の入射信号sの信号分離を行う。なお、信号分離部5は信号分離手段を構成している。図1の例では、信号分離装置の構成要素である白色化部1、時間差相互相関計算部2、ピーク列抽出部3、周期推定部4及び信号分離部5のそれぞれが専用のハードウェア例えば、CPUを実装している半導体集積回路、あるいは、ワンチップマイコンなどで構成されているものを想定しているが、信号分離装置がコンピュータで構成されていてもよい。信号分離装置がコンピュータで構成される場合、白色化部1、時間差相互相関計算部2、ピーク列抽出部3、周期推定部4及び信号分離部5の処理内容を記述しているプログラムをコンピュータのメモリに格納し、当該コンピュータのCPUが当該メモリに格納されているプログラムを実行するようにすればよい。図2はこの発明の実施の形態1による信号分離装置の処理内容を示すフローチャートである。図3はこの発明の実施の形態1による信号分離装置の周波数推定部4を示す構成図である。図3において、第1周期設定部11はピーク列抽出部3により特定された複数の時間差τのうち、第1番目の時間差τを第1番目の周期τ1に設定する処理を実施する。周期格納部12は設定済みの周期τpp=1,2,,Pを格納する記録媒体である。次ピーク時間差選択部13はピーク列抽出部3により特定された複数の時間差τのうち、第n番目の時間差τを選択する処理を実施する。ピーク時間差/周期除算部14は次ピーク時間差選択部13により選択された第n番目の時間差τを周期格納部12により格納されている設定済みの周期τpでそれぞれ除算し第n番目の時間差τ/τ1,第n番目の時間差τ/τ2,,第n番目の時間差τ/τP、その除算したときの余りrpp=1,2,,Pを同一列判定部15に出力する処理を実施する。同一列判定部15はピーク時間差/周期除算部14から出力された余りrpが零、あるいは、余りrpが十分に小さい余りが所定値εより小さい場合、第n番目の時間差τが設定済みの周期τpの列に属していると判定して、第n番目の時間差τを設定済みの周期τpで現れる時間差のグループに分類する。一方、ピーク時間差/周期除算部14から出力された全ての余りrpが十分に大きい余りが所定値εより大きい場合、第n番目の時間差τが何れの周期τpの列にも属していないと判定して、第n番目の時間差τを第P1番目の周期τp1に設定する処理を実施する。処理終了判定部16はピーク列抽出部3により特定された全ての時間差τについて、同一列判定部15の判定処理が終了しているか否かを判定し、判定処理が終了していなければ、現時点で未だ選択していない時間差τの選択を次ピーク時間差選択部13に指示する処理を実施する。次に動作について説明する。図1の信号分離装置は、概略的には、時間差τの時間差相互相関行列Rzτの非対角項である白色化受信信号zの時間差相互相関値が大きい時間差τを選択し、その時間差τをTDに設定するものである。ただし、理想的な周期信号の場合、時間差相互相関行列Rzτには、各入射信号sの周期でピークが立つため、それぞれのピーク列から1つずつ時間差τを選択するようにする。即ち、所定の範囲内にある各々の時間差τ毎に、白色化受信信号zの時間差相互相関値を計算し、その計算結果の中で、ピークを有する時間差相互相関値に対応する時間差τを選択する。複数の時間差τを選択すると、それらの時間差τの列から各ピーク列に対応する周期τpを推定する。以下、図1の信号分離装置の処理内容を具体的に説明する。信号受信部図示せぬは、複数のアンテナが混合信号信号源が異なる複数の入射信号sが混合している信号を受信すると図12を参照、その混合信号に対するA/D変換処理等の受信処理を実施して、ディジタルの混合信号である受信信号xを白色化部1に出力する。白色化部1は、ディジタルの受信信号xを受けると、下記の式11に示すように、その受信信号xを白色化し、白色化後の受信信号である白色化受信信号zを時間差相互相関計算部2及び信号分離部5に出力するステップST1。受信信号xの白色化の処理自体は公知の技術であるため詳細な説明を省略する。ただし、Mは白色化用の行列である。なお、白色化部1は、信号分離部5により分離される入射信号sの振幅が不定であるため、受信信号xを事前に白色化しており、白色化後の受信信号である白色化受信信号zに対して、後述する独立成分分析であるTDを実施する場合、白色化していない受信信号xに対してTDを実施する場合よりも、TDにおける学習が高速化されることが期待される。時間差相互相関計算部2は、白色化部1から白色化受信信号zを受けると、その白色化受信信号zにおけるピーク間の時間差τを変数とする白色化受信信号zの時間差相互相関行列Rzτを生成するステップST2。例えば、入射信号数が2である場合、時間差τにおける白色化受信信号zの時間差相互相関行列Rzτは、下記の式12のように表される。式12の1,2要素と2,1要素に時間差相互相関が現れているため、これらを最小化するような分離行列Wを計算する。これは式12の対角化を意味している。なお、Eは期待値または平均値を意味する数学記号であり、Hは行列の転置を意味し、は共役複素数を意味する。時間差相互相関計算部2は、白色化受信信号zの時間差相互相関行列Rzτを生成すると、下記の式13に示すように、所定の範囲内の任意の時間差τを時間差相互相関行列Rzτに順次代入することで、各時間差τの相互相関値として、時間差相互相関の二乗和評価関数値Jτを算出するステップST3。ただし、off.は行列の非対角項を意味する。また、Rzj1j2チルダは、行列Rzチルダの第j1,j2要素であることを意味し、j1≠j2は非対角項を意味している。ここでは、所定の範囲内の任意の時間差τを時間差相互相関行列Rzτに代入しているが、所定の範囲は、例えば、受信信号xの長さの半分程度が考えられる。時間差相互相関行列Rzτに代入する時間差τは離散的な値であり、例えば、受信信号xの長さが2000であれば、0~1000の間の離散的な値である。図4は時間差相互相関計算部2による二乗和評価関数値Jτの算出例を示す説明図であり、各時間差τの二乗和評価関数値Jτが算出されている。入射信号sの周期に対応する時間差τの二乗和評価関数値Jτは大きなピークとなるが、入射信号sの周期に対応していない時間差τの二乗和評価関数値Jτにも小さなピークが立つことがあり、その小さなピークと本来の大きなピークを区別する必要がある。そこで、ピーク列抽出部3は、時間差相互相関計算部2が各時間差τの二乗和評価関数値Jτを算出すると、下記の式14に示すように、各時間差τの二乗和評価関数値Jτの中で、隣接している時間差τ-1,τ1の二乗和評価関数値Jτ-1,Jτ1より大きく、かつ、所定の閾値thresholdより大きい二乗和評価関数値Jτを検出して、その二乗和評価関数値Jτに対応する時間差τを特定するステップST4。閾値thresholdの設定方法は様々な方法が考えられるが、ここでは、下記の式15に示すように、隣接している時間差τi-1,τi1の二乗和評価関数値Jτi-1,Jτi1より大きい全ての二乗和評価関数値Jτiの平均値としている。ただし、mean.は平均値の計算を意味している。図5はピーク列抽出部3により検出された二乗和評価関数値Jτに対応する時間差τを示す説明図である。図5の例では、14個の二乗和評価関数値Jτが検出されて、14個の時間差τが特定されている。周期推定部4は、ピーク列抽出部3が複数の二乗和評価関数値Jτを検出して、複数の二乗和評価関数値Jτに対応する時間差τを特定すると、複数の時間差τの中で、同じ周期で現れる時間差τを同一のグループに分類し、各グループに属している時間差τの周期τpを推定するステップST5。p=1,2,,Pであり、Pは分類されたグループ数である。以下、周期推定部4の処理内容を具体的に説明する。図6は周期推定部4による周期推定方法の概要を示す説明図である。周期推定部4では、周期的なピークの時間差は、周期の整数倍の時間差に存在する性質を利用し、最初のピークの時間差から順番に時間差の割り算を行い、余りが小さければ、同一周期のピーク列に存在し、余りが大きければ、異なる周期のピーク列に属するとみなすことで、それぞれの周期を検出する方法を採用している。周期推定部4の第1周期設定部11は、ピーク列抽出部3が複数の二乗和評価関数値Jτを検出して、複数の二乗和評価関数値Jτに対応する時間差τを特定すると、複数の時間差τの中で、第1番目の時間差τを第1番目の周期τ1に設定する。なお、第1周期設定部11により設定された第1番目の周期τ1は周期格納部12により格納される。次ピーク時間差選択部13は、ピーク列抽出部3により特定された複数の時間差τのうち、第n番目の時間差τを選択して、その時間差τをピーク時間差/周期除算部14に出力する。n=1,2,,Nであり、Nはピーク列抽出部3により特定された時間差τの個数である。ピーク時間差/周期除算部14は、次ピーク時間差選択部13から第n番目の時間差τを受けると、第n番目の時間差τを周期の候補τcとして、その周期の候補τcを周期格納部12により格納されている設定済みの周期τpで除算し、その除算したときの余りrpを同一列判定部15に出力する。即ち、ピーク時間差/周期除算部14は、modτc,τpを演算することにより、除算時の余りrpを同一列判定部15に出力する。なお、moda,bはaをbで割った余りを返す関数である。初期の段階では、第1番目の周期τ1だけが周期格納部12に格納されているため、modτc,τ1だけを演算して余りr1を出力するが、その後、P-1個の周期τ2,τ3,,τPが設定されて、P個の周期τ1,τ2,,τPが周期格納部12に格納されていれば、P個のmodτc,τ1,modτc,τ2,,modτc,τPを演算して、P個の余りr1,r2,,rPを出力する。同一列判定部15は、ピーク時間差/周期除算部14から余りrpを受けると、その余りrpと所定値εを比較する。所定値εは周期τpより小さい零に近い値であり、例えば、2や3などの値が設定される。なお、ピーク時間差/周期除算部14からP個の余りr1,r2,,rPを受けていれば、P個の余りr1,r2,,rPと所定値εをそれぞれ比較する。同一列判定部15は、余りrpが零である場合、あるいは、余りrpが所定値εより小さい場合0<rp<ε、第n番目の時間差τが設定済みの周期τpの列に属していると判定して、第n番目の時間差τを設定済みの周期τpで現れる時間差のグループに分類する。例えば、ピーク時間差/周期除算部14からP個の余りr1,r2,,rPを受けているとき、余りr2だけが所定値εより小さければ、第n番目の時間差τを設定済みの周期τ2で現れる時間差のグループに分類する。同一列判定部15は、ピーク時間差/周期除算部14から出力された全ての余りrpが所定値εより大きい場合P個の余りr1,r2,,rPを受けていれば、余りr1,r2,,rPの全てが所定値εより大きい場合、第n番目の時間差τが何れの周期τpの列にも属していないと判定して、新たに第n番目の時間差τを第P1番目の周期τP1に設定する。なお、同一列判定部15により設定された第P1番目の周期τP1は周期格納部12により格納される。処理終了判定部16は、ピーク列抽出部3により特定された全ての時間差τについて、同一列判定部15の判定処理が終了しているか否かを判定し、判定処理が終了していなければ、現時点で未だ選択していない時間差τの選択を次ピーク時間差選択部13に指示する。一方、判定処理が終了していれば、設定済みの周期τpの出力を周期格納部12に指示する。なお、M個の周期τpが設定されていれば、周期格納部12からM個の周期τ1,τ2,,τMが信号分離部5に出力される。ここでは、ピーク時間差/周期除算部14が、modτc,τpを演算するものを示したが、modτc,τpの関数から返される余りrpは正の値だけである。しかし、ピーク出現誤差は正負両方に出るため、modτc,τpを演算する他に、τp-modτc,τpも演算して、双方の演算から得られる余りrpを所定値εと比較する検定を行うようにしてもよい。なお、検定式は、下記の式16の通りである。また、ここでは、処理終了判定部16が、ピーク列抽出部3により特定された全ての時間差τについて、同一列判定部15の判定処理が終了しているか否かを判定するものを示したが、所望の数だけ周期τpが設定されたか否かを判定し、所望の数だけ周期τpが設定されていれば、ピーク列抽出部3により特定された全ての時間差τについての判定処理が終了していなくても、設定済みの周期τpの出力を周期格納部12に指示するようにしてもよい。信号分離部5は、周期推定部4からM個の周期τ1,τ2,,τMを受けると、受信信号xに混信している各入射信号sの周期が、M個の周期τ1,τ2,,τMのいずれかと一致すると仮定して、白色化部1から出力された白色化受信信号zに混信している複数の入射信号sの信号分離を行うステップST6。即ち、信号分離部5は、TDにおける時間差パラメータとして、周期推定部4から出力されたM個の周期τ1,τ2,,τMを順番に設定し、その設定した周期τpにおける白色化受信信号zの時間差相互相関値を最小化することで、その白色化受信信号zから当該周期τpを有する入射信号sを信号分離する。ここでは、周期推定部4がM個の周期τ1,τ2,,τMを信号分離部5に出力し、信号分離部5がTDにおける時間差パラメータとして、M個の周期τ1,τ2,,τMを順番に設定するものを示したが、TDの処理負荷を軽減するために、M個の周期τ1,τ2,,τMの中の或る周期が、他の周期の整数倍である場合、周期推定部4が、どちらか一方の周期だけを信号分離部5に出力して、他方の周期を信号分離部5に出力しないようにしてもよい。例えば、周期τ4が周期τ1の整数倍である場合、周期τ4又は周期τ1の一方だけを信号分離部5に出力して、他方の周期を信号分離部5に出力しないようにする。この場合、TDにおける時間差パラメータとして、周期τ1が設定されれば、周期τ4は設定されず、周期τ4が設定されれば、周期τ1は設定されない。以下、図1の信号分離装置による入射信号の信号分離のシミュレーション結果について説明する。ここでは、周期が30の入射信号と、周期が74の入射信号と、周期が77の入射信号との3波が混信している受信信号周期性ガウス信号が受信される場合のシミュレーション結果について説明する。図7は3波の入射信号の時間差自己相関を示す説明図である。図7より明らかなように、各入射信号は周期的にピークを有している。ただし、入射信号の周期は事前に分からないため未知である。図8は白色化受信信号の時間差相互相関を示す説明図であり、図9はTDにおける時間差パラメータとして設定される時間差ττ=30,74,77を示す説明図である。また、図10は実施の形態1による信号分離装置及び従来の信号分離装置における入射信号の分離性能のシミュレーション結果を示す説明図である。図8から明らかなように、白色化受信信号zの時間差相互相関には、各入射信号sの周期でピークが現れる。このシミュレーションでは、図9に示すように、TDにおける時間差パラメータとして、時間差ττ=30,74,77が設定されており、3波の入射信号の周期30,74,77と一致している。この実施の形態1による信号分離装置によれば、図10に示すように、従来の信号分離装置時間差τをランダムに選択する方法を用いる信号分離装置よりも、入射信号の分離性能が約20dB改善している。以上で明らかなように、この実施の形態1によれば、信号源が異なる複数の入射信号sが混信している受信信号xを白色化する白色化部1と、白色化部1により白色化された受信信号である白色化受信信号zにおけるピーク間の時間差τを変数とする白色化受信信号zの時間差相互相関行列Rzτを生成し、所定の範囲内の任意の時間差τを時間差相互相関行列Rzτに順次代入して、各時間差τの相互相関値である二乗和評価関数値Jτを算出する時間差相互相関計算部2と、時間差相互相関計算部2により算出された各時間差τの二乗和評価関数値Jτの中で、隣接している時間差τ-1,τ1の二乗和評価関数値Jτ-1,Jτ1より大きく、かつ、所定の閾値thresholdより大きい二乗和評価関数値Jτを検出して、その二乗和評価関数値Jτに対応する時間差τを特定するピーク列抽出部3と、ピーク列抽出部3により特定された複数の時間差τの中で、同じ周期で現れる時間差τを同一のグループに分類し、各グループに属している時間差τの周期τpを推定する周期推定部4とを設け、信号分離部5が、入射信号sの周期が周期推定部4により推定された周期τpと一致すると仮定して、その白色化受信信号zに混信している複数の入射信号sの信号分離を行うように構成したので、事前に入射信号sの周期を知っていなくても、適正な時間差τを自動的に選択して、受信信号xに混信している複数の入射信号sを信号分離することができる効果を奏する。また、この実施の形態1によれば、あるグループに属している時間差τの周期が、他のグループに属している時間差τの周期の整数倍である場合、周期推定部4が、どちらか一方のグループに属している時間差τの周期だけを推定結果として信号分離部5に出力するように構成したので、無駄なTDの実施を省くことができるようになり、その結果、信号分離性能の劣化を招くことなく、信号分離部5におけるTDの処理負荷を軽減することができる効果を奏する。また、この実施の形態1によれば、周期推定部4が、ピーク列抽出部3により特定された複数の時間差τの中から任意の時間差τを選択し、その時間差τを設定済みの周期τpで除算した場合の余りrpが所定値εより小さければ、その時間差τを周期τpで現れる時間差τのグループに分類し、その余りrpが所定値εより大きければ、その選択した時間差τを周期τp以外の周期で現れる時間差τのグループに分類するように構成したので、同じ周期で現れる時間差τを同一のグループに分類することができる効果を奏する。実施の形態2.上記実施の形態1では、周期推定部4が、ピーク列抽出部3により特定された複数の時間差τの中から任意の時間差τを選択し、その時間差τを設定済みの周期τpで除算した場合の余りrpが所定値εより小さければ、その時間差τを周期τpで現れる時間差τのグループに分類し、その余りrpが所定値εより大きければ、その選択した時間差τを周期τp以外の周期で現れる時間差τのグループに分類するものを示したが、図11に示すように、周期推定部4がPRI変換部4aを内蔵し、PRI変換部4aが、ピーク列抽出部3により特定された複数の時間差τをPRI変換することで、複数の時間差τの中で、同じ周期で現れる時間差τを特定して、その時間差τのグループ分け同じ周期で現れる時間差τは同じグループに分類を行うようにしてもよく、上記実施の形態1と同様の効果を奏することができる。なお、PRI変換は、公知の技術であるため詳細な説明は省略するが例えば、特公昭62-26603号公報を参照、異なる周期を持つピーク列が混在しているピーク列から、それぞれのピーク列を検出して、それぞれのピーク列の周期を検出するアルゴリズムである。なお、本願発明はその発明の範囲内において、各実施の形態の自由な組み合わせ、あるいは各実施の形態の任意の構成要素の変形、もしくは各実施の形態において任意の構成要素の省略が可能である。"