JP2006238971A - Cycle calculator and program - Google Patents

Cycle calculator and program Download PDF

Info

Publication number
JP2006238971A
JP2006238971A JP2005055222A JP2005055222A JP2006238971A JP 2006238971 A JP2006238971 A JP 2006238971A JP 2005055222 A JP2005055222 A JP 2005055222A JP 2005055222 A JP2005055222 A JP 2005055222A JP 2006238971 A JP2006238971 A JP 2006238971A
Authority
JP
Japan
Prior art keywords
period
sensor
measured
reference period
pulse wave
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
Application number
JP2005055222A
Other languages
Japanese (ja)
Other versions
JP4326486B2 (en
Inventor
Yoshinori Tanabe
義憲 田邉
Hirokazu Genno
広和 源野
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.)
Sanyo Electric Co Ltd
Original Assignee
Sanyo Electric Co 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 Sanyo Electric Co Ltd filed Critical Sanyo Electric Co Ltd
Priority to JP2005055222A priority Critical patent/JP4326486B2/en
Publication of JP2006238971A publication Critical patent/JP2006238971A/en
Application granted granted Critical
Publication of JP4326486B2 publication Critical patent/JP4326486B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To calculate appropriately a standard cycle that corresponds to heart rate or pulse rate of a subject. <P>SOLUTION: A cycle calculator 1 is provided which is arranged on a subject detachably. A pulse wave sensor 10 measures a standard cycle, and an acceleration sensor 90 measures kinetic amount of the subject that correlates with the standard cycle. A CPU 70 corrects the standard cycle measured by the pulse wave sensor 10 using the kinetic amount measured by the acceleration sensor. <P>COPYRIGHT: (C)2006,JPO&NCIPI

Description

本発明は、基準周期を算出する周期算出装置及びプログラムに関する。   The present invention relates to a period calculation device and a program for calculating a reference period.

従来から、血管に光を照射し、当該血管の拍動に応じた反射光又は透過光の光量変化に基づいて脈波を特定し、特定した脈波の周期(以下では単に脈波周期と称する)を算出する周期算出装置などが提案されている。この周期算出装置では、脈波周期がn個測定され、測定された脈波周期が大きさ順に並べられて、その並べられた順で中央にくる脈波周期(以下では中央周期と称する)が用いられることにより、脈拍数が算出されている(例えば、特許文献1参照)。これにより、測定されたn個の脈波周期の中からノイズによる影響の少ない脈波周期が用いられるため、当該脈波周期に対応する脈拍数が高精度に算出されるように見える。
特開2002−28138号公報
Conventionally, a blood vessel is irradiated with light, a pulse wave is identified based on a change in the amount of reflected or transmitted light according to the pulsation of the blood vessel, and the cycle of the identified pulse wave (hereinafter simply referred to as a pulse wave cycle). ) Has been proposed. In this cycle calculation device, n pulse wave cycles are measured, the measured pulse wave cycles are arranged in order of magnitude, and a pulse wave cycle (hereinafter referred to as a central cycle) that comes to the center in the arranged order. The pulse rate is calculated by being used (for example, refer patent document 1). As a result, a pulse wave period that is less influenced by noise is used among the n pulse wave periods measured, and thus it appears that the pulse rate corresponding to the pulse wave period is calculated with high accuracy.
JP 2002-28138 A

しかしながら、上記周期算出装置では、本来の脈波周期の間に1又は複数ノイズ(例えば、図13に示すP2−P4の脈波周期の間であるP3)が含まれることにより、当該本来の脈波周期(例えば、図13に示すP2−P4の脈波周期)が分断された場合には、当該本来の脈波周期が上記中央周期にはなり得ず、当該本来の脈波周期が脈拍数の算出に用いられなくなるため、高精度に脈拍数が算出され難かった。   However, in the period calculation device, one or a plurality of noises (for example, P3 which is between P2 and P4 pulse wave periods shown in FIG. 13) are included in the original pulse wave period. When the wave period (for example, P2-P4 pulse wave period shown in FIG. 13) is divided, the original pulse wave period cannot be the central period, and the original pulse wave period is the pulse rate. It is difficult to calculate the pulse rate with high accuracy.

また、脈波がフィルタ回路を通過すると、一部の脈波(例えば、図13に示すP4-5の脈波,点線部分の脈波)が排除されることがあるため、当該一部の脈波が脈拍数の算出に用いられず、高精度に脈拍数が算出され難かった。   Further, when the pulse wave passes through the filter circuit, a part of the pulse wave (for example, the pulse wave of P4-5 shown in FIG. 13 or the pulse wave of the dotted line portion) may be excluded. Waves were not used to calculate the pulse rate, and it was difficult to calculate the pulse rate with high accuracy.

一方、出願人らは、基準周期を設定し、該基準周期に対する上限値及び下限値の範囲内である特定範囲に、該基準周期の次の周期である次周期が属している場合には、該次周期を脈波周期として算出する装置を提案している。これにより、該次周期が当該特定範囲に属しない場合には、当該次周期が脈波周期として算出されないこととなる。   On the other hand, when the applicant sets a reference period and the next period, which is the next period of the reference period, belongs to the specific range that is within the range of the upper limit value and the lower limit value for the reference period, An apparatus for calculating the next period as a pulse wave period is proposed. Thereby, when the next cycle does not belong to the specific range, the next cycle is not calculated as a pulse wave cycle.

従って、上記次周期が上記特定範囲に属しない短い周期(例えば、図13に示すP2−P3;P3はノイズによる脈波)や、上記次周期が上記特定範囲を超える長い周期(例えば、図11に示すP4-5;本来のP4-5はフィルタ回路などにより除去されている)が脈波周期とされないため、適切な脈拍数が算出されるように見える。   Therefore, a short period (for example, P2-P3; P3 shown in FIG. 13 is a pulse wave due to noise) in which the next period does not belong to the specific range, or a long period (for example, FIG. 11 in which the next period exceeds the specific range). P4-5; original P4-5 is removed by a filter circuit or the like) is not set as a pulse wave period, so that it seems that an appropriate pulse rate is calculated.

ところが、上記装置は、上記基準周期に対する上限値及び下限値の範囲内である特定範囲に上記次周期が属しているか否かにより、当該次周期を脈波周期として算出している。このため、上記装置は、上記基準周期が適切でなければ上記次周期を脈波周期として適切に算出することができなかった。   However, the apparatus calculates the next cycle as a pulse wave cycle depending on whether or not the next cycle belongs to a specific range that is within a range of an upper limit value and a lower limit value with respect to the reference cycle. For this reason, the apparatus cannot appropriately calculate the next period as a pulse wave period unless the reference period is appropriate.

具体的には、上記装置は、フィルタ回路を備えているため、本来存在しなければならない脈波(例えば、図13に示すP4-5)を除去することがある。このため、上記装置は、本来の脈波(例えば、図13に示すP4-5)が除かれると、その本来の脈波が除かれた状態での各脈波間(例えば、図13に示すP4−P5間)の周期を基準周期としなければならず、その基準周期を用いることにより上記次周期を誤算出することがあった。   Specifically, since the apparatus includes a filter circuit, a pulse wave (for example, P4-5 shown in FIG. 13) that should originally exist may be removed. For this reason, when the original pulse wave (for example, P4-5 shown in FIG. 13) is removed, the above-described apparatus can detect the interval between the pulse waves in the state where the original pulse wave is removed (for example, P4 shown in FIG. 13). -P5) must be a reference period, and the next period may be erroneously calculated by using the reference period.

これにより、従来から脈拍数が算出される基となる上記基準周期を適切にすることのできる技術の開発が望まれていた。   Accordingly, it has been desired to develop a technique that can appropriately set the reference period that is a basis for calculating the pulse rate.

そこで、本発明は上記点に鑑みてなされたものであり、基準周期を適切に算出することができる周期算出装置及びプログラムを提供することを目的とする。   Therefore, the present invention has been made in view of the above points, and an object thereof is to provide a period calculation device and a program that can appropriately calculate a reference period.

本発明は、上記課題を解決するために、被験者に着脱可能に装着される周期算出装置であって、第1センサにより測定される基準周期と相関のある被験者の運動量を測定する第2センサと、第2センサにより測定された運動量の値を用いて、第1センサにより測定された基準周期を補正する補正部とを備えることを特徴とする。   In order to solve the above-mentioned problem, the present invention is a period calculation device that is detachably attached to a subject, and a second sensor that measures a subject's momentum correlated with a reference period measured by the first sensor; And a correction unit that corrects the reference period measured by the first sensor using the value of the momentum measured by the second sensor.

かかる本発明の特徴によれば、測定された基準周期が本来の基準周期(以下では、本来基準周期と称する)とは大きくずれていても、第2センサにより測定された運動量の値を用いて、測定された基準周期を本来基準周期に近づけることができる。   According to such a feature of the present invention, even if the measured reference period is greatly deviated from the original reference period (hereinafter referred to as the original reference period), the momentum value measured by the second sensor is used. The measured reference period can be brought close to the original reference period.

上記発明においては、補正値は、第1センサにより測定された基準周期が所定範囲内に属しない場合には、第2センサにより測定された運動量の値を用いて基準周期を補正してもよい。   In the above invention, the correction value may be corrected using the momentum value measured by the second sensor when the reference period measured by the first sensor does not fall within the predetermined range. .

上記発明においては、第1センサは被験者の脈波の振幅の1つのピークから次のピークまでの間隔である基準周期を測定するものであり、基準周期に対する下限及び上限の間である特定範囲を設定する範囲設定手段と、第1センサにより測定された基準周期に対して次の周期である次周期が特定範囲に属する場合には、次周期を脈波周期として算出する周期算出部が備えられてもよい。   In the said invention, a 1st sensor measures the reference | standard period which is the space | interval from the one peak of the amplitude of a test subject's pulse wave to the next peak, and the specific range which is between the minimum and the upper limit with respect to a reference | standard period. A range setting means for setting and a period calculation unit for calculating the next period as a pulse wave period when the next period, which is the next period relative to the reference period measured by the first sensor, belongs to the specific range. May be.

上記発明においては、基準周期は被験者の脈拍数又は心拍数に対応する周期であり、運動量は互いに直交する3軸方向に作用する加速度の合成値であり、補正部は、多数の被験者の前記加速度の合成値と該多数の被験者の脈拍数又は心拍数との関係により算出された近似式と、前記第2センサにより測定された前記加速度の合成値とを用いて、前記第1センサにより測定された前記基準周期を補正してもよい。   In the above invention, the reference period is the pulse rate of the subject or a period corresponding to the heart rate, the momentum is a composite value of accelerations acting in three axis directions orthogonal to each other, and the correction unit is configured to measure the accelerations of a large number of subjects. Using the approximate expression calculated from the relationship between the combined value of the multiple and the pulse rate or heart rate of the large number of subjects, and the combined value of the acceleration measured by the second sensor. The reference period may be corrected.

本発明の特徴によれば、基準周期を適切に算出することができる周期算出装置及びプログラムを提供することができる。   According to the characteristics of the present invention, it is possible to provide a period calculation device and a program that can appropriately calculate a reference period.

[第1実施形態]
(周期算出装置の基本構成)
本発明における周期算出装置1について図面を参照しながら説明する。図1は、本実施形態における周期算出装置1の内部ブロックを示す図である。本実施形態における周期算出装置1は、脈波センサ10と、検出部20と、バンドパスフィルタ部30と、報知部40と、RAM50と、ROM60と、CPU70とを備えている。
[First Embodiment]
(Basic configuration of period calculation device)
The period calculation device 1 according to the present invention will be described with reference to the drawings. FIG. 1 is a diagram showing an internal block of a period calculation device 1 in the present embodiment. The period calculation device 1 in the present embodiment includes a pulse wave sensor 10, a detection unit 20, a bandpass filter unit 30, a notification unit 40, a RAM 50, a ROM 60, and a CPU 70.

脈波センサ10は、光を照射すると共に、当該光による皮膚などからの反射光を受光して受光量に応じた脈波を測定する。検出部20は、脈波センサ10により測定された脈波を検出する。バンドパスフィルタ部30は、検出部20により検出された脈波のうちの特定の周波数成分の脈波のみを抽出するものである。   The pulse wave sensor 10 irradiates light, receives light reflected from the skin or the like by the light, and measures a pulse wave corresponding to the amount of received light. The detection unit 20 detects the pulse wave measured by the pulse wave sensor 10. The band pass filter unit 30 extracts only a pulse wave having a specific frequency component from the pulse wave detected by the detection unit 20.

報知部40は、後述するCPU70により脈波周期が算出されない場合にはその旨を報知する。また、報知部40は、脈波周期、脈波周期に対応する脈拍数等を出力するものであり、例えば、液晶表示装置、スピーカなどから構成されている。RAM50は、一時的にデータを記憶するものである。ROM60は、周期算出装置1を動作させるためのプログラムなどを記憶するものである。   If the pulse wave cycle is not calculated by the CPU 70 described later, the notification unit 40 notifies that fact. Moreover, the alerting | reporting part 40 outputs the pulse rate etc. corresponding to a pulse wave period and a pulse wave period, for example, is comprised from a liquid crystal display device, a speaker, etc. The RAM 50 temporarily stores data. The ROM 60 stores a program for operating the cycle calculation device 1 and the like.

前記CPU70は、ROM60に記憶されたプログラムに基づいて処理を実行するものであり、ピーク間隔特定部71と、次ピーク間隔特定部72と、範囲設定部73と、周期算出部74とを備えている。   The CPU 70 executes processing based on a program stored in the ROM 60, and includes a peak interval specifying unit 71, a next peak interval specifying unit 72, a range setting unit 73, and a period calculating unit 74. Yes.

ピーク間隔特定部71は、脈波の振幅の1つのピークから次のピークまでの間隔をピーク間隔(後述する図3に示すt1,t2…;基準周期の次の周期である次周期)として特定する。次ピーク間隔特定部72は、ピーク間隔特定部71により特定されたピーク間隔(例えば、後述する図3(c)に示すt1)とそれに続くピーク間隔である連続ピーク間隔(例えば、後述する図3(c)に示すt2)とを1つのピーク間隔(例えば、後述する図3(c)に示すt1+t2、基準周期の次の周期である次周期)として特定する。   The peak interval specifying unit 71 specifies an interval from one peak to the next peak of the amplitude of the pulse wave as a peak interval (t1, t2,... Shown in FIG. 3 to be described later; the next cycle that is the next cycle of the reference cycle). To do. The next peak interval specifying unit 72 is a peak interval specified by the peak interval specifying unit 71 (for example, t1 shown in FIG. 3C described later) and a continuous peak interval (for example, FIG. (T2) shown in (c) is specified as one peak interval (for example, t1 + t2 shown in FIG. 3 (c) described later, the next cycle that is the next cycle of the reference cycle).

範囲設定部73は、基準周期T(後述する図3(a)参照)に対する下限及び上限の間である特定範囲TA(後述する図3(a)参照)を設定する。この基準周期Tには、予め設定された脈波周期、直前のピーク間隔の脈波周期、各ピーク間隔の平均の脈波周期などが挙げられる。本実施形態における基準周期Tは、所定のピーク間隔に対して直前に算出された脈波周期として説明する。   The range setting unit 73 sets a specific range TA (see FIG. 3A described later) that is between a lower limit and an upper limit with respect to the reference period T (see FIG. 3A described later). Examples of the reference period T include a preset pulse wave period, a pulse wave period at the immediately preceding peak interval, an average pulse wave period at each peak interval, and the like. The reference period T in the present embodiment will be described as a pulse wave period calculated immediately before a predetermined peak interval.

周期算出部74は、ピーク間隔特定部71により特定されたピーク間隔(例えば、後述する図3(b)に示すt1)が特定範囲TAに属する場合には、当該ピーク間隔を脈波周期として算出する。また、周期算出部74は、次ピーク間隔特定部72により特定されたピーク間隔(例えば、後述する図3(c)に示すt1+t2)が特定範囲TAに属する場合には、当該ピーク間隔を脈波周期として算出する。   When the peak interval specified by the peak interval specifying unit 71 (for example, t1 shown in FIG. 3B described later) belongs to the specific range TA, the cycle calculating unit 74 calculates the peak interval as a pulse wave cycle. To do. Further, when the peak interval specified by the next peak interval specifying unit 72 (for example, t1 + t2 shown in FIG. 3C described later) belongs to the specific range TA, the period calculating unit 74 converts the peak interval into a pulse wave. Calculated as a period.

さらに、周期算出部74は、ピーク間隔特定部71又は次ピーク間隔特定部72により特定された複数のピーク間隔(例えば、後述する図3(d)に示すt1+t2,t1+t2+t3)が特定範囲TAに属する場合には、複数のピーク間隔のうちの基準周期Tに最も近いピーク間隔(例えば、t1+t2)を脈波周期として算出する。   Further, in the period calculation unit 74, a plurality of peak intervals (for example, t1 + t2, t1 + t2 + t3 shown in FIG. 3D described later) specified by the peak interval specifying unit 71 or the next peak interval specifying unit 72 belong to the specific range TA. In this case, the peak interval (for example, t1 + t2) closest to the reference cycle T among the plurality of peak intervals is calculated as the pulse wave cycle.

また、周期算出部74は、ピーク間隔特定部71又は次ピーク間隔特定部72により特定された複数のピーク間隔(例えば、後述する図4(c)に示すt1,t1+t2)が特定範囲TAに属しない場合には、複数のピーク間隔のうちの最初のピーク間隔(例えば、t1)の終点を基準周期Tの開始点として設定する(例えば、後述する図4(d)及び(g))。   In addition, the period calculation unit 74 includes a plurality of peak intervals specified by the peak interval specifying unit 71 or the next peak interval specifying unit 72 (for example, t1, t1 + t2 shown in FIG. 4C described later) belonging to the specific range TA. If not, the end point of the first peak interval (for example, t1) among the plurality of peak intervals is set as the start point of the reference period T (for example, FIGS. 4D and 4G described later).

(周期算出装置の動作)
以下において、周期算出装置の動作である周期算出方法について、図面を参照しながら説明する。図2は、本実施形態における周期算出方法を示す図である。図2に示すように、ステップ1において、検出部20は、脈波センサ10から出力された信号を脈波として検出する。その後、バンドパスフィルタ部30は、検出部20により検出された脈波のうちの特定の周波数成分のみを抽出する。
(Operation of period calculation device)
Hereinafter, a period calculation method that is an operation of the period calculation device will be described with reference to the drawings. FIG. 2 is a diagram illustrating a period calculation method in the present embodiment. As shown in FIG. 2, in step 1, the detection unit 20 detects a signal output from the pulse wave sensor 10 as a pulse wave. Thereafter, the band pass filter unit 30 extracts only a specific frequency component from the pulse wave detected by the detection unit 20.

ステップ2において、CPU70は、脈波の振幅の1つのピークから次のピークまでの間隔をピーク間隔(後述する図3に示すt1,t2…)として特定する。ステップ3において、CPU70は、特定したピーク間隔を用いて周期算出方法を行う。ステップ4において、CPU70は、脈波周期が算出されたか否か確認する。また、CPU70は、脈波周期が算出された場合にはステップ5の処理に移り、脈波周期が算出されていない場合にはステップ6の処理に移る。   In step 2, the CPU 70 specifies an interval from one peak of the amplitude of the pulse wave to the next peak as a peak interval (t1, t2,... Shown in FIG. 3 described later). In step 3, the CPU 70 performs a period calculation method using the specified peak interval. In step 4, the CPU 70 confirms whether or not the pulse wave cycle has been calculated. Further, the CPU 70 proceeds to the process of step 5 when the pulse wave period is calculated, and proceeds to the process of step 6 when the pulse wave period is not calculated.

ステップ5において、CPU70は、算出された脈波周期に対応する脈拍数を出力するように指示するコマンドを報知部40に出力する。この報知部40は、CPU70により出力されたコマンドに基づいて脈拍数などを表示する。ステップ6において、CPU70は、脈波周期がn回以上連続して算出されていない場合にはステップ7の処理に移り、脈波周期がn回以内算出されていない場合にはステップ1の処理に戻る。   In step 5, the CPU 70 outputs a command for instructing to output the pulse rate corresponding to the calculated pulse wave cycle to the notification unit 40. The notification unit 40 displays the pulse rate and the like based on the command output by the CPU 70. In step 6, the CPU 70 proceeds to the process of step 7 when the pulse wave cycle has not been calculated continuously n times or more, and proceeds to the process of step 1 when the pulse wave cycle has not been calculated within n times. Return.

ステップ7において、CPU70は、処理がエラーであることを報知するように指示するコマンドを報知部40に出力する。この報知部40は、CPU70により出力されたコマンドに基づいて特定情報(例えば、脈波センサ10を適切な位置に付けて下さいなどの文字)を表示する。   In step 7, the CPU 70 outputs to the notification unit 40 a command for instructing that the process is an error. The notification unit 40 displays specific information (for example, characters such as “please attach the pulse wave sensor 10 to an appropriate position”) based on the command output by the CPU 70.

ステップ8において、CPU70は、終了条件(例えば、脈波周期の算出を終了させるための操作が受け付けられた場合など)が成立したか否か確認する。また、CPU70は、終了条件が成立した場合には本処理を終了し、終了条件が成立していない場合にはステップ1の処理に戻る。   In step 8, the CPU 70 confirms whether an end condition (for example, an operation for ending the calculation of the pulse wave period is accepted) is satisfied. Further, the CPU 70 ends the process when the end condition is satisfied, and returns to the process of step 1 when the end condition is not satisfied.

次に、上述したステップ3の周期算出処理について説明する。図3及び図4は、ステップ3の周期算出処理により脈波周期が算出されるまでの様子を示す図である。この周期算出処理はパターン1及びパターン2に大きく分けられる。以下ではパターン1及びパターン2について詳細に説明する。   Next, the period calculation process in step 3 described above will be described. 3 and 4 are diagrams showing a state until the pulse wave cycle is calculated by the cycle calculation processing in step 3. This period calculation process is roughly divided into pattern 1 and pattern 2. Hereinafter, pattern 1 and pattern 2 will be described in detail.

(A)パターン1について
図3(a)に示すように、先ず、範囲設定部73は、基準周期である前回の脈波周期T及びその脈波周期Tに対する下限Tmin及び上限Tmaxの間である特定範囲TAを設定する。そして、最初のピーク間隔t1が特定範囲TAに属する場合には、周期算出部74は、ピーク間隔t1を脈波周期として算出する。
(A) Pattern 1 As shown in FIG. 3A, first, the range setting unit 73 is between the previous pulse wave period T, which is the reference period, and the lower limit Tmin and the upper limit Tmax with respect to the pulse wave period T. A specific range TA is set. When the first peak interval t1 belongs to the specific range TA, the cycle calculating unit 74 calculates the peak interval t1 as a pulse wave cycle.

一方、図3(b)に示すように、最初のピーク間隔t1(又はt1’)が特定範囲TAに属しない場合には、周期算出部74は、ピーク間隔t1(又はt1’)を脈波周期として算出しない。   On the other hand, as shown in FIG. 3B, when the first peak interval t1 (or t1 ′) does not belong to the specific range TA, the period calculation unit 74 uses the peak interval t1 (or t1 ′) as a pulse wave. Not calculated as a period.

この場合には、図3(c)に示すように、次ピーク間隔特定部72は、ピーク間隔t1とそれに続くピーク間隔である連続ピーク間隔t2とをピーク間隔(t1+t2)として特定する。そして、図3(c)に示すように、ピーク間隔(t1+t2)が特定範囲TAに属する場合には、周期算出部74は、ピーク間隔(t1+t2)を脈波周期として算出する。   In this case, as shown in FIG. 3C, the next peak interval specifying unit 72 specifies the peak interval t1 and the continuous peak interval t2, which is the subsequent peak interval, as the peak interval (t1 + t2). Then, as shown in FIG. 3C, when the peak interval (t1 + t2) belongs to the specific range TA, the cycle calculation unit 74 calculates the peak interval (t1 + t2) as a pulse wave cycle.

但し、図3(d)に示すように、ピーク間隔(t1+t2)のみならず、ピーク間隔(t1+t2)とそれに続く連続ピーク間隔t3との総和であるピーク間隔(t1+t2+t3)も特定範囲TAに属する場合には、周期算出部74は、最も基準周期Tに近いピーク間隔(t1+t2)を脈波周期として算出する。   However, as shown in FIG. 3D, not only the peak interval (t1 + t2) but also the peak interval (t1 + t2 + t3) which is the sum of the peak interval (t1 + t2) and the subsequent continuous peak interval t3 belongs to the specific range TA. The period calculation unit 74 calculates the peak interval (t1 + t2) closest to the reference period T as the pulse wave period.

(B)パターン2について
図4(a)に示すように、先ず、範囲設定部73は、基準周期である前回の脈波周期T及びその脈波周期Tに対する下限Tmin及び上限Tmaxの間である特定範囲TAを設定する。そして、図4(b)及び(c)に示すように、最初のピーク間隔t1が特定範囲TAに属しない場合には、次ピーク間隔特定部72は、ピーク間隔t1とそれに続くピーク間隔である連続ピーク間隔t2とをピーク間隔(t1+t2)として特定する。
(B) Pattern 2 As shown in FIG. 4A, first, the range setting unit 73 is between the previous pulse wave period T, which is the reference period, and the lower limit Tmin and the upper limit Tmax with respect to the pulse wave period T. A specific range TA is set. As shown in FIGS. 4B and 4C, when the first peak interval t1 does not belong to the specific range TA, the next peak interval specifying unit 72 is the peak interval t1 and the subsequent peak interval. The continuous peak interval t2 is specified as the peak interval (t1 + t2).

その後、図4(c)及び(d)に示すように、最初のピーク間隔t1のみならず、次に特定されたピーク間隔(t1+t2)も特定範囲TAに属しない場合には、周期算出部74は、最初のピーク間隔t1の終点を基準周期Tの開始点に変更する(図4(g)参照)。したがって、最初のピーク間隔t1が削除されることとなる。   Thereafter, as shown in FIGS. 4C and 4D, when not only the first peak interval t1 but also the next specified peak interval (t1 + t2) does not belong to the specific range TA, the period calculation unit 74 Changes the end point of the first peak interval t1 to the start point of the reference period T (see FIG. 4G). Therefore, the first peak interval t1 is deleted.

なお、本実施形態では、次に特定されたピーク間隔(t1+t2)が上限Tmaxを超える場合には、周期算出部74は、最初のピーク間隔t1の終点を基準周期Tの開始点に変更する(図4(d)参照)。   In the present embodiment, when the next specified peak interval (t1 + t2) exceeds the upper limit Tmax, the cycle calculating unit 74 changes the end point of the first peak interval t1 to the start point of the reference cycle T ( (Refer FIG.4 (d)).

そして、図4(d)に示すように、ピーク間隔t2が特定範囲TAに属しない場合には、次ピーク間隔特定部72は、ピーク間隔t2とそれに続くピーク間隔である連続ピーク間隔t3とをピーク間隔(t2+t3)として特定する。その後、図4(e)及び(f)に示すように、ピーク間隔(t2+t3)が特定範囲TAに属する場合には、周期算出部74は、ピーク間隔(t2+t3)を脈波周期として算出する。   As shown in FIG. 4D, when the peak interval t2 does not belong to the specific range TA, the next peak interval specifying unit 72 sets the peak interval t2 and the continuous peak interval t3 that is the subsequent peak interval. It is specified as the peak interval (t2 + t3). Thereafter, as shown in FIGS. 4E and 4F, when the peak interval (t2 + t3) belongs to the specific range TA, the period calculation unit 74 calculates the peak interval (t2 + t3) as a pulse wave period.

ここで、周期算出部74は、上記パターン1又は上記パターン2を用いて算出した複数の脈波周期の平均値(秒)を算出し、60をその平均値で除算することにより1分間における脈拍数(bpm)を算出する。   Here, the period calculation unit 74 calculates an average value (seconds) of a plurality of pulse wave periods calculated using the pattern 1 or the pattern 2, and divides 60 by the average value to thereby calculate the pulse in one minute. The number (bpm) is calculated.

(作用及び効果)
このような本発明の特徴によれば、ピーク間隔(例えば、図3に示すt1)とそれに続く連続ピーク間隔(例えば、図3に示すt2)との総和であるピーク間隔(t1+t2)が特定範囲内TAに属する場合には、当該ピーク間隔(t1+t2)を脈波周期として算出する。これにより、特定範囲内TAに属しない程短いピーク間隔、すなわち本来の脈波の1つのピークとノイズのピークとの間隔である可能性の高いピーク間隔(t1)を脈波周期として算出せずに、そのピーク間隔t1とそれに続く連続間隔t2との総和である本来のピーク間隔(t1+t2)を脈波周期として算出するため、脈波にノイズが含まれても適切な脈波周期を算出することができる。
(Function and effect)
According to such a feature of the present invention, the peak interval (t1 + t2), which is the sum of the peak interval (eg, t1 shown in FIG. 3) and the subsequent continuous peak interval (eg, t2 shown in FIG. 3), is in a specific range. When belonging to the inner TA, the peak interval (t1 + t2) is calculated as the pulse wave period. As a result, a peak interval that is so short that it does not belong to the specific range TA, that is, a peak interval (t1) that is likely to be the interval between one peak of the original pulse wave and the noise peak is not calculated as the pulse wave period. In addition, since the original peak interval (t1 + t2), which is the sum of the peak interval t1 and the subsequent continuous interval t2, is calculated as the pulse wave cycle, an appropriate pulse wave cycle is calculated even if noise is included in the pulse wave. be able to.

さらに、バンドパスフィルタ部30などにより排除された脈波の前後に存在する各ピークの間隔であるピーク間隔(例えば、図3に示すt1’)が本来のピーク間隔よりも広くなることがある。この場合には、当該広いピーク間隔(t1’)が特定範囲TAを超える場合には、その広いピーク間隔を脈波周期として算出しないため、一部の脈波が排除されても適切な脈波周期を算出することができる。   Furthermore, the peak interval (for example, t1 ′ shown in FIG. 3) that is the interval between the peaks existing before and after the pulse wave eliminated by the bandpass filter unit 30 or the like may be wider than the original peak interval. In this case, when the wide peak interval (t1 ′) exceeds the specific range TA, the wide peak interval is not calculated as the pulse wave period. The period can be calculated.

また、複数のピーク間隔(例えば、図3に示すt1+t2,t1+t2+t3)が特定範囲TAに属しても、当該複数のピーク間隔のうちの基準周期Tに最も近いピーク間隔(例えば、図3に示すt1+t2)を脈波周期として算出するため、より本来の周期に近いピーク間隔(t1+t2)を脈波周期として算出することができる。   Further, even if a plurality of peak intervals (for example, t1 + t2, t1 + t2 + t3 shown in FIG. 3) belong to the specific range TA, a peak interval closest to the reference period T (for example, t1 + t2 shown in FIG. 3) among the plurality of peak intervals. ) Is calculated as the pulse wave period, the peak interval (t1 + t2) closer to the original period can be calculated as the pulse wave period.

さらに、複数のピーク間隔(例えば、図4に示すt1,t1+t2)が特定範囲TAに属しない場合には、複数のピーク間隔のうちの最初のピーク間隔(例えば、図4に示すt1)の終点を基準周期Tの開始点とする。これにより、当該最初のピーク間隔(t1)を除いてその次のピーク間隔(t2)が特定範囲に属するか否かについて再び判定するため、より適切な脈波周期を算出することができる。   Furthermore, when a plurality of peak intervals (for example, t1, t1 + t2 shown in FIG. 4) do not belong to the specific range TA, the end points of the first peak interval (for example, t1 shown in FIG. 4) among the plurality of peak intervals. Is the starting point of the reference period T. Accordingly, since it is determined again whether or not the next peak interval (t2) belongs to the specific range except the first peak interval (t1), a more appropriate pulse wave cycle can be calculated.

また、脈波周期が算出されない場合にはその旨が報知されることにより、測定者は、脈波を測定するセンサの位置をより適切な位置に取付け直すことができる。   Further, when the pulse wave period is not calculated, the fact is notified so that the measurer can reattach the position of the sensor that measures the pulse wave to a more appropriate position.

[第2実施形態]
(周期算出装置の基本構成)
第1実施形態及び第2実施形態では、基準周期T(図3及び図4参照)は、前回の脈波周期Tであるという点で共通するが、第2実施形態では、基準周期Tは、さらに被験者に作用する運動量に応じて補正されるという点で異なる。以下では、第1実施形態と異なる点のみを説明し、共通する部分の説明については省略する。
[Second Embodiment]
(Basic configuration of period calculation device)
In the first embodiment and the second embodiment, the reference period T (see FIGS. 3 and 4) is common in that it is the previous pulse wave period T. In the second embodiment, the reference period T is Furthermore, it differs in that it is corrected according to the momentum acting on the subject. Below, only a different point from 1st Embodiment is demonstrated, and it abbreviate | omits about description of a common part.

図5は、本実施形態における周期算出装置1の内部ブロックを示す図である。なお、図5に示す周期算出装置1に含まれている各部は、図1に示す周期算出装置1に含まれている同名の各部と同様の機能を有するため、詳細な説明は省略する。   FIG. 5 is a diagram showing an internal block of the period calculation device 1 in the present embodiment. Note that each unit included in the period calculation device 1 illustrated in FIG. 5 has the same function as each unit of the same name included in the cycle calculation device 1 illustrated in FIG.

図5に示すように、周期算出装置1は、上記第1実施形態の構成に加えて、補正部75と、加速度センサ80と、検出部90とを備えている。   As shown in FIG. 5, the period calculation device 1 includes a correction unit 75, an acceleration sensor 80, and a detection unit 90 in addition to the configuration of the first embodiment.

加速度センサ80は、脈波センサ10、検出部20及びCPU70により測定される周期(基準周期T)と相関のある被験者の運動量を測定する第2センサを構成する。この運動量は、被験者に作用する加速度、被験者の歩数、被験者から吐き出されたガスの分析結果などでもよい。本実施形態における運動量は、被験者に作用する加速度であるものとする。なお、脈波センサ10、検出部20及びCPU70は第1センサを構成する。また、加速度センサ80は、ホルター心電計に備えられてもよい。   The acceleration sensor 80 constitutes a second sensor that measures the amount of exercise of the subject having a correlation with the period (reference period T) measured by the pulse wave sensor 10, the detection unit 20, and the CPU 70. This momentum may be the acceleration acting on the subject, the number of steps of the subject, the analysis result of the gas exhaled from the subject, and the like. The amount of exercise in the present embodiment is assumed to be acceleration acting on the subject. The pulse wave sensor 10, the detection unit 20, and the CPU 70 constitute a first sensor. The acceleration sensor 80 may be provided in a Holter electrocardiograph.

検出部90は、加速度センサ80により測定された加速度の値を検出する。   The detection unit 90 detects the value of acceleration measured by the acceleration sensor 80.

CPU70は、加速度センサ80により測定された加速度の値を用いて、脈波センサ10により測定された基準周期を補正する補正部75を備えている。例えば、補正部75は、多数の被験者の加速度の値と該多数の被験者の脈拍数又は心拍数との関係により算出された近似式と、加速度センサ80により測定された加速度の合成値とを用いて、脈波センサ10により測定された基準周期を補正する。   The CPU 70 includes a correction unit 75 that corrects the reference period measured by the pulse wave sensor 10 using the acceleration value measured by the acceleration sensor 80. For example, the correction unit 75 uses an approximate expression calculated from the relationship between the acceleration values of a large number of subjects and the pulse rate or heart rate of the large numbers of subjects, and a combined value of accelerations measured by the acceleration sensor 80. Thus, the reference period measured by the pulse wave sensor 10 is corrected.

また、補正値77は、脈波センサ10により測定された基準周期が所定範囲内に属しない場合には、加速度センサ80により測定された加速度の値を用いて基準周期を補正してもよい。   The correction value 77 may be corrected using the acceleration value measured by the acceleration sensor 80 when the reference period measured by the pulse wave sensor 10 does not fall within the predetermined range.

図6(a)は、本実施形態における脈波センサ10及び加速度センサ80が被験者に装着されている状態を示す図である。図6(b)は、被験者に作用する加速度の方向を示す図である。   FIG. 6A is a diagram illustrating a state where the pulse wave sensor 10 and the acceleration sensor 80 according to the present embodiment are attached to a subject. FIG. 6B is a diagram showing the direction of acceleration acting on the subject.

図6(a)に示すように、脈波センサ10は被験者の耳に装着され、加速度センサ80は被験者の胴体部分に装着される。なお、加速度センサ80は、脈波センサ10と共に被験者の耳などに装着されてもよい。   As shown in FIG. 6A, the pulse wave sensor 10 is attached to the ear of the subject, and the acceleration sensor 80 is attached to the trunk of the subject. The acceleration sensor 80 may be attached to the subject's ear together with the pulse wave sensor 10.

また、図6(b)に示すように、加速度センサ80が被験者に装着されることにより、加速度センサ80は、互いに直交する3軸の方向に加速度を測定する。すなわち、加速度センサ80は、被験者の進退方向に沿うx軸の加速度ix,被験者の左右方向に沿うy軸の加速度iy,被験者の鉛直方向に沿うz軸の加速度izを測定する。これらの加速度ix,加速度iy及び加速度izの合成値はitotalとする。   As shown in FIG. 6B, when the acceleration sensor 80 is attached to the subject, the acceleration sensor 80 measures acceleration in the directions of three axes orthogonal to each other. That is, the acceleration sensor 80 measures the x-axis acceleration ix along the subject's forward / backward direction, the y-axis acceleration iy along the subject's left-right direction, and the z-axis acceleration iz along the subject's vertical direction. The combined value of the acceleration ix, acceleration iy, and acceleration iz is itotal.

(周期算出装置の動作)
次に、周期算出装置の動作である基準周期算出処理1及び基準周期算出処理2について、図面を参照しながら説明する。以下では、基準周期算出処理1及び基準周期算出処理2の順に説明する。
(Operation of period calculation device)
Next, reference cycle calculation processing 1 and reference cycle calculation processing 2 which are operations of the cycle calculation apparatus will be described with reference to the drawings. Below, it demonstrates in order of the reference period calculation process 1 and the reference period calculation process 2. FIG.

(1)基準周期算出処理1
図7は、本実施形態における基準周期算出処理1を示す図である。図7に示すように、ステップ11において、加速度センサ80は、互いに直交する3軸に作用する被験者の加速度(ix,iy,iz)を測定する。
(1) Reference cycle calculation process 1
FIG. 7 is a diagram showing the reference cycle calculation process 1 in the present embodiment. As shown in FIG. 7, in step 11, the acceleration sensor 80 measures the acceleration (ix, iy, iz) of the subject acting on three axes orthogonal to each other.

ステップ12において、補正部75は、測定された加速度(ix,iy,iz)に基づいて、当該加速度(ix,iy,iz)の合成値である合成加速度の平均値を合成加速度
itotalとして算出する。すなわち、所定時間間隔における複数の合成加速度(ここではm個)が存在する場合に、その所定時間間隔における複数の合成加速度の平均値を合成加速度itotalとすると、合成加速度itotalは以下の式1で表現することができる。

Figure 2006238971
In step 12, the correction unit 75 calculates the average value of the combined acceleration, which is the combined value of the accelerations (ix, iy, iz), as the combined acceleration itotal based on the measured accelerations (ix, iy, iz). . That is, when there are a plurality of synthetic accelerations (m in this case) at a predetermined time interval, and the average value of the plurality of synthetic accelerations at the predetermined time interval is the synthetic acceleration itotal, the synthetic acceleration itotal is expressed by the following formula 1. Can be expressed.
Figure 2006238971

ステップ13において、補正部75は、予め設定された合成加速度−脈拍数特性(後述する式2の近似式)と、上記算出した合成加速度itotalとに基づいて、該合成加速度itotalに対する仮基準脈拍数を算出する。この仮基準脈拍数は、上記基準周期Tに対応する前回の脈拍数を最適値に補正するための仮の脈拍数である。この仮基準脈拍数は以下の式2の近似式で表現することができる。   In step 13, the correction unit 75 determines a provisional reference pulse rate for the combined acceleration itotal based on a previously set combined acceleration-pulse rate characteristic (approximate expression of equation 2 described later) and the calculated combined acceleration itotal. Is calculated. The temporary reference pulse rate is a temporary pulse rate for correcting the previous pulse rate corresponding to the reference period T to an optimum value. This provisional reference pulse rate can be expressed by an approximate expression of Expression 2 below.

仮基準脈拍数=a・itotal+b (a,b;定数)…式2
ステップ14において、補正部75は、前回の脈拍数と仮基準脈拍数とに基づいて当該前回の脈拍数を補正し、補正後の脈拍数を基準脈拍数とする。この基準脈拍数は以下の式3で表現することができる。
Temporary reference pulse rate = a · itotal + b (a, b; constant)... Formula 2
In step 14, the correction unit 75 corrects the previous pulse rate based on the previous pulse rate and the provisional reference pulse rate, and sets the corrected pulse rate as the reference pulse rate. This reference pulse rate can be expressed by Equation 3 below.

基準脈拍数=(A・前回の脈拍数+B・仮基準脈拍数)/(A+B) (A,B;定数)…式3
なお、上記Aは2、Bは1であることが好ましい。また、補正部75は、上記式3より前回の脈拍数と仮基準脈拍数とを按分することにより、基準脈拍数を算出しているが、これに限定されずに、重回帰分析、PLS解析を用いて、上記前回の脈拍数、上記合成加速度itotal、被験者の歩数、被験者により吐き出されたガス成分の分析結果に基づいて基準脈拍数を算出してもよい。
Reference pulse rate = (A · Previous pulse rate + B · Temporary reference pulse rate) / (A + B) (A, B; constant) ... Equation 3
In addition, it is preferable that A is 2 and B is 1. Further, the correction unit 75 calculates the reference pulse rate by dividing the previous pulse rate and the provisional reference pulse rate according to the above equation 3, but the present invention is not limited to this, and multiple regression analysis and PLS analysis are performed. May be used to calculate the reference pulse rate based on the previous pulse rate, the combined acceleration itotal, the number of steps of the subject, and the analysis result of the gas component exhaled by the subject.

ステップ15において、補正部75は、算出した基準脈拍数に基づいて基準周期Tを算出する。この基準周期Tは以下の式4で表現することができる。なお、基準脈拍数はbpmの単位であり、基準周期Tは秒の単位であるものとする。   In step 15, the correction unit 75 calculates the reference period T based on the calculated reference pulse rate. This reference period T can be expressed by the following Equation 4. It is assumed that the reference pulse rate is a unit of bpm and the reference period T is a unit of second.

基準周期T=60/基準脈拍数 …式4
ここで、第1実施形態における周期算出処理においてCPU70は、上記式4により算出された基準周期Tを前回の脈波周期(図3及び図4参照)として設定し、その脈波周期の次の脈波周期(次周期)を算出する。この設定された基準周期Tにより次の脈波周期が算出されるまでの過程は、上述した周期算出処理(図3及び図4参照)と同様であるため、詳細な説明は省略する。
Reference period T = 60 / reference pulse rate (Formula 4)
Here, in the cycle calculation process in the first embodiment, the CPU 70 sets the reference cycle T calculated by the above equation 4 as the previous pulse wave cycle (see FIG. 3 and FIG. 4), and next to the pulse wave cycle. The pulse wave cycle (next cycle) is calculated. Since the process until the next pulse wave period is calculated based on the set reference period T is the same as the above-described period calculation process (see FIGS. 3 and 4), detailed description thereof is omitted.

次に、図8を参照しながら、上記S13において用いられる合成加速度−脈拍数特性(近似式)を算出するまでの処理について説明する。図8は、合成加速度−脈拍数特性(近似式)を算出するまでの処理を示す図である。   Next, with reference to FIG. 8, a process until calculation of the combined acceleration-pulse rate characteristic (approximate expression) used in S13 will be described. FIG. 8 is a diagram illustrating a process until calculating a combined acceleration-pulse rate characteristic (approximate expression).

S21において、周期算出装置1は、所定数(n)の被験者に作用する加速度の合成値である合成加速度と、その合成加速度が発揮されている所定数(n)の被験者の脈拍数との入力を促す。なお、所定数(n)は、後述する最小二乗法により近似式が算出可能な数である。   In S <b> 21, the period calculation device 1 inputs a composite acceleration that is a composite value of accelerations acting on a predetermined number (n) of subjects and a pulse number of the predetermined number (n) of subjects who are exhibiting the combined acceleration. Prompt. The predetermined number (n) is a number for which an approximate expression can be calculated by the least square method described later.

S22において、周期算出装置1は、入力された所定数(n)の被験者の合成加速度とそれに対する所定数(n)の被験者の脈拍数との関係により、最小二乗法を用いて近似式(上記式2参照)を求める。ここで、図9は、所定数(n)の被験者に作用する合成加速度itotal−脈拍数特性との関係を示す図である。本実施形態では、仮基準脈拍数は次の近似式(上記式2と同様)により表現することができる。   In S22, the period calculation device 1 uses the least squares method to calculate an approximate expression (the above-described formula) based on the relationship between the input composite acceleration of the predetermined number (n) of subjects and the pulse number of the predetermined number (n) of subjects. Equation 2) is obtained. Here, FIG. 9 is a diagram showing the relationship between the combined acceleration itotal-pulse rate characteristics acting on a predetermined number (n) of subjects. In the present embodiment, the provisional reference pulse rate can be expressed by the following approximate expression (similar to Expression 2 above).

仮基準脈拍数=a・itotal+b
=0.0004・itotal+68.666
次に、図10を参照しながら、時間と基準周期T(図7に示すS15参照)に対する脈拍数との関係について説明する。図10は、第2実施形態の処理が実行された場合の時間に対する心拍数及び脈拍数の関係を示す図である。図10に示すように、脈拍数は、心拍数と略同様の値を示しており、心拍数から大きくずれていないことが分かる。
Temporary reference pulse rate = a · itotal + b
= 0.0004 · itotal + 68.666
Next, the relationship between time and the pulse rate with respect to the reference period T (see S15 shown in FIG. 7) will be described with reference to FIG. FIG. 10 is a diagram illustrating the relationship between the heart rate and the pulse rate with respect to time when the process of the second embodiment is executed. As shown in FIG. 10, the pulse rate shows substantially the same value as the heart rate, and it can be seen that it does not deviate significantly from the heart rate.

具体的には、500s付近では脈拍数が160であるが、その後の525s付近では脈拍数が125となり、脈拍数が急激に減少している。これは、525s付近において本来の脈波がフィルタ回路などにより除去されてしまい、脈拍数が誤算出されたためである。   Specifically, the pulse rate is 160 in the vicinity of 500 s, but the pulse rate is 125 in the vicinity of 525 s thereafter, and the pulse rate decreases rapidly. This is because the original pulse wave was removed by a filter circuit or the like in the vicinity of 525 s, and the pulse rate was erroneously calculated.

ところが、500s以降においても上記基準周期Tが上記合成加速度itotalなどにより適宜補正(上記式1乃至式4参照)される。このため、525s付近において脈拍数が誤算出されても、その後において適宜修正された基準周期Tにより適切な脈拍数が算出される。   However, even after 500 s, the reference period T is appropriately corrected (see the above formulas 1 to 4) by the combined acceleration itotal and the like. For this reason, even if the pulse rate is erroneously calculated in the vicinity of 525 s, an appropriate pulse rate is calculated based on the reference period T that is appropriately corrected thereafter.

ここで、図11は、第2実施形態の処理が実行されない場合の時間に対する心拍数及び脈拍数の関係を示す図である。図11に示すように、0〜485sにおいては脈拍数が心拍数と同様の数値を示すが、それ以降においては、フィルタ回路の影響により脈拍数が心拍数から大きくずれることがある。   Here, FIG. 11 is a diagram illustrating a relationship between the heart rate and the pulse rate with respect to time when the process of the second embodiment is not executed. As shown in FIG. 11, the pulse rate shows the same value as the heart rate at 0 to 485 s, but after that, the pulse rate may deviate greatly from the heart rate due to the influence of the filter circuit.

例えば、550s付近に対する脈拍数80の周期がその前の485s付近に対する脈拍数160の周期に対して2倍、3倍となることがある。つまり、フィルタ回路などの影響により550s付近において本来の脈波(例えば、図13に示すP4-5)が除去され、550s付近の周期がその前の480s付近の周期に対して2倍、3倍となることがある。 このため、550s付近の周期が基準周期Tとして誤算出されると、それ以降においては、誤算出された基準周期Tがそのまま用いられるため、その基準周期Tにより脈拍数も誤算出され続ける。本発明では、図10で説明したように、基準周期Tが誤算出されても、上記合成加速度itotalなどにより最適な基準周期Tに適宜補正(上記式1乃至式4参照)されるため、その最適な基準周期Tにより適切な脈拍数が算出される。   For example, the period of the pulse rate 80 for the vicinity of 550 s may be twice or three times the period of the pulse rate 160 for the vicinity of 485 s. That is, the original pulse wave (for example, P4-5 shown in FIG. 13) is removed in the vicinity of 550 s due to the influence of the filter circuit or the like, and the period in the vicinity of 550 s is doubled or tripled in comparison with the previous period in the vicinity of 480 s. It may become. For this reason, if a cycle in the vicinity of 550 s is erroneously calculated as the reference cycle T, thereafter, the erroneously calculated reference cycle T is used as it is, and therefore the pulse rate continues to be erroneously calculated by the reference cycle T. In the present invention, as described with reference to FIG. 10, even if the reference period T is erroneously calculated, it is appropriately corrected to the optimum reference period T by the synthetic acceleration itotal and the like (see the above formulas 1 to 4). An appropriate pulse rate is calculated based on the optimum reference period T.

なお、周期算出装置1は、脈拍数を算出しているが、同様の処理により心拍数も当然に算出することができる。   Note that the cycle calculating device 1 calculates the pulse rate, but the heart rate can also be naturally calculated by the same processing.

なお、周期算出装置1は、基準周期Tに対応する前回の脈拍数と被験者に作用する合成加速度itotalとを用いて当該基準周期Tを補正しているが、これに限定されずに、次に示す変更を加えてもよい。すなわち、周期算出装置1は、加速度センサ80に替えて、被験者の歩数を計数する歩数計、被験者により吐き出されたガス成分を分析する分析機器にしてもよい。この周期算出装置1は、基準周期Tに対応する前回の脈拍数と、歩数計により計数された被験者の歩数又は分析機器により分析されたガス成分とを用いて当該基準周期Tを補正することができる。   The period calculating device 1 corrects the reference period T by using the previous pulse rate corresponding to the reference period T and the combined acceleration itotal acting on the subject. Changes shown may be made. That is, the period calculation device 1 may be replaced with the acceleration sensor 80 and may be a pedometer that counts the number of steps of the subject or an analysis device that analyzes the gas component exhaled by the subject. This period calculation device 1 can correct the reference period T using the previous pulse rate corresponding to the reference period T and the number of steps of the subject counted by the pedometer or the gas component analyzed by the analysis device. it can.

なお、周期算出装置1は、合成加速度−脈拍数特性である近似式(式2参照)を予め設定しているが、これに限定されずに、次に示す変更を加えてもよい。具体的には、周期算出装置1は、S12において算出した合成加速度itotalと、その合成加速度itotalに対する脈拍数の値とを順次記憶する。   In addition, although the period calculation apparatus 1 presets the approximate expression (refer Formula 2) which is a composite acceleration-pulse rate characteristic, it is not limited to this, You may add the change shown next. Specifically, the cycle calculating device 1 sequentially stores the combined acceleration itotal calculated in S12 and the pulse rate value for the combined acceleration itotal.

そして、周期算出装置1は、その合成加速度itotalと脈拍数の値とが所定期間において所定数分記憶されたことを条件に、記憶された略全ての合成加速度itotalと脈拍数の値とを用いて近似式を算出し、予め設定されている近似式から、算出した近似式に変更する。   Then, the cycle calculating apparatus 1 uses almost all of the stored composite acceleration itotal and the pulse rate value on condition that the combined acceleration itotal and the pulse rate value are stored for a predetermined number of times in a predetermined period. The approximate expression is calculated, and the approximate expression set in advance is changed to the calculated approximate expression.

この場合には、周期算出装置1が、より多くの合成加速度itotalと脈拍数の値とに基づいて近似式を順次更新することにより、より適切な近似式(式2参照)を算出することができる。このため、周期算出装置1は、その近似式を用いて被験者の合成加速度itotalに対する仮基準脈拍数をより正確に算出することができ、結果的には基準脈拍数(式3参照)及び基準周期T(式4参照)もより正確に算出することができる。   In this case, the cycle calculation device 1 may calculate a more appropriate approximate expression (see Expression 2) by sequentially updating the approximate expression based on more synthetic acceleration itotal and the value of the pulse rate. it can. For this reason, the period calculation device 1 can more accurately calculate the provisional reference pulse rate for the synthetic acceleration itotal of the subject using the approximate expression, and as a result, the reference pulse rate (see Expression 3) and the reference period T (see Equation 4) can also be calculated more accurately.

(2)基準周期算出処理2
周期算出装置1は、上記基準周期算出処理1においては被験者に作用する合成加速度itotalを用いて常に基準周期Tを算出(図7参照)しているが、基準周期算出処理2においては被験者に作用する合成加速度itotalを用いて常に基準周期Tを算出(図7参照)しない点で相違する。
(2) Reference cycle calculation process 2
The period calculation device 1 always calculates the reference period T using the combined acceleration itotal acting on the subject in the reference period calculation process 1 (see FIG. 7), but the reference period calculation process 2 acts on the subject. The difference is that the reference cycle T is not always calculated using the resultant acceleration itotal (see FIG. 7).

図12は、本実施形態における基準周期算出処理2を示す図である。図12に示すように、S21において、周期算出装置1は、基準周期Tが所定範囲内に属するか否か判定する。この所定範囲には、予め設定された範囲、基準周期Tの前の周期に対して下限及び上限の間である範囲などが挙げられる。また、周期算出装置1は、この判定がYESである場合にはS22の処理に移り、NOである場合には本処理を終了する。   FIG. 12 is a diagram illustrating the reference cycle calculation process 2 in the present embodiment. As shown in FIG. 12, in S21, the cycle calculation device 1 determines whether or not the reference cycle T belongs to a predetermined range. Examples of the predetermined range include a preset range and a range that is between the lower limit and the upper limit with respect to the cycle before the reference cycle T. Moreover, the period calculation apparatus 1 moves to the process of S22 when this determination is YES, and ends this process when it is NO.

S22乃至S26までの処理は、図7に示すS11乃至S15までの処理と同様であるため、詳細な説明は省略する。   Since the processing from S22 to S26 is the same as the processing from S11 to S15 shown in FIG. 7, detailed description thereof is omitted.

この場合には、周期算出装置1は、周期算出装置1は、基準周期Tが所定範囲内に属しない場合にのみ、加速度センサ80により測定された加速度の値を用いて当該基準周期Tを修正する。このため、周期算出装置1が基準周期Tを逐一補正する処理を実行しないため、周期算出装置1は周期算出装置1における処理負荷を軽減することができる。   In this case, the cycle calculation device 1 corrects the reference cycle T using the acceleration value measured by the acceleration sensor 80 only when the reference cycle T does not belong to the predetermined range. To do. For this reason, since the cycle calculation device 1 does not execute the process of correcting the reference cycle T one by one, the cycle calculation device 1 can reduce the processing load on the cycle calculation device 1.

なお、上記周期算出装置において動作するプログラムは、上記CPU70と同様の機能を有する。なお、プログラムは、記録媒体に記録されてもよい。この記録媒体は、ハードディスク、フレキシブルディスク、コンパクトディスク、ICチップ、カセットテープなどが挙げられる。   Note that the program that operates in the period calculation device has the same function as the CPU 70. The program may be recorded on a recording medium. Examples of the recording medium include a hard disk, a flexible disk, a compact disk, an IC chip, and a cassette tape.

第1実施形態における周期算出装置を示す概略構成図である。It is a schematic block diagram which shows the period calculation apparatus in 1st Embodiment. 第1実施形態における周期算出方法を示す図である。It is a figure which shows the period calculation method in 1st Embodiment. 第1実施形態における周期算出処理を示す図である(その1)。It is a figure which shows the period calculation process in 1st Embodiment (the 1). 第1実施形態における周期算出処理を示す図である(その2)。It is a figure which shows the period calculation process in 1st Embodiment (the 2). 第2実施形態における周期算出装置を示す概略構成図である。It is a schematic block diagram which shows the period calculation apparatus in 2nd Embodiment. 第2実施形態における被験者に作用する加速度の方向を示す図である。It is a figure which shows the direction of the acceleration which acts on the test subject in 2nd Embodiment. 第2実施形態における基準周期算出処理1を示す図である。It is a figure which shows the reference period calculation process 1 in 2nd Embodiment. 第2実施形態における合成加速度−脈拍数特性・算出処理を示す図である。It is a figure which shows the synthetic acceleration-pulse rate characteristic and calculation process in 2nd Embodiment. 第2実施形態における合成加速度に対する脈拍数の関係を示す図である。It is a figure which shows the relationship of the pulse rate with respect to the synthetic | combination acceleration in 2nd Embodiment. 第2実施形態における時間に対する心拍数及び脈拍数の関係を示す図である。It is a figure which shows the relationship of the heart rate and the pulse rate with respect to time in 2nd Embodiment. 第2実施形態における基準周期算出処理2を示す図である。It is a figure which shows the reference period calculation process 2 in 2nd Embodiment. 従来における時間に対する心拍数及び脈拍数の関係を示す図である。It is a figure which shows the relationship between the heart rate and the pulse rate with respect to the time in the past. 従来における脈波の一例を示す図である。It is a figure which shows an example of the conventional pulse wave.

符号の説明Explanation of symbols

1…周期算出装置、10…脈波センサ、20…検出部、30…バンドパスフィルタ部、40…報知部、50…RAM、60…ROM、70…CPU、71…ピーク間隔特定部、72…次ピーク間隔特定部、73…範囲設定部、74…周期算出部、77…補正部、80…加速度センサ、90…検出部   DESCRIPTION OF SYMBOLS 1 ... Period calculation apparatus, 10 ... Pulse wave sensor, 20 ... Detection part, 30 ... Band pass filter part, 40 ... Notification part, 50 ... RAM, 60 ... ROM, 70 ... CPU, 71 ... Peak interval specific | specification part, 72 ... Next peak interval specifying unit, 73 ... range setting unit, 74 ... period calculation unit, 77 ... correction unit, 80 ... acceleration sensor, 90 ... detection unit

Claims (8)

被験者に着脱可能に装着される周期算出装置であって、
第1センサにより測定される基準周期と相関のある被験者の運動量を測定する第2センサと、
前記第2センサにより測定された運動量の値を用いて、前記第1センサにより測定された前記基準周期を補正する補正部と
を備えることを特徴とする周期算出装置。
A cycle calculation device that is detachably attached to a subject,
A second sensor for measuring a subject's momentum correlated with a reference period measured by the first sensor;
A period calculation apparatus comprising: a correction unit that corrects the reference period measured by the first sensor using a value of momentum measured by the second sensor.
前記基準周期は、被験者の脈拍数又は心拍数に対応する周期であり、前記運動量は、互いに直交する3軸方向に作用する加速度の合成値であり、
前記補正部は、多数の被験者の前記加速度の合成値と該多数の被験者の脈拍数又は心拍数との関係により算出された近似式と、前記第2センサにより測定された前記加速度の合成値とを用いて、前記第1センサにより測定された前記基準周期を補正することを特徴とする請求項1に記載の周期算出装置。
The reference period is a period corresponding to the pulse rate or heart rate of the subject, and the momentum is a composite value of accelerations acting in three axis directions orthogonal to each other,
The correction unit includes an approximate expression calculated based on a relationship between a combined value of the accelerations of a large number of subjects and a pulse rate or a heart rate of the large numbers of subjects, and a combined value of the accelerations measured by the second sensor. The period calculation device according to claim 1, wherein the reference period measured by the first sensor is corrected by using the first and second sensors.
前記補正値は、前記第1センサにより測定された前記基準周期が所定範囲内に属しない場合には、前記第2センサにより測定された運動量の値を用いて前記基準周期を補正し、前記第1センサにより測定された前記基準周期が所定範囲内に属する場合には、前記第2センサにより測定された運動量の値を用いて前記基準周期を補正することを特徴とする請求項1に記載の周期算出装置。   When the reference period measured by the first sensor does not fall within a predetermined range, the correction value corrects the reference period using a value of momentum measured by the second sensor, and 2. The reference period according to claim 1, wherein when the reference period measured by one sensor belongs to a predetermined range, the reference period is corrected using a value of momentum measured by the second sensor. Period calculation device. 前記第1センサは、被験者の脈波の振幅の1つのピークから次のピークまでの間隔である前記基準周期を測定するものであり、
前記基準周期に対する下限及び上限の間である特定範囲を設定する範囲設定手段と、
前記第1センサにより測定された前記基準周期に対して次の周期である次周期が前記特定範囲に属する場合には、該次周期を脈波周期として算出する周期算出部と
を備えることを特徴とする請求項1又は請求項3のいずれかに記載の周期算出装置。
The first sensor measures the reference period which is an interval from one peak of the amplitude of the pulse wave of the subject to the next peak,
Range setting means for setting a specific range between a lower limit and an upper limit with respect to the reference period;
A period calculation unit that calculates the next period as a pulse wave period when a next period that is a next period with respect to the reference period measured by the first sensor belongs to the specific range; The period calculation device according to claim 1 or 3.
被験者に着脱可能に装着され、第1センサにより測定される基準周期と相関のある被験者の運動量を測定する第2センサを備える周期算出装置において動作するプログラムであって、
コンピュータを、
前記第2センサにより測定された運動量の値を用いて、前記第1センサにより測定された前記基準周期を補正する補正部として機能させることを特徴とするプログラム。
A program that is detachably attached to a subject and that operates in a period calculation device that includes a second sensor that measures a subject's momentum correlated with a reference period measured by a first sensor,
Computer
A program that functions as a correction unit that corrects the reference period measured by the first sensor using the value of the momentum measured by the second sensor.
前記基準周期は、被験者の脈拍数又は心拍数に対応する周期であり、前記運動量は、互いに直交する3軸方向に作用する加速度の合成値であり、
前記補正部は、多数の被験者の前記加速度の合成値と該多数の被験者の脈拍数又は心拍数との関係により算出された近似式と、前記第2センサにより測定された前記加速度の合成値とを用いて、前記第1センサにより測定された前記基準周期を補正することを特徴とする請求項5に記載のプログラム。
The reference period is a period corresponding to the pulse rate or heart rate of the subject, and the momentum is a composite value of accelerations acting in three axis directions orthogonal to each other,
The correction unit includes an approximate expression calculated based on a relationship between a combined value of the accelerations of a large number of subjects and a pulse rate or a heart rate of the large numbers of subjects, and a combined value of the accelerations measured by the second sensor. The program according to claim 5, wherein the reference period measured by the first sensor is corrected by using.
前記補正値は、前記第1センサにより測定された前記基準周期が所定範囲内に属しない場合には、前記第2センサにより測定された運動量の値を用いて前記基準周期を補正し、前記第1センサにより測定された前記基準周期が所定範囲内に属する場合には、前記第2センサにより測定された運動量の値を用いて前記基準周期を補正することを特徴とする請求項5に記載のプログラム。   When the reference period measured by the first sensor does not fall within a predetermined range, the correction value corrects the reference period using a value of momentum measured by the second sensor, and The said reference period is correct | amended using the value of the momentum measured by the said 2nd sensor, when the said reference period measured by 1 sensor belongs in the predetermined range. program. 前記第1センサは、被験者の脈波の振幅の1つのピークから次のピークまでの間隔である前記基準周期を測定するものであり、
前記基準周期に対する下限及び上限の間である特定範囲を設定する範囲設定手段と、
前記第1センサにより測定された前記基準周期に対して次の周期である次周期が前記特定範囲に属する場合には、該次周期を脈波周期として算出する周期算出部と
を備えることを特徴とする請求項5又は請求項7のいずれかに記載のプログラム。



The first sensor measures the reference period which is an interval from one peak of the amplitude of the pulse wave of the subject to the next peak,
Range setting means for setting a specific range between a lower limit and an upper limit with respect to the reference period;
A period calculation unit that calculates the next period as a pulse wave period when a next period that is a next period with respect to the reference period measured by the first sensor belongs to the specific range; The program according to any one of claims 5 and 7.



JP2005055222A 2005-02-28 2005-02-28 Period calculation device and program Expired - Fee Related JP4326486B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005055222A JP4326486B2 (en) 2005-02-28 2005-02-28 Period calculation device and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005055222A JP4326486B2 (en) 2005-02-28 2005-02-28 Period calculation device and program

Publications (2)

Publication Number Publication Date
JP2006238971A true JP2006238971A (en) 2006-09-14
JP4326486B2 JP4326486B2 (en) 2009-09-09

Family

ID=37045969

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005055222A Expired - Fee Related JP4326486B2 (en) 2005-02-28 2005-02-28 Period calculation device and program

Country Status (1)

Country Link
JP (1) JP4326486B2 (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009213636A (en) * 2008-03-10 2009-09-24 Denso Corp State estimation device
JP2010131209A (en) * 2008-12-04 2010-06-17 Japan Atomic Energy Agency Ear plug type individual heat illness alarming device
JP2010234623A (en) * 2009-03-31 2010-10-21 Brother Ind Ltd Liquid drop ejector
CN106880351A (en) * 2015-12-15 2017-06-23 德州仪器公司 Reduce the artifact of exercise induced in photo-plethysmographic (PPG) signal
WO2018096632A1 (en) * 2016-11-24 2018-05-31 オリンパス株式会社 Data processing device, computer readable medium, data processing method, and program
US10172564B2 (en) 2016-11-24 2019-01-08 Olympus Corporation Apparatus, computer-readable medium, and method for detecting biological data of target patient from attachable sensor attached to target patient
JP2019201888A (en) * 2018-05-23 2019-11-28 倉敷紡績株式会社 Physical condition evaluation system
KR20200058908A (en) * 2018-11-20 2020-05-28 전자부품연구원 Motion Analysis Method for Fitness Posture Correction
JP2022087093A (en) * 2017-02-07 2022-06-09 浩士 中村 Alarm emission system

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009213636A (en) * 2008-03-10 2009-09-24 Denso Corp State estimation device
JP2010131209A (en) * 2008-12-04 2010-06-17 Japan Atomic Energy Agency Ear plug type individual heat illness alarming device
JP2010234623A (en) * 2009-03-31 2010-10-21 Brother Ind Ltd Liquid drop ejector
CN106880351A (en) * 2015-12-15 2017-06-23 德州仪器公司 Reduce the artifact of exercise induced in photo-plethysmographic (PPG) signal
WO2018096632A1 (en) * 2016-11-24 2018-05-31 オリンパス株式会社 Data processing device, computer readable medium, data processing method, and program
US10172564B2 (en) 2016-11-24 2019-01-08 Olympus Corporation Apparatus, computer-readable medium, and method for detecting biological data of target patient from attachable sensor attached to target patient
JP2022087093A (en) * 2017-02-07 2022-06-09 浩士 中村 Alarm emission system
JP7297190B2 (en) 2017-02-07 2023-06-26 浩士 中村 alarm system
JP2019201888A (en) * 2018-05-23 2019-11-28 倉敷紡績株式会社 Physical condition evaluation system
JP7074566B2 (en) 2018-05-23 2022-05-24 倉敷紡績株式会社 Physical condition evaluation system
KR20200058908A (en) * 2018-11-20 2020-05-28 전자부품연구원 Motion Analysis Method for Fitness Posture Correction
KR102123743B1 (en) * 2018-11-20 2020-06-17 전자부품연구원 Motion Analysis Method for Fitness Posture Correction

Also Published As

Publication number Publication date
JP4326486B2 (en) 2009-09-09

Similar Documents

Publication Publication Date Title
JP4326486B2 (en) Period calculation device and program
JP4855721B2 (en) Blood pressure measurement device
US8897864B2 (en) Heart rate meter and method for removing noise of heart beat waveform
US20160345930A1 (en) Blood pressure measurement device and blood pressure measurement method
US20060224074A1 (en) Heartbeat measuring apparatus
CN104114084A (en) Method and system for optical blood pressure monitoring
EP1932473A1 (en) Cough detecting device and cough detecting method
JP6211679B2 (en) Heart rate detector
US20160038117A1 (en) Ultrasonic blood pressure measuring device and blood pressure measuring method
JP6181576B2 (en) Hemodynamic measurement apparatus and hemodynamic measurement method
JP6692027B2 (en) Blood vessel elasticity index value measuring device, blood pressure measuring device, blood vessel elasticity index value measuring method and blood pressure measuring method
CN104027109A (en) Atrial fibrillation analyzer and program
US20150238169A1 (en) Ultrasonic measurement apparatus and ultrasonic measurement method
US10905381B2 (en) Blood pressure correction information generating device, blood pressure measurement device and blood pressure correction information generating method
US20200008690A1 (en) Blood pressure data processing apparatus, blood pressure data processing method, and program
CN110418600A (en) Blood pressure measuring device, method and program
CN104994780A (en) Method for measuring a physiological parameter, such as a biological rhythm, on the basis of at least two sensors
KR20190081650A (en) Apparatus and method for measuring bio information
JP2008079688A (en) Blood pressure measuring instrument
JP6536038B2 (en) Period estimation apparatus, period estimation method and program
CN110811589A (en) Pulse wave evaluation device and pulse wave evaluation method
US7344502B2 (en) Pulse wave measuring apparatus of high measurement accuracy
JP2016055093A (en) Blood pressure measurement apparatus and blood pressure measurement method
WO2020213431A1 (en) Biological information monitoring system, bed system, and biological information monitoring method
JP5092885B2 (en) Electronic blood pressure monitor

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20060725

RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20080118

RD13 Notification of appointment of power of sub attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7433

Effective date: 20080201

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20090421

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: 20090512

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: 20090609

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120619

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120619

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120619

Year of fee payment: 3

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120619

Year of fee payment: 3

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: 20120619

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130619

Year of fee payment: 4

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313113

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313113

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees