JP4401014B2 - 脈波解析装置 - Google Patents
脈波解析装置 Download PDFInfo
- Publication number
- JP4401014B2 JP4401014B2 JP2000262600A JP2000262600A JP4401014B2 JP 4401014 B2 JP4401014 B2 JP 4401014B2 JP 2000262600 A JP2000262600 A JP 2000262600A JP 2000262600 A JP2000262600 A JP 2000262600A JP 4401014 B2 JP4401014 B2 JP 4401014B2
- Authority
- JP
- Japan
- Prior art keywords
- pulse wave
- subband
- wave data
- analysis
- wavelet transform
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Description
【発明の属する技術分野】
本発明は、被検者の脈波を解析する脈波解析装置に関する。
【0002】
【従来の技術】
近年、健康管理を目的として、日常生活や、ジョギング等の運動時における心臓の拍動数(心拍数)をモニターするニーズが高まっている。
心拍数を検出するには、胸部に電極を取り付け、心拍に伴って発生する活動電位を計測して心電図を描き、この心電図に描かれる波形のピーク間隔時間から算出する方法が一般に知られている。
【0003】
しかし、電極を胸部に張り付ける必要があり、手間がかかるので、計測した脈波から簡便に脈拍数を検出する方法が代用として行われている。
脈波は、心拍によって起きる動脈内の圧力変化が、末梢動脈に波動として伝わったものであり、例えば、光学式脈波センサによって検出することができる。
光学式脈波センサは、血液中のヘモグロビンの光吸収特性を利用して、末梢動脈の血液の波動的な流量変化を計測するものであり、人体の適所(指、腕、こめかみ等)に簡便に装着することができ、被験者に大きな負担をかけない。
【0004】
【発明が解決しようとする課題】
上記従来の脈拍数検出装置は、以下の課題を有する。
図4に示す様に、心拍数、脈拍数は、それぞれ、60を心電波形と脈波波形の振幅のピーク間隔時間(秒)で割って求められる。
通常、心電波形と脈波波形の振幅のピーク位置とは同期しており、心拍数と脈拍数とは一致する。
しかし、被験者は、日常生活や運動を行っているので、時々、脈波センサを装着した計測部位に体動が生じる。
体動が生じると、末梢動脈の血流が乱れ、心拍とは無関係な脈波の振幅のピークが発生するので、心拍数と脈拍数とが一致しなくなり、脈拍数を心拍数の代用として利用できなくなる。
【0005】
また、心拍と無関係な脈波の振幅のピークの周波数が、心拍に同期する脈波の振幅のピークの周波数に近いので、心拍と無関係な脈波をノイズ除去用のフィルタで除去できない。
【0006】
一方、特開平7- 299044号公報には、運動ノイズ検出センサを設け、脈波検出センサから得られる脈拍信号からピーク検出信号とその時間を求め、ピーク検出信号のうち運動ノイズによる信号を除去し、他のピーク検出信号から脈拍信号の周期を求め、該周期から脈拍数を求める脈拍検出装置が記載されている。しかし、運動ノイズ検出センサでは検知できない、血流の乱れに伴うノイズが脈拍信号に含まれている場合には有効に機能しない。
【0007】
本発明の目的は、被検者が体動しても、脈波を精度良く解析することができる脈波解析装置の提供にある。
【0008】
【課題を解決するための手段】
脈波計測手段は、被検者の適所に装着した脈波センサが出力する脈波を計測する。
ウェーブレット変換手段は、脈波計測手段が計測した脈波データを複数のサブバンドにウェーブレット変換し、脈波データが細かく周波数成分に分解される。
逆ウェーブレット変換手段は、選定されたサブバンドを逆ウェーブレット変換する。
解析手段は、再構成された脈波データを解析する。
脈波解析装置は、脈波が整形されてピークが明確であるので、被検者が体動しても、脈波を精度良く解析することができる。
【0009】
また、脈波解析装置は、サブバンドを選定するサブバンド選定手段とを備える。
サブバンド選定手段は、安静状態にある被検者から所定時間のあいだ計測した脈波データから解析指標X1を算出しておき(s11、s12)、さらに該脈波データを複数のサブバンドにウェーブレット変換し(s13)、複数のサブバンドから任意の1つのサブバンドを取り除いた残りのサブバンド群を逆ウェーブレット変換して再構成した脈波データから解析指標X2を求め(s15、s16)、解析指標X1と解析指標X2との差異が設定値以上になった場合に、取り除いたサブバンドを採用する(s17、s18)ことを繰り返すことによりサブバンドを選定することを特徴とする(図2参照)。
【0010】
なお、差異は、脈波の立ち上がりの最大傾斜点間の距離の差の統計値、または脈波のピーク間距離の差の統計値である。
【0011】
選定処理が終了すると、被検者の適所に装着した脈波センサが出力する脈波を所定時間、脈波計測手段が計測する(ステップS4)。
つぎに、この計測した脈波データをウェーブレット変換手段が複数のサブバンドにウェーブレット変換する(ステップS5)。ウェーブレット変換することにより、脈波データは細かく周波数成分に分解される。なお、その分解されたデータ係数の集合を一つのサブバンドと呼ぶ。
【0012】
逆ウェーブレット変換手段は、選定された複数のサブバンドを用いて逆ウェーブレット変換する(ステップS6)。これにより、脈波データが再構成される。再構成された脈波データを解析手段が各種解析処理を行う(ステップS7)。解析処理は、例えば、脈拍数の算出や脈拍間隔の揺らぎ解析等である。
【0013】
【発明の実施の形態】
本発明の一実施例を図2、図3、図5、および図6に基づいて説明する。図3に示す様に、脈波解析装置Aは、被検者に装着する脈波センサ1と、脈波センサ1を駆動するための駆動回路2と、脈波を計測して解析するデータ処理部3と、表示器(図示せず)とを有する。
【0014】
本実施例では、脈波センサ1は光学式であり、窓11を有するセンサユニット12内に、発光素子13と受光素子14とを配設したものであり、例えば、被検者の指に装着する。
【0015】
発光素子13から被検者の皮膚15に向かって光130が照射されて人体内部16に入り、毛細血管17に当たった光の一部が、毛細血管17中を流れる血液中のヘモグロビンに吸収される。また、残りの光が毛細血管17で反射して散乱し、その一部が受光素子14に入光する。
【0016】
この際、血液の流動により、毛細血管17中のヘモグロビンの量が波動的に変化するので、ヘモグロビンに吸収される光も波動的に変化する。
その結果、毛細血管17で反射して受光素子14で検出される受光量が変化し、その受光素子14の受光量が変化し、センサ出力100(脈波)の電圧変化に現れる。
【0017】
駆動回路2は、発光素子13に駆動用の電力を供給する回路であり、駆動用の電力が供給されると発光素子13が発光する。
データ処理部3は、センサ出力100(脈波)を増幅する検出回路31と、増幅された脈波101(アナログ信号)をデジタル信号102に変換するA/Dコンバータ32と、デジタル信号102を処理するマイクロコンピュータ33とを有する。
【0018】
マイクロコンピュータ33には、デジタル信号102を図2に示すフローチャートに基づいて処理するプログラム(脈波計測手段、ウェーブレット変換手段、逆ウェーブレット変換手段、解析手段、サブバンド選定手段に相当)が組み込まれている。
【0019】
つぎに、マイクロコンピュータ33の作動および脈拍解析装置Aの作用効果を図2に示すフローチャートに基づいて説明する。
マイクロコンピュータ33は、脈波の計測・解析に先立って後述するサブバンドの選定処理(ステップs11〜s19)をステップs1で行う。
【0020】
サブバンドの選定処理が終了すると、ステップs2でデジタル信号102を所定時間、取り込み脈波データを得る。
つぎに、ステップs3で、この脈波データをN個のサブバンドにウェーブレット変換する。ウェーブレット変換することにより、脈波データは細かく周波数成分に分解される。なお、その分解されたデータ係数の集合を一つのサブバンドと呼ぶ。
【0021】
後述するサブバンドの選定処理(ステップs11〜s19)により選定された複数のサブバンドを用いて、ステップs4で脈波データを逆ウェーブレット変換する。これにより、脈波データが再構成される。
ステップs5で、再構成された脈波データを用いて各種解析処理を行う。解析処理は、例えば、脈拍数の算出や脈拍間隔の揺らぎ解析等である。
【0022】
サブバンドの選定処理(ステップs11〜s19)について
ステップs11で、安静状態にある被検者から所定時間のあいだ脈波101をサンプリングしてデジタル信号102を所定時間、取り込み脈波データを得る。ステップs12で、この脈波データから解析指標X1を算出する。解析指標X1は、例えば、脈拍数や脈拍間隔である。
【0023】
ステップs13で、安静状態にある被検者から所定時間のあいだ脈波101をサンプリングしてデジタル信号102を所定時間、取り込み、得られた脈波データをN個のサブバンドにウェーブレット変換する。また、変数Iを0とする。ウェーブレット変換することにより、脈波データが細かく周波数成分に分解される。なお、分解されたデータ係数の集合を一つのサブバンドと呼ぶ。
【0024】
ステップs14で変数Iを1〜Nまで変化させる。
ステップs15で、I個目のサブバンドを除いた残りのサブバンドを逆ウェーブレット変換する。
【0025】
ステップs16で、逆ウェーブレット変換により再構成された所定時間の脈波データから解析指標X2を算出する。
ステップs17で、|X1−X2|≧設定値であるか否かを判別し、|X1−X2|≧設定値である場合(YES)にはステップs18に進み、|X1−X2|<設定値である場合(NO)にはステップs14に戻る。
ステップs18でI番目のサブバンドを選定する。
ステップs19でI≧Nであるか否かを判別し、I≧Nである場合(YES)にはステップs2に進み、I<Nである場合(NO)にはステップs14に戻る。
【0026】
I番目のサブバンドが解析指標に大きく影響を与えている場合には、|X1−X2|≧設定値になる。つまり、I番目のサブバンドには解析指標に必要な情報が含まれているといえる。
このようにして、サブバンドを選定すれば、解析しようとしている指標に係る情報(例えば、脈拍数の解析に必要な情報)のみを抽出できる。
【0027】
脈波解析装置Aは、サブバンドの選定処理(ステップs11〜s18)により選定された複数のサブバンドを用いて、脈波データを逆ウェーブレット変換して脈波データを再構成している。
このため、被検者の体動により血流に乱れが生じて、図5の(a)に示す脈波データ(原波形)が計測されても、図5の(b)に示す波形(再構成した波形)に整形されてピークが明確となるので脈波データを精度良く解析することができる。
【0028】
本発明は、上記実施例以外に、つぎの実施態様を含む。
a.図6の(a)は、安静状態にある被検者から脈波101を計測して、デジタル信号102に変換して得られた脈波データ(安静時脈波データの原波形)である。
図6の(b)は、その脈波データ(安静時脈波データの原波形)からI番目のサブバンドを除いた残りのサブバンドを逆ウェーブレット変換して再構成した脈波データ(安静時脈波データの原波形から、I番目のサブバンドを除いて再構成した脈波データ)である。
【0029】
或る解析指数に基づいて、この二つの波形を評価した場合、その差が大きければ、I番目のサブバンドは、その解析指標にとって重要であるので、I番目のサブバンドを選定する。
ここで、解析指標を図6に示すY1…、X1…とする。
Y1…は、安静時脈波データの立ち上がりの最大傾斜点間の距離(時間の単位)である。
X1…は、安静時脈波データの原波形から、I番目のサブバンドを除いて再構成した脈波データの立ち上がりの最大傾斜点間の距離(時間の単位)である。
【0030】
これら、Y1…、X1…に基づいてサブバンドの重要性を判断する。例えば、|Y1−X1|、|Y2−X2|、|Y3−X3|……の合計、平均値、二乗平均値、分散、標準偏差といった統計値である。
また、上記最大傾斜点を、脈波データのピーク間距離としても良い。
【0031】
b.脈波センサは、指以外に、被検者の耳、足、胴体、首等に装着しても良い。また、装着する脈波センサは、光学式以外に、超音波式、圧力式、ドップラー式等でも良い。
【0032】
c.上記実施例の脈波センサ1の発光素子13は、LEDが使用でき、受光素子14はフォトダイオードやフォトトランジスタが使用できる。また、受光素子14の、抵抗値変化、電流値変化、または電圧値変化に基づいて脈波を検出すれば良い。
【図面の簡単な説明】
【図1】請求項1を説明するためのフローチャートである。
【図2】本発明の一実施例に係る脈波解析装置の作動を示すフローチャートである。
【図3】本発明の一実施例に係る脈波解析装置の説明図である。
【図4】心電波形と脈波波形とを比較したグラフである。
【図5】(a)は被検者の体動により血流に乱れが生じている場合の脈波データ(原波形)であり、(b)はその脈波データを再構成した場合の波形(再構成した波形)である。
【図6】(a)は安静状態にある被検者から得られた脈波データ(安静時脈波データの原波形)であり、(b)はその脈波データからI番目のサブバンドを除いた残りのサブバンドを逆ウェーブレット変換して再構成した脈波データ(安静時脈波の原波形から、I番目のサブバンドを除いて再構成した脈波データ)である。
【符号の説明】
A 脈波解析装置
1 脈波センサ
33 マイクロコンピュータ(脈波計測手段、ウェーブレット変換手段、逆ウェーブレット変換手段、解析手段、サブバンド選定手段)
X1、X2 解析指標(差異)
Claims (3)
- 被検者の適所に装着した脈波センサが出力する脈波を計測する脈波計測手段と、
脈波データを複数のサブバンドにウェーブレット変換するウェーブレット変換手段と、
選定されたサブバンドを逆ウェーブレット変換する逆ウェーブレット変換手段と、
再構成された脈波データを解析する解析手段と、
サブバンドを選定するサブバンド選定手段とを備え、
該サブバンド選定手段は、
安静状態にある前記被検者から所定時間のあいだ計測した脈波データから解析指標X1を算出しておき、さらに該脈波データを複数のサブバンドにウェーブレット変換し、複数のサブバンドから任意の1つのサブバンドを取り除いた残りのサブバンド群を逆ウェーブレット変換して再構成した脈波データから解析指標X2を求め、
解析指標X1と解析指標X2との差異が設定値以上になった場合に、取り除いたサブバンドを採用することを繰り返すことによりサブバンドを選定することを特徴とする脈波解析装置。 - 前記差異は、前記脈波の立ち上がりの最大傾斜点間の距離の差の統計値であることを特徴とする請求項1記載の脈波解析装置。
- 前記差異は、脈波のピーク間距離の差の統計値であることを特徴とする請求項1記載の脈波解析装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000262600A JP4401014B2 (ja) | 2000-08-31 | 2000-08-31 | 脈波解析装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000262600A JP4401014B2 (ja) | 2000-08-31 | 2000-08-31 | 脈波解析装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2002065622A JP2002065622A (ja) | 2002-03-05 |
JP4401014B2 true JP4401014B2 (ja) | 2010-01-20 |
Family
ID=18750258
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000262600A Expired - Fee Related JP4401014B2 (ja) | 2000-08-31 | 2000-08-31 | 脈波解析装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4401014B2 (ja) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4052110B2 (ja) * | 2002-12-05 | 2008-02-27 | 株式会社デンソー | 振動波判定装置 |
US7099714B2 (en) * | 2003-03-31 | 2006-08-29 | Medtronic, Inc. | Biomedical signal denoising techniques |
CN100457023C (zh) * | 2006-09-22 | 2009-02-04 | 北京工业大学 | 一种基于提升小波的脉搏信号处理方法 |
JP6206201B2 (ja) * | 2014-01-17 | 2017-10-04 | トヨタ自動車株式会社 | ステータ通電加熱装置 |
-
2000
- 2000-08-31 JP JP2000262600A patent/JP4401014B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2002065622A (ja) | 2002-03-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9737266B2 (en) | Systems and methods for determining respiration information from a photoplethysmograph | |
US6905470B2 (en) | Apparatus and method for detecting heartbeat using PPG | |
US6024700A (en) | System and method for detecting a thought and generating a control instruction in response thereto | |
US6393311B1 (en) | Method, apparatus and system for removing motion artifacts from measurements of bodily parameters | |
EP2757944B1 (en) | Systems and methods for determining respiration information from a photoplethysmograph | |
US8880576B2 (en) | Systems and methods for determining respiration information from a photoplethysmograph | |
US8655436B2 (en) | Heart rate meter and heart beat detecting method | |
US9402554B2 (en) | Systems and methods for determining respiration information from a photoplethysmograph | |
US9693709B2 (en) | Systems and methods for determining respiration information from a photoplethysmograph | |
JP6190466B2 (ja) | 生体信号測定器及び接触状態推定方法 | |
CN103690152B (zh) | 一种基于脉搏解析的动脉弹性功能评估装置 | |
US20130191035A1 (en) | Method and system for detection and rejection of motion/noise artifacts in physiological measurements | |
JP2004261366A (ja) | 生体状態検出装置及びセンサ並びに生体状態検出方法 | |
EP2854620A1 (en) | Narrow band feature extraction from cardiac signals | |
US20150065896A1 (en) | Pulsation detecting device, electronic apparatus, and program | |
US20140257124A1 (en) | Atrial fibrillation analyzer and program | |
JP2008132012A (ja) | 脈波検出装置 | |
US20130172686A1 (en) | Systems and methods for determining physiological information using autocorrelation with gaps | |
Pu et al. | Novel tailoring algorithm for abrupt motion artifact removal in photoplethysmogram signals | |
US20100074482A1 (en) | Method for detecting bio signal features in the presence of noise | |
JP4401014B2 (ja) | 脈波解析装置 | |
JP2011194217A (ja) | 振動測定法におけるアーティファクトの周波数スペクトルの使用法 | |
JP4122003B2 (ja) | 心拍や呼吸等の生体信号の抽出法及び装置 | |
JP4731031B2 (ja) | 睡眠解析装置及びプログラム並びに記録媒体 | |
Ram et al. | Use of multi-scale principal component analysis for motion artifact reduction of PPG signals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20061025 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20090715 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090721 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090916 |
|
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: 20091027 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20091027 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121106 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131106 Year of fee payment: 4 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |