(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2022-04-27
(45)【発行日】2022-05-11
(54)【発明の名称】心拍検出システム、心拍検出方法
(51)【国際特許分類】
A61B 5/11 20060101AFI20220428BHJP
A61B 5/0245 20060101ALI20220428BHJP
【ＦＩ】
A61B5/11 110
A61B5/0245 Z
A61B5/11 ZDM
【請求項の数】
10
(21)【出願番号】P 2018014373
(22)【出願日】2018-01-31
(65)【公開番号】P2019129996
(43)【公開日】2019-08-08
【審査請求日】2021-01-22
【新規性喪失の例外の表示】特許法第３０条第２項適用 ・発行日：平成２９年８月２９日 刊行物：信学技報、ｖｏｌ．１００ ｎｏ．Ｂ－２０－１０、ｐｐ．３９３、２０１７年８月発行、「確率的勾配降下法を用いたドップラーレーダによるロバストな心拍検出」・開催日：平成２９年９月１２日～１５日 集会名、開催場所：２０１７年電子情報通信学会ソサイエティ大会、東京都市大学世田谷キャンパス（東京都世田谷区玉堤１－２８－１）
【国等の委託研究の成果に係る記載事項】（出願人による申告）平成２９年度、国立研究開発法人科学技術振興機構、研究成果展開事業 センター・オブ・イノベーション（ＣＯＩ）プログラム、ＣＯＩ拠点「真の社会イノベーションを実現する革新的「健やか力」創造拠点」委託事業、産業技術力強化法第１９条の適用を受ける特許出願
(73)【特許権者】
【識別番号】521369068
【氏名又は名称】データソリューションズ株式会社
(74)【代理人】
【識別番号】100079108
【弁理士】
【氏名又は名称】稲葉 良幸
(74)【代理人】
【識別番号】100109346
【弁理士】
【氏名又は名称】大貫 敏史
(74)【代理人】
【識別番号】100117189
【弁理士】
【氏名又は名称】江口 昭彦
(74)【代理人】
【識別番号】100134120
【弁理士】
【氏名又は名称】内藤 和彦
(72)【発明者】
【氏名】大槻 知明
(72)【発明者】
【氏名】葉 臣
【審査官】▲瀬▼戸井 綾菜
(56)【参考文献】
【文献】欧州特許出願公開第３１８２４１３（ＥＰ，Ａ１）
【文献】特開２００５－１６０６４０（ＪＰ，Ａ）
【文献】国際公開第２０１７／１４０６６３（ＷＯ，Ａ１）
【文献】特表２０１７－５１９５４８（ＪＰ，Ａ）
【文献】特開２０１０－１４２４５６（ＪＰ，Ａ）
【文献】特開２００６－５５５０４（ＪＰ，Ａ）
【文献】特開２０１７－２２５５６０（ＪＰ，Ａ）
【文献】米国特許出願公開第２０１０／０２４５０９１（ＵＳ，Ａ１）
【文献】国際公開第２０１８／００３１３９（ＷＯ，Ａ１）
【文献】茂木 瑛梨子 他，スペクトログラムを用いたドップラーレーダによる心拍検出，電子情報通信学会技術研究報告，2017年01月，Vol.116, No.407，pp.25-30
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ６１Ｂ ５／０６－５／２２
Ａ６１Ｂ ５／０２－５／０３
Ａ６１Ｂ １０／００
(57)【特許請求の範囲】
【請求項1】
被験者からの反射波を受信して時間領域のＩ信号及びＱ信号を取得するドップラーセンサと、時間領域の前記Ｉ信号及び前記Ｑ信号を周波数領域のＩ信号及びＱ信号に変換し、スパーススペクトルを再構成するスパーススペクトル再構成部と、
前記スパーススペクトルのピークを追跡し、心拍数を推定するピーク追跡部と、
を有し、
前記スパーススペクトル再構成部は、確率的勾配降下法に属するアルゴリズムを用いて心拍のスパーススペクトルを再構成する、
心拍検出システム。
【請求項2】
前記ピーク追跡部は、２つの隣接する心拍数を推定する際に、前回の推定値を中心とした所定の探索範囲内で今回の推定値を選択し、推定値の変化が所定の閾値を超える場合は、前回の推定値を出力する請求項１に記載の心拍検出システム。
【請求項3】
時間領域の前記Ｉ信号及び前記Ｑ信号に含まれる心拍に対する雑音主成分を除去し、心拍成分を抽出する特異スペクトル分析部を有し、前記スパーススペクトル再構成部は、前記特異スペクトル分析部により雑音主成分を除去された時間領域の前記Ｉ信号及び前記Ｑ信号を周波数領域のＩ信号及び前記Ｑ信号に変換する請求項１又は２に記載の心拍検出システム。
【請求項4】
時間領域の前記Ｉ信号及び前記Ｑ信号に含まれる心拍に対する非周期雑音を抑える時間差分処理部を有し、前記スパーススペクトル再構成部は、前記特異スペクトル分析部により雑音主成分を除去され、前記時間差分処理部により非周期雑音を抑えられた時間領域の前記Ｉ信号及び前記Ｑ信号を周波数領域のＩ信号及び前記Ｑ信号に変換する請求項３に記載の心拍検出システム。
【請求項5】
前記確率的勾配降下法に属するアルゴリズムは、ＺＡ－ＳＬＭＳアルゴリズムである請求項１乃至４の何れか一項に記載の心拍検出システム。
【請求項6】
前記確率的勾配降下法に属するアルゴリズムは、ＺＡ－ＬＭＳアルゴリズムである請求項１乃至４の何れか一項に記載の心拍検出システム。
【請求項7】
前記確率的勾配降下法に属するアルゴリズムは、ＦＯＣＵＳＳアルゴリズムである請求項１乃至４の何れか一項に記載の心拍検出システム。
【請求項8】
前記スパーススペクトル再構成部は、Ｉ信号のスペクトルピークとＱ信号のスペクトルピークのうち、ピーク値が高い方を採用する、請求項１乃至７の何れか一項に記載の心拍検出システム。
【請求項9】
サンプリングされた信号をｄ（ｎ）、入力信号をｘ（ｎ）、フィルタ係数をｈ（ｎ）として、下記式で与えられる回帰誤差ｅ（ｎ）を最小化する、
ｅ（ｎ）＝ｄ（ｎ）－ｘ T （ｎ）ｈ（ｎ）
請求項５又は６に記載の心拍検出システム。
【請求項10】
ドップラーセンサで被験者からの反射波を受信して時間領域のＩ信号及びＱ信号を取得するステップと、時間領域の前記Ｉ信号及び前記Ｑ信号を周波数領域のＩ信号及び前記Ｑ信号に変換し、スパーススペクトルを再構成するステップと、
前記スパーススペクトルのピークを追跡し、心拍数を推定するステップと、
を有し、
前記スパーススペクトルを再構成することにおいて、確率的勾配降下法に属するアルゴリズムを用いて心拍のスパーススペクトルを再構成する、
心拍検出方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、心拍検出システム、心拍検出方法に関する。
【背景技術】
【0002】
近年、被験者の健康状態や精神状態を把握するために、例えばドップラーセンサを用いて非接触で心拍数や心拍間隔を検出する心拍検出システムが検討されている。ドップラーセンサを用いた心拍検出システムにおいて、心拍検出精度を向上することは重要な課題であり、そのために様々な提案がなされている。
【0003】
様々な提案のうちの１つとして、例えば、スペクトログラム法が挙げられる。スペクトログラム法は、スペクトログラム上の心拍に反応していると考えられる周波数から心拍成分を含む周波数帯域を決定し、周波数帯域内の心拍による周波数の振幅を積分し、積分値の時間変化から心拍を検出する手法である。
【0004】
スペクトログラム法では、被験者が着座静止時には優れた心拍検出精度が得られるが、タイピング時や運転時等の様々な体動を伴う動作時には、十分な心拍検出精度が得られなかった。
【先行技術文献】
【非特許文献】
【0005】
【文献】茂木，大槻，"スペクトログラムを用いたドップラーレーダによる心拍検出"，電子情報通信学会ＡＳＮ研究会２０１７．
【発明の概要】
【発明が解決しようとする課題】
【0006】
本発明は、上記の点に鑑みてなされたものであり、心拍検出精度を向上することを目的とする。
【課題を解決するための手段】
【0007】
本心拍検出システムは、被験者からの反射波を受信して時間領域のＩ信号及びＱ信号を取得するドップラーセンサと、時間領域の前記Ｉ信号及び前記Ｑ信号を周波数領域のＩ信号及びＱ信号に変換するスパーススペクトル再構成部と、を有し、前記スパーススペクトル再構成部は、確率的勾配降下法に属するアルゴリズムを用いて心拍のスパーススペクトルを再構成することを要件とする。
【発明の効果】
【0008】
開示の技術によれば、心拍検出精度を向上することができる。
【図面の簡単な説明】
【0009】
【図1】第１の実施の形態に係る心拍検出システムの概略構成を例示する図である。
【図2】ドップラーセンサで得た信号波形の一例である。
【図3】第１の実施の形態に係る信号処理部のハードウェアブロックを例示する図である。
【図4】第１の実施の形態に係る信号処理部の機能ブロックを例示する図である。
【図5】第１の実施の形態に係る心拍検出システムの動作を示すフローチャートの一例である。
【図6】所定のステップで得た信号波形の一例である。
【図7】ＳＳＲに用いる適応フィルタについて説明する図である。
【図8】ＳＳＲ問題について説明する図である。
【図9】ＺＡ－ＳＬＭＳアルゴリズムについて説明する図である。
【図10】実施例の代表的な測定結果を示す図（その１）である。
【図11】実施例の代表的な測定結果を示す図（その２）である。
【図12】実施例の代表的な測定結果を示す図（その３）である。
【図13】実施例の代表的な測定結果を示す図（その４）である。
【図14】５人の被験者について着座静止している場合の心拍数の評価結果をまとめた図である。
【図15】５人の被験者について着座静止している場合の心拍間隔の評価結果をまとめた図である。
【図16】５人の被験者についてタイピングしている場合の心拍数の評価結果をまとめた図である。
【図17】５人の被験者についてタイピングしている場合の心拍間隔の評価結果をまとめた図である。
【発明を実施するための形態】
【0010】
以下、図面を参照して、実施の形態の説明を行う。なお、各図面において、同一構成部分には同一符号を付し、重複した説明を省略する場合がある。
【0011】
【0012】
ドップラーセンサ１０は、ドップラー効果による送信信号と受信信号の周波数シフトを観測することで、観測対象（被験者）の動きを検出するセンサである。本実施の形態では、一例として、送信波として無変調連続波（ＣＷ：Continuous Wave）を用いる。
【0013】
ドップラーセンサ１０は被験者の近傍に配置され、被験者の心臓又はその近傍で反射された信号（反射波）を受信し、観測対象の動きにより発生した信号を得ることができる。被験者としては、例えば、ＰＣでタイピングする作業者、車両の運転者等が挙げられる。
【0014】
図２は、ドップラーセンサで得た信号波形の一例である。図２に示す信号は、送信信号と受信信号の間の周波数シフトを表わす信号を時間の関数として取得したものであり、送信信号と同相（In-phase）成分であるＩ信号、及び直交位相（Quadrature）成分であるＱ信号で構成される。
【0015】
心拍や呼吸、瞬き、体動（身体による動作）等の身体の表面の動きはドップラーセンサ１０により観測できる。心拍検出を行う際には、ハイパスフィルタやバンドパスフィルタ等を用いて適宜ノイズ除去を行うことが好ましい。
【0016】
図１に戻り、信号処理部２０は、ドップラーセンサ１０の出力信号に基づいて、被験者の心拍数や心拍間隔を検出する。信号処理部２０は、適宜、ドップラーセンサ１０で受信した信号のＩ信号及びＱ信号をそのまま利用したり、Ｉ信号及びＱ信号に基づいて各種信号（振幅、位相、それらの積分値等）を生成したりすることができる。
【0017】
図３は、第１の実施の形態に係る信号処理部のハードウェアブロックを例示する図である。図３を参照するに、信号処理部２０は、ＣＰＵ２１と、ＲＯＭ２２と、ＲＡＭ２３と、Ｉ／Ｆ２４と、バスライン２５とを有している。ＣＰＵ２１、ＲＯＭ２２、ＲＡＭ２３、及びＩ／Ｆ２４は、バスライン２５を介して相互に接続されている。
【0018】
ＣＰＵ２１は、信号処理部２０の各機能を制御する。記憶手段であるＲＯＭ２２は、ＣＰＵ２１が信号処理部２０の各機能を制御するために実行するプログラムや、各種情報を記憶している。記憶手段であるＲＡＭ２３は、ＣＰＵ２１のワークエリア等として使用される。又、ＲＡＭ２３は、所定の情報を一時的に記憶することができる。Ｉ／Ｆ２４は、心拍検出システム１を他の機器等と接続するためのインターフェイスである。心拍検出システム１は、Ｉ／Ｆ２４を介して、外部ネットワーク等と接続されてもよい。
【0019】
但し、信号処理部２０の一部又は全部は、ハードウェアのみにより実現されてもよい。ハードウェアの一例としては、ＡＳＩＣ（Application Specific Integrated Circuit)、ＤＳＰ（Digital Signal Processor）、ＦＰＧＡ（Field Programmable Gate Array）等が挙げられる。又、信号処理部２０は、物理的に複数の装置等により構成されてもよい。
【0020】
図４は、第１の実施の形態に係る信号処理部の機能ブロックを例示する図である。図４を参照するに、信号処理部２０は、機能ブロックとして、信号取得部２０１と、特異スペクトル分析部２０２と、時間差分処理部２０３と、スパーススペクトル再構成部２０４と、ピーク追跡部２０５とを有している。各機能ブロックの具体的な機能については、図５等の説明の中で後述する。
【0021】
【0022】
まず、ステップＳ１１では、信号取得部２０１は、ドップラーセンサ１０で被験者からの反射波を受信して信号を取得する（信号取得ステップ）。ここで取得する信号は、例えば、図２に示すようなＩ信号及びＱ信号で構成されている。図２に示すＩ信号及びＱ信号はサンプリングされた信号であり、サンプリングレートは、例えば、１０００Ｈｚである。
【0023】
次に、ステップＳ１２では、信号取得部２０１は、ステップＳ１１で受信したＩ信号及びＱ信号のダウンサンプリングを行う。ダウンサンプリングのサンプリングレートは、平常時の心拍（１～１．６Ｈｚ程度）の２倍（３．２Ｈｚ）よりも十分に高ければ任意の値としてよく、例えば、１２５Ｈｚ程度とすることができる。
【0024】
例えば、ステップＳ１１で取得するドップラーセンサ１０のサンプリングレートが１０００Ｈｚであり、Ｉ信号及びＱ信号を取得するタイムウィンドウ（観測時間）が８秒であるとすると、サンプリング数は８０００となる。これを１２５Ｈｚのサンプリングレート（周波数帯域：０～１２５Ｈｚ）でダウンサンプリングすることにより、サンプリング数は１０００となり、信号処理における計算量を大幅に低減することができる。
【0025】
次に、ステップＳ１３では、ステップＳ１２でダウンサンプリングを行ったＩ信号及びＱ信号のノイズ成分の除去を行う。ノイズ成分の除去は、例えば、バンドパスフィルタを用いてハードウェア的に行うことができる。又、信号処理部２０内においてディジタル信号処理（ディジタルフィルタ等）により行ってもよい。一般に、平常時の心拍は１～１．６Ｈｚ程度の範囲で変化するので、バンドパスフィルタの通過帯域は、例えば、０．４Ｈｚ～５Ｈｚ程度とすることができる。
【0026】
図６（ａ）に信号波形の一例を示す。図６（ａ）に示す信号波形は、タイムウィンドウを８秒とした場合に、ステップＳ１２でサンプリングレート１０００ＨｚのＩ信号及びＱ信号を１２５Ｈｚでダウンサンプリングし、ステップＳ１３でノイズ成分の除去を行ったドップラー信号である（つまり、図６（ａ）は８秒間のデータである）。
【0027】
次に、ステップＳ１４では、特異スペクトル分析部２０２は、特異スペクトル分析（ＳＳＡ：Singular Spectrum Analysis）を行う。ＳＳＡは、時間領域のＩ信号及びＱ信号を振動的な成分とノイズ成分とに分解し、時間領域のＩ信号及びＱ信号に含まれる呼吸や体動による雑音主成分（心拍に対する雑音主成分）を除去し、心拍成分を抽出するものである。なお、ステップＳ１４では、必要に応じて、１回前の処理値を用いることができる。例えば、１回前の処理値からの範囲を制限して今回の処理値を求めることにより、精度が向上する場合がある。
【0028】
次に、ステップＳ１５では、時間差分処理部２０３は、時間領域のＩ信号及びＱ信号に含まれる心拍に対する非周期雑音を抑える時間差分処理を行う。具体的には、ステップＳ１４で特異スペクトル分析を行って呼吸や体動による雑音主成分を除去したサンプリング信号ｙ＝［ｙ１、ｙ２、・・・、ｙＭ］を求め（例えば、Ｍ＝１０００）、サンプリング信号ｙの１次差分ｙ'＝［ｙ２－ｙ１、ｙ３－ｙ２、・・・、ｙＭ－ｙＭ－１］を求める。時間差分処理を行うことで、ステップＳ１６で再構成するスペクトルのピークを実際の周期的な心拍のピークに近づかせることができる。
【0029】
図６（ｂ）に信号波形の一例を示す。図６（ｂ）に示す信号波形は、図６（ａ）のＩ信号及びＱ信号に対しステップＳ１４で特異スペクトル分析を行って呼吸や体動による雑音主成分を除去したサンプリング信号ｙを求め（Ｍ＝１０００）、サンプリング信号ｙの１次差分ｙ'を求めたものである。
【0030】
次に、ステップＳ１６では、スパーススペクトル再構成部２０４は、時間領域のＩ信号及びＱ信号を周波数領域のＩ信号及びＱ信号に変換するスパーススペクトル再構成（ＳＳＲ：Sparse Spectrum Reconstruction）を行う（スパーススペクトル再構成ステップ）。スパーススペクトル再構成部２０４は、確率的勾配降下法（ＳＧＤ：Stochastic Gradient Descent）に属するアルゴリズムを用いて心拍のスパーススペクトルを再構成する。本実施の形態では、一例として、確率的勾配降下法に属するＺＡ－ＬＭＳ（Zero Attracting Least Mean Square）アルゴリズムをＳＳＲに適用し、呼吸や体動による雑音を除去する例を示す。
【0031】
なお、スパーススペクトルとは、単に心拍ピークがあり、ほかの周波数に対応するエネルギーがほとんどゼロという（心拍）スペクトルである。
【0032】
図６（ｃ）に信号波形の一例を示す。図６（ｃ）に示す信号波形は、図６（ｂ）のＩ信号及びＱ信号に対しＺＡ－ＬＭＳアルゴリズムによりスパーススペクトル再構成を行った後のＩ信号である。図６（ｃ）において、ＧＴ（Ground Truth：真値）は心電計で測定した心拍数（図６（ｄ）に対応）、ＥＳＴ（Estimation：推定値）はスパーススペクトル再構成を行ったスペクトラムのピーク周波数である。図６（ｃ）ではＥＳＴは約１．２Ｈｚであり、即ち心臓が１分間に７２回拍動することである（７２ＢＰＭ）。図６（ｃ）の例では、この推定値（７２ＢＰＭ）とＧＴとのずれ（絶対値誤差）は、０．７８ＢＰＭである。
【0033】
【0034】
ここで、ＺＡ－ＬＭＳアルゴリズムについて詳しく説明する。図７は、ＳＳＲに用いる適応フィルタについて説明する図である。図７に示すように、未知システム５０１に入った入力信号ｘ（ｎ）にインパルス応答ｈｏｒｉをかけ、内積ｘＴ（ｎ）ｈｏｒｉを得る。ｘＴ（ｎ）ｈｏｒｉは、付加的な雑音ｚ（ｎ）に干渉されて目的信号ｄ（ｎ）＝ｘＴ（ｎ）ｈｏｒｉ＋ｚ（ｎ）となり、適応フィルタ５０２に入力される。
【0035】
目的信号ｄ（ｎ）から、入力信号ｘ（ｎ）とフィルタ係数ｈ（ｎ）との内積ｘＴ（ｎ）ｈ（ｎ）を減算することで、回帰誤差ｅ（ｎ）が得られる。適応アルゴリズムを用い、回帰誤差ｅ（ｎ）を最小化するように確率的勾配降下法に基づいてフィルタ係数ｈ（ｎ）の回帰更新を行ってｈｏｒｉを推定することができる。
【0036】
ＳＳＲ問題とは、ドップラー信号に対する少数のサンプリング信号ｙにより、フーリエ観測行列Φを用いて高解像度の心拍スペクトルを復元することである。本実施の形態では、ＳＳＲ問題と適応フィルタ５０２のパラメータとの対応関係に基づき、図７に示す適応フィルタ５０２を用い、図８に示すＳＳＲ問題への応用を行う。
【0037】
具体的には、図８に示すフーリエ観測ベクトルφｊを入力信号ｘ（ｎ）、スペクトル未知解ｓ（ｎ）をフィルタ係数ｈ（ｎ）、サンプリング信号の要素ｙｊを目的信号ｄ（ｎ）に対応させる。なお、図８において、５０３は時間領域のドップラー信号を模式的に示したものであり、５０４は５０３をフーリエ変換した周波数領域のドップラー信号を模式的に示したものである。又、ＨＲ（Heart Rate）は心拍数に対応するピークを模式的に示したものである。
【0038】
ＺＡ－ＬＭＳアルゴリズムの場合、インパルス応答ｈｏｒｉが２変数を含むことを仮定する。回帰毎に一つの入力信号ｘ（ｎ）を確率的に用い、２次の回帰誤差ｅ（ｎ）による目的関数の勾配∇ｆを逆方向にし、ステップサイズμをかけ、係数修正を行う。係数修正により、フィルタ係数ｈ（ｎ）がｈｏｒｉに少しずつ近づいていく。
【0039】
以上がＺＡ－ＬＭＳアルゴリズムの説明であるが、非ガウス雑音環境では、回帰誤差が不安定に更新する可能性が高く、ＺＡ－ＬＭＳアルゴリズムでは理想的な復号結果を達成することが困難な場合がある。そこで、インパルス雑音に有効に対応するため、図９に示すように、回帰誤差を一定の範囲内に制限し、突然のインパルス外乱に対しても誤差が急増しないことを実現するＺＡ－ＳＬＭＳ（Zero Attracting Sign Least Mean Square）アルゴリズムを適用することがより好ましい。
【0040】
図９において、破線の円及び実線の矢印は回帰誤差を制限した勾配降下を示し、一点鎖線の矢印はインパルス雑音に干渉された勾配降下を示している。又、Ｓ（ｎ）はスペクトル未知解を、Ｓｏｒｉは心拍スペクトルを示している。このように、ＺＡ－ＬＭＳアルゴリズムにスパース制約を導入することにより、スペクトル未知解Ｓ（ｎ）は、より実際のスパースな心拍スペクトルＳｏｒｉに近づくことができる。
【0041】
ＺＡ－ＬＭＳの目的関数は、式（１）で示すことができる。
【0042】
【数1】
【0043】
【数2】
【0044】
【数3】
【0045】
【数4】
【0046】
【数5】
【0047】
ピーク追跡部２０５は、選択操作及び確認操作を行う。具体的には、ピーク追跡部２０５は、前回の推定値を中心とした所定の探索範囲内（例えば、前後１０ＢＰＭの範囲内）で、スパーススペクトルがピークとなる周波数を今回の推定値として選択する（選択操作）。次に、ピーク追跡部２０５は、隣接する推定値は大きく変化しないため、推定値の変化が所定の閾値（例えば、閾値θ＝６ＢＰＭ）を超える場合は、前回の推定値を出力する（確認操作）。
【0048】
例えば、図６（ｃ）の信号波形を例にとると、前回の推定値が７２ＢＰＭであった場合、ピーク追跡部２０５は、例えば、６２ＢＰＭ～８２ＢＰＭの範囲内で今回の推定値を選択する。そして、例えば、推定値の変化が６ＢＰＭ以下であれば今回の推定値を出力し、推定値の変化が６ＢＰＭを超える場合には前回の推定値を出力する。
【0049】
以上の各ステップにより、雑音を除去し、心拍検出精度を向上することができる。なお、Ｉ信号のスペクトルピークとＱ信号のスペクトルピークのうち、ピーク値が高い方のデータを採用することができる。
【0050】
このように、心拍検出システム１では、非接触センサであるドップラーセンサで取得したＩ信号及びＱ信号に対して、スパーススペクトル再構成を行う。この際、スパーススペクトル再構成に、確率的勾配降下法に属するアルゴリズムを用いることにより、呼吸や体動による雑音を高精度で除去することができる。その結果、着座静止時のみならず、タイピング時や運転時等の様々な体動を伴う動作時においても、心拍検出精度を向上することができる。ここで、心拍検出は、心拍数の検出や心拍間隔（ＲＲＩ：R-R Interval）の検出を含むものである。
【0051】
なお、確率的勾配降下法に属するアルゴリズムには、雑音に高いロバスト性を持つＺＡ－ＬＭＳアルゴリズムを用いることが好ましく、雑音に更に高いロバスト性を持つＺＡ－ＳＬＭＳアルゴリズムを用いることがより好ましい。但し、これらには限定されず、以下の第１の実施の形態の変形例で示すように、ＦＯＣＵＳＳ（Focal Underdetermined System Solver）アルゴリズム等を用いても構わない。
【0052】
又、スパーススペクトル再構成を行う前に、特異スペクトル分析を行って呼吸や体動による雑音主成分を除去し、心拍成分を抽出すると、心拍検出精度をより向上することができる。又、特異スペクトル分析の後に、時間差分処理を行って再構成するスペクトルのピークを実際の周期的な心拍のピークに近づかせることで、心拍検出精度を更に向上することができる。
【0053】
〈第１の実施の形態の変形例〉
第１の実施の形態の変形例では、図５のステップＳ１６のスパーススペクトル再構成に、第１の実施の形態とは異なるアルゴリズムを適用する例を示す。なお、第１の実施の形態の変形例において、既に説明した実施の形態と同一構成部についての説明は省略する場合がある。
第１の実施の形態の変形例では、図５のステップＳ１６のスパーススペクトル再構成に、第１の実施の形態とは異なるアルゴリズムを適用する例を示す。なお、第１の実施の形態の変形例において、既に説明した実施の形態と同一構成部についての説明は省略する場合がある。
【0054】
第１の実施の形態の変形例では、図５のステップＳ１６のスパーススペクトル再構成に、確率的勾配降下法に属するＦＯＣＵＳＳアルゴリズムを適用し、呼吸や体動による雑音を高精度で除去する。
【0055】
ＦＯＣＵＳＳアルゴリズムは、近似凸最適化に属する非凸ノルムアルゴリズムである。ここでは、圧縮センシング（ＣＳ）技術に基づき、Ｍ×Ｎフーリエ観測行列Φ∈ＣＭ×Ｎを用いたＳＳＲアルゴリズムを適用する。Φは式（６）で示すことができる。
【0056】
【数6】
【0057】
ＳＳＲ問題は式（７）の方程式で表示できる。
【0058】
【数7】
【0059】
【数8】
【0060】
次に、加重最小ノルム解を求める。最小ノルム解と同じ原理に基づき、||q||（ｑ＝Ｗ＋ｘ、ＷはＮ×Ｎ正方行列）を最適化目標として最小にし、式（９）を再定義する。
【0061】
【数9】
【0062】
【数10】
【0063】
【数11】
【0064】
【数12】
【0065】
【数13】
【0066】
【数14】
【0067】
【数15】
【0068】
【数16】
【0069】
【数17】
【0070】
【数18】
【0071】
【数19】
【0072】
【数20】
【0073】
このように、スパーススペクトル再構成に用いる確率的勾配降下法に属するアルゴリズムとして、ＺＡ－ＬＭＳアルゴリズムやＺＡ－ＳＬＭＳアルゴリズムに代えて、ＦＯＣＵＳＳアルゴリズムを用いても、呼吸や体動による雑音を高精度で除去することができる。
【0074】
但し、ＦＯＣＵＳＳアルゴリズムは、復元結果が雑音の影響を受けやすいため、雑音に高いロバスト性を持つＺＡ－ＬＭＳアルゴリズムやＺＡ－ＳＬＭＳアルゴリズムを用いる方が、呼吸や体動による雑音をより高精度で除去することができる。
【0075】
［実施例］
実施例では、心拍検出システム１を用いて心拍検出の実験を行った。実験の諸元を表１に示す。
実施例では、心拍検出システム１を用いて心拍検出の実験を行った。実験の諸元を表１に示す。
【0076】
【表1】
【0077】
図１０～図１３に、代表的な測定結果を示す。図１０は、５人の被験者のうち被験者３が着座静止している場合の心拍数の評価結果である。図１１は、５人の被験者のうち被験者２がタイピングしている場合の心拍数の評価結果である。図１２は、５人の被験者のうち被験者３が着座静止している場合の心拍間隔の評価結果である。図１３は、５人の被験者のうち被験者２がタイピングしている場合の心拍間隔の評価結果である。
【0078】
なお、図１０及び図１１において、ＲＭ－ＦＯＣＵＳＳアルゴリズム、ＺＡ－ＬＭＳアルゴリズム、及びＺＡ－ＳＬＭＳアルゴリズムを用いて測定した各データは、タイムウィンドウを８秒、ずらし時間を２秒として、図６（ａ）～図６（ｃ）を参照して説明した方法によりスパーススペクトル再構成を行ったスペクトラムのピークの周波数（ＥＳＴ）から推定した心拍数をタイムウィンドウごとにプロットしたものである。又、図１０（着席静止時）及び図１１（タイピング時）に示したデータに基づいて心拍数、図１２（着席静止時）及び図１３（タイピング時）に示したデータに基づいて心拍間隔をプロットしたものである。
【0079】
図１４は、５人の被験者について着座静止している場合の心拍数の評価結果をまとめた図であり、基準値に対する絶対値誤差を比較している（単位はＢＰＭ）。図１４に示すように、各アルゴリズムについて、着座静止している場合の心拍数の基準値に対する絶対値誤差の平均値を比較すると、比較例であるスペクトログラム法に比べて、実施例である３つのアルゴリズムは、絶対値誤差の平均値が大幅に低減されており、心拍検出精度を向上できることが確認された。中でも、ＺＡ－ＳＬＭＳアルゴリズムを用いた場合に、絶対値誤差の平均値が最も小さくなっており、ＺＡ－ＳＬＭＳアルゴリズムは心拍検出精度を向上するための好適なアルゴリズムであることが確認された。
【0080】
図１５は、５人の被験者について着座静止している場合の心拍間隔の評価結果をまとめた図であり、ＲＭＳＥ（Root Mean Square Error）を比較している（単位はｍｓｅｃ）。図１５に示すように、各アルゴリズムについて、着座静止している場合の心拍間隔のＲＭＳＥの平均値を比較すると、比較例であるスペクトログラム法に比べて、実施例である３つのアルゴリズムは、ＲＭＳＥの平均値が大幅に低減されており、心拍検出精度を向上できることが確認された。中でも、ＺＡ－ＬＭＳ及びＺＡ－ＳＬＭＳアルゴリズムを用いた場合に、ＲＭＳＥの平均値が最も小さくなっており、ＺＡ－ＬＭＳ及びＺＡ－ＳＬＭＳアルゴリズムは心拍検出精度を向上するための好適なアルゴリズムであることが確認された。
【0081】
図１６は、５人の被験者についてタイピングしている場合の心拍数の評価結果をまとめた図であり、基準値に対する絶対値誤差を比較している（単位はＢＰＭ）。図１６に示すように、各アルゴリズムについて、タイピングしている場合の心拍数の基準値に対する絶対値誤差の平均値を比較すると、比較例であるスペクトログラム法に比べて、実施例である３つのアルゴリズムは、絶対値誤差の平均値が大幅に低減されており、心拍検出精度を向上できることが確認された。中でも、ＺＡ－ＳＬＭＳアルゴリズムを用いた場合に、絶対値誤差の平均値が最も小さくなっており、ＺＡ－ＳＬＭＳアルゴリズムは心拍検出精度を向上するための好適なアルゴリズムであることが確認された。
【0082】
図１７は、５人の被験者についてタイピングしている場合の心拍間隔の評価結果をまとめた図であり、ＲＭＳＥを比較している（単位はｍｓｅｃ）。図１７に示すように、各アルゴリズムについて、タイピングしている場合の心拍間隔のＲＭＳＥの平均値を比較すると、比較例であるスペクトログラム法に比べて、実施例である３つのアルゴリズムは、ＲＭＳＥの平均値が大幅に低減されており、心拍検出精度を向上できることが確認された。中でも、ＺＡ－ＳＬＭＳアルゴリズムを用いた場合に、ＲＭＳＥの平均値が最も小さくなっており、ＺＡ－ＳＬＭＳアルゴリズムは心拍検出精度を向上するための好適なアルゴリズムであることが確認された。
【0083】
以上、好ましい実施の形態等について詳説したが、上述した実施の形態等に制限されることはなく、特許請求の範囲に記載された範囲を逸脱することなく、上述した実施の形態等に種々の変形及び置換を加えることができる。
【符号の説明】
【0084】
１ 心拍検出システム
１０ ドップラーセンサ
２０ 信号処理部
２１ ＣＰＵ
２２ ＲＯＭ
２３ ＲＡＭ
２４ Ｉ／Ｆ
２５ バスライン
２０１ 信号取得部
２０２ 特異スペクトル分析部
２０３ 時間差分処理部
２０４ スパーススペクトル再構成部
２０５ ピーク追跡部
１０ ドップラーセンサ
２０ 信号処理部
２１ ＣＰＵ
２２ ＲＯＭ
２３ ＲＡＭ
２４ Ｉ／Ｆ
２５ バスライン
２０１ 信号取得部
２０２ 特異スペクトル分析部
２０３ 時間差分処理部
２０４ スパーススペクトル再構成部
２０５ ピーク追跡部
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】