JP2020049134A - 血圧推定装置および血圧推定プログラム - Google Patents

Abstract

【課題】特別な装置や煩雑な校正手順を必要とせず、循環動態（血圧、心拍出量、及び全末梢血管抵抗の関係性）の理論的な背景を基盤に、合理的な精度で血圧を推定することが可能な血圧推定装置および血圧推定プログラムの提供。【解決手段】光電容積脈波測定により心拍数（ＨＲ；heart rate）および修正規準化脈波容積（ｍＮＰＶ；modified normalized pulse volume）を得る測定手段２と、測定手段により得られた心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）から血圧（ＢＰ；blood pressure）をＢＰ＝ｅｘｐ［ａ（ｌｎ ＨＲ）＋ｂ（ｌｎ ｍＮＰＶ）＋ｃ］（但し、ａ，ｂ，ｃは補正係数）により算出する演算手段３とを含む。【選択図】図１

Description

本発明は、医療やヘルスケア等の分野において血圧を推定するための血圧推定装置および血圧推定プログラムに関する。

超高齢社会の進行に伴い高血圧患者が増加し続けており、予防および治療のための血圧管理が重要となっている。しかしながら、現状では一般的に血圧値の参考とされるのは外来血圧または家庭血圧であり、加えて年に１回ないしは数回の健康診断で測定するのみに留まっている人の割合が大きい。これでは短期的な血圧変動等の重要な指標を見逃してしまうことになる。

従来、最も利用されている血圧測定方法は、上腕または手首にカフを巻いて測定するカフ式血圧計を使用する方法（カフ振動法）である。しかしながら、この血圧測定方法では、カフ適用部分を圧迫するため、被験者に不快感を与えてしまうという問題がある。また、カフ付属の専用の血圧測定装置（カフ式血圧計）が必要になるため、血圧を日常的に測定するには向いていないと言える。そこで、カフを使用せずに血圧を測定するカフレス血圧測定方法として、例えば特許文献１〜３に記載のものが開発されている。

特許文献１に記載のカフレス血圧測定方法は、心電図および光電容積脈波（ＰＰＧ；Photo-plethysmogram）を測定するものであり、被験者に装着された生体信号検出センサによって得られる生体信号（心電図と脈波）から脈波伝搬時間（ＰＴＴ；pulse transit time）を算出し、血管に関する特性を表す２つの定数パラメータ（Ｂ１，Ｂ２）を用いて収縮期血圧を算出する。
ここで、定数パラメータ（Ｂ１，Ｂ２）は、あらかじめ被験者の血圧を変動させながら取得した参照データセットを用いて決定しておく（キャリブレーション）必要がある。

特許文献２に記載のカフレス血圧測定方法は、光電容積脈波（ＰＰＧ）の周波数特性を使用するものであり、生体の動脈における脈波から得られる心拍数情報から、低周波信号成分ＬＦと高周波信号成分ＨＦとの比を算出し、交感神経の活動状態および副交感神経の活動状態を判定し、血圧推定に利用する。

特許文献３に記載のカフレス血圧測定方法は、光電容積脈波（ＰＰＧ）のみを測定するものであり、脈波を測定した後、各波高値、波高値の比、脈波立ち上がり点間の時間、各波相互の時間間隔、脈波の積分値、脈拍数などの特徴量を計算し、同様に速度脈波や加速度脈波についても特徴量を計算する。このカフレス血圧測定方法は、これらを用いて血圧の基準値との関係を現場で徐々に学習させ、徐々に精度を向上させる手法である。

特許第５９８４０８８号公報 特許第３２１３２７８号公報 特開平１０−２９５６５７号公報

上記特許文献１，２に記載のカフレス血圧測定方法では、合理的な精度で血圧を推定することを可能にしているが、心電図と光電容積脈波を測定する特別な装置が必要になってしまうことに加え、前述のカフ式血圧計を用いた煩雑な校正手順を必要とするため、実用的とは言い難い。

また、上記特許文献３に記載のカフレス血圧測定方法では、血圧推定のプロセスにおいて根拠とする生理学的背景がはっきりしておらず、現時点では使用の限界や適合性を見積もることが難しいという問題がある。

そこで、本発明においては、特別な装置や煩雑な校正手順を必要とせず、循環動態（血圧、心拍出量、及び全末梢血管抵抗の関係性）の理論的な背景を基盤に、合理的な精度で血圧を推定することが可能な血圧推定装置および血圧推定プログラムを提供することを目的とする。

本発明の血圧推定装置は、光電容積脈波測定により心拍数（ＨＲ；heart rate）および修正規準化脈波容積（ｍＮＰＶ；modified normalized pulse volume）を得る測定手段と、測定手段により得られた心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）から血圧（ＢＰ；blood pressure）を
ＢＰ＝ｅｘｐ［ａ（ｌｎ ＨＲ）＋ｂ（ｌｎ ｍＮＰＶ）＋ｃ］
（但し、ａ，ｂ，ｃは補正係数）
により算出する演算手段とを含むものである。

循環生理学の知見（循環器系をおおざっぱに線形系と考える）によれば、血行力学は概してオームの法則を満たしており、オームの法則「電圧（Ｖ）＝電流（Ｉ）×抵抗（Ｒ）」になぞらえれば、心臓血管系において血圧（ＢＰ）は、心拍出量（ＣＯ；cardiac output）と全末梢血管抵抗（ＴＰＲ；total peripheral resistance）の積により次のように求められる。
ＢＰ＝ＣＯ×ＴＰＲ
この方程式の両辺に対数変換（ｌｎ）を施すと次のようになる。
ｌｎ ＢＰ＝ｌｎ ＣＯ＋ｌｎ ＴＰＲ

ここで、心拍出量（ＣＯ）と心拍数（ＨＲ）は、両方ともβ−アドレナリン作動性交感神経活動の影響を受ける点で、心拍数（ＨＲ）に関連する。同様に、全末梢血管抵抗（ＴＰＲ）と修正規準化脈波容積（ｍＮＰＶ）は両方ともα−アドレナリン作動性の交感神経活動の影響を受けるという点で、指の血管の緊張度（収縮度）を反映する測定値である。したがって、血圧（ＢＰ）は単純な線形多項式によって推定することができる。

すなわち、
ｌｎ ＢＰ＝ａ×ｌｎ ＨＲ＋ｂ×ｌｎ ｍＮＰＶ＋ｃ
（但し、ａ，ｂ，ｃは補正係数）
これを変形して、
ＢＰ＝ｅｘｐ［ａ（ｌｎ ＨＲ）＋ｂ（ｌｎ ｍＮＰＶ）＋ｃ］
ここで、心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）の両方は、光電容積脈波を使用して測定できる生理学的変数である。

すなわち、血圧推定に必要なものは、光電容積脈波測定から得られる心拍数（ＨＲ）と修正規準化脈波容積（ｍＮＰＶ）だけであるため、光電容積脈波による心拍数（ＨＲ）と修正規準化脈波容積（ｍＮＰＶ）を測定できさえすれば、手軽に血圧を推定することができ、本発明の血圧推定装置によれば、光電容積脈波測定により得られた心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）から、前述の循環動態の理論的な背景（ＢＰ＝ＣＯ×ＴＰＲ）を基盤に合理的な精度で血圧（ＢＰ）を推定することができる。

ここで、血圧（ＢＰ）として、平均動脈圧（ＭＡＰ；mean arterial pressure）を算出する場合、カフを用いた上腕測定による平均動脈圧（ＭＡＰ）を独立変数とし、心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）を従属変数として、
ＭＡＰ＝ｅｘｐ［ａ_MAP（ｌｎ ＨＲ）＋ｂ_MAP（ｌｎ ｍＮＰＶ）＋ｃ_MAP］
を重回帰分析して得られた補正係数ａ_MAP，ｂ_MAP，ｃ_MAPを用いることにより、合理的な精度で平均動脈圧（ＭＡＰ）を推定することができる。

さらに、収縮期血圧（ＳＢＰ；systolic ＢＰ）は心拍数（ＨＲ）および心拍出量（ＣＯ）と密接に関連し、拡張期血圧（ＤＢＰ；diastolic ＢＰ）は全末梢血管抵抗（ＴＰＲ）および修正規準化脈波容積（ｍＮＰＶ）と密接に関連していることが知られており、収縮期血圧（ＳＢＰ）は心拍数（ＨＲ）のより高い寄与により、また、拡張期血圧（ＤＢＰ）は修正規準化脈波容積（ｍＮＰＶ）のより高い寄与により推定することができる。すなわち、平均動脈圧（ＭＡＰ）、収縮期血圧（ＳＢＰ）および拡張期血圧（ＤＢＰ）のそれぞれに関して、補正係数ａ，ｂ，ｃを個別に用意することで、上記の方程式１つで３つの血圧値を推定することが可能となる。

また、本発明の血圧推定プログラムは、光電容積脈波測定により心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）を得る測定手段と、測定手段により得られた心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）から血圧（ＢＰ）を、
ＢＰ＝ｅｘｐ［ａ（ｌｎ ＨＲ）＋ｂ（ｌｎ ｍＮＰＶ）＋ｃ］
（但し、ａ，ｂ，ｃは補正係数）
により算出する演算手段としてコンピュータを機能させるためのものである。このプログラムを実行したコンピュータによれば、上記本発明の血圧推定装置と同様の作用、効果を奏することができる。

本発明によれば、ユーザはカフによる圧迫という負荷から解放され、かつ特別な装置も煩雑な校正手順も必要とすることなく、スマートフォン（小型携帯端末機器）をはじめとする日常的な機器を利用することで、容易に血圧を推定することが可能となる。

本発明の実施の形態における血圧推定装置の機能ブロック図である。 図１の血圧推定装置のハードウェア構成を示すブロック図である。 図２の血圧推定装置の使用時の様子を示す図であって、（Ａ）は血圧推定装置の外観を示す図、（Ｂ）は血圧測定時の持ち方を示す図である。 本実施形態における血圧推定装置による血圧推定のフローチャートである。 実験手順を示す説明図である。 実験結果を示す図である。

図１は本発明の実施の形態における血圧推定装置の機能ブロック図、図２は図１の血圧推定装置のハードウェア構成を示すブロック図、図３は図２の血圧推定装置の使用時の様子を示す図であって、（Ａ）は血圧推定装置の外観を示す図、（Ｂ）は血圧測定時の持ち方を示す図である。

図１に示すように、本発明の実施の形態における血圧推定装置１は、光電容積脈波測定により心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）を得る測定手段２と、測定手段２により得られた心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）から血圧（ＢＰ）を算出する演算手段３と、演算手段３により算出した血圧（ＢＰ）を出力する出力手段４とを有する。

血圧推定装置１は、図２に示すように、光源１０、駆動回路１１、光検出器１２、フィルタ１３、Ａ／Ｄ変換器１４、ＣＰＵ１５、ＲＡＭ１６、ＲＯＭ１７、通信手段１８、表示器１９や電源２０等を備える。血圧推定装置１は、例えば図３に示すスマートフォン２１により構成される。スマートフォン２１は、光源１０としてＬＥＤ（フラッシュ用白色ＬＥＤ）２２を備えている。また、光検出器１２としてカメラ（ＣＭＯＳカメラ）２３を備えている。

スマートフォン２１はコンピュータの一種であり、ＲＯＭ１７に記憶されたＯＳをＣＰＵ１５によって実行する。また、スマートフォン２１は、ＯＳ上で上記測定手段２、演算手段３および出力手段４として機能させるための血圧推定プログラムを実行することにより、本実施形態における血圧推定装置１として機能する。血圧推定プログラムはＲＡＭ１６に記憶され、ＣＰＵ１５によって実行される。

血圧測定時には、スマートフォン２１を片手で持ち、１本の指でＬＥＤ２２とカメラ２３の両方を覆うようにする。測定手段２は、駆動回路１１を制御することでＬＥＤ２２（光源１０）から光電容積脈波を取得するための光を発し、カメラ２３（光検出器１２）を制御して反射光画像を撮影し、フィルタ１３およびＡ／Ｄ変換器１４を介して入力する。スマートフォン２１は、１本の指でＬＥＤ２２とカメラ２３の両方を覆って生体によって多重散乱された反射光を撮影するため、ＬＥＤ２２とカメラ２３の位置は比較的近くにあることが必要である。

測定手段２は、この入力した反射光画像に基づいて光電容積脈波を取得し、心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）を得る。この光電容積脈波測定により心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）を得る方法については公知（K. Matsumura and T. Yamakoshi, iPhysioMeter: A new approach for measuring heart rate and normalized pulse volume using only a smartphone. Behavior Research Methods, 45(4), 1272-1278, 2013. doi:10.3758/s13428-012-0312-z）であり、その詳細な説明は省略する。

演算手段３は、測定手段２により得られた心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）から血圧（ＢＰ）を
ＢＰ＝ｅｘｐ［ａ（ｌｎ ＨＲ）＋ｂ（ｌｎ ｍＮＰＶ）＋ｃ］
により算出する。

ここで、ａ，ｂ，ｃは補正係数であり、予めカフを用いた上腕測定による血圧（ＢＰ）を独立変数とし、心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）を従属変数として、
ＢＰ＝ｅｘｐ［ａ（ｌｎ ＨＲ）＋ｂ（ｌｎ ｍＮＰＶ）＋ｃ］
を重回帰分析して得られたものを用いる。

例えば、血圧（ＢＰ）として平均動脈圧（ＭＡＰ）を算出する場合、カフを用いた上腕測定による平均動脈圧（ＭＡＰ）を独立変数とし、心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）を従属変数として、
ＭＡＰ＝ｅｘｐ［ａ_MAP（ｌｎ ＨＲ）＋ｂ_MAP（ｌｎ ｍＮＰＶ）＋ｃ_MAP］
を重回帰分析して得られた補正係数ａ_MAP，ｂ_MAP，ｃ_MAPを用いる。

また、血圧（ＢＰ）として収縮期血圧（ＳＢＰ）を算出する場合には、カフを用いた上腕測定による収縮期血圧（ＳＢＰ）を独立変数とし、心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）を従属変数として、
ＳＢＰ＝ｅｘｐ［ａ_SBP（ｌｎ ＨＲ）＋ｂ_SBP（ｌｎ ｍＮＰＶ）＋ｃ_SBP］
を重回帰分析して得られた補正係数ａ_SBP，ｂ_SBP，ｃ_SBPを用いる。

さらに、血圧（ＢＰ）として拡張期血圧（ＤＢＰ）を算出する場合には、カフを用いた上腕測定による拡張期血圧（ＤＢＰ）を独立変数とし、心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）を従属変数として、
ＤＢＰ＝ｅｘｐ［ａ_DBP（ｌｎ ＨＲ）＋ｂ_DBP（ｌｎ ｍＮＰＶ）＋ｃ_DBP］
を重回帰分析して得られた補正係数ａ_DBP，ｂ_DBP，ｃ_DBPを用いる。

出力手段４は、この演算手段３により算出した平均動脈圧（ＭＡＰ）、収縮期血圧（ＳＢＰ）および拡張期血圧（ＤＢＰ）の各血圧（ＢＰ）を、通信手段１８や表示器１９により出力する。

図４は本実施形態における血圧推定装置１による血圧推定のフローチャートである。
図４に示すように、まず測定手段２の初期化を行い（Ｓ１０１）、測定手段２により光電容積脈波測定を行う（Ｓ１０２）。撮影した画像の明るさ判定を行い（Ｓ１０３）、明るさが規定範囲を超えているかどうか判定し（Ｓ１０４）、超えている場合にはカメラ２３（光検出器１２）のゲインを変更する（Ｓ１０５）。

そして、測定手段２により心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）を算出し（Ｓ１０６）、演算手段３により平均動脈圧（ＭＡＰ）、収縮期血圧（ＳＢＰ）および拡張期血圧（ＤＢＰ）の各血圧（ＢＰ）を推定する（Ｓ１０７）。演算手段３の結果は表示器１９に表示したり、通信手段１８により外部へ出力したりする（Ｓ１０８）。

この血圧推定装置１では、ユーザはカフによる圧迫という負荷から解放され、かつ特別な装置も煩雑な校正手順も必要とすることなく、スマートフォン２１をはじめとする日常的な機器を利用することで、容易に血圧を推定することが可能である。

なお、スマートフォン２１以外に、例えば、光電容積脈波測定を行うことが可能な以下の機器を利用することが可能である。
・補聴器、イヤホン（耳周辺に固定することができ、比較的安定した光電容積脈波測定が可能）
・ポインティングデバイスとしてのマウス（指が当たる部分にセンサを実装することで、比較的良好な光電容積脈波測定が可能）
・指輪
・スマートウォッチ（常時身に着けるものであり、無意識的な光電容積脈波測定が可能）

健常成人１３名（女性６名、男性７名、平均年齢２０．６歳）により実験を行った。実験に用いた計測装置は、上記血圧推定装置１として動作するｉＰｈｏｎｅ ６ｓ（アップル社製）と、従来のカフ式血圧計としての上腕式自動血圧計（ＤＳ−Ｓ１０（日本精密測定器株式会社製））である。計測部位は、ｉＰｈｏｎｅ ６ｓが左手第２指、カフ式血圧計が右上腕部である。

図５は実験手順を示している。図５に示すように、計測装置の取付け完了後、同時に計測を行い、その後、３分間の安静状態と、３分間の課題実施を行い、カフ式血圧計により９０秒ごとに計４回血圧を測定した。なお、血圧推定装置１（ｉＰｈｏｎｅ ６ｓ）では連続計測を行った。課題実施は、血圧を変動させるために実施するものであり、５０００から１３を引いていく暗算課題（能動的ストレッサー課題）を行った。

この実験によりカフ式血圧計から得られた平均動脈圧（ＭＡＰ）、収縮期血圧（ＳＢＰ）および拡張期血圧（ＤＢＰ）をそれぞれ独立変数とし、血圧推定装置１（ｉＰｈｏｎｅ ６ｓ）により得られた心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）を従属変数として前述の重回帰分析を行い、平均動脈圧（ＭＡＰ）、収縮期血圧（ＳＢＰ）および拡張期血圧（ＤＢＰ）の補正係数ａ_MAP，ａ_SBP，ａ_DBP，ｂ_MAP，b_SBP，b_DBP，c_MAP，c_SBP，c_DBPをそれぞれ得た。

図６（Ａ）は平均動脈圧（ＭＡＰ）、収縮期血圧（ＳＢＰ）および拡張期血圧（ＤＢＰ）の血圧推定装置１（ｉＰｈｏｎｅ ６ｓ）だけを用いて推定した値と従来のカフ式血圧計で測定した値の散布図（Ｎ＝４９）を示している。上列の実線は回帰直線を表しており、その式をｒ値（Pearsonの相関係数）とともに各散布図に示している。下列は対応するBland-Altman plots（誤差分析）であり、実線および破線はそれぞれＭ（平均）およびＭ＋ＳＤ（標準偏差）の範囲を示している。平均＝（推定血圧＋上腕血圧）／２、差＝上腕血圧−推定血圧である。

図６（Ｂ）は平均動脈圧（ＭＡＰ）、収縮期血圧（ＳＢＰ）および拡張期血圧（ＤＢＰ）を近赤外光（８１０ ｎｍ）の従来式光電容積脈波（ＰＰＧ）を用いて推定した値と従来のカフ式血圧計で測定した値の散布図（Ｎ＝５１）を示している。上列と下列の精度評価図は上と同じである。

図６（Ａ）から分かるように、血圧推定装置１では、構成を伴い，かつ理論的根拠の乏しい図６（Ｂ）のＰＷＶベースの手法｛｜ｒ｜s＞＞０．７４（X. Ding et al.: Pulse transit time based continuous cuffless blood pressure estimation: A new extension and a comprehensive evaluation. Scientific Reports, 7, 11554, 2017. doi:10.1038/s41598-017-11507-3），｜ｒ｜ｓ＜０．６７（R. A. Payne, et al.: Pulse transit time measured from the ECG: an unreliable marker of beat-to-beat blood pressure. Journal of Applied Physiology, 100, 136-141, 2006. doi:10.1152/japplphysiol.00657.2005）｝と比較しても、遜色ない結果を得ることができる。

本発明の血圧推定装置および血圧推定プログラムは、医療やヘルスケア等の分野において血圧を推定するための装置およびプログラムとして有用である。

１ 血圧推定装置
２ 測定手段
３ 演算手段
４ 出力手段
１０ 光源
１１ 駆動回路
１２ 光検出器
１３ フィルタ
１４ Ａ／Ｄ変換器
１５ ＣＰＵ
１６ ＲＡＭ
１７ ＲＯＭ
１８ 通信手段
１９ 表示器
２０ 電源
２１ スマートフォン
２２ ＬＥＤ
２３ カメラ

Claims (5)

1. 光電容積脈波測定により心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）を得る測定手段と、
   前記測定手段により得られた心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）から血圧（ＢＰ）を
   ＢＰ＝ｅｘｐ［ａ（ｌｎ ＨＲ）＋ｂ（ｌｎ ｍＮＰＶ）＋ｃ］
   （但し、ａ，ｂ，ｃは補正係数）
   により算出する演算手段と
   を含む血圧推定装置。
2. 前記演算手段は、カフを用いた上腕測定による平均動脈圧（ＭＡＰ）を独立変数とし、前記心拍数（ＨＲ）および前記修正規準化脈波容積（ｍＮＰＶ）を従属変数として、
   ＭＡＰ＝ｅｘｐ［ａ_MAP（ｌｎ ＨＲ）＋ｂ_MAP（ｌｎ ｍＮＰＶ）＋ｃ_MAP］
   を重回帰分析して得られた補正係数ａ_MAP，ｂ_MAP，ｃ_MAPを用いることにより、前記血圧（ＢＰ）としての平均動脈圧（ＭＡＰ）を算出するものである
   請求項１記載の血圧推定装置。
3. 前記演算手段は、カフを用いた上腕測定による収縮期血圧（ＳＢＰ）を独立変数とし、前記心拍数（ＨＲ）および前記修正規準化脈波容積（ｍＮＰＶ）を従属変数として、
   ＳＢＰ＝ｅｘｐ［ａ_SBP（ｌｎ ＨＲ）＋ｂ_SBP（ｌｎ ｍＮＰＶ）＋ｃ_SBP］
   を重回帰分析して得られた補正係数ａ_SBP，ｂ_SBP，ｃ_SBPを用いることにより、前記血圧（ＢＰ）としての収縮期血圧（ＳＢＰ）を算出するものである
   請求項１または２に記載の血圧推定装置。
4. 前記演算手段は、カフを用いた上腕測定による拡張期血圧（ＤＢＰ）を独立変数とし、前記心拍数（ＨＲ）および前記修正規準化脈波容積（ｍＮＰＶ）を従属変数として、
   ＤＢＰ＝ｅｘｐ［ａ_DBP（ｌｎ ＨＲ）＋ｂ_DBP（ｌｎ ｍＮＰＶ）＋ｃ_DBP］
   を重回帰分析して得られた補正係数ａ_DBP，ｂ_DBP，ｃ_DBPを用いることにより、前記血圧（ＢＰ）としての拡張期血圧（ＤＢＰ）を算出するものである
   請求項１から３のいずれか１項に記載の血圧推定装置。
5. 光電容積脈波測定により心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）を得る測定手段と、
   前記測定手段により得られた心拍数（ＨＲ）および修正規準化脈波容積（ｍＮＰＶ）から血圧（ＢＰ）を、
   ＢＰ＝ｅｘｐ［ａ（ｌｎ ＨＲ）＋ｂ（ｌｎ ｍＮＰＶ）＋ｃ］
   （但し、ａ，ｂ，ｃは補正係数）
   により算出する演算手段と
   してコンピュータを機能させるための血圧推定プログラム。