(57)【要約】【課題】生体の体内時計の指示値（例えば時刻）を実時間的かつ高精度に特定する。【解決手段】第１生体情報および第２生体情報の各々の周期的な変動と時刻との関係を表す基準リズム情報を利用して、生体における前記第１生体情報および前記第２生体情報の各々の測定値から、前記生体の体内時計の指示値を特定する解析処理部を具備する生体解析装置。【選択図】図２
2019055087.tif 000002
【特許請求の範囲】
【請求項１】
第１生体情報および第２生体情報の各々の周期的な変動と時刻との関係を表す基準リズム情報を利用して、生体における前記第１生体情報および前記第２生体情報の各々の測定値から、前記生体の体内時計の指示値を特定する解析処理部
を具備する生体解析装置。
【請求項２】
前記基準リズム情報は、前記第１生体情報に対応する第１軸と前記第２生体情報に対応する第２軸とが設定された座標空間において前記第１生体情報と前記第２生体情報との関係を表す環状の基準リズム曲線を規定し、
前記解析処理部は、前記座標空間内において前記第１生体情報および前記第２生体情報の各々の測定値が示す測定点と、前記基準リズム曲線の中心点とを通過する直線が、前記基準リズム曲線に交差する地点に対応する時刻に応じて、前記体内時計の指示値を特定する
請求項１の生体解析装置。
【請求項３】
前記基準リズム情報は、前記第１生体情報に対応する第１軸と前記第２生体情報に対応する第２軸とが設定された座標空間において前記第１生体情報と前記第２生体情報との関係を表す環状の基準リズム曲線を規定し、
前記解析処理部は、前記座標空間内において前記第１生体情報および前記第２生体情報の各々の測定値が示す測定点の位置ベクトルの角度に応じて、前記体内時計の指示値を特定する
請求項１の生体解析装置。
【請求項４】
前記第１生体情報と前記第２生体情報とは、時間的な変動の位相が相違する
請求項１から請求項３の何れかの生体解析装置。
【請求項５】
前記第１生体情報は、前記生体の脈拍数、酸素飽和度および血圧の何れかであり、
前記第２生体情報は、前記生体の深部温度、皮膚温度、脈拍間隔および脈拍変動の何れかである
請求項４の生体解析装置。
【請求項６】
第１生体情報および第２生体情報の各々の周期的な変動と時刻との関係を表す基準リズム情報を利用して、生体における前記第１生体情報および前記第２生体情報の各々の測定値から、前記生体の体内時計の指示値を特定する
生体解析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、生体の状態を解析するための技術に関する。
【背景技術】
【０００２】
脈拍数や体温等の各種の生体情報は、例えば日単位または月単位で周期的に変動する。特許文献１から特許文献５には、生体情報の周期的な変動（すなわち生体リズム）を解析する技術が提案されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１７−１８１９４号公報
【特許文献２】特開２０１５−２９６０５号公報
【特許文献３】特開２０１２−２４３３２号公報
【特許文献４】特開２０１０−１５８２６７号公報
【特許文献５】特開平５−３８７８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
生体リズムは体内時計（生物時計）により調整される。したがって、体内時計の時刻は、生体状態（例えば健康状態や疾病状態）を評価するうえで重要な指標である。しかし、特許文献１から特許文献５の技術は、生体リズムを解析する技術に過ぎず、体内時計の時刻までは特定できない。
【０００５】
また、特許文献１から特許文献５の技術では、生体情報の複数の測定値の時系列を三角関数により近似することで生体リズムを特定する。したがって、測定時点に近い短時間内の測定値のみからでは生体リズムを正確に特定できない。１種類の生体情報のみからでは、実際には生体リズムを正確に特定できないという問題もある。以上の事情を考慮して、本発明は、生体の体内時計の指示値（例えば時刻）を実時間的かつ高精度に特定することを目的とする。
【課題を解決するための手段】
【０００６】
以上の課題を解決するために、本発明の好適な態様に係る生体解析装置は、第１生体情報および第２生体情報の各々の周期的な変動と時刻との関係を表す基準リズム情報を利用して、生体における前記第１生体情報および前記第２生体情報の各々の測定値から、前記生体の体内時計の指示値を特定する解析処理部を具備する。以上の態様では、第１生体情報および第２生体情報の各々の周期的な変動と時刻との関係を表す基準リズム情報を利用して、第１生体情報および第２生体情報の各々の測定値から体内時計の指示値が特定される。したがって、体内時計の時刻を実時間的かつ高精度に特定することが可能である。
【０００７】
本発明の好適な態様において、前記基準リズム情報は、前記第１生体情報に対応する第１軸と前記第２生体情報に対応する第２軸とが設定された座標空間において前記第１生体情報と前記第２生体情報との関係を表す環状の基準リズム曲線を規定し、前記解析処理部は、前記座標空間内において前記第１生体情報および前記第２生体情報の各々の測定値が示す測定点と、前記基準リズム曲線の中心点とを通過する直線が、前記基準リズム曲線に交差する地点に対応する時刻に応じて、前記体内時計の指示値を特定する。以上の態様では、座標空間内の測定点と基準リズム曲線の中心点とを通過する直線が基準リズム曲線に交差する地点に対応する時刻に応じて、体内時計の指示値が特定される。したがって、体内時計の指示値を簡便な処理により特定できるという利点がある。
【０００８】
本発明の好適な態様において、前記基準リズム情報は、前記第１生体情報に対応する第１軸と前記第２生体情報に対応する第２軸とが設定された座標空間において前記第１生体情報と前記第２生体情報との関係を表す環状の基準リズム曲線を規定し、前記解析処理部は、前記座標空間内において前記第１生体情報および前記第２生体情報の各々の測定値が示す測定点の位置ベクトルの角度に応じて、前記体内時計の指示値を特定する。以上の態様では、座標空間内の測定点と基準リズム曲線の中心点とを通過する直線の角度に応じて体内時計の指示値が特定される。したがって、体内時計の指示値を簡便な処理により特定できるという利点がある。
【０００９】
本発明の好適な態様において、前記第１生体情報と前記第２生体情報とは、時間的な変動の位相が相違する。以上の態様では、第１生体情報と第２生体情報とで時間的な変動の位相が相違するから、第１生体情報と第２生体情報とが略同等の位相で変動する場合と比較して、体内時計の時刻を高精度に特定することが可能である。
【００１０】
本発明の好適な態様において、前記第１生体情報は、前記生体の脈拍数、酸素飽和度および血圧の何れかであり、前記第２生体情報は、前記生体の深部温度、皮膚温度、脈拍間隔および脈拍変動の何れかである。以上の態様では、第１生体情報の時間的な変動と第２生体情報の時間的な変動との位相差が概略的にはπ／２となる。したがって、第１生体情報と第２生体情報との関係を表す環状の基準リズム曲線を利用して、体内時計の指示値を高精度に特定できるという利点がある。
【００１１】
本発明の好適な態様に係る生体解析方法は、第１生体情報および第２生体情報の各々の周期的な変動と時刻との関係を表す基準リズム情報を利用して、生体における前記第１生体情報および前記第２生体情報の各々の測定値から、前記生体の体内時計の指示値を特定する。以上の態様では、第１生体情報および第２生体情報の各々の周期的な変動と時刻との関係を表す基準リズム情報を利用して、第１生体情報および第２生体情報の各々の測定値から体内時計の指示値が特定される。したがって、体内時計の時刻を実時間的かつ高精度に特定することが可能である。
【図面の簡単な説明】
【００１２】
【図１】本発明の第１実施形態に係る生体解析装置の外観図である。
【図２】生体解析装置の機能に着目した構成図である。
【図３】生体情報の時間的な変化を示すグラフである。
【図４】基準リズム曲線を利用して体内時計時刻を特定する方法の説明図である。
【図５】基準リズム情報の模式図である。
【図６】基準リズム情報の他例の模式図である。
【図７】第２実施形態において基準リズム曲線を利用して体内時計時刻を特定する方法の説明図である。
【図８】第２実施形態における基準リズム情報の模式図である。
【図９】第２実施形態の解析処理部が体内時計時刻を特定する処理のフローチャートである。
【図１０】第３実施形態の解析処理部が体内時計時刻を特定する処理のフローチャートである。
【図１１】第３実施形態において基準リズム曲線を利用して体内時計時刻を特定する方法の説明図である。
【図１２】第３実施形態の効果を示すグラフである。
【図１３】実製品の構成図である。
【図１４】変形例における表示装置の構成図である。
【図１５】変形例における生体解析装置の構成図である。
【図１６】変形例における生体解析装置の構成図である。
【図１７】変形例における生体解析装置の構成図である。
【発明を実施するための形態】
【００１３】
＜第１実施形態＞
図１は、第１実施形態における生体解析装置１００の外観図である。生体解析装置１００は、生体に装着される可搬型の電子機器であり、生体の体内時計の時刻（以下「体内時計時刻」という）を特定する。第１実施形態の体内時計時刻は、日周期（サーカディアンリズム）の体内時計において現在の生体の状態が何れの時点（位相）に相当するかを示す指示値である。第１実施形態では、生体の手首に装着される腕時計型の生体解析装置１００を例示する。ただし、生体解析装置１００が装着される部位は手首に限定されない。例えば、頸部、上腕部、足首、胴体（胸部，腰部）、手足の甲部等、生体の任意の部位に装着可能な形態に生体解析装置１００は構成される。
【００１４】
図１に例示される通り、生体解析装置１００は、時針１４１と分針１４２と指針１４３とを含む表示装置１４を具備する。時針１４１および分針１４２は現在時刻を表示し、指針１４３は生体の体内時計時刻を表示する。体内時計時刻は、例えば時針１４１が示す時間を基準とした相対的な時間として表示される。なお、体内時計時刻を表示する構成は以上の例示に限定されない。例えば、液晶表示パネルや有機ＥＬ（Electroluminescence）表示パネル等の表示器を表示装置１４として体内時計時刻を表示することも可能である。また、酸素飽和度や血圧等の各種の生体情報を表示装置１４が表示してもよい。
【００１５】
図２は、生体解析装置１００の機能的な構成を例示するブロック図である。図２に例示される通り、第１実施形態の生体解析装置１００は、処理装置１１と記憶装置１２と検出装置１３と表示装置１４とを具備する。
【００１６】
検出装置１３は、生体の状態を検出するためのセンサーである。第１実施形態の検出装置１３は、図２に例示される通り、光センサー１３１と熱センサー１３２とを含んで構成される。光センサー１３１は、生体の状態を光学的に検出するセンサーである。例えば、光センサー１３１は、生体に光（例えば近赤外光）を照射する発光素子と、生体から到来する光を受光する受光素子とを具備する。熱センサー１３２は、生体の温度を検出するセンサーであり、例えばサーミスターや熱電対等の測温体を含んで構成される。
【００１７】
処理装置１１は、ＣＰＵ（Central Processing Unit）またはＦＰＧＡ（Field-Programmable Gate Array）等の演算回路であり、生体解析装置１００の全体を制御する。記憶装置１２は、例えば不揮発性の半導体メモリーで構成され、処理装置１１が実行するプログラムと処理装置１１が使用する各種のデータとを記憶する。第１実施形態の処理装置１１は、記憶装置１２に記憶されたプログラムを実行することで、生体の状態を解析するための複数の機能（測定処理部２１および解析処理部２２）を実現する。なお、処理装置１１の機能を複数の集積回路に分散した構成、または、処理装置１１の一部または全部の機能を専用の電子回路で実現した構成も採用され得る。また、図２では処理装置１１と記憶装置１２とを別個の要素として図示したが、記憶装置１２を内包する処理装置１１を例えばＡＳＩＣ（Application Specific Integrated Circuit）等により実現することも可能である。
【００１８】
測定処理部２１は、検出装置１３による検出結果から生体情報を測定する。例えば、光センサー１３１から供給される検出信号について公知の演算処理を実行することで、測定処理部２１は、光電脈波、容積脈波、脈拍数（心拍数）、脈拍間隔（心拍間隔）、脈拍変動（心拍変動）、血流速、血流量、血液灌流量、血管抵抗、血圧（拡張期血圧／収縮期血圧）および酸素飽和度（ＳｐＯ２）等の生体情報を測定することが可能である。また、熱センサー１３２から供給される検出信号について公知の演算処理を実行することで、測定処理部２１は、熱流束（単位面積内の熱流）、皮膚温度および深部温度等の生体情報を測定することが可能である。皮膚温度および深部温度は、例えば手首について測定される。深部温度は、例えば熱流束と皮膚温度とから特定される。
【００１９】
図３は、複数種の生体情報（脈拍数，酸素飽和度，皮膚温度，深部温度，熱流束）の時間的な変化のグラフである。図３には、実測値と近似値とが併記されている。近似値は、公知のコサイナー法を利用して実測値を三角関数で近似した数値である。
【００２０】
図３に例示される通り、各生体情報は、例えば日単位で周期的に変動する。生体情報の時間的な変動の位相は生体情報の種類毎に相違し得る。例えば、脈拍数や熱流束が極大となる時刻（以下「頂点位相時刻」という）は15:30前後であるのに対し、深部温度の頂点位相時刻は概ね21:00〜24:00である。また、酸素飽和度の頂点位相時刻は概ね2:30頃である。すなわち、脈拍数の頂点位相時刻が夕方の時間帯であるのに対し、深部温度の頂点位相時刻は夜間から深夜の時間帯である。
【００２１】
以上のように生体情報の位相が相違する理由としては、交感神経と副交感神経とからなる自律神経の機能の相違が想定される。一般的に脈拍数や熱流束等の代謝に関する機能は日中（15:00頃）に上昇して夜間は低下する。自律神経は代謝機能の変動リズムを制御し、脈拍数が低下する夜間は副交感神経が優位になる。他方、手首等の末端部位の温度は血流により変動する。自律神経は末端部位の血流を制御することで放熱を促進し、体温を調整すると推測される。以上の説明から理解される通り、各生体情報における位相の相違は、心臓等の臓器や深部等の体幹部位に対する自律神経の影響と、皮膚血流等の末端部位に対する自律神経の影響との間に時間差があることが要因であると推測される。
【００２２】
以下の説明では、時間的な変動の位相が相違する関係にある生体情報Ｘ（第１生体情報の例示）と生体情報Ｙ（第２生体情報の例示）とに着目する。具体的には、生体情報Ｘと生体情報Ｙとの位相差はπ/２である。例えば、生体情報Ｘは、脈拍数、酸素飽和度および血圧の何れかであり、生体情報Ｙは、深部温度、皮膚温度、脈拍間隔および脈拍変動の何れかである。以下の説明では、生体情報Ｘが脈拍数であり、生体情報Ｙが深部温度である場合を便宜的に例示する。
【００２３】
図２の測定処理部２１は、生体情報Ｘの測定値ｘと生体情報Ｙの測定値ｙとを測定する。測定処理部２１による測定値ｘおよび測定値ｙの測定は、所定長の単位期間毎に反復される。例えば、１秒毎に算定された数値を６０秒毎に平均して複数の地点にわたり移動平均を算定することで、生体情報Ｘの測定値ｘと生体情報Ｙの測定値ｙとが確定する。すなわち、単位期間毎に測定値ｘと測定値ｙとの組が算定される。なお、検出装置１３から出力される検出信号のノイズを除去するフィルターを測定処理部２１の前段に設置してもよい。
【００２４】
図２の解析処理部２２は、測定処理部２１から順次に供給される生体情報Ｘの測定値ｘと生体情報Ｙの測定値ｙとから生体の体内時計時刻を特定する。解析処理部２２による体内時計時刻の特定は、所定の周期で（例えば測定値ｘおよび測定値ｙの測定毎に）順次に実行される。
【００２５】
解析処理部２２による体内時計時刻の特定には基準リズム情報Ｄが利用される。基準リズム情報Ｄは、記憶装置１２に記憶されたデータであり、生体情報Ｘおよび生体情報Ｙの各々の周期的な変動と時刻との関係を規定する。基準リズム情報Ｄについて以下に詳述する。
【００２６】
生体情報における日単位の周期的な変動（以下「サーカディアンリズム（circadian rhythm）」という）は、以下の数式(1)および数式(2)で表現される。数式(1)は生体情報Ｘのサーカディアンリズムであり、数式(2)は生体情報Ｙのサーカディアンリズムである。
【数１】
2019055087.tif 000003
【００２７】
数式(1)および数式(2)の記号Ｍxおよび記号Ｍyはメサー（mesor）であり、周期的な変動における最大値と最小値との中央値に相当する。記号Ａxおよび記号Ａyは変動の振幅を意味し、記号Φxおよび記号Φyは変動の位相を意味する。記号Ｔは時刻（時間単位）である。
【００２８】
生体解析装置１００の利用者について測定された多数の測定値ｘの時系列が数式(1)で近似され、当該利用者の多数の測定値ｙの時系列が数式(2)で近似されるように、メサー（Ｍx,Ｍy）と振幅（Ａx,Ａy）と位相（Φx,Φy）とが事前に確定される。すなわち、生体情報の測定値を三角関数で近似（フィッティング）するコサイナー法により、生体情報Ｘのサーカディアンリズム（数式(1)）と生体情報Ｙのサーカディアンリズム（数式(2)）とが確定される。なお、以上の説明では生体解析装置１００の利用者の測定値（ｘ,ｙ）から生体情報Ｘおよび生体情報Ｙの各々のサーカディアンリズムを確定したが、例えば特定の母集団に属する生体について事前に観測された生体情報の統計値をサーカディアンリズムの確定に利用してもよい。また、サーカディアンリズムを確定する方法はコサイナー法に限定されない。例えばコンスタントルーチン法、自己相関法等の公知の生体リズム測定方法が任意に採用される。
【００２９】
図４に例示される通り、生体情報Ｘに対応するＸ軸（第１軸の例示）と生体情報Ｙに対応するＹ軸（第２軸の例示）とが設定された座標空間（以下「バイタル空間」という）を想定する。数式(1)で表現された生体情報Ｘのサーカディアンリズムと数式(2)で表現された生体情報Ｙのサーカディアンリズムとをバイタル空間に展開すると、図４に例示される通り、生体情報Ｘと生体情報Ｙとの関係を表す曲線（以下「基準リズム曲線」という）Ｃが規定される。
【００３０】
基準リズム曲線Ｃは、バイタル空間内の特定の地点Ｑc（座標（Ｘ0,Ｙ0））を中心点とする環状の曲線である。中心点ＱcのＸ座標は生体情報ＸのメサーＭxに相当し（Ｘ0＝Ｍx）、中心点ＱcのＹ座標は生体情報ＹのメサーＭyに相当する（Ｙ0＝Ｍy）。ただし、所定の時間にわたる生体情報Ｘの測定値ｘの代表値（例えば平均値や中央値）を中心点ＱcのＸ座標として利用し、生体情報Ｙの測定値ｙの代表値を中心点ＱcのＹ座標として利用してもよい。図４に例示される通り、基準リズム曲線Ｃ上の任意の１個の地点（Ｘ,Ｙ）は特定の時刻ｔに対応する。すなわち、基準リズム曲線Ｃを例えば反時計回りに辿ることは時刻ｔの進行に相当する。
【００３１】
記憶装置１２に記憶された基準リズム情報Ｄは、以上に説明した基準リズム曲線Ｃを規定する。具体的には、基準リズム情報Ｄは、図５に例示される通り、基準リズム曲線Ｃ上の生体情報Ｘ（脈拍数）および生体情報Ｙ（深部温度）の各数値を時刻ｔに対応させたデータテーブルである。なお、数式(1)および数式(2)における各変数（Ｍx,Ｍy,Ａx,Ａy,Φx,Φy）の数値が確定すれば基準リズム曲線Ｃは規定される。したがって、図６に例示される通り、メサー（Ｍx,Ｍy）と振幅（Ａx,Ａy）と位相（Φx,Φy）と頂点位相時刻とを生体情報Ｘおよび生体情報Ｙの各々について登録したデータテーブルを基準リズム情報Ｄとして記憶装置１２に記憶してもよい。以上の説明から理解される通り、第１実施形態の基準リズム情報Ｄは、生体情報Ｘおよび生体情報Ｙの各々の周期的な変動と時刻ｔとの関係を表すデータである。
【００３２】
図２の解析処理部２２は、以上に説明した基準リズム情報Ｄを利用して、測定処理部２１が測定した生体情報Ｘの測定値ｘと生体情報Ｙの測定値ｙとから生体の体内時計時刻を特定する。測定処理部２１が測定値ｘおよび測定値ｙを測定するたびに（すなわち単位期間毎に）、解析処理部２２は体内時計時刻を特定する。図４に例示される通り、生体情報Ｘの測定値ｘと生体情報Ｙの測定値ｙとは、バイタル空間内の任意の１個の地点（以下「測定点」という）Ｐに対応する。すなわち、測定値ｘは測定点ＰのＸ軸上の座標値であり、測定値ｙは測定点ＰのＹ軸上の座標値である。第１実施形態の解析処理部２２は、測定値ｘおよび測定値ｙに対応する測定点Ｐと、基準リズム曲線Ｃの中心点Ｑcとを通過する直線Ｌが、基準リズム曲線Ｃに交差する地点Ｑに対応する時刻ｔを、体内時計時刻として特定する。
【００３３】
直線Ｌは以下の数式(3)で表現される。
【数２】
2019055087.tif 000004
解析処理部２２は、数式(1)と数式(2)と数式(3)とを含む未知数３の連立方程式を解くことで体内時計時刻を算定する。なお、基準リズム曲線Ｃと直線Ｌとは２点で交差するため、連立方程式から２個の時刻ｔが特定される。解析処理部２２は、２個の時刻ｔのうち、基準リズム曲線Ｃ上で測定点Ｐに近い時点に対応する時刻ｔを、体内時計時刻として確定する。解析処理部２２が特定した体内時計時刻が、図１を参照して前述した通り、表示装置１４により表示される。
【００３４】
以上に説明した通り、第１実施形態では、生体情報Ｘおよび生体情報Ｙの各々の周期的な変動と時刻ｔとの関係を表す基準リズム情報Ｄを利用して、生体情報Ｘの測定値ｘと生体情報Ｙの測定値ｙとから体内時計時刻が特定される。すなわち、基準リズム情報Ｄが記憶装置１２に記憶された状態では、生体情報Ｘの測定値ｘと生体情報Ｙの測定値ｙとの１組さえ測定されれば体内時計時刻を特定でき、複数の測定値ｘの時系列や複数の測定値ｙの時系列は不要である。したがって、体内時計の時刻を実時間的に特定できるという利点がある。また、第１実施形態では、生体情報Ｘおよび生体情報Ｙの双方の測定値（ｘ,ｙ）が利用されるから、１種類の生体情報のみから体内時計時刻を特定する構成と比較して、体内時計の時刻を高精度に特定できるという利点もある。
【００３５】
第１実施形態では、バイタル空間内の測定点Ｐと基準リズム曲線Ｃの中心点Ｑcとを通過する直線Ｌが基準リズム曲線Ｃと交差する地点Ｑに対応する時刻ｔに応じて体内時計時刻が特定される。したがって、体内時計時刻を簡便な処理により特定できるという利点がある。
【００３６】
＜第２実施形態＞
本発明の第２実施形態について説明する。なお、以下に例示する各形態において作用または機能が第１実施形態と同様である要素については、第１実施形態の説明で使用した符号を流用して各々の詳細な説明を適宜に省略する。
【００３７】
第１実施形態では、生体情報Ｘおよび生体情報Ｙによりバイタル空間に基準リズム曲線Ｃを規定した。第２実施形態では、所定の範囲内の数値に正規化された生体情報Ｘnおよび生体情報Ｙnによりバイタル空間に基準リズム曲線Ｃnを規定する。
【００３８】
数式(1)で表現される生体情報Ｘは、以下の数式(4)で表現される正規化により、所定の範囲（−１≦Ｘn≦１）内の生体情報Ｘnに変換される。同様に、数式(2)で表現される生体情報Ｙは、以下の数式(5)で表現される正規化により、所定の範囲（−１≦Ｙn≦１）内の生体情報Ｙnに変換される。
【数３】
2019055087.tif 000005
【００３９】
図７に例示される通り、正規化後の生体情報Ｘnに対応するＸ軸と正規化後の生体情報Ｙnに対応するＹ軸とが規定されたバイタル空間を想定する。生体情報Ｘnのサーカディアンリズムと生体情報Ｙnのサーカディアンリズムとをバイタル空間に展開すると、図７に例示される通り、生体情報Ｘnと生体情報Ｙnとの関係を表す基準リズム曲線Ｃnが規定される。生体情報Ｘと生体情報Ｙとの位相差がπ/２に近似または一致する場合、基準リズム曲線Ｃnは単位円に相当する。
【００４０】
第２実施形態の基準リズム情報Ｄは、図８に例示される通り、基準リズム曲線Ｃn上の複数の地点Ｑnの各々について角度Θと時刻ｔとが登録されたデータテーブルである。図７に例示される通り、基準リズム曲線Ｃn上の任意の１個の地点Ｑnに対応する角度Θは、当該地点Ｑnに対応する位置ベクトルＶq（Ｖq＝（Ｘn,Ｙn））がＸ軸に対してなす角度（偏角）である。角度Θは、基準リズム曲線Ｃnの中心点Ｑc（すなわちバイタル空間の原点）と地点Ｑnとを通過する直線がＸ軸に対してなす角度とも換言される。以上の説明から理解される通り、基準リズム曲線Ｃn上の任意の地点Ｑnに対応する角度Θは、以下の数式(6)で表現される。
【数４】
2019055087.tif 000006
【００４１】
図９は、第２実施形態の解析処理部２２が体内時計時刻を特定する処理のフローチャートである。測定処理部２１が測定値ｘおよび測定値ｙを測定するたびに図９の処理が実行される。
【００４２】
図９の処理を開始すると、解析処理部２２は、測定処理部２１が測定した測定値ｘおよび測定値ｙを所定の範囲内の数値に正規化する（Ｓa1）。具体的には、解析処理部２２は、以下の数式(7)で表現される正規化により、生体情報Ｘの測定値ｘを所定の範囲（−１≦ｘn≦１）内の測定値ｘnに変換する。同様に、解析処理部２２は、以下の数式(8)で表現される正規化により、生体情報Ｙの測定値ｙを所定の範囲（−１≦ｙn≦１）内の測定値ｙnに変換する。添字ｎは正規化（normalization）を意味する。
【数５】
2019055087.tif 000007
【００４３】
解析処理部２２は、正規化後の測定値ｘnおよび測定値ｙnが示す測定点Ｐnの角度θを特定する（Ｓa2）。図７に例示される通り、角度θは、測定点Ｐnに対応する位置ベクトルＶp（Ｖp＝（ｘn,ｙn））がＸ軸に対してなす角度（偏角）である。角度θは、基準リズム曲線Ｃnの中心点Ｑc（すなわちバイタル空間の原点）と測定点Ｐnとを通過する直線がＸ軸に対してなす角度とも換言される。具体的には、第２実施形態の解析処理部２２は、正規化後の測定値ｘnおよび測定値ｙnについて以下の数式(9)の演算を実行することで角度θを算定する。
【数６】
2019055087.tif 000008
【００４４】
解析処理部２２は、測定点Ｐnに対応する位置ベクトルＶpの角度θから体内時計時刻を特定する（Ｓa3）。具体的には、解析処理部２２は、基準リズム情報Ｄに含まれる複数の角度Θのうち位置ベクトルＶpの角度θに近似または一致する角度Θを特定し、基準リズム情報Ｄにおいて当該角度Θに対応する時刻ｔを体内時計時刻として特定する。
【００４５】
第２実施形態においても第１実施形態と同様の効果が実現される。また、第２実施形態では、測定点Ｐnに対応する位置ベクトルＶpの角度θに応じて体内時計時刻が特定される。したがって、体内時計時刻を簡便な処理により特定できるという利点がある。
【００４６】
＜第３実施形態＞
図１０は、第３実施形態の解析処理部２２が体内時計時刻を特定する処理のフローチャートである。測定処理部２１が測定値ｘおよび測定値ｙを測定するたびに図１０の処理が実行される。図１０の処理を開始すると、解析処理部２２は、第２実施形態と同様に、測定処理部２１が測定した測定値ｘおよび測定値ｙを所定の範囲内に正規化する（Ｓb1）。すなわち、測定値ｘを正規化した測定値ｘnと測定値ｙを正規化した測定値ｙnとが算定される。
【００４７】
解析処理部２２は、図１１に例示される通り、位置ベクトルＶqと位置ベクトルＶpとがなす角度Δθの絶対値|Δθ|を算定する（Ｓb2）。位置ベクトルＶqは、基準リズム曲線Ｃにおいて任意の時刻（以下「実時刻」という）Ｔに対応する地点Ｑnの位置ベクトルである。具体的には、解析処理部２２は、位置ベクトルＶqと位置ベクトルＶpとの内積を含む以下の数式(10)の演算により角度Δθの絶対値|Δθ|を算定する。実時刻Ｔは、例えば現在時刻である。角度Δθは、実時刻Ｔを基準とした体内時計時刻に対応する。数式(10)の分子は、位置ベクトルＶqと位置ベクトルＶpとの内積である。
【数７】
2019055087.tif 000009
【００４８】
数式(10)により算定された角度Δθの絶対値|Δθ|は、基準リズム曲線Ｃnのうち実時刻Ｔに対応する地点Ｑnの位置ベクトルＶqと測定点Ｐnの位置ベクトルＶpとがなす角度である。ただし、測定点Ｐnの位置ベクトルＶpが位置ベクトルＶqに対して正側および負側の何れに位置するのかは、数式(10)の演算のみでは確定できない。位置ベクトルＶpが位置ベクトルＶqに対して正側および負側の何れに位置するのかを判別するために、解析処理部２２は、位置ベクトルＶpと位置ベクトルＶqとの外積crossを算定する（Ｓb3）。外積crossは以下の数式(11)で表現される。
【数８】
2019055087.tif 000010
【００４９】
解析処理部２２は、外積crossの符号により、角度Δθの符号を以下の数式(12)のように判別する（Ｓb4）。
【数９】
2019055087.tif 000011
【００５０】
解析処理部２２は、以上の処理で算定した角度Δθを含む以下の数式(13)の演算により、体内時計時刻τを算定する（Ｓb5）。
【数１０】
2019055087.tif 000012
以上の説明から理解される通り、第３実施形態においても第２実施形態と同様に、生体情報Ｘの測定値ｘnと生体情報Ｙの測定値ｙnとが示す測定点Ｐnの位置ベクトルＶpの角度θに応じて体内時計時刻τが特定される。
【００５１】
第３実施形態においても第１実施形態および第２実施形態と同様の効果が実現される。また、第３実施形態では、測定点Ｐnの位置ベクトルＶpと基準リズム曲線Ｃn上の地点Ｑnの位置ベクトルＶqとの内積および外積により体内時計時刻τが特定される。したがって、第１実施形態や第２実施形態と比較して、体内時計時刻を簡便な処理により特定できるという利点がある。
【００５２】
図１２は、第３実施形態により特定される体内時計時刻の誤差のグラフである。図１２には、脈拍数を利用したコサイナー法で算定された体内時計時刻の誤差（対比例１）と、深部温度を利用したコサイナー法で算定された体内時計時刻の誤差（対比例２）とが図示されている。第３実施形態によれば、対比例１および対比例２と比較して体内時計時刻の誤差が全体的に小さいことが確認できる。実際に誤差の標準偏差を算定すると、対比例１では６．１７［時間］であり、対比例２では６．１３［時間］であるのに対し、第３実施形態では４．４０［時間］である。すなわち、第３実施形態によれば、コサイナー法を利用した場合と比較して１．７時間ほど測定誤差を低減できる。
【００５３】
＜第４実施形態＞
第２実施形態および第３実施形態の例示から理解される通り、生体情報Ｘと生体情報Ｙとの位相差がπ/２に近似または一致する場合、基準リズム曲線Ｃnは円環状となる。他方、生体情報Ｘと生体情報Ｙとの位相差がπ/２から乖離する場合、基準リズム曲線Ｃnは楕円状または直線状となる。したがって、生体情報Ｘと生体情報Ｙとの位相差がπ/２から乖離する場合には、バイタル空間内での角度θを利用した体内時計時刻の測定精度が低下する可能性がある。第４実施形態は、生体情報Ｘと生体情報Ｙとの位相差がπ/２から乖離する場合でも体内時計時刻を高精度に特定するための形態である。
【００５４】
正規化後の生体情報Ｘnおよび生体情報Ｙnを加法定理で展開すると、以下の数式(14)および数式(15)が導出される。
【数１１】
2019055087.tif 000013
【００５５】
係数αおよび係数βを利用した生体情報Ｘnと生体情報Ｙnとの線形和（Ｘn＋αＹn,Ｘn＋βＹn）を、以下の数式(16)および数式(17)のように立式する。
【数１２】
2019055087.tif 000014
【００５６】
数式(16)の正弦成分と数式(17)の余弦成分とがゼロとなるように、以下の数式(18)および数式(19)の条件が成立すると仮定する。
【数１３】
2019055087.tif 000015
【００５７】
数式(18)および数式(19)から以下の数式(20)および数式(21)が導出される。
【数１４】
2019055087.tif 000016
【００５８】
数式(20)を数式(16)に適用することで以下の数式(22)が導出され、数式(21)を数式(17)に適用することで以下の数式(23)が導出される。
【数１５】
2019055087.tif 000017
【００５９】
数式(22)を正規化することで数式(24)の関数Ａnが導出され、数式(23)を正規化することで数式(25)の関数Ｂnが導出される。
【数１６】
2019055087.tif 000018
【００６０】
関数Ａnと関数Ｂnとの間の位相差はπ/２である。数式(24)の関数Ａnを示すＸ軸と数式(25)の関数Ｂnを示すＹ軸とが規定されたバイタル空間を想定すると、関数Ａnと関数Ｂnとの関係を表す円環状の基準リズム曲線Ｃが規定される。
【００６１】
測定処理部２１が算定した測定値ｘおよび測定値ｙについて、以上と同様の処理を実行することで、測定値ｘを補正した補正値ａnを表現する数式(26)と、測定値ｙを補正した補正値ｂnを表現する数式(27)とが導出される。
【数１７】
2019055087.tif 000019
【００６２】
第２実施形態または第３実施形態における基準リズム曲線Ｃを、数式(24)および数式(25)で規定される基準リズム曲線Ｃに置換し、第２実施形態または第３実施形態における測定値ｘnおよび測定値ｙnを補正値ａnおよび補正値ｂnに置換することで、体内時計時刻を特定することが可能である。以上の説明から理解される通り、第４実施形態によれば、生体情報Ｘと生体情報Ｙとの位相差がπ/２から乖離した関係にある場合でも、体内時計時刻を高精度に特定できるという利点がある。
【００６３】
＜補足＞
第１実施形態から第４実施形態では、相異なる２種類の生体情報（生体情報Ｘおよび生体情報Ｙ）を利用して体内時計時刻を特定する構成（以下「構成Ａ」という）を採用する。体内時計時刻を特定および表示する実際の生体解析装置（以下「実製品」という）が、構成Ａを採用しているか否かを判別するための方法を以下に説明する。
【００６４】
図１３に例示される通り、第１センサー５１と第２センサー５２と処理装置５３と表示装置５４とを具備する実製品５０を想定する。第１センサー５１および第２センサー５２は、相異なる物理量（例えば光量または温度）を検出する検出器である。処理装置５３は体内時計時刻を特定し、表示装置５４は、処理装置５３が特定した体内時計時刻を表示する。
【００６５】
第１に、第２センサー５２を無効状態に設定したうえで、体内時計時刻の測定を実製品５０に実行させる。無効状態とは、例えば第２センサー５２が目的の物理量を検出できない状態（例えば遮光または断熱された状態）、または、第２センサー５２が処理装置５３から電気的に絶縁された状態である。以上の状態では、第１センサー５１から供給される検出信号のみが、処理装置５３による体内時計時刻の測定に有効に利用される。表示装置５４は、処理装置５３が測定した体内時計時刻τ1を表示する。
【００６６】
第２に、第２センサー５２を有効状態に設定したうえで、体内時計時刻の測定を実製品５０に実行させる。以上の状態では、第１センサー５１および第２センサー５２の双方から供給される検出信号が、処理装置５３による体内時計時刻の測定に利用される。表示装置５４は、処理装置５３が測定した体内時計時刻τ2を表示する。
【００６７】
実製品５０が体内時計時刻の測定に第１センサー５１および第２センサー５２の双方の検出信号を利用している場合、第２センサー５２を無効状態として測定された体内時計時刻τ1と、第２センサー５２を有効状態として測定された体内時計時刻τ2とは相違する。したがって、体内時計時刻τ1と体内時計時刻τ2とが乖離する場合には、実製品５０が構成Ａを採用している可能性が充分に高い。具体的には、体内時計時刻τ1と体内時計時刻τ2との差異|τ1−τ2|が１０分以上である場合には、実製品５０が構成Ａを採用している可能性が高い。
【００６８】
次に、構成Ａを採用していることが確認された製品（以下「採用品」という）を想定する。前述の通り、第１センサー５１および第２センサー５２の双方が有効状態に設定された実製品５０により測定された体内時計時刻τ2が、採用品により測定された体内時計時刻τ3に近似する場合にも、実製品５０が構成Ａを採用している可能性が高い。具体的には、体内時計時刻τ1と体内時計時刻τ3との差異|τ1−τ3|が１０分以上である場合には、実製品５０が構成Ａを採用している可能性が高い。
【００６９】
＜変形例＞
以上に例示した各形態は多様に変形され得る。具体的な変形の態様を以下に例示する。以下の例示から任意に選択された２以上の態様を適宜に併合することも可能である。
【００７０】
（１）測定値ｘおよび測定値ｙを正規化する方法は前述の例示（数式(7)および数式(8)）に限定されない。例えば、以下の数式(28)の演算で測定値ｘを測定値ｘn'に正規化し、以下の数式(29)の演算で測定値ｙを測定値ｙn'に正規化してもよい。測定値ｘn'および測定値ｙn'を利用する構成では、位置ベクトルＶpが単位ベクトルとなる。
【数１８】
2019055087.tif 000020
【００７１】
（２）前述の各形態では日単位の体内時計時刻を例示したが、月周期（インフラディアンリズム）の体内時計において生体の状態が何れの時点に相当するかを示す指示値（以下「体内時計日」という）を、前述の各形態と同様の構成により特定することも可能である。インフラディアンリズムの典型例は、２８日から３０日程度の月経周期である。
【００７２】
具体的には、生体情報Ｘの１日内の代表値（例えば平均値）を測定値ｘとし、生体情報Ｙの１日内の代表値（例えば平均値）を測定値ｙとして、月周期にわたる測定値ｘおよび測定値ｙの時系列が生成される。複数の測定値ｘの時系列を三角関数により近似することで生体情報Ｘのインフラディアンリズムが確定され、複数の測定値ｙの時系列を三角関数により近似することで生体情報Ｙのインフラディアンリズムが確定される。生体情報Ｘのインフラディアンリズムと生体情報Ｙのインフラディアンリズムとをバイタル空間に展開することで、月周期の基準リズム曲線Ｃが規定される。月周期の基準リズム曲線Ｃは、生体情報Ｘおよび生体情報Ｙの各々の月周期の変動と日付との関係を表す。なお、月周期の基準リズム曲線Ｃの中心点Ｑcは、例えば、日周期を対象とした数式(1)のメサーＭxおよび数式(2)のメサーＭyで規定される。
【００７３】
解析処理部２２は、以上の方法で規定された基準リズム曲線Ｃを利用して、生体情報Ｘの測定値ｘ（１日内の代表値）と生体情報Ｙの測定値ｙ（１日内の代表値）とに応じた体内時計日を特定する。正規化後の測定値ｘnおよび測定値ｙnに対応する基準リズム曲線Ｃnを体内時計日の特定に利用してもよい。基準リズム曲線Ｃ（または基準リズム曲線Ｃn）を利用する方法は第１実施形態から第４実施形態と同様である。
【００７４】
図１４は、体内時計日を表示する表示装置１４の構成例である。図１４の表示装置１４は、図１と同様の指針（時針１４１，分針１４２，指針１４３）に加えて月齢表示盤１４４と指針１４５とを具備する。月齢表示盤１４４（ムーンフェイズ）は、月周期（２８日〜３０日）で１周する円板であり、月齢の変化（新月→満月→新月）を図形的に表示する。指針１４５は、月齢表示盤１４４に連動して回転し、体内時計日を表示する。月齢表示盤１４４および指針１４５を視認することで、利用者は月経やホルモンの状態を確認または予測することが可能である。なお、体内時計日は個人性の高い情報であるから、例えば利用者による操作に応じて表示／非表示を切替える構成が好適である。
【００７５】
（３）前述の各形態では、単体の機器として構成された生体解析装置１００を例示したが、以下の例示の通り、生体解析装置１００の複数の要素は相互に別体の装置として実現され得る。
【００７６】
前述の各形態では、検出装置１３を具備する生体解析装置１００を例示したが、図１５に例示される通り、検出装置１３を生体解析装置１００とは別体とした構成も想定される。検出装置１３は、例えば生体の特定の部位（例えば手首や上腕等）に装着される可搬型の光学センサーモジュールである。生体解析装置１００は、例えば携帯電話機またはスマートフォン等の情報端末で実現される。腕時計型の情報端末で生体解析装置１００を実現してもよい。検出装置１３が生成した検出信号が有線または無線により生体解析装置１００に送信される。生体解析装置１００の処理装置１１は、検出信号から体内時計時刻を算定して表示装置１４に表示させる。以上の説明から理解される通り、検出装置１３は生体解析装置１００から省略され得る。
【００７７】
前述の各形態では、表示装置１４を具備する生体解析装置１００を例示したが、図１６に例示される通り、表示装置１４を生体解析装置１００とは別体とした構成も想定される。生体解析装置１００の処理装置１１は、前述の各形態と同様の方法で特定した体内時計時刻の表示データを、有線または無線により表示装置１４に送信する。表示装置１４は、生体解析装置１００から受信した表示データを利用して体内時計時刻を表示する。なお、表示装置１４は、専用の表示機器であってもよいが、例えば、携帯電話機もしくはスマートフォン等の情報端末、または、被験者が携帯可能な腕時計型の情報端末に搭載されてもよい。以上の説明から理解される通り、表示装置１４は生体解析装置１００から省略され得る。
【００７８】
図１７に例示される通り、検出装置１３および表示装置１４を生体解析装置１００（処理装置１１）とは別体とした構成も想定される。なお、検出装置１３と生体解析装置１００とを別体とした構成において、測定処理部２１を検出装置１３に搭載することも可能である。測定処理部２１が測定した測定値ｘおよび測定値ｙが有線または無線により検出装置１３から生体解析装置１００に送信される。以上の説明から理解される通り、測定処理部２１は生体解析装置１００から省略され得る。
【００７９】
（４）前述の各形態では体内時計時刻を例示したが、生体解析装置１００が特定する情報は、体内時計時刻に限定されない。例えば、体内時計時刻と基準時刻（例えば現在時刻）との時間差を生体解析装置１００が特定してもよい。以上の説明から理解される通り、生体解析装置１００が特定する情報は、「体内時計の指示値」として包括的に表現される。「体内時計の指示値」の典型例は、第１実施形態から第４実施形態に例示した体内時計時刻であるが、月周期の体内時計における前述の体内時計日や、体内時計時刻と基準時刻との時間差も、「体内時計の指示値」の概念に包含される。
【００８０】
（５）前述の各形態では、体内時計時刻を表示装置１４に表示したが、体内時計時刻の表示は必須ではない。例えば、生体解析装置１００が特定した体内時計時刻を通信網から他の通信装置に送信してもよい。また、生体解析装置１００の記憶装置１２や生体解析装置１００に着脱可能な可搬型の記録媒体に体内時計時刻を格納してもよい。
【００８１】
（６）前述の各形態に係る生体解析装置１００は、前述の例示の通り、処理装置１１とプログラムとの協働により実現される。本発明の好適な態様に係るプログラムは、コンピューターが読取可能な記録媒体に格納された形態で提供されてコンピューターにインストールされ得る。また、配信サーバーが具備する記録媒体に格納されたプログラムを、通信網を介した配信の形態でコンピューターに提供することも可能である。記録媒体は、例えば非一過性（non-transitory）の記録媒体であり、ＣＤ-ＲＯＭ等の光学式記録媒体（光ディスク）が好例であるが、半導体記録媒体または磁気記録媒体等の公知の任意の形式の記録媒体を包含し得る。なお、非一過性の記録媒体とは、一過性の伝搬信号（transitory, propagating signal）を除く任意の記録媒体を含み、揮発性の記録媒体を除外するものではない。
【符号の説明】
【００８２】
１００…生体解析装置、１１…処理装置、１２…記憶装置、１３…検出装置、１３１…光センサー、１３２…熱センサー、１４…表示装置、２１…測定処理部、２２…解析処理部。
【図１】
2019055087.tif 000021
【図２】
2019055087.tif 000022
【図３】
2019055087.tif 000023
【図４】
2019055087.tif 000024
【図５】
2019055087.tif 000025
【図６】
2019055087.tif 000026
【図７】
2019055087.tif 000027
【図８】
2019055087.tif 000028
【図９】
2019055087.tif 000029
【図１０】
2019055087.tif 000030
【図１１】
2019055087.tif 000031
【図１２】
2019055087.tif 000032
【図１３】
2019055087.tif 000033
【図１４】
2019055087.tif 000034
【図１５】
2019055087.tif 000035
【図１６】
2019055087.tif 000036
【図１７】
2019055087.tif 000037