(57)【要約】【課題】心機能の健全性を評価する指標を、短時間且つ簡便な測定により取得する。【解決手段】被測定者の測定部位から脈波を検出するセンサ（３１）を有した脈波検出装置と、脈波検出装置とデータを遣り取りする装置（３０）とを備えるシステム（１００）であって、装置は、１２時間未満の予め定められた長さの時間においてセンサにより連続して検出される脈波のデータを時系列に収集し、収集される前記時系列の脈波データから、予め定められた周波数帯域のパワースペクトル密度を算出し、予め定められた周波数帯域におけるパワースペクトル密度を代表する値を、心機能の健全性を評価する指標として算出する。【選択図】図１
2019213820.tif 000002
【特許請求の範囲】
【請求項１】
被測定者の測定部位から脈波を検出するセンサを有した脈波検出装置と、
前記脈波検出装置とデータを遣り取りする装置とを備えるシステムであって、
前記装置は、
１２時間未満の予め定められた長さの時間において前記センサにより連続して検出されるデータを時系列に収集する収集手段と、
収集される前記時系列の脈波のデータから、予め定められた周波数帯域のパワースペクトル密度を算出する周波数解析手段を含み、
前記装置は、
前記予め定められた周波数帯域における前記パワースペクトル密度を代表する値を、心機能の健全性を評価する指標として算出する、システム。
【請求項２】
前記予め定められた長さの時間は、５分以上６０分以下の長さの時間を含む、請求項１に記載のシステム。
【請求項３】
前記周波数解析手段は、自己回帰移動平均モデルの演算に従い、収集された前記時系列の脈波のデータから、前記予め定められた周波数帯域のパワースペクトル密度を算出する、請求項１または２に記載のシステム。
【請求項４】
前記周波数解析手段は、収集された前記時系列の脈波のデータから、フーリエ変換により、前記予め定められた周波数帯域のパワースペクトル密度を算出する、請求項１または２に記載のシステム。
【請求項５】
前記予め定められた周波数帯域は、生体の圧受容器反射が血圧変動に作用する周波数帯域に相当する、請求項１から４のいずれか１項に記載のシステム。
【請求項６】
前記予め定められた周波数帯域は、０．０１Ｈｚ〜０．１Ｈｚの帯域を含む、請求項１から５のいずれか１項に記載のシステム。
【請求項７】
前記代表する値は、前記予め定められた周波数帯域における前記パワースペクトル密度の大きさを示す値を含む、請求項１から６のいずれか１項に記載のシステム。
【請求項８】
前記大きさを示す値は、前記予め定められた周波数帯域における前記パワースペクトル密度の積分値を含む、請求項７に記載のシステム。
【請求項９】
前記装置は、さらに、評価部を備え、
前記評価部は、
前記心機能の健全性を評価する指標が示す値を対応の閾値と比較し、比較の結果に基き、当該指標の値を前記心機能のリスクの段階のいずれかに分類する、請求項１から８のいずれか１項に記載のシステム。
【請求項１０】
前記装置は、さらに、
前記予め定められた周波数帯域におけるパワースペクトル密度の変位を、前記日内の血圧変動性を評価する指標として算出する、請求項９に記載のシステム。
【請求項１１】
前記評価部は、さらに、
前記日内の血圧変動性を評価する指標が示す値を対応の閾値と比較し、比較の結果に基き、当該指標の値を循環器のリスクの段階のいずれかに分類する、請求項１０に記載のシステム。
【請求項１２】
前記装置は、さらに、
収集される時系列の脈波のデータから血圧を算出する、請求項９から１１のいずれか１項に記載のシステム。
【請求項１３】
前記評価部は、さらに、
前記血圧を対応の閾値と比較し、比較の結果に基き、当該血圧を循環器のリスクの段階のいずれかに分類する、請求項１２に記載のシステム。
【請求項１４】
前記装置は、さらに、
前記リスクの段階を示す情報を、表示デバイスに表示させる、請求項９から１３のいずれか１項に記載のシステム。
【請求項１５】
前記装置は、さらに、
前記リスクの段階に関連付けられた治療方針を、表示デバイスに表示させる、請求項９から１４のいずれか１項に記載のシステム。
【請求項１６】
１２時間未満の予め定められた長さの時間において測定部位からセンサにより連続して検出される脈波のデータを時系列に収集する手段と、
収集される時系列の脈波のデータから、パワースペクトル密度を算出する周波数解析手段と、を備え、
予め定められた周波数帯域における前記パワースペクトル密度を代表する値を、心機能の健全性を評価する指標として算出する、装置。
【請求項１７】
プログラムであって、前記プログラムはコンピュータに、
１２時間未満の予め定められた長さの時間において測定部位からセンサにより連続して検出される脈波のデータを時系列に収集するステップと、
収集される時系列の脈波のデータから、パワースペクトル密度を算出するステップと、
予め定められた周波数帯域における前記パワースペクトル密度を代表する値を、心機能の健全性を評価する指標として算出するステップを実行させる、プログラム。
【請求項１８】
被測定者の測定部位から脈波を検出するセンサを有した脈波検出装置と、
前記脈波検出装置とデータを遣り取りする装置とを備えるシステムであって、
前記装置は、
１２時間未満の予め定められた長さの時間において前記センサにより連続して検出されるデータを時系列に収集する収集手段と、
収集される前記時系列の脈波のデータから、予め定められた周波数帯域の脈波のデータを抽出する手段を含み、
前記装置は、
抽出された前記脈波のデータが示す血圧のばらつきの大きさを、心機能の健全性を評価する指標として算出する、システム。
【請求項１９】
１２時間未満の予め定められた長さの時間において測定部位からセンサにより連続して検出される脈波のデータを時系列に収集する手段と、
収集される前記時系列の脈波のデータから、予め定められた周波数帯域の脈波のデータを抽出する手段を含み、
抽出された前記脈波のデータが示す血圧のばらつきの大きさを、心機能の健全性を評価する指標として算出する、装置。
【請求項２０】
プログラムであって、前記プログラムはコンピュータに、
１２時間未満の予め定められた長さの時間において測定部位からセンサにより連続して検出される脈波のデータを時系列に収集するステップと、
収集される前記時系列の脈波のデータから、予め定められた周波数帯域の脈波のデータを抽出するステップと、
抽出された前記脈波のデータが示す血圧のばらつきの大きさを、心機能の健全性を評価する指標として算出するステップと、を実行させる、プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この開示は、心機能の健全性の評価を支援するシステム、装置およびプログラムに関する。
【背景技術】
【０００２】
従来の高血圧診療では、カフを用いた血圧計等により診察室や家庭でスポット的に計測した血圧値を基準に、高血圧に関する評価がなされてきた。また、血圧は一拍ごとに変動し、その変動自体も循環器リスク指標の一つであることが知られている。したがって単回の血圧測定では血圧の変動性を評価できないばかりか、計測した血圧値は血圧変動の影響によりばらつきが生じる（測定頻度が血圧変動頻度よりも遅いために原理的に生じるエリアスと呼ばれるバラツキ）ため、高血圧の診断にはある程度長期間の血圧値を平均して評価する必要がある。つまり、患者の日常の行動（排泄、食事、喫煙、運動、会話、就寝等）による血圧の変動を考慮した評価が必要とされている。
【０００３】
このような背景のもと、ＡＢＰＭ（Guidelines for the clinical use of 24 hour ambulatory blood pressure monitoring）を用いて日内の行動による血圧変動を評価することが実施されている。ＡＢＰＭでは、１５〜３０分に１回の測定で２４時間に亘り血圧測定が実施される。ＡＢＰＭにより測定された日内の血圧値（収縮期血圧、拡張期血圧等）の標準偏差ＳＤ（standard deviation）または変動係数ＣＶ（Coefficient of Variation）から、血圧変動が大きいと評価された場合、心血管イベントのリスクは高いと判断される。
【０００４】
非特許文献１において、発明者らは、生体の圧受容器反射（baroreflex）が血圧変動を抑制するように作用することに着目している。すなわち、非特許文献１は、圧受容器反射の感度（baroreflex gain）が強力に作用する周波数帯域が０．０１〜０．１Ｈｚでは血圧の変動は、主に循環器のリスクを表すとの知見のもとで、１２時間に亘って測定された血圧波形のＰＳＤ（パワースペクトル密度：Power Spectral Density）の０．０１〜０．１Ｈｚにおける変位（ＰＳＤの傾き）は、圧受容体反射のゲインの算出に用いることができることを開示する。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】Hiroshi Mannoji, Keita Saku, Takuya Nishikawa, Takeshi Tohyama, Yasuhiro Oga, Kiyokazu Abe, Takuya Kishi, Hiroyuki Tsutsui, and Kenji Sunagawa “Identification of open-loop baroreflex gain from continuous arterial pressure monitoring”，[Online]，1 Apr 2017，Abstract Number:848.11，FASEB Journal, [2018年4月26日検索]，インターネット〈URL：https://www.fasebj.org/doi/abs/10.1096/fasebj.31.1_supplement.848.11〉
【発明の概要】
【発明が解決しようとする課題】
【０００６】
高血圧患者において、安静時の血圧、血圧変動および心機能を統合的に評価すること心臓または血管の機能不全を防止することが望まれている。しかし、心機能の検査は、エコー検査に代表されるように、血圧測定とは別の検査を実施する必要があった。また、エコー検査のための特別な装置が必要であり、検査に長時間を要していた。したがって、心機能の健全性を評価するための指標を、短時間且つ簡便な測定により、また、他の指標（血圧等）とともに取得できる技術の提供が望まれていた。
【０００７】
本開示は、心機能の健全性を評価する指標を、短時間且つ簡便な測定により、他の指標とともに取得できることを一つの目的としている。
【課題を解決するための手段】
【０００８】
本開示の一例に係るシステムは、被測定者の測定部位から脈波を検出するセンサを有した脈波検出装置と、脈波検出装置とデータを遣り取りする装置とを備えるシステムであって、装置は、１２時間未満の予め定められた長さの時間においてセンサにより連続して検出されるデータを時系列に収集する収集手段と、収集される時系列の脈波のデータから、予め定められた周波数帯域のパワースペクトル密度を算出する周波数解析手段を含み、装置は、予め定められた周波数帯域におけるパワースペクトル密度を代表する値を、心機能の健全性を評価する指標として算出する。
【０００９】
この開示によれば、測定部位から検出される脈波から、心機能の健全性を評価するための指標を算出することができるから、特別な装置（エコー検査装置等）を必要とせずに当該指標の取得が可能となる。また、測定部位から脈波を検出する期間を、１２時間未満の予め定められた長さの時間という、短い期間にすることができる。また、脈波が表す他の指標（例えば血圧）とともに、心機能の健全性を評価するための指標を取得することが可能になる。
【００１０】
上述の開示において、予め定められた長さの時間は、５分以上６０分以下の長さの時間を含む。
【００１１】
この開示によれば、予め定められた長さの時間を５分以上６０分以下にすることができる。
【００１２】
上述の開示において、周波数解析手段は、自己回帰移動平均モデルの演算に従い、収集された時系列の脈波のデータから、予め定められた周波数帯域のパワースペクトル密度を算出する。
【００１３】
この開示によれば、周波数解析に、自己回帰移動平均モデルを用いて、予め定められた周波数帯域のパワースペクトル密度を算出することができる。
【００１４】
上述の開示において、周波数解析手段は、収集された時系列の脈波のデータから、フーリエ変換により、予め定められた周波数帯域のパワースペクトル密度を算出する。
【００１５】
この開示によれば、周波数解析に、フーリエ変換を用いて、予め定められた周波数帯域のパワースペクトル密度を算出することができる。
【００１６】
上述の開示において、予め定められた周波数帯域は、生体の圧受容器反射が血圧変動に作用する周波数帯域に相当する。
【００１７】
この開示によれば、圧受容器反射が作用する周波数帯域におけるパワースペクトル密度から心機能の健全性を評価する指標を算出することができる。これにより、心機能の健全性を評価するための指標として、圧受容器反射を考慮した指標の算出が可能となる。
【００１８】
上述の開示において、予め定められた周波数帯域は、０．０１Ｈｚ〜０．１Ｈｚの帯域を含む。
【００１９】
この開示によれば、０．０１Ｈｚ〜０．１Ｈｚの帯域におけるパワースペクトル密度から心機能の健全性を評価する指標を算出することができる。
【００２０】
上述の開示において、代表する値は、予め定められた周波数帯域におけるパワースペクトル密度の大きさを示す値を含む。
【００２１】
この開示によれば、予め定められた周波数帯域におけるパワースペクトル密度の大きさを示す値を算出することで、心機能の健全性を評価する指標を得ることができる。
【００２２】
上述の開示において、代表する値は、予め定められた周波数帯域におけるパワースペクトル密度の積分値を含む。
【００２３】
この開示によれば、予め定められた周波数帯域におけるパワースペクトル密度の積分値を算出することで、心機能の健全性を評価する指標を得ることができる。
【００２４】
上述の開示において装置は、さらに、評価部を備え、評価部は、心機能の健全性を評価する指標が示す値を対応の閾値と比較し、比較の結果に基き、当該指標の値を心機能のリスクの段階のいずれかに分類する。
【００２５】
この開示によれば、心機能を評価する指標から、被測定者における心機能のリスクを段階的に分類することができる。
【００２６】
上述の開示において、装置は、さらに、予め定められた周波数帯域におけるパワースペクトル密度の変位を、日内の血圧変動性を評価する指標として算出する。
【００２７】
この開示によれば、予め定められた周波数帯域におけるパワースペクトル密度から、日内の血圧変動性を評価する指標も算出することができる。
【００２８】
上述の開示において、評価部は、さらに、日内の血圧変動性を評価する指標が示す値を対応の閾値と比較し、比較の結果に基き、当該指標の値を循環器のリスクの段階のいずれかに分類する。
【００２９】
この開示によれば、日内の血圧変動性を評価する指標から、被測定者における循環器のリスクを段階的に分類することができる。
【００３０】
上述の開示において、装置は、さらに、収集される時系列の脈波のデータから血圧を算出する。
【００３１】
この開示によれば、被測定者の血圧を算出することもできる。
上述の開示において、評価部は、さらに、血圧を対応の閾値と比較し、比較の結果に基き、当該血圧を循環器のリスクの段階のいずれかに分類する。
【００３２】
この開示によれば、血圧から、被測定者における循環器のリスクを段階的に分類することができる。
【００３３】
上述の開示において、装置は、さらに、リスクの段階を示す情報を、表示デバイスに表示させる。
【００３４】
この開示によれば、例えば、分類されたリスクの段階を視覚的に提示することができる。
【００３５】
上述の開示において、装置は、さらに、リスクの段階に関連付けられた治療方針を、表示デバイスに表示させる。
【００３６】
この開示によれば、例えば、被測定者の心機能または循環器についてのリスクの段階に関連付けられた治療方針を視覚的に提示することができる。
【００３７】
本開示の一例に係る装置は、１２時間未満の予め定められた長さの時間において測定部位からセンサにより連続して検出される脈波のデータを時系列に収集する手段と、収集される時系列の脈波のデータから、パワースペクトル密度を算出する周波数解析手段と、を備え、予め定められた周波数帯域における前記パワースペクトル密度を代表する値を、心機能の健全性を評価する指標として算出する。
【００３８】
この開示によれば、測定部位から検出される脈波から、心機能の健全性を評価するための指標を算出することができる。また、測定部位から脈波を検出する期間を、１２時間未満の予め定められた長さの時間という、短い期間にすることができる。また、脈波が表す他の指標（例えば血圧）とともに、心機能の健全性を評価するための指標を取得することが可能になる。
【００３９】
本開示の一例に係るプログラムはコンピュータに、１２時間未満の予め定められた長さの時間において測定部位からセンサにより連続して検出される脈波のデータを時系列に収集するステップと、収集される時系列の脈波のデータから、パワースペクトル密度を算出するステップと、予め定められた周波数帯域における前記パワースペクトル密度を代表する値を、心機能の健全性を評価する指標として算出するステップを実行させる。
【００４０】
この開示によれば、プログラムが実行されることにより、測定部位から検出される脈波から、心機能の健全性を評価するための指標を算出することができる。また、測定部位から脈波を検出する期間を、１２時間未満の予め定められた長さの時間という、短い期間にすることができる。また、脈波が表す他の指標（例えば血圧）とともに、心機能の健全性を評価するための指標を取得することが可能になる。
【００４１】
この開示に一例に係るシステムは、被測定者の測定部位から脈波を検出するセンサを有した脈波検出装置と、脈波検出装置とデータを遣り取りする装置とを備えるシステムであって、装置は、１２時間未満の予め定められた長さの時間においてセンサにより連続して検出されるデータを時系列に収集する収集手段と、収集される時系列の脈波のデータから、予め定められた周波数帯域の脈波のデータを抽出する手段を含み、装置は、抽出された脈波のデータが示す血圧のばらつきの大きさを、心機能の健全性を評価する指標として算出する。
【００４２】
この開示によれば、測定部位から検出される脈波から、心機能の健全性を評価するための指標を算出することができる。また、測定部位から脈波を検出する期間を、１２時間未満の予め定められた長さの時間という、短い期間にすることができる。また、脈波が表す他の指標（例えば血圧）とともに、心機能の健全性を評価するための指標を取得することが可能になる。
【００４３】
この開示の一例に係る装置は、１２時間未満の予め定められた長さの時間において測定部位からセンサにより連続して検出される脈波のデータを時系列に収集する手段と、収集される時系列の脈波のデータから、予め定められた周波数帯域の脈波のデータを抽出する手段を含み、抽出された脈波のデータが示す血圧のばらつきの大きさを、心機能の健全性を評価する指標として算出する。
【００４４】
この開示によれば、測定部位から検出される脈波から、心機能の健全性を評価するための指標を算出することができる。また、測定部位から脈波を検出する期間を、１２時間未満の予め定められた長さの時間という、短い期間にすることができる。また、脈波が表す他の指標（例えば血圧）とともに、心機能の健全性を評価するための指標を取得することが可能になる。
【００４５】
この開示の一例に係るプログラムはコンピュータに、１２時間未満の予め定められた長さの時間において測定部位からセンサにより連続して検出される脈波のデータを時系列に収集するステップと、収集される時系列の脈波のデータから、予め定められた周波数帯域の脈波のデータを抽出するステップと、抽出された脈波のデータが示す血圧のばらつきの大きさを、心機能の健全性を評価する指標として算出するステップと、を実行させる。
【００４６】
この開示によれば、プログラムが実行されると、測定部位から検出される脈波から、心機能の健全性を評価するための指標を算出することができる。また、測定部位から脈波を検出する期間を、１２時間未満の予め定められた長さの時間という、短い期間にすることができる。また、脈波が表す他の指標（例えば血圧）とともに、心機能の健全性を評価するための指標を取得することが可能になる。
【発明の効果】
【００４７】
本開示によれば、測定部位から検出される脈波から、心機能の健全性を評価するための指標を算出することができる。また、測定部位から脈波を検出する期間を、１２時間未満の予め定められた長さの時間という、短い期間にすることができる。また、脈波が表す他の指標（例えば血圧等）とともに、当該心機能の健全性を評価するための指標を取得することが可能になる。
【図面の簡単な説明】
【００４８】
【図１】本発明の実施の形態に係るシステムの構成の一例を模式的に示す図である。
【図２】本発明の実施の形態に係る血圧計１が左手首９０に装着された状態を模式的に示す図である。
【図３】本発明の実施の形態に係る血圧計の制御系のハードウェア構成の一例を示すブロック図である。
【図４】本発明の実施の形態に係る脈波計測時の圧力センサに対する押圧レベルの変化の一例を表すグラフである。
【図５】本発明の実施の形態に係る血圧変動性の評価に関する発明者らの知見を説明する図である。
【図６】本発明の実施の形態に係る血圧変動性の評価に関する発明者らの知見を説明する図である。
【図７】本発明の実施の形態に係るに全体的な処理の一例を示すフローチャートである。
【図８】本発明の実施の形態に係る算出と評価の処理の一例を示すフローチャートである。
【図９】本発明の実施の形態に係る自己回帰移動平均モデルを表す演算式の一例を示す図である。
【図１０】本発明の実施の形態に係る周波数解析処理の一例を示すフローチャートである。
【図１１】本発明の実施の形態に係る血圧の波形を拡大して模式的に示す図である。
【図１２】本発明の実施の形態に係るセグメント分割の一例を模式的に示す図である。
【図１３】本発明の実施の形態に係る血圧変動性を示す指標の算出の処理の一例を示すフローチャートである。
【図１４】本発明の実施の形態に係る周波数系列のＰＳＤを模式的に示す図である。
【図１５】本発明の実施の形態に係るΔＰＳＤとＳＤの相関性と所定長時間の関係性を表す実験結果を示す図である。
【図１６】本発明の実施の形態に係る心機能の健全性を示す指標の算出の処理の一例を示すフローチャートである。
【図１７】本発明の実施の形態に係る周波数系列のＰＳＤにおける指標ＨＩを模式的に示す図である。
【図１８】本発明の実施の形態に係る指標の評価基準の一例をテーブル形式で示す図である。
【図１９】本発明の実施の形態に係る表示画面の一例を模式的に示す図である。
【図２０】本発明の実施の形態に係る表示画面の他の例を模式的に示す図である。
【図２１】本発明の実施の形態に係る端末装置１０の構成の一例を示す図である。
【図２２】本発明の実施の形態に係る装置３０の構成の一例を示す図である。
【図２３】本発明の実施の形態に係る装置３０が備えるモジュールの一例を模式的に示す図である。
【図２４】本発明の実施の形態に係る心機能の健全性を示す指標の算出の処理の他の例を示すフローチャートである。
【図２５】図２４の処理における心機能の健全性を示す指標の算出を説明する図である。
【発明を実施するための形態】
【００４９】
本発明の実施の形態について、図面を参照しながら詳細に説明する。なお、図中の同一または相当部分については、同一符号を付してその説明は繰返さない。
【００５０】
＜Ａ．適用例＞
図１は、本発明の実施の形態に係るシステム１００の構成の一例を模式的に示す図である。図１を参照して、本発明が適用される場面の一例について説明する。図１を参照して、システム１００は、被測定者の測定部位から脈波を検出する圧力センサ３１を有した血圧計１と、血圧計１とデータを遣り取りする端末装置１０Ａ，１０Ｂおよび装置３０とを備える。装置３０は、例えば医家が操作するコンピュータであって、１２時間未満の予め定められた長さの時間において圧力センサ３１により連続して測定された時系列の脈波のデータを収集し、収集された脈波のデータを処理する。このデータ処理は、予め定められた長さの時間で連続して検出される脈波の時系列のデータに基づき、心機能の健全性を評価するための指標を算出する処理を含む。「心機能の健全性を評価するための指標」は、後述する心機能の健全性を推定できる指標ＨＩに相当する。
【００５１】
上記の予め定められた長さの時間は、発明者らが、実験に基づき得られた５分以上６０分以下の長さの時間を含む。
【００５２】
このように、心機能の健全性を推定できる指標を取得するために、従来は２４時間に亘って測定される脈波から算出する必要があったが、システム１００によれば、１２時間未満の予め定められた長さの時間で連続して測定される脈波を用いて当該指標を算出することができる。これにより、脈波測定のために被測定者を長時間拘束する必要がなくなる。５分以上６０分以下（例えば３０分程度）の期間の連続して検出される脈波を用いることで、本開示では、短時間の測定で指標を算出（推定）することが可能になる。
【００５３】
データ処理における指標算出では、収集された時系列の脈波のデータから、周波数解析により、予め定められた周波数帯域におけるパワースペクトル密度（ＰＳＤ）を算出する。そして、予め定められた周波数帯域におけるパワースペクトル密度の値を、心機能の健全性を評価するための指標として算出する。この値は、例えば、予め定められた周波数帯域におけるパワースペクトル密度の積分値を含む。
【００５４】
上記の予め定められた周波数は、生体の圧受容器反射が血圧変動に作用する周波数帯域である、０．０１〜０．１Ｈｚに相当する。したがって、圧受容器反射が作用する周波数帯域におけるパワースペクトル密度を用いて心機能の健全性を推定できる指標を算出することができるから、全ての帯域のパワースペクトル密度を算出する必要はなく、演算量を低減することができる。
【００５５】
装置３０は、心機能の健全性を評価する指標が示す値を対応の閾値と比較し、比較の結果に基き、当該指標の値を心機能のリスクの段階のいずれかに分類し、分類されたリスクの段階を表示する。また、分類されたリスクの段階に関連付けられた治療方針を表示する。これにより、例えば装置３０を操作する医家に対し、治療方針の判断を支援する情報を提供することができる。
【００５６】
装置３０は、さらに、他の指標を算出してもよい。例えば、予め定められた周波数帯域におけるパワースペクトル密度の変位を、日内の血圧変動性を評価する指標として算出する。または、例えば、予め定められた周波数帯域における脈波のデータから血圧（収縮期血圧、拡張期血圧等）を算出する。装置３０は、これら各指標の値も対応の閾値と比較し、比較の結果に基き、当該指標の値を循環器のリスクの段階のいずれかに分類し、分類されたリスクの段階を表示してもよい。さらに装置３０は、上記の各リスクの段階を組合わせて、その組合わせに関連付けられた治療方針を、表示してもよい。
【００５７】
このように、測定部位から検出される脈波から、心機能の健全性を評価するための指標を、脈波が表す他の指標（例えば血圧、血圧変動性等の指標）とともに取得することが可能になる。
【００５８】
上記のシステム１００では、医家の装置３０が、心機能の健全性を評価するための指標を算出するデータ処理を実施するが、データ処理は、端末装置１０Ａ，１０Ｂ、または血圧計１により実施してもよい。また、端末装置１０Ａ，１０Ｂは、被測定者が携帯するスマートフォン等の端末であってよい。例えば、血圧計１が測定した脈波のデータを、最寄りの端末装置１０Ａ，１０Ｂに転送し、端末装置１０Ａ，１０Ｂは受信した脈波のデータから、当該指標を算出する。これにより、被測定者は、血圧計１で血圧を測定しながら、携帯する端末装置１０Ａ，１０Ｂで、リアルタイムに心機能の健全性を評価するための指標を算出させることができる。
【００５９】
また、上記のシステム１００では、被測定者の測定部位から脈波を検出する装置は、血圧計１としているが、血圧波形に同期した脈波を検出する装置であれば、血圧計１に限定されない。例えば、容積脈波を検出する装置、例えば動脈に含まれる酸素の飽和度を測定するパルスオキシメーター（pulse oximeter）であってよい。その場合は、脈波として、パルスオキシメーターは光電式で測定する容積脈波であってよい。また、血圧計１は、図１では携帯型であるが、携帯型に限定されず据え置き型であってもよい。
【００６０】
以下、本開示のより具体的な応用例として、本実施の形態に係るシステム１００のより詳細な構成および処理について説明する。
【００６１】
＜Ｂ．システムの構成＞
再び図１を参照して、システム１００は、被測定者（ユーザともいう）が使用する端末装置１０Ａ，１０Ｂと、血圧計１と、装置３０と、ネットワーク４１，４３とを含む。以下では、端末装置１０Ａと端末装置１０Ｂを「端末装置１０」と総称する場合がある。
【００６２】
血圧計１は、携帯型であって、例えば本体とカフ（腕帯）とが一体となった腕時計型で提供されるが、腕時計型に限定されない携帯型であってもよい。本体とカフが別体であってもよい。また、携帯型に限定されず据え置き型であってもよい。
【００６３】
本実施の形態では、血圧計１は、腕時計のように長時間腕に装着して、２４時間連続して脈波を１拍ごとに測定する機能や、常時装着して測定開始ボタンを押すだけで測定できる機能などを有する。なお、血圧計１は、血圧とは異なる種類（歩数等）の生体情報を測定する機能を備えてもよい。
【００６４】
端末装置１０は、たとえば、タッチパネルを備えるスマートフォンである。以下では、スマートフォンを「端末装置１０」の代表例として説明を行なう。ただし、端末装置１０は、折り畳み式携帯電話、タブレット端末装置、ＰＣ（personal  computer）、ＰＤＡ（Personal  Data  Assistance）などのような他の装置であってもよい。
【００６５】
ネットワーク４１は、端末装置１０Ａと、端末装置１０Ｂと、装置３０とを互いに接続するために、インターネット、移動体端末通信網などの各種ネットワークを含む。端末装置１０Ｂと血圧計１とを接続するためのネットワーク４３は、典型的にはＵＳＢ（Universal Serial Bus）シリアル通信が採用され得るが、ここでは、近距離無線通信方式を採用しており、例えばＢＬＥ（Bluetooth（登録商標）  low  energy）による通信方式を採用する。他の通信方式としては、ＬＡＮ（Local Area Network）等が採用されてもよい。ネットワーク４３は、これらに限らず、各種の有線または無線の通信方式を採用してもよい。
【００６６】
装置３０は、端末装置１０から送信されるデータを受信し、受信したデータを格納し、また受信したデータを処理する。装置３０は、例えば、医家が操作可能なコンピュータに相当する。装置３０によるデータ処理の詳細は後述する。
【００６７】
血圧計１は、測定した血圧のデータを、端末装置１０を介して、装置３０に送信する。血圧計１が装置３０と通信する場合、端末装置１０は当該通信を中継する中継器に相当する。なお、血圧計１は、端末装置１０を介さずに、装置３０と直接通信してもよい。
【００６８】
＜Ｃ．血圧計の構成＞
図２は、本発明の実施の形態に係る血圧計１が左手首９０に装着された状態（以下、「装着状態」とも称する）を模式的に示す図である。図２では、左手首９０の長手方向に対して垂直な断面が模式的に示される。本実施の形態では、左手首９０が血圧の測定部位となる。なお、「測定部位」は、動脈が通っている部位であればよい。測定部位は、例えば、手首、上腕などの上肢であってもよいし、足首、大腿などの下肢であってもよい。
【００６９】
図２を参照して、血圧計１は、本体７０とベルト２０を備える。ベルト２０は、帯状体２３と、カフ２１とを含む。
【００７０】
カフ２１は、帯状体２３の内周面に沿って取り付けられる。カフ２１は、伸縮可能な２枚のポリウレタンシートからなる流体袋として構成されている。流体袋は、流体を収容可能な袋状の部材であればよい。「流体」は、液体と気体の両方を含み、例えば、水、空気などを用いることができる。
【００７１】
本体７０は、ベルト２０と一体に設けられる。なお、ベルト２０と本体７０とを別々に形成し、ベルト２０に対して本体７０を係合部材（例えば、ヒンジ）を介して、一体に取り付ける構成でもよい。本実施の形態では、本体７０が配置された部位は、装着状態において左手首９０の背側面（手の甲側の面）９０ｂに対応する。図２中には、左手首９０内で掌側面（手の平側の面）９０ａ近傍を通る橈骨動脈９１が示されている。
【００７２】
本体７０の頂面（測定部位から最も遠い側の面）には、後述のディスプレイ５０が設けられる。本体７０の側面は、ユーザからの指示を入力するための後述の操作部５２が設けられる。
【００７３】
帯状体２３は、厚さ方向に関して可撓性を有し、周方向（長手方向）に関して非伸縮性を有するプラスチック材料から構成される。
【００７４】
血圧計１（ベルト２０）は、全体として略環状に構成されるとともに、本体７０の底面とベルト２０の端部とが、図示しないバックル等によって図２中の矢印Ａ方向に開閉可能に構成される。
【００７５】
ユーザは、血圧計１を装着する場合、環状のベルト２０に左手を通す。装着時は、左手首９０の周りのベルト２０の角度位置を調節して、左手首９０を通る橈骨動脈９１上（測定部位）にカフ２１の圧力センサ３１を位置させる。これにより、圧力センサ３１は、左手首９０の掌側面９０ａのうち橈骨動脈９１に対応する部分に当接する状態となる。このようにして、ユーザは血圧計１（ベルト２０）を左手首９０に装着する。
【００７６】
（ｃ１．ハードウェア構成）
図３は、本発明の実施の形態に係る血圧計１の制御系のハードウェア構成の一例を示すブロック図である。図４を参照して、本体７０は、血圧計１を制御するためのＣＰＵ(Central Processing Unit)１００と、ディスプレイ５０と、メモリ５１と、操作部５２と、電池５３と、通信Ｉ／Ｆ（インターフェィスの略）５９と、外部の記憶媒体１０７を読み書きするＲ／Ｗ（Read/Writer）１０６を含む。また、本体７０は、さらに、ポンプ３２と、弁３３とを含む。
【００７７】
さらに、本体７０は、各部に電力を供給する電池５３、圧力センサ３１からの出力を周波数に変換する発振回路３１０と、ポンプ３２を駆動するポンプ駆動回路３２０とを含む。圧力センサ３１は、カフ２１に一体的に備えられる。
【００７８】
ディスプレイ５０は、例えば、有機ＥＬ(Electro Luminescence)ディスプレイで構成され、ＣＰＵ１００からの制御信号に従って、血圧測定結果などの血圧に関する情報、その他の情報を表示する。
【００７９】
操作部５２は、例えば、プッシュ式スイッチで構成され、ユーザによる血圧測定開始または停止の指示に応じた操作信号をＣＰＵ１００に出力する。
【００８０】
メモリ５１は、血圧計１を制御するためのプログラム、血圧計１を制御するために用いられるデータ、血圧測定の結果のデータ等を記憶する。また、メモリ５１は、プログラムが実行されるときのワークメモリ等として用いられる。
【００８１】
ＣＰＵ１００は、メモリ５１に記憶された制御プログラムを実行することにより、例えば、トノメトリ法による血圧測定を実行する場合、ＣＰＵ１００は、操作部５２からの血圧測定開始の指示に応じて、圧力センサ３１からの信号に基づいて、ポンプ３２（および弁３３）を駆動する制御を行なう。
【００８２】
ＣＰＵ１００は、トノメトリ法による血圧測定を実行する場合、操作部５２からの血圧測定開始の指示に応じて、カフ２１内の空気を排出させるために弁３３を駆動する制御を行なう。また、ＣＰＵ１００は、圧力センサ３１からの信号に基づいて、血圧波形を導出する制御を行なう。
【００８３】
通信Ｉ／Ｆ５９は、ネットワーク４３を介して通信するための送受信回路を含む。
ポンプ３２および弁３３は、エア配管３９ａを介してカフ２１に接続されている。圧力センサ３１はカフ２１に接続されて、カフ２１内の圧力を検出する。
【００８４】
ポンプ３２は、例えば、圧電ポンプで構成される。ポンプ３２は、カフ２１内の圧力（カフ圧）を加圧するために、エア配管３９ａを通してカフ２１に加圧用の流体としての空気を供給する。
【００８５】
弁３３は、ポンプ３２に搭載され、ポンプ３２のオン／オフに伴って開閉が制御されるように構成される。具体的には、弁３３は、ポンプ３２がオンされると閉じて、カフ２１内に空気を封入する一方、ポンプ３２がオフされると開いて、カフ２１の空気をエア配管３９ａを通して大気中へ排出させる。
【００８６】
弁３３は、逆止弁としても動作し、排出されるエアが逆流することはない。ポンプ駆動回路３２０は、ポンプ３２をＣＰＵ１００から与えられる制御信号に基づいて駆動する。
【００８７】
圧力センサ３１は、例えば、ピエゾ抵抗式圧力センサであり、カフ２１に接続されている。圧力センサ３１は、ベルト２０（カフ２１）の圧力、例えば、大気圧を基準（ゼロ）とした圧力を検出して時系列の信号として出力する。
【００８８】
発振回路３１０は、圧力センサ３１からのピエゾ抵抗効果による電気抵抗の変化に基づく電気信号値に応じた周波数を有する周波数信号をＣＰＵ１００に出力する。圧力センサ３１の出力は、カフ２１内の圧力を制御するため、および、血圧値（収縮期血圧（ＳＹＳ：Systolic Blood Pressure）と拡張期血圧（ＤＩＡ：Diastolic Blood Pressure）等を含む。）を算出するために用いられる。
【００８９】
＜Ｄ．脈波の検出処理＞
本発明の実施の形態では、血圧測定方法として、例えばトノメトリ法を採用する。トノメトリ法は、血管を皮膚上から圧迫し、圧脈波（以下、脈波という）を非侵襲的に計測できるとのメリットを有する。なお、脈波を検出する方法はトノメトリ法に限定されない。
【００９０】
図４は、本発明の実施の形態に係る脈波計測時の圧力センサ３１に対する押圧レベルの変化の一例を表すグラフである。該グラフの縦軸に押圧の圧力値（単位はｍｍＨｇ）が採られて、横軸に時間経過が採られている。図４を参照して、時間Ｔ１にて計測開始されると動脈９１直上に設置された圧力センサ３１はカフ２１の内圧の上昇によって生体に押圧される。カフ２１の内圧が徐々に高まるにつれて、圧力センサ３１は動脈を橈骨に押しつける状態となり、動脈は押しつぶされていく。圧力センサ３１が最適圧で生体に押圧されているとき、動脈は扁平した状態（トノメトリ状態）となり、圧力センサ３１は動脈内圧の変動に起因した脈波を正確に計測することができる。さらにカフ２１による押圧レベルを上げると、動脈はさらに押しつぶされ、最後に閉塞する。以上の過程にておいては圧力センサ３１および発振回路３１０により検出される周波数信号から脈波の波形変化を検出し、検出した波形変化に基づきＣＰＵ１００が最適圧を判定している。計測では、最適圧以上まで押圧レベルを上昇させたときの時間Ｔ２から、押圧レベルは最適圧まで減圧され、ＣＰＵ１００が最適圧にまで達したと判断した時の時間Ｔ３以降は、最適圧が保持される。最適圧を保持している間に圧力センサ３１により検出される脈波は発振回路３１０により周波数信号に変換されて、ＣＰＵ１００に出力される。ＣＰＵ１００は、周波数信号を血圧の波形を示す脈波データに変換（校正）しメモリ５１に格納する、または脈波データを通信Ｉ／Ｆ５９を介して時系列に装置３０に送信する。
【００９１】
＜Ｅ．知見＞
図５および図６は、本発明の実施の形態に係る血圧変動性の評価に関する発明者らの知見を説明する図である。発明者らは、実験により、被測定者が安静時の１２時間未満の所定長時間で連続して測定される脈波（血圧波形）を、周波数解析して得られる、圧受容器反射（baroreflex）が作用する０．０１〜０．１ＨｚのΔＰＳＤ（すなわちＰＳＤの変位（傾き）の大きさ）は、被測定者の血圧をＡＢＰＭ等により２４時間測定して得られる日内（２４時間）の血圧の標準偏差ＳＤ（standard deviation）、すなわち日内の血圧変動と、高い相関性を有することを発見した（図５参照）。図５は、発明者らの実験結果を示しており、ΔＰＳＤが大きいほど日内の血圧の標準偏差ＳＤ（変動）が大きく、ΔＰＳＤが小さいほど日内の血圧の標準偏差ＳＤ（変動）が小さいことを示している。このように、０．０１〜０．１ＨｚのΔＰＳＤは、日内の血圧変動性を推定することができる指標として使えるとの知見を得た。
【００９２】
さらに、発明者らは、上記の０．０１〜０．１ＨｚにおけるＰＳＤから、被測定者の心機能を評価することができるとの知見を得た。すなわち、血圧は心拍出と血管抵抗から成る。心機能が低下すると心拍出の変化量は低下し、圧受容器反射に拘わらず、血圧は変動できなくなる。発明者らは、この現象が０．０１〜０．１ＨｚにおけるＰＳＤの低下として現れることを発見した（図６参照）。このように発明者らは、血圧の波形（すなわち脈波）を周波数解析して得られるＰＳＤの値は心機能を評価することができる指標として使えるとの知見も得た。
【００９３】
発明者らは、これらの知見に基づき、本明細書で開示する血圧変動性または心機能の健全性を評価する技術を提案する。
【００９４】
＜Ｆ．装置の構成＞
図２１は、本発明の実施の形態に係る端末装置１０の構成の一例を示す図である。図２２は、本発明の実施の形態に係る装置３０の構成の一例を示す図である。図２１を参照して、端末装置１０は、ＣＰＵ１５２、タイマ１５３、メモリ１５４、端末装置１０に対するユーザの操作を受付ける操作部１５６、ディスプレイ１５８、アンテナ１６２を介した無線のための通信部１６０、メモリインターフェイス（Ｉ／Ｆ）１６４、通信インターフェイス（Ｉ／Ｆ）１６６、音声出力のためのスピーカ１６８および音声入力のためのマイク１７０とを含む。
【００９５】
ＣＰＵ１５２は、メモリ１５４に記憶されたプログラムを実行することにより、各部を制御する。
【００９６】
メモリ１５４の記憶領域は、ＲＡＭ（Random  Access  Memory）、ＲＯＭ（Read-Only  Memory）、フラッシュメモリ、ハードディスク装置などにより構成される。メモリ１５４には、ＣＰＵ１５２によって実行されるプログラム、またはＣＰＵ１５２によって用いられるデータなどが格納される。
【００９７】
操作部１５６は、端末装置１０に対する操作入力を受付ける。典型的には、操作部１５６は、タッチパネルを含んで実現される。タッチパネルは、ディスプレイ１５８上に設けられる。また、操作部１５６は、スイッチ・ボタンなどを含んでもよい。
【００９８】
通信部１６０は、アンテナ１６２を介して移動体通信網に接続し無線通信のための信号を送受信する。これにより、端末装置１０は、他の通信装置（たとえば、装置３０、他の端末装置１０）との通信が可能となる。
【００９９】
メモリインターフェイス１６４は、外部の記憶媒体１６５からデータを読出す。ＣＰＵ１５２は、メモリインターフェイス１６４を介して記憶媒体１６５に格納されているデータを読出して、当該データをメモリ１５４に格納する。ＣＰＵ１５２は、メモリ１５４からデータを読出して、メモリインターフェイス１６４を介して当該データを外部の記憶媒体１６５に格納する。
【０１００】
記憶媒体１６５は、ＣＤ（Compact  Disc）、ＤＶＤ（Digital  Versatile  Disk）、ＢＤ（Blu-ray（登録商標）  Disc）、ＵＳＢ（Universal  Serial  Bus）メモリ、ＳＤ（Secure  Digital）メモリカードなどの不揮発的にプログラムまたはデータを格納する媒体を含む。
【０１０１】
通信Ｉ／Ｆ（インターフェイスの略）１６６は、血圧計１との間でデータをやり取りするための通信インターフェイスであり、アダプタやコネクタなどによって実現される。
【０１０２】
図２２を参照して、装置３０は、ＣＰＵ５０２を含むコントローラ５０１、情報を出力するための出力デバイス５０３、および装置３０に対するユーザ操作を受付けるためのボタン・スイッチなどを含む操作デバイス５０４を備える。装置３０は、さらに、ネットワーク４１を介して端末装置１０などと通信するための通信Ｉ／Ｆ５０６、プログラムおよびデータを格納するＲＯＭ、ＲＡＭなどを含む記憶媒体であるストレージ５０５、データベース５１０などの各種データを格納するためのＨＤＤ（Hard  Disc  Drive）５０７、外部の記憶媒体５０９が脱着自在に装着されるドライバ５０８を備える。
【０１０３】
出力デバイス５０３は、ディスプレイ５１１、プリンタ等を含む。ドライバ５０８は、装着された記憶媒体５０９からデータまたはプログラムを読出し、また、装着された記憶媒体５０９にデータを格納する。記憶媒体５０９から読出されたデータまたはプログラムは、ストレージ５０５に格納される。記憶媒体５０９は、ＣＤ（Compact  Disc）、ＤＶＤ（Digital  Versatile  Disk）、ＢＤ（Blu-ray（登録商標）  Disc）、ＵＳＢ（Universal  Serial  Bus）メモリ、ＳＤ（Secure  Digital）メモリカードなどの不揮発的にプログラムまたはデータを格納する媒体を含む。
【０１０４】
＜Ｇ．機能構成＞
図２３は、本発明の実施の形態に係る装置３０が備えるモジュールの一例を模式的に示す図である。図２３を参照して、装置３０は、血圧計１から受信することにより脈波データを収集する収集部６００、血圧の変動を示す指標を算出するために、脈波データから各種指標を算出する指標算出部６０１、血圧を算出する血圧算出部６０４、および指標を評価する評価部６０５を備える。指標算出部６０１は血圧計１から受信した脈波データを用いた周波数解析を実施する周波数解析部６０３を含む。
【０１０５】
モジュールがプログラムである場合、ＣＰＵ５０２がストレージ５０５または記憶媒体５０９等からプログラムを読出し、ＣＰＵ５０２が読出されたプログラムを実行することにより、収集部６００、指標算出部６０１および評価部６０５が実現される。なお、収集部６００または指標算出部６０１または評価部６０５は、プログラムモジュールから構成されるものに限定されず、ＦＰＧＡ（Field-Programmable Gate Array）またはＡＳＩＣ（Application Specific Integrated Circuit）等の回路モジュールを含んで実現されてもよい。
【０１０６】
＜Ｈ．全体フローチャート＞
図７は、本発明の実施の形態に係るに全体的な処理の一例を示すフローチャートである。図７では、血圧計１と端末装置１０間の通信および端末装置１０と装置３０間の通信のシーケンスも示される。図７に示すように、本実施の形態では、血圧計１と装置３０との両者間の通信は、端末装置１０を中継するとしている。なお、両者は中継器を介さずに直接通信するとしてもよい。
【０１０７】
図７を参照して、まず被測定者は血圧計１を測定部位に装着する。血圧計１のＣＰＵ１００は、操作部５２等から血圧測定開始の指示を受付けると、図４に従うカフ圧の変化がなされるように各部を制御する。本実施の形態では、被測定者は、安静状態で血圧測定を実施する。
【０１０８】
カフ圧の制御により図４の時間Ｔ３以降の最適圧が保持される期間において、ＣＰＵ１００は、測定部位における圧力センサ３１により１拍毎に連続して検出される脈波を表す周波数信号を、発振回路３１０から受付けて、受付けた周波数信号を時系列の脈波データに変換しメモリ５１に格納する。本実施の形態では、この周波数信号は、圧力センサ３１からの出力を例えば５ミリ秒毎にサンプリングした信号からなる。
【０１０９】
ＣＰＵ１００は、上記の最適圧が保持される期間における予め定められた長さの時間（以下、所定長時間という）内において周波数信号を受付けて、周波数信号をデジタルデータに変換することで脈波データを収集する（ステップＳ１）。ＣＰＵ１００は、脈波データを収集開始してから所定長時間が経過したか否かに基づき、収集が完了したか否かを判断する（ステップＳ３）。ＣＰＵ１００が、収集は完了していないと判断すると（ステップＳ３でＮＯ）、ステップＳ１の処理が継続するが、収集は完了したと判断すると（ステップＳ３でＹＥＳ）、送信処理（ステップＳ４）を実施する。
【０１１０】
本実施の形態では、所定長時間は、１２時間未満の長さの時間であり、より典型的には、５分以上（例えば、３０分等）６０分以下の時間を含み得る。
【０１１１】
ステップＳ４では、ＣＰＵ１００は、メモリ５１から所定長時間の長さに相当する時系列の脈波データを読出し、通信Ｉ／Ｆ５９を介して、端末装置１０を中継して、装置３０に送信する（ステップＳ４およびＳ５）。装置３０の収集部６００は、血圧計１から所定長時間の時系列の脈波データを受信することにより、脈波データを収集する（ステップＳ６）。なお、中継用の端末装置１０は血圧計１の最寄りに位置する端末であって、典型的には、被測定者が所持する端末装置１０を含む。
【０１１２】
本実施の形態では、端末装置１０は、所定長時間の脈波データの収集が完了したときに、収集された脈波データを装置３０に送信開始するが、送信のタイミングは脈波データの収集完了時に限定されない。例えば、血圧計１は、脈波の測定と、脈波データを装置３０に送信する処理とを並行して実施してもよい。
【０１１３】
装置３０のＣＰＵ５０２は、受信した所定長時間の脈波データを時系列にメモリ５１等に格納する。ＣＰＵ５０２は、所定長時間の脈波データを処理する。具体的には、指標算出部６０１は、所定長時間の脈波データから、日内（２４時間）の血圧変動性を推定できる指標と心機能の健全性を推定できる指標とを算出し、安静時の血圧を算出し、評価部６０５は指標等を評価する（ステップＳ７）。
【０１１４】
装置３０のＣＰＵ５０２は、ステップＳ７の処理の結果を、端末装置１０を中継し、血圧計１に送信する（ステップＳ８）。また、ＣＰＵ５０２は、ステップＳ７の処理の結果を、出力デバイス５０３を介し表示，印字等し、また、被測定者に関連付けてデータベース５１０に格納する（ステップＳ１３）。
【０１１５】
端末装置１０のＣＰＵ１５２は、装置３０から処理の結果を受信すると、受信データを中継するか否かを判断する（ステップＳ９）。この判断は、例えばメモリ１５４に設定されたフラグに基づきなされる。
【０１１６】
ＣＰＵ１５２は、中継すると判断すると（ステップＳ９でＹＥＳ）、装置３０からの処理の結果をそのまま血圧計１に送信する。血圧計１のＣＰＵ１００は、端末装置１０から処理の結果を受信すると（ステップＳ１０）、受信した処理の結果を、ディスプレイ５０を介して表示し、また、メモリ５１に格納する（ステップＳ１２）。
【０１１７】
一方、ＣＰＵ１５２は、中継しないと判断すると（ステップＳ９でＮＯ）、装置３０から受信した処理の結果を、ディスプレイ１５８を介し表示し、また、メモリ１５４等に格納する（ステップＳ１１）。
【０１１８】
以上により、血圧計１で所定長時間の脈波が測定された場合、被測定者は、日内の血圧変動性の指標と、指標を用いた評価とが提示される。また、血圧計１が脈波の測定と、脈波データの装置３０への送信とを並行して実施する場合には、脈波の測定時に、リアルタイムに上記の指標と評価の算出がなされて、その結果が被測定者に提示される。
【０１１９】
上記のステップＳ９で用いるフラグは、例えば、血圧計１の機種に従い設定され得る。具体的には、機種によっては血圧計１のディスプレイ５０の画面サイズが制限される。その場合は、中継を許可しない（中継を実施しない）値をフラグに設定する。これにより、画面サイズがより大きい端末装置１０のディスプレイ１５８に処理の結果を表示させることができる。なお、フラグの設定方法は、血圧計１の機種に基づく方法に限定されない。
【０１２０】
＜Ｉ．指標等の算出＞
図８は、本発明の実施の形態に係る算出と評価の処理の一例を示すフローチャートである。図８を参照して、図７のステップＳ７における算出と評価の処理を説明する。まず、指標算出部６０１は、所定長時間の脈波データを用いて日内（２４時間）の血圧変動性を評価するための指標としてΔＰＳＤと、心機能の健全性を推定できる指標ＨＩとして後述する積分値を算出する（ステップＳ２５、Ｓ２７、Ｓ２９）。
【０１２１】
（ｉ１．指標の算出）
指標ΔＰＳＤと積分値の算出を説明する。
【０１２２】
（ｉ１-１．ＰＳＤの算出例）
図９は、本発明の実施の形態に係る自己回帰移動平均モデルを表す演算式の一例を示す図である。図８のステップＳ２５において、周波数解析部６０３は、ΔＰＳＤを算出するために、所定時間長の脈波データである時系列データを、例えば、図９の自己回帰移動平均（ＡＲＭＡ：Autoregressive Moving Average)モデルを用いて解析し、解析によりＰＳＤを算出する。
【０１２３】
図９を参照して、自己回帰移動平均モデルＦ１は、血圧変動の時系列ｙｔを、時点tの値ｙtが過去のｙの値ｙt-ｉ（i=１,２，・・・p）と過去のホワイト・ノイズεｔを用いた関数として表されるモデルである。また、自己回帰移動平均モデルＦ１の定数Ｃ、およびパラメータγiとパラメータθiは、血圧変動の制御を行なう生体の生理学的機構に依存した時定数等の値を示す。
【０１２４】
（ｉ１-２．ＰＳＤの他の算出例）
ＰＳＤの他の算出例を説明する。図１０は、本発明の実施の形態に係る周波数解析処理の一例を示すフローチャートである。図１１は、本発明の実施の形態に係る血圧の波形を拡大して模式的に示す図である。図１２は、本発明の実施の形態に係るセグメント分割の一例を模式的に示す図である。
【０１２５】
図８のステップＳ２５において周波数解析部６０３は、図１０の処理に従い、所定時間長の脈波データ（時系列データ）から、フーリエ変換を用いて、各周波数帯域のＰＳＤを算出する。この脈波データは、図１１の示すように拍動の大きさに相当する振幅を有した血圧の波形を示す。
【０１２６】
図１０を参照して、フーリエ変換を適用するために、周波数解析部６０３は、所定時間長の脈波データ（例えば、図１１の血圧波形のデータ）を、予め定められた時間長さのセグメント単位に分割する（ステップＴ１）。この分割は、例えば図１２に示されるように、各セグメントがｎ分間の長さを有し、且つ隣接するセグメントどうしがｍ分間の重複を許容するようにして実施される。例えば、ｎ分は２分であり、ｍ分は１分である。
【０１２７】
周波数解析部６０３は、各セグメントの脈波データに前処理を施す（ステップＴ３）。具体的には、周波数解析部６０３は、各セグメントの脈波データに予め定められた関数を適用することにより、脈波データから、線形トレンドの成分を除去する線形デトレンド処理を実施し、その後、フーリエ変換におけるデータ端での不連続性を回避するため、各セグメントの脈波データに予め定められた窓関数を適用する。このような線形デトレンド処理により、測定された脈波データから、評価対象よりも十分に遅い周波数成分の影響が除去される。
【０１２８】
次に、周波数解析部６０３は、前処理がなされた各セグメントの脈波データに基づきＰＳＤを算出（推定）する（ステップＴ５）。具体的には、周波数解析部６０３は、各セグメントに高速フーリエ変換を実施する。これにより、各セグメントの脈波データから、周波数系列のデータを取得する。周波数解析部６０３は、全てのセグメントの周波数系列のデータを加算し、平均化することにより、ＰＳＤを推定（算出）する。
【０１２９】
なお、本実施の形態では、ＰＳＤの推定の方法として、Welch法やマルチテーパー法などのノンパラメトリック法でもよいし、自己回帰モデルまたはＡＲＭＡモデルといったパラメトリック法でもよい。
【０１３０】
（ｉ１-３．ΔＰＳＤの算出）
図１３は、本発明の実施の形態に係る血圧変動性を示す指標の算出の処理の一例を示すフローチャートである。図１４は、本発明の実施の形態に係る周波数系列のＰＳＤを模式的に示す図である。図１４のグラフでは、横軸の目盛は周波数（Ｈｚ）のオーダー（１０のべき乗）を取り、縦軸の目盛はＰＳＤのオーダー（１０のべき乗）を取る。図１５は、本発明の実施の形態に係るΔＰＳＤとＳＤの相関性と所定長時間の関係性を表す実験結果を示す図である。
【０１３１】
図８のステップＳ２７において指標算出部６０１は、ΔＰＳＤを算出する。具体的には、指標算出部６０１は、上記の算出された周波数系列のＰＳＤ（例えば、図１４参照）のうちから、圧受容器反射が作用する周波数帯域０．１Ｈｚ〜０．０１ＨｚのＰＳＤを抽出する。そして、指標算出部６０１は、抽出した帯域成分に基づき、０．１ＨｚのＰＳＤに対する０．０１ＨｚのＰＳＤの比であるΔＰＳＤを算出する（図１３のステップＲ１）。ΔＰＳＤは、後述の図１５に示すように、日内の血圧の標準偏差ＳＤと相関するから、ΔＰＳＤは日内の血圧変動性を推定できる指標として用いることができる。
【０１３２】
なお、日内の血圧変動性を推定できる指標は、周波数帯域０．１Ｈｚ〜０．０１ＨｚのＰＳＤの変位を表す値であればよく、上記のΔＰＳＤに限定されない。例えば、当該変位は、図１４に示すような、ＰＳＤと周波数に対応のＸ軸，Ｙ軸で規定される平面上において、０．１Ｈｚのパワースペクトル密度（ＰＳＤ）と０．０１Ｈｚのパワースペクトル密度（ＰＳＤ）の両方の値を結ぶ直線Ｌの傾きの大きさであってもよい。
【０１３３】
（ｉ１-４．指標ＨＩの算出例）
図１６は、本発明の実施の形態に係る心機能の健全性の指標の算出の処理の一例を示すフローチャートである。図１７は、本発明の実施の形態に係る周波数系列のＰＳＤにおける積分値を模式的に示す図である。
【０１３４】
図８のステップＳ２９において指標算出部６０１は、心機能の健全性の指標ＨＩを算出する。具体的には、指標算出部６０１は、上記の算出された周波数系列のＰＳＤ（例えば、図１７参照）のうちから、圧受容器反射が作用する周波数帯域０．１Ｈｚ〜０．０１Ｈｚに着目し、０．１Ｈｚ〜０．０１Ｈｚの帯域のＰＳＤを積分することにより積分値を算出する（図１６のステップＰ１）。
【０１３５】
図１７の積分値は、被測定者の心機能の健全性を推定できる指標ＨＩとして有効である。具体的には、ＰＳＤは、心臓の拍動の大きさを表すから、心臓の血液を送り出すポンプとしての機能が低下するとＰＳＤは小さくなる。また、圧受容器反射が強く作用する周波数帯域（０．１Ｈｚ〜０．０１Ｈｚの帯域）は、心機能の健全性を把握し易い帯域でもある。したがって、０．１Ｈｚ〜０．０１Ｈｚの帯域におけるＰＳＤの積分値は、心機能の健全性を推定する指標ＨＩとして有効である。なお、心機能の健全性を推定する指標ＨＩは、０．１Ｈｚ〜０．０１Ｈｚの帯域におけるＰＳＤの大きさを表す値であればよく、積分値に限定されず、例えば、平均値、中央値等であってもよい。
【０１３６】
（ｉ１-５．指標ＨＩの他の算出例）
図２４は、本発明の実施の形態に係る心機能の健全性を示す指標ＨＩの算出の処理の他の例を示すフローチャートである。図２５は、図２４の処理における心機能の健全性を示す指標ＨＩの算出を説明する図である。図２５では、血圧計１により測定される時系列の脈波データ２５１と、脈波データ２５１にバンドパスフィルタ処理を施して得られる血圧の波形を示すデータＡＰを模式的に示す。
【０１３７】
図２４と図２５を参照して、指標ＨＩの他の算出例を説明する。指標ＨＩの算出は、上記に述べたＰＳＤの圧受容器反射が作用する周波数帯域（０．１Ｈｚ〜０．０１Ｈｚの帯域）の積分値を算出する方法に限定されない。例えば、指数算出部６０１は、図２４に示すように、血圧計１から収集した時系列の脈波データ（例えば、図２５の脈波データ２５１）を、バンドパスフィルタ処理を施して、血圧の波形を示すデータＡＰを抽出する（ステップＰ２）。このバンドパスフィルタ処理の特性は、圧受容器反射が作用する周波数帯域に相当する脈波データ２５２のみを通過させるようなカットオフ周波数を有する。
【０１３８】
指標算出部６０１は、データＡＰが示す血圧のばらつきの大きさ（例えば、分散、標準偏差、２乗平均平方根等）を指標ＨＩとして算出する（ステップＰ３）。
【０１３９】
図２４の処理によれば、周波数領域での処理は必要ないので、指標ＨＩを簡単に算出することができる。
【０１４０】
（ｉ１-６．安静時血圧の算出）
血圧算出部６０４は、所定長時間の脈波データが示す血圧波形から、安静時血圧を算出する（図８のステップＳ３１）。この安静時血圧は、絶対的な血圧値を示す。具体的には、血圧算出部６０４は、所定長時間の脈波データから脈波包絡線を示すデータを生成し、この脈波包絡線データから血圧を算出する。このような血圧の算出は、公知の方法に従うので、説明は繰り返さない。
【０１４１】
＜Ｊ．ΔＰＳＤとＳＤ間の相関性と所定長時間に関する実験＞
発明者らは、実験により、ΔＰＳＤとＳＤ間に相関性があることと、相関性と所定長時間との間の関係性を確認した。図１５は、ラットにカテーテルを挿入し逐次血圧を測定して得られた実験結果を示すグラフである。図１５のグラフの縦軸には相関係数Ｒをとり、横軸には経過時間をとっている。相関係数Ｒは、実験により所定長時間で収集した血圧波形のデータからＡＲＭＡＳモデルＦ１で算出した血圧変動性を示す指標(ΔＰＳＤ)と、実験により測定した２４時間（日内）の平均血圧の標準偏差ＳＤとの相関の程度を表す。ただし、相関係数Ｒは、最大１．０である。
【０１４２】
図１５の実験結果によれば、血圧変動性の指標として所定長時間を０分から６０分の間で変化させた場合に、所定長時間(６０分以下)の測定値によるΔＰＳＤは、日内の平均血圧の標準偏差ＳＤと高い相関性を有することが確認できた。さらに、所定長時間は、５分以上であれば、相関の程度を示す係数Ｒは０．８以上と高い値を示すことも確認できた。このような実験結果から、発明者らは、図７のステップＳ１において、血圧を測定する所定長時間は、５分以上６０分以下であればよいとの知見を得た。
【０１４３】
＜Ｋ．評価＞
図１８は、本発明の実施の形態に係る指標の評価基準の一例をテーブル形式で示す図である。図１８を参照して、テーブル３００は、例えば装置３０のストレージ５０５またはＨＤＤ５０７に格納される。テーブル３００には、測定結果の組３０１のそれぞれに関連付けて心機能または循環器の治療方針を示すデータ３０２が登録されている。測定結果に組３０１は、血圧の絶対値ＢＰ、ΔＰＳＤおよび指標ＨＩとからなる組を示す。
【０１４４】
図８のステップＳ３３において、評価部６０５は、心機能（心血管）または循環器のリスクを評価する場合に、指標算出部６０１から出力されるΔＰＳＤおよび指標ＨＩ、および血圧算出部６０４から出力される血圧の絶対値ＢＰの各値を対応の閾値と比較し、比較の結果に基づき当該指標の値を心機能または循環器のリスクの段階のいずれかに分類する。具体的には、血圧の絶対値ＢＰについて「高い」か否か、日内の血圧変動性ΔＰＳＤについて「正常」または「高い」、および心機能の評価の指標ＨＩについて「正常」または「低い」とリスクの段階のいずれかに分類することにより評価する。
【０１４５】
また、評価部６０５は測定結果の組３０１の評価の組合せ（「高い」、「正常」、「正常」）、（「高い」、「高い」、「正常」）および（「高い」、「高い」、「低い」）に基づき、テーブル３００から関連付けされた治療方針のデータ３０２を読出す。このように、評価部６０５は、評価結果として、血圧の絶対値ＢＰ、日内の血圧変動性ΔＰＳＤ、および心機能の評価の指標ＨＩと、治療方針のデータ３０２を出力する。データ３０２が示す治療方針は、ＨＨＳ(Heart Health Score：心臓健康指数)の評価を定期的に受けることを推奨するメッセージを含んでいる。
【０１４６】
治療方針のデータ３０２は、例えば図７のステップＳ１３において、ディスプレイ５１１に表示される。これにより、医家に対し、治療方針の判断を支援する情報を提供することができる。
【０１４７】
＜Ｌ．表示画面例＞
図１９は、本発明の実施の形態に係る表示画面の一例を模式的に示す図である。図１９の画面は血圧波形の解析によるＨＨＳに関する情報を表示する。図１９を参照して、例えば装置３０のＣＰＵ５０２はディスプレイ５１１に、ＨＨＳ１９１、個別指標１９２、被測定者の血圧波形１９３およびＰＳＤの波形１９４を表示させる。個別指標１９２は、血圧の絶対値ＢＰ、日内の血圧変動性ΔＰＳＤおよび心機能の評価指標ＨＩを、すなわち心機能または循環器のリスクの大きさを、例えば３段階の色分け（例えば、赤、緑、青）で示される扇状のピクトグラムを含む。ＨＨＳ１９１は、個別指標１９２が示す血圧の絶対値ＢＰ、日内の血圧変動性ΔＰＳＤおよび心機能の評価指標ＨＩが分類された各段階の値のトータル値をスコアとして示す。ＰＳＤの波形１９４のうち斜線部分の波形は、圧受容器反射が作用する周波数帯域における波形を示している。
【０１４８】
医家は、図１９の画面から、直感的に被測定者（患者）の循環器または心機能（心血管）のリスクについての評価を把握することができる。
【０１４９】
図２０は、本発明の実施の形態に係る表示画面の他の例を模式的に示す図である。図２０の画面は血圧波形の情報を表示する。図２０を参照して、装置３０のＣＰＵ５０２は、ディスプレイ５１１に、血圧計１が測定した脈波による時系列の血圧波形２０１と、血圧波形２０１を部分的に拡大した波形２０２を表示させる。血圧計１が脈波の測定に並行して脈波データを装置３０に送信する場合は、リアルタイムに血圧波形２０１（２０２）が表示される。
【０１５０】
なお、図１９または図２０の情報の表示先は、装置３０のディスプレイ５１１に限定されず、端末装置１０のディスプレイ１５８または血圧計１のディスプレイ５０であってもよい。また、図１９または図２０の情報を、複数枚の画面に分割して表示してもよい。
【０１５１】
このように、本実施の形態では、被測定者の血圧を非侵襲的に１２時間よりも短い所定長時間において連続測定することで、血圧絶対値ＢＰと血圧変動性ΔＰＳＤ（標準偏差ＳＤ）を同時に評価することができる。また、それら血圧変動性を評価する指標に心機能の評価の指標ＨＩを加えることで、医家は、診療の効率化とともにこれまでに見逃された心機能のリスクの抽出や心不全にまで進行した高リスク，高血圧の患者の最適な患者管理をすることが期待できる。
【０１５２】
＜Ｍ．プログラム＞
上述したフローチャートを用いて説明した処理は、プログラムとして提供することもできる。このようなプログラムは、血圧計１のメモリ５１、端末装置１０のメモリ１５４、装置３０のストレージ５０５等の記憶媒体に格納されており、ＣＰＵが記憶媒体から当該プログラムを読出し、プログラムの命令コードを実行することにより処理が実現される。
【０１５３】
血圧計１の場合、このプログラムは、通信Ｉ／Ｆ５９を介して、通信回線を経由し、外部の情報処理装置（端末装置１０など）からメモリ５１ダウンロードされることにより供給され得る。または、カードＲ／Ｗ１０６を介して記憶媒体１０７からメモリ５１にロードされる。
【０１５４】
端末装置１０の場合、このプログラムは、通信Ｉ／Ｆ１６６を介して、通信回線を経由し、外部の情報処理装置（装置３０など）からメモリ１５４にダウンロードされる。または、メモリＩ／Ｆ１６４を介して記憶媒体１６５からメモリ１５４にロードされる。
【０１５５】
装置３０の場合、このプログラムは、通信Ｉ／Ｆ５０６を介して、通信回線を経由し、外部の情報処理装置からストレージ５０５にダウンロードされる。または、ドライバ５０８を介して記憶媒体５０９からストレージ５０５にロードされる。
【０１５６】
＜Ｎ．変形例＞
（ｎ１．変形例１）
血圧変動性を評価する指標ΔＰＳＤを算出するために用いるデータとして、本実施の形態では、血圧の波形データを用いたが、血圧の波形に限定されない。つまり、ΔＰＳＤは相対値であるから、ＰＳＤの絶対値は血圧変動性の評価に影響しない。したがって、血圧波形に関連する値であれば、血圧変動性を評価する指標ΔＰＳＤの算出に用いることができる。例えば、血圧計１において血圧の波形に校正する前の波形のデータ、すなわち圧力センサ３１から出力されるトノメトリ圧の変化を示す波形データ、または容積脈波の検出波形のデータであってもよい。容積脈波検出器には、例えば光電式脈波センサが含まれる。
【０１５７】
（ｎ２．変形例２）
上記の実施の形態では、収集部６００、指標算出部６０１、血圧算出部６０４および評価部６０５の全てを装置３０に備えたが、これらの少なくとも１つは、血圧計１または端末装置１０に備えられても良い。例えば、血圧計１が、収集部６００、指標算出部６０１、血圧算出部６０４および評価部６０５を備えても良い。また、血圧計１は血圧算出部６０４を備えて、端末装置１０または装置３０が収集部６００、指標算出部６０１または評価部６０５を備えても良い。
【０１５８】
（ｎ３．変形例３）
上記の実施の形態における脈波の測定（検出）方法であるトノメトリ法は、圧力センサ３１を経皮的に橈骨動脈９１に押し当てるだけであり装着がシンプルである利点がある。なお、脈波の測定方法はトノメトリ法に限定されない。例えば、容積補償法は、容積の動的な制御が必要でありシステムは本実施の形態の血圧計１に比較すると大型化し易いが、脈波の測定方法として、容積補償法を用いてもよい。
【０１５９】
（ｎ４．変形例４）
上記の実施の形態では、「時系列の脈波データ」は、脈波（血圧波形）の時間的な変化を、連続的に（または一定間隔をおいて不連続に）測定して得られた脈波のデータの系列（一連の脈波データ）を示す。したがって、本実施の形態では、１拍毎または５ミリ秒毎等を検出の間隔として連続的に脈波を検出（測定）することで時系列の脈波データが取得される。なお、脈波を「連続的」に検出するとの概念は、上記で述べた１拍毎または５ミリ秒毎等の略等間隔で脈波を検出する態様に限定されない。例えば、ノイズが含まれる部分（期間）のカットしながら連続的となるように脈波を検出する態様、また、例えば信頼性がない脈波を除きながら連続的となるように脈波を検出する態様、例えば１拍飛ばし（または２拍以上を飛ばしながら）で連続的となるように脈波を検出する態様等も、脈波を「連続的」に検出するとの概念に含まれ得る。
【０１６０】
＜Ｏ．付記＞
上述したような本実施の形態は、以下のような技術思想を含む。
【０１６１】
［構成１］
被測定者の測定部位から脈波を検出するセンサ（３１）を有した脈波検出装置（１）と、
前記脈波検出装置とデータを遣り取りする装置（３０）とを備えるシステム（１００）であって、
前記装置は、
１２時間未満の予め定められた長さの時間において前記センサにより連続して検出されるデータを時系列に収集する収集手段（６００）と、
収集される前記時系列の脈波のデータから、予め定められた周波数帯域のパワースペクトル密度を算出する周波数解析手段（６０１）を含み、
前記装置は、
前記予め定められた周波数帯域における前記パワースペクトル密度を代表する値を、心機能の健全性を評価する指標として算出する、システム。
【０１６２】
[構成２]
前記予め定められた長さの時間は、５分以上６０分以下の長さの時間を含む、構成１に記載のシステム。
【０１６３】
[構成３]
前記周波数解析手段は、自己回帰移動平均モデルの演算に従い、収集された前記時系列の脈波のデータから、前記予め定められた周波数帯域のパワースペクトル密度を算出する、構成１または２に記載のシステム。
【０１６４】
[構成４]
前記周波数解析手段は、収集された前記時系列の脈波のデータから、フーリエ変換により、前記予め定められた周波数帯域のパワースペクトル密度を算出する、構成１または２に記載のシステム。
【０１６５】
[構成５]
前記予め定められた周波数帯域は、生体の圧受容器反射が血圧変動に作用する周波数帯域に相当する、構成１から４のいずれか１に記載のシステム。
【０１６６】
[構成６]
前記予め定められた周波数帯域は、０．０１Ｈｚ〜０．１Ｈｚの帯域を含む、構成１から５のいずれか１に記載のシステム。
【０１６７】
[構成７]
前記代表する値は、前記予め定められた周波数帯域における前記パワースペクトル密度の大きさを示す値を含む、構成１から６のいずれか１に記載のシステム。
【０１６８】
[構成８]
前記大きさを示す値は、前記予め定められた周波数帯域における前記パワースペクトル密度の積分値を含む、構成７に記載のシステム。
【０１６９】
[構成９]
前記装置は、さらに、評価部（６０５）を備え、
前記評価部は、
前記心機能の健全性を評価する指標が示す値を対応の閾値と比較し、比較の結果に基き、当該指標の値を前記心機能のリスクの段階のいずれかに分類する、構成１から８のいずれか１に記載のシステム。
【０１７０】
[構成１０]
前記装置は、さらに、
前記予め定められた周波数帯域におけるパワースペクトル密度の変位を、前記日内の血圧変動性を評価する指標として算出する、構成９に記載のシステム。
【０１７１】
[構成１１]
前記評価部は、さらに、
前記日内の血圧変動性を評価する指標が示す値を対応の閾値と比較し、比較の結果に基き、当該指標の値を循環器のリスクの段階のいずれかに分類する、構成１０に記載のシステム。
【０１７２】
[構成１２]
前記装置は、さらに、
収集される時系列の脈波のデータから血圧を算出する、構成９から１１のいずれか１に記載にシステム。
【０１７３】
[構成１３]
前記評価部は、さらに、
前記血圧を対応の閾値と比較し、比較の結果に基き、当該血圧を循環器のリスクの段階のいずれかに分類する、構成１２に記載のシステム。
【０１７４】
[構成１４]
前記装置は、さらに、
前記リスクの段階を示す情報を、表示デバイスに表示させる、構成９から１３のいずれか１に記載のシステム。
【０１７５】
[構成１５]
前記装置は、さらに、
前記リスクの段階に関連付けられた治療方針を、表示デバイスに表示させる、構成９から１４のいずれか１に記載にシステム。
【０１７６】
[構成１６]
１２時間未満の予め定められた長さの時間において測定部位からセンサ（３１）により連続して検出される脈波のデータを時系列に収集する手段（６００）と、
収集される時系列の脈波のデータから、パワースペクトル密度を算出する周波数解析手段（６０１）と、を備え、
予め定められた周波数帯域における前記パワースペクトル密度を代表する値を、心機能の健全性を評価する指標として算出する、装置。
【０１７７】
[構成１７]
プログラムであって、前記プログラムはコンピュータに、
１２時間未満の予め定められた長さの時間において測定部位からセンサ（３１）により連続して検出される脈波のデータを時系列に収集するステップと、
収集される時系列の脈波のデータから、パワースペクトル密度を算出するステップと、
予め定められた周波数帯域における前記パワースペクトル密度を代表する値を、心機能の健全性を評価する指標として算出するステップを実行させる、プログラム。
【０１７８】
[構成１８]
被測定者の測定部位から脈波を検出するセンサ（３１）を有した脈波検出装置（１）と、
前記脈波検出装置とデータを遣り取りする装置（３０）とを備えるシステム（１００）であって、
前記装置は、
１２時間未満の予め定められた長さの時間において前記センサにより連続して検出されるデータを時系列に収集する収集手段（６００）と、
収集される前記時系列の脈波のデータから、予め定められた周波数帯域の脈波のデータを抽出する手段を含み、
前記装置は、
抽出された前記脈波のデータが示す血圧のばらつきの大きさを、心機能の健全性を評価する指標として算出する、システム（１００）。
【０１７９】
[構成１９]
１２時間未満の予め定められた長さの時間において測定部位からセンサにより連続して検出される脈波のデータを時系列に収集する手段（６００）と、
収集される前記時系列の脈波のデータから、予め定められた周波数帯域の脈波のデータを抽出する手段を含み、
抽出された前記脈波のデータが示す血圧のばらつきの大きさを、心機能の健全性を評価する指標として算出する、装置（３０）。
【０１８０】
[構成２０]
プログラムであって、前記プログラムはコンピュータに、
１２時間未満の予め定められた長さの時間において測定部位からセンサにより連続して検出される脈波のデータを時系列に収集するステップと、
収集される前記時系列の脈波のデータから、予め定められた周波数帯域の脈波のデータを抽出するステップと、
抽出された前記脈波のデータが示す血圧のばらつきの大きさを、心機能の健全性を評価する指標として算出するステップと、を実行させる、プログラム。
【０１８１】
今回開示された実施の形態はすべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記した説明ではなく、特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
【符号の説明】
【０１８２】
１  血圧計、１０，１０Ａ，１０Ｂ  端末装置、２０  ベルト、２１  カフ、３０  装置、３１  圧力センサ、３２  ポンプ、３３  弁、４１，４３  ネットワーク、５０，１５８，５１１  ディスプレイ、７０  本体、９０  左手首、９１  橈骨動脈、１００  システム、１０７，１６５，５０９  記憶媒体、１５３  タイマ、１９２  個別指標、１９３，２０１  血圧波形、１９４，２０２  波形、３００  テーブル、３１０  発振回路、３２０  ポンプ駆動回路、５０１  コントローラ、６００  収集部、６０１  指標算出部、６０３  周波数解析部、６０４  血圧算出部、６０５  評価部、ＢＰ  血圧絶対値、Ｆ１  自己回帰移動平均モデル、ＨＩ  指標、Ｒ  相関係数、ＳＤ  標準偏差、ＰＳＤ  パワースペクトル密度、ΔＰＳＤ  ＰＳＤの変位。
【図１】
2019213820.tif 000003
【図２】
2019213820.tif 000004
【図３】
2019213820.tif 000005
【図４】
2019213820.tif 000006
【図５】
2019213820.tif 000007
【図６】
2019213820.tif 000008
【図７】
2019213820.tif 000009
【図８】
2019213820.tif 000010
【図９】
2019213820.tif 000011
【図１０】
2019213820.tif 000012
【図１１】
2019213820.tif 000013
【図１２】
2019213820.tif 000014
【図１３】
2019213820.tif 000015
【図１４】
2019213820.tif 000016
【図１５】
2019213820.tif 000017
【図１６】
2019213820.tif 000018
【図１７】
2019213820.tif 000019
【図１８】
2019213820.tif 000020
【図１９】
2019213820.tif 000021
【図２０】
2019213820.tif 000022
【図２１】
2019213820.tif 000023
【図２２】
2019213820.tif 000024
【図２３】
2019213820.tif 000025
【図２４】
2019213820.tif 000026
【図２５】
2019213820.tif 000027