(57)【特許請求の範囲】
【請求項１】
心臓から発生する心音を検出して該心音を表す心音信号を出力する心音マイクと、
該心音マイクから出力される心音信号のうち、第１音を表す心音信号を時間および周波数について解析する時間周波数解析手段と、
該時間周波数解析手段により解析された第１音の時間周波数解析波形について、第１音の初期期間における主成分の周波数帯域よりも高周波側に設定された第１周波数の信号強度が予め定められた閾値を越えた時間に基づいて、大動脈弁の開放時期を決定する大動脈弁開放時期決定手段と
を、含むことを特徴とする心音検出装置。
【請求項２】
心臓から発生する心音を検出して該心音を表す心音信号を出力する心音マイクと、
該心音マイクから出力される心音信号のうち、第１音を表す心音信号を時間および周波数について解析する時間周波数解析する時間周波数解析手段と、
該時間周波数解析手段により解析された第１音の時間周波数解析波形について、第１音の初期期間における主成分の周波数帯域内の信号強度が極小となった時間に基づいて、大動脈弁の開放時期を決定する大動脈弁開放時期決定手段と
を、含むことを特徴とする心音検出装置。
【請求項３】
生体の心臓の収縮開始から血液が拍出されるまでの前駆出期間を測定する前駆出期間測定装置であって、
前記生体の心電誘導波形を検出する心電誘導装置と、
前記請求項１または２記載の心音検出装置と、
前記心電誘導装置により心電誘導波形の心室筋の興奮を表す波形が検出された時点から、前記心音検出装置の大動脈弁開放時期決定手段により決定された大動脈弁の開放時期までの時間差を前駆出期間として算出する前駆出期間算出手段と
を含むことを特徴とする前駆出期間測定装置。
【請求項４】
生体の動脈内を伝播する脈波の伝播速度に関連する脈波伝播速度情報を測定するための脈波伝播速度情報測定装置であって、
前記請求項１または２記載の心音検出装置と、
前記生体の一部に装着され、該生体の動脈を通して伝播する脈波を検出する脈波検出装置と、
前記心音検出装置の大動脈弁開放時期決定手段により決定された大動脈弁の開放時期と前記脈波検出装置により脈波の立ち上がり点が検出された時点との時間差に基づいて、前記脈波伝播速度情報を算出する脈波伝播速度情報算出手段と
を含むことを特徴とする脈波伝播速度情報測定装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、大動脈弁の開放時期を決定することができる心音検出装置、および、その心音検出装置を利用した前駆出期間測定装置・脈波伝播速度情報測定装置に関するものである。
【０００２】
【従来の技術】
心音の第１音と左心系（すなわち、左心室、左心房、大動脈）の内圧との関係は以下のようになっている。まず、心室筋が収縮を始めると、左心室内圧は急速に上昇し、左心房内圧を上回ったときに僧帽弁が閉鎖し、その僧帽弁の閉鎖によって発生する音が第１音の初期部分を構成する。左心室内圧がさらに急上昇を続け、大動脈圧を越えると、大動脈弁が開き、左心室から大動脈への血液の流出（駆出）が始まる。この大動脈弁の開放によって生じる音が第１音の後半部分を構成する。
【０００３】
【発明が解決しようとする課題】
大動脈弁の開放によって発生する音は、僧帽弁の閉鎖によって発生する音が消失しないうちに発生することから、第１音には、僧帽弁の閉鎖によって発生する音に、大動脈弁の開放によって発生する音が重畳している。さらに、第１音には、右心系の音（たとえば、僧帽弁とほぼ同時に閉鎖する三尖弁の閉鎖によって発生する音など）が含まれ、また、生体内において発生する体内ノイズ等の影響もあることから、第１音は複雑な形状であり、第１音から大動脈弁の開放時期を決定することは困難であった。
【０００４】
本発明は以上のような事情を背景として為されたものであり、その目的とするところは、大動脈弁の開放時期を正確に決定することができる心音検出装置を提供することにある。
【０００５】
【課題を解決するための第１の手段】
上記目的を達成するための第１発明の要旨とするところは、心音検出装置において、(a) 心臓から発生する心音を検出してその心音を表す心音信号を出力する心音マイクと、(b) その心音マイクから出力される心音信号のうち、第１音を表す心音信号を時間および周波数について解析する時間周波数解析する時間周波数解析手段と、(c-1) その時間周波数解析手段により解析された第１音の時間周波数解析波形について、第１音の初期期間における主成分の周波数帯域よりも高周波側に設定された第１周波数の信号強度が予め定められた閾値を越えた時間に基づいて、大動脈弁の開放時期を決定する大動脈弁開放時期決定手段とを、含むことにある。
【０００６】
【第１発明の効果】
大動脈弁の開放時に発生する音は、第１音の初期期間において主成分を構成する僧帽弁および三尖弁の閉鎖音よりも遅れて発生し、第１音の後半部分の一部を構成することから、このように、大動脈弁開放時期決定手段により、時間周波数解析手段によって得られた時間周波数解析波形について、第１音の初期期間における主成分の周波数帯域よりも高周波側に設定された第１周波数が予め設定された閾値を越えたことに基づいて大動脈弁の開放時期が決定されれば、大動脈弁の開放時期を正確に決定することができる。
【０００７】
【課題を解決するための第２の手段】
上記目的を達成するための第２発明の要旨とするところは、心音検出装置において、(a) 心臓から発生する心音を検出してその心音を表す心音信号を出力する心音マイクと、(b) その心音マイクから出力される心音信号のうち、第１音を表す心音信号を時間および周波数について解析する時間周波数解析する時間周波数解析手段と、(c-2) その時間周波数解析手段により解析された第１音の時間周波数解析波形について、第１音の初期期間における主成分の周波数帯域内の信号強度が極小となった時間に基づいて、大動脈弁の開放時期を決定する大動脈弁開放時期決定手段とを、含むことにある。
【０００８】
【第２発明の効果】
第１音の初期期間における主成分の周波数域内の信号強度は一旦減少していくが、大動脈弁の開放によって再び強度が増加するので、このようにすれば、大動脈弁の開放時期を正確に決定することができる。
【０００９】
【発明の他の態様】
ここで、前記心音検出装置は、前駆出期間測定装置に利用できる。すなわち、その前駆出期間測定装置は、生体の心臓の収縮開始から血液が拍出されるまでの前駆出期間を測定する前駆出期間測定装置であって、前記生体の心電誘導波形を検出する心電誘導装置と、前記心音検出装置と、前記心電誘導装置により心電誘導波形の心室筋の興奮を表す波形が検出された時点から前記心音検出装置の大動脈弁開放手段により決定された大動脈弁の開放時期までの時間差を前駆出期間として算出する前駆出期間算出手段とを含むことを特徴とするものである。このようにすれば、前記心音検出装置の大動脈弁開放時期決定手段により大動脈弁の開放時期が正確に決定され、前駆出期間算出手段により、心電誘導波形の心室筋の興奮を表す波形が検出された時点から前記心音検出装置により決定された大動脈弁の開放時期までの時間差が前駆出期間として算出されるので、正確な前駆出期間が算出できる。
【００１０】
また、前記心音検出装置は、脈波伝播時間または脈波伝播速度などの脈波伝播速度情報を測定する脈波伝播速度情報測定装置にも利用できる。すなわち、その脈波伝播速度情報測定装置は、生体の動脈内を伝播する脈波の伝播速度に関連する脈波伝播速度情報を測定するための脈波伝播速度情報測定装置であって、前記心音検出装置と、前記生体の一部に装着され、その生体の動脈を通して伝播する脈波を検出する脈波検出装置と、前記心音検出装置の大動脈弁開放時期決定手段により決定された大動脈弁の開放時期と前記脈波検出装置により脈波の立ち上がり点が検出された時点との時間差に基づいて、前記脈波伝播速度情報を算出する脈波伝播速度情報算出手段とを含むことを特徴とするものでる。このようにすれば、前記心音検出装置の大動脈弁開放時期決定手段によって大動脈弁の開放時期が正確に決定され、脈波伝播速度情報算出手段では、その正確に決定された大動脈弁の開放時期と、脈波検出装置により脈波の立ち上がり点が検出された時点との時間差に基づいて脈波伝播速度情報が算出される。上記大動脈弁の開放時期と脈波の立ち上がり点が検出された時点との時間差には前駆出期間が含まれないことから、その時間差に基づいて正確な脈波伝播速度情報が算出できる。
【００１１】
因みに、従来、脈波伝播速度情報を測定する場合、上流側において心拍同期波を検出するために、心電誘導波形を検出する心電誘導装置や心音マイクが用いられることがあり、その場合、心電誘導波形の心室筋の興奮を表す信号（すなわちＱ波、Ｒ波、またはＳ波）が検出された時点、或いは心音の第１音の開始点が第１基準点として用いられることがある。しかし、その第１基準点と、下流側において脈波検出装置により検出される第２基準点との時間差には、心室筋が興奮を開始してから実際に大動脈へ血液が駆出されるまでの前駆出期間ＰＥＰが含まれる。そのため、脈波伝播時間ＤＴを算出するためには、上記時間差から前駆出期間ＰＥＰを差し引く必要があるが、前駆出期間ＰＥＰを測定することは容易でなかったことから、前駆出期間ＰＥＰは予め決定された一定値とされることが多かった。しかしながら、実際には、前駆出期間ＰＥＰは心筋の状態によって変動するので、上記従来の方法によって算出される脈波伝播速度情報は精度が不十分だったのである。なお、前駆出期間ＰＥＰは、心筋の前負荷或いは後負荷等によって変動することから、心機能を評価する指標となる。
【００１２】
【発明の好適な実施の形態】
以下、本発明の一実施例を図面に基づいて詳細に説明する。図１は、本発明が適用され、心音検出装置および前駆出期間測定装置としての機能を備えた脈波伝播速度情報測定装置１０の構成を説明するブロック図である。
【００１３】
図１において、マイクロホン１２は、心音マイクとして機能するものであり、被測定者１４の胸部上に図示しない粘着テープ等により固着される。そして、マイクロホン１２の図示しない内部に備えられている圧電素子において、被測定者１４の心臓から発生する心音等が電気信号すなわち心音信号ＳＨに変換される。心音信号増幅器１６には、心音の高音成分をよく記録するためにエネルギーの大きい低音成分を弱める図示しない４種類のフィルタが備えられ、マイクロホン１２から出力される心音信号ＳＨが増幅され、且つ、ろ波された後に、Ａ／Ｄ変換器１８を介して電子制御装置２０へ供給される。
【００１４】
１組の電極２２は、心筋の活動電位を示す心電信号ＳＥを出力するために、生体の所定部位に装着される。本実施例においては、II誘導による心電図を測定するために、被測定者１４の右手首と左足首とにそれぞれ装着されている。電極２２から出力された心電信号ＳＥは、所定の増幅器を備えた心電誘導装置２４により増幅された後、Ａ／Ｄ変換器２６を介し前記電子制御装置２０へ供給される。
【００１５】
光電脈波センサ２８は、毛細血管を含む末梢細動脈へ伝播した脈波を検出する脈波検出装置として機能し、たとえば脈拍検出などに用いるものと同様に構成されており、生体の一部（たとえば指尖部）を収容可能な図示しないハウジング内に、図示しない発光素子および受光素子を備えている。その発光素子は、ヘモグロビンによって反射可能な波長帯の赤色光或いは赤外光、好ましくは酸素飽和度によって影響を受けない８００ｎｍ程度の波長を生体の表皮に向かって照射する光源であり、受光素子は、表皮内からの散乱光を検出し、毛細血管内の血液容積に対応する光電脈波信号ＳＭを出力する。光電脈波センサ２８から出力された光電脈波信号ＳＭは、Ａ／Ｄ変換器３１を介して電子制御装置２０へ供給される。この光電脈波信号ＳＭは、一拍毎に脈動する信号であって、表皮内の毛細血管内のヘモグロビンの量すなわち血液量に対応している。
【００１６】
さらに、電子制御装置２０には、押しボタン２９から起動信号ＳＳが供給されるようになっている。
【００１７】
上記電子制御装置２０は、ＣＰＵ３０，ＲＯＭ３２，ＲＡＭ３４，および図示しないＩ／Ｏポート等を備えた所謂マイクロコンピュータにて構成されており、ＣＰＵ３０は、ＲＯＭ３２に予め記憶されたプログラムに従ってＲＡＭ３４の記憶機能を利用しつつ信号処理を実行することにより、心音信号ＳＨが表す心音波形を時間周波数解析して大動脈弁の開放時期を決定し、さらに、その大動脈弁の開放時期に加えて、心電信号ＳＥおよび光電脈波信号ＳＭに基づいて、前駆出期間ＰＥＰおよび脈波伝播速度情報を算出し、その算出した前駆出期間ＰＥＰおよび脈波伝播速度情報を表示器３６の表示画面に表示させる。
【００１８】
図２は、上記脈波伝播速度情報測定装置１０における電子制御装置２０の制御機能の要部を説明する機能ブロック線図である。時間周波数解析手段４０は、マイクロホン１２から供給される心音信号ＳＨが表す心音のうち、第１音を含む部分を切り出し、その切り出した心音信号ＳＨを時間周波数解析する。時間周波数解析とは、信号を時間と周波数の両面から同時に解析して周波数解析値の時間変化を示すものであり、たとえば、ウェーブレット変換や、所定の時間区分毎の高速フーリエ変換（ＦＦＴ）等の手法が用いられる。
【００１９】
上記ウェーブレット変換とは、図３にその一例が示されるウェーブレット関数ψ（ｔ）を、時間軸方向に平行移動させる移動変数ｂと、ウェーブレット関数ψ（ｔ）が表す波形の時間軸方向の大きさを伸縮させる伸縮変数ａとの関数として、そのウェーブレット関数ψ（（ｔ−ｂ）／ａ）と切り出した心音信号ＳＨを表す関数ｆ（ｔ）との積を時間ｔについて積分して得られる上記ａおよびｂの関数として定義される。すなわち、下記式１のように定義される。なお、上記ウェーブレット関数ψ（（ｔ−ｂ）／ａ）において、伸縮変数ａに対応してψ（ｔ）の幅がａ倍になることから、１／ａが周波数に対応し、移動変数ｂに対応してψ（ｔ）が時間軸方向に平行移動することから、ｂは時間に対応する。
【００２０】
【式１】
0004352558.tif 000002
【００２１】
図４、図５は、上記式１のウェーブレット変換式の意味を説明するための図であり、図４（Ａ）は、パラメータａ、ｂを適当に選ぶことにより、ウェーブレット関数ψ（（ｔ−ｂ）／ａ）がある関数ｇ（ｔ）の一部分に略一致している状態を示し、図５（Ａ）は、ウェーブレット関数ψ（（ｔ−ｂ）／ａ）がある関数ｈ（ｔ）の一部分を近似していない状態を示している。そして、図４（Ｂ）および図５（Ｂ）は、図４（Ａ）および図５（Ａ）の場合におけるウェーブレット関数ψ（（ｔ−ｂ）／ａ）と関数ｇ（ｔ）またはｈ（ｔ）との積をそれぞれ示す図である。図４（Ｂ）に示されるように、ウェーブレット関数ψ（（ｔ−ｂ）／ａ）が関数ｇ（ｔ）の一部分に略一致している場合は、ψ（（ｔ−ｂ）／ａ）とｇ（ｔ）との積は符号の変化がないので、積分値は大きくなる。しかし、図５（Ｂ）に示されるように、ウェーブレット関数ψ（（ｔ−ｂ）／ａ）が関数ｈ（ｔ）の一部分を近似していない場合は、ψ（（ｔ−ｂ）／ａ）とｈ（ｔ）との積はｔの変化とともに激しく符号が変化するので、積分値は小さくなる。従って、上記式１は、パラメータａ、ｂを変更することにより、ウェーブレット関数ψ（（ｔ−ｂ）／ａ）が、切り出した心音信号ＳＨを表す関数ｆ（ｔ）の一部分に似ている場合に大きい値を示し、切り出した心音信号ＳＨを表す関数ｆ（ｔ）の一部分に似ていない場合は小さい値を示す。
【００２２】
上記式１において、周波数に対応する伸縮変数ａおよび時間に対応する移動変数ｂを逐次変化させていき、その都度、ウェーブレット変換を行なうと、時間、周波数、および信号強度の三次元の時間周波数解析波形が得られる。図６（ａ）は、切り出した心音信号ＳＨの一例を示し、図６（ｂ）は、図６（ａ）の心音信号ＳＨをウェーブレット変換して得た時間周波数解析波形を、等高線図（等強度線図）として表したグラフである。
【００２３】
また、高速フーリエ変換は、所定の時間区分内に得られた信号（この信号は、時間と信号強度の二次元の信号である）を周波数および信号強度の二次元の信号に変換するものであるので、上記時間区分を移動させ、その都度、高速フーリエ変換を行なうと、前述のウェーブレット変換と同様に、時間、周波数、および信号強度の三次元の時間周波数解析波形が得られる。
【００２４】
大動脈弁開放時期決定手段４２は、上記時間周波数解析手段４０により解析された第１音の周波数解析値の時間変化すなわち第１音の時間周波数解析波形に基づいて、たとえば第１音の初期期間における主成分の周波数帯域ＲＦ（以下、初期主周波数帯域という）よりも高周波側に設定された第１周波数ｆ1 の信号強度が予め設定された閾値ＴＨを越えた時間ｔに基づいて大動脈弁の開放時期を決定する。
心音の第１音は、僧帽弁および三尖弁の閉鎖によって開始することから、第１音は、当初、僧帽弁および三尖弁の閉鎖音により構成され、大動脈弁の開放音は含まれていない。僧帽弁および三尖弁の閉鎖後、等容収縮期を経過すると大動脈弁が開放するが、第１音の初期期間とは、第１音に大動脈弁の開放音が含まれるようになる前の期間、すなわち、上記等容収縮期よりも短い期間として予め定められた期間を意味する。また、第１音の初期期間の主成分は、僧帽弁および三尖弁の閉鎖音の主成分であり、上記初期主周波数帯域ＲＦは、たとえば、初期期間において、信号強度が予め設定された所定値を越えた周波数の最低値ＲＦmin から最高値ＲＦmax までとされる。
そして、大動脈弁の開放により、初期期間には存在しなかった高周波成分が第１音に現れることが実験により判明している。従って、上記第１周波数ｆ1 は、たとえば、第１音の初期主周波数帯域ＲＦの最高周波数ＲＦmax に所定値を加えた周波数、或いは、その最高周波数ＲＦmax に１以上の所定値を掛けた周波数に決定される。または、大動脈弁の開放に特徴的な周波数として予め実験に基づいて決定されてもよい。
また、上記閾値ＴＨは、たとえば、第１音に含まれる体内ノイズや外来ノイズの信号レベルよりも大きい範囲で可及的に低い値に設定される。
図６（ｂ）の大動脈弁開放時期ｔO は、この大動脈弁開放時期決定手段４２により決定された大動脈弁開放時期ｔO の一例であり、第１周波数ｆ1 が予め実験に基づいて決定された場合である。
【００２５】
前駆出期間算出手段４４は、心電誘導装置２４により心電誘導波形の心室筋の興奮を表す波形（Ｑ波、Ｒ波、Ｓ波）が検出された時点から、前記大動脈弁開放時期決定手段４２により決定された大動脈弁の開放時期までの時間差を前駆出期間ＰＥＰとして算出し、その算出した前駆出期間ＰＥＰを表示器３６に表示する。図７には、心室筋の興奮を表す波形としてＲ波が用いられた場合の前駆出期間ＰＥＰが示してある。なお、図７の大動脈起始部脈波形は説明の便宜上示してあるが、本実施例においては測定されない。
【００２６】
脈波伝播速度情報算出手段４６は、図７にも図示するように、大動脈弁開放時期決定手段４２により決定された大動脈弁の開放時期ｔO と前記光電脈波センサ２８により検出された光電脈波の立ち上がり点までの時間差すなわち脈波伝播時間ＤＴ(sec) を算出する脈波伝播時間算出手段を含み、予め設定された式２から、その脈波伝播時間算出手段により算出された脈波伝播時間ＤＴに基づいて脈波伝播速度ＰＷＶ(m/sec) を算出し、その算出した脈波伝播速度ＰＷＶを表示器３６に表示する。なお、式２において、Ｌは、大動脈起始部から光電脈波センサ２８の装着位置までの伝播距離(m) であり、予め実験に基づいて決定されたものである。上記大動脈弁の開放時期ｔO は、実際に大動脈へ血液が駆出される時点を示すものであるので、前駆出期間ＰＥＰが含まれていない正確な脈波伝播時間ＤＴが得られ、また、その脈波伝播時間ＤＴに基づいて算出される脈波伝播速度ＰＷＶも正確なものとなる。
【００２７】
（式２）    ＰＷＶ＝Ｌ／ＤＴ
【００２８】
図８は、図２の機能ブロック線図に示した電子制御装置２０の制御作動をさらに具体的に説明するためのフローチャートである。本ルーチンは、押しボタン２９が操作され、起動信号ＳＳが供給された場合に実行される。
【００２９】
図８において、ステップＳ１（以下、ステップを省略する。）では、マイクロホン１２から供給される心音信号ＳＨ、心電誘導装置２４から供給される心電信号ＳＥ、および光電脈波センサ２８から供給される光電脈波信号ＳＭが読み込まれる。
【００３０】
続くＳ２では、上記心音信号ＳＨ、心電信号ＳＥ、および光電脈波信号ＳＭが一拍分読み込まれたか否かが、たとえば光電脈波の立ち上がり点に基づいて判断される。この判断が否定された場合には、前記Ｓ１が繰り返し実行されることにより、心音信号ＳＨ、心電信号ＳＥ、および光電脈波信号ＳＭの読み込みが継続される。
【００３１】
一方、上記Ｓ２の判断が肯定された場合には、時間周波数解析手段４０に対応するＳ３乃至Ｓ４が実行される。まず、Ｓ３において、心音信号ＳＨから第１音を十分に含む期間（たとえば１５０msec）が切り出され、続いて、Ｓ４において、上記Ｓ３で切り出された心音信号ＳＨがウェーブレット変換される。
【００３２】
そして、続く大動脈弁開放時期決定手段４２に対応するＳ５では、上記Ｓ４でウェーブレット変換して得られた時間周波数解析波形を用いて、その時間周波数解析波形の、大動脈弁の開放音に特徴的な周波数として予め設定された第１周波数ｆ1 における信号強度が、予め設定された閾値ＴＨを越えた時点が大動脈弁の開放時期ｔO に決定される。
【００３３】
続くＳ６では、前記Ｓ１において読み込まれた心電信号ＳＥから、Ｑ波が発生した時間が決定され、さらに、続く前駆出期間算出手段４４に対応するＳ７では、上記Ｓ６において決定された心電誘導波形のＱ波発生時間から、前記Ｓ５において決定された大動脈弁の開放時期ｔO までの時間差が前駆出期間ＰＥＰとして算出される。
【００３４】
続いて、Ｓ８では、前記Ｓ１において読み込まれた光電脈波信号ＳＭから、光電脈波の立ち上がり点の時間が決定され、続いて、脈波伝播速度情報算出手段４６に対応するＳ９乃至Ｓ１０が実行される。まず、Ｓ９では、前記Ｓ５において決定された大動脈弁の開放時期ｔO から、上記Ｓ８で決定された光電脈波の立ち上がり点までの時間差が脈波伝播時間ＤＴとして算出される。そして、続くＳ１０では、その脈波伝播時間ＤＴが前記式２に代入されることにより、脈波伝播速度ＰＷＶが算出される。
【００３５】
続くＳ１１では、前記Ｓ７で算出された前駆出期間ＰＥＰ、および上記Ｓ１１で算出された脈波伝播速度ＰＷＶが表示器３６に表示される。
【００３６】
上述のように、本実施例によれば、大動脈弁開放時期決定手段４２（Ｓ５）により、時間周波数解析手段４０（Ｓ３乃至Ｓ４）によってウェーブレット変換されて得られた時間周波数解析波形について、初期主周波数帯域ＲＦよりも高周波側に設定された第１周波数ｆ1 が予め設定された閾値ＴＨを越えた時点が大動脈弁の開放時期が決定されるので、大動脈弁の開放時期を正確に決定することができる。
【００３７】
また、本実施例によれば、大動脈弁開放時期決定手段４２（Ｓ５）により大動脈弁の開放時期ｔO が正確に決定され、前駆出期間算出手段４４（Ｓ７）により、心電誘導波形のＱ波が検出された時点から大動脈弁開放時期決定手段４２（Ｓ５）により決定された大動脈弁の開放時期ｔo までの時間差が前駆出期間ＰＥＰとして算出されるので、正確な前駆出期間ＰＥＰが算出できる。
【００３８】
また、本実施例によれば、大動脈弁開放時期決定手段４２（Ｓ５）によって大動脈弁の開放時期ｔo が正確に決定され、脈波伝播速度情報算出手段４６（Ｓ９乃至Ｓ１０）では、その正確に決定された大動脈弁の開放時期ｔo と、光電脈波センサ２８により光電脈波の立ち上がり点が検出された時点との時間差に基づいて、脈波伝播時間ＤＴおよび脈波伝播速度ＰＷＶが算出される。上記大動脈弁の開放時期ｔO と光電脈波の立ち上がり点が検出された時点との時間差には前駆出期間ＰＥＰが含まれないことから、その時間差に基づいて正確な脈波伝播時間ＤＴおよび脈波伝播速度ＰＷＶが算出できる。
【００３９】
以上、本発明の一実施例を図面に基づいて詳細に説明したが、本発明はその他の態様においても適用される。
【００４０】
たとえば、大動脈弁開放時期決定手段は、前記周波数解析手段４０により解析された第１音の時間周波数解析波形について、初期主周波数帯域ＲＦ内の信号強度が極小となった時間に基づいて、大動脈弁の開放時期を決定してもよい。
前述したように、第１音の初期期間の主成分は僧帽弁および三尖弁の閉鎖音であることから、僧帽弁および三尖弁の閉鎖が完了すると、初期主周波数帯域ＲＦの信号強度は低下するが、その後、大動脈弁が開放されると、初期主周波数帯域ＲＦの信号強度も増加するので、初期主周波数帯域ＲＦ内の信号強度が極小となった時間に基づいて大動脈弁の開放時期を決定できるのである。
たとえば、初期主周波数帯域ＲＦの一部或いは全部の周波数帯域の信号強度を平均し、その平均が極小となる時間に基づいて大動脈弁の開放時期ｔo を決定してもよいし、その初期主周波数帯域ＲＦ内の一点（以下、第２周波数ｆ2 という）の信号強度が極小となる時間に基づいて大動脈弁の開放時期ｔO を決定してもよい。図６（ｂ）の大動脈弁開放時期ｔO は、図６（ｂ）に示す周波数帯域ＲＦの中間値を第２周波数ｆ2 として、その第２周波数ｆ2 に基づいて決定したものでもある。すなわち、このようにしても、前述の第１の手段により決定される場合と略同じ時間が大動脈弁の開放時期ｔo に決定される。
また、前述の実施例のように、時間周波数解析手段４０によって時間周波数解析された時間周波数解析波形について、第１周波数ｆ1 が予め設定された閾値ＴＨを越えた時点と、上記のように、時間周波数解析波形の初期主周波数帯域ＲＦ内の信号強度が極小となった時点とをそれぞれ決定し、それらの中間点を大動脈弁の開放時期ｔ0 に決定してもよい。
【００４１】
また、前述の実施例では、脈波検出装置として、生体の指尖部に装着される光電脈波センサ２８が用いられていたが、生体の所定部位を押圧し、その部位に伝播する圧脈波を検出する形式の圧脈波センサ、生体の所定部位（たとえば上腕）に装着される圧迫帯内の圧力の変動を検出する形式の圧脈波センサ、オキシメータ用の光電脈波検出プローブ、指に装着された電極を介してインピーダンス変化を検出するインピーダンス脈波検出装置などが脈波検出装置として用いられてもよい。
【００４２】
また、前述の実施例では、脈波伝播速度情報算出手段４６（Ｓ９）は、大動脈弁の開放時期ｔO から光電脈波の立ち上がり点までの時間差を脈波伝播時間ＤＴとして算出していたが、心電誘導波形の心室筋の興奮を表す波形が検出された時間から光電脈波の立ち上がり点までの時間差を算出し、その時間差から、前駆出期間算出手段４４（Ｓ７）で算出された前駆出期間ＰＥＰを引くことにより脈波伝播時間ＤＴを算出してもよい。
【００４３】
また、前述の実施例において、時間周波数解析手段４０（Ｓ３乃至Ｓ４）では、心音信号ＳＨのうち第１音を含む部分を切り出すことにより、心音信号ＳＨの一部を時間周波数解析していたが、マイクロホン１２から供給される全ての心音信号ＳＨを時間周波数解析してもよい。
【００４４】
また、前述の実施例の脈波伝播速度情報測定装置１０には、前駆出期間ＰＥＰを算出する前駆出期間算出手段４４（Ｓ７）が設けられていたが、前駆出期間算出手段４４（Ｓ７）は設けられなくてもよい。
【００４５】
また、前述の実施例の脈波伝播速度情報測定装置１０において、脈波伝播速度情報算出手段４６（Ｓ９乃至Ｓ１０）が設けられず、前駆出期間ＰＥＰを算出する前駆出期間測定装置として用いられてもよい。
【００４６】
なお、本発明はその主旨を逸脱しない範囲においてその他種々の変更が加えられ得るものである。
【図面の簡単な説明】
【図１】本発明が適用され、心音検出装置および前駆出期間測定装置としての機能を備えた脈波伝播速度情報測定装置の構成を説明するブロック図である。
【図２】図１の脈波伝播速度情報測定装置における電子制御装置の制御機能の要部を説明する機能ブロック線図である。
【図３】ウェーブレット関数の一例を示す図である。
【図４】ウェーブレット関数ψ（（ｔ−ｂ）／ａ）がある関数ｇ（ｔ）の一部分に略一致している状態、およびそのときの２つの関数の積を示す図である。
【図５】ウェーブレット関数ψ（（ｔ−ｂ）／ａ）がある関数ｈ（ｔ）の一部分を近似していない状態、およびそのときの２つの関数の積を示す図である。
【図６】（ａ）は、図２の時間周波数解析手段において切り出した心音信号ＳＨの一例を示す図である。（ｂ）は、（ａ）の心音信号ＳＨをウェーブレット変換して得た時間周波数解析波形を、等高線図として表したグラフである。
【図７】心室筋の興奮を表す波形としてＲ波が用いられた場合の前駆出期間ＰＥＰを例示する図である。
【図８】図２の機能ブロック線図に示した電子制御装置の制御作動をさらに具体的に説明するためのフローチャートである。
【符号の説明】
１０：脈波伝播速度情報測定装置（心音検出装置、前駆出期間測定装置）
１２：マイクロホン（心音マイク）
２４：心電誘導装置
２８：光電脈波センサ（脈波検出装置）
４０：時間周波数解析手段
４２：大動脈弁開放時期決定手段
４４：前駆出期間算出手段
４６：脈波伝播速度情報算出手段
【図１】
0004352558.tif 000003
【図２】
0004352558.tif 000004
【図３】
0004352558.tif 000005
【図４】
0004352558.tif 000006
【図５】
0004352558.tif 000007
【図６】
0004352558.tif 000008
【図７】
0004352558.tif 000009
【図８】
0004352558.tif 000010