(57)【要約】【課題】  対象者の気分の変動量を推定することができる、気分推定システムを提供すること。【解決手段】  気分推定システム１が備える気分変動量推定部６７は、Ｘ軸が安静状態であるときにおける生体情報でありＹ軸が非安静状態であるときにおける生体情報であるＸ−Ｙ平面上において、Ｘ成分が安静状態であるときにおける対象者の生体情報の平均値でありＹ成分が対象者の生体情報の現在値（生体情報現在値）である座標点Ａと、安静状態であるときにおける生体情報と非安静状態であるときにおける生体情報との相関関係を表す相関線との距離である乖離距離Ｌに基づいて、現在時点における対象者の気分の変動量を推定する。【選択図】  図８Ｂ
2019017946.tif 000002
【特許請求の範囲】
【請求項１】
心拍情報を含む対象者の現在時点における生体情報であって、気分の変動に影響する生体情報の現在値である生体情報現在値を取得する生体情報現在値取得部と、
前記生体情報現在値取得部が取得した対象者の前記生体情報現在値に基づいて、対象者の気分状態が安静状態であるか非安静状態であるかを判定する気分状態判定部と、
前記気分状態判定部により対象者の気分状態が非安静状態であると判定されたときに、現在時点における対象者の気分の変動量を推定する気分変動量推定部と、
を備え、
前記気分変動量推定部は、Ｘ軸が安静状態であるときにおける生体情報でありＹ軸が非安静状態であるときにおける生体情報であるＸ−Ｙ平面上において、Ｘ成分が安静状態であるときにおける対象者の生体情報の平均値でありＹ成分が前記生体情報現在値である座標点と、安静状態であるときにおける生体情報と非安静状態であるときにおける生体情報との相関関係を表す相関線との距離である乖離距離に基づいて、現在時点における対象者の気分の変動量を推定する、
を備える、気分推定システム。
【請求項２】
請求項１に記載の気分推定システムにおいて、
前記気分変動量推定部は、現在時点における対象者の気分のうち変動している気分の変動量を推定する、気分推定システム。
【請求項３】
請求項１又は２に記載の気分推定システムにおいて、
前記気分変動量推定部は、前記乖離距離が大きいほど、現在時点における対象者の気分の変動量が大きいと推定する。気分推定システム。
【請求項４】
請求項１乃至３のいずれか１項に記載の気分推定システムにおいて、
対象者の現在の行動状態を推定する行動状態推定部を備え、
前記気分変動量推定部は、前記Ｘ−Ｙ平面上に表される前記相関線からの前記座標点の乖離方向及び、前記行動状態推定部が推定した対象者の現在の行動状態に基づいて、現在時点における対象者の気分のうち変動している気分を特定する、気分推定システム。
【請求項５】
請求項４に記載の気分推定システムにおいて、
対象者のストレスの強弱及び心理的負担の大小を推定するストレス推定部をさらに備える、気分推定システム。
【請求項６】
請求項５に記載の気分推定システムにおいて、
前記ストレス推定部は、安静状態であるときにおける対象者のＲＲＩデータに基づいて求められるローレンツプロットである安静時ローレンツプロットのＬＰ面積及び非安静状態であるときにおける対象者のＲＲＩデータに基づいて求められるローレンツプロットである非安静時ローレンツプロットのＬＰ面積に基づいて、対象者のストレスの強弱を推定し、前記安静時ローレンツプロットのＭ値及び前記非安静時ローレンツプロットのＭ値に基づいて、対象者の心理的負担の大小を推定する、気分推定システム。
【請求項７】
請求項５又は６に記載の気分推定システムにおいて、
前記気分変動量推定部が推定した現在時点において変動している気分及びその変動量、前記ストレス推定部が推定した対象者のストレスの強弱及び心理的負担、前記体動状態推定部が推定した対象者の現在時点における行動状態に基づいて、対象者の行動状態、変動している気分及びその変動量、ストレスの強弱及び心理的負担の関連付けを行う統合判断部を備える、心理的刺激推定システム。
【請求項８】
請求項１乃至７のいずれか１項に記載の気分推定システムにおいて、
前記気分変動量推定部が現在時点における対象者の気分の変動量を推定する場合に用いる前記生体情報が、ＲＲＩ分散である、気分推定システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、気分推定システムに関する。
【背景技術】
【０００２】
対象者の気分、或いは、心理状態、を推定するシステムが開発されている。例えば、特許文献１は、対象者の顔に持続的に表れる主表情の種類の変化、及び、対象者の顔に瞬間的に表れる微表情の有無を検出し、検出された主表情の種類の変化及び微表情の有無に基づいて、対象者の心理状態を推定するシステムを開示する。また、特許文献２は、対象者の生体データを収集し、収集した生体データに対して設定された推定区間内における心理状態の変化を推定する方法を開示する。また、特許文献３は、時針及び分針を有する腕時計式の装置の時針及び分針により表される時間軸の周に沿って、対象者の生理学的パラメータから得られた指標を表示する装置を開示する。また、特許文献４は、生体情報に基づいて推定された心理状態と、体動情報又は位置情報に基づいて推定された活動状態とに基づいて、対象者の嗜好を判定し、嗜好の判定に基づいて、視聴が推奨される推奨コンテンツを決定するように構成された画像表示システムを開示する
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１６−１４９０６３号公報
【特許文献２】特開２０１４−１４４０５２号公報
【特許文献３】特許第５０８１２２３号明細書
【特許文献４】特開２０１６−９２５１０号公報
【発明の概要】
【０００４】
（発明が解決しようとする課題）
上記特許文献記載の技術によれば、対象者の気分或いは心理状態を推定することができるものの、対象者の気分の変動量まで推定することができない。そこで、本発明は、対象者の気分の変動量を推定することができる、気分推定システムを提供することを目的とする。
【０００５】
（課題を解決するための手段）
本発明に係る気分推定システムは、心拍情報を含む対象者の現在時点における生体情報であって、気分の変動に影響する生体情報の現在値である生体情報現在値を取得する生体情報現在値取得部（６１）と、生体情報現在値取得部が取得した対象者の生体情報現在値に基づいて、対象者の気分状態が安静状態であるか非安静状態であるかを判定する気分状態判定部と（６２）、気分状態判定部により対象者の気分状態が非安静状態であると判定されたときに、現在時点における対象者の気分の変動量を推定する気分変動量推定部（６７）と、を備える。そして、気分変動量推定部は、Ｘ軸が安静状態であるときにおける生体情報でありＹ軸が非安静状態であるときにおける生体情報であるＸ−Ｙ平面上において、Ｘ成分が安静状態であるときにおける対象者の生体情報の平均値でありＹ成分が生体情報現在値である座標点と、安静状態であるときにおける生体情報と非安静状態であるときにおける生体情報との相関関係を表す相関線との距離である乖離距離に基づいて、現在時点における対象者の気分の変動量を推定する。
【０００６】
本発明に係る気分推定システムによれば、安静時の生体情報と非安静時の生体情報との相関線に対して、対象者の現在の生体情報（生体情報現在値）がどの程度乖離しているかに基づいて、現在時点における対象者の気分の変動量が定量的に推定される。このように、本発明によれば、対象者の気分の変動量を推定することができる気分推定システムを提供することができる。
【０００７】
上記発明において、「気分」とは、「感情」或いは「情動」と言い換えることもでき、人が物事に対して抱く気持である。「気分」は、複数の種類に分類され得る。例えば、「好き嫌い」、「快適感」、「焦燥感」、「集中」、「眠気」、「喜怒哀楽」等は、気分の一種であり、それぞれが、気分を構成する。
【０００８】
また、上記発明において、「安静状態」とは、気分を表す状態が発現していない状態を言い、「非安静状態」とは、気分を表す状態が発現している状態を言う。また、上記発明における「相関関係」は、安静状態であるときにおける生体情報と非安静状態であるときにおける生体情報との関係であって複数の対象者に対して適用され得る一般的な相関関係である。この相関関係は、複数の対象者（被験者）から得られる安静状態であるときにおける生体情報と非安静状態であるときにおける生体情報に基づいて構築される。
【０００９】
また、上記発明において、「気分の変動量」とは、ある種類の気分に対して、その気分が発現したときにおけるその気分についての安静時からの気持の強さ（刺激の強さ）の変化量を意味する。例えば、「好き嫌い」という気分に対して「嫌悪感の度合い」、「焦燥感」という気分に対して「焦りの度合い」が、気分の変動量（刺激の強さ）である。
【００１０】
また、上記発明において、「生体情報現在値」とは、心拍数、心拍間隔（ＲＲＩ）、皮膚温度、呼吸数、発汗量（皮膚電位）等の、生体現象によって体内から発せられる信号（生体信号）の現在値及び、生体信号の統計値（平均、分散）の現在値を表す。なお「生体信号の統計値の現在値」は、現在時点に至るまでの微小時間内に取得された生体信号データを統計処理することにより算出され得る。
【００１１】
気分変動量推定部は、現在時点における対象者の気分のうち変動している気分の変動量を推定するのがよい。これによれば、対象者の現在時点において変動している気分の変動量を推定することができる。ここで、「変動している気分」とは、非安静状態である気分を言う
【００１２】
気分変動量推定部は、上記乖離距離が大きいほど、現在時点における対象者の気分の変動量が大きいと推定するのがよい。上記相関線からの対象者の生体情報現在値を表す座標点の乖離距離が大きいほど、その生体情報に影響する気分がより大きく変動していると考えられる。従って、乖離距離が大きいほど気分の変動量が大きいと推定することにより、気分の変動量を精度良く推定することができる。
【００１３】
また、本発明に係る気分推定システムは、対象者の現在の行動状態を推定する行動状態推定部（６６）を備えるとよい。そして、気分変動量推定部は、Ｘ−Ｙ平面上に表される相関式からの対象者の生体情報現在値を表す座標点の乖離方向及び、行動状態推定部が推定した対象者の現在の行動状態に基づいて、現在時点における対象者の気分のうち変動している気分を特定するのがよい。これによれば、変動している気分を精度良く特定することができる。
【００１４】
また、本発明に係る気分推定システムは、対象者のストレスの強弱及び心理的負担の大小を推定するストレス推定部（６８）をさらに備えるとよい。これによれば、対象者の気分の変動量のみならず、気分の変動に伴って対象者が受けるストレスの強弱及び心理的負担の大きさをも推定することができる。
【００１５】
この場合、ストレス推定部は、安静状態であるときにおける対象者のＲＲＩデータに基づいて求められるローレンツプロットである安静時ローレンツプロットのＬＰ面積及び非安静状態であるときにおける対象者のＲＲＩデータに基づいて求められるローレンツプロットである非安静時ローレンツプロットのＬＰ面積に基づいて、対象者のストレスの強弱を推定し、安静時ローレンツプロットのＭ値及び非安静時ローレンツプロットのＭ値に基づいて、対象者の心理的負担の大小を推定するのがよい。
【００１６】
これによれば、安静時ローレンツプロットのＬＰ面積及びＭ値と、非安静時ローレンツプロットのＬＰ面積及びＭ値を比較することにより、精度良く、ストレスの強弱及び心理的負担の大小を推定することができる。
【００１７】
また、本発明に係る気分推定システムは、気分変動量推定部が推定した現在時点において変動している気分及びその変動量、ストレス推定部が推定した対象者のストレスの強弱及び心理的負担、体動状態推定部が推定した対象者の現在時点における行動状態に基づいて、対象者の行動状態、変動している気分及びその変動量、ストレスの強弱及び心理的負担の関連付けを行う統合判断部（６９）を備えるのがよい。これによれば、対象者の行動状態、そのときに変動している気分及びその変動量、ストレスの強弱及び心理的負担の大きさが、どのように関連しているかを推定することができる。
【００１８】
また、気分変動量推定部が現在時点における対象者の気分の変動量を推定する場合に用いる生体情報が、ＲＲＩ分散であるとよい。ここで、ＲＲＩ分散とは、心拍間隔（Ｒ−Ｒ  ｉｎｔｅｒｖａｌ）の分散である。安静状態であるときにおけるＲＲＩ分散と非安静状態であるときにおけるＲＲＩ分散の相関式に対する対象者のＲＲＩ分散の現在値を表す座標点の乖離距離は、対象者の気分の変動量を精度良く表す。従って、ＲＲＩ分散についての乖離距離に基づいて、気分の変動量をより精度良く推定することができる。なお、ＲＲＩ分散の現在値は、現在時点に至るまでの微小時間の間に取得される対象者の複数のＲＲＩデータに基づいて算出することができる。
【図面の簡単な説明】
【００１９】
【図１】本実施形態に係る気分推定システムの概略構成を示す図である。
【図２】サーバの機能構成を示すブロック図である。
【図３】気分状態判定部実行する気分状態判定処理ルーチンの流れを示すフローチャートである。
【図４】焦燥感に関する発生確率ｐの推移を時系列的に示すグラフの一例である。
【図５】安静時における心拍数の平均値と非安静時における心拍数の平均値との相関関係を表すグラフの一例である。
【図６】安静時における皮膚温度の平均値と非安静時における皮膚温度の平均値との相関関係を表すグラフの一例である。
【図７】行動状態推定部が実行する行動状態推定処理ルーチンの流れを示すフローチャートである。
【図８Ａ】気分変動量推定部が実行する気分変動量推定処理ルーチンの流れを示すフローチャートである。
【図８Ｂ】図８Ａに続くフローチャートである。
【図９】Ｘ軸が安静時における生体情報でありＹ軸が非安静時における生体情報であるＸ−Ｙ座標平面上に、相関式Ｆ（Ｘ）と、座標点Ａとが示されたグラフの一例である。
【図１０Ａ】ストレス推定部が実行するストレス推定処理ルーチンの流れを示すフローチャートである。
【図１０Ｂ】図１０Ａに続くフローチャートである。
【図１１】ローレンツプロットの一例を示す図である。
【図１２】統合処理部が実行する統合処理ルーチンの流れを示すフローチャートである。
【図１３】環境制御部が実行する環境制御処理ルーチンの流れを示すフローチャートである。
【発明を実施するための形態】
【００２０】
以下、本発明の実施形態について図面を参照して説明する。図１は、本実施形態に係る気分推定システムの概略構成を示す図である。図１に示すように、本実施形態に係る気分推定システム１は、生体信号検出センサ２と、体動検出センサ３と、情報端末４と、中継器（モデム）５と、サーバ６とを備える。
【００２１】
生体信号検出センサ２は、少なくとも心拍数或いはＲＲＩ等の心拍信号を含み、対象者の様々な気分に影響すると考えられる生体信号を検出する。生体信号検出センサ２は、対象者の複数の生体信号を個々に検出する複数のセンサにより構成されるセンサ群であってもよい。センサ群を構成する複数のセンサとして、対象者の心拍信号（心拍間隔ＲＲＩ、心拍数等）を検出する心拍信号検出センサ、対象者の脈波信号を検出する脈波信号検出センサ、対象者の皮膚電位（ここで、皮膚電位の大きさは、発汗量に影響する）を検出する皮膚電位検出センサ、呼吸信号（呼吸数、呼吸間隔等）を検出する呼吸信号検出センサ、体温（皮膚温度）を検出する体温検出センサ、血圧検出センサ等を例示することができる。また、体動検出センサ３は、対象者の動きを加速度信号として検出するとともに、検出した加速度信号を出力する。生体信号検出センサ２及び体動検出センサ３は、ウェアラブルセンサとして対象者が身に着けているとよい。生体信号検出センサ２及び体動検出センサ３は、有線通信手段或いは無線通信手段により、中継器５と通信可能に構成されており、検出した信号を、中継器５に送信する。
【００２２】
情報端末４は、対象者が携帯しており、無線通信手段によりサーバ６と送受信可能に構成される。情報端末４は、サーバ６から送信された情報を表示する機能を有する。情報端末４として、例えばスマートフォンを例示することができる。
【００２３】
中継器５は、生体信号検出センサ２が出力した生体信号及び体動検出センサ３が出力した加速度信号を逐次的に受信する。また、中継器５は、サーバ６と無線通信により送受信可能に構成される。中継器５は、生体信号検出センサ２から受信した対象者の生体信号及び体動検出センサ３から受信した対象者の加速度信号を、無線通信によってサーバ６に送信する。また、中継器５は、複数の環境制御機器７と通信可能に構成される。環境制御機器７は、対象者の周辺環境（温度、明るさ等）を制御することができる機器である。環境制御機器７として、空調装置、照明、音響機器等を例示することができる。
【００２４】
サーバ６は、データの送受信、データの記憶・保管、データの解析・判定・統合処理、等を行うことができるハードウェアである。サーバ６は、例えばクラウドコンピューティング上に設置することができる。
【００２５】
また、本実施形態に係る気分推定システム１は、複数の対象者についての気分を推定することができるように構成される。従って、気分推定システム１は、図１に示すように、複数の対象者（図１においては対象者Ａ、対象者Ｂ、対象者Ｃ）のそれぞれに対して割り当てられた、生体信号検出センサ２、体動検出センサ３、情報端末４、中継器５を有することになる。
【００２６】
図２は、サーバ６の機能構成を示すブロック図である。図２に示すように、サーバ６は、現在値取得部６１と、気分状態判定部６２と、情報記憶部６３と、統計解析部６４と、相関関係導出部６５と、行動状態推定部６６と、気分変動量推定部６７と、ストレス推定部６８と、統合処理部６９と、環境制御部７０と、を有する。
【００２７】
現在値取得部６１は、生体信号検出センサ２及び体動検出センサ３が検出した複数の対象者の直近の生体信号及び加速度信号を、中継器５を介して入力する。そして、入力された信号に基づいて、対象者の現在時点における生体情報（生体情報現在値）及び体動情報（体動情報現在値）を取得する。現在値取得部６１が、本発明の生体情報現在値取得部に相当する。
【００２８】
ここで、生体信号検出センサ２が検出する信号には、心拍信号が含まれる。また、生体信号検出センサ２が検出する信号は、気分の変動に影響する生体信号である。従って、現在値取得部６１は、心拍情報を含む対象者の現在時点における生体情報であって、気分の変動に影響する生体情報現在値を取得することになる。また、現在値取得部６１が取得する生体情報現在値には、生体信号の現在値に加え、生体信号の統計値（平均及び分散）の現在値が含まれる。例えば、対象者のＲＲＩの平均値の現在値や分散（ＲＲＩ分散）の現在値も、生体情報現在値に含まれる。この場合において、生体信号の平均値の現在値は、現在時点に至るまでの微小時間（例えば１分間）の間に入力した生体信号データを平均することにより取得することができ、生体信号の分散の現在値は、現在時点に至るまでの微小時間の間に入力した生体信号データの平均値及びその平均値と各生体信号データとの偏差に基づいて取得することができる
【００２９】
気分状態判定部６２は、対象者の気分が安静状態であるか非安静状態であるかを判定する。ここで、「安静状態」とは、複数種類の気分の全てについて、安静な状態（その気分が発現されていない状態）を言う。例えば、複数種類の気分が、「焦燥感」、「集中」「好き嫌い」であるとすると、焦燥感及び集中度が低く、且つ、好きという感情も嫌いという感情もないような状態が、「安静状態」である。一方、「非安静状態」とは、上記した複数種類の気分の少なくとも一つの気分が安静でない状態（その気分が発現されている状態）を言う。
【００３０】
気分状態判定部６２は、気分状態判定処理を実行することにより、対象者の現在の気分状態が安静状態であるか非安静状態であるかを判定する。図３は、気分状態判定部６２が、気分状態を判定するために実行する気分状態判定処理ルーチンの流れを示すフローチャートである。図３に示すルーチンが起動すると、気分状態判定部６２は、まず、図３のＳ１０１にて、現在値取得部６１が取得した対象者の生体情報現在値を読み出す。次いで、気分状態判定部６２は、対象者の個々の気分が安静状態であるという事象の発生確率ｐを、全ての種類の気分について演算する。
【００３１】
ある種の気分が安静状態であるか非安静状態であるかは、少なくとも心拍情報を含む特定の生体情報の大きさに現れる。従って、心拍情報を含む特定の生体情報から、その気分が安静状態であるか非安静状態であるかを推定することができる。本実施形態においては、ロジスティック回帰式を用いて、気分が安静状態であるか非安静状態であるかが判定される。この場合、気分状態判定部６２は、Ｓ１０２にて、気分が安静状態であるという事象の発生確率ｐの対数オッズ（ｌｎ（ｐ／（１−ｐ））を従属変数、少なくとも心拍情報を含む生体情報ｘｋを独立変数とした以下の（１）式に示すようなロジスティック回帰式
【数１】
2019017946.tif 000003
を用いたロジスティック回帰分析により、発生確率ｐを演算する。ここで、ｋは次元数と呼ばれるデータの種類数を示し、ｉはデータ数を示す。また、生体情報ｘｋについての最適な回帰係数αｋは、以下の（２）式に示す目的関数Ｅを最大化する最尤法により求められる。
【数２】
2019017946.tif 000004
ここで、（２）式において、ｔは確率の正解値である。
【００３２】
また、（１）式を変形することにより、以下の（３）式が導出される。
【数３】
2019017946.tif 000005
【００３３】
最尤法により求められる回帰係数と心拍情報を含む生体情報とを（３）式に代入することにより、発生確率ｐを演算することができる。一般的には、得られた発生確率ｐに１００を乗じて、百分率により発生確率ｐが表現される。
【００３４】
また、上述したように、気分には、複数の種別が存在する。従って、上記した発生確率ｐは、それぞれの気分に対応して求めることができる。例えば、「焦燥感」という気分に対しては、以下の（４）式に示すロジスティック回帰式が与えられる。
【数４】
2019017946.tif 000006
【００３５】
（４）式からわかるように、「焦燥感」は、生体情報としての皮膚電位、ＲＲＩ（心拍間隔）、ＲＲＩ分散（心拍間隔の分散）に影響する。よって、気分状態判定部６２は、現在値取得部６１から読み出した対象者の皮膚電位の現在値、ＲＲＩの現在値、及び、ＲＲＩ分散の現在値に基づいて、「焦燥感」についての気分が安静状態であるという発生確率ｐを求めることができる。
【００３６】
図４は、焦燥感に関する発生確率ｐの推移を時系列的に示すグラフの一例である。図４の横軸が、焦燥感に関する発生確率ｐの演算開始からの経過時間（ｓｅｃ．）であり、縦軸が発生確率ｐ（％）である。図４に示すように、経過時間が２４１秒以前の区間においては、発生確率ｐが原則的に５０％以上であり、経過時間が２４１秒以降の区間においては、発生確率ｐが原則的に５０％未満である。つまり、２４１秒の時点を境に、発生確率ｐが変化していることがわかる。
【００３７】
気分状態判定部６２は、Ｓ１０２にて上記した手法に基づいて対象者の現在の各気分が安静状態であるという事象の発生確率ｐを演算した後に、Ｓ１０３にて、演算した発生確率ｐが、全ての気分について５０％以上であるか否かを判断する。全ての気分についての発生確率ｐが５０％以上である場合（Ｓ１０３：Ｙｅｓ）、気分状態判定部６２は、Ｓ１０４に処理を進めて、対象者の現在の気分状態が安静状態であると判定し、次いで、Ｓ１０６に処理を進めて、Ｓ１０１にて読み出した生体情報現在値を安静時生体情報として情報記憶部６３に出力する。その後、気分状態判定部６２は、このルーチンを終了する。一方、全ての気分についての発生確率ｐが５０％以上ではない場合（Ｓ１０３：Ｎｏ）、気分状態判定部６２は、Ｓ１０５に処理を進める。なお、Ｓ１０３の判定結果がＮｏとなる場合は、全ての気分についての発生確率ｐが５０％未満である場合のみならず、ある種類の気分が５０％以上でありそれ以外の気分についての発生確率ｐが５０％未満である場合を含む。
【００３８】
Ｓ１０５では、気分状態判定部６２は、対象者の現在の気分状態が非安静状態であると判定し、次いで、Ｓ１０７に処理を進めて、Ｓ１０１にて読み出した生体情報現在値を非安静時生体情報として情報記憶部６３に出力する。その後、気分状態判定部６２は、このルーチンを終了する。
【００３９】
気分状態判定部６２が上記した気分状態判定処理を微小時間間隔で繰り返し実行することにより、対象者の現在の気分の状態が安静状態であるか非安静状態であるかが逐次的に判定される。
【００４０】
情報記憶部６３は、気分状態判定部６２から安静時生体情報を入力した場合、入力した生体情報を、安静時における生体情報として、対象者ごとに記憶する。また、情報記憶部６３は、気分状態判定部６２から非安静時生体情報を入力した場合、入力した生体情報を、非安静時における生体情報として、対象者ごとに記憶する。情報記憶部６３により、対象者ごとに、安静状態であるときにおける生体情報及び非安静状態であるときにおける生体情報が、累積的に記憶される。
【００４１】
統計解析部６４は、情報記憶部６３に記憶されているそれぞれの対象者の安静時生体情報の統計値である平均値及び分散、並びに、非安静時生体情報の統計値である平均値及び分散を、対象者ごとに、且つ、気分の種類ごとに、演算する。なお、統計解析部６４が演算する各生体情報の平均値及び分散は、現在値取得部６１が取得する生体情報現在値に含まれる生体信号の平均値の現在値或いは分散の現在値とは異なる。現在値取得部６１が取得した生体信号の平均値の現在値のうち安静時に取得した生体信号の平均値の現在値（生体情報現在値）の平均値が、統計解析部６４にて演算される安静時生体情報の平均値と言え、現在値取得部６１が取得した生体信号の平均値の現在値のうち非安静時に取得した生体信号の平均値の現在値（生体情報現在値）の平均値が、統計解析部６４にて演算される非安静時生体情報の平均値と言える。また、現在値取得部６１が取得した生体信号の分散の現在値のうち安静時に取得した生体信号の分散の現在値（生体情報現在値）の平均値が、統計解析部６４にて演算される安静時生体情報の分散と言え、現在値取得部６１が取得した生体信号の分散の現在値のうち非安静時に取得した生体信号の分散の現在値（生体情報現在値）の平均値が、統計解析部６４にて演算される非安静時生体情報の分散と言える。統計解析部６４が演算した安静時生体情報の平均値及び分散、並びに非安静時生体情報の平均値及び分散は、相関関係導出部６５に出力される。
【００４２】
相関関係導出部６５は、気分状態が安静状態であるとき（安静時）における生体情報と、気分状態が非安静状態であるとき（非安静時）における生体情報との相関関係を表す相関式を導出する。本実施形態では、相関関係導出部６５は、生体情報の平均値及び分散に関しての相関式を導出する。この場合、相関関係導出部６５は、複数の対象者についての安静時生体情報の平均値及び非安静時生体情報の平均値を統計解析部６４から入力するとともに、入力した複数の対象者についてのこれらの平均値に基づいて、安静時における生体情報の平均値と非安静時における生体情報の平均値との相関関係を表す相関式を導出する。また、相関関係導出部６５は、複数の対象者についての安静時生体情報の分散及び非安静時生体情報の分散を統計解析部６４から入力するとともに、入力した複数の対象者についてのこれらの分散に基づいて、安静時における生体情報の分散と非安静時における生体情報の分散との相関関係を表す相関式を導出する。相関関係は、例えば、最小二乗法を用いて導出することができる。
【００４３】
なお、相関式の導出に用いるデータ数が少ない場合、十分に精度の高い相関式を導出することができない。この点に関し、本実施形態においては、相関関係導出部６５には、予め複数の対象者に対して実験により得られた、安静時における生体情報の統計値（平均値、分散）と非安静時における生体情報の統計値（平均値、分散）との相関関係を表す相関式が記憶されている。そして、相関関係導出部６５は、統計解析部６４から入力した安静時生体情報の平均値及び分散並びに非安静時生体情報の平均値及び分散を用いて、予め記憶されている相関式を補正することにより、相関式を導出する。これにより、相関式の精度が向上する。
【００４４】
図５は、相関関係導出部６５が導出した、安静時における心拍数の平均値と非安静時における心拍数の平均値との相関関係を表すグラフの一例であり、図６は、相関関係導出部６５が導出した、安静時における皮膚温度の平均値と非安静時における皮膚温度の平均値との相関関係を表すグラフの一例である。図５及び図６において、横軸が安静時における生体情報（心拍数、皮膚温度）の平均値であり、縦軸が非安静時（タスク時）における生体情報（心拍数、皮膚温度）の平均値である。また、図５及び図６のグラフ中にプロットされた点が、一人の対象者についての安静時における生体情報（心拍数、皮膚温度）の平均値と非安静時における生体情報（心拍数、皮膚温度）の平均値との関係を表す。図５及び図６からわかるように、安静時における生体情報の平均値と非安静時における生体情報の平均値は、強い相関を有する。
【００４５】
また、図５及び図６において、グラフ中にプロットされた点のうち、菱形で示される点は、「好き嫌い」という気分についての、安静時における生体情報（心拍数、皮膚温度）の平均値と非安静時における生体情報（心拍数、皮膚温度）の平均値との関係を示し、四角で示される点は、「焦燥感」という気分についての、安静時における生体情報の平均値と非安静時における生体情報の平均値との関係を示し、三角で示された点は、「集中」という気分についての、安静時における生体情報の平均値と非安静時における生体情報の平均値との関係を示す。図５及び図６からわかるように、いずれの種類の気分においても、特定の生体情報（例えば心拍数、或いは皮膚温度）に関して、安静時における生体情報の平均値と非安静時における生体情報の平均値との相関は、同じような傾向を示す。
【００４６】
行動状態推定部６６は、対象者の行動状態を推定する。ここにおいて、行動状態とは、対象者の動作態様を表す。一例として、行動状態として、「眠っている」、「座っている」、「立っている」、「歩いている」、「走っている」といった、動作態様を例示することができる。
【００４７】
行動状態推定部６６は、行動状態推定処理を実行することにより、対象者の現在時点における行動状態を推定する。図７は、行動状態推定部６６が実行する行動状態推定処理ルーチンの流れを示すフローチャートである。このルーチンが起動すると、行動状態推定部６６は、まず、図７のＳ２０１にて、現在時点に至るまでの所定の微小時間の間に体動検出センサ３が検出した対象者の加速度信号の現在値を、現在値取得部６１から読み出す。次いで、行動状態推定部６６は、入力した加速度信号データに基づいて、加速度の平均値αを算出する（Ｓ２０２）。
【００４８】
続いて、行動状態推定部６６は、加速度の平均値αが、第一閾値加速度α１未満であるか否かを判断する（Ｓ２０３）。第一閾値加速度α１は、加速度の平均値αがそれ未満であるかそれ以上であるかによって対象者が臥位状態（寝ている状態）であるか座位状態（座っている状態）であるかを区別することができる程度の加速度の閾値として、予め定められる。加速度の平均値αが第一閾値加速度α１未満である場合（Ｓ２０３：Ｙｅｓ）、行動状態推定部６６は、Ｓ２０７に処理を進め、対象者の現在の行動状態を臥位状態に設定する。その後、行動状態推定部６６は、このルーチンを終了する。
【００４９】
また、Ｓ２０３にて、加速度の平均値αが、第一閾値加速度α１以上であると判断した場合（Ｓ２０３：Ｎｏ）、行動状態推定部６６は、Ｓ２０４に処理を進め、加速度の平均値αが、第一閾値加速度α１よりも大きい第二閾値加速度α２未満であるか否かを判断する。第二閾値加速度α２は、加速度の平均値αがそれ未満であるかそれ以上であるかによって対象者が座位状態であるか起立状態（立っている状態）であるかを区別することができる程度の加速度の閾値として、予め定められる。加速度の平均値αが第二閾値加速度α２未満である場合（Ｓ２０４：Ｙｅｓ）、行動状態推定部６６は、Ｓ２０８に処理を進め、対象者の行動状態を座位状態に設定する。その後、行動状態推定部６６は、このルーチンを終了する。
【００５０】
また、Ｓ２０４にて、加速度の平均値αが、第二閾値加速度α２以上であると判断した場合（Ｓ２０４：Ｎｏ）、行動状態推定部６６は、Ｓ２０５に処理を進め、加速度の平均値αが、第二閾値加速度α２よりも大きい第三閾値加速度α３未満であるか否かを判断する。第三閾値加速度α３は、加速度の平均値αがそれ未満であるかそれ以上であるかによって対象者が起立状態であるか歩行状態（歩いている状態）であるかを区別することができる程度の加速度の閾値として、予め定められる。加速度の平均値αが第三閾値加速度α３未満である場合（Ｓ２０５：Ｙｅｓ）、行動状態推定部６６は、Ｓ２０９に処理を進め、対象者の行動状態を起立状態に設定する。その後、行動状態推定部６６は、このルーチンを終了する。
【００５１】
また、Ｓ２０５にて、加速度の平均値αが、第三閾値加速度α３以上であると判断した場合（Ｓ２０５：Ｎｏ）、行動状態推定部６６は、Ｓ２０６に処理を進め、加速度の平均値αが、第三閾値加速度α３よりも大きい第四閾値加速度α４未満であるか否かを判断する。第四閾値加速度α４は、加速度の平均値αがそれ未満であるかそれ以上であるかによって対象者が歩行状態であるか走行状態（走っている状態）であるかを区別することができる程度の加速度の閾値として、予め定められる。加速度の平均値αが第四閾値加速度α４未満である場合（Ｓ２０６：Ｙｅｓ）、行動状態推定部６６は、Ｓ２１０に処理を進め、対象者の行動状態を歩行状態に設定する。その後、行動状態推定部６６は、このルーチンを終了する。
【００５２】
また、Ｓ２０６にて、加速度の平均値αが、第四閾値加速度α４以上であると判断した場合（Ｓ２０６：Ｎｏ）、行動状態推定部６６は、Ｓ２１１に処理を進める。そして、Ｓ２１１にて、対象者の行動状態を、走行状態に設定する。その後、行動状態推定部６６は、このルーチンを終了する。
【００５３】
行動状態推定部６６が上記した行動状態判定処理を実行して、対象者の体動（加速度）に基づいて対象者の行動状態を、臥位状態、座位状態、起立状態、歩行状態、走行状態、のいずれかに設定することにより、対象者の現在時点における行動状態が推定される。なお、閾値加速度を更に細分化することにより、より詳細に対象者の行動状態を推定することもできる。例えば、対象者が座位状態で筆記作業をしている状態、対象者が歩行状態で軽作業をしている状態等も、推定することができる。
【００５４】
気分変動量推定部６７は、例えば外部の刺激により変動した気分の種類を特定し、特定した気分の変動量を推定する。気分変動量推定部６７は、変動した気分の種類の特定及びその変動量を推定するにあたり、気分変動量推定処理を実行する。図８Ａ、図８Ｂは、気分変動量推定部６７が実行する気分変動量推定処理ルーチンの流れを示すフローチャートである。
【００５５】
図８Ａ，図８Ｂに示す気分変動量推定処理ルーチンが起動すると、気分変動量推定部６７は、まず、図８ＡのＳ３０１にて、気分状態判定部６２によって対象者の現在の気分状態が非安静状態であると判断されているか否かを判断する。対象者の現在の気分状態が安静状態であると判断されている場合（Ｓ３０１：Ｎｏ）、気分変動量推定部６７はこのルーチンを終了する。
【００５６】
一方、対象者の気分状態が非安静状態であると判断されている場合（Ｓ３０１：Ｙｅｓ）、気分変動量推定部６７は、Ｓ３０２に処理を進めて、現在値取得部６１が取得した対象者の生体情報現在値βｎｏｗを読み出す。
【００５７】
次いで、気分変動量推定部６７は、対象者の安静時における生体情報βの統計値βｓｔを、統計解析部６４から読み出す。ここで、統計値βｓｔは、対象者の安静時における生体情報βの平均値又は分散である。
【００５８】
続いて、気分変動量推定部６７は、Ｘ成分が、対象者が安静状態であるときにおける生体情報βの統計値βｓｔであり、Ｙ成分が生体情報現在値βｎｏｗである座標点Ａ（βｓｔ、βｎｏｗ）を、設定する（Ｓ３０４）。ここで、生体情報現在値βｎｏｗが対象者の生体信号（例えば心拍数）の現在値である場合、統計値βｓｔは、対象者の安静時における生体信号（例えば心拍数）の平均値である。一方、生体情報現在値βｎｏｗが対象者の生体信号（例えばＲＲＩ）の分散（例えばＲＲＩ分散）の現在値である場合、統計値βｓｔは、対象者の安静時における生体信号の分散である。ここで、上述したように、統計解析部６４が演算する安静時生体情報の分散は、現在値取得部６１が取得する安静時における生体信号の分散の現在値の平均値と言える。つまり、Ｙ成分が分散の現在値である場合でも、座標点ＡのＸ成分は、対象者の安静時における生体情報の平均値（生体情報が分散である場合は分散の平均値）である。言い換えれば、座標点Ａは、Ｘ成分が対象者の安静時における生体情報の平均値であり、Ｙ成分が対象者の生体情報現在値である。
【００５９】
次に、気分変動量推定部６７は、安静時における生体情報βの統計値と非安静時における生体情報βの統計値との相関式Ｆ（Ｘ）を、相関関係導出部６５から読み出す（Ｓ３０５）。ここで、生体情報現在値βｎｏｗが対象者の生体信号の現在値である場合、読み出される相関式は、安静時における生体情報βの平均値と非安静時における生体情報βの平均値との相関関係を表す相関式である。一方、生体情報現在値βｎｏｗが対象者の生体信号の分散（例えばＲＲＩ分散）の現在値である場合、読み出される相関式は、安静時における生体情報βの分散と非安静時における生体情報βの分散との相関関係を表す相関式である。相関式Ｆ（Ｘ）は、以下の（５）式のように一次関数として表されるものとする。
【数５】
2019017946.tif 000007
ここで、Ｘは、安静時における生体情報の統計値（平均又は分散）であり、Ｆ（Ｘ）は非安静時における生体情報βの統計値（平均又は分散）である。
【００６０】
続いて、気分変動量推定部６７は、生体情報現在値βｎｏｗが、Ｆ（βｓｔ）よりも大きいか否かを判断する（Ｓ３０６）。生体情報現在値βｎｏｗがＦ（βｓｔ）よりも大きい場合（Ｓ３０６：Ｙｅｓ）、気分変動量推定部６７は、Ｓ３０８に処理を進めて、相関式Ｆ（Ｘ）により表される相関線から生体情報現在値βｎｏｗが乖離している方向（生体情報現在値βｎｏｗの乖離方向）をプラスに設定する。一方、生体情報現在値βｎｏｗがＦ（βｓｔ）以下である場合（Ｓ３０６：Ｎｏ）、気分変動量推定部６７は、Ｓ３０７に処理を進めて、生体情報現在値βｎｏｗがＦ（βｓｔ）よりも小さいか否かを判断する。生体情報現在値βｎｏｗがＦ（βｓｔ）よりも小さい場合（Ｓ３０７：Ｙｅｓ）、気分変動量推定部６７は、Ｓ３０９に処理を進めて、生体情報現在値βｎｏｗの乖離方向をマイナスに設定する。
【００６１】
図９は、Ｘ軸（横軸）が安静時における生体情報βでありＹ軸（縦軸）が非安静時における生体情報βであるＸ−Ｙ平面上に、安静時における生体情報βの統計値と非安静時における生体情報βの統計値との関係を表す相関式Ｆ（Ｘ）により表される相関線と、座標点Ａとが示されたグラフの一例である。図９に示すように、座標点Ａが相関式Ｆ（Ｘ）を表す直線よりも上方に位置している場合、生体情報現在値βｎｏｗがＦ（βｓｔ）よりも大きい。従って、この場合、生体情報現在値βｎｏｗの乖離方向がプラスに設定される。一方、座標点Ａが相関式Ｆ（Ｘ）を表す直線よりも下側に位置している場合、生体情報現在値βｎｏｗがＦ（βｓｔ）よりも小さい。従って、この場合、生体情報現在値βｎｏｗの乖離方向がマイナスに設定される。このように、座標点Ａが相関式Ｆ（Ｘ）を表す相関線よりも上側に位置するか下側に位置するかによって、乖離方向がプラスであるかマイナスであるかが決定される。
【００６２】
なお、図８のＳ３０７の判定結果がＮｏである場合、生体情報現在値βｎｏｗがＦ（βｓｔ）に等しい。この場合、気分変動量推定部６７は、このルーチンを終了する。
【００６３】
Ｓ３０８またはＳ３０９にて、生体情報現在値βｎｏｗの乖離方向をプラス或いはマイナスに設定した後に、気分変動量推定部６７は、図８ＢのＳ３１０に処理を進める。Ｓ３１０では、気分変動量推定部６７は、行動状態推定部６６が推定（設定）した対象者の現在時点における行動状態を読み出す。次いで、気分変動量推定部６７は、対象者の生体情報現在値βｎｏｗ及び対象者の現在時点における行動状態に基づいて、対象者の気分のうち変動している気分を推定（特定）する（Ｓ３１１）。
【００６４】
気分の変動は、上記（４）式からわかるように、特定の生体情報に影響する。また、気分の種類によって、気分の変動時における生体情報の変化の方向が異なる。例えば、「焦燥感」という気分と、「集中」という気分について言えば、一般的には、「焦燥感」が大きい場合には心拍数が増加するが、「集中」が大きい（集中力が高い）場合には心拍数が低下する。従って、対象者の現在時点における心拍数（心拍数現在値）βｎｏｗが、対象者の非安静時における心拍数の平均値（Ｆ（βｓｔ））よりも大きい場合（つまり心拍数現在値βｎｏｗの乖離方向がプラスである場合）、「焦燥感」という気分が変動している確率が高い。一方、対象者の心拍数現在値βｎｏｗが、対象者の非安静時における心拍数の平均値（Ｆ（βｓｔ））未満である場合（つまり心拍数現在値βｎｏｗの乖離方向がマイナスである場合）、「集中」という気分が変動している（すなわち集中している）確率が高い。
【００６５】
また、複数種類の気分について、それらの気分の変動時に、ある生体情報現在値βｎｏｗの乖離方向が同じである場合も存在する。この場合、対象者の生体情報現在値βｎｏｗの乖離方向がプラス方向（又はマイナス方向）であっても、それが、どの気分の変動を表しているのか区別できない可能性がある。このような場合に、対象者の現在の行動状態を用いて、どのような気分が変動しているのかが推定される。例えば、生体情報現在値βｎｏｗの乖離方向から、「焦燥感」又は「好き嫌い」という気分が変動していると判断される場合において、対象者の現在の行動状態が、「臥位状態」或いは「座位状態」であれば、例えば対象者が座った状態で不快な映像を視聴していることにより、視聴しているコンテンツに対して嫌悪感を抱いているというようなことが想定され得る。この場合、「好き嫌い」という気分が変動していると推定される。これに対し、対象者の現在の行動状態が、「歩行状態」或いは「走行状態」であれば、例えば対象者が焦っていてじっとしていられないということが想定され得る。この場合、「焦燥感」という気分が変動していると推定される。このように、気分変動量推定部６７は、Ｓ３１１にて、対象者の生体情報現在値βｎｏｗの乖離方向及び対象者の現在の行動状態に基づいて、変動している気分を推定（特定）する。
【００６６】
次に、気分変動量推定部６７は、Ｓ３１２に処理を進める。そして、Ｓ３１２以降の処理を実行することにより、変動している気分の変動量を推定する。この場合、気分変動量推定部６７は、Ｓ３１２にて、図９に示すような、Ｘ軸が安静時における生体情報βでありＹ軸が非安静時における生体情報βであるＸ−Ｙ平面上において、相関式Ｆ（Ｘ）により表される相関線と座標点Ａとの距離を、乖離距離Ｌとして算出する。乖離距離Ｌは、相関式Ｆ（Ｘ）により表される相関線と座標点Ａとの位置関係から算出することができる。
【００６７】
乖離距離Ｌは、相関式Ｆ（Ｘ）により表される相関線からの、対象者の生体情報現在値βｎｏｗの乖離量を表す。この乖離距離Ｌが大きいほど、対象者の生体情報現在値βｎｏｗは通常の気分の変動よりも大きく変動していると言え、それ故に、対象者の気分の変動量は大きいと言える。また、この乖離距離Ｌが小さい場合、対象者の生体情報現在値βｎｏｗは通常の気分の変動と同程度に変動していると言え、それ故に、対象者の気分の変動量は小さいと言える。つまり、乖離距離Ｌは、変動している気分の変動量を表す。従って、乖離距離Ｌに基づいて、気分の変動量を推定することができる。
【００６８】
なお、発明者の研究では、生体情報のうち、ＲＲＩ分散についての乖離距離Ｌが、変動している気分の変動量をより精度良く表すことが実験により確かめられた。従って、乖離距離Ｌは、Ｘ軸が安静状態であるときにおけるＲＲＩ分散でありＹ軸が非安静状態であるときにおけるＲＲＩ分散であるＸ−Ｙ平面上において、Ｘ成分が安静状態であるときにおける対象者のＲＲＩ分散でありＹ成分が現在時点における対象者のＲＲＩ分散である座標点Ａと、安静状態であるときにおけるＲＲＩ分散と非安静状態であるときにおけるＲＲＩ分散との相関式により表される相関線との距離であるとよい。つまり、気分変動量推定部６７が変動している気分の変動量を推定する場合に用いる生体情報は、ＲＲＩ分散であるのがよい。
【００６９】
気分変動量推定部６７は、Ｓ３１２にて乖離距離Ｌを算出した後に、Ｓ３１３に処理を進めて、乖離距離Ｌに基づいて気分変動量Ｅを算出する。この場合、乖離距離Ｌが大きいほど気分変動量Ｅが大きくなるように、気分変動量Ｅが算出される。気分変動量Ｅは、距離Ｌに比例する値として算出されていてもよい。或は、気分変動量Ｅは、最も気分の変動が大きいと考えられる場合に得られる距離Ｌｍａｘに対する、Ｓ３１２で算出された距離Ｌの比率の百分率として算出されていてもよい。この場合、距離Ｌｍａｘは、全ての対象者に対して得られる最大距離であってもよいし、或は、その対象者個人に対して得られる最大距離であってもよい。Ｓ３１３にて気分変動量Ｅを算出した後に、気分変動量推定部６７はこのルーチンを終了する。
【００７０】
気分変動量推定部６７が上記した気分推定処理を実行することにより、対象者についての変動している気分を特定し且つ特定した気分の変動量（気分変動量Ｅ）を推定することができる。なお、変動している気分は、気分状態判定部６２がＳ１０３にて演算する各気分についての発生確率ｐ及び、対象者の現在の行動状態に基づいても、推定することができる。
【００７１】
ストレス推定部６８は、ストレス推定処理を実行することにより、対象者のストレスの強弱及び心理的負荷の大小を推定する。図１０Ａ及び図１０Ｂは、ストレス推定部６８が実行するストレス推定処理ルーチンの流れを示すフローチャートである。このルーチンが起動すると、ストレス推定部６８は、まず、図１０ＡのＳ４０１にて、安静時における対象者のＲＲＩの平均値及び分散を、統計解析部６４から読み出す。次いで、ストレス推定部６８は、Ｓ４０１にて読み出した平均値及び分散により表される正規分布関数に対して乱数列を用いた演算を実行して、安静時におけるＲＲＩの正規分布ダミーデータを作成する（Ｓ４０２）。
【００７２】
続いて、ストレス推定部６８は、安静時におけるＲＲＩの平均値と非安静時におけるＲＲＩの平均値との相関式を、相関関係導出部６５から読み出し（Ｓ４０３）、読み出した相関式に安静時における対象者のＲＲＩの平均値を代入することにより、非安静時における対象者のＲＲＩの平均値を算出する（Ｓ４０４）。次に、ストレス推定部６８は、安静時におけるＲＲＩの分散と非安静時におけるＲＲＩの分散との相関式を相関関係導出部６５から読み出し（Ｓ４０５）、読み出した相関式に安静時における対象者のＲＲＩの分散を代入することにより、非安静時における対象者のＲＲＩの分散を算出する（Ｓ４０６）。続いて、ストレス推定部６８は、Ｓ４０４にて算出した非安静時における対象者のＲＲＩの平均値及びＳ４０６にて算出した非安静時における対象者のＲＲＩの分散により表される正規分布関数に対して乱数列を用いた演算を実行して、非安静時におけるＲＲＩの正規分布ダミーデータを作成する（Ｓ４０７）。以上の処理により、安静時におけるＲＲＩのダミーデータ及び非安静時におけるＲＲＩのダミーデータが作成される。
【００７３】
次に、ストレス推定部６８は、安静時におけるＲＲＩのダミーデータを用いて、安静時におけるＲＲＩのローレンツプロット（安静時ローレンツプロット）を作成する（Ｓ４０８）。ローレンツプロットは、Ｘ成分が乱数列のｎ番目に求められたＲＲＩデータ（ＲＲＩ（ｎ））でありＹ成分が乱数列の（ｎ＋１）番目に求められたＲＲＩデータ（ＲＲＩ（ｎ＋１））である座標点（ＲＲＩ（ｎ），ＲＲＩ（ｎ＋１））を、Ｘ−Ｙ平面上にプロットすることにより作成される。
【００７４】
図１１に、ローレンツプロットの一例を示す。図１１に示すように、ローレンツプロットを構成する点群により表される図形は、一般的に、Ｙ＝Ｘの直線を中心として広がりを持つ形状である。ローレンツプロットを構成する点群により表される形状の面積は、ＬＰ面積と呼ばれる。また、ローレンツプロットを構成する各点をＹ＝Ｘ軸とＹ＝−Ｘ軸に投影後、Ｙ＝Ｘ軸において、原点（０，０）からの距離の平均値は、Ｍ値と呼ばれる。ＬＰ面積は、外部の刺激に対する人の反応に応じて変化する。例えば、ＬＰ面積が大きい人は、変化や刺激に対して感受性が高く、それによる気分の変動を自分で抑制し難いタイプに分類される。また、ＬＰ面積が小さい人は、変化や刺激に対して自分で気分の変動を抑えることができるタイプに分類される。従って、ＬＰ面積に基づいて、人のタイプ分けを実施することができる。
【００７５】
ストレス推定部６８は、上記のようにして安静時ローレンツプロットを作成した後に、安静時ローレンツプロットのＬＰ面積（ＬＰ１）及びＭ値（Ｍ１）を算出する（Ｓ４０９）。なお、ＬＰ面積は、ローレンツプロットを表す画像から得ることもできるし、また、ローレンツプロットの形状が略楕円であるような場合には、楕円の面積の公式（長径×短径×円周率／４）から求めることもできる。Ｍ値は、ローレンツプロットの形状が略楕円であるような場合にはその楕円の中心座標の原点からの距離として求めることもできるし、或いは、ローレンツプロットを構成する各点をＹ＝Ｘ軸とＹ＝−Ｘ軸に投影後、Ｙ＝Ｘ軸において、原点（０，０）からの距離の平均値により求めることもできる。
【００７６】
次いで、ストレス推定部６８は、非安静時におけるＲＲＩの正規分布ダミーデータを用いて、非安静時におけるローレンツプロット（非安静時ローレンツプロット）を作成し（Ｓ４１０）、非安静時ローレンツプロットについてのＬＰ面積（ＬＰ２）及びＭ値（Ｍ２）を算出する（Ｓ４１１）。
【００７７】
次に、ストレス推定部６８は、安静時ローレンツプロットのＭ値（Ｍ１）が、非安静時ローレンツプロットのＭ値（Ｍ２）よりも小さいか否かを判断する（Ｓ４１２）。Ｍ１がＭ２未満である場合（Ｓ４１２：Ｙｅｓ）、ストレス推定部６８は、気分変動時における対象者の心理的なストレスが小さいと判断する。この場合、ストレス推定部６８は、ストレス量Ｓを、「小」に設定する（Ｓ４１３）。一方、Ｍ１がＭ２以上である場合（Ｓ４１２：Ｎｏ）、ストレス推定部６８は、気分変動時における対象者の心理的なストレスが大きいと判断する。この場合、ストレス推定部６８は、ストレス量Ｓを「大」に設定する（Ｓ４１４）。このように、安静時と非安静時におけるローレンツプロットのＭ値の大きさの比較によって、心理的なストレスの大小が判断される。
【００７８】
Ｓ４１３又はＳ４１４にてストレス量Ｓを設定した後に、ストレス推定部６８は、Ｓ４１５に処理を進める。Ｓ４１５では、ストレス推定部６８は、安静時ローレンツプロットのＬＰ面積（ＬＰ１）が、非安静時ローレンツプロットのＬＰ面積（ＬＰ２）未満であるか否かを判断する（Ｓ４１５）。ＬＰ１がＬＰ２未満である場合（Ｓ４１５：Ｙｅｓ）、ストレス推定部６８は、気分変動時における対象者のストレスの変動が大きく、心理的な負担が大きいと判断する。この場合、ストレス推定部６８は、心理的負担Ｖを「大」に設定する（Ｓ４１６）。その後、ストレス推定部６８はこのルーチンを終了する。一方、ＬＰ１がＬＰ２以上である場合（Ｓ４１５：Ｎｏ）、ストレス推定部６８は、気分変動時における対象者のストレスの変動が小さく、心理的な負担が小さいと判断する。この場合、ストレス推定部６８は、心理的負担Ｖを「小」に設定する（Ｓ４１７）。その後、ストレス推定部６８はこのルーチンを終了する。
【００７９】
このように、ストレス推定部６８は、ストレス推定処理を実行して、安静時におけるＲＲＩのローレンツプロットのＬＰ面積及びＭ値と、非安静時におけるＲＲＩのローレンツプロットのＬＰ面積及びＭ値とを比較する。そして、その比較の結果に応じて、対象者の気分が変動した場合におけるストレスの強弱及び心理的負担の大小を推定する。
【００８０】
統合処理部６９は、上記した各機能が各処理を実行した結果に基づいて、各結果の関連付けを行う。この場合、統合処理部６９は、統合処理を実行する。
【００８１】
図１２は、統合処理部６９が実行する統合処理ルーチンの流れを示すフローチャートである。このルーチンが起動すると、統合処理部６９は、まず、図１２のＳ５０１にて、行動状態推定部６６が推定した対象者の現在時点における行動状態を取得する。次いで、統合処理部６９は、気分変動量推定部６７が特定した変動している気分、及び、変動している気分の変動量（気分変動量）Ｅを取得する（Ｓ５０２）。続いて、統合処理部６９は、ストレス推定部６８が推定したストレス量Ｓ及び心理的負担Ｖを取得する（Ｓ５０３）。
【００８２】
その後、統合処理部６９は、Ｓ５０４にて、関連付け処理を実施する。この関連付け処理では、統合処理部６９は、Ｓ５０１にて取得した行動状態、Ｓ５０２にて取得した気分及び気分変動量Ｅ、Ｓ５０３にて取得したストレス量Ｓ及び心理的負担Ｖ、に基づいて、対象者がどのような行動状態であり、且つ、そのときにおける対象者の気分、及び、対象者が感じているストレス、心理的負担についての関連付けを行う。
【００８３】
例えば、対象者が、座っている状態で不快な映像を視聴している場面を想定する。この場合、対象者の行動状態は座位状態である。また、このとき、「嫌悪感」という気分が変動し、その気分の変動量が大きく、さらに、ストレス量Ｓ及び心理的負担Ｖが「大」であると推定されたとする。この場合、対象者は、不快な映像の視聴により「嫌悪感」を強く感じ、且つ、それによりストレス及び心理的負担が大きいというように、対象者の気分の状態とストレスの強弱及び心理的負担の大小を、関連付けることができる。
【００８４】
また、対象者が、忘れ物をして自宅に走って戻っている状態を想定する。この場合、対象者の行動状態は走行状態である。また、このとき、「焦燥感」という気分が変動し、気分変動量Ｅが大きいが、ストレス量Ｓ及び心理的負担Ｖは「小」であると推定されたとする。この場合、対象者は、忘れ物をしたことにより「焦り」を強く感じているものの、ストレス及び心理的負担はさほど大きくないといった、対象者の気分の状態とストレス及び心理的負担の大小を、関連付けることができる。
【００８５】
統合処理部６９は、Ｓ５０４にて、関連付け処理を実施した後に、その結果を出力する。出力された結果は、サーバ６から対象者の情報端末４に送信される。これにより、対象者の情報端末４に、その対象者の現在の行動状態、気分状態、ストレス、心理的負担が、関連付けて表示される。対象者は、こうして表示された情報に基づいて、自らの心理状態を把握することができる。なお、統合処理部６９は、Ｓ５０４の処理によって対象者の心理的負担が大きいと判断する場合には、Ｓ５０５にて、注意を喚起するような情報（例えば、「心理的負担が増大しています」というような文字情報）、或いは、心理的負担を軽減するようなアドバイス情報（例えば、「気分転換したらどうですか」というような文字情報）を合わせて送信することができる。その後、統合処理部６９は、このルーチンを終了する。
【００８６】
統合処理部６９が上記した統合判断処理を実行することにより、対象者の現在時点における行動状態、変動している気分及びその気分についての気分変動量、ストレス及び心理的負担が、関連付けられる。そして、関連付けられたこれらの情報から、対象者がどのような行動状態であるときにどの程度の気分の変動が発生しているか、及び、その時点において対象者が心理的に負担を感じているか、といったことを、推定することができる。
【００８７】
環境制御部７０は、環境制御処理を実行することにより、統合処理部６９が関連付けた結果（行動状態、気分状態、ストレスの強弱、心理的負担の大小の関連付け）に基づいて、対象者の周辺の環境を制御する。図１３は、環境制御部７０が実行する環境制御処理ルーチンの流れを示すフローチャートである。図１３に示すルーチンが起動すると、環境制御部７０は、まず図１３のＳ６０１にて、対象者について統合処理部６９がＳ５０４にて判断した結果（行動状態、気分状態及びその気分の変動量（気分変動量Ｅ）、ストレス量Ｓ、心理的負担Ｖの関連付け）を入力する。次いで、環境制御部７０は、Ｓ６０２にて、対象者の気分の変動量（気分変動量Ｅ）が閾値変動量Ｅｔｈよりも大きいか否かを判断する。閾値変動量Ｅｔｈは、対象者の気分変動量Ｅがそれよりも大きい場合に生理的になんらかの影響を及ぼす可能性がある値として予め定められる。
【００８８】
気分変動量Ｅが閾値変動量Ｅｔｈ以下である場合（Ｓ６０２：Ｎｏ）、環境制御部７０はこのルーチンを終了する。一方、気分変動量Ｅが閾値変動量Ｅｔｈよりも大きい場合（Ｓ６０２：Ｙｅｓ）、環境制御部７０はＳ６０３に処理を進める。Ｓ６０３では、環境制御部７０は、対象者の周囲環境を制御することができる環境制御機器７のうち、対象者の気分の変動量を小さくすることができるような環境制御機器を選択する。
【００８９】
例えば、対象者が「緊張」しているような気分状態である場合（「緊張」という気分の変動量が大きい場合）、心身を落ち着かせるような音楽や映像を視聴させることにより、緊張状態が緩和する。従って、この場合、環境制御部７０は、Ｓ６０３にて、楽曲或は映像を再生することができる機器を選択する。
【００９０】
次に、環境制御部７０は、対象者の心理的負担Ｖが「大」に設定されているか否かを判断する（Ｓ６０４）。心理的負担Ｖが「大」に設定されていない場合（Ｓ６０４：Ｎｏ）、環境制御部７０は、このルーチンを終了する。一方、心理的負担Ｖが「大」に設定されている場合（Ｓ６０４：Ｙｅｓ）、環境制御部７０はＳ６０５に処理を進めて、Ｓ６０３にて選択した環境制御機器を動作させるため指令信号を中継器５に送信する。中継器５は、環境制御部７０（サーバ６）から指令信号を入力すると、環境制御部７０がＳ６０３にて選択した環境制御機器が動作するように、環境制御機器に動作指令信号を出力する。これにより選択された環境制御機器が動作する。その後、環境制御部７０はこのルーチンを終了する。
【００９１】
環境制御部７０が上記した環境制御処理を実行した場合、対象者の気分が大きく変動していて、且つ、心理的負担が大きいときに、変動した気分の変動量を小さくして心理的負担を軽減するように環境制御機器が動作する。これにより、対象者の気分状態を安静状態にさせることができる。
【００９２】
以上、本発明の実施形態について説明したが、本発明は、上記実施形態に限定されるべきものではない。例えば、上記実施形態では、統合処理部６９が統合処理の結果を対象者の情報端末４に自動的に送信する例を示したが、対象者が情報端末４から送信要求をサーバ６に送信したときに、統合処理部６９が統合処理の結果を情報端末４に送信するように構成してもよい。つまり、対象者の要求に応じて統合処理の結果が対象者に送信されるように構成されていてもよい。また、本発明に係る気分推定システムは、対象者の周囲の環境状況（温度、照度等）を検出する環境状況検出センサを備えていてもよい。この場合、環境制御部７０は、環境制御処理にて、環境状況をも考慮して環境制御機器を選択することができる。また、本発明に係る気分推定システムは、ストレス推定部が作成する各対象者のローレンツプロットに基づいて、或いは各対象者のＲＲＩ分散に基づいて、安静時における生体情報と非安静時における生体情報との関連性に関する対象者のタイプ分けを実施してもよい。こうしてタイプ分けをして各対象者の安静時及び非安静時の生体情報の関連性の傾向を把握しておくことにより、より正確に、個々の対象者の気分の変動量を精度良く推定することができる。また、本発明に係る気分推定システムは、対象者の安静時における気分の変動状態と非安静時における気分の変動状態とを過去の履歴から学習する学習機能を備えていても良い。これによれば、学習機能を備えることにより、対象者の個々に対して精度よく、気分の変動量を推定することができる。このように、本発明は、その趣旨を逸脱しない限りにおいて変形可能である。
【符号の説明】
【００９３】
１…気分推定システム、２…生体信号検出センサ、３…体動検出センサ、４…情報端末、５…中継器、６…サーバ、６１…現在値取得部（生体情報現在値取得部）、６２…気分状態判定部、６３…情報記憶部、６４…統計解析部、６５…相関関係導出部、６６…行動状態推定部、６７…気分変動量推定部、６８…ストレス推定部、６９…統合処理部、７０…環境制御部、７…環境制御機器、Ｅ…気分変動量、Ｌ…乖離距離、Ｓ…ストレス量、Ｖ…心理的負担
【図１】
2019017946.tif 000008
【図２】
2019017946.tif 000009
【図３】
2019017946.tif 000010
【図４】
2019017946.tif 000011
【図５】
2019017946.tif 000012
【図６】
2019017946.tif 000013
【図７】
2019017946.tif 000014
【図８Ａ】
2019017946.tif 000015
【図８Ｂ】
2019017946.tif 000016
【図９】
2019017946.tif 000017
【図１０Ａ】
2019017946.tif 000018
【図１０Ｂ】
2019017946.tif 000019
【図１１】
2019017946.tif 000020
【図１２】
2019017946.tif 000021
【図１３】
2019017946.tif 000022