(57)【要約】【課題】脈波信号を処理して血圧を精度良く演算する。【解決手段】脈波信号処理装置は、脈波信号に基づく特徴量を抽出する特徴量抽出部と、前記特徴量を保持する特徴量保持部と、前記特徴量に基づき血圧を演算する血圧演算部とを備える。前記血圧演算部は、過去に前記特徴量保持部に保持された前記特徴量から算出される値と、前記特徴量抽出部によって抽出された現在の前記特徴量から算出される値とから血圧を演算する。【選択図】図１
2014008200.tif 000002
【特許請求の範囲】
【請求項１】
脈波信号に基づく特徴量を抽出する特徴量抽出部と、
前記特徴量を保持する特徴量保持部と、
前記特徴量に基づき血圧を演算する血圧演算部とを備え、
前記血圧演算部は、過去に前記特徴量保持部に保持された前記特徴量から算出される値と、前記特徴量抽出部によって抽出された現在の前記特徴量から算出される値とから血圧を演算する脈波信号処理装置。
【請求項２】
前記特徴量は、前記脈波信号を２回微分した加速度脈波信号に関する特徴量を含み、
前記過去に前記特徴量保持部に保持された前記特徴量から算出される値は、前記加速度脈波信号の極値を用いて算出される請求項１に記載の脈波信号処理装置。
【請求項３】
前記過去に前記特徴量保持部に保持された前記特徴量から算出される値は、前記脈波信号が１周期の最小値になるタイミングの直前の前記加速度脈波信号の極小値、前記脈波信号が１周期の最小値になるタイミングの直前の前記加速度脈波信号の極大値、及び前記脈波信号が１周期の最小値になるタイミングの直後の前記加速度脈波信号の極大値を用いて算出される請求項２に記載の脈波信号処理装置。
【請求項４】
前記特徴量保持部は、生体毎に前記特徴量を保持する請求項１〜３のいずれか１項に記載の脈波信号処理装置。
【請求項５】
生体の血流循環により生じる脈波を検出する脈波検出部と、
前記脈波検出部から出力される脈波信号を処理する請求項１〜４のいずれか１項に記載の脈波信号処理装置とを備える血圧測定装置。
【請求項６】
コンピュータを、
脈波信号に基づく特徴量を抽出する特徴量抽出部、
前記特徴量を保持する特徴量保持部、及び
前記特徴量に基づき血圧を演算する血圧演算部として機能させるための脈波信号処理プログラムであって、
前記血圧演算部は、過去に前記特徴量保持部に保持された前記特徴量から算出される値と、前記特徴量抽出部によって抽出された現在の前記特徴量から算出される値とから血圧を演算する脈波信号処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、脈波信号を処理して血圧を演算する脈波信号処理装置及びそれを備える血圧測定装置及び脈波信号処理プログラムに関する。
【背景技術】
【０００２】
従来より、生体の血流循環により生じる脈波を用いて血圧を測定する技術が提案されている（例えば特許文献１〜３参照）。
【０００３】
特許文献１には、脈波に含まれる進行波成分と反射波成分を抽出することで脈波伝播時間を算出し、脈波伝播時間から血圧を算出する方法が開示されている。特許文献１で開示されている血圧の算出式は、次の通りである。
血圧＝α×脈波伝播時間＋β
【０００４】
なお、特許文献１には、係数α、βは固定値、もしくは、他の方法（カフによる測定や観血的な測定）で得られた実測値と、対応する時点の脈波伝播時間とを用いて都度最適な値となるよう更新するように構成することも可能であると記載されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００７−７０５５号公報
【特許文献２】特開平１０−２９５６５７号公報
【特許文献３】特開平１１−１５５８２６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、α、βが固定値である場合は、血圧の測定対象である生体毎に係数α、βの最適値が異なるため、血圧を精度良く算出することはできない。また、係数α、βの値を最適値に更新するには、他の方法で血圧を測定しなければならない。
【０００７】
図９Ａ及び図９Ｂは、複数の被験者に対して、本発明者が試行錯誤して係数α、βを最適な固定値になるように調整した場合の、脈波伝播時間の一次関数式から算出した血圧とカフを用いる血圧計で計測した血圧との相関を示す図である。図９Ａは係数α、βをそれぞれ−０．２、１６４とした場合の最高血圧に関する相関図であり、図９Ｂは係数α、βをそれぞれ−０．２、１２３とした場合の最低血圧に関する相関図である。図９Ａ及び図９Ｂから、測定対象である被験者が複数である場合、被験者毎に係数α、βの最適値が異なるため、脈波伝播時間の一次関数式から算出した血圧とカフを用いる血圧計で計測した血圧との相関は低く、脈波伝播時間の一次関数式から算出した血圧は精度が良くないことが分かる。
【０００８】
また、特許文献２及び３では、脈波信号から一拍毎の脈波波形を複数個抽出して基線を合わせて平均し、平均的な脈波信号を求める工夫や部位毎に演算された血圧の平均値を演算する工夫等が行われているが、血圧を精度良く算出する算出式に関する開示はない。
【０００９】
本発明は、上記の状況に鑑み、脈波信号を処理して血圧を精度良く演算する脈波信号処理装置及びそれを備える血圧測定装置及び脈波信号処理プログラムを提供することを目的とする。
【課題を解決するための手段】
【００１０】
上記目的を達成するために本発明に係る脈波信号処理装置は、脈波信号に基づく特徴量を抽出する特徴量抽出部と、前記特徴量を保持する特徴量保持部と、前記特徴量に基づき血圧を演算する血圧演算部とを備え、前記血圧演算部は、過去に前記特徴量保持部に保持された前記特徴量から算出される値と、前記特徴量抽出部によって抽出された現在の前記特徴量から算出される値とから血圧を演算する。なお、本明細書において、「現在の前記特徴量」とは、最新の脈波測定によって得られる脈波信号に基づいて抽出される特徴量を意味している。また、複数の生体の血圧演算が可能な仕様にする場合には、最新の脈波測定に関して、複数の生体全部の範囲で「最新」を定義するのではなく、生体毎の範囲で「最新」を定義する。すなわち、ｎ（ｎは２以上の自然数）個の生体それぞれに対して少なくとも一回は脈波測定が行われた場合、最新の脈波測定はｎ個存在することになり、各最新の脈波測定が各生体に一対一対応する。
【００１１】
上記目的を達成するために本発明に係る血圧測定装置は、生体の血流循環により生じる脈波を検出する脈波検出部と、前記脈波検出部から出力される脈波信号を処理する上記構成の脈波信号処理装置とを備える。
【００１２】
上記目的を達成するための脈波信号処理方法は、脈波信号に基づく特徴量を抽出する特徴量抽出ステップと、前記特徴量を保持する特徴量保持ステップと、前記特徴量に基づき血圧を演算する血圧演算ステップとを備え、前記血圧演算ステップは、過去に前記特徴量保持ステップにおいて保持された前記特徴量から算出される値と、前記特徴量抽出ステップにおいて抽出された現在の前記特徴量から算出される値とから血圧を演算するステップである。
【００１３】
上記目的を達成するために本発明に係る脈波信号処理プログラムは、コンピュータを、脈波信号に基づく特徴量を抽出する特徴量抽出部、前記特徴量を保持する特徴量保持部、及び前記特徴量に基づき血圧を演算する血圧演算部として機能させるための脈波信号処理プログラムであって、前記血圧演算部は、過去に前記特徴量保持部に保持された前記特徴量から算出される値と、前記特徴量抽出部によって抽出された現在の前記特徴量から算出される値とから血圧を演算する。
【発明の効果】
【００１４】
本発明によると、過去に保持された特徴量から算出される値と、現在の特徴量から算出される値とから血圧を演算するので、現在の特徴量のみから血圧を演算する従来の技術に比べて、血圧を精度良く演算することができる。
【図面の簡単な説明】
【００１５】
【図１】本発明の一実施形態に係る血圧測定装置を示すブロック図である。
【図２】脈波信号、速度脈波信号、及び加速度脈波信号の各波形例を示すタイムチャートである。
【図３Ａ】標準血圧に関連する指標の平均値とカフを用いる血圧計で計測した最高血圧の平均値との相関を示す図である。
【図３Ｂ】標準血圧に関連する指標の平均値とカフを用いる血圧計で計測した最低圧の平均値との相関を示す図である。
【図４Ａ】血圧演算部が算出した現在の最高血圧とカフを用いる血圧計で計測した最高血圧との相関を示す図である。
【図４Ｂ】血圧演算部が算出した現在の最低血圧とカフを用いる血圧計で計測した最低血圧との相関を示す図である。
【図５】図１に示す血圧測定装置が実施する血圧測定処理方法の一例を示すフローチャートである。
【図６】スマートフォンの一例を示す斜視図である。
【図７】スマートフォンの他の例を示す斜視図である。
【図８】本発明の一実施形態に係る血圧測定装置の変形例を示すブロック図である。
【図９Ａ】脈波伝播時間の一次関数式から算出した最高血圧とカフを用いる血圧計で計測した最高血圧との相関を示す図である。
【図９Ｂ】脈波伝播時間の一次関数式から算出した最低血圧とカフを用いる血圧計で計測した最低血圧との相関を示す図である。
【発明を実施するための形態】
【００１６】
本発明の実施形態について図面を参照して以下に説明する。
【００１７】
＜＜血圧測定装置の構成＞＞
図１は、本発明の一実施形態に係る血圧測定装置を示すブロック図である。図１に示す血圧測定装置は、ＬＥＤ等の発光部１と、フォトトランジスタ等の受光部２と、Ａ／Ｄ変換部３と、加速度脈波信号算出部４と、特徴量第１抽出部５と、特徴量第２抽出部６と、不揮発性メモリ７と、血圧演算部８とを備えている。発光部１、受光部２、及びＡ／Ｄ変換部３が本発明の「脈波検出部」の一例に該当し、加速度脈波信号算出部４、特徴量第１抽出部５、及び特徴量第２抽出部６が本発明の「特徴量抽出部」の一例に該当し、不揮発性メモリ７が本発明の「特徴量保持部」の一例に該当し、血圧演算部８が本発明の「血圧演算部」の一例に該当し、本発明の「脈波信号処理装置」は、「特徴量抽出部」、「特徴量保持部」及び「血圧演算部」によって構成される。なお、Ａ／Ｄ変換部３を本発明の「脈波検出部」の一例の構成要素から除いて、本発明の「脈波信号処理装置」の一例がＡ／Ｄ変換部３を備えるようにしてもよい。
【００１８】
発光部１及び受光部２はフォトリフレクタを構成している。当該フォトリフレクタの発光面及び受光面の上に人の指１００が置かれると、発光部１から人の指１００内に照射された光の反射光が受光部２によって受光され、光電変換される。血液中のヘモグロビンが光を吸収するため、上記反射光の光量は、血管の膨張収縮の度合いに応じて変動する。したがって、心拍に同期したアナログ電気信号である脈波信号が受光部２から出力される。
【００１９】
なお、発光部１及び受光部２の配置を図１の場合とは逆にしてもよい。すなわち、発光部１を指先側に配置し、受光部２の指の根元側に配置してもよい。また、図１の場合とは異なり、指差し方向に垂直な方向に発光部１及び受光部２が並ぶようにしてもよい。この場合、発光部１を指の爪側に配置し、受光部２の指の腹側に配置してもよく、これとは逆に、発光部１を指の腹側に配置し、受光部２の指の爪側に配置してもよい。本実施形態では反射光を利用して脈波信号を生成する光電脈波検出部を用いたが、透過光を利用して脈波信号を生成する光電脈波検出部を用いても構わない。また、光電脈波検出部以外の脈波検出部を用いることも可能である。また、脈波の測定箇所は、指に限らず、生体の血流循環により生じる脈波を皮膚を通して検出できる箇所であればよく、例えば耳たぶ等であってもよい。
【００２０】
Ａ／Ｄ変換部３は、受光部２から送られてくるアナログ脈波信号をデジタル脈波信号に変換してから、加速度脈波信号算出部４及び特徴量第１抽出部５に出力する。加速度脈波信号算出部４は、脈波信号を２回微分した加速度脈波信号を生成して特徴量第２抽出部６に出力する。
【００２１】
特徴量第１抽出部５は、脈波信号に関する特徴量を抽出して不揮発性メモリ７及び血圧演算部８に出力する。特徴量第２抽出部６は、加速度脈波信号に関する特徴量を抽出して不揮発性メモリ７及び血圧演算部８に出力する。
【００２２】
不揮発性メモリ７は、特徴量第１抽出部５から送られてくる脈波信号に関する特徴量と、特徴量第２抽出部６から送られてくる加速度脈波信号に関する特徴量を記憶する。後述する標準血圧は生体毎に異なるので、図１に示す血圧測定装置を複数の生体の血圧測定が可能な仕様にする場合には、不揮発性メモリ７が生体毎に特徴量を記憶するようにし、血圧測定対象の生体に応じた特徴量が不揮発性メモリ７から読み出されるようにすればよい。
【００２３】
不揮発性メモリ７は、血圧測定装置内に固定的に設けられている記憶媒体であってもよく、血圧測定装置本体に対して着脱可能な記憶媒体であってもよい。また、不揮発性メモリ７は血圧測定装置本体から物理的に離れた場所に設けられ、血圧測定装置本体側と不揮発性メモリ７側との双方にデータ通信部を設けてデータ通信部を介して不揮発性メモリ７への特徴量の書き込み及び不揮発性メモリ７からの特徴量の読み出しを行うようにしてもよい。
【００２４】
血圧演算部８は、過去に不揮発性メモリ７に記憶された特徴量を不揮発性メモリ７から読み出し、過去に不揮発性メモリ７に記憶された脈波信号に関する特徴量及び過去に不揮発性メモリ７に記憶された加速度脈波信号に関する特徴量から算出される値と、特徴量第１抽出部５によって抽出された現在の脈波信号に関する特徴量及び特徴量第２抽出部６によって抽出された現在の加速度脈波信号に関する特徴量から算出される値とから血圧を演算する。
【００２５】
血圧演算部８が不揮発性メモリ７から読み出す過去に不揮発性メモリ７に記憶された特徴量は、過去に不揮発性メモリ７に記憶された特徴量の全てであっても一部であっても良い。一部を読み出す場合、例えば、所定の期間（例えば、直近の１ヶ月間、前年同月同日から直近の１ヶ月間）の特徴量を読み出しても良く、また、例えば、同一種類の他の特徴量と比較して著しく値が異なる特徴量を除いて読み出しても良い。
【００２６】
＜＜血圧演算部が実施する演算の内容＞＞
ここで、脈波信号、脈波信号を１回微分した速度脈波信号、及び脈波信号を２回微分した加速度脈波信号の各波形例を図２に示す。なお、図２中の太実線は脈波信号の波形を示しており、図２中の細実線は速度脈波信号の波形を示しており、図２中の一点鎖線は加速度脈波信号の波形を示している。また、図２の左側の縦軸は脈波信号の値（電圧値）を示しており、図２の右側の縦軸は速度脈波信号、加速度脈波信号の相対値を示している。
【００２７】
以下、図２を参照にして、血圧演算部８が実施する演算の内容をより詳細に説明する。
【００２８】
＜反射波伝搬時間の算出＞
血圧演算部８は、血流の反射波が生体の２点間を伝搬するのに要する時間（反射波伝搬時間）を算出する。
【００２９】
例えば、脈波信号の１周期の最小点ｅでの時間ｔ＿ｅの直後の加速度脈波信号の極大点ｄでの時間ｔ＿ｄから脈波信号の１周期の最小点ｅでの時間ｔ＿ｅを引くことで、血圧演算部８が反射波伝搬時間ｔを算出するとよい。すなわち、下記（１）式により反射波伝搬時間ｔを算出するとよい。
ｔ＝ｔ＿ｄ−ｔ＿ｅ    ・・・（１）
【００３０】
また、例えば、脈波信号の１周期の最小点ｅでの時間ｔ＿ｅの直後の加速度脈波信号の極小点ｃでの時間ｔ＿ｃから、脈波信号の１周期の最小点ｅでの時間ｔ＿ｅの直前の加速度脈波信号の極小点ａでの時間ｔ＿ａを引くことで、血圧演算部８が反射波伝搬時間ｔを算出してもよい。すなわち、下記（２）式により反射波伝搬時間ｔを算出してもよい。
ｔ＝ｔ＿ｃ−ｔ＿ａ    ・・・（２）
【００３１】
また、例えば、脈波信号の１周期の最小点ｅでの時間ｔ＿ｅの直後の加速度脈波信号の極大点ｄでの時間ｔ＿ｄから、脈波信号の１周期の最小点ｅでの時間ｔ＿ｅの直前の加速度脈波信号の極大点ｂでの時間ｔ＿ｂを引くことで、血圧演算部８が反射波伝搬時間を算出してもよい。すなわち、下記（３）式により反射波伝搬時間ｔを算出してもよい。
ｔ＝ｔ＿ｄ−ｔ＿ｂ    ・・・（３）
【００３２】
＜脈波周期の算出＞
血圧演算部８は、脈波周期を算出する。
【００３３】
例えば、脈波信号の１周期の最大点ｆでの時間ｔ＿ｆを２回連続で取得し、後の時間ｔ＿ｆから前の時間ｔ＿ｆを引くことで、血圧演算部８が脈波周期Ｔを算出するとよい。
【００３４】
また、例えば、脈波信号の１周期の最小点ｅでの時間ｔ＿ｅを２回連続で取得し、後の時間ｔ＿ｅから前の時間ｔ＿ｅを引くことで、血圧演算部８が脈波周期Ｔを算出してもよい。同様に、加速度脈波信号の１周期の最大点を２回連続で取得して脈波周期Ｔを算出してもよく、加速度脈波信号の１周期の最小点を２回連続で取得して脈波周期Ｔを算出してもよい。同様に、速度脈波信号の１周期の最大点を２回連続で取得して脈波周期Ｔを算出することもでき、速度脈波信号の１周期の最小点を２回連続で取得して脈波周期Ｔを算出することもできる。このように速度脈波信号に関する特徴量を用いる場合、脈波信号処理装置が脈波信号から速度脈波信号を算出する速度脈波信号算出部を備えるようにすればよい。
【００３５】
＜現在の血圧の算出＞
血圧演算部８は、血圧の測定対象である生体毎に個体差がある生体の標準血圧と、現在の血圧が標準血圧からどれだけ高いか或いは低いかを示す血圧変動とを合わせて、現在の血圧を算出する。以下に示すように、標準血圧は、過去に不揮発性メモリ７に記憶された特徴量から算出することができ、血圧変動は、過去に不揮発性メモリ７に記憶された特徴量と現在の特徴量とから算出することができる。
【００３６】
＜標準血圧の算出＞
血圧演算部８は、血圧の測定対象である生体毎に個体差がある生体の標準血圧を算出する。例えば、加速度脈波信号の極値を用いて、標準血圧を算出するようにしてもよい。標準血圧を算出するために用いる加速度脈波信号の極値としては、例えば、血管の硬さを表すと考えられている加速度脈波信号の３つの極値、すなわち脈波信号の１周期の最小点ｅでの時間ｔ＿ｅの直前の加速度脈波信号の極小点ａでの極小値ｙ＿ａ、脈波信号の１周期の最小点ｅでの時間ｔ＿ｅの直前の加速度脈波信号の極大点ｂでの極大値ｙ＿ｂ、及び脈波信号の１周期の最小点ｅでの時間ｔ＿ｅの直後の加速度脈波信号の極大点ｄでの極大値ｙ＿ｄを用いるようにしてもよい。さらに、血圧演算部８は、反射波伝搬時間ｔも用いて、標準血圧を算出するようにしてもよい。
【００３７】
例えば、極大値ｙ＿ｄから極大値ｙ＿ｂを引いて得られる第１の値を極小値ｙ＿ａで除して第２の値を得て、反射波伝搬時間ｔに係数αを掛けて得られる第３の値を第２の値から引いて指標Ｗを得た場合、すなわち、下記（４）式により指標Ｗを算出した場合、図３Ａ及び図３Ｂから分かるように、指標Ｗの平均値（例えば１日１回特定の人の脈波を測定して算出した指標Ｗの１ヶ月間の平均値）とカフを用いる血圧計で計測した血圧の平均値（例えば１日１回特定の人の血圧を測定した結果の１ヶ月間の平均値）とに高い相関がある。なお、図３Ａ及び図３Ｂは、反射波伝搬時間ｔの単位をｍｓｅｃとし、係数αを０．００１とした時の値を示している。
【数１】
2014008200.tif 000003
【００３８】
図３Ａは、指標Ｗの平均値ＷAVEとカフを用いる血圧計で計測した最高血圧の平均値ＳＢＰAVEとの相関を示す図であり、図中の一次関数は下記（５）式で表される。
ＳＢＰAVE＝２１０ＷAVE＋２０  ・・・（５）
【００３９】
図３Ｂは、指標Ｗの平均値ＷAVEとカフを用いる血圧計で計測した最低血圧の平均値ＤＢＰAVEとの相関を示す図であり、図中の一次関数は下記（６）式で表される。
ＤＢＰAVE＝１７０ＷAVE  ・・・（６）
【００４０】
したがって、血圧演算部８は、指標Ｗの平均値ＷAVEを用いて標準血圧を算出するようにしてもよい。
【００４１】
＜血圧変動の算出＞
血圧演算部８は、現在の血圧が標準血圧からどれだけ高いか或いは低いかを示す血液変動を算出する。例えば、過去の反射波伝搬時間の平均値及び現在の反射波伝搬時間を用いて血液変動を算出するようにしてもよい。さらに、過去の脈波周期の平均値及び現在の脈波周期の平均値も用いて血液変動を算出するようにしてもよい。
【００４２】
例えば、現在の反射波伝搬時間と現在の脈波周期に係数βを掛けた値とを加算した現在の加算値から、反射波伝搬時間と脈波周期に係数βを掛けた値とを加算した加算値の過去の平均値を引いて得られる値に係数γを掛けて得られる値を血圧変動としてもよい。この血圧変動の算出例と上記（５）式及び上記（６）式と組み合わせて血圧演算部８が血圧を演算する場合、現在の最高血圧ＳＢＰCURは下記（７）式で表され、現在の最低血圧ＤＢＰCURは下記（８）式で表される。
ＳＢＰCUR＝２１０ＷAVE＋２０＋γ×（現在のＶ−Ｖの過去の平均値）・・・（７）
ただし、Ｖ＝ｔ＋βＴ
ＤＢＰCUR＝１７０ＷAVE＋γ×（現在のＶ−Ｖの過去の平均値）・・・（８）
ただし、Ｖ＝ｔ＋βＴ
【００４３】
図４Ａは、血圧演算部８が上記（７）式に基づいて演算した現在の最高血圧ＳＢＰCURとカフを用いる血圧計で計測した最高血圧との相関を示す図である。図４Ｂは、血圧演算部８が上記（８）式に基づいて演算した現在の最低血圧ＤＢＰCURとカフを用いる血圧計で計測した最低血圧との相関を示す図である。
【００４４】
脈波伝播時間の一次関数式から算出した血圧すなわち特許文献１の血圧の算出式から算出されていた血圧とカフを用いる血圧計で計測した血圧との相関が低い（図９Ａ及び図９Ｂ参照）のに対して、血圧演算部８が上記（７）式及び上記（８）式に基づいて演算した現在の血圧（ＳＢＰCUR、ＤＢＰCUR）とカフを用いる血圧計で計測した血圧との相関は高い（図４Ａ及び図４Ｂ参照）。なお、図４Ａ及び図４Ｂは、反射波伝搬時間ｔと脈波周期Ｔの単位をｍｓｅｃとし、係数βを０．０１とし、係数γを−０．２とした時の値である。
【００４５】
このように本発明は、過去に保持された特徴量から算出される値と、現在の特徴量から算出される値とから血圧を演算するので、現在の特徴量のみから血圧を演算する従来の技術に比べて、血圧を精度良く演算することができる。
【００４６】
＜＜血圧測定装置が実施する血圧測定処理方法＞＞
次に、図１に示す血圧測定装置が実施する血圧測定処理方法の一例について図５のフローチャートを参照して説明する。例えば、図１に示す血圧測定装置に設けられている操作部（不図示）が、血圧測定の開始を指示するユーザ操作を受け付けると、図５のフローチャートの動作が開始される。なお、図５のステップＳ２０〜Ｓ６０が脈波信号処理方法の一例に該当する。
【００４７】
ステップＳ１０では、発光部１、受光部２、及びＡ／Ｄ変換部３がデジタル信号である脈波信号を生成する。ステップＳ１０によって生成される脈波信号は、数秒から数十秒程度の期間の信号であることが好ましい。
【００４８】
ステップＳ１０に続くステップＳ２０では、加速度脈波信号算出部４が、脈波信号を２回微分処理して加速度脈波信号を算出する。
【００４９】
ステップＳ２０に続くステップＳ３０では、特徴量第１抽出部５が、脈波信号に関する特徴量（例えば、図２中に示されている時間ｔ＿ｅや２回連続の時間ｔ＿ｆなど）を抽出し、特徴量第２抽出部６が、加速度脈波信号に関する特徴量（例えば、図２中に示されている時間ｔ＿ｄや極大値ｙ＿ｄ、極大値ｙ＿ｂ、極小値ｙ＿ａなど）を抽出する。ステップＳ１０によって生成された脈波信号の期間が長く、同一種類の特徴量が複数個抽出できる場合は、任意の１個のみを抽出してもよいが、複数個抽出したのちに平均を求め、その平均を抽出した特徴量としてもよい。複数個抽出したのちに平均を求める場合、時間に関しては周期毎に所定の基準時点（時間ｔ＿ｅ）を設け、その所定の基準時点からの相対時間を用いるとよい。
【００５０】
ステップＳ３０に続くステップＳ４０では、血圧演算部８が、過去（現在の血圧測定開始よりも前）に不揮発性メモリ７に記憶された特徴量を読み出す。
【００５１】
ステップＳ４０に続くステップＳ５０では、特徴量第１抽出部５及び特徴量第２抽出部６によって抽出された現在の特徴量が不揮発性メモリ７に書き込まれる。なお、ステップＳ４０とステップＳ５０との実行順序は逆であっても構わない。また、ステップＳ４０とステップＳ５０とが並列して実行されても構わない。
【００５２】
ステップＳ５０に続くステップＳ６０では、血圧演算部８が、不揮発性メモリ７から読み出した特徴量から算出される値（例えば、上述した指標Ｗの平均値ＷAVEやＶの過去の平均値など）と、現在の特徴量から算出される値（例えば、上述した現在のＶ）とから血圧を演算する。そして、演算された血圧の値は、例えば図１に示す血圧測定装置に設けられている表示部（不図示）に表示させるとよい。
【００５３】
＜＜脈波信号処理プログラム＞＞
図１に示す血圧測定装置は、専用装置によって実現してもよく、脈波検出部を備えるコンピュータ或いは脈波検出部から脈波信号を受け取ることができるコンピュータと、本発明の一実施形態に係る脈波信号処理プログラムとによって実現してもよい。
【００５４】
本発明の一実施形態に係る脈波信号処理プログラムは、電気通信回線を介して、或いは、コンピュータ読み取り可能な記憶媒体（例えばＣＤ−ＲＯＭ）を介して、脈波検出部を備えるコンピュータ或いは脈波検出部から脈波信号を受け取ることができるコンピュータに供給されてインストールされる。そして、インストールされた本発明の一実施形態に係る脈波信号処理プログラムが実行されると、当該コンピュータが図１に示す血圧測定装置の加速度脈波信号算出部４、特徴量第１抽出部５、特徴量第２抽出部６、不揮発性メモリ７、及び血圧演算部８として機能する。
【００５５】
脈波検出部を備えるコンピュータの例としては、図６に示すようにフォトリフレクタを備え、フォトリフレクタの発光面９及び受光面１０が筐体表面に露出しているスマートフォン１１を挙げることができる。また、脈波検出部から脈波信号を受け取ることができるコンピュータの例としては、例えば図７に示すようなＩ／Ｏインターフェース１２を備え、図１に示す血圧測定装置の発光部１、受光部２、及びＡ／Ｄ変換部３に相当する脈波検出ユニット１３のコネクタ部１４から出力される脈波信号をＩ／Ｏインターフェース１２によって受け取ることができるスマートフォン１５を挙げることができる。
【００５６】
当然の事ながら、スマートフォン以外の携帯型コンピュータやデスクトップパーソナルコンピュータなどの据置型コンピュータも、脈波検出部を備えるコンピュータ或いは脈波検出部から脈波信号を受け取ることができるコンピュータとなり得る。
【００５７】
＜＜その他＞＞
以上、本発明に係る実施形態について説明したが、本発明の範囲はこれに限定されるものではなく、発明の主旨を逸脱しない範囲で種々の変更を加えて実行することができる。例えば、加速度脈波信号算出部４をアナログ回路とし、図８に示すように、加速度脈波信号算出部４が受光部２から送られてくるアナログ脈波信号を用いてアナログ加速度脈波信号を生成し、Ａ／Ｄ変換部１６が加速度脈波信号算出部４から送られてくるアナログ加速度脈波信号をデジタル加速度脈波信号に変換してから特徴量第２抽出部６に出力するようにしてもよい。なお、Ａ／Ｄ変換部３とＡ／Ｄ変換部１６は同一のハードウェア構成にすることができるので、Ａ／Ｄ変換部３とＡ／Ｄ変換部１６を共用化してもよい。
【符号の説明】
【００５８】
１  発光部
２  受光部
３、１６  Ａ／Ｄ変換部
４  特徴量第１抽出部
５  加速度脈波信号算出部
６  特徴量第２抽出部
７  不揮発性メモリ
８  血圧演算部
９  フォトリフレクタの発光面
１０  フォトリフレクタの受光面
１１、１５  スマートフォン
１２  Ｉ／Ｏインターフェース
１３  脈波検出ユニット
１４  脈波検出ユニットのコネクタ部
１００  人の指
【図１】
2014008200.tif 000004
【図２】
2014008200.tif 000005
【図３Ａ】
2014008200.tif 000006
【図３Ｂ】
2014008200.tif 000007
【図４Ａ】
2014008200.tif 000008
【図４Ｂ】
2014008200.tif 000009
【図５】
2014008200.tif 000010
【図６】
2014008200.tif 000011
【図７】
2014008200.tif 000012
【図８】
2014008200.tif 000013
【図９Ａ】
2014008200.tif 000014
【図９Ｂ】
2014008200.tif 000015