JP5812265B2 - Autonomic nerve state evaluation system - Google Patents

Autonomic nerve state evaluation system Download PDF

Info

Publication number
JP5812265B2
JP5812265B2 JP2011158788A JP2011158788A JP5812265B2 JP 5812265 B2 JP5812265 B2 JP 5812265B2 JP 2011158788 A JP2011158788 A JP 2011158788A JP 2011158788 A JP2011158788 A JP 2011158788A JP 5812265 B2 JP5812265 B2 JP 5812265B2
Authority
JP
Japan
Prior art keywords
pitch frequency
frequency band
index value
voice data
cem
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
JP2011158788A
Other languages
Japanese (ja)
Other versions
JP2013022180A (en
Inventor
格一 塩見
格一 塩見
あすか 原田
あすか 原田
Original Assignee
国立研究開発法人 電子航法研究所
国立研究開発法人 電子航法研究所
格一 塩見
格一 塩見
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 国立研究開発法人 電子航法研究所, 国立研究開発法人 電子航法研究所, 格一 塩見, 格一 塩見 filed Critical 国立研究開発法人 電子航法研究所
Priority to JP2011158788A priority Critical patent/JP5812265B2/en
Publication of JP2013022180A publication Critical patent/JP2013022180A/en
Application granted granted Critical
Publication of JP5812265B2 publication Critical patent/JP5812265B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Description

この発明は、人の自律神経系の状態を評価するシステムに関する。   The present invention relates to a system for evaluating the state of a human autonomic nervous system.

人の自律神経の状態が、心拍や脈拍、発汗状態、眼球の動き、瞳孔の大きさなどに表れることは、従来から知られている。そして、心拍数、脈拍数、発汗状態、眼球の動き、瞳孔の大きさを検出して、その人の自律神経の状態を評価することが行われている(特許文献1〜3)。
また、自律神経のうち、交感神経が優位になったときには、心拍数や脈拍数、発汗量、眼球の動きが多くなり、反対に副交感神経が優位になったときにはこれらの値が少なくなることも知られている。
そして、交感神経と副交感神経のどちらが優位な状態になっているのかという自律神経系の状態を定量的な情報として評価する方法として、例えば、心拍変動スペクトル解析が行われ、医療分野で応用されている(非特許文献1)。
It has been conventionally known that the state of human autonomic nerves appears in heartbeat, pulse, sweating, eyeball movement, pupil size, and the like. Then, the heart rate, the pulse rate, the sweating state, the movement of the eyeball, and the size of the pupil are detected to evaluate the state of the person's autonomic nerve (Patent Documents 1 to 3).
Also, among the autonomic nerves, when the sympathetic nerve is dominant, the heart rate, the pulse rate, the amount of sweating, and the movement of the eyeball increase. On the other hand, when the parasympathetic nerve is dominant, these values may decrease. Are known.
And, as a method for evaluating the state of the autonomic nervous system, which is the dominant state of sympathetic nerve or parasympathetic nerve, as quantitative information, for example, heart rate variability spectrum analysis is performed and applied in the medical field (Non-Patent Document 1).

この方法は、被検者の心電図を測定し、そのR−R間隔を抽出して解析するものである。具体的には、上記R−Rの間隔のフーリエ変換によるパワースペクトルをとる。
このパワースペクトルにおいて、高周波成分は副交感神経活動のみを反映し、低周波数成分は副交感神活動と交感神経活動の双方を反映している。従って、これら高周波成分と低周波成分のスペクトルの大きさによって、副交感神経系と交感神経系のいずれが優位な状態なのかを知ることができる。
このように自律神経系の状態を知ることは、作業中の緊張や疲労度、特定の病気の患者の状態を推測するために利用される。
In this method, an electrocardiogram of a subject is measured, and the RR interval is extracted and analyzed. Specifically, a power spectrum is obtained by Fourier transform of the RR interval.
In this power spectrum, the high frequency component reflects only parasympathetic activity, and the low frequency component reflects both parasympathetic activity and sympathetic activity. Therefore, it is possible to know which of the parasympathetic nervous system and the sympathetic nervous system is dominant by the magnitude of the spectrum of these high-frequency components and low-frequency components.
Knowing the state of the autonomic nervous system in this way is used to estimate the tension and fatigue during work and the state of a patient with a specific disease.

特開平11−342120号公報JP 11-342120 A 特開2005−034470号公報JP 2005-034470 A 特開平07−171104号公報Japanese Patent Laid-Open No. 07-171104

谷 明博、鍵谷 俊文、堀 正二,“心拍変動のスペクトル解析による自律神経機能評価”、CLINICIAN、エーザイ株式会社、1992年12月1日、Vol.39 NO.416、P.57-63Akihiro Tani, Toshifumi Kagiya, Shoji Hori, “Evaluation of Autonomic Nervous Function by Spectral Analysis of Heart Rate Variability”, CLINICIAN, Eisai Co., Ltd., December 1, 1992, Vol.39 NO.416, P.57-63 塩見格一、「発話分析から考える脳機能モデル」、感性工学研究論文集、日本、日本感性工学学会、2004年2月、第4巻、1号、p.3−12Seiichi Shiomi, “Brain function model considered from speech analysis”, Kansei Engineering Research Papers, Japan, Japan Society for Kansei Engineering, February 2004, Vol. 4, No. 1, p. 3-12 M.Sano、Y.Sawada、“Measurement of theLyapunov Spectrum from a Chaotic Time Series”、The American Physical Society、アメリカ、1985年、Vol.55NO.10、P.1082-1085M.Sano, Y.Sawada, “Measurement of the Lyapunov Spectrum from a Chaotic Time Series”, The American Physical Society, USA, 1985, Vol. 55NO.10, P.1082-1085

上記心拍変動スペクトル解析を行ためには、被験者の心電図が必要である。そのため、被験者に心電図計を取り付けなければならない。
また、特許文献1、2に記載された測定装置も、被験者に皮膚抵抗や、血圧を計測する装置を取り付ける必要がある。
このように、被験者にその身体情報を計測するための計測装置を取り付ける必要がある。そのため、被験者が、何らかの作業をしているときには、上記計測装置が邪魔になる。そのため、作業中の自律神経系の状態を正確に測定することは難しかった。
また、上記心拍変動スペクト解析を行うためには、解析対照の心電データを取り込んでおく必要があり、被験者の作業中に解析結果を得ることも難しかった。
To cormorants row the heart rate variability spectrum analysis, it is necessary electrocardiogram of the subject. Therefore, an electrocardiograph must be attached to the subject.
In addition, the measuring devices described in Patent Documents 1 and 2 also need to attach a device for measuring skin resistance and blood pressure to the subject.
Thus, it is necessary to attach a measuring device for measuring the physical information to the subject. Therefore, when the subject is doing some work, the measurement device becomes an obstacle. Therefore, it is difficult to accurately measure the state of the autonomic nervous system during work.
Further, in order to perform the HRV spectrum analysis, must incorporate an electrocardiograph data analysis control was also difficult to obtain an analysis result while working subject.

一方、特許文献3の自律神経機能検査装置は非接触で瞳孔を測定するものであるが、瞳孔の大きさを光学的に正確に測定するためには、被験者に対する光の方向が重要である。この方法も、被験者が測定装置を意識せずに行っている作業中に、正確な測定をすることは困難である。
この発明の目的は、被験者に測定装置を意識させずに、作業中の自律神経系の状態を測定することができる自律神経の状態評価システムを提供することである。
On the other hand, the autonomic nervous function test apparatus of Patent Document 3 measures the pupil in a non-contact manner, but in order to optically accurately measure the size of the pupil, the direction of light with respect to the subject is important. Also in this method, it is difficult to perform an accurate measurement while the subject is performing work without being aware of the measuring device.
An object of the present invention is to provide an autonomic nerve state evaluation system capable of measuring the state of the autonomic nervous system during work without making the subject aware of the measuring device.

第1の発明は、音声入力手段と、この音声入力手段から入力した音声データをデジタル音声データに変換するA/D変換手段と、デジタル音声データのピッチ周波数帯域を取り出すフィルターと、このフィルターで取り出したピッチ周波数帯のデジタル音声データを処理し、脈拍変動率と相関性を有するカオス論的指標値を算出し、それを出力する信号処理手段とを備えた点を特徴とする。 The first invention is a voice input means, an A / D conversion means for converting voice data inputted from the voice input means into digital voice data, a filter for taking out the pitch frequency band of the digital voice data, and taking out by this filter. And a signal processing means for processing the digital audio data in the pitch frequency band, calculating a chaotic index value having a correlation with the pulse fluctuation rate, and outputting it.

第2の発明は、上記第1の発明を前提とし、上記デジタル音声データを記憶する音声データ記憶手段と、このデータ記憶手段が記憶しているデジタル音声データからピッチ周波数を算出するとともに、このピッチ周波数に基づいてピッチ周波数帯域を算出するピッチ周波数帯域算出手段とを備え、上記フィルターは、上記音声データ記憶手段が記憶しているデジタル音声データから、上記ピッチ周波数帯域算出手段が算出したピッチ周波数帯域の音声データを取り出す点を特徴とする。   The second invention is based on the first invention, and calculates the pitch frequency from the voice data storage means for storing the digital voice data, the digital voice data stored in the data storage means, and the pitch. Pitch frequency band calculating means for calculating a pitch frequency band based on the frequency, and the filter is a pitch frequency band calculated by the pitch frequency band calculating means from the digital audio data stored in the audio data storage means. It is characterized in that the voice data is extracted.

第1、第2の発明の自律神経の状態評価システムによれば、被験者の音声データを取得するだけで、そのときの自律神経の状態を測定することができる。このシステムは、従来の測定方法のように、他の作業の妨げになるような計測装置を身体に取り付ける必要がないので、通常の作業状態における自律神経の状態を測定できる。   According to the autonomic nerve state evaluation system of the first and second inventions, the state of the autonomic nerve at that time can be measured only by acquiring the voice data of the subject. This system can measure the state of the autonomic nerve in a normal working state because there is no need to attach a measuring device that interferes with other work to the body unlike the conventional measuring method.

第2の発明では、被験者ごとのピッチ周波数帯域を算出してフィルターを設定することができるので、正確にピッチ周波数帯域の音声データを取り出すことができる。その結果として、自律神経の状態を正確に測定することができる。
特に、ピッチ周波数帯域を取り出すフィルターを通す音声データに基づいて算出したピッチ周波数帯域を上記フィルターに設定するようにすれば、一人の被験者のピッチ周波数の変動にも対応してピッチ周波数帯域の音声データをより正確に取り出すことができるようになる。
In the second invention, the pitch frequency band for each subject can be calculated and the filter can be set, so that the voice data in the pitch frequency band can be accurately extracted. As a result, the state of the autonomic nerve can be accurately measured.
In particular, if the pitch frequency band calculated based on the voice data that passes through the filter for extracting the pitch frequency band is set in the filter, the voice data in the pitch frequency band corresponding to the variation of the pitch frequency of one subject. Can be taken out more accurately.

図1はこの発明の第1実施形態のシステム構成図である。FIG. 1 is a system configuration diagram of a first embodiment of the present invention. 図2は、第1実施形態のシステムを用いて算出した被験者の座位と立位におけるピッチ周波数帯域の音声データの指標値CEMの変動率と対応人数との関係を示したグラフである。FIG. 2 is a graph showing the relationship between the change rate of the index value CEM of the voice data in the pitch frequency band in the sitting position and the standing position of the subject calculated using the system of the first embodiment, and the corresponding number of persons. 図3は、第1実施形態のシステムを用いて算出した、車両運転中における被験者のピッチ周波数帯域の音声データのカオス論的指数値CEM及び心拍変動率と走行速度との関係を示したグラフである。FIG. 3 is a graph showing the relationship between the chaotic theoretical index value CEM and the heart rate variability and the running speed of the voice data in the pitch frequency band of the subject while driving the vehicle, calculated using the system of the first embodiment. is there. 図4はこの発明の第2実施形態のシステム構成図である。FIG. 4 is a system configuration diagram of the second embodiment of the present invention.

図1に示す第1実施形態は、自律神経の状態を評価すべき、被験者の声を収集するためのこの発明の音声入力手段であるマイクロフォン1と、このマイクロフォン1から入力されたアナログ音声データをデジタル変換するA/D変換器2と、このA/D変換器2でデジタル変換されたデジタル音声データを記憶するデータ記憶部6と、このデータ記憶部6に接続したピッチ周波数帯域通過フィルター3及びピッチ周波数算出部7と、上記ピッチ周波数通過フィルター3を通過したデータを処理するデータ処理部4と、データ処理部4の処理結果を出力する出力部5とを備えている。   In the first embodiment shown in FIG. 1, a microphone 1 which is a voice input means of the present invention for collecting a subject's voice for which the state of an autonomic nerve is to be evaluated, and analog voice data input from the microphone 1 are used. A / D converter 2 for digital conversion, a data storage unit 6 for storing digital audio data digitally converted by the A / D converter 2, a pitch frequency bandpass filter 3 connected to the data storage unit 6, and A pitch frequency calculation unit 7, a data processing unit 4 that processes data that has passed through the pitch frequency pass filter 3, and an output unit 5 that outputs a processing result of the data processing unit 4 are provided.

上記ピッチ周波数帯域通過フィルター3は、音声データのうちピッチ周波数帯域の音声データを取り出すためのフィルターである。ピッチ周波数とは、被験者の声帯の振動によって決まる周波数のことで、声帯から***までの声道の形状による影響を除いた情報である。
なお、上記ピッチ周波数算出部7は、マイクロフォン1から取り込んだ全帯域のデジタル音声データからピッチ周波数を算出し、そのピッチ周波数に基づいた帯域を上記ピッチ周波数帯域通過フィルター3に設定する機能を備えている。
The pitch frequency band pass filter 3 is a filter for extracting audio data in the pitch frequency band from the audio data. The pitch frequency is a frequency determined by vibration of the subject's vocal cords, and is information excluding the influence of the shape of the vocal tract from the vocal cords to the lips.
The pitch frequency calculation unit 7 has a function of calculating a pitch frequency from digital audio data of the entire band captured from the microphone 1 and setting a band based on the pitch frequency in the pitch frequency band pass filter 3. Yes.

この第1実施形態では、上記ピッチ周波数算出部7が上記データ記憶部6を介して入力したサンプリング時間当たりの音声データをケプストラム分析してピッチ周波数を算出し、そのピッチ周波数に基づいた帯域をピッチ周波数帯域通過フィルター3に設定する。上記ピッチ周波数帯域通過フィルター3は、ピッチ周波数算出部7によって設定されたピッチ周波数に対応したピッチ周波数帯域の音声データを通過させてデータ処理部4へ入力する。
なお、上記ピッチ周波数帯域の音声データを取り出すフィルターとしては、ピッチ周波数帯域のデータのみを通過させるものや、ピッチ周波数帯域以外のデータを減衰させるものなどどのような原理を利用したものでもよい。
In the first embodiment, the pitch frequency calculation unit 7 performs cepstrum analysis on the audio data per sampling time input via the data storage unit 6 to calculate the pitch frequency, and the band based on the pitch frequency is pitched. Set to frequency band pass filter 3. The pitch frequency band pass filter 3 passes audio data in a pitch frequency band corresponding to the pitch frequency set by the pitch frequency calculation unit 7 and inputs the audio data to the data processing unit 4.
The filter for extracting the audio data in the pitch frequency band may use any principle such as a filter that passes only the data in the pitch frequency band or a filter that attenuates data other than the pitch frequency band.

そして、この第1実施形態では、上記ピッチ周波数帯域通過フィルター3には、上記ケプストラム分析の対象とした音声データそのものを入力するようにしている。そのため、ピッチ周波数算出部7によってピッチ周波数を算出し、上記ピッチ周波数帯域通過フィルター3にピッチ周波数帯域が設定されるまでは、ピッチ周波数帯域フィルター3に音声データを入力しないように、データ記憶部6が音声データを保持している。   In the first embodiment, the pitch frequency bandpass filter 3 is inputted with the audio data itself as the object of the cepstrum analysis. Therefore, the data storage unit 6 prevents the voice data from being input to the pitch frequency band filter 3 until the pitch frequency is calculated by the pitch frequency calculation unit 7 and the pitch frequency band is set in the pitch frequency band pass filter 3. Holds voice data.

また、データ処理部4は、上記ピッチ周波数帯域通過フィルター3を通過した、すなわちピッチ周波数帯域の音声データを処理して、カオス論的指標値を算出する機能を備えている。
上記カオス論的指標値とは、カオス性を有する時系列データから作成したアトラクタの広がりを表す指標値である。例えば、非特許文献3に記載された佐野・澤田のアルゴリズムにより算出される最大リアプノフ指数や、非特許文献2で定義されたSiCECAアルゴリズムで算出されるマクロ脳活性度指標値などである。
Further, the data processing unit 4 has a function of calculating the chaotic index value by processing the voice data that has passed through the pitch frequency band pass filter 3, that is, the pitch frequency band.
The chaotic index value is an index value representing the spread of attractors created from time-series data having chaotic properties. For example, the maximum Lyapunov exponent calculated by the Sano / Sawada algorithm described in Non-Patent Document 3, the macro brain activity index value calculated by the SiCECA algorithm defined in Non-Patent Document 2, and the like.

上記データ処理部4は、上記ピッチ周波数帯域通過フィルター3を介して入力されたピッチ周波数帯域の音声データから上記最大リアプノフ指数やSiCECAアルゴリズムで定義される上記指標値CEMを算出する機能を備えている。
この第1実施形態では、上記データ処理部4が、非特許文献2に記載された上記SiCECAアルゴリズムによってマクロ脳活性度指数を算出する計算手順に基づいてピッチ周波数帯域の音声データのカオス論的指標値CEMを算出するようにしている。
そして、データ処理部4が算出した上記指標値CEMを出力部5へ出力する。
The data processing unit 4 has a function of calculating the maximum Lyapunov exponent and the index value CEM defined by the SiCECA algorithm from the voice data of the pitch frequency band input via the pitch frequency band pass filter 3. .
In the first embodiment, the data processing unit 4 performs a chaotic index of speech data in the pitch frequency band based on a calculation procedure for calculating a macro brain activity index using the SiCECA algorithm described in Non-Patent Document 2. The value CEM is calculated.
Then, the index value CEM calculated by the data processing unit 4 is output to the output unit 5.

上記のようにして被験者のピッチ周波数帯域の音声データから算出したカオス論的指標値である指標値CEMが、被験者の自律神経系の状態を表しているかどうかを確認するため、実験1を行った。
〔実験1〕
50人の被験者それぞれについて、図1に示すシステムを用い、座った状態で発声した音声データと、立った状態で発声した音声データを取得し、これら音声データのピッチ周波数帯域のデータに基づくCEM値を算出した。さらに、座った状態の音声データに基づくCEM値(座位CEM)と、立った状態の音声データに基づくCEM値(立位CEM)を算出し、CEMの変化率として座位CEMに対する立位CEMを演算した。
Experiment 1 was performed in order to confirm whether the index value CEM, which is a chaotic index value calculated from the voice data in the pitch frequency band of the subject as described above, represents the state of the subject's autonomic nervous system. .
[Experiment 1]
For each of the 50 subjects, the system shown in FIG. 1 is used to acquire voice data uttered while sitting and voice data uttered while standing, and CEM values based on the pitch frequency band data of these voice data Was calculated. Furthermore, the CEM value (sitting CEM) based on the voice data in the sitting state and the CEM value (standing CEM) based on the voice data in the standing state are calculated, and the standing CEM with respect to the sitting CEM is calculated as the change rate of the CEM. did.

この実験1の結果を図2のグラフに示す。
図2は、CEM値の変動率に対する被験者の人数を表したグラフであり、ほとんどの被験者の立位CEM/座位CEMは1未満となっている。
このことから、この実施形態のデータ処理部5で算出するピッチ周波数帯域の音声データの指標値CEMは、座っているときよりも立った時に、小さくなることが分かった。
The result of Experiment 1 is shown in the graph of FIG.
FIG. 2 is a graph showing the number of subjects with respect to the variation rate of the CEM value. Most subjects' standing CEM / sitting CEM is less than 1.
From this, it was found that the index value CEM of the sound data in the pitch frequency band calculated by the data processing unit 5 of this embodiment is smaller when standing than when sitting.

一方で、一般的に人は座っているときよりも立っているときの方が、約10%心拍数が高くなることが確認されている。つまり、立位のときの方が座位のときよりも、交感神経が優位となるということである。
これらのことから、自律神経系の状態が変化すると、上記データ処理部4で算出した指標値CEM変化し、しかも交感神経が優位となると指標値CEMが小さくなり、副交感神経が優位になると指標値CEMが大きくなることが分かる。
つまり、この実験1の結果から、自律神経の変化と上記指標値CEMとには相関があることが確認できた。
On the other hand, it has been confirmed that a person generally has a higher heart rate of about 10% when standing than when sitting. In other words, the sympathetic nerve is superior in the standing position than in the sitting position.
Therefore, when the state of the autonomic nervous system changes, the index value CEM calculated by the data processing unit 4 changes, and when the sympathetic nerve becomes dominant, the index value CEM decreases, and when the parasympathetic nerve becomes dominant, the index value It can be seen that the CEM increases.
That is, from the result of Experiment 1, it was confirmed that there was a correlation between the change in the autonomic nerve and the index value CEM.

〔実験2〕
テストコースにおける走行実験を行い、走行中の音声データについて上記指標値CEMを求め、この指標値CEMに対応する時点の脈拍変動率と対比した。
具体的な実験方法は次のとおりである。
被験者は、走行テストを始める前に2分間の安静状態を保ち、その間の安静時脈拍数Aを計測する。
その後、テストコースを決められた速度で所定の距離走行し、走行中に走行時脈拍数Bを計測する。
運転終了後、被験者は安静にして2分間の安静時脈拍数Cを計測する。
[Experiment 2]
A running experiment in a test course was performed, the index value CEM was obtained for the running voice data, and compared with the pulse fluctuation rate at the time corresponding to the index value CEM.
The specific experimental method is as follows.
The subject keeps a resting state for 2 minutes before starting the running test, and measures the resting pulse rate A during that time.
Thereafter, the test course is traveled at a predetermined speed for a predetermined distance, and the running pulse rate B is measured during travel.
After the driving is completed, the subject is rested and the resting pulse rate C for 2 minutes is measured.

なお、車両の走行中には、被験者である運転士の音声データを、図1に示すマイクロフォン1から収集し、この第1実施形態の評価システムによって、上記走行中に収集した音声データに基づいて上記指標値CEMを算出する。この指標値CEMはピッチ周波数帯域の音声データに基づくものである。
さらに、上記安静時脈拍数の平均値(A+C)/2と、走行時脈拍数Bの比率を、脈拍変動率として算出する。
以上の手順を、走行速度、60〔km/H〕、80〔km/H〕、100〔km/H〕、120〔km/H〕の4条件で実行する。
その結果を図3のグラフに示す。
While the vehicle is running, voice data of the driver, who is the subject, is collected from the microphone 1 shown in FIG. 1 and based on the voice data collected during the running by the evaluation system of the first embodiment. The index value CEM is calculated. This index value CEM is based on audio data in the pitch frequency band.
Furthermore, the ratio of the average value (A + C) / 2 of the resting pulse rate and the running pulse rate B is calculated as the pulse fluctuation rate.
The above procedure is executed under four conditions of traveling speed, 60 [km / H], 80 [km / H], 100 [km / H], and 120 [km / H].
The result is shown in the graph of FIG.

図3は、走行速度ごとの走行中の上記指標値CEMを示した太線のグラフvと、走行速度ごとの上記脈拍変動率を示した細線のグラフpである。
これらのグラフから、上記指標値CEMと上記脈拍変動率とは相関があるように見える。実際に、これらの相関係数を算出すと、−0.95となり、上記指標値CEMと脈拍変動率は高い負の相関性を示すことが分かった。
また、脈拍が自律神経の状態に依存することは知られている。
このように、第1実施形態の評価システムから出力される上記指標値CEMが、脈拍変動率と高い相関性を示すということから、上記指標値CEMが自律神経の状態を示すものであることが確認できた。
従って、この実施形態のシステムによれば、被験者の身体に計測装置などを取り付けることなく、音声データから自律神経系の指標値CEMを出力でき、自律神経の状態を評価することができる。
FIG. 3 is a thick line graph v indicating the index value CEM during traveling for each traveling speed, and a thin line graph p indicating the pulse fluctuation rate for each traveling speed.
From these graphs, it appears to be a correlation between the index value CEM and the pulse variation rate. Indeed, if you calculate these correlation coefficients, becomes -0.95, the index value CEM and pulse variation rate was found to exhibit a high negative correlation.
It is also known that the pulse depends on the state of the autonomic nerve.
Thus, the index value CEM outputted from the evaluation system of the first embodiment, since they exhibit pulse variation rate and high correlation, that the index value CEM is indicative of the state of the autonomic nervous It could be confirmed.
Therefore, according to the system of this embodiment, the autonomic nervous system index value CEM can be output from the voice data without attaching a measuring device or the like to the body of the subject, and the state of the autonomic nerve can be evaluated.

なお、図3において点線で示したグラフsは上記実験2と同様にして収集した音声データを、上記ピッチ周波数帯域通過フィルター3を通さないでそのままデータ処理部4に入力して得た指標値を示したものである。この指標値も、上記指標値CEMと同じ計算手順によって算出したものである。上記ピッチ周波数帯域通過フィルター3を通していないため、上記自律神経系の指標値CEMとは異なるが、図3には全帯域のCEMと表示している。
そして、この全帯域の指標値と脈拍変動率との相関係数は+0.55となり、これらの相関性は低いことが分かった。
つまり、被験者の全帯域の音声データではなく、ピッチ周波数帯域の音声データ、すなわち声帯により生成されるデータが、自律神経の状態を表すということである。
Note that a graph s indicated by a dotted line in FIG. 3 shows an index value obtained by inputting the voice data collected in the same manner as in the experiment 2 to the data processing unit 4 without passing through the pitch frequency bandpass filter 3. It is shown. This index value is also calculated by the same calculation procedure as the index value CEM. Since the pitch frequency band pass filter 3 is not passed through, the index value CEM of the autonomic nervous system is different, but in FIG.
The correlation coefficient between the index values of the entire band and the pulse fluctuation rate was +0.55, and it was found that these correlations were low.
That is, the voice data in the pitch frequency band, that is, the data generated by the vocal cords, represents the state of the autonomic nerve, not the voice data in the entire band of the subject.

なお、上記相関係数は、次の式(1)によって算出した値である。

Figure 0005812265
The correlation coefficient is a value calculated by the following equation (1).
Figure 0005812265

上記第1実施形態では、被験者の音声データをもとにピッチ周波数帯域通過フィルター3を設定し、しかも、ピッチ周波数帯域フィルター3を通過させる音声データと、そのフィルターに設定したピッチ周波数帯域を算出するためのケプストラム分析の対象とが一致するようにしている。そのため、一人の被験者のピッチ周波数が測定中に変動したとしても、上記自律神経系の指標値CEMを算出するためのピッチ周波数帯域の音声データを正確に取り出すことができる。
ただし、被験者の音声データに基づいて算出したピッチ周波数帯域を、一旦上記フィルター3に設定したら、その設定をそのまま利用するようにしてもよい。この場合には、個人のピッチ周波数の微妙な変動には対応できないが、一般的なピッチ周波数から外れる声帯をもった被験者のピッチ周波数帯域の音声データを正確に取り出して上記指標値CEMを正確に算出することができる。
In the first embodiment, the pitch frequency band pass filter 3 is set based on the voice data of the subject, and the voice data that passes the pitch frequency band filter 3 and the pitch frequency band set for the filter are calculated. Therefore, the target of cepstrum analysis is matched. Therefore, even if the pitch frequency of one subject fluctuates during measurement, the voice data in the pitch frequency band for calculating the autonomic nervous system index value CEM can be accurately extracted.
However, once the pitch frequency band calculated based on the voice data of the subject is set in the filter 3, the setting may be used as it is. In this case, although it cannot cope with subtle fluctuations in the pitch frequency of an individual, the index value CEM is accurately obtained by accurately extracting the voice data of the subject's pitch frequency band having a vocal cord deviating from the general pitch frequency. Can be calculated.

また、上記第1実施形態では、ピッチ周波数算出部7がケプストラム分析によってピッチ周波数を算出するようにしているが、ピッチ周波数の算出方法はこれに限らない。
例えば、上記SiCECAアルゴリズムにおける埋め込み点が埋め込み空間において相互に接近する周期の逆数を、ピッチ周波数として算出することもできる。また、その他の計算方法を用いてもよい。
ただし、ピッチ周波数帯域として、ほとんどの成人男性あるいは女性が含まれる範囲を設定することは可能である。例えば、平均的な成人男性は、100[Hz]〜150[Hz]、成人女性は200[Hz]〜300[Hz]、子供はその中間である。
そこで、被験者が成人男性、成人女性、あるいは子供であるかによって、上記周波数範囲より多少広い周波数範囲を、予め、上記ピッチ周波数帯域通過フィルター3に固定的に設定しておくこともできる。
In the first embodiment, the pitch frequency calculation unit 7 calculates the pitch frequency by cepstrum analysis, but the pitch frequency calculation method is not limited to this.
For example, the reciprocal of the period in which the embedding points in the SiCECA algorithm approach each other in the embedding space can be calculated as the pitch frequency. Further, other calculation methods may be used.
However, it is possible to set a range including most adult men or women as the pitch frequency band. For example, the average adult male is 100 [Hz] to 150 [Hz], the adult female is 200 [Hz] to 300 [Hz], and the child is intermediate.
Therefore, depending on whether the subject is an adult male, an adult female, or a child, a frequency range slightly wider than the above frequency range can be fixedly set in the pitch frequency bandpass filter 3 in advance.

このように、ピッチ周波数帯域通過フィルター3に、予めピッチ周波数帯域を設定した第2実施形態のシステム構成図を図4に示す。
この第2実施形態のシステムは、データ記憶部6及ピッチ周波数算出部7を備えず、A/D変換器2の出力が直接ピッチ周波数帯通過フィルター3に入力される点が図1に示す第1実施形態と異なる。但し、その他の構成は第1実施形態と同じである。
第1実施形態と同じ構成要素には同じ符号を用い、同じ符号を用いた同じ構成要素は第1実施形態と同じ機能を有する。
FIG. 4 shows a system configuration diagram of the second embodiment in which the pitch frequency band is set in advance in the pitch frequency band pass filter 3 as described above.
The system according to the second embodiment does not include the data storage unit 6 and the pitch frequency calculation unit 7, and the output of the A / D converter 2 is directly input to the pitch frequency band pass filter 3 as shown in FIG. Different from one embodiment. However, other configurations are the same as those in the first embodiment.
The same reference numerals are used for the same constituent elements as those in the first embodiment, and the same constituent elements using the same reference numerals have the same functions as those in the first embodiment.

すなわち、この第2実施形態のシステムも、マイクロフォン1から入力された音声データをA/D変換器2でデジタルデータにして、ピッチ周波数帯域フィルター3を通して取り出されたピッチ周波数帯域の音声データに基づいて、データ処理部がカオス論的指標値の算出手順によって自律神経系の指標値CEMを出力する。
出力された指標値CEMから、交感神経と副交感神経もいずれが優位な状態なのかという自律神経の状態を評価することができる。そのため、この第2実施形態の評価システムでも、被験者の身体には脈拍計や血圧計などの装置を取り付けることなく、音声データのみから簡単に自律神経の状態を評価できることになる。
なお、この第2実施形態においてピッチ周波数帯域通過フィルター3に設定するピッチ周波数帯域は、一般人のピッチ周波数に基づいて設定してもよいし、特定の被験者の音声データから別のシステムによって算出したピッチ周波数に基づいて設定してもよい。
That is, in the system of the second embodiment, the audio data input from the microphone 1 is converted into digital data by the A / D converter 2 and is based on the audio data in the pitch frequency band extracted through the pitch frequency band filter 3. The data processing unit outputs the index value CEM of the autonomic nervous system according to the calculation procedure of the chaotic index value.
From the output index value CEM, it is possible to evaluate the state of the autonomic nerve, which is the dominant state of both the sympathetic nerve and the parasympathetic nerve. Therefore, even in the evaluation system of the second embodiment, the state of the autonomic nerve can be easily evaluated only from the voice data without attaching a device such as a pulse meter or a blood pressure monitor to the body of the subject.
Note that the pitch frequency band set in the pitch frequency band pass filter 3 in the second embodiment may be set based on the pitch frequency of a general person, or a pitch calculated by another system from voice data of a specific subject. You may set based on a frequency.

また、上記第1、第2実施形態では、上記データ処理部4がカオス論的指標値として、SiCECAアルゴリズムに基づいた自律神経系の指標値CEMを算出するようにしているが、自律神経系の指標値としてのカオス論的指標値は、上記CEMに限らない。
例えば、佐野・澤田のアルゴリズムに基づいた最大リアプノフ指数などでもよい。
要するに、被験者のピッチ周波数帯域の音声データのカオス論的揺らぎを評価できる指標値ならば、いずれも自律神経系の指標値とすることができる。
In the first and second embodiments, the data processing unit 4 calculates the index value CEM of the autonomic nervous system based on the SiCECA algorithm as the chaotic index value. The chaotic theoretical index value as the index value is not limited to the CEM.
For example, the maximum Lyapunov exponent based on the Sano / Sawada algorithm may be used.
In short, any index value can be used as an index value of the autonomic nervous system as long as it is an index value that can evaluate the chaotic fluctuation of the voice data in the pitch frequency band of the subject.

この発明の自律神経の状態評価システムは、運転など作業中の自律神経の状態評価に適している。   The autonomic nerve state evaluation system according to the present invention is suitable for autonomic nerve state evaluation during work such as driving.

1 マイクロフォン
2 A/D変換器
3 ピッチ周波数帯域通過フィルター
4 データ処理部
5 出力部
6 データ記憶部
7 ピッチ周波数算出部
DESCRIPTION OF SYMBOLS 1 Microphone 2 A / D converter 3 Pitch frequency bandpass filter 4 Data processing part 5 Output part 6 Data storage part 7 Pitch frequency calculation part

Claims (2)

音声入力手段と、この音声入力手段から入力した音声データをデジタル音声データに変換するA/D変換手段と、デジタル音声データのピッチ周波数帯域を取り出すフィルターと、このフィルターで取り出したピッチ周波数帯のデジタル音声データを処理し、脈拍変動率と相関性を有するカオス論的指標値を算出し、それを出力する信号処理手段とを備えた自律神経の状態評価システム。 Audio input means, A / D conversion means for converting audio data input from the audio input means into digital audio data, a filter for extracting the pitch frequency band of the digital audio data, and a digital of the pitch frequency band extracted by this filter An autonomic nerve state evaluation system comprising: signal processing means for processing voice data , calculating a chaotic index value having a correlation with a pulse fluctuation rate, and outputting the calculated chaotic index value. 上記デジタル音声データを記憶する音声データ記憶手段と、このデータ記憶手段が記憶しているデジタル音声データからピッチ周波数を算出するとともに、このピッチ周波数に基づいてピッチ周波数帯域を算出するピッチ周波数帯域算出手段とを備え、上記フィルターは、上記音声データ記憶手段が記憶しているデジタル音声データから、上記ピッチ周波数帯域算出手段が算出したピッチ周波数帯域の音声データを取り出す請求項1に記載の自律神経の状態評価システム。   Audio data storage means for storing the digital audio data, and pitch frequency band calculation means for calculating the pitch frequency from the digital audio data stored in the data storage means and calculating the pitch frequency band based on the pitch frequency 2. The state of the autonomic nerve according to claim 1, wherein the filter extracts voice data in the pitch frequency band calculated by the pitch frequency band calculating means from the digital voice data stored in the voice data storage means. Evaluation system.
JP2011158788A 2011-07-20 2011-07-20 Autonomic nerve state evaluation system Active JP5812265B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2011158788A JP5812265B2 (en) 2011-07-20 2011-07-20 Autonomic nerve state evaluation system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2011158788A JP5812265B2 (en) 2011-07-20 2011-07-20 Autonomic nerve state evaluation system

Publications (2)

Publication Number Publication Date
JP2013022180A JP2013022180A (en) 2013-02-04
JP5812265B2 true JP5812265B2 (en) 2015-11-11

Family

ID=47781200

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011158788A Active JP5812265B2 (en) 2011-07-20 2011-07-20 Autonomic nerve state evaluation system

Country Status (1)

Country Link
JP (1) JP5812265B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6193826B2 (en) * 2014-08-29 2017-09-06 日本電信電話株式会社 Pulsation sensation information generation apparatus, pulsation sensation information generation method, distribution system, and program
JP7403129B2 (en) * 2018-05-23 2023-12-22 パナソニックIpマネジメント株式会社 Eating and swallowing function evaluation method, program, eating and swallowing function evaluation device, and eating and swallowing function evaluation system

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3764663B2 (en) * 2001-09-14 2006-04-12 独立行政法人電子航法研究所 Psychosomatic diagnosis system
JP3512398B2 (en) * 2001-09-25 2004-03-29 独立行政法人電子航法研究所 Audio processing device
JP2003280684A (en) * 2002-03-25 2003-10-02 Electronic Navigation Research Institute Chaologic diagnosis sensitizing apparatus
JP2010054568A (en) * 2008-08-26 2010-03-11 Oki Electric Ind Co Ltd Emotional identification device, method and program

Also Published As

Publication number Publication date
JP2013022180A (en) 2013-02-04

Similar Documents

Publication Publication Date Title
KR101656611B1 (en) Method for obtaining oxygen desaturation index using unconstrained measurement of bio-signals
KR100954817B1 (en) System and method for testing blood vessel's health and stress through signal analysis of pluse wave
JP4469746B2 (en) Heart rate measuring device and method of operating heart rate measuring device
JP5475658B2 (en) Pain determination device, pain determination program, and control method for pain determination device
US20030055348A1 (en) Apparatus for detecting sleep apnea using electrocardiogram signals
JP6579890B2 (en) Fatigue meter
JP5132003B2 (en) Blood pressure estimation device and blood pressure estimation method
WO2014152565A1 (en) Wavelet analysis in neuro diagnostics
JP6513005B2 (en) Fatigue meter
JP6602248B2 (en) Electronic blood pressure monitor
JPWO2009150744A1 (en) Sleep state monitoring device, monitoring system, and computer program
JP6842214B2 (en) Emotion estimator
CN114974566B (en) Cognitive function assessment method and system
WO2023087512A1 (en) Depression objective screening system, method, and device
JP5812265B2 (en) Autonomic nerve state evaluation system
JP2017064364A (en) Biological state estimation apparatus, biological state estimation method, and computer program
JP4238321B2 (en) Mental stress assessment device
EP3375368A1 (en) Respiration estimating method and device
EP2666408A1 (en) Physiological adaptability system with multiple sensors
WO2019142720A1 (en) Buckle, state of alertness determination system, and state of alertness determination method
JP6557489B2 (en) Biological state estimation device and computer program
CN107374616B (en) Autonomic nerve activity judgment method and device and electrocardiogram monitoring device
KR102053329B1 (en) Method for analyzing stress using biological and exercising signals
JP6742620B2 (en) Swallowing diagnostic device and program
KR101487106B1 (en) Apparatus and method for determining pulse conditions

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20140714

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20140714

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150313

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150407

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150520

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20150908

R150 Certificate of patent or registration of utility model

Ref document number: 5812265

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313115

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250