(57)【要約】      （修正有）【課題】２つの標本化データの相互相関のピーク位置を求めるピーク位置探索方法およびピーク位置探索装置等に関し、２つの標本化データの相互相関のピーク位置を高精度に求める。【解決手段】２つの標本化データをそれぞれ離散フーリエ変換し２つのスペクトルを得る。次に一方を複素共役スペクトルとして他方のスペクトルとの積の逆離散フーリエ変換を行って相互相関関数を求め、この相互相関関数のΔｔごとの標本値から相互相関のピーク位置を求める。【選択図】図６
2008282268.tif 000002
【特許請求の範囲】
【請求項１】
２つの標本化データの相互相関のピーク位置を求めるピーク位置探索方法において、
所定の標本化間隔
【数１】
2008282268.tif 000141
で標本化された２つの標本化データ
【数２】
2008282268.tif 000142
を取得するデータ取得ステップと、
前記データ取得ステップで取得した２つの標本化データそれぞれに
【数３】
2008282268.tif 000143
個のゼロ点を付加した２つの標本化データ
【数４】
2008282268.tif 000144
のそれぞれについて離散フーリエ変換することにより、２つのスペクトル
【数５】
2008282268.tif 000145
を算出するスペクトル算出ステップと、
前記スペクトル算出ステップで算出された２つのスペクトルのうちの一方のスペクトル
【数６】
2008282268.tif 000146
の複素共役と他方のスペクトル
【数７】
2008282268.tif 000147
との積の逆離散フーリエ変換をすることにより、相互相関関数
【数８】
2008282268.tif 000148
を求める相互相関関数算出ステップと、
前記標本化間隔
【数９】
2008282268.tif 000149
ごとの標本化点
【数１０】
2008282268.tif 000150
の途中の点を含む任意点
【数１１】
2008282268.tif 000151
の相互相関関数
【数１２】
2008282268.tif 000152
を表わす式
【数１３】
2008282268.tif 000153
但し、
【数１４】
2008282268.tif 000154
は内挿関数であって、
【数１５】
2008282268.tif 000155
で表わされる、
に従って、ピーク位置を求めるピーク位置探索ステップとを有することを特徴とするピーク位置探索方法。
【請求項２】
前記ピーク位置探索ステップが、前記相互相関関数算出ステップで求められた標本化間隔
【数１６】
2008282268.tif 000156
ごとの相互相関関数から求められる、誤差を含む第１のピーク位置近傍の第１区間を該標本化間隔
【数１７】
2008282268.tif 000157
よりも細かい間隔
【数１８】
2008282268.tif 000158
で細分化した各細分点の相互相関関数を求める第１細分化ステップと、該第１細分化ステップで求められた間隔
【数１９】
2008282268.tif 000159
ごとの相互相関関数から求められる、誤差を含む第２のピーク位置近傍の、前記第１区間よりも狭い第２区間を該間隔
【数２０】
2008282268.tif 000160
よりも細かい間隔
【数２１】
2008282268.tif 000161
で細分化した各細分点の相互相関関数を求める第２細分化ステップと、以下同様にして所望精度のピーク位置が探索されるまで、第ｎ−１細分化ステップ（ｎ＝２，３，…）で求められた間隔
【数２２】
2008282268.tif 000162
ごとの相互相関関数から求められる、誤差を含む第ｎ−１のピーク位置近傍の第ｎ−１区間よりも狭い第ｎ区間を該間隔
【数２３】
2008282268.tif 000163
よりも細かい間隔
【数２４】
2008282268.tif 000164
で細分化した各細分点の相互相関関数を求める第ｎ細分化ステップとを有することを特徴とする請求項１記載のピーク位置探索方法。
【請求項３】
前記標本化間隔
【数２５】
2008282268.tif 000165
が標本化時間間隔
【数２６】
2008282268.tif 000166
であり、前記データ取得ステップは、該標本化時間間隔
【数２７】
2008282268.tif 000167
で標本化された２つの標本化時系列データ
【数２８】
2008282268.tif 000168
但し、
【数２９】
2008282268.tif 000169
での信号値
【数３０】
2008282268.tif 000170
を
【数３１】
2008282268.tif 000171
などと略記した、
を取得するステップであることを特徴とする請求項１又は２記載のピーク位置探索方法。
【請求項４】
２つの標本化データの相互相関のピーク位置を求めるピーク位置探索装置において、
所定の標本化間隔
【数３２】
2008282268.tif 000172
で標本化された２つの標本化データ
【数３３】
2008282268.tif 000173
を取得するデータ取得部と、
前記データ取得ステップで取得した２つの標本化データそれぞれに
【数３４】
2008282268.tif 000174
個のゼロ点を付加した２つの標本化データ
【数３５】
2008282268.tif 000175
のそれぞれについて離散フーリエ変換することにより、２つのスペクトル
【数３６】
2008282268.tif 000176
を算出するスペクトル算出部と、
前記スペクトル算出ステップで算出された２つのスペクトルのうちの一方のスペクトル
【数３７】
2008282268.tif 000177
の複素共役と他方のスペクトル
【数３８】
2008282268.tif 000178
との積の逆離散フーリエ変換をすることにより、相互相関関数
【数３９】
2008282268.tif 000179
を求める相互相関関数算出部と、
前記標本化間隔
【数４０】
2008282268.tif 000180
ごとの標本化点
【数４１】
2008282268.tif 000181
の途中の点を含む任意点
【数４２】
2008282268.tif 000182
の相互相関関数
【数４３】
2008282268.tif 000183
を表わす式
【数４４】
2008282268.tif 000184
但し、
【数４５】
2008282268.tif 000185
は内挿関数であって、
【数４６】
2008282268.tif 000186
で表わされる、
に従って、ピーク位置を求めるピーク位置探索部とを備えたことを特徴とするピーク位置探索装置。
【請求項５】
前記ピーク位置探索部は、前記相互相関関数算出ステップで求められた標本化間隔
【数４７】
2008282268.tif 000187
ごとの相互相関関数から求められる、誤差を含む第１のピーク位置近傍の第１区間を該標本化間隔
【数４８】
2008282268.tif 000188
よりも細かい間隔
【数４９】
2008282268.tif 000189
で細分化した各細分点の相互相関関数を求める第１細分化ステップと、該第１細分化ステップで求められた間隔
【数５０】
2008282268.tif 000190
ごとの相互相関関数から求められる、誤差を含む第２のピーク位置近傍の前記第１区間よりも狭い第２区間を該間隔
【数５１】
2008282268.tif 000191
よりも細かい間隔
【数５２】
2008282268.tif 000192
で細分化した各細分点の相互相関関数を求める第２細分化ステップと、以下同様にして所望精度のピーク位置が探索されるまで、第ｎ−１細分化ステップ（ｎ＝２，３，…）で求められた間隔
【数５３】
2008282268.tif 000193
ごとの相互相関関数から求められる、誤差を含む第ｎ−１のピーク位置近傍の第ｎ−１区間よりも狭い第ｎ区間を該間隔
【数５４】
2008282268.tif 000194
よりも細かい間隔
【数５５】
2008282268.tif 000195
で細分化した各細分点の相互相関関数を求める第ｎ細分化ステップとを経ることにより所望精度のピーク位置を求めるものであることを特徴とする請求項４記載のピーク位置探索装置。
【請求項６】
前記標本化間隔
【数５６】
2008282268.tif 000196
が標本化時間間隔
【数５７】
2008282268.tif 000197
であり、前記データ取得部は、該標本化時間間隔
【数５８】
2008282268.tif 000198
で標本化された２つの標本化時系列データ
【数５９】
2008282268.tif 000199
但し、
【数６０】
2008282268.tif 000200
での信号値
【数６１】
2008282268.tif 000201
を
【数６２】
2008282268.tif 000202
などと略記した、
を取得するものであることを特徴とする請求項４又は５記載のピーク位置探索装置。
【請求項７】
前記データ取得部の前段に、アナログ信号を受け取り該アナログ信号を標本化時間間隔
【数６３】
2008282268.tif 000203
で標本化してデジタルの標本化時系列データを生成するアナログ・デジタル変換部を備えたことを特徴とする請求項６記載のピーク位置探索装置。
【請求項８】
前記アナログ・デジタル変換部のさらに前段に、エイリアジングノイズ防止用のアンチエイリアジングフィルタを備えたことを特徴とする請求項７記載のピーク位置探索装置。
【請求項９】
コンピュータ内で実行され該コンピュータを、２つの標本化データの相互相関のピーク位置を求めるピーク位置探索装置として動作させるピーク位置探索プログラムであって、
前記コンピュータを、
所定の標本化間隔
【数６４】
2008282268.tif 000204
で標本化された２つの標本化データ
【数６５】
2008282268.tif 000205
を取得するデータ取得部と、
前記データ取得ステップで取得した２つの標本化データそれぞれに
【数６６】
2008282268.tif 000206
個のゼロ点を付加した２つの標本化データ
【数６７】
2008282268.tif 000207
のそれぞれについて離散フーリエ変換することにより、２つのスペクトル
【数６８】
2008282268.tif 000208
を算出するスペクトル算出部と、
前記スペクトル算出ステップで算出された２つのスペクトルのうちの一方のスペクトル
【数６９】
2008282268.tif 000209
の複素共役と他方のスペクトル
【数７０】
2008282268.tif 000210
との積の逆離散フーリエ変換をすることにより、相互相関関数
【数７１】
2008282268.tif 000211
を求める相互相関関数算出部と、
前記標本化間隔
【数７２】
2008282268.tif 000212
ごとの標本化点
【数７３】
2008282268.tif 000213
の途中の点を含む任意点
【数７４】
2008282268.tif 000214
の相互相関関数
【数７５】
2008282268.tif 000215
を表わす式
【数７６】
2008282268.tif 000216
但し、
【数７７】
2008282268.tif 000217
は内挿関数であって、
【数７８】
2008282268.tif 000218
で表わされる、
に従って、ピーク位置を求めるピーク位置探索部とを備えたピーク位置探索装置として動作させることを特徴とするピーク位置探索プログラム。
【請求項１０】
前記ピーク位置探索部は、前記相互相関関数算出ステップで求められた標本化間隔
【数７９】
2008282268.tif 000219
ごとの相互相関関数から求められる、誤差を含む第１のピーク位置近傍の第１区間を該標本化間隔
【数８０】
2008282268.tif 000220
よりも細かい間隔
【数８１】
2008282268.tif 000221
で細分化した各細分点の相互相関関数を求める第１細分化ステップと、該第１細分化ステップで求められた間隔
【数８２】
2008282268.tif 000222
ごとの相互相関関数から求められる、誤差を含む第２のピーク位置近傍の前記第１区間よりも狭い第２区間を該間隔
【数８３】
2008282268.tif 000223
よりも細かい間隔
【数８４】
2008282268.tif 000224
で細分化した各細分点の相互相関関数を求める第２細分化ステップと、以下同様にして所望精度のピーク位置が探索されるまで、第ｎ−１細分化ステップ（ｎ＝２，３，…）で求められた間隔
【数８５】
2008282268.tif 000225
ごとの相互相関関数から求められる、誤差を含む第ｎ−１ピーク位置近傍の第ｎ−１区間よりも狭い第ｎ区間を該間隔
【数８６】
2008282268.tif 000226
よりも細かい間隔
【数８７】
2008282268.tif 000227
で細分化した各細分点の相互相関関数を求める第ｎ細分化ステップとを経ることにより所望精度のピーク位置を求めるものであることを特徴とする請求項９記載のピーク位置探索プログラム。
【請求項１１】
前記標本化間隔
【数８８】
2008282268.tif 000228
が標本化時間間隔
【数８９】
2008282268.tif 000229
であり、前記データ取得部は、該標本化時間間隔
【数９０】
2008282268.tif 000230
で標本化された２つの標本化時系列データ
【数９１】
2008282268.tif 000231
但し、
【数９２】
2008282268.tif 000232
での信号値
【数９３】
2008282268.tif 000233
を
【数９４】
2008282268.tif 000234
などと略記した、
を取得するものであることを特徴とする請求項９又は１０記載のピーク位置探索プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、２つの標本化データの相互相関のピーク位置を求めるピーク位置探索方法およびピーク位置探索装置、並びに、コンピュータ内で実行されてそのコンピュータをピーク位置探索装置として動作させるピーク位置検索プログラムに関する。
【背景技術】
【０００２】
相互相関法により遅延時間を求めることは、計測の分野で広範囲に利用されている技術である。このことから相互相関法による遅延時間の推定は基本的なテーマであると言えるが、標本化データを用いる場合には、その標本化時間間隔により分解能に制限が出て来る。こういった問題に対し、従来から、相互相関関数のピーク付近に２次関数をフィットさせてピーク位置を推定する方法や、標本化時間間隔を再設定し相互相関関数を計算し直し、ピーク位置を推定する方法が採られてきた（非特許文献１，２参照）。このうち、前者の方法はピーク位置の誤差量が推定しにくい近似的方法の感が拭い切れないし、また、後者の方法は、あくまでも、標本化時間間隔で精度が定まってしまい、精度上の問題を無くそうとすれば無限小の標本化時間間隔を要する、すなわち、無限個の標本化データを必要とすることから、実用上は不可能というジレンマに陥ることになる。他方、周期関数の標本化定理を利用することにより、１周期分の有限個の標本化データから任意の時刻での信号値を正確に求める手法が提案されている（特許文献１）。しかしながらこの特許文献１の技術を相互相関に応用した技術は何ら提案されていない。例えば非特許文献２は、特許文献１が公開されてから５年近く経た後の文献であるにも拘らず、精度向上を目的として標本化時間間隔を狭めることのみに終始している。
【非特許文献１】Ｓ．Ｄ．Ｓｔｅａｒｎｓ  ａｎｄ  Ｒ．Ａ．Ｄａｖｉｄ，“Ｓｉｇａｎａｌ  Ｐｒｏｃｅｓｓｏｎｇ  Ａｌｇｏｒｉｔｈｍｓ，”（Ｐｒｅｎｔｉｃｅ−Ｈａｌｌ，Ｉｎｃ．，ＥｎｇｌｅｗｏｏｄＣｌｉｆｆｓ，Ｎｅｗ  Ｊｅｒｓｅｙ，１９８８）ｐｐ．２１２−２２３
【非特許文献２】牧野康一、横田考俊、山本貢平、岡田恭明、吉久光一、“相互相関法を利用した自動車走行騒音の指向性の測定，”日本音響学会講演論文集  ２００７年３月，ｐｐ．８２７−８２８
【特許文献１】特開２００２−２７８９４８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
本発明は、上記事情に鑑み、２つの標本化データの相互相関のピーク位置を高精度に求めることのできるピーク位置探索方法およびピーク位置探索装置、並びに、コンピュータ内で実行されて、そのコンピュータを、ピーク位置を高精度に求めることのできるピーク位置探索装置として動作させることのできるピーク位置探索プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００４】
上記目的を達成する本発明のピーク位置探索方法は、
２つの標本化データの相互相関のピーク位置を求めるピーク位置探索方法であって、
所定の標本化間隔
【０００５】
【数１】
2008282268.tif 000003
【０００６】
で標本化された２つの標本化データ
【０００７】
【数２】
2008282268.tif 000004
【０００８】
を取得するデータ取得ステップと、
前記データ取得ステップで取得した２つの標本化データそれぞれに
【０００９】
【数３】
2008282268.tif 000005
【００１０】
個のゼロ点を付加した２つの標本化データ
【００１１】
【数４】
2008282268.tif 000006
【００１２】
のそれぞれについて離散フーリエ変換することにより、２つのスペクトル
【００１３】
【数５】
2008282268.tif 000007
【００１４】
を算出するスペクトル算出ステップと、
前記スペクトル算出ステップで算出された２つのスペクトルのうちの一方のスペクトル      
【００１５】
【数６】
2008282268.tif 000008
【００１６】
の複素共役と他方のスペクトル
【００１７】
【数７】
2008282268.tif 000009
【００１８】
との積の逆離散フーリエ変換をすることにより、相互相関関数
【００１９】
【数８】
2008282268.tif 000010
【００２０】
を求める相互相関関数算出ステップと、
前記標本化間隔
【００２１】
【数９】
2008282268.tif 000011
【００２２】
ごとの標本化点
【００２３】
【数１０】
2008282268.tif 000012
【００２４】
の途中の点を含む任意点
【００２５】
【数１１】
2008282268.tif 000013
【００２６】
の相互相関関数
【００２７】
【数１２】
2008282268.tif 000014
【００２８】
を表わす式
【００２９】
【数１３】
2008282268.tif 000015
【００３０】
但し、
【００３１】
【数１４】
2008282268.tif 000016
【００３２】
は内挿関数であって、
【００３３】
【数１５】
2008282268.tif 000017
【００３４】
で表わされる、
に従って、ピーク位置を求めるピーク位置探索ステップとを有することを特徴とする。
【００３５】
本発明のピーク位置探索方法は、いわば、前掲の特許文献１の技術を相互相関のピーク位置探索に応用したものであるが、過去には、特許文献１の技術を相互相関やそれに類似した技術に応用した例は存在しない。本発明者は、特許文献１の技術を相互相関に適用することによってピーク位置を高精度に探索することができるという点に想到し、本発明の完成に至ったのである。
【００３６】
ここで、上記本発明のピーク位置探索方法において、上記ピーク位置探索ステップが、相互相関関数算出ステップで求められた標本化間隔
【００３７】
【数１６】
2008282268.tif 000018
【００３８】
ごとの相互相関関数から求められる、誤差を含む第１のピーク位置近傍の第１区間を該標本化間隔
【００３９】
【数１７】
2008282268.tif 000019
【００４０】
よりも細かい間隔
【００４１】
【数１８】
2008282268.tif 000020
【００４２】
で細分化した各細分点の相互相関関数を求める第１細分化ステップと、その第１細分化ステップで求められた間隔
【００４３】
【数１９】
2008282268.tif 000021
【００４４】
ごとの相互相関関数から求められる、誤差を含む第２のピーク位置近傍の、上記第１区間よりも狭い第２区間を間隔
【００４５】
【数２０】
2008282268.tif 000022
【００４６】
よりも細かい間隔
【００４７】
【数２１】
2008282268.tif 000023
【００４８】
で細分化した各細分点の相互相関関数を求める第２細分化ステップと、以下同様にして所望精度のピーク位置が探索されるまで、第ｎ−１細分化ステップ（ｎ＝２，３，…）で求められた間隔
【００４９】
【数２２】
2008282268.tif 000024
【００５０】
ごとの相互相関関数から求められる、誤差を含むピーク位置近傍の、第ｎ−１区間よりも狭い第ｎ区間を該間隔
【００５１】
【数２３】
2008282268.tif 000025
【００５２】
よりも細かい間隔
【００５３】
【数２４】
2008282268.tif 000026
【００５４】
で細分化した各細分点の相互相関関数を求める第ｎ細分化ステップとを有することが好ましい。
【００５５】
この手法を採用することにより、ピーク位置を所望の精度で求めるための演算量が減り高速化が図られる。
【００５６】
また、本発明のピーク位置探索方法において、上記標本化間隔
【００５７】
【数２５】
2008282268.tif 000027
【００５８】
が標本化時間間隔
【００５９】
【数２６】
2008282268.tif 000028
【００６０】
であり、前記データ取得ステップは、標本化時間間隔
【００６１】
【数２７】
2008282268.tif 000029
【００６２】
で標本化された２つの標本化時系列データ
【００６３】
【数２８】
2008282268.tif 000030
【００６４】
但し、
【００６５】
【数２９】
2008282268.tif 000031
【００６６】
での信号値
【００６７】
【数３０】
2008282268.tif 000032
【００６８】
を
【００６９】
【数３１】
2008282268.tif 000033
【００７０】
などと略記した、
を取得するステップであってもよい。
【００７１】
但し、本発明のピーク位置探索方法は、時間を変数とする標本化データのみでなく、空間位置を変数とする標本化データにも同様に適用することができる。
【００７２】
また、上記目的を達成する本発明のピーク位置探索装置は、２つの標本化データの相互相関のピーク位置を求めるピーク位置探索装置であって、
所定の標本化間隔
【００７３】
【数３２】
2008282268.tif 000034
【００７４】
で標本化された２つの標本化データ
【００７５】
【数３３】
2008282268.tif 000035
【００７６】
を取得するデータ取得部と、
データ取得ステップで取得した２つの標本化データそれぞれに
【００７７】
【数３４】
2008282268.tif 000036
【００７８】
個のゼロ点を付加した２つの標本化データ
【００７９】
【数３５】
2008282268.tif 000037
【００８０】
のそれぞれについて離散フーリエ変換をすることにより、２つのスペクトル
【００８１】
【数３６】
2008282268.tif 000038
【００８２】
を算出するスペクトル算出部と、
スペクトル算出ステップで算出された２つのスペクトルのうちの一方のスペクトル
【００８３】
【数３７】
2008282268.tif 000039
【００８４】
の複素共役と他方のスペクトル
【００８５】
【数３８】
2008282268.tif 000040
【００８６】
との積の逆離散フーリエ変換することにより、相互相関関数
【００８７】
【数３９】
2008282268.tif 000041
【００８８】
を求める相互相関関数算出部と、
標本化間隔
【００８９】
【数４０】
2008282268.tif 000042
【００９０】
ごとの標本化点
【００９１】
【数４１】
2008282268.tif 000043
【００９２】
の途中の点を含む任意点
【００９３】
【数４２】
2008282268.tif 000044
【００９４】
の相互相関関数
【００９５】
【数４３】
2008282268.tif 000045
【００９６】
を表わす式
【００９７】
【数４４】
2008282268.tif 000046
【００９８】
但し、
【００９９】
【数４５】
2008282268.tif 000047
【０１００】
は内挿関数であって、
【０１０１】
【数４６】
2008282268.tif 000048
【０１０２】
で表わされる、
に従って、ピーク位置を求めるピーク位置探索部とを備えたことを特徴とする。
【０１０３】
ここで、本発明のピーク位置探索装置においても、上記ピーク位置探索部が、相互相関関数算出ステップで求められた標本化間隔
【０１０４】
【数４７】
2008282268.tif 000049
【０１０５】
ごとの相互相関関数から求められる、誤差を含む第１のピーク位置近傍の第１区間を該標本化間隔
【０１０６】
【数４８】
2008282268.tif 000050
【０１０７】
よりも細かい間隔
【０１０８】
【数４９】
2008282268.tif 000051
【０１０９】
で細分化した各細分点の相互相関関数を求める第１細分化ステップと、その第１細分化ステップで求められた間隔
【０１１０】
【数５０】
2008282268.tif 000052
【０１１１】
ごとの相互相関関数から求められる、誤差を含む第２のピーク位置近傍の、第１区間よりも狭い第２区間を間隔
【０１１２】
【数５１】
2008282268.tif 000053
【０１１３】
よりも細かい間隔
【０１１４】
【数５２】
2008282268.tif 000054
【０１１５】
で細分化した各細分点の相互相関関数を求める第２細分化ステップと、以下同様にして所望精度のピーク位置が探索されるまで、第ｎ−１細分化ステップ（ｎ＝２，３，…）で求められた間隔
【０１１６】
【数５３】
2008282268.tif 000055
【０１１７】
ごとの相互相関関数から求められる、誤差を含むピーク位置近傍の、第ｎ−１のピーク位置近傍の第ｎ−１区間よりも狭い第ｎ区間を該間隔
【０１１８】
【数５４】
2008282268.tif 000056
【０１１９】
よりも細かい間隔
【０１２０】
【数５５】
2008282268.tif 000057
【０１２１】
で細分化した各細分点の相互相関関数を求める第ｎ細分化ステップとを経ることにより所望精度のピーク位置を求めるものであってもよい。
【０１２２】
また、本発明のピーク位置探索装置において、上記標本化間隔
【０１２３】
【数５６】
2008282268.tif 000058
【０１２４】
が標本化時間間隔
【０１２５】
【数５７】
2008282268.tif 000059
【０１２６】
であり、前記データ取得部は、該標本化時間間隔
【０１２７】
【数５８】
2008282268.tif 000060
【０１２８】
で標本化された２つの標本化時系列データ
【０１２９】
【数５９】
2008282268.tif 000061
【０１３０】
但し、
【０１３１】
【数６０】
2008282268.tif 000062
【０１３２】
での信号値
【０１３３】
【数６１】
2008282268.tif 000063
【０１３４】
を
【０１３５】
【数６２】
2008282268.tif 000064
【０１３６】
などと略記した、
を取得するものであってもよい。
【０１３７】
さらに、本発明のピーク位置探索装置において、上記データ取得部の前段に、アナログ信号を受け取りそのアナログ信号を標本化時間間隔
【０１３８】
【数６３】
2008282268.tif 000065
【０１３９】
で標本化してデジタルの標本化時系列データを生成するアナログ・デジタル変換部を備えたことが好ましい。
【０１４０】
本発明のピーク位置探索装置は、そのピーク位置探索装置の外部でアナログ信号からデジタルの標本化時系列データに変換されてその変換後の標本化時系列データを取得するものであってもよいが、本発明のピーク位置探索装置がアナログ・デジタル変換部を備え、外部からアナログ信号を受け取り、自分の内部で標本化時系列データに変換するものであってもよい。
【０１４１】
さらに本発明のピーク位置探索装置において、上記アナログ・デジタル変換部のさらに前段に、エイリアジングノイズ防止用のアンチエイリアジングフィルタを備えることが好ましい。
【０１４２】
こうすることにより、自分自身でアンチエイリアジングノイズの混入の防止をコントロールすることができる。
【０１４３】
また、上記目的を達成する本発明のピーク位置探索プログラムは、コンピュータ内で実行されコンピュータを、２つの標本化データの相互相関のピーク位置を求めるピーク位置探索装置として動作させるピーク位置探索プログラムであって、
上記コンピュータを、
所定の標本化間隔
【０１４４】
【数６４】
2008282268.tif 000066
【０１４５】
で標本化された２つの標本化データ
【０１４６】
【数６５】
2008282268.tif 000067
【０１４７】
を取得するデータ取得部と、
データ取得ステップで取得した２つの標本化データそれぞれに
【０１４８】
【数６６】
2008282268.tif 000068
【０１４９】
個のゼロ点を付加した２つの標本化データ
【０１５０】
【数６７】
2008282268.tif 000069
【０１５１】
のそれぞれについて離散フーリエ変換することにより、２つのスペクトル
【０１５２】
【数６８】
2008282268.tif 000070
【０１５３】
を算出するスペクトル算出部と、
スペクトル算出ステップで算出された２つのスペクトルのうちの一方のスペクトル
【０１５４】
【数６９】
2008282268.tif 000071
【０１５５】
の複素共役と他方のスペクトル
【０１５６】
【数７０】
2008282268.tif 000072
【０１５７】
との積の逆離散フーリエ変換をすることにより、相互相関関数
【０１５８】
【数７１】
2008282268.tif 000073
【０１５９】
を求める相互相関関数算出部と、
標本化間隔
【０１６０】
【数７２】
2008282268.tif 000074
【０１６１】
ごとの標本化点
【０１６２】
【数７３】
2008282268.tif 000075
【０１６３】
の途中の点を含む任意点
【０１６４】
【数７４】
2008282268.tif 000076
【０１６５】
の相互相関関数
【０１６６】
【数７５】
2008282268.tif 000077
【０１６７】
を表わす式
【０１６８】
【数７６】
2008282268.tif 000078
【０１６９】
但し、
【０１７０】
【数７７】
2008282268.tif 000079
【０１７１】
は内挿関数であって、
【０１７２】
【数７８】
2008282268.tif 000080
【０１７３】
で表わされる。
に従って、ピーク位置を所望の精度で求めるピーク位置探索部とを備えたピーク位置探索装置として動作させることを特徴とする。
【０１７４】
ここで、本発明のピーク位置探索プログラムにおいても、上記ピーク位置探索部は、相互相関関数算出ステップで求められた標本化間隔
【０１７５】
【数７９】
2008282268.tif 000081
【０１７６】
ごとの相互相関関数から求められる、誤差を含む第１のピーク位置近傍の第１区間を該標本化間隔
【０１７７】
【数８０】
2008282268.tif 000082
【０１７８】
よりも細かい間隔
【０１７９】
【数８１】
2008282268.tif 000083
【０１８０】
で細分化した各細分点の相互相関関数を求める第１細分化ステップと、その第１細分化ステップで求められた間隔
【０１８１】
【数８２】
2008282268.tif 000084
【０１８２】
ごとの相互相関関数から求められる、誤差を含む第２のピーク位置近傍の、第１区間よりも狭い第２区間を間隔
【０１８３】
【数８３】
2008282268.tif 000085
【０１８４】
よりも細かい間隔
【０１８５】
【数８４】
2008282268.tif 000086
【０１８６】
で細分化した各細分点の相互相関関数を求める第２細分化ステップと、以下同様にして所望精度のピーク位置が探索されるまで、第ｎ−１細分化ステップ（ｎ＝２，３，…）で求められた間隔
【０１８７】
【数８５】
2008282268.tif 000087
【０１８８】
ごとの相互相関関数から求められる、誤差を含む第ｎ−１ピーク位置近傍の、第ｎ−１区間よりも狭い第ｎ区間を該間隔
【０１８９】
【数８６】
2008282268.tif 000088
【０１９０】
よりも細かい間隔
【０１９１】
【数８７】
2008282268.tif 000089
【０１９２】
で細分化した各細分点の相互相関関数を求める第ｎ細分化ステップとを経ることにより所望精度のピーク位置を求めるものであること。
【０１９３】
さらに、本発明のピーク位置探索プログラムにおいて、上記標本化間隔
【０１９４】
【数８８】
2008282268.tif 000090
【０１９５】
が標本化時間間隔
【０１９６】
【数８９】
2008282268.tif 000091
【０１９７】
であり、上記データ取得部は、該標本化時間間隔
【０１９８】
【数９０】
2008282268.tif 000092
【０１９９】
で標本化された２つの標本化時系列データ
【０２００】
【数９１】
2008282268.tif 000093
【０２０１】
但し、
【０２０２】
【数９２】
2008282268.tif 000094
【０２０３】
での信号値
【０２０４】
【数９３】
2008282268.tif 000095
【０２０５】
を
【０２０６】
【数９４】
2008282268.tif 000096
【０２０７】
などと略記した、
を取得するものであってもよい。
【発明の効果】
【０２０８】
以上説明したとおり、本発明によれば、相互相関のピーク位置を高精度に探索することができる。
【発明を実施するための最良の形態】
【０２０９】
以下、本発明の実施形態について説明する。ここでは、本発明を移動速度の測定に適用した例を挙げて説明する。
【０２１０】
図１は、本発明のピーク位置探索装置の一実施形態として動作するピーク位置探索用コンピュータの外観斜視図である。本発明の一実施形態としてのピーク位置探索装置は、このピーク位置探索用コンピュータ１００のハードウェアとその内部で実行されるソフトウェアとの組合せにより実現されている。
【０２１１】
このピーク位置探索用コンピュータ１００は、ＣＰＵ、ＲＡＭメモリ、磁気ディスク、通信用ボード等を内蔵した本体１０１、本体からの指示によりその表示画面１０２ａ上に画面表示を行なうディスプレイ１０２、このピーク位置探索用コンピュータ内に、オペレータの指示や文字情報を入力するためのキーボード１０３、表示画面上の任意の位置を指定することによりその位置に表示されているアイコン等に応じた指示を入力するマウス１０４を備えている。
【０２１２】
本体１０１には、ＣＤ−ＲＯＭ１０５（図２参照）が取り出し自在に装填され、装填されたＣＤ−ＲＯＭ１０５をドライブするＣＤ−ＲＯＭドライブも内蔵されている。
【０２１３】
ここでは、ＣＤ−ＲＯＭ１０５に、本発明の一実施形態としてのピーク位置探索プログラムが記憶されており、そのＣＤ−ＲＯＭ１０５が本体１０１内に装填され、ＣＤ−ＲＯＭドライブによりそのＣＤ−ＲＯＭ１０５に記憶されたピーク位置探索プログラムがそのピーク位置探索演算用コンピュータ１００の磁気ディスク内にインストールされる。ピーク位置探索用コンピュータ１００の磁気ディスク内にインストールされたピーク位置探索プログラムが起動されると、このピーク位置探索用コンピュータ１００は、本発明のピーク位置探索装置の一実施形態として動作する。
【０２１４】
本実施形態では、このピーク位置探索用コンピュータ１００には、２つの光電センサ１１，１２が接続されている。これら２つの光電センサ１１，１２は、それぞれ、図１に矢印方向に移動する測定対象物１０に向けて光を照射しその反射光を受光して受光信号をピーク位置探索用コンピュータ１００に入力するものであり、互いに距離ｄだけ離れた位置に置かれている。ピーク位置探索用コンピュータ１００では、２つの光電センサ１１，１２で得られた２つの受光信号の相互相関のピーク位置が求められ、そのピーク位置
【０２１５】
【数９５】
2008282268.tif 000097
【０２１６】
と距離
【０２１７】
【数９６】
2008282268.tif 000098
【０２１８】
とから測定対象物１０の移動速度
【０２１９】
【数９７】
2008282268.tif 000099
【０２２０】
が、
【０２２１】
【数９８】
2008282268.tif 000100
【０２２２】
として計測される。
【０２２３】
２つの光電センサ１１，１２は極めて正確に距離ｄ
だけ離れた位置に置かれており、以下では、ピーク位置探索用コンピュータ１００内での演算によるピーク位置
【０２２４】
【数９９】
2008282268.tif 000101
【０２２５】
を高精度に求める技術について説明する。
【０２２６】
図２は、図１に示すピーク位置探索用コンピュータ１００のハードウェア構成図である。
【０２２７】
このハードウェア構成図には、中央演算処理装置（ＣＰＵ）１１１、ＲＡＭ１１２、磁気ディスクコントローラ１１３、ＣＤ−ＲＯＭドライブ１１５、マウスコントローラ１１６、キーボードコントローラ１１７、ディスプレイコントローラ１１８、およびＡ／Ｄ変換ボード１２０が示されており、それらはバス１１０で相互に接続されている。
【０２２８】
Ａ／Ｄ変換ボード１２０上には、フィルタ１２１およびＡ／Ｄ変換器１２２が２チャンネル分搭載されており、このＡ／Ｄ変換ボード１２０は、図１に示すピーク位置探索用コンピュータ１００の本体１０１内に装填されている。
【０２２９】
このＡ／Ｄ変換ボード１２０上のＡ／Ｄ変換器１２２は、入力されたアナログの受光信号を所定の標本化時間間隔（サンプリング間隔）
【０２３０】
【数１００】
2008282268.tif 000102
【０２３１】
で標本化（サンプリング）してディジタルの標本化時系列データを生成するものであり、またこのＡ／Ｄ変換ボード１２０に搭載されたフィルタ１２１は、Ａ／Ｄ変換器１２２により得られるディジタルの標本化時系列データにエイリアジングノイズ（折り返し雑音）が混入するのを防止するための、低域通過型のアンチエイリアジングフィルタである。
【０２３２】
ＣＤ−ＲＯＭドライブ１１５は、図１を参照して説明したように、ＣＤ−ＲＯＭ１０５が装填され、装填されたＣＤ−ＲＯＭ１０５をアクセスするものである。
【０２３３】
また、図２には、磁気ディスクコントローラ１１３によりアクセスされる磁気ディスク１１４、マウスコントローラ１１６により制御されるマウス１０４、キーボードコントローラ１１７により制御されるキーボード１０３、およびディスプレイコントローラ１１８により制御されるディスプレイ１０２も示されている。
【０２３４】
図３は、Ａ／Ｄ変換器１２２の入力側のアナログ受光信号と出力側の標本化時系列データの模式図である。
【０２３５】
図３（Ａ），（Ｂ）は、それぞれ、光電センサ１１，１２から得られアンチエイリアジングフィルタ通過した後のアナログ信号波形であり、図１に示す距離ｄ（既知）と測定対象面１０の移動速度Ｖ（ここでは未知数）とで決まる分だけ位相がずれた波形となっている。
【０２３６】
Ａ／Ｄ変換器１２２では、これら２つのアナログの受光信号が標本化時間間隔（サンプリング間隔）
【０２３７】
【数１０１】
2008282268.tif 000103
【０２３８】
で標本化（サンプリング）されて、○印に相当する各時刻の信号の大きさを表わす数値が並んだ２つの標本化時系列データが生成される。
【０２３９】
図４は、本発明の一実施形態としてのピーク位置探索プログラムの構成を示す図である。
【０２４０】
ここでは、このピーク位置探索プログラム２００は、ＣＤ−ＲＯＭ１０５に記憶されており、このピーク位置探索プログラム２００は、データ取得部２１０、スペクトル算出部２２０、相互相関関数算出部２３０、およびピーク位置探索部２４０から構成されている。
【０２４１】
このピーク位置探索プログラム２００を構成する各部の作用については後述する。
【０２４２】
図５は、本発明の一実施形態としてのピーク位置探索方法の構成を示すフローチャートである。
【０２４３】
このピーク位置探索方法は、データ取得ステップ（ステップａ）、スペクトル算出ステップ（ステップｂ）、相互相関関数算出ステップ（ステップｃ）、およびピーク位置探索ステップ（ステップｄ）とから構成されている。
【０２４４】
図６は、ピーク位置探索ステップ（ステップｄ）の詳細フローを示す図であり、このピーク位置検索ステップ（ステップｄ）では、細分化（ステップｄ１）が所望精度に達するまで（ステップｄ２）、繰り返される。
【０２４５】
図５、図６に示すピーク位置探索方法は、本実施形態では、図４に示すピーク位置探索プログラムが図１、図２に示すピーク位置探索用コンピュータ１００にインストールされて実行されることにより、そのピーク位置探索プログラムにより実施される。
【０２４６】
このピーク位置探索方法の各ステップの作用については、以下に説明する図７のピーク位置探索装置の説明の際に、図４のピーク位置探索プログラムの各部の作用も合わせて一緒に説明する。
【０２４７】
図７は、図１、図２に示すピーク位置探索用コンピュータに図４に示すピーク位置探索プログラムがインストールされて実行されることにより実現された、本発明のピーク位置探索装置の一実施形態の機能ブロック図である。
【０２４８】
この図５に示すピーク位置探索装置３００には、アンチエイリアジングフィルタ３０１、アナログ・デジタル変換部３０２、データ取得部３１０、スペクトル算出部３２０、相互相関関数算出部３３０、およびピーク位置探索部３４０から構成されている。
【０２４９】
アンチエイリアジングフィルタ３０１およびアナログ・デジタル変換部３０２は、図２に示すＡ／Ｄ変換ボード１２０上のフィルタ１２１およびＡ／Ｄ変換器１２２に相当する。それ以外のデータ取得部３１０、スペクトル算出部３２０、相互相関関数算出部３３０、およびピーク位置探索部３４０は、ピーク位置探索用コンピュータ１００のＣＰＵ１１１（図１、図２参照）で実行されるピーク位置探索プログラムによってピーク位置探索用コンピュータ１００内に実現される機能である。図４に示すピーク位置探索プログラム２００および図５に示すピーク位置探索方法との関係では、図７のピーク位置探索装置３００のデータ取得部３１０は、図４のピーク位置探索プログラム２００のデータ取得部２１０および図５のピーク位置探索方法のデータ取得ステップ（ステップａ）に対応し、図７のピーク位置探索装置３００のスペクトル算出部３２０は、図４のピーク位置探索プログラム２００のスペクトル算出部２２０および図５のピーク位置探索方法のスペクトル算出ステップ（ステップｂ）に対応し、図７のピーク位置探索装置３００の相互相関関数算出部３３０は、図４のピーク位置探索プログラム２００の相互相関関数算出部２３０および図５のピーク位置探索方法の相互相関関数算出ステップ（ステップｃ）に対応し、図７のピーク位置探索装置３００のピーク位置探索部３４０は、図４のピーク位置探索プログラム２００のピーク位置探索部２４０および図５のピーク位置探索方法のピーク位置探索ステップ（ステップｄ）に対応する。ただし、図７のピーク位置探索装置３００のデータ取得部３１０、スペクトル算出部３２０および相互相関関数算出部３４０は、それぞれがハードウェアとソフトウェアとの複合からなるものであるのに対し、図４のピーク位置探索プログラム２００の、それらに対応する要素はアプリケーションソフトウェアのみからなるものである。また図５のピーク位置探索方法は、本実施形態では、図４のピーク位置探索プログラムが、図１、図２のピーク位置探索用コンピュータ１００にインストールされてそのピーク位置探索用コンピュータ１００内で実行されることにより、そのピーク位置探索用コンピュータ１００内で実施される。
【０２５０】
以下、この図７のピーク位置探索装置３００の各部の作用を説明することで、図４のピーク位置探索プログラム２００および図５のピーク位置探索方法の説明を兼ねることとする。
【０２５１】
図７のピーク位置探索装置のデータ取得部３１０は、図２に示すＡ／Ｄ変換ボード１２０を制御して、そのＡ／Ｄ変換ボード１２０におけるフィルタリングおよびＡ／Ｄ変換により得られたデジタルの標本化時系列データを、図４のピーク位置探索プログラム２００内に取り込む役割りを成すものであり、ハードウェア上は、図４のピーク位置探索プログラム２００を構成するデータ取得部２１０を実行するＣＰＵ１１１等がこれに相当する。
【０２５２】
このデータ取得部３１０では、上記のようにして、図１に示す２つの光電センサ１１，１２で生成された２つの受光信号に対応する、それぞれがＮ点からなる２つの標本化時系列データ
【０２５３】
【数１０２】
2008282268.tif 000104
【０２５４】
が得られたものとする。ただし、
【０２５５】
【数１０３】
2008282268.tif 000105
【０２５６】
での信号値
【０２５７】
【数１０４】
2008282268.tif 000106
【０２５８】
を
【０２５９】
【数１０５】
2008282268.tif 000107
【０２６０】
などと略記した。
【０２６１】
図７のスペクトル算出部３２０では、上記の２つの標本化時系列データそれぞれのスペクトルが求められる。
【０２６２】
ここでは、先ず、相互相関関数計算後の回り込みノイズの混入を防ぐために、それぞれ標本化時系列データにＮ個のゼロ点を付け加えた、
【０２６３】
【数１０６】
2008282268.tif 000108
【０２６４】
が生成され、それらにＦＦＴ演算処理が形成されて、２つのスペクトル
【０２６５】
【数１０７】
2008282268.tif 000109
【０２６６】
が求められる。
【０２６７】
図７の相互相関関数算出部３３０では、それら２つのスペクトルのうちの一方のスペクトル
【０２６８】
【数１０８】
2008282268.tif 000110
【０２６９】
の複素共役と、もう一方のスペクトル
【０２７０】
【数１０９】
2008282268.tif 000111
【０２７１】
との積の逆フーリエ変換が行なわれ、相互相関関数
【０２７２】
【数１１０】
2008282268.tif 000112
【０２７３】
が、
【０２７４】
【数１１１】
2008282268.tif 000113
【０２７５】
として、算出される。従って、時系列
【０２７６】
【数１１２】
2008282268.tif 000114
【０２７７】
が時系列
【０２７８】
【数１１３】
2008282268.tif 000115
【０２７９】
を、標本化時間間隔
【０２８０】
【数１１４】
2008282268.tif 000116
【０２８１】
の整数倍
【０２８２】
【数１１５】
2008282268.tif 000117
【０２８３】
だけ時間をずらした時系列
【０２８４】
【数１１６】
2008282268.tif 000118
【０２８５】
である場合には、相互相関関数は、
【０２８６】
【数１１７】
2008282268.tif 000119
【０２８７】
となり、
【０２８８】
【数１１８】
2008282268.tif 000120
【０２８９】
で相互相関関数は最大値を取る。しかし、時系列
【０２９０】
【数１１９】
2008282268.tif 000121
【０２９１】
が時系列
【０２９２】
【数１２０】
2008282268.tif 000122
【０２９３】
を、標本化時間間隔
【０２９４】
【数１２１】
2008282268.tif 000123
【０２９５】
の任意実数倍
【０２９６】
【数１２２】
2008282268.tif 000124
【０２９７】
だけ時間をずらした時系列
【０２９８】
【数１２３】
2008282268.tif 000125
【０２９９】
である場合には、相互相関関数のピーク位置は、
【０３００】
【数１２４】
2008282268.tif 000126
【０３０１】
が整数の値を取る位置には存在せず、遅延時間の推定値としては、最大、標本化時間間隔の半分
【０３０２】
【数１２５】
2008282268.tif 000127
【０３０３】
程度の誤差を含んだものとなる。
【０３０４】
そこで、図７のピーク位置探索装置３００のピーク位置探索部３４０では、以下のようにしてピーク位置の正確な探索が行なわれる。
【０３０５】
上記の相互相関関数は、元の時系列
【０３０６】
【数１２６】
2008282268.tif 000128
【０３０７】
が帯域制限されているため、相関関数に対する計算式（Ｄ）からも判るように、相互相関関数も帯域制限されている。したがって、任意のラグタイム
【０３０８】
【数１２７】
2008282268.tif 000129
【０３０９】
に於ける相互相関関数の値は、
【０３１０】
【数１２８】
2008282268.tif 000130
【０３１１】
【数１２９】
2008282268.tif 000131
【０３１２】
で与えられる。この式に従って、相互相関関数のピーク位置を求めるには、まず、最初の離散データから得られる相互相関関数のピーク付近の位置３点を求め、この区間を細分し、（Ａ）式を用いて細分点での相互相関関数数値を求める。更に、得られた相互相関関数のピーク付近の位置３点を求め、この区間を更に細分し、（Ａ）式を用いて再細分点での相互相関関数を求める。こういったことを所望の時間精度に達するまで反復することにより、ピーク位置を所望の高精度に求めることができる。
【０３１３】
図８は、ピーク位置探索部３４０における繰り返し演算の説明図である。
【０３１４】
図８（Ａ）は、相互相関算出部３３０で求められた相互相関関数のピーク位置近傍を模式的に示す図である。
【０３１５】
○印を付した、間隔
【０３１６】
【数１３０】
2008282268.tif 000132
【０３１７】
ごとの離散的な各点の相互相関関数が求められている。それらの点を結ぶ線は便宜的に付したものである。
【０３１８】
ここでは、図８（Ａ）の相互相関関数から、誤差を含んだ第１のピーク位置ａ１を挟む両側ａ２，ａ３の第１区間が認識され、今度は図８（Ｂ）に示すように、（Ａ）式によりその第１区間内を、間隔
【０３１９】
【数１３１】
2008282268.tif 000133
【０３２０】
よりも小さい間隔
【０３２１】
【数１３２】
2008282268.tif 000134
【０３２２】
で細分した各点について相互相関関数が求められる。尚、図８（Ｂ）の縦軸、横軸のスケールは、図８（Ａ）と比べ拡大されている。
【０３２３】
このときの相互相関関数は、単なる補間演算とは違い、理論上極めて正確に算出される。
【０３２４】
図８（Ｂ）でも、誤差を含んだ第２のピーク位置ｂ１を間に置いた両側ｂ２（＝ａ１），ｂ３の第２区間が認識され、今度は、図８（Ｃ）に示すように、（Ａ）式により、その第２区間内を間隔
【０３２５】
【数１３３】
2008282268.tif 000135
【０３２６】
よりもさらに小さい
【０３２７】
【数１３４】
2008282268.tif 000136
【０３２８】
で細分し、その細分した各点について相互相関関数が求められる。
【０３２９】
尚、図８（Ｃ）は図８（Ｂ）と比べ、縦軸、横軸のスケールがさらに拡大されている。
【０３３０】
さらに、図８（Ｃ）における第３のピーク位置ｃ１を間に置いた両側ｃ１，ｃ３＝ｂ１の第３区間が認識され、その第３区間がさらに細分され、（Ａ）式に従って、その細分された各点の相互相関関数が求められる。
【０３３１】
この細分化の演算が、必要な時間精度が得られるまで繰り返されて、ピーク位置が極めて高精度に求められる。
【０３３２】
ここでは、図８（Ｃ）の段階で必要な時間精度が得られたものとすると、その段階でのピーク位置
【０３３３】
【数１３５】
2008282268.tif 000137
【０３３４】
を用いて、図１に示す測定対象面１０の移動速度
【０３３５】
【数１３６】
2008282268.tif 000138
【０３３６】
が、
【０３３７】
【数１３７】
2008282268.tif 000139
【０３３８】
により極めて高精度に求められる。
【０３３９】
尚、ここでは、本発明を、図１に示す測定対象面１０の移動速度
【０３４０】
【数１３８】
2008282268.tif 000140
【０３４１】
の測定に適用した例を示したが、相互相関のピーク位置を求めて何らかの物理量を計測することは極めて広範な分野で採用されており、本発明は、計測の分野を問わず、相互相関のピーク位置を求めて何らかの物理量を計測する分野に広く適用することができる。
【図面の簡単な説明】
【０３４２】
【図１】本発明のピーク位置探索装置の一実施形態として動作するピーク位置探索用コンピュータの外観斜視図である。
【図２】図１に示すピーク位置探索用コンピュータのハードウェア構成図である。
【図３】Ａ／Ｄ変換器の入力側のアナログ受光信号と出力側の標本化時系列データの模式図である。
【図４】本発明の一実施形態としてのピーク位置探索プログラムの構成を示す図である。
【図５】本発明の一実施形態としてのピーク位置探索方法の構成を示すフローチャートである。
【図６】ピーク位置探索ステップの詳細フローを示す図である。
【図７】本発明のピーク位置探索装置の一実施形態の機能ブロック図である。
【図８】ピーク位置探索部における繰り返し演算の説明図である。
【符号の説明】
【０３４３】
１０    測定対象物
１１，１２    光電センサ
１００    ピーク位置探索用コンピュータ
１０１    本体
１０２    ディスプレイ
１０３    キーボード
１０４    マウス
１０５    ＣＤ−ＲＯＭ
１１０    バス
１１１    ＣＰＵ
１２０    Ａ／Ｄ変換ボード
１２１    フィルタ
１２２    Ａ／Ｄ変換器
２００    ピーク位置探索プログラム
２１０    データ取得部
２２０    スペクトル算出部
２３０    相互相関関数算出部
２４０    ピーク位置探索部
３００    ピーク位置探索装置
３０１    アンチエイリアジングフィルタ
３０２    アナログ・デジタル変換部
３１０    データ取得部
３２０    スペクトル算出部
３３０    相互相関関数算出部
３４０    ピーク位置探索部
【図１】
2008282268.tif 000235
【図２】
2008282268.tif 000236
【図３】
2008282268.tif 000237
【図４】
2008282268.tif 000238
【図５】
2008282268.tif 000239
【図６】
2008282268.tif 000240
【図７】
2008282268.tif 000241
【図８】
2008282268.tif 000242