ホーム > 特許ランキング > ヘルスセンシング株式会社
知財求人 - 知財ポータルサイト「IP Force」
特許6358768睡眠状態測定装置、δ波推定方法、位相コヒーレンス算出装置及びストレス状態測定装置
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】6358768
(24)【登録日】2018年6月29日
(45)【発行日】2018年7月18日
(54)【発明の名称】睡眠状態測定装置、δ波推定方法、位相コヒーレンス算出装置及びストレス状態測定装置
(51)【国際特許分類】
A61B 5/16 20060101AFI20180709BHJP
A61B 5/113 20060101ALI20180709BHJP
【ＦＩ】
A61B5/16 130
A61B5/16ZDM
A61B5/113
【請求項の数】20
【全頁数】42
(21)【出願番号】特願2017-545768(P2017-545768)
(86)(22)【出願日】2017年2月15日
(86)【国際出願番号】JP2017005556
(87)【国際公開番号】WO2017141976
(87)【国際公開日】20170824
【審査請求日】2017年9月4日
(31)【優先権主張番号】特願2016-26008(P2016-26008)
(32)【優先日】2016年2月15日
(33)【優先権主張国】JP
【早期審査対象出願】
(73)【特許権者】
【識別番号】514118321
【氏名又は名称】ヘルスセンシング株式会社
(74)【代理人】
【識別番号】100123858
【弁理士】
【氏名又は名称】磯田 志郎
(72)【発明者】
【氏名】鐘ヶ江 正巳
(72)【発明者】
【氏名】藤田 努
(72)【発明者】
【氏名】新関 久一
【審査官】
九鬼 一慶
(56)【参考文献】
【文献】
特開平０６−２６１８６５（ＪＰ，Ａ）
【文献】
特表２０１１−５２１６７４（ＪＰ，Ａ）
【文献】
特表２０１１−５１７９８２（ＪＰ，Ａ）
【文献】
山形大学工学部，ヘルスセンシング（株）と山形大学の共同研究がＪＳＴの研究成果最適展開支援プログラム（Ａ−ＳＴＥＰ）に採択，プレス発表資料［オンライン］，日本，２０１５年 ２月１７日，１頁，［検索日 2018.2.1］，インターネット：＜http://www2.yz.yamagata-u.ac.jp/kaiken/150217press.pdf＞
(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】
非特許文献２には、心拍間隔の変動と呼吸パターンの瞬時位相差の位相コヒーレンスが心理ストレスと相関しており、位相コヒーレンスを心理ストレスの指標とが評価できる可能性があると開示されている。非特許文献２では、胸部に貼付した電極により心電図を取得し、胸部及び腹部に巻いて装着した呼吸努力ベルトで呼吸パターンを実測し、心電図から心拍間隔の変動（呼吸性洞性不整脈）を算出し、さらに心拍間隔の変動の瞬時位相を算出し、また、呼吸パターンの瞬時位相を算出し、これらから時刻ｔにおける心拍と呼吸リズムの瞬時位相差Ψ（ｔ）jを求め、得られた瞬時位相差を使用して、式（１）により位相コヒーレンスを算出することが記載されている。式（１）のＮはサンプリングされたデータ数であり、Ｎ個平均して求める。【0006】
【数1】
【0007】
位相コヒーレンスは、０〜１の値をとり、心拍間隔の変動と呼吸パターンの瞬時位相差が一定の関係に近いほど１に近づき、瞬時位相差がランダムになるほど０に近づく。安静リラックス時には位相コヒーレンスが１に近く，精神ストレスがかかると位相コヒーレンスが低下するので，心理ストレスを正規化した指標で推定することが可能となる。【0008】
ところで、ヒトの睡眠はノンレム睡眠とレム睡眠に大別され、睡眠中に約９０分周期で繰り返すと言われている。客観的に眠りの質や量（睡眠周期や睡眠深度）を評価する為には、睡眠ポリグラフィが必要であり、医療機関での計測や医療機関から機器を借りてきて計測が行われるが、脳波や心電図、筋電図、眼電図などの複数のセンサが必要であり、取り扱いが煩雑である。また、複数のセンサを身体の各所に固定する必要があり、センサによる拘束で動きが制限されていた。睡眠中の無意識の動作によってセンサが外れる可能性もあり、センサが外れると信号が検出できず、状態を把握できない点でも問題があった。さらに、装置が大型であるため、個人で評価できるような指標ではなく、睡眠の質を個人が把握できる状況には至っていない。このため、個人が簡便にストレス状態や睡眠の質を把握できる生体計測手法とその装置が求められている。【0009】
また、高齢者の孤独死や介護施設で不幸にして亡くなった場合、事故死なのか病死なのか判断がつかない場合が多い。このような事態を防ぐために、日常の身体活動を拘束せずに無侵襲で容易に計測できる生体情報計測機器の登場が待ち望まれている。特許文献２には、ベッドパッド又はマットレスの上部や下部に、人の身体から発生する振動を検出する圧電センサを設置し、非拘束タイプの圧電センサを用いて身体振動信号を検出することが開示されている。さらに、圧電センサによって検出した身体振動信号を差動信号増幅アンプによって増幅した後、対象とする人に起因する身体振動である心臓の拍動（心拍数換算：３０〜２４０回／分、周波数帯域換算：０．５〜４Ｈｚ）による振動、肺呼吸活動（呼吸数換算：６０回／分以下、周波数帯域換算：１Ｈｚ以下）による振動、鼾に起因する鼾振動を分離フィルタ機能を用いて抽出することが開示されている。特許文献２では、心臓の拍動に起因する拍動振動、肺呼吸活動に起因する肺呼吸振動または鼾振動が有りの状態が所定の存在継続時間以上を超えることで、人が所定場所に存在していること、及び拍動振動、肺呼吸振動または鼾振動が無い状態が所定の不在継続時間以上を超えることで、人が所定場所に不在であると判定している。【0010】
特許文献２の人存在不在検出装置では、１個の圧電センサを用いて身体振動信号を取得し、身体振動信号を周波数フィルタにより分離して呼吸、心拍、鼾、体動の４つの信号を得ていた。周波数フィルタは、コンデンサや抵抗及びオペアンプ等で構成されたローパスフィルタ（ＬＰＦ）やハイパスフィルタ（ＨＰＦ）のアナログフィルタ、又は身体振動信号をＡ／Ｄコンバータでデジタル信号に変換し数値化されたデータをもとにＣＰＵ（中央処理装置）の演算処理にてフィルタリングを行なうデジタルフィルタの何れか一方又は双方で構成することが開示されている。【先行技術文献】
【特許文献】
【0011】
【特許文献1】特開２００６−３２５７６６号公報
【特許文献2】特開２０１３−２１０３６７号公報
【非特許文献】
【0012】
【非特許文献1】Gary G Berntson, J. Thomas Bigger Jr. et al., “Heart rate variability: Origin, methods, and interpretive caveats”. Psychophysiology (USA) 34巻, p.623-648, 1997年
【非特許文献2】Niizeki K and Saitoh T., “Incoherent oscillations of respiratory sinus arrhythmia during acute mental stress in humans”. American Journal of Physiology, Heart and Circulatory Physiology (USA) 302巻, p.359-367, 2012年
【発明の概要】
【発明が解決しようとする課題】
【0013】
ストレス評価とも関連するが、近年、睡眠障害が多発しており、自己管理の一環として簡易に睡眠状態を測定、観察又は評価できる手段が望まれている。睡眠の質や量（睡眠周期や睡眠深度）を客観的に評価するためには、医療機関における計測や大型の専用装置が必要であり、一般的な通常生活において測定することは困難である。また、従来の装置は、脳波や心電図、筋電図、眼電図などを取得するための複数のセンサを必要とし、複数のセンサを身体の各所に固定する必要がある。また、センサによる拘束で動きが制限されるため、できるだけセンサの数を減らしたいとの要望もあった。睡眠中の無意識の動作によってセンサが外れる可能性もあることから、取り付けが容易なセンサで測定できることも望まれていた。さらに、複数のセンサを身体の各所に固定する煩雑さを緩和することも求められていた。このように、個人が簡便にストレス状態や睡眠状態を測定できる装置が求められている。【0014】
心理ストレスの指標となる情報を取得する装置は、ウェアラブルであり、できるだけ身体の動きに対して邪魔にならないような構成であることが求められる。この点で、特許文献１に記載の指輪型の生体信号測定装置は、手指に装着できるので好ましい。しかしながら、かかる指輪型の生体信号測定装置は、ＬＥＤ発光部を有する光学センサによって、血中のヘモグロビンによる光の吸収度を測定し、心拍に同期する光の吸収度の変動を測定することによって心拍を測定するものであるため、ＬＥＤ発光部などを動作させるための消費電力が大きいという問題がある。さらに、かかる光学センサを用いた生体信号測定装置は、ＬＥＤ発光部からの光が血管に照射される必要があり、血管から光がずれると生体信号が測定できなくなるという問題があり、生体信号測定装置の装着位置が極めて限定され、装着後のずれも防止する必要があった。また、特許文献１に記載の指輪型の生体信号測定装置によれば、脈波から脈拍を測定できるとされているが、呼吸に関する情報については測定できなかった。【0015】
別の課題として、非特許文献１に記載されたストレス評価法は、上記のとおり、ストレス負荷がない状態でも指標が変化すること、年齢の依存性やフィットネスの影響など個人差も大きく定量的に評価できないという問題があった。また、非特許文献１に記載されたストレス評価法では、心電図から抽出された心拍リズムゆらぎをストレス指標としているが、かかる指標は、心拍リズムゆらぎ呼吸周波数や一回換気量の影響を受け、ストレス負荷がない状態でも変化するため、汎化的な評価法になっていない。また、年齢の依存性やフィットネスの影響など個人差も大きく定量的に評価できないという問題があった。【0016】
非特許文献２に記載された方法では、胸部に貼付した電極により心電図を取得し、胸部及び腹部に巻いて装着した呼吸努力ベルトで呼吸パターンを実測し、心拍と呼吸パターンから位相コヒーレンスを算出し、位相コヒーレンスに基づいてストレス評価が可能であるとされている。しかしながら、非特許文献２に記載された方法では、心拍を検出するための心電図用のセンサと呼吸パターンを検出するためのセンサが必要であった。特に、呼吸パターンを検出するためのセンサとして、口と鼻を覆うマスクを装着して呼吸量を測定するセンサの場合は、日常生活、運動中又は睡眠中において測定するのが困難である上に、マスク装着によってストレスが生じる虞があった。よって、より負担がなく、ストレスフリーな状態でストレスを評価可能な手法の開発が望まれていた。【0017】
また、別の課題として、特許文献２に記載されているように、１つの圧電センサから取得した身体振動信号から周波数フィルタを用いて心拍情報や呼吸情報を取得する場合、抽出される呼吸情報や心拍情報として、できるだけ精度の高い情報を得ることが望まれていた。特に、上記のようなストレス評価にも使用可能な精度で心拍情報又は呼吸情報を複数の情報が混在した信号から間接的に心拍間隔の変動や呼吸パターンを取得する方法が求められている。【0018】
本発明は、上記の課題の少なくとも一部を解決することを目的とするものである。【課題を解決するための手段】
【0019】
前記課題を解決するために、本発明の睡眠状態測定装置は、動物の睡眠中に取得した心拍間隔の変動の瞬時位相と、前記心拍間隔の変動と同一時系列における前記動物の呼吸パターンの瞬時位相との間の瞬時位相差に基づいて位相コヒーレンスを算出する位相コヒーレンス算出手段を含む。さらに、睡眠状態測定装置は、心電図を計測するセンサと、前記センサで計測した動物の心電図から呼吸パターンに関する信号を抽出する呼吸波形抽出手段とを有していてもよい。また、睡眠状態測定装置は、振動を計測するセンサと、前記センサで計測した信号から心拍間隔の変動を算出する心拍間隔算出手段と、前記センサで計測した信号から呼吸パターンに関する信号を抽出する呼吸波形抽出手段とを含んでいてもよい。上記睡眠状態測定装置は、前記センサのサンプリング周波数が１００Ｈｚ以上であることが好ましい。さらに、前記センサは、動物の体に装着可能なウェアラブルセンサであることが好ましく、前記ウェアラブルセンサは、ヒトの四肢又は頭部に装着する装着部に実装されていることがより好ましい。さらに、上記睡眠状態測定装置は、前記位相コヒーレンス算出手段で算出した位相コヒーレンスに基づいて、睡眠状態等を判定する判定機能を有していてもよい。前記判定機能は、さらに前記呼吸パターンから睡眠中の呼吸状態を判定してもよい。【0020】
また、本発明の睡眠状態を測定する方法は、動物の睡眠中に取得した同一時系列における心拍間隔の変動と呼吸周波数との間の瞬時位相差の位相コヒーレンスに基づいて睡眠状態を測定する。さらに、上記睡眠状態を測定する方法において、前記動物の心電図を計測し、計測した心電図から呼吸パターンに関する信号を取得してもよい。また、上記睡眠状態を測定する方法において、前記動物の心弾動図波形又は生体振動信号を計測し、計測した心弾動図波形又は生体振動信号から心拍間隔の変動と呼吸パターンに関する信号を取得してもよく、さらに、前記動物の心弾動図波形又は生体振動信号から心弾動の伝達特性を推定し、前記伝達特性の逆伝達関数を前記動物の心弾動図波形又は生体振動信号に用いて心拍波形を取得することが好ましい。【0021】
また、本発明の位相コヒーレンス算出装置は、少なくとも動物の心拍に関する情報及び呼吸に関する情報の両方を含む生体情報を取得する生体情報取得手段と、前記生体情報から呼吸パターンを抽出する呼吸波形抽出手段と、前記生体情報から心拍間隔の変動を算出する心拍間隔算出手段と、前記呼吸パターンと前記心拍間隔の変動との間の瞬時位相差の位相コヒーレンスを算出する位相コヒーレンス算出手段とを含む。上記位相コヒーレンス算出装置において、前記生体情報は心電図であってもよいし、前記生体情報は心弾動図波形又は生体振動信号であってもよい。上記位相コヒーレンス算出装置において、前記生体情報は１００Ｈｚ以上のサンプリング周波数で取得されることが好ましい。さらに、前記生体情報取得手段は、動物の体に装着可能なウェアラブルセンサを含むことが好ましく、前記ウェアラブルセンサは、ヒトの四肢又は頭部に装着する装着部に実装されていることがより好ましい。【0022】
また、本発明のストレス測定装置は、上記位相コヒーレンス算出装置を備えたものである。また、本発明の睡眠状態測定装置は、上記位相コヒーレンス算出装置を備えたものである。【0023】
また、本発明の心拍波形抽出方法は、動物の心弾動図波形又は生体振動信号から心拍波形を抽出する方法であって、前記心弾動図波形又は生体振動信号から心弾動の伝達特性を推定し、前記伝達特性の逆伝達関数を前記心弾動図波形又は生体振動信号に用いて心拍波形を取得する。上記心拍波形抽出方法において、前記心弾動図波形又は生体振動信号を一回の心拍が含まれるように波形断片を切り出し、複数の前記波形断片を重ね合わせて平均化して心弾動の伝達特性を推定することが好ましい。【0024】
また、本発明の心拍間隔算出方法は、動物の心弾動図波形又は生体振動信号から心拍間隔を算出する方法であって、前記心弾動図波形又は生体振動信号に由来する信号に対し、心拍の周波数よりも高い周波数を下限周波数とするハイパスフィルタを通過させ、前記ハイパスフィルタを通過させた後の信号について絶対値をとる処理を行うことを特徴とする。【0025】
また、本発明の信号処理方法は、動物の心拍に関する情報を含んだ生体情報又は呼吸に関する情報を含んだ生体情報の信号処理方法であって、前記生体情報に由来する信号のパワースペクトルから上限周波数又は下限周波数を求め、前記上限周波数又は下限周波数を遮断周波数とするフィルタを通過させる処理を含むことを特徴とする。【発明の効果】
【0026】
本発明の睡眠状態測定装置及び方法は、心拍間隔の変動の瞬時位相と呼吸パターンの瞬時位相との瞬時位相差の位相コヒーレンスが、δ波の振幅と相関関係があることを利用して睡眠状態を測定することができる。位相コヒーレンスは、呼吸周波数の影響を受けにくく、より正確に睡眠状態を測定することができる。また、一つのセンサでも睡眠状態測定装置を実現することが可能であり、使用者に負担が少なくストレスを感じにくい装置を提供できる。特に振動を計測するセンサの場合は動物を拘束せずに計測することが可能であり、より簡便な装置を提供できる。位相コヒーレンスは、リアルタイムに測定することができ、ノンレム睡眠時間の発生時点やノンレム―レム睡眠のリズム周期を測定することができ、睡眠状態を的確に把握することができる。また、同時に呼吸パターンを測定又は抽出しているので、睡眠時の無呼吸状態を検出することも可能である。【0027】
本発明の位相コヒーレンス算出装置は、少なくとも動物の心拍に関する情報及び呼吸に関する情報の両方を含む生体情報から、心拍間隔の変動と呼吸パターンを取得して、位相コヒーレンスを算出するので、一つのセンサで実現され、使用者に負担が少なくストレスを感じにくい装置を提供できる。特に振動を計測するセンサの場合は動物を拘束せずに計測することが可能であり、より簡便な装置を提供できる。位相コヒーレンスを測定することで、呼吸周波数に依存せずにリアルタイムで被測定者の心理ストレス状態や睡眠状態を確認することが可能となる。【0028】
また、上記装置において、動物の体に装着可能なウェアラブルセンサ、特にヒトの四肢又は頭部、例えば、ヒトの手指に装着する装着部に実装されるウェアラブルセンサとしたので、被験者にとって負担が少なくストレスを感じにくく、ヒトの動きを拘束することがなく、より簡便な装置を実現できる。また、心拍に関する情報及び呼吸に関する情報の両方を含む生体振動信号を測定することが可能であり、心拍及び呼吸に関する生体情報をリアルタイムで算出することができる。特に、心拍間隔の変動の瞬時位相と呼吸パターンの瞬時位相との瞬時位相差の位相コヒーレンスもリアルタイムで算出することができ、呼吸周波数に依存せずにリアルタイムで被験者の心理ストレス状態や睡眠状態を確認することが可能となる。【図面の簡単な説明】
【0029】
【図1】睡眠中の脳波におけるδ波の振幅及び睡眠中の位相コヒーレンスの時間的変化を示す図。
【図2】睡眠状態測定装置の概略ブロック図。
【図3】睡眠状態測定装置の一実施形態。
【図4】（Ａ）は、心電図波形（上）と抽出した呼吸パターン（下実線）及び呼吸パターンの実測値（下点線）であり、（Ｂ）は、心拍間隔の変動を示すグラフであり、（Ｃ）は、実線が心拍間隔の変動の瞬時位相（上実線）及び呼吸パターンの瞬時位相（上点線）と算出した位相コヒーレンスλ（下）である。
【図5】呼吸数を１分間に８、１０、１２、１５、１８、２０、２４回で変化させた時の心電図から抽出した推定呼吸周波数と、実測呼吸周波数の相関関係を示す図。
【図6】呼吸数を１分間に８、１０、１２、１５、１８、２０、２４回で変化させた時の推定呼吸周波数を用いて算出した位相コヒーレンスλecgと実測呼吸周波数を用いて算出した位相コヒーレンスλの相関関係を示す図。
【図7】位相コヒーレンス算出装置の概略ブロック図。
【図8】（Ａ）は、心弾動図波形（上）と心弾動図波形から抽出された模擬心電図波形（下）であり、（Ｂ）は、心弾動図波形から求めた平均の伝達特性。
【図9】心電図波形を取得する際のサンプリング周波数として、１ｋＨｚ、５００Ｈｚ、２００Ｈｚ、１００Ｈｚ、５０Ｈｚと変化させた場合の位相コヒーレンスを比較した図。
【図10】（Ａ）は、安静時における心拍（実線）と呼吸（点線）の瞬時位相であり、（Ｂ）は、安静時の心拍と呼吸の瞬時位相のリサージュ図である。（Ｃ）は、暗算課題時（ストレス状態）における心拍（実線）と呼吸（点線）の瞬時位相であり、（Ｄ）は、暗算課題時の心拍と呼吸の瞬時位相のリサージュ図であり、（Ｅ）及び（Ｆ）は、安静状態と暗算課題時の呼吸性不整脈（ＲＳＡ）及び位相コヒーレンスの変化を示す図。
【図11】自発的に呼吸周波数を１５回／分から２５回／分に変化させた場合の心拍間隔（ＲＲＩ）、呼吸性不整脈の振幅（ＡRSA）、位相コヒーレンス（λ）及び呼吸周波数（ｆR）を示す図。
【図12】本発明の生体振動信号測定装置の概略ブロック図。
【図13】シート型圧電センサの構造を示す図。
【図14】実施例で検出、算出した信号波形を示す図。
【図15】心電図及び生体振動信号由来の信号のパワースペクトル密度（ＰＳＤ）を示す図。
【図16】第１の実施形態の指輪型の生体振動信号測定装置の構造を示す概略図。
【図17】第２の実施形態の腕輪型の生体振動信号測定装置の構造を示す概略図。
【図18】振動センサ（圧電センサ）の一例を示す図。
【図19】（Ａ）は、被験者の手指で検出した生体振動信号の原信号、（Ｂ）は、その処理後信号の波形である。
【図20】（Ａ）は、被験者の手首で検出した生体振動信号の原信号、（Ｂ）は、その処理後信号の波形である。
【図21】（Ａ）（Ｂ）（Ｃ）は、各々、生体振動信号から得られた心拍動パルス、心電図から得られた心拍動パルス、実測した呼吸パターンの例である。
【図22】仰臥位姿勢の被験者の人差し指付根で検出した生体振動信号から算出した各種信号を示す図。
【図23】座位姿勢の被験者の人差し指付根で検出した生体振動信号から算出した各種信号を示す図。
【図24】立位姿勢の被験者の人差し指付根で検出した生体振動信号から算出した各種信号を示す図。
【図25】（Ａ）は、生体振動信号から推定した心拍間隔ＲＲＩ、（Ｂ）は、心拍間隔のパワースペクトル波形である。
【図26】生体振動信号から求めた位相コヒーレンスλと心拍間隔のパワースペクトル成分との関係を示す図。
【図27】座位姿勢の被験者の上腕部で検出した生体振動信号から算出した各種信号を示す図。
【図28】立位姿勢の被験者のこめかみで検出した生体振動信号から算出した各種信号を示す図。
【発明を実施するための形態】
【0030】
［睡眠状態測定装置及び方法］睡眠の状態は、脳波と眼球運動のパターンによって分類されている。脳波は脳の神経活動に伴う電位差によって発生する電流の時間的変化を記録した脳電図（Electroencephalogram、「ＥＥＧ」）で計測される。覚醒状態において、安静時にはα波（周波数：８〜１３Ｈｚ、振幅：約５０μＶ）が主成分であるが、精神活動中はα波が抑制され、振幅の小さいβ波（周波数：１４〜Ｈｚ、振幅：多くは約３０μＶ以下）が現れる。浅い睡眠の段階ではα波が徐々に減少し、θ波（周波数：４〜８Ｈｚ）が現れ、深い睡眠中はδ波（周波数：１〜４Ｈｚ）が現れる。α波よりも周波数の速いβ波を速波と呼ぶこともある。また、α波よりも周波数の遅い脳波を徐波と呼ぶこともあり、θ波やδ波は徐波の一種である。
【0031】
睡眠は、大きくはレム睡眠とノンレム睡眠に分類されるが、ノンレム睡眠は、さらにステージＩからＩＶまでの４段階に分類される。ステージＩ（まどろみ期、入眠期）は、うとうとした状態であり、α波のリズムが失われ、徐々に平坦化する。ステージＩＩは、ステージＩよりは眠りが深い状態であり、睡眠紡錘波（ｓｐｉｎｄｌｅ）とＫ複合波が出現する。ステージＩＩＩは、かなり深い睡眠であり、δ波が２０％以上５０％未満を占める。ステージＩＶは、δ波が５０％以上を占める段階であり、最も深い睡眠状態である。ノンレム睡眠のステージＩＩＩ及びＩＶは、脳波が徐波を示すことから徐波睡眠とも呼ばれる。一方、レム睡眠は、覚醒時に似た低振幅速波パターンの脳波を示すのに、覚醒させるために徐波睡眠よりも強い刺激を必要とする深い睡眠状態であり、急速眼球運動を伴う。睡眠状態は、通常、約９０分の周期（睡眠周期）で変動し、ステージＩから徐々に睡眠が深くなり、その後、一旦睡眠が浅くなり、レム睡眠に移行する。ただし、睡眠時間にも影響し、睡眠の前半ではノンレム睡眠が優勢であり、後半はレム睡眠が優勢となることが多い。上記のように、脳波における周波数１〜４Ｈｚのδ波は、徐波睡眠において観察されることから、δ波を睡眠の深さの指標とすることができる。【0032】
本発明者らは、睡眠中の心拍間隔の変動と呼吸パターンの瞬時位相差の位相コヒーレンス（以下、「位相コヒーレンス」という）が、睡眠中の脳波におけるδ波と相関していることを発見し、位相コヒーレンスを計測することにより、睡眠状態を測定できることを発明したのである。図１は、睡眠中の検体の脳波におけるδ波の振幅（δ）と位相コヒーレンス（λ）の時間的変化を示すものである。図１から、δ波の振幅が大きくなると位相コヒーレンス（λ）が１に近づき、δ波の振幅が小さくなると位相コヒーレンス（λ）が０に近づく傾向があり、両者の相関関係が確認された。被験者１１名による計測では、λとδ波との相関係数の平均は０．５３であり、その標準偏差は０．１０であった。また，λの変化はδ波の変化に平均で１１．６分先行していた。【0033】
位相コヒーレンスは、２つの信号間の位相差のばらつきの程度を示すものであり、本発明においては、呼吸由来の心拍間隔の変動の位相差と呼吸パターンの位相差とのばらつきの程度を示す。心拍間隔は、呼吸性不整脈（ＲＳＡ）により、息を吸ったとき（吸気時）に短くなり、息を吐いたとき（呼気時）に長くなることから、心拍間隔は呼吸の影響により呼吸パターンと類似の周期で変動する。そして、上記の発見によれば、睡眠の深度が深い場合（δ波が優勢の場合）は、心拍間隔の変動と呼吸パターンの位相差のばらつきは少なく（位相コヒーレンスが１に近い）、睡眠の深度が浅くなるとばらつきが大きくなる（位相コヒーレンスが０に近い）。【0034】
位相コヒーレンスは、同一時系列における心拍間隔の変動の瞬時位相と呼吸パターンの瞬時位相との瞬時位相差を各心拍間隔の情報及び呼吸パターンから算出することが可能である。各心拍間隔の情報は、心拍の情報を含んだ生体情報から取得することができる。心拍の情報を含んだ生体情報とは、例えば、心拍に伴う活動電位の時間的変化を計測した心電図、心拍に伴う振動の時間的変化を計測した心弾動図波形、動物の振動（心弾動を含む）の時間的変化を計測した生体振動信号等が利用可能である。また、呼吸パターンは、呼吸パターンの情報を含んだ生体情報から取得することもできる。呼吸パターンの情報を含んだ生体情報とは、例えば、呼吸による空気の流れを計測した実測値、呼吸に伴う胸郭インピーダンス変化を計測した実測値、呼吸による温度変化を計測した実測値、呼吸運動に伴う腹部の動きを計測した実測値、心電図、心弾動図波形、生体振動信号等が利用可能である。心電図は、心拍に伴う心筋活動電位の時間的変化を体表面から計測するものであるが、呼吸に伴う体表面電位の時間的変化も含むものであり、心拍の情報及び呼吸パターンの情報の両方を含む生体情報である。また、心弾動図波形は、心拍に伴う振動の時間的変化を計測するものであるが、呼吸に伴う振動の時間的変化も含むものであり、心拍の情報及び呼吸パターンの情報の両方を含む生体情報である。また、生体振動信号は、動物の振動を計測するものであり、かかる振動には、心拍に伴う振動、呼吸に伴う振動も含まれており、心拍の情報及び呼吸パターンの情報の両方を含む生体情報である。心電図、心弾動図波形又は生体振動信号の一つ又は複数を組み合わせて、位相コヒーレンスの算出に必要なレベルで心拍の情報又は呼吸パターンの情報を抽出し、抽出した心拍の情報又は呼吸パターンの情報に基づいて位相コヒーレンスの算出に利用することができる。【0035】
位相コヒーレンスは、呼吸に伴う心拍間隔の変動の瞬時位相ψh（ｔ）と呼吸パターンの瞬時位相ψr（ｔ）とを算出し、これらの差（瞬時位相差）を算出し、算出された瞬時位相差を使用して算出される。心拍間隔の変動の瞬時位相ψh（ｔ）は、心拍間隔のデータから呼吸に伴う心拍間隔変動の時間的変化（S（ｔ））を算出し、心拍間隔変動の時間的変化（S（ｔ））を下記式（２）で表されるヒルベルト変換により解析信号にすることによって算出できる。なお、式（２）及び式（３）のＨ［…］はヒルベルト変数であり、Ｐ．Ｖ．はコーシーの主値を意味する。【0036】
【数2】
【0037】
また、呼吸パターンの瞬時位相ψr（ｔ）は、呼吸パターンの情報から、呼吸パターン時間的変化（Ｒ（ｔ））を算出し、呼吸パターンの時間的変化（Ｒ（ｔ））を下記式（３）で表されるヒルベルト変換することによって算出できる。【0038】
【数3】
【0039】
瞬時位相差Ψ（ｔ）は、式（２）及び式（３）によって算出された心拍間隔の変動の瞬時位相ψh（ｔ）及び呼吸パターンの瞬時位相ψr（ｔ）を用いて、次の式（４）で算出できる。Ψ（ｔ）＝ψh（ｔ）−ψr（ｔ）+２ｎπ （４）
ここで、ｎは−π≦Ψ≦πとなる適当な整数である。
【0040】
そして、これらから時刻ｔk における位相コヒーレンスを下記式（１）から算出することができる。式（１）のＮはサンプリングされたデータ数であり、Ｎ個平均して求める。【0041】
【数4】
【0042】
位相コヒーレンスを利用した睡眠状態測定装置１は、図２に示すように、少なくとも情報取得部２及び情報処理部３を備えている。さらに睡眠状態測定装置１は、操作部４、出力部５及び記憶部６を備えていてもよい。【0043】
情報取得部２は、位相コヒーレンスの算出に必要な情報を取得するものであり、動物を計測するためのセンサ及びセンサの情報を有線又は無線で入力する入力部を含む構成であってもよいし、すでに計測済みの情報が記録された他の記録媒体からの情報を有線又は無線で入力可能な入力部を含む構成であってもよい。すなわち、情報取得部２は、少なくとも情報を入力する入力部を備えており、場合によっては入力部と有線又は無線で接続された生体情報を計測するためのセンサを備えていてもよい。センサで生体情報を計測する場合、サンプリング周波数は１００Ｈｚ以上であることが好ましい。【0044】
情報処理部３は、入力された情報を処理するものであり、例えば、コンピュータのＣＰＵ（中央処理装置）の演算処理機能を利用することができる。また、情報処理の中には、デジタル回路ではなくアナログ回路で実現することも可能である。例えば、情報処理として周波数フィルタを行う場合は、コンデンサや抵抗及びオペアンプ等で構成されたローパスフィルタ（ＬＰＦ）やハイパスフィルタ（ＨＰＦ）のアナログフィルタで実現してもよいし、ＣＰＵの演算処理機能によってフィルタリングを行なうデジタルフィルタで実現してもよい。情報処理部３は、情報処理の種類に応じて、デジタル回路とアナログ回路の両方を含んでいてもよいし、入力される情報がアナログであれば、アナログ−デジタル変換回路によってデジタル信号に変換してもよい。情報処理部３は、入力される情報によって必要となる機能又は処理が異なるが、少なくとも心拍間隔の変動の瞬時位相と呼吸パターンの瞬時位相との瞬時位相差を用いて位相コヒーレンスを算出する位相コヒーレンス算出機能を有している。なお、位相コヒーレンス算出機能を有する装置を「位相コヒーレンス算出装置」と呼び、本発明の睡眠状態測定装置１も位相コヒーレンス算出装置の一種である。【0045】
位相コヒーレンスを算出するためには、例えば、次のような情報を入力して算出することができる。Ａ）情報取得部２に、同一時系列における心拍間隔の変動の瞬時位相と呼吸パターンの瞬時位相との瞬時位相差を入力し、情報処理部３が位相コヒーレンス算出機能によって入力された瞬時位相差を用いて位相コヒーレンスを算出する。Ｂ）情報取得部２に、同一時系列における心拍間隔の変動の瞬時位相及び呼吸パターンの瞬時位相を入力し、情報処理部３が、両者の瞬時位相差を算出する瞬時位相差算出機能を有し、瞬時位相差算出機能によって両者の瞬時位相差を算出し、算出された瞬時位相差を用いて位相コヒーレンス算出機能によって位相コヒーレンスを算出する。Ｃ）情報取得部２に、同一時系列における心拍間隔の変動及び呼吸パターンを入力し、情報処理部３が、瞬時位相算出機能を有し、瞬時位相算出機能によって心拍間隔の変動の瞬時位相と呼吸パターンの瞬時位相を算出し、算出された瞬時位相を用いて瞬時位相差算出機能及び位相コヒーレンス算出機能によって位相コヒーレンスを算出する。Ｄ）情報取得部２に、心拍に関する情報を含んだ生体情報及び呼吸に関する情報を含んだ生体情報を入力し、情報処理部３が、心拍間隔算出機能及び呼吸パターン算出機能を有し、心拍間隔算出機能によって心拍の情報を含んだ生体情報から心拍間隔の変動を算出し、呼吸パターン算出機能によって呼吸パターンの情報を含んだ生体情報から呼吸パターンを算出し、その後、上記Ｃ）と同様の処理を行う。Ｅ）情報取得部２に、心拍に関する情報及び呼吸に関する情報の両方を含む生体情報を入力し、情報処理部３が、かかる生体情報から心拍の情報又は呼吸パターンの情報を検出又は抽出する機能を有し、検出又は抽出した心拍の情報又は呼吸パターンの情報を用いてその後の処理をしてもよい。【0046】
さらに、情報処理部３は、算出した位相コヒーレンスに基づいて、睡眠状態等を判定する判定機能を有していてもよい。判定機能としては、例えば、算出した位相コヒーレンスを閾値と比較して、閾値よりも大きい場合は深い睡眠であると判定し、小さい場合は浅い睡眠と判定してもよいし、位相コヒーレンスが閾値よりも大きい値だった時間で睡眠の品質を評価してもよいし、位相コヒーレンスの変動の周期によって睡眠の品質を評価してもよい。閾値は、予め定めた数値であってもよいし、計測対象の過去に算出した位相コヒーレンスの数値から特定してもよいし、複数の数値を設定し、段階的に睡眠の質を評価してもよい。さらに、情報処理部３の判定機能によって、睡眠中の呼吸状態を判定することもできる。例えば、呼吸パターンの取得方法にもよるが、いずれの取得方法であっても、中枢性睡眠時無呼吸（脳の呼吸中枢の異常等により呼吸運動が停止して起こる無呼吸）については、呼吸運動が停止するので呼吸パターンが検出できなくなるので、無呼吸状態であることを判定できる。さらに、呼吸による空気の流れを計測して呼吸パターンを取得した場合は、呼吸の気流が停止して呼吸パターンが検出されなくなるので、中枢性睡眠時無呼吸だけではなく、閉塞性睡眠時無呼吸（呼吸運動はあるが気道の閉塞等による無呼吸）も判定できる。【0047】
操作部４は、使用者が睡眠状態測定装置１を操作するためのスイッチ、タッチパネル、ボタン、つまみ、キーボード、マウス、音声入力用マイク等の操作端子が設けられている。また操作部４には、操作内容等を表示するディスプレイが設けられていてもよい。出力部５は、算出した位相コヒーレンスを出力してもよいし、位相コヒーレンス以外の生体情報を出力してもよいし、判定機能で判定した睡眠状態等を出力してもよい。出力部５としては、結果を画像で表示するディスプレイ、結果を紙で出力するプリンター、結果を音声で出力するスピーカー、結果を電子情報で出力する有線又は無線の出力端子などを使用することができる。なお、出力部５としてのディスプレイを操作部４におけるタッチパネルや操作内容等を表示するディスプレイと兼用させる構成であってもよい。記憶部６は、情報取得部２で取得した情報や、情報処理部３で算出した結果、判定機能で判定した結果などを記憶することができる。【0048】
図３は、睡眠状態測定装置１の一例である。睡眠状態測定装置１は、センサ２１、アナログ−デジタル変換回路３１、心拍抽出手段３２、呼吸波形抽出手段３３、心拍間隔算出手段３４、ヒルベルト変換フィルタ３５、３６、瞬時位相差算出手段３７、位相コヒーレンス算出手段３８、操作ボタン４１、タッチパネル４２、音声入力用マイク４３、表示ディスプレイ４４、無線通信手段４５、スピーカー４６、記録装置４７を有している。【0049】
センサ２１は、心拍に関する情報を含んだ生体情報及び呼吸に関する情報を含んだ生体情報を検出するものである。例えば、心拍に関する情報を含んだ生体情報を検出するものとして心電図計測用センサ又は振動を計測するセンサなどがあり、呼吸に関する情報を含んだ生体情報を検出するものとして心電図計測用センサ、振動を計測するセンサ又は呼吸センサなどがある。図３においては１つであるが、複数種類のセンサを含んでいてもよい。心電図計測用センサ又は振動を計測するセンサは、心拍に関する生体情報も呼吸に関する生体情報も含む信号を検出できるものであり、１つのセンサで両方検出することが可能であるが、何れかの生体情報を別のセンサで検出してもよい。心電図計測用センサの場合は、ディスポーザブル電極を用いて専用の電子回路を生体の胸部に貼付して計測することが好ましく、電極によって心電図波形が計測される。導出法は単極誘導または双極誘導でも良い。【0050】
また、振動を計測するセンサは、接触式でも非接触式でもよく、接触型の振動を計測するセンサの場合は、動物に直接又は間接的に接触させて配置することによって、心弾動図波形又は生体振動信号を検出することができる。心弾動図波形又は生体振動信号を検出するための接触型の振動を計測するセンサは、振動を発生する種々の生物に直接又は近傍に配置され、生物からの振動を検出し電気信号として出力できれば足りる。振動を計測するセンサとしては、圧電センサとしてピエゾ素子が好適に用いられるが、その他のセンサ、例えば高分子圧電体（ポリオレフィン系材料）を用いてもよい。ピエゾ素子の素材としては、例えば、多孔性ポリプロピレンエレクトレットフィルム（ＥｌｅｃｔｒｏＭｅｃｈａｎｉｃａｌ Ｆｉｌｍ（ＥＭＦＩ））、またはＰＶＤＦ（ポリフッ化ビニリデンフィルム）、またはポリフッ化ビニリデンと三フッ化エチレン共重合体（Ｐ（ＶＤＦ−ＴｒＦＥ））、又はポリフッ化ビニリデンと四フッ化エチレン共重合体（Ｐ（ＶＤＦ−ＴＦＥ））を用いてもよい。圧電センサとしては、フィルム状であることが好ましい。さらに、圧電センサの場合、動物を拘束せずに心弾動図波形又は生体振動信号を取得することが可能であり、よりストレスフリーで測定できるので好ましい。ただし、圧電センサは、リストバンド、ベルト、腕時計、指輪、ヘッドバンド等に取り付けて、動物に装着してウェアラブルセンサとして利用することもできる。また、その他の種類の振動を計測するセンサとして、例えば、高感度の加速度センサを用いて、腕時計、携帯端末のように体と接触させて、あるいはベッド、椅子等の一部に加速度センサを設置して心弾動図波形又は生体振動信号を取得してもよいし、チューブ内の空気圧又は液体圧の変化を圧力センサ等で検知して、心弾動図波形又は生体振動信号を取得してもよい。さらに、振動を計測するセンサとして、マイクロ波等を用いた信号受発信に伴って非接触で心弾動図波形又は生体振動信号を取得できる非接触式のセンサを利用してもよい。マイクロ波としてはマイクロ波ドップラーセンサ、ＵＷＢ（ウルトラワイドバンド）インパルスの反射遅延時間を測定し、対象物との距離を測定する受信波による心弾動図波形又は生体振動信号、マイクロ波以外の電磁波を用いて得られた心弾動図波形又は生体振動信号、ＬＥＤ光を使った反射又は透過光から得られる心弾動図波形又は生体振動信号、さらには、超音波の反射波から得られる心弾動図波形又は生体振動信号を使用してもよい。これらのマイクロ波等を用いたセンサは、小型化が可能であり、非接触かつ非拘束で信号を取得でき、遠隔から信号を取得できる。なお、加速度センサも小型化が可能である。【0051】
呼吸センサの場合は、例えば、呼吸による空気の流れを計測した実測値、呼吸に伴う胸郭インピーダンス変化を計測した実測値、呼吸による温度変化を計測した実測値、呼吸運動に伴う腹部の動きを計測した実測値等によって呼吸パターンの時間変化が計測される。センサ２１で検出した信号は、有線又は無線で睡眠状態測定装置１のアナログ−デジタル変換回路３１に入力される。図４（Ａ）の上側のグラフは、心電図計測用センサで計測された心電図波形である。【0052】
アナログ−デジタル変換回路３１は、センサ２１からのアナログ信号をデジタル信号に変換する回路である。センサ２１内にアナログ−デジタル変換回路３１を設けてもよいし、センサ２１がデジタル信号を検出できる場合には設けなくてもよい。また、センサ２１からのアナログ信号をフィルタリング等の処理を行った後にアナログ−デジタル変換回路３１によってデジタル信号に変換してもよい。【0053】
心拍抽出手段３２は、センサ２１で検出した信号から心拍に関する信号を抽出する手段であり、センサの種類又は入力される信号に応じて適宜適当な処理が選択される。心電図波形や心弾動図波形が入力された場合、通常、心電図波形や心弾動図波形には呼吸の影響を受けているため、呼吸の成分を取り除くための処理を行うことが好ましいが、心拍間隔の算出に問題がなければ心拍抽出手段を使用しなくてもよい。また、生体振動信号が入力された場合、通常、生体振動信号には、心臓の拍動による心弾動だけではなく、呼吸による振動や、体動、発声、外部環境等に基づく振動も含まれる場合があり、これらのノイズを除去する処理を行うことが好ましい。かかる処理としては、例えば、心電図波形や振動信号の強度をｎ乗（ｎは２以上の整数であり、ｎが奇数の場合は絶対値を取る）して強調処理した後、バンドパスフィルタ（ＢＰＦ）を通過させてもよい。心拍抽出手段３２のＢＰＦは、通過域の下限周波数が０．５Ｈｚ以上、０．６Ｈｚ以上、０．７Ｈｚ以上、０．８Ｈｚ以上、０．９Ｈｚ又は１Ｈｚ以上であることが好ましく、上限周波数が１０Ｈｚ以下、８Ｈｚ以下、６Ｈｚ以下、５Ｈｚ以下、３Ｈｚ以下であることが好ましく、これらの下限周波数の何れかと上限周波数の何れかを組み合わせた通過域を持つことが好ましい。心拍抽出手段３２の下限周波数が、呼吸波形抽出手段３３の上限周波数と同じであってもよいし、呼吸波形抽出手段３３の上限周波数よりも低く、心拍抽出手段３２の通過域の一部が呼吸波形抽出手段３３の通過域と重畳していてもよい。また、心拍間隔抽出方法として、取得した心弾動図波形又は生体振動信号からフィルタの上限周波数又は下限周波数を求めることが好ましく、さらに好ましくは、定期的又は不定期に、取得した心弾動図波形又は生体振動信号からフィルタの上限周波数又は下限周波数を求めることが好ましい。例えば、心弾動図波形又は生体振動信号もしくはこれらの信号に前処理（例えば、ノイズ除去、強調処理等）したもの（心弾動図波形又は生体振動信号に由来する信号）について、パワースペクトルを求め、０．５Ｈｚ以上から密度を検索して最初のピークを同定し、そのピークが所定の閾値（例えばピークの半値幅）まで低下する低周波側および／又は高周波側の周波数の帯域を通過周波数としても良い。パワースペクトルは、例えばフーリエ変換することにより求めることができる。このように、取得した心弾動図波形又は生体振動信号から求めた上限周波数又は下限周波数のフィルタを用いて信号処理を行うことにより、取得した生体に特有の生体情報や取得時の体勢、体調、環境等の条件が反映され、個人差や取得時の条件に対応したフィルタを設定することができ、リアルタイムで位相コヒーレンスを算出できた。なお、センサ２１の一部に呼吸センサを使用していた場合には、呼吸センサからの信号は、心拍抽出手段３２に入力する必要はない。【0054】
さらに、心弾動図波形から心拍間隔を算出する方法として、心弾動図波形を心拍の周波数よりも高い周波数を下限周波数とするハイパスフィルタ（ＨＰＦ）を通過させ、ＨＰＦ後の信号について絶対値をとることにより、利用したＨＰＦを通過した信号の包絡線信号から、心弾動図波形の各心拍動の山を得ることができ、そのピーク値または心拍動の山の開始時点から心拍間隔を求めることができる。通常の心拍の周波数は最大でも３Ｈｚ程度であるが、このハイパスフィルタの下限周波数は、５Ｈｚ以上であることが好ましく、１０Ｈｚ、２０Ｈｚ、３０Ｈｚ、４０Ｈｚであってもよい。この信号処理方法で心拍間隔の変動を算出して得られた位相コヒーレンス（λ）の値は、心電図波形から求めた位相コヒーレンス（λ）の値に非常に近い値を得ることができた。また、ＨＰＦ後の信号について絶対値をとった信号について、上記のパワースペクトルから求めた通過周波数のＢＰＦ（またはＬＰＦ）を利用することがより好ましい。さらに、ＨＰＦを通過させる前に、ノイズ除去の前処理等を行ってもよい。生体振動信号（心弾動図波形を含む）は心拍動に伴う振動波形を含むものであるが、呼吸運動に伴う振動成分が心拍動成分と重畳すると波形が安定せず、原波形からＲＲＩに相当する一拍毎の拍動間隔を得るのは従来難しかった。本手法では心拍動の基本周波数成分と重畳した呼吸周波数成分を予め除去した後、心拍動由来の高周波振動成分から拍動間隔を求めることにより、より正確な拍動間隔を得ることができ、その結果、心拍間隔の変動である呼吸性不整脈の検出も正確になり、算出した位相コヒーレンスは心電図と実測呼吸から求めたものとほぼ一致する。【0055】
呼吸波形抽出手段３３は、センサ２１で検出した信号から呼吸パターンに関する信号を抽出する手段であり、センサの種類又は入力される信号に応じて適宜適当な処理が選択される。センサ２１の一部に呼吸センサを使用し、呼吸センサによって呼吸パターンが実測される場合には、呼吸波形抽出手段３３を設けなくてもよいし、呼吸波形抽出手段３３でノイズとなる信号を除去する処理を行ってもよい。心電図波形、心弾動図波形又は生体振動信号から呼吸パターンに関する信号を抽出する場合には、かかる処理としては、例えば、心電図波形や振動信号の強度をｎ乗（ｎは２以上の整数であり、ｎが奇数の場合は絶対値を取る）して強調処理した後、０．５Ｈｚ以下の周波数範囲の通過域を有するローパスフィルタ（ＬＰＦ）を通過させてもよい。呼吸波形抽出手段３３のＬＰＦの遮断周波数は、０．３、０．４、０．６、０．７Ｈｚ、０．８Ｈｚであってもよい。また、呼吸波形抽出手段３３の遮断周波数は、心拍抽出手段３２の下限周波数と同じであってもよいし、下限周波数よりも高くして通過域の一部が重畳していてもよい。また、呼吸波形抽出方法として、取得した心電図波形、心弾動図波形又は生体振動信号からフィルタの上限周波数又は下限周波数を求めることが好ましく、さらに好ましくは、定期的又は不定期に、取得した心弾動図波形又は生体振動信号からフィルタの上限周波数又は下限周波数を求めることが好ましい。心電図波形、心弾動図波形又は生体振動信号もしくはこれらの信号に前処理（例えばノイズ除去のフィルタや強調処理等）したものについて、パワースペクトルを求め、低周波側からパワースペクトル密度を検索して最初のピークを同定し、そのピークが所定の閾値（例えばピークの半値幅）まで低下する高周波側の周波数を遮断周波数としても良い。このように、取得した心電図波形、心弾動図波形又は生体振動信号から求めた上限周波数又は下限周波数のフィルタを用いて信号処理を行うことにより、取得した生体に特有の生体情報や取得時の体勢、体調、環境等の条件が反映され、個人差や取得時の条件に対応したフィルタを設定することができ、リアルタイムで位相コヒーレンスを算出できた。また、ＬＰＦの代わりにＢＰＦを通過させても良く、この場合、ＢＰＦの下限周波数は十分低い周波数であれば足り、例えば０．１Ｈｚに設定してもよい。図４（Ａ）の下側のグラフは、呼吸パターンを示すものであり、実線は心電図波形から抽出した呼吸パターンであり、点線は呼吸による空気の流れを計測した実測値である。図４（Ａ）の下側のグラフより、心電図波形から抽出した呼吸パターンでも、周期は実測値と一致していることが確認できる。【0056】
図５は、８名の被験者で呼吸数を１分間に８、１０、１２、１５、１８、２０、２４回で変化させた時の心電図から抽出した推定呼吸周波数と、実測呼吸周波数の相関関係を示した図である。点線は呼吸周波数の９５％信頼区間を示している。呼吸周波数が約０．４Ｈｚを超えると９５％信頼区間がｉｄｅｎｔｉｔｙ ｌｉｎｅから外れ、過小評価されるが、０．４Ｈｚ未満では１０％以内の精度で呼吸周波数を推定できる。図６は、図５の推定呼吸周波数を用いて算出した位相コヒーレンスλecgと実測呼吸周波数を用いて算出した位相コヒーレンスλの相関関係を示した図である。呼吸周波数が０．３３Ｈｚ（２０回／分）までは、９割以上の精度で位相コヒーレンスを求められた。【0057】
心拍間隔算出手段３４は、心拍抽出手段３２からの信号が入力され、心拍の間隔を算出する。心拍の間隔は、例えば心電図のＰ波、Ｒ波、Ｔ波又はＵ波の間隔、特にＲ波が鋭いピークを有するので、Ｒ波から次のＲ波までの間隔を計測することが好ましい。心弾動図波形や生体振動信号から抽出された心拍に関する信号の場合も、鋭いピークのＲ波に相当する波形の間隔を計測することが好ましい。図４（Ｂ）は、図４（Ａ）の心電図から算出した心拍間隔の変動を示すグラフであり、縦軸が心拍間隔（ｍｓ）、横軸が時間（ｓ）である。図４（Ｂ）から、心拍間隔が一定の周期で変動していることが確認できる。なお、呼吸性不整脈（ＲＳＡ）の振幅も、心理ストレス等を評価する指標の一つとして利用可能であるが、後述するように、呼吸周波数によっても呼吸性不整脈（ＲＳＡ）の振幅が変化するので、位相コヒーレンスによる評価と組み合わせて補助的又は追加的に評価するのが好ましい。また、心拍間隔算出手段３４は、心拍に関する情報を含んだ生体情報から直接心拍の間隔を算出してもよい。この場合、心拍抽出手段３２の機能を含む心拍間隔算出手段３４であってもよいし、信号処理方法によっては、心拍抽出手段３２を必要とせず、心拍に関する情報を含んだ生体情報から直接心拍の間隔を算出できる。【0058】
ヒルベルト変換フィルタ３５は、心拍間隔の変動について瞬時位相と瞬時振幅を出力するものであり、ヒルベルト変換フィルタ３６は、呼吸パターンについて瞬時位相と瞬時振幅を出力するものである。ヒルベルト変換は、アナログ回路で９０度位相差分波器を実現しても良いし，有限インパルス応答型のデジタルフィルタで構成しても良い。ヒルベルト変換した信号と実信号を加えて解析信号を得て、解析信号の実部と虚部の比から瞬時位相を求めることができる。図４（Ｃ）の上側のグラフは、実線が心拍間隔の変動の瞬時位相であり、点線が呼吸パターンの瞬時位相であり、縦軸が位相（ラジアン）であり、横軸が時間（ｓ）である。【0059】
瞬時位相差算出手段３７は、心拍間隔の変動の瞬時位相と呼吸パターンの瞬時位相との位相差（瞬時位相差）を算出し、結果を位相コヒーレンス算出手段３８に出力する。位相コヒーレンス算出手段３８では、上記のとおり、心拍間隔の変動の瞬時位相と呼吸パターンの瞬時位相との瞬時位相差を用いて位相コヒーレンスを算出する。位相コヒーレンスを求める際のデータは最低でも呼吸１周期の窓長で計算する。図４（Ｃ）の下側のグラフは、算出した位相コヒーレンスλである。図４（Ｃ）では、常に位相コヒーレンスが１に近く、比較的ばらつきが少ない状態であることが確認できる。睡眠状態測定装置１は、さらに算出した位相コヒーレンスλから睡眠状態を判定する機能を有していてもよい。【0060】
操作ボタン４１、タッチパネル４２、音声入力用マイク４３は、使用者が睡眠状態測定装置１を操作するための入力手段であり、睡眠状態測定装置１を作動させたり、必要な情報を出力させたりすることができる。表示ディスプレイ４４、スピーカー４６は、心拍、呼吸パターン、呼吸性不整脈、位相コヒーレンスλや、位相コヒーレンスλから推定される睡眠状態などを出力する出力手段として利用することができる。無線通信手段４５は、算出した位相コヒーレンスλや睡眠状態の出力手段に利用してもよいし、センサ２０からの信号を入力する入力手段として使用してもよい。または、音声で睡眠状態などを出力してもよい。記録装置４７は、入力された情報、各種手段のプログラム、計測結果などが記録される。【0061】
睡眠状態測定装置１は、携帯端末（たとえば、携帯電話、スマートフォン等）とセンサで実現することもできる。センサ、例えば心電図センサにＡ／Ｄ変換回路と無線通信機能を設け、センサで検出した信号をＡ／Ｄ変換回路でデジタル信号に変換し、デジタル信号を無線通信機能によって携帯端末に送信する。無線通信機能としては、例えばＢｌｕｅｔｏｏｔｈ（登録商標）、Ｗｉ−ｆｉ（登録商標）などを利用することが好ましい。【0062】
本発明の睡眠状態測定装置１は、心拍間隔の変動の瞬時位相と呼吸パターンの瞬時位相との瞬時位相差の位相コヒーレンスが、δ波の振幅と相関関係があることを利用して睡眠状態を測定することができる。位相コヒーレンスは、後述するように、呼吸周波数の影響を受けにくく、より正確に睡眠状態を測定することができる。また、一つのセンサでも睡眠状態測定装置１を実現することが可能であり、使用者に負担が少なくストレスを感じにくい装置を提供できる。特に振動を計測するセンサの場合は動物を拘束せずに計測することが可能であり、より簡便な装置を提供できる。位相コヒーレンスは、リアルタイムに測定することができ、ノンレム睡眠時間の発生時点やノンレム―レム睡眠のリズム周期を測定することができ、睡眠状態を的確に把握することができる。また、同時に呼吸パターンを測定又は抽出しているので、睡眠時の無呼吸状態を検出することも可能である。さらに、被測定者の呼吸性不整脈の大きさも測定できるので、睡眠状態について補助的な判断指標とすることも可能である。【0063】
［位相コヒーレンス算出装置及び方法］位相コヒーレンスは、睡眠状態だけではなく、心理ストレスの評価にも利用することができる。特に、本実施形態で説明する位相コヒーレンス算出装置は、少なくとも動物の心拍に関する情報及び呼吸に関する情報の両方を含む生体情報から、心拍間隔の変動と呼吸パターンを取得して、位相コヒーレンスを算出する。
【0064】
図７は、位相コヒーレンス算出装置１１の概略ブロック図である。位相コヒーレンス算出装置１１は、少なくとも生体情報取得手段１２と、呼吸波形抽出手段１３と、心拍間隔算出手段１４と、位相コヒーレンス算出手段１５とを含んでいる。さらに、位相コヒーレンス算出装置１１は、図３に記載されたその他の手段を具備していてもよい。本実施形態における位相コヒーレンス算出装置１１は、生体情報取得手段１２が、少なくとも動物の心拍に関する情報及び呼吸に関する情報の両方を含む生体情報を取得するものであり、かかる生体情報として、例えば、心電図、心弾動図波形又は生体振動信号が挙げられる。生体情報取得手段１２は、動物を計測するためのセンサ及びセンサの情報を有線又は無線で入力する入力部を含む構成であってもよいし、すでに計測済みの情報が記録された他の記録媒体からの情報を有線又は無線で入力可能な入力部を含む構成であってもよい。すなわち、情報取得部２は、少なくとも情報を入力する入力部を備えており、場合によっては入力部と有線又は無線で接続された生体情報を計測するためのセンサを備えていてもよい。センサで生体情報を計測する場合、サンプリング周波数は１００Ｈｚ以上であることが好ましい。【0065】
心電図から呼吸パターンを抽出し、心拍間隔の変動と呼吸パターンを取得できることは上記睡眠状態測定装置１の説明において詳述したとおりである。心弾動図波形又は生体振動信号の場合、振動を計測するセンサは、生体に直接接触させてもよいが、生体から振動が伝搬する部材（床、ベッド、椅子、机、衣服、靴、絨毯、シーツ、カバーなどを含む）に接触させてもよいし、非接触式の場合は生体から離間して配置してもよい。接触式の振動を計測するセンサが設置される部材は、生体から複数の部材を介在させたものでもよい。例えば、床の上、ベッドマットの上又は下、椅子の座板や背もたれの表面又は裏面、机の天板の表面又は裏面等に振動を計測するセンサを設置又は埋設してもよい。さらに、ベッドの脚において、床と接触する部位に設けられた保護部材の一つに、振動を計測するセンサを設けてもよいし、ベッドのフレーム、ヘッドボード又はサイドレール等に設けてもよいし、椅子の脚、ひじ掛け、フレーム等に設けてもよいし、机の脚、貫、幕板等に設けてもよい。また、トイレの便座又は便器に振動を計測するセンサを設けてもよく、便座又は便器の表面、裏側又は内部に振動を計測するセンサを配置することができる。例えば、便座の裏面における便器との接触箇所に設けられた緩衝部又は便器の表面における便座との接触箇所に振動を計測するセンサ部を配置してもよい。なお、振動を計測するセンサとしては、上記のとおり、圧電センサ、加速度センサ、圧力センサ、非接触式のセンサ等を使用することができる。【0066】
心弾動図波形又は生体振動信号の場合、心拍に関する情報と、呼吸に関する情報をそれぞれ抽出する必要がある。ローパスフィルタ（ＬＰＦ）、バンドパスフィルタ（ＢＰＦ）、ハイパスフィルタ（ＨＰＦ）によって周波数で心拍に関する情報と、呼吸に関する情報を抽出することができる。より正確に心弾動図波形又は生体振動信号から心拍を求めるために、心弾動の伝達特性を算出し、逆伝達関数を推定することで心弾動図波形又は生体振動信号から心電図に相当する波形を求めることができる。逆伝達関数は、予め被測定者の心電図と心弾動とを測定し、伝達特性を調査してもよいが、心弾動図波形又は生体振動信号のみから伝達特性を推定することも可能である。心弾動図波形から低周波成分を除去した後にWiener フィルタなどを適用し、逆伝達関数と元の心弾動図波形を重畳積分して模擬心電図波形を得ることができる。また、心弾動図波形から心拍間隔を算出する他の方法として、心弾動図波形を心拍の周波数よりも高い周波数を下限周波数とするハイパスフィルタ（ＨＰＦ）を通過させ、ＨＰＦ後の信号について絶対値をとることにより、利用したＨＰＦを通過した信号の包絡線信号から、心弾動図波形の各心拍動の山を得ることができ、そのピーク値または心拍動の山の開始時点から心拍間隔を求めてもよい。【0067】
図８（Ａ）は、上側のグラフが心弾動図波形であり、下側のグラフが心弾動図波形から抽出された模擬心電図波形である。まず、心弾動図波形においても、鋭いピークが周期的に生じているので、そのピークの時刻（Ｔ）を求め、その時刻が心電図から得られるＲ波時刻と仮に推定し、心弾動図波形の鋭いピークを起点として、それに伴う心弾動図波形を所定の窓長（長くとも次のピークまで）の波形断片で切り出してまとめる。まとめた際にばらつきの大きい波形断片は、時刻（Ｔ）がＲ波時刻と一致していないとして排除することが好ましい。図８（Ｂ）の実線は、１００個の波形断片を重ね合わせた際の平均であり、点線は標準偏差である。この実線を心電図と心弾動間の平均伝達特性とし、かかる伝達特性の逆伝達関数を用いて心弾動図波形から模擬心電図波形を得ることができた。このように、心弾動の伝達特性を算出し、逆伝達関数を推定することで心弾動図波形又は生体振動信号から心電図に相当する波形を求めることでより正確に抽出することができる。特に、本発明の位相コヒーレンス算出装置及び睡眠状態測定装置は、心拍間隔の変動や、呼吸パターンの抽出を行うものであり、より正確に心電図に相当する波形を求めることは重要である。【0068】
また、心弾動図波形から呼吸に関する信号は、ローパスフィルタ（ＬＰＦ）を通過させることで抽出したり、模擬心電図波形から呼吸由来の振幅変調を抽出したりすることで抽出できる。【0069】
センサで生体情報を計測する場合、サンプリング周波数は１００Ｈｚ以上であることが好ましい。図９は、心電図波形を取得する際のサンプリング周波数として、１ｋＨｚ、５００Ｈｚ、２００Ｈｚ、１００Ｈｚ、５０Ｈｚと変化させた場合の位相コヒーレンスを比較した図である。１ｋＨｚでサンプリングした心電図の位相コヒーレンスと比較した二乗平均平方根誤差（Root Mean Square Error:ＲＭＳＥ）は、５００Ｈｚが０．０２８、２００Ｈｚが０．０３９、１００Ｈｚが０．０４５、５０Ｈｚが０．１０９であった。このように、１００Ｈｚ以上のサンプリング周波数であれば、十分精度の高い位相コヒーレンスを得ることができる。【0070】
図１０は、安静時と暗算課題時（ストレス状態）の心拍間隔の変動と呼吸パターンの瞬時位相の関係を示すものである。図１０（Ａ）は、安静時における心拍（実線）と呼吸（点線）の瞬時位相であり、（Ｂ）は、安静時の心拍と呼吸の瞬時位相のリサージュ図である。図１０（Ｃ）は、暗算課題時（ストレス状態）における心拍（実線）と呼吸（点線）の瞬時位相であり、（Ｄ）は、暗算課題時の心拍と呼吸の瞬時位相のリサージュ図である。図１０（Ｅ）は、安静状態から暗算課題を課したとき、および暗算課題終了後の呼吸性不整脈（ＲＳＡ）の変化である。図１０（Ｆ）は、安静状態から暗算課題を課したとき，および暗算課題終了後の位相コヒーレンスの変化であり、点線は、呼吸流速による実測の呼吸パターンを用いて算出された位相コヒーレンスであり、実線は、心電図から算出した呼吸パターンを用いて算出された位相コヒーレンスである。安静時の位相コヒーレンスは、０．６９±０．１２（９５％信頼区間：０．６３〜０．７５）であり、暗算課題時の位相コヒーレンスは、０．４５±０．１７（９５％信頼区間：０．４１〜０．４９）であり、安静時に比べて暗算課題時は有意に位相コヒーレンスが低下していた。図１０から、瞬時位相差は、安静時では非常に安定しており、暗算課題などの精神ストレスを課すと呼吸性不整脈（ＲＳＡ）の大きさが減弱するだけではなく、位相差が乱れることが確認できる。これは呼吸中枢により生成される呼吸振動子と自律神経支配下にある心拍振動子の協調関係がストレスにより攪乱されることを意味している。【0071】
図１１は、自発的に呼吸周波数を変化させた場合の心拍間隔（ＲＲＩ）、呼吸性不整脈の振幅（ＡＲＳＡ）、位相コヒーレンス（λ）及び呼吸周波数（ｆR）を示す図である。最初、呼吸数は１分当たり１５回であったが、３００ｓから、意識的に１分当たり２５回に増やした。呼吸数が増えると呼吸性不整脈の振幅（ＡRSA）は小さくなるが、位相コヒーレンスはほとんど変化しない。呼吸性不整脈は自律神経の副交感神経活動に由来する生理現象であり、その振幅は副交感神経活動の緊張度を反映すると言われているが、図１１のように、副交感神経活動の緊張度を変化させるのではなく、呼吸数を変化させただけでも呼吸性不整脈の振幅に影響するため、必ずしも呼吸性不整脈から自律神経活動を推定することはできなかった。この点、位相コヒーレンスは呼吸数を変化させても変化しないので、位相コヒーレンスから自律神経活動を推定することでより正確に自律神経活動を推定できる。このように，位相コヒーレンスを測定することで、呼吸周波数に依存せずにリアルタイムで被測定者の心理ストレス状態や睡眠状態を確認することが可能となる。さらに、被測定者の呼吸性不整脈の大きさも測定できるので、心理ストレス状態や睡眠状態について補助的な判断指標とすることができる。さらに、心拍間隔の変動も情報として取得することができるので、心拍間隔の変動の時間変化を周波数分析（フーリエ変換）することにより、遅いゆらぎの成分（ＬＦ）及び速いゆらぎの成分（ＨＦ）を求め、ＬＨ／ＨＦを交感神経活動指標とし、ＨＦを副交感神経活動指標として検出することも可能である。例えば、心拍間隔の変動の時間変化をフーリエ変換したパワースペクトルの０．０４〜０．１５Ｈｚの成分をＬＦとし、０．１５〜０．４Ｈｚの成分をＨＦとして算出してもよい。これらの指標は、従来、心電図から得られたＲ波の心拍間隔の変動解析として用いられていたが、本発明においては、心電図のみではなく、心弾動図波形又は生体振動信号から得られた心拍間隔の変動においても、解析手法として用いることができ、交感神経・副交感神経活動の調節障害を調べることが可能となった。なお、本発明の睡眠状態測定装置において心電図を利用した場合に、これらの指標を同時に検出してもよい。【0072】
動物の心拍に関する情報を含んだ生体情報又は呼吸に関する情報を含んだ生体情報の信号処理方法として、生体情報に由来する信号のパワースペクトルから上限周波数又は下限周波数を求め、前記上限周波数又は下限周波数を遮断周波数とするフィルタを通過させる処理を含むことが好ましい。生体情報に由来する信号のパワースペクトルから上限周波数又は下限周波数を求めることにより、取得した生体に特有の生体情報や、生体情報を取得する際の体勢、体調、環境等の条件が反映され、個人差や取得時の条件に対応したフィルタを設定することができる。さらに、体勢、体調、環境等の条件は常に変化するので、定期的又は不定期に、フィルタの上限周波数又は下限周波数を更新することが好ましい。【0073】
心拍に関する情報を含んだ生体情報とは、例えば、心電図、心弾動図波形、動物の振動（心弾動を含む）の時間的変化を計測した生体振動信号等が利用可能であり、呼吸に関する情報を含んだ生体情報とは、例えば、呼吸による空気の流れを計測した実測値、呼吸に伴う胸郭インピーダンス変化を計測した実測値、呼吸による温度変化を計測した実測値、呼吸運動に伴う腹部の動きを計測した実測値、心電図、心弾動図波形、生体振動信号等が利用可能である。生体情報に由来する信号とは、これらの心拍に関する情報を含んだ生体情報又は呼吸に関する情報を含んだ生体情報それ自体だけではなく、生体情報に前処理（例えば、ノイズ除去、強調処理等）したものを含む。かかる信号処理方法は、例えば、上記睡眠状態測定装置において、心拍間隔算出機能によって心拍の情報を含んだ生体情報から心拍間隔の変動を算出したり、呼吸パターン算出機能によって呼吸パターンの情報を含んだ生体情報から呼吸パターンを算出したりする際に利用してもよいし、位相コヒーレンス算出装置において、動物の心拍に関する情報及び呼吸に関する情報の両方を含む生体情報から、心拍間隔の変動を算出したり、呼吸パターンを算出したりする際に利用してもよい。【0074】
［生体振動信号測定装置（ウェアラブルセンサ）］図１２は、本発明の生体振動信号測定装置１０１の概略ブロック図である。生体振動信号測定装置１０１は、ヒトに装着可能な装着部に少なくとも振動センサ１０２を備え、必要に応じて、情報処理手段１０３、通信手段１０４、電力供給手段１０５、記憶手段１０６、表示出力手段１０７、操作手段１０８などの一つ又は複数を備えていてもよい。装着部は、動物の体に装着可能であり、例えば、ヒト又は動物の四肢又は頭部に装着可能なウェアラブルなものが好ましい。ヒト又は動物の上肢としては、手指、手首、腕などに装着することが好ましく、例えば、指輪、腕輪、指サック、リストバンドなどにセンサを実装することができる。また、ヒト又は動物の下肢としては、腿、脛、足首に装着することが好ましく、例えば、バンド、靴下、スパッツなどにセンサを実装してもよい。また、ヒト又は動物の頭部においても、首、こめかみ、耳などに装着することが好ましく、例えば、ヘッドバンド、ネクタイ、ネックレス、ピアス、イアリングなどにセンサを実装してもよい。
【0075】
振動センサ１０２は、振動を計測するセンサであり、心拍に関する情報と呼吸に関する情報を含んだ生体振動信号を取得し、装着部を介してヒトに装着可能とされている。図１２においては、１つの振動センサ１０２が例示されているが、複数の振動センサを含んでもよいし、必要に応じて、他の種類のセンサ（例えば、光学センサ、温度センサ）を振動センサ１０２と合わせて利用してもよい。振動センサ１０２で生体振動信号を計測する場合、サンプリング周波数は１００Ｈｚ以上であることが好ましい。振動センサ１０２は、ヒトに直接又は間接的に接触させて配置することによって、生体振動信号（心弾動図波形を含んでもよい）を検出することができる。振動センサ１０２としては、上記のとおり、ピエゾ素子、高分子圧電体などを用いることができる。振動センサ１０２で検出された生体振動信号は、装置内のバス回路又は通信手段１０４を介して情報処理手段１０３等に伝達される。【0076】
情報処理手段１０３は、入力された生体振動信号を処理する手段であり、例えば、電子回路や、ＣＰＵ（中央処理装置）の演算処理機能を利用することができる。情報処理手段１０３は、生体振動信号測定装置１０１に設けられていることが好ましいが、生体振動信号測定装置１０１とは別に、通信手段１０４を介して生体振動信号測定装置１０１から生体振動信号を伝達可能に設けられていてもよい。ＣＰＵの演算処理機能によれば、例えば、デジタルフィルタを構成し、周波数フィルタリングを実現することもできる。また、情報処理手段１０３は、デジタル回路ではなくアナログ回路で実現することも可能である。例えば、コンデンサや抵抗及びオペアンプ等で構成されたローパスフィルタ（ＬＰＦ）やハイパスフィルタ（ＨＰＦ）などのアナログフィルタによって、周波数フィルタリングを実現してもよい。また、入力される生体振動信号がアナログ信号であれば、アナログ−デジタル変換回路によってデジタル信号に変換してもよい。【0077】
情報処理手段１０３は、少なくとも、心拍間隔算出手段１３１、呼吸波形抽出手段１３２を含み、さらに位相コヒーレンス算出手段１３３を含むことが好ましい。なお、位相コヒーレンス算出機能を有する装置を「位相コヒーレンス算出装置」ということもある。本発明の生体振動信号測定装置１０１は、位相コヒーレンス算出機能を備える場合、位相コヒーレンス算出装置の一種となる。心拍間隔算出手段１３１、呼吸波形抽出手段１３２、位相コヒーレンス算出手段１３３の具体的な処理については、上記睡眠状態測定装置１及び位相コヒーレンス算出装置１１の説明において詳述したとおりである。【0078】
通信手段１０４は、有線又は無線通信を介して各種信号を受け渡す機能を有する。通信手段１０４は、振動センサ１０２に接続された配線、ケーブルであってもよい。無線の通信手段１０４は、例えば、振動センサ１０２が取得した生体振動信号を情報処理手段１０３、記憶手段１０６、表示出力手段１０７、外部装置（図示せず）などに送信してもよいし、情報処理手段１０３が算出した位相コヒーレンスの情報を記憶手段１０６、表示出力手段１０７、外部装置（図示せず）などに送信してもよいし、記憶手段１０６に格納された生体振動信号又は位相コヒーレンスを、情報処理手段１０３、表示出力手段１０７などに送信してもよい。また、通信手段１０４は、操作手段１０８を介して使用者から入力された情報を、情報処理手段１０３、記憶手段１０６、表示出力手段１０７などに送信してもよい。通信手段１０４として無線の場合は、例えば、Ｂｌｕｅｔｏｏｔｈ（登録商標）、Ｗｉ−ｆｉ（登録商標）、近接場型の近距離無線通信（ＮＦＣ：Ｎｅａｒ ｆｉｅｌｄ ｒａｄｉｏ ｃｏｍｍｕｎｉｃａｔｉｏｎ）などを利用することが好ましい。なお、通信手段１０４は、生体振動信号測定装置１０１の態様によっては、必ずしも双方向の通信としなくてもよい。【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】
図１７に示すように、ヒトの手首には、小指側に尺骨動脈１１３が走行し、親指側に橈骨動脈１１４が走行している。振動センサの一つ１０２（１）は、尺骨動脈１１３上に置かれ、他の振動センサ２（２）は、橈骨動脈１１４上に置かれる。このため、腕輪型の生体振動信号測定装置１０１Ｄは、かかる二本の動脈から手首の生体振動信号を取得することができる。【0114】
なお、図１７に示す腕輪型の生体振動信号測定装置１０１Ｄは、情報処理手段１０３を含むように構成されているが、これに限定されない。振動センサ１０２が取得した生体振動信号を記憶する記憶手段を設け、適当な通信手段（例えば、Ｂｌｕｅｔｏｏｔｈ（登録商標）など）を介して外部の情報処理手段に受け渡すように構成してもよいし、表示手段を設けて生体振動信号又は生体情報を表示してもよい。また、図１７では、二つの振動センサ１０２が示されているが、単なる例示であって、これに限定されない。振動センサ１０２は、長い帯状に構成され、かかる振動センサ１０２が手首の全周又は全周の一部を覆うように取り付けられてもよいし、いずれか一方の振動センサ１０２だけを設けてもよい。また、情報処理手段１０３が腕輪型の生体振動信号測定装置に設けられる場合、情報処理手段１０３は、適当な通信手段などを介して信号処理の結果を外部の装置（例えば、携帯端末、コンピュータなど）に送信してもよい。信号処理の結果は、外部の装置の表示装置に表示されてもよい。振動センサ１０２として、圧電センサ（以下、圧電センサ１０２と記載することもある）を用いることが好ましい。【0115】
腕輪型の生体振動信号測定装置によれば、手首の生体振動信号を二箇所で測定するので、二つの生体信号から、心拍間隔、呼吸性不整脈、位相コヒーレンスを算出することができ、睡眠状態又はストレス状態をより高い精度で評価することができる。ただし、振動センサを使用する場合、手首を一定の強さで締め付ける必要があり、ヒトに違和感を与える場合もある。また、締め付けが足りず、振動センサが皮膚に密着しない場合は、振動センサと手首表面との間に隙間ができ、取得した生体振動信号に雑音が混ざる場合がある。さらに、手首においては、振動センサが尺骨動脈１１３又は橈骨動脈１１４の少なくとも一方の上に配置されないと、振動を正確に測定することは難しく、振動センサ１０２の配置が限定されていた。この点、指輪型の生体振動信号測定装置は、圧電センサを実装するのに好適であり、手指の周りを一定の強さで締め付けたとしても、ヒトに違和感を与えることが少なく、雑音も生じにくいので好ましい。【0116】
図１８は、図１２、図１６、図１７における振動センサ１０２の一例を示す図であり、シート厚を断面視したものである。振動センサ１０２は、シート状の振動センサ素材１２１を挟んで上下に信号電極層１２２及び接地電極層１２３を有し、信号電極層１２２及び接地電極層１２３にはそれぞれ上部引出電極１２４及び下部引出電極１２５が接続される。振動センサ１０２の少なくとも一部は、被覆カバー１２０で保護されている。振動センサ素材１２１としては、フッ素系の有機薄膜強誘電体材料であるポリフッ化ビニリデン（ＰＶＤＦ）を使用することができる。圧電センサ１０２は、例えば、ピエゾ素子であり、機械的な力（僅かな力）が印加されると、振動センサ素材１２１に起電力が発生し、振動センサ素材１２１に蓄積された電荷を電流電圧変換し、電気信号として取り出すことができる。【0117】
図１６に示す指輪部材１１０が導電性材料で構成されている場合、圧電センサ１０２の上部引出電極１２４及び下部引出電極１２５は、指輪部材１１０に接続され、圧電センサ１０２が取得した生体振動信号は、指輪部材１１０を介して情報処理手段１０３などに入力されてもよい。また、上部引出電極１２４及び下部引出電極１２５は、図１６又は図１７に示す情報処理手段１０３、通信手段１０４、記憶手段１０６などの入力端子に直接接続されてもよい。【0118】
［実施例４］本実施例では、指輪型の生体振動信号測定装置を被験者の人差し指付根、人差し指腹、薬指付根に装着し、手指の生体振動信号を測定し、手指の生体振動信号から心拍に関する情報を分離抽出した。振動センサは、コンデンサ型の圧電素子であり、金属電極（２０ｍｍφの電極）に１５ｍｍφの圧電体を挟んだ３層構造であり、皮膚と接触する側の金属板は円形ゴムで保護されており、ゴムを介して間接的に皮膚に接触し、振動を検出した。
【0119】
同様に、腕輪型の生体振動信号測定装置を被験者の手首に装着し、振動センサを手首中央（図１７の尺骨動脈１１３と橈骨動脈１１４の間）に配置した場合と、橈骨動脈１１４上に配置した場合に手首から得られる生体振動信号を測定した。【0120】
心臓から押し出された血液による拍動は、動脈の壁や血液を伝わって手足の末梢まで届くため、手指や手首の血管は、血液の流入によって生じる容積変化が生じ、血管運動反応を測定することによって、心臓や呼吸の動きを捉えることができる。かかる血管運動反応は、手指又は手首の表面を通じて振動波として圧電センサに伝わり、圧電センサでμＶオーダーの起電力が発生する。この信号には、心拍、呼吸等の所望する情報に加えて、障害となる雑音信号も含まれるので、その後の信号処理アルゴリズム（電子回路及びソフトウェア）を用いて心拍、呼吸等の情報を分離抽出した。また、比較のため、被験者の胸部に心電図用電極を貼付し、単極導出により心電図（ＥＣＧ）を同時に計測するともに、熱線型呼吸流速計により呼吸パターンの実測値を同時に計測した。【0121】
図１９（Ａ）は、左側から順にそれぞれ、被験者の人差し指付根、人差し指腹、薬指付根で検出した生体振動信号の原信号である。図１９（Ｂ）は、その処理後信号（全波整流積分後の心拍由来の振動成分波形を微分処理した後）の波形である。図１９（Ｂ）に示されるように、処理後信号は、ヒトの手指のいずれの部位であれ、同程度の強度で現れることが分かる。【0122】
図２０（Ａ）は、左側から順に、被験者の手首中央で検出した生体振動信号の原信号、手首親指側の橈骨動脈で検出した生体振動信号の原信号である。図２０（Ｂ）は、その処理後信号（全波整流積分後の心拍由来の振動成分波形を微分処理した後）の波形である。図２０（Ｂ）に示されるように、手首中央前の処理後信号は、橈骨動脈の処理後信号より振幅が小さく、手首の部位によって、取得可能な生体振動信号の強度が大きく異なることが分かる。【0123】
図２１（Ａ）は、手指の生体振動信号（信号処理後の波形）から得られた心拍動パルスの例を示す。図２１（Ｂ）は、比較のため測定した心電図から得られた心拍動パルスであり、図２１（Ｃ）は、比較のため測定した実測した呼吸パターン（呼吸流速波形）である。【0124】
［実施例５］本実施例では、指輪型の生体振動信号測定装置を被験者の人差し指付根に装着し、被験者を仰臥位姿勢、座位姿勢、立位姿勢の状態にして、手指の生体振動信号を測定し、手指の生体振動信号から心拍に関する情報と呼吸に関する情報を分離抽出し、位相コヒーレンスを算出した。振動センサは、実施例４と同様であり、ゴムを介して間接的に皮膚に接触し、振動を検出した。また、比較のため、被験者の胸部に心電図用電極を貼付し、単極導出により心電図（ＥＣＧ）を同時に計測するともに、熱線型呼吸流速計により呼吸パターンの実測値を同時に計測した。
【0125】
手指及び手首の生体振動信号、心電図、呼吸流速波形は、１００Ｈｚ〜１０００Ｈｚでサンプリングして保存した。生体振動信号は、心拍検出用と呼吸検出用の特定のデジタルフィルタを用いて信号処理を行った。生体振動信号に高周波バンドパスフィルタ（ハイパスフィルタでもよい）をかけた後、時定数０．１〜０．２秒で全波整流積分を行い心拍由来の振動成分のみを抽出し、この波形のピークを求めるため、全波整流積分後の心拍由来の振動成分波形を微分し、閾値を設定してピーク（心拍動パルス）を検出し、そのピークの間隔から拍動間隔を求めた。ここでは、心電図から求めた心拍間隔をＲＲＩ、生体振動信号から求めた心拍間隔をＢＢＩと称する。また生体振動信号に低周波のバンドパスフィルタ（ローパスフィルタでもよい）をかけ呼吸由来の振動成分を抽出し、呼吸パターンを推定した。ＲＲＩならびにＢＢＩをスプライン補間により１０Ｈｚで再サンプリングした。熱線型呼吸流速計により実測した呼吸パターンと生体振動信号から推定した呼吸パターンも１０Ｈｚで再サンプリングした。呼吸性不整脈と呼吸パターンをヒルベルト変換し、解析信号から瞬時位相を求め、その位相差Ψから位相コヒーレンスλを算出した。位相コヒーレンスは１０秒の計算窓で５秒ずつシフトさせながら求めた。【0126】
図２２は、「仰臥位姿勢」の被験者の人差し指付根で検出した生体振動信号から得られた各種信号αと心電図又は実測の呼吸パターンから得られた各種信号βとを比較した図である。図２２（Ａ）は、上側に人差し指付根の生体振動信号から得られた心拍動パルス（波形α）、下側に心電図から得られた心拍動パルス（波形β）を示した。図２２（Ｂ）は、人差し指付根の生体振動信号から推定した呼吸パターン（波形α）と、実測した呼吸パターン（波形β）とを示す。図２２（Ｃ）は、人差し指付根の生体振動信号から推定した心拍間隔ＢＢＩ（〇印）と、心電図から得られた心拍間隔（実線）とを示す。図２２（Ｄ）は、人差し指付根の生体振動信号から推定した呼吸性不整脈の振幅（波形α）と、心電図から得られた呼吸性不整脈の振幅（波形β）とを示す。図２２（Ｅ）は、人差し指付根の生体振動信号から推定した位相コヒーレンス（波形α）と、心電図及び実測した呼吸パターンから算出した位相コヒーレンス（波形β）とを示す。【0127】
図２２（Ａ）に示すように、手指の生体振動信号から得られた心拍動パルス（波形α）は、実測した心電図から得られた心拍動パルス（波形β）と近似していた。図２２（Ｂ）に示すように、手指の生体振動信号から推定した呼吸パターン（波形α）は、実測した呼吸パターン（波形β）と近似していた。また、図２２（Ｃ）に示すように、手指の生体振動信号から推定した心拍間隔ＢＢＩ（〇印）は、心電図から同時計測で得られた心拍間隔ＲＲＩ（実線）と酷似しており、信号処理によって、手指の生体振動信号から心拍間隔の変動を算出することができた。図２２（Ｄ）に示すように、手指の生体振動信号から推定した呼吸性不整脈（波形α）の振幅は、心電図から得られた呼吸性不整脈の振幅（波形β）とほぼ一致していた。さらに、図２２（Ｅ）に示すように、手指の生体振動信号のみから算出した位相コヒーレンス（波形α）は、心電図及び実測した呼吸パターンから算出した位相コヒーレンス（波形β）に近似していた。これらのことから、手指の生体振動信号を測定するだけで、心拍間隔の変動（呼吸性不整脈）、呼吸パターン及び位相コヒーレンスを算出することが可能であった。【0128】
図２３は、「座位姿勢」の被験者の人差し指付根で検出した生体振動信号から得られた各種信号αと心電図又は実測の呼吸パターンから得られた各種信号βとを比較した図である。図２３（Ａ）から図２３（Ｅ）に示す信号は、図２２（Ａ）から図２２（Ｅ）に示すものと同様である。【0129】
図２３（Ａ）に示すように、手指の生体振動信号から得られた心拍動パルス（波形α）は、実測した心電図から得られた心拍動パルス（波形β）と近似していた。図２３（Ｂ）に示すように、手指の生体振動信号から推定した呼吸パターン（波形α）は、若干の位相のずれはあるものの実測した呼吸パターン（波形β）と近似していた。また、図２３（Ｃ）に示すように、手指の生体振動信号から推定した心拍間隔ＢＢＩ（〇印）は、心電図から同時計測で得られた心拍間隔ＲＲＩ（実線）と酷似していた。図２３（Ｄ）に示すように、手指の生体振動信号から推定した呼吸性不整脈の振幅（波形α）は、心電図から得られた呼吸性不整脈の振幅（波形β）とおおむね一致していた。さらに、図２３（Ｅ）に示すように、手指の生体振動信号のみから算出した位相コヒーレンス（波形α）は、心電図及び実測した呼吸パターンから算出した位相コヒーレンス（波形β）と類似の変動を有する。【0130】
図２４は、「立位姿勢」の被験者の人差し指付根で検出した生体振動信号から得られた各種信号αと心電図又は実測の呼吸パターンから得られた各種信号βとを比較した図である。図２４（Ａ）から図２４（Ｅ）に示す信号の種類は、図２２（Ａ）から図２２（Ｅ）に示すものと同様である。【0131】
図２４（Ａ）に示すように、手指の生体振動信号から得られた心拍動パルス（波形α）は、実測した心電図から得られた心拍動パルス（波形β）と近似していた。図２４（Ｂ）に示すように、手指の生体振動信号から推定した呼吸パターン（波形α）は、若干の位相のずれはあるものの実測した呼吸パターン（波形β）と近似していた。また、図２４（Ｃ）に示すように、手指の生体振動信号から推定した心拍間隔ＢＢＩ（〇印）は、心電図から同時計測で得られた心拍間隔ＲＲＩ（実線）と酷似していた。図２４（Ｄ）に示すように、手指の生体振動信号から推定した呼吸性不整脈の振幅（波形α）は、心電図から得られた呼吸性不整脈の振幅（波形β）とおおむね一致していた。さらに、図２４（Ｅ）に示すように、手指の生体振動信号のみから算出した位相コヒーレンス（波形α）は、心電図及び実測した呼吸パターンから算出した位相コヒーレンス（波形β）と類似の変動を有する。【0132】
図２５（Ａ）は、「仰臥位姿勢」の被験者の人差し指付根で検出した生体振動信号（図２２）から推定した心拍間隔ＢＢＩ（上段）とそのパワースペクトル波形（下段）である。図２５（Ｂ）は、「座位姿勢」の被験者の人差し指付根で検出した生体振動信号（図２３）から推定した心拍間隔ＢＢＩ（上段）とそのパワースペクトル波形（下段）である。図２５（Ｃ）は、「立位姿勢」の被験者の人差し指付根で検出した生体振動信号（図２４）から推定した心拍間隔ＢＢＩ（上段）とそのパワースペクトル波形（下段）である。図２５（Ａ）乃至図２５（Ｃ）の下段に示されるように、仰臥位姿勢、座位姿勢、立位姿勢の順に高周波成分（ＨＦ成分）が減弱していることがわかる。【0133】
図２６は、仰臥位姿勢、座位姿勢、立位姿勢の被験者の人差し指付根で検出した生体振動信号から求めた位相コヒーレンスλと心拍間隔のパワースペクトルとの関係を示す図である。図２６（Ａ）は、仰臥位姿勢、座位姿勢、立位姿勢における位相コヒーレンスλと、図２５下段のパワースペクトル波形におけるＨＦ成分のパワー（ｍｓ2）との関係を示す。ＨＦ成分は周波数０．１５Ｈｚ〜０．４Ｈｚまでを積分して求めたものである。図２６（Ｂ）は、仰臥位姿勢、座位姿勢、立位姿勢における位相コヒーレンスλと正規化されたＨＦ成分（ｎＨＦ）のパワー（％）との関係を示す。ｎＨＦはＨＦ成分のパワーを平均のＢＢＩで除したものである。図２６（Ｃ）は、仰臥位姿勢、座位姿勢、立位姿勢における位相コヒーレンスλと低周波成分（ＬＦ成分）のパワー／高周波成分（ＨＦ成分）のパワーの比との関係を示す【0134】
一般に、高周波成分（ＨＦ成分）は呼吸によって生ずる副交感神経活動によって影響を受け、低周波成分（ＬＦ成分）は交感神経と副交感神経活動によって影響を受ける。このため、各姿勢におけるＨＦ成分のパワーは、副交感神経機能の指標となり、各姿勢におけるＬＦ成分のパワー／ＨＦ成分のパワーは、交感神経機能の指標となる。【0135】
図２６（Ａ）及び図２６（Ｂ）に示されるように、位相コヒーレンスλは、ＨＦ成分のパワーに対して正の相関を示すことが分かった。一方、図２６（Ｃ）に示されるように、位相コヒーレンスλは、ＬＦ成分のパワー／ＨＦ成分のパワーに対して負の相関を示すことが分かった。ＨＦやｎＨＦは心臓の副交感神経活動を反映すると言われている。また，ＬＦ／ＨＦは交感神経活動を反映すると言われている。したがって，位相コヒーレンスλは心臓自律神経活動のバランスを表しており，０であれば交感神経活動の亢進を，１であれば副交感神経活動の亢進を間接的に示すことになる。図２６に示すように、仰臥位姿勢では、位相コヒーレンスλが大きく、ＨＦ成分及びｎＨＦ成分が大きく、ＬＦ／ＨＦは小さくなっており、座位姿勢では仰臥位姿勢に比べて位相コヒーレンスλが小さくなり、ＨＦ成分及びｎＨＦ成分も小さく、ＬＦ／ＨＦは大きくなり、立位姿勢では位相コヒーレンスλ、ＨＦ成分及びｎＨＦ成分がさらに小さくなり、ＬＦ／ＨＦはさらに大きくなっている。このように、手指の生体振動信号から得られた生体情報は、自律神経の指標として利用可能である。【0136】
［実施例６］本実施例では、腕輪型の生体振動信号測定装置を座位姿勢の被験者（ヒト）の上腕部（肘と肩の間）に４０μｍ厚の薄型ＰＶＤＦフィルムの両面のほぼ全面に１μｍ厚のカーボン電極を被着し、さらにその両面を２５μｍ厚のＰＥＴフィルムでラミネート加工することによりＰＶＤＦフィルムを保護し、さらに電磁シールドを目的として、その両面をアルミニウム薄膜（アルミニウム薄膜とＰＥＴの重ね膜）で覆って最終的なセンサ（幅８ｃｍ，長さ１２ｃｍ、厚さ０．５ｍｍ）を構成した。本シートセンサを上腕部の半円周分を巻き付け装着して、上腕の生体振動信号（脈波）を測定し、上腕の生体振動信号（脈波）から心拍に関する情報と呼吸に関する情報を分離抽出し、位相コヒーレンスを算出した。上腕部に振動センサを内側に実装したバンドを巻き付けて装着した。振動センサは、厚さ数ミクロン程度の薄い有機薄膜シート状フィルムを介して間接的に皮膚に接触し、振動を検出した。また、比較のため、被験者の胸部に心電図用電極を貼付し、単極導出により心電図（ＥＣＧ）を同時に計測した。
【0137】
図２７は、上腕部で検出した生体振動信号から得られた各種信号（実線、×）と心電図から得られた各種信号（点線、〇）とを比較した図である。図２７（Ａ）は、上腕部の生体振動信号から得られた心拍動パルス（実線）と、心電図から得られた心拍動パルス（点線）を示した。図２７（Ｂ）は、上腕部の生体振動信号のピークから推定した心拍間隔ＢＢＩ（×）と、心電図のピークから得られた心拍間隔ＲＲＩ（〇）とを示す。図２７（Ｃ）は、上腕部の生体振動信号のピーク値間隔の変動から推定した呼吸パターン（実線）と、心電図のピーク値間隔（心拍間隔）の変動から推定した呼吸パターン（点線）とを示す。図２７（Ｄ）は、上腕部の生体振動信号から推定した呼吸性不整脈の振幅（実線）と、心電図から得られた呼吸性不整脈の振幅（点線）とを示す。図２７（Ｅ）は、上腕部の生体振動信号から推定した位相コヒーレンス（実線）と、心電図から算出した位相コヒーレンス（点線）とを示す。【0138】
図２７（Ａ）に示すように、上腕部の生体振動信号から得られた心拍動パルス（実線）は、実測した心電図から得られた心拍動パルス（点線）と近似していた。図２７（Ｂ）に示すように、上腕部の生体振動信号から推定した心拍間隔ＢＢＩ（×）は、心電図から同時計測で得られた心拍間隔ＲＲＩ（〇）と酷似しており、信号処理によって、上腕部の生体振動信号から心拍間隔の変動を算出することができた。また、図２７（Ｃ）に示すように、上腕部の生体振動信号から推定した呼吸パターン（実線）と、心電図から推定した呼吸パターン（点線）とでは波形が若干異なっており、同じアルゴリズムによる呼吸パターンの算出は本実施例では適切ではなかった。図２７（Ｄ）に示すように、上腕部の生体振動信号から推定した呼吸性不整脈（実線）の振幅は、心電図から得られた呼吸性不整脈の振幅（点線）とほぼ一致していた。さらに、図２７（Ｅ）に示すように、上腕部の生体振動信号のみから算出した位相コヒーレンス（実線）と心電図から算出した位相コヒーレンス（点線）とは、部分的に大きなずれがあるものの、相関している部分もあり、上腕の生体振動信号から傾向を把握することは可能である。これらのことから、上腕の生体振動信号を測定するだけで、少なくとも心拍間隔の変動（呼吸性不整脈）、呼吸パターン及び位相コヒーレンスを算出することが可能であった。【0139】
［実施例７］本実施例では、実施例６と同じように、矩形の薄型センサ（幅５ｃｍ、長さ５ｃｍ程度）を製作し、それをゴム系のバンドに内蔵させたヘッドバンド型の生体振動信号測定装置を座位姿勢の被験者（ヒト）の頭部（こめかみ）に装着して、こめかみの生体振動信号を測定し、こめかみの生体振動信号から心拍に関する情報と呼吸に関する情報を分離抽出し、位相コヒーレンスを算出した。頭部に振動センサを内側に実装したヘッドバンドを巻き付けて装着した。振動センサは、ゴムを介して間接的にこめかみの皮膚に接触し、振動を検出した。また、比較のため、被験者の胸部に心電図用電極を貼付し、単極導出により心電図（ＥＣＧ）を同時に計測した。
【0140】
図２８は、こめかみで検出した生体振動信号から得られた各種信号（実線、×）と心電図から得られた各種信号（点線、〇）とを比較した図である。図２８（Ａ）は、こめかみの生体振動信号から得られた心拍動パルス（実線）と、心電図から得られた心拍動パルス（点線）を示した。図２８（Ｂ）は、こめかみの生体振動信号のピークから推定した心拍間隔ＢＢＩ（×）と、心電図のピークから得られた心拍間隔ＲＲＩ（〇）とを示す。図２８（Ｃ）は、こめかみの生体振動信号のピーク値間隔の変動から推定した呼吸パターン（実線）と、心電図のピーク値間隔の変動から推定した呼吸パターン（点線）とを示す。図２８（Ｄ）は、こめかみの生体振動信号から推定した呼吸性不整脈の振幅（実線）と、心電図から得られた呼吸性不整脈の振幅（点線）とを示す。図２８（Ｅ）は、こめかみの生体振動信号から推定した位相コヒーレンス（実線）と、心電図から算出した位相コヒーレンス（点線）とを示す。【0141】
図２８（Ａ）に示すように、こめかみの生体振動信号から得られた心拍動パルス（実線）は、実測した心電図から得られた心拍動パルス（点線）と近似していた。図２８（Ｂ）に示すように、こめかみの生体振動信号から推定した心拍間隔ＢＢＩ（×）は、心電図から同時計測で得られた心拍間隔ＲＲＩ（〇）と酷似しており、信号処理によって、こめかみの生体振動信号から心拍間隔の変動を算出することができた。また、図２８（Ｃ）に示すように、こめかみの生体振動信号から推定した呼吸パターン（実線）と、心電図から推定した呼吸パターン（点線）とでは波形が若干異なっており、同じアルゴリズムによる呼吸パターンの算出は本実施例では適切ではなかった。図２８（Ｄ）に示すように、こめかみの生体振動信号から推定した呼吸性不整脈（実線）の振幅は、心電図から得られた呼吸性不整脈の振幅（点線）とほぼ一致していた。さらに、図２８（Ｅ）に示すように、こめかみの生体振動信号のみから算出した位相コヒーレンス（実線）と心電図から算出した位相コヒーレンス（点線）とは、部分的にずれがあるものの、相関している部分もあり、こめかみの生体振動信号から傾向を把握することは可能である。これらのことから、こめかみの生体振動信号を測定するだけで、少なくとも心拍間隔の変動（呼吸性不整脈）、呼吸パターン及び位相コヒーレンスを算出することが可能であった。【0142】
本発明の一形態の指輪型、腕輪型又はバンド型の生体振動信号測定装置によれば、心拍間隔の変動の瞬時位相と呼吸パターンの瞬時位相との瞬時位相差の位相コヒーレンスが、δ波の振幅と相関関係があることを利用して睡眠状態を測定することができる。位相コヒーレンスは、呼吸周波数の影響を受けにくいため、より正確に睡眠状態を測定することができる。また、本発明の一形態の指輪型、腕輪型又はバンド型の生体振動信号測定装置によれば、被験者に負担が少なく、ストレスを感じにくく、より簡便な睡眠状態測定装置を提供できる。位相コヒーレンスは、リアルタイムに測定することができ、ノンレム睡眠時間の発生時点やノンレム―レム睡眠のリズム周期を測定することができ、睡眠状態を的確に把握することができる。また、同時に呼吸パターンを測定又は抽出しているので、睡眠時の無呼吸状態を検出することも可能である。さらに、被験者の呼吸性不整脈の大きさも測定できるので、睡眠状態について補助的な判断指標とすることも可能である。【0143】
さらに、本発明の一形態の指輪型、腕輪型又はバンド型の生体振動信号測定装置を備えた睡眠状態測定装置を構成する場合、睡眠状態測定装置は、算出した位相コヒーレンスに基づいて、睡眠状態等を判定する判定機能を有していてもよい。判定機能としては、例えば、算出した位相コヒーレンスを閾値と比較して、閾値よりも大きい場合は深い睡眠であると判定し、小さい場合は浅い睡眠と判定してもよいし、位相コヒーレンスが閾値よりも大きい値だった時間で睡眠の品質を評価してもよいし、位相コヒーレンスの変動の周期によって睡眠の品質を評価してもよい。閾値は、予め定めた数値であってもよいし、計測対象の過去に算出した位相コヒーレンスの数値から特定してもよいし、複数の数値を設定し、段階的に睡眠の質を評価してもよい。さらに、生体振動信号測定装置が備える情報処理手段が判定機能を有し、睡眠中の呼吸状態を判定することもできる。例えば、呼吸パターンの取得方法にもよるが、いずれの取得方法であっても、中枢性睡眠時無呼吸（脳の呼吸中枢の異常等により呼吸運動が停止して起こる無呼吸）については、呼吸運動が停止するので呼吸パターンが検出できなくなるので、無呼吸状態であることを判定できる。さらに、呼吸による空気の流れを計測して呼吸パターンを取得した場合は、呼吸の気流が停止して呼吸パターンが検出されなくなるので、中枢性睡眠時無呼吸だけではなく、閉塞性睡眠時無呼吸（呼吸運動はあるが気道の閉塞等による無呼吸）も判定できる。生体振動信号測定装置の表示出力手段は、判定機能で判定した睡眠状態等の結果を出力してもよい。【0144】
さらに、位相コヒーレンスは、睡眠状態だけではなく、心理ストレスの評価にも利用することができるため、本発明の位相コヒーレンス算出機能を備えた生体振動信号測定装置は、ストレス状態測定装置を構成することもできる。【0145】
本発明の生体振動信号測定装置は、さまざまな装飾品電子機器等に組み込んで使用することができる。例えば、腕輪、腕時計、指輪に振動を計測するセンサを組み込み、使用者のストレス状態を計測してもよい。また、生体振動信号測定装置で取得した生体振動信号や、得られた生体情報を通信手段を介して、携帯端末、コンピュータ等に位相コヒーレンス算出装置を組み込み、日常生活の種々の場面において心理ストレスを評価したり、車、電車、飛行機等に送信し、運転手の状態を管理することもできる。本発明の生体振動信号測定装置は、円柱形やヒト動物の体の一部に巻き付けたり、貼り付けたりすることにより、脈波等の生体信号を取得することができる。本装置に、Ｂｌｕｅｔｏｏｔｈ（登録商標）やＺｉｇｂｅｅ（登録商標）等の通信機能とバッテリー機能を組み込むことができる。この結果、本装置をトンネル工事や地下等で作業する作業員の体の一部に貼り付けることにより、作業員から作業中のリアルタイムで脈波等の生体振動信号を取得することによる作業員の健康管理にも活用することができる。また、女性の下着ブラジャーの一部に本装置を組み込んで、女性の心拍・呼吸等の生体信号を取得し、女性の健康管理に活用することもできるし、妊婦の腹帯の中に本センサを挿入して、妊婦の健康管理や、胎児の心拍や体動等の振動信号等の生体情報を取得することもできる。【0146】
本発明の生体振動信号測定装置は、ヒトからの生体情報取得に限らず、モルモット等の実験動物、犬、猫、兎、ハムスター等のペット用動物、牛、馬、羊等の家畜等の動物にも適用できる。適用箇所としては、首輪、腕輪、足環、耳輪等、上記動物の体の一部に装着して使うことが出来る。【符号の説明】
【0147】
１ 睡眠状態測定装置２ 情報取得部
３ 情報処理部
４ 操作部
５ 出力部
６ 記憶部