測定対象血管の力学的機能の評価方法、測定対象血管の力学的機能評価装置、測定対象血管の力学的機能の評価プログラム及び記憶媒体
【課題】
測定対象血管が人の左右対称的な位置にそれぞれ位置する血管から取得される振動のデータに基づいて、前記血管の振動に関する固有振動数を抽出することにより、該固有振動数を血管の力学的指標として、前記血管の動脈硬化の評価を行うことができる測定対象血管の力学的機能の評価方法及び評価装置を提供する。
【解決手段】
評価装置１０では、マイク１１Ｌ，１１Ｒ、振動センサ１２Ｌ，１２Ｒが左右の浅側頭動脈７０Ｌ，７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの血流による振動を同時に計測する。評価装置１０の計算処理装置１４は、計測して得られた振動のデータを周波数解析して、各血管の固有振動数を生成する。評価装置１０は、固有振動数とスティフネスとの相関関係、及び前記固有振動数と予想破裂血圧との相関関係の両者の相関関係に基づいて、抽出した固有振動数に対応するスティフネス及び予想破裂血圧を求める。
測定対象血管が人の左右対称的な位置にそれぞれ位置する血管から取得される振動のデータに基づいて、前記血管の振動に関する固有振動数を抽出することにより、該固有振動数を血管の力学的指標として、前記血管の動脈硬化の評価を行うことができる測定対象血管の力学的機能の評価方法及び評価装置を提供する。
【解決手段】
評価装置１０では、マイク１１Ｌ，１１Ｒ、振動センサ１２Ｌ，１２Ｒが左右の浅側頭動脈７０Ｌ，７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの血流による振動を同時に計測する。評価装置１０の計算処理装置１４は、計測して得られた振動のデータを周波数解析して、各血管の固有振動数を生成する。評価装置１０は、固有振動数とスティフネスとの相関関係、及び前記固有振動数と予想破裂血圧との相関関係の両者の相関関係に基づいて、抽出した固有振動数に対応するスティフネス及び予想破裂血圧を求める。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、測定対象血管の力学的機能の評価方法、測定対象血管の力学的機能の評価装置、測定対象血管の力学的機能の評価プログラム及び記憶媒体に関する。
【背景技術】
【０００２】
生活習慣の多様化、世界的高年齢化の進行や食生活の多様化により近年では若年層も含め動脈硬化が進行し、該症状に関連した心筋梗塞、脳梗塞、脳出血など血管系疾患による長期治療患者が増加している。特に、動脈硬化は、自覚症状がほとんどないまま進行するため、動脈硬化を予防するためには定期的な診察が必要である。
【０００３】
現在の頭部血管機能検査は、Ｘ線ＣＴ、ＮＭＲ等の画像化法により形状異常、血管の閉塞状況等を画像化して検出することが主流であり、相当に症状が進行した状態での検査である。これらの装置は、光学な装置であり、しかも大型であるため、設置されている施設が限られている。
【０００４】
また、これらの検査に加えて、超音波エコーを使用して、血管に関して、血管画像、血圧、脈波を計測する装置の技術が公知である（特許文献１〜３参照）。
下記の特許文献１乃至６は、本出願人が調査した出願時の技術水準を示すものである。特許文献１には、医用超音波診断装置で、頸動脈を画像化して、頸動脈径の心拍変動を測定することにより、動脈硬化の指標を得ることが提案されている。特許文献２では、医用超音波診断装置で得られる２次元画像情報により、頸動脈等の動脈を３次元画像にする技術が提案されている。特許文献３では、頸動脈のエコー動画像に基づいて頸動脈の弾性率を計測し、該弾性率と破裂圧力との関係に基づいて頸動脈強度解析を行う生体頸動脈強度解析システムが提案されている。特許文献４では、耳に装着した検出器で血管に流れる脈波を検出して、該脈波に基づき、頸動脈を含めた首からの上の血管における動脈硬化を評価するための指数を測定する装置が提案されている。特許文献５，６では、浅側頭動脈又はその周辺部の血圧測定のためのカフと、血管に流れる血液の脈波を検出するための脈波検出手段を備え、前記カフと脈波検出手段を、耳に取付けることが開示されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特許第３８８２０８４号
【特許文献２】特開２００５−２７０３５１号公報
【特許文献３】特開２００８−７３０８７号公報
【特許文献４】特開２００６−１０２２５１号公報
【特許文献５】特開２００７−６１１８５号公報
【特許文献６】特開２００７−２５９９５７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
なお、特許文献１乃至特許文献３は、医用超音波診断装置を使用するものであり、この装置は、特別な装置であって、設置されている施設が限られている。
特許文献４は、外耳及びその周辺部を測定部位として脈波を検出し、該脈波に基づき評価指数を算出し、該評価指数を提示するようにしている。又、特許文献５，６は、浅側頭動脈の血圧を測定するものである。特許文献４乃至特許文献６は、耳又はその周辺部に脈波、又は血圧を検出するようにしているが、測定対象血管が人の左右対称的な位置にそれぞれ位置する血管の力学的機能（例えば、スティフネス、強度）をどのように評価するかについては開示されていない。
【０００７】
人の頭部血管の力学的機能は、年齢の進行とともに低下しており、定期的に計測する技術の確立が望まれている。特に人の頭部血管は、一般的に左右対称に位置しており、その左右対称に位置した血管の力学的機能は、少年期、青年期、老年期の経過とともに変化する。この左右対称に位置した血管の力学的機能が同時に評価されれば、頭部血管の動脈硬化の予防に役立つことが期待される。
【０００８】
本発明の目的は、測定対象血管が人の左右対称的な位置にそれぞれ位置する血管から取得される振動のデータに基づいて、前記血管の振動に関する固有振動数を抽出することにより、該固有振動数を血管の力学的指標として、前記血管の動脈硬化の評価を行うことができる測定対象血管の力学的機能の評価方法及び評価装置を提供することにある。
【０００９】
又、本発明の他の目的は、測定対象血管が人の左右対称的な位置にそれぞれ位置する血管から取得される振動のデータに基づいて、前記血管の振動に関する固有振動数を抽出することにより、該固有振動数を血管の力学的指標として、前記血管の動脈硬化の評価を行うことができる測定対象血管の力学的機能の評価プログラム及び記憶媒体を提供することにある。
【課題を解決するための手段】
【００１０】
上記目的を達成するために、請求項１に記載の発明は、測定対象血管が人の左右対称的な位置にそれぞれ位置する血管であって、血流による各測定対象血管の振動を同時に計測する第１ステップと、計測して得られた前記振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する第２ステップと、前記測定対象血管における固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係のうち、少なくともいずれか一方の相関関係に基づいて、第２ステップで抽出した固有振動数に対応するスティフネス及び予想破裂血圧のうち少なくともいずれか一方を求める第３ステップを含むことを特徴とする測定対象血管の力学的機能の評価方法を要旨とするものである。
【００１１】
請求項１の発明の測定対象血管の力学的機能の評価方法は、人の左右対称的に位置する血管を測定対象とする。人の左右対称的に位置する血管は、例えば頭部血管の左右浅側頭動脈、或いは左右頸動脈があり、これらの血管の力学的機能は、年齢とともに物理的劣化が不均一に進行する。前記力学的機能は、血管の弾性率（スティフネス）、又は、予想破裂血圧に関係する。このような左右対称的に位置する血管の振動を同時計測し、計測して得られた前記振動のデータを周波数解析すると、測定対象血管の振動の固有振動数が得られる。このことにより、左右の測定対象血管の力学的機能の評価が可能となる。
【００１２】
請求項２の発明は、請求項１において、前記第２ステップでは、さらに、前記振動のデータに基づいて、前記各測定対象血管の振動をリサージュ波形で表示することを特徴とする。
【００１３】
請求項３の発明は、請求項１又は請求項２において、前記測定対象血管が人の左右対称的な位置にそれぞれ位置する血管が、左右の浅側頭動脈であることを特徴とする。
請求項４の発明は、測定対象血管が人の左右対称的な位置にそれぞれ位置する血管であって、血流による各測定対象血管の振動を同時に計測する振動検出手段と、前記計測して得られた前記振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する指標抽出手段と、前記測定対象血管における固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係のうち、少なくともいずれか一方の相関関係に基づいて、前記指標抽出手段が抽出した固有振動数に対応するスティフネス及び予想破裂血圧のうち少なくともいずれか一方を求める手段とを備えることを特徴とする測定対象血管の力学的機能評価装置を要旨とするものである。
【００１４】
請求項５の発明は、請求項４において、前記各測定対象血管の振動のデータに基づいて、リサージュ波形で表示する表示手段を備えることを特徴とする。
請求項６の発明は、コンピュータを、人の左右対称的な位置にそれぞれ位置する各測定対象血管の血流による振動のデータであって、同時に計測して得られた前記振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する指標抽出手段と、前記測定対象血管における固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係のうち、少なくともいずれか一方の相関関係に基づいて、前記指標抽出手段が抽出した固有振動数に対応するスティフネス及び予想破裂血圧のうち少なくともいずれか一方を求める手段として機能させることを特徴とする測定対象血管の力学的機能の評価プログラムを要旨とするものである。
【００１５】
請求項７の発明は、請求項６において、コンピュータを、前記各測定対象血管の振動のデータに基づいて、リサージュ波形で表示する表示手段として機能させることを特徴とする。
【００１６】
請求項８の発明は、請求項６又は請求項７に記載の測定対象血管の力学的機能の評価プログラムを記憶した記憶媒体を要旨とするものである。
【発明の効果】
【００１７】
請求項１の発明によれば、測定対象血管が人の左右対称的な位置にそれぞれ位置する血管から取得される振動のデータに基づいて、前記血管の振動に関する固有振動数を抽出することにより、該固有振動数を血管の力学的指標として、前記血管の動脈硬化の評価を行うことができる測定対象血管の力学的機能の評価方法を提供できる。又、請求項１の発明によれば、超音波装置、Ｘ線ＣＴ、ＮＭＲと異なり、特別な大型の装置は必要でなく、血管の振動を検出する装置を使用することができ、短時間で、その評価を行うことができる。
【００１８】
請求項２の発明によれば、測定対象血管の振動のリサージュ波形で表示されるため、左右の血管の物理量の差を容易に評価することができる測定対象血管の力学的機能の評価方法を提供できる。
【００１９】
請求項３の発明によれば、左右の浅側頭動脈に対して、請求項１又は請求項２の効果を奏することができる。
請求項４の発明によれば、測定対象血管が人の左右対称的な位置にそれぞれ位置する血管から取得される振動のデータに基づいて、前記血管の振動に関する固有振動数を抽出することにより、該固有振動数を血管の力学的指標として、前記血管の動脈硬化の評価を行うことができる測定対象血管の力学的機能の評価装置を提供できる。又、請求項４の発明によれば、超音波装置、Ｘ線ＣＴ、ＮＭＲと異なり、特別な大型の装置は必要でなく、血管の振動を検出する装置を使用することができ、短時間で、その評価を行うことができる。
【００２０】
請求項５の発明によれば、測定対象血管の振動のリサージュ波形で表示されるため、左右の血管の物理量（測定対象血管の振動の位相、振幅、角振動数）の差を容易に評価することができる測定対象血管の力学的機能の評価装置を提供できる。
【００２１】
請求項６の発明によれば、コンピュータを指標抽出手段及び指標抽出手段が抽出した固有振動数に対応するスティフネス及び予想破裂血圧のうち少なくともいずれか一方を求める手段として機能させることができ、測定対象血管の力学的機能の評価装置の構築を容易に行うことができる評価プログラムを提供できる。
【００２２】
請求項７の発明によれば、請求項６の効果に加えて、コンピュータを、前記各測定対象血管の振動のデータに基づいて、リサージュ波形で表示する表示手段として機能させることができるため、測定対象血管の力学的機能の評価装置の構築を容易に行うことができる評価プログラムを提供できる。
【００２３】
請求項８の発明によれば、請求項６、又は請求項７の評価プログラムを記憶することにより、コンピュータを請求項６又は請求項７の手段に容易に機能させることができる。
【図面の簡単な説明】
【００２４】
【図１】（ａ）は、マイク１１Ｌ，１１Ｌの取付け状態を示す説明図、（ｂ）は振動センサの取付け状態を示す説明図、（ｃ）は測定対象血管の力学的機能の評価装置の概略図。
【図２】測定対象血管の血管位置と測定結果の例を示す説明図。
【図３】ＣＰＵが実行する評価プログラムのフローチャート。
【図４】左浅側頭動脈のマイクで測定された拍動４サイクル分の検出信号のチャート。
【図５】右浅側頭動脈のマイクで測定された拍動４サイクル分の検出信号のチャート。
【図６】図４、図５の１サイクル（１周期）分の信号を同じ時間軸で示したチャート。
【図７】図６に示した１サイクル（１周期）分の信号（左右の血管）をフーリエ級数式（ｎ＝１０）で、ディスプレイ２５に表示する際のチャート。
【図８】（ａ）は、２１歳男性の左浅側頭動脈の振動を測定した電位−周波数解析図、（ｂ）は、６１歳男性の左浅側頭動脈の振動を測定した電位−周波数解析図。
【図９】心電と胸部振動音の同時計測のチャート。
【図１０】（ａ）は計測装置１０によって計測された左右の血管の振動波形図、（ｂ）は周波数解析された左右の血管の振動の周波数を示すチャート、（ｃ）はリサージュ波形図、（ｄ）は固有振動数とスティフネスＥｔｈとの関係、固有振動数と予想破裂血圧の関係を示す表示例の説明図。
【図１１】（ａ）は位相差に関するリサージュ波形図、（ｂ）は振幅に関するリサージュ波形図、（ｃ）は角振動数に関するリサージュ波形図。
【図１２】固有振動数とスティフネスＥｔｈとの関係、固有振動数と予想破裂血圧の関係を示す表示例の説明図。
【図１３】（ａ）は、Ｉ音のリサージュ波形図、（ｂ）は、同じくII音のリサージュ波形図、（ｃ）はＩ音とII音の同時のリサージュ波形図。
【図１４】青年者の左右血管のリサージュ波形図。
【図１５】青年者の左右血管のリサージュ波形図。
【図１６】老年者の左右血管のリサージュ波形図。
【図１７】フーリエ級数式化後の左右血管のリサージュ波形図。
【発明を実施するための形態】
【００２５】
以下、本発明の測定対象血管の力学的機能評価装置（以下、単に評価装置という）及び評価方法を具体化した一実施形態を図１〜１７を参照して説明する。
図１（ａ）〜（ｃ）に示すように、評価装置１０は、振動検出手段としての一対のマイク１１Ｌ，１１Ｒ、振動検出手段としての振動センサ１２Ｒ，１２Ｌ、心電図を取得するための電極１３、前記マイク１１Ｌ，１１Ｒ、振動センサ１２Ｌ，１２Ｒ、及び計算処理装置１４を備えている。
【００２６】
マイク１１Ｌ、１１Ｒは、人の左右の耳甲介腔５０Ｌ，５０Ｒにそれぞれ着脱自在に挿入可能な円筒状のケース１１内に、マイク部１１ｍが収納されて固定されている。又、各ケース１１の先端面、又は、先端の周囲には開口１１ａが形成されている。開口１１ａのケース１１に対する設けられる位置は限定されるものではないが、開口１１ａは、音源に対して、すなわち、右の耳であれば右浅側頭動脈７０Ｒに対して、左の耳であれば、左浅側頭動脈７０Ｌに対して、指向が最適となるように各鼓膜に対して垂直、或いは平行となるように設けられていることが好ましい。なお、開口１１ａの音源（振動源）に対する位置がマイク１１Ｌ，１１Ｒの測定対象血管、すなわち、左浅側頭動脈７０Ｌ、右浅側頭動脈７０Ｒからの血管音が最大（すなわち、音圧が最大）となる最適な位置となるように、耳甲介腔５０Ｌ、５０Ｒに対する挿入量の調整、或いは、各ケース１１の軸心の回りの回転角度調整を行う。
【００２７】
マイク１１Ｌ、１１Ｒからの検出信号は増幅器１５，１６により増幅され、Ａ／Ｄ１７，１８によりアナログ信号からデシタル信号に変換されて計算処理装置１４に入力される。
【００２８】
振動センサ１２Ｌ，１２Ｒは、マイク１１Ｌ，１１Ｒが取付けされた人の首１００において、左頸動脈８０Ｌ，右頸動脈８０Ｒの振動が最適に測定できる皮膚表面に粘着テープ等により取着される。なお、振動センサ１２Ｌ，１２Ｒは、首の皮膚表面に対して皮膚表面に密着させることが好ましい。なお、振動センサ１２Ｌ，１２Ｒを皮膚表面に密着させない場合においても、振動センサ１２Ｌ，１２Ｒを皮膚表面に近接して配置すれば、左頸動脈８０Ｌ，右頸動脈８０Ｒの振動が皮膚表面を介して空気の振動として計測できるため、これでもよい。
【００２９】
振動センサ１２Ｌ，１２Ｒからの検出信号は増幅器１９，２０により増幅され、Ａ／Ｄ２１，２２によりアナログ信号からデシタル信号に変換されて計算処理装置１４に入力される。
【００３０】
心電図を取得するための電極１３は、公知であるため、図１（ｃ）では説明の便宜上、１つのみ代表として図示されているが、複数の電極を含む。電極１３が検出した信号（検出信号）は、増幅器２３により増幅され、Ａ／Ｄ２４によりアナログ信号からデシタル信号に変換されて計算処理装置１４に入力される。
【００３１】
計算処理装置１４は、コンピュータからなり、ディスプレイ２５、記憶装置２６、及び通信部２７を備えている。計算処理装置１４は、ＣＰＵ１４ａ、記憶媒体としてのＲＯＭ１４ｂ，ＲＡＭ１４ｃを備えている。ＲＯＭ１４ｂには評価プログラムが格納されている。ＲＡＭ１４ｃは、前記評価プログラムを実行処理する際の作業用メモリである。
【００３２】
記憶装置２６は、前記プログラムを実行処理して得られた結果の各種データを格納する。ディスプレイ２５は、ＣＰＵ１４ａにより、表示制御されて、各種データをグラフ、チャート等により表示する。
【００３３】
前記計算処理装置１４は、指標抽出手段に相当する。又、ディスプレイ２５を備える計算処理装置１４は、表示手段に相当し、記憶装置２６を有する計算処理装置１４は記憶手段に相当する。
【００３４】
計算処理装置１４は、通信部２７を備えている。通信部２７は、有線通信又は無線通信が可能となっており、外部装置であるコンピュータ３０に対して、該コンピュータ３０が備えている通信部３２を介して通信が可能である。コンピュータ３０は、例えば、病歴等を扱うデータセンタに備えられており、評価装置１０で測定された各種データは、通信部２７を介して、コンピュータ３０に送信され、このデータセンタが備えるデータベースに格納されるようにされている。
【００３５】
（実施形態の作用）
次に上記のように構成された評価装置１０において、ＣＰＵ１４ａが実行する評価プログラムの処理を図３を参照して説明する。なお、マイク１１Ｌ、１１Ｒは、予め左右の耳甲介腔５０Ｌ，５０Ｒに装着されるとともに振動センサ１２Ｌ、１２Ｒは、左頸動脈８０Ｌ，右頸動脈８０Ｒの振動が最適に測定できる首１００の皮膚表面に取着されているものとする。又、電極１３は、心電図を得るために測定対象者の胸部等に、予めつけられているものとする。なお、心電図（ＥＣＧ：electrocardiography）は、マイク１１Ｌ、１１Ｒ、振動センサ１２Ｌ，１２Ｒの検出信号を時系列表示する際の基準時間軸として利用する。
【００３６】
図３は、ＣＰＵ１４ａが実行する評価プログラムのフローチャートである。
なお、評価プログラムを実行するに当たり、予め測定対象者の氏名、年齢、性別、測定部位、測定日時等の基礎項目は図示しないキーボードにより入力されているものとする。
【００３７】
ステップ（以下、ステップをＳという）Ｓ１０において、同時に計測した、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ、左頸動脈８０Ｌ，右頸動脈８０Ｒ、電極１３からの検出信号を取得する。この検出時間は数秒でよい。
【００３８】
Ｓ２０では、ＣＰＵ１４ａは計測された各信号の１周期分を取出す。なお、取出す周期は、複数周期分でもよいが、重なり等が生ずるため、１周期分が好ましい。
Ｓ３０及びＳ４０では、フーリエ変換による周波数解析を行う。具体的には、Ｓ３０では、ＣＰＵ１４ａは、Ｓ２０で取出した左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ、左頸動脈８０Ｌ，右頸動脈８０Ｒの各信号を、ディスプレイ２５に下記のフーリエ級数式で表示（ｎ＝１０〜２０項程度まで）する。
【００３９】
ｙ（ｔ）＝A0 ＋ A1ｃｏｓ（ω1ｔ＋φ1） ＋ A2ｃｏｓ（2ω2ｔ＋φ2） ＋ A3ｃｏｓ（3ω3ｔ＋φ3） ＋・・・ ＋ Anｃｏｓ(nωnt＋φn)
ここで、Ａｎは振幅、ωn＝２πｆｎ、ｆは心拍の基本周波数、ｎ＝１，２，３・・・であり、ωは角振動数、φｎは位相差である。又、Ａ０は直流成分である。なお、Ｓ３０における表示の具体例については、後述する。
【００４０】
Ｓ４０では、ＣＰＵ１４ａは、前記フーリエ級数式の各係数（Ａ０，Ａｎ，φｎ）を決定する。
Ｓ５０では、ＣＰＵ１４ａは、左右の血管について、Ｓ４０で決定した各係数を比較することにより、左右の血管（左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ、左頸動脈８０Ｌ，右頸動脈８０Ｒ）毎に取得された周波数の中から、それぞれピーク周波数を固有振動数として抽出する。詳細には、振動の周波数解析を行い、周波数が一番多い帯域の中心（ピーク周波数）が固有振動数として求められている。
【００４１】
ここで、固有振動数ｐは、振動体の質量Ｍ、バネ定数ｋとすると、
ｐ＝√（ｋ／Ｍ）
の関係がある。バネ定数ｋは力／変形量（単位はＮ／ｍ，ニュートン／メータ）である。又、バネ定数ｋは形状により異なり、同じ材料でも、パイプ状、棒状、コイルのように巻いたものでは異なる。従って、本実施形態において、測定対象血管の固有振動数を求めることは、バネ定数に相関するものを求めていることになる。これは、
スティフネスの定義＝応力／ひずみ＝力の大きさ／変形量
と等価となる。
【００４２】
又、合わせて、ＣＰＵ１４ａは左右の血管について、Ｓ４０で決定した各係数を比較することにより、左右の血管（左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ、左頸動脈８０Ｌ，右頸動脈８０Ｒ）の振幅Ａｎ、角振動数、位相差φｎをデータ化する。ここで、固有振動数、振幅Ａｎ、位相、角振動数は、振動に関する物理量に相当する。
【００４３】
Ｓ６０では、動脈硬化の評価をデータ化し、今回の測定日時等の前記基礎項目と関連づけて記憶装置２６に格納する。前記動脈硬化の評価（すなわち、血管の力学的評価）とは、各血管毎に抽出された前記固有振動数と対応するスティフネスＥｔｈを記憶装置２６に格納された固有振動数とスティフネスとの相関マップを参照して、又、各血管毎に抽出された前記固有振動数と対応する予想破裂血圧を固有振動数と予想破裂血圧との相関マップを参照して、それぞれ算出することである。
【００４４】
なお、これらの相関マップは、予め試験データに基づいて作成され、記憶装置２６に格納されている。
図１０（ｄ）は、測定対象血管の力学的機能の指標となる固有振動数と相関するスティフネスとの相関マップ、及び固有振動数と相関する予想破裂血圧との相関マップである。図１０（ｄ）において、縦軸は振動数、横軸はスティフネスＥｔｈ、及び予想破裂血圧である。従って、Ｓ５０で、各血管の固有振動数が抽出されると、前記相関マップを参照して、抽出された固有振動数と、対応するスティフネスＥｔｈ、及び予想破裂血圧が得られる。
【００４５】
なお、Ｓ６０において、さらに、振幅Ａｎ、位相差φｎについての定量的評価を行ってもよい。この場合の定量的評価は、振幅Ａｎ、位相差φｎの大きさによる評価であって、振幅Ａｎ、位相差φｎ毎に、予めテーブル化されて記憶装置２６に格納されており、該テーブルをＣＰＵ１４ａが参照することにより定量的評価を行う。
【００４６】
Ｓ７０では、ＣＰＵ１４ａは、左右の血管（左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ、左頸動脈８０Ｌ，右頸動脈８０Ｒ）の固有振動数、振幅Ａｎ、位相、及び前記データ化した結果をディスプレイ２５に表示する。
【００４７】
又、Ｓ７０において、ディスプレイ２５における前記データ化した結果の表示には、同時計測による左右の血管における振動のリサージュ波形の表示を含み、さらに、各血管の固有振動数から、前記相関マップに基づいて得られる各血管のスティフネス、及び各血管の予想破裂血圧の表示を含む。
【００４８】
Ｓ７０の処理が終了すると、この評価プログラムを終了する。
なお、これらの表示の具体例は後述する。
上記のフローチャートの流れに沿って、本実施形態の評価方法及び表示の具体例について説明する。
【００４９】
図１０（ａ）に示すように、Ｓ１０で取得した例えば、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒの振動波形は、周波数解析された結果、図１０（ｂ）に示すように、Ｓ５０で、各血管毎にピーク周波数を固有振動数ＰHz（Ｒ），ＰHz（Ｌ）として抽出される。
【００５０】
Ｓ２０〜Ｓ６０を経て、図１０（ｃ）に示すように、Ｓ７０において、同時計測による左右の血管における振動のリサージュ波形をディスプレイ２５に表示する。図１０（ｃ）は、左右の血管の振動のデータに基づいてリサージュ波形を得たときの例を示したものである。図１０（ｃ）において、点線は老年者、実線は青年者の場合の例を示している。
【００５１】
本実施形態では、縦軸に右Ａの血管の位相（又は、振幅、振動数）、横軸に左Ｂの血管の位相（又は、振幅、振動数）を取る。
リサージュ波形について説明する。
【００５２】
図１１（ａ）〜（ｃ）には、リサージュ波形がそれぞれ示されている。
ここで右Ａの血管の振動を、ａ０・sin(ω0+φ0)とし、左Ｂの血管の振動をａ1・sin(ω1+φ1)で表わすことができたとする。このとき、図１１（ａ）の位相に関するリサージュ波形は、振幅ａ0，ａ1が同じ、及び角振動数ω0，ω1が同じとしたとき、位相φ0と位相φ1の差（位相差）が大きいと、同図において、破線で示すようになる。又、位相のリサージュ波形は、振幅ａ0，ａ1が同じ、及び角振動数ω0，ω1が同じとしたとき、位相φ0と位相φ1の差（位相差）が小さいと実線で示すようになる。
【００５３】
このことから、位相差が小さい場合は、左右の血管の物理量はほぼ等価であることを意味する。又、位相差が大きい場合は、片方の血管が伝播速度大であって、スティフネスが大きいことを示している。
【００５４】
又、図１１（ｂ）の振幅に関するリサージュ波形は、位相φ0，φ1が同じ、及び角振動数ω0，ω1が同じとしたとき、振幅ａ0，ａ1比（差）が大きいと、同図において、破線で示すようになる。又、振幅のリサージュ波形は、位相φ0，φ1が同じ、及び角振動数ω0，ω1が同じとしたとき、振幅ａ0，ａ1比（差）が同じであると、実線で示すようになる。
【００５５】
このことから、振幅差（比）が小である場合は、左右の血管の物理量はほぼ等価であることを意味する。振幅差（比）が大である場合は、片方の血管が減衰係数が小であり、スティフネスが大きいことを示している。
【００５６】
又、図１１（ｃ）の角振動数に関するリサージュ波形は、振幅ａ0，ａ1が同じ、及び位相φ0，φ1が同じとしたとき、角振動数ω0，ω1比が大きいと、同図において、破線で示すようになる。又、角振動数のリサージュ波形は、振幅ａ0，ａ1が同じ、及び位相φ0，φ1が同じとしたとき、角振動数ω0，ω1比が小であると、実線で示すようになる。
【００５７】
このことから、角振動数比が小である場合は、左右の血管の物理量はほぼ等価であることを意味する。角振動数比が大である場合は、片方の血管の固有振動数が大であり、スティフネスが大きいことを示している。
【００５８】
ここで、さらに、上記の評価プログラムを実行して、左右の血管の測定した際の実際の例を挙げて説明する。
図２は、１〜５ｃｈ（チャンネル）で、測定した血管、及び、心電図である。図２の１ｃｈ及び２ｃｈは、測定対象者の左耳部の左浅側頭動脈７０Ｌ及び右浅側頭動脈７０Ｒで測定した振動音が示されている。又、図２の３ｃｈ及び４ｃｈは、測定対象者の左頸動脈８０Ｌ，右頸動脈８０Ｒで測定した振動が示されている。なお、本測定例では、振動センサ１２Ｌ，１２Ｒを皮膚表面に密着させず、振動センサ１２Ｌ，１２Ｒを皮膚表面に近接配置して、左頸動脈８０Ｌ，右頸動脈８０Ｒの振動を皮膚表面を介して空気の振動を計測した。又、図２において、５ｃｈは、心電図での測定を示す。前記フローチャートのＳ１０において、実際の測定時には、図２で示されるチャートが得られることになる。なお、図２で示されている各チャートは、横軸が時間、縦軸が音圧（又は振動）に比例するパワー（電位）である。
【００５９】
図４、図５は、Ｓ１０において、左血管（左浅側頭動脈）、及び、右血管（右浅側頭動脈）の、マイク１１Ｌ，１１Ｒで測定された拍動４サイクル（周期）分の検出信号のチャートである。図４、図５において横軸は時間、縦軸は音圧である。なお、図９は、Ｓ１０で取得される心電図を示している。なお、説明の便宜上、参考のために、この心電図には前記心電図と同時に計測した胸部振動音も図示している。同図において、胸部振動音のＩ音とは、第１心音のことで、急速な心室収縮期（三尖弁、僧帽弁が閉じて圧力上昇の期間）に血液が心室壁に擦れて生じる音で緊張音とも呼ばれる。同図において、胸部振動音のII音とは、心室収縮期の終わりに大動脈弁、肺動脈弁が閉じる時の音である。なお、これらのＩ音、及びII音は拍動１サイクル中に生ずる音である。
【００６０】
Ｓ２０では、図４、５で示されている信号から１サイクル（１周期）分を取り出す。例えば、図６は、図４、図５の１サイクル（１周期）分の信号を同じ時間軸で示したものである。
【００６１】
図７は、図６に示した１サイクル（１周期）分の信号（左右の血管）をＳ３０でのフーリエ級数式（ｎ＝１０迄）で、ディスプレイ２５に表示したものである。すなわち、図７は、後述する表１で示される係数（振幅Ａｎ，位相差φｎ）で数式化されたものであり、これらの数字により左右血管の振幅、位相差の比較が容易となる。
【００６２】
図８（ａ）は、測定対象者が、２１歳の男性の左耳の左浅側頭動脈の振動音を測定し、その振動音をＳ３０，Ｓ４０で処理して得られた周波数（横軸（対数軸））と電位（縦軸）を図示したものである。同図において、ピーク周波数が固有振動数となる。同図においては、固有振動数は、３１．８Ｈｚである。
【００６３】
図８（ｂ）は、測定対象者が、６１歳の男性の左耳の左浅側頭動脈の振動音を測定し、その振動音をＳ３０，Ｓ４０で処理して得られた周波数（横軸（対数軸））と電位（縦軸）を図示したものである。同図において、ピーク周波数が固有振動数となる。同図においては、固有振動数は、４８Ｈｚである。
【００６４】
ここで、それぞれの左右の血管の固有振動数が得られた場合、青年期では、左右の血管の振動の物理量の差はゼロであり、動脈硬化進行により左右の血管の振動の物理量の差が大きくなるため、上記図８（ａ）では左耳の固有振動数が得られているが、右耳の固有振動数も同様にして得られた。この２１歳の男性の左右の血管の固有振動数の差は、同じく６１歳の男性の左右血管の固有振動数の差よりも小さく、従って、６１歳の男性の左右血管の固有振動数の差が大きいため、動脈硬化が進行していることが分かる。
【００６５】
【表１】
【００６６】
表１は、Ｓ３０，Ｓ４０において得られた左右の血管（左浅側頭動脈、右浅側頭動脈）音の振幅Ａｎ、左右の血管の位相差（ラジアン及び角度）で示したものである。なお、ｎは１０迄としている。
【００６７】
ここでは、ｎ＝５，６，７当たりに強い成分ｃｎを持ち、その時の位相差は、６０−８０（°）となる。これは、右血管と左血管よりも、６０−８０（°）ほど遅れることを意味する。この場合、動脈硬化の進行により、左右血管の振幅差、位相差は大きくなる傾向がある。
【００６８】
図１２は、図８（ｂ）で説明した測定対象者が、６１歳の男性の場合、前記固有振動数が４８Ｈｚであるため、相関マップを参照して得られる予想破裂血圧は６５０ｍｍＨｇであり、４８Ｈｚに対応するスティフネスＥｔｈが得られる。なお、図１２では、説明の便宜上、図１０（ｃ）と異なり、相関マップは、縦軸がスティフネスＥｔｈ、及び予想破裂血圧、横軸が振動数である。
【００６９】
ここで、図９で説明した胸部振動音のリサージュ波形を参考にして図示する。図１３（ａ）は、Ｉ音のリサージュ波形であり、図１３（ｂ）は、同じくII音のリサージュ波形、図１３（ｃ）はＩ音とII音の同時のリサージュ波形である。
【００７０】
図１４は、測定対象者が、２１歳（青年者）の男性の左浅側頭動脈、及び右浅側頭動脈の振動のリサージュ波形であり、説明の便宜上、拍動４サイクル分を表示したものである。図１５は、測定対象者が、２６歳（青年者）の男性の左浅側頭動脈、及び右浅側頭動脈の振動のリサージュ波形であり、説明の便宜上、拍動４サイクル分を表示したものである。図１６は、測定対象者が、６１歳（老年者）の男性の左浅側頭動脈、及び右浅側頭動脈の振動のリサージュ波形であり、説明の便宜上、拍動４サイクル分を表示したものである。
【００７１】
上記のように、左右血管の振動のリサージュ波形を取ると、動脈硬化が少ない青年者と、動脈硬化が大である老年者とに大きくリサージュ波形が異なることが分かる。
なお、図１７は、参考例として、前記表１の例、即ち、フーリエ数式化後に、リサージュ波形としたものである。フーリエ数式化してリサージュ波形とした図１７と、数式化しないでリサージュ波形とした図１３（ｃ）ものと等価で有ることが分かる。
【００７２】
本実施形態によって発揮される効果について、以下に記載する。
（１） 本実施形態の測定対象血管の力学的機能の評価方法は、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒに血流による各測定対象血管の振動を同時に計測する第１ステップ（Ｓ１０）と、計測して得られた前記振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する第２ステップ（Ｓ５０）とを備える。
【００７３】
又、本実施形態の測定対象血管の力学的機能の評価方法は、さらに、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒにおける固有振動数とスティフネスＥｔｈとの相関関係、及び前記固有振動数と予想破裂血圧との相関関係の両者の相関関係に基づいて、Ｓ５０で抽出した固有振動数に対応するスティフネスＥｔｈ及び予想破裂血圧を求める第３ステップ（Ｓ６０）を備える。
【００７４】
この結果、本実施形態の評価方法は、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒから取得される振動のデータに基づいて、前記血管の振動に関する固有振動数を抽出することにより、該固有振動数を血管の力学的指標として、前記血管の動脈硬化の評価を行うことができる測定対象血管の力学的機能の評価方法を提供できる。又、本実施形態の評価方法によれば、超音波装置、Ｘ線ＣＴ、ＮＭＲと異なり、特別な大型の装置は必要でなく、血管の振動を検出する装置を使用することができ、短時間で、その評価を行うことができる。
【００７５】
（２） 本実施形態の測定対象血管の力学的機能の評価方法は、Ｓ５０（第２ステップ）では、さらに、振動のデータに基づいて、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの振動をリサージュ波形で表示する。この結果、本実施形態の測定対象血管の力学的機能の評価方法によれば、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの振動のリサージュ波形で表示されるため、左右の血管の物理量の差を容易に評価することができる。
【００７６】
（３） 本実施形態の評価装置１０は、マイク１１Ｌ，１１Ｒ、振動センサ１２Ｌ，１２Ｒは振動検出手段として、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの血流による振動を同時に計測する。又、評価装置１０の計算処理装置１４は、指標抽出手段として、計測して得られた振動のデータを周波数解析して、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの力学的機能の指標となる固有振動数を抽出する。
【００７７】
そして、計算処理装置１４は、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒにおける固有振動数とスティフネスＥｔｈとの相関マップ（相関関係）、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関マップ（相関関係）に基づいて、前記指標抽出手段が抽出した固有振動数に対応するスティフネスＥｔｈ及び予想破裂血圧を求める。
【００７８】
本実施形態の評価装置１０によれば、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒから取得される振動のデータに基づいて、前記血管の振動に関する固有振動数を抽出することにより、該固有振動数を血管の力学的指標として、前記血管の動脈硬化の評価を行うことができる測定対象血管の力学的機能の評価装置を提供できる。
【００７９】
本実施形態の評価装置１０によれば、超音波装置、Ｘ線ＣＴ、ＮＭＲと異なり、特別な大型の装置は必要でなく、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの振動を検出する装置を使用することができ、短時間で、その評価を行うことができる。
【００８０】
（４） 本実施形態の評価装置１０の計算処理装置１４は、表示手段として、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの振動のデータに基づいて、リサージュ波形で表示する。この結果、評価装置１０によれば、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの振動のリサージュ波形で表示されるため、左右の血管の物理量（測定対象血管の振動の位相、振幅、角振動数）の差を容易に評価することができる。
【００８１】
（５） 本実施形態の評価プログラムは、コンピュータである計算処理装置１４を、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの血流による振動のデータであって、同時に計測して得られた振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する指標抽出手段として機能させる。
【００８２】
又、本実施形態の評価プログラムは、計算処理装置１４を、固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係に基づいて、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒにおける振動の固有振動数に対応するスティフネス及び予想破裂血圧を求める手段として機能させる。
【００８３】
この結果、本実施形態の評価プログラムによれば、コンピュータを指標抽出手段及び指標抽出手段が抽出した固有振動数に対応するスティフネス及び予想破裂血圧を求める手段として機能させることができ、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの力学的機能の評価装置の構築を容易に行うことができる評価プログラムを提供できる。
【００８４】
（６） 本実施形態の評価プログラムは、ディスプレイ２５を備える計算処理装置１４を、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの振動のデータに基づいて、リサージュ波形で表示する表示手段として機能させる。このため、本実施形態の評価プログラムによれば、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの力学的機能の評価装置の構築を容易に行うことができる評価プログラムを提供できる。
【００８５】
（７） 本実施形態のＲＯＭ１４ｂは、記憶媒体として前記評価プログラムを記憶する。この結果、本実施形態のＲＯＭ１４ｂによれば、計算処理装置１４を、指標抽出手段、表示手段、及び、固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係に基づいて、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒにおける振動の固有振動数に対応するスティフネス及び予想破裂血圧を求める手段として容易に機能させることができる。
【００８６】
なお、前記実施形態を次のように変更して構成することもできる。
○ 前記実施形態では、左右の血管として、左浅側頭動脈７０Ｌ、右浅側頭動脈７０Ｒ、左頸動脈８０Ｌ、右頸動脈８０Ｒ、及び心電図の各種検出信号を取得するようにしたが、下記のように変更してもよい。
【００８７】
左浅側頭動脈７０Ｌ、右浅側頭動脈７０Ｒ、及び心電図の各種検出信号を取得してもよい。この場合、左右の浅側頭動脈に対して、前記実施形態の（１）、（２）の作用効果を奏することができる。すなわち、この場合は、前記Ｓ５０では、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒの固有振動数が取得される。又、Ｓ６０では、左浅側頭動脈７０Ｌ，右浅側頭動脈７０ＲのスティフネスＥｔｈ及び予想破裂血圧が得られる。
【００８８】
或いは、左頸動脈８０Ｌ、右頸動脈８０Ｒ、及び心電図の各種検出信号を取得するようにしてもよい。この場合、左右の浅側頭動脈に対して、前記実施形態の（１）、（２）の作用効果を奏することができる。
【００８９】
この場合は、左頸動脈８０Ｌ、右頸動脈８０Ｒに対して、前記実施形態の（１）、（２）の作用効果を奏することができる。すなわち、この場合は、前記Ｓ５０では、左浅側左頸動脈８０Ｌ、右頸動脈８０Ｒの固有振動数が取得される。又、Ｓ６０では、左頸動脈８０Ｌ、右頸動脈８０ＲのスティフネスＥｔｈ及び予想破裂血圧が得られる。
【００９０】
○ 前記実施形態では、固有振動数に基づいて、スティフネスＥｔｈ及び予想破裂血圧を求めるようにしたが、スティフネスＥｔｈ及び予想破裂血圧のいずれか一方のみを求めるようにしてもよい。
【００９１】
○ 前記実施形態では、マイク部１１ｍを収納するケース１１は、円筒状としたが、イヤホーン形状としてもよく、耳甲介腔５０Ｌ，５０Ｒに挿入できる形状であればよい。
○ 前記実施形態では、ＲＯＭ１４ｂを評価プログラムを記憶する記憶媒体としたが、記憶装置２６を評価プログラムを記憶する記憶媒体にしてもよい。或いは、ＣＤや、ＤＶＤ、ＵＳＢメモリ等の記憶媒体に評価プログラムを格納して、ＣＤやＤＶＤ等の記憶媒体を読込みするドライバ装置等にて読込みして評価プログラムを実行するようにしてもよい。
【００９２】
○ 前記実施形態において、マイク１１Ｌ，１１Ｒに代えて振動計としてもよい。
○ 前記実施形態において、振動センサ１２Ｌ，１２Ｒに代えてマイクとしてもよい。
○ 前記実施形態では、振動に関する物理量としては、固有振動数、位相、振幅としたが、振動の減衰率としてもよい、減衰率は、動脈硬化が進行すると、青年期の血管よりも老年期の血管の方が、弾性率が悪くなるため早く減衰する傾向がある。従って、左右血管の振動の減衰率を測定（抽出）し、差を算出（抽出）するようにしてもよい。
【００９３】
○ Ｓ６０において、データ化した結果、及び左右の血管の固有振動数、振幅Ａｎ等をディスプレイ２５に表示する際、これらの左右の血管の差をディスプレイ２５に表示するようにしてもよい。
【符号の説明】
【００９４】
１０…評価装置、１１Ｌ，１１Ｒ…マイク（振動検出手段）、
１２Ｌ，１２Ｒ…振動センサ（振動検出手段）、
１３…電極、１４…計算処理装置（指標抽出手段、表示手段）、
１４ａ…ＲＯＭ（記憶媒体）、２６…記憶装置。
【技術分野】
【０００１】
本発明は、測定対象血管の力学的機能の評価方法、測定対象血管の力学的機能の評価装置、測定対象血管の力学的機能の評価プログラム及び記憶媒体に関する。
【背景技術】
【０００２】
生活習慣の多様化、世界的高年齢化の進行や食生活の多様化により近年では若年層も含め動脈硬化が進行し、該症状に関連した心筋梗塞、脳梗塞、脳出血など血管系疾患による長期治療患者が増加している。特に、動脈硬化は、自覚症状がほとんどないまま進行するため、動脈硬化を予防するためには定期的な診察が必要である。
【０００３】
現在の頭部血管機能検査は、Ｘ線ＣＴ、ＮＭＲ等の画像化法により形状異常、血管の閉塞状況等を画像化して検出することが主流であり、相当に症状が進行した状態での検査である。これらの装置は、光学な装置であり、しかも大型であるため、設置されている施設が限られている。
【０００４】
また、これらの検査に加えて、超音波エコーを使用して、血管に関して、血管画像、血圧、脈波を計測する装置の技術が公知である（特許文献１〜３参照）。
下記の特許文献１乃至６は、本出願人が調査した出願時の技術水準を示すものである。特許文献１には、医用超音波診断装置で、頸動脈を画像化して、頸動脈径の心拍変動を測定することにより、動脈硬化の指標を得ることが提案されている。特許文献２では、医用超音波診断装置で得られる２次元画像情報により、頸動脈等の動脈を３次元画像にする技術が提案されている。特許文献３では、頸動脈のエコー動画像に基づいて頸動脈の弾性率を計測し、該弾性率と破裂圧力との関係に基づいて頸動脈強度解析を行う生体頸動脈強度解析システムが提案されている。特許文献４では、耳に装着した検出器で血管に流れる脈波を検出して、該脈波に基づき、頸動脈を含めた首からの上の血管における動脈硬化を評価するための指数を測定する装置が提案されている。特許文献５，６では、浅側頭動脈又はその周辺部の血圧測定のためのカフと、血管に流れる血液の脈波を検出するための脈波検出手段を備え、前記カフと脈波検出手段を、耳に取付けることが開示されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特許第３８８２０８４号
【特許文献２】特開２００５−２７０３５１号公報
【特許文献３】特開２００８−７３０８７号公報
【特許文献４】特開２００６−１０２２５１号公報
【特許文献５】特開２００７−６１１８５号公報
【特許文献６】特開２００７−２５９９５７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
なお、特許文献１乃至特許文献３は、医用超音波診断装置を使用するものであり、この装置は、特別な装置であって、設置されている施設が限られている。
特許文献４は、外耳及びその周辺部を測定部位として脈波を検出し、該脈波に基づき評価指数を算出し、該評価指数を提示するようにしている。又、特許文献５，６は、浅側頭動脈の血圧を測定するものである。特許文献４乃至特許文献６は、耳又はその周辺部に脈波、又は血圧を検出するようにしているが、測定対象血管が人の左右対称的な位置にそれぞれ位置する血管の力学的機能（例えば、スティフネス、強度）をどのように評価するかについては開示されていない。
【０００７】
人の頭部血管の力学的機能は、年齢の進行とともに低下しており、定期的に計測する技術の確立が望まれている。特に人の頭部血管は、一般的に左右対称に位置しており、その左右対称に位置した血管の力学的機能は、少年期、青年期、老年期の経過とともに変化する。この左右対称に位置した血管の力学的機能が同時に評価されれば、頭部血管の動脈硬化の予防に役立つことが期待される。
【０００８】
本発明の目的は、測定対象血管が人の左右対称的な位置にそれぞれ位置する血管から取得される振動のデータに基づいて、前記血管の振動に関する固有振動数を抽出することにより、該固有振動数を血管の力学的指標として、前記血管の動脈硬化の評価を行うことができる測定対象血管の力学的機能の評価方法及び評価装置を提供することにある。
【０００９】
又、本発明の他の目的は、測定対象血管が人の左右対称的な位置にそれぞれ位置する血管から取得される振動のデータに基づいて、前記血管の振動に関する固有振動数を抽出することにより、該固有振動数を血管の力学的指標として、前記血管の動脈硬化の評価を行うことができる測定対象血管の力学的機能の評価プログラム及び記憶媒体を提供することにある。
【課題を解決するための手段】
【００１０】
上記目的を達成するために、請求項１に記載の発明は、測定対象血管が人の左右対称的な位置にそれぞれ位置する血管であって、血流による各測定対象血管の振動を同時に計測する第１ステップと、計測して得られた前記振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する第２ステップと、前記測定対象血管における固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係のうち、少なくともいずれか一方の相関関係に基づいて、第２ステップで抽出した固有振動数に対応するスティフネス及び予想破裂血圧のうち少なくともいずれか一方を求める第３ステップを含むことを特徴とする測定対象血管の力学的機能の評価方法を要旨とするものである。
【００１１】
請求項１の発明の測定対象血管の力学的機能の評価方法は、人の左右対称的に位置する血管を測定対象とする。人の左右対称的に位置する血管は、例えば頭部血管の左右浅側頭動脈、或いは左右頸動脈があり、これらの血管の力学的機能は、年齢とともに物理的劣化が不均一に進行する。前記力学的機能は、血管の弾性率（スティフネス）、又は、予想破裂血圧に関係する。このような左右対称的に位置する血管の振動を同時計測し、計測して得られた前記振動のデータを周波数解析すると、測定対象血管の振動の固有振動数が得られる。このことにより、左右の測定対象血管の力学的機能の評価が可能となる。
【００１２】
請求項２の発明は、請求項１において、前記第２ステップでは、さらに、前記振動のデータに基づいて、前記各測定対象血管の振動をリサージュ波形で表示することを特徴とする。
【００１３】
請求項３の発明は、請求項１又は請求項２において、前記測定対象血管が人の左右対称的な位置にそれぞれ位置する血管が、左右の浅側頭動脈であることを特徴とする。
請求項４の発明は、測定対象血管が人の左右対称的な位置にそれぞれ位置する血管であって、血流による各測定対象血管の振動を同時に計測する振動検出手段と、前記計測して得られた前記振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する指標抽出手段と、前記測定対象血管における固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係のうち、少なくともいずれか一方の相関関係に基づいて、前記指標抽出手段が抽出した固有振動数に対応するスティフネス及び予想破裂血圧のうち少なくともいずれか一方を求める手段とを備えることを特徴とする測定対象血管の力学的機能評価装置を要旨とするものである。
【００１４】
請求項５の発明は、請求項４において、前記各測定対象血管の振動のデータに基づいて、リサージュ波形で表示する表示手段を備えることを特徴とする。
請求項６の発明は、コンピュータを、人の左右対称的な位置にそれぞれ位置する各測定対象血管の血流による振動のデータであって、同時に計測して得られた前記振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する指標抽出手段と、前記測定対象血管における固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係のうち、少なくともいずれか一方の相関関係に基づいて、前記指標抽出手段が抽出した固有振動数に対応するスティフネス及び予想破裂血圧のうち少なくともいずれか一方を求める手段として機能させることを特徴とする測定対象血管の力学的機能の評価プログラムを要旨とするものである。
【００１５】
請求項７の発明は、請求項６において、コンピュータを、前記各測定対象血管の振動のデータに基づいて、リサージュ波形で表示する表示手段として機能させることを特徴とする。
【００１６】
請求項８の発明は、請求項６又は請求項７に記載の測定対象血管の力学的機能の評価プログラムを記憶した記憶媒体を要旨とするものである。
【発明の効果】
【００１７】
請求項１の発明によれば、測定対象血管が人の左右対称的な位置にそれぞれ位置する血管から取得される振動のデータに基づいて、前記血管の振動に関する固有振動数を抽出することにより、該固有振動数を血管の力学的指標として、前記血管の動脈硬化の評価を行うことができる測定対象血管の力学的機能の評価方法を提供できる。又、請求項１の発明によれば、超音波装置、Ｘ線ＣＴ、ＮＭＲと異なり、特別な大型の装置は必要でなく、血管の振動を検出する装置を使用することができ、短時間で、その評価を行うことができる。
【００１８】
請求項２の発明によれば、測定対象血管の振動のリサージュ波形で表示されるため、左右の血管の物理量の差を容易に評価することができる測定対象血管の力学的機能の評価方法を提供できる。
【００１９】
請求項３の発明によれば、左右の浅側頭動脈に対して、請求項１又は請求項２の効果を奏することができる。
請求項４の発明によれば、測定対象血管が人の左右対称的な位置にそれぞれ位置する血管から取得される振動のデータに基づいて、前記血管の振動に関する固有振動数を抽出することにより、該固有振動数を血管の力学的指標として、前記血管の動脈硬化の評価を行うことができる測定対象血管の力学的機能の評価装置を提供できる。又、請求項４の発明によれば、超音波装置、Ｘ線ＣＴ、ＮＭＲと異なり、特別な大型の装置は必要でなく、血管の振動を検出する装置を使用することができ、短時間で、その評価を行うことができる。
【００２０】
請求項５の発明によれば、測定対象血管の振動のリサージュ波形で表示されるため、左右の血管の物理量（測定対象血管の振動の位相、振幅、角振動数）の差を容易に評価することができる測定対象血管の力学的機能の評価装置を提供できる。
【００２１】
請求項６の発明によれば、コンピュータを指標抽出手段及び指標抽出手段が抽出した固有振動数に対応するスティフネス及び予想破裂血圧のうち少なくともいずれか一方を求める手段として機能させることができ、測定対象血管の力学的機能の評価装置の構築を容易に行うことができる評価プログラムを提供できる。
【００２２】
請求項７の発明によれば、請求項６の効果に加えて、コンピュータを、前記各測定対象血管の振動のデータに基づいて、リサージュ波形で表示する表示手段として機能させることができるため、測定対象血管の力学的機能の評価装置の構築を容易に行うことができる評価プログラムを提供できる。
【００２３】
請求項８の発明によれば、請求項６、又は請求項７の評価プログラムを記憶することにより、コンピュータを請求項６又は請求項７の手段に容易に機能させることができる。
【図面の簡単な説明】
【００２４】
【図１】（ａ）は、マイク１１Ｌ，１１Ｌの取付け状態を示す説明図、（ｂ）は振動センサの取付け状態を示す説明図、（ｃ）は測定対象血管の力学的機能の評価装置の概略図。
【図２】測定対象血管の血管位置と測定結果の例を示す説明図。
【図３】ＣＰＵが実行する評価プログラムのフローチャート。
【図４】左浅側頭動脈のマイクで測定された拍動４サイクル分の検出信号のチャート。
【図５】右浅側頭動脈のマイクで測定された拍動４サイクル分の検出信号のチャート。
【図６】図４、図５の１サイクル（１周期）分の信号を同じ時間軸で示したチャート。
【図７】図６に示した１サイクル（１周期）分の信号（左右の血管）をフーリエ級数式（ｎ＝１０）で、ディスプレイ２５に表示する際のチャート。
【図８】（ａ）は、２１歳男性の左浅側頭動脈の振動を測定した電位−周波数解析図、（ｂ）は、６１歳男性の左浅側頭動脈の振動を測定した電位−周波数解析図。
【図９】心電と胸部振動音の同時計測のチャート。
【図１０】（ａ）は計測装置１０によって計測された左右の血管の振動波形図、（ｂ）は周波数解析された左右の血管の振動の周波数を示すチャート、（ｃ）はリサージュ波形図、（ｄ）は固有振動数とスティフネスＥｔｈとの関係、固有振動数と予想破裂血圧の関係を示す表示例の説明図。
【図１１】（ａ）は位相差に関するリサージュ波形図、（ｂ）は振幅に関するリサージュ波形図、（ｃ）は角振動数に関するリサージュ波形図。
【図１２】固有振動数とスティフネスＥｔｈとの関係、固有振動数と予想破裂血圧の関係を示す表示例の説明図。
【図１３】（ａ）は、Ｉ音のリサージュ波形図、（ｂ）は、同じくII音のリサージュ波形図、（ｃ）はＩ音とII音の同時のリサージュ波形図。
【図１４】青年者の左右血管のリサージュ波形図。
【図１５】青年者の左右血管のリサージュ波形図。
【図１６】老年者の左右血管のリサージュ波形図。
【図１７】フーリエ級数式化後の左右血管のリサージュ波形図。
【発明を実施するための形態】
【００２５】
以下、本発明の測定対象血管の力学的機能評価装置（以下、単に評価装置という）及び評価方法を具体化した一実施形態を図１〜１７を参照して説明する。
図１（ａ）〜（ｃ）に示すように、評価装置１０は、振動検出手段としての一対のマイク１１Ｌ，１１Ｒ、振動検出手段としての振動センサ１２Ｒ，１２Ｌ、心電図を取得するための電極１３、前記マイク１１Ｌ，１１Ｒ、振動センサ１２Ｌ，１２Ｒ、及び計算処理装置１４を備えている。
【００２６】
マイク１１Ｌ、１１Ｒは、人の左右の耳甲介腔５０Ｌ，５０Ｒにそれぞれ着脱自在に挿入可能な円筒状のケース１１内に、マイク部１１ｍが収納されて固定されている。又、各ケース１１の先端面、又は、先端の周囲には開口１１ａが形成されている。開口１１ａのケース１１に対する設けられる位置は限定されるものではないが、開口１１ａは、音源に対して、すなわち、右の耳であれば右浅側頭動脈７０Ｒに対して、左の耳であれば、左浅側頭動脈７０Ｌに対して、指向が最適となるように各鼓膜に対して垂直、或いは平行となるように設けられていることが好ましい。なお、開口１１ａの音源（振動源）に対する位置がマイク１１Ｌ，１１Ｒの測定対象血管、すなわち、左浅側頭動脈７０Ｌ、右浅側頭動脈７０Ｒからの血管音が最大（すなわち、音圧が最大）となる最適な位置となるように、耳甲介腔５０Ｌ、５０Ｒに対する挿入量の調整、或いは、各ケース１１の軸心の回りの回転角度調整を行う。
【００２７】
マイク１１Ｌ、１１Ｒからの検出信号は増幅器１５，１６により増幅され、Ａ／Ｄ１７，１８によりアナログ信号からデシタル信号に変換されて計算処理装置１４に入力される。
【００２８】
振動センサ１２Ｌ，１２Ｒは、マイク１１Ｌ，１１Ｒが取付けされた人の首１００において、左頸動脈８０Ｌ，右頸動脈８０Ｒの振動が最適に測定できる皮膚表面に粘着テープ等により取着される。なお、振動センサ１２Ｌ，１２Ｒは、首の皮膚表面に対して皮膚表面に密着させることが好ましい。なお、振動センサ１２Ｌ，１２Ｒを皮膚表面に密着させない場合においても、振動センサ１２Ｌ，１２Ｒを皮膚表面に近接して配置すれば、左頸動脈８０Ｌ，右頸動脈８０Ｒの振動が皮膚表面を介して空気の振動として計測できるため、これでもよい。
【００２９】
振動センサ１２Ｌ，１２Ｒからの検出信号は増幅器１９，２０により増幅され、Ａ／Ｄ２１，２２によりアナログ信号からデシタル信号に変換されて計算処理装置１４に入力される。
【００３０】
心電図を取得するための電極１３は、公知であるため、図１（ｃ）では説明の便宜上、１つのみ代表として図示されているが、複数の電極を含む。電極１３が検出した信号（検出信号）は、増幅器２３により増幅され、Ａ／Ｄ２４によりアナログ信号からデシタル信号に変換されて計算処理装置１４に入力される。
【００３１】
計算処理装置１４は、コンピュータからなり、ディスプレイ２５、記憶装置２６、及び通信部２７を備えている。計算処理装置１４は、ＣＰＵ１４ａ、記憶媒体としてのＲＯＭ１４ｂ，ＲＡＭ１４ｃを備えている。ＲＯＭ１４ｂには評価プログラムが格納されている。ＲＡＭ１４ｃは、前記評価プログラムを実行処理する際の作業用メモリである。
【００３２】
記憶装置２６は、前記プログラムを実行処理して得られた結果の各種データを格納する。ディスプレイ２５は、ＣＰＵ１４ａにより、表示制御されて、各種データをグラフ、チャート等により表示する。
【００３３】
前記計算処理装置１４は、指標抽出手段に相当する。又、ディスプレイ２５を備える計算処理装置１４は、表示手段に相当し、記憶装置２６を有する計算処理装置１４は記憶手段に相当する。
【００３４】
計算処理装置１４は、通信部２７を備えている。通信部２７は、有線通信又は無線通信が可能となっており、外部装置であるコンピュータ３０に対して、該コンピュータ３０が備えている通信部３２を介して通信が可能である。コンピュータ３０は、例えば、病歴等を扱うデータセンタに備えられており、評価装置１０で測定された各種データは、通信部２７を介して、コンピュータ３０に送信され、このデータセンタが備えるデータベースに格納されるようにされている。
【００３５】
（実施形態の作用）
次に上記のように構成された評価装置１０において、ＣＰＵ１４ａが実行する評価プログラムの処理を図３を参照して説明する。なお、マイク１１Ｌ、１１Ｒは、予め左右の耳甲介腔５０Ｌ，５０Ｒに装着されるとともに振動センサ１２Ｌ、１２Ｒは、左頸動脈８０Ｌ，右頸動脈８０Ｒの振動が最適に測定できる首１００の皮膚表面に取着されているものとする。又、電極１３は、心電図を得るために測定対象者の胸部等に、予めつけられているものとする。なお、心電図（ＥＣＧ：electrocardiography）は、マイク１１Ｌ、１１Ｒ、振動センサ１２Ｌ，１２Ｒの検出信号を時系列表示する際の基準時間軸として利用する。
【００３６】
図３は、ＣＰＵ１４ａが実行する評価プログラムのフローチャートである。
なお、評価プログラムを実行するに当たり、予め測定対象者の氏名、年齢、性別、測定部位、測定日時等の基礎項目は図示しないキーボードにより入力されているものとする。
【００３７】
ステップ（以下、ステップをＳという）Ｓ１０において、同時に計測した、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ、左頸動脈８０Ｌ，右頸動脈８０Ｒ、電極１３からの検出信号を取得する。この検出時間は数秒でよい。
【００３８】
Ｓ２０では、ＣＰＵ１４ａは計測された各信号の１周期分を取出す。なお、取出す周期は、複数周期分でもよいが、重なり等が生ずるため、１周期分が好ましい。
Ｓ３０及びＳ４０では、フーリエ変換による周波数解析を行う。具体的には、Ｓ３０では、ＣＰＵ１４ａは、Ｓ２０で取出した左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ、左頸動脈８０Ｌ，右頸動脈８０Ｒの各信号を、ディスプレイ２５に下記のフーリエ級数式で表示（ｎ＝１０〜２０項程度まで）する。
【００３９】
ｙ（ｔ）＝A0 ＋ A1ｃｏｓ（ω1ｔ＋φ1） ＋ A2ｃｏｓ（2ω2ｔ＋φ2） ＋ A3ｃｏｓ（3ω3ｔ＋φ3） ＋・・・ ＋ Anｃｏｓ(nωnt＋φn)
ここで、Ａｎは振幅、ωn＝２πｆｎ、ｆは心拍の基本周波数、ｎ＝１，２，３・・・であり、ωは角振動数、φｎは位相差である。又、Ａ０は直流成分である。なお、Ｓ３０における表示の具体例については、後述する。
【００４０】
Ｓ４０では、ＣＰＵ１４ａは、前記フーリエ級数式の各係数（Ａ０，Ａｎ，φｎ）を決定する。
Ｓ５０では、ＣＰＵ１４ａは、左右の血管について、Ｓ４０で決定した各係数を比較することにより、左右の血管（左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ、左頸動脈８０Ｌ，右頸動脈８０Ｒ）毎に取得された周波数の中から、それぞれピーク周波数を固有振動数として抽出する。詳細には、振動の周波数解析を行い、周波数が一番多い帯域の中心（ピーク周波数）が固有振動数として求められている。
【００４１】
ここで、固有振動数ｐは、振動体の質量Ｍ、バネ定数ｋとすると、
ｐ＝√（ｋ／Ｍ）
の関係がある。バネ定数ｋは力／変形量（単位はＮ／ｍ，ニュートン／メータ）である。又、バネ定数ｋは形状により異なり、同じ材料でも、パイプ状、棒状、コイルのように巻いたものでは異なる。従って、本実施形態において、測定対象血管の固有振動数を求めることは、バネ定数に相関するものを求めていることになる。これは、
スティフネスの定義＝応力／ひずみ＝力の大きさ／変形量
と等価となる。
【００４２】
又、合わせて、ＣＰＵ１４ａは左右の血管について、Ｓ４０で決定した各係数を比較することにより、左右の血管（左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ、左頸動脈８０Ｌ，右頸動脈８０Ｒ）の振幅Ａｎ、角振動数、位相差φｎをデータ化する。ここで、固有振動数、振幅Ａｎ、位相、角振動数は、振動に関する物理量に相当する。
【００４３】
Ｓ６０では、動脈硬化の評価をデータ化し、今回の測定日時等の前記基礎項目と関連づけて記憶装置２６に格納する。前記動脈硬化の評価（すなわち、血管の力学的評価）とは、各血管毎に抽出された前記固有振動数と対応するスティフネスＥｔｈを記憶装置２６に格納された固有振動数とスティフネスとの相関マップを参照して、又、各血管毎に抽出された前記固有振動数と対応する予想破裂血圧を固有振動数と予想破裂血圧との相関マップを参照して、それぞれ算出することである。
【００４４】
なお、これらの相関マップは、予め試験データに基づいて作成され、記憶装置２６に格納されている。
図１０（ｄ）は、測定対象血管の力学的機能の指標となる固有振動数と相関するスティフネスとの相関マップ、及び固有振動数と相関する予想破裂血圧との相関マップである。図１０（ｄ）において、縦軸は振動数、横軸はスティフネスＥｔｈ、及び予想破裂血圧である。従って、Ｓ５０で、各血管の固有振動数が抽出されると、前記相関マップを参照して、抽出された固有振動数と、対応するスティフネスＥｔｈ、及び予想破裂血圧が得られる。
【００４５】
なお、Ｓ６０において、さらに、振幅Ａｎ、位相差φｎについての定量的評価を行ってもよい。この場合の定量的評価は、振幅Ａｎ、位相差φｎの大きさによる評価であって、振幅Ａｎ、位相差φｎ毎に、予めテーブル化されて記憶装置２６に格納されており、該テーブルをＣＰＵ１４ａが参照することにより定量的評価を行う。
【００４６】
Ｓ７０では、ＣＰＵ１４ａは、左右の血管（左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ、左頸動脈８０Ｌ，右頸動脈８０Ｒ）の固有振動数、振幅Ａｎ、位相、及び前記データ化した結果をディスプレイ２５に表示する。
【００４７】
又、Ｓ７０において、ディスプレイ２５における前記データ化した結果の表示には、同時計測による左右の血管における振動のリサージュ波形の表示を含み、さらに、各血管の固有振動数から、前記相関マップに基づいて得られる各血管のスティフネス、及び各血管の予想破裂血圧の表示を含む。
【００４８】
Ｓ７０の処理が終了すると、この評価プログラムを終了する。
なお、これらの表示の具体例は後述する。
上記のフローチャートの流れに沿って、本実施形態の評価方法及び表示の具体例について説明する。
【００４９】
図１０（ａ）に示すように、Ｓ１０で取得した例えば、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒの振動波形は、周波数解析された結果、図１０（ｂ）に示すように、Ｓ５０で、各血管毎にピーク周波数を固有振動数ＰHz（Ｒ），ＰHz（Ｌ）として抽出される。
【００５０】
Ｓ２０〜Ｓ６０を経て、図１０（ｃ）に示すように、Ｓ７０において、同時計測による左右の血管における振動のリサージュ波形をディスプレイ２５に表示する。図１０（ｃ）は、左右の血管の振動のデータに基づいてリサージュ波形を得たときの例を示したものである。図１０（ｃ）において、点線は老年者、実線は青年者の場合の例を示している。
【００５１】
本実施形態では、縦軸に右Ａの血管の位相（又は、振幅、振動数）、横軸に左Ｂの血管の位相（又は、振幅、振動数）を取る。
リサージュ波形について説明する。
【００５２】
図１１（ａ）〜（ｃ）には、リサージュ波形がそれぞれ示されている。
ここで右Ａの血管の振動を、ａ０・sin(ω0+φ0)とし、左Ｂの血管の振動をａ1・sin(ω1+φ1)で表わすことができたとする。このとき、図１１（ａ）の位相に関するリサージュ波形は、振幅ａ0，ａ1が同じ、及び角振動数ω0，ω1が同じとしたとき、位相φ0と位相φ1の差（位相差）が大きいと、同図において、破線で示すようになる。又、位相のリサージュ波形は、振幅ａ0，ａ1が同じ、及び角振動数ω0，ω1が同じとしたとき、位相φ0と位相φ1の差（位相差）が小さいと実線で示すようになる。
【００５３】
このことから、位相差が小さい場合は、左右の血管の物理量はほぼ等価であることを意味する。又、位相差が大きい場合は、片方の血管が伝播速度大であって、スティフネスが大きいことを示している。
【００５４】
又、図１１（ｂ）の振幅に関するリサージュ波形は、位相φ0，φ1が同じ、及び角振動数ω0，ω1が同じとしたとき、振幅ａ0，ａ1比（差）が大きいと、同図において、破線で示すようになる。又、振幅のリサージュ波形は、位相φ0，φ1が同じ、及び角振動数ω0，ω1が同じとしたとき、振幅ａ0，ａ1比（差）が同じであると、実線で示すようになる。
【００５５】
このことから、振幅差（比）が小である場合は、左右の血管の物理量はほぼ等価であることを意味する。振幅差（比）が大である場合は、片方の血管が減衰係数が小であり、スティフネスが大きいことを示している。
【００５６】
又、図１１（ｃ）の角振動数に関するリサージュ波形は、振幅ａ0，ａ1が同じ、及び位相φ0，φ1が同じとしたとき、角振動数ω0，ω1比が大きいと、同図において、破線で示すようになる。又、角振動数のリサージュ波形は、振幅ａ0，ａ1が同じ、及び位相φ0，φ1が同じとしたとき、角振動数ω0，ω1比が小であると、実線で示すようになる。
【００５７】
このことから、角振動数比が小である場合は、左右の血管の物理量はほぼ等価であることを意味する。角振動数比が大である場合は、片方の血管の固有振動数が大であり、スティフネスが大きいことを示している。
【００５８】
ここで、さらに、上記の評価プログラムを実行して、左右の血管の測定した際の実際の例を挙げて説明する。
図２は、１〜５ｃｈ（チャンネル）で、測定した血管、及び、心電図である。図２の１ｃｈ及び２ｃｈは、測定対象者の左耳部の左浅側頭動脈７０Ｌ及び右浅側頭動脈７０Ｒで測定した振動音が示されている。又、図２の３ｃｈ及び４ｃｈは、測定対象者の左頸動脈８０Ｌ，右頸動脈８０Ｒで測定した振動が示されている。なお、本測定例では、振動センサ１２Ｌ，１２Ｒを皮膚表面に密着させず、振動センサ１２Ｌ，１２Ｒを皮膚表面に近接配置して、左頸動脈８０Ｌ，右頸動脈８０Ｒの振動を皮膚表面を介して空気の振動を計測した。又、図２において、５ｃｈは、心電図での測定を示す。前記フローチャートのＳ１０において、実際の測定時には、図２で示されるチャートが得られることになる。なお、図２で示されている各チャートは、横軸が時間、縦軸が音圧（又は振動）に比例するパワー（電位）である。
【００５９】
図４、図５は、Ｓ１０において、左血管（左浅側頭動脈）、及び、右血管（右浅側頭動脈）の、マイク１１Ｌ，１１Ｒで測定された拍動４サイクル（周期）分の検出信号のチャートである。図４、図５において横軸は時間、縦軸は音圧である。なお、図９は、Ｓ１０で取得される心電図を示している。なお、説明の便宜上、参考のために、この心電図には前記心電図と同時に計測した胸部振動音も図示している。同図において、胸部振動音のＩ音とは、第１心音のことで、急速な心室収縮期（三尖弁、僧帽弁が閉じて圧力上昇の期間）に血液が心室壁に擦れて生じる音で緊張音とも呼ばれる。同図において、胸部振動音のII音とは、心室収縮期の終わりに大動脈弁、肺動脈弁が閉じる時の音である。なお、これらのＩ音、及びII音は拍動１サイクル中に生ずる音である。
【００６０】
Ｓ２０では、図４、５で示されている信号から１サイクル（１周期）分を取り出す。例えば、図６は、図４、図５の１サイクル（１周期）分の信号を同じ時間軸で示したものである。
【００６１】
図７は、図６に示した１サイクル（１周期）分の信号（左右の血管）をＳ３０でのフーリエ級数式（ｎ＝１０迄）で、ディスプレイ２５に表示したものである。すなわち、図７は、後述する表１で示される係数（振幅Ａｎ，位相差φｎ）で数式化されたものであり、これらの数字により左右血管の振幅、位相差の比較が容易となる。
【００６２】
図８（ａ）は、測定対象者が、２１歳の男性の左耳の左浅側頭動脈の振動音を測定し、その振動音をＳ３０，Ｓ４０で処理して得られた周波数（横軸（対数軸））と電位（縦軸）を図示したものである。同図において、ピーク周波数が固有振動数となる。同図においては、固有振動数は、３１．８Ｈｚである。
【００６３】
図８（ｂ）は、測定対象者が、６１歳の男性の左耳の左浅側頭動脈の振動音を測定し、その振動音をＳ３０，Ｓ４０で処理して得られた周波数（横軸（対数軸））と電位（縦軸）を図示したものである。同図において、ピーク周波数が固有振動数となる。同図においては、固有振動数は、４８Ｈｚである。
【００６４】
ここで、それぞれの左右の血管の固有振動数が得られた場合、青年期では、左右の血管の振動の物理量の差はゼロであり、動脈硬化進行により左右の血管の振動の物理量の差が大きくなるため、上記図８（ａ）では左耳の固有振動数が得られているが、右耳の固有振動数も同様にして得られた。この２１歳の男性の左右の血管の固有振動数の差は、同じく６１歳の男性の左右血管の固有振動数の差よりも小さく、従って、６１歳の男性の左右血管の固有振動数の差が大きいため、動脈硬化が進行していることが分かる。
【００６５】
【表１】
【００６６】
表１は、Ｓ３０，Ｓ４０において得られた左右の血管（左浅側頭動脈、右浅側頭動脈）音の振幅Ａｎ、左右の血管の位相差（ラジアン及び角度）で示したものである。なお、ｎは１０迄としている。
【００６７】
ここでは、ｎ＝５，６，７当たりに強い成分ｃｎを持ち、その時の位相差は、６０−８０（°）となる。これは、右血管と左血管よりも、６０−８０（°）ほど遅れることを意味する。この場合、動脈硬化の進行により、左右血管の振幅差、位相差は大きくなる傾向がある。
【００６８】
図１２は、図８（ｂ）で説明した測定対象者が、６１歳の男性の場合、前記固有振動数が４８Ｈｚであるため、相関マップを参照して得られる予想破裂血圧は６５０ｍｍＨｇであり、４８Ｈｚに対応するスティフネスＥｔｈが得られる。なお、図１２では、説明の便宜上、図１０（ｃ）と異なり、相関マップは、縦軸がスティフネスＥｔｈ、及び予想破裂血圧、横軸が振動数である。
【００６９】
ここで、図９で説明した胸部振動音のリサージュ波形を参考にして図示する。図１３（ａ）は、Ｉ音のリサージュ波形であり、図１３（ｂ）は、同じくII音のリサージュ波形、図１３（ｃ）はＩ音とII音の同時のリサージュ波形である。
【００７０】
図１４は、測定対象者が、２１歳（青年者）の男性の左浅側頭動脈、及び右浅側頭動脈の振動のリサージュ波形であり、説明の便宜上、拍動４サイクル分を表示したものである。図１５は、測定対象者が、２６歳（青年者）の男性の左浅側頭動脈、及び右浅側頭動脈の振動のリサージュ波形であり、説明の便宜上、拍動４サイクル分を表示したものである。図１６は、測定対象者が、６１歳（老年者）の男性の左浅側頭動脈、及び右浅側頭動脈の振動のリサージュ波形であり、説明の便宜上、拍動４サイクル分を表示したものである。
【００７１】
上記のように、左右血管の振動のリサージュ波形を取ると、動脈硬化が少ない青年者と、動脈硬化が大である老年者とに大きくリサージュ波形が異なることが分かる。
なお、図１７は、参考例として、前記表１の例、即ち、フーリエ数式化後に、リサージュ波形としたものである。フーリエ数式化してリサージュ波形とした図１７と、数式化しないでリサージュ波形とした図１３（ｃ）ものと等価で有ることが分かる。
【００７２】
本実施形態によって発揮される効果について、以下に記載する。
（１） 本実施形態の測定対象血管の力学的機能の評価方法は、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒに血流による各測定対象血管の振動を同時に計測する第１ステップ（Ｓ１０）と、計測して得られた前記振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する第２ステップ（Ｓ５０）とを備える。
【００７３】
又、本実施形態の測定対象血管の力学的機能の評価方法は、さらに、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒにおける固有振動数とスティフネスＥｔｈとの相関関係、及び前記固有振動数と予想破裂血圧との相関関係の両者の相関関係に基づいて、Ｓ５０で抽出した固有振動数に対応するスティフネスＥｔｈ及び予想破裂血圧を求める第３ステップ（Ｓ６０）を備える。
【００７４】
この結果、本実施形態の評価方法は、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒから取得される振動のデータに基づいて、前記血管の振動に関する固有振動数を抽出することにより、該固有振動数を血管の力学的指標として、前記血管の動脈硬化の評価を行うことができる測定対象血管の力学的機能の評価方法を提供できる。又、本実施形態の評価方法によれば、超音波装置、Ｘ線ＣＴ、ＮＭＲと異なり、特別な大型の装置は必要でなく、血管の振動を検出する装置を使用することができ、短時間で、その評価を行うことができる。
【００７５】
（２） 本実施形態の測定対象血管の力学的機能の評価方法は、Ｓ５０（第２ステップ）では、さらに、振動のデータに基づいて、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの振動をリサージュ波形で表示する。この結果、本実施形態の測定対象血管の力学的機能の評価方法によれば、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの振動のリサージュ波形で表示されるため、左右の血管の物理量の差を容易に評価することができる。
【００７６】
（３） 本実施形態の評価装置１０は、マイク１１Ｌ，１１Ｒ、振動センサ１２Ｌ，１２Ｒは振動検出手段として、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの血流による振動を同時に計測する。又、評価装置１０の計算処理装置１４は、指標抽出手段として、計測して得られた振動のデータを周波数解析して、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの力学的機能の指標となる固有振動数を抽出する。
【００７７】
そして、計算処理装置１４は、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒにおける固有振動数とスティフネスＥｔｈとの相関マップ（相関関係）、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関マップ（相関関係）に基づいて、前記指標抽出手段が抽出した固有振動数に対応するスティフネスＥｔｈ及び予想破裂血圧を求める。
【００７８】
本実施形態の評価装置１０によれば、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒから取得される振動のデータに基づいて、前記血管の振動に関する固有振動数を抽出することにより、該固有振動数を血管の力学的指標として、前記血管の動脈硬化の評価を行うことができる測定対象血管の力学的機能の評価装置を提供できる。
【００７９】
本実施形態の評価装置１０によれば、超音波装置、Ｘ線ＣＴ、ＮＭＲと異なり、特別な大型の装置は必要でなく、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの振動を検出する装置を使用することができ、短時間で、その評価を行うことができる。
【００８０】
（４） 本実施形態の評価装置１０の計算処理装置１４は、表示手段として、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの振動のデータに基づいて、リサージュ波形で表示する。この結果、評価装置１０によれば、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの振動のリサージュ波形で表示されるため、左右の血管の物理量（測定対象血管の振動の位相、振幅、角振動数）の差を容易に評価することができる。
【００８１】
（５） 本実施形態の評価プログラムは、コンピュータである計算処理装置１４を、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの血流による振動のデータであって、同時に計測して得られた振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する指標抽出手段として機能させる。
【００８２】
又、本実施形態の評価プログラムは、計算処理装置１４を、固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係に基づいて、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒにおける振動の固有振動数に対応するスティフネス及び予想破裂血圧を求める手段として機能させる。
【００８３】
この結果、本実施形態の評価プログラムによれば、コンピュータを指標抽出手段及び指標抽出手段が抽出した固有振動数に対応するスティフネス及び予想破裂血圧を求める手段として機能させることができ、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの力学的機能の評価装置の構築を容易に行うことができる評価プログラムを提供できる。
【００８４】
（６） 本実施形態の評価プログラムは、ディスプレイ２５を備える計算処理装置１４を、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの振動のデータに基づいて、リサージュ波形で表示する表示手段として機能させる。このため、本実施形態の評価プログラムによれば、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒの力学的機能の評価装置の構築を容易に行うことができる評価プログラムを提供できる。
【００８５】
（７） 本実施形態のＲＯＭ１４ｂは、記憶媒体として前記評価プログラムを記憶する。この結果、本実施形態のＲＯＭ１４ｂによれば、計算処理装置１４を、指標抽出手段、表示手段、及び、固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係に基づいて、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒ，左頸動脈８０Ｌ，右頸動脈８０Ｒにおける振動の固有振動数に対応するスティフネス及び予想破裂血圧を求める手段として容易に機能させることができる。
【００８６】
なお、前記実施形態を次のように変更して構成することもできる。
○ 前記実施形態では、左右の血管として、左浅側頭動脈７０Ｌ、右浅側頭動脈７０Ｒ、左頸動脈８０Ｌ、右頸動脈８０Ｒ、及び心電図の各種検出信号を取得するようにしたが、下記のように変更してもよい。
【００８７】
左浅側頭動脈７０Ｌ、右浅側頭動脈７０Ｒ、及び心電図の各種検出信号を取得してもよい。この場合、左右の浅側頭動脈に対して、前記実施形態の（１）、（２）の作用効果を奏することができる。すなわち、この場合は、前記Ｓ５０では、左浅側頭動脈７０Ｌ，右浅側頭動脈７０Ｒの固有振動数が取得される。又、Ｓ６０では、左浅側頭動脈７０Ｌ，右浅側頭動脈７０ＲのスティフネスＥｔｈ及び予想破裂血圧が得られる。
【００８８】
或いは、左頸動脈８０Ｌ、右頸動脈８０Ｒ、及び心電図の各種検出信号を取得するようにしてもよい。この場合、左右の浅側頭動脈に対して、前記実施形態の（１）、（２）の作用効果を奏することができる。
【００８９】
この場合は、左頸動脈８０Ｌ、右頸動脈８０Ｒに対して、前記実施形態の（１）、（２）の作用効果を奏することができる。すなわち、この場合は、前記Ｓ５０では、左浅側左頸動脈８０Ｌ、右頸動脈８０Ｒの固有振動数が取得される。又、Ｓ６０では、左頸動脈８０Ｌ、右頸動脈８０ＲのスティフネスＥｔｈ及び予想破裂血圧が得られる。
【００９０】
○ 前記実施形態では、固有振動数に基づいて、スティフネスＥｔｈ及び予想破裂血圧を求めるようにしたが、スティフネスＥｔｈ及び予想破裂血圧のいずれか一方のみを求めるようにしてもよい。
【００９１】
○ 前記実施形態では、マイク部１１ｍを収納するケース１１は、円筒状としたが、イヤホーン形状としてもよく、耳甲介腔５０Ｌ，５０Ｒに挿入できる形状であればよい。
○ 前記実施形態では、ＲＯＭ１４ｂを評価プログラムを記憶する記憶媒体としたが、記憶装置２６を評価プログラムを記憶する記憶媒体にしてもよい。或いは、ＣＤや、ＤＶＤ、ＵＳＢメモリ等の記憶媒体に評価プログラムを格納して、ＣＤやＤＶＤ等の記憶媒体を読込みするドライバ装置等にて読込みして評価プログラムを実行するようにしてもよい。
【００９２】
○ 前記実施形態において、マイク１１Ｌ，１１Ｒに代えて振動計としてもよい。
○ 前記実施形態において、振動センサ１２Ｌ，１２Ｒに代えてマイクとしてもよい。
○ 前記実施形態では、振動に関する物理量としては、固有振動数、位相、振幅としたが、振動の減衰率としてもよい、減衰率は、動脈硬化が進行すると、青年期の血管よりも老年期の血管の方が、弾性率が悪くなるため早く減衰する傾向がある。従って、左右血管の振動の減衰率を測定（抽出）し、差を算出（抽出）するようにしてもよい。
【００９３】
○ Ｓ６０において、データ化した結果、及び左右の血管の固有振動数、振幅Ａｎ等をディスプレイ２５に表示する際、これらの左右の血管の差をディスプレイ２５に表示するようにしてもよい。
【符号の説明】
【００９４】
１０…評価装置、１１Ｌ，１１Ｒ…マイク（振動検出手段）、
１２Ｌ，１２Ｒ…振動センサ（振動検出手段）、
１３…電極、１４…計算処理装置（指標抽出手段、表示手段）、
１４ａ…ＲＯＭ（記憶媒体）、２６…記憶装置。
【特許請求の範囲】
【請求項１】
測定対象血管が人の左右対称的な位置にそれぞれ位置する血管であって、血流による各測定対象血管の振動を同時に計測する第１ステップと、
計測して得られた前記振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する第２ステップと、
前記測定対象血管における固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係のうち、少なくともいずれか一方の相関関係に基づいて、第２ステップで抽出した固有振動数に対応するスティフネス及び予想破裂血圧のうち少なくともいずれか一方を求める第３ステップを含むことを特徴とする測定対象血管の力学的機能の評価方法。
【請求項２】
前記第２ステップでは、さらに、前記振動のデータに基づいて、前記各測定対象血管の振動をリサージュ波形で表示することを特徴とする請求項１に記載の測定対象血管の力学的機能の評価方法。
【請求項３】
前記測定対象血管が人の左右対称的な位置にそれぞれ位置する血管が、左右の浅側頭動脈であることを特徴とする請求項１又は請求項２に記載の測定対象血管の力学的機能の評価方法。
【請求項４】
測定対象血管が人の左右対称的な位置にそれぞれ位置する血管であって、血流による各測定対象血管の振動を同時に計測する振動検出手段と、
前記計測して得られた前記振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する指標抽出手段と、
前記測定対象血管における固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係のうち、少なくともいずれか一方の相関関係に基づいて、前記指標抽出手段が抽出した固有振動数に対応するスティフネス及び予想破裂血圧のうち少なくともいずれか一方を求める手段とを備えることを特徴とする測定対象血管の力学的機能評価装置。
【請求項５】
前記各測定対象血管の振動のデータに基づいて、リサージュ波形で表示する表示手段を備えることを特徴とする請求項４に記載の測定対象血管の力学的機能評価装置。
【請求項６】
コンピュータを、
人の左右対称的な位置にそれぞれ位置する各測定対象血管の血流による振動のデータであって、同時に計測して得られた前記振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する指標抽出手段と、
前記測定対象血管における固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係のうち、少なくともいずれか一方の相関関係に基づいて、前記指標抽出手段が抽出した固有振動数に対応するスティフネス及び予想破裂血圧のうち少なくともいずれか一方を求める手段として機能させることを特徴とする測定対象血管の力学的機能の評価プログラム。
【請求項７】
請求項６において、
コンピュータを、
前記各測定対象血管の振動のデータに基づいて、リサージュ波形で表示する表示手段として機能させることを特徴とする測定対象血管の力学的機能の評価プログラム。
【請求項８】
請求項６又は請求項７に記載の測定対象血管の力学的機能の評価プログラムを記憶した記憶媒体。
【請求項１】
測定対象血管が人の左右対称的な位置にそれぞれ位置する血管であって、血流による各測定対象血管の振動を同時に計測する第１ステップと、
計測して得られた前記振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する第２ステップと、
前記測定対象血管における固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係のうち、少なくともいずれか一方の相関関係に基づいて、第２ステップで抽出した固有振動数に対応するスティフネス及び予想破裂血圧のうち少なくともいずれか一方を求める第３ステップを含むことを特徴とする測定対象血管の力学的機能の評価方法。
【請求項２】
前記第２ステップでは、さらに、前記振動のデータに基づいて、前記各測定対象血管の振動をリサージュ波形で表示することを特徴とする請求項１に記載の測定対象血管の力学的機能の評価方法。
【請求項３】
前記測定対象血管が人の左右対称的な位置にそれぞれ位置する血管が、左右の浅側頭動脈であることを特徴とする請求項１又は請求項２に記載の測定対象血管の力学的機能の評価方法。
【請求項４】
測定対象血管が人の左右対称的な位置にそれぞれ位置する血管であって、血流による各測定対象血管の振動を同時に計測する振動検出手段と、
前記計測して得られた前記振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する指標抽出手段と、
前記測定対象血管における固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係のうち、少なくともいずれか一方の相関関係に基づいて、前記指標抽出手段が抽出した固有振動数に対応するスティフネス及び予想破裂血圧のうち少なくともいずれか一方を求める手段とを備えることを特徴とする測定対象血管の力学的機能評価装置。
【請求項５】
前記各測定対象血管の振動のデータに基づいて、リサージュ波形で表示する表示手段を備えることを特徴とする請求項４に記載の測定対象血管の力学的機能評価装置。
【請求項６】
コンピュータを、
人の左右対称的な位置にそれぞれ位置する各測定対象血管の血流による振動のデータであって、同時に計測して得られた前記振動のデータを周波数解析して、前記各測定対象血管の力学的機能の指標となる固有振動数を抽出する指標抽出手段と、
前記測定対象血管における固有振動数とスティフネスとの相関関係、及び、前記測定対象血管における固有振動数と予想破裂血圧との相関関係のうち、少なくともいずれか一方の相関関係に基づいて、前記指標抽出手段が抽出した固有振動数に対応するスティフネス及び予想破裂血圧のうち少なくともいずれか一方を求める手段として機能させることを特徴とする測定対象血管の力学的機能の評価プログラム。
【請求項７】
請求項６において、
コンピュータを、
前記各測定対象血管の振動のデータに基づいて、リサージュ波形で表示する表示手段として機能させることを特徴とする測定対象血管の力学的機能の評価プログラム。
【請求項８】
請求項６又は請求項７に記載の測定対象血管の力学的機能の評価プログラムを記憶した記憶媒体。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図１６】
【図１７】
【図９】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図１６】
【図１７】
【図９】
【公開番号】特開２０１０−１８７９２８（Ｐ２０１０−１８７９２８Ａ）
【公開日】平成２２年９月２日（２０１０．９．２）
【国際特許分類】
【出願番号】特願２００９−３５４８８（Ｐ２００９−３５４８８）
【出願日】平成２１年２月１８日（２００９．２．１８）
【国等の委託研究の成果に係る記載事項】（出願人による申告）平成１６〜２０年度、文部科学省、地域科学技術振興施策、委託研究（知的クラスター創成事業、岐阜・大垣地域ロボティック先端医療クラスター）、産業技術力強化法第１９条の適用を受ける特許出願
【出願人】（３０４０１９３９９）国立大学法人岐阜大学 (289)
【Ｆターム（参考）】
【公開日】平成２２年９月２日（２０１０．９．２）
【国際特許分類】
【出願日】平成２１年２月１８日（２００９．２．１８）
【国等の委託研究の成果に係る記載事項】（出願人による申告）平成１６〜２０年度、文部科学省、地域科学技術振興施策、委託研究（知的クラスター創成事業、岐阜・大垣地域ロボティック先端医療クラスター）、産業技術力強化法第１９条の適用を受ける特許出願
【出願人】（３０４０１９３９９）国立大学法人岐阜大学 (289)
【Ｆターム（参考）】
[ Back to top ]