JP7160460B2 - Pulse rate calculation method, pulse rate calculation device, recording medium - Google Patents

Pulse rate calculation method, pulse rate calculation device, recording medium Download PDF

Info

Publication number
JP7160460B2
JP7160460B2 JP2021524805A JP2021524805A JP7160460B2 JP 7160460 B2 JP7160460 B2 JP 7160460B2 JP 2021524805 A JP2021524805 A JP 2021524805A JP 2021524805 A JP2021524805 A JP 2021524805A JP 7160460 B2 JP7160460 B2 JP 7160460B2
Authority
JP
Japan
Prior art keywords
pulse rate
pulse
peak
rate calculation
low
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.)
Active
Application number
JP2021524805A
Other languages
Japanese (ja)
Other versions
JPWO2020246370A1 (en
Inventor
重樹 金井
秀一 佐藤
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.)
NEC Solutions Innovators Ltd
Original Assignee
NEC Solutions Innovators 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 NEC Solutions Innovators Ltd filed Critical NEC Solutions Innovators Ltd
Publication of JPWO2020246370A1 publication Critical patent/JPWO2020246370A1/en
Application granted granted Critical
Publication of JP7160460B2 publication Critical patent/JP7160460B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/0245Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals

Description

本発明は、脈拍数算出方法、脈拍数算出装置、記録媒体に関する。 The present invention relates to a pulse rate calculation method, a pulse rate calculation device, and a recording medium.

身体に装着したウェアラブルデバイスを用いて脈拍を測定することにより脈拍数を算出することがある。このようなウェアラブルデバイスを用いた測定を行う場合、体動などによる外乱ノイズが問題となることが知られている。 The pulse rate may be calculated by measuring the pulse using a wearable device worn on the body. It is known that disturbance noise due to body movement or the like poses a problem when performing measurements using such wearable devices.

上記のようなノイズに対する技術として、例えば、特許文献1がある。特許文献1には、脈波センサと体動センサと脈波信号フィルタリング部とフィルタ係数設定部とを有する拍動検出装置が記載されている。特許文献1によると、フィルタ係数設定部は、体動信号に基づいて体動変化が所定の閾値を超えて増大したことを検出すると、適応フィルタの係数を所定値に設定する。また、脈波信号フィルタリング部は、体動信号に基づいて適応フィルタを生成して脈波信号中のノイズ信号を抽出するとともに、脈波信号からノイズを除去した拍動信号を出力する。 For example, Patent Literature 1 is known as a technique for dealing with noise as described above. Patent Literature 1 describes a pulse detection device having a pulse wave sensor, a body motion sensor, a pulse wave signal filtering section, and a filter coefficient setting section. According to Patent Document 1, a filter coefficient setting unit sets a coefficient of an adaptive filter to a predetermined value when it detects that a change in body motion has increased beyond a predetermined threshold based on a body motion signal. The pulse wave signal filtering unit generates an adaptive filter based on the body motion signal to extract a noise signal in the pulse wave signal, and outputs a pulse signal obtained by removing noise from the pulse wave signal.

特開2010-172645号公報JP 2010-172645 A

特許文献1に記載されている技術の場合、ノイズを除去するためには、脈波センサの他に体動センサを有することが必要になる。そのため、体動センサを有さない場合にノイズを除去することができず、容易な構成でノイズの影響を低減させることが難しい、という課題が生じていた。 In the case of the technique described in Patent Document 1, in order to remove noise, it is necessary to have a body motion sensor in addition to the pulse wave sensor. Therefore, the noise cannot be removed without the body motion sensor, and it is difficult to reduce the influence of the noise with a simple configuration.

また、体動等の外乱ノイズを含んだ脈波信号から脈拍数を推定する方法としては、高速フーリエ変換等による周波数領域でノイズ除去する方法などもある。しかしながら、フーリエ変換において時間分解能と周波数分解能はトレードオフである。そのため、フーリエ変換を行う場合、心拍変動解析に必要な時間分解能を保持しつつ、精度よく脈拍数を推定するための周波数分解能を得ることは難しかった。なお、自律神経機能の評価等に用いられる心拍変動解析には、例えば、脈波信号のサンプリングレートが約25Hz程度必要だという報告がある。 Also, as a method of estimating the pulse rate from a pulse wave signal containing disturbance noise such as body motion, there is a method of removing noise in the frequency domain by fast Fourier transform or the like. However, there is a trade-off between time resolution and frequency resolution in the Fourier transform. Therefore, when Fourier transform is performed, it is difficult to obtain frequency resolution for accurately estimating the pulse rate while maintaining the time resolution necessary for heart rate variability analysis. There is a report that, for example, heart rate variability analysis used for evaluation of autonomic nerve function requires a pulse wave signal sampling rate of about 25 Hz.

このように、ウェアラブルデバイスを用いて脈拍数を算出する際に、フーリエ変換に依らず、容易な構成でノイズの影響を低減させることが難しい、という課題が生じていた。そこで、本発明は、ウェアラブルデバイスを用いて脈拍数を算出する際に、フーリエ変換に依らず、容易な構成でノイズの影響を低減させることが難しい、という課題を解決する脈拍数算出方法、脈拍数算出装置、記録媒体を提供することにある。 As described above, when calculating the pulse rate using a wearable device, there is a problem that it is difficult to reduce the influence of noise with a simple configuration without relying on the Fourier transform. Therefore, the present invention provides a pulse rate calculation method that solves the problem that it is difficult to reduce the influence of noise with a simple configuration without relying on the Fourier transform when calculating the pulse rate using a wearable device. An object of the present invention is to provide a number calculation device and a recording medium.

かかる目的を達成するため本発明の一形態である脈拍数算出方法は、
測定した脈波信号に基づいて脈拍数を算出する脈拍数算出装置が行う脈拍数算出方法であって、
減衰させる周波数の異なる複数種類のローパスフィルタを前記脈波信号に適用して、適用した結果に基づいて脈拍数を算出する
という構成をとる。
In order to achieve such an object, a pulse rate calculation method, which is one aspect of the present invention, comprises:
A pulse rate calculation method performed by a pulse rate calculation device that calculates a pulse rate based on a measured pulse wave signal,
A plurality of types of low-pass filters with different attenuation frequencies are applied to the pulse wave signal, and the pulse rate is calculated based on the application result.

また、本発明の他の形態である脈拍数算出装置は、
測定した脈波信号に基づいて脈拍数を算出する脈拍数算出装置であって、
減衰させる周波数の異なる複数種類のローパスフィルタを前記脈波信号に適用して、適用した結果に基づいて脈拍数を算出する脈拍数算出部を有する
という構成をとる。
In addition, a pulse rate calculation device, which is another embodiment of the present invention,
A pulse rate calculator for calculating a pulse rate based on a measured pulse wave signal,
and a pulse rate calculator that applies a plurality of types of low-pass filters with different frequencies to attenuate to the pulse wave signal and calculates the pulse rate based on the result of the application.

また、本発明の他の形態である記録媒体は、
測定した脈波信号に基づいて脈拍数を算出する脈拍数算出装置に、
減衰させる周波数の異なる複数種類のローパスフィルタを前記脈波信号に適用して、適用した結果に基づいて脈拍数を算出する脈拍数算出部を実現するためのプログラムを記録した、コンピュータが読み取り可能な記録媒体である。
A recording medium according to another aspect of the present invention is
A pulse rate calculator that calculates the pulse rate based on the measured pulse wave signal,
A computer readable program recording a program for realizing a pulse rate calculation unit that applies a plurality of types of low-pass filters with different frequencies to attenuate to the pulse wave signal and calculates the pulse rate based on the application results. A recording medium.

本発明は、以上のように構成されることにより、ウェアラブルデバイスを用いて脈拍数を算出する際に、フーリエ変換に依らず、容易な構成でノイズの影響を低減させることが難しい、という課題を解決する脈拍数算出方法、脈拍数算出装置、記録媒体を提供することが可能となる。 By being configured as described above, the present invention solves the problem that it is difficult to reduce the influence of noise with a simple configuration without relying on the Fourier transform when calculating the pulse rate using a wearable device. It is possible to provide a pulse rate calculation method, a pulse rate calculation device, and a recording medium that solve the problem.

本発明の第1の実施形態におけるウェアラブルデバイスの構成の一例を示すブロック図である。1 is a block diagram showing an example of a configuration of a wearable device according to a first embodiment of the invention; FIG. 脈波センサの処理の一例を説明するための図である。It is a figure for demonstrating an example of a process of a pulse wave sensor. 図1で示す処理部の構成の一例を示すブロック図である。2 is a block diagram showing an example of a configuration of a processing unit shown in FIG. 1; FIG. 基線変動を除去する処理の一例を説明するための図である。FIG. 5 is a diagram for explaining an example of processing for removing baseline fluctuations; 基線変動を除去する処理の一例を説明するための図である。FIG. 5 is a diagram for explaining an example of processing for removing baseline fluctuations; ピーク値を検出する処理の一例を説明するための図である。It is a figure for demonstrating an example of the process which detects a peak value. ピーク値を検出する際に用いるピーク間隔の一例を説明するための図である。FIG. 4 is a diagram for explaining an example of peak intervals used when detecting peak values; FIG. フィードバック処理の一例を説明するための図である。FIG. 4 is a diagram for explaining an example of feedback processing; FIG. 本発明の第1の実施形態における処理部の動作の一例を示すフローチャートである。4 is a flow chart showing an example of the operation of a processing unit according to the first embodiment of the present invention; ピーク値を検出する処理を行う際の動作の一例を示すフローチャートである。7 is a flow chart showing an example of an operation when performing processing for detecting a peak value; 脈拍数を算出する処理を行う際の動作の一例を示すフローチャートである。4 is a flow chart showing an example of an operation when performing processing for calculating a pulse rate; フィードバック処理を行う際の動作の一例を示すフローチャートである。7 is a flow chart showing an example of an operation when feedback processing is performed; 本発明の第2の実施形態における脈拍数算出装置の構成の一例を示すブロック図である。FIG. 4 is a block diagram showing an example of the configuration of a pulse rate calculation device according to a second embodiment of the present invention; FIG.

[第1の実施形態]
本発明の第1の実施形態を図1から図12までを参照して説明する。図1は、ウェアラブルデバイス100の構成の一例を示すブロック図である。図2は、脈波センサ210の処理の一例を説明するための図である。図3は、処理部220の構成の一例を示すブロック図である。図4、図5は、基線変動を除去する処理の一例を説明するための図である。図6は、ピーク値を検出する処理の一例を説明するための図である。図7は、ピーク値を検出する際に用いるピーク間隔の一例を説明するための図である。図8は、フィードバック処理の一例を説明するための図である。図9は、処理部220の動作の一例を示すフローチャートである。図10は、ピーク値を検出する処理を行う際の動作の一例を示すフローチャートである。図11は、脈拍数を算出する処理を行う際の動作の一例を示すフローチャートである。図12は、フィードバック処理を行う際の動作の一例を示すフローチャートである。
[First embodiment]
A first embodiment of the present invention will be described with reference to FIGS. 1 to 12. FIG. FIG. 1 is a block diagram showing an example of the configuration of the wearable device 100. As shown in FIG. FIG. 2 is a diagram for explaining an example of processing of pulse wave sensor 210. As shown in FIG. FIG. 3 is a block diagram showing an example of the configuration of the processing section 220. As shown in FIG. 4 and 5 are diagrams for explaining an example of processing for removing baseline fluctuations. FIG. 6 is a diagram for explaining an example of processing for detecting peak values. FIG. 7 is a diagram for explaining an example of peak intervals used when detecting peak values. FIG. 8 is a diagram for explaining an example of feedback processing. FIG. 9 is a flow chart showing an example of the operation of the processing unit 220. As shown in FIG. FIG. 10 is a flow chart showing an example of the operation when performing processing for detecting a peak value. FIG. 11 is a flowchart showing an example of the operation when performing processing for calculating the pulse rate. FIG. 12 is a flow chart showing an example of the operation when feedback processing is performed.

本発明の第1の実施形態においては、脈波センサ210により取得した脈波信号を用いて、1分あたりの脈拍の数である脈拍数を算出するウェアラブルデバイス100について説明する。後述するように、本実施形態におけるウェアラブルデバイス100では、取得した脈波信号にハイパスフィルタ(HPF:High-pass filter)を適用して基線変動を除去した後、減衰させる周波数の異なる2種類のローパスフィルタ(LPS:Low-pass filter)を適用する。そして、ウェアラブルデバイス100は、2種類のローパスフィルタを適用した結果に基づいて脈拍のピークを検出する。その後、ウェアラブルデバイス100は、ピークの検出結果に基づいて求められるピーク間隔から脈拍数を算出する。 In the first embodiment of the present invention, the wearable device 100 that calculates the pulse rate, which is the number of pulses per minute, using the pulse wave signal acquired by the pulse wave sensor 210 will be described. As will be described later, in the wearable device 100 of the present embodiment, after applying a high-pass filter (HPF: High-pass filter) to the acquired pulse wave signal to remove baseline fluctuations, two types of low-pass with different frequencies to be attenuated Apply a filter (LPS: Low-pass filter). Then, wearable device 100 detects a pulse peak based on the result of applying two types of low-pass filters. After that, the wearable device 100 calculates the pulse rate from the peak interval obtained based on the peak detection result.

ウェアラブルデバイス100は、身体に装着して脈波信号を測定する、腕時計型などの情報処理装置である。図1は、ウェアラブルデバイス100の構成の一例を示している。図1を参照すると、ウェアラブルデバイス100には、脈波センサ210と処理部220とを含むセンサボード200が内蔵されている。 The wearable device 100 is an information processing device such as a wristwatch that is worn on the body and measures a pulse wave signal. FIG. 1 shows an example of the configuration of the wearable device 100. As shown in FIG. Referring to FIG. 1 , wearable device 100 incorporates sensor board 200 including pulse wave sensor 210 and processing unit 220 .

脈波センサ210は、例えば、光電式のセンサである。図2は、脈波センサ210による測定の一例を示している。図2で示すように、例えば、脈波センサ210は、緑色のLEDの光を身体に向けて照射する。そして、脈波センサ210は、ウェアラブルデバイス100を装着した身体から反射した光を計測することで、脈波信号を測定する。例えば、以上のように、脈波センサ210は、光電容量脈波方式(PPG:Photoplethysmography)を用いて脈波信号を測定可能なよう構成されている。 Pulse wave sensor 210 is, for example, a photoelectric sensor. FIG. 2 shows an example of measurement by the pulse wave sensor 210. As shown in FIG. As shown in FIG. 2, for example, the pulse wave sensor 210 irradiates the body with green LED light. Pulse wave sensor 210 measures the pulse wave signal by measuring the light reflected from the body on which wearable device 100 is worn. For example, as described above, the pulse wave sensor 210 is configured to be able to measure a pulse wave signal using a photoelectric capacitive pulse wave method (PPG: Photoplethysmography).

なお、自律神経機能評価などに用いられる心拍変動解析には、1秒当たりの測定回数を示す脈波信号のサンプリングレートが約25Hz程度必要であるという報告がある。そのため、本実施形態における脈波信号のサンプリングレートは、25Hz以上あることが望ましい。 It is reported that heart rate variability analysis used for autonomic nerve function evaluation requires a pulse wave signal sampling rate of about 25 Hz, which indicates the number of measurements per second. Therefore, it is desirable that the sampling rate of the pulse wave signal in this embodiment is 25 Hz or more.

処理部220は、脈波センサ210が測定した脈波信号に基づいて脈拍数を算出するマイコンなどの情報処理部(脈拍数算出装置)である。図3は、処理部220の構成の一例を示している。図3を参照すると、処理部220は、例えば、基線変動除去部221とピーク検出部222と脈拍数算出部223とフィードバック部224とを有している。 The processing unit 220 is an information processing unit (pulse rate calculation device) such as a microcomputer that calculates the pulse rate based on the pulse wave signal measured by the pulse wave sensor 210 . FIG. 3 shows an example of the configuration of the processing unit 220. As shown in FIG. Referring to FIG. 3, the processing section 220 has, for example, a baseline fluctuation removal section 221, a peak detection section 222, a pulse rate calculation section 223, and a feedback section 224. FIG.

例えば、上述した各処理部は、ピーク検出を行うピーク検出回路や脈拍数を算出してフィードバックを行う脈拍数推定回路などの各種論理回路、ハードウェアにより実現することが出来る。上述した各処理部は、例えば、記憶装置に格納されたプログラムをCPU(Central Processing Unit)などの演算装置が実行することにより実現されても構わない。 For example, each of the processing units described above can be realized by various logic circuits and hardware, such as a peak detection circuit that performs peak detection and a pulse rate estimation circuit that calculates and feeds back the pulse rate. Each of the processing units described above may be implemented by executing a program stored in a storage device by an arithmetic device such as a CPU (Central Processing Unit).

基線変動除去部221は、脈波センサ210が測定した脈波信号にハイパスフィルタを適用する。これにより、基線変動除去部221は、脈波センサ210が測定した脈波信号から基線変動を除去する。 The baseline fluctuation remover 221 applies a high-pass filter to the pulse wave signal measured by the pulse wave sensor 210 . Thereby, the baseline fluctuation removing unit 221 removes the baseline fluctuation from the pulse wave signal measured by the pulse wave sensor 210 .

例えば、図4で示すように、脈波センサ210が測定した脈波信号では、ウェアラブルデバイス100を装着した装着者が動いていることなどにより、基線が揺らいでいることがある。そこで、基線変動除去部221は、例えば、予め定められた0.5Hzなどのハイパスフィルタを適用する。これにより、図5で示すように、脈波信号の基線を0に揃えたフラットな状態にすることが出来る。 For example, as shown in FIG. 4, in the pulse wave signal measured by the pulse wave sensor 210, the baseline may fluctuate due to the wearer wearing the wearable device 100 moving. Therefore, the baseline fluctuation removing unit 221 applies a predetermined high-pass filter such as 0.5 Hz, for example. As a result, as shown in FIG. 5, a flat state in which the baseline of the pulse wave signal is aligned to 0 can be obtained.

なお、基線変動除去部221が減衰させる周波数の値は、信号の種類などに応じて適宜変更可能なよう構成して構わない。つまり、基線変動除去部221が用いるハイパスフィルタは、0.5Hzのハイパスフィルタに限定されない。 It should be noted that the value of the frequency attenuated by the baseline fluctuation removing section 221 may be configured to be appropriately changeable according to the type of signal. That is, the high-pass filter used by the baseline fluctuation removing section 221 is not limited to the 0.5 Hz high-pass filter.

ピーク検出部222は、減衰させる周波数の異なる2種類のローパスフィルタを適用した後、適用した結果に基づいて脈拍のピークを検出する。 The peak detection unit 222 applies two types of low-pass filters with different attenuation frequencies, and then detects a pulse peak based on the application results.

例えば、ピーク検出部222は、基線変動除去部221により基線変動を除去した脈波信号を2つにコピーする。そして、ピーク検出部222は、コピーした2つのうちの一方に3Hzのローパスフィルタを適用するとともに、コピーした2つのうちの他方に10Hzのローパスフィルタを適用する。 For example, the peak detector 222 makes two copies of the pulse wave signal from which the baseline fluctuation has been removed by the baseline fluctuation remover 221 . Then, the peak detection unit 222 applies a 3 Hz low-pass filter to one of the two copies, and applies a 10 Hz low-pass filter to the other of the two copies.

ピーク検出部222は、3Hzのローパスフィルタを適用した脈波信号から後述する推定範囲内に存在するピークを検出する。また、3Hzのローパスフィルタを適用した脈波信号において推定範囲内にピークが存在しない場合、ピーク検出部222は、10Hzのローパスフィルタを適用した脈波信号から、同様の推定範囲内に存在するピークを検出する。このように、ピーク検出部222は、3Hzのローパスフィルタを適用した脈波信号と10Hzのローパスフィルタを適用した脈波信号とを用いて、ピークを検出する。なお、ピーク検出部222によるピークの検出は、脈波信号の接線の傾きの変化を検出するなど既知の方法を用いて行って構わない。 A peak detector 222 detects a peak present within an estimated range described later from the pulse wave signal to which the 3 Hz low-pass filter is applied. Also, if the pulse wave signal to which the 3 Hz low-pass filter is applied does not have a peak within the estimation range, the peak detection unit 222 detects the peak present within the same estimation range from the pulse wave signal to which the 10 Hz low-pass filter is applied. to detect In this way, the peak detector 222 detects peaks using the pulse wave signal to which the 3 Hz low-pass filter is applied and the pulse wave signal to which the 10 Hz low-pass filter is applied. Note that the peak detection by the peak detection unit 222 may be performed using a known method such as detecting a change in the slope of the tangent line of the pulse wave signal.

ここで、図6を参照して、推定範囲について詳細に説明する。図6は、推定範囲について説明するための図である。図6を参照すると、ピーク検出部222は、基準となる値である脈拍数(推定用)に基づいて、次のピークを推定する範囲となる推定範囲を設定する。例えば、ピーク検出部222は、脈拍数(推定用)に所定の値を加算した場合の、前回のピークから次のピークまでのピーク間隔を算出する。また、ピーク検出部222は、脈拍数(推定用)から所定の値を減算した場合のピーク間隔を算出する。そして、ピーク検出部222は、脈拍数(推定用)に所定の値を加算した場合のピーク間隔と脈拍数(推定用)から所定の値を減算した場合のピーク間隔との間を推定範囲として設定する。 Now, with reference to FIG. 6, the estimation range will be described in detail. FIG. 6 is a diagram for explaining the estimated range. Referring to FIG. 6, the peak detection unit 222 sets an estimation range, which is a range for estimating the next peak, based on the pulse rate (for estimation), which is a reference value. For example, the peak detector 222 calculates the peak interval from the previous peak to the next peak when a predetermined value is added to the pulse rate (for estimation). The peak detector 222 also calculates the peak interval when a predetermined value is subtracted from the pulse rate (for estimation). Then, the peak detection unit 222 defines an estimation range between the peak interval when a predetermined value is added to the pulse rate (for estimation) and the peak interval when a predetermined value is subtracted from the pulse rate (for estimation). set.

また、ピーク検出部222は、脈拍数(推定用)から基準値を算出することが出来る。基準値は、脈拍数(推定用)である場合のピーク間隔を示す値である。 Also, the peak detector 222 can calculate the reference value from the pulse rate (for estimation). A reference value is a value indicating a peak interval in the case of a pulse rate (for estimation).

例えば、ピーク間隔をサンプリングしたデータの数で表すものとすると、ピーク間隔は、(サンプリングレート[Hz])/(脈拍数(推定用)に必要に応じて加算、または、減算した値)×60[秒]で求めることが出来る。例えば、加算する値が30bpm、減算する値が25bpmとして予め定められており、脈拍数(推定用)として初期値である80bpmを用いるとする。また、サンプリングレートが110Hzであるとする。この場合、ピーク検出部222は、110/(80+30)×60=60を算出するとともに、110/(80-25)×60=120を算出する。そして、ピーク検出部222は、前回のピークからサンプルデータの数が60~120の間を推定範囲として設定する。また、ピーク検出部222は、基準値として、110/80×60=82.5を算出することが出来る。 For example, if the peak interval is represented by the number of sampled data, the peak interval is (sampling rate [Hz]) / (pulse rate (for estimation) added or subtracted as necessary) x 60 It can be obtained in [seconds]. For example, it is assumed that the value to be added is 30 bpm and the value to be subtracted is 25 bpm, and the initial value of 80 bpm is used as the pulse rate (for estimation). Also assume that the sampling rate is 110 Hz. In this case, the peak detector 222 calculates 110/(80+30)×60=60 and 110/(80−25)×60=120. Then, the peak detection unit 222 sets the estimated range to be between 60 and 120 pieces of sample data from the previous peak. Also, the peak detector 222 can calculate 110/80×60=82.5 as the reference value.

なお、ノイズなどの影響により、推定範囲内に複数のピークが検出される場合がある。この場合、ピーク検出部222は、検出された複数のピークのうち、上述した基準値により近いピークを採用するよう構成することが出来る。 A plurality of peaks may be detected within the estimation range due to the influence of noise or the like. In this case, the peak detector 222 can be configured to adopt a peak closer to the above-described reference value among the plurality of detected peaks.

また、10Hzのローパスフィルタを適用した脈波信号において、推定範囲内にピークが存在しない場合もある。この場合、ピーク検出部222は、ピーク検出不可と判断する。そして、ピーク検出部222は、例えば、基準値の位置にピークが存在するものとして扱うことが出来る。 Also, there are cases where the pulse wave signal to which the 10 Hz low-pass filter is applied does not have a peak within the estimation range. In this case, the peak detector 222 determines that the peak cannot be detected. Then, the peak detection unit 222 can handle, for example, that a peak exists at the position of the reference value.

以上をまとめると、例えば、図7で示すようになる。図7は、ピーク検出部222によるピーク検出処理の一例を示している。図7を参照すると、ピーク検出部222は、3Hzのローパスフィルタを適用した脈波信号を用いて、前回のピークから推定範囲内に存在するピークを検出する。3Hzのローパスフィルタを適用した脈波信号において、推定範囲内にピークが存在する場合、ピーク検出部222は、10Hzのローパスフィルタを適用した脈波信号の確認を行わない。また、3Hzのローパスフィルタを適用した脈波信号において推定範囲内にピークを検出できない場合、ピーク検出部222は、10Hzのローパスフィルタを適用した脈波信号を用いて、同様の推定範囲内のピークを検出する。また、10Hzのローパスフィルタを適用した脈波信号においても推定範囲内にピークが存在しない場合、ピーク検出部222は、基準値の位置にピークが存在するものとして扱う。例えば、以上のような処理を繰り返すことにより、ピーク検出部222は、図7で示すように、3Hzのローパスフィルタを適用した脈波信号と10Hzのローパスフィルタを適用した脈波信号とに基づくピークの検出を行う。 If the above is put together, it will come to show in FIG. 7, for example. FIG. 7 shows an example of peak detection processing by the peak detector 222. As shown in FIG. Referring to FIG. 7, the peak detector 222 detects a peak existing within an estimated range from the previous peak using the pulse wave signal to which a 3 Hz low-pass filter is applied. If the pulse wave signal to which the 3 Hz low-pass filter is applied has a peak within the estimation range, the peak detector 222 does not check the pulse wave signal to which the 10 Hz low-pass filter is applied. Also, if the pulse wave signal to which the 3 Hz low-pass filter is applied cannot detect a peak within the estimation range, the peak detection unit 222 uses the pulse wave signal to which the 10 Hz low-pass filter is applied to detect the peak within the same estimation range. to detect Also, if the pulse wave signal to which the 10 Hz low-pass filter is applied does not have a peak within the estimation range, the peak detector 222 treats it as having a peak at the position of the reference value. For example, by repeating the above processing, the peak detection unit 222 detects a peak based on the pulse wave signal to which a 3 Hz low-pass filter is applied and the pulse wave signal to which a 10 Hz low-pass filter is applied, as shown in FIG. detection.

なお、脈拍数(推定用)の初期値、加算する値、減算する値は、それぞれ例示した以外の値であっても構わない。また、脈拍数(推定用)は後述するフィードバック処理により適宜変更されることになる。つまり、基準となる値である脈拍数(推定用)は、脈拍数の算出結果に基づいて適宜変更される。 Note that the initial value of the pulse rate (for estimation), the value to be added, and the value to be subtracted may be values other than those exemplified. In addition, the pulse rate (for estimation) is appropriately changed by feedback processing, which will be described later. In other words, the pulse rate (for estimation), which is a reference value, is appropriately changed based on the calculation result of the pulse rate.

また、ピーク検出部222が減衰させる周波数の値は、信号の種類などに応じて適宜変更可能なよう構成して構わない。つまり、ピーク検出部222が用いるローパスフィルタは、3Hzや10Hzのローパスフィルタに限定されない。 Further, the value of the frequency attenuated by the peak detection section 222 may be configured to be appropriately changeable according to the type of signal or the like. That is, the low-pass filter used by the peak detector 222 is not limited to the low-pass filters of 3 Hz and 10 Hz.

脈拍数算出部223は、ピーク検出部222が検出したピークに基づいて脈拍数を算出する。例えば、脈拍数算出部223は、前回のピークとピーク検出部222が検出したピークとの間のピーク間隔から、脈拍数を算出する。 Pulse rate calculator 223 calculates the pulse rate based on the peak detected by peak detector 222 . For example, the pulse rate calculator 223 calculates the pulse rate from the peak interval between the previous peak and the peak detected by the peak detector 222 .

例えば、脈拍数算出部223は、前回のピークとピーク検出部222が検出したピークとの間のサンプリングデータの個数を計測することで、ピーク間隔を算出する。そして、脈拍数算出部223は、算出したピーク間隔に基づいて脈拍数を算出する。具体的には、例えば、脈拍数算出部223は、(サンプリングレート[Hz])/(ピーク間隔)×60[秒]を計算することで、脈拍数を算出する。例えば、サンプリングレートが110Hz、ピーク間隔が110であるとすると、脈拍数算出部223は、110/110×60=60を脈拍数として算出する。 For example, the pulse rate calculator 223 calculates the peak interval by counting the number of pieces of sampling data between the previous peak and the peak detected by the peak detector 222 . Then, the pulse rate calculator 223 calculates the pulse rate based on the calculated peak interval. Specifically, for example, the pulse rate calculator 223 calculates the pulse rate by calculating (sampling rate [Hz])/(peak interval)×60 [seconds]. For example, if the sampling rate is 110 Hz and the peak interval is 110, the pulse rate calculator 223 calculates 110/110×60=60 as the pulse rate.

なお、3Hzのローパスフィルタを適用した脈波信号と10Hzのローパスフィルタを適用した脈波信号のいずれにおいても推定範囲内にピークが存在しない場合、脈拍数算出部223は、新たな脈拍数の算出を行わず、前回算出した脈拍数を採用するよう構成しても構わない。 If neither the pulse wave signal to which the 3 Hz low-pass filter is applied nor the pulse wave signal to which the 10 Hz low-pass filter is applied does not have a peak within the estimation range, the pulse rate calculation unit 223 calculates a new pulse rate. is not performed, and the previously calculated pulse rate may be used.

フィードバック部224は、脈拍数算出部223が算出した脈拍数に基づくフィードバック処理を行う。例えば、フィードバック部224は、脈拍数算出部223が算出した脈拍数に基づいて、推定範囲を設定する際に用いた脈拍数(推定用)の変更を行う。 The feedback section 224 performs feedback processing based on the pulse rate calculated by the pulse rate calculation section 223 . For example, the feedback unit 224 changes the pulse rate (for estimation) used when setting the estimation range based on the pulse rate calculated by the pulse rate calculation unit 223 .

例えば、フィードバック部224は、脈拍数算出部223が算出した脈拍数が前回算出した脈拍数に所定値を加算した値より大きいか、又は、脈拍数算出部223が算出した脈拍数が前回算出した脈拍数から所定値を減算した値より小さいか、確認する。そして、脈拍数算出部223が算出した脈拍数が前回算出した脈拍数に所定値を加算した値より大きい場合、フィードバック部224は、脈拍数(推定用)に予め定められた値を加算して新たな脈拍数(推定用)とする。また、脈拍数算出部223が算出した脈拍数が前回算出した脈拍数から所定値を減算した値より小さい場合、フィードバック部224は、脈拍数(推定用)から予め定められた値を減算して新たな脈拍数(推定用)とする。 For example, the feedback unit 224 determines whether the pulse rate calculated by the pulse rate calculation unit 223 is greater than the value obtained by adding a predetermined value to the pulse rate calculated last time, or the pulse rate calculated by the pulse rate calculation unit 223 is greater than the pulse rate calculated last time. Check if it is less than the value obtained by subtracting a predetermined value from the pulse rate. If the pulse rate calculated by the pulse rate calculator 223 is greater than the value obtained by adding a predetermined value to the previously calculated pulse rate, the feedback section 224 adds a predetermined value to the pulse rate (for estimation). New pulse rate (for estimation). If the pulse rate calculated by the pulse rate calculator 223 is smaller than the value obtained by subtracting a predetermined value from the previously calculated pulse rate, the feedback section 224 subtracts a predetermined value from the pulse rate (for estimation). New pulse rate (for estimation).

図8は、上述したフィードバック部224によるフィードバック処理の一例を示している。図8を参照すると、例えば、脈拍数算出部223が算出した脈拍数が前回算出した脈拍数に所定値である5bpmを加算した値より大きい場合、フィードバック部224は、脈拍数(推定用)に予め定められた値である3bpmを加算して新たな脈拍数(推定用)とする。また、脈拍数算出部223が算出した脈拍数が前回算出した脈拍数から所定値である5bpmを減算した値より小さい場合、フィードバック部224は、脈拍数(推定用)から予め定められた値である3bpmを減算して新たな脈拍数(推定用)とする。なお、脈拍数算出部223が算出した脈拍数が前回算出した脈拍数から所定値である5bpmを減算した値以上前回算出した脈拍数に所定値である5bpmを加算した値以下である場合、フィードバック部224は、前回の脈拍数(推定用)をそのまま新たな脈拍数(推定用)とすることが出来る。 FIG. 8 shows an example of feedback processing by the feedback section 224 described above. Referring to FIG. 8, for example, when the pulse rate calculated by the pulse rate calculation unit 223 is greater than the value obtained by adding a predetermined value of 5 bpm to the previously calculated pulse rate, the feedback unit 224 changes the pulse rate (for estimation) to A new pulse rate (for estimation) is obtained by adding a predetermined value of 3 bpm. Further, when the pulse rate calculated by the pulse rate calculation unit 223 is smaller than the value obtained by subtracting the predetermined value of 5 bpm from the pulse rate calculated last time, the feedback unit 224 calculates the pulse rate (for estimation) with a predetermined value. A certain 3 bpm is subtracted to obtain a new pulse rate (for estimation). If the pulse rate calculated by the pulse rate calculation unit 223 is greater than or equal to the value obtained by subtracting a predetermined value of 5 bpm from the previously calculated pulse rate and less than or equal to the value obtained by adding the predetermined value of 5 bpm to the previously calculated pulse rate, feedback is performed. The unit 224 can use the previous pulse rate (for estimation) as a new pulse rate (for estimation).

例えば、以上のように、フィードバック部224は、脈拍数算出部223が算出した脈拍数に基づいて、推定範囲を設定する際に用いる脈拍数(推定用)に対するフィードバック処理を行う。なお、フィードバック部224は、脈拍数算出部223が算出した脈拍数を脈拍数(推定用)とするなど、上記例示した以外の方法でフィードバック処理を行うよう構成しても構わない。 For example, as described above, the feedback unit 224 performs feedback processing on the pulse rate (for estimation) used when setting the estimation range, based on the pulse rate calculated by the pulse rate calculation unit 223 . Note that the feedback unit 224 may be configured to perform feedback processing by a method other than the above example, such as using the pulse rate calculated by the pulse rate calculation unit 223 as the pulse rate (for estimation).

なお、上述した各種所定値、予め定められた値は、それぞれ例示である。フィードバック部224がフィードバック処理を行う際に用いる値は、例示した以外のものであっても構わない。 Note that the various predetermined values and predetermined values described above are examples. Values used when the feedback unit 224 performs feedback processing may be values other than those illustrated.

以上が、処理部220の構成の一例である。続いて、図9から図12までを参照して、処理部220の動作の一例について説明する。まず、図9を参照して、処理部220が脈拍数を算出する際の全体の動作の流れについて説明する。 The above is an example of the configuration of the processing unit 220 . Next, an example of the operation of the processing unit 220 will be described with reference to FIGS. 9 to 12. FIG. First, referring to FIG. 9, the flow of the overall operation when the processing unit 220 calculates the pulse rate will be described.

図9を参照すると、基線変動除去部221は、脈波センサ210が測定した脈波信号にハイパスフィルタを適用する。これにより、基線変動除去部221は、脈波センサ210が測定した脈波信号から基線変動を除去する(ステップS101)。 Referring to FIG. 9 , baseline fluctuation remover 221 applies a high-pass filter to the pulse wave signal measured by pulse wave sensor 210 . Thereby, the baseline fluctuation removing unit 221 removes the baseline fluctuation from the pulse wave signal measured by the pulse wave sensor 210 (step S101).

ピーク検出部222は、減衰させる周波数の異なる2種類のローパスフィルタを適用する(ステップS202)。例えば、ピーク検出部222は、基線変動除去部221により基線変動を除去した脈波信号を2つにコピーする。そして、ピーク検出部222は、コピーした2つのうちの一方に3Hzのローパスフィルタを適用するとともに、コピーした2つのうちの他方に10Hzのローパスフィルタを適用する。 The peak detector 222 applies two types of low-pass filters with different attenuation frequencies (step S202). For example, the peak detector 222 makes two copies of the pulse wave signal from which the baseline fluctuation has been removed by the baseline fluctuation remover 221 . Then, the peak detection unit 222 applies a 3 Hz low-pass filter to one of the two copies, and applies a 10 Hz low-pass filter to the other of the two copies.

ピーク検出部222は、3Hzのローパスフィルタを適用した脈波信号と10Hzのローパスフィルタを適用した脈波信号とを用いて、ピークを検出する(ステップS103)。例えば、ピーク検出部222は、まず、3Hzのローパスフィルタを適用した脈波信号を用いて、前回のピークから推定範囲内に存在するピークを検出しようとする。また、3Hzのローパスフィルタを適用した脈波信号において、前回のピークから推定範囲内に次のピークが存在しない場合、ピーク検出部222は、10Hzのローパスフィルタを適用した脈波信号を用いて、前回のピークから推定範囲内に存在するピークを検出しようとする。 The peak detector 222 detects peaks using the pulse wave signal to which the 3 Hz low-pass filter is applied and the pulse wave signal to which the 10 Hz low-pass filter is applied (step S103). For example, the peak detection unit 222 first tries to detect a peak existing within the estimated range from the previous peak using the pulse wave signal to which a 3 Hz low-pass filter is applied. Also, in the pulse wave signal to which the 3 Hz low-pass filter is applied, if the next peak does not exist within the estimation range from the previous peak, the peak detection unit 222 uses the pulse wave signal to which the 10 Hz low-pass filter is applied, Attempt to find peaks that are within the estimated range from the previous peak.

脈拍数算出部223は、前回のピークとピーク検出部222が検出したピークとの間のピーク間隔から、脈拍数を算出する(ステップS104)。例えば、脈拍数算出部223は、(サンプリングレート[Hz])/(ピーク間隔)×60[秒]を計算することで、脈拍数を算出する。 The pulse rate calculator 223 calculates the pulse rate from the peak interval between the previous peak and the peak detected by the peak detector 222 (step S104). For example, the pulse rate calculator 223 calculates the pulse rate by calculating (sampling rate [Hz])/(peak interval)×60 [seconds].

フィードバック部224は、脈拍数算出部223が算出した脈拍数に基づくフィードバック処理を行う(ステップS105)。例えば、フィードバック部224は、脈拍数算出部223が算出した脈拍数に基づいて、推定範囲を設定する際に用いた脈拍数(推定用)の変更を行う。 The feedback unit 224 performs feedback processing based on the pulse rate calculated by the pulse rate calculation unit 223 (step S105). For example, the feedback unit 224 changes the pulse rate (for estimation) used when setting the estimation range based on the pulse rate calculated by the pulse rate calculation unit 223 .

処理部220は、推定範囲を用いたピークの検出、ピーク間隔に応じた脈拍数の算出、脈拍数(推定用)を変更するフィードバック処理、の各処理を繰り返すことで、リアルタイムで脈拍数の算出処理を行う。続いて、図10を参照して、ステップS103の処理についてより詳細に説明する。 The processing unit 220 repeats each process of detecting a peak using the estimation range, calculating the pulse rate according to the peak interval, and feedback processing for changing the pulse rate (for estimation), thereby calculating the pulse rate in real time. process. Next, with reference to FIG. 10, the process of step S103 will be described in more detail.

図10を参照すると、ピーク検出部222は、3Hzのローパスフィルタを適用した脈波信号から推定範囲内に存在するピークを検出する(ステップS201)。 Referring to FIG. 10, the peak detector 222 detects a peak existing within the estimation range from the pulse wave signal to which the 3 Hz low-pass filter is applied (step S201).

3Hzのローパスフィルタを適用した脈波信号において、推定範囲内にピークが存在する場合(ステップS202、Yes)、ピーク検出部222は、ピークの検出を終了する。一方、推定範囲内にピークが存在しない場合(ステップS202、No)、ピーク検出部222は、10Hzのローパスフィルタを適用した脈波信号から、同様の推定範囲内に存在するピークを検出する(ステップS203)。 If the pulse wave signal to which the 3 Hz low-pass filter is applied has a peak within the estimation range (step S202, Yes), the peak detection unit 222 terminates peak detection. On the other hand, if there is no peak within the estimated range (step S202, No), the peak detection unit 222 detects a peak within the same estimated range from the pulse wave signal to which the 10 Hz low-pass filter is applied (step S203).

10Hzのローパスフィルタを適用した脈波信号において推定範囲内にピークが存在する場合(ステップS204、Yes)、ピーク検出部222は、ピークの検出を終了する。一方、推定範囲内にピークが存在しない場合(ステップS204、No)、ピーク検出部222は、ピーク検出不可と判断する(ステップS205)。この場合、ピーク検出部222は、例えば、基準値の位置にピークが存在するものとして扱うよう構成することが出来る。 If the pulse wave signal to which the 10 Hz low-pass filter is applied has a peak within the estimation range (step S204, Yes), the peak detection unit 222 terminates peak detection. On the other hand, if there is no peak within the estimated range (step S204, No), the peak detector 222 determines that the peak cannot be detected (step S205). In this case, the peak detection unit 222 can be configured to handle, for example, a peak at the position of the reference value.

続いて、図11を参照して、ステップS104の処理についてより詳細に説明する。 Next, with reference to FIG. 11, the process of step S104 will be described in more detail.

図11を参照すると、ピーク検出部222によりピークの検出が出来ていた場合(ステップS301、Yes)、脈拍数算出部223は、前回のピークとピーク検出部222が検出したピークとの間のピーク間隔を算出する(ステップS302)。例えば、脈拍数算出部223は、ピーク間のサンプリングデータの個数を計測することで、ピーク間隔を算出する。 Referring to FIG. 11, when the peak detection unit 222 has successfully detected a peak (step S301, Yes), the pulse rate calculation unit 223 detects the peak between the previous peak and the peak detected by the peak detection unit 222. Intervals are calculated (step S302). For example, the pulse rate calculator 223 calculates the peak interval by counting the number of sampling data between peaks.

脈拍数算出部223は、算出したピーク間隔を用いて、脈拍数を算出する。例えば、脈拍数算出部223は、(サンプリングレート[Hz])/(ピーク間隔)×60[秒]を計算することで、脈拍数を算出する(ステップS303)。 The pulse rate calculator 223 calculates the pulse rate using the calculated peak interval. For example, the pulse rate calculator 223 calculates the pulse rate by calculating (sampling rate [Hz])/(peak interval)×60 [seconds] (step S303).

一方、ピーク検出部222によるピークの検出が出来ていなかった場合(ステップS301、No)、脈拍数算出部223は、前回算出した脈拍数を採用する(ステップS304)。 On the other hand, if the peak detection unit 222 could not detect a peak (step S301, No), the pulse rate calculation unit 223 adopts the pulse rate calculated last time (step S304).

続いて、図12を参照して、ステップS105の処理についてより詳細に説明する。 Next, with reference to FIG. 12, the process of step S105 will be described in more detail.

図12を参照すると、脈拍数算出部223が算出した脈拍数が前回算出した脈拍数に所定値である5bpmを加算した値より大きい場合(ステップS401、Yes)、フィードバック部224は、脈拍数(推定用)に予め定められた値である3bpmを加算する(ステップS402)。そして、フィードバック部204は、脈拍数(推定用)に予め定められた値である3bpmを加算した値を、新たな脈拍数(推定用)として設定する(ステップS406)。 Referring to FIG. 12, when the pulse rate calculated by the pulse rate calculator 223 is greater than the value obtained by adding a predetermined value of 5 bpm to the previously calculated pulse rate (step S401, Yes), the feedback section 224 outputs the pulse rate ( 3 bpm, which is a predetermined value, is added to (for estimation) (step S402). Then, the feedback unit 204 sets a value obtained by adding a predetermined value of 3 bpm to the pulse rate (for estimation) as a new pulse rate (for estimation) (step S406).

脈拍数算出部223が算出した脈拍数が前回算出した脈拍数に所定値である5bpmを加算した値以下である場合(ステップS401、No)、フィードバック部224は、脈拍数算出部223が算出した脈拍数が前回算出した脈拍数から所定値である5bpmを減算した値より小さいか否か確認する。 If the pulse rate calculated by the pulse rate calculation unit 223 is equal to or less than the value obtained by adding the predetermined value of 5 bpm to the previously calculated pulse rate (step S401, No), the feedback unit 224 outputs the pulse rate calculated by the pulse rate calculation unit 223. It is checked whether or not the pulse rate is smaller than a value obtained by subtracting a predetermined value of 5 bpm from the previously calculated pulse rate.

脈拍数算出部223が算出した脈拍数が前回算出した脈拍数から所定値である5bpmを減算した値より小さい場合(ステップS403、Yes)、フィードバック部224は、脈拍数(推定用)から予め定められた値である3bpmを減算する(ステップS404)。そして、フィードバック部204は、脈拍数(推定用)から予め定められた値である3bpmを減算した値を、新たな脈拍数(推定用)として設定する(ステップS406)。一方、脈拍数算出部223が算出した脈拍数が前回算出した脈拍数から所定値である5bpmを減算した値以上である場合(ステップS403、No)、フィードバック部224は、前回の脈拍数(推定用)をそのまま新たな脈拍数(推定用)として(ステップS405)設定する(ステップS406)。 If the pulse rate calculated by the pulse rate calculation unit 223 is smaller than the value obtained by subtracting the predetermined value of 5 bpm from the pulse rate calculated last time (step S403, Yes), the feedback unit 224 determines a predetermined value from the pulse rate (for estimation). 3 bpm, which is the obtained value, is subtracted (step S404). The feedback unit 204 then sets a value obtained by subtracting a predetermined value of 3 bpm from the pulse rate (for estimation) as a new pulse rate (for estimation) (step S406). On the other hand, when the pulse rate calculated by the pulse rate calculation unit 223 is equal to or greater than the value obtained by subtracting the predetermined value of 5 bpm from the pulse rate calculated last time (step S403, No), the feedback unit 224 returns the previous pulse rate (estimated (for estimation) is set as it is (step S405) as a new pulse rate (for estimation) (step S406).

このように、処理部220は、ピーク検出部222と脈拍数算出部223とを有している。このような構成により、ピーク検出部222は、2種類のローパスフィルタを適用した結果に基づいて脈拍のピークを検出することが出来る。また、脈拍数算出部223は、ピーク検出部222が2種類のローパスフィルタを用いて検出したピークに基づいて、脈拍数を算出することが出来る。その結果、より確実にピークを検出することが可能となり、より精度よく脈拍数を算出することが可能となる。つまり、上記構成によると、体動センサなどを有さなくても、また、フーリエ変換を用いることなく、容易な構成によりノイズの影響を低減させて精度よく脈拍数を算出することが可能となる。 Thus, the processing section 220 has a peak detection section 222 and a pulse rate calculation section 223 . With such a configuration, the peak detector 222 can detect a pulse peak based on the results of applying two types of low-pass filters. The pulse rate calculator 223 can also calculate the pulse rate based on the peaks detected by the peak detector 222 using two types of low-pass filters. As a result, the peak can be detected more reliably, and the pulse rate can be calculated more accurately. That is, according to the above configuration, it is possible to reduce the influence of noise and accurately calculate the pulse rate with a simple configuration without using a body motion sensor or the like and without using Fourier transform. .

また、ピーク検出部222は、推定範囲に基づいてピークを検出するよう構成されている。このような構成により、より精度よくピークを検出することが可能となる。また、推定範囲は、フィードバック部224によるフィードバック処理により更新される脈拍数(推定用)に基づいて設定されるよう構成されている。このような構成により、より精度よくピークを検出することが可能となる。 Also, the peak detector 222 is configured to detect a peak based on the estimated range. With such a configuration, peaks can be detected with higher accuracy. Also, the estimation range is configured to be set based on the pulse rate (for estimation) updated by feedback processing by the feedback section 224 . With such a configuration, peaks can be detected with higher accuracy.

なお、本実施形態においては、ピーク検出部222が2種類のローパスフィルタを適用する場合について例示した。しかしながら、ピーク検出部222が適用するローパスフィルタの数は、2種類に限定されない。例えば、ピーク検出部222は、3種類以上の複数のローパスフィルタを適用して、ピークを検出するよう構成しても構わない。 In addition, in this embodiment, the peak detection part 222 illustrated the case where two types of low-pass filters are applied. However, the number of low-pass filters applied by peak detector 222 is not limited to two. For example, the peak detection unit 222 may be configured to detect peaks by applying a plurality of low-pass filters of three or more types.

また、本実施形態においては、脈拍数算出部223が脈拍数を算出するごとにフィードバック部224によりフィードバック処理が行われる場合について例示した。しかしながら、フィードバック部224によるフィードバック処理は、必ずしも毎回行うよう構成しなくても構わない。 Moreover, in the present embodiment, the case where the feedback processing is performed by the feedback section 224 each time the pulse rate calculation section 223 calculates the pulse rate has been exemplified. However, the feedback processing by the feedback unit 224 does not necessarily have to be performed every time.

[第2の実施形態]
次に、図13を参照して、本発明の第2の実施形態について説明する。第2の実施形態では、脈拍数算出装置30の構成の概要について説明する。
[Second embodiment]
A second embodiment of the present invention will now be described with reference to FIG. In the second embodiment, an overview of the configuration of the pulse rate calculator 30 will be described.

脈拍数算出装置30は、測定した脈波信号に基づいて脈拍数を算出する情報処理装置(または、回路装置)である。図13を参照すると、脈拍数算出装置30は、例えば、脈拍数算出部31を有している。 The pulse rate calculator 30 is an information processing device (or circuit device) that calculates the pulse rate based on the measured pulse wave signal. Referring to FIG. 13, the pulse rate calculator 30 has a pulse rate calculator 31, for example.

上述した処理部は、例えば、論理回路などのハードウェアにより実現することが出来る。上述した処理部は、例えば、記憶装置に格納されたプログラムをCPUなどの演算装置が実行することにより実現されても構わない。 The processing unit described above can be realized by hardware such as a logic circuit, for example. The processing unit described above may be implemented by, for example, executing a program stored in a storage device by an arithmetic device such as a CPU.

脈拍数算出部31は、脈波信号を取得すると、減衰させる周波数の異なる複数種類のローパスフィルタを取得した脈波信号に適用する。そして、脈拍数算出部31は、適用した結果に基づいて、脈拍数を算出する。 After acquiring the pulse wave signal, the pulse rate calculator 31 applies a plurality of types of low-pass filters with different frequencies to attenuate to the acquired pulse wave signal. Then, the pulse rate calculator 31 calculates the pulse rate based on the result of the application.

このように、脈拍数算出装置30は、脈拍数算出部31を有している。このような構成により、脈拍数算出部31は、減衰させる周波数の異なる複数種類のローパスフィルタを脈波信号に適用した結果に基づいて、脈拍数を算出することが出来る。その結果、より確実にピークを検出することが可能となり、より精度よく脈拍数を算出することが可能となる。これにより、フーリエ変換を用いることなく、容易な構成によりノイズの影響を低減させて精度よく脈拍数を算出することが可能となる。 Thus, the pulse rate calculator 30 has the pulse rate calculator 31 . With such a configuration, the pulse rate calculator 31 can calculate the pulse rate based on the result of applying a plurality of types of low-pass filters with different attenuation frequencies to the pulse wave signal. As a result, the peak can be detected more reliably, and the pulse rate can be calculated more accurately. This makes it possible to reduce the influence of noise and accurately calculate the pulse rate with a simple configuration without using Fourier transform.

また、上述した脈拍数算出装置30は、当該脈拍数算出装置30に所定のプログラムが組み込まれることで実現できる。具体的に、本発明の他の形態であるプログラムは、取得した脈波信号に基づいて脈拍数を算出する脈拍数算出装置30に、減衰させる周波数の異なる複数種類のローパスフィルタを脈波信号に適用して、適用した結果に基づいて脈拍数を算出する脈拍数算出部31を実現するためのプログラムである。 Further, the pulse rate calculation device 30 described above can be realized by incorporating a predetermined program into the pulse rate calculation device 30 . Specifically, the program, which is another embodiment of the present invention, applies a plurality of types of low-pass filters having different frequencies to attenuate to the pulse wave signal in the pulse rate calculation device 30 that calculates the pulse rate based on the acquired pulse wave signal. It is a program for realizing a pulse rate calculation unit 31 that applies and calculates the pulse rate based on the result of application.

また、上述した脈拍数算出装置30により実行される脈拍数算出方法は、取得した脈波信号に基づいて脈拍数を算出する脈拍数算出装置が行う脈拍数算出方法であって、減衰させる周波数の異なる複数種類のローパスフィルタを脈波信号に適用して、適用した結果に基づいて脈拍数を算出する、という方法である。 Further, the pulse rate calculation method executed by the pulse rate calculation device 30 described above is a pulse rate calculation method performed by the pulse rate calculation device that calculates the pulse rate based on the acquired pulse wave signal. In this method, different types of low-pass filters are applied to the pulse wave signal, and the pulse rate is calculated based on the results of the application.

上述した構成を有する、プログラム、又は、脈拍数算出方法、の発明であっても、上記脈拍数算出装置30と同様の作用・効果を有するために、上述した本発明の目的を達成することが出来る。 Even the invention of the program or the pulse rate calculation method having the configuration described above can achieve the above-described object of the present invention because it has the same action and effect as the pulse rate calculation device 30 described above. I can.

<付記>
上記実施形態の一部又は全部は、以下の付記のようにも記載されうる。以下、本発明における脈拍数算出方法などの概略を説明する。但し、本発明は、以下の構成に限定されない。
<Appendix>
Some or all of the above embodiments may also be described as the following appendices. An outline of the pulse rate calculation method and the like according to the present invention will be described below. However, the present invention is not limited to the following configurations.

(付記1)
取得した脈波信号に基づいて脈拍数を算出する脈拍数算出装置が行う脈拍数算出方法であって、
減衰させる周波数の異なる複数種類のローパスフィルタを前記脈波信号に適用して、適用した結果に基づいて脈拍数を算出する
脈拍数算出方法。
(付記2)
付記1に記載の脈拍数算出方法であって、
複数種類のローパスフィルタを適用した結果に基づいて検出した脈拍のピークに基づいて脈拍数を算出する
脈拍数算出方法。
(付記3)
付記2に記載の脈拍数算出方法であって、
基準となる値に基づいて推定範囲を設定し、前回のピークから前記推定範囲内にある脈拍のピークを検出する
脈拍数算出方法。
(付記4)
付記3に記載の脈拍数算出方法であって、
前記基準となる値に所定の値を加算した場合のピーク間隔と、前記基準となる値から所定の値を減算した場合のピーク間隔と、を算出し、算出した結果に基づいて前記推定範囲を設定する
脈拍数算出方法。
(付記5)
付記3又は付記4に記載の脈拍数算出方法であって、
前記基準となる値は、脈拍数の算出結果に基づいて変更される値である
脈拍数算出方法。
(付記6)
付記3から付記5までのいずれか1項に記載の脈拍数算出方法であって、
あるローパスフィルタを適用した結果において前記推定範囲内に脈拍のピークが存在しない場合、他のローパスフィルタを適用した結果において前記推定範囲内に脈拍のピークが存在するか確認する
脈拍数算出方法。
(付記7)
付記2から付記6までのいずれか1項に記載の脈拍数算出方法であって、
検出した脈拍のピークと、前回の脈拍のピークと、の間のピーク間隔に基づいて脈拍数を算出する
脈拍数算出方法。
(付記8)
付記1から付記7までのいずれか1項に記載の脈拍数算出方法であって、
前記脈波信号に所定のハイパスフィルタを適用した後、減衰させる周波数の異なる複数種類のローパスフィルタを適用する
脈拍数算出方法。
(付記9)
取得した脈波信号に基づいて脈拍数を算出する脈拍数算出装置であって、
減衰させる周波数の異なる複数種類のローパスフィルタを前記脈波信号に適用して、適用した結果に基づいて脈拍数を算出する脈拍数算出部を有する
脈拍数算出装置。
(付記10)
取得した脈波信号に基づいて脈拍数を算出する脈拍数算出装置に、
減衰させる周波数の異なる複数種類のローパスフィルタを前記脈波信号に適用して、適用した結果に基づいて脈拍数を算出する脈拍数算出部を実現するためのプログラム。
(Appendix 1)
A pulse rate calculation method performed by a pulse rate calculation device that calculates a pulse rate based on an acquired pulse wave signal,
A pulse rate calculation method, comprising: applying a plurality of types of low-pass filters with different frequencies to attenuate to the pulse wave signal, and calculating the pulse rate based on the results of the application.
(Appendix 2)
The pulse rate calculation method according to Supplementary Note 1,
A pulse rate calculation method for calculating a pulse rate based on a pulse peak detected based on a result of applying a plurality of types of low-pass filters.
(Appendix 3)
The pulse rate calculation method according to appendix 2,
A pulse rate calculation method, comprising: setting an estimation range based on a reference value; and detecting a pulse peak within the estimation range from a previous peak.
(Appendix 4)
The pulse rate calculation method according to appendix 3,
Calculating a peak interval when a predetermined value is added to the reference value and a peak interval when a predetermined value is subtracted from the reference value, and determining the estimated range based on the calculated result Sets the pulse rate calculation method.
(Appendix 5)
The pulse rate calculation method according to Supplementary Note 3 or Supplementary Note 4,
The pulse rate calculation method, wherein the reference value is a value that is changed based on a pulse rate calculation result.
(Appendix 6)
The pulse rate calculation method according to any one of appendices 3 to 5,
A pulse rate calculation method for determining whether a peak pulse exists within the estimated range in a result obtained by applying another low-pass filter when a pulse peak does not exist in the estimated range in a result obtained by applying a certain low-pass filter.
(Appendix 7)
The pulse rate calculation method according to any one of appendices 2 to 6,
A pulse rate calculation method for calculating a pulse rate based on a peak interval between a detected pulse peak and a previous pulse peak.
(Appendix 8)
The pulse rate calculation method according to any one of appendices 1 to 7,
A pulse rate calculation method comprising: applying a predetermined high-pass filter to the pulse wave signal, and then applying a plurality of types of low-pass filters having different frequencies to attenuate.
(Appendix 9)
A pulse rate calculation device that calculates a pulse rate based on the acquired pulse wave signal,
A pulse rate calculation device, comprising: a pulse rate calculator that applies a plurality of types of low-pass filters with different frequencies to attenuate to the pulse wave signal and calculates a pulse rate based on the result of the application.
(Appendix 10)
A pulse rate calculator that calculates the pulse rate based on the acquired pulse wave signal,
A program for implementing a pulse rate calculator that applies a plurality of types of low-pass filters with different frequencies to attenuate to the pulse wave signal and calculates the pulse rate based on the results of the application.

なお、上記各実施形態及び付記において記載したプログラムは、記憶装置に記憶されていたり、コンピュータが読み取り可能な記録媒体に記録されていたりする。例えば、記録媒体は、フレキシブルディスク、光ディスク、光磁気ディスク、及び、半導体メモリ等の可搬性を有する媒体である。 It should be noted that the programs described in each of the above embodiments and supplementary notes are stored in a storage device or recorded in a computer-readable recording medium. For example, the recording medium is a portable medium such as a flexible disk, an optical disk, a magneto-optical disk, and a semiconductor memory.

以上、上記各実施形態を参照して本願発明を説明したが、本願発明は、上述した実施形態に限定されるものではない。本願発明の構成や詳細には、本願発明の範囲内で当業者が理解しうる様々な変更をすることが出来る。 Although the present invention has been described with reference to the above-described embodiments, the present invention is not limited to the above-described embodiments. Various changes that can be understood by those skilled in the art can be made to the configuration and details of the present invention within the scope of the present invention.

なお、本発明は、日本国にて2019年6月5日に特許出願された特願2019-104970の特許出願に基づく優先権主張の利益を享受するものであり、当該特許出願に記載された内容は、全て本明細書に含まれるものとする。 In addition, the present invention enjoys the benefit of claiming priority based on the patent application of Japanese Patent Application No. 2019-104970 filed on June 5, 2019 in Japan, and is described in the patent application. The contents are hereby incorporated by reference in their entirety.

100 ウェアラブルデバイス
200 センサボード
210 脈波センサ
220 処理部
221 基線変動除去部
222 ピーク検出部
223 脈拍数算出部
224 フィードバック部
30 脈拍数算出装置
31 脈拍数算出部
100 wearable device 200 sensor board 210 pulse wave sensor 220 processing unit 221 baseline fluctuation removal unit 222 peak detection unit 223 pulse rate calculation unit 224 feedback unit 30 pulse rate calculation device 31 pulse rate calculation unit

Claims (8)

取得した脈波信号に基づいて脈拍数を算出する脈拍数算出装置が行う脈拍数算出方法であって、
減衰させる周波数の異なる複数種類のローパスフィルタを前記脈波信号に適用して、複数種類のローパスフィルタを適用した結果に基づいて検出した脈拍のピークに基づいて脈拍数を算出し、
脈拍のピークを検出する際は、基準となる値に基づいて推定範囲を設定し、前回のピークから前記推定範囲内にある脈拍のピークを検出し、
前記基準となる値は、脈拍数の算出結果に基づいて変更される値である
脈拍数算出方法。
A pulse rate calculation method performed by a pulse rate calculation device that calculates a pulse rate based on an acquired pulse wave signal,
Applying a plurality of types of low-pass filters with different frequencies to be attenuated to the pulse wave signal, and calculating a pulse rate based on the pulse peak detected based on the results of applying the plurality of types of low-pass filters,
When detecting a pulse peak, an estimation range is set based on a reference value, and a pulse peak within the estimation range from the previous peak is detected,
The reference value is a value that is changed based on the pulse rate calculation result.
Pulse rate calculation method.
請求項に記載の脈拍数算出方法であって、
前記基準となる値に所定の値を加算した場合のピーク間隔と、前記基準となる値から所定の値を減算した場合のピーク間隔と、を算出し、算出した結果に基づいて前記推定範囲を設定する
脈拍数算出方法。
The pulse rate calculation method according to claim 1 ,
Calculating a peak interval when a predetermined value is added to the reference value and a peak interval when a predetermined value is subtracted from the reference value, and determining the estimated range based on the calculated result Sets the pulse rate calculation method.
請求項1または請求項2に記載の脈拍数算出方法であって、
あるローパスフィルタを適用した結果において前記推定範囲内に脈拍のピークが存在しない場合、他のローパスフィルタを適用した結果において前記推定範囲内に脈拍のピークが存在するか確認する
脈拍数算出方法。
The pulse rate calculation method according to claim 1 or claim 2 ,
A pulse rate calculation method for determining whether a peak pulse exists within the estimated range in a result obtained by applying another low-pass filter when a pulse peak does not exist in the estimated range in a result obtained by applying a certain low-pass filter.
請求項1から請求項3までのいずれか1項に記載の脈拍数算出方法であって、
検出した脈拍のピークと、前回の脈拍のピークと、の間のピーク間隔に基づいて脈拍数を算出する
脈拍数算出方法。
The pulse rate calculation method according to any one of claims 1 to 3 ,
A pulse rate calculation method for calculating a pulse rate based on a peak interval between a detected pulse peak and a previous pulse peak.
請求項1から請求項までのいずれか1項に記載の脈拍数算出方法であって、
前記脈波信号に所定のハイパスフィルタを適用した後、減衰させる周波数の異なる複数種類のローパスフィルタを適用する
脈拍数算出方法。
The pulse rate calculation method according to any one of claims 1 to 4 ,
A pulse rate calculation method comprising: applying a predetermined high-pass filter to the pulse wave signal, and then applying a plurality of types of low-pass filters having different frequencies to attenuate.
取得した脈波信号に基づいて脈拍数を算出する脈拍数算出装置であって、
減衰させる周波数の異なる複数種類のローパスフィルタを前記脈波信号に適用して、複数種類のローパスフィルタを適用した結果に基づいて検出した脈拍のピークに基づいて脈拍数を算出する脈拍数算出部を有し、
前記脈拍数算出部は、脈拍のピークを検出する際は、基準となる値に基づいて推定範囲を設定し、前回のピークから前記推定範囲内にある脈拍のピークを検出し、前記基準となる値は、脈拍数の算出結果に基づいて変更される値である
脈拍数算出装置。
A pulse rate calculation device that calculates a pulse rate based on the acquired pulse wave signal,
a pulse rate calculation unit that applies a plurality of types of low-pass filters with different frequencies to attenuate to the pulse wave signal, and calculates a pulse rate based on a pulse peak detected based on the results of applying the plurality of types of low-pass filters ; have
When detecting a pulse peak, the pulse rate calculator sets an estimation range based on a reference value, detects a pulse peak within the estimation range from the previous peak, and uses it as the reference. value is the value that changes based on the pulse rate calculation result
Pulse rate calculator.
取得した脈波信号に基づいて脈拍数を算出する脈拍数算出装置に、
減衰させる周波数の異なる複数種類のローパスフィルタを前記脈波信号に適用して、複数種類のローパスフィルタを適用した結果に基づいて検出した脈拍のピークに基づいて脈拍数を算出し、
脈拍のピークを検出する際は、基準となる値に基づいて推定範囲を設定し、前回のピークから前記推定範囲内にある脈拍のピークを検出し、
前記基準となる値は、脈拍数の算出結果に基づいて変更される値である
処理を実現するためのプログラム
A pulse rate calculator that calculates the pulse rate based on the acquired pulse wave signal,
Applying a plurality of types of low-pass filters with different frequencies to be attenuated to the pulse wave signal, and calculating a pulse rate based on the pulse peak detected based on the results of applying the plurality of types of low-pass filters,
When detecting a pulse peak, an estimation range is set based on a reference value, and a pulse peak within the estimation range from the previous peak is detected,
The reference value is a value that is changed based on the pulse rate calculation result.
A program for realizing processing .
取得した脈波信号に基づいて脈拍数を算出する脈拍数算出装置が行う脈拍数算出方法であって、
減衰させる周波数の異なる複数種類のローパスフィルタを前記脈波信号に適用して、複数種類のローパスフィルタを適用した結果に基づいて検出した脈拍のピークに基づいて脈拍数を算出し、
脈拍のピークを検出する際は、基準となる値に基づいて推定範囲を設定し、前回のピークから前記推定範囲内にある脈拍のピークを検出し、
複数種類のローパスフィルタを前記脈波信号に適用した結果に基づいて脈拍のピークを検出する際は、あるローパスフィルタを適用した結果において前記推定範囲内に脈拍のピークが存在しない場合に、他のローパスフィルタを適用した結果において前記推定範囲内に脈拍のピークが存在するか確認する
脈拍数算出方法。
A pulse rate calculation method performed by a pulse rate calculation device that calculates a pulse rate based on an acquired pulse wave signal,
Applying a plurality of types of low-pass filters with different frequencies to be attenuated to the pulse wave signal, and calculating a pulse rate based on the pulse peak detected based on the results of applying the plurality of types of low-pass filters,
When detecting a pulse peak, an estimation range is set based on a reference value, and a pulse peak within the estimation range from the previous peak is detected,
When detecting a pulse peak based on the results of applying a plurality of types of low-pass filters to the pulse wave signal, if there is no pulse peak within the estimated range in the result of applying a certain low-pass filter, other A pulse rate calculation method for confirming whether a pulse peak exists within the estimated range in the result of applying a low-pass filter .
JP2021524805A 2019-06-05 2020-05-28 Pulse rate calculation method, pulse rate calculation device, recording medium Active JP7160460B2 (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2019104970 2019-06-05
JP2019104970 2019-06-05
PCT/JP2020/021210 WO2020246370A1 (en) 2019-06-05 2020-05-28 Pulse rate calculating method, pulse rate calculating device, and recording medium

Publications (2)

Publication Number Publication Date
JPWO2020246370A1 JPWO2020246370A1 (en) 2021-12-02
JP7160460B2 true JP7160460B2 (en) 2022-10-25

Family

ID=73652204

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021524805A Active JP7160460B2 (en) 2019-06-05 2020-05-28 Pulse rate calculation method, pulse rate calculation device, recording medium

Country Status (2)

Country Link
JP (1) JP7160460B2 (en)
WO (1) WO2020246370A1 (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006271731A (en) 2005-03-29 2006-10-12 Toshiba Corp Heartbeat measuring apparatus and heartbeat measuring method
JP2009189562A (en) 2008-02-14 2009-08-27 Citizen Holdings Co Ltd Pulse wave measuring apparatus
US20160324477A1 (en) 2015-05-08 2016-11-10 Texas Instruments Incorporated Accuracy of heart rate estimation from photoplethysmographic (ppg) signals
JP2017213123A (en) 2016-05-31 2017-12-07 日本電信電話株式会社 Blood flow rate measurement device and blood flow rate measurement method
WO2018066342A1 (en) 2016-10-07 2018-04-12 株式会社村田製作所 Handheld pulse wave measurement device

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006271731A (en) 2005-03-29 2006-10-12 Toshiba Corp Heartbeat measuring apparatus and heartbeat measuring method
JP2009189562A (en) 2008-02-14 2009-08-27 Citizen Holdings Co Ltd Pulse wave measuring apparatus
US20160324477A1 (en) 2015-05-08 2016-11-10 Texas Instruments Incorporated Accuracy of heart rate estimation from photoplethysmographic (ppg) signals
JP2017213123A (en) 2016-05-31 2017-12-07 日本電信電話株式会社 Blood flow rate measurement device and blood flow rate measurement method
WO2018066342A1 (en) 2016-10-07 2018-04-12 株式会社村田製作所 Handheld pulse wave measurement device

Also Published As

Publication number Publication date
JPWO2020246370A1 (en) 2021-12-02
WO2020246370A1 (en) 2020-12-10

Similar Documents

Publication Publication Date Title
KR102532764B1 (en) Apparatus and method for estimating biophysiological rates
US7455643B1 (en) Continuous non-invasive blood pressure measurement apparatus and methods providing automatic recalibration
Castells-Rufas et al. Simple real-time QRS detector with the MaMeMi filter
US8655436B2 (en) Heart rate meter and heart beat detecting method
WO2016057781A1 (en) Method and apparatus for non-contact fast vital sign acquisition based on radar signal
WO2007032226A1 (en) Heart rate meter and method for removing noise of heart beat waveform
US20130006150A1 (en) Bruxism detection device and bruxism detection method
US10135424B2 (en) Digital filter with confidence input
TWI609671B (en) Signal detection method
CN107743048B (en) Signal processing system for removing OTDR noise
Li et al. A UV-visible absorption spectrum denoising method based on EEMD and an improved universal threshold filter
JP7160460B2 (en) Pulse rate calculation method, pulse rate calculation device, recording medium
CN111920429B (en) Mental stress detection method and device and electronic equipment
CN107405086B (en) Measurement device, measurement method, and program
TW201801671A (en) Method of heart rate detection
TWI593238B (en) Impulsive noise detection circuit and method thereof
JP7428605B2 (en) Respiratory heart rate measuring device and respiratory heart rate measuring program
JP2013183845A (en) Pulsation detector, electronic device and program
TWI488127B (en) Method and device for detection of touch points in touch panel
JP2013208311A (en) Pulsation detector, electronic equipment, and program
JP2013208312A (en) Signal analyzer, electronic equipment and program
JP5854551B2 (en) Real-time frequency analysis method
CN110403593B (en) Heart rate detection method
JP6152806B2 (en) Biological information measurement method
JP2013202076A (en) Pulsation detection device, electronic equipment and program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210727

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220628

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220824

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20221006

R151 Written notification of patent or utility model registration

Ref document number: 7160460

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151