(57)【特許請求の範囲】
【請求項１】
血管機能検査装置であって、
検査側の上腕に装着され、駆血及び開放を行う加圧部と、
検査側と対照側の各指尖に装着され、又は検査側と対照側の各手首に装着され、駆血前後の脈波を測定するプローブと、
前記加圧部及び前記プローブを制御する制御部と、
前記制御部に指令を与え、前記制御部から前記脈波に係るデータを取得する処理部と、を備え、
前記処理部が、駆血前後の前記脈波に係るデータに基づいて平均振幅値の計算及びカオス解析を行い、心血管疾患と血管内皮機能障害の判定を行うための特徴量を定量化した数値を算出し、
前記処理部が、前記定量化した数値に基づき、検査側における駆血後の数値を（Ａ）と、同じく駆血前の数値を（Ｂ）と、対照側における駆血後の数値を（Ｃ）と、同じく駆血前の数値を（Ｄ）とするとき、検査側駆血前後比（Ａ／Ｂ）、対照側駆血前後比（Ｃ／Ｄ）及び前記検査側駆血前後比の前記対照側駆血前後比に対する検査側／対照側比（（Ａ／Ｂ）／（Ｃ／Ｄ））を求めることを特徴とする血管機能検査装置。
【請求項２】
前記処理部が、前記カオス解析として、前記脈波の振幅値及びその変化率に基づいて駆血前後の血管ゆらぎを示すカオスアトラクターを算出し、前記算出されたカオスアトラクターについて、リアプノフ指数、外周面積及び短径と長径の固有ベクトルの少なくとも１つを定量化して前記特徴量を定量化した数値を取得することを特徴とする請求項１に記載の血管機能検査装置。
【請求項３】
前記プローブが、光電式又は圧力式によって前記脈波を測定することを特徴とする請求項１又は２に記載の血管機能検査装置。
【請求項４】
前記プローブが光電式であって、測定部位に最適圧力を与えるための空気袋と、前記測定部位に光を照射する発光部と、前記測定部位からの反射光を受光する受光部と、を備えることを特徴とする請求項３に記載の血管機能検査装置。
【請求項５】
前記プローブが、検査側の第１指尖に装着され、細い血管に係る前記脈波を測定する第１プローブを含むことを特徴とする請求項１から４のいずれか１項に記載の血管機能検査装置。
【請求項６】
前記プローブが、検査側の第１手首に装着され、太い血管に係る前記脈波を測定する第２プローブをさらに含むことを特徴とする請求項５に記載の血管機能検査装置。
【請求項７】
前記プローブが、対照側の第２指尖に装着され、細い血管に係る前記脈波を測定する第３プローブと、対照側の第２手首に装着され、太い血管に係る前記脈波を測定する第４プローブと、をさらに含むことを特徴とする請求項６に記載の血管機能検査装置。
【請求項８】
前記処理部が、検査側と対照側の前記指尖の脈波と、検査側と対照側の前記手首の脈波との両方について前記検査側駆血前後比、前記対照側駆血前後比及び前記検査側／対照側比を求めることを特徴とする請求項１から６のいずれか１項に記載の血管機能検査装置。
【請求項９】
前記心血管疾患ごとに既知のカオスアトラクター群を記憶する記憶部をさらに備え、
前記処理部が、前記算出されたカオスアトラクターを前記既知のカオスアトラクター群とパターンマッチングにより対比することを特徴とする請求項２から８のいずれか１項に記載の血管機能検査装置。
【請求項１０】
前記処理部によって処理された結果を表示する表示部をさらに備えることを特徴とする請求項１から９のいずれか１項に記載の血管機能検査装置。
【請求項１１】
血管機能検査方法であって、
検査側である第１指尖及び第１手首、並びに対照側である第２指尖及び第２手首に予備的に圧力を加え、それぞれの脈波を測定する第１ステップと、
前記圧力が被験者にとって最適圧力であるかどうかを判定し、ＹＥＳならば前記最適圧力を維持する第２ステップと、
測定された駆血前の各脈波に基づいて振幅値を計算し、前記第１指尖、前記第１手首、前記第２指尖及び前記第２手首に係る駆血前のカオスアトラクターを算出する第３ステップと、
検査側の前記第１指尖及び前記第１手首に対し、所定時間の間、駆血を施す第４ステップと、
測定された駆血後の各脈波に基づいて振幅値を計算し、前記第１指尖、前記第１手首、前記第２指尖及び前記第２手首に係る駆血後のカオスアトラクターを算出する第５ステップと、
前記第５ステップにおいて計算された振幅値及び算出されたカオスアトラクターについて、検査側における駆血後の数値を（Ａ）と、同じく駆血前の数値を（Ｂ）と、対照側における駆血後の数値を（Ｃ）と、同じく駆血前の数値を（Ｄ）とするとき、検査側駆血前後比（Ａ／Ｂ）、対照側駆血前後比（Ｃ／Ｄ）及び前記検査側駆血前後比の前記対照側駆血前後比に対する検査側／対照側比（（Ａ／Ｂ）／（Ｃ／Ｄ））を求める第６ステップと、
前記検査側駆血前後比、前記対照側駆血前後比及び前記検査側／対照側比をそれぞれの基準値と対比する第７ステップと、
前記第７ステップにおける対比の結果から心血管疾患と血管内皮機能障害に関する情報を出力する第８ステップと、を備えることを特徴とする血管機能検査方法。
【請求項１２】
前記第５ステップの後に、前記算出されたカオスアトラクターを心血管疾患ごとに既知のカオスアトラクター群とパターンマッチングにより対比する第９ステップをさらに備えることを特徴とする請求項１１に記載の血管機能検査方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、血管機能検査装置に関する。
【背景技術】
【０００２】
重大な心血管疾患については、これに先立ち、血管機能障害が先行することが知られている。したがって、血管機能は、血管リスクを最も早く把握できる指標であり、心血管疾患に至る前に、投薬や生活習慣の変更によって回復可能である。つまり、血管機能の低下を早期に発見できれば、回復不可能な疾患の状態へ移る前に予防するができることになる。
【０００３】
血管機能検査は、侵襲的方途と非侵襲的方法とに大別される。侵襲的方法としては、カテーテルによる血管造影法がある。しかし、これは入院検査であって高コスト、高リスクであり、非侵襲的方法が好ましく用いられる。非侵襲的方法としては、（１）エコー画像による血流依存性血管拡張反応（ＦＭＤ；Ｆｌｏｗ−Ｍｅｄｉａｔｅｄ  Ｄｉｌａｔｉｏｎ）検査がある。手動操作でエコー探触子を血管に押し当て、血管径を測定するものであるが、再現性が低く臨床応用に限界がある。（２）また、末梢動脈圧計側（ＰＡＴ；Ｐｅｒｉｐｈｅｒａｌ  Ａｒｔｅｒｉｓｌ  Ｔｏｎｏｍｅｔｒｙ）検査がある。上腕駆血開放後の血管拡張反応を指尖の末梢脈拍動の血流量の変化として検出するものであるが、検査用プローブが使い捨てのためコストが高い。また、上腕の駆血は血圧計を用いて手動で行うため、操作者の負担が大きい。以下、（１）ＦＭＤ検査と、（２）ＰＡＴについて、説明する。
【０００４】
（１）ＦＭＤ検査とは、血管がアセチルコリン受容体からのシグナル伝達により、アルギニンから合成された一酸化窒素（ＮＯ）を介して拡張することや、血管が機械的刺激や「ずり応力（ｓｈｅａｒ  ｓｔｒｅｓｓ）」にも応答して血管内皮細胞よりＮＯが産生されて拡張することを応用したものである。近年、駆血開放後に血管内の血流増加により強力なずり応力が生じ、その後、一過性に血管が拡張する現象が明らかにされ、ＦＭＤと命名された。すなわち、ＦＭＤはＮＯ依存性の一過性血管拡張反応と理解されている。ＦＭＤにより血管内皮細胞から産生されるＮＯ等の血管拡張物質に反応する血管拡張性が診断できる。ＦＭＤを利用した検査方法としては上腕動脈の血管超音波検査で血管径変化を直接計測する方法を用いる。
【０００５】
しかし、ＦＭＤの上腕動脈拡張率を超音波で測定する機器（ＢＡＵＳ：Ｂｒａｃｈｉａｌ  Ａｒｔｅｒｙ  Ｕｌｔｒａｓｏｕｎｄ）は、超音波による血管内径側定を原理とするため、内径描出が一定とならず、個人差が大きく再現性も乏しい。また、超音波エコー検査者の技量に依存するところが大きく超音波検査機器が高額（数千万円）である問題があった。また、測定時間も１５分間と比較的長いという問題があった。
【０００６】
（２）そこで、最近、ＦＭＤ検査に関して、上腕動脈拡張率をエコーで測定するというデータの不安定と高額なエコー機器を使用しない上腕駆血解放後の血管拡張反応を指尖の毛細血管の圧変化でとらえるＰＡＴ技術による新しいＦＭＤ測定法が開発されている。例えば、非特許文献１に、上腕動脈駆血解放後の血管拡張反応に対して、指尖の血流量変動を圧変化として検出する方法が開示されている。また、ＰＡＴ技術に関して特許文献１〜３が存在する。ＰＡＴでも動脈硬化の判定が可能である報告もある（三上正俊、鍵谷昭文、小澤禎治：加圧脈波による動脈硬化診断。日本臨床生理学会雑誌３５：９−１８，２００５）。
【０００７】
しかし、ＰＡＴ技術によりこれまでに商品化されている機器（商品名「エンドパッド２０００」。なお、「ＥＮＤＯ  ＰＡＴ」は登録商標）は、上腕の駆血部分が手動になっており、加圧が一定でなく不安定で、駆血できているかの判断は、検者の波形による主観的判断に依存し、手動で適宣、再加圧している。指尖脈圧の駆血前後の指尖血流量を脈圧の比率（駆血後／駆血前）で計算するため、変化率が大きく出る。また、指尖血管圧測定プローブは、指尖を覆うプローブ内膜への依存性が高く、１回限りの「使い捨て」であり、ランニングコストが非常に高い問題を有している。また、測定時間も１５分間と、ＦＭＤ検査のそれから短縮されていないという問題があった。
【０００８】
さらには次のような課題もあった。一般に、指尖脈波の振幅は、指尖部をセンサーに圧着させる圧力によって変化する。最適な圧力でセンサーに圧着させると、大きな振幅のきれいな脈波をセンシングすることができるが、圧力が弱すぎると、振幅が小さくノイズ成分が大きくなり、圧力が強すぎると、逆に血流が阻害されて、振幅が小さくなる。最適な圧力は、指の太さや、血管の硬さにより異なるので、バネ式の一般のプローブでは、圧力が弱すぎたり、強すぎたりするために、脈波が精度良く取れない場合がある。また、バネ式のプローブでは、１５分間の測定を続けて指先を圧迫し続けるために、脈波の振幅が減衰してしまうこともある。指尖を内膜による覆うプローブも、多少の空気の漏れがあるため、時間経過に連れて脈波の振幅の減衰がある。
【０００９】
さらにまた、ＦＭＤやＰＡＴでは、測定結果の精度や再現性が低いという課題があった。血管の充血反応原理においては、図１に示すように、血流のフローによるずり力Ｐ１はＮＯを誘発し、ＮＯは、血管を構成する平滑筋の拡張を生じさせるのに対し、血管壁圧Ｐ２の変化は非ＮＯ依存要素を誘発し、この非ＮＯ依存要素は、血管を構成する平滑筋の変化（収縮・膨張）につながる。ＦＭＤやＰＡＴでは、非ＮＯ依存要素による血管壁圧変化の影響を避けることが困難という課題があることと、ずり力Ｐ１及び血管壁圧Ｐ２による複雑で非均一な血管変形から、測定結果の精度や再現性を高らしめるためには、非ＮＯ依存要素による血管壁圧の変化に影響されない測定手段が望まれていた。またさらに、太い血管の測定とするＦＭＤは太い血管に関わる血管内皮機能の検査が得意とする。一方、細い血管の測定とするＰＡＴは細い血管に関わる血管内皮機能の検査が得意とする。血管機能検査の汎用性が限定されていた。そこで、太い血管と細い血管の両方ともの血管機能の検査でき、更に同時観測による両者の割合を定量化できることにより、血管機能障害の病因を突き止められるような測定手段が望まれていた。
【先行技術文献】
【特許文献】
【００１０】
【特許文献１】特表２０００−５１５７８９号公報
【特許文献２】特表２００３−５２７１４９号公報
【特許文献３】特表２００４−５２８０５２号公報
【非特許文献】
【００１１】
【非特許文献１】Bonetti PO,Pumper GM,Higano ST,et al:Non-invasive Identification ofpatients with early colonary atherosclerosis by assessment of digitasl reacyivehyperemia.J Am Coll cardiol 44:2137-2141,2004
【発明の概要】
【発明が解決しようとする課題】
【００１２】
本発明は、上記の課題に鑑みなされたものであり、血管内皮機能測定に適用できる高い再現性と凡用性に優れ、測定時間の短縮化、測定に要するコストの低減化を図った指尖及び手首（橈骨）の脈波解析装置の提供を目的とする。
【課題を解決するための手段】
【００１３】
本発明は、上記目的を達成するために、以下の構成によって把握される。
（１）本発明の第１の観点は、血管機能検査装置であって、検査側の上腕に装着され、駆血及び開放を行う加圧部と、指尖又は手首に装着され、駆血前後の脈波を測定するプローブと、前記加圧部及び前記プローブを制御する制御部と、前記制御部に指令を与え、前記制御部から前記脈波に係るデータを取得する処理部と、を備え、前記処理部が、前記脈波に係るデータに基づいて平均振幅値の計算及びカオス解析を行い、心血管疾患と血管内皮機能障害の判定を行うために前記脈波の特徴量を定量化することを特徴とする。
【００１４】
（２）上記（１）の構成において、前記処理部は、前記カオス解析として、前記脈波の振幅値及びその変化率に基づいて駆血前後の血管ゆらぎを示すカオスアトラクターを算出し、前記算出されたカオスアトラクターについて、リアプノフ指数、外周面積及び短径と長径の固有ベクトルの少なくとも１つを定量化してもよい。
【００１５】
（３）上記（１）又は（２）の構成において、前記プローブは、光電式又は圧力式によって前記脈波を測定してもよい。
【００１６】
（４）上記（３）の構成において、前記プローブは、光電式であって、測定部位に最適圧力を与えるための空気袋と、前記測定部位に光を照射する発光部と、前記測定部位からの反射光を受光する受光部と、を備えてもよい。
【００１７】
（５）上記（１）から（４）のいずれか１つの構成において、前記プローブは、検査側の第１指尖に装着され、細い血管に係る前記脈波を測定する第１プローブを含んでもよい。
【００１８】
（６）上記（５）の構成において、前記プローブは、検査側の第１手首に装着され、太い血管に係る前記脈波を測定する第２プローブをさらに含んでもよい。
【００１９】
（７）上記（６）の構成において、前記プローブは、対照側の第２指尖に装着され、細い血管に係る前記脈波を測定する第３プローブと、対照側の第２手首に装着され、太い血管に係る前記脈波を測定する第４プローブと、をさらに含んでもよい。
【００２０】
（８）上記（７）の構成において、前記処理部が、前記脈波について、検査側における駆血後の数値を（Ａ）と、同じく駆血前の数値を（Ｂ）と、対照側における駆血後の数値を（Ｃ）と、同じく駆血前の数値を（Ｄ）とするとき、検査側駆血前後比（Ａ／Ｂ）、対照側駆血前後比（Ｃ／Ｄ）及び前記検査側駆血前後比の前記対照側駆血前後比に対する検査側／対照側比（（Ａ／Ｂ）／（Ｃ／Ｄ））を求めてもよい。
【００２１】
（９）上記（２）から（８）のいずれか１つの構成において、前記心血管疾患ごとに既知のカオスアトラクター群を記憶する記憶部をさらに備え、前記処理部は、前記算出されたカオスアトラクターを前記既知のカオスアトラクター群とパターンマッチングにより対比してもよい。
【００２２】
（１０）上記（１）から（９）のいずれか１つの構成において、前記処理部によって処理された結果を表示する表示部をさらに備えてもよい。
【００２３】
（１１）本発明の第２の観点は、血管機能検査方法であって、検査側である第１指尖及び第１手首、並びに対照側である第２指尖及び第２手首に予備的に圧力を加え、それぞれの脈波を測定する第１ステップと、前記圧力が被験者にとって最適圧力であるかどうかを判定し、ＹＥＳならば前記最適圧力を維持する第２ステップと、測定された駆血前の各脈波に基づいて振幅値を計算し、前記第１指尖、前記第１手首、前記第２指尖及び前記第２手首に係る駆血前のカオスアトラクターを算出する第３ステップと、検査側の前記第１指尖及び前記第１手首に対し、所定時間の間、駆血を施す第４ステップと、測定された駆血後の各脈波に基づいて振幅値を計算し、前記第１指尖、前記第１手首、前記第２指尖及び前記第２手首に係る駆血後のカオスアトラクターを算出する第５ステップと、前記第５ステップにおいて計算された振幅値及び算出されたカオスアトラクターについて、検査側における駆血後の数値を（Ａ）と、同じく駆血前の数値を（Ｂ）と、対照側における駆血後の数値を（Ｃ）と、同じく駆血前の数値を（Ｄ）とするとき、検査側駆血前後比（Ａ／Ｂ）、対照側駆血前後比（Ｃ／Ｄ）及び前記検査側駆血前後比の前記対照側駆血前後比に対する検査側／対照側比（（Ａ／Ｂ）／（Ｃ／Ｄ））を求める第６ステップと、前記検査側駆血前後比、前記対照側駆血前後比及び前記検査側／対照側比をそれぞれの基準値と対比する第７ステップと、前記第７ステップにおける対比の結果から心血管疾患と血管内皮機能障害を判定する第８ステップと、を備えることを特徴とする。
【００２４】
（１２）上記１１）の構成において、前記第５ステップの後に、前記算出されたカオスアトラクターを心血管疾患ごとに既知のカオスアトラクター群とパターンマッチングにより対比する第９ステップをさらに備えてもよい。
【発明の効果】
【００２５】
本発明によれば、血管内皮機能測定に適用できる高い再現性と凡用性に優れ、測定時間の短縮化、測定に要するコストの低減化を図った指尖及び手首（橈骨）の脈波解析装置を提供することができる。
【図面の簡単な説明】
【００２６】
【図１】血管の充血反応原理を示す図である。
【図２】本実施形態に係る血管機能検査装置の構成例を示す図である。
【図３】図２のうち、プローブ及びカフについて拡大して示す図である。
【図４】本実施形態に係る血管機能検査装置の制御例を示す図である。
【図５】本実施形態に係る血管機能検査装置による測定、解析の手順を示すフローチャートである。
【図６】脈波の振幅、プローブを指尖又は手首に圧着させる装着圧力、血管のコンプライアンスＣとの関係を示す図である。
【図７】検査側と対照側における脈波の測定データ（振幅）の例を示す図である。
【図８】駆血前後のカオスアトラクターの例を示す図である。
【図９】カオスアトラクターのパターンと、血管内皮障害／変形と、病因との対応例を示す図である。
【図１０】カオスアトラクターについてのターケンスの手法を示す第１の図である。
【図１１】カオスアトラクターについてのターケンスの手法を示す第２の図である。
【図１２】カオスアトラクターのリアプノフ指数を考え方を示す図である。
【図１３】カオス解析結果を表示した図であって、心筋梗塞、脳梗塞、冠状動脈硬化などの病因を有する被験者の解析結果の例を示す図である。
【図１４】カオス解析結果を表示した図であって、糖尿病、桃源病、末梢疾患などの病因を有する被験者の解析結果の例を示す図である。
【図１５】本実施形態とエンドパット２０００との対比試験の結果を示す図である。
【発明を実施するための形態】
【００２７】
＜血管機能検査装置１の全体構成＞
本発明の実施形態に係る血管機能検査装置１について、図を参照しつつ説明する。図２は、本実施形態の構成例を示しており、検査側である第１上肢Ａ１の近位側すなわち上腕には駆血するためのカフ１０（加圧部）が装着され、検査側である第１上肢Ａ１及び対照側である第２上肢Ａ２の遠位側すなわち両手指尖及び両手首にはそれぞれの脈波を測定する４つのプローブ２０が、それぞれ装着されている。プローブ２０は、第１上肢Ａ１の第１指尖Ｆ１に装着された第１プローブ２１と、第１上肢Ａ１の第１手首Ｗ１に装着された第２プローブ２２と、第２上肢Ａ２の第２指尖Ｆ２に装着された第３プローブ２３と、第２上肢Ａ２の第２手首Ｗ２に装着された第４プローブ２４とを含んでいる。
【００２８】
なお、以下の説明において、検査側を駆血側と、対照側を非駆血側と称することもある。また、添付の図面において、検査側を「ｔｅｓｔ」と、対照側を「ｃｏｎｔｒｏｌ」と表示することもある。
【００２９】
各プローブは空気チューブ２１ａ、２２ａ、２３ａ、２４ａを介して、カフ１０は空気チューブ１０ａを介して、これらを制御する装置本体３０（制御部）に接続されており、装置本体３０は、プローブ２０により取得されたデータを解析し、表示するコンピュータ４０（処理部）に接続されている。なお、各空気チューブには、接続された要素同士間の信号を送受信するための配線が施されている。また、コンピュータ４０は、例えばスマートホンやタブレットなどの携帯情報端末としてもよく、装置本体３０との接続は、有線でも無線でもいずれでも構わない。コンピュータ４０は、データベース４１（記憶部）及びディスプレイ４２（表示部）を備えているが、これらはコンピュータ４０の一部であっても外付けであってもいずれでも構わない。
【００３０】
図３は、カフ１０の装着と、プローブ２０の装着とを拡大して示す図であり、装置本体３０及びコンピュータ４０との接続は省略されている。プローブ２０と装置本体３０及びコンピュータ４０との関係は、図４を参照して後述する。図３に示すように、両手指尖にはそれぞれ第１プローブ２１及び第３プローブ２３が、両手首には第２プローブ２２及び第４プローブ２４が装着されている。第１プローブ２１と第３プローブ２３、第２プローブ２２と第４プローブ２４はそれぞれ同じ構成（ただし、第１プローブ２１と第２プローブ２２は検査側、第３プローブ２３と第４プローブ２４は対照側）であるので、第１プローブ２１と第２プローブ２２に代表させて、以下説明する。
【００３１】
第１プローブ２１は、検査側の第１指尖Ｆ１に装着され、その背側に位置する空気袋２０１と、第１指尖Ｆ１の掌側に向けて光を照射するＬＥＤ２０２（発光部）と、ＬＥＤ２０２と併設され血管からの反射光を受光する受光センサ２０３（受光部）とを備えている。第２プローブ２２は、検査側の第１手首Ｗ１に装着され、その橈骨の掌側に位置する空気袋２０１と、第１手首Ｗ１の橈骨の掌側に向けて光を照射するＬＥＤ２０２（発光部）と、ＬＥＤ２０２と併設され血管からの反射光を受光する受光センサ２０３（受光部）とを備えている。
【００３２】
それぞれの空気袋２０１は、装置本体３０の内部又は外部に設けられている圧力ポンプ２０４及び排気弁２０５によって加圧又は減圧され、その圧力は圧力センサ２０６によって測定される。空気袋２０１は、脈波を測定するに際して被験者にとっての最適圧力を測定部位に付与するものであり、圧力ポンプ２０４、排気弁２０５及び圧力センサ２０６は装置本体３０によって制御される。詳しくは、図４及び図５を参照して後述する。
【００３３】
それぞれのＬＥＤ２０２としては、例えば赤外線発光ダイオードを好適に用いることができ、受光センサ２０３としては、フォトダイオードやフォトトランジスタを好適に用いることができる。ＬＥＤ２０２と受光センサ２０３は第１プローブ２１又は第２プローブ２２の内部で、それらの光軸が交差するように配置されている。ＬＥＤ２０２から出た光（赤外線）は、皮膚を通して血管に照射されるが、血流中のヘモグロビンにて一部の波長帯（近赤外線）が吸収され、その反射光量は、心拍による血管の容量変動（すなわち、容量脈波）に応じて、変化することとなる。受光センサ２０３は、この変化を捉え、電圧変化として出力する。
【００３４】
なお、ここでは、第１プローブ２１又は第２プローブ２２を、ＬＥＤ２０２及び受光センサ２０３を有する光電式の構成によって容量脈波を測定するものとして説明したが、これに限ることなく、圧力式の構成によって圧脈波を測定するようにしてもよい。圧力式の場合は、ＬＥＤ２０２及び受光センサ２０３に代えて、圧力センサ（圧力素子）により、血管の圧力変動を捕捉する。圧力式の場合には、正確な測定のため、圧力センサを動脈直上部に正しく押圧することが肝要である。
【００３５】
駆血するためのカフ１０としては、公知の自動血圧計などに採用されているものを用いることができ、第１プローブ２１又は第２プローブ２２とは別に、圧力ポンプ１０４又は排気弁１０５によって加圧又は減圧され、その圧力は圧力センサ１０６によって測定される。圧力ポンプ１０４、排気弁１０５及び圧力センサ１０６は装置本体３０によって制御される。カフ１０への加圧と開放は全自動化されており、加圧の知能化制御を行われる。カフ１０の加圧はコンピュータ４０によって自動監視されて完全駆血できるまで制御され、過加圧や欠加圧が避けられる。
【００３６】
＜プローブ２０による測定＞
次に、プローブ２０による脈波の測定について、図４及び図５を参照しつつ、第１プローブ２１に代表させて説明する。第３プローブ２３については同様である。第２プローブ２２及び第４プローブについても基本的には同様であるが、測定部位が異なることに伴って相違する点は後述する。
【００３７】
第１プローブ２１の空気袋２０１は、図４に示すように、空気チューブ２１ａを介して、圧力ポンプ２０４、排気弁２０５及び圧力センサ２０６に接続されている。空気袋２０１の圧力は、圧力センサ２０６による検知結果に応じてコンピュータ４０が装置本体３０の制御部３１に与える圧力制御コマンドに基づき、制御部３１が圧力ポンプ２０４と排気弁２０５を制御することによって、調整される。ＬＥＤ２０２は、コンピュータ４０が制御部３１に与える脈波測定コマンドにより制御され、第１指尖Ｆ１に光を照射し、反射光は受光センサ２０３に受光される。
【００３８】
装置本体３０には、制御部３１として、オートゲイン回路、直流成分抽出回路、Ａ／Ｄ（アナログ／デジタル）変換回路、ＵＳＢコントローラ、アナログフィルター、増幅回路、微分回路（２回微分）が設けられている。制御部３１は、コンピュータ４０から受け取った圧力制御コマンドによって空気袋２０１を制御するが、圧力制御コマンドとしては、加圧スピード、加圧開始、現在圧力維持、減圧スピード、減圧開始及び圧力開放が含まれる。また、制御部３１は、コンピュータ４０から受け取った脈波測定コマンドによってＬＥＤ２０２及び受光センサ２０３を制御するが、脈波制御コマンドとしては、測定開始及び測定終了が含まれる。
【００３９】
制御部３１は、圧力センサ２０６及び受光センサ２０３が取得した信号を、圧力時系列データ、加速度脈波時系列デジタルデータ、脈波直流成分時系列デジタルデータ、脈波時系列デジタルデータとして、コンピュータ４０に送信する。なお、加速度脈波は、元の波形である容量脈波を微分回路で２回微分することにより得られる。
【００４０】
これらのデジタルデータを受けたコンピュータ４０は、脈波の振幅値とその変化率に基づいてカオス解析を行い、血管ゆらぎとしてのカオスアトラクターを算出する。そして、そのゆらぎ指数としてリアプノフ指数を計算するとともに、エントロピーの計算、ピーク値（振幅の上限側のピーク値）の検出、高速フーリエ変換（ＦＦＴ）による心拍変動解析（ＨＲＶ解析；Ｈｅａｒｔ  Ｒａｔｅ  Ｖａｒｉａｂｉｌｉｔｙ）を行う。カオス解析については、項を改めて、後述する。
【００４１】
＜測定及び解析の手順＞
上述した構成に基づき、図５に示す手順に従って測定及び解析を行う。まず、本格測定を開始する前に、第１指尖Ｆ１の掌側に受光センサ２０３を圧着させる最適圧力を求めるため、予備測定を行う（ステップＳ１０１）。例えば、第１プローブ２１の空気袋２０１への加圧を、２０ｍｍＨｇから１２０ｍｍＨｇの範囲で変化させながら約１分間加圧し、被験者ごとに最適圧力を見出す（ステップＳ１０２）。
【００４２】
ここで、最適圧力とは、前述したように、脈波の振幅は、測定部位である第１指尖Ｆ１の掌側を受光センサ２０３に圧着させる圧力によって変化することから、圧力が弱すぎて振幅が小さくノイズ成分が大きくなることを避け、また、逆に圧力が強すぎて血流が阻害され振幅が小さくなることを避けるような圧力を指す。また、測定するたびに指尖や手首に加える装着圧力を最適圧力に設定することは、測定結果のバラツキの要素を吸収し、優れた再現性に寄与する。
【００４３】
この点について図示すると、図６に示すようになる。図６は、横軸に時間（秒）をとり、縦軸に脈波の振幅（一点鎖線で表示）、第１プローブ２１の装着圧力（点線で表示）及び血管のコンプライアンスＣ（破線で表示）をとったものであり、装着圧力の上昇につれて、コンプライアンスＣがまず極地点をとる最大コンプライアンス圧力が出現し、続いて脈波の振幅が極地点をとる最大振幅圧力が出現する。そこで、ノイズの影響を最小限に抑え、かつ、圧力が強すぎて血流が阻害されないように、再現性の高いデータを得るには、このコンプライアンスＣが最大になる最大コンプライアンス圧力又は振幅が最大になる最大振幅圧力を、プローブ２０の最適圧力として採用し、本格測定の続行中、維持することが望ましい。
【００４４】
このうち、装着圧力が相対的に低い最大コンプライアンス圧力を用いると、より変化しやすく鋭敏なデータの変化を入手し易い。このようなことから、プローブ２０のうち、指尖における細い血管（末梢血管）の脈波を測定する第１プローブ２１（第３プローブ２３も同じ）としては、最大コンプライアンス圧力を最適圧力とすることが好ましい。ここで、コンプライアンスＣは、血管壁の伸展性であって、血管の容積変化／壁圧変化で表されるが、コンプライアンスＣの極値点（コンプライアンスＣの変化がゼロ）における装着圧力を最適圧力として設定、維持して脈波を測定することにより、非ＮＯ依存要素の影響を除外できるとともに、最大感度での測定が実現する。すなわち、血管内皮機能を反映したＮＯ（一酸化窒素）による充血反応のみによって測定されので、高精度で血管内皮機能を評価することができる。
【００４５】
これに対し、装着圧力が相対的に高い最大振幅圧力を用いると、よりノイズを小さくしたデータを入手し易い。このようなことから、プローブ２０のうち、手首における太い血管（橈骨動脈）の脈波を測定する第２プローブ２２や第４プローブ２４としては、最大振幅圧力を最適圧力とすることが好ましい。
【００４６】
このようにプローブ２０ごとに最適圧力を維持したまま、本格測定を開始し、駆血前の測定を所定の時間行う。測定時間は、概ね２分間である。受光センサ２０３によって測定された電圧変化は、増幅回路にて増幅され、Ａ／Ｄ変換後にコンピュータ４０に送信されて振幅値として計算され記憶されるとともに、表示部に表示される。また、コンピュータ４０は、計算された振幅値とその変化率を元にカオス解析によってカオスアトラクターを算出する（ステップＳ１０３）。カオス解析については、項を改めて、後述する。
【００４７】
次に、検査側をカフ１０で加圧し（ステップＳ１０４）、駆血を所定の時間行う（ステップＳ１０５）。駆血時間は、概ね３分間である。完全な駆血に達したかどうかを判断し、達成していない場合は駆血が達成されるまで加圧を継続する。カフ１０は、図３に示したように、圧力ポンプ１０４、減圧弁１０５及び圧力センサ１０６によって加圧、開放が施される。駆血が完全であるかどうかの判断方法は、第１プローブ２１と第３プローブ２３、又は第２プローブ２２と第４プローブ２４の間で脈波の相関係数を計算し、相関係数が閾値よりも小さい時に駆血されていると判断する。駆血達成後は、圧力の低下を監視、制御しつつ、所定の駆血時間が経過するまで設定圧力を自動的に維持する。
【００４８】
その後、カフ１０を開放して血流を再開し（ステップＳ１０６）、駆血後の測定を所定の時間行う。測定時間は概ね１分間である。これにより、本実施形態における測定時間の合計は、概ね５分間となる。受光センサ２０３によって測定された電圧変化は、増幅回路にて増幅され、Ａ／Ｄ変換後にコンピュータ４０に送信されて振幅値として計算され記憶されるとともに、表示部に表示される。また、コンピュータ４０は、計算された振幅値とその変化率を元にカオス解析によってカオスアトラクターを算出する（ステップＳ１０７）。カオス解析については、項を改めて、後述する。
【００４９】
図７は、プローブ２０のうち、上段に検査側の第１プローブ２１を、下段に対照側の第３プローブ２３によって測定された脈波の振幅を示している。横軸は時間を、縦軸は脈波の振幅値を示す。検査側では駆血により一定時間、脈波が消失している。図中、Ａは検査側の駆血後の測定区間、Ｂは同じく検査側の駆血前の測定区間、Ｃは対照側の駆血後の測定区間、Ｄは同じく対照側の駆血前の測定区間を表している。被験者の血管拡張率は、得られた駆血前後の振幅値から算出された平均振幅値に基づき、次のように設定される。すなわち、駆血側であって検査側の駆血後の平均振幅値を（Ａ）、同じく駆血前の平均振幅値を（Ｂ）、非駆血側であって対照側の駆血後の平均振幅値を（Ｃ）、同じく駆血前の平均振幅値を（Ｄ）とするとき、検査側の駆血前後比（Ａ／Ｂ）と対照側の駆血前後比（Ｃ／Ｄ）を求め、前者（Ａ／Ｂ）の後者（Ｃ／Ｄ）に対する検査側／対照側比（（Ａ／Ｂ）／（Ｃ／Ｄ））を、血管拡張率として求める（ステップＳ１０８）。
【００５０】
この検査側／対照側比について、ＦＭＤとエンドパット２０００との関係において次のようにいうことができる。すなわち、ＦＭＤとエンドパット２０００は検査側及び対照側の各駆血前後比を求めて判定しているが、本実施形態では、最終的に検査側／対照側比によって表されるので、異常な結果は少なくなる。また、前述した再現性について、本実施形態、ＦＭＤ、エンドパット２０００それぞれの方法で被験者１人を対象に３日間各複数回の測定を行ったところ、変動係数として、この順に、４．７％、２２％、１８％という結果が得られている。ここで、検査側／対照側比は、リスクファクター（加齢、血管硬化）に対してよい逆相関関係があることが判明しており、加齢については、検査側／対照側比＝−０．００１９４＊年齢＋２．４０１９、相関係数ｒ＝０．４６、血管硬化については、検査側／対照側比＝−０．００１１２＊血管硬化（血管年齢）＋１．８７２８、相関係数ｒ＝０．５７である。
【００５１】
上述した検査側の駆血前後比（Ａ／Ｂ）、対照側の駆血前後比（Ｃ／Ｄ）及び検査側／対照側比（（Ａ／Ｂ）／（Ｃ／Ｄ））は、平均振幅値から得られる血管拡張率のみならず、次項以降に説明するように、カオス解析によって振幅値及びその変化率から血管ゆらぎとして算出されたカオスアトラクターについて、得られたリアプノフ指数の変化に対してもゆらぎ指数の変化として適用される。さらに、カオスアトラクターをパターン（様子）として可視化し、コンピュータ４０内外のデータベース４１に記憶されている心血管疾患に係る既知のカオスアトラクター群とパターンマッチングにより対比する。このように、各比率の基準値との比較により（ステップＳ１０９）、及び／又は、カオスアトラクターのパターンマッチングにより、血管機能の状態を判定し（ステップＳ１１０）、ひいては心血管疾患を判定する。
【００５２】
＜カオス解析＞
本実施形態では、上述した血管機能検査装置１において、測定された脈波の振幅値とその変化率に基づいて独自のカオス解析技術を駆使し、血管ゆらぎすなわち血管拡張のときに生ずる僅かな変化を分析する。図８に示すように、カオス解析によるカオスアトラクターを血管ゆらぎとして把握し、血管拡張のパターン（様子）を示し、血流のフローによるずり力Ｐ１と血管壁圧Ｐ２による複雑で非均一な血管変形を見えることができ、更に血管機能障害の病因を見極めることを可能とするものである。
【００５３】
図８の上段は、検査側の第１指尖Ｆ１に係る駆血前（血管拡張前）及び駆血後（血管拡張後）のカオスアトラクターを、中段は、第１指尖Ｆ１の振幅を、下段は、対照側の第２指尖Ｆ２の振幅を、それぞれ示している。ここでは、図８の中段に示したように、検査側の脈波の振幅が駆血前に比べて駆血後において大きくなっている例を示しており、これに対応して、駆血後のカオスアトラクターは、駆血前のカオスアトラクターよりもサイズ（面積）が大きく、かつ、各ループ間の幅が広くなっている。駆血前後のカオストラクターを対比することによって、血管ゆらぎの変化を把握することができる。
【００５４】
さらに、駆血後（血管拡張後）のカオスアトラクターのパターンを、心血管疾患ごとに対応してデータベース４１に記憶されている既知のカオスアトラクターとパターンマッチングによって対比することにより、心血管疾患の病因を判定することができる。図９は、駆血後（血管拡張後）のカオスアトラクターのパターンと、血管内皮障害／変形と、病因との対応例を示したものであり、上段は心筋梗塞を、下段は桃源病を示している。
【００５５】
カオス解析について、以下に詳しく説明する。カオスアトラクターとは、一般に、ある力学系がそこに向かって時間発展をする集合のことである。その力学系において、カオスアトラクターに十分近い点から運動するとき、そのカオスアトラクターに十分近いままであり続ける。カオスアトラクターに含まれる軌道は、そのカオスアトラクターの内部にとどまり続けること以外に制限はなく、周期的であったり、カオス的であったりする。本実施形態では、横軸を脈波の振幅値とし、縦軸をその振幅値の変化率（振幅速度）として、プロットした集合である。カオスアトラクターの形状と大きさは、血管集団の有効な血管集団形状変化とサイズとを非線形的に捉えることができる。定量的には、以下の２つの事項を把握する。
【００５６】
１つ目は、形状変化を示すゆらぎ指数として、リアプノフ指数を把握する。リアプノフ指数は、後述するように、力学系においてごく接近した軌道が離れていく度合いを表す量である。このリアプノフ指数は、不安定さ（バラツギさ）を示す数値であり、その値が大きければ不安定でバラツキが多く、小さければ安定であることを示す。ここでいう不安定（バラツキ）とはその力学系に小さな力を加えたときに、大きな状態の変化となって現れることを意味する。
【００５７】
２つ目は、大きさを示す面積として、カオスアトラクターの外周面積を捉える。外周面積は、カオスアトラクターの動いた範囲が全部囲めるような、横軸（ｘ軸）、縦軸（ｙ軸）に平行に辺をとったときの四角形の面積である。算出式は、以下の式で表される。
Ｅｎｖ．Ａｒｅａ＝（ｍａｘ（ｘ）−ｍｉｎ（ｘ））＊（ｍａｘ（ｙ）−ｍｉｎ（ｙ））
【００５８】
検査側の第１指尖Ｆ１の第１プローブ２１及び第１手首Ｗ１の第２プローブ２２を対象としたとき、上述の２つの数値を用いて、血管内皮機能の総合指数Ｉｎｄｅｘを、以下のように表すことができる。
Ｉｎｄｅｘ＝Ｃｓ＊｛Ａ１＊第１指尖の平均振幅値（又は、カオスアトラクターの外周面積）＋Ａ２＊第１指尖のゆらぎ指数｝＋Ｃｌ＊｛Ａ３＊第１手首の平均振幅値（又は、カオスアトラクターの外周面積）＋Ａ４＊第１手首のゆらぎ指数｝
ここで、Ｃｓ：第１指尖への寄与度、Ｃｌ：第１手首への寄与度、Ａ１〜Ａ４：それぞれの依存係数である。
【００５９】
＜＜カオスアトラクターの可視化＞＞
次に、脈波などの生体信号からカオスアトラクターを可視化する方法について、説明する。時系列波形の生体信号からカオスアトラクターの形状を描く方法は、ターケンスの手法と呼ばれ、Ｎ点の状態変数、例えば時系列波形Ｓｔ＝｛Ｓ１，Ｓ２，…，ＳＮ｝をｄ個に復元することを考える。まず、遅れ時間τを用いてｄ個の状態変数の組Ｐｉをつくる。
Ｐｉ＝｛Ｓｉ，Ｓｉ＋τ，Ｓｉ＋２τ，…，Ｓｉ＋（ｄ−１）τ｝
仮に復元する状態変数の個数を３とすれば、その組Ｐｉは次のようになる。
Ｐｉ＝｛Ｓｉ，Ｓｉ＋τ，Ｓｉ＋２τ｝
【００６０】
この操作を１≦ｉ＜Ｎ−（ｄ−１）τの範囲で繰り返すことにより、力学系の構造を近似的に復元することができる。さらに、Ｐｉをｄ次元の状態空間座標とみなし、各座標を結ぶことにより力学系の運動をカオスアトラクターとして表現することができる。ターケンスの手法のイメージを図１０及び図１１に示す。Ｓｉ，Ｓｉ＋τ，Ｓｉ＋２τをそれぞれＸ軸座標、Ｙ軸座標、Ｚ軸座標と考え、図１０の（ｘ１，ｙ１，ｚ１）、（ｘ２，ｙ２，ｚ２）、（ｘ３，ｙ３，ｚ３）を、図１１の左図に示したＸ，Ｙ，Ｚの３軸をもつ３次元状態空間に埋め込む。これを１≦ｉ＜Ｎ−（ｄ−１）τの範囲で繰り返すと図１１の右図のようになる。この手法において１つの状態変数をｄ個に復元することを埋め込みといい、遅れ時間τを埋め込み遅延時間、復元する状態変数の個数ｄ（これは状態空間の次元と考えることもできる）を埋め込み次元と呼ぶ。各々のパラメータは自己相関やフラクタル次元を用いて最適値を求める必要がある。
【００６１】
＜＜リアプノフ指数＞＞
リアプノフ指数は、カオスアトラクターの中に軌道を追跡して、次のように得ることができる。初期時刻で互いに似た状態は、超球サイズの中に似た位置（距離が近い）にある。図１２に示すように、この距離をｒ（０）とし、この軌道を発展して時刻ｔでの距離（軌道の差、離れる度合い）をｒ（ｔ）とすると、
ｒ（ｔ）＝ｒ（０）＊ｅｘｐ（λｔ）
と表せ、このλをリアプノフ指数と呼ぶ。Ｒｏｓｅｎｓｔｅｉｎ（１９９３）のアルゴリズムにより、リアプノフ指数の推定は次のように行われる。近傍軌道間の初期変位をＤ、ｔ秒後の変位をｄ（ｔ）としたとき、リアプノフ指数λは、
ｄ（ｔ）＝Ｄ＊ｅｘｐ（λｔ）
で定義されるので、上式の両辺に対数変換を施し
ｌｏｇ（ｄ（ｔ））＝λｔ＋ｌｏｇ（Ｄ）
により、近傍軌道間変位−時間関数の線形近似直線の傾きとして、リアプノフ指数λを推定する。
【００６２】
＜＜カオスアトラクターのパターンマッチング＞＞
カオスアトラクターのパターンマッチングについては、文字認識などで用いられているテンプレートマッチングという手法を応用して、パターンの類似度を計算する。例えば、計算ファイルのカオスアトラクターと、基本データファイルのカオスアトラクターをそれぞれ、２次元の場合は、１００×１００のセルに分割し、両方のカオスアトラクターのセルに点が存在する場合の数を数える。３次元の場合は、「平面」から「立体」へ展開し、１００×１００×１００のセルに分割し、両方のカオスアトラクターのセルに点が存在する場合の数を数える。
【００６３】
ここで、２つのカオスアトラクターの類似度の計算は、それぞれのカオスアトラクターのみに点が存在するセルの数をＴ１、両方のカオスアトラクターのセルに共通して点が存在するセルの数をＴ２とすると、両者のパターン類似度は、以下の式で表される。
パターン類似度＝（Ｔ２／Ｔ１）×１００
【００６４】
＜＜カオスアトラクターの特徴量を示す他の例＞＞
カオスアトラクターの特徴量を判断する方法は、上述したもののほかに、例えば次のようなものとしてもよい。カオスアトラクターのそれぞれのループの短径Ｄと長径Ｌの特徴長さを計算し、これによって診断を高精度に行うことができる。ベクトルに対して線形変換をし、向きが同じで大きさが変わったベクトルが算出され、ＤとＬを固有ベクトルといい、その値（長さ）を「固有値」と呼ぶ。そして、もっとも大きい固有値に対応する固有ベクトルのことをＬで、２番目に大きい固有値に対応する固有ベクトルをＤとしたとき、
Ｑ＝Ｄ／Ｌ
をもって診断することが可能であり、例えば、Ｑ＝０．７５以上であれば軽度故障、Ｑ＝０．５０−０．７５であれば中度故障、Ｑ＝０．５０以下であれば重度故障などと、診断できる。
【００６５】
＜＜カオス解析結果の表示＞＞
図１３及び図１４は、カオス解析によって得られた結果をディスプレイ４２に表示した例を示している。これらの図では、左側に、上段から下段に向かって、検査側の第１指尖Ｆ１の脈波の振幅、対照側の第２指尖Ｆ２の脈波の振幅、検査側の第１手首Ｗ１の脈波の振幅、対照側の第２手首Ｗ２の脈波の振幅が表示され、中央に、それぞれに対応する駆血前後のカオスアトラクターが表示されている。右側の上段には、解析結果として、指尖及び手首それぞれについて、平均振幅値、最大振幅値−最小振幅値、リアプノフ指数に係る検査側の駆血前後比、対照側の駆血前後比、検査側／対照側比が表示されている。右側の下段には、指尖及び手首それぞれについて、駆血前の測定開始時間、駆血前のベースライン継続時間、駆血後の測定開始時間、駆血後の関心領域（ＲＯＩ）継続時間が表示される。
【００６６】
図１３は、心筋梗塞、脳梗塞、冠状動脈硬化などの病因を有する被験者の解析結果を示しており、この場合、太い血管すなわち手首の橈骨動脈における解析結果（血管内皮機能低下）を有効に用いることができる。これに対し、図１４は、糖尿病、桃源病、末梢疾患などの病因を有する被験者の解析結果を示しており、この場合、細い血管すなわち指尖の末梢血管における解析結果（血管内皮機能低下）を有効に用いることができる。
【００６７】
＜エンドパッド２０００との対比＞
最後に、本実施形態とエンドパッド２０００との対比について、説明する。本実施形態においては、前述したとおり、両指尖及び両手首のプローブ２０は、エンドパット２０００のプローブが１回限りの使い捨てであるのに対し、複数回使用可能とするために、プローブ内膜に依存性のないＬＥＤ／受光センサ式を採用し、最適圧力で指尖又は手首を押さえつける構造としている。これにより、ランニングコストが抑制され、半永久的に使用が可能であって安価なプローブを提供でき、コスト低減が図られる。
【００６８】
また、エンドパット２０００と本実施形態との対比試験を行ったところ、両者の検査結果の間に高い相関が確認できた。具体的には、健常者を被験者として、エンドパット２０００と、本実施形態を反映したプロトタイプ機器による同時測定を行い、両者の血管内皮機能評価指数（ＲＨＩ；Ｒｅａｃｔｉｖｅ  Ｈｙｐｅｒｅｍｉａ  Ｉｎｄｅｘ）を解析した。その結果、図１５に示すように、相関係数Ｒ＝０．６８（有意確率Ｐ＜０．０５）、推定誤差０．０３６の十分な相関が得られ、有意性が認められた。すなわち、本プロトタイプ機器が普及型機器として、輸入品のエンドパット２０００に代わり得ることが確認された。
【００６９】
＜本実施形態の変形例＞
上述した実施形態は、プローブ２０が、検査側の第１指尖Ｆ１に装着される第１プローブ２１、同じく検査側の第１手首Ｗ１に装着される第２プローブ２２、対照側の第２指尖Ｆ２に装着される第３プローブ２３、同じく対照側の第２手首Ｗ２に装着される第４プローブ２４が含まれる場合すなわち４点測定の場合について説明した。しかし、プローブ２０の設定については、その検査目的や被験者の状況などに応じて、４つのプローブ２０を装着した４点測定でなくとも、例えば、１つのプローブ２０を装着した１点測定や、２つのプローブ２０を装着した２点測定などにおいても、測定や解析の構成及び手順について、本発明の要旨を反映させることができる。
【００７０】
１点測定としては、例えば、一方の上肢が不自由な被験者の場合、他方の上肢を検査側として、細い血管（末梢血管）のみを測定することで十分であれば、検査側の第１指尖Ｆ１に装着される第１プローブ２１のみとしたり、太い血管（橈骨動脈）のみを測定することで十分であれば、検査側の第１手首Ｗ１に装着される第２プローブ２２のみとしたりする場合が挙げられる。
【００７１】
２点測定としては、やはり、一方の上肢が不自由な被験者の場合、他方の上肢を検査側として、検査側の第１指尖Ｆ１に装着される第１プローブ２１と検査側の第１手首Ｗ１に装着される第２プローブ２２の組合せとする場合が挙げられる。また、両上肢が自由な被験者の場合、細い血管（末梢血管）のみを測定することで十分であれば、検査側の第１指尖Ｆ１に装着される第１プローブ２１と対照側の第２指尖Ｆ２に装着される第３プローブ２３の組合せとすればよく、太い血管（橈骨動脈）のみを測定することで十分であれば、検査側の第１手首Ｗ１に装着される第２プローブ２２と対照側の第２手首Ｗ２に装着される第４プローブ２４の組合せとすればよい。
【００７２】
１点測定や対照側にプローブ２０を装着しない２点測定の場合は、脈波の振幅値やリアプノフ指数について、対照側の駆血前後比や検査側／対照側比を算出できないが、そのような場合であっても、検査側における駆血前後比や駆血前後のカオスアトラクターによる対比、既知の心血管疾患ごとの既知のカオスアトラクターとの対比は可能であるので、問題なく本発明の要旨を適用できる。
【００７３】
＜本実施形態の効果＞
本実施形態の効果を整理すると、次のとおりである。
（１）測定時間は、ＦＭＤやエンドパット２０００の１５分から概ね５分に大幅に短縮化できる。
（２）検査側の上腕に巻くカフは、加圧と開放の全自動化、加圧の知能化制御を行い、カフの加圧はコンピュータによって自動監視されて完全駆血できるまで制御され、過加圧や欠加圧が避けられる。
（３）指尖や手首に加える装着圧力を、コンプライアンスの極値点（コンプライアンスの変化がゼロ）における装着圧力を最適圧力として設定、維持して脈波を測定することにより、非ＮＯ依存要素の影響を除外できるとともに、最大感度での測定が実現する。すなわち、血管内皮機能を反映したＮＯ（一酸化窒素）による充血反応のみによって測定されので、高精度で血管内皮機能を評価することができる。
（４）ＦＭＤとエンドパット２０００よりも、血管内皮機能は充血による振幅の変化（検査側／対照側比）によって表され、異常な結果は少なくなる。
（５）測定するたびに指尖や手首に加える装着圧力は最適圧力に設定されることで、測定結果のバラツキの要素を吸収でき、再現性に優れる。
（６）検査側／対照側比とリスクファクター（加齢、血管硬化）との間で得られた逆相関関係を活用できる。
（７）ＦＭＤとエンドパット２０００の欠点を克服し、安価でプローブを提供でき、コスト低減が図られ、全自動測定と解析でだれでも使いやすく、一般臨床検査に使える医療機器として提供できる。
（８）太い血管と細い血管の両方ともの血管機能の検査でき、更に同時観測による両者の割合を定量化できることにより、広範に血管機能障害の検査が可能となる。
（９）血管ゆらぎによる血管拡張のパターン（複雑な様子）を把握でき、血管機能障害による疾患の病因を見極められる。
【００７４】
このように、本発明は、上述した具体的な実施形態に限定されるものではなく、種々の変更を行ったものも含まれるものであり、そのことは、当業者にとって特許請求の範囲の記載から明らかである。
【産業上の利用可能性】
【００７５】
本発明は、生体における非浸襲的な直接血管機能を検査できる方法であるため、メタボリック検診事業等に活用できる可能性がある。動脈硬化発症前診断として、血管病変早期発見に加え、メタボロックシンドロームで要指導となった住民の健康管理事業にも貢献できる。また、スタチン系薬剤・抗血小板薬の動脈硬化に対する生体内薬効判定が可能となる。
【符号の説明】
【００７６】
１…血管機能検査装置
１０…カフ（加圧部）
１０４…圧力ポンプ
１０５…排気弁
１０６…圧力センサ
２０…プローブ
２０１…空気袋
２０２…ＬＥＤ（発光部）
２０３…受光センサ（受光部）
２０４…圧力ポンプ
２０５…排気弁
２０６…圧力センサ
２１…第１プローブ
２１ａ…空気チューブ
２２…第２プローブ
２２ａ…空気チューブ
２３…第３プローブ
２３ａ…空気チューブ
２４…第４プローブ
２４ａ…空気チューブ
３０…装置本体
３１…制御部
４０…コンピュータ（処理部）
４１…データベース（記憶部）
４２…ディスプレイ（表示部）
【図１】
0006620333.tif 000002
【図２】
0006620333.tif 000003
【図３】
0006620333.tif 000004
【図４】
0006620333.tif 000005
【図５】
0006620333.tif 000006
【図６】
0006620333.tif 000007
【図７】
0006620333.tif 000008
【図８】
0006620333.tif 000009
【図９】
0006620333.tif 000010
【図１０】
0006620333.tif 000011
【図１１】
0006620333.tif 000012
【図１２】
0006620333.tif 000013
【図１３】
0006620333.tif 000014
【図１４】
0006620333.tif 000015
【図１５】
0006620333.tif 000016