(57)【要約】      （修正有）【課題】心拍変動解析要件を満たさない周波数において、信頼度の高い自律神経活性度を算出する装置を提供する。【解決手段】予め定められた周波数を用いて生体信号を取得する生体信号取得部１０１と、生体信号取得部１０１により取得した生体信号から低周波成分と高周波成分を算出する周波数成分算出部１０２と、周波数成分算出部１０２により算出された低周波成分及び高周波成分の誤差を、予め定められた周波数を用いて補正する誤差補正部１０３と、誤差補正部１０３により誤差が補正された低周波成分及び高周波成分に基づいて自律神経活性度を算出する自律神経活性度算出部１０４と、を備える。【選択図】図２
2016043191.tif 000002
【特許請求の範囲】
【請求項１】
予め定められた周波数を用いて生体信号を取得する生体信号取得手段と、
前記生体信号取得手段により取得した生体信号から低周波成分と高周波成分を算出する周波数成分算出手段と、
前記周波数成分算出手段により算出された前記低周波成分及び前記高周波成分の誤差を、前記予め定められた周波数を用いて補正する誤差補正手段と、
前記誤差補正手段により誤差が補正された前記低周波成分及び前記高周波成分に基づいて自律神経活性度を算出する自律神経活性度算出手段と、
を備えることを特徴とする生体信号解析装置。
【請求項２】
前記誤差補正手段は、前記周波数成分算出手段により算出された前記低周波成分から前記低周波成分の誤差成分を差し引くことにより前記低周波成分の誤差を補正することを特徴とする請求項１記載の生体信号解析装置。
【請求項３】
前記誤差補正手段は、前記周波数成分算出手段により算出された前記高周波成分から前記高周波成分の誤差成分を差し引くことにより前記高周波成分の誤差を補正することを特徴とする請求項１又は２記載の生体信号解析装置。
【請求項４】
前記誤差補正手段は、前記予め定められた周波数をｆｓと、前記低周波成分の誤差成分をＬＦｄｄとするとき、ＬＦｄｄ＝３４７０／ｆｓ2［ｍｓｅｃ2］に基づいて前記低周波成分の誤差成分を算出することを特徴とする請求項１から３の何れか１項に記載の生体信号解析装置。
【請求項５】
前記誤差補正手段は、前記予め定められた周波数をｆｓと、前記高周波成分の誤差成分をＨＦｄｄとするとき、ＨＦｄｄ＝４７５３５／ｆｓ2［ｍｓｅｃ2］に基づいて前記高周波成分の誤差成分を算出することを特徴とする請求項１から４の何れか１項に記載の生体信号解析装置。
【請求項６】
予め定められた周波数を用いて生体信号を取得する生体信号取得装置と情報処理装置とがネットワークを介して接続された生体信号解析システムであって、
前記情報処理装置は、
前記生体信号取得装置により取得された生体信号を受け付ける受付手段と、
前記受付手段により受け付けられた生体信号から低周波成分と高周波成分を算出する周波数成分算出手段と、
前記周波数成分算出手段により算出された前記低周波成分及び前記高周波成分の誤差を、前記予め定められた周波数を用いて補正する誤差補正手段と、
前記誤差補正手段により誤差が補正された前記低周波成分及び前記高周波成分に基づいて自律神経活性度を算出する自律神経活性度算出手段と、
を備えることを特徴とする生体信号解析システム。
【請求項７】
予め定められた周波数を用いて生体信号を取得するステップと、
取得した前記生体信号から低周波成分と高周波成分を算出するステップと、
算出された前記低周波成分及び前記高周波成分の誤差を、前記予め定められた周波数を用いて補正するステップと、
誤差が補正された前記低周波成分及び前記高周波成分に基づいて自律神経活性度を算出するステップと、
を備えることを特徴とする生体信号解析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、生体信号解析装置、生体信号解析システム、及び、生体信号解析方法に関する。
【背景技術】
【０００２】
従来から、心拍のピーク間隔より生成した心拍変動波形の解析を行うことにより、心拍変動性指標（以下、「ＨＲＶ」という。）が得られることが既に知られている。このＨＲＶは、交感・副交感神経系の活動水準を反映する指標として知られている。
【０００３】
そして、ＨＲＶの周波数成分中における低周波（以下、「ＬＦ」という。）成分、高周波（以下、「ＨＦ」という。）成分のピークパワー値等を検討することで精神的作業負荷、心的作業負荷、覚醒度、緊張度の定量化が行える可能性が示唆されている。ＨＲＶは、医学的に検討価値の高い指標として知られている。なお、一般的に、ＬＦ成分は、０．０４〜０．１５サイクル／ビートの周波数成分をいい、ＨＦ成分は、０．１５〜０．４０サイクル／ビートの周波数成分をいう。
【０００４】
即ち、ＬＦ成分とは交換神経系を反映する血圧変動性の成分が低周波であることから重要となっている成分であり、そのパワーは精神的緊張の増大、起立性の刺激、例えば姿勢の変化などにより増大するものとして認識されている。また、ＨＦ成分とは呼吸変動性の成分が高周波であることからこのように称されるものであり、安静状態や睡眠中に高い値を示し緊張度の増大により消失傾向に向かうことが知られているものである。
【０００５】
ＨＲＶとしては、一般的に心電図から取得した心電波形を用いるが、指や耳たぶの容積脈波を用いることもできる。容積脈波は例えば皮膚血管の膨張・収縮を皮膚表面から光学的／電気的／機械的に捕らえたものである。特許文献１には、容積脈波の２次微分波すなわち加速度脈波のピーク間隔が心拍のピーク間隔と高い相関があることが記載されている。
【０００６】
容積脈波を光学的に捕らえる手法の一つとして、例えばＲＧＢカメラ等のイメージセンサにより容積脈波を取得するフォトプレチスモグラフィが知られている。また、測定対象である人の指尖や顔を定点撮影することで、血管のヘモグロビン吸収ピークがＧ信号帯域[５００〜６００ｎｍ]にあることから、Ｇ信号帯域の時系列変化が一定の周期、すなわち容積脈波として現れることが既に知られている。
【０００７】
また、非特許文献１には、スマートフォンのカメラで取得した指尖容積脈波から心拍変動を解析することが記載されており、これは心電計以外のデバイスで自律神経活性度を測定できることを示している。
【発明の概要】
【発明が解決しようとする課題】
【０００８】
ここで、非特許文献２にはＨＲＶガイドラインが記載されており、その内の心拍変動解析要件に関し、生体信号のＡ／Ｄ変換のサンプリング周波数として、２５０Ｈｚ以上を推奨することが記載されている。
【０００９】
他方、特許文献２には、低サンプリング周波数で心拍変動を推定する目的で、スペクトル解析とバンドパスフィルタを利用した心拍変動推定方法が開示されている。
【００１０】
しかし、非特許文献１におけるスマートフォンのカメラのサンプリング周波数は２４．９９９Ｈｚと記載されており、例えば非特許文献１に記載された心拍変動解析要件である２５０Ｈｚ以上を満たしていない。心拍変動解析要件を満たさない周波数で心拍変動解析を実施する場合、自律神経活性度を示すＬＦ／ＨＦにおいて測定誤差が発生する。しかし、どの程度の測定誤差が発生するのか分からないため、測定結果の信頼度に影響するという問題がある。
【００１１】
本発明は、このような実情に鑑みてなされたものであって、心拍変動解析要件を満たさない周波数において、信頼度の高い自律神経活性度を算出することを目的とする。
【課題を解決するための手段】
【００１２】
上記の課題を解決するため、本発明の生体信号解析装置は、予め定められた周波数を用いて生体信号を取得する生体信号取得手段と、生体信号取得手段により取得した生体信号から低周波成分と高周波成分を算出する周波数成分算出手段と、周波数成分算出手段により算出された低周波成分及び高周波成分の誤差を、予め定められた周波数を用いて補正する誤差補正手段と、誤差補正手段により誤差が補正された低周波成分及び高周波成分に基づいて自律神経活性度を算出する自律神経活性度算出手段と、を備えることを特徴とする。
【発明の効果】
【００１３】
本発明によれば、心拍変動解析要件を満たさない周波数において、信頼度の高い自律神経活性度を算出することが可能となる。
【図面の簡単な説明】
【００１４】
【図１】本発明の実施形態における生体信号解析装置のハードウェア構成図である。
【図２】本発明の実施形態における生体信号解析装置に係る機能ブロック図である。
【図３】本発明の実施形態における脈拍変動の算出について説明する模式図である。
【図４】本発明の実施形態における周波数解析の一例について説明する模式図である。
【図５】本発明の実施形態における誤差成分の算出について説明する模式図である。
【図６】本発明の実施形態における生体信号解析処理手順について説明するフローチャートである。
【図７】図６におけるＬＦ及びＨＦの算出手順の詳細について説明するフローチャートである。
【図８】本発明の実施形態を適用した一例について説明する模式図である。
【発明を実施するための形態】
【００１５】
本発明の実施形態の生体信号解析装置に関し以下図面を用いて説明するが、本発明の趣旨を越えない限り、何ら本実施形態に限定されるものではない。なお、各図中、同一又は相当する部分には同一の符号を付しており、その重複説明は適宜に簡略化乃至省略する。
【００１６】
また、本実施形態においては、生体信号として、人の指先をイメージセンサ等で撮影することで血管の容積変化を捉えた容積脈波を取得しているが、指先に限らず、顔を撮影してもよい。また、イメージセンサに限らず、光センサによって容積脈波を取得してもよい。
【００１７】
さらに、本実施形態における生体信号解析処理機能を備えた生体信号解析装置として、スマートフォンを用いて説明するが、スマートフォンに限らず、情報処理機能を備えた携帯電話やＰＤＡ、タブレットＰＣ、ノートＰＣ等の情報処理装置であってもよい。
【００１８】
また、本実施形態においては、生体信号解析装置としてのスマートフォンが容積脈波を取得する生体信号取得手段としてのセンサ等を備えるものとして説明する。他方、センサ等の生体信号取得装置をスマートフォン等とは別体とし、生体信号取得装置により取得された容積脈波をスマートフォンが受け付ける受付手段としてのインタフェースを備えるシステム構成であってもよい。この場合、生体信号取得装置としてのセンサ等を備えた、例えば体組成計等を用いてもよい。
【００１９】
本実施形態の生体信号解析装置のハードウェア構成の概略について図１を参照して説明する。本実施形態の生体信号解析装置１は、撮影部１１と、ＣＰＵ１２と、ＲＡＭ１３と、ＲＯＭ１４と、ＨＤＤ１５と、表示部１６と、通信部１７を備える。
【００２０】
撮影部１１は、例えば容積脈波を取得するカメラ等のイメージセンサである。なお、撮影部１１に代えて、上述した光センサを用いて容積脈波を取得してもよい。
【００２１】
ＣＰＵ１２は、後述する本実施形態における生体信号解析処理を実行する制御手段である。
【００２２】
ＲＡＭ１３は、後述する本実施形態における生体信号解析処理のための処理データを一時的に保存するため等に活用される。ＲＯＭ１４は、例えばＣＰＵ１２において実行されるプログラム等を格納する。ＨＤＤ１５は、例えば撮影部１１を用いて取得された容積脈波等を保存する大容量記憶装置である。なお、ＨＤＤ１５に代えて、ＤＶＤ等の記憶媒体を読み込むＤＶＤ−ＲＯＭ等のデバイスであってもよい。また、ＵＳＢ等の外部記憶媒体であってもよい。
【００２３】
さらに、取得された容積脈波等を後述する通信部１７によってネットワークを介して所定のネットワークサーバへ送信し、当該サーバに容積脈波を格納し、必要なときに、当該サーバから容積脈波を取得することとしてもよい。
【００２４】
表示部１６は、例えば本実施形態における生体信号解析処理結果等を表示するディスプレイ等である。通信部１７は、外部との間でデータの送受信を行うＬＡＮやＷＡＮ等の通信インタフェースである。なお、通信は有線でも無線でも何れであってもよい。
【００２５】
なお、図１は、本実施形態の生体信号解析装置１が備える構成として、最小限の構成を示したものであるが、図示されないハードウェアとして、例えばスピーカ等の音声出力部等を備えていてもよい。また、本実施形態において、例えば表示部１６や通信部１７は必須の構成でなくてもよい。例えば、処理結果を表示部１６によって表示することに代えて、音声出力部を用いて音声で通知する場合、表示部１６はなくてもよい。
【００２６】
次に、本実施形態における生体信号解析装置１の機能ブロックについて図２を参照して説明する。本実施形態の生体信号解析装置１は、機能ブロックとして、生体信号取得部１０１と、周波数成分算出部１０２と、誤差補正部１０３と、自律神経活性度算出部１０４を備える。これらの機能ブロックは、図１に示したＣＰＵ１２によって実現される。
【００２７】
生体信号取得部１０１は、予め定められた周波数を用いて生体信号を取得する生体信号取得手段であり、図１に示したイメージセンサ等の撮影部１１に相当する。予め定められた周波数とは、例えば撮影部１１による撮影結果に含まれるＧプレーンの一定の領域の平均値をサンプリングして時系列変化としての容積脈波を取得するためのサンプリング周波数ｆｓである。
【００２８】
なお、生体信号取得部１０１は、取得した容積脈波のピーク鮮鋭化のため容積脈波の２次微分派、すなわち加速度脈波を算出してもよい。
【００２９】
周波数成分算出部１０２は、生体信号取得部１０１により取得した生体信号から低周波成分（以下、「ＬＦ」という。）と高周波成分（以下、「ＨＦ」という。）を算出する周波数成分算出手段である。例えば、周波数成分算出部１０２は、生体信号取得部１０１により取得した容積脈波の脈拍変動から周波数解析を行い、ＬＦ、ＨＦを算出する。なお、ＬＦは交感神経活性度を示し、ＨＦは副交感神経活性度を示すものである。
【００３０】
誤差補正部１０３は、周波数成分算出部１０２により算出されたＬＦ及びＨＦの誤差を、サンプリング周波数ｆｓを用いて補正する誤差補正手段である。誤差補正処理の詳細については後述する。
【００３１】
自律神経活性度算出部１０４は、誤差補正部１０３により誤差が補正されたＬＦ及びＨＦに基づいて自律神経活性度（以下、「ＡＮＦ」という。）を算出する自律神経活性度算出手段である。自律神経活性度ＡＮＦはＬＦ／ＨＦにより算出することができる。
【００３２】
本実施形態における脈拍変動の算出について図３を参照して説明する。なお、図３（ｂ）は図３（ｃ）の破線Ａで囲まれた部分を拡大したものである。脈拍変動の算出は、次のようにして行う。まず、図３（ａ）に示すように、例えばＲ１波と次のＲ２波との間隔をとってＲ-Ｒ間隔を測定する。次に、図３（ｂ）及び図３（ｃ）に示すように、測定したＲ-Ｒ間隔データを後方のＲ波の時間的位置にプロットし、これを補間する。その後、図３（ｂ）の実線で示すように、等間隔で再サンプリングしたデータを作成する。
【００３３】
次に、本実施形態における周波数解析の一例について図４を参照して説明する。図４は、横軸に周波数φ［Ｈｚ］をとり、縦軸にパワースペクトル密度［ｍｓｅｃ2］をとったグラフである。図４に示すように、ＬＦは、０．０４〜０．１５Ｈｚのパワースペクトル成分の積分値であり、ＨＦは、０．１５〜０．４０Ｈｚのパワースペクトル成分の積分値である。
【００３４】
本実施形態における誤差補正処理の詳細について説明する。誤差補正部１０３は、周波数成分算出部１０２により算出されたＬＦからＬＦの誤差成分を差し引くことによりＬＦの誤差を補正する。補正後のＬＦ成分を「ＬＦｔ」、周波数成分算出部１０２により算出されたＬＦ成分を「ＬＦｍ」、ＬＦの誤差成分を「ＬＦｄｄ」とすると、以下のように定義される。
ＬＦｔ＝ＬＦｍ−ＬＦｄｄ＝ＬＦｔ−３４７０／ｆｓ2・・・［１］
【００３５】
また、誤差補正部１０３は、周波数成分算出部１０２により算出されたＨＦからＨＦの誤差成分を差し引くことによりＨＦの誤差を補正する。補正後のＨＦ成分を「ＨＦｔ」、周波数成分算出部１０２により算出されたＨＦ成分を「ＨＦｍ」、ＨＦの誤差成分を「ＨＦｄｄ」とすると、以下のように定義される。
ＨＦｔ＝ＨＦｍ−ＨＦｄｄ＝ＨＦｔ−４７５３５／ｆｓ2・・・［２］
【００３６】
上記のように、ＨＦｄｄ及びＬＦｄｄは、サンプリング周波数ｆｓが一意に決まるため、定数となることが分かる。
【００３７】
ＨＦｄｄ及びＬＦｄｄの算出について図５を参照して説明する。ここで、サンプリング周波数ｆｓでスペクトル分析を行った各スペクトル（φ）の成分は以下のように表記することができる。
Ｓｍ（φ）＝Ｓｔ（φ）＋Ｓｄｄ（φ）
【００３８】
Ｓｍ（φ）はサンプリング周波数ｆｓで観測されるスペクトル（φ）の振幅成分であり、Ｓｔ（φ）は真のスペクトル（φ）の振幅成分であり、Ｓｄｄ（φ）はサンプリング周波数ｆｓによるスペクトル（φ）の測定誤差成分である。Ｓｄｄ（φ）は、サンプリング周波数ｆｓにより一意に決まり、次式のように定義される。
Ｓｄｄ（φ）＝１／６ｆｓ2［１−ｃｏｓ（２πφ）］
【００３９】
ＨＦｄｄ及びＬＦｄｄの算出について説明する。ＬＦ成分は、周波数φが０．０４から０．１５までの振幅成分の積算値である。振幅成分の積算値をＬＦｍ、真のＬＦの振幅成分をＬＦｔ、サンプリング周波数ｆｓによるＬＦスペクトル領域の測定誤差成分をＬＦｄｄとすると、以下のように定義できる。
ＬＦｍ＝ＬＦｔ＋ＬＦｄｄ
【００４０】
他方、ＨＦ成分は、周波数φが０．１５から０．４０までの振幅成分の積算値である。振幅成分の積算値をＨＦｍ、真のＨＦの振幅成分をＨＦｔ、サンプリング周波数ｆｓによるＨＦスペクトル領域の測定誤差成分をＨＦｄｄとすると、以下のように定義できる。
ＨＦｍ＝ＨＦｔ＋ＨＦｄｄ
【００４１】
ＬＦｄｄはＳｄｄ（φ）を用いて以下のように定義できる。ここでは、φは０．０４〜０．１５である。
ＬＦｄｄ＝∫Ｓｄｄ（φ）＝３４７０／ｆｓ2
【００４２】
他方、ＨＦｄｄはＳｄｄ（φ）を用いて以下のように定義できる。ここでは、φは０．１５〜０．４０である。
ＨＦｄｄ＝∫Ｓｄｄ（φ）＝４７５３５／ｆｓ2
【００４３】
すなわち、上述した式［１］及び式［２］で示したように、補正後のＬＦ成分ＬＦｔ、補正後のＨＦ成分ＨＦｔを算出することができる。
【００４４】
上述を踏まえ、本実施形態の生体信号解析処理手順について図６を参照して説明する。まず、不図示の入力部を介し、初期設定として、検査が開始される前等に、患者情報が生体信号解析装置としてのスマートフォン等に入力される（ステップＳ１）。患者情報とは、例えば、被測定者の氏名、年齢、ＩＤ番号、性別、糖尿病や血管障害等の既往歴等である。
【００４５】
次に、生体信号取得部１０１は、指先等から脈拍信号として容積脈波を取得する（ステップＳ２）。
【００４６】
そして、周波数成分算出部１０２は、生体信号取得部１０１により取得された容積脈波からＬＦ及びＨＦを算出する（ステップＳ３）。本ステップの詳細については後述する。
【００４７】
誤差補正部１０３は、サンプリング周波数ｆｓを用いてＬＦ及びＨＦの誤差補正を実行する（ステップＳ４）。
【００４８】
その後、自律神経活性度算出部１０４は、誤差補正部１０３により補正されたＬＦ、ＨＦから自律神経活性度ＡＮＦを算出する（ステップＳ５）。
【００４９】
そして、自律神経活性度算出部１０４により算出されたＡＮＦが表示部１６に表示され、あるいは通信部１７により外部装置に送信され、あるいはＨＤＤ１５等の記憶部に記憶される（ステップＳ６）。
【００５０】
ステップＳ３の周波数成分算出部１０２によるＬＦ及びＨＦの算出手順の詳細について図７を参照して説明する。周波数成分算出部１０２は、図３を用いて説明したように脈拍変動を算出する（ステップＳ３１）。
【００５１】
次に、周波数成分算出部１０２は、図４を用いて説明したように周波数解析処理を行い（ステップＳ３２）、ＬＦ算出処理（ステップＳ３３）、ＨＦ算出処理（ステップＳ３４）を行う。
【００５２】
本実施形態を適用した一例について図８を参照して説明する。ここでは、例えばカメラのサンプリング周波数ｆｓが６０Ｈｚであるスマートフォンでカメラの前に指先を当てて自律神経活性度を算出する例を用いることとする。図８に示すように、スマートフォン２００の裏面に設けられたカメラ２１０によりユーザの指先を撮影し、スマートフォン２００による解析処理結果２２１がスマートフォン２００の表面に設けられたディスプレイ２２０に表示される。
【００５３】
そうすると、ＬＦｄｄ及びＨＦｄｄは以下のように算出される。
ＬＦｄｄ＝３４７０／６０2＝０．９６３［ｍｓｅｃ2］
ＨＦｄｄ＝４７５３５／６０2＝１３．２０４［ｍｓｅｃ2］
【００５４】
そして、周波数成分算出部１０２により、ＬＦｍ及びＨＦｍが以下のように算出されたとする。
ＬＦｍ＝８２４．８８５［ｍｓｅｃ2］
ＨＦｍ＝５３０．９５６［ｍｓｅｃ2］
【００５５】
そうすると、補正後のＬＦ成分であるＬＦｔ及び補正後のＨＦ成分であるＨＦｔは上述した式［１］及び式［２］を用いて以下のように算出される。
ＬＦｔ＝８２４．８８５−０．９６３＝８２３．９２９［ｍｓｅｃ2］
ＨＦｔ＝５３０．９５６−１３．２０４＝５１７．７５２［ｍｓｅｃ2］
【００５６】
その結果、補正後の自律神経活性度ＡＮＦは、以下のように算出される。
ＡＮＦ＝ＬＦｔ／ＨＦｔ＝１．５９１
そして、例えば、算出されたＡＮＦが解析処理結果２２１として、スマートフォン２００のディスプレイ２２０上に表示される。
【００５７】
なお、上述する各実施の形態は、本発明の好適な実施の形態であり、本発明の要旨を逸脱しない範囲内において種々変更実施が可能である。例えば、上述した本実施形態の生体信号解析装置における各処理を、ハードウェア、又は、ソフトウェア、あるいは、両者の複合構成を用いて実行することも可能である。
【００５８】
なお、ソフトウェアを用いて処理を実行する場合には、処理シーケンスを記録したプログラムを、専用のハードウェアに組み込まれているコンピュータ内のメモリにインストールして実行させることが可能である。あるいは、各種処理が実行可能な汎用コンピュータにプログラムをインストールして実行させることが可能である。
【符号の説明】
【００５９】
１    生体信号解析装置
１１    撮影部
１２    ＣＰＵ
１３    ＲＡＭ
１４    ＲＯＭ
１５    ＨＤＤ
１６    表示部
１７    通信部
１０１    生体信号取得部
１０２    周波数成分算出部
１０３    誤差補正部
１０４    自律神経活性度算出部
２００    スマートフォン
２１０    カメラ
２２０    ディスプレイ
２２１    解析処理結果
【先行技術文献】
【特許文献】
【００６０】
【特許文献１】特開２０１０−５１８２２号公報
【特許文献２】特開２０００−１６６８８９号公報
【非特許文献】
【００６１】
【非特許文献１】Christopher G. Scully, Student Member, IEEE, Jinseok Lee. Joseph Meyer. Alexander M. Gorbach, Domhnull Granquist-Fraser, Member, IEEE, Yitzhak Mendelson, Member, IEEE, and Ki H. Chon, Senior Member, IEEE "Physiological Parameter Monitoring from Optical Recordings with a Mobile Phone",［online］, IEEE Trans Biomed Eng. (2012 February); 59(2): 303-306,［平成26年8月25日検索］、インターネット〈URL:http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3476722/pdf/nihms410059.pdf〉
【非特許文献２】Task Force of The European Society of Cardiology and The North AmericanSociety of Pacing and Electrophysiology (Membership of the Task Force listed in the Appendix) “Heart rate variability: Standards of measurement, physiological interpretation, and clinical use",［online］, European Heart Journal (1996) 17, 354-381,［平成26年8月25日検索］、インターネット〈URL: http://eurheartj.oxfordjournals.org/content/17/3/354.full.pdf+html〉
【図１】
2016043191.tif 000003
【図２】
2016043191.tif 000004
【図３】
2016043191.tif 000005
【図４】
2016043191.tif 000006
【図５】
2016043191.tif 000007
【図６】
2016043191.tif 000008
【図７】
2016043191.tif 000009
【図８】
2016043191.tif 000010