(57)【要約】      （修正有）【課題】非拘束、リアルタイムかつ連続的に血圧を推定することのできる血圧推定装置を提供する。【解決手段】血圧推定装置１００は、人の背部から採取される生体信号を解析し、心室充満期から等容性収縮期の血流変動によって生じる生体内の振動（生体内振動）を捉え、生体内振動の様子を示す指標のゆらぎの様子（ゆらぎに関する指標（ゆらぎ指標））を捉える。ゆらぎ指標は、人の状態、特に、血圧に関する情報と関連性がある。乗物のシート、事務用や家具用などの各種椅子、あるいは、ベッドなどの寝具等に取り付けられた、非拘束で生体信号を測定できる生体信号測定装置に、人が着席あるいは仰臥して背部を当接すると、容易に血圧の推定、特に、正常域血圧（拡張期血圧が９０ｍｍＨｇ未満、収縮期血圧が１４０ｍｍＨｇ未満）を示す領域に含まれるか高血圧を示す領域に含まれるかを推定できる。【選択図】図２
2019122502.tif 000002
【特許請求の範囲】
【請求項１】
人の背部に当接され、前記背部の体表面を介して伝播する生体信号を非拘束で捉える生体信号測定装置から前記生体信号を受信し、受信した前記生体信号を解析する生体信号処理手段を有しており、
前記生体信号処理手段が、
前記生体信号の時系列波形を所定の周波数帯域でフィルタリングし、心周期を顕在化させたフィルタ処理波形を得るフィルタリング手段と、
前記生体信号の測定と同時に測定した心電計から得られる心電図波形データに前記フィルタ処理波形を照合し、前記フィルタ処理波形における心室充満期から等容性収縮期における時系列波形の範囲の波形成分を特定し、心室充満期から等容性収縮期における血流変動に起因して生じる振動に関する指標を求める血流変動指標算出手段と、
前記血流変動に起因して生じる振動に関する指標の時系列波形を求め、その時系列波形の変動の仕方を示すゆらぎに関する指標を求めるゆらぎ指標算出手段と、
前記ゆらぎに関する指標を、予め記憶部に記憶された前記ゆらぎに関する指標と血圧との関係を示す血圧推定用相関データを用いて、前記人の血圧を推定する推定手段と
を有することを特徴とする血圧推定装置。
【請求項２】
前記生体信号処理手段は、連続的に受信した前記生体信号を用い、前記人の血圧を連続的に推定可能である請求項１記載の血圧推定装置。
【請求項３】
前記ゆらぎ指標算出手段が、
前記血流変動に起因して生じる振動に関する指標の時系列波形を求めると共に、その時系列波形を周波数解析し、パワースペクトルと周波数の両対数グラフに周波数解析結果であるゆらぎ解析図を得るゆらぎ解析手段と、
前記ゆらぎ解析図の回帰直線の傾きを求め、前記ゆらぎ解析図の回帰直線の傾きを前記ゆらぎに関する指標とするゆらぎ解析図傾き算出手段と
を有する請求項１又は２記載の血圧推定装置。
【請求項４】
前記ゆらぎ解析図傾き算出手段は、
前記ゆらぎ解析図に関する回帰直線の傾きとして、ＶＬＦからＬＦの範囲に属する所定の周波数帯の回帰直線の傾きを求め、その傾きを前記ゆらぎに関する指標とする手段である請求項３記載の血圧推定装置。
【請求項５】
前記血圧推定用相関データは、前記ＶＬＦからＬＦの範囲に属する所定の周波数帯の回帰直線の傾きが、正常域血圧では−１に近い傾向を示し、高血圧では＋１に近い傾向を示す相関データとして形成されている請求項４記載の血圧推定装置。
【請求項６】
前記血流変動指標算出手段が、
前記生体信号中、前記心室充満期から等容性収縮期に対応する時系列波形の範囲において２つの波形成分を一組以上特定し、前記２つの波形成分の各全振幅を用いて散布図を作成する散布図作成手段と、
前記散布図にプロットされた点群の回帰直線の傾きを求め、前記点群の回帰直線の傾きを前記血流変動に起因して生じる振動に関する指標とする散布図傾き算出手段と
を有する請求項１〜５のいずれか１に記載の血圧推定装置。
【請求項７】
前記散布図作成手段が、前記心室充満期の心房収縮時における前記２つの波形成分の各全振幅を用いた第１の散布図と、前記第１の散布図に用いた前記２つの波形成分より後のタイミング以降であって等容性収縮期に移行する房室弁閉鎖時付近のタイミングに対応する前記２つの波形成分の各全振幅を用いた第２の散布図とを作成する手段であり、
前記ゆらぎ指標算出手段が前記第１の散布図を用いて求めたゆらぎ指標から、前記推定手段が拡張期血圧を推定し、
前記ゆらぎ指標算出手段が前記第２の散布図を用いて求めたゆらぎ指標から、前記推定手段が収縮期血圧を推定する請求項６記載の血圧推定装置。
【請求項８】
人の背部に当接され、前記背部の体表面を介して伝播する生体信号を非拘束で捉える生体信号測定装置から前記生体信号を受信し、
前記生体信号の時系列波形を所定の周波数帯域でフィルタリングし、心周期を顕在化させたフィルタ処理波形を得て、
前記生体信号の測定と同時に測定した心電計から得られる心電図波形データに前記フィルタ処理波形を照合し、前記フィルタ処理波形における心室充満期から等容性収縮期における時系列波形の範囲の波形成分を特定し、心室充満期から等容性収縮期における血流変動に起因して生じる振動に関する指標を求め、
前記血流変動に起因して生じる振動に関する指標の時系列波形を求めた後、その時系列波形の変動の仕方を示すゆらぎに関する指標を求め、
前記ゆらぎに関する指標を、予め記憶部に記憶された前記ゆらぎに関する指標と血圧との関係を示す血圧推定用相関データを用いて、前記人の血圧を推定することを特徴とする血圧推定方法。
【請求項９】
前記生体信号測定装置から、前記生体信号を連続的に受信し、前記人の血圧を連続的に推定する請求項８記載の血圧推定方法。
【請求項１０】
前記生体信号測定装置の生体信号検出部を、前記人の背部において、鎖骨の位置に対応する部位から剣状突起の位置に対応する部位までの範囲に設けて前記生体信号を捉える請求項８又は９記載の血圧推定方法。
【請求項１１】
人の背部に当接され、前記背部の体表面を介して伝播する生体信号を非拘束で捉える生体信号測定装置から前記生体信号を受信し、受信した前記生体信号を解析する生体信号処理手順を実行させ、コンピュータを血圧推定装置として機能させるコンピュータプログラムであって、
前記生体信号処理手順として、
前記生体信号の時系列波形を所定の周波数帯域でフィルタリングし、心周期を顕在化させたフィルタ処理波形を得る手順と、
前記生体信号の測定と同時に測定した心電計から得られる心電図波形データに前記フィルタ処理波形を照合し、前記フィルタ処理波形における心室充満期から等容性収縮期における時系列波形の範囲の波形成分を特定し、心室充満期から等容性収縮期における血流変動に起因して生じる振動に関する指標を求める手順と、
前記血流変動に起因して生じる振動に関する指標の時系列波形を求め、その時系列波形の変動の仕方を示すゆらぎに関する指標を求める手順と、
前記ゆらぎに関する指標を、予め記憶部に記憶された前記ゆらぎに関する指標と血圧との関係を示す血圧推定用相関データを用いて、前記人の血圧を推定する手順と
を前記コンピュータに実行させるコンピュータプログラム。
【請求項１２】
前記生体信号処理手順は、連続的に受信した前記生体信号を用い、前記コンピュータに、前記人の血圧の推定を連続的に実行させる請求項１１記載のコンピュータプログラム。
【請求項１３】
前記ゆらぎに関する指標を求める手順では、
前記血流変動に起因して生じる振動に関する指標の時系列波形を求めると共に、その時系列波形を周波数解析し、パワースペクトルと周波数の両対数グラフに周波数解析結果であるゆらぎ解析図を得る手順と、
前記ゆらぎ解析図の回帰直線の傾きを求め、前記ゆらぎ解析図の回帰直線の傾きを前記ゆらぎに関する指標とする手順と
を前記コンピュータに実行させる請求項１１又は１２記載のコンピュータプログラム。
【請求項１４】
前記ゆらぎ解析図に関する回帰直線の傾きとして、ＶＬＦからＬＦの範囲に属する所定の周波数帯の回帰直線の傾きを求め、その傾きを前記ゆらぎに関する指標とする請求項１３記載のコンピュータプログラム。
【請求項１５】
血流変動に起因して生じる振動に関する指標を求める手順では、
前記生体信号中、前記心室充満期から等容性収縮期に対応する時系列波形の範囲において２つの波形成分を一組以上特定し、前記２つの波形成分の各全振幅を用いて散布図を作成する手順と、
前記散布図にプロットされた点群の回帰直線の傾きを求め、前記点群の回帰直線の傾きを前記血流変動に起因して生じる振動に関する指標とする手順と
を前記コンピュータに実行させる請求項１１〜１４のいずれか１に記載のコンピュータプログラム。
【請求項１６】
前記散布図を作成する手順として、前記心室充満期の心房収縮時における前記２つの波形成分の各全振幅を用いた第１の散布図と、前記第１の散布図に用いた前記２つの波形成分より後のタイミング以降であって等容性収縮期に移行する房室弁閉鎖時付近のタイミングに対応する前記２つの波形成分の各全振幅を用いた第２の散布図とを作成する手順が実行され、
前記ゆらぎに関する指標を求める手順では、前記第１の散布図を用いて求めたゆらぎ指標から拡張期血圧を推定し、前記第２の散布図を用いて求めたゆらぎ指標から、収縮期血圧を推定する請求項１５記載のコンピュータプログラム。
【請求項１７】
人の背部に当接され、前記背部の体表面を介して伝播する生体信号を非拘束で捉える生体信号測定装置から前記生体信号を受信し、受信した前記生体信号を解析する生体信号処理手順を実行させ、コンピュータを血圧推定装置として機能させる請求項１１〜１６のいずれか１に記載のコンピュータプログラムが記録されたコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、人の背部から得られる生体信号を用いて、血圧を推定する技術に関する。
【背景技術】
【０００２】
本発明者らは、特許文献１〜４等において、人の上体の背部の体表面に生じる振動を非拘束で捉え、その振動を解析して人の状態を推定する技術を提案している。人の上体の背部の体表面に生じる振動は、心臓と大動脈等の生体内の振動が伝播したものであり、心房及び心室の収縮期及び拡張期の情報や、循環の補助ポンプとなる血管壁の弾力情報及び反射波の情報を含んでいる。
【０００３】
特許文献１では、体表面を介して伝播する振動（生体信号）から抽出した１Ｈｚ近傍の背部体表脈波（Ａｏｒｔｉｃ  Ｐｕｌｓｅ  Ｗａｖｅ（ＡＰＷ））の時系列波形に所定の時間幅を適用してスライド計算を行って周波数傾きの時系列波形を求め、その変化の傾向から、例えば、振幅が増幅傾向にあるか、減衰傾向にあるかなどによって生体状態の推定を行っている。また、生体信号を周波数解析し、予め定めたＵＬＦ帯域（極低周波帯域）からＶＬＦ帯域（超低周波帯域）に属する機能調整信号、疲労受容信号及び活動調整信号に相当する各周波数のパワースペクトルを求め、各パワースペクトルの時系列変化から人の状態を判定することも開示している。
【０００４】
特許文献２〜４では、恒常性維持機能レベルを判定する手段を開示している。恒常性維持機能レベルを判定する手段は、周波数傾き時系列波形の微分波形の正負、周波数傾き時系列波形を積分した積分波形の正負、ゼロクロス法を利用した周波数傾き時系列波形とピーク検出法を利用した周波数傾き時系列波形をそれぞれ絶対値処理して得られた各周波数傾き時系列波形の絶対値等のうち、少なくとも１つ以上を用いて判定する。これらの組み合わせにより、恒常性維持機能のレベルがいずれに該当するかを求める。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１１−１６７３６２号公報
【特許文献２】ＷＯ２０１１／０４６１７８号公報
【特許文献３】特開２０１４−１１７４２５号公報
【特許文献４】特開２０１４−２２３２７１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
上記した技術によれば、生体調節機能に関する人の状態、すなわち、入眠予兆現象、切迫睡眠現象、覚低走行状態、恒常性維持機能レベル、初期疲労状態、気分の状態などを捉えることができる。
その一方、高血圧の人や心疾患等を患っている人は、健常者と比較すると、血管壁の状態、心臓の動き方等に起因して、血液の流れにも影響が生じる。従来、高血圧か否かを判断する場合、カフを腕に巻き付けるタイプの血圧計によって血圧を測定するのが一般的である。しかしながら、これらの測定には、人体を拘束しかつ安定した状態とした上で、各装置を取り付けて実施する必要がある。そのため、人体の拘束等を極力減らし、血圧をより簡易に測定できる技術の開発が望まれている。
【０００７】
本発明は上記に鑑みなされたものであり、人の背部の体表面を介して伝播する生体信号を非拘束で捉えることができる生体信号測定装置を用い、取得した生体信号を解析し、非拘束で血圧の推定を行う技術を提供することを課題とする。
【課題を解決するための手段】
【０００８】
上記課題を解決するため、本発明者は、酸素の摂取量と血圧の相関性から、静脈還流の変動及び心筋収縮力の変動に伴う血流量の変動が血圧変動を生むことに着目した。また、左心室の収縮は血液の体循環の原動力となるが、右心室は左心室を包むように存在し、左心室は胸腔内の正中面左方に深く隠れていると共に、右心室の前面にある心膜は胸骨心膜靱帯で前胸壁と結合している。そのため、前胸壁から得られる心臓の情報は限られ、特に、体循環の原動力である左心室の情報は前胸壁からは計測しがたい。左心室の情報は、前胸壁よりも、背部の体表面に伝播される振動（本明細書において、単に「振動」という場合には、音（音波）も含む）に多く含まれている。そこで、本発明者は、背部の体表面に伝播される振動（生体信号）から、血流変動の状態を知ることで、血圧を推定できると考え、本発明を完成するに至った。
【０００９】
すなわち、本発明の血圧推定装置は、人の背部に当接され、前記背部の体表面を介して伝播する生体信号を非拘束で捉える生体信号測定装置から前記生体信号を受信し、受信した前記生体信号を解析する生体信号処理手段を有しており、前記生体信号処理手段が、  前記生体信号の時系列波形を所定の周波数帯域でフィルタリングし、心周期を顕在化させたフィルタ処理波形を得るフィルタリング手段と、前記生体信号の測定と同時に測定した心電計から得られる心電図波形データに前記フィルタ処理波形を照合し、前記フィルタ処理波形における心室充満期から等容性収縮期における時系列波形の範囲の波形成分を特定し、心室充満期から等容性収縮期における血流変動に起因して生じる振動に関する指標を求める血流変動指標算出手段と、前記血流変動に起因して生じる振動に関する指標の時系列波形を求め、その時系列波形の変動の仕方を示すゆらぎに関する指標を求めるゆらぎ指標算出手段と、前記ゆらぎに関する指標を、予め記憶部に記憶された前記ゆらぎに関する指標と血圧との関係を示す血圧推定用相関データを用いて、前記人の血圧を推定する推定手段とを有することを特徴とする。
【００１０】
前記生体信号処理手段は、連続的に受信した前記生体信号を用い、前記人の血圧を連続的に推定可能である。
前記ゆらぎ指標算出手段が、前記血流変動に起因して生じる振動に関する指標の時系列波形を求めると共に、その時系列波形を周波数解析し、パワースペクトルと周波数の両対数グラフに周波数解析結果であるゆらぎ解析図を得るゆらぎ解析手段と、前記ゆらぎ解析図の回帰直線の傾きを求め、前記ゆらぎ解析図の回帰直線の傾きを前記ゆらぎに関する指標とするゆらぎ解析図傾き算出手段とを有することが好ましい。
前記ゆらぎ解析図傾き算出手段は、前記ゆらぎ解析図に関する回帰直線の傾きとして、ＶＬＦからＬＦの範囲に属する所定の周波数帯の回帰直線の傾きを求め、その傾きを前記ゆらぎに関する指標とする手段であることが好ましい。
前記血圧推定用相関データは、前記ＶＬＦからＬＦの範囲に属する所定の周波数帯の回帰直線の傾きが、正常域血圧では−１に近い傾向を示し、高血圧では＋１に近い傾向を示す相関データとして形成されていることが好ましい。
前記血流変動指標算出手段が、前記生体信号中、前記心室充満期から等容性収縮期に対応する時系列波形の範囲において２つの波形成分を一組以上特定し、前記２つの波形成分の各全振幅を用いて散布図を作成する散布図作成手段と、前記散布図にプロットされた点群の回帰直線の傾きを求め、前記点群の回帰直線の傾きを前記血流変動に起因して生じる振動に関する指標とする散布図傾き算出手段とを有することが好ましい。
【００１１】
前記散布図作成手段が、前記心室充満期の心房収縮時における前記２つの波形成分の各全振幅を用いた第１の散布図と、前記第１の散布図に用いた前記２つの波形成分より後のタイミング以降であって等容性収縮期に移行する房室弁閉鎖時付近のタイミングに対応する前記２つの波形成分の各全振幅を用いた第２の散布図とを作成する手段であり、前記ゆらぎ指標算出手段が前記第１の散布図を用いて求めたゆらぎ指標から、前記推定手段が拡張期血圧を推定し、前記ゆらぎ指標算出手段が前記第２の散布図を用いて求めたゆらぎ指標から、前記推定手段が収縮期血圧を推定する構成であることが好ましい。
【００１２】
本発明の血圧推定方法は、人の背部に当接され、前記背部の体表面を介して伝播する生体信号を非拘束で捉える生体信号測定装置から前記生体信号を受信し、前記生体信号の時系列波形を所定の周波数帯域でフィルタリングし、心周期を顕在化させたフィルタ処理波形を得て、前記生体信号の測定と同時に測定した心電計から得られる心電図波形データに前記フィルタ処理波形を照合し、前記フィルタ処理波形における心室充満期から等容性収縮期における時系列波形の範囲の波形成分を特定し、心室充満期から等容性収縮期における血流変動に起因して生じる振動に関する指標を求め、前記血流変動に起因して生じる振動に関する指標の時系列波形を求めた後、その時系列波形の変動の仕方を示すゆらぎに関する指標を求め、前記ゆらぎに関する指標を、予め記憶部に記憶された前記ゆらぎに関する指標と血圧との関係を示す血圧推定用相関データを用いて、前記人の血圧を推定することを特徴とする。
【００１３】
前記生体信号測定装置から、前記生体信号を連続的に受信し、前記人の血圧を連続的に推定することが好ましい。
前記生体信号測定装置の生体信号検出部を、前記人の背部において、鎖骨の位置に対応する部位から剣状突起の位置に対応する部位までの範囲に設けて前記生体信号を捉えることが好ましい。
【００１４】
本発明のコンピュータプログラムは、人の背部に当接され、前記背部の体表面を介して伝播する生体信号を非拘束で捉える生体信号測定装置から前記生体信号を受信し、受信した前記生体信号を解析する生体信号処理手順を実行させ、コンピュータを血圧推定装置として機能させるコンピュータプログラムであって、前記生体信号処理手順として、前記生体信号の時系列波形を所定の周波数帯域でフィルタリングし、心周期を顕在化させたフィルタ処理波形を得る手順と、前記生体信号の測定と同時に測定した心電計から得られる心電図波形データに前記フィルタ処理波形を照合し、前記フィルタ処理波形における心室充満期から等容性収縮期における時系列波形の範囲の波形成分を特定し、心室充満期から等容性収縮期における血流変動に起因して生じる振動に関する指標を求める手順と、前記血流変動に起因して生じる振動に関する指標の時系列波形を求め、その時系列波形の変動の仕方を示すゆらぎに関する指標を求める手順と、前記ゆらぎに関する指標を、予め記憶部に記憶された前記ゆらぎに関する指標と血圧との関係を示す血圧推定用相関データを用いて、前記人の血圧を推定する手順とを前記コンピュータに実行させることを特徴とする。
【００１５】
前記生体信号処理手順は、連続的に受信した前記生体信号を用い、前記コンピュータに、前記人の血圧の推定を連続的に実行させることが好ましい。
前記ゆらぎに関する指標を求める手順では、前記血流変動に起因して生じる振動に関する指標の時系列波形を求めると共に、その時系列波形を周波数解析し、パワースペクトルと周波数の両対数グラフに周波数解析結果であるゆらぎ解析図を得る手順と、前記ゆらぎ解析図の回帰直線の傾きを求め、前記ゆらぎ解析図の回帰直線の傾きを前記ゆらぎに関する指標とする手順とを前記コンピュータに実行させることが好ましい。
前記ゆらぎ解析図に関する回帰直線の傾きとして、ＶＬＦからＬＦの範囲に属する所定の周波数帯の回帰直線の傾きを求め、その傾きを前記ゆらぎに関する指標とすることが好ましい。
【００１６】
血流変動に起因して生じる振動に関する指標を求める手順では、前記生体信号中、前記心室充満期から等容性収縮期に対応する時系列波形の範囲において２つの波形成分を一組以上特定し、前記２つの波形成分の各全振幅を用いて散布図を作成する手順と、前記散布図にプロットされた点群の回帰直線の傾きを求め、前記点群の回帰直線の傾きを前記血流変動に起因して生じる振動に関する指標とする手順とを前記コンピュータに実行させることが好ましい。
前記散布図を作成する手順として、前記心室充満期の心房収縮時における前記２つの波形成分の各全振幅を用いた第１の散布図と、前記第１の散布図に用いた前記２つの波形成分より後のタイミング以降であって等容性収縮期に移行する房室弁閉鎖時付近のタイミングに対応する前記２つの波形成分の各全振幅を用いた第２の散布図とを作成する手順が実行され、前記ゆらぎに関する指標を求める手順では、前記第１の散布図を用いて求めたゆらぎ指標から拡張期血圧を推定し、前記第２の散布図を用いて求めたゆらぎ指標から、収縮期血圧を推定することが好ましい。
【００１７】
本発明のコンピュータ読み取り可能な記録媒体は、人の背部に当接され、前記背部の体表面を介して伝播する生体信号を非拘束で捉える生体信号測定装置から前記生体信号を受信し、受信した前記生体信号を解析する生体信号処理手順を実行させ、コンピュータを血圧推定装置として機能させる上記コンピュータプログラムが記録されていることを特徴とする。
【００１８】
ここで、本発明によって血圧が推定できる原理を説明する。
人の心拍出量は、身体の代謝需要に応じて規定され、心拍出量は、心拍数と一回拍出量の積に等しい。心拍数は心臓血管中枢を中心とする自律神経系支配と体液性支配で制御されている。拍出量は、室の前負荷に影響を受け、一回拍出量は心収縮力と前負荷・後負荷によって決まる。血流量の変動の様子は、ベルヌーイ定理を用いた数理モデルで説明することができる。ベルヌーイの定理は圧力エネルギー、運動エネルギー、位置エネルギーからなり、水柱の高さに置き換えて圧力水頭（pressure head）、速度水頭（velocity head）、位置水頭（elevation head）と呼ばれ、各水頭の総和が全水頭（total head）となる。摩擦・曲がり・出入り口で失われるエネルギーは、速度水頭に比例する損失水頭（head loss）であり、流量変動の支配方程式は、前負荷・後負荷の血液の水頭と流量の関数で求められる。
【００１９】
本発明者は、これらを応用し、上記の体循環の原動力である左心室の血流量の変動の仕方、血流制御の在り方を知る方法として、損失水頭の変動から知ることを考えた。心臓から背部体表面に伝わる振動（音）の変動情報は心筋の収縮力の変動による血流量の変動より作られる。血流量の変動は、心筋の拡張能、収縮能より作られる。心拍出量と全末梢血管抵抗の積は血圧となり、心拍数と心収縮力の変動が血圧変動を生む。血圧の勾配によって、血液は流れる。収縮と弛緩により拍動が生じ、拡張期と収縮期が交互にある心臓の拍動により、動脈圧は変化する。つまり、この動脈圧の変動を捉えるためのパラメータが、左心室の振動（音）変動情報の中に隠されているということになる。そして心周期のどの時相の変動を見ればよいのか、それが重要なポイントになる。心室周期は四つの期に分けて考えられる。心室充満期、等容性収縮期、拍出期、等容性弛緩期である。心室充満期は、初期の急速充満期、ゆるやかな充満期、そして最後の心房収縮による充満期がある。心房収縮は心電図Ｐ波に始まる。圧格差に基づく自然の心室充満に対して、左心房はブースター機能となり、併せて左室の前負荷となる。この静脈還流に基づく前負荷により、左室の拡張末期圧が上昇する。心房収縮時以降の血流変動は左心拡張期圧変動と心拍数の変動に相関することは知られている。ここに拡張期血圧を知る手掛かりは心室充満期、特に心房収縮による充満期にあると考えられる。
【００２０】
血流変動の定常性は平均値、相関関数（ＰＳＤ）などの統計量が測定時刻によらず一定であるという性質にある。平均値は心拍数であり、心拍数の変動からなる時系列波形のＰＳＤによって自律神経系活動が定量され、心拍ゆらぎと呼ばれる。心拍ゆらぎは、１／ｆβゆらぎと呼ばれる数値で表される。心拍数の時系列データが定常であれば、ＰＳＤの対数表示が周波数の対数表示に反比例する。スペクトル指数β＝１、１／ｆゆらぎと呼ばれるものであり、βが１より大きい場合、生理学的意味においては非定常と定義される。言い換えると、一般的に定常状態と呼ばれる安静時では、心拍ゆらぎは、全体として１／ｆゆらぎとなり、フラクタル・ブラウン運動を示す。時折の血圧上昇や降圧剤は、ゆらぎを高止まりにする。高止まりとなる周期成分が任意の時間幅で混入すると、心拍数と心収縮力に変動が生じ、血圧が変化する。ストレスは血圧を変化させ、血流変動をもたらす。血流変動は、各周期成分のＰＳＤに大変動と停滞をもたらし、時刻に依存する時系列信号を生み出す。大変動と停滞は、血流変動のゆらぎに非定常性をもたらす。
【００２１】
生体に対する外乱が存在しない状態、すなわちストレスフリーの状態では、心拍数や血流量のゆらぎは１／ｆβ型となり、正常血圧者は、β＝１、１／ｆβは−１になる。高血圧患者は、自律神経系支配と体液性支配が安静状態でも現れ、フラクタル・ブラウン運動に交感神経支配による周期成分が含まれるようになる。すなわち、高血圧患者は、心拍数や血流量のゆらぎを示す１／ｆβは全体としておおよそ−１となるが、そのうち交感神経支配の状態を示すＬＦ帯域を含む領域のゆらぎに高止まりが現れ、その領域において＋成分が現れる。 
【００２２】
また、降圧剤投与による人為的な血圧低下は、反射性交感神経緊張による心収縮性の増大により、血圧変動は抑制される。また生体の反射性調節が働かない場合には、心室と動脈系の力学的相互関係の最適結合に従って決定される血圧が維持される。自律神経系が作動する状況では、生体はあらかじめ設定した血圧を維持することを優先して、血液循環を調節する。すなわち、心血管系の整合は自律神経系により大きく修飾されるが、心臓および血管の交感神経応答性の違いにより変化し、その整合状況が維持されるわけではない。人為的な降圧剤投与は、心血管系の整合に対して、強い影響力を発揮し、降圧剤による制御は、血流変動のゆらぎが１／ｆβの一部に線形部分、すなわち、＋成分が現れるように推移していくことが考えられる。 
【００２３】
また、左心拡張期圧の支配方程式は、心拍数の自律神経支配と体液性支配に相関するパラメータであるＶＬＦからＬＦ帯域までのゆらぎの積分強度が関与している。血流変動並びに正常化因子は、心臓に帰ってきた血流によって引き起こされる物理現象から心機能を判断でき、心室の前負荷が鍵となり、時相は心房収縮から心室収縮に至る境界にあるとする。そこで、解析対象の時相は、等容性収縮期の直前直後となる。 
【００２４】
交感神経は、β受容体刺激により促進性Ｇ蛋白を介して心拍数増加、心収縮力増強などを生じさせる。これら交感神経緊張に伴い、臨床的に観察される現象は、約１５秒でピークに達するといわれている。一方、交感神経刺激に対する効果の発現は、心拍数に対して、一拍ごとの動的な調整はできないことも知られている。ここに安静状態での心拍変動と心収縮力のゆらぎを観察することで、β受容体刺激による心拍数増加、心収縮力増強による血圧変動を捉えられる可能性があると考えられる。血圧変動は、血流変動のゆらぎと相関し、血流変動は背部体表音の周波数帯域上に現れる。
【００２５】
ところで、外側の心膜と内側の心膜の間にある心膜腔には数ｍｌの液体がたまっている。この心膜腔液は心臓のポンプ作用に影響を与える。成人の心臓の重さは２００〜３００ｇで心音の周波数は２５〜４５Ｈｚ、ないし５０Ｈｚといわれている。心音の周波数が心臓を浮動させる構造体の固有振動数に一致するとして心臓浮動構造体の動的ばね定数を計算すると、約２．５ｋｇ／ｍｍ前後となる。この計算結果を正として考えると、心臓浮動構造体は７０Ｈｚ以上の高周波に対して除振機能を有し、低周波帯域は剛体として考えられ、固有振動数以下の２５Ｈｚ前後の低周波数帯域の振動・音響波形は線形・弱非線形情報となる。 従って、共振帯域を除く２５Ｈｚ以下の振動波形を計測対象とする。
【００２６】
体循環の支配因子である自律神経支配とホルモンによる体液性支配は、心拍変動の周波数帯域のゆらぎの積分強度で識別され、０．１５〜０．４ＨｚのＨＦ帯域と０．０４〜０．１５ＨｚのＬＦ帯域のゆらぎの積分強度は、副交感間神経および交感神経の各活動状態を表し、０．００３〜０．０４ＨｚのＶＬＦ帯域の積分強度は体液性支配により発現するといわれている。そして、経験的に得られたＶＬＦの亢進は、視床下部室傍核のコルチコトロピン分泌ホルモン（ＣＲＨ）により活動を亢進される副腎皮質刺激ホルモン（ＡＣＴＨ）と副交感神経およびエピネフリンの亢進、そして交感神経とアンジオテンシンの抑制が条件となり、心拍数増大に寄与する。またＣＲＨ−ＡＣＴＨ系と交感神経がともに賦活するのは耐ストレス反応である。なお、夜間はメラトニン分泌による交感神経系の強力な抑制により、ＣＲＨ−ＡＣＴＨ系の亢進と副交感神経の賦活化にもかかわらず心拍数および血圧は低値に維持されることも知られている。
【００２７】
血流変動は、一回拍出量と心拍数の変動で表され、計測状態を安静状態かつ短時間で規定すると一回拍出量の変動は安定し、心拍数の変動が血流変動の支配因子として高い寄与率を持つ。心拍数の変動は、ＶＬＦとＬＦが支配因子となるため、計測のための時間幅をＶＬＦから求めると、３６０秒間となり、１／３６０は０．００３ＨｚでＶＬＦの帯域の周波数となる。左心室の拡張能と収縮能が血圧の変動に対して依存性を持つと考えられ、血流変動は、質量の変動と速度の変動が支配因子で、質量の変動は血液量の変動として現れ、速度の変動は心拍変動として現れる。力積となる血圧変動は、質量の変動と速度の変動が支配因子で、結果として２つの因子が絡み合って生じるのが時系列波形の振幅変動（ゆらぎ）となる。ここに、弁を管路とした数理モデルを検討する。 
【００２８】
図２４（ａ），（ｂ）は、心室充満期から等容性収縮期で弁が解放していくまでの心房と心室の容積変化と血液の動きを示したものである。ｘ は心房・僧帽弁・心室を移動する血液の線座標を示す。僧帽弁の入り口を０とし、出口をｘ＝ｌとする。添え字は入り口で１、出口で２とする。僧帽弁の断面積をＡ、単位体積の血液の重量をγ、重力加速度をｇ、流速をｖ、血液の水頭をｈ、圧力をＰとする。 
【００２９】
僧帽弁入り口ｘ＝０のときの血液の圧力はＰ１＝γｈ１ 
となり、 
時間Δｔの間に流量Ｑ の血液を押し込むのに要する仕事ΔＷ１は、
ΔＷ１＝Ｐ１ＱΔｔ＝γｈ１ＱΔｔ・・・（１）
となる。 
【００３０】
同様に僧帽弁出口端でなされる仕事は、
ΔＷ２＝γｈ２ＱΔｔ ・・・（２）
となる。 
【００３１】
流体の圧力損失は流量Ｑの関数であるからｆ（Ｑ）とする。この損失のために僧帽弁内で、時間Δｔの間に費やされるエネルギーΔＬは、
ΔＬ＝γｆ（Ｑ）ＱΔｔ・・・（３）
となる。 
【００３２】
僧帽弁内の移動量ｄｘに対応する運動エネルギーは、
ｄＫ＝（Ｖ２／２ｇ）×Ａγｄｘ ・・・（４）
で表される。
Ａは血液の面積で、Ａγは血液の重量である。 
【００３３】
従って、
ｄｋ＝（１／２ｇＡ）×γＱ２ｄｘ・・・（５）
となる。
【００３４】
ゆえに血液の持つ運動エネルギーＫは、
Ｋ＝γＱ２ｍ／２・・・（６）
で表される。 
但し、ｍは僧帽弁の性状により定まる定数で、弁定数と呼ぶ。 
【００３５】
微小時間Δｔ間の僧帽弁内の血液の運動エネルギーの変化ΔＫを求める。
ΔＫ＝ｍγＱ（ｄＱ／Δｔ）・・・（７）
となる。
【００３６】
ΔＫ＝ΔＷ１−ΔＷ２−ΔＬ・・・（８）
なる関係より、
ｍ（ｄＱ／ｄｔ）＝ｈ１−ｈ２−ｆ（Ｑ）・・・（９）
が求められる。 
【００３７】
図２４（ａ）に示すｈ１は前負荷である心房収縮による圧力水頭、ｈ２は後負荷で心室拡張期圧による圧力水頭である。
心室拡張気圧の変化は、流量Ｑの時間的変化が与えられれば、求めることができる。 
（９）式は整理すると、 
ｈ１−ｈ２＝ｍ（ｄＱ／ｄｔ）＋ｆ（Ｑ）・・・（１０）
で表され、 
中心静脈圧、心房圧は、心房収縮時直前では、ほぼ一定と考え、かつ水頭損失はＩＶが小さいため、ｆ（Ｑ）≒０とすると、 
ｈ１−ｈ２＝ｍ（ｄＱ／ｄｔ）・・・（１１）
となる。 
【００３８】
流量はＱ０から０に急速に、かつ一様に変化するとし、僧房弁閉鎖に要する時間はＴｍで表すと、心室拡張圧の上昇Δ（ｈ１−ｈ２）は、 
Δ（ｈ１−ｈ２）＝ｍＱ０／Ｔｍ＝ｌ／ｇＡ×Ｑ０／Ｔｍ＝ｌｖ０／ｇＴｍ．．．（１２）となる。 
【００３９】
ここに、水頭の上昇は急速充満期の流速が大きいほど、僧房弁を閉じる時間が短いほど、水撃作用は大きくなり、大きなＩ音が生じることがわかる。ｍは僧帽弁・大動脈弁の導管としての機能によって決まる定数である。ｈ１は前負荷を表し、心房収縮の静脈還流圧と心房収縮能で決定される心房圧である。後負荷は心室内圧である。ここにｈ１−ｈ２は、心室まで帰ってくる血液を押し出すのに必要な心房の収縮能に相関し、上腕の拡張期血圧とも相関性がある。また左房内圧と左室内圧の差がｈ１−ｈ２で、この差が等しくなる時がＲ波の出現ポイントで僧房弁の閉鎖の時である。 
ｆ（Ｑ）＝ｍ（ｄＱ／ｄｔ）・・・（１３）
Δｆ（Ｑ）∝ｄＰｍａｘ／ｄｔ・・・（１４）
となり、収縮期血圧も同様にしてＲ波近傍の血液量の変動と相関性が取れる可能性が示唆される。
【００４０】
ここで、高血圧患者と正常被験者の差について検討してみる。心室は拡張しきると、静脈圧によって血液がさらに心室内に流入する。充満期の最後の時相で、左房収縮でさらなる血液が左室内に押し込まれる。２０〜３０歳代の成人で安静状態では、心房収縮による血液の流入量は全体の１０〜２０％に過ぎないが、加齢とともにその比率は上昇し、４６％前後まで上昇することは知られている。左房収縮による血液充満は、高血圧患者では、その比率が上昇すると考えられる。なぜなら、高血圧患者は２０〜３０歳代の若年成人時に比べて比較的心拍数は増加傾向にある。心拍数が増加すると受動的な血液充満のための時間が短くなる。この充満期の最後の時相における流量は、左房収縮による左室充満が拡張期血圧を決定する際の重要なファクターになるとし、このファクターは、交感神経刺激により変動する。この変動の仕方は、神経伝達物質であるノルアドレナリン、アドレナリンの分泌の影響を受け、ノルアドレナリンは心筋細胞膜のβ１受容体と結合し、アドレナリンも心筋β１受容体に作用する。アドレナリンとノルアドレナリンは、カテコールアミンと呼ばれ、β１受容体を活性化させ、心拍数の増加、収縮力の増加を引き起こす。
【００４１】
次に、式（１０）中のｄＱ／ｄｔを代数方程式に変換するために、血流量の変動を捉える図式解法としてのローレンツプロット法の利用について説明する。心電図を用いて心室周期の任意の時相を限定して、任意の時相に相当する時系列波形の振幅変動を、ローレンツプロット法により隣り合う振幅変動ａ１，ａ２のプロット群を作図する。振幅ａ１，ａ２の集合体から最小二乗法により傾きｔａｎθ（以下、場合により、この傾きｔａｎθを「Ｆｒａｃｔａｌ  Ａｎｇｌｅ（ＦＡ）」と呼ぶ）を求める。この傾きｔａｎθの時間変動を時系列波形で表し、時系列波形の周波数解析結果を両対数表示して、周波数とＰＳＤに関する傾きを求める。具体的には、計測時間幅は前述の３６０秒間とし、３０秒間の傾きを求め、３０秒間の９０％ラップ、すなわち３秒毎にｔａｎθのゆらぎをプロットし、ｔａｎθの変動時系列波形を求める。ｔａｎθの変動時系列波形を周波数解析し、両対数表示図からＬＦ帯域の傾き（この傾きを「Ｆｒａｃｔａｌ  Ｓｌｏｐｅ（ＦＳ）」と呼ぶ）を求める。このＬＦ帯域の傾きＦＳが「＋」の値を取るか、「−」の値を取るかで、β１受容体の活性化の有無を判定できる。「＋」の値を取る場合は周期成分が強く強い制御が行われている状態で、「−」の値を取る場合は、フラクタル・ブラウン運動で１／ｆのゆらぎとなり、正常被験者が安静状態の定常状態において恒常性維持機能を発現している状態を示す。 
【００４２】
以上より、左心血流変動の全体の傾きは、ＶＬＦ〜ＬＦ帯域に現れ、また、交感神経活動が亢進した場合の全体の傾きと様相の異なる傾きもこれらの帯域に現れる。そして、交感神経活動が亢進した場合に発現する傾きは、「＋」と「−」の傾きがあり、特に「＋」の値を示す場合は心筋β１受容体に作用して心拍数は周波数成分に、血圧についてはＦＳに、周期性をもたらし、ゆらぎは線形化する。自律神経支配と体液性支配の二重支配構造をもつ心拍数（ＨＲ）制御は、心電図のＲＲＩの時系列波形の周波数解析結果のＰＳＤ−Ｈｚの両対数表示図から求められる。この心電図ＲＲＩの周波数解析によるゆらぎには、β１受容体の活性化による心筋の収縮力の増加による血流量の変動から生じるゆらぎ成分は含まれていない。心筋の収縮力変動による血流変動は、後負荷の心室の拡張圧の変動からも影響を受ける。したがって心室の拡張圧の変動は、血流変動のゆらぎと相関する。
【００４３】
心房収縮力は、静脈還流による血圧変動を超える圧力を生み出し、心室拡張圧、さらに心室拡張圧と平衡状態にある大動脈の拡張圧を超える時に、僧房弁閉鎖となる。僧房弁閉鎖直前の血流量変動のゆらぎと平均上腕拡張期血圧および平均上腕収縮期血圧に関する代数方程式を、ＦＳ（±ＬＦ）をパラメータにして表すと、
Δ（ｈ１−ｈ２）＝ｌ／ｇＡ・Ｑ（±ＬＦ）／Ｔｍ＝ｌ／ｇＡ・Ｆｆｓ（±ＬＦ）・・・（１５）
となる。 
【００４４】
したがって、 拡張期血圧（Ｄ．Ｂ．Ｐ）は、 
Ｈ（Ｄ．Ｂ．Ｐ）＝Δ（ｈ１−ｈ２）＋Ｈ０＝ｌ／ｇＡ・Ｑ・Ｆｆｓ（±ＬＦ）＋Ｈ０・・・（１６）
となる。 
【００４５】
ここに、Ｈ０はＦＳが０のときなので、安静状態では、高血圧と正常域血圧の境界の値となる。これは、正常被験者は、安静状態では、１／ｆゆらぎは−１となるためである。一方高血圧患者は、拡張期血圧が高値にあるため高心拍数、高心収縮力と相関関係にある正の傾きをもつＬＦ値を用いることになる。 
【００４６】
上記の考え方を心室収縮期に適応させると収縮期血圧の推定が可能になる。約０．３５秒の心室収縮期は等容性収縮期とそれより長い心室拍出期で構成される。図２４（ｂ）に示すように、等容性収縮期に左室圧が上昇し、大動脈圧を上回ると、大動脈弁が開き、左室から大動脈への血液の駆出が始まる。左室から拍出される血流は末梢の血管に流れ去る血流より大きなエネルギー量をもつため、弾性血管を拡張させ、一時的に留まることにより、動脈圧が上昇し、収縮期血圧に達することになる。そして動脈圧は左室圧を上回るが、慣性力により左室から大動脈への血液の流れは続く。血管の弾性の低下は血管拡張によるエネルギー損失が少なくし、末梢血管からの反射波により収縮期血圧をより高くする。そして急速拍出期の時間も短くなることから、収縮期血圧を高くする。そのため、高血圧者のＩＩ音は、大動脈弁閉鎖の時間短縮効果による水撃作用との合わせ技でＩ音よりも振幅が大きくなる可能性が示唆される。ここで式（９）を、心室収縮期の時相に合わせた式に変換する。
【００４７】
前負荷ｈ１がｈ３となり、時相は心房収縮から心室収縮にかわる。後負荷ｈ２はｈ４となる。Δｈ４が心室の収縮能を示し、収縮期血圧に関係するパラメータとなる。時相は、心電図Ｒ波の終わり近くから等容性収縮期の時間帯となる。ｄＱ／ｄｔは等容性収縮期の心臓の収縮性の指標としてよく用いられるｄＰ／ｄｔｍａｘと相関し、初期駆出血流変動は左室収縮期圧変動と相関する。動脈の収縮期血圧は、心拍変動と等容性収縮期の心室収縮能と良く相関し、ここでもＬＦをパラメータにした式が成立する。 
【００４８】
したがって、 収縮期血圧（Ｓ．Ｂ．Ｐ）は、ＦＳを因子とした関数Ｆｆｓを含む式で表され、 
Ｈ（Ｓ．Ｂ．Ｐ）＝Δｈ４＋Ｈ０＝ｌ／ｇＡ・Ｑ・Ｆｆｓ（±ＬＦ）＋Ｈ０・・・（１７）
となる。
【発明の効果】
【００４９】
本発明は、人の背部から採取される生体信号を解析し、心室充満期から等容性収縮期の血流変動によって生じる生体内の振動（生体内振動）を捉え、さらにこの生体内振動の様子を示す指標のゆらぎの様子（ゆらぎに関する指標（ゆらぎ指標））を捉える構成である。このゆらぎ指標は、上記のように人の状態、特に、血圧に関する情報と関連性がある。よって、本発明によれば、乗物のシート、事務用や家具用などの各種椅子、あるいは、ベッドなどの寝具等に取り付けられた、非拘束で生体信号を測定できる生体信号測定装置に、人が着席あるいは仰臥して背部を当該生体信号測定装置に当接するだけで、容易に血圧の推定、特に、正常域血圧（拡張期血圧が９０ｍｍＨｇ未満、収縮期血圧が１４０ｍｍＨｇ未満）を示す領域に含まれるか高血圧を示す領域に含まれるかの推定を行うことができる。また、カフを用いる血圧測定では、１回測定後、次の測定までに再度上腕を圧迫しなければならず、その間数分程度を要し、連続的な血圧測定を行うことは困難であったが、本発明によれば、非拘束であるため、背部から次々に捉えられる生体信号を用いて、連続的（ここでいう「連続的」には、カフを用いる場合よりも短い間隔で間欠的に測定する場合も含む）に血圧を推定することができる。
【図面の簡単な説明】
【００５０】
【図１】図１（ａ）は、本発明の一の実施形態において用いた背部体表脈波を測定する生体信号測定装置の一例を示した分解図であり、図１（ｂ）は、その要部断面図である。
【図２】図２は、本発明の一の実施形態に係る血圧推定装置の構成を模式的に示した図である。
【図３】図３（ａ）は、心電図波形を示した図であり、図３（ｂ）は、心音波を示した図であり、図３（ｃ）は、上記実施形態の生体信号測定装置から得られる生体信号の時系列波形をフィルタリングしたフィルタ処理波形（ＲＣＷ）を示した図である。
【図４】図４は、図３に示した心電図波形、心音波及びフィルタ処理波形（ＲＣＷ）の所定時間帯における拡大図である。
【図５】図５（ａ）は、心房収縮の際の血流変動を反映する波形成分に関する散布図であり、図５（ｂ）は、房室弁閉鎖の直前及び直後の心室充満期から等容性収縮期に移行する際の血流変動を反映する波形成分に関する散布図である。
【図６】図６（ａ），（ｂ）は、散布図傾き算出手段の機能を説明するための図である。
【図７】図７（ａ）〜（ｃ）は、ゆらぎ解析手段の機能を説明するための図である。
【図８】図８は、血圧推定用相関データを作成する手順を示したフローチャートである。
【図９】図９は、高血圧被験者４名の散布図及び散布図傾きを示した図である。
【図１０】図１０は、正常高値血圧被験者４名の散布図及び散布図傾きを示した図である。
【図１１】図１１は、至適血圧被験者８名の散布図及び散布図傾きを示した図である。
【図１２】図１２は、正常血圧被験者２名の散布図及び散布図傾きを示した図である。
【図１３】図１３は、低血圧被験者２名の散布図及び散布図傾きを示した図である。
【図１４】図１４は、心房収縮の際の血流変動を反映するａｎ（１），ａｎ（２）（ａ１ａ２）の波形成分を用いた散布図、散布図傾き、ゆらぎ解析図、並びに、それらから求めたフラクタルスロープ（ＦＳ）の例をまとめて示した図である。
【図１５】図１５は、心室充満期から等容性収縮期に移行する際の血流変動を反映するａｎ（２），ａｎ（３）（ａ２ａ３）の波形成分を用いた散布図、散布図傾き、ゆらぎ解析図、並びに、それらから求めたフラクタルスロープ（ＦＳ）の例をまとめて示した図である。
【図１６】図１６は、本実施形態の血圧推定用相関データの例を示した図である。
【図１７】図１７は、フラクタルスロープ（ＦＳ）の値の平均値、上腕血圧値の平均値を用いて求めた本実施形態の血圧推定用相関データの例を示した図である。
【図１８】図１８（ａ）は、図１６の血圧推定用相関データを再掲した図であり、図１８（ｂ）は、上腕血圧計の収縮期血圧及び拡張期血圧とフラクタルスロープ（ＦＳ）の対応関係を図１８（ａ）とは反対にして作成した比較データである。
【図１９】図１９は、血圧の推定過程を説明するためのフローチャートである。
【図２０】図２０は、血圧の推定過程を説明するためのフローチャートであり、図１９の続きである。
【図２１】図２１（ａ），（ｂ）は、フラクタルスロープを求める手法を説明するための図である。
【図２２】図２２（ａ），（ｂ）は、上記実施形態により推定した血圧の推定値と上腕血圧の実測値とを比較したグラフである。
【図２３】図２３（ａ），（ｂ）は、フラクタルスロープを求める回帰直線の傾きを図２２の場合と異ならせて求めた血圧の推定値と上腕血圧の実測値とを比較したグラフである。
【図２４】図２４（ａ），（ｂ）は、心室充満期から等容性収縮期で弁が解放していくまでの心房と心室の容積変化と血液の動きを示した図である。
【発明を実施するための形態】
【００５１】
以下、図面に示した本発明の実施形態に基づき、本発明をさらに詳細に説明する。本発明においては人の背部の体表面を伝播する生体信号（背部体表脈波）を用いる。この背部体表脈波は、上記のように、心房や心室の収縮期や拡張期において心房や心室内の血液の流れによる各弁や心筋内壁に衝突する際の振動、大動脈を流れる際に血管壁を押圧することによって生じる振動（音として伝わる振動を含む）が体表面に伝達されたものである。従って、背部体表脈波には、本発明においてターゲットとする心室充満期における生体内の振動、すなわち心室充満期から等容性収縮期の血流変動によって生じる生体内の振動（生体内振動）の情報が含まれる。
【００５２】
背部体表脈波を採取するための生体信号測定装置は、好ましくは、（株）デルタツーリング製の居眠り運転警告装置（スリープバスター（登録商標））で使用されている生体信号測定装置１を用いる。図１は生体信号測定装置１の概略構成を示したものである。この生体信号測定装置１は、図示しない電気的な配線等に加え、生体信号検出部１０を有してなり、該生体信号検出部１０が、測定用の椅子、ベッド、あるいは、乗物の運転席のシート等に組み込んで使用することができ、腕、手、指等を拘束することなく生体信号を採取できる。
【００５３】
生体信号検出部１０は、図１（ａ），（ｂ）に示したように、それぞれ所定の幅及び長さを有する略長方形の第一層部材１１、第二層部材１２及び第三層部材１３が積層された三層構造で構成されている。第一層部材１１は、三次元立体編物等からなり、生体信号の検出対象である人体側に位置させて用いられ、人体の背部の体表面を介しての生体信号は第一層部材１１にまず伝播される。第二層部材１２は、第一層部材１１から伝播される微弱な背部体表脈波を共鳴現象又はうなり現象によって強調させる共鳴層として機能し、ビーズ発泡体等からなるベース部材１２１、固有振動子の機能を果たす三次元立体編物１２２、膜振動を生じるフィルム１２３を有して構成される。ベース部材１２１の中心を挟んだ対称位置に、２つの配置孔１２１ａ，１２１ａが形成され、この配置孔１２１ａ，１２１ａに固有振動子の機能を果たす三次元立体編物１２２、１２２が配置されている。フィルム１２３，１２３は、この固有振動子の機能を果たす三次元立体編物１２２，１２２の各露出面を覆うように、第二層部材１２の両面に積層される。そして、一方の三次元立体編物１２２とフィルム１２３との間に背部体表脈波による振動（音）を検出するマイクロフォンセンサ１４が配設されている。なお、微弱な振動（音）である背部体表脈波を検出できるものであればマイクロフォンセンサ１４に限定されず他のセンサを用いてもよい。また、必要に応じて、２つの三次元立体編物１２２，１２２のそれぞれに対応してセンサを配置することもできる。第三層部材１３は、第二層部材１２を介して第一層部材１１の反対側に積層され、外部からの振動入力を低減する。第三層部材１３は、好ましくは１００Ｈｚを超える高周波の外部振動を除振する機能を有する。第三層部材１３としては、このようなフィルタリング機能を果たすために、第一層部材１１と同様に三次元立体編物を用いることが好ましい。
【００５４】
なお、三次元立体編物は、グランド編地の編目密度、グランド糸の太さや素材、連結糸の配設密度、連結糸の太さや素材などを調整することにより、必要な除振性能をもたせることができる。また、第一層部材１１として用いた三次元立体編物、第二層部材１２の固有振動子として用いた三次元立体編物１２２のいずれよりも、連結糸の配設密度を高くした三次元立体編物を第三層部材１３として用いることにより、所定以上の高周波振動を伝わり難くしている。
【００５５】
生体信号検出部１０は、上記のように第一層部材１１、第二層部材１２及び第三層部材１３を備えてなるが、背部の体表面を介して伝播される生体信号、特に、左心室の拡張、左心室の収縮に伴う大動脈への血液の流出による生体内の振動情報を含む背部体表脈波を検出するため、上記の椅子、ベッド等において、鎖骨に対応する高さから剣状突起に対応する高さまでの範囲に対応するに設けられることが好ましい。これは、本発明者が生体信号検出部１０の取り付け位置を種々変化させた中で、この範囲における検出感度が高かったことに基づくものである。但し、「心臓血液駆出に伴う胸部インピーダンス変化検出による電気的インピーダンス心拍出量計測用至適スポット電極配置」、生体医工学、４６巻６号（２００８年１２月）という研究報告によれば、心拍出量（ｃａｒｄｉａｃ  ｏｕｔｐｕｔ（ｃｏ））を電気的インピーダンス法を用いて計測するに当たって、電極を胸部前面に配置する場合、鎖骨から剣状突起に対応する範囲では、心房及び心室間の血液量変化の影響が大きく、動脈の血液量変化を捉えるには不十分で、鎖骨レベル、剣状突起レベルに電極を配置することが望ましいとする一方で、電極を背部に配置した場合には、鎖骨対応部位と剣状突起対応部位間では、一方から他方にずらしていくにしたがって、動脈の血液量が一定の率で変化していくと共に、心房及び心室間の血液量も一定の率で変化していくことをデータとして示している。このことは、背部においては、鎖骨対応部位と剣状突起対応部位間に生体信号検出部１０を当接すると、特に、左心室の拡張期における血液量の変化と、左心室が収縮して大動脈が駆出する際の大動脈における血液量の変化の両方を感度良く捉えることができることを示す間接的な裏付けとも言える。
【００５６】
次に、本実施形態の血圧推定装置１００の構成について図２に基づいて説明する。血圧推定装置１００は、生体信号処理手段２００を有して構成される。血圧推定装置１００は、コンピュータ（パーソナルコンピュータ、機器に組み込まれるマイクロコンピュータ等も含む）から構成され、生体信号処理手段２００として機能する生体信号処理手順を実行させるコンピュータプログラムが記憶部（当該コンピュータ（生体信号測定装置１００）の内蔵のハードディスク等の記録媒体のほか、リムーバブルの各種記録媒体、通信手段で接続された他のコンピュータの記録媒体等も含む）に記憶されている。また、生体信号処理手段２００として機能するコンピュータプログラムは、フィルタリング手段２１０、血流変動指標算出手段２２０、ゆらぎ指標算出手段２３０、推定手段２４０として機能し、フィルタリング手順、血流変動指標算出手順、ゆらぎ指標算出手順、推定手順をコンピュータに実行させる。また、フィルタリング手順、血流変動指標算出手順、ゆらぎ指標算出手順、推定手順を実現するコンピュータプログラムが組み込まれた１以上の記憶回路を有する電子回路で実現することもできる。
【００５７】
また、コンピュータプログラムは、記録媒体に記憶させて提供することができる。コンピュータプログラムを記憶した記録媒体は、非一過性の記録媒体であっても良い。非一過性の記録媒体は特に限定されないが、例えば フレキシブルディスク、ハードディスク、ＣＤ−ＲＯＭ、ＭＯ（光磁気ディスク）、ＤＶＤ−ＲＯＭ、メモリカードなどの記録媒体が挙げられる。また、通信回線を通じてコンピュータプログラムをコンピュータに伝送してインストールすることも可能である。
【００５８】
フィルタリング手段２１０は、生体信号測定装置１の生体信号検出部１０に組み込まれたセンサ１４から得られる生体信号の時系列波形（搬送波に生体情報がのり、生体信号検出部１０の第二層部材１２において共鳴現象等により強調された波形（以下、「センサ波形」））をフィルタリングする。フィルタリング手段２１０は、例えば、中心周波数２０Ｈｚ近傍のバンドパスフィルタ、好ましくは、１０〜３０Ｈｚの周波数帯域のバンドパスフィルタをかける手段であり、より好ましくは、高周波成分を除くため、これにさらに５０Ｈｚと３５Ｈｚのノッチフィルタをかける手段を有しており、これらのフィルタリングにより、センサ波形が、１０〜３０Ｈｚのフィルタ処理波形（図３（ｃ）の波形（以下、「ＲＣＷ」という）に変換される。心拍数は約１〜１．５Ｈｚ前後が標準的な範囲であるが、ＲＣＷは、図３（ｃ）に示したように、約１秒周期で相対的に大きな全振幅の波形成分が出現していることが捉えられていることから、心周期が顕在化される。この点は、図３（ａ）の心電図の波形、図３（ｂ）の心音図の波形と比較して明らかである。しかしながら、ＲＣＷには、心周期のリズムだけでなく、血流変動に伴う生体内の振動（音）のエネルギーが含まれている。
【００５９】
血流変動指標算出手段２２０は、ＲＣＷ中、心室充満期から等容性収縮期における時系列波形の範囲の波形成分を用いて、心室充満期から等容性収縮期における血流変動に起因して生体内で生じる振動に関する指標を求める手段である。血流変動指標算出手段２２０は、図２に示したように、散布図作成手段（散布図作成回路）２２１と、散布図傾き算出手段（散布図傾き算出回路）２２２とを有している。
【００６０】
散布図作成手段２２１は、ＲＣＷから上記の波形成分を２つ特定すると共に、この２つの波形成分の各全振幅（ａｎ（ｉ），ａｎ（ｉ＋１））を用いて散布図を作成する。心室充満期の心房収縮は、房室弁閉鎖前の約０．２秒間（例えば、安静時の人の平均的な１心周期が０．９秒程度（心拍数６７回）の場合で約０．１７秒間）である。ＲＣＷ中、この間の波形成分を用いる。房室弁の閉鎖は、図３（ａ）に示した心電図の波形のＲ波の直後であり、それを境として心房収縮から等容性収縮期に移行するため、心電図のＲ波のタイミングを基準とした波形成分を選択する。図４は、図３（ａ）〜（ｃ）における約１０１秒〜約１０１．８秒の範囲の拡大図である。図４に示したように、まず、心電図のＲ波の直前に位置する房室弁閉鎖の直前の極点Ｂｎ（０）を求める。次に、極点Ｂｎ（０）の一つ前に位置する極点をＢｎ（１）とし、Ｂｎ（１）を基準として、その直前の波形成分の下側極点から上側極点までの縦軸に沿った振幅（すなわち全振幅（ｐ−ｐ））を「ａｎ（１）」とし、その直後の波形成分の全振幅を「ａｎ（２）」とする。ここで特定した全振幅：ａｎ（１），ａｎ（２）の波形成分はＲ波の直前０．１秒以内の範囲であり、この２つの波形成分の各全振幅は、心房収縮の際の血流変動を反映している。
【００６１】
一方、Ｒ波の直後の極点、すなわち、極点Ｂｎ（０）の直後に位置する極点をＢｎ（２）とし、極点Ｂｎ（２）を基準として、その直前の波形成分の全振幅を「ａｎ（２）」とし、その直後の全振幅を「ａｎ（３）」とする。これらは、心電図のＲ波を挟んだ直前及び直後の波形成分であり、心室充満期から等容性収縮期に移行する際の血流変動を反映している。
【００６２】
次に、散布図作成手段２２１は、このようにして求めた隣接する２つの波形成分の各全振幅（ａｎ（ｉ），ａｎ（ｉ＋１））の値を、一方の値を縦軸に、他方の値を横軸にとった座標系にプロットし、散布図（ローレンツプロット図）を作成する（図５参照）。この散布図にプロットされた点が、２つの波形成分の各全振幅の比を示す。
【００６３】
ローレンツプロット法を用いることにより、Δｔの間に流れる血液の変化量を求めることができる。すなわち、心房の収縮能により、血液に運動エネルギーが与えられ、この血流の振動系は増幅する自由振動になる。自由振動系は、一般的に粘性減衰と固体摩擦の二つの減衰が作用するが、ローレンツプロット法を用いて、不規則なゆらぎを見せる血流量の中で規則性を見出すことができる。
【００６４】
上記のように、散布図を作成するための２つの波形成分を特定するためには、房室弁閉鎖直前の心電図波形に出現するＲ波を特定する必要がある。本実施形態では、図２に示したように、心電計３００からの心電図波形データを血圧推定装置１００で受信し、その情報を血流変動指標算出手段２２０の散布図作成手段２２１で上記のように用いるようになっている。なお、心電計３００としては、Ｒ波のタイミングを特定できるものであればその種類は全く限定されるものではない。よって、本実施形態では、被験者の背部に生体信号測定装置１の生体信号検出部１０を配置すると共に、心電計３００の電極を被験者の胸部に装着してデータを採取する。
【００６５】
図５（ａ）は、上記のようにして求めた心房収縮の際の血流変動を反映する２つの波形成分の各全振幅ａｎ（１）、ａｎ（２）をプロットした散布図（ローレンツプロット図）であり、図５（ｂ）は、Ｒ波の出現直後である心室充満期から等容性収縮期に移行する際の血流変動を反映する２つの波形成分の各全振幅ａｎ（２）、ａｎ（３）をプロットした散布図である。いずれも、ＲＣＷの３０秒間のデータを一つの座標系にプロットしたものである（例えば心拍数約６０回／分の場合で約３０点プロットされる）。散布図作成手段２２１は、このような散布図をＲＣＷについて時系列に求めていくが、その際、ＲＣＷの時間幅を例えば９０％ラップで重ねてスライドさせて求めていく。これにより、本例では全振幅の比を示す散布図が３秒間隔で作成されることになる。なお、この時間幅やスライド計算のラップ率はあくまで一例であり、これに限定されるものではない。但し、設定する時間幅が長すぎたり、ラップ率が低すぎたりすると、心房収縮の状態の把握精度が低下するため、時間幅は、１０〜６０秒程度の範囲とし、ラップ率は、７０〜９５％程度の範囲とすることが好ましい。
【００６６】
散布図傾き算出手段２２２は、散布図にプロットされた点群について、最小二乗法により回帰直線Ａを引き、さらにその回帰直線Ａの傾きを求める。なお、この傾きは、横軸又は縦軸に対する傾きをそのまま用いることもできるが、図６（ｂ）に示したように、任意の基準線（総時間プロット基準線）に対する傾き（総時間プロット基準線と３０秒間の回帰直線Ａとの間の角度）を、散布図傾き算出手段２２２から出力される傾き（θｉ）とすることが好ましい。また、回帰直線Ａの傾きｔａｎθはθ（度）に近似していると見なして差し支えないため本実施形態では、θ（度）のゆらぎの時系列波形を計算に用いる。
【００６７】
総時間プロット基準線は、ＲＣＷについて、上記の各散布図を作成する時間幅よりも数倍以上長い時間（例えば計測総時間）について全振幅比を求め、その長い時間に関する散布図に引いた線である。具体的には、図６（ａ）に示したように、例えば計測総時間３６０秒間のデータの中で、ボックスカウンティング法により出現回数の少ないデータ（例えば２０％以下）は棄却し、全振幅の小さい波形のプロットが集まっている点（プロット群の下端）を起点として、プロットの重心を通る絶対最大長の線を総時間プロット基準線とする。プロットの重心は、Ｘ座標（ａｎ（ｉ＋１））、Ｙ座標（ａｎ（ｉ））の平均値である。
【００６８】
ゆらぎ指標算出手段２３０は、血流変動指標算出手段２２０の散布図傾き算出手段２２２により求められた生体内で生じる振動に関する指標（血流変動指標）の時系列波形を求めると共に、その時系列波形のゆらぎに関する指標を求める手段である。ゆらぎ指標算出手段２３０は、図２に示したように、ゆらぎ解析手段２３１とゆらぎ解析図傾き算出手段２３２とを有している。
【００６９】
ゆらぎ解析手段２３１は、まず、散布図傾き算出手段２２２により求められた点群の回帰直線の傾き（θｉ）の値（血流変動指標）を時系列にプロットする。上記の例では、点群の回帰直線の傾き（θｉ）の値（血流変動指標）が３秒間隔で求められるため、図７（ａ）に示したように、当該回帰直線の傾き（θｉ）の値を３秒間隔でプロットし、図７（ｂ）に示した時系列波形を改めて構成する。なお、図７（ｂ）では、回帰直線の傾き（θｉ）のそのままの値をプロットした時系列波形（全体の平均値をゼロ点にして引き直した時系列波形でもよい）を「θｉ時系列波形」とし、「θｉ時系列波形」に０．０８Ｈｚのローパスフィルターをかけた時系列波形を「ＬＰＦ処理波形」として示している。上記の例の場合、３秒に１回のプロットになるため、０．０８Ｈｚを超える高周波成分は信頼性が低いと考えられる。この点群の回帰直線の傾き（θｉ）の値（血流変動指標）の時系列波形は、ある状態から次の状態へと移り変わることに対応するためのゆらぎがある。そこで、ゆらぎ解析手段２３１は、血流変動指標の時系列波形について周波数解析（ＦＦＴ）を行って、パワースペクトル密度と周波数の両対数グラフに、周波数解析結果を表示する。なお、周波数解析の対象とする回帰直線の傾き（θｉ）の値（血流変動指標）の時系列波形は、図７（ｂ）における「θｉ時系列波形」よりも「ＬＰＦ処理波形」を用いることが好ましい。この周波数解析により、ゆらぎを支配する周波数がいずれであるかがわかる（図７（ｃ）参照）。本明細書では、この周波数解析の結果を「ゆらぎ解析図」と称する。
【００７０】
ゆらぎ解析図傾き算出手段２３２は、ゆらぎ解析手段２３１で得られたゆらぎ解析図の回帰直線Ｂを求めると共に、このゆらぎ解析図の回帰直線Ｂの傾き（本願では、このゆらぎ解析図の回帰直線Ｂの傾きを「フラクタルスロープ（ＦＳ）」と称する）を算出する。ゆらぎ解析図の回帰直線Ｂは、本実施形態では、図７（ｃ）に示したように、パワースペクトル密度の対数が周波数の対数に全体的には反比例するため、パワースペクトル密度の対数目盛りを縦軸に周波数の対数目盛りを横軸にとると右肩下がりになる。例えば、４５度で右肩下がりのゆらぎ解析図の回帰直線Ｂの傾き（フラクタルスロープ）は「−１（すなわち１／ｆゆらぎの傾き）」で示すことができる。
【００７１】
上記のように、左心血流変動の全体の傾きは、ＶＬＦからＬＦ帯域に現れ、特に交感神経活動が亢進した場合は全体の傾きと様相の異なる傾きがこれらの帯域に現れ、また、これらの帯域の傾きが「＋」の値を取るか、「−」の値を取るかで、高血圧患者に特有のβ１受容体の活性化の有無を判定できる。さらに、０．００３〜０．０８Ｈｚ付近が血圧調節ホルモン（レニン・アンジオテンシン・アルドステロン系）の働きを反映する周波数帯であり、Ｍａｙｅｒ波と呼ばれる０．１Ｈｚの振動（交感神経性血管収縮神経の興奮レベルの周期的振動）も血圧に影響を及ぼしている。これらのことを勘案すると、ゆらぎ指標であるゆらぎ解析図の回帰直線Ｂの傾き（フラクタルスロープ）は、ＶＬＦ帯域（超低周波帯域：０．００３３〜０．０４Ｈｚ）からＬＦ帯域（低周波帯域：０．０４〜０．１５Ｈｚ）において設定される所定の周波数帯で求めることが好ましい。所定の周波数帯は、本実施形態の場合、上記のように血流変動指標の時系列波形を作成する際に０．０８Ｈｚでローパスフィルターをかけているため、まず上限値は０．０８Ｈｚ付近に設定することが好ましい。一方、下限値は、０．０１Ｈｚ以下の情報には血流量の変動情報が顕著には含まれていないことから、０．０１〜０．０４Ｈｚの間で設定することが好ましい。例えば、心拍数が早くなることにより血圧が低下傾向となる場合があり、その際には、血流量変動のゆらぎが低周波側に現れる傾向となる場合もある。よって、所定の周波数帯は、例えば、０．０３〜０．０８Ｈｚと一律に設定することも可能であるが、血圧を推定する被験者毎に設定することが好ましい。
血圧を推定する被験者毎に設定する方法として、本実施形態では次の方法を用いた。すなわち、後述するように、血圧の推定に当たっては、被験者のゆらぎ指標であるフラクタルスロープ（ＦＳ）を、血圧推定用相関データに照合して行われる。よって、上腕血圧値とフラクタルスロープ（ＦＳ）とを対応させて血圧推定用相関データ（血圧推定用相関式）を予め求めて記憶部に設定しておく必要がある。従って、この血圧推定相関データの作成には、各被験者のゆらぎ解析図毎に、被験者の心拍数等を考慮して個別に、フラクタルスロープ（ＦＳ）を求めるための適切な周波数帯が設定されている。そこで、これらのゆらぎ解析図を教師データとして用い、これから血圧を推定する被験者のゆらぎ解析図を各教師データに照合し、一致度の高い教師データを抽出し、その教師データに用いられている周波数帯を、これから血圧を推定する被験者のゆらぎ解析図において用い、フラクタルスロープ（ＦＳ）を算出することができる。このような手法を用いることで、血圧推定用相関データを作成した後は、所定の周波数帯を被験者毎に自動で設定することができる。なお、新たに求められたデータは、新たな教師データとして蓄積し、その新たな教師データも含めて血圧推定用相関データ（血圧推定用相関式）を改めて求める、というステップを繰り返すことで、血圧の推定精度を高めていくことができる。
【００７２】
ゆらぎ解析図傾き算出手段２３２により出力されるゆらぎ指標としてのフラクタルスロープ（ＦＳ）は、後述する推定手段２４０において血圧の推定に用いられるが、そのためには予め、血圧とフラクタルスロープ（ＦＳ）との相関データ（血圧推定用相関データ）を作成しておく必要がある。従って、本実施形態において生体信号処理手段２００は、ゆらぎ解析図傾き算出手段２３２により出力されるゆらぎ指標としてのフラクタルスロープを用いて、血圧推定用相関データを作成する相関データ作成手段２６０を有している（図２）。この相関データ作成手段２６０では、次のような手順が実行されて相関データを作成し、記憶部に記憶させる。
【００７３】
なお、血圧推定用相関データを作成する場合、散布図傾き算出手段２２２により求められた点群の回帰直線の傾き（θｉ）の値（血流変動指標）として、拡張期血圧の推定には、心房収縮の際の血流変動を反映する２つの波形成分の各全振幅ａｎ（１）、ａｎ（２）をプロットした散布図から得られるものを採用し、収縮期血圧の推定には、房室弁閉鎖の直前及び直後の心室充満期から等容性収縮期に移行する際の血流変動を反映する２つの波形成分の各全振幅ａｎ（２）、ａｎ（３）をプロットした散布図から得られるものを採用する。
【００７４】
そして、図８に示したように、ゆらぎ指標算出手段２３０により回帰直線Ｂの傾き（フラクタルスロープ（ＦＳ））を算出したならば（Ｓ１１０１）、相関データ作成手段２６０がこの情報を受け取る（Ｓ１１０３）。また、ほぼ同じタイミングで血圧計により測定した被験者の上腕血圧のデータ（Ｓ１１０２）も受信する（Ｓ１１０３）。なお、上腕血圧の測定方法は限定されるものではなく、上腕にカフを巻き付けるタイプの血圧計（上腕式血圧計）によって測定してもよいし、手首式血圧計等により測定してもよいことはもちろんである。このようにして、血圧推定用相関データ（血圧／ｆｒａｃｔａｌ  ｓｌｏｐｅ  ＭＡＰ）が作成され（Ｓ１１０４）、このデータが記憶部に記憶される。後述する図１６あるいは図１７が血圧推定用相関データの例であり、図１６及び図１７ともに、横軸にＬＦ帯域の回帰直線Ｂの傾き（フラクタルスロープ（ＦＳ））をとり、縦軸に上腕血圧データをとり、拡張期血圧推定用相関データと収縮期血圧推定用相関データをそれぞれ示す二本の相関線（相関式）がそれぞれ求められている。
【００７５】
推定手段２４０は、ゆらぎ解析図傾き算出手段２３２により求められたフラクタルスロープの値をゆらぎ指標として用い、血圧の推定を行う。すなわち、評価対象者のフラクタルスロープ（ＦＳ）の値を、記憶部に記憶された上記の血圧推定用相関データに照合して血圧を推定する（図２参照）。例えば、ａｎ（１）（図中及び以下においては単に「ａ１」と表記）及びａｎ（２）（図中及び以下においては単に「ａ２」と表記）を用いて求めたゆらぎ指標であるフラクタルスロープ（ＦＳ）が「−０．２」と出力された場合、図１６では、ｙ＝１９．０３９ｘ＋８４．３６１の式に代入され、約８１ｍｍＨｇと求められ、図１７を用いた場合には、ｙ＝２０．４１４ｘ＋８５．０４３の式に代入され、約８１ｍｍＨｇと推定される。また、ａ２及びａｎ（３）（図中及び以下においては単に「ａ３」と表記）を用いて求めたゆらぎ指標であるフラクタルスロープ（ＦＳ）が「−１」と出力された場合、図１６では、ｙ＝１９．８１３ｘ＋１４２．５の式に代入され、約１２３ｍｍＨｇと推定され、図１７に照らすと、ｙ＝１９．００６ｘ＋１４１．９３の式に代入され、約１２３ｍｍＨｇと推定される。ａ１及びａ２を用いたフラクタルスロープ（ＦＳ）、並びに、ａ２及びａ３を用いたフラクタルスロープ（ＦＳ）のいずれも、本実施形態によれば、次々に求めることができるため、被験者の血圧を、非拘束でありながら連続的に推定することができる。
【００７６】
ここで、図９〜図１５を用いて図１６及び図１７に示した血圧推定用相関データの作成過程を説明する。まず、図９〜図１３は、代表的な被験者２０人（被験者（Ｓｕｂ．）Ｎｏ．０１〜２０）の背部体表脈波をフィルタリング手段２１０によりフィルタリングした後、血流量変動指標算出手段２２０の散布図作成手段２２１によって作成した散布図（ローレンツプロット図）を示している。なお、これらの被験者から採油した背部体表脈波は安静状態のものである。これらは、いずれも、ａ１及びａ２を用いて求めた点（図の凡例中、「ａ１ａ２」と表記）とａ２及びａ３を用いて求めた点（図の凡例中、「ａ２ａ３」と表記）の両方がプロットされている。また、「ａ１ａ２」点、「ａ２ａ３」点の各点群について、散布図傾き算出手段２２２によって求めた回帰直線Ａ（それぞれ、「ａ１ａ２近似線」、「ａ２ａ３近似線」）が示されている。そして、散布図傾き算出手段２２２は、これらａ１ａ２近似線及びａ２ａ３近似線の傾き（θ１）を求める。
【００７７】
なお、図９〜図１５の分析に当たっては、日本高血圧学会の高血圧ガイドライン（２０１４年）の血圧分類に基づいて、安静状態における各被験者の上腕血圧値から、次のように血圧の種類を分類した。まず、収縮期血圧１４０ｍｍＨｇ以上かつ／又は拡張期血圧９０ｍｍＨｇ以上の場合に「高血圧」と分類し、「高血圧」の条件未満の場合に「正常域血圧」と分類した。また、「正常域血圧」について、収縮期血圧１３０〜１３９ｍｍＨｇかつ／又は拡張期血圧：８５〜８９ｍｍＨｇの場合を「正常高値血圧」とし、収縮期血圧１２０〜１２９ｍｍＨｇかつ／又は拡張期血圧：８０〜８４ｍｍＨｇの場合を「正常血圧」とし、収縮期血圧１２０ｍｍＨｇ未満かつ拡張期血圧８０ｍｍＨｇ未満の場合を「至適血圧」とした。「低血圧」については上記ガイドラインにおいて規定がないため、世界保健機関（ＷＨＯ）に規定されている収縮期血圧１００ｍｍＨｇ以下、拡張期血圧６０ｍｍＨｇ以下の基準を採用した。
【００７８】
次に、ゆらぎ指標算出手段２３０のゆらぎ解析手段２３１が、散布図傾き算出手段２２２により求めたａ１ａ２近似線及びａ２ａ３近似線の各傾き（θ１）の時系列波形を求め（図７（ｂ）参照）、周波数解析（ＦＦＴ）を行って、パワースペクトル密度と周波数の両対数グラフからなるゆらぎ解析図を求める（図７（ｃ）参照）。ゆらぎ解析図傾き算出手段２３２は、上記のようにＶＬＦからＬＦ帯域において、各被験者の心拍数等を参照して所定の周波数帯を設定し、ゆらぎ解析図の回帰直線Ｂの傾きであるフラクタルスロープ（ＦＳ）を求める。
【００７９】
次に、相関データ作成手段２６０により、拡張期血圧に対応するａ１ａ２に関しては図１４に示したようなフラクタルスロープ（ＦＳ）が得られ、収縮期血圧に対応するａ２ａ３に関しては図１５に示したようなフラクタルスロープが得られる。図１４及び図１５中、中央列（「（２）」で表示した列）が上記のゆらぎ解析図に相当し、そこにおける所定の周波数帯のフラクタルスロープ（ＦＳ）を求めている。なお、図１４及び図１５中、左列（「（１）」で表示した列）は散布図（ローレンツプロット図）を示し、右列（「（３）」で表示した列は、心電計３００から得られた心拍数の時系列波形のゆらぎを解析したスペクトルである。心拍数の時系列波形のゆらぎを解析したスペクトルのフラクタルスロープは血圧との有意な相関は見られなかった。
（２）の中央列を見ると、例えば、正常高値血圧の範囲に分類される血圧を示す被験者Ｎｏ．０５は、図１４のａ１ａ２に基づいたものは「−０．０５」のフラクタルスロープ（ＦＳ）が得られ、図１５のａ２ａ３に基づいたものは、「−０．７８」のフラクタルスロープ（ＦＳ）が得られた。この際の被験者Ｎｏ．０５の上腕血圧は、拡張期血圧が８３、収縮期血圧が１３３であった。そして、図１６において、ａ１ａ２に関するフラクタルスロープ（ＦＳ）の値「−０．０５」に、拡張期血圧８３を対応させてプロットし、ａ２ａ３に関するフラクタルスロープ（ＦＳ）の値「−０．７８」に、収縮期血圧１３３を対応させてプロットする。
【００８０】
同様に、至適血圧の範囲に分類される血圧を示す被験者Ｎｏ．１１は、図１４のａ１ａ２に基づいたものは「−０．５３」のフラクタルスロープ（ＦＳ）が得られ、図１５のａ２ａ３に基づいたものは、「−１．３０」のフラクタルスロープ（ＦＳ）が得られた。この際の被験者Ｎｏ．１１の上腕血圧は、拡張期血圧が６５、収縮期血圧が１０８であった。そして、図１６において、ａ１ａ２に関するフラクタルスロープ（ＦＳ）の値「−０．５３」に、拡張期血圧６５を対応させてプロットし、ａ２ａ３に関するフラクタルスロープ（ＦＳ）の値「−１．３０」に、収縮期血圧１０８を対応させてプロットする。
さらに、高血圧の範囲に分類される血圧を示す被験者Ｎｏ．０１は、図１４のａ１ａ２に基づいたものは「＋１．８１」のフラクタルスロープ（ＦＳ）が得られ、図１５のａ２ａ３に基づいたものは、「＋０．５４」のフラクタルスロープ（ＦＳ）が得られた。この際の被験者Ｎｏ．０１の上腕血圧は、拡張期血圧が１１６、収縮期血圧が１６４であった。そして、図１６において、ａ１ａ２に関するフラクタルスロープ（ＦＳ）の値「＋１．８１」に、拡張期血圧１１６を対応させてプロットし、ａ２ａ３に関するフラクタルスロープ（ＦＳ）の値「＋０．５４」に、収縮期血圧１６４を対応させてプロットする。
【００８１】
その他の被験者についても同様の処理を行い、図１６にプロットしていく。このようにしてプロットしていくと、高血圧に分類される被験者のフラクタルスロープ（ＦＳ）の値はａ１ａ２及びａ２ａ３のいずれに関してもほぼ０以上となり、正常域血圧（正常高値血圧、正常血圧、至適血圧、低血圧）の被験者のフラクタルスロープ（ＦＳ）の値はａ１ａ２及びａ２ａ３のいずれに関してもほぼ０以下となった。そして、図１６にプロットされたａ１ａ２に基づいた点群に引いた回帰直線はｙ＝１９．０３９ｘ＋８４．３６１となり、決定係数は０．７５１１で高い相関性が示された。また、フラクタルスロープ（ＦＳ）値が「０」のときのｙ切片の値が「８４．３６１」であり、拡張期血圧の高血圧と正常域血圧との境界である９０ｍｍＨｇに近い値であると共に、正常高値血圧と正常血圧との境界である８５ｍｍＨｇとほぼ同じ値であった。この回帰直線が最終的に求める拡張期血圧推定用相関データとなる。
【００８２】
一方、ａ２ａ３に基づいた点群に引いた回帰直線はｙ＝１９．８８２ｘ＋１４２．４７となり、決定係数は０．９１０１で高い相関性が示された。また、フラクタルスロープ（ＦＳ）値が「０」のときのｙ切片の値が「１４２．４７」であり、収縮期血圧の高血圧と正常域血圧との境界である１４０ｍｍＨｇに近い値であった。この回帰直線が最終的に求める収縮期血圧推定用相関データとなる。なお、図１６にプロットされたデータは、ａ１ａ２に関しては被験者２５名分の１９２データであり、ａ２ａ３に関しては被験者３４名分の１４６データである。
【００８３】
図１７は、図１６にプロットしたデータのうち、複数回計測を行った被験者のａ１ａ２及びａ２ａ３に関する各フラクタルスロープ（ＦＳ）の値の平均値、上腕血圧計により測定された拡張期血圧の平均値と収縮期血圧の平均値を用いて作成した血圧推定用相関データである。ａ１ａ２を用いた拡張期血圧推定用相関データの回帰直線はｙ＝２０．４１４ｘ＋８５．０４３であり、その決定係数は０．９２４６であった。ａ２ａ３を用いた収縮期血圧推定用相関データの回帰直線はｙ＝１９．２２２ｘ＋１４２．０８であり、その決定係数は０．９６１５であった。いずれも、図１６で求めた回帰直線よりも相関性が高くなっており、血圧推定用相関データの作成に当たっては、被験者の複数の測定値の平均値を用いることが好ましい。
【００８４】
ここで、図１８（ａ）は、図１８（ｂ）との比較のために、図１６の血圧推定用相関データを再掲したものである。図１８（ｂ）は、図１８（ａ）とは逆に、ａ１ａ２に基づくフラクタルスロープ（ＦＳ）の値に対して、上腕血圧計により測定した収縮期血圧を対応させてプロットし、ａ２ａ３に基づくフラクタルスロープ（ＦＳ）の値に対して、上腕血圧計により測定した拡張期血圧を対応させてプロットした各点群を示したものである。図１８（ｂ）では、ａ１ａ２に基づくフラクタルスロープ（ＦＳ）の値に収縮期血圧を対比させた点群の回帰直線が、ｙ＝２２．４９９ｘ＋１３１．３４、決定係数が０．５２４４と、図１８（ａ）と比較して相関性が大きく劣っていた。また、ａ２ａ３に基づくフラクタルスロープ（ＦＳ）の値に拡張期血圧を対比させた点群の回帰直線が、ｙ＝１３．２７２ｘ＋９０．０３７、決定係数が０．７３０９であった。決定係数の値が若干低くなっていると共に、ｙ切片が９０を若干上回っていた。また、図１８（ａ）では収縮期血圧と拡張期血圧の各回帰直線の傾きがほぼ同じであるのに対し、図１８（ｂ）の手法では、収縮期血圧と拡張期血圧の各回帰直線の傾きが大きく異なっていることがわかる。換言すれば、図１８（ｂ）の手法の場合、収縮期血圧と拡張期血圧との差がほとんどないデータ構成が存在することになり、両者を明確に区別して推定するための相関データとしての信頼性が乏しいことがわかる。よって、本実施形態のように、拡張期血圧推定用にはａ１ａ２に基づくフラクタルスロープ（ＦＳ）の値を用い、収縮期血圧推定用にはａ２ａ３に基づくフラクタルスロープ（ＦＳ）の値を用いることが適切である。
【００８５】
次に、図１９及び図２０のフローチャートに基づき、血圧の推定過程を具体的に説明する。まず、血圧を推定する対象となる人（血圧推定対象者）は、背部体表脈波を検出する生体信号検出装置１の生体信号検出部１０が上記のように背部に設けられた測定用の椅子等に着席する。また、血圧推定対象者には心電計３００の電極を胸部に装着し、心電図波形データが血圧推定装置１００に送信されるように設定する。
【００８６】
測定が開始されると、血圧推定装置１００は、血圧推定対象者の背部体表脈波（センサ波形）を生体信号測定装置１により受信する（Ｓ１９０１）。すると、血圧推定装置１００におけるフィルタリング手段２１０が、受信したセンサ波形をフィルタリングし（Ｓ１９０２）、ＲＣＷを求める（Ｓ１９０３）。次に、血流量変動指標算出手段２２０の散布図作成手段２２１は、図４に示したように、心電計３００から得た心電図波形データにおいてＲ波のタイミングを特定し（Ｓ１９０４）、それに基づいて上記の各全振幅ａ１，ａ２，ａ３を求める（Ｓ１９０５）。次に、散布図作成手段２２１は、ａ１及びａ２を用いた散布図（ローレンツプロット）、ａ２及びａ３を用いた散布図（ローレンツプロット）を作成する（Ｓ１９０６，Ｓ１９０７）。次に、散布図傾き算出手段２２２により、各散布図の回帰直線の傾きθｉを求め、ゆらぎ指標算出手段２３０のゆらぎ解析手段２３１がθｉの時系列波形を拡張期用、収縮期用のそれぞれについて作成する（Ｓ１９０８）。次に、ゆらぎ解析手段２３１は、周波数解析（スペクトル解析）を行い（Ｓ１９０９）、θｉ時系列波形のスペクトルであるゆらぎ解析図を拡張期用、収縮期用のそれぞれについて作成する（Ｓ１９１０）。
【００８７】
次に、ゆらぎ解析図傾き算出手段２３２が、Ｓ１９１０で求めた血圧推定対象者の拡張期用、収縮期用の各ゆらぎ解析図（スペクトル）を、記憶部のデータベースにアクセスし（図２参照）、ゆらぎ解析図（スペクトル）の教師データに照合する（Ｓ１９１１）。各教師データは、上記のように、被験者毎に応じて適切な回帰直線Ｂの傾き（フラクタルスロープ）を求めるために、ＶＬＦからＬＦに属する範囲でそれぞれの特性に合わせた所定の周波数帯が設定されている。そこで、Ｓ１９１１では、血圧推定対象者のＳ１９１０で求めたゆらぎ解析図（拡張期用、収縮期用）を、各教師データに設定された周波数帯の範囲において照合し、正規化相互相関関数を用いて一致度の高い教師データを抽出する。次に、ゆらぎ解析図傾き算出手段２３２が、抽出された教師データに設定されている周波数帯を、Ｓ１９１０で求めた血圧推定対象者のゆらぎ解析図（拡張期用、収縮期用）に適用し（Ｓ２００１）、その範囲における回帰直線Ｂの傾きをフラクタルスロープ（ＦＳ）として求める（Ｓ２００２、Ｓ２００３）。図２１は、その一例を示しており、図２１（ｂ）が血圧推定対象者のゆらぎ解析図（スペクトル）である。図２１（ａ）が教師データの一例であり、太線がこの教師データの回帰直線Ｂを示し、この教師データの回帰直線Ｂを引いた周波数帯と同じ周波数帯においてゆらぎ解析図（スペクトル）の一致度が比較される。教師データに設定された周波数帯は、上記のように教師データ毎に異なるため、この一致度を全ての教師データについて比較し、一致度の高い教師データに設定された周波数帯と同じ周波数帯を図２１（ｂ）の血圧推定対象者のゆらぎ解析図に設定し、その周波数帯について回帰直線Ｂを引いてその傾きをフラクタルスロープとして求める。
なお、上記の教師データとの一致度の比較法はあくまで一例であり、一致度を比較するに当たって、例えば、ＶＬＦからＬＦに属するゆらぎ解析図（スペクトル）の全体について一致度を判定し、それによって一致度が高いと判定された教師データに設定された周波数帯を血圧推定対象者のゆらぎ解析図に利用するようにすることもできる。
【００８８】
次に、推定手段２４０が動作して記憶部に記憶されている図１６又は図１７に示された血圧推定用相関データにアクセスし（Ｓ２００４）、拡張期用、収縮期用の各血圧推定用相関式に各フラクタルスロープ（ＦＳ）を代入する。これにより血圧推定対象者の拡張期血圧、収縮期血圧を求めることができる（Ｓ２００５，Ｓ２００６）。
【００８９】
図２２（ａ）は、図１９〜図２０に示したステップに従って推定した複数の血圧推定対象者（２５名、４９データ）の収縮期血圧と上腕血圧の実測値との比較を示した図である。図２２（ａ）は、血圧推定用相関データとして、図１６の全事例を用いて求めたものを利用して血圧を推定し、図２２（ｂ）は、血圧推定用相関データとして、図１７の平均値を利用して求めたものを利用して血圧を推定したものである。いずれも、傾きが１に近いと共に、決定係数も高く、本実施形態の手法により推定される血圧が、上腕血圧と極めて高い相関を有することがわかる。
【００９０】
一方、図２３（ａ），（ｂ）は、血圧推定対象者（２７名、７５データ）について、教師データを参照して設定した所定の周波数帯で傾きを算出する際（図２０のＳ２００１，Ｓ２００２）、その周波数帯におけるゆらぎ解析図の波形成分全体に関する回帰直線の傾き（図において単に「傾き」と表記）ではなく、図２１（ｂ）に示したように当該波形成分における図の上側のピーク点のみに対する回帰直線の傾き（図において「ｐｅａｋ傾き」と表記）を求めて推定した収縮期血圧と上腕血圧の実測値との比較を示した図である。なお、図２３（ａ）は、図１６の全事例を用いた血圧推定用相関データを利用し、図２３（ｂ）は、図１７の平均値を用いた血圧推定用相関データを利用したものである。図２３（ａ），（ｂ）から明らかなように、この場合の推定値は実測値との相関が極めて低かった。よって、上記のようにゆらぎ解析図における所定の周波数帯の波形成分に関し、その全体の回帰直線Ｂの傾きをフラクタルスロープ（ＦＳ）として採用することが好ましい。
【００９１】
本実施形態によれば、上記のように、複数人の被験者の上腕血圧とフラクタルスロープ（ＦＳ）との関係で血圧推定用相関データ（拡張期血圧推定用、収縮血圧推定用）を作成して記憶しておけば、その後は、血圧推定対象者のフラクタルスロープ（ＦＳ）を求めるだけで、この相関データに照らして、血圧を、非拘束、リアルタイムで推定でき、かつ連続的に推定することができる。
【００９２】
よって、例えば、生体信号測定装置１の生体信号検出部１０をベッド、椅子の背部その他の人体支持手段に設けておけば、被験者の血圧をいつでも推定することができる。従って、被験者の血圧の急変も察知することができ、睡眠中の血圧の変動も知ることができる。もちろん、これらの推定した血圧情報は、血圧推定装置１００に内蔵若しくは接続した通信手段（図示せず）により、モニタ、警告灯などの任意の表示装置にその旨を知らせたり、管理者等の端末装置に通知したりすることもできる。
【符号の説明】
【００９３】
１  生体信号測定装置
１０  生体信号検出部
１１  コアパッド
１２  スペーサパッド
１４  センサ
１００  血圧推定装置
２００  生体信号処理手段
２１０  フィルタリング手段
２２０  血流変動指標算出手段
２２１  散布図作成手段
２２２  散布図傾き算出手段
２３０  ゆらぎ指標算出手段
２３１  ゆらぎ解析手段
２３２  ゆらぎ解析図傾き算出手段
２４０  推定手段
【図１】
2019122502.tif 000003
【図２】
2019122502.tif 000004
【図３】
2019122502.tif 000005
【図４】
2019122502.tif 000006
【図５】
2019122502.tif 000007
【図６】
2019122502.tif 000008
【図７】
2019122502.tif 000009
【図８】
2019122502.tif 000010
【図９】
2019122502.tif 000011
【図１０】
2019122502.tif 000012
【図１１】
2019122502.tif 000013
【図１２】
2019122502.tif 000014
【図１３】
2019122502.tif 000015
【図１４】
2019122502.tif 000016
【図１５】
2019122502.tif 000017
【図１６】
2019122502.tif 000018
【図１７】
2019122502.tif 000019
【図１８】
2019122502.tif 000020
【図１９】
2019122502.tif 000021
【図２０】
2019122502.tif 000022
【図２１】
2019122502.tif 000023
【図２２】
2019122502.tif 000024
【図２３】
2019122502.tif 000025
【図２４】
2019122502.tif 000026