JP4611206B2 - 疲労度測定装置、疲労検出装置及びコンピュータプログラム - Google Patents
疲労度測定装置、疲労検出装置及びコンピュータプログラム Download PDFInfo
- Publication number
- JP4611206B2 JP4611206B2 JP2005515039A JP2005515039A JP4611206B2 JP 4611206 B2 JP4611206 B2 JP 4611206B2 JP 2005515039 A JP2005515039 A JP 2005515039A JP 2005515039 A JP2005515039 A JP 2005515039A JP 4611206 B2 JP4611206 B2 JP 4611206B2
- Authority
- JP
- Japan
- Prior art keywords
- peak value
- fatigue
- value
- slope
- lyapunov exponent
- 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
Links
- 238000004590 computer program Methods 0.000 title claims description 14
- 238000004364 calculation method Methods 0.000 claims description 94
- 238000004458 analytical method Methods 0.000 claims description 50
- 238000001514 detection method Methods 0.000 claims description 50
- 230000008859 change Effects 0.000 claims description 24
- 238000000034 method Methods 0.000 claims description 17
- 238000012545 processing Methods 0.000 claims description 9
- 238000009499 grossing Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims description 5
- 238000005259 measurement Methods 0.000 abstract description 4
- 238000011002 quantification Methods 0.000 abstract 1
- 206010016256 fatigue Diseases 0.000 description 195
- 230000036544 posture Effects 0.000 description 47
- 230000007958 sleep Effects 0.000 description 47
- 238000002474 experimental method Methods 0.000 description 35
- 238000010586 diagram Methods 0.000 description 26
- 239000006185 dispersion Substances 0.000 description 22
- 238000012360 testing method Methods 0.000 description 17
- 238000011156 evaluation Methods 0.000 description 14
- 230000002093 peripheral effect Effects 0.000 description 13
- 230000007423 decrease Effects 0.000 description 12
- 230000001953 sensory effect Effects 0.000 description 10
- 238000012423 maintenance Methods 0.000 description 9
- 210000003205 muscle Anatomy 0.000 description 9
- 238000001228 spectrum Methods 0.000 description 8
- 208000008035 Back Pain Diseases 0.000 description 7
- 230000003387 muscular Effects 0.000 description 7
- 238000005070 sampling Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 230000003340 mental effect Effects 0.000 description 6
- 208000008930 Low Back Pain Diseases 0.000 description 5
- 230000007704 transition Effects 0.000 description 5
- 206010049565 Muscle fatigue Diseases 0.000 description 4
- 208000024891 symptom Diseases 0.000 description 4
- 206010041349 Somnolence Diseases 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 235000019577 caloric intake Nutrition 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 230000033001 locomotion Effects 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 230000001737 promoting effect Effects 0.000 description 3
- 230000002829 reductive effect Effects 0.000 description 3
- 230000004620 sleep latency Effects 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 241000723438 Cercidiphyllum japonicum Species 0.000 description 2
- JOYRKODLDBILNP-UHFFFAOYSA-N Ethyl urethane Chemical compound CCOC(N)=O JOYRKODLDBILNP-UHFFFAOYSA-N 0.000 description 2
- 208000002193 Pain Diseases 0.000 description 2
- 230000003542 behavioural effect Effects 0.000 description 2
- 230000037149 energy metabolism Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000011514 reflex Effects 0.000 description 2
- 230000033764 rhythmic process Effects 0.000 description 2
- 230000035899 viability Effects 0.000 description 2
- 239000000055 Corticotropin-Releasing Hormone Substances 0.000 description 1
- 102000012289 Corticotropin-Releasing Hormone Human genes 0.000 description 1
- 108010022152 Corticotropin-Releasing Hormone Proteins 0.000 description 1
- 208000019914 Mental Fatigue Diseases 0.000 description 1
- 208000032140 Sleepiness Diseases 0.000 description 1
- 210000003489 abdominal muscle Anatomy 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000037007 arousal Effects 0.000 description 1
- 230000000386 athletic effect Effects 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000009910 autonomic response Effects 0.000 description 1
- 230000031018 biological processes and functions Effects 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 230000035080 detection of muscle activity involved in regulation of muscle adaptation Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 210000003041 ligament Anatomy 0.000 description 1
- 210000004705 lumbosacral region Anatomy 0.000 description 1
- 231100000862 numbness Toxicity 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000003836 peripheral circulation Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 230000037321 sleepiness Effects 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 208000011580 syndromic disease Diseases 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/16—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
- A61B5/18—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state for vehicle drivers or machine operators
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/4538—Evaluating a particular part of the muscoloskeletal system or a particular medical condition
- A61B5/4561—Evaluating static posture, e.g. undesirable back curvature
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6887—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient mounted on external non-worn devices, e.g. non-medical devices
-
- A—HUMAN NECESSITIES
- A63—SPORTS; GAMES; AMUSEMENTS
- A63B—APPARATUS FOR PHYSICAL TRAINING, GYMNASTICS, SWIMMING, CLIMBING, OR FENCING; BALL GAMES; TRAINING EQUIPMENT
- A63B71/00—Games or sports accessories not covered in groups A63B1/00 - A63B69/00
- A63B71/06—Indicating or scoring devices for games or players, or for other sports activities
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/021—Measuring pressure in heart or blood vessels
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/4519—Muscles
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7239—Details of waveform analysis using differentiation including higher order derivatives
-
- A—HUMAN NECESSITIES
- A63—SPORTS; GAMES; AMUSEMENTS
- A63B—APPARATUS FOR PHYSICAL TRAINING, GYMNASTICS, SWIMMING, CLIMBING, OR FENCING; BALL GAMES; TRAINING EQUIPMENT
- A63B2230/00—Measuring physiological parameters of the user
- A63B2230/08—Measuring physiological parameters of the user other bio-electrical signals
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Animal Behavior & Ethology (AREA)
- Physical Education & Sports Medicine (AREA)
- Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Psychiatry (AREA)
- Educational Technology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Dentistry (AREA)
- Developmental Disabilities (AREA)
- Rheumatology (AREA)
- Hospice & Palliative Care (AREA)
- Child & Adolescent Psychology (AREA)
- Psychology (AREA)
- Social Psychology (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
- Two-Way Televisions, Distribution Of Moving Picture Or The Like (AREA)
- Control Of Vending Devices And Auxiliary Devices For Vending Devices (AREA)
Description
本発明は、人の疲労度を定量的に測定するための疲労度検出装置、疲労状態を検出する疲労検出装置及びそれらに用いるコンピュータプログラムに関する。
【背景技術】
人の生体の状態、例えば、活性状態(覚醒状態)であるか、睡眠状態であるかを検出するには、従来、脳波を測定し、その脳波パターンを解析することにより行われている。しかしながら、脳波を測定するには、被検者の頭部に脳波電極や眼電位電極を取り付ける必要があるなど、人の通常動作を制約する環境下で行わなければならず、例えば、自動車、電車などの各種輸送機器の運転時における生体状態を運転者に負担をかけずに評価することは困難である。
一方、運転中の運転者の生体状態(心身状態)を監視することは、近年、事故予防策として注目されており、例えば、特許文献1(特開平9−308614号公報)、特許文献2(特開平10−146321号公報)には、心拍又は脈拍を用いて生体状態を監視する技術が提案されている。特許文献1及び2に開示の技術によれば、脳波測定用の大がかりな装置の頭部への装着が不要で、簡易に運転者の生休状態を評価できる。
特許文献1及び2に開示された装置は、いずれも、心拍又は脈拍についてのカオス指標を算出し、そのカオス指標から運転者の心身状態を判定するものである。具体的には、カオス指標の一つとして、心拍又は脈拍のリアプノフ指数を求め、その時間的変化において、一定時間以上のリアプノフ指数の減少を示した場合に、休憩が必要な程度にストレス負荷が生じ、居眠りに至る直前状態(眠気を感じた状態)となっていると判定する構成である。生体信号のカオス指標によって生体の状態を客観的に診断できることは、既に特開平4−208136号公報により報告されており、特許文献1及び2に開示された装置によって居眠りに至る直前状態を検出することも可能であると推定される。
しかしながら、特許文献1及び2は、測定した心拍又は脈拍を、リアプノフ指数などのカオス指標のみによって処理し、リアプノフ指数の減少変化から精神的に安定した状態に至る過程を検出し、眠気を催すような疲労状態が生じているか否かを判定することを開示しているに過ぎず、疲労度を定量的に捉える試みはなされていない。
また、一般に、体力には、生命を維持するために関わる生存性の体力と、生存性の体力を背景とする活動性の体力がある。生存性の体力は、生命、健康を維持する能力に相当し、防衛体力とも呼ばれ、活動性の体力は体を動かす行動体力であり、一般的には運動能力として理解されている。行動体力を支える機能には、エネルギー発生システム、エネルギー供給システム、及びエネルギー制御システムがある。エネルギー発生システムは、筋肉系の働き具合で、筋力、持久力などの筋肉疲労の元となる機能である。エネルギー供給システムは、酸素摂取量や心拍数から求められる呼吸・循環系の機能であり、エネルギー制御システムは、敏捷性、協調性、平衡性及び適応性の機能である。従って、エネルギー供給システムの状態から肉体的負担を判定でき、エネルギー制御システムの状態から精神的負担を判定できる。そして、これらの判定から、結果として、筋肉疲労の元となるエネルギー発生システムの活動状態を把握できる。
ここで、生体信号データのリアプノフ指数を求めることにより、上記エネルギー制御システムの状態、すなわち、精神的負担の状態を把握でき、生体信号の周期のピーク値から求められる抵抗力(本明細書では、「パワー値」と呼ぶ)を測定処理することにより、上記エネルギー供給システムの状態、すなわち、肉体的負担の状態を把握できる。換言すれば、生体に種々のストレッサーが加わった際に生じる汎適応症候群と呼ばれる生体反応により、ACTH放出ホルモンが多彩な生物作用を生じさせるが、抵抗力とは、かかる作用により現れる心拍出量増加などの自律神経反応や攻撃性の増加などの、外部からの侵入・破壊に対して戦う力であり、それにより、エネルギーを消費し、カロリーを放出する。従って、消費カロリーの低下を検知することにより疲労度を判定することが有効である。本発明においては、かかる抵抗力を生体信号の周期のピーク値から求め、パワー値として定義するものである。
一方、疲労には、末梢性疲労、中枢性疲労があり、例えば、自動車シートに長時間着座した状態では、シートの状態や体調などにより個人差はあるものの、通常、前半は末梢性疲労が優勢で、後半は中枢性疲労が優勢となる。また、疲労に至る過程では、消費カロリーが高いまま、精神的にリラックスし、やがて消費カロリーが低下していく場合と、精神的に刺激を受けている中で消費カロリーが低下し、その後精神的にリラックスし、さらに消費カロリーが低下する場合とがある。すなわち、いずれの疲労も、消費カロリーの低下を伴う点で共通しており、この点からも疲労を捉えるにはパワー値を分析することが有効と考えられる。しかるに、従来は、専らリアプノフ指数の分析によって精神的負担の面から、睡眠に至る直前状態の疲労を検出するものであり、パワー値の分析が考慮されていないため、特に、末梢性疲労の検出には不向きであった。
また、特許文献1及び2におけるリアプノフ指数及び心拍数の値は、時系列変化を捉えているものの、15分あるいは30分ごとの値を用いている。従って、運転時のモニタリングに必要な実質的にリアルタイムの状態変化を見ることはできない。
本発明は、上記に鑑みなされたものであり、疲労度の定量化を図ることができると共に、末梢性疲労、中枢性疲労を問わずに疲労信号を検出でき、特に運転者の疲労度の測定及び疲労の検出に適した疲労度測定装置、疲労検出装置及びコンピュータプログラムを提供することを課題とする。
【発明の開示】
上記した課題を解決するため、請求項1記載の本発明では、指尖容積脈波を測定する生体信号測定器により採取された生体信号データである指尖容積脈波の原波形の各周期のピーク値を検出する生体信号ピーク値検出手段と、
前記生体信号ピーク値検出手段により得られる各ピーク値から、所定時間範囲ごとに上限側のピーク値と下限側のピーク値との差を算出し、この差をパワー値として設定するパワー値算出手段と、
前記パワー値の所定時間範囲における時間軸に対する傾きを、前記所定時間に対して所定のラップ率で所定回数スライド計算して求めるパワー値傾き算出手段と、
前記パワー値傾き算出手段によりスライド計算して得られたパワー値の傾きの時系列信号を絶対値処理して、積分値を算出し、得られた積分値を疲労度として求める疲労度算出手段と
を具備することを特徴とする疲労度測定装置を提供する。
請求項2記載の本発明では、前記生体信号ピーク値検出手段が、前記生体信号データを平滑化微分し、波形の変動幅に対して所定の閾値で上限側のピーク値と下限側のピーク値を求める手段であることを特徴とする請求項1記載の疲労度測定装置を提供する。
請求項3記載の本発明では、前記パワー値算出手段が、前記生体信号データの所定時間範囲における、上限側のピーク値の平均値と下限側のピーク値の平均値との差をパワー値として算出する手段であることを特徴とする請求項1記載の疲労度測定装置を提供する。
請求項4記載の本発明では、前記パワー値算出手段が、前記生体信号データの所定時間範囲における、上限側のピーク値の平均値と下限側のピーク値の平均値との差の二乗値をパワー値として算出する手段であることを特徴とする請求項3記載の疲労度測定装置を提供する。
請求項5記載の本発明では、前記パワー値傾き算出手段におけるスライド計算に用いる時間間隔が180秒であり、ラップ率が90%であることを特徴とする請求項1記載の疲労度測定装置を提供する。
請求項6記載の本発明では、さらに、前記生体信号データをカオス解析して最大リアプノフ指数を算出する最大リアプノフ指数算出手段と、
算出された最大リアプノフ指数の時系列変化波形の各周期のピーク値を検出する最大リアプノフ指数ピーク値検出手段と、
前記パワー値の傾きに加え、前記最大リアプノフ指数ピーク値検出手段により得られる最大リアプノフ指数の各ピーク値の所定時間範囲における時間軸に対する傾きを、前記所定時間に対して所定のラップ率で所定回数スライド計算して求める最大リアプノフ指数傾き算出手段と、
前記パワー値傾き算出手段によりスライド計算して得られたパワー値の傾き及び前記最大リアプノフ指数傾き算出手段によりスライド計算して得られた最大リアプノフ指数の傾きが、時系列信号の中で略180度の位相差を安定して示す時点を疲労信号の発生点として判定する比較・判定手段と
を備えることを特徴とする請求項1記載の疲労度測定装置を提供する。
請求項7記載の本発明では、前記最大リアプノフ指数ピーク値検出手段が、最大リアプノフ指数の時系列変化波形を平滑化微分し、波形の変動幅に対して所定の閾値で上限側のピーク値と下限側のピーク値を求める手段であることを特徴とする請求項6記載の疲労度測定装置を提供する。
請求項8記載の本発明では、前記最大リアプノフ指数の傾き算出手段におけるスライド計算に用いる時間間隔が180秒であり、ラップ率が90%であることを特徴とする請求項6記載の疲労度測定装置を提供する。
請求項9記載の本発明では、指尖容積脈波を測定する生体信号測定器により採取された生体信号データである指尖容積脈波の原波形の各周期のピーク値を検出する生体信号ピーク値検出手段と、
前記生体信号ピーク値検出手段により得られる各ピーク値から、所定時間範囲ごとに上限側のピーク値と下限側のピーク値との差を算出し、この差をパワー値として設定するパワー値算出手段と、
前記パワー値の所定時間範囲における時間軸に対する傾きを、前記所定時間に対して所定のラップ率で所定回数スライド計算して求めるパワー値傾き算出手段と、
前記生体信号データをカオス解析して最大リアプノフ指数を算出する最大リアプノフ指数算出手段と、
算出された最大リアプノフ指数の時系列変化波形の各周期のピーク値を検出する最大リアプノフ指数ピーク値検出手段と、
前記最大リアプノフ指数ピーク値検出手段により得られる最大リアプノフ指数の各ピーク値の所定時間範囲における時間軸に対する傾きを、前記所定時間に対して所定のラップ率で所定回数スライド計算して求める最大リアプノフ指数傾き算出手段と、
前記パワー値傾き算出手段によりスライド計算して得られたパワー値の傾き及び前記最大リアプノフ指数傾き算出手段によりスライド計算して得られた最大リアプノフ指数の傾きが、時系列信号の中で略180度の位相差を安定して示す時点を疲労信号の発生点として判定する比較・判定手段と
を備えることを特徴とする疲労検出装置を提供する。
請求項10記載の本発明では、前記生体信号ピーク値検出手段が、前記生体信号データを平滑化微分し、波形の変動幅に対して所定の閾値で上限側のピーク値と下限側のピーク値を求める手段であり、前記最大リアプノフ指数ピーク値検出手段が、最大リアプノフ指数の時系列変化波形を平滑化微分し、波形の変動幅に対して所定の閾値で上限側のピーク値と下限側のピーク値を求める手段であることを特徴とする請求項9記載の疲労検出装置を提供する。
請求項11記載の本発明では、前記パワー値傾き算出手段及び前記最大リアプノフ指数傾き算出手段におけるスライド計算に用いる時間間隔が180秒であり、ラップ率が90%であることを特徴とする請求項9記載の疲労検出装置を提供する。
請求項12記載の本発明では、人の指尖容積脈波を測定する生体信号測定器により採取された生体信号データである指尖容積脈波の原波形を解析して疲労度を測定するプロセスをコンピュータに実行させるためのコンピュータプログラムであって、
前記生体信号データの原波形の各周期のピーク値を検出する生体信号ピーク値検出ステップと、
前記生休信号ピーク値検出ステップにより得られる各ピーク値から、所定時間範囲ごとに上限側のピーク値と下限側のピーク値との差を算出し、この差をパワー値として設定するパワー値算出ステップと、
前記パワー値の所定時間範囲における時間軸に対する傾きを、前記所定時間に対して所定のラップ率で所定回数スライド計算して求めるパワー値傾き算出ステップと、
前記パワー値傾き算出ステップによりスライド計算して得られたパワー値の傾きの時系列信号を絶対値処理して、積分値を算出し、得られた積分値を疲労度として求める疲労度算出ステップと
をコンピュータに実行させるコンピュータプログラムを提供する。
請求項13記載の本発明では、人の指尖容積脈波を測定する生体信号測定器により採取された生体信号データである指尖容積脈波の原波形を解析して疲労を検出するプロセスをコンピュータに実行させるためのコンピュータプログラムであって、
前記生体信号データの原波形の各周期のピーク値を検出する生体信号ピーク値検出ステップと、
前記生体信号ピーク値検出ステップにより得られる各ピーク値から、所定時間範囲ごとに上限側のピーク値と下限側のピーク値との差を算出し、この差をパワー値として設定するパワー値算出ステップと、
前記パワー値の所定時間範囲における時間軸に対する傾きを、前記所定時間に対して所定のラップ率で所定回数スライド計算して求めるパワー値傾き算出ステップと、
前記生体信号データをカオス解析して最大リアプノフ指数を算出する最大リアプノフ指数算出ステップと、
算出された最大リアプノフ指数の時系列変化波形の各周期のピーク値を検出する最大リアプノフ指数ピーク値検出ステップと、
前記最大リアプノフ指数ピーク値検出ステップにより得られる最大リアプノフ指数の各ピーク値の所定時間範囲における時間軸に対する傾きを、前記所定時間に対して所定のラップ率で所定回数スライド計算して求める最大リアプノフ指数傾き算出ステップと、
前記パワー値傾き算出ステップによりスライド計算して得られたパワー値の傾き及び前記最大リアプノフ指数傾き算出ステップによりスライド計算して得られた最大リアプノフ指数の傾きが、時系列信号の中で略180度の位相差を安定して示す時点を疲労信号の発生点として判定する比較・判定ステップと
をコンピュータに実行させるコンピュータプログラムを提供する。
【発明の効果】
本発明の疲労度測定装置及びコンピュータプログラムは、生体信号データの原波形の各周期のピーク値を検出する生体信号ピーク値検出手段と、生体信号ピーク値検出手段により得られる各ピーク値から、所定時間範囲ごとに上限側のピーク値と下限側のピーク値との差を算出し、この差をパワー値として設定するパワー値算出手段と、パワー値の傾きを求めるパワー値傾き算出手段とを備え、パワー値の傾きの時系列信号を絶対値処理して積分値を算出し、この積分値を疲労度として求める構成である。この結果、人の疲労度の定量化を図ることができた。
また、本発明の疲労検出装置及びコンピュータプログラムは、パワー値の傾きを算出する構成に加え、最大リアプノフ指数の時間軸に対する傾きを求める最大リアプノフ指数傾き算出手段を有する構成とすることにより、疲労信号を検出することができる。また、比較・判定手段により、出現した疲労信号に対応する疲労の種類を判定することができる。
【図面の簡単な説明】
図1は、本発明の一の実施の形態にかかる疲労度測定装置の構成を示すブロック図である。
図2は、スライド計算の手法を説明するための図である。
図3(a)〜(e)は、30分間の睡眠実験において、最適な傾き計算を行うため、サンプリング時間を異ならせた場合のパワー値の傾き及び最大リアプノフ指数の傾きを示す図である。図3(f)はその波高係数を示す図である。
図4(a)〜(d)は、30分間の睡眠実験において、最適な傾き計算を行うため、スライドラップ率を時間を異ならせた場合の/パワー値の傾き及び最大リアプノフ指数の傾きを示す図である。図4(e)は波高係数を示す図である。
図5(a)は、図3(a)〜(e)の周波数分析結果を示す図であり、図5(b)は、図4(a)〜(d)の周波数分析結果を示す図である。
図6は、180分の睡眠実験で得られたパワー値の傾き及び最大リアプノフ指数の傾きについての周波数分析結果を示す図であり、(a)はスライドラップ率を90%として、サンプリング時間を変化させた場合、(b)はサンプリング時間を180秒として、スライドラップ率を変化させた場合を示す。
図7は、短時間着座実験により得られたパワー値の傾きと最大リアプノフ指数の傾きの30分間の時系列変化を示す図であり、(a)は、猫背姿勢で着座した際のデータを、(b)は、同じシートに強制姿勢で着座した際のデータを示す図である。
図8(a),(b)は、図7(a),(b)の猫背姿勢、強制姿勢の周波数分析結果を示す図である。
図9は、3時間静的着座実験により得られたパワー値の傾き及び最大リアプノフ指数の傾きを示す図であり、(a)は体圧分散型シートに着座した場合を、(b)は姿勢維持型シートに着座した場合を示す図である。
図10(a)は、図9(a)の体圧分散型シートの周波数分析結果を示す図であり、図10(b)は、図9(b)の姿勢維持型シートの周波数分析結果を示す図である。
図11(a)は、図9(a)の体圧分散型シート及び図9(b)の姿勢維持型シートの各疲労度を、パワー値の傾きの時系列信号を絶対値処理し、積分値として算出したデータを示す図であり、図11(b)は、体圧分散型シートを基準として姿勢維持型シートを評価した疲労曲線を示す図である。
図12は、ランダム励振下での3時間着座実験により得られたパワー値の傾き及び最大リアプノフ指数の傾きを示す図であり、(a)は体圧分散型シートに着座した場合を、(b)は姿勢維持型シートに着座した場合を示す図である。
図13(a)は、図12(a)の体圧分散型シートの周波数分析結果を示す図であり、図13(b)は、図12(b)の姿勢維持型シートの周波数分析結果を示す図である。
図14(a)は、図12(a)の体圧分散型シート及び図12(b)の姿勢維持型シートの各疲労度を、パワー値の傾きの時系列信号を絶対値処理し、積分値として算出したデータを示す図であり、図14(b)は、体圧分散型シートの疲労度(計算値)及び官能評価値を、姿勢維持型シートを基準として評価した疲労曲線である。
図15は、着座(座位)車両走行実験により得られた運転席の被験者のパワー値の傾き及び最大リアプノフ指数の傾きを示す図である。
図16は、図15の周波数分析結果を示す図である。
図17は、着座(座位)車両走行実験により得られた助手席の被験者のパワー値の傾き及び最大リアプノフ指数の傾きを示す図である。
図18は、図17の周波数分析結果を示す図である。
図19は、疲労度の計算値を示した図であり、図19(a)は運転席被験者の疲労曲線を、図19(b)は助手席被験者の疲労曲線を、図19(c)は両者の疲労曲線を重ね合わせたものをそれぞれ示す。
図20は、車両走行実験による入眠予兆信号、疲労信号及び睡眠信号の波高係数を比較した図である。
図21は、腰痛促進短時間着座実験におけるパワー値の傾き、最大リアプノフ指数の傾きの周波数分析、及び疲労曲線を示した図であり、(a)は、筋肉質被験者のデータを、(b)は腰痛被験者のデータを、(c)は痩せ型被験者のデータを示す
【発明を実施するための最良の形態】
以下、図面に示した実施形態に基づき本発明をさらに詳細に説明する。図1は、本発明の一の実施形態に係る疲労度測定装置1のブロック図である。この図に示したように、本実施形態の疲労度測定装置1は、生体信号測定器10から採取された生体信号データを受信し、所定の分析処理を実行する。
生体信号測定器10は、脈波、心拍などの生体信号を採取できるものであればよいが、指尖容積脈波などの末梢循環の状態を観察するものが好ましい。指尖容積脈波の測定器としては、例えば、赤外線発光ダイオードとフォトトランジスタを備え、指に装着して測定するものを用いることができる。また、例えば、自動車や電車の運転席に着座した人の生体信号を検知する場合には、運転席のシートバックやシートクッションに取り付けられる圧力センサを備え、圧力値変化から脈波を検知するもの用いることができる。なお、この場合、着座時において、人に異物感を感じさせない必要があり、圧力センサとしては、例えば、フィルム状の圧電素子を用い、これを、シートバック、シートクッションの表面などに貼着して取り付けることが好ましい。
疲労度測定装置1は、生体信号測定器10により採取された生体信号データを受信する受信手段を備え、プログラムとしての、最大リアプノフ指数算出手段(最大リアプノフ指数算出ステップ)21と、最大リアプノフ指数ピーク値検出手段(最大リアプノフ指数ピーク値検出ステップ)22と、生体信号ピーク値検出手段(生体信号ピーク値検出ステップ)23と、パワー値算出手段(パワー値算出ステップ)24と、傾き算出手段(傾き算出ステップ)25と、比較・判定手段(比較・判定ステップ)26と、疲労度算出手段(疲労度算出ステップ)27とを備えて構成される。
最大リアプノフ指数は、カオス指標の一つであり、カオスの初期値依存性の程度を指数で示した数値で、カオスアトラクタが描く軌道のうち、近接した2本の軌道間の距離が、時間経過に伴って離れていく度合いを示す量である。具体的には、生体信号測定器10により採取した生体信号データを、最大リアプノフ指数算出手段(最大リアプノフ指数算出ステップ)21により、まず、時間遅れ法によって生体信号(例えば、指尖容積脈波)の時系列信号を状態空間に再構成する。脈波の時系列の遅延時間は、50msで、埋め込み次元はFNN(False Near Neighbors)法を用いると、次元3のときFNNはほぼ零になり、次元4で完全に零になったことから、最適な埋め込み次元を4次元とした。ここで、得られた連続的なデータ計算値に対し、30秒のスライディングウインドウを用いて、図2に示すスライド計算を行い、リアプノフ指数を数値化した。リアプノフ指数の中の最大リアプノフ指数の値を1秒ごとにプロットし、最大リアプノフ指数の時系列データを算定した。次いで、傾き算出手段25において、最大リアプノフ指数の時系列データに対し、図2に示すスライド計算を行うが、詳細は後述する。
本実施形態の最大リアプノフ指数ピーク値検出手段(最大リアプノフ指数ピーク値検出ステップ)22は、上記により計算される最大リアプノフ指数の時系列変化波形の各周期のピーク値を検出する。具体的には、上記により計算される最大リアプノフ指数を、SavitzkyとGolayによる平滑化微分法により平滑化し、上限側のピーク値と下限側のピーク値(ボトム値)を検出する。波形の変動幅に対して所定の閾値で、好ましくは波形の変動幅の70%を閾値として平滑化を併用した一次微分波形で、微分値=0を各ピーク値とする。また、平滑化微分法による平滑化により、ノイズの影響を小さくできる。
生体信号ピーク値検出手段(生体信号ピーク値検出ステップ)23は、生体信号測定器10により得られた生体信号データの原波形の各周期のピーク値を検出するものである。具体的には、生体信号データをSavitzkyとGolayにより平滑化微分し、波形の変動幅に対して所定の閾値で、好ましくは、波形の変動幅の70%を閾値として検出を行い上限側のピーク値と下限側のピーク値(ボトム値)を求める。
パワー値算出手段(パワー値算出ステップ)24においては、生体信号ピーク値検出手段23により得られた生体信号データの各ピーク値を、予め設定した所定の時間範囲ごと、例えば、5秒(s)ごとに切り分け、その時間範囲の中で上限側ピーク値と下限側ピーク値の平均値を求め、それらの差をパワー値として求める。但し、変化量を強調するために、本実施形態では、上記の所定時間範囲における上限側ピーク値の平均値と下限側ピーク値の平均値との差を二乗してパワー値としている。
傾き算出手段(傾き算出ステップ)25は、最大リアプノフ指数の傾き算出手段(最大リアプノフ指数の傾き算出ステップ)と、パワー値の傾き算出手段(パワー値の傾き算出ステップ)とを備えている。最大リアプノフ指数の傾き算出手段は、最大リアプノフ指数ピーク値検出手段22により得られた最大リアプノフ指数の各ピーク値の所定時間範囲における時間軸に対する傾きを、パワー値の傾き算出手段は、パワー値算出手段24により得られたパワー値の所定時間範囲における時間軸に対する傾きを、それぞれ前記所定時間に対して所定のラップ率で所定回数スライド計算して求める(図2参照)。スライド計算は、次のようにして行う。
例えば、T秒(s)間における傾きを、スライドラップ率90%で求める場合には、まず、0(s)〜T(s)間における最大リアプノフ指数のピーク値、及びパワー値の時間軸に対する傾きを、最小二乗近似により求める。次いで、
スライド計算(1):T/10(s)〜T+T/10(s)間、
スライド計算(2):2×T/10(s)〜T+2×T/10(s)間、
スライド計算(n):n×T/10(s)〜T+n×T/10(s)間
における各傾きを最小二乗近似により求めていく。
ここで、最大リアプノフ指数とパワー値の時間領域における特徴を大域的に把握するためには、スライド計算を行う際のサンプリング時間間隔(T秒間)は180秒間が最適であり、スライドラップ率は90%が最適である。これは、数名の被験者について、同じ環境下で30分間の睡眠実験を行い、指尖容積脈波を採取して分析した結果から得られたものである。図3〜図5は、その一例を示すものである。
図3(a)〜(e)は、それぞれ、傾き計算のためのサンプリング時間間隔を、60秒、120秒、180秒、240秒、300秒とし、スライドラップ率を90%に統一して示した最大リアプノフ指数の傾き及びパワー値の傾きであり、図5(a)はその周波数分析結果を示す。なお、図中、aは入眠予兆信号の振幅を示し、bは入眠予兆が現れてから睡眠に遷移する遷移状態の信号の振幅を示し、cは睡眠に入ったときの睡眠信号の振幅を示す。
そして、それぞれにおいて、各兆候信号(ここでは、入眠予兆信号a、遷移状態信号b、睡眠信号c)の時系列信号から、傾きの離散的信号の波高係数:Cf=Xp/Xs(但し、Xpは、兆候信号の最大振幅を示し、Xsは兆候信号が生じる前又は後の定常状態の信号の振幅を示す)を求め、その値から傾きの特徴が最も感度良く現れる条件を求めた。その結果が図3(f)であり、この図から、傾き計算するための時間間隔として180秒が最も感度が良いことがわかる。中央値を180秒に設定した理由は、疲労による筋肉活動の指令の発射頻度は筋肉内の末梢性反射機構が多くになうことによる。つまり、疲労による上位中枢の興奮性の減衰と末梢性の抑制性反射機構の関与により筋肉活動の指令が減少するが、血流を正常に戻すと180秒間で中枢性の興奮水準は回復することと関連していると予測される。
一方、スライドラップ率は、サンプリング時間間隔180秒間の場合において、70%から95%までを計算した。70%未満については時系列信号が疎になるため省略した。その結果が図4(a)〜(d)であり、図5(b)はその周波数分析結果を示す。この図から、スライドラップ率が90%と95%のときがノイズが小さかったが、その波高係数を示した図4(e)のグラフを参照すると、スライドラップ率90%の場合が最も感度が高かった。このことから、兆候信号a,b,cを明確に拾える時間間隔180秒とスライドラップ率90%が、適切な情報抽出のために最も好ましい条件である。
なお、上記した結果は、実験時間30分の場合であるが、実験時間を180分とした場合にも、図6(a),(b)に示したように、兆候信号a,b,cは、時間間隔180秒とスライドラップ率90%とした場合に顕著に特徴を抽出できた。
以上の結果から、傾き算出手段25においては、180秒間に対する最大リアプノフ指数の変化率と180秒間に対するパワー値の変化率を最小二乗法により求め、次いで、18秒後を起点として180秒間に対する傾きを最小二乗法により求めていくことが最適である。
ここで、循環系の生体信号の周波数帯は10Hz以下の周波数帯に集中する。呼吸で0.25〜0.33Hz、心拍数で0.83〜1.17Hz、脈波で0.5〜10Hzである。従来の脈波分析は、脈波の波形の型による分析で血管の硬さ、血液粘度などの情報を得て、10Hz以上の周波数帯のノイズに関してはローパスフィルタを設けることで対応してきている。しかし、10Hz以下の周波数帯のノイズの混入による影響を抑えることは難しく、脈波分析の採取現場は限られてきた。一方、自動車や体動が生じる環境の中で採取される指尖容積脈波は、一般的には不規則な振動源によって加振される振動であり、不規則振動源下でのノイズの影響を小さくしないと、運転者の疲労度を検出するに当たって実用的ではない。ところが、上記のような処理を行い、パワー値の傾き及び最大リアプノフ指数の傾きの大域的な傾向変化を把握する構成とすることで、ノイズの影響を小さく抑えている。
すなわち、指尖容積脈波等の生体信号の原波形の最大リアプノフ指数及びパワー値の微係数(傾き)を、複数回のスライド計算を行って求めるという手法により、ノイズの混入を防ぎ、低周波ゆらぎを精度良く捉えることができる。
比較・判定手段(比較・判定ステップ)26は、傾き算出手段25により得られたパワー値の傾きと最大リアプノフ指数のピーク値の傾きとを比較し、疲労信号の出現を判定する。疲労信号か否かは、後述する試験結果から、パワー値の傾き及び最大リアプノフ指数の傾きとを、同じグラフの時間軸上にプロットしていった場合に、その時系列信号の中で、両者が略180度の位相差(逆位相)を安定して示す特徴的な信号群が出現している時点をもって判定する。後述の試験において、このような特徴的な信号群が出現した時点が、被験者の自己申告コメント、観察者のコメント、及びビデオ撮影記録とを総合的に比較した場合に、疲労や眠気を生じていると判断できるタイミングと一致していた。
なお、略180度の位相差を安定して示す範囲には、疲労兆候を示す疲労信号と、入眠兆候を示す入眠予兆信号とが混在しているが、後述の試験結果から、パワー値の傾きにおいて、入眠時(入眠直後以降の睡眠時)の睡眠信号を基準とした場合に、その直前において出現した略180度位相差を伴う上記特徴的な信号が、睡眠信号の2倍以上の大きな振幅でその後の変化が顕著に小さくなっているのに対し、他のタイミングで出現した略180度位相差を伴う特徴的な信号では、その前後信号に対するパワー値の傾きの振幅差が比較的小さかった。換言すれば、入眠予兆信号は、睡眠に抵抗しながら現れてくるため、その前後に存在する時系列信号とは大きく異なるが、末梢性疲労、中枢性疲労によって生じる疲労信号は、うなりに似た大きな周期性の中に存在する。そこで、前者を入眠予兆信号として判定し、後者を疲労信号として判定して区別した。なお、入眠予兆信号は、疲労の終着点である入眠潜時状態を示す信号でもあり、疲労を示す信号の一種でもある。従って、本実施形態によれば、広い意味での疲労に関する信号が、疲労信号であるか、入眠予兆信号であるかという種類を判定することができる。
比較・判定手段(比較・判定ステップ)26は、疲労状態判定手段(ステップ)を備えている。これは、時系列で出現する上記パワー値の傾きの変化と最大リアプノフ指数の傾きとを周波数分析するものであるが、後述の試験結果から、最大リアプノフ指数の傾きのパワースペクトラムが大きい場合に、中枢性疲労優位状態と判定し、パワー値の傾きのパワースペクトラムが大きい場合に、末梢性疲労優位状態と判定する。これにより、抽出した疲労信号を生じた原因が、末梢性疲労によるものか、中枢性疲労によるのものであるか判定することができる。
疲労度算出手段(疲労度算出ステップ)27は、上記傾き算出手段25におけるパワー値の傾き算出手段により得られたパワー値の傾きの時系列信号を絶対値処理して、積分値を算出することにより、エネルギー代謝量を推定し、該積分値を疲労度(疲労の進行度合い)として算出する構成である。エネルギー代謝量と疲労度とは連動するからである。これにより、疲労度を客観的に把握でき、上記した比較・判定手段26により疲労信号や入眠予兆信号を検出した時点とを総合評価し、例えば、運転者に何らかの注意を喚起すべきか否かを判定することができる。また、疲労度を定量化できるため、その数値により、運転者は、自ら休憩が必要か否かを自覚することができる。
なお、上記した実施形態においては、最大リアプノフ指数算出手段(最大リアプノフ指数算出ステップ)、最大リアプノフ指数ピーク値検出手段(最大リアプノフ指数ピーク値検出ステップ)と、生体信号ピーク値検出手段(生体信号ピーク値検出ステップ)と、パワー値算出手段(パワー値算出ステップ)と、傾き算出手段(傾き算出ステップ)のほかに、比較・判定手段(比較・判定ステップ)及び疲労度算出手段(疲労度算出ステップ)を備えている。従って、上記した実施形態は、本発明の疲労度測定装置と疲労検出装置の両方の機能を備える。疲労度を測定するに当たっては、定量的な疲労度の値のほかに、疲労信号が生じたタイミングを検出する構成とすることが好ましく、それにより、運転者に対する警報装置の発生タイミング等を効果的に制御できる。しかしながら、疲労度算出手段を備えていない構成においても、比較・判定手段を備えることにより、疲労検出装置を形成することも可能である。
また、上記した最大リアプノフ指数算出手段(最大リアプノフ指数算出ステップ)、最大リアプノフ指数ピーク値検出手段(最大リアプノフ指数ピーク値検出ステップ)と、生体信号ピーク値検出手段(生体信号ピーク値検出ステップ)と、パワー値算出手段(パワー値算出ステップ)と、傾き算出手段(傾き算出ステップ)と、比較・判定手段(比較・判定ステップ)、疲労度算出手段(疲労度算出ステップ)等を含んで構成される本発明のコンピュータプログラムは、記録媒体へ記憶させて提供することができる。「記録媒体」とは、それ自身では空間を占有し得ないプログラムを担持することができる媒体であり、例えば、フレキシブルディスク、ハードディスク、CD−ROM、MO(光磁気ディスク)、DVD−ROMなどである。また、本発明に係るプログラムをインストールしたコンピュータから、通信回線を通じて他のコンピュータへ伝送することも可能である。また、汎用的な端末装置に対して、上記のプログラムをプリインストール、あるいはダウンロードすることで、本発明の疲労度測定装置又は疲労検出装置を形成することは、もちろん可能である。
(試験例)
(短時間着座による疲労実験及びその結果)
男性被験者13名により、5〜30分間の短時間着座実験を行った。図7は、そのうちの一人のパワー値の傾きと最大リアプノフ指数の傾きの30分間の時系列変化を示す。図7(a)は、当該被験者が猫背姿勢で着座した際のデータを、図7(b)は、同じシートに当該被験者が強制姿勢で着座した際のデータを示す。なお、猫背姿勢とは、筋力を使用せず、靱帯で姿勢を支持することをいい、強制姿勢とは、胸を張って筋力を使用して姿勢を支持することをいう。一般に猫背姿勢は、筋疲労が少ない状態であるため、中枢性疲労が優位であり、強制姿勢は、末梢性疲労が優位である。
図7の時系列信号を考察すると、図7(a)の中枢性疲労優位の猫背姿勢では、パワー値の傾きの変化率が徐々に減少し、最大リアプノフ指数の傾きもゆらぎながら減少傾向にあり、疲労信号の発生も実験の最初に生じ、後半は実験にも慣れ、精神的・肉体的にリラックスしていると読みとれる。一方、図7(b)の末梢性疲労優位の強制姿勢の場合には、パワー値、最大リアプノフ指数の各傾きは、共に増大傾向にあり、疲労度が非線形で上昇していると考えられ、緊張と姿勢維持のために肉体的・精神的疲労が生じている。従って、比較・判定手段(ステップ)26においては、パワー値及び最大リアプノフ指数の各傾きが前者のような傾向を示した場合に、中枢性疲労が優位であると判定するように設定し、後者のような傾向を示した場合に、末梢性疲労が優位であると判定するように設定することが好ましい。
図8(a),(b)は、図7(a),(b)の猫背姿勢、強制姿勢の周波数分析結果を示す。この図から、中枢性疲労優位状態では、最大リアプノフ指数の傾きのパワースペクトラムが大きく、末梢性疲労は、パワー値の傾きのパワースペクトラムが大きくなる傾向にあることがわかる。従って、比較・判定手段(ステップ)26においては、かかる周波数分析により、いずれの傾きのパワースペクトラムが大きいかを比較することで、中枢性疲労優位か、末梢性疲労優位かを判定する構成とすることがより好ましい。
(長時間着座による疲労実験)
3時間静的着座実験と、米国のミシガン州でワゴン車を用いて採取した1.3Hzにおいて振幅のP−P値で2.0Gの衝撃性振動が発生する突起乗り越えを含むランダム励振下での3時間着座実験を行った。被験者は、いずれも入眠潜時10〜15分の20〜30歳代の男性3名であり、時間帯は午後1時〜4時の間であった。
実験には、腰椎と座骨結節の支持性が高く、最終安定姿勢を維持し易い硬めのクッション性を有する姿勢維持性を重視したバケット型の自動車用シート(姿勢維持型シート)、上記の姿勢維持型を基本としつつも、やわらかめのクッション性を有し、接触面積を広げ体圧のピーク値を低減させた体圧分散型の自動車用シート(体圧分散型シート)を使用した。
(3時間静的着座実験の結果)
図9に、被験者のうちの一人のパワー値の傾きと最大リアプノフ指数の傾きの3時間時系列変化を示す。図9(a)は、体圧分散型シートに着座した場合であり、(b)は、姿勢維持型シートに着座した場合である。
図9(a)におけるT−F1、T−F2、T−F3、T−Fs、図9(b)におけるS−F1、S−F2、S−F3、S−F4、S−Fsが上記した特徴的な信号が現れた信号群を示す。このうち、T−FsとS−Fsは、他の信号と比較し、その絶対値が小さく、その後のパワー値の傾きの変化が小さくなっていることから、入眠予兆信号であると特定できる。そして、入眠予兆信号後の140分で瞬眠が生じていたと判定できる。これは、入眠予兆信号が疲労曲線でいえば危険領域にあり、信号の振幅が小さくなっているため波高係数は急激に大きくなるが、その絶対値は小さくなるものと推察される。一方、T−F1、T−F2、T−F3、S−F1、S−F2、S−F3、S−F4は、振幅は増大傾向にあり、絶対値が大きい。従って、これらの信号は、入眠予兆信号と異なる特性を有しているため、疲労信号として判定する。また、体圧分散型シートに見られるように、疲労信号の振幅が減少し、リズムが安定化していき、定常状態に近づく場合は、上記した図7の結果を参照すると、中枢性疲労が優位になったものと判定できる。
図10(a)は、図9(a)の体圧分散型シートのパワー値の傾きの周波数分析と最大リアプノフ指数の傾きの周波数分析を示す図であり、図10(b)は、図9(b)の姿勢維持型シートのパワー値の傾きの周波数分析と最大リアプノフ指数の傾きの周波数分析を示す図である。
そして、図11(a)は、図9(a)の体圧分散型シート及び図9(b)の姿勢維持型シートの各疲労度を、パワー値の傾きの時系列信号を絶対値処理し、積分値として算出したものである。また、疲労度(計算値)に合わせて、官能評価値をそれぞれグラフ化して示している。なお、官能評価値は、ボルグの指標に基づいて示したものである(「シート感性品質評価」(自動車技術会 学術講演会前刷集、No91−99、21−24、2002)参照)。また、図11(b)は、体圧分散型シートを基準として姿勢維持型シートを定性評価したものである。
図11(a),(b)から、本発明の疲労度測定装置により求めた疲労度(計算値)により示される疲労曲線と官能評価値に基づく疲労曲線とは、定性的傾向がよく似ていることが認められ、本発明の疲労度測定装置により測定された疲労度が、定量的かつ客観的であることを示すことができた。
また、図10(a),(b)の周波数分析及び図11(a),(b)の疲労曲線を考察すると、体圧分散型シートは、姿勢支持性も高く、かつ姿勢変化に対応しやすいシートで、リラックスした状態で着座し続けることができることがわかる。また、最大リアプノフ指数の傾きの周波数分析で疲労信号が増大傾向にあり、中枢性疲労による瞬眠が実験開始後約140分で訪れているが(図9(a)参照)、パワー値の傾きの周波数分析と疲労曲線の各段階がほぼ線形で推移していることから、疲労感の少ないシートであることがわかる。
一方、姿勢維持型シートは、被験者の筋肉が姿勢を維持し、30分までは非常に良好な着座感を示しているが、筋肉疲労が生じる30分過ぎから急激に疲労感が生じ、断続的に疲労信号が生じている。つまり、被験者は実験開始後140分に、このシートでも瞬眠が訪れるが、疲労信号の振幅が増大傾向にあることから、肉体的疲労が原因の瞬眠で、瞬眠と体動の後に疲労度の軽減があることが、120〜150分間の疲労曲線の形から推定できる(図9(b)、図11参照)。また、図10のパワー値の傾きの周波数分析から、体圧分散型シートに比べてエネルギーを消費していることがわかる。最大リアプノフ指数の傾きの周波数分析からは、疲労信号のゆらぎの大きさは増大傾向にあるため、緊張していき精神的刺激の高い状態であると考えられる。
本発明は、このように、疲労度を算出して得られる疲労曲線、及びパワー値の傾き及び最大リアプノフ指数の傾きの周波数分析から、シートの客観的な評価に使用することもできる。
(ランダム励振下での3時間着座実験の結果)
図12は、被験者のうちの一人のパワー値の傾きと最大リアプノフ指数の傾きの3時間時系列変化を示す図である。図12(a)は、体圧分散型シートに着座した場合であり、(b)は、姿勢維持型シートに着座した場合である。
図13(a)は、図12(a)の体圧分散型シートのパワー値の傾きの周波数分析と最大リアプノフ指数の傾きの周波数分析を示す図であり、図13(b)は、図12(b)の姿勢維持型シートのパワー値の傾きの周波数分析と最大リアプノフ指数の傾きの周波数分析を示す図である。
図14(a)は、図12(a)の体圧分散型シート及び図12(b)の姿勢維持型シートの各疲労度を、パワー値の傾きの時系列信号を絶対値処理し、積分値として算出したものである。また、疲労度(計算値)に合わせて、官能評価値をそれぞれグラフ化して示している。図14(b)は、体圧分散型シートの疲労度(計算値)及び官能評価値を、姿勢維持型シートを基準として評価した疲労曲線である。
図14(a)の疲労曲線において、体圧分散型シートは150分まで、姿勢維持型シートは120分まで、計算値で示される疲労度と官能評価値がほぼ一致する。120分の以降の官能評価値との差は、上下振動から誘発されるシートバックからの腰部へのたたきから生じる腰痛の影響によるもので、150分以降の官能評価値との差はしびれによるものと考えられる。この点、図14(b)を見ても、計算値と官能評価値の定性的傾向が150分までは良く一致するが、150分以降は乖離の傾向を示している。また、図12を見ると、パワー値の傾きと最大リアプノフ指数の傾きの変化は、姿勢維持型シートの場合には極端に押さえ込まれている。これは、姿勢維持型シート上の被験者が、痛みに耐えながら着座していたためと予測される。なお、時系列信号の中で、最大リアプノフ指数の傾きが振り切れている信号は、体動が生じていることを示す。
また、図13(a),(b)の周波数分析から、体圧分散型シートは、実験時間内では、疲労度が線形に変化し、適当なゆらぎの中で着座できており、大きな疲労促進は見あたらない。一方、姿勢維持型シートは、実験の早い段階から筋力を使用している。そして、体動により、別の筋肉を使用し、持ち直しを図ろうとするが、姿勢維持のための筋力は収束方向にあり、急速に疲労感が増大していることが読みとれる。
(着座(座位)車両走行実験)
セダン型車両を使用し、運転席と助手席において実施した。運転席の被験者は入眠潜時10〜15分の40歳代の男性であり、助手席の被験者は入眠潜時10〜15分の30歳代の女性である。そして、名神自動車道の桂川PA−山陽自動車道の吉備SA間、吉備SA−小谷SA間において、午後1時半〜午後5時までの3.5時間分の指尖容積脈波を採取した。なお、運転席の被験者、助手席の被験者共に、実験走行開始前に2時間の運転を実施している。
(着座(座位)車両走行実験の結果)
図15は、運転席の被験者のパワー値の傾きと最大リアプノフ指数の傾きの時系列信号を示したものであり、図16は、図15のパワー値の傾きと最大リアプノフ指数の傾きの周波数分析を示したものである。この結果から、運転席の被験者は、桂川PAから吉備SAを走行した100分間は覚醒と疲労が周期的に生じているが、中程度の覚醒状態を継続していることがわかる。さらに、75分から90分の間で肉体的疲労が生じているが、緊張により適応力を上げている。この点は、図16の周波数分析結果の最大リアプノフ指数の傾きのパワースペクトラムが大きくなっていることからも推定できる。
そして、吉備SAから小谷SA間の走行で、運転再開後40分で二回目の疲労が生じている。この疲労は、周波数分析結果に示されるように、パワー値の傾きのパワースペクトラムの増大はなく、図19(a)の矢印aで示す疲労曲線でもわかるように、桂川PAから小谷SAを走行していたときよりも疲労度が下回っている。すなわち、これらの結果から中枢性疲労が優位と推定され、この中枢性疲労には緊張により適応している。なお、この中枢性の疲労も10分後には回復している。
上記の点は、視察、自己申告とも一致した。
また、図16の周波数分析から、桂川PA−吉備SA間は、パワー値の傾きが増大傾向にあり、比較的リラックスしながらも、体力を使った運転であることがわかる。一方、吉備SA−小谷SA間は、中間で疲労が生じているが、緊張感を上げて疲労感に対応していることがわかる。つまり、この運転席の被験者は、体力を使った運転で、疲労が生じると精神力でリズムを作り、疲労に対処していると推定される。
図17は、助手席で着座した被験者のパワー値の傾きと最大リアプノフ指数の傾きの時系列信号を示したものである。この被験者は、生体信号を採取する前の2時間の走行により疲労状態にあった。これは、図19(b)の疲労曲線の矢印bで示されるように、桂川PAを出発して5分後に疲労曲線が急激な立ち上がりを示すことから推定される。そして、90〜100分の間に、入眠予兆信号と瞬眠が生じている。これは、図18の周波数分析結果から、パワー値と最大リアプノフ指数の各傾きのパワースペクトラムが共に高い水準を示していることから、蓄積された肉体疲労に中枢性疲労が加味されて生じたものと推測される。入眠予兆信号が生じて約5分後に瞬眠が生じ一時的に回復し、吉備SAで休憩に入った。再び走行した後、40分後に入眠予兆信号の第2波が訪れている。そして、その5分後に瞬眠が生じ、その後、覚醒状態に移行している。時系列信号のパワー値の傾きの振幅の増減と逆位相状態を保持する時間の割合から、この被験者は総じて疲労状態にあることがわかる。なお、吉備SA−小谷SA間は瞬眠のあと適応力が上昇し、帰着による期待感から精神力が高揚し、疲労感が持ち直している。この点は、図18の周波数分析からもわかる。また、視察、被験者のコメントとも一致していた。
一方、図19は、疲労度の計算値を示す。図19(a)は運転席被験者の疲労曲線であり、図19(b)は助手席被験者の疲労曲線であり、図19(c)は両者の疲労曲線を重ね合わせて示したものである。この疲労曲線から、男性と女性の体力差及び運転しているか否か、あるいは帰着の期待感などの精神的条件差による疲労度の遷移傾向がよく捉えられていることがわかる。図20は、車両走行実験による入眠予兆信号、疲労信号及び睡眠信号の波高係数を比較したものである。波高係数は、外的刺激を受ける環境であるため、比較的高い水準にある。
また、図19の計算による疲労度の推移は、被験者の感覚及びコメントと良く一致しており、本発明による疲労度の計算が、実車走行においても有効であることがわかった。
(腰痛促進短時間着座実験)
シートフレームの撓みを小さくし、座部に厚さ30mmのチップウレタンを敷き、背部の腰椎支持部に脱着可能で厚さ調節可能なウレタンをおいたシートを用い、体格差の異なる3名の男性により実施した。被験者の特徴は、筋肉質、腰痛持ち、痩せ型である。実験においては、背もたれ無しの自由姿勢、立位の脊柱カーブに準じた姿勢(通常運転姿勢)、腰部の張り出しを強調した姿勢の3種類で姿勢を変化させて行った。
(腰痛促進短時間着座実験の結果)
図21は、本実験におけるパワー値の傾き、最大リアプノフ指数の傾きの周波数分析、及び疲労曲線を示したものである。図21(a)は、筋肉質被験者のデータを、(b)は腰痛被験者のデータを、(c)は痩せ型被験者のデータを示す。この図から、腰痛持ちの被験者は、他の2人と異なり、早い時期に疲労度が上昇し、背もたれのない場合には痛みに耐えているのがわかる。また、筋肉質と痩せ型の被験者は、シートの姿勢支持の方法の差に拘わらず、比較的安定した状態を示す。特に、筋肉質の被験者は、その背筋、腹筋で姿勢を維持できているのが分かる。なお、疲労度は体力により、疲労の進行度に差が生じていた。これらの点は、被験者のコメントとよく一致していた。
従って、この実験においても、本発明による疲労の検出及び疲労度の測定が、客観的かつ定量的に示すことができることがわかる。
【産業上の利用可能性】
本発明では、上記のように、疲労度を定量化し、中枢性疲労のみではなく肉体的疲労も合わせて、客観的に算出することが可能である。また、疲労度を算出する過程で示されるパワー値の傾き及び最大リアプノフ指数の傾きの時系列信号の状態、さらには周波数分析の状態から、疲労信号、入眠予兆信号が生じたタイミングやそれらの種別も判定することができる。従って、これらを用いて、所定の疲労度に至った場合、あるいは、疲労信号や入眠予兆信号を検知した場合に、任意の警告装置を機能させるようにすることもできる。また、疲労度が客観的に表示され、疲労信号や入眠予兆信号も検知できるため、人が着座するシート、あるいは寝具などの性能評価に用いることもできる。同時に、その際の人の心身状態を把握する診断装置として用いることもできる。さらに、生体信号データに代え、例えば、機械装置の微振動を検出する測定器を用いることにより、機械装置の故障診断などを行うことに、本発明を応用することも考えられる。
Claims (13)
- 指尖容積脈波を測定する生体信号測定器により採取された生体信号データである指尖容積脈波の原波形の各周期のピーク値を検出する生体信号ピーク値検出手段と、
前記生体信号ピーク値検出手段により得られる各ピーク値から、所定時間範囲ごとに上限側のピーク値と下限側のピーク値との差を算出し、この差をパワー値として設定するパワー値算出手段と、
前記パワー値の所定時間範囲における時間軸に対する傾きを、前記所定時間に対して所定のラップ率で所定回数スライド計算して求めるパワー値傾き算出手段と、
前記パワー値傾き算出手段によりスライド計算して得られたパワー値の傾きの時系列信号を絶対値処理して、積分値を算出し、得られた積分値を疲労度として求める疲労度算出手段と
を具備することを特徴とする疲労度測定装置。 - 前記生体信号ピーク値検出手段が、前記生体信号データを平滑化微分し、波形の変動幅に対して所定の閾値で上限側のピーク値と下限側のピーク値を求める手段であることを特徴とする請求項1記載の疲労度測定装置。
- 前記パワー値算出手段が、前記生体信号データの所定時間範囲における、上限側のピーク値の平均値と下限側のピーク値の平均値との差をパワー値として算出する手段であることを特徴とする請求項1記載の疲労度測定装置。
- 前記パワー値算出手段が、前記生体信号データの所定時間範囲における、上限側のピーク値の平均値と下限側のピーク値の平均値との差の二乗値をパワー値として算出する手段であることを特徴とする請求項3記載の疲労度測定装置。
- 前記パワー値傾き算出手段におけるスライド計算に用いる時間間隔が180秒であり、ラップ率が90%であることを特徴とする請求項1記載の疲労度測定装置。
- さらに、前記生体信号データをカオス解析して最大リアプノフ指数を算出する最大リアプノフ指数算出手段と、
算出された最大リアプノフ指数の時系列変化波形の各周期のピーク値を検出する最大リアプノフ指数ピーク値検出手段と、
前記パワー値の傾きに加え、前記最大リアプノフ指数ピーク値検出手段により得られる最大リアプノフ指数の各ピーク値の所定時間範囲における時間軸に対する傾きを、前記所定時間に対して所定のラップ率で所定回数スライド計算して求める最大リアプノフ指数傾き算出手段と、
前記パワー値傾き算出手段によりスライド計算して得られたパワー値の傾き及び前記最大リアプノフ指数傾き算出手段によりスライド計算して得られた最大リアプノフ指数の傾きが、時系列信号の中で略180度の位相差を安定して示す時点を疲労信号の発生点として判定する比較・判定手段と
を備えることを特徴とする請求項1記載の疲労度測定装置。 - 前記最大リアプノフ指数ピーク値検出手段が、最大リアプノフ指数の時系列変化波形を平滑化微分し、波形の変動幅に対して所定の閾値で上限側のピーク値と下限側のピーク値を求める手段であることを特徴とする請求項6記載の疲労度測定装置。
- 前記最大リアプノフ指数の傾き算出手段におけるスライド計算に用いる時間間隔が180秒であり、ラップ率が90%であることを特徴とする請求項6記載の疲労度測定装置。
- 指尖容積脈波を測定する生体信号測定器により採取された生体信号データである指尖容積脈波の原波形の各周期のピーク値を検出する生体信号ピーク値検出手段と、
前記生体信号ピーク値検出手段により得られる各ピーク値から、所定時間範囲ごとに上限側のピーク値と下限側のピーク値との差を算出し、この差をパワー値として設定するパワー値算出手段と、
前記パワー値の所定時間範囲における時間軸に対する傾きを、前記所定時間に対して所定のラップ率で所定回数スライド計算して求めるパワー値傾き算出手段と、
前記生体信号データをカオス解析して最大リアプノフ指数を算出する最大リアプノフ指数算出手段と、
算出された最大リアプノフ指数の時系列変化波形の各周期のピーク値を検出する最大リアプノフ指数ピーク値検出手段と、
前記最大リアプノフ指数ピーク値検出手段により得られる最大リアプノフ指数の各ピーク値の所定時間範囲における時間軸に対する傾きを、前記所定時間に対して所定のラップ率で所定回数スライド計算して求める最大リアプノフ指数傾き算出手段と、
前記パワー値傾き算出手段によりスライド計算して得られたパワー値の傾き及び前記最大リアプノフ指数傾き算出手段によりスライド計算して得られた最大リアプノフ指数の傾きが、時系列信号の中で略180度の位相差を安定して示す時点を疲労信号の発生点として判定する比較・判定手段と
を備えることを特徴とする疲労検出装置。 - 前記生体信号ピーク値検出手段が、前記生体信号データを平滑化微分し、波形の変動幅に対して所定の閾値で上限側のピーク値と下限側のピーク値を求める手段であり、前記最大リアプノフ指数ピーク値検出手段が、最大リアプノフ指数の時系列変化波形を平滑化微分し、波形の変動幅に対して所定の閾値で上限側のピーク値と下限側のピーク値を求める手段であることを特徴とする請求項9記載の疲労検出装置。
- 前記パワー値傾き算出手段及び前記最大リアプノフ指数傾き算出手段におけるスライド計算に用いる時間間隔が180秒であり、ラップ率が90%であることを特徴とする請求項9記載の疲労検出装置。
- 人の指尖容積脈波を測定する生体信号測定器により採取された生体信号データである指尖容積脈波の原波形を解析して疲労度を測定するプロセスをコンピュータに実行させるためのコンピュータプログラムであって、
前記生体信号データの原波形の各周期のピーク値を検出する生体信号ピーク値検出ステップと、
前記生休信号ピーク値検出ステップにより得られる各ピーク値から、所定時間範囲ごとに上限側のピーク値と下限側のピーク値との差を算出し、この差をパワー値として設定するパワー値算出ステップと、
前記パワー値の所定時間範囲における時間軸に対する傾きを、前記所定時間に対して所定のラップ率で所定回数スライド計算して求めるパワー値傾き算出ステップと、
前記パワー値傾き算出ステップによりスライド計算して得られたパワー値の傾きの時系列信号を絶対値処理して、積分値を算出し、得られた積分値を疲労度として求める疲労度算出ステップと
をコンピュータに実行させるコンピュータプログラム。 - 人の指尖容積脈波を測定する生体信号測定器により採取された生体信号データである指尖容積脈波の原波形を解析して疲労を検出するプロセスをコンピュータに実行させるためのコンピュータプログラムであって、
前記生体信号データの原波形の各周期のピーク値を検出する生体信号ピーク値検出ステップと、
前記生体信号ピーク値検出ステップにより得られる各ピーク値から、所定時間範囲ごとに上限側のピーク値と下限側のピーク値との差を算出し、この差をパワー値として設定するパワー値算出ステップと、
前記パワー値の所定時間範囲における時間軸に対する傾きを、前記所定時間に対して所定のラップ率で所定回数スライド計算して求めるパワー値傾き算出ステップと、
前記生体信号データをカオス解析して最大リアプノフ指数を算出する最大リアプノフ指数算出ステップと、
算出された最大リアプノフ指数の時系列変化波形の各周期のピーク値を検出する最大リアプノフ指数ピーク値検出ステップと、
前記最大リアプノフ指数ピーク値検出ステップにより得られる最大リアプノフ指数の各ピーク値の所定時間範囲における時間軸に対する傾きを、前記所定時間に対して所定のラップ率で所定回数スライド計算して求める最大リアプノフ指数傾き算出ステップと、
前記パワー値傾き算出ステップによりスライド計算して得られたパワー値の傾き及び前記最大リアプノフ指数傾き算出ステップによりスライド計算して得られた最大リアプノフ指数の傾きが、時系列信号の中で略180度の位相差を安定して示す時点を疲労信号の発生点として判定する比較・判定ステップと
をコンピュータに実行させるコンピュータプログラム。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003363902 | 2003-10-23 | ||
JP2003363902 | 2003-10-23 | ||
PCT/JP2004/016058 WO2005039415A1 (ja) | 2003-10-23 | 2004-10-22 | 疲労度測定装置、疲労検出装置及びコンピュータプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2005039415A1 JPWO2005039415A1 (ja) | 2007-02-15 |
JP4611206B2 true JP4611206B2 (ja) | 2011-01-12 |
Family
ID=34510078
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2005515039A Expired - Fee Related JP4611206B2 (ja) | 2003-10-23 | 2004-10-22 | 疲労度測定装置、疲労検出装置及びコンピュータプログラム |
Country Status (7)
Country | Link |
---|---|
US (1) | US7470231B2 (ja) |
EP (1) | EP1685798B1 (ja) |
JP (1) | JP4611206B2 (ja) |
CN (1) | CN100411589C (ja) |
AT (1) | ATE499878T1 (ja) |
DE (1) | DE602004031659D1 (ja) |
WO (1) | WO2005039415A1 (ja) |
Families Citing this family (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4636861B2 (ja) * | 2004-11-25 | 2011-02-23 | 株式会社デルタツーリング | 快適感評価装置 |
JP4832914B2 (ja) * | 2006-02-08 | 2011-12-07 | 株式会社デルタツーリング | 筋疲労評価装置 |
RU2309668C1 (ru) | 2006-02-20 | 2007-11-10 | Александр Сергеевич Парфенов | Способ неинвазивного определения функции эндотелия и устройство для его осуществления |
JP5044230B2 (ja) * | 2007-02-14 | 2012-10-10 | 株式会社デルタツーリング | 生体信号分析装置、シート及び生体信号分析方法 |
JP2008195301A (ja) * | 2007-02-14 | 2008-08-28 | Delta Tooling Co Ltd | シート |
JP5179798B2 (ja) | 2007-05-22 | 2013-04-10 | 株式会社デルタツーリング | 座席構造 |
JP5022530B2 (ja) * | 2007-07-20 | 2012-09-12 | 株式会社デルタツーリング | 疲労度演算装置及びコンピュータプログラム |
JP5086834B2 (ja) * | 2008-02-20 | 2012-11-28 | 株式会社デルタツーリング | 疲労解析装置及びコンピュータプログラム |
US8945028B2 (en) * | 2008-05-20 | 2015-02-03 | Ekso Bionics, Inc. | Device and method for decreasing energy consumption of a person by use of a lower extremity exoskeleton |
CA2724062C (en) * | 2008-05-20 | 2019-11-26 | The Regents Of The University Of California | Device and method for decreasing oxygen consumption of a person during steady walking by use of a load-carrying exoskeleton |
CN102098986B (zh) * | 2008-07-23 | 2015-09-09 | 伯克利仿生技术公司 | 外骨骼和控制该外骨骼的迈步腿的方法 |
EP2329769B1 (en) * | 2008-08-19 | 2015-03-25 | Delta Tooling Co., Ltd. | Biometric signal measuring device and organism condition analyzing system |
EP2422700B1 (en) | 2009-04-25 | 2019-06-05 | Delta Tooling Co., Ltd. | Device and computer program for analyzing the state of a living body |
US8057400B2 (en) | 2009-05-12 | 2011-11-15 | Angiologix, Inc. | System and method of measuring changes in arterial volume of a limb segment |
JP5444898B2 (ja) * | 2009-07-09 | 2014-03-19 | アイシン精機株式会社 | 状態検出装置、状態検出方法及びプログラム |
JP5704651B2 (ja) | 2009-10-14 | 2015-04-22 | 株式会社デルタツーリング | 生体状態推定装置、生体状態推定システム及びコンピュータプログラム |
US9030294B2 (en) | 2010-09-20 | 2015-05-12 | Pulsar Informatics, Inc. | Systems and methods for collecting biometrically verified actigraphy data |
JP5733499B2 (ja) * | 2010-10-29 | 2015-06-10 | 株式会社デルタツーリング | 生体状態推定装置及びコンピュータプログラム |
US20140073936A1 (en) * | 2012-09-11 | 2014-03-13 | Nellcor Puritan Bennett Llc | Methods and systems for conditioning physiological information based on a stability function |
US9751534B2 (en) | 2013-03-15 | 2017-09-05 | Honda Motor Co., Ltd. | System and method for responding to driver state |
NZ630770A (en) | 2013-10-09 | 2016-03-31 | Resmed Sensor Technologies Ltd | Fatigue monitoring and management system |
US9814934B2 (en) * | 2014-06-11 | 2017-11-14 | Brian Alexander Mabrey | Baseline attenuated muscle (BAM) method |
CN104720799A (zh) * | 2015-04-03 | 2015-06-24 | 上海帝仪科技有限公司 | 基于低频脑电信号的疲劳检测方法及系统 |
CN104825174B (zh) * | 2015-04-17 | 2017-11-07 | 深圳还是威健康科技有限公司 | 一种疲劳状态的检测方法及终端 |
CN105030244B (zh) * | 2015-06-29 | 2018-05-11 | 杭州镜之镜科技有限公司 | 一种眨眼的检测方法和检测系统 |
CN105105732A (zh) * | 2015-08-14 | 2015-12-02 | 姚丽峰 | 实现脉搏心率检测以及情绪和疲劳程度监测的系统及方法 |
CN105467807B (zh) * | 2015-11-19 | 2018-10-30 | 西安交通大学 | 面向全息展示的多通道人机交互工效测评系统及测评方法 |
JP6216406B1 (ja) * | 2016-04-27 | 2017-10-18 | 株式会社デンソーアイティーラボラトリ | ドライバ緊張度判定装置及びドライバ緊張度判定方法 |
WO2018014953A1 (en) * | 2016-07-20 | 2018-01-25 | Toyota Motor Europe | Control device, system and method for determining a comfort level of a driver |
JP6676499B2 (ja) * | 2016-08-12 | 2020-04-08 | オムロンヘルスケア株式会社 | 疲労度判定装置、疲労度判定方法、疲労度判定プログラム、及び、生体情報測定装置 |
US9795322B1 (en) | 2016-10-14 | 2017-10-24 | Right Posture Pte. Ltd. | Methods and systems for monitoring posture with alerts and analytics generated by a smart seat cover |
WO2018074371A1 (ja) * | 2016-10-21 | 2018-04-26 | シチズン時計株式会社 | 検出装置 |
CN107007290B (zh) * | 2017-03-27 | 2019-09-20 | 广州视源电子科技股份有限公司 | 基于时域及相空间的脑电放松度识别方法及装置 |
CN108464831A (zh) * | 2018-04-19 | 2018-08-31 | 福州大学 | 一种穿戴式肌肉疲劳检测的装置及方法 |
WO2020043855A1 (en) * | 2018-08-29 | 2020-03-05 | Sleep Advice Technologies S.R.L. | Photoplethysmography based detection of transitions between awake, drowsiness, and sleep phases of a subject |
JP7217576B2 (ja) * | 2018-11-05 | 2023-02-03 | 株式会社安藤・間 | ドライバーの状態推定方法及び装置 |
CN110477897B (zh) * | 2019-08-15 | 2022-04-01 | 中国体育国际经济技术合作有限公司 | 一种体能测试方法及系统 |
US20220161748A1 (en) * | 2020-11-26 | 2022-05-26 | Pixart Imaging Inc. | Multiple detection system |
CN112716469B (zh) * | 2020-12-29 | 2022-07-19 | 厦门大学 | 一种基于指尖视频的实时心率提取方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09308614A (ja) * | 1996-05-21 | 1997-12-02 | Matsushita Electric Ind Co Ltd | 運転者モニター装置とこれを用いた安全装置 |
JP2004344613A (ja) * | 2003-05-21 | 2004-12-09 | Delta Tooling Co Ltd | 運転席用座席システム及び覚醒状態復帰装置 |
JP2004344612A (ja) * | 2003-05-21 | 2004-12-09 | Delta Tooling Co Ltd | 生体評価システム、コンピュータプログラム及び記録媒体 |
Family Cites Families (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CH681990A5 (ja) | 1989-04-14 | 1993-06-30 | Matec Holding | |
JPH0339236A (ja) | 1989-07-07 | 1991-02-20 | Asahi Chem Ind Co Ltd | 視力補正用器具の製造方法 |
JPH0349987A (ja) | 1989-07-18 | 1991-03-04 | Mitsubishi Paper Mills Ltd | 感熱記録材料 |
JPH0370327A (ja) | 1989-08-10 | 1991-03-26 | Nec Corp | 秘話端末装置 |
JPH0370346A (ja) | 1989-08-10 | 1991-03-26 | Nec Corp | データ入力システム |
JPH0717077B2 (ja) | 1989-09-04 | 1995-03-01 | 三菱電機株式会社 | プリンタの用紙給排路構造 |
JPH069546B2 (ja) * | 1990-11-30 | 1994-02-09 | 株式会社コンピュータコンビニエンス | 体表面から採取した脈波及び/又は心拍を用いる診断装置 |
DE4203084A1 (de) | 1992-02-04 | 1993-08-05 | Bayer Ag | Substituierte aminosaeureamide |
JP3047712B2 (ja) | 1993-11-19 | 2000-06-05 | セイコーエプソン株式会社 | 脈波診断装置 |
JP3070346B2 (ja) * | 1993-08-09 | 2000-07-31 | セイコーエプソン株式会社 | ストレス評価装置および生理的年齢評価装置 |
JP3092369B2 (ja) | 1993-01-07 | 2000-09-25 | セイコーエプソン株式会社 | 脈波解析装置 |
JP3049987B2 (ja) | 1993-03-11 | 2000-06-05 | セイコーエプソン株式会社 | 脈波解析装置 |
JP3070327B2 (ja) | 1993-03-08 | 2000-07-31 | セイコーエプソン株式会社 | 診断装置 |
JP3039236B2 (ja) | 1993-11-30 | 2000-05-08 | セイコーエプソン株式会社 | 脈波解析装置 |
WO1994015526A1 (en) | 1993-01-07 | 1994-07-21 | Seiko Epson Corporation | Pulse wave analyzer, and diagnosis apparatus using the same |
JP3011358B2 (ja) * | 1994-11-16 | 2000-02-21 | パイオニア株式会社 | 心拍測定装置 |
JP3811995B2 (ja) * | 1996-07-03 | 2006-08-23 | 松下電器産業株式会社 | 生体モニタリング装置とそれを備えた携帯用保護具 |
JP3596198B2 (ja) * | 1996-11-20 | 2004-12-02 | 松下電器産業株式会社 | 運転者監視装置 |
US6483484B1 (en) * | 1998-12-18 | 2002-11-19 | Semiconductor Energy Laboratory Co., Ltd. | Goggle type display system |
JP4064028B2 (ja) * | 2000-01-05 | 2008-03-19 | 株式会社タニタ | 身体疲労度判定装置 |
US20020156392A1 (en) * | 2001-03-06 | 2002-10-24 | Mitsubishi Chemical Corporation | Method and apparatus for inspecting biological rhythms |
JP2002272708A (ja) * | 2001-03-21 | 2002-09-24 | Mitsuba Corp | 体調判定方法および体調判定装置 |
US6954700B2 (en) * | 2001-04-25 | 2005-10-11 | Computer Convenience Inc. | Efficacy of biological state and action affecting biological state, judging apparatus, judging system, judging program and recording medium holding the program |
JP2003061921A (ja) * | 2001-08-30 | 2003-03-04 | Mitsuba Corp | 体調判別方法および体調判別装置 |
WO2005000119A1 (ja) * | 2003-06-27 | 2005-01-06 | Soiken Inc. | 疲労度評価方法、疲労度評価装置、およびデータベース |
-
2004
- 2004-10-22 WO PCT/JP2004/016058 patent/WO2005039415A1/ja active Application Filing
- 2004-10-22 DE DE602004031659T patent/DE602004031659D1/de active Active
- 2004-10-22 CN CNB2004800312169A patent/CN100411589C/zh not_active Expired - Fee Related
- 2004-10-22 JP JP2005515039A patent/JP4611206B2/ja not_active Expired - Fee Related
- 2004-10-22 AT AT04793166T patent/ATE499878T1/de not_active IP Right Cessation
- 2004-10-22 EP EP04793166A patent/EP1685798B1/en not_active Not-in-force
- 2004-10-22 US US10/576,173 patent/US7470231B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09308614A (ja) * | 1996-05-21 | 1997-12-02 | Matsushita Electric Ind Co Ltd | 運転者モニター装置とこれを用いた安全装置 |
JP2004344613A (ja) * | 2003-05-21 | 2004-12-09 | Delta Tooling Co Ltd | 運転席用座席システム及び覚醒状態復帰装置 |
JP2004344612A (ja) * | 2003-05-21 | 2004-12-09 | Delta Tooling Co Ltd | 生体評価システム、コンピュータプログラム及び記録媒体 |
Also Published As
Publication number | Publication date |
---|---|
US7470231B2 (en) | 2008-12-30 |
DE602004031659D1 (de) | 2011-04-14 |
ATE499878T1 (de) | 2011-03-15 |
US20070078351A1 (en) | 2007-04-05 |
CN1870939A (zh) | 2006-11-29 |
JPWO2005039415A1 (ja) | 2007-02-15 |
EP1685798A1 (en) | 2006-08-02 |
WO2005039415A1 (ja) | 2005-05-06 |
CN100411589C (zh) | 2008-08-20 |
EP1685798A4 (en) | 2009-07-01 |
EP1685798B1 (en) | 2011-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4611206B2 (ja) | 疲労度測定装置、疲労検出装置及びコンピュータプログラム | |
US7250029B2 (en) | Human condition evaluation system, computer program, and computer-readable record medium | |
JP4637098B2 (ja) | 負荷体状態判定装置、乗物用シート及びコンピュータプログラム | |
JP5553303B2 (ja) | 生体状態推定装置及びコンピュータプログラム | |
JP5704651B2 (ja) | 生体状態推定装置、生体状態推定システム及びコンピュータプログラム | |
JP6118097B2 (ja) | 運転時生体状態判定装置及びコンピュータプログラム | |
KR100851383B1 (ko) | 심전도를 이용한 운전자 스트레스지수 제공시스템 및 방법 | |
WO2012057331A1 (ja) | 生体状態推定装置及びコンピュータプログラム | |
JP2004344613A (ja) | 運転席用座席システム及び覚醒状態復帰装置 | |
JP5086834B2 (ja) | 疲労解析装置及びコンピュータプログラム | |
JP5036383B2 (ja) | 睡眠状態判定装置、睡眠状態判定方法及びコンピュータプログラム | |
JP4636861B2 (ja) | 快適感評価装置 | |
JP4832914B2 (ja) | 筋疲労評価装置 | |
JP5327584B2 (ja) | 生体状態分析装置、コンピュータプログラム及び記録媒体 | |
JP4502712B2 (ja) | 心身状態判定システム | |
JP2017064364A (ja) | 生体状態推定装置、生体状態推定方法及びコンピュータプログラム | |
JPWO2006057313A1 (ja) | 生体信号検出装置 | |
JP6209395B2 (ja) | 生体状態推定装置及びコンピュータプログラム | |
JP6588035B2 (ja) | 生体状態分析装置及びコンピュータプログラム | |
WO2017057022A1 (ja) | 生体状態推定装置、生体状態推定方法及びコンピュータプログラム | |
JP4636860B2 (ja) | 人体支持構造の評価装置 | |
JP2016112144A (ja) | 生体状態分析装置及びコンピュータプログラム | |
Hao et al. | Development of Real-Time Driver's Health Detection System by Using Smart Handle |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070928 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100629 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100828 |
|
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: 20100917 |
|
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: 20101013 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131022 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4611206 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |