JP3741226B2

JP3741226B2 - 脈拍計測装置

Info

Publication number: JP3741226B2
Application number: JP35115895A
Authority: JP
Inventors: 祐子石川; 耕二小熊; 適江
Original assignee: Tanita Corp
Current assignee: Tanita Corp
Priority date: 1995-12-27
Filing date: 1995-12-27
Publication date: 2006-02-01
Anticipated expiration: 2015-12-27
Also published as: JPH09173308A

Description

【０００１】
【発明の属する技術分野】
静止状態のみでなく、体動のある時の脈拍をも計測する脈拍計測装置に関する。
【０００２】
【従来の技術】
従来より図１に示す如くの回路構成により、脈拍センサからの微弱信号を増幅すると共にバンドパスフィルタにより脈拍の周波数域以外のノイズを除去し、その出力値が一定値以上であれば信号を出力するコンパレータにより１、０のデジタル値に変換し、演算装置は信号が現れたときの時間を測定し、それを信号の発生時間として順次記憶してその信号発生の周期から脈拍を計算していた。
【０００３】
他の方法として、センサから得られた波形をコンパレータにより２値化するのではなく、波形のピーク電圧の時点の時間を得て信号発生の時間とする等、センサから得られた波形をＡ／Ｄ変換してその結果を直接記憶して計算処理していた。
【０００４】
しかし静止状態のみでなく、長時間装着して体動のあるときの脈拍をも計測する必要がある場合、非観血の脈拍センサは、体表組織を通して心拍状態を観測しなければ成らず、そのため体動によりセンサの出力は大きく変動し、脈拍による出力値を越えてしまう場合が多く、従来の技術による計測では正確な脈拍数を得ることは出来なかった。
【０００５】
【発明が解決しようとする課題】
本発明の課題は、従来の技術では処理できない、体動に依る、脈拍とは関係のない信号を脈拍信号と分別し、脈拍信号と見なせる信号を得ることが出来なかった期間の脈拍数を推測・補正し長時間帯の正確な脈拍数を得ることにある。
【０００６】
【課題を解決するための手段】
本発明の一つの観点によれば、脈拍を得るためのセンサを有し、センサの出力を処理し、個々の脈波の発生時間を算出するためのアナログ回路、演算装置および記憶装置を有し、脈波の発生時間を順次記憶し、その周期から脈波周期を算出する計算手段を有する脈拍計測装置において、前記計算手段は、一時的に脈波と見なせるパルスがないと判断された場合、脈波と見なせるパルスを得られなかった期間を脈波周期の代表値で除した値に最も近い整数Nを求め、脈波と見なせるパルスを得られなかった期間を前記整数Ｎで除した推定脈波周期を、該脈波と見なせるパルスを得られなかった期間の脈波周期と見なし、その期間の脈拍数を算出決定することを特徴とする脈拍計測装置が提供される。
本発明の一つの実施の形態によれば、前記計算手段は、一時的に脈波と見なせるパルスがないと判断された場合、脈波と見なせるパルスを得られなかった期間を脈波周期の代表値で除した値に最も近い整数Nを求め、前記推定脈波周期の許容上限値×（Ｎ−１）がパルスを得られなかった期間以上であるか、または、推定脈波周期の許容下限値×（Ｎ＋１）が、パルスを得られなかった期間以下となる場合は、確定的な値を得られないとする。
本発明の別の実施の形態によれば、前記脈波周期の代表値を、脈波が得られなかった期間を挟まない最近数回の脈波周期の平均値とする。
本発明のさらに別の実施の形態によれば、前記推定脈波周期の許容限界値を、脈波周期の代表値から脈波周期の代表値の一定率を増減した値とする。
本発明のさらに別の実施の形態によれば、前記推定脈波周期の許容限界値を、直前の脈波周期に直前の脈波周期と前記脈波周期の代表値との差に関連した値で補正した値とする。
本発明のさらに別の実施の形態によれば、確定的な値を得られないとされた場合、継続的に信号を計測し続け、その複数の信号のそれぞれの時刻を一定時間保持しておき、その信号の内３つ以上の信号が、脈波周期の許容下限値以上許容上限値以下で、周期的になったものを脈波信号とみなして脈波周期を算出し、前記確定的な値を得られないとされた以前の脈波周期との平均値を、確定的な値が得られないとされた期間の脈波周期の代表値として決定する。
【０００７】
【実施の形態】
脈拍は急激な変化は無く、過去の脈波周期から、次の周期の範囲は推測されることを利用し、信号が現れる時間間隔が脈波とみなせるか否かを検知パルスの周期に関連して判断し、センサから得られた波形を脈波に依るものと、体動その他の要因に依る物とを分別し、脈波に依る信号と確認できたものだけをカウントし、信号間隔が許容範囲より大きい場合は脈拍信号が欠落したとして有るべく脈拍数を推測する。
本発明の実施の形態をフローチャート図４で示す。脈波周期の代表値をＴtyp 、許容下限値をＴmin、上限値をTmaxとおき、Ｔxから何番目の信号時間について計算するかを指定する変数ｎを１とし、脈波周期の代表値との差を計算するための変数Ｄを十分大きな値Ｄmaxとしておく。
Ｔx+nと、Ｔxとの差をとった周期時間Ｔwを計算し、周期時間Ｔwが許容下限値Ｔminを越えるまで、ｎ＝１番目から、順に計算する。
こうしてＴwがＴmin以下である信号を除外し、さらにＴwがＴmax以下である信号を求め、脈波の可能性がある候補とする。この候補について、代表値との誤差｜Ｔw−Ｔtyp｜を求め、誤差がＤより小であればＤを｜Ｔw−Ｔtyp｜に更新し、ｎ＝ｎ＋１として、次の信号を計算する。
【０００８】
ＴwがＴmax以下である信号について順次誤差を求め、その値が直前の信号の誤差より小さければ順次計算を進め、誤差が直前の信号の誤差より大きくなったら、直前の信号が最も代表値に近いことが解る。その場合はｎ＝ｎ−１として直前の信号に戻し、Ｔx+n−Ｔxを脈波周期として、表示・記録する。
そしてＸ＝Ｘ＋ｎとして、新たに求まった信号を基準として計算を繰り返す。
ＴwがＴmaxを越えた場合は、前の信号が許容上限、下限値に入っていたかを調べる。これはＤ＝Ｄmax でなければ、誤差の最小値を求めている途中であると解り、直前の信号の誤差が最小値と確定するので、上と同様に処理する。
Ｄ＝Ｄmaxであった場合、ＴxとＴx+n との間に、許容範囲に入った信号がまだないことが解る。つまり、脈波信号がノイズによって欠落していると考えられる。この時の信号を脈波信号と仮定し、ＴxとＴx+nの間が何周期分（Ｎ）にあたるか代表値Ｔtypとの差が最小になるＮの値を求める。
【０００９】
求まった周期Ｔw／Ｎが、脈波周期の許容誤差範囲に入るかを調べ、範囲を超えている場合は、信号が脈波によるものではなかったとみなし、次の信号について同様にしてＮを求める。Ｔw／Ｎが範囲に入った場合は、これを脈波周期として、表示、記録し、Ｘ＝Ｘ＋ｎとして、新たに求まった信号を基準として計算を繰り返す。
【００１０】
この方法は、断続後の信号が確実に脈波と認められる必要がある。その信頼性を上げるために、ＴwがＴmaxを越え、ＴxとＴx+nとの間に、許容範囲に入った信号がない場合、つまり脈波の断続が認められた場合、それ以降の信号を所定時間（数秒）記憶しておき、その信号の内、所定回数（数回）の信号が周期的になっているものを、脈波信号とみなし、上記の計算を行う。
【００１１】
正確な脈拍を記録するためには、上記の方法で得た脈拍値が正しいか否かを判定する必要がある場合がある。
この場合、断続した場合の脈波周期の誤差、または許容上限値、下限値を求めておけば、
パルスを得られなかった期間 ≦ 脈波周期の上限値 × （Ｎ−１）
または
パルスを得られなかった期間 ≧ 脈波周期の下限値 × （Ｎ＋１）
となる場合は、確定的な値を得られないとすることにより、信頼性のある値か否かを判定できる。
【実施例】
【００１２】
図３の場合を、最近５拍の平均値を脈波周期の代表値とし、１秒（Ｔtyp ＝１）、推定脈波周期の許容限界を脈波周期の代表値±２０％（Ｔmax＝１．２Ｔtyp，Ｔmin＝０．８Ｔtpy）、Ｄmax を１００として前記実施の形態の方法で測定した場合を検証する。
第１拍目のＴwは１秒であるからＴw＞Ｔminで有り、Ｔw＞Ｔmax で無いので図４のフローチャートの▲１▼→▲２▼→▲３▼を経て○10に行く。ここで、Ｔw−Ｔtypは０で有るからＤ＜｜Ｔw−Ｔtyp｜では無く、○11でＤ＝０に更新し▲９▼を経て▲１▼に戻る。
【００１３】
第２拍目のＴwは１．５（秒）であるから、Ｔw＞Ｔmaxで有るから▲４▼に行くが、Ｄを０にしているためＤ＝Ｄmaxでは無く○12でｎ＝２をｎ−１＝１に戻し、○13で脈波周期を第１拍目の１秒に決定し、脈拍値を計算し表示、記憶して、○16 でＸを１として最初に戻る。
【００１４】
ここでの第１拍目の周期は１．５−１＝０．５（秒）でＴw＜Ｔminであるから▲９▼でｎを２として▲１▼に戻る。
第２拍目の周期は４−１＝３（秒）であるからＴw＞Ｔmax、Ｄ＝Ｄmax であるから▲１▼、▲２▼、▲３▼、▲４▼、▲５▼、▲６▼を経て▲７▼に行く。
｜Ｔw／Ｎ−Ｔtyp｜＝｜３／１−１｜＝２＜１００＝Ｄmax であるから○17でＤを２に、▲６▼でＮを２に変更し再度▲７▼を実行する。
｜Ｔw／Ｎ−Ｔtyp｜＝｜３／２−１｜＝０．５＜２＝Ｄとなり再度同じ経路を経て、Ｄには０．５を、Ｎには３を適用して▲７▼に戻る。
｜Ｔw／Ｎ−Ｔtyp｜＝｜３／３−１｜＝０＜０．５＝Ｄであり更に○17、▲６▼、▲７▼を繰り返す。この時はＤ＝０、Ｎ＝４である。
｜Ｔw／Ｎ−Ｔtyp｜＝｜３／４−１｜＝０．２５０＞０となり▲８▼に行く。
Ｔw／（Ｎ−１）＝３／（４−１）＝１であるからＴmin＜Ｔw／（Ｎ−１）＜Ｔmaxとなり○18に行きＴw／（Ｎ−１）を脈波周期としてこの間の脈拍値を計算する。
【００１５】
この間のパルスを得られなかった期間は３秒であり計算により得られた脈拍数Ｎは３であるが、この数値が信頼できるか否かを確認する。
脈波周期の上限値×（Ｎ−１）＝１．２×２＝２．４＜３
脈拍周期の下限値×（Ｎ＋１）＝０．８×４＝３．２＞３
でパルスを得られなかった期間のの脈拍数は３と決定される。
即ち２拍とすると、３／２＝１．５秒となり許容範囲の上限値を超え、４拍とすると、３／４＝０．７５秒となり許容範囲のとなり下限値以下となり、３拍の時の１秒だけが脈拍周期の許容範囲に入ることが判る。
【００１６】
１０秒間、脈波と見なせるパルスを得られなかった場合は、先の計算と同様に実行され、Ｔtypに最も近い脈拍周期となるＮ＝１０が得られるが、
脈波周期の上限値×（Ｎ−１）＝１．２× ９＝１０．８＞１０
脈拍周期の下限値×（Ｎ＋１）＝０．８×１１＝８．８＜１０
であり、確定した値が得られないことになる。
１０秒間、脈波と見なせるパルスを得られない場合、
脈波周期９拍分と仮定すると１０／９＝１．１０
脈波周期１０拍分と仮定すると１０／１０＝１．００
脈波周期１１拍分と仮定すると１０／１１＝０．９１
脈波周期１２拍分と仮定すると１０／１２＝０．８３
代表値には１０拍が最も近い物であるが、９拍、１１拍及び１２拍の何れもが許容限界に入り、確定することは出来ない。
【００１７】
この場合信頼性のある脈拍数が得られないとしてこの期間の脈拍数の決定は保留し、それ以降の信号を３秒以上記憶しておきその信号の内３つ以上の信号が周期的になった物を脈波信号と見なし計算してこの間の脈波周期の平均を脈波周期の代表値とし、断続期間中の脈波周期は連続して変化したとして、断続前の脈波周期の代表値とで断続期間中の脈波周期を補間して決定する。
【００１８】
【発明の効果】
一時的に脈波と見なせるパルスがないと判断された場合、脈波と見なせるパルスを得られなかった期間を脈波周期の代表値で除した値に最も近い整数Nを求め、脈波と見なせるパルスを得られなかった期間を前記整数Ｎで除した推定脈波周期を、該脈波と見なせるパルスを得られなかった期間の脈波周期と見なし、その期間の脈拍数を算出するので、一時的に脈波と見なせるパルスがない時でも、脈拍数を算出できる。
前記推定脈波周期の許容上限値×（Ｎ−１）がパルスを得られなかった期間以上であるか、または、推定脈波周期の許容下限値×（Ｎ＋１）が、パルスを得られなかった期間以下となる場合は、確定的な値を得られないとするので、脈拍数の計算値の信頼性をあげることが出来る。
通常の脈波周期は狭い範囲しか取り得ないこと（５秒以上（１２拍／分以下）や０．２秒以下（３００拍／分以上）に成ることはない）、及び脈波は通常、急激には変化しないものであること利用した、直前の脈波に関連して随時脈波周期のあり得る範囲を更新する本発明は、確実に脈波による信号では無い信号を除外する事が出来、ノイズによる脈波信号の欠落を推測し補正することが出来、体動により変化する脈拍を正確にカウントすることが出来る。
【図面の簡単な説明】
【図１】従来の脈拍計の回路構成例
【図２】ノイズがない場合の脈波とコンパレータ出力を示す図
【図３】ノイズで出力波形が乱れた場合の出力例
【図４】本発明の基本形態のフローチャート

Claims

脈拍を得るためのセンサを有し、センサの出力を処理し、個々の脈波の発生時間を算出するためのアナログ回路、演算装置および記憶装置を有し、脈波の発生時間を順次記憶し、その周期から脈波周期を算出する計算手段を有する脈拍計測装置において、前記計算手段は、一時的に脈波と見なせるパルスがないと判断された場合、脈波と見なせるパルスを得られなかった期間を脈波周期の代表値で除した値に最も近い整数 N を求め、脈波と見なせるパルスを得られなかった期間を前記整数Ｎで除した推定脈波周期を、該脈波と見なせるパルスを得られなかった期間の脈波周期と見なし、その期間の脈拍数を算出決定することを特徴とする脈拍計測装置。
前記計算手段は、一時的に脈波と見なせるパルスがないと判断された場合、脈波と見なせるパルスを得られなかった期間を脈波周期の代表値で除した値に最も近い整数 N を求め、前記推定脈波周期の許容上限値×（Ｎ−１）がパルスを得られなかった期間以上であるか、または、推定脈波周期の許容下限値×（Ｎ＋１）が、パルスを得られなかった期間以下となる場合は、確定的な値を得られないとすることを特徴とする請求項１に記載の脈拍計測装置。
前記脈波周期の代表値を、脈波が得られなかった期間を挟まない最近数回の脈波周期の平均値とすることを特徴とする請求項１および請求項２に記載の脈拍計測装置。
前記推定脈波周期の許容限界値を、脈波周期の代表値から脈波周期の代表値の一定率を増減した値とすることを特徴とする請求項２に記載の脈拍計測装置。
前記推定脈波周期の許容限界値を、直前の脈波周期に直前の脈波周期と前記脈波周期の代表値との差に関連した値で補正した値とすることを特徴とする請求項２に記載の脈拍計測装置。
確定的な値を得られないとされた場合、継続的に信号を計測し続け、その複数の信号のそれぞれの時刻を一定時間保持しておき、その信号の内３つ以上の信号が、脈波周期の許容下限値以上許容上限値以下で、周期的になったものを脈波信号とみなして脈波周期を算出し、前記確定的な値を得られないとされた以前の脈波周期との平均値を、確定的な値が得られないとされた期間の脈波周期の代表値として決定する請求項２記載の脈拍計測装置。