JPS63270027A - R波検出装置 - Google Patents
R波検出装置Info
- Publication number
- JPS63270027A JPS63270027A JP62105387A JP10538787A JPS63270027A JP S63270027 A JPS63270027 A JP S63270027A JP 62105387 A JP62105387 A JP 62105387A JP 10538787 A JP10538787 A JP 10538787A JP S63270027 A JPS63270027 A JP S63270027A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- change
- rate
- sampling data
- synchronization signal
- 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.)
- Pending
Links
- 230000008859 change Effects 0.000 claims abstract description 56
- 238000005070 sampling Methods 0.000 claims abstract description 41
- 230000001360 synchronised effect Effects 0.000 claims abstract description 20
- 238000001514 detection method Methods 0.000 claims description 9
- 230000036982 action potential Effects 0.000 claims description 7
- 238000003745 diagnosis Methods 0.000 abstract description 3
- 230000000737 periodic effect Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000005236 sound signal Effects 0.000 description 2
- 206010060965 Arterial stenosis Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 210000001765 aortic valve Anatomy 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
技術分野
本発明は、心電センサから出力される心電信号に基づい
て心電図のR波に同期した同期信号を高い精度で検出す
るR波検出装置に関するものである。
て心電図のR波に同期した同期信号を高い精度で検出す
るR波検出装置に関するものである。
従来技術
心臓の活動電位の変化を表す心電図は、患者の健康状態
を表す重要な指標であり、医学的な判断の基礎情報の一
つである。なかでも、最も検出し易いR波は、心電図全
体の波形を識別・解析するために用いられるのみならず
、心音や脈波、コロトコフ(KOROTOKOFF)音
等の採込み範囲を定める際のタイミング信号などとして
、幅広く利用されている。このようなことから、心臓の
活動電位を検出する心電センサから出力される心電信号
に基づいて、上記R波に同期した同期信号を検出するこ
とが行われている。
を表す重要な指標であり、医学的な判断の基礎情報の一
つである。なかでも、最も検出し易いR波は、心電図全
体の波形を識別・解析するために用いられるのみならず
、心音や脈波、コロトコフ(KOROTOKOFF)音
等の採込み範囲を定める際のタイミング信号などとして
、幅広く利用されている。このようなことから、心臓の
活動電位を検出する心電センサから出力される心電信号
に基づいて、上記R波に同期した同期信号を検出するこ
とが行われている。
発明が解決しようとする問題点
ところで、かかる従来のR波検出装置は、一般に、心電
センサから出力される心電信号を所定のサンプリング周
期で読み込み、そのサンプリングデータの変化率からR
波に同期した同期信号を決定していたため、必ずしも充
分な検出精度が得られないという問題があった。すなわ
ち、心臓の活動電位には、通常、商用交流電源に起因す
る周期性のノイズが含まれているところから、心電セン
サから出力される心電信号はその商用交流電源の周波数
に同期して振動する一方、上記サンプリング周期はこの
商用交流電源の周波数の周期よりも極めて短いため、上
記周期性ノイズに起因してサンプリングデータの変化率
、更には同期信号にばらつきを生じてしまうのである。
センサから出力される心電信号を所定のサンプリング周
期で読み込み、そのサンプリングデータの変化率からR
波に同期した同期信号を決定していたため、必ずしも充
分な検出精度が得られないという問題があった。すなわ
ち、心臓の活動電位には、通常、商用交流電源に起因す
る周期性のノイズが含まれているところから、心電セン
サから出力される心電信号はその商用交流電源の周波数
に同期して振動する一方、上記サンプリング周期はこの
商用交流電源の周波数の周期よりも極めて短いため、上
記周期性ノイズに起因してサンプリングデータの変化率
、更には同期信号にばらつきを生じてしまうのである。
問題点を解決するための手段
本発明は以上の事情を背景として為されたものであり、
その目的とするところは、心電信号に混入する商用交流
電源の周期性ノイズに拘らずR波に同期した同期信号を
高い精度で検出することにある。
その目的とするところは、心電信号に混入する商用交流
電源の周期性ノイズに拘らずR波に同期した同期信号を
高い精度で検出することにある。
そして、かかる目的を達成するため、本発明は、心臓の
活動電位を検出する心電センサから出力される心電信号
に基づいて心電図のR波に同期した同期信号を検出する
装置であって、(a)前記心電信号を予め定められたサ
ンプリング周期で読み込む読込み手段と、(blその読
込み手段によって読み込まれたサンプリングデータのう
ち、商用交流電源の周波数に対応する周期で読み込まれ
た二つのサンプリングデータの変化率を算出する算出手
段と、(C1その算出手段によって算出された前記変化
率に基づいて前記サンプリングデータの中から前記同期
信号を決定する同期信号決定手段とを有することを特徴
とする。
活動電位を検出する心電センサから出力される心電信号
に基づいて心電図のR波に同期した同期信号を検出する
装置であって、(a)前記心電信号を予め定められたサ
ンプリング周期で読み込む読込み手段と、(blその読
込み手段によって読み込まれたサンプリングデータのう
ち、商用交流電源の周波数に対応する周期で読み込まれ
た二つのサンプリングデータの変化率を算出する算出手
段と、(C1その算出手段によって算出された前記変化
率に基づいて前記サンプリングデータの中から前記同期
信号を決定する同期信号決定手段とを有することを特徴
とする。
作用および発明の効果
このようなR波検出装置においては、読込み手段により
所定のサンプリング周期で読み込まれたサンプリングデ
ータのうち、商用交流電源の周波数に対応する周期で読
み込まれた二つのサンプリングデータ、すなわち商用交
流電源の周期変化に対する位相が互いに一致するサンプ
リングデータから、その変化率が算出手段によって算出
される。
所定のサンプリング周期で読み込まれたサンプリングデ
ータのうち、商用交流電源の周波数に対応する周期で読
み込まれた二つのサンプリングデータ、すなわち商用交
流電源の周期変化に対する位相が互いに一致するサンプ
リングデータから、その変化率が算出手段によって算出
される。
このため、心電センサから出力される心電信号に商用交
流電源に起因する周期性ノイズが混入していても、上記
変化率がその周期性ノイズに影響されることはない。そ
して、この変化率に基づいてR波に同期した同期信号が
同期信号決定手段により決定されるため、その決定され
た同期信号はR波に対して常に高い精度で同期させられ
ることとなる。
流電源に起因する周期性ノイズが混入していても、上記
変化率がその周期性ノイズに影響されることはない。そ
して、この変化率に基づいてR波に同期した同期信号が
同期信号決定手段により決定されるため、その決定され
た同期信号はR波に対して常に高い精度で同期させられ
ることとなる。
ここで、前記読込み手段を、前記商用交流電源の周波数
に対応する周期よりも小さいサンプリング周期で前記心
電信号を読み込むように構成する一方、前記算出手段を
、そのサンプリング周期で前記心電信号が読み込まれる
毎に前記変化率を算出するように構成し、その算出され
た全ての変化率に基づいて同期信号を決定させるように
すれば、その同期信号を一層高い精度で検出することが
できるが、読込み手段により商用交流電源の周波数に対
応する周期と同じ周期でサンプリングデータを読み込ま
せるようにしても、商用交流電源の周期性ノイズに起因
する同期信号のばらつきを防止するという本発明の効果
は得られる。
に対応する周期よりも小さいサンプリング周期で前記心
電信号を読み込むように構成する一方、前記算出手段を
、そのサンプリング周期で前記心電信号が読み込まれる
毎に前記変化率を算出するように構成し、その算出され
た全ての変化率に基づいて同期信号を決定させるように
すれば、その同期信号を一層高い精度で検出することが
できるが、読込み手段により商用交流電源の周波数に対
応する周期と同じ周期でサンプリングデータを読み込ま
せるようにしても、商用交流電源の周期性ノイズに起因
する同期信号のばらつきを防止するという本発明の効果
は得られる。
また、前記同期信号決定手段は、例えば、前記変化率の
絶対値が予め求められた閾値よりも大きい場合に、その
変化率が正のときには前記二つのサンプリングデータの
うち後から読み込まれたサンプリングデータを前記同期
信号として決定し、その変化率が負のときには前記二つ
のサンプリングデータのうち先に読み込まれたサンプリ
ングデータを前記同期信号として決定するように構成さ
れる。
絶対値が予め求められた閾値よりも大きい場合に、その
変化率が正のときには前記二つのサンプリングデータの
うち後から読み込まれたサンプリングデータを前記同期
信号として決定し、その変化率が負のときには前記二つ
のサンプリングデータのうち先に読み込まれたサンプリ
ングデータを前記同期信号として決定するように構成さ
れる。
なお、R波に同期した同期信号とは、心電図のR波のピ
ークに対して一定の関係にある信号のことで、必ずしも
R波のピークに対応する信号を意味するものではない。
ークに対して一定の関係にある信号のことで、必ずしも
R波のピークに対応する信号を意味するものではない。
実施例
以下、本発明の一実施例を図面に基づいて詳細に説明す
る。
る。
第1図において、IOは心臓の活動電位を検出する心電
センサとしての心電計で、生体の胸部等に接着させられ
る複数の電極12を備えている。
センサとしての心電計で、生体の胸部等に接着させられ
る複数の電極12を備えている。
心電計10からは上記活動電位に対応する心電信号SE
が出力され、この心電信号SEはA/Dコンバータ14
においてデジタル信号SEDに変換された後、タイミン
グ信号STに従ってCPUl6に読み込まれる。CPU
16には、押釦PBおよびクロック信号源18からそれ
ぞれ起動信号SAおよび所定周波数のパルス信号spが
供給されるとともに、図示は省略するが、心臓から発生
する心音を表す心音信号が心音計等から供給されるよう
になっている。
が出力され、この心電信号SEはA/Dコンバータ14
においてデジタル信号SEDに変換された後、タイミン
グ信号STに従ってCPUl6に読み込まれる。CPU
16には、押釦PBおよびクロック信号源18からそれ
ぞれ起動信号SAおよび所定周波数のパルス信号spが
供給されるとともに、図示は省略するが、心臓から発生
する心音を表す心音信号が心音計等から供給されるよう
になっている。
上記CPUI 6はRAM20およびROM22と共に
マイクロコンピュータを構成しており、RAM20の一
時記憶機能を利用しつつROM22に予め記憶されたプ
ログラムに従って信号処理を行う。第2図は、かかる信
号処理ロジックを示すフローチャートの一例であり、以
下、このフローチャートに従って本実施例の作動を説明
する。
マイクロコンピュータを構成しており、RAM20の一
時記憶機能を利用しつつROM22に予め記憶されたプ
ログラムに従って信号処理を行う。第2図は、かかる信
号処理ロジックを示すフローチャートの一例であり、以
下、このフローチャートに従って本実施例の作動を説明
する。
先ず、図示しない電源スィッチが投入されると図示しな
い初期化ステップを経てステップs1が実行され、押6
0PBが押圧操作されたが否が、換言すれば起動信号S
Aが供給されているが否がか判断される。起動信号SA
が供給されると、次にステップS2が実行され、クロッ
ク信号源18がら供給されているパルス信号SPを計数
するタイマの計数内容tがリセフトされるとともに、続
いてステップS3が実行される。
い初期化ステップを経てステップs1が実行され、押6
0PBが押圧操作されたが否が、換言すれば起動信号S
Aが供給されているが否がか判断される。起動信号SA
が供給されると、次にステップS2が実行され、クロッ
ク信号源18がら供給されているパルス信号SPを計数
するタイマの計数内容tがリセフトされるとともに、続
いてステップS3が実行される。
ステップS3においては、予め定められた一定のサンプ
リング周期、本実施例では約0.4.44 m秒(周波
数で2250Hz)でタイミング信号STをA/Dコン
バータ14へ出力することにより、そのサンプリング周
期で前記デジタル信号SEDを読み込み、次にステップ
S4が実行される。このステップS4は、ステップS3
において読み込まれたデジタル信号SED、すなわちサ
ンプリングデータの変化率Siを次式(1)に従って算
出する。
リング周期、本実施例では約0.4.44 m秒(周波
数で2250Hz)でタイミング信号STをA/Dコン
バータ14へ出力することにより、そのサンプリング周
期で前記デジタル信号SEDを読み込み、次にステップ
S4が実行される。このステップS4は、ステップS3
において読み込まれたデジタル信号SED、すなわちサ
ンプリングデータの変化率Siを次式(1)に従って算
出する。
ここで、f (i+37)、f (i)は、何れも
読み込まれたデジタル信号S E D <、、+ の信
号値f(n)であり、変化率Stは、前記サンプリング
周期の37倍の周期で読み込まれた二つのデジタル信号
S E D (11およびS E D (i+371の
信号値f(i)およびf (i+37)に基づいて算
出される。これは、心臓の活動電位には、通常、商用交
流電源に起因する周期性のノイズが含まれているため、
その商用交流電源の周波数に対応した周期で読み込まれ
たデジタル信号SEDの信号値を比較してその変化率S
tを算出するためで、37は、商用交流電源の周期T(
この実施例では1760秒)を前記サンプリング周期(
0,444m秒)で割算することによって求められた値
である。これにより、デジタル信号S E D l、、
) の信号値f (n)が第3図に示されているよう
に商用交流電源の周期Tで変動しても、変化率Siを算
出するための基礎となるデジタル信号S E D t=
、の信号値f (i)とデジタル信号S E D li
。37〉の信号値f (i+37)とはその周期Tに
対する位相が互いに一致しているため、変化率Siが商
用交流電源による周期性ノイズに影響されることはない
。なお、この段階においては未だデジタル信号SEDは
一つしかないため、前記ステップS3が繰り返されて3
8個のデジタル信号S E D +++ 〜S E D
(1111が読み込まれることにより、初めて最初の
デジタル信号5EDLl、の信号値r (1)と38番
目のデジタル信号S E D (3111の信号値f(
38)とに基づいて変化率S、が算出される。
読み込まれたデジタル信号S E D <、、+ の信
号値f(n)であり、変化率Stは、前記サンプリング
周期の37倍の周期で読み込まれた二つのデジタル信号
S E D (11およびS E D (i+371の
信号値f(i)およびf (i+37)に基づいて算
出される。これは、心臓の活動電位には、通常、商用交
流電源に起因する周期性のノイズが含まれているため、
その商用交流電源の周波数に対応した周期で読み込まれ
たデジタル信号SEDの信号値を比較してその変化率S
tを算出するためで、37は、商用交流電源の周期T(
この実施例では1760秒)を前記サンプリング周期(
0,444m秒)で割算することによって求められた値
である。これにより、デジタル信号S E D l、、
) の信号値f (n)が第3図に示されているよう
に商用交流電源の周期Tで変動しても、変化率Siを算
出するための基礎となるデジタル信号S E D t=
、の信号値f (i)とデジタル信号S E D li
。37〉の信号値f (i+37)とはその周期Tに
対する位相が互いに一致しているため、変化率Siが商
用交流電源による周期性ノイズに影響されることはない
。なお、この段階においては未だデジタル信号SEDは
一つしかないため、前記ステップS3が繰り返されて3
8個のデジタル信号S E D +++ 〜S E D
(1111が読み込まれることにより、初めて最初の
デジタル信号5EDLl、の信号値r (1)と38番
目のデジタル信号S E D (3111の信号値f(
38)とに基づいて変化率S、が算出される。
このようにして最初の変化率S1が算出されると、次に
ステップS5が実行されてその絶対値IS11を記憶し
た後、ステップS6において前記タイマの計数内容tが
1.8秒を経過したか否かが判断される。最初の変化率
S、を求めるのに必要な時間は約1/60秒であるため
、この段階での判断はNOであり、ステップ83以下が
再び実行されて39番目のデジタル信号S E D (
391が読み込まれるとともに、その信号値f(39)
と2番目に読み込まれたデジタル信号S E D tz
、の信号値「 (2)とから変化率S2が算出され、そ
の絶対値1s21が記憶される。この時の所要時間は前
記サンプリング周期と略同じで約0.444m秒である
。したがって、タイマの計数内容tが1.8秒を経過す
ることはなく、引き続きステップ83以下の実行が繰り
返されてその都度変化率Siが算出され、その絶対値1
si lが記憶される。
ステップS5が実行されてその絶対値IS11を記憶し
た後、ステップS6において前記タイマの計数内容tが
1.8秒を経過したか否かが判断される。最初の変化率
S、を求めるのに必要な時間は約1/60秒であるため
、この段階での判断はNOであり、ステップ83以下が
再び実行されて39番目のデジタル信号S E D (
391が読み込まれるとともに、その信号値f(39)
と2番目に読み込まれたデジタル信号S E D tz
、の信号値「 (2)とから変化率S2が算出され、そ
の絶対値1s21が記憶される。この時の所要時間は前
記サンプリング周期と略同じで約0.444m秒である
。したがって、タイマの計数内容tが1.8秒を経過す
ることはなく、引き続きステップ83以下の実行が繰り
返されてその都度変化率Siが算出され、その絶対値1
si lが記憶される。
その後、タイマの計数内容tが1.8秒を経過するとス
テップS7が実行され、前記ステップS5において記憶
された多数の変化率l511の中からその値が最も大き
い最大変化率max l S i lを求める。そして
、続くステップS8においては、次式(2)に従ってそ
の最大変化率maxlsilに0゜8を掛算することに
より、閾値Cを算出する。
テップS7が実行され、前記ステップS5において記憶
された多数の変化率l511の中からその値が最も大き
い最大変化率max l S i lを求める。そして
、続くステップS8においては、次式(2)に従ってそ
の最大変化率maxlsilに0゜8を掛算することに
より、閾値Cを算出する。
C= maxl S i l Xo、8
・・・(2)この閾値Cは、その後に心電計10から
出力される心電信号SEに基づいて、心電図のR波に同
期したデジタル信号SED、すなわち同期信号SE D
(1)を検出するための基準となるものである。
・・・(2)この閾値Cは、その後に心電計10から
出力される心電信号SEに基づいて、心電図のR波に同
期したデジタル信号SED、すなわち同期信号SE D
(1)を検出するための基準となるものである。
この閾値Cを求めるために、前記ステップ83〜S5を
1.8秒以上繰り返したのは、心臓の拍動周期は通常7
0回/分以上であるため、1.8秒あれば個体差に関係
なくR波に対応する心電信号SEが供給されるからであ
る。また、変化率1si lに基づいて閾値Cを決定し
たのは、心電図のR波は、第4図に示されているように
デジタル信号SE D 、、、の信号値f (n)が
最も急激に変化する部分だからであり、最大変化率ma
xlsilに0゜8を掛算したのは、R波に同期した同
期信号5ED(8) を確実に検出するためである。
1.8秒以上繰り返したのは、心臓の拍動周期は通常7
0回/分以上であるため、1.8秒あれば個体差に関係
なくR波に対応する心電信号SEが供給されるからであ
る。また、変化率1si lに基づいて閾値Cを決定し
たのは、心電図のR波は、第4図に示されているように
デジタル信号SE D 、、、の信号値f (n)が
最も急激に変化する部分だからであり、最大変化率ma
xlsilに0゜8を掛算したのは、R波に同期した同
期信号5ED(8) を確実に検出するためである。
このようにして閾値Cが決定されると、次にステップS
9および310が実行され、前記ステップS3.S4と
同様にして、新たにデジタル信号SEDが前述したサン
プリング周期(0,444m秒)で読み込まれるととも
に、その読み込まれたデジタル信号SEDの変化率Si
が前記(1)式に従って算出される。そして、ステップ
311においては、算出された変化率の絶対値1si
lが前記閾値Cよりも大きいか否かが判断され、小さい
場合には上記ステップS9および310が繰り返される
が、大きい場合には次のステップS12が実行されて心
電図のR波に同期した同期信号5ED(R)が決定され
る。すなわち、このステップS12においては、先ず、
変化率Siの正負が判断され、変化率Stが正の場合に
は後から読み込まれたデジタル信号S E D (i。
9および310が実行され、前記ステップS3.S4と
同様にして、新たにデジタル信号SEDが前述したサン
プリング周期(0,444m秒)で読み込まれるととも
に、その読み込まれたデジタル信号SEDの変化率Si
が前記(1)式に従って算出される。そして、ステップ
311においては、算出された変化率の絶対値1si
lが前記閾値Cよりも大きいか否かが判断され、小さい
場合には上記ステップS9および310が繰り返される
が、大きい場合には次のステップS12が実行されて心
電図のR波に同期した同期信号5ED(R)が決定され
る。すなわち、このステップS12においては、先ず、
変化率Siの正負が判断され、変化率Stが正の場合に
は後から読み込まれたデジタル信号S E D (i。
1.)を同期信号S E D unとして決定する一方
、変化率Siが負の場合には先に読み込まれたデジタル
信号S E D t+、を同期信号S E D (11
1として決定する。本実施例では、これ等の一連の信号
処理ロジックのうち上記ステップS9を実行する部分が
読込み手段に相当し、ステップSIOを実行する部分が
算出手段に相当し、ステップSllおよびS12を実行
する部分が同期信号決定手段に相当する。
、変化率Siが負の場合には先に読み込まれたデジタル
信号S E D t+、を同期信号S E D (11
1として決定する。本実施例では、これ等の一連の信号
処理ロジックのうち上記ステップS9を実行する部分が
読込み手段に相当し、ステップSIOを実行する部分が
算出手段に相当し、ステップSllおよびS12を実行
する部分が同期信号決定手段に相当する。
なお、心電図のQ波のピークからS波のピークまでの時
間は通常100m秒程度で、変化率Stを算出するため
の基礎となるデジタル信号5EDO)およびS E D
(i+3?)の読込み周期(1/60秒’=17m秒
)に比較して充分に大きいため、かかる変化率Siに基
づいてR波に同期する同期信号5ED(。を決定するこ
とは充分に可能なのである。また、本実施例の閾値Cは
、前記(2)式から明らかなように最大変化率−,1S
ilに0.8を掛算して求められたもので、比較的小さ
いため、上記同期信号S E D 、、、 は、第4図
に示されているように、通常はデジタル信号S E D
Li+ が心電図のQ波のピークに対応する場合のデ
ジタル信号S E D (i+3’7)に決定される。
間は通常100m秒程度で、変化率Stを算出するため
の基礎となるデジタル信号5EDO)およびS E D
(i+3?)の読込み周期(1/60秒’=17m秒
)に比較して充分に大きいため、かかる変化率Siに基
づいてR波に同期する同期信号5ED(。を決定するこ
とは充分に可能なのである。また、本実施例の閾値Cは
、前記(2)式から明らかなように最大変化率−,1S
ilに0.8を掛算して求められたもので、比較的小さ
いため、上記同期信号S E D 、、、 は、第4図
に示されているように、通常はデジタル信号S E D
Li+ が心電図のQ波のピークに対応する場合のデ
ジタル信号S E D (i+3’7)に決定される。
そして、このようにしてR波に同期した同期信号S E
D (11が決定されると、次にステップS13の心
音データ採込みルーチンが実行され、同期信号5EDf
。を基準として例えば心音の中の大動脈弁閉鎖に起因す
る第1音を決定するとともに、その第1音を基準として
予め定められた一定時間だけ心音を表す心音信号が採り
込まれる。その後、ステップS14では、ステップS1
3において採り込まれた心音データの数が予め定められ
た一定数αを超えたか否かが判断され、心音データの数
が一定数αに達するまで上記ステップ89〜S13が繰
り返される。これにより、心臓の一拍毎に同期信号S
B D (R)を基準として第1音が決定されるととも
に、その第1音に基づいて心音データが採り込まれ、そ
のデータ数がαに達すると、以後、その心音データに基
づいて例えば冠状動脈の狭窄異常等の自動診断が行われ
る。
D (11が決定されると、次にステップS13の心
音データ採込みルーチンが実行され、同期信号5EDf
。を基準として例えば心音の中の大動脈弁閉鎖に起因す
る第1音を決定するとともに、その第1音を基準として
予め定められた一定時間だけ心音を表す心音信号が採り
込まれる。その後、ステップS14では、ステップS1
3において採り込まれた心音データの数が予め定められ
た一定数αを超えたか否かが判断され、心音データの数
が一定数αに達するまで上記ステップ89〜S13が繰
り返される。これにより、心臓の一拍毎に同期信号S
B D (R)を基準として第1音が決定されるととも
に、その第1音に基づいて心音データが採り込まれ、そ
のデータ数がαに達すると、以後、その心音データに基
づいて例えば冠状動脈の狭窄異常等の自動診断が行われ
る。
ここで、上記第1音を決定するための基準となる同期信
号S E D (Ill は、予め定められたサンプリ
ング周期で読み込まれた多数のデジタル信号SEDい、
のうち、商用交流電源の周期Tで読み込まれてその商用
交流電源の周期変化に対する位相が互いに一致する二つ
のデジタル信号S E D (i+およびS E D
(io37)の信号値f (i)およびf (i+3
7)の変化率Siに基づいて決定されているため、心電
計10から出力される心電信号SE、更にはデジタル信
号S E D (−1に商用交流電源に起因する周期性
ノイズが混入していても、上記変化率Si、更にはその
変化率Siに基づいて決定される同期信号S E D
+u がその周期性ノイズに影響されることはなく、R
波に対して常に高い精度で同期させられる。したがって
、かかる同期信号S E D 、、、を基準として決定
される第1音、更にはその第1音に基づいて採り込まれ
る心音データの信輔性が高くなり、その心音データに基
づいてより正確な自動診断が行われるようになるのであ
る。
号S E D (Ill は、予め定められたサンプリ
ング周期で読み込まれた多数のデジタル信号SEDい、
のうち、商用交流電源の周期Tで読み込まれてその商用
交流電源の周期変化に対する位相が互いに一致する二つ
のデジタル信号S E D (i+およびS E D
(io37)の信号値f (i)およびf (i+3
7)の変化率Siに基づいて決定されているため、心電
計10から出力される心電信号SE、更にはデジタル信
号S E D (−1に商用交流電源に起因する周期性
ノイズが混入していても、上記変化率Si、更にはその
変化率Siに基づいて決定される同期信号S E D
+u がその周期性ノイズに影響されることはなく、R
波に対して常に高い精度で同期させられる。したがって
、かかる同期信号S E D 、、、を基準として決定
される第1音、更にはその第1音に基づいて採り込まれ
る心音データの信輔性が高くなり、その心音データに基
づいてより正確な自動診断が行われるようになるのであ
る。
以上、本発明の一実施例を図面に基づいて詳細に説明し
たが、本発明は他の態様で実施することもできる。
たが、本発明は他の態様で実施することもできる。
例えば、前記実施例では予め最大変化率l5i1に基づ
いて閾値Cが求められ、その閾値Cを基準として同期信
号S E D (11を決定するようになっているが、
変化率Siが正となる場合のみから閾値C1更には同期
信号S E D (Ill を決定するようにしたり、
或いは変化率Siが負となる場合のみから閾値C1更に
は同期信号S E D n+ を決定するようにしたり
することもできる。なお、最大変化率lsi lに対す
る閾値Cの割合を適宜変更し得ることは勿論である。
いて閾値Cが求められ、その閾値Cを基準として同期信
号S E D (11を決定するようになっているが、
変化率Siが正となる場合のみから閾値C1更には同期
信号S E D (Ill を決定するようにしたり、
或いは変化率Siが負となる場合のみから閾値C1更に
は同期信号S E D n+ を決定するようにしたり
することもできる。なお、最大変化率lsi lに対す
る閾値Cの割合を適宜変更し得ることは勿論である。
また、変化率StはR波のピークを境として急激に減少
するため、この変化率Siの変化からR波のピークに対
応する、或いはR波のピークに対・ して一定の関係を有する同期信号S E D (R1を
決定するようにすることもできる。その場合には、閾値
Cを予め求めておく必要はない。
するため、この変化率Siの変化からR波のピークに対
応する、或いはR波のピークに対・ して一定の関係を有する同期信号S E D (R1を
決定するようにすることもできる。その場合には、閾値
Cを予め求めておく必要はない。
また、前記実施例では同期信号S E D (Illを
基準として心音の第1音を決定し、その第1音に基づい
て心音データを採り込む場合について説明したが、脈波
やコロトコフ音等を採り込む場合、或いは心電図の波形
の識別・解析等を行う場合にも本発明は同様に適用され
得る。
基準として心音の第1音を決定し、その第1音に基づい
て心音データを採り込む場合について説明したが、脈波
やコロトコフ音等を採り込む場合、或いは心電図の波形
の識別・解析等を行う場合にも本発明は同様に適用され
得る。
その他−々例示はしないが、本発明はその精神を逸脱す
ることなく当業者の知識に基づいて種々の変更、改良を
加えた態様で実施することができる。
ることなく当業者の知識に基づいて種々の変更、改良を
加えた態様で実施することができる。
第1図は本発明の一実施例であるR波検出装置の構成を
説明するブロック線図である。第2図は第1図の装置の
作動を説明するフローチャートである。第3図はデジタ
ル信号に混入した商用交流電源の周期性ノイズを説明す
る図である。第4図はデジタル信号の信号値の変化を示
す図である。 10:心電計(心電センサ) 16 :CPU 20 :RAM22:RO
M 5E:心電信号 SED:デジタル信号(サンプリングデータ)SED(
R1:同期信号 Si:変化率ステップS9:読込み手
段 ステップS10:算出手段 ステップSll、S12:同朋信号決定手段出願人
コーリン電子株式会社 第3図 第4図
説明するブロック線図である。第2図は第1図の装置の
作動を説明するフローチャートである。第3図はデジタ
ル信号に混入した商用交流電源の周期性ノイズを説明す
る図である。第4図はデジタル信号の信号値の変化を示
す図である。 10:心電計(心電センサ) 16 :CPU 20 :RAM22:RO
M 5E:心電信号 SED:デジタル信号(サンプリングデータ)SED(
R1:同期信号 Si:変化率ステップS9:読込み手
段 ステップS10:算出手段 ステップSll、S12:同朋信号決定手段出願人
コーリン電子株式会社 第3図 第4図
Claims (3)
- (1)心臓の活動電位を検出する心電センサから出力さ
れる心電信号に基づいて心電図のR波に同期した同期信
号を検出する装置であって、 前記心電信号を予め定められたサンプリング周期で読み
込む読込み手段と、 該読込み手段によって読み込まれたサンプリングデータ
のうち、商用交流電源の周波数に対応する周期で読み込
まれた二つのサンプリングデータの変化率を算出する算
出手段と、 該算出手段によって算出された前記変化率に基づいて前
記サンプリングデータの中から前記同期信号を決定する
同期信号決定手段と を有することを特徴とするR波検出装置。 - (2)前記読込み手段は、前記商用交流電源の周波数に
対応する周期よりも小さいサンプリング周期で前記心電
信号を読み込むものであり、前記算出手段は、該サンプ
リング周期で前記心電信号が読み込まれる毎に前記変化
率を算出するものである特許請求の範囲第1項に記載の
R波検出装置。 - (3)前記同期信号決定手段は、前記変化率の絶対値が
予め求められた閾値よりも大きい場合に、該変化率が正
のときには前記二つのサンプリングデータのうち後から
読み込まれたサンプリングデータを前記同期信号として
決定し、該変化率が負のときには前記二つのサンプリン
グデータのうち先に読み込まれたサンプリングデータを
前記同期信号として決定するものである特許請求の範囲
第1項または第2項に記載のR波検出装置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP62105387A JPS63270027A (ja) | 1987-04-28 | 1987-04-28 | R波検出装置 |
| FR8805471A FR2614441B1 (fr) | 1987-04-24 | 1988-04-25 | Appareil d'enregistrement d'images en couleurs |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP62105387A JPS63270027A (ja) | 1987-04-28 | 1987-04-28 | R波検出装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPS63270027A true JPS63270027A (ja) | 1988-11-08 |
Family
ID=14406244
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP62105387A Pending JPS63270027A (ja) | 1987-04-24 | 1987-04-28 | R波検出装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS63270027A (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20030070479A (ko) * | 2002-02-25 | 2003-08-30 | 주식회사 이디지털메드 | 심장박동 감시 시스템 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61245024A (ja) * | 1985-04-23 | 1986-10-31 | Yamatake Honeywell Co Ltd | 電磁流量計 |
-
1987
- 1987-04-28 JP JP62105387A patent/JPS63270027A/ja active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61245024A (ja) * | 1985-04-23 | 1986-10-31 | Yamatake Honeywell Co Ltd | 電磁流量計 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20030070479A (ko) * | 2002-02-25 | 2003-08-30 | 주식회사 이디지털메드 | 심장박동 감시 시스템 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4446873A (en) | Method and apparatus for detecting heart sounds | |
| US4649930A (en) | Apparatus for beat buffering techniques varified by arrhythmias detection for stopaction frames of cardiac function | |
| US4549552A (en) | Heart sound detector and cardiac cycle data are combined for diagnostic reliability | |
| JP3559284B2 (ja) | 仮想トリガを使用する心拍同期化用パルス酸素計 | |
| US4548204A (en) | Apparatus for monitoring cardiac activity via ECG and heart sound signals | |
| US4546777A (en) | Heart sound detector and synchronization for diagnostics | |
| US5033472A (en) | Method of and apparatus for analyzing propagation of arterial pulse waves through the circulatory system | |
| US7137955B2 (en) | Methods and systems for distal recording of phonocardiographic signals | |
| US4993421A (en) | Cardiac monitoring system | |
| CN115023182B (zh) | 用于同步光学光电描记(ppg)信号与心电图(ecg)信号的计算机实现的方法 | |
| JPS6244265A (ja) | 患者の呼吸監視装置 | |
| US3878832A (en) | Method and apparatus for detecting and quantifying cardiovascular murmurs and the like | |
| US3902479A (en) | Method and apparatus for heartbeat rate monitoring | |
| US3581735A (en) | Heartbeat frequency determining apparatus and method | |
| JP2011212364A (ja) | 心音測定装置 | |
| JP3415796B2 (ja) | 周期性生体情報測定装置 | |
| JP2001198094A (ja) | 脈拍数検出装置 | |
| US20260053418A1 (en) | Heart sound segmentation method, apparatus, electronic device, and readable storage medium | |
| CN100448395C (zh) | 血压测量装置 | |
| JP5203210B2 (ja) | パルスオキシメトリーのための信号処理 | |
| JPS63270027A (ja) | R波検出装置 | |
| JPS63290544A (ja) | 狭窄音検出装置 | |
| JP2000093399A (ja) | 心音計 | |
| JP2798737B2 (ja) | 心拍性同期信号採取装置 | |
| FI129789B (en) | A method, apparatus and computer program product for determining pulse transmission time |