(57)【要約】      （修正有）【課題】ネットワークにおけるパケット転送の遅延を測定するための最適な間隔を自動的に決定する。【解決手段】遅延時間測定処理部５１は、ネットワークにおけるパケット転送の遅延時間を測定して、遅延時間の時系列データを生成し、ＦＦＴ処理部５２は、時系列データに対して、高速フーリエ変換を行い、周波数の時系列データを生成し、平滑化処理部５３は、周波数の時系列データに対して、局所線形回帰による平滑化を行い、平滑した周波数の時系列データを生成し、周波数ν＊算出処理部５４は、平滑化した周波数の時系列データを用いて、周波数ν＊と周波数ν（ν＊＜ν）の区間内のパワーの最大と最小の差が全体のパワーの最大Ｐｍａｘと最小Ｐｍｉｎの差より小さくなる最小の周波数ν＊を求め、測定間隔決定処理部５５は、最小の周波数ν＊の逆数を、ネットワークにおけるパケット転送の遅延を測定するための最適な間隔として決定する。【選択図】図５
2010283662.tif 000002
【特許請求の範囲】
【請求項１】
プログラムされたコンピュータ処理によって、
ネットワークにおけるパケット転送の遅延を測定するための最適な間隔を決定するパケット転送遅延測定周期決定システムであって、
プログラムされたコンピュータ処理を実行する手段として、
上記ネットワークにおけるパケット転送の遅延時間を測定し、遅延時間の時系列データを生成して記憶装置に格納する遅延時間測定処理手段と、
上記記憶装置から上記遅延時間の時系列データを読み出し、当該時系列データに対して、高速フーリエ変換を行い、周波数の時系列データを生成して記憶装置に格納するＦＦＴ処理手段と、
上記記憶装置から上記周波数の時系列データを読み出し、該周波数の時系列データに対して、局所線形回帰による平滑化を行い、平滑した周波数の時系列データを記憶装置に格納する平滑化処理手段と、
上記記憶装置から上記平滑化した周波数の時系列データを読み出し、
周波数ν＊と周波数ν（ν＊＜ν）の区間内のパワーの最大と最小の差が全体のパワーの最大Ｐｍａｘと最小Ｐｍｉｎの差より小さくなる最小の周波数ν＊を求めて記憶装置に格納する周波数ν＊算出処理手段と、
上記記憶装置から上記最小の周波数ν＊を読み出し、該最小の周波数ν＊の逆数を、上記ネットワークにおけるパケット転送の遅延を測定するための最適な間隔として決定する測定間隔決定処理手段と
を有することを特徴とするパケット転送遅延測定周期決定システム。
【請求項２】
請求項１に記載のパケット転送遅延測定周期決定システムであって、
上記周波数ν＊算出処理手段は、
下記ステップの手順を実行することで、上記最小の周波数ν＊を求めることを特徴とするパケット転送遅延測定周期決定システム。
＜手順＞
ステップＳ１：νｔｍｐ＝νｍａｘとする。
ステップＳ２：νｔｍｐ＜＝ν＜＝νｍａｘにおいてＰ（ν）＜Ｐｍｉｎ×（Ｐｍａｘ÷Ｐｍｉｎ）×Ｐｔｒを満たすか否かを、以下の手順で判定する。
ステップＳ２−１：  νｔｍｐ＝νとする。
ステップＳ２−２：  Ｐ（ν）＜Ｐｍｉｎ×（Ｐｍａｘ÷Ｐｍｉｎ）×Ｐｔｒを満たすか否かを判定し、満たしていれば、ステップＳ２−３の処理手順に進み、満たしていなければ、ステップＳ２の処理においては条件を満たさないとしてステップＳ３の処理手順に進む。
ステップＳ２−３：  ν×２＜＝νｍａｘであれば、「ν×２」を新たなνとして、ステップＳ２−２の処理手順に戻り、そうでなければ、ステップＳ２の処理においては条件を満たしたものとしてステップＳ３の処理手順に進む。
ステップＳ３：  ステップＳ２での処理手順において条件を満たしているとの判定であれば、「νｔｍｐ÷２」を新たなνｔｍｐとして、ステップＳ２の処理手順を実行し、満たしていなければ、その時点の「νｔｍｐ×２」を最小の周波数ν＊とする。
【請求項３】
プログラムされたコンピュータ処理により、ネットワークにおけるパケット転送の遅延を測定するための最適な間隔を決定するコンピュータ装置によるパケット転送遅延測定周期決定方法であって、
プログラムされたコンピュータ処理実行手順として、
請求項１もしくは請求項２のいずれかに記載のパケット転送遅延測定周期決定システムにおける各手順を実行する手段を具備したことを特徴とするパケット転送遅延測定周期決定方法。
【請求項４】
コンピュータを、請求項１もしくは請求項２のいずれかに記載のパケット転送遅延測定周期決定システムにおける各手段として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、パケットネットワークの品質管理技術に係り、特に、ＩＰ（Ｉｎｔｅｒｎｅｔ  Ｐｒｏｔｏｃｏｌ）ネットワークにおけるパケットの転送品質の測定処理の効率化と精度の向上を図るのに好適な技術に関するものである。
【背景技術】
【０００２】
ＩＰネットワークにおけるパケット転送品質には、例えば非特許文献１，２に記載のように、パケット損失率（ＩＰＬＲ）、平均パケット転送遅延時間（ＩＰＴＤ）や、パケット転送遅延ゆらぎ（ＩＰＤＶ）などの尺度がある。
【０００３】
これらのパケット転送品質は、図１に示すような構成で、次のようにして測定する。
【０００４】
（１）測定対象のネットワーク（図中「ｔａｒｇｅｔ  ＮＷ」と記載）を通過する個々のパケットに対し、入る地点（ＭＰ１）と出る地点（ＭＰ２）において損失の有無と、通過時刻の差分による遅延時間を記録する。
【０００５】
（２）遅延時間を複数測定し、損失の有無をＩＰＬＲ、平均値をＩＰＴＤ，（９９．９％−最小値）をＩＰＤＶとする。
【０００６】
通常、ＩＰネットワークの品質管理においては、図２に示すようにして、上述の測定および統計処理を一定時間毎に実行し、処理結果と予め定められた閾値とを比較することにより、異常を知らせるアラームを発生させたり、記録、保存して事後の検討において利用したりする。
【０００７】
このようなＩＰネットワークの品質管理における測定間隔は、従来は、管理者の経験により決定されており、ネットワーク品質変動の特性を加味したものではなかった。
【０００８】
例えば、従来の人手によるＩＰネットワークの品質管理における測定間隔では、ＶｏＩＰ（Ｖｏｉｃｅ  ｏｖｅｒ  Ｉｎｔｅｒｎｅｔ  Ｐｒｏｔｏｃｏｌ）の保留時間を元にして３分間隔で測定する、または、機器の性能上の制約により５分間隔で測定する、などとして決定されていた。
【０００９】
しかし、本来、測定間隔は、その目的である、品質の把握のために適切な間隔で実行されるべきである。例えば、測定間隔は、短時間で繰り返すと変動が詳細に把握できる。しかし、過度に細かい場合、大きな変動を把握できない上に、データ保存のため多くのリソースが必要となる、などの問題がある。
【００１０】
広く利用されているＩＰネットワークでは、長期依存性が観測されているが、この条件のもとでは、ネットワークの遅延時間には、長時間で観測すると大きな変動はあっても、短時間の特性はランダムに近似できるという特性が、例えば非特許文献３，４等において報告されている。
【００１１】
図３においては、コンピュータシミュレーションにより生成したパケット転送時間を時系列で示している。これらに、周期性を分析する技術である、ピリオドグラムを適用すると、図４に示すようになる。
【００１２】
図４において、横軸は周波数（Ｆｒｅｑｕｅｎｃｙ）、縦軸はパワー（Ｐｏｗｅｒ）を表している。このグラフでは、水平になっていると各周期成分が均等に混合されており、元の測定値がランダムに変化していることを表している。
【００１３】
周波数の高い短周期の領域（図４のグラフにおける右側）では、グラフが水平になっており、周波数の低い長周期の領域（図４のグラフにおける左側）では、グラフに傾きがある。よって、長い周期の変動はあるが、短時間で見ると、ランダムに変動していることとなる。
【００１４】
ネットワーク管理においては、その境目（図４における丸印の部分）にあたる波長を測定間隔とすれば、測定間隔より長い時間の変動はデータとして蓄積し、測定周期より短い時間の変動は統計処理により丸めることとなる、すなわち、ランダムな変動は丸めつつ、品質の変動を把握することが可能となる。
【００１５】
しかし、従来、図４において、丸印をつけた部分にあたる波長を求める際、グラフを図示し、人間が判断するしかなかった。
【先行技術文献】
【非特許文献】
【００１６】
【非特許文献１】ＩＴＵ−Ｔ  Ｒｅｃｏｍｍｅｎｄａｔｉｏｎ  Ｙ．１５４０，“Ｉｎｔｅｒｎｅｔ  ｐｒｏｔｏｃｏｌ  ｄａｔａ  ｃｏｍｍｕｎｉｃａｔｉｏｎ  ｓｅｒｖｉｃｅ−ＩＰ  ｐａｃｋｅｔ  ｔｒａｎｓｆｅｒ  ａｎｄ  ａｖａｉｌａｂｉｌｉｔｙ  ｐｅｒｆｏｒｍａｎｃｅ  ｐａｒａｍｅｔｅｒｓ”，ＩＴＵ−Ｔ  Ｒｅｃｏｍｍｅｎｄａｔｉｏｎ  Ｙ．１５４０，２００７．
【非特許文献２】ＩＴＵ−Ｔ  Ｒｅｃｏｍｍｅｎｄａｔｉｏｎ  Ｙ．１５４１，“Ｉｎｔｅｒｎｅｔ  ｐｒｏｔｏｃｏｌ  ｃｏｍｍｕｎｉｃａｔｉｏｎ  ｓｅｒｖｉｃｅ−ＩＰ  ｐｅｒｆｏｒｍａｎｃｅ  ａｎｄ  ａｖａｉｌａｂｉｌｉｔｙ  ｏｂｊｅｃｔｉｖｅｓ  ａｎｄ  ａｌｌｏｃａｔｉｏｎｓ”，ＩＴＵ−Ｔ  Ｒｅｃｏｍｍｅｎｄａｔｉｏｎ  Ｙ．１５４１，２００６．
【非特許文献３】織克典，奥山隆文，“パケット転送遅延の周波数分析とネットワーク管理への応用”，信学技報，ＣＱ２００８−５２，ｐｐ．３９−４４，ＩＥＩＣＥ，ＩＥＩＣＥ，２００８．
【非特許文献４】Ｈｉｒｏｋｉ  Ｆｕｒｕｙａ，Ｈａｊｉｍｅ  Ｎａｋａｍｕｒａ，Ｓｈｉｎｉｃｈｉ  Ｎｏｍｏｔｏ，ａｎｄ  Ｔｅｔｓｕｙａ  Ｔａｋｉｎｅ，“Ｌｏｃａｌ  Ｐｏｉｓｓｏｎ  ｐｒｏｐｅｒｔｙ  ｏｆ  ａｇｇｒｅｇａｔｅｄ  ｉｐ  ｔｒａｆｆｉｃ”，ＩＥＩＣＥ  Ｔｒａｎｓ．Ｃｏｍｍｕｎ．，ｖｏｌ．８６，ｎｏ．８，ｐｐ．２３６８−２３７６，２００３．
【発明の概要】
【発明が解決しようとする課題】
【００１７】
解決しようとする問題点は、ＩＰネットワークの品質管理のための測定間隔は、基本的に、短時間で繰り返すと変動が詳細に把握できるが、過度に細かい場合、大きな変動を把握できない上に、データ保存のため多くのリソースが必要となる、などの問題があるが、従来、ＩＰネットワークの品質管理のための測定間隔は、管理者の経験により決定されており、ネットワーク品質変動の特性を加味したものではなく、品質の把握のための適切な測定間隔が決定されていなかった点である。
【００１８】
本発明の目的は、これら従来技術の課題を解決し、ＩＰネットワークにおけるパケット転送遅延の測定処理の効率化と測定精度の向上を図り、ＩＰネットワークの品質管理の精度を向上させることである。
【課題を解決するための手段】
【００１９】
上記目的を達成するため、本発明では、プログラムされたコンピュータ処理により、ネットワークにおけるパケット転送の遅延を測定するための最適な間隔を自動的に決定する。すなわち、まず、ネットワークにおけるパケット転送の遅延時間を測定して、遅延時間の時系列データを生成し、次に、この時系列データに対して、高速フーリエ変換を行い、周波数の時系列データを生成し、さらに、この周波数の時系列データに対して、局所線形回帰による平滑化を行い、平滑した周波数の時系列データを生成し、そして、平滑化した周波数の時系列データを用いて、周波数ν＊と周波数ν（ν＊＜ν）の区間内のパワーの最大と最小の差が全体のパワーの最大Ｐｍａｘと最小Ｐｍｉｎの差より小さくなる最小の周波数ν＊を求め、この最小の周波数ν＊の逆数を、ネットワークにおけるパケット転送の遅延を測定するための最適な間隔として決定する。
【発明の効果】
【００２０】
本発明によれば、ＩＰネットワークにおけるパケット転送の遅延を測定するための最適な間隔（測定間隔）を自動的に決定することができ、ＩＰネットワークにおけるパケット転送遅延の測定処理の効率化と測定精度の向上を図ることができ、ＩＰネットワークの品質管理精度を向上させることが可能となる。
【図面の簡単な説明】
【００２１】
【図１】ＩＰネットワークの品質測定に係るシステム構成例を示す説明図である。
【図２】図１におけるシステムでの測定値に対する統計処理手順例を示す説明図である。
【図３】コンピュータシミュレーションにおいて生成した時系列のパケット転送時間データを示す説明図である。
【図４】図３におけるパケット転送時間データのピリオドグラフを示す説明図である。
【図５】本発明に係るパケット転送遅延測定周期決定システムの構成例を示すブロック図である。
【図６】図５におけるパケット転送遅延測定周期決定システムの本発明に係る処理手順例を示すフローチャートである。
【発明を実施するための形態】
【００２２】
以下、図を用いて本発明を実施するための最良の形態例を説明する。図５においては、は、本発明に係るパケット転送遅延測定周期決定システムの構成例を示しており、本例のパケット転送遅延測定周期決定システム５０は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）や主メモリ等を具備したコンピュータ装置からなり、プログラムされたコンピュータ処理を実行する機能として、遅延時間測定処理部５１、ＦＦＴ処理部５２、平滑化処理部５３、周波数ν＊算出処理部５４、測定間隔決定処理部５５を具備している。尚、メモリ５１ａ〜５５ａはメインメモリやＨＤＤ等の外部記憶装置における記憶領域を利用している。
【００２３】
このような構成により、パケット転送遅延測定周期決定システム５０は、図１におけるＩＰネットワーク（ｔａｒｇｅｔ  ＮＷ）におけるパケット転送の遅延を測定するための最適な間隔を自動的に決定する。
【００２４】
すなわち、遅延時間測定処理部５１は、図１におけるパケットの入る地点ＭＰ１と出る地点ＭＰ２におけるパケットの損失の有無と通過時刻の差分による遅延時間を測定し、その測定結果としての、図３に示すような遅延時間の時系列データを生成してメモリ５１ａに格納する。
【００２５】
ＦＦＴ処理部５２は、メモリ５１ａに格納された遅延時間の時系列データを読み出し、当該時系列データに対して、例えば「森下巌，小畑秀文，「信号処理」，コロナ社，１９８２．」等において記載の高速フーリエ変換（ＦＦＴ：Ｆａｓｔ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍ）を行い、図４に例示する細かい実線相当のデータを求め、当該データをメモリ５２ａに格納する。
【００２６】
ＦＦＴ処理部５２で求められるデータは、図４に例示する実線相当のグラフで示すように、上下に激しく変動しており、水平領域を数値的に算出するのが困難である。
【００２７】
そこで、本例では、平滑化処理部５３において、メモリ５２ａに格納されたデータに対して、例えば、「Ｊ．Ｍ．チェンバース，Ｔ．Ｍ．ヘイスティ，「Ｓと統計モデル−データ科学の新しい波−」，共立出版株式会社，１９９９．」等において記載のＬＯＥＳＳ（局所線形回帰）による平滑化を行い、図４に例示する灰色太線相当のデータを得て、メモリ５３ａに格納する。
【００２８】
このようにして平滑化処理部５３において求めた平滑化後のデータも、厳密な意味で水平になっているわけではないため、周波数ν＊算出処理部５４により、特定の周波数（ν＊）より周波数の大きい区間（ν＊＜ν）を選択した時、その区間内のパワーの最大、最小の差が全体のパワーの最大Ｐｍａｘ、最小Ｐｍｉｎの差と比較して、相対的に小さくなるよう、区間を選択する。
【００２９】
そのため、相対的なパワーの小ささを表すパラメータＰｔｒを決めておく。そして、周波数νにおけるパワーをＰ（ν）と表すと、周波数ν＊算出処理部５４では、ν＊＜νにおいて、常にＰ（ν）＜Ｐｍｉｎ×（Ｐｍａｘ÷Ｐｍｉｎ）×Ｐｔｒを満たす最小の周波数ν＊を求め、メモリ５４ａに格納する。
【００３０】
図４に示す例においては、Ｐｔｒ＝１０１÷３の時、丸印を付けた点の周波数が最小の周波数ν＊に相当する。
【００３１】
そして、測定間隔決定処理部５５は、周波数ν＊算出処理部５４が求めた区間の最小周波数ν＊の逆数を、測定間隔として決定し、メモリ５５ａに格納する。
【００３２】
このようにして、本例のパケット転送遅延測定周期決定システム５０によれば、ＩＰネットワークにおけるパケット転送の遅延を測定するための最適な間隔を自動的に決定することができ、ＩＰネットワークにおけるパケット転送遅延の測定処理の効率化と測定精度の向上を図ることができ、ＩＰネットワークの品質管理精度を向上させることが可能となる。
【００３３】
次に、図６を用いて、このようなパケット転送遅延測定周期決定システム５０による測定間隔を求める処理手順例を説明する。
【００３４】
以下の前提１〜５のもとで、測定間隔を求める。
【００３５】
前提１：  測定間隔を得るため、事前にサンプル測定を行い、時間間隔１÷νｍａｘ毎の遅延データが得られている。
【００３６】
前提２：  遅延データに損失はないものとする。または、損失があった場合、適切な遅延時間に置き換える（前述の非特許文献１では、３秒が推奨されている）
【００３７】
前提３：  遅延時間の時系列は、ピリオドグラム分析、すなわち、フーリエ変換が可能である。
【００３８】
前提４：  遅延時間のピリオドグラムに対して、局所線形回帰分析が可能である。
【００３９】
前提５：  パラメータＰｔｒと，周波数ν＊を求める際の粒度ｉが決められている。
【００４０】
ステップＳ６０１においては、図１におけるパケットの入る地点ＭＰ１と出る地点ＭＰ２におけるパケットの損失の有無と通過時刻の差分による遅延時間を測定し、その測定結果としての、図３に示すような遅延時間の時系列データを生成してメモリ６０に格納する。
【００４１】
ステップＳ６０２においては、ステップＳ６０１で生成した遅延時間の時系列データ（遅延データ）に対し、ＦＦＴを適用し、ピリオドグラムを作成する。
【００４２】
ステップＳ６０３においては、ステップＳ６０２で作成したピリオドグラムに対してＬＯＥＳＳ（局所線形回帰）を適用する。
【００４３】
ステップＳ６０４においては、周波数νと対応するＬＯＥＳＳ適用後のパワーＰ（ν）、および、その最大Ｐｍａｘ、最小Ｐｍｉｎ、に対して、周波数ν＊＜周波数νにおいて、常にＰ（ν）＜Ｐｍｉｎ×（Ｐｍａｘ÷Ｐｍｉｎ）×Ｐｔｒを満たす最小の周波数ν＊を求める。
【００４４】
ステップＳ６０５においては、「１÷ν＊」を最適な測定間隔として出力する。ただし、「ν＊＞νｍａｘ」の場合は、事前に行ったサンプル測定からは、最適な測定間隔は導出できないこととする。
【００４５】
尚、ステップＳ６０４においては、下記ステップの手順で、常にＰ（ν）＜Ｐｍｉｎ×（Ｐｍａｘ÷Ｐｍｉｎ）×Ｐｔｒを満たす最小の周波数ν＊を求める。
【００４６】
ステップＳ１：νｔｍｐ＝νｍａｘとする。
ステップＳ２：νｔｍｐ＜＝ν＜＝νｍａｘにおいてＰ（ν）＜Ｐｍｉｎ×（Ｐｍａｘ÷Ｐｍｉｎ）×Ｐｔｒを満たすか否かを、以下の手順で判定する。
【００４７】
ステップＳ２−１：  νｔｍｐ＝νとする。
ステップＳ２−２：  Ｐ（ν）＜Ｐｍｉｎ×（Ｐｍａｘ÷Ｐｍｉｎ）×Ｐｔｒを満たすか否かを判定する。満たしていれば、ステップＳ２−３の処理手順に進み、満たしていなければ、ステップＳ２の処理においては条件を満たさないとしてステップＳ３の処理手順に進む。
【００４８】
ステップＳ２−３：  ν×２＜＝νｍａｘであれば、「ν×２」を新たなνとして、ステップＳ２−２の処理手順に戻り、そうでなければ、ステップＳ２の処理においては条件を満たしたものとしてステップＳ３の処理手順に進む。
【００４９】
ステップＳ３：  ステップＳ２での処理手順において条件を満たしているとの判定であれば、「νｔｍｐ÷２」を新たなνｔｍｐとして、ステップＳ２の処理手順を実行し、満たしていなければ、その時点の「νｔｍｐ×２」を最小の周波数ν＊とする。
【００５０】
以上、図１〜図６を用いて説明したように、本例では、プログラムされたコンピュータ処理により、ネットワークにおけるパケット転送の遅延を測定するための最適な間隔を自動的に決定する。すなわち、まず、遅延時間測定処理部５１により、ネットワークにおけるパケット転送の遅延時間を測定して、遅延時間の時系列データを生成し、次に、ＦＦＴ処理部５２により、この時系列データに対して、高速フーリエ変換を行い、周波数の時系列データを生成し、さらに、平滑化処理部５３により、この周波数の時系列データに対して、局所線形回帰による平滑化を行い、平滑した周波数の時系列データを生成し、そして、周波数ν＊算出処理部５４により、平滑化した周波数の時系列データを用いて、周波数ν＊と周波数ν（ν＊＜ν）の区間内のパワーの最大と最小の差が全体のパワーの最大Ｐｍａｘと最小Ｐｍｉｎの差より小さくなる最小の周波数ν＊を求め、測定間隔決定処理部５５において、この最小の周波数ν＊の逆数を、ネットワークにおけるパケット転送の遅延を測定するための最適な間隔として決定する。
【００５１】
このように、本例によれば、ＩＰネットワークにおけるパケット転送の遅延を測定するための最適な間隔を自動的に決定することができ、ＩＰネットワークにおけるパケット転送遅延の測定処理の効率化と測定精度の向上を図ることができ、ＩＰネットワークの品質管理精度を向上させることが可能となる。
【００５２】
尚、本発明は、図１〜図６を用いて説明した例に限定されるものではなく、その要旨を逸脱しない範囲において種々変更可能である。例えば、本例では、パケット転送遅延測定周期決定システムを、図１における測定対象のネットワーク（ｔａｒｇｅｔ  ＮＷ）における個々のパケットが入る地点（ＭＰ１）と出る地点（ＭＰ２）以外に箇所に設けた構成で説明しているが、パケット転送遅延測定周期決定システムを地点ＭＰ１，２のいずれかに設けた構成としても良い。
【００５３】
また、本例のパケット転送遅延測定周期決定システムを実装するコンピュータ装置に関しては、ＣＰＵと主メモリの他に、表示装置、入力装置、外部記憶装置を具備し、光ディスク駆動装置等を介してＣＤ−ＲＯＭ等の記憶媒体に記録されたプログラムやデータを外部記憶装置内にインストールした後、この外部記憶装置から主メモリに読み込みＣＰＵで処理することにより、各処理部の機能を実行する構成としても良い。尚、記録媒体としては、ＦＤ（Ｆｌｅｘｉｂｌｅ  Ｄｉｓｋ）等を記録媒体として用いることでも良い。また、プログラムのインストールに関しても、通信装置を介してネットワーク経由でプログラムをダウンロードしてインストールすることでも良い。
【符号の説明】
【００５４】
５１：遅延時間測定処理部、５２：ＦＦＴ処理部、５３：平滑化処理部、５４：周波数ν＊算出処理部、５５：測定間隔決定処理部、５１ａ〜５５ａ，６０：メモリ。
【図１】
2010283662.tif 000003
【図２】
2010283662.tif 000004
【図３】
2010283662.tif 000005
【図４】
2010283662.tif 000006
【図５】
2010283662.tif 000007
【図６】
2010283662.tif 000008