(57)【特許請求の範囲】
【請求項１】
人体に装着して脈拍を計測する脈拍計において、
脈波センサを有し、脈波検出信号を出力する脈波検出部と、
加速度センサを有し、静脈血の挙動に影響を与える体動に対応する加速度を検出して体動検出信号として出力する体動検出部と、
前記体動検出信号に基づいて前記脈波検出信号に含まれる体動成分を除去する体動成分除去部と、
前記体動成分除去後の前記脈波検出信号に基づいて脈拍数を算出する脈拍数算出部と、を備え、
前記加速度センサは、装着者の腕の末梢方向に沿った軸をＸ軸とし、前記Ｘ軸に垂直で手の甲に垂直な軸をＺ軸とし、前記Ｘ軸および前記Ｚ軸に垂直な軸をＹ軸とし、前記Ｘ軸、前記Ｙ軸、前記Ｚ軸のそれぞれに沿った方向の加速度を検出する３軸加速度センサであり、
前記体動成分除去部は、前記体動検出信号に対応する、前記Ｘ軸に沿った方向の加速度成分であるＸ軸加速度成分、前記Ｙ軸に沿った方向の加速度成分であるＹ軸加速度成分および前記Ｚ軸に沿った方向の加速度成分であるＺ軸加速度成分に基づいて前記体動成分を生成する体動成分生成部を備えたことを特徴とする脈拍計。
【請求項２】
請求項１記載の脈拍計において、
前記加速度センサは、前記脈波センサの近傍に配置されていることを特徴とする脈拍計。
【請求項３】
請求項１記載の脈拍計において、
前記加速度センサは、前記脈波センサに略積層状態で配置されていることを特徴とする脈拍計。
【請求項４】
請求項１記載の脈拍計において、
前記体動成分除去部は、前記Ｙ軸加速度成分および前記Ｚ軸加速度成分をベクトルとみなし、これらのベクトルの合成ベクトルである２軸加速度合成成分と、前記Ｘ軸加速度成分と、に基づいて前記体動成分を生成する体動成分生成部を備えたことを特徴とする脈拍
計。
【請求項５】
請求項１記載の脈拍計において、
前記体動成分除去部は、前記Ｘ軸加速度成分、前記Ｙ軸加速度成分および前記Ｚ軸加速度成分の３つの加速度成分をベクトルとみなし、これらのベクトルの合成ベクトルである３軸加速度合成成分に基づいて前記体動成分を生成する体動成分生成部を備えたことを特徴とする脈拍計。
【請求項６】
請求項１記載の脈拍計において、
前記Ｘ軸加速度成分、前記Ｙ軸加速度成分および前記Ｚ軸加速度成分のうち少なくともいずれか一の加速度成分について重み付けを行うことを特徴とする脈拍計。
【請求項７】
請求項１記載の脈拍計において、
前記体動成分除去部は、前記Ｘ軸加速度成分、前記Ｙ軸加速度成分および前記Ｚ軸加速度成分に基づいて適応フィルタ係数を生成するフィルタ係数生成部と、
今回の前記脈波検出信号から前回の前記脈波検出信号に前記適応フィルタ係数を適用して得られる体動成分を除去する除去処理部と、
を備えたことを特徴とする脈拍計。
【請求項８】
請求項１記載の脈拍計において、
前記体動成分除去部は、前記体動成分として含まれる所定の低周波数帯域成分を所定の疑似低周波信号を用いて前記脈波検出信号から除去することを特徴とする脈拍計。
【請求項９】
請求項１記載の脈拍計において、
前記体動成分除去部は、前記体動成分として含まれる所定の低周波数帯域成分を所定の疑似低周波信号を用いて前記脈波検出信号から除去すべく、前記疑似低周波成分に基づいて適応フィルタ係数を生成するフィルタ係数生成部と、
今回の前記脈波検出信号から前回の前記脈波検出信号に前記適応フィルタ係数を適用して得られる体動成分を除去する除去処理部と、
を備えたことを特徴とする脈拍計。
【請求項１０】
請求項１記載の脈拍計において、
前記脈波検出信号に含まれる体動成分からピッチあるいは歩数を検出する体動情報検出部を備えたことを特徴とする脈拍計。
【請求項１１】
腕に装着される腕時計型情報機器において、
脈波センサを有し、脈波検出信号を出力する脈波検出部と、
装着者の腕の末梢方向に沿った軸をＸ軸とし、前記Ｘ軸に垂直で手の甲に垂直な軸をＺ軸とし、前記Ｘ軸および前記Ｚ軸に垂直な軸をＹ軸とし、前記Ｘ軸、前記Ｙ軸、前記Ｚ軸のそれぞれに沿った方向の加速度を検出する３軸加速度センサを有し、前記Ｘ軸に沿った方向の加速度成分であるＸ軸加速度成分、前記Ｙ軸に沿った方向の加速度成分であるＹ軸加速度成分および前記Ｚ軸に沿った方向の加速度成分であるＺ軸加速度成分に基づいて前記体動成分を生成する体動成分生成部と、
前記脈波検出信号から生成した前記体動成分を除去する体動成分除去部と、
前記体動成分除去後の前記脈波検出信号に基づいて脈拍数を算出する脈拍数算出部と、
前記脈拍数を表示する表示部と、
を備えたことを特徴とする腕時計型情報機器。
【請求項１２】
脈波センサを有し、脈波検出データを出力する脈波検出部と、装着者の腕の末梢方向に沿った軸をＸ軸とし、前記Ｘ軸に垂直で手の甲に垂直な軸をＺ軸とし、前記Ｘ軸および前
記Ｚ軸に垂直な軸をＹ軸とし、前記Ｘ軸、前記Ｙ軸、前記Ｚ軸のそれぞれに沿った方向の加速度を検出する３軸加速度センサを有し、静脈血の挙動に影響を与える体動に対応する加速度を検出して体動検出データとして出力する体動検出部と、を備え、人体に装着して脈拍を計測する脈拍計をコンピュータにより制御する制御プログラムにおいて、
前記体動検出データに対応する、前記Ｘ軸に沿った方向の加速度成分であるＸ軸加速度成分データ、前記Ｙ軸に沿った方向の加速度成分であるＹ軸加速度成分データおよび前記Ｚ軸に沿った方向の加速度成分であるＺ軸加速度成分データに基づいて体動成分データを生成させ、
前記体動成分データに基づいて前記脈波検出データに対応する脈波成分に含まれる体動成分を除去させ、
前記体動成分除去後の前記脈波成分に基づいて脈拍数を算出させる、
ことを特徴とする制御プログラム。
【請求項１３】
請求項１２記載の制御プログラムを記録したことを特徴とするコンピュータ読取可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、脈拍計、脈拍計の制御方法、腕時計型情報機器、制御プログラムおよび記録媒体に係り、特に人体の腕に装着して歩行中あるいは走行中の脈拍を測定するのに好適な脈拍計、脈拍計の制御方法、腕時計型情報機器、制御プログラムおよび記録媒体に関する。
【背景技術】
【０００２】
従来より、人体の一部に装着して歩行中あるいは走行中の脈拍を測定する脈波計が知られている。
例えば、特許文献１には、腕時計型の脈拍計が開示されている。
上記特許文献１に開示されている脈拍計は、加速度センサにより検出した体動信号の周波数分析結果に基づいて体動信号の高調波成分のすべてに相当する周波数成分を脈波信号の周波数分析結果から除去し、体動信号の高調波成分が除去された脈波信号の周波数分析結果のなかから最大のパワーを有する周波数成分を抽出し、当該抽出した周波数成分に基づいて脈拍数を算出するという構成を採っていた。
【特許文献１】特許第２８１６９４４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
上記従来の脈拍計においては、体動成分の検出を加速度センサにより行っていたため、脈拍センサ信号に含まれる生体内部で発生する体動成分を全て把握できているわけではなく、体動成分の除去が不完全となるという恐れがあった。
従来は体動成分を完全に把握できないため、脈拍センサ信号に含まれる体動成分を除去するために周波数分析結果の高調波成分の特徴を利用して体動信号を特定し、特定した体動信号を除去して脈波信号を抽出しているため、体動が周期的な特性を有していない場合には、体動成分を除去することができず、ひいては、脈拍を正しく求める事ができないという問題点があった。
そこで、本発明の目的は、脈拍センサ信号に含まれる体動成分をより正確に把握することにより、体動成分が周期的な特性を有していない場合であっても、脈波成分から生体内部で発生する体動成分を確実に除去して正確に脈拍数算出を行うことができる脈拍計、腕時計型情報機器、制御プログラムおよび記録媒体を提供することにある。
【課題を解決するための手段】
【０００４】
上記課題を解決するため、人体に装着して脈拍を計測する脈拍計は、脈波センサを有し、脈波検出信号を出力する脈波検出部と、加速度センサを有し、静脈血の挙動に影響を与える体動に対応する加速度を検出して体動検出信号として出力する体動検出部と、前記体動検出信号に基づいて前記脈波検出信号に含まれる体動成分を除去する体動成分除去部と、前記体動成分除去後の前記脈波検出信号に基づいて脈拍数を算出する脈拍数算出部と、を備え、前記加速度センサは、装着者の腕の末梢方向に沿った軸をＸ軸とし、前記Ｘ軸に垂直で手の甲に垂直な軸をＺ軸とし、前記Ｘ軸および前記Ｚ軸に垂直な軸をＹ軸とし、前記Ｘ軸、前記Ｙ軸、前記Ｚ軸のそれぞれに沿った方向の加速度を検出する３軸加速度センサであり、前記体動成分除去部は、前記体動検出信号に対応する、前記Ｘ軸に沿った方向の加速度成分であるＸ軸加速度成分、前記Ｙ軸に沿った方向の加速度成分であるＹ軸加速度成分および前記Ｚ軸に沿った方向の加速度成分であるＺ軸加速度成分に基づいて前記体動成分を生成する体動成分生成部を備えたことを特徴としている。
上記構成によれば、脈波検出部は、脈波検出信号を体動成分検出除去部に出力する。
一方、体動成分検出除去部の体動成分生成部は、体動検出信号に対応する、Ｘ軸に沿った方向の加速度成分であるＸ軸加速度成分、Ｙ軸に沿った方向の加速度成分であるＹ軸加速度成分およびＺ軸に沿った方向の加速度成分であるＺ軸加速度成分に基づいて体動成分を生成する。
これにより、体動成分除去部は、脈波検出信号に含まれる体動成分を除去する
これらにより、生体情報計測部は、体動成分除去後の脈波検出信号に基づいて生体情報を計測する。
【０００５】
この場合において、前記加速度センサは、前記脈波センサの近傍に配置されているようにしてもよい。
さらに、前記加速度センサは、前記脈波センサに略積層状態で配置されているようにしてもよい。
さらにまた、前記体動成分除去部は、前記Ｘ軸に沿った方向の加速度成分であるＸ軸加速度成分、前記Ｙ軸に沿った方向の加速度成分であるＹ軸加速度成分および前記Ｚ軸に沿った方向の加速度成分であるＺ軸加速度成分に基づいて前記体動成分を生成する体動成分生成部を備えるようにしてもよい。
【０００６】
また、前記体動成分除去部は、前記Ｙ軸加速度成分および前記Ｚ軸加速度成分をベクトルとみなし、これらのベクトルの合成ベクトルである２軸加速度合成成分と、前記Ｘ軸加速度成分と、に基づいて前記体動成分を生成する体動成分生成部を備えるようにしてもよ。
さらに、前記体動成分除去部は、前記Ｘ軸加速度成分、前記Ｙ軸加速度成分および前記Ｚ軸加速度成分の３つの加速度成分をベクトルとみなし、これらのベクトルの合成ベクトルである３軸加速度合成成分に基づいて前記体動成分を生成する体動成分生成部を備えるようにしてもよい。
さらにまた、前記Ｘ軸加速度成分、前記Ｙ軸加速度成分および前記Ｚ軸加速度成分のうち少なくともいずれか一の加速度成分について重み付けを行うようにしてもよい。
また、前記体動成分除去部は、前記Ｘ軸加速度成分、前記Ｙ軸加速度成分および前記Ｚ軸加速度成分に基づいて適応フィルタ係数を生成するフィルタ係数生成部と、今回の前記脈波検出信号から前回の前記脈波検出信号に前記適応フィルタ係数を適用して得られる体動成分を除去する除去処理部と、を備えるようにしてもよい。
さらに、前記体動成分除去部は、前記体動成分として含まれる所定の低周波数帯域成分を所定の疑似低周波信号を用いて前記脈波検出信号から除去するようにしてもよい。
【０００７】
さらにまた、前記体動成分除去部は、前記体動成分として含まれる所定の低周波数帯域成分を所定の疑似低周波信号を用いて前記脈波検出信号から除去すべく、前記疑似低周波成分に基づいて適応フィルタ係数を生成するフィルタ係数生成部と、今回の前記脈波検出信号から前回の前記脈波検出信号に前記適応フィルタ係数を適用して得られる体動成分を除去する除去処理部と、を備えるようにしてもよい。
また、前記脈波検出信号に含まれる体動成分からピッチあるいは歩数を検出する体動情報検出部を備えるようにしてもよい。
【００１０】
また、腕に装着される腕時計型情報機器は、脈波センサを有し、脈波検出信号を出力する脈波検出部と、前記加速度センサは、装着者の腕の末梢方向に沿った軸をＸ軸とし、前記Ｘ軸に垂直で手の甲に垂直な軸をＺ軸とし、前記Ｘ軸および前記Ｚ軸に垂直な軸をＹ軸とし、前記Ｘ軸、前記Ｙ軸、前記Ｚ軸のそれぞれに沿った方向の加速度を検出する３軸加速度センサを有し、前記Ｘ軸に沿った方向の加速度成分であるＸ軸加速度成分、前記Ｙ軸に沿った方向の加速度成分であるＹ軸加速度成分および前記Ｚ軸に沿った方向の加速度成分であるＺ軸加速度成分に基づいて前記体動成分を生成する体動成分生成部と、前記脈波検出信号から生成した前記体動成分を除去する体動成分除去部と、前記体動成分除去後の前記脈波検出信号に基づいて脈拍数を算出する脈拍数算出部と、前記脈拍数を表示する表示部と、を備えたことを特徴としている。
上記構成によれば、脈波検出部は、体動成分除去部に脈波検出信号を出力する。
体動成分生成部は、Ｘ軸に沿った方向の加速度成分であるＸ軸加速度成分、Ｙ軸に沿った方向の加速度成分であるＹ軸加速度成分およびＺ軸に沿った方向の加速度成分であるＺ軸加速度成分に基づいて体動成分を生成する。
これにより、体動成分除去部は、脈波検出信号から体動成分を除去し、脈拍数算出部は、体動成分除去後の脈波検出信号に基づいて脈拍数を算出し、表示部は、算出された脈拍数を表示する。
【００１１】
また、脈波センサを有し、脈波検出データを出力する脈波検出部と、装着者の腕の末梢方向に沿った軸をＸ軸とし、前記Ｘ軸に垂直で手の甲に垂直な軸をＺ軸とし、前記Ｘ軸および前記Ｚ軸に垂直な軸をＹ軸とし、前記Ｘ軸、前記Ｙ軸、前記Ｚ軸のそれぞれに沿った方向の加速度を検出する３軸加速度センサを有し、静脈血の挙動に影響を与える体動に対応する加速度を検出して体動検出データとして出力する体動検出部と、を備え、人体に装着して脈拍を計測する脈拍計をコンピュータにより制御する制御プログラムは、前記体動検出データに対応する、前記Ｘ軸に沿った方向の加速度成分であるＸ軸加速度成分データ、前記Ｙ軸に沿った方向の加速度成分であるＹ軸加速度成分データおよび前記Ｚ軸に沿った方向の加速度成分であるＺ軸加速度成分データに基づいて体動成分データを生成させ、前記体動成分データに基づいて前記脈波検出データに対応する脈波成分に含まれる体動成分を除去させ、前記体動成分除去後の前記脈波成分に基づいて脈拍数を算出させる、
ことを特徴としている。
さらに上記制御プログラムをコンピュータ読取可能な記録媒体に記録するようにしてもよい。
【発明の効果】
【００１２】
本発明によれば、脈波検出信号に含まれる体動成分を確実に検出し、体動成分除去後の脈波検出信号に基づいて脈拍数を算出するので、生体内部で発生する体動成分、特に静脈血に起因する体動成分を確実に除去し、脈波成分に対応するスペクトルを容易に特定して、脈拍検出精度を向上させることができる。
【発明を実施するための最良の形態】
【００１３】
次に本発明の好適な実施の形態について図面を参照して説明する。
図１は、実施形態の脈拍測定装置の装着状態の説明図である。
脈拍測定装置１０は、ユーザの腕１１に装着されて用いられ、装置本体（時計ケース）１０Ａと、装置本体１０Ａを腕に装着するためのリストバンド１０Ｂと、を備えている。
【００１４】
図２は、実施形態の脈拍測定装置の断面図である。
リストバンド１０Ｂを手首に巻きつけて脈拍測定装置１０を装着すると、装置本体１０Ａの裏面側が手首の甲に密着する。
そこで、装置本体１０Ａの裏面側には３軸（Ｘ軸、Ｙ軸、Ｚ軸）加速度センサ１２および脈波センサ１３が設けられている。この場合において、３軸加速度線センサ１２は、体動センサとして機能している。
図２に示すように、脈波センサユニット１３は、脈波の検出用光を射出するＬＥＤ１３Ａと、人体により反射された検出用光を受光するＰＤ（Photo Detector）１３Ｂと、ＬＥＤ１３ＡおよびＰＤ１３Ｂを保護するとともに、ＬＥＤ１３Ａの照射光、生体を介して得られる反射光を透過してＰＤ１３Ｂに入射させるための透明ガラス１３Ｃと、を備えている。ここで、透明ガラス１３Ｃは、装置本体１０Ａを構成する裏蓋１４により固定されている。
この脈波センサユニット１３の構成によれば、ＬＥＤ１３Ａからの光が透明ガラス５５を介して手首の甲に照射され、その反射光がフォトダイオード１３Ｂに受光されるようになっている。
【００１５】
装置本体１０Ａの表面側には、現在時刻や日付に加えて、脈波センサ１３の検出結果に基づく脈拍数ＨＲなどの生体情報も表示する液晶表示装置１５が設けられている。
また、装置本体１０Ａの内部には、メイン基板１６の上側にＣＰＵ等の各種ＩＣ回路が設けられており、これによってデータ処理回路１７が構成される。
また、メイン基板１６の裏面側には電池１８が設けられており、電池１８から３軸加速度センサ、脈波センサ１３、液晶表示装置１５およびメイン基板１６電源が供給されるようになっている。
３軸加速度センサ１２および脈波センサ１３とメイン基板１６とは、ヒートシール１９によって接続されている。これにより、ヒートシール１９形成される配線によって、メイン基板１６から電源が３軸加速度センサ１２および脈波センサ１３に供給される。
この結果、３軸加速度センサ１２からメイン基板１６に加速度検出信号が供給される用になっている。また、脈波センサ１３からメイン基板１６に脈波検出信号が供給されるようになっている。
【００１６】
データ処理回路１７は、加速度検出信号および脈波検出信号にＦＦＴ処理を施し、それらの処理結果を解析することにより、脈拍数ＨＲを算出している。なお、装置本体１０Ａの外側面には、図１に示すように、時刻合わせや表示モードの切換などを行うためのボタンスイッチ２０Ａ、２０Ｂ、２０Ｃ、２０Ｄ、２０Ｅが設けられている。
図３は実施形態の脈拍測定装置の概要構成ブロック図である。
脈拍測定装置１０は、大別すると、上述した３軸加速度センサ１２、脈波センサ１３および液晶表示装置１５の他に、脈波信号増幅回路２１と、加速度信号増幅回路２２と、Ａ／Ｄ変換回路２３と、ＭＰＵ２４と、ＲＡＭ２５と、ＲＯＭ２６と、を備えている。
３軸加速度センサ１２は、図１あるいは図２に示すＸ軸方向の加速度を検出するＸ軸加速度センサ１２Ｘと、Ｙ軸方向の加速度を検出するＹ軸加速度センサ１２Ｙと、Ｚ軸方向の加速度を検出するＺ軸加速度センサ１２Ｚと、を備えている。
脈波信号増幅回路２１は、脈波センサ１３から出力される脈波検出信号を所定増幅率で増幅して増幅脈波検出信号としてＡ／Ｄ変換回路２３に出力する。
加速度信号増幅回路２２は、３軸加速度センサ１２から出力されるＸ軸加速度検出信号、Ｙ軸加速度検出信号およびＺ軸加速度検出信号をそれぞれ所定増幅率で増幅して増幅Ｘ軸加速度検出信号、増幅Ｙ軸加速度検出信号および増幅Ｚ軸加速度検出信号としてＡ／Ｄ変換回路２３に出力する。
【００１７】
Ａ／Ｄ変換回路２３は、入力された増幅脈波検出信号、増幅Ｘ軸加速度検出信号、増幅Ｙ軸加速度検出信号および増幅Ｚ軸加速度検出信号および増幅圧力検出信号をそれぞれ別個にアナログ／ディジタル変換して脈波検出データ、Ｘ軸加速度検出データＫｘ、Ｙ軸加速度検出データＫｙおよびＺ軸加速度検出データＫｚとしてＭＰＵ２４へ出力する。
ＭＰＵ２４は、脈波検出データ、Ｘ軸加速度検出データＫｘ、Ｙ軸加速度検出データＫｙおよびＺ軸加速度検出データＫｚをＲＡＭ２５に格納するとともに、ＲＯＭ２６に格納された制御プログラムに基づいて脈拍数を算出し、表示装置１５に表示することとなる。
より詳細には、ＭＰＵ２４は、ＲＡＭ２５に格納した脈波検出データ並びにＸ軸加速度検出データＫｘ、Ｙ軸加速度検出データＫｙおよびＺ軸加速度検出データＫｚに基づいて得られる体動検出データを時系列順に並べ、対応するサンプリングタイミング毎に脈波検出データと体動検出検出データとの両者の差である残差データを求める。
そして残差データの周波数分析（ＦＦＴ：Fast fourier Transformation）を行って、脈波の高調波成分を抽出し、その周波数から脈拍数を算出することとなる。
【００１８】
次により具体的な脈拍数算出処理について説明する。
まず、実施形態の具体的説明に先立ち、実施形態の動作原理を説明する。
脈波を検出するための脈波センサの出力には、脈波成分に加えて様々な体動成分が含まれている。この体動成分は、脈拍被測定者であるユーザの運動（歩行・走行動作、腕の振りなど）に起因する、生体内部の変化により発生するものであることが知られている。
ところで、体動成分を検出するセンサとして３軸加速度センサを用いた場合、特に末梢方向、すなわち、Ｘ軸方向の体動成分の影響が大きいことがわかっているが、他の２軸（Ｙ軸およびＺ軸）方向の体動成分についても無視することができない。
そこで、発明者らは、同一の体動成分を発生させた場合の２軸方向の加速度をベクトルとして、２軸の加速度ベクトルの合成ベクトルの変化量と脈波センサの出力に含まれる体動成分量（ストローク成分量）との関係を調べた。
【００１９】
ここで、３軸加速度センサ１２の構成について詳細に検討する。
図４は、３軸（Ｘ、Ｙ、Ｚ軸）加速度センサ１２を用いたときの後述するＸ軸方向の加速度と脈波センサの出力信号に含まれる体動成分（ストローク成分）との関係説明図である。
図５は、３軸加速度センサ１２を用いたときの後述するＹ軸方向の加速度と脈波センサの出力信号に含まれる体動成分（ストローク成分）との関係説明図である。
【００２０】
図６は、３軸加速度センサ１２を用いたときの後述するＺ軸方向の加速度と脈波センサの出力信号に含まれる体動成分（ストローク成分）との関係説明図である。
図７は３軸の説明図である。
図７に示すように、Ｘ軸は、図１に示した末梢方向に延在する軸であり、Ｙ軸は、手のひらを平面上に載置した場合にＸ軸に垂直で平面上に存在する軸であり、Ｚ軸は、手のひらを載置した平面に垂直な軸である。
図４ないし図６に示したように、脈波センサの出力信号に含まれる体動成分は、Ｘ軸方向の成分が支配的であることがわかる。従って、加速度センサとしては、Ｘ軸方向、すなわち、図１に示した末梢方向の加速度のみを検出可能な１軸加速度センサを用いれば、脈波センサで検出される体動成分を推定することが可能となる。
これは、換言すれば、２軸の加速度ベクトルの合成ベクトル変化量を検出できれば、脈波センサの出力に含まれている静脈血の影響量を推し量ることができるということである。
そこで、本実施形態においては、静脈起因の体動成分を３軸加速度センサ１２により検出するとともに、これらの検出出力を脈波センサの出力から所定割合で差し引くことにより、静脈血の影響を除去した信号に基づいて正確に脈拍数を検出するようにしている。
図８は、実施形態の適応フィルタの一例の概要構成ブロック図である。
適応フィルタ３０は、大別するとフィルタ係数生成部３１と、合成部３２と、を備えている。
フィルタ係数生成部３１の係数コントロール部３１Ａは、体動成分除去部として機能しており、合成部３２が前回出力したフィルタ適用後のデータに基づいて適応フィルタ係数ｈを生成する。
【００２１】
そしてフィルタ係数生成部３１は、入力された体動成分検出信号であるＸ軸加速度データＫｘ、Ｙ軸加速度データＫｙとＺ軸加速度データＫｚの合成データである合成加速度ベクトルデータ（＝ｙ）および疑似低周波信号（＝ｚ）に係数コントロール部３１Ａが生成した適用フィルタ係数ｈを適用して、体動除去データｈ（ｘ）、ｈ（ｙ）、ｈ（ｚ）をそれぞれ生成し、合成部３２に出力する。
合成部３２は、除去処理部として機能しており、前回抽出した脈波検出データ（＝脈波成分＋体動成分）と体動除去データｈ（ｘ）、ｈ（ｙ）、ｈ（ｚ）を合成し、今回の脈波検出データに含まれる体動成分を実質的に除去（減算）して、脈波成分ｅ（ｎ）を抽出する。
【００２２】
ここで、疑似低周波信号を用いる理由について説明する。
発明者らの実験によれば、脈波検出データから体動除去データｈ（ｘ）、ｈ（ｙ）を除去したとしても、得られる脈波成分には低周波変動成分が残り正確に脈拍数が求まらない場合があった。
これは、呼吸や神経の影響と考えられるが、これを検出して影響を除去するためには大がかりなシステムが必要となり、携帯型の脈拍測定装置を実現することは不可能になってしまう。
そこで、発明者らは、体動検出センサである３軸加速度センサの出力信号に当該低周波変動成分に相当する周波数の疑似低周波信号を掛け合わせた上で適応フィルタを施すことにより、その影響を除去することができた。
この場合において、疑似低周波信号は、周波数分析の際に所定の周波数分布を持って低周波変動成分を除去する必要があることとと、その周波数帯域が０．５Ｈｚ以下であるという観点から、０．５Ｈｚ以下の三角波もしくは矩形波であることが望ましい。なお、その周波数帯域および波形形状は実際に含まれる低周波変動成分に応じて適宜変更が可能である。
【００２３】
次に本実施形態の具体的な脈拍数算出処理について説明する。
図９は、Ｘ軸加速度センサ１２Ｘから出力されるＸ軸加速度検出信号に対応するＸ軸加速度データＫｘを時系列順に並べたグラフである。
また、図１０は、図９に対応するＸ軸加速度検出データＫｘにＦＦＴを施して得られた周波数分析結果である。
図１１は、Ｙ軸加速度センサ１２Ｙから出力されるＹ軸加速度検出信号に対応するＹ軸加速度データＫｙを時系列順に並べたグラフである。
また、図１２は、図１１に対応するＹ軸加速度検出データＫｙにＦＦＴを施して得られた周波数分析結果である。
図１３は、Ｚ軸加速度センサ１２Ｚから出力されるＺ軸加速度検出信号に対応するＺ軸加速度データＫｚを時系列順に並べたグラフである。
また、図１４は、図１３に対応するＺ軸加速度検出データＫｚにＦＦＴを施して得られた周波数分析結果である。
図９と、図１１および図１３とを比較すると、体動成分としては、Ｘ軸加速度成分の影響がＹ軸加速度成分およびＺ軸加速度成分より影響が大きいことがわかる。
そこで、発明者らは、測定精度を維持しつつ、処理の簡略化を図るべく上述したようにＹ軸加速度成分およびＺ軸加速度成分を一体として取り扱うこととし、２軸の加速度ベクトルの合成ベクトル変化量を検出することとした。
図１５は、Ｙ軸加速度センサ１２Ｙから出力されるＹ軸加速度検出信号に対応するＹ軸加速度データＫｙと、Ｚ軸加速度センサ１２Ｚから出力されるＺ軸加速度検出信号に対応するＺ軸加速度データＫｚと、をそれぞれベクトルとして扱い、その合成ベクトルとして得られる合成加速度ベクトルデータを時系列順に並べたグラフである。
また、図１６は、図１５に対応する合成加速度ベクトルデータ（＝√（Ｋｙ2 ＋Ｋｚ2 ））すなわち、２軸加速度合成成分にＦＦＴを施して得られた周波数分析結果である。
図１７は、あらかじめ設定した疑似低周波信号（三角波を使用）を時系列順に並べたグラフである。
また、図１８は、図１７に対応する疑似低周波信号にＦＦＴを施して得られた周波数分析結果である。
【００２４】
図１８をみてわかるとおり、その周波数はおよそ０．５Ｈｚ以下であり、所定の周波数分布を持っている。
図１９は、脈波検出データの一例を時系列順に並べてグラフ化したものである。
また、図２０は、図１９に対応する脈波検出データにＦＦＴを施して得られた周波数分析結果である。
まず、ＭＰＵ２４は、ＲＡＭ２５に格納した脈波検出データ、Ｘ軸加速度検出データ、Ｙ軸加速度検出データおよびＺ軸加速度検出データを順次読み出し、一つのサンプリングタイミングにおける脈波検出データを合成部３２に出力する。
これと並行して、ＭＰＵ２４は、合成部３２に出力した脈波検出データに対応するＸ軸加速度検出データＫｘ、Ｙ軸加速度検出データＫｙおよびＺ軸加速度検出データＫｚをフィルタ係数生成部３１に出力する。
【００２５】
これによりフィルタ係数生成部３１の係数コントロール部３１Ａは、合成部３２が前回出力したフィルタ適用後のデータに基づいて適応フィルタ係数ｈを生成する。
そしてフィルタ係数生成部３１は、入力された体動成分検出信号であるＸ軸加速度データＫｘ（＝ｘ）、Ｙ軸加速度データＫｙとＺ軸加速度データＫｚの合成データである合成加速度ベクトルデータ（＝ｙ）および疑似低周波信号（＝ｚ）に係数コントロール部３１Ａの制御下でそれぞれ適用フィルタ係数ｈを適用して、体動除去データｈ（ｘ）、ｈ（ｙ）、ｈ（ｚ）をそれぞれ生成し、合成部３２に出力する。
これにより合成部３２は、今回の脈波データと体動除去データｈ（ｘ）、ｈ（ｙ）、ｈ（ｚ）とを合成して、今回の脈波検出データに含まれる体動成分を実質的に除去（減算）して、脈波成分を抽出して適応フィルタ適用後のデータである残差データｅ（ｎ）を出力する。
【００２６】
図２１は、図１７の脈波検出信号に対して図９の増幅Ｘ軸加速度検出信号、図１５の合成加速度ベクトル信号および図１７の疑似低周波信号に適応フィルタを適用して得られた信号を合成して得られた残差データを時系列順に並べてグラフ化したものである。
次にＭＰＵ２４は、残差データに対しＦＦＴを施す。
図２２は、図２１の残差データにＦＦＴを施して得られた周波数分析結果である。
これにより、得られた周波数分析結果は、実質的に脈波センサの出力信号（脈波成分＋体動成分）から、静脈起因の体動成分を除去したもの、すなわち、主として脈波成分に対応する脈波データとなる。
ここで、比較のため、疑似低周波信号を用いなかった場合に得られる脈波データについて説明する。
図２３は、図１７の脈波検出信号に対して図９の増幅Ｘ軸加速度検出信号および図１５の合成加速度ベクトル信号に適応フィルタを適用して得られた信号を合成して得られた残差データを時系列順に並べてグラフ化したものである。
また、図２４は、図２３の残差データにＦＦＴを施して得られた周波数分析結果である。
【００２７】
図２２と図２４とを比較することにより、本実施形態の構成によれば、低周波変動成分の低減が図れ、ひいては、脈拍数検出における低周波変動成分の影響を容易に除去できることが容易にわかる。
さらにＭＰＵ２４は、主として脈波成分を含む得られた脈波データから最大の周波数成分を脈拍スペクトルとして、その周波数から脈拍数を算出する。
そしてＭＰＵ２４は、液晶表示装置１５に脈拍数を表示することとなる。
以上の説明のように、本実施形態によれば、体動センサとして機能する３軸加速度センサ１２および脈波センサ１３を用いて、さらに疑似低周波信号を用いることにより、生体内部で発生する体動成分の主要因である静脈の変動を確実に検出して把握できる。このため、体動成分を確実に除去でき、正確な脈波成分検出、ひいては、正確な脈拍数の測定が行える。
【００２８】
［１．１］第１変形例
以上の説明においては、Ｙ軸加速度データＫｙとＺ軸加速度データＫｚの合成データである合成加速度ベクトルデータ（＝√（Ｋｙ2 ＋Ｋｚ2 ））を用いる場合の実施形態であったが、本第１変形例は、Ｘ軸加速度データ、Ｙ軸加速度データおよびＺ軸加速度データの３つの加速度データを合成した合成加速度ベクトルデータ（＝√（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2 ））、すなわち、３軸加速度合成成分を用いる場合の実施形態である。
図２５は、実施形態の第１変形例の適応フィルタの一例の概要構成ブロック図である。
適応フィルタ４０は、大別すると、フィルタ係数生成部４１と、積算部４２と、合成部４３と、を備えている。
フィルタ係数生成部４１は、体動成分除去部として機能しており、合成部４３が前回出力したフィルタ適用後のデータに基づいて適応フィルタ係数ｈを生成する。
これと並行して、積算部４２は、Ｘ軸加速度データ、Ｙ軸加速度データおよびＺ軸加速度データの３つの加速度データを合成した合成加速度ベクトルデータ（＝√（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2 ））に対し、あらかじめ設定した疑似低周波信号を掛け合わせてフィルタ係数生成部４１に出力する。
この結果、フィルタ係数生成部４１は、積算部４２の出力に生成した適用フィルタ係数ｈを適用して、体動除去データｈ（√（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2））を生成し、合成部４３に出力する。
合成部４３は、除去処理部として機能しており、前回抽出した脈波検出データ（＝脈波成分＋体動成分）と体動除去データｈ（√（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2））を合成し、今回の脈波検出データに含まれる体動成分を実質的に除去（減算）して、脈波成分ｅ（ｎ）を抽出する。
【００２９】
次に本第１変形例の具体的な脈拍数算出処理について説明する。
図２６は、Ｘ軸加速度センサ１２Ｘから出力されるＸ軸加速度検出信号に対応するＸ軸加速度データＫｘを時系列順に並べたグラフである。
また、図２７は、図２６に対応するＸ軸加速度検出データＫｘにＦＦＴを施して得られた周波数分析結果である。
図２８は、Ｙ軸加速度センサ１２Ｙから出力されるＹ軸加速度検出信号に対応するＹ軸加速度データＫｙを時系列順に並べたグラフである。
また、図２９は、図２８に対応するＹ軸加速度検出データＫｙにＦＦＴを施して得られた周波数分析結果である。
【００３０】
図３０は、Ｚ軸加速度センサ１２Ｚから出力されるＺ軸加速度検出信号に対応するＺ軸加速度データＫｚを時系列順に並べたグラフである。
また、図３１は、図３０に対応するＺ軸加速度検出データＫｚにＦＦＴを施して得られた周波数分析結果である。
図３２は、Ｘ軸加速度センサ１２Ｘから出力されるＸ軸加速度データＫｚと、Ｙ軸加速度センサ１２Ｙから出力されるＹ軸加速度検出信号に対応するＹ軸加速度データＫｙと、Ｚ軸加速度センサ１２Ｚから出力されるＺ軸加速度検出信号に対応するＺ軸加速度データＫｚと、をそれぞれベクトルとして扱い、その合成ベクトルとして得られる合成加速度ベクトルデータ（＝√（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2 ））を時系列順に並べたグラフである。
また、図３３は、図３２に対応する合成加速度ベクトルデータ（＝√（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2 ））にＦＦＴを施して得られた周波数分析結果である。
図３４は、脈波検出データの一例を時系列順に並べてグラフ化したものである。
また、図３５は、図３４に対応する脈波検出データにＦＦＴを施して得られた周波数分析結果である。
【００３１】
まず、ＭＰＵ２４は、ＲＡＭ２５に格納した脈波検出データ、Ｘ軸加速度検出データ、Ｙ軸加速度検出データおよびＺ軸加速度検出データを順次読み出し、一つのサンプリングタイミングにおける脈波検出データを合成部４３に出力する。
これと並行して、ＭＰＵ２４は、合成部４３に出力した脈波検出データに対応するＸ軸加速度検出データＫｘ、Ｙ軸加速度検出データＫｙおよびＺ軸加速度検出データＫｚを積算部４２に出力する。
積算部４２は、Ｘ軸加速度データ、Ｙ軸加速度データおよびＺ軸加速度データの３つの加速度データを合成した合成加速度ベクトルデータ（＝√（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2 ））に対し、図１７および図１８に示したような疑似低周波信号を掛け合わせてフィルタ係数生成部４１に出力する。
これによりフィルタ係数生成部４１は、合成部４３が前回出力したフィルタ適用後のデータに基づいて適応フィルタ係数ｈを生成する。
そしてフィルタ係数生成部４１は、入力された合成加速度ベクトルデータ（＝√（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2 ））に適用フィルタ係数ｈを適用して、体動除去データｈ（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2 ）を生成し、合成部４３に出力する。
これにより合成部４３は、今回の脈波データと体動除去データｈ（√（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2 ）とを合成して、今回の脈波検出データに含まれる体動成分を実質的に除去（減算）して、脈波成分を抽出して適応フィルタ適用後のデータである残差データｅ（ｎ）を出力する。
【００３２】
図３６は、図３４の脈波検出データに対して図３２の合成加速度ベクトルデータおよび図１７の疑似低周波信号に適応フィルタを適用して得られたデータを合成して得られた残差データを時系列順に並べてグラフ化したものである。
次にＭＰＵ２４は、残差データに対しＦＦＴを施す。
図３７は、図３６の残差データにＦＦＴを施して得られた周波数分析結果である。
これにより、得られた周波数分析結果は、実施形態の場合と比較して低周波数領域（＜０．５Ｈｚ）に脈波成分とは無関係なスペクトルが残っているものの、脈波成分の周波数帯域（２Ｈｚ〜２．５Ｈｚ）に影響を与えるものではなく、実質的に脈波センサの出力信号（脈波成分＋体動成分）から、静脈起因の体動成分を除去したもの、すなわち、主として脈波成分に対応する脈波データとなる。
【００３３】
［１．２］第２変形例
以上の説明においては、疑似低周波信号を処理に用いる場合であったが、本第２変形例は、処理および装置構成の簡略化のため、疑似低周波信号を処理に用いない場合の変形例である。
図３８は、実施形態の第２変形例の適応フィルタの一例の概要構成ブロック図である。
適応フィルタ５０は、大別するとフィルタ係数生成部５１と、合成部５２と、を備えている。
【００３４】
フィルタ係数生成部５１の係数コントロール部５１Ａは、体動成分除去部として機能しており、合成部５２が前回出力したフィルタ適用後のデータに基づいて適応フィルタ係数ｈを生成する。
そしてフィルタ係数生成部５１は、入力された体動成分検出信号であるＸ軸加速度データＫｘ、Ｙ軸加速度データＫｙおよびＺ軸加速度データｚに係数コントロール部５１Ａが生成した適用フィルタ係数ｈを適用して、体動除去データｈ（ｘ）、ｈ（ｙ）、ｈ（ｚ）をそれぞれ生成し、合成部５２に出力する。
合成部５２は、除去処理部として機能しており、前回抽出した脈波検出データ（＝脈波成分＋体動成分）と体動除去データｈ（ｘ）、ｈ（ｙ）、ｈ（ｚ）を合成し、今回の脈波検出データに含まれる体動成分を実質的に除去（減算）して、脈波成分ｅ（ｎ）を抽出する。
次に具体的な処理データの一例を説明する。
図３９は、脈波検出データの一例を時系列順に並べてグラフ化したものである。
また、図４０は、図３９に対応する脈波検出データにＦＦＴを施して得られた周波数分析結果である。
図４１は、図３４の脈波検出データに対して図３２の合成加速度ベクトル信号および図１７の疑似低周波信号に適応フィルタを適用して得られた信号を合成して得られた残差データを時系列順に並べてグラフ化したものである。
図４２は、図４１の残差データにＦＦＴを施して得られた周波数分析結果である。
【００３５】
ＭＰＵ２４が残差データｅ（ｎ）に対しＦＦＴを施すことにより、得られた周波数分析結果は、図３７に示すように、実施形態の場合と同様に実質的に脈波センサの出力信号（脈波成分＋体動成分）から、静脈起因の体動成分を除去したもの、すなわち、主として脈波成分に対応する脈波データとなる。
［１．３］第３変形例
本第３変形例は、上述した第１変形例のさらなる変形例であり、第１変形例において、疑似低周波信号を用いない構成を採用する変形例である。
図４３は、実施形態の第３変形例の適応フィルタの一例の概要構成ブロック図である。
適応フィルタ６０は、大別すると、フィルタ係数生成部６１と、合成部６２と、を備えている。
フィルタ係数生成部６１は、体動成分除去部として機能しており、合成部６２が前回出力したフィルタ適用後のデータに基づいて適応フィルタ係数ｈを生成する。
さらにフィルタ係数生成部６１は、Ｘ軸加速度データＫｘ、Ｙ軸加速度データＫｙおよびＺ軸加速度データＫｚの３つの加速度データを合成した合成加速度ベクトルデータ（＝√（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2 ））に生成した適用フィルタ係数ｈを適用して、体動除去データｈ（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2 ）を生成し、合成部６２に出力する。
合成部６２は、除去処理部としあて機能しており、前回抽出した脈波検出データ（＝脈波成分＋体動成分）と体動除去データｈ（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2 ）を合成し、今回の脈波検出データに含まれる体動成分を実質的に除去（減算）して、脈波成分ｅ（ｎ）を抽出する。
本第３変形例によれば、第１変形例と同様の効果を得ることができるとともに、疑似低周波信号を用いていないため、装置構成および処理の簡略化をよりいっそう図ることができる。
【００３６】
［１．４］第４変形例
本第４変形例は、実施形態において、疑似低周波信号を用いない構成を採用する変形例である。
図４４は、実施形態の適応フィルタの一例の概要構成ブロック図である。
適応フィルタ７０は、大別するとフィルタ係数生成部７１と、合成部７２と、を備えている。
フィルタ係数生成部７１の係数コントロール部７１Ａは、体動成分除去部として機能しており、合成部７２が前回出力したフィルタ適用後のデータに基づいて適応フィルタ係数ｈを生成する。
【００３７】
そしてフィルタ係数生成部７１は、入力された体動成分検出信号であるＸ軸加速度データＫｘ（＝ｘ）およびＹ軸加速度データＫｙとＺ軸加速度データＫｚの合成データである合成加速度ベクトルデータ（＝ｙ）に係数コントロール部３１Ａが生成した適用フィルタ係数ｈを適用して、体動除去データｈ（ｘ）、ｈ（ｙ）をそれぞれ生成し、合成部７２に出力する。
合成部７２は、除去処理部として機能しており、前回抽出した脈波検出データ（＝脈波成分＋体動成分）と体動除去データｈ（ｘ）、ｈ（ｙ）を合成し、今回の脈波検出データに含まれる体動成分を実質的に除去（減算）して、脈波成分ｅ（ｎ）を抽出する。
本第４変形例によれば、実施形態と同様の効果を得ることができるとともに、疑似低周波信号を用いていないため、装置構成および処理の簡略化をよりいっそう図ることができる。
【００３８】
［１．５］第５変形例
以上の説明においては、Ｘ軸加速度データＫｘ、Ｙ軸加速度データＫｙおよびＺ軸加速度データＫｚの３つの加速度データを合成した合成加速度ベクトルデータ（＝√（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2 ））あるいはＹ軸加速度データＫｙおよびＺ軸加速度データＫｚの２つの加速度データを合成した合成加速度ベクトルデータ（＝√（Ｋｙ2 ＋Ｋｚ2 ））を算出するに際して、何ら重み付けを行っていなかったが、各合成加速度ベクトルデータの元となる加速度データに適宜重み付けを行うように構成することも可能である。
例えば、Ｘ軸加速度データＫｘ、Ｙ軸加速度データＫｙおよびＺ軸加速度データＫｚの３つの加速度データにより合成加速度ベクトルデータを求めるに際し、次式を用いても良い。
√（ａ・Ｋｘ2 ＋ｂ・Ｋｙ2 ＋ｃＫｚ2）
ここで、ａ＞ｂ≧ｃ＞０とする。
また、合成加速度ベクトルデータを用いない場合であっても、Ｘ軸加速度データＫｘ、Ｙ軸加速度データＫｙおよびＺ軸加速度データＫｚに同様に適宜重み付けを行って適応フィルタ係数を適用するようにしてもよい。
さらに、疑似低周波信号についても重み付けを行うようにしてもよい。
さらにまた、以上の説明においては、３軸加速度センサを腕に取り付ける場合について説明したが、指元や指先に設けるようにすることも可能である。
［１．６］第６変形例
以上の説明においては、制御用プログラムがＲＯＭ２６内にあらかじめ記憶されている場合について説明したが、各種磁気ディスク、光ディスク、メモリカードなどの記録媒体に制御用プログラムをあらかじめ記録し、これらの記録媒体から読み込み、インストールするように構成することも可能である。また、通信インターフェースを設け、インターネット、ＬＡＮなどのネットワークを介して制御用プログラムをダウンロードし、インストールして実行するように構成することも可能である。
【図面の簡単な説明】
【００３９】
【図１】実施形態の脈拍測定装置の装着状態の説明図である。
【図２】実施形態の脈拍測定装置の断面図である。
【図３】実施形態の脈拍測定装置の概要構成ブロック図である。
【図４】３軸（Ｘ、Ｙ、Ｚ軸）加速度センサ１２を用いたときのＸ軸方向の加速度と脈波センサの出力信号に含まれる体動成分との関係説明図である。
【図５】３軸加速度センサ１２を用いたときの後述するＹ軸方向の加速度と脈波センサの出力信号に含まれる体動成分との関係説明図である。
【図６】３軸加速度センサ１２を用いたときの後述するＺ軸方向の加速度と脈波センサの出力信号に含まれる体動成分（ストローク成分）との関係説明図である。
【図７】３軸の説明図である。
【図８】実施形態の適応フィルタの一例の概要構成ブロック図である。
【図９】Ｘ軸加速度センサ１２Ｘから出力されるＸ軸加速度検出信号に対応するＸ軸加速度データＫｘを時系列順に並べたグラフである。
【図１０】図９に対応するＸ軸加速度検出データＫｘにＦＦＴを施して得られた周波数分析結果である。
【図１１】Ｙ軸加速度センサ１２Ｙから出力されるＹ軸加速度検出信号に対応するＹ軸加速度データＫｙを時系列順に並べたグラフである。
【図１２】図１１に対応するＹ軸加速度検出データＫｙにＦＦＴを施して得られた周波数分析結果である。
【図１３】Ｚ軸加速度センサ１２Ｚから出力されるＺ軸加速度検出信号に対応するＺ軸加速度データＫｚを時系列順に並べたグラフである。
【図１４】図１３に対応するＺ軸加速度検出データＫｚにＦＦＴを施して得られた周波数分析結果である。
【図１５】Ｙ軸加速度センサ１２Ｙから出力されるＹ軸加速度検出信号に対応するＹ軸加速度データＫｙと、Ｚ軸加速度センサ１２Ｚから出力されるＺ軸加速度検出信号に対応するＺ軸加速度データＫｚと、をそれぞれベクトルとして扱い、その合成ベクトルとして得られる合成加速度ベクトルデータを時系列順に並べたグラフである。
【図１６】図１５に対応する合成加速度ベクトルデータ（＝√（Ｋｙ2 ＋Ｋｚ2 ））にＦＦＴを施して得られた周波数分析結果である。
【図１７】あらかじめ設定した疑似低周波信号（三角波を使用）を時系列順に並べたグラフである。
【図１８】図１７に対応する疑似低周波信号にＦＦＴを施して得られた周波数分析結果である。
【図１９】脈波検出データの一例を時系列順に並べてグラフ化したものである。
【図２０】図１７に対応する脈波検出データにＦＦＴを施して得られた周波数分析結果である。
【図２１】図１７の脈波検出信号に対して図９の増幅Ｘ軸加速度検出信号、図１５の合成加速度ベクトル信号および図１７の疑似低周波信号に適応フィルタを適用して得られた信号を合成して得られた残差データを時系列順に並べてグラフ化したものである。
【図２２】図２１の残差データにＦＦＴを施して得られた周波数分析結果である。
【図２３】図１７の脈波検出信号に対して図９の増幅Ｘ軸加速度検出信号および図１５の合成加速度ベクトル信号に適応フィルタを適用して得られた信号を合成して得られた残差データを時系列順に並べてグラフ化したものである。
【図２４】図２３の残差データにＦＦＴを施して得られた周波数分析結果である。
【図２５】実施形態の第１変形例の適応フィルタの一例の概要構成ブロック図である。
【図２６】Ｘ軸加速度データＫｘを時系列順に並べたグラフである。
【図２７】図２６に対応するＸ軸加速度検出データＫｘにＦＦＴを施して得られた周波数分析結果である。
【図２８】Ｙ軸加速度データＫｙを時系列順に並べたグラフである。
【図２９】図２８に対応するＹ軸加速度検出データＫｙにＦＦＴを施して得られた周波数分析結果である。
【図３０】Ｚ軸加速度データＫｚを時系列順に並べたグラフである。
【図３１】図３０に対応するＺ軸加速度検出データＫｚにＦＦＴを施して得られた周波数分析結果である。
【図３２】合成加速度ベクトルデータ（＝√（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2 ））を時系列順に並べたグラフである。
【図３３】合成加速度ベクトルデータ（＝√（Ｋｘ2 ＋Ｋｙ2 ＋Ｋｚ2 ））にＦＦＴを施して得られた周波数分析結果である。
【図３４】脈波検出データの一例を時系列順に並べてグラフ化したものである。
【図３５】図３４に対応する脈波検出データにＦＦＴを施して得られた周波数分析結果である。
【図３６】図３６は、図３４の脈波検出データに対して図３２の合成加速度ベクトルデータおよび図１７の疑似低周波信号に適応フィルタを適用して得られたデータを合成して得られた残差データを時系列順に並べてグラフ化したものである。
【図３７】図３６の残差データにＦＦＴを施して得られた周波数分析結果である。
【図３８】実施形態の第２変形例の適応フィルタの一例の概要構成ブロック図である。
【図３９】脈波検出データの一例を時系列順に並べてグラフ化したものである。
【図４０】図３９に対応する脈波検出データにＦＦＴを施して得られた周波数分析結果である。
【図４１】図３４の脈波検出データに対して図３２の合成加速度ベクトル信号および図１７の疑似低周波信号に適応フィルタを適用して得られた信号を合成して得られた残差データを時系列順に並べてグラフ化したものである。
【図４２】図４１の残差データにＦＦＴを施して得られた周波数分析結果である。
【図４３】実施形態の第３変形例の適応フィルタの一例の概要構成ブロック図である。
【図４４】実施形態の適応フィルタの一例の概要構成ブロック図である。
【符号の説明】
【００４０】
１０…脈拍測定装置、１０Ａ…装置本体（時計ケース）、１０Ｂ…リストバンド１０Ｂ、１２…３軸（Ｘ軸、Ｙ軸、Ｚ軸）加速度センサ、１２Ｘ…Ｘ軸加速度センサ、１２Ｙ…ｙ軸加速度センサ、１２Ｚ…Ｚ軸加速度センサ、１３…脈波センサ、１３Ａ…ＬＥＤ、１３Ｂ…ＰＤ、１３Ｃ…透明ガラス、１５…液晶表示装置、２１…脈波信号増幅回路、２２…体動信号増幅回路、２３…Ａ／Ｄ変換回路、２４…ＭＰＵ、２５…ＲＡＭ、２６…ＲＯＭ、２７…表示装置、３０、４０、５０、６０、７０…適応フィルタ、３１、４１、５１、６１、７１…適応フィルタ係数生成部、３１Ａ、５１Ａ、７１Ａ…係数コントロール部、３２、４３、５２…合成部、４２…積算部。
【図１】
0003726832.tif 000002
【図２】
0003726832.tif 000003
【図３】
0003726832.tif 000004
【図４】
0003726832.tif 000005
【図５】
0003726832.tif 000006
【図６】
0003726832.tif 000007
【図７】
0003726832.tif 000008
【図８】
0003726832.tif 000009
【図９】
0003726832.tif 000010
【図１０】
0003726832.tif 000011
【図１１】
0003726832.tif 000012
【図１２】
0003726832.tif 000013
【図１３】
0003726832.tif 000014
【図１４】
0003726832.tif 000015
【図１５】
0003726832.tif 000016
【図１６】
0003726832.tif 000017
【図１７】
0003726832.tif 000018
【図１８】
0003726832.tif 000019
【図１９】
0003726832.tif 000020
【図２０】
0003726832.tif 000021
【図２１】
0003726832.tif 000022
【図２２】
0003726832.tif 000023
【図２３】
0003726832.tif 000024
【図２４】
0003726832.tif 000025
【図２５】
0003726832.tif 000026
【図２６】
0003726832.tif 000027
【図２７】
0003726832.tif 000028
【図２８】
0003726832.tif 000029
【図２９】
0003726832.tif 000030
【図３０】
0003726832.tif 000031
【図３１】
0003726832.tif 000032
【図３２】
0003726832.tif 000033
【図３３】
0003726832.tif 000034
【図３４】
0003726832.tif 000035
【図３５】
0003726832.tif 000036
【図３６】
0003726832.tif 000037
【図３７】
0003726832.tif 000038
【図３８】
0003726832.tif 000039
【図３９】
0003726832.tif 000040
【図４０】
0003726832.tif 000041
【図４１】
0003726832.tif 000042
【図４２】
0003726832.tif 000043
【図４３】
0003726832.tif 000044
【図４４】
0003726832.tif 000045