ホーム > 特許ランキング > 株式会社デルタツーリング
知財求人 - 知財ポータルサイト「IP Force」
特許5704651生体状態推定装置、生体状態推定システム及びコンピュータプログラム
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】5704651
(24)【登録日】2015年3月6日
(45)【発行日】2015年4月22日
(54)【発明の名称】生体状態推定装置、生体状態推定システム及びコンピュータプログラム
(51)【国際特許分類】
A61B 5/18 20060101AFI20150402BHJP
B60K 28/06 20060101ALN20150402BHJP
【ＦＩ】
A61B5/18
!B60K28/06 A
【請求項の数】20
【全頁数】41
(21)【出願番号】特願2011-536174(P2011-536174)
(86)(22)【出願日】2010年10月14日
(86)【国際出願番号】JP2010068068
(87)【国際公開番号】WO2011046178
(87)【国際公開日】20110421
【審査請求日】2013年10月7日
(31)【優先権主張番号】特願2009-237802(P2009-237802)
(32)【優先日】2009年10月14日
(33)【優先権主張国】JP
(73)【特許権者】
【識別番号】594176202
【氏名又は名称】株式会社デルタツーリング
(74)【代理人】
【識別番号】100101742
【弁理士】
【氏名又は名称】麦島 隆
(72)【発明者】
【氏名】藤田 悦則
(72)【発明者】
【氏名】小倉 由美
(72)【発明者】
【氏名】前田 慎一郎
(72)【発明者】
【氏名】落合 直輝
【審査官】
宮川 哲伸
(56)【参考文献】
【文献】
国際公開第２００９／１０４４６０（ＷＯ，Ａ１）
【文献】
国際公開第２００８／１４３２４９（ＷＯ，Ａ１）
【文献】
国際公開第２００５／０９２１９３（ＷＯ，Ａ１）
【文献】
特開２００９−０２２６１０（ＪＰ，Ａ）
(58)【調査した分野】（Int.Cl.，ＤＢ名）
Ａ６１Ｂ ５／００
Ａ６１Ｂ ５／１８
Ｂ６０Ｋ ２８／０６
(57)【特許請求の範囲】
【請求項1】
生体信号測定手段により人の上体から採取した生体信号を用いて、人の状態を推定する生体状態推定装置であって、
前記生体信号を分析して所定時点における人の恒常性維持機能レベルを複数の段階に区分して求める恒常性維持機能レベル演算手段と、
前記恒常性維持機能レベルの各段階を縦軸に時間を横軸にとり、前記恒常性維持機能レベル演算手段により求められる恒常性維持機能レベルを時系列にプロットし、恒常性維持機能レベルのゆらぎ度合いを表示手段にグラフとして表示させる出力手段と
を具備することを特徴とする生体状態推定装置。
前記生体信号を分析して所定時点における人の恒常性維持機能レベルを複数の段階に区分して求める恒常性維持機能レベル演算手段と、
前記恒常性維持機能レベルの各段階を縦軸に時間を横軸にとり、前記恒常性維持機能レベル演算手段により求められる恒常性維持機能レベルを時系列にプロットし、恒常性維持機能レベルのゆらぎ度合いを表示手段にグラフとして表示させる出力手段と
を具備することを特徴とする生体状態推定装置。
【請求項2】
前記恒常性維持機能レベル演算手段は、恒常性維持機能レベルを高活性状態から機能低下状態まで３〜１０段階の範囲で区分して求め、前記出力手段は、高活性状態を最上部とし、機能低下状態を最下部として縦軸に３〜１０段階の範囲に区分して表示させる請求項１記載の生体状態推定装置。
【請求項3】
前記恒常性維持機能レベル演算手段は、恒常性維持機能レベルを高活性状態から機能低下状態まで５段階に区分して求め、前記出力手段は、高活性状態を最上部とし、機能低下状態を最下部として縦軸に５段階に区分して表示させる請求項２記載の生体状態推定装置。
【請求項4】
前記出力手段により表示手段に表示されるグラフから、人の状態の変遷を判定する判定手段をさらに有する請求項１〜３のいずれか１に記載の生体状態推定装置。
【請求項5】
前記判定手段は、前記出力手段により表示手段に表示されるグラフに、恒常性維持機能レベルの２段階以上の急低下が所定回数以上生じた場合に異常状態と判定する異常状態判定手段を有する請求項４記載の生体状態推定装置。
【請求項6】
前記恒常性維持機能レベル演算手段は、
前記生体信号測定手段により得られる生体信号の周波数を求める周波数演算手段と、
前記周波数演算手段により得られた前記生体信号の周波数から、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列波形を求める周波数傾き時系列解析演算手段と、
前記周波数傾き時系列解析演算手段により求めた前記周波数傾きの時系列波形を微分する微分手段と、
前記周波数傾き時系列解析演算手段により求めた前記周波数傾きの時系列波形を積分する積分手段と、
前記周波数傾き時系列解析演算手段により求めた前記周波数傾きの時系列波形の増減から方形波を求める方形波算出手段と、
任意に設定した第１時間帯における前記周波数傾き時系列波形と前記第１時間帯よりも後の第２時間帯における前記周波数傾き時系列波形と間の記述関数及び記述関数振幅値を求める記述関数算出手段と、
前記周波数演算手段により得られる、前記生体信号の時系列波形の極大値を用いた生体信号の周波数の時系列波形と、前記生体信号の時系列波形の正負が切り替わるゼロクロス地点を用いた生体信号の周波数の時系列波形とを用いて、前記周波数傾き時系列解析演算手段により求めた、それぞれの周波数傾きの時系列波形を絶対値処理する絶対値処理手段と、
前記周波数傾き時系列解析演算手段により求めた前記周波数傾き、前記微分手段により求めた微分値、前記積分手段により求めた積分値、前記方形波算出手段により求めた方形波の正負、前記記述関数算出手段により求めた記述関数振幅値、及び前記絶対値処理手段により求めた周波数傾き時系列波形の２つの絶対値のうち、少なくとも１つ以上を用いて恒常性維持機能レベルの段階を求める恒常性維持機能段階算出手段と
を具備してなる請求項１〜５のいずれか１に記載の生体状態推定装置。
前記生体信号測定手段により得られる生体信号の周波数を求める周波数演算手段と、
前記周波数演算手段により得られた前記生体信号の周波数から、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列波形を求める周波数傾き時系列解析演算手段と、
前記周波数傾き時系列解析演算手段により求めた前記周波数傾きの時系列波形を微分する微分手段と、
前記周波数傾き時系列解析演算手段により求めた前記周波数傾きの時系列波形を積分する積分手段と、
前記周波数傾き時系列解析演算手段により求めた前記周波数傾きの時系列波形の増減から方形波を求める方形波算出手段と、
任意に設定した第１時間帯における前記周波数傾き時系列波形と前記第１時間帯よりも後の第２時間帯における前記周波数傾き時系列波形と間の記述関数及び記述関数振幅値を求める記述関数算出手段と、
前記周波数演算手段により得られる、前記生体信号の時系列波形の極大値を用いた生体信号の周波数の時系列波形と、前記生体信号の時系列波形の正負が切り替わるゼロクロス地点を用いた生体信号の周波数の時系列波形とを用いて、前記周波数傾き時系列解析演算手段により求めた、それぞれの周波数傾きの時系列波形を絶対値処理する絶対値処理手段と、
前記周波数傾き時系列解析演算手段により求めた前記周波数傾き、前記微分手段により求めた微分値、前記積分手段により求めた積分値、前記方形波算出手段により求めた方形波の正負、前記記述関数算出手段により求めた記述関数振幅値、及び前記絶対値処理手段により求めた周波数傾き時系列波形の２つの絶対値のうち、少なくとも１つ以上を用いて恒常性維持機能レベルの段階を求める恒常性維持機能段階算出手段と
を具備してなる請求項１〜５のいずれか１に記載の生体状態推定装置。
【請求項7】
前記恒常性維持機能レベル演算手段により処理される前に、前記生体信号測定手段により得られる生体信号を、所定の周波数帯域でフィルタリングするフィルタリング手段をさらに有し、
前記フィルタリング手段が、前記生体信号測定手段により測定した静的環境下における生体信号の平均周波数を求め、該平均周波数を基準値として、フィルタリングする周波数帯域を設定する構成である請求項１〜６のいずれか１に記載の生体状態推定装置。
前記フィルタリング手段が、前記生体信号測定手段により測定した静的環境下における生体信号の平均周波数を求め、該平均周波数を基準値として、フィルタリングする周波数帯域を設定する構成である請求項１〜６のいずれか１に記載の生体状態推定装置。
【請求項8】
前記生体信号測定手段は、人の背部に当接され、心臓の動きによる大動脈の揺れに伴って圧力が変動するエアパックを備えてなり、前記恒常性維持機能レベル演算手段は、前記エアパックの圧力変動を分析対象の前記生体信号として用いる請求項１〜７のいずれか１に記載の生体状態推定装置。
【請求項9】
前記エアパックが、乗物用シートのシートバック部に装着されている請求項８記載の生体状態推定装置。
【請求項10】
前記表示手段が車載用モニタである請求項１〜９のいずれか１に記載の生体状態推定装置。
【請求項11】
人の上体から生体信号を採取する生体信号測定手段と、
前記生体信号測定手段により採取した生体信号を用いて、人の状態を推定する生体状態推定装置とを備えた生体状態推定システムであって、
前記生体状態推定装置は、
前記生体信号を分析して所定時点における人の恒常性維持機能レベルを複数の段階に区分して求める恒常性維持機能レベル演算手段と、
前記恒常性維持機能レベルの各段階を縦軸に時間を横軸にとり、前記恒常性維持機能レベル演算手段により求められる恒常性維持機能レベルを時系列にプロットし、恒常性維持機能レベルのゆらぎ度合いを表示手段にグラフとして表示させる出力手段と
を具備することを特徴とする生体状態推定システム。
前記生体信号測定手段により採取した生体信号を用いて、人の状態を推定する生体状態推定装置とを備えた生体状態推定システムであって、
前記生体状態推定装置は、
前記生体信号を分析して所定時点における人の恒常性維持機能レベルを複数の段階に区分して求める恒常性維持機能レベル演算手段と、
前記恒常性維持機能レベルの各段階を縦軸に時間を横軸にとり、前記恒常性維持機能レベル演算手段により求められる恒常性維持機能レベルを時系列にプロットし、恒常性維持機能レベルのゆらぎ度合いを表示手段にグラフとして表示させる出力手段と
を具備することを特徴とする生体状態推定システム。
【請求項12】
前記生体信号測定手段は、人の背部に当接され、心臓の動きによる大動脈の揺れに伴って圧力が変動するエアパックを備えてなり、前記恒常性維持機能レベル演算手段は、前記エアパックの圧力変動を分析対象の前記生体信号として用いる請求項１１記載の生体状態推定システム。
【請求項13】
前記エアパックが、乗物用シートのシートバック部に装着されている請求項１２記載の生体状態推定システム。
【請求項14】
前記表示手段が車載用モニタである請求項１０〜１３のいずれか１に記載の生体状態推定システム。
【請求項15】
生体信号測定手段により人の上体から採取した生体信号を用いて、人の状態を推定する生体状態推定装置に、
前記生体信号を分析して所定時点における人の恒常性維持機能レベルを複数の段階に区分して求める恒常性維持機能レベル演算ステップと、
前記恒常性維持機能レベルの各段階を縦軸に時間を横軸にとり、前記恒常性維持機能レベル演算ステップにより求められる恒常性維持機能レベルを時系列にプロットして表示手段にグラフとして表示させる出力ステップと
を実行させることを特徴とするコンピュータプログラム。
前記生体信号を分析して所定時点における人の恒常性維持機能レベルを複数の段階に区分して求める恒常性維持機能レベル演算ステップと、
前記恒常性維持機能レベルの各段階を縦軸に時間を横軸にとり、前記恒常性維持機能レベル演算ステップにより求められる恒常性維持機能レベルを時系列にプロットして表示手段にグラフとして表示させる出力ステップと
を実行させることを特徴とするコンピュータプログラム。
【請求項16】
前記出力ステップにより表示手段に表示されるグラフが折れ線グラフである請求項１５記載のコンピュータプログラム。
【請求項17】
前記恒常性維持機能レベル演算ステップは、恒常性維持機能レベルを高活性状態から機能低下状態まで５段階に区分して求め、前記出力ステップは、高活性状態を最上部とし、機能低下状態を最下部として縦軸に５段階に区分して表示させる請求項１５又は１６記載のコンピュータプログラム。
【請求項18】
前記出力ステップにより表示手段に表示されるグラフから、人の状態の変遷を判定する判定ステップをさらに実行させる請求項１５〜１７のいずれか１に記載のコンピュータプログラム。
【請求項19】
前記判定ステップにおいて、前記出力ステップにより表示手段に表示されるグラフに、恒常性維持機能レベルの２段階以上の低下が所定回数以上生じた場合に異常状態と判定する異常状態判定ステップを実行させる請求項１８記載のコンピュータプログラム。
【請求項20】
前記恒常性維持機能レベル演算ステップとして、
前記生体信号測定手段により得られる生体信号の周波数を求める周波数演算ステップと、
前記周波数演算ステップにより得られた前記生体信号の周波数から、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列波形を求める周波数傾き時系列解析演算ステップと、
前記周波数傾き時系列解析演算ステップにより求めた前記周波数傾きの時系列波形を微分する微分ステップと、
前記周波数傾き時系列解析演算ステップにより求めた前記周波数傾きの時系列波形を積分する積分ステップと、
前記周波数傾き時系列解析演算ステップにより求めた前記周波数傾きの時系列波形の増減から方形波を求める方形波算出ステップと、
任意に設定した第１時間帯における前記周波数傾き時系列波形と前記第１時間帯よりも後の第２時間帯における前記周波数傾き時系列波形と間の記述関数及び記述関数振幅値を求める記述関数算出ステップと、
前記周波数演算ステップにより得られる、前記生体信号の時系列波形の極大値を用いた生体信号の周波数の時系列波形と、前記生体信号の時系列波形の正負が切り替わるゼロクロス地点を用いた生体信号の周波数の時系列波形とを用いて、前記周波数傾き時系列解析演算ステップにより求めた、それぞれの周波数傾きの時系列波形を絶対値処理する絶対値処理ステップと、
前記周波数傾き時系列解析演算ステップにより求めた前記周波数傾き、前記微分手段により求めた微分値、前記積分手段により求めた積分値、前記方形波算出手段により求めた方形波の正負、前記記述関数算出ステップにより求めた記述関数振幅値、及び前記絶対値処理ステップにより求めた周波数傾き時系列波形の２つの絶対値のうち、少なくとも１つ以上を用いて恒常性維持機能レベルの段階を求める恒常性維持機能段階算出ステップと
を実行させる請求項１５〜１９のいずれか１に記載のコンピュータプログラム。
前記生体信号測定手段により得られる生体信号の周波数を求める周波数演算ステップと、
前記周波数演算ステップにより得られた前記生体信号の周波数から、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列波形を求める周波数傾き時系列解析演算ステップと、
前記周波数傾き時系列解析演算ステップにより求めた前記周波数傾きの時系列波形を微分する微分ステップと、
前記周波数傾き時系列解析演算ステップにより求めた前記周波数傾きの時系列波形を積分する積分ステップと、
前記周波数傾き時系列解析演算ステップにより求めた前記周波数傾きの時系列波形の増減から方形波を求める方形波算出ステップと、
任意に設定した第１時間帯における前記周波数傾き時系列波形と前記第１時間帯よりも後の第２時間帯における前記周波数傾き時系列波形と間の記述関数及び記述関数振幅値を求める記述関数算出ステップと、
前記周波数演算ステップにより得られる、前記生体信号の時系列波形の極大値を用いた生体信号の周波数の時系列波形と、前記生体信号の時系列波形の正負が切り替わるゼロクロス地点を用いた生体信号の周波数の時系列波形とを用いて、前記周波数傾き時系列解析演算ステップにより求めた、それぞれの周波数傾きの時系列波形を絶対値処理する絶対値処理ステップと、
前記周波数傾き時系列解析演算ステップにより求めた前記周波数傾き、前記微分手段により求めた微分値、前記積分手段により求めた積分値、前記方形波算出手段により求めた方形波の正負、前記記述関数算出ステップにより求めた記述関数振幅値、及び前記絶対値処理ステップにより求めた周波数傾き時系列波形の２つの絶対値のうち、少なくとも１つ以上を用いて恒常性維持機能レベルの段階を求める恒常性維持機能段階算出ステップと
を実行させる請求項１５〜１９のいずれか１に記載のコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、人の上体から得られる生体信号の時系列波形を用いて生体の状態を推定する技術に関する。【背景技術】
【0002】
運転中の運転者の生体状態を監視することは、近年、事故予防策等として注目されている。本出願人は、特許文献１〜３において、シートクッション部に圧力センサを配置し、臀部脈波を採取して分析し、入眠予兆現象を判定する手法を開示している。【0003】
具体的には、脈波の時系列波形を、それぞれ、ＳａｖｉｔｚｋｙとＧｏｌａｙによる平滑化微分法により、極大値と極小値を求める。そして、５秒ごとに極大値と極小値を切り分け、それぞれの平均値を求める。求めた極大値と極小値のそれぞれの平均値の差の二乗をパワー値とし、このパワー値を５秒ごとにプロットし、パワー値の時系列波形を作る。この時系列波形からパワー値の大域的な変化を読み取るために、ある時間窓Ｔｗ（１８０秒）について最小二乗法でパワー値の傾きを求める。次に、オーバーラップ時間Ｔｌ（１６２秒）で次の時間窓Ｔｗを同様に計算して結果をプロットする。この計算（移動計算）を順次繰り返してパワー値の傾きの時系列波形を得る。一方、脈波の時系列波形をカオス解析して最大リアプノフ指数を求め、上記と同様に、平滑化微分によって極大値を求め、移動計算することにより最大リアプノフ指数の傾きの時系列波形を得る。【0004】
そして、パワー値の傾きの時系列波形と最大リアプノフ指数の傾きの時系列波形が逆位相となっており、さらには、パワー値の傾きの時系列波形で低周波、大振幅の波形が生じている波形を、入眠予兆を示す特徴的な信号と判定し、その後に振幅が小さくなったポイントを入眠点と判定している。【0005】
また、特許文献４として、内部に三次元立体編物を挿入した空気袋（エアパック）を備え、このエアパックを人の背部に対応する部位に配置し、エアパックの空気圧変動を測定し、得られた空気圧変動の時系列データから人の生体信号を検出し、人の生体の状態を分析するシステムを開示している。また、非特許文献１及び２においても、腰腸肋筋に沿うようにエアパックセンサを配置して人の生体信号を検出する試みを報告している。このエアパックの空気圧変動は、心臓の動きに伴う下行大動脈の揺れによるものであり、特許文献１及び２の臀部脈波を利用する場合よりも、心臓の動きにより近い状態変化を捉えることができる。【0006】
また、出願人は、特許文献５として、パワー値の傾きの時系列信号を絶対値処理して積分値を算出し、その積分値を疲労度として求める技術を提案している。さらに、出願人は、特許文献６として、パワー値の時系列データを用いて主として肉体的な疲労度を求める手法にプラスして、最大リアプノフ指数の時系列データを用いた交感神経活動による疲労の代償作用を考慮することで、精神的疲労も加味した、より官能評価に近い疲労度を求める技術を提案している。【先行技術文献】
【特許文献】
【0007】
【特許文献1】特開２００４−３４４６１２号公報
【特許文献2】特開２００４−３４４６１３号公報
【特許文献3】ＷＯ２００５／０９２１９３Ａ１公報
【特許文献4】特開２００７−９００３２号公報
【特許文献5】ＷＯ２００５／０３９４１５Ａ１公報
【特許文献6】特開２００９−２９５３８４号公報
【非特許文献】
【0008】
【非特許文献1】「非侵襲型センサによって測定された生体ゆらぎ信号の疲労と入眠予知への応用」、落合直輝（外６名）、第３９回日本人間工学会 中国・四国支部大会 講演論文集、平成１８年１１月２５日発行、発行所：日本人間工学会 中国・四国支部事務局
【非特許文献2】「非侵襲生体信号センシング機能を有する車両用シートの試作」、前田慎一郎（外４名）、第３９回日本人間工学会 中国・四国支部大会 講演論文集、平成１８年１１月２５日発行、発行所：日本人間工学会 中国・四国支部事務局
【発明の概要】
【発明が解決しようとする課題】
【0009】
特許文献１〜４及び非特許文献１〜２の技術は、上記したように、パワー値の傾きの時系列波形と最大リアプノフ指数の傾きの時系列波形が逆位相となり、かつ、パワー値の傾きの時系列波形で低周波、大振幅の波形が生じた時点をもって入眠予兆現象と捉えている。この手法によれば、人がシートに着席すれば、入眠予兆現象を捉えることができる。また、特許文献５，６によれば、人の疲労度を求めることができる。【0010】
しかし、入眠予兆現象が生じたか否かを波形の変化から判断するには熟練を要し、一般のドライバが波形変化を見て入眠予兆を判断することは困難である。このため、上記特許文献では、時系列波形の逆位相等が生じた時点で、何らかの警告装置（例えば、ブザー音、シートバックの振動など）を作動させて覚知させる手段を採用している。また、特許文献５，６の疲労度は、時間と共に疲労度が右肩上がりに増加していくグラフで判断するものである。従って、この場合も、例えば、このグラフを自動車に搭載したモニタに表示させたとして、各ドライバがそのグラフを見て自らの疲労の程度を把握するのは容易ではない。このグラフによる判定に慣れた者が見て判断すれば別であるが、そうでなければ、疲労度がある数値以上に至った時点等、何らかの変曲点において、音等の警告によりドライバに覚知させるしかない。【0011】
また、従来の手法は、いずれも人の代謝変化に起因する状態を主として判定するものであった。つまり、時間経過と共に一方的に肉体が疲れていく過程を主として見ているものであった。しかし、緊張状態では、精神疲労の寄与率が高く、リラックス状態では肉体疲労が疲労の進行を支配すると考えられる。ところが、従来の手法では、脳機能の活性化による疲労の抑制が考慮されていない。一般的な運転環境下では、疲労の進行度合いに拘わらず、脳機能の活性化が生じる。この脳機能の活性化には創造性や好奇心などがあり、これらは覚醒誘導をもたらす。特許文献６のように、肉体的な疲労度を求める手法にプラスして、最大リアプノフ指数の時系列データを用いた交感神経活動による疲労の代償作用を考慮し、精神的疲労も加味した疲労度を求めているが、この場合でも、あくまで、一方的に増加していく疲労の過程、交感神経の代償作用を出力しているに過ぎず、脳機能が活性化しているか否かまでは捉えていない。【0012】
本発明は上記に鑑みなされたものであり、自律神経の活動だけでなく脳機能の活性化状態も捉えることができると共に、表示手段を見るだけで人の恒常性維持機能レベルを視覚的に容易に把握でき、それにより疲労の進行度合いも把握できる生体状態推定装置、生体状態推定システム及びコンピュータプログラムを提供することを課題とする。【課題を解決するための手段】
【0013】
上記課題を解決するため、本発明の生体状態推定装置は、生体信号測定手段により人の上体から採取した生体信号を用いて、人の状態を推定する生体状態推定装置であって、前記生体信号を分析して所定時点における人の恒常性維持機能レベルを複数の段階に区分して求める恒常性維持機能レベル演算手段と、前記恒常性維持機能レベルの各段階を縦軸に時間を横軸にとり、前記恒常性維持機能レベル演算手段により求められる恒常性維持機能レベルを時系列にプロットし、恒常性維持機能レベルのゆらぎ度合いを表示手段にグラフとして表示させる出力手段とを具備することを特徴とする。【0014】
前記恒常性維持機能レベル演算手段は、恒常性維持機能レベルを高活性状態から機能低下状態まで３〜１０段階の範囲で区分して求め、前記出力手段は、高活性状態を最上部とし、機能低下状態を最下部として縦軸に３〜１０段階の範囲に区分して表示させることが好ましい。前記恒常性維持機能レベル演算手段は、恒常性維持機能レベルを高活性状態から機能低下状態まで５段階に区分して求め、前記出力手段は、高活性状態を最上部とし、機能低下状態を最下部として縦軸に５段階に区分して表示させることが好ましい。【0015】
前記出力手段により表示手段に表示されるグラフから、人の状態の変遷を判定する判定手段をさらに有することが好ましい。前記判定手段は、前記出力手段により表示手段に表示されるグラフに、恒常性維持機能レベルの２段階以上の急低下が所定回数以上生じた場合に異常状態と判定する異常状態判定手段を有することが好ましい。【0016】
前記恒常性維持機能レベル演算手段は、前記生体信号測定手段により得られる生体信号の周波数を求める周波数演算手段と、前記周波数演算手段により得られた前記生体信号の周波数から、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列波形を求める周波数傾き時系列解析演算手段と、前記周波数傾き時系列解析演算手段により求めた前記周波数傾きの時系列波形を微分する微分手段と、前記周波数傾き時系列解析演算手段により求めた前記周波数傾きの時系列波形を積分する積分手段と、前記周波数傾き時系列解析演算手段により求めた前記周波数傾きの時系列波形の増減から方形波を求める方形波算出手段と、任意に設定した第１時間帯における前記周波数傾き時系列波形と前記第１時間帯よりも後の第２時間帯における前記周波数傾き時系列波形と間の記述関数及び記述関数振幅値を求める記述関数算出手段と、前記周波数演算手段により得られる、前記生体信号の時系列波形の極大値を用いた生体信号の周波数の時系列波形と、前記生体信号の時系列波形の正負が切り替わるゼロクロス地点を用いた生体信号の周波数の時系列波形とを用いて、前記周波数傾き時系列解析演算手段により求めた、それぞれの周波数傾きの時系列波形を絶対値処理する絶対値処理手段と、前記周波数傾き時系列解析演算手段により求めた前記周波数傾き、前記微分手段により求めた微分値、前記積分手段により求めた積分値、前記方形波算出手段により求めた方形波の正負、前記記述関数算出手段により求めた記述関数振幅値、及び前記絶対値処理手段により求めた周波数傾き時系列波形の２つの絶対値のうち、少なくとも１つ以上を用いて恒常性維持機能レベルの段階を求める恒常性維持機能段階算出手段とを具備してなることが好ましい。【0017】
前記恒常性維持機能レベル演算手段により処理される前に、前記生体信号測定手段により得られる生体信号を、所定の周波数帯域でフィルタリングするフィルタリング手段をさらに有し、前記フィルタリング手段が、前記生体信号測定手段により測定した静的環境下における生体信号の平均周波数を求め、該平均周波数を基準値として、フィルタリングする周波数帯域を設定する構成であることが好ましい。【0018】
前記生体信号測定手段は、人の背部に当接され、心臓の動きによる大動脈の揺れに伴って圧力が変動するエアパックを備えてなり、前記恒常性維持機能レベル演算手段は、前記エアパックの圧力変動を分析対象の前記生体信号として用いることが好ましい。乗物に適用する場合、前記エアパックが、乗物用シートのシートバック部に装着されていることが好ましく、前記表示手段として車載用モニタを用いることができる。【0019】
また、本発明の生体状態推定システムは、人の上体から生体信号を採取する生体信号測定手段と、前記生体信号測定手段により採取した生体信号を用いて、人の状態を推定する生体状態推定装置とを備えた生体状態推定システムであって、前記生体状態推定装置は、前記生体信号を分析して所定時点における人の恒常性維持機能レベルを複数の段階に区分して求める恒常性維持機能レベル演算手段と、前記恒常性維持機能レベルの各段階を縦軸に時間を横軸にとり、前記恒常性維持機能レベル演算手段により求められる恒常性維持機能レベルを時系列にプロットし、恒常性維持機能レベルのゆらぎ度合いを表示手段にグラフとして表示させる出力手段とを具備することを特徴とする。【0020】
前記生体信号測定手段は、人の背部に当接され、心臓の動きによる大動脈の揺れに伴って圧力が変動するエアパックを備えてなり、前記恒常性維持機能レベル演算手段は、前記エアパックの圧力変動を分析対象の前記生体信号として用いることが好ましい。乗物に適用する場合、前記エアパックが、乗物用シートのシートバック部に装着されていることが好ましく、前記表示手段として車載用モニタを用いることができる。【0021】
また、本発明のコンピュータプログラムは、生体信号測定手段により人の上体から採取した生体信号を用いて、人の状態を推定する生体状態推定装置に、前記生体信号を分析して所定時点における人の恒常性維持機能レベルを複数の段階に区分して求める恒常性維持機能レベル演算ステップと、前記恒常性維持機能レベルの各段階を縦軸に時間を横軸にとり、前記恒常性維持機能レベル演算ステップにより求められる恒常性維持機能レベルを時系列にプロットして表示手段にグラフとして表示させる出力ステップとを実行させることを特徴とする。【0022】
前記出力ステップにより表示手段に表示されるグラフが折れ線グラフであることが好ましい。前記恒常性維持機能レベル演算ステップは、恒常性維持機能レベルを高活性状態から機能低下状態まで５段階に区分して求め、前記出力ステップは、高活性状態を最上部とし、機能低下状態を最下部として縦軸に５段階に区分して表示させることが好ましい。【0023】
前記出力ステップにより表示手段に表示されるグラフから、人の状態の変遷を判定する判定ステップをさらに実行させることが好ましい。前記判定ステップにおいて、前記出力ステップにより表示手段に表示されるグラフに、恒常性維持機能レベルの２段階以上の低下が所定回数以上生じた場合に異常状態と判定する異常状態判定ステップを実行させることが好ましい。【0024】
前記恒常性維持機能レベル演算ステップとして、前記生体信号測定手段により得られる生体信号の周波数を求める周波数演算ステップと、前記周波数演算ステップにより得られた前記生体信号の周波数から、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列波形を求める周波数傾き時系列解析演算ステップと、前記周波数傾き時系列解析演算ステップにより求めた前記周波数傾きの時系列波形を微分する微分ステップと、前記周波数傾き時系列解析演算ステップにより求めた前記周波数傾きの時系列波形を積分する積分ステップと、前記周波数傾き時系列解析演算ステップにより求めた前記周波数傾きの時系列波形の増減から方形波を求める方形波算出ステップと、任意に設定した第１時間帯における前記周波数傾き時系列波形と前記第１時間帯よりも後の第２時間帯における前記周波数傾き時系列波形と間の記述関数及び記述関数振幅値を求める記述関数算出ステップと、前記周波数演算ステップにより得られる、前記生体信号の時系列波形の極大値を用いた生体信号の周波数の時系列波形と、前記生体信号の時系列波形の正負が切り替わるゼロクロス地点を用いた生体信号の周波数の時系列波形とを用いて、前記周波数傾き時系列解析演算ステップにより求めた、それぞれの周波数傾きの時系列波形を絶対値処理する絶対値処理ステップと、前記周波数傾き時系列解析演算ステップにより求めた前記周波数傾き、前記微分手段により求めた微分値、前記積分手段により求めた積分値、前記方形波算出手段により求めた方形波の正負、前記記述関数算出ステップにより求めた記述関数振幅値、及び前記絶対値処理ステップにより求めた周波数傾き時系列波形の２つの絶対値のうち、少なくとも１つ以上を用いて恒常性維持機能レベルの段階を求める恒常性維持機能段階算出ステップとを実行させることが好ましい。
【発明の効果】
【0025】
本発明は、人の上体から採取した生体信号を用いて、恒常性維持機能レベルを複数の段階に区分して求め、恒常性維持機能レベルの各段階を時間軸に従って時系列にプロットして表示手段にグラフとして表示させる構成であり、好ましくは、グラフの縦軸の最上部を高活性状態として最下部を機能低下状態として表示する。これにより、自律神経による恒常性を維持しようとするゆらぎの様子と、脳による刺激が加えら得た様子を周期関数として捉えることができる。【0026】
また、このように恒常性維持機能レベルをグラフで時系列に表示した場合、人はそれを視認すると、恒常性維持機能の劣化状態（疲労の進行過程）を熟練を要することなく容易に把握することができる。恒常性維持機能の変化であるこのグラフを見るという行為は、自分では主観的及び客観的に把握できない状態を知りたいという高次の欲求を満たす行為であるため、それにより脳の活性化が生じ、入眠に入りにくくすることができる。いわば、このグラフを見るそのこと自体が、覚醒誘導をもたらすものである。従って、本発明は脳の活性化を促す効果が高く、従来のブザー音などにより聴覚を刺激する方法や、シートバックの振動などのような触覚・圧覚を刺激する方法等と比較して高い覚醒誘導効果が期待できる。【0027】
また、本発明において恒常性維持機能レベルを求める生体信号は、人の上体の生体反応を利用して検出されるものであるため、心臓の動きに伴う大動脈の揺れを捕捉することができる。これにより、生命維持機能である恒常性維持機能を見ることができ、外乱に対する適応力を見ている指尖容積脈波を用いて人の状態を推定する技術と比較し、より正確に人の状態を把握することができる。【0028】
人の状態変化のサイクルには、覚醒と睡眠の各状態とその二つの状態への遷移状態に大別される。人は、高活性状態から睡眠への遷移状態いわゆる疲労代償期を経由して、疲労している状態すなわち機能低下状態に遷移する。その後、回復期の遷移状態を経由して高活性状態に戻る。これらの遷移状態において入眠予兆信号、覚醒予兆信号が発現する。また、人体は、恒常性を維持するために、血液の循環を利用して代謝で使われるエネルギーの変化を最小限にする。人体は血液の循環動態を変えるために、心拍の周波数（心拍数）と血管径を変える。従って、心拍の周波数変動と血管径の変化を検出することで、代謝で使われるエネルギーの遷移を同定できる。この遷移の様子が捉えられると、高活性状態から機能低下状態への遷移が捉えられ、入眠への移行を予測できる。この状態変化は、自律神経系によって引き起こされる。この変化を末梢系の反応を通して捉えると外的刺激への適応状態が分かるが、中枢系の反応すなわち心臓の動きに伴う大動脈の揺れを捕捉すれば疲労の質と量の同定が可能となる。【0029】
また、グラフの進行過程の中で、恒常性維持機能レベルの２段階以上の低下が生じるということは、エネルギーレベルの急激な変化を示すものである。本発明のようなグラフを用いることにより、そのような急低下を視認できるため、人は急激な体調変化を意識し、自ら高活性状態に誘導しようと意識したり、積極的に休息をとることを意識したりすることができる。また、上記のような急激な変化が所定回数以上生じた場合には異常状態（入眠予兆信号の出現）であると自動的に判定する手段を設けることがより好ましい。異常状態と判定された場合には、警告音を鳴らすなど、聴覚や圧覚の刺激等と併用することができる。【図面の簡単な説明】
【0030】
【図9】図９は、生体信号測定手段により検出した脈波（心拍変動）のピーク値又はゼロクロス地点を用いて、周波数変動時系列波形、周波数変動時系列波形の基線、周波数変動の傾き時系列である周波数傾き時系列波形、及び積分曲線を求める方法を説明するための図である。
【図10】図１０は上記実施形態に係る生体状態推定装置において用いたフィルタリング手段を施した波形と施さない波形とを示したものであり、（ａ）は自動車のアイドリング状態で行った測定結果を、（ｂ）は走行状態で行った測定結果を示した図である。
【発明を実施するための形態】
【0031】
以下、図面に示した本発明の実施形態に基づき、本発明をさらに詳細に説明する。図１は、本実施形態に係る生体状態推定装置６０の分析対象である生体信号、ここでは心臓の動きに伴う背部の大動脈の揺れによる圧力変動を採取する生体信号測定手段１を組み込んだ自動車用のシート５００の外観を示した図である。なお、生体信号測定手段１と生体状態推定装置６０とにより、生体状態推定システムが構成される。生体信号測定手段１は、シートバック部５１０に組み込まれて用いられる。ここで、生体信号測定手段１によって採取される圧力変動の検知信号には、ノイズ信号がより少ないことが望ましい。そこで、本実施形態の生体信号測定手段１は、以下に説明するように、自動車の走行中等の振動環境下においても、センサの出力信号自体に含まれるノイズ信号を少なくできるような工夫がなされている。【0032】
生体信号測定手段１は、エアパックユニット１００と、ビーズ発泡樹脂弾性部材２０とを有して構成されている。エアパックユニット１００は、収容体１５と、該収容体１５に収容した２つのエアパック１０を備えて構成される。各エアパック１０は、３つの小空気袋１１１が縦方向に連接されている一方、そのそれぞれは空気の流通がないように形成されている。各小空気袋１１１内には、復元力付与部材としての三次元立体編物１１２が配置されている。【0033】
本実施形態では、このようなエアパック１０が左右に配置される。左右に配置することにより、着座者の背への当たりが左右均等になり、違和感を感じにくくなる。また、左右のエアパック１０，１０のいずれか一方を構成するいずれかの小空気袋１１１にセンサ取付チューブ１１１ａが設けられ、その内側に空気圧変動を測定するセンサ１１１ｂが固定されている。なお、センサ取付チューブ１１１ａは密閉されている。小空気袋１１１は、このような生体信号による空気圧変動に敏感に反応させるために、大きさは、幅４０〜１００ｍｍ、長さ１２０〜２００ｍｍの範囲が好ましい。小空気袋１１１の素材は限定されるものではないが、例えば、ポリウレタンエラストマー（例えば、シーダム株式会社製、品番「ＤＵＳ６０５−ＣＤＲ」）からなるシートを用いて形成することができる。センサ１１１ｂとしては、小空気袋１１１内の空気圧を測定できるものであればよく、例えば、コンデンサ型マイクロフォンセンサを用いることができる。【0034】
小空気袋１１１を３つ連接した場合の全体の大きさとしては、自動車のシート５００のシートバック部５１０に用いる場合、幅４０〜１００ｍｍ、全長４００〜６００ｍｍの範囲とすることが好ましい。長さが短い場合、シートバック部５１０において、着座者が、腰部付近の一部分のみに異物感を感じるため、４００ｍｍ以上の長さとして、できるだけ、着座者の背全体に対応させることが好ましい。【0035】
空気圧変動を検出するセンサ１１１ｂは、本実施形態では、着座者の左側に配置されるエアパック１０を構成する中央の小空気袋１１１に設けている。この小空気袋１１１の位置は、着座者の背部の大動脈（特に、「下行大動脈」）の心臓の動きに伴って生じる揺れを検知可能な領域に相当する。このような背部の大動脈の動きを検知可能な領域は、着座者の体格により一律ではないが、身長１５８ｃｍの日本人女性から身長１８５ｃｍの日本人男性までの様々な体格の被験者２０名で測定したところ、該小空気袋１１１（幅６０ｍｍ、長さ１６０ｍｍ）をシートバック部５１０の中心寄りの側縁と下縁の交差部Ｐ(図２及び図３参照）が、シートクッション部５２０の上面からシートバック部５１０の表面に沿った長さＬ：２２０ｍｍ、シートバック部５１０の中心からの距離Ｍ：８０ｍｍとなるように設定したところ、上記全ての被験者において大動脈の動きによる圧力変動を検知できた。小空気袋１１１の大きさが、幅４０〜１００ｍｍ、長さ１２０〜２００ｍｍの範囲の場合、交差部Ｐの位置を、シートクッション部５２０の上面からシートバック部５１０の表面に沿った長さで１５０〜２８０ｍｍ、シートバック部５１０の中心から６０〜１２０ｍｍの範囲に設定することが好ましい。【0036】
上記した２つのエアパック１０をシートバック部５１０において容易に所定の位置に設定できるようにユニット化しておくことが好ましい。従って、図２〜図４に示したような収容体１５にエアパック１０を装着したエアパックユニット１００として構成とすることが好ましい。収容体１５は、両側にエアパック１０を収容する袋状のエアパック収容部１５１を有し、２つのエアパック収容部１５１間に接続部１５２を有している。【0037】
２つのエアパック収容部１５１には、それぞれエアパック１０が挿入される。また、エアパック収容部１５１には、エアパック１０とほぼ同じ大きさの三次元立体編物４０を、エアパック１０の裏側エアパック１２の背面側に重ねて挿入することが好ましい（図３（ｄ）参照）。三次元立体編物４０を配置することにより、エアパック１０が該三次元立体編物４０によっていわば浮くように支持されるため、シートバック部５１０からの外部振動がエアパック内部の圧力変動として伝わりにくくなる。すなわち、三次元立体編物４０を配置することにより、高周波小振幅の外部振動が入力された場合には、三次元立体編物４０のパイルと空気圧の変動から、エアパック内にバネ定数の低い、バネ・マス・ダンパ系が作られる。そして、それが三次元立体編物４０を内蔵したエアパック１０において、低・高周波入力に対するフィルタ（ローパスフィルタ・ハイパスフィルタ）として作用し、該外部振動を減衰する。【0038】
接続部１５２は、２つのエアパック部１５１を所定間隔をおいて支持できるものであればよく、幅６０〜１２０ｍｍ程度で形成される。接続部１５２も、袋状に形成し、その内部に三次元立体編物４５を挿入することが好ましい(図３（ｄ）及び図４参照）。これにより、該接続部１５２を通じて入力される振動も、該三次元立体編物４５を挿入することにより効果的に除振でき、センサ１１１ｂを備えたエアパック１０への外部振動の伝達を抑制できる。【0039】
なお、上記したように、小空気袋１１１は、例えば、ポリウレタンエラストマー（例えば、シーダム株式会社製、品番「ＤＵＳ６０５−ＣＤＲ」）からなるシートを用いて形成されるが、収容体１５も同じ素材を用いて形成することが好ましい。また、小空気袋１１１、エアパック収容部１５１及び接続部１５２内に装填される各三次元立体編物は、例えば、特開２００２−３３１６０３号公報に開示されているように、互いに離間して配置された一対のグランド編地と、該一対のグランド編地間を往復して両者を結合する多数の連結糸とを有する立体的な三次元構造となった編地である。【0040】
一方のグランド編地は、例えば、単繊維を撚った糸から、ウェール方向及びコース方向のいずれの方向にも連続したフラットな編地組織（細目）によって形成され、他方のグランド編地は、例えば、短繊維を撚った糸から、ハニカム状（六角形）のメッシュを有する編み目構造に形成されている。もちろん、この編地組織は任意であり、細目組織やハニカム状以外の編地組織を採用することもできるし、両者とも細目組織を採用するなど、その組み合わせも任意である。連結糸は、一方のグランド編地と他方のグランド編地とが所定の間隔を保持するように、２つのグランド編地間に編み込んだものである。【0041】
ビーズ発泡樹脂弾性部材２０は、シートバック部５１０の表皮部材とエアパック１０を収容した収容体１５（エアパックユニット１００）との間に配設され、２つのエアパック１０の全長に相当する長さを有し、２つのエアパック１０の頂部間の長さに相当する幅を有している。従って、長さが４００〜６００ｍｍ、幅が２５０〜３５０ｍｍ程度の大きさのものを用いることが好ましい。これにより、２つのエアパック１０が共に覆われるため、２つのエアパック１０の凹凸を感じにくくなる。【0042】
ビーズ発泡樹脂弾性部材２０は、平板状に形成されたビーズ発泡体と、その外面に貼着される被覆材とから構成されている。ビーズ発泡体としては、ポリスチレン、ポリプロピレン及びポリエチレンのいずれか少なくとも一つを含む樹脂のビーズ法による発泡成形体が用いられる。なお、発泡倍率は任意であり限定されるものではない。被覆材は、ビーズ発泡体の外面に接着により貼着され、高い伸度と回復率を有する素材であり、好ましくは、伸度２００％以上、１００％伸長時の回復率が８０％以上である弾性繊維不織布が用いられる。例えば、特開２００７−９２２１７号公報に開示された熱可塑性エラストマー弾性繊維が相互に溶融接着された不織布を用いることができる。具体的には、ＫＢセーレン（株）製、商品名「エスパンシオーネ」を用いることができる。【0043】
また、ビーズ発泡樹脂弾性部材２０の外面を覆う被覆材として、例えば、熱可塑性ポリエステルからなる不織布を用いることもできる。具体的には、帝人（株）製のポリエチレンナフタレート（ＰＥＮ）繊維（１１００ｄｔｅｘ）から形成した２軸織物（縦：２０本／ｉｎｃｈ、横：２０本／ｉｎｃｈ）を用いることができる。【0044】
ビーズ発泡樹脂弾性部材２０を構成するビーズ発泡体としては、厚さ約５〜６ｍｍ程度のものを用いることができ、その外面に、厚さ約１ｍｍ以下の上記した弾性繊維不織布や熱可塑性ポリエステルからなる不織布を貼着して形成される。なお、本実施形態では、ビーズ発泡樹脂弾性部材２０の表皮部材５１１に対向する面とその反対面に、上記のエスパンシオーネ（商品名）を貼着している。これにより、生体信号の伝達性が向上する。【0045】
本実施形態において人体支持手段を構成するシート５００のシートバック部５１０は、表皮部材５１１と該表皮部材５１１の背面側に配設されるクッション支持部材５１２とを備えてなり、該表皮部材５１１とクッション支持部材５１２との間にエアパック１０を保持した収容体１５（エアパックユニット１００）とビーズ発泡樹脂弾性部材２０が組み込まれる。この際、クッション支持部材５１２側にまずエアパック１０を保持した収容体１５（エアパックユニット１００）が配置され、その表面側にビーズ発泡樹脂弾性部材２０が配置された上で、表皮部材５１１により被覆される。なお、クッション支持部材５１２は、例えば、三次元立体編物をシートバック部５１０の左右一対のサイドフレームの後端縁間に張って形成することもできるし、合成樹脂板から形成することもできる。表皮部材５１１は、例えば、三次元立体編物、合成皮革、皮革、あるいはこれらの積層体などを左右一対のサイドフレームの前縁間に張って設けることができる。【0046】
このように、本実施形態においては、表皮部材５１１の裏面側に所定の大きさのビーズ発泡樹脂弾性部材２０が積層して配置され、さらにその後方に左右一対のエアパック１０を保持した収容体１５（エアパックユニット１００）が配置される構成であるため、着座者が背にエアパック１０の凹凸を感じることがなくなり、生体信号を測定するためのエアパック１０を有する構成でありながら、座り心地が向上する。なお、上記した説明では、ビーズ発泡樹脂弾性部材２０を１枚用いたのみであるが、複数枚重ねて配置することも可能である。【0047】
なお、上記した生体信号測定手段１は、エアパックユニット１００をシート５００のシートバック部５１０に組み込んでいるが、シートバック部５１０の表面に後付で取り付けるシート用クッションに組み込むようにしてもよい。また、エアパック１０を形成する小空気袋１１１を３つ連接しているが、２つであってもよいし、１つであってもよい。また、背部の大動脈の揺動を捉えられる位置に配置して用いることができる限り、例えば、１つの空気袋からなるエアパック１０のみから生体信号測定手段１を構成することもできる。また、上記したサイズは小空気袋１１１を３つ連接して用いた場合の適切な例であり、エアパック１０のサイズ、エアパックユニット１００のサイズは、背部の大動脈の揺動を捉えられる限り全く限定されるものではないことはもちろんである。【0048】
次に、生体状態推定装置６０の構成について図６に基づいて説明する。生体状態推定装置６０には、ハードディスク等の記憶部に設定されるコンピュータプログラムからなる恒常性維持機能レベル演算手段（恒常性維持機能レベル演算ステップ）６１と、出力手段（出力ステップ）６２とを有して構成される。なお、コンピュータプログラムは、フレキシブルディスク、ハードディスク、ＣＤ−ＲＯＭ、ＭＯ（光磁気ディスク）、ＤＶＤ−ＲＯＭ、メモリカードなどの記録媒体へ記憶させて提供することもできるし、通信回線を通じて伝送することも可能である。【0049】
恒常性維持機能レベル演算手段（恒常性維持機能レベル演算ステップ）６１は、生体信号測定手段１により検出された背部の大動脈の動きによる圧力変動の信号（以下、場合により「エアパック信号」という）の時系列波形から人の状態を分析し、所定時点における人の恒常性維持機能レベルを複数の段階に区分して求める手段である。人の状態は、上記したように、大きく高活性状態と機能低下状態とに分けられると共に、両者間は、さらに定常状態が含まれ、これら各状態間に遷移状態が存在する。但し、本実施形態では、高活性状態から機能低下状態までを５段階の恒常性維持機能レベルとして区分した。つまり、覚醒時における最も元気な疲労を感じていない状態（高活性状態）から、疲労に陥った状態（機能低下状態）までの間を５段階の恒常性維持機能レベルに区分した。本実施形態では５段階に区分したが、これに限定されるものではない。但し、２段階では、異常状態（入眠予兆信号）を検知する恒常性維持機能レベルの２段階以上の急激な低下を示すことができず、１０段階を超えると細分過ぎて視覚的に把握しにくくなるため、３〜１０段階の範囲で区分することが好ましい。より好ましくは５段階である。なお、恒常性維持機能レベル演算手段６１による恒常性維持機能レベルの求め方の詳細については後述する。【0050】
出力手段（出力ステップ）６２は、恒常性維持機能レベルの各段階を縦軸として、時間を横軸としてなるプロットエリアを表示手段６５に表示する。そして、このプロットエリアに、恒常性維持機能レベル演算手段６１により求められる恒常性維持機能レベルを算出時間に対応して時系列にプロットし、折れ線グラフとして表示する。恒常性維持機能レベルを時系列にプロットしていくことにより、自律神経による恒常性を維持しようとするゆらぎと脳による賦活化の様子の周期関数が表される。【0051】
縦軸の恒常性維持機能レベルは、図１１に示したように、本実施形態では５段階に区分されて表示される。具体的には、高活性状態（覚醒時の中でも最も元気な疲労を感じていない状態）が最上部に表示され、機能低下状態が最下部に表示され、両者の間に３段階の恒常性維持機能段階が表示される。ここでは、中央部（上から３つ目）を定常状態とし、上から２つ目を高活性状態と定常状態との間の状態、下から２つ目を定常状態と機能低下状態と間の状態とし、そのイメージを人の顔の様子を模した図で表示している。なお、説明の便宜上、最も元気な状態である最上部を恒常性維持機能段階１（高活性状態（調子の良い状態：「イイ調子」等と表示））とし、以下、下方に向かって順に、恒常性維持機能段階２（比較的良い状態：「大丈夫」等と表示）、恒常性維持機能段階３（定常（普通）状態）、恒常性維持機能段階４（やや機能低下状態：「運転注意」等と表示）、恒常性維持機能段階５（機能低下状態：「疲労状態」等と表示））と表現する。【0052】
図１１（ａ）〜（ｅ）は、ドライバの恒常性維持機能レベルを求め、出力手段６２により表示手段６５に折れ線グラフで表示した例を示したものである。これは、複数の被験者の試験データをまとめてその傾向を示したグラフである。なお、表示手段６５としては、自動車の運転席から見やすい位置に搭載した車載用モニタを用いることが好ましい。【0053】
図１１（ａ）の例は、基本的には、最も元気な状態である恒常性維持機能段階１の範囲でグラフが描かれていく一方で、５〜８分間毎に、恒常性維持機能段階２に一時的に移行し、再び恒常性維持機能段階１に戻ることを繰り返している。仮に人が長時間に亘って恒常性維持機能段階１を維持し続けるとすれば、これは緊張感を伴って元気な状態を維持しているものであるが、緊張して集中力の高い状態、あるいはリラックスして集中力が高く元気で快適な場合には、自律神経の働きによる恒常性維持機能により恒常性維持機能段階１と恒常性維持機能段階２との間をこのように揺らぐ波形となる。この間で、ゆらぎの回数が多いときはリラックス状態を示し、ゆらぎの少ないときは緊張状態を示す。【0054】
図１１（ｂ）の例は、恒常性維持機能段階１から恒常性維持機能段階２に移行し、恒常性維持機能段階３、恒常性維持機能段階２、恒常性維持機能段階３、恒常性維持機能段階４、恒常性維持機能段階３、恒常性維持機能段階２というように各恒常性維持機能段階を順に行き来するゆらぎを伴いながら、徐々に疲労が増していく波形である。つまり、肉体疲労の進行により恒常性維持機能が順次低下し、やがては睡眠を伴う休息が必要となるケースである。いわば、極端な状態変化を伴わない通常の疲労の進行過程を示したものである。【0055】
図１１（ｃ）の例は、運転スタート時の恒常性維持機能段階１から、恒常性維持機能段階２、恒常性維持機能段階３、恒常性維持機能段階５と一気に疲労が進行したケースである。恒常性維持機能段階５まで進行したため、交感神経の働き（代償作用）によって強制的に恒常性維持機能段階３まで戻す機能が働くが、その後、恒常性維持機能段階５と恒常性維持機能段階３との間で極めて短時間で変化している。これは、眠気と安静状態との間の漫然とした状態が生じていることを示すものであり、ヒューマンエラーが生じやすく、運転に注意が必要な状態と判定できる。【0056】
図１１（ｄ）の例は、１５分過ぎに交感神経代償作用によって恒常性維持機能段階３から恒常性維持機能段階１に戻しているものの、３０分過ぎからはゆらぎがほとんど生じず恒常性維持機能段階１から恒常性維持機能段階５へとほぼ一方的に進行している。これは恒常性維持が難しくなっていることを示すものであり、早期の休憩が必要な状態と判定できる。【0057】
図１１（ｅ）の例は、高活性状態である恒常性維持機能段階１の範囲でグラフが描かれていくが、ほとんどゆらぎがない。これは、過緊張して集中力の高い状態が継続し、ゆらぎが少なくなり、エアパックがゆらぎを検出していない状態、つまり脳の賦活化によって過緊張している状態を示すものである。ところが、図１１（ｅ）のグラフ中、２５分付近、４０分付近、４７分付近において、数秒から数十秒で恒常性維持機能レベルが２段階以上（この例では、恒常性維持機能段階１から恒常性維持機能段階４までの３段階）一気に低下して一気に元に戻る尖鋭な突出ラインが生じている。この尖鋭な突出ラインは過緊張の継続に対する反動であり、過緊張状態の中での疲労の進行に通じる。この実験において、尖鋭な突出ラインが複数回生じた後に、ドライバにマイクロスリープが生じているが、これは生体の恒常性維持機能の著しい低下と考えられることから、この尖鋭な突出ラインはマイクロスリープに至る直前の信号である入眠予兆を示すものと定義した。なお、入眠予兆信号のより正確な判定については、この恒常性維持機能レベルの急低下だけでなく、後述するように、方形波の正負、ゼロクロス法により求めた絶対値の優位性を考慮することが好ましい。よって、このような突出ラインが複数回生じた場合には、それを入眠予兆信号と判定し、速やかな休憩が必要な状態と判定する。【0058】
このように、本実施形態によれば、恒常性維持機能レベルを視覚的に明確に捉えることができる。すなわち、ドライバは、例えば、図１１（ａ）のグラフを視覚的に捉えた場合には、順調に運転できる状態であることを理解し、図１１（ｃ）のグラフを視覚的に捉えた場合には、運転に注意しつつ比較的速やかに休憩した方がよい状態であることを理解し、図１１（ｅ）のような突出ラインを視認した場合には、入眠予兆信号の可能性があることから、比較的早期に休憩を取る方がよいことを理解できる。つまり、このような表示手段６５へのグラフによる表示は、短時間ではあるが、ドライバがそれを見て自らの状態を把握することで、ドライバの脳を刺激し賦活化が生じ、覚醒方向に誘導して、早めの休憩を促すことができるものである。【0059】
本実施形態では、ドライバが表示手段６５に表示されたグラフを見ることによって自らの状態がどのような状態であるかを容易に把握できるものであり、それにより脳への刺激を促すものであるが、休憩が必要な状態か、入眠予兆信号が出現したかなどを自動的に判定する判定手段（判定ステップ）６３を組み込むこともできる。【0060】
判定手段６３は、上記した恒常性維持機能段階のグラフにける基線の動きを監視し、例えば、図１１（ｃ）、（ｄ）のように、恒常性維持機能段階３〜５の範囲に基線が移動して、疲労の回復が見られないようになった状態に至ったならば、音、振動などによる警告を発するようにすることができる。本実施形態は、恒常性維持機能段階を視覚的に把握することで脳の活性化が生じ覚醒方向に誘導されるものであるが、判定手段６３によるいわば監視システムを含めることで、音、振動などの警告を併用することもできる。特に、図１１（ｅ）のような入眠予兆信号が出現した場合、例えば、３回出現した段階で、それを異常状態と判定し、より大きな警告音などを発する異常状態判定手段を設けた構成とすると、居眠り運転をより確実に抑制するのに役立つ。【0061】
ここで、上記出力手段６２は、人の恒常性維持機能レベルを時系列に折れ線グラフで表示させるものであるが、表示手段６５に、上記の折れ線グラフとは異なる表示を異なる画面で出力する別画面出力手段６２１を備える構成とすることができる。図１２（ａ）〜（ｅ）は、この別画面出力手段６２１により表示手段６５に出力された画像を示したものである。【0062】
別画面出力手段６２１により表示される画面には、図１３に拡大して示した代表例のように、この実施形態では、恒常性維持機能段階１〜５に対応したイメージを中性的な人の顔の様子（上記折れ線グラフの縦軸の顔のイメージ図と同じ）を示した第１画像部６２２と時系列のグラフからなる第２画像部６２３とが表示される。第１画像部６２２は、中心に顔イメージ表示部６２２ａと、その下部の心拍数表示部６２２ｂと、その上部の文字表示部６２２ｃと、右側に表示された入眠予兆信号表示部６２２ｄとを有している。顔イメージ表示部６２２ａには、恒常性維持機能レベル演算手段６１により求められる恒常性維持機能段階に応じて、別画面出力手段６２１が各恒常性維持機能段階に対応する顔のイメージを表示する。また、文字表示部６２２ｃには、別画面出力手段６２１が各恒常性維持機能段階に応じて設定した文字、例えば、「イイ調子」、「大丈夫」、「普通状態」、「運転注意」、「疲労状態」などを表示する。また、入眠予兆信号表示部６２２ｄは、５つの表示ランプを有しており、図１１（ｅ）に示したような入眠予兆信号を生じた場合に、その出現個数が各表示ランプに点灯するようになっている。これにより、人は、表示ランプの点灯を見て、休憩が必要なタイミングを自覚することができる。【0063】
第２画像部６２３には、その上部表示部６２３ａに、生体信号測定手段１により検出された背部の大動脈の動きによる圧力変動の信号（エアパック信号）の時系列波形のグラフが表示され、下部表示部６２３ｂに、周波数の傾きの時系列波形のグラフが表示される。【0064】
別画面出力手段６２１により表示される第１画像部６２２、第２画像部６２３の各画像部は、これらの表示に限定されるものではないが、特に、第１画像部６２２の顔イメージ表示部６２２ａ、文字表示部６２２ｃ及び入眠予兆信号表示部６２２ｄをドライバが視認することにより、現在の状況を容易に把握できる。但し、第１画像部６２２の表示は、ある瞬間における状態を示すだけであるため、画面を見て、自らの状態変化を経時的に把握し、将来の自分の状態変化を予測することで、休憩の有無を強く意識しやすいことから、運転中は、図１１に示したような折れ線グラフの画面をドライバが視認できるようにすることが好ましい。【0065】
なお、始動時において、図１１の折れ線グラフの画面を表示させるか、図１２の第１画像部６２２，第２画像部６２３の画面を表示させるかは任意に設定できるが、本実施形態では、初期画面として、図１２及び図１３に示した、第１画像部６２２，第２画像部６２３の画面を表示させ、「ＳＴＡＲＴ」ボタンを押すと、測定が開始され、恒常性維持機能レベル演算手段６１による演算が始まるようにしている。そして、「履歴」ボタンを押すと、図１１の折れ線グラフの画面が表示されるように設定している。すなわち、その瞬間における変化を図１２及び図１３に示した第１画像部６２２，第２画像部６２３からなる画面により表示し、過去の履歴を含んだ時系列の変化を図１１に示した折れ線グラフの画面により表示する。また、本実施形態では、運転中、図１１の折れ線グラフの画面が表示されるようにするために、「ＳＴＡＲＴ」、「履歴」の各ボタンを操作するようにすることで手指の動きによって脳の活性化が図られるという利点もある。また、本実施形態によれば、図１１の折れ線グラフの画面、図１２の第１画像部６２２（特に心拍数の部分）、第２画像部６２３の合計３つの画面を視認することになるが、３つの画面を見ることで、３つの画面から総合的に自らの状態を把握しようとするため、そのことも脳の活性化につながる。３つより少ない場合には脳の活性化効果が低く、３つより多い場合には煩雑になりすぎる。なお、「ＳＴＡＲＴ」、「履歴」の各ボタン操作は運転の支障にならないように、ハンドル付近に設けておくなどすることが好ましいことはもちろんである。【0066】
次に、恒常性維持機能レベル演算手段（恒常性維持機能レベル演算ステップ）６１に設定した具体的な演算手法について説明する。すなわち、恒常性維持機能レベル演算手段６１は、周波数演算手段（周波数演算ステップ）６１１、周波数傾き時系列解析演算手段（周波数傾き時系列解析演算ステップ）６１２、周波数変動時系列解析演算手段（周波数変動時系列解析演算ステップ）６１３、微分手段（微分ステップ）６１４、積分手段（積分ステップ）６１５、方形波算出手段（方形波算出ステップ）６１６、記述関数算出手段（記述関数算出ステップ）６１７、絶対値処理手段（絶対値処理ステップ）６１８、恒常性維持機能段階算出手段（恒常性維持機能段階算出ステップ）６１９とを備えてなる。【0067】
周波数演算手段６１１は、生体信号測定手段１から得られるエアパック信号の時系列データ（好ましくは後述のようにフィルタリング処理された所定の周波数領域のエアパック信号フィルタ波形）における周波数の時系列波形を求める。この時系列波形の求め方には２種類あり、第１の手法は、副交感神経機能の指標として用いられるＨＦに関するもので、第２の手法は、交感神経機能の指標として用いられるＬＦ／ＨＦ（なお、ＬＦは０．０５−０．２０Ｈｚの周波数成分、ＨＦは０．２０−０．３５Ｈｚの周波数成分である）に関するものである。【0068】
第１の手法は、エアパック信号の時系列波形を平滑化微分して極大値（ピーク）を用いて時系列波形を求める方法（以下、「ピーク検出法」という）である。ピーク検出法は、基本は、ＨＦの機能に相当する時系列波形であり、ＬＦ／ＨＦの代償機能に相当する波形が重畳されることもある時系列波形である。例えば、ＳａｖｉｔｚｋｙとＧｏｌａｙによる平滑化微分法により極大値を求める。次に、例えば５秒ごとに極大値を切り分け、その５秒間に含まれる時系列波形の極大値（波形の山側頂部）間の時間間隔の逆数を個別周波数ｆとして求め、その５秒間における個別周波数ｆの平均値を当該５秒間の周波数Ｆの値として採用する（図９の[1]のステップ）。そして、この５秒毎に得られる周波数Ｆをプロットすることにより、周波数の時系列波形を求める（図９の[2]のステップ）。【0069】
第２の手法は、エアパック信号の時系列波形において、正から負に切り替わる地点（以下、「ゼロクロス地点」という）を用いて時系列波形を求める方法（以下、「ゼロクロス法」という）である。このゼロクロス法は、脈波の周波数の基本成分を捉えるもので、ＬＦ／ＨＦに相当するものである。この方法では、まず、ゼロクロス地点を求めたならば、それを例えば５秒毎に切り分け、その５秒間に含まれる時系列波形のゼロクロス地点間の時間間隔の逆数を個別周波数ｆとして求め、その５秒間における個別周波数ｆの平均値を当該５秒間の周波数Ｆの値として採用する（図９の[1]のステップ）。そして、この５秒毎に得られる周波数Ｆをプロットすることにより、周波数の時系列波形を求める（図９の[2]のステップ）。【0070】
周波数傾き時系列解析演算手段（周波数傾き時系列解析演算ステップ）６１２は、周波数演算手段６１１によって、ピーク検出法又はゼロクロス法を用いて得られたエアパック信号の周波数の時系列波形から、所定の時間幅の時間窓を設定し、時間窓毎に最小二乗法により該エアパック信号の周波数の傾きを求め、その時系列波形を出力する構成である。周波数傾き時系列解析演算手段６１２により得られる周波数傾き時系列波形は、交感神経及び副交感神経の発現のバランスを捉えた生体のゆらぎを示すものとして出力される。具体的には、まず、ある時間窓Ｔｗ１における周波数の傾きを最小二乗法により求めてプロットする（図９の[3],[5]のステップ）。次に、オーバーラップ時間Ｔｌ（図９の[6]のステップ）で次の時間窓Ｔｗ２を設定し、この時間窓Ｔｗ２における周波数の傾きを同様に最小二乗法により求めてプロットする。この計算（移動計算）を順次繰り返し、エアパック信号の周波数の傾きの時系列変化を周波数傾き時系列波形として出力する（図９の[8]のステップ）。なお、時間窓Ｔｗの時間幅は１８０秒に設定することが好ましく、オーバーラップ時間Ｔｌは１６２秒に設定することが好ましい。これは、本出願人による上記特許文献３（ＷＯ２００５／０９２１９３Ａ１公報）において示したように、時間窓Ｔｗの時間幅及びオーバーラップ時間Ｔｌを種々変更して行った睡眠実験から、特徴的な信号波形が最も感度よく出現する値として選択されたものである。【0071】
周波数変動時系列解析演算手段（周波数変動時系列解析演算ステップ）６１３は、周波数演算手段６１１により得られたエアパック信号の周波数の時系列波形（図９の[2]のステップ）に、所定の時間幅の時間窓（好ましくは１８０秒）を設定し、周波数の平均値を求める（図９の[3],[4]のステップ）。次に、所定のオーバーラップ時間（好ましくは１６２秒）で設定した所定の時間窓（好ましくは１８０秒）毎にエアパック信号の周波数の平均値を求める移動計算を行い、プロットする。そして、時間窓毎にプロットされた周波数の平均値の時系列変化を周波数変動時系列波形として出力する（図９の[7]のステップ）。そして、周波数傾き時系列波形と周波数変動時系列波形とを併せて出力すると図９[9]のステップのようになる。なお、周波数変動時系列波形から心拍数が求められる。【0072】
ここで、ピーク検出法は、生体信号の波形の乱れ、波形の状態を検出しているのに対し、ゼロクロス法は周波数成分を検出している。すなわち、両者が一致している時は、周期関数に近似していることを意味し、乖離している時は、波形に複数の成分が重畳されていることを意味している。波形の乱れは、低周波成分に高周波成分が重畳されたり、低周波成分にさらに低周波成分が重畳されたりしている状態を示し、振幅の増大や減少となり、交感神経の亢進や低下に対応し、これが過緊張・我慢、放心、リラックス、慢然、賦活化の各状態に関係する。従って、ピーク検出法及びゼロクロス法のそれぞれにより求めた周波数傾き時系列波形を絶対値処理して両者を比較することにより、交感神経系が亢進した過緊張状態や我慢状態にあるか、副交感神経優位のリラックス状態にあるかの大雑把な状態判定を行うことができる。その一方、周波数傾き時系列波形から記述関数を求める場合、並びに、周波数変動時系列波形から心拍数の増減等を求める場合には、ゼロクロス法により得られる周波数の時系列波形を用いることが好ましい。これは、記述関数の算出目的が、波形の変化度合いを見ることで、活動代謝にあるか、安静代謝の中での変化の状態にあるか、あるいは、各代謝の中での活性度合・機能の低下度合を知ることにあるためである。【0073】
微分手段（微分ステップ）６１４は、周波数傾き時系列解析演算手段６１２により求めた周波数傾きの時系列波形を微分し、周波数傾き時系列波形の瞬間的な変動率を見る。【0074】
積分手段（積分ステップ）６１５は、周波数傾き時系列解析演算手段６１２により求めた前記周波数傾きの時系列波形を積分し、周波数傾き時系列波形の時間経過による累積を見る。【0075】
方形波算出手段（方形波算出ステップ）６１６は、周波数傾き時系列解析演算手段６１２により求めた周波数傾きの時系列波形のピークを求め、包絡線を作成し、包絡線の傾きから増減傾向を求めて方形波を描く。すなわち、包絡線の傾きが増加傾向の区間はプラス、減少傾向の区間はマイナスとして方形波を描く。【0076】
記述関数算出手段（記述関数算出ステップ）６１７は、任意に設定した第１時間帯における周波数傾き時系列波形を高速フーリエ変換すると共に、第１時間帯よりも後の第２時間帯における周波数傾き時系列波形を高速フーリエ変換して、「第２時間帯における高速フーリエ変換（ａ２＋ｉ・ｂ２・・・（ａ２は実数部、ｂ２は虚数部）／第１時間帯における高速フーリエ変換（ａ１＋ｉ・ｂ１・・・（ａ１は実数部、ｂ１は虚数部）」により両者間の記述関数（等価伝達関数）を求める手段である。つまり、第１時間帯における周波数傾き時系列波形を入力関数とし、第２時間帯における周波数傾き時系列波形を出力関数として、両者間の記述関数を求めることにより、基本の状態が相殺され、変化の様子が残り、第１時間帯と第２時間帯における人の恒常性維持機能の変化を捉えるものである。なお、第１時間帯及び第２時間帯は、各時間帯における周波数傾き時系列波形の関数を求めることができる時間幅であればよく限定されるものではなく、例えば、半周期分から３周期分（例えば約３分間から約２０分間）の範囲で設定できる。統計的に所定の範囲に設定してもよいし、個人毎に第１時間帯及び第２時間帯の時間幅を設定できるようにしてもよい。【0077】
得られた記述関数は「（ａ１ａ２＋ｂ１ｂ２）／（ａ１２＋ｂ１２）＋ｉ（ａ１ｂ２−ａ２ｂ１）／（ａ１２＋ｂ１２）」で表される。そして、「（ａ１ａ２＋ｂ１ｂ２）／（ａ１２＋ｂ１２）」＝Ａ、「ｉ（ａ１ｂ２−ａ２ｂ１）／（ａ１２＋ｂ１２）」＝Ｂとおき、「Ａ２＋Ｂ２」の平方根を記述関数振幅値として求める。そして、この記述関数振幅値により、２つの時間帯における状態変化を捉える。値の設定は任意であり個人毎に設定することもできるが、統計的に所定の値に設定するようにしてもよい。すなわち、所定の値あるいは値の範囲にある場合に、安静状態、漫然状態、交感神経代償作用の働いている状態等という設定を行うことができる。【0078】
絶対値処理手段（絶対値処理ステップ）６１８は、周波数傾き時系列解析演算手段６１２により得られる周波数傾き時系列波形を絶対値処理して出力する。このとき、絶対値処理手段６１８は、ピーク検出法を用いた場合とゼロクロス法を用いた場合の双方で得られた周波数傾き時系列波形を絶対値処理して出力する。具体的には、周波数演算手段６１１が、ピーク検出法を用いて生体信号の周波数の時系列波形を求めると共に、ゼロクロス法を用いて生体信号の周波数の時系列波形を求める。そして、周波数傾き時系列解析演算手段６１２が、ピーク検出法及びゼロクロス法による各周波数の時系列波形のそれぞれについて、周波数傾き時系列波形を求め、各周波数傾き時系列波形について絶対値処理手段６１８が絶対値処理する（図９[10]のステップ）。【0079】
上記のように、ピーク検出法は低周波成分に高周波成分が付加された生体信号の波形の乱れを検出しているのに対し、ゼロクロス法は生体信号の低周波成分を検出しており、波形の乱れは少ない。すなわち、両者が一致している時は、重畳成分が少ないことを意味し、乖離している時は、重畳成分が多く含まれていることを意味している。従って、絶対値処理手段６１８により得られた２つのデータを比較することで、交感神経系が亢進した過緊張状態あるいは我慢状態にあるか、交感神経の低下により副交感神経優位となるか、交感神経の低下を伴わずに副交感神経が優位となりリラックス状態に誘導されるか、あるいは交感神経・副交感神経共に低下することによる機能低下状態又は休息必要状態にあるかの状態判定を行うことができる。すなわち、絶対値処理手段６１８により、ゼロクロス法を用いて算出した絶対値がピーク検出法を用いて算出した絶対値よりも高い場合には交感神経亢進の状態であり、逆の場合は副交感神経優位な状態であると判定できる（図８（ａ）参照）。【0080】
恒常性維持機能段階算出手段（恒常性維持機能段階算出ステップ）６１９は、周波数傾き時系列解析演算手段６１２により求めた周波数傾き、微分手段６１４により求めた微分値、積分手段６１５により求めた積分値、方形波算出手段６１７により求めた方形波の正負、記述関数算出手段６１８により求めた記述関数振幅値、及び絶対値処理手段６１９により求めた周波数傾き時系列波形の２つの絶対値のうち、少なくとも１つ以上を用いて恒常性維持機能レベルの段階を求める。これらの組み合わせにより、上記の恒常性維持機能段階１〜５のいずれに該当するかを求める。例えば、周波数傾きと積分値を用いて、所定以上の場合に「恒常性維持機能段階１」と判定し、あるいは、微分値が所定位置以下であって、方形波が「負」、記述関数振幅値が所定値以上、かつ、２つの絶対値のうちの「ピーク優位」の場合に「恒常性維持機能段階４」と判定するように設定できる。これらの組み合わせ、判定の際の閾値等は限定されるものではなく、複数の被験者のデータを統計処理して決定することもできるし、個人毎に設定することもできる。【0081】
なお、恒常性維持機能段階算出手段６１９による状態判定は、各指標を次のように組み合わせて行うことが好ましい。【0082】
（１）恒常性維持機能段階１の判定：・ゼロクロス法による周波数傾き時系列波形が正の場合、
又は、
・ゼロクロス法による周波数傾き時系列波形の積分波形が正の場合
これらの波形が正の場合、人の交感神経活動が優位となっているため、高活性状態と判定するのに適しているからである。
【0083】
（２）恒常性維持機能段階２の判定：・ゼロクロス法による周波数傾き時系列波形が負の場合において、今回の傾きと前回の傾きとの比が所定値以上で、この比の出現が所定回数以下の場合、
又は、
・積分波形が所定の値の範囲の場合
周波数傾き時系列波形が負の場合であっても、所定値以上であればあるいは積分波形が所定の値の範囲内であれば、ある程度の活性状態が維持されているからである。
【0084】
（３）恒常性維持機能段階３の判定：・ゼロクロス法による周波数傾き時系列波形が負の場合において、今回の傾きと前回の傾きとの比が所定値以上で、この比の出現が所定回数以上の場合、
又は、
・積分波形が所定の値の範囲（恒常性維持機能段階２の場合よりも低い範囲）の場合、
又は、
・微分波形が所定の値以下の場合
これらの範囲で交感神経と副交感神経のバランスが良い状態を示すからである。
【0085】
（４）恒常性維持機能段階４の判定：・微分波形が所定の値以下（恒常性維持機能段階３の場合よりも低い値）の場合、
又は、
・積分波形が所定の値の範囲（恒常性維持機能段階３の場合よりも低い範囲）の場合、
又は、
・傾き時系列波形の記述関数振幅値が所定値以上であって、方形波が負、かつ、ピーク検出法で算出した絶対値がゼロクロス法による絶対値よりも優位の場合
これらの範囲で、副交感神経活動が徐々に優位になる傾向を示すからである。
【0086】
（５）恒常性維持機能段階５の判定：・微分波形が所定の値以下（恒常性維持機能段階４の場合よりも低い値）の場合、
又は、
・積分波形が所定の値以下（恒常性維持機能段階４の場合の下限値よりも低い値）の場合、
又は、
・傾き時系列波形の記述関数振幅値が所定値以上（恒常性維持機能段階４の場合よりも高い値）であって、方形波が負、かつ、ピーク検出法で算出した絶対値がゼロクロス法による絶対値よりも優位の場合
これらの範囲で、副交感神経活動が優位となり、機能低下によって疲労を感じるからである。
【0087】
（６）入眠予兆の判定：・ゼロクロス法で算出した絶対値がピーク検出法による絶対値よりも優位の場合であって、方形波が正、かつ、恒常性維持機能段階を示すグラフが２段階以上下がる場合
ゆらぎの範囲を逸脱する大きな動きであり、このような現象が生じた後にマイクロスリープが出現するからである。
【0088】
ここで、本出願人が従来行っている周波数傾き時系列波形と周波数変動時系列波形とを用いた判定は、入眠予兆信号という人の代謝変化に起因する状態を主として判定するものであった。つまり、時間経過と共に一方的に肉体が疲れていく過程を主として見ているものであった。しかしながら、本実施形態の恒常性維持機能段階算出手段６１９により表示される図１１（ａ）〜（ｅ）のような折れ線グラフによる表示は、上記したように、周波数傾き、微分値、積分値、方形波の正負、記述関数振幅値、及び絶対値の各指標を複数組み合わせて行われるため、代謝作用による人の状態変化だけでなく、各状態における恒常性維持機能の復元性能、減衰性能が示唆され、かつ、脳の刺激に伴う賦活化及び状態変化も併せて見ることができる。そのメカニズムを説明すると次の通りである。【0089】
まず、睡眠状態は代謝のエネルギーが最少レベルにあり、平常・運動状態ではエネルギーレベルが高位の状態にあり、安静状態はその中間に位置する。リラックス状態での定常運転状態では安静状態に近いエネルギーレベルになる。一方、睡眠との戦い、周辺状況への配慮・注意、危険運転などでは心拍数が上昇し、平常・運動状態に近いものとなる。各状態における恒常性の維持は、個人の身体・精神能力によって決まるゆらぎ幅を逸脱しないように調整が行われている。【0090】
人の状態変化は、低次元カオスによって、特徴抽出が可能であることが知られている。睡眠状態での恒常性維持機能のゆらぎはカオティックに変動し、平常・活動状態ではランダムあるいは周期的（以下、非カオス的と呼ぶ）に変動する。この非カオス的な変動を示している状態では、外的ストレスなどの環境変化への適用性が乏しくなる。ドライバーが車を運転している状態で、カオティックな変動を示している状況の下では、ドライバーは外部振動を知覚し難くなり、呼吸や心拍のゆらぎに影響を受け難い。すなわちカオティックな変動は恒常性維持のためのエネルギー消費量を最少限にし、疲労の進行を抑制する。この状態のドライバーはリラックス状態に誘導され、交感神経の過度な緊張が軽減される。一方、非カオス的な変動状態ではドライバーに外部振動が伝わりやすく、呼吸や心拍に影響を与え、各器官に無駄な力が入り、の抵抗疲労感を生む。この疲労感がエネルギー消費量を増大させ、疲労の進行の勾配が大きくなる。ここで人は、交感神経優位な状態となる。【0091】
人体の恒常性維持機能をばね−マス−ダッシュポット系で置き換えて考えてみた。ばね系はエネルギーを蓄え力に変換する機能で、これを交感神経系とした。一方、ダッシュポット系はエネルギーを減衰させる機能で、これを副交感神経系とした。マスは、生命維持機能である心循環器系とした。生命維持機能の恒常性維持を受け持つ中枢系と、外界からの刺激を緩和する末梢系の恒常性維持機能の２つの機能で生体の制御を行うこととした。これら２つの機能の特性は非線形性を含むため、記述関数として表される。背部に置かれたエアパックセンサは、主に中枢系に近いところの状態を捕捉する。エアパックセンサで捉えられた信号は、心臓の動きによる大動脈の揺れである。【0092】
人は恒常性の維持を自律神経系で制御する。副交感神経は、大きな変動は起こさない。従って、副交感神経優位は、交感神経の機能低下で生じる。つまり、心拍数の安定性は交感神経と副交感神経のバランスで制御される。心拍数の過度の上昇は副交感神経で抑制される（結果として迷走迷走神経反射が起こり、徐脈になる。）。人は交感神経と副交感神経のバランスのとれた自然体の状態、交感神経が亢進した緊張状態、交感神経も副交感神経も機能低下する疲労状態、副交感神経優位な回復期にあたるリラックス状態があり、これらの組み合わせが変化に対応する柔軟性を産み、恒常性を維持する助けとなる。【0093】
これらのことから、人の状態を推定していくには、自律神経系の変化を逆問題として解いていくことであり、生体のゆらぎの微分情報、交感神経と副交感神経のゆらぎの出現度合いから、まず緊張状態かリラックス状態かを判定する（図８（ａ））、次にゆらぎ波形の特徴のある変化を見て、活性状態、安静状態、機能低下状態あるいはそれらの遷移状態にあるかを判断する（図８（ｂ））。次いで、過去の状態の組み合わせから、その後の状態を推定する（図８（ｃ））。過去の状態の組み合わせとは、交感神経の亢進を表すゆらぎ波形の突出の状態とその出現回数、疲労の急激な進行（恒常性維持機能の急激な変化）の状態、その急激変化の回数、心拍の周波数の変動の傾向を表す方形波の情報である。これが入眠予兆現象などの推定に役立つ。【0094】
一方、精神疲労と肉体疲労の寄与率に応じた疲労の仕方を考えた場合、従来、疲労の自然な進行として、疲労感を感じていない段階、疲労感を交感神経の代償作用で補い疲労感を感じていない状態、疲労感を覚え、ヒューマンエラーが発生する段階の３つの段階があると言われている。【0095】
しかし、緊張状態では、精神疲労の寄与率が高く、リラックス状態では肉体疲労が疲労の進行を支配すると考えられる。ところが、従来の疲労の進行の考え方では、脳機能の活性化による疲労の抑制が考慮されていない。一般的な運転環境下では、疲労の進行度合いに拘わらず、脳機能の活性化が生じる。この脳機能の活性化には創造性や好奇心などがあり、これらは覚醒誘導をもたらす。【0096】
このことは、図１１（ａ）〜（ｅ）に示したように、生体ゆらぎ（代謝作用）による増減だけでなく、恒常性維持機能段階が急に元に戻る過程が生じたり、時間経過によっても恒常性維持機能段階がほとんど進行しない過程が見られることがわかる。つまり、単に、人の肉体が疲れている様を示しているのではなく、疲労に対する許容量の拡大過程をも示すグラフになっていると考えられる。疲労に対する許容量の拡大は、脳から指令によるもの、つまり脳の賦活化に起因する活性化によるものであり、このグラフは、代謝作用による人の状態変化だけでなく、脳の作用による人の状態（特に、「心のありさま（情態）」）の変化をも捉えているものと言える。従って、本発明において、「人（生体）の状態」の推定には、代謝作用に起因した状態変化だけでなく、恒常性維持機能のゆらぎの状態を通して、その性能劣化の様子から疲労の進行を推断すること、また、脳の賦活化に起因する再活性による状態変化を推定することも含む。【0097】
上記生体状態推定装置６０による状態推定過程を図７のフローチャートに基づき説明すると、まず、生体信号測定手段１により背部の大動脈の動きによるエアパック信号を得る（Ｓ１０１）。なお、Ｓ１０１のエアパック信号は、エアパックの圧力変動の信号そのものではなく、所定のフィルタリング（生体信号以外のノイズ信号を低減するアナログフィルタリング処理）を施した後の信号である。このエアパック信号を周波数演算手段６１１により処理することもできるが、より精度を高めるために、このエアパック信号にさらにフィルタリング手段を施し（Ｓ１０２）、解析に必要な周波数領域の信号波形（エアパック信号デジタルフィルタ波形）を得ている（Ｓ１０３）。このさらなるフィルタリング手段（以下、「デジタルフィルタリング手段」という）は、自動車の走行時において車体フロアから入力される振動（ノイズ）をできるだけ除去するものである。そのため、このデジタルフィルタリング手段は、静的環境下（なお、本明細書でいう「静的環境下」には、エンジンが始動しておらず全く振動の入力がない測定環境のほか、実質的な静的環境下（すなわち静的相当環境下）と言えるアイドリング状態（走行時のように路面凹凸に起因する振動が入力されない状態）での測定環境下も含む意味である）において、数秒間から数十秒間、生体信号測定手段１により生体信号を測定して、その平均周波数を求め、該平均周波数を基準値としてフィルタリングする周波数帯域（下限周波数と上限周波数）を設定する構成である。これにより、走行時に路面から車体フロアを介して入力される振動波形の混入が少なくなり、人の恒常性維持機能段階が同じとした場合における、上記デジタルフィルタリング手段を適用した静的環境下において得られるエアパック信号フィルタ波形と走行時に得られるエアパック信号デジタルフィルタ波形とがほぼ同じ波形となる。【0098】
なお、デジタルフィルタリング手段において基準となる平均周波数の求め方は限定されるものではないが、本実施形態では、生体信号測定手段１により得られた生体信号そのものの波形（すなわち圧力波形）を微分（好ましくは２階微分）し、微分波形に全波整流を適用し、さらにローパスフィルターを適用する。ローパスフィルターのフィルタリング周波数は、背部大動脈の動きによる生体信号の一般的な周波数である１〜３Ｈｚの範囲で適宜に設定される。ローパスフィルターを適用した波形からピークを求め、さらにピークの時間間隔から平均周波数を求める。この手法によるフィルタリングを施して得られた波形を図１０（ａ），（ｂ）に示す。図１０（ａ）は静的相当環境下であるアイドリング状態におけるエアパック信号波形とエアパック信号デジタルフィルタ波形との比較であり、図１０（ｂ）は運転時（走行時）におけるエアパック信号波形とエアパック信号デジタルフィルタ波形との比較である。この図から、エアパック信号波形は、運転時（走行時）においては静的相当環境下であるアイドリング状態と異なり、上下に突出する乱れた波形を示しているのに対し、エアパック信号デジタルフィルタ波形の場合は、静的アイドリング状態と運転時（走行時）とにおいてほぼ同じ波形となっている。従って、走行前のアイドリング状態（静的相当環境下）で予め生体信号を測定し、それを基準としてフィルタリングすることは高い精度の測定結果を得るために非常に好ましい。【0099】
Ｓ１０３において得られたエアパック信号フィルタ波形は、周波数演算手段６１１によって処理され、該エアパック信号フィルタ波形の周波数の時系列波形がゼロクロス法及びピーク検出法により求められる（Ｓ１０４）。次に、ゼロクロス法によるエアパック信号フィルタ波形の周波数の時系列波形を用いて、周波数傾き時系列波形、その微分波形（微分値）、積分波形（積分値）、方形波、記述関数（記述関数振幅値）を求めると共に、ゼロクロス法の周波数傾き時系列波形とピーク検出法による周波数傾き時系列波形を絶対値処理した波形を求める（Ｓ１０５）。次に、恒常性維持機能段階算出手段６１９により恒常性維持機能段階を判定し（Ｓ１０６）、出力手段６２により表示手段６５に出力する（Ｓ１０７）。【0100】
（試験例１）本実施形態に係るエアパック１０を備えた生体信号測定手段１により検出された背部の大動脈の動きによる圧力変動の信号（エアパック信号）の有効性を検証する試験として、心電図との比較を行った。
【0101】
被験者は２０歳代〜３０歳代の計４名(男性２名、女性２名)である。測定項目は心電図と、人体各部位に当接したエアパックから得られた圧力変動の信号（エアパック信号）である。エアパック信号は、上記実施形態の自動車用シート５００とベッド用マットレスを使用して、頸部、上腕、橈骨、大腿から採取した。実験は座位、仰臥位にてそれぞれ３分間計測して行った。【0102】
図１４は、３０歳代女性被験者の心電図とエアパック信号の比較結果であり、図１４（ａ）はエアパック信号と心電図の原波形の一部を示す。心電図とエアパック信号の原波形を比較すると心電図のピーク値（Ｒ波）はエアパック信号の極小値とほぼ一致していることがわかる。図１４（ｂ）はエアパック信号と心電図から算出した心拍周波数時系列波形を比較した結果を示す。心拍周波数は座位、仰臥位共に心電図と背部、頸部で波形の一致がみられる。他の３名も同様の結果であり、特に心電図と背部エアパック信号は安定した一致をみせた。また、仰臥位の状態においても、体動により一時的に乱れても、その後心電図との高い相関を示し、乱れたまま崩れていくことはなかった。【0103】
背部エアパック信号は心臓に最も近い部位から得た生体信号であり、心電図との相関性も高いことから、背部に設けられたエアパックセンサは心臓の動きによる大動脈の揺れを捕捉していると考えられる。また、心電図とエアパック信号の原波形の極性反転は、心臓が収縮している収縮期の動脈血管の拡張と心臓が弛緩して、血液の流入により拡張している拡張期の動脈血管の収縮の関係に相当するものと考えられる。図１４（ｂ）上のＡ，Ｂ，Ｃ部に示す心電図とエアパック信号の大きなずれは体動によるものである。つまりエアパック信号は、体動による変化も敏感に捉えることができると考えられる。そして、拘束性の高い自動車用シートの被験者の微小な体動は図１４（ａ）のＡ，Ｂ，Ｃ部のように捉えることはできないが、大きな運転動作による体動（体を前に大きくせり出した場合など）は図１４（ｂ）のＡ，Ｂ，Ｃ部のように捉えることができる。これらの体動は、人が通常とは異なる不安定状態、あるいは異常状態になった時に生じることから、エアパックセンサにより信号を検知できなくなった途端に何らかの異常事態であると判定することもできる。【0104】
（試験例２）−解析例１−
（試験条件）
５０歳代の男性被験者Ａ、４０歳代の女性被験者Ｂ及び４０歳代の男性被験者Ｂの３名を、広島県東広島市から富山県黒部市まで自動車で３日間かけて往復させた。１泊目は黒部市に宿泊し、２泊目は京都市に宿泊した。１日目は、東広島市から黒部市に向かい、２日目は、黒部市から京都市に向かい、３日目は、京都市から大阪市に向かって打ち合わせを行った後に東広島市に向かった。また、３名の被験者は中途で適宜に運転を交代しながら試験した。試験に用いた車両の運転席には、本実施形態のエアパック１０をシートバック部５１０に内蔵させたシート５００を搭載した。結果を図１５〜図３０に示す。各図において（ａ）は、ゼロクロス法を用いた周波数傾き時系列解析演算手段６１２から得られた周波数傾き時系列波形（図では「エアパック脈波周波数傾き」と表示）と、周波数変動時系列解析演算手段６１３から得られた周波数変動波形（図では「エアパック脈波周波数変動」と表示」）を示した。（ｂ）は、（ａ）の周波数傾き時系列波形を用いて、記述関数算出手段６１７により得られた記述関数（伝達関数）の推移を示し、（ｃ）は、（ａ）の周波数傾き時系列波形を用いて、方形波算出手段６１６から得られた方形波の推移を示す。（ｄ）は、（ａ）〜（ｃ），（ｆ）の各データ等を用いて、恒常性維持機能段階算出手段６１９により算出された恒常性維持機能段階の推移（判定結果）を示したものである。（ｅ）は、絶対値処理手段６１８により得られたピーク検出法による絶対値とゼロクロス法による絶対値とをそれぞれ積分したグラフである。（ｆ）は、絶対値処理手段６１８により得られたピーク検出法による絶対値とゼロクロス法による絶対値との時系列の推移を示したものでる。（ｇ）は、周波数傾き時系列波形の積分波形を示したものである。
【0105】
（考察）主として、各図の（ｄ）に示した判定結果と各被験者の実際の様子とを基に説明する。
（１）男性被験者Ａ 黒部方面１（往路）（図１５）
計測開始〜２０００秒前後までは疲労感はあるものの、比較的調子の良い状態が続いている。疲労感の軽減は生体ゆらぎによるもので、交感神経の代償作用は生じていないと思われる。その後入眠予兆現象が発現し、疲労が進行して行った。そして、３０００秒強で眠気が発生し運転を交代した。
【0106】
（２）女性被験者Ｂ 黒部方面２（往路）（図１６）１８００〜３０００秒間は調子の良い状態が続いた。それ以外は混雑（追い越し車線をのろのろ走行の車に車線をふさがれて）でイライラ運転であった。４２００秒強で眠気の発生が認められる。
【0107】
（３）男性被験者Ａ 黒部方面３（往路）（図１７）最初の６００秒間は調子の良い状態が続いたが、その後渋滞が発生し疲労状態に入った。この間の疲労は精神的なものが大きく、５１００秒強で最初の入眠予兆信号が発現した。その後は交感神経の代償作用が機能し、９６００秒強で二回目の入眠予兆信号が生じた。渋滞のイライラと夜間運転と緊張による疲労感で、２回の入眠予兆信号の発現で夕食をとった。
【0108】
（４）男性被験者Ａ 黒部方面４（往路）（図１８）夕食後元気が出たが、その元気は長続きせずに１２００秒強で眠気に対抗する交感代償作用が発現した。その後は２０００秒までは交感神経の亢進で元気な状態が続いた。途中交感神経のレベルが下がって副交感優位な状態になり、１回目の入眠予兆信号が生じた。２３００秒以降は緊張感の低下と共に、自然な疲労の進行が認められた。
【0109】
（５）女性被験者Ｂ 黒部方面５（往路）（図１９）計測開始時から交感代償作用が発現し、緊張を伴った疲労状態にある。１０００秒強から１８５０秒の間で、リラックス状態で脳機能の活性が生じたと思われる区間があり、その後疲労が回復傾向になり、やや疲労感があるものの高活性状態に誘導されている。これは脳機能の活性がリラックス状態に誘導し、疲労の回復を促したと思われる。
【0110】
（６）男性被験者Ａ 黒部方面６（往路）（図２０）助手席における睡眠の効果があり１２００秒強までの間、疲労感を感じていない状態がある。その後自然な覚醒度の低下があるが、これはサーカディアンリズムの低下に連動するものと思われる。
【0111】
（７）男性被験者Ａ 京都方面１（復路）（図２１）開始より７００秒間は、前日の疲れの影響がある。その後は、道に迷ったために色々と考えている。この思慮による効果で脳機能の活性化が生じたものと思われる。脳機能の活性が収まると、反動で疲労が進行していった。
【0112】
（８）男性被験者Ａ 京都方面２（復路）（図２２）打ち合わせ後の疲労感が生体ゆらぎを大きくし、一気に疲労状態に陥っている。男性被験者Ａは空腹感を覚えていたが、１１００秒付近でパン屋を発見したことに起因し、一気に覚醒度が上昇している。
【0113】
（９）男性被験者Ａ 京都方面３（復路）（図２３）パンを食することで注意散漫になり、運転注意が表示された。その後リラックス状態の中で緩やかに疲労感が軽減されていった。途中交感神経のレベルが下がり、二度ほど眠気が生じた。その後生体ゆらぎが生じ再度覚醒方向に振れていった。
【0115】
（１１）男性被験者Ａ 京都方面５（復路）（図２５）当初から疲れがあり、ゆらぎながら恒常性を維持しようとしている。機能低下が認められ交感代償作用が２４００秒前後から発現した。その後、眠気が生じたが、交感代償作用により持ち直した。その後、緊張が解けてリラックス状態になることにより、疲労感が出てきたものと思われる。
【0116】
（１２）男性被験者Ａ 京都方面６（復路）（図２６）休憩後運転を再開したが休憩中仕事上の電話を受けていた。運転開始当初は、休憩前の状態にすばやく戻っていった。その後は、電話内容に起因すると思われるイライラが続いている。４６００秒以降は交感神経が亢進したままで、急速な機能低下と眠気が各一回生じた。
【0117】
（１３）男性被験者Ａ 京都方面７（復路）（図２７）休憩後運転を再開した。休憩後、緊張と交感神経の亢進が解けて、一気に疲労感生じた。その疲労感の解消にリラックス状態の中で栄養ドリンクを飲むことで一時的に回復したが、その後は生体ゆらぎの幅が大きくなり、疲労感を軽減していった。栄養ドリンクはリラックス方向に誘導しながらの生体ゆらぎによる回復を促すと考えられる。
【0118】
（１４）男性被験者Ｂ 東広島方面１（復路）（図２８）疲労の蓄積は無い。揺らぎをうまく使いながら、緊張状態の中で疲労感に対処している。２５５０秒を超えたところで眠気が発生している。その後の機能低下は走行車線をふさぐ低速走行の車にイライラさせられることで生じた。
【0119】
（１５）男性被験者Ｂ 東広島方面２（復路）（図２９）食事と休息で疲労感が解消され、元気な状態になった。緊張感の反動で、眠気と大きな揺らぎが生じたが、交感代償作用による緊張状態で疲労感を感じていない状態にもどった。
【0120】
（１６）男性被験者Ｂ 東広島方面３（復路）（図３０）１２００秒を過ぎたところで急激に疲労感が出てきた。交感代償作用で一時的にもどるが、長続きしない。最後は到着に対する期待感で緊張状態になり、疲労感が小さくなった。
【0121】
上記のことから、各被験者の実際の様子と、各図の（ｄ）に示した判定結果から判断できる様子とがよく一致していると言える。【0122】
（試験例３）−解析例２−
（試験条件）
上記実施形態のエアパック１０を内蔵したクッションを大型トラックの運転席のシートバック部に取り付けて、夜間の長距離運転によるドライバの状態推定に関する実験を行った。実験は東京−大阪間の約５００ｋｍの高速道路の往復路で実施した。午後１０時頃東京を出発し、中間地点付近のサービスエリアで、１時間程度の休憩を取り、その後、大阪へ向けて実験を再開した。午前５時頃大阪に到着後、６〜８時間程度の睡眠を取り、午後１０頃に東京へ向けて復路の実験を行った。復路も往路同様に、中間地点付近のサービスエリアで、１時間程度の休憩を取った。被験者は３０歳代〜４０歳代の健康な男性計４名である。計測に使用した生体信号は、背部から採れる圧力波形である。
【0123】
（考察）図３１（ａ）〜（ｄ）は、被験者Ａの東京−大阪間の往復で行われた疲労実験の結果を示す。縦軸は恒常性維持機能段階１〜５(図では正常（高活性状態）〜異常（機能低下状態）)を示し、横軸は時間軸を示す。図３１（ａ）は往路の被験者Ａの東京−美合間の恒常性維持機能の変化の様子を示す。ゆらぎの中心軸は、恒常性維持機能段階２ないし３にあり、それを基準にしてゆらぎが生じている。また、恒常性維持機能段階１の占める割合が全行程の５２．５％になり、交感神経優位な緊張状態で頑張って運転していることが分かる。
【0124】
図３１（ｂ）は、復路の被験者Ａの大阪−小笠間の疲労進行度合を示す。図３１（ａ）と同様に図３１（ｂ）の被験者Ａは、恒常性維持機能段階２ないし３を基準にしてゆらぎが生じている。異なるのは、恒常性維持機能段階１にある時間帯が少なくなる事である。確実に疲労の進行が早くなっており、長時間睡眠による完全な疲労回復に至っていないものと思われる。【0125】
図３１（ｃ）は往路で休息後の被験者Ａの美合−大阪間の疲労進行度合を示す。この行程では、休息前の状態とほぼ同様の様相を示した。したがって、休息は疲労回復に効果があるものと思われる。図３１（ｄ）は復路で休息後の被験者Ａの小笠−東京間の恒常性維持機能の変化の様子を示す。この行程では、図３１（ｃ）の往路とは異なる現象を示すことから、休息は一時的な疲労回復に効果があるものの、蓄積された疲労には改善効果が無いことが伺われる。これら４つの事例から、往路で一時的であるが休息を挟む事によって恒常性維持機能の活性化が期待できることが示唆される。復路の休息は、疲労の蓄積により、改善効果は認められないものの、疲労の抑制効果は有るものと思われる。他の３名についても同様な結果が得られた。【0126】
なお、時間経過と共に疲労が進行している場合であっても、恒常性維持機能レベルを時系列にプロットして示されるグラフが、生体のゆらぎに対応したゆらぎの範囲である場合（その時点のゆらぎのほぼ中心位置に対して、恒常性維持機能レベルで１段階から２段階の範囲）には、特に問題はない。そこで、その範囲であるか否かを視覚的に把握することを容易にするため、図３２（ａ），（ｂ）に示したように、ゆらぎの範囲を明示するようにすることが好ましい。図３２（ａ）は中心に対するゆらぎの範囲を上下に伸びる矢印で示したものであり、図３２（ｂ）は、上限を示すＡラインと下限を示すＢラインとにより取り囲んで示したものである。【0127】
また、図３２（ｃ）に示したように、入眠予兆信号が生じたことを時系列のグラフ中に明確に表示するため、恒常性維持機能段階１〜５の下に、「入眠予兆」の段階を表示する構成とすることもできる。グラフの急低下が生じ、判定手段６３の異常状態判定手段により入眠予兆と判定された場合、そのときのグラフの下端を最下段に表示された「入眠予兆」段階まで引き延ばして表示する。これにより、入眠予兆信号が生じたタイミングを視覚的に容易に把握することができる。【0128】
なお、生体信号測定手段としては、上記したエアパック１０を用いたものに限らず、図３３に示したものを用いることもできる。図３３に示した生体信号測定手段２００は、三次元立体編物２１０、三次元立体編物支持部材２１５、フィルム２１６、板状発泡体２２１，２２２、振動センサ２３０を有して構成される。【0129】
三次元立体編物２１０は、図１等に示した生体信号測定手段１と同様のものを用いることができる。三次元立体編物２１０は、厚み方向の荷重−たわみ特性が、測定板上に載置して直径３０ｍｍ又は直径９８ｍｍの加圧板で加圧した際に、荷重１００Ｎまでの範囲で、人の臀部の筋肉の荷重−たわみ特性に近似したバネ定数を備えることが好ましい。具体的には直径３０ｍｍの加圧板で加圧した際の当該バネ定数が０．１〜５Ｎ／ｍｍの範囲、又は、直径９８ｍｍの加圧板で加圧した際の当該バネ定数が１〜１０Ｎ／ｍｍであるものを用いることが好ましい。人の臀部の筋肉の荷重−たわみ特性に近似していることにより、三次元立体編物と筋肉とが釣り合い、生体信号が伝播されると、三次元立体編物が人の筋肉と同様の振動を生じることになり、生体信号を大きく減衰させることなく伝播できる。【0130】
板状発泡体２２１，２２２は、ビーズ発泡体により構成することが好ましい。ビーズ発泡体としては、例えば、ポリスチレン、ポリプロピレン及びポリエチレンのいずれか少なくとも一つを含む樹脂のビーズ法による発泡成形体が用いることができる。ビーズ発泡体からなる板状発泡体２２１，２２２は、個々の微細なビーズを構成している発泡により形成された球状の樹脂膜の特性により、微小な振幅を伴う生体信号を膜振動（横波）として伝播する。この膜振動（横波）が三次元立体編物に弦振動として伝わり、これらの膜振動（横波）と弦振動が重畳され、生体信号は、膜振動（横波）と弦振動が重畳されることによって増幅された機械振動として、後述する振動センサ２３０により検出される。従って、生体信号の検出が容易になる。【0131】
板状発泡体２２１，２２２をビーズ発泡体から構成する場合、発泡倍率は２５〜５０倍の範囲で、厚さがビーズの平均直径以下に形成されていることが好ましい。例えば、３０倍発泡のビーズの平均直径が４〜６ｍｍ程度の場合では、板状発泡体２２１，２２２の厚さは３〜５ｍｍ程度にスライスカットする。これにより、板状発泡体２２１，２２２に柔らかな弾性が付与され、振幅の小さな振動に共振し、フィルム上を伝わる横波に減衰が生じにくくなる。なお、板状発泡体２２１，２２２は、本実施形態のように、三次元立体編物２１０を挟んで両側に配置されていても良いが、いずれか片側、好ましくは、シートバック側のみに配置した構成とすることもできる。【0132】
ここで、三次元立体編物２１０は、幅４０〜１００ｍｍ、長さ１００〜３００ｍｍの範囲の短冊状のものが用いられる。この大きさのものだと、三次元立体編物２１０に予備圧縮（連結糸に張力が発生する状態）を生じやすくなり、人と三次元立体編物２１０との間で平衡状態が作りやすい。本実施形態では、人が背部が当接した際の違和感軽減のため、脊柱に対応する部位を挟んで対象に２枚配設するようにしている。三次元立体編物２１０を簡単に所定位置に配置するようにするため、三次元立体編物２１０は三次元立体編物支持部材２１５に支持させた構成とすることが好ましい。三次元立体編物支持部材２１５は、板状に成形され、脊柱に対応する部位を挟んで対称位置に、縦長の配置用貫通孔２１５ａ，２１５ａが２つ形成されている。三次元立体編物支持部材２１５は、上記板状発泡体２２１，２２２と同様に、板状に形成されたビーズ発泡体から構成することが好ましい。三次元立体編物支持部材２１５をビーズ発泡体から構成する場合の好ましい発泡倍率、厚さの範囲は上記板状発泡体２２１，２２２と同様である。但し、生体信号により膜振動（横波）をより顕著に起こさせるためには、三次元立体編物２１０，２１０の上下に積層される板状発泡体２２１，２２２の厚さが、三次元立体編物支持部材２１５の厚さよりも薄いことが好ましい。【0133】
三次元立体編物支持部材２１５に形成した配置用貫通孔２１５ａ，２１５ａに、２つの三次元立体編物２１０，２１０を挿入配置した状態で、三次元立体編物２１０，２１０の表側及び裏側にフィルム２１６，２１６を積層する。なお、配置用貫通孔２１５ａ，２１５ａの形成位置（すなわち、三次元立体編物２１０，２１０の配設位置）は、心房と大動脈（特に、「下行大動脈」）の拍出に伴う動きによって生じる揺れ及び大動脈弁の動き（心部揺動波）を検知可能な領域に相当する位置とすることが好ましい。この結果、三次元立体編物２１０，２１０は、上下面が板状発泡体２２１，２２２によりサンドイッチされ、周縁部が三次元立体編物支持部材２１５によって取り囲まれており、板状発泡体２２１，２２２及び三次元立体編物支持部材２１５が共振箱（共鳴箱）の機能を果たす。【0134】
また、三次元立体編物支持部材２１５よりも、三次元立体編物２１０，２１０の方が厚いものを用いることが好ましい。つまり、三次元立体編物２１０，２１０を配置用貫通孔２１５ａ，２１５ａに配置した場合には、三次元立体編物２１０，２１０の表面及び裏面が、該配置用貫通孔２１５ａ，２１５ａよりも突出するような厚さ関係とする。これにより、フィルム２１６，２１６の周縁部を配置用貫通孔２１５ａ，２１５ａの周縁部に貼着すると、三次元立体編物２１０，２１０は厚み方向に押圧されるため、フィルム２１６，２１６の反力による張力が発生し、該フィルム２１６，２１６に固体振動（膜振動（横波））が生じやすくなる。一方、三次元立体編物２１０，２１０にも予備圧縮が生じ、三次元立体編物の厚さ形態を保持する連結糸にも反力による張力が生じて弦振動が生じやすくなる。なお、フィルム２１６，２１６は、三次元立体編物２１０，２１０の表側及び裏側の両側に設けることが好ましいが、いずれか少なくとも一方に設けた構成とすることも可能である。フィルム２１６，２１６としては、例えば、ポリウレタンエラストマーからなるプラスチックフィルム(例えば、シーダム株式会社製、品番「ＤＵＳ６０５−ＣＤＲ」）等を用いることができる。【0135】
振動センサ２３０は、上記したフィルム２１６，２１６を積層する前に、いずれか一方の三次元立体編物２１０に固着して配設される。三次元立体編物２１０は一対のグランド編地と連結糸とから構成されるが、各連結糸の弦振動がグランド編地との節点を介してフィルム２１６，２１６及び板状発泡体２２１，２２２に伝達されるため、振動センサ２３０は感知部２３０ａを三次元立体編物２１０の表面（グランド編地の表面）に固着することが好ましい。振動センサ２３０としては、マイクロフォンセンサ、中でも、コンデンサ型マイクロフォンセンサを用いることが好ましい。本実施形態では、マイクロフォンセンサを配置した部位（すなわち、三次元立体編物２１０を配置した配置用貫通孔２１５ａ）の密閉性を考慮する必要がないため、マイクロフォンセンサのリード線の配線は容易に行うことができる。生体信号によって生じる人の筋肉を介した体表面の振動は、三次元立体編物２１０だけでなく、板状発泡体２２１，２２２、フィルム２１６にも伝播され、それらが振動（弦振動、膜振動（横波））して減衰を防止しつつ重畳されて増幅する。よって、振動センサ２３０は、三次元立体編物２１０に限らず、振動伝達経路を構成する板状発泡体２２１，２２２及びフィルム２１６に、その感知部２３０ａを固定することもできる。【0136】
生体信号測定手段２００としては、図３３に示したものに限定されず、例えば、図３４（ａ）に示したように、２つの三次元立体編物２１０，２１０の両方を覆うことのできる大きさのフィルム２１７を少なくとも一方に用いても良い。また、図３４（ｂ）に示したように、略長方形の三次元立体編物を両端縁から中央部に向かって折り曲げ、重なり合った部分の中央部を縫い合わせたランバーサポート２１８を配置するようにしてもよい。ランバーサポート２１８は、三次元立体編物支持部材２１５に面ファスナ等を用いて固定される。ランバーサポート２１８をこのようにして設けることにより、狭いスペース中でストローク感を高めるのに寄与する。【0137】
上記した生体信号測定手段２００は、例えば、図３５に示したように、自動車用シート１０００のシートバックフレーム１１００に被覆される表皮１２００の内側に配置される。なお、配置作業を容易にするため、生体信号測定手段２００を構成する三次元立体編物２１０、三次元立体編物支持部材２１５、フィルム２１６、板状発泡体２２１，２２２、振動センサ２３０等は予めユニット化しておくことが好ましい。【0138】
上記した生体信号測定手段２００によれば、生体信号により、筋肉の荷重−たわみ特性に近似する荷重−たわみ特性をもつ板状発泡体２２１，２２２やフィルム２１６に膜振動（横波）が生じると共に、人の筋肉の荷重−たわみ特性に近似した荷重−たわみ特性を有する三次元立体編物２１０に弦振動が生じる。そして、三次元立体編物２１０の弦振動は再びフィルム２１６等の膜振動（横波）に影響を与え、これらの振動が重畳して作用する。その結果、生体信号に伴って体表面から入力される振動は、減衰することなく弦振動と膜振動（横波）との重畳によって増幅された固体振動として直接振動センサ２３０により検出されることになる。【0139】
図１等に示したエアパック１０内の空気圧変動を検出する生体信号測定手段１の場合、体積と圧力が反比例関係にあるため、密閉袋の体積を小さくしないと圧力変動を検出しにくい。これに対し、図３３〜図３５に示した生体信号測定手段２００によれば、空気圧変動ではなく、上記のように、機械的増幅デバイス（三次元立体編物２１０、板状発泡体２２１，２２２、フィルム２１６又はフィルム２１７）に伝播される増幅された固体振動を検出するものであるため、その容積（体積）が検出感度の観点から制限されることはほとんどなく、心部揺動波という振幅の小さな振動を感度良く検出できる。このため、多様な体格を有する人に対応できる。従って、図３３〜図３５に示した生体信号測定手段２００は、乗物用シートのように、多様な体格を有する人が利用し、さらに多様な外部振動が入力される環境下においても感度良く生体信号を検出できる。【産業上の利用可能性】
【0140】
本発明は、ベッドなどの寝具に生体信号測定手段を配置し、背部の大動脈の揺動を捉え、上記した生体状態推定装置により解析して、人の状態推定を行うことに適用することもできる。これにより、寝ている人（特に、病人、介護を要する人）の健康状態を表示手段のモニタに示される画面により容易に把握することができる。同様に、家庭内に配置される椅子、事務用椅子等に生体信号測定手段を配置して状態推定を行うことに適用することもでき、それにより、日常生活における健康状態を容易に把握することができる。【符号の説明】
【0141】
１ 生体信号測定手段１０ エアパック
１５ 収容体
１００ エアパックユニット
２０ ビーズ発泡樹脂弾性部材
４０，４５ 三次元立体編物
５００ シート
５１０ シートバック部
５１１ 表皮部材
５１２ クッション支持部材
５２０ シートクッション部
６０ 生体状態推定装置
６１ 恒常性維持機能レベル演算手段
６１１ 周波数演算手段（周波数演算ステップ）
６１２ 周波数傾き時系列解析演算手段（周波数傾き時系列解析演算ステップ）
６１３ 周波数変動時系列解析演算手段（周波数変動時系列解析演算ステップ）
６１４ 微分手段（微分ステップ）
６１５ 積分手段（積分ステップ）
６１６ 方形波算出手段（方形波算出ステップ）
６１７ 記述関数算出手段（記述関数算出ステップ）
６１８ 絶対値処理手段（絶対値処理ステップ）
６１９ 恒常性維持機能段階算出手段（恒常性維持機能段階算出ステップ）
６２ 出力手段
６３ 判定手段
６５ 表示手段