(57)【特許請求の範囲】
【請求項１】
ビート信号の二次元平面における分布に基づいて前記ビート信号を一次元の第１候補信号に変換する第１変換部と、
前記ビート信号の二次元平面における位置の時間変化に基づいて前記ビート信号を一次元の第２候補信号に変換する第２変換部と、
前記第１変換部から出力された前記第１候補信号および前記第２変換部から出力された前記第２候補信号に基づいて前記第１候補信号または前記第２候補信号の少なくともいずれかを選択することにより一次元信号を決定する信号決定部と、
を備え、
前記ビート信号は、物体の運動に応じてドップラーセンサにより検出され、
前記二次元平面は、ＩＱ平面である、振動状態推定装置。
【請求項２】
前記信号決定部は、前記第１候補信号の特徴量と前記第２候補信号の特徴量とを比較することにより決定される評価値に基づいて前記第１候補信号または前記第２候補信号の少なくともいずれかを選択する、請求項１に記載の振動状態推定装置。
【請求項３】
前記信号決定部は、前記第１候補信号の周波数と前記第２候補信号の周波数とを比較することにより上記評価値を決定する、請求項２に記載の振動状態推定装置。
【請求項４】
前記信号決定部は、変換前の前記ビート信号の二次元平面における分布に基づいて推定される前記分布の中心を基準とする前記ビート信号の回転角度の変化量に応じて決定される評価値に基づいて前記第１候補信号または前記第２候補信号の少なくともいずれかを選択する、請求項１に記載の振動状態推定装置。
【請求項５】
前記信号決定部は、前記第１候補信号および前記第２候補信号の双方を選択する場合、前記第１候補信号および前記第２候補信号を、前記評価値を用いて重み付けすることにより、前記一次元信号を決定する、請求項２〜４のいずれか１項に記載の振動状態推定装置。
【請求項６】
前記信号決定部は、前記第１候補信号もしくは前記第２候補信号のうち一時刻前に選択された一方の候補信号と異なる候補信号を選択する場合、または、前記第１候補信号および前記第２候補信号の双方を選択する場合、前記第１候補信号または前記第２候補信号のうち一方の候補信号の位相を、他方の候補信号の位相に整合させる、請求項１〜５のいずれか１項に記載の振動状態推定装置。
【請求項７】
前記第１変換部は、前記ビート信号を表現する二次元ベクトルと、前記ビート信号の共分散行列に対する最大固有値に対応した固有ベクトルとの内積を算出することにより、前記ビート信号を前記第１候補信号に変換する、請求項１〜６のいずれか１項に記載の振動状態推定装置。
【請求項８】
前記第２変換部は、前記ビート信号の二次元平面における分布に基づいて推定される前記分布の中心からの前記ビート信号の位置までの距離と、前記中心を基準とする前記ビート信号の回転角度の変化量との積を算出することにより、前記ビート信号を前記第２候補信号に変換する、請求項１〜７のいずれか１項に記載の振動状態推定装置。
【請求項９】
前記信号決定部により決定された前記一次元信号の周波数を推定する周波数推定部と、
前記周波数推定部により推定された前記周波数に基づいて前記一次元信号の基準位置を推定する基準位置推定部と、
前記基準位置推定部により推定された一の基準位置と連続する他の基準位置との間隔を前記一次元信号の周期として算出する周期算出部と、
をさらに備える、請求項１〜８のいずれか１項に記載の振動状態推定装置。
【請求項１０】
前記ビート信号は、生体の呼吸に伴う振動に応じて前記ドップラーセンサにより検出されたビート信号である、請求項１〜９のいずれか１項に記載の振動状態推定装置。
【請求項１１】
ビート信号の二次元平面における分布に基づいて前記ビート信号を一次元の第１候補信号に変換するステップと、
前記ビート信号の二次元平面における位置の時間変化に基づいて前記ビート信号を一次元の第２候補信号に変換するステップと、
変換された前記第１候補信号および変換された前記第２候補信号に基づいて前記第１候補信号または前記第２候補信号の少なくともいずれかを選択することにより一次元信号を決定するステップと、
を含み、
前記ビート信号は、物体の運動に応じてドップラーセンサにより検出され、
前記二次元平面は、ＩＱ平面である、振動状態推定方法。
【請求項１２】
コンピュータを、
ビート信号の二次元平面における分布に基づいて前記ビート信号を一次元の第１候補信号に変換する第１変換部と、
前記ビート信号の二次元平面における位置の時間変化に基づいて前記ビート信号を一次元の第２候補信号に変換する第２変換部と、
前記第１変換部から出力された前記第１候補信号および前記第２変換部から出力された前記第２候補信号に基づいて前記第１候補信号または前記第２候補信号の少なくともいずれかを選択することにより一次元信号を決定する信号決定部と、
として機能させるためのプログラムであり、
前記ビート信号は、物体の運動に応じてドップラーセンサにより検出され、
前記二次元平面は、ＩＱ平面である、プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、振動状態推定装置、振動状態推定方法、およびプログラムに関する。
【背景技術】
【０００２】
近年、被験者の健康状態を判定するために、被験者の呼吸などの振動を非接触で検出し、検出された振動状態を推定する技術が開発されている。例えば、下記特許文献１には、被験者が使用する寝具に備えつけられたエアパッドを介して被験者の振動を計測し、計測された信号から被験者の振動の周期を算出する手段について開示されている。また、下記特許文献２には、ドップラーセンサにより取得された被験者の振動を示すＩＱ信号（Ｉ：Ｉｎ−ｐｈａｓｅ、Ｑ：Ｑｕａｄｒａｔｕｒｅ−ｐｈａｓｅ）について、ＩＱ平面と呼ばれる二次元平面上における信号のノルムから振動の１周期ごとの波形の境界を特定することにより振動の周期を算出する技術が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００６−２６３０３２号公報
【特許文献２】特開２０１４−０１４７０８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかし、上記特許文献１においては、振動波形に含まれる微細な揺らぎをローパスフィルタにより減少させているため、微細な揺らぎを含む振動の波形は正確に抽出されておらず、推定される振動の周期の精度は高くない。また、上記特許文献２においては、振動を示すＩＱ信号にノイズが重畳している場合、または、振動が連続的で切れ目なく生じている場合、振動の周期の区切りとなる振動の停止区間が検出されにくくなるため、振動の周期を推定することが困難となる。
【０００５】
そこで、本発明は、上記問題に鑑みてなされたものであり、本発明の目的とするところは、振動の周期を推定するための振動の波形を精度高く抽出することが可能な、新規かつ改良された振動状態推定装置を提供することにある。
【課題を解決するための手段】
【０００６】
上記課題を解決するために、本発明のある観点によれば、ビート信号の二次元平面における分布に基づいて上記ビート信号を一次元の第１候補信号に変換する第１変換部と、上記ビート信号の二次元平面における位置の時間変化に基づいて上記ビート信号を一次元の第２候補信号に変換する第２変換部と、上記第１変換部から出力された上記第１候補信号および上記第２変換部から出力された上記第２候補信号に基づいて上記第１候補信号または上記第２候補信号の少なくともいずれかを選択することにより一次元信号を決定する信号決定部と、を備え、上記ビート信号は、物体の運動に応じてドップラーセンサにより検出され、上記二次元平面は、ＩＱ平面である、振動状態推定装置が提供される。
【０００７】
上記信号決定部は、上記第１候補信号の特徴量と上記第２候補信号の特徴量とを比較することにより決定される評価値に基づいて上記第１候補信号または上記第２候補信号の少なくともいずれかを選択してもよい。
【０００８】
上記信号決定部は、上記第１候補信号の周波数と上記第２候補信号の周波数とを比較することにより上記評価値を決定してもよい。
【０００９】
上記信号決定部は、変換前の上記ビート信号の二次元平面における分布に基づいて推定される上記分布の中心を基準とする上記ビート信号の回転角度の変化量に応じて決定される評価値に基づいて上記第１候補信号または上記第２候補信号の少なくともいずれかを選択してもよい。
【００１０】
上記信号決定部は、上記第１候補信号および上記第２候補信号の双方を選択する場合、上記第１候補信号および上記第２候補信号を、上記評価値を用いて重み付けすることにより、上記一次元信号を決定してもよい。
【００１１】
上記信号決定部は、上記第１候補信号もしくは上記第２候補信号のうち一時刻前に選択された一方の候補信号と異なる候補信号を選択する場合、または、上記第１候補信号および上記第２候補信号の双方を選択する場合、上記第１候補信号または上記第２候補信号のうち一方の候補信号の位相を、他方の候補信号の位相に整合させてもよい。
【００１２】
上記第１変換部は、上記ビート信号を表現する二次元ベクトルと、上記ビート信号の共分散行列に対する最大固有値に対応した固有ベクトルとの内積を算出することにより、上記ビート信号を上記第１候補信号に変換してもよい。
【００１３】
上記第２変換部は、上記ビート信号の二次元平面における分布に基づいて推定される上記分布の中心からの上記ビート信号の位置までの距離と、上記中心を基準とする上記ビート信号の回転角度の変化量との積を算出することにより、上記ビート信号を上記第２候補信号に変換してもよい。
【００１４】
上記振動状態推定装置は、上記信号決定部により決定された上記一次元信号の周波数を推定する周波数推定部と、上記周波数推定部により推定された上記周波数に基づいて上記一次元信号の基準位置を推定する基準位置推定部と、上記基準位置推定部により推定された一の基準位置と連続する他の基準位置との間隔を上記一次元信号の周期として算出する周期算出部と、をさらに備えてもよい。
【００１５】
上記ビート信号は、生体の呼吸に伴う振動に応じて上記ドップラーセンサにより検出されたビート信号であってもよい。
【００１６】
また、上記課題を解決するために、本発明の別の観点によれば、ビート信号の二次元平面における分布に基づいて上記ビート信号を一次元の第１候補信号に変換するステップと、上記ビート信号の二次元平面における位置の時間変化に基づいて上記ビート信号を一次元の第２候補信号に変換するステップと、変換された上記第１候補信号および変換された上記第２候補信号に基づいて上記第１候補信号または上記第２候補信号の少なくともいずれかを選択することにより一次元信号を決定するステップと、を含み、上記ビート信号は、物体の運動に応じてドップラーセンサにより検出され、上記二次元平面は、ＩＱ平面である、振動状態推定方法が提供される。
【００１７】
また、上記課題を解決するために、本発明の別の観点によれば、コンピュータを、ビート信号の二次元平面における分布に基づいて上記ビート信号を一次元の第１候補信号に変換する第１変換部と、上記ビート信号の二次元平面における位置の時間変化に基づいて上記ビート信号を一次元の第２候補信号に変換する第２変換部と、上記第１変換部から出力された上記第１候補信号および上記第２変換部から出力された上記第２候補信号に基づいて上記第１候補信号または上記第２候補信号の少なくともいずれかを選択することにより一次元信号を決定する信号決定部と、として機能させるためのプログラムであり、上記ビート信号は、物体の運動に応じてドップラーセンサにより検出され、上記二次元平面は、ＩＱ平面である、プログラムが提供される。
【発明の効果】
【００１８】
以上説明したように本発明によれば、振動の周期を推定するための振動の波形を精度高く抽出することが可能である。
【図面の簡単な説明】
【００１９】
【図１】本発明の一実施形態に係る振動状態推定システムの概要を示す図である。
【図２】本発明の一実施形態に係る振動状態推定装置の構成例を示すブロック図である。
【図３】本発明の一実施形態に係る信号変換部の構成例を示すブロック図である。
【図４】ビート信号のＩＱ平面における軌跡の一例を示す図である。
【図５】ビート信号のＩＱ平面における軌跡の一例を示す図である。
【図６】本発明の一実施形態に係る第１変換部により出力された第１候補信号の一例を示す図である。
【図７】ビート信号のＩＱ平面における軌跡について本発明の一実施形態に係る第２変換部による変換処理を実施する例を示す図である。
【図８】本発明の一実施形態に係る信号決定部における一次元信号の決定処理フローの一例を示す図である。
【図９】本発明の一実施形態に係る信号決定部により算出された指標値の分布を示す図である。
【図１０】本発明の一実施形態に係る振動状態推定装置の動作例を示すフローチャートである。
【図１１】本発明の実施形態に係る振動状態推定装置のハードウェア構成例を示すブロック図である。
【発明を実施するための形態】
【００２０】
以下に添付図面を参照しながら、本発明の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。
【００２１】
＜１．振動状態推定システムの構成例＞
図１は、本発明の一実施形態に係る振動状態推定システム１の概要を示す図である。図１を参照すると、振動状態推定システム１は、ドップラーセンサ２、および振動状態推定装置１０を有する。
【００２２】
ドップラーセンサ２は、図１に示すように、例えば部屋の天井等に設置され、検出エリアである部屋の屋内に向けて、光、電磁波、または音波等の放射波を放射し、物体（例えば図１においては人物Ｐ）により反射された反射波を受信する。このとき、反射波の周波数は、物体の振動等の運動により生じるドップラー効果により、放射波の周波数から変化する。ドップラーセンサ２は、放射波の周波数と反射波の周波数との差分の周波数を有するビート信号を生成する。ドップラーセンサ２は直交検波方式を搭載してもよく、その場合、ドップラーセンサ２は、余弦波成分（Ｉ成分）および正弦波成分（Ｑ成分）の２種類のビート信号を生成する。ドップラーセンサ２は、生成したビート信号を振動状態推定装置１０に出力する。
【００２３】
なお、ドップラーセンサ２は、物体の振動が検出可能であれば任意の位置に設置可能である。また、図１に示した例では、ドップラーセンサ２は放射波の送信部と反射波の受信部とが一体となって構成されているが、ドップラーセンサ２は、送信部と受信部とを分離した構成により実現されてもよい。また、ドップラーセンサ２が放射する放射波は、物体の振動によりドップラー効果を生じさせることが可能であれば、任意の周波数帯の波であってもよい。例えば、呼吸の振動をドップラーセンサが検知する場合、２．４ＧＨｚ〜１０．５ＧＨｚ帯などの、呼吸の振動の変位と比較して波長が大きい波が用いられることが多い。しかし、本実施形態においては、呼吸の振動の幅広い変位を対象として、かつ波形に含まれる微細なゆらぎを検出するために、２４ＧＨｚ帯などの準ミリ波やミリ波の帯域の波が用いられることが好ましい。
【００２４】
振動状態推定装置１０は、ドップラーセンサ２から出力されたビート信号から、被験者（例えば人物Ｐ）の振動状態を推定する装置である。振動状態推定装置１０は、例えば、ネットワーク上の一または複数の情報処理装置によって実現されてもよい。より具体的には、振動状態推定装置１０は、サーバやＰＣ（Ｐｅｒｓｏｎａｌ  Ｃｏｍｐｕｔｅｒ）等により実現されてもよい。ここで、本明細書において振動状態とは、例えば物体の振動波形から得られる振幅や周期を意味する。例えば、振動状態推定装置１０は、人物Ｐの呼吸等の振動から得られるビート信号から、人物Ｐの振動状態を推定することが可能である。振動状態推定装置１０は、例えば、不図示の通信部を介して、推定した振動状態を有線または無線により外部装置等に送信することが可能である。これにより、例えば、上記外部装置において、振動状態について解析を行うことが可能となる。振動状態の解析により、例えば、人物Ｐの健康状態等を把握することが可能となる。なお、振動状態の解析については、上記のような外部装置ではなく、振動状態推定装置１０本体において行うことも可能である。
【００２５】
ここで、人物Ｐの振動から得られるビート信号には、人物Ｐの様々な運動に起因する信号が含まれ得る。例えば、人物Ｐの呼吸の振動のみならず、人物Ｐの体動や心拍などの振動もビート信号に含まれ得る。さらに、ドップラーセンサ２が生成したビート信号には、周囲の環境に起因するノイズ等も含まれ得る。そのため、例えば呼吸の振動を検出対象とする場合、呼吸の波形情報のみをビート信号から検出することが望ましい。
【００２６】
例えば、特開２００６−２６３０３２号公報においては、呼吸の周期を検出するために、振動信号に対して、呼吸の周波数帯に相当する０．１５〜０．３０Ｈｚのバンドパスフィルタによるフィルタ処理が施されている。しかしその場合、呼吸の１周期ごとに微細な揺らぎが生じている場合において、その揺らぎがバンドパスフィルタにより低減されてしまうため、微細な揺らぎを伴う振動を表現する波形を抽出することが困難となる。また、上記公報に開示されている技術をビート信号に対して適用させる場合、振動の振幅の大きさによっては、詳細は後述するが、振動波形に高周波成分が重畳されてしまう。この場合、一回の呼吸において複数のピークを有する波形が抽出されるため、呼吸による振動の周期を誤って推定してしまう可能性がある。
【００２７】
また、特開２０１４−０１４７０８号公報においては、ＩＱ平面上におけるビート信号のノルムを算出し、ノルムが略ゼロとなる時刻を呼吸の吸い込みと吐き出しの切り替わり時として検出することにより、呼吸の周期が推定されている。しかし、例えばビート信号にノイズが多く重畳している場合や、呼吸の動作が切れ目なく連続的である場合、ノルムが略ゼロとなる時刻の検出が困難となる。そのため、呼吸の周期を誤って推定してしまう可能性がある。
【００２８】
そこで、上記事情を一着眼点にして、振動状態推定装置１０を創作するに至った。本発明の一実施形態に係る振動状態推定装置１０は、例えば、呼吸などの振動の周期を推定するための振動の波形を精度高く抽出することが可能である。以下、一実施形態に係る振動状態推定装置１０の構成についてより詳細に説明する。
【００２９】
＜２．振動状態推定装置の構成例＞
図２は、本発明の一実施形態に係る振動状態推定装置１０の構成例を示すブロック図である。図２を参照すると、振動状態推定装置１０は、ビート信号取得部１１０、フィルタ部１２０、信号変換部１３０、周波数推定部１４０、基準位置推定部１５０、および周期算出部１６０を含む。図示された振動状態推定装置１０は、取得したビート信号を一次元信号に変換し、変換した一次元信号の周波数を推定し、推定された周波数からビート信号の周期を算出するための一次元信号における基準位置を推定する。これにより、微細な揺らぎを含む呼吸の振動を抽出し、さらに微細な揺らぎを含む呼吸の振動の周期を一回の呼吸ごとに推定することが可能となる。これにより、例えば推定された振動の周期のばらつきから、被験者のストレス状態などを検知することが可能である。以下、各構成の機能について説明する。
【００３０】
（ビート信号取得部）
ビート信号取得部１１０は、ドップラーセンサ２から出力されるビート信号Ｄ（ｔ）を取得する。ドップラーセンサ２から出力されるビート信号Ｄ（ｔ）はＩ成分とＱ成分の２波の成分を有し、振幅をＡ（ｔ）、波長をλ、時刻ｔにおけるドップラーセンサ２と対象物体（例えば図１に示された人物Ｐ）との距離をｄ（ｔ）、初期位相をφ０、直流成分をＯ、ノイズ成分をｗとすると、ビート信号Ｄ（ｔ）は、下記の数式１のように表現される。
【００３１】
【数１】
0006477199.tif 000002
【００３２】
ビート信号取得部１１０は、取得したビート信号Ｄ（ｔ）をフィルタ部１２０に出力する。なお、詳しくは後述するが、他の実施形態において振動状態推定装置がフィルタ部１２０を有さない場合は、ビート信号取得部１１０はビート信号Ｄ（ｔ）を信号変換部１３０に出力してもよい。また、ビート信号取得部１１０は、取得したビート信号Ｄ（ｔ）を時系列に不図示の記憶部に記憶してもよい。
【００３３】
（フィルタ部）
フィルタ部１２０は、ビート信号取得部１１０が取得したビート信号Ｄ（ｔ）に含まれる直流成分Ｏなどの低周波成分や、ノイズ成分ｗなどを減少または除去するフィルタ処理を実施し、フィルタ処理されたビート信号Ｄ（ｔ）を信号変換部１３０に出力する。これにより、例えば、ビート信号Ｄ（ｔ）の振幅が微弱である場合において、ビート信号Ｄ（ｔ）を不図示の増幅器等により増幅させることが可能となる。フィルタ部１２０は、例えば、ローパスフィルタ、ハイパスフィルタ、バンドパスフィルタ、ＩＩＲフィルタなどの多様なフィルタを用いることができ、またこれらのフィルタの組み合わせを用いてもよい。また、フィルタ部１２０は、例えば、人物Ｐの体動や心拍など、呼吸とは関係のない振動を含む成分を低減させてもよい。なお、本実施形態ではフィルタ部１２０によりビート信号Ｄ（ｔ）に含まれる直流成分Ｏやノイズ成分ｗ等が低減されるが、他の実施形態においてビート信号Ｄ（ｔ）についてのフィルタ処理が不要であると判断された場合、振動状態推定装置１０は必ずしもフィルタ部１２０を含まなくてもよい。
【００３４】
（信号変換部）
信号変換部１３０は、ビート信号取得部１１０またはフィルタ部１２０から出力されたビート信号Ｄ（ｔ）を一次元信号ｒ（ｔ）に変換する。より具体的には、信号変換部１３０は、少なくとも２つの変換手段によりビート信号Ｄ（ｔ）から変換された一次元の候補信号のうち、少なくともいずれかを用いて決定された一次元信号ｒ（ｔ）を出力する。Ｉ成分とＱ成分を有する二次元信号であるビート信号Ｄ（ｔ）を一次元信号ｒ（ｔ）に変換することにより、呼吸の振動の周期を推定するための処理が容易になる。なお、信号変換部１３０の構成および処理については後述する。
【００３５】
（周波数推定部）
周波数推定部１４０は、信号変換部１３０から出力された一次元信号ｒ（ｔ）の有する周波数ｆＥ（ｔ）を推定する。例えば、周波数推定部１４０は、呼吸周期の大きな揺らぎに合わせるために、一次元信号ｒ（ｔ）と参照信号との位相差、および位相差の時間変化に基づいて一次元信号ｒ（ｔ）の周波数を推定する、第１周波数推定部を有してもよい。また、周波数推定部１４０は、第１周波数推定部により推定された推定周波数に対応する周期に相当する区間長を有する比較信号を一次元信号ｒ（ｔ）から切り出し、比較信号と上記の区間長を有する一次元信号ｒ（ｔ）との相関係数を算出することにより一次元信号ｒ（ｔ）の周波数を推定する、第２周波数推定部を有してもよい。以下、第１周波数推定部により推定された一次元信号ｒ（ｔ）の周波数を第１推定周波数ｆＡ（ｔ）と、また、第２周波数推定部により推定された一次元信号ｒ（ｔ）の周波数を第２推定周波数ｆＢ（ｔ）と呼称する。
【００３６】
−第１周波数推定部
第１周波数推定部は、時刻ｔｋにおける一次元信号ｒ（ｔk）と参照信号との位相差ψ（ｔk）から、第１推定周波数ｆＡ（ｔk）を算出する。第１周波数推定部は、まず、一次元信号ｒ（ｔｋ）と、参照信号とを乗算し、乗算された信号の低周波成分を、ローパスフィルタ等を用いて抽出することにより、一次元信号ｒ（ｔk）と参照信号との位相差ψ（ｔk）を算出する。ここでは参照信号として、位相変化による周期推定を容易とするため、位相とその時間変化をパラメータにもつ正弦波、または余弦波が用いられる。その場合、参照信号は、時刻ｔk−１に推定された第１周波数ｆＡ（ｔk−１）を有する正弦波、あるいは余弦波の信号である。そして、第１周波数推定部は、位相差ψ（ｔk）と、時刻ｔk−１に算出された位相差ψ（ｔk−１）との変化分を算出する。位相差の時間変化分ψ（ｔk）−ψ（ｔk−１）は、一次元信号ｒ（ｔ）の周波数の変化によって生じたものと仮定して、第１周波数推定部は、参照信号の周波数ｆＡ（ｔk−１）に、位相差の時間変化分ψ（ｔk）−ψ（ｔk−１）に相当する周波数を加算して、第１推定周波数ｆＡ（ｔk）を出力する。
【００３７】
なお、上記位相差の時間変化分が所定の閾値以上であった場合、周波数のフィードバック値が大きく見積もられやすいため、第１推定周波数ｆＡ（ｔk）が振動を伴って収束しないといった問題がある。そのため、参照信号の周波数ｆＡ（ｔk−１）に加算する周波数の値は、位相差の時間変化分ψ（ｔk）−ψ（ｔk−１）に相当する周波数に、１未満の係数を乗算した値でもよい。これにより、比較的短い時間で第１推定周波数ｆＡ（ｔｋ）が収束することが可能である。
【００３８】
また、算出された第１推定周波数ｆＡ（ｔk）が所定の周波数帯域に存在しない場合、参照信号の周波数ｆＡ（ｔk−１）を、第１推定周波数ｆＡ（ｔk）として出力してもよい。ここで所定の周波数帯域とは、例えば、本実施形態においては、呼吸運動の周期に相当する周波数帯域を意味する。
【００３９】
−第２周波数推定部
第２周波数推定部は、入力された一次元信号ｒ（ｔ）の相関分析を行うことにより、第１周波数推定部により得られた第１推定周波数ｆＡ（ｔk）より局所的な変化において精度の高い第２推定周波数ｆＢ（ｔk）を推定する。まず、第２周波数推定部は、一次元信号ｒ（ｔ）から、時刻ｔk−ＴＡから時刻ｔkまでの区間を切り出し、これを比較信号として扱う。ＴＡは、先に推定された第１推定周波数ｆＡ（ｔk−１）に対応する周期である。次に、第２周波数推定部は、切り出された比較信号を用いて、相関係数を算出する。相関係数の算出範囲は、ｔk−２ＴＡからｔk−ＴＡの範囲を含むことが好ましい。そして第２周波数推定部は、上記の算出範囲において、区間長ＴＡを有する一次元信号ｒ（ｔ）と、比較信号との相関係数Ｒ（τ）が最大となるタイムラグτｍａｘを算出する。このとき、第２周波数推定部は、τ＝τｍａｘのときの相関係数Ｒ（τｍａｘ）を不図示の記憶部に記憶してもよい。第２周波数推定部は、相関係数算出範囲において算出される相関係数が最大となるときのタイムラグτｍａｘから、第２推定周波数ｆＢ（ｔｋ）を算出する。例えば、第２推定周波数ｆＢ（ｔｋ）は、タイムラグτｍａｘの逆数であってもよい。
【００４０】
−判定部
周波数推定部１４０は、上述した２つの推定された周波数について、いずれの周波数を用いるかを判定する判定部を有してもよい。例えば、判定部は、２つの推定された周波数の差が所定の閾値以内であるかに基づいて、いずれの周波数を推定周波数ｆＥ（ｔ）として用いるかを判定してもよい。また、第２周波数推定部で算出されたタイムラグτｍａｘの値に基づいて、第１推定周波数ｆＡ（ｔk）と第２推定周波数ｆＢ（ｔｋ）のいずれの周波数を用いるかを判定してもよい。
【００４１】
以上より、周波数推定部１４０は、おおよその周波数を推定できる第１周波数推定部と、イレギュラーな変動には追従することが困難であるが、波形の局所的な変化を検出することにより周波数を推定する第２周波数推定部とのどちらを用いるかを判定部により判定することにより、微細な揺らぎを多く含む波形の周波数を、１回の呼吸ごとに高い精度で推定することができる。周波数推定部１４０により推定された推定周波数ｆＥ（ｔ）は、基準位置推定部１５０に出力される。
【００４２】
（基準位置推定部）
基準位置推定部１５０は、推定周波数ｆＥ（ｔ）に基づいて一次元信号ｒ（ｔ）の基準位置を推定する。ここで基準位置とは、本明細書においては、一次元信号ｒ（ｔ）の１周期ごとの波形について、各波形の類似する波形のうち、ピーク位置などの共通する特徴を持った位置を意味する。一次元信号ｒ（ｔ）は時系列信号であるため、基準位置は時刻により特定される。また、基準位置は、例えば波形のボトム位置でもよいし、ｒ（ｔ）＝０となる位置でもよい。
【００４３】
基準位置推定部１５０は、まず、推定時刻ｔｋを終点とする推定周波数ｆＥ（ｔk）に対応する周期ＴＥを区間長とする基準信号を、一次元信号ｒ（ｔ）から切り出す。基準位置推定部１５０は、切り出された基準信号の基準位置を推定する。例えば、基準位置推定部１５０は、基準信号と同一の周期ＴＥを有する所定の初期位相（本実施形態においては、初期位相を０とする）の余弦波または正弦波を生成し、上記余弦波または正弦波と基準信号との位相差を算出してもよい。ここで算出された位相差と時刻ｔｋから、一次元信号ｒ（ｔ）の基準信号区間内における基準位置を推定することが可能である。
【００４４】
なお、基準位置推定部１５０は、推定された基準位置の誤差を補正するために、上記基準信号の区間内において複数回推定された基準位置の分布から、統計手段を用いて基準位置を確定してもよい。例えば、基準位置推定部１５０は、基準位置の分布から得られるモード（最頻値）を、基準位置として確定してもよい。
【００４５】
以上より、基準位置推定部１５０は、一次元信号ｒ（ｔ）の基準位置を推定することにより、微細な揺らぎを含む波形においても、周期を算出するための情報を導出することが可能となる。基準位置推定部１５０は、推定または確定した基準位置を周期算出部１６０に出力する。
【００４６】
（周期算出部）
周期算出部１６０は、基準位置推定部１５０から出力された複数の基準位置から、ビート信号Ｄ（ｔ）の周期を算出する。具体的には、周期算出部１６０は、連続する基準位置の差分を、ビート信号Ｄ（ｔ）の１周期として算出してもよい。これにより、微細な揺らぎを含む波形により表現される呼吸振動の周期を、１回の呼吸ごとに算出することが可能である。
【００４７】
＜３．信号変換部の構成例＞
図３は、本発明の一実施形態に係る信号変換部１３０の構成例を示すブロック図である。図３を参照すると、信号変換部１３０は、第１変換部１３１と、第２変換部１３２と、信号決定部１３３とを含む。以下、それぞれの構成要素について説明する。
【００４８】
（第１変換部）
第１変換部１３１は、ビート信号取得部１１０、またはフィルタ部１２０から取得したビート信号Ｄ（ｔ）について、ビート信号Ｄ（ｔ）のＩＱ平面における分布に基づいて、ビート信号Ｄ（ｔ）を一次元の第１候補信号Ｒ１（ｔ）に変換し、変換した第１候補信号Ｒ１（ｔ）を信号決定部１３３に出力する。例えば、第１変換部１３１は、ビート信号Ｄ（ｔ）を表現する二次元ベクトルと、ビート信号Ｄ（ｔ）の共分散行列の最大固有値に対応する固有ベクトルとの内積を算出することにより、ビート信号Ｄ（ｔ）を第１候補信号Ｒ１（ｔ）に変換してもよい。つまり、第１変換部１３１は、ビート信号Ｄ（ｔ）を最大固有値に対応する主成分方向に射影した第１候補信号Ｒ１（ｔ）を得ることができる。具体的には、第１候補信号Ｒ１（ｔ）は、下記の数式２により算出されてもよい。
【００４９】
【数２】
0006477199.tif 000003
【００５０】
ここで、ｂ（ｔ）はビート信号Ｄ（ｔ）のＩ成分とＱ成分を表現する二次元ベクトル、ｐ（ｔ）は、ビート信号Ｄ（ｔ）の共分散行列の最大固有値に対応する固有ベクトルである。なお、ビート信号Ｄ（ｔ）を表現する二次元ベクトルの原点は、例えばビート信号Ｄ（ｔ）の分布に基づいて推定されてもよい。例えば、過去数秒〜数十秒におけるビート信号Ｄ（ｔ）の軌跡の、ＩＱ平面の原点からの最大位置および最小位置の中点が、原点であってもよい。また、最小二乗法等によりビート信号Ｄ（ｔ）の軌跡を円や楕円として近似し、その近似した図形の中心点が原点であってもよい。また、フィルタ部１２０においてハイパスフィルタやバンドパスフィルタを用いてフィルタ処理されたビート信号はＩＱ平面の原点に近接するので、ＩＱ平面の原点が二次元ベクトルの原点であってもよい。
【００５１】
図４は、ビート信号Ｄ（ｔ）のＩＱ平面における軌跡の一例を示す図である。図４に示されるグラフＧ１を参照すると、ＩＱ平面上においてビート信号Ｄ（ｔ）の軌跡Ｔｒ１が示されている。例えば、呼吸の振動の変位がドップラーセンサ２から放射される放射波の波長の１／１０程度であるなど比較的小さい場合、ビート信号Ｄ（ｔ）は、軌跡Ｔｒ１のような扁平な８の字の分布となることが多い。このような場合、例えば軌跡Ｔｒ１について主成分分析を行った場合、軸Ｐｃ１が示す第１主成分に対応する最大固有値は、他の主成分と比較して大きな値を示す。第１変換部１３１は、ビート信号Ｄ（ｔ）を第１主成分に射影することにより、第１候補信号Ｒ１（ｔ）を得ることができる。
【００５２】
なお、本実施形態に係る第１変換部１３１は、軌跡Ｔｒ１について主成分分析等を行うことにより得られる固有値および固有ベクトルを用いてビート信号Ｄ（ｔ）を第１候補信号Ｒ１（ｔ）に変換する。しかし、例えば、軌跡Ｔｒ１のような扁平な８の字の分布について、上記扁平方向の軸方向を抽出し、上記軸方向における信号の動きを抽出することが可能であれば、第１変換部１３１は主成分分析以外の方法を用いてビート信号Ｄ（ｔ）を第１候補信号Ｒ１（ｔ）に変換してもよい。
【００５３】
（課題）
しかし、呼吸の振動の変位が比較的大きい場合、上述した第１変換部１３１によりビート信号Ｄ（ｔ）から変換された第１候補信号Ｒ１（ｔ）の波形は、呼吸の振動を必ずしも正確に表現するとは限らない。図５は、ビート信号Ｄ（ｔ）のＩＱ平面における軌跡の一例を示す図である。図５に示されるグラフＧ２ａを参照すると、ＩＱ平面上においてビート信号Ｄ（ｔ）の軌跡Ｔｒ２が示されている。例えば、呼吸の振動の変位がドップラーセンサ２から放射される放射波の波長の１／２以上であるなど比較的大きい場合、軌跡Ｔｒ２のように、円弧を描くような分布となることが多い。この軌跡Ｔｒ２により表現されるビート信号Ｄ（ｔ）の分布において主成分分析を行った場合、最大固有値の値は低くなり、主成分方向を適切に定めることが困難である。また、仮に主成分方向を示す軸が適切に得られた場合であっても、第１変換部１３１によって得られた第１候補信号Ｒ１（ｔ）は、本来の呼吸周期よりも短い周期を有する高周波の信号を含む形状となり得る。
【００５４】
図６は、本発明の一実施形態に係る第１変換部１３１により出力された第１候補信号Ｒ１（ｔ）の一例を示す図である。図６を参照すると、第１候補信号Ｒ１（ｔ）の時間変化を示す波Ｗ１が示されている。このとき、波Ｗ１は、高周波成分による歪みＷ１１〜Ｗ１４を含んでいる。これは、ビート信号Ｄ（ｔ）の移動方向が呼吸の切り替わりや揺らぎにより反転したり、ビート信号Ｄ（ｔ）が１回の呼吸の間に円弧を何周も描くように移動することにより、主成分方向に射影される第１候補信号Ｒ１（ｔ）の波Ｗ１が、呼吸の周期よりも短い周期（高周波）の振動を含んでしまうためである。そのため、呼吸の振動の変位が、例えば、ドップラーセンサ２の放射波の波長の１／２以上である場合、第１変換部１３１により変換された第１候補信号Ｒ１（ｔ）の波形は、実際の呼吸の振動を表す波形（図６に示した例では、波Ｗ１が波形Ｗ２１〜Ｗ２４を含む波形）とは異なり、高周波成分による歪みを含んでしまう。その結果、例えば、一次元信号ｒ（ｔ）の波形が高周波の信号を含む場合、波形から振動の周期を判定することは困難となり得る。
【００５５】
そこで、信号変換部１３０は、ビート信号Ｄ（ｔ）の位置の時間変化を反映させてビート信号Ｄ（ｔ）を一次元の第２候補信号Ｒ２（ｔ）に変換する第２変換部１３２をさらに含む。第２変換部１３２は、第１変換部１３１では正確な抽出が困難である変位の大きい呼吸の振動を抽出することができる。具体的には、第２変換部１３２は、ビート信号Ｄ（ｔ）の分布ではなくビート信号Ｄ（ｔ）の逐次的な変化を追いかけることにより一次元の信号を得ることが可能である。
【００５６】
（第２変換部）
第２変換部１３２は、ビート信号取得部１１０、またはフィルタ部１２０から取得したビート信号Ｄ（ｔ）について、ビート信号Ｄ（ｔ）のＩＱ平面における位置の時間変化に基づいて、ビート信号Ｄ（ｔ）を一次元の第２候補信号Ｒ２（ｔ）に変換し、変換した第２候補信号Ｒ２（ｔ）を信号決定部１３３に出力する。例えば、第２変換部１３２は、ビート信号Ｄ（ｔ）のＩＱ平面における分布に基づいて推定される分布の中心からのビート信号Ｄ（ｔ）の位置までの距離（振幅に相当）と、分布の中心を基準とするビート信号Ｄ（ｔ）の回転角度の変化量との積を算出することにより、ビート信号Ｄ（ｔ）を第２候補信号Ｒ２（ｔ）に変換してもよい。例えば、Ｒ２（ｔ）は、下記の数式３における連続関数（本実施形態においては、シグモイド関数）を利用した関数を用いて算出されてもよい。
【００５７】
【数３】
0006477199.tif 000004
【００５８】
ここで、Ａｍｐ（ｔ）はビート信号Ｄ（ｔ）の振幅、θ（ｔ）はビート信号Ｄ（ｔ）のＩＱ平面上における角度であり、θ’（ｔ）は、角度θ（ｔ）の単位時間の角度の変化量、または角度θ（ｔ）の時間微分値である。また、ａ、ｂ、およびｃはいずれも定数であり、自由に設定可能である。
【００５９】
図７は、ビート信号Ｄ（ｔ）のＩＱ平面における軌跡について第２変換部１３２による変換処理を実施する例を示す図である。図７に示されるグラフＧ２ｂを参照すると、軌跡Ｔｒ２上におけるビート信号Ｄ（ｔ）について、振幅Ａｍｐ（ｔ）および角度θ（ｔ）が示されている。振幅Ａｍｐ（ｔ）および角度θ（ｔ）の変化量θ’（ｔ）の積は、図７で示されている領域Ａｖ１の面積速度に相当する。この面積速度は、呼吸の振動の単位時間当たりの大きさに相当する。また、変化量θ’（ｔ）は、呼吸の振動の方向（息の吸い込み、吐き出し）により正負の値を取るため、呼吸の振動の方向を表現できる。そのため、振幅Ａｍｐ（ｔ）と変化量θ’（ｔ）の積により、呼吸運動の状態を的確に表現することが可能となる。
【００６０】
なお、振幅Ａｍｐ（ｔ）は、ビート信号Ｄ（ｔ）のＩＱ平面における分布に基づいて推定される分布の中心からのビート信号Ｄ（ｔ）の位置までの距離に相当し、また、角度θ（ｔ）は、分布の中心を基準とするビート信号Ｄ（ｔ）の回転角度に相当する。ここで、分布の中心とは、例えば図７に示した中心点Ｃ１を指す。ここで、ビート信号Ｄ（ｔ）の分布の中心は、ビート信号Ｄ（ｔ）の分布に基づいて推定されてもよい。例えば、過去数秒〜数十秒におけるビート信号Ｄ（ｔ）の分布について、ＩＱ平面の原点からの最大位置と最小位置との中間位置や、分布の中央に相当する位置が、ビート信号Ｄ（ｔ）の分布の中心であってもよい。また、最小二乗法等によりビート信号Ｄ（ｔ）の軌跡を円や楕円として近似し、その近似した図形の中心点が分布の中心であってもよい。また、フィルタ部１２０においてハイパスフィルタやバンドパスフィルタを用いてフィルタ処理されたビート信号Ｄ（ｔ）はＩＱ平面の原点に近接するように分布するので、ＩＱ平面の原点がビート信号Ｄ（ｔ）の分布の中心であってもよい。
【００６１】
（信号決定部）
信号決定部１３３は、第１変換部１３１から出力された第１候補信号Ｒ１（ｔ）および第２変換部１３２から出力された第２候補信号Ｒ２（ｔ）に基づいて第１候補信号Ｒ１（ｔ）または第２候補信号Ｒ２（ｔ）の少なくともいずれかを選択することにより、周波数推定部１４０に出力する一次元信号ｒ（ｔ）を決定する。例えば、信号決定部１３３は、第１候補信号Ｒ１（ｔ）または第２候補信号Ｒ２（ｔ）の一方を一次元信号ｒ（ｔ）と決定してもよい。また、信号決定部１３３は、第１候補信号Ｒ１（ｔ）と第２候補信号Ｒ２（ｔ）とを統合させた信号を一次元信号ｒ（ｔ）と決定してもよい。
【００６２】
図４に示した呼吸の振動が小さい場合のビート信号Ｄ（ｔ）の軌跡Ｔｒ１においては、分布の中心からの角度θ（ｔ）の変化量θ’（ｔ）は、図５に示したビート信号Ｄ（ｔ）の軌跡Ｔｒ２の場合と比較して、ノイズや量子化誤差の影響により精度が低くなる。つまり、呼吸の振動が小さい場合は、第２候補信号Ｒ２（ｔ）よりも、第１候補信号Ｒ１（ｔ）の方が、呼吸の振動の波形を正確に表現し得る。一方で、上述したように、図５に示した呼吸の振動が大きい場合のビート信号Ｄ（ｔ）の軌跡Ｔｒ２においては、第１候補信号Ｒ１（ｔ）よりも、第２候補信号Ｒ２（ｔ）の方が、呼吸の振動の波形を正確に表現し得る。そのため、信号決定部１３３は、第１候補信号Ｒ１（ｔ）または第２候補信号Ｒ２（ｔ）についてそれぞれ比較し、比較結果に基づいて出力する一次元信号ｒ（ｔ）を決定する。
【００６３】
ここで、信号決定部１３３は、一次元信号ｒ（ｔ）を決定するためのパラメータとして、評価値ｓを決定する。評価値ｓは、呼吸の振動の変位に応じて呼吸の振動を表現する波形の精度が変化する第１候補信号Ｒ１（ｔ）および第２候補信号Ｒ２（ｔ）のどちらを優先的に用いるべきかを評価する値である。つまり、評価値ｓを算出し、評価値ｓに基づいて一次元信号ｒ（ｔ）を決定することにより、呼吸の振動の変位にとらわれず、呼吸の振動がより正確に反映された波形を生成することが可能となる。評価値ｓは、例えば、０〜１の間で決定する値であってもよい。この場合、例えば、評価値ｓが０に近ければ第１候補信号Ｒ１（ｔ）が、評価値ｓが１に近ければ第２候補信号Ｒ２（ｔ）が、他方の候補信号よりも呼吸の振動の波形をより正確に表現し得るということを意味する。また、本実施形態においては、評価値ｓは単一であるが、他の実施形態においてはかかる例に限定されない。例えば、評価値ｓは、第１候補信号Ｒ１（ｔ）および第２候補信号Ｒ２（ｔ）についてそれぞれ算出されてもよく、信号決定部１３３は、各候補信号の評価値ｓに基づいて出力する一次元信号ｒ（ｔ）を決定してもよい。また、評価値ｓの数値範囲についても、０〜１の範囲に限定されず、任意の範囲であってもよい。
【００６４】
以下、信号決定部１３３における一次元信号ｒ（ｔ）の決定処理フローについて、評価値ｓの決定方法を含めて説明する。
【００６５】
＜４．信号決定処理フロー＞
図８は、本発明の一実施形態に係る信号決定部１３３における一次元信号ｒ（ｔ）の決定処理フローの一例を示す図である。本実施形態に係る信号決定部１３３は、まず、評価値ｓを決定するための指標値を算出し、その後、信号決定部１３３は、算出した指標値を用いて、評価値ｓを決定する。その後、信号決定部１３３は、決定された評価値ｓを用いて一次元信号ｒ（ｔ）を決定する。
【００６６】
［４−１．指標値の算出］
まず、信号決定部１３３は、第１候補信号Ｒ１（ｔ）および第２候補信号Ｒ２（ｔ）の少なくともいずれかを選択するための評価値ｓを決定するための指標値を算出する処理を行う（Ｓ２０１）。信号決定部１３３は、例えば、以下の２つの手段の少なくともいずれかを用いて算出する。
【００６７】
（周波数による指標値の算出）
信号決定部１３３は、例えば、第１候補信号Ｒ１（ｔ）の特徴量と第２候補信号Ｒ２（ｔ）の特徴量とを比較することにより決定される評価値ｓに基づいて、少なくともいずれかの候補信号を選択することが可能である。この特徴量とは、例えば、本実施形態においては、各候補信号の周波数であってもよい。つまり、信号決定部１３３は、第１候補信号Ｒ１（ｔ）および第２候補信号Ｒ２（ｔ）のそれぞれについて周波数を算出し、周波数の各々を比較することにより評価値ｓを決定してもよい。上述したように、呼吸の振動の変位がドップラーセンサ２から放射される放射波の波長の１／２以上であるなど比較的大きい場合、図６で示したように、第１候補信号Ｒ１（ｔ）は高周波成分を含むような波形を形成し得る。その場合、第１候補信号Ｒ１（ｔ）の周波数は、第２候補信号Ｒ２（ｔ）の周波数と比較して、高い値を取り得る。一方で、第１候補信号Ｒ１（ｔ）が高周波成分をあまり含まない場合は、各周波数はほぼ同じ値を取り得る。したがって、信号決定部１３３は、各候補信号の周波数を比較することにより、第１候補信号Ｒ１（ｔ）または第２候補信号Ｒ２（ｔ）の一方又は双方を選択するための評価値ｓを決定することが可能である。
【００６８】
ここで、信号決定部１３３は、評価値ｓを決定するための指標値ｘ１（ｔ）を、各候補信号の周波数に基づいて算出してもよい。例えば、信号決定部１３３は、第１候補信号Ｒ１（ｔ）の周波数をｆｒ１（ｔ）、第２候補信号Ｒ２（ｔ）の周波数をｆｒ２（ｔ）として、以下の数式４を用いて、評価値ｓを決定するための指標値ｘ１（ｔ）を算出してもよい。
【００６９】
【数４】
0006477199.tif 000005
【００７０】
ここで、ｋは定数であり、０．０〜１．０程度の値が好ましい。例えばｋ＝１とした場合において、呼吸の振動の変位が比較的大きい場合、第１候補信号Ｒ１（ｔ）に高周波成分が含まれる場合、指標値ｘ１（ｔ）は１．０よりも大きい値を取り得る。一方で、第１候補信号Ｒ１（ｔ）が高周波成分をあまり含まない場合は、指標値ｘ１（ｔ）はほぼ１．０に近い値を取り得る。
【００７１】
なお、各候補信号の周波数は、例えば、各候補信号についてヒルベルト変換を用いることにより得られる複素成分を有する解析信号の位相の変化から算出されてもよい。一般に、ヒルベルト変換は、入力信号について、入力信号を実部とし、π／２の位相遅延器を経由した信号を虚部とする複素成分を有する解析信号に変換することができる。これにより、例えば、連続する２時刻における解析信号の位相の時間変化を抽出することで、入力信号である候補信号の周波数を算出することが可能である。なお、本実施形態において各候補信号の周波数はヒルベルト変換により得られる解析信号により算出されるが、上記周波数の算出方法についてはかかる例に限定されない。例えば、ＤＦＴ（Ｄｉｓｃｒｅｔｅ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍ：離散フーリエ変換）やウェーブレット解析などを用いて各候補信号のピークとなる周波数成分を算出することも可能である。ただし、上述したように、各候補信号は微細な揺らぎを含む一次元の波形を有しているため、各候補信号をフーリエ変換等により周波数領域に変換すると、上記のような微細な揺らぎを表現する周波数成分に関する情報が希薄化することが考えられる。そのため、各候補信号について、ヒルベルト変換など、時間領域において解析する手法を用いて周波数を算出することが好ましい。
【００７２】
（角度変化量による指標値の算出）
また、信号決定部１３３は、変換前のビート信号Ｄ（ｔ）のＩＱ平面における分布に基づいて推定される分布の中心を基準とするビート信号Ｄ（ｔ）の回転角度の変化量に応じて評価値ｓを決定してもよい。ここで、分布の中心、および回転角度は、第２変換部１３２において用いられた第２候補信号Ｒ２（ｔ）の算出方法における分布の中心、および回転角度θ（ｔ）を意味する。例えば、呼吸の振動の変位が比較的小さい場合、回転角度θ（ｔ）の変化量θ’（ｔ）の変化も比較的小さくなる。一方で、呼吸の振動の変位が比較的大きい場合、θ’（ｔ）の変化も比較的大きくなる。そのため、θ’（ｔ）の大きさや偏差に基づいて呼吸の振動の変位を把握することが可能である。つまり、θ’（ｔ）の変化について評価することにより、第１候補信号Ｒ１（ｔ）または第２候補信号Ｒ２（ｔ）の一方又は双方を選択するための評価値ｓを決定することが可能である。
【００７３】
例えば、信号決定部１３３は、以下の数式５を用いて、角度の変化量θ’（ｔ）の偏差に基づいて、評価値ｓを決定するための指標値ｘ２（ｔ）を算出してもよい。
【００７４】
【数５】
0006477199.tif 000006
【００７５】
ここで、Ｌは任意の遅延時間、Ｅは期待値を意味する。例えば、呼吸の振動が比較的小さい場合は、指標値ｘ２（ｔ）の値は小さくなる。一方で、呼吸の振動が比較的大きい場合は、指標値ｘ２（ｔ）の値は大きくなる。
【００７６】
なお、呼吸の振動が比較的小さい場合は、ビート信号Ｄ（ｔ）の回転角度が小さいため、変化量θ’（ｔ）が異常値を示す場合も多い。そのため、例えば信号決定部１３３は、変化量θ’（ｔ）の分布についてメディアンフィルタ、移動平均による平滑化処理、または外れ値の除去処理等を実施することにより、呼吸の振動が比較的小さい場合におけるθ’（ｔ）の値を補正してもよい。これにより、変化量θ’（ｔ）による各候補信号の選択をより適切に実施することができる。また、変化量θ’（ｔ）の値は、呼吸の振動が小さい場合は小さくなり、呼吸の振動が大きい場合は大きくなる。そのため、指標値ｘ２（ｔ）の値は、例えば、θ’（ｔ）と所定の閾値とを比較することにより算出されてもよいし、また、変化量θ’（ｔ）の分布から決定されるマハラノビス距離に基づいて算出されてもよい。
【００７７】
（その他の指標値）
なお、本実施形態においては、各候補信号の周波数、およびビート信号Ｄ（ｔ）の角度の変化量に基づいて指標値ｘ１（ｔ）およびｘ２（ｔ）がそれぞれ算出されたが、本発明はかかる例に限定されず、多様な方法により指標値が算出されてもよい。例えば、信号決定部１３３は、各候補信号の振幅に含まれる特徴に基づいて指標値を算出してもよい。より具体的に説明すると、信号決定部１３３は、第１候補信号Ｒ１（ｔ）と第２候補信号Ｒ２（ｔ）の同時刻における振幅を比較し、複数の時刻における比較結果の分布と上記時刻における比較結果とに基づいて評価値ｓを決定するための指標値が算出されてもよい。これにより、一方の候補信号の振幅に異常があった際に、すぐに他方の候補信号を用いるように評価値ｓの値を決定することが可能である。また、信号決定部１３３は、振幅に含まれる特徴に基づいて指標値を算出する際に、双方の候補信号の振幅を正規化してもよい。これにより、双方の候補信号の振幅の大きさ以外の特徴を抽出が容易となり、評価値ｓを決定するための要因を増やすことが可能となる。
【００７８】
［４−２．評価値の決定］
信号決定部１３３は、算出された指標値ｘ１（ｔ）またはｘ２（ｔ）の少なくともいずれかの分布に基づいて、評価値ｓを決定する処理を行う（Ｓ２０２）。例えば、信号決定部１３３は、算出された指標値ｘ１（ｔ）またはｘ２（ｔ）の少なくともいずれかと、所定の閾値とを比較することにより、評価値ｓを決定してもよい。また、信号決定部１３３は、指標値ｘ１（ｔ）またはｘ２（ｔ）の少なくともいずれかについて、機械学習の手法により最適化されたモデルを用いて評価値ｓを決定してもよい。例えば、信号決定部１３３は、マハラノビス距離や線形判別関数などの判別分析や、ＲＶＭ（Ｒｅｌｅｖａｎｃｅ  Ｖｅｃｔｏｒ  Ｍａｃｈｉｎｅ：関連ベクトルマシン）、スパースロジスティック回帰、およびベイズロジスティック回帰等のロジスティック回帰、並びに、Ｄｅｅｐ  Ｌｅａｒｎｉｎｇ（深層学習）等を含むニューラルネットワークなどを用いて評価値ｓを決定してもよい。以下、評価値ｓの決定手法の例について説明する。
【００７９】
（指標値の分布による評価）
例えば、信号決定部１３３は、算出された指標値ｘ１（ｔ）またはｘ２（ｔ）の少なくともいずれかが所定の閾値以上であるか否かに応じて、評価値ｓを０または１のいずれかに決定してもよい。また、信号決定部１３３は、指標値ｘ１（ｔ）またはｘ２（ｔ）の少なくともいずれかの値と上記の所定の閾値との差に基づいて、評価値ｓを０〜１の間において決定してもよい。
【００８０】
図９は、本発明の一実施形態に係る信号決定部１３３により算出された指標値の分布を示す図である。横軸は指標値ｘ１（ｔ）を、縦軸は指標値ｘ２（ｔ）を示す。なお、指標値ｘ１（ｔ）は数式４により算出された値であり（ｋ＝１）、指標値ｘ２（ｔ）は数式５により算出された値である。また、図中の△は、第１候補信号Ｒ１（ｔ）に高周波成分が含まれている場合を示し、図中の○は、第１候補信号Ｒ１（ｔ）に高周波成分が含まれていない場合を示している。図９を参照すると、指標値ｘ１（ｔ）およびｘ２（ｔ）がともに比較的高い値を示す場合（境界Ｂ１の上部にプロットされた場合）、高周波成分を含む波形が多くなっている。この場合、例えば、境界Ｂ１を判別境界として、判別境界Ｂ１よりも上部にプロットされている指標値を示す場合、高周波成分を含む可能性の低い第２候補信号Ｒ２（ｔ）が選択されるように評価値ｓが決定されてもよい。具体的には、信号決定部１３３は、判別境界Ｂ１よりも上部にプロットされた場合は評価値ｓを１と決定し、判別境界Ｂ１よりも下部にプロットされた場合は評価値ｓを０と決定してもよい。これにより、評価値ｓを容易に決定することが可能となる。また、信号決定部１３３は、指標値ｘ１（ｔ）およびｘ２（ｔ）についてプロットされた点と判別境界Ｂ１との距離から評価値ｓを決定してもよい。
【００８１】
また、信号決定部１３３は、プロットされた点と判別境界Ｂ１との距離に基づき、任意の関数を用いて評価値ｓを算出してもよい。この場合、信号決定部１３３は、例えば、判別境界Ｂ１の評価値ｓを０．５として、破線Ｂ２上の評価値ｓを０．７５、破線Ｂ３上の評価値ｓを０．２５となるような関数を用いて評価値ｓを算出してもよい。これにより、評価値ｓをより適切に決定することが可能である。また、図９に示した例では、単純な線形の判別境界Ｂ１に基づいて評価値ｓが算出され得るが、例えば、ベイズの識別規則やＳＶＭ（Ｓｕｐｐｏｒｔ  Ｖｅｃｔｏｒ  Ｍａｃｈｉｎｅ：サポートベクトルマシン）、カーネルトリック等を用いて、判別境界が設定されてもよい。これにより、より適切な評価値ｓを決定することが可能となる。
【００８２】
（ロジスティック回帰による評価）
また、信号決定部１３３は、ロジスティック回帰により構築した確率的識別モデルから、評価値ｓを決定してもよい。ここで、評価値ｓは数式６のシグモイド関数により定義され、評価値ｓを算出するための確率的識別モデルは、数式７に示された尤度ｐ（ｙ｜ｗ）で与えられるロジスティック回帰により構築される。なお、数式６および７において、σ（ｘ）はシグモイド関数を意味する。
【００８３】
【数６】
0006477199.tif 000007
【００８４】
【数７】
0006477199.tif 000008
【００８５】
ただし、ｙは教師データ、ｗは学習パラメータ、ｘｎは指標値が格納されるベクトルであり、これらは以下の数式８〜１０により定義される。
【００８６】
【数８】
0006477199.tif 000009
【００８７】
【数９】
0006477199.tif 000010
【００８８】
【数１０】
0006477199.tif 000011
【００８９】
評価値ｓは、数式７において尤度ｐ（ｙ｜ｗ）が最大となる場合の学習パラメータｗ、および指標値が格納されるベクトルｘｎに基づいて、数式６を用いて算出される。ロジスティック回帰を用いることにより、第１候補信号Ｒ１（ｔ）と第２候補信号Ｒ２（ｔ）のどちらが適切な波形であるかを決定する評価値ｓを、指標値の分布に基づいて算出することが可能である。なお、数式１０においてはベクトルｘｎに格納される指標値はｘ１（ｔ）およびｘ２（ｔ）の双方であったが、指標値ｘ１（ｔ）またはｘ２（ｔ）のいずれかのみがベクトルｘｎに格納されてもよい。また、周波数や角度変化量以外の多次元の指標値をベクトルｘｎに追加してもよい。これにより、評価値ｓを様々なパラメータから決定することが可能である。
【００９０】
（ニューラルネットワークによる評価）
信号決定部１３３は、ニューラルネットワークを用いて、評価値ｓを算出してもよい。例えば、信号決定部１３３は、上記のような指標値を入力するための入力層、第１候補信号Ｒ１（ｔ）および第２候補信号Ｒ２（ｔ）を選択するための出力値（選定出力値）を出力するための出力層、および中間層からなる３層のフィードフォワード型ニューラルネットワークを構成してもよい。中間層は、例えば複数のノードにより構成されてもよく、また、Ｄｅｅｐ  Ｌｅａｒｎｉｎｇのように複数の階層により構成されてもよい。ここで、ニューラルネットワークにより評価値ｓを決定するために、信号決定部１３３は、例えば、誤差逆伝播法等によりニューラルネットワークを学習させる必要がある。例えば、信号決定部１３３は、教師信号付きの指標値を入力層に入力し、出力された選定出力値と教師信号とを比較し、比較の結果生じた誤差を中間層にフィードバックすることにより、最適な選定出力値が出力されるようにニューラルネットワークを学習させる。これにより、中間層の各ノードにおいて演算された選定出力値の精度を向上させることが可能である。
【００９１】
そして、信号決定部１３３は、学習されたニューラルネットワークを用いて評価値ｓを決定する。信号決定部１３３は、例えば、入力層に各指標値を与え、その後出力層から出力された第１候補信号Ｒ１（ｔ）および第２候補信号Ｒ２（ｔ）の選定出力値に基づいて評価値ｓを決定してもよい。例えば、評価値ｓは、第２候補信号Ｒ２（ｔ）の選定出力値をシグモイド関数により変換された値としてもよい。また、評価値ｓは、各候補信号の選定出力値の合計に対する第２候補信号Ｒ２（ｔ）の選定出力値の比率としてもよい。これにより、第１候補信号Ｒ１（ｔ）および第２候補信号Ｒ２（ｔ）のどちらが適切な波形であるかを決定する評価値ｓを、ニューラルネットワークから直接算出することが可能である。また、ニューラルネットワークを用いることにより、周波数や角度変化量以外の多次元の指標値も入力層に入力することが可能である。これにより、評価値ｓを様々なパラメータから決定することが可能である。
【００９２】
（バイアスの設定）
また、第１候補信号Ｒ１（ｔ）または第２候補信号Ｒ２（ｔ）のいずれかを優先的に選択できるように、評価値ｓの決定プロセスにおいて、予めバイアスが加えられてもよい。例えば、呼吸の振動の変位が放射波の波長の１／３〜１／２程度である場合、第１候補信号Ｒ１（ｔ）および第２候補信号Ｒ２（ｔ）は、いずれも呼吸の振動の波形を同程度の精度で表現し得る。この場合、呼吸の振動の方向（吸う方向および吐く方向）を反映する第２候補信号Ｒ２（ｔ）を優先的に用いることが望ましい。そのため、例えば、信号決定部１３３は、評価値ｓの決定の際に、第２候補信号Ｒ２（ｔ）に重みを与えるようにバイアスを加えておいてもよい。この場合、具体的には、上記のようなロジスティック回帰における確率的識別モデルやニューラルネットワークの中間層において、第２候補信号Ｒ２（ｔ）に関する値が第１候補信号Ｒ１（ｔ）に関する値と比較して高い値が出力されるように調整されてもよい。これにより、第１候補信号Ｒ１（ｔ）および第２候補信号Ｒ２（ｔ）がいずれも呼吸の振動の波形を同程度の精度で表現している場合において、呼吸の振動の方向を反映した第２候補信号Ｒ２（ｔ）を優先的に選択することが可能となる。
【００９３】
［４−３．一次元信号の決定］
次に、信号決定部１３３は、評価値ｓに基づいて、出力すべき一次元信号ｒ（ｔ）を決定する処理を行う（Ｓ２０３）。まず、信号決定部１３３は、評価値ｓに基づいて、第１候補信号Ｒ１（ｔ）または第２候補信号Ｒ２（ｔ）の一方、または双方を選択する。そして、信号決定部１３３は必要な場合に応じて、複数の候補信号の位相を整合させる。そして、信号決定部１３３は、候補信号の選択時に双方の候補信号を選択した場合、双方の候補信号を統合する処理を行う。以下、一次元信号ｒ（ｔ）の決定処理について説明する。
【００９４】
（評価値の利用）
まず、信号決定部１３３は、評価値ｓに基づいて第１候補信号Ｒ１（ｔ）または第２候補信号Ｒ２（ｔ）の一方、または双方を選択する。例えば、信号決定部１３３は、評価値ｓが所定の閾値未満である場合に第１候補信号Ｒ１（ｔ）を、評価値ｓが所定の閾値以上である場合に第２候補信号Ｒ２（ｔ）を一次元信号ｒ（ｔ）として決定してもよい。また、信号決定部１３３は、双方の候補信号を統合することにより一次元信号ｒ（ｔ）を決定する場合、後述するように、評価値ｓを各候補信号の重みとして用いてもよい。
【００９５】
（位相の整合）
信号決定部１３３は、例えば、一次元信号ｒ（ｔ）として一時刻前に一方の候補信号を選択したあとに他方の候補信号に切り替える場合、または、双方の候補信号を統合する場合、双方の候補信号の位相を整合させる必要がある。例えば、第１候補信号Ｒ１（ｔ）の位相φ１と第２候補信号Ｒ２（ｔ）の位相φ２は必ずしも一致しているとは限らない。そのため、信号決定部１３３は、一方の候補信号の位相を他方の候補信号の位相に整合させる。これにより、一次元信号ｒ（ｔ）の信号の切り替え処理や、双方の候補信号の統合処理において、各候補信号の位相の不整合による波形の乱れを防ぐことができる。なお、本実施形態では、第１候補信号Ｒ１（ｔ）の位相を第２候補信号Ｒ２（ｔ）の位相に整合させる処理が実施される。以下、位相の整合処理について説明する。
【００９６】
信号決定部１３３は、位相の整合処理の前に、各候補信号の振幅を揃えるための正規化を行う。例えば、各候補信号における出力値の標準偏差、最大値と最小値の差、または絶対値の最大値等から、各候補信号の振幅が算出されてもよい。そして、算出された振幅に基づいて線形変換等により各候補信号の振幅が正規化されてもよい。また、一方の候補信号の振幅を他方の候補信号の振幅に合わせてもよい。この場合、信号決定部１３３は、呼吸の１周期よりも長い時間の信号を候補信号の各々から切り出し、切り出した信号から各候補信号の振幅を算出してもよい。例えば、呼吸の１周期は約３秒であるので、約３秒〜１０秒程度の信号を切り出すことが好ましい。これにより、呼吸の変動による振幅の変動の影響を含んだうえで各候補信号の振幅を算出することが可能となる。なお、本実施形態では双方の候補信号の振幅について正規化を行うが、後処理において振幅の影響を問題としない場合は、振幅の正規化が行われなくてもよい。また、本実施形態では信号決定部１３３が正規化を行うが、ビート信号取得部１１０やフィルタ部１２０、第１変換部１３１や第２変換部１３２がビート信号Ｄ（ｔ）や第１候補信号Ｒ１（ｔ）、第２候補信号Ｒ２（ｔ）について正規化を行ってもよい。また、信号決定部１３３は、正規化された各候補信号の振幅を元に戻す処理を行ってもよい。
【００９７】
各候補信号の正規化処理が実施された場合、各候補信号はそれぞれ以下の数式１１および数式１２のように表現される。
【００９８】
【数１１】
0006477199.tif 000012
【００９９】
【数１２】
0006477199.tif 000013
【０１００】
ここで、本実施形態においては、上述のとおり、信号決定部１３３は、第１候補信号Ｒ１（ｔ）の位相を第２候補信号Ｒ２（ｔ）の位相に整合させる処理を行う。信号決定部１３３は、ヒルベルト変換を用いて第１候補信号Ｒ１（ｔ）を解析信号Ｒｅｆ（ｔ）に変換する（数式１３）。
【０１０１】
【数１３】
0006477199.tif 000014
【０１０２】
次に、信号決定部１３３は、Ｒ２（ｔ）とＲｅｆ（ｔ）を乗算し、位相差Δφを抽出する（数式１４）。
【０１０３】
【数１４】
0006477199.tif 000015
【０１０４】
ここで、Δｆ＝ｆｒ１−ｆｒ２、Δφ＝φ１−φ２であり、また、ＦＬＰＦ（ｘ）は、ｆｒ１＋ｆｒ２により表現される周波数成分を低減するローパスフィルタを表す式である。このローパスフィルタは、遮断特性が急峻であり、例えば０．１Ｈｚ以上の周波数成分を除去するフィルタであってもよい。これにより、Ｒ２（ｔ）とＲｅｆ（ｔ）を乗算させたことにより生じるｆｒ１＋ｆｒ２の成分を除去することが可能である。また、Δｆは各候補信号の周波数であり、ヒルベルト変換等によりΔｆを算出することは可能である。しかし、各候補信号の周波数は近い値を有するので、信号決定部１３３は、Δｆを０とみなしてもよい。これにより、Ｐｈ（ｔ）を、各候補信号の位相差を示す信号とみなすことができる（数式１５）。
【０１０５】
【数１５】
0006477199.tif 000016
【０１０６】
信号決定部１３３は、先に得たＲｅｆ（ｔ）と上記のＰｈ（ｔ）の複素共役を乗算することにより、第１候補信号Ｒ１（ｔ）の周波数成分と第２候補信号Ｒ２（ｔ）の位相を有する複素信号Ｒｑ（ｔ）を得ることができる（数式１６）。そして、信号決定部１３３は、複素信号Ｒｑ（ｔ）の実部を、第２候補信号Ｒ２（ｔ）の位相を有する第１候補信号Ｒ１ａ（ｔ）として得ることができる（数式１７）。
【０１０７】
【数１６】
0006477199.tif 000017
【０１０８】
【数１７】
0006477199.tif 000018
【０１０９】
以上の位相の整合処理により、信号決定部１３３は、候補信号の切り替え時や候補信号の統合時において、各候補信号の位相を揃えた状態で処理を行うことが可能となる。これにより連続的な一次元信号ｒ（ｔ）を出力することが可能となる。なお、例えば切り替え時において位相を整合させた後に一次元信号ｒ（ｔ）を出力した場合、揺らぎを多く伴う波形においては、波形が歪な形に変形する可能性がある。そのため、信号決定部１３３は、例えば位相を整合させたのちに、上記位相に緩やかに変化するよう、指数平滑化平均等を用いて処理してもよい。また、一方の候補信号の位相が整合のために変化した場合、時間方向に対して上記候補信号の波形が変化することが考えられる。そのため、信号決定部１３３は、遅延処理等を行うことにより、緩やかに波形を変化させながら、上記候補信号の位相を他方の候補信号の位相に合わせこんでもよい。
【０１１０】
（統合処理）
信号決定部１３３は、双方の候補信号を用いて一次元信号ｒ（ｔ）を決定する場合、位相を整合させたうえで、双方の候補信号を統合する処理を実施してもよい。例えば、信号決定部１３３は、以下の数式１８のように評価値ｓを重みとして用いて一次元信号ｒ（ｔ）を決定することができる。
【０１１１】
【数１８】
0006477199.tif 000019
【０１１２】
数式１８で示された一次元信号ｒ（ｔ）は、位相補正された第１候補信号Ｒ１ａ（ｔ）と第２候補信号Ｒ２（ｔ）の加重和として表現される。これにより、信号決定部１３３は、評価値ｓの値に追随して一次元信号ｒ（ｔ）の値を決定することが可能である。
【０１１３】
以上のような一次元信号ｒ（ｔ）の決定処理によって、呼吸の振動の変位に応じて変動する第１候補信号Ｒ１（ｔ）および第２候補信号Ｒ２（ｔ）のうち、呼吸の振動の波形をより適切に表現している信号を選択することができる。これにより、呼吸の振動の変位の大きさに関わらず、呼吸の振動をより正確に反映させた一次元信号ｒ（ｔ）を出力することが可能である。そのため、後処理である周期推定処理において、微細な揺らぎを含み、振動の変位が大きく変動するような呼吸の振動の周期をより正確に推定することが可能となる。
【０１１４】
＜５．振動状態推定装置の動作例＞
次に、本発明の一実施形態に係る振動状態推定装置１０の動作例について説明する。図１０は、本発明の一実施形態に係る振動状態推定装置１０の動作例を示すフローチャートである。
【０１１５】
まず、ビート信号取得部１１０は、ドップラーセンサ２から出力されたビート信号Ｄ（ｔ）を取得する（Ｓ３０１）。そして、フィルタ部１２０は、取得されたビート信号Ｄ（ｔ）から直流成分などの低周波成分を低減または除去するフィルタ処理を行う（Ｓ３０２）。
【０１１６】
次に、信号変換部１３０は、フィルタ処理されたビート信号Ｄ（ｔ）を一次元信号ｒ（ｔ）に変換する（Ｓ３０３）。そして、周波数推定部１４０は、変換された一次元信号ｒ（ｔ）の有する周波数を、一または複数の手段を用いて推定する（Ｓ３０４）。次に、基準位置推定部１５０は、推定された周波数を用いて、一次元信号ｒ（ｔ）の基準位置を推定する（Ｓ３０５）。そして、周期算出部１６０は、一次元信号ｒ（ｔ）における連続する基準位置の間隔から、呼吸の周期を算出する（Ｓ３０６）。周期の算出後、振動状態推定装置１０は、引き続き振動の周期算出を実施するか否かを、ビート信号Ｄ（ｔ）の取得の有無により、または、振動状態推定装置１０のユーザの入力等により判定する（Ｓ３０７）。引き続き振動の周期算出が実施される場合（ＹＥＳ）、再度ステップＳ３０１のビート信号取得処理が実施される。一方、振動の周期算出が実施されない場合（ＮＯ）、振動状態推定装置１０は動作を終了する。
【０１１７】
＜６．ハードウェア構成例＞
以上、本発明の一実施形態に係る振動状態推定装置１０の動作例について説明した。上述した振動状態推定装置１０の情報処理は、ソフトウェアと、振動状態推定装置１０との協働により実現される。以下では、本発明の実施形態に係る振動状態推定装置１０のハードウェア構成を説明する。
【０１１８】
図１１は、本発明の実施形態に係る振動状態推定装置１０のハードウェア構成例を示すブロック図である。図１１を参照すると、振動状態推定装置１０は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）９０１と、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）９０２と、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）９０３と、ホストバス９０４と、を備える。また、振動状態推定装置１０は、ブリッジ９０５と、外部バス９０６と、インタフェース９０７と、入力装置９０８と、出力装置９０９と、ストレージ装置９１０と、ドライブ９１１と、ネットワークインタフェース９１２と、を備える。
【０１１９】
ＣＰＵ９０１は、演算処理装置および制御装置として機能し、各種プログラムに従って振動状態推定装置１０内の動作全般を制御する。また、ＣＰＵ９０１は、マイクロプロセッサであってもよい。なお、ＣＰＵ９０１は、振動状態推定装置１０内の動作全般またはその一部を制御する。例えば、ＲＯＭ９０２は、ＣＰＵ９０１が使用するプログラムや演算パラメータ等を記憶する。ＲＡＭ９０３は、ＣＰＵ９０１の実行において使用するプログラムや、その実行において適宜変化するパラメータ等を一時記憶する。これらはＣＰＵバスなどから構成されるホストバス９０４により相互に接続されている。
【０１２０】
ホストバス９０４は、ブリッジ９０５を介して、ＰＣＩ（Ｐｅｒｉｐｈｅｒａｌ  Ｃｏｍｐｏｎｅｎｔ  Ｉｎｔｅｒｃｏｎｎｅｃｔ／Ｉｎｔｅｒｆａｃｅ）バスなどの外部バス９０６に接続されている。なお、必ずしもホストバス９０４、ブリッジ９０５および外部バス９０６を分離構成する必要はなく、１つのバスにこれらの機能を実装してもよい。
【０１２１】
入力装置９０８は、マウス、キーボード、タッチパネル、ボタン、マイクロフォン、スイッチおよびレバーなどユーザが情報を入力するための入力手段と、ユーザによる入力に基づいて入力信号を生成し、ＣＰＵ９０１に出力する入力制御回路などから構成されている。振動状態推定装置１０のユーザは、当該入力装置９０８を操作することにより、振動状態推定装置１０に対して各種のデータを入力したり処理動作を指示したりすることができる。
【０１２２】
出力装置９０９は、例えば、ＣＲＴディスプレイ装置、液晶ディスプレイ（ＬＣＤ）装置、ＯＬＥＤ装置およびランプなどの表示装置を含む。さらに、出力装置９０９は、スピーカ及びヘッドホンなどの音声出力装置を含む。出力装置９０９は、例えば、再生されたコンテンツを出力する。具体的には、表示装置は再生された映像データ等の各種情報をテキストまたはイメージで表示する。一方、音声出力装置は、再生された音声データや表示装置に表示されたテキストデータ等を音声に変換して出力する。
【０１２３】
ストレージ装置９１０は、本発明の実施形態に係る振動状態推定装置１０におけるデータ格納用の装置である。ストレージ装置９１０は、記憶媒体、記憶媒体にデータを記録する記録装置、記憶媒体からデータを読み出す読み出し装置および記憶媒体に記憶されたデータを削除する削除装置などを含んでも良い。ストレージ装置は、例えば、ＨＤＤ（Ｈａｒｄ  Ｄｉｓｃ  Ｄｒｉｖｅ）やＳＳＤ（Ｓｏｌｉｄ  Ｓｔａｔｅ  Ｄｒｉｖｅ）で構成される。このストレージ装置９１０は、ＣＰＵ９０１が実行するプログラムや各種データを格納する。なお、ストレージ装置９１０は、不図示の記憶部の機能を実現する。
【０１２４】
ドライブ９１１は、記憶媒体用リーダライタであり、振動状態推定装置１０に内蔵、あるいは外付けされる。ドライブ９１１は、装着されている磁気ディスク、光ディスク、光磁気ディスク、または半導体メモリ等のリムーバブル記憶媒体９６に記録されている情報を読みだして、ＲＡＭ９０３に出力する。また、ドライブ９１１は、リムーバブル記憶媒体９６に情報を書き込むこともできる。
【０１２５】
ネットワークインタフェース９１２は、例えば、他の装置に接続するための通信デバイス等で構成された通信インタフェースである。また、ネットワークインタフェース９１２は、無線ＬＡＮ（Ｌｏｃａｌ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）対応通信装置、ＬＴＥ（Ｌｏｎｇ  Ｔｅｒｍ  Ｅｖｏｌｕｔｉｏｎ）対応通信装置、またはブルートゥース通信装置であってもよい。また、ネットワークインタフェース９１２は、有線による通信を行うワイヤー通信装置であってもよい。なお、ネットワークインタフェース９１２は、不図示の通信部の機能を実現する。
【０１２６】
以上、振動状態推定装置１０のハードウェア構成の一例を示した。上記の各構成要素は、汎用的な部材を用いて構成されていてもよいし、各構成要素の機能に特化したハードウェアにより構成されていてもよい。かかる構成は、実施する時々の技術レベルに応じて適宜変更され得る。
【０１２７】
＜７．まとめ＞
ここまで、図１〜図１１を用いて、本発明の一実施形態について説明した。本発明の一実施形態によれば、振動状態推定装置１０は、取得したビート信号を複数の信号変換手段を用いて一次元の候補信号に変換しつつ、変換された複数の候補信号のいずれが呼吸の振動をより正確に表現しているかを評価し、評価結果に基づいて出力すべき一次元信号を決定する。これにより、呼吸の振動の変動が大きい場合であっても、呼吸の振動をより正確に反映させた一次元信号を抽出することが可能である。したがって、後処理である振動の周期推定処理において、実際の呼吸の振動に近い波形に基づいて振動の周期を推定することが可能となるので、振動の周期をより高い精度で推定することができる。
【０１２８】
以上、添付図面を参照しながら本発明の好適な実施形態について詳細に説明したが、本発明はかかる例に限定されない。本発明の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本発明の技術的範囲に属するものと了解される。
【０１２９】
例えば、上記実施形態では、振動状態推定装置１０は、ドップラーセンサ２から出力されるビート信号をリアルタイムで取得し、並行して解析する動作を実施していたが、本発明はかかる例に限定されない。例えば、他の実施形態に係る振動状態推定装置は、一度取得した信号を記憶部に保存し、推定処理を実施したい際に記憶部から再度信号を取り出して解析を行ってもよい。これにより、過去に取得したビート信号についてまとめて解析を行うことが可能となる。
【０１３０】
また、上記実施形態では、被験者がドップラーセンサ２による検知が可能なエリアに存在していることが想定されていたが、本発明はかかる例に限定されない。例えば、被験者が椅子やベッドなどの位置にいる場合にのみドップラーセンサが被験者を検知することができるように、ドップラーセンサが設置されてもよい。これにより、解析対象となる被験者を限定することが可能となる。また、振動状態推定装置は、被験者の有無や被験者による歩行などの大きな活動を検知する手段を組み合わせてもよい。これにより、被験者が睡眠しているなど限られた状態にある場合にのみ推定処理を実施することが可能となる。
【０１３１】
また、上記実施形態では、生体による呼吸の振動を解析対象とする例を説明したが、本発明はかかる例に限定されない。例えば、本発明は、生体の心拍、脈拍、およびけいれん運動、並びに、機械の振動およびクランク等の往復運動など、多様な振動運動に対して適用可能である。この場合、ドップラーセンサの放射波の波長を解析対象の振動の変位に応じて適切に設定することにより、あらゆる振動の微細な揺らぎを含む一次元信号を抽出することが可能である。
【０１３２】
また、上記実施形態では、ドップラーセンサ２と振動状態推定装置１０とが振動状態推定システム１を構成していたが、本発明は係る例に限定されない。例えば、ドップラーセンサ２と振動状態推定装置１０とが一体となった振動状態推定装置または振動状態推定システムが提供されてもよい。
【０１３３】
また、本明細書の振動状態推定装置１０および信号変換部１３０の処理における各ステップは、必ずしもフローチャートとして記載された順序に沿って時系列に処理する必要はない。例えば、振動状態推定装置１０または信号変換部１３０の処理における各ステップは、フローチャートとして記載した順序と異なる順序で処理されても、並列的に処理されてもよい。
【０１３４】
また、振動状態推定装置１０に内蔵されるＣＰＵ、ＲＯＭおよびＲＡＭなどのハードウェアを、上記の振動状態推定装置１０の各構成と同等の機能を発揮させるためのコンピュータプログラムも作成可能である。また、該コンピュータプログラムを記憶させた記憶媒体も提供される。
【符号の説明】
【０１３５】
１      振動状態推定システム
２      ドップラーセンサ
１０    振動状態推定装置
１１０  ビート信号取得部
１２０  フィルタ部
１３０  信号変換部
１３１  第１変換部
１３２  第２変換部
１３３  信号決定部
１４０  周波数推定部
１５０  基準位置推定部
１６０  周期算出部
【図１】
0006477199.tif 000020
【図２】
0006477199.tif 000021
【図３】
0006477199.tif 000022
【図４】
0006477199.tif 000023
【図５】
0006477199.tif 000024
【図６】
0006477199.tif 000025
【図７】
0006477199.tif 000026
【図８】
0006477199.tif 000027
【図９】
0006477199.tif 000028
【図１０】
0006477199.tif 000029
【図１１】
0006477199.tif 000030