JP5760876B2 - Atrial fibrillation determination device, atrial fibrillation determination method and program - Google Patents

Atrial fibrillation determination device, atrial fibrillation determination method and program Download PDF

Info

Publication number
JP5760876B2
JP5760876B2 JP2011194904A JP2011194904A JP5760876B2 JP 5760876 B2 JP5760876 B2 JP 5760876B2 JP 2011194904 A JP2011194904 A JP 2011194904A JP 2011194904 A JP2011194904 A JP 2011194904A JP 5760876 B2 JP5760876 B2 JP 5760876B2
Authority
JP
Japan
Prior art keywords
atrial fibrillation
waveform signal
power
determination
unit
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
JP2011194904A
Other languages
Japanese (ja)
Other versions
JP2013055982A (en
JP2013055982A5 (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.)
Seiko Epson Corp
Original Assignee
Seiko Epson Corp
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 Seiko Epson Corp filed Critical Seiko Epson Corp
Priority to JP2011194904A priority Critical patent/JP5760876B2/en
Priority to US13/604,064 priority patent/US20130060154A1/en
Publication of JP2013055982A publication Critical patent/JP2013055982A/en
Publication of JP2013055982A5 publication Critical patent/JP2013055982A5/ja
Application granted granted Critical
Publication of JP5760876B2 publication Critical patent/JP5760876B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は、心房細動を判定する技術に関する。   The present invention relates to a technique for determining atrial fibrillation.

心臓疾患に関する医療分野において、心房細動を判定する技術がある。特許文献1には、1拍毎の心電図から得られるRR間隔を計測し、その標準偏差と度数分布とに基づいて、心房細動を判定する技術が開示されている。非特許文献1には、心房細動はRR間隔が不規則であり、心房細動の心拍の周波数解析をすると1/fβ成分が存在し、そのゆらぎのため白色雑音状になることが記載されている。   There is a technique for determining atrial fibrillation in the medical field related to heart disease. Patent Document 1 discloses a technique for measuring an RR interval obtained from an electrocardiogram for each beat and determining atrial fibrillation based on the standard deviation and frequency distribution. Non-Patent Document 1 describes that atrial fibrillation has an irregular RR interval, and a frequency analysis of the heartbeat of atrial fibrillation has a 1 / fβ component, resulting in white noise due to the fluctuation. ing.

特開2009−89883号公報JP 2009-89883 A

Hayano J, Yamasaki F, Sakata S, Okada A, Mukai S, Fujinami T 「Spectral characteristics of ventricular response to atrial fibrillation.」 Am J Physiol 1997; 273 : H2811-H2816Hayano J, Yamasaki F, Sakata S, Okada A, Mukai S, Fujinami T `` Spectral characteristics of ventricular response to atrial fibrillation. '' Am J Physiol 1997; 273: H2811-H2816

上述した特許文献1および非特許文献1においては、正確な心房細動の判定には、1拍毎のRR間隔を正確に計測する必要がある。このRR間隔の計測にあたっては、心電を測定して得られる心電図の波形信号から計測することも可能であるが、脈波を測定して得られる脈波信号から計測することも可能である。   In Patent Document 1 and Non-Patent Document 1 described above, it is necessary to accurately measure the RR interval for each beat for accurate determination of atrial fibrillation. In measuring the RR interval, it is possible to measure from the waveform signal of the electrocardiogram obtained by measuring the electrocardiogram, but it is also possible to measure from the pulse wave signal obtained by measuring the pulse wave.

ところが、脈波を測定する場合には、測定中に被験者が自由に動きまわることができる場合が多いことから、体動ノイズの影響が脈波信号に含まれやすい。心電を測定する場合においても、脈波を測定する場合と比べて程度の差はあるものの、体動ノイズの影響が心電図の波形信号に含まれる場合がある。このように体動ノイズの影響を受けた場合には、1拍毎のRR間隔を正確に計測することは非常に困難である。
そのため、特許文献1および非特許文献1に開示された技術のように、1拍毎の正確なRR間隔を計測することを前提としている場合には、体動ノイズの影響が含まれる信号を用いて心房細動の判定を行うことはできなかった。
However, when measuring a pulse wave, the subject can often move freely during the measurement, so the influence of body motion noise is likely to be included in the pulse wave signal. Even when measuring electrocardiogram, there is a case where the influence of body motion noise is included in the waveform signal of the electrocardiogram, although there is a difference in degree compared to the case of measuring pulse wave. Thus, when affected by body motion noise, it is very difficult to accurately measure the RR interval for each beat.
Therefore, when it is assumed that an accurate RR interval for each beat is measured as in the techniques disclosed in Patent Document 1 and Non-Patent Document 1, a signal including the influence of body motion noise is used. Therefore, atrial fibrillation could not be determined.

本発明は、上述した事情に鑑みてなされたものであり、その目的の一つは、脈波信号および心電図の波形信号などのRR間隔が計測可能な信号において体動ノイズの影響が含まれていても、その信号から心房細動を判定することにある。   The present invention has been made in view of the above-described circumstances, and one of the purposes thereof includes the influence of body motion noise in a signal capable of measuring an RR interval such as a pulse wave signal and an electrocardiogram waveform signal. However, it is to determine atrial fibrillation from the signal.

上述の課題を解決するため、本発明は、心電または脈波の検出結果を示す検出波形信号を取得する取得部と、前記取得された検出波形信号の周波数解析により得られる各フレームのスペクトルに基づいて、フレーム毎にRR間隔に相当するパラメーターを算出し、当該パラメーターの時間変化を示すRR波形信号を算出するRR間隔算出部と、前記RR波形信号における予め決められた周波数帯域のパワーの時間変化を算出するパワー算出部と、前記算出されたパワーが特定の条件を満たすか否かを判定し、判定結果に応じて心房細動の有無に関する情報を出力する判定部とを具備することを特徴とする心房細動判定装置を提供する。
この心房細動判定装置によれば、脈波信号および心電図の波形信号などのRR間隔が計測可能な信号において体動ノイズの影響が含まれていても、その信号から心房細動を判定することができる。
In order to solve the above-described problem, the present invention provides an acquisition unit that acquires a detection waveform signal indicating a detection result of an electrocardiogram or a pulse wave, and a spectrum of each frame obtained by frequency analysis of the acquired detection waveform signal. A parameter corresponding to the RR interval for each frame, an RR interval calculation unit for calculating an RR waveform signal indicating a time change of the parameter, and a time of power in a predetermined frequency band in the RR waveform signal A power calculation unit that calculates a change, and a determination unit that determines whether or not the calculated power satisfies a specific condition and outputs information on the presence or absence of atrial fibrillation according to the determination result. A characteristic atrial fibrillation determination device is provided.
According to this atrial fibrillation determination device, even if the influence of body motion noise is included in a signal capable of measuring an RR interval such as a pulse wave signal and an electrocardiogram waveform signal, atrial fibrillation is determined from the signal. Can do.

また、別の好ましい態様において、前記周波数帯域の最低周波数は、前記フレームの時間の逆数以上であることを特徴とする。
この心房細動判定装置によれば、心房細動の判定の精度を向上させることができる。
In another preferred aspect, the lowest frequency of the frequency band is equal to or greater than the reciprocal of the time of the frame.
According to this atrial fibrillation determination device, the accuracy of determination of atrial fibrillation can be improved.

また、別の好ましい態様において、検出対象者の前記心電または前記脈波を検出する検出部と、前記判定部によって出力された情報に基づいてユーザーに報知する報知部とを具備し、前記取得部は、前記検出された結果に応じて得られる検出波形信号を取得することを特徴とする。
この心房細動判定装置によれば、検出対象者がリアルタイムに心房細動の判定結果を確認することができる。
Moreover, in another preferable aspect, the detection unit includes a detection unit that detects the electrocardiogram or the pulse wave of a detection target person, and a notification unit that notifies a user based on information output by the determination unit, and the acquisition The unit acquires a detected waveform signal obtained according to the detected result.
According to this atrial fibrillation determination device, the person to be detected can confirm the determination result of atrial fibrillation in real time.

また、別の好ましい態様において、前記取得部は、前記検出波形信号に対して、体動ノイズ成分を低減するフィルター処理を施し、前記検出波形信号として出力するノイズ低減部を備えることを特徴とする。
この心房細動判定装置によれば、心房細動の判定の精度を向上させつつ、検出対象者がリアルタイムに心房細動の判定結果を確認することができる。
In another preferred aspect, the acquisition unit includes a noise reduction unit that performs a filtering process to reduce a body movement noise component on the detected waveform signal and outputs the detected waveform signal as the detected waveform signal. .
According to this atrial fibrillation determination device, the detection subject can confirm the determination result of atrial fibrillation in real time while improving the accuracy of determination of atrial fibrillation.

また、本発明は、心電または脈波の検出結果を示す検出波形信号を取得する取得ステップと、前記取得された検出波形信号の周波数解析により得られる各フレームのスペクトルに基づいて、フレーム毎にRR間隔に相当するパラメーターを算出し、当該パラメーターの時間変化を示すRR波形信号を算出するRR間隔算出ステップと、前記RR波形信号における予め決められた周波数帯域のパワーの時間変化を算出するパワー算出ステップと、前記算出されたパワーが特定の条件を満たすか否かを判定し、判定結果に応じて心房細動の有無に関する情報を出力する判定ステップとを具備することを特徴とする心房細動判定方法を提供する。
この心房細動判定方法によれば、脈波信号および心電図の波形信号などのRR間隔が計測可能な信号において体動ノイズの影響が含まれていても、その信号から心房細動を判定することができる。
Further, the present invention provides an acquisition step for acquiring a detection waveform signal indicating a detection result of an electrocardiogram or a pulse wave, and a frame for each frame based on a spectrum of each frame obtained by frequency analysis of the acquired detection waveform signal. An RR interval calculating step for calculating a parameter corresponding to the RR interval, calculating an RR waveform signal indicating a time change of the parameter, and a power calculation for calculating a time change of power in a predetermined frequency band in the RR waveform signal. Atrial fibrillation comprising: a step; and a determination step of determining whether or not the calculated power satisfies a specific condition and outputting information on the presence or absence of atrial fibrillation according to the determination result A determination method is provided.
According to this atrial fibrillation determination method, even if the influence of body motion noise is included in a signal capable of measuring an RR interval such as a pulse wave signal and an electrocardiogram waveform signal, atrial fibrillation is determined from the signal. Can do.

また、本発明は、コンピューターを、心電または脈波の検出結果を示す検出波形信号を取得する取得部と、前記取得された検出波形信号の周波数解析により得られる各フレームのスペクトルに基づいて、フレーム毎にRR間隔に相当するパラメーターを算出し、当該パラメーターの時間変化を示すRR波形信号を算出するRR間隔算出部と、前記RR波形信号における予め決められた周波数帯域のパワーの時間変化を算出するパワー算出部と、前記算出されたパワーが特定の条件を満たすか否かを判定し、判定結果に応じて心房細動の有無に関する情報を出力する判定部として機能させるためのプログラムを提供する。
このプログラムによれば、脈波信号および心電図の波形信号などのRR間隔が計測可能な信号において体動ノイズの影響が含まれていても、その信号から心房細動を判定することができる。
Further, the present invention provides a computer based on an acquisition unit that acquires a detection waveform signal indicating a detection result of an electrocardiogram or a pulse wave, and a spectrum of each frame obtained by frequency analysis of the acquired detection waveform signal. A parameter corresponding to the RR interval is calculated for each frame, an RR interval calculation unit that calculates an RR waveform signal indicating a time change of the parameter, and a time change of power in a predetermined frequency band in the RR waveform signal is calculated. And a program for determining whether or not the calculated power satisfies a specific condition and outputting information regarding the presence or absence of atrial fibrillation according to the determination result. .
According to this program, even if the influence of body motion noise is included in a signal capable of measuring an RR interval such as a pulse wave signal and an electrocardiogram waveform signal, atrial fibrillation can be determined from the signal.

実施形態における脈波測定装置の外観を説明する図である。It is a figure explaining the external appearance of the pulse wave measuring device in an embodiment. 実施形態における脈波測定装置の構成を説明する図である。It is a figure explaining the composition of the pulse wave measuring device in an embodiment. 実施形態における心房細動判定装置の機能構成を説明する図である。It is a figure explaining the functional structure of the atrial fibrillation determination apparatus in embodiment. 検出波形信号の周波数解析を行うときのフレームを説明する図である。It is a figure explaining a frame when performing frequency analysis of a detection waveform signal. RR間隔算出部における周波数特性を説明する図である。It is a figure explaining the frequency characteristic in a RR space | interval calculation part. RR波形信号の周波数解析を行うときのフレームを説明する図である。It is a figure explaining the frame when performing the frequency analysis of a RR waveform signal. 心房細動として判定された期間を説明する図である。It is a figure explaining the period determined as atrial fibrillation. 実施形態における心房細動判定処理を説明するフローチャートである。It is a flowchart explaining the atrial fibrillation determination process in embodiment. 非特許文献1における技術を用いた心房細動の判定方法を説明する図である。It is a figure explaining the determination method of atrial fibrillation using the technique in a nonpatent literature 1. FIG. RR間隔算出部の処理を施してから非特許文献1における技術を用いた心房細動の判定方法に適用した場合を説明する図である。It is a figure explaining the case where it applies to the determination method of the atrial fibrillation using the technique in a nonpatent literature 1 after performing the process of a RR space | interval calculation part.

<実施形態>
[脈波測定装置1の構成]
図1は、実施形態における脈波測定装置1の外観を説明する図である。本発明の実施形態における脈波測定装置1は、図1(a)に示すように、ユーザーである検出対象者の手1000における手首部分(腕)に腕時計のようにして装着される装置本体10と、検出部位に装着されて脈波を検出する脈波検出部20とを有する。装置本体10と脈波検出部20とは、ケーブル30により接続されている。ケーブル30は、脈波検出部20から出力される脈波信号(以下、検出波形信号Lという)を装置本体10に供給する一方、装置本体10からの電力を脈波検出部20に供給する。
<Embodiment>
[Configuration of Pulse Wave Measuring Device 1]
Drawing 1 is a figure explaining the appearance of pulse wave measuring device 1 in an embodiment. As shown in FIG. 1A, a pulse wave measuring device 1 according to an embodiment of the present invention is a device main body 10 that is worn like a wristwatch on a wrist portion (arm) of a user's 1000 hand to be detected. And a pulse wave detector 20 that is attached to the detection site and detects a pulse wave. The apparatus main body 10 and the pulse wave detection unit 20 are connected by a cable 30. The cable 30 supplies a pulse wave signal (hereinafter referred to as a detection waveform signal L) output from the pulse wave detection unit 20 to the apparatus main body 10, and supplies power from the apparatus main body 10 to the pulse wave detection unit 20.

装置本体10には、リストバンド50が取り付けられている。装置本体10は、リストバンド50がユーザーの腕に巻き付けられることにより腕に装着される。装置本体10には、操作部14および表示部15が設けられている。操作部14は、ユーザーが脈波測定装置1に機能選択の指示などを入力するためのボタンスイッチなどの操作子である。操作部14には、表示部15上に設けられたタッチセンサーなどが含まれていてもよい。表示部15は、液晶ディスプレイ、有機ELディスプレイなどの表示デバイスである。   A wristband 50 is attached to the apparatus main body 10. The apparatus main body 10 is attached to the arm by the wristband 50 being wound around the user's arm. The apparatus main body 10 is provided with an operation unit 14 and a display unit 15. The operation unit 14 is an operator such as a button switch for the user to input a function selection instruction or the like to the pulse wave measurement device 1. The operation unit 14 may include a touch sensor provided on the display unit 15. The display unit 15 is a display device such as a liquid crystal display or an organic EL display.

図1(b)に示すように、脈波検出部20が装着される検出部位は、この例においては、手1000における人差し指の根元から第2指関節までの間の一部であるものとするが、脈波を検出できる部位であればどの部位であってもよい。脈波検出部20は、固定バンド40によって固定されることにより、検出部位に装着される。このとき、固定バンド40は、脈波検出部20を覆った状態であり、脈波検出部20の受光部には固定バンド40の外側からの光が到達しないように遮光する構成である。   As shown in FIG. 1B, in this example, the detection site where the pulse wave detection unit 20 is worn is a part from the base of the index finger to the second finger joint in the hand 1000. However, any site may be used as long as it can detect a pulse wave. The pulse wave detector 20 is attached to the detection site by being fixed by the fixing band 40. At this time, the fixed band 40 is in a state of covering the pulse wave detection unit 20 and is configured to shield the light from the outside of the fixed band 40 from reaching the light receiving unit of the pulse wave detection unit 20.

脈波検出部20は、以下のようにして脈波を検出し、検出結果を示す検出波形信号Lを出力する。脈波検出部20は、発光部(例えば、緑色LED(Light Emitting Diode))と受光部とを有する。脈波検出部20は、装置本体10からケーブル30を介して供給された電力に応じた光を発光部から照射する。脈波検出部20は、発光部からの光のうち、毛細血管中のヘモグロビンによって反射した光を受光部により受光し、受光レベルに応じた信号を、検出波形信号Lとして、ケーブル30を介して装置本体10に供給する。   The pulse wave detector 20 detects a pulse wave as follows and outputs a detection waveform signal L indicating the detection result. The pulse wave detection unit 20 includes a light emitting unit (for example, a green LED (Light Emitting Diode)) and a light receiving unit. The pulse wave detection unit 20 emits light corresponding to the power supplied from the apparatus main body 10 via the cable 30 from the light emitting unit. The pulse wave detection unit 20 receives light reflected by hemoglobin in the capillary blood vessel from the light emission unit, and receives a signal corresponding to the received light level as a detection waveform signal L via the cable 30. It is supplied to the apparatus main body 10.

図2は、実施形態における脈波測定装置1の構成を説明する図である。脈波測定装置1は、CPU(Central Processing Unit)11、RAM(Random Access Memory)12、ROM(Read Only Memory)13、操作部14、表示部15、発振回路16、計時回路17、A/D変換回路18、増幅回路19を有する装置本体10と、脈波検出部20とを有する。増幅回路19および脈波検出部20を除く各構成は、バスを介して接続されている。   FIG. 2 is a diagram illustrating the configuration of the pulse wave measurement device 1 according to the embodiment. The pulse wave measuring apparatus 1 includes a CPU (Central Processing Unit) 11, a RAM (Random Access Memory) 12, a ROM (Read Only Memory) 13, an operation unit 14, a display unit 15, an oscillation circuit 16, a timing circuit 17, and an A / D. The apparatus main body 10 having a conversion circuit 18 and an amplification circuit 19 and a pulse wave detection unit 20 are included. The components other than the amplifier circuit 19 and the pulse wave detector 20 are connected via a bus.

CPU11は、ROM13に記憶されている制御プログラムにしたがって、各部の制御およびデータの転送などを行う。RAM12は、、検出波形信号Lなどの生体情報、およびCPU11における制御プログラムの実行中に発生する各種データを一時記憶する。CPU11は制御プログラムを実行することにより、心房細動判定機能を実現し、脈波測定装置1を心房細動判定装置として機能させる。なお、CPU11は、制御プログラムを実行することにより、心房細動判定機能以外の様々な機能を実現するようにしてもよい。これらの機能は、例えば、ユーザーが操作部14を操作することにより実現されるようにすればよい。   The CPU 11 controls each unit and transfers data according to a control program stored in the ROM 13. The RAM 12 temporarily stores biological information such as the detected waveform signal L and various data generated during execution of the control program in the CPU 11. By executing the control program, the CPU 11 realizes an atrial fibrillation determination function and causes the pulse wave measurement device 1 to function as an atrial fibrillation determination device. The CPU 11 may realize various functions other than the atrial fibrillation determination function by executing a control program. These functions may be realized by, for example, the user operating the operation unit 14.

操作部14は、上述したように、ユーザーの指示を脈波測定装置1に入力するためのボタンスイッチなどを有する。操作部14は、ユーザーによって操作されると、CPU11に対して、操作内容を示す操作信号を出力する。
表示部15は、上述したように、液晶ディスプレイ、有機ELディスプレイなどの表示デバイスを有し、CPU11によって表示内容が制御される。この表示内容は、例えば、時刻表示、各種メニュー画面、脈波測定結果、心房細動の判定結果などを示す各種画像である。
発振回路16は、制御の基礎となるクロック信号をCPU11に供給する。
計時回路17は、CPU11の制御により、時間を計測する。
As described above, the operation unit 14 includes a button switch for inputting a user instruction to the pulse wave measuring device 1. When operated by the user, the operation unit 14 outputs an operation signal indicating the operation content to the CPU 11.
As described above, the display unit 15 includes a display device such as a liquid crystal display or an organic EL display, and the display content is controlled by the CPU 11. The display contents are various images showing time display, various menu screens, pulse wave measurement results, atrial fibrillation determination results, and the like.
The oscillation circuit 16 supplies the CPU 11 with a clock signal that is the basis of control.
The timer circuit 17 measures time under the control of the CPU 11.

増幅回路19は、脈波検出部20からケーブル30を介して供給された検出波形信号Lを増幅する。増幅時のゲインは、CPU11の制御により設定される。
A/D変換回路18は、増幅回路19において増幅されたアナログ信号の検出波形信号Lを、デジタル信号に変換する。この例においては、サンプリング周波数は100Hzであり、脈波から得られるRR間隔に比べて十分高い周波数となっている。また、この例においては、量子化は10ビットで行われる。なお、サンプリング周波数、量子化ビットについては、必要とする精度に応じて異なる値に決められていてもよい。
続いて、CPU11によって実現される心房細動判定装置の機能構成(心房細動判定機能)について説明する。
The amplification circuit 19 amplifies the detection waveform signal L supplied from the pulse wave detection unit 20 via the cable 30. The gain at the time of amplification is set by the control of the CPU 11.
The A / D conversion circuit 18 converts the detected waveform signal L of the analog signal amplified by the amplification circuit 19 into a digital signal. In this example, the sampling frequency is 100 Hz, which is sufficiently higher than the RR interval obtained from the pulse wave. In this example, quantization is performed with 10 bits. Note that the sampling frequency and the quantization bit may be set to different values depending on the required accuracy.
Subsequently, the functional configuration (atrial fibrillation determination function) of the atrial fibrillation determination apparatus realized by the CPU 11 will be described.

[機能構成]
図3は、実施形態における心房細動判定装置100の機能構成を説明する図である。心房細動判定装置100は、ノイズ低減部111、RR間隔算出部112、パワー算出部113、判定部114、および表示制御部115を有するとともに、各種データの記憶領域となる検出波形信号記憶領域121、RR波形信号記憶領域122、およびパワー波形信号記憶領域123の各機能構成により実現される。
[Function configuration]
FIG. 3 is a diagram illustrating a functional configuration of the atrial fibrillation determination device 100 according to the embodiment. The atrial fibrillation determination device 100 includes a noise reduction unit 111, an RR interval calculation unit 112, a power calculation unit 113, a determination unit 114, and a display control unit 115, and a detected waveform signal storage region 121 that is a storage region for various data. , The RR waveform signal storage area 122 and the power waveform signal storage area 123 are realized by each functional configuration.

検出波形信号記憶領域121は、A/D変換回路18によってデジタル信号に変換された検出波形信号Lが記憶されるRAM12上に設けられた領域である。
ノイズ低減部111は、検出波形信号記憶領域121に記憶された検出波形信号Lから、RR間隔に相当する周波数帯域以外の体動ノイズ成分を低減するフィルター処理を施して出力する。フィルター処理としては、例えば、ハイパスフィルター、バンドパスフィルター、適応フィルターなどによる処理である。ノイズ低減部111において体動ノイズ成分が低減された検出波形信号Lについては、一旦RAM12において記憶されるようにしてもよい。検出波形信号記憶領域121およびノイズ低減部111は、RR間隔算出部112において周波数解析に用いられる検出波形信号Lを取得する取得部として機能する。
The detection waveform signal storage area 121 is an area provided on the RAM 12 in which the detection waveform signal L converted into a digital signal by the A / D conversion circuit 18 is stored.
The noise reduction unit 111 performs a filtering process for reducing body motion noise components other than the frequency band corresponding to the RR interval from the detection waveform signal L stored in the detection waveform signal storage area 121 and outputs the result. The filter process is, for example, a process using a high pass filter, a band pass filter, an adaptive filter, or the like. The detected waveform signal L in which the body movement noise component is reduced in the noise reduction unit 111 may be temporarily stored in the RAM 12. The detection waveform signal storage area 121 and the noise reduction unit 111 function as an acquisition unit that acquires the detection waveform signal L used for frequency analysis in the RR interval calculation unit 112.

なお、この処理においては体動ノイズ成分が低減され、その影響が検出波形信号Lから減少することにはなるが、背景技術として示した技術(特許文献1、非特許文献1)において精密な心房細動の判定ができるほど、正確なRR間隔を計測することができるまでには至らない。   In this process, the body motion noise component is reduced, and the influence thereof is reduced from the detected waveform signal L. However, the precise atrium in the technique shown as the background art (Patent Document 1, Non-Patent Document 1). It is not possible to measure an accurate RR interval to the extent that fibrillation can be determined.

RR間隔算出部112は、ノイズ低減部111において体動ノイズ成分が低減された検出波形信号Lについて、サンプリング毎にフレームを切り出し、短時間での周波数解析(STFT(Short-Time Fourier transform)解析)により周波数スペクトルを算出する。そして、RR間隔算出部112は、算出した周波数スペクトルに基づいて、フレーム毎にRR間隔に相当するパラメータを算出し、このパラメータの時間変化を示すRR波形信号FRRを、RAM12上に設けられた領域のRR波形信号記憶領域122に記憶する。   The RR interval calculation unit 112 cuts out a frame for each sampling of the detected waveform signal L from which the body motion noise component has been reduced in the noise reduction unit 111, and performs frequency analysis in a short time (STFT (Short-Time Fourier transform) analysis). To calculate the frequency spectrum. Then, the RR interval calculation unit 112 calculates a parameter corresponding to the RR interval for each frame based on the calculated frequency spectrum, and an RR waveform signal FRR indicating a time change of the parameter is an area provided on the RAM 12. Is stored in the RR waveform signal storage area 122.

算出されるパラメータは、この例においては、フレーム内におけるRR間隔の平均を示す値(平均脈波RR間隔)であり、例えば、周波数スペクトルの最大ピークとなる周波数である。したがって、RR波形信号FRRは、平均脈波RR間隔の時間変化を示すことになる。RR間隔算出部112における処理により、ノイズ低減部111において体動ノイズが完全に除去されなくても、RR波形信号FRRに含まれる体動ノイズの影響を大幅に低減することができる。   In this example, the calculated parameter is a value (average pulse wave RR interval) indicating an average of the RR intervals in the frame, and is, for example, a frequency that is the maximum peak of the frequency spectrum. Therefore, the RR waveform signal FRR indicates a time change of the average pulse wave RR interval. Even if the body motion noise is not completely removed by the noise reduction unit 111, the influence of the body motion noise included in the RR waveform signal FRR can be significantly reduced by the processing in the RR interval calculation unit 112.

図4は、検出波形信号Lの周波数解析を行うときのフレームを説明する図である。図4に示す波形は検出波形信号Lの波形の例である。図4に示すように、各フレームの時間は、この例においては4秒であり、1秒毎にサンプリングされ周波数解析が行われる。すなわち、各フレームは1秒ずつずれて設定され、次のフレームと3秒間オーバーラップしている。このように、サンプリングタイミング、フレームが設定されているため、平均脈波RR間隔はRR間隔の4秒間における平均値であり、RR波形信号FRRは1秒毎の平均脈波RR間隔の変化を示すものとなる。   FIG. 4 is a diagram for explaining a frame when performing frequency analysis of the detected waveform signal L. The waveform shown in FIG. 4 is an example of the waveform of the detected waveform signal L. As shown in FIG. 4, the time of each frame is 4 seconds in this example, and sampling is performed every second to perform frequency analysis. That is, each frame is set to be shifted by 1 second and overlaps with the next frame for 3 seconds. Thus, since the sampling timing and the frame are set, the average pulse wave RR interval is an average value for 4 seconds of the RR interval, and the RR waveform signal FRR indicates a change in the average pulse wave RR interval per second. It will be a thing.

図5は、RR間隔算出部112における周波数特性を説明する図である。RR間隔算出部112において、上述のように設定したフレームで周波数解析を行うことは、移動平均処理における周波数特性が重畳することと同等である。図5に示す周波数特性は、フレームの時間の4秒に相当する0.25Hzおよびその整数倍の周波数において谷が生じ、山の頂点を結ぶような全体的な傾向として、高周波数になるほどレベルが低くなる、すなわち、マイナスの傾きを持つ周波数特性となっている。フレームの時間が長くなるほど、この傾きは急になる。一方、フレームの時間が短くなると、傾きが「0」に近づくことになるが、検出波形信号Lにおける体動ノイズ成分の残存が多くなってしまう。そのため、フレーム時間は1秒以上5秒以下、望ましくは2秒以上4秒以下とするとよい。   FIG. 5 is a diagram for explaining frequency characteristics in the RR interval calculation unit 112. In the RR interval calculation unit 112, performing the frequency analysis with the frame set as described above is equivalent to superimposing the frequency characteristics in the moving average process. In the frequency characteristics shown in FIG. 5, valleys occur at a frequency of 0.25 Hz corresponding to 4 seconds of the frame time and an integer multiple thereof, and as a general tendency to connect the peaks, the level increases as the frequency increases. The frequency characteristic becomes low, that is, has a negative slope. This slope becomes steeper as the frame time increases. On the other hand, when the frame time is shortened, the inclination approaches “0”, but the body motion noise component remains in the detected waveform signal L. Therefore, the frame time is 1 to 5 seconds, preferably 2 to 4 seconds.

ここで、心房細動判定装置100の機能構成の説明を中断し、体動ノイズがほとんど存在しない理想的な心電図の波形信号から得られたRR間隔(心電RR間隔)を用いて、背景技術として示した非特許文献1に開示された技術により心房細動を判定する場合について説明する。また、この説明と対比して、この心電図の波形信号に対して、RR間隔算出部112における処理を施した波形信号(平均心電RR間隔)を用いて、非特許文献1に開示された技術によっては心房細動を判定することができないことについて説明する。   Here, the description of the functional configuration of the atrial fibrillation determination device 100 is interrupted, and the RR interval (electrocardiogram RR interval) obtained from the waveform signal of an ideal electrocardiogram having almost no body motion noise is used. A case where atrial fibrillation is determined by the technique disclosed in Non-Patent Document 1 shown as will be described. Further, in contrast to this description, a technique disclosed in Non-Patent Document 1 using a waveform signal (average electrocardiogram RR interval) obtained by performing processing in the RR interval calculation unit 112 on the waveform signal of the electrocardiogram. Explains that atrial fibrillation cannot be determined in some cases.

図9は、非特許文献1における技術を用いた心房細動の判定方法を説明する図である。図9は、心電RR間隔の変動を示す波形信号について、480秒間を1フレームとし、そのフレームにおいて0.01Hzから0.2Hzの帯域で周波数解析を行い、ピーク周波数とパワーとを対数変換して表したグラフである。図9(a)は、心房細動を発症していないときの心電RR間隔を用いた場合を示し、図9(b)は、心房細動を発症しているときの心電RR間隔を用いた場合を示している。図中の直線は、プロットされたデータから得られる1次回帰直線を示す。これらのグラフから、1次回帰直線の傾きβと相関係数γとを算出すると、以下に示す結果となる。   FIG. 9 is a diagram for explaining a method for determining atrial fibrillation using the technique in Non-Patent Document 1. FIG. 9 shows a waveform signal indicating fluctuations in the electrocardiogram RR interval. One frame is 480 seconds, and frequency analysis is performed in a band from 0.01 Hz to 0.2 Hz in that frame, and the peak frequency and power are logarithmically converted. It is the graph which expressed. FIG. 9 (a) shows a case where an electrocardiogram RR interval when no atrial fibrillation is developed, and FIG. 9 (b) shows an electrocardiogram RR interval when atrial fibrillation is developed. The case where it is used is shown. The straight line in the figure represents a linear regression line obtained from the plotted data. When the slope β and the correlation coefficient γ of the primary regression line are calculated from these graphs, the following results are obtained.

図9(a)に示す心房細動を発症していない場合においては、γ=−0.72、β=−1.29となる。また、図9(b)に示す心房細動を発症している場合においては、γ=−0.07、β=−0.13となる。このように、心房細動を発症している場合には、相関がなくなって白色雑音状になり、また、傾きβが「0」に近づくことがわかる。
一方、図9に示す場合と同じ心電RR間隔の変動を示す波形信号について、上述したRR間隔算出部112における処理を施した波形信号(平均心電RR間隔)を用いて、同様の解析を行った場合について、図10を用いて説明する。
In the case where the atrial fibrillation shown in FIG. 9A has not occurred, γ = −0.72 and β = −1.29. In the case where atrial fibrillation shown in FIG. 9B has developed, γ = −0.07 and β = −0.13. Thus, when atrial fibrillation has developed, it becomes clear that the correlation is lost and white noise occurs, and the slope β approaches “0”.
On the other hand, the same analysis is performed using the waveform signal (average electrocardiogram RR interval) subjected to the processing in the above-described RR interval calculation unit 112 with respect to the waveform signal showing the same variation of the electrocardiogram RR interval as shown in FIG. The case where it performed is demonstrated using FIG.

図10は、RR間隔算出部112の処理を施してから非特許文献1における技術を用いた心房細動の判定方法に適用した場合を説明する図である。図10においても、平均心電RR間隔の変動を示す波形信号について、480秒間を1フレームとし、そのフレームにおいて0.01Hzから0.2Hzの帯域で周波数解析を行い、ピーク周波数とパワーとを対数変換して表したグラフである。図10(a)は、心房細動を発症していないときの心電RR間隔を用いた場合を示し、図10(b)は、心房細動を発症しているときの心電RR間隔を用いた場合を示している。図中の直線は、プロットされたデータから得られる1次回帰直線を示す。   FIG. 10 is a diagram illustrating a case where the processing of the RR interval calculation unit 112 is performed and then applied to the determination method for atrial fibrillation using the technique in Non-Patent Document 1. In FIG. 10 as well, the waveform signal indicating the variation of the average electrocardiogram RR interval is set to one frame for 480 seconds, and frequency analysis is performed in the band from 0.01 Hz to 0.2 Hz in that frame, and the peak frequency and power are logarithmically. It is the graph expressed by converting. FIG. 10 (a) shows a case where an electrocardiogram RR interval when no atrial fibrillation is developed, and FIG. 10 (b) shows an electrocardiogram RR interval when atrial fibrillation develops. The case where it is used is shown. The straight line in the figure represents a linear regression line obtained from the plotted data.

図10(a)に示す心房細動を発症していない場合においては、γ=−0.68、β=−1.40となる。また、図10(b)に示す心房細動を発症している場合においては、γ=−0.41、β=−1.02となる。このように、非特許文献1に示す技術を適用した場合、心電RR間隔を用いたときには、図9に示すように心房細動の有無による違いが明確である。一方で、体動ノイズの影響を低減するために、平均心電RR間隔を用いたときには、図10に示すように、心房細動の発症の有無によるγ、βの有意差がなく、このデータからは心房細動の有無の判定が困難である。   In the case where atrial fibrillation shown in FIG. 10A has not developed, γ = −0.68 and β = −1.40. In the case where atrial fibrillation shown in FIG. 10B has developed, γ = −0.41 and β = −1.02. Thus, when the technique shown in Non-Patent Document 1 is applied, when the electrocardiogram RR interval is used, the difference depending on the presence or absence of atrial fibrillation is clear as shown in FIG. On the other hand, when the average electrocardiogram RR interval is used to reduce the influence of body motion noise, there is no significant difference in γ and β depending on the presence or absence of the onset of atrial fibrillation as shown in FIG. Therefore, it is difficult to determine the presence or absence of atrial fibrillation.

これは、上述した図5に示す周波数特性が重畳して現れ、傾きβがマイナス側に変化することが原因である。したがって、体動ノイズの影響を低減するために、RR間隔算出部112における処理を行うと、非特許文献1における技術を適用して心房細動の有無を判定することはできない。
一方で、心房細動が発症すると、図9(b)に示すように相関がなくなって白色雑音状になるが、図9(a)と図9(b)とを比較したときに、高周波数帯域側においてパワーが増大しているとも言える。本発明においては、このパワーの増大を指標として用いる構成である。
This is because the frequency characteristics shown in FIG. 5 described above are superimposed and the slope β changes to the negative side. Therefore, if processing in the RR interval calculation unit 112 is performed in order to reduce the influence of body motion noise, it is not possible to determine the presence or absence of atrial fibrillation by applying the technique in Non-Patent Document 1.
On the other hand, when atrial fibrillation develops, the correlation disappears and white noise appears as shown in FIG. 9B. However, when comparing FIG. 9A and FIG. It can be said that the power is increasing on the band side. In the present invention, this increase in power is used as an index.

例えば、図9、図10に示す例において、0.2Hz近傍の周波数帯域についてパワーを比較すると、以下のようになる。まず、心電RR間隔を用いた場合、心房細動を発症していないときのパワーが「1.59」、心房細動を発症しているときのパワーが「4.97」となる。また、平均心電RR間隔を用いた場合、心房細動を発症していないときのパワーが「0.05」、心房細動を発症しているときのパワーが「0.30」となる。このように、心房細動を発症しているときには、発症していないときに比べて、この周波数帯域のパワーが数倍に増加して有意差が見られる。したがって、この周波数帯域におけるパワーの変化は、心房細動の有無を判定する指標となり得る。   For example, in the example shown in FIGS. 9 and 10, the power is compared in the frequency band near 0.2 Hz as follows. First, when the electrocardiogram RR interval is used, the power when no atrial fibrillation is developed is “1.59”, and the power when atrial fibrillation is developed is “4.97”. When the average electrocardiogram RR interval is used, the power when atrial fibrillation does not develop is “0.05”, and the power when atrial fibrillation is developed is “0.30”. Thus, when atrial fibrillation has developed, the power in this frequency band increases several times and a significant difference is seen compared to when it does not. Therefore, the change in power in this frequency band can be an index for determining the presence or absence of atrial fibrillation.

図3に戻って説明を続ける。パワー算出部113は、RR波形信号記憶領域122に記憶されたRR波形信号FRRについて、短時間での周波数解析(STFT解析)を行い、得られる周波数スペクトルに基づいて、一部の周波数帯域(以下、算出周波数帯域という)のパワー(以下、帯域パワーという)を算出する。パワー算出部113は、算出した帯域パワーの時間変化を示すパワー波形信号Paを、RAM12上に設けられた領域のパワー波形信号記憶領域123に記憶する。   Returning to FIG. 3, the description will be continued. The power calculation unit 113 performs a short-time frequency analysis (STFT analysis) on the RR waveform signal FRR stored in the RR waveform signal storage area 122, and based on the obtained frequency spectrum, a part of frequency bands (hereinafter referred to as a frequency band). , Calculated frequency band) (hereinafter referred to as band power). The power calculation unit 113 stores the power waveform signal Pa indicating the time variation of the calculated band power in the power waveform signal storage area 123 of the area provided on the RAM 12.

図6は、RR波形信号FRRの周波数解析を行うときのフレームを説明する図である。図6に示す波形はRR波形信号FRRの波形の例である。図6に示すように、各フレームの時間は、この例においては120秒であり、60秒毎にサンプリングされ周波数解析が行われる。すなわち、各フレームは60秒ずつずれて設定され、次のフレームと60秒間オーバーラップしている。   FIG. 6 is a diagram for explaining a frame when performing frequency analysis of the RR waveform signal FRR. The waveform shown in FIG. 6 is an example of the waveform of the RR waveform signal FRR. As shown in FIG. 6, the time of each frame is 120 seconds in this example, and sampling is performed every 60 seconds for frequency analysis. That is, each frame is set so as to be shifted by 60 seconds and overlaps with the next frame for 60 seconds.

また、パワー算出部113において帯域パワーが算出される上述の算出周波数帯域は、予め決められ、この例においては、0.25Hzから0.5Hzまでの帯域であるものとする。これは、図5に示す周波数特性の2つの谷(0.25Hz、0.5Hzの谷)の間として決められている。これは、谷の部分におけるパワーは抑えられてしまうことから心房細動の有無の判定にはほとんど寄与しないため、心房細動の有無の判定に寄与する部分を中心に算出周波数帯域が決められている。すなわち、0.3Hzから0.45Hzまでの帯域といったように、算出周波数帯域は、さらに、周波数特性の谷の部分が除かれて、山の部分だけが含まれるように、その範囲が狭くなるように設定されていてもよい。   In addition, the above-described calculated frequency band in which the band power is calculated in the power calculation unit 113 is determined in advance, and in this example, is a band from 0.25 Hz to 0.5 Hz. This is determined as between two valleys (0.25 Hz and 0.5 Hz valleys) of the frequency characteristics shown in FIG. This is because the power in the valley portion is suppressed, so it hardly contributes to the determination of the presence or absence of atrial fibrillation, so the calculated frequency band is determined around the portion that contributes to the determination of the presence or absence of atrial fibrillation. Yes. That is, the range of the calculated frequency band, such as the band from 0.3 Hz to 0.45 Hz, is further narrowed so that only the mountain part is included by excluding the valley part of the frequency characteristics. May be set.

ここで、算出周波数帯域の最低周波数(下限)と最高周波数(上限)とは、この例においては、RR間隔算出部112における周波数特性、すなわちRR間隔算出部112における周波数解析で用いるフレームの時間に応じて決められていた。一方で、上下限の周波数のいずれか一方、または双方の周波数は、必ずしもフレームの時間に応じて決められていなくてもよい。   Here, the lowest frequency (lower limit) and the highest frequency (upper limit) of the calculated frequency band are the frequency characteristics in the RR interval calculation unit 112, that is, the frame time used in the frequency analysis in the RR interval calculation unit 112 in this example. It was decided accordingly. On the other hand, one or both of the upper and lower limit frequencies may not necessarily be determined according to the frame time.

算出周波数帯域の最低周波数は、図9、図10に示すように、パワーの変化が明確となる0.1Hz以上、望ましくは0.2Hz以上であるとよい。このとき、最低周波数は、上述したように、RR間隔算出部112における周波数解析で用いるフレームの時間の逆数以上であることがさらに望ましい。
また、算出周波数帯域の最高周波数は、ナイキスト周波数による影響を考え、RR間隔算出部112における周波数解析のサンプリング周波数の1/2以下であることが望ましい。このとき、最高周波数は、上述したように、RR間隔算出部112における周波数解析で用いるフレームの時間の逆数の2倍以下であることがさらに望ましい。
As shown in FIGS. 9 and 10, the lowest frequency in the calculated frequency band is 0.1 Hz or higher, preferably 0.2 Hz or higher, in which the power change becomes clear. At this time, as described above, it is more desirable that the minimum frequency is equal to or greater than the reciprocal of the frame time used in the frequency analysis in the RR interval calculation unit 112.
The maximum frequency in the calculated frequency band is preferably less than or equal to ½ of the sampling frequency of the frequency analysis in the RR interval calculation unit 112 in consideration of the influence of the Nyquist frequency. At this time, as described above, it is more desirable that the maximum frequency is not more than twice the reciprocal of the frame time used in the frequency analysis in the RR interval calculation unit 112.

図3に戻って説明を続ける。判定部114は、パワー波形信号記憶領域123に記憶されたパワー波形信号Paについて、特定の判定条件を満たすか否かを判定し、判定結果に応じた情報を出力する。特定の判定条件とは、この例においては、判定時点から過去30分の観測時間に対して、パワー波形信号Paの値(帯域パワー)が、予め決められた閾値Pth(例えば「1」)を超えた時間が50%以上の割合となることである。すなわち、この特定の判定条件を満たす場合には、その判定時点において心房細動であると判定する。   Returning to FIG. 3, the description will be continued. The determination unit 114 determines whether or not a specific determination condition is satisfied for the power waveform signal Pa stored in the power waveform signal storage area 123, and outputs information according to the determination result. In this example, the specific determination condition is that the value (band power) of the power waveform signal Pa is a predetermined threshold Pth (for example, “1”) for the observation time of the past 30 minutes from the determination time. The excess time is 50% or more. That is, when this specific determination condition is satisfied, it is determined that atrial fibrillation is present at the determination time.

なお、特定の判定条件は、様々に設定可能であり、パワー波形信号Paの値が一定時間にわたって閾値Pthを超えた場合に心房細動であると判定されるようにしてもよい。また、30分の観測時間に対して、閾値Pthを超えた時間が50%未満であっても、閾値Pthを超えた回数が一定回数を超えていることを条件として、心房細動であると判定されるようにしてもよいし、30分間の帯域パワーの積分値が一定の値を超えていれば、心房細動であると判定されるようにしてもよい。
また、リアルタイムの判定で無ければ、上記の過去30分の観測時間ではなく、例えば、判定時点の前後15分ずつの観測時間としてもよい。
The specific determination condition can be set in various ways, and may be determined as atrial fibrillation when the value of the power waveform signal Pa exceeds the threshold value Pth for a certain period of time. Moreover, even if the time exceeding the threshold Pth is less than 50% with respect to the observation time of 30 minutes, the number of times exceeding the threshold Pth exceeds a certain number of times, so that it is atrial fibrillation. It may be determined, or if the integrated value of the band power for 30 minutes exceeds a certain value, it may be determined as atrial fibrillation.
If the determination is not in real time, the observation time may be 15 minutes before and after the determination time instead of the observation time of the past 30 minutes.

図7は、心房細動として判定された期間を説明する図である。図7に示す波形は、パワー波形信号Paの波形の例であり、24時間にわたって脈波を検出して得られた波形である。期間Ta、Tb、Tcは、判定部114によって心房細動であると判定された期間を示している。この期間Ta、Tb、Tcにおいては、図7に示すように、パワー波形信号Paの値が、大きくなっている。   FIG. 7 is a diagram for explaining a period determined as atrial fibrillation. The waveform shown in FIG. 7 is an example of the waveform of the power waveform signal Pa, and is a waveform obtained by detecting a pulse wave over 24 hours. Periods Ta, Tb, and Tc indicate periods determined by the determination unit 114 as atrial fibrillation. In the periods Ta, Tb, and Tc, as shown in FIG. 7, the value of the power waveform signal Pa is large.

図3に戻って説明を続ける。判定部114は、心房細動であると判定すると、その判定結果を示す情報を表示制御部115に出力する。判定部114から出力される情報は、例えば、心房細動であることの判定を示す情報など、心房細動の有無に関する情報であればよい。情報表示制御部115は、判定部114から出力された情報に基づいて、表示部15の表示内容を制御し、心房細動であると判定されたことを示す画像を表示させる。ユーザーは、この表示内容を見ることにより、心房細動であるか否かを確認することができる。なお、この表示内容は、リアルタイムに心房細動の判定結果を示す表示であってもよいし、心房細動であると判定された期間を示す表示であってもよい。
以上が、心房細動判定装置100の機能構成の説明である。続いて、心房細動判定装置100の動作(心房細動判定処理)について、図8を用いて説明する。
Returning to FIG. 3, the description will be continued. If the determination unit 114 determines that it is atrial fibrillation, the determination unit 114 outputs information indicating the determination result to the display control unit 115. The information output from the determination unit 114 may be information regarding the presence or absence of atrial fibrillation, such as information indicating determination of atrial fibrillation, for example. The information display control unit 115 controls the display content of the display unit 15 based on the information output from the determination unit 114, and displays an image indicating that it is determined as atrial fibrillation. The user can confirm whether or not it is atrial fibrillation by viewing this display content. The display content may be a display indicating the determination result of atrial fibrillation in real time or a display indicating a period determined to be atrial fibrillation.
The functional configuration of the atrial fibrillation determination device 100 has been described above. Next, the operation (atrial fibrillation determination process) of the atrial fibrillation determination apparatus 100 will be described with reference to FIG.

[心房細動判定処理]
図8は、実施形態における心房細動判定処理を説明するフローチャートである。まずユーザーが操作部14を操作し、心房細動の判定処理を開始する指示が入力されると、CPU11は、図8に示すフローを開始する。CPU11は、ユーザーが操作部14を操作して、判定処理を終了する指示が入力されているか否かを判定する(ステップS110)。CPU11は、判定処理を終了する指示が入力されている場合(ステップS110;Yes)には、心房細動の判定処理を終了する。
[Atrial fibrillation determination processing]
FIG. 8 is a flowchart illustrating atrial fibrillation determination processing in the embodiment. First, when the user operates the operation unit 14 and an instruction to start the determination process of atrial fibrillation is input, the CPU 11 starts the flow shown in FIG. The CPU 11 determines whether or not the user operates the operation unit 14 to input an instruction to end the determination process (step S110). CPU11 complete | finishes the determination process of atrial fibrillation, when the instruction | indication which complete | finishes the determination process is input (step S110; Yes).

CPU11は、判定処理を終了する指示が入力されていない場合(ステップS110;No)には、脈波検出部20において脈波を検出させて検出波形信号Lを測定し(ステップS120)、ノイズ低減部111により体動ノイズ低減処理を行う(ステップS130)。このとき、CPU11は、検出波形信号LをRAM12の検出波形信号記憶領域121に記憶に記憶するが、体動ノイズ低減処理を行った検出波形信号Lを記憶するようにしてもよい。   When the instruction to end the determination process is not input (step S110; No), the CPU 11 causes the pulse wave detector 20 to detect the pulse wave and measure the detected waveform signal L (step S120), thereby reducing noise. The body motion noise reduction process is performed by the unit 111 (step S130). At this time, the CPU 11 stores the detection waveform signal L in the detection waveform signal storage area 121 of the RAM 12, but may store the detection waveform signal L subjected to the body movement noise reduction processing.

CPU11は、体動ノイズ低減処理を行った波形信号がRAM12に1フレーム蓄積されたか否かを判定する(ステップS140)。CPU110は、1フレームの蓄積がされていない場合(ステップS140;No)には、ステップS110に戻って処理を続ける。一方、1フレームの蓄積がされた場合(ステップS140;Yes)には、CPU11は、RR間隔算出部112により平均脈波RR間隔を算出する(ステップS210)。   The CPU 11 determines whether or not the waveform signal subjected to the body movement noise reduction process is stored in the RAM 12 for one frame (step S140). If one frame has not been accumulated (step S140; No), the CPU 110 returns to step S110 and continues processing. On the other hand, when one frame is accumulated (step S140; Yes), the CPU 11 calculates the average pulse wave RR interval by the RR interval calculation unit 112 (step S210).

CPU11は、RR間隔算出部112により算出した平均脈波RR間隔をRR波形信号記憶領域122に記憶する(ステップS220)。この記憶領域に記憶された平均脈波RR間隔の時間変化はRR波形信号FRRとなる。
CPU11は、RR波形信号記憶領域122に記憶されたRR波形信号FRRが1フレーム蓄積されたか否かを判定する(ステップS230)。CPU110は、1フレームの蓄積がされていない場合(ステップS230;No)には、ステップS110に戻って処理を続ける。一方、1フレームの蓄積がされた場合(ステップS230;Yes)には、CPU11は、パワー算出部113により帯域パワーを算出する(ステップS310)。
The CPU 11 stores the average pulse wave RR interval calculated by the RR interval calculation unit 112 in the RR waveform signal storage area 122 (step S220). The time change of the average pulse wave RR interval stored in this storage area becomes the RR waveform signal FRR.
The CPU 11 determines whether or not one frame of the RR waveform signal FRR stored in the RR waveform signal storage area 122 has been accumulated (step S230). When one frame is not accumulated (step S230; No), the CPU 110 returns to step S110 and continues the process. On the other hand, when one frame is accumulated (step S230; Yes), the CPU 11 calculates the band power by the power calculation unit 113 (step S310).

CPU11は、パワー算出部113により算出した帯域パワーをパワー波形信号記憶領域123に記憶する(ステップS320)。この記憶領域に記憶された帯域パワーの時間変化はパワー波形信号Paとなる。
CPU11は、記憶されているパワー波形信号Paを参照し、判定部114により、過去30分における帯域パワーの50%以上が閾値Pthを超えていることという心房細動の判定条件を満たすか否かを判定する(ステップS410)。CPU11は、判定条件を満たしていないと判定した場合(ステップS410;No)には、ステップ110に戻って処理を続ける。
一方、CPU11は、判定条件を満たしていると判定した場合(ステップS410;Yes)には、表示制御部115により心房細動であることの判定結果を表示部15に表示させ(ステップS420)、ステップS110に戻って処理を続ける。
The CPU 11 stores the band power calculated by the power calculation unit 113 in the power waveform signal storage area 123 (step S320). The time change of the band power stored in this storage area becomes the power waveform signal Pa.
The CPU 11 refers to the stored power waveform signal Pa, and determines whether or not the determination unit 114 satisfies the determination condition of atrial fibrillation that 50% or more of the band power in the past 30 minutes exceeds the threshold value Pth. Is determined (step S410). When the CPU 11 determines that the determination condition is not satisfied (step S410; No), the CPU 11 returns to step 110 and continues the processing.
On the other hand, if the CPU 11 determines that the determination condition is satisfied (step S410; Yes), the display control unit 115 causes the display unit 15 to display a determination result indicating that atrial fibrillation is present (step S420). It returns to step S110 and continues a process.

なお、CPU11は、ステップS110からステップS140の処理をステップS140の判定にかかわらず繰り返し実行してもよい。この場合、CPU11は、ステップS140においてYesとなる度に、ステップS110からステップS140の処理と並行して、ステップS210以降の処理を実行するようにしてもよい。このとき、ステップS230におけるNoの場合、または、ステップS410におけるNoの場合には、並行して実行したステップS210以降の処理を終了するようにすればよい。
以上が、心房細動判定処理の説明である。
Note that the CPU 11 may repeatedly execute the processing from step S110 to step S140 regardless of the determination in step S140. In this case, whenever it becomes Yes in step S140, you may make it perform the process after step S210 in parallel with the process of step S110 to step S140. At this time, in the case of No in step S230 or in the case of No in step S410, the processes after step S210 executed in parallel may be terminated.
The above is the description of the atrial fibrillation determination process.

このように、本発明の実施形態における脈波測定装置1においては、1拍毎の脈波RR間隔の代わりに平均脈波RR間隔を計測することにより体動ノイズの影響を低減しつつ、心房細動の判定を行うことができる。   As described above, in the pulse wave measuring device 1 according to the embodiment of the present invention, the influence of body motion noise is reduced by measuring the average pulse wave RR interval instead of the pulse wave RR interval for each beat, and the atrium. Fibrosis can be determined.

<変形例>
以上、本発明の実施形態について説明したが、本発明は以下のように、さまざまな態様で実施可能である。
[変形例1]
上述した実施形態においては、検出波形信号Lは、脈波検出部20において脈波を検出した結果を示す信号であったが、心電を検出した結果として得られる波形信号であってもよい。すなわち、RR間隔に相当するパラメーターが取得可能な波形信号であればよい。
<Modification>
As mentioned above, although embodiment of this invention was described, this invention can be implemented in various aspects as follows.
[Modification 1]
In the embodiment described above, the detection waveform signal L is a signal indicating a result of detecting a pulse wave in the pulse wave detection unit 20, but may be a waveform signal obtained as a result of detecting an electrocardiogram. That is, any waveform signal that can acquire a parameter corresponding to the RR interval may be used.

[変形例2]
上述した実施形態においては、心房細動判定装置100の機能構成としてノイズ低減部111が設けられていたが、必ずしも設けられていなくてもよい。この場合には、RR間隔算出部112は、検出波形信号記憶領域121から周波数解析を行う検出波形信号Lを取得すればよい。
[Modification 2]
In the above-described embodiment, the noise reduction unit 111 is provided as a functional configuration of the atrial fibrillation determination device 100. However, the noise reduction unit 111 is not necessarily provided. In this case, the RR interval calculation unit 112 may acquire the detection waveform signal L for performing frequency analysis from the detection waveform signal storage area 121.

[変形例3]
上述した実施形態においては、心房細動判定装置100は、脈波測定装置1において実現されていたがが、パーソナルコンピューターなどの情報処理装置において実現されるようにしてもよい。この場合には、情報処理装置は、予め測定された検出波形信号Lを外部装置から取得して、検出波形信号記憶領域121に記憶するようにすればよい。そして、情報処理装置は、この検出波形信号Lを心房細動判定処理により解析して、心房細動の有無を判定すればよい。
[Modification 3]
In the embodiment described above, the atrial fibrillation determination device 100 is realized in the pulse wave measurement device 1, but may be realized in an information processing device such as a personal computer. In this case, the information processing apparatus may acquire the detection waveform signal L measured in advance from the external device and store it in the detection waveform signal storage area 121. Then, the information processing apparatus may analyze the detected waveform signal L by atrial fibrillation determination processing to determine the presence or absence of atrial fibrillation.

[変形例4]
上述した実施形態においては、装置本体10と脈波検出部20とは、ケーブル30により有線で接続されていたが、無線により接続されていてもよい。この場合には、装置本体10と脈波検出部20とは、脈波検出部20の制御に必要な制御信号および脈波検出部20において生成される検出波形信号Lなどの各種信号を無線通信によりやり取りすればよい。また、装置本体10および脈波検出部20のそれぞれにおいて電力を供給可能な電池などの構成を有するようにすればよい。
[Modification 4]
In the above-described embodiment, the apparatus main body 10 and the pulse wave detection unit 20 are connected by wire with the cable 30, but may be connected wirelessly. In this case, the apparatus main body 10 and the pulse wave detection unit 20 wirelessly communicate various signals such as a control signal necessary for controlling the pulse wave detection unit 20 and a detection waveform signal L generated in the pulse wave detection unit 20. You can communicate with each other. Moreover, what is necessary is just to make it have the structure of the battery etc. which can supply electric power in each of the apparatus main body 10 and the pulse-wave detection part 20.

[変形例5]
上述した実施形態においては、心房細動の判定結果が表示部15に表示され、ユーザーに報知されていたが、音、振動などにより報知されるようにしてもよい。例えば、音でユーザーに報知する場合には、スピーカー、および判定部114からの情報に基づいて、スピーカーの放音内容を制御する音制御部を設ければよい。また、振動でユーザーに報知する場合には、振動アクチュエーター、および判定部114からの情報に基づいて、振動アクチュエーターの振動内容を制御する振動制御部を設ければよい。このように、実施形態における表示制御部115および表示部15は、心房細動の判定結果に応じてユーザーに報知する報知部として概念することもできる。
[Modification 5]
In the embodiment described above, the determination result of atrial fibrillation is displayed on the display unit 15 and notified to the user, but may be notified by sound, vibration, or the like. For example, when notifying the user with sound, a sound control unit for controlling the sound emission content of the speaker may be provided based on information from the speaker and the determination unit 114. When notifying the user by vibration, a vibration control unit that controls the vibration content of the vibration actuator based on the information from the vibration actuator and the determination unit 114 may be provided. Thus, the display control unit 115 and the display unit 15 in the embodiment can also be conceptualized as a notification unit that notifies the user according to the determination result of atrial fibrillation.

[変形例6]
上述した実施形態における制御プログラムは、磁気記録媒体(磁気テープ、磁気ディスクなど)、光記録媒体(光ディスクなど)、光磁気記録媒体、半導体メモリなどのコンピュータ読み取り可能な記録媒体に記憶した状態で提供し得る。また、脈波測定装置1は、各プログラムをネットワーク経由でダウンロードしてもよい。
[Modification 6]
The control program in the above-described embodiment is provided in a state stored in a computer-readable recording medium such as a magnetic recording medium (magnetic tape, magnetic disk, etc.), an optical recording medium (optical disk, etc.), a magneto-optical recording medium, or a semiconductor memory. Can do. Moreover, the pulse wave measuring apparatus 1 may download each program via a network.

1…脈波測定装置、10…装置本体、11…CPU、12…RAM、13…ROM、14…操作部、15…表示部、16…発振回路、17…計時回路、18…A/D変換回路、19…増幅回路、20…脈波検出部、30…ケーブル、40…固定バンド、50…リストバンド、100…心房細動判定装置、111…ノイズ低減部、112…RR間隔算出部、113…パワー算出部、114…判定部、115…表示制御部、121…検出波形信号記憶領域、122…RR波形信号記憶領域、123…パワー波形信号記憶領域、1000…手 DESCRIPTION OF SYMBOLS 1 ... Pulse wave measuring device, 10 ... Apparatus main body, 11 ... CPU, 12 ... RAM, 13 ... ROM, 14 ... Operation part, 15 ... Display part, 16 ... Oscillation circuit, 17 ... Time-counting circuit, 18 ... A / D conversion Circuit: 19 ... Amplifier circuit, 20 ... Pulse wave detector, 30 ... Cable, 40 ... Fixed band, 50 ... Wristband, 100 ... Atrial fibrillation determination device, 111 ... Noise reduction unit, 112 ... RR interval calculation unit, 113 DESCRIPTION OF SYMBOLS ... Power calculation part, 114 ... Determination part, 115 ... Display control part, 121 ... Detection waveform signal storage area, 122 ... RR waveform signal storage area, 123 ... Power waveform signal storage area, 1000 ... Hand

Claims (6)

心電または脈波の検出結果を示す検出信号を取得する取得部と、
得された前記検出信号の周波数解析により得られるフレームのスペクトルに基づいて、RR間隔に相当するパラメーターを算出し、前記パラメーターの時間変化を示すRR波形信号を算出するRR間隔算出部と、
前記RR波形信号予め決められた周波数帯域におけるパワーの時間変化を算出するパワー算出部と、
出された前記パワーの時間変化が特定の条件を満たすか否かに応じて心房細動の有無に関する情報を出力する判定部と
を具備することを特徴とする心房細動判定装置。
An acquisition unit for acquiring a detection signal indicating a detection result of an electrocardiogram or a pulse wave;
Based on the spectrum of the frames obtained by the frequency analysis of the obtained have been the detection signal taken to calculate the parameter corresponding to the RR interval, and RR interval calculation unit for calculating the RR waveform signal indicating the time change of the parameter,
A power calculation section for calculating a time change of the power in a predetermined frequency band of the RR waveform signal,
Atrial fibrillation determination apparatus characterized by time variation of the calculated issued said power; and a determination unit for outputting information about the presence of atrial fibrillation in accordance with whether or not a particular condition is satisfied.
前記周波数帯域の最低周波数は、前記フレームの時間の逆数以上である
ことを特徴とする請求項1に記載の心房細動判定装置。
The atrial fibrillation determination apparatus according to claim 1, wherein the lowest frequency of the frequency band is equal to or greater than the reciprocal of the time of the frame.
検出対象者の前記心電または前記脈波を検出する検出部と、
前記判定部によって出力された情報に基づいてユーザーに報知する報知部と
を具備し、
前記取得部は、検出された結果に応じて得られる検出波形信号を取得する
ことを特徴とする請求項1または請求項2に記載の心房細動判定装置。
A detection unit for detecting the electrocardiogram or the pulse wave of the person to be detected;
A notification unit for notifying a user based on the information output by the determination unit;
The atrial fibrillation determination device according to claim 1, wherein the acquisition unit acquires a detection waveform signal obtained according to a detected result.
前記取得部は、前記検出波形信号に対して、体動ノイズ成分を低減するフィルター処理を施し、前記検出波形信号として出力するノイズ低減部を備える
ことを特徴とする請求項3に記載の心房細動判定装置。
The atrial cell according to claim 3, wherein the acquisition unit includes a noise reduction unit that performs a filtering process to reduce a body movement noise component on the detection waveform signal and outputs the detection waveform signal. Motion determination device.
心電または脈波の検出結果を示す検出信号を取得する取得ステップと、
得された前記検出信号の周波数解析により得られるフレームのスペクトルに基づいて、RR間隔に相当するパラメーターを算出し、前記パラメーターの時間変化を示すRR波形信号を算出するRR間隔算出ステップと、
前記RR波形信号予め決められた周波数帯域におけるパワーの時間変化を算出するパワー算出ステップと、
出された前記パワーの時間変化が特定の条件を満たすか否かに応じて心房細動の有無に関する情報を出力する判定ステップと
を具備することを特徴とする心房細動判定方法。
An acquisition step of acquiring a detection signal indicating a detection result of an electrocardiogram or a pulse wave;
Based on the spectrum of the frames obtained by the frequency analysis of the obtained have been the detection signal taken to calculate the parameter corresponding to the RR interval, and RR interval calculation step of calculating the RR waveform signal indicating the time change of the parameter,
A power calculating step of calculating a time change of the power in a predetermined frequency band of the RR waveform signal,
Atrial fibrillation determination method characterized in that the time variation of the calculated issued said power; and a determination step of outputting information about the presence of atrial fibrillation in accordance with whether or not a particular condition is satisfied.
コンピューター
心電または脈波の検出結果を示す検出信号を取得
得された前記検出信号の周波数解析により得られるフレームのスペクトルに基づいて、RR間隔に相当するパラメーターを算出し、前記パラメーターの時間変化を示すRR波形信号を算出
前記RR波形信号予め決められた周波数帯域におけるパワーの時間変化を算出する
出された前記パワーの時間変化が特定の条件を満たすか否かに応じて心房細動の有無に関する情報を出力させる
プログラム。
On your computer,
Acquires a detection signal indicating the detection result of the electrocardiogram or pulse wave,
Based on the spectrum of the frames obtained by the frequency analysis of the obtained have been the detection signal taken to calculate the parameter corresponding to the RR interval, and calculates the RR waveform signal indicating the time change of the parameter,
It calculates the time change of the power in a predetermined frequency band of the RR waveform signal,
Program time change of calculated issued said power causes outputs information on the presence or absence of atrial fibrillation in accordance with whether or not a particular condition is satisfied.
JP2011194904A 2011-09-07 2011-09-07 Atrial fibrillation determination device, atrial fibrillation determination method and program Active JP5760876B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2011194904A JP5760876B2 (en) 2011-09-07 2011-09-07 Atrial fibrillation determination device, atrial fibrillation determination method and program
US13/604,064 US20130060154A1 (en) 2011-09-07 2012-09-05 Atrial fibrillation decision apparatus, and method and program for deciding presence of atrial fibrillation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2011194904A JP5760876B2 (en) 2011-09-07 2011-09-07 Atrial fibrillation determination device, atrial fibrillation determination method and program

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2015116573A Division JP6098673B2 (en) 2015-06-09 2015-06-09 Atrial fibrillation determination device, operation method and program for atrial fibrillation determination device

Publications (3)

Publication Number Publication Date
JP2013055982A JP2013055982A (en) 2013-03-28
JP2013055982A5 JP2013055982A5 (en) 2014-09-25
JP5760876B2 true JP5760876B2 (en) 2015-08-12

Family

ID=48132439

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011194904A Active JP5760876B2 (en) 2011-09-07 2011-09-07 Atrial fibrillation determination device, atrial fibrillation determination method and program

Country Status (1)

Country Link
JP (1) JP5760876B2 (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6060563B2 (en) * 2012-08-24 2017-01-18 セイコーエプソン株式会社 Atrial fibrillation determination device, atrial fibrillation determination method and program
JP2014171589A (en) 2013-03-07 2014-09-22 Seiko Epson Corp Atrial fibrillation analyzation equipment and program
JP2014171660A (en) 2013-03-08 2014-09-22 Seiko Epson Corp Atrial fibrillation analyzation equipment, atrial fibrillation analysis system, atrial fibrillation analysis method and program
JP2016047093A (en) 2014-08-27 2016-04-07 セイコーエプソン株式会社 Biological information analysis system, biological information processing system and biological information analysis device
JP6150825B2 (en) 2015-02-05 2017-06-21 ユニオンツール株式会社 Atrial fibrillation detection system
CN106037702A (en) 2015-04-17 2016-10-26 精工爱普生株式会社 Biological information processing system, biological information processing device, terminal device, method for generating analysis result information, and biological information processing method
CN106037701A (en) 2015-04-17 2016-10-26 精工爱普生株式会社 Biological information processing system, biological information processing device, and method for generating analysis result information
CN117100255B (en) * 2023-10-25 2024-01-23 四川大学华西医院 Method for judging fall prevention based on neural network model and related products

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3319140B2 (en) * 1994-04-19 2002-08-26 いすゞ自動車株式会社 Heart rate variability waveform analysis method and apparatus
JPH08280637A (en) * 1995-04-19 1996-10-29 Nissan Motor Co Ltd Mental activity judging apparatus
WO1997038626A1 (en) * 1996-04-17 1997-10-23 Seiko Epson Corporation Arrhythmia detector
JP3729143B2 (en) * 2002-03-14 2005-12-21 株式会社デンソー Pulse wave measuring device
US7092751B2 (en) * 2003-09-23 2006-08-15 Instrumentarium Corp. Detection of atrial arrhythmia
JP4753881B2 (en) * 2004-01-16 2011-08-24 コンプメディクス リミテッド Apparatus and signal processing method for monitoring, detecting and classifying sleep disordered breathing derived from electrocardiogram
JP2009089883A (en) * 2007-10-09 2009-04-30 Konica Minolta Sensing Inc Atrial fibrillation detector, system and method

Also Published As

Publication number Publication date
JP2013055982A (en) 2013-03-28

Similar Documents

Publication Publication Date Title
JP5760876B2 (en) Atrial fibrillation determination device, atrial fibrillation determination method and program
US9504401B2 (en) Atrial fibrillation analyzer and program
JP4855721B2 (en) Blood pressure measurement device
US20130060154A1 (en) Atrial fibrillation decision apparatus, and method and program for deciding presence of atrial fibrillation
JP5098721B2 (en) Blood pressure measurement device, blood pressure derivation program, and blood pressure derivation method
WO2013132844A1 (en) Pulse monitor and program
JP6060563B2 (en) Atrial fibrillation determination device, atrial fibrillation determination method and program
JP2014171660A (en) Atrial fibrillation analyzation equipment, atrial fibrillation analysis system, atrial fibrillation analysis method and program
JP6679602B2 (en) Device for detecting heart rate and heart rate variability
JP4892616B2 (en) Biological information detection device
WO2007032226A1 (en) Heart rate meter and method for removing noise of heart beat waveform
JP6036122B2 (en) Pulse wave measuring device and program
KR20120113530A (en) Method and apparatus for detecting peak from biological signal
WO2018168805A1 (en) Blood pressure measurement device, method, and program
JP6522327B2 (en) Pulse wave analyzer
JP4459713B2 (en) Biological information detection device
US20190320926A1 (en) Heartbeat detection method and heartbeat detection device
JP5127526B2 (en) Pulse wave measuring device and pulse wave measuring method
JP5998516B2 (en) Pulsation detection device, electronic device and program
JP5582051B2 (en) Pulse wave measuring device and program
JP6098673B2 (en) Atrial fibrillation determination device, operation method and program for atrial fibrillation determination device
JP2014132934A (en) Measuring apparatus and measuring method
TWI655930B (en) Measuring device, measuring method and recording medium
JP2022160608A (en) Biological signal processing device and control method of the same
JP2013183845A (en) Pulsation detector, electronic device and program

Legal Events

Date Code Title Description
A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140813

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20140813

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150424

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20150525

R150 Certificate of patent (=grant) or registration of utility model

Ref document number: 5760876

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350