JP4568825B2 - Calculation method of parasympathetic nerve activity index - Google Patents

Calculation method of parasympathetic nerve activity index Download PDF

Info

Publication number
JP4568825B2
JP4568825B2 JP2004166421A JP2004166421A JP4568825B2 JP 4568825 B2 JP4568825 B2 JP 4568825B2 JP 2004166421 A JP2004166421 A JP 2004166421A JP 2004166421 A JP2004166421 A JP 2004166421A JP 4568825 B2 JP4568825 B2 JP 4568825B2
Authority
JP
Japan
Prior art keywords
heartbeat
frequency
nerve activity
time
parasympathetic nerve
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.)
Expired - Fee Related
Application number
JP2004166421A
Other languages
Japanese (ja)
Other versions
JP2005319256A (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.)
Toyama Prefecture
Kanazawa Institute of Technology (KIT)
Original Assignee
Toyama Prefecture
Kanazawa Institute of Technology (KIT)
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 Toyama Prefecture, Kanazawa Institute of Technology (KIT) filed Critical Toyama Prefecture
Priority to JP2004166421A priority Critical patent/JP4568825B2/en
Publication of JP2005319256A publication Critical patent/JP2005319256A/en
Application granted granted Critical
Publication of JP4568825B2 publication Critical patent/JP4568825B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Description

心臓に関わる自律神経系の活動評価のうち、特に副交感神経活動指標の算出方法に関する。 Among the activities evaluation of the autonomic nervous system related to the heart, about the particular method of calculating the parasympathetic activity index.

従来の副交感神経活動指標の算出方法は、身体が定常的な状態にある時に、心電図のR波のピーク時間間隔を求め、この間隔から心拍周波数を算出する。この心拍周波数のペクトル分析を行い、このスペクトルの高周波成分(HF成分)の多寡により副交感神経活動を評価している。この高周波成分とは、呼吸周波数に等しい周波数成分であり、呼吸性洞性不整脈による影響であることが知られている。副交感神経が活発になると、呼吸性洞性不整脈(RSA)が顕著になり、心拍周波数が呼吸の影響を強く受けるようになる。 In the conventional parasympathetic activity index calculation method, when the body is in a steady state, the peak time interval of the R wave of the electrocardiogram is obtained, and the heartbeat frequency is calculated from this interval. A spectrum analysis of the heartbeat frequency is performed, and parasympathetic nerve activity is evaluated based on the number of high-frequency components (HF components) of this spectrum. This high-frequency component is a frequency component equal to the respiratory frequency and is known to be an effect of respiratory sinus arrhythmia. When parasympathetic nerves become active, respiratory sinus arrhythmia (RSA) becomes prominent, and the heart rate frequency is strongly influenced by respiration.

これに対して、交感神経活動は、比較的ゆっくりと働くために、この評価には、長時間の安定的な状態が必要である。副交感神経は、素早く働く性質があるが、医学や生理学の分野では、この交感神経活動と共に評価されることが多いため、ほとんど、長時間の安定的な状態でしか議論されなかった。このため、短時間に変動する様子を調べるための副交感神経活動の定量的な指標の算出方法はなかった。 On the other hand, since sympathetic nerve activity works relatively slowly, this evaluation requires a stable state for a long time. Although parasympathetic nerves have the property of working quickly, in the field of medicine and physiology, since they are often evaluated together with this sympathetic nerve activity, they were mostly discussed only in a stable state for a long time. For this reason, there has been no method for calculating a quantitative index of parasympathetic nerve activity for examining changes in a short time.

副交感神経活動は、本来、素早く働く自律神経系の働きであり、その働きを調べるには、短い時間単位で評価する必要がある。しかしながら、これまでは、心拍周波数の揺らぎをスペクトル分析法で評価する方法しかなく、有効なスペクトル解析を行うために、比較的長い時間における平均値的な評価しか出来なかった。
この発明は、上記背景技術の問題点に鑑み点されたもので、一呼吸ごとに、副交感神経活動指標を求めてその評価を可能にする副交感神経活動指標の算出方法を提供することを目的とする
Parasympathetic nerve activity is inherently a function of the autonomic nervous system that works quickly, and in order to examine its function, it is necessary to evaluate it in short time units. However, until now, there was only a method for evaluating fluctuations in the heartbeat frequency by a spectrum analysis method, and only an average evaluation over a relatively long time could be performed in order to perform an effective spectrum analysis.
The present invention has been made in view of the problems of the background art described above, and an object thereof is to provide a method for calculating a parasympathetic nerve activity index that allows a parasympathetic nerve activity index to be obtained and evaluated for each breath. To do .

心電図から、一心拍毎にR波のピーク時間を読みとり、このR波のピーク時間から一心拍ごとの心拍周期を求め、この心拍周期から心拍周波数を計算する。このとき、心拍周波数は一心拍ごとに定まるため、時間的に飛び飛びの値となる。この飛び飛びの値の間を補間して、任意の時間における心拍周波数を求める。また、呼吸周期も同時に計測し、一呼吸周期に相当する時間帯を切り出して、この時間帯の心拍周波数の揺らぎとこの呼吸周期との関係を次のように調べる。
この切り出した時間帯は、そのときの一呼吸に要した時間で、その逆数は呼吸周波数である。この時間帯における心拍周波数の変動成分のうち、この呼吸周波数と等しい周波数成分のパワーを計算する。そして、計算されたパワーを一呼吸時間における副交感神経活動指標とする。
From the electrocardiogram, the peak time of the R wave is read for each heartbeat, the heartbeat cycle for each heartbeat is obtained from the peak time of the R wave, and the heartbeat frequency is calculated from this heartbeat cycle . At this time, since the heartbeat frequency is determined for each heartbeat, it becomes a jump value in time. By interpolating between the jump values, the heart rate frequency at an arbitrary time is obtained. The respiratory cycle is also measured at the same time, a time zone corresponding to one respiratory cycle is cut out, and the relationship between the fluctuation of the heart rate frequency in this time zone and this respiratory cycle is examined as follows.
This cut-out time zone is the time required for one breath at that time, and its reciprocal is the breathing frequency. Of the fluctuation components of the heart rate frequency in this time zone, the power of the frequency component equal to the respiration frequency is calculated. The calculated power is used as an index of parasympathetic nerve activity during one breathing time.

この発明の副交感神経活動指標の算出方法は、副交感神経活動は身体の特定の器官に一種のリラックス状態を与える働きがあることを利用して、身体のリラックス状態の評価を可能とするものである。また、この神経活動は迅速に働く特徴があり、この活動を一呼吸毎に評価できれば、健康器具等を使用したとき身体へのリラックス効果の評価が可能となる。この他、人にリラックス感を与えることを目的とした商品の効能について、副交感神経活動から見た評価が可能となる。 The parasympathetic nerve activity index calculation method according to the present invention makes it possible to evaluate the relaxed state of the body by utilizing that parasympathetic nerve activity has a function of giving a kind of relaxed state to a specific organ of the body. . In addition, the neural activity has quickly work features, if you can evaluate this activity every single breath, it is possible to evaluate the relaxing effect on the body when using health equipment and the like. In addition, it is possible to evaluate the efficacy of products aimed at giving people a sense of relaxation from the viewpoint of parasympathetic activity.

本発明は、一呼吸毎に副交感神経活動を評価することができることから、使用時の副交感神経活動の状況を計測することができる。また、これらのデータを蓄積することにより、副交感神経活動の年齢的な特性や、個人特有のパターンなどが明らかになり、これを踏まえた適切なサービスをフィードバックさせることができる。
具体的な例としては、呼吸開始時間を計測する装置器具を有し、R波のピーク時間を何らかの方法で計測できる装置を有する製品を開発することにより、一呼吸ごとに副交感神経活動を評価する機能を持たせることが出来る。副交感神経活動は、一種のリラックス状態を評価する指標でもあるので、このような製品を利用している時のリラックス状態を評価して、この時の心身の状態に適した身体刺激を与える機能を持つ製品を開発することができる。このとき、利用者の年齢や病歴などを事前に考慮して装置を調整する必要がある。
Since this invention can evaluate a parasympathetic nerve activity for every respiration, the condition of the parasympathetic nerve activity at the time of use can be measured. Also, by accumulating these data, the age-related characteristics of parasympathetic nerve activity, individual-specific patterns, etc. are clarified, and appropriate services based on this can be fed back.
As a specific example, parasympathetic nerve activity is evaluated for each breath by developing a product having a device that measures the respiratory start time and having a device that can measure the peak time of the R wave by some method. Can have a function. Parasympathetic nerve activity is also an index for evaluating a kind of relaxed state, so it has a function that evaluates the relaxed state when using such a product and gives physical stimulation suitable for the state of mind and body at this time. You can develop the products you have. At this time, it is necessary to adjust the apparatus in consideration of the age and medical history of the user in advance.

以下、本発明を実施例によって具体的に説明する。呼吸性洞性不整脈は吸気中に心拍周波数が多くなり、呼気中に心拍周波数が小さくなる現象である。
図1は心電図と呼吸の計測結果を示したものである。呼吸開始は、吸気からはじまり、グラフに表示された呼吸曲線の値が極小となった時間であり、吸気は次の呼吸曲線が極大となったところで終了する。この時点から呼気が開始され、次の曲線が極小値になる時刻に呼気が終了する。これが、1呼吸のサイクルである。
図1から分かるように、吸気時のR波は、ピーク同士の間隔が狭く、呼気時のR波は、ピーク同士の間隔は広くなっている。この現象が呼吸性洞性不整脈である。この現象を定量的に算出するために、R波のピーク同士の中点心拍周波数プロットする。心拍周波数は、R波のピーク同士の間隔から算出される。この心拍周波数は黒点で表したように飛び飛びの値となるが、図2のように、この黒点の間にスプライン補間を行うと、心拍周波数の変動を、連続的に表すことができる。
Hereinafter, the present invention will be specifically described by way of examples. Respiratory sinus arrhythmia is a phenomenon in which the heart rate frequency increases during inspiration and the heart rate frequency decreases during expiration.
FIG. 1 shows electrocardiogram and respiration measurement results. The start of respiration starts from inspiration and is the time when the value of the respiration curve displayed on the graph becomes a minimum, and inspiration ends when the next respiration curve becomes maximum. Expiration starts from this point, and expires at the time when the next curve reaches a minimum value. This is a cycle of one breath.
As can be seen from Figure 1, R-wave during inspiration, narrow interval between the peaks, R-wave during expiration, the interval between the peaks is wider. This phenomenon is respiratory sinus arrhythmia. In order to quantitatively calculate this phenomenon , the heart rate frequency is plotted at the midpoint between the peaks of the R wave. The heartbeat frequency is calculated from the interval between the peaks of the R wave. The heartbeat frequency is a jump value as represented by black dots. However, if spline interpolation is performed between the black dots as shown in FIG. 2, fluctuations in the heartbeat frequency can be represented continuously.

次に、この求めた心拍周波数と呼吸との相関性を次のように計算する。心拍周波数成分は一呼吸の間にも変化している。この変化している心拍周波数は、色々な周波数成分を持っている。このうち、この一呼吸に対応する呼吸周期の逆数を呼吸周波数として計算し、そのパワーを心拍周波数と呼吸の相関性を示す値とする。図3に示すように、このように計算された相関性を副交感神経活動の強さとする。 Next, the correlation between the obtained heartbeat frequency and respiration is calculated as follows. The heart rate frequency component changes during one breath. This changing heart rate frequency has various frequency components. Among them, the inverse of the breathing cycle that corresponds to the single breath was calculated as the breathing frequency, a value indicating the correlation between respiration and the power and the heartbeat frequency. As shown in FIG. 3, the correlation calculated in this way is defined as the strength of parasympathetic nerve activity.

今日、健康効果を持つ製品や、ストレス状態におかれた人に何らかの癒しをもたらす製品が注目されている。しかしながら、健康状態は、年齢や病歴などの個人の特性が関わる部分が大きい。また、ストレス状態に対する対処法の種類は、個人の特性のほか、その時の生理的な状態によって微妙に異なる。
この発明は、人間の生理学的な知見をデータベース化して、これとコンピューター制御技術とを組み合わせることにより、個人の特性や、使用するときの状態に合わせた個別的なサービスを提供できる可能性がある。
Today, products with health benefits and products that bring some kind of healing to people in stress are drawing attention. However, the health condition is largely related to individual characteristics such as age and medical history. Moreover, the kind of the coping method with respect to a stress state is slightly different according to the physiological state at that time besides an individual characteristic.
This invention has a possibility of providing individual services according to individual characteristics and conditions when used by creating a database of human physiological knowledge and combining it with computer control technology. .

呼吸性洞性不整脈(RSA)の内容を示すグラフである。呼吸性洞性不整脈とは、吸気時に、心拍が早くなり、呼気時に心拍が遅くなる生理現象のことである。 It is a graph which shows the content of respiratory sinus arrhythmia (RSA). Respiratory sinus arrhythmia is a physiological phenomenon in which the heart beats faster during inspiration and slows during exhalation. 心電図から、心拍周波数を計算する方法を示すグラフである。隣合うR波のピークから、心拍周期を算出して、このR波の中間地点の心拍周波数を計算する。これらの飛び飛びの値をスプライン補間して、任意の時間での心拍周波数を求める。 It is a graph which shows the method of calculating a heart rate frequency from an electrocardiogram. The heartbeat period is calculated from the peak of the adjacent R wave, and the heartbeat frequency at the midpoint of this R wave is calculated. These jump values are spline interpolated to obtain the heart rate frequency at an arbitrary time. 副交感神経を評価する方法を示すグラフである。心拍周波数と呼気量比曲線との関係から、呼吸性洞性不整脈を定量的に評価している。 It is a graph which shows the method of evaluating a parasympathetic nerve. Respiratory sinus arrhythmia is quantitatively evaluated from the relationship between heart rate frequency and expiratory volume ratio curve.

Claims (1)

一心拍毎の心拍周期を求め、この心拍周期から心拍周波数を算出し、この心拍周波数の値の時間軸方向の間を補間して、任意の時間における心拍周波数を求め、同時に呼吸周期も計測し、一呼吸周期に相当する時間帯毎に、この時間帯の前記補間した心拍周波数の値の変動を示す波形について、前記呼吸周期から求めた呼吸周波数と等しい周波数成分のパワーを計算し、計算されたパワーを一呼吸の時間帯における副交感神経活動指標とすることを特徴とする副交感神経活動指標の算出方法。 Obtain the heartbeat cycle for each heartbeat, calculate the heartbeat frequency from this heartbeat cycle, interpolate between the heartbeat frequency values in the time axis direction, find the heartbeat frequency at any time, and simultaneously measure the respiratory cycle For each time zone corresponding to one breathing cycle, the power of the frequency component equal to the breathing frequency obtained from the breathing cycle is calculated for the waveform indicating the fluctuation of the interpolated heartbeat frequency value in this time zone. The parasympathetic nerve activity index is calculated by using the measured power as a parasympathetic nerve activity index in a breathing time zone .
JP2004166421A 2004-05-07 2004-05-07 Calculation method of parasympathetic nerve activity index Expired - Fee Related JP4568825B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2004166421A JP4568825B2 (en) 2004-05-07 2004-05-07 Calculation method of parasympathetic nerve activity index

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2004166421A JP4568825B2 (en) 2004-05-07 2004-05-07 Calculation method of parasympathetic nerve activity index

Publications (2)

Publication Number Publication Date
JP2005319256A JP2005319256A (en) 2005-11-17
JP4568825B2 true JP4568825B2 (en) 2010-10-27

Family

ID=35466963

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004166421A Expired - Fee Related JP4568825B2 (en) 2004-05-07 2004-05-07 Calculation method of parasympathetic nerve activity index

Country Status (1)

Country Link
JP (1) JP4568825B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8852116B2 (en) 2012-05-08 2014-10-07 Gm3 Co., Ltd. Psychiatric symptom and psychiatric disorder onset risk evaluation method and evaluator using heart rate variability index
EP2862509A1 (en) 2013-10-16 2015-04-22 GM3 Co., Ltd. Psychiatric symptom and psychiatric disorder onset risk evaluation method and evaluator using heart rate variability index

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4357503B2 (en) 2006-06-28 2009-11-04 株式会社東芝 Biological information measuring device, biological information measuring method, and biological information measuring program
JP4993980B2 (en) * 2006-09-15 2012-08-08 Gac株式会社 Apparatus and method capable of outputting expiration time
DE102006060819A1 (en) * 2006-12-21 2008-07-03 Fresenius Medical Care Deutschland Gmbh Patient's respiration rate determining method, involves determining momentary respiration rates, where weights of individual respiration rates depend on difference between respective respiration rates and estimated value
US20110054277A1 (en) * 2008-05-09 2011-03-03 Koninklijke Philips Electronics N.V. Contactless respiration monitoring of a patient and optical sensor for a photoplethysmography measurement
US9414770B2 (en) 2010-12-29 2016-08-16 Biosense Webster (Israel) Ltd. Respiratory effect reduction in catheter position sensing
JP5929020B2 (en) * 2011-07-04 2016-06-01 株式会社豊田中央研究所 Consciousness state estimation device and program
JP6815069B2 (en) * 2015-07-28 2021-01-20 日本光電工業株式会社 A computer-readable storage medium in which a waveform analysis method, a waveform analysis device, a waveform analysis program, and a waveform analysis program are stored.
JP2020174756A (en) * 2019-04-15 2020-10-29 公益財団法人鉄道総合技術研究所 Autonomic nerve activity monitoring method and system thereof

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08117198A (en) * 1994-10-25 1996-05-14 Isuzu Motors Ltd Apparatus for decision of physiological and device therefor state

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08117198A (en) * 1994-10-25 1996-05-14 Isuzu Motors Ltd Apparatus for decision of physiological and device therefor state

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8852116B2 (en) 2012-05-08 2014-10-07 Gm3 Co., Ltd. Psychiatric symptom and psychiatric disorder onset risk evaluation method and evaluator using heart rate variability index
EP2862509A1 (en) 2013-10-16 2015-04-22 GM3 Co., Ltd. Psychiatric symptom and psychiatric disorder onset risk evaluation method and evaluator using heart rate variability index

Also Published As

Publication number Publication date
JP2005319256A (en) 2005-11-17

Similar Documents

Publication Publication Date Title
Pinna et al. Effect of paced breathing on ventilatory and cardiovascular variability parameters during short-term investigations of autonomic function
Guenette et al. Sex differences in exercise-induced diaphragmatic fatigue in endurance-trained athletes
Al-Ani et al. Changes in RR variability before and after endurance training measured by power spectral analysis and by the effect of isometric muscle contraction
Sarang et al. Effects of two yoga based relaxation techniques on heart rate variability (HRV).
EP1507474B1 (en) Procedure for deriving reliable information on respiratory activity from heart period measurement
EP1495715B1 (en) A method and apparatus based on combination of three phsysiological parameters for assessment of analgesia during anesthesia or sedation
Russell et al. Inclusion of a rest period in diaphragmatic breathing increases high frequency heart rate variability: Implications for behavioral therapy
Katayama et al. Inspiratory muscle fatigue increases sympathetic vasomotor outflow and blood pressure during submaximal exercise
EP1273265B1 (en) Monitoring a condition of a patient under anaesthesia or sedation
US20120065480A1 (en) Stress monitor system and method
CA2604741C (en) Methods and devices for relieving stress
Ahuja et al. GSR and HRV: its application in clinical diagnosis
US9486148B2 (en) Method and apparatus for estimating depth of anethesia by using ECG signal
US10980489B2 (en) Method and device for quantifying a respiratory sinus arrhythmia and use of said type of method or said type of device
Archiza et al. Temporal characteristics of exercise-induced diaphragmatic fatigue
JP4568825B2 (en) Calculation method of parasympathetic nerve activity index
Stankovski et al. Time-frequency methods and voluntary ramped-frequency breathing: a powerful combination for exploration of human neurophysiological mechanisms
Sakakibara et al. Efficacy of paced breathing at the low-frequency peak on heart rate variability and baroreflex sensitivity
Blasco-Lafarga et al. Cardiovascular and autonomic responses to a maximal exercise test in elite youngsters
Matos-Santos et al. Cardiovascular responses to resistance exercise performed with large and small muscle mass
De Abreu et al. Effects of inspiratory muscle-training intensity on cardiovascular control in amateur cyclists
Stergiopoulos et al. Second ventilatory threshold assessed by heart rate variability in a multiple shuttle run test
Fisher et al. A method for more accurate determination of resonance frequency of the cardiovascular system, and evaluation of a program to perform it
Martins de Abreu et al. Effects of inspiratory muscle-training intensity on cardiovascular control in amateur cyclists
Chua et al. Combined electrocardiogram and photoplethysmogram measurements as an indicator of objective sleepiness

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20070328

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20070424

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20070328

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20070424

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100323

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100524

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

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20100713

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

Free format text: PAYMENT UNTIL: 20130820

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees