JPH021227A - 周期測定装置 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 ［産業上の利用分野］ この発明は、生体信号、特に胎児の心拍信号を表わす電
気信号の周期を自己相関方式で測定する周期測定装置に
関する。

［従来の技術］ 自己相関方式とは、心拍信号を適当なサンプリング周期
をもってサンプリングし、サンプリングされたデータに
基いて心拍信号の自己相関関数を計算し、計算された自
己相関関数から心拍信号を測定する方式である。自己相
関関数とは、心拍信号の成る時刻における波数がその時
刻から成る時間だけずれた時刻における波形とどれだけ
類似性を有しているかを示すものである。換言するなら
ば、心拍信号の繰返し波形の類似度を示すものである。

このことを第１図を参照して述べると、［周期がＴなら
ばその周期Ｔで繰返される部分Ｍ、を、時間（周期）Ｔ
だけ時間軸上で移動させると後続する次の部分Ｍ２に最
も精度高く重なるＪというように表現することができる
。

ところで生体信号をｆ　（ｔ）であられすと、自己相関
関数φ（て）は、 で求めることができる。

測定する信号をサンプリングして得られたデータをｆ　
　（ｋ　）　　（ｋ　＝　１　、　２、−ｎ　）とする
と、上式（１）は ｎは１サンプリングサイクルにおける、電気信号であら
れされる心拍信号のサンプリングの回数であり、ｋはサ
ンプリング序数である。

■サンプリングサイクルとは、ｎ回のサンプリングによ
って位相差変数での成る値における１つの自己相関関数
値を計算する過程をいう。

（２）式は展開すると次のようになる。

φ　（τ）　　＝−（ｆ（１）ｆ（１＋で）　＋ｆ　（
２）　ｆ　（２十て）＋・・・÷ｆ　（ｎ）　ｆ　（ｎ
十で））　　　　　・・・（３）すなわち位相差でだけ
ずれている時刻における二つのデータの積の和によって
あられされる。

（１）式において、Ｔは信号の周期を示している。

（１）、（２）、（３）式において、では、としてあら
れされる。

心拍信号についてのある時刻とその時刻からある時間だ
けずれたある時刻までの時間をあられしている。すなわ
ちでは生体信号ｆ　（ｔ）に位相差を与える変数であり
、信号の一周期と考えられる範囲で変化するものである
。

ところで一般的に胎児の心拍信号の、自己相関方式によ
る周期測定の場合について考えてみると、まず所定のサ
ンプリング周期をもって心拍信号をサンプリングするこ
とからはじめる。

そして胎児心拍信号の周期は臨床実験により知られてい
るように極めて広くほぼ３００ｍ５ないし１５００ｍｓ
の範囲にある。したがって従来、測定に際してはでを３
００ｍ５ないし１５００ｍｓの範囲で変化させていた。

実際にはサンプリング方式ではサンプリング周期をＴ、
とするとて／Ｔｓをてとして使用するため、てとして使
用するため、ては３００／ＴＳないし１５００／Ｔｓの
範囲で変化させることになる。この範囲で求めた自己相
関関数は、てが心拍信号の周期Ｔおよびその整数倍の時
間２Ｔ、３Ｔ、・・・の時にピークを有するので、周期
Ｔに相当するピークを検出することによって心拍信号の
周期を求めることができる。

［発明が解決しようとする問題点］ しかしながら広範囲にわたって自己相関関数を計算する
と、信号処理に要する時間が長くなり、実時間処理が強
く望まれる周期測定においては好ましいことではない。

さらに、また広範囲に測定することによって雑音に影響
されるおそれも生じてくる。

この発明は上述のような事情に鑑みなされたものであっ
て、その目的は、雑音に影響されるおそれがなく、かつ
ほぼ実時間で処理することのできる周期測定装置を提供
することである。

また、本発明の他の目的は、各サンプリング時点でピー
クを検出すると共に、不要な計算を省き、その計算結果
を記憶する相関メモリの容量を縮小した周期測定装置を
提案することである。

［問題点を解決するための手段］ 本発明は、胎児の心拍信号の予測最大変化数がほぼ±１
５ＢＰＭ以内であることに注目して、心拍信号の自己相
関関数計算範囲を単位時間当りの最新心拍数上予測最大
心拍変化数に対応する周期の範囲に限定した。

この発明によれば、心拍信号を表わす電気信号を所定の
サンプリング周期でサンプリングするサンプリング手段
と、各サンプリング毎に、前記サンプリング手段によっ
て得られたデータを用いて、前記心拍信号についての自
己相関関数を、前記電気信号に位相差を与える位相差変
数を所定変化範囲に亙り順次変化させて、計算する自己
相関関数計算手段と、位相差変数の前記所定変化範囲に
亘って計算された自己相関関数からピークを検出するピ
ーク検出手段と、該ピーク検出手段の検出したピークの
自己相関関数から心拍信号の周期を計算する周期計算手
段と、前記自己相関関数計算手段が計算する際の位相差
変数の前記所定変化範囲を、前記周期計算手段の計算し
た周期に対応する心拍数を単位時間当りの最新心拍数と
すると、（単位時間当りの最新心拍数±予測最大心拍変
化数）に対応する周期の範囲に設定する計算範囲設定手
段とを具備してなる周期測定装置が提供される。

心拍信号を胎児の心拍信号とした場合には、予測最大心
拍変化数は通常±１５ＰＭＢであるがそれ以上の変化を
生じる場合もあるので若干の余裕をみて±２０ＢＰＭと
するのが好ましい。

しかしながら通常の最大値であるほぼ±１５ＢＰＭに設
定しても、実用上十分な精度を得ることが可能であり、
このようにすれば計算時間をさらに短縮することが可能
であるから特に高い精度を要求されない場合には最大変
化範囲として±１５８ＰＭを選択してもよい。

さらに、前記サンプリング手段の前記サンプリング周期
を前記心拍信号の周期の変化に対応させて変更する手段
を備えることが好ましい。

また、前記サンプリング周期の各変化段階相互間におけ
る変化割合を一定比率することが好ましい。

［実施例］ 以下この発明の実施例を第２図ないし第４図を参照しな
がら説明する。

第２図は、この発明による周期測定装置を、胎児の心拍
信号周期測定に使用した場合について説明するためのも
のであり、横軸は心拍信号周期を示しており、矢印で示
す範囲は、心拍信号周期における自己相関関数計算範囲
、すなわち位相差変数てを変化させる範囲を示している
ものである。

胎児の心拍信号に限らず一般に生体信号の自己相関関数
計算範囲は測定データの精度を実質的に低下せしめない
程度内でできるだけ狭く設定することが実時間処理の観
点から望まれるところである。すなわち測定結果に実質
的に影響を与える信号範囲のみを計算範囲と定め、この
範囲内のみについて演算処理することにより測定データ
の精度を実質的に低下せしめることなく実時間で、処理
することが望まれる。

また計算範囲を不必要に広く設定した場合には、雑音に
影響される可能性も生起される。この観点からも自己相
関関数計算範囲を、測定結果に実質的に影響を与える範
囲で計算するように制御することが望まれるところであ
る。

本発明者は、このような考えを背景として、胎児の単位
時間当りの心拍数の最大変化数が前述のように成る予測
し得る変化数範囲にあることから、位相差変数での変化
範囲をすなわち自己相関関数計算範囲を、（単位時間当
りの最新心拍数上予測最大心拍変化数）に対応する周期
の範囲に制御する制御手段を設けた周期測定装置を案出
したものである。すなわちこの位相差変数てを上記時間
範囲内で変化させてこの自己相関関数を計算することに
より、測定精度を実質的に低下せしめることなく、はぼ
実時間で処理することを可能にしたものである。第２図
を参照して説明すると、横軸は胎児の心拍信号周期を示
しており、矢印で示す範囲は自己相関関数計算範囲を示
しているものであり、３００ｍ５ないし１５００ｍｓの
周期域における計算範囲を矢印範囲で示すように（１分
間当りの最新心拍数±２０ＢＰＭ）に対応する時間の範
囲としている。換言するならば、位相差変数での変化さ
せる範囲を上述の時間範囲に限定している。この例に見
られるように１分間当りの予測最大心拍変化数は、臨床
実験で得られるデータの最大値±１５８ＰＭに対して若
干の余裕をみて多少太き目に例えば±２０ＢＰＭ程度を
選定することが好ましい。これによって計算もれによる
測定データの精度の低下を防止することができる。

第２図に示す実施例では、上述のように自己相関関数の
計数範囲を換言するならば位相差変数τの変化範囲を、
算出される周期に実質的に影響を及ぼす範囲にのみ限定
して計算するように制御したことにより、実質的に意味
のない多量のデータの演算処理が不要となり、相関方式
周期測定方式の実用上の観点から強く望まれている実時
間処理に大きく寄与すると共に、雑音に影響される可能
性も大巾に低減させることができる。

第２図に示すような、相関方式による周期測定方式は、
例えば第３図に示すような構成の周期測定装置によって
達成される。

トランスジューサ２は例えば婦人の腹部Ｗに配置されて
、胎児の心拍信号を検出し、対応する電気信号を発生す
る。トランスジューサ２の出力の心拍信号を表わす電気
信号は、それに接続された前処理回路３によって波形成
形がなされた後、サンプリング回路４において設定され
たサンプリング周期でサンプリングされかつデジタル信
号の形態にアナログ−デジタル変換（Ａ−Ｄ変換）され
る。サンプリングされたデータはサンプリング回路４に
接続されているデータメモリ６に記憶される。データメ
モリ６は複数のシフトレジスタで構成されているもので
あり、常に最新のＮ個例えば２５６個のデータを記憶し
ている。データメモリ６には乗算器８が接続されており
、乗算器８には加算器１０が接続されている。乗算器８
および加算器１０は、データメモリ６に格納されている
データに基いて実質的に（３）式に示す自己相関関数計
算を行い、その結果を、加算器１０に接続されている相
関メモリ１２に格納する。したがって乗算器８および加
算器１０は心拍信号の自己相関関数計算回路と考えるこ
とができる。

相関メモリ１２にはピーク検出器１４が接続されており
、ピーク検出器１４は相関メモリ１２に貯えられた自己
相関関数データからピークを検出する。

ピーク検出器１４には周期計算回路１６が接続されてお
り、周期計算回路１６はピーク検出回路１４からピーク
検出信号を受けて、心拍信号の周期を求める。周期Ｔは
、得られたピークの時間軸上の位置によって決定される
位相差変数γの値に対応した時間すなわちＴ＝γＸｔ　
（ｔ、はサンプリング周期）であられされる。周期計算
回路１６には心拍数計算回路１８が接続されており、心
拍数計算回路１８は周期計算回路１６からの、周期を示
す信号に基づいてその心拍数を計算する。心拍数計算回
路１８は制御回路２０に接続されている。制御回路２０
には、例えば発光ダイオード（ＬＥＤ）を具備して成る
表示器２２が接続されている。表示器２２は制御回路２
０を介して心拍数計算回路１８から出力される信号に基
づいて心拍信号の心拍数を発光表示する。なおこの時、
制御回路２０には心拍数計算回路１８からの信号が雑音
成分を含んでいる場合あるいはプローブはずれが生じた
ような場合、心拍数計算回路１８からの信号が表示器２
２へ入力しないように制御して誤まった心拍数の表示を
防止するような補助手段を備えているとよい。しかし、
これはこの発明と直接関係はないので、詳細な説明は省
略する。制御回路２０は、さらに、位相差γの変化範囲
すなわち自己相関関数の計算範囲を設定する計算範囲設
定回路２４に接続されている。計算範囲設定回路２４は
、さらに乗算器８と加算器１０とに接続されている。制
御回路２ｏにはさらに基準レベル検出器２６が接続され
ており、基準レベル検出器２６はサンプリング回路４に
接続されている。

上述のような構成において計算範囲設定回路２４は、心
拍数計算回路１８からの心拍数をあられす信号を受けて
、（その心拍数±２０ＢＰＭ）に対応する時間を計算す
る。そして、乗算器８に、この時間範囲内における位相
差変数γについての信号を出力する。なおこの場合、制
御回路２０は、心拍数計算回路１８からの信号が雑音成
分を含んでいるような場合、あるいはプローブはずれが
生じたような場合に心拍拍数計算回路１８からの信号が
計算範囲設定回路２４に入力されるのを阻止するように
制御する。乗算器８は、計算範囲設定回路２４から入力
されてくる信号が示す位相差変数γの値だけ離れた二つ
のサンプリングデータをデータメモリ６から読出し、乗
算する。計算範囲設定回路２４はさらに加算器１０にタ
イミング信号を出力し、加算器１０はタイミング信号に
応答して乗算器８で計算されたデータと相関メモリ１２
からの対応する位相差変数γについてのそれ以前の自己
相関数計算値とを読出して計算する。加算結果は再び相
関メモリ１２の所定のアドレスに格納される。このよう
な演算をデータのサンプリング毎に行い、心拍信号の自
己相関関数を相関メモリ１２に格納する。

このようにして相関メモリ１２に格納された自己相関関
数から前述のようにピーク検出器１４かピークを検出し
、周期計算回路１６で、このピークの時間軸上の位置に
より決定される位相差変数γの値に対応する時間であら
れされる心拍信号の周期を計算し、心拍数計算回路１８
でその周期から心拍数を計算する。

制御回路２０はさらに適当な時間間隔で基準レベル検出
器２６に信号を出力する。基準レベル検出器２６は制御
回路２０からの信号を受けて、サンプリングされたデー
タに符号付けする場合の最適な基準レベル（ゼロレベル
）を検出するためのものである。詳述すると、サンプリ
ングされたデータに符号付けする際データの正負のバラ
ンスが正確にとれているほど自己相関関数曲線を周期性
を明確に表すものであり、基準レベル検出器２６はその
ために設けられていてサンプリングの際データの最大値
、最小値、あるいは平均値を検出して基準レベルの最適
値を求めるものである。

第３図に示す実施例において、計算範囲設定回路２４に
より位相差変数γの変化範囲すなわち自己相関関数範囲
を、（１分間当りの最新心拍数±２０　Ｂ　ＰＭ）に対
応する時間の範囲に制御することにより、実質的に意味
のない多量のデータをサンプリングして不必要に計算時
間を増大せしめてしまうこともなく、またデータ精度の
実質的な低下を招くことのない周期測定が得られる。

ところで、実時間処理を望む観点からは、はぼ３００ｍ
５ないし１５００ｍｓなる胎児の心拍信号周期の全域に
わたって−様な一定サンプリング周期をもつサンプリン
グすることは望ましいことではない。それは周期の短い
心拍信号領域においてはサンプリング周期を短く設定し
て密なデータ検出を行うことが高精度な測定を達成する
観点から望ましいが、一方周期の長い心拍信号領域にお
いては時間の変化に対して信号変化はそれほと急激では
ないのでサンプリング周期を長く設定しても測定データ
の精度を実質的に低減せしめることにはならず、むしろ
サンプリング周期を短い周期の心拍信号のサンプリング
周期と同じ周期に設定した場合にはデータが実質的に不
要に多量にサンプリングされ、演算回数が無意味に増大
し、実時間測定の大きな妨げとなる。さらには雑音によ
る影響を受ける可能性もある。

このような観点から、自己相関関数計算範囲の限定に加
えて、第４図に示すように、心拍信号の周期の変化に対
応させて段階的にサンプリング周期を変化させ、実質的
に意味のないデータの演算処理をなくすことはさらに好
ましいことである。

心拍信号の周期の変化に対応させて段階的にサンプリン
グ周期を変化させていく他の根拠は、次のようなことで
ある。すなわち、周期は心拍数に逆比例するので、例え
ば心拍数が低くなれば周期は広がっていく。このため上
２０ＢＰ対応する時間の範囲は広がっていくので、サン
プリング周期も長くとらなければならなくなる。このよ
うに心拍信号周期の変化に対応して上２０ＢＰ で、心拍信号周期の変化に対応させてサンプリング周期
の変化させていくことは好ましいことである。

そのための具体的手段としては、例えば広範囲な心拍信
号領域をいくつかの領域に分割し、各領域の心拍信号の
大きさにしたがってそれぞれ対応した大きさのサンプリ
ング周期に定め、心拍数の高い領域、すなわち、心拍信
号周期の小さな領域には短いサンプリング周期を設定し
、一方心拍数の低い領域、すなわち心拍信号周期の大き
な領域には長いサンプリング周期を設定する。

サンプリング周期の設定の一例を第４図を参照して説明
すると、心拍信号周期域に二つの閾値Ｔ　Ｈ　＋　、　
Ｔ　Ｈ　２を定めて胎児心拍信号周期域を三つの領域、
Ｉ，ＴＩ，ＩＩＩに分割し、各領域毎に対応する異なる
サンプリング周期を設定する。

閾値ＴＨ＋　、ＴＨ２としては、それぞれ、例えば６０
０ｍ５，１０１００Ｏが定められ、この場合には、領域
Ｉ，ＩＩ，ＩＩＩの範囲はそれぞれ３００−６００ｍｓ
，６００−１　０００ｍ５，１　０００−１５００ｍｓ
となる。

これらの領域Ｉ　、　　ＩＩ，　ＩＩＩにおけるサンプ
リング周期をそれぞれＴ，−　Ｉ，Ｔ，−ＩＩ，Ｔ．−
Ｉｎとすると、各サンプリング周期の関係は次のように
する必要がある。

Ｔ，　−　Ｉ　＜Ｔ．　−　ＩＩ　＜Ｔ．　−ＩＩＩサ
ンプリング周期Ｔ．−Ｉ・Ｔ．−ＩＩの設定については
、領域Ｉ　、　ＩＩ，　ＩＩＩの分割の形態によっても
異なるが、領域Ｉ　、　ＩＩ，　ＩＩＩを上述の例のよ
うに３００ｍｓ−６００ｍｓ，６００ｍｓ−１０００ｍ
ｓ１　０００ｍ５　−　１　５００ｍ５と設定した場合
には、サンプリング周？ＣＩＩＴ．−　Ｉ，Ｔ．−ＩＩ
は例えばそれぞれ５ｍｓ，　７．　５ｍｓ．　　１　１
．　２５ｍ５と定めることができる。

なお、領域の変更が生じた場合、以前の領域における測
定で得られているサンプリングデータな新たな領域に設
定されているサンプリング周期に対応する周期のデータ
に補正して用いる場合に、補正演算を容易にするために
、隣接する領域相互間におけるサンプリング周期の変更
割合は一定比率とすることが望ましい。特に、この一定
比率は、例えば３／２．４／３等のように分数比であら
れされる一定比率とすることが好ましい。

なお、サンプリング周期の変更領域の数は任意に設定で
きるが、やたらに多くすることは繁雑になるばかりで好
ましいことではない。測定対象。

精度、計算速度の短縮化等を考慮して、例えば実施例に
示したように３個程度の領域に定めることが適当である
。

心拍信号の全周期域を例えば三つの領域に区分し、心拍
信号の周期の変化に対応させて領域を適当に変更するた
めに、第３図に示す実施例では領域設定回路２８が設け
られている。領域設定回路２８は、制御回路２０．サン
プリング回路４８周期計算回路１６に接続されている。

領域設定回路２８は制御回路２０からの領域変更の指示
信号を受けて領域の変更を行う。制御回路２０は心拍数
計算回路１８からの心拍数を示す信号を受け、その心拍
数に対応する周期を計算し、その周期の属する領域を指
示する信号を出力する。したがって、制御回路２０は計
算して得られた心拍信号の周期がその時設定されている
領域における周期範囲を越えた場合、その周期が属する
周期範囲の新たな領域を指示する信号を領域設定回路２
８に出力する。例えば今周期範囲が３００ｍｓ　−６０
０ｍｓと定められている領域■が設定されており領域Ｉ
にて測定がなされている場合において、心拍数計算回路
１８から得られた信号の心拍数に対応する周期が例えば
５９０ｍ５から６１０ｍ５になったような場合、測定領
域を例えばサンプリング周期５ｍｓである領域■から例
えば６００ｍｓ−１０００ｍｓの周期範囲を定めている
領域ＩＩに変更指示する信号を出力する。領域設定回路
２４はこの変更指示信号を受けて、サンプリング回路４
にサンプリング周期変更信号を出力し、サンプリング回
路４におけるサンプリング周期を、領域ＩＩに予め設定
されているサンプリング周期例えば７．５ｍｓに変更す
る。このように、測定された心拍数に対応する周期が設
定領域において予め定められている周期範囲を越えると
、領域の変更が行われ、サンプリング周期が新たな領域
において予め設定されている周期に変更される。

領域設定回路２８はまた設定された領域において定めら
れているサンプリング周期を示す信号を周期計算回路１
６に出力する。周期計算回路１６は、ピーク検出器１４
からのピーク検出信号を受けて、このピーク時間軸上の
位置によって決定される位相差変数γの値に対応する時
間であられされる心拍信号の周期Ｔ＝ｔＸγ、（て、は
サンプリング周期）を求める。

以上のようにして領域を変更してサンプリング周期の変
更が行なわれると共に、心拍信号の周期が計算される。

［発明の効果］ 以上述べたようにこの発明によれば、心拍信号の周期測
定において、検出された最新の周期に従って、次の自己
相関関数の計算の位相差変数の変化範囲を、算出される
周期に実質的に影響を与える範囲、例えば（単位時間当
りの最新心拍数±２０　Ｂ　ＰＭ）に対応する周期の範
囲になるよう制御したことにより、実質的に意味のない
多量のデータのサンプリングや無駄な計算及び相関メモ
リの記憶容量を省略することが可能になり、雑音に影響
されるおそれがない他、計算時間の短縮により実質的に
ほぼ実時間で処理することのできる周期測定装置が提供
される。

さらに、この発明による自己相関関数の計算範囲の限定
に加えて、心拍信号の周期の変化に対応させてサンプリ
ング周期を変化させることにより、データ精度を実質的
に低下させることなく、計算時間を一層短縮せしめたも
のとし、かつサンプリング周期の変化を段階的に一定比
率をもって変化させることにより古いデータを補正して
そのまま新しいデータとして使用することができそれに
よって連続測定を可能とし、はぼ実時間で処理される周
期測定装置が提供される。

【図面の簡単な説明】

第１図は、自己相関方式による周期測定を説明するため
に用いた心拍信号周期図、 第２図は、この発明による周期測定装置を胎児の心拍信
号周期の測定に適用した場合について説明するための図
、 第３図は、この発明の周期測定装置の一実施例の構成を
ブロックダイヤグラムの形態で概略的に示す図、 第４図は、この発明の周期測定において心拍信号の周期
の変化に対応して段階的にサンプリング周期を変化せし
める方式を説明するための図である。 図中、２・・・トランスジューサ、３・・・前処理回路
、４・・・サンプリング回路、６・・・データメモリ、
８・・・乗算器、１２・・・相関メモリ、１４・・・ピ
ーク検出器、１６・・・周期計算回路、１８・・・心拍
数計算回路、２０・・・制御回路、２２・・・表示器、
２４・・・計算範囲設定回路、２６・・・基準レベル検
出器、２８・・・領域設定回路、Ｔ−Ｉ、Ｔ、−ＩＩ。 ＩＩＩ −ＩＩＩ・・・サンプリング周期である。 特許

Claims (5)

【特許請求の範囲】
1. （１）心拍信号を表わす電気信号を所定のサンプリング
   周期でサンプリングするサンプリング手段と、 各サンプリング毎に、前記サンプリング手段によつて得
   られたデータを用いて心拍信号についての自己相関関数
   を、前記電気信号に位相差を与える位相差変数を所定変
   化範囲に亙り順次変化させて計算する自己相関関数計算
   手段と、 位相差変数の前記所定変化範囲に亘つて計算された自己
   相関関数からピークを検出するピーク検出手段と、 該ピーク検出手段の検出したピークの自己相関関数から
   心拍信号の周期を計算する周期計算手段と、 前記自己相関関数計算手段が計算する際の位相差変数の
   前記所定変化範囲を、前記周期計算手段の計算した周期
   に対応する心拍数を単位時間当りの最新心拍数とすると
   、（単位時間当りの最新心拍数±予測最大心拍変化数）
   に対応する周期の範囲に設定する計算範囲設定手段とを
   具備してなる周期測定装置。
2. （２）前記心拍信号は胎児の心拍信号であり、前記予測
   最大心拍変化数は２０ＢＰＭであることを特徴とする特
   許請求の範囲第１項に記載の周期測定装置。
3. （３）前記心拍信号は胎児の心拍信号であり、前記最大
   予測心拍変化数は１５ＢＰＭであることを特徴とする特
   許請求の範囲第１項に記載の周期測定装置。
4. （４）前記サンプリング手段の前記サンプリング周期を
   前記心拍信号の変化に対応させて変更する手段を備えて
   いることを特徴とする特許請求の範囲第１項に記載の周
   期測定装置。
5. （５）前記サンプリング周期の各変化段階相互間におけ
   る変化割合を一定比率とした特許請求の範囲第４項に記
   載の周期測定装置。