運動時の脈拍数，酸素飽和度測定方法及び測定装置
【課題】 被測定者の運動等による体動に起因する雑音の周期が一致或いは近接し、その振幅が大きい場合にも、正確に脈拍数及び酸素飽和度の測定が可能な脈拍数及び酸素飽和度測定方法及び測定装置を提供する。
【解決手段】 生体組織を透過又は反射した２種類の脈波信号を信号成分と雑音成分に分離するステップと、前記分離された信号成分と雑音成分の周波数スペクトルを演算するステップと、前記信号成分と雑音成分の周波数スペクトルに基づいて、体動の有無を判定するステップと、前記体動の有無の判定結果に応じて、脈拍周波数を決定するステップと前記決定された脈拍周波数から脈拍数を演算するステップと、前記決定された脈拍周波数に対応する前記２種類の脈波信号のスペクトル比から酸素飽和度を演算するステップとを含むことを特徴とする運動時の脈拍数及び酸素飽和度測定方法。
【解決手段】 生体組織を透過又は反射した２種類の脈波信号を信号成分と雑音成分に分離するステップと、前記分離された信号成分と雑音成分の周波数スペクトルを演算するステップと、前記信号成分と雑音成分の周波数スペクトルに基づいて、体動の有無を判定するステップと、前記体動の有無の判定結果に応じて、脈拍周波数を決定するステップと前記決定された脈拍周波数から脈拍数を演算するステップと、前記決定された脈拍周波数に対応する前記２種類の脈波信号のスペクトル比から酸素飽和度を演算するステップとを含むことを特徴とする運動時の脈拍数及び酸素飽和度測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、運動時に装着して脈拍数や酸素飽和度（ＳｐＯ２）の測定に適したパルスオキシメータであって、脈拍周期と運動による体動に起因する雑音の周期が一致或いは近接した場合にも適用可能な信号処理に関する。
【背景技術】
【０００２】
一つの媒体からほぼ同時に抽出された２つの信号から信号成分と雑音成分に分離する方
法には様々な方法が提案されている。
それらは、一般的には周波数領域や時間領域による処理が行われている。
医療現場でも、光電脈波計と言われる脈波波形や脈拍数を測定する装置、血液に含まれ
る吸光物質の濃度測定として、酸素飽和度ＳｐＯ2の測定装置、一酸化炭素ヘモグロビン
やＭｅｔヘモグロビン等の特殊ヘモグロビン濃度の測定装置、注入色素濃度の測定装置な
どがパルスフォトメータとして知られている。
中でも酸素飽和度ＳｐＯ2測定装置を特にパルスオキシメータと呼んでいる。
【０００３】
パルスフォトメータの原理は、対象物質への吸光性が異なる複数の波長の光を生体組織
に透過又は反射させ、その透過光又は反射光の光量を連続的に測定することで得られる脈
波データ信号から対象物質の濃度を求めるものである。
そしてその脈波データに雑音が混入すると、正しい濃度の計算が出来ず、誤処置につな
がる危険が生じる。
従来、パルスフォトメータにおいても雑音を低減するために周波数帯域を分割して信号
成分に着目し、２つの信号の相関を取るなどの方法が提案されている。
【０００４】
本出願人は、異なる２つの波長の光を生体組織に照射して透過光から得られる
２つの脈波信号のそれぞれの大きさを縦軸、横軸としてグラフを描き、その回帰直線を求
め、その回帰直線の傾きに基づいて、動脈血中の酸素飽和度ないし吸光物質濃度を求める
ことを提案している。（特許文献１）
この発明により、測定精度を高め、低消費電力化することができた。
しかし、各波長の脈波信号についての多くのサンプリングデータを用いて回帰直線ないしその傾きを求めるためには、なお多くの計算処理を要していた。
【０００５】
更に、本出願人は、周波数解析を用いてはいるが、その解析においては従来技術のように脈波信号そのものを抽出するのではなく、脈波信号の基本周波数を求め、さらには精度を高めるためにその高調波周波数を用いたフィルタを用いて脈波信号をフィルタリングする方法を提案している。（特許文献２）
また、信号分離法を使い信号とノイズを分離する方法も提案している。（特許文献３、
特許文献４）
【０００６】
しかし、いずれの方法においても、振幅比で脈波の１０倍もあるような大きな体動性ノ
イズが混入した場合、パルスレートおよび動脈血酸素飽和度を計算することは困難であり
、更なる改善が望まれていた。
【０００７】
上記改善の例として、同一の媒体からほぼ同時に抽出される２つの同種の信号を処理して共通の信号成分を抽出する計算処理負担を軽減した信号処理方法も本出願人が出願している。（特許文献５参照）
上記特許文献５では、、上記信号処理方法を適用して、前記媒体の体動によるノイズが脈波信号に生じた場合であっても、対象物質の濃度を精度よく測定することにある。
また、体動によるノイズが脈波データ信号に生じた場合であっても、脈波信号からノイ
ズを除去し、精度よく脈拍を求めている。
【特許文献１】特許第３２７０９１７号公報
【特許文献２】特開２００３−１３５４３４号公報
【特許文献３】特開２００５−９５５８１号公報
【特許文献４】特開２００５−２４５５７４号公報
【特許文献５】特開２００７−８３０２１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００８】
しかし、上記特許文献１−５に記載の技術でも、被測定者の運動等による体動に起因する雑音の周波数が脈拍の周波数に一致或いは近接し、その振幅が大きい場合には、脈拍数の正確な測定ができなかった。
【０００９】
本発明の課題は、被測定媒体からほぼ同時に抽出された２つの信号（IR,R）からそれぞれ信号成分と雑音成分に分離した後、それぞれの周波数スペクトルを求めて、求めた信号成分と雑音成分の周波数スペクトルのそれぞれの基本波及び高調波のピーク周波数の関係から、脈拍数の演算に採用する周波数を決定して脈拍数を演算する。更に、脈拍数の演算に決定された周波数に基づいて酸素飽和度（ＳｐＯ２）を演算することを特徴とする脈拍数及び酸素飽和度の測定方法及び測定装置を提供することにある。
【課題を解決するための手段】
【００１０】
前記課題を解決するために、本発明の運動時の脈拍数及び酸素飽和度測定方法は、
生体組織を透過又は反射した２種類の脈波信号を信号成分と雑音成分に分離するステップと、
前記分離された信号成分と雑音成分の周波数スペクトルを演算するステップと、
前記信号成分と雑音成分の周波数スペクトルに基づいて、
体動の有無を判定するステップと、
前記体動の有無の判定結果に応じて、脈拍周波数を決定するステップと、
前記決定された脈拍周波数から脈拍数を演算するステップと、
前記決定された脈拍周波数に対応する前記２種類の脈波信号のスペクトル比から酸素飽和度を演算するステップとを含むことを特徴とする。（請求項１）
【００１１】
また、前記脈拍周波数の決定のステップは、体動が無い場合には、前記信号成分の基本周波数を前記脈拍周波数として決定するステップを含むことを特徴とする。（請求項２）
【００１２】
また、前記脈拍周波数の決定のステップは、体動が有る場合に、信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以上の場合に、前記信号成分周波数スペクトルの最大振幅ピーク周波数を前記脈拍周波数として決定するステップを含むことを特徴とする。（請求項３）
【００１３】
また、前記脈拍周波数の決定のステップは、更に、前記信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以下で、前記信号成分周波数スペクトルの最大振幅ピーク周波数と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以下の場合に、前記信号成分周波数スペクトルの最大振幅ピーク周波数を前記脈拍周波数として決定するステップを含むことを特徴とする。（請求項４）
【００１４】
また、前記脈拍周波数の決定のステップは、更に、前記信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以下で、前記信号成分周波数スペクトルの最大振幅ピーク周波数と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以上で、信号成分周波数スペクトルの２番目の振幅のピーク周波数 と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以下の場合に、前記信号成分周波数スペクトルの２番目の振幅のピーク周波数を前記脈拍周波数として決定するステップを含むことを特徴とする。（請求項５）
【００１５】
また、前記脈拍周波数の決定のステップは、更に、前記信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以下で、前記信号成分周波数スペクトルの最大振幅ピーク周波数と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以上で、前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２と前時間窓脈拍周波数との差と、前記前時間窓脈拍周波数との比が所定％以下の場合には、前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２を脈拍周波数として決定するステップを含むことを特徴とする。（請求項６）
【００１６】
また、本発明の運動時の脈拍数及び酸素飽和度測定装置は、生体組織を透過又は反射した２種類の脈波信号を信号成分と雑音成分に分離する手段と、
前記分離された信号成分と雑音成分の周波数スペクトルを演算する手段と、
前記信号成分と雑音成分の周波数スペクトルに基づいて、体動の有無を判定する手段と、
前記体動の有無の判断結果に応じて、前記信号成分と雑音成分の周波数スペクトルに基づいて、脈拍周波数を決定する手段と、
前記決定された脈拍周波数から脈拍数を演算する手段と、
前記決定された脈拍周波数に対応する前記２種類の脈波信号のスペクトル比から酸素飽和度を演算する手段とを備えることを特徴とする。（請求項７）
【発明の効果】
【００１７】
請求項１〜７に記載の本発明の脈拍数及び酸素飽和度の測定方法及び測定装置によれば、運動時の体動によるノイズ（雑音）の影響を排除して正確な脈拍数の測定が可能になる。
また、運動時に体動による動脈血以外の組織や静脈血の動きがある場合には、その動きによる周波数帯域の信号を用いないので、正確な酸素飽和度の測定が可能になる。
また、体動ない安静時には、高い時間分解能の脈拍計測が可能となる従来のゼロクロス法等の測定方法（装置）を併用して、選択して使用することができる。その場合には、脈拍数が低い安静時の計測に十分な時間分解能が得られる。また、自律神経の検査のための脈拍数のゆらぎ計測に必要な時間分解能を得られる。
【発明を実施するための最良の形態】
【００１８】
本発明の脈拍数及び酸素飽和度の測定方法及び測定装置におけるデータ処理の全体的な構成を図１の処理フローを用いて説明する。
・同一の被測定媒体からほぼ同時に抽出される２つの同種の信号（IR,R）を、信号成分と雑音成分に分離する。（ステップS1）
・ステップS1で分離した信号成分と、雑音成分に対してFFT処理等により、それぞれの周波数スペクトルを算出する。（ステップS2）
・体動の有無の判定及び後述の判定手法によって脈拍周波数算出する。（ステップS3）
・ステップS3で算出した脈拍周波数に基づいて、脈拍数を算出する。（ステップS4）
・ステップS3で算出して脈拍周波数におけるRとIRの周波数スペクトルを算出する。（ステップS5）
・ステップS5で算出した脈拍周波数におけるRとIRの周波数スペクトル比を抽出する。（ステップS6）
・ステップS6のRとIRの周波数スペクトル比から酸素飽和度（SpO2）を算出する。（ステップS7）
【００１９】
本発明の脈拍周波数の決定判断は、以下の性質に基づいている。
被測定媒体から測定された脈拍周波数スペクトルには、基本波の他に第２高調波、第３高調波等が含まれ、それぞれ基本波周波数の整数倍のところに現れて、通常、その振幅はほぼ基本波、第２高調波、第３高調波の順に小さくなる性質がある。
また、体動が無い時には、動脈血が脈動して厚みが変化し、体動がある時は動脈血と静脈血の厚みが同時に変化する。
【００２０】
本発明は、以下の３つの観点から解析（分類）されている。
(a) 体動で主として動いた部分は、動脈血か？ 静脈血か？
信号成分は動脈血により生成されるので、体動で主として動いた成分が動脈血の場合は、その影響は信号成分に現れる。
また、体動で主として動いた成分が静脈血の場合は、その影響は雑音成分に現れる。
(b) 体動の周波数は脈拍の周波数と同じか？ 隣り合っているか？ 異なっているか？
体動と脈拍の周波数が同じか、隣り合っている場合には、信号成分のスペクトルと体動によるスペクトルとは分離できず、単峰性となる。
(c) 体動の強度は強いか？弱いか？
強い体動の時には、動脈血、静脈血の動きが大きくなるので、体動によるスペクトルの振幅は大きくなる。
弱い体動の時には、体動のスペクトルの振幅は小さくなる。
【００２１】
以下に、被測定媒体からほぼ同時に抽出された２つの信号（IR,R）からそれぞれ信号成分と雑音成分に分離した後、上記の３つの観点から、それぞれの周波数スペクトルを求めた例で、体動があった場合の種々の例について図２の表１を用いて説明する。
【００２２】
図２の表１は、体動があるときの信号成分周波数スペクトルと雑音成分周波数スペクトルとの関係を示すものである。
図２の表１では、脈拍によるスペクトルを実線で示し体動によるスペクトルを破線で示している。
【００２３】
また、脈拍によるスペクトルと体動によるスペクトルが同じ周波数の場合は、信号成分と雑音成分が分離できず、両者が合計されたスペクトルとして表示されている。（(1)，(4)，（7），(10)）
【００２４】
また、脈拍によるスペクトルと体動によるスペクトルが極めて近い隣り合った周波数S1とSn周波数の場合は、完全には分離できずFFT処理しても単峰となって、区別ができず、S1単独では認識できない。（(2)，(5)，（8），(11)）
【００２５】
また、体動の強弱は動脈血及び静脈血共に動きが大きくなるので、強い体動作の時は体動によるスペクトル振幅は大きくなり、体動が弱い体動の時は体動のスペクトルは小さくなる。
【００２６】
以下に、図２の表１における１２の例について、信号成分周波数スペクトルのS1が、脈拍の周波数と判断できる理由を説明する。
【００２７】
(6)，(9)，(12)，(5)，(8)，(11)では、信号成分周波数スペクトルの最大振幅ピーク周波数S1が雑音成分周波数スペクトルの最大振幅ピーク周波数Nと一致していないので、S1が脈拍の周波数と判断できる。
【００２８】
(1)，(4)，(7)，(10)では、信号成分周波数スペクトルの最大振幅ピーク周波数S1が雑音成分周波数スペクトルの最大振幅ピーク周波数Nと一致するが、S1の２倍の周波数ピークS2があるので、S1が脈拍の周波数と判断できる。
【００２９】
(3)では、信号成分周波数スペクトルの最大振幅ピーク周波数Snが雑音成分周波数スペクトルの最大振幅ピーク周波数Nと一致するが、Snの２倍の周波数ピークがないので、
Snは脈拍の周波数ではなく、体動の周波数であることが判断できる。
信号成分周波数スペクトルの第２の振幅ピーク周波数S1には２倍の周波数ピークS2があるのでS1が脈拍の周波数と判断できる。
【００３０】
(2)では、信号成分周波数スペクトルの最大振幅ピーク周波数Snが雑音成分周波数スペクトルの最大振幅ピーク周波数Nと一致するが、Snの２倍の周波数ピークがないので、
Snは脈拍の周波数ではなく、体動の周波数であることが判断できる。
また、信号成分周波数スペクトルの脈拍によるピーク周波数S1はSnと区別がつかず、S1を見つけることはできない。
そこで、信号成分周波数スペクトルの第２の振幅ピーク周波数S2の1/2倍の周波数を先行する時間窓（前回の判断）で決定した脈拍の周波数と比較して、その変化率が一定の範囲内であれば、S2の1/2倍の周波数を脈拍の周波数として決定できる。
ここで、一定の範囲とは、時間窓を切り替える時間内に脈拍数が生理的に最大変化し得る変化率の大きさ（例えば、10％）である。
【００３１】
上記判断例を表に纏めたものが、図３の表２の処理方法に示す。
図３の表２の優先順位１，２の例は、体動が無い場合であるので、前記図２の表１の(1)〜(12)には含まれてはいない。
【００３２】
本発明の脈拍数（PR）の詳細な演算手順を図４のフローチャートを用いて説明する。
この演算手順を「運動アルゴリズムで決定された脈拍の周波数決定方法」と称する。
このフローチャートにおける判断手法は、脈拍数及び酸素飽和度の測定装置に搭載されるコンピュータによる自動判断処理に適している。
【００３３】
・被測定媒体からほぼ同時に抽出された２つの信号（IR,R）を信号成分と雑音成分に分離して、FFT処理等によって、それぞれの周波数スペクトルを求める。（ステップS11）
・fN（体動による雑音成分スペクトル）があるか否かの判断をする。（ステップS12）
・ステップS12の判断がNoの場合には、信号成分スペクトルの周波数であるfS1，fS2，fS3が整数倍で、そのピークがfS1＞fS2＞fS3で、その差が±２０％以内であるか否かの判断をする。（ステップS13）
【００３４】
このステップS13の判断は、図３の表２の優先順位１のFFTパターンに相当する場合であって、基本波(fS1)と第２高調波（fS2）及び第３高調波(fs3)のピークがfS1＞fS2＞fS3の条件を満たしている。
【００３５】
・ステップS13の判断がYesの場合には、脈拍周波数としてfS1を採用して、脈拍数（PR）をPR=fS1×６０で算出する。（ステップS19）
【００３６】
・ステップS13の判断がNoの場合には、信号成分スペクトルの周波数であるfS1，fS2，fS3が整数倍で、そのピークがfS1＞fS3＞fS2で、その差が±２０％以内であるか否かの判断をする。（ステップS14）
【００３７】
このステップS14の判断は、図３の表２の優先順位２のFFTパターンに相当する場合であって、基本波(fS1)と第２高調波（fS3）及び第３高調波(fs2)のピークがfS1＞fS2＞fS3の条件を満たしている。
【００３８】
・ステップS14の判断がYesの場合には、脈拍周波数としてfS1を採用して、脈拍数（PR）をPR=fS1×６０で算出する。（ステップS19）
【００３９】
・ステップS14の判断がNo又はステップS12の判断がYesの場合には、
|fS1-fN|＞Δf1(Hz)であるか否かの判断をする。（ステップS15）
ここで、fS1は、信号成分周波数スペクトルの最大振幅ピーク周波数、fNは、雑音成分周波数スペクトルの最大振幅ピーク周波数、Δf1の値は、例えば、０．２５（Hz)である。
【００４０】
このステップS15の判断は、図３の表２の優先順位３のFFTパターンに相当する場合であって、図２の表１の（５）、（６）、（８）、（９）、（１１）、（１２）、における基本波の信号成分S1の周波数(fS1)と雑音成分Nの周波数（fN）との差が、所定値（０．２５Hz）以上の場合である。
【００４１】
・ステップS15の判断がYesの場合には、脈拍周波数としてfS1を採用して、脈拍数（PR）をPR=fS1×６０で算出する。（ステップS19）
【００４２】
・ステップS15の判断がNoの場合には、
|fS1-fN|＜Δf1(Hz)及び |fS1-fS2/2| Δf2(Hz)であるか否かの判断をする。（ステップS16）
ここで、fS1は、信号成分周波数スペクトルの最大振幅ピーク周波数、fNは、雑音成分周波数スペクトルの最大振幅ピーク周波数、fS2は、信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数、Δf1及びΔf2の値は、例えば、共に０．２５（Hz)である。
【００４３】
このステップS16の判断は、図３の表２の優先順位４のFFTパターンに相当する場合であって、図２の表１の（１）、（４）、（７）、（１０）、における基本波の信号成分(S1)の周波数(fS1)と雑音成分Nの周波数(fN)との差が、所定値（０．２５Hz）以下の場合であり、且つ、基本波の信号成分S1の周波数(fS1)と第２高調波S2の周波数(fS2)の1/2の周波数との差が、所定値（０．２５Hz）以下の場合である。
【００４４】
・ステップS16の判断がYesの場合には、脈拍周波数としてfS1を採用して、脈拍数（PR）をPR=fS1×６０で算出する。（ステップS19）
【００４５】
・ステップS16の判断がNoの場合には、
|fSn-fN|＜Δf1(Hz)及び |fSn-fS2/2|＞Δf2(Hz) AND
|fSn-fS3/2|＞Δf2(Hz)及び|fS1-fS2/2|＜Δf2(Hz)か否かの判断をする。（ステップS17）
ここで、fSnは、信号成分周波数スペクトルの最大振幅ピーク周波数、fN1は、雑音成分周波数スペクトルの最大振幅ピーク周波数、fS2及びfS3は、信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数、fS1は、信号成分周波数スペクトルの２番目の振幅のピーク周波数、Δf1及びΔf2の値は、例えば、共に０．２５（Hz)である。
【００４６】
このステップS17の判断は、図３の表２の優先順位５のFFTパターンに相当する場合であって、図２の表１の（３）における信号成分周波数スペクトルに現れた雑音（体動）Snの周波数(fSn)と雑音成分Nの周波数(fN)との差が、所定値（０．２５Hz）以下の場合であり、信号成分の周波数スペクトルに現れた雑音（体動）Snの周波数(fSn)と第２高調波S2の周波数(fS2)の1/2の周波数との差が、所定値（０．２５Hz）以上の場合である。
また、信号成分の周波数スペクトルに現れた雑音（体動）Snの周波数（fSn）と第３高調波S3の周波数(fS3)の1/2の周波数との差が、所定値（０．２５Hz）以上の場合であり、信号成分の基本波S1の周波数(fS1)と第２高調波S2の周波数(fS2)の1/2の周波数との差が、所定値（０．２５Hz）以下の場合である。
【００４７】
・ステップS17の判断がYesの場合には、脈拍周波数としてfS1を採用して、脈拍数（PR）をPR=fS1×６０で算出する。（ステップS19）
【００４８】
・ステップS16の判断がNoの場合には、
|fSn-fN|＜Δf1(Hz)及び |fSn-fS2/2|＞Δf2(Hz) AND
|fSn-fS3/2|＜Δf2(Hz)か否かの判断及び |fS2/2-前時間窓脈拍f|／前時間窓脈拍f＜ｒ（％）否かの判断をする。（ステップS18）
ここで、fSnは、信号成分周波数スペクトルの最大振幅ピーク周波数、fNは、雑音成分周波数スペクトルの最大振幅ピーク周波数、fS2及びfS3は、信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数、Δf1及びΔf2の値は、例えば、共に０．２５（Hz)であり、ｒの値は１０（％）である。
【００４９】
このステップS18の判断は、図３の表２の優先順位６のFFTパターンに相当する場合であって、図２の表１の（２）における信号成分周波数スペクトルに現れた雑音（体動）Snの周波数(fSn)と雑音成分Nの周波数(fN)との差が、所定値（０．２５Hz）以下の場合であり、信号成分の周波数スペクトルに現れた雑音（体動）Snの周波数(fSn)と第２高調波S2の周波数(fS2)の1/2の周波数との差が、所定値（０．２５Hz）以上の場合である。
また、信号成分の周波数スペクトルに現れた雑音（体動）Snの周波数(fSn)と第３高調波S3の周波数(fS3)の1/2の周波数との差が、所定値（０．２５Hz）以上の場合であり、第２高調波S2の周波数(fS2)の1/2と前時間窓脈拍f との差と前時間窓脈拍ｆとの比が所定％（１０％）以下の場合である。
【００５０】
・ステップS18の判断がYesの場合には、脈拍周波数としてfS2/2を採用して、脈拍数（PR）をPR=fS2/2×６０で算出する。（ステップS20）
・ステップS18の判断がNoの場合には、脈拍周波数を決定せずに、ステップS12に戻る。
【００５１】
図４のフローチャートでは、脈拍数を算出するまでが開示されているが、酸素飽和度（SpO2）は、脈拍数を算出するのに使用した脈拍周波数における被測定媒体からほぼ同時に抽出された２つの信号（IR,R）の比を抽出することによって得られる。
【００５２】
本発明の「運動アルゴリズムで決定された脈拍の周波数決定方法」と、従来の「ゼロクロス法」脈拍数の算出とを併用した例を図５に示す。
図５では、「運動アルゴリズムで決定された脈拍の周波数決定方法」による測定装置１と「ゼロクロス法」による測定装置２とを併用して、ＰＲの選択手段３で運動時には、１を選択し、安静時には２を選択する。
【００５３】
この選択によって、運動時に体動による動脈血以外の組織や静脈血の動きがある場合には、その動きによる周波数帯域の信号を用いないので、正確な酸素飽和度の測定が可能になる。
また、体動ない安静時には、高い時間分解能の脈拍計測が可能となる従来のゼロクロス法等の測定方法（装置）を併用して、選択して使用することができる。その場合には、脈拍数が低い安静時の計測に十分な時間分解能が得られる。また、自律神経の検査のための脈拍数のゆらぎ計測に必要な時間分解能を得られる。
【図面の簡単な説明】
【００５４】
【図１】本発明の脈拍数及び酸素飽和度の測定方法及び測定装置におけるデータ処理の全体的な構成を示すフローチャートである。
【図２】被測定媒体からほぼ同時に抽出された２つの信号（IR,R）からそれぞれ信号成分と雑音成分に分離した後、それぞれの周波数スペクトルを求めた例で、体動があった場合の種々の例を示す図である。
【図３】図２の表に示した例における脈拍数を算出するための判断例を表に纏めたものである。
【図４】本発明の脈拍数（PR）の詳細な算出手順を示すフローチャートである。
【図５】「運動アルゴリズムで決定された脈拍の周波数決定方法」と「ゼロクロス法」による脈拍数の算出を併用した例を示す図である。
【符号の説明】
【００５５】
１：「運動アルゴリズムで決定された脈拍の周波数決定方法」による脈拍数及び酸素飽和度測定装置
２：「ゼロクロス法」による脈拍数及び酸素飽和度測定装置
３：選択手段
【技術分野】
【０００１】
本発明は、運動時に装着して脈拍数や酸素飽和度（ＳｐＯ２）の測定に適したパルスオキシメータであって、脈拍周期と運動による体動に起因する雑音の周期が一致或いは近接した場合にも適用可能な信号処理に関する。
【背景技術】
【０００２】
一つの媒体からほぼ同時に抽出された２つの信号から信号成分と雑音成分に分離する方
法には様々な方法が提案されている。
それらは、一般的には周波数領域や時間領域による処理が行われている。
医療現場でも、光電脈波計と言われる脈波波形や脈拍数を測定する装置、血液に含まれ
る吸光物質の濃度測定として、酸素飽和度ＳｐＯ2の測定装置、一酸化炭素ヘモグロビン
やＭｅｔヘモグロビン等の特殊ヘモグロビン濃度の測定装置、注入色素濃度の測定装置な
どがパルスフォトメータとして知られている。
中でも酸素飽和度ＳｐＯ2測定装置を特にパルスオキシメータと呼んでいる。
【０００３】
パルスフォトメータの原理は、対象物質への吸光性が異なる複数の波長の光を生体組織
に透過又は反射させ、その透過光又は反射光の光量を連続的に測定することで得られる脈
波データ信号から対象物質の濃度を求めるものである。
そしてその脈波データに雑音が混入すると、正しい濃度の計算が出来ず、誤処置につな
がる危険が生じる。
従来、パルスフォトメータにおいても雑音を低減するために周波数帯域を分割して信号
成分に着目し、２つの信号の相関を取るなどの方法が提案されている。
【０００４】
本出願人は、異なる２つの波長の光を生体組織に照射して透過光から得られる
２つの脈波信号のそれぞれの大きさを縦軸、横軸としてグラフを描き、その回帰直線を求
め、その回帰直線の傾きに基づいて、動脈血中の酸素飽和度ないし吸光物質濃度を求める
ことを提案している。（特許文献１）
この発明により、測定精度を高め、低消費電力化することができた。
しかし、各波長の脈波信号についての多くのサンプリングデータを用いて回帰直線ないしその傾きを求めるためには、なお多くの計算処理を要していた。
【０００５】
更に、本出願人は、周波数解析を用いてはいるが、その解析においては従来技術のように脈波信号そのものを抽出するのではなく、脈波信号の基本周波数を求め、さらには精度を高めるためにその高調波周波数を用いたフィルタを用いて脈波信号をフィルタリングする方法を提案している。（特許文献２）
また、信号分離法を使い信号とノイズを分離する方法も提案している。（特許文献３、
特許文献４）
【０００６】
しかし、いずれの方法においても、振幅比で脈波の１０倍もあるような大きな体動性ノ
イズが混入した場合、パルスレートおよび動脈血酸素飽和度を計算することは困難であり
、更なる改善が望まれていた。
【０００７】
上記改善の例として、同一の媒体からほぼ同時に抽出される２つの同種の信号を処理して共通の信号成分を抽出する計算処理負担を軽減した信号処理方法も本出願人が出願している。（特許文献５参照）
上記特許文献５では、、上記信号処理方法を適用して、前記媒体の体動によるノイズが脈波信号に生じた場合であっても、対象物質の濃度を精度よく測定することにある。
また、体動によるノイズが脈波データ信号に生じた場合であっても、脈波信号からノイ
ズを除去し、精度よく脈拍を求めている。
【特許文献１】特許第３２７０９１７号公報
【特許文献２】特開２００３−１３５４３４号公報
【特許文献３】特開２００５−９５５８１号公報
【特許文献４】特開２００５−２４５５７４号公報
【特許文献５】特開２００７−８３０２１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００８】
しかし、上記特許文献１−５に記載の技術でも、被測定者の運動等による体動に起因する雑音の周波数が脈拍の周波数に一致或いは近接し、その振幅が大きい場合には、脈拍数の正確な測定ができなかった。
【０００９】
本発明の課題は、被測定媒体からほぼ同時に抽出された２つの信号（IR,R）からそれぞれ信号成分と雑音成分に分離した後、それぞれの周波数スペクトルを求めて、求めた信号成分と雑音成分の周波数スペクトルのそれぞれの基本波及び高調波のピーク周波数の関係から、脈拍数の演算に採用する周波数を決定して脈拍数を演算する。更に、脈拍数の演算に決定された周波数に基づいて酸素飽和度（ＳｐＯ２）を演算することを特徴とする脈拍数及び酸素飽和度の測定方法及び測定装置を提供することにある。
【課題を解決するための手段】
【００１０】
前記課題を解決するために、本発明の運動時の脈拍数及び酸素飽和度測定方法は、
生体組織を透過又は反射した２種類の脈波信号を信号成分と雑音成分に分離するステップと、
前記分離された信号成分と雑音成分の周波数スペクトルを演算するステップと、
前記信号成分と雑音成分の周波数スペクトルに基づいて、
体動の有無を判定するステップと、
前記体動の有無の判定結果に応じて、脈拍周波数を決定するステップと、
前記決定された脈拍周波数から脈拍数を演算するステップと、
前記決定された脈拍周波数に対応する前記２種類の脈波信号のスペクトル比から酸素飽和度を演算するステップとを含むことを特徴とする。（請求項１）
【００１１】
また、前記脈拍周波数の決定のステップは、体動が無い場合には、前記信号成分の基本周波数を前記脈拍周波数として決定するステップを含むことを特徴とする。（請求項２）
【００１２】
また、前記脈拍周波数の決定のステップは、体動が有る場合に、信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以上の場合に、前記信号成分周波数スペクトルの最大振幅ピーク周波数を前記脈拍周波数として決定するステップを含むことを特徴とする。（請求項３）
【００１３】
また、前記脈拍周波数の決定のステップは、更に、前記信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以下で、前記信号成分周波数スペクトルの最大振幅ピーク周波数と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以下の場合に、前記信号成分周波数スペクトルの最大振幅ピーク周波数を前記脈拍周波数として決定するステップを含むことを特徴とする。（請求項４）
【００１４】
また、前記脈拍周波数の決定のステップは、更に、前記信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以下で、前記信号成分周波数スペクトルの最大振幅ピーク周波数と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以上で、信号成分周波数スペクトルの２番目の振幅のピーク周波数 と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以下の場合に、前記信号成分周波数スペクトルの２番目の振幅のピーク周波数を前記脈拍周波数として決定するステップを含むことを特徴とする。（請求項５）
【００１５】
また、前記脈拍周波数の決定のステップは、更に、前記信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以下で、前記信号成分周波数スペクトルの最大振幅ピーク周波数と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以上で、前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２と前時間窓脈拍周波数との差と、前記前時間窓脈拍周波数との比が所定％以下の場合には、前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２を脈拍周波数として決定するステップを含むことを特徴とする。（請求項６）
【００１６】
また、本発明の運動時の脈拍数及び酸素飽和度測定装置は、生体組織を透過又は反射した２種類の脈波信号を信号成分と雑音成分に分離する手段と、
前記分離された信号成分と雑音成分の周波数スペクトルを演算する手段と、
前記信号成分と雑音成分の周波数スペクトルに基づいて、体動の有無を判定する手段と、
前記体動の有無の判断結果に応じて、前記信号成分と雑音成分の周波数スペクトルに基づいて、脈拍周波数を決定する手段と、
前記決定された脈拍周波数から脈拍数を演算する手段と、
前記決定された脈拍周波数に対応する前記２種類の脈波信号のスペクトル比から酸素飽和度を演算する手段とを備えることを特徴とする。（請求項７）
【発明の効果】
【００１７】
請求項１〜７に記載の本発明の脈拍数及び酸素飽和度の測定方法及び測定装置によれば、運動時の体動によるノイズ（雑音）の影響を排除して正確な脈拍数の測定が可能になる。
また、運動時に体動による動脈血以外の組織や静脈血の動きがある場合には、その動きによる周波数帯域の信号を用いないので、正確な酸素飽和度の測定が可能になる。
また、体動ない安静時には、高い時間分解能の脈拍計測が可能となる従来のゼロクロス法等の測定方法（装置）を併用して、選択して使用することができる。その場合には、脈拍数が低い安静時の計測に十分な時間分解能が得られる。また、自律神経の検査のための脈拍数のゆらぎ計測に必要な時間分解能を得られる。
【発明を実施するための最良の形態】
【００１８】
本発明の脈拍数及び酸素飽和度の測定方法及び測定装置におけるデータ処理の全体的な構成を図１の処理フローを用いて説明する。
・同一の被測定媒体からほぼ同時に抽出される２つの同種の信号（IR,R）を、信号成分と雑音成分に分離する。（ステップS1）
・ステップS1で分離した信号成分と、雑音成分に対してFFT処理等により、それぞれの周波数スペクトルを算出する。（ステップS2）
・体動の有無の判定及び後述の判定手法によって脈拍周波数算出する。（ステップS3）
・ステップS3で算出した脈拍周波数に基づいて、脈拍数を算出する。（ステップS4）
・ステップS3で算出して脈拍周波数におけるRとIRの周波数スペクトルを算出する。（ステップS5）
・ステップS5で算出した脈拍周波数におけるRとIRの周波数スペクトル比を抽出する。（ステップS6）
・ステップS6のRとIRの周波数スペクトル比から酸素飽和度（SpO2）を算出する。（ステップS7）
【００１９】
本発明の脈拍周波数の決定判断は、以下の性質に基づいている。
被測定媒体から測定された脈拍周波数スペクトルには、基本波の他に第２高調波、第３高調波等が含まれ、それぞれ基本波周波数の整数倍のところに現れて、通常、その振幅はほぼ基本波、第２高調波、第３高調波の順に小さくなる性質がある。
また、体動が無い時には、動脈血が脈動して厚みが変化し、体動がある時は動脈血と静脈血の厚みが同時に変化する。
【００２０】
本発明は、以下の３つの観点から解析（分類）されている。
(a) 体動で主として動いた部分は、動脈血か？ 静脈血か？
信号成分は動脈血により生成されるので、体動で主として動いた成分が動脈血の場合は、その影響は信号成分に現れる。
また、体動で主として動いた成分が静脈血の場合は、その影響は雑音成分に現れる。
(b) 体動の周波数は脈拍の周波数と同じか？ 隣り合っているか？ 異なっているか？
体動と脈拍の周波数が同じか、隣り合っている場合には、信号成分のスペクトルと体動によるスペクトルとは分離できず、単峰性となる。
(c) 体動の強度は強いか？弱いか？
強い体動の時には、動脈血、静脈血の動きが大きくなるので、体動によるスペクトルの振幅は大きくなる。
弱い体動の時には、体動のスペクトルの振幅は小さくなる。
【００２１】
以下に、被測定媒体からほぼ同時に抽出された２つの信号（IR,R）からそれぞれ信号成分と雑音成分に分離した後、上記の３つの観点から、それぞれの周波数スペクトルを求めた例で、体動があった場合の種々の例について図２の表１を用いて説明する。
【００２２】
図２の表１は、体動があるときの信号成分周波数スペクトルと雑音成分周波数スペクトルとの関係を示すものである。
図２の表１では、脈拍によるスペクトルを実線で示し体動によるスペクトルを破線で示している。
【００２３】
また、脈拍によるスペクトルと体動によるスペクトルが同じ周波数の場合は、信号成分と雑音成分が分離できず、両者が合計されたスペクトルとして表示されている。（(1)，(4)，（7），(10)）
【００２４】
また、脈拍によるスペクトルと体動によるスペクトルが極めて近い隣り合った周波数S1とSn周波数の場合は、完全には分離できずFFT処理しても単峰となって、区別ができず、S1単独では認識できない。（(2)，(5)，（8），(11)）
【００２５】
また、体動の強弱は動脈血及び静脈血共に動きが大きくなるので、強い体動作の時は体動によるスペクトル振幅は大きくなり、体動が弱い体動の時は体動のスペクトルは小さくなる。
【００２６】
以下に、図２の表１における１２の例について、信号成分周波数スペクトルのS1が、脈拍の周波数と判断できる理由を説明する。
【００２７】
(6)，(9)，(12)，(5)，(8)，(11)では、信号成分周波数スペクトルの最大振幅ピーク周波数S1が雑音成分周波数スペクトルの最大振幅ピーク周波数Nと一致していないので、S1が脈拍の周波数と判断できる。
【００２８】
(1)，(4)，(7)，(10)では、信号成分周波数スペクトルの最大振幅ピーク周波数S1が雑音成分周波数スペクトルの最大振幅ピーク周波数Nと一致するが、S1の２倍の周波数ピークS2があるので、S1が脈拍の周波数と判断できる。
【００２９】
(3)では、信号成分周波数スペクトルの最大振幅ピーク周波数Snが雑音成分周波数スペクトルの最大振幅ピーク周波数Nと一致するが、Snの２倍の周波数ピークがないので、
Snは脈拍の周波数ではなく、体動の周波数であることが判断できる。
信号成分周波数スペクトルの第２の振幅ピーク周波数S1には２倍の周波数ピークS2があるのでS1が脈拍の周波数と判断できる。
【００３０】
(2)では、信号成分周波数スペクトルの最大振幅ピーク周波数Snが雑音成分周波数スペクトルの最大振幅ピーク周波数Nと一致するが、Snの２倍の周波数ピークがないので、
Snは脈拍の周波数ではなく、体動の周波数であることが判断できる。
また、信号成分周波数スペクトルの脈拍によるピーク周波数S1はSnと区別がつかず、S1を見つけることはできない。
そこで、信号成分周波数スペクトルの第２の振幅ピーク周波数S2の1/2倍の周波数を先行する時間窓（前回の判断）で決定した脈拍の周波数と比較して、その変化率が一定の範囲内であれば、S2の1/2倍の周波数を脈拍の周波数として決定できる。
ここで、一定の範囲とは、時間窓を切り替える時間内に脈拍数が生理的に最大変化し得る変化率の大きさ（例えば、10％）である。
【００３１】
上記判断例を表に纏めたものが、図３の表２の処理方法に示す。
図３の表２の優先順位１，２の例は、体動が無い場合であるので、前記図２の表１の(1)〜(12)には含まれてはいない。
【００３２】
本発明の脈拍数（PR）の詳細な演算手順を図４のフローチャートを用いて説明する。
この演算手順を「運動アルゴリズムで決定された脈拍の周波数決定方法」と称する。
このフローチャートにおける判断手法は、脈拍数及び酸素飽和度の測定装置に搭載されるコンピュータによる自動判断処理に適している。
【００３３】
・被測定媒体からほぼ同時に抽出された２つの信号（IR,R）を信号成分と雑音成分に分離して、FFT処理等によって、それぞれの周波数スペクトルを求める。（ステップS11）
・fN（体動による雑音成分スペクトル）があるか否かの判断をする。（ステップS12）
・ステップS12の判断がNoの場合には、信号成分スペクトルの周波数であるfS1，fS2，fS3が整数倍で、そのピークがfS1＞fS2＞fS3で、その差が±２０％以内であるか否かの判断をする。（ステップS13）
【００３４】
このステップS13の判断は、図３の表２の優先順位１のFFTパターンに相当する場合であって、基本波(fS1)と第２高調波（fS2）及び第３高調波(fs3)のピークがfS1＞fS2＞fS3の条件を満たしている。
【００３５】
・ステップS13の判断がYesの場合には、脈拍周波数としてfS1を採用して、脈拍数（PR）をPR=fS1×６０で算出する。（ステップS19）
【００３６】
・ステップS13の判断がNoの場合には、信号成分スペクトルの周波数であるfS1，fS2，fS3が整数倍で、そのピークがfS1＞fS3＞fS2で、その差が±２０％以内であるか否かの判断をする。（ステップS14）
【００３７】
このステップS14の判断は、図３の表２の優先順位２のFFTパターンに相当する場合であって、基本波(fS1)と第２高調波（fS3）及び第３高調波(fs2)のピークがfS1＞fS2＞fS3の条件を満たしている。
【００３８】
・ステップS14の判断がYesの場合には、脈拍周波数としてfS1を採用して、脈拍数（PR）をPR=fS1×６０で算出する。（ステップS19）
【００３９】
・ステップS14の判断がNo又はステップS12の判断がYesの場合には、
|fS1-fN|＞Δf1(Hz)であるか否かの判断をする。（ステップS15）
ここで、fS1は、信号成分周波数スペクトルの最大振幅ピーク周波数、fNは、雑音成分周波数スペクトルの最大振幅ピーク周波数、Δf1の値は、例えば、０．２５（Hz)である。
【００４０】
このステップS15の判断は、図３の表２の優先順位３のFFTパターンに相当する場合であって、図２の表１の（５）、（６）、（８）、（９）、（１１）、（１２）、における基本波の信号成分S1の周波数(fS1)と雑音成分Nの周波数（fN）との差が、所定値（０．２５Hz）以上の場合である。
【００４１】
・ステップS15の判断がYesの場合には、脈拍周波数としてfS1を採用して、脈拍数（PR）をPR=fS1×６０で算出する。（ステップS19）
【００４２】
・ステップS15の判断がNoの場合には、
|fS1-fN|＜Δf1(Hz)及び |fS1-fS2/2| Δf2(Hz)であるか否かの判断をする。（ステップS16）
ここで、fS1は、信号成分周波数スペクトルの最大振幅ピーク周波数、fNは、雑音成分周波数スペクトルの最大振幅ピーク周波数、fS2は、信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数、Δf1及びΔf2の値は、例えば、共に０．２５（Hz)である。
【００４３】
このステップS16の判断は、図３の表２の優先順位４のFFTパターンに相当する場合であって、図２の表１の（１）、（４）、（７）、（１０）、における基本波の信号成分(S1)の周波数(fS1)と雑音成分Nの周波数(fN)との差が、所定値（０．２５Hz）以下の場合であり、且つ、基本波の信号成分S1の周波数(fS1)と第２高調波S2の周波数(fS2)の1/2の周波数との差が、所定値（０．２５Hz）以下の場合である。
【００４４】
・ステップS16の判断がYesの場合には、脈拍周波数としてfS1を採用して、脈拍数（PR）をPR=fS1×６０で算出する。（ステップS19）
【００４５】
・ステップS16の判断がNoの場合には、
|fSn-fN|＜Δf1(Hz)及び |fSn-fS2/2|＞Δf2(Hz) AND
|fSn-fS3/2|＞Δf2(Hz)及び|fS1-fS2/2|＜Δf2(Hz)か否かの判断をする。（ステップS17）
ここで、fSnは、信号成分周波数スペクトルの最大振幅ピーク周波数、fN1は、雑音成分周波数スペクトルの最大振幅ピーク周波数、fS2及びfS3は、信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数、fS1は、信号成分周波数スペクトルの２番目の振幅のピーク周波数、Δf1及びΔf2の値は、例えば、共に０．２５（Hz)である。
【００４６】
このステップS17の判断は、図３の表２の優先順位５のFFTパターンに相当する場合であって、図２の表１の（３）における信号成分周波数スペクトルに現れた雑音（体動）Snの周波数(fSn)と雑音成分Nの周波数(fN)との差が、所定値（０．２５Hz）以下の場合であり、信号成分の周波数スペクトルに現れた雑音（体動）Snの周波数(fSn)と第２高調波S2の周波数(fS2)の1/2の周波数との差が、所定値（０．２５Hz）以上の場合である。
また、信号成分の周波数スペクトルに現れた雑音（体動）Snの周波数（fSn）と第３高調波S3の周波数(fS3)の1/2の周波数との差が、所定値（０．２５Hz）以上の場合であり、信号成分の基本波S1の周波数(fS1)と第２高調波S2の周波数(fS2)の1/2の周波数との差が、所定値（０．２５Hz）以下の場合である。
【００４７】
・ステップS17の判断がYesの場合には、脈拍周波数としてfS1を採用して、脈拍数（PR）をPR=fS1×６０で算出する。（ステップS19）
【００４８】
・ステップS16の判断がNoの場合には、
|fSn-fN|＜Δf1(Hz)及び |fSn-fS2/2|＞Δf2(Hz) AND
|fSn-fS3/2|＜Δf2(Hz)か否かの判断及び |fS2/2-前時間窓脈拍f|／前時間窓脈拍f＜ｒ（％）否かの判断をする。（ステップS18）
ここで、fSnは、信号成分周波数スペクトルの最大振幅ピーク周波数、fNは、雑音成分周波数スペクトルの最大振幅ピーク周波数、fS2及びfS3は、信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数、Δf1及びΔf2の値は、例えば、共に０．２５（Hz)であり、ｒの値は１０（％）である。
【００４９】
このステップS18の判断は、図３の表２の優先順位６のFFTパターンに相当する場合であって、図２の表１の（２）における信号成分周波数スペクトルに現れた雑音（体動）Snの周波数(fSn)と雑音成分Nの周波数(fN)との差が、所定値（０．２５Hz）以下の場合であり、信号成分の周波数スペクトルに現れた雑音（体動）Snの周波数(fSn)と第２高調波S2の周波数(fS2)の1/2の周波数との差が、所定値（０．２５Hz）以上の場合である。
また、信号成分の周波数スペクトルに現れた雑音（体動）Snの周波数(fSn)と第３高調波S3の周波数(fS3)の1/2の周波数との差が、所定値（０．２５Hz）以上の場合であり、第２高調波S2の周波数(fS2)の1/2と前時間窓脈拍f との差と前時間窓脈拍ｆとの比が所定％（１０％）以下の場合である。
【００５０】
・ステップS18の判断がYesの場合には、脈拍周波数としてfS2/2を採用して、脈拍数（PR）をPR=fS2/2×６０で算出する。（ステップS20）
・ステップS18の判断がNoの場合には、脈拍周波数を決定せずに、ステップS12に戻る。
【００５１】
図４のフローチャートでは、脈拍数を算出するまでが開示されているが、酸素飽和度（SpO2）は、脈拍数を算出するのに使用した脈拍周波数における被測定媒体からほぼ同時に抽出された２つの信号（IR,R）の比を抽出することによって得られる。
【００５２】
本発明の「運動アルゴリズムで決定された脈拍の周波数決定方法」と、従来の「ゼロクロス法」脈拍数の算出とを併用した例を図５に示す。
図５では、「運動アルゴリズムで決定された脈拍の周波数決定方法」による測定装置１と「ゼロクロス法」による測定装置２とを併用して、ＰＲの選択手段３で運動時には、１を選択し、安静時には２を選択する。
【００５３】
この選択によって、運動時に体動による動脈血以外の組織や静脈血の動きがある場合には、その動きによる周波数帯域の信号を用いないので、正確な酸素飽和度の測定が可能になる。
また、体動ない安静時には、高い時間分解能の脈拍計測が可能となる従来のゼロクロス法等の測定方法（装置）を併用して、選択して使用することができる。その場合には、脈拍数が低い安静時の計測に十分な時間分解能が得られる。また、自律神経の検査のための脈拍数のゆらぎ計測に必要な時間分解能を得られる。
【図面の簡単な説明】
【００５４】
【図１】本発明の脈拍数及び酸素飽和度の測定方法及び測定装置におけるデータ処理の全体的な構成を示すフローチャートである。
【図２】被測定媒体からほぼ同時に抽出された２つの信号（IR,R）からそれぞれ信号成分と雑音成分に分離した後、それぞれの周波数スペクトルを求めた例で、体動があった場合の種々の例を示す図である。
【図３】図２の表に示した例における脈拍数を算出するための判断例を表に纏めたものである。
【図４】本発明の脈拍数（PR）の詳細な算出手順を示すフローチャートである。
【図５】「運動アルゴリズムで決定された脈拍の周波数決定方法」と「ゼロクロス法」による脈拍数の算出を併用した例を示す図である。
【符号の説明】
【００５５】
１：「運動アルゴリズムで決定された脈拍の周波数決定方法」による脈拍数及び酸素飽和度測定装置
２：「ゼロクロス法」による脈拍数及び酸素飽和度測定装置
３：選択手段
【特許請求の範囲】
【請求項１】
生体組織を透過又は反射した２種類の脈波信号を信号成分と雑音成分に分離するステップと、
前記分離された信号成分と雑音成分の周波数スペクトルを演算するステップと、
前記信号成分と雑音成分の周波数スペクトルに基づいて、
体動の有無を判定するステップと、
前記体動の有無の判定結果に応じて、脈拍周波数を決定するステップと、
前記決定された脈拍周波数から脈拍数を演算するステップと、
前記決定された脈拍周波数に対応する前記２種類の脈波信号のスペクトル比から酸素飽和度を演算するステップと、
を含むことを特徴とする運動時の脈拍数及び酸素飽和度測定方法。
【請求項２】
前記脈拍周波数の決定のステップは、
体動が無い場合には、前記信号成分の基本周波数を前記脈拍周波数として決定するステップを、
含むことを特徴とする請求項１に記載の運動時の脈拍数及び酸素飽和度測定方法。
【請求項３】
前記脈拍周波数の決定のステップは、
体動が有る場合に、
信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以上の場合に、前記信号成分周波数スペクトルの最大振幅ピーク周波数を前記脈拍周波数として決定するステップを、
含むことを特徴とする請求項２に記載の運動時の脈拍数及び酸素飽和度測定方法。
【請求項４】
前記脈拍周波数の決定のステップは、
更に、前記信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以下で、前記信号成分周波数スペクトルの最大振幅ピーク周波数と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以下の場合に、前記信号成分周波数スペクトルの最大振幅ピーク周波数を前記脈拍周波数として決定するステップを、
含むことを特徴とする請求項３に記載の運動時の脈拍数及び酸素飽和度測定方法。
【請求項５】
前記脈拍周波数の決定のステップは、
更に、前記信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以下で、前記信号成分周波数スペクトルの最大振幅ピーク周波数と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以上で、信号成分周波数スペクトルの２番目の振幅のピーク周波数 と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以下の場合に、前記信号成分周波数スペクトルの２番目の振幅のピーク周波数を前記脈拍周波数として決定するステップを、
含むことを特徴とする請求項４に記載の運動時の脈拍数及び酸素飽和度測定方法。
【請求項６】
前記脈拍周波数の決定のステップは、
更に、前記信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以下で、前記信号成分周波数スペクトルの最大振幅ピーク周波数と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以上で、
前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２と前時間窓脈拍周波数との差と、前記前時間窓脈拍周波数との比が所定％以下の場合には、前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２を脈拍周波数として決定するステップを、
含むことを特徴とする請求項５に記載の運動時の脈拍数及び酸素飽和度測定方法。
【請求項７】
生体組織を透過又は反射した２種類の脈波信号を信号成分と雑音成分に分離する手段
前記分離された信号成分と雑音成分の周波数スペクトルを演算する手段と、
前記信号成分と雑音成分の周波数スペクトルに基づいて、体動の有無を判定する手段と、
前記体動の有無の判断結果に応じて、前記信号成分と雑音成分の周波数スペクトルに基づいて、脈拍周波数を決定する手段と、
前記決定された脈拍周波数から脈拍数を演算する手段と、
前記決定された脈拍周波数に対応する前記２種類の脈波信号のスペクトル比から酸素飽和度を演算する手段と、
を備えることを特徴とする運動時の脈拍数及び酸素飽和度測定装置。
【請求項１】
生体組織を透過又は反射した２種類の脈波信号を信号成分と雑音成分に分離するステップと、
前記分離された信号成分と雑音成分の周波数スペクトルを演算するステップと、
前記信号成分と雑音成分の周波数スペクトルに基づいて、
体動の有無を判定するステップと、
前記体動の有無の判定結果に応じて、脈拍周波数を決定するステップと、
前記決定された脈拍周波数から脈拍数を演算するステップと、
前記決定された脈拍周波数に対応する前記２種類の脈波信号のスペクトル比から酸素飽和度を演算するステップと、
を含むことを特徴とする運動時の脈拍数及び酸素飽和度測定方法。
【請求項２】
前記脈拍周波数の決定のステップは、
体動が無い場合には、前記信号成分の基本周波数を前記脈拍周波数として決定するステップを、
含むことを特徴とする請求項１に記載の運動時の脈拍数及び酸素飽和度測定方法。
【請求項３】
前記脈拍周波数の決定のステップは、
体動が有る場合に、
信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以上の場合に、前記信号成分周波数スペクトルの最大振幅ピーク周波数を前記脈拍周波数として決定するステップを、
含むことを特徴とする請求項２に記載の運動時の脈拍数及び酸素飽和度測定方法。
【請求項４】
前記脈拍周波数の決定のステップは、
更に、前記信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以下で、前記信号成分周波数スペクトルの最大振幅ピーク周波数と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以下の場合に、前記信号成分周波数スペクトルの最大振幅ピーク周波数を前記脈拍周波数として決定するステップを、
含むことを特徴とする請求項３に記載の運動時の脈拍数及び酸素飽和度測定方法。
【請求項５】
前記脈拍周波数の決定のステップは、
更に、前記信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以下で、前記信号成分周波数スペクトルの最大振幅ピーク周波数と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以上で、信号成分周波数スペクトルの２番目の振幅のピーク周波数 と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以下の場合に、前記信号成分周波数スペクトルの２番目の振幅のピーク周波数を前記脈拍周波数として決定するステップを、
含むことを特徴とする請求項４に記載の運動時の脈拍数及び酸素飽和度測定方法。
【請求項６】
前記脈拍周波数の決定のステップは、
更に、前記信号成分周波数スペクトルの最大振幅ピーク周波数と雑音成分周波数スペクトルの最大振幅ピーク周波数との差が所定値以下で、前記信号成分周波数スペクトルの最大振幅ピーク周波数と前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２との差が所定値以上で、
前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２と前時間窓脈拍周波数との差と、前記前時間窓脈拍周波数との比が所定％以下の場合には、前記信号成分周波数スペクトルの最大振幅ピーク周波数よりも高い周波数のピーク周波数の１／２を脈拍周波数として決定するステップを、
含むことを特徴とする請求項５に記載の運動時の脈拍数及び酸素飽和度測定方法。
【請求項７】
生体組織を透過又は反射した２種類の脈波信号を信号成分と雑音成分に分離する手段
前記分離された信号成分と雑音成分の周波数スペクトルを演算する手段と、
前記信号成分と雑音成分の周波数スペクトルに基づいて、体動の有無を判定する手段と、
前記体動の有無の判断結果に応じて、前記信号成分と雑音成分の周波数スペクトルに基づいて、脈拍周波数を決定する手段と、
前記決定された脈拍周波数から脈拍数を演算する手段と、
前記決定された脈拍周波数に対応する前記２種類の脈波信号のスペクトル比から酸素飽和度を演算する手段と、
を備えることを特徴とする運動時の脈拍数及び酸素飽和度測定装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図２】
【図３】
【図４】
【図５】
【公開番号】特開２００９−２２４８４（Ｐ２００９−２２４８４Ａ）
【公開日】平成２１年２月５日（２００９．２．５）
【国際特許分類】
【出願番号】特願２００７−１８７８８４（Ｐ２００７−１８７８８４）
【出願日】平成１９年７月１９日（２００７．７．１９）
【出願人】（０００２３０９６２）日本光電工業株式会社 (179)
【Ｆターム（参考）】
【公開日】平成２１年２月５日（２００９．２．５）
【国際特許分類】
【出願日】平成１９年７月１９日（２００７．７．１９）
【出願人】（０００２３０９６２）日本光電工業株式会社 (179)
【Ｆターム（参考）】
[ Back to top ]