ホーム > 特許ランキング > 公立大学法人北九州市立大学
知財求人 - 知財ポータルサイト「IP Force」
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2023037552
(43)【公開日】2023-03-15
(54)【発明の名称】生体データ取得装置、及び生体データ取得方法
(51)【国際特許分類】
A61B 5/11 20060101AFI20230308BHJP
【ＦＩ】
A61B5/11 110
【審査請求】未請求
【請求項の数】7
【出願形態】ＯＬ
(21)【出願番号】P 2022011676
(22)【出願日】2022-01-28
(31)【優先権主張番号】17/466,549
(32)【優先日】2021-09-03
(33)【優先権主張国・地域又は機関】US
(71)【出願人】
【識別番号】506060258
【氏名又は名称】公立大学法人北九州市立大学
(74)【代理人】
【識別番号】100088155
【弁理士】
【氏名又は名称】長谷川 芳樹
(74)【代理人】
【識別番号】100145012
【弁理士】
【氏名又は名称】石坂 泰紀
(74)【代理人】
【識別番号】100212026
【弁理士】
【氏名又は名称】中村 真生
(72)【発明者】
【氏名】梶原 昭博
【テーマコード（参考）】
4C038
【Ｆターム（参考）】
4C038VA04
4C038VB33
(57)【要約】
【課題】生体データを非接触で精度良く検知する。
【解決手段】生体データ取得装置は、第１～第Ｎの距離変動データから、対応強度が時系列に沿って並べられた時間変動データを生成する第１の生成部と、時間変動データに基づいて、生体の検知部位における生体データを生成するように構成された第２の生成部と、抑圧処理を行うことで第ｎの距離変動データを生成する信号強度算出部とを備える。第１～第Ｊの受信信号データのうちの第ｊの受信信号データは、検知領域を距離方向に分割して得られる第１～第Ｋのレンジビンにおける信号強度ａｊ１～ａｊＫを含む。抑圧処理は、第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉを周波数変換して得られた結果に基づいて係数ｖｉを算出することと、第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと係数ｖｉとを乗算することで、第ｎの距離変動データを算出することと、を含む。
【選択図】図２２
【解決手段】生体データ取得装置は、第１～第Ｎの距離変動データから、対応強度が時系列に沿って並べられた時間変動データを生成する第１の生成部と、時間変動データに基づいて、生体の検知部位における生体データを生成するように構成された第２の生成部と、抑圧処理を行うことで第ｎの距離変動データを生成する信号強度算出部とを備える。第１～第Ｊの受信信号データのうちの第ｊの受信信号データは、検知領域を距離方向に分割して得られる第１～第Ｋのレンジビンにおける信号強度ａｊ１～ａｊＫを含む。抑圧処理は、第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉを周波数変換して得られた結果に基づいて係数ｖｉを算出することと、第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと係数ｖｉとを乗算することで、第ｎの距離変動データを算出することと、を含む。
【選択図】図２２
【特許請求の範囲】
【請求項1】
第１～第Ｎ（Ｎは２以上の自然数）の距離変動データが時系列順に並べられた時系列データを記憶するように構成された記憶部であって、前記第１～第Ｎの距離変動データはそれぞれ、互いに異なる時刻において、広帯域又は超広帯域の無線電波が生体において反射された反射波に基づいて得られる、距離に対する信号強度の変化を示すデータである、記憶部と、
前記第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち前記生体の所定の検知部位に対応する対応強度を取得する強度取得処理を前記第１～第Ｎの距離変動データに対して行うことにより、前記対応強度が時系列に沿って並べられた時間変動データを生成するように構成された第１の生成部と、
前記時間変動データに基づいて、前記生体の検知部位における生体データを生成するように構成された第２の生成部と、
第１～第Ｊ（Ｊは２以上の自然数）の受信信号データに対して抑圧処理を行うことで、前記第ｎの距離変動データを生成するように構成された信号強度算出部とを備え、
前記第１～第Ｊの受信信号データは、フレーム周期内の時系列順に並ぶ異なる複数の時刻において、前記生体が存在する検知領域にて前記無線電波が反射された反射波に基づいて得られ、
前記第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、前記第ｊの受信信号データは、前記検知領域を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含み、
前記抑圧処理は、
前記第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは、１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、
前記第ｊの受信信号データについて、前記第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと前記係数ｖｉとを乗算することで、前記第ｎの距離変動データを算出することと、を含む、生体データ取得装置。
前記第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち前記生体の所定の検知部位に対応する対応強度を取得する強度取得処理を前記第１～第Ｎの距離変動データに対して行うことにより、前記対応強度が時系列に沿って並べられた時間変動データを生成するように構成された第１の生成部と、
前記時間変動データに基づいて、前記生体の検知部位における生体データを生成するように構成された第２の生成部と、
第１～第Ｊ（Ｊは２以上の自然数）の受信信号データに対して抑圧処理を行うことで、前記第ｎの距離変動データを生成するように構成された信号強度算出部とを備え、
前記第１～第Ｊの受信信号データは、フレーム周期内の時系列順に並ぶ異なる複数の時刻において、前記生体が存在する検知領域にて前記無線電波が反射された反射波に基づいて得られ、
前記第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、前記第ｊの受信信号データは、前記検知領域を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含み、
前記抑圧処理は、
前記第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは、１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、
前記第ｊの受信信号データについて、前記第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと前記係数ｖｉとを乗算することで、前記第ｎの距離変動データを算出することと、を含む、生体データ取得装置。
【請求項2】
第１～第Ｎ（Ｎは２以上の自然数）の距離変動データが時系列順に並べられた時系列データを記憶するように構成された記憶部であって、前記第１～第Ｎの距離変動データはそれぞれ、互いに異なる時刻において、広帯域又は超広帯域の無線電波が生体において反射された反射波に基づいて得られる、距離に対する信号強度の変化を示すデータである、記憶部と、
前記第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち前記生体の所定の検知部位に対応する対応強度を取得する強度取得処理を前記第１～第Ｎの距離変動データに対して行うことにより、前記対応強度が時系列に沿って並べられた時間変動データを生成するように構成された第１の生成部と、
前記時間変動データに基づいて、前記生体の検知部位における生体データを生成するように構成された第２の生成部と、
第１～第Ｊ（Ｊは２以上の自然数）の受信信号データに対して抑圧処理を行うことで、前記第ｎの距離変動データを生成するように構成された信号強度算出部とを備え、
前記第１～第Ｊの受信信号データは、フレーム周期内の時系列順に並ぶ異なる複数の時刻において、前記生体が存在する検知領域にて前記無線電波が反射された反射波に基づいて得られ、
前記第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、前記第ｊの受信信号データは、前記検知領域を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含み、
前記抑圧処理は、
前記第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、
前記第１～第Ｊの受信信号データのうちの複数の受信信号データを積分して積分信号データを生成することであって、前記積分信号データは、前記第１～第Ｋのレンジビンにおける信号強度ｂ１～ｂＫを含むことと、
前記積分信号データについて前記第１～第Ｋのレンジビンのそれぞれにおいて信号強度ｂｉと前記係数ｖｉとを乗算して、前記第ｎの距離変動データを算出することとを含む、生体データ取得装置。
前記第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち前記生体の所定の検知部位に対応する対応強度を取得する強度取得処理を前記第１～第Ｎの距離変動データに対して行うことにより、前記対応強度が時系列に沿って並べられた時間変動データを生成するように構成された第１の生成部と、
前記時間変動データに基づいて、前記生体の検知部位における生体データを生成するように構成された第２の生成部と、
第１～第Ｊ（Ｊは２以上の自然数）の受信信号データに対して抑圧処理を行うことで、前記第ｎの距離変動データを生成するように構成された信号強度算出部とを備え、
前記第１～第Ｊの受信信号データは、フレーム周期内の時系列順に並ぶ異なる複数の時刻において、前記生体が存在する検知領域にて前記無線電波が反射された反射波に基づいて得られ、
前記第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、前記第ｊの受信信号データは、前記検知領域を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含み、
前記抑圧処理は、
前記第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、
前記第１～第Ｊの受信信号データのうちの複数の受信信号データを積分して積分信号データを生成することであって、前記積分信号データは、前記第１～第Ｋのレンジビンにおける信号強度ｂ１～ｂＫを含むことと、
前記積分信号データについて前記第１～第Ｋのレンジビンのそれぞれにおいて信号強度ｂｉと前記係数ｖｉとを乗算して、前記第ｎの距離変動データを算出することとを含む、生体データ取得装置。
【請求項3】
前記信号強度算出部は、前記抑圧処理を実行することの前にフレーム設定処理を実行し、
前記フレーム設定処理は、前記生体の動作の速度を示す動作情報を取得することと、前記動作情報が示す速度に基づいて前記フレーム周期の長さを設定することとを含む、請求項１又は２に記載の装置。
前記フレーム設定処理は、前記生体の動作の速度を示す動作情報を取得することと、前記動作情報が示す速度に基づいて前記フレーム周期の長さを設定することとを含む、請求項１又は２に記載の装置。
【請求項4】
第１～第Ｎ（Ｎは２以上の自然数）の距離変動データが時系列順に並べられた時系列データを記憶するように構成された記憶部であって、前記第１～第Ｎの距離変動データはそれぞれ、互いに異なる時刻において、広帯域又は超広帯域の無線電波が生体において反射された反射波に基づいて得られる、距離に対する信号強度の変化を示すデータである、記憶部と、
前記第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち前記生体の所定の検知部位に対応する対応強度を取得する強度取得処理を前記第１～第Ｎの距離変動データに対して行うことにより、前記対応強度が時系列に沿って並べられた時間変動データを生成するように構成された第１の生成部と、
前記時間変動データに基づいて、前記生体の検知部位における生体データを生成するように構成された第２の生成部と、
第１のフレーム周期と、当該第１のフレーム周期とは異なる第２のフレーム周期とのそれぞれにおいて、第１～第Ｊ（Ｊは２以上の自然数）の受信信号データを取得するように構成された信号強度算出部とを備え、
前記第１～第Ｊの受信信号データは、１つのフレーム周期内の時系列順に並ぶ異なる複数の時刻において、前記生体が存在する検知領域にて前記無線電波が反射された反射波に基づいて得られ、
前記第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、前記第ｊの受信信号データは、前記検知領域を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含み、
前記信号強度算出部は、所定の抑圧処理を行うことで、前記第１のフレーム周期における前記第ｎの距離変動データを生成するように構成されており、
前記抑圧処理は、
前記第１のフレーム周期及び前記第２のフレーム周期のそれぞれにおいて、前記第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは、１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、
前記第１～第Ｋのレンジビンの全てについて、前記第１のフレーム周期と前記第２のフレーム周期との間での前記係数ｖｉの変化の程度が小さくなるにつれて値が小さくなるように、前記係数ｖｉの変化の程度に基づいて係数ｃｉを算出することと、
前記第１のフレーム周期の前記第ｊの受信信号データについて、前記第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと前記係数ｃｉとを乗算することで、前記第ｎの距離変動データを算出することとを含む、生体データ取得装置。
前記第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち前記生体の所定の検知部位に対応する対応強度を取得する強度取得処理を前記第１～第Ｎの距離変動データに対して行うことにより、前記対応強度が時系列に沿って並べられた時間変動データを生成するように構成された第１の生成部と、
前記時間変動データに基づいて、前記生体の検知部位における生体データを生成するように構成された第２の生成部と、
第１のフレーム周期と、当該第１のフレーム周期とは異なる第２のフレーム周期とのそれぞれにおいて、第１～第Ｊ（Ｊは２以上の自然数）の受信信号データを取得するように構成された信号強度算出部とを備え、
前記第１～第Ｊの受信信号データは、１つのフレーム周期内の時系列順に並ぶ異なる複数の時刻において、前記生体が存在する検知領域にて前記無線電波が反射された反射波に基づいて得られ、
前記第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、前記第ｊの受信信号データは、前記検知領域を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含み、
前記信号強度算出部は、所定の抑圧処理を行うことで、前記第１のフレーム周期における前記第ｎの距離変動データを生成するように構成されており、
前記抑圧処理は、
前記第１のフレーム周期及び前記第２のフレーム周期のそれぞれにおいて、前記第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは、１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、
前記第１～第Ｋのレンジビンの全てについて、前記第１のフレーム周期と前記第２のフレーム周期との間での前記係数ｖｉの変化の程度が小さくなるにつれて値が小さくなるように、前記係数ｖｉの変化の程度に基づいて係数ｃｉを算出することと、
前記第１のフレーム周期の前記第ｊの受信信号データについて、前記第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと前記係数ｃｉとを乗算することで、前記第ｎの距離変動データを算出することとを含む、生体データ取得装置。
【請求項5】
第１～第Ｎ（Ｎは２以上の自然数）の距離変動データが時系列順に並べられた時系列データを生成することであって、前記第１～第Ｎの距離変動データはそれぞれ、互いに異なる時刻において、広帯域又は超広帯域の無線電波が生体において反射された反射波に基づいて得られる、距離に対する信号強度の変化を示すデータである、前記時系列データを生成することと、
前記第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち前記生体の所定の検知部位に対応する対応強度を取得する強度取得処理を前記第１～第Ｎの距離変動データに対して行うことにより、前記対応強度が時系列に沿って並べられた時間変動データを生成することと、
前記時間変動データに基づいて、前記生体の検知部位における生体データを生成することと、
第１～第Ｊ（Ｊは２以上の自然数）の受信信号データに対して抑圧処理を行うことで、前記第ｎの距離変動データを生成することとを含み、
前記第１～第Ｊの受信信号データは、フレーム周期内の時系列順に並ぶ異なる複数の時刻において、前記生体が存在する検知領域にて前記無線電波が反射された反射波に基づいて得られ、
前記第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、前記第ｊの受信信号データは、前記検知領域を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含み、
前記抑圧処理は、
前記第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、
前記第ｊの受信信号データについて、前記第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと前記係数ｖｉとを乗算することで、前記第ｎの距離変動データを算出することと、を含む、生体データ取得方法。
前記第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち前記生体の所定の検知部位に対応する対応強度を取得する強度取得処理を前記第１～第Ｎの距離変動データに対して行うことにより、前記対応強度が時系列に沿って並べられた時間変動データを生成することと、
前記時間変動データに基づいて、前記生体の検知部位における生体データを生成することと、
第１～第Ｊ（Ｊは２以上の自然数）の受信信号データに対して抑圧処理を行うことで、前記第ｎの距離変動データを生成することとを含み、
前記第１～第Ｊの受信信号データは、フレーム周期内の時系列順に並ぶ異なる複数の時刻において、前記生体が存在する検知領域にて前記無線電波が反射された反射波に基づいて得られ、
前記第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、前記第ｊの受信信号データは、前記検知領域を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含み、
前記抑圧処理は、
前記第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、
前記第ｊの受信信号データについて、前記第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと前記係数ｖｉとを乗算することで、前記第ｎの距離変動データを算出することと、を含む、生体データ取得方法。
【請求項6】
第１～第Ｎ（Ｎは２以上の自然数）の距離変動データが時系列順に並べられた時系列データを生成することであって、前記第１～第Ｎの距離変動データはそれぞれ、互いに異なる時刻において、広帯域又は超広帯域の無線電波が生体において反射された反射波に基づいて得られる、距離に対する信号強度の変化を示すデータである、前記時系列データを生成することと、
前記第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち前記生体の所定の検知部位に対応する対応強度を取得する強度取得処理を前記第１～第Ｎの距離変動データに対して行うことにより、前記対応強度が時系列に沿って並べられた時間変動データを生成することと、
前記時間変動データに基づいて、前記生体の検知部位における生体データを生成することと、
第１～第Ｊ（Ｊは２以上の自然数）の受信信号データに対して抑圧処理を行うことで、前記第ｎの距離変動データを生成することとを含み、
前記第１～第Ｊの受信信号データは、フレーム周期内の時系列順に並ぶ異なる複数の時刻において、前記生体が存在する検知領域にて前記無線電波が反射された反射波に基づいて得られ、
前記第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、前記第ｊの受信信号データは、前記検知領域を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含み、
前記抑圧処理は、
前記第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、
前記第１～第Ｊの受信信号データのうちの複数の受信信号データを積分して積分信号データを生成することであって、前記積分信号データは、前記第１～第Ｋのレンジビンにおける信号強度ｂ１～ｂＫを含むことと、
前記積分信号データについて前記第１～第Ｋのレンジビンのそれぞれにおいて信号強度ｂｉと前記係数ｖｉとを乗算して、前記第ｎの距離変動データを算出することとを含む、生体データ取得方法。
前記第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち前記生体の所定の検知部位に対応する対応強度を取得する強度取得処理を前記第１～第Ｎの距離変動データに対して行うことにより、前記対応強度が時系列に沿って並べられた時間変動データを生成することと、
前記時間変動データに基づいて、前記生体の検知部位における生体データを生成することと、
第１～第Ｊ（Ｊは２以上の自然数）の受信信号データに対して抑圧処理を行うことで、前記第ｎの距離変動データを生成することとを含み、
前記第１～第Ｊの受信信号データは、フレーム周期内の時系列順に並ぶ異なる複数の時刻において、前記生体が存在する検知領域にて前記無線電波が反射された反射波に基づいて得られ、
前記第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、前記第ｊの受信信号データは、前記検知領域を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含み、
前記抑圧処理は、
前記第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、
前記第１～第Ｊの受信信号データのうちの複数の受信信号データを積分して積分信号データを生成することであって、前記積分信号データは、前記第１～第Ｋのレンジビンにおける信号強度ｂ１～ｂＫを含むことと、
前記積分信号データについて前記第１～第Ｋのレンジビンのそれぞれにおいて信号強度ｂｉと前記係数ｖｉとを乗算して、前記第ｎの距離変動データを算出することとを含む、生体データ取得方法。
【請求項7】
第１～第Ｎ（Ｎは２以上の自然数）の距離変動データが時系列順に並べられた時系列データを生成することであって、前記第１～第Ｎの距離変動データはそれぞれ、互いに異なる時刻において、広帯域又は超広帯域の無線電波が生体において反射された反射波に基づいて得られる、距離に対する信号強度の変化を示すデータである、前記時系列データを生成することと、
前記第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち前記生体の所定の検知部位に対応する対応強度を取得する強度取得処理を前記第１～第Ｎの距離変動データに対して行うことにより、前記対応強度が時系列に沿って並べられた時間変動データを生成することと、
前記時間変動データに基づいて、前記生体の検知部位における生体データを生成することと、
第１のフレーム周期と、当該第１のフレーム周期とは異なる第２のフレーム周期とのそれぞれにおいて第１～第Ｊ（Ｊは２以上の自然数）の受信信号データを取得することと、 所定の抑圧処理を行うことで、前記第１のフレーム周期における前記第ｎの距離変動データを生成することとを含み、
前記第１～第Ｊの受信信号データは、１つのフレーム周期内の時系列順に並ぶ異なる複数の時刻において、前記生体が存在する検知領域にて前記無線電波が反射された反射波に基づいて得られ、
前記第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、前記第ｊの受信信号データは、前記検知領域を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含み、
前記抑圧処理は、
前記第１のフレーム周期及び前記第２のフレーム周期のそれぞれにおいて、前記第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは、１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、
前記第１～第Ｋのレンジビンの全てについて、前記第１のフレーム周期と前記第２のフレーム周期との間での前記係数ｖｉの変化の程度が小さくなるにつれて値が小さくなるように、前記係数ｖｉの変化の程度に基づいて係数ｃｉを算出することと、
前記第１のフレーム周期の前記第ｊの受信信号データについて、前記第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと前記係数ｃｉとを乗算することで、前記第ｎの距離変動データを算出することとを含む、生体データ取得方法。
前記第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち前記生体の所定の検知部位に対応する対応強度を取得する強度取得処理を前記第１～第Ｎの距離変動データに対して行うことにより、前記対応強度が時系列に沿って並べられた時間変動データを生成することと、
前記時間変動データに基づいて、前記生体の検知部位における生体データを生成することと、
第１のフレーム周期と、当該第１のフレーム周期とは異なる第２のフレーム周期とのそれぞれにおいて第１～第Ｊ（Ｊは２以上の自然数）の受信信号データを取得することと、 所定の抑圧処理を行うことで、前記第１のフレーム周期における前記第ｎの距離変動データを生成することとを含み、
前記第１～第Ｊの受信信号データは、１つのフレーム周期内の時系列順に並ぶ異なる複数の時刻において、前記生体が存在する検知領域にて前記無線電波が反射された反射波に基づいて得られ、
前記第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、前記第ｊの受信信号データは、前記検知領域を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含み、
前記抑圧処理は、
前記第１のフレーム周期及び前記第２のフレーム周期のそれぞれにおいて、前記第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは、１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、
前記第１～第Ｋのレンジビンの全てについて、前記第１のフレーム周期と前記第２のフレーム周期との間での前記係数ｖｉの変化の程度が小さくなるにつれて値が小さくなるように、前記係数ｖｉの変化の程度に基づいて係数ｃｉを算出することと、
前記第１のフレーム周期の前記第ｊの受信信号データについて、前記第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと前記係数ｃｉとを乗算することで、前記第ｎの距離変動データを算出することとを含む、生体データ取得方法。
【発明の詳細な説明】
【技術分野】
【0001】
本開示は、生体データ取得装置、及び生体データ取得方法に関する。
【背景技術】
【0002】
特許文献１は、ユーザの心拍に関する心拍データを測定し、ユーザの心拍データと識別情報とを送信対象端末に送信するウェアラブルデバイスを開示している。このウェアラブルデバイスは、リストバンド型又は腕時計型の測定機器であり、ウェアラブルデバイスを装着したユーザの心拍データが測定される。
【先行技術文献】
【特許文献】
【0003】
【特許文献1】特開２０１７－１６７８０２号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
特許文献１に記載の測定機器を用いてユーザの心拍データを測定する場合、測定機器をユーザに接触させた状態で測定を行う必要がある。このため、測定に際して測定機器の着脱に伴う手間が発生し、測定機器の装着によりユーザに不快感を与えてしまう懸念がある。そのため、電波センサ等を用いて心拍データ等の生体データを非接触で測定するための技術が検討されている。ところが、電波センサ等によって非接触で生体データを測定しようとする場合、生体データの測定対象である生体（例えば、人体、動物など）が動いていると、生体からの反射波に基づいて生体データを取得する際に影響が生じうる。
【0005】
そこで、本開示は、生体データを非接触で精度良く検知することが可能な生体データ取得装置、及び生体データ取得方法を説明する。
【課題を解決するための手段】
【0006】
本開示の１つの観点に係る生体データ取得装置は、第１～第Ｎ（Ｎは２以上の自然数）の距離変動データが時系列順に並べられた時系列データを記憶するように構成された記憶部であって、第１～第Ｎの距離変動データはそれぞれ、互いに異なる時刻において、広帯域又は超広帯域の無線電波が生体において反射された反射波に基づいて得られる、距離に対する信号強度の変化を示すデータである、記憶部と、第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち生体の所定の検知部位に対応する対応強度を取得する強度取得処理を第１～第Ｎの距離変動データに対して行うことにより、対応強度が時系列に沿って並べられた時間変動データを生成するように構成された第１の生成部と、時間変動データに基づいて、生体の検知部位における生体データを生成するように構成された第２の生成部と、第１～第Ｊ（Ｊは２以上の自然数）の受信信号データに対して抑圧処理を行うことで、第ｎの距離変動データを生成するように構成された信号強度算出部とを備える。第１～第Ｊの受信信号データは、フレーム周期内の時系列順に並ぶ異なる複数の時刻において、生体が存在する検知領域にて上記無線電波が反射された反射波に基づいて得られる。第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、第ｊの受信信号データは、検知領域を距離方向に分割して得られる第１～第Ｋのレンジビンにおける信号強度ａｊ１～ａｊＫを含む。抑圧処理は、第１～第Ｋ（Ｋは２以上の自然数）のレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、第ｊの受信信号データについて、第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと係数ｖｉとを乗算することで、第ｎの距離変動データを算出することと、を含む。。
【発明の効果】
【0007】
本開示に係る生体データ取得装置、及び生体データ取得方法によれば、生体データを非接触で精度良く検知することが可能となる。
【図面の簡単な説明】
【0008】
【発明を実施するための形態】
【0009】
以下に、本開示に係る実施形態の一例について、図面を参照しつつより詳細に説明する。以下の説明において、同一要素又は同一機能を有する要素には同一符号を用いることとし、重複する説明は省略する。
【0010】
［生体データ取得システム］
生体データ取得システム１の構成について、図１～図３を参照して説明する。生体データ取得システム１は、送受信装置１０と、生体データ取得装置２０とを備える。生体データ取得システム１は、例えば室内又は屋外の生体Ｈの検知部位における呼吸や拍動などの各生体データを取得するように構成されている。室内は、例えば、リビング、オフィス、寝室、トイレ、浴室、乗物などの内部を含む。
生体データ取得システム１の構成について、図１～図３を参照して説明する。生体データ取得システム１は、送受信装置１０と、生体データ取得装置２０とを備える。生体データ取得システム１は、例えば室内又は屋外の生体Ｈの検知部位における呼吸や拍動などの各生体データを取得するように構成されている。室内は、例えば、リビング、オフィス、寝室、トイレ、浴室、乗物などの内部を含む。
【0011】
検知対象となる生体Ｈとしては、例えば人間、動物などが挙げられるが、本明細書では生体Ｈが人間であるものとして説明する。生体Ｈは、動作していてもよく、椅子等に着座した状態で作業していてもよく、静止していてもよい。なお、生体Ｈは静止していても、無意識下においてわずかに体が動いていることがある。本明細書では、生体Ｈによる意図的な動作、無意識下でのわずかな動作（微体動）、及び外部要因に由来して生体Ｈに生ずる動作（振動等）を含めて生体Ｈの「動き」として説明する。生体Ｈによる意図的な動作は、生体の全体の動作、及び生体の一部分の動作を含む。生体の全体の動作は、例えば、歩行等を含む。生体の一部分の動作は、例えば、立位、座位、及び臥位状態での生体の一部分の動作を含む。検知部位は、例えば生体Ｈの心臓部分を含む部位（胸部又は背中）であってもよい。
【0012】
送受信装置１０は、生体データを検知する対象領域として予め設定された検知領域Ｒに向けて、広帯域又は超広帯域（ＵＷＢ：Ultra Wide Band）の無線電波を発信するように構成されている。検知領域Ｒは、例えば生体Ｈの検知部位に無線電波が到達するように設定されている。送受信装置１０は、例えば室内の壁に設けられていてもよい。送受信装置１０は、図２に示されるように、発信機１１（発信部）と、受信機１２（受信部）と、通信機１３（通信部）と、制御部１４と、バス１５とを含む。
【0013】
発信機１１は、制御部１４の指示に基づいて、検知領域Ｒに向けて広帯域又は超広帯域の無線電波を発信するように構成されている。発信機１１のアンテナ水平面における指向角は、例えば、４０°～８０°程度であってもよいし、６０°～８０°程度であってもよい。発信機１１のアンテナ垂直面における指向角は、例えば、３０°～４０°程度であってもよいし、３０°以下であってもよい。これらの指向角の範囲に対応して検知領域Ｒが設定される。
【0014】
本明細書において「広帯域」とは、周波数帯域幅が１００ＭＨｚ以上で且つ５００ＭＨｚ以下の場合をいうものとする。そのため、発信機１１が発信する広帯域の無線電波の周波数帯域幅は、例えば、１００ＭＨｚ以上であってもよいし、３００ＭＨｚ以上であってもよい。「超広帯域」とは、周波数帯域幅が５００ＭＨｚを超える場合をいうものとする。そのため、発信機１１が発信する超広帯域の無線電波の周波数帯域幅は、例えば、３ＧＨｚ以上であってもよいが、電波法及びコストパフォーマンスに鑑みて４ＧＨｚ以下であってもよい。広帯域又は超広帯域の無線電波を用いる場合、無線電波の発信出力の電力スペクトルを小さくすることができ、生体Ｈに対する無線電波の影響を小さくすることができる。以下では、特に説明がない限り、「無線電波」は広帯域又は超広帯域の無線電波を意味する。
【0015】
受信機１２は、発信機１１から発信された無線電波の反射波を受信可能に構成されている。発信機１１において広帯域又は超広帯域の無線電波を用いているので、受信機１２は、無線電波の反射波を反射パス長に対応した時間軸上で分離して受信することができる。すなわち、受信機１２は、検知領域Ｒの距離方向において複数に区切られたレンジビン（周波数帯域幅で決まる）ごとに、信号を抽出することが可能である。各レンジビンの幅は、例えば１０ｃｍ以下であってもよいし、３ｃｍ～４ｃｍ程度であってもよい。検知領域Ｒにおけるレンジビンの数は、レンジビンの幅及び検知領域Ｒの大きさに応じて適宜変更可能である。
【0016】
受信機１２によるデータ取得の時間間隔（サンプリング周期）は、適宜設定しうる。生体データ取得装置２０による拍動成分の抽出（詳しくは後述する）を容易にするために、受信機１２は、生体の拍動に対応する周波数を基準として、反射波をオーバサンプリングしてもよい。受信機１２は、例えば、拍動に応じた周波数範囲の上限値（例えば１．６Ｈｚ）の２倍より大きいサンプリング周波数（例えば上限値の整数倍）にて反射波をオーバサンプリングしてもよい。なお、受信機１２以外でオーバサンプリングが行われてもよい。
【0017】
通信機１３は、生体データ取得装置２０の通信機２３と通信可能に構成されている。通信機１３の通信機２３との通信方式は特に限定されず、例えば、無線通信であってもよい。無線通信の例として、ＬＴＥ（Long Term Evolution）、ＬＴＥ－Ａ（LTE-Advanced）、ＳＵＰＥＲ３Ｇ、ＩＭＴ－Ａｄｖａｎｃｅｄ、４Ｇ、５Ｇ、ＦＲＡ（Future Radio Access）、Ｗ－ＣＤＭＡ（登録商標）、ＧＳＭ（登録商標）、ＣＤＭＡ２０００、ＵＭＢ（Ultra Mobile Broadband）、ＩＥＥＥ ８０２．１１（Ｗｉ－Ｆｉ）、ＩＥＥＥ ８０２．１６（ＷｉＭＡＸ）、ＩＥＥＥ ８０２．２０、ＵＷＢ、Ｂｌｕｅｔｏｏｔｈ（登録商標）、その他の通信方式が用いられてもよい。
【0018】
制御部１４は、バス１５を介して、発信機１１、受信機１２及び通信機１３との間で信号の送受信を行い、これらの動作を制御するように構成されている。制御部１４は、例えば、プロセッサ１６と、メモリ１７と、ストレージ１８とを含む。
【0019】
プロセッサ１６は、メモリ１７、ストレージ１８等のハードウェアに所定のソフトウェア（プログラム）を読み込むと所定の演算を行い、発信機１１からの無線電波の発信、受信機１２が受信した反射波の受信、通信機１３による通信、メモリ１７及びストレージ１８におけるデータの読み出し又は書き込みを実行する。これにより、送受信装置１０における各機能が実現される。
【0020】
生体データ取得装置２０は、送受信装置１０により受信された反射波（受信データ）に基づいて、生体Ｈの生体データを取得するように構成されている。例えば、生体Ｈが呼吸をすると、それに応じて胸部が膨張収縮する。そのため、受信データには生体Ｈの呼吸に関する情報が含まれる。また、生体Ｈの心臓の拍動に応じて心臓が膨張収縮する。そのため、受信データには生体Ｈの拍動に関する情報も含まれる。
【0021】
生体データ取得装置２０は、呼吸に関する情報及び拍動に関する情報の少なくとも一方を生体データとして取得する。生体データ取得装置２０は、図１に示されるように、送受信装置１０と離間して配置されていてもよい。生体データ取得装置２０は、例えば、生体Ｈが存在する室内とは異なる場所に設けられてもよい。生体データ取得装置２０に含まれる各要素を収容する筐体は、送受信装置１０に含まれる各要素を収容する筐体と別体であってもよい。
【0022】
生体データ取得装置２０は、図３に示されるように、通信機２３と、制御部２４と、バス２５とを備える。通信機２３は、上述したように、送受信装置１０の通信機１３と通信可能に構成されている。生体データ取得装置２０は、通信機２３が内蔵されたコンピュータ装置又はパーソナルコンピュータにより構成されていてもよい。
【0023】
制御部２４は、バス２５を介して、通信機２３との間で信号の送受信を行い、これの動作を制御するように構成されている。制御部２４は、例えば、プロセッサ２６と、メモリ２７（記憶部）と、ストレージ２８と（記憶部）とを含む。
【0024】
プロセッサ２６は、メモリ２７、ストレージ２８等のハードウェアに所定のソフトウェア（プログラム）を読み込むと所定の演算を行い、送受信装置１０を介して取得した反射波の解析、通信機２３による通信、メモリ２７及びストレージ２８におけるデータの読み出し又は書き込みを実行する。これにより、生体データ取得装置２０における各機能が実現される。
【0025】
続いて、送受信装置１０及び生体データ取得装置２０における各機能について、図４～図９を参照して説明する。送受信装置１０の制御部１４は、図４に示されるように、機能ブロックとして、送受信処理部３１と、記憶部３２とを含む。
【0026】
送受信処理部３１は、発信機１１、受信機１２及び通信機１３の間で信号を送受信する機能を有する。具体的には、送受信処理部３１は、発信機１１に指示信号を送信して、発信機１１から無線電波を発信させる。送受信処理部３１は、受信機１２から反射波のデータ（反射波に応じた受信データ）を受信し、当該データを通信機２３に出力する。送受信処理部３１は、通信機１３に指示信号を送信して、受信データを生体データ取得装置２０に送信させる。
【0027】
記憶部３２は、種々のデータを記憶する機能を有する。記憶部３２が記憶するデータとしては、例えば、読み出したプログラム、発信機１１の動作設定データ、受信機１２が受信した反射波に関する受信データ等が挙げられる。
【0028】
生体データ取得装置２０の制御部２４は、図４に示されるように、機能ブロックとして、データ取得部４１と、信号強度算出部４２と、記憶部４３と、データ補正部４４（補正部）と、時間変動データ生成部４５（第１の生成部）と、生体データ生成部４６（第２の生成部）とを含む。
【0029】
データ取得部４１は、通信機１３及び通信機２３を介して、反射波に応じた受信データを取得する機能を有する。データ取得部４１は、取得した受信データを信号強度算出部４２に出力する。
【0030】
信号強度算出部４２は、データ取得部４１から出力された受信データに基づいて、各レンジビンのそれぞれにおける信号強度を算出する機能を有する。換言すると、信号強度算出部４２は、所定の時刻における、距離に対する信号強度の変化（レンジビン毎の信号強度の変化）を示す距離変動データを算出する。信号強度算出部４２は、受信データに対して、生体Ｈが検知領域Ｒに存在しない場合の無線電波の反射に基づいた参照受信データを用いて差分処理を施した上で、距離変動データを算出してもよい。この差分処理は、例えば、人または動いている物体（動体）を検知したり、それらの距離や動線を検出したり、マルチパスによる影響を低減するために行われる。
【0031】
信号強度算出部４２によって算出された信号強度は、記憶部４３に出力される。以下では、サンプリング周期ごとに、時刻ｔ１，ｔ２，・・・，ｔＮ（Ｎは２以上の自然数）における距離変動データをそれぞれ、「第１の距離変動データ」、「第２の距離変動データ」、・・・、「第Ｎの距離変動データ」などと称し、第１～第Ｎの距離変動データのうち時刻ｔｎ（ｎは１～Ｎの自然数）における距離変動データを「第ｎの距離変動データ」と称する。信号強度算出部４２は、無線電波が生体Ｈから反射された反射波に基づいて、互いに異なる時刻ｔ１～ｔＮにおける第１～第Ｎの距離変動データをそれぞれ時間順に算出してもよい。
【0032】
記憶部４３は、種々のデータを記憶する機能を有する。記憶部４３が記憶するデータとしては、例えば、読み出したプログラム、信号強度算出部４２が算出した複数の距離変動データ、データ補正部４４が補正した時系列データ（以下、「補正時系列データ」という。）、時間変動データ生成部４５が生成した時間変動データ、生体データ生成部４６が生成した生体データ（例えば、心拍間隔データ）、距離変動データの蓄積単位を示す設定値（以下、「設定時間」という。）等が挙げられる。設定時間は、メモリの容量に応じて予め定められてもよい。設定時間は、例えば、２０秒～３分程度の範囲内で定められてもよいし、３分以上であってもよい。
【0033】
記憶部４３は、設定時間において信号強度算出部４２により算出された第１～第Ｎの距離変動データ（所定時間分の第１～第Ｎの距離変動データ）が時系列順に並べられた時系列データを記憶する。時系列データは、距離（レンジ）軸、時間軸、及び信号強度軸を含む３次元グラフにより示される。第ｎの距離変動データは、距離成分、時間成分、及び信号強度成分を含んでいる。図５（ａ）は、時系列データの一例を示す。図５（ａ）に示される時系列データは、６０秒間において得られた距離変動データである。この時系列データは、被験者（生体Ｈ）が送受信装置１０に正対した測定環境において、送受信装置１０から約１．３ｍだけ離れた位置に静止した状態で測定を開始し、約３０秒後から送受信装置１０に向かって徐々に近づいていった場合のデータである。そのため、被験者の胸部及び心臓（検知部位）からの反射波の信号強度が、距離変動データにおけるピーク付近に現れており、当該ピークでの距離成分が徐々に小さくなっている。
【0034】
データ補正部４４は、記憶部４３に記憶されている時系列データを補正する機能を有する。データ補正部４４は、例えば、第ｎの距離変動データにおける検知部位に対応する距離成分が基準値に近づくように第ｎの距離変動データの距離成分を補正する補正処理を行ってもよい。基準値は、例えば、生体データの取得ごとに、任意に設定される値であってもよい。データ補正部４４は、第１の距離変動データに含まれるピークでの距離成分の位置を基準値に設定してもよい。データ補正部４４は、当該補正処理を第１～第Ｎの距離変動データに対して行うことにより、補正された時系列データ（補正時系列データ）を得てもよい。データ補正部４４は、補正処理として、例えば、第１～第Ｎの距離変動データを分割した複数のブロックそれぞれにおいて、生体Ｈの動きを近似（一次補正）し、検知部位に対応する距離成分を補正（二次補正）してもよい。そして、データ補正部４４は、複数のブロック間での距離成分を補正（三次補正）してもよい。データ補正部４４は、補正時系列データを時間変動データ生成部４５に出力する。
【0035】
ここで、補正処理の例について具体的に説明する。この補正処理は、生体Ｈが動いていた場合であっても、検知部位からの反射に基づいた信号（信号強度の時間変動）を抽出するために行われる。まず、データ補正部４４は、第１～第Ｎの距離変動データにおいてピーク（散乱点）をそれぞれ算出する。このときに算出する「ピーク」は、信号強度が所定値よりも大きく、且つ極大となる点であってもよい。
【0036】
ところで、マルチパス等の外乱の影響により、１つの距離変動データに複数のピークが算出される場合がある。この場合には、検知部位からの反射に基づくピーク以外のピークが時系列データには含まれうる。ここで、図６（ａ）に、外乱等に基づくピークを含む複数のピークｐｎが距離軸－時間軸の平面に白丸印（○印）でプロットされたグラフの一例を示す。以下では、図６（ａ）に例示される複数のピークｐｎを用いて、外乱等に由来するピークを含む距離変動データを補正するための補正処理を説明する。
【0037】
まず、データ補正部４４は、時系列データが時間軸において所定の時間間隔（例えば数秒～数十秒間隔）で分割された複数のブロックごとに、複数のピークｐｎに対して近似線Ｃを描く（生成処理）。以下では、時系列データを時系列に沿って順に分割して得られる複数のブロックを、「第１のブロックＢ１」、「第２のブロックＢ２」、・・・、「第ＭのブロックＢＭ」（Ｍは、２以上の自然数）とそれぞれ称する場合がある。第１～第ＭのブロックＢ１～ＢＭのうち１つのブロックを「第ｍのブロックＢｍ」（ｍは、１～Ｍの自然数）と称する場合がある。
【0038】
近似線Ｃは、１つのブロック内（所定の時間間隔）での生体Ｈの移動を推定した動線軌跡を意味する。近似線は、近似直線、あるいは二次又は多項近似を用いた多項近似曲線であってもよい。人の動作は、１秒以内であれば、ほぼ直線の動線軌跡となり、数秒程度であっても動線軌跡を二次曲線で近似することが可能である。
【0039】
図６（ａ）の例では、１つの距離変動データごとに近似線Ｃを中心として左右に２つずつ目盛が示されている。１つの目盛間隔は１つのレンジビンに対応している。すなわち、図６（ａ）の目盛は、近似線Ｃを中心として距離軸の正側及び負側にそれぞれ２つずつレンジビンを描いたものである。データ補正部４４は、これら４つのレンジビンから外れた位置にあるピークを以降の処理において除外してもよい。レンジビンの設定数は、特に限定されない。例えば、近似線Ｃを中心として距離軸の正側及び負側にそれぞれ１つずつであってもよいし、ｘ個（ｘは３以上の自然数）ずつであってもよい。あるいは、近似線Ｃを中心として距離軸の正側に設定されるレンジビンの数は、近似線Ｃを中心として距離軸の負側に設定されるレンジビンの数よりも少なくてもよいし、多くてもよいし、同一であってもよい。
【0040】
一般的に、近似線Ｃ付近のピークと胸部付近からの反射波とは一致することが多いが、生体Ｈが送受信装置１０に対して斜め方向を向いていると、近似線Ｃ付近のピークと胸部付近からの反射波とが異なることもある。そのため、データ補正部４４は、近似線Ｃを基準とした所定の範囲（以下、「散乱範囲」という。）において、検知部位からの反射点を検出してもよい。例えば、生体Ｈが人間であり、その横幅が４０～５０ｃｍ程度と仮定し、１ＧＨｚの無線電波が用いられる場合、図６（ａ）に例示されるように、散乱範囲は近似線Ｃから±２レンジビンに設定されてもよい。なお、この場合に、レンジビンの幅は１５ｃｍ程度であってもよい。
【0041】
次に、データ補正部４４は、図６（ｂ）に示されるように、近似線Ｃが第１のブロックＢ１での基準値（ここでは、ｘ＝０の位置）と一致するように、ブロックＢ１内に含まれる距離変動データの距離成分を距離軸に沿って平行移動する（一次補正：ブロック補正処理）。ブロックＢ１での基準値は、任意に設定される値である。次に、データ補正部４４は、散乱範囲において、レンジビンごとに信号強度の時間変動を算出して、当該時間変動が最大値を示すレンジビンを選択してもよい。データ補正部４４は、例えば、同じレンジビン同士で、第ｎの距離変動データと第ｎ＋１の距離変動データとの間で信号強度の差分を求めることにより、レンジビンごとに信号強度の時間変動を算出してもよい。
【0042】
図７（ａ）は、レンジビンにおける信号強度の時間変動が最も大きい点を検知部位からの反射点ｐｔとして検出した結果の一例を示す。図７（ａ）において、反射点ｐｔは黒丸印（●印）で示されている。次に、データ補正部４４は、図７（ｂ）に示されるように、検出された反射点ｐｔの距離成分が基準値と一致するように、ブロックＢ１内に含まれる距離変動データの距離成分を距離軸に沿って平行移動する（二次補正）。外乱等を含みうるピークｐｎではなく、レンジビンにおける信号強度の時間変動が最も大きい点を反射点ｐｔとして検出しているので、図７（ａ）又は図７（ｂ）に示されるように、時刻によって、ピークｐｎが反射点ｐｔと一致する場合もあり、ピークｐｎが反射点ｐｔと一致しない場合もある。以上の生成処理、一次補正及び二次補正が行われることで、マルチパス等の外乱に由来するピークが除去されうる。
【0043】
データ補正部４４は、他のブロックについても上記と同様に、ブロックごとに同様の生成処理、一次補正及び二次補正を行う。以上のように、データ補正部４４は、第ｍのブロックＢｍにおいて、第１～第Ｎの距離変動データのうち第ｍのブロックＢｍに含まれる複数の距離変動データにおける信号強度のピークｐｎでの距離成分に基づいて近似線Ｃを生成する生成処理と、ブロックＢｍについて定められる基準値に当該近似線Ｃが近づくように複数の距離変動データの距離成分を補正する処理（一次補正及び二次補正）とを行う。データ補正部４４は、生成処理、一次補正及び二次補正を、第１～第Ｍのブロックそれぞれに対して行う。
【0044】
ブロックごとの生成処理、一次補正及び二次補正では、ブロックごとに任意に定められる基準値が用いられてもよい。例えば、ブロック内において最初の距離変動データにおけるピークの位置が、当該ブロックでの基準値に設定されてもよい。このため、複数のブロックそれぞれにて生成処理、一次補正及び二次補正を行った後では、複数のブロック間において基準値の位置が互いにずれている。そのため、データ補正部４４は、複数のブロック間の基準値の位置が一致するように、それぞれのブロック内に含まれる距離変動データの距離成分を距離軸に沿って平行移動する（三次補正）。データ補正部４４は、三次補正において、最初のブロックＢ１における基準値を全てのブロックの基準値として採用してもよい。
【0045】
以上のように、生成処理、一次補正、二次補正及び三次補正を経ることにより、図５（ｂ）に示されるように、生体Ｈが基準値で静止しているような補正時系列データが得られる。なお、図６（ａ）等では、説明を簡単にするために基準値の位置をｘ＝０としているが、図５（ａ）及び図５（ｂ）に示される例では、最初のブロックでの基準値の位置が１．３ｍ付近であり、１．３ｍ付近にピーク及び検知部位からの反射点が揃えられている。
【0046】
時間変動データ生成部４５は、第ｎの距離変動データにおける信号強度のうち生体Ｈの検知部位に対応する対応強度を取得する強度取得処理を行う機能を有する。時間変動データ生成部４５は、例えば、所定時間分の第１～第Ｎの距離変動データに対して強度取得処理を行うことにより、対応強度が時系列に沿って並べられた時間変動データを生成する。時間変動データ生成部４５は、例えば、補正時系列データ内の基準値（一例として、第１のブロックＢ１についての基準値）において時系列に沿って並べられた信号強度の集合を、時間変動データとして得てもよい。時間変動データは、検知部位からの反射波に応じた信号強度の時間変動を示している。図８（ａ）に示される波形「ｓ」は、時間変動データ生成部４５により生成された時間変動データの一例を示す。データ補正部４４による補正処理を経ることで、時間変動データでは、生体Ｈの動き（動作又は微体動）による影響が低減されている。時間変動データには、生体Ｈの検知部位における生体情報（呼吸及び拍動）が含まれている。
【0047】
生体データ生成部４６は、生体Ｈの検知部位における生体データを生成する機能を有する。生体データ生成部４６は、例えば、時間変動データに対して多重解像度解析を行うことにより、生体Ｈの検知部位における生体データを生成してもよい。生体データ生成部４６は、時間変動データに対する多重解像度解析において、当該時間変動データを高周波成分（高解像度成分）と低周波成分（低解像度成分）とに分解するウェーブレット変換を繰り返す。多重解像度解析における解像度レベルは、例えば、５～１０程度であってもよいし、７程度であってもよい。生体データ生成部４６は、例えば、基底関数となるウェーブレット関数として、Haarウェーブレット、Daubechiesウェーブレット、Symletウェーブレット、又はCoifletウェーブレットを用いてウェーブレット変換を行ってもよい。
【0048】
ここで、時間変動データに対する多重解像度解析の例について具体的に説明する。図８（ｂ）～図８（ｉ）は、図８（ａ）に示される時間変動データを、多重解像度解析により分解した複数（ここでは８個）の波形（以下、「分解波形」という。）の一例を示す。生体データ生成部４６は、まず、時間変動データに対して１回目のウェーブレット変換を行うことで、高周波成分と低周波成分とに分解する。この時に分解された高周波成分は、図８（ｉ）に示される波形「ｄ１」である。生体データ生成部４６は、１回目のウェーブレット変換において分解された低周波成分に対して、２回目のウェーブレットを行う。この時に分解された高周波成分は、図８（ｈ）に示される波形「ｄ２」である。
【0049】
生体データ生成部４６は、３回目以降のウェーブレット変換を同様に繰り返す。この例では、時間変動データに対して、７回のウェーブレット変換を行うことにより、当該時間変動データが８個のレベル（周波数帯）に分解されている。なお、図８（ｇ）、図８（ｆ）、図８（ｅ）、図８（ｄ）及び図８（ｃ）に示される波形「ｄ３」～「ｄ７」はそれぞれ、３～７回目のウェーブレット変換により分解された高周波成分である。図８（ｂ）に示される波形「ａ７」は、７回目のウェーブレット変換により分解された低周波成分である。
【0050】
生体データ生成部４６は、複数の分解波形から拍動信号及び呼吸信号の少なくとも一方を抽出する。拍動信号は、時間変動データに含まれる生体Ｈの拍動に応じた成分である。呼吸信号は、時間変動データに含まれる生体Ｈの呼吸に応じた成分である。生体データ生成部４６は、例えば、複数の分解波形のうち、呼吸に応じた周波数範囲（例えば、０．５Ｈｚ以下）を有する分解波形を、呼吸信号として抽出してもよい。生体データ生成部４６は、複数の分解波形のうち、拍動に応じた周波数範囲（例えば０．８Ｈｚ～１．６Ｈｚ）を有する分解波形を、拍動信号として抽出してもよい。
【0051】
図８（ｂ）～図８（ｉ）に示される例では、生体データ生成部４６は、図８（ｂ）に示される波形「ａ７」を呼吸信号として抽出してもよい。生体データ生成部４６は、図８（ｃ）に示される波形「ｄ７」を拍動信号として抽出してもよい。生体データ生成部４６は、複数の分解波形（例えば、図８（ｃ）及び図８（ｄ）に示される２つの波形）の周波数範囲が、拍動（呼吸）に応じた周波数範囲に対応する場合には、当該複数の分解波形を再合成することで、再合成された波形を拍動信号（呼吸信号）として抽出してもよい。再合成する際には、生体データ生成部４６は、当該複数の分解波形同士を加算することで、再合成された波形を生成してもよい。なお、生体データ生成部４６は、多重解像度解析に代えて独立成分分析法または移動平均差分法を用いて、時間変動データから呼吸信号と拍動信号とを分離してもよい。
【0052】
生体データ生成部４６は、複数の分解波形から拍動信号を抽出する際に、拍動とは異なる生体Ｈの動きに対応する信号（以下、「非拍動信号」という。）を除外してもよい。拍動とは異なる生体Ｈの動きは、例えば、生体Ｈ自身に由来する拍動以外の周期的な動作、及び外部要因に由来して生体Ｈに生ずる周期的な動作を含む。前者の一例は、生体Ｈの呼吸動作を含む。後者の一例は、生体Ｈが搭乗している乗物の運転に伴い乗物から生体Ｈに作用する振動を含む。生体データ生成部４６は、拍動信号を抽出する際に、複数の分解波形のうちの非拍動信号を除外してもよい。生体データ生成部４６は、非拍動信号が除外された後の残りの複数の分解波形を再合成することで得られる信号を、拍動信号として抽出してもよい。生体データ生成部４６は、例えば、生体Ｈの呼吸動作に対応する図８（ｂ）に示される波形「ａ７」を除外したうえで、波形「ａ７」が除外された後の残りの波形「ｄ１」～「ｄ７」を再合成することで得られる信号を、拍動信号として抽出してもよい。
【0053】
生体データ生成部４６は、拍動信号に含まれる極大となる複数の点において互いに隣り合う点同士の時間間隔を算出していくことで、図９に示されるように、心拍間隔の時間変動（Ｒ－Ｒ間隔：ＲＲＩ）を示す心拍間隔データを生体データとして生成してもよい。生体データ生成部４６は、拍動信号から心拍間隔データを生成する際に、心拍間隔を高精度に算出するために、拍動信号に対してスペクトル解析を行ってもよい。生体データ生成部４６は、例えば、自己回帰モデル（ＡＲ：autoregressive model）を用いたスペクトル解析を拍動信号に対して行ってもよい。ＡＲモデルを用いたスペクトル解析では、生体データ生成部４６は、バーグ法によりＡＲモデルのパラメータ（係数）を求めてもよく、赤池情報量基準（ＡＩＣ：Akaike's Information Criterion）によりＡＲモデルの次数を決定してもよい。
【0054】
生体データ生成部４６は、ＡＲモデルを用いたスペクトル解析によりノイズ成分が除去された拍動信号から、心拍間隔データを算出してもよい。生体データ生成部４６は、このように算出された心拍間隔データに基づいて、単位時間あたりの心拍の数（心拍数：heart rate，平均心拍：average heart rate）を生体データとして求めてよい。生体データ生成部４６は、上記と同様に、単位時間（例えば１分間）あたりの呼吸の数（呼吸数：respiratory rate）を生体データとして求めてもよい。なお、生体データ生成部４６は、ＡＲモデルに代えて、最大エントロピー法（ＭＥＭ：Maximum Entropy Method）又はメムカルク（MemCalc）解析システムを用いてスペクトル解析を行ってもよい。
【0055】
生体データ取得装置２０は、生体データ生成部４６により生成された生体データ（例えば、心拍間隔データ）を外部に出力する機能を有していてもよい。生体データ取得システム１は、例えば、報知機（不図示）を備えていてもよく、生体データ取得装置２０は、当該報知機に生体データを出力してもよい。生体データ取得装置２０は、送受信装置１０及び生体データ取得装置２０とは別の場所（建物）に設けられている報知機に通信機２３を介してデータを送信してもよい。生体データ取得装置２０は、取得した生体データを記憶部４３に記憶してもよい。
【0056】
［生体データ取得方法］
続いて、図１０を参照して、生体データ取得方法について説明する。なお、ここでは、生体データとして生体Ｈの拍動に関する情報（心拍間隔の時間変動）を取得する場合を例として説明する。
続いて、図１０を参照して、生体データ取得方法について説明する。なお、ここでは、生体データとして生体Ｈの拍動に関する情報（心拍間隔の時間変動）を取得する場合を例として説明する。
【0057】
まず、送受信装置１０（制御部１４）の送受信処理部３１は、発信機１１に指示して、発信機１１から無線電波を検知領域Ｒに向けて発信させる。次に、送受信処理部３１は、受信機１２を介して無線電波の反射波に基づいた受信データを受信して、当該受信データを通信機１３に出力する。そして、送受信処理部３１は、通信機１３に指示して、受信データを生体データ取得装置２０の通信機２３に送信させる。
【0058】
次に、生体データ取得装置２０（制御部２４）のデータ取得部４１は、通信機２３に送信された受信データを、時刻ｔ１における受信データとして取得する。データ取得部４１は、取得した時刻ｔ１の受信データを信号強度算出部４２に出力する。
【0059】
次に、信号強度算出部４２は、時刻ｔ１の受信データに基づいて、距離に対する信号強度の変動を示す第１の距離変動データを算出する（図１０のステップＳ１）。信号強度算出部４２は、受信データに対して、生体Ｈが検知領域Ｒに存在しない場合の無線電波の反射に基づいた参照受信データを用いて差分処理を施した上で、第１の距離変動データを算出してもよい。信号強度算出部４２は、算出した第１の距離変動データを記憶部４３に出力する。
【0060】
次に、記憶部４３は、第１の距離変動データを記憶する（図１０のステップＳ２）。そして、制御部２４は、初期時刻から設定時間が経過したかどうかを判断する（図１０のステップＳ３）。設定時間が経過していないと判断された場合（ステップＳ３：ＮＯ）、制御部２４は、ステップＳ１～Ｓ３を繰り返す。これにより、設定時間において、第１～第Ｎの距離変動データが記憶部４３に記憶される。
【0061】
設定時間が経過したと判断された場合（ステップＳ３：ＹＥＳ）、データ補正部４４は、記憶部４３に記憶されている時系列データ（第１～第Ｎの距離変動データ）を補正する（図１０のステップＳ４）。データ補正部４４は、例えば、第ｎの距離変動データにおける検知部位に対応する距離成分（検知部位からの反射点）が基準値に近づくように第ｎの距離変動データの距離成分を補正する補正処理を行う。データ補正部４４は、当該補正処理を第１～第Ｎの距離変動データに対して行うことにより、補正された時系列データ（補正時系列データ）を生成する（図５（ａ）及び図５（ｂ）参照）。データ補正部４４は、図６（ａ）～図７（ｂ）に示される例のように、一次補正、二次補正、及び三次補正を行うことで、補正時系列データを生成してもよい。データ補正部４４は、補正時系列データを時間変動データ生成部４５に出力する。
【0062】
次に、時間変動データ生成部４５は、補正時系列データに基づいて、時間変動データを生成する（図１０のステップＳ５）。時間変動データ生成部４５は、例えば、第ｎの距離変動データにおける信号強度のうち生体Ｈの検知部位に対応する対応強度を取得する強度取得処理を行う。時間変動データ生成部４５は、第１～第Ｎの距離変動データに対して強度取得処理を行うことにより、対応強度が時系列に沿って並べられた時間変動データを生成する。時間変動データ生成部４５は、補正時系列データにおいて、生体データの取得処理ごとに設定される基準値に沿って並べられた信号強度の集合を、時間変動データとして得てもよい（図５（ｂ）及び図８（ａ）参照）。時間変動データ生成部４５は、生成した時間変動データを生体データ生成部４６に出力する。
【0063】
次に、生体データ生成部４６は、時間変動データに対して多重解像度解析を行う（図１０のステップＳ６）。生体データ生成部４６は、時間変動データに対して、高周波成分と低周波成分とに分解するウェーブレット変換を繰り返すことで、互いに異なる周波数範囲を有する複数の分解波形（図８（ｂ）～図８（ｉ）参照）を生成する。
【0064】
次に、生体データ生成部４６は、複数の分解波形から拍動成分を抽出する（図１０のステップＳ７）。生体データ生成部４６は、例えば、複数の分解波形のうち、拍動に応じた周波数範囲（例えば０．８Ｈｚ～１．６Ｈｚ）を有する分解波形を、拍動成分を示す拍動信号として抽出してもよい。生体データ生成部４６は、複数の分解波形の周波数範囲が、拍動に応じた周波数範囲に対応する場合には、当該複数の分解波形を再合成することで、再合成された波形を拍動信号として抽出してもよい。
【0065】
次に、生体データ生成部４６は、拍動信号に基づいて心拍間隔データを生成する（ステップＳ８）。生体データ生成部４６は、例えば、拍動信号に含まれる極大となる複数の点において互いに隣り合う点同士の時間間隔を算出していくことで、心拍間隔データを生体データとして生成してもよい（図９参照）。生体データ生成部４６は、拍動信号から心拍間隔データを生成する際に、心拍間隔を高精度に算出するために、拍動信号に対してスペクトル解析を行ってもよい。生体データ生成部４６は、例えば、ＡＲモデルを用いたスペクトル解析によりノイズ成分が除去された拍動信号から、心拍間隔データを算出してもよい。
【0066】
生体データ取得装置２０は、生体データ生成部４６により生成された心拍間隔データを外部（例えば外部の報知機）に出力してもよい。生体データ取得装置２０は、生体データ生成部４６により生成された心拍間隔データを記憶部４３に記憶してもよい。生体データ取得装置２０（制御部２４）は、ステップＳ３において最初の設定時間終了後に、ステップＳ１～Ｓ３を繰り返すことで、設定時間ごとに時系列データ（第１～第Ｎの距離変動データ）を繰り返し取得してもよい。生体データ取得装置２０は、設定時間ごとに行うステップＳ４～Ｓ８の一連の処理を繰り返すことで、設定時間ごとに心拍間隔データを算出してもよい。一の設定時間で行うステップＳ４～Ｓ８の一連の処理が、他の設定時間で行うステップＳ１～Ｓ３の処理と並行して行われてもよい。
【0067】
［作用］
以上の例によれば、設定時間において反射波から得られる一定数の第１～第Ｎの距離変動データが時系列データとして記憶部４３に蓄積されている。そして、時系列データから、検知部位に対応する信号強度が時系列に沿って並べられた時間変動データが得られる。時間変動データに含まれる信号強度は、生体Ｈの検知部位に対応しているので、検知対象である生体Ｈが動いていても、検知部位に対応した信号強度に基づいて生体データが得られる。その結果、非接触にて生体データを精度良く検知することが可能となる。
以上の例によれば、設定時間において反射波から得られる一定数の第１～第Ｎの距離変動データが時系列データとして記憶部４３に蓄積されている。そして、時系列データから、検知部位に対応する信号強度が時系列に沿って並べられた時間変動データが得られる。時間変動データに含まれる信号強度は、生体Ｈの検知部位に対応しているので、検知対象である生体Ｈが動いていても、検知部位に対応した信号強度に基づいて生体データが得られる。その結果、非接触にて生体データを精度良く検知することが可能となる。
【0068】
検知対象である生体Ｈが動作している条件下で動きを考慮せずに時間変動データを得る場合、ある基準時間における生体Ｈからの反射位置を特定し、当該位置を測定位置（基準値）として、測定位置における信号強度の時間変動を抽出することになる。この場合、生体Ｈの動作により測定位置において検知部位の情報が含まれていない期間があるので、動きによって生体情報を正確に検知し難い。これに対して、上記の例では、第ｎの距離変動データにおける検知部位に対応する距離成分（レンジビン）が基準値に近づくように第ｎの距離変動データの距離成分を補正する補正処理を第１～第Ｎの距離変動データに対して行うことにより、補正時系列データが得られる。そして、時間変動データ生成部４５は、強度取得処理を実行する際に、補正時系列データから基準値における信号強度を対応強度として取得する。この場合、生体Ｈの動きが信号強度のピークでの距離成分の変動に現れ、当該ピーク付近での検知部位に対応する距離成分を基準値に近づけるように補正が行われることで、時系列データにおいて生体Ｈによる動きが補正される。そのため、生体Ｈの動きによる影響が低減された信号強度の時間変動データが得られ、当該時間変動データに基づいて生体データが取得される。その結果、生体Ｈが動いていたとしても、生体の生体データを精度良く検知することが可能となる。
【0069】
以上の例によれば、補正処理は、第１～第Ｎの距離変動データのうち所定時間内に含まれる第ｐ～第ｑ（ｐ及びｑはそれぞれ、１≦ｐ＜ｑ≦Ｎを満たす自然数）の距離変動データにおける信号強度のピークでの距離成分に基づいて得られる近似線が基準値に近づくように第ｐ～第ｑの距離変動データの距離成分を補正することを含んでいる。生体Ｈの動線軌跡は、距離変動データに含まれるピークでの距離成分の近似線によって表すことができる。これにより、距離変動データにおいて、マルチパス等の外乱によるピークが含まれていても、当該ピークは近似線から離れているので除外することができる。そのため、生体データの高精度な検知をより確実に行うことが可能となる。
【0070】
以上の例によれば、補正処理は、時系列データを時系列に沿って順に分割することで得られる第１～第Ｍのブロックのうち第ｍのブロックＢｍにおいて、第１～第Ｎの距離変動データのうち第ｍのブロックに含まれる複数の距離変動データにおける信号強度のピークでの距離成分に基づいて得られる近似線が、第１～第ＭのブロックＢ１～ＢＭのそれぞれについて定められる任意の基準値に近づくように複数の距離変動データの距離成分を補正するブロック補正処理を、第１～第ＭのブロックＢ１～ＢＭに対して行うことを含んでいる。この場合、ブロックごとに、生体Ｈの動線軌跡を表す近似線を算出するので、近似線を生体Ｈの動線軌跡により近づけることができる。そのため、生体データのより高精度な検知を行うことが可能となる。
【0071】
以上の例によれば、生体データ生成部４６は、時間変動データに対して多重解像度解析を行うことで、生体Ｈの検知部位における生体データを生成するように構成されている。この場合、多重解像度解析により、時間変動データから生体データに対応する周波数範囲の波形が抽出される。そのため、時間変動データに生体データ以外の成分が含まれていても、生体データを取得することが可能となる。
【0072】
以上の例によれば、生体データ生成部４６は、時間変動データに対して多重解像度解析を行うことで、時間変動データを複数の分解波形に分解し、複数の分解波形のうち生体Ｈの拍動に対応する拍動信号に基づいて生体Ｈの拍動に関する拍動データを生成するように構成されている。この場合、多重解像度解析により、時間変動データから拍動に対応する周波数範囲の波形が抽出される。そのため、時間変動データに生体Ｈの拍動以外の成分が含まれていても、拍動データを取得することが可能となる。
【0073】
以上の例によれば、生体データ生成部４６は、ＡＲモデルを用いて拍動信号をスペクトル解析することにより、生体Ｈの心拍間隔の時間変動を示す心拍間隔データを生成するように構成されている。この場合、拍動信号に対してＡＲモデルを用いたスペクトル解析が行われることで、拍動信号に含まれるノイズが除去されるので、心拍間隔を精度良く算出することができる。
【0074】
以上の例によれば、受信機１２は、生体の拍動に対応する周波数を基準として、反射波をオーバサンプリングしうる。すなわち、受信機１２から出力された受信データに基づいて信号強度算出部４２が算出する第１～第Ｎの距離変動データは、生体の拍動に対応する周波数を基準として、反射波をオーバサンプリングすることで得られたデータである。この場合、サンプリング周期が小さくなり、より多くのデータが得られるので、時間変動データから拍動信号を抽出することが容易である。
【0075】
以上の例によれば、生体データ生成部４６は、複数の分解波形のうち、拍動とは異なる生体Ｈの動きに対応する信号を除外したうえで、心拍間隔データを生成するように構成されている。この場合、拍動とは異なる生体Ｈの動きによる情報が時間変動データから除外されるので、拍動データを精度良く検知することができる。
【0076】
以上の例によれば、反射波に応じた受信データを生体データ取得装置２０に無線通信を用いて送信するように構成された通信機１３が備えられる。無線電波の発信及び反射波の受信を行う送受信装置１０に搭載可能な制御部（メモリ）の容量が小さい場合があるが、別体の生体データ取得装置２０内でデータの蓄積と処理を行うことで、大容量のデータの取扱いができる。そのため、長時間にわたり生体データの検知を継続することが可能となる。
【0077】
以上に説明した実施形態は一例であり、特許請求の範囲及びその要旨を逸脱しない範囲で種々の変形を上記の実施形態に加えてもよい。
【0078】
（１）生体データ取得システム１において、送受信装置１０（通信機１３）と生体データ取得装置２０（通信機２３）との間は、有線で通信可能に接続されていてもよい。
【0079】
（２）生体データ取得装置２０は、送受信装置１０と同じ建物内に設けられてもよく、送受信装置１０とは異なる建物に取り付けられていてもよい。当該異なる建物としては、例えば、生体Ｈと別居している家族の家、健康管理サービスを提供する施設、救急医療サービスを提供する施設（日本において消防署）、病院などが挙げられる。これらの場合、生体Ｈとは離れた場所にいても、生体Ｈの状態を知ることができる。そのため、生体Ｈが単身生活者の場合であっても、生体Ｈの状態を見守ることが可能となる。生体データ取得装置２０は、取得した生体データから、生体Ｈの状態（異常の有無、健康状態）等を推定する機能を有していてもよい。
【0080】
（３）生体データ取得システム１ではなく、図１１に示されるように、生体データ取得装置２０が、発信機１１及び受信機１２を更に有していてもよい。換言すると、生体データを取得する装置（生体データ取得装置２０）内に、無線電波の送受を行う発信機１１及び受信機１２が内蔵されていてもよい。図１１に示される生体データ取得装置２０は、制御部２４Ａと、発信機１１と、受信機１２とを有する点において、図３に示される生体データ取得装置２０と相違する。制御部２４Ａは、発信機１１及び受信機１２を制御する機能を有しており、図１２に示されるように、データ取得部４１に代えて送受信処理部４８を有する点において制御部２４と相違する。送受信処理部４８は、送受信処理部３１と同様の機能を有していてもよい。図１１に示される生体データ取得装置２０を備える生体データ取得システム１では、無線電波の発信及び受信を行うための機能が、生体データ取得装置２０とは別体の送受信装置１０ではなく、生体データ取得装置２０自体に設けられているので、システムを簡略化することが可能となる。
【0081】
（４）生体データ取得システム１において、発信機１１及び受信機１２が別々の装置として構成されてもよい。つまり、生体データ取得システム１は、送受信装置１０に代えて、発信機１１及び一の制御部を有する発信装置と、受信機１２及び他の制御部を有する受信装置とを備えていてもよい。これらの発信装置、受信装置、及び生体データ取得装置２０は、互いに無線又は有線により通信可能に接続されていてもよい。発信機１１、受信機１２及び送受信処理部４８において、送信及び受信の双方に複数のアンテナ素子が配置されるＭＩＭＯ（Multiple Input Multiple Output）レーダ技術が用いられてもよい。
【0082】
（５）生体データ取得装置２０は、複数の生体Ｈについて、生体データ（心拍間隔データ）をそれぞれ取得してもよい。検知領域Ｒに複数の生体Ｈ（例えば二人の検知対象者）が存在する場合がある。例えば、一の検知対象者は送受信装置１０から遠ざかり、他の検知対象者は送受信装置１０に近づく場合を考えると、距離変動データから得られる複数のピークｐｎ（散乱点）は、図１３に示されるような２つの散乱点群Ｃｌ１，Ｃｌ２（信号強度と距離の情報）として示される。
【0083】
データ補正部４４は、１つのブロック内に複数の散乱点群が存在すること（１つのブロック内に複数人が存在すること）を検知した場合に、それぞれの散乱点群をクラスタ化してもよい。データ補正部４４は、クラスタ毎に一次補正、二次補正、及び三次補正を行って、クラスタ化された散乱点ごとに補正時系列データを１つずつ生成してもよい。時間変動データ生成部４５は、クラスタ化された散乱点ごとに得られた補正時系列データから時間変動データをそれぞれ生成してもよい。生体データ生成部４６は、クラスタ化された散乱点ごとに得られた時間変動データから生体データをそれぞれ生成してもよい。一の検知対象者に対する生体データの取得処理が完了した後に、他の検知対象者に対する生体データの取得処理が行われてもよい。
【0084】
図１４～図１８を参照して、検知領域Ｒ内に一の生体Ｈ１と別の生体Ｈ２が存在する場合の生体データの取得処理の一例を説明する。検知領域Ｒ内に生体Ｈ１，Ｈ２が存在する場合、信号強度算出部４２は、無線電波が生体Ｈ１及び生体Ｈ２において反射された反射波に基づいて、時刻ｔ１～ｔＮにおける第１～第Ｎの距離変動データをそれぞれ時間順に算出してもよい。無線電波が生体Ｈ１及び生体Ｈ２において反射された反射波は、生体Ｈ１のみから反射された反射波と、生体Ｈ２のみから反射された反射波とを含む。記憶部４３は、生体Ｈ１及び生体Ｈ２において反射された反射波に基づいて得られる第１～第Ｎの距離変動データを含む時系列データを記憶してもよい。距離変動データの蓄積単位を示す設定時間が、上述した例に加えて、数秒～数十秒程度の範囲に設定されてもよい。設定時間の一例は、１秒～１８０秒程度であってもよく、３秒～１２０秒程度であってもよく、４秒～６０秒程度であってもよい。設定時間の他の例は、１秒よりも短い時間（例えば、数ミリ秒から数百ミリ秒程度）であってもよい。
【0085】
図１４に示される時系列データは、約４秒間において得られる複数の距離変動データを含んでいる。当該時系列データは、被験者である生体Ｈ１及び生体Ｈ２が送受信装置１０に正対した測定環境において得られたデータである。より詳しくは、当該時系列データは、生体Ｈ１が送受信装置１０から約１．９ｍだけ離れた位置に静止し、且つ生体Ｈ２が送受信装置１０から約２．２ｍだけ離れた位置に静止した状態で測定が開始され、測定開始の約２秒後から、生体Ｈ１が送受信装置１０に徐々に近づくと共に、生体Ｈ２が送受信装置１０から徐々に遠ざかったときに取得されたデータである。
【0086】
データ補正部４４は、記憶部４３が記憶する時系列データ（所定時間分の第１～第Ｎの距離変動データ）に基づいて、生体Ｈ１について補正処理と、生体Ｈ２についての補正処理とをそれぞれ行ってもよい。データ補正部４４は、検出された生体の数に応じて、複数のピークｐｎ（散乱点）を公知の手法によりクラスタ化してもよい。データ補正部４４は、例えば、レンジビンごとの信号強度の変化に基づいて、各生体（生体Ｈ１及び生体Ｈ２それぞれ）の動作を追跡することで、複数のピークｐｎを生体ごとにクラスタ化してもよい。
【0087】
データ補正部４４は、例えば、第ｎの距離変動データにおける生体Ｈ１の検知部位に対応する距離成分が基準値に近づくように第ｎの距離変動データの距離成分を補正する上述の補正処理を行ってもよい。図１５（ａ）は、生体Ｈ１についてクラスタ化された複数のピークｐｎに基づき補正処理が行われて得られた補正時系列データを示している。なお、図１５（ｂ）に示されるように、二次補正が行われることで、観測時間によっては、ピークｐｎと一致しない反射点ｐｔが基準値（破線で示す位置）に揃えられる。
【0088】
時間変動データ生成部４５は、第ｎの距離変動データにおける信号強度のうち生体Ｈ１の検知部位に対応する対応強度を取得する強度取得処理を行ってもよい。時間変動データ生成部４５は、第１～第Ｎの距離変動データに対して、生体Ｈ１についての強度取得処理を行うことにより、対応強度が時系列に沿って並べられた時間変動データを生成してもよい。時間変動データ生成部４５は、例えば、生体Ｈ１についての補正時系列データ内の基準値において時系列に沿って並べられた信号強度の集合を、時間変動データとして得てもよい。
【0089】
図１６（ａ）は、図１４に示される補正前の時系列データにおいて、基準値（例えば、距離が１．９ｍ）での信号強度の時間変化を示している。図１６（ｂ）は、図１５（ａ）に示される補正後の時系列データにおいて、基準値での信号強度の時間変化を示している。図１６（ａ）に示す信号強度の時間変化では、測定開始の約２秒後から、生体Ｈ１の検知部位ではなく他の反射位置からの反射波の信号強度が得られている。図１６（ｂ）に示す信号強度の時間変化では、時系列データにおいて補正が行われているので、生体Ｈ１の検知部位からの反射波の信号強度が連続して得られている。
【0090】
データ補正部４４は、生体Ｈ１についての補正処理とは別に、第ｎの距離変動データにおける生体Ｈ２の検知部位に対応する距離成分が別の基準値に近づくように第ｎの距離変動データの距離成分を補正する補正処理を行ってもよい。図１７は、生体Ｈ２についてクラスタ化された複数のピークｐｎに基づき補正処理が行われて得られた補正時系列データを示している。
【0091】
時間変動データ生成部４５は、生体Ｈ１についての強度取得処理とは別に、第ｎの距離変動データにおける信号強度のうち生体Ｈ２の検知部位に対応する対応強度を取得する強度取得処理を行ってもよい。時間変動データ生成部４５は、第１～第Ｎの距離変動データに対して、生体Ｈ２についての強度取得処理を行うことにより、対応強度が時系列に沿って並べられた時間変動データを生成してもよい。時間変動データ生成部４５は、例えば、生体Ｈ２についての補正時系列データ内の基準値において時系列に沿って並べられた信号強度の集合を、時間変動データとして得てもよい。
【0092】
図１８（ａ）は、図１４に示される補正前の時系列データにおいて、基準値（例えば、距離が２．２ｍ）での信号強度の時間変化を示している。図１８（ｂ）は、図１７に示される補正後の時系列データにおいて、基準値での信号強度の時間変化を示している。図１８（ａ）に示す信号強度の時間変化では、測定開始の約２秒後から、生体Ｈ２の検知部位ではなく他の反射位置からの反射波の信号強度が得られている。図１８（ｂ）に示す信号強度の時間変化では、時系列データにおいて補正が行われているので、生体Ｈ２の検知部位からの反射波の信号強度が連続して得られている。なお、生体Ｈ２は、送受信装置１０から徐々に遠ざかっているので、後半部分において信号強度が全般的に低下している。
【0093】
生体データ生成部４６は、生体Ｈ１についての時間変動データに基づいて、生体Ｈ１の検知部位における生体データを生成してもよい。生体データ生成部４６は、例えば、生体Ｈ１についての時間変動データに対して多重解像度解析を行うことにより、生体Ｈ１の検知部位における生体データを生成してもよい。生体データ生成部４６は、生体Ｈ２についての時間変動データに基づいて、生体Ｈ２の検知部位における生体データを生成してもよい。生体データ生成部４６は、例えば、生体Ｈ２についての時間変動データに対して多重解像度解析を行うことにより、生体Ｈ２の検知部位における生体データを生成してもよい。
【0094】
以上の例の生体データ取得装置２０によれば、検知領域Ｒに複数の生体が存在していても、生体ごとに生体データを監視することが可能となる。
【0095】
（６）検知対象の生体Ｈを人間としたときに、生体Ｈが同じ場所に留まって動作している場合（例えば、机での作業又は自動車運転時等）がある。この場合、データ取得部４１により取得される時系列データは、図６（ｂ）に示される状態に近くなる。このような場合、データ補正部４４は、ブロックごとの近似線の算出（一次補正）を行うことなく、二次補正及び三次補正を行ってもよい。動きが限定的な場合において、データ補正部４４は、散乱範囲に含まれるレンジビンの数を増加してもよい。
【0096】
（７）生体データ生成部４６は、上述した多重解像度解析に代えて、時間変動データに対して高速フーリエ変換（FFT：Fast Fourier Transform）又はウェーブレット変換を行うことにより、平均呼吸数及び平均心拍数の少なくとも一方を生体データとして算出（生成）してもよい。拍動信号は呼吸信号に比べて微弱だが、例えば、観測時間を数分程度と長く設定することで、動いている生体Ｈの平均呼吸数及び心拍数（呼吸スペクトル及び拍動スペクトル）を推定することが可能となる。
【0097】
（８）データ補正部４４は、二次補正を行った後に、ブロックごとに、第ｎの距離変動データにおける信号強度のピークでの距離成分と当該ブロックでの基準値との間の偏差に基づく誤差を算出してもよい。データ補正部４４は、算出した誤差が所定の閾値よりも小さいかどうかを判断してもよい。データ補正部４４は、当該誤差が所定の閾値よりも小さくなるまで、生成処理、一次補正及び二次補正を繰り返し行ってもよい。データ補正部４４は、例えば、上記誤差として、第１～第Ｎの距離変動データそれぞれにおけるピークでの距離成分と、基準値との間の平均平方二乗誤差（Root Mean Square Error）を算出してもよい。データ補正部４４は、ブロックごとの上記誤差が所定の閾値よりも小さくなった場合に、三次補正を行ってもよい。基準値は、例えば、生体Ｈの現実の位置であってもよい。
【0098】
（９）データ補正部４４は、誤差の算出を行うことなく生成処理、一次補正及び二次補正を行い、三次補正を行った後に、１つのブロック（例えば、第１のブロックＢ１）の基準値を基準として、上記と同様に誤差を算出してもよい。図１９に示される生体データ取得方法は、ステップＳ４１が更に実行される点において、図１０に示される取得方法と相違する。データ補正部４４は、ステップＳ４の時系列データの補正を行った後に、上記誤差が閾値εよりも小さいかどうかを判断する（ステップＳ４１）。ステップＳ４１において、上記誤差が閾値ε以上であると判断された場合（ステップＳ４１：ＮＯ）、データ補正部４４は、ステップＳ４の補正処理を再度実行する。ステップＳ４１において、上記誤差が閾値εよりも小さいと判断された場合（ステップＳ４１：ＹＥＳ）、制御部２４は、ステップＳ５～Ｓ８の一連の処理を実行する。
【0099】
（１０）生体データ取得システム１は乗物に設けられていてもよい。生体データ取得システム１は、乗物内に設定された検知領域Ｒに存在する生体Ｈを検知対象として生体データを取得してもよい。乗物の一例は、自動車、自動二輪車（オートバイ）、鉄道車両、船、飛行機、及びヘリコプターを含む。乗物内の検知対象は、例えば、当該乗物の操縦者（運転者）及び同乗者を含む。図２０には、生体データ取得システム１が設けられた乗物６０が例示されている。
【0100】
図２０に示されるように、生体データ取得システム１は、乗物６０の車室内に設けられてもよい。例えば、生体データ取得システム１の送受信装置１０は、検知領域Ｒが乗物６０の運転席を含むように配置されてもよい。一例として、送受信装置１０は、図２０に示されるように、フロントガラス５１に設けられてもよい。送受信装置１０は、フロントガラス５１に代えて、ダッシュボート５２、ハンドル部材５３、ルームミラー５４、及びルーフパネル５５のうちのいずれかの箇所に設けられてもよい。あるいは、送受信装置１０は、運転席に近い側のピラー５６（フロントガラス５１の横に位置するＡピラー）に設けられてもよい。生体データ取得システム１の生体データ取得装置２０も、乗物６０内に配置されてもよい。発信機１１及び受信機１２を有する生体データ取得装置２０（図１１参照）が、上記のいずれかの箇所に配置されてもよい。
【0101】
生体データ取得システム１が、乗物６０に設けられた場合でも、生体データ取得装置２０において、図１０又は図１９に示されるデータ取得方法の一連の処理が実行されてもよい。この場合、ステップＳ１～Ｓ３が繰り返し実行されて得られる第１～第Ｎの距離変動データはそれぞれ、無線電波が乗物６０内の生体Ｈにおいて反射された反射波に基づくデータである。乗物６０内の生体Ｈからの反射波に基づく時系列データに対して、ステップＳ４，Ｓ５の処理が施されることで時間変動データが生成されてもよい。
【0102】
生体データ生成部４６は、時間変動データに対して、上述の例（７回のウェーブレット変換）とは異なり８回のウェーブレット変換を行うことにより、複数の分解波形を生成してもよい。生体データ生成部４６は、例えば、拍動信号を抽出する際に、呼吸信号に加えて、乗物６０に由来して生体Ｈに加わる振動に対応する非拍動信号を除外してもよい。生体データ生成部４６は、乗物６０に由来する振動の固有周波数が含まれる周波数範囲を有する分解波形を除外してもよい。生体データ生成部４６は、呼吸信号及び乗物６０に由来する非拍動信号を除外した後の残りの複数の分解波形を再合成することで、生体Ｈの拍動に応じた拍動信号を抽出してもよい。生体データ生成部４６は、当該拍動信号に対して、上述と同様に、ＡＲモデルを用いたスペクトル解析を行うことで、生体データ（例えば、心拍間隔データ）を生成してもよい。
【0103】
図２１には、乗物６０内の生体Ｈから取得された心拍間隔データの一例が示されている。図２１に示される心拍間隔データは、道路上を走行させている乗物６０（自動車）内の運転者の心拍変動（Ｈｚ）を測定した結果であり、図２１には２回の測定を行った結果が示されている。なお、心拍変動（Ｈｚ）は、運転者の心拍間隔（ｓ）の逆数に対応する。図２１において、「測定値１」は、７回のウェーブレット変換を行い、呼吸信号を除いて得られた１回目の測定値であり、「測定値２」は、８回のウェーブレット変換を行い、呼吸信号と乗物６０に由来する非拍動信号とを除外して得られた２回目の測定値である。「参照値１」及び「参照値２」は、それぞれ測定値１及び測定値２に対応しており、比較のために接触式のセンサで得られた計測値である。参照値１及び参照値２は、測定値１及び測定値２をそれぞれ取得する際に、運転者の指先にＰＰＧ（光学式心拍センサ）を取り付けて計測されている。なお、送受信装置１０が運転席に近い側のピラー５６に取り付けられて、生体データの測定が行われている。
【0104】
図２１に示される測定結果から、乗物６０内においても、図９に示される結果と同様に、生体Ｈの心拍間隔データが取得されることがわかる。また、測定値１と参照値１とを比較した場合の両者の相関は０．９５以下であり、測定値２と参照値２とを比較した場合の両者の相関は０．９６～０．９７であった。このことから、生体Ｈに接触させて計測した参照値を真値とした場合に、測定値１に比べて測定値２に含まれる誤差が小さいことがわかる。すなわち、乗物６０に由来する非拍動信号を除外することで、乗物６０内の生体Ｈの生体データがより精度良く検出されている。
【0105】
以上の例の生体データ取得装置２０を有する乗物６０によれば、乗物６０内の１又は複数の生体に関する生体データを監視することが可能となる。生体データ生成部４６が、乗物６０に由来する非拍動信号を除外したうえで、拍動データを生成するので、乗物６０に由来する情報が時系列データに含まれていても、拍動データを精度良く検知することができる。なお、上述の補正処理が行われることで、運転者の運転時における動作（例えば、左右方向の確認動作、又は各種の操作に伴う動作）があっても、精度良く生体データを検知できる。
【0106】
（１１）生体データ取得装置２０は、上記生体以外の物体からの不要な反射波の影響を抑制するための処理を実行することで、上記第１～第Ｎの距離変動データを含む上記時系列データを生成する。不要な反射波は、例えば、壁からの反射波、及び、什器（例えば、机、又は椅子）からの反射波を含む。以下では、不要な反射波の影響を低減するための処理を「抑圧処理」と称する。図２２には、上記抑圧処理を説明するための模式図が示されている。
【0107】
発信機１１及び受信機１２は、ＦＭ（Frequency Modulation）レーダを構成してもよい。発信機１１は、上記無線電波である複数の送信波ＴＷを検知領域Ｒに向けて順に発信するように構成されていてもよい。発信機１１は、複数の送信波ＴＷを時系列順に並ぶ異なるタイミングでそれぞれ発信してもよい。発信機１１は、一定の時間間隔で複数の送信波ＴＷそれぞれを発信してもよい。発信機１１は、例えば、０．１秒間において、４０個～１４０個の送信波ＴＷを発信するように構成されている。
【0108】
受信機１２は、複数の送信波ＴＷそれぞれが検知領域Ｒにおいて反射されて得られる複数の受信波ＲＷを受信するように構成されている。複数の受信波ＲＷは、上述した無線電波の反射波である。発信機１１及び受信機１２は、図２２に示されるように、１つの送信波ＴＷの送信と、１つの受信波ＲＷの受信とを交互に繰り返す。受信機１２は、検知領域Ｒを距離方向に沿って分割して得られる第１～第Ｋのレンジビン（上記複数のレンジビン）に含まれる１つのレンジビンごとに、信号を抽出してもよい。Ｋは、２以上の自然数であり、複数のレンジビンの数を示す。図２２では、第１～第Ｋのレンジビンそれぞれの距離（位置）が、ｘ１，ｘ２，ｘ３…，ｘＫによって示されている。
【0109】
上記抑圧処理では、任意に設定されるフレーム周期単位で（１つのフレーム周期ごとに）、第ｎの距離変動データが生成される。１つのフレーム周期の長さＴは、０．０１秒間～０．２秒間程度であってもよい。１つのフレーム周期の長さＴが０．１秒間に設定されている場合、４０個～１４０個の送信波ＴＷに対応する４０個～１４０個の受信波ＲＷに基づいて、第ｎの距離変動データが生成されてもよい。図２２に示される例では、１つのフレーム周期内において、送信波ＴＷ１，ＴＷ２，・・・，ＴＷＪが送信され、受信波ＲＷ１，ＲＷ２，・・・，ＲＷＪが受信される。Ｊは、２以上の自然数であり、例えば、フレーム周期の長さＴが０．１秒間である場合に、Ｊは、４０～１４０であってもよい。
【0110】
信号強度算出部４２は、１つの受信波ＲＷに応じた受信データに基づいて、第１～第Ｋのレンジビンそれぞれにおける信号強度を算出してもよい。信号強度算出部４２は、１つのフレーム周期内で受信する複数の受信波ＲＷに応じた受信データに基づいて、第１～第Ｊの受信信号データを生成（算出）してもよい。信号強度算出部４２は、受信波ＲＷｊに応じた受信データに基づいて、第１～第Ｊの受信信号データのうちの第ｊの受信信号データを生成してもよく、ｊは、１～Ｊの自然数である。
【0111】
第１～第Ｊの受信信号データは、１つのフレーム周期内の時系列順に並ぶ異なる時刻において検知領域Ｒで反射された反射波に基づいて得られる。上記反射波は、広帯域の無線電波又は超広帯域の無線電波である。第ｊの受信信号データ（第１～第Ｊの受信信号データそれぞれ）は、距離に対する信号強度の変化を示すデータである。信号強度は、複素数（又は、ＩＱ信号）によって示されてもよい。第ｊの受信信号データは、第１～第Ｋのレンジビンにおける信号強度ａｊ１～ａｊＫを含む。信号強度ａｙｚのうち「ｙ」は第ｙの受信信号データの信号強度であることを示しており、信号強度ａｙｚのうち「ｚ」は第ｚのレンジビンでの信号強度であることを示している。Ｋは、２以上の自然数であり、ｙは、１～Ｊの自然数であり、ｚは、１～Ｋの自然数である。第ｊの受信信号データの第１のレンジビンにおける信号強度がａｊ１によって表され、第ｊの受信信号データの第２のレンジビンにおける信号強度がａｊ２によって表される。第ｊの受信信号データの第３のレンジビン～第（Ｋ－１）のレンジビンにおける信号強度がそれぞれａｊ３～ａｊ（Ｋ－１）によって表される。第ｊの受信信号データの第Ｋのレンジビンにおける信号強度がａｊＫによって表される。
【0112】
第１の受信信号データは、第１～第Ｋのレンジビンにおける信号強度ａ１１，ａ１２，・・・，ａ１Ｋを含む。第２の受信信号データは、第１～第Ｋのレンジビンにおける信号強度ａ２１，ａ２２，・・・，ａ２Ｋを含む。第Ｊの受信信号データは、第１～第Ｋのレンジビンにおける信号強度ａＪ１，ａＪ２，・・・，ａＪＫを含む。記憶部４３は、信号強度算出部４２によって生成された第１～第Ｊの受信信号データを記憶してもよい。信号強度算出部４２は、第１～第Ｊの受信信号データに対して上記抑圧処理を行うことで、第ｎの距離変動データを生成する。
【0113】
上記抑圧処理は、係数算出処理と、データ生成処理とを含んでもよい。信号強度算出部４２は、係数算出処理として、信号強度ａ１ｉ～ａＪｉを周波数変換して得られた結果に基づいて係数ｖｉを算出する処理を実行してもよく、ｉは１～Ｋの自然数である。信号強度算出部４２は、信号強度ａ１ｉ～ａＪｉを含む時間データを周波数変換してもよい。信号強度算出部４２は、係数算出処理を実行する際に、第１～第Ｋのレンジビンに含まれる１つのレンジビンごとに、第１～第Ｊの受信信号データの当該１つのレンジビンにおける複数の信号強度を含む時間データを周波数変換してもよい。信号強度算出部４２は、第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉを周波数変換して得られた結果に基づいて係数ｖｉを算出する。
【0114】
信号強度算出部４２は、第１のレンジビンにおいて、信号強度ａ１１～ａＪ１を含む時間データを周波数変換して得られた結果に基づいて係数ｖ１を算出してもよい。信号強度算出部４２は、第２のレンジビンにおいて、信号強度ａ１２～ａＪ２を含む時間データを周波数変換して得られた結果に基づいて係数ｖ２を算出してもよい。信号強度算出部４２は、第３～第Ｋのレンジビンそれぞれにおいても、第１のレンジビンと同様に係数ｖ３～ｖＫを算出してもよい。信号強度算出部４２は、上記１つのレンジビンごとに、周波数変換により得られた結果に基づいて係数を算出してもよい。上記時間データは、時間軸と信号強度軸を含むデータであり、１つのレンジビンでの信号強度の時間変化を示している。
【0115】
信号強度算出部４２は、高速フーリエ変換（Fast Fourier Transform）を用いて信号強度ａ１ｉ～ａＪｉを周波数変換してもよい。信号強度算出部４２は、上記１つのレンジビンごとに、複数の信号強度を周波数変換して得られる周波数スペクトルにおいて最も大きい値を有する周波数成分の周波数を抽出してもよい。信号強度算出部４２は、第１～第Ｋのレンジビンにおいてそれぞれ抽出された複数の周波数を正規化することで、上記係数を算出してもよい。信号強度算出部４２は、上記複数の周波数のうちの最も大きい周波数が１となるように正規化を行ってもよい。
【0116】
信号強度算出部４２は、データ生成処理として、第１～第Ｊの受信信号データのいずれか１つにおける複数の信号強度に係数ｖｉを乗算することで、第ｎの距離変動データを算出する処理を実行してもよい。信号強度算出部４２は、第ｊの受信信号データ（例えば、第１の受信信号データ）について、第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉに係数ｖｉを乗算することで、第ｎの距離変動データを算出する。第ｎの距離変動データの第１のレンジビンにおける信号強度は、信号強度ａｊ１に係数ｖ１を乗算することで得られる。第ｎの距離変動データの第２のレンジビンにおける信号強度は、信号強度ａｊ２に係数ｖ２を乗算することで得られる。第ｎの距離変動データの第３～第（Ｋ－１）のレンジビンにおける信号強度は、信号強度ａｊ３～ａｊ（Ｋ－１）に係数ｖ３～ｖ（Ｋ－１）を乗算することでそれぞれ得られる。第ｎの距離変動データの第Ｋのレンジビンにおける信号強度は、信号強度ａｊＫに係数ｖＫを乗算することで得られる。信号強度算出部４２は、第１の受信信号データについて、第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａ１ｉに係数ｖｉを乗算することで、第ｎの距離変動データを算出してもよい。
【0117】
信号強度算出部４２は、１つのフレーム周期において第ｎの距離変動データを算出した後に、次のフレーム周期において得られる第１～第Ｊの受信信号データに対して上記抑圧処理を実行することで、第（ｎ＋１）の距離変動データを算出してもよい。図２３は、第１～第Ｎの距離変動データを含む時系列データを生成する手順の一例を示すフローチャートである。最初に、信号強度算出部４２は、受信波ＲＷ１に応じた受信データに基づいて、第１～第Ｋのレンジビンそれぞれにおける信号強度ａ１１～ａ１Ｋを算出することで、距離に対する信号強度の変化を示す第１の受信信号データを取得する（図２３のステップＳ１１）。
【0118】
次に、記憶部４３は、ステップＳ１１で取得された第１の受信信号データを記憶する（図２３のステップＳ１２）。そして、制御部２４は、１回目のステップＳ１１の開始時刻から１つのフレーム周期が経過したかどうかを判断する（図２３のステップＳ１３）。１つのフレーム周期を経過していないと判断された場合（ステップＳ１３：ＮＯ）、制御部２４は、ステップＳ１１～Ｓ１３を繰り返す。２回目のステップＳ１１では、信号強度算出部４２が、第２の受信信号データを生成し、Ｊ回目のステップＳ１１では、信号強度算出部４２が、第Ｊの受信信号データを生成する。これにより、１つのフレーム周期において、第１～第Ｊの受信信号データを含むデータ群が、記憶部４３に記憶される。
【0119】
１つのフレーム周期が経過したと判断された場合（ステップＳ１３：ＹＥＳ）、信号強度算出部４２は、第１～第Ｋのレンジビンにおいて係数ｖ１～ｖＫを算出する（図２３のステップＳ１４）。信号強度算出部４２は、第１～第Ｋのレンジビンの全てそれぞれについて、信号強度ａ１ｉ～ａＪｉを含む時間データを周波数変換して得られた結果に基づいて、係数ｖｉを算出してもよい。そして、信号強度算出部４２は、第１の受信信号データの信号強度ａ１ｉに係数ｖｉを乗算することで、第１の距離変動データを生成してもよい（図２３のステップＳ１５）。
【0120】
次に、記憶部４３は、ステップＳ１５で生成した第１の距離変動データを記憶する（図２３のステップＳ１６）。そして、制御部２４は、初期時刻から上記設定時間が経過したかどうかを判断する（図２３のステップＳ１７）。上記設定時間を経過していないと判断された場合（ステップＳ１７：ＮＯ）、制御部２４が実行する処理は、ステップＳ１１に戻る。以降、ステップＳ１１～Ｓ１３での繰り返し処理とステップＳ１４～Ｓ１７とを含む一連の処理が繰り返されることで、第２～第Ｎの距離変動データが生成及び記憶される。
【0121】
上記設定時間が経過したと判断された場合（ステップＳ１７：ＹＥＳ）、１つの設定時間において時系列データを取得する処理が終了する。そして、制御部２４は、図１０に示されるステップＳ４～Ｓ８の一連の処理を実行することで、抑圧処理が施されたうえで取得された時系列データに基づいて、心拍間隔データ等の生体データを生成する。
【0122】
上記抑圧処理では、１つのレンジビンにおいて、時系列に並ぶ複数の信号強度を含む時間データを周波数変換して得られた結果に基づき係数が算出される。検知対象物である生体以外の壁等が存在するレンジビンでは、複数の信号強度の時間変動が小さいので、時間データを周波数変換すると低い周波数成分が得られる。動作している生体が存在するレンジビンでは、生体の動作（動き）により複数の信号強度の時間変動が比較的大きくなり、時間データを周波数変換すると高い周波数成分が得られる。そのため、検知対象外の物体が存在するレンジビンでの係数を小さくし、生体が存在するレンジビンでの係数を大きくすることができる。その結果、第ｎの距離変動データにおいて、検知対象外の物体からの反射波に基づく信号強度と、生体からの反射波に基づく信号強度との差を大きくすることができる。従って、生体データのより高精度な検知を行うことが可能となる。
【0123】
図２４（ａ）に示される時系列データは、約６０秒間において得られる複数の距離変動データを含んでおり、上記抑圧処理に代えて、後述の積分処理が行われることで得られている。図２４（ｂ）に示される時系列データは、約６０秒間において得られる複数の距離変動データを含んでおり、上記抑圧処理が行われることで得られている。図２４（ａ）に示される時系列データ及び図２４（ｂ）に示される時系列データの間では、距離変動データの算出方法が異なるが、同じ受信データに基づいて複数の距離変動データが算出されている。
【0124】
上記受信データは、被験者である生体Ｈ１及び生体Ｈ２が、検知領域Ｒにおいて歩行している測定環境において得られたデータである。上記測定環境では、生体Ｈ１が、送受信装置１０から約１．２ｍだけ離れた位置と、約４．５ｍだけ離れた位置との間で往復するように歩行している。上記測定環境では、生体Ｈ２が、送受信装置１０から約４．０ｍだけ離れた位置と、約５．３ｍだけ離れた位置との間で往復するように歩行している。上記積分処理では、１つのレンジビンごとに、第１～第Ｊの受信信号データの当該１つのレンジビンにおける複数の信号強度を含む時間データを積分することで得られる積分値に基づき、第ｎの距離変動データが生成されている。
【0125】
図２４（ａ）に示される時系列データでは、送受信装置１０から約３．８ｍだけ離れた位置、約４．７ｍだけ離れた位置、及び、約５．４ｍだけ離れた位置それぞれからの不要な反射波に基づく信号が含まれている。図２４（ａ）に示される時系列データでは、検知対象外の物体からの反射波に基づく信号強度と、生体Ｈ１又は生体Ｈ２からの反射波に基づく信号強度との差が小さい。図２４（ｂ）に示される時系列データでは、上記不要な反射波に基づく信号の強度が、ほぼゼロである。そのため、検知対象外の物体からの反射波に基づく信号強度と、生体Ｈ１又は生体Ｈ２からの反射波に基づく信号強度との差が大きい。
【0126】
生体Ｈ１又は生体Ｈ２が少なくとも一部の期間で静止していても、生体の動作には無意識下の動作が含まれるので、生体に基づく複数の信号強度の時間変動が、固定物に基づく信号強度の時間変動に比べて大きくなり得る。そのため、歩行等の意図的な動作をしている場合と同様に、抑圧処理が実行されることで、検知対象外の物体からの反射波に基づく信号強度と、生体Ｈ１又は生体Ｈ２からの反射波に基づく信号強度との差が大きくなると考えられる。
【0127】
壁等が存在するレンジビンでの信号強度は雑音と同程度で小さく、周波数変換すると、０Ｈｚ付近において周波数成分が最も大きくなる。動作している生体が存在するレンジビンでは、レンジビンでの時間データを周波数変換すると、０Ｈｚ以外の周波数成分が最も大きくなる結果が得られる。上述の例では、壁等からの反射波に基づく係数ｖがほぼゼロとなるが、抑圧処理において、生体Ｈの動きよりも小さいレンジビンでの係数ｖを０に設定してもよい。例えば、周波数変換して得られる成分が最も大きい周波数が所定の閾値以下となるレンジビンでの係数ｖを０に設定し、所定の閾値よりも大きいレンジビンでの係数ｖを１に設定してもよい。
【0128】
（１２）上記積分処理が実行されることで、一部のノイズ成分が低減される。信号強度算出部４２は、抑圧処理を実行する際に、係数算出処理とデータ生成処理とに加えて、積分処理を行う。信号強度算出部４２は、積分処理を実行する際に、第１～第Ｊの受信信号データのうちの複数の受信信号データを積分して積分信号データを生成する。積分信号データは、第１～第Ｋのレンジビンにおける信号強度ｂ１～ｂＫを含む。積分信号データの第１のレンジビンにおける信号強度がｂ１によって表され、積分信号データの第２のレンジビンにおける信号強度がｂ２によって表される。積分信号データの第Ｋのレンジビンにおける信号強度がｂＫによって表される。
【0129】
信号強度算出部４２は、第１～第ｒの受信信号データを積分して積分信号データを生成してもよく、ｒは、２～Ｊの自然数であり、ｒはＪよりも小さくてもよい。信号強度算出部４２は、レンジビンごとに、第１～第ｒの受信信号データの当該レンジビンにおける複数の信号強度を積分することで得られる積分値に基づいて、信号強度ｂ１～ｂＫを算出してもよい。信号強度算出部４２は、第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａｒｉを積分することで得られる積分値に基づいて、信号強度ｂｉを算出してもよい。
【0130】
信号強度算出部４２は、第１のレンジビンにおいて、信号強度ａ１１～ａｒ１を積分することで得られる積分値に基づいて、信号強度ｂ１を算出してもよい。信号強度算出部４２は、第２のレンジビンにおいて、信号強度ａ１２～ａｒ２を積分することで得られる積分値に基づいて、信号強度ｂ２を算出してもよい。信号強度算出部４２は、第３～（Ｋ－１）のレンジビンにおいて、第１のレンジビンと同様に、係数ｂ３～ｂ（Ｋ－１）を算出してもよい。信号強度算出部４２は、第Ｋのレンジビンにおいて、信号強度ａ１Ｋ～ａｒＫを積分することで得られる積分値に基づいて、信号強度ｂＫを算出してもよい。
【0131】
信号強度算出部４２は、上記データ生成処理を実行する際に、積分信号データについて第１～第Ｋのレンジビンのそれぞれにおいて信号強度ｂｉと係数ｖｉとを乗算して、第ｎの距離変動データを算出してもよい。第ｎの距離変動データの第１のレンジビンにおける信号強度は、信号強度b１に係数ｖ１を乗算することで得られてもよい。第ｎの距離変動データの第２のレンジビンにおける信号強度は、信号強度ｂ２に係数ｖ２を乗算することで得られてもよい。第ｎの距離変動データの第Ｋのレンジビンにおける信号強度は、信号強度ｂＫに係数ｖＫを乗算することで得られてもよい。
【0132】
以上の例では、積分処理によって一部の受信信号データのＳＮ比を向上させたうえで、信号強度と係数ｖ１～ｖＫとの乗算が行われる。これにより、第ｎの距離変動データにおいて、検知対象外の物体からの反射波に基づく信号強度と、生体Ｈ１又は生体Ｈ２からの反射波に基づく信号強度との差を大きくすることができる。従って、生体データのより高精度な検知を行うことが可能となる。
【0133】
（１３）検知領域Ｒが室内に設定される場合、１つのフレーム周期内での生体の動きを無視できる場合が多い。しかしながら、生体Ｈの動作が速くなり、１つのフレーム周期内において生体の動きが無視できなくなると、生体データの検出精度が低下するおそれがある。信号強度算出部４２は、抑圧処理を実行することの前にフレーム設定処理を実行してもよい。上記フレーム設定処理は、生体Ｈの動作の速度を示す動作情報を取得することと、動作情報が示す速度に基づいてフレーム周期の長さを設定することとを含んでもよい。
【0134】
信号強度算出部４２は、生体データ取得装置２０以外の装置によって得られた上記動作情報を、その装置から取得してもよい。生体データ取得装置２０が、上記動作情報を算出してもよい。信号強度算出部４２は、生体Ｈの動作の速度が、所定の閾値以下である場合には、フレーム周期の長さを標準値（例えば、０．１秒間）に設定してもよい。信号強度算出部４２は、生体Ｈの動作の速度が上記閾値よりも大きい場合に、フレーム周期の長さを標準値よりも小さい値（例えば、０．０１秒間）に設定してもよい。信号強度算出部４２は、生体Ｈの動作の速度と上記閾値との差によって、フレーム周期の長さの調節量を変化させてもよい。この例では、生体Ｈの動作が速くなった場合に、フレーム周期を短くすることができる。フレーム周期が短いと第１～第Ｎの距離変動データに生体の動きが反映されやすい。従って、生体データのより高精度な検知を行うことが可能となる。
【0135】
（１４）抑圧処理において、影響が抑制される不要な反射波は、壁及び什器等の動作しない物体（以下、「固定物」という。）からの反射波に代えて又は加えて、動きが一定である物体からの反射波を含んでもよい。動きが一定である物体は、周期的な動作を行う。動きが一定である物体としては、例えば、室内に設置された扇風機又はエアコンに含まれるファン及びルーバー（可動羽）が挙げる。動きが一定である物体が存在するレンジビンでは、信号強度ａ１ｉ～ａＪｉを周波数変換して得られる結果において、その動きに応じた周波数に大きな周波数成分（成分のピーク）が表れる。動きが一定である物体に由来する周波数成分は、生体Ｈに由来する周波数成分よりも大きな値になる場合もある。
【0136】
信号強度算出部４２は、第１～第Ｎのフレーム周期それぞれにおいて第１～第Ｊの受信信号データを取得した後に、所定の抑圧処理を実行してもよい。この抑圧処理は、係数算出処理と、係数補正処理と、データ生成処理とを含む。係数算出処理では、信号強度算出部４２が、第１～第Ｎのフレーム周期それぞれにおいて、第１～第Ｋのレンジビンについて上記係数ｖｉを算出する。係数補正処理では、信号強度算出部４２が、第１～第Ｋのレンジビンに含まれるレンジビンごとに、第１～第Ｎのフレーム周期それぞれで算出される上記係数ｖｉの変化の程度を求める。
【0137】
そして、係数補正処理では、信号強度算出部４２が、第１～第Ｎのフレーム周期のうちの第ｎのフレーム周期について、レンジビンごとに、係数ｖｉの変化の程度が小さくなるにつれて値が小さくなるように、係数ｖｉの変化の程度に基づいて係数ｃｉ（補正後の係数）を算出する。係数補正処理では、係数ｖｉの変化の程度が相対的に小さいレンジビンでは、係数ｃｉが小さい値となり、係数ｖｉの変化の程度が相対的に大きいレンジビンでは、係数ｃｉが大きい値となるように、係数ｃｉが算出される。信号強度算出部４２は、係数ｖｉの変化の程度が小さくなるにつれて値が小さくなるように、第ｎのフレーム周期における係数ｖｉの値を補正することで、係数ｃｉ（補正後の係数）を算出してもよい。
【0138】
データ生成処理では、信号強度算出部４２が、第ｎのフレーム周期における第ｊの受信信号データ（例えば、第１の受信信号データ）について、第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと係数ｃｉとを乗算して、第ｎの距離変動データを算出する。動きが一定である物体が存在するレンジビン、及び、固定物が存在するレンジビンでは、複数のフレーム周期で周波数変換された結果（周波数成分のピークが表れる周波数）の変化の程度が小さい。固定物が存在するレンジビンでは、周波数がゼロ付近に成分のピークが表れ、動きが一定である物体が存在するレンジビンでは、その動きに応じた周波数に成分のピークが表れる。これに対して、生体Ｈが存在するレンジビンでは、生体Ｈの動きが一定ではないために、複数のフレーム周期において、周波数変換された結果の変化の程度が相対的に大きい。そのため、上述の係数補正処理を実行することで、生体Ｈが存在するレンジビンでの係数に比べて、動きが一定である物体又は固定物が存在するレンジビンでの係数が小さくなる。
【0139】
図２５は、第１～第Ｎの距離変動データを含む時系列データを生成する手順の一例を示すフローチャートである。最初に、制御部２４（信号強度算出部４２及び記憶部４３）が、図２３に示されるフローチャートと同様に、ステップＳ１１～Ｓ１４を実行する。これにより、１つのフレーム周期において、第１～第Ｋのレンジビンの全てについて係数ｖｉが算出される。
【0140】
ステップＳ１４の実行後、制御部２４は、初期時刻から上記設定時間が経過したかどうかを判断する（図２５のステップＳ１７）。上記設定時間を経過していないと判断された場合（ステップＳ１７：ＮＯ）、制御部２４が実行する処理は、ステップＳ１１に戻り、上記ステップＳ１１～Ｓ１４の一連の処理が繰り返される。これにより、１つのフレーム周期（第１のフレーム周期）と、そのフレーム周期とは異なる別のフレーム周期（第２のフレーム周期）とを含む第１～第Ｎのフレーム周期それぞれにおいて第１～第Ｊの受信信号データの取得及び記憶が行われ、第１～第Ｎのフレーム周期それぞれにおいて係数ｖｉが算出される。
【0141】
上記設定時間が経過したと判断された場合（ステップＳ１７：ＹＥＳ）、制御部２４が実行する処理は、ステップＳ２１に進む。ステップＳ２１では、制御部２４の信号強度算出部４２が、第１～第Ｋのレンジビンに含まれるレンジビンごとに、第１～第Ｎのフレーム周期それぞれで算出される上記係数ｖｉの変化の程度を求める。信号強度算出部４２は、例えば、レンジビンごとに、１つのフレーム周期と次のフレーム周期との間での係数ｖｉの変化量を算出することを第Ｎのフレーム周期まで繰り返して、その変化量の総計（累計）を変化の程度を示す指標として算出する。
【0142】
そして、信号強度算出部４２は、レンジビンごとに、係数ｖｉの変化の程度が小さくなるにつれて値が小さくなるように、係数ｖｉの変化量の総計に基づいて係数ｃｉ（補正後の係数）を算出する。この場合、信号強度算出部４２は、第１～第Ｋのレンジビンにおいて係数ｃ１～ｃＫを算出する。信号強度算出部４２は、第１～第Ｋのレンジビンにおいて、上記変化量の総計を正規化して得られる値を係数ｃｉ（係数ｃ１～ｃＫ）として算出してもよい。又は、信号強度算出部４２は、上記変化量の総計を所定の閾値と比較し、上記変化量の総計が閾値よりも小さくなるレンジビンでの係数ｃｉをゼロに設定し、上記変化量の総計が閾値以上となるレンジビンでの係数ｃｉを１に設定してもよい。
【0143】
次に、信号強度算出部４２は、第ｎのフレーム周期における第ｊの受信信号データについて、レンジビンごとに、信号強度ａｊｉに係数ｃｉを乗算することで第ｎの距離変動データを算出する処理を、第１～第Ｎのフレーム周期それぞれにおいて実行する。そして、記憶部４３が、第１～第Ｎのフレーム周期における第１～第Ｎの距離変動データを記憶する（図２５のステップＳ２２）。なお、このフローチャートの例では、Ｎ個のフレーム周期での係数ｖｉの変化の程度が求められて係数ｃｉが算出されているが、第１～第Ｎのフレーム周期のうちの任意に選択される複数（２以上）のフレーム周期において係数ｖｉの変化の程度が求められて、係数ｃｉが算出されてもよい。例えば、時間的に連続した複数のフレーム周期での係数ｖｉの変化の程度が求められてもよく、時間的に連続していない複数のフレーム周期での係数ｖｉの変化の程度が求められてもよい。
【0144】
ステップＳ２２の実行後、制御部２４は、図１０に示されるステップＳ４～Ｓ８の一連の処理を実行することで、係数補正処理を含む抑圧処理が施されたうえで取得された時系列データに基づいて、心拍間隔データ等の生体データを生成する。
【0145】
例１．本開示の１つの例に係る生体データ取得装置は、第１～第Ｎ（Ｎは２以上の自然数）の距離変動データが時系列順に並べられた時系列データを記憶するように構成された記憶部（４３）であって、第１～第Ｎの距離変動データはそれぞれ、互いに異なる時刻（時刻ｔ１～ｔＮ）において、広帯域又は超広帯域の無線電波が生体（Ｈ）において反射された反射波に基づいて得られる、距離に対する信号強度の変化を示すデータである、記憶部（４３）と、第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち生体（Ｈ）の所定の検知部位に対応する対応強度を取得する強度取得処理を第１～第Ｎの距離変動データに対して行うことにより、対応強度が時系列に沿って並べられた時間変動データを生成するように構成された第１の生成部（４５）と、時間変動データに基づいて、生体（Ｈ）の検知部位における生体データを生成するように構成された第２の生成部（４６）と、第１～第Ｊ（Ｊは２以上の自然数）の受信信号データに対して抑圧処理を行うことで、前記第ｎの距離変動データを生成するように構成された信号強度算出部（４２）とを備える。第１～第Ｊの受信信号データは、フレーム周期内の時系列順に並ぶ異なる複数の時刻において、生体（Ｈ）が存在する検知領域（Ｒ）にて無線電波が反射された反射波に基づいて得られる。第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、第ｊの受信信号データは、検知領域（Ｒ）を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含む。抑圧処理は、第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、第ｊの受信信号データについて、前記第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと前記係数ｖｉとを乗算することで、前記第ｎの距離変動データを算出することと、を含む。この場合、所定期間において反射波から得られる一定数の距離変動データが時系列データとして記憶部に蓄積されている。そして、時系列データから、検知部位に対応する信号強度が時系列に沿って並べられた時間変動データが得られる。時間変動データに含まれる信号強度は、生体の検知部位に対応しているので、検知対象である生体が動いていても、検知部位に対応した信号強度に基づいて生体データが得られる。また、上記抑圧処理が行われることで、第ｎの距離変動データにおいて、検知対象外の物体からの反射波に基づく信号強度と、生体からの反射波に基づく信号強度との差を大きくすることができる。その結果、非接触にて生体データを精度良く検知することが可能となる。
【0146】
例２．本開示の１つの例に係る生体データ取得装置は、第１～第Ｎ（Ｎは２以上の自然数）の距離変動データが時系列順に並べられた時系列データを記憶するように構成された記憶部（４３）であって、第１～第Ｎの距離変動データはそれぞれ、時系列順に並ぶ異なる時刻（時刻ｔ１～ｔＮ）において、広帯域又は超広帯域の無線電波が生体（Ｈ）において反射された反射波に基づいて得られる、距離に対する信号強度の変化を示すデータである、記憶部（４３）と、第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち生体（Ｈ）の所定の検知部位に対応する対応強度を取得する強度取得処理を第１～第Ｎの距離変動データに対して行うことにより、対応強度が時系列に沿って並べられた時間変動データを生成するように構成された第１の生成部（４５）と、時間変動データに基づいて、生体（Ｈ）の検知部位における生体データを生成するように構成された第２の生成部（４６）と、第１～第Ｊの受信信号データに対して抑圧処理を行うことで、前記第ｎの距離変動データを生成するように構成された信号強度算出部（４２）とを備える。第１～第Ｊ（Ｊは２以上の自然数）の受信信号データは、フレーム周期内の時系列順に並ぶ異なる複数の時刻において、生体が存在する検知領域にて無線電波が反射した反射波に基づいて得られる。第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、第ｊの受信信号データは、検知領域を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含む。抑圧処理は、第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することことと、第１～第Ｊの受信信号データのうちの複数の受信信号データを積分して積分信号データを生成することであって、積分信号データは、第１～第Ｋのレンジビンにおける信号強度ｂ１～ｂＫを含むことと、積分信号データについて第１～第Ｋのレンジビンのそれぞれにおいて信号強度ｂｉと係数ｖｉとを乗算して、第ｎの距離変動データを算出することとを含む。この場合、例１の生体データ取得装置と同様に、検知対象である生体が動いていても、検知部位に対応した信号強度に基づいて生体データが得られる。また、積分信号データの算出を含む上記抑圧処理が行われることで、第ｎの距離変動データにおいて、検知対象外の物体からの反射波に基づく信号強度と、生体からの反射波に基づく信号強度との差をより大きくすることができる。従って、生体データのより高精度な検知を行うことが可能となる。
【0147】
例３．例１又は例２の生体データ取得装置において、信号強度算出部（４２）は、上記抑圧処理を実行することの前にフレーム設定処理を実行してもよい。フレーム設定処理は、生体（Ｈ）の動作の速度を示す動作情報を取得することと、動作情報が示す速度に基づいてフレーム周期の長さを設定することとを含んでもよい。この場合、生体の動作が速くなると、フレーム周期を短くすることができる。フレーム周期が短いと第１～第Ｎの距離変動データに生体の動きが反映されやすい。従って、生体データのより高精度な検知を行うことが可能となる。
【0148】
例４．本開示の１つの例に係る生体データ取得装置は、第１～第Ｎ（Ｎは２以上の自然数）の距離変動データが時系列順に並べられた時系列データを記憶するように構成された記憶部（４３）であって、第１～第Ｎの距離変動データはそれぞれ、時系列順に並ぶ異なる時刻（時刻ｔ１～ｔＮ）において、広帯域又は超広帯域の無線電波が生体（Ｈ）において反射された反射波に基づいて得られる、距離に対する信号強度の変化を示すデータである、記憶部（４３）と、第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち生体（Ｈ）の所定の検知部位に対応する対応強度を取得する強度取得処理を第１～第Ｎの距離変動データに対して行うことにより、対応強度が時系列に沿って並べられた時間変動データを生成するように構成された第１の生成部（４５）と、時間変動データに基づいて、生体（Ｈ）の検知部位における生体データを生成するように構成された第２の生成部（４６）と、第１のフレーム周期と、当該第１のフレーム周期とは異なる第２のフレーム周期とのそれぞれにおいて第１～第Ｊ（Ｊは２以上の自然数）の受信信号データを取得するように構成された信号強度算出部（４２）とを備える。第１～第Ｊの受信信号データは、１つのフレーム周期内の時系列順に並ぶ異なる複数の時刻において、生体（Ｈ）が存在する検知領域（Ｒ）にて無線電波が反射された反射波に基づいて得られる。第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、第ｊの受信信号データは、検知領域（Ｒ）を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含む。信号強度算出部（４２）は、所定の抑圧処理を行うことで、第１のフレーム周期における第ｎの距離変動データを生成するように構成されている。抑圧処理は、第１のフレーム周期及び第２のフレーム周期のそれぞれにおいて、第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは、１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、第１～第Ｋのレンジビンの全てについて、第１のフレーム周期と第２のフレーム周期との間での係数ｖｉの変化の程度が小さくなるにつれて値が小さくなるように、係数ｖｉの変化の程度に基づいて係数ｃｉを算出することと、第１のフレーム周期の第ｊの受信信号データについて、第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと係数ｃｉとを乗算することで、第ｎの距離変動データを算出することとを含む。この場合、例１の生体データ取得装置と同様に、検知対象である生体が動いていても、検知部位に対応した信号強度に基づいて生体データが得られる。また、係数ｃｉの算出を含む上記抑圧処理が行われることで、第ｎの距離変動データにおいて、検知対象外であり動作しない物体又は動きが一定である物体からの反射波に基づく信号強度と、生体からの反射波に基づく信号強度との差を大きくすることができる。従って、生体データのより高精度な検知を行うことが可能となる。
【0149】
例５．本開示の１つの例に係る生体データ取得方法は、第１～第Ｎ（Ｎは２以上の自然数）の距離変動データが時系列順に並べられた時系列データを生成することであって、第１～第Ｎの距離変動データはそれぞれ、互いに異なる時刻において、広帯域又は超広帯域の無線電波が生体において反射された反射波に基づいて得られる、距離に対する信号強度の変化を示すデータである、時系列データを生成することと、第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち生体（Ｈ）の所定の検知部位に対応する対応強度を取得する強度取得処理を第１～第Ｎの距離変動データに対して行うことにより、対応強度が時系列に沿って並べられた時間変動データを生成することと、時間変動データに基づいて、生体（Ｈ）の検知部位における生体データを生成することと、第１～第Ｊ（Ｊは２以上の自然数）の受信信号データに対して抑圧処理を行うことで、第ｎの距離変動データを生成することとを含む。第１～第Ｊの受信信号データは、フレーム周期内の時系列順に並ぶ異なる複数の時刻において、生体（Ｈ）が存在する検知領域（Ｒ）にて無線電波が反射された反射波に基づいて得られる。第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、第ｊの受信信号データは、検知領域（Ｒ）を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含む。抑圧処理は、第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、第ｊの受信信号データについて、第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと係数ｖｉとを乗算することで、第ｎの距離変動データを算出することと、を含む。この場合、例１の生体データ取得装置と同様の作用効果が奏される。
【0150】
例６．本開示の１つの例に係る生体データ取得方法は、第１～第Ｎ（Ｎは２以上の自然数）の距離変動データが時系列順に並べられた時系列データを生成することであって、第１～第Ｎの距離変動データはそれぞれ、時系列順に並ぶ異なる時刻において、広帯域又は超広帯域の無線電波が生体において反射された反射波に基づいて得られる、距離に対する信号強度の変化を示すデータである、時系列データを生成することと、第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち生体（Ｈ）の所定の検知部位に対応する対応強度を取得する強度取得処理を第１～第Ｎの距離変動データに対して行うことにより、対応強度が時系列に沿って並べられた時間変動データを生成することと、時間変動データに基づいて、生体（Ｈ）の検知部位における生体データを生成することと、第１～第Ｊ（Ｊは２以上の自然数）の受信信号データに対して抑圧処理を行うことで、第ｎの距離変動データを生成することとを含む。第１～第Ｊの受信信号データは、フレーム周期内の時系列順に並ぶ異なる複数の時刻において、生体（Ｈ）が存在する検知領域（Ｒ）にて無線電波が反射した反射波に基づいて得られる。第１～第Ｊの受信信号データのうちの第ｊの受信信号データは、距離に対する信号強度の変化を示すデータであり、第ｊ（ｊは１～Ｊの自然数）の受信信号データは、検知領域（Ｒ）を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含む。抑圧処理は、第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、第１～第Ｊの受信信号データのうちの複数の受信信号データを積分して積分信号データを生成することであって、積分信号データは、第１～第Ｋのレンジビンにおける信号強度ｂ１～ｂＫを含むことと、積分信号データについて第１～第Ｋのレンジビンのそれぞれにおいて信号強度ｂｉと係数ｖｉとを乗算して、第ｎの距離変動データを算出することとを含む。この場合、例２の生体データ取得装置と同様の作用効果が奏される。
【0151】
例７．本開示の１つの例に係る生体データ取得方法は、第１～第Ｎ（Ｎは２以上の自然数）の距離変動データが時系列順に並べられた時系列データを生成することであって、第１～第Ｎの距離変動データはそれぞれ、時系列順に並ぶ異なる時刻において、広帯域又は超広帯域の無線電波が生体において反射された反射波に基づいて得られる、距離に対する信号強度の変化を示すデータである、時系列データを生成することと、第１～第Ｎの距離変動データのうち第ｎ（ｎは１～Ｎの自然数）の距離変動データにおける信号強度のうち生体（Ｈ）の所定の検知部位に対応する対応強度を取得する強度取得処理を第１～第Ｎの距離変動データに対して行うことにより、対応強度が時系列に沿って並べられた時間変動データを生成することと、時間変動データに基づいて、生体（Ｈ）の検知部位における生体データを生成することと、第１のフレーム周期と、当該第１のフレーム周期とは異なる第２のフレーム周期とのそれぞれにおいて第１～第Ｊ（Ｊは２以上の自然数）の受信信号データを取得することと、所定の抑圧処理を行うことで、第１のフレーム周期における第ｎの距離変動データを生成することとを含む。第１～第Ｊの受信信号データは、１つのフレーム周期内の時系列順に並ぶ異なる複数の時刻において、生体（Ｈ）が存在する検知領域（Ｒ）にて無線電波が反射された反射波に基づいて得られる。第１～第Ｊの受信信号データのうちの第ｊ（ｊは１～Ｊの自然数）の受信信号データは、距離に対する信号強度の変化を示すデータであり、第ｊの受信信号データは、検知領域（Ｒ）を距離方向に分割して得られる第１～第Ｋ（Ｋは２以上の自然数）のレンジビンにおける信号強度ａｊ１～ａｊＫを含む。抑圧処理は、第１のフレーム周期及び第２のフレーム周期のそれぞれにおいて、第１～第Ｋのレンジビンの全てについて、信号強度ａ１ｉ～ａＪｉ（ｉは、１～Ｋの自然数）を周波数変換して得られた結果に基づいて係数ｖｉを算出することと、第１～第Ｋのレンジビンの全てについて、第１のフレーム周期と第２のフレーム周期との間での係数ｖｉの変化の程度が小さくなるにつれて値が小さくなるように、係数ｖｉの変化の程度に基づいて係数ｃｉを算出することと、第１のフレーム周期の第ｊの受信信号データについて、第１～第Ｋのレンジビンのそれぞれにおいて信号強度ａｊｉと係数ｃｉとを乗算することで、第ｎの距離変動データを算出することとを含む。この場合、例４の生体データ取得装置と同様の作用効果が奏される。
【符号の説明】
【0152】
１…生体データ取得システム、１０…送受信装置、１１…発信機、１２…受信機、１３…通信機、２０…生体データ取得装置、２４…制御部、４２…信号強度算出部、４３…記憶部、４４…データ補正部、４５…時間変動データ生成部、４６…生体データ生成部、６０…乗物。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図25】