JPH10258039A - 脈拍計 - Google Patents
脈拍計Info
- Publication number
- JPH10258039A JPH10258039A JP9063486A JP6348697A JPH10258039A JP H10258039 A JPH10258039 A JP H10258039A JP 9063486 A JP9063486 A JP 9063486A JP 6348697 A JP6348697 A JP 6348697A JP H10258039 A JPH10258039 A JP H10258039A
- Authority
- JP
- Japan
- Prior art keywords
- data
- pulse
- pulse wave
- body movement
- frequency
- 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
Links
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
ックはDB2'、DB3'、DB4'およびDB1になる。こ
こで、データブロックDB2'〜DB4'には体動成分が重
畳しており、データブロックDB1には体動成分が重畳
していない(d)。このため、期間T6において、バッ
ファ制御信号は、データブロックDB2'〜DB4'を0デ
ータに置換して出力し、これに続いてデータブロックD
B1を出力するようにバッファを制御する(e)。ここ
で、0データは脈波データの中心値である。したがっ
て、体動成分が重畳したデータブロックDB2'〜DB4'
を周波数成分の解析に影響を与えないデータに置換する
ことができる。
Description
管理や健康管理などを行うための脈拍計に関するもので
ある。特に、安静時および運動時のいずれの状態におい
ても、高い精度をもって脈拍数を計測するための信号処
理技術に関するものである。
計測できれば、利用者の運動量の管理や健康管理(危険
防止)を行うことができることから、腕などに装着した
まま脈拍数の計測を行うことができる携帯用の脈拍計が
案出されている。かかる携帯用の脈拍計では、光学セン
サなどを用いて脈波信号を計測し、この脈波信号から脈
拍に相当する信号を抽出して脈拍数を求めている。但
し、ジョギング中に計測した脈波信号には、その体動に
起因する信号成分も含まれている。このため、脈波成分
から体動成分を除去する必要がある。
開発されている。まず、特開平7−227383号公報
に開示されている脈拍計にあっては、体動センサによっ
て検出された体動信号に基づいて体動の有無を判定す
る。体動が無いと判定された場合には、脈波センサによ
って検出された脈波信号を矩形波に波形整形しこの矩形
波に基づいて脈拍数を算出する。一方、体動が有りと判
定された場合には、脈波信号にFFT処理を施した解析
結果から体動信号にFFT処理を施した解析結果を減算
して体動周波数成分が除去された脈波周波数成分を求
め、この演算結果に基づいて脈拍数を算出する。この脈
拍計にあっては、体動が無い場合には、FFT処理を行
う必要がないので消費電力を低減することができる。
ている脈拍計にあっては、体動センサによって検出され
た体動信号に基づいて体動の有無を判定する。体動が無
いと判定された場合には、脈波センサによって検出され
た脈波信号にFFT処理を施して周波数解析を行いこの
解析結果に基づいて脈拍数を算出する。一方、体動が有
ると判定された場合には、脈波信号にFFT処理を施し
た解析結果から体動信号にFFT処理を施した解析結果
を減算して体動周波数成分が除去された脈波周波数成分
を求め、この演算結果に基づいて脈拍数を算出する。こ
の脈拍計にあっては、脈波信号を矩形波に波形整形を行
う回路を省略できるので、構成を簡易することができ
る。
来の脈拍計においては、FFT(Fast Fourier Transfo
rm:高速フーリエ変換)処理を行っているが、FFT処
理した結果得られる最も低い周波数は、分析時間の逆数
で決定される。すなわち、FFT処理は、一定のサンプ
ル数を分析単位として行われる。
から体動無しに切り替わった場合、体動が無くなった後
半の期間においては、本来、脈波信号の周波数解析結果
に基づいて脈拍数の算出を行うべきであるにもかかわら
ず、脈波の解析結果から体動の解析結果を減算した脈波
周波数成分に基づいて脈拍数の算出が行われる。この場
合には、正確な脈拍数の検出が、行われない場合があっ
た。このように従来の脈拍計にあっては、利用者が運動
をやめた場合であっても、すぐには脈拍数の算出方法を
切り替えることができず、応答性に問題があった。
本来、体動が無いにも拘わらず脈波成分から体動成分が
差し引かれ、正確な脈波を検出することができない。こ
のような異常状態から通常状態に戻った場合には、すぐ
に正確な脈波を検出することが望ましい。しかし、上述
したように脈拍数の算出方法の切り替えには時間がかか
るので、通常状態に復帰した後にも、不正確な脈拍数を
検出してしまうといった問題があった。
のであり、応答性を改善した脈拍計を提供することを目
的とする。
め、請求項1に記載の発明にあっては、一定のサンプル
数に対応するデータブロック毎に脈拍を計測する脈拍計
において、生体の体動を検出して体動データを出力する
体動検出手段と、前記体動データに基づいて体動の有無
を前記データブロック毎に判定する判定手段と、前記生
体の脈動を検出して脈波データを出力する脈波検出手段
と、前記脈波データを記憶するとともに、前記判定手段
の判定結果が体動有りから体動無しに切り替わった後の
所定期間においては、体動有りと判定された期間中に記
憶した前記脈波データを周波数解析に影響を与えないデ
ータに置換して出力し、前記所定期間以外の期間におい
ては記憶した前記脈波データを出力する置換手段と、前
記置換手段の出力データに対して複数の前記データブロ
ック単位で周波数解析を施し、次の周波数解析を1デー
タブロックずらした前記置換手段の出力データに対して
施すことにより、前記データブロック毎に解析結果を得
る第1の演算手段と、前記体動データに対して複数の前
記データブロック単位で周波数解析を施し、次の周波数
解析を1データブロックずらした前記体動データに対し
て施すことにより、前記データブロック毎に解析結果を
得る第2の演算手段と、前記判定手段によって体動有り
と判定された場合には、前記第1の演算手段の周波数解
析に基づいて脈拍に対応する周波数を抽出し、前記判定
手段によって体動無しと判定された場合には、前記第1
の演算手段と前記第2の演算手段の周波数解析結果に基
づいて、脈拍に対応する周波数を抽出する脈波成分抽出
手段と、前記脈波成分抽出手段によって抽出された前記
脈拍の周波数から脈拍数を算出する脈拍数算出手段と、
前記脈拍数を表示する表示手段とを備えたことを特徴と
する。
前記第2の演算手段の周波数解析結果から、前記生体の
体動のピッチを算出するピッチ算出手段と、前記判定手
段によって体動有りと判定された場合には前記ピッチを
表示するように前記表示手段を制御し、前記判定手段に
よって体動無しと判定された場合には、前記ピッチをに
表示しないように前記表示手段を制御する表示制御手段
とを備えたことを特徴とする。
前記第1の演算手段および前記第2の演算手段で行う周
波数解析処理は高速フーリエ変換処理であり、前記デー
タブロックのサンプル数は、高速フーリエ変換演算対象
データに対して、1/4もしくは1/2であることを特
徴とする。
機能ブロック図である。図において、脈拍計fは、一定
のサンプル数に対応するデータブロック毎に脈拍を計測
するようになっている。f1は体動検出手段であって、
生体の体動を検出して体動データを出力する。体動検出
手段f1は、例えば加速度センサによって構成される。
また、f2は判定手段であって、体動検出手段f1によ
って検出される体動データに基づいて体動の有無をデー
タブロック毎に判定する。
の脈動を検出して脈波データを出力する。脈波検出手段
f3としては各種のものがあるが、例えば、光電式脈波
センサが好適である。また、f4は置換手段であって、
脈波データを記憶するとともに、判定手段f2の判定結
果が体動有りから体動無しに切り替わった後の所定期間
においては、体動有りと判定された期間中に記憶した脈
波データを周波数解析に影響を与えないデータに置換し
て出力し、所定期間以外の期間においては記憶した脈波
データを出力する。また、f5は第1の演算手段であっ
て、置換手段f4の出力データに対して複数のデータブ
ロック単位で周波数解析を施し、次の周波数解析を1デ
ータブロックずらした置換手段f4の出力データに対し
て施すことにより、データブロック毎に解析結果を得
る。これにより、第1の演算手段f2は、体動有りから
体動無しに切り替わった後の所定期間において、実質的
に体動無しに対応する脈波データに基づいて周波数解析
を行うことができる。
動データに対して複数の前記データブロック単位で周波
数解析を施し、次の周波数解析を1データブロックずら
した前記体動データに対して施すことにより、前記デー
タブロック毎に解析結果を得る。また、f7は脈波成分
抽出手段であって、判定手段f2によって体動有りと判
定された場合には、第1の演算手段f5の周波数解析に
基づいて脈拍に対応する周波数を抽出し、判定手段f2
によって体動無しと判定された場合には、第1の演算手
段f5と第2の演算手段f7の周波数解析結果に基づい
て、脈拍に対応する周波数を抽出する。また、f8は脈
拍数検出手段であって、脈波成分抽出手段f7によって
抽出された前記脈拍の周波数から脈拍数を算出する。
出手段f8によって検出された脈拍数を表示する。ま
た、f10はピッチ算出手段であって、第2の演算手段
f6の周波数解析結果から、前記生体の体動のピッチを
算出する。これにより、ランニング時のピッチを求める
ことができる。また、f11は表示制御手段であって、
判定手段f2によって体動有りと判定された場合にはピ
ッチを表示するように表示手段f9を制御し、判定手段
f2によって体動無しと判定された場合には、ピッチを
表示しないように表示手段f9制御する。
説明図である。図2において、本例の脈拍計1(携帯用
脈波計測装置)は、腕時計構造を有する装置本体10
と、この装置本体10に接続されるケーブル20と、こ
のケーブル20の先端側に設けられた脈波検出用センサ
ユニット30(脈波信号検出用センサ)とから大略構成
されている。ケーブル20の先端側にはコネクタピース
80が構成されており、このコネクタピース80は、装
置本体10の6時の側に構成されているコネクタ部70
に対して着脱自在である。装置本体10には、腕時計に
おける12時方向から腕に巻きついてその6時方向で固
定されるリストバンド12が設けられ、このリストバン
ド12によって、装置本体10は、腕に着脱自在であ
る。脈波検出用センサユニット30は、センサ固定用バ
ンド40によって遮光されながら人差し指の根本に装着
される。このように、脈波検出用センサユニット30を
指の根本に装着すると、ケーブル20が短くて済むの
で、ケーブル20は、ランニング中に邪魔にならない。
また、掌から指先までの体温の分布を計測すると、寒い
ときには、指先の温度が著しく低下するのに対し、指の
根本の温度は比較的低下しない。従って、指の根本に脈
波検出用センサユニット30を装着すれば、寒い日に屋
外でランニングしたときでも、脈拍数などを正確に計測
できる。
ス11(本体ケース)を備えており、この時計ケース1
1の表面側には、現在時刻や日付に加えて、走行時や歩
行時のピッチ、および脈拍数などの脈波情報などを表示
するELバックライト付きの液晶表示装置13が構成さ
れている。液晶表示装置13には、セグメント表示領域
の他、ドット表示領域が構成されており、ドット表示領
域では、各種の情報をグラフィック表示可能である。
91を利用して、体の動きを体動信号TSとして検出す
る体動検出用センサ装置90が内蔵されている。また、
時計ケース11の内部には、脈波検出用センサユニット
30が計測した脈波信号MSに基づいて脈拍数の変化な
どを求めるとともに、それを液晶表示装置13に表示す
るために、各種の制御やデータ処理を行うマイクロコン
ピュータなどからなる制御部が構成されている。制御部
には計時回路も構成されており、通常時刻、ラップタイ
ム、スプリットタイムなども液晶表示装置13に表示で
きるようになっている。また、時計ケース11の外周部
には、時刻合わせや表示モードの切換などの外部操作を
行うためのボタンスイッチ111〜115が構成されて
いる。
LED、フォトトランジスタなどから構成される。LE
Dから照射された光は血液によって反射され、その反射
光がフォトトランジスタによって受光されるようになっ
ている。また、LED31としては、InGaN系(イ
ンジウム−ガリウム−窒素系)の青色LEDが好適であ
る。青色LEDの発光スペクトルは、例えば450nm
に発光ピークを有し、その発光波長域は、350nmか
ら600nmまでの範囲にある。この場合には、かかる
発光特性を有するLEDに対応させてフォトトランジス
タとして、GaAsP系(ガリウム−砒素−リン系)の
フォトトランジスタを用いればよい。このフォトトラン
ジスタの受光波長領域は、例えば、主要感度領域が30
0nmから600nmまでの範囲にあって、300nm
以下にも感度領域がある。このような青色LEDとフォ
トトランジスタとを組み合わせると、その重なり領域で
ある300nmから600nmまでの波長領域におい
て、脈波が検出される。この場合には、以下の利点があ
る。
が700nm以下の光は、指の組織を透過しにくい傾向
があるため、外光がセンサ固定用バンドで覆われていな
い指の部分に照射されても、指の組織を介してフォトト
ランジスタ32まで到達せず、検出に影響を与えない波
長領域の光のみがフォトトランジスタ32に達する。一
方、300nmより低波長領域の光は、皮膚表面でほと
んど吸収されるので、受光波長領域を700nm以下と
しても、実質的な受光波長領域は、300nm〜700
nmとなる。したがって、指を大掛かりに覆わなくと
も、外光の影響を抑圧することができる。また、血液中
のヘモグロビンは、波長が300nmから700nmま
での光に対する吸光係数が大きく、波長が880nmの
光に対する吸光係数に比して数倍〜約100倍以上大き
い。したがって、この例のように、ヘモグロビンの吸光
特性に合わせて、吸光特性が大きい波長領域(300n
mから700nm)の光を検出光として用いると、その
検出値は、血量変化に応じて感度よく変化するので、血
量変化に基づく脈波信号MSのS/N比を高めることが
できる。
図3は本実施形態に係わる脈拍計のブロック図である。
30によって検出された脈波信号MSは、アンプ201
によって所定のレベルの増幅された後、A/D変換器2
02を介してデジタル信号に変換され脈波データMDと
してバッファ203に供給され、そこに格納されるよう
になっている。この脈波データMDは1サンプル当たり
8ビットで表され、0レベルを中心として±127レベ
ルを指示する。
検出された体動信号TSは、アンプ204によって所定
のレベルの増幅された後、A/D変換器205を介して
デジタル信号に変換され体動データTDとして出力され
る。この体動データTDは1サンプル当たり8ビットで
表される。なお、この例において、A/D変換器20
2,205のサンプリング周波数は、8Hzに設定され
ており、16秒当たり128サンプルが得られるように
なっている。
を順次格納し、後述するバッファ制御信号BCSに基づ
いて、所定期間内に供給された脈波データMDあるい
は、その一部を0データに置換したデータを脈波信号用
FFT回路207に出力する。
ッファ203から出力される脈波データMDに対してF
FT処理を施し、脈波信号MSの各周波数脈成分を示す
脈波解析データMKDを生成する。このFFT処理を図
4を参照して説明する。図において、時刻t0は脈波の
計測開始時刻であり、DB1,DB2,…DB5,…は脈
波データMDを32サンプル単位でまとめたデータブロ
ックである。
タブロック)を対象として演算が行われる。仮に、12
8サンプル毎に解析処理を行うと、16秒(=128*
1/8)毎に脈拍数が更新される。この場合、更新間隔
を短くするためにFFT処理の対象となるサンプル数を
減らすと、必要な周波数帯域での解析結果を得ることが
できない。そこで、この例にあっては、1回の演算毎に
対象となるサンプルを1データブロック単位でずらし、
直前の128サンプルに対してFFT処理を行うように
している。これにより、データブロック(32サンプ
ル)毎に解析結果を出力して、脈拍数の更新を4秒毎に
行うことが可能となる。このように、現在の演算対象と
するデータを前回の演算対象となったデータと一部重複
させ、順次、演算対象となるデータをずらしながら行う
FFT処理を、以下、シフトFFT処理と称する。
脈波信号用FFT回路207に第1のデータブロックD
B1が取り込まれた後に行われる。この場合、脈波デー
タMDとしては計測開始時刻t0の直後にサンプルされ
た32サンプルしかないので、データブロックDB1に
0データを3データブロックDB(96サンプル)付加
したものを演算対象データとしている。なお、0データ
は、脈波データMDの中心値であるので、0データを付
加することによって、脈波信号MSに含まれていない周
波数スペクトルが演算結果に表れるといったことはな
い。
データを1データブロック分ずらして行われる。この場
合、データブロックDB1,DB2に0データを2データ
ブロックDB(96サンプル)付加したものを演算対象
データとする。この後、演算の対象となるデータを順次
ずらして、FFT処理を継続する。これにより、脈波の
計測開始から16秒経過して脈波データMDを128サ
ンプル得なくとも、4秒経過した時点で周波数解析を行
うことができるので、脈波計測開始から脈拍数の表示ま
での時間を大幅に短縮することが可能となる。
08は、体動データTDに対して、FFT処理を施し、
体動信号TSの各周波数成分を指示する体動解析データ
TKDを生成する。この場合も、上述した脈波データM
Dと同様に、シフトFFT処理が行われる。すなわち、
体動信号用FFT回路208は、まず、直前の128サ
ンプルに対してFFT処理を行ない、次の処理では32
サンプルずらした128サンプルに対してFFT処理を
行う。
KDのみから脈拍に相当する周波数を抽出し、その結果
を脈波数算出手段210に出力する(第1の抽出方
法)。また、脈波抽出手段209は、脈波解析データM
KDと体動解析データTKDとを比較して、脈波解析デ
ータMKDの示す各周波数成分から体動解析データTK
Dの示す各周波数成分を差し引いて、その結果に基づい
て脈拍に相当する周波数を抽出し、これを脈拍数算出手
段210に出力することも可能である(第2の抽出方
法)。この抽出方法の切換は、制御手段206から供給
される抽出切換信号TKSによって制御される。
に基づいて、体動の有無、すなわち利用者が運動状態に
あるのか安静状態にあるのかを1データブロック(32
サンプル)毎に判別する。具体的には、体動データTD
の正のピーク値から負のピーク値を減算して振幅値PPを
求め、これを予め設定された基準値REFと比較して、振
幅値PPが基準値REFを1データブロック期間連続して上
回る場合に体動有りと判定し、一方、振幅値PPが基準値
REFを1データブロック期間中のいずれかのサンプルで
下回る場合には体動無しと判定している。この場合、基
準値REFは、体動の有無を判別できるように実験等によ
って予め定められるが、この例にあっては、「5」に設
定している。そして、この判定結果に基づいて、脈波抽
出方法を指示する抽出切換信号TKS、バッファ203
から出力するデータを制御するバッファ制御信号BC
S、およびピッチ表示の有無を制御するピッチ表示制御
信号PCSが各々生成される。
抽出手段209に対して、体動無しと判定された場合に
は第1の抽出方法で脈拍に相当する周波数を算出するよ
う指示し、一方、体動有りと判定した場合には第2の抽
出方法で脈拍に相当する周波数を算出するよう指示す
る。また、バッファ制御信号BCSは、バッファ203
に対して、体動の状態が有りから無しに変化した後、一
定の期間、脈波データMDの一部を0データに置換して
出力するように指示する。また、ピッチ表示制御信号P
CSは、ピッチ検出手段211に対して、体動が有る場
合にピッチの値を示すピッチデータを出力し、体動が無
い場合にピッチデータを出力しないように指示する。
抽出手段209によって抽出された周波数に基づいて脈
拍数を算出し、これを示す脈拍数データを出力する。こ
の場合、抽出される周波数fは、脈波信号MSの基本波
成分に対応するので、脈拍数算出手段210は「60/
f」を演算によって求め、脈拍数を特定する。
データTKDに基づいて、ランニング中のピッチを指示
するピッチデータを算出する。この場合、体動解析デー
タTKDは、32サンプル(4秒)毎に生成されるの
で、ピッチ算出手段211によって算出されるピッチデ
ータも4秒毎に生成される。但し、ピッチデータの表示
部212への出力は、ピッチ表示制御信号PCSによっ
て制御されるため、制御手段206によって体動無しと
判定されると直ちにピッチデータの出力が停止され、ピ
ッチは表示されなくなる。これにより、ランニングを終
了した後に、ピッチが表示されるといったことがなく
る。
づいて脈拍数を液晶表示装置13に表示させるととも
に、ピッチデータに基づいてピッチを液晶表示装置13
に表示させる。これにより、利用者は、脈拍数やピッチ
といった生体の情報を認識することができる。
めのタイミングチャートである。図5(a)は体動デー
タTDの振幅値PPの変化についてその一例を示すグラフ
である。この例にあっては、時刻tにおいて、振幅値PP
が基準値REFを下回る。したがって、制御手段206
は、期間T1〜期間T5では体動有りと判定し、期間T
6〜期間T8では体動無しと判定する。なお、期間T5
で体動有りと判定されるのは、時刻ttに達した時点で
始めて制御手段206は体動が無かったと検知するから
である。また、図5(b)は、A/D変換器202から
出力される脈波データMDである。また、体動有りと判
定される期間中に検出されたデータブロックに
は、「’」を付してある。
号BCSは直前の128サンプルに相当するデータブロ
ックDB1'〜DB4'を出力するようにバッファ203を
制御する。これにより、図5(c)に示すデータがバッ
ファ203から出力される。また、期間T5において
は、体動有りと判定されるので、抽出切換信号TKS
は、第2の抽出方法で脈拍に対応した周波数を抽出する
ように脈波抽出手段209を制御する。この場合、脈波
成分抽出手段209は、データブロックDB1'〜DB4'
に基づいて生成された脈波解析データMKDと体動解析
データTKDとを比較し、この結果に基づいて脈拍に対
応する周波数を特定する。
ブロックは、図5(d)に示すDB2'、DB3'、DB4'
およびDB1になる。ここで、データブロックDB2'〜
DB4'には体動成分が重畳しており、データブロックD
B1には体動成分が重畳していない。したがって、デー
タブロックDB1に着目すれば、脈波成分抽出手段20
9は、脈波解析データMKDにのみ基づいて脈拍に対応
する周波数を抽出するのが望ましい(第1の抽出方
法)。しかし、FFT処理の対象となるデータブロック
DB2'〜DB4'は、体動有りに対応するものであるた
め、第1の抽出方法では、正確に脈拍に対応する周波数
を特定することができない。
御信号BCSは、データブロックDB2'〜DB4'を0デ
ータに置換して出力し、これに続いてデータブロックD
B1を出力するようにバッファ203を制御する。図5
(e)はバッファ203から出力されるデータである。
ここで、0データは脈波データMDの中心値であるか
ら、体動成分が重畳したデータブロックDB2'〜DB4'
は周波数成分の解析に影響を与えないデータに置換され
る。これにより、脈波信号用FFT回路207は、デー
タブロックDB1に基づいて周波数解析を行うことがで
きる。この際、抽出切換信号TKSは、第1の抽出方法
で脈拍に対応した周波数を抽出するように脈波抽出手段
209を制御する。この場合、脈波成分抽出手段209
は、脈波解析データMKDにのみ基づいて、脈拍に対応
する周波数を特定する。
り具体的に説明する。例えば、期間T2〜期間T5にお
ける脈波データMDの変化が図6に示すものである場
合、この脈波データMDにFFT処理を施すと、図7に
示す周波数スペクトラムが得られる。この例では、比較
的低域に体動成分が存在している。このように体動成分
が重畳した周波数解析結果から脈拍に対応する周波数を
特定することは、困難である。図8に、期間T2〜期間
T4において0データに置換した脈波データMDの変化
を示す。この脈波データMDは、図5(e)に示すもの
に相当し、これにFFT処理を施すと、図9に示す周波
数スペクトラムが得られる。この場合の周波数解析結果
には、体動成分が表れていないので、周波数解析結果か
ら脈拍に対応する周波数を特定することが可能となる。
すなわち、0データに置換することにより、実質的に体
動成分がないデータにのみ基づいて周波数解析が行われ
る。
Tの演算処理の対象となるデータブロックは、本来、図
5(f)に示すDB3'、DB4'、DB1、DB2である
が、バッファ制御信号BCSは、図5(g)に示すよう
にデータブロックDB3'、DB4'を0データに置換した
データを出力し、これに続いてデータブロックDB1、
DB2を出力するようにバッファ203を制御する。ま
た、この際、抽出切換信号TKSは、第1の抽出方法で
脈拍に対応した周波数を抽出するように脈波抽出手段2
09を制御する。また、期間T8において、バッファ制
御信号BCSは、図5(h)に示すデータブロックDB
4'、DB1、DB2、DB3のうち、データブロックDB
4'を0データに置換して、図5(i)に示すデータを出
力するようにバッファ203を制御する。このように、
判定結果が体動有りから体動無しに切り替わった直後の
3ブロック期間には、体動成分が重畳したブロックデー
タが存在するため、これらのデータを0データに置換す
る処理がバッファ203において行われる。この意味に
おいて、バッファ203は脈波データの一部を周波数解
析結果に影響を与えないデータへ置換する置換手段とし
て機能する。
号BCSは、図5(j)に示すデータブロックDB1〜
DB4を出力するようにバッファ203を制御する。デ
ータブロックDB1〜DB4は、体動成分が重畳していな
い脈波データMDであり、脈波成分抽出手段209は、
脈波解析データMKDにのみ基づいて脈拍に対応する周
波数を特定する。
ら体動無しに切り替わった後の所定期間において、体動
有りと判定された期間中に記憶した脈波データMDを周
波数解析に影響を与えないデータに置換して出力するの
で、当該期間において脈波信号用FFT回路207は、
実質的に、体動無しの期間に検出された脈波データMD
にのみ基づいて脈波の周波数解析を行うことができる。
この結果、例えば、ランニングを終えた直後のように生
体の状態が体動有りから体動無し変化した場合であって
も、正確な脈拍を直ちに表示させることができる。ま
た、体動信号TSにノイズが突発的に重畳して異常状態
に陥った場合であっても通常状態に戻ると、正確な脈拍
を直ちに表示させることができる。
理について説明する。図10は脈波抽出手段の動作を示
すフローチャートである。ステップS1において、脈波
抽出手段209は、抽出切換信号TKSが第2の抽出方
法を指示するか否かを判定する。第2の抽出方法を指示
する場合は、運動時であって、以下の手順に従って、脈
拍に対応する周波数を特定する。まず、ステップS2で
は、体動解析データTKDの最大周波数成分fmを特定
し、その1/2の周波数にある一定値TH以上の体動成
分があるかどうか判定する(ステップS2)。一定値T
H以上の体動成分がある場合、fmは第2高調波として
特定する(ステップS3)。一方、無い場合には、fm
の1/3の周波数に、ある一定値以上の体動成分がある
か否かを判定する(ステップS4)。ある一定値TH以
上の体動成分が存在するならば、fmを第3高調波とし
て特定する(ステップS5)。一方、無かった場合には
fmを基本波として特定する(ステップS6)。
何番目(変数HMCできていされている。)の高調波で
あるかを求め、ステップS7で基本波を求めるためにf
mを除する数値(変数HMC)を決定する。そこでステ
ップS7では、体動の基本波を求める。
おいては、脈波の周波数解析の結果の大きな線スペクト
ル順にその周波数と体動周波数との比較を行い、その周
波数が体動信号の基本波、第2高調波、第3高調波と一
致するかどうかを判定する。すなわち、脈波検出用セン
サユニット30によって検出された脈波信号MSを周波
数解析した結果と、体動検出用センサ装置90によって
検出された体動信号TSを周波数解析した結果との間
で、互いに重なる周波数があるかどうかを判定するもの
である。
数の基本波との比較を行い、ステップS10においては
体動周波数の第2高調波との比較を行い、ステップS1
1においては体動周波数の第3高調波との比較を行う。
これを、検出されたすべての脈波の周波数成分について
繰り返して行い、もし、仮に一致する周波数が存在する
場合は、この周波数成分を除去する。但し、脈波周波数
成分のうち最大のレベルを有する周波数成分のみを用い
て判定してもよい。これは、脈波の基本波のレベルが、
通常一番大きいからである。この処理を行うことで、ス
テップS12において、体動成分と一致しない最大脈波
成分fnを抽出することができる。
法を指示する場合には、ステップS1の判定結果はNO
となり、ステップS13に進んで、脈波解析データMK
Dのうち最大のものを脈波周波数成分fnとして特定す
る。この場合は、安静時の抽出方法に対応する。
サンプルを1データブロックとし、4データブロックに
対してFFT処理を行ったが(データブロックのサンプ
ル数は、高速フーリエ変換演算対象データに対して、1
/4)、64サンプルを1データブロックとしてもよ
い。この場合、FFT処理の演算の対象を128サンプ
ル(2データブロック)とし、1回の解析毎に1データ
ブロックずらすシフトFFT処理を行ってもよい(デー
タブロックのサンプル数は、高速フーリエ変換演算対象
データに対して、1/4)。
項によれば、生体の状態が体動有りから体動無しに切り
替わった場合に、体動有りの期間に検出された脈波デー
タについては周波数解析の結果に影響を与えないデータ
に置換するから、実質的に体動無しの期間に検出された
脈波データに基づいて脈拍数を特定することができ、こ
の結果、脈拍数表示の応答性を改善することができる。
また、体動無しの場合には、生体の体動ピッチを表示し
ないので、利用者が戸惑うこともない。
ク図である。
示す説明図である。
る。
である。
ためのタイミングチャートである。
一例を示す図である。
した周波数スペクトラムを示す図である。
データMDの変化を示す図である。
した周波数スペクトラムを示す図である。
すフローチャートである。
Claims (3)
- 【請求項1】 一定のサンプル数に対応するデータブロ
ック毎に脈拍を計測する脈拍計において、 生体の体動を検出して体動データを出力する体動検出手
段と、 前記体動データに基づいて体動の有無を前記データブロ
ック毎に判定する判定手段と、 前記生体の脈動を検出して脈波データを出力する脈波検
出手段と、 前記脈波データを記憶するとともに、前記判定手段の判
定結果が体動有りから体動無しに切り替わった後の所定
期間においては、体動有りと判定された期間中に記憶し
た前記脈波データを周波数解析に影響を与えないデータ
に置換して出力し、前記所定期間以外の期間においては
記憶した前記脈波データを出力する置換手段と、 前記置換手段の出力データに対して複数の前記データブ
ロック単位で周波数解析を施し、次の周波数解析を1デ
ータブロックずらした前記置換手段の出力データに対し
て施すことにより、前記データブロック毎に解析結果を
得る第1の演算手段と、 前記体動データに対して複数の前記データブロック単位
で周波数解析を施し、次の周波数解析を1データブロッ
クずらした前記体動データに対して施すことにより、前
記データブロック毎に解析結果を得る第2の演算手段
と、 前記判定手段によって体動有りと判定された場合には、
前記第1の演算手段の周波数解析に基づいて脈拍に対応
する周波数を抽出し、前記判定手段によって体動無しと
判定された場合には、前記第1の演算手段と前記第2の
演算手段の周波数解析結果に基づいて、脈拍に対応する
周波数を抽出する脈波成分抽出手段と、 前記脈波成分抽出手段によって抽出された前記脈拍の周
波数から脈拍数を算出する脈拍数算出手段と前記脈拍数
を表示する表示手段とを備えたことを特徴とする脈拍
計。 - 【請求項2】 前記第2の演算手段の周波数解析結果か
ら、前記生体の体動のピッチを算出するピッチ算出手段
と、 前記判定手段によって体動有りと判定された場合には前
記ピッチを表示するように前記表示手段を制御し、前記
判定手段によって体動無しと判定された場合には、前記
ピッチを表示しないように前記表示手段を制御する表示
制御手段とを備えたことを特徴とする請求項1に記載の
脈拍計。 - 【請求項3】 前記第1の演算手段および前記第2の演
算手段で行う周波数解析処理は高速フーリエ変換処理で
あり、前記データブロックのサンプル数は、高速フーリ
エ変換演算対象データに対して、1/4もしくは1/2
であることを特徴とする請求項1に記載の脈拍計。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP06348697A JP3564255B2 (ja) | 1997-03-17 | 1997-03-17 | 脈拍計 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP06348697A JP3564255B2 (ja) | 1997-03-17 | 1997-03-17 | 脈拍計 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH10258039A true JPH10258039A (ja) | 1998-09-29 |
| JP3564255B2 JP3564255B2 (ja) | 2004-09-08 |
Family
ID=13230634
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP06348697A Expired - Lifetime JP3564255B2 (ja) | 1997-03-17 | 1997-03-17 | 脈拍計 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3564255B2 (ja) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005323906A (ja) * | 2004-05-17 | 2005-11-24 | Seiko Instruments Inc | 生体情報計測装置及び生体情報計測方法 |
| US7195596B2 (en) | 2002-10-04 | 2007-03-27 | Seiko Instruments Inc. | Pulse wave detecting apparatus and fourier transform process apparatus |
| JP2015506746A (ja) * | 2012-01-16 | 2015-03-05 | ヴァレンセル, インコーポレイテッドValencell, Inc. | 慣性律動による生理的測定エラーの軽減 |
| WO2015129557A1 (ja) * | 2014-02-25 | 2015-09-03 | ローム株式会社 | 心拍数検出装置 |
| JP2016146933A (ja) * | 2015-02-12 | 2016-08-18 | ルネサスエレクトロニクス株式会社 | 脈拍計、周波数解析装置及び脈拍測定方法 |
| JP2017189288A (ja) * | 2016-04-12 | 2017-10-19 | ソニーモバイルコミュニケーションズ株式会社 | 検出装置、検出方法、およびプログラム |
| US9993204B2 (en) | 2013-01-09 | 2018-06-12 | Valencell, Inc. | Cadence detection based on inertial harmonics |
| US10390762B2 (en) | 2012-01-16 | 2019-08-27 | Valencell, Inc. | Physiological metric estimation rise and fall limiting |
| JP2019180823A (ja) * | 2018-04-10 | 2019-10-24 | セイコーエプソン株式会社 | 生体解析装置および生体解析方法 |
| JP2019209041A (ja) * | 2018-06-08 | 2019-12-12 | 株式会社デンソー | 脈波形解析方法、及びプログラム |
-
1997
- 1997-03-17 JP JP06348697A patent/JP3564255B2/ja not_active Expired - Lifetime
Cited By (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7195596B2 (en) | 2002-10-04 | 2007-03-27 | Seiko Instruments Inc. | Pulse wave detecting apparatus and fourier transform process apparatus |
| JP2005323906A (ja) * | 2004-05-17 | 2005-11-24 | Seiko Instruments Inc | 生体情報計測装置及び生体情報計測方法 |
| US10542896B2 (en) | 2012-01-16 | 2020-01-28 | Valencell, Inc. | Reduction of physiological metric error due to inertial cadence |
| JP2015506746A (ja) * | 2012-01-16 | 2015-03-05 | ヴァレンセル, インコーポレイテッドValencell, Inc. | 慣性律動による生理的測定エラーの軽減 |
| US11350884B2 (en) | 2012-01-16 | 2022-06-07 | Valencell, Inc. | Physiological metric estimation rise and fall limiting |
| US10631740B2 (en) | 2012-01-16 | 2020-04-28 | Valencell, Inc. | Reduction of physiological metric error due to inertial cadence |
| US10349844B2 (en) | 2012-01-16 | 2019-07-16 | Valencell, Inc. | Reduction of physiological metric error due to inertial cadence |
| US10390762B2 (en) | 2012-01-16 | 2019-08-27 | Valencell, Inc. | Physiological metric estimation rise and fall limiting |
| US11363987B2 (en) | 2013-01-09 | 2022-06-21 | Valencell, Inc. | Cadence detection based on inertial harmonics |
| US9993204B2 (en) | 2013-01-09 | 2018-06-12 | Valencell, Inc. | Cadence detection based on inertial harmonics |
| WO2015129557A1 (ja) * | 2014-02-25 | 2015-09-03 | ローム株式会社 | 心拍数検出装置 |
| JP2016146933A (ja) * | 2015-02-12 | 2016-08-18 | ルネサスエレクトロニクス株式会社 | 脈拍計、周波数解析装置及び脈拍測定方法 |
| US10568581B2 (en) | 2015-02-12 | 2020-02-25 | Renesas Electronics Corporation | Pulsimeter, frequency analysis device, and pulse measurement method |
| KR20160099489A (ko) | 2015-02-12 | 2016-08-22 | 르네사스 일렉트로닉스 가부시키가이샤 | 맥박계, 주파수 해석 장치 및 맥박 측정 방법 |
| US20160235371A1 (en) * | 2015-02-12 | 2016-08-18 | Renesas Electronics Corporation | Pulsimeter, frequency analysis device, and pulse measurement method |
| JP2017189288A (ja) * | 2016-04-12 | 2017-10-19 | ソニーモバイルコミュニケーションズ株式会社 | 検出装置、検出方法、およびプログラム |
| US10849570B2 (en) | 2016-04-12 | 2020-12-01 | Sony Corporation | Detection apparatus, detection method, and program |
| JP2019180823A (ja) * | 2018-04-10 | 2019-10-24 | セイコーエプソン株式会社 | 生体解析装置および生体解析方法 |
| JP2019209041A (ja) * | 2018-06-08 | 2019-12-12 | 株式会社デンソー | 脈波形解析方法、及びプログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3564255B2 (ja) | 2004-09-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3605216B2 (ja) | 脈拍計 | |
| CN102008299B (zh) | 利用具有两种波长的光波测量脉搏的方法和设备 | |
| US4869262A (en) | Device for displaying blood pressure | |
| WO1998041142A1 (en) | Pulsometer and pulsation indication method | |
| US7144375B2 (en) | Pulsimeter, control method for pulsimeter, wristwatch information device, control program, and recording medium | |
| RU2675399C2 (ru) | Система мониторинга частоты сердечных сокращений | |
| JP3610148B2 (ja) | 周期・周波数計測装置 | |
| JPH10258038A (ja) | 脈波検出装置および脈拍計 | |
| JP3564255B2 (ja) | 脈拍計 | |
| JPH11276448A (ja) | 信号抽出装置および信号抽出方法 | |
| EP0181067B1 (en) | Device for displaying a blood pressure value | |
| JP5605253B2 (ja) | 脈波計測装置、およびプログラム | |
| JP5582051B2 (ja) | 脈波計測装置、およびプログラム | |
| JPH10179528A (ja) | 脈波解析装置 | |
| JPH10234685A (ja) | 計測装置 | |
| JP5012866B2 (ja) | 生体情報処理装置、および、生体情報処理方法 | |
| JP3454989B2 (ja) | 計測データに対する補正方法 | |
| EP2497423A1 (en) | Biometric apparatus | |
| JP3584140B2 (ja) | 計測装置 | |
| JP2004113821A (ja) | 運動処方支援装置 | |
| JPH1071137A (ja) | ストレス度表示装置及びストレス度表示方法 | |
| JP2011212383A (ja) | 生体情報処理装置 | |
| JP2012095795A (ja) | 脈波解析法 | |
| JP2009195589A (ja) | 生体情報処理装置及びその制御方法 | |
| JPH10258042A (ja) | 健康管理指標測定装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20040210 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040409 |
|
| 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: 20040518 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040607 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080611 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090611 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090611 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100611 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110611 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120611 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130611 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140611 Year of fee payment: 10 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| EXPY | Cancellation because of completion of term |