(57)【特許請求の範囲】
【請求項１】
生体の心拍波形を検出する心拍波形検出部と、
前記心拍波形から心拍を検出するための信号処理部とを備え、
前記信号処理部は、
前記心拍波形に対して信号処理を行って心拍信号を形成する心拍信号形成処理部と、
前記心拍信号形成処理部により形成された心拍信号から心拍を検出する心拍検出処理部とを備え、
前記心拍信号形成処理部は、前記心拍波形の主成分以外の周波数帯域にあるノイズ成分を除去するフィルタ部と、前記フィルタ部を通過した前記心拍波形の周波数成分を信号増幅する信号増幅部とを有し、当該信号増幅部の信号増幅特性は、心拍波形の低域側の周波数成分を低増幅率で増幅し、心拍波形の高域側の周波数成分を高増幅率で増幅することを特徴とする心拍計。
【請求項２】
生体の心拍波形を検出する心拍波形検出部と、
前記心拍波形から心拍を検出するための信号処理部とを備え、
前記信号処理部は、
前記心拍波形に対して信号処理を行って心拍信号を形成する心拍信号形成処理部と、
前記心拍信号形成処理部により形成された心拍信号から心拍を検出する心拍検出処理部とを備え、
前記心拍信号形成処理部は、前記心拍波形の主成分以外の周波数帯域にあるノイズ成分を除去するフィルタ部と、前記フィルタ部を通過した前記心拍波形の周波数成分を微分特性で信号増幅する信号増幅部とを有することを特徴とする心拍計。
【請求項３】
前記心拍信号形成処理部は、
前記心拍波形から所定狭帯域の周波数成分を通過させて外乱ノイズ成分又は体動ノイズ成分あるいは、外乱ノイズ成分と体動ノイズ成分の両方を除去するバンドパスフィルタを備え、
前記信号増幅部は、前記バンドパスフィルタを通過した周波数成分を信号増幅することを特徴とする請求項２に記載の心拍計。
【請求項４】
前記信号増幅部は、
前記微分特性を２階以上の微分特性で前記心拍波形の周波数成分を信号増幅することを特徴とする請求項２又は３に記載の心拍計。
【請求項５】
前記心拍検出処理部は、
前記心拍信号の波形スロープ成分を検出し、所定の波形スロープ成分であるとき心拍を検出することを特徴とする請求項１から４の何れか一つに記載の心拍計。
【請求項６】
前記心拍検出処理部は、
前記心拍信号の波高値と傾きを検出条件として、心拍信号に付加された波形歪みを検出する波形歪み検出部を備え、この波形歪み検出部で検出した波形歪みを心拍として検出することを特徴とする請求項１から５の何れか一つに記載の心拍計。
【請求項７】
生体の心拍波形を検出する心拍波形検出工程と、
前記心拍波形の主成分以外の周波数帯域にあるノイズ成分を除去するフィルタ処理と、前記フィルタ処理後の前記心拍波形の低域側の周波数成分を低増幅率で増幅し、心拍波形の高域側の周波数成分を高増幅率で増幅する信号増幅処理とを施して心拍信号を形成する心拍信号形成工程と、
前記心拍信号形成工程において心拍信号に付加された波形歪みを検出する波形歪み検出工程と、
前記波形歪みを心拍として検出する心拍検出工程とを備えることを特徴とする心拍検出方法。
【請求項８】
生体の心拍波形を検出する心拍波形検出工程と、
前記心拍波形の主成分以外の周波数帯域にあるノイズ成分を除去するフィルタ処理と、前記フィルタ処理後の周波数成分を微分特性で信号増幅する信号増幅処理とを施して心拍信号を形成する心拍信号形成工程と、
前記心拍信号形成工程において心拍信号に付加された波形歪みを検出する波形歪み検出工程と、
前記波形歪みの検出により心拍波形の基本波又は２倍高調波あるいは、基本波と２倍高調波の両方を心拍として検出する心拍検出工程とを備えることを特徴とする心拍検出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、心拍計および心拍検出方法に関し、特に心拍波形に含まれる体動ノイズや外乱ノイズのノイズ除去に関する。
【背景技術】
【０００２】
生体の心拍を測定する心拍計は従来より種々提案されている。発光素子から身体に光を照射し、その反射光または透過光を受光素子で検出し、受光信号を脈拍信号に変換して脈拍を検出する脈拍検出回路として、例えば、特許文献１が知られている。【０００３】
このような心拍計では、雑音に対して脈拍数が安定して表示されることが求められており、脈拍数表示の安定性を高める提案されている（例えば、特許文献２、特許文献３、特許文献４）。【０００４】
特許文献２は、雑音に起因するパルスが比較的幅のせまいパルスであることに着目し、脈波検出回路と脈波演算手段との間に、脈波検出回路が出力する脈波信号のパルス幅を評価する脈波パルス評価手段を設け、正規の脈波信号と認められる信号のみを脈波演算手段に伝達することによって、脈拍表示の安定性を得ようとするものである。【０００５】
また、特許文献３には、加速度センサとウエーブレット変換を用いて体動を除去するものであり、加速度センサによって、体動を加速度として体動波形を検出し、この体動波形にウエーブレット変換を施して周波数領域毎の体動解析データを生成し、また、生体の検出部位から検出した脈波波形にウエーブレット変換を施して周波数領域毎の脈波解析データを生成し、脈波解析データから体動解析データを減算して脈拍を検出することが開示されている。【０００６】
また、心拍情報を光学的に取り出すフォトプレチスムグラフィー（PPG）において、PPG信号から高周波成分の雑音をウエーブレット変換を用いて除去する点については、特許文献４にも示されている。【０００７】
【特許文献１】特公昭６１−２９７３０号公報
【特許文献２】特公平４−７９２５０号公報
【特許文献３】特開平１１−９５６４号公報
【特許文献４】特開２００３−３１０５６２号公報
【非特許文献１】順天堂大学名誉教授  石河利寛著、「健康・体力のための運動生理学」、杏林書院、2000年4月, P.75-78
【発明の開示】
【発明が解決しようとする課題】
【０００８】
心拍計による心拍数の測定は、例えば非運動時や運動時というように様々な測定条件で行われる。このように異なる測定条件下では、心拍センサで検出される心拍波形は、基本波形に限らずこの基本波形に重畳されるノイズ成分についても、その周波数や波高値などの信号特性が異なる。【０００９】
図２４は心拍センサ等の心拍検出部で検出される心拍波形の一例を示している。心拍波形の主成分は基本波と高調波である。心拍波形は、この主成分に加えてノイズ成分が重畳され、多数の周波数成分を含んでいる。【００１０】
例えば、基本波形は基本周波数成分の他に２倍高調波成分をはじめとする高調波成分を含み、ノイズ成分は、心拍センサや伝送系中に侵入する電気的な高周波ノイズである外乱ノイズや、心拍センサを装着した非測定者が運動した際に心拍センサの装着位置ずれ等で発生する微動ノイズである体動ノイズを含んでいる。【００１１】
図２５は図２４の心拍波形に含まれるノイズ成分の一例であり、図２４（ａ）は外乱ノイズの一例を示し、図２５（ｂ）は体動ノイズの一例を示している。【００１２】
これら基本波や高調波の周波数や波高値は、被測定者が心拍計を使用する状態に応じて変化し、例えば、非運動時、運動開始時、運動中等の各状態で異なる特性を示す。【００１３】
図２６，図２７，図２８は、それぞれ非運動時，運動開始時，及び運動中における周波数特性を示している。心拍波形は、上記したように、各状態において基本波や２倍高調波等の高調波の主成分に加えて外乱ノイズや体動ノイズを含み、運動中には、基本波及び２倍高調波の周波数は非運動時よりも高くなる他、低い周波数成分が増加する（図２８参照）。また、運動開始時には、基本波や２倍高調波等の高調波の主成分と比較して、高周波成分が増加する（図２７参照）。【００１４】
ここで、人体の心拍数は、おおよそ20〜220bpm（一分当たりの心拍数）であり、このときの基本周波数は0.33（＝20bpm/60sec）〜3.67（＝220bpm/60sec）Hzの周波数帯域となり、２倍高調波までの高調波を考慮すると0.33（＝20bpm/60sec）〜7.34（3.67×2）Hzの周波数帯域となる。なお、心拍数について記載する文献としては、例えば、非特許文献１がある。この非特許文献１には、運動時の最大心拍数を、運動時の最大心拍数＝220−年齢という関係式で表しており、220bpmを上限としている。他方、下限について、非特許文献１には、最も低い実測例として、安静時心拍数が28bpmのスキー選手が確認されている旨の記載がある。【００１５】
一般的に、非運動時は例えば心拍数が80bpm以下の低い値であって、基本周波数は1.34Hz以下でその振幅も小さい。非特許文献１には、非運動時において、測定した人の90％以上が心拍数50〜80bpmの範囲に入っているという実測データが記載されている。これに対して、運動中には心拍数は80bpmよりも高い値となると共に基本周波数も高くなり、また、運動による血流増加に伴った信号振幅も例えば1.5倍以上に増加する。【００１６】
心拍波形に含まれるノイズ成分を除去するには、一般的には、所定の周波数幅を通過帯域とするバンドパスフィルタ等によって、主成分が含まれる周波数帯域の信号成分を抽出することが考えられる。【００１７】
図２９は非運動時の例（基本周波数を1.34Hzとする）を示し、図３０は運動中の例（基本周波数を3.0Hzとする）を示し、0.33（＝20bpm/60sec）〜3.67（＝220bpm/60sec）Hzの通過帯域を有するバンドパスフィルタによりノイズ除去を行う例を示している。【００１８】
上記したように、人体の心拍数は0.33（＝20bpm/60sec）〜3.67（＝220bpm/60sec）Hzを主成分の周波数帯域と考えることができるため、図２９，図３０に示すように、心拍波形の周波数成分（図２９（ａ），図３０（ａ））を0.33（＝20bpm/60sec）〜3.67（＝220bpm/60sec）Hzの通過帯域のバンドパスフィルタに通すことで、心拍波形の主成分である基本周波数やその２倍高調波を抽出することができる。例えば、非運動時には基本周波数（1.34Hz）とその２倍高調波（2.67Hz）を抽出することができ（図２９（ｃ））、運動時には基本周波数（3.0Hz）を抽出することができる（図３０（ｃ））。【００１９】
しかしながら、ノイズ成分の周波数帯域は3〜30Hzであるため、心拍波形の主成分の周波数帯域とノイズ成分の周波数帯域とで、周波数帯域に重なる部分（例えば、3.0Hz〜3.67Hz）が存在する。そのため、この周波数帯域の重なり部分にノイズ成分がある場合には、このノイズ成分を除去することができない（例えば、図２９（ｃ）中のノイズ成分）という問題がある。【００２０】
このように、心拍波形に含まれるノイズ成分の除去を、主成分が含まれる周波数帯域の信号成分を抽出し、その他の周波数帯域の信号成分を除去する信号処理によって行う場合には、帯域内に存在するノイズ成分を除去することができないという問題がある。この問題は、上記したバンドパスフィルタによる処理に限らず、前記したウエーブレット変換を用いるといった周波数解析による手法においても同様と考えられる。【００２１】
また、図３０に示すように、運動中には低い周波数のノイズ成分が増えるため、バンドパスフィルタの通過帯域内に発生したノイズ成分を除去することができない（図３０（ｃ））。また、この運動中では基本周波数が高くなることで２倍高調波はバンドパスフィルタの通過帯域外となって除去され、主成分は基本周波数のみとなる。このように、運動中では、ノイズ成分が残り、さらに、２倍高調波が除去されて主成分は基本周波数のみとなるため、Ｓ／Ｎ比が低下するという問題がある。【００２２】
上記のように、測定条件が変化すると、心拍波形に含まれる基本波や高調波やノイズ成分の周波数や波高値などの信号特性が変化するため、ノイズ成分を除去して心拍波形に含まれる主成分を有効に検出することが困難となり、心拍の正確な検出が困難となるという問題がある。【００２３】
そこで、本発明は従来の問題を解決し、測定条件が変化し、心拍波形に含まれる基本波や高調波やノイズ成分の周波数や波高値などの信号特性が変化した場合であっても、心拍を正確に検出することを目的とする。
【課題を解決するための手段】
【００２４】
本発明の心拍計は、生体の心拍波形を検出する心拍波形検出部と、心拍波形から心拍を検出するための信号処理部とを備える。【００２５】
信号処理部は、心拍波形に対して信号処理を施して心拍信号を形成する心拍信号形成処理部と、この心拍信号形成処理部により形成された心拍信号から心拍を検出する心拍検出処理部とを備える。【００２６】
本発明の心拍信号形成処理部は、所定の周波数特性を付加する信号処理を行うことによって、周波数成分の信号強度を増幅させた心拍信号を形成する。この周波数成分の信号強度の増幅は、心拍波形の主成分を増加させ、ノイズ成分を低減させるものである。【００２７】
心拍信号形成処理部は、心拍波形の周波数成分を信号増幅する信号増幅部を有し、この信号増幅部は、心拍波形の低域側の周波数成分を低増幅率で増幅し、心拍波形の高域側の周波数成分を高増幅率で増幅する信号増幅特性とすることによって、心拍波形に所定の周波数特性を付加し、心拍波形の主成分を増加させ、ノイズ成分を低減させる。この信号処理としては、例えば、心拍波形の周波数成分を微分特性で信号増幅する。【００２８】
また、心拍信号形成処理部は帯域通過のフィルタ処理を含むことができ、これによって心拍波形の所定狭帯域の信号に対して信号処理を行い、心拍波形の主成分以外の周波数帯域にあるノイズ成分を除去する。【００２９】
この心拍信号形成処理部で行うノイズ成分の除去及び信号強度の増幅の信号処理は、心拍信号に波形歪みを付加する。本発明による波形歪みは、心拍波形の主成分を増加し、ノイズ成分を低減する。これにより、波形歪みは心拍と強い相関で関係付けられる。したがって、波形歪みは心拍情報を良好なＳ／Ｎ比で備えることになる。【００３０】
本発明の心拍検出処理部は、心拍信号形成処理部により心拍信号に付加された波形歪みを検出することで、ノイズの影響を低減して、誤検出することなく心拍を検出することができる。【００３１】
なお、心拍検出部は心拍センサと検出回路で構成することができる。心拍センサは、例えば光学式の反射センサあるいは透過センサを用いることができ、検出回路は、心拍センサが光信号を出力する場合には電気信号に変換し、必要に応じて信号増幅を行う。【００３２】
本発明の心拍信号形成処理部は、心拍信号が備える周波数内の所定の周波数帯域に対して信号処理を施す狭帯域信号処理を施す。【００３３】
心拍信号形成処理部は、心拍波形の周波数成分から所定狭帯域の周波数成分を通過させて外乱ノイズ成分又は体動ノイズ成分あるいは、外乱ノイズ成分と体動ノイズ成分の両方を除去するバンドパスフィルタと、このバンドパスフィルタを通過した周波数成分を、高周波側を高い増幅率で増幅する信号増幅部とを有する。【００３４】
バンドパスフィルタは、心拍波形の主成分以外の周波数帯域にあるノイズ成分を除去する。バンドパスフィルタの通過帯域は、心拍の第１の状態における最高心拍波形の２倍高調波の周波数以上の周波数を高域側の遮断周波数とし、また、心拍の第１の状態における最低心拍波形の基本周波数以下の周波数を低域側の遮断周波数とする周波数範囲とする。【００３５】
信号増幅部は、例えば微分特性による増幅によって、高周波側の周波数成分を強調して増幅する。微分特性は、１階微分特性や２階微分特性の他、さらに高次の微分特性とすることができるが、心拍波形のＳ／Ｎ比や、信号処理部を形成するハードウエア構成やその費用対効果を考慮すると、４０dBの増幅率が得られる２階微分特性が適当である。この２階微分特性等の微分特性によって、心拍信号の高周波側の周波数成分は強調されて増幅される。【００３６】
信号増幅部は、心拍の第２の状態における最高心拍波形の基本周波数の付近の増幅度を最も高くする増幅特性を備え、この周波数付近の成分をより強調して増幅する。これによって、ノイズ成分については低い増幅率で増幅し、心拍の基本周波数や２倍高調波の周波数の成分については高い増幅率で増幅し、ノイズの影響を低減する。【００３７】
また、信号増幅部の増幅特性は、最も高い増幅度の周波数よりも高い周波数域側において、増幅率が漸次低下する周波数帯域を備える。これにより、２倍高調波の周波数が変動した場合についても、信号成分を取りこぼすことなく有効に取得することができ、また、高次高調波成分についても取り込むことができる。【００３８】
したがって、本発明の心拍信号形成処理部は、心拍の第１の状態における最低心拍波形の基本周波数以下の周波数を低域側の遮断周波数とし、心拍の第１の状態における最高心拍波形の２倍高調波以上の周波数を高域側の遮断周波数とする周波数範囲と、第２の状態における最高心拍波形の基本周波数の周波数付近の増幅度を最も高くする増幅特性とを組み合わせた周波数特性を備える。【００３９】
上記する周波数特性において、周波数範囲は第１の心拍状態で定まる周波数によって周波数範囲を設定するものであり、心拍の第１の状態を心拍数が少なく基本周波数が低い状態とし、心拍の第２の状態を心拍数が多く基本周波数が高い状態とする。また、心拍波形の主成分は、基本周波数と２倍高調波とする。【００４０】
ここで、第１の状態において取り得る周波数には幅があり、測定対象や測定条件によって心拍波形の基本周波数及び２倍高調波は所定の周波数幅内で変動する。【００４１】
この周波数幅内において、最低の心拍波形の基本周波数以下の周波数を低域側の遮断周波数とすることによって、主成分を含む低周波数の内で想定し得る最も低い周波数を低周側波数の目安として設定することができる。【００４２】
また、同じく基本周波数及び２倍高調波の周波数幅内において、最高の心拍波形の２倍高調波の周波数以上を高域側の遮断周波数とすることによって、主成分を含む高周波数の内で想定し得る最も高い周波数を高域側周波数の目安として設定することができる。【００４３】
したがって、この通過帯域は、この最低の心拍波形の基本周波数や最高の心拍波形の２倍高調波の周波数を基準とする帯域に幅を持たせた周波数帯域である。【００４４】
上記のように周波数範囲を設定することで、心拍波形の周波数成分からこの周波数範囲の周波数成分を検出すれば、その中には必ずその心拍波形の基本周波数あるいは２倍高調波の主成分が含まれることになるため、検出した主成分から心拍信号を求めることが可能となる。【００４５】
また、上記する周波数特性において、増幅特性は高周波側の周波数成分を、低周波数側の周波数成分より大きく増幅するものであり、心拍の第２の状態における最高心拍波形の基本周波数付近の増幅度を最も高く設定する。心拍の第２の状態を、単位時間当たりの心拍数が多く、基本周波数が高い状態とすることによって、主成分を含む高周波数の内で想定し得る最も高い基本周波数成分について高い増幅度で増幅することができる。【００４６】
上記のように増幅特性を設定することで、第２の状態において基本周波数が高周波数となり、２倍高調波が前記した周波数範囲から外れ、周波数範囲内に含まれる心拍波形の主成分が基本周波数のみとなった場合であっても、基本周波数を確実に増幅することができると共に、その信号を他の周波数の信号と比較して大きく増幅することができるため、心拍情報のＳ／Ｎ比を高めることができ、検出した主成分から心拍信号を求めることが可能となる。【００４７】
心拍の第１の状態は非運動状態とし、また、心拍の第２の状態は運動状態とすることができる。【００４８】
前記した心拍信号形成処理部では高周波数側の周波数成分を強調して増幅する処理を行うため、形成された心拍信号はスロープ成分が強調される。スロープ成分は、心拍波形の主成分に基づいて形成されるものであるから、スロープ部分の個数と心拍数とは対応するものとなる。本発明の心拍検出処理部は、このスロープ成分を検出することによって、心拍を検出する。【００４９】
本発明の心拍検出処理部の一態様は、心拍信号の波形スロープ成分を検出し、所定の波形スロープ成分であるとき心拍を検出する。なお、スロープ成分は、心拍信号波形の立ち上がり期間、あるいは立ち下がり期間における信号変化の何れかを用いてもよく、また、両方の信号変化を用いてもよい。この場合の検出数は心拍数の２倍となるため、検出数を１／２とすることで心拍数を求める。【００５０】
本発明の心拍検出処理部の他の態様は、心拍信号の波高値と傾きを検出条件として、心拍信号に付加された波形歪みを検出する波形歪み検出部を備える。この波形歪み検出部は、波形歪みを心拍として検出する。【００５１】
ここで、検出条件は、心拍信号が所定波高値以上であり、かつ、心拍信号の傾きが所定期間において単純増加又は単純減少あるいは、単純増加と単純減少の両方である。【００５２】
所定波高値は、心拍の非運動時における最高心拍波形の２倍高調波の周波数付近での周波数成分とすることができる。通常、心拍波形の基本周波数の波高値は、２倍高調波の波高値よりも大きく、また、運動中の波高値は非運動時の波高値よりも大きくなる。また、一般に、心拍波形の２倍高調波の波高値はノイズ成分の波高値よりも大きい場合が多い。【００５３】
そこで、非運動時における２倍高調波の周波数付近での周波数成分よりも大きな波高値を検出条件とすることによって、心拍波形の主成分に有無を判定することができる。【００５４】
また、所定期間は、心拍波形の周期の１／４以下の周期幅とすることができる。心拍波形の主成分は、少なくとも基本周波数と２倍高調波を含む場合、その周期特性は長くとも心拍波形の周期の１／４と想定されるため、心拍の主成分は、少なくとも心拍波形の周期の１／４以下の周期幅に渡って単調に増加あるいは単調に減少する。この周期幅は、前記したスロープ部分の幅に相当する。【００５５】
ここで、心拍信号が心拍波形をサンプリングすることで得られる場合には、心拍波形の周期の１／４の時間幅をサンプリング周期で除算した値に相当する個数の連続するサンプリング点を検出し、そのサンプリング点の値が心拍の非運動時における最高心拍波形の２倍高調波の周波数付近での周波数成分よりも大きく、かつ、単調に増加あるいは単調に減少する場合に、波形歪みを心拍として検出する。【００５６】
本発明は、上記した心拍計の態様の他に、心拍波形のノイズを除去し、心拍を検出する方法の態様を含む。【００５７】
本発明の心拍検出方法は、生体の心拍波形を検出する心拍波形検出工程と、心拍波形の低域側の周波数成分を低増幅率で増幅し、心拍波形の高域側の周波数成分を高増幅率で増幅する信号増幅処理を施して心拍信号を形成する心拍信号形成工程と、心拍信号形成工程において心拍信号に付加された波形歪みを検出する波形歪み検出工程と、波形歪みを心拍として検出する心拍検出工程とを備える。【００５８】
また、心拍信号形成工程は、心拍波形の周波数成分を微分特性で信号増幅する信号増幅処理とすることができ、心拍検出工程は、波形歪みの検出により心拍波形の基本波又は２倍高調波あるいは、基本波と２倍高調波の両方を心拍として検出することができる。【００５９】
ここで、心拍信号形成工程の一態様は、心拍波形の周波数成分から所定狭帯域の周波数成分を通過させて外乱ノイズ成分又は体動ノイズ成分あるいは、外乱ノイズ成分と体動ノイズ成分の両方を除去する帯域通過処理と、狭帯域内の周波数成分を、高周波側を高い増幅度とする２階微分特性等の微分特性で信号増幅する信号増幅処理とを含む狭帯域信号処理の２つの処理を含む。【００６０】
また、心拍信号形成工程の別の態様は、帯域通過処理は、心拍の第１の状態における最低心拍波形の基本周波数以下の周波数を低域側の遮断周波数とし、心拍の第１の状態における最高心拍波形の２倍高調波の周波数以上の周波数を高域側の遮断周波数とする周波数範囲に制限することで行う。また、信号増幅処理は、第２の状態における最高心拍波形の基本周波数付近の増幅度を最も高くする増幅特性で信号増幅することで行う。これらの処理を組み合わせることで、周波数範囲において第２の状態における最高心拍波形の基本周波数付近の心拍波形を強調して信号増幅する。【００６１】
帯域通過処理は、非運動状態において測定して最低心拍波形の基本周波数を低域側周波数とし、運動状態において測定した最高心拍波形の基本周波数を高域側周波数とする周波数範囲に制限する。【００６２】
本発明の波形歪み検出工程の一態様は、心拍信号の波形スロープ成分を検出し、所定の波形スロープ成分を抽出することによりノイズ成分を除去して心拍のみを検出する。【００６３】
スロープ成分は、心拍信号波形の立ち上がり期間あるいは立ち下がり期間における信号変化、もしくはその両方期間における信号変化としてもよい。【００６４】
本発明の波形歪み検出工程の他の態様は、心拍信号の波高値と傾きを検出条件として心拍信号に付加された波形歪みを検出し、検出した波形歪みを心拍として検出する。【００６５】
ここで、検出条件は、前記したように、心拍信号が所定波高値以上であり、かつ、所定期間において単純増加又は単純減少あるいは、単純増加と単純減少の両方とし、所定波高値は、心拍の非運動状態における最高心拍波形の２倍高調波の周波数付近での周波数成分とする。
【発明の効果】
【００６６】
本発明によれば、測定条件が変化し、心拍波形に含まれる基本波や高調波やノイズ成分の周波数や波高値などの信号特性が変化した場合であっても、心拍を正確に検出することができる。
【図面の簡単な説明】
【００６７】
【図１】本発明の概略を説明するための図である。
【図２】心拍波形、基本周波数及び２倍高調波の概略を説明するための図である。
【図３】本発明の心拍計の概略構成を説明するための図である。
【図４】心拍センサの一構成例を説明するための概略断面図である。
【図５】本発明の心拍計の別の構成を説明するための図である。
【図６】心拍センサの別の構成例を説明するための概略断面図である。
【図７】バンドパスフィルタ及び２階微分特性処理による狭帯域信号処理の詳細説明図である。
【図８】心拍の第１の状態（非運動時）におけるバンドパスフィルタ及び２階微分特性による動作を説明するための図である。
【図９】心拍の第２の状態（運動中）におけるバンドパスフィルタ及び２階微分特性による動作を説明するための図である。
【図１０】非運動時における処理前の波形及び周波数成分を説明するための図である。
【図１１】非運動時における処理後の波形及び周波数成分を説明するための図である。
【図１２】運動開始直後における処理前の波形及び周波数成分を説明するための図である。
【図１３】運動開始直後における処理後の波形及び周波数成分を説明するための図である。
【図１４】運動中における処理前の波形及び周波数成分を説明するための図である。
【図１５】運動中における処理後の波形及び周波数成分を説明するための図である。
【図１６】狭帯域信号処理の一構成例を説明するための図である。
【図１７】本発明の信号増幅の増幅特性を説明するための図である。
【図１８】本発明の波形歪み検出（スロープ検出）を説明するためのフローチャートである。
【図１９】周期とスロープ部分との関係を説明するための図である。
【図２０】基本周波数によるは波形歪みを説明するための図である。
【図２１】サンプリング点によるスロープ検出を説明するための図である。
【図２２】波形歪み検出の実例図を用いて説明するための図である。
【図２３】波形歪み検出部の一構成例を説明するための図である。
【図２４】心拍センサ等の心拍検出部で検出される心拍波形の一例を示す図である。
【図２５】図２４の心拍波形に含まれるノイズ成分の一例を示す図である。
【図２６】非運動時における周波数特性を示す図である。
【図２７】運動開始時における周波数特性を示す図である。
【図２８】運動中における周波数特性を示す図である。
【図２９】非運動時の例（基本周波数を1.34Hzとする）を示す図である。
【図３０】運動中の例（基本周波数を3.0Hzとする）を示す図である。
【符号の説明】
【００６８】
１  心拍計
２  心拍センサ
２ａ  発光素子部
２ｂ  受光素子部
２ｃ  遮光部
２Ａ  触覚センサ
３  検出回路
３ａ  受光回路部
３ｂ  増幅回路部
３ｃ  Ａ／Ｄ変換部
４  信号処理部
５  狭帯域信号処理部
５a  バンドパスフィルタ
５ｂ  ２階微分特性
６  スロープ検出部
６ａ  波高値
６ｂ  傾斜
７  心拍計数部
８  心拍報知部
９  発光回路部
１０  生体
１１  血管
１２  振動波
１３  生体組織
１４  皮膚
２０  心拍波形
２１  基本周波数
２２  ２倍高調波
５０  遅延器
５１  係数器
５２  加算器
６０  記憶手段
６１  第１の比較手段
６２  第２の比較手段
６３  カウンタ
【発明を実施するための最良の形態】
【００６９】
以下、本発明の心拍計及び心拍波形からノイズ成分を除去して心拍を検出する手順について図を用いて詳細に説明する。【００７０】
はじめに、本発明の概略について図１を用いて説明する。図１において、本発明の心拍計１は、生体の心拍波形を検出する心拍波形検出部Ａと、検出した心拍波形を信号処理した心拍を検出する信号処理部４とを備える。ここで、信号処理部４は、心拍信号形成処理部Ｂと心拍検出処理部Ｃを備える。信号処理部４で検出した心拍は、心拍数検出部Ｄで心拍数を検出する。ここで、心拍数検出部Ｄは、心拍数計数部７において心拍数を計数し、心拍数報知部８において、計数した心拍数を表示、送信、記録等により報知する。【００７１】
なお、前記した心拍波形検出部Ａは、例えば、心拍センサ２と、心拍センサ２の出力から検出信号を取得する検出回路３とにより構成することができる。心拍センサ２は、例えば、光学式センサを用いることができ、検出回路３は光信号等の心拍センサ２で取得された出力を電気信号に変換し、必要に応じて信号増幅したり、Ａ／Ｄ変換によりデジタル信号に変換することができる。【００７２】
信号処理部４は、心拍波形の所定狭帯域の信号に対してフィルタ処理及び信号増幅処理を行って心拍信号を形成する心拍信号形成処理部Ｂと、この心拍信号形成処理部Ｂによる所定の周波数特性の付加によって、心拍信号に付加された波形歪みから心拍を検出する心拍検出処理部Ｃとを備える。【００７３】
心拍信号形成処理部Ｂは、狭帯域信号処理部５によって構成することができる。狭帯域信号処理部５は、２つの処理内容５ａ，５ｂを含んでいる。【００７４】
第１の処理内容５ａは、心拍波形の主成分以外の周波数帯域にあるノイズ成分を除去するものであり、取り込んだ心拍波形の周波数帯域から所定狭帯域の信号を抽出することで行うことができる。心拍波形は、主成分として基本周波数とそのｎ倍高調波の周波数成分を含んでいる。ｎ倍高調波は“ｎ”の数を多くとるほど心拍波形に近づくが、一般にｎ倍高調波の波高値は“ｎ”が大きいほど小さいため、心拍数を検出することを目的とする場合などの比較的に低い精度で十分な場合には、“ｎ”として“２”を用いた２倍高調波で十分である。【００７５】
したがって、ここでは、心拍波形の主成分は基本周波数と２倍高調波の周波数成分とする。図２は心拍波形、基本周波数及び２倍高調波の概略を説明するための図である。図２（ａ）は心拍波形２０の一例であり、図２（ｂ）は心拍波形２０に含まれる基本周波数２１の周波数成分と２倍高調波２２の周波数成分を示している。なお、基本周波数や２倍高調波はフーリエ解析等の周波数解析で求めることができる。【００７６】
周波数軸上において、心拍波形の主成分である基本周波数と２倍高調波の周波数成分の他にノイズによる周波数成分が存在する。【００７７】
ノイズとしては、前記したように心拍センサや回路上で侵入する外乱ノイズや、心拍センサが被測定位置からずれることで発生する位置計体動ノイズがある。ノイズに含まれる周波数はノイズ源等によって様々であり、心拍波形の主成分が存在する周波数範囲の内外にわたって現れる。【００７８】
狭帯域信号処理部５は、心拍波形の主成分である基本周波数や２倍高調波の周波数成分が現れる周波数範囲を設定し、この周波数範囲を用いて心拍波形を狭帯域化して当該周波数範囲内の信号のみの狭帯域信号を生成する。この心拍波形の狭帯域化により、周波数範囲の外側の周波数を持つノイズ成分を除去することができる。【００７９】
このノイズ成分を除去する第１の処理内容５ａは、例えば、バンドパスフィルタで実現することができる。【００８０】
次に、狭帯域信号処理部５が持つ第２の処理内容５ｂは、周波数成分の信号強度を増幅させ、これによって心拍信号を形成するものである。この信号処理は、心拍波形に所定の周波数特性を付加する信号処理を施すことによって行うことができる。この所定の周波数特性を付加する信号処理は、心拍波形の主成分を増加させ、ノイズ成分を低減させるものである。【００８１】
前記した狭帯域信号処理によってノイズ成分を除去した場合であっても、基本周波数と２倍高調波の周波数成分の周波数範囲内で現れるノイズ成分は除去されずに残存することになり、心拍信号のＳ／Ｎ比を低下させる要因となる。【００８２】
狭帯域信号処理部５の第２の処理は、このような、同じ周波数範囲内に存在する心拍波形の主成分とノイズ成分について、主成分の波高値を信号増幅させることによってノイズ成分との信号強度差を大きくし、これによって心拍信号のＳ／Ｎ比を向上させるものである。ここでは、２倍高調波の周波数成分を信号増幅する例を示している。【００８３】
この主成分の波高値を信号増幅させる第２の処理内容５ｂは、例えば、主に２倍高調波を信号増幅する微分特性処理で実現することができる。微分特性処理は、例えば２階微分特性処理とすることができる。なお、２階微分特性処理は、後述するように微分特性処理の一例であり、また微分特性処理も、高い周波数域において高い増幅率で信号増幅する一例である。【００８４】
この狭帯域信号処理５で行うノイズ成分の除去及び信号強度増幅は、心拍信号に波形歪みを付加することになる（図１中のＥ）。狭帯域信号処理５で行う信号強度増幅は、主成分の周波数成分、特に基本周波数及び２倍高調波を非線形に信号増幅するため、この信号増幅で得られた心拍信号に含まれる周波数成分の強度比率は、当初の心拍波形に含まれている周波数成分の強度比率とは異なることになる。このことは、心拍信号に波形歪みを付加することを意味している。【００８５】
この波形歪みは心拍波形を当初の信号形状から歪ませることになるが、心拍波形の主成分が増加され、ノイズ成分が低減されるため、波形歪みは心拍波形よりもより強く、心拍と強い相関で関係付けられることになる。したがって、狭帯域信号処理５により得られる波形歪みを検出することによって、心拍情報を良好なＳ／Ｎ比で求めることになる。【００８６】
なお、バンドパスフィルタ及び２階微分特性処理による狭帯域信号処理の詳細例については、図７を用いて後述する。【００８７】
心拍検出処理部Ｃは、波形歪み検出部６によって構成することができる。波形歪み検出部６は、例えば心拍信号中に含まれる所定スロープ成分を検出し、このスロープと心拍と対応付けることによって心拍を検出する。スロープ成分は、心拍信号波形の立ち上がり期間での信号変化から求めても、立ち下がり期間での信号変化から求めても良く、あるいは、その両期間での信号変化から求めても良い。【００８８】
波形歪み検出部６は、心拍信号の波高値６ａと傾き６ｂを検出条件として心拍信号に付加された波形歪みを検出する。【００８９】
検出条件は、心拍信号が所定波高値以上であり、かつ、心拍信号の傾きが所定期間において単純増加又は単純減少あるいは、単純増加と単純減少の両方である。【００９０】
所定波高値は、例えば、心拍の非運動時における最高心拍波形の２倍高調波の周波数付近での周波数成分とすることができる。【００９１】
心拍波形について見ると、通常、基本周波数の波高値は２倍高調波の波高値よりも大きく、また、運動中の波高値は非運動時の波高値よりも大きい、また、心拍波形の２倍高調波の波高値はノイズ成分の波高値よりも大きい場合が多い。【００９２】
そこで、波高値を第１の条件として心拍の主成分を検出する。非運動時における２倍高調波の周波数付近での周波数成分の波高値を閾値とすることで、この波高値よりも小さい周波数成分には心拍の主成分は含まれず、この波高値よりも大きい周波数成分には心拍の主成分は含まれている可能性がある。【００９３】
ノイズ成分の大きさによっては前記した波高値を超える場合があるため、上記の波高値による心拍波形の主成分の判定のみによって心拍の主成分を検出することはできない。【００９４】
そこで、所定期間における傾きを第２の条件とすることで心拍の主成分を検出する。心拍波形の主成分が少なくとも基本周波数と２倍高調波を含む場合には、その周期特性は長くとも心拍波形の周期の１／４と想定される。このとき、心拍の主成分は、少なくとも心拍波形の周期の１／４以下の周期幅に渡って単調に増加あるいは単調に減少する。【００９５】
したがって、心拍波形の周期の１／４以下の周期幅の所定期間内において、信号が単調に増加あるいは単調に減少することを傾きの条件として、心拍の検出を行う。【００９６】
心拍信号がサンプリングデータである場合には、心拍波形の周期の１／４の時間幅をサンプリング周期で除算した値に相当する個数の連続するサンプリング点を求め、このサンプリング点の値が心拍の非運動時における最高心拍波形の２倍高調波の周波数付近での周波数成分よりも大きく、かつ、単調に増加あるいは単調に減少する場合に、心拍として検出する。【００９７】
次に、本発明の心拍計の概略構成について図３を用いて説明する。なお、ここでは、心拍センサ２として光学式センサの例について示している。【００９８】
心拍計１は、生体１０から心拍情報を取得する心拍センサ２と、心拍センサ２の出力から検出信号を形成する検出回路３と、検出回路３からの検出信号を信号処理して心拍を検出する信号処理部４と、心拍を計数する心拍数計数部７と、計数した心拍数を報知する心拍数報知部８を備える。【００９９】
心拍センサ２は、発光回路９によって駆動された生体１０に光を照射する発光素子部２ａと、生体１０で散乱あるいは反射、又は透過した光を受光する受光素子部２ｂを備える。図４は心拍センサ２の一構成例を説明するための概略断面図であり、生体１０に対して光を照射し、反射した光を検出する構成例を示している。発光素子部２ａと受光素子部２ｂは遮光板２ｃを挟んで対向し、照射点（図示していない）に対して対称の位置に設置されている。ここで、遮光板２ｃは発光素子部２ａから受光素子部２ｂへ直接に入射する光を阻止している。【０１００】
発光素子部２ａから生体１０に照射された光は、生体１０内の組織や血管１１内の血液で散乱し、再び生体１０の外に出射される。生体１０から出射する光の強度は、血流に応じて変動する。本発明の心拍計１は、この光の強度の血流に応じて変動する光強度の変化に基づいて心拍を検出するものである。【０１０１】
検出回路３は、受光素子部２ｂから得られる光信号を受光し、電気信号の検出信号に変換する受光回路部３ａと、検出信号を信号増幅する増幅回路部３ｂと、デジタル信号に変換するＡ／Ｄ変換部３ｃを備える。【０１０２】
信号処理部４は、検出回路３からの検出信号を、前記したバンドパスフィルタや２階微分特性等によってノイズ除去や２倍高調波の信号増幅を行う狭帯域信号処理部５、狭帯域信号処理部５で生成された心拍信号の波形歪みを形成して心拍を検出する心拍検出部（スロープ検出部）６を備える。【０１０３】
心拍数計数部７は信号処理部４で求めた心拍を計数し、心拍数報知部８は心拍数計数部７で計数した心拍数を報知する。なお、報知は、心拍数の表示、記録、送信等を含むものである。【０１０４】
図５は、図３で示した心拍計の別の構成を説明するための図である。この構成例では、心拍センサ２として触覚センサ２Ａの例について示している。【０１０５】
心拍計１は、前記図３で示したように心拍センサ２と、検出回路３と、信号処理部４と、心拍計数部７と、報知部８を備え、心拍センサ２である触覚センサ２Ａの検出信号は、検出回路３中の増幅回路部３ｂで信号増幅され、Ａ／Ｄ変換部３ｃでデジタル信号に変換される。【０１０６】
この構成例では、心拍センサ２として触覚センサ２Ａを備える。ここで、この触覚センサ２Ａは、生体から発せられる振動を検出するセンサを総称しており、例えば、生体中の動脈血管の脈波を検出することで、心拍を検出する。【０１０７】
図６は触覚センサ２Ａの一構成例を説明するための概略断面図である。生体内の動脈血管１１は、血管内を流れる血液の変動に応じて、脈拍に同期して振動する。この動脈血管の振動は、生体の組織１３内を振動波１２となって伝搬する。触覚センサ２Ａは、生体の例えば皮膚表面に触覚して取り付けられ、生体の組織１３内を伝搬した振動波１２を検出する。この振動波の検出は、圧力変化や振動変化として検出される。【０１０８】
触覚センサ２Ａは、振動波を検出する態様に応じた各種センサを用いることができる。例えば、触覚センサ２Ａとして圧力センサを用いた場合には、振動波を圧力変化として検出する。また、触覚センサ２Ａとして振動センサを用いた場合には、振動波を振動変化として検出する。振動変化としては、例えば、振幅や周波数等の変化を検出する。【０１０９】
触覚センサ２Ａは、生体の組織１３内を伝搬した振動波を、皮膚を介して検出する。そのため、触覚センサ２Ａは生体の脈拍を検出する測定部位の近傍に設けられる。触覚センサ２Ａは、測定部位の近傍の皮膚１４に接触させることで検出感度を高めることができ、また、触覚センサ２Ａを皮膚１４に押し当てることによって、より検出感度を高めることができる。【０１１０】
この構成例の触覚センサ２Ａは、血流に応じて変動する振動波を検出し、検出回路３は検出信号を増幅回路部３ｂで信号増幅し、Ａ／Ｄ変換部３ｃでデジタル信号に変換する。【０１１１】
なお、信号処理部４内の狭帯域信号処理部５、波形歪検出部６、および心拍数計数部７の構成や信号処理、また、心拍数報知部８の構成や信号処理、及び各部の作用動作は、前記図３を用いた説明と同様であるため、ここでの説明は省略する。【０１１２】
なお、心拍センサとして光学センサを用いた場合には、心拍信号を検出する部位を狭い範囲で特定することができ、特定位置での心拍状態を検出することができる。また、心拍センサとして触覚センサを用いた場合には、心拍信号を広範囲から取得することができ、心拍センサの取り付けの位置決めに高い精度を不要とすることができ、また、使用中の位置ずれによる検出不良を緩和することができる。【０１１３】
次に、心拍信号形成処理部Ｂにおけるノイズ成分除去５ａと２倍高調波の信号増幅５ｂについて、それぞれバンドパスフィルタと２階微分特性の例を図７を用いて説明する。【０１１４】
図７（ａ）はバンドパスフィルタの通過帯域を示し、図７（ｂ）は２階微分特性の周波数特性を示し、図７（ｃ）はバンドパスフィルタと２階微分特性の周波数特性とを組み合わせて得られる周波数特性を示している。【０１１５】
バンドパスフィルタは、心拍波形の周波数成分から所定狭帯域の周波数成分を通過させて外乱ノイズ成分又は体動ノイズ成分あるいは、外乱ノイズと体動ノイズの両方を除去する帯域通過フィルタである。【０１１６】
図７（ａ）に示すように、バンドパスフィルタは、低域側周波数Ｆ１を低域側の遮断周波数とし、高域側周波数Ｆ２を高域側の遮断周波数とする狭帯域を通過帯域として設定する。【０１１７】
ここで、バンドパスフィルタの低域側周波数Ｆ１は、心拍の第１の状態における最低心拍波形の基本周波数以下の周波数で設定し、高域側周波数Ｆ２は心拍の第１の状態における最高心拍波形の２倍高調波以上の周波数で設定することができる。ここで、心拍の第１の状態は、被測定者が運動していない非運動時とすることができる。この非運動時における心拍数は個人差があるが、例えば、20bpm〜80bpmとすることができる。なお、［bpm］は１分当たりの心拍数を表している。この心拍数はそれぞれ0.33Hzと1.34Hzに相当するため、低域側周波数Ｆ１は0.33Hzとなり、高域側周波数Ｆ２は2.68 Hz（1.34Hz×２）となる。【０１１８】
また、２階微分特性は高周波側を高い増幅度で信号増幅する。２階微分特性は、図７（ｂ）に示すように、低周波数側から高周波側周波数Ｆ３に向かって増幅度が直線的に増加する。なお、高周波側周波数Ｆ３は２階微分の微分上限値を設定する。【０１１９】
ここで、２階微分特性の高周波側周波数Ｆ３は、心拍の第２の状態における最高心拍波形の基本周波数とし、設定することができる。ここで、心拍の第２の状態は、被測定者が運動中とすることができる。この運動中における心拍数は個人差があるが、例えば、220bpmとすることができる。この心拍数は3.67Hzに相当する。【０１２０】
なお、低周波数側では、バンドパスフィルタの低域側周波数Ｆ１付近の増幅度がほぼ１となるように設定している。【０１２１】
図７（ｃ）は、図７（ａ）に示すバンドパスフィルタと図７（ｂ）に示す２階微分特性の周波数特性とを組み合わせて得られる周波数特性を示している。低域側周波数Ｆ１、高域側周波数Ｆ２、及び高周波側周波数Ｆ３の周波数の設定例を前記した数値例に設定した場合には、周波数は、Ｆ１、Ｆ２，Ｆ３の順に高くなるため、低域側周波数Ｆ１付近の増幅度をほぼ１とし、この低域側周波数Ｆ１付近から高域側周波数Ｆ２に向かって直線的に増幅度が上昇して、高域側周波数Ｆ２で最高の増幅度となり、高域側周波数Ｆ２から高周波側周波数Ｆ３に向かって増幅度が減少し、高周波側周波数Ｆ３付近の周波数を超える周波数範囲では増幅度はほぼ１以下となる。また、低域側周波数Ｆ１付近より低い周波数範囲においても増幅度はほぼ１以下となる。【０１２２】
ここで、運動時におけるノイズ成分（体動ノイズや外乱ノイズのノイズ成分）は、およそ3〜30Hzであるため、上記した数値例の周波数特性と重ねると、図７（ｃ）に示すように、運動時におけるノイズ成分の大部分は、バンドパスフィルタの通過帯域外となって除去される。運動時におけるノイズ成分の内の3〜3.67Hz近傍のノイズ成分は、２階微分特性と重なるため、除去されない。この周波数範囲に現れるノイズ成分については、前記したように、本発明の波形歪み検出（スロープ検出）によって、心拍波形の主成分を強調して増幅させることでノイズ成分との波高値に差を設けることによって、Ｓ／Ｎ比を向上させて、心拍検出の精度を高める。【０１２３】
なお、図７に示す周波数特性は前記した数値例に基づくものであるため、異なる数値例を用いた場合には異なる周波数特性を示す。例えば、最高の増幅度を示す周波数部分も多少ずれることになるが、低周波数側から高周波数側に向かって増幅度が増加する傾向は共通である。【０１２４】
なお、高周波側周波数Ｆ３よりも高い周波数では、増幅率を急峻に下げずに漸次低下するように設定する。これによって、心拍信号の基準周波数や２倍高調波が高周波側周波数Ｆ３よりも高い周波数にずれた場合であっても、切り捨てることなく心拍検出に寄与させることができる。また、３倍高調波を心拍検出に寄与させることもできる。【０１２５】
以下、図８を用いて心拍の第１の状態（非運動時）におけるバンドパスフィルタ及び２階微分特性による動作を説明し、図９を用いて心拍の第２の状態（運動中）におけるバンドパスフィルタ及び２階微分特性による動作を説明する。なお、以下に示す各周波数成分の大きさの比率は説明上から任意に設定したものであり、実際の状態を示すものではない。【０１２６】
はじめに、図８を用いて心拍の第１の状態（非運動時）におけるバンドパスフィルタ及び２階微分特性による動作を説明する。【０１２７】
図８（ａ）は、非運動時における検出信号（心拍波形）の周波数成分を表し、基本周波数と２倍高調波が共に、低域側周波数Ｆ１と高域側周波数Ｆ２との間にあるものとする。【０１２８】
図８（ｃ）の信号は、上記した図８（ａ）の周波数成分を図８（ｂ）に示すバンドパスフィルタに通して得られた周波数成分を示している。バンドパスフィルタによって、高域側周波数Ｆ２より高周波数にあるノイズ成分は除去される。【０１２９】
図８（ｅ）の信号は、バンドパスフィルタ通過後の信号を、図８（ｄ）に示す２階微分特性によって信号増幅して得られた周波数成分を示している。２階微分特性によって、高域側周波数Ｆ３付近の周波数成分はより大きく増幅される。図８（ｄ）では、２倍高調波の周波数成分及びノイズ成分が強調されて増幅される。【０１３０】
次に、図９を用いて心拍の第２の状態（運動中）におけるバンドパスフィルタ及び２階微分特性による動作を説明する【０１３１】
図９（ａ）は、運動中における検出信号（心拍波形）の周波数成分を表し、基本周波数と２倍高調波が共に、高域側周波数Ｆ２よりも高周波数側にあるものとする。なお、基本周波数は２階微分特性の高周波数側周波数Ｆ３よりも低周波数とする。【０１３２】
図９（ｃ）の信号は、上記した図９（ａ）の周波数成分を図９（ｂ）に示すバンドパスフィルタに通して得られた周波数成分を示している。バンドパスフィルタによって、高域側周波数Ｆ２より高周波数にあるノイズ成分及び２倍高調波は除去される。【０１３３】
図９（ｅ）の信号は、バンドパスフィルタ通過後の信号を、図９（ｄ）に示す２階微分特性によって信号増幅して得られた周波数成分を示している。２階微分特性によって、高域側周波数Ｆ３付近の周波数成分はより大きく増幅される。図９（ｄ）では、基本周波数の周波数成分及びノイズ成分が強調されて増幅される。【０１３４】
以下の図１０〜図１５に心拍信号形成処理部Ａ（狭帯域信号処理）による信号処理例を示す。図１０，図１１は、非運動時における処理前と処理後を示し、図１２，図１３は運動開始直後における処理前と処理後を示し、図１４，図１５は運動中における処理前と処理後を示し、それぞれ（ａ）は波形を示し、（ｂ）は周波数成分を示している。【０１３５】
はじめに、非運動時の場合について説明する。【０１３６】
図１０（ａ）に示すように、非運動時における狭帯域信号処理前の元波形では、血流が低レベルにあるため心拍波形の振幅が低く、また、非運動であるためノイズのレベルも低い。図１１（ａ）は、この波形に対して狭帯域信号処理を施して得られる波形を示している。この波形は、狭帯域信号処理による波形歪みによって元波形から変形するが、Ｓ／Ｎ比は向上する。なお、図１１（ａ）中の実線は高域側周波数Ｆ２の曲率成分を示し、これ以上の波形部分は心拍の主成分を表している。【０１３７】
また、図１０（ｂ）と図１１（ｂ）の周波数成分を比較すると、狭帯域信号処理によりＦ１〜Ｆ３の周波数範囲の信号レベルが高められている。【０１３８】
次に、運動開始直後の場合について説明する。
図１２（ａ）に示すように、運動開始直後における狭帯域信号処理前の元波形では、血流が低レベルにあるため心拍波形の振幅が低いが、運動が開始されたためノイズのレベルは高くなるため、Ｓ／Ｎ比は悪化する。【０１３９】
図１３（ａ）は、この波形に対して狭帯域信号処理を施して得られる波形を示している。この波形は、狭帯域信号処理により、波形歪みによって元波形から変形するが、ノイズよりも２倍高調波が強調されて増幅されＳ／Ｎ比は向上する。なお、図１３（ａ）中の実線は高域側周波数Ｆ２の曲率成分を示し、これ以上の波形部分は心拍の主成分を表している。【０１４０】
また、図１２（ｂ）と図１３（ｂ）の周波数成分を比較すると、狭帯域信号処理によりＦ１〜Ｆ３の周波数範囲の信号レベルが高められている。【０１４１】
次に、運動中の場合について説明する。【０１４２】
図１４（ａ）に示すように、運動中における狭帯域信号処理前の元波形では、心拍及び血流が増加し、心拍波形の振幅が高まる。【０１４３】
図１５（ａ）は、この波形に対して狭帯域信号処理を施して得られる波形を示している。この波形は、狭帯域信号処理により高周波数側の信号が増幅される。このとき、高周波数側にノイズ成分がある場合には、このノイズ成分の増幅によって、ピーク点がばらつくとった波形歪みが発生する。【０１４４】
また、図１４（ｂ）と図１５（ｂ）の周波数成分を比較すると、狭帯域信号処理によりＦ１〜Ｆ３の周波数範囲のノイズを含む信号レベルが高められている。【０１４５】
次に、図１６を用いて狭帯域信号処理の一構成例を説明する。なお、図１６に示す構成例は、入力信号をサンプリング信号とし、デジタルフィルタで構成する場合を示している。【０１４６】
一般に、デジタルフィルタは、遅延器５０と係数器５１と加算器５２によって構成される。【０１４７】
図１６（ａ）は本発明の狭帯域信号処理を構成するバンドパスフィルタの構成例を示し、図１６（ｂ）は本発明の狭帯域信号処理を構成する２階微分特性の構成例を示している。【０１４８】
図１６（ａ）に示すバンドパスフィルタは、ＦＩＲ型（Finite Impulse Response：有限インパルス応答）の例を示し、狭帯域の通過帯域の周波数Ｆ１，Ｆ２は、サンプリング周期と、係数器５１に定める係数h0〜hmで設定する。【０１４９】
また、図１６（ｂ）に示す２階微分特性は、サンプリング周期と、後退差分の遅延器の設定個数によって、微分上限値Ｆ３を設定する。【０１５０】
上記した例では、信号増幅の増幅特性において２階微分特性の場合を示しているが、２階微分特性に限らず１階微分あるいは、より高次の微分特性とする他、微分特性以外の増幅特性としてもよい。【０１５１】
図１７は、本発明の信号増幅の増幅特性を説明するための図である。図１７（ａ）〜（ｃ）は微分特性による増幅特性の例であり、図１７（ｄ）はステップ状の増幅特性の例である。【０１５２】
図１７（ａ）は１階微分の微分特性を示し、図１７（ｂ）は２階微分の微分特性を示し、図１７（ｃ）は３階微分の微分特性を示している。各微分特性は微分の回数により増幅率が異なり、１階微分では２０dBの増幅率となり、２階微分では４０dBの増幅率となり、３階微分では６０dBの増幅率となる。なお、各dBの数値は、所定の周波数幅に対する増幅率で表している。【０１５３】
何れの階数の微分特性とするかは任意に設定することができるが、ここでは、検出される心拍信号のＳ／Ｎ比の程度や、信号処理の回路構成や、費用対効果等から、例えば４０dBの増幅率が得られる２階微分を用いた例を示している。【０１５４】
また、高い周波数域で高い増幅率を実現する増幅特性として、微分特性の他、図１７（ｄ）に示すようにステップ状の増幅率を設定してもよい。図１７（ｄ）では、周波数Ｆ１付近の増幅率を“１”とし、高い周波数Ｆ３付近の増幅率を高く設定し、この間をステップ状に変化させる例を示している。なお、ステップの段数は任意に設定することができる。【０１５５】
次に、本発明の心拍検出処理部Ｃである波形歪み検出（スロープ検出）について、図１８のフローチャート、図１９の周期とスロープ部分との関係を説明するための図、図２０の基本周波数による波形歪みを説明するための図、図２１のサンプリング点によるスロープ検出を説明するための図、図２２の波形歪み検出の実例図を用いて説明する。【０１５６】
図１８に示す複数のフローチャートでは、波形歪み検出（スロープ検出）を、検出波形について所定周期でサンプリングして得られるサンプリング値を用いて行う(S11)。【０１５７】
図１９において、心拍の周期をＴとしたとき、波形歪み検出は心拍波形の周期Ｔの１／４の時間幅で行う。この時間幅Ｔ／４内で波形歪み検出を行うには、サンプリング周期ｔとしたとき、時間幅Ｔ／４をサンプリング周期ｔで除算した値（（Ｔ／４）／ｔ）に相当する個数の連続するサンプリング点（Ｓ１〜Ｓ５）を検出し、このサンプリング点（Ｓ１〜Ｓ５）の値が心拍の非運動時における最高心拍波形の２倍高調波の周波数付近での周波数成分よりも大きく、かつ、単調に増加あるいは単調に減少する場合に、波形歪みを心拍として検出する。【０１５８】
例えば、運動時の最大心拍Ｆ３を例えば、3.367Hzとしたとき、この周期の１／４以下の時間幅（68msec）でスロープ検出を行う。ここで、サンプリング周期を12msecとすると、5.67サンプリング回数以下のサンプリングで検出される。なお、このサンプリング回数以上である場合には、１／４周期を超えてしまうため、スロープ検出を行うことができなくなる。【０１５９】
そこで、ここでは、サンプリング回数として５回を設定する。S1で求めたサンプリング値を、少なくともスロープの傾斜を求める際に必要な個数分求め、順次S1〜S5に格納する。なお、各サンプリング値は、サンプリングが行われる毎に順次入れ替え、古いサンプリング値から順から廃棄する(S12)。【０１６０】
波形歪みを検出する第１の条件として、サンプリング点（Ｓ１〜Ｓ５）の値が設定値以上であるかを判定する。設定値としては、心拍の非運動時における最高心拍波形の２倍高調波の周波数付近での周波数成分とすることができ、前記したように、この設定値以上であれば、ノイズ成分を除去することができる。(S13)。【０１６１】
次に、波形歪みを検出する第２の条件として、サンプリング点（Ｓ１〜Ｓ５）の傾斜が所定期間で連続して単純増加あるいは単純減少であるかを判定する。【０１６２】
図２０はこの第１の条件及び第２の条件による判定を示している。図２０（ａ）において、Ｈは第１の条件である設定値を示している。【０１６３】
図２０（ａ）の場合には、サンプリング値Ｓ１〜Ｓ５の全てが第１の条件である設定値Ｈ以上であり、サンプリング値Ｓ１〜Ｓ５がＳ１＜Ｓ２＜Ｓ３＜Ｓ４＜Ｓ５の関係にあり、第２の条件を満足している。これより、このスロープ部分によって一心拍と判定する。【０１６４】
図２０（ｂ）の場合には、サンプリング値Ｓ１〜Ｓ５がＳ１＜Ｓ２＜Ｓ３＜Ｓ４＜Ｓ５の関係にあって第２の条件を満足しているが、サンプリング値Ｓ１は設定値Ｈ以下であるため、このスロープ部分による一心拍の判定は行わない。【０１６５】
また、図２０（ｃ）の場合には、サンプリング値Ｓ１〜Ｓ５の全てが設定値Ｈ以上であって第１の条件を満たすが、サンプリング値Ｓ５はＳ４よりも小さく、Ｓ１＜Ｓ２＜Ｓ３＜Ｓ４＜Ｓ５の第２の条件を満足していないため、このスロープ部分による一心拍の判定は行わない(S14)。【０１６６】
S13とS14の第１、第２の条件の両条件を満たした場合には、このスロープ部分によって一心拍があったと判定する(S15)。この心拍数を順次計数し（S16）、計数した心拍数を報知する。報知としては、表示するほか、記録媒体へのデータ記録や、他の装置へのデータ送信を含む。【０１６７】
図２１は、心拍波形の基本周波数の大きさによる２倍高調波の現れ方を説明している。図２１（ａ）は、心拍波形の基本周波数が高域側周波数Ｆ２以下の場合（1.5Hz）を示している。このときには、主にＦ２（２倍高調波）の曲率成分（2.68Hz）Ｅが現れる。このＥについて、スロープの傾きを判定することで心拍を判定する。【０１６８】
また、図２１（ｂ）は、心拍波形の基本周波数が高域側周波数Ｆ２と微分上限値Ｆ３との間の場合（2.68Hz）を示している。このときには、高調波がカットされるため、Ｆ２（２倍高調波）の曲率成分は零となり、基本周波数の正弦波そのものの曲率成分Ｆが現れる。このＦについて、スロープの傾きを判定することで心拍を判定する。【０１６９】
図２２は波形歪み検出の実例図を示している。ここでは、心拍の基本周波数としてはＦ３の3.67Hzが取り得るため、右肩上がりの傾きを検出するために、運動時最大心拍Ｆ３の周期の１／４以下の時間幅（68msec以下）の期間でスロープ検出を行う。また、ここでは、サンプリング値の最低値として２を設定している（図２２中のＨ）。【０１７０】
図２３は、波形歪み検出部の一構成例を説明するための図である。【０１７１】
図２３に示す構成では、入力したサンプリング値Ｓ１〜Ｓ５を順次記憶する記憶手段６０と、この記憶手段６０に記憶されるサンプリング値Ｓ１〜Ｓ５を設定値と比較する第１の比較手段６１と、第１の比較手段６１の比較結果と記憶手段６０に記憶されるサンプリング値Ｓ１〜Ｓ５を入力し、サンプリング値Ｓ１〜Ｓ５が設定値よりも大きい場合に、サンプリング値Ｓ１〜Ｓ５の大きさの順を判定する第２の比較手段６２と、第２の比較手段６２の判定結果に基づいて心拍を計数するカウンタ６３を備える。
【産業上の利用可能性】
【０１７２】
本発明の心拍波形の検出方法は、心拍計に適用する他、身体機能測定装置等の心拍を一つのデータとして身体機能を測定する装置に適応することができる。
【図１】
0004657300.tif 000002
【図２】
0004657300.tif 000003
【図３】
0004657300.tif 000004
【図４】
0004657300.tif 000005
【図５】
0004657300.tif 000006
【図６】
0004657300.tif 000007
【図７】
0004657300.tif 000008
【図８】
0004657300.tif 000009
【図９】
0004657300.tif 000010
【図１０】
0004657300.tif 000011
【図１１】
0004657300.tif 000012
【図１２】
0004657300.tif 000013
【図１３】
0004657300.tif 000014
【図１４】
0004657300.tif 000015
【図１５】
0004657300.tif 000016
【図１６】
0004657300.tif 000017
【図１７】
0004657300.tif 000018
【図１８】
0004657300.tif 000019
【図１９】
0004657300.tif 000020
【図２０】
0004657300.tif 000021
【図２１】
0004657300.tif 000022
【図２２】
0004657300.tif 000023
【図２３】
0004657300.tif 000024
【図２４】
0004657300.tif 000025
【図２５】
0004657300.tif 000026
【図２６】
0004657300.tif 000027
【図２７】
0004657300.tif 000028
【図２８】
0004657300.tif 000029
【図２９】
0004657300.tif 000030
【図３０】
0004657300.tif 000031