(57)【特許請求の範囲】
【請求項１】
生体信号測定装置により測定した人の生体信号を分析し、生体状態を推定する生体状態推定装置であって、
前記生体信号が前記生体信号測定装置により背部から採取される音・振動情報であり、
前記背部から採取される音・振動情報を分析して、脳機能、自律神経機能、肉体・精神疲労又は感覚との関連性の高いゆらぎに起因する指標を含む、生体調節機能要素の状態の変動を示す複数の指標を、それぞれ予め設定された所定の判定時間毎に算出してその時系列変化を求める生体調節機能要素判定手段と、
前記生体調節機能要素判定手段における前記各生体調節機能要素の各判定時間よりも長く設定される基礎的体調推定時間に対応した、前記人の基礎的な体調を推定する基礎的体調推定手段と、
前記基礎的体調推定手段によって推定される前記人の基礎的な体調のレベルを、前記所定の基礎的体調推定時間の経過毎に出力する基礎的体調出力手段と
を有し、
前記基礎的体調推定手段は、
前記生体調節機能要素判定手段により求められた前記各生体調節機能要素に関する各時系列変化のうち、優先順位の高い前記生体調節機能要素の時系列変化が、所定の基準を満たす場合に、前記所定の基礎的体調推定時間における前記人の基礎的な体調を、所定のレベルと推定する第１基礎的体調推定手段と、
前記第１基礎的体調推定手段における推定対象とならない場合に、前記第１基礎的体調推定手段において用いた前記生体調節機能要素よりも優先順位の低い他の生体調節機能要素に関する前記時系列変化を用い、前記人の基礎的な体調を、所定の基準に基づいて分類される所定のレベルと推定する第２基礎的体調推定手段と
を有し、
前記基礎的体調出力手段が、前記第１基礎的体調推定手段又は前記第２基礎的体調推定手段によって推定される前記人の基礎的な体調のレベルを出力することを特徴とする生体状態推定装置。
【請求項２】
前記第１基礎的体調推定手段は、前記生体調節機能要素判定手段において判定された複数の前記生体調節機能要素のうち、前記自律神経機能への関連性の高い指標又は前記肉体・精神疲労への関連性の高い指標の時系列変化を用い、その時系列変化が所定の基準を満たす場合に、前記人の基礎的な体調のレベルを、「良好」又は「不良」と推定する請求項１記載の生体状態推定装置。
【請求項３】
前記第１基礎的体調推定手段は、前記自律神経機能への関連性の高い指標の時系列変化が、所定の基準を満たす場合に、前記人の基礎的な体調を「不良」と推定し、前記肉体・精神疲労への関連性の高い指標の時系列変化が、所定の基準を満たす場合に、前記人の基礎的な体調を「良好」と推定する手段である請求項２記載の生体状態推定装置。
【請求項４】
前記第２基礎的体調推定手段は、前記第１基礎的体調推定手段において「良好」又は「不良」と推定されない場合に、前記生体調節機能要素判定手段において判定された複数の前記生体調節機能要素のうち、前記感覚への関連性の高い指標の時系列変化から、前記人の基礎的な体調のレベルを、「良好」、「不良」又はそれらの「中間状態」のいずれかと推定する請求項２又は３記載の生体状態推定装置。
【請求項５】
コンピュータを用いて、生体信号測定装置により測定した人の生体信号を分析し、生体状態を推定する生体状態推定方法であって、
前記生体信号が前記生体信号測定装置により背部から採取される音・振動情報であり、
前記背部から採取される音・振動情報を分析して、脳機能、自律神経機能、肉体・精神疲労又は感覚との関連性の高いゆらぎに起因する指標を含む、生体調節機能要素の状態の変動を示す複数の指標を、それぞれ予め設定された所定の判定時間毎に算出してその時系列変化を求める生体調節機能要素判定手順と、
前記生体調節機能要素判定手順における前記各生体調節機能要素の各判定時間よりも長く設定される基礎的体調推定時間に対応した、前記人の基礎的な体調を推定する基礎的体調推定手順と、
前記基礎的体調推定手順によって推定される前記人の基礎的な体調のレベルを、前記所定の基礎的体調推定時間の経過毎に出力する基礎的体調出力手順と
を有し、
前記基礎的体調推定手順は、
前記生体調節機能要素判定手順により求められた前記各生体調節機能要素に関する各時系列変化のうち、優先順位の高い前記生体調節機能要素の時系列変化が、所定の基準を満たす場合に、前記所定の基礎的体調推定時間における前記人の基礎的な体調を、所定のレベルと推定する第１基礎的体調推定手順と、
前記第１基礎的体調推定手順における推定対象とならない場合に、前記第１基礎的体調推定手順において用いた前記生体調節機能要素よりも優先順位の低い他の生体調節機能要素に関する前記時系列変化を用い、前記人の基礎的な体調を、所定の基準に基づいて分類される所定のレベルと推定する第２基礎的体調推定手順と
を有し、
前記基礎的体調出力手順が、前記第１基礎的体調推定手順又は前記第２基礎的体調推定手順によって推定される前記人の基礎的な体調のレベルを出力することを特徴とする生体状態推定方法。
【請求項６】
生体状態推定装置としてのコンピュータに、
生体信号測定装置により測定した人の生体信号を分析し、生体状態を推定する生体状態推定手順を実行させるコンピュータプログラムであって、
前記生体信号が前記生体信号測定装置により背部から採取される音・振動情報であり、
前記生体状態推定手順として、
前記背部から採取される音・振動情報を分析して、脳機能、自律神経機能、肉体・精神疲労又は感覚との関連性の高いゆらぎに起因する指標を含む、生体調節機能要素の状態の変動を示す複数の指標を、それぞれ予め設定された所定の判定時間毎に算出してその時系列変化を求める生体調節機能要素判定手順と、
前記生体調節機能要素判定手順における前記各生体調節機能要素の各判定時間よりも長く設定される基礎的体調推定時間に対応した、前記人の基礎的な体調を推定する基礎的体調推定手順と、
前記基礎的体調推定手順によって推定される前記人の基礎的な体調のレベルを、前記所定の基礎的体調推定時間の経過毎に出力する基礎的体調出力手順と
を実行させ、
前記基礎的体調推定手順は、
前記生体調節機能要素判定手順により求められた前記各生体調節機能要素に関する各時系列変化のうち、優先順位の高い前記生体調節機能要素の時系列変化が、所定の基準を満たす場合に、前記所定の基礎的体調推定時間における前記人の基礎的な体調を、所定のレベルと推定する第１基礎的体調推定手順と、
前記第１基礎的体調推定手順における推定対象とならない場合に、前記第１基礎的体調推定手順において用いた前記生体調節機能要素よりも優先順位の低い他の生体調節機能要素に関する前記時系列変化を用い、前記人の基礎的な体調を、所定の基準に基づいて分類される所定のレベルと推定する第２基礎的体調推定手順と
を実行し、
前記基礎的体調出力手順が、前記第１基礎的体調推定手順又は前記第２基礎的体調推定手順によって判定される前記人の基礎的な体調のレベルを出力することを特徴とするコンピュータプログラム。
【請求項７】
前記第１基礎的体調推定手順は、前記生体調節機能要素判定手順において判定された複数の前記生体調節機能要素のうち、前記自律神経機能への関連性の高い指標又は前記肉体・精神疲労への関連性の高い指標の時系列変化を用い、その時系列変化が所定の基準を満たす場合に、前記人の基礎的な体調のレベルを、「良好」又は「不良」と推定する請求項６記載のコンピュータプログラム。
【請求項８】
前記第２基礎的体調推定手順は、前記第１基礎的体調推定手順において「良好」又は「不良」と推定されない場合に、前記生体調節機能要素判定手順において判定された複数の前記生体調節機能要素のうち、前記感覚への関連性の高い指標の時系列変化から、前記人の基礎的な体調のレベルを、「良好」、「不良」又はそれらの「中間状態」のいずれかと推定する請求項７記載のコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、生体信号測定装置により測定した人の生体信号に基づき、人の基礎的な体調を推定する生体状態推定装置、生体状態推定方法及びコンピュータプログラムに関する。
【背景技術】
【０００２】
本発明者らは、特許文献１〜５等において、人の上体の中で背部の体表面に生じる振動を生体信号測定装置により検出し、人の状態を解析する技術を提案している。人の上体背部から検出される心臓と大動脈の運動から生じる音・振動情報は、心臓と大動脈の運動から生じる圧力振動であり、心室の収縮期及び拡張期の情報と、循環の補助ポンプとなる血管壁の弾力情報及び反射波の情報を含んでいる。すなわち、心臓と大動脈の運動から背部表面に生じる１Ｈｚ近傍の背部体表脈波（Ａｏｒｔｉｃ  Ｐｕｌｓｅ  Ｗａｖｅ（ＡＰＷ））を含む）や、心拍に伴って背部側に伝わる音（「疑似心音」（本明細書では胸部側から採取される心臓の音である心音に対して、背部側で採取される心臓の音を「疑似心音」とする））の情報を含んでいる。そして、心拍変動に伴う信号波形は交感神経系及び副交感神経系の神経活動情報を含み、大動脈の揺動に伴う信号波形には交感神経活動の情報を含んでいる。
【０００３】
特許文献１では、採取した生体信号（音・振動情報）から抽出した１Ｈｚ近傍の背部体表脈波（ＡＰＷ）の時系列波形に所定の時間幅を適用してスライド計算を行って周波数傾きの時系列波形を求め、その変化の傾向から、例えば、振幅が増幅傾向にあるか、減衰傾向にあるかなどによって生体状態の推定を行っている。また、生体信号を周波数解析し、予め定めたＵＬＦ帯域（極低周波帯域）からＶＬＦ帯域（超低周波帯域）に属する機能調整信号、疲労受容信号及び活動調整信号に相当する各周波数のパワースペクトルを求め、各パワースペクトルの時系列変化から人の状態を判定することも開示している。疲労受容信号は、通常の活動状態における疲労の進行度合いを示すため、これに併せて、機能調整信号や活動調整信号のパワースペクトルの優勢度合いを比較することにより、人の状態（交感神経優位の状態、副交感神経優位の状態など）を判定することができる。また、これら３つの信号に相当する周波数成分のパワースペクトルの値の合計を１００とした際の各周波成分の分布率を時系列に求め、その分布率の時系列変化を利用して人の状態を判定することも開示している。
【０００４】
特許文献２では、生体状態の定量化手法として、生体状態を体調マップ及び感覚マップとして表示する技術を提案している。これは、上記したＡＰＷを周波数分析し、対象となる解析区間について、解析波形を両対数軸表示に表し、その解析波形を低周波帯域、中周波帯域、高周波帯域に分け、区分けした解析波形の傾きと、全体の解析波形の形とから一定の基準に基づいて解析波形の点数化を行い、それを座標軸にプロットしたものである。体調マップは、自律神経系の制御の様子を交感神経と副交感神経のバランスとして見たものであり、感覚マップは、体調マップに心拍変動の変化の様子を重畳させたものである。
【０００５】
特許文献３〜５では、恒常性維持機能レベルを判定する手段を開示している。恒常性維持機能レベル判定する手段は、周波数傾き時系列波形の微分波形の正負、周波数傾き時系列波形を積分した積分波形の正負、ゼロクロス法を利用した周波数傾き時系列波形とピーク検出法を利用した周波数傾き時系列波形をそれぞれ絶対値処理して得られた各周波数傾き時系列波形の絶対値等のうち、少なくとも１つ以上を用いて判定する。これらの組み合わせにより、恒常性維持機能のレベルがいずれに該当するかを求める。例えば、周波数傾きと積分値を用いて、所定以上の場合に「恒常性維持機能レベル１」と判定し、あるいは、微分値が所定位置以下であって、かつ、２つの絶対値のうちの「ピーク優位」の場合に「恒常性維持機能レベル４」と判定するように設定できる。これらの組み合わせ、判定の際の閾値等は多数の被験者のデータを統計処理して決定している。また、恒常性維持機能のレベルは、例えば、５〜７段階に分け、恒常性維持機能の優れる場合（集中度合いの高い場合等）から、恒常性維持機能が劣る場合（過緊張状態の場合、脇見運転等による集中力の低下等）を判定する。モニタに表示するに当たっては、５〜７段階のレベルを文字で表示したり、あるいは、中間レベル（普通の状態）以上の場合には、一括して恒常性維持機能が優れる場合と判定し、それよりも下の場合には一括しての恒常性維持機能が劣る場合と判定し、それぞれについて、モニタに異なる色彩表示がなされるように設定したりすることも開示されている。
【０００６】
非特許文献１では、指尖容積脈波情報に関し、交感神経の情報を反映するパワー値の周波数傾き時系列波形を求め、それを絶対値処理した積分値を疲労度として時系列にプロットし、これにより疲労曲線を描き、筋疲労を捉える技術が開示されている。非特許文献２では、エアパックセンサを用いて人の背部から取得した生体信号を同様の手法で演算処理して疲労曲線を描き、筋疲労を捉える技術が開示されている。すなわち、交感神経の情報を反映したパワー値の周波数傾き時系列波形（ＡＰＷの場合にはゼロクロス法による周波数傾き時系列波形）を用いることによって筋疲労の状態を把握することができる。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０１１−１６７３６２号公報
【特許文献２】特開２０１２−２３９４８０号公報
【特許文献３】ＷＯ２０１１／０４６１７８号公報
【特許文献４】特開２０１４−１１７４２５号公報
【特許文献５】特開２０１４−２２３２７１号公報
【非特許文献】
【０００８】
【非特許文献１】原著「指尖容積脈波情報を用いた長時間着座疲労の簡易評価法の開発」藤田悦則（外８名）、人間工学  Ｖｏｌ．４０、Ｎｏ．５（’０４）
【非特許文献２】「非侵襲型センサによって測定された生体ゆらぎ信号の疲労と入眠予知への応用」、落合直輝（外６名）、第３９回日本人間工学会  中国・四国支部大会  講演論文集、平成１８年１１月２５日発行、発行所：日本人間工学会  中国・四国支部事務局
【発明の概要】
【発明が解決しようとする課題】
【０００９】
上記した技術は、いずれも、生体調節機能に関してゆらぎに起因する各要素を分析して人の状態を判定するものであるが、生体信号に対する演算処理がそれぞれ異なり、出力される判定結果として、入眠予兆のタイミングを判定したり、疲労の度合いを判定したり、恒常性維持機能レベルの変化を判定したり、それぞれの目的に応じたものとなっている。しかし、これらは、いずれも別々に出力される。特許文献４では、入眠予兆現象、切迫睡眠現象、覚低走行状態、恒常性維持機能レベル、初期疲労状態、気分判定など、生体調節機能のゆらぎに起因する各要素に関する複数の指標の時系列変化を１台の装置で判定し、それらを１つのモニタに出力する技術も開示しているが、いずれにしても、各指標の時系列変化を個別に判定していることに変わりない。
【００１０】
特許文献４に開示の技術のように様々な指標を求めることは、人の状態、特に運転者の状態をより正確に把握するために非常に重要であることはもとよりであるが、上記従来の技術においては、自律神経系に関する指標を基本的に用いて、それに対応する人の状態をより正確に出力することに重点がおかれているため、判定対象となる生体信号の抽出時間（判定時間）が比較的短く、数十秒から数分に１回といった頻度で判定結果を出力している。しかし、例えば、運転中、運転者あるいはその情報を遠隔的に監視する管理者において、３０分、１時間といったある程度長い時間に亘る基礎的な体調（その時間帯における支配的、代表的ないしは平均的な体調）を概略的に把握したいという要請もある。つまり、例えば、１時間運転したが、その間、基本的に、良好な体調のもとで運転していたと概略的に推定できるのか、体調不良の状態で恒常性を維持する生体調節機能に助けられて運転していたと概略的に推定できるのか、運転者自ら把握ができれば、早めに休憩をとるなどの自己判断を促すことができる。
【００１１】
本発明は、上記した点に鑑みなされたものであり、比較的短時間で次々と判定される生体調節機能別のゆらぎに起因する各要素（生体調節機能要素）の時系列変化の情報を複数用いると共に、これを所定の条件下で処理し、より長時間に設定された基礎的体調推定時間における人の基礎的な体調を出力し、運転者等による概略的な体調把握を容易にすることができる生体状態推定装置、生体状態推定方法及びコンピュータプログラムを提供することを課題とする。
【課題を解決するための手段】
【００１２】
上記課題を解決するため、本発明者が鋭意検討し、次の点に着目して本発明を完成するに至った。すなわち、人の上体背部から検出される心臓と大動脈の運動から生じる音・振動情報、特に、それらのうちの１Ｈｚ近傍の背部体表脈波（ＡＰＷ）は、血管の弾性情報や反射波の情報等を含んでいる。このため、背部体表脈波（ＡＰＷ）を解析して後述の周波数傾き時系列波形を求めることにより、生体調節機能要素である生体の総合的なゆらぎの情報を求めることができる。このゆらぎの情報の中で周波数毎の変動の様子を捉えると、周波数毎の分布率が脳波（θ波、α波、β波）のゆらぎの要素が反映された変動の仕方をとるため、その時系列波形を用いて解析することにより、脳波のどの周波数帯域が支配的なゆらぎなのかといったゆらぎの情報を捉えることができる。また、ゼロクロス法による周波数傾き時系列波形が、自律神経系の支配するところにある一方で、脳波のゆらぎを反映しているところは、ゼロクロス法による周波数傾き時系列波形を周波数解析し、その中で、０．００１７Ｈｚに代表される機能調整信号、０．００３５Ｈｚに代表される疲労受容信号、０．００５３Ｈｚに代表される活動調整信号の３点の周波数成分のパワースペクトル比で示されるものであり、いわば周波数解析したパワースペクトルの形を表す３点の周波数成分の分布率は、その急変する部位を自律神経反応というよりも内分泌系の機能の発現を示す部位として捉えられると考えられる。また、ゼロクロス法により求めた周波数傾き時系列波形は、絶対値処理することにより交感神経の発現の度合いを示し、ピーク検出法により求めた周波数傾き時系列波形は副交感神経の発現の度合いを示す。よって、これらを用いることで生体調節機能の発現の様子をより詳しく捉えることができる。例えば、ゼロクロス法により求めた周波数傾き時系列波形を絶対値処理し、これを積分することで人の疲労の度合いを示す疲労曲線が求められ、筋疲労の状態を把握できる。また、ゼロクロス法を用いた各周波数傾き時系列波形の微分波形の正負、ゼロクロス法又はピーク検出法の各周波数傾き時系列波形の絶対値等のうち、少なくとも１つ以上を用いることにより、恒常性維持機能レベルの変動の様子を捉えることができる。疲労曲線、恒常性維持機能レベルの時系列波形も、周波数傾き時系列波形から派生したものであり、自律神経系、脳機能等のゆらぎの情報を反映した指標となる。これに加え、肉体・精神疲労への関連性の高い指標（体調マップ、感覚マップ）、及び、感覚への関連性の高い指標（恒常性維持機能レベルの注意、警告に相当するレベルの頻度から求められる疲労として意識する感覚あるいは倦怠感等が生じている頻度）も求める。これらは、脳、自律神経系、内分泌系の各機能のゆらぎの様子を複数の観点から捉えているものであり、これらの複数種類のゆらぎの指標を後述のように統計的に処理することで、分析対象の人の基礎的な体調がどのような状態であるかを推定できると考えられる。
【００１３】
すなわち、本発明の生体状態推定装置は、生体信号測定装置により測定した人の生体信号を分析し、生体状態を推定する生体状態推定装置であって、前記生体信号を分析して、脳機能、自律神経機能、肉体・精神疲労又は感覚との関連性の高いゆらぎに起因する指標を含む、生体調節機能要素の状態の変動を示す複数の指標を、それぞれ予め設定された所定の判定時間毎に算出してその時系列変化を求める生体調節機能要素判定手段と、前記生体調節機能要素判定手段における前記各生体調節機能要素の各判定時間よりも長く設定される基礎的体調推定時間に対応した、前記人の基礎的な体調を推定する基礎的体調推定手段と、前記基礎的体調推定手段によって推定される前記人の基礎的な体調のレベルを、前記所定の基礎的体調推定時間の経過毎に出力する基礎的体調出力手段とを有し、前記基礎的体調推定手段が、前記生体調節機能要素判定手段によって求められる前記各生体調節機能要素の状態に関する各時系列変化を、予め設定された優先順位に従って分析し、所定の基準に照らして、前記所定の基礎的体調推定時間における前記人の基礎的な体調を推定する手段であることを特徴とする。
【００１４】
前記基礎的体調推定手段は、前記生体調節機能要素判定手段により求められた前記各生体調節機能要素に関する各時系列変化のうち、優先順位の高い前記生体調節機能要素の時系列変化が、所定の基準を満たす場合に、前記所定の基礎的体調推定時間における前記人の基礎的な体調を、所定のレベルと推定する第１基礎的体調推定手段と、前記第１基礎的体調推定手段における推定対象とならない場合に、前記第１基礎的体調推定手段において用いた前記生体調節機能要素よりも優先順位の低い他の生体調節機能要素に関する前記時系列変化を用い、前記人の基礎的な体調を、所定の基準に基づいて分類される所定のレベルと推定する第２基礎的体調推定手段とを有し、前記基礎的体調出力手段が、前記第１基礎的体調推定手段又は前記第２基礎的体調推定手段によって推定される前記人の基礎的な体調のレベルを出力する構成とすることが好ましい。
【００１５】
前記第１基礎的体調推定手段は、前記生体調節機能要素判定手段において判定された複数の前記生体調節機能要素のうち、前記自律神経機能への関連性の高い指標又は前記肉体・精神疲労への関連性の高い指標の時系列変化を用い、その時系列変化が所定の基準を満たす場合に、前記人の基礎的な体調のレベルを、「良好」又は「不良」と推定する構成とすることが好ましい。
前記第１基礎的体調推定手段は、前記自律神経機能への関連性の高い指標の時系列変化が、所定の基準を満たす場合に、前記人の基礎的な体調を「不良」と推定し、前記肉体・精神疲労への関連性の高い指標の時系列変化が、所定の基準を満たす場合に、前記人の基礎的な体調を「良好」と推定する手段であることが好ましい。
【００１６】
前記第２基礎的体調推定手段は、前記第１基礎的体調推定手段において「良好」又は「不良」と推定されない場合に、前記生体調節機能要素判定手段において判定された複数の前記生体調節機能要素のうち、前記感覚への関連性の高い指標の時系列変化から、前記人の基礎的な体調のレベルを、「良好」、「不良」又はそれらの「中間状態」のいずれかと推定する構成であることが好ましい。
【００１７】
本発明の生体状態推定方法は、コンピュータを用いて、生体信号測定装置により測定した人の生体信号を分析し、生体状態を推定する生体状態推定方法であって、前記生体信号を分析して、脳機能、自律神経機能、肉体・精神疲労又は感覚との関連性の高いゆらぎに起因する指標を含む、生体調節機能要素の状態の変動を示す複数の指標を、それぞれ予め設定された所定の判定時間毎に算出してその時系列変化を求める生体調節機能要素判定手順と、前記生体調節機能要素判定手順における前記各生体調節機能要素の各判定時間よりも長く設定される基礎的体調推定時間に対応した、前記人の基礎的な体調を推定する基礎的体調推定手順と、前記基礎的体調推定手順によって推定される前記人の基礎的な体調のレベルを、前記所定の基礎的体調推定時間の経過毎に出力する基礎的体調出力手順とを有し、前記基礎的体調推定手順が、前記生体調節機能要素判定手順によって求められる前記各生体調節機能要素の状態に関する各時系列変化を、予め設定された優先順位に従って分析し、所定の基準に照らして、前記所定の基礎的体調推定時間における前記人の基礎的な体調を推定する手順であることを特徴とする。
【００１８】
前記生体状態推定方法の前記基礎的体調推定手順は、前記生体調節機能要素判定手順により求められた前記各生体調節機能要素に関する各時系列変化のうち、優先順位の高い前記生体調節機能要素の時系列変化が、所定の基準を満たす場合に、前記所定の基礎的体調推定時間における前記人の基礎的な体調を、所定のレベルと推定する第１基礎的体調推定手順と、前記第１基礎的体調推定手順における推定対象とならない場合に、前記第１基礎的体調推定手順において用いた前記生体調節機能要素よりも優先順位の低い他の生体調節機能要素に関する前記時系列変化を用い、前記人の基礎的な体調を、所定の基準に基づいて分類される所定のレベルと推定する第２基礎的体調推定手順とを有し、前記基礎的体調出力手順が、前記第１基礎的体調推定手順又は前記第２基礎的体調推定手順によって推定される前記人の基礎的な体調のレベルを出力する構成であることが好ましい。
【００１９】
本発明のコンピュータプログラムは、生体状態推定装置としてのコンピュータに、生体信号測定装置により測定した人の生体信号を分析し、生体状態を推定する生体状態推定手順を実行させるコンピュータプログラムであって、前記生体状態推定手順として、前記生体信号を分析して、脳機能、自律神経機能、肉体・精神疲労又は感覚との関連性の高いゆらぎに起因する指標を含む、生体調節機能要素の状態の変動を示す複数の指標を、それぞれ予め設定された所定の判定時間毎に算出してその時系列変化を求める生体調節機能要素判定手順と、前記生体調節機能要素判定手順における前記各生体調節機能要素の各判定時間よりも長く設定される基礎的体調推定時間に対応した、前記人の基礎的な体調を推定する基礎的体調推定手順と、前記基礎的体調推定手順によって推定される前記人の基礎的な体調のレベルを、前記所定の基礎的体調推定時間の経過毎に出力する基礎的体調出力手順とを実行させ、前記基礎的体調推定手順が、前記生体調節機能要素判定手順によって求められる前記各生体調節機能要素の状態に関する各時系列変化を、予め設定された優先順位に従って分析し、所定の基準に照らして、前記所定の基礎的体調推定時間における前記人の基礎的な体調を推定する手順であることを特徴とする。
【００２０】
前記コンピュータプログラムの前記基礎的体調推定手順は、前記生体調節機能要素判定手順により求められた前記各生体調節機能要素に関する各時系列変化のうち、優先順位の高い前記生体調節機能要素の時系列変化が、所定の基準を満たす場合に、前記所定の基礎的体調推定時間における前記人の基礎的な体調を、所定のレベルと推定する第１基礎的体調推定手順と、前記第１基礎的体調推定手順における推定対象とならない場合に、前記第１基礎的体調推定手順において用いた前記生体調節機能要素よりも優先順位の低い他の生体調節機能要素に関する前記時系列変化を用い、前記人の基礎的な体調を、所定の基準に基づいて分類される所定のレベルと推定する第２基礎的体調推定手順とを実行し、前記基礎的体調出力手順が、前記第１基礎的体調推定手順又は前記第２基礎的体調推定手順によって判定される前記人の基礎的な体調のレベルを出力する構成であることが好ましい。
【００２１】
前記コンピュータプログラムの前記第１基礎的体調推定手順は、前記生体調節機能要素判定手順において判定された複数の前記生体調節機能要素のうち、前記自律神経機能への関連性の高い指標又は前記肉体・精神疲労への関連性の高い指標の時系列変化を用い、その時系列変化が所定の基準を満たす場合に、前記人の基礎的な体調のレベルを、「良好」又は「不良」と推定する構成であることが好ましい。
前記コンピュータプログラムの前記第２基礎的体調推定手順は、前記第１基礎的体調推定手順において「良好」又は「不良」と推定されない場合に、前記生体調節機能要素判定手順において判定された複数の前記生体調節機能要素のうち、前記感覚への関連性の高い指標の時系列変化から、前記人の基礎的な体調のレベルを、「良好」、「不良」又はそれらの「中間状態」のいずれかと推定する構成であることが好ましい。
【発明の効果】
【００２２】
本発明によれば、脳機能、自律神経機能、肉体・精神疲労又は感覚との関連性の高いゆらぎに起因する指標を含む、生体調節機能要素の状態の変動を示す複数の指標を、それぞれ予め設定された所定の判定時間毎に算出し、その時系列変化を求め、各時系列変化を優先順位に従って分析することで、各生体調節機能要素の判定時間よりも長く設定される基礎的体調推定時間に対応した人の調節システムの基本となるゆらぎ性能から体調を求めることができる。すなわち、所定の基礎的体調推定時間内において、二律背反性のある複数の生体調節機能要素毎に求められる時系列変化に基づき、基本的な体調を推定するため、各ゆらぎ性能毎に優先順位を統計的に求め、その優先順位に従った生体調節機能要素の指標を用いて、基礎的体調推定時間における基礎的な体調のレベル、例えば「良好」、「不良」を推定し、さらに、この推定に使用した生体調節機能によって基礎的な体調の分類ができない場合に、優先順位の低い他の生体調節機能に関する時系列変化を用いて、調節システムの外乱に対応する性能に対応する基礎的な体調のレベル、例えば、「良好」、「不良」、それらの「中間状態」などを推定して出力する。
【００２３】
従って、本発明では、体調レベルに応じて関与する生体の調節機能要素（生体調節機能要素）の判定結果が別々に出力されるだけでなく、各生体調節機能要素の判定結果を利用して、各生体調節機能要素の判定結果とは別に基礎的な体調をゆらぎの持つ恒常性維持機能から推定して出力しており、運転者や管理者等は、過去３０分や１時間等における基礎的な体調を一見して把握することができる。そのため、この情報が積み重なることで、運転者等は自らが休憩をとった方がよいなどの自己判断を行うことが可能となる。もちろん、従来の個別の生体調節機能要素の判定結果に基づいた警告と併用ないしはそれらの警告を優先することが好ましいが、基礎的な体調を定期的に出力することで、休憩判断その他の必要な行為（管理者からの警告など）をより適切に行うことができる。
【図面の簡単な説明】
【００２４】
【図１】図１（ａ）は、本発明の一の実施形態において用いた背部体表脈波を測定する生体信号測定装置の一例を示した分解図であり、図１（ｂ）は、その要部断面図である。
【図２】図２は、本発明の一の実施形態に係る生体状態推定装置の構成を模式的に示した図である。
【図３】図３（ａ）は、周波数傾き時系列波形演算手段により求められるゼロクロス法、ピーク検出法による周波数傾き時系列波形の一例を示した図であり、図３（ｂ）は、振幅の変動を判別するために図３（ａ）をスムージング処理した波形である。
【図４】図４は、分布率演算手段により求められる分布率の時系列波形の一例を示した図である。
【図５】図５は、疲労曲線演算手段により求められる疲労曲線の一例を示した図である。
【図６】図６は、恒常性維持機能レベル演算手段により求められる恒常性維持機能レベルの時系列波形の一例を示した図である。
【図７】図７（ａ）は、体調マップ演算手段により求められる体調マップの一例を示した図であり、図７（ｂ）は、感覚マップ演算手段により求められる感覚マップの一例を示した図である。
【図８】図８は、基礎的体調推定手段による基礎的体調の推定ステップを説明するためのフローチャートである。
【図９】図９（ａ），（ｂ）は、実験例２における被験者Ａ及びＢのある月における実験結果を示した図である。
【図１０】図１０は、図９（ｂ）に示した被験者Ｂの１４日と１５日における生体状態推定装置による体調推定結果の時系列データと、その前後の運行日における運行開始時刻を示した図である。
【図１１】図１１は、実験例２の全被験者の主観評価と生体状態推定装置による体調推定の「良」の割合との相関を正規化した分布図である。
【発明を実施するための形態】
【００２５】
以下、図面に示した本発明の実施形態に基づき、本発明をさらに詳細に説明する。本発明において採取する生体信号は、例えば、指尖容積脈波、背部から採取される音・振動情報（以下、「背部音・振動情報」）等が挙げられるが、好ましくは、背部音・振動情報である。背部音・振動情報には、上記のように、人の上体背部から検出される心臓と大動脈の運動から生じる音・振動情報であり、心室の収縮期及び拡張期の情報と、血液循環の補助ポンプとなる血管壁の弾性情報及び血圧による弾性情報並びに反射波の情報、すなわち、背部体表脈波（ＡＰＷ）や疑似心音情報を含んでいる。また、心拍変動に伴う信号波形は交感神経系及び副交感神経系の神経活動情報（交感神経の代償作用を含んだ副交感神経系の活動情報）を含み、大動脈の揺動に伴う信号波形には交感神経活動の情報や内分泌系の情報を含んでいるため、異なる観点から生体調節機能要素を判定するのに適している。
【００２６】
生体信号を採取するための生体信号測定装置は、指尖容積脈波であれば指尖容積脈波計を用いることができる。背部音・振動情報であれば、例えば、圧力センサを用いることも可能であるが、好ましくは、（株）デルタツーリング製の居眠り運転警告装置（スリープバスター（登録商標））で使用されている生体信号測定装置１を用いる。図１は生体信号測定装置１の概略構成を示したものである。この生体信号測定装置１は、乗物の運転席に組み込んで使用することができ、手指を拘束することなく生体信号を採取できる。
【００２７】
生体信号測定装置１を簡単に説明すると、図１（ａ），（ｂ）に示したように、上層側から順に、第一層１１、第二層１２及び第三層１３が積層された三層構造からなり、三次元立体編物等からなる第一層１１を生体信号の検出対象である人体側に位置させて用いられる。従って、人体の体幹背部からの生体信号、特に、心室、心房、大血管の振動に伴って発生する生体音（体幹直接音ないしは生体音響信号）を含む心臓・血管系の音・振動情報（背部体表脈波（ＡＰＷを含む））は、生体信号入力系である第一層１１にまず伝播される。第二層１２は、第一層１１から伝播される生体信号、特に心臓・血管系の音・振動を共鳴現象又はうなり現象によって強調させる共鳴層として機能し、ビーズ発泡体等からなる筐体１２１、固有振動子の機能を果たす三次元立体編物１２２、膜振動を生じるフィルム１２３を有して構成される。第二層１２内において、マイクロフォンセンサ１４が配設され、音・振動情報を検出する。第三層１３は、第二層１２を介して第一層１１の反対側に積層され、外部からの音・振動入力を低減する。
【００２８】
次に、本実施形態の生体状態推定装置１００の構成について図２に基づいて説明する。生体状態推定装置１００は、生体調節機能要素判定手段２００、基礎的体調推定手段３００及び基礎的体調出力手段４００を有して構成されている。生体状態推定装置１００は、コンピュータ（マイクロコンピュータ等も含む）から構成され、生体調節機能要素判定手段２００、基礎的体調推定手段３００及び基礎的体調出力手段４００として機能する、生体状態推定手順を実施する生体調節機能要素判定手順、基礎的体調推定手順及び基礎的体調出力手順をコンピュータに実行させるコンピュータプログラムが記憶部に設定されている。
【００２９】
なお、コンピュータプログラムは、記録媒体に記憶させてもよい。この記録媒体を用いれば、例えば上記コンピュータに上記プログラムをインストールすることができる。ここで、上記プログラムを記憶した記録媒体は、非一過性の記録媒体であっても良い。非一過性の記録媒体は特に限定されないが、例えば フレキシブルディスク、ハードディスク、ＣＤ−ＲＯＭ、ＭＯ（光磁気ディスク）、ＤＶＤ−ＲＯＭ、メモリカードなどの記録媒体が挙げられる。また、通信回線を通じて上記プログラムを上記コンピュータに伝送してインストールすることも可能である。
【００３０】
生体調節機能要素判定手段２００は、本実施形態では、上記の生体信号測定装置１により測定された生体信号である背部音・振動情報を分析し、人の基礎的な体調の推定に用いる生体調節機能要素に関する複数種類の指標を、それぞれ予め設定された所定の判定時間毎に算出してその時系列変化を求める。なお、分析対象の生体信号は、指尖容積脈波等であってもよく、本発明がそれを排除するものではないが、背部音・振動情報が好ましいことは上記のとおりである。
【００３１】
生体調節機能要素判定手段２００において判定される複数種類の生体調節機能要素は限定されるものではないが、少なくとも、脳機能や自律神経機能への関連性の高い指標、肉体・精神疲労への関連性の高い指標、及び、感覚への関連性の高い指標を含むものであることが好ましい。これらは、人の恒常性維持機能に影響を与える脳波のゆらぎの変動の様子、あるいは、体温調節機能に代表される生体調節機能が仕事をする様子を示す指標であり、体調により各調節機能に与える影響が大きいためである。
【００３２】
脳機能や自律神経機能への関連性の高い指標としては、例えば、採取した生体信号を処理して得られる周波数傾きの時系列波形、上記従来技術の項で説明した３つの信号の分布率の時系列波形、疲労曲線、恒常性維持機能レベルの判定の時系列の変動が挙げられる。周波数傾きの時系列波形は、恒常性維持機能の調節作用のベースにあるものはゆらぎを示すものであり、そのゆらぎは二律背反性のある機能のバランスをうまく調整し、人の自律神経機能との関連性を特に高く示している。これは統計的な手法による裏付けがなされていることである。分布率から求められた各周波数帯域の時系列波形は、ゆらぎのリズムに間接的に関与する脳波の種類（θ波、α波、β波）に対応し、人の脳機能及び自律神経機能に加え、内分泌系の調節機能との関連性を高く示している。恒常性は、内分泌系、自律神経系など様々な調節システムによって保たれるため、そのレベルの変動は、脳、自律神経系及び内分泌系のゆらぎによる調節性能とも深く関連している。
【００３３】
生体調節機能要素判定手段２００は、上記の脳機能、自律神経機能及び内分泌系のゆらぎの変動の仕方を捉える指標を求める演算手段として、周波数傾きの時系列波形を求める周波数傾き時系列波形演算手段２１０、分布率を求める分布率演算手段２２０、疲労曲線を求める疲労曲線演算手段２３０、及び恒常性維持機能レベルを求める恒常性維持機能レベル演算手段２４０とを有している。
【００３４】
周波数傾き時系列波形演算手段２１０は、生体信号測定装置１のセンサ１４から得られる背部音・振動情報をフィルタリング処理した１Ｈｚ近傍の背部体表脈波（ＡＰＷ）の時系列波形から周波数の時系列波形を求めた後、周波数の時系列波形をスライド計算して周波数傾き時系列波形を求める（図３（ａ），（ｂ）参照）。周波数傾き時系列波形演算手段２１０は、本発明者らによる上記特許文献１等に開示されているように、背部体表脈波（ＡＰＷ）の時系列波形において、正から負に切り替わる点（ゼロクロス点）を用いる手法（ゼロクロス法）と、背部体表脈波（ＡＰＷ）の時系列波形を平滑化微分して極大値（ピーク点）を用いて時系列波形を求める方法（ピーク検出法）の２つの方法がある。
【００３５】
ゼロクロス法では、ゼロクロス点を求めたならば、それを例えば５秒毎に切り分け、その５秒間に含まれる時系列波形のゼロクロス点間の時間間隔の逆数を個別周波数ｆとして求め、その５秒間における個別周波数ｆの平均値を当該５秒間の周波数Ｆの値として採用する。そして、この５秒毎に得られる周波数Ｆを時系列にプロットすることにより、周波数の変動の時系列波形を求める。
【００３６】
ピーク検出法では、背部体表脈波（ＡＰＷ）の上記時系列波形を、例えば、ＳａｖｉｔｚｋｙとＧｏｌａｙによる平滑化微分法により極大値を求める。次に、例えば５秒ごとに極大値を切り分け、その５秒間に含まれる時系列波形の極大値間の時間間隔の逆数を個別周波数ｆとして求め、その５秒間における個別周波数ｆの平均値を当該５秒間の周波数Ｆの値として採用する。そして、この５秒毎に得られる周波数Ｆを時系列にプロットすることにより、周波数の変動の時系列波形を求める。
【００３７】
周波数傾き時系列波形演算手段２１０は、ゼロクロス法又はピーク検出法により求められた周波数の変動の時系列波形から、所定のオーバーラップ時間（例えば１８秒）で所定の時間幅（例えば１８０秒）の時間窓を設定し、時間窓毎に最小二乗法により周波数の傾きを求め、その傾きの時系列波形を出力する。このスライド計算を順次繰り返し、ＡＰＷの周波数の傾きの時系列変化を周波数傾き時系列波形として出力する。
【００３８】
背部体表脈波（ＡＰＷ）は、中枢系である心臓の制御の様子を主として含む生体信号、すなわち、動脈の交感神経支配の様子、並びに、交感神経系と副交感神経系の出現情報を含む生体信号であり、ゼロクロス法により求めた周波数傾き時系列波形（図３（ａ），（ｂ）において「０ｘ」と表示した波形）は、心臓の制御の状態により関連しており、交感神経の出現状態を反映しているが、ピーク検出法により求めた周波数傾き時系列波形（図３（ａ），（ｂ）において「Ｐｅａｋ」と表示した波形）は、心拍変動により関連している。従って、自律神経機能の状態をより明確に把握するためには、ゼロクロス法を用いて求めた周波数傾き時系列波形を用いることが好ましい。
【００３９】
交感神経の活動は、血管弾性や血管径に影響を与え、さらに、血管壁からの反射波の影響が、人の背部から検出される音・振動情報に含まれる疑似心音情報（背部から検出されるため、心臓から背部表面までの間の筋肉、皮膚等により２０Ｈｚ近傍の信号として検出される）の疑似Ｉ音（心音Ｉ音に相当）と疑似ＩＩ音（心音ＩＩ音に相当）の間の波形成分に重畳される。これが、ゼロクロス法におけるゼロクロス点間の幅と、ピーク検出法におけるピーク点間の幅とを異ならせる理由であり、ゼロクロス法では反射波の影響を受けた周期となっている。よって、ゼロクロス法による周波数傾き時系列波形を見ることで交感神経の情報を捉えることができる。
一方、ピーク値は上記のように心拍変動の情報を反映しているが、心拍変動は主に副交感神経によって制御されている。そのため、ピーク値を見ると副交感神経の情報を捉えることができる。
【００４０】
周波数傾き時系列波形演算手段２１０により得られるゼロクロス法による周波数傾き時系列波形は、睡眠前の所定のタイミングで眠気に対する抵抗として生じる交感神経活動の一時的亢進に伴って振幅が拡大し、長周期化する傾向を示した場合に、入眠予兆現象の指標と捉えられることが知られている（特許文献４参照）。また、入眠予兆現象を示す波形が出現した後、波形が収束傾向を示し、その後、より長周期の大きな変動ゆらぎを示すと、その長周期のゆらぎを示し始めたポイントが、入眠直前の切迫睡眠現象を示す指標と捉えられることが知られている。
【００４１】
分布率演算手段２２０は、まず、周波数傾き時系列波形演算手段２１０から得られる周波数傾き時系列波形をそれぞれ周波数分析して、心循環系のゆらぎの特性が切り替わる周波数である上記の０．００３３Ｈｚよりも低い周波数の機能調整信号、機能調整信号よりも高い周波数の疲労受容信号、及び疲労受容信号よりも高い周波数の活動調整信号に相当するＵＬＦ帯域からＶＬＦ帯域に属する各周波数成分を抜き出す。次に、これらの周波数成分のそれぞれの分布率を時系列に求める。すなわち、３つの周波数成分のパワースペクトルの値の合計を１とした際の各周波数成分の割合を分布率として時系列に求める（図４参照）。
【００４２】
本実施形態では、図４に示したように、機能調整信号として０．００１７Ｈｚの周波数成分を用い、疲労受容信号として０．００３５Ｈｚの周波数成分を用い、活動調整信号として０．００５３Ｈｚの周波数成分を用いている。心疾患の一つである心房細動において、心・循環系のゆらぎの特性が切り替わる周波数は、０．００３３Ｈｚと言われており、０．００３３Ｈｚ近傍のゆらぎの変化を捉えることで、自律神経の活動、恒常性維持に関する情報が得られる。また、０．００３３Ｈｚ近傍以下と０．００５３Ｈｚ近傍の周波数帯は、主に体温調節に関連するもので、０．０１〜０．０４Ｈｚの周波数帯は自律神経の制御に関連するものと言われている。そして、本発明者らが実際に、生体信号に内在するこれら低周波のゆらぎを算出する周波数傾き時系列波形を求め、それを周波数解析したところ、０．００３３Ｈｚよりも低周波の０．００１７Ｈｚ、０．００３３Ｈｚ近傍の０．００３５Ｈｚを中心とする周波数帯のゆらぎと、さらにこれら２つ以外に、０．００５３Ｈｚを中心とする周波数帯のゆらぎがあることが確認できた。但し、各信号の周波数成分は個人差等により調整することも可能であり、機能調整信号は０．００３３Ｈｚ未満の範囲で好ましくは０．００１〜０．００２７Ｈｚの範囲で、疲労受容信号は０．００２〜０．００５２Ｈｚの範囲で、活動調整信号は０．００４〜０．００７Ｈｚの範囲で調整して用いることができる。
【００４３】
分布率演算手段２２０により求められる分布率の時系列変化は、特許文献２に示されているように、例えば、０．００１７Ｈｚの分布率が急低下し、かつ０．００５３Ｈｚの分布率が急上昇する変化を示す時点を切迫睡眠現象の出現時点と捉えることができる。
【００４４】
疲労曲線演算手段２３０は、本発明者らの特開２００９−２２６１０号公報に開示されている手段であり、ゼロクロス法による求めた周波数傾き時系列波形を絶対値処理して積分値を算出し、この積分値を疲労度として所定の判定時間毎に求めて、時間に対応してプロットし、図５に示したような疲労曲線を求める手段である。筋活動は、筋肉の収縮又は弛緩であり、交感神経の情報を反映しているゼロクロス法による周波数傾き時系列波形の積分情報である疲労曲線は筋活動との相関性が高い（非特許文献１参照）。よって、疲労曲線では、その傾きが所定以上変動するポイントが特異点を示しており、各特異点は、増大する疲労に対応して、筋活動が生じたことを示すポイントや血流量が増大したポイントを示している。
【００４５】
恒常性維持機能レベル判定手段２４０は、特許文献３に開示の技術に基づくものであり、周波数傾き時系列波形演算手段２１０により得られるゼロクロス法を用いた各周波数傾き時系列波形の微分波形の正負、周波数傾き時系列波形を積分した積分波形の正負、ゼロクロス法を利用した周波数傾き時系列波形とピーク検出法を利用した周波数傾き時系列波形をそれぞれ絶対値処理して得られた各周波数傾き時系列波形の絶対値等のうち、少なくとも１つ以上を用いて判定する。これらの組み合わせにより、恒常性維持機能のレベルがいずれに該当するかを求める。例えば、周波数傾きと積分値を用いて、所定以上の場合に「恒常性維持機能レベル１」と判定し、あるいは、微分値が所定位置以下であって、かつ、２つの絶対値のうちの「ピーク優位」の場合に「恒常性維持機能レベル４」と判定するように設定できる。そして、例えば、上記の条件を様々に組み合わせ、人の状態との相関をとり、レベル１〜３と判定される場合を、普通から良好な状態、レベル４〜６と判定される場合を、注意の必要な状態と判定する。また、入眠予兆や切迫睡眠の兆候が生じているなどと判定された場合には、直ちに警告を要するレベルとして、それぞれの状態によりレベル７〜１１といった指標を付与する。株式会社デルタツーリング製、商品名「スリープバスター」では、恒常性維持機能レベル判定手段２４０による判定結果が、例えば、図６に示したように表示されるように設定されている。
【００４６】
肉体・精神疲労への関連性の高い指標としては、特許文献２に開示された指標である体調マップ及び感覚マップを用いることができる。これらは、ゆらぎの変動の仕方をグラフ化したもので、人の肉体・精神疲労との関連性を高く示している。
【００４７】
そのため、本実施形態の生体調節機能要素判定手段２００は、さらに体調マップ演算手段２５０及び感覚マップ演算手段２６０を有している。生体信号測定装置１から取得した背部音・振動情報から得られる背部体表脈波（ＡＰＷ）を周波数分析し、対象となる解析区間について、解析波形を両対数軸表示に表し、その解析波形を低周波帯域、中周波帯域、高周波帯域に分け、区分けした解析波形の傾きと、全体の解析波形の形とから一定の基準に基づいて解析波形の点数化を行い、それを座標軸にプロットしたものである。体調マップは、自律神経系の制御の様子を交感神経と副交感神経のバランスとして見たものであり、感覚マップは、体調マップに心拍変動の変化の様子を重畳させたものである。
【００４８】
具体的には、体調マップ演算手段２５０は、背部体表脈波を周波数解析した解析波形について、所定周期領域毎に回帰直線をまず求める。次に、周期領域毎に求められる各回帰直線を、その傾きに基づいて領域得点を付与すると共に、隣接する周波数領域における回帰直線間のパワースペクトル密度の値の較差及び回帰直線間の傾きの違いに基づき、各回帰直線全体における分岐現象を示す折れ点数を求め、その折れ点数に基づいた形状得点を付与し、領域得点及び形状得点の少なくとも一方を用いて、各解析波形についての判定基準点を求める。領域得点としては、各領域における各回帰直線の傾きを略水平状態、上向き及び下向きの３つに分け、例えば略水平状態の得点を基準として、上向きの場合と下向きの場合とで得点を増減させて得点を付与する。形状得点としては、折れ点数が少ないほど高得点を付与する。
【００４９】
判定基準点を求める際には、ゼロクロス法により求めた周波数傾き時系列波形を用いて第１の判定基準点を求め、ピーク検出法により求めた周波数傾き時系列波形を用いて第２の判定基準点を求める。そして、第１の判定基準点に基づく指標を一方の軸に、第２の判定基準点に基づく指標を他方の軸にとって、座標点をプロットし、図７（ａ）に例示したような体調マップが作成される。体調マップでは、座標点同士を結んだ座標時系列変化線が、１／ｆの傾きに近似した変化傾向であると判定された場合には快適と判定され、上下方向に変化していると判定された場合には不快と判定される。図７（ａ）は、座標原点に合わせずに複数の座標点を結んでいるが、時間的に異なる２点の変化傾向を見る場合、１点目を座標原点に合わせて、２点目が第４象限にプロットされると、この生体調節機能要素に関しては「良好」ということになり、判断がより容易になる。
【００５０】
感覚マップ演算手段２６０は、心拍変動に関連するピーク検出法を用いた周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に周波数の平均値を求める移動計算を行い、時間窓毎に得られる周波数の平均値の時系列変化を周波数変動時系列波形として求め、さらに、ゼロクロス法を用いた周波数の時系列波形から求められる機能点に対応する指標を一方の軸にとると共に、ピーク検出法により求められる上記の周波数変動時系列波形の所定の時間幅における変化量に対応する指標を他方の軸にとり、機能点と変化量とから求められる座標の時系列変化を求めていく手段である。図７（ｂ）がこのようにして求めた感覚マップの一例である。図７（ｂ）では、座標原点に合わせずに複数の座標点を結んだものであるが、時間的に異なる２点の変化傾向を見る場合、１点目を座標原点に合わせて、２点目をプロットすると、両者間の離隔距離及び離隔方向が判断しやすくなる。
【００５１】
なお、機能点は、比較対象の前後２つの時間範囲における解析波形の判定基準点間において、次式：
機能点＝後時間範囲の判定基準点＋（後時間範囲の判定基準点−前時間範囲の判定基準点）×ｎ、（但し、ｎは補正係数）
により求められる。
【００５２】
感覚への関連性の高い指標としては、上記の恒常性維持機能レベル判定手段２４０により求められる恒常性維持機能レベルの時系列変化のうち、例えば、周波数傾きと積分値を用いて、普通から良好といえるレベルの指標（上記の例では、レベル１〜３）、注意を要するレベルの指標（上記の例では、レベル４〜６）を用いてそれらの頻出頻度を用いて判定できる。恒常性維持機能レベルは、上記のように自律神経機能の状態と高く関連しているが、体調、基礎的な体力、あるいは動機付けにより、疲労に対して交感神経代償作用が発現した際、疲労感を感じるときと感じないときがある。従って、疲労に対する交感代償作用と基礎的な体調は、それを疲労として感じる感覚との関連性が高い。なお、ここでいう感覚とは、倦怠感あるいは覚低状態を伴う喪失感に似た感覚のことである。
【００５３】
基礎的体調推定手段３００は、上記の生体調節機能要素判定手段２００において求められる各生体調節機能要素のゆらぎ性能に関する各時系列変化から、所定の基準に照らして分析対象の人の基礎的な体調（基礎的体調）を推定する手段である。生体調節機能要素判定手段２００においては、上記のように、生体調節機能要素のゆらぎ性能に関する時系列変化が複数種類得られるように設定されているが、この複数種類得られる各時系列変化は、所定の判定時間毎に得られる。例えば、周波数傾き時系列波形演算手段２１０は、生体信号測定装置１からのデータを取得した後、最初の演算結果が出力されるまで数分かかるが、その後は、例えば、１８秒ごとに得られ、それにより時系列変化が求められる。分布率演算手段２２０により得られる分布率、疲労曲線演算手段２３０により得られる疲労度、及び恒常性維持機能レベル演算手段２４０により得られる恒常性維持機能レベルも最初の演算結果が出力されるまで数分かかり、その後、例えば１８秒毎に得られ、それぞれ時系列変化が求められる。体調マップ演算手段２５０及び感覚マップ演算手段２６０によりそれぞれ得られる演算結果は、最初は２０〜３０分かかるが、２点目はその約十数分後、３点目以降は数分毎に得られる。これに対し、基礎的体調推定手段３００は、各生体調節機能要素におけるこれらの各判定時間よりも長い時間（基礎的体調推定時間）について、基礎的体調を推定する。
【００５４】
例えば、運転者が体調を把握する場合、各生体調節機能要素のゆらぎ性能の判定結果を個別にモニタに表示させることでももちろん可能であるが、この場合には、データを分析する運転者の分析能力、判定能力によりバラツキが生じやすい。従って、機械的に判定することによってバラツキを吸収し、それらの組み合わせを必要に応じて寄与率を考慮するなどして統計処理することにより、基礎的な体調を推定することが好ましい。すなわち、比較的短時間毎に出力されるデータを統計的に組み合わせることにより、現在の生体の基礎的な制御能（体調）を推定する方法である。各生体調節機能要素の判定結果は、逐次変化する運転者の最新の心身状態を把握するのに適しているが、運転という作業を行っている最中にこれらの判定結果を逐次出力されても、運転者自身がその状態を分析して基礎的な体調を自覚することは難しい場合がある。そこで、基礎的体調推定手段３００は、より長い時間に設定された基礎的体調推定時間における基礎的体調を推定するようにしたものである。基礎的体調推定時間は、１５分、３０分、６０分等、任意に設定できるが、頻度が多すぎ、状態の変化が多すぎると結果的に運転者が体調を把握しにくいことにつながるため、２０〜４０分程度の時間とすることが好ましい。
【００５５】
基礎的体調推定手段３００は、各生体調節機能要素の判定結果を用いて基礎的体調を判定する。具体的には、所定の長い時間に設定された基礎的体調推定時間における体調をより反映している生体調節機能要素の判定結果を優先して用いる。そのため、本実施形態では、基礎的体調推定手段３００を構成させるコンピュータプログラムである基礎的体調推定手順に、各生体調節機能要素の判定結果を利用する優先順位が設定されている。基礎的体調推定手段３００は、より詳細には、図２に示したように、第１基礎的体調推定手段３１０と第２基礎的体調推定手段３２０とを有し、図８のフローチャートに示したように、まず、第１基礎的体調推定手段３１０において基礎的体調の推定を行い（Ｓ１１０，Ｓ１２０）、次に、第２基礎的体調推定手段３２０における基礎的体調の推定（Ｓ１３０）を行うように設定されている。
【００５６】
第１基礎的体調推定手段３１０は、優先順位の設定された生体調節機能要素のうち、優先順位の高い生体調節機能要素の時系列変化が、所定の基準を満たすか否かを求め、所定のレベルを満たす場合に、基礎的体調を当該所定のレベルと推定する。第２基礎的体調推定手段３２０は、第１基礎的体調推定手段において判定対象とならない場合、すなわち、第１基礎的体調推定手段３１０における所定の基準を満たさない場合に、第１基礎的体調推定手段３１０において用いられた生体調節機能要素よりも優先順位の低い他の生体調節機能要素に関する時系列変化を用い、所定の基準に基づいて基礎的体調を推定する。従って、第２基礎的体調推定手段３２０が基礎的体調の推定の最終手段となるため、第２基礎的体調推定手段３２０は、判定対象のデータをいずれかの基礎的体調に必ず分類するように設定している。
【００５７】
基礎的体調推定手段３００において推定対象となる生体調節機能要素判定手段２００によって得られる各生体調節機能要素の優先順位としては、後述の実験例のように多数の事例を分析して統計的に、必要に応じてそれらに寄与率も考慮して決定することが好ましい。それにより、より推定精度を向上できる。本実施形態においては、後述の実験例における統計処理に基づいて、自律神経機能への関連性の高い指標又は肉体・精神疲労への関連性の高い指標の時系列変化のいずれかの優先順位を高く設定し、次の優先順位に、感覚への関連性の高い指標の時系列変化となるように設定している。従って、第１基礎的体調推定手段３１０においては、自律神経機能への関連性の高い指標の時系列変化が、所定の条件を満たすか否か、肉体・精神疲労への関連性の高い指標の時系列変化が、所定の条件を満たすか否かを判別するように設定している。また、自律神経機能への関連性の高い指標の時系列変化が、所定の条件を満たす場合に、基礎的体調を「不良」と推定し、肉体・精神疲労への関連性の高い指標の時系列変化が、所定の条件を満たす場合に、基礎的体調を「良好」と推定するように設定している。
【００５８】
脳機能、自律神経機能あるいは内分泌系の調節機能のゆらぎに基づいた恒常性維持機能に関する指標である上記の周波数傾き時系列波形、分布率、疲労曲線（疲労度）、及び恒常性維持機能レベルは、入眠予兆、切迫睡眠、覚低走行状態など、疲労の蓄積の結果生じる兆候を判別しやすい指標である。また、脳機能によって調節されている恒常性維持機能のゆらぎは、その周波数帯域の差により、内分泌系など、支配される調節システムを異にするが、上記の中でも分布率は、これらの調節システムの急変時、減衰時、増大時がよく反映される指標である。そこで、これらを用いて運転等の作業を行う上で注意や警告を要する状況が基礎的体調推定時間において顕著に現れたり、その頻度が所定以上だったりする場合に、基礎的体調を「不良」と推定する。肉体・精神疲労への関連性の高い指標である体調マップ・感覚マップは、快調、快適に感じている場合の指標を顕著に判別しやすい。そこで、この指標が快調、快適を示す条件の場合に、基礎的体調を「良好」と推定する。第１基礎的体調推定手段３１０における基礎的体調の推定において、脳機能、自律神経機能及び内分泌系の調節機能への関連性の高い指標と、肉体・精神疲労への関連性の高い指標とのいずれを優先して用いるかは制限されるものではないが、肉体・精神疲労の変調も自律神経との関わりが基本的に存在するため、本実施形態のように、脳機能や自律神経機能等への関連性の高い指標を用いた推定を実施し、次に、肉体・精神疲労への関連性の高い指標を用いた推定を実施することが好ましい。
【００５９】
脳機能、自律神経機能及び内分泌系のホルモン分泌調節機能への関連性の高い指標は、本実施形態では上記のように、周波数傾き時系列波形、分布率、疲労曲線（疲労度）、及び恒常性維持機能レベルの４つある。このうち、一つにおいて、入眠予兆等の兆候を所定以上検出した場合に、基礎的体調を「不良」と推定するように設定することも可能であるが、複数の指標において所定の兆候を検出した場合に、基礎的体調を「不良」と推定することは信頼度を高めるため好ましい。
【００６０】
そこで、本実施形態の第１基礎的体調推定手段３１０は、周波数傾き時系列波形演算手段２１０から求められる周波数傾き時系列波形、分布率演算手段２１０から求められる分布率の時系列波形、疲労曲線演算手段２３０から求められる疲労曲線（疲労度の時系列波形）、及び恒常性維持機能レベル演算手段２４０から求められる恒常性維持機能レベルのうち、３つ以上の指標が所定の基準を満たす場合（図８のＳ１１０で「Ｙｅｓ」と判定された場合）に「不良」と推定するように設定している。
【００６１】
本実施形態において「不良」と推定する所定の基準は、次のように設定している。
（ａ）周波数傾き時系列波形演算手段２１０から求められる指標
ゼロクロス法を用いた周波数傾き時系列波形において、振幅変化を比較し、複数回（通常、２〜４回の範囲で設定）連続で１つ前の振幅の９〜６割未満に変化する収束箇所が生じた場合（交感神経活動が低下し、眠気に抵抗できない状態に陥ったことを推定する指標）
（ｂ）分布率演算手段２２０から求められる指標
ゼロクロス法を用いた周波数傾き時系列波形の分布率の時系列変化において、所定時間の範囲（通常、６０〜１２０秒間の範囲で設定）で、０．００１７Ｈｚの分布率が急減（通常、減少率１５％以上で設定）し、その間に０．００５３Ｈｚの分布率が急増（通常、増加率１５％以上で設定）した場合（入眠予兆現象の出現を推定する指標）
（ｃ）疲労曲線演算手段２３０から求められる指標
所定時間（通常、３〜１０分の範囲で設定）の間における、ピーク検出法を用いた疲労曲線（ピーク検出法を用いた周波数傾き時系列波形の絶対値の積算の時系列波形）の傾きが、ゼロクロス法を用いた疲労曲線（ゼロクロス法を用いた周波数傾き時系列波形の絶対値の積算の時系列波形）の傾きよりも大きく変化する箇所が１箇所以上存在し、かつ、所定時間経過時に、ピーク検出法を用いた疲労曲線が所定の値以上に至った場合（副交感神経活動が極端に優位な状態であることを推定する指標）
（ｄ）恒常性維持機能レベル演算手段２４０から求められる指標
１８秒毎に得られる恒常性維持機能レベルのうち、普通レベルよりは低いレベル、上記の例ではレベル４〜６という注意判定が数回から十数回以上出現する場合（副交感神経活動が優位な状態と推定されるときに出現する指標）、あるいは、警告を要するレベル、上記の例ではレベル７〜１１の警告判定が数回以上出現する場合（交感神経活動の急激な亢進や極端な低下などが推定されるときに出現する指標）
【００６２】
また、本実施形態の第１基礎的体調推定手段３１０は、上記の（ａ）〜（ｄ）の指標のうち３つ以上において「不良」と推定されないデータ（図８のＳ１１０で「Ｎｏ」と判定されたデータ）に関し、体調マップ演算手段２５０及び感覚マップ演算手段２６０の指標を用いて所定の基準を満たす「良好」に相当するか否かを推定する（図８のＳ１２０）。
（ｅ）基礎的体調「良好」と判定される場合の指標
基礎的体調「良好」に相当する基準として、本実施形態では、体調マップ演算手段２５０から求められる時系列変化が、一つ手前の演算結果が出力されるポイント（上記のように、１点目、２点目は所定の時間経過後に出力されるが、３点目以降は数分毎に出力される）を座標原点に合わせた際に、次のポイントが第４象限にプロットされ、かつ、感覚マップ演算手段２６０から求められる時系列変化が、同じく一つ手前のポイントを座標原点に合わせた際に、Ｘ軸方向に所定以上離隔してプロットされる場合に、「良好」と推定するように設定している（図８のＳ１２０で「Ｙｅｓ」の場合）。
【００６３】
なお、基礎的体調が「不良」と推定される（ａ）〜（ｄ）の判定基準及び基礎的体調が「良好」と推定される（ｅ）の判定基準は、後述の多数の事例の統計的分析に基づくものであるが、これに限定されるものではない。例えば、個人毎にデータを蓄積して、個人毎に統計的に条件を設定するようにしてもよい。
【００６４】
第２基礎的体調推定手段３２０は、推定対象のデータが、第１基礎的体調推定手段３１０において「不良」、「良好」の各基準のいずれも満たさない場合（図８のＳ１１０で「Ｎｏ」と判定され、かつ、Ｓ１２０で「Ｎｏ」と判定された場合）に実行される。第２基礎的体調推定手段３２０は、恒常性維持機能レベル判定手段２４０により求められる恒常性維持機能レベルの時系列変化のうち、交感神経活動が優位で普通から良好といえるレベルの指標（上記の例では、レベル１〜３）と、副交感神経活動が優位で注意を要するレベルの指標（上記の例では、レベル４〜６）の境界付近のレベルの出現頻度を比較する。但し、レベルの１段階の違いでは、状態の違いは小さいため、２段階以上違うレベルで比較することが好ましい。本実施形態では、普通から良好といえるレベルの指標のうちの真ん中のレベル２の指標と、注意し始める必要のあるレベル４の指標の出現頻度の割合を比較している。基本的には、交感神経活動が優位で良好状態を示すレベル２の出現頻度が高く、副交感神経活動が優位で注意状態を示すレベル４の出現頻度が低い場合には「良好」と推定でき、出現頻度が逆の関係の場合には「不良」と推定できるが、第２基礎的体調推定手段３２０の分析対象となるデータは、第１基礎的体調推定手段３１０において明確に「良好」、「不良」と推定されなかったものであるため、いずれにも分類しにくいデータも想定される。そこで、本発明では、後述の試験例において、多数の事例を分析し、ベイズ推定の手法により、「良好」、「不良」及びそれらの「中間状態」に分類する基準を設定している。詳細は後述する。
【００６５】
基礎的体調出力手段４００は、本実施形態では、第１基礎的体調推定手段３１０の上記条件を満たす場合（図８のＳ１１０で「Ｙｅｓ」と判定された場合、又は、Ｓ１２０で「Ｙｅｓ」と判定された場合）、その推定結果として基礎的体調のレベルを「不良」（図８のＳ１１１）又は「良好」（図８のＳ１２１）と出力する。第１基礎的体調推定手段３１０の条件を満たさない場合（図８のＳ１１０で「Ｎｏ」と判定され、かつ、Ｓ１２０で「Ｎｏ」と判定された場合）には、第２基礎的体調推定手段３２０の推定結果である「良好」、「不良」それらの「中間状態」を出力する（図８のＳ１３１）。基礎的体調出力手段４００は、基礎的体調の推定結果である「良好」、「不良」又は「中間状態」を、人が認識可能な媒体を通じて出力する。例えば、本実施形態の生体信号測定装置１及び生体状態推定装置１００が、自動車に搭載され、運転者の状態を検知するものであれば、車載モニタに、基礎的体調推定時間の経過毎に推定結果を表示させる。表示方法は、文字でもよいし、図柄等であってもよい。図柄として、運転者が視認しやすいように、例えば、「良好」の場合には、天気の「晴れ」の記号を用い、「不良」の場合には天気の「雨」の記号を用い、「中間状態」の場合には天気の「曇」の記号を用いることができる。そのほか、複数種類のキャラクタで表示してもよいし、これらを併合して表示してもよい。また、車載スピーカを介して音声により出力したりすることもできる。
【００６６】
本実施形態によれば、例えば、運転者の状態を解析する場合、各生体調節機能要素の短時間毎の変化に加え、所定の基礎的体調推定時間毎に、運転者の基礎的な体調を推定して出力することができる。そのため、運転者は、自らが例えばここ３０分間どのような状態で運転しているかを概略的に自覚しやすく、例えば基礎的体調として「不良」が連続して推定されるようであれば、速やかに休憩するなどの決断を促しやすい。また、車載されている生体状態推定装置１００と管理者側のコンピュータとを通信手段を介してつなげておくことにより、運転者の個別の各生体調節機能要素だけでなく、基礎的な体調もリアルタイムで把握でき、管理者側の判断の助けになる。
【００６７】
（実験例１）
次に、基礎的体調推定手段３００を構成する第１基礎的体調推定手段３１０の推定基準、第２基礎的体調推定手段３２０の推定基準の設定に関する実験結果を説明する。
【００６８】
（実験方法）
生体信号測定装置として生体信号測定装置１が運転席のシートバックに装着され、生体信号測定装置１により測定される背部音・振動情報のデータが入力されるコンピュータである生体状態推定装置１００が搭載されている自動車に被験者を着座させて走行させた。なお、生体状態推定装置１００は、株式会社デルタツーリング製、商品名「スリープバスター」を用いた。
【００６９】
被験者は６０歳代の健康な日本人男性であり、約４０分間の連続運転を異なる日に複数回実施した。被験者には、運転終了後、運転中の体調を自己申告させた。自己申告は、「良好」、「やや良好」、「普通」、「やや不良」、「不良」の５段階で行わせた。自己申告のうち、「やや良好」、「普通」、「やや不良」という曖昧なケースは除外し、被験者が明確に「良好」、「不良」と判定した５５例の走行データについて、本実施形態の手法による基礎的体調の推定結果とを比較した。
【００７０】
Ｉ．第１基礎的体調推定手段３１０による基礎的体調の推定
（１）周波数傾き時系列波形演算手段２１０から求められる指標が上記（ａ）の条件を満たすか否かの検討
（ａ）の条件は、２回連続で１つ前の振幅の９割未満に変化する波形の収束箇所が１回以上生じているか否かに設定し、被験者の主観と対比したカイ二乗検定を行った。２×２クロステーブルを次表に示すが、ｐ＝０．２４、正答率は６７％であった。
【００７１】
【表１】
0006666705.tif 000002
【００７２】
（２）分布率演算手段２２０から求められる指標が上記（ｂ）の条件を満たすか否かの検討
（ｂ）の条件は、９０秒間で、０．００１７Ｈｚの分布率が２０％以上減少し、その間に０．００５３Ｈｚの分布率が２０％以上増加する場合に設定し、被験者の主観と対比したカイ二乗検定を行った。２×２クロステーブルを次表に示すが、ｐ＝０．１１、正答率は６９％であった。
【００７３】
【表２】
0006666705.tif 000003
【００７４】
（３）疲労曲線演算手段２３０から求められる指標が上記（ｃ）の条件を満たすか否かの検討
（ｃ）の条件は、６分間の範囲で、ピーク検出法を用いた疲労曲線の傾きが、ゼロクロス法を用いた疲労曲線の傾きの１．５倍以上に変化する箇所が１箇所以上存在し、かつ、ピーク検出法を用いた疲労曲線が３０分経過時に０．０３より大きくなっている場合に設定し、被験者の主観と対比したカイ二乗検定を行った。２×２クロステーブルを次表に示すが、ｐ＝０．１９、正答率は６４％であった。
【００７５】
【表３】
0006666705.tif 000004
【００７６】
（４）恒常性維持機能レベル演算手段２４０から求められる指標が上記（ｄ）の条件を満たすか否かの検討
（ｄ）の条件は、６分間で上記の例でレベル４〜６という判定が１０回以上、又は、レベル７〜１１の判定が７回以上出現する場合に設定し、被験者の主観と対比したカイ二乗検定を行った。２×２クロステーブルを次表に示すが、ｐ＝０．０１８、正答率は７１％であった。
【００７７】
【表４】
0006666705.tif 000005
【００７８】
５５例の走行データのうち、上記（ａ）〜（ｄ）の条件を満たす数をまとめると次表のとおりである。
【００７９】
【表５】
0006666705.tif 000006
【００８０】
表５から明らかなように、（ａ）〜（ｄ）の条件のうち３つ以上満たす場合の７例は、いずれも被験者の主観においても「不良」と評価されている。従って、第１基礎的体調推定手段３１０においては、周波数傾き時系列波形演算手段２１０から求められる周波数傾き時系列波形、分布率演算手段２１０から求められる分布率の時系列波形、疲労曲線演算手段２３０から求められる疲労曲線（疲労度の時系列波形）、及び恒常性維持機能レベル演算手段２４０から求められる恒常性維持機能レベルのうち、３つ以上の指標がそれぞれ上記の（ａ）〜（ｄ）の条件を満たす場合に「不良」と推定するように設定することが好ましい。
【００８１】
（５）基礎的体調「良好」と判定される上記（ｅ）の条件を満たすか否かの検討
５５例の走行データの全てについて、体調マップ演算手段２５０により、測定開始から最初の判定結果を座標原点に合わせた際に２番目の判定結果の座標点が第４象限にプロットされるか否かを調べた。同じく、５５例の走行データの全てについて、感覚マップ演算手段２６０により、測定開始から最初の判定結果を座標原点に合わせた際に２番目の判定結果の座標点がＸ軸の正方向に所定目盛り以上（本実験例では、主観評価「良好」の３８例中多くの事例（２３例）の離隔距離が目盛り５以上であったため、当該目盛り５以上に設定）離隔してプロットされるか否かを調べた。体調マップ演算手段２５０及び感覚マップ演算手段２６０の演算結果が、両方とも、上記条件を満たす場合が、第１基礎的体調推定手段３１０において基礎的体調「良好」という推定結果となるが、被験者の主観と対比したカイ二乗検定の結果は、次表のとおりである。ｐ＝０．０９５、正答率は５６％であった。
【００８２】
【表６】
0006666705.tif 000007
【００８３】
上記表６の検討の際は、５５例全ての走行データを用いているが、図８のフローチャートに示したように、第１基礎的体調推定手段３１０は、まず、上記（ａ）〜（ｄ）の条件のうち３つ以上を満足する場合について基礎的体調を「不良」と推定し（Ｓ１１０）、「不良」と推定されない場合について、体調マップ及び感覚マップによる「良好」に相当するか否かの判断を行う（Ｓ１２０）。そこで、上記（ａ）〜（ｄ）の条件のうち３つ以上を満足する７例を除いた４８例についての比較検討したところ次の結果が得られた。この表より、７例を除いた場合には、ｐ＝０．０３２、正答率５６％であった。ｐ＜０．０５となっており、「不良」推定の条件を満たさない場合に、「良好」推定を行うことが好ましい。
【００８４】
【表７】
0006666705.tif 000008
【００８５】
ＩＩ．第２基礎的体調推定手段３２０による基礎的体調の推定
生体状態推定装置１００の恒常性維持機能レベル演算手段２４０によって求められる恒常性維持機能レベルのうち、レベル１〜６の合計出現数を１として、交感神経活動が優位な状態を示すレベル２及び副交感神経活動が優位な状態を示すレベル４の出現割合を求めた。これらを５５例の全てについて行い、ベイズ推定を用いて次の確率表を作成した。
【００８６】
【表８】
0006666705.tif 000009
【００８７】
表８の確率表を用いて第２基礎的体調推定手段３２０により推定される基礎的体調は、次表のとおりである。
【００８８】
【表９】
0006666705.tif 000010
【００８９】
上記表８及び表９の検討の際は、確率表作成のため、５５例全ての走行データを用いているが、図８のフローチャートに示したように、第２基礎的体調推定手段３２０は、第１基礎的体調推定手段３１０において基礎的体調が推定されない場合に体調推定を行う。そこで、第１基礎的体調推定手段３１０において「良好」又は「不良」と推定された２６例を除いた２９例について比較検討した結果が次表のとおりである。
【００９０】
【表１０】
0006666705.tif 000011
【００９１】
本実施形態の生体状態推定装置１００は、図８のフローチャートに示したように、推定対象のデータについて、第１基礎的体調推定手段３１０において、まず、「不良」の状態を推定し（Ｓ１１０）、「不良」と推定されない場合に、「良好」の状態を推定し（Ｓ１２０）、そのいずれとも推定されない場合に、第２基礎的体調推定手段において、表８に示したベイズ推定による確率表に基づき、「良好」、「不良」、それらの「中間状態」（上記表では「中間」と表示）のいずれかであるかを推定するが（Ｓ１３０）、本実験例では、被験者が明確に「良好」又は「不良」と評価した５５例を用いている。
【００９２】
従って、これら５５例中、生体状態推定装置１００において、最終的に「中間」と推定されたデータは、被験者の主観評価とは合致しない。そこで、この「中間」のデータ（１１例）を除くと共に、第１基礎的体調推定手段３１０において「不良」と推定されたデータ（不良：７例）、その「不良」のデータを除いた中で、「良好」と推定されたデータ（良好：１９例）、並びに、第１基礎的体調推定手段３１０において「不良」、「良好」と推定されたデータを除き、第２基礎的体調推定手段３２０において「不良」、「良好」と推定されたデータ（不良：７例、良好：１１例）について、すなわち、重複のないデータについて最終的に基礎的体調出力手段４００によって出力される推定結果を被験者の主観評価と比較してまとめると次表のとおりになった。
【００９３】
【表１１】
0006666705.tif 000012
【００９４】
表１１によれば、ｐ＝１．５８×１０−７で、正答率９１％であり、本実施形態の生体状態推定装置１０による基礎的体調の推定結果が、高い確率で被験者の主観評価と一致し、本実施形態の基礎的体調の推定結果の信頼度が高いことがわかる。
【００９５】
なお、５５例中、第１基礎的体調推定手段３１０において「良好」又は「不良」と推定されたデータを除いて、第２基礎的体調推定手段３２０において「中間」と推定された１１例のデータは、被験者の主観評価とは合致しないが、これは、被験者の基本的な体力や思考傾向によって左右されるものと考えられる。すなわち、本実施権例の被験者の場合、ポジティブ思考の傾向が強く、体力もあるため、１１例中の多くは、主観的には「良好」と自己申告したものと考えられる。
【００９６】
（実験例２）
次に、運輸事業者のトラックの運転席に、生体信号測定装置１として、（株）デルタツーリング製、スリープバスター（登録商標）を装着し、職業運転手の生体信号を採取して、実験例１で示した各推定基準が設定された生体状態推定装置１００により、図８に示したフローチャートに従った分析を実施し、被験者の３０分毎の体調を、「良」、「中間」、「不良」の３段階で推定した。被験者である職業運転手は、２０歳代から５０歳代までの男性７名で、全被験者の総運行回数７１回を対象に分析を行った。また、各被験者には、運行当日の自身の主観的な体調を、運行終了後、「絶不調」（レベル：１）から「絶好調」（レベル：５）までの５段階で申告させた。
【００９７】
なお、「運行回数」は、各被験者の勤務日１日（夜勤など、カレンダー上の日付を跨ぐ場合もある）における運行開始時点から当該勤務日の運行終了時点までを１回（一運行）とカウントし、その中には、数時間連続の運転やその間の休憩時間等も含んでいる。一方、生体状態推定装置１００による推定結果は約３０分に１回出力される。従って、一運行において、１時間以上の連続運転が含まれる場合には、生体状態推定装置１００により複数回の推定結果が得られる。その一方、各被験者の運行当日の主観評価は、一運行終了時点における１回の申告である。
【００９８】
図９（ａ），（ｂ）は被験者Ａ及びＢのある月の結果を示し、生体状態推定装置１００による一運行中の体調推定の全判定数のうち、「良」と推定された割合（各グラフの右側縦軸）と、被験者の５段階の当日運行終了後の主観的な体調評価（各グラフの左側縦軸）を運行日別にプロットして時系列に示したものである。
【００９９】
図９（ａ），（ｂ）から、いずれの被験者のデータも、生体状態推定装置１００により「良」と判定された割合の変化の増減傾向と、主観評価の変化の傾向とが近似していることがわかる。表１２は、被験者Ａの体調推定の結果及び主観評価を、前日（前回の運行）との比較で示したもので、前日との比較で、体調推定「良」の割合が増加している場合及び主観評価がよくなっている場合を「＋」で表示し、体調推定「良」の割合が減少している場合及び主観評価が悪くなっている場合を「−」で表示し、いずれも前日と同じ場合には「±」で表示した。また、増減傾向が一致する場合を丸印で示した。
【０１００】
【表１２】
0006666705.tif 000013
【０１０１】
表１２によれば、主観評価が前日と同じ場合を除いて、主観と体調推定が近似した傾向になることが示唆される。
【０１０２】
表１３は、表１２と同様の前日（前回の運行）との比較を行った全被験者の結果をまとめて示したものであり、７名の被験者全員の増減傾向の一致率は０．７１で、フィッシャーの正確確率検定の結果、ｐ＝０．０３２となり、生体状態推定装置１００による体調推定の「良」の割合の増減と主観評価には有意な関連が認められた。また、体調推定の「良」の割合が「増加」したときに主観評価が「悪い」となったデータ数は有意に少なく、体調推定「良」の割合が「減少」したときに主観評価が「悪い」となったデータ数は有意に多い結果となった。その一方、主観評価が、前日の運行と同一であるとした回答数が７１運行中３２件、「普通（レベル：３）」との回答が７１運行中３８件あり、日々の体調の変化を厳密に自覚することが比較的難しいことも示唆された。
【０１０３】
【表１３】
0006666705.tif 000014
【０１０４】
ここで、図９（ｂ）の被験者Ｂのデータにおいて、主観評価で１（絶不調）と回答した１４日と１５日の運行に着目する。図１０は、被験者Ｂの１４日と１５日における生体状態推定装置１００による体調推定結果の時系列データと、８日から１６日までの運行開始時刻を示している。１４日の運行は、運行開始時が１５時４８分で運行終了時が１５日午前３時過ぎであり、１５日の運行は、運行開始時が１８時４６分で運行終了時が１６日午前３時過ぎであった。
【０１０５】
１４日は、体調推定結果が「良」と「不良」を行き来し、終盤には「不良」の推定が連続している。１４日は、被験者Ｂの運行後のアンケートで軽い眠気が運行終盤に生じた旨を回答していることから、１１日までの早朝・昼勤務から、１４日以降、午後・夜勤務に変わったこと、並びに、１２日、１３日が休暇であったことにより、体の仕事モードや夜勤モードへの順応、切り替わりが十分ではなく、眠気や体調の不良を感じたものと推定される。１５日は、前半においてほぼ「不良」との推定であるが、後半には「良」の推定も出ていることから、１５日の主観評価が１４日のレベル：１（絶不調）からレベル：２（不調）へと変化したことに対応している。従って、本実施形態の生体状態推定装置１００による体調推定結果は、このように、被験者の主観評価に近い状態を推定していると言えると共に、約３０分に１回出力することで、被験者が自覚しにくい体調の変化を反映できることも示唆される。
【０１０６】
図１１は、全被験者の主観評価と生体状態推定装置１００による体調推定の「良」の割合との相関を正規化した分布図であり、被験者毎に分布の線形近似を求めて細線で示している。図１１において４５度斜めに引いた破線より右側に位置するデータは、主観評価に比べて体調推定結果「良」の割合が高いことを示し、破線より左側に位置するデータは、主観評価に比べて体調推定結果「良」の割合が低いことを示している。また、各被験者の線形近似の直線は、傾きが４５度に近いほど主観評価と体調推定結果との間に正の相関が高く、傾き０度に近づくほど相関性が低いことを示し、傾きがマイナスの場合には負の相関が高いことを示す。主観評価が５段階に対し、生体状態推定装置１００による体調推定が「良」、「中間」、「不良」の３つであることから、主観評価に対して体調推定結果が左右方向へ広がる傾向があるものの、概ね正の傾きを示しており、図１１からも、生体状態推定装置１００による体調推定と主観評価との間に一定程度の相関があることが示唆される。
【０１０７】
なお、上記した実験例では、車載のコンピュータである生体状態推定装置１００の設定されたコンピュータプログラム及びその記憶部に記憶されたデータを用いて解析を行っているが、これらのプログラムや各データは、運行管理者のコンピュータに設定し、上記と同様な解析を行うことができることはもちろんである。この場合、解析対象のデータを、車載の生体状態推定装置１００から通信回線（無線等）を介して運行管理者のコンピュータに送信し、運行管理者のコンピュータにおいてもリアルタイムに解析できるようにしてもよいし、また、運転業務終了後に、車載の生体状態推定装置１００に記憶されているデータを取り出し、運行管理者のコンピュータにおいて事後解析するようにしてもよい。運行管理者は、このようなデータを収集することで、運転者毎に、運転時の状況を把握できると共に、より適切な運転をするためのアドバイスの付与などにも活用できる。
【符号の説明】
【０１０８】
１  生体信号測定装置
１１  コアパッド
１２  スペーサパッド
１３  センサ
１００  生体状態推定装置
２００  生体調節機能要素判定手段
２１０  周波数傾き時系列波形演算手段
２２０  分布率演算手段
２３０  疲労曲線演算手段
２４０  恒常性維持機能レベル演算手段
２５０  体調マップ演算手段
２６０  感覚マップ演算手段
３００  基礎的体調推定手段
３１０  第１基礎的体調推定手段
３２０  第２基礎的体調推定手段
【図１】
0006666705.tif 000015
【図２】
0006666705.tif 000016
【図３】
0006666705.tif 000017
【図４】
0006666705.tif 000018
【図５】
0006666705.tif 000019
【図６】
0006666705.tif 000020
【図７】
0006666705.tif 000021
【図８】
0006666705.tif 000022
【図９】
0006666705.tif 000023
【図１０】
0006666705.tif 000024
【図１１】
0006666705.tif 000025