JP2014054447A - 脈拍測定装置、ならびに、脈拍測定方法および脈拍測定プログラム - Google Patents

脈拍測定装置、ならびに、脈拍測定方法および脈拍測定プログラム Download PDF

Info

Publication number
JP2014054447A
JP2014054447A JP2012201910A JP2012201910A JP2014054447A JP 2014054447 A JP2014054447 A JP 2014054447A JP 2012201910 A JP2012201910 A JP 2012201910A JP 2012201910 A JP2012201910 A JP 2012201910A JP 2014054447 A JP2014054447 A JP 2014054447A
Authority
JP
Japan
Prior art keywords
pulse
frequency
search
pulse wave
intensity
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2012201910A
Other languages
English (en)
Other versions
JP5987578B2 (ja
JP2014054447A5 (ja
Inventor
Kenji Fujii
健司 藤井
Tatsuya Kobayashi
達矢 小林
Toshihiko Ogura
敏彦 小椋
Yukiya Sawanoi
幸哉 澤野井
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Omron Healthcare Co Ltd
Original Assignee
Omron Healthcare Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Omron Healthcare Co Ltd filed Critical Omron Healthcare Co Ltd
Priority to JP2012201910A priority Critical patent/JP5987578B2/ja
Priority to PCT/JP2013/070750 priority patent/WO2014041913A1/ja
Priority to CN201380045775.4A priority patent/CN104602595B/zh
Priority to DE112013004467.6T priority patent/DE112013004467T5/de
Publication of JP2014054447A publication Critical patent/JP2014054447A/ja
Priority to US14/641,546 priority patent/US20150173627A1/en
Publication of JP2014054447A5 publication Critical patent/JP2014054447A5/ja
Application granted granted Critical
Publication of JP5987578B2 publication Critical patent/JP5987578B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02416Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/1118Determining activity level
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6802Sensor mounted on worn items
    • A61B5/681Wristwatch-type devices
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A61B5/7207Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis

Abstract

【課題】被測定者が安静状態になくても、その被測定者の脈拍数を正しく測定できる脈拍測定装置を提供すること。
【解決手段】脈波を脈波センサ(15)によって検知して脈波信号を取得するデータ取得部(31)と、動きを体動センサ(33)によって検知して運動強度信号を取得する運動強度取得部(31)と、脈波信号を記憶する記憶部(32)と、時間領域の脈波信号を周波数領域に変換し脈波信号の周波数スペクトルを求める周波数変換部(31)と、その周波数軸において強度ピークを探索するための探索周波数範囲を設定する探索範囲設定部(31)と、探索周波数範囲において強度ピークを抽出するピーク抽出部(31)と、抽出された強度ピークの周波数に応じて被測定者の脈拍数を求める脈拍数算出部(31)と、を備え、探索範囲設定部は、運動強度信号が示す運動強度に応じて、探索周波数範囲を変更することを特徴とする脈拍測定装置。
【選択図】図2

Description

この発明は脈拍測定装置に関し、より詳しくは、被測定者の血管の脈動を検知して脈拍数を測定する脈拍測定装置に関する。
また、この発明は脈拍測定方法、および、脈拍測定プログラムに関し、より詳しくは、被測定者の血管の脈動を検知して脈拍数を測定するための脈拍測定方法、および、脈拍測定プログラムに関する。
従来、被測定者の脈拍を測定するための装置としては、心電センサが取り付けられたベルトを被測定者の胸部に巻き回して被測定者の心臓の拍動を心電的に測定することにより、被測定者の脈拍数(心拍数)を測定する装置がある。
また、上述の装置が被測定者の心拍を心電的に検知するのに対し、被測定者の血管の脈動を非心電的に検知することにより脈拍数を測定する装置もある。
後者の装置としては、たとえば、光電センサによって被測定者の皮下の血管の脈動を光電的に検知することにより、被測定者の脈拍数を測定する装置がある(たとえば、特許文献1(特開平10−234684号公報)を参照)。
そのような後者の装置においては、被測定者の皮下の血管の脈動を表す信号(脈波信号)を取得し、該脈波信号の時間変動の周期性にもとづいて脈拍数を測定している。
特開平10−234684号公報
しかしながら、非心電的に、たとえば光電的に、被測定者の皮下の血管の脈動を検知することにより被測定者の脈拍数を測定する手法を採る装置では、被測定者が例えば運動しているとき、被測定者の脈拍数を正しく測定することに困難が伴う。
その理由は、測定時に被測定者が運動を行うと、運動により血管に加速度が生じ、それによって血流に乱れが生じるからである。乱れは、脈波信号に外乱成分として重畳される。そのために、脈波信号から脈動に起因した時間変動の周期を抽出することが困難になる。
また、被測定者が運動を行うことにより、被測定者の身体部位に取り付けられたセンサ手段にも加速度が生じ、センサ手段が、身体部位に対して位置ずれを起こしたり、一時的にではあっても、センサ手段が身体部位から乖離したりするような事象が発生する。これら事象もまた、脈波信号に外乱成分として重畳されてしまう。そのような事象も、脈波信号から脈動に起因した時間変動の周期を抽出することを困難にする一因である。
脈波信号において、血管の脈動に起因した信号強度の変動と、上述のような外乱成分に起因した信号強度の変動とを正しく区別することは極めて困難である。そのため、非心電的に、たとえば光電的に被測定者の皮下の血管の脈動を検知することにより被測定者の脈拍数を測定する手法を採る場合、脈波信号に上述の外乱成分が重畳されることを避けるため、被測定者は、測定の間、安静状態を保つ必要があった。
このことは脈拍測定装置の利便性や、測定条件、測定環境の多様性を制限する要因であった。
そこで、この発明の課題は、被測定者が安静状態になくても、その被測定者の脈拍数を正しく測定可能な脈拍測定装置を提供することにある。
また、この発明の課題は、被測定者が安静状態になくても、その被測定者の脈拍数を正しく測定可能な脈拍測定方法を提供すること、および、そのような脈拍測定方法をコンピュータに実行させることができる脈拍測定プログラムを提供することにある。
上記課題を解決するため、この発明の脈拍測定装置は、
被測定者の脈波を脈波センサによって検知して脈拍を表す脈波信号を取得するデータ取得部と、
前記被測定者の動きを体動センサによって検知して前記被測定者が行っている運動の強度を表す運動強度信号を取得する運動強度取得部と、
前記脈波信号を記憶する記憶部と、
前記記憶部に記憶された時間領域の前記脈波信号を周波数領域に変換し、前記脈波信号の周波数スペクトルを求める周波数変換部と、
前記周波数スペクトルの周波数軸において、強度ピークを探索するための探索周波数範囲を設定する探索範囲設定部と、
前記周波数スペクトルの前記設定された探索周波数範囲において強度ピークを抽出するピーク抽出部と、
抽出された前記強度ピークの周波数に応じて前記被測定者の脈拍数を求める脈拍数算出部と、
を備え、
前記探索範囲設定部は、前記運動強度信号が示す運動強度に応じて、前記探索周波数範囲を変更することを特徴とする。
なお、本明細書で、データ取得部は、脈波センサから脈波信号を直接取得してもよいし、それに代えて、脈波センサから脈波信号をサーバ(記憶部を有する)等に一旦記憶させ、そのサーバ等から取得(間接的取得)してもよい。また、運動強度取得部は、体動センサから運動強度信号を直接取得してもよいし、それに代えて、体動センサから運動強度信号をサーバ(記憶部を有する)等に一旦記憶させ、そのサーバ等から取得(間接的取得)してもよい。
また、「脈拍数」とは、単位時間当たりの脈拍の数(例えば、毎分当たりの脈拍の数であるビート・パー・ミニッツ(BPM))を指す。
この発明の脈拍測定装置では、データ取得部は、被測定者の脈波を脈波センサによって検知して脈拍を表す脈波信号を取得する。運動強度取得部は、被測定者の動きを体動センサによって検知して被測定者が行っている運動の強度を表す運動強度信号を取得する。記憶部は、脈波信号を記憶する。周波数変換部は、記憶部に記憶された時間領域の脈波信号を周波数領域に変換して、脈波信号の周波数スペクトルを求める。探索範囲設定部は、上記周波数スペクトルの周波数軸において、強度ピークを探索するための探索周波数範囲を設定する。ピーク抽出部は、周波数スペクトルの設定された上記探索周波数範囲において強度ピークを抽出する。脈拍数算出部は、抽出された強度ピークの周波数に応じて被測定者の脈拍数を求める。そして、探索範囲設定部は、上記の運動強度信号が示す運動強度に応じて、探索周波数範囲を変更する。
ここで、探索範囲設定部が、周波数スペクトルの周波数軸において、強度ピークを探索するための探索周波数範囲を設定することは、被測定者がする運動に由来する周波数成分、高調波成分を、ピーク抽出部による強度ピーク抽出の周波数範囲から除外する、という意味がある。また、被測定者がする運動により、自身の脈拍にも影響が及ぶ。たとえば、被測定者が強い運動を行えば、その脈拍数は増加する傾向を示す。また、被測定者の運動の強度が低下すれば、その脈拍数は減少する傾向を示す。よって、運動強度信号が示す運動強度に応じて脈拍数の変動の傾向を予測し探索周波数範囲を変更することによって、脈波信号のうち被測定者の血管の脈動に由来する周波数成分(の基本周波数成分)が探索周波数範囲に含まれることが確保される。したがって、被測定者が安静状態になくても、その被測定者の脈拍数を正しく算出することができるようになる。
一実施形態の脈拍測定装置では、前記周波数変換部、前記運動強度取得部、前記探索範囲設定部、前記ピーク抽出部、および、前記脈拍数算出部は、予め定められた周期で処理を繰り返し、第1の周期において前記脈拍数算出部が前記被測定者の脈拍数として第1の値を算出した場合には、前記探索範囲設定部は、前記第1の値に対し予め定められた比率範囲に含まれる値を、前記第1の周期に続く第2の周期のための前記探索周波数範囲に設定することを特徴とする。
この形態の脈拍測定装置では、第1の周期において被測定者の脈拍数として第1の値を算出した場合、第1の周期に続く第2の周期のための探索周波数範囲として、第1の値に対し予め定められた比率範囲に含まれる値を設定する。よって、第2の周期においても、探索周波数範囲に被測定者の血管の脈動に由来する周波数成分(の基本周波数成分)が探索周波数範囲に含まれることがより確実になる。したがって、被測定者が安静状態になくても、その被測定者の脈拍数を正しく算出することができるようになる。
一実施形態の脈拍測定装置では、前記運動強度取得部が第3の周期に続く第4の周期において取得した運動強度が前記第3の周期において取得した運動強度よりも大きい場合には、前記探索範囲設定部は、前記第3の周期のための前記探索周波数範囲よりも高周波数側にシフトされた周波数範囲を前記第4の周期のための前記探索周波数範囲に設定することを特徴とする。
この形態の脈拍測定装置では、第4の周期において取得した運動強度が第3の周期において取得した運動強度よりも大きい場合、第3の周期のための探索周波数範囲よりも高周波数側にシフトされた周波数範囲が第4の周期のための探索周波数範囲に設定される。そうすることで、第4の周期においても、探索周波数範囲に被測定者の血管の脈動に由来する周波数成分(の基本周波数成分)が探索周波数範囲に含まれることがより確実になる。したがって、被測定者の運動強度が変化した場合にも、その脈拍数を正しく算出することができるようになる。
一実施形態の脈拍測定装置では、前記探索範囲設定部は、前記第3の周期のための前記探索周波数範囲の周波数幅と同じ幅を有するように、前記第4の周期のための前記探索周波数範囲を設定することを特徴とする。
なお、本明細書で「探索周波数範囲の周波数幅」とは、探索周波数範囲の上限にあたる周波数と、下限にあたる周波数との差の絶対値を指す。ここでの周波数の単位は、BPM等でよい。
この形態の脈拍測定装置では、装置がする処理の負荷が軽減される。
この発明の脈拍数の測定方法は、
脈拍測定装置がする被測定者の脈拍数の測定方法であって、
脈波センサによって、前記被測定者の脈拍を表す脈波信号を取得するデータ取得ステップと、
前記脈波信号を記憶部に記憶する記憶ステップと、
前記記憶部に記憶された時間領域の前記脈波信号を周波数領域に変換し、前記脈波信号の周波数スペクトルを求める周波数変換ステップと、
体動センサによって、前記被測定者が行っている運動の強度を表す運動強度信号を取得する運動強度取得ステップと、
前記周波数スペクトルの周波数軸において、強度ピークを探索するための探索周波数範囲を設定する探索範囲設定ステップと、
前記周波数スペクトルの前記設定された探索周波数範囲において強度ピークを抽出するピーク抽出ステップと、
抽出された前記強度ピークの周波数に応じて前記被測定者の脈拍数を求める脈拍数算出ステップと、を有し、
前記探索範囲設定ステップは、前記運動強度信号が示す運動強度に応じて、前記探索周波数範囲を変更するステップを含むことを特徴とする。
この発明の脈拍数の測定方法によれば、運動強度信号が示す運動強度に応じて探索周波数範囲を変更することによって、脈波信号のうち被測定者の血管の脈動に由来する周波数成分(基本周波数成分)が探索周波数範囲に含まれることが確保される。したがって、被測定者が安静状態になくても、その被測定者の脈拍数を正しく算出することができるようになる。
この発明の脈拍数測定用コンピュータ・プログラムは、上述の脈拍数の測定方法をコンピュータに実行させるためのプログラムである。
この発明の脈拍数測定用コンピュータ・プログラムによれば、コンピュータに上述の脈拍測定方法を実行させることができる。
以上より明らかなように、この発明の脈拍測定装置および脈拍数の測定方法によれば、運動強度信号が示す運動強度に応じて探索周波数範囲を変更することによって、脈波信号のうち被測定者の血管の脈動に由来する周波数成分(基本周波数成分)が探索周波数範囲に含まれることが確保される。したがって、被測定者が安静状態になくても、その被測定者の脈拍数を正しく算出することができるようになる。
また、この発明の脈拍数測定用コンピュータ・プログラムによれば、コンピュータに上述の脈拍数の測定方法を実行させることができる。
この発明の一実施形態の脈拍測定装置の外観の模式的な斜視図である。 この発明の一実施形態の脈拍測定装置の模式的な断面図である。 前記脈拍測定装置の機能的な構成を示すブロック図である。 前記脈拍測定装置の、脈波信号を測定するための脈波センサ部の回路構成を例示する図である。 前記脈拍測定装置の動作フローを示す図である。 脈波信号(時間領域)の一例を示す図である。 脈波信号(時間領域)のAC成分の一例を示す図である。 脈波信号AC成分(周波数領域)の一例を示す図である。 (a)被測定者の運動強度の時間変化の例を示す図である。(b)脈拍数算出タイミングと、各タイミングにおいて脈拍数算出に用いられる脈波信号AC成分の時間範囲との関係を示す図である。 運動強度に応じて変更される探索周波数範囲の例を示す図である。
以下、この発明の実施の形態を、図面を参照しながら詳細に説明する。
図1Aおよび図1Bは、一実施形態の脈拍測定装置の構成を模式的に示している。図1Aは、一実施形態の脈拍測定装置の外観の模式的な斜視図であり、図1Bは、同脈拍測定装置の模式的な断面図である。なお、説明の都合上、図示しない被測定部位の側を本体10の「下面側」、被測定部位の反対側を本体10の「上面側」としている。
この脈拍測定装置1は、本体10と、バンド20と、を含む。図1Aに示すように、脈拍測定装置1は、腕時計のように、バンド20を被測定者の測定部位3(たとえば、手首)に巻き回すことにより、本体を被測定者の手首に固定することができる。
脈拍測定装置1の本体10は、被測定者の被測定部位(図示しない)に対して密着して配置されて被測定部位との接触面を形成する下面13と、当該下面13の反対側に位置する上面11と、を有する。本体10は、下面13に沿った面方向に関して、そのサイズが小さく構成されたくびれ形状wを有する(図1B)。
脈拍測定装置1の本体10は、下面13の側に配置されて被測定者の脈拍を測定する脈波センサとしての測定部15と、上面11の側に配置されて測定部15によって測定された脈拍に関する情報を表示する表示部14と、を備える。下面13の側に配置された測定部15は、測定光(例えば赤外光又は近赤外光)を発光する発光ダイオードのような発光素子16と、フォトダイオード又はフォトトランジスタのような受光素子17と、を備える光学式センサである。発光素子16は、被測定部位に向けて或る発光強度で光を照射する発光部として働く。また、受光素子17は、被測定部位からの反射光又は透過光を受光する受光部として働く。
本体10が被測定部位3に密着して配置された状態で、発光素子16から発せられた測定光(例えば赤外光又は近赤外光)を被測定部位の皮下にある血管(たとえば、動脈)に照射すると、動脈を流れる赤血球によって照射光が反射され、この反射光が受光素子17で受光される。受光素子17で受光される反射光の光量は、動脈の脈動に応じて変化する。したがって、当該測定部15により、脈波情報を検出して脈拍数を計測することができる。なお、図1では、測定部15が下面13に接するように配置されているが、測定部15が本体10の内部に配置されるとともに、本体10の内部に配置された測定部15と本体10の下面13と連通する空間部を備える構成であってもよい。また、図1Aおよび図1Bに示した脈拍測定装置1は、測定部15が発光素子16と発光素子16の近傍に配置された受光素子17とから構成されて、被測定部位3からの反射光を検出するタイプのものを例示しているが、測定部15が発光素子16と発光素子16に対して対向配置された受光素子17とから構成されて、被測定部位3を透過した透過光を検出するタイプとすることもできる。
この脈拍測定装置1では、脈波センサとして光電式センサからなる測定部15を備えているので、簡単な構成で、脈拍を含む脈波情報を精度よく検知することができる。
表示部14が、本体10の上面11の側すなわち頭部13に配置されている。表示部14は、表示画面(例えば、LCD(Liquid Crystal Display)またはEL(Electroluminescence)ディスプレイなど)を含む。表示部14は、被測定者の脈拍に関する情報(例えば、脈拍数)等を表示画面に表示する。当該表示画面の制御は、表示制御部として機能する制御部31(CPU)(後述)によって行われる。
本体10を被測定者の被測定部位3に取り付けるためのバンド20は、本体10を密着保持するための本体保持部21と、被測定部位を取り巻くための取り巻き部25と、を有する。
本体保持部21には、本体10のくびれ形状wの外形サイズと略一致するように開口部が形成され、これにより、くびれ形状wの部分にて、本体10とバンド20とが係合している。
本体保持部21の一端部には、略矩形形状に屈折されたバックル部材22が取り付けられている。バックル部材22の穴23を通して、取り巻き部25の端部24が被測定部位3から外向きに挿通され、折り返されている。
取り巻き部25のうち端部24以外の部分には、外側面(被測定部位3に接する内側面と反対の面)に長手方向に延びる長めの雌側面ファスナーが設けられており、雌側面ファスナーは、端部24に取り付けられた雄側面ファスナー26と着脱自在に係合している。
このようにして、バンド20によって、本体10が被測定部位3に対して密着して保持されている。
図2は、脈拍測定装置1の機能的なブロック構成を示している。この脈拍測定装置1の本体10は、制御部(CPU)31と、記憶部32と、表示部14と、操作部34と、脈波センサ部15と、体動センサ部33と、を含む。脈拍測定装置1は、さらに、図示しない通信部を含んでもよい。その場合、脈拍測定装置1は、図示しない外部の装置との間で、データ通信を行うことができる。
制御部31は、CPU(Central Processing Unit;中央演算処理装置)およびその補助回路を含み、脈拍測定装置1を構成する各部を制御し、記憶部32に記憶されたプログラムおよびデータに従って各種の処理を実行する。すなわち、制御部(CPU)31は、操作部34、脈波センサ部15、体動センサ部33、および、図示しない通信部から入力されたデータを処理し、処理したデータを、記憶部32に記憶させたり、表示部14で表示させたり、通信部から出力させたりする。
記憶部32は、制御部(CPU)31でプログラムを実行するために必要な作業領域として用いられるRAM(Random Access Memory)と、制御部(CPU)31で実行するための基本的なプログラムを記憶するためのROM(Read Only Memory)と、を含む。また、記憶部32の記憶領域を補助するための補助記憶装置の記憶媒体として、半導体メモリ(メモリカード、SSD(Solid State Drive))などを用いることができる。この記憶部32は、被測定者毎に、脈波センサ部15によって検知された被測定者の脈拍を表す脈波信号(特にそのAC成分)を時系列で格納することができる。
操作部34は、例えば、脈拍測定装置1の電源をON又はOFFするために操作される電源スイッチと、被測定者毎の測定結果を記憶部32に保存するためにいずれの被測定者であるか、あるいは、どのような測定を行うか、を選択するために操作される操作スイッチと、を備える。なお、操作部34は、本体10の上面11の側あるいは側面12に設置することができる。
このように、脈拍測定装置1は、単独の装置として構成することが可能である。だが、図示しない通信部を備えることによって、ネットワーク上でも使用可能である。
その通信部は、有線又は無線のネットワークを介して、制御部(CPU)31によって生成されたデータや記憶部32に格納されていたデータをサーバへ送信したり、サーバの制御部(図示しない)によって生成されたデータやサーバの記憶部(図示しない)に格納されていたデータを受信したりするために用いられる。ここで、サーバとあるのは、通常のサーバに加えて、例えば、パーソナルコンピュータのような据え置き型端末、あるいは、携帯電話やスマートフォンやPDA(パーソナル・デジタル・アシスタンツ)やタブレット(tablet)、および、テレビ等のAV機器のリモートコントローラのような携帯型端末、ならびに、テレビ等のAV機器に内蔵されたコンピュータを含む広い概念を意味している。
なお、脈拍測定装置1の各部へ、図示しない電源から、操作部34の電源スイッチに対するユーザの操作に応じて電力が供給される。
図3は、この脈拍測定装置1の脈波センサ部15の回路構成を例示している。この脈波センサ部15は、CPU31の制御の下で動作することにより脈波センサ部15の動作を制御する脈波センサコントローラ41を備える。
脈波センサコントローラ41は、パルス駆動回路42を制御して発光素子16をパルス駆動させる。つまり、パルス駆動回路42は、脈波センサコントローラ41から供給された駆動パルスに応じてnpn形のトランジスタがスイッチングされることにより、発光素子16の発光状態(周波数とデューティ)を制御する。
また、脈波センサコントローラ41は、発光強度制御回路43を制御して発光素子16の発光強度(すなわち駆動電流)を制御する。つまり、発光強度制御回路43は、CPU31に制御された脈波センサコントローラ41からの発光強度制御信号に応じて可変抵抗の抵抗値が変更されることにより該抵抗値によって規定される駆動電流で発光素子16を駆動して、発光素子16の発光強度を制御する。すなわち、発光素子16を流れる駆動電流が大きくなるほど、発光素子16の発光強度(すなわち発光光量)が大きくなる。
受光素子17は、受光した光の強さに応じた光電出力を出力する。脈波センサコントローラ41は上述のようにして発光素子16を制御するとともに、受光感度調整回路44を制御して受光素子17の受光感度(すなわち光電出力のゲイン)を制御する。受光感度調整回路44は、CPU31に制御された脈波センサコントローラ41からの光電出力制御信号に応じて、その可変抵抗の抵抗値を増減させることで、受光素子17からの光電出力(図5Aにおける脈波DC成分PDC)の大きさを調整する。
なお、ここでは、受光素子17からの光電出力を、脈波DC成分PDCと呼んでいる。実際には、受光素子17から出力される光電出力は一定レベル(DC成分)にAC成分が重畳された脈流であるが、その脈動の大きさは光電出力の大きさに較べて極めて小さいため、ここでは、受光素子17からの光電出力を、脈波DC成分PDCとも呼ぶこととする。
受光素子17からの光電出力(図5Aにおける脈波DC成分PDC)は、二手に分岐されて、一方がバンドパスフィルタ(BPF)45に入力され、他方がA/D変換回路(DC成分用ADC)47Dに入力される。
BPF45は、受光素子17からの光電出力PDCからAC成分を取り出す作用を有し、増幅器46は、BPF45からの出力を増幅する作用を有する。BPF45の通過帯域は、ヒトの一般的な脈拍数の範囲(30BPM〜300BPM)に対応した周波数帯域(0.5Hz〜5Hz)を含めばよい。増幅器46からは、光電出力PDCのAC成分(図5Aにおける脈波AC成分PS(t))が出力され、該出力は、A/D変換回路(AC成分用ADC)47Aに入力される。
受光素子17から出力された光電信号PDCは、A/D変換器47Dを経て、アナログ信号からデジタル信号に変換され、ADC47Aの出力からは脈波AC成分PS(t)のデジタル信号がCPU31に入力される。脈波AC成分PS(t)のデジタル信号は、後述するように、被測定者の脈拍数の算出に用いられる。ADC47Dの出力からは、光電信号(脈波DC成分PDC)がCPU31に入力されて、発光強度を制御するためのパラメータ等の演算処理に用いられる。
なお、本例では、ADC47A(AC成分用ADC)およびADC変換回路47D(DC成分用ADC)から出力されるデジタル信号がCPU31へ入力されるが、ADC47A,47Dは、CPU31に内蔵された態様であってもよい。
体動センサ部33は、加速度センサ48を備える。加速度センサ48は、被測定部位に作用する加速度の大きさを測定し、測定結果を増幅器49へ出力する。増幅器49の出力は、A/D変換回路(ADC)50へ入力され、ADC50からは、加速度の情報を含んだデジタル信号がCPU31へ入力される。ここでは、加速度センサ48に作用する加速度の大きさは、被測定者が行う運動の強度とよく対応していると考え、加速度センサ48の出力を、被測定者が行っている運動の強度を表す運動強度信号として利用する。
この脈拍測定装置1は、全体として、図4に示す脈拍測定方法のフローに従って動作する。
概略的には、脈拍測定装置1は、まず、測定開始時に、安静状態にある被測定者の脈拍数(安静時脈拍数)を算出する。そして、脈拍測定装置1は、次回の測定周期においては、安静時脈拍数にもとづいて、周波数領域で表された脈波信号(より具体的には、脈波AC成分)におけるスペクトル強度のピークを探索するべき周波数範囲(探索周波数範囲)を決定し、探索周波数範囲に存在するスペクトル強度のピークを抽出し、抽出された強度ピークの周波数にもとづいて被測定者の脈拍数を算出する。それ以降の測定周期では、脈拍測定装置1は、体動センサ部から出力される被測定者が行っている運動の強度を表す運動強度信号に応じて探索周波数範囲を前回の測定で用いた探索周波数範囲からずらして当該範囲内でスペクトル強度のピークを抽出することにより、前回の測定周期において算出された脈拍数からの脈拍数の変化を追跡するようにして今回の測定周期における脈拍数を算出する。
i) まず、ステップS1に示すように、安静状態での脈拍数を測定するため、CPU31は、体動センサ部33から出力される運動強度信号に基づいて被測定者が安静状態にあるか否か、を判断する。CPU31が、被測定者が安静状態にあると判断した場合(ステップS1における「YES」)、処理は、ステップS2に進む。そうでない場合、CPU31は、ステップS1を予め設定された周期で繰り返す。なお、ステップS1においては、CPU31は、脈波センサ部から取得した脈波信号(脈波AC成分PS(t))の周波数スペクトルを求め、スペクトル強度分布の形状から被測定者が安静状態にあるか否かを判断してもよい。
ii) 次に、ステップS2に示すように、CPU31は、脈波センサ部15から被測定者の脈拍を表す安静時の脈波信号(脈波AC成分PS(t))を取得するデータ取得部として動作する。より具体的には、データ取得部として動作するCPU31は、光電信号PDCに含まれるAC成分PS(t)を取得する。(図5Aおよび図5B参照。)
図5Aは、受光素子17から出力される光電信号(脈波DC成分PDC)の一例を示す図である。図5Aでは、横軸が時間(秒)を表し、縦軸が脈波DC成分PDCの強度(単位は任意)を表す。光電信号(脈波DC成分PDC)は、上述のように、微小なAC成分を含んだ脈流である。すなわち、脈波DC成分PDCは、組織や滞留している血液などに吸収及び散乱された光からに起因する、周期的に変動しない一定レベルの成分(直流成分)に、生体の脈動(すなわち血液の脈波)を反映して周期的に変動する成分(交流成分)PS(t)が重畳された脈流として出力される。なお、通常は、一定レベルの成分(直流成分)の大きさに比して、周期的に変動する脈波AC成分PS(t)の大きさ(振幅)は2桁程度小さい。このため、光電信号(脈波DC成分PDC)から脈波AC成分PS(t)を取り出して、データとして取り扱えるように増幅するのが望ましい。本例では、増幅器46はオペアンプを含み、CPU31の制御の下、入力抵抗と帰還抵抗との抵抗比を調整することによって、脈波AC成分の増幅ゲインを制御する。増幅器46から出力された脈波AC成分PS(t)は、ADC47Aを経て、デジタル信号の脈波AC成分PS(t)となってCPU31に入力される。
図5Bは、CPU31に入力される脈波AC成分PS(t)の波形を例示している。なお、図5Bでは、横軸が時間(秒)を表し、縦軸が脈波AC成分PS(t)の強度(単位は任意)を表す。脈波AC成分PS(t)は、生体の脈動(すなわち血液の脈波)に応じて周期的に変化している。つまり、脈波AC成分PS(t)は、被測定者の脈拍を表す脈波信号である。この脈波AC成分PS(t)は、図2中に示した記憶部32に、時系列で記憶される。
iii) 次に、図4のステップS3に示すように、CPU31は、記憶部32に記憶された時間領域の安静時の脈波信号(脈波AC成分PS(t))を周波数領域に変換し、脈波信号(脈波AC成分PS(t))の周波数スペクトル(PS(f))を求める周波数変換部として動作する。より具体的には、周波数変換部として動作するCPU31は、記憶部32に記憶された時間領域の安静時の脈波信号(脈波AC成分PS(t))を周波数領域に変換して、安静時の脈波AC成分の周波数スペクトルPS(f)を求める。この例では、CPU31は、周波数変換として動作することにより、安静時の脈波信号(脈波AC成分PS(t))について高速フーリエ変換(FFT)を行う。CPU31は、図7(b)に例示するように、時系列で記憶部32に記憶された安静時の脈波AC成分PS(t)の予め定められた或る長さの期間Td(たとえば、16秒間、8秒間、4秒間等)に含まれる安静時のAC成分PS(t)の周波数スペクトルPS(f)を求める。
図6は、周波数領域に変換された安静時の脈波AC成分PS(f)の一例を示す図である。図6では、横軸が脈拍数(単位は、BPM(30BPMが、0.5Hzに相当する。))を表し、縦軸がスペクトル強度(単位は任意)を表す。この例では、周波数領域に変換された安静時のAC成分PS(f)には、およそ60BPMに大きなピークが見られる。およそ120BPM、およそ180BPMには、その高調波成分が現れている。
iv) 次に、図4のステップS4に示すように、CPU31は、周波数スペクトルの設定された探索周波数範囲において強度ピークを抽出するピーク抽出部として動作する。測定開始時には(安静状態にある被測定者の脈拍数(安静時脈拍数)を求める場合には)、探索周波数範囲は、全周波数範囲(たとえば、30BPM〜300BPM、すなわち、0.5Hz〜5Hz)としてよい。図6の例では、CPU31は、およそ60BPMにおいて周波数スペクトルPS(f)の強度ピークを抽出する。CPU31は、およそ120BPM、およそ180BPMにおける比較的小さな強度ピークについては、およそ60BPMに現れた強度ピークの高調波成分であるとみなして捨てる。次に、CPU31は、抽出された強度ピークの周波数に応じて被測定者の安静時の脈拍数を求める脈拍数算出部として動作することにより、抽出された強度ピークの周波数(図6の場合は、1Hz)に応じて被測定者の安静時脈拍数はおよそ60BPMであると判断する。
v) 次に、図4のステップS5に示すように、CPU31は、前記周波数スペクトルの周波数軸において、強度ピークを探索するための探索周波数範囲を設定する探索範囲設定部として動作する。具体的には、探索範囲設定部として動作するCPU31は、前回の測定において算出された脈拍数(ここでは、安静時の脈拍数(およそ60BPM))に対し予め定められた比率範囲(たとえば、プラスマイナス20%以内)に含まれる値を、次回の測定周期のための前記探索周波数範囲に設定する。たとえば、CPU31は、前回の測定周期において算出された脈拍数(安静時脈拍数)に対してプラスマイナス20%以内にある値の範囲を次回の測定周期のための探索周波数範囲に設定する。図6のように前回の測定周期において算出された脈拍数が60BPMであれば、48BPM〜72BPMの範囲が次回の測定周期のための探索周波数範囲として設定される。
以降、図4のステップS6〜ステップS13にわたる処理ループは、測定開始時から数えて2回目以降の脈拍数測定にかかる処理フローである。ステップS6からステップS13に至る一連の処理が、1回の脈拍数測定で実施される。この一連の処理は、測定終了まで、所定の測定周期で(たとえば、5秒間隔(図7における時間間隔Ts))で実施される。脈拍測定装置1は、測定開始時から数えて2回目以降の脈拍数測定においては、体動センサ部33から出力される運動強度信号にもとづいて、必要に応じて探索周波数範囲を前回の探索周波数範囲からずらすように変化させ、当該探索周波数範囲においてスペクトル強度のピークを抽出し、脈拍数を算出する。
vi) 図4のステップS6に示すように、CPU31は、運動強度取得部として動作し、体動センサ部33から被測定者が行っている運動の強度を表す運動強度信号を取得する。
vii) 次に、ステップS7に示すように、探索範囲設定部として動作するCPU31は、前回の測定周期における被測定者の運動強度と今回の測定周期における被測定者の運動強度とを、運動強度信号に基づいて比較し、今回の測定周期における運動強度が前回の測定周期における運動強度よりも増大したか、変化しなかったか、または、低減したか、を判断する。
図7(a)は、運動強度の時間変化例(3例)と、測定周期との関係を示す図である。横軸は、時刻であり、縦軸は、運動強度信号に基づいて決定された、被測定者の運動の強度である。なお、ここでの運動強度は、体動センサ部33(加速度センサ48)の出力する加速度の各時刻における値でよい。あるいは、運動強度は、加速度センサ48の出力を予め定めた時間間隔にわたって積分することによって得られる値であってもよいし、その他の予め定める計算方法によって体動センサ部33が出力する運動強度信号を処理して得られる値でもよい。たとえば、体動センサ部33(加速度センサ48)の出力から被測定者の歩行ピッチ(走行ピッチ)を求め、該ピッチを運動強度としてもよい。
第1の運動強度時間変化例WLaは、前回の測定周期における運動強度よりも今回の測定周期の運動強度が増大した場合を示す例である。第1の運動強度時間変化例WLaでは、前回の測定周期(時刻t1)において運動強度はla1であったのに対し、今回の測定周期(時刻t2)において運動強度はla2(la2:la2>la1)である。このような場合、CPU31は、図4のステップS7において、今周期の運動強度は直前周期の運動強度から増大するように変化した、と判断する(ステップS7における「YES」)。したがって、処理は、ステップS8へ進む。
第2の運動強度時間変化例WLbは、前回の測定周期における運動強度と今回の測定周期の運動強度との間で変化がなかった場合を示す例である。第2の運動強度時間変化例WLbでは、前回の測定周期(時刻t1)において運動強度はlb1であり、今回の測定周期(時刻t2)においても運動強度はlb2(lb2:lb2=lb1)である。このような場合、CPU31は、図4のステップS7において、今周期の運動強度は直前周期の運動強度から変化しなかった、と判断する(ステップS7における「NO」)。したがって、処理は、ステップS9へ進む。
第3の運動強度時間変化例WLcは、前回の測定周期における運動強度よりも今回の測定周期の運動強度が低減した場合を示す例である。第3の運動強度時間変化例WLcでは、前回の測定周期(時刻t1)において運動強度はlc1であったのに対し、今回の測定周期(時刻t2)において運動強度はlc2(lc2:lc2<lc1)である。このような場合、CPU31は、図4のステップS7において、今周期の運動強度は直前周期の運動強度から低減するように変化した、と判断する(ステップS7における「YES」)。したがって、処理は、ステップS8へ進む。
viii) 図4のステップS8に示すように、CPU31は、探索範囲設定部として動作して、今回の測定周期における運動強度が、前回の測定周期における運動強度よりも大きい場合(図7(a)の運動強度WLaのような場合)、探索周波数範囲を、前回の探索周波数範囲よりも高周波数側(高BPM側)にずらす(シフトさせる)。
逆に、CPU31は、今回の測定周期における運動強度が、前回の測定周期における運動強度よりも小さい場合(図7(a)の運動強度WLcのような場合)、探索周波数範囲を、前回の探索周波数範囲よりも低周波数側(低BPM側)にずらす(シフトさせる)。
ステップS9に示すように、CPU31は、今回の測定周期における運動強度が、前回の測定周期における運動強度から変化していない場合(図7(a)の運動強度WLbのような場合)、探索周波数範囲を、前回の探索周波数範囲から変化させない。
図8は、図4のステップS8およびステップS9による探索周波数範囲の変更(維持)の様子を表す図である。横軸は、脈拍数(BPM)を表し、縦軸は、スペクトル強度(単位は任意)を表す。図7(b)に示す時間領域の脈波信号(脈波AC成分PS(t))のうち、t=t2−Tdからt=t2までの脈波信号PS(t)を周波数領域に変換したものが、図8における周波数スペクトルPS(f)である。
図8の例では、前回の測定周期の探索周波数範囲は、周波数範囲SR1である。周波数範囲SR1は、下限周波数fL1および上限周波数fH1で規定される周波数幅(fH1−fL1)の周波数範囲であって、前回の脈拍測定では当該範囲において強度ピークの抽出が行われたものとする。
図4のステップS8において、CPU31は、今回の測定周期における運動強度が、前回の測定周期における運動強度よりも大きい場合(図7(a)の運動強度WLaのような場合)、探索周波数範囲を、前回の探索周波数範囲SR1よりも高周波数側(高BPM側)の周波数範囲SR2aにずらす(シフトさせる)。これにより、今回の測定周期における探索周波数範囲は、下限周波数fL2a(fL2a=fL1+dPb)および上限周波数fH2a(fH2a=fH1+dPt)で規定される周波数幅(fH2a−fL2a)の周波数範囲となる。ここで、dPb=dPtとしてよく、その場合、今回の探索周波数範囲の周波数幅は、前回の探索周波数範囲の周波数幅と同一となる。また、CPU31は、今回の測定周期における運動強度と前回の測定周期における運動強度との差(図7(a)における運動強度の差(la2−la1)が大きければ大きいほどに、探索周波数範囲のシフト量(図8におけるdPtおよびdPb)を大きくしてもよい。そうすることで、運動負荷の増大に伴って増加であろう脈拍数をより確実に追跡することが可能となる。
逆に、CPU31は、今回の測定周期における運動強度が、前回の測定周期における運動強度よりも小さい場合(図7(a)の運動強度WLcのような場合)、探索周波数範囲を、前回の探索周波数範囲SR1よりも低周波数側(低BPM側)の周波数範囲SR2cにずらす(シフトさせる)。これにより、今回の測定周期における探索周波数範囲は、下限周波数fL2c(fL2c=fL1−dMb)および上限周波数fH2c(fH2c=fH1−dMt)で規定される周波数幅(fH2c−fL2c)の周波数範囲となる。ここでも、dMb=dMtとしてよく、その場合、今回の探索周波数範囲の周波数幅は、前回の探索周波数範囲の周波数幅と同一となる。また、CPU31は、今回の測定周期における運動強度と前回の測定周期における運動強度との差(図7(a)における運動強度の差(lc1−lc2)が大きければ大きいほどに、探索周波数範囲のシフト量(図8におけるdMtおよびdMb)を大きくしてもよい。そうすることで、運動負荷の減少に伴って低下するであろう脈拍数をより確実に追跡することが可能となる。
図4のステップS9においては、CPU31は、今回の測定周期における運動強度が、前回の測定周期における運動強度から変化していない場合(図7(a)の運動強度WLbのような場合)、探索周波数範囲を、前回の探索周波数範囲SR1から変化させない。
ix) 図4のステップS10に示すように、CPU31は、データ取得部として動作することにより、今回の測定周期にための脈波信号(脈波AC成分PS(t))の時系列データを記憶部32から取得する。たとえば、CPU31は、図7(b)のような脈波信号(脈波AC成分PS(t))の時系列データが記憶部32に記憶されている場合、時刻t2−Tdから今回の測定周期の時刻t2までの脈波信号(脈波AC成分PS(t))の時系列データを記憶部32から取得する。
x) 次に、図4のステップS11に示すように、CPU31は、周波数変換部として動作することにより、記憶部32に記憶された時間領域の脈波信号(脈波AC成分PS(t))を周波数領域に変換し、脈波信号(脈波AC成分PS(t))の周波数スペクトル(PS(f))を求める。たとえば、CPU31は、ステップS10で取得した予め定められた期間Tdの脈波信号(脈波AC成分PS(t))の時系列データについて高速フーリエ変換(FFT)を行い、図8に示すような脈波信号の周波数スペクトルPS(f)を導出する。
xi) そして、図4のステップS12に示すように、CPU31は、ピーク抽出部として動作することにより、ステップS8またはステップS9で設定された今回の測定周期の探索周波数範囲(SR2a、SR2b、または、SR2c)において、周波数スペクトルの強度ピーク(極大点)を抽出する。次に、CPU31は、脈拍数算出部として動作することにより、抽出された強度ピークの周波数に応じて被測定者の安静時脈拍数を算出する。
xii) ステップS13において、CPU31は、脈拍測定を終了するか否かを判断し、脈拍測定を継続すると判断した場合には、ステップS6に戻って次回の測定周期のための処理を行う。
以上のように、一実施形態における脈拍測定装置1は、被測定者が行っている運動の強度に基づいて脈拍の変動の傾向を予測し、予測した脈拍変動の方向性を考慮して、前回の探索周波数範囲を高周波数側あるいは低周波数側にずらしたり、前回の探索周波数範囲をそのまま維持したりして、脈拍に起因したスペクトル強度のピークを周波数領域の脈波信号から抽出する。そうすることで、たとえ被測定者が運動を行うことにより外乱成分が脈波信号に重畳されるような場合であっても、外乱成分に起因したスペクトル強度のピークを脈拍に起因したスペクトル強度のピークであると誤認識することがなくなり(少なくとも、誤認識する頻度が低減され)、被測定者が安静状態になくても、その被測定者の脈拍数を正しく測定することが可能となっている。
一実施形態における脈拍測定装置1は、非心電的に取得された脈波信号の周波数スペクトル強度分布に基づいて被測定者の脈拍数を算出する脈拍測定装置である。非心電的とは、たとえば、光電方式を指すが、これに限定されない。非心電的な方法には、光電方式のほか、圧電方式等が含まれる。
一実施形態における脈拍測定装置1は、脈波信号として、光電出力PDCのうち被測定者の脈拍数として想定する範囲(30BPM〜300BPM)の周期で変動する成分を取りだして利用している。しかしながら、光電出力PDCをそのまま脈波信号として利用してもよい。
また、上述の脈拍測定の方法をコンピュータに実行させるためのプログラムとして構築してもよい。
また、そのようなプログラム(脈拍測定プログラム)をCD−ROMなどのコンピュータ読み取り可能な記録媒体に記録して配布できるようにしても良い。上記脈拍測定プログラムを汎用コンピュータにインストールすることで、汎用コンピュータによって上記脈拍測定方法を実行することが可能である。
また、記憶部32に記憶されているプログラムを、メモリその他の非一時的なコンピュータ読み取り可能な記録媒体(メモリ、ハードデイスクドライブ、光ディスクなど)にエンコードしておき、汎用コンピュータに上述の脈拍測定方法を実行させてもよい。また、プログラムは、インターネット等を通じて配信されてもよい。
また、上の例では、CPU31は、周波数領域への変換として高速フーリエ変換(FFT)を行ったが、これに限られるものではない。時間領域の光電信号PDCを周波数領域に変換できるものであれば、他の変換方式を採用してもよい。
また、CPU31として、上記脈拍測定の方法を実行する専用のハードウェアロジック回路を用いてもよい。つまり、データ取得部、運動強度取得部、探索範囲設定部、ピーク抽出部、脈拍数算出部の少なくともいずれか1つは、専用のハードウェア回路によって実現されてもよい。
また、上の例では、図4のステップS1において被測定者が安静状態にあると判断されたとき、図4のステップS4において、脈波信号の周波数スペクトルに含まれた最大の強度ピークが示す周波数を、被測定者の安静状態の脈拍数として求めた。しかしながら、これに限られるものではない。脈波信号(脈波AC成分PS(t))のピークまたは谷の数をカウントして、脈波信号(脈波AC成分PS(t))の変動の繰り返し数から1分あたりの変動回数を求め、それに基づいて被測定者の安静状態の脈拍数を求めてもよい。
以上の実施形態は例示であり、この発明の範囲から離れることなく様々な変形が可能である。
1 脈拍測定装置
10 本体
15 脈波センサ部
16 発光素子
17 受光素子
31 CPU
32 記憶部
33 体動センサ部
54 発光素子
56 受光素子

Claims (6)

  1. 被測定者の脈波を脈波センサによって検知して脈拍を表す脈波信号を取得するデータ取得部と、
    前記被測定者の動きを体動センサによって検知して前記被測定者が行っている運動の強度を表す運動強度信号を取得する運動強度取得部と、
    前記脈波信号を記憶する記憶部と、
    前記記憶部に記憶された時間領域の前記脈波信号を周波数領域に変換し、前記脈波信号の周波数スペクトルを求める周波数変換部と、
    前記周波数スペクトルの周波数軸において、強度ピークを探索するための探索周波数範囲を設定する探索範囲設定部と、
    前記周波数スペクトルの前記設定された探索周波数範囲において強度ピークを抽出するピーク抽出部と、
    抽出された前記強度ピークの周波数に応じて前記被測定者の脈拍数を求める脈拍数算出部と、を有し、
    前記探索範囲設定部は、前記運動強度信号が示す運動強度に応じて、前記探索周波数範囲を変更することを特徴とする、脈拍測定装置。
  2. 請求項1に記載の脈拍測定装置において、
    前記周波数変換部、前記運動強度取得部、前記探索範囲設定部、前記ピーク抽出部、および、前記脈拍数算出部は、予め定められた周期で処理を繰り返し、
    第1の周期において前記脈拍数算出部が前記被測定者の脈拍数として第1の値を算出した場合には、前記探索範囲設定部は、前記第1の値に対し予め定められた比率範囲に含まれる値を、前記第1の周期に続く第2の周期のための前記探索周波数範囲に設定することを特徴とする、脈拍測定装置。
  3. 請求項1または2に記載の脈拍測定装置において、
    前記運動強度取得部が第3の周期に続く第4の周期において取得した運動強度が前記第3の周期において取得した運動強度よりも大きい場合には、前記探索範囲設定部は、前記第3の周期のための前記探索周波数範囲よりも高周波数側にシフトされた周波数範囲を前記第4の周期のための前記探索周波数範囲に設定することを特徴とする、脈拍測定装置。
  4. 請求項3に記載の脈拍測定装置において、
    前記探索範囲設定部は、前記第3の周期のための前記探索周波数範囲の周波数幅と同じ幅を有するように、前記第4の周期のための前記探索周波数範囲を設定することを特徴とする、脈拍測定装置。
  5. 脈拍測定装置がする被測定者の脈拍数の測定方法であって、
    脈波センサによって、前記被測定者の脈拍を表す脈波信号を取得するデータ取得ステップと、
    前記脈波信号を記憶部に記憶する記憶ステップと、
    前記記憶部に記憶された時間領域の前記脈波信号を周波数領域に変換し、前記脈波信号の周波数スペクトルを求める周波数変換ステップと、
    体動センサによって、前記被測定者が行っている運動の強度を表す運動強度信号を取得する運動強度取得ステップと、
    前記周波数スペクトルの周波数軸において、強度ピークを探索するための探索周波数範囲を設定する探索範囲設定ステップと、
    前記周波数スペクトルの前記設定された探索周波数範囲において強度ピークを抽出するピーク抽出ステップと、
    抽出された前記強度ピークの周波数に応じて前記被測定者の脈拍数を求める脈拍数算出ステップと、を有し、
    前記探索範囲設定ステップは、前記運動強度信号が示す運動強度に応じて、前記探索周波数範囲を変更するステップを含むことを特徴とする、脈拍数の測定方法。
  6. 脈拍測定装置のコンピュータが実行可能な脈拍数測定用コンピュータ・プログラムであって、
    前記コンピュータに、
    脈波センサによって、被測定者の脈拍を表す脈波信号を取得するデータ取得ステップと、
    前記脈波信号を記憶部に記憶する記憶ステップと、
    前記記憶部に記憶された時間領域の前記脈波信号を周波数領域に変換し、前記脈波信号の周波数スペクトルを求める周波数変換ステップと、
    体動センサによって、前記被測定者が行っている運動の強度を表す運動強度信号を取得する運動強度取得ステップと、
    前記周波数スペクトルの周波数軸において、強度ピークを探索するための探索周波数範囲を設定する探索範囲設定ステップと、
    前記周波数スペクトルの前記設定された探索周波数範囲において強度ピークを抽出するピーク抽出ステップと、
    抽出された前記強度ピークの周波数に応じて前記被測定者の脈拍数を求める脈拍数算出ステップと、を実行させる脈拍測定プログラムであって、
    前記探索範囲設定ステップは、前記運動強度信号が示す運動強度に応じて、前記探索周波数範囲を変更するステップを含むことを特徴とする、脈拍数測定用コンピュータ・プログラム。
JP2012201910A 2012-09-13 2012-09-13 脈拍測定装置、ならびに、脈拍測定方法および脈拍測定プログラム Active JP5987578B2 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP2012201910A JP5987578B2 (ja) 2012-09-13 2012-09-13 脈拍測定装置、ならびに、脈拍測定方法および脈拍測定プログラム
PCT/JP2013/070750 WO2014041913A1 (ja) 2012-09-13 2013-07-31 脈拍測定装置、ならびに、脈拍測定方法および脈拍測定プログラム
CN201380045775.4A CN104602595B (zh) 2012-09-13 2013-07-31 脉搏测定装置和脉搏测定方法
DE112013004467.6T DE112013004467T5 (de) 2012-09-13 2013-07-31 Pulsmessvorrichtung, Pulsmessverfahren und Pulsmessprogramm
US14/641,546 US20150173627A1 (en) 2012-09-13 2015-03-09 Pulse measurement device, pulse measurement method, and pulse measurement program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2012201910A JP5987578B2 (ja) 2012-09-13 2012-09-13 脈拍測定装置、ならびに、脈拍測定方法および脈拍測定プログラム

Publications (3)

Publication Number Publication Date
JP2014054447A true JP2014054447A (ja) 2014-03-27
JP2014054447A5 JP2014054447A5 (ja) 2015-09-17
JP5987578B2 JP5987578B2 (ja) 2016-09-07

Family

ID=50278041

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012201910A Active JP5987578B2 (ja) 2012-09-13 2012-09-13 脈拍測定装置、ならびに、脈拍測定方法および脈拍測定プログラム

Country Status (5)

Country Link
US (1) US20150173627A1 (ja)
JP (1) JP5987578B2 (ja)
CN (1) CN104602595B (ja)
DE (1) DE112013004467T5 (ja)
WO (1) WO2014041913A1 (ja)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017136165A (ja) * 2016-02-02 2017-08-10 富士通株式会社 センサ情報処理装置、センサユニット、及び、センサ情報処理プログラム
JP2017136164A (ja) * 2016-02-02 2017-08-10 富士通株式会社 センサ情報処理装置、センサユニット、及び、センサ情報処理プログラム
JP2017136163A (ja) * 2016-02-02 2017-08-10 富士通株式会社 センサ情報処理装置、センサユニット、及び、センサ情報処理プログラム
JP2018040750A (ja) * 2016-09-09 2018-03-15 フロンティアシステム株式会社 マーキング装置
KR20180067348A (ko) * 2016-12-12 2018-06-20 경북대학교 산학협력단 Ppg를 이용한 맥박수 추정방법 및 장치
JP2018138138A (ja) * 2017-02-24 2018-09-06 富士通株式会社 心拍数推定方法および心拍数推定装置
CN108685569A (zh) * 2017-03-30 2018-10-23 瑞萨电子株式会社 脉搏测量设备、脉搏测量方法以及非暂态计算机可读介质
JP7171193B2 (ja) 2018-01-16 2022-11-15 フクダ電子株式会社 生体信号処理装置およびその制御方法

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5991100B2 (ja) 2012-09-13 2016-09-14 オムロンヘルスケア株式会社 脈拍測定装置、脈拍測定方法、および脈拍測定プログラム
US9814400B1 (en) * 2015-05-26 2017-11-14 Verily Life Sciences Llc Method for improving accuracy of pulse rate estimation
WO2016187847A1 (zh) * 2015-05-27 2016-12-01 深圳市长桑技术有限公司 一种信号获取方法与系统
JP2017113191A (ja) * 2015-12-22 2017-06-29 富士通株式会社 電子機器及び脈拍数算出プログラム
CN111493846B (zh) * 2019-01-31 2023-03-21 深圳市爱都科技有限公司 一种血氧饱和度和心率值的计算方法和穿戴设备
CN109924960A (zh) * 2019-01-31 2019-06-25 深圳市爱都科技有限公司 一种血氧饱和度、心率值和压力等级的计算方法和穿戴设备
CN111329481A (zh) * 2020-03-03 2020-06-26 中国科学院深圳先进技术研究院 生理参数确定方法、装置、生理参数检测设备及介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08289876A (ja) * 1995-02-20 1996-11-05 Seiko Epson Corp 脈拍計
JPH10258038A (ja) * 1997-03-17 1998-09-29 Seiko Epson Corp 脈波検出装置および脈拍計
JP2005160640A (ja) * 2003-12-01 2005-06-23 Denso Corp 生体状態検出装置
JP2007330431A (ja) * 2006-06-14 2007-12-27 Mitsuba Corp 生体情報判定システム及び生体情報判定方法並びに生体情報判定プログラム
JP2012170701A (ja) * 2011-02-23 2012-09-10 Seiko Epson Corp 拍動検出装置

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101039617A (zh) * 2004-10-15 2007-09-19 普尔塞特拉瑟技术有限公司 用于生理脉冲测量的光学输入信号的运动消除
JP5181477B2 (ja) * 2007-01-11 2013-04-10 ヤマハ株式会社 フィットネス運動状態表示装置
JP5139106B2 (ja) * 2008-02-12 2013-02-06 株式会社東芝 脈波間隔計測装置及び計測方法
JP2009297106A (ja) * 2008-06-10 2009-12-24 Univ Fukuoka 心音測定装置、至適運動強度測定装置、心音測定方法および心音測定プログラム
CN101991410B (zh) * 2009-08-31 2012-09-19 深圳市理邦精密仪器股份有限公司 一种脉率搜索和计算方法
US8718980B2 (en) * 2009-09-11 2014-05-06 Qualcomm Incorporated Method and apparatus for artifacts mitigation with multiple wireless sensors
CN102270264B (zh) * 2010-06-04 2014-05-21 中国科学院深圳先进技术研究院 生理信号质量评估系统及方法
CN102028457B (zh) * 2010-11-24 2012-10-03 北京麦邦光电仪器有限公司 脉率测量方法及指环式脉率测量仪
US9307943B2 (en) * 2011-07-01 2016-04-12 Seiko Epson Corporation Biological information processing device

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08289876A (ja) * 1995-02-20 1996-11-05 Seiko Epson Corp 脈拍計
JPH10258038A (ja) * 1997-03-17 1998-09-29 Seiko Epson Corp 脈波検出装置および脈拍計
JP2005160640A (ja) * 2003-12-01 2005-06-23 Denso Corp 生体状態検出装置
JP2007330431A (ja) * 2006-06-14 2007-12-27 Mitsuba Corp 生体情報判定システム及び生体情報判定方法並びに生体情報判定プログラム
JP2012170701A (ja) * 2011-02-23 2012-09-10 Seiko Epson Corp 拍動検出装置

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017136165A (ja) * 2016-02-02 2017-08-10 富士通株式会社 センサ情報処理装置、センサユニット、及び、センサ情報処理プログラム
JP2017136164A (ja) * 2016-02-02 2017-08-10 富士通株式会社 センサ情報処理装置、センサユニット、及び、センサ情報処理プログラム
JP2017136163A (ja) * 2016-02-02 2017-08-10 富士通株式会社 センサ情報処理装置、センサユニット、及び、センサ情報処理プログラム
US10667758B2 (en) 2016-02-02 2020-06-02 Fujitsu Limited Sensor information processing apparatus
JP2018040750A (ja) * 2016-09-09 2018-03-15 フロンティアシステム株式会社 マーキング装置
KR20180067348A (ko) * 2016-12-12 2018-06-20 경북대학교 산학협력단 Ppg를 이용한 맥박수 추정방법 및 장치
KR102059685B1 (ko) * 2016-12-12 2019-12-26 경북대학교 산학협력단 Ppg를 이용한 맥박수 추정방법 및 장치
JP2018138138A (ja) * 2017-02-24 2018-09-06 富士通株式会社 心拍数推定方法および心拍数推定装置
CN108685569A (zh) * 2017-03-30 2018-10-23 瑞萨电子株式会社 脉搏测量设备、脉搏测量方法以及非暂态计算机可读介质
JP2018166883A (ja) * 2017-03-30 2018-11-01 ルネサスエレクトロニクス株式会社 脈拍計測装置、脈拍計測方法、及びプログラム
JP7171193B2 (ja) 2018-01-16 2022-11-15 フクダ電子株式会社 生体信号処理装置およびその制御方法

Also Published As

Publication number Publication date
CN104602595A (zh) 2015-05-06
US20150173627A1 (en) 2015-06-25
WO2014041913A1 (ja) 2014-03-20
CN104602595B (zh) 2016-09-28
DE112013004467T5 (de) 2015-05-28
JP5987578B2 (ja) 2016-09-07

Similar Documents

Publication Publication Date Title
JP5987578B2 (ja) 脈拍測定装置、ならびに、脈拍測定方法および脈拍測定プログラム
JP2014054447A5 (ja)
JP6079075B2 (ja) 脈拍測定装置、ならびに、脈拍測定方法および脈拍測定プログラム
JP5991100B2 (ja) 脈拍測定装置、脈拍測定方法、および脈拍測定プログラム
JP5855004B2 (ja) 腕部装着式血圧計
US20190254590A1 (en) Method and apparatus for providing biofeedback during meditation exercise
JP6126220B2 (ja) 生体状態推定装置
US7566306B2 (en) Biological information processing apparatus and operation program product for the same
JP6406624B2 (ja) 生体情報測定装置およびこの装置を備える生体情報測定システム
EP2821008A1 (en) Bioinformation processing apparatus and signal processing method
JP5562805B2 (ja) 脈拍数測定方法及び血中酸素飽和度測定方法
JP2016016203A (ja) 生体情報検出装置
JP5760876B2 (ja) 心房細動判定装置、心房細動判定方法およびプログラム
JP2014171589A (ja) 心房細動解析装置およびプログラム
JP6060563B2 (ja) 心房細動判定装置、心房細動判定方法およびプログラム
JP5958222B2 (ja) 生体情報測定装置及び該装置におけるパラメータ設定方法
JP2018149183A (ja) 血圧データ処理装置、血圧データ処理方法、およびプログラム
JP2016000149A (ja) 脈拍測定装置、脈拍測定方法、および、脈拍測定プログラム
JP2016000153A (ja) 脈拍測定装置、脈拍測定方法、および、脈拍測定プログラム
JP2014212796A (ja) 脈波計測装置
JP2022160608A (ja) 生体信号処理装置およびその制御方法
JP2019069236A (ja) 血圧測定装置、腕時計端末、及び血圧測定方法
JP6098673B2 (ja) 心房細動判定装置、心房細動判定装置の作動方法およびプログラム
US20210038162A1 (en) Method, device, and program for estimating pulse rate

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20150730

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150803

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20160712

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160725

R150 Certificate of patent or registration of utility model

Ref document number: 5987578

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150