JP5892678B2 - 生体状態推定装置及びコンピュータプログラム - Google Patents
生体状態推定装置及びコンピュータプログラム Download PDFInfo
- Publication number
- JP5892678B2 JP5892678B2 JP2011108909A JP2011108909A JP5892678B2 JP 5892678 B2 JP5892678 B2 JP 5892678B2 JP 2011108909 A JP2011108909 A JP 2011108909A JP 2011108909 A JP2011108909 A JP 2011108909A JP 5892678 B2 JP5892678 B2 JP 5892678B2
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- region
- time
- reference point
- time series
- 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
Links
- 238000004590 computer program Methods 0.000 title claims description 21
- 230000008859 change Effects 0.000 claims description 186
- 238000004458 analytical method Methods 0.000 claims description 129
- 238000000034 method Methods 0.000 claims description 122
- 238000004364 calculation method Methods 0.000 claims description 116
- 238000001514 detection method Methods 0.000 claims description 101
- 230000033001 locomotion Effects 0.000 claims description 44
- 238000005259 measurement Methods 0.000 claims description 35
- 230000003044 adaptive effect Effects 0.000 claims description 23
- 230000000737 periodic effect Effects 0.000 claims description 22
- 230000006978 adaptation Effects 0.000 claims description 21
- 230000003595 spectral effect Effects 0.000 claims description 19
- 238000012731 temporal analysis Methods 0.000 claims description 17
- 238000000700 time series analysis Methods 0.000 claims description 17
- 230000003247 decreasing effect Effects 0.000 claims description 12
- 230000001965 increasing effect Effects 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims 1
- 238000012360 testing method Methods 0.000 description 114
- 239000004744 fabric Substances 0.000 description 103
- 239000007787 solid Substances 0.000 description 74
- 238000010586 diagram Methods 0.000 description 72
- 239000006260 foam Substances 0.000 description 38
- 230000002889 sympathetic effect Effects 0.000 description 33
- 230000006870 function Effects 0.000 description 28
- 210000005036 nerve Anatomy 0.000 description 23
- 210000005037 parasympathetic nerve Anatomy 0.000 description 20
- 230000003068 static effect Effects 0.000 description 19
- -1 polyethylene terephthalate Polymers 0.000 description 18
- 239000011324 bead Substances 0.000 description 16
- 241000282414 Homo sapiens Species 0.000 description 15
- 230000006698 induction Effects 0.000 description 14
- 238000009826 distribution Methods 0.000 description 13
- 239000012528 membrane Substances 0.000 description 13
- 210000003205 muscle Anatomy 0.000 description 13
- 206010062519 Poor quality sleep Diseases 0.000 description 12
- 201000010099 disease Diseases 0.000 description 12
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 12
- 239000000835 fiber Substances 0.000 description 12
- 238000011156 evaluation Methods 0.000 description 11
- 230000001953 sensory effect Effects 0.000 description 11
- 230000007423 decrease Effects 0.000 description 10
- 230000008035 nerve activity Effects 0.000 description 10
- 229920000139 polyethylene terephthalate Polymers 0.000 description 9
- 239000005020 polyethylene terephthalate Substances 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 210000000709 aorta Anatomy 0.000 description 8
- 206010022000 influenza Diseases 0.000 description 8
- 230000001734 parasympathetic effect Effects 0.000 description 8
- 210000003403 autonomic nervous system Anatomy 0.000 description 7
- 230000013632 homeostatic process Effects 0.000 description 7
- 239000000463 material Substances 0.000 description 7
- 230000002093 peripheral effect Effects 0.000 description 7
- 206010041349 Somnolence Diseases 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 6
- 210000002837 heart atrium Anatomy 0.000 description 6
- 230000029058 respiratory gaseous exchange Effects 0.000 description 6
- 208000024891 symptom Diseases 0.000 description 6
- 208000032140 Sleepiness Diseases 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 5
- 230000003321 amplification Effects 0.000 description 5
- 230000001746 atrial effect Effects 0.000 description 5
- 230000006835 compression Effects 0.000 description 5
- 238000007906 compression Methods 0.000 description 5
- 238000003199 nucleic acid amplification method Methods 0.000 description 5
- 229920002215 polytrimethylene terephthalate Polymers 0.000 description 5
- 230000002040 relaxant effect Effects 0.000 description 5
- 230000037321 sleepiness Effects 0.000 description 5
- 230000007704 transition Effects 0.000 description 5
- 230000037007 arousal Effects 0.000 description 4
- 238000009499 grossing Methods 0.000 description 4
- 230000001939 inductive effect Effects 0.000 description 4
- 238000012423 maintenance Methods 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 4
- 230000002861 ventricular Effects 0.000 description 4
- 210000001217 buttock Anatomy 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 3
- 238000005187 foaming Methods 0.000 description 3
- 230000036541 health Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000036651 mood Effects 0.000 description 3
- 230000010355 oscillation Effects 0.000 description 3
- 230000008447 perception Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000000644 propagated effect Effects 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 230000002441 reversible effect Effects 0.000 description 3
- 230000035807 sensation Effects 0.000 description 3
- 210000000689 upper leg Anatomy 0.000 description 3
- 241001465754 Metazoa Species 0.000 description 2
- 239000004698 Polyethylene Substances 0.000 description 2
- 239000004743 Polypropylene Substances 0.000 description 2
- 210000001015 abdomen Anatomy 0.000 description 2
- 230000016571 aggressive behavior Effects 0.000 description 2
- 210000002376 aorta thoracic Anatomy 0.000 description 2
- 210000001765 aortic valve Anatomy 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 2
- 230000000994 depressogenic effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 206010025482 malaise Diseases 0.000 description 2
- 239000004745 nonwoven fabric Substances 0.000 description 2
- 229920000573 polyethylene Polymers 0.000 description 2
- 229920001155 polypropylene Polymers 0.000 description 2
- 229920005989 resin Polymers 0.000 description 2
- 239000011347 resin Substances 0.000 description 2
- 206010003658 Atrial Fibrillation Diseases 0.000 description 1
- 206010016352 Feeling of relaxation Diseases 0.000 description 1
- 241000282412 Homo Species 0.000 description 1
- 125000002066 L-histidyl group Chemical group [H]N1C([H])=NC(C([H])([H])[C@](C(=O)[*])([H])N([H])[H])=C1[H] 0.000 description 1
- 239000004793 Polystyrene Substances 0.000 description 1
- 206010039897 Sedation Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 210000000467 autonomic pathway Anatomy 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- ZPUCINDJVBIVPJ-LJISPDSOSA-N cocaine Chemical compound O([C@H]1C[C@@H]2CC[C@@H](N2C)[C@H]1C(=O)OC)C(=O)C1=CC=CC=C1 ZPUCINDJVBIVPJ-LJISPDSOSA-N 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 210000004177 elastic tissue Anatomy 0.000 description 1
- 230000008713 feedback mechanism Effects 0.000 description 1
- 230000003862 health status Effects 0.000 description 1
- 208000019622 heart disease Diseases 0.000 description 1
- 230000003284 homeostatic effect Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000007794 irritation Effects 0.000 description 1
- 238000005304 joining Methods 0.000 description 1
- 238000007562 laser obscuration time method Methods 0.000 description 1
- WABPQHHGFIMREM-UHFFFAOYSA-N lead(0) Chemical compound [Pb] WABPQHHGFIMREM-UHFFFAOYSA-N 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 230000003340 mental effect Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000474 nursing effect Effects 0.000 description 1
- 210000001002 parasympathetic nervous system Anatomy 0.000 description 1
- 239000002985 plastic film Substances 0.000 description 1
- 229920006255 plastic film Polymers 0.000 description 1
- 229920000728 polyester Polymers 0.000 description 1
- 239000011112 polyethylene naphthalate Substances 0.000 description 1
- 229920002223 polystyrene Polymers 0.000 description 1
- 229920003225 polyurethane elastomer Polymers 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 230000036280 sedation Effects 0.000 description 1
- 210000000697 sensory organ Anatomy 0.000 description 1
- 230000003860 sleep quality Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 210000002820 sympathetic nervous system Anatomy 0.000 description 1
- 229920003002 synthetic resin Polymers 0.000 description 1
- 239000000057 synthetic resin Substances 0.000 description 1
- KKEYFWRCBNTPAC-UHFFFAOYSA-L terephthalate(2-) Chemical compound [O-]C(=O)C1=CC=C(C([O-])=O)C=C1 KKEYFWRCBNTPAC-UHFFFAOYSA-L 0.000 description 1
- 229920001169 thermoplastic Polymers 0.000 description 1
- 239000004416 thermosoftening plastic Substances 0.000 description 1
- 239000002759 woven fabric Substances 0.000 description 1
Description
前記状態推定手段は、前記ゼロクロス検出手段を用いた周波数の時系列波形から得られる第1の判定基準点と、前記ピーク検出手段を用いた周波数の時系列波形から得られる第2の判定基準点とを用い、前記第1の判定基準点に基づく指標を一方の軸に、前記第2の判定基準点に基づく指標を他方の軸にとり、第1の判定基準点と第2の判定基準点とから求められる座標の時系列変化を求め、生体の状態を推定する構成であることが好ましい。
活性・適応領域及び耐性・適応領域間である場合に、体調良好と推定し、
活性・抵抗領域及び耐性・適応領域間である場合に、通常状態と推定し、
耐性・抵抗領域及び活性・適応領域間である場合に、体調の急変のおそれのある状態と推定する構成とすることが好ましい。
さらに、ピーク検出手段を用いた周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に周波数の平均値を求める移動計算を行い、時間窓毎に得られる周波数の平均値の時系列変化を周波数変動時系列波形として出力する周波数変動演算手段を有し、前記状態推定手段は、前記ゼロクロス検出手段を用いた周波数の時系列波形から求められる前記機能点に対応する指標を一方の軸にとると共に、前記周波数変動演算手段により求められる周波数変動時系列波形の所定の時間幅における変化量に対応する指標を他方の軸にとり、前記機能点と前記変化量とから求められる座標の時系列変化を求め、感覚に関する生体の状態を推定する構成とすることが好ましい。
機能点=後時間範囲の判定基準点+(後時間範囲の判定基準点−前時間範囲の判定基準点)×n、(但し、nは補正係数)
により求められる機能点を時系列に求め、機能点の時系列変化から、生体の状態を推定する構成であることが好ましい。
前記第1の判定基準点に基づく指標を一方の軸に、前記第2の判定基準点に基づく指標を他方の軸にとり、
第1の判定基準点と第2の判定基準点とから求められる座標の時系列変化を求め、生体の状態を推定することが好ましい。
活性・適応領域及び耐性・適応領域間である場合に、体調良好と推定し、
活性・抵抗領域及び耐性・適応領域間である場合に、通常状態と推定し、
耐性・抵抗領域及び活性・適応領域間である場合に、体調の急変のおそれのある状態と推定することが好ましい。
さらに、ピーク検出手順を用いた周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に周波数の平均値を求める移動計算を行い、時間窓毎に得られる周波数の平均値の時系列変化を周波数変動時系列波形として出力する周波数変動演算手順を有し、前記状態推定手順は、前記ゼロクロス検出手順を用いた周波数の時系列波形から求められる前記機能点に対応する指標を一方の軸にとると共に、前記周波数変動演算手順により求められる周波数変動時系列波形の所定の時間幅における変化量に対応する指標を他方の軸にとり、前記機能点と前記変化量とから求められる座標の時系列変化を求め、感覚に関する生体の状態を推定することが好ましい。
材質:
表側のグランド編地・・・300デシテックス/288fのポリエチレンテレフタレート繊維仮撚加工糸と700デシテックス/192fのポリエチレンテレフタレート繊維仮撚加工糸との撚り糸
裏側のグランド編地・・・450デシテックス/108fのポリエチレンテレフタレート繊維仮撚加工糸と350デシテックス/1fのポリトリメチレンテレフタレートモノフィラメントとの組み合わせ
連結糸・・・・・・・・・350デシテックス/1fのポリトリメチレンテレフタレートモノフィラメント
材質:
グランド編地(縦糸)・・・600デシテックス/192fのポリエチレンテレフタレート繊維仮撚加工糸
グランド編地(横糸)・・・300デシテックス/72fのポリエチレンテレフタレート繊維仮撚加工糸
連結糸・・・・・・・・・800デシテックス/1fのポリエチレンテレフタレートモノフィラメント
材質:
表側のグランド編地・・・450デシテックス/108fのポリエチレンテレフタレート繊維仮撚加工糸の2本の撚り糸
裏側のグランド編地・・・450デシテックス/108fのポリエチレンテレフタレート繊維仮撚加工糸の2本の撚り糸
連結糸・・・・・・・・・350デシテックス/1fのポリトリメチレンテレフタレートモノフィラメント
材質:
表側のグランド編地・・・450デシテックス/144fのポリエチレンテレフタレート繊維仮撚加工糸の2本の撚り糸
裏側のグランド編地・・・450デシテックス/144fのポリエチレンテレフタレート繊維仮撚加工糸と350デシテックス/1fのポリトリメチレンテレフタレートモノフィラメントとの組み合わせ
連結糸・・・・・・・・・350デシテックス/1fのポリトリメチレンテレフタレートモノフィラメント
比較対象の前後2つの時間範囲における解析波形の判定基準点間において、次式:
機能点=後時間範囲の判定基準点+(後時間範囲の判定基準点−前時間範囲の判定基準点)×n、(但し、nは補正係数)
により求められる機能点を時系列に求めていく手段である。
交感神経系の活動状態を示すゼロクロス検出手段を用いた周波数の時系列波形から得られる第1の判定基準点と、副交感神経系の活動状態を示すピーク検出手段を用いた周波数の時系列波形から得られる第2の判定基準点とを用い、第1の判定基準点に基づく指標を一方の軸に、第2の判定基準点に基づく指標を他方の軸にとり、第1の判定基準点と第2の判定基準点とから求められる全身状態の変化の様子を座標の時系列変化として求める手段である。
なお、図17(b)及び後述の各試験例の判定結果(図16、図23、図32、図41、図52)において示した座標系では、横軸のプラス側を「適応」、マイナス側を「抵抗」と設定し、縦軸のプラス側を「活性」、マイナス側を「耐性」と設定しているが、これを、心理学上の用語で表現すると次のようになる。すなわち、交感神経活動に関連する横軸は、プラス方向に値が大きくなるほど「プラスの「発揮」大」となり、マイナス方向に値が大きくなるほど「マイナスの「発揮」大」となる。副交感神経活動に関連する縦軸は、プラス方向に値が大きくなるほど「プラスの「求められている」大」となり、マイナス方向に値が大きくなるほど「マイナスの「求められている」大」となる。そして、座標原点(0,0)に近づくほど、「発揮」及び「求められている」はいずれも小さくなる。このように心理学上の用語で設定すると、例えば、「ハツラツ領域」は「プラスの「発揮」大」でかつ「プラスの「求められている」大」の領域となり、「ハツラツ」あるいは「集中・高揚・緊張・喜び」として表現した領域の意味合いにマッチする。
図1に示した生体信号測定手段1を、(株)デルタツーリング製、シートクッションである商品名「ツインランバー」の背部の裏側に積層し、自動車用シートに取り付け、被験者を着座させ、座位姿勢での心房や心室及び大動脈の揺動による生体信号、いわゆる体表脈波(以下、「心部揺動波」というが、「APW」と略記する場合もある)を採取した。なお、生体信号測定手段1を構成する板状発泡体21,22及び三次元立体編物支持部材15は、ビーズの平均直径が約5mmで、厚さ3mmにスライスカットしたビーズ発泡体を用いた。三次元立体編物10は、住江織物(株)製、製品番号:49011Dで、厚さ10mmのものであった。フィルム16は、シーダム株式会社製、品番「DUS605−CDR」を用いた。被験者は、70歳代の健康な男性である。また、上記自動車用シートを助手席に搭載して、停止状態の静的条件下で測定する試験と、その後自動車を走行させながら測定した実車走行試験とを行った。
図1に示した生体信号測定手段1を、(株)デルタツーリング製、シートクッションである商品名「ツインランバー」の背部の裏側に積層し、自動車用シートに取り付け、被験者を着座させ、座位姿勢(覚醒誘導姿勢(図では「覚醒姿勢」と表示))で心房や心室及び大動脈の揺動による生体信号(以下、「心部揺動波」というが、「APW」と略記する場合もある)を採取した。また、シートバックを通常の座位姿勢(覚醒誘導姿勢)よりも大腿部と脊柱の開度を広げた座位姿勢(リラックス姿勢(図では「寝姿勢」と表示))でも生体信号(APW)を採取した。生体信号測定手段1のその他の構成は、試験例1のものと同じである、また、被験者は40歳代の健康な男性である。なお、図25及び図26は、試験中に測定した脳波、指尖容積脈波、心拍による解析結果を示したものであり、図25は覚醒誘導姿勢の結果を、図26はリラックス姿勢の結果を示し、覚醒誘導姿勢では試験中、被験者は覚醒状態であり、リラックス姿勢では、試験中盤と終盤において睡眠状態であった。
ゼロクロス検出手段を用いた第1判定基準点は、図18(a)では3点、図18(b)では6点、図18(c)では3点、図18(d)では3点であった。ピーク検出手段を用いた第2判定基準点は、図19(a)では4点、図19(b)では4点、図19(c)では3点、図19(d)では4点であった。
ゼロクロス検出手段を用いた第1判定基準点は、図20(a)では4点、図20(b)では5点、図20(c)では9点、図20(d)では6点であった。ピーク検出手段を用いた第2判定基準点は、図21(a)では3点、図21(b)では8点、図21(c)では7点、図21(d)では4点であった。
20歳代男性被験者について、試験例2と同じ条件で、通常の座位姿勢(覚醒誘導姿勢(図では「覚醒姿勢」と表示))及び大腿部と脊柱の開度を広げた座位姿勢(リラックス姿勢(図では「寝姿勢」と表示))で生体信号(APW)を採取した。なお、図34及び図35は、試験中に測定した脳波、指尖容積脈波、心拍による解析結果を示したものであり、図34は覚醒誘導姿勢の結果を、図35はリラックス姿勢の結果を示すが、これらの図から、被験者は、いずれも前半は若干うとうとしている傾向が見られるが、後半は覚醒していることがわかる。
ゼロクロス検出手段を用いた第1判定基準点は、図27(a)では6点、図27(b)では4点、図27(c)では5点、図27(d)では6点、図27(e)では8点であった。ピーク検出手段を用いた第2判定基準点は、図28(a)では4点、図28(b)では4点、図28(c)では7点、図28(d)では7点、図28(e)では4点であった。
ゼロクロス検出手段を用いた第1判定基準点は、図29(a)では7点、図29(b)では4点、図29(c)では2点、図29(d)では8点、図29(e)では3点であった。ピーク検出手段を用いた第2判定基準点は、図30(a)では6点、図30(b)では4点、図30(c)では2点、図30(d)では1点、図30(e)では4点であった。
図36〜図46は試験例3と同じ20歳代男性被験者のデータである。このうち、図36及び図37は、インフルエンザ発症時のデータであるが、これは、試験例3の図27及び図28と同じデータである。比較しやすくするために再掲している。図38及び図39は、インフルエンザが治癒し、十分に体力が回復した後(試験例3の測定日から2週間経過後)に、覚醒誘導姿勢で測定した健常時のデータである。
40歳代男性被験者について、試験例2,3と同じ条件で、通常の座位姿勢(覚醒誘導姿勢(図では「覚醒姿勢」と表示))及び大腿部と脊柱の開度を広げた座位姿勢(リラックス姿勢(図では「寝姿勢」と表示))でも生体信号(APW)を採取した。なお、図54及び図55は、試験中に測定した脳波、指尖容積脈波、心拍による解析結果を示したものであり、図54は覚醒誘導姿勢の結果を、図55はリラックス姿勢の結果を示すが、これらの図から、被験者は、いずれも試験中浅い眠りに陥った状態であった。
ゼロクロス検出手段を用いた第1判定基準点は、図47(a)では7点、図47(b)では6点、図47(c)では6点、図47(d)では6点、図47(e)では4点であった。ピーク検出手段を用いた第2判定基準点は、図48(a)では8点、図48(b)では3点、図48(c)では3点、図48(d)では1点、図48(e)では3点であった。
ゼロクロス検出手段を用いた第1判定基準点は、図49(a)では5点、図49(b)では4点、図49(c)では4点、図49(d)では4点、図49(e)では4点であった。ピーク検出手段を用いた第2判定基準点は、図50(a)では7点、図50(b)では3点、図50(c)では3点、図50(d)では6点、図50(e)では4点であった。
図1に示した生体信号測定手段1を、(株)デルタツーリング製、シートクッションである商品名「ツインランバー」の背部の裏側に積層し、トラックの運転席に取り付け、運転中の被験者の心房や心室及び大動脈の揺動による生体信号(以下、「心部揺動波」というが、「APW」と略記する場合もある)を採取した。被験者数は、延べ153名であった。
次に、人がどう感じているのかという感覚の変化、知覚感覚器の変化の様子を主として捉える方法について説明する。ここでは、試験例1のデータを用い、状態推定手段660により、次のような処理を行った。その結果が、図62(c)、図63(c)、図64(c)及び図65(c)である。これらの図は、図16に示された結果を加味すると共に、ゼロクロス検出手段を周波数解析して得られた機能点(図62〜図65の(a))と、ピーク検出手段により得られた周波数変動の時系列波形(図62〜図65の(b))から求められるグラフの変化量(傾き)とを用いて示したものである。
図62〜図65が、図16の静的状態、「小谷SA〜与島PA」間の走行時、「与島PA〜山越」間の走行時、「山越〜石鎚山SA」間の走行時にそれぞれ対応する。
10 三次元立体編物
15 三次元立体編物支持部材
15a 配置用貫通孔
16 フィルム
21,22 板状発泡体
30 振動センサ
100 シート
110 シートバックフレーム
120 表皮
60 生体状態推定装置
610 周波数演算手段
620 周波数傾き時系列解析演算手段
630 周波数解析手段
640 回帰直線演算手段
650 判定基準点算出手段
660 状態推定手段
Claims (26)
- 生体信号測定手段により採取した生体信号を用いて、生体の状態を推定する生体状態推定装置であって、
前記生体信号測定手段により得られる所定の測定時間における生体信号の時系列波形から、周波数の時系列波形を求める周波数演算手段と、
前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列変化を周波数傾き時系列波形として出力する周波数傾き時系列解析演算手段と、
前記周波数傾き時系列解析演算手段から得られる所定時間範囲における周波数傾き時系列波形を周波数解析し、パワースペクトル密度と周波数との関係を示す解析波形を所定時間範囲毎に出力する周波数解析手段と、
前記周波数解析手段により出力される各解析波形について、所定周期領域毎に回帰直線を求める回帰直線演算手段と、
前記周期領域毎に求められる各回帰直線を、その傾きに基づいて領域得点を付与すると共に、隣接する周波数領域における回帰直線間のパワースペクトル密度の値の較差及び回帰直線間の傾きの違いに基づき、各回帰直線全体における分岐現象を示す折れ点数を求め、その折れ点数に基づいた形状得点を付与し、前記領域得点及び形状得点の少なくとも一方を用いて、各解析波形についての判定基準点を求める判定基準点算出手段と、
前記判定基準点算出手段により求められる前記各解析波形の判定基準点の時系列の変化を基に、生体の状態を推定する状態推定手段と
を具備することを特徴とする生体状態推定装置。 - 前記回帰直線演算手段は、分析対象の解析波形を、長周期領域、中周期領域及び短周期領域に分けて、前記回帰直線を求める請求項1記載の生体状態推定装置。
- 前記回帰直線演算手段は、前記長周期領域においては、前記長周期領域の中心周波数を境界としてULF領域とVLF領域との2つに分けてそれぞれについて回帰直線を求め、ULF領域の回帰直線及びVLF領域の回帰直線の各傾きの積が所定値以下か否かを判定し、所定値以下の場合には前記ULF領域及びVLF領域の各回帰直線を採用し、所定値を上回る場合には前記長周期領域の全体における前記回帰直線を採用する請求項2記載の生体状態推定装置。
- 前記判定基準点算出手段は、前記領域得点として、前記各領域における各回帰直線の傾きを略水平状態、上向き及び下向きの3つに分け、略水平状態の得点を基準として、上向きの場合と下向きの場合とで得点を増減させる構成である請求項1〜3のいずれか1に記載の生体状態推定装置。
- 前記判定基準点算出手段は、前記形状得点として、前記折れ点数が少ないほど高得点を付与する構成である請求項1〜4のいずれか1に記載の生体状態推定装置。
- 前記判定基準点算出手段は、前記折れ点数を、隣接する周期領域の2つの回帰直線間において、パワースペクトル密度の値の較差が所定以上の場合、及び、隣接する周期領域の2つの回帰直線間において、パワースペクトル密度の値の較差が所定以内であって2つの回帰直線の傾きの角度の違いが予め設定した所定角度以上の場合に、それぞれ折れ点としてカウントする請求項1〜5のいずれか1に記載の生体状態推定装置。
- 前記状態推定手段は、比較対象の前後2つの時間範囲における解析波形の判定基準点間において、次式:
機能点=後時間範囲の判定基準点+(後時間範囲の判定基準点−前時間範囲の判定基準点)×n、(但し、nは補正係数)
により求められる機能点を時系列に求め、機能点の時系列変化から、生体の状態を推定する請求項1〜6のいずれか1に記載の生体状態推定装置。 - 前記周波数演算手段は、前記生体信号の時系列波形におけるゼロクロス点を用いて周波数の時系列波形を求めるゼロクロス検出手段と、前記生体信号の時系列波形のピーク点を用いて周波数の時系列波形を求めるピーク検出手段とのいずれか少なくとも一方の手段を備える請求項1〜7のいずれか1に記載の生体状態推定装置。
- 前記状態推定手段は、前記ゼロクロス検出手段を用いた周波数の時系列波形から得られる第1の判定基準点と、前記ピーク検出手段を用いた周波数の時系列波形から得られる第2の判定基準点とを用い、
前記第1の判定基準点に基づく指標を一方の軸に、前記第2の判定基準点に基づく指標を他方の軸にとり、
第1の判定基準点と第2の判定基準点とから求められる座標の時系列変化を求め、生体の状態を推定する請求項8記載の生体状態推定装置。 - 前記状態推定手段は、前記座標同士を結んだ座標時系列変化線が、1/fの傾きに近似した変化傾向であると判定された場合には快適と判定し、上下方向に変化していると判定された場合には不快と判定する請求項9記載の生体状態推定装置。
- 前記状態推定手段は、前記第1の判定基準点に基づく指標を一方の軸に、前記第2の判定基準点に基づく指標を他方の軸にとった座標系を象限毎に、活性・適応領域、活性・抵抗領域、耐性・抵抗領域、耐性・適応領域に区分し、異なる測定時間において求められた複数の前記座標時系列変化線同士を比較した場合に、比較対象の前記座標時系列変化線の全体の移動方向により、体調を推定する手段を有する請求項10記載の生体状態推定装置。
- 前記座標時系列変化線の全体の主な移動方向が、
活性・適応領域及び耐性・適応領域間である場合に、体調良好と推定し、
活性・抵抗領域及び耐性・適応領域間である場合に、通常状態と推定し、
耐性・抵抗領域及び活性・適応領域間である場合に、体調の急変のおそれのある状態と推定する請求項10記載の生体状態推定装置。 - 前記周波数演算手段は、前記生体信号の時系列波形におけるゼロクロス点を用いて周波数の時系列波形を求めるゼロクロス検出手段と、前記生体信号の時系列波形のピーク点を用いて周波数の時系列波形を求めるピーク検出手段とを備えており、
前記周波数演算手段におけるピーク検出手段を用いた周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に周波数の平均値を求める移動計算を行い、時間窓毎に得られる周波数の平均値の時系列変化を周波数変動時系列波形として出力する周波数変動演算手段をさらに有し、
前記状態推定手段が、前記ゼロクロス検出手段を用いた周波数の時系列波形から求められる前記機能点に対応する指標を一方の軸にとると共に、前記周波数変動演算手段により求められる周波数変動時系列波形の所定の時間幅における変化量に対応する指標を他方の軸にとり、
前記機能点と前記変化量とから求められる座標の時系列変化を求め、感覚に関する生体の状態を推定する請求項7記載の生体状態推定装置。 - 生体信号測定手段により採取した生体信号を用いて、生体の状態を推定する生体状態推定装置に設定されるコンピュータプログラムであって、
前記生体信号測定手段により得られる所定の測定時間における生体信号の時系列波形から、周波数の時系列波形を求める周波数演算手順と、
前記周波数演算手順により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列変化を周波数傾き時系列波形として出力する周波数傾き時系列解析演算手順と、
前記周波数傾き時系列解析演算手順から得られる所定時間範囲における周波数傾き時系列波形を周波数解析し、パワースペクトル密度と周波数との関係を示す解析波形を所定時間範囲毎に出力する周波数解析手順と、
前記周波数解析手順により出力される各解析波形について、所定周期領域毎に回帰直線を求める回帰直線演算手順と、
前記周期領域毎に求められる各回帰直線を、その傾きに基づいて領域得点を付与すると共に、隣接する周波数領域における回帰直線間のパワースペクトル密度の値の較差及び回帰直線間の傾きの違いに基づき、各回帰直線全体における分岐現象を示す折れ点数を求め、その折れ点数に基づいた形状得点を付与し、前記領域得点及び形状得点の少なくとも一方を用いて、各解析波形についての判定基準点を求める判定基準点算出手順と、
前記判定基準点算出手順により求められる前記各解析波形の判定基準点の時系列の変化を基に、生体の状態を推定する状態推定手順と
をコンピュータに実行させるコンピュータプログラム。 - 前記回帰直線演算手順は、分析対象の解析波形を、長周期領域、中周期領域及び短周期領域に分けて、前記回帰直線を求める請求項14記載のコンピュータプログラム。
- 前記回帰直線演算手順は、前記長周期領域においては、前記長周期領域の中心周波数を境界としてULF領域とVLF領域との2つに分けてそれぞれについて回帰直線を求め、ULF領域の回帰直線及びVLF領域の回帰直線の各傾きの積が所定値以下か否かを判定し、所定値以下の場合には前記ULF領域及びVLF領域の各回帰直線を採用し、所定値を上回る場合には前記長周期領域の全体における前記回帰直線を採用する請求項15記載のコンピュータプログラム。
- 前記判定基準点算出手順は、前記領域得点として、前記各領域における各回帰直線の傾きを略水平状態、上向き及び下向きの3つに分け、略水平状態の得点を基準として、上向きの場合と下向きの場合とで得点を増減させる請求項14〜16のいずれか1に記載のコンピュータプログラム。
- 前記判定基準点算出手順は、前記形状得点として、前記折れ点数が少ないほど高得点を付与する請求項14〜17のいずれか1に記載のコンピュータプログラム。
- 前記判定基準点算出手順は、前記折れ点数を、隣接する周期領域の2つの回帰直線間において、パワースペクトル密度の値の較差が所定以上の場合、及び、隣接する周期領域の2つの回帰直線間において、パワースペクトル密度の値の較差が所定以内であって2つの回帰直線の傾きの角度の違いが予め設定した所定角度以上の場合に、それぞれ折れ点としてカウントする請求項14〜18のいずれか1に記載のコンピュータプログラム。
- 前記状態推定手順は、比較対象の前後2つの時間範囲における解析波形の判定基準点間において、次式:
機能点=後時間範囲の判定基準点+(後時間範囲の判定基準点−前時間範囲の判定基準点)×n、(但し、nは補正係数)
により求められる機能点を時系列に求め、機能点の時系列変化から、生体の状態を推定する請求項14〜19のいずれか1に記載のコンピュータプログラム。 - 前記周波数演算手順は、前記生体信号の時系列波形におけるゼロクロス点を用いて周波数の時系列波形を求めるゼロクロス検出手順と、前記生体信号の時系列波形のピーク点を用いて周波数の時系列波形を求めるピーク検出手順とのいずれか少なくとも一方の手順を備える請求項14〜20のいずれか1に記載のコンピュータプログラム。
- 前記状態推定手順は、前記ゼロクロス検出手順を用いた周波数の時系列波形から得られる第1の判定基準点と、前記ピーク検出手順を用いた周波数の時系列波形から得られる第2の判定基準点とを用い、
前記第1の判定基準点に基づく指標を一方の軸に、前記第2の判定基準点に基づく指標を他方の軸にとり、
第1の判定基準点と第2の判定基準点とから求められる座標の時系列変化を求め、生体の状態を推定する請求項21記載のコンピュータプログラム。 - 前記状態推定手順は、前記座標同士を結んだ座標時系列変化線が、1/fの傾きに近似した変化傾向であると判定された場合には快適と判定し、上下方向に変化していると判定された場合には不快と判定する請求項22記載のコンピュータプログラム。
- 前記状態推定手順は、前記第1の判定基準点に基づく指標を一方の軸に、前記第2の判定基準点に基づく指標を他方の軸にとった座標系を象限毎に、活性・適応領域、活性・抵抗領域、耐性・抵抗領域、耐性・適応領域に区分し、異なる測定時間において求められた複数の前記座標時系列変化線同士を比較した場合に、比較対象の前記座標時系列変化線の全体の移動方向により、体調を推定する請求項23記載のコンピュータプログラム。
- 前記座標時系列変化線の全体の主な移動方向が、
活性・適応領域及び耐性・適応領域間である場合に、体調良好と推定し、
活性・抵抗領域及び耐性・適応領域間である場合に、通常状態と推定し、
耐性・抵抗領域及び活性・適応領域間である場合に、体調の急変のおそれのある状態と推定する請求項23記載のコンピュータプログラム。 - 前記周波数演算手順は、前記生体信号の時系列波形におけるゼロクロス点を用いて周波数の時系列波形を求めるゼロクロス検出手順と、前記生体信号の時系列波形のピーク点を用いて周波数の時系列波形を求めるピーク検出手順とを備えており、
前記周波数演算手順におけるピーク検出手順を用いた周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に周波数の平均値を求める移動計算を行い、時間窓毎に得られる周波数の平均値の時系列変化を周波数変動時系列波形として出力する周波数変動演算手順をさらに有し、
前記状態推定手順が、前記ゼロクロス検出手順を用いた周波数の時系列波形から求められる前記機能点に対応する指標を一方の軸にとると共に、前記周波数変動演算手順により求められる周波数変動時系列波形の所定の時間幅における変化量に対応する指標を他方の軸にとり、
前記機能点と前記変化量とから求められる座標の時系列変化を求め、感覚に関する生体の状態を推定する請求項20記載のコンピュータプログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011108909A JP5892678B2 (ja) | 2011-05-14 | 2011-05-14 | 生体状態推定装置及びコンピュータプログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011108909A JP5892678B2 (ja) | 2011-05-14 | 2011-05-14 | 生体状態推定装置及びコンピュータプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2012239480A JP2012239480A (ja) | 2012-12-10 |
JP5892678B2 true JP5892678B2 (ja) | 2016-03-23 |
Family
ID=47461855
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011108909A Active JP5892678B2 (ja) | 2011-05-14 | 2011-05-14 | 生体状態推定装置及びコンピュータプログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5892678B2 (ja) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6118097B2 (ja) * | 2012-12-14 | 2017-04-19 | 株式会社デルタツーリング | 運転時生体状態判定装置及びコンピュータプログラム |
JP2016112144A (ja) * | 2014-12-12 | 2016-06-23 | 株式会社デルタツーリング | 生体状態分析装置及びコンピュータプログラム |
JP6666705B2 (ja) * | 2015-09-28 | 2020-03-18 | デルタ工業株式会社 | 生体状態推定装置、生体状態推定方法及びコンピュータプログラム |
WO2017057022A1 (ja) * | 2015-09-28 | 2017-04-06 | デルタ工業株式会社 | 生体状態推定装置、生体状態推定方法及びコンピュータプログラム |
US20180360315A1 (en) | 2015-12-12 | 2018-12-20 | Delta Kogyo Co., Ltd. | Biological state estimation device, biological state estimation method, computer program, and recording medium |
WO2017099257A1 (ja) | 2015-12-12 | 2017-06-15 | デルタ工業株式会社 | 生体状態推定装置、生体状態推定方法、コンピュータプログラム及び記録媒体 |
JP6836264B2 (ja) * | 2016-12-09 | 2021-02-24 | デルタ工業株式会社 | 生体状態推定装置、生体状態推定方法、コンピュータプログラム及び記録媒体 |
JP6836265B2 (ja) * | 2016-12-09 | 2021-02-24 | デルタ工業株式会社 | 生体状態推定装置、生体状態推定方法、コンピュータプログラム及び記録媒体 |
JP7217576B2 (ja) * | 2018-11-05 | 2023-02-03 | 株式会社安藤・間 | ドライバーの状態推定方法及び装置 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2816799B2 (ja) * | 1993-08-10 | 1998-10-27 | 株式会社日立製作所 | センシング装置と環境制御機器及び環境制御方法 |
JP3310498B2 (ja) * | 1994-09-02 | 2002-08-05 | 独立行政法人産業技術総合研究所 | 生体情報解析装置および生体情報解析方法 |
JP4347621B2 (ja) * | 2003-05-21 | 2009-10-21 | 株式会社デルタツーリング | 生体評価システム、コンピュータプログラム及び記録媒体 |
JP4247055B2 (ja) * | 2003-05-21 | 2009-04-02 | 株式会社デルタツーリング | 運転席用座席システム |
JP4637098B2 (ja) * | 2004-03-25 | 2011-02-23 | 株式会社デルタツーリング | 負荷体状態判定装置、乗物用シート及びコンピュータプログラム |
JP4636861B2 (ja) * | 2004-11-25 | 2011-02-23 | 株式会社デルタツーリング | 快適感評価装置 |
JP4959178B2 (ja) * | 2005-02-28 | 2012-06-20 | 株式会社デルタツーリング | クッション材及び圧力変動検出装置 |
JP2007022443A (ja) * | 2005-07-20 | 2007-02-01 | Toyota Motor Corp | 車両用内装材 |
JP5044230B2 (ja) * | 2007-02-14 | 2012-10-10 | 株式会社デルタツーリング | 生体信号分析装置、シート及び生体信号分析方法 |
JP5582478B2 (ja) * | 2009-07-16 | 2014-09-03 | 株式会社デルタツーリング | 生体状態推定装置及びコンピュータプログラム |
-
2011
- 2011-05-14 JP JP2011108909A patent/JP5892678B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
JP2012239480A (ja) | 2012-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5892678B2 (ja) | 生体状態推定装置及びコンピュータプログラム | |
JP5553303B2 (ja) | 生体状態推定装置及びコンピュータプログラム | |
JP5733499B2 (ja) | 生体状態推定装置及びコンピュータプログラム | |
JP6460560B2 (ja) | 音・振動情報収集機構及び音・振動情報センシングシステム | |
JP5704651B2 (ja) | 生体状態推定装置、生体状態推定システム及びコンピュータプログラム | |
JP5044230B2 (ja) | 生体信号分析装置、シート及び生体信号分析方法 | |
JP6118097B2 (ja) | 運転時生体状態判定装置及びコンピュータプログラム | |
JP5977938B2 (ja) | 生体状態推定装置及びコンピュータプログラム | |
JP5429946B2 (ja) | 生体状態分析装置及びコンピュータプログラム | |
JP5669287B2 (ja) | 飲酒検知システム及びコンピュータプログラム | |
US20060155175A1 (en) | Biological sensor and support system using the same | |
JP5582478B2 (ja) | 生体状態推定装置及びコンピュータプログラム | |
JP5751475B2 (ja) | 生体状態推定装置及びコンピュータプログラム | |
JP5327584B2 (ja) | 生体状態分析装置、コンピュータプログラム及び記録媒体 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20140510 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20141226 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20150108 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20150306 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20150805 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20150817 |
|
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: 20160204 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20160220 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5892678 Country of ref document: JP 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 |