JP2001333890A - Method of computing pulse rate - Google Patents

Method of computing pulse rate

Info

Publication number
JP2001333890A
JP2001333890A JP2000157258A JP2000157258A JP2001333890A JP 2001333890 A JP2001333890 A JP 2001333890A JP 2000157258 A JP2000157258 A JP 2000157258A JP 2000157258 A JP2000157258 A JP 2000157258A JP 2001333890 A JP2001333890 A JP 2001333890A
Authority
JP
Japan
Prior art keywords
pulse
period
pulse wave
value
cycle
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.)
Withdrawn
Application number
JP2000157258A
Other languages
Japanese (ja)
Inventor
Hiroyuki Ibe
浩行 井邊
Kaoru Fukuyori
薫 福頼
Akihiro Michimori
章弘 道盛
Hiroshi Hagiwara
啓 萩原
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Panasonic Electric Works Co Ltd
Original Assignee
Matsushita Electric Works Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Matsushita Electric Works Ltd filed Critical Matsushita Electric Works Ltd
Priority to JP2000157258A priority Critical patent/JP2001333890A/en
Publication of JP2001333890A publication Critical patent/JP2001333890A/en
Withdrawn legal-status Critical Current

Links

Landscapes

  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

PROBLEM TO BE SOLVED: To figure out a representing value for a period of pulse without being affected by noise. SOLUTION: When pulse rate is figured out by obtaining a period of pulse from output by a pulse sensor which detects pulse, the pulse rate of a period, when there is no signals considered as a pulse, is presumed by using a representing value for a period of the pulse which is calculated from a period of the pulse obtained before. An auto-correlation function for the pulse signals during a pre-determined period is figured out. When the first peak value going into the positive direction is greater than a pre-determined threshold value, the period of the pulse during the period is taken as a representing value for the period of the pulse. Influence caused by noise is avoided by calculating the representing value for the period of the pulse using the auto-correlation function.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は体動などによるノイ
ズが混入した際に脈拍数を算出することができる脈拍数
の算出方法に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a pulse rate calculating method capable of calculating a pulse rate when noise due to body motion or the like is mixed.

【0002】[0002]

【従来の技術】脈拍数を非観血で行うにあたり、脈拍セ
ンサの出力から脈波とみなすことができる信号を抽出し
て、1拍毎の脈波の発生時刻を順次記憶し、この発生時
刻から脈波の周期を算出して、該周期から脈拍数を求め
ることが行われている。
2. Description of the Related Art When a pulse rate is measured in a non-invasive manner, a signal which can be regarded as a pulse wave is extracted from the output of a pulse sensor, and the time of occurrence of each pulse is sequentially stored. Calculating the pulse wave period from the above, and obtaining the pulse rate from the period.

【0003】この場合、脈波を確実に検出することがで
きていれば問題は無いのであるが、実際には体動による
ノイズによって、一時的に脈波とみなせる信号が得られ
なかったりすることがあり、この時には存在しているは
ずの脈波を補完するかたちで周期を算出しなくては、正
確な脈拍数を得ることができない。
In this case, there is no problem as long as the pulse wave can be reliably detected. However, in practice, a signal that can be temporarily regarded as a pulse wave cannot be obtained due to noise due to body movement. In this case, an accurate pulse rate cannot be obtained unless the cycle is calculated in a manner that complements the pulse wave that should exist.

【0004】ここにおいて、特開平9−173308号
公報には、一時的に脈波とみなせる信号が得られなかっ
た期間の脈拍数を推測する方法として、それまでの脈波
周期を基に次の脈波周期の範囲を限定し、脈波である可
能性がある信号までの周期をそれまでの脈波周期に基づ
いて求めた脈波周期代表値から推測することが示されて
いる。
Here, Japanese Patent Application Laid-Open No. 9-173308 discloses a method for estimating a pulse rate during a period in which a signal that can be temporarily regarded as a pulse wave is not obtained, based on the pulse wave period up to that time. It is shown that the range of the pulse wave cycle is limited, and the cycle up to a signal that may be a pulse wave is estimated from a representative pulse wave cycle value obtained based on the previous pulse wave cycle.

【0005】また、上記公報には、脈波とみなせる信号
が得られなかった期間を除くとともにこの期間を挟まな
い最近数回の脈波周期の平均値を脈波周期代表値として
採用することが示されている。また、ここでは脈波とみ
なせる信号が得られなかった期間は、脈波周期代表値か
ら一定率を増減した領域に入らない脈波周期が1つ以上
続いた場合としている。
In the above publication, the average value of several recent pulse wave periods excluding a period in which a signal that can be regarded as a pulse wave is not obtained and not including this period is adopted as a pulse wave period representative value. It is shown. Here, the period during which a signal that can be regarded as a pulse wave is not obtained is a case in which one or more pulse wave periods that do not enter a region where the constant rate is increased or decreased from the pulse wave period representative value continue.

【0006】[0006]

【発明が解決しようとする課題】この場合、脈波周期代
表値が正しければ、これを根拠に算出した脈拍数もほぼ
正確なものとなるが、上記脈波周期代表値に関して上記
公報に示されたものは次の問題を有している。
In this case, if the representative value of the pulse wave cycle is correct, the pulse rate calculated on the basis thereof will be almost accurate. Have the following problems:

【0007】すなわち、体動によって脈波信号にノイズ
が混入した時、常に脈波周期代表値に一定率を増減した
領域に入らない状態が生じるわけではなく、脈波周期代
表値に一定率を増減した領域に入り続けることがあり、
結果的に脈波周期が徐々に長くなる場合がある。図4の
M1〜M4はこの場合を示しており、M3及びM4は本
来の脈波ではなく、ノイズによるものである。
That is, when noise is mixed in the pulse wave signal due to body movement, a state in which the pulse wave cycle representative value does not always enter a region where the constant rate is increased or decreased does not always occur. May continue to enter the increased or decreased area,
As a result, the pulse wave cycle may gradually increase. M1 to M4 in FIG. 4 show this case, and M3 and M4 are not due to the original pulse wave but due to noise.

【0008】このような状態が生じると、本来の脈波で
ある信号M5〜M7が再度得られた時には、脈拍数算定
のために算出している脈波周期を代表する値が、真の脈
波周期からかけ離れた値をとっているために、図4に示
す例では実際上、M6の脈波が無視されてしまう場合が
生じるものであり、この時には正確な脈拍数を得ること
ができない。つまりは、脈波信号と見なせる信号を得る
ことが出来なかった期間の脈拍数を推測するにあたり、
その推測の根拠となる脈波周期代表値そのものがノイズ
による影響を受けてしまっていることから、正確な脈拍
数を得ることができなくなる場合が存在するということ
である。
When such a state occurs, when the signals M5 to M7, which are the original pulse waves, are obtained again, the value representative of the pulse wave period calculated for calculating the pulse rate becomes the true pulse wave. In the example shown in FIG. 4, the pulse wave of M6 may be ignored in practice because the value is far from the wave period. At this time, an accurate pulse rate cannot be obtained. In other words, in estimating the pulse rate during the period when a signal that can be regarded as a pulse wave signal could not be obtained,
Since the pulse wave cycle representative value itself, which is the basis for the estimation, has been affected by noise, there is a case where an accurate pulse rate cannot be obtained.

【0009】本発明はこのような点に鑑みなされたもの
であって、その目的とするところはノイズの影響を受け
ることなく脈波周期代表値の算出を行うことができ、こ
のために正確な脈拍数を算出することができる脈拍数の
算出方法を提供することにある。
SUMMARY OF THE INVENTION The present invention has been made in view of the above circumstances, and a purpose thereof is to calculate a pulse wave cycle representative value without being affected by noise. An object of the present invention is to provide a pulse rate calculation method capable of calculating a pulse rate.

【0010】[0010]

【課題を解決するための手段】しかして本発明は、脈波
を検出する脈拍センサの出力から脈波周期を求めて脈拍
数を算出するに際して、一時的に脈波とみなせる信号が
ないと判断された期間の脈拍数を、それ以前に求めてお
いた脈波周期から算出した脈波周期代表値を利用して推
定する脈拍数の算出方法において、予め定めた期間の脈
波信号の自己相関関数を算出し、その一番目のプラス方
向のピーク値が予め定めた閾値より大きい場合に、上記
期間内の脈波周期をもって前記脈波周期代表値とするこ
とに第1の特徴を有しており、予め定めた期間の脈波信
号の自己相関関数を算出し、その一番目のマイナス方向
のピーク値が予め定めた閾値より小さく且つその一番目
のプラス方向のピーク値が予め定めた閾値より大きい場
合に、上記期間内の脈波周期をもって前記脈波周期代表
値とすることに第2の特徴を有している。
According to the present invention, when a pulse rate is calculated by obtaining a pulse wave period from an output of a pulse sensor for detecting a pulse wave, it is determined that there is no signal which can be temporarily regarded as a pulse wave. The pulse rate during the specified period is estimated using the pulse wave cycle representative value calculated from the pulse wave cycle determined before that, and the autocorrelation of the pulse wave signal during the predetermined period is calculated. A first feature is that a function is calculated, and when the first peak value in the positive direction is larger than a predetermined threshold value, the pulse wave cycle within the period is used as the pulse wave cycle representative value. The autocorrelation function of the pulse wave signal for a predetermined period is calculated, and the first peak value in the minus direction is smaller than a predetermined threshold value, and the first peak value in the positive direction is smaller than a predetermined threshold value. If greater, within the above period With the pulse wave period and has a second feature in that with the pulse wave period representative value.

【0011】自己相関関数を利用して脈波周期代表値を
算出するために、ノイズの影響を防ぐことができるもの
であり、誤った脈波周期代表値を基準として誤った脈拍
数を算出してしまう事態が生じるのを防ぐことができ
る。
Since the pulse wave cycle representative value is calculated using the autocorrelation function, the influence of noise can be prevented, and an erroneous pulse rate is calculated based on an erroneous pulse wave cycle representative value. Can be prevented from occurring.

【0012】[0012]

【発明の実施の形態】以下本発明を実施の形態の一例に
基づいて詳述する。なお、脈拍数の測定のための器材や
基本となる脈拍数の算出方法に関しては、前記公報に記
載のものと基本的に同じであることから説明は省略し、
一時的に脈波とみなせる信号がないと判断された期間の
脈拍数の算出根拠となる脈波周期代表値の算出方法につ
いてのみ説明すると、図1に示すように、まず脈波周期
代表値の初期値を予め定めておく。この初期値には、例
えば1分間の脈拍数が70拍と予想されるので、これに
基づいた値を設定する。そして、予め定めた時間分の脈
波信号データを読み込み、自己相関関数を算出し、次に
脈波周期代表値を基準として必要に応じて脈波周期の補
間を行う。さらに、その自己相関関数のプラス側の一番
目のピーク値(図2のロ参照)が予め定めた値より大の
場合、脈波信号にノイズは混入していないと判断し、そ
の期間の脈波信号の脈波周期を代表値として採用して脈
波周期の代表値の書き換えを行い、次の脈波信号データ
の読み込みを行う。
DESCRIPTION OF THE PREFERRED EMBODIMENTS The present invention will be described below in detail based on an embodiment. The equipment for measuring the pulse rate and the method of calculating the basic pulse rate are not described because they are basically the same as those described in the above publication,
Only a method of calculating a pulse wave cycle representative value serving as a basis for calculating a pulse rate during a period in which it is determined that there is no signal that can be temporarily regarded as a pulse wave will be described. First, as shown in FIG. Initial values are determined in advance. As the initial value, for example, the pulse rate per minute is expected to be 70, so a value based on this is set. Then, the pulse wave signal data for a predetermined time is read, an autocorrelation function is calculated, and then the pulse wave cycle is interpolated as necessary based on the pulse wave cycle representative value. Further, when the first peak value on the plus side of the autocorrelation function (see FIG. 2B) is larger than a predetermined value, it is determined that noise is not mixed in the pulse wave signal, and the pulse wave signal during that period is determined. The representative value of the pulse wave cycle is rewritten using the pulse wave cycle of the wave signal as a representative value, and the next pulse wave signal data is read.

【0013】なお、1回に読み込む脈波信号データは自
己相関関数を算出して脈波信号の1周期分周辺にピーク
がたつのに十分な長さ、例えば3秒以上とする。さら
に、自己相関関数のピーク値を算出するために、自己相
関係数を十分に細かい時間間隔で算出しておく。
The pulse wave signal data to be read at one time is calculated to have an autocorrelation function and has a length sufficient for a peak around one period of the pulse wave signal, for example, 3 seconds or more. Furthermore, in order to calculate the peak value of the autocorrelation function, the autocorrelation coefficient is calculated at sufficiently small time intervals.

【0014】図3に他例を示す。基本的には上記のもの
と同じであるが、自己相関関数から脈波周期代表値の評
価を行うにあたり、そのプラス側の一番目のピーク値
(図2のロ参照)が予め定められた値より大でありさら
にその自己相関関数のマイナス側の一番目のピーク値
(図2のイ参照)が予め定められた値より小である場合
に、脈波信号にはノイズは混入していないと判断し、そ
の期間の脈波信号の脈波周期を代表値として採用して脈
波周期の代表値の書き換えを行い、次の脈波信号データ
の読み込みを行う。この場合、図1に示したものに比し
て、ノイズの混入のおそれをより確実に排除することが
できる。
FIG. 3 shows another example. Basically the same as the above, but when evaluating the representative value of the pulse wave period from the autocorrelation function, the first peak value on the plus side (see FIG. 2B) is a predetermined value. If it is larger and the first peak value on the minus side of the autocorrelation function (see a in FIG. 2) is smaller than a predetermined value, it is determined that the pulse wave signal contains no noise. Judgment is made, the pulse wave cycle of the pulse wave signal in that period is adopted as a representative value, the representative value of the pulse wave cycle is rewritten, and the next pulse wave signal data is read. In this case, the possibility of noise mixing can be more reliably eliminated than that shown in FIG.

【0015】[0015]

【発明の効果】以上のように本発明においては、一時的
に脈波とみなせる信号が得られなかった期間の脈拍数を
予め求めた脈波周期代表値を基に算出する算出方法にお
いて、自己相関関数を用いて脈波周期代表値を決めてい
ることから、ノイズの影響を受けていない脈波周期代表
値を得ることができるものであり、このために正確な脈
拍数を求めることができる。
As described above, according to the present invention, in the calculation method for calculating the pulse rate in a period in which a signal that can be temporarily regarded as a pulse wave was not obtained based on the pulse wave cycle representative value obtained in advance, Since the pulse wave cycle representative value is determined using the correlation function, a pulse wave cycle representative value that is not affected by noise can be obtained, and therefore, an accurate pulse rate can be obtained. .

【0016】また、自己相関関数を用いるにあたり、プ
ラス側のピーク値だけでなくマイナス側のピーク値も参
照して評価を行うことで、脈波周期代表値をより正確な
ものとして算出することができ、これに伴ってより正確
な脈拍数の算出を行うことができる。
In using the autocorrelation function, the evaluation is performed with reference to not only the peak value on the plus side but also the peak value on the minus side, so that the pulse wave cycle representative value can be calculated as more accurate. Accordingly, a more accurate calculation of the pulse rate can be performed.

【図面の簡単な説明】[Brief description of the drawings]

【図1】本発明の実施の形態の一例における要部のフロ
ーチャートである。
FIG. 1 is a flowchart of a main part according to an example of an embodiment of the present invention.

【図2】自己相関関数におけるピーク値の説明図であ
る。
FIG. 2 is an explanatory diagram of a peak value in an autocorrelation function.

【図3】他例における要部のフローチャートである。FIG. 3 is a flowchart of a main part in another example.

【図4】ノイズの混入している脈波信号の一例を示すタ
イムチャートである。
FIG. 4 is a time chart showing an example of a pulse wave signal mixed with noise.

【図5】ノイズの混ざっていない脈波信号の一例を示す
タイムチャートである。
FIG. 5 is a time chart showing an example of a pulse wave signal in which noise is not mixed.

───────────────────────────────────────────────────── フロントページの続き (72)発明者 道盛 章弘 大阪府門真市大字門真1048番地松下電工株 式会社内 (72)発明者 萩原 啓 大阪府門真市大字門真1048番地松下電工株 式会社内 Fターム(参考) 4C017 AA02 AA10 BC16 FF05  ──────────────────────────────────────────────────続 き Continuing on the front page (72) Inventor Akihiro Michimori 1048 Kadoma Kadoma, Osaka Prefecture Matsushita Electric Works, Ltd. F term (reference) 4C017 AA02 AA10 BC16 FF05

Claims (2)

【特許請求の範囲】[Claims] 【請求項1】 脈波を検出する脈拍センサの出力から脈
波周期を求めて脈拍数を算出するに際して、一時的に脈
波とみなせる信号がないと判断された期間の脈拍数を、
それ以前に求めておいた脈波周期から算出した脈波周期
代表値を利用して推定する脈拍数の算出方法において、
予め定めた期間の脈波信号の自己相関関数を算出し、そ
の一番目のプラス方向のピーク値が予め定めた閾値より
大きい場合に、上記期間内の脈波周期をもって前記脈波
周期代表値とすることを特徴とする脈拍数の算出方法。
When calculating a pulse rate by obtaining a pulse cycle from an output of a pulse sensor for detecting a pulse wave, the pulse rate during a period in which it is determined that there is no signal that can be temporarily regarded as a pulse wave is calculated by:
In a method of calculating a pulse rate to be estimated using a pulse wave cycle representative value calculated from a pulse wave cycle obtained before that,
Calculate the autocorrelation function of the pulse wave signal for a predetermined period, and when the first positive peak value is larger than a predetermined threshold, the pulse wave period representative value with the pulse wave period within the period. A method of calculating a pulse rate.
【請求項2】 脈波を検出する脈拍センサの出力から脈
波周期を求めて脈拍数を算出するに際して、一時的に脈
波とみなせる信号がないと判断された期間の脈拍数を、
それ以前に求めておいた脈波周期から算出した脈波周期
代表値を利用して推定する脈拍数の算出方法において、
予め定めた期間の脈波信号の自己相関関数を算出し、そ
の一番目のマイナス方向のピーク値が予め定めた閾値よ
り小さく且つその一番目のプラス方向のピーク値が予め
定めた閾値より大きい場合に、上記期間内の脈波周期を
もって前記脈波周期代表値とすることを特徴とする脈拍
数の算出方法。
2. When calculating a pulse rate by obtaining a pulse cycle from an output of a pulse sensor for detecting a pulse wave, a pulse rate during a period in which it is determined that there is no signal that can be temporarily regarded as a pulse wave is calculated by:
In a method of calculating a pulse rate to be estimated using a pulse wave cycle representative value calculated from a pulse wave cycle obtained before that,
Calculate the autocorrelation function of the pulse wave signal for a predetermined period, and when the first peak value in the negative direction is smaller than the predetermined threshold value and the first peak value in the positive direction is larger than the predetermined threshold value A pulse wave period within the period, the pulse wave period representative value being used as the pulse wave period representative value.
JP2000157258A 2000-05-26 2000-05-26 Method of computing pulse rate Withdrawn JP2001333890A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2000157258A JP2001333890A (en) 2000-05-26 2000-05-26 Method of computing pulse rate

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2000157258A JP2001333890A (en) 2000-05-26 2000-05-26 Method of computing pulse rate

Publications (1)

Publication Number Publication Date
JP2001333890A true JP2001333890A (en) 2001-12-04

Family

ID=18661919

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2000157258A Withdrawn JP2001333890A (en) 2000-05-26 2000-05-26 Method of computing pulse rate

Country Status (1)

Country Link
JP (1) JP2001333890A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012130391A (en) * 2010-12-20 2012-07-12 Oki Electric Industry Co Ltd Heart rate detection device, method, and program
CN104706337A (en) * 2015-02-11 2015-06-17 华东师范大学 Automatic pulse wave crest and trough detection method
CN111803041A (en) * 2020-06-19 2020-10-23 浙江澍源智能技术有限公司 Human body pulse peak prediction method and device

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012130391A (en) * 2010-12-20 2012-07-12 Oki Electric Industry Co Ltd Heart rate detection device, method, and program
CN104706337A (en) * 2015-02-11 2015-06-17 华东师范大学 Automatic pulse wave crest and trough detection method
CN111803041A (en) * 2020-06-19 2020-10-23 浙江澍源智能技术有限公司 Human body pulse peak prediction method and device
CN111803041B (en) * 2020-06-19 2024-05-24 浙江澍源智能技术有限公司 Human pulse peak prediction method and device

Similar Documents

Publication Publication Date Title
JPH06148003A (en) Ultrasonic temperature measuring equipment
JP3432210B2 (en) Flow measurement device
EP0341445B1 (en) Method of and apparatus for measuring revolution speed
US7805277B2 (en) Step number measuring apparatus
KR101276973B1 (en) Pulse frequency measurement method and apparatus
JPH0240173B2 (en)
JP2001333890A (en) Method of computing pulse rate
JP4476667B2 (en) Pedometer
JPS624971B2 (en)
EP0290352B1 (en) Electronic clinical thermometer
JPH09173308A (en) Correcting method for case temporarily incapable of getting pulse wave
JP4353682B2 (en) Flow rate measuring method and apparatus
US20020039027A1 (en) Dynamically periodical detecting method and detector
JP2699832B2 (en) Signal detection processing circuit
JPH08336502A (en) Method and device for heartbeat interval detection
JPH08233842A (en) Speed detector and detection of abnormality thereof
JP2008032461A (en) Apparatus and program for detecting state change
JP2005017235A (en) Pulse width measuring device and method
JPH0377067A (en) Rotating speed measuring instrument
JP3824043B2 (en) Encoder system
JPS59214920A (en) Calculating method of pulse frequency
JPH11118651A (en) Digital pressure gage
JP2569071Y2 (en) Frequency voltage converter
JPH021566A (en) Cycle measuring instrument
JPH11183554A (en) Fault point locating device

Legal Events

Date Code Title Description
A300 Withdrawal of application because of no request for examination

Free format text: JAPANESE INTERMEDIATE CODE: A300

Effective date: 20070807