JP6098261B2 - 信号処理装置、脈波測定装置および信号処理方法 - Google Patents

Description

本発明は、生体の脈波を測定する装置に関する。

ノイズ成分を除去するための演算処理を行う脈波測定装置が開発されている。例えば、特許文献１は、発光部と受光部とを複数設けて測定した複数の脈波信号を演算処理し、ノイズ成分を除去する方法が記載されている。

特開２００５−３２３９０６号公報

しかし、特許文献１に開示された技術は、第１受光部で脈波を含む成分を測定し、第２受光部で“生体情報を含まない”とされる成分を測定する。そしてこの技術は、両者を正規化した後に第１受光部で得た測定結果から第２受光部で得た測定結果を減算している。しかし、第２受光部が第１受光部における脈波以外の成分を測定するには、第２受光部を第１受光部の近傍に設置する必要があるため、第２受光部で測定される成分を、完全に“生体情報を含まない”ものとすることは困難である。そして、第２受光部に脈波が含まれる場合に減算処理によって脈波が消滅してしまう場合がある。さらに、任意の周波数で正規化する処理を行う際に、脈波が現れている周波数を選択してしまうと、この脈波が消滅する場合がある。

本発明は、脈波を含む複数の信号からノイズを除去する技術を提供する。

上述した課題を解決するため、本発明に係る信号処理装置は、生体の脈波を測定する第１測定部から、当該脈波を示す第１信号を取得する第１取得部と、前記第１測定部と異なる感度で前記生体の脈波を測定する第２測定部から、当該脈波を示す第２信号を取得する第２取得部と、前記第１信号の第１スペクトルと第２信号の第２スペクトルから、前記第１測定部と前記第２測定部との感度の比を推算する推算部と、前記推算部により推算された比を用いて、前記第１スペクトルと前記第２スペクトルに含まれるノイズを相殺するように当該第１スペクトルから当該第２スペクトルを減算する減算部とを備えることを特徴とする。
この構成によれば、脈波を含む複数の信号からノイズを除去することができる。

好ましくは、前記推算部は、予め決められた周波数帯における前記第１スペクトルの積分値と、前記周波数帯における前記第２スペクトルの積分値との比を前記感度の比として推算するとよい。
また、上述の態様において、前記周波数帯は、前記脈波に相当する周波数帯であることが望ましい。
また、上述の態様において、前記周波数帯は、０．５Ｈｚ以上３．５Ｈｚ以下の範囲であることが望ましい。
これらの構成によれば、脈波を含む複数の信号からノイズを除去することができる。

また、好ましくは、前記推算部は、前記第１スペクトルにおける所定の周波数帯のスペクトル強度と、前記第２スペクトルにおける当該周波数帯のスペクトル強度との比を前記感度の比として推算するとよい。
この構成によれば、予めノイズが発生する周波数帯が分かっている場合に、高精度でノイズを除去することができる。

また、好ましくは、前記推算部は、前記第１スペクトルと前記第２スペクトルのうち、いずれか一方のスペクトルにおける閾値以上のスペクトル強度を示す周波数帯のスペクトル強度と、他方のスペクトルにおける当該周波数帯のスペクトル強度との比を前記感度の比として推算するとよい。
この構成によれば、脈波よりもノイズの方が必ず強く測定される場合に、高精度でノイズを除去することができる。

また、好ましくは、前記推算部は、前記第１スペクトルと前記第２スペクトルのうち、いずれか一方のスペクトルにおいて強い順に選ばれたスペクトル強度を示す周波数帯のスペクトル強度と、他方のスペクトルにおける当該周波数帯のスペクトル強度との比を前記感度の比として推算するとよい。
この構成によれば、ノイズが発生する周波数帯が不明であるが、脈波よりもノイズの方が必ず強く測定される場合に、高精度でノイズを除去することができる。

また、上述の態様において、前記第１測定部と前記第２測定部との感度の比として基準となる基準値を記憶する記憶部と、前記推算部により推算された比と、前記記憶部に記憶された基準値とに基づいて、前記第１スペクトルおよび前記第２スペクトルに決められた割合を超えるノイズが含まれるか否かを判定する判定部とを備え、前記減算部は、前記判定部によって前記第１スペクトルと前記第２スペクトルに前記割合を超えるノイズが含まれないと判定された場合に、減算を行わないことが望ましい。
この構成によれば、不要な演算処理を行わなくて済む。

また、本発明に係る脈波測定装置は、生体の脈波を測定する第１測定部と、前記第１測定部と異なる感度で前記生体の脈波を測定する第２測定部と、上述の信号処理装置とを備えることを特徴とする。
この構成によれば、脈波を含む複数の信号からノイズを除去して、その脈波を測定することができる。

また、本発明に係る信号処理方法は、第１取得部が、生体の脈波を測定する第１測定部から、当該脈波を示す第１信号を取得する第１取得ステップと、第２取得部が、前記第１測定部と異なる感度で前記生体の脈波を測定する第２測定部から、当該脈波を示す第２信号を取得する第２取得ステップと、推算部が、前記第１信号の第１スペクトルと前記第２信号の前記第２スペクトルから、前記第１測定部と前記第２測定部との感度の比を推算する推算ステップと、減算部が、前記推算ステップにおいて推算された比を用いて、前記第１スペクトルと前記第２のスペクトルに含まれるノイズを相殺するように当該第１スペクトルから当該第２スペクトルを減算する減算ステップとを備える。
この構成によれば、脈波を含む複数の信号からノイズを除去することができる。

脈波測定装置の外観を示す図である。 脈波測定装置の構成を示すブロック図である。 測定部の各構成の配置を示した図である。 制御部の機能的構成を示す図である。 脈波測定装置の動作を示すフロー図である。 脈波測定装置によるノイズ成分の除去の状態を示す図である。

１．実施形態
１−１．全体構成
図１は、脈波測定装置１の外観を示す図である。脈波測定装置１は、リストバンド２によって使用者の手首に固定される腕時計のような構造を有している。脈波測定装置１のおもて面には液晶パネルなどからなる表示面１４１（後述）が設けられており、脈波測定装置１の側面には使用者が押して操作をするボタンスイッチなどの操作子１５１（後述）が設けられている。脈波測定装置１のうら面は使用者の手首に接触する。

図２は、脈波測定装置１の構成を示すブロック図である。制御部１１は、ＣＰＵ（Central Processing Unit）、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）を有し、ＣＰＵがＲＯＭや記憶部１２に記憶されているコンピュータプログラム（以下、単にプログラムという）を読み出して実行することにより脈波測定装置１の各部を制御する。

記憶部１２は、ソリッドステートドライブ（ＳＳＤ：Solid State Drive）などの大容量の記憶手段であり、ＣＰＵに読み込まれるプログラムを記憶する。また、記憶部１２は、ノイズの割合が、所定の閾値よりも下回る状態（以下、基準状態という）で予め測定された脈波に関する情報である基準情報１２１を記憶している。基準情報１２１は、例えば、基準状態で測定した２種類の脈波信号に対し、それぞれパワー積分値を算出して得られたその各積分値の比などである。

表示部１４は、液晶などを利用した表示面１４１を備え、制御部１１からの指示に応じてこの表示面１４１に画像を表示させる。
操作部１５は各種の指示をするためのボタンスイッチなどの操作子１５１を備えており、ユーザによる操作を受け付けてその操作内容に応じた信号を制御部１１に供給する。なお、操作子１５１は、表示面１４１に重ねられた透明なタッチパネルを含んでいてもよい。

測定部１３は、第１測定部１３１、第２測定部１３２、増幅部１３３、およびＡ／Ｄ変換部１３４を有する。第１測定部１３１は、生体の脈波を測定してその脈波を示す第１信号を出力する構成である。具体的に第１測定部１３１は、第１発光部１３１１と第１受光部１３１２とを有しており、これらが生体（この場合、手首の皮膚面）側に接触するように配置される。

第２測定部１３２は、第１測定部１３１と異なる感度で生体の脈波を測定してその脈波を示す第２信号を出力する構成である。ここで「感度」とは、ノイズ感度であり、脈波成分（脈波による信号の強度）に対するノイズ成分（ノイズによる信号の強度）の割合を表したものである。ここでは、第２測定部１３２の感度の方が、第１測定部１３１の感度よりも高くなるように構成されている。
具体的に第２測定部１３２は、第２発光部１３２１と第２受光部１３２２とを有しており、これらが生体側に接触するように配置される。

増幅部１３３は、第１測定部１３１および第２測定部１３２によってそれぞれ出力された信号を増幅するアンプである。Ａ／Ｄ変換部１３４は、増幅部１３３によって増幅されたアナログ信号をデジタル信号に変換する。

図３は、測定部１３の各構成の配置を示した図である。図３（ａ）に示すように、第１発光部１３１１、第１受光部１３１２、第２発光部１３２１および第２受光部１３２２はいずれも脈波測定装置１において、使用者の手首の皮膚面側に接触するように配置されている。なお、図３（ｂ）に示すように、第１発光部１３１１が第２発光部１３２１を兼ねていてもよい。

第１発光部１３１１は、図示しない電源から供給された電流によって、その電流に応じた光量の光を生体組織に向けて照射する。第１受光部１３１２は、第１発光部１３１１が照射した光のうち、生体組織によって反射した光を受光し、その受光強度に応じた信号を第１信号として出力する。この反射光には、様々なものが含まれ得るが、このうち、血管中のヘモグロビンによって反射した光は、脈波を示すものである。一方、例えば使用者の体の動き（体動）があると反射光は影響され、脈波に無関係なノイズを含む場合がある。

第２発光部１３２１は、供給された電流に応じた光量の光を生体組織に向けて照射する。第２受光部１３２２は、第２発光部１３２１が照射した光のうち、生体組織によって反射した光を受光し、その受光強度に応じた信号を第２信号として出力する。第２測定部１３２は、第１測定部１３１と感度が異なるよう、例えば、以下の（１）〜（３）に示すように構成されている。なお、これらの構成は組み合わされていてもよい。

（１）第１発光部１３１１と第１受光部１３１２との距離と異なる距離で第２発光部１３２１と第２受光部１３２２とを配置する。
（２）第２発光部１３２１が第１発光部１３１１と異なる波長の光を照射する。
（３）第１発光部１３１１または第１受光部１３１２を生体に押しつける圧力とは異なる圧力で、第２発光部１３２１または第２受光部１３２２を生体に押しつける。

例えば、発光部と受光部との距離が異なると、生体への光の浸透深度が変化する。距離が長いほど浸透深度は深くなり、生体に深さ依存のノイズ源がある場合には発光部と受光部との距離の違いによって、第１測定部１３１および第２測定部１３２は、異なる感度でそれぞれ脈波を測定する。

また、照射する光の波長によって血管中のヘモグロビンの吸光係数が異なる。特に酸化ヘモグロビンと還元ヘモグロビンとでは吸光係数が異なり、波長をいずれかの吸光係数に合うように調整することで動脈血か静脈血のいずれか一方の影響が強い脈波信号が得られる。ノイズ源が動脈あるいは静脈のいずれか一方にある場合、照射光の波長を異ならせることで第１測定部１３１および第２測定部１３２は、異なる感度でそれぞれ脈波を測定する。

また、装着面への圧力を変えることにより生体組織が押しつぶされる量や領域は変化する。特に真皮層は毛細血管が多く柔らかい組織であるため低い圧力で押しつぶされる。したがって、真皮層の毛細血管にノイズ源がある場合、真皮層が押しつぶされる大きさの違いに応じて、第１測定部１３１および第２測定部１３２は、異なる感度でそれぞれ脈波を測定する。

１−２．制御部の機能的構成
図４は、制御部１１の機能的構成を示す図である。第１測定部１３１により出力された第１信号、および第２測定部１３２により出力された第２信号は、それぞれ増幅部１３３を経て増幅され、Ａ／Ｄ変換部１３４によってデジタル信号に変換されてから制御部１１に供給される。つまり、制御部１１は、生体の脈波を測定する第１測定部１３１から、その脈波を示す第１信号を取得する第１取得部として機能する。また、制御部１１は、第１測定部１３１と異なる感度で生体の脈波を測定する第２測定部１３２から、その脈波を示す第２信号を取得する第２取得部として機能する。

また、制御部１１は、フレーム分割部１１１、スペクトル算出部１１２、感度比推算部１１３、減算係数算出部１１４、ノイズ判定部１１５、減算フィルタ算出部１１６、スペクトル減算部１１７、波形算出部１１８、および波形合成部１１９として機能する。そして、制御部１１は、これらの機能を備えることにより、取得した第１信号および第２信号に含まれるノイズを除去する処理を行う信号処理装置として機能する。

フレーム分割部１１１は、第１測定部１３１および第２測定部１３２によりそれぞれ出力された第１信号および第２信号を、それぞれ所定の時間間隔ごとに分割する。ここで分割された信号をフレームという。なお、ここでフレーム分割部１１１は、後述するスペクトル算出部１１２により周波数成分の解析がし易くなるように、分割した各フレームに対して窓関数をかける処理を行う。窓関数としては、例えば以下の式（１）で与えられるハニング窓関数ω（ｎ）が用いられる。

ここで、Ｎは１フレームあたりのサンプル数であり、ｎはフレーム内におけるサンプル位置を表す。

スペクトル算出部１１２は、フレーム分割部１１１によって分割された各信号を高速フーリエ変換などのアルゴリズムを用いて処理することによりそれぞれスペクトル情報に変換する。このスペクトル情報により、各信号に含まれる周波数成分の強度分布が得られる。

感度比推算部１１３は、第１信号および第２信号からそれぞれ変換されたスペクトル情報を用いて、予め決められた周波数帯である「対象帯域」について各スペクトルを積分する。そして、感度比推算部１１３は、積分により得られた各積分値に基づいて、ノイズ感度比ＮＲを推算する。ここで対象帯域とは、脈波測定装置１が想定している脈拍（例えば、４０〜２００拍／分）に相当する周波数帯の成分であり、具体的には例えば０．６７Ｈｚから３．３３Ｈｚまでである。なお、対象帯域は、０．５Ｈｚ以上３．５Ｈｚ以下の範囲であることが望ましい。

ノイズ感度比ＮＲとは、第１測定部１３１の感度に対する第２測定部１３２の感度の比を推算した値である。ノイズ感度比ＮＲは、様々な方法によって推算されるが、ここでは積分値を用いた推算方法を用いる。

第１信号のスペクトル（以下、第１スペクトルという）に含まれる対象帯域の積分値をＰ₁、第２信号のスペクトル（以下、第２スペクトルという）に含まれる対象帯域の積分値をＰ₂とすると、ノイズ感度比ＮＲは以下の式（２）で表される。すなわち、感度比推算部１１３は、式（２）に沿ってノイズ感度比ＮＲを推算する。

すなわち、感度比推算部１１３は、第１スペクトルおよび第２スペクトルに基づいて、第１測定部１３１と第２測定部１３２との感度の比を推算する推算部として機能する。

減算係数算出部１１４は、感度比推算部１１３によって推算されたノイズ感度比ＮＲを用いて減算係数αを算出する。減算係数αはノイズ感度比ＮＲの逆数であるため、以下の式（３）で表される。すなわち、減算係数算出部１１４は、式（３）に沿って減算係数αを算出する。

ノイズ判定部１１５は、第１スペクトルおよび第２スペクトルに決められた割合を超えるノイズが含まれるか否かを判定する。具体的には、ノイズ判定部１１５は、記憶部１２に記憶された基準情報１２１を参照し、減算係数算出部１１４で算出した減算係数αを検査することで上述した割合を超えるノイズが含まれるか否かを判定する。

すなわち、基準情報１２１を記憶する記憶部１２は、第１測定部１３１と第２測定部１３２との感度の比として基準となる基準値を記憶する記憶部として機能する。
また、ノイズ判定部１１５は、推算部（感度比推算部１１３）により推算された比（ノイズ感度比ＮＲ）と、記憶部１２に記憶された基準値（基準情報１２１）とに基づいて、第１スペクトルおよび第２スペクトルに決められた割合を超えるノイズが含まれるか否かを判定する判定部として機能する。

なお、基準情報１２１には、予め基準状態で測定された第１信号および第２信号の各スペクトルの積分値、またはこれらの積分値の比である脈波係数βが記述されている。ノイズ判定部１１５は、例えば、減算係数算出部１１４が算出した減算係数αを脈波係数βで除算し、その値である（α／β）が予め定められた下限値Ｌを上回り、且つ、予め定められた上限値Ｈを下回るか否かを判定する。そして、（α／β）が下限値Ｌを上回り、且つ、上限値Ｈを下回ると判定した場合に、ノイズ判定部１１５は、第１スペクトルおよび第２スペクトルに決められた割合を超えるノイズが含まれないと判定する。

減算フィルタ算出部１１６は、減算係数算出部１１４により算出された減算係数α、第１スペクトルであるＸ₁（ω）、および第２スペクトルであるＸ₂（ω）によりスペクトル減算のフィルタ係数Ｈ（ω）を、次の式（４）に沿って算出する。

スペクトル減算部１１７は、第１スペクトルおよび第２スペクトルに減算フィルタ算出部１１６が算出したフィルタ係数を適用することで、各スペクトルの減算を行う。具体的に、Ｘ₁（ω）から減算係数αをかけたＸ₂（ω）を減算する処理において、減算後のスペクトルをＳ（ω）とすれば、Ｓ（ω）は次の式（５）によって表される。

ここで、ｘ（ｔ）とは、第１測定部１３１において測定される第１信号である。つまり、減算処理は、第１信号をフィルタＨ（ω）でフィルタリングすることと同じである。この減算処理の結果、各スペクトルに含まれるノイズが相殺されたスペクトルＳ（ω）が得られる。すなわち、減算係数算出部１１４、減算フィルタ算出部１１６、およびスペクトル減算部１１７は、推算部（感度比推算部１１３）により推算された比（ノイズ感度比ＮＲ）を用いて、各スペクトルに含まれるノイズを相殺するように各スペクトルのうち一方から他方を減算する減算部として機能する。なお、この場合、減算係数算出部１１４、減算フィルタ算出部１１６、およびスペクトル減算部１１７は、第１スペクトルから第２スペクトルを減算することで、各スペクトルに含まれるノイズを相殺する。また、ここでいう「ノイズを相殺するように減算する」の意味は、各スペクトルに含まれていたノイズが完全に無くなる場合に限らず、それらのノイズが減算の前に比べて低減する場合も含む。

ただし、ノイズ判定部１１５においてスペクトルに決められた割合を超えるノイズが含まれていないと判定された場合、スペクトル減算部１１７は上述した減算処理を行わず、上記の各スペクトルのうち、いずれか一方をそのまま出力する。なお、出力する方のスペクトルは、予め決められていてもよいし、例えば、相対的に感度の低い方が選ばれてもよい。

波形算出部１１８は、スペクトル減算部１１７の出力を逆フーリエ変換処理することで時間領域の波形を算出する。
波形合成部１１９は、波形算出部１１８によって出力される波形を合成して出力する。

１−３．動作
図５は、脈波測定装置１の動作を示すフロー図である。脈波測定装置１の制御部１１は、測定部１３から第１信号および第２信号を受け取ると、これら各信号をフレームに分割し（ステップＳ０１）、分割した各フレームに対して窓関数をかける処理（窓関数処理）を行う（ステップＳ０２）。そして制御部１１は、窓関数処理を施された各フレームに対して高速フーリエ変換などを行ってスペクトルを算出し（ステップＳ０３）、各スペクトルの対象帯域をそれぞれ積分して（ステップＳ０４）、得られた積分値からノイズ感度比ＮＲを推算する（ステップＳ０５）。また制御部１１は、推算したノイズ感度比ＮＲに基づき、減算係数αを算出する（ステップＳ０６）。

次に、制御部１１は、基準情報１２１から読み出される脈波係数βと、ステップＳ０６において算出した減算係数αとに基づいて、第１スペクトルおよび第２スペクトルに決められた割合を超えるノイズが含まれるか否かを判定する（ステップＳ０７）。ノイズが含まれると判定した場合（ステップＳ０７；ＹＥＳ）、制御部１１は、算出された減算係数αによりスペクトル減算のフィルタ係数を算出し（ステップＳ０８）、算出したフィルタ係数を適用することで、各スペクトルの減算を行う（ステップＳ０９）。そして、制御部１１は、減算されたスペクトルを用いて逆フーリエ変換を行って時間信号を算出し（ステップＳ１０）、さらにこの時間信号を重複加算した合成波形を算出する（ステップＳ１１）。一方、ノイズが含まれないと判定した場合（ステップＳ０７；ＮＯ）、制御部１１は、ステップＳ０１で分離したフレームをそのまま出力する（ステップＳ１２）。

図６は、脈波測定装置１によるノイズ成分の除去の状態を示す図である。図６（ａ）に示す実線は第１測定部１３１により出力された第１信号を変換して得られた第１スペクトルであり、破線は第２測定部１３２により出力された第２信号を変換して得られた第２スペクトルである。脈波測定装置１が、第１測定部１３１と第２測定部１３２との感度の比を算出し、その比を用いて減算処理を行うと、図６（ｂ）に示すスペクトルが得られる。これは第１信号および第２信号にそれぞれ含まれていたノイズが除去され、対象帯域が抽出された信号である。

以上、説明した通り、脈波測定装置１は、２つの測定部の感度が異なることを利用して、それぞれが出力する信号に含まれるノイズを除去するので、従来の技術に比べて脈波を高精度に測定することができる。

２．変形例
以上が実施形態の説明であるが、この実施形態の内容は以下のように変形し得る。また、以下の変形例を組み合わせてもよい。

２−１．変形例１
上述した実施形態において、感度比推算部１１３は、第１スペクトルの積分値と、第２スペクトルの積分値との比を、ノイズ感度比ＮＲとして推算していたが、ノイズ感度比ＮＲの推算方法はこれに限られない。感度比推算部１１３は、第１スペクトルおよび第２スペクトルのうち、いずれか一方における閾値以上のスペクトル強度を示す周波数帯のスペクトル強度と、他方のスペクトルにおけるその周波数帯のスペクトル強度との比を感度の比として推算してもよい。

例えば、感度比推算部１１３は、第１スペクトルのうち、決められた閾値以上のスペクトル強度（以下、第１スペクトル強度という）を示す周波数帯を特定し、第２スペクトルにおけるその周波数帯のスペクトル強度（以下、第２スペクトル強度という）を特定する。そして、感度比推算部１１３は、第１スペクトル強度に対する第２スペクトル強度の比をノイズ感度比ＮＲとして推算してもよい。

脈波よりもノイズの方が必ず強く測定される場合、強いスペクトル強度ほど、脈波でなくノイズである可能性が高い。上述したように感度比推算部１１３が一方のスペクトルにおいて閾値以上のスペクトル強度を示す周波数帯を特定し、他方におけるその周波数帯のスペクトル強度を特定してこれらスペクトル強度の比を算出することで、真のノイズがどの周波数帯に発生しているか不明な場合であっても、ノイズ感度比ＮＲを高精度で推算することができる。

なお、決められた閾値以上の値を示す第１スペクトル強度が複数ある場合、これらのそれぞれに対応する第２スペクトル強度を特定し、それぞれの組について比を算出して、これら複数の比に基づいてノイズ感度比ＮＲを推算してもよい。例えば、複数の比の相加平均値をノイズ感度比ＮＲとしてもよい。

２−２．変形例２
また、感度比推算部１１３は、第１スペクトルおよび第２スペクトルのうち、いずれか一方のスペクトルにおいて強い順に選ばれたスペクトル強度を示す周波数帯のスペクトル強度と、他方のスペクトルにおけるその周波数帯のスペクトル強度との比を感度の比として推算してもよい。

例えば、感度比推算部１１３は、第１スペクトルのうち、強度の強い順に決められた数の第１スペクトル強度を選び、これら第１スペクトル強度を示す周波数帯を特定する。そして、感度比推算部１１３は、第２スペクトルにおけるその周波数帯のスペクトル強度を第２スペクトル強度として特定する。そして、感度比推算部１１３は、第１スペクトル強度に対する第２スペクトル強度の比をノイズ感度比ＮＲとして推算してもよい。この構成によっても、ノイズ感度比ＮＲを高精度で推算することができる。

２−３．変形例３
また、感度比推算部１１３は、第１スペクトルにおける所定の周波数帯のスペクトル強度と、第２スペクトルにおけるその周波数帯のスペクトル強度との比を感度の比として推算してもよい。例えば、感度比推算部１１３は、第１スペクトルおよび第２スペクトルにおいて、３．６Ｈｚから４．０Ｈｚまでの範囲などのように、上述した対象帯域から外れた帯域（所定の周波数帯）におけるスペクトル強度を、それぞれ特定する。そして、感度比推算部１１３は、特定したこれらスペクトル強度の比をノイズ感度比ＮＲとして推算してもよい。

この構成によれば、予めノイズが発生する周波数帯が分かっている場合に、ノイズ感度比ＮＲを高精度で推算することができる。

２−４．変形例４
上述した実施形態において、ノイズ判定部１１５は、基準情報１２１を参照し、減算係数算出部１１４で算出した減算係数αを検査することで第１スペクトルおよび第２スペクトルに決められた割合を超えるノイズが含まれるか否かを判定していたが、制御部１１は、ノイズ判定部１１５として機能しなくてもよい。この場合、記憶部１２は基準情報１２１を記憶しなくてもよい。そして、この場合、減算フィルタ算出部１１６は、減算係数算出部１１４が算出した減算係数αがいかなる値であっても、スペクトル減算のフィルタ係数を算出し、スペクトル減算部１１７は、このフィルタ係数を第１スペクトルおよび第２スペクトルに適用することで、各スペクトルの減算を行えばよい。

１…脈波測定装置、１１…制御部（信号処理装置）、１１１…フレーム分割部、１１２…スペクトル算出部、１１３…感度比推算部、１１４…減算係数算出部、１１５…ノイズ判定部、１１６…減算フィルタ算出部、１１７…スペクトル減算部、１１８…波形算出部、１１９…波形合成部、１２…記憶部、１２１…基準情報、１３…測定部、１３１…第１測定部、１３１１…第１発光部、１３１２…第１受光部、１３２…第２測定部、１３２１…第２発光部、１３２２…第２受光部、１３３…増幅部、１３４…Ａ／Ｄ変換部、１４…表示部、１４１…表示面、１５…操作部、１５１…操作子、２…リストバンド。

Claims (10)

1. 生体の脈波を測定する第１測定部から、前記脈波を示す第１信号を取得する第１取得部と、
   前記第１測定部と異なる感度で前記生体の脈波を測定する第２測定部から、前記脈波を示す第２信号を取得する第２取得部と、
   前記第１信号の第１スペクトルと前記第２信号の第２スペクトルから、予め決められた周波数帯における前記第１スペクトルの積分値と、前記周波数帯における前記第２スペクトルの積分値との比を、前記第１測定部と前記第２測定部との感度の比をとして推算する推算部と、
   前記推算部により推算された比を用いて、前記第１スペクトルと前記第２スペクトルに含まれるノイズを相殺するように当該第１スペクトルから当該第２スペクトルを減算する減算部と、
   を備えることを特徴とする信号処理装置。
2. 前記周波数帯は、脈拍に相当する周波数帯である
   ことを特徴とする請求項１に記載の信号処理装置。
3. 前記周波数帯は、０．５Ｈｚ以上３．５Ｈｚ以下の範囲である
   ことを特徴とする請求項２に記載の信号処理装置。
4. 生体の脈波を測定する第１測定部から、当該脈波を示す第１信号を取得する第１取得部と、
   前記第１測定部と異なる感度で前記生体の脈波を測定する第２測定部から、当該脈波を示す第２信号を取得する第２取得部と、
   前記第１信号の第１スペクトルと第２信号の第２スペクトルのうち、いずれか一方のスペクトルにおける閾値以上のスペクトル強度を示す周波数帯のスペクトル強度と、他方のスペクトルにおける当該周波数帯のスペクトル強度との比を前記第１測定部と前記第２測定部との感度の比として推算する推算部と、
   前記推算部により推算された比を用いて、前記第１スペクトルと前記第２スペクトルに含まれるノイズを相殺するように当該第１スペクトルから当該第２スペクトルを減算する減算部と、
   を備えることを特徴とする信号処理装置。
5. 生体の脈波を測定する第１測定部から、当該脈波を示す第１信号を取得する第１取得部と、
   前記第１測定部と異なる感度で前記生体の脈波を測定する第２測定部から、当該脈波を示す第２信号を取得する第２取得部と、
   前記第１信号の第１スペクトルと第２信号の第２スペクトルのうち、いずれか一方のスペクトルにおいて強い順に選ばれたスペクトル強度を示す周波数帯のスペクトル強度と、他方のスペクトルにおける当該周波数帯のスペクトル強度との比を前記第１測定部と前記第２測定部との感度の比として推算する推算部と、
   前記推算部により推算された比を用いて、前記第１スペクトルと前記第２スペクトルに含まれるノイズを相殺するように当該第１スペクトルから当該第２スペクトルを減算する減算部と、
   を備えることを特徴とする信号処理装置。
6. 前記第１測定部と前記第２測定部との感度の比として基準となる基準値を記憶する記憶部と、
   前記推算部により推算された比と、前記記憶部に記憶された基準値とに基づいて、前記第１スペクトルおよび前記第２スペクトルに決められた割合を超えるノイズが含まれるか否かを判定する判定部とを備え、
   前記減算部は、前記判定部によって前記第１スペクトルと前記第２スペクトルに前記割合を超えるノイズが含まれないと判定された場合に、減算を行わない
   ことを特徴とする請求項１から５のいずれか１項に記載の信号処理装置。
7. 生体の脈波を測定する第１測定部と、
   前記第１測定部と異なる感度で前記生体の脈波を測定する第２測定部と、
   請求項１から６のいずれか１項に記載の信号処理装置と
   を備えることを特徴とする脈波測定装置。
8. 第１取得部が、生体の脈波を測定する第１測定部から、前記脈波を示す第１信号を取得する第１取得ステップと、
   第２取得部が、前記第１測定部と異なる感度で前記生体の脈波を測定する第２測定部から、前記脈波を示す第２信号を取得する第２取得ステップと、
   推算部が、前記第１信号の第１スペクトルと前記第２信号の第２スペクトルから、予め決められた周波数帯における前記第１スペクトルの積分値と、前記周波数帯における前記第２スペクトルの積分値との比を、前記第１測定部と前記第２測定部との感度の比として推算する推算ステップと、
   減算部が、前記推算ステップにおいて推算された比を用いて、前記第１スペクトルと前記第２のスペクトルに含まれるノイズを相殺するように当該第１スペクトルから当該第２スペクトルを減算する減算ステップと、
   を備える信号処理方法。
9. 第１取得部が、生体の脈波を測定する第１測定部から、前記脈波を示す第１信号を取得する第１取得ステップと、
   第２取得部が、前記第１測定部と異なる感度で前記生体の脈波を測定する第２測定部から、前記脈波を示す第２信号を取得する第２取得ステップと、
   推算部が、前記第１信号の第１スペクトルと第２信号の第２スペクトルのうち、いずれか一方のスペクトルにおける閾値以上のスペクトル強度を示す周波数帯のスペクトル強度と、他方のスペクトルにおける当該周波数帯のスペクトル強度との比を前記第１測定部と前記第２測定部との感度の比として推算する推算ステップと、
   減算部が、前記推算ステップにおいて推算された比を用いて、前記第１スペクトルと前記第２のスペクトルに含まれるノイズを相殺するように当該第１スペクトルから当該第２スペクトルを減算する減算ステップと、
   を備える信号処理方法。
10. 第１取得部が、生体の脈波を測定する第１測定部から、前記脈波を示す第１信号を取得する第１取得ステップと、
    第２取得部が、前記第１測定部と異なる感度で前記生体の脈波を測定する第２測定部から、前記脈波を示す第２信号を取得する第２取得ステップと、
    推算部が、前記第１信号の第１スペクトルと第２信号の第２スペクトルのうち、いずれか一方のスペクトルにおいて強い順に選ばれたスペクトル強度を示す周波数帯のスペクトル強度と、他方のスペクトルにおける当該周波数帯のスペクトル強度との比を前記第１測定部と前記第２測定部との感度の比として推算する推算ステップと、
    減算部が、前記推算ステップにおいて推算された比を用いて、前記第１スペクトルと前記第２のスペクトルに含まれるノイズを相殺するように当該第１スペクトルから当該第２スペクトルを減算する減算ステップと、
    を備える信号処理方法。

Images