JP2004358022A - 脈波解析法、この解析法に基づく自律神経機能評価法及び自律神経機能評価装置 - Google Patents
脈波解析法、この解析法に基づく自律神経機能評価法及び自律神経機能評価装置 Download PDFInfo
- Publication number
- JP2004358022A JP2004358022A JP2003161709A JP2003161709A JP2004358022A JP 2004358022 A JP2004358022 A JP 2004358022A JP 2003161709 A JP2003161709 A JP 2003161709A JP 2003161709 A JP2003161709 A JP 2003161709A JP 2004358022 A JP2004358022 A JP 2004358022A
- Authority
- JP
- Japan
- Prior art keywords
- pulse wave
- interval
- finger
- light
- standard
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
【課題】心電図に基づくのではなく、加速度脈波を用いて自律神経機能などを評価するための加速度脈波の解析法、この解析法を利用する自律神経機能評価方法及び自律神経機能評価装置の提供。
【解決手段】所定の時間連続して計測した脈波の波形を２次微分して算出した加速度脈波の連続する波形からａ−ａ間隔を求め、このａ−ａ間隔の変動を心電図のＲ−Ｒ間隔の変動に対応する間隔とする。この解析法に基づいて自律神経機能を評価する。
【選択図】 図３
【解決手段】所定の時間連続して計測した脈波の波形を２次微分して算出した加速度脈波の連続する波形からａ−ａ間隔を求め、このａ−ａ間隔の変動を心電図のＲ−Ｒ間隔の変動に対応する間隔とする。この解析法に基づいて自律神経機能を評価する。
【選択図】 図３
Description
【０００１】
【発明の属する技術分野】
本発明は、脈波解析法、この解析法を利用した自律神経機能評価法及び自律神経機能評価装置に関する。
【０００２】
【従来の技術】
従来から、心拍変動性による自律神経機能を定量的指標で評価するために、非侵襲的で簡便に実施できる方法として、もっぱら心電計による計測に基づく評価法が実施されている。この心電図記録による自律神経機能判定は、１００拍、１５０拍、２分間で判定されることが多い。この場合、心電図のＲ−Ｒ間隔の変動係数を求めて、年齢別の「正常値」を参考に臨床的に診断するのが一般的である。ただし、いわゆる「正常値」は普遍的な数値ではなく、文献上の数値を参考にした上で、各医療機関で独自に決めている。
心拍変動性は加齢の影響を受け、その変動係数は加齢によって減少する。これは、これまで報告された心拍変動係数をまとめて示す表１から明らかである。
【０００３】
（表１）
【０００４】
上記心拍変動係数の「正常値」に性差はないといわれている。また、昼間の心拍変動係数は低く、夜間の心拍変動係数は高いという日内変動を示すが、日中変動の幅は２％以内であって有意差はなく判定に影響はないことが報告されている（例えば、非特許文献１参照）。また、心拍変動性の低下は、自律神経機能の異常、なかでも副交感神経機能低下を反映するものと考えられている。この自律神経機能の異常とは、交感神経機能と副交感神経機能とのバランスの異常である。
【０００５】
例えば、糖尿病患者では、しびれ、アキレス腱反射の減弱又は消失、振動覚障害などの体性神経障害とともに、消化管機能異常、排尿障害、起立性低血圧などの自律神経障害が多くみられるが、このような自律神経症状のない患者でも、心電図Ｒ−Ｒ間隔の心拍変動係数は減少していることが多い（例えば、非特許文献２参照）。また、一部の脳梗塞患者や、神経疾患の患者でも心拍変動係数の低下がみられる（例えば、非特許文献３）。このような自律神経機能異常は循環動態の不安定化を招き、合併症のリスクを高める。
また、脈波による心拍変動の計測も実験的には行われており、例えば、検出した容積脈波を微分した速度脈波を解析して心電図のＲ−Ｒ間隔に相当する情報を取り出す脈波解析法が知られている（例えば、特許文献１参照）。しかし、脈波による心拍変動の計測は未だ医療の現場では実用化されていない。
【０００６】
脈圧とは収縮期血圧と拡張期血圧との差であるが、これは１回拍出量と動脈系のコンプライアンスとの関数である。脈圧の波形は大動脈から末梢の動脈に行くにつれて変容していく。これは種々の部位において投射波と反射波とが合成され、共鳴が起こるからと考えられ、変容の程度は血管の性状又は特性の影響の総和とみることができる。大動脈と末梢動脈とにおける動脈内圧を比較すると、末梢動脈における平均動脈圧及び拡張期血圧は大動脈の場合よりも低い。しかし、収縮期血圧は、投射波と反射波との合成により、末梢動脈の方が高くなるため、脈圧は大きくなるという現象が起こる。それ故、末梢脈波の波形ほど中心脈波よりも凹凸が大きく波形判別がし易いという利点がある。それでも原波形は基線が安定せず変曲点の認識が困難な場合がある。
【０００７】
そこで、指先で測定する末梢脈波である容積脈波の原波形を２回微分して得られる二次微分波である加速度脈波が、より評価に適した波形として研究や臨床で利用されている。現在、多く使われている脈波計は光電式指先容積脈波計である。この脈波計の原理は、指先にヘモグロビンに吸光特異性のある波長の光を当てて、吸収光、又は反射光から血管内の血流の容積変化を求め、波形を得る方法に基づいている。容積脈波と圧脈波とは反映しているものは異なるが、波形のもつ意味は同等である。
上記末梢脈波の加速度脈波が生体のどのような生理活動（自律神経機能など）を表現しているのかについては、平均波形の変化と自律神経機能との対応の研究については知られているが、個々の成分波の変化と自律神経機能との対応を検討した研究については知られていない。
【０００８】
加速度脈波は心臓の収縮期の波形であり、ａ〜ｅ波の５つの成分波をもつ（図１）。これらの成分波は、生体の条件により、また、加齢に応じて一定の変化をする。頂点ａを有する成分波（ａ波）は原波形の立ち上がり部分の波形に相当し、頂点ｅを有する成分波（ｅ波）は原波形の収縮終期のノッチ部分に相当する。それゆえ、期外収縮時や頻脈時にはａ−ｅ間隔は短縮し、各成分波も変形する。波形の基線より上を正、下を負の象限としたとき、ａ波は基線より常に上に位置する陽性波であり、ｂ波は基線より常に下に位置する陰性波であり、ｃ波、ｄ波、ｅ波は生体の条件により陽性又は陰性に変化する成分波である。
【０００９】
また、加速度脈波の波形は、加齢によって変化していくので、生理的変化と病的変化の判別も難しい。測定器械、測定環境、年齢などによって異なる波形を正しく評価するには、どのような条件にも適応できるような共通の波形評価法を確立しなければならない。すなわち、測定環境が同一であれば、どの器械にも、どの年齢層にも共通して使用できる共通の基準を作成し、それに基づいて評価することが必要であるが、いまだ有用な評価法は確立されていない。
【００１０】
【非特許文献１】
景山茂他：心電図Ｒ−Ｒ間隔変動、内科、５５： ２４２−２４６（１９８５）
【非特許文献２】
Ｂｅｒｎａｒｄｉ Ｌ．， Ｒｉｃｏｒｄｉ Ｌ．， Ｌａｚｚａｒｉ Ｐ． ｅｔ ａｌ：Ｉｍｐａｉｒｅｄ ｃｉｒｃａｄｉａｎ ｍｏｄｕｌａｔｉｏｎ ｏｆ ｓｉｍｐａｔｈｅｔｉｃ ａｃｔｉｖｉｔｙ ｉｎ ｄｉａｂｅｔｅｓ， Ｃｉｒｃｕｌａｔｉｏｎ， １９９２； ８６： １４４３−１４５２
【非特許文献３】
宮田正和、松浦達雄：ＣＶＲ−Ｒ、心身医療、Ｖｏｌ． ３， Ｎｏ． ８ １９９１
【特許文献１】
特開２００１−７０２６５号公報（特許請求の範囲）
【００１１】
【発明が解決しようとする課題】
自律神経活動評価としての心拍変動測定の基礎は、心臓の拍動の原因である心電パルスの発生間隔の変動である。特に、心室収縮の開始を表すＲ波の間隔（Ｒ−Ｒ間隔）を評価する。しかし、心電計による計測は簡便ではない。特に、体温計や血圧計による計測の場合と異なり、患者が家庭で測定して経過を見るということはできない。従って、上記のような心電図のＲ−Ｒ間隔の変動係数を求めて、自律神経機能を評価する簡便な方法の開発が求められている。
また、心電図のＲ−Ｒ間隔の変動係数の代わりに、容積脈波や速度脈波を利用して自律神経機能を評価しようとしても、基線変動や体動ノイズなどによる脈の誤診や検出ミスがあり、特に脈が小さい場合、判定が困難であるという問題もある。
【００１２】
そこで、本発明の課題は、心電図に基づくのではなく、加速度脈波を用いて自律神経機能などを評価するための加速度脈波の解析法、この解析法を利用する自律神経機能評価方法及び自律神経機能評価装置を提供することにある。
【００１３】
【課題を解決するための手段】
本発明の脈波解析法は、生体の脈波から心電図のＲ−Ｒ間隔の変動に対応する情報を求める脈波解析法であって、所定の時間連続して計測した脈波の波形を２次微分して加速度脈波を算出し、この加速度脈波の連続する波形からａ−ａ間隔を求め、このａ−ａ間隔の変動を心電図のＲ−Ｒ間隔の変動に対応する間隔とすることを特徴とする。
【００１４】
本発明の脈波解析法はまた、生体の脈波から心電図のＲ−Ｒ間隔の変動に対応する情報を求める脈波解析法であって、過去の加速度脈波、好ましくは過去数拍から数十拍の加速度脈波から求めた基準（標準）となる波形を標準加速度脈波とし、この標準加速度脈波と測定加速度脈波との波形の類似度を評価して個別加速度脈波を決定し、この個別加速度脈波の連続する波形からａ−ａ間隔を求め、このａ−ａ間隔の変動を心電図のＲ−Ｒ間隔の変動に対応する間隔とすることを特徴とする。
上記波形の類似度の評価を、測定加速度脈波のａ波候補の波高値が標準加速度脈波のａ波の波高値と同一若しくは両者の中間値になるように規格化した波形を用いて行い、また、標準加速度脈波のａ波と測定加速度脈波のａ波候補との時間位置を合わせ、規格化した測定加速度脈波と標準加速度脈波との波高差を求めて行うことを特徴とする。
【００１５】
上記波形の類似度の評価を、個別加速度脈波のａ波判定条件として、測定加速度脈波をそのａ波候補の波高値と標準加速度脈波のａ波波高値との中間値若しくはａ波波高値で規格化した個別加速度脈波と標準加速度脈波との波高差の積分値又は多重積分値を用いて算出した類似度に基づいて行うことを特徴とする。
上記脈波解析法において、隣り合うａ−ａ間隔の変化が一定範囲を外れるものを除外したａ−ａ間隔について、これらの値により標準ａ−ａ間隔を構成して、ａ−ａ間隔の異常を判定することを特徴とする。
【００１６】
上記標準ａ−ａ間隔として、過去のデータの平均若しくは重み付き平均を用い、異常判定がでたデータに対して、未来データを含む平均若しくは重み付き平均で再判定することを特徴とする。以下、「平均若しくは重み付き平均」を簡略的に「（重み付き）平均」と言うこともある。
上記脈波解析法において、評価している時点のａ−ａ間隔が、過去の平均的なａ−ａ間隔、又は未来よりも過去のａ−ａ間隔データの重みを大きくした重み付き推移平均、又は直前のａ−ａ間隔を基準として、一定の基準値を超えるとき、その時点の直後から順次、隣り合うａ−ａ間隔の変化が一定範囲内に入るものを数点選び出し、これらの値により標準ａ−ａ間隔を更新して、ａ−ａ間隔の異常を判定することを特徴とする。
【００１７】
上記脈波解析法において、生体の脈波を、発光部と受光部とを有する反射型脈波センサであって、発光部が、受光部より指の動脈血流の下流側に、その上面が受光部の上面よりも突出しかつ所定の距離だけ指の腹部を載置する床面よりも突出するように配置され、そして指の先端部を発光部よりも指動脈血流のさらに下流側に装着するための空間が床面の先端部分に設けられている脈波センサに指先を装着して測定することを特徴とする。
脈波の測定は、上記脈波センサの発光部がその上面が指の腹部を載置する床面より０．１〜１．５ｍｍ突出するように、発光部が配置されている脈波センサを用いて行われることを特徴とする。
【００１８】
上記脈波の測定は、脈波センサの受光部の上面が発光部の上面と同一レベルか若しくはそれより下になるように、また、所定の距離だけ床面より低くなるように、受光部が配置されている脈波センサを用いて行われることを特徴とする。
上記脈波の測定は、脈波センサの指の先端部を装着する空間の床面と対向する面に押圧材が取付けられ、この押圧材で発光部よりも指動脈血流のさらに下流側の指の先端部の上面を押さえるように構成されている脈波センサを用いて行われることを特徴とする。
上記押圧材は、この押圧材で発光部よりも指動脈血流のさらに下流側の指の先端部の上面を５０〜２００ｇＷの圧力で押さえるように構成されていることを特徴とする。
【００１９】
上記脈波の測定は、脈波センサの発光部の側面が、内面が照射光に対する反射特性を有する遮光壁で囲まれ、そして遮光壁上端が指の腹部を載置する床面から０．１〜１．５ｍｍ突出するように構成されている脈波センサを用いて行われることを特徴とする。
上記脈波の測定は、脈波センサの発光部及び受光部の上方に、赤外線透過性の窓部を、窓部上面が指の腹部を載置する床面よりも所定の距離だけ突出するように配置してある上記脈波センサを用いて行われることを特徴とする。
上記窓部の上面が指の腹部を載置する床面より０．１〜０．５ｍｍ突出するように配置されていることを特徴とする。
【００２０】
また、本発明の自律神経機能評価法は、上記脈波解析法に基づき、連続するａ−ａ間隔について、隣り合うａ−ａ間隔の最初のａ−ａ間隔が標準ａ−ａ間隔より正常範囲を外れて短い場合で、かつ、前回のａ−ａ間隔と今回のａ−ａ間隔との和が標準ａ−ａ間隔の１倍若しくは２倍を基準にしてａ−ａ間隔の変動の分布から定まる所定の範囲に入るとき、期外収縮と判定することを特徴とする。
この自律神経機能評価法はまた、ａ−ａ間隔の統計量に基づいて自律神経機能を評価する方法であって、上記脈波解析法に基づき、ａ−ａ間隔が正常範囲を外れる場合はこれを除外してａ−ａ間隔の統計量を算出し、これに基づいて自律神経機能を評価することを特徴とする。
【００２１】
この自律神経機能評価法はまた、ａ−ａ間隔の時間変動を周波数解析して自律神経機能を評価する方法であって、隣り合うａ−ａ間隔について、最初のａ−ａ間隔が、上記脈波解析法に基づき、標準ａ−ａ間隔より正常範囲を外れて短い場合で、かつ、最初のａ−ａ間隔と次のａ−ａ間隔との和が、標準ａ−ａ間隔の１倍若しくは２倍を基準にして所定の範囲（例えば、０．７５〜１．２５若しくは１．５〜２．５倍）に入る時、１倍を基準とする場合には、最初のａ−ａ間隔を最初のａ−ａ間隔と次のａ−ａ間隔との両者の和に補正して次のａ−ａ間隔を削除し、また、２倍を基準とする場合には、それぞれ二つのａ−ａ間隔を両者の和の二分の一に補正して、周波数解析することにより自律神経機能を評価することを特徴とする。
【００２２】
上記自律神経機能評価法において、ａ−ａ間隔の変動性を各年代別標準参考値を参照しながら評価することにより、自律神経機能を評価することを特徴とする。
上記自律神経機能評価法において、上記脈波解析法に基づき、長時間連続測定や期外収縮波の検出を行うことにより、また、期外収縮除去後の推定変動係数を算出することにより、糖尿病、神経障害、脳血管疾患、冠動脈疾患、喘息、更年期障害などを含めて、合併症リスクの予測や、治療効果判定、自己健康管理を行うことを特徴とする。
【００２３】
さらに、本発明の加速度脈波の変動解析による自律神経機能評価装置は、生体の脈波を検出して脈波の大きさに応じた信号を出力する脈波計測手段と、脈波計測手段により得られた脈波の波形を２次微分して算出された加速度脈波の波形からその波形パラメータを解析する波形パラメータ解析手段とを備え、この波形パラメータ解析手段が、所定の時間連続して計測した脈波の波形を２次微分して加速度脈波を算出する手段と、この加速度脈波の連続する波形から心電図のＲ−Ｒ間隔の変動に対応するａ−ａ間隔の変動を求める手段とを有することを特徴とする。
上記自律神経機能評価装置は、さらに、隣り合うａ−ａ間隔の変化が一定範囲を外れるものを除外したａ−ａ間隔について、これらの値により標準ａ−ａ間隔を構成して、ａ−ａ間隔の異常を判定する手段を有することを特徴とする。
【００２４】
【発明の実施の形態】
本発明によれば、加速度脈波を変動解析し、この変動解析に基づいて自律神経機能などを評価し、各種医療診断や健康診断を行うことができる。この変動解析は、生体の指先に取り付けて使用される脈波センサの計測・検出結果に基づき行われる。この脈波センサからなる脈波測定装置は、所定の加速度脈波情報からａ−ａ間隔を抽出するための脈波解析プログラムが組み込まれたマイクロコンピュータを有している。
【００２５】
図２に脈波センサ（（株）ユメディカ社製のＡＲＴＥＴＴ（登録商標））を用いて測定・記録した容積脈波、速度脈波（一次微分波）、加速度脈波（二次微分波）を示す。この加速度脈波は、図１に示すようにａ波、ｂ波、ｃ波、ｄ波、ｅ波の５つの成分波をもつ。これら５つの成分波で構成される加速度脈波は心臓収縮期の波であり、これらの成分波は、生体の状態により、また、加齢に応じて一定の変化をしていく。ただし、ｅ波は収縮と拡張の境界波である。なお、この脈波センサを用いて測定することにより、体動などによるノイズの発生を極力押さえた連続波形データを得ることができる。
【００２６】
図３に、心電図（ＥＣＧ）と加速度脈波（ＡＰＧ）とを同時測定した波形を示す。図３から明らかなように、心電図Ｒ波の出現に一定の時間的な遅れ（α）を伴って、加速度脈波のａ波が出現している。このとき、加速度脈波のａ−ａ間隔（ＡＰＧ／ａａ間隔（Ｔａａ））は、心電図のＲ−Ｒ間隔（ＥＣＧ／ＲＲ間隔）にほぼ一致する（図４）。
【００２７】
［加速度脈波ａ−ａ間隔の計測と変動係数算出］
本発明者らの知見によれば、血管の老化の程度を評価するための加速度脈波波形評価には、１０秒から２０秒程度までの測定期間において容積脈波が安定している期間の１０拍程度が検出できれば充分であった。しかしながら、自律神経機能評価のためには、１００拍以上、場合によっては１５分以上の連続測定が必要であって、脈拍間隔の１ｍｓｅｃ以内の精度での連続測定が要求される。従って、一拍毎の個別加速度脈波の精度の高い検出と、ａ波の時間位置検出精度が問題となる。
【００２８】
上記したように、従来の自律神経活動評価における心拍変動測定の基礎は、心臓の拍動の原因である心電パルスの発生間隔の変動、特に、心室収縮の開始を表すＲ波の間隔（Ｒ−Ｒ間隔）の変動を測定し、評価することにある。このＲ−Ｒ間隔を、他の心拍検出手段で代用することも様々に行われている。例えば、心音や、頸動脈圧脈波、指尖容積脈波などで心拍数若しくは脈拍数を測定することが行われている。従来の方法では、一般的には、これらの拍動を電気信号に変換したもの（一般的には電圧）が一定レベルを超えるときにパルスを発生させて、１分間の拍数をカウントしているに過ぎない。
しかるに、数％程度の変動を計測しなければならない自律神経機能の評価のためには、ミリ秒程度の精度で個々の拍動の間隔を計測する必要があるので、パルスの発生位置が問題となる。
【００２９】
容積脈波は、基線が変動し、安定した検出が困難である。また、時間位置精度もでない。気温が低い場合や体調の影響により、脈波の振幅が非常に小さい場合があり、この場合、Ｓ／Ｎ（信号／雑音）比が小さくなり、脈の検出そのものができない場合がある。速度脈波や加速度脈波の検出も試みられているが、Ｓ／Ｎ比が悪い場合には、ノイズのピークが多く発生し、検出すべきピークの同定が困難である。
【００３０】
そして、この計測された心拍間隔が心電図のＲ−Ｒ間隔と一致していれば問題はないが、このことは必ずしも自明ではない。例えば、指尖容積脈波の極小点や極大点を拍動の時間位置として検出することが考えられるが、これらは、基線変動の影響を受けるし、個々の心拍の駆出力などの生理的な変動の影響も受けやすい。そして、検出する心拍変動係数が小さく、心電図のＲ−Ｒ間隔が１％以内となる自律神経障害を有効に検出するためには、少なくとも心電図のＲ−Ｒ間隔と数ミリ秒以内の精度での対応付けができている必要がある。しかし、従来の方法では、このような精度を達成することは困難である。
【００３１】
また、実用的な応用において、パーソナルコンピュータによる自動解析の出力結果が信頼できるものであるためには、異常なａ−ａ間隔を検出したときに、誤検出（ノイズ若しくはａ波以外のピークをａ波として検出した場合）、又は検出ミス（実際に存在するａ波を検出し損なった場合）によってａ−ａ間隔を間違って計測したのか、あるいは実際のａ−ａ間隔を検出しているのかが問題となる。
平均波形を算出して波形の特徴量だけを求める場合には、確実に個々の脈波を検出して検出できたものだけで平均を取ればよいが、連続するａ−ａ間隔の変動を検出し、ａ−ａ間隔の異常、正常の判定をするためには、個々の脈波の検出精度を飛躍的に高める必要がある。実際上は、測定波形とａ波検出位置とを同時に表示して確認するのが最善であり、如何に解析ソフトの解析精度が向上しようとも、視覚的な確認の必要性は残ると考えられるが、誤検出や検出ミスはできうる限り排除できるように構成すべきである。
【００３２】
以下、個別加速度脈波の検出・解析について、図５のフロー図を参照して説明する。
１．標準加速度脈波パラメータの初期値算出
測定開始時に、類似度判定のための標準加速度脈波とそのパラメータ及びａ−ａ間隔異常判定のためのａ−ａ間隔変動幅の初期値を求める必要がある。このため、波形をモニターしながら容積脈波の基線の変動及びノイズが少ないことを確認して測定を始め、容積脈波のピークやボトム若しくは立ち上がりを検出して脈波を切り出して、その時間間隔からａ−ａ間隔とその変動幅を求めると共に、平均波形を求め、これから標準加速度脈波とａ波波高値の初期値を求める。または、速度脈波や加速度脈波が安定していることを確認して、同様の処理をする。
【００３３】
以後、上で算出した標準値を用いて、測定開始時点から、改めて次の解析を行う。
２．ａ波候補の抽出
（１）測定加速度脈波の極大点を求め、求めた極大値とａ波の標準波高値とを比較して、一定範囲（例えば、６０〜１８０％）に入るものをａ波候補とする。６０％未満のものはノイズ若しくはｅ波を検出しているかもしれず、１８０％を超えるものはノイズの可能性が高いという問題がある。
【００３４】
３．個別加速度脈波のａ波判定
上記したように、個別加速度脈波のａ波判定条件として、ａ波の波高値で規格化した個別加速度脈波と標準加速度脈波との波高差の積分値若しくは多重（二重）積分値を用いて類似度を算出し、判定している。
（１）測定加速度脈波の波形を、標準加速度脈波のａ波の波高値若しくはこのａ波の波高値とａ波候補の波高値との中間値で規格化する。この規格化は、測定加速度脈波の振幅に［（標準加速度脈波のａ波の波高値）＋ｑ＊（ａ波候補の波高値）］／［（１＋ｑ）＊（ａ波候補の波高値）］（ｑは０以上の実数値であり、実用的には１程度である）を乗じて行う。
【００３５】
（２）測定加速度脈波のａ波候補点と標準加速度脈波のａ波との時間（ｔ）位置を合わせ（Ｔａとする）、ａ波の波高値で規格化した測定加速度脈波と標準加速度脈波との波高差を求める。
（３）Ｔａ−ｔ１〜Ｔａ＋ｔ２の範囲で、求めた波高差の絶対値を積分して類似度とする。ここで、ｔ１は、例えば０．１秒とし、ｔ２は、例えばａ−ｅ間隔程度の時間とすればよい。この範囲を広く取りすぎると、期外収縮などの不整脈で、ｅ波直後に次ぎのａ波が出現した場合など、拡張期の標準波形との差が大きく出て、類似度が小さくなり、脈が検出できなくなる。
この場合、Ｔａ−ｔ１〜Ｔａ＋ｔ２の範囲で、波高差を積分したものの絶対値を更に積分（二重積分）して類似度としてもよい。
【００３６】
（４）類似度が高いとき、ａ波候補をａ波として確定する。類似度が低いときは、ａ波候補をａ波と判定せず、ノイズと判定して、再度ａ波候補の抽出を行う。
（５）今回確定したａ波と前回のａ波との時間差を求めて、これをａ−ａ間隔Ｔａａとして算出する。
【００３７】
上記した加速度脈波の波高差の積分値若しくは多重積分値を用いて類似度を評価する根拠は以下の通りである。
（１）一般的に、加速度脈波は基線変動に強いので、波高差で評価しても、容積脈波、速度脈波の基線変動の影響を受けないというメリットがある。しかし、微分波形は一般に高周波ノイズの影響を受けやすく、波高差での評価ではスパイク状のノイズには対応できない。
【００３８】
（２）加速度脈波の積分は基本的には速度脈波なので、速度脈波の波高差を類似度とすることが考えられる。しかしながら、この場合、速度脈波の基線変動が大きいと、判定ミスが発生する。これは、速度脈波の波形がまったく同形でも、基線変化分だけの波高差が発生するからである。一方、速度脈波の微分波形である加速度脈波の値は、速度脈波の波形の勾配であり、積分値は速度脈波の波高値の変化を表し、基線変動の影響を受けない。そこで、速度脈波に基線変動があっても、変動の影響が除去される微分波形（加速度脈波）から積分値を求めて、これから類似度を算出して評価することにより、安定した個別脈波の判定が可能となる。
【００３９】
（３）速度脈波にスパイク状のノイズが無視できない場合には、上記第３（３）項に示した二重積分による類似度評価が有効である。類似度を求めるための積分の変形として、特に限られるわけではなく、上記した絶対値の積分以外に、波形データを二乗したものを積分し、これの平方根を算出してもよい。
（４）なお、速度脈波の波高差を用いて類似度を評価する場合の上記問題点（基線変動が大きい）を解決するため、速度脈波を脈拍周期で積分した値を速度脈波の基線として補正することが考えられるが、不整脈がある場合にはうまくいかない。
【００４０】
次に、上記した加速度脈波の波高差の積分（差積分）による類似度判定について、図６を参照して詳細に説明する。図６の横軸は時間（ｍｓｅｃ）であり、縦軸は波高値である。
測定加速度脈波Ａ１と標準加速度脈波Ａ０との差Ａ１−Ａ０のゼロ点Ｔ（ｎ）では、標準加速度に対応する速度脈波（標準速度脈波）Ｖ０の勾配と、測定加速度脈波に対応する速度脈波（測定速度脈波）Ｖ１の勾配とは一致しており、このときのＶ０に対するＶ１の波高差（Ｖ１−Ｖ０の値）をΔＶ（ｎ）とする（図の場合、この値は負になっている）。差Ａ１−Ａ０が正の期間では、測定加速度脈波は標準加速度脈波に対して相対的に上昇し、負の期間では下降する。例えば、正の期間（Ｔ（０）からＴ（１）までの期間）における差分Ａ１−Ａ０の積分値は、速度脈波の波高値の上昇ΔＶ（１）−ΔＶ（０）を表しており、負の期間（Ｔ（１）からＴ（２）までの期間）におけるＡ１−Ａ０の積分値は、速度脈波の波高値の下降ΔＶ（２）−ΔＶ（１）を表している。
【００４１】
一般的に、速度脈波は基線変動により全体のレベルが変化していることがある。図６に示した場合では、測定速度脈波Ｖ１の基線は低くなっており、このような場合、速度脈波の波高値で脈拍を検出すると検出ミスが発生することがある。標準速度脈波Ｖ０に対する測定速度脈波Ｖ１の波高差の変化は、このような基線変動の影響を受けない。従って、差分Ａ１−Ａ０の絶対値の積分は、速度脈波の基線変動の影響を受けずに、標準速度脈波Ｖ０に対する測定速度脈波Ｖ１の形状差の目安とすることができる。
【００４２】
しかしながら、波高差で評価する場合には、高周波ノイズやパルス状のノイズがあると、類似度の評価値が悪くなる。この問題は、差分Ａ１−Ａ０（絶対値ではない）の積分を行い、これの絶対値の積分から類似度を評価することにより解決できる。この場合は、基線変動の影響を除いて容積脈波の形状差を評価していることになる。
また、リアルタイム性を重視する場合、計算量削減のために上記第３（２）〜（４）項と第３（５）項とは順番を逆にしてもよい。つまり、ａ波候補と前回ａ波との時間間隔が正常範囲の時、類似度を評価せずにａ波と判定し、ａ−ａ間隔が異常なときのみ、類似度を評価することもできる。
【００４３】
４．ａ−ａ間隔の異常値判定法
測定中にａ−ａ間隔の異常が生じた場合、以下のようにしてその異常が判定される。
隣り合うａ−ａ間隔Ｔａａの変化ＤＴａａ、すなわち（ａ−ａ間隔変動）ＤＴａａ［ｎ］＝Ｔａａ［ｎ］−Ｔａａ［ｎ−１］が、一定範囲（例えば、ＤＴａａの標準偏差Ｓｄｖ＿ＤＴａａの２倍以内）を外れるＴａａを除外したＴａａについて、これらの値により標準ａ−ａ間隔を構成して、ａ−ａ間隔の異常を判定する。この標準ａ−ａ間隔は判定対象データ前後のデータの推移平均を用いることができるが、リアルタイム性考慮する場合、過去データの（重み付き）平均を用い、異常判定がでたデータに対して、未来データを含む（重み付き）平均で再判定すればよい。
【００４４】
リアルタイム性を考慮して、アルゴリズムを構成する場合は、評価している時点でのａ−ａ間隔が、過去の平均的なａ−ａ間隔、若しくは直前のａ−ａ間隔（標準ａ−ａ間隔）を基準として、一定の基準値（例えば、４＊Ｓｄｖ＿ＤＴａａ）を超えるとき、その時点の直後から順次、隣り合うａ−ａ間隔の変化ＤＴａａが一定範囲（例えば、２＊Ｓｄｖ＿ＤＴａａ）内に入るＴａａを数点選び出し、これらの値により標準ａ−ａ間隔を更新して、ａ−ａ間隔の異常を判定してもよい。
【００４５】
図７に基づいて、ａ−ａ間隔再判定処理の流れを説明する。
▲１▼ｔ＝０までのデータを用いて求めた標準ａ−ａ間隔をＳｔ＿Ｔａａ（補正前）とする。この標準ａ−ａ間隔の算出式は、例えば、後出（段落番号［００５４］）の重み付き平均による更新の際の算式による。
▲２▼ｔ＝１の時のＴａａ（Ｔａａ（１））をＳｔ＿Ｔａａ（補正前）と比較して異常と判定する。
▲３▼ＤＴａａ（２）＝Ｔａａ（２）−Ｔａａ（１）を評価して、Ｔａａ（２）を平均値算出データとする。この評価は、前述のように、ＤＴａａの標準偏差Ｓｄｖ＿ＤＴａａから定めた閾値よりも、評価しているＤＴａａが小さければ平均値算出データとし、大きければ平均値算出データに含めないものとする。
【００４６】
▲４▼同様にして、ＤＴａａ（３）、ＤＴａａ（４）、ＤＴａａ（５）を評価し、Ｔａａ（３）、Ｔａａ（４）は平均値算出データから除外し、Ｔａａ（５）を平均値算出データとする。
▲５▼Ｓｔ＿Ｔａａ（補正前）とＴａａ（２）、Ｔａａ（５）の平均値を求めてＳｔ＿Ｔａａ（補正後）とする。
▲６▼Ｔａａ（１）をＳｔ＿Ｔａａ（補正後）と比較し、▲２▼で一旦異常と判定した結果を、正常と再判定する。
【００４７】
図８は、実測データに本アルゴリズムを適用した場合の標準ａ−ａ間隔の推移を示す。図には、Ｔａａの変化（前回Ｔａａとの差）も示してあるが、これから分かるように、前回Ｔａａの差のみでは異常Ｔａａを判定することができない。
心拍間隔の変化は、カオス的であるともフラクタル的な変化パターンを示すとも言われており、非常に複雑であって、単一の判定基準で判定することが困難である。図８に示すようなステップ状の変化も正常な変化として解析データに含める必要がある。
一般によく用いられる単一の基準を適用する場合の問題点を次に説明する。
正常範囲を標準ａ−ａ間隔の±１５％（８５％〜１１５％）のような一定範囲とすることが考えられるが、ａ−ａ間隔の変動は個人差が大きいので、変動が大きい被験者の場合には、正常な変動を異常とする検出ミスが発生することが考えられ、逆に変動が小さい被験者の場合には、異常な変動を正常とする誤検出が発生することが考えられるので、好ましくない。
【００４８】
ａ−ａ間隔の平均値を基準にして、ａ−ａ間隔の分布全体から正常ａ−ａ間隔の範囲を定めた場合、上記と異なり、次のような問題点がある。ゆっくりした大きな変動において、例えばＴａａが長いときに変動範囲内の短いＴａａが発生すると、異常ａ−ａ間隔と判定すべきところが正常ａ−ａ間隔と判定されることになる（図９）。
また、判定対象になっているａ−ａ間隔測定値Ｔａａ［ｎ］の前後数拍の平均値を標準ａ−ａ間隔とすることの問題点、すなわち、推移平均から求めた標準ａ−ａ間隔からの差のみでは異常Ｔａａが判定できない理由は次の通りである（図１０）。
【００４９】
呼吸周期の変動幅が小さい場合に、大きなゆっくりした変化があるときで、呼吸周期の変動幅を超える階段状の変化があるとき、変化点前後の中間値が標準ａ−ａ間隔となる。このとき、この変化を正常と見なすのであれば、正常範囲の限界ＴｈｒＤＤＴａａはこの階段状の変化の半分以上にする必要がある。ところがこの場合、単発のａ−ａ間隔異常でこのＴｈｒＤＤＴａａ以内のものは、検出できないことになる。
また、図１０に示すように、必ずしも異常とはいえないのに異常と判定される場合がある。
【００５０】
図１１に、ａ−ａ間隔異常値判定に関し、上記で求めた標準ａ−ａ間隔と推移平均法で求めた標準ａ−ａ間隔とを、実測データについて比較したグラフを示す。この図から明らかなように、推移平均では、データ列に階段状のギャップがあるような場合に、正常データを異常データと判定するケースが発生するが、本法により判定ミス（誤判定）を回避することができる。
【００５１】
５．異常データの処理
（１）ａ−ａ間隔の統計解析（変動解析）においては、異常と判定されたａ−ａ間隔は全て除外して解析する。
（２）スペクトル解析においては、以下の判定条件により期外収縮と判定された場合は、補正して解析を行う。その他の異常データがある場合は、異常データが含まれないよう、解析区間を分割して解析すべきである。
【００５２】
隣り合うａ−ａ間隔の最初のａ−ａ間隔が、標準ａ−ａ間隔より正常範囲を外れて短い場合で、前回ａ−ａ間隔と今回ａ−ａ間隔との和が、標準ａ−ａ間隔の１倍若しくは２倍を基準にして所定の範囲（例えば、０．７５〜１．２５若しくは１．５〜２．５倍）に入るとき、期外収縮と定める。
上記所定の範囲は、望ましくは隣り合うａ−ａ間隔の差の変動の分布から定める。例えば、隣り合うａ−ａ間隔の差の標準偏差をＳｄｖ＿ＤＴａａとして、±ｎ＊Ｓｄｖ＿ＤＴａａ（ｎの値は前者の場合、４〜５、後者の場合８〜１０程度が望ましい）を所定の範囲とする。
【００５３】
６．標準波形及び標準パラメータの更新
前回測定のａ波Ｔａ［ｎ−１］と今回測定のａ波Ｔａ［ｎ］との差Ｔａａ［ｎ］＝Ｔａ［ｎ］−Ｔａ［ｎ−１］が、上記したように定めた正常範囲に入るとき、標準ａ−ａ間隔（Ｓｔ＿Ｔａａ）、ａ−ａ間隔の変動幅（Ｓｄｖ＿ＤＴａａ）、標準ａ波波高値（Ｓｔ＿ＡＣＣａ）及び標準加速度脈波波形を更新する。
【００５４】
この更新は、例えば、平均値若しくは重み付き平均値による。重み付き平均による更新は次のような算式よってもよい。
Ｓｔ＿Ｔａａ（今回）＝（（ｋ−１）＊Ｓｔ＿Ｔａａ（前回）＋Ｔａａ（正常測定値））／ｋこの場合、ｋ＝１のときは、今回測定のａ−ａ間隔を次回の標準ａ−ａ間隔とすることに相当する。
【００５５】
標準加速度脈波については、前回ａ波Ｔａ［ｎ−１］を基準とする前回加速度脈波で、標準加速度脈波の波形を更新する。正常範囲に入らないときは更新しない。
上記のようにして個別加速度脈波を検出・解析すれば、期外収縮波やその他の異常値を除去して推定変動係数を算出できるので、期外収縮が存在する事実と変動係数評価とを別にでき、臨床上有用である。
【００５６】
［加速度脈波のａ−ａ間隔の変動係数の標準参考値］
次に、加速度脈波のａ−ａ間隔の変動係数の標準参考値、糖尿病患者のａ−ａ間隔の変動係数、加速度脈波による自律神経機能評価の価値について説明する。
図４に示すように、心電図のＲ−Ｒ間隔と加速度脈波のａ−ａ間隔とは一致している。通常、心電図は仰臥位で測定するのに対して、加速度脈波は座位で測定する。仰臥位と座位のふたつの姿勢においては、自律神経バランスの状態が異なるため、同一人であっても、それぞれの変動係数が異なることが考えられる。そこで、健康人１２１人について、加速度脈波を２分間記録して、ａ−ａ間隔の変動係数（ＣＶａａ％）の平均値と標準偏差とを年令別に算定した。年代別の変動係数を図１２に示す。また、各年代における標準的な変動係数の範囲を、［平均値−標準偏差］〜［平均値＋標準偏差］として図１３に示す。
これらの図１２及び１３から明らかなように、加速度脈波のａ−ａ間隔の変動係数は、加齢と共に減少するが、その標準偏差はほぼ同じ範囲内に入ることが分かる。
【００５７】
［糖尿病患者のａ−ａ間隔の変動係数］
上記加速度脈波のａ−ａ間隔の変動係数について、糖尿病患者と健康人とを比較した。糖尿病患者２６人の加速度脈波を２分間記録して、上記したようにしてａ−ａ間隔の変動係数を算出した。図１４に健康人（Ｎｏｒｍａｌ）と糖尿病患者（ＤＭ）との変動係数を比較して示す。図１４から明らかなように、４０代及び５０代において糖尿病患者の変動係数の低下が見られた。
【００５８】
［加速度脈波による自律神経機能評価法の価値］
上記したように、加速度脈波ａ−ａ間隔の変動と心電図Ｒ−Ｒ間隔の変動とは対応する。従って、加速度脈波のａ−ａ間隔の変動性の評価は心電図による心拍変動性評価と同等のものであるといえる。加速度脈波のａ−ａ間隔の変動係数は、心電図のＲ−Ｒ間隔の変動係数よりも、年齢により０．２〜１．５程度高くでる可能性があるが、これは、座位測定のために仰臥位よりは呼吸の影響が大きいためであると考えられる。そこで、ａ−ａ間隔の変動係数は、その年代別標準参考値を参照しながら心拍変動性を判定することにより、自律神経機能評価に用いることができる。
【００５９】
上記システムでは、長時間連続測定や期外収縮波の検出も可能である。また、このシステムは、期外収縮除去後の推定変動係数を算出できるので、臨床上応用可能の範囲が広がるという利点がある。さらに、加速度脈波測定は着衣の着脱もなく、座位にて指先の脈波を測定するだけであるから、心電図の場合よりもより簡便である。さらにまた、加速度脈波測定装置は心電計に比べて安価である。従って、糖尿病、神経障害、脳血管疾患、冠動脈疾患、喘息、更年期障害など多方面において、合併症リスクの予測や、治療効果判定、自己管理などに幅広く用いることができる。
【００６０】
上記加速度脈波を測定するための脈波センサについて以下説明する。
本発明で用いることができる脈波センサは、発光部と受光部とを有し、発光部と受光部とにより被験者の指の脈波を測定する反射型の脈波センサであって、波形がかわってもａ波の検出ができるものであれば特に制限はなく、従来の脈波計測装置を用いることができる。例えば、本発明者らが開発した脈波センサであって、発光部を、受光部より指の動脈血流の下流側に、その上面が受光部の上面よりも突出しかつ所定の距離だけ指の腹部を載置する床面よりも突出するように配置し、そして指の先端部を発光部よりも指動脈血流のさらに下流側に装着するための空間が床面の先端部分に設けられているようなものであってもよい。この場合、このように構成することにより、指の発光部に対する密着性が良好になると共に、指動脈血流の下流側が圧迫されても、上流側が圧迫されるのに比べて脈波への影響は少なく、再現性よく脈波情報を得ることができる。また、指と発光部との接触面積は発光部の上面の面積と等しくなり、指にかかる接触面積が小さいので、波形の変化はほとんどない。
【００６１】
上記発光部は、その上面が指の腹部を載置する床面より通常０．１〜１．５ｍｍ程度、好ましくは０．２〜１．０ｍｍ程度、さらに好ましくは０．３〜０．５ｍｍ程度突出するように配置されてもよい。発光部をこのように配置すると、指腹部の皮膚面が発光部上面から下に被さるので、測定データに対する外乱光や漏洩光や反射光の影響を小さくすることができると共に、被験者が指を装着する際に、その突出部に触れることによりセンサ位置を認識して、指を所定の位置に載置し易いという利点がある。しかし、０．１ｍｍ未満であると、センサ位置を確認し難いので、指先を所定位置に載置し難いと共に、測定データに対する反射光の影響が大きくなる。また、１．５ｍｍを超えると、指皮膚面が床面から浮くため、不安定な装着状態になると共に、指を載置した時の指に対する圧力により波形の変形が生じて再現性が悪くなり、測定される脈波データにバラツキが生じて、正確な脈波情報を得ることが難しくなる。
【００６２】
上記受光部は、その上面が指の腹部を載置する床面と同じレベルになるように配置されているか、又は所定の距離だけ床面より低くなるように配置されていることが好ましく、これにより発光部に対する指の密着性がより良好になる。
上記指の先端部を装着する空間の床面と対向する面に押圧材を取付け、押圧材で発光部よりも指動脈血流のさらに下流側の指の先端部の上面を押さえるように構成することが好ましい。被験者は、脈波測定中に意識的・無意識的に指先に力を入れる場合がある。この場合に、力を抜くように指導すると、被験者の指の形状によっては、センサへの密着性が悪くなることがある。力を入れる場合も、力を抜く場合も、指の小さな動きによるノイズが発生する。押圧材を設けることにより、ノイズが減少し、測定データの再現性が増し、正確な脈波情報を得ることができる。
【００６３】
上記したように、押圧部位を指の先端部上面の小面積に限るように構成してあるので、押圧部位の上流側の指動脈血流部分の脈波を再現性よく正確に測定することができる。
上記発光部の側面は、発光部から指内に照射される光が外部に漏れるのを阻止しかつ指の腹部表面からの反射光を阻止するために、遮光壁で囲まれていることが好ましい。
【００６４】
上記発光部は、内面が照射光に対する反射特性を有する遮光壁内部に配置され、遮光壁の上端が、指の腹部を載置する床面から通常０．１〜１．５ｍｍ程度、好ましくは０．２〜１．０ｍｍ程度、さらに好ましくは０．３〜０．５ｍｍ程度突出するように構成され、指の腹部がこの上端に載置されて遮光壁上端全面を覆うようにすることが好ましい。遮光壁の上端をこのような範囲内になるように突出させると、指腹部の皮膚面が発光部上面から下に被さるので、測定データに対する外乱光や漏洩光や反射光の影響を小さくすることができると共に、被験者が指を装着する際に、遮光壁の突出位置に触れることによりセンサ位置を認識して、指を所定の位置に載置し易いという利点がある。
【００６５】
しかし、０．１ｍｍ未満であると、遮光壁位置を確認し難いので、指先を所定位置に載置し難いと共に、発光部からの照射光や指腹部表面からの反射光が漏れ易くなり、測定データに対する反射光の影響が大きくなる。また、１．５ｍｍを超えると、指皮膚面が床面から浮くため、不安定な装着状態になると共に、指を載置した時の指に対する圧力により波形の変形が生じて再現性が悪くなり、測定される脈波データにバラツキが生じて、正確な脈波情報を得ることが難しくなる。
上記押圧材で指の先端部の上面を押さえる際の圧力を通常５０〜２００ｇＷ、望ましくは７０〜１５０ｇＷになるように構成することが好ましい。
【００６６】
また、本発明で用いることができる脈波センサは、発光部と受光部とを有し、発光部と受光部とにより被験者の指の脈波を測定する反射型の脈波センサにおいて、発光部を、受光部より指の動脈血流の下流側に、その上面が受光部の上面よりも突出するように配置し、また、発光部及び受光部の上方に、赤外線透過性の窓部を、窓部上面が指の腹部を載置する床面よりも所定の距離だけ突出するように配置し、そして指の先端部を発光部よりも指動脈血流のさらに下流側に装着するための空間が床面の先端部分に設けられているものであっても良い。
【００６７】
上記窓部の上面は、指の腹部を載置する床面より通常０．１〜１．５ｍｍ程度、望ましくは０．２〜１．０ｍｍ程度、更に望ましくは０．３〜０．５ｍｍ程度突出するように配置されることが好ましい。
上記押圧材で窓部より先の指の先端部の上面を押さえる際の圧力を通常５０〜２００ｇＷ、望ましくは７０〜１５０ｇＷになるように構成することが好ましい。５０ｇＷ未満であると、脈波の振幅が小さく、測定時の振動や被験者の指のふるえによるノイズが入りやすく、また、被験者の緊張などにより測定波形が変形しやすい。また、２００ｇＷを超えると測定波形が極端に変形する。なお、指全体を加圧すると測定波形が変形し、再現性が低下する。
【００６８】
以下、図面を参照して本発明で用いることができる脈波センサについて図面を参照して詳細に説明する。
図１５は、反射型脈波センサの要部である指装着部の断面構造を示す図（図１５（ａ））、並びに発光部及び受光部の近傍を拡大して示す図（図１５（ｂ））であり、指を装着した状態で示す。
【００６９】
この反射型脈波センサは、発光部及び受光部を有し、被験者の指を装着して脈波の測定ができるものであって、蓋部を構成する開閉自在の合成樹脂製の上方部分と、指腹部を載置することができるように構成された合成樹脂製の床部分とを有してなる。上方部分はその内面が指の外形に沿った形状をしていてもよく、また、床部分はその床面が平坦であっても又は外乱光を遮断するために指の付け根側がやや高く又は低くなるように傾斜した形状になっていてもよい。以下述べるように、上方部分の先端部には押圧材が設けられ、床面に載置した指の動脈血流の発光部よりもさらに下流側の指の先端部分を押さえて固定できるようになっており、また、床部分には発光部及び受光部が所定の位置に配置されている。この押圧材は、指の先端部分を押さえて固定できるものであればよく、例えば、クッション材や、バネ材のような板材などでよく、所定の圧力で指の先端部分を押さえることができるように構成されている。また、このセンサには、反射光の電流／電圧変換回路、増幅器が設けられており、このセンサをパーソナルコンピュータに接続すれば、センサからの出力に基づき正確な脈波情報を得ることができる。
【００７０】
この脈波センサの場合、指装着部に指を差し込み、指の先端部分の腹部に赤外線などの光を当てると、毛細血管中のヘモグロビン（赤血球）が光の一部を吸収し、光の反射量が変化する（血液量が多い部分は光の反射量が少なくなる）。この微妙な光の反射量の変化を検出し、検出された反射光を電流から電圧へ変換し、増幅器に送信し、増幅された信号電圧をパーソナルコンピュータを利用してＡＤ変換して出力し、脈波情報として活用する。
【００７１】
脈波センサの要部を構成する指装着部には、図１５（ａ）及び（ｂ）に示すように、発光ダイオード（ＬＥＤ）などの半導体発光素子からなる発光部１が、フォトダイオード（ＰＤ）などの半導体受光素子からなる受光部２よりも被験者の指３の動脈血流の下流側に配置される。指内の発光部１からの照射光１ａの経路を見ると、発光部分の光束は、指内を進むにつれて拡散して拡がる。このため、発光部１からの入射光量変化による受光部２の光量変化は大きく、受光部２の位置変化による受光する拡散光の光量変化は小さい。従って、発光部１を指に密着するようにすることが必要となる。しかし、密着性を良くすることは、指に余分の圧力を加えることにつながる。そこで、本発明では、発光部１を受光部２より指動脈血流の下流側に配置し、指に余分の圧力が加わらないようにする。
【００７２】
また、発光部１は、その上面が受光部２の上面より所定の距離だけ突出する（すなわち、高くなる）ように配置される。すなわち、発光部１の高さＨ１が受光部２の高さＨ２より所定の距離だけ高くなるように構成する。指装着部の先端部分には、赤外ＬＥＤなどからなる発光部１から照射される光１ａの経路よりも指動脈血流のさらに下流側に空間４が設けられ、指３の先端部がこの空間内に載置され得るように構成される。
【００７３】
指装着部の指腹部を載置する面は、指載置床面５として構成される。床面５には発光部１及び受光部２が所定の位置に設けられ、床面の先端部分は傾斜して立ち上がり、指の先端が適切に納まるように構成される。この指装着部において、発光部１の配置された位置より動脈血流の下流側であって、床面と対向する面に押圧材６が設けられる。この押圧材により、装着された指の先端部分（好ましくは、爪部）を軽く押さえ、装着された指が動かないように固定する。このように構成することにより、被験者の意識的・無意識的な指の小さな動きが少なくなって、ノイズ発生が減少し、その結果、測定波形の変化も少なくなる。なお、押圧材により動脈血流の下流側を圧迫しても、脈波への影響は小さい。
【００７４】
発光部１からの照射光１ａが指の皮膚表面で反射して受光部２に入り込むと、この反射光がノイズとなり、受光部２に入る受光量が変動する。このために、正確な脈波を測定することができなくなる。また、照射光１ａが脈波センサの外部に漏れると、照射光の効率が下がり、かつ、受光部が受け取る反射光２ａの光量が減少して正確な脈波を測定することが困難になる。そのため、本発明では、余分な反射光や漏洩光を阻止するために、発光部１の側面を遮光壁７で囲むことが好ましい。
【００７５】
この遮光壁７は、反射光や漏洩光をなくすような形状であれば、その形状は問わないが、例えば、発光部１の外周形状に沿った円筒状などの形状が好ましい。装着された指は、この遮光壁の上面に７ａの点で密着し、固定される。遮光壁７は、その受光部２側が黒色であってもよく、また、その内面が鏡面であってもよい。遮光壁の材料としては、赤外線を遮る性質を有する材料であれば特に制限はなく、例えば、赤外線を実質的に透過させないポリプロピレン系樹脂やＡＢＳ系樹脂などの熱可塑性樹脂など、又はこれらに黒色塗装などの表面処理を施したものをあげることができる。
【００７６】
発光部１の上面に赤外線透過性の絶縁体キャップ８を設け、発光部１と指３とが直接接触しないようにする。これは、発光部の通電部に影響を及ぼさないようにするためと、発光部表面の汚れを清拭しないで済むようにするためである。絶縁体キャップ８の外形は、例えば、発光体１の上方部分の形状に沿った円筒形状などの形状であればよい。絶縁体キャップ８の上面を凹レンズで構成すれば、射出光の指向性をさらに広げることができる。この絶縁体キャップの材料としては、赤外線に対して透光性が高い赤外線透過性材料であれば特に制限はなく、例えば、アクリル樹脂、ポリエチレン樹脂、ポリカーボネート樹脂、ポリスチレン樹脂などをあげることができる。また、受光部２と指３とが直接接触して指に圧力がかからないように、受光部２と指３との間に隙間を設けるような構造とすることが好ましい。
【００７７】
発光部１の発光素子と受光部２の受光素子との指向性の影響について図１６に示す。図１６（ａ）に示すように、発光部１の発光素子と受光部２の受光素子との配置が従来の指向性が強い配置の場合、発光部１の発光ダイオードの光軸がずれると、有効な検出領域もずれる。しかし、本発明におけるように、指向性の弱い発光部１の発光素子と受光部２の受光素子とを接近して配置すれば、図１６（ｂ）に示すように、発光ダイオードの光軸のずれに対する有効検出領域のずれは相対的に小さい。そのため、得られる脈波情報は正確となる。
発光部１からの照射光の射出角（半値角）αを通常５０度以上、好ましくは５０〜８５度、より好ましくは５０〜８０度とすることにより、有効検出領域のずれは相対的に少なくなる。５０度未満であると、有効検出領域のずれが大きくなり、正確な脈波データを得ることが困難になる。
【００７８】
上記脈波センサにおいて、発光部１と受光部２との距離が長くなる程、加速度脈波の波形であるａ波の振幅は小さくなり、ノイズ成分が発生し易くなって、測定波形の変形が大きくなる傾向がある。また、この距離が長い程、圧力の影響がある指部位の脈波を測定することになり、測定波形が変形し易い。そのため、発光部と受光部との距離を、所定の距離、例えば、通常８ｍｍ以内、好ましくは６ｍｍ以内に設定すれば、加速度脈波のａ波の振幅及びｂ波とａ波との比（ｂ／ａ）は適切な範囲内に納まる。この場合、光軸のずれも少なく、有効検出領域のずれも少なく、また、波形は変形し難い。なお、この距離が上記範囲を外れた動脈上流側の指部位では、動脈血管が膨らんで、ｂ／ａが小さい（絶対値が大きい）状態になり、また、下流側の指部位では、鬱血状態となって、ｂ／ａが大きい（絶対値が小さい）状態になる。また、発光部と受光部との距離の下限は、特に制限はなく、発光部と受光部との物理的な大きさや脈波センサの大きさなどに依って所望により設定できる最低距離であればよい。例えば、２〜３ｍｍ程度に設定してもよい。
【００７９】
また、絶縁体キャップが脱落しないようにして、脈波センサ本体の取り扱い性を向上させるため、図１７に示すように、絶縁体キャップ１４の下方部分につば部１４ａを設けた構造としてもよい。図１７において、１１は発光部、１１ａは発光部からの照射光、１２は受光部、１３は遮光壁、１５は指載置床面を示す。発光部１１、受光部１２、遮光壁１３、床面１５などの配置位置関係については、図１５に示す場合と同様である。また、遮光壁１３及び絶縁体キャップ１４の材料も、図１５に示す遮光壁７及び絶縁体キャップ８の材料と同様である。さらに、絶縁体キャップ１４の上面を凹レンズで構成すれば、射出光の指向性をさらに広げることができる。
【００８０】
上記した脈波センサでは、受光部は、指に圧力がかからないように、その上面が指装着部の床面と同一高さ又はそれより下になるように配置される。これにより、受光部に入射する光の割合が最も多い受光部上面位置にあたる指部分が圧迫されないようになる。例えば、受光部を脈波センサの指載置床面より１ｍｍ程度低くなるように配置すればよい。
図１８は、さらに別の反射型脈波センサの要部である指装着部の断面構造を示す図（図１８（ａ））、並びに発光部及び受光部の近傍を拡大して示す図（図１８（ｂ））であり、指を装着した状態で示す。図１８において、図１５と同じ構成要素については同じ参照符号で示す。以下、図１５の場合と異なる構成要素について説明する。
【００８１】
図１８（ａ）及び（ｂ）に示す実施の形態によれば、発光部１及び受光部２の上方に赤外線透過性の窓部９が設けられている。この窓部９は、その上面が被験者の指の腹部を載置する床面５よりも所定の距離（０．１ｍｍ以上、例えば約０．３５ｍｍ）だけ高くなるように配置されている。窓部９を、図１８（ｂ）に示すように、発光部１及び受光部２を覆うように床面の縁端部上に載置して固定しても、又は縁端部の切欠き上に載置して固定してもよく、その配置の仕方に制限はない。窓部９を設けることにより、発光部１及び受光部２と被験者の指とが直接接触しないようになる。これにより、通電部に影響が及ぶことはなく、また、発光部及び受光部表面の汚れを清拭しないで済み、メンテナンスが楽になる。
【００８２】
窓部９の外形は特に制限されず、例えば、厚さ０．５ｍｍ程度の板状でよい。この窓部９の上面を凹レンズで構成すれば、射出光の指向性をさらに広げることができる。この窓部の材料としては、赤外線に対して透光性が高い赤外線透過性材料であれば特に制限はなく、例えば、アクリル樹脂、ポリエチレン樹脂、ポリカーボネート樹脂、ポリスチレン樹脂などをあげることができる。
【００８３】
上記脈波センサによれば、発光部を、受光部より指の動脈血流の下流側に、その上面が受光部の上面よりも所定の距離突出しかつ所定の距離だけ指の腹部を載置する床面よりも突出するように配置し、そして指の先端部が発光部よりも指動脈血流のさらに下流側に装着され得るように、床面の先端部分に空間を設けてあるので、指の発光部に対する密着性が良好になると共に、指動脈血流の下流側が圧迫されても、上流側が圧迫されるのに比べて脈波への影響は少なく、再現性よく正確な脈波情報を得ることができる。
また、発光部及び受光部の上方に赤外線透過性の窓部を配置する場合、窓部の上面が指の腹部を載置する床面よりも所定の距離だけ突出するように配置してあるので、測定波形への影響は少なく、再現性よく正確な脈波情報を得ることができる。
【００８４】
【発明の効果】
本発明によれば、所定の時間連続して計測した生体脈波の波形を２次微分して算出した加速度脈波の連続する波形からａ−ａ間隔を求めて、このａ−ａ間隔の変動が心電図のＲ−Ｒ間隔の変動に対応することから、このａ−ａ間隔の変動性の評価を心電図による心拍変動性評価と同等のものとして、自律神経機能評価などに用いることができるという効果を奏する。
本発明によれば、長時間連続測定や期外収縮波の検出も可能である。
【００８５】
また、これまでの心電計を用いる方法ではできなかった期外収縮補正後のａ−ａ間隔変動のパワースペクトルが算出できるので、臨床上応用可能の範囲が広がるという利点がある。 さらに、加速度脈波測定は着衣の着脱もなく、座位にて指先の脈波を測定するだけであるから、心電図の場合よりもより簡便である。
さらにまた、加速度脈波測定装置は心電計に比べて安価である。
従って、本発明によれば、糖尿病、神経障害、脳血管疾患、冠動脈疾患、喘息、更年期障害など多方面において、合併症リスクの予測や、治療効果判定、自己管理などにも幅広く用いることができる。
【図面の簡単な説明】
【図１】ａ、ｂ、ｃ、ｄ、ｅ波の５つの成分波をもつ加速度脈波の標準波形でり、心臓の収縮期の波形図。
【図２】脈波センサを用いて測定・記録した容積脈波、速度脈波（一次微分波）、加速度脈波（二次微分波）の波形図。
【図３】心電図（ＥＣＧ）と加速度脈波（ＡＰＧ）とを同時測定した波形を示す波形図。
【図４】同時測定した加速度脈波のａ−ａ間隔と心電図のＲ−Ｒ間隔との関係を示すグラフ。
【図５】加速度脈波解析のフロー図。
【図６】加速度脈波の波高差の積分（差積分）による類似度判定について説明するためのグラフ。
【図７】ａ−ａ間隔再判定処理の流れを説明するためのグラフ。
【図８】ａ−ａ間隔の異常値判定を説明するための、拍数に対するａ−ａ間隔（ｍｓｅｃ）の関係を示すグラフ。
【図９】Ｔａａの分布全体から正常範囲を決定すると異常ａ−ａ間隔が検出されないケースはあることを説明するための、拍数に対するａ−ａ間隔（ｍｓｅｃ）の関係を示すグラフ。
【図１０】拍数に対するａ−ａ間隔（ｍｓｅｃ）の関係を示すグラフであって、単純推移平均ではうまくいかないケースを説明するためのグラフ。
【図１１】ａ−ａ間隔の異常値判定に関し、推移平均と標準ａ−ａ間隔とを比較して、拍数に対するａ−ａ間隔（ｍｓｅｃ）の関係を示すグラフ。
【図１２】健康人についての加速度脈波のａ−ａ間隔の変動係数（ＣＶａａ％）の平均値と標準偏差との関係を示すグラフ。
【図１３】図１２に示す各年代における標準的な変動係数の範囲を、［平均値−標準偏差］〜［平均値＋標準偏差］として示すグラフ。
【図１４】糖尿病患者の加速度脈波のａ−ａ間隔の変動係数を健康人の場合と比較して示すグラフ。
【図１５】本発明で用いることができる脈波センサの構造を示す断面図であり、（ａ）はこの脈波センサの要部である指装着部の模式的構造を示す断面図、（ｂ）は図１５（ａ）の発光部及び受光部の近傍の拡大断面図。
【図１６】光の指向性を示す模式図であり、（ａ）は従来技術における発光素子と受光素子との指向性の影響について示す図、（ｂ）は本発明における発光素子と受光素子との指向性の影響について示す図。
【図１７】本発明で用いることができる別の脈波センサにおいてつば部を有する絶縁体キャップの周辺の模式的構造を示す断面図。
【図１８】本発明で用いることがができる更に別の脈波センサの構造を示す断面図であり、（ａ）はこの脈波センサの要部である指装着部の模式的構造を示す断面図、（ｂ）は（ａ）の発光部及び受光部の近傍の拡大断面図。
【符号の説明】
１ 発光部 ２ 受光部
２ａ 反射光 Ｈ１ 発光部の高さ
Ｈ２ 受光部の高さ ３ 指
４ 空間 ５ 指載置床面
６ 押圧材 ７ 遮光壁
７ａ 指と遮光壁との密着点 ９ 窓部
１５ 指載置床面 α 照射光の射出角（半値角）
【発明の属する技術分野】
本発明は、脈波解析法、この解析法を利用した自律神経機能評価法及び自律神経機能評価装置に関する。
【０００２】
【従来の技術】
従来から、心拍変動性による自律神経機能を定量的指標で評価するために、非侵襲的で簡便に実施できる方法として、もっぱら心電計による計測に基づく評価法が実施されている。この心電図記録による自律神経機能判定は、１００拍、１５０拍、２分間で判定されることが多い。この場合、心電図のＲ−Ｒ間隔の変動係数を求めて、年齢別の「正常値」を参考に臨床的に診断するのが一般的である。ただし、いわゆる「正常値」は普遍的な数値ではなく、文献上の数値を参考にした上で、各医療機関で独自に決めている。
心拍変動性は加齢の影響を受け、その変動係数は加齢によって減少する。これは、これまで報告された心拍変動係数をまとめて示す表１から明らかである。
【０００３】
（表１）
上記心拍変動係数の「正常値」に性差はないといわれている。また、昼間の心拍変動係数は低く、夜間の心拍変動係数は高いという日内変動を示すが、日中変動の幅は２％以内であって有意差はなく判定に影響はないことが報告されている（例えば、非特許文献１参照）。また、心拍変動性の低下は、自律神経機能の異常、なかでも副交感神経機能低下を反映するものと考えられている。この自律神経機能の異常とは、交感神経機能と副交感神経機能とのバランスの異常である。
【０００５】
例えば、糖尿病患者では、しびれ、アキレス腱反射の減弱又は消失、振動覚障害などの体性神経障害とともに、消化管機能異常、排尿障害、起立性低血圧などの自律神経障害が多くみられるが、このような自律神経症状のない患者でも、心電図Ｒ−Ｒ間隔の心拍変動係数は減少していることが多い（例えば、非特許文献２参照）。また、一部の脳梗塞患者や、神経疾患の患者でも心拍変動係数の低下がみられる（例えば、非特許文献３）。このような自律神経機能異常は循環動態の不安定化を招き、合併症のリスクを高める。
また、脈波による心拍変動の計測も実験的には行われており、例えば、検出した容積脈波を微分した速度脈波を解析して心電図のＲ−Ｒ間隔に相当する情報を取り出す脈波解析法が知られている（例えば、特許文献１参照）。しかし、脈波による心拍変動の計測は未だ医療の現場では実用化されていない。
【０００６】
脈圧とは収縮期血圧と拡張期血圧との差であるが、これは１回拍出量と動脈系のコンプライアンスとの関数である。脈圧の波形は大動脈から末梢の動脈に行くにつれて変容していく。これは種々の部位において投射波と反射波とが合成され、共鳴が起こるからと考えられ、変容の程度は血管の性状又は特性の影響の総和とみることができる。大動脈と末梢動脈とにおける動脈内圧を比較すると、末梢動脈における平均動脈圧及び拡張期血圧は大動脈の場合よりも低い。しかし、収縮期血圧は、投射波と反射波との合成により、末梢動脈の方が高くなるため、脈圧は大きくなるという現象が起こる。それ故、末梢脈波の波形ほど中心脈波よりも凹凸が大きく波形判別がし易いという利点がある。それでも原波形は基線が安定せず変曲点の認識が困難な場合がある。
【０００７】
そこで、指先で測定する末梢脈波である容積脈波の原波形を２回微分して得られる二次微分波である加速度脈波が、より評価に適した波形として研究や臨床で利用されている。現在、多く使われている脈波計は光電式指先容積脈波計である。この脈波計の原理は、指先にヘモグロビンに吸光特異性のある波長の光を当てて、吸収光、又は反射光から血管内の血流の容積変化を求め、波形を得る方法に基づいている。容積脈波と圧脈波とは反映しているものは異なるが、波形のもつ意味は同等である。
上記末梢脈波の加速度脈波が生体のどのような生理活動（自律神経機能など）を表現しているのかについては、平均波形の変化と自律神経機能との対応の研究については知られているが、個々の成分波の変化と自律神経機能との対応を検討した研究については知られていない。
【０００８】
加速度脈波は心臓の収縮期の波形であり、ａ〜ｅ波の５つの成分波をもつ（図１）。これらの成分波は、生体の条件により、また、加齢に応じて一定の変化をする。頂点ａを有する成分波（ａ波）は原波形の立ち上がり部分の波形に相当し、頂点ｅを有する成分波（ｅ波）は原波形の収縮終期のノッチ部分に相当する。それゆえ、期外収縮時や頻脈時にはａ−ｅ間隔は短縮し、各成分波も変形する。波形の基線より上を正、下を負の象限としたとき、ａ波は基線より常に上に位置する陽性波であり、ｂ波は基線より常に下に位置する陰性波であり、ｃ波、ｄ波、ｅ波は生体の条件により陽性又は陰性に変化する成分波である。
【０００９】
また、加速度脈波の波形は、加齢によって変化していくので、生理的変化と病的変化の判別も難しい。測定器械、測定環境、年齢などによって異なる波形を正しく評価するには、どのような条件にも適応できるような共通の波形評価法を確立しなければならない。すなわち、測定環境が同一であれば、どの器械にも、どの年齢層にも共通して使用できる共通の基準を作成し、それに基づいて評価することが必要であるが、いまだ有用な評価法は確立されていない。
【００１０】
【非特許文献１】
景山茂他：心電図Ｒ−Ｒ間隔変動、内科、５５： ２４２−２４６（１９８５）
【非特許文献２】
Ｂｅｒｎａｒｄｉ Ｌ．， Ｒｉｃｏｒｄｉ Ｌ．， Ｌａｚｚａｒｉ Ｐ． ｅｔ ａｌ：Ｉｍｐａｉｒｅｄ ｃｉｒｃａｄｉａｎ ｍｏｄｕｌａｔｉｏｎ ｏｆ ｓｉｍｐａｔｈｅｔｉｃ ａｃｔｉｖｉｔｙ ｉｎ ｄｉａｂｅｔｅｓ， Ｃｉｒｃｕｌａｔｉｏｎ， １９９２； ８６： １４４３−１４５２
【非特許文献３】
宮田正和、松浦達雄：ＣＶＲ−Ｒ、心身医療、Ｖｏｌ． ３， Ｎｏ． ８ １９９１
【特許文献１】
特開２００１−７０２６５号公報（特許請求の範囲）
【００１１】
【発明が解決しようとする課題】
自律神経活動評価としての心拍変動測定の基礎は、心臓の拍動の原因である心電パルスの発生間隔の変動である。特に、心室収縮の開始を表すＲ波の間隔（Ｒ−Ｒ間隔）を評価する。しかし、心電計による計測は簡便ではない。特に、体温計や血圧計による計測の場合と異なり、患者が家庭で測定して経過を見るということはできない。従って、上記のような心電図のＲ−Ｒ間隔の変動係数を求めて、自律神経機能を評価する簡便な方法の開発が求められている。
また、心電図のＲ−Ｒ間隔の変動係数の代わりに、容積脈波や速度脈波を利用して自律神経機能を評価しようとしても、基線変動や体動ノイズなどによる脈の誤診や検出ミスがあり、特に脈が小さい場合、判定が困難であるという問題もある。
【００１２】
そこで、本発明の課題は、心電図に基づくのではなく、加速度脈波を用いて自律神経機能などを評価するための加速度脈波の解析法、この解析法を利用する自律神経機能評価方法及び自律神経機能評価装置を提供することにある。
【００１３】
【課題を解決するための手段】
本発明の脈波解析法は、生体の脈波から心電図のＲ−Ｒ間隔の変動に対応する情報を求める脈波解析法であって、所定の時間連続して計測した脈波の波形を２次微分して加速度脈波を算出し、この加速度脈波の連続する波形からａ−ａ間隔を求め、このａ−ａ間隔の変動を心電図のＲ−Ｒ間隔の変動に対応する間隔とすることを特徴とする。
【００１４】
本発明の脈波解析法はまた、生体の脈波から心電図のＲ−Ｒ間隔の変動に対応する情報を求める脈波解析法であって、過去の加速度脈波、好ましくは過去数拍から数十拍の加速度脈波から求めた基準（標準）となる波形を標準加速度脈波とし、この標準加速度脈波と測定加速度脈波との波形の類似度を評価して個別加速度脈波を決定し、この個別加速度脈波の連続する波形からａ−ａ間隔を求め、このａ−ａ間隔の変動を心電図のＲ−Ｒ間隔の変動に対応する間隔とすることを特徴とする。
上記波形の類似度の評価を、測定加速度脈波のａ波候補の波高値が標準加速度脈波のａ波の波高値と同一若しくは両者の中間値になるように規格化した波形を用いて行い、また、標準加速度脈波のａ波と測定加速度脈波のａ波候補との時間位置を合わせ、規格化した測定加速度脈波と標準加速度脈波との波高差を求めて行うことを特徴とする。
【００１５】
上記波形の類似度の評価を、個別加速度脈波のａ波判定条件として、測定加速度脈波をそのａ波候補の波高値と標準加速度脈波のａ波波高値との中間値若しくはａ波波高値で規格化した個別加速度脈波と標準加速度脈波との波高差の積分値又は多重積分値を用いて算出した類似度に基づいて行うことを特徴とする。
上記脈波解析法において、隣り合うａ−ａ間隔の変化が一定範囲を外れるものを除外したａ−ａ間隔について、これらの値により標準ａ−ａ間隔を構成して、ａ−ａ間隔の異常を判定することを特徴とする。
【００１６】
上記標準ａ−ａ間隔として、過去のデータの平均若しくは重み付き平均を用い、異常判定がでたデータに対して、未来データを含む平均若しくは重み付き平均で再判定することを特徴とする。以下、「平均若しくは重み付き平均」を簡略的に「（重み付き）平均」と言うこともある。
上記脈波解析法において、評価している時点のａ−ａ間隔が、過去の平均的なａ−ａ間隔、又は未来よりも過去のａ−ａ間隔データの重みを大きくした重み付き推移平均、又は直前のａ−ａ間隔を基準として、一定の基準値を超えるとき、その時点の直後から順次、隣り合うａ−ａ間隔の変化が一定範囲内に入るものを数点選び出し、これらの値により標準ａ−ａ間隔を更新して、ａ−ａ間隔の異常を判定することを特徴とする。
【００１７】
上記脈波解析法において、生体の脈波を、発光部と受光部とを有する反射型脈波センサであって、発光部が、受光部より指の動脈血流の下流側に、その上面が受光部の上面よりも突出しかつ所定の距離だけ指の腹部を載置する床面よりも突出するように配置され、そして指の先端部を発光部よりも指動脈血流のさらに下流側に装着するための空間が床面の先端部分に設けられている脈波センサに指先を装着して測定することを特徴とする。
脈波の測定は、上記脈波センサの発光部がその上面が指の腹部を載置する床面より０．１〜１．５ｍｍ突出するように、発光部が配置されている脈波センサを用いて行われることを特徴とする。
【００１８】
上記脈波の測定は、脈波センサの受光部の上面が発光部の上面と同一レベルか若しくはそれより下になるように、また、所定の距離だけ床面より低くなるように、受光部が配置されている脈波センサを用いて行われることを特徴とする。
上記脈波の測定は、脈波センサの指の先端部を装着する空間の床面と対向する面に押圧材が取付けられ、この押圧材で発光部よりも指動脈血流のさらに下流側の指の先端部の上面を押さえるように構成されている脈波センサを用いて行われることを特徴とする。
上記押圧材は、この押圧材で発光部よりも指動脈血流のさらに下流側の指の先端部の上面を５０〜２００ｇＷの圧力で押さえるように構成されていることを特徴とする。
【００１９】
上記脈波の測定は、脈波センサの発光部の側面が、内面が照射光に対する反射特性を有する遮光壁で囲まれ、そして遮光壁上端が指の腹部を載置する床面から０．１〜１．５ｍｍ突出するように構成されている脈波センサを用いて行われることを特徴とする。
上記脈波の測定は、脈波センサの発光部及び受光部の上方に、赤外線透過性の窓部を、窓部上面が指の腹部を載置する床面よりも所定の距離だけ突出するように配置してある上記脈波センサを用いて行われることを特徴とする。
上記窓部の上面が指の腹部を載置する床面より０．１〜０．５ｍｍ突出するように配置されていることを特徴とする。
【００２０】
また、本発明の自律神経機能評価法は、上記脈波解析法に基づき、連続するａ−ａ間隔について、隣り合うａ−ａ間隔の最初のａ−ａ間隔が標準ａ−ａ間隔より正常範囲を外れて短い場合で、かつ、前回のａ−ａ間隔と今回のａ−ａ間隔との和が標準ａ−ａ間隔の１倍若しくは２倍を基準にしてａ−ａ間隔の変動の分布から定まる所定の範囲に入るとき、期外収縮と判定することを特徴とする。
この自律神経機能評価法はまた、ａ−ａ間隔の統計量に基づいて自律神経機能を評価する方法であって、上記脈波解析法に基づき、ａ−ａ間隔が正常範囲を外れる場合はこれを除外してａ−ａ間隔の統計量を算出し、これに基づいて自律神経機能を評価することを特徴とする。
【００２１】
この自律神経機能評価法はまた、ａ−ａ間隔の時間変動を周波数解析して自律神経機能を評価する方法であって、隣り合うａ−ａ間隔について、最初のａ−ａ間隔が、上記脈波解析法に基づき、標準ａ−ａ間隔より正常範囲を外れて短い場合で、かつ、最初のａ−ａ間隔と次のａ−ａ間隔との和が、標準ａ−ａ間隔の１倍若しくは２倍を基準にして所定の範囲（例えば、０．７５〜１．２５若しくは１．５〜２．５倍）に入る時、１倍を基準とする場合には、最初のａ−ａ間隔を最初のａ−ａ間隔と次のａ−ａ間隔との両者の和に補正して次のａ−ａ間隔を削除し、また、２倍を基準とする場合には、それぞれ二つのａ−ａ間隔を両者の和の二分の一に補正して、周波数解析することにより自律神経機能を評価することを特徴とする。
【００２２】
上記自律神経機能評価法において、ａ−ａ間隔の変動性を各年代別標準参考値を参照しながら評価することにより、自律神経機能を評価することを特徴とする。
上記自律神経機能評価法において、上記脈波解析法に基づき、長時間連続測定や期外収縮波の検出を行うことにより、また、期外収縮除去後の推定変動係数を算出することにより、糖尿病、神経障害、脳血管疾患、冠動脈疾患、喘息、更年期障害などを含めて、合併症リスクの予測や、治療効果判定、自己健康管理を行うことを特徴とする。
【００２３】
さらに、本発明の加速度脈波の変動解析による自律神経機能評価装置は、生体の脈波を検出して脈波の大きさに応じた信号を出力する脈波計測手段と、脈波計測手段により得られた脈波の波形を２次微分して算出された加速度脈波の波形からその波形パラメータを解析する波形パラメータ解析手段とを備え、この波形パラメータ解析手段が、所定の時間連続して計測した脈波の波形を２次微分して加速度脈波を算出する手段と、この加速度脈波の連続する波形から心電図のＲ−Ｒ間隔の変動に対応するａ−ａ間隔の変動を求める手段とを有することを特徴とする。
上記自律神経機能評価装置は、さらに、隣り合うａ−ａ間隔の変化が一定範囲を外れるものを除外したａ−ａ間隔について、これらの値により標準ａ−ａ間隔を構成して、ａ−ａ間隔の異常を判定する手段を有することを特徴とする。
【００２４】
【発明の実施の形態】
本発明によれば、加速度脈波を変動解析し、この変動解析に基づいて自律神経機能などを評価し、各種医療診断や健康診断を行うことができる。この変動解析は、生体の指先に取り付けて使用される脈波センサの計測・検出結果に基づき行われる。この脈波センサからなる脈波測定装置は、所定の加速度脈波情報からａ−ａ間隔を抽出するための脈波解析プログラムが組み込まれたマイクロコンピュータを有している。
【００２５】
図２に脈波センサ（（株）ユメディカ社製のＡＲＴＥＴＴ（登録商標））を用いて測定・記録した容積脈波、速度脈波（一次微分波）、加速度脈波（二次微分波）を示す。この加速度脈波は、図１に示すようにａ波、ｂ波、ｃ波、ｄ波、ｅ波の５つの成分波をもつ。これら５つの成分波で構成される加速度脈波は心臓収縮期の波であり、これらの成分波は、生体の状態により、また、加齢に応じて一定の変化をしていく。ただし、ｅ波は収縮と拡張の境界波である。なお、この脈波センサを用いて測定することにより、体動などによるノイズの発生を極力押さえた連続波形データを得ることができる。
【００２６】
図３に、心電図（ＥＣＧ）と加速度脈波（ＡＰＧ）とを同時測定した波形を示す。図３から明らかなように、心電図Ｒ波の出現に一定の時間的な遅れ（α）を伴って、加速度脈波のａ波が出現している。このとき、加速度脈波のａ−ａ間隔（ＡＰＧ／ａａ間隔（Ｔａａ））は、心電図のＲ−Ｒ間隔（ＥＣＧ／ＲＲ間隔）にほぼ一致する（図４）。
【００２７】
［加速度脈波ａ−ａ間隔の計測と変動係数算出］
本発明者らの知見によれば、血管の老化の程度を評価するための加速度脈波波形評価には、１０秒から２０秒程度までの測定期間において容積脈波が安定している期間の１０拍程度が検出できれば充分であった。しかしながら、自律神経機能評価のためには、１００拍以上、場合によっては１５分以上の連続測定が必要であって、脈拍間隔の１ｍｓｅｃ以内の精度での連続測定が要求される。従って、一拍毎の個別加速度脈波の精度の高い検出と、ａ波の時間位置検出精度が問題となる。
【００２８】
上記したように、従来の自律神経活動評価における心拍変動測定の基礎は、心臓の拍動の原因である心電パルスの発生間隔の変動、特に、心室収縮の開始を表すＲ波の間隔（Ｒ−Ｒ間隔）の変動を測定し、評価することにある。このＲ−Ｒ間隔を、他の心拍検出手段で代用することも様々に行われている。例えば、心音や、頸動脈圧脈波、指尖容積脈波などで心拍数若しくは脈拍数を測定することが行われている。従来の方法では、一般的には、これらの拍動を電気信号に変換したもの（一般的には電圧）が一定レベルを超えるときにパルスを発生させて、１分間の拍数をカウントしているに過ぎない。
しかるに、数％程度の変動を計測しなければならない自律神経機能の評価のためには、ミリ秒程度の精度で個々の拍動の間隔を計測する必要があるので、パルスの発生位置が問題となる。
【００２９】
容積脈波は、基線が変動し、安定した検出が困難である。また、時間位置精度もでない。気温が低い場合や体調の影響により、脈波の振幅が非常に小さい場合があり、この場合、Ｓ／Ｎ（信号／雑音）比が小さくなり、脈の検出そのものができない場合がある。速度脈波や加速度脈波の検出も試みられているが、Ｓ／Ｎ比が悪い場合には、ノイズのピークが多く発生し、検出すべきピークの同定が困難である。
【００３０】
そして、この計測された心拍間隔が心電図のＲ−Ｒ間隔と一致していれば問題はないが、このことは必ずしも自明ではない。例えば、指尖容積脈波の極小点や極大点を拍動の時間位置として検出することが考えられるが、これらは、基線変動の影響を受けるし、個々の心拍の駆出力などの生理的な変動の影響も受けやすい。そして、検出する心拍変動係数が小さく、心電図のＲ−Ｒ間隔が１％以内となる自律神経障害を有効に検出するためには、少なくとも心電図のＲ−Ｒ間隔と数ミリ秒以内の精度での対応付けができている必要がある。しかし、従来の方法では、このような精度を達成することは困難である。
【００３１】
また、実用的な応用において、パーソナルコンピュータによる自動解析の出力結果が信頼できるものであるためには、異常なａ−ａ間隔を検出したときに、誤検出（ノイズ若しくはａ波以外のピークをａ波として検出した場合）、又は検出ミス（実際に存在するａ波を検出し損なった場合）によってａ−ａ間隔を間違って計測したのか、あるいは実際のａ−ａ間隔を検出しているのかが問題となる。
平均波形を算出して波形の特徴量だけを求める場合には、確実に個々の脈波を検出して検出できたものだけで平均を取ればよいが、連続するａ−ａ間隔の変動を検出し、ａ−ａ間隔の異常、正常の判定をするためには、個々の脈波の検出精度を飛躍的に高める必要がある。実際上は、測定波形とａ波検出位置とを同時に表示して確認するのが最善であり、如何に解析ソフトの解析精度が向上しようとも、視覚的な確認の必要性は残ると考えられるが、誤検出や検出ミスはできうる限り排除できるように構成すべきである。
【００３２】
以下、個別加速度脈波の検出・解析について、図５のフロー図を参照して説明する。
１．標準加速度脈波パラメータの初期値算出
測定開始時に、類似度判定のための標準加速度脈波とそのパラメータ及びａ−ａ間隔異常判定のためのａ−ａ間隔変動幅の初期値を求める必要がある。このため、波形をモニターしながら容積脈波の基線の変動及びノイズが少ないことを確認して測定を始め、容積脈波のピークやボトム若しくは立ち上がりを検出して脈波を切り出して、その時間間隔からａ−ａ間隔とその変動幅を求めると共に、平均波形を求め、これから標準加速度脈波とａ波波高値の初期値を求める。または、速度脈波や加速度脈波が安定していることを確認して、同様の処理をする。
【００３３】
以後、上で算出した標準値を用いて、測定開始時点から、改めて次の解析を行う。
２．ａ波候補の抽出
（１）測定加速度脈波の極大点を求め、求めた極大値とａ波の標準波高値とを比較して、一定範囲（例えば、６０〜１８０％）に入るものをａ波候補とする。６０％未満のものはノイズ若しくはｅ波を検出しているかもしれず、１８０％を超えるものはノイズの可能性が高いという問題がある。
【００３４】
３．個別加速度脈波のａ波判定
上記したように、個別加速度脈波のａ波判定条件として、ａ波の波高値で規格化した個別加速度脈波と標準加速度脈波との波高差の積分値若しくは多重（二重）積分値を用いて類似度を算出し、判定している。
（１）測定加速度脈波の波形を、標準加速度脈波のａ波の波高値若しくはこのａ波の波高値とａ波候補の波高値との中間値で規格化する。この規格化は、測定加速度脈波の振幅に［（標準加速度脈波のａ波の波高値）＋ｑ＊（ａ波候補の波高値）］／［（１＋ｑ）＊（ａ波候補の波高値）］（ｑは０以上の実数値であり、実用的には１程度である）を乗じて行う。
【００３５】
（２）測定加速度脈波のａ波候補点と標準加速度脈波のａ波との時間（ｔ）位置を合わせ（Ｔａとする）、ａ波の波高値で規格化した測定加速度脈波と標準加速度脈波との波高差を求める。
（３）Ｔａ−ｔ１〜Ｔａ＋ｔ２の範囲で、求めた波高差の絶対値を積分して類似度とする。ここで、ｔ１は、例えば０．１秒とし、ｔ２は、例えばａ−ｅ間隔程度の時間とすればよい。この範囲を広く取りすぎると、期外収縮などの不整脈で、ｅ波直後に次ぎのａ波が出現した場合など、拡張期の標準波形との差が大きく出て、類似度が小さくなり、脈が検出できなくなる。
この場合、Ｔａ−ｔ１〜Ｔａ＋ｔ２の範囲で、波高差を積分したものの絶対値を更に積分（二重積分）して類似度としてもよい。
【００３６】
（４）類似度が高いとき、ａ波候補をａ波として確定する。類似度が低いときは、ａ波候補をａ波と判定せず、ノイズと判定して、再度ａ波候補の抽出を行う。
（５）今回確定したａ波と前回のａ波との時間差を求めて、これをａ−ａ間隔Ｔａａとして算出する。
【００３７】
上記した加速度脈波の波高差の積分値若しくは多重積分値を用いて類似度を評価する根拠は以下の通りである。
（１）一般的に、加速度脈波は基線変動に強いので、波高差で評価しても、容積脈波、速度脈波の基線変動の影響を受けないというメリットがある。しかし、微分波形は一般に高周波ノイズの影響を受けやすく、波高差での評価ではスパイク状のノイズには対応できない。
【００３８】
（２）加速度脈波の積分は基本的には速度脈波なので、速度脈波の波高差を類似度とすることが考えられる。しかしながら、この場合、速度脈波の基線変動が大きいと、判定ミスが発生する。これは、速度脈波の波形がまったく同形でも、基線変化分だけの波高差が発生するからである。一方、速度脈波の微分波形である加速度脈波の値は、速度脈波の波形の勾配であり、積分値は速度脈波の波高値の変化を表し、基線変動の影響を受けない。そこで、速度脈波に基線変動があっても、変動の影響が除去される微分波形（加速度脈波）から積分値を求めて、これから類似度を算出して評価することにより、安定した個別脈波の判定が可能となる。
【００３９】
（３）速度脈波にスパイク状のノイズが無視できない場合には、上記第３（３）項に示した二重積分による類似度評価が有効である。類似度を求めるための積分の変形として、特に限られるわけではなく、上記した絶対値の積分以外に、波形データを二乗したものを積分し、これの平方根を算出してもよい。
（４）なお、速度脈波の波高差を用いて類似度を評価する場合の上記問題点（基線変動が大きい）を解決するため、速度脈波を脈拍周期で積分した値を速度脈波の基線として補正することが考えられるが、不整脈がある場合にはうまくいかない。
【００４０】
次に、上記した加速度脈波の波高差の積分（差積分）による類似度判定について、図６を参照して詳細に説明する。図６の横軸は時間（ｍｓｅｃ）であり、縦軸は波高値である。
測定加速度脈波Ａ１と標準加速度脈波Ａ０との差Ａ１−Ａ０のゼロ点Ｔ（ｎ）では、標準加速度に対応する速度脈波（標準速度脈波）Ｖ０の勾配と、測定加速度脈波に対応する速度脈波（測定速度脈波）Ｖ１の勾配とは一致しており、このときのＶ０に対するＶ１の波高差（Ｖ１−Ｖ０の値）をΔＶ（ｎ）とする（図の場合、この値は負になっている）。差Ａ１−Ａ０が正の期間では、測定加速度脈波は標準加速度脈波に対して相対的に上昇し、負の期間では下降する。例えば、正の期間（Ｔ（０）からＴ（１）までの期間）における差分Ａ１−Ａ０の積分値は、速度脈波の波高値の上昇ΔＶ（１）−ΔＶ（０）を表しており、負の期間（Ｔ（１）からＴ（２）までの期間）におけるＡ１−Ａ０の積分値は、速度脈波の波高値の下降ΔＶ（２）−ΔＶ（１）を表している。
【００４１】
一般的に、速度脈波は基線変動により全体のレベルが変化していることがある。図６に示した場合では、測定速度脈波Ｖ１の基線は低くなっており、このような場合、速度脈波の波高値で脈拍を検出すると検出ミスが発生することがある。標準速度脈波Ｖ０に対する測定速度脈波Ｖ１の波高差の変化は、このような基線変動の影響を受けない。従って、差分Ａ１−Ａ０の絶対値の積分は、速度脈波の基線変動の影響を受けずに、標準速度脈波Ｖ０に対する測定速度脈波Ｖ１の形状差の目安とすることができる。
【００４２】
しかしながら、波高差で評価する場合には、高周波ノイズやパルス状のノイズがあると、類似度の評価値が悪くなる。この問題は、差分Ａ１−Ａ０（絶対値ではない）の積分を行い、これの絶対値の積分から類似度を評価することにより解決できる。この場合は、基線変動の影響を除いて容積脈波の形状差を評価していることになる。
また、リアルタイム性を重視する場合、計算量削減のために上記第３（２）〜（４）項と第３（５）項とは順番を逆にしてもよい。つまり、ａ波候補と前回ａ波との時間間隔が正常範囲の時、類似度を評価せずにａ波と判定し、ａ−ａ間隔が異常なときのみ、類似度を評価することもできる。
【００４３】
４．ａ−ａ間隔の異常値判定法
測定中にａ−ａ間隔の異常が生じた場合、以下のようにしてその異常が判定される。
隣り合うａ−ａ間隔Ｔａａの変化ＤＴａａ、すなわち（ａ−ａ間隔変動）ＤＴａａ［ｎ］＝Ｔａａ［ｎ］−Ｔａａ［ｎ−１］が、一定範囲（例えば、ＤＴａａの標準偏差Ｓｄｖ＿ＤＴａａの２倍以内）を外れるＴａａを除外したＴａａについて、これらの値により標準ａ−ａ間隔を構成して、ａ−ａ間隔の異常を判定する。この標準ａ−ａ間隔は判定対象データ前後のデータの推移平均を用いることができるが、リアルタイム性考慮する場合、過去データの（重み付き）平均を用い、異常判定がでたデータに対して、未来データを含む（重み付き）平均で再判定すればよい。
【００４４】
リアルタイム性を考慮して、アルゴリズムを構成する場合は、評価している時点でのａ−ａ間隔が、過去の平均的なａ−ａ間隔、若しくは直前のａ−ａ間隔（標準ａ−ａ間隔）を基準として、一定の基準値（例えば、４＊Ｓｄｖ＿ＤＴａａ）を超えるとき、その時点の直後から順次、隣り合うａ−ａ間隔の変化ＤＴａａが一定範囲（例えば、２＊Ｓｄｖ＿ＤＴａａ）内に入るＴａａを数点選び出し、これらの値により標準ａ−ａ間隔を更新して、ａ−ａ間隔の異常を判定してもよい。
【００４５】
図７に基づいて、ａ−ａ間隔再判定処理の流れを説明する。
▲１▼ｔ＝０までのデータを用いて求めた標準ａ−ａ間隔をＳｔ＿Ｔａａ（補正前）とする。この標準ａ−ａ間隔の算出式は、例えば、後出（段落番号［００５４］）の重み付き平均による更新の際の算式による。
▲２▼ｔ＝１の時のＴａａ（Ｔａａ（１））をＳｔ＿Ｔａａ（補正前）と比較して異常と判定する。
▲３▼ＤＴａａ（２）＝Ｔａａ（２）−Ｔａａ（１）を評価して、Ｔａａ（２）を平均値算出データとする。この評価は、前述のように、ＤＴａａの標準偏差Ｓｄｖ＿ＤＴａａから定めた閾値よりも、評価しているＤＴａａが小さければ平均値算出データとし、大きければ平均値算出データに含めないものとする。
【００４６】
▲４▼同様にして、ＤＴａａ（３）、ＤＴａａ（４）、ＤＴａａ（５）を評価し、Ｔａａ（３）、Ｔａａ（４）は平均値算出データから除外し、Ｔａａ（５）を平均値算出データとする。
▲５▼Ｓｔ＿Ｔａａ（補正前）とＴａａ（２）、Ｔａａ（５）の平均値を求めてＳｔ＿Ｔａａ（補正後）とする。
▲６▼Ｔａａ（１）をＳｔ＿Ｔａａ（補正後）と比較し、▲２▼で一旦異常と判定した結果を、正常と再判定する。
【００４７】
図８は、実測データに本アルゴリズムを適用した場合の標準ａ−ａ間隔の推移を示す。図には、Ｔａａの変化（前回Ｔａａとの差）も示してあるが、これから分かるように、前回Ｔａａの差のみでは異常Ｔａａを判定することができない。
心拍間隔の変化は、カオス的であるともフラクタル的な変化パターンを示すとも言われており、非常に複雑であって、単一の判定基準で判定することが困難である。図８に示すようなステップ状の変化も正常な変化として解析データに含める必要がある。
一般によく用いられる単一の基準を適用する場合の問題点を次に説明する。
正常範囲を標準ａ−ａ間隔の±１５％（８５％〜１１５％）のような一定範囲とすることが考えられるが、ａ−ａ間隔の変動は個人差が大きいので、変動が大きい被験者の場合には、正常な変動を異常とする検出ミスが発生することが考えられ、逆に変動が小さい被験者の場合には、異常な変動を正常とする誤検出が発生することが考えられるので、好ましくない。
【００４８】
ａ−ａ間隔の平均値を基準にして、ａ−ａ間隔の分布全体から正常ａ−ａ間隔の範囲を定めた場合、上記と異なり、次のような問題点がある。ゆっくりした大きな変動において、例えばＴａａが長いときに変動範囲内の短いＴａａが発生すると、異常ａ−ａ間隔と判定すべきところが正常ａ−ａ間隔と判定されることになる（図９）。
また、判定対象になっているａ−ａ間隔測定値Ｔａａ［ｎ］の前後数拍の平均値を標準ａ−ａ間隔とすることの問題点、すなわち、推移平均から求めた標準ａ−ａ間隔からの差のみでは異常Ｔａａが判定できない理由は次の通りである（図１０）。
【００４９】
呼吸周期の変動幅が小さい場合に、大きなゆっくりした変化があるときで、呼吸周期の変動幅を超える階段状の変化があるとき、変化点前後の中間値が標準ａ−ａ間隔となる。このとき、この変化を正常と見なすのであれば、正常範囲の限界ＴｈｒＤＤＴａａはこの階段状の変化の半分以上にする必要がある。ところがこの場合、単発のａ−ａ間隔異常でこのＴｈｒＤＤＴａａ以内のものは、検出できないことになる。
また、図１０に示すように、必ずしも異常とはいえないのに異常と判定される場合がある。
【００５０】
図１１に、ａ−ａ間隔異常値判定に関し、上記で求めた標準ａ−ａ間隔と推移平均法で求めた標準ａ−ａ間隔とを、実測データについて比較したグラフを示す。この図から明らかなように、推移平均では、データ列に階段状のギャップがあるような場合に、正常データを異常データと判定するケースが発生するが、本法により判定ミス（誤判定）を回避することができる。
【００５１】
５．異常データの処理
（１）ａ−ａ間隔の統計解析（変動解析）においては、異常と判定されたａ−ａ間隔は全て除外して解析する。
（２）スペクトル解析においては、以下の判定条件により期外収縮と判定された場合は、補正して解析を行う。その他の異常データがある場合は、異常データが含まれないよう、解析区間を分割して解析すべきである。
【００５２】
隣り合うａ−ａ間隔の最初のａ−ａ間隔が、標準ａ−ａ間隔より正常範囲を外れて短い場合で、前回ａ−ａ間隔と今回ａ−ａ間隔との和が、標準ａ−ａ間隔の１倍若しくは２倍を基準にして所定の範囲（例えば、０．７５〜１．２５若しくは１．５〜２．５倍）に入るとき、期外収縮と定める。
上記所定の範囲は、望ましくは隣り合うａ−ａ間隔の差の変動の分布から定める。例えば、隣り合うａ−ａ間隔の差の標準偏差をＳｄｖ＿ＤＴａａとして、±ｎ＊Ｓｄｖ＿ＤＴａａ（ｎの値は前者の場合、４〜５、後者の場合８〜１０程度が望ましい）を所定の範囲とする。
【００５３】
６．標準波形及び標準パラメータの更新
前回測定のａ波Ｔａ［ｎ−１］と今回測定のａ波Ｔａ［ｎ］との差Ｔａａ［ｎ］＝Ｔａ［ｎ］−Ｔａ［ｎ−１］が、上記したように定めた正常範囲に入るとき、標準ａ−ａ間隔（Ｓｔ＿Ｔａａ）、ａ−ａ間隔の変動幅（Ｓｄｖ＿ＤＴａａ）、標準ａ波波高値（Ｓｔ＿ＡＣＣａ）及び標準加速度脈波波形を更新する。
【００５４】
この更新は、例えば、平均値若しくは重み付き平均値による。重み付き平均による更新は次のような算式よってもよい。
Ｓｔ＿Ｔａａ（今回）＝（（ｋ−１）＊Ｓｔ＿Ｔａａ（前回）＋Ｔａａ（正常測定値））／ｋこの場合、ｋ＝１のときは、今回測定のａ−ａ間隔を次回の標準ａ−ａ間隔とすることに相当する。
【００５５】
標準加速度脈波については、前回ａ波Ｔａ［ｎ−１］を基準とする前回加速度脈波で、標準加速度脈波の波形を更新する。正常範囲に入らないときは更新しない。
上記のようにして個別加速度脈波を検出・解析すれば、期外収縮波やその他の異常値を除去して推定変動係数を算出できるので、期外収縮が存在する事実と変動係数評価とを別にでき、臨床上有用である。
【００５６】
［加速度脈波のａ−ａ間隔の変動係数の標準参考値］
次に、加速度脈波のａ−ａ間隔の変動係数の標準参考値、糖尿病患者のａ−ａ間隔の変動係数、加速度脈波による自律神経機能評価の価値について説明する。
図４に示すように、心電図のＲ−Ｒ間隔と加速度脈波のａ−ａ間隔とは一致している。通常、心電図は仰臥位で測定するのに対して、加速度脈波は座位で測定する。仰臥位と座位のふたつの姿勢においては、自律神経バランスの状態が異なるため、同一人であっても、それぞれの変動係数が異なることが考えられる。そこで、健康人１２１人について、加速度脈波を２分間記録して、ａ−ａ間隔の変動係数（ＣＶａａ％）の平均値と標準偏差とを年令別に算定した。年代別の変動係数を図１２に示す。また、各年代における標準的な変動係数の範囲を、［平均値−標準偏差］〜［平均値＋標準偏差］として図１３に示す。
これらの図１２及び１３から明らかなように、加速度脈波のａ−ａ間隔の変動係数は、加齢と共に減少するが、その標準偏差はほぼ同じ範囲内に入ることが分かる。
【００５７】
［糖尿病患者のａ−ａ間隔の変動係数］
上記加速度脈波のａ−ａ間隔の変動係数について、糖尿病患者と健康人とを比較した。糖尿病患者２６人の加速度脈波を２分間記録して、上記したようにしてａ−ａ間隔の変動係数を算出した。図１４に健康人（Ｎｏｒｍａｌ）と糖尿病患者（ＤＭ）との変動係数を比較して示す。図１４から明らかなように、４０代及び５０代において糖尿病患者の変動係数の低下が見られた。
【００５８】
［加速度脈波による自律神経機能評価法の価値］
上記したように、加速度脈波ａ−ａ間隔の変動と心電図Ｒ−Ｒ間隔の変動とは対応する。従って、加速度脈波のａ−ａ間隔の変動性の評価は心電図による心拍変動性評価と同等のものであるといえる。加速度脈波のａ−ａ間隔の変動係数は、心電図のＲ−Ｒ間隔の変動係数よりも、年齢により０．２〜１．５程度高くでる可能性があるが、これは、座位測定のために仰臥位よりは呼吸の影響が大きいためであると考えられる。そこで、ａ−ａ間隔の変動係数は、その年代別標準参考値を参照しながら心拍変動性を判定することにより、自律神経機能評価に用いることができる。
【００５９】
上記システムでは、長時間連続測定や期外収縮波の検出も可能である。また、このシステムは、期外収縮除去後の推定変動係数を算出できるので、臨床上応用可能の範囲が広がるという利点がある。さらに、加速度脈波測定は着衣の着脱もなく、座位にて指先の脈波を測定するだけであるから、心電図の場合よりもより簡便である。さらにまた、加速度脈波測定装置は心電計に比べて安価である。従って、糖尿病、神経障害、脳血管疾患、冠動脈疾患、喘息、更年期障害など多方面において、合併症リスクの予測や、治療効果判定、自己管理などに幅広く用いることができる。
【００６０】
上記加速度脈波を測定するための脈波センサについて以下説明する。
本発明で用いることができる脈波センサは、発光部と受光部とを有し、発光部と受光部とにより被験者の指の脈波を測定する反射型の脈波センサであって、波形がかわってもａ波の検出ができるものであれば特に制限はなく、従来の脈波計測装置を用いることができる。例えば、本発明者らが開発した脈波センサであって、発光部を、受光部より指の動脈血流の下流側に、その上面が受光部の上面よりも突出しかつ所定の距離だけ指の腹部を載置する床面よりも突出するように配置し、そして指の先端部を発光部よりも指動脈血流のさらに下流側に装着するための空間が床面の先端部分に設けられているようなものであってもよい。この場合、このように構成することにより、指の発光部に対する密着性が良好になると共に、指動脈血流の下流側が圧迫されても、上流側が圧迫されるのに比べて脈波への影響は少なく、再現性よく脈波情報を得ることができる。また、指と発光部との接触面積は発光部の上面の面積と等しくなり、指にかかる接触面積が小さいので、波形の変化はほとんどない。
【００６１】
上記発光部は、その上面が指の腹部を載置する床面より通常０．１〜１．５ｍｍ程度、好ましくは０．２〜１．０ｍｍ程度、さらに好ましくは０．３〜０．５ｍｍ程度突出するように配置されてもよい。発光部をこのように配置すると、指腹部の皮膚面が発光部上面から下に被さるので、測定データに対する外乱光や漏洩光や反射光の影響を小さくすることができると共に、被験者が指を装着する際に、その突出部に触れることによりセンサ位置を認識して、指を所定の位置に載置し易いという利点がある。しかし、０．１ｍｍ未満であると、センサ位置を確認し難いので、指先を所定位置に載置し難いと共に、測定データに対する反射光の影響が大きくなる。また、１．５ｍｍを超えると、指皮膚面が床面から浮くため、不安定な装着状態になると共に、指を載置した時の指に対する圧力により波形の変形が生じて再現性が悪くなり、測定される脈波データにバラツキが生じて、正確な脈波情報を得ることが難しくなる。
【００６２】
上記受光部は、その上面が指の腹部を載置する床面と同じレベルになるように配置されているか、又は所定の距離だけ床面より低くなるように配置されていることが好ましく、これにより発光部に対する指の密着性がより良好になる。
上記指の先端部を装着する空間の床面と対向する面に押圧材を取付け、押圧材で発光部よりも指動脈血流のさらに下流側の指の先端部の上面を押さえるように構成することが好ましい。被験者は、脈波測定中に意識的・無意識的に指先に力を入れる場合がある。この場合に、力を抜くように指導すると、被験者の指の形状によっては、センサへの密着性が悪くなることがある。力を入れる場合も、力を抜く場合も、指の小さな動きによるノイズが発生する。押圧材を設けることにより、ノイズが減少し、測定データの再現性が増し、正確な脈波情報を得ることができる。
【００６３】
上記したように、押圧部位を指の先端部上面の小面積に限るように構成してあるので、押圧部位の上流側の指動脈血流部分の脈波を再現性よく正確に測定することができる。
上記発光部の側面は、発光部から指内に照射される光が外部に漏れるのを阻止しかつ指の腹部表面からの反射光を阻止するために、遮光壁で囲まれていることが好ましい。
【００６４】
上記発光部は、内面が照射光に対する反射特性を有する遮光壁内部に配置され、遮光壁の上端が、指の腹部を載置する床面から通常０．１〜１．５ｍｍ程度、好ましくは０．２〜１．０ｍｍ程度、さらに好ましくは０．３〜０．５ｍｍ程度突出するように構成され、指の腹部がこの上端に載置されて遮光壁上端全面を覆うようにすることが好ましい。遮光壁の上端をこのような範囲内になるように突出させると、指腹部の皮膚面が発光部上面から下に被さるので、測定データに対する外乱光や漏洩光や反射光の影響を小さくすることができると共に、被験者が指を装着する際に、遮光壁の突出位置に触れることによりセンサ位置を認識して、指を所定の位置に載置し易いという利点がある。
【００６５】
しかし、０．１ｍｍ未満であると、遮光壁位置を確認し難いので、指先を所定位置に載置し難いと共に、発光部からの照射光や指腹部表面からの反射光が漏れ易くなり、測定データに対する反射光の影響が大きくなる。また、１．５ｍｍを超えると、指皮膚面が床面から浮くため、不安定な装着状態になると共に、指を載置した時の指に対する圧力により波形の変形が生じて再現性が悪くなり、測定される脈波データにバラツキが生じて、正確な脈波情報を得ることが難しくなる。
上記押圧材で指の先端部の上面を押さえる際の圧力を通常５０〜２００ｇＷ、望ましくは７０〜１５０ｇＷになるように構成することが好ましい。
【００６６】
また、本発明で用いることができる脈波センサは、発光部と受光部とを有し、発光部と受光部とにより被験者の指の脈波を測定する反射型の脈波センサにおいて、発光部を、受光部より指の動脈血流の下流側に、その上面が受光部の上面よりも突出するように配置し、また、発光部及び受光部の上方に、赤外線透過性の窓部を、窓部上面が指の腹部を載置する床面よりも所定の距離だけ突出するように配置し、そして指の先端部を発光部よりも指動脈血流のさらに下流側に装着するための空間が床面の先端部分に設けられているものであっても良い。
【００６７】
上記窓部の上面は、指の腹部を載置する床面より通常０．１〜１．５ｍｍ程度、望ましくは０．２〜１．０ｍｍ程度、更に望ましくは０．３〜０．５ｍｍ程度突出するように配置されることが好ましい。
上記押圧材で窓部より先の指の先端部の上面を押さえる際の圧力を通常５０〜２００ｇＷ、望ましくは７０〜１５０ｇＷになるように構成することが好ましい。５０ｇＷ未満であると、脈波の振幅が小さく、測定時の振動や被験者の指のふるえによるノイズが入りやすく、また、被験者の緊張などにより測定波形が変形しやすい。また、２００ｇＷを超えると測定波形が極端に変形する。なお、指全体を加圧すると測定波形が変形し、再現性が低下する。
【００６８】
以下、図面を参照して本発明で用いることができる脈波センサについて図面を参照して詳細に説明する。
図１５は、反射型脈波センサの要部である指装着部の断面構造を示す図（図１５（ａ））、並びに発光部及び受光部の近傍を拡大して示す図（図１５（ｂ））であり、指を装着した状態で示す。
【００６９】
この反射型脈波センサは、発光部及び受光部を有し、被験者の指を装着して脈波の測定ができるものであって、蓋部を構成する開閉自在の合成樹脂製の上方部分と、指腹部を載置することができるように構成された合成樹脂製の床部分とを有してなる。上方部分はその内面が指の外形に沿った形状をしていてもよく、また、床部分はその床面が平坦であっても又は外乱光を遮断するために指の付け根側がやや高く又は低くなるように傾斜した形状になっていてもよい。以下述べるように、上方部分の先端部には押圧材が設けられ、床面に載置した指の動脈血流の発光部よりもさらに下流側の指の先端部分を押さえて固定できるようになっており、また、床部分には発光部及び受光部が所定の位置に配置されている。この押圧材は、指の先端部分を押さえて固定できるものであればよく、例えば、クッション材や、バネ材のような板材などでよく、所定の圧力で指の先端部分を押さえることができるように構成されている。また、このセンサには、反射光の電流／電圧変換回路、増幅器が設けられており、このセンサをパーソナルコンピュータに接続すれば、センサからの出力に基づき正確な脈波情報を得ることができる。
【００７０】
この脈波センサの場合、指装着部に指を差し込み、指の先端部分の腹部に赤外線などの光を当てると、毛細血管中のヘモグロビン（赤血球）が光の一部を吸収し、光の反射量が変化する（血液量が多い部分は光の反射量が少なくなる）。この微妙な光の反射量の変化を検出し、検出された反射光を電流から電圧へ変換し、増幅器に送信し、増幅された信号電圧をパーソナルコンピュータを利用してＡＤ変換して出力し、脈波情報として活用する。
【００７１】
脈波センサの要部を構成する指装着部には、図１５（ａ）及び（ｂ）に示すように、発光ダイオード（ＬＥＤ）などの半導体発光素子からなる発光部１が、フォトダイオード（ＰＤ）などの半導体受光素子からなる受光部２よりも被験者の指３の動脈血流の下流側に配置される。指内の発光部１からの照射光１ａの経路を見ると、発光部分の光束は、指内を進むにつれて拡散して拡がる。このため、発光部１からの入射光量変化による受光部２の光量変化は大きく、受光部２の位置変化による受光する拡散光の光量変化は小さい。従って、発光部１を指に密着するようにすることが必要となる。しかし、密着性を良くすることは、指に余分の圧力を加えることにつながる。そこで、本発明では、発光部１を受光部２より指動脈血流の下流側に配置し、指に余分の圧力が加わらないようにする。
【００７２】
また、発光部１は、その上面が受光部２の上面より所定の距離だけ突出する（すなわち、高くなる）ように配置される。すなわち、発光部１の高さＨ１が受光部２の高さＨ２より所定の距離だけ高くなるように構成する。指装着部の先端部分には、赤外ＬＥＤなどからなる発光部１から照射される光１ａの経路よりも指動脈血流のさらに下流側に空間４が設けられ、指３の先端部がこの空間内に載置され得るように構成される。
【００７３】
指装着部の指腹部を載置する面は、指載置床面５として構成される。床面５には発光部１及び受光部２が所定の位置に設けられ、床面の先端部分は傾斜して立ち上がり、指の先端が適切に納まるように構成される。この指装着部において、発光部１の配置された位置より動脈血流の下流側であって、床面と対向する面に押圧材６が設けられる。この押圧材により、装着された指の先端部分（好ましくは、爪部）を軽く押さえ、装着された指が動かないように固定する。このように構成することにより、被験者の意識的・無意識的な指の小さな動きが少なくなって、ノイズ発生が減少し、その結果、測定波形の変化も少なくなる。なお、押圧材により動脈血流の下流側を圧迫しても、脈波への影響は小さい。
【００７４】
発光部１からの照射光１ａが指の皮膚表面で反射して受光部２に入り込むと、この反射光がノイズとなり、受光部２に入る受光量が変動する。このために、正確な脈波を測定することができなくなる。また、照射光１ａが脈波センサの外部に漏れると、照射光の効率が下がり、かつ、受光部が受け取る反射光２ａの光量が減少して正確な脈波を測定することが困難になる。そのため、本発明では、余分な反射光や漏洩光を阻止するために、発光部１の側面を遮光壁７で囲むことが好ましい。
【００７５】
この遮光壁７は、反射光や漏洩光をなくすような形状であれば、その形状は問わないが、例えば、発光部１の外周形状に沿った円筒状などの形状が好ましい。装着された指は、この遮光壁の上面に７ａの点で密着し、固定される。遮光壁７は、その受光部２側が黒色であってもよく、また、その内面が鏡面であってもよい。遮光壁の材料としては、赤外線を遮る性質を有する材料であれば特に制限はなく、例えば、赤外線を実質的に透過させないポリプロピレン系樹脂やＡＢＳ系樹脂などの熱可塑性樹脂など、又はこれらに黒色塗装などの表面処理を施したものをあげることができる。
【００７６】
発光部１の上面に赤外線透過性の絶縁体キャップ８を設け、発光部１と指３とが直接接触しないようにする。これは、発光部の通電部に影響を及ぼさないようにするためと、発光部表面の汚れを清拭しないで済むようにするためである。絶縁体キャップ８の外形は、例えば、発光体１の上方部分の形状に沿った円筒形状などの形状であればよい。絶縁体キャップ８の上面を凹レンズで構成すれば、射出光の指向性をさらに広げることができる。この絶縁体キャップの材料としては、赤外線に対して透光性が高い赤外線透過性材料であれば特に制限はなく、例えば、アクリル樹脂、ポリエチレン樹脂、ポリカーボネート樹脂、ポリスチレン樹脂などをあげることができる。また、受光部２と指３とが直接接触して指に圧力がかからないように、受光部２と指３との間に隙間を設けるような構造とすることが好ましい。
【００７７】
発光部１の発光素子と受光部２の受光素子との指向性の影響について図１６に示す。図１６（ａ）に示すように、発光部１の発光素子と受光部２の受光素子との配置が従来の指向性が強い配置の場合、発光部１の発光ダイオードの光軸がずれると、有効な検出領域もずれる。しかし、本発明におけるように、指向性の弱い発光部１の発光素子と受光部２の受光素子とを接近して配置すれば、図１６（ｂ）に示すように、発光ダイオードの光軸のずれに対する有効検出領域のずれは相対的に小さい。そのため、得られる脈波情報は正確となる。
発光部１からの照射光の射出角（半値角）αを通常５０度以上、好ましくは５０〜８５度、より好ましくは５０〜８０度とすることにより、有効検出領域のずれは相対的に少なくなる。５０度未満であると、有効検出領域のずれが大きくなり、正確な脈波データを得ることが困難になる。
【００７８】
上記脈波センサにおいて、発光部１と受光部２との距離が長くなる程、加速度脈波の波形であるａ波の振幅は小さくなり、ノイズ成分が発生し易くなって、測定波形の変形が大きくなる傾向がある。また、この距離が長い程、圧力の影響がある指部位の脈波を測定することになり、測定波形が変形し易い。そのため、発光部と受光部との距離を、所定の距離、例えば、通常８ｍｍ以内、好ましくは６ｍｍ以内に設定すれば、加速度脈波のａ波の振幅及びｂ波とａ波との比（ｂ／ａ）は適切な範囲内に納まる。この場合、光軸のずれも少なく、有効検出領域のずれも少なく、また、波形は変形し難い。なお、この距離が上記範囲を外れた動脈上流側の指部位では、動脈血管が膨らんで、ｂ／ａが小さい（絶対値が大きい）状態になり、また、下流側の指部位では、鬱血状態となって、ｂ／ａが大きい（絶対値が小さい）状態になる。また、発光部と受光部との距離の下限は、特に制限はなく、発光部と受光部との物理的な大きさや脈波センサの大きさなどに依って所望により設定できる最低距離であればよい。例えば、２〜３ｍｍ程度に設定してもよい。
【００７９】
また、絶縁体キャップが脱落しないようにして、脈波センサ本体の取り扱い性を向上させるため、図１７に示すように、絶縁体キャップ１４の下方部分につば部１４ａを設けた構造としてもよい。図１７において、１１は発光部、１１ａは発光部からの照射光、１２は受光部、１３は遮光壁、１５は指載置床面を示す。発光部１１、受光部１２、遮光壁１３、床面１５などの配置位置関係については、図１５に示す場合と同様である。また、遮光壁１３及び絶縁体キャップ１４の材料も、図１５に示す遮光壁７及び絶縁体キャップ８の材料と同様である。さらに、絶縁体キャップ１４の上面を凹レンズで構成すれば、射出光の指向性をさらに広げることができる。
【００８０】
上記した脈波センサでは、受光部は、指に圧力がかからないように、その上面が指装着部の床面と同一高さ又はそれより下になるように配置される。これにより、受光部に入射する光の割合が最も多い受光部上面位置にあたる指部分が圧迫されないようになる。例えば、受光部を脈波センサの指載置床面より１ｍｍ程度低くなるように配置すればよい。
図１８は、さらに別の反射型脈波センサの要部である指装着部の断面構造を示す図（図１８（ａ））、並びに発光部及び受光部の近傍を拡大して示す図（図１８（ｂ））であり、指を装着した状態で示す。図１８において、図１５と同じ構成要素については同じ参照符号で示す。以下、図１５の場合と異なる構成要素について説明する。
【００８１】
図１８（ａ）及び（ｂ）に示す実施の形態によれば、発光部１及び受光部２の上方に赤外線透過性の窓部９が設けられている。この窓部９は、その上面が被験者の指の腹部を載置する床面５よりも所定の距離（０．１ｍｍ以上、例えば約０．３５ｍｍ）だけ高くなるように配置されている。窓部９を、図１８（ｂ）に示すように、発光部１及び受光部２を覆うように床面の縁端部上に載置して固定しても、又は縁端部の切欠き上に載置して固定してもよく、その配置の仕方に制限はない。窓部９を設けることにより、発光部１及び受光部２と被験者の指とが直接接触しないようになる。これにより、通電部に影響が及ぶことはなく、また、発光部及び受光部表面の汚れを清拭しないで済み、メンテナンスが楽になる。
【００８２】
窓部９の外形は特に制限されず、例えば、厚さ０．５ｍｍ程度の板状でよい。この窓部９の上面を凹レンズで構成すれば、射出光の指向性をさらに広げることができる。この窓部の材料としては、赤外線に対して透光性が高い赤外線透過性材料であれば特に制限はなく、例えば、アクリル樹脂、ポリエチレン樹脂、ポリカーボネート樹脂、ポリスチレン樹脂などをあげることができる。
【００８３】
上記脈波センサによれば、発光部を、受光部より指の動脈血流の下流側に、その上面が受光部の上面よりも所定の距離突出しかつ所定の距離だけ指の腹部を載置する床面よりも突出するように配置し、そして指の先端部が発光部よりも指動脈血流のさらに下流側に装着され得るように、床面の先端部分に空間を設けてあるので、指の発光部に対する密着性が良好になると共に、指動脈血流の下流側が圧迫されても、上流側が圧迫されるのに比べて脈波への影響は少なく、再現性よく正確な脈波情報を得ることができる。
また、発光部及び受光部の上方に赤外線透過性の窓部を配置する場合、窓部の上面が指の腹部を載置する床面よりも所定の距離だけ突出するように配置してあるので、測定波形への影響は少なく、再現性よく正確な脈波情報を得ることができる。
【００８４】
【発明の効果】
本発明によれば、所定の時間連続して計測した生体脈波の波形を２次微分して算出した加速度脈波の連続する波形からａ−ａ間隔を求めて、このａ−ａ間隔の変動が心電図のＲ−Ｒ間隔の変動に対応することから、このａ−ａ間隔の変動性の評価を心電図による心拍変動性評価と同等のものとして、自律神経機能評価などに用いることができるという効果を奏する。
本発明によれば、長時間連続測定や期外収縮波の検出も可能である。
【００８５】
また、これまでの心電計を用いる方法ではできなかった期外収縮補正後のａ−ａ間隔変動のパワースペクトルが算出できるので、臨床上応用可能の範囲が広がるという利点がある。 さらに、加速度脈波測定は着衣の着脱もなく、座位にて指先の脈波を測定するだけであるから、心電図の場合よりもより簡便である。
さらにまた、加速度脈波測定装置は心電計に比べて安価である。
従って、本発明によれば、糖尿病、神経障害、脳血管疾患、冠動脈疾患、喘息、更年期障害など多方面において、合併症リスクの予測や、治療効果判定、自己管理などにも幅広く用いることができる。
【図面の簡単な説明】
【図１】ａ、ｂ、ｃ、ｄ、ｅ波の５つの成分波をもつ加速度脈波の標準波形でり、心臓の収縮期の波形図。
【図２】脈波センサを用いて測定・記録した容積脈波、速度脈波（一次微分波）、加速度脈波（二次微分波）の波形図。
【図３】心電図（ＥＣＧ）と加速度脈波（ＡＰＧ）とを同時測定した波形を示す波形図。
【図４】同時測定した加速度脈波のａ−ａ間隔と心電図のＲ−Ｒ間隔との関係を示すグラフ。
【図５】加速度脈波解析のフロー図。
【図６】加速度脈波の波高差の積分（差積分）による類似度判定について説明するためのグラフ。
【図７】ａ−ａ間隔再判定処理の流れを説明するためのグラフ。
【図８】ａ−ａ間隔の異常値判定を説明するための、拍数に対するａ−ａ間隔（ｍｓｅｃ）の関係を示すグラフ。
【図９】Ｔａａの分布全体から正常範囲を決定すると異常ａ−ａ間隔が検出されないケースはあることを説明するための、拍数に対するａ−ａ間隔（ｍｓｅｃ）の関係を示すグラフ。
【図１０】拍数に対するａ−ａ間隔（ｍｓｅｃ）の関係を示すグラフであって、単純推移平均ではうまくいかないケースを説明するためのグラフ。
【図１１】ａ−ａ間隔の異常値判定に関し、推移平均と標準ａ−ａ間隔とを比較して、拍数に対するａ−ａ間隔（ｍｓｅｃ）の関係を示すグラフ。
【図１２】健康人についての加速度脈波のａ−ａ間隔の変動係数（ＣＶａａ％）の平均値と標準偏差との関係を示すグラフ。
【図１３】図１２に示す各年代における標準的な変動係数の範囲を、［平均値−標準偏差］〜［平均値＋標準偏差］として示すグラフ。
【図１４】糖尿病患者の加速度脈波のａ−ａ間隔の変動係数を健康人の場合と比較して示すグラフ。
【図１５】本発明で用いることができる脈波センサの構造を示す断面図であり、（ａ）はこの脈波センサの要部である指装着部の模式的構造を示す断面図、（ｂ）は図１５（ａ）の発光部及び受光部の近傍の拡大断面図。
【図１６】光の指向性を示す模式図であり、（ａ）は従来技術における発光素子と受光素子との指向性の影響について示す図、（ｂ）は本発明における発光素子と受光素子との指向性の影響について示す図。
【図１７】本発明で用いることができる別の脈波センサにおいてつば部を有する絶縁体キャップの周辺の模式的構造を示す断面図。
【図１８】本発明で用いることがができる更に別の脈波センサの構造を示す断面図であり、（ａ）はこの脈波センサの要部である指装着部の模式的構造を示す断面図、（ｂ）は（ａ）の発光部及び受光部の近傍の拡大断面図。
【符号の説明】
１ 発光部 ２ 受光部
２ａ 反射光 Ｈ１ 発光部の高さ
Ｈ２ 受光部の高さ ３ 指
４ 空間 ５ 指載置床面
６ 押圧材 ７ 遮光壁
７ａ 指と遮光壁との密着点 ９ 窓部
１５ 指載置床面 α 照射光の射出角（半値角）
Claims (22)
- 生体の脈波から心電図のＲ−Ｒ間隔の変動に対応する情報を求める脈波解析法であって、所定の時間連続して計測した脈波の波形を２次微分して加速度脈波を算出し、この加速度脈波の連続する波形からａ−ａ間隔を求め、このａ−ａ間隔の変動を該心電図のＲ−Ｒ間隔の変動に対応する間隔とすることを特徴とする脈波解析法。
- 生体の脈波から心電図のＲ−Ｒ間隔の変動に対応する情報を求める脈波解析法であって、過去の加速度脈波から求めた基準となる波形を標準加速度脈波とし、この標準加速度脈波と測定加速度脈波との波形の類似度を評価して個別加速度脈波を決定し、この個別加速度脈波の連続する波形からａ−ａ間隔を求め、このａ−ａ間隔の変動を該心電図のＲ−Ｒ間隔の変動に対応する間隔とすることを特徴とする脈波解析法。
- 請求項２において、波形の類似度の評価を、測定加速度脈波のａ波候補の波高値が標準加速度脈波のａ波の波高値と同一若しくは両者の中間値になるように規格化した波形を用いて行い、また、標準加速度脈波のａ波と測定加速度脈波のａ波候補との時間位置を合わせ、規格化した測定加速度脈波と標準加速度脈波との波高差を求めて行うことを特徴とする脈波解析法。
- 請求項２又は３において、波形の類似度の評価を、個別加速度脈波のａ波判定条件として、測定加速度脈波をそのａ波候補の波高値と標準加速度脈波のａ波波高値との中間値で若しくは該ａ波波高値で規格化した個別加速度脈波と標準加速度脈波との波高差の積分値又は多重積分値を用いて算出した類似度に基づいて行うことを特徴とする脈波解析法。
- 請求項１〜４のいずれかにおいて、隣り合うａ−ａ間隔の変化が一定範囲を外れるものを除外したａ−ａ間隔について、これらの値により標準ａ−ａ間隔を構成して、ａ−ａ間隔の異常を判定することを特徴とする脈波解析法。
- 請求項５において、標準ａ−ａ間隔として、過去のデータの平均若しくは重み付き平均を用い、異常判定がでたデータに対して、未来データを含む平均若しくは重み付き平均で再判定することを特徴とする脈波解析法。
- 請求項１〜４のいずれかにおいて、評価している時点のａ−ａ間隔が、過去の平均的なａ−ａ間隔、又は未来よりも過去のａ−ａ間隔データの重みを大きくした重み付き推移平均、又は直前のａ−ａ間隔を基準として、一定の基準値を超えるとき、その時点の直後から順次、隣り合うａ−ａ間隔の変化が一定範囲内に入るものを数点選び出し、これらの値により標準ａ−ａ間隔を更新して、ａ−ａ間隔の異常を判定することを特徴とする脈波解析法。
- 請求項１〜７のいずれかにおいて、生体の脈波を、発光部と受光部とを有する反射型脈波センサであって、該発光部が、該受光部より指の動脈血流の下流側に、その上面が該受光部の上面よりも突出しかつ所定の距離だけ指の腹部を載置する床面よりも突出するように配置され、そして指の先端部を該発光部よりも指動脈血流のさらに下流側に装着するための空間が該床面の先端部分に設けられている脈波センサに指先を装着して測定することを特徴とする脈波解析法。
- 請求項８において、該脈波センサの発光部の上面が指の腹部を載置する床面より０．１〜１．５ｍｍ突出するように、発光部が配置されている脈波センサに指先を装着して脈波を測定することを特徴とする脈波解析法。
- 請求項８又は９において、該脈波センサの受光部の上面が発光部の上面と同一レベルか若しくはそれより下になるように、また、所定の距離だけ該床面より低くなるように、受光部が配置されている脈波センサに指先を装着して脈波を測定することを特徴とする脈波解析法。
- 請求項８〜１０のいずれかにおいて、該脈波センサの指の先端部を装着する空間の床面と対向する面に押圧材が取付けられ、該押圧材で発光部よりも指動脈血流のさらに下流側の指の先端部の上面を押さえるように構成されている脈波センサに指先を装着して脈波を測定することを特徴とする脈波解析法。
- 請求項１１において、該押圧材は、この押圧材で発光部よりも指動脈血流のさらに下流側の指の先端部の上面を５０〜２００ｇＷの圧力で押さえるように構成されている脈波センサに指先を装着して脈波を測定することを特徴とする脈波解析法。
- 請求項８〜１２のいずれかにおいて、該脈波センサの発光部の側面は、内面が照射光に対する反射特性を有する遮光壁で囲まれ、そして該遮光壁上端が指の腹部を載置する床面から０．１〜１．５ｍｍ突出するように構成されている脈波センサに指先を装着して脈波を測定することを特徴とする脈波解析法。
- 請求項８〜１３のいずれかにおいて、該脈波センサの発光部及び受光部の上方に、赤外線透過性の窓部を、窓部上面が指の腹部を載置する床面よりも所定の距離だけ突出するように配置してある脈波センサに指先を装着して測定することを特徴とする脈波解析法。
- 請求項１４において、該窓部の上面が指の腹部を載置する床面より０．１〜０．５ｍｍ突出するように配置されている脈波センサに指先を装着して測定することを特徴とする脈波解析法。
- 請求項１〜１５のいずれかに記載の脈波解析法に基づき、連続するａ−ａ間隔について、隣り合うａ−ａ間隔の最初のａ−ａ間隔が標準ａ−ａ間隔より正常範囲を外れて短い場合で、かつ、前回のａ−ａ間隔と今回のａ−ａ間隔との和が標準ａ−ａ間隔の１倍若しくは２倍を基準にしてａ−ａ間隔の変動の分布から定まる所定の範囲に入るとき、期外収縮と判定することを特徴とする自律神経機能評価法。
- ａ−ａ間隔の統計量に基づいて自律神経機能を評価する方法であって、請求項１〜１５のいずれかに記載の脈波解析法に基づき、ａ−ａ間隔が正常範囲を外れる場合はこれを除外してａ−ａ間隔の統計量を算出し、これに基づいて自律神経機能を評価することを特徴とする自律神経機能評価法。
- ａ−ａ間隔の時間変動を周波数解析して自律神経機能を評価する方法であって、隣り合うａ−ａ間隔について、最初のａ−ａ間隔が、請求項１〜１５のいずれかに記載の脈波解析法に基づき、標準ａ−ａ間隔より正常範囲を外れて短い場合で、かつ、最初のａ−ａ間隔と次のａ−ａ間隔との和が、標準ａ−ａ間隔の１倍若しくは２倍を基準にして所定の範囲に入る時、該１倍を基準とする場合には、最初のａ−ａ間隔を最初のａ−ａ間隔と次のａ−ａ間隔との両者の和に補正して次のａ−ａ間隔を削除し、また、該２倍を基準とする場合には、それぞれ二つのａ−ａ間隔を該両者の和の二分の一に補正して、周波数解析することにより自律神経機能を評価することを特徴とする自律神経機能評価法。
- 請求項１６〜１８のいずれかにおいて、ａ−ａ間隔の変動性を各年代別標準参考値を参照しながら評価することにより、自律神経機能を評価することを特徴とする自律神経機能評価法。
- 請求項１〜１５のいずれかに記載の脈波解析法に基づき、長時間連続測定や期外収縮波の検出を行うことにより、また、期外収縮除去後の推定変動係数を算出することにより、糖尿病、神経障害、脳血管疾患、冠動脈疾患、喘息、更年期障害などを含めて、合併症リスクの予測や、治療効果判定、自己健康管理を行うことを特徴とする自律神経機能評価法。
- 生体の脈波を検出して脈波の大きさに応じた信号を出力する脈波計測手段と、該脈波計測手段により得られた脈波の波形を２次微分して算出された加速度脈波の波形からその波形パラメータを解析する波形パラメータ解析手段とを備え、該波形パラメータ解析手段が、所定の時間連続して計測した脈波の波形を２次微分して加速度脈波を算出する手段と、この加速度脈波の連続する波形から心電図のＲ−Ｒ間隔の変動に対応するａ−ａ間隔の変動を求める手段とを有することを特徴とする加速度脈波の変動解析による自律神経機能評価装置。
- 請求項２１において、さらに、隣り合うａ−ａ間隔の変化が一定範囲を外れるものを除外したａ−ａ間隔について、これらの値により標準ａ−ａ間隔を構成して、ａ−ａ間隔の異常を判定する手段を有することを特徴とする自律神経機能評価装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003161709A JP4625886B2 (ja) | 2003-06-06 | 2003-06-06 | 脈波解析法及び自律神経機能評価装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003161709A JP4625886B2 (ja) | 2003-06-06 | 2003-06-06 | 脈波解析法及び自律神経機能評価装置 |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2009278051A Division JP5080550B2 (ja) | 2009-12-07 | 2009-12-07 | 自律神経機能評価装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2004358022A true JP2004358022A (ja) | 2004-12-24 |
| JP4625886B2 JP4625886B2 (ja) | 2011-02-02 |
Family
ID=34054051
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2003161709A Expired - Fee Related JP4625886B2 (ja) | 2003-06-06 | 2003-06-06 | 脈波解析法及び自律神経機能評価装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4625886B2 (ja) |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006271731A (ja) * | 2005-03-29 | 2006-10-12 | Toshiba Corp | 心拍計測装置及び心拍計測方法 |
| JP2007125366A (ja) * | 2005-10-06 | 2007-05-24 | Konica Minolta Sensing Inc | 脈波データ解析方法、システム、プログラム |
| JP2007203041A (ja) * | 2006-01-31 | 2007-08-16 | Up Management Gmbh & Co Med-Systems Kg | 心肺相互作用を用いる患者の血液動態状態の評価用装置 |
| JP2007236807A (ja) * | 2006-03-10 | 2007-09-20 | Haruko Takada | 運動耐容能評価装置 |
| JP2008094787A (ja) * | 2006-10-13 | 2008-04-24 | Paq Kk | 落花生処理物を用いた交感神経興奮剤 |
| JP2008237574A (ja) * | 2007-03-27 | 2008-10-09 | Toshiba Corp | 自律神経指標を計測する装置、方法および生体情報を検出する装置 |
| EP2085025A1 (en) | 2008-01-31 | 2009-08-05 | Nihon Kohden Corporation | Autonomic nerve activity measuring apparatus and autonomic nerve activity measuring method |
| JP2009189416A (ja) * | 2008-02-12 | 2009-08-27 | Toshiba Corp | 脈波間隔計測装置及び計測方法 |
| JP2010162069A (ja) * | 2009-01-13 | 2010-07-29 | Nihon Univ | 非接触生体情報収集による年齢推定装置、年齢推定方法及び年齢推定プログラム |
| WO2010103817A1 (ja) * | 2009-03-10 | 2010-09-16 | 株式会社カオテック研究所 | 自律神経バランス演算装置またはその方法 |
| JP2012147879A (ja) * | 2011-01-18 | 2012-08-09 | Hiro Kagaku Kenkyusho:Kk | 自律神経機能年齢の判定システム及び判定方法 |
| JP2015097932A (ja) * | 2009-05-12 | 2015-05-28 | アンジオロジックス インコーポレイテッド | 肢部分の動脈容積における変化の計測システムおよび記憶媒体 |
| JP2015215351A (ja) * | 2014-05-09 | 2015-12-03 | 三星電子株式会社Ｓａｍｓｕｎｇ Ｅｌｅｃｔｒｏｎｉｃｓ Ｃｏ．，Ｌｔｄ． | 分光センサ、及びそれを採用した分光器 |
| US9289177B2 (en) | 2011-01-20 | 2016-03-22 | Nitto Denko Corporation | Sensing device, a method of preparing a sensing device and a personal mobile sensing system |
| JP2016165333A (ja) * | 2015-03-09 | 2016-09-15 | 株式会社デンソー | 脈波検出装置 |
| WO2018194205A1 (ko) * | 2017-04-21 | 2018-10-25 | 조선대학교 산학협력단 | 광전용적맥파를 이용한 뇌혈관 질환 및 협착 분석 장치 및 방법 |
| JP2020008590A (ja) * | 2014-05-09 | 2020-01-16 | 三星電子株式会社Ｓａｍｓｕｎｇ Ｅｌｅｃｔｒｏｎｉｃｓ Ｃｏ．，Ｌｔｄ． | 分光センサ、及びそれを採用した分光器 |
| WO2021125387A1 (ko) * | 2019-12-17 | 2021-06-24 | 전자부품연구원 | 당뇨 추정 시스템 및 그 방법 |
| CN113100727A (zh) * | 2021-05-12 | 2021-07-13 | 深圳市通久电子有限公司 | 实时分析识别脉搏波峰的方法 |
| CN113598724A (zh) * | 2021-09-03 | 2021-11-05 | 上海市高血压研究所 | 用于评估心血管自主神经功能的装置 |
| JP2021536307A (ja) * | 2018-09-06 | 2021-12-27 | ヴァンダービルト ユニバーシティ | 被験体を評価するための非侵襲的静脈波形分析 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5113692U (ja) * | 1974-07-18 | 1976-01-31 | ||
| JPS56143140A (en) * | 1980-04-04 | 1981-11-07 | Matsushita Electric Works Ltd | Pulse detecting sensor |
| JPS5828603U (ja) * | 1981-08-14 | 1983-02-24 | 松下電工株式会社 | 脈拍センサ |
| JPS62145601U (ja) * | 1986-03-11 | 1987-09-14 | ||
| JPH1043150A (ja) * | 1996-08-06 | 1998-02-17 | Omron Corp | 脈波検出装置 |
| JP2003290162A (ja) * | 2002-04-03 | 2003-10-14 | U-Medica Inc | 脈波計測装置及びノイズ成分除去方法 |
| JP2004000316A (ja) * | 2002-04-26 | 2004-01-08 | U-Medica Inc | 脈波計測・解析装置及び脈波の平均波形算出・提示方法 |
| JP2004000467A (ja) * | 2002-03-15 | 2004-01-08 | U-Medica Inc | 脈波センサ |
-
2003
- 2003-06-06 JP JP2003161709A patent/JP4625886B2/ja not_active Expired - Fee Related
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5113692U (ja) * | 1974-07-18 | 1976-01-31 | ||
| JPS56143140A (en) * | 1980-04-04 | 1981-11-07 | Matsushita Electric Works Ltd | Pulse detecting sensor |
| JPS5828603U (ja) * | 1981-08-14 | 1983-02-24 | 松下電工株式会社 | 脈拍センサ |
| JPS62145601U (ja) * | 1986-03-11 | 1987-09-14 | ||
| JPH1043150A (ja) * | 1996-08-06 | 1998-02-17 | Omron Corp | 脈波検出装置 |
| JP2004000467A (ja) * | 2002-03-15 | 2004-01-08 | U-Medica Inc | 脈波センサ |
| JP2003290162A (ja) * | 2002-04-03 | 2003-10-14 | U-Medica Inc | 脈波計測装置及びノイズ成分除去方法 |
| JP2004000316A (ja) * | 2002-04-26 | 2004-01-08 | U-Medica Inc | 脈波計測・解析装置及び脈波の平均波形算出・提示方法 |
Non-Patent Citations (1)
| Title |
|---|
| 高田晴子、鷲見尚美: "「加速度脈波ａ−ａ間隔変動による自律神経評価法の検討」", 教育医学, vol. 第４８巻第１号, JPN6009033036, 8 July 2002 (2002-07-08), JP, pages 166 - 167, ISSN: 0001360532 * |
Cited By (32)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006271731A (ja) * | 2005-03-29 | 2006-10-12 | Toshiba Corp | 心拍計測装置及び心拍計測方法 |
| JP2007125366A (ja) * | 2005-10-06 | 2007-05-24 | Konica Minolta Sensing Inc | 脈波データ解析方法、システム、プログラム |
| JP2007203041A (ja) * | 2006-01-31 | 2007-08-16 | Up Management Gmbh & Co Med-Systems Kg | 心肺相互作用を用いる患者の血液動態状態の評価用装置 |
| JP2007236807A (ja) * | 2006-03-10 | 2007-09-20 | Haruko Takada | 運動耐容能評価装置 |
| JP2008094787A (ja) * | 2006-10-13 | 2008-04-24 | Paq Kk | 落花生処理物を用いた交感神経興奮剤 |
| US8157732B2 (en) | 2007-03-27 | 2012-04-17 | Kabushiki Kaisha Toshiba | Method and apparatus for measuring autonomic-nervous index and apparatus for detecting biological information |
| JP2008237574A (ja) * | 2007-03-27 | 2008-10-09 | Toshiba Corp | 自律神経指標を計測する装置、方法および生体情報を検出する装置 |
| EP2085025A1 (en) | 2008-01-31 | 2009-08-05 | Nihon Kohden Corporation | Autonomic nerve activity measuring apparatus and autonomic nerve activity measuring method |
| US8382679B2 (en) | 2008-01-31 | 2013-02-26 | Nihon Kohden Corporation | Autonomic nerve activity measuring apparatus and autonomic nerve activity measuring method |
| JP2009189416A (ja) * | 2008-02-12 | 2009-08-27 | Toshiba Corp | 脈波間隔計測装置及び計測方法 |
| US8608663B2 (en) | 2008-02-12 | 2013-12-17 | Kabushiki Kaisha Toshiba | Pulse wave period measuring apparatus, pulse wave period measuring method, heartbeat period measuring method, and recording medium |
| JP2010162069A (ja) * | 2009-01-13 | 2010-07-29 | Nihon Univ | 非接触生体情報収集による年齢推定装置、年齢推定方法及び年齢推定プログラム |
| WO2010103817A1 (ja) * | 2009-03-10 | 2010-09-16 | 株式会社カオテック研究所 | 自律神経バランス演算装置またはその方法 |
| CN102209493A (zh) * | 2009-03-10 | 2011-10-05 | 混沌技术研究所 | 植物神经平衡运算装置或其方法 |
| JP5566997B2 (ja) * | 2009-03-10 | 2014-08-06 | 株式会社カオテック研究所 | 自律神経バランス演算装置およびその方法 |
| JP2015097932A (ja) * | 2009-05-12 | 2015-05-28 | アンジオロジックス インコーポレイテッド | 肢部分の動脈容積における変化の計測システムおよび記憶媒体 |
| JP2012147879A (ja) * | 2011-01-18 | 2012-08-09 | Hiro Kagaku Kenkyusho:Kk | 自律神経機能年齢の判定システム及び判定方法 |
| US9289177B2 (en) | 2011-01-20 | 2016-03-22 | Nitto Denko Corporation | Sensing device, a method of preparing a sensing device and a personal mobile sensing system |
| JP2015215351A (ja) * | 2014-05-09 | 2015-12-03 | 三星電子株式会社Ｓａｍｓｕｎｇ Ｅｌｅｃｔｒｏｎｉｃｓ Ｃｏ．，Ｌｔｄ． | 分光センサ、及びそれを採用した分光器 |
| JP2020008590A (ja) * | 2014-05-09 | 2020-01-16 | 三星電子株式会社Ｓａｍｓｕｎｇ Ｅｌｅｃｔｒｏｎｉｃｓ Ｃｏ．，Ｌｔｄ． | 分光センサ、及びそれを採用した分光器 |
| JP2016165333A (ja) * | 2015-03-09 | 2016-09-15 | 株式会社デンソー | 脈波検出装置 |
| WO2018194205A1 (ko) * | 2017-04-21 | 2018-10-25 | 조선대학교 산학협력단 | 광전용적맥파를 이용한 뇌혈관 질환 및 협착 분석 장치 및 방법 |
| KR20180118472A (ko) * | 2017-04-21 | 2018-10-31 | 조선대학교산학협력단 | 광전용적맥파를 이용한 뇌혈관 질환 및 협착 분석 장치 및 방법 |
| KR101994788B1 (ko) * | 2017-04-21 | 2019-07-01 | 조선대학교 산학협력단 | 광전용적맥파를 이용한 뇌혈관 질환 및 협착 분석 장치 및 방법 |
| US11259711B2 (en) | 2017-04-21 | 2022-03-01 | Industry-Academic Cooperation Foundation Chosun University | Device and method for analyzing cerebrovascular disease and stenosis by using photoplethysmography |
| JP2021536307A (ja) * | 2018-09-06 | 2021-12-27 | ヴァンダービルト ユニバーシティ | 被験体を評価するための非侵襲的静脈波形分析 |
| KR20210077246A (ko) * | 2019-12-17 | 2021-06-25 | 한국전자기술연구원 | 당뇨 추정 시스템 및 그 방법 |
| WO2021125387A1 (ko) * | 2019-12-17 | 2021-06-24 | 전자부품연구원 | 당뇨 추정 시스템 및 그 방법 |
| KR102371699B1 (ko) | 2019-12-17 | 2022-03-08 | 한국전자기술연구원 | 당뇨 추정 시스템 및 그 방법 |
| CN113100727A (zh) * | 2021-05-12 | 2021-07-13 | 深圳市通久电子有限公司 | 实时分析识别脉搏波峰的方法 |
| CN113100727B (zh) * | 2021-05-12 | 2023-09-19 | 深圳市通久电子有限公司 | 实时分析识别脉搏波峰的方法 |
| CN113598724A (zh) * | 2021-09-03 | 2021-11-05 | 上海市高血压研究所 | 用于评估心血管自主神经功能的装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP4625886B2 (ja) | 2011-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4625886B2 (ja) | 脈波解析法及び自律神経機能評価装置 | |
| Charlton et al. | Wearable photoplethysmography for cardiovascular monitoring | |
| US20060224073A1 (en) | Integrated physiological signal assessing device | |
| US20210244302A1 (en) | Methods to estimate the blood pressure and the arterial stiffness based on photoplethysmographic (ppg) signals | |
| KR100552681B1 (ko) | 수면 무호흡 진단 장치 및 방법 | |
| US11160459B2 (en) | Monitoring health status of people suffering from chronic diseases | |
| US6293915B1 (en) | Pulse wave examination apparatus, blood pressure monitor, pulse waveform monitor, and pharmacological action monitor | |
| US7044917B2 (en) | Method and apparatus for measuring pulsus paradoxus | |
| JP5080550B2 (ja) | 自律神経機能評価装置 | |
| US20230060335A1 (en) | Blood pressure and autoregulation monitoring | |
| JP2006158974A (ja) | 一体型生理学的信号評価装置 | |
| US7628760B2 (en) | Circulation monitoring system and method | |
| Peltokangas et al. | Monitoring arterial pulse waves with synchronous body sensor network | |
| TW200425873A (en) | System for measurement and analysis of vasodilatation | |
| Scardulla et al. | Photoplethysmograhic sensors, potential and limitations: Is it time for regulation? A comprehensive review | |
| US10327648B2 (en) | Blood vessel mechanical signal analysis | |
| KR20190105421A (ko) | 광혈류 측정기 기반의 인체착용형 혈압 측정장치 및 혈압 측정방법 | |
| JP4347338B2 (ja) | 血管年齢の評価方法 | |
| CN210697620U (zh) | 一种全自动穿戴式心血管健康诊断手环 | |
| JP7133576B2 (ja) | 位相差法による連続血圧測定システム | |
| WO2021249850A1 (en) | Wearable device | |
| Freithaler et al. | Smartphone-Based Blood Pressure Monitoring via the Oscillometric Finger Pressing Method: Analysis of Oscillation Width Variations Can Improve Diastolic Pressure Computation | |
| Huotari et al. | Arterial pulse wave analysis based on PPG and EMFi measurements | |
| WO2023162756A1 (ja) | 血行動態推定方法 | |
| 양승만 | Non-invasive Continuous Arterial Blood Pressure Monitoring System: Large Population Models and Single Chest-worn Device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060605 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20070830 |
|
| RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20070830 |
|
| RD05 | Notification of revocation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7425 Effective date: 20071001 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20071005 |
|
| RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20071005 |
|
| RD05 | Notification of revocation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7425 Effective date: 20071005 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20071005 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090701 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090831 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20091007 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20091207 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100113 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100315 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100407 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100409 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20100818 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20100915 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20101005 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131119 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 4625886 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |