JPH10258038A - 脈波検出装置および脈拍計 - Google Patents
脈波検出装置および脈拍計Info
- Publication number
- JPH10258038A JPH10258038A JP9063487A JP6348797A JPH10258038A JP H10258038 A JPH10258038 A JP H10258038A JP 9063487 A JP9063487 A JP 9063487A JP 6348797 A JP6348797 A JP 6348797A JP H10258038 A JPH10258038 A JP H10258038A
- Authority
- JP
- Japan
- Prior art keywords
- pulse wave
- pulse
- wave signal
- signal
- filter
- 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
- 238000001514 detection method Methods 0.000 claims abstract description 95
- 238000004364 calculation method Methods 0.000 claims description 40
- 238000001914 filtration Methods 0.000 claims description 7
- 230000007423 decrease Effects 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 37
- 239000011295 pitch Substances 0.000 description 49
- 238000012545 processing Methods 0.000 description 45
- 238000001228 spectrum Methods 0.000 description 42
- 238000000034 method Methods 0.000 description 33
- 230000008569 process Effects 0.000 description 18
- 230000008859 change Effects 0.000 description 17
- 238000010586 diagram Methods 0.000 description 13
- 239000004973 liquid crystal related substance Substances 0.000 description 12
- 230000002238 attenuated effect Effects 0.000 description 11
- 238000000605 extraction Methods 0.000 description 6
- 238000003672 processing method Methods 0.000 description 6
- 210000000707 wrist Anatomy 0.000 description 6
- 230000003321 amplification Effects 0.000 description 5
- 238000003199 nucleic acid amplification method Methods 0.000 description 5
- 239000008280 blood Substances 0.000 description 4
- 210000004369 blood Anatomy 0.000 description 4
- 239000000284 extract Substances 0.000 description 3
- 102000001554 Hemoglobins Human genes 0.000 description 2
- 108010054147 Hemoglobins Proteins 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 239000011521 glass Substances 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 206010005746 Blood pressure fluctuation Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000006793 arrhythmia Effects 0.000 description 1
- 206010003119 arrhythmia Diseases 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000000295 emission spectrum Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- IJGRMHOSHXDMSA-UHFFFAOYSA-N nitrogen Substances N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 1
- 229910052698 phosphorus Inorganic materials 0.000 description 1
- 239000011574 phosphorus Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
- 230000000284 resting effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 210000002303 tibia Anatomy 0.000 description 1
- 210000000623 ulna Anatomy 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
- A61B5/721—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts using a separate sensor to detect motion or using motion information derived from signals other than the physiological signal to be measured
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, 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/024—Detecting, measuring or recording pulse rate or heart rate
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, 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/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02438—Detecting, measuring or recording pulse rate or heart rate with portable devices, e.g. worn by the patient
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, 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/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
- A61B5/02427—Details of sensor
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, 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/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/0245—Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
- A61B5/02455—Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals provided with high/low alarm devices
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Cardiology (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physiology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
波信号を得るとともに、当該脈波信号に基づいて高い精
度で脈拍数を求める。 【解決手段】 脈波検出用センサユニット30からの脈
波信号はバッファ503に一時格納される。インパルス
ノイズ検出手段505によって、バッファ503上の脈
波信号中にインパルスノイズが検出されると、第1のデ
ジタルフィルタ506の通過帯域は、前回の脈拍数に相
当する周波数を中心にした山型となり、バッファ503
から出力される脈波信号中のインパルスノイズが低減さ
れる。その後、第2のデジタルフィルタ507および第
3のデジタルフィルタ508によって、脈波信号中の全
体的なノイズおよび体動成分が低減された後に、周波数
解析部509によって周波数解析され、この解析結果か
ら脈拍数が算出される。
Description
脈波検出装置、および当該脈波検出装置を用いた脈拍計
に関する。
・告知する脈拍計が市販されている。この種の脈拍計
は、人体の被測定部位に近接配置された脈波検出用セン
サからの出力信号(脈波信号)を用いて脈拍数を算出す
る。脈拍数の算出方法としては、以下に述べる矩形波処
理による方法および周波数分析による方法が知られてい
る。
し、この矩形波の周期を測定して脈拍数(周期の逆数に
比例)を算出する。すなわち、脈波信号のレベルの変動
をそのまま時間領域で調べることによって脈拍数を算出
しており、少ない演算量かつ小規模な回路構成で脈拍数
を算出することができるという特徴がある。
その結果として得られるスペクトラムにおける最大レベ
ルのスペクトルを抽出し、当該スペクトルの周波数から
脈拍数を算出する。すなわち、脈波信号に対して周波数
領域でのレベル比較を行って脈拍数を算出している。な
お、周波数解析方法としては、通常、FFTが採用され
る。
波成分以外のノイズが重畳することがある。ここで重畳
するノイズは、定常ノイズとは限らないので、入力段に
アナログフィルタを設けただけでは、その影響を十分に
取り除くことができない。そこで、本出願人は、以下に
述べるノイズの影響を低減する処理〜を提案してい
る。
イズの総称であり、このインパルスノイズが重畳した脈
波信号の一例を図11に示す。なお、図11(a)には
時間領域における脈波信号が、図11(b)には当該脈
波信号にFFT処理を施して得られるスペクトラムが示
されている。図から明らかなように、インパルスノイズ
が重畳したために、図11(a)中の期間t1〜t2に
おいて脈波信号が異常に変化し、図11(b)において
は、脈波の基本波を表すスペクトルSPよりも高レベル
のスペクトルが存在している。前述したように、周波数
分析による方法では、最も高レベルのスペクトルに基づ
いて脈拍数を算出するので、図に示すようなインパルス
ノイズが重畳した脈波信号に対して周波数分析を行って
も正確な脈拍数を算出することはできない。
ノイズの発生原因となる事象の有無を監視し、その監視
結果に基づいて、インパルスノイズが重畳する虞がある
場合には、脈波信号中のインパルスノイズを含む区間
(例えば図11(a)中の期間t1〜t2)にダミー信
号を挿入した後に周波数分析を行う装置を提案している
(詳しくは、特願平7−273238号の特許願に添付
した明細書および図面参照)。この装置によれば、例え
ば、図12(a)に示すように、インパルスノイズが重
畳していた期間t1〜t2に値が0のダミー信号が挿入
されるので、図12(b)に示すように、FFT処理の
結果として得られるスペクトラムにおいて、脈波の基本
波を表すスペクトルSPが最も高レベルのスペクトルと
なる。なお、図12から明らかなように、上述した装置
は周波数分析を前提としている。
出値から大きくずれる可能性は低い。このことを利用し
て、前回の検出値に一定の係数を掛けて今回の検出値の
適正な範囲(ウィンドウ)を設定し、範囲外の検出値が
得られた場合にはこれをノイズによる異常値として除去
する処理がウィンドウ処理である。なお、脈拍数に対す
るウィンドウの追従性が悪いと、図13に示すように、
運動開始時(t1)などの脈拍数が急激に変化し得ると
きには追従できず、その結果、検出値が正常値であって
も異常値として除去されてしまうという現象が発生す
る。しかも、この現象は、ウィンドウが正しく補正され
るまで継続する。そこで、本出願人はウィンドウの追従
性を向上させる技術を提案している(詳しくは、特願平
8−24511号の特許願に添付した明細書および図面
参照)。
被測定部位に近接配置されるので、使用者が運動してい
る場合には、その体動成分が脈波信号に重畳してしま
う。ここで、このような体動成分が重畳した脈波信号に
FFT処理を施して得られるスペクトラムを図14
(a)に例示する。この図に示す例では、図中左側のス
ペクトル群が脈波成分、右側のスペクトル群が体動成分
であり、いずれの群のスペクトルも同等程度のレベルと
なっている。もちろん、図14(a)は一例にすぎず、
場合によっては、体動成分のスペクトル中に最も高レベ
ルのスペクトルが存在することも有り得る。したがっ
て、体動成分が重畳した脈波信号に対して周波数分析を
行っても正確な脈拍数を算出することはできない。
設け、図14(a)に示すスペクトラムから、当該体動
検出用センサの出力信号(体動信号)にFFT処理を施
して得たスペクトラム(図14(b))を減算し、図1
4(c)に示すような脈波成分のみからなるスペクトラ
ムを得た後に、最も高レベルのスペクトルを選択する装
置を提案している(詳しくは、特開平7−273238
号公報参照)。この装置によれば、図から明らかなよう
に、選択されるスペクトルは、脈波の基本波を表すスペ
クトルSPとなる。なお、スペクトラムの減算を行うこ
とから明らかなように、上述した装置は周波数分析を前
提としている。
ルスノイズ除去処理では、インパルスノイズの発生原因
となる事象、例えば、腕時計型に適用した場合にはバッ
クライトの点灯やアラームの鳴鐘などを予め設定してお
き、設定された事象発生時に脈波信号にダミー信号を挿
入している。しかしながら、一般的なインパルスノイズ
は自装置の内部状態に無関係に発生することもあり、そ
れらの全ての事象を検出することは極めて困難である。
一方、事象の発生に関わらずにダミー信号を挿入するよ
うにすれば、インパルスノイズを完全に除去することが
できるが、本来の脈波成分も全て除去されてしまう。す
なわち、の処理だけでは本来の脈波成分への影響を最
小に抑制しつつ全てのインパルスノイズを除去すること
は困難である。
ウ外の検出値は除去されるので、ウィンドウの変化が脈
拍の変化に追従できない場合(例えば不整脈が発生した
場合)には、正当な検出値が異常値として除去されてし
まう。
クトラムの減算を行うようにしているが、実際には脈波
検出用センサの出力中の体動成分(図14(a)中右
側)と、体動検出用センサの出力中の体動成分(図14
(b))とが完全には一致しないことがあり、このよう
な場合には、前者から後者を減算しても体動成分を完全
に除去することはできない。
もので、脈波波形からノイズ成分を適切に除去した脈波
信号を得ることができる脈波検出装置を提供することを
第1の目的とし、当該脈波検出装置を用いて脈拍数を高
い精度で求めることができる脈拍計を提供することを第
2の目的としている。
ために、請求項1記載の発明は、脈波を検出し脈波信号
を出力する脈波検出用センサと、前記脈波検出用センサ
から出力される脈波信号をフィルタリングして出力する
可変特性のフィルタと、前記フィルタでフィルタリング
された脈波信号に基づいて脈拍数を算出する脈拍数算出
手段と、前記脈拍数算出手段で算出された脈拍数に基づ
いて前記フィルタの特性を設定する特性設定手段とを具
備することを特徴としている。
波を検出し脈波信号を出力する脈波検出用センサと、前
記脈波検出用センサから出力される脈波信号をフィルタ
リングして出力する可変特性のフィルタと、体動を検出
し体動信号を出力する体動検出用センサと、前記体動検
出用センサから出力される体動信号に基づいて体動のピ
ッチを算出するピッチ算出手段と、前記ピッチ算出手段
で算出されたピッチに基づいて前記フィルタの特性を設
定する前記特性設定手段とを具備することを特徴として
いる。
請求項1記載のものにおいて、前記脈波センサから出力
される脈波信号を一時格納してから前記フィルタへ出力
するバッファと、前記バッファに一時格納された脈波信
号からインパルスノイズを検出するインパルスノイズ検
出手段とを具備し、前記特性設定手段は、前記インパル
スノイズ検出手段の検出結果をも加味して前記フィルタ
の特性を設定することを特徴としている。
求項1記載のものにおいて、前記フィルタは、脈波信号
の通過レベルが基準周波数から脈波の基本波の下限周波
数および上限周波数にかけて緩やかに低くなるよう構成
されており、前記特性設定手段は、前記基準周波数を設
定することで前記フィルタの特性を設定することを特徴
としている。
請求項1または3あるいは4記載のものにおいて、前記
特性設定手段は、前記フィルタに設定する特性を前記脈
拍数算出手段で算出された脈拍数の変動状況に応じて変
更することを特徴としている。また、請求項6記載の脈
波検出装置は、請求項2記載の脈波検出装置において、
前記特性設定手段は、前記フィルタに設定する特性を前
記ピッチ算出手段で算出されたピッチの変動状況に応じ
て変更することを特徴としている。
1、3、4または5記載の脈波検出装置を備え、前記脈
拍数算出手段で算出された脈拍数を告知する告知手段を
具備することを特徴としている。また、請求項8記載の
脈拍計は、請求項2または6記載の脈波検出装置を備
え、前記フィルタでフィルタリングされた脈波信号に基
づいて脈拍数を算出する脈拍数算出手段と、前記脈拍数
算出手段で算出された脈拍数を告知する告知手段とを具
備することを特徴としている。
の実施形態について説明する。なお、本実施形態による
脈拍計は、一般的なデジタル腕時計としての機能を備え
ており、時計モードと脈拍計モードとを切り替えて使用
される。
を有する装置本体10と、この装置本体10に接続され
るケーブル20と、このケーブル20の先端側に設けら
れた脈波検出用センサユニット(脈波検出用センサ)3
0とから大略構成されている。装置本体10には、腕時
計における12時方向から腕に巻きついてその6時方向
で固定されるリストバンド12が設けられ、このリスト
バンド12によって、装置本体10は、腕に着脱自在に
なっている。脈波検出用センサユニット30は、センサ
固定用バンド40によって遮光されながら人差し指の根
元に装着されている。脈波検出用センサユニット30を
指の根元に装着する理由としては、ケーブル20の短縮
化による運動時の障害除去と外気の影響の低減が挙げら
れる。
などを外した状態で示す平面図、図3は、本脈拍計を3
時の方向から眺めた側面図である。図2において、装置
本体10は、樹脂製の時計ケース11(本体ケース)を
備えており、この時計ケース11の表面側には、現在時
刻や日付に加えて、脈拍数などの脈波情報などを表示す
る液晶表示装置13(表示装置)が構成されている。こ
の液晶表示装置13には、表示面の左上側に位置する第
1のセグメント表示領域131、右上側に位置する第2
のセグメント表示領域132、右下側に位置する第3の
セグメント表示領域133、及び左下側に位置するドッ
ト表示領域134が設けられており、ドット表示領域1
34には各種の情報をグラフィック表示することができ
る。
ンサユニット30による検出結果(脈波信号)に基づい
て脈拍数の変化などを液晶表示装置13で表示するため
に、この表示装置に対する制御や検出信号に対する信号
処理などを行う制御部5が設けられている。制御部5に
は、計時回路も構成されているため、通常時刻、ラップ
タイム、スプリットタイムなども液晶表示装置13に表
示可能である。また、時計ケース11の外周部及び表面
部には、時刻合わせ、モードの切換、ラップタイムや脈
波情報の計測を開始する旨の外部操作などを行うための
ボタンスイッチ111〜117が設けられている。
ース11に内蔵されているボタン形の電池59であり、
ケーブル20は、電池59から脈波検出用センサユニッ
ト30に電力を供給するとともに、脈波検出用センサユ
ニット30の検出結果を時計ケース11の制御部5に入
力している。なお、脈拍計1では、手首の屈曲の自由度
の確保と転倒時における手の甲の保護のために、装置本
体10の大型化の方向を3時及び9時の方向とし、さら
に腕時計における9時の方向に大きな張出部分101を
有するよう構成し、加えてリストバンド12を3時の方
向側に偏った位置で接続するようにしている。
池59に対して9時の方向には、ブザー用(報知音発生
用)の扁平な圧電素子58が配置されている。電池59
は、圧電素子58に比較して重いため、装置本体10の
重心位置は、3時の方向に偏るが、この偏った側にリス
トバンド12を接続して装着状態を安定化している。さ
らに、電池59と圧電素子58とを面方向に配置したこ
とによって本体10の薄型化を、裏面部119に電池蓋
118を設けたことによって電池59の交換の容易化を
実現している。
ストバンド12の端部に取り付けられた止め軸121を
保持するための連結部105が形成されている。時計ケ
ース11の6時の方向には、腕に巻かれたリストバンド
12か長さ方向の途中位置で折り返されるとともに、こ
の途中位置を保持するための留め具122が取り付けら
れる受け部106が形成されている。
部119から受け部106に至る部分は、時計ケース1
1と一体に成形されて裏面部119に対して約115゜
の角度をなす回転止め部108になっている。すなわ
ち、リストバンド12によって装置本体10を右の手首
L(腕)の上面部L1(手の甲の側)に位置するように
装着したとき、時計ケース11の裏面部119は、手首
Lの上面部L1に密着する一方、回転止め部108は、
榛骨Rのある側面部L2に当接する。この状態で、本体
10の裏面部119は、榛骨Rと尺骨Uを跨ぐ感じにあ
る一方、回転止め部108と裏面部119との屈曲部分
109から回転止め部108にかけては、榛骨Rに当接
する感じになる。このように、回転止め部108と裏面
部119とは、約115゜という解剖学的に理想的な角
度をなしているため、本体10を矢印Aの方向に、ま
た、装置本体10を矢印Bの方向に回そうとしても、本
体10はそれ以上不必要にずれない。また、裏面部11
9及び回転止め部108によって腕の回りの片側2カ所
で装置本体10の回転を規制するだけであるため、腕が
細くても、裏面部119及び回転止め部108は確実に
腕に接するので、回転止め効果が確実に得られる一方、
腕が太くても窮屈な感じがない。
の断面図である。図4において、脈波検出用センサユニ
ット30は、そのケース体としてのセンサ枠36の裏側
に蓋302が被されることによって、内側に部品収納空
間300が構成されている。部品収納空間300の内部
には、回路基板35が配置されている。回路基板35に
は、LED31、フォトトランジスタ32、その他の電
子部品が実装されている。脈波検出用センサユニット3
0には、ブッシユ393によってケーブル20の端部が
固定され、ケーブル20の各配線は、各回路基板35の
パターン上にはんだ付けされている。ここで、脈波検出
用センサユニット30は、ケーブル20が指の根元割か
ら装置本体10の側に引き出されるようにして指に取り
付けられる。従って、LED31及びフォトトランジス
タ32は、指の長さ方向に沿って配列されることにな
り、そのうち、LED31は指の先端側に位置し、フォ
トトランジスタ32は指の根元の方に位置するので、外
光がフォトトランジスタ32に届き難くなっている。
サ枠36の上面部分(実質的な脈波信号検出部)にガラ
ス板からなる透光板34によって光透過窓が形成され、
この透光板34に対して、LED31及びフォトトラン
ジスタ32は、それぞれ発光面及び受光面を透光板34
の方に向けている。このため、透光板34の外側表面3
41(指表面との接触面/センサ面)に指表面を密着さ
せると、LED31は、指表面の側に向けて光を発する
ことになり、フォトトランジスタ32は、LED31が
発した光のうち指の側から反射してくる光を受光するこ
とになる。ここで、透光板34の外側表面341と指表
面との密着性を高めるために、透光板34の外側表面3
41は、その周囲部分361から突出するよう構成され
ている。
GaN系(インジウム−ガリウム−窒素系)の青色LE
Dを用いており、その発光スペクトルは、450nmに
発光ピークを有し、その発光波長領域は、350nmか
ら600nmまでの範囲にある。かかる発光特性を有す
るLED31に対応させて、本実施形態では、フォトト
ランジスタ32として、GaAsP系(ガリウム−砒素
−リン系)のフォトトランジスタを用いてあり、その素
子自身の受光波長領域は、主要感度領域が300nmか
ら600nmまでの範囲にあって、300nm以下にも
感度領域がある。
ット30を、センサ固定用バンド40によって指の根元
に装着し、この状態で、LED31から指に向けて光を
照射すると、この光が血管に届いて血液中のヘモグロビ
ンによって光の一部が吸収され、一部が反射する。指
(血管)から反射してきた光は、フォトトランジスタ3
2によって受光され、その受光圧変化が血圧変化(血液
の脈波)に対応する。すなわち、血量が多いときには、
反射光が弱くなる一方、血圧が少なくなると、反射光が
強くなるので、反射光強度の変化を検出すれば、脈拍数
などを計測することができる。
mから600nmまでの範囲にあるLED31と、受光
波長領域が300nmから600nmまでの範囲のフォ
トトランジスタ32とを用いてあり、その重なり領域で
ある約300nmから約600nmまでの波長領域、す
なわち、約700nm以下の波長領域における検出結果
に基づいて脈波情報を表示する。かかる脈波検出用セン
サユニット30を用いれば、外光が指の露出部分にあた
っても、外光に含まれる光のうち波長領域が700nm
以下の光は、指を導光体としてフォトトランジスタ32
(受光部)にまで到達しない。外光に含まれる波長領域
が700nm以下の光は、指を透過しにくい傾向にある
ため、外光かセンサ固定用バンド40で覆われていない
指の部分に照射されても、点線Xで示すように、指を通
ってフォトトランジスタ32まで届かないからである。
また、血液中のヘモグロビンは、波長が300nmから
700nmまでの光に対する吸光係数が大きいので、約
700nm以下の波長領域の光を利用することにより、
S/N比の高い脈波信号を得ることができる。
ニットとの接続構造 図1及び図3に示したように、本体10の6時の方向に
おいて、回転止め部108として延設されている部分の
表面側には、コネクタ部70が構成され、そこには、ケ
ーブル20の端部に構成されたコネクタピース80を着
脱できるようになっている。従って、コネクタピース8
0をコネクタ部70から外せば、脈波計1を通常の腕時
計やストップウオッチとして用いることができる(この
際には所定のコネクタカバーを装着してコネクタ部70
を保護する)。
が、時計機能に関する部分については周知であるので、
当該部分を除いた脈拍計機能に関する部分についてのみ
説明する。
いる制御部5の要部の構成を示すブロック図であり、こ
の図に示すように、脈拍計の本体内部には加速度センサ
などの体動検出用センサ101が設けられており、制御
部5は、当該体動検出用センサ101の検出結果(体動
信号)と脈波検出用センサ30の検出結果(脈波信号)
に基づいて、脈拍数と運動のピッチとを求めるための構
成を備えている。
ンサユニット30の検出結果(脈波信号)を増幅して出
力する脈波信号増幅回路、502は脈波信号増幅回路5
01の出力(アナログ電圧信号)を所定ビットのデジタ
ル信号(例えば、−127〜127)に変換して出力す
るA/D変換器、503はA/D変換器502の出力を
一時格納して出力するバッファである。なお、バッファ
503の容量は、後段の周波数分析での検出時間(例え
ば16秒)に対応して適宜設定されるものであり、ここ
では、4秒相当の検出値を格納可能であるものとする。
なお、バッファ503上の4秒相当の検出値は、後述す
るインパルスノイズの検出処理の終了後すぐに出力され
るので、この4秒が装置全体の遅延時間に影響を与える
ことはほとんどない。
り、バッファ503に格納されている検出値がインパル
スノイズによる影響を受けているか否かを判定し、判定
結果を出力する。上記判定の手法は任意であるが、ここ
では、バッファ503上の全検出値数に対する所定の範
囲外の検出値の割合が所定のしきい値を超過した場合に
“1”を、それ以外の場合に“0”を判定結果として出
力するものとする。
ィルタであり、バッファ503から出力される脈波信号
に対して、順に、インパルスノイズの除去処理、疑似ウ
ィンドウ処理、体動成分除去処理を施す。各デジタルフ
ィルタとしては、例えばFIRフィルタが好適に用いら
れる。なお、各デジタルフィルタの処理の内容について
は後に詳述する。また、509は第3のデジタルフィル
タ508から出力される脈波信号に対して周波数解析
(例えばFFT処理)を施し、その結果(スペクトラ
ム)を出力する周波数解析部であり、図示を略すが、デ
ジタルフィルタ508からの脈波信号を所定の検出時間
(例えば16秒)だけ記憶可能な記憶手段を備えてい
る。
検出結果(体動信号)を増幅して出力する体動信号増幅
回路、511は体動信号増幅回路510の出力(アナロ
グ電圧信号)を所定ビットのデジタル信号(例えば、−
127〜127)に変換して出力するA/D変換器、5
12はA/D変換器511の出力を一時格納して出力す
るバッファであり、バッファ503と同一の容量を有す
る。
周波数解析部であり、バッファ512から出力される体
動信号に対して周波数解析(例えばFFT処理)を施
し、その結果(スペクトラム)を出力する。なお、通
常、各周波数解析部509,513の記憶手段は同一の
メモリ(例えばRAM)上に設けられる。
波信号を矩形波に変換して出力する脈波矩形波処理手
段、515は体動信号増幅回路510からの体動信号の
振幅が所定値(例えば50[mV])を超過した場合に
体動検出信号を出力する体動信号検出手段、516は各
周波数解析部509,513の出力結果から脈拍に対応
する周波数および体動のピッチに対応する周波数を抽出
して出力する脈波・体動成分抽出手段、517は脈波矩
形波処理手段514からの矩形波信号と、脈波・体動成
分抽出手段516からの各周波数とを入力し、脈拍数お
よびピッチの算出のための信号を出力する演算方法切替
手段であり、体動信号検出手段515からの体動検出信
号に基づいて、後段へ出力する信号を切り替える。な
お、体動のピッチに対する周波数(体動の基本波の周波
数)を求める処理については、動作の説明において詳述
する。
手段515から体動無し信号を受け取ると、脈波矩形波
処理手段514からの矩形波信号と脈波・体動成分抽出
手段516で抽出された体動のピッチに対応した周波数
とを後段へ出力するとともに、A/D変換器502、バ
ッファ503、インパルスノイズ検出手段505、第1
〜第3のデジタルフィルタ506〜508、周波数解析
部509の動作を停止させる。一方、体動信号検出手段
515から体動有り信号を受け取ると、A/D変換器5
02、バッファ503、インパルスノイズ検出手段50
5、第1〜第3のデジタルフィルタ506〜508、周
波数解析部509の動作を開始させ、演算方法切替手段
517は、以後、脈波・体動成分抽出手段516からの
各周波数を後段へ出力する。
介して供給される脈波の基本波の周波数または矩形波信
号に基づいて脈拍数を算出する脈拍数算出手段、519
は演算方法切替手段517を介して供給される体動の基
本波の周波数に基づいて体動のピッチを算出するピッチ
算出手段であり、それぞれ、算出した脈拍数/ピッチを
液晶表示装置13へ供給する。なお、周波数からの脈拍
数/ピッチの算出は、周波数に定数“60”を乗算する
ことで実現され、矩形波信号からの脈拍数/ピッチの算
出は、矩形波の周期を測定し、測定値の逆数に定数“6
0”を乗算することで実現される。
算出結果(脈拍数)に基づいて第1のデジタルフィルタ
506に設定する係数を算出して出力する第1の係数算
出手段(特性設定手段)であり、第1のデジタルフィル
タ506の特性は、インパルスノイズ検出手段505の
判定結果が“1”の場合、すなわちインパルスノイズが
検出された場合にのみ第1の係数算出手段520で算出
された係数で表される特性となり、その他の場合にはス
ルー(全領域通過)となる。なお、第1の係数算出手段
520が出力する係数には、基準となる周波数(基準周
波数)fb1と遮断周波数fl,fhとが含まれてお
り、基準周波数fb1は直前に検出された脈拍数に基づ
いて算出され、遮断周波数fl,fhは、脈拍数の変動
に応じて設定される。例えば、脈拍数の変動が大きく、
脈拍推移の安定度が低い場合には基準周波数fb1との
間隔が広くなるように、脈拍数の変動が小さく、脈拍推
移の安定度が高い場合には基準周波数fb1との間隔が
狭くなるように遮断周波数fl,fhが設定される。な
お、基準周波数fb1との間隔が一定となるように遮断
周波数fl,fhを設定する態様も考えられる。
する係数値で特定される第1のデジタルフィルタの特性
の一例を図6(a)に示す。この図に示すように、イン
パルスノイズ検出手段505においてインパルスノイズ
が検出された時に第1のデジタルフィルタ506に設定
される特性は、基準周波数fb1の通過レベルが最も高
い山型となる。すなわち、第1のデジタルフィルタ50
6に入力された信号は、周波数領域において、基準周波
数fb1近傍では減衰率が低く、基準周波数fb1から
遠ざかるにつれて減衰率が増大し、ある程度以上遠ざか
って遮断周波数fl,fhを超えると遮断される。
fb1が2.3[Hz]の場合、すなわち直前までの脈
拍数が140[拍/分]程度の場合には、インパルスノ
イズ検出時における第1のデジタルフィルタ506の出
力値が、5標本点の重み付け平均値 XN = (−1・XN−1・XN+1+4・XN+2−1・XN+
3−1・XN+4)/8(N=0,1,2,3,4,…,N
−1) となるよう設定されている。
果(脈拍数)に基づいて第2のデジタルフィルタ507
に設定する係数を算出して出力する第2の係数算出手段
(特性設定手段)であり、第2のデジタルフィルタ50
7の特性は、例外なく、第2の係数算出手段521で算
出された係数で表される特性となる。ここで、第2の係
数算出手段521が出力する係数値で特定される第2の
デジタルフィルタ507の特性の一例を図6(b)に示
す。この図に示すように、第2のデジタルフィルタ50
7に設定される特性は、基準周波数fb1の減衰率が最
も低く、裾野に向かって緩やかに減衰率が増大していく
山型となる。
の脈拍数を所定の数だけ記憶する記憶手段を備え、当該
記憶手段に記憶した脈拍数の変動に応じた係数を第2の
デジタルフィルタ507へ供給し、特性の稜線の傾きを
変化させる。例えば、脈拍数の変動が小さい場合には上
記稜線の傾きを急にし、同変動が大きい場合には上記稜
線の傾きを緩やかにする。このようにすることによっ
て、脈波信号中から脈波成分を大幅に低減してしまうと
いった事態を回避しつつ、脈波成分以外のノイズ成分を
低減することができる。
果(ピッチ)に基づいて第3のデジタルフィルタ508
に設定する係数を算出して出力する第3の係数算出手段
(特性設定手段)であり、第3のデジタルフィルタ50
8の特性は、例外なく、第3の係数算出手段522で算
出された係数で表される特性となる。ここで、第3の係
数算出手段522が出力する係数値で特定される第3の
デジタルフィルタ508の特性の一例を図6(c)に示
す。この図に示すように、第3のデジタルフィルタ50
8に設定される特性は、基準周波数fb2(体動の基本
波の周波数)の倍周波数においてくびれた特性、すなわ
ち基準周波数fb2の倍周波数の成分(体動成分)を遮
断する特性となる。
のピッチを所定の数だけ記憶する記憶手段を備え、当該
記憶手段に記憶したピッチの変動に応じた係数を第3の
デジタルフィルタ508へ供給し、くびれの強弱を変化
させる。例えば、ピッチの変動が小さい場合には、上記
くびれを強くし、基本周波数fb2の倍周波数の成分を
大幅に減衰させ、同変動が大きい場合には、上記くびれ
を弱くし、基本周波数fb2の倍周波数の成分であって
もあまり減衰しないようにする。このようにすることに
よって、脈波信号中から脈波成分を大幅に低減してしま
うといった事態を回避しつつ、脈波成分以外のノイズ成
分(体動成分)を大幅に低減することができる。
13は、脈拍数およびピッチを所定の領域に数値表示す
るが、数値に応じた絵文字や、脈拍数およびピッチの変
動の様子を表すグラフ等を表示するようにしてもよい。
また、第3のデジタルフィルタ508から出力される脈
波信号で表される波形(脈波波形)をスイープ表示する
ようにしてもよい。
構成要素は、実際にはCPU、ROM、RAMなどの回
路要素によってソフトウェア的に実現されているが、実
現手法は適宜設計事項であるので、いずれの構成要素が
どの回路要素によって実現されているか等の説明は省略
する。なお、省スペース化および低コスト化の観点か
ら、時計モードおよび脈拍計モードにおいて極力、上記
回路要素を兼用するよう設計することが望ましい。
計としての動作については周知であるので、その説明を
省略する。なお、所定のボタンを押下することにより、
時計モードと脈拍計モードとが切り替わる。
ャートであり、この図に示すように、脈拍計モードで
は、脈波検出用センサユニット30および体動検出用セ
ンサ101は、常時、脈波信号および体動信号を出力す
る。体動検出用センサ101から出力された体動信号
は、体動信号増幅回路510によって増幅され、体動信
号検出手段515へ供給される。その結果、脈拍数の演
算方法の切り替え基準となる体動検出信号が体動信号検
出手段515から演算方法切替手段517へ供給され
る。演算方法切替手段517では、体動信号検出手段5
15から供給された体動検出信号に基づいて脈拍数の演
算方法(周波数分析方法/矩形波処理方法)を切り替え
る(ステップS401,402)。
15から体動検出信号が供給されている場合、および体
動検出信号が供給されていない時間が所定時間未満であ
る場合には「体動有り」と判定し、脈波・体動成分抽出
手段516からの出力をそのまま後段の脈拍数算出手段
518およびピッチ算出手段519へ供給することによ
って、演算方法を周波数分析方法に切り替える。逆に、
上記体動検出信号が供給されていない時間が所定時間以
上となった場合には「体動無し」と判定し、脈波矩形波
処理手段514からの矩形波信号および脈波・体動成分
抽出手段516からの体動のピッチに対する周波数を後
段の脈拍数算出手段518およびピッチ算出手段519
へ供給することによって、演算方法を矩形波処理方法に
切り替える。なお、省電力の観点から、各方法において
不要な回路要素の動作あるいは同回路要素への給電は停
止される。
うに、脈波検出用センサユニット30によって検出され
た脈波信号は、脈波信号増幅回路501によって増幅さ
れ、脈波矩形波処理手段514によって矩形波信号に変
換され、演算方法切替手段517を介して脈拍数算出手
段518へ供給される。一方、体動検出用センサ101
によって検出された体動信号は、体動信号増幅回路51
0によって増幅され、A/D変換器511でデジタル信
号に変換され、バッファ512に一時格納される。バッ
ファ512から出力された体動信号は周波数解析部52
2で周波数解析処理(例えば検出時間が16秒のFFT
処理)を施され、当該周波数解析結果から脈波・体動成
分抽出手段516によって体動の基本波の周波数が抽出
される。
ラム)から、脈波・体動成分抽出手段516が脈波およ
び体動の基本波の周波数を求める処理について、図8を
参照して説明する。図8は体動信号の周波数解析結果の
一例を示す図であり、この図に示すように、一般には、
体動の基本波(腕の振りの基本波)の周波数成分より
も、体動の第2高調波の周波数成分の方が高レベルとな
る。体動の基本波は、腕の振り出し、引き戻しを一周期
とする振り子運動に相当し、通常走行において腕の振り
を滑らかな振り子運動にするのは難しいので、この成分
のレベルは低めとなるが、体動の第2高調波成分は、右
足をステップした時と左足をステップした時に均等に発
生する上下動や、腕の振り出し、引き戻しのそれぞれの
瞬間に加わる加速度に相当しているので、この成分のレ
ベルは高めとなるからである。
特徴的に得易いということになる。ところで、通常走行
では、2〜4[Hz]の範囲であれば、走行ペースの速
い遅いを考えても第2高調波が出現する領域をカバーす
ることができる。したがって、この領域に限定した上で
特徴的な第2高調波成分を抽出することで検出精度を上
げることができる。
の周波数分析結果から最も高レベルのスペクトルの周波
数(fs)を求め、次に、fsの1/2の周波数領域
に、所定のしきい値以上のレベルのスペクトルが存在す
るか否かを判定する。「存在する」と判定した場合に
は、fsを体動の第2高調波の周波数とし、fs/2を
体動の基本波の周波数として特定する。なお、「存在し
ない」と判定した場合には、fsを第3高調波の周波数
と仮定し、fs/3の周波数領域に所定のしきい値以上
のレベルのスペクトルが存在するか否かを判定し、存在
する場合にはfs/3を体動の基本波の周波数として特
定し、存在しない場合にはfsを体動の基本波の周波数
として特定する。なお、第3高調波まで考慮するのは、
体動の基本波の存在し得る範囲として、2〜4[Hz]
の範囲を想定しているためである。こうして特定された
体動の基本波の周波数は演算方法切替手段517を介し
てピッチ算出手段519へ供給される。以上が図7のス
テップS407の処理である。
分抽出手段516を介して供給された矩形波信号の周期
(波間値)を求め、この周期の逆数(周波数)に60を
乗じた値を脈拍数とし、ピッチ算出手段519では、脈
波・体動成分抽出手段516を介して供給された周波数
に60を乗じた値をピッチとする(ステップS40
8)。なお、各算出手段518は算出した脈拍数を係数
算出手段520,521へ、各算出手段519は算出し
たピッチを係数算出手段522へ供給するよう構成され
ているが、矩形波処理方法が選択されている間は、この
供給処理は停止される(もちろん、いかなる場合にも供
給処理を続行する態様としてもよい)。各算出手段51
8,519で算出された脈波数およびピッチは液晶表示
装置13へ供給されて表示され、使用者に視認される
(ステップS409)。もちろん、脈拍数やピッチを音
声で知らせるなど、視覚以外の感覚に訴えて告知する態
様としてもよい。
検出用センサユニット30によって検出された脈波信号
は、脈波信号増幅回路501によって増幅され、A/D
変換器502でデジタル信号(例えば−127〜127
の整数値をとる)に変換され(ステップS403)、バ
ッファ503に一時格納される。なお、体動信号につい
ては前述の「B−2:矩形波処理方法」における処理と
同一の処理が行われ、体動の基本波の周波数が得られる
ので、その説明を省略する。
波信号に対して、第1〜第3のデジタルフィルタ506
〜508によってフィルタ処理が施され、当該脈波信号
中のノイズ成分が低減または除去される(ステップS4
04)。以下、このフィルタ処理について順に説明す
る。バッファ503に格納された検出値(例えば4秒間
分の脈波信号)がインパルスノイズによる影響を受けて
いない場合、具体的には、バッファ503上の全検出値
数に対する所定の範囲外の検出値の割合が所定のしきい
値以下の場合には、インパルスノイズ検出手段505か
ら、インパルスノイズが検出されなかった旨の信号(例
えば、値が“0”の信号)が、バッファ503から第1
のデジタルフィルタ506へ上記検出値が出力されるタ
イミングで出力される。第1のデジタルフィルタ506
には、直前に算出された脈拍数に基づいて第1の係数算
出手段520によって算出された係数が供給されている
が、インパルスノイズ検出手段505からインパルスノ
イズが検出されなかった旨の信号が供給されているの
で、その特性は全領域スルーとなる。したがって、脈波
信号はそのまま第2のデジタルフィルタ507へ供給さ
れる。
がインパルスノイズによる影響を受けている場合、具体
的には、バッファ503上の全検出値数に対する所定の
範囲外の検出値の割合が所定のしきい値を超えている場
合には、インパルスノイズ検出手段505から、インパ
ルスノイズが検出された旨の信号(例えば、値が“1”
の信号)が、バッファ503から第1のデジタルフィル
タ506へ上記検出値が出力されるタイミングで出力さ
れる。この場合には、第1のデジタルフィルタ506の
特性は、図6(a)に示すような特性となり、予想され
る脈波の周波数成分以外の成分は減衰または遮断され
る。
よるフィルタ処理前の脈波信号の波形の例を図9(a)
に、そのFFT処理結果を図9(b)に、同脈波信号に
第1のデジタルフィルタ506によるフィルタ処理を施
した後の脈波信号の波形の例を図10(a)に、そのF
FT処理結果を図10(b)に示す。これらの図に示す
ように、第1のデジタルフィルタ506により、低周波
成分(インパルスノイズ)が大幅に減衰される。すなわ
ち、第1のデジタルフィルタ506によって、脈波信号
中のインパルスノイズ成分が大幅に低減または除去され
る。
に算出された脈拍数に基づいて第2の係数算出手段52
1によって算出された係数が供給されており、その特性
は図6(b)に示すようなものになっている。すなわ
ち、脈波の周波数から離れるにつれて減衰の程度が増す
ような特性となっている。したがって、インパルスノイ
ズ以外のノイズ成分もここで減衰される。ところで、安
静時のように脈拍数が急激には変化しない場合には、予
想される脈波の周波数から僅かでも離れた周波数もノイ
ズ成分と考えられる。
動するような場合には、予想される脈波の周波数からあ
る程度離れた周波数も本来の脈波成分を含んでいる可能
性がある。そこで、第2の係数算出手段521は、脈拍
数の変化が小さい場合には図6(b)の稜線の傾きを急
にし、脈拍数の変化が大きい場合には同傾きを緩やかに
して本来の脈波成分をあまり減衰させずに、ノイズ成分
を減衰させるようにしている。
に算出されたピッチに基づいて第3の係数算出手段52
2によって算出された係数が供給されており、その特性
は図6(c)に示すようなものになっている。前述した
ように、基準周波数fb2は体動の基本波の周波数であ
るので、第3のデジタルフィルタ508によって、脈波
信号中の体動の基本波成分および高調波成分が減衰され
る。なお、第2の係数算出手段521が第2のデジタル
フィルタ507に対して行うのと同様に、第3の係数算
出手段522は、ピッチの変化が小さい場合には図6
(c)のくびれを強くし、ピッチの変化が大きい場合に
は同くびれを弱くするので、ピッチが急激には変化しな
い場合には、体動成分が大幅に減衰・除去され、運動開
始時のようにピッチが急激に変動するような場合には、
前回のピッチと今回のピッチが大幅に異なる可能性があ
るので、体動成分の減衰は小幅に留められる。こうして
整形された脈波信号は、周波数解析部509によって所
定の周波数解析処理(本実施形態では検出時間が16秒
のFFT処理)を施される(ステップS405)。
ム)から、脈波・体動成分抽出手段516が脈波および
体動の基本波の周波数を特定する(ステップS40
6)。ここで、脈波の基本波の周波数を求める処理につ
いて説明する。脈波・体動成分抽出手段516は、ま
ず、脈波信号の周波数解析結果からレベルの高い順にス
ペクトルを選択し、逐次、選択したスペクトルの周波数
と体動の基本波の周波数(例えばfs/2)および高調
波の周波数(例えばfs,3fs/2)を比較し、一致
しない場合に当該スペクトルの周波数を脈波の基本波の
周波数として特定する。こうして特定された脈波の基本
波の周波数は演算方法切替手段517を介して脈拍数算
出手段518へ供給される。
タ509によって、脈波信号中の体動の基本波および高
調波の成分は大幅に減衰あるいは除去されているので、
単に、脈波信号の周波数解析結果において最も高レベル
のスペクトルの周波数を脈波の基本波の周波数として特
定するようにしてもよい。ただし、この場合、ピッチが
急激に変化するような状況下では、第3のデジタルフィ
ルタによる体動の基本波および高調波の成分の減衰量は
小幅にとどまるので、残留している体動成分のスペクト
ルを脈波の基本波のスペクトルとして選択してしまう可
能性がある。よって、ピッチの変化の程度に応じて処理
を変えるのが望ましい。
脈拍数算出手段518では、当該周波数から脈拍数が算
出され、この脈拍数が液晶表示装置13、第1の係数算
出手段520、および第2の係数算出手段521へ供給
される。また、体動の基本波の周波数が供給されたピッ
チ算出手段519では、当該周波数からピッチが算出さ
れ、このピッチが液晶表示装置13および第3の係数算
出手段522へ供給される(ステップS408)。な
お、脈拍数算出手段518から第1の係数算出手段52
0および第2の係数算出手段521へ供給する値、ピッ
チ算出手段519から第3の係数算出手段522へ供給
される値を、演算方法切替手段517を介して供給され
る各周波数とし、各係数算出手段520〜522を、当
該各周波数に基づいて係数を算出するように設計しても
よい。なお、液晶表示装置13に供給された脈拍数およ
びピッチは、対応する領域に表示され、使用者に告知さ
れる(ステップS409)。
ば、第1〜第3のデジタルフィルタ506〜508によ
って脈波信号中のノイズ成分を減衰または除去すること
ができるので、後段の周波数分析による脈拍数の検出処
理の精度を向上させることができる。また、脈波信号自
体を整形するので、例えば、脈波波形自体を表示するよ
うな態様にあっては、より正確な脈波波形を表示するこ
ともできる。
は、インパルスノイズが発生していないときには、全領
域スルーとなるので、インパルスノイズが重畳していな
い脈波信号をフィルタリングせずに済む。さらに、フィ
ルタリングによって主に減衰または除去されるのはイン
パルスノイズ成分であるので、インパルスノイズの重畳
部分における本来の脈波成分を大幅に減衰または除去し
てしまう虞がない。
は、除去される周波数成分は存在せず、あくまでも減衰
にとどめているので、脈拍数が突発的に変化しても、本
来の脈波成分を除去してしまう虞がない。さらに、第3
のデジタルフィルタによって体動の基本波および高調波
の周波数成分が大幅に減衰または除去されるので、後段
の処理を軽減することができる。また、第2および第3
のデジタルフィルタにあっては、脈拍数またはピッチの
変化の程度に応じて特性を変更することができるので、
より精度の高い脈波波形、脈拍数を得ることができる。
タを同時に使用しているが、いずれか1つまたは2つを
のみを使用する態様としてもよい。また、第2および第
3のデジタルフィルタを合成し、1つのフィルタとして
実現するようにしてもよい。さらに、バッファおよびイ
ンパルスノイズ検出手段を省略し、第1のデジタルフィ
ルタによるフィルタリングを全ての脈波信号に対して施
すような態様も考えられる。この場合には、第1〜第3
のデジタルフィルタを1つのフィルタとして実現可能で
ある。また、実装態様は、ネックレスや眼鏡など、腕時
計に限定されるものではない。もちろん、脈拍計単独の
装置としてもよいし、脈波を検出するための装置として
もよい。
ば、脈波検出用センサから出力される脈波信号をフィル
タリングして出力する可変特性のフィルタに、フィルタ
リング後の脈波信号に基づいて算出された脈拍数に基づ
いた特性、あるいは体動検出用センサから出力される体
動信号に基づいて算出された体動のピッチに基づいた特
性を設定することにより、脈波信号中のノイズ成分を低
減することができる。したがって、より高い精度で脈波
を検出することができる。
パルスノイズを検出した脈波信号に対して選択的に上記
フィルタリングを行うことにより、本来の脈波成分を大
幅に低減または除去することなく、インパルスノイズ成
分を低減または除去することができる。また、通過レベ
ルが基準周波数(前回の脈拍数に応じた周波数)から脈
波の基本波の下限周波数および上限周波数にかけて緩や
かに低くなるように上記フィルタの特性を設定すること
により、本来の脈波成分を大幅に低減または除去するこ
となく、ノイズ成分を全体的に低減することができる。
さらに、上記フィルタに設定する特性を脈拍数またはピ
ッチの変動状況に応じて変更することにより、より高い
精度で脈波を検出することができる。また、検出した脈
波に基づいた脈拍数を告知することにより、高い精度で
脈拍数を測定する脈拍計を実現することができる。
を示す図である。
ルなどを外した状態で示す平面図である。
る。
断面図である。
5の要部の構成を示すブロック図である。
を示す図であり、(a)は第1のデジタルフィルタの特
性、(b)は第2のデジタルフィルタの特性、(c)は
第3のデジタルフィルタの特性を示す。
ャートである。
求める処理を説明するための図である。
によるフィルタ処理前の脈波信号の波形の例を、(b)
は(a)のFFT処理結果を示す図である。
ルフィルタによってフィルタ処理した後の脈波信号の波
形の例を、(b)は(a)のFFT処理結果を示す図で
ある。
例を示す図であり、(a)は時間領域における脈波信号
を、(b)は(a)の脈波信号にFFT処理を施して得
られるスペクトラムを示す。
説明するための図であり、(a)は時間領域における脈
波信号を、(b)は(a)の脈波信号にFFT処理を施
して得られるスペクトラムを示す。
めの図である。
ための図であり、(a)は体動成分が重畳した脈波信号
にFFT処理を施して得られるスペクトラムを、(b)
は体動信号にFFT処理を施して得たスペクトラムを、
(c)は(a)から(b)を減算した結果のスペクトラ
ムを示す。
ト、101…体動検出用センサ、501…脈波信号増幅
回路、502,511…A/D変換器、503,512
…バッファ、505…インパルスノイズ検出手段、50
6…第1のデジタルフィルタ、507…第2のデジタル
フィルタ、508…第3のデジタルフィルタ、509,
513…周波数解析部、510…体動信号増幅回路、5
14…脈波矩形波処理手段、515…体動信号検出手
段、516…脈波・体動成分抽出手段、517…演算方
法切替手段、518…脈拍数算出手段、519…ピッチ
算出手段、520…第1の係数算出手段、521…第2
の係数算出手段、522…第3の係数算出手段
Claims (8)
- 【請求項1】 脈波を検出し脈波信号を出力する脈波検
出用センサと、 前記脈波検出用センサから出力される脈波信号をフィル
タリングして出力する可変特性のフィルタと、 前記フィルタでフィルタリングされた脈波信号に基づい
て脈拍数を算出する脈拍数算出手段と、 前記脈拍数算出手段で算出された脈拍数に基づいて前記
フィルタの特性を設定する特性設定手段とを具備するこ
とを特徴とする脈波検出装置。 - 【請求項2】 脈波を検出し脈波信号を出力する脈波検
出用センサと、 前記脈波検出用センサから出力される脈波信号をフィル
タリングして出力する可変特性のフィルタと、 体動を検出し体動信号を出力する体動検出用センサと、 前記体動検出用センサから出力される体動信号に基づい
て体動のピッチを算出するピッチ算出手段と、 前記ピッチ算出手段で算出されたピッチに基づいて前記
フィルタの特性を設定する前記特性設定手段とを具備す
ることを特徴とする脈波検出装置。 - 【請求項3】 前記脈波センサから出力される脈波信号
を一時格納してから前記フィルタへ出力するバッファ
と、 前記バッファに一時格納された脈波信号からインパルス
ノイズを検出するインパルスノイズ検出手段とを具備
し、 前記特性設定手段は、前記インパルスノイズ検出手段の
検出結果をも加味して前記フィルタの特性を設定するこ
とを特徴とする請求項1記載の脈波検出装置。 - 【請求項4】 前記フィルタは、脈波信号の通過レベル
が基準周波数から脈波の基本波の下限周波数および上限
周波数にかけて緩やかに低くなるよう構成されており、 前記特性設定手段は、前記基準周波数を設定することで
前記フィルタの特性を設定することを特徴とする請求項
1記載の脈波検出装置。 - 【請求項5】 前記特性設定手段は、前記フィルタに設
定する特性を前記脈拍数算出手段で算出された脈拍数の
変動状況に応じて変更することを特徴とする請求項1ま
たは3あるいは4記載の脈波検出装置。 - 【請求項6】 前記特性設定手段は、前記フィルタに設
定する特性を前記ピッチ算出手段で算出されたピッチの
変動状況に応じて変更することを特徴とする請求項2記
載の脈波検出装置。 - 【請求項7】 請求項1、3、4または5記載の脈波検
出装置を備え、 前記脈拍数算出手段で算出された脈拍数を告知する告知
手段を具備することを特徴とする脈拍計。 - 【請求項8】 請求項2または6記載の脈波検出装置を
備え、 前記フィルタでフィルタリングされた脈波信号に基づい
て脈拍数を算出する脈拍数算出手段と、 前記脈拍数算出手段で算出された脈拍数を告知する告知
手段とを具備することを特徴とする脈拍計。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP06348797A JP3584143B2 (ja) | 1997-03-17 | 1997-03-17 | 脈波検出装置および脈拍計 |
PCT/JP1998/001128 WO1998041143A1 (fr) | 1997-03-17 | 1998-03-17 | Detecteur d'onde d'impulsion et pulsometre |
EP98907274A EP0922433A4 (en) | 1997-03-17 | 1998-03-17 | PULSE WAVE DETECTOR AND PULSOMETER |
US09/180,341 US6155983A (en) | 1997-03-17 | 1998-03-17 | Pulse wave detecting device and pulse measurer |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP06348797A JP3584143B2 (ja) | 1997-03-17 | 1997-03-17 | 脈波検出装置および脈拍計 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH10258038A true JPH10258038A (ja) | 1998-09-29 |
JP3584143B2 JP3584143B2 (ja) | 2004-11-04 |
Family
ID=13230666
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP06348797A Expired - Lifetime JP3584143B2 (ja) | 1997-03-17 | 1997-03-17 | 脈波検出装置および脈拍計 |
Country Status (4)
Country | Link |
---|---|
US (1) | US6155983A (ja) |
EP (1) | EP0922433A4 (ja) |
JP (1) | JP3584143B2 (ja) |
WO (1) | WO1998041143A1 (ja) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006512990A (ja) * | 2003-01-13 | 2006-04-20 | ネルコアー ピューリタン ベネット インコーポレイテッド | 信号品質に基づく予め設定されたフィルタパラメータの選択 |
JP2006271731A (ja) * | 2005-03-29 | 2006-10-12 | Toshiba Corp | 心拍計測装置及び心拍計測方法 |
JP2007143623A (ja) * | 2005-11-24 | 2007-06-14 | Seiko Instruments Inc | 生体情報測定装置 |
JP2010004972A (ja) * | 2008-06-25 | 2010-01-14 | Fujitsu Ltd | 光電脈波計測装置および光電脈波計測用プログラム |
KR100963253B1 (ko) | 2008-03-10 | 2010-06-10 | (주)에이치쓰리시스템 | 맥박 신호에서 잡음을 제거하기 위한 방법 및 장치와기록매체 |
JP2012065707A (ja) * | 2010-09-21 | 2012-04-05 | Rohm Co Ltd | 測定システム |
JP2013059439A (ja) * | 2011-09-13 | 2013-04-04 | Nippon Koden Corp | 生体信号測定装置 |
JP2014054447A (ja) * | 2012-09-13 | 2014-03-27 | Omron Healthcare Co Ltd | 脈拍測定装置、ならびに、脈拍測定方法および脈拍測定プログラム |
US8795184B2 (en) | 2010-07-12 | 2014-08-05 | Rohm Co., Ltd. | Wireless plethysmogram sensor unit, a processing unit for plethysmogram and a plethysmogram system |
JP2015205187A (ja) * | 2015-06-09 | 2015-11-19 | セイコーエプソン株式会社 | 心房細動判定装置、心房細動判定方法およびプログラム |
JP2021027995A (ja) * | 2020-10-15 | 2021-02-25 | 日本電気株式会社 | 判定システム、判定方法、および判定プログラム |
US11660025B2 (en) | 2017-07-11 | 2023-05-30 | Nec Corporation | Identification device, identification method, and recording medium with recorded identification program |
Families Citing this family (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7279318B1 (en) * | 1999-06-09 | 2007-10-09 | Hybrid Systems Limited | Modification of biological elements |
ATE481922T1 (de) * | 2001-09-28 | 2010-10-15 | Csem Ct Suisse Electronique | Verfahren und vorrichtung zur pulsmessung |
US7341561B2 (en) * | 2003-05-30 | 2008-03-11 | Casio Computer Co., Ltd. | Wrist-worn high-accuracy pulsation measuring apparatus |
US20070070800A1 (en) * | 2005-09-27 | 2007-03-29 | Nathalie Virag | Externally worn vasovagal syncope detection device |
EP1908402B1 (fr) * | 2006-10-06 | 2015-12-09 | ETA SA Manufacture Horlogère Suisse | Méthode et dispositif intégré de mesure d'une pulsation cardiaque |
EP1908401A1 (fr) * | 2006-10-06 | 2008-04-09 | ETA SA Manufacture Horlogère Suisse | Méthode et dispositif de mesure d'une pulsation cardiaque lors de la pratique d'un sport rythmique |
US8157730B2 (en) | 2006-12-19 | 2012-04-17 | Valencell, Inc. | Physiological and environmental monitoring systems and methods |
US8652040B2 (en) | 2006-12-19 | 2014-02-18 | Valencell, Inc. | Telemetric apparatus for health and environmental monitoring |
EP2175773A2 (en) | 2007-08-07 | 2010-04-21 | Council of Scientific&Industrial Research | A non-invasive device nadi tarangini useful for quantitative detection of arterial nadi pulse waveform |
US8251903B2 (en) | 2007-10-25 | 2012-08-28 | Valencell, Inc. | Noninvasive physiological analysis using excitation-sensor modules and related devices and methods |
KR20110069084A (ko) | 2008-09-19 | 2011-06-22 | 센서즈 포 메드슨 앤드 사이언스 인코포레이티드 | 광학 센서 조립체 |
US8700111B2 (en) | 2009-02-25 | 2014-04-15 | Valencell, Inc. | Light-guiding devices and monitoring devices incorporating same |
US8788002B2 (en) | 2009-02-25 | 2014-07-22 | Valencell, Inc. | Light-guiding devices and monitoring devices incorporating same |
KR100997444B1 (ko) * | 2009-11-17 | 2010-11-30 | (주)에이치쓰리시스템 | 광용적맥파 측정기 |
US8888701B2 (en) | 2011-01-27 | 2014-11-18 | Valencell, Inc. | Apparatus and methods for monitoring physiological data during environmental interference |
JP5682369B2 (ja) * | 2011-02-23 | 2015-03-11 | セイコーエプソン株式会社 | 拍動検出装置 |
EP3222210B1 (en) | 2011-08-02 | 2024-09-25 | Yukka Magic LLC | Systems and methods for variable filter adjustment by heart rate metric feedback |
JP6003470B2 (ja) * | 2012-09-25 | 2016-10-05 | オムロンヘルスケア株式会社 | 血圧測定装置、脈波検出方法 |
GB2520487A (en) * | 2013-11-19 | 2015-05-27 | Eduardo Mangieri | Blood-flow sensor apparatus |
CN103767694B (zh) * | 2014-01-06 | 2015-07-08 | 西安交通大学 | 一种准确提取袖带压力震荡波的方法 |
WO2015131065A1 (en) | 2014-02-28 | 2015-09-03 | Valencell, Inc. | Method and apparatus for generating assessments using physical activity and biometric parameters |
JP6464004B2 (ja) | 2015-03-19 | 2019-02-06 | 株式会社東芝 | 測定装置、測定方法およびプログラム |
US10945618B2 (en) | 2015-10-23 | 2021-03-16 | Valencell, Inc. | Physiological monitoring devices and methods for noise reduction in physiological signals based on subject activity type |
WO2017070463A1 (en) | 2015-10-23 | 2017-04-27 | Valencell, Inc. | Physiological monitoring devices and methods that identify subject activity type |
WO2018009736A1 (en) | 2016-07-08 | 2018-01-11 | Valencell, Inc. | Motion-dependent averaging for physiological metric estimating systems and methods |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5825102U (ja) * | 1981-08-13 | 1983-02-17 | 松下電工株式会社 | 脈拍検出回路 |
JPH01288230A (ja) * | 1988-05-14 | 1989-11-20 | Matsushita Electric Works Ltd | 脈拍検出装置 |
DE4105457A1 (de) * | 1991-02-21 | 1992-09-03 | Seca Gmbh | Verfahren und vorrichtung zur nichtinvasiven kontinuierlichen blutdruckmessung bei menschen |
JP3422128B2 (ja) * | 1994-11-15 | 2003-06-30 | オムロン株式会社 | 血圧計測装置 |
JP3605216B2 (ja) * | 1995-02-20 | 2004-12-22 | セイコーエプソン株式会社 | 脈拍計 |
JP3075169B2 (ja) * | 1996-02-29 | 2000-08-07 | ヤマハ株式会社 | 半導体記憶装置 |
US5868679A (en) * | 1996-11-14 | 1999-02-09 | Colin Corporation | Blood-pressure monitor apparatus |
-
1997
- 1997-03-17 JP JP06348797A patent/JP3584143B2/ja not_active Expired - Lifetime
-
1998
- 1998-03-17 EP EP98907274A patent/EP0922433A4/en not_active Withdrawn
- 1998-03-17 WO PCT/JP1998/001128 patent/WO1998041143A1/ja not_active Application Discontinuation
- 1998-03-17 US US09/180,341 patent/US6155983A/en not_active Expired - Lifetime
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006512990A (ja) * | 2003-01-13 | 2006-04-20 | ネルコアー ピューリタン ベネット インコーポレイテッド | 信号品質に基づく予め設定されたフィルタパラメータの選択 |
JP2010227602A (ja) * | 2003-01-13 | 2010-10-14 | Nellcor Puritan Bennett Llc | 信号品質に基づく予め設定されたフィルタパラメータの選択 |
JP2006271731A (ja) * | 2005-03-29 | 2006-10-12 | Toshiba Corp | 心拍計測装置及び心拍計測方法 |
JP2007143623A (ja) * | 2005-11-24 | 2007-06-14 | Seiko Instruments Inc | 生体情報測定装置 |
KR100963253B1 (ko) | 2008-03-10 | 2010-06-10 | (주)에이치쓰리시스템 | 맥박 신호에서 잡음을 제거하기 위한 방법 및 장치와기록매체 |
JP2010004972A (ja) * | 2008-06-25 | 2010-01-14 | Fujitsu Ltd | 光電脈波計測装置および光電脈波計測用プログラム |
US8795184B2 (en) | 2010-07-12 | 2014-08-05 | Rohm Co., Ltd. | Wireless plethysmogram sensor unit, a processing unit for plethysmogram and a plethysmogram system |
JP2012065707A (ja) * | 2010-09-21 | 2012-04-05 | Rohm Co Ltd | 測定システム |
JP2013059439A (ja) * | 2011-09-13 | 2013-04-04 | Nippon Koden Corp | 生体信号測定装置 |
US9125605B2 (en) | 2011-09-13 | 2015-09-08 | Nihon Kohden Corporation | Biological signal measuring apparatus |
JP2014054447A (ja) * | 2012-09-13 | 2014-03-27 | Omron Healthcare Co Ltd | 脈拍測定装置、ならびに、脈拍測定方法および脈拍測定プログラム |
JP2015205187A (ja) * | 2015-06-09 | 2015-11-19 | セイコーエプソン株式会社 | 心房細動判定装置、心房細動判定方法およびプログラム |
US11660025B2 (en) | 2017-07-11 | 2023-05-30 | Nec Corporation | Identification device, identification method, and recording medium with recorded identification program |
JP2021027995A (ja) * | 2020-10-15 | 2021-02-25 | 日本電気株式会社 | 判定システム、判定方法、および判定プログラム |
Also Published As
Publication number | Publication date |
---|---|
JP3584143B2 (ja) | 2004-11-04 |
WO1998041143A1 (fr) | 1998-09-24 |
EP0922433A4 (en) | 2001-11-14 |
EP0922433A1 (en) | 1999-06-16 |
US6155983A (en) | 2000-12-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3584143B2 (ja) | 脈波検出装置および脈拍計 | |
EP0729726B1 (en) | Pulse rate meter | |
JP5369726B2 (ja) | 拍動検出装置、および拍動検出方法 | |
EP0805414B1 (en) | Pace measuring apparatus, electronic appliance and pace measuring method | |
JP3726832B2 (ja) | 脈拍計、腕時計型情報機器、制御プログラムおよび記録媒体 | |
JP2816944B2 (ja) | 脈拍計 | |
JPH10258040A (ja) | 脈拍計 | |
EP2859842A1 (en) | Reflection type photodetection apparatus, and biological information measuring apparatus | |
JP2004097371A (ja) | 脈拍計、脈拍計の制御方法、時計型情報機器、制御プログラムおよび記録媒体 | |
JP6115330B2 (ja) | 生体情報測定装置、および生体情報測定方法 | |
JP5605253B2 (ja) | 脈波計測装置、およびプログラム | |
JP5582051B2 (ja) | 脈波計測装置、およびプログラム | |
JP5012866B2 (ja) | 生体情報処理装置、および、生体情報処理方法 | |
JP3564255B2 (ja) | 脈拍計 | |
EP0803733B1 (en) | Measuring instrument, portable electronic equipment and measuring method | |
JP2005218595A (ja) | 精神ストレス評価方法および装置 | |
JP3584140B2 (ja) | 計測装置 | |
JPH1071137A (ja) | ストレス度表示装置及びストレス度表示方法 | |
JP2004113821A (ja) | 運動処方支援装置 | |
JPH09276254A (ja) | 作業疲労度観測システム | |
JP2011087838A (ja) | 拍動検出装置及び拍動検出方法 | |
JP5655721B2 (ja) | 生体情報処理装置及び生体情報処理方法 | |
JP2009207713A (ja) | 生体情報処理装置、生体情報処理方法および制御プログラム | |
JP3558428B2 (ja) | 生体情報計測装置 | |
JPH0824232A (ja) | リラックス度表示装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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: 20040706 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040802 |
|
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: 20080806 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090806 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100806 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100806 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110806 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120806 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130806 Year of fee payment: 9 |
|
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 |