JP7433554B2 - バイタル測定装置、バイタル測定方法及びバイタル測定システム - Google Patents

Description

本開示は、バイタル測定装置、バイタル測定方法及びバイタル測定システムに関するものである。

被測定者のバイタルを測定する測定装置がある（特許文献１を参照）。当該測定装置は、室内に存在している被測定者に向けてマイクロ波を送信したのち、被測定者による反射後のマイクロ波である反射波を受信し、反射波の受信信号を出力するアンテナを備えている。また、当該測定装置は、アンテナから出力された受信信号に含まれている同相信号を抽出する第１の混合器と、アンテナから出力された受信信号に含まれている直交信号を抽出する第２の混合器と、当該同相信号と当該直交信号とを含む複素信号を生成する信号処理装置とを備えている。
アンテナから送信されたマイクロ波が、被測定者だけでなく、室内の壁に反射されたような場合、当該信号処理装置により生成される複素信号は、被測定者による反射波に係る複素信号と、壁による反射波に係る複素信号とが重ね合わされた信号となる。壁は、静止物であるため、壁による反射波に係る複素信号の位相は、一定である。一方、被測定者の胸部は、呼吸に伴って往復運動するため、被測定者による反射波に係る複素信号の位相は、時間の経過に伴って変化する。当該信号処理装置は、生成した複素信号の中から、位相が一定の複素信号を除去することで、被測定者による反射波に係る複素信号を取得し、被測定者による反射波に係る複素信号から、被測定者のバイタルを検出する。

特開２０２０－１５７０００号公報

特許文献１に開示されている測定装置では、複数の被測定者が室内に存在していれば、それぞれの被測定者のバイタルに係る複素信号が互いに重ね合わされる。信号処理装置は、生成した複素信号の中から、位相が一定の複素信号を除去したとしても、それぞれの被測定者による反射波に係る複素信号を分離することができない。このため、当該信号処理装置は、それぞれの被測定者のバイタルを検出することができないという課題があった。

本開示は、上記のような課題を解決するためになされたもので、被測定者が複数の場合であっても、それぞれの被測定者のバイタルを検出することができるバイタル測定装置及びバイタル測定方法を得ることを目的とする。

本開示に係るバイタル測定装置は、対象物体による反射波を受信するアンテナから、反射波の受信信号を取得する信号取得部と、信号取得部により取得された受信信号を時間方向にフーリエ変換するフーリエ変換部とを備えている。また、バイタル測定装置は、フーリエ変換部によるフーリエ変換後の信号を用いて、バイタル測定装置からの複数の距離ビンのそれぞれに対応する複素電力の２次元方位マップを算出するマップ算出部と、マップ算出部により算出されたそれぞれの距離ビンに対応する複素電力の２次元方位マップの時間変化に基づいて、対象物体に含まれる被測定者が存在している位置を特定し、被測定者が存在している位置についての複素電力の位相の時間変化から、被測定者のバイタルを推定するバイタル推定部とを備えている。

本開示によれば、被測定者が複数の場合であっても、それぞれの被測定者のバイタルを検出することができる。

実施の形態１に係るバイタル測定装置２０を含むバイタル測定システムを示す構成図である。 実施の形態１に係るバイタル測定装置２０のハードウェアを示すハードウェア構成図である。 バイタル測定装置２０が、ソフトウェア又はファームウェア等によって実現される場合のコンピュータのハードウェア構成図である。 バイタル測定装置２０の処理手順であるバイタル測定方法を示すフローチャートである。 それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、信号発生器１２ａにより生成されるアップチャープの信号Ｔｘ（１）～Ｔｘ（Ｎ）を示す説明図である。 或る送信サイクルｃにおけるＮ個のフーリエ変換後の信号Ｓ（ｒ，ｇ，ｈ，ｃ）を示す説明図である。 それぞれの距離ビンｒ（ｒ＝１，・・・，Ｒ）に対応する複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の２次元方位マップを示す説明図である。 同一の距離ビンｒに属し、かつ、同一の２次元方位（ａｚ，ｅｌ）に属している複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）における位相の時間変化を示す説明図である。 存在位置特定部２５により算出される位相変化信号θ（ｒ，ｃ）を示す説明図である。 呼吸スペクトルＳ_ＲＲ（ｒ，ｓｆ）の一例を示す説明図である。 スカログラムの一例を示す説明図である。 スロータイム方向の２次元スペクトルＷ（ｒ，ｆ，ｓｆ）の一例を示す説明図である。 心拍数推定部２８による最大比合成処理を示す説明図である。 心拍数推定部２８により得られる心拍スペクトルＳ_ＨＲ（ｒ，ｓｆ）の一例を示す説明図である。 実施の形態２に係るバイタル測定装置２０を含むバイタル測定システムを示す構成図である。 実施の形態２に係るバイタル測定装置２０のハードウェアを示すハードウェア構成図である。 信号発生器１２ａにより生成されるヒット数Ｑのアップチャープの信号Ｔｘ（１）～Ｔｘ（Ｎ）を示す説明図である。 図１８Ａは、アンテナ１１－１～１１－ＮからＴＤＭ（Ｔｉｍｅ Ｄｉｖｉｓｉｏｎ Ｍｕｌｔｉｐｌｅｘｉｎｇ）方式で送信波が放射されたのち、信号取得部２１により取得される受信データＳ（ｔ，ｇ，ｈ，ｃ）を示す説明図、図１８Ｂは、ＡＤオフセット補正後の受信データＳ’（ｔ，ｇ，ｈ，ｃ）を示す説明図、図１８Ｃは、移動体による反射波に係る信号抑圧後の受信データＳ_Ｙ（ｔ，ｇ，ｈ，ｃ）を示す説明図である。 実施の形態３に係るバイタル測定装置２０を含むバイタル測定システムを示す構成図である。 実施の形態３に係るバイタル測定装置２０のハードウェアを示すハードウェア構成図である。

以下、本開示をより詳細に説明するために、本開示を実施するための形態について、添付の図面に従って説明する。

実施の形態１．
図１は、実施の形態１に係るバイタル測定装置２０を含むバイタル測定システムを示す構成図である。
図２は、実施の形態１に係るバイタル測定装置２０のハードウェアを示すハードウェア構成図である。
図１に示すバイタル測定システムは、センサ１０及びバイタル測定装置２０を備えている。
センサ１０は、Ｎ個のアンテナ１１－１～１１－Ｎ、信号送信部１２、Ｎ個のサーキュレータ１３－１～１３－Ｎ及びＮ個の信号受信部１４－１～１４－Ｎを備えている。Ｎは、２以上の整数である。
アンテナ１１－１～１１－Ｎのそれぞれは、送受信アンテナである。
図１に示すバイタル測定システムでは、ＭＩＭＯ（Ｍｕｌｔｉｐｌｅ－Ｉｎｐｕｔ Ｍｕｌｔｉｐｌｅ－Ｏｕｔｐｕｔ）によって受信信号の分解能を高めるために、センサ１０が、Ｎ個のアンテナ１１－１～１１－Ｎを備えている。そして、アンテナ１１－１～１１－Ｎのそれぞれが、送信アンテナと受信アンテナとを兼ねている。しかし、これは一例に過ぎず、センサ１０が、送信アンテナと受信アンテナとを別々に備えていてもよい。

センサ１０から送信波が放射されるときは、アンテナ１１－１～１１－Ｎの中から、送信波を放射する１つのアンテナが選択される。
アンテナ１１－１～１１－Ｎにおける送信波の送信順序は、決まっている。例えば、アンテナ１１－１、アンテナ１１－２、・・・、アンテナ１１－Ｎの順番に送信順序が決まっている。ただし、これは一例に過ぎず、例えば、アンテナ１１－Ｎ、・・・、アンテナ１１－２、アンテナ１１－１の順番に送信順序が決まっていてもよい。
アンテナ１１－ｎ（ｎ＝１，・・・，Ｎ）は、サーキュレータ１３－ｎから出力された送信信号に係る送信波を対象物体が存在している空間に放射する。アンテナ１１－ｎから放射された送信波は、対象物体によって反射される。対象物体の中には、空間に存在している１人以上の被測定者ｋのほか、空間を形成している部屋の壁、又は、空間に存在している机等の静止物が含まれる。ｋ＝１，・・・，Ｋであり、Ｋは、１以上の整数である。
アンテナ１１－１～１１－Ｎのそれぞれは、対象物体による反射波を受信し、反射波の受信信号をサーキュレータ１３－ｎに出力する。

信号送信部１２は、信号発生器１２ａ及び出力先選択部１２ｂを備えている。
信号送信部１２は、Ｎ個のアンテナ１１－１～１１－Ｎの中から、送信波を放射させる１つのアンテナ１１－ｎを順番に選択する。
信号送信部１２は、選択したアンテナ１１－ｎから送信波を空間に放射させるために、選択したアンテナ１１－ｎと接続されているサーキュレータ１３－ｎに送信信号を出力する。
信号発生器１２ａは、例えば、時間の経過に伴って周波数が変化する送信信号、又は、パルスの送信信号を発生させる。時間の経過に伴って周波数が変化する送信信号としては、例えば、アップチャープの信号、又は、ダウンチャープの信号がある。
信号発生器１２ａは、送信信号を出力先選択部１２ｂに出力する。
出力先選択部１２ｂは、Ｎ個のサーキュレータ１３－１～１３－Ｎの中で、次に送信波を放射させる順番のアンテナ１１－ｎと接続されているサーキュレータ１３－ｎに対して、信号発生器１２ａにより発生された送信信号を出力する。

サーキュレータ１３－ｎ（ｎ＝１，・・・，Ｎ）は、出力先選択部１２ｂから出力された送信信号をアンテナ１１－ｎに出力する。
また、サーキュレータ１３－ｎは、アンテナ１１－ｎから出力された受信信号を信号受信部１４－ｎに出力する。

信号受信部１４－ｎ（ｎ＝１，・・・，Ｎ）は、サーキュレータ１３－ｎから出力された受信信号に対する受信処理を実施する。受信処理としては、例えば、受信信号の周波数をダウンコンバートする処理のほか、周波数変換後の受信信号をアナログ信号からデジタル信号に変換する処理がある。
信号受信部１４－ｎは、デジタル信号である受信データＳ（ｔ，ｇ，ｈ，ｃ）をバイタル測定装置２０に出力する。
ｔは、アンテナ１１－ｎによる反射波の受信時刻である。ｇは、送信波を放射したアンテナ１１－ｎを識別する変数であり、ｇ＝１，・・・，Ｎである。ｈは、反射波を受信したアンテナ１１－ｎを識別する変数であり、ｈ＝１，・・・，Ｎである。ｃは、アンテナ１１－１～１１－Ｎによる送信波の送信サイクルを識別する変数であり、ｃ＝１，・・・，Ｃである。Ｃは、２以上の整数である。

バイタル測定装置２０は、信号取得部２１、フーリエ変換部２２、マップ算出部２３及びバイタル推定部２４を備えている。
信号取得部２１は、例えば、図２に示す信号取得回路３１によって実現される。
信号取得部２１は、信号受信部１４－１～１４－Ｎのそれぞれから、反射波の受信信号として、受信データＳ（ｔ，ｇ，ｈ，ｃ）を取得し、受信データＳ（ｔ，ｇ，ｈ，ｃ）をフーリエ変換部２２に出力する。

フーリエ変換部２２は、例えば、図２に示すフーリエ変換回路３２によって実現される。
フーリエ変換部２２は、信号取得部２１によって、信号受信部１４－１～１４－Ｎからの受信データＳ（ｔ，ｇ，ｈ，ｃ）が取得される毎に、それぞれの受信データＳ（ｔ，ｇ，ｈ，ｃ）を時間方向にフーリエ変換する。フーリエ変換としては、例えば、高速フーリエ変換（ＦＦＴ：Ｆａｓｔ Ｆｏｕｒｉｅｒ Ｔｒａｎｓｆｏｒｍ）、又は、離散フーリエ変換（ＤＦＴ：Ｄｉｓｃｒｅｔｅ Ｆｏｕｒｉｅｒ Ｔｒａｎｓｆｏｒｍ）がある。
フーリエ変換部２２は、それぞれのフーリエ変換後の信号Ｓ（ｒ，ｇ，ｈ，ｃ）をマップ算出部２３に出力する。ｒは、バイタル測定装置２０からの距離ビンを識別する変数である。ｒ＝１，・・・，Ｒである。Ｒは、２以上の整数である。

マップ算出部２３は、例えば、図２に示すマップ算出回路３３によって実現される。
マップ算出部２３は、フーリエ変換部２２によって、それぞれの受信データＳ（ｔ，ｇ，ｈ，ｃ）がフーリエ変換される毎に、それぞれのフーリエ変換後の信号Ｓ（ｒ，ｇ，ｈ，ｃ）を用いて、バイタル測定装置２０からのそれぞれの距離ビンｒ（ｒ＝１，・・・，Ｒ）に対応する複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の２次元方位マップを算出する。
マップ算出部２３は、それぞれの距離ビンｒに対応する複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の２次元方位マップをバイタル推定部２４に出力する。

バイタル推定部２４は、例えば、図２に示すバイタル推定回路３４によって実現される。
バイタル推定部２４は、存在位置特定部２５及びバイタル推定処理部２６を備えている。
バイタル推定部２４は、マップ算出部２３によって、それぞれの距離ビンｒ（ｒ＝１，・・・，Ｒ）に対応する複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の２次元方位マップが算出される毎に、２次元方位マップを取得する。
バイタル推定部２４は、それぞれの距離ビンｒに対応する複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の２次元方位マップの時間変化に基づいて、対象物体に含まれる被測定者ｋ（ｋ＝１，・・・，Ｋ）が存在している位置を特定する。
バイタル推定部２４は、それぞれの被測定者ｋが存在している位置についての複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の位相の時間変化から、それぞれの被測定者ｋのバイタルを推定する。

存在位置特定部２５は、マップ算出部２３によって、それぞれの距離ビンｒ（ｒ＝１，・・・，Ｒ）に対応する複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の２次元方位マップが算出される毎に、２次元方位マップを取得する。
存在位置特定部２５は、それぞれの距離ビンｒに対応する複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の２次元方位マップの時間変化に基づいて、それぞれの被測定者ｋ（ｋ＝１，・・・，Ｋ）が存在している位置を特定する。
存在位置特定部２５は、それぞれの被測定者ｋが存在している位置についての複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の位相の時間変化を示す位相変化信号θ（ｒ，ｃ）をバイタル推定処理部２６に出力する。

バイタル推定処理部２６は、呼吸数推定部２７及び心拍数推定部２８を備えている。
バイタル推定処理部２６は、存在位置特定部２５から出力されたそれぞれの被測定者ｋ（ｋ＝１，・・・，Ｋ）についての位相変化信号θ（ｒ，ｃ）に基づいて、それぞれの被測定者ｋのバイタルを推定する。
呼吸数推定部２７は、それぞれの被測定者ｋについての位相変化信号θ（ｒ，ｃ）に基づいて、それぞれの被測定者ｋの呼吸数ＲＲを推定する。
心拍数推定部２８は、それぞれの被測定者ｋについての位相変化信号θ（ｒ，ｃ）に基づいて、それぞれの被測定者ｋの心拍数ＨＲを推定する。

図１では、バイタル測定装置２０の構成要素である信号取得部２１、フーリエ変換部２２、マップ算出部２３及びバイタル推定部２４のそれぞれが、図２に示すような専用のハードウェアによって実現されるものを想定している。即ち、バイタル測定装置２０が、信号取得回路３１、フーリエ変換回路３２、マップ算出回路３３及びバイタル推定回路３４によって実現されるものを想定している。
信号取得回路３１、フーリエ変換回路３２、マップ算出回路３３及びバイタル推定回路３４のそれぞれは、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）、ＦＰＧＡ（Ｆｉｅｌｄ－Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）、又は、これらを組み合わせたものが該当する。

バイタル測定装置２０の構成要素は、専用のハードウェアによって実現されるものに限るものではなく、バイタル測定装置２０が、ソフトウェア、ファームウェア、又は、ソフトウェアとファームウェアとの組み合わせによって実現されるものであってもよい。
ソフトウェア又はファームウェアは、プログラムとして、コンピュータのメモリに格納される。コンピュータは、プログラムを実行するハードウェアを意味し、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、中央処理装置、処理装置、演算装置、マイクロプロセッサ、マイクロコンピュータ、プロセッサ、あるいは、ＤＳＰ（Ｄｉｇｉｔａｌ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｏｒ）が該当する。

図３は、バイタル測定装置２０が、ソフトウェア又はファームウェア等によって実現される場合のコンピュータのハードウェア構成図である。
バイタル測定装置２０が、ソフトウェア又はファームウェア等によって実現される場合、信号取得部２１、フーリエ変換部２２、マップ算出部２３及びバイタル推定部２４におけるそれぞれの処理手順をコンピュータに実行させるためのプログラムがメモリ４１に格納される。そして、コンピュータのプロセッサ４２がメモリ４１に格納されているプログラムを実行する。

また、図２では、バイタル測定装置２０の構成要素のそれぞれが専用のハードウェアによって実現される例を示し、図３では、バイタル測定装置２０がソフトウェア又はファームウェア等によって実現される例を示している。しかし、これは一例に過ぎず、バイタル測定装置２０における一部の構成要素が専用のハードウェアによって実現され、残りの構成要素がソフトウェア又はファームウェア等によって実現されるものであってもよい。

次に、図１に示すバイタル測定システムの動作について説明する。
図４は、バイタル測定装置２０の処理手順であるバイタル測定方法を示すフローチャートである。
信号送信部１２の信号発生器１２ａは、送信信号として、例えば、周波数変調連続波（ＦＭＣＷ：Ｆｒｅｑｕｅｎｃｙ Ｍｏｄｕｌａｔｅｄ Ｃｏｎｔｉｎｕｏｕｓ Ｗａｖｅ）方式、又は、高速チャープ変調（ＦＣＭ：Ｆａｓｔ－Ｃｈｉｒｐ Ｍｏｄｕｌａｔｉｏｎ）方式に従って、アップチャープの信号、又は、ダウンチャープの信号を生成する。
図１に示すバイタル測定システムでは、信号発生器１２ａが、図５に示すように、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、アンテナ１１－１～１１－Ｎの数分だけ、アップチャープの信号Ｔｘ（ｎ）を繰り返し生成している。したがって、信号発生器１２ａは、全部で、Ｎ×Ｃ個のアップチャープの信号Ｔｘ（ｎ）を生成している。
図５は、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、信号発生器１２ａにより生成されるアップチャープの信号Ｔｘ（１）～Ｔｘ（Ｎ）を示す説明図である。図５の例では、Ｎ＝３である。
信号発生器１２ａは、生成したアップチャープの信号Ｔｘ（ｎ）（ｎ＝１，・・・，Ｎ）の周波数を高周波数帯の信号にアップコンバートし、アップコンバート後の信号を送信信号Ｔｘ’（ｎ）として、出力先選択部１２ｂに出力する。高周波数帯としては、例えば、３０～３００ＧＨｚ程度のミリ波帯である。

出力先選択部１２ｂは、信号発生器１２ａから、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、送信信号Ｔｘ’（ｎ）を取得する。
出力先選択部１２ｂは、Ｎ個のサーキュレータ１３－１～１３－Ｎの中で、次に送信波を放射する順番のアンテナ１１－ｎと接続されているサーキュレータ１３－ｎに対して、送信信号Ｔｘ’（ｎ）を出力する。
出力先選択部１２ｂは、例えば、次に送信波を放射する順番のアンテナ１１－ｎがアンテナ１１－１であれば、送信信号Ｔｘ’（１）をサーキュレータ１３－１に出力し、次に送信波を放射する順番のアンテナ１１－ｎがアンテナ１１－２であれば、送信信号Ｔｘ’（２）をサーキュレータ１３－２に出力する。
また、出力先選択部１２ｂは、例えば、次に送信波を放射する順番のアンテナ１１－ｎがアンテナ１１－Ｎであれば、送信信号Ｔｘ’（Ｎ）をサーキュレータ１３－Ｎに出力する。

サーキュレータ１３－ｎ（ｎ＝１，・・・，Ｎ）は、出力先選択部１２ｂから送信信号Ｔｘ’（ｎ）を受けると、送信信号Ｔｘ’（ｎ）をアンテナ１１－ｎに出力する。
アンテナ１１－ｎ（ｎ＝１，・・・，Ｎ）は、サーキュレータ１３－ｎから送信信号Ｔｘ’（ｎ）を受けると、送信信号Ｔｘ’（ｎ）に係る電波である送信波を対象物体が存在している空間に放射する。アンテナ１１－ｎから放射された送信波は、対象物体によって反射される。即ち、アンテナ１１－ｎから放射された送信波は、それぞれの被測定者ｋによって反射されるほか、空間を形成している壁等によって反射される。
それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、Ｎ個のアンテナ１１－１～１１－Ｎの中の１つのアンテナが、送信波を順番にＮ回空間に放射するため、アンテナ１１－１～１１－Ｎのそれぞれは、それぞれの送信サイクルｃにおいて、反射波をＮ回受信する。
アンテナ１１－ｎ（ｎ＝１，・・・，Ｎ）は、それぞれの反射波の受信信号をサーキュレータ１３－ｎに出力する。
サーキュレータ１３－ｎ（ｎ＝１，・・・，Ｎ）は、アンテナ１１－ｎから出力された受信信号を信号受信部１４－ｎに出力する。

信号受信部１４－ｎ（ｎ＝１，・・・，Ｎ）は、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、サーキュレータ１３－ｎから出力されたＮ個の受信信号のそれぞれに対する受信処理を実施する。
信号受信部１４－ｎは、受信信号に対する受信処理として、例えば、受信信号の周波数を中間周波数帯の周波数にダウンコンバートする処理のほか、周波数変換後の受信信号をアナログ信号からデジタル信号に変換する処理を行う。
信号受信部１４－ｎは、それぞれの送信サイクルｃにおいて、Ｎ個のデジタル信号である受信データＳ（ｔ，ｇ，ｈ，ｃ）のそれぞれをバイタル測定装置２０に出力する。
これにより、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、信号受信部１４－１～１４－Ｎから、全部でＮ×Ｎ（＝Ｇ×Ｈ）個の受信データＳ（ｔ，ｇ，ｈ，ｃ）がバイタル測定装置２０に与えられる。

バイタル測定装置２０の信号取得部２１は、信号受信部１４－１～１４－Ｎから、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、Ｎ×Ｎ個の受信データＳ（ｔ，ｇ，ｈ，ｃ）を取得する（図４のステップＳＴ１）。
信号取得部２１は、それぞれの送信サイクルｃにおいて、Ｎ×Ｎ個の受信データＳ（ｔ，ｇ，ｈ，ｃ）をフーリエ変換部２２に出力する。

フーリエ変換部２２は、信号取得部２１から、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、Ｎ×Ｎ個の受信データＳ（ｔ，ｇ，ｈ，ｃ）を取得する。
フーリエ変換部２２は、それぞれの送信サイクルｃにおいて、Ｎ×Ｎ個の受信データＳ（ｔ，ｇ，ｈ，ｃ）の中から、反射波を受信したそれぞれのアンテナ１１－ｎ（ｎ＝１，・・・，Ｎ）に係るＮ個の受信データＳ（ｔ，ｇ，ｈ＝ｎ，ｃ）の取り出しを行う。
フーリエ変換部２２は、それぞれの送信サイクルｃにおいて、それぞれのアンテナ１１－ｎに係るＮ個の受信データＳ（ｔ，ｇ，ｈ＝ｎ，ｃ）のそれぞれを時間方向にフーリエ変換する（図４のステップＳＴ２）。
それぞれの送信サイクルｃにおいて、それぞれのアンテナ１１－ｎに係るＮ個の受信データＳ（ｔ，ｇ，ｈ＝ｎ，ｃ）のそれぞれが、フーリエ変換部２２によってフーリエ変換されることで、Ｎ個のフーリエ変換後の信号Ｓ（ｒ，ｇ，ｈ，ｃ）が生成される。
フーリエ変換部２２によるフーリエ変換後の信号Ｓ（ｒ，ｇ，ｈ，ｃ）は、図６に示すように、距離ビンｒに対応する複素電力を示す信号である。
図６は、或る送信サイクルｃにおけるＮ個のフーリエ変換後の信号Ｓ（ｒ，ｇ，ｈ，ｃ）を示す説明図である。
図６において、横軸は、距離ビンを示し、縦軸は、複素電力を示している。
フーリエ変換部２２は、それぞれの送信サイクルｃにおいて、Ｎ個のフーリエ変換後の信号Ｓ（ｒ，ｇ，ｈ，ｃ）をマップ算出部２３に出力する。

マップ算出部２３は、フーリエ変換部２２から、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、Ｎ個のフーリエ変換後の信号Ｓ（ｒ，ｇ，ｈ，ｃ）を取得する。
マップ算出部２３は、それぞれの送信サイクルｃにおいて、図７に示すように、それぞれのフーリエ変換後の信号Ｓ（ｒ，ｇ，ｈ，ｃ）から、それぞれの距離ビンｒ（ｒ＝１，・・・，Ｒ）に対応する複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の２次元方位マップを算出する（図４のステップＳＴ３）。
図７は、それぞれの距離ビンｒ（ｒ＝１，・・・，Ｒ）に対応する複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の２次元方位マップを示す説明図である。
図７において、横軸は、アジマス方向を示し、縦軸は、エレベーション方向を示している。
具体的には、マップ算出部２３は、それぞれのフーリエ変換後の信号Ｓ（ｒ，ｇ，ｈ，ｃ）について、距離ビンｒ毎に、ＤＢＦ（Ｄｉｇｉｔａｌ Ｂｅａｍ Ｆｏｒｍｉｎｇ）による２次元測角処理を行う。マップ算出部２３が２次元測角処理を行うことで、２次元方位（ａｚ，ｅｌ）についての複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）を示す２次元方位マップが得られる。ここでの２次元測角処理は、アジマス方向とエレベーション方向との２次元測角処理である。また、マップ算出部２３は、測角手法として、ＤＢＦを用いている。しかし、これは一例に過ぎず、マップ算出部２３は、Ｃａｐｏｎ法等の他の測角手法を用いて、２次元測角処理を行うようにしてもよい。なお、２次元測角処理自体は、公知の技術であるため詳細な説明を省略する。
それぞれの送信サイクルｃにおいて、マップ算出部２３が、Ｎ個のフーリエ変換後の信号Ｓ（ｒ，ｇ，ｈ，ｃ）を取得しているため、距離ビンｒの数が例えばＲ＝１０００個であれば、Ｎ個のフーリエ変換後の信号Ｓ（ｒ，ｇ，ｈ，ｃ）から、Ｎ×１０００個の２次元方位マップが算出される。
マップ算出部２３は、それぞれの距離ビンｒに対応するＮ個の２次元方位マップの合成処理を行う。それぞれの距離ビンｒに対応するＮ個の２次元方位マップの合成処理としては、例えば、Ｎ個の２次元方位マップに含まれている同一の２次元方位同士の複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の加算処理、あるいは、同一の２次元方位同士の複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の平均処理がある。
マップ算出部２３は、それぞれの送信サイクルｃにおいて、それぞれの距離ビンｒに対応する合成処理後の２次元方位マップをバイタル推定部２４に出力する。

バイタル推定部２４の存在位置特定部２５は、マップ算出部２３から、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、それぞれの距離ビンｒに対応する合成処理後の２次元方位マップを取得する。
存在位置特定部２５は、それぞれの距離ビンｒに対応する合成処理後の２次元方位マップの時間変化に基づいて、対象物体に含まれる被測定者ｋ（ｋ＝１，・・・，Ｋ）が存在している位置を特定する（図４のステップＳＴ４）。
以下、存在位置特定部２５による位置の特定処理を具体的に説明する。

被測定者ｋ（ｋ＝１，・・・，Ｋ）による反射波には、被測定者ｋの呼吸を示す呼吸信号と、被測定者ｋの心拍を示す心拍信号とが重畳されているものの、複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の大部分は、呼吸信号の電力である。呼吸信号は、被測定者ｋの胸部における往復運動によって生じる位相変動を示す信号である。このため、同一の距離ビンｒに属し、かつ、同一の２次元方位（ａｚ，ｅｌ）に属している複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）は、複素信号空間において、図８に示すように、時間の経過に伴って円周上の往復運動を行う。
図８は、同一の距離ビンｒに属し、かつ、同一の２次元方位（ａｚ，ｅｌ）に属している複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）における位相の時間変化を示す説明図である。図８では、複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）が、複素信号空間において、円周上の往復運動を行っている。
図８において、●は、送信サイクルｃ＝１～Ｃにおける複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の位相θ（ｒ，ａｚ，ｅｌ，ｃ）を示し、Ｃ個の●がプロットされている。

存在位置特定部２５は、それぞれの距離ビンｒ（ｒ＝１，・・・，Ｒ）に対応する２次元方位（ａｚ，ｅｌ）についての複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）における位相の時間変化が描くフィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）を算出する。
具体的には、存在位置特定部２５は、送信サイクルｃ（ｃ＝１，・・・，Ｃ）の複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）を示す複素データがｓ_ｃであるとして、以下の式（１）に示す評価関数Ｊ（α，β）が最小になるα，βを求めることで、最適なフィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）を算出する。αは、フィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）の中心であり、βは、フィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）の半径である。ｓ_ｃの下付きのｃは、送信サイクルを示す変数である。
存在位置特定部２５によるフィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）の算出数は、Ｒ×ＡＺ×ＥＬである。ＡＺは、複素データｓ_ｃにおけるアジマス方向の分解能に対応するセル数、ＥＬは、複素データｓ_ｃにおけるエレベーション方向の分解能に対応するセル数である。

評価関数Ｊ（α，β）は、式変形によって、以下の式（２）に示すような行列形式で表すことができる。

式（２）において、＊は、複素共役を表す数学記号である。

存在位置特定部２５は、以下の式（３）に示すように、最小２乗法を用いて、式（２）に示すパラメータベクトルｐを推定する。

式（３）において、＋は、疑似逆行列を表す数学記号である。ｐの文字の上に表示されている“～”は、パラメータベクトルｐの推定結果であることを示す記号である。

パラメータベクトルｐの推定結果が、以下の式（４）のように表されるとすれば、存在位置特定部２５は、以下の式（５）及び式（６）に示すように、フィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）の中心αと、フィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）の半径βとを求めることができる。フィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）の中心αと、フィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）の半径βとが求まれば、フィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）が算出されたことになる。

式（５）及び式（６）において、Ｒｅ（□）は、複素数である□の実部であることを表し、Ｉｍ（□）は、複素数である□の虚部であることを表している。
ここでは、存在位置特定部２５が、Ｋａｓａ ｆｉｔと呼ばれる円フィッティング手法を用いて、フィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）を算出している。しかし、これは一例に過ぎず、存在位置特定部２５が、Ｐｒａｔｔ ｆｉｔ、Ｔａｕｂｉｎ ｆｉｔ、又は、Ｈｙｐｅｒ ｆｉｔ等の円フィッティング手法を用いて、フィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）を算出するようにしてもよい。

次に、存在位置特定部２５は、以下の式（７）に示すように、フィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）の誤差であるフィッティング誤差Ｅ（ａｚ，ｅｌ）を算出する。

式（７）において、Ｓ（ｒ，ａｚ，ｅｌ，ｃ）は、２次元方位スペクトルであり、２次元方位スペクトルは、送信サイクルｃ（ｃ＝１，・・・，Ｃ）における、距離ビンｒに対応する２次元方位（ａｚ，ｅｌ）についての複素電力ＣＰ（ｒ，ａｚ，ｅｌ）に相当する。
α（ｒ，ａｚ，ｅｌ）は、フィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）の中心であり、β（ｒ，ａｚ，ｅｌ）は、フィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）の半径である。

被測定者ｋの呼吸を示す呼吸信号と、被測定者ｋの心拍を示す心拍信号とが反射波に重畳されていれば、フィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）の半径βが大きくなり、かつ、フィッティング誤差Ｅ（ａｚ，ｅｌ）が小さくなる。一方、反射波が、壁等の静止物による反射波であれば、フィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）の半径βが小さくなり、かつ、フィッティング誤差Ｅ（ａｚ，ｅｌ）が大きくなる。
存在位置特定部２５は、算出したＲ×ＡＺ×ＥＬ個のフィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）の半径βと、第１の閾値Ｔｈ_１とを比較する。第１の閾値Ｔｈ_１は、存在位置特定部２５の内部メモリに格納されていてもよいし、バイタル測定装置２０の外部から与えられるものであってもよい。
存在位置特定部２５は、算出したＲ×ＡＺ×ＥＬ個のフィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）の中で、半径βが第１の閾値Ｔｈ_１以上となるフィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）を探索する。

次に、存在位置特定部２５は、半径βが第１の閾値Ｔｈ_１以上となる１つ以上のフィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）のフィッティング誤差Ｅ（ａｚ，ｅｌ）と、第２の閾値Ｔｈ_２とを比較する。第２の閾値Ｔｈ_２は、存在位置特定部２５の内部メモリに格納されていてもよいし、バイタル測定装置２０の外部から与えられるものであってもよい。
存在位置特定部２５は、半径βが第１の閾値Ｔｈ_１以上となる１つ以上のフィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）の中で、フィッティング誤差Ｅ（ａｚ，ｅｌ）が第２の閾値Ｔｈ_２以下となるフィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）を探索する。

存在位置特定部２５は、それぞれの被測定者ｋ（ｋ＝１，・・・，Ｋ）が存在している位置として、半径βが第１の閾値Ｔｈ_１以上となり、かつ、フィッティング誤差Ｅ（ａｚ，ｅｌ）が第２の閾値Ｔｈ_２以下となるフィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）に係る距離ビンｒ、アジマス方向ａｚ及びエレベーション方向ｅｌのそれぞれを特定する。
存在位置特定部２５は、送信サイクルｃ＝１～Ｃにおける、それぞれの被測定者ｋが存在している位置についての複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）から、図９に示すような、複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）における位相の時間変化を示す位相変化信号θ（ｒ，ｃ）を算出する。
存在位置特定部２５は、被測定者ｋが存在している位置についての複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）から位相変化信号θ（ｒ，ｃ）を算出する方法として、例えば、ＡＤ（Ａｒｃｔａｎｇｅｎｔ Ｄｅｍｏｄｕｌａｔｉｏｎ）法、又は、ＣＳＤ（Ｃｏｍｐｌｅｘ Ｓｉｇｎａｌ Ｄｅｍｏｄｕｌａｔｉｏｎ）法を用いることができる。

図９は、存在位置特定部２５により算出される位相変化信号θ（ｒ，ｃ）を示す説明図である。
図９において、横軸は、送信サイクルｃに対応する時間、縦軸は、距離ビンｒにおける複素電力ＣＰ（ｒ，ａｚ，ｅｌ，ｃ）の位相［ｒａｄ］を示している。
存在位置特定部２５は、それぞれの被測定者ｋについての位相変化信号θ（ｒ，ｃ）を呼吸数推定部２７及び心拍数推定部２８のそれぞれに出力する。

呼吸数推定部２７は、存在位置特定部２５から、それぞれの被測定者ｋ（ｋ＝１，・・・，Ｋ）についての位相変化信号θ（ｒ，ｃ）を取得する。
呼吸数推定部２７は、それぞれの被測定者ｋについての位相変化信号θ（ｒ，ｃ）をフーリエ変換することで、それぞれの被測定者ｋの呼吸数ＲＲを推定する（図４のステップＳＴ５）。
以下、呼吸数推定部２７による呼吸数ＲＲの推定処理を具体的に説明する。

被測定者ｋの呼吸を示す呼吸信号の波形は、概ね正弦波である。このため、位相変化信号θ（ｒ，ｃ）の中の正弦波的な変動は、被測定者ｋの呼吸によるものである。
呼吸数推定部２７は、それぞれの被測定者ｋについての位相変化信号θ（ｒ，ｃ）をスロータイム方向にフーリエ変換することで、図１０に示すような呼吸スペクトルＳ_ＲＲ（ｒ，ｓｆ）を得る。スロータイムは、送信波の送信時刻である。呼吸スペクトルＳ_ＲＲ（ｒ，ｓｆ）は、位相変化信号θ（ｒ，ｃ）のフーリエ変換結果であり、ｓｆは、スロータイム方向の周波数である。
図１０は、呼吸スペクトルＳ_ＲＲ（ｒ，ｓｆ）の一例を示す説明図である。
図１０において、横軸は、呼吸数ＲＲ［ｂｐｍ］、縦軸は、呼吸スペクトル［ｄＢ］である。

呼吸数推定部２７は、それぞれの被測定者ｋについての呼吸スペクトルＳ_ＲＲ（ｒ，ｓｆ）の中で、第３の閾値Ｔｈ_３以上の呼吸スペクトルＳ_ＲＲを探索する。第３の閾値Ｔｈ_３は、呼吸数推定部２７の内部メモリに格納されていてもよいし、バイタル測定装置２０の外部から与えられるものであってもよい。呼吸スペクトルＳ_ＲＲ（ｒ，ｓｆ）の中には、第３の閾値Ｔｈ_３以上になる呼吸スペクトルＳ_ＲＲが１つ以上存在している。
呼吸数推定部２７は、第３の閾値Ｔｈ_３以上の呼吸スペクトルＳ_ＲＲに対応する呼吸数ＲＲの中で、第４の閾値Ｔｈ_４以上の呼吸数ＲＲを探索する。第４の閾値Ｔｈ_４は、呼吸数推定部２７の内部メモリに格納されていてもよいし、バイタル測定装置２０の外部から与えられるものであってもよい。
呼吸数推定部２７は、第３の閾値Ｔｈ_３以上の呼吸スペクトルＳ_ＲＲに対応する呼吸数ＲＲの中で、第４の閾値Ｔｈ_４以上の呼吸数ＲＲが、被測定者ｋの呼吸数ＲＲであると推定する。

ここでは、呼吸数推定部２７が、呼吸スペクトルＳ_ＲＲ（ｒ，ｓｆ）の中で、第３の閾値Ｔｈ_３以上の呼吸スペクトルＳ_ＲＲであり、かつ、当該呼吸スペクトルＳ_ＲＲに対応する呼吸数ＲＲが第４の閾値Ｔｈ_４以上であれば、当該呼吸数ＲＲが被測定者ｋの呼吸数ＲＲであると推定している。しかし、これは一例に過ぎず、呼吸数推定部２７が、第３の閾値Ｔｈ_３以上の呼吸スペクトルＳ_ＲＲの中で、最大の呼吸スペクトルＳ_ＲＲに対応する呼吸数ＲＲが、被測定者ｋの呼吸数ＲＲであると推定するようにしてもよい。

心拍数推定部２８は、存在位置特定部２５から、それぞれの被測定者ｋ（ｋ＝１，・・・，Ｋ）についての位相変化信号θ（ｒ，ｃ）を取得する。
心拍数推定部２８は、それぞれの被測定者ｋについての位相変化信号θ（ｒ，ｃ）に基づいて、被測定者ｋの心拍数ＨＲを推定する（図４のステップＳＴ６）。
以下、心拍数推定部２８による心拍数ＨＲの推定処理を具体的に説明する。

まず、心拍数推定部２８は、それぞれの被測定者ｋについての位相変化信号θ（ｒ，ｃ）を図示せぬハイパスフィルタ（ＨＰＦ：Ｈｉｇｈ－Ｐａｓｓ Ｆｉｌｔｅｒ）に通すことで、位相変化信号θ（ｒ，ｃ）に重畳されている呼吸信号を除去する。
次に、心拍数推定部２８は、それぞれの呼吸信号除去後の位相変化信号θ’（ｒ，ｃ）を連続ウェーブレット変換することで、図１１に示すようなスカログラムを得る。スカログラムは、呼吸信号除去後の位相変化信号θ’（ｒ，ｃ）の連続ウェーブレット変換結果である。
図１１は、スカログラムの一例を示す説明図である。
図１１において、横軸は、送信サイクルｃであり、縦軸は、周波数［Ｈｚ］である。
スカログラムには、図１１に示すように、心拍信号が生じている位置に縞状が現れる。

心拍数推定部２８は、スカログラムの絶対値をスロータイム方向にフーリエ変換することで、図１２に示すような、スロータイム方向の２次元スペクトルＷ（ｒ，ｆ，ｓｆ）を得る。ｆは、周波数成分である。
図１２は、スロータイム方向の２次元スペクトルＷ（ｒ，ｆ，ｓｆ）の一例を示す説明図である。
図１２において、横軸は、心拍数ＨＲ［ｂｐｍ］、縦軸は、周波数［Ｈｚ］である。
２次元スペクトルＷ（ｒ，ｆ，ｓｆ）には、図１２に示すように、複数の周波数成分ｆに心拍信号が分散している。心拍信号が分散している高周波領域のスペクトルを最大比合成することで、心拍信号が強調されている心拍スペクトルＳ_ＨＲ（ｒ，ｓｆ）が得られる。

心拍数推定部２８は、以下の式（８）に示すように、それぞれの被測定者ｋについての２次元スペクトルＷ（ｒ，ｆ，ｓｆ）に含まれている高周波成分を示す行列Ｕを定義する。

（９）
式（８）において、ｆ_Ｌは、高周波成分の下限周波数、ｆ_Ｈは、高周波成分の上限周波数である。

次に、心拍数推定部２８は、図１３に示すように、行列Ｕの相関行列Ｒ_ｘｘ＝ＵＵ^Ｈを生成し、相関行列Ｒ_ｘｘを固有値分解することで、第１の固有ベクトルｕ_１を得る。Ｈは、エルミート転置を示す数学記号である。
心拍数推定部２８は、以下の式（１０）に示すように、第１の固有ベクトルｕ_１が示す方向に、行列Ｕを射影することで最大比合成を行い、図１４に示すような心拍スペクトルＳ_ＨＲ（ｒ，ｓｆ）を得る。

図１３は、心拍数推定部２８による最大比合成処理を示す説明図である。
図１３において、横軸は、心拍数ＨＲ［ｂｐｍ］、縦軸は、周波数［Ｈｚ］である。
図１４は、心拍数推定部２８により得られる心拍スペクトルＳ_ＨＲ（ｒ，ｓｆ）の一例を示す説明図である。
図１４において、横軸は、心拍数ＨＲ［ｂｐｍ］、縦軸は、心拍スペクトル［ｄＢ］である。

心拍数推定部２８は、それぞれの被測定者ｋについての心拍スペクトルＳ_ＨＲ（ｒ，ｓｆ）の中で、第５の閾値Ｔｈ_５以上の心拍スペクトルＳ_ＨＲを探索する。第５の閾値Ｔｈ_５は、心拍数推定部２８の内部メモリに格納されていてもよいし、バイタル測定装置２０の外部から与えられるものであってもよい。被測定者ｋが存在していれば、心拍スペクトルＳ_ＨＲが第５の閾値Ｔｈ_５以上になる。
心拍数推定部２８は、第５の閾値Ｔｈ_５以上の心拍スペクトルＳ_ＨＲに対応する心拍数ＨＲが、被測定者ｋの心拍数ＨＲであると推定する。

以上の実施の形態１では、対象物体による反射波を受信するアンテナ１１－１～１１－Ｎから、反射波の受信信号を取得する信号取得部２１と、信号取得部２１により取得された受信信号を時間方向にフーリエ変換するフーリエ変換部２２とを備えるように、バイタル測定装置２０を構成した。また、バイタル測定装置２０は、フーリエ変換部２２によるフーリエ変換後の信号を用いて、バイタル測定装置２０からの複数の距離ビンのそれぞれに対応する複素電力の２次元方位マップを算出するマップ算出部２３と、マップ算出部２３により算出されたそれぞれの距離ビンに対応する複素電力の２次元方位マップの時間変化に基づいて、対象物体に含まれる被測定者が存在している位置を特定し、被測定者が存在している位置についての複素電力の位相の時間変化から、被測定者のバイタルを推定するバイタル推定部２４とを備えている。したがって、バイタル測定装置２０は、被測定者が複数の場合であっても、それぞれの被測定者のバイタルを検出することができる。

図１に示すバイタル測定装置２０では、存在位置特定部２５が、それぞれの被測定者ｋが存在している位置として、半径βが第１の閾値Ｔｈ_１以上となり、かつ、フィッティング誤差Ｅ（ａｚ，ｅｌ）が第２の閾値Ｔｈ_２以下となるフィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）に係る距離ビンｒ、アジマス方向ａｚ及びエレベーション方向ｅｌのそれぞれを特定している。しかし、これは一例に過ぎず、存在位置特定部２５が、それぞれの被測定者ｋが存在している位置として、半径βが第１の閾値Ｔｈ_１以上となるフィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）に係る距離ビンｒ、アジマス方向ａｚ及びエレベーション方向ｅｌのそれぞれを特定するようにしてもよい。
また、存在位置特定部２５が、それぞれの被測定者ｋが存在している位置として、フィッティング誤差Ｅ（ａｚ，ｅｌ）が第２の閾値Ｔｈ_２以下となるフィッティング円Ｃｉｒｃ（ｒ，ａｚ，ｅｌ）に係る距離ビンｒ、アジマス方向ａｚ及びエレベーション方向ｅｌのそれぞれを特定するようにしてもよい。

実施の形態２．
実施の形態２では、信号取得部２１により取得された受信信号に含まれている、移動体による反射波に係る信号を抑圧する信号抑圧部２９を備えているバイタル測定装置２０について説明する。

図１５は、実施の形態２に係るバイタル測定装置２０を含むバイタル測定システムを示す構成図である。図１５において、図１と同一符号は同一又は相当部分を示すので説明を省略する。
図１６は、実施の形態２に係るバイタル測定装置２０のハードウェアを示すハードウェア構成図である。図１６において、図２と同一符号は同一又は相当部分を示すので説明を省略する。

図１５に示すバイタル測定装置２０は、信号取得部２１、信号抑圧部２９、フーリエ変換部２２、マップ算出部２３及びバイタル推定部２４を備えている。
信号抑圧部２９は、例えば、図１６に示す信号抑圧回路３５によって実現される。
信号抑圧部２９は、信号取得部２１によって、信号受信部１４－１～１４－Ｎからの受信データＳ（ｔ，ｇ，ｈ，ｃ）が取得される毎に、それぞれの受信データＳ（ｔ，ｇ，ｈ，ｃ）に含まれている、移動体による反射波に係る信号を抑圧する。
信号抑圧部２９は、それぞれの信号抑圧後の受信データＳ’（ｔ，ｇ，ｈ，ｃ）をフーリエ変換部２２に出力する。

図１５では、バイタル測定装置２０の構成要素である信号取得部２１、信号抑圧部２９、フーリエ変換部２２、マップ算出部２３及びバイタル推定部２４のそれぞれが、図１６に示すような専用のハードウェアによって実現されるものを想定している。即ち、バイタル測定装置２０が、信号取得回路３１、信号抑圧回路３５、フーリエ変換回路３２、マップ算出回路３３及びバイタル推定回路３４によって実現されるものを想定している。
信号取得回路３１、信号抑圧回路３５、フーリエ変換回路３２、マップ算出回路３３及びバイタル推定回路３４のそれぞれは、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ＡＳＩＣ、ＦＰＧＡ、又は、これらを組み合わせたものが該当する。

バイタル測定装置２０の構成要素は、専用のハードウェアによって実現されるものに限るものではなく、バイタル測定装置２０が、ソフトウェア、ファームウェア、又は、ソフトウェアとファームウェアとの組み合わせによって実現されるものであってもよい。
バイタル測定装置２０が、ソフトウェア又はファームウェア等によって実現される場合、信号取得部２１、信号抑圧部２９、フーリエ変換部２２、マップ算出部２３及びバイタル推定部２４におけるそれぞれの処理手順をコンピュータに実行させるためのプログラムが図３に示すメモリ４１に格納される。そして、図３に示すプロセッサ４２がメモリ４１に格納されているプログラムを実行する。

また、図１６では、バイタル測定装置２０の構成要素のそれぞれが専用のハードウェアによって実現される例を示し、図３では、バイタル測定装置２０がソフトウェア又はファームウェア等によって実現される例を示している。しかし、これは一例に過ぎず、バイタル測定装置２０における一部の構成要素が専用のハードウェアによって実現され、残りの構成要素がソフトウェア又はファームウェア等によって実現されるものであってもよい。

次に、図１５に示すバイタル測定システムの動作について説明する。
図１に示すバイタル測定システムでは、信号発生器１２ａが、図５に示すように、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、アップチャープの信号Ｔｘ（１）～Ｔｘ（Ｎ）を生成している。
図１５に示すバイタル測定システムでは、信号発生器１２ａが、図１７に示すように、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、信号Ｔｘ（１）～Ｔｘ（Ｎ）をＱ回繰り返し生成する。即ち、信号発生器１２ａは、それぞれの送信サイクルｃにおいて、ヒットｑ（ｑ＝１，・・・，Ｑ）の信号Ｔｘ（１）～Ｔｘ（Ｎ）を生成する。Ｑは、２以上の整数である。即ち、信号発生器１２ａは、それぞれの送信サイクルｃにおいて、ヒット数Ｑの信号Ｔｘ（１）～Ｔｘ（Ｎ）を生成する。
図１７は、信号発生器１２ａにより生成されるヒット数Ｑのアップチャープの信号Ｔｘ（１）～Ｔｘ（Ｎ）を示す説明図である。図１７の例では、Ｎ＝３である。ヒット数Ｑは、例えば１６である。
このため、アンテナ１１－ｎ（ｎ＝１，・・・，Ｎ）は、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、Ｑ個の送信信号Ｔｘ’（ｎ）に係る送信波のそれぞれを対象物体が存在している空間に放射する。

アンテナ１１－１～１１－Ｎのそれぞれは、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、Ｎ×Ｑ個の反射波を受信する。
アンテナ１１－ｎから送信波が放射された空間に、対象物体として、被測定者ｋ（ｋ＝１，・・・，Ｋ）のほかに、移動体が存在していれば、移動体による反射波についてもアンテナ１１－１～１１－Ｎに入射される。
被測定者ｋは、例えば、ベッドに寝ている概ね静止状態の人間、又は、椅子等に座っている概ね静止状態の人間である。移動体は、送信波が放射された空間内を移動している人間等である。
アンテナ１１－ｎ（ｎ＝１，・・・，Ｎ）は、それぞれの送信サイクルｃにおいて、対象物体による反射波を受信する。即ち、アンテナ１１－ｎ（ｎ＝１，・・・，Ｎ）は、それぞれの送信サイクルｃにおいて、被測定者ｋ及び移動体のそれぞれによる反射波を受信する。
アンテナ１１－ｎ（ｎ＝１，・・・，Ｎ）は、それぞれの送信サイクルｃにおいて、被測定者ｋ及び移動体のそれぞれによる反射波の信号を含むＮ×Ｑ個の受信信号のそれぞれをサーキュレータ１３－ｎに出力する。
サーキュレータ１３－ｎは、それぞれの送信サイクルｃにおいて、アンテナ１１－ｎから出力されたＮ×Ｑ個の受信信号を信号受信部１４－ｎに出力する。

信号受信部１４－ｎ（ｎ＝１，・・・，Ｎ）は、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、サーキュレータ１３－ｎから出力されたＮ×Ｑ個の受信信号のそれぞれに対する受信処理を実施する。
信号受信部１４－ｎは、それぞれの送信サイクルｃにおいて、Ｎ×Ｑ個の受信データＳ（ｔ，ｇ，ｈ，ｃ）のそれぞれをバイタル測定装置２０に出力する。
これにより、送信サイクルｃ＝１～Ｃにおいて、信号受信部１４－１～１４－Ｎから、全部でＮ×Ｑ×Ｃ個の受信データＳ（ｔ，ｇ，ｈ，ｃ）がバイタル測定装置２０に与えられる。

バイタル測定装置２０の信号取得部２１は、信号受信部１４－１～１４－Ｎから、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、Ｎ×Ｑ個の受信データＳ（ｔ，ｇ，ｈ，ｃ）を取得する。
信号取得部２１は、それぞれの送信サイクルｃにおいて、Ｎ×Ｑ個の受信データＳ（ｔ，ｇ，ｈ，ｃ）のそれぞれを信号抑圧部２９に出力する。

信号抑圧部２９は、信号取得部２１から、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、Ｎ×Ｑ個の受信データＳ（ｔ，ｇ，ｈ，ｃ）を取得する。
信号抑圧部２９は、それぞれの送信サイクルｃにおいて、Ｎ×Ｑ個の受信データＳ（ｔ，ｇ，ｈ，ｃ）のそれぞれに含まれている、移動体による反射波に係る信号を抑圧する。
信号抑圧部２９は、それぞれの送信サイクルｃにおいて、信号抑圧後の信号として、Ｎ×Ｑ個の受信データＳ_Ｙ（ｔ，ｇ，ｈ，ｃ）をフーリエ変換部２２に出力する。
以下、信号抑圧部２９による信号抑圧処理を具体的に説明する。

図１８Ａは、アンテナ１１－１～１１－ＮからＴＤＭ（Ｔｉｍｅ Ｄｉｖｉｓｉｏｎ Ｍｕｌｔｉｐｌｅｘｉｎｇ）方式で送信波が放射されたのち、信号取得部２１により取得される受信データＳ（ｔ，ｇ，ｈ，ｃ）を示している。
図１８Ａにおいて、横軸は、時間を示し、縦軸は、振幅を示している。ＭＩＭＯ（Ｍｕｌｔｉｐｌｅ－Ｉｎｐｕｔ Ｍｕｌｔｉｐｌｅ－Ｏｕｔｐｕｔ）による仮想チャネル数は、Ｇ×Ｈである。
受信データＳ（ｔ，ｇ，ｈ，ｃ）には、図１８Ａに示すように、受信系ハードウェアに起因するＤＣ（Ｄｉｒｅｃｔ Ｃｕｒｒｅｎｔ）オフセット成分が重畳されている。
信号取得部２１から信号抑圧部２９に与えられる受信データＳ（ｔ，ｇ，ｈ，ｃ）の数は、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、Ｎ×Ｑ個である。

信号抑圧部２９は、図１８Ｂに示すように、仮想チャネル毎に、受信データＳ（ｔ，ｇ，ｈ，ｃ）に重畳されているＤＣオフセット成分を除去するＡＤ（Ａｎａｌｏｇ－ｔｏ－Ｄｉｇｉｔａｌ）オフセット補正を行う。ＡＤオフセット補正自体は、公知の技術であるため詳細な説明を省略する。
図１８Ｂは、ＡＤオフセット補正後の受信データＳ’（ｔ，ｇ，ｈ，ｃ）を示している。
図１８Ｂにおいて、横軸は、時間を示し、縦軸は、振幅を示している。

信号抑圧部２９は、図１８Ｃに示すように、それぞれの送信サイクルｃにおいて、仮想チャネル毎に、Ｇ×Ｈ個のＡＤオフセット補正後の受信データＳ’（ｔ，ｇ，ｈ，ｃ）を平均化する。
信号抑圧部２９によって、受信データＳ’（ｔ，ｇ，ｈ，ｃ）が平均化されることで、受信データＳ’（ｔ，ｇ，ｈ，ｃ）に含まれている、移動体による反射波に係る信号が抑圧される。
信号抑圧部２９は、それぞれの送信サイクルｃにおいて、Ｎ個の受信データＳ（ｔ，ｇ，ｈ，ｃ）として、Ｎ個の信号抑圧後の受信データＳ_Ｙ（ｔ，ｇ，ｈ，ｃ）をフーリエ変換部２２に出力する。
図１８Ｃは、移動体による反射波に係る信号抑圧後の受信データＳ_Ｙ（ｔ，ｇ，ｈ，ｃ）を示している。
図１８Ｃにおいて、横軸は、時間を示し、縦軸は、振幅を示している。
フーリエ変換部２２以降の処理は、図１に示すバイタル測定装置２０と同様であるため、説明を省略する。

以上の実施の形態２では、信号取得部２１により取得された受信信号に含まれている、移動体による反射波に係る信号を抑圧する信号抑圧部２９を備え、フーリエ変換部２２が、信号抑圧部２９による信号抑圧後の受信信号を時間方向にフーリエ変換するように、図１５に示すバイタル測定装置２０を構成した。したがって、図１５に示すバイタル測定装置２０は、図１に示すバイタル測定装置２０と同様に、被測定者が複数の場合であっても、それぞれの被測定者のバイタルを検出することができるほか、空間内に移動体が存在していても、被測定者のバイタルを検出することができる。

実施の形態３．
実施の形態３では、信号取得部２１により取得された受信信号に基づいて、移動体を検出する移動体検出部３０を備えるバイタル測定装置２０について説明する。

図１９は、実施の形態３に係るバイタル測定装置２０を含むバイタル測定システムを示す構成図である。図１９において、図１５と同一符号は同一又は相当部分を示すので説明を省略する。
図２０は、実施の形態３に係るバイタル測定装置２０のハードウェアを示すハードウェア構成図である。図２０において、図１６と同一符号は同一又は相当部分を示すので説明を省略する。

図１９に示すバイタル測定装置２０は、信号取得部２１、信号抑圧部２９、フーリエ変換部２２、マップ算出部２３、バイタル推定部２４及び移動体検出部３０を備えている。
移動体検出部３０は、例えば、図２０に示す移動体検出回路３６によって実現される。
移動体検出部３０は、信号取得部２１により取得された受信データＳ（ｔ，ｇ，ｈ，ｃ）に含まれている、移動体による反射波に係る信号を抽出する。
移動体検出部３０は、移動体による反射波に係る信号に基づいて、移動体を検出する。

図１９では、バイタル測定装置２０の構成要素である信号取得部２１、信号抑圧部２９、フーリエ変換部２２、マップ算出部２３、バイタル推定部２４及び移動体検出部３０のそれぞれが、図２０に示すような専用のハードウェアによって実現されるものを想定している。即ち、バイタル測定装置２０が、信号取得回路３１、信号抑圧回路３５、フーリエ変換回路３２、マップ算出回路３３、バイタル推定回路３４及び移動体検出回路３６によって実現されるものを想定している。
信号取得回路３１、信号抑圧回路３５、フーリエ変換回路３２、マップ算出回路３３、バイタル推定回路３４及び移動体検出回路３６のそれぞれは、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ＡＳＩＣ、ＦＰＧＡ、又は、これらを組み合わせたものが該当する。

バイタル測定装置２０の構成要素は、専用のハードウェアによって実現されるものに限るものではなく、バイタル測定装置２０が、ソフトウェア、ファームウェア、又は、ソフトウェアとファームウェアとの組み合わせによって実現されるものであってもよい。
バイタル測定装置２０が、ソフトウェア又はファームウェア等によって実現される場合、信号取得部２１、信号抑圧部２９、フーリエ変換部２２、マップ算出部２３、バイタル推定部２４及び移動体検出部３０におけるそれぞれの処理手順をコンピュータに実行させるためのプログラムが図３に示すメモリ４１に格納される。そして、図３に示すプロセッサ４２がメモリ４１に格納されているプログラムを実行する。

また、図２０では、バイタル測定装置２０の構成要素のそれぞれが専用のハードウェアによって実現される例を示し、図３では、バイタル測定装置２０がソフトウェア又はファームウェア等によって実現される例を示している。しかし、これは一例に過ぎず、バイタル測定装置２０における一部の構成要素が専用のハードウェアによって実現され、残りの構成要素がソフトウェア又はファームウェア等によって実現されるものであってもよい。

次に、図１９に示すバイタル測定システムの動作について説明する。移動体検出部３０以外は、図１５に示すバイタル測定システムと同様である。このため、ここでは、移動体検出部３０の動作のみを説明する。

移動体検出部３０は、信号取得部２１から、それぞれの送信サイクルｃ（ｃ＝１，・・・，Ｃ）において、Ｎ×Ｃ個の受信データＳ（ｔ，ｇ，ｈ，ｃ）を取得する。
移動体検出部３０は、Ｎ×Ｃ個の受信データＳ（ｔ，ｇ，ｈ，ｃ）に基づいて、空間に存在している移動体の追尾処理を実施することで、時刻の経過に伴って変化する、移動体の存在位置を推定する。移動体の追尾処理自体は、公知の技術であるため詳細な説明を省略する。
例えば、空間に存在している移動体の数が１つであれば、１つの移動体の存在位置を推定する。例えば、２つの移動体として、第１の移動体と第２の移動体とが空間に存在していれば、第１の移動体の存在位置と、第２の移動体の存在位置とを推定する。
移動体検出部３０は、移動体の存在位置の推定結果に基づいて、移動体が存在している２次元方位を検出する。２次元方位は、バイタル測定装置２０から移動体を見たアジマス方向と、バイタル測定装置２０から移動体を見たエレベーション方向とである。

ここでは、移動体検出部３０が、空間に存在している移動体の追尾処理を実施することで、移動体の存在位置を推定している。しかし、これは一例に過ぎず、移動体検出部３０が、Ｎ×Ｃ個の受信データＳ（ｔ，ｇ，ｈ，ｃ）をヒット方向にフーリエ変換することで、移動体のドップラー周波数ｆｄを算出するようにしてもよい。移動体のドップラー周波数ｆｄを算出することで、移動体の速度を算出することができる。

以上の実施の形態３では、信号取得部２１により取得された受信信号に基づいて、移動体を検出する移動体検出部３０を備えるように、図１９に示すバイタル測定装置２０を構成した。したがって、図１９に示すバイタル測定装置２０は、図１５に示すバイタル測定装置２０と同様に、被測定者が複数の場合であっても、それぞれの被測定者のバイタルを検出することができるほか、移動体を検出することができる。

なお、本開示は、各実施の形態の自由な組み合わせ、あるいは各実施の形態の任意の構成要素の変形、もしくは各実施の形態において任意の構成要素の省略が可能である。

本開示は、バイタル測定装置、バイタル測定方法及びバイタル測定システムに適している。

１０ センサ、１１－１～１１－Ｎ アンテナ、１２ 信号送信部、１２ａ 信号発生器、１２ｂ 出力先選択部、１３－１～１３－Ｎ サーキュレータ、１４－１～１４－Ｎ 信号受信部、２０ バイタル測定装置、２１ 信号取得部、２２ フーリエ変換部、２３ マップ算出部、２４ バイタル推定部、２５ 存在位置特定部、２６ バイタル推定処理部、２７ 呼吸数推定部、２８ 心拍数推定部、２９ 信号抑圧部、３０ 移動体検出部、３１ 信号取得回路、３２ フーリエ変換回路、３３ マップ算出回路、３４ バイタル推定回路、３５ 信号抑圧回路、３６ 移動体検出回路、４１ メモリ、４２ プロセッサ。

Claims (10)

1. 対象物体による反射波を受信するアンテナから、前記反射波の受信信号を取得する信号取得部と、
   前記信号取得部により取得された受信信号を時間方向にフーリエ変換するフーリエ変換部と、
   前記フーリエ変換部によるフーリエ変換後の信号を用いて、バイタル測定装置からの複数の距離ビンのそれぞれに対応する複素電力の２次元方位マップを算出するマップ算出部と、
   前記マップ算出部により算出されたそれぞれの距離ビンに対応する複素電力の２次元方位マップの時間変化に基づいて、前記対象物体に含まれる被測定者が存在している位置を特定し、前記被測定者が存在している位置についての複素電力の位相の時間変化から、前記被測定者のバイタルを推定するバイタル推定部と
   を備えたバイタル測定装置。
2. 前記対象物体の中に、１人以上の被測定者が含まれており、
   前記バイタル推定部は、
   前記マップ算出部により算出されたそれぞれの距離ビンに対応する複素電力の２次元方位マップの時間変化に基づいて、それぞれの被測定者が存在している位置を特定し、それぞれの被測定者が存在している位置についての複素電力の位相の時間変化を示す位相変化信号を出力する存在位置特定部と、
   前記存在位置特定部から出力されたそれぞれの位相変化信号に基づいて、それぞれの被測定者のバイタルを推定するバイタル推定処理部とを備えていることを特徴とする請求項１記載のバイタル測定装置。
3. 前記存在位置特定部は、
   前記マップ算出部により算出された複数の２次元方位マップのそれぞれに含まれている複数の複素電力の中から、距離ビンが互いに同一であり、かつ、２次元方位が互いに同一である複素電力の組を複数抽出し、それぞれの組に含まれている複数の複素電力についての位相の時間変化を描くフィッティング円を算出し、被測定者が存在している位置として、算出した複数のフィッティング円の中で、第１の閾値以上の半径を有するフィッティング円に係る組に含まれている複素電力の位置を特定することを特徴とする請求項２記載のバイタル測定装置。
4. 前記存在位置特定部は、
   前記マップ算出部により算出された複数の２次元方位マップのそれぞれに含まれている複数の複素電力の中から、距離ビンが互いに同一であり、かつ、２次元方位が互いに同一である複素電力の組を複数抽出し、それぞれの組に含まれている複数の複素電力についての位相の時間変化を描くフィッティング円を算出し、被測定者が存在している位置として、算出した複数のフィッティング円の中で、フィッティング円の誤差であるフィッティング誤差が第２の閾値以下であるフィッティング円に係る組に含まれている複素電力の位置を特定することを特徴とする請求項２記載のバイタル測定装置。
5. 前記存在位置特定部は、
   前記マップ算出部により算出された複数の２次元方位マップのそれぞれに含まれている複数の複素電力の中から、距離ビンが互いに同一であり、かつ、２次元方位が互いに同一である複素電力の組を複数抽出し、それぞれの組に含まれている複数の複素電力についての位相の時間変化を描くフィッティング円を算出し、被測定者が存在している位置として、算出した複数のフィッティング円の中で、第１の閾値以上の半径を有し、かつ、フィッティング円の誤差であるフィッティング誤差が第２の閾値以下であるフィッティング円に係る組に含まれている複素電力の位置を特定することを特徴とする請求項２記載のバイタル測定装置。
6. 前記バイタル推定処理部は、
   前記存在位置特定部から出力されたそれぞれの位相変化信号に基づいて、それぞれの被測定者の呼吸数を推定する呼吸数推定部と、
   前記存在位置特定部から出力されたそれぞれの位相変化信号に基づいて、それぞれの被測定者の心拍数を推定する心拍数推定部とを備えていることを特徴とする請求項２記載のバイタル測定装置。
7. 前記信号取得部により取得された受信信号に含まれている、移動体による反射波に係る信号を抑圧する信号抑圧部を備え、
   前記フーリエ変換部は、
   前記信号抑圧部による信号抑圧後の受信信号を時間方向にフーリエ変換することを特徴とする請求項１記載のバイタル測定装置。
8. 前記信号取得部により取得された受信信号に基づいて、前記移動体を検出する移動体検出部を備えたことを特徴とする請求項７記載のバイタル測定装置。
9. 信号取得部が、対象物体による反射波を受信するアンテナから、前記反射波の受信信号を取得し、
   フーリエ変換部が、前記信号取得部により取得された受信信号を時間方向にフーリエ変換し、
   マップ算出部が、前記フーリエ変換部によるフーリエ変換後の信号を用いて、バイタル測定装置からの複数の距離ビンのそれぞれに対応する複素電力の２次元方位マップを算出し、
   バイタル推定部が、前記マップ算出部により算出されたそれぞれの距離ビンに対応する複素電力の２次元方位マップの時間変化に基づいて、前記対象物体に含まれる被測定者が存在している位置を特定し、前記被測定者が存在している位置についての複素電力の位相の時間変化から、前記被測定者のバイタルを推定する
   バイタル測定方法。
10. アンテナと、
    前記アンテナから、対象物体が存在している空間に送信波を放射させる信号送信部と、
    前記対象物体による反射後の送信波である反射波を受信する前記アンテナから、前記反射波の受信信号を取得する信号取得部と、
    前記信号取得部により取得された受信信号を時間方向にフーリエ変換するフーリエ変換部と、
    前記フーリエ変換部によるフーリエ変換後の信号を用いて、バイタル測定装置からの複数の距離ビンのそれぞれに対応する複素電力の２次元方位マップを算出するマップ算出部と、
    前記マップ算出部により算出されたそれぞれの距離ビンに対応する複素電力の２次元方位マップの時間変化に基づいて、前記対象物体に含まれる被測定者が存在している位置を特定し、前記被測定者が存在している位置についての複素電力の位相の時間変化から、前記被測定者のバイタルを推定するバイタル推定部と
    を備えたバイタル測定システム。

Images