(57)【要約】  心拍数検出装置は、生体の脈波に応じた脈波信号を取得する脈波信号取得部と、脈波信号に基づき生体の心拍数を導出する演算部と、を備える。演算部は、前回導出した心拍数を内包し且つ所定の大きさを持つ数値範囲内で今回の心拍数を導出する際、体動信号取得部で取得される体動信号に応じて前記数値範囲を可変設定する。或いは、演算部は、脈波信号の内、体動ノイズが重畳する周波数帯域の信号成分を減衰させることで補正脈波信号を生成し、補正脈波信号に基づいて心拍数を導出する。或いは、演算部は、脈波信号の周波数解析により得られる脈波パワースペクトルについて、各周波数に対するパワースペクトル密度で重み付けされた周波数平均値を算出することにより、心拍数を導出する。
2015129557.tif 000002
【特許請求の範囲】
【請求項１】
生体の脈波に応じた脈波信号を取得する脈波信号取得部と、
前記生体の動きに応じた体動信号を取得する体動信号取得部と、
順次取得される前記脈波信号及び前記体動信号に基づき前記生体の心拍数を順次導出する演算部と、を備え、
前記演算部は、前回導出した心拍数を内包し且つ所定の大きさを持つ数値範囲内で今回の心拍数を導出する際、前記体動信号に応じて前記数値範囲を可変設定する、
ことを特徴とする心拍数検出装置。
【請求項２】
前記演算部は、前記体動信号に基づき、前記生体の動き状態を、前記動きの無い状態を含む第１動き状態又は前記第１動き状態よりも前記動きの大きい第２動き状態に分類判定し、
今回の心拍数を導出する際において前記動き状態が前記第２動き状態であると判定されたとき、前記数値範囲の中心を、前回導出した心拍数よりも大きくする、
ことを特徴とする請求項１に記載の心拍数検出装置。
【請求項３】
今回の心拍数を導出する際において前記動き状態が前記第１動き状態であると判定されたとき、前記演算部は、前回導出した心拍数を前記数値範囲の中心に設定する、
ことを特徴とする請求項２に記載の心拍数検出装置。
【請求項４】
前記動き状態が前記第２動き状態から前記第１動き状態に切り替わったとき、前記演算部は、その切り替わり後の所定の期間において、前記数値範囲の中心を、前回導出した心拍数よりも小さくする、
ことを特徴とする請求項２または請求項３に記載の心拍数検出装置。
【請求項５】
前記動き状態が前記第２動き状態から前記第１動き状態に切り替わったとき、前記演算部は、その切り替わり後の所定の期間において、前記数値範囲に依存せずに前記心拍数を導出する、
ことを特徴とする請求項２または請求項３に記載の心拍数検出装置。
【請求項６】
前記演算部は、前記体動信号のパワースペクトルにおける複数のパワースペクトル密度の総和に基づき、前記生体の動き状態に対する前記分類判定を行う、
ことを特徴とする請求項２〜請求項５のいずれか一項に記載の心拍数検出装置。
【請求項７】
前記演算部は、前記総和が所定値以上であるとき、前記動き状態が前記第２動き状態に属すると判定する一方、前記総和が前記所定値未満であるとき、前記動き状態が前記第１動き状態に属すると判定する、
ことを特徴とする請求項６に記載の心拍数検出装置。
【請求項８】
前記演算部は、今回の心拍数を導出するための脈波信号に基づくパワースペクトルにおいて、前記数値範囲に対応する帯域内でパワースペクトル密度が最大となる周波数を同定することにより、今回の心拍数を導出する、
ことを特徴とする請求項１〜請求項７のいずれか一項に記載の心拍数検出装置。
【請求項９】
前記演算部は、今回の心拍数を導出する際、
前記体動信号のパワースペクトルにおいて、パワースペクトル密度が所定の閾値以上となっている帯域を含む減衰対象帯域を設定し、
前記脈波信号取得部にて取得された脈波信号のパワースペクトルにおいて、前記減衰対象帯域内のパワースペクトル密度を減衰させることで補正パワースペクトルを導出し、
前記補正パワースペクトルにおいて、前記数値範囲に対応する帯域内でパワースペクトル密度が最大となる周波数を同定することにより、今回の心拍数を導出する、
ことを特徴とする請求項８に記載の心拍数検出装置。
【請求項１０】
前記演算部は、時間領域上の信号として取得された前記脈波信号及び前記体動信号の夫々を周波数領域上の信号に変換する時間／周波数変換処理を介して前記心拍数を導出し、前記時間／周波数変換処理における周波数分解能を、前記脈波信号及び前記体動信号間で同じにする、
ことを特徴とする請求項１〜請求項９のいずれか一項に記載の心拍数検出装置。
【請求項１１】
前記演算部は、前回導出した心拍数を内包し且つ所定の大きさを持つ数値範囲内で今回の心拍数を導出する際、前記体動信号のベクトル値と各軸における振幅変化の総和を各々算出し、その算出結果に応じて前記数値範囲を可変的に設定することを特徴とする請求項１に記載の心拍数検出装置。
【請求項１２】
前記演算部は、前記数値範囲を設定する際、前回導出した心拍数と所定の閾値との比較結果も併せて参照することを特徴とする請求項１１に記載の心拍数検出装置。
【請求項１３】
前記演算部は、前記脈波信号の周波数解析により得られる脈波パワースペクトルについて、前記数値範囲に対応する帯域内で、各周波数に対するパワースペクトル密度で重み付けされた周波数平均値を算出することにより、前記心拍数を導出することを特徴とする請求項１１または請求項１２に記載の心拍数検出装置。
【請求項１４】
前記演算部は、前記脈波信号の周波数解析により得られる脈波パワースペクトルについて、前記数値範囲に対応する帯域内で、パワースペクトル密度が最大となる周波数を同定することにより、前記心拍数を導出することを特徴とする請求項１１または請求項１２に記載の心拍数検出装置。
【請求項１５】
前記演算部は、前記心拍数を導出するに際して、同一条件で前記脈波信号と前記体動信号の周波数解析を行うことにより、脈波パワースペクトルと体動パワースペクトルを各々算出することを特徴とする請求項１１〜請求項１４のいずれか一項に記載の心拍数検出装置。
【請求項１６】
前記演算部は、前記体動パワースペクトルのピーク位置を解析し、その解析結果に応じて体動ノイズ成分を減衰するように前記脈波パワースペクトルを補正することを特徴とする請求項１５に記載の心拍数検出装置。
【請求項１７】
前記演算部は、複数検出されたピーク位置のうちパワースペクトル密度の上位順に所定数のピーク位置だけを抽出して前記脈波パワースペクトルの補正に反映させることを特徴とする請求項１６に記載の心拍数検出装置。
【請求項１８】
前記演算部は、前記ピーク位置を中心とした前後の周波数帯域を体動ノイズ周波数帯域として設定し、前記体動ノイズ周波数帯域内の脈波パワースペクトルに対して所定の補正係数を乗算することを特徴とする請求項１６または請求項１７に記載の心拍数検出装置。
【請求項１９】
前記演算部は、前記補正係数を周波数に応じて可変的に設定することを特徴とする請求項１８に記載の心拍数検出装置。
【請求項２０】
生体の脈波に応じた脈波信号を取得する脈波信号取得部と、
前記脈波信号に基づき前記生体の心拍数を導出する演算部と、を備え、
前記演算部は、前記脈波信号の内、体動ノイズが重畳する周波数帯域の信号成分を減衰させることで補正脈波信号を生成し、前記補正脈波信号に基づいて前記心拍数を導出することを特徴とする心拍数検出装置。
【請求項２１】
前記生体の動きに応じた体動信号を取得する体動信号取得部をさらに備え、
前記演算部は、
前記体動信号のパワースペクトルにおいて、パワースペクトル密度が所定の閾値以上となっている帯域を含む減衰対象帯域を設定し、
前記脈波信号取得部にて取得された脈波信号のパワースペクトルにおいて、前記減衰対象帯域内のパワースペクトル密度を減衰させることで補正パワースペクトルを導出し、
前記補正パワースペクトルに基づき、前記心拍数を導出する、
ことを特徴とする請求項２０に記載の心拍数検出装置。
【請求項２２】
前記脈波信号取得部にて順次取得される脈波信号を原脈波信号とし、
前記演算部は、前記原脈波信号の内、所定周波数未満の信号成分を減衰させることで前記補正脈波信号を生成し、前記補正脈波信号に基づいて前記心拍数を導出する、
ことを特徴とする請求項２０に記載の心拍数検出装置。
【請求項２３】
前記演算部は、前記脈波信号取得部にて順次取得される前記原脈波信号の信号値を用いて前記心拍数を順次導出し、今回の心拍数を、前回導出した心拍数を内包し且つ所定の大きさを持った数値範囲内で導出する、
ことを特徴とする請求項２２に記載の心拍数検出装置。
【請求項２４】
前記演算部は、今回の心拍数を導出するための原脈波信号のパワースペクトルにおいて、前記所定周波数未満のパワースペクトル密度を減衰させることで前記補正脈波信号のパワースペクトルを生成し、前記補正脈波信号のパワースペクトルにおいて、前記数値範囲に対応する帯域内でパワースペクトル密度が最大となる周波数を同定することにより、今回の心拍数を導出する、
ことを特徴とする請求項２３に記載の心拍数検出装置。
【請求項２５】
前記生体の動きに応じた体動信号を取得する体動信号取得部をさらに備え、
前記演算部は、前記心拍数を導出するに際して、同一条件で前記脈波信号と前記体動信号の周波数解析を行うことにより、脈波パワースペクトルと体動パワースペクトルを各々算出した上で、前記体動パワースペクトルのピーク位置を解析し、その解析結果に応じて体動ノイズ成分を減衰するように前記脈波パワースペクトルを補正することを特徴とする請求項２０に記載の心拍数検出装置。
【請求項２６】
生体の脈波に応じた脈波信号を取得する脈波信号取得部と、
前記脈波信号に基づき前記生体の心拍数を導出する演算部と、を備え、
前記演算部は、前記脈波信号の周波数解析により得られる脈波パワースペクトルについて、各周波数に対するパワースペクトル密度で重み付けされた周波数平均値を算出することにより、前記心拍数を導出することを特徴とする心拍数検出装置。
【請求項２７】
前記脈波信号取得部は、前記生体に光を照射する発光部及び前記発光部からの光に基づき前記生体を透過又は反射した光を受光する受光部を有し、前記受光部の受光結果に基づき前記脈波信号を取得する
ことを特徴とする請求項１〜請求項２６のいずれか一項に記載の心拍数検出装置。
【請求項２８】
前記受光部は、前記生体を透過又は反射した光を各々に受光する複数の受光素子を有し、
前記発光部は、前記複数の受光素子の間に配置される
ことを特徴とする請求項２７に記載の心拍数検出装置。
【請求項２９】
前記心拍数の導出対象となる解析周波数帯域は、０．６Ｈｚ〜４．０Ｈｚに限定されていることを特徴とする請求項１〜請求項２８のいずれか一項に記載の心拍数検出装置。
【請求項３０】
前記脈波信号取得部は、前記生体の手首に装着されるものであり、前記体動信号によって示される前記生体の動きは、前記生体の前腕の長手方向の動きを含むことを特徴とする請求項１〜請求項２９のいずれか一項に記載の心拍数検出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、心拍数検出装置に関する。
【背景技術】
【０００２】
生体に光を照射して生体を透過した光の受光信号に基づき、生体の脈波信号を取得する脈波センサが知られており（例えば特許文献１参照）、脈波信号から生体の心拍数（ハートレート値）等を求めることができる。【０００３】
脈波信号から生体の心拍数を算出する方法としてピーク解析方法及びＦＦＴ解析方法が存在する。ピーク解析方法では、脈波信号の時系列波形において、ピーク（極大値）の位置を特定し、隣接するピークの時間間隔から心拍数を算出する。ＦＦＴ解析方法では、ＦＦＴ（高速フーリエ変換）によって脈波信号の時系列データから脈波信号のパワースペクトルを算出し、パワースペクトル密度が最大となる周波数を同定することにより心拍数を算出する。【０００４】
また、心拍数を順次算出するシステムにおいて、今回の心拍数を算出する際、前回算出した心拍数を中心とする数値範囲を設定し、当該数値範囲内に限定して今回の心拍数を同定する帯域制限技術も提案されている（非特許文献１参照）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平０５-１６１６１５号公報
【非特許文献】
【０００６】
【非特許文献１】巻口、他３名，「スマートフォンの音声入出力端子をインターフェースとする脈波測定装置の実装と評価」，情報処理学会インタラクション２０１２，情報処理学会，平成２４年３月，ｐ．５９３−５９８
【発明の概要】
【発明が解決しようとする課題】
【０００７】
生体の心拍数の短時間における変化量は限られており、今回の検出心拍数が前回の検出心拍数と大きくかけ離れているとは考えにくい。上述の帯域制限技術は、この知見を利用したものであり、当該技術により体動等によるノイズの影響を低減することができる（心拍数検出値がノイズに基づく異常値になることが抑制される）。但し、帯域制限の条件が不適切であると正確な心拍数検出が難しくなるため、どのような条件で帯域制限をかけるかが重要となる。心拍数検出技術において、正確な心拍数検出が有益であることは言うまでもない。【０００８】
正確な心拍数検出の妨げとなるノイズの要因の一つは、生体の動き（すなわち体動）である。体動等によるノイズを低減することができれば、正確な心拍数検出が可能となる。【０００９】
そこで、本発明は、正確な心拍数検出を可能とする心拍数検出装置を提供することを目的とする。
【課題を解決するための手段】
【００１０】
本明細書中に開示されている心拍数検出装置は、生体の脈波に応じた脈波信号を取得する脈波信号取得部と、前記生体の動きに応じた体動信号を取得する体動信号取得部と、順次取得される前記脈波信号及び前記体動信号に基づき前記生体の心拍数を順次導出する演算部と、を備え、前記演算部は、前回導出した心拍数を内包し且つ所定の大きさを持つ数値範囲内で今回の心拍数を導出する際、前記体動信号に応じて前記数値範囲を可変設定する構成とされている。【００１１】
また、本明細書中に開示されている心拍数検出装置は、生体の脈波に応じた脈波信号を取得する脈波信号取得部と、前記脈波信号に基づき前記生体の心拍数を導出する演算部とを備え、前記演算部は、前記脈波信号の内、体動ノイズが重畳する周波数帯域の信号成分を減衰させることで補正脈波信号を生成し、前記補正脈波信号に基づいて前記心拍数を導出する構成とされている。【００１２】
また、本明細書中に開示されている心拍数検出装置は、生体の脈波に応じた脈波信号を取得する脈波信号取得部と、前記脈波信号に基づき前記生体の心拍数を導出する演算部とを備え、前記演算部は、前記脈波信号の周波数解析により得られる脈波パワースペクトルについて、各周波数に対するパワースペクトル密度で重み付けされた周波数平均値を算出することにより、前記心拍数を導出する構成とされている。【００１３】
また、本明細書中に開示されている心拍数検出方法は、生体の脈波に応じた脈波信号を取得する脈波信号取得ステップと、前記生体の動きに応じた体動信号を取得する体動信号取得ステップと、順次取得される前記脈波信号及び前記体動信号に基づき前記生体の心拍数を順次導出する演算ステップと、を備え、前記演算ステップでは、前回導出した心拍数を内包し且つ所定の大きさを持つ数値範囲内で今回の心拍数を導出する際、前記体動信号に応じて前記数値範囲を可変設定する構成とされている。【００１４】
また、本明細書中に開示されている心拍数検出方法は、生体の脈波に応じた脈波信号を取得する脈波信号取得ステップと、前記脈波信号に基づき前記生体の心拍数を導出する演算ステップと、を備え、前記演算ステップでは、前記脈波信号の内、体動ノイズが重畳する周波数帯域の信号成分を減衰させることで補正脈波信号を生成し、前記補正脈波信号に基づいて前記心拍数を導出する構成とされている。【００１５】
また、本明細書中に開示されている心拍数検出方法は、生体の脈波に応じた脈波信号を取得する脈波信号取得ステップと、前記脈波信号に基づき前記生体の心拍数を導出する演算ステップと、を備え、前記演算ステップでは、前記脈波信号の周波数解析により得られる脈波パワースペクトルについて、各周波数に対するパワースペクトル密度で重み付けされた周波数平均値を算出することにより、前記心拍数を導出する構成とされている。【００１６】
なお、本発明に関連する上記以外の特徴、要素、ステップ、利点、及び、特性については、以下に続く詳細な説明やこれに関する添付の図面によって、さらに明らかとなる。
【発明の効果】
【００１７】
本発明によれば、正確な心拍数検出を可能とする心拍数検出装置及び心拍数検出方法を提供することが可能である。
【図面の簡単な説明】
【００１８】
【図１】手首での脈波測定の原理を説明するための模式図である。
【図２】生体内における光の減衰量が時間的に変化する様子を示す波形図である。
【図３】本発明の第１実施形態に係る脈波センサのブロック図である。
【図４】本発明の第１実施形態に係る光センサ部の回路例と、光センサ部、フィルタ部及び制御部の関係とを示す図である。
【図５】本発明の第１実施形態に係る脈波ＦＦＴ解析の解析条件を示す図である。
【図６】本発明の第１実施形態に係る心拍数の導出フローチャートである。
【図７】脈波データが補完される様子を示した図である。
【図８】図６における脈波ＦＦＴ解析（Ａ１）の詳細フローチャートである。
【図９】図６における脈波ＦＦＴ解析（Ａ２）の詳細フローチャートである。
【図１０】本発明の第１実施形態に係る検索対象帯域（帯域制限技術）を示す図である。
【図１１】図６における脈波ＦＦＴ解析（Ａ３）の詳細フローチャートである。
【図１２】本発明の第１実施形態に係り、減衰補正前後の脈波信号のパワースペクトルを示す図である。
【図１３】本発明の第１実施形態の実験による心拍数測定結果を示す図である。
【図１４】本発明の第１実施形態の実験による心拍数測定結果を示す図である。
【図１５】本発明の第２実施形態に係り、生体と、本体ユニット及びベルトと、Ｘ、Ｙ及びＺ軸との関係を示す図である。
【図１６】本発明の第２実施形態に係る制御部の内部ブロック図である。
【図１７】本発明の第２実施形態に係る心拍数の導出フローチャートである。
【図１８】本発明の第２実施形態に係り、加速度信号に基づく解析動作のフローチャートである。
【図１９】図１７における脈波ＦＦＴ解析（Ｂ３）の詳細フローチャートである。
【図２０】本発明の第２実施形態に係るパワースペクトルの補正例を示す図である。
【図２１】本発明の第２実施形態に係る検索対象帯域の設定内容を示す図である。
【図２２】本発明の第２実施形態の変形技術に係り、帯域制限の一時停止方法を説明するための図である。
【図２３】本発明の第２実施形態の実験による心拍数測定結果を示す図である。
【図２４】本発明の第２実施形態の実験による心拍数測定結果を示す図である。
【図２５】本発明の第３実施形態に係る光センサ部の構成を示す図である。
【図２６】本発明の第４実施形態に係る解析動作のフローチャートである。
【図２７】ピーク位置検出に基づく体動ノイズキャンセル動作の一例を示す図である。
【図２８】ベクトル値Ｖと振幅変化総和Ｄの定義式を示す図である。
【図２９】体動モードの一例を示すテーブルである。
【図３０】ＰＳＤ最大周波数に基づくＨＲ導出時の問題点を説明するための図である。
【図３１】周波数平均値Ａａｖｅの算出式を示す図である。
【図３２】本発明の第４実施形態の実験による心拍数測定結果を示す図である。
【発明を実施するための形態】
【００１９】
以下、本発明の実施形態の例を、図面を参照して具体的に説明する。参照される各図において、同一の部分には同一の符号を付し、同一の部分に関する重複する説明を原則として省略する。尚、本明細書では、記述の簡略化上、情報、信号、物理量、状態量又は部材等を参照する記号又は符号を記すことによって、該記号又は符号に対応する情報、信号、物理量、状態量又は部材等の名称を省略又は略記することがある。また、後述の任意のフローチャートにおいて、任意の複数のステップにおける複数の処理は、処理内容に矛盾が生じない範囲で、任意に実行順序を変更できる又は並列に実行できる。
【００２０】
＜＜脈波測定の原理＞＞
図１は、手首での脈波測定の原理を説明するための模式図であり、図２は、生体内における光の減衰量（吸光度）が時間的に変化する様子を示す波形図である。後述の各実施形態では、以下の原理に基づく脈波測定が行われる。
【００２１】
容積脈波法による脈波測定では、例えば、図１に示すように、測定窓に押し当てられた生体の一部（図１では手首）に向けて発光部（ＬＥＤ（Light Emitting Diode）など）から光が照射され、体内を透過して体外に出てくる光の強度が受光部（フォトダイオードやフォトトランジスタなど）で検出される。ここで、図２に示したように、生体組織や静脈血（脱酸素化ヘモグロビンＨｂ）による光の減衰量（吸光度）は一定であるが、動脈血（酸素化ヘモグロビンＨｂＯ2）による光の減衰量（吸光度）は拍動によって時間的に変動する。従って、可視領域から近赤外領域にある「生体の窓」（光が生体を透過しやすい波長領域）を利用し、受光部の受光結果に基づき末梢動脈の吸光度変化を測定することで、非侵襲で容積脈波を測定することができる。
【００２２】
尚、図１では、図示の便宜上、発光部及び受光部を有する脈波センサを手首の背側（外側）に装着した様子が描写されているが、脈波センサの装着位置についてはこれに限定されるものではなく、手首の腹側（内側）であってもよいし、他の部位（指先、指の第３関節、額、眉間、鼻先、頬、眼下、こめかみ、耳たぶ、耳穴など）であってもよい。
【００２３】
また、心臓及び自立神経の支配を受けている脈波は、常に一定の挙動を示すものではなく、被験者の状態によって様々な変化（揺らぎ）を生じるものである。従って、脈波の変化（揺らぎ）を解析することにより、被験者の様々な身体情報を得ることができる。例えば、心拍数からは、被験者の運動能力や緊張度などを知ることができ、心拍変動からは、被験者の疲労度、快眠度、及び、ストレスの大きさなどを知ることができる。また、脈波を時間軸で２回微分することにより得られる加速度脈波からは、被験者の血管年齢や動脈硬化度などを知ることができる。
【００２４】
＜＜第１実施形態＞＞
本発明の第１実施形態を説明する。図３には、第１実施形態に係る脈波センサ１のブロック図が、生体２及び外部機器ＥＥと共に示されている。脈波センサ１は、本体ユニット１０と、本体ユニット１０の両端部に取り付けられて生体２（具体的には手首）に巻き回されるベルト２０と、を備えた腕輪構造（腕時計型構造）を有する。ベルト２０の素材として、皮革、金属、樹脂などを用いることができる。本明細書において、被験者とは生体２を有する人間を指す。
【００２５】
本体ユニット１０には、光センサ部１１、フィルタ部１２、制御部１３、表示部１４、通信部１５、電源部１６及び体動センサ部１７が設けられる。
【００２６】
光センサ部１１は、本体ユニット１０の裏面（即ち、生体２と対向する側の面）に設けられており、発光部１１Ａ及び受光部１１Ｂを有する。光センサ部１１では、ＬＥＤ等の発光部１１Ａから生体２に光を照射し、照射した光の内、生体２内を透過した光の強度を受光部１１Ｂで検出することにより、脈波信号を取得する。尚、ここでは、発光部１１Ａと受光部１１Ｂが生体２に対していずれも同じ側に設けられた反射型構成（図１の実線矢印を参照）が採用されているが、発光部１１Ａと受光部１１Ｂが生体２を挟んで互いに反対側に設けられる透過型構成（図１の破線矢印を参照）が採用されても良い。反射型構成では、厳密に表現すると、発光部１１Ａから生体２に光を照射され、照射された光の内、生体２の一部を透過してから生体２内部で反射した光が受光部１１Ｂにて受光される。一方、透過型構成では、発光部１１Ａから生体２の一面に対して光を照射され、照射された光の内、生体２を貫通するように透過して生体２の反対面（上記一面に対する反対面）から出た光が受光部１１Ｂにて受光される。まとめると、受光部１１Ｂは、発光部１１Ａからの光に基づき生体２を透過又は反射した光を受光するものであって良い。
【００２７】
フィルタ部１２は、光センサ部１１の出力信号にフィルタ処理及び増幅処理を施して制御部１３に伝達する。制御部１３は、脈波センサ１全体の動作を統括的に制御するほか、フィルタ部１２の出力信号に各種の信号処理を施すことにより、脈波に関する種々な情報（以下、脈波情報と呼ぶ）を取得する。脈波情報は、例えば、心拍数、脈波の揺らぎ、心拍変動及び加速度脈波を含む。ＭＰＵ（micro processing unit）等を用いて、制御部１３を形成すると良い。表示部１４は、本体ユニット１０の表面（生体２と対向しない側の面）に設けられた液晶表示パネル等から成り、脈波情報や日付や時間に関する情報等を含む表示情報を表示する。表示部１４は、腕時計の文字盤面に相当すると言える。通信部１５は、脈波センサ１と異なる外部機器ＥＥとの間で任意の通信を行う。通信部１５は、脈波情報を含む脈波センサ１の測定データを、外部機器ＥＥに無線又は有線で送信できる。外部機器ＥＥは、パーソナルコンピュータや携帯電話機などの任意の電子機器であり、ネットワーク網を介して通信部１５と接続されていても良い。電源部１６は、バッテリとＤＣ／ＤＣコンバータを含み、バッテリからの入力電圧を所望の出力電圧に変換して脈波センサ１の各部に供給する。体動センサ部１７は、第１実施形態では利用されないので、後述の他の実施形態で説明する。
【００２８】
図４に、光センサ部１１の回路例と共に、光センサ部１１、フィルタ部１２及び制御部１３の関係を示す。制御部１３は、Ａ／Ｄ変換部３１、演算部３２及びメモリ３３を有する。受光部１１Ｂとしてのフォトトランジスタには、上記透過した光の強度に応じた電流が流れ、当該電流の大きさに応じた電圧信号Ｐａｏがフィルタ部１２に入力される。フィルタ部１２は、電圧信号Ｐａｏにフィルタ処理及び増幅処理を施して得られるアナログ電圧信号ＰａをＡ／Ｄ変換部３１に入力する。Ａ／Ｄ変換部３１（例えば、１０ビット、サンプリングレート：８Ｈｚ）は、アナログ電圧信号Ｐａをデジタル電圧信号Ｐｄに変換して出力する。信号Ｐａｏ、Ｐａ及びＰｄの夫々は、生体２の脈波を示す脈波信号の一種である。演算部３２は、脈波信号Ｐｄに対して様々な演算処理を施す（詳細は後述）。メモリ３３は、演算部３２にて実行される処理を規定するプログラムの保存用のプログラムメモリ、及び、演算部３２にて利用又は算出される各種データを一時記憶するデータメモリを有する。
【００２９】
［脈波ＦＦＴ解析の解析条件］
制御部１３は、心拍数の導出に当たり、脈波信号に対して、ＦＦＴ（高速フーリエ変換）を用いた解析である脈波ＦＦＴ解析を実行する。図５の（ａ）欄〜（ｄ）欄を参照して脈波ＦＦＴ解析の解析条件を示す。脈波ＦＦＴ解析の解析条件には、サンプリング周波数、解析データ数、実行間隔及び周波数分解能が含まれ、更に、解析帯域を定める解析下限周波数ｆLL及び解析上限周波数ｆHHが含まれる（図５の（ａ）欄〜（ｄ）欄参照）。
【００３０】
説明の明確化のため、時刻を表す記号として整数ｔを導入する。任意の整数ｉに関し、“ｔ＝ｉ＋１”である時刻は“ｔ＝ｉ”である時刻よりも単位時間Δｔだけ遅い。Ａ／Ｄ変換部によるアナログ信号のデジタル信号への変換をサンプリングと呼ぶ。Ａ／Ｄ変換部３１のサンプリング周波数は“１／Δｔ”である（図５の（ａ）欄参照）。時刻ｔ＝ｉのサンプリングにて取得される脈波信号Ｐｄのデジタル値を“Ｐｄ［ｉ］”にて表す。即ち、時刻ｔ＝ｉにて、第ｉ番目の脈波信号値である値Ｐｄ［ｉ］が得られるものとする。以下では、Ｐｄ［ｉ］にて値が表される脈波信号を、脈波データと呼ぶことがある。
【００３１】
図５の（ｂ）欄に示す如く、パラメータＰＲAは、１回の脈波ＦＦＴ解析に用いるデータ数（即ち、脈波データの個数）である解析データ数を示す。脈波ＦＦＴ解析では、その時点で得られている最新のＰＲA個の脈波データを解析対象データ列として用いてＦＦＴを行う。従って例えば、時刻ｔ＝ｉでは、Ｐｄ［ｉ−ＰＲA＋１］〜Ｐｄ［ｉ］が解析対象データ列を形成する。
【００３２】
図５の（ｃ）欄に示す如く、脈波ＦＦＴ解析は、ＰＲB個分の脈波データが得られるごとに実行される。従って例えば、Ｐｄ［ｉ−ＰＲA＋１］〜Ｐｄ［ｉ］を解析対象データ列として用いた第ｊ回目の脈波ＦＦＴ解析が実行された後、Ｐｄ［ｉ＋１］〜Ｐｄ［ｉ＋ＰＲB］が取得されると、Ｐｄ［ｉ−ＰＲA＋１＋ＰＲB］〜Ｐｄ［ｉ＋ＰＲB］を解析対象データ列として用いた第（ｊ＋１）回目の脈波ＦＦＴ解析が実行される。ＰＲBは、脈波ＦＦＴ解析の実行間隔を表している。例えば、実行間隔を１秒相当に設定する場合、ＰＲBはサンプリング周波数１／Δｔと一致する。
【００３３】
脈波ＦＦＴ解析により、図５の（ｄ）欄に示すような脈波信号のパワースペクトルが得られる。パワースペクトルは、第１〜第ｍ要素帯域の夫々におけるパワースペクトル密度（以下ＰＳＤ値という）を、情報として含んでいる。ｍは２以上の所定の整数である。任意の整数ｉに関し、第ｉ要素帯域の上限周波数は、第（ｉ＋１）要素帯域の下限周波数と一致している。第１要素帯域の中心周波数及び第ｍ要素帯域の中心周波数は、脈波ＦＦＴ解析における解析下限周波数ｆLL及び解析上限周波数ｆHHである。第１〜第ｍ要素帯域において要素帯域の帯域幅は共通であり、１つの要素帯域の帯域幅が周波数分解能Ｒに相当する。
【００３４】
周波数ｆLL及びｆHHは、心拍数がとりうる値を考慮して設定され、それぞれ、例えば、０．６Ｈｚ（ヘルツ）、４．０Ｈｚである。周波数ｆLLを０．６Ｈｚ以外にすることもできるし、周波数ｆHHを４．０Ｈｚ以外にすることもできる（但し、勿論、ｆLL＜ｆHH）。サンプリング周波数は、解析上限周波数ｆHHの２倍以上の任意の周波数に設定され、例えば、８、１６、３２、６４又は１２８Ｈｚであって良い。解析データ数ＰＲAは、任意の所定値とされ、例えば、１２８、２５６、５１２、１０２８又は２０４８であって良い。実行間隔ＰＲBも任意の所定値とされ、例えば、０．５秒相当、１秒相当、２秒相当とされる。周波数分解能Ｒも任意の所定値とされる。ここでは、周波数分解能Ｒは、０．０３１２５Ｈｚであるとする。０．０３１２５Ｈｚは、心拍数の単位で表すと、１．８７５ｂｐｍ（beat per minute）に相当する。
【００３５】
サンプリング周波数、解析データ数ＰＲA及び周波数分解能の増大によって、心拍数の検出精度（検出分解能を含む）の向上が期待されるが、一方で演算負荷が重くなる。必要な精度と演算負荷のバランスを考慮して解析条件を決定すると良い。実行間隔ＰＲBも同様である。
【００３６】
［心拍数の導出方法］
図６は、制御部１３による、心拍数を導出するためのフローチャートである。まず、ステップＳ１１において、Ａ／Ｄ変換部３１は、所定のサンプリング周波数にて脈波データの取り込み（即ちサンプリング）を開始する。順次得られる脈波データはメモリ３３に保存される。ステップＳ１１に続くステップＳ１２において、制御部１３は、これまでに取り込まれた脈波データの個数ｎを、所定の第１閾値（ＰＲA−Ｑ）、及び、解析データ数と一致する第２閾値ＰＲAと比較する。Ｑは所定の正の整数（例えば４０）である。第１不等式“ＰＲA−Ｑ≦ｎ＜ＰＲA”の成立時にはステップＳ１２からステップＳ１３に進み、第２不等式“ＰＲA≦ｎ”の成立時にはステップＳ１２からステップＳ１９に進む。不等式“ｎ＜ＰＲA−Ｑ”の成立時には、データ数が少なすぎるため、第１不等式が成立するまで脈波データの取り込みを継続しつつステップＳ１２で待機する。
【００３７】
ステップＳ１３の段階で取得されている脈波データは、Ｐｄ［１］〜Ｐｄ［ｎ］であり、それらの総個数は解析データ数ＰＲAより少ない。故に、時刻ｔにて実行されるステップＳ１３において、演算部３２は、時刻ｔ＝ｎ＋１〜時刻ｔ＝ＰＲAにおける脈波データＰｄ［ｎ＋１］〜Ｐｄ［ＰＲA］として、所定の補完用データを補完する。図７に、この補完の様子を示す。補完用データは、脈波信号の中心電位を持つ。例えば、脈波信号が、０Ｖ〜３Ｖの範囲内で１．５Ｖを中心に変動する脈動信号（想定脈波強度：５０ｍＶ〜５００ｍＶ（ｔｙｐ．２００ｍＶ））である場合、補完用データは１．５Ｖに相当する固定値を持つ。尚、Ｐｄ［１］〜Ｐｄ［ｎ］の平均値を補完用データとして用いても良いし、上記中心電位と異なる電位に相当する固定値を補完用データとして用いても構わない。ステップＳ１３により、実測によって得られた脈波データＰｄ［１］〜Ｐｄ［ｎ］と、補完された脈波データＰｄ［ｎ＋１］〜Ｐｄ［ＰＲA］とから成る解析対象データ列が形成される。
【００３８】
ステップＳ１３に続くステップＳ１４において、演算部３２は、解析対象データ列が存在する区間だけ０以外の値を持つ窓関数（例えばハニング関数）を脈波データに対して掛け合わせる窓関数処理を実行する。その後、ステップＳ１５において、演算部３２は、これから行うべき脈波ＦＦＴ解析が、何回目の脈波ＦＦＴ解析であるのかを確認する。そして、これから行うべき脈波ＦＦＴ解析が、１回目の脈波ＦＦＴ解析である場合には、ステップＳ１６に進んで脈波ＦＦＴ解析Ａ１を実行し、２回目以降の脈波ＦＦＴ解析である場合には、ステップＳ１７に進んで脈波ＦＦＴ解析Ａ２を実行する。脈波ＦＦＴ解析Ａ１及びＡ２は、ステップＳ１４の窓関数処理を経た脈波データ（即ち、脈波データ中の解析対象データ列）に対して実行される。
【００３９】
ステップＳ１６の脈波ＦＦＴ解析Ａ１又はステップＳ１７の脈波ＦＦＴ解析Ａ２を終えると、ステップＳ１８へ進む。ステップＳ１８において、Ａ／Ｄ変換部３１により、新たにＰＲB個分の脈波データが取り込まれ、ステップＳ１２に戻る。新たに取り込まれた脈波データはメモリ３３に保存される。
【００４０】
一方、“ＰＲA≦ｎ”の場合、既に取得されている脈波データの総個数は解析データ数ＰＲA以上であるので、上述のような補完は不要である。ステップＳ１９において、演算部３２は、上述の窓関数処理と同様の窓関数処理を実行する。ステップＳ１９における解析対象データ列は、実測によって得られた脈波データＰｄ[ｎ−ＰＲA＋１]〜Ｐｄ[ｎ]から成る（図５の（ｂ）欄参照）。その後、ステップＳ２０において、ステップＳ１９の窓関数処理を経た脈波データ（即ち、脈波データ中の解析対象データ列）に対し脈波ＦＦＴ解析Ａ３を実行する。
【００４１】
ステップＳ２０の脈波ＦＦＴ解析Ａ３を終えると、ステップＳ２１へ進む。ステップＳ２１において、Ａ／Ｄ変換部３１により、新たにＰＲB個分の脈波データが取り込まれた後、ステップＳ１９に戻る。新たに取り込まれた脈波データはメモリ３３に保存される。この際、メモリ３３に保存されている脈波データの内、より古くから保存されているＰＲB個分の脈波データをメモリ３３から削除しても良い。
【００４２】
脈波ＦＦＴ解析の解析条件は、脈波ＦＦＴ解析Ａ１、Ａ２及びＡ３間で共通である。脈波ＦＦＴ解析Ａ１、Ａ２及びＡ３において、周波数ｆLL及びｆHHは、夫々、例えば、０．６３５Ｈｚ及び４．０Ｈｚである。
【００４３】
――脈波ＦＦＴ解析Ａ１（図８）――
図８を参照し、脈波ＦＦＴ解析Ａ１について説明する。脈波ＦＦＴ解析Ａ１は、ステップＳ３１〜Ｓ３３の処理から成る。ステップＳ３１において、演算部３２は、窓関数処理を経た脈波データ（即ち、脈波データ中の解析対象データ列）に対してＦＦＴを適用することで、脈波信号のパワースペクトルを算出する。脈波信号のパワースペクトルは、上述したように、脈波信号の第１〜第ｍ要素帯域の夫々におけるＰＳＤ値を情報として含んでいる。続くステップＳ３２において、演算部３２は、ステップＳ３１で得られた第１〜第ｍ要素帯域におけるｍ個のＰＳＤ値の内、最大のＰＳＤ値に対応する要素帯域の周波数を、ＰＳＤ最大周波数Ｆ［ｃｙ］として検索及び特定する。要素帯域の周波数は、要素帯域に属する周波数（例えば、中心周波数）を指す。
【００４４】
記号ｃｙは、今回実行する脈波ＦＦＴ解析の回数、即ち、今回実行する脈波ＦＦＴ解析が第何回目の脈波ＦＦＴ解析であるのかを示す。回数ｃｙは、脈波ＦＦＴ解析Ａ１、Ａ２及びＡ３の全実行回数である。従って例えば、脈波ＦＦＴ解析Ａ２及びＡ３を未だ１回も行っておらず、１回目の脈波ＦＦＴ解析として脈波ＦＦＴ解析Ａ１を行う時には、ｃｙ＝１である。また例えば、脈波ＦＦＴ解析Ａ１を１回目の脈波ＦＦＴ解析として行った後、初めて脈波ＦＦＴ解析Ａ２を行う時には、ｃｙ＝２である。また例えば、脈波ＦＦＴ解析Ａ１を１回目の脈波ＦＦＴ解析として行った後に脈波ＦＦＴ解析Ａ２を３０回繰り返し実行し、その後、初めて脈波ＦＦＴ解析Ａ３を行う時には、ｃｙ＝３２である。
【００４５】
ステップＳ３２に続くステップＳ３３において、演算部３２は、周波数Ｆ［ｃｙ］を、“ＨＲ［ｙ］＝Ｆ［ｃｙ］×６０”に従ってハートレート値ＨＲ［ｃｙ］に変換し、Ｆ［ｃｙ］及びＨＲ［ｃｙ］をメモリ３３に保存する。ハートレート値は心拍数の検出値を指し、ハートレート値の単位は、心拍数と同様、“ｂｐｍ（beat per minute）”である。
【００４６】
――脈波ＦＦＴ解析Ａ２（図９）――
図９を参照し、脈波ＦＦＴ解析Ａ２について説明する。脈波ＦＦＴ解析Ａ２は、ステップＳ４１〜Ｓ４４の処理から成る。ステップＳ４１において、演算部３２は、窓関数処理を経た最新の脈波データ（即ち、脈波データ中の解析対象データ列）に対してＦＦＴを適用することで、脈波信号のパワースペクトルを算出する。
【００４７】
続くステップＳ４２において、演算部３２は、前回の脈波ＦＦＴ解析にて特定された周波数Ｆ［ｃｙ−１］を元に検索対象帯域を設定する。周波数Ｆ［ｃｙ］の導出に用いる検索対象帯域は、周波数Ｆ［ｃｙ−１］を内包し且つ所定の帯域幅を有する。ここでは、図１０に示す如く、周波数“Ｆ［ｃｙ−１］−７Ｒ”及び“Ｆ［ｃｙ−１］＋７Ｒ”を、夫々、検索対象帯域の下限及び上限に設定するものとする（Ｒは周波数分解能を表す；図５の（ｄ）欄参照）。そして、演算部３２は、最新の脈波データによるパワースペクトルを参照し、周波数Ｆ［ｃｙ−１］を基準に設定した検索対象帯域内で、ＰＳＤ値が最大となる要素帯域の周波数をＰＳＤ最大周波数Ｆ［ｃｙ］として検索及び特定する。つまり、周波数Ｆ［ｃｙ−１］を基準に帯域制限をかけた状態でＰＳＤ最大周波数Ｆ［ｃｙ］を検索する。
【００４８】
ステップＳ４２に続くステップＳ４３において、演算部３２は、周波数Ｆ［ｃｙ］を、“ＨＲ［ｙ］＝Ｆ［ｃｙ］×６０”に従ってハートレート値ＨＲ［ｃｙ］に変換し、Ｆ［ｃｙ］及びＨＲ［ｃｙ］をメモリ３３に保存する。更にその後、ステップＳ４４において、演算部３２は、ハートレート値の移動平均値ＨＲAVE［ｃｙ］を求めてメモリ３３に保存する。移動平均値ＨＲAVE［ｃｙ］は、直近の計ＡＶＥＮＵＭMAX個のハートレート値ＨＲ［ｃｙ−ＡＶＥＮＵＭMAX＋１］〜ＨＲ［ｃｙ］を移動平均することで得られる。移動平均は、単純移動平均でも良いし、加重移動平均でも良い。ＡＶＥＮＵＭMAXは２以上の任意の整数である。移動平均の対象の個数がＡＶＥＮＵＭMAX未満である場合には、ハートレート値ＨＲ［１］〜ＨＲ［ｃｙ］を移動平均することで移動平均値ＨＲAVE［ｃｙ］を得ればよい。
【００４９】
――脈波ＦＦＴ解析Ａ３（図１１）――
図１１を参照し、脈波ＦＦＴ解析Ａ３について説明する。脈波ＦＦＴ解析Ａ３は、ステップＳ５１〜Ｓ５５の処理から成る。ステップＳ５１において、演算部３２は、窓関数処理を経た最新の脈波データ（即ち、脈波データ中の解析対象データ列）に対してＦＦＴを適用することで、脈波信号のパワースペクトルを算出する。
【００５０】
ステップＳ５２において、演算部３２は、ステップＳ５１で算出された脈波信号のパワースペクトルに対し減衰補正を施すことで補正パワースペクトル（換言すれば、補正脈波信号のパワースペクトル）を導出する。説明の明確化のため、減衰補正前のパワースペクトルを原パワースペクトルと呼ぶ。この減衰補正では、原パワースペクトルにおいて、所定の周波数ＦTH未満の各要素帯域のＰＳＤ値を所定の減衰率（例えば１／１００）にて減衰させる。この減衰補正の様子を図１２の（ａ）欄に示す。パワースペクトル３１０及び３１０ｃは、夫々、原パワースペクトル及び補正パワースペクトルの例である（図１２の（ａ）欄の説明は後にも設けられる）。
【００５１】
続くステップＳ５３において、演算部３２は、前回の脈波ＦＦＴ解析にて特定された周波数Ｆ［ｃｙ−１］を元に検索対象帯域を設定する。検索対象帯域の設定方法は、上述したものと同様である。そして、演算部３２は、最新の脈波データに基づく補正パワースペクトルにおいて、周波数Ｆ［ｃｙ−１］を基準に設定した検索対象帯域内で、ＰＳＤ値が最大となる要素帯域の周波数をＰＳＤ最大周波数Ｆ［ｃｙ］として検索及び特定する。
【００５２】
その後、ステップＳ５４において、演算部３２は、周波数Ｆ［ｃｙ］を、“ＨＲ［ｙ］＝Ｆ［ｃｙ］×６０”に従ってハートレート値ＨＲ［ｃｙ］に変換し、Ｆ［ｃｙ］及びＨＲ［ｃｙ］をメモリ３３に保存する。更にその後、ステップＳ５５において、演算部３２は、ハートレート値の移動平均値ＨＲAVE［ｃｙ］を求めてメモリ３３に保存する。移動平均値ＨＲAVE［ｃｙ］の導出方法は上述した通りである。
【００５３】
制御部１３は、通信部１５を用いて、表示部１４及び外部機器ＥＥに脈波情報を出力することができ、順次導出されるハートレート値の移動平均値ＨＲAVE［ｃｙ］を出力脈波情報に含めることができる。通常は、最新の移動平均値ＨＲAVE［ｃｙ］を、生体２の心拍数の検出値として表示部１４及び外部機器ＥＥに出力すれば良い。但し、移動平均前のハートレート値ＨＲ［ｃｙ］や周波数Ｆ［ｃｙ］を出力脈波情報に含めることも可能である。ＨＲ［ｃｙ］もＨＲAVE［ｃｙ］も、脈波センサ１にて導出される心拍数の検出値である。尚、ＨＲAVE［ｃｙ］が導出されたとき、過去に算出された一部又は全部の移動平均値をメモリ３３から削除するようにしても良い。例えば、ＨＲAVE［ｃｙ］が導出されたとき、ＨＲAVE［ｃｙ］を含む、直近１分間分のハートレート値の移動平均値をメモリ３３に残存させておく一方で、それより過去に導出されたハートレート値の移動平均値をメモリ３３から削除しても良い。
【００５４】
［減衰補正の意義］
このように、演算部３２は、光センサ部１１の測定により取得された原脈波信号（Ｐａｏ、Ｐａ又はＰｄ）の内、所定周波数ＦTH未満の信号成分を減衰補正によって減衰させることで補正脈波信号を生成し、補正脈波信号に基づいて心拍数（ＨＲ［ｃｙ］又はＨＲAVE［ｃｙ］）を導出する。原脈波信号は、減衰補正前の脈波信号を指す。当該減衰補正を、時間領域上の原脈波信号に対して実行することも可能ではあるが、上述の動作例では、当該減衰補正を周波数領域上で実行している。即ち、原パワースペクトルに対して当該減衰補正を行うことで、補正脈波信号のパワースペクトルである補正パワースペクトルを導出している（図１１のステップＳ５２参照）。
【００５５】
周波数ＦTHは、生体２の体動の周期の期待値に基づき設定される。ステップＳ５２の減衰補正は、脈波信号に混入している体動の信号成分の低減を目指すものである。一方、生体２の体動（体の動き）の周期は、一般に、短くても６０／８０秒〜６０／１００秒程度である。故に、周波数ＦTHを、例えば８０／６０Ｈｚ〜１００／６０Ｈｚの間の周波数から選ぶと良く、典型的には１．５Ｈｚとしても良い。１．５Ｈｚの周波数ＦTHをｂｐｍの単位で表すと９０ｂｐｍである。減衰補正における減衰率として、例えば、１／１０〜１／１００の範囲内の数値を選ぶと良い。
【００５６】
既に参照した図１２の（ａ）欄のパワースペクトル３１０は、被験者が、ジョギング等の活動を行っている状態で観測される原パワースペクトルの例である。原パワースペクトル３１０及び補正パワースペクトル３１０ｃの夫々において、ＰＳＤ値のピーク（極大値）ＰＫN及びＰＫBは、夫々、体動による信号成分及び真の脈波の信号成分によるものである。活動中における真の心拍数は値（ＦTH×６０）を超えることが高確度で期待されるため、減衰補正によって、真の心拍数検出にとってノイズとなるピークＰＫNの値だけが効果的に減衰されることになる。つまり、上述の減衰補正により、生体２が活動していても精度良く心拍数を導出することが可能になる。
【００５７】
一方、図１２の（ｂ）欄のパワースペクトル３２０及び３２０ｃは、夫々、被験者が静止している状態における原パワースペクトル及び補正パワースペクトルの例である。被験者が静止しているときも、真の心拍数によるＰＳＤ値が原パワースペクトル３２０においてピークＰＫBとして現れる。また、被験者が静止しておれば、ピークＰＫBの周波数が周波数ＦTH未満となって上記減衰補正によりピークＰＫBの値が減衰されることもある。しかしながら、被験者静止時には、体動によるピークＰＫNが存在しないため（つまりノイズが十分に小さいため）、ピークＰＫBの値が減衰されたとしても、正しく心拍数を検出することができる。
【００５８】
尚、上述の動作例では、減衰補正を脈波ＦＦＴ解析Ａ３においてのみ実行しているが、当該減衰補正を脈波ＦＦＴ解析Ａ２においても実行可能である（この場合、図６のステップＳ１７にて脈波ＦＦＴ解析Ａ３が実行されることになる）。
【００５９】
［検索対象帯域の設定（帯域制限）の意義］
また、演算部３２は、光センサ部１１の測定で順次取得される原脈波信号の信号値を用いて心拍数の順次導出を行うが、各回の心拍数を、帯域制限を適用した状態で導出する。帯域制限では、今回の心拍数（ＨＲ［ｃｙ］）を、前回導出した心拍数（ＨＲ［ｃｙ−１］）を内包し且つ所定の大きさを持った数値範囲内で検索し導出する。上述の動作例では、当該探索を“Ｈｚ”を単位とした周波数領域上で実行しており、図１０の検索対象帯域の単位を“ｂｐｍ”に変換したものが当該数値範囲に相当する（図９のステップＳ４２、図１１のステップＳ５３も参照）。
【００６０】
より具体的には、今回の心拍数（ＨＲ［ｃｙ］）を導出するための原脈波信号のパワースペクトル（即ち、ＨＲ［ｃｙ］の導出に用いる解析対象データ列のパワースペクトル）から補正パワースペクトルを生成した後、当該補正パワースペクトルにおいて上記数値範囲内に対応する検索対象帯域を設定し、検索対象帯域内でＰＳＤ値が最大となる周波数をＦ［ｃｙ］として同定することにより今回の心拍数（ＨＲ［ｃｙ］）を導出する。生体２の心拍数の短時間における変化量は限られており、今回の検出心拍数が前回の検出心拍数と大きくかけ離れているとは考えにくい。そこで、上記の如く検索対象帯域内で最大ＰＳＤ値を探索するようにする。これにより、心拍数検出値（即ちハートレート値）がノイズに基づく異常値になることが抑制される。換言すれば、心拍数検出に対するノイズの影響を低減することができる。
【００６１】
検索対象帯域の大きさは任意である。１秒間における人間の心拍数の変化量の最大値は１５ｂｐｍ程度と考えられるため、脈波ＦＦＴ解析の実行間隔が１秒相当である場合には、今回の心拍数（ＨＲ［ｃｙ］）を算出する際の検索対象帯域を、例えば、ｂｐｍを単位として、“ＨＲ［ｃｙ−１］−１５”以上且つ“ＨＲ［ｃｙ−１］＋１５”以下の帯域に設定すると良い。
【００６２】
［実験結果など］
図１３及び図１４に、被験者を用いた実験結果の例を示す。波形３４１〜３４３を含むグラフ３４０（図１３の（ａ）欄）、波形３５１〜３５３を含むグラフ３５０（図１３の（ｂ）欄）、波形３６１〜３６３を含むグラフ３６０（図１４の（ａ）欄）、波形３７１〜３７３を含むグラフ３７０（図１４の（ｂ）欄）は、夫々、実験Ｅ1A、Ｅ1B、Ｅ1C、Ｅ1Dの結果を表している。各実験では、被験者が静止している２つの区間中に一定区間だけ被験者を歩行又は走行させた。被験者の歩行又は走行速度は、実験Ｅ1A、Ｅ1B、Ｅ1C、Ｅ1Dにおいて、夫々、時速６、８、１０、１２ｋｍである。破線波形３４１、３５１、３６１、３７１は、市販の胸ベルト装着型心拍計を用いた、真値とみなせる心拍数の測定結果を表す。実線波形３４２、３５２、３６２、３７２は、脈波センサ１による心拍数の測定結果（ＨＲAVE［ｃｙ］の時系列データ）を表す。一点鎖線波形３４３、３５３、３６３、３７３は、従来のピーク解析方法を用いたときの測定結果である。尚、グラフ３４０において、破線波形３４１及び実線波形３４２を区別して示すべく、それらを実際の波形位置から若干ずらして示している。グラフ３５０〜３７０においても同様であり、後述のグラフ４４０〜４７０（図２３及び図２４）においても同様である。
【００６３】
時速６ｋｍの歩行時には、従来のピーク解析方法でも胸ベルト装着型心拍計と相関のある測定結果が得られた。しかし、時速８ｋｍ以上の走行時において、従来のピーク解析方法による測定結果は、体動ノイズの影響を受けて真値（胸ベルト装着型心拍計の測定結果）から大きくかけ離れている。被験者が安静にしている状態又はそれに近い状態では、体動等によるノイズが十分に小さい脈波信号波形が得られるため、ピーク解析方法でも、正確な心拍数検出が可能である。但し、被験者の活動時など、体動等によるノイズが脈波信号波形に多く含まれるようになると、取得波形が乱れて、図１３及び図１４に示す如く正確な心拍数検出が難しくなる。これに対し、第１実施形態の脈波センサ１を用いると、様々な活動環境下で、胸ベルト装着型心拍計と相関のある測定結果が得られていることが分かる。
【００６４】
尚、上述の動作例では、脈波ＦＦＴ解析の実行回数ｃｙが１回目のときにだけ脈波ＦＦＴ解析Ａ１を実行しているが、脈波ＦＦＴ解析Ａ１を複数回実行するようにしても良い。例えば、図６において、回数ｃｙが１、２及び３回目である時の夫々においてステップＳ１５からステップＳ１６に進み、回数ｃｙが４回目以降である状態でステップＳ１５に至ったときにステップＳ１７へ進むようにしても良い。
【００６５】
この場合、１〜３回目のＰＳＤ最大周波数Ｆ［１］〜Ｆ［３］は脈波ＦＦＴ解析Ａ１で求められ、４回目のＰＳＤ最大周波数Ｆ［４］は脈波ＦＦＴ解析Ａ２で求められる。この際、周波数Ｆ［４］を求めるときの検索対象帯域の中心を、周波数Ｆ［１］〜Ｆ［３］に基づき定めればよい。例えば、Ｆ［１］〜Ｆ［３］の平均又は加重平均を、周波数Ｆ［４］の導出用の検索対象帯域の中心に設定しても良い。或いは例えば、Ｆ［２］とＦ［３］が互いに一致又は近似している一方で、Ｆ［１］とＦ［２］又はＦ［３］とが大きく相違している場合には、Ｆ［１］の誤差が大きいと判断して、Ｆ［２］及びＦ［３］の平均を、周波数Ｆ［４］の導出用の検索対象帯域の中心に設定しても良い。本実施形態では、検索対象帯域の設定による帯域制限を介して各回の心拍数を検出するため、一度、不適切な検索対象帯域を設定してしまうと、心拍数の真値検出が困難な状況が長時間継続する。故に、帯域制限をかけない脈波ＦＦＴ解析Ａ１を複数回実施し、複数回の脈波ＦＦＴ解析Ａ１の結果を統計的に処理して検索対象帯域を設定することにより、検索対象帯域の信頼性が高まる。
【００６６】
＜＜第２実施形態＞＞
本発明の第２実施形態を説明する。第２実施形態及び後述の第３実施形態は第１実施形態を基礎とする実施形態であり、第２及び第３実施形態において特に述べない事項に関しては、矛盾の無い限り、第１実施形態の記載が第２及び第３実施形態にも適用される。矛盾の無い限り、第１〜第３実施形態の内、任意の複数の実施形態を組み合わせても良い。第２実施形態では、上述の脈波センサ１による他の心拍数導出方法を説明する。第１及び第２実施形態間で脈波センサ１の構成は同じである。
【００６７】
第２実施形態では、体動センサ部１７が有益に機能する（図３参照）。体動センサ部１７は、生体２の動き（即ち体動）を検出することで、生体２の動きに応じた体動信号を生成及び出力する。ここでは、体動センサ部１７が、生体２の加速度を検出する加速度センサ（例えばサンプリングレート：８Ｈｚ)から成るものとする。故に、以下では、体動信号を加速度信号と呼ぶ。加速度信号は生体２の加速度を示す。体動センサ部１７内の加速度センサは、互いに直交するＸ、Ｙ及びＺ軸方向の生体２の加速度を個別に検出できても良い。この場合、体動センサ部１７は、Ｘ軸、Ｙ軸、Ｚ軸方向における生体２の加速度の検出結果を示すＸ軸、Ｙ軸、Ｚ軸方向の加速度信号を生成及び出力する。但し、体動センサ部１７によって検出及び出力される、生体２の動きを示す物理量は、加速度以外の物理量（速度、角速度、角加速度など）でもよい。
【００６８】
図１５の（ａ）欄及び（ｂ）欄には、生体２と、本体ユニット１０及びベルト２０と、Ｘ、Ｙ及びＺ軸との関係を表している。上述したように、脈波センサ１は、腕輪構造を有し、生体２の手首にベルト２０を巻回して固定され、これに伴って、手首における手の甲側の面又は手の平側の面に、本体ユニット１０の裏面である装着面が接触固定される。Ｘ軸は、生体２の前腕の長手方向（即ち、生体２の前腕が伸びる方向）に平行である。生体２の手首にベルト２０が巻回固定されるのであるから、ベルト２０と加速度センサとの位置関係から脈波センサ１においてＸ軸方向を特定できる。上記装着面はＸ軸及びＹ軸に平行であって且つＺ軸に直交する。
【００６９】
体動センサ部１７からは、アナログの加速度信号Ｍａが出力され、図１６に示す如く、制御部１３内のＡ／Ｄ変換部３６にてアナログの加速度信号Ｍａがデジタルの加速度信号Ｍｄに変換される。Ａ／Ｄ変換部３６は、体動センサ部１７内に備えられていても良い。体動センサ部１７が、複数の軸方向における複数の加速度信号Ｍａを出力する際、Ａ／Ｄ変換部３６は、複数の加速度信号Ｍａの夫々を加速度信号Ｍｄに変換する。体動センサ部１７の出力信号にフィルタ処理及び増幅処理を施した信号が加速度信号Ｍａであっても良い。
【００７０】
第２実施形態では、後述されるように加速度信号から生体２の動き状態を検出して心拍数の検出精度向上を図る。一方、生体２がジョギング等をするとき、脈波センサ１を装着した前腕の振動の主成分はＸ軸に平行であるため、Ｘ軸方向の加速度信号に体動の主成分が含まれる。このため、体動センサ部１７は、少なくともＸ軸方向の加速度信号を取得することが好ましい。生体２がジョギング等の活動をする際、前腕振動に対応するＸ軸方向の加速度を検出することで生体２の動き状態を良好に把握できるからである。勿論、体動センサ部１７は、Ｘ軸、Ｙ軸、Ｚ軸方向の加速度信号の内、２つ以上の加速度信号を取得しても良いし、それら全ての加速度信号を取得しても良い。以下の説明において、加速度信号とは、体動センサ部１７がＸ軸方向の加速度信号のみを取得する場合にあってはＸ軸方向の加速度信号を指し、体動センサ部１７がＸ、Ｙ軸、Ｚ軸方向の加速度信号の内の２つ以上の加速度信号を取得する場合にあっては当該２つ以上の加速度信号の合成ベクトル信号を指す。
【００７１】
制御部１３は、脈波信号だけでなく、加速度信号に対してもＦＦＴ解析を行う。加速度信号に対するＦＦＴ解析を加速度ＦＦＴ解析と呼ぶ。加速度ＦＦＴ解析の解析条件にも、脈波ＦＦＴ解析と同様に、サンプリング周波数、解析データ数、実行間隔、周波数分解能、解析下限周波数ｆLL及び解析上限周波数ｆLLが含まれるが（図５の（ａ）欄〜（ｄ）欄を参照）、加速度ＦＦＴ解析の解析条件は脈波ＦＦＴ解析の解析条件と同じとされる。第１実施形態で述べた脈波ＦＦＴ解析の解析条件の説明文が加速度ＦＦＴ解析の解析条件にも適用され、当該適用の際、第１実施形態に記載の脈波ＦＦＴ解析、Ａ／Ｄ変換部３１、脈波信号、脈波データ、記号“Ｐｄ”が、夫々、加速度ＦＦＴ解析、Ａ／Ｄ変換部３６、加速度信号、加速度データ、記号“Ｍｄ”に読み替えられる。従って例えば、時刻ｔ＝ｉのサンプリングにて取得される加速度信号Ｍｄのデジタル値は“Ｍｄ［ｉ］”にて表される。つまり、時刻ｔ＝ｉにて、第ｉ番目の加速度信号値、即ちＭｄ［ｉ］が得られる。Ｐｄ［ｉ］にて値が表される加速度信号を、加速度データとも呼ぶ。
【００７２】
［心拍数の導出方法］
図１７は、制御部１３による、心拍数を導出するためのフローチャートである。まず、ステップＳ１１１において、Ａ／Ｄ変換部３１及び３６は、所定のサンプリング周波数にて脈波データ及び加速度データの取り込み（即ちサンプリング）を開始する。順次得られる脈波データ及び加速度データはメモリ３３に保存される。ステップＳ１１１に続くステップＳ１１２において、制御部１３は、これまでに取り込まれた脈波データの個数ｎを、所定の第１閾値（ＰＲA−Ｑ）、及び、解析データ数と一致する第２閾値ＰＲAと比較する（図５の（ｂ）欄を参照）。第１不等式“ＰＲA−Ｑ≦ｎ＜ＰＲA”の成立時にはステップＳ１１２からステップＳ１１３に進み、第２不等式“ＰＲA≦ｎ”の成立時にはステップＳ１１２からステップＳ１１９に進む。不等式“ｎ＜ＰＲA−Ｑ”の成立時には、データ数が少なすぎるため、第１不等式が成立するまで脈波データ及び加速度データの取り込みを継続しつつステップＳ１１２で待機する。
【００７３】
ステップＳ１１３の段階で取得されている脈波データは、Ｐｄ［１］〜Ｐｄ［ｎ］であり、それらの総個数は解析データ数ＰＲAより少ない。故に、時刻ｔにて実行されるステップＳ１１３において、演算部３２は、時刻ｔ＝ｎ＋１〜時刻ｔ＝ＰＲAにおける脈波データＰｄ［ｎ＋１］〜Ｐｄ［ＰＲA］として、所定の補完用データを補完する。この補完の方法は第１実施形態のそれと同様である。ステップＳ１１３では、実測によって得られた脈波データＰｄ［１］〜Ｐｄ［ｎ］と、補完された脈波データＰｄ［ｎ＋１］〜Ｐｄ［ＰＲA］とから成る解析対象データ列が形成される。
【００７４】
ステップＳ１１３に続くステップＳ１１４において、演算部３２は、解析対象データ列が存在する区間だけ０以外の値を持つ窓関数（例えばハニング関数）を脈波データに対して掛け合わせる窓関数処理を実行する。その後、ステップＳ１１５において、演算部３２は、これから行うべき脈波ＦＦＴ解析が、何回目の脈波ＦＦＴ解析であるのかを確認する。そして、これから行うべき脈波ＦＦＴ解析が、１回目の脈波ＦＦＴ解析である場合には、ステップＳ１１６に進んで脈波ＦＦＴ解析Ｂ１を実行し、２回目以降の脈波ＦＦＴ解析である場合には、ステップＳ１１７に進んで脈波ＦＦＴ解析Ｂ２を実行する。脈波ＦＦＴ解析Ｂ１及びＢ２は、ステップＳ１１４の窓関数処理を経た脈波データ（即ち、脈波データ中の解析対象データ列）に対して実行される。
【００７５】
ステップＳ１１６の脈波ＦＦＴ解析Ｂ１又はステップＳ１１７の脈波ＦＦＴ解析Ｂ２を終えると、ステップＳ１１８へ進む。ステップＳ１１８において、Ａ／Ｄ変換部３１及び３６により、新たにＰＲB個分の脈波データ及びＰＲB個分の加速度データが取り込まれ、ステップＳ１１２に戻る。新たに取り込まれた脈波データ及び加速度データはメモリ３３に保存される。
【００７６】
一方、“ＰＲA≦ｎ”の場合、既に取得されている脈波データの総個数は解析データ数ＰＲA以上であるので、上述のような補完は不要である。ステップＳ１１９において、演算部３２は窓関数処理を実行する。ステップＳ１１９の窓関数処理では、解析対象データ列が存在する区間だけ０以外の値を持つ窓関数（例えばハニング関数）を脈波データに対して掛け合わせる処理に加えて、同様の窓関数を加速度データに対しても掛け合わせる。時刻ｔ＝ｎにおいては、脈波信号の解析対象データ列が脈波データＰｄ［ｎ−ＰＲA＋１］〜Ｐｄ［ｎ］から成り（図５の（ｂ）欄を参照）、このとき、加速度データＭｄ［ｎ−ＰＲA＋１］〜Ｍｄ［ｎ］が存在する区間だけ０以外の値を持つ窓関数が加速度データに掛け合わされる。
【００７７】
その後、ステップＳ１２０において、演算部３２は、ステップＳ１１９の窓関数処理を経た加速度データに対し加速度ＦＦＴ解析を実行し、その結果を利用して体動モード判定及び補正係数設定を行う（詳細は後述）。そして、ステップＳ１２１において、演算部３２は、ステップＳ１２０の結果を利用しつつ、ステップＳ１１９の窓関数処理を経た脈波データ（即ち、脈波データ中の解析対象データ列）に対し脈波ＦＦＴ解析Ｂ３を実行する（詳細は後述）。
【００７８】
ステップＳ１２１の脈波ＦＦＴ解析Ｂ３を終えると、ステップＳ１２２へ進む。ステップＳ１２２において、Ａ／Ｄ変換部３１及び３６により、新たにＰＲB個分の脈波データ及びＰＲB個分の加速度データが取り込まれた後、ステップＳ１１９に戻る。新たに取り込まれた脈波データ及び加速度データはメモリ３３に保存される。この際、メモリ３３に保存されている脈波データの内、より古くから保存されているＰＲB個分の脈波データをメモリ３３から削除しても良い。加速度データについても同様である。
【００７９】
脈波ＦＦＴ解析の解析条件は、脈波ＦＦＴ解析Ｂ１、Ｂ２及びＢ３間で共通である。脈波ＦＦＴ解析Ｂ１、Ｂ２及びＢ３において、周波数ｆLL及びｆHHは、夫々、例えば、０．６３５Ｈｚ及び４．０Ｈｚである。
【００８０】
――脈波ＦＦＴ解析Ｂ１、Ｂ２――
脈波ＦＦＴ解析Ｂ１、Ｂ２は、夫々、第１実施形態における脈波ＦＦＴ解析Ａ１、Ａ２と同じものである。尚、第２実施形態においては、“ｃｙ”は、脈波ＦＦＴ解析Ｂ１、Ｂ２及びＢ３の全実行回数を表すものと解される。
【００８１】
――加速度ＦＦＴ解析を介した体動モード判定＆補正係数設定――
加速度ＦＦＴ解析を伴うステップＳ１２０の処理を説明する。図１８を参照する。ステップＳ１２０の処理は、ステップＳ１３１〜Ｓ１３４の処理を含む。まずステップＳ１３１において、演算部３２は、窓関数処理を経た加速度データに対してＦＦＴを適用することで、加速度信号のパワースペクトルを算出する。加速度信号のパワースペクトルは、加速度信号の第１〜第ｍ要素帯域の夫々におけるＰＳＤ値を情報として含んでいる。続くステップＳ１３２において、演算部３２は、ステップＳ１３１で求めた複数のＰＳＤ値の総和ＰＳＤSUMを求める。総和ＰＳＤSUMは、加速度信号の第１〜第ｍ要素帯域における計ｍ個のＰＳＤ値の総和であっても良いし、加速度信号の特定帯域における複数のＰＳＤ値の総和であっても良い。特定帯域は、第１〜第ｍ要素帯域の全帯域の一部である。生体２の動き方にも依存するが、生体２の動きが激しくなるほど総和ＰＳＤSUMは増大する。
【００８２】
そしてステップＳ１３３において、演算部３２は、総和ＰＳＤSUMに基づき体動モードｍｏｄｅを設定する。具体的には、総和ＰＳＤSUMが所定の閾値ＴＨPSDSUMより小さいとき、演算部３２は、生体２が静止している又は生体２の動きが小さいと判断して“ｍｏｄｅ”に１を設定し、総和ＰＳＤSUMが当該閾値ＴＨPSDSUM以上であるとき、生体２が活動していると判断して“ｍｏｄｅ”に２を設定する。
【００８３】
一方、ステップＳ１３４において、演算部３２は、要素帯域ごとに加速度データのＰＳＤ値を所定の閾値ＴＨPSDと比較し、要素帯域ごとに、ＰＳＤ値が閾値ＴＨPSD以上なら補正係数に第１補正値を設定する一方でＰＳＤ値が閾値ＴＨPSD未満なら補正係数に第２補正値を設定する。第１及び第２補正値は共に任意の正の所定値であるが、第１補正値は第２補正値よりも小さい。第２補正値は１より大きくても構わない。例えば、第２補正値が１．０であるとき、第１補正値を、０．１〜０．０１の範囲内から選べば良い。ここでは、第１補正値は０．１であり、第２補正値は１．０であるとする。第ｉ要素帯域に対して設定された補正係数を“ＫＣi”にて表す。
【００８４】
尚、補正係数の設定を複数の要素帯域ごとに行うようにしても良い。例えば、２つの要素帯域ごとに補正係数の設定を行う場合、加速度データにおける第ｉ及び第（ｉ＋１）要素帯域のＰＳＤ値の合計を所定の閾値ＴＨPSDと比較し、その合計が閾値ＴＨPSD以上なら補正係数ＫＣi及びＫＣi+1に第１補正値（０．１）を設定し、そうでないなら、補正係数ＫＣi及びＫＣi+1に第２補正値（１．０）を設定すれば良い。
【００８５】
――脈波ＦＦＴ解析Ｂ３（図１９）――
図１９を参照し、脈波ＦＦＴ解析Ｂ３について説明する。脈波ＦＦＴ解析Ｂ３は、ステップＳ１４１〜Ｓ１４６の処理から成る。ステップＳ１４１において、演算部３２は、窓関数処理を経た最新の脈波データ（即ち、脈波データ中の解析対象データ列）に対してＦＦＴを適用することで、脈波信号のパワースペクトルを算出する。
【００８６】
次にステップＳ１４２において、演算部３２は、ステップＳ１３４（図１８参照）で設定された補正係数ＫＣ1〜ＫＣmを用いて、ステップＳ１４１で算出したパワースペクトルを補正することで、補正パワースペクトルを算出する。第１実施形態と同様、脈波信号に関する補正前のパワースペクトルを原パワースペクトルと呼ぶ。但し、第２実施形態における補正パワースペクトルは、第１実施形態と異なり、補正係数に基づき算出されたものを指す。補正係数に基づく補正では、要素帯域ごとに、対応する補正係数を原パワースペクトルのＰＳＤ値に乗じることで補正パワースペクトルを得る。つまり、任意の整数ｉに関し、補正パワースペクトルの第ｉ要素帯域のＰＳＤ値は、原パワースペクトルの第ｉ要素帯域のＰＳＤ値と補正係数ＫＣiとの積である。このような補正により、体動の加速度が大きい帯域、即ち、体動ノイズが多く混入していると予想される帯域において、脈波信号のＰＳＤ値が減衰せしめられる。この減衰補正の様子を図２０に示す。パワースペクトル４１０は加速度信号のパワースペクトルの例であり、パワースペクトル４２０及び４２０ｃは、夫々、原パワースペクトル及び補正パワースペクトルの例である（図２０を参照した説明は後にも設けられる）。
【００８７】
次にステップＳ１４３において、演算部３２は、体動モードｍｏｄｅに基づき検索対象帯域を設定する。検索対象帯域の意義は第１実施形態で述べた通りであり、周波数Ｆ［ｃｙ］の導出に用いる検索対象帯域は、周波数Ｆ［ｃｙ−１］を内包し且つ所定の帯域幅を有する。ｍｏｄｅが“１”である時と“２”である時とで、検索対象帯域は互いに異なる。具体的には、図２１の（ａ）欄及び（ｂ）欄に示す如く、“ｍｏｄｅ”が静止モードに対応する“１”であるとき、周波数“Ｆ［ｃｙ−１］−７Ｒ”及び“Ｆ［ｃｙ−１］＋７Ｒ”を、夫々、検索対象帯域の下限及び上限に設定する一方、“ｍｏｄｅ”が活動モードに相当する“２”であるとき、周波数“Ｆ［ｃｙ−１］−１Ｒ”及び“Ｆ［ｃｙ−１］＋１３Ｒ”を、夫々、検索対象帯域の下限及び上限に設定する。このような数値は、勿論例示であり、様々に変更できる。但し、ｍｏｄｅ＝２であるとき、周波数Ｆ［ｃｙ］を求めるための検索対象帯域の中心は、周波数Ｆ［ｃｙ−１］よりも大きい。ｍｏｄｅ＝１であるとき、周波数Ｆ［ｃｙ］を求めるための検索対象帯域の中心は、周波数Ｆ［ｃｙ−１］と一致していて良い。
【００８８】
続くステップＳ１４４にて、演算部３２は、最新の脈波データに基づく補正パワースペクトルにおいて、周波数Ｆ［ｃｙ−１］を基準に設定した検索対象帯域内で、ＰＳＤ値が最大となる要素帯域の周波数をＰＳＤ最大周波数Ｆ［ｃｙ］として検索及び特定する。
【００８９】
その後、ステップＳ１４５において、演算部３２は、周波数Ｆ［ｃｙ］を、“ＨＲ［ｙ］＝Ｆ［ｃｙ］×６０”に従ってハートレート値ＨＲ［ｃｙ］に変換し、Ｆ［ｃｙ］及びＨＲ［ｃｙ］をメモリ３３に保存する。更にその後、ステップＳ１４６において、演算部３２は、ハートレート値の移動平均値ＨＲAVE［ｃｙ］を求めてメモリ３３に保存する。移動平均値ＨＲAVE［ｃｙ］の導出方法は上述した通りである。
【００９０】
［体動モード判定に応じた、検索対象帯域可変設定の意義］
第１実施形態と同様、第２実施形態でも、各回の心拍数を、帯域制限を適用した状態で導出する。つまり、今回の心拍数（ＨＲ［ｃｙ］）を、前回導出した心拍数（ＨＲ［ｃｙ−１］）を内包し且つ所定の大きさを持った数値範囲内で検索し導出する。上述の動作例では、当該探索を“Ｈｚ”を単位とした周波数領域上で実行しており、図１９の検索対象帯域の単位を“ｂｐｍ”に変換したものが当該数値範囲に相当する。より具体的には、今回の心拍数（ＨＲ［ｃｙ］）を導出するための脈波信号に基づくパワースペクトルにおいて、上記数値範囲内に対応する検索対象帯域を設定し、検索対象帯域内でＰＳＤ値が最大となる周波数をＦ［ｃｙ］として同定することにより今回の心拍数（ＨＲ［ｃｙ］）を導出する。生体２の心拍数の短時間における変化量は限られており、今回の心拍数が前回の心拍数と大きくかけ離れているとは考えにくい。そこで、上記の如く検索対象帯域内で最大ＰＳＤ値を探索するようにする。これにより、心拍数導出に対するノイズの影響を低減することができる。
【００９１】
但し、第２実施形態では、上記数値範囲（即ち検索対象帯域）を、加速度信号に基づく体動モードｍｏｄｅの設定を介して、動的に可変設定する。“ｍｏｄｅ”を１に設定することは、生体２の動き状態が動き状態ＭＳAに属すると分類判定することに相当し、“ｍｏｄｅ”を２に設定することは、生体２の動き状態が動き状態ＭＳBに属すると分類判定することに相当する。動き状態ＭＳAは、生体２の動きが無い状態を含み、生体２が比較的小さく動いている状態を含みうる。動き状態ＭＳBは、生体２が動いている状態を含み、少なくとも、動き状態ＭＳBにおける生体２の動きの大きさは動き状態ＭＳAにおける生体２の動きの大きさよりも大きい。生体２の動きの大きさは、例えば、上記の総和ＰＳＤSUMとして評価される。
【００９２】
そして、今回の心拍数（ＨＲ［ｃｙ］）を導出する際に動き状態が動き状態ＭＳBであると判定されたとき、即ち“ｍｏｄｅ＝２”であるとき、今回の心拍数（ＨＲ［ｃｙ］）を導出するための検索対象帯域の中心を、前回導出した心拍数の周波数（Ｆ［ｃｙ−１］）よりも大きくするオフセット処理を行う（図２１の（ｂ）欄を参照）。“ｍｏｄｅ＝２”であるとき、生体２の活動により心拍数が増加傾向にあることが予想されるため、検索対象帯域を高域側にシフトすることにより、より妥当な心拍数検出が見込める（真の心拍数の周波数が検索対象帯域外に存在するような状況の発生が回避される）。尚、上述のオフセット処理は、前回導出した心拍数（ＨＲ［ｃｙ−１］）が所定心拍数以下の場合に限って実施されるものであっても良い。この所定心拍数は例えば１２０ｂｐｍである（勿論、１２０ｂｐｍ以外にすることも可能である）。
【００９３】
一方、“ｍｏｄｅ＝１”であるとき、心拍数の変動は少ないことが予想されるため、今回の心拍数（ＨＲ［ｃｙ］）を導出するための検索対象帯域の中心を、前回導出した心拍数の周波数（Ｆ［ｃｙ−１］）と一致させる良い（図２１の（ａ）欄を参照）。
【００９４】
但し、第１時刻において動き状態が動き状態ＭＳBであると判定され、その後の第２時刻において動き状態が動き状態ＭＳBから動き状態ＭＳAに切り替わったとき、演算部３２は、第２時刻を起点とし且つ所定の時間長を有する特定期間を設けるようにしても良い。特定期間では、今回の心拍数（ＨＲ［ｃｙ］）を導出するための検索対象帯域の中心が、前回導出した心拍数の周波数（Ｆ［ｃｙ−１］）よりも小さくされる。例えば、特定期間においては、周波数“Ｆ［ｃｙ−１］−８Ｒ”及び“Ｆ［ｃｙ−１］＋４Ｒ”を、夫々、検索対象帯域の下限及び上限に設定して良い。特定期間中は、生体２の活動停止に伴って心拍数が減少傾向にあることが予想されるためであり、検索対象帯域を低域側にシフトすることにより、より妥当な心拍数検出が見込める。特定期間は、第２時刻における心拍数ＨＲ［ｃｙ］が所定値（例えば、１１０ｂｐｍ）以上である場合に限って設定されるものであっても良い。心拍数ＨＲ［ｃｙ］が比較的少ない状態において体動が停止したとしても（即ち例えば、ほんの少し運動してから止まったような状態では）、時間経過に伴う心拍数低下は殆ど生じないことが予想されるからである。
【００９５】
尚、生体２の動き状態が動き状態ＭＳBに属するとき、生体２の動き状態を更に複数の動き状態に細分化しても良い。例えば、演算部３２は、“ＴＨPSDSUM＜ＰＳＤSUM≦ＴＨ２PSDSUM”の成立時には、生体２の動き状態が動き状態ＭＳB1に属すると判定してｍｏｄｅに“２”を設定する一方で、“ＴＨ２PSDSUM＜ＰＳＤSUM”の成立時には、生体２の動き状態が動き状態ＭＳB2に属すると判定してｍｏｄｅに“３”を設定する（但し、所定の閾値ＴＨ２PSDSUMは閾値ＴＨPSDSUMより大きい）。動き状態ＭＳB1及びＭＳB2は共に動き状態ＭＳBに属するが、動き状態ＭＳB2の方が動き状態ＭＳB1よりも生体の動きの大きさが大きい。そして、“ｍｏｄｅ＝３”であるとき、今回の心拍数（ＨＲ［ｃｙ］）を導出するための検索対象帯域の中心を、“ｍｏｄｅ＝２”であるときのそれよりも高くすると良い。
【００９６】
［帯域制限の一時停止］
検索対象帯域を用いた帯域制限は心拍数の検出精度向上に寄与するものの、一度、不適切な検索対象帯域を設定してしまうと、心拍数の真値検出が困難な状況が長時間継続するおそれもある。不適切な検索対象帯域の設定は生体２の活動中に起こりやすい。一方、生体２の静止中には、図１２の（ｂ）欄に示すような体動ノイズの少ない脈波信号が得られるため、帯域制限を設けなくても、各回の脈波信号に基づき十分に信頼性の高い心拍数検出が可能である。
【００９７】
そこで、第１時刻において動き状態が動き状態ＭＳBであると判定され、その後の第２時刻において動き状態が動き状態ＭＳBから動き状態ＭＳAに切り替わったとき、図２２に示す如く、演算部３２は、第２時刻を起点とし且つ所定の時間長を有する停止期間を設けるようにしても良い（但し、第１時刻では帯域制限が実行されていたものとする）。停止期間では、上述の帯域制限の実行が停止され、図８の脈波ＦＦＴ解析Ａ１と同じ脈波ＦＦＴ解析Ｂ１を実行することで、ＰＳＤ最大周波数Ｆ［ｃｙ］及びハートレート値ＨＲ［ｃｙ］を求める（つまり、検索対象帯域に依存せずにＦ［ｃｙ］及びＨＲ［ｃｙ］が求められる）。そして、停止期間の終了後は、停止期間中に求めたＰＳＤ最大周波数を用いて、検索対象帯域の設定を介した帯域制限を再開する。これにより、生体２の活動停止を契機に検索対象帯域が見直され、不適切な検索対象帯域の設定が長時間継続することが抑制される。
【００９８】
単純には例えば、停止期間の時間長が脈波ＦＦＴ解析の実行間隔（図５の（ｃ）欄を参照）と同じである場合において、停止期間中にＰＳＤ最大周波数Ｆ［ｃｙ−１］として周波数ＦSTOP1が求められた場合、停止期間終了直後の検索対象帯域の中心を周波数ＦSTOP1に設定し、その設定状態で脈波ＦＦＴ解析Ｂ２又はＢ３の繰り返し実行を再開すれば良い。停止期間の時間長が脈波ＦＦＴ解析の実行間隔の２倍以上である場合、停止期間中に求めた複数のＰＳＤ最大周波数に基づき停止期間終了直後の検索対象帯域の中心を設定し、その設定状態で脈波ＦＦＴ解析Ｂ２又はＢ３の繰り返し実行を再開すれば良い。例えば、複数のＰＳＤ最大周波数の平均を、停止期間終了直後の検索対象帯域の中心に設定して良い。
【００９９】
或いは例えば、以下のような方法を用いても良い。即ち、停止期間直前に求められたＰＳＤ最大周波数が周波数ＦA0である場合において、停止期間中に求められた１回目のＰＳＤ最大周波数が周波数ＦA1であるとき、演算部３２は、周波数ＦA0及びＦA1の差ＤＩＦ01を所定の差分閾値と比較する。そして、その差ＤＩＦ01が差分閾値未満であるならば、停止期間を直ちに終了して、周波数ＦA0若しくはＦA1又はそれらの平均値を停止期間終了直後の検索対象帯域の中心に設定し、その設定状態で脈波ＦＦＴ解析Ｂ２又はＢ３の繰り返し実行を再開する。差ＤＩＦ01が差分閾値以上であるならば、停止期間を終了させずに、停止期間中における２回目のＰＳＤ最大周波数を周波数ＦA2として求めて、周波数ＦA1及びＦA2の差ＤＩＦ12を所定の差分閾値と比較する。差ＤＩＦ12が差分閾値未満であるならば、停止期間を直ちに終了して、周波数ＦA1若しくはＦA2又はそれらの平均値を停止期間終了直後の検索対象帯域の中心に設定し、その設定状態で脈波ＦＦＴ解析Ｂ２又はＢ３の繰り返し実行を再開する。差ＤＩＦ12が差分閾値以上になることは考えにくいが、仮に、差ＤＩＦ12が差分閾値以上であったならば、停止期間中に連続して求めた２つのＰＳＤ最大周波数の差又は停止期間中に連続して求めた２以上のＰＳＤ最大周波数の最大差が差分閾値以下になるまで停止期間を継続し、差分閾値以下の差又は最大差が得られた時点で停止期間を終了すると良い。この場合、停止期間の終了直前に求められた１つのＰＳＤ最大周波数、又は、停止期間の終了直前に求められた２以上のＰＳＤ最大周波数の平均値を、停止期間終了直後の検索対象帯域の中心に設定し、その設定状態で脈波ＦＦＴ解析Ｂ２又はＢ３の繰り返し実行を再開すれば良い。
【０１００】
［補正係数によるパワースペクトル補正の意義］
また、第２実施形態では、加速度信号のパワースペクトルにおいて、ＰＳＤ値が閾値ＴＨPSD以上になっている要素帯域に対して第２補正値を持つ補正係数が設定される（図１８のステップＳ１３４参照）。加速度信号のパワースペクトルにおいて、ＰＳＤ値が閾値ＴＨPSD以上になっている各要素帯域は減衰対象帯域に属すると考えられる。演算部３２は、図１８のステップＳ１３４にて減衰対象帯域を特定した後、脈波信号の原パワースペクトルにおいて減衰対象帯域内のＰＳＤ値を減衰させることで補正パワースペクトルを得ている（図１９のステップＳ１４２参照）。そして、補正パワースペクトルにおいて、上記数値範囲内に対応する検索対象帯域を設定し、検索対象帯域内でＰＳＤ値が最大となる周波数をＦ［ｃｙ］として同定することにより今回の心拍数（ＨＲ［ｃｙ］）を導出する。
【０１０１】
減衰対象帯域は体動の加速度が大きい帯域、即ち、脈波信号に体動ノイズが多く混入していると予想される帯域である。上記のような補正係数を用いた減衰補正により、図２０に示す如く、脈波信号に混入した体動ノイズが減少せしめられて補正パワースペクトル４２０ｃが得られる。補正パワースペクトル４２０ｃを用いてＰＳＤ最大周波数Ｆ［ｃｙ］を同定するようにすれば、体動ノイズの影響が軽減された精度の高い心拍数検出が可能である。
【０１０２】
尚、加速度信号の各要素帯域のＰＳＤ値に応じて補正係数を３段階以上に設定しても良い。つまり例えば、図１８のステップＳ１３４において、演算部３２は、要素帯域ごとに、ＰＳＤ値が所定閾値ＴＨPSD以上なら補正係数に所定値ＶＡＬ1（例えば０．１）を設定し、ＰＳＤ値が閾値ＴＨPSD未満であるが所定閾値ＴＨ２PSD以上なら補正係数に所定値ＶＡＬ2（例えば０．５）を設定し、ＰＳＤ値が閾値ＴＨ２PSD未満なら補正係数に所定値ＶＡＬ3（例えば１．０）を設定するようにしても良い。ここで、ＴＨPSD＞ＴＨ２PSD、且つ、ＶＡＬ1＜ＶＡＬ2＜ＶＡＬ3、である。
【０１０３】
また、上述の動作例では、加速度データを脈波ＦＦＴ解析Ｂ３においてのみ利用しているが、加速度データを脈波ＦＦＴ解析Ｂ２にも利用するようにしても良い。この場合、図１７のステップＳ１１３及びＳ１１４にて、脈波データと同様の補完処理及び窓関数処理を加速度データに対しても実行し、ステップＳ１１７にて、脈波ＦＦＴ解析Ｂ２の代わりにステップＳ１２０及びＳ１２１の処理を実行すればよい。
【０１０４】
［実験結果など］
図２３及び図２４に、被験者を用いた実験結果の例を示す。波形４４１〜４４３を含むグラフ４４０（図２３の（ａ）欄）、波形４５１〜４５３を含むグラフ４５０（図２３の（ｂ）欄）、波形４６１〜４６３を含むグラフ４６０（図２４の（ａ）欄）、波形４７１〜４７３を含むグラフ４７０（図２４の（ｂ）欄）は、夫々、実験Ｅ2A、Ｅ2B、Ｅ2C、Ｅ2Dの結果を表している。各実験では、被験者が静止している２つの区間中に一定区間だけ被験者を歩行又は走行させた。被験者の歩行又は走行速度は、実験Ｅ2A、Ｅ2B、Ｅ2C、Ｅ2Dにおいて、夫々、時速６、８、１０、１２ｋｍである。破線波形４４１、４５１、４６１、４７１は、市販の胸ベルト装着型心拍計を用いた、真値とみなせる心拍数の測定結果を表す。実線波形４４２、４５２、４６２、４７２は、脈波センサ１による心拍数の測定結果（ＨＲAVE［ｃｙ］の時系列データ）を表す。一点鎖線波形４４３、４５３、４６３、４７３は、従来のピーク解析方法を用いたときの測定結果である。
【０１０５】
図２３及び図２４に対応する実験では、歩行時においてすら、従来のピーク解析方法による測定結果が、体動ノイズの影響を受けて真値（胸ベルト装着型心拍計の測定結果）から大きくかけ離れている。これに対し、第２実施形態の脈波センサ１を用いると、様々な活動環境下で、胸ベルト装着型心拍計と相関のある測定結果が得られている。第１実施形態の方法を用いた場合、被験者によっては、脈波センサ１の測定心拍数と胸ベルト装着型心拍計との相関が薄れる可能性があるが、第２実施形態の方法を用いれば被験者に依存せず当該相関をより強固に確保できる。
【０１０６】
尚、上述の動作例では、脈波ＦＦＴ解析の実行回数ｃｙが１回目のときにだけ脈波ＦＦＴ解析Ｂ１を実行しているが、脈波ＦＦＴ解析Ｂ１を複数回実行するようにしても良い。例えば、図１７において、回数ｃｙが１、２及び３回目である時の夫々においてステップＳ１１５からステップＳ１１６に進み、回数ｃｙが４回目以降である状態でステップＳ１１５に至ったときにステップＳ１１７へ進むようにしても良い。この場合、１〜３回目のＰＳＤ最大周波数Ｆ［１］〜Ｆ［３］は脈波ＦＦＴ解析Ｂ１で求められ、４回目のＰＳＤ最大周波数Ｆ［４］は脈波ＦＦＴ解析Ｂ２で求められる。この際、周波数Ｆ［４］を求めるときの検索対象帯域の中心を、周波数Ｆ［１］〜Ｆ［３］に基づき定めればよい（当該定め方は第１実施形態で述べたものと同様）。
【０１０７】
加速度ＦＦＴ解析の解析条件は全て脈波ＦＦＴ解析の解析条件と同じであって良いが、それらの解析条件の内、一部（例えば、サンプリング周波数）は互いに異なっていても良い。但し、それらの間で周波数分解能が異なっていると、要素帯域ごとの補正係数による補正（図１９のステップＳ１４２）を行い難くなる。故に、加速度ＦＦＴ解析の解析条件と脈波ＦＦＴ解析の解析条件との間で、少なくとも周波数分解能は互いに同じにしておくと良い。周知の如く、ＦＦＴは、時間領域上の信号（即ち時間軸上のデータ列）を周波数領域上の信号に変換する時間／周波数変換処理の一種である。本発明において、時間／周波数変換処理は、ＦＦＴに属するものでなくても構わない。
【０１０８】
上述したように、Ａ／Ｄ変換部３６は加速度センサに内包されるものであって良い。この場合において、加速度センサ内のＡ／Ｄ変換部３６のサンプリング周波数が加速度ＦＦＴ解析の解析条件におけるサンプリング周波数より大きい場合、演算部３２は、公知の周波数変換処理（間引き等）を利用して加速度ＦＦＴ解析用の加速度データを得れば良い。
【０１０９】
＜＜第３実施形態＞＞
本発明の第３実施形態を説明する。第３実施形態で述べる技術は第１及び第２実施形態に適用される。
【０１１０】
発光部１１Ａは、生体２に対して光を照射する発光素子（ＬＥＤ等）を有し、発光部１１Ａに含まれる発光素子の個数は１つでも良いし、２以上でも良い。受光部１１Ｂは、発光部１１Ａが照射した光の内、生体２内を透過した光を受光する受光素子（フォトダイオード又はフォトトランジスタ）を有し、受光部１１Ｂに含まれる受光素子の個数は１つでも良いし、２以上でも良い。
【０１１１】
望ましい光センサ部１１の形態の１つを、図２５に示す。図２５の光センサ部１１は、生体２と接触する装着面５０と、装着面５０に固定された１つの発光素子５１及び４つの受光素子６１〜６４と、を有する。発光素子５１は、緑色の光を出力するＬＥＤであって発光部１１Ａを形成する。受光素子６１〜６４はフォトトランジスタであって受光部１１Ｂを形成する。発光素子５１は、受光素子６１及び６３の中心に位置すると共に、受光素子６２及び６４の中心に位置する。受光素子６１及び６３の各中心を結ぶ線と、受光素子６２及び６４の各中心を結ぶ線は、互いに直交する。発光素子５１の中心と受光素子６１〜６４の中心との間の距離は互いに等しい。受光素子６１〜６４の夫々は、発光部１１Ａからの光に基づき生体２を透過又は反射した光を受光するものであって良い。
【０１１２】
受光素子６１が受光した光の量に応じた電流が受光素子６１に流れる。受光素子６２〜６４についても同様である。このとき、光センサ部１１は、受光素子６１〜６４に流れる電流の合計電流に基づき脈波信号を形成できる。このため、光センサ部１１では、発光素子及び受光素子が夫々１つの光センサ部に比べて、信号強度を高めることができる。図２５の構成から受光素子６２及び６４を削除することもできるし、図２５の構成に対して更に他の受光素子を追加することもできる。何れにせよ、受光部１１Ｂにおける複数の受光素子の間に、発光素子５１を配置することが好ましい。
【０１１３】
＜＜第４実施形態＞＞
先述の第２実施形態では、加速度信号に基づくパワースペクトル（以下では、加速度パワースペクトルと略称する）を算出し、その算出結果に応じて脈波信号のパワースペクトル（以下では、脈波パワースペクトルと略称する）を補正すると共に、心拍数ＨＲの検索対象帯域（帯域制限範囲）を可変制御することにより、安静時だけでなく活動時においても心拍数ＨＲを精度良く導出することができる旨を提案した。
【０１１４】
なお、第２実施形態のアルゴリズムは、歩行時、ジョギング時、或いは、ランニング時など、脈波センサ部１１が装着されている腕を比較的大きく振るような活動時（腕振運動時）の心拍数を検出する手法として有効である。
【０１１５】
一方、本願の発明者は、さらなる鋭意研究の結果、リストロール時、リストカール時、キータイピンク時、或いは、エルゴメータ使用時など、あまり腕を振らない活動時（非腕振運動時）には、第２実施形態のアルゴリズムによる心拍数検出精度がやや低下するという知見を得た。
【０１１６】
そこで、第４実施形態では、あまり腕を振らない活動時においても、心拍数ＨＲを精度良く導出することのできる解析手法について提案する。
【０１１７】
［フローチャート］
図２６は、本発明の第４実施形態に係る解析動作の概略（全体像）を説明するためのフローチャートである。本実施形態の解析動作においても、先の第２実施形態と同じく、光センサ部１１から得られる時系列データ（脈波データ）と体動センサ部１７（例えば３軸加速度センサ）から得られる時系列データ（加速度データ）を用いて、安静時及び活動時の心拍数ＨＲを算出するアルゴリズムが採用されている。なお、本フローの動作主体は、基本的に演算部３２である。
【０１１８】
フローが開始されると、ステップＳ２０１において、脈波データ及び加速度データの蓄積が行われる。ここで蓄積された脈波データ及び加速度データは、ステップＳ２０２におけるＦＦＴ解析に供される。また、加速度データは、ステップＳ２０６におけるベクトル値Ｖ及び振幅変化総和Ｄの算出処理にも供される。なお、ステップＳ２０２とステップＳ２０６は、それぞれ並列に処理が進められる。
【０１１９】
ステップＳ２０２では、同一の条件で脈波データ及び加速度データのＦＦＴ解析が行われて、脈波パワースペクトルと加速度パワースペクトルが各々算出される。ＦＦＴ解析については、先に説明した通りであるので重複した説明は割愛する。ここで算出された脈波パワースペクトルは、ステップＳ２０５における補正脈波パワースペクトルの生成処理に供される。一方、加速度パワースペクトルは、ステップＳ２０３におけるピーク位置検出処理に供される。
【０１２０】
ステップＳ２０３では、加速度パワースペクトルのピーク位置が検出される（詳細は後述）。その後、フローはステップＳ２０４に進められる。
【０１２１】
ステップＳ２０４では、加速度パワースペクトルのピーク位置に応じて体動ノイズ成分を減衰するための補正データが生成される（詳細は後述）。ここで生成された補正データは、ステップＳ２０５における補正脈波パワースペクトルの生成処理に供される。
【０１２２】
ステップＳ２０５では、脈波パワースペクトルに対する補正データの乗算処理が行われて補正脈波パワースペクトルが生成される（詳細は後述）。ここで生成された補正脈波パワースペクトルは、ステップＳ２０８における帯域制限処理に供される。
【０１２３】
一方、ステップＳ２０６では、加速度データのベクトル値Ｖと振幅変化総和Ｄが算出される（詳細は後述）。その後、フローはステップＳ２０７に進められる。
【０１２４】
ステップＳ２０７では、加速度データのベクトル値Ｖと振幅変化総和Ｄ（ないしは前回の心拍数ＨＲ）に応じて体動モードが判定される（詳細は後述）。ここでの判定結果は、ステップＳ２０８における帯域制限処理に供される。
【０１２５】
ステップＳ２０８では、補正脈波パワースペクトルに対して体動モードに応じた帯域制限が掛けられる（詳細は後述）。その後、フローはステップＳ２０９に進められる。
【０１２６】
ステップＳ２０９では、帯域制限の掛けられた補正脈波パワースペクトルから心拍数ＨＲが導出される（詳細は後述）。
【０１２７】
以下では、上記各ステップの動作について個別具体的に説明していく。
【０１２８】
［ピーク位置検出に基づく体動ノイズキャンセル動作］
図２７は、加速度パワースペクトルのピーク位置検出に基づく体動ノイズキャンセル動作（図２６のステップＳ２０３〜Ｓ２０５に相当）の一例を示す図である。なお、（ａ）欄には加速度パワースペクトル、（ｂ）欄には体動ノイズ周波数帯域を設定するための帯域フラグα、（ｃ）欄には脈波パワースペクトルの減衰率を設定するための補正係数β、（ｄ）欄には脈波パワースペクトルに乗算される補正データγ、及び、（ｅ）欄には脈波パワースペクトル（破線）と補正脈波パワースペクトル（実線）が各々描写されている。
【０１２９】
（ａ）欄で示すように、演算部３２は、加速度パワースペクトル（体動パワースペクトルに相当）のピーク位置を解析する。ピーク位置の解析手法としては、隣接周波数間で加速度パワースペクトルの差分値を順次算出し、その変極点（差分値が正から負に切り替わる極大点）を探索すればよい。上記の差分演算に先立ち、加速度パワースペクトルの平滑化処理を行うことも可能である。なお、ピーク位置に相当する周波数には「１」、ピーク位置に相当しない周波数には「０」がフラグ付けされる。
【０１３０】
本図の例では、加速度パワースペクトルに７つのピーク位置が検出されている。このように、多数のピーク位置が検出された場合、それら全てを脈波パワースペクトルの補正に反映させると、脈波信号を不必要に減衰させてしまい、却って心拍数ＨＲの検出精度を悪化させてしまうおそれがある。
【０１３１】
そこで、演算部３２は、複数検出されたピーク位置のうちパワースペクトル密度の上位順に所定数のピーク位置（本図の例では、丸印を付した３つのピーク位置）だけを抽出して脈波パワースペクトルの補正に反映させる。このような選定処理を行うことにより、多数のピーク位置が検出された場合でも、脈波パワースペクトルを適切に補正することが可能となる。
【０１３２】
ピーク位置の選定後、演算部３２は、ピーク位置を中心とした前後の周波数帯域を体動ノイズ周波数帯域として設定する。具体的には、（ｂ）欄で示すように、演算部３２は、選定した３つのピーク位置を中心に各々±ｘＲ（例えばｘ＝３、５、７、９）の周波数範囲で帯域フラグαを「１」とし、その余の周波数範囲で帯域フラグαを「０」とする。
【０１３３】
先述の第２実施形態では、加速度パワースペクトル密度が所定の閾値よりも高い周波数帯域を体動ノイズ周波数帯域として設定し、同帯域内に属する脈波パワースペクトルを一律的に減衰させていた。そのため、加速度パワースペクトルの分布状態（ピーク幅）によって体動ノイズ周波数帯域が広くなったり狭くなったりするので、脈波パワースペクトルを適切に補正することができないおそれがあった。
【０１３４】
一方、本実施形態であれば、体動ノイズ周波数帯域を常に一定幅（±ｘＲ）に固定することができる。従って、加速度パワースペクトルの分布状態（ピーク幅）に依らず、脈波パワースペクトルを適切に補正することが可能となる。
【０１３５】
体動ノイズ周波数帯域の設定後、演算部３２は、帯域フラグαと補正係数βを用いて補正データγを生成する。具体的には、（ｄ）欄で示したように、体動ノイズ周波数帯域内（α＝１）ではγ＝β（＝α×β）となり、体動ノイズ周波数帯域外（α＝０）ではγ＝１（減衰なし）に固定される。
【０１３６】
また、本実施形態では、体動ノイズ周波数帯域に属する脈波パワースペクトルを一律的に減衰するのではなく、（ｃ）欄で示したように、脈波パワースペクトルの減衰率を設定するための係数βが周波数に応じた変動値として用意されている。
【０１３７】
本図の例では、周波数ｆ１以下の周波数範囲で補正係数βが「β１」に固定されて、周波数ｆ１〜ｆ２の周波数範囲で補正係数βが「β１」から「β２」まで線形的に変化されて、周波数ｆ２以上の周波数範囲で補正係数βが「β２」に固定される（例えば、ｆ１＝１２０ｂｐｍ相当、ｆ２＝１８０ｂｐｍ相当、β１＝０．０１、β２＝０．１）。なお、補正係数βの設定例についてはこれに限定されるものではない。
【０１３８】
このように、補正係数βを周波数に応じて可変的に設定することにより、体動ノイズが重畳しやすい周波数帯域と、体動ノイズが重畳しにくい周波数帯域のそれぞれにおいて、脈波パワースペクトルの減衰量を最適に調整することが可能となる。
【０１３９】
補正データの生成後、演算部３２は、脈波パワースペクトルに補正データγを乗算する（体動ノイズ周波数帯域に属する脈波パワースペクトルの各周波数成分に対して各々に対応した補正係数βを乗算する）ことにより、補正脈波パワースペクトルを生成する。
【０１４０】
上記一連の体動ノイズキャンセル動作により、（ｅ）欄で示したように、体動ノイズ成分（破線）を取り除いた補正脈波パワースペクトルが得られるので、安静時だけでなく活動時においても心拍数ＨＲを精度良く導出することが可能となる。
【０１４１】
［加速度データのベクトル値Ｖと振幅変化総和Ｄに基づく体動モード判定動作］
次に、加速度データのベクトル値Ｖと振幅変化総和Ｄに基づく体動モード判定動作（図２６のステップＳ２０６〜Ｓ２０７に相当）について詳述する。
【０１４２】
図２８は、ベクトル値Ｖと振幅変化総和Ｄの定義式を示す図である。時刻ｔに体動センサ部１７で得られる３軸加速度データを各々Ｘ（ｔ）、Ｙ（ｔ）、Ｚ（ｔ）とした場合、加速度データのベクトル値Ｖ及び振幅変化総和Ｄ（＝各軸における振幅変化の総和）は、それぞれ、次の（１）式及び（２）式により算出される。
【０１４３】
【数１】
2015129557.tif 000003
【０１４４】
【数２】
2015129557.tif 000004
【０１４５】
先述の第２実施形態では、加速度パワースペクトルの総和が所定値を超えているか否かに応じて体動モードを２パターンに分類し、体動モードに応じた帯域制限処理を行っていた（図１８のステップＳ１３３や図１９のステップＳ１４３を参照）。
【０１４６】
一方、本実施形態では、上記（１）式及び（２）式を用いて算出される加速度データのベクトル値Ｖ及び振幅変化総和Ｄ（ないしは前回の心拍数ＨＲ）に応じて体動モードが６パターンに分類されており、体動モードに応じた帯域制限処理がよりきめ細かく実施される。以下、図面を参照しながら詳細に説明する。
【０１４７】
図２９は、体動モードの一例を示すテーブルであり、左から順に、体動モード、帯域制限（前回心拍数ＨＲ±分解能）、想定活動状態、閾値設定（ベクトル値Ｖ、振幅変化総和Ｄ、及び、前回心拍数ＨＲ）、並びに、備考が記載されている。なお、静止中の被験者に対して１Ｇが掛かっている状態では、Ｖ＝１０００が算出されるものとする。
【０１４８】
体動モードとしては、ｍｏｄｅ０〜ｍｏｄｅ５の６パターンが用意されている。ｍｏｄｅ０は、心拍数ＨＲの計測開始直後に設定される体動モードである。ｍｏｄｅ０では、被験者が連続静止状態であるという想定の下、帯域制限が解除されて全帯域が心拍数ＨＲの検索対象とされる。なお、Ｖ＜１０５０、かつ、Ｄ＜２０が満たされているときには、ｍｏｄｅ０と判定される（前回心拍数ＨＲは不問）。ただし、帯域制限の解除後に心拍数ＨＲが導出された時点で体動モードはｍｏｄｅ１に切り替えられる。
【０１４９】
ｍｏｄｅ１では、被験者が軽微な動きを行っているという想定の下、心拍数ＨＲの検出対象として、前回心拍数ＨＲ±３Ｒの帯域制限が掛けられる。ｍｏｄｅ１での帯域制限については、図２１の（ａ）欄の±７Ｒを±３Ｒとして読み替えて参照することにより、容易に理解することができる。なお、Ｖ＜１０５０、Ｄ＜２０、かつ、ＨＲ＜１２０が満たされているときには、ｍｏｄｅ１と判定される。
【０１５０】
ただし、ｍｏｄｅ１が所定時間（例えば１６秒間）に亘って連続すると、体動モードがｍｏｄｅ０に切り替えられる。言い換えると、被験者が安静状態となる度に、体動モードがｍｏｄｅ０に戻されて帯域制限が解除される。このような構成とすることにより、何らかの不具合により不適切な帯域制限が掛けられた場合であっても、被験者が安静状態に戻れば、誤った帯域制限も自動的に解除されるので、運動再開時には心拍数ＨＲを正しく検出することが可能となる。
【０１５１】
ｍｏｄｅ２では、被験者が活動を開始したという想定の下、心拍数ＨＲの検出対象として、前回心拍数ＨＲ±７Ｒ（オフセット＋６Ｒあり）の帯域制限が掛けられる。ｍｏｄｅ２での帯域制限については、図２１の（ｂ）欄をそのまま参照することにより、容易に理解することができる。なお、１１５０＜Ｖ、かつ、ＨＲ＜１２０が満たされているときには、ｍｏｄｅ２と判定される（振幅変化総和Ｄは不問）。ただし、ｍｏｄｅ２は、心拍数ＨＲの急上昇が想定される活動開始直後の４秒間にのみ有効であり、当該期間の経過後は他の体動モードに移行される。
【０１５２】
ｍｏｄｅ３では、被験者が活動中（腕振運動中）であるという想定の下、心拍数ＨＲの検出対象として、前回心拍数ＨＲ±５Ｒの帯域制限が掛けられる。ｍｏｄｅ３での帯域制限については、図２１の（ａ）欄の±７Ｒを±５Ｒとして読み替えて参照することにより、容易に理解することができる。なお、１０５０＜Ｖ、または、１２０＜ＨＲが満たされているときには、ｍｏｄｅ３と判定される（振幅変化総和Ｄは不問）。
【０１５３】
ｍｏｄｅ４では、被験者が活動中（手首運動などの非腕振運動中）であるという想定の下、当該体動モードへの突入時における心拍数ＨＲ±５Ｒの帯域制限が掛けられる。すなわち、ｍｏｄｅ１〜ｍｏｄｅ３の帯域制限中心周波数は、いずれも前回心拍数ＨＲを基準とした可変値であるのに対して、ｍｏｄｅ４の帯域制限中心周波数は、ｍｏｄｅ４突入時の心拍数ＨＲを基準とした固定値とされている。このような体動モードを用意しておくことにより、腕振運動時だけでなく非腕振運動時においても精度良く心拍数ＨＲを導出することが可能となる。なお、Ｖ＜１０５０、かつ、８０＜Ｄが満たされているときには、ｍｏｄｅ４と判定される（前回心拍数ＨＲは不問）。
【０１５４】
ｍｏｄｅ５では、被験者が軽微な動き（クールダウン）を行っているという想定の下、心拍数ＨＲの検出対象として、前回心拍数ＨＲ±５Ｒの帯域制限が掛けられる。ｍｏｄｅ５での帯域制限については、先のｍｏｄｅ３と同様、図２１の（ａ）欄の±７Ｒを±５Ｒとして読み替えて参照することにより、容易に理解することができる。なお、Ｖ＜１０５０、かつ、２０＜Ｄ＜８０が満たされているときには、ｍｏｄｅ５と判定される（前回心拍数ＨＲは不問）。ｍｏｄｅ５の帯域制限中心周波数は、ｍｏｄｅ１〜ｍｏｄｅ３と同様、前回心拍数ＨＲを基準とした可変値である。
【０１５５】
上記の体動モード切替制御（特にｍｏｄｅ４やｍｏｄｅ５の創設）は、非腕振運動時にはベクトル値Ｖがあまり変化せずに振幅変化総和Ｄが比較的大きく変化する、という本願発明者の新たな知見に基づいて初めて成し得るものである。
【０１５６】
なお、本図の例では、加速度データのベクトル値Ｖと振幅変化総和Ｄ（ないしは前回心拍数ＨＲ）に応じて、６パターンの体動モードを切り替える説明を行ったが、体動モードのパターン数についてはこれに限定されるものではなく適宜増減することも可能である。また、体動モードの切替に際して加速度データのベクトル方向を併せて参照してもよい。
【０１５７】
［周波数平均値に基づく心拍数ＨＲの導出アルゴリズム］
次に、帯域制限の掛けられた補正脈波パワースペクトルから心拍数ＨＲを導出するための新規アルゴリズム（図２６のステップＳ２０９に相当）について詳述する。
【０１５８】
先述の第２実施形態では、補正脈波パワースペクトルのパワースペクトル密度が最大となる周波数（以下ではＰＳＤ最大周波数と呼ぶ）を検索し、そのＰＳＤ最大周波数を心拍数ＨＲに変換していた。もちろん、上記アルゴリズムは、図２６のステップＳ２０９にも適用することが可能である。しかしながら、帯域制限範囲内に補正脈波パワースペクトルのピークが複数存在する場合には、上記アルゴリズムによる心拍数ＨＲの導出に不具合を生じるおそれがあった。
【０１５９】
図３０は、ＰＳＤ最大周波数に基づくＨＲ導出時の問題点を説明するための図である。帯域制限範囲に補正脈波パワースペクトルのピークが複数存在していた場合であっても、（ａ）欄で示すように、最大ピークのパワースペクトル密度が突出しているのであれば、ＰＳＤ最大周波数に基づく心拍数ＨＲの導出に特段の問題は生じない。
【０１６０】
しかしながら、（ｂ）欄で示すように、同程度のパワースペクトル密度を持つ複数のピークが互いに離れた周波数に現れている場合には、いずれか一方のピークを選択してＰＳＤ最大周波数としても、心拍数ＨＲを正しく導出することができないおそれがある。
【０１６１】
また、（ｃ）欄で示すように、帯域制限範囲に突出したピークが現れない場合には、僅かなパワースペクトル密度差に基づいてＰＳＤ最大周波数を選択せざるを得ず、心拍数ＨＲの導出精度が悪化してしまう。
【０１６２】
そこで、本願発明者は、上記の問題点に鑑み、心拍数ＨＲの導出アルゴリズムを全面的に刷新するという決断を下した。より具体的に述べると、演算部３２は、帯域制限の掛けられた補正脈波パワースペクトルから心拍数ＨＲを導出する際、補正脈波パワースペクトルについて、帯域制限範囲内で各周波数に対するパワースペクトル密度で重み付けされた周波数平均値Ａａｖｅを算出し、その算出結果から心拍数ＨＲを導出する。
【０１６３】
図３１は、周波数平均値Ａａｖｅの算出式を示す図であり、左から順に、補正脈波パワースペクトルの周波数成分Ａ（ｋ）、パワースペクトル密度Ｂ（ｋ）、及び、両者の乗算値Ｃ（ｋ）（＝Ａ（ｋ）×Ｂ（ｋ））が描写されている。なお、ｋ＝１、２、…、ｍであり、ｍは帯域制限幅±ｘＲに応じて変化する可変値（ｍ＝２ｘ＋１）である。例えば、ｘ＝３のときにはｍ＝７となる。
【０１６４】
ここで、周波数平均値Ａａｖｅは、次の（３）式により算出される。
【０１６５】
【数３】
2015129557.tif 000005
【０１６６】
この周波数平均値Ａａｖｅに基づいて心拍数ＨＲを導出することにより、補正脈波パワースペクトルに複数のピークが存在している場合であっても、心拍数ＨＲを正しく導出することが可能となる。
【０１６７】
［実験結果］
図３２は、本発明の第４実施形態の実験による心拍数測定結果を示す図である。なお、（ａ）欄には手首運動（静止→リストロール→静止→指運動→静止→リストカール→静止→歩行→静止）を行っているときの測定結果が示されており、（ｂ）欄にはエルゴメータ運動（静止→エルゴメータ使用→静止）を行っているときの測定結果が示されている。
【０１６８】
また、各欄の実線は第４実施形態の測定結果を示しており、一点鎖線は第２実施形態の測定結果を示している。一方、各欄の破線は、胸ベルト装着型心拍計の測定結果を比較参照値（正しい測定結果に相当）として示している。
【０１６９】
本図から明らかなように、第２実施形態では手首運動時やエルゴメータ運動時に測定された心拍数ＨＲが比較参照値から乖離してしまっている。一方、第４実施形態ではそのような乖離は生じておらず、手首運動時やエルゴメータ運動時においても心拍数ＨＲを正しく測定することができている。
【０１７０】
＜＜変形等＞＞
本発明の実施形態は、特許請求の範囲に示された技術的思想の範囲内において、適宜、種々の変更が可能である。以上の実施形態は、あくまでも、本発明の実施形態の例であって、本発明ないし各構成要件の用語の意義は、以上の実施形態に記載されたものに制限されるものではない。上述の説明文中に示した具体的な数値は、単なる例示であって、当然の如く、それらを様々な数値に変更することができる。上述の任意の式において、矛盾が生じない範囲で、不等号“≧”又は“≦”を夫々不等号“＞”又は“＜”に置き換えても良いし、その逆も可能である。上述の実施形態に適用可能な注釈事項として、以下に、注釈１〜注釈５を記す。各注釈に記載した内容は、矛盾なき限り、任意に組み合わせることが可能である。
【０１７１】
［注釈１］
上述の脈波センサ１は、脈波信号取得部、体動信号取得部及び演算部３２を有して、生体２の心拍数を検出及び導出する心拍数検出装置を内包している。心拍数の検出又は導出は、心拍数の検出値であるハートレート値の導出と同義である。脈波信号取得部は光センサ部１１を含む。フィルタ部１２及びＡ／Ｄ変換部３１も脈波信号取得部の構成要素に含まれる、と考えても良い。体動信号取得部は体動センサ部１７を含む。Ａ／Ｄ変換部３６、及び、Ａ／Ｄ変換部３６と体動センサ部１７との間に設けられうるフィルタ部（不図示）も、体動信号取得部の構成要素に含まれる、と考えても良い。
【０１７２】
［注釈２］
図４の演算部３２及びメモリ３３は、脈波センサ１ではなく、外部機器ＥＥ（図３参照）に設けられていても良い。つまり、演算部３２及びメモリ３３における上述の動作を外部機器ＥＥ内で実現しても良い。演算部３２及びメモリ３３における上述の動作の内、一部を脈波センサ１内で実行し、残りを外部機器ＥＥ内で実行しても良い。
【０１７３】
［注釈３］
発光部１１Ａにおける発光素子の発光は、常時発光でも良いし、パルス発光でも良い。パルス発光では、発光素子が発光している状態と発光素子が発光していない状態とが交互に訪れる。
【０１７４】
［注釈４］
脈波センサ１又は心拍数検出装置である対象装置を、集積回路等のハードウェア、或いは、ハードウェアとソフトウェアの組み合わせによって構成することができる。対象装置にて実現される機能の全部又は一部である任意の特定の機能をプログラムとして記述して、該プログラムを対象装置に搭載可能なフラッシュメモリに保存しておいても良い。そして、該プログラムをプログラム実行装置（例えば、対象装置に搭載可能なマイクロコンピュータ）上で実行することによって、その特定の機能を実現するようにしてもよい。上記プログラムは任意の記録媒体に記憶及び固定されうる。上記プログラムを記憶及び固定する記録媒体は対象装置と異なる機器（サーバ機器等）に搭載又は接続されても良い。
【０１７５】
［注釈５］
上述の各実施形態では、光センサ部１１を用いた光学的方式にて脈波信号を検出及び取得しているが、本発明はこれに限定されない。つまり本発明において、脈波を示す脈波信号は、生体２の血管内の圧力変化を検出する圧電方式で検出及び取得されたものでも良いし、心電図を計測できる心電計により検出及び取得されたものであっても良い（この場合、心電図の波形が脈波信号の波形に相当する）。
【符号の説明】
【０１７６】
１  脈波センサ
２  生体
１０  本体ユニット
１１  光センサ部
１１Ａ  発光部
１１Ｂ  受光部
１２  フィルタ部
１３  制御部
１７  体動センサ部
２０  ベルト
３１、３６  Ａ／Ｄ変換部
３２  演算部
３３  メモリ
５１  発光素子
６１〜６４  受光素子
【図１】
2015129557.tif 000006
【図２】
2015129557.tif 000007
【図３】
2015129557.tif 000008
【図４】
2015129557.tif 000009
【図５】
2015129557.tif 000010
【図６】
2015129557.tif 000011
【図７】
2015129557.tif 000012
【図８】
2015129557.tif 000013
【図９】
2015129557.tif 000014
【図１０】
2015129557.tif 000015
【図１１】
2015129557.tif 000016
【図１２】
2015129557.tif 000017
【図１３】
2015129557.tif 000018
【図１４】
2015129557.tif 000019
【図１５】
2015129557.tif 000020
【図１６】
2015129557.tif 000021
【図１７】
2015129557.tif 000022
【図１８】
2015129557.tif 000023
【図１９】
2015129557.tif 000024
【図２０】
2015129557.tif 000025
【図２１】
2015129557.tif 000026
【図２２】
2015129557.tif 000027
【図２３】
2015129557.tif 000028
【図２４】
2015129557.tif 000029
【図２５】
2015129557.tif 000030
【図２６】
2015129557.tif 000031
【図２７】
2015129557.tif 000032
【図２８】
2015129557.tif 000033
【図３０】
2015129557.tif 000034
【図３１】
2015129557.tif 000035
【図３２】
2015129557.tif 000036
【図２９】
2015129557.tif 000037
【国際調査報告】
2015129557.tif 000038
2015129557.tif 000039
2015129557.tif 000040
2015129557.tif 000041
2015129557.tif 000042
2015129557.tif 000043
2015129557.tif 000044
2015129557.tif 000045