(57)【特許請求の範囲】
【請求項１】
生体の体表面に電磁波を送信し、上記生体の体表面で反射される反射波を受信信号として受信する信号受送信手段と、
上記受信信号を直交検波して、Ｉ信号とＱ信号を生成する直交検波手段と、
上記生成されたＩ信号とＱ信号とを時系列に順次取得するＩＱ信号取得手段と、
所定の時点において取得したＩ信号とＱ信号のＩＱ平面上の速度ベクトルの一定期間での方向変化に基づいて、上記生体の呼吸の変化を示す呼吸信号を生体情報として生成する信号処理手段とを備え、
上記信号処理手段は、上記取得したＩＱ平面上の速度ベクトルの一定期間での方向変化を、所定の時点でのＩＱ平面上の速度ベクトルと、一定期間過去において取得したＩ信号とＱ信号のＩＱ平面上の速度ベクトルのなす角の余弦として算出することを特徴とする生体状態取得装置。
【請求項２】
上記信号処理手段は、上記ＩＱ平面上の速度ベクトルの一定期間での方向変化を検出するときに、上記呼吸信号に対して互いに異なる第１と第２のしきい値を設定して、上記呼吸信号が上記第１と第２の両方のしきい値を同じ方向で通過したことを検出することを特徴とする請求項１記載の生体状態取得装置。
【請求項３】
上記生成された呼吸信号に対して所定のヒステリシス特性を用いてフィルタするヒステリシスフィルタをさらに備えたことを特徴とする請求項１又は２記載の生体状態取得装置。
【請求項４】
生体の体表面に電磁波を送信し、上記生体の体表面で反射される反射波を受信信号として受信する信号受送信手段と、
上記受信信号を直交検波して、Ｉ信号とＱ信号を生成する直交検波手段と、
上記生成されたＩ信号とＱ信号とを時系列に順次取得するＩＱ信号取得手段と、
所定の時点において取得したＩ信号とＱ信号のＩＱ平面上の速度ベクトルの一定期間での方向変化に基づいて、上記生体の呼吸の変化を示す呼吸信号を生体情報として生成する信号処理手段とを備え、
上記信号処理手段は、上記ＩＱ平面上の速度ベクトルの一定期間での方向変化を検出するときに、上記呼吸信号に対して互いに異なる第１と第２のしきい値を設定して、上記呼吸信号が上記第１と第２の両方のしきい値を同じ方向で通過したことを検出することを特徴とする生体状態取得装置。
【請求項５】
生体の体表面に電磁波を送信し、上記生体の体表面で反射される反射波を受信信号として受信する信号受送信手段と、
上記受信信号を直交検波して、Ｉ信号とＱ信号を生成する直交検波手段と、
上記生成されたＩ信号とＱ信号とを時系列に順次取得するＩＱ信号取得手段と、
所定の時点において取得したＩ信号とＱ信号のＩＱ平面上の速度ベクトルの一定期間での方向変化に基づいて、上記生体の呼吸の変化を示す呼吸信号を生体情報として生成する信号処理手段と、
上記生成された呼吸信号に対して所定のヒステリシス特性を用いてフィルタするヒステリシスフィルタとを備えたことを特徴とする生体状態取得装置。
【請求項６】
上記生成されたＩ信号とＱ信号とに基づいて、速度ベクトルの変化方向を判定することで上記生体の呼気と吸気とを判定する呼吸判定手段をさらに備えたことを特徴とする請求項１〜５のうちのいずれか１つに記載の生体状態取得装置。
【請求項７】
コンピュータを、
請求項１〜６のうちのいずれか１つに記載のＩＱ信号取得手段と、
請求項１〜６のうちのいずれか１つに記載の信号処理手段と、
請求項６記載の呼吸判定手段として機能させるための生体状態取得プログラム。
【請求項８】
請求項１〜６のうちのいずれか１つに記載の生体状態取得装置と、
上記生体状態取得装置で取得された生体状態に基づいて機器本体の運転を制御する制御手段を備えたことを特徴とする機器。
【請求項９】
請求項１〜６のうちのいずれか１つに記載の生体状態取得装置と、
空間を空調する空調手段と、
上記生体状態取得装置で取得された生体情報に基づいて空調手段を制御する制御手段を備えたことを特徴とする空調機器。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば人間などの生体の状態（以下、生体状態という。）を取得する生体状態取得装置、生体状態取得プログラム、生体状態取得装置を備えた機器及び空調機器に関する。
【背景技術】
【０００２】
従来、非接触でかつ非拘束で睡眠中の人の生体信号を取得し、状態判定を行う技術が提案されている。例えば特許文献１に開示された生体状態取得装置は、睡眠中の生体の体表面に電磁波を送信し反射波を受信する送受信部と、反射波の直交検波を行ってＩ信号とＱ信号とからなるＩＱ信号を生成する検波部と、ＩＱ平面上での信号の軌跡から生体の状態を判定する信号処理部とを備えて構成される。ここで、信号処理部はＩＱ平面上においてＩＱ信号の速度ベクトルのノルムを算出することで呼吸信号を得て、この増減を、しきい値を設定してカウントすることで、呼吸頻度を算出する。
【０００３】
このような方法は、受信信号をアナログ処理する際に生じる信号の劣化（ＩＱ平面上での原点ずれなど）がある場合でも、周波数解析を用いずに簡易に呼吸頻度を検出することを特徴としている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１１−１５８８７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、ＩＱ平面上で得られるＩＱ信号の速度ベクトルのノルムは、電磁波の送受信部と反射面の距離及び反射強度などに影響を受ける。このため、この従来例では、速度ベクトルのノルム増減をカウントするためにしきい値を調整する必要があった。
【０００６】
本発明の目的は以上の問題点を解決し、従来技術に比較して簡単な方法で呼吸信号を得ることができる生体状態取得装置を提供することにある。
【課題を解決するための手段】
【０００７】
本発明の一態様に係る生体状態取得装置は、
生体の体表面に電磁波を送信し、上記生体の体表面で反射される反射波を受信信号として受信する信号受送信手段と、
上記受信信号を直交検波して、Ｉ信号とＱ信号を生成する直交検波手段と、
上記生成されたＩ信号とＱ信号とを時系列に順次取得するＩＱ信号取得手段と、
所定の時点において取得したＩ信号とＱ信号のＩＱ平面上の速度ベクトルの一定期間での方向変化に基づいて、上記生体の呼吸の変化を示す呼吸信号を生体情報として生成する信号処理手段とを備え、
上記信号処理手段は、上記取得したＩＱ平面上の速度ベクトルの一定期間での方向変化を、所定の時点でのＩＱ平面上の速度ベクトルと、一定期間過去において取得したＩ信号とＱ信号のＩＱ平面上の速度ベクトルのなす角の余弦として算出することを特徴とする。
【発明の効果】
【０００８】
本発明によれば、呼吸信号として固定された範囲内で変動する時系列信号を得ることができる。これによりしきい値の設定が容易になり呼吸頻度を精度よく算出することができる。
【図面の簡単な説明】
【０００９】
【図１】本発明の実施の形態１に係る生体状態取得装置の構成を示すブロック図である。
【図２】図１の生体状態取得装置において取得されたＩＱ信号の挙動を示すグラフである。
【図３】図１の生体状態取得装置においてＩＱ信号のＩＱ平面上での速度ベクトルの算出方法を説明するＩＱ信号のグラフである。
【図４】図１の生体状態取得装置において速度ベクトルの変化を算出する方法を示すＩＱ信号のグラフである。
【図５】図１の生体状態取得装置において呼吸曲線の算出及び呼吸状態の判定例を示す呼吸信号の波形図である。
【図６】図１の生体状態取得装置において呼吸頻度の利用例を示す、呼吸頻度に対する保湿量を示すグラフである。
【図７】本発明の実施の形態２に係る生体状態取得装置の構成を示すブロック図である。
【図８】図７の生体状態取得装置における呼吸曲線の補正例を示す呼吸信号の波形図である。
【図９】図７の生体状態取得装置において用いるヒステリシスフィルタ１２の入出力電圧特性を示すグラフである。
【図１０】図７の生体状態取得装置における、ヒステリシスフィルタ１２によるフィルタ処理後の呼吸曲線の一例を示す呼吸信号の波形図である。
【図１１】本発明の実施の形態３に係る生体状態取得装置の構成を示すブロック図である。
【図１２】図１１の生体状態取得装置における、速度ベクトルの変化方向の判定方法を示すＩＱ信号のグラフである。
【図１３】図１１の生体状態取得装置において呼吸信号、並びに呼気及び吸気判定曲線を用いた呼気／吸気の判定方法を示すグラフである。
【発明を実施するための形態】
【００１０】
以下、本発明に係る実施形態について図面を参照して説明する。なお、以下の各実施形態において、同様の構成要素については同一の符号を付している。
【００１１】
実施の形態１．
図１は本発明の実施の形態１に係る生体状態取得装置の構成を示すブロック図である。図１において、本実施の形態にかかる生体状態取得装置は例えばヒト又は動物である生体２０の生体状態を取得する装置であって、本実施の形態では、生体状態として呼吸頻度を検出し、呼吸頻度は単位時間当たりの呼吸回数を表すものとする。
【００１２】
図１において、生体状態取得装置は、ドップラレーダセンサ１と、直交検波回路２と、バンドパスフィルタ３と、ＡＤ変換器４と、演算装置５と、表示装置１１とを備えて構成される。ここで、演算装置５は、ＩＱ信号取得回路６と、信号処理回路７と、呼吸判定回路１０とを備える。信号処理回路７は、速度ベクトル変化検出回路８と、呼吸信号生成回路９とを備える。
【００１３】
ドップラレーダセンサ１は信号送受信手段として設けられ、ドップラレーダセンサ１のアンテナは生体２０の体表面（好ましくは、胸表面）に向けられており、生体２０の体表面に向けて電磁波が照射され、当該体表面で反射された電磁波である反射波を受信信号として受信できるように設置されている。直交検波回路２は受信信号を直交検波することでＩ信号とＱ信号からなるＩＱ信号を生成してバンドパスフィルタ３に出力する。バンドパスフィルタ３は、予め呼吸によるドップラシフト成分以外を除去するように設定された周波数通過特性を有し、フィルタ後のＩＱ信号をＡＤ変換器４に出力する。ＡＤ変換器４は入力されるアナログＩＱ信号をＡＤ変換値のデジタルＩＱ信号に変換して演算装置５に出力する。
【００１４】
演算装置５のＩＱ信号取得回路６は信号バッファメモリを有し、デジタルＩＱ信号を時系列で順次、所定の周期で離散的に受信して所定期間のデジタルＩＱ信号を当該信号バッファメモリに格納して信号処理回路７における信号処理に用いる。速度ベクトル変化検出回路８は、任意の時点におけるＩＱ信号のＩＱ平面上の速度ベクトルと、一定期間過去において取得したＩＱ信号のＩＱ平面上の速度ベクトルのなす角の余弦を算出する。また、呼吸信号生成回路９は算出された余弦を呼吸信号として呼吸判定回路１０に出力する。呼吸判定回路１０は呼吸信号の増減から呼吸を判定して呼吸の頻度（以下、呼吸頻度という。）を算出し、算出結果を表示装置１１に表示し、もしくはプリンタなどの出力装置に出力する。
【００１５】
次いで、速度ベクトル変化検出回路８で行われる速度ベクトルの算出について以下説明する。
【００１６】
図２は図１の生体状態取得装置においてヒトに対して取得されたＩＱ信号の挙動を示すグラフである。図２に示すように、ＡＤ変換器４から出力されたＩ信号及びＱ信号の組であるＩＱ信号の座標をＩＱ平面上にプロットすると、生体２０の一例である例えばヒトが呼吸するときの体表面の動きに応じて図示するような軌跡を描く。ここで、ＩＱ信号の座標は、体表面がドップラレーダセンサ１に近づく場合は反時計回り、体表面がドップラレーダセンサ１ら遠ざかる場合は時計回りに軌跡を描く。なお、ドップラレーダセンサ１は、生体の胸郭の膨張収縮を精度良く捉えられる位置に設置されることが望ましい。
【００１７】
呼吸時の体表面の動作は、吸い込み時は、吸い込み開始後、徐々に速度が上昇してピークに達する。そして、呼吸の端点（呼気及び吸気の切り替わり時）に向かって今度は速度が下降し、呼吸の端点で略ゼロとなる。また、呼吸の端点前後では速度ベクトルが逆方向になる。
【００１８】
図３は図１の生体状態取得装置においてＩＱ信号のＩＱ平面上での速度ベクトルの算出方法を説明するＩＱ信号のグラフである。図３において、時刻ｔにおけるＩＱ平面上での座標がＰ（ｔ）とし、時刻ｔから一定時間（Δｔ）過去である時刻ｔ−Δｔにおける座標がＰ（ｔ−Δｔ）であった場合、時刻ｔにおける速度ベクトルは次式として算出する。
【００１９】
Ｖ（ｔ）＝Ｐ（ｔ）−Ｐ（ｔ−Δｔ）
【００２０】
また、時刻ｔにおける速度ベクトルがＶ（ｔ）であった場合、Ｖ（ｔ）と一定時間（Δｔｖ）過去である時刻ｔ−Δｔｖにおける速度ベクトルＶ（ｔ−Δｔｖ）がなす角の余弦を算出することで、速度ベクトルＶ（ｔ）のＶ（ｔ−Δｔｖ）からの方向の変化を求めることができる。
【００２１】
図４は図１の生体状態取得装置において速度ベクトルの変化を算出する方法を示すＩＱ信号のグラフである。すなわち、図４はヒトの呼吸動作時のＩＱ信号のＩＱ平面上での速度ベクトルの推移を模擬したものである。図４において、Ｖ（ｔ１）、Ｖ（ｔ１−Δｔｖ）、Ｖ（ｔ２−Δｔｖ）は、吸気時の速度ベクトルを表しており、Ｖ（ｔ２）、Ｖ（ｔ３−Δｔｖ）、Ｖ（ｔ３）は、呼気時の速度ベクトルを表している。なお、○は吸気から呼気への切り替えを示し、□は呼気から吸気への切り替えを示す。
【００２２】
図４において、吸気時の速度ベクトルＶ（ｔ１）が、Ｖ（ｔ１−Δｔｖ）となす角はπ／２より小さく、これらの余弦は＋１に近い値となる。呼気中の速度ベクトルＶ（ｔ３）が、Ｖ（ｔ３−Δｔｖ）となす角はπ／２よりも小さく、これらの余弦も＋１に近い値となる。また、吸気から呼気に切り替わる時点では、速度ベクトルＶ（ｔ２）は、Ｖ（ｔ２−Δｔｖ）とは逆方向になるため、その余弦は−１に近い値となる。
【００２３】
速度ベクトルＶ（ｔ）と一定時間過去である時刻ｔ−Δｔｖにおける速度ベクトルＶ（ｔ−Δｔｖ）がなす角の余弦を算出し、この時系列変化を呼吸信号と定義すると、呼吸動作を繰り返す時の呼吸信号は図５のように推移する。
【００２４】
図５は図１の生体状態取得装置において呼吸曲線の算出及び呼吸状態の判定例を示す呼吸信号の波形図である。なお、本実施の形態では呼吸信号生成回路９は速度ベクトル変化検出回路８が出力した値を、そのまま呼吸信号として出力している。なお、ＡＤ変換器４におけるＩＱ信号のサンプリング周期は、例えば５０サンプル／秒であり、差分時間Δｔ、Δｔｖは両方とも例えば１秒程度であるが、必要な計測精度や、要求されるコストなどを勘案して、適宜設定することができる。以上の方法による算出を行えば、一回の呼吸動作が長く、ＩＱ平面上で数周回転するような場合であっても、変動が±１の間に収まる呼吸信号を得ることができる。
【００２５】
なお、得られた呼吸信号から単位時間当たりの呼吸頻度を算出する場合は、しきい値の調整は不要で、しきい値をゼロとして交点をカウントすればよい。具体的には一回の呼吸につき４回ゼロを通過するので、単位時間当たりにカウントした交点の数を４で割った数を呼吸頻度とすればよい。
【００２６】
以上の実施の形態によれば、呼吸信号として固定された範囲内で変動する時系列信号を得ることができる。これによりしきい値の設定が容易になり呼吸頻度を精度よく算出することができる。また、呼吸信号の変動範囲を既知の限定した範囲内に収めることができ、呼気及び吸気の切り替えを検出するためのしきい値調整が不要となる。すなわち、しきい値を角度とした場合は必ずπ／２を通過するのでこれをしきい値とすればよい。
【００２７】
さらに、方向変化を任意の時点でのＩＱ平面上の速度ベクトルと、一定期間過去において取得したＩＱ平面上の速度ベクトルのなす角の余弦として算出しているので、呼吸信号の変動範囲を±１の範囲内に収めることができ、呼気及び吸気の切り替えを検出するためのしきい値調整が不要となる。
【００２８】
なお、演算装置５から出力された呼吸頻度は別の演算装置に入力され、機器制御などに用いてもよい。
【００２９】
図６は図１の生体状態取得装置において呼吸頻度の利用例を示す、呼吸頻度に対する保湿量を示すグラフである。例えば、睡眠中の呼吸を補助するための加湿器において、図６に示すように加湿空気放出量（保湿量）を呼吸頻度に応じて逐次調整する制御が考えられる。この制御により、吸頻度が高く鼻腔の粘膜が乾燥しやすい場合に重点的に水分を供給することで、空気に対する過剰な水分供給を抑制することができる。図６では、呼吸頻度と保湿量を比例関係で表わしているが、この限りではなく、呼吸頻度が増加した場合に保湿量も増加する関係であればよい。従って、結果として、加湿器の設置された部屋における結露発生を防ぐことができ、また加湿器に組み込まれる給水タンクの容量を低減できる。
【００３０】
実施の形態２．
図７は本発明の実施の形態２に係る生体状態取得装置の構成を示すブロック図である。図７において、実施の形態２に係る生体状態取得装置は、図１の生体状態取得装置に比較して、以下の点が異なる。
（１）呼吸判定回路１０は、後述するように２個のしきい値Ａ，Ｂ（図８）を設定して呼吸の判定を行う。
（２）呼吸信号生成回路９と呼吸判定回路１０との間に、ヒステリシスフィルタ１２を挿入して設けた。
【００３１】
図８は図７の生体状態取得装置における呼吸曲線の補正例を示す呼吸信号の波形図である。また、図９は図７の生体状態取得装置において用いるヒステリシスフィルタ１２の入出力電圧特性を示すグラフである。さらに、図１０は図７の生体状態取得装置における、ヒステリシスフィルタ１２によるフィルタ処理後の呼吸曲線の一例を示す呼吸信号の波形図である。
【００３２】
図８に示すように、呼吸信号と周波数の近いノイズが重畳した場合、これをバンドパスフィルタ３などで取り除くことは難しく、実施の形態１で示した方法で呼吸頻度を算出すると、実際よりも多く呼吸をカウントしてしまう。このような場合、図８に示すように２個のしきい値Ａ，Ｂ（例えば呼吸信号の最大値と０との間のプラスのしきい値Ａ、及び最小値と０との間のマイナスのしきい値Ｂ）を設定して、呼吸信号が両方のしきい値Ａ，Ｂを同じ方向に通過したことを検出することで呼吸頻度を算出してもよい。
【００３３】
図８にはしきい値Ａとしきい値Ｂがあり、呼吸信号がこれらを通過した点を矢印で示してある。しきい値Ａを下から上に向かって通過した場合はＵＡ、上から下に向かって通過した場合はＤＡ、しきい値ＢについてはそれぞれＵＢ、ＤＢとしてそれぞれの矢印近くに示した。ＵＢ、ＵＡと連続して検出した場合、ＤＡ、ＤＢと連続して検出した場合のみをカウントすると、×で示した部分はカウントされない。
【００３４】
このような処理は、例えば図９のようなヒステリシス特性を持ったヒステリシスフィルタ１２を用いることで実現する。図９のヒステリシス特性を有するヒステリシスフィルタ１２を適用することによって、図８で示した呼吸信号は図１０に示すような波形になり、しきい値Ａ及びしきい値Ｂと比較することで、一呼吸を従来技術に比較して高精度で検出することができる。
【００３５】
以上の本実施の形態では、例えば図８に両側矢印線で示した区間を一呼吸として検出することで、呼吸頻度を従来技術に比較して高精度で算出することができる。また、２個のしきい値Ａ，Ｂを設定して、呼吸信号が両方のしきい値を同じ方向に通過したことを検出することで、呼吸信号と近い周波数をもつノイズが乗じた場合であっても、正確な呼吸数カウントが可能となる。
【００３６】
実施の形態３．
図１１は本発明の実施の形態３に係る生体状態取得装置の構成を示すブロック図である。図１１において、実施の形態３に係る生体状態取得装置は、図１の呼吸判定回路１０に代えて、さらにＩＱ信号のＩＱ平面上の速度ベクトルの変化方向を判定する呼吸判定回路１０Ａを備えたことを特徴としている。
【００３７】
呼気と吸気の判別を行う場合には、ＩＱ信号が呼吸端から離れている時点での速度ベクトルの変化方向（回転角度）を算出することが有効である。速度ベクトルの変化方向は、任意の時点での速度ベクトルと過去の速度ベクトルの外積を算出することで実現できる。
【００３８】
図１２は図１１の生体状態取得装置における、速度ベクトルの変化方向の判定方法を示すＩＱ信号のグラフである。
【００３９】
図１２において吸気時の速度ベクトルはＶ（ｔ１）と表わされているが、速度ベクトルＶ（ｔ１−Δｔｖ）から比較して反時計まわりに回転する（回転角度をθ１とする）。逆に呼気時の速度ベクトルはＶ（ｔ３）と表わされているが、速度ベクトルＶ（ｔ３−Δｔｖ）から比較して時計まわりに回転する（回転角度をθ３）。例えば、外積Ｖ（ｔ１）×Ｖ（ｔ１−Δｔｖ）は次式の３次元ベクトルとなる。
【００４０】
【数１】
0006465827.tif 000002
【００４１】
ここで、呼吸の端点から離れている場合に、式（１）の｜Ｖ（ｔ１）｜｜Ｖ（ｔ１−Δｔ）｜ｓｉｎθ１は正の値となる。また、呼気時の外積Ｖ（ｔ３）×Ｖ（ｔ３−Δｔｖ）も同様に計算すると、呼吸の端点から離れている場合は｜Ｖ（ｔ３）｜｜Ｖ（ｔ３−Δｔ）｜ｓｉｎθ３は負の値となる。これらを呼気及び吸気判定曲線とする。ここで、呼吸の端点から離れていることは、実施の形態１及び２で算出した、呼吸曲線に対して所定のしきい値を設定し、呼吸の端点が上記しきい値よりも大きいことにより判定できる。
【００４２】
図１３は図１１の生体状態取得装置において呼吸信号、並びに呼気及び吸気判定曲線を用いた呼気／吸気の判定方法を示すグラフである。図１３では、図５で示した呼吸信号と、当該呼吸信号に対してしきい値Ｃを設定し、呼吸信号がしきい値Ｃよりも大きい場合に、上述した呼気及び吸気判定曲線を同じ時間軸で表した。図１３のグラフ下には判定結果を示している。この例では、左から呼気、吸気、呼気、吸気と順次判定されている。
【００４３】
以上説明したように本実施の形態によれば、ＩＱ信号の速度ベクトルの変化方向を判定することで、呼気又は吸気を従来技術に比較して高精度で判定できる。
【００４４】
変形例．
以上の実施形態においては、演算装置５内の各回路６〜１０を例えばハードウエアで構成しているが、本発明はこれに限らず、演算装置を１つのデジタル計算機などのコンピュータで構成し、各演算をコンピュータにより実行されるソフトウェアプログラムで構成してもよい。なお、ＩＱ信号取得回路６では、ＩＱ信号を一時的に格納するメモリを備える。
【００４５】
以上の実施の形態１〜３の演算装置５から出力された呼吸信号と呼気及び吸気判定結果は、別の制御装置に入力されて、機器制御などに用いることができる。例えば、睡眠中の呼吸を補助するための加湿器において、加湿空気の放出を吸気のタイミングに合わせて行うことが可能である。結果として、空気に対する過剰な水分供給を抑制し、加湿器の設置された部屋における結露を防ぐことができ、また加湿器に組み込まれる給水タンクの容量を低減できる。
【００４６】
なお、生体状態取得装置を備えた機器としては、生体状態取得装置で取得された生体状態に基づいて機器本体の運転を制御する制御手段を備え、例えば、目覚まし時計や空気調和機の他、例えば照明、芳香機能など５感を刺激する機能を有する機器、もしくは、ＴＶや音楽プレーヤーなどのＡＶ機器、湯たんぽ、加湿器、除湿機、空気清浄機などの空調機器であって、空調機器は、生体状態取得装置と、空間を空調する空調手段と、上記生体状態取得装置で取得された生体情報に基づいて空調手段を制御する制御手段を備えた空調機器であってもよい。
【産業上の利用可能性】
【００４７】
以上詳述したように、本発明によれば、呼吸信号として固定された範囲内で変動する時系列信号を得ることができる。これによりしきい値の設定が容易になり呼吸頻度を精度よく算出することができる。
【符号の説明】
【００４８】
１  ドップラレーダセンサ、２  直交検波回路、３  バンドパスフィルタ、４  ＡＤ変換器、５  演算装置、６  ＩＱ信号取得回路、７  信号処理回路、８  速度ベクトル変化検出回路、９  呼吸信号生成回路、１０，１０Ａ  呼吸判定回路、１１  表示装置、１２  ヒステリシスフィルタ、２０  生体。
【図１】
0006465827.tif 000003
【図２】
0006465827.tif 000004
【図３】
0006465827.tif 000005
【図４】
0006465827.tif 000006
【図５】
0006465827.tif 000007
【図６】
0006465827.tif 000008
【図７】
0006465827.tif 000009
【図８】
0006465827.tif 000010
【図９】
0006465827.tif 000011
【図１０】
0006465827.tif 000012
【図１１】
0006465827.tif 000013
【図１２】
0006465827.tif 000014
【図１３】
0006465827.tif 000015