血圧測定装置及びプログラム並びに記録媒体
【課題】簡易な方法で高い測定精度が得られる血圧測定装置を提供すること。
【解決手段】ステップ140では、心電及び脈波の同時計測を行う。なお、脈波しか使用しない場合は、脈波のみの計測を行う。ステップ150では、測定によって得られた心電信号及び脈波信号の解析を行って、血圧の算出に用いたり容積AIの補正に用いる特徴量などを算出する。ステップ160では、心電信号及び脈波信号の妥当性を確認する。ステップ170では、容積AIの補正を行う。ステップ180では、例えば前記式(1)を利用し、各ステップの処理にて得られた脈波伝播時間(PTT)と補正容積AIを用いて、血圧を算出する。ステップ190では、算出した血圧を、ディスプレイ15に表示したり、スピーカ17で報知する。
【解決手段】ステップ140では、心電及び脈波の同時計測を行う。なお、脈波しか使用しない場合は、脈波のみの計測を行う。ステップ150では、測定によって得られた心電信号及び脈波信号の解析を行って、血圧の算出に用いたり容積AIの補正に用いる特徴量などを算出する。ステップ160では、心電信号及び脈波信号の妥当性を確認する。ステップ170では、容積AIの補正を行う。ステップ180では、例えば前記式(1)を利用し、各ステップの処理にて得られた脈波伝播時間(PTT)と補正容積AIを用いて、血圧を算出する。ステップ190では、算出した血圧を、ディスプレイ15に表示したり、スピーカ17で報知する。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、例えば車両に搭載されて、カフ等を用いることなくドライバなどの血圧の測定が可能な血圧測定装置及びプログラム並びに記録媒体に関する。
【背景技術】
【0002】
従来より、血圧を測定する場合には、聴診法においてもオシロメトリック法においても、カフによる締め付けが必要であり、装置が大型化したり、被測定者に苦痛があったり、連続的な血圧測定ができないなどの問題があった。
【0003】
これに対して、近年では、脈波伝播速度や脈波特徴量を用いて血圧を算出する測定法が検討されている(特許文献1〜9参照)。例えば、容積脈波であれば、光を用いて測定できるため、装置の小型化を実現できるととともに、被測定者の苦痛を解消でき、更には、連続的な血圧測定が可能になる。
【特許文献1】特開平10−295656号公報
【特許文献2】特開平10−295657号公報
【特許文献3】特開平11−318837号公報
【特許文献4】特表2001−504362号公報
【特許文献5】特表2003−555号公報
【特許文献6】特表2003−527149号公報
【特許文献7】特開2001−8907号公報
【特許文献8】特開2006−263354号公報
【特許文献9】特開2006−6897号公報
【発明の開示】
【発明が解決しようとする課題】
【0004】
しかしながら、いずれの引用文献に記載の技術の場合も、基本的に、測定精度が足りないという問題があった。また、この測定精度を高めるためには、個人適合を行うことが考えられるが、その場合には、測定機器に体重などの個人特徴量(個人特徴情報)を自分で入力する必要があり、手間がかかるという問題もあった。
【0005】
本発明は、上記点に鑑み、簡易な方法で高い測定精度が得られる血圧測定装置及びプログラム並びに記録媒体を提供することを目的とする。
【課題を解決するための手段】
【0006】
(1)請求項1の発明は、生体から検出された心電信号及び脈波信号のうち、少なくとも前記脈波信号から得られる生体情報を用いて、血圧を算出する血圧測定装置であって、中枢血管の血管弾性を表す中枢血管特徴量と、末梢血管の血管弾性を表す末梢血管特徴量とを用いて、血圧を算出することを特徴とする。
【0007】
本発明では、カフ等を用いることなく、心電信号や脈波信号などを用いて血圧を求めるので、測定者の負担が少ないという利点がある。
特に、本発明では、血圧を算出する際には、中枢血管の血管弾性を表す中枢血管特徴量と、末梢血管の血管弾性を表す末梢血管特徴量とを用いるので、従来の様に、体重等をマニュアルで入力しなくとも、高い測定精度が得られる。
【0008】
つまり、後述する実験例からも明らかな様に、血圧は大きな血管である中枢血管の弾性と小さな血管である末梢血管の弾性に大きく影響されるので、各血管の弾性を表す指標である中枢血管特徴量(中枢血管の血管弾性に対応した値(相関関係がある値))と末梢血管特徴量(末梢血管の血管弾性に対応した値(相関関係がある値))とを用いることにより、血管全体の状態(弾性の状態)を精度良く把握でき、これらを用いることで血圧の推定精度を向上できる。
【0009】
これにより、測定精度が足りないという課題を改善でき、カフによる締め付けのない高精度な血圧計を提供可能となる。
なお、ここで中枢血管とは、大動脈などの太い動脈を表し、末梢血管とは、上腕、トウ骨動脈などの動脈、細動脈などを意味する。
【0010】
さらに、中枢血管特徴量、末梢血管特徴量とも、相関があり、完全に切り分けられるものではなく、中枢血管の情報を多く含む特徴量、末梢血管の情報を多く含む特徴量という意味である。
【0011】
(2)請求項2の発明では、前記脈波信号は、容積脈波を表す信号であることを特徴とする。
本発明は、脈波信号を例示したものである。容積脈波は、血流量等に応じて変化する血管の容積の変化を示す脈波であり、通常は光学式の脈波センサを用いて測定することができ、そのため簡易な計測が可能となる。なお、容積脈波を補正することにより、圧脈波(血管の圧力の変化に対応する脈波)に近い精度が得られる。
【0012】
(3)請求項3の発明では、前記中枢血管特徴量は、前記心電信号と前記脈波信号との時間差(PTT)、又は脈波波形により得られる進行波から反射波までの時間、若しくはその値(即ち、PTT、又は進行波から反射波までの時間)を用いて求められるもの(例えば前記時間を速度で示した値など)であることを特徴とする。
【0013】
本発明は、中枢血管特徴量を例示したものである。
ここで、心電信号と脈波信号との時間差は、図1に示す様に、脈波伝播時間(PTT)であり、心電図のR波と手指で測定した脈波との時差として求めることができる。
【0014】
また、脈波波形により得られる進行波(心臓の収縮時に発生する駆出波)から反射波(駆出波が末梢血管や動脈の分岐部で発生する反射波)までの時間は、図2に示す様に、脈波における反射波(P2)と進行波(P1)との時間差ΔTとして求めることができる。
【0015】
(4)請求項4の発明では、前記末梢血管特徴量は、脈波のAI、又は加速度脈波の特徴量のd、若しくはその値(即ち脈波のAI、又はd)を用いて求められるものであることを特徴とする。
【0016】
本発明は、末梢血管特徴量を例示したものである。
ここで、脈波のAI(Augmentation Index)は、図2に示す様に、脈波における反射波(P2)と進行波(P1)の大きさの比(P2/P1)であり、脈波(容積脈波又は圧脈波)から求めることができる。
【0017】
また、加速度脈波の特徴量としては、時間軸に沿った波のピークを示すa波(収縮初期陽性波)、b波(収縮初期陰性波)、c波(収縮中期再上昇波)、d波(収縮後期再下降波)、e波(拡張初期陽性波)、f波(拡張初期陰性波)があり、この加速度脈波の特徴量a〜fは、各波のピークの大きさ(振幅)を示している。
【0018】
なお、脈波信号を一階微分することにより、速度脈波を求めることができ、脈波信号を二階微分することにより、加速度脈波を求めることができる。
(5)請求項5の発明では、前記容積脈波から得られる容積AIを、前記加速度脈波の特徴量のd、e、fのうち1種を用いて補正することを特徴とする。
【0019】
本発明者等の研究によれば、容積AI(容積脈波から求めたAI)を、加速度脈波の特徴量のd、e、fのうち1種を用いて補正することにより、圧AI(圧脈波から求めたAI)に近い値が得られることが確認されている。
【0020】
よって、この容積AIを補正した補正容積AI(単に補正AIとも記す)を用いて 血圧を算出することにより、精度良く血圧を求めることができる。
(6)請求項6の発明では、前記容積脈波から得られる容積AIを、前記加速度脈波の特徴量のe−f点間における脈波の傾きと、心拍数とを用いて補正することを特徴とする。 本発明者等の研究によれば、容積AIを、加速度脈波の特徴量のe−f点間における脈波の傾きと心拍数(HR)を用いて補正することにより、圧AIに近い値が得られることが確認されている。
【0021】
よって、この容積AIを補正した補正容積AIを用いて 血圧を算出することにより、精度良く血圧を求めることができる。
ここで、加速度脈波の特徴量のe−f点間における脈波の傾きとは、大動脈弁閉鎖によるノッチ(notch)からノッチ直後の後隆波のピークまでの傾きであり、図2に示す様に、速度脈波の特徴量のe1(振幅)に相当する値である。なお、このe1とは、加速度脈波e−f点の間にある速度脈波のピーク値である。また、後隆波にピークが見られない場合には、ノッチから所定時間経過後までの傾きを採用してよい。
【0022】
尚、心拍数の代わりに心拍間隔RRIを用いてよい。
(7)請求項7の発明では、前記心電信号及び前記脈波信号以外の生体情報も加味して、血圧を算出することを特徴とする。
【0023】
これにより、一層血圧の算出精度が向上する。
(8)請求項8の発明では、前記心電信号及び前記脈波信号以外の生体情報は、身体特徴量及び生体のインピーダンスのうち少なくとも1種であることを特徴とする。
【0024】
本発明は、脈波信号以外の生体情報を例示したものである。例えば身長、体重、年齢、性別の様な身体特徴情報は、血流量(従って血圧)に相関関係があるので、この身体特徴情報を用いることにより、精度良く血圧を算出することができる。また、インピーダンスは、体重や体脂肪(従って血流量)に相関関係があるので、精度の高い血圧測定が可能である。
【0025】
(9)請求項9の発明は、生体から検出された心電信号及び脈波信号のうち、少なくとも前記脈波信号から得られる生体情報を用いて、血圧を算出する血圧測定装置であって、脈波波形の立ち下がりの状態を示す立ち下がり情報を用いて、血圧を算出することを特徴とする。
【0026】
脈波波形の立ち下がり情報は、血管弾性の情報を含むものである。つまり、脈波波形の立ち下がりの部分は、脈波の立ち上がりの部分と異なり、心拍出波形の影響を受けづらいため、血管弾性情報を心臓の影響を少なくして抽出できる。この血管弾性情報を用いることにより、血圧を精度良く算出できる。
【0027】
ここで、立ち下がりの範囲とは、大動脈弁閉鎖によるノッチ直後の脈波(後隆波)のピークからの立ち下がりを示しており、ノッチ直後にピークが無い場合には、例えば4階微分のゼロクロス点に対応する脈波の位置からの立ち下がりを採用できる。
【0028】
なお、立ち下がり情報としては、複数の情報を用いると、血圧の算出精度が向上するので、好適である。
(10)請求項10の発明では、前記立ち下がり情報として、脈波波形の立ち下がりの最大傾きを用いることを特徴とする。
【0029】
本発明は、好ましい立ち下がり情報を例示したものである。本発明者等の研究によれば、脈波波形の立ち下がりの最大傾きを用いることにより、血圧を精度良く算出できることが確認されている。
【0030】
(11)請求項11の発明では、前記立ち下がり情報として、加速度脈波の特徴量のe点以降における脈波波形の最大値の8〜2割の範囲の傾きを用いることを特徴とする。
本発明は、好ましい立ち下がり情報を例示したものである。本発明者等の研究によれば、図3に示す様に、加速度脈波の特徴量のe点以降における脈波波形の最大値の8〜2割の範囲の傾きを用いることにより、血圧を精度良く算出できることが確認されている。
【0031】
(12)請求項12の発明では、(前記請求項9〜11のいずれかに記載の血圧測定装置に加えて)更に、前記請求項1〜8のいずれかに記載の血圧測定装置の構成を備えたことを特徴とする。
【0032】
これにより、血圧の算出精度を一層高めることができる。
(13)請求項13の発明は、前記請求項1〜12のいずれかに記載の血圧測定装置の機能を実現させるためのプログラムである。
【0033】
よって、マイクロコンピュータ等でこのプログラムを実施することにより、上述した生体情報などを用い、所定の演算式を使用して、血圧を算出することが可能である。
(14)請求項14の発明は、前記請求項13に記載のプログラムを記録したコンピュータ読取可能な記録媒体である。
【0034】
よって、この記録媒体に記録したプログラムを用いて、上述したように、血圧を算出することができる。
【発明を実施するための最良の形態】
【0035】
次に、本発明の最良の形態(実施形態)について説明する。
[第1実施形態]
本実施形態では、自動車に搭載されてドライバの血圧を測定する血圧測定装置について説明する。
【0036】
a)まず、本実施形態の血圧測定装置のシステム構成を説明する。
図4に示す様に、本実施形態では、自動車に搭載された血圧測定装置1と、ドライバ等に情報の報知を行う報知装置3と、手動にてデータの入力を行うマニュアル入力部5と、ステアリング7に取り付けられた脈波センサ9と、ステアリング7に取り付けられた一対の電極11、13とを備えている。
【0037】
前記血圧測定装置1は、周知のマイクロコンピュータを中心とした電子制御装置であり、脈波センサ9や電極11、13等からの信号に基づいて、血圧の算出や報知装置3の制御等を行う。
【0038】
前記報知装置3は、血圧などを表示する液晶等のディスプレイ15や、その内容を音声等で出力するスピーカ17により構成される。
前記マニュアル入力部5は、マニュアルにて、体重、身長、年齢、性別等の個人別の身体特徴情報が入力可能な例えばキーボードやテンキーやリモコンなどの入力装置である。なお、ディスプレイ15の表示画面をタッチパネルとして、データを入力するようにしてもよい。
【0039】
前記脈波センサ9は、周知の発光素子(LED)や受光素子(PD)を備えた光学式のセンサであり、例えばドライバの指先等に光を照射し、その反射波を利用して脈波(容積脈波)を検出することができる。従って、後述する血圧の算出に用いる脈波信号は、この脈波センサ9から得ることができる。
【0040】
前記一対の電極11、13は、ドライバの左右の手にそれぞれ接触するように、ステアリング9の左右に配置されている。この電極11、13は、心電信号を得るために電圧を印加する心電計の電極として用いられる。ここでは、両電極11、13と血圧測定装置1が、心電計の機能を果たす。従って、後述する血圧の算出に用いる心電信号は、この電極11、13を利用して得ることができる。
【0041】
なお、一対の電極11、13間に電流を流しそれによってインピーダンスが得られる。
b)次に、本実施形態の血圧測定装置等の機能をブロック図にて更に詳しく説明する。
図5に示す様に、本実施形態の血圧測定装置1には、インピーダンス計測部21と、心電/インピーダンス計測切替部23と、心電信号取得部25と、脈波信号取得部27と、心電信号解析部31と、心電及び脈波信号解析部33と、脈波信号解析部35と、血圧演算部37とを備えている。
【0042】
このうち、インピーダンス計測部21は、電極11、13間に微少電流(例えば高周波及び低周波の複数の周波数の交流電流)を流してインピーダンスを測定するものである。 心電信号取得部25は、心臓の活動に伴う電気的興奮を、電極11、13間の電位差(心電信号)として計測するものである。
【0043】
心電/インピーダンス計測切替部23は、インピーダンス計測部21と心電信号取得部25との機能を自動的に切り替えるものである。
脈波信号取得部27は脈波センサ9を駆動して脈波信号を取得するものである。
【0044】
心電信号解析部31は、心電信号を解析して、例えば周知のQT信号等を算出するものである。
心電及び脈波信号解析部33は、前記図1に示す様に、心電信号と脈波信号を用いて、心電信号に対する脈波信号の遅れ時間である脈波伝播時間(PTT)を求めるものである。
【0045】
脈波信号解析部35は、前記図2に示す様に、脈波信号を解析して、1階微分(速度脈波)、2階微分(加速度脈波)、3階微分、4階微分を行うとともに、各微分における特徴量(例えば速度脈波のa1〜k1、加速度脈波のa〜f等)の算出、AI(容積AI)の算出、容積AIの補正などを行うものである。
【0046】
血圧演算部37は、PTT、速度脈波や加速度脈波などの特徴量、容積AI等を用い、所定の演算式に基づいて血圧を算出するものである。
c)次に、血圧の算出に用いる演算式について説明する。
【0047】
(i)血圧の算出式
ここでは、中枢血管特徴量(PTT、脈波面積)と末梢血管特徴量(容積AI)とを用いた演算式、例えば下記式(1)〜(5)のいずれかの演算式を用いて、血圧(推定血圧:EBP)を求める。
・PTTを使用する場合
EBP=α・PTT+β・容積AI+γ ・・(1)
EBP=α・PTT+β・容積AI+γ・(脈波特徴量)+・・+ ・・(2)
EBP=α・PTT+β・容積AI+γ・d+δ・e+ε・f
+ζ・HR+・・+ ・・(3)
・PTTを使用しない場合
EBP=α・脈波面積+β・容積AI+γ ・・(4)
・より精度の良い推定式
EBP=α・PTT+β・容積AI+γ・(脈波特徴量)+δ・(身体特徴量)
+ε・(インピーダンス)+ζ ・・(5)
前記各式で、α、β、γ、δ、ε、ζは係数である。脈波特徴量は、例えばa1、b1など(図2参照)である。HRは脈拍であり、脈拍に代えて脈拍間隔(RRI)を用いることができる。脈波面積とは、脈波の基線をゼロとした場合の波形に囲まれる部分の面積である。身体特徴量は、身長、体重、性別、年齢などである。
【0048】
以下、各式について説明する。
前記式(1)では、推定血圧の算出に、PTTと容積AIを用いる。このうち、PTTは中枢血管の血管弾性を示すパラメータであり、容積AIは末梢血管の血管弾性を示すパラメータであるので、両パラメータを用いることにより、中枢、末梢系の情報を取り入れられているため、精度良く血圧を推定することができる。
【0049】
前記式(2)では、推定血圧の算出に、PTTと容積AIに加えて、脈波特徴量を用いる。よって、前記式(1)よりも精度良く血圧を推定することができる。ここで、脈波特徴量を加味すると推定精度が向上するのは、特徴量の非線形性の効果、及び心臓の拍出の影響等を取り入れられているからである。
【0050】
前記式(3)では、推定血圧の算出に、PTTと容積AIに加えて、加速度脈波のd、e、f、及びHRを用いる。よって、前記式(1)よりも精度良く血圧を推定することができる。ここで、加速度脈波のd、e、fとHRを加味すると推定精度が向上するのは、後述の容積AIを圧AIに補正する効果、及び心臓の拍出の影響の取り入れがあるからである。
【0051】
前記式(4)では、推定血圧の算出に、脈波面積(S)と容積AIを用いている。このうち、Sは中枢血管の血管弾性を示すパラメータであると考えられ(実験により推定)、容積AIは末梢血管の血管弾性を示すパラメータであるので、両パラメータを用いることにより、マニュアルにて体重を入力しなくても、精度良く血圧を推定することができる。
【0052】
前記式(5)では、推定血圧の算出に、PTTと容積AIに加えて、脈波特徴量、身体特徴量、インピーダンスを用いている。よって、前記式(2)よりも精度良く血圧を推定することができる。なお、身長や体重などの身体特徴量は、血流量と相関関係があり、また、インピーダンスは体重や体脂肪と相関関係がある(よって血流量とも相関関係がある)ので、これらを加味することにより、血圧の推定精度が向上すると考えられる。
【0053】
(ii)容積AIの補正式
容積AIはその測定方法の特性上、圧脈波により求めた圧AIよりは精度が低い傾向にある。従って、ここでは、容積AIを下記式(6)又は(7)を用いて補正する。
補正容積AI=容積AI+α・HR+β・ノッチ−後隆波傾き ・・・(6)
補正容積AI=容積AI+α・d+β・e+γ・f ・・・(7)
前記各式で、α、β、γは係数である。また、HRは心拍数である(心拍間隔RRIを用いてもよい)。ノッチ−後隆波傾きとは、上述したノッチ(notch)からノッチ直後の後隆波のピークまでの傾き(即ち加速度脈波の特徴量のe−f点間における脈波の傾き=速度脈波の特徴量のg1)である。d、e、fは、上述した加速度脈波の特徴量である。
【0054】
以下、各式について説明する。
前記式(6)では、容積AIの補正に、HRとノッチ−後隆波傾きを用いているので、圧AIに近い容積AIとすることができる。なお、HRとノッチ−後隆波傾きを用いると良いことは、実験的に検証した。
【0055】
前記式(7)では、容積AIの補正に、加速度脈波の特徴量のd、e、fを用いているので、圧AIに近い容積AIとすることができる。なお、d、e、fを用いると良いことは、実験的に検証した。
【0056】
本実施形態では、前記式(6)、(7)より補正した容積AIを、前記式(1)〜(5)に用いることにより、血圧の推定精度が一層向上する。
d)次に、血圧検出装置1にて実施される制御処理について説明する。
【0057】
図6のフローチャートに示す様に、ステップ(S)100では、身体特徴情報を利用するか否かを判定する。例えばディスプレイ15に、身体特徴情報を利用するかどうかを表示し、ドライバの入力を促すようにしても良い。ここで肯定判断されるとステップ110に進み、一方否定判断されるとステップ120に進む。
【0058】
ステップ110では、身体特徴情報の入力が選択されたので、ドライバ等により、例えばタッチパネルによって体重等の入力がされた場合には、そのデータを取得する。
ステップ120では、インピーダンスの計測を行うか否かを判定する。例えばディスプレイ15に、インピーダンスを計測するか否かを表示し、ドライバの入力を促すようにしても良い。ここで肯定判断されるとステップ130に進み、一方否定判断されるとステップ140に進む。
【0059】
ステップ130では、例えばタッチパネルによって体重等の入力がされた場合には、そのデータを取得する。
続くステップ135では、インピーダンスの計測を行う。
【0060】
ステップ140では、心電及び脈波の同時計測を行う。なお、脈波しか使用しない場合は、脈波のみの計測を行う。
続くステップ150では、測定によって得られた心電信号及び脈波信号の解析を行って、血圧の算出や容積AIの補正に用いる演算値(特徴量など)を算出する。
【0061】
例えば、前記式(1)を用いて血圧を算出する場合には、 脈波伝播時間(PTT)と容積AIが必要であるので、ここでは、心電信号と脈波信号とを用いて脈波伝播時間(PTT)を算出するとともに、P1、P2を用いて容積AIを算出する。なお、他の式を用いる場合には、各式に応じて、血圧の算出に必要な特徴量などの値を算出する。
【0062】
続くステップ160では、心電信号及び脈波信号の妥当性を確認する。つまり、心電信号や脈波信号が測定値として妥当な値であるか否かを判定する。ここで、妥当な値であると判断されるとステップ170に進み、妥当な値でないと判断されると、再度心電及び脈波の測定を行うために、ステップ140に戻る。
【0063】
続くステップ170では、容積AIの補正が必要な場合には、ここで容積AIの補正を行う。
例えば前記式(6)を用いて容積AIを補正する場合には、心拍数(HR)とノッチ−後隆波傾きとを用いて、補正容積AIを算出する。
【0064】
ステップ180では、例えば前記式(1)を利用し、上述した各ステップの処理にて得られた脈波伝播時間(PTT)と補正容積AIを用いて、血圧を算出する。
続くステップ190では、算出した血圧を、ディスプレイ15に表示したり、スピーカ17で報知し、一旦本処理を終了する。
【0065】
d)次に、本実施形態による効果を確認した実験例について説明する。
図7(a)は、横軸にトノメトリ血圧により得られたトウ骨動脈圧脈波を解析して求めた圧AIをとり、縦軸に光学式脈波センサにより求めた脈波信号から算出した容積AI(未補正)をとったものである。
【0066】
この図7(a)から、圧AIと(補正前の)容積AIとの相関係数(r)は0.82であり、圧AIと容積AIとの相関関係が小さいことが分かる。よって、容積AIを用いて血圧を算出する場合には、圧AIを用いた場合より精度が低くなる。
【0067】
図7(b)は、横軸に圧AIをとり、縦軸に前記式(6)により補正した補正容積AI(補正AI)をとったものである。
この図7(b)から、圧AIと補正容積AIとの相関係数は0.93であり、圧AIと補正容積AIとの相関関係が大きいことが分かる。よって、圧AIに代えて補正容積AIを用いることにより、精度の高い血圧を算出できることが分かる。
【0068】
図8(a)は、横軸にカフの締め付けを行うオシロメトリック式血圧計により測定した実血圧をとり、縦軸にPTTのみを用いて算出した推定血圧(即ち前記式(1)でPTT項のみを用いて推定血圧を算出したもの)をとったものである。
【0069】
この図8(a)から、実血圧と推定血圧との相関が少なく、PTTのみを用いて血圧を推定した場合には、精度の高い血圧を求めることができないことが分かる。
図8(b)は、前記式(1)を用いて推定血圧を求め、実血圧との相関をとったものである。この場合には、PTTと補正容積AIを用いて血圧を推定するので、推定血圧の精度が高いことがわかる。
【0070】
図8(c)は、前記式(1)と式(7)を用いて推定血圧を求め、実血圧との相関をとったものである。この場合には、PTTと補正容積AIとd、e、fを用いて血圧を推定するので、推定血圧の精度が一層高いことがわかる。
【0071】
図9は、前記式(5)と式(7)を用いて推定血圧を求め、実血圧との相関をとったものである。この場合には、PTTと補正容積AIと脈波特徴量とインピーダンスと身体特徴量を用いて血圧を推定しているので、推定血圧の精度がより一層高いことがわかる。
【0072】
e)この様に、本実施形態では、例えば前記式(1)及び式(6)を利用し、心電信号や脈波信号から得られた例えば脈波伝播時間(PTT)や補正容積AIを用いて血圧を算出するので、簡易な手法で精度の高い血圧測定が可能である。
【0073】
なお、血圧の算出に前記式(1)〜(5)を用いる場合は、補正しない容積AIを用いることも可能であるが、補正容積AIを用いた方が、血圧の算出精度が高まるので好適である。
[第2実施形態]
次に、第2実施形態を説明するが、前記第1実施形態と同様な内容の説明は省略する。
【0074】
本実施形態は、ハード構成は前記第1実施形態と同様であるので、使用する演算式や演算処理の手順について説明する。
a)まず、本実施形態における血圧の算出に用いる演算式について説明する。
【0075】
(i)血圧の算出式
ここでは、脈波(詳しくは後隆波)の立ち下がり情報を用いた演算式、例えば下記式(8)〜(16)のいずれかの演算式を用いて、血圧(推定血圧:EBP)を求める。
・立ち下がりの傾きのみを使用する場合
EBP=α・e1 ・・(8)
EBP=α・e1+β・f1 ・・(9)
EBP=α・(0.8−0.6)の傾き+β・(0.4−0.2) ・・(10)
EBP=α・最大傾き ・・(11)
・脈波及びPTTを使用する場合
EBP=α・e1+β・PTT ・・(12)
EBP=α・e1+β・f1+γ・PTT ・・(13)
EBP=α・e1+β・f1+γ・PTT+δ・(脈波特徴量)+・ ・・(14)
・C1、C2を使用する場合
EBP=α・PTT+β・容積AI+γ・C1+δ・C2+ε ・・(15)
・より精度の良い推定式
EBP=α・e1+β・f1+γ・PTT+δ・(脈波特徴量)
+ε・(身体特徴量)+ζ・(インピーダンス)+η ・・(16)
前記各式で、α、β、γ、δ、ε、ζ、ηは係数である。脈波特徴量は、例えばa1、d、AIである。e1は、図2に示す様に、1階微分の脈波の立ち下がりの傾きの最小値であり、f1は、2階微分のf点後での、4階微分のゼロクロスに対応する点での脈波の立ち下がりの傾斜である。身体特徴量は、身長、体重、性別、年齢などである。C1、C2は、脈波の立下りを指数関数で近似することにより求められる指数であり、指数関数での近似により、精度良く情報の抽出が可能となる(Daniel A Duprez,Determinateion of Radial Artery Pulse Wave Analysis in Asymptomatic Individuals,American Journal of Hypertension 2004;17:647−653参照)。
【0076】
以下、各式について説明する。
前記式(8)では、推定血圧の算出に、脈波の立ち下がり情報であるe1を用いるので、精度良く血圧を推定することができる。ここで、e1を用いると推定精度が高いのは、心臓系の影響の少ない血管弾性情報を用いることができているからである。
【0077】
前記式(9)では、推定血圧の算出に、脈波の立ち下がり情報であるe1とf1を用いるので、e1だけを用いる場合よりも、精度良く血圧を推定することができる。ここで、f1を加味すると測定精度が向上するのは、大動脈及び小動脈の血管情報を併せて考慮できているからである。
【0078】
前記式(10)では、推定血圧の算出に、図3に示す様に、後隆波の傾きの一部を使用する(ここでは、後隆波のピークを1とし、それを20%づつ5等分し、各区分における傾きを求める)ので、精度良く血圧を推定することができる。なお、これらの傾きの平均値を用いてもよい。
【0079】
前記式(11)では、推定血圧の算出に、後隆波の最大傾きを用いるので、精度良く血圧を推定することができる。なお、図3の様に区分した領域での傾きの最大値を、最大傾きとしてもよい。ここで、最大傾きを加味すると測定精度が向上するのは、e1とほぼ同義であり、ノイズの多い環境下では、微分して求めるよりも安定して求められるため、精度が向上するからである。
【0080】
前記式(12)では、推定血圧の算出に、脈波の立ち下がり情報であるe1とPTTとを用いるので、精度良く血圧を推定することができる。
前記式(13)では、推定血圧の算出に、脈波の立ち下がり情報であるe1、f1とPTTとを用いるので、より精度良く血圧を推定することができる。
【0081】
前記式(14)では、推定血圧の算出に、脈波の立ち下がり情報であるe1、f1とPTTと脈波特徴量を用いるので、より一層精度良く血圧を推定することができる。
前記式(15)では、推定血圧の算出に、脈波の立ち下がり情報であるC1、C2とPTTと容積脈AI(補正容積AIでも良い)を用いるので、精度良く血圧を推定することができる。
【0082】
前記式(16)では、推定血圧の算出に、脈波の立ち下がり情報であるe1、f1とPTTと脈波特徴量と身体特徴量とインピーダンスとを用いるので、一層精度良く血圧を推定することができる。
【0083】
また前記式(14)〜(16)において、脈波特徴量として、AIもしくはd等を使用すると、請求項12に示すように、請求項1〜8による血圧推定精度向上も同時に行うことになる。
【0084】
この場合は、PTT及びAI、dにより、中枢及び末梢の弾性、e1及びf1により、同様に中枢及び末梢の弾性による血圧推定精度の向上を図るが、AI、dは心拍出波形の影響も受けるため、e1、f1を併せて用いることで、中枢、末梢、心臓系を考慮した血圧推定が実施可能となり、推定精度が大きく向上する。
【0085】
b)次に、血圧検出装置にて実施される制御処理について説明する。
図10のフローチャートに示す様に、ステップ(S)200では、身体特徴情報を利用するか否かを判定する。ここで肯定判断されるとステップ210に進み、一方否定判断されるとステップ220に進む。
【0086】
ステップ210では、例えばタッチパネルによって体重等の入力がされた場合には、そのデータを取得する。
ステップ220では、インピーダンスの計測を行うか否かを判定する。ここで肯定判断されるとステップ230に進み、一方否定判断されるとステップ240に進む。
【0087】
ステップ230では、例えばタッチパネルによって体重等の入力がされた場合には、そのデータを取得する。
続くステップ235では、インピーダンスの計測を行う。
【0088】
ステップ240では、心電及び脈波の同時計測を行う。なお、脈波しか使用しない場合は、脈波のみの計測を行う。
続くステップ250では、測定によって得られた心電信号及び脈波信号の解析を行って、血圧の算出に必要な特徴量などを求める。
【0089】
続くステップ260では、心電信号及び脈波信号の妥当性を確認する。ここで、妥当な値であると判断されるとステップ270に進み、妥当な値でないと判断されると、ステップ240に戻る。
【0090】
続くステップ270では、前記各式(8)〜(16)を利用して、血圧を算出する。
続くステップ280では、算出した血圧を報知し、一旦本処理を終了する。
c)次に、本実施形態による効果を確認した実験例について説明する。
【0091】
図11(a)は、横軸に実血圧をとり、縦軸にPTTのみを用いて算出した推定血圧(即ち前記式(3)でPTTの項のみを用いて推定血圧を算出したもの)をとったものである。
【0092】
この図11(a)から、PTTのみを用いて血圧を推定した場合には、精度の高い血圧を求めることができないことが分かる。
図11(b)は、前記式(12)を用いて推定血圧を求め、実血圧との相関をとったものである。この場合には、e1とPTTと用いて血圧を推定しているので、推定血圧の精度が高いことがわかる。
【0093】
図11(c)は、前記式(13)を用いて推定血圧を求め、実血圧との相関をとったものである。この場合には、e1、f1とPTTを用いて血圧を推定しているので、推定血圧の精度が一層高いことがわかる。
【0094】
この様に、本実施形態では、脈波の立ち下がり情報を用いて血圧を推定するので、簡易な手法で精度の高い血圧測定が可能である。
尚、本発明は前記実施形態になんら限定されるものではなく、本発明を逸脱しない範囲において種々の態様で実施しうることはいうまでもない。
【0095】
例えば、上述した血圧測定装置の機能は、コンピュータのプログラムにより実行される処理により実現することができ、このプログラムは、記録媒体に記録できる。
つまり、上述した様なプログラムをコンピュータシステムにて実現する機能は、例えば、コンピュータシステム側で起動するプログラムとして備えることができる。このようなプログラムの場合、例えば、フレキシブルディスク、光磁気ディスク、CD−ROM、DVD、ハードディスク等のコンピュータ読み取り可能な記録媒体に記録し、必要に応じてコンピュータシステムにロードして起動することにより用いることができる。この他、ROMやバックアップRAM等をコンピュータ読み取り可能な記録媒体として前記プログラムを記録しておき、このROMあるいはバックアップRAM等をコンピュータシステムに組み込んで用いても良い。
【図面の簡単な説明】
【0096】
【図1】心電信号と脈波信号を示すグラフである。
【図2】脈波信号とそれを微分した信号を示すグラフである。
【図3】後隆波の傾きを示すグラフである。
【図4】第1実施形態の車室内に配置された血圧測定装置の概略構成を示す説明図である。
【図5】血圧測定装置の機能を示すブロック図である。
【図6】第1実施形態の制御処理を示すフローチャートである。
【図7】(a)圧AIと容積AIとの関係を示すグラフ、(b)圧AIと補正容積AIとの関係を示すグラフである。
【図8】カフを用いた実測値と血圧測定装置による演算値との相関関係を示すグラフである。
【図9】カフを用いた実測値と血圧測定装置による演算値との相関関係を示すグラフである。
【図10】第2実施形態の制御処理を示すフローチャートである。
【図11】カフを用いた実測値と血圧測定装置による演算値との相関関係を示すグラフである。
【符号の説明】
【0097】
1…血圧測定装置
3…報知装置
5…マニュアル入力部
7…ステアリング
9…脈波センサ
11、13…電極
【技術分野】
【0001】
本発明は、例えば車両に搭載されて、カフ等を用いることなくドライバなどの血圧の測定が可能な血圧測定装置及びプログラム並びに記録媒体に関する。
【背景技術】
【0002】
従来より、血圧を測定する場合には、聴診法においてもオシロメトリック法においても、カフによる締め付けが必要であり、装置が大型化したり、被測定者に苦痛があったり、連続的な血圧測定ができないなどの問題があった。
【0003】
これに対して、近年では、脈波伝播速度や脈波特徴量を用いて血圧を算出する測定法が検討されている(特許文献1〜9参照)。例えば、容積脈波であれば、光を用いて測定できるため、装置の小型化を実現できるととともに、被測定者の苦痛を解消でき、更には、連続的な血圧測定が可能になる。
【特許文献1】特開平10−295656号公報
【特許文献2】特開平10−295657号公報
【特許文献3】特開平11−318837号公報
【特許文献4】特表2001−504362号公報
【特許文献5】特表2003−555号公報
【特許文献6】特表2003−527149号公報
【特許文献7】特開2001−8907号公報
【特許文献8】特開2006−263354号公報
【特許文献9】特開2006−6897号公報
【発明の開示】
【発明が解決しようとする課題】
【0004】
しかしながら、いずれの引用文献に記載の技術の場合も、基本的に、測定精度が足りないという問題があった。また、この測定精度を高めるためには、個人適合を行うことが考えられるが、その場合には、測定機器に体重などの個人特徴量(個人特徴情報)を自分で入力する必要があり、手間がかかるという問題もあった。
【0005】
本発明は、上記点に鑑み、簡易な方法で高い測定精度が得られる血圧測定装置及びプログラム並びに記録媒体を提供することを目的とする。
【課題を解決するための手段】
【0006】
(1)請求項1の発明は、生体から検出された心電信号及び脈波信号のうち、少なくとも前記脈波信号から得られる生体情報を用いて、血圧を算出する血圧測定装置であって、中枢血管の血管弾性を表す中枢血管特徴量と、末梢血管の血管弾性を表す末梢血管特徴量とを用いて、血圧を算出することを特徴とする。
【0007】
本発明では、カフ等を用いることなく、心電信号や脈波信号などを用いて血圧を求めるので、測定者の負担が少ないという利点がある。
特に、本発明では、血圧を算出する際には、中枢血管の血管弾性を表す中枢血管特徴量と、末梢血管の血管弾性を表す末梢血管特徴量とを用いるので、従来の様に、体重等をマニュアルで入力しなくとも、高い測定精度が得られる。
【0008】
つまり、後述する実験例からも明らかな様に、血圧は大きな血管である中枢血管の弾性と小さな血管である末梢血管の弾性に大きく影響されるので、各血管の弾性を表す指標である中枢血管特徴量(中枢血管の血管弾性に対応した値(相関関係がある値))と末梢血管特徴量(末梢血管の血管弾性に対応した値(相関関係がある値))とを用いることにより、血管全体の状態(弾性の状態)を精度良く把握でき、これらを用いることで血圧の推定精度を向上できる。
【0009】
これにより、測定精度が足りないという課題を改善でき、カフによる締め付けのない高精度な血圧計を提供可能となる。
なお、ここで中枢血管とは、大動脈などの太い動脈を表し、末梢血管とは、上腕、トウ骨動脈などの動脈、細動脈などを意味する。
【0010】
さらに、中枢血管特徴量、末梢血管特徴量とも、相関があり、完全に切り分けられるものではなく、中枢血管の情報を多く含む特徴量、末梢血管の情報を多く含む特徴量という意味である。
【0011】
(2)請求項2の発明では、前記脈波信号は、容積脈波を表す信号であることを特徴とする。
本発明は、脈波信号を例示したものである。容積脈波は、血流量等に応じて変化する血管の容積の変化を示す脈波であり、通常は光学式の脈波センサを用いて測定することができ、そのため簡易な計測が可能となる。なお、容積脈波を補正することにより、圧脈波(血管の圧力の変化に対応する脈波)に近い精度が得られる。
【0012】
(3)請求項3の発明では、前記中枢血管特徴量は、前記心電信号と前記脈波信号との時間差(PTT)、又は脈波波形により得られる進行波から反射波までの時間、若しくはその値(即ち、PTT、又は進行波から反射波までの時間)を用いて求められるもの(例えば前記時間を速度で示した値など)であることを特徴とする。
【0013】
本発明は、中枢血管特徴量を例示したものである。
ここで、心電信号と脈波信号との時間差は、図1に示す様に、脈波伝播時間(PTT)であり、心電図のR波と手指で測定した脈波との時差として求めることができる。
【0014】
また、脈波波形により得られる進行波(心臓の収縮時に発生する駆出波)から反射波(駆出波が末梢血管や動脈の分岐部で発生する反射波)までの時間は、図2に示す様に、脈波における反射波(P2)と進行波(P1)との時間差ΔTとして求めることができる。
【0015】
(4)請求項4の発明では、前記末梢血管特徴量は、脈波のAI、又は加速度脈波の特徴量のd、若しくはその値(即ち脈波のAI、又はd)を用いて求められるものであることを特徴とする。
【0016】
本発明は、末梢血管特徴量を例示したものである。
ここで、脈波のAI(Augmentation Index)は、図2に示す様に、脈波における反射波(P2)と進行波(P1)の大きさの比(P2/P1)であり、脈波(容積脈波又は圧脈波)から求めることができる。
【0017】
また、加速度脈波の特徴量としては、時間軸に沿った波のピークを示すa波(収縮初期陽性波)、b波(収縮初期陰性波)、c波(収縮中期再上昇波)、d波(収縮後期再下降波)、e波(拡張初期陽性波)、f波(拡張初期陰性波)があり、この加速度脈波の特徴量a〜fは、各波のピークの大きさ(振幅)を示している。
【0018】
なお、脈波信号を一階微分することにより、速度脈波を求めることができ、脈波信号を二階微分することにより、加速度脈波を求めることができる。
(5)請求項5の発明では、前記容積脈波から得られる容積AIを、前記加速度脈波の特徴量のd、e、fのうち1種を用いて補正することを特徴とする。
【0019】
本発明者等の研究によれば、容積AI(容積脈波から求めたAI)を、加速度脈波の特徴量のd、e、fのうち1種を用いて補正することにより、圧AI(圧脈波から求めたAI)に近い値が得られることが確認されている。
【0020】
よって、この容積AIを補正した補正容積AI(単に補正AIとも記す)を用いて 血圧を算出することにより、精度良く血圧を求めることができる。
(6)請求項6の発明では、前記容積脈波から得られる容積AIを、前記加速度脈波の特徴量のe−f点間における脈波の傾きと、心拍数とを用いて補正することを特徴とする。 本発明者等の研究によれば、容積AIを、加速度脈波の特徴量のe−f点間における脈波の傾きと心拍数(HR)を用いて補正することにより、圧AIに近い値が得られることが確認されている。
【0021】
よって、この容積AIを補正した補正容積AIを用いて 血圧を算出することにより、精度良く血圧を求めることができる。
ここで、加速度脈波の特徴量のe−f点間における脈波の傾きとは、大動脈弁閉鎖によるノッチ(notch)からノッチ直後の後隆波のピークまでの傾きであり、図2に示す様に、速度脈波の特徴量のe1(振幅)に相当する値である。なお、このe1とは、加速度脈波e−f点の間にある速度脈波のピーク値である。また、後隆波にピークが見られない場合には、ノッチから所定時間経過後までの傾きを採用してよい。
【0022】
尚、心拍数の代わりに心拍間隔RRIを用いてよい。
(7)請求項7の発明では、前記心電信号及び前記脈波信号以外の生体情報も加味して、血圧を算出することを特徴とする。
【0023】
これにより、一層血圧の算出精度が向上する。
(8)請求項8の発明では、前記心電信号及び前記脈波信号以外の生体情報は、身体特徴量及び生体のインピーダンスのうち少なくとも1種であることを特徴とする。
【0024】
本発明は、脈波信号以外の生体情報を例示したものである。例えば身長、体重、年齢、性別の様な身体特徴情報は、血流量(従って血圧)に相関関係があるので、この身体特徴情報を用いることにより、精度良く血圧を算出することができる。また、インピーダンスは、体重や体脂肪(従って血流量)に相関関係があるので、精度の高い血圧測定が可能である。
【0025】
(9)請求項9の発明は、生体から検出された心電信号及び脈波信号のうち、少なくとも前記脈波信号から得られる生体情報を用いて、血圧を算出する血圧測定装置であって、脈波波形の立ち下がりの状態を示す立ち下がり情報を用いて、血圧を算出することを特徴とする。
【0026】
脈波波形の立ち下がり情報は、血管弾性の情報を含むものである。つまり、脈波波形の立ち下がりの部分は、脈波の立ち上がりの部分と異なり、心拍出波形の影響を受けづらいため、血管弾性情報を心臓の影響を少なくして抽出できる。この血管弾性情報を用いることにより、血圧を精度良く算出できる。
【0027】
ここで、立ち下がりの範囲とは、大動脈弁閉鎖によるノッチ直後の脈波(後隆波)のピークからの立ち下がりを示しており、ノッチ直後にピークが無い場合には、例えば4階微分のゼロクロス点に対応する脈波の位置からの立ち下がりを採用できる。
【0028】
なお、立ち下がり情報としては、複数の情報を用いると、血圧の算出精度が向上するので、好適である。
(10)請求項10の発明では、前記立ち下がり情報として、脈波波形の立ち下がりの最大傾きを用いることを特徴とする。
【0029】
本発明は、好ましい立ち下がり情報を例示したものである。本発明者等の研究によれば、脈波波形の立ち下がりの最大傾きを用いることにより、血圧を精度良く算出できることが確認されている。
【0030】
(11)請求項11の発明では、前記立ち下がり情報として、加速度脈波の特徴量のe点以降における脈波波形の最大値の8〜2割の範囲の傾きを用いることを特徴とする。
本発明は、好ましい立ち下がり情報を例示したものである。本発明者等の研究によれば、図3に示す様に、加速度脈波の特徴量のe点以降における脈波波形の最大値の8〜2割の範囲の傾きを用いることにより、血圧を精度良く算出できることが確認されている。
【0031】
(12)請求項12の発明では、(前記請求項9〜11のいずれかに記載の血圧測定装置に加えて)更に、前記請求項1〜8のいずれかに記載の血圧測定装置の構成を備えたことを特徴とする。
【0032】
これにより、血圧の算出精度を一層高めることができる。
(13)請求項13の発明は、前記請求項1〜12のいずれかに記載の血圧測定装置の機能を実現させるためのプログラムである。
【0033】
よって、マイクロコンピュータ等でこのプログラムを実施することにより、上述した生体情報などを用い、所定の演算式を使用して、血圧を算出することが可能である。
(14)請求項14の発明は、前記請求項13に記載のプログラムを記録したコンピュータ読取可能な記録媒体である。
【0034】
よって、この記録媒体に記録したプログラムを用いて、上述したように、血圧を算出することができる。
【発明を実施するための最良の形態】
【0035】
次に、本発明の最良の形態(実施形態)について説明する。
[第1実施形態]
本実施形態では、自動車に搭載されてドライバの血圧を測定する血圧測定装置について説明する。
【0036】
a)まず、本実施形態の血圧測定装置のシステム構成を説明する。
図4に示す様に、本実施形態では、自動車に搭載された血圧測定装置1と、ドライバ等に情報の報知を行う報知装置3と、手動にてデータの入力を行うマニュアル入力部5と、ステアリング7に取り付けられた脈波センサ9と、ステアリング7に取り付けられた一対の電極11、13とを備えている。
【0037】
前記血圧測定装置1は、周知のマイクロコンピュータを中心とした電子制御装置であり、脈波センサ9や電極11、13等からの信号に基づいて、血圧の算出や報知装置3の制御等を行う。
【0038】
前記報知装置3は、血圧などを表示する液晶等のディスプレイ15や、その内容を音声等で出力するスピーカ17により構成される。
前記マニュアル入力部5は、マニュアルにて、体重、身長、年齢、性別等の個人別の身体特徴情報が入力可能な例えばキーボードやテンキーやリモコンなどの入力装置である。なお、ディスプレイ15の表示画面をタッチパネルとして、データを入力するようにしてもよい。
【0039】
前記脈波センサ9は、周知の発光素子(LED)や受光素子(PD)を備えた光学式のセンサであり、例えばドライバの指先等に光を照射し、その反射波を利用して脈波(容積脈波)を検出することができる。従って、後述する血圧の算出に用いる脈波信号は、この脈波センサ9から得ることができる。
【0040】
前記一対の電極11、13は、ドライバの左右の手にそれぞれ接触するように、ステアリング9の左右に配置されている。この電極11、13は、心電信号を得るために電圧を印加する心電計の電極として用いられる。ここでは、両電極11、13と血圧測定装置1が、心電計の機能を果たす。従って、後述する血圧の算出に用いる心電信号は、この電極11、13を利用して得ることができる。
【0041】
なお、一対の電極11、13間に電流を流しそれによってインピーダンスが得られる。
b)次に、本実施形態の血圧測定装置等の機能をブロック図にて更に詳しく説明する。
図5に示す様に、本実施形態の血圧測定装置1には、インピーダンス計測部21と、心電/インピーダンス計測切替部23と、心電信号取得部25と、脈波信号取得部27と、心電信号解析部31と、心電及び脈波信号解析部33と、脈波信号解析部35と、血圧演算部37とを備えている。
【0042】
このうち、インピーダンス計測部21は、電極11、13間に微少電流(例えば高周波及び低周波の複数の周波数の交流電流)を流してインピーダンスを測定するものである。 心電信号取得部25は、心臓の活動に伴う電気的興奮を、電極11、13間の電位差(心電信号)として計測するものである。
【0043】
心電/インピーダンス計測切替部23は、インピーダンス計測部21と心電信号取得部25との機能を自動的に切り替えるものである。
脈波信号取得部27は脈波センサ9を駆動して脈波信号を取得するものである。
【0044】
心電信号解析部31は、心電信号を解析して、例えば周知のQT信号等を算出するものである。
心電及び脈波信号解析部33は、前記図1に示す様に、心電信号と脈波信号を用いて、心電信号に対する脈波信号の遅れ時間である脈波伝播時間(PTT)を求めるものである。
【0045】
脈波信号解析部35は、前記図2に示す様に、脈波信号を解析して、1階微分(速度脈波)、2階微分(加速度脈波)、3階微分、4階微分を行うとともに、各微分における特徴量(例えば速度脈波のa1〜k1、加速度脈波のa〜f等)の算出、AI(容積AI)の算出、容積AIの補正などを行うものである。
【0046】
血圧演算部37は、PTT、速度脈波や加速度脈波などの特徴量、容積AI等を用い、所定の演算式に基づいて血圧を算出するものである。
c)次に、血圧の算出に用いる演算式について説明する。
【0047】
(i)血圧の算出式
ここでは、中枢血管特徴量(PTT、脈波面積)と末梢血管特徴量(容積AI)とを用いた演算式、例えば下記式(1)〜(5)のいずれかの演算式を用いて、血圧(推定血圧:EBP)を求める。
・PTTを使用する場合
EBP=α・PTT+β・容積AI+γ ・・(1)
EBP=α・PTT+β・容積AI+γ・(脈波特徴量)+・・+ ・・(2)
EBP=α・PTT+β・容積AI+γ・d+δ・e+ε・f
+ζ・HR+・・+ ・・(3)
・PTTを使用しない場合
EBP=α・脈波面積+β・容積AI+γ ・・(4)
・より精度の良い推定式
EBP=α・PTT+β・容積AI+γ・(脈波特徴量)+δ・(身体特徴量)
+ε・(インピーダンス)+ζ ・・(5)
前記各式で、α、β、γ、δ、ε、ζは係数である。脈波特徴量は、例えばa1、b1など(図2参照)である。HRは脈拍であり、脈拍に代えて脈拍間隔(RRI)を用いることができる。脈波面積とは、脈波の基線をゼロとした場合の波形に囲まれる部分の面積である。身体特徴量は、身長、体重、性別、年齢などである。
【0048】
以下、各式について説明する。
前記式(1)では、推定血圧の算出に、PTTと容積AIを用いる。このうち、PTTは中枢血管の血管弾性を示すパラメータであり、容積AIは末梢血管の血管弾性を示すパラメータであるので、両パラメータを用いることにより、中枢、末梢系の情報を取り入れられているため、精度良く血圧を推定することができる。
【0049】
前記式(2)では、推定血圧の算出に、PTTと容積AIに加えて、脈波特徴量を用いる。よって、前記式(1)よりも精度良く血圧を推定することができる。ここで、脈波特徴量を加味すると推定精度が向上するのは、特徴量の非線形性の効果、及び心臓の拍出の影響等を取り入れられているからである。
【0050】
前記式(3)では、推定血圧の算出に、PTTと容積AIに加えて、加速度脈波のd、e、f、及びHRを用いる。よって、前記式(1)よりも精度良く血圧を推定することができる。ここで、加速度脈波のd、e、fとHRを加味すると推定精度が向上するのは、後述の容積AIを圧AIに補正する効果、及び心臓の拍出の影響の取り入れがあるからである。
【0051】
前記式(4)では、推定血圧の算出に、脈波面積(S)と容積AIを用いている。このうち、Sは中枢血管の血管弾性を示すパラメータであると考えられ(実験により推定)、容積AIは末梢血管の血管弾性を示すパラメータであるので、両パラメータを用いることにより、マニュアルにて体重を入力しなくても、精度良く血圧を推定することができる。
【0052】
前記式(5)では、推定血圧の算出に、PTTと容積AIに加えて、脈波特徴量、身体特徴量、インピーダンスを用いている。よって、前記式(2)よりも精度良く血圧を推定することができる。なお、身長や体重などの身体特徴量は、血流量と相関関係があり、また、インピーダンスは体重や体脂肪と相関関係がある(よって血流量とも相関関係がある)ので、これらを加味することにより、血圧の推定精度が向上すると考えられる。
【0053】
(ii)容積AIの補正式
容積AIはその測定方法の特性上、圧脈波により求めた圧AIよりは精度が低い傾向にある。従って、ここでは、容積AIを下記式(6)又は(7)を用いて補正する。
補正容積AI=容積AI+α・HR+β・ノッチ−後隆波傾き ・・・(6)
補正容積AI=容積AI+α・d+β・e+γ・f ・・・(7)
前記各式で、α、β、γは係数である。また、HRは心拍数である(心拍間隔RRIを用いてもよい)。ノッチ−後隆波傾きとは、上述したノッチ(notch)からノッチ直後の後隆波のピークまでの傾き(即ち加速度脈波の特徴量のe−f点間における脈波の傾き=速度脈波の特徴量のg1)である。d、e、fは、上述した加速度脈波の特徴量である。
【0054】
以下、各式について説明する。
前記式(6)では、容積AIの補正に、HRとノッチ−後隆波傾きを用いているので、圧AIに近い容積AIとすることができる。なお、HRとノッチ−後隆波傾きを用いると良いことは、実験的に検証した。
【0055】
前記式(7)では、容積AIの補正に、加速度脈波の特徴量のd、e、fを用いているので、圧AIに近い容積AIとすることができる。なお、d、e、fを用いると良いことは、実験的に検証した。
【0056】
本実施形態では、前記式(6)、(7)より補正した容積AIを、前記式(1)〜(5)に用いることにより、血圧の推定精度が一層向上する。
d)次に、血圧検出装置1にて実施される制御処理について説明する。
【0057】
図6のフローチャートに示す様に、ステップ(S)100では、身体特徴情報を利用するか否かを判定する。例えばディスプレイ15に、身体特徴情報を利用するかどうかを表示し、ドライバの入力を促すようにしても良い。ここで肯定判断されるとステップ110に進み、一方否定判断されるとステップ120に進む。
【0058】
ステップ110では、身体特徴情報の入力が選択されたので、ドライバ等により、例えばタッチパネルによって体重等の入力がされた場合には、そのデータを取得する。
ステップ120では、インピーダンスの計測を行うか否かを判定する。例えばディスプレイ15に、インピーダンスを計測するか否かを表示し、ドライバの入力を促すようにしても良い。ここで肯定判断されるとステップ130に進み、一方否定判断されるとステップ140に進む。
【0059】
ステップ130では、例えばタッチパネルによって体重等の入力がされた場合には、そのデータを取得する。
続くステップ135では、インピーダンスの計測を行う。
【0060】
ステップ140では、心電及び脈波の同時計測を行う。なお、脈波しか使用しない場合は、脈波のみの計測を行う。
続くステップ150では、測定によって得られた心電信号及び脈波信号の解析を行って、血圧の算出や容積AIの補正に用いる演算値(特徴量など)を算出する。
【0061】
例えば、前記式(1)を用いて血圧を算出する場合には、 脈波伝播時間(PTT)と容積AIが必要であるので、ここでは、心電信号と脈波信号とを用いて脈波伝播時間(PTT)を算出するとともに、P1、P2を用いて容積AIを算出する。なお、他の式を用いる場合には、各式に応じて、血圧の算出に必要な特徴量などの値を算出する。
【0062】
続くステップ160では、心電信号及び脈波信号の妥当性を確認する。つまり、心電信号や脈波信号が測定値として妥当な値であるか否かを判定する。ここで、妥当な値であると判断されるとステップ170に進み、妥当な値でないと判断されると、再度心電及び脈波の測定を行うために、ステップ140に戻る。
【0063】
続くステップ170では、容積AIの補正が必要な場合には、ここで容積AIの補正を行う。
例えば前記式(6)を用いて容積AIを補正する場合には、心拍数(HR)とノッチ−後隆波傾きとを用いて、補正容積AIを算出する。
【0064】
ステップ180では、例えば前記式(1)を利用し、上述した各ステップの処理にて得られた脈波伝播時間(PTT)と補正容積AIを用いて、血圧を算出する。
続くステップ190では、算出した血圧を、ディスプレイ15に表示したり、スピーカ17で報知し、一旦本処理を終了する。
【0065】
d)次に、本実施形態による効果を確認した実験例について説明する。
図7(a)は、横軸にトノメトリ血圧により得られたトウ骨動脈圧脈波を解析して求めた圧AIをとり、縦軸に光学式脈波センサにより求めた脈波信号から算出した容積AI(未補正)をとったものである。
【0066】
この図7(a)から、圧AIと(補正前の)容積AIとの相関係数(r)は0.82であり、圧AIと容積AIとの相関関係が小さいことが分かる。よって、容積AIを用いて血圧を算出する場合には、圧AIを用いた場合より精度が低くなる。
【0067】
図7(b)は、横軸に圧AIをとり、縦軸に前記式(6)により補正した補正容積AI(補正AI)をとったものである。
この図7(b)から、圧AIと補正容積AIとの相関係数は0.93であり、圧AIと補正容積AIとの相関関係が大きいことが分かる。よって、圧AIに代えて補正容積AIを用いることにより、精度の高い血圧を算出できることが分かる。
【0068】
図8(a)は、横軸にカフの締め付けを行うオシロメトリック式血圧計により測定した実血圧をとり、縦軸にPTTのみを用いて算出した推定血圧(即ち前記式(1)でPTT項のみを用いて推定血圧を算出したもの)をとったものである。
【0069】
この図8(a)から、実血圧と推定血圧との相関が少なく、PTTのみを用いて血圧を推定した場合には、精度の高い血圧を求めることができないことが分かる。
図8(b)は、前記式(1)を用いて推定血圧を求め、実血圧との相関をとったものである。この場合には、PTTと補正容積AIを用いて血圧を推定するので、推定血圧の精度が高いことがわかる。
【0070】
図8(c)は、前記式(1)と式(7)を用いて推定血圧を求め、実血圧との相関をとったものである。この場合には、PTTと補正容積AIとd、e、fを用いて血圧を推定するので、推定血圧の精度が一層高いことがわかる。
【0071】
図9は、前記式(5)と式(7)を用いて推定血圧を求め、実血圧との相関をとったものである。この場合には、PTTと補正容積AIと脈波特徴量とインピーダンスと身体特徴量を用いて血圧を推定しているので、推定血圧の精度がより一層高いことがわかる。
【0072】
e)この様に、本実施形態では、例えば前記式(1)及び式(6)を利用し、心電信号や脈波信号から得られた例えば脈波伝播時間(PTT)や補正容積AIを用いて血圧を算出するので、簡易な手法で精度の高い血圧測定が可能である。
【0073】
なお、血圧の算出に前記式(1)〜(5)を用いる場合は、補正しない容積AIを用いることも可能であるが、補正容積AIを用いた方が、血圧の算出精度が高まるので好適である。
[第2実施形態]
次に、第2実施形態を説明するが、前記第1実施形態と同様な内容の説明は省略する。
【0074】
本実施形態は、ハード構成は前記第1実施形態と同様であるので、使用する演算式や演算処理の手順について説明する。
a)まず、本実施形態における血圧の算出に用いる演算式について説明する。
【0075】
(i)血圧の算出式
ここでは、脈波(詳しくは後隆波)の立ち下がり情報を用いた演算式、例えば下記式(8)〜(16)のいずれかの演算式を用いて、血圧(推定血圧:EBP)を求める。
・立ち下がりの傾きのみを使用する場合
EBP=α・e1 ・・(8)
EBP=α・e1+β・f1 ・・(9)
EBP=α・(0.8−0.6)の傾き+β・(0.4−0.2) ・・(10)
EBP=α・最大傾き ・・(11)
・脈波及びPTTを使用する場合
EBP=α・e1+β・PTT ・・(12)
EBP=α・e1+β・f1+γ・PTT ・・(13)
EBP=α・e1+β・f1+γ・PTT+δ・(脈波特徴量)+・ ・・(14)
・C1、C2を使用する場合
EBP=α・PTT+β・容積AI+γ・C1+δ・C2+ε ・・(15)
・より精度の良い推定式
EBP=α・e1+β・f1+γ・PTT+δ・(脈波特徴量)
+ε・(身体特徴量)+ζ・(インピーダンス)+η ・・(16)
前記各式で、α、β、γ、δ、ε、ζ、ηは係数である。脈波特徴量は、例えばa1、d、AIである。e1は、図2に示す様に、1階微分の脈波の立ち下がりの傾きの最小値であり、f1は、2階微分のf点後での、4階微分のゼロクロスに対応する点での脈波の立ち下がりの傾斜である。身体特徴量は、身長、体重、性別、年齢などである。C1、C2は、脈波の立下りを指数関数で近似することにより求められる指数であり、指数関数での近似により、精度良く情報の抽出が可能となる(Daniel A Duprez,Determinateion of Radial Artery Pulse Wave Analysis in Asymptomatic Individuals,American Journal of Hypertension 2004;17:647−653参照)。
【0076】
以下、各式について説明する。
前記式(8)では、推定血圧の算出に、脈波の立ち下がり情報であるe1を用いるので、精度良く血圧を推定することができる。ここで、e1を用いると推定精度が高いのは、心臓系の影響の少ない血管弾性情報を用いることができているからである。
【0077】
前記式(9)では、推定血圧の算出に、脈波の立ち下がり情報であるe1とf1を用いるので、e1だけを用いる場合よりも、精度良く血圧を推定することができる。ここで、f1を加味すると測定精度が向上するのは、大動脈及び小動脈の血管情報を併せて考慮できているからである。
【0078】
前記式(10)では、推定血圧の算出に、図3に示す様に、後隆波の傾きの一部を使用する(ここでは、後隆波のピークを1とし、それを20%づつ5等分し、各区分における傾きを求める)ので、精度良く血圧を推定することができる。なお、これらの傾きの平均値を用いてもよい。
【0079】
前記式(11)では、推定血圧の算出に、後隆波の最大傾きを用いるので、精度良く血圧を推定することができる。なお、図3の様に区分した領域での傾きの最大値を、最大傾きとしてもよい。ここで、最大傾きを加味すると測定精度が向上するのは、e1とほぼ同義であり、ノイズの多い環境下では、微分して求めるよりも安定して求められるため、精度が向上するからである。
【0080】
前記式(12)では、推定血圧の算出に、脈波の立ち下がり情報であるe1とPTTとを用いるので、精度良く血圧を推定することができる。
前記式(13)では、推定血圧の算出に、脈波の立ち下がり情報であるe1、f1とPTTとを用いるので、より精度良く血圧を推定することができる。
【0081】
前記式(14)では、推定血圧の算出に、脈波の立ち下がり情報であるe1、f1とPTTと脈波特徴量を用いるので、より一層精度良く血圧を推定することができる。
前記式(15)では、推定血圧の算出に、脈波の立ち下がり情報であるC1、C2とPTTと容積脈AI(補正容積AIでも良い)を用いるので、精度良く血圧を推定することができる。
【0082】
前記式(16)では、推定血圧の算出に、脈波の立ち下がり情報であるe1、f1とPTTと脈波特徴量と身体特徴量とインピーダンスとを用いるので、一層精度良く血圧を推定することができる。
【0083】
また前記式(14)〜(16)において、脈波特徴量として、AIもしくはd等を使用すると、請求項12に示すように、請求項1〜8による血圧推定精度向上も同時に行うことになる。
【0084】
この場合は、PTT及びAI、dにより、中枢及び末梢の弾性、e1及びf1により、同様に中枢及び末梢の弾性による血圧推定精度の向上を図るが、AI、dは心拍出波形の影響も受けるため、e1、f1を併せて用いることで、中枢、末梢、心臓系を考慮した血圧推定が実施可能となり、推定精度が大きく向上する。
【0085】
b)次に、血圧検出装置にて実施される制御処理について説明する。
図10のフローチャートに示す様に、ステップ(S)200では、身体特徴情報を利用するか否かを判定する。ここで肯定判断されるとステップ210に進み、一方否定判断されるとステップ220に進む。
【0086】
ステップ210では、例えばタッチパネルによって体重等の入力がされた場合には、そのデータを取得する。
ステップ220では、インピーダンスの計測を行うか否かを判定する。ここで肯定判断されるとステップ230に進み、一方否定判断されるとステップ240に進む。
【0087】
ステップ230では、例えばタッチパネルによって体重等の入力がされた場合には、そのデータを取得する。
続くステップ235では、インピーダンスの計測を行う。
【0088】
ステップ240では、心電及び脈波の同時計測を行う。なお、脈波しか使用しない場合は、脈波のみの計測を行う。
続くステップ250では、測定によって得られた心電信号及び脈波信号の解析を行って、血圧の算出に必要な特徴量などを求める。
【0089】
続くステップ260では、心電信号及び脈波信号の妥当性を確認する。ここで、妥当な値であると判断されるとステップ270に進み、妥当な値でないと判断されると、ステップ240に戻る。
【0090】
続くステップ270では、前記各式(8)〜(16)を利用して、血圧を算出する。
続くステップ280では、算出した血圧を報知し、一旦本処理を終了する。
c)次に、本実施形態による効果を確認した実験例について説明する。
【0091】
図11(a)は、横軸に実血圧をとり、縦軸にPTTのみを用いて算出した推定血圧(即ち前記式(3)でPTTの項のみを用いて推定血圧を算出したもの)をとったものである。
【0092】
この図11(a)から、PTTのみを用いて血圧を推定した場合には、精度の高い血圧を求めることができないことが分かる。
図11(b)は、前記式(12)を用いて推定血圧を求め、実血圧との相関をとったものである。この場合には、e1とPTTと用いて血圧を推定しているので、推定血圧の精度が高いことがわかる。
【0093】
図11(c)は、前記式(13)を用いて推定血圧を求め、実血圧との相関をとったものである。この場合には、e1、f1とPTTを用いて血圧を推定しているので、推定血圧の精度が一層高いことがわかる。
【0094】
この様に、本実施形態では、脈波の立ち下がり情報を用いて血圧を推定するので、簡易な手法で精度の高い血圧測定が可能である。
尚、本発明は前記実施形態になんら限定されるものではなく、本発明を逸脱しない範囲において種々の態様で実施しうることはいうまでもない。
【0095】
例えば、上述した血圧測定装置の機能は、コンピュータのプログラムにより実行される処理により実現することができ、このプログラムは、記録媒体に記録できる。
つまり、上述した様なプログラムをコンピュータシステムにて実現する機能は、例えば、コンピュータシステム側で起動するプログラムとして備えることができる。このようなプログラムの場合、例えば、フレキシブルディスク、光磁気ディスク、CD−ROM、DVD、ハードディスク等のコンピュータ読み取り可能な記録媒体に記録し、必要に応じてコンピュータシステムにロードして起動することにより用いることができる。この他、ROMやバックアップRAM等をコンピュータ読み取り可能な記録媒体として前記プログラムを記録しておき、このROMあるいはバックアップRAM等をコンピュータシステムに組み込んで用いても良い。
【図面の簡単な説明】
【0096】
【図1】心電信号と脈波信号を示すグラフである。
【図2】脈波信号とそれを微分した信号を示すグラフである。
【図3】後隆波の傾きを示すグラフである。
【図4】第1実施形態の車室内に配置された血圧測定装置の概略構成を示す説明図である。
【図5】血圧測定装置の機能を示すブロック図である。
【図6】第1実施形態の制御処理を示すフローチャートである。
【図7】(a)圧AIと容積AIとの関係を示すグラフ、(b)圧AIと補正容積AIとの関係を示すグラフである。
【図8】カフを用いた実測値と血圧測定装置による演算値との相関関係を示すグラフである。
【図9】カフを用いた実測値と血圧測定装置による演算値との相関関係を示すグラフである。
【図10】第2実施形態の制御処理を示すフローチャートである。
【図11】カフを用いた実測値と血圧測定装置による演算値との相関関係を示すグラフである。
【符号の説明】
【0097】
1…血圧測定装置
3…報知装置
5…マニュアル入力部
7…ステアリング
9…脈波センサ
11、13…電極
【特許請求の範囲】
【請求項1】
生体から検出された心電信号及び脈波信号のうち、少なくとも前記脈波信号から得られる生体情報を用いて、血圧を算出する血圧測定装置であって、
中枢血管の血管弾性を表す中枢血管特徴量と、末梢血管の血管弾性を表す末梢血管特徴量とを用いて、血圧を算出することを特徴とする血圧測定装置。
【請求項2】
前記脈波信号は、容積脈波を表す信号であることを特徴とする請求項1に記載の血圧測定装置。
【請求項3】
前記中枢血管特徴量は、前記心電信号と前記脈波信号との時間差、又は脈波波形により得られる進行波から反射波までの時間、若しくはその値を用いて求められるものであることを特徴とする請求項1又は2に記載の血圧測定装置。
【請求項4】
前記末梢血管特徴量は、脈波のAI、又は加速度脈波の特徴量のd、若しくはその値を用いて求められるものであることを特徴とする請求項1〜3のいずれかに記載の血圧測定装置。
【請求項5】
前記容積脈波から得られる容積AIを、前記加速度脈波の特徴量のd、e、fのうち1種を用いて補正することを特徴とする請求項4に記載の血圧測定装置。
【請求項6】
前記容積脈波から得られる容積AIを、前記加速度脈波の特徴量のe−f点間における脈波の傾きと、心拍数とを用いて補正することを特徴とする請求項4に記載の血圧測定装置。
【請求項7】
前記心電信号及び前記脈波信号以外の生体情報も加味して、血圧を算出することを特徴とする請求項1〜6のいずれかに記載の血圧測定装置。
【請求項8】
前記心電信号及び前記脈波信号以外の生体情報は、身体特徴量及び生体のインピーダンスのうち少なくとも1種であることを特徴とする請求項7に記載の血圧測定装置。
【請求項9】
生体から検出された心電信号及び脈波信号のうち、少なくとも前記脈波信号から得られる生体情報を用いて、血圧を算出する血圧測定装置であって、
脈波波形の立ち下がりの状態を示す立ち下がり情報を用いて、血圧を算出することを特徴とする血圧測定装置。
【請求項10】
前記立ち下がり情報として、脈波波形の立ち下がりの最大傾きを用いることを特徴とする請求項9に記載の血圧測定装置。
【請求項11】
前記立ち下がり情報として、加速度脈波の特徴量のe点以降における脈波波形の最大値の8〜2割の範囲の傾きを用いることを特徴とする請求項9に記載の血圧測定装置。
【請求項12】
更に、前記請求項1〜8のいずれかに記載の血圧測定装置の構成を備えたことを特徴とする請求項9〜11のいずれかに記載の血圧測定装置。
【請求項13】
前記請求項1〜12のいずれかに記載の血圧測定装置の機能を実現させるためのプログラム。
【請求項14】
前記請求項13に記載のプログラムを記録したコンピュータ読取可能な記録媒体。
【請求項1】
生体から検出された心電信号及び脈波信号のうち、少なくとも前記脈波信号から得られる生体情報を用いて、血圧を算出する血圧測定装置であって、
中枢血管の血管弾性を表す中枢血管特徴量と、末梢血管の血管弾性を表す末梢血管特徴量とを用いて、血圧を算出することを特徴とする血圧測定装置。
【請求項2】
前記脈波信号は、容積脈波を表す信号であることを特徴とする請求項1に記載の血圧測定装置。
【請求項3】
前記中枢血管特徴量は、前記心電信号と前記脈波信号との時間差、又は脈波波形により得られる進行波から反射波までの時間、若しくはその値を用いて求められるものであることを特徴とする請求項1又は2に記載の血圧測定装置。
【請求項4】
前記末梢血管特徴量は、脈波のAI、又は加速度脈波の特徴量のd、若しくはその値を用いて求められるものであることを特徴とする請求項1〜3のいずれかに記載の血圧測定装置。
【請求項5】
前記容積脈波から得られる容積AIを、前記加速度脈波の特徴量のd、e、fのうち1種を用いて補正することを特徴とする請求項4に記載の血圧測定装置。
【請求項6】
前記容積脈波から得られる容積AIを、前記加速度脈波の特徴量のe−f点間における脈波の傾きと、心拍数とを用いて補正することを特徴とする請求項4に記載の血圧測定装置。
【請求項7】
前記心電信号及び前記脈波信号以外の生体情報も加味して、血圧を算出することを特徴とする請求項1〜6のいずれかに記載の血圧測定装置。
【請求項8】
前記心電信号及び前記脈波信号以外の生体情報は、身体特徴量及び生体のインピーダンスのうち少なくとも1種であることを特徴とする請求項7に記載の血圧測定装置。
【請求項9】
生体から検出された心電信号及び脈波信号のうち、少なくとも前記脈波信号から得られる生体情報を用いて、血圧を算出する血圧測定装置であって、
脈波波形の立ち下がりの状態を示す立ち下がり情報を用いて、血圧を算出することを特徴とする血圧測定装置。
【請求項10】
前記立ち下がり情報として、脈波波形の立ち下がりの最大傾きを用いることを特徴とする請求項9に記載の血圧測定装置。
【請求項11】
前記立ち下がり情報として、加速度脈波の特徴量のe点以降における脈波波形の最大値の8〜2割の範囲の傾きを用いることを特徴とする請求項9に記載の血圧測定装置。
【請求項12】
更に、前記請求項1〜8のいずれかに記載の血圧測定装置の構成を備えたことを特徴とする請求項9〜11のいずれかに記載の血圧測定装置。
【請求項13】
前記請求項1〜12のいずれかに記載の血圧測定装置の機能を実現させるためのプログラム。
【請求項14】
前記請求項13に記載のプログラムを記録したコンピュータ読取可能な記録媒体。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【公開番号】特開2008−302127(P2008−302127A)
【公開日】平成20年12月18日(2008.12.18)
【国際特許分類】
【出願番号】特願2007−154086(P2007−154086)
【出願日】平成19年6月11日(2007.6.11)
【出願人】(000004260)株式会社デンソー (27,639)
【出願人】(504147243)国立大学法人 岡山大学 (444)
【Fターム(参考)】
【公開日】平成20年12月18日(2008.12.18)
【国際特許分類】
【出願日】平成19年6月11日(2007.6.11)
【出願人】(000004260)株式会社デンソー (27,639)
【出願人】(504147243)国立大学法人 岡山大学 (444)
【Fターム(参考)】
[ Back to top ]