(57)【要約】      （修正有）【課題】被験者の呼吸数の推定値を、高い信頼性をもって出力することのできる生体状態モニタリングシステムを提供する。【解決手段】ベッド上の被験者の荷重を検出する荷重検出器と、荷重検出器の検出値の時間的変動に基づいて、被験者に体動が生じているか否かを判定する体動判定部と、荷重検出器の検出値の時間的変動に基づいて、被験者の呼吸数の推定値を逐次求めて出力する呼吸数推定部と、被験者の体動の持続時間に基づいて、体動がトゥイッチであるか否かを判定するトゥイッチ判定部とを備える。トゥイッチ判定部が体動がトゥイッチではないと判定した場合には、呼吸数推定部は、呼吸数の推定値の出力を停止し、又は体動の直前の、体動判定部が被験者に体動が生じていないと判定した期間において最後に求めた第１推定値を出力する。体動がトゥイッチであると判定した場合には、呼吸数推定部は、新たに求める呼吸数の推定値を順次出力する。【選択図】図９
2019187941.tif 000002
【特許請求の範囲】
【請求項１】
ベッド上の被験者の生体状態をモニタする生体状態モニタリングシステムであって、
前記ベッド上の被験者の荷重を検出する少なくとも１つの荷重検出器と、
前記荷重検出器の検出値の時間的変動に基づいて、被験者に体動が生じているか否かを判定する体動判定部と、
前記荷重検出器の検出値の時間的変動に基づいて、被験者の呼吸数の推定値を逐次求めて出力する呼吸数推定部と、
被験者の体動の持続時間に基づいて、該体動がトゥイッチであるか否かを判定するトゥイッチ判定部とを備え、
前記トゥイッチ判定部が前記体動がトゥイッチではないと判定した場合には、前記呼吸数推定部は、前記呼吸数の推定値の出力を停止し、又は前記体動の直前の、前記体動判定部が被験者に体動が生じていないと判定した期間において最後に求めた第１推定値を出力し、前記トゥイッチ判定部が前記体動がトゥイッチであると判定した場合には、前記呼吸数推定部は、新たに求める呼吸数の推定値を順次出力する生体状態モニタリングシステム。
【請求項２】
前記トゥイッチ判定部が前記体動がトゥイッチではないと判定した場合には、前記呼吸数推定部は第１推定値を出力する請求項１に記載の生体状態モニタリングシステム。
【請求項３】
前記トゥイッチ判定部は、前記体動判定部の判定結果が体動有りである期間の長さと所定の閾値との比較に基づいて、被験者の体動がトゥイッチであるか否かを判定する請求項１又は２に記載の生体状態モニタリングシステム。
【請求項４】
ベッド上の被験者の生体状態をモニタする生体状態モニタリングシステムであって、
前記ベッド上の被験者の荷重を検出する少なくとも１つの荷重検出器と、
前記荷重検出器の検出値の時間的変動に基づいて、被験者に体動が生じているか否かを判定する体動判定部と、
前記荷重検出器の検出値の時間的変動に基づいて、被験者の呼吸数の推定値を逐次求めて出力する呼吸数推定部とを備え、
前記体動判定部が被験者に体動が生じていると判定する期間中は、前記呼吸数推定部は、前記呼吸数の推定値の出力を停止し、又は前記期間の直前の、前記体動判定部が被験者に体動が生じていないと判定した期間において最後に求めた第１推定値を出力する生体状態モニタリングシステム。
【請求項５】
前記体動判定部が被験者に体動が生じていると判定する期間中は、前記呼吸数推定部は第１推定値を出力する請求項４に記載の生体状態モニタリングシステム。
【請求項６】
更に、前記呼吸数推定部が出力した前記被験者の呼吸数の推定値を表示する表示部を備える請求項１〜５のいずれか一項に記載の生体状態モニタリングシステム。
【請求項７】
ベッドと、
請求項１〜６のいずれか一項に記載の生体状態モニタリングシステムとを備えるベッドシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、荷重検出器の検出値に基づいてベッド上の被験者の生体状態をモニタする生体状態モニタリングシステムに関する。
【背景技術】
【０００２】
医療や介護の分野において、荷重検出器を介してベッド上の被験者の荷重を検出し、検出した荷重に基づいて被験者の生体状態をモニタすることが提案されている。具体的には例えば、検出した荷重に基づいて被験者の呼吸数を推定することが提案されている。
【０００３】
特許文献１は、ベッドの脚部の下に配置した荷重センサの計測値を周波数解析して、被験者の呼吸数を検出することを開示している。特許文献２は、ベッドの脚部の下に荷重検出器を配置してベッド上の被検生体の重心の移動を求め、この重心の移動に基づいて被検生体の呼吸運動と心拍運動とを求めることを開示している。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第４８８３３８０号
【特許文献２】特公昭６１−２４０１０号
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明は、被験者の呼吸数の推定値を、高い信頼性をもって出力することのできる生体状態モニタリングシステムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明の第１の態様に従えば、
ベッド上の被験者の生体状態をモニタする生体状態モニタリングシステムであって、
前記ベッド上の被験者の荷重を検出する少なくとも１つの荷重検出器と、
前記荷重検出器の検出値の時間的変動に基づいて、被験者に体動が生じているか否かを判定する体動判定部と、
前記荷重検出器の検出値の時間的変動に基づいて、被験者の呼吸数の推定値を逐次求めて出力する呼吸数推定部と、
被験者の体動の持続時間に基づいて、該体動がトゥイッチであるか否かを判定するトゥイッチ判定部とを備え、
前記トゥイッチ判定部が前記体動がトゥイッチではないと判定した場合には、前記呼吸数推定部は、前記呼吸数の推定値の出力を停止し、又は前記体動の直前の、前記体動判定部が被験者に体動が生じていないと判定した期間において最後に求めた第１推定値を出力し、前記トゥイッチ判定部が前記体動がトゥイッチであると判定した場合には、前記呼吸数推定部は、新たに求める呼吸数の推定値を順次出力する生体状態モニタリングシステムが提供される。
【０００７】
第１の態様の生体状態モニタリングシステムにおいて、前記トゥイッチ判定部が前記体動がトゥイッチではないと判定した場合には、前記呼吸数推定部は第１推定値を出力してもよい。
【０００８】
第１の態様の生体状態モニタリングシステムにおいて、前記トゥイッチ判定部は、前記体動判定部の判定結果が体動有りである期間の長さと所定の閾値との比較に基づいて、被験者の体動がトゥイッチであるか否かを判定してもよい。
【０００９】
本発明の第２の態様に従えば、
ベッド上の被験者の生体状態をモニタする生体状態モニタリングシステムであって、
前記ベッド上の被験者の荷重を検出する少なくとも１つの荷重検出器と、
前記荷重検出器の検出値の時間的変動に基づいて、被験者に体動が生じているか否かを判定する体動判定部と、
前記荷重検出器の検出値の時間的変動に基づいて、被験者の呼吸数の推定値を逐次求めて出力する呼吸数推定部とを備え、
前記体動判定部が被験者に体動が生じていると判定する期間中は、前記呼吸数推定部は、前記呼吸数の推定値の出力を停止し、又は前記期間の直前の、前記体動判定部が被験者に体動が生じていないと判定した期間において最後に求めた第１推定値を出力する生体状態モニタリングシステムが提供される。
【００１０】
第２の態様の生体状態モニタリングシステムにおいて、前記体動判定部が被験者に体動が生じていると判定する期間中は、前記呼吸数推定部は第１推定値を出力してもよい。
【００１１】
第１の態様又は第２の態様の生体状態モニタリングシステムは、更に、前記呼吸数推定部が出力した前記被験者の呼吸数の推定値を表示する表示部を備えてもよい。
【００１２】
本発明の第３の態様に従えば、
ベッドと、
第１の態様又は第２の態様の生体状態モニタリングシステムとを備えるベッドシステムが提供される。
【発明の効果】
【００１３】
本発明の生体状態モニタリングシステムは、被験者の呼吸数の推定値を、高い信頼性をもって出力することができる。
【図面の簡単な説明】
【００１４】
【図１】図１は、本発明の第１実施形態に係る生体状態モニタリングシステムの構成を示すブロック図である。
【図２】図２は、荷重検出器のベッドに対する配置を示す説明図である。
【図３】図３は、本発明の第１実施形態の生体状態モニタリングシステムを用いた呼吸数の推定方法を示すフローチャートである。
【図４】図４は、荷重検出器により検出された荷重値の変動の様子を、被験者が呼吸のみを行っている安静期間と、被験者が体動を行っている体動期間の両方について示す概略的なグラフである。
【図５】図５（ａ）は、被験者の重心が、被験者の呼吸に応じて被験者の体軸方向に振動する様子を概念的に示す説明図である。図５（ｂ）は、被験者の呼吸に応じた被験者の重心の振動に基づいて描画される呼吸波形の一例を示すグラフである。
【図６】図６は、体動判定部における体動判定の具体的方法の一例を説明するための説明図である。
【図７】図７は、被験者の呼吸数を呼吸波形に基づいて推定する方法の一例を説明するための説明図である。
【図８】図８は、呼吸数推定部により行われる呼吸数推定工程の具体例を説明するための説明図である。
【図９】図９は、本発明の第２実施形態に係る生体状態モニタリングシステムの構成を示すブロック図である。
【図１０】図１０は、本発明の第２実施形態の生体状態モニタリングシステムを用いた呼吸数の推定方法を示すフローチャートである。
【図１１】図１１は、第２実施形態の呼吸数推定工程中の、被験者の体動がトゥイッチか否かに基づいて出力する推定値の内容を異ならせる工程を示すフローチャートである。
【図１２】図１２（ａ）、図１２（ｂ）は、呼吸数推定部がトゥイッチ判定部の判定結果を参照して出力する呼吸数の推定値の調整を行う工程について説明するための説明図である。
【図１３】図１３は、変形例に係るベッドシステムの全体構成を示すブロック図である。
【発明を実施するための形態】
【００１５】
＜第１実施形態＞
本発明の第１実施形態の生体状態モニタリングシステム１００（図１）について、これをベッドＢＤ（図２）と共に使用して、ベッドＢＤ上の被験者Ｓの呼吸数を推定（決定）する場合を例として説明する。
【００１６】
図１に示す通り、本実施形態の生体状態モニタリングシステム１００は、荷重検出部１、制御部３、記憶部４を主に有する。荷重検出部１と制御部３とは、Ａ／Ｄ変換部２を介して接続されている。制御部３には更に表示部５、報知部６、入力部７が接続されている。
【００１７】
荷重検出部１は、４つの荷重検出器１１、１２、１３、１４を備える。荷重検出器１１、１２、１３、１４のそれぞれは、例えばビーム形のロードセルを用いて荷重を検出する荷重検出器である。このような荷重検出器は例えば、特許第４８２９０２０号や特許第４００２９０５号に記載されている。荷重検出器１１、１２、１３、１４はそれぞれ、配線又は無線によりＡ／Ｄ変換部２に接続されている。
【００１８】
図２に示す通り、荷重検出部１の４つの荷重検出器１１〜１４は、被験者Ｓが使用するベッドＢＤの四隅の脚ＢＬ１、ＢＬ２、ＢＬ３、ＢＬ４の下端部に取り付けられたキャスターＣ１、Ｃ２、Ｃ３、Ｃ４の下にそれぞれ配置される。
【００１９】
Ａ／Ｄ変換部２は、荷重検出部１からのアナログ信号をデジタル信号に変換するＡ／Ｄ変換器を備え、荷重検出部１と制御部３にそれぞれ配線又は無線で接続されている。
【００２０】
制御部３は、専用又は汎用のコンピュータであり、内部に標準偏差算出部３１、呼吸波形描画部（呼吸波形算出部、呼吸波形取得部）３２、体動判定部３３、呼吸数推定部（呼吸数決定部、呼吸数算出部）３４が構築されている。
【００２１】
記憶部４は、生体状態モニタリングシステム１００において使用されるデータを記憶する記憶装置であり、例えばハードディスク（磁気ディスク）を用いることができる。
【００２２】
表示部５は、制御部３から出力される情報を生体状態モニタリングシステム１００の使用者に表示する液晶モニター等のモニターである。
【００２３】
報知部６は、制御部３からの情報に基づいて所定の報知を聴覚的に行う装置、例えばスピーカを備える。
【００２４】
入力部７は、制御部３に対して所定の入力を行うためのインターフェイスであり、キーボード及びマウスにし得る。
【００２５】
このような生体状態モニタリングシステム１００を使用して、ベッド上の被験者の呼吸数を推定（決定）する動作について説明する。
【００２６】
生体状態モニタリングシステム１００を使用した被験者の呼吸数の推定は、図３のフローチャートに示す通り、被験者Ｓの荷重を検出する荷重検出工程Ｓ１１と、検出した荷重の変動の程度を示す標準偏差を算出する標準偏差算出工程Ｓ１２と、検出した荷重に基づいて被験者の呼吸波形を描画する呼吸波形描画工程Ｓ１３と、標準偏差算出工程Ｓ１２で求めた標準偏差及び呼吸波形描画工程Ｓ１３で描いた呼吸波形の振幅を用いて被験者の体動判定を行う体動判定工程Ｓ１４と、呼吸波形描画工程Ｓ１３で描いた呼吸波形に基づいて被験者の呼吸数を推定して出力する呼吸数推定工程Ｓ１５と、出力された呼吸数を表示する表示工程Ｓ１６とを含む。
【００２７】
［荷重検出工程］
荷重検出工程Ｓ１１では、荷重検出器１１、１２、１３、１４を用いてベッドＢＤ上の被験者Ｓの荷重を検出する。ベッドＢＤ上の被験者Ｓの荷重は、ベッドＢＤの四隅の脚ＢＬ１〜ＢＬ４の下に配置された荷重検出器１１〜１４に分散して付与され、これらによって分散して検出される。
【００２８】
荷重検出器１１〜１４はそれぞれ、荷重（荷重変化）を検出してアナログ信号としてＡ／Ｄ変換部２に出力する。Ａ／Ｄ変換部２は、サンプリング周期を例えば５ミリ秒として、アナログ信号をデジタル信号に変換し、デジタル信号（以下「荷重信号」）として制御部３に出力する。以下では、荷重検出器１１、１２、１３、１４から出力されたアナログ信号をＡ／Ｄ変換部２においてデジタル変換して得られる荷重信号を、それぞれ荷重信号ｓ１、ｓ２、ｓ３、ｓ４と呼ぶ。
【００２９】
［標準偏差算出工程］
標準偏差算出工程Ｓ１２では、標準偏差算出部３１が、荷重信号ｓ１、ｓ２、ｓ３、ｓ４の各々について、その所定のサンプリング期間（一例として５秒間）に含まれるサンプリング値の標準偏差（移動標準偏差）σ１、σ２、σ３、σ４を算出する。算出は、例えば過去５秒間をサンプリング期間として常時行われ得る。
【００３０】
標準偏差は、サンプリング値のばらつきの大きさを表わすため、図４に示すように、ベッドＢＤ上の被験者Ｓが安静にしており、荷重信号ｓ１〜ｓ４の変動の量が小さい期間ＰＥＲ１においては標準偏差σ１〜σ４も小さくなる。一方で被験者Ｓが身体を動かしており（被験者Ｓに体動が生じており）、荷重信号ｓ１〜ｓ４の変動の量が大きい期間ＰＥＲ２においては標準偏差σ１〜σ４も大きくなる。
【００３１】
したがって、被験者Ｓに体動が生じている期間においては、被験者Ｓに体動が生じていない期間に比べて標準偏差σ１〜σ４の値が大きくなる。
【００３２】
本明細書及び本発明において「体動」とは、被験者の頭部、胴部（体幹）、四肢の移動を意味する。呼吸や心拍等に伴う臓器、血管等の移動は体動には含まれない。体動は、一例として、被験者Ｓの胴部（体幹）の移動を伴う大きな体動と、被験者の四肢や頭部の移動のみを伴う小さな体動とに分類し得る。大きな体動の一例は、寝返りや起き上がり等であり、小さな体動の一例は、睡眠中の手足や頭部の移動等である。被験者に心拍や呼吸、体動が生じた場合には、これに応じて荷重検出器１１〜１４からの荷重信号ｓ１〜ｓ４が変動する。その変動の量は、被験者Ｓの心拍に応じた変動の量、被験者Ｓの呼吸に応じた変動の量、被験者Ｓの小さな体動に応じた変動の量、被験者Ｓの大きな体動に応じた変動の量の順に大きくなる。
【００３３】
なお、本明細書及び本発明に記載する被験者の体動判定において、荷重信号ｓ１〜ｓ４の、被験者Ｓの心拍に応じた変動の量は無視できるほどに小さい。そのため、本明細書及び本発明において、「被験者が呼吸のみを行う」や、荷重値、荷重信号が「呼吸のみに応じて変動する」とは、被験者が体動を行っていないことや、荷重値、荷重信号が体動に応じた変動を示していないことを意味し、被験者に心拍が生じていないことや、荷重値、荷重信号が心拍に応じた変動を含んでいないことを意味するものではない。
【００３４】
［呼吸波形描画工程］
呼吸波形描画工程Ｓ１３では、呼吸波形描画部（呼吸波形算出部、呼吸波形取得部）３２が、荷重信号ｓ１〜ｓ４に基づいて被験者Ｓの呼吸波形を描画する。
【００３５】
人間の呼吸は、胸郭及び横隔膜を移動させて、肺を膨張及び収縮させることにより行われる。ここで吸気時、すなわち肺が膨張する時には横隔膜は下方に下がり、内臓も下方に移動する。一方で呼気時、すなわち肺が収縮する時には横隔膜は上方に上がり、内臓も上方に移動する。本件の出願人に付与された特許第６１０５７０３号の明細書に記載されている通り、この内臓移動に伴って重心Ｇはわずかに移動し、その移動方向は背骨の延在方向（体軸方向）にほぼ沿っている。
【００３６】
本発明及び本明細書において「呼吸波形」とは、被験者の呼吸に応じて被験者の体軸方向に振動する被験者の重心の振動の様子を、時間軸に展開して示す波形を意味する。呼吸波形の１周期は、被験者の１回の呼吸（呼気及び吸気）に対応する。呼吸波形の振幅は、被験者の体格や呼吸の深さの影響を受ける。具体的には例えば、被験者が大柄であったり、被験者が深い呼吸を行った場合には振幅は大きくなり、被験者が小柄であったり、被験者が浅い呼吸を行った場合には振幅は小さくなる。
【００３７】
呼吸波形描画部３２は、具体的には次のようにして呼吸波形の描画を行う。
【００３８】
呼吸波形描画部３２は、まず、荷重検出部１からの荷重信号ｓ１〜ｓ４に基づいて、各サンプリング時刻ごとに被験者Ｓの重心Ｇの位置を算出する。被験者Ｓの重心Ｇは、図５（ａ）に示すように、被験者Ｓの呼吸に応じて、被験者Ｓの体軸ＳＡの方向に振動している。
【００３９】
呼吸波形描画部３２は、次いで、体軸ＳＡの方向を縦軸、時間軸を横軸として、各時間における重心Ｇの位置を体軸ＳＡに投影した位置と、重心Ｇの呼吸に応じた振動の振動中心との間の距離を縦軸にプロットすることにより、呼吸波形ＢＷ（図５（ｂ））を描画する。
【００４０】
なお、呼吸波形描画部３２は、必ずしも実際に呼吸波形を描画する必要はなく、呼吸波形を示すデータを求めるのみでもよい。
【００４１】
［体動判定工程］
体動判定工程Ｓ１４では、体動判定部３３が、標準偏差算出工程Ｓ１２で算出した標準偏差σ１〜σ４と、呼吸波形描画工程Ｓ１３で描画した呼吸波形ＢＷの振幅とを用いて、被験者Ｓに体動が生じているか否かの判定を行う。
【００４２】
判定は、具体的には例えば、次のように行う。
【００４３】
体動判定部３３は、まず、呼吸波形描画工程Ｓ１３において描画された呼吸波形ＢＷについてピーク検出を行い、最新のピークｐｎとその１つ前のピークｐｎ−１との間の期間における最小値と、最新のピークｐｎとの間の差を、呼吸波形ＢＷの最新の振幅Ａｎ（図５（ｂ））として算出する。そして、振幅Ａｎと、振幅Ａｎよりも以前に算出された振幅Ａｎ−１、Ａｎ−２、・・・の単純平均値である平均振幅ＡＡＶｎを求める。平均振幅ＡＡＶｎの算出に用いる振幅Ａｎの数は任意であるが、一例として５秒程度のサンプリング期間に得られる数とし得る。
【００４４】
なお、振幅Ａｎ、Ａｎ−１、Ａｎ−２、・・・、及び平均振幅ＡＡＶｎの算出に関わる呼吸波形ＢＷが、安静期間（被験者が体動を生ずることなく呼吸のみを行っている期間）に得られた呼吸波形であることを確認するために（即ち、体動に応じた重心移動に起因する誤差を含んでいないことを確認するために）、標準偏差σ１〜σ４の単純平均値σＡＶ、或いは標準偏差σ１〜σ４のいずれか１つを、閾値σ０と比較してもよい。閾値σ０は単純平均値σＡＶ、或いは標準偏差σ１〜σ４のいずれか１つがこの値未満であれば確実に被験者Ｓに体動は生じていないと判定し得るほど小さな値とする（反対に、単純平均値σＡＶ、或いは標準偏差σ１〜σ４のいずれか１つがこの値以上であっても、被験者Ｓに体動は生じていない場合があり得る）。
【００４５】
次に、体動判定部３３は、下記の式（１）により、規格化標準偏差σｓ１〜σｓ４を求める。
【数１】
2019187941.tif 000003
【００４６】
そして、算出した規格化標準偏差σｓ１〜σｓ４の単純平均値σｓＡＶを求め、単純平均値σｓＡＶと所定の閾値σｓＴＨとの比較に基づいて、被験者Ｓに体動が生じているか否かを判定する。具体的には例えば、単純平均値σｓＡＶが所定の閾値σｓＴＨ以上であれば被験者Ｓに体動が生じていると判定し、単純平均値σｓＡＶが所定の閾値σｓＴＨ未満であれば被験者Ｓに体動が生じていないと判定する。
【００４７】
被験者Ｓに体動が生じていないと判定されている期間については、体動判定部３３は、呼吸波形描画部３２において取得される呼吸波形ＢＷに基づいて、呼吸波形ＢＷの１周期ごとに、最新の振幅Ａｎを算出する。そして、算出した最新の振幅Ａｎの値を用いて新たに算出した平均振幅ＡＡＶｎと、各サンプリング時刻ごとに（一例として５ミリ秒ごとに）算出される標準偏差σ１〜σ４の値とを用いて、（式１）により、逐次、規格化標準偏差σｓ１〜σｓ４を算出し、単純平均値σｓＡＶを算出する。
【００４８】
被験者Ｓの体動の有無の判定は、このようにして逐次算出される単純平均値σｓＡＶを、所定の閾値σｓＴＨと逐次比較することによりなされる。
【００４９】
判定の結果、被験者Ｓに体動が生じたと判定された後は、体動判定部３３は、平均振幅ＡＡＶｎの更新を停止し、その時点で算出されている平均振幅ＡＡＶｎを用いた単純平均値σｓＡＶの算出を続ける。そして算出された単純平均値σｓＡＶと所定の閾値σｓＴＨとの比較に基づいて、被験者Ｓの体動の有無の判定を続ける。これは、被験者Ｓに体動が生じている期間においては、呼吸波形ＢＷの振幅及び周期が体動の影響により大きく変動し、新たな振幅Ａｎ、及び平均振幅ＡＡＶｎを誤差なく算出することが困難となるためである。
【００５０】
次に、再度、被験者Ｓに体動がないと判定した場合には、体動判定部３３は、呼吸波形描画部３２において取得される呼吸波形ＢＷに基づいて、最新の振幅Ａｎの値を算出し、これを用いて新たに平均振幅ＡＡＶｎを算出する。そして、最新の平均振幅ＡＡＶｎを用いた単純平均値σｓＡＶの算出を行い、算出された単純平均値σｓＡＶと所定の閾値σｓＴＨとの比較に基づいて、被験者Ｓの体動の有無の判定を続ける。
【００５１】
この工程の具体例を、図６を用いて説明する。
【００５２】
図６に示すように、時刻Ｔ１において呼吸波形ＢＷのピークＰ１、最新の振幅Ａ１が得られた場合、体動判定部３３は、時刻Ｔ１以降の期間については、振幅Ａ１とその直前に得られた振幅Ａ０、Ａ−１・・・（いずれも不図示）の単純平均値である平均振幅ＡＡＶ１を用いて、逐次規格化標準偏差σｓ１〜σｓ４、及び単純平均値σｓＡＶ（ここでは、σｓＡＶ（ＡＡＶ１）と記載する）を算出し、単純平均値σｓＡＶ（ＡＡＶ１）と閾値σｓＴＨとの比較に基づき、被験者Ｓの体動の有無を判定する。
【００５３】
次に、時刻Ｔ２において呼吸波形ＢＷのピークＰ２、最新の振幅Ａ２が得られた場合、体動判定部３３は、時刻Ｔ２以降の期間については、振幅Ａ２を用いて算出した平均振幅ＡＡＶ２を用いて、逐次規格化標準偏差σｓ１〜σｓ４、及び単純平均値σｓＡＶ（ここでは、σｓＡＶ（ＡＡＶ２）と記載する）を算出し、単純平均値σｓＡＶ（ＡＡＶ２）と閾値σｓＴＨとの比較に基づき、被験者Ｓの体動の有無を判定する。
【００５４】
時刻Ｔ２１において、被験者Ｓに体動が生じた場合、体動判定部３３は、単純平均値σｓＡＶ（ＡＡＶ２）が閾値σｓＴＨ以上となることに基づいて、被験者Ｓに体動が生じたと判定する。そしてその後、単純平均値σｓＡＶ（ＡＡＶ２）が閾値σｓＴＨ以上である期間については、単純平均値σｓＡＶ（ＡＡＶ２）を用いた体動判定を継続する。
【００５５】
時刻Ｔ２２において、被験者Ｓの体動が収まった場合、体動判定部３３は、単純平均値σｓＡＶ（ＡＡＶ２）が閾値σｓＴＨ未満となることに基づいて、被験者Ｓに体動は生じていないと判定する。その後、再び呼吸波形ＢＷの最新の振幅Ａｎが得られ、最新の平均振幅ＡＡＶｎが算出されるまで、単純平均値σｓＡＶ（ＡＡＶ２）を用いた体動判定を継続する。
【００５６】
体動判定部３３が、呼吸波形の振幅Ａｎ（平均振幅ＡＡＶｎ）を用いて標準偏差σ１〜σ４の値の規格化を行うのは次の理由による。
【００５７】
上述した通り、標準偏差σ１〜σ４の値は、一般に被験者Ｓに体動が生じている期間に大きくなる。そのため、標準偏差σ１〜σ４の値を所定の閾値と比較することにより被験者Ｓに体動が生じているか否かの判定を行うことが考えられる。
【００５８】
しかしながら、本発明の発明者の知見によれば、標準偏差σ１〜σ４の値は、体格の大きい被験者においては、呼吸に応じた臓器の移動量も大きいため、体動を生じることなく呼吸のみを行っている期間においても比較的大きな値となり得る。また、被験者の体格が大きくない場合においても、例えば当該被験者が深呼吸を行った際には、比較的大きな値となり得る。標準偏差σ１〜σ４の値を所定の閾値と比較することにより体動判定を行えば、このような場合に、被験者に体動が生じていないにもかかわらず、被験者に体動が生じたとの誤判定がなされ得る。
【００５９】
一方で、呼吸波形に着目すると、呼吸波形の振幅は、上記の通り、被験者の体格や呼吸の深さの影響を受け、被験者が大柄であったり、被験者が深い呼吸を行った場合には振幅は大きくなり、被験者が小柄であったり、被験者が浅い呼吸を行った場合には振幅は小さくなる。
【００６０】
即ち、標準偏差σ１〜σ４の値を呼吸波形ＢＷの平均振幅ＡＡＶｎで除して、規格化を行うことにより、被験者Ｓの体格や呼吸の深さが標準偏差σ１〜σ４の値に及ぼす影響を軽減（補償、補正）することができる。そして、このような規格化により得られた規格化標準偏差σｓ１〜σｓ４を用いた体動判定を行うことにより、体動判定の精度を高めることができる。
【００６１】
［呼吸数推定工程］
呼吸数推定工程Ｓ１５においては、呼吸数推定部（呼吸数決定部）３４が、呼吸波形描画工程Ｓ１３において描画した被験者Ｓの呼吸波形ＢＷに基づいて被験者Ｓの呼吸数を推定（決定）し、推定（決定）した呼吸数を表示部５に出力する。
【００６２】
呼吸数推定部３４は、具体的には、次の方法により、呼吸波形ＢＷに基づく呼吸数の推定を行う。
【００６３】
呼吸数推定部３４は、呼吸波形描画工程Ｓ１３において描かれた呼吸波形ＢＷに対して、所定の周期でピーク検出を行う。そして新たなピークが検出された場合には、当該ピークが検出された時刻と、そこから所定期間（一例として２５秒間。以下、適宜「推定用サンプリング期間」と呼ぶ）だけ遡った時刻との間に含まれる複数のピークについて、隣接する２つのピークの間の距離（以下、ピーク間距離ｔと呼ぶ）をそれぞれ求める。
【００６４】
図７に、２５秒間の推定用サンプリング期間に含まれる呼吸波形ＢＷの一例、この呼吸波形ＢＷが有するピークｐ１〜ｐ６、及びピークｐ１〜ｐ６の間のピーク間距離ｔ１〜ｔ５を示す。
【００６５】
次に、呼吸数推定部３４は、次の式２により、複数のピーク間距離ｔの平均値ｔＡＶを求める。
【数２】
2019187941.tif 000004
ここで、Ｎは推定用サンプリング期間に含まれる呼吸波形ＢＷが有するピークの数である。この平均値ｔＡＶは、推定用サンプリング期間における被験者Ｓの呼吸の周期の平均値にほぼ相当する。
【００６６】
その後、呼吸数推定部３４は、次の式３により、被験者Ｓの１分間の呼吸数の推定値Ｒを算出する。そして算出した推定値Ｒを、推定した呼吸数として表示部５に出力する。
【数３】
2019187941.tif 000005
【００６７】
また、本実施形態の呼吸数推定部３４は、被験者Ｓの体動の有無に応じて、出力する呼吸数の値を調整する。以下、この調整について説明する。
【００６８】
呼吸数推定部３４は、体動判定部３３が被験者Ｓに体動が生じていないと判定している期間（安静期間）、及び体動判定部３３が被験者Ｓに体動が生じていると判定している期間（体動期間）のいずれにおいても、呼吸波形ＢＷのピークが検出される度に、上記の方法に従い、式２、式３を用いて呼吸数の推定値Ｒを逐次算出する。
【００６９】
そして、安静期間（後述する、体動期間から安静期間に変移した後の所定期間を除く）においては、算出した推定値Ｒを、推定した呼吸数として逐次表示部５に出力する。
【００７０】
一方、呼吸数推定部３４は、体動期間においては、上記の方法に従い式２、式３を用いて逐次算出する呼吸数の推定値Ｒに代えて、安静期間から体動期間に変移する直前に算出した推定値（第１推定値）Ｒ、即ちその直前の安静期間の最後に算出した推定値Ｒを、推定した呼吸数として表示部５に出力する。
【００７１】
また、呼吸数推定部３４は、安静期間から体動期間に変移し、その後、再び安静期間に変移した後も、安静期間への変移時から所定期間が経過するまでは、上記の方法に従い式２、式３を用いて逐次算出する呼吸数の推定値Ｒに代えて、直前の体動期間に出力していた値、即ち安静期間から体動期間に変移する直前に算出した推定値Ｒを、推定した呼吸数として表示部５に出力する。
【００７２】
その後、安静期間の呼吸波形ＢＷのピークが２つ検出された後に、当該２つのピーク間のピーク間距離ｔを求め、このピーク間距離ｔの値を平均値ｔＡＶの値として用いて、式３より呼吸数の推定値Ｒを算出する。そして体動期間及びその後の所定期間において固定的に出力していた値、即ち安静期間から体動期間に変移する直前に算出した推定値Ｒに代えて、新たに算出した推定値Ｒを、推定した呼吸数として表示部５に出力する。
【００７３】
同様に、呼吸数推定部３４は、呼吸波形ＢＷのピークが３つ、４つと検出された時点で、当該３つ、４つのピークの内の隣り合う２つのピーク間のピーク間距離ｔをそれぞれ求め、これらのピーク間距離ｔの単純平均値を平均値ｔＡＶの値として用いて、式３より呼吸数の推定値Ｒを算出する。そして算出した値を、推定した呼吸数として表示部５に出力する。その後、体動期間から安静期間に変化してから推定用サンプリング期間（ここでは２５秒間）が経過した後は、上記の式２を用いた平均値ｔＡＶの算出を再開し、算出した平均値ｔＡＶを上記の式３に適用して推定値Ｒを算出し、算出した値を、推定した呼吸数として表示部５に出力する。
【００７４】
この工程の具体例を、図８を用いて説明する。
【００７５】
図８に示すように、安静期間に含まれる時刻Ｔ４において呼吸波形ＢＷのピークＰ４が検出された後は、呼吸数推定部３４は、時刻Ｔ４を終端とする２５秒間の推定用サンプリング期間に求められた呼吸波形ＢＷに基づいて、上記の式２、式３より呼吸数の推定値Ｒを算出し、算出した値を推定した呼吸数として出力する。ここでは呼吸数Ｒ４が出力されたものとする。同様に、時刻Ｔ５、Ｔ６において呼吸波形ＢＷのピークＰ５、Ｐ６が検出された後は、呼吸数推定部３４は、時刻Ｔ５、Ｔ６を終端とする推定用サンプリング期間に求められた呼吸波形ＢＷに基づいて、上記の式２、式３より呼吸数の推定値Ｒを算出し、算出した値を推定した呼吸数として出力する。ここでは呼吸数Ｒ５、Ｒ６が出力されたものとする。
【００７６】
次に、時刻Ｔ６１において被験者Ｓに体動が生じた場合、呼吸数推定部３４は、体動判定部３３が被験者Ｓに体動が生じたと判定したことに基づいて、出力する呼吸数の推定値を、上記の式２、式３を用いて新たに算出する推定値Ｒから、直前に算出した推定値（呼吸数Ｒ６に等しい）に切り換える。そして、体動判定部３３が被験者Ｓに体動が生じていると判定している期間、即ち体動期間においては、推定した呼吸数として、呼吸数Ｒ６を出力し続ける。
【００７７】
また、呼吸数推定部３４は、時刻Ｔ６２において被験者Ｓの体動が収まり、体動判定部３３が被験者Ｓに体動は生じていないと判定した後（即ち安静期間に戻った後）も、所定期間の間、上記の式２、式３を用いて新たに算出する推定値Ｒではなく、時刻Ｔ６１の直前に算出した推定値（呼吸数Ｒ６に等しい）を、推定した呼吸数として出力し続ける。
【００７８】
その後、呼吸数推定部３４は、時刻Ｔ８において被験者Ｓの体動が収まった後の２つ目のピークＰ８が検出された後に、被験者Ｓの体動が収まった後の１つ目のピークＰ７とピークＰ８との間のピーク間距離ｔ８を求める。そして、ピーク間距離ｔ８を平均値ｔＡＶの値として用いて、式３より呼吸数の推定値Ｒを算出し、時刻Ｔ６１以降固定的に出力していた値（呼吸数Ｒ６）に代えて、新たに算出した推定値Ｒを推定した呼吸数として出力する。ここでは呼吸数Ｒ８が出力されたものとする。
【００７９】
呼吸数推定部３４は、時刻Ｔ９において被験者Ｓの体動が収まった後の３つ目のピークＰ９が検出された後に、被験者Ｓの体動が収まった後の２つ目のピークＰ８とピークＰ９との間のピーク間距離ｔ９を求める。そして、ピーク間距離ｔ８とピーク間距離ｔ９の単純平均値を平均値ｔＡＶの値として用いて、式３より呼吸数の推定値Ｒを算出し、算出した値を推定した呼吸数として出力する。ここでは呼吸数Ｒ９が出力されたものとする。
【００８０】
同様に、呼吸数推定部３４は、時刻Ｔ１０、Ｔ１１において被験者Ｓの体動が収まった後の４つ目のピークＰ１０、５つ目のピークＰ１１が検出された後に、ピークＰ９とピークＰ１０との間のピーク間距離ｔ１０、ピークＰ１０とピークＰ１１との間のピーク間距離ｔ１１を求める。そして、ピーク間距離ｔ８〜ｔ１１の内、既に得られているものの単純平均値を平均値ｔＡＶの値として用いて、式３より呼吸数の推定値Ｒを算出し、算出した値を推定した呼吸数として出力する。ここでは呼吸数Ｒ１０、Ｒ１１が出力されたものとする。
【００８１】
その後、呼吸数推定部３４は、時刻Ｔ６２から推定用サンプリング期間に相当する２５秒が経過した後は、上記の式２を用いた平均値ｔＡＶの算出を再開し、算出した平均値ｔＡＶを上記の式３に適用して呼吸数の推定値Ｒを算出する。そして、算出した値を、推定した呼吸数として表示部５に出力する。
【００８２】
呼吸数推定部３４が、体動期間、及びその後の安静期間の初期の所定期間において、当該期間に算出する推定値Ｒに代えて、安静期間から体動期間に変化する直前に算出した呼吸数の推定値Ｒを、被験者Ｓの呼吸数の推定値として出力する理由は次の通りである。
【００８３】
図８に示す通り、体動期間においても呼吸波形ＢＷにピークは存在し、これが検出される。しかしながら図８に示される波形から明らかなように、呼吸波形ＢＷは被験者Ｓの体動の影響を受けて乱れており、呼吸波形ＢＷの周期も変化している。そのため、体動期間の呼吸波形ＢＷのピークに基づいて求められたピーク間距離は、呼吸の周期を反映しない値となる。推定用サンプリング期間に体動期間が含まれている場合には、式２で求める平均値ｔＡＶの値が、このような呼吸の周期を反映しないピーク間距離の影響を受けた値となるため、式３を用いて算出される呼吸数の推定値Ｒも、被験者の実際の呼吸数とは乖離した値となっている可能性が高い。
【００８４】
そこで、呼吸数推定部３４は、体動期間が推定用サンプリング期間に含まれてしまう体動期間、及び体動期間の後の安静期間の初期の所定期間においては、当該期間の推定値として、新たに算出する推定値Ｒではなく、安静期間から体動期間に変化する直前に算出した推定値Ｒを出力する。これにより、生体状態モニタリングシステム１００の使用者に提示する呼吸数の推定値の信頼性が維持される。
【００８５】
なお、呼吸数推定部３４は、体動期間中、及び体動期間から安静期間に変移した後の所定期間中は、呼吸波形ＢＷのピーク検出や、ピーク間距離の算出、式２、式３を用いた推定値Ｒの算出を停止してもよい。
【００８６】
本実施形態では、安静期間の呼吸波形ＢＷのピークが２つ検出された時点で、これらのピークの間のピーク間距離を用いた推定値Ｒの算出を行い、算出された値を出力している。この方法によれば、体動期間から安静期間に変移した後、推定用サンプリング期間に相当する時間が経過するのを待つことなく、体動期間における呼吸波形ＢＷの乱れの影響を除いた呼吸数の推定値Ｒの算出及び出力を再開することができる。しかしながら、体動終了後、推定用サンプリング期間に略等しい時間が経過するまで、体動期間に変化する直前に算出された推定値Ｒを出力し続けてもよい。
【００８７】
［表示工程］
表示工程Ｓ１６においては、制御部３が、呼吸数推定工程Ｓ１５において出力された被験者Ｓの呼吸数（推定値）を、表示部５に表示する。また表示工程Ｓ１６では、表示部５を用いた表示に加えて、又はこれに代えて、報知部６を用いた報知を行っても良い。この場合は例えば、被験者Ｓの呼吸数が所定の設定値に至った場合に報知音を発し、生体状態モニタリングシステム１００の使用者である看護師や介護士等に報せる。
【００８８】
本実施形態の生体状態モニタリングシステム１００の効果を以下にまとめる。
【００８９】
本実施形態の生体状態モニタリングシステム１００の備える呼吸数推定部３４は、被験者Ｓに体動があった場合に、体動中及び体動終了後の所定期間の間、体動発生直前に算出した推定値を、当該期間の呼吸数の推定値として出力する。したがって、本実施形態の生体状態モニタリングシステム１００は、被験者Ｓに体動が生じた場合も、体動の影響を抑制して、被験者の呼吸数の推定値を高い信頼性をもって出力することができる。
【００９０】
本実施形態の生体状態モニタリングシステム１００の備える呼吸数推定部３４は、体動期間から安静期間に変移した後、呼吸波形ＢＷのピークが２つ検出された時点で、これらのピークの間のピーク間距離を用いた呼吸数の推定値の算出を再開する。したがって、本実施形態の生体状態モニタリングシステム１００は、体動終了後、より早期に、表示部５に表示される呼吸数が体動直前に推定された値に固定された状態を解除することができる。
【００９１】
本実施形態の生体状態モニタリングシステム１００は、ベッドＢＤの脚ＢＬ１〜ＢＬ４の下に配置した荷重検出器１１〜１４を用いて被験者Ｓの体動の有無を判定している。したがって、被験者Ｓの身体に計測装置を取り付ける必要がなく、被験者Ｓに不快感や違和感を与えることがない。
【００９２】
＜第２実施形態＞
次に、本発明の第２実施形態の生体状態モニタリングシステム２００について説明する。
【００９３】
上記の第１実施形態の生体状態モニタリングシステム１００においては、呼吸数推定部３４は、被験者Ｓの体動により呼吸波形に乱れが生じている体動期間等において、出力する被験者Ｓの呼吸数の推定値を、体動期間の直前に算出された推定値（第１推定値）に固定する。生体状態モニタリングシステム１００は、こうすることにより、呼吸波形の乱れの影響で呼吸数の推定値の信頼性が低下することを防止している。
【００９４】
ここで、レム睡眠中の人間には、トゥイッチ（ｔｗｉｔｃｈ）と呼ばれる筋肉の攣縮運動（ピクッとした動き）が生じることが知られている。トゥイッチも、体動の一種である。本発明の発明者は、トゥイッチが生体状態のモニタリングに与える影響を研究し、被験者にトゥイッチが生じても、トゥイッチにより呼吸波形に大きな乱れが生じることはないとの知見を得た。これは、持続時間（継続時間）の短い体動であるトゥイッチの影響が、呼吸波形を取得するための計算処理やフィルタリング処理において除去又は軽減されるためであると考えられる。一方で、本発明の発明者の知見によれば、体動がトゥイッチである場合でも、体動判定部は体動が生じたとの判定を行う。
【００９５】
第２実施形態の生体状態モニタリングシステム２００は、本発明の発明者の上記知見に基づくものであり、体動判定部が被験者に体動が生じたと判定した場合であっても、当該体動がトゥイッチである場合には、第１実施形態の生体状態モニタリングシステム１００の呼吸数推定部３４が行う推定値の固定を行わないように構成されている。これにより、第２実施形態の生体状態モニタリングシステムは、被験者の呼吸数の推定値をより一層高い信頼性をもって出力することができる。
【００９６】
第２実施形態の生体状態モニタリングシステム２００（図９）は、被験者の体動がトゥイッチであるか否かを判定するトゥイッチ判定部３５を備える点を除いては、第１実施形態の生体状態モニタリングシステム１００と同一の構成を有する。第１、第２実施形態の生体状態モニタリングシステム１００、２００において共通する構成については説明を省略する。トゥイッチ判定部３５の動作については後述する。
【００９７】
生体状態モニタリングシステム２００を使用した被験者の呼吸数の推定は、図１０のフローチャートに示す通り、被験者Ｓの荷重を検出する荷重検出工程Ｓ２１、検出した荷重の変動の程度を示す標準偏差を算出する標準偏差算出工程Ｓ２２、検出した荷重に基づいて被験者の呼吸波形を描画する呼吸波形描画工程Ｓ２３、標準偏差算出工程Ｓ２２で求めた標準偏差及び呼吸波形描画工程Ｓ２３で描いた呼吸波形の振幅を用いて被験者の体動判定を行う体動判定工程Ｓ２４、被験者Ｓの体動がトゥイッチであるか否かを判定するトゥイッチ判定工程Ｓ２５、呼吸波形描画工程Ｓ２３で描いた呼吸波形に基づいて被験者の呼吸数を推定して出力する呼吸数推定工程Ｓ２６、及び出力された呼吸数を表示する表示工程Ｓ２７を含む。
【００９８】
荷重検出工程Ｓ２１、標準偏差算出工程Ｓ２２、呼吸波形描画工程Ｓ２３、体動判定工程Ｓ２４、表示工程Ｓ２７は、それぞれ、第１実施形態の生体状態モニタリングシステム１００が実行する荷重検出工程Ｓ１１、標準偏差算出工程Ｓ１２、呼吸波形描画工程Ｓ１３、体動判定工程Ｓ１４、表示工程Ｓ１６と実質的に同一であるため説明を省略する。
【００９９】
［トゥイッチ判定工程、呼吸数推定工程］
トゥイッチ判定工程Ｓ２５においては、トゥイッチ判定部３５が、体動判定部３３の判定結果に基づいて、被験者Ｓに生じている体動がトゥイッチであるか否かを判定する。判定は、具体的には次の方法により行われる。
【０１００】
トゥイッチ判定部３５は、体動判定部３３から体動判定の結果を受け取り、判定結果が体動無しから体動有りに変化した時点で、時間計測を開始する。そして、体動判定部３３の判定結果が体動無しから体動有りに変化した時点の後、判定結果が再び体動無しに戻るまでの長さ（持続時間）を計測により求め、求めた持続時間を所定の閾値Ｄｔｈと比較する。
【０１０１】
比較の結果、ある体動の持続時間が閾値Ｄｔｈよりも短い場合には、当該体動は、持続時間の短いトゥイッチであると判定する。一方で、ある体動の持続時間が閾値Ｄｔｈよりも長い場合には、当該体動は持続時間の長い体動であり、トゥイッチではないと判定する。
【０１０２】
判定に用いる閾値Ｄｔｈは、トゥイッチの一般的な持続時間に基づいて適宜決定し得る。また、トゥイッチの影響によりベッド等に生じる振動が収束するまでの時間も考慮して閾値Ｄｔｈの決定を行っても良く、システムを実装する上でのその他の要因を考慮して閾値Ｄｔｈの決定を行っても良い。閾値Ｄｔｈは、具体的には例えば、０．５秒〜６秒程度とし得る。
【０１０３】
本明細書及び本発明においては、持続時間が所定時間よりも短い体動、換言すれば、その発生が体動判定に影響を及ぼす期間が所定期間（一例として０．５秒〜６秒）よりも短い体動をトゥイッチとみなす。
【０１０４】
呼吸数推定工程Ｓ２６においては、呼吸数推定部３４が、呼吸波形描画部３２が描画する呼吸波形ＢＷ、体動判定部３３の判定結果、及びトゥイッチ判定部３５の判定結果に基づいて、被験者Ｓの呼吸数の推定値を算出し、出力する。
【０１０５】
呼吸数推定部３４は、呼吸波形描画部３２から呼吸波形ＢＷ、体動判定部３３から体動判定の結果を受け取り、これらを所定時間Ｘだけ遅延させて用いる。このような遅延を行う理由は後述する。
【０１０６】
図１１のフローチャートに示す通り、呼吸数推定部３４は、まず、体動判定部３３の判定結果に基づいて、被験者Ｓに体動が生じているか否かを判定する（図１１。Ｓ２６１）。そして、被験者Ｓに体動が生じていない場合（Ｓ２６１：Ｎｏ）は、式２、式３を用いて逐次新たに算出する推定値Ｒを、推定した呼吸数として順次表示部５に出力する（図１１。Ｓ２６３）。
【０１０７】
呼吸数推定部３４は、被験者Ｓに体動が生じている場合（Ｓ２６１：Ｙｅｓ）には、第１実施形態の呼吸数推定工程Ｓ１５とは異なり、次の処理を行う。
【０１０８】
呼吸数推定部３４は、体動判定部３３の判定結果が体動無しから体動有りに変化した時点で、まず、トゥイッチ判定部３５の判定結果に基づいて、当該体動がトゥイッチであるか否かを判定する（図１１。Ｓ２６２）。
【０１０９】
そして、当該体動がトゥイッチではない場合（Ｓ２６２：Ｎｏ）には、第１実施形態の呼吸数推定工程Ｓ１５において被験者Ｓに体動が生じた場合と同様に、式２、式３を用いて逐次算出する呼吸数の推定値Ｒに代えて、安静期間から体動期間に変移する直前に算出した推定値（第１推定値）Ｒを、推定した呼吸数として表示部５に出力する（図１１。Ｓ２６４）。続いて、体動期間から安静期間に移行した後の処理も、第１実施形態の呼吸数推定工程Ｓ１５と同様である。
【０１１０】
一方、当該体動がトゥイッチである場合（Ｓ２６２：Ｙｅｓ）には、呼吸数推定部３４は、安静期間と同様に、式２、式３を用いて新たに算出される推定値Ｒを、推定した呼吸数として逐次表示部５に出力する（図１１。Ｓ２６３）。すなわち、呼吸数推定部３４は、体動判定部３３が被験者Ｓに体動が生じていると判定した場合であっても、トゥイッチ判定部３５により当該体動がトゥイッチであると判定された場合には、当該体動が生じている期間を体動期間とはみなさず、安静期間と同様の処理を行うように構成されている。
【０１１１】
なお、呼吸数推定部３４は、上記の通り、呼吸波形描画部３２から受け取った呼吸波形ＢＷ、及び体動判定部３３から受け取った体動判定の結果を、所定時間Ｘだけ遅延させて用いている。したがって、所定時間Ｘを閾値Ｄｔｈよりも大きな値とすれば、呼吸数推定部３４における処理上で体動判定部３３の判定結果が体動無しから体動有りに変化する時点においては、トゥイッチ判定部３５は、既に、当該体動がトゥイッチであるか否かの判定を完了していることになる。したがって、呼吸数推定部３４は、上記の通り、体動判定部３３の判定結果が体動無しから体動有りに変化した時点で、当該体動がトゥイッチであるか否かを判定（確認）することができる。
【０１１２】
この工程の具体例を、図１２（ａ）、図１２（ｂ）を用いて説明する。なお、図１２（ａ）の下側のグラフ、及び図１２（ｂ）において、縦軸の「１」は体動有りを、「０」は体動無しを意味する。
【０１１３】
図１２（ｂ）に示すように、体動判定部３３が、時刻Ｔ１３に、被験者Ｓに体動が生じたと判定した場合、トゥイッチ判定部３５は、当該時点から、体動判定部３３の判定結果が再び体動無しに戻る時刻Ｔ１４までの持続時間（継続時間）を計測し、求めた持続時間を閾値Ｄｔｈと比較する。ここでは、閾値Ｄｔｈが３秒、時刻Ｔ１３から時刻Ｔ１４までの持続時間が２秒であり、トゥイッチ判定部３５は、時刻Ｔ１３から時刻Ｔ１４の間に生じた体動は、トゥイッチであると判定したものとする。
【０１１４】
その後、図１２（ａ）の下側のグラフに示すように、時刻Ｔ１３から所定時間Ｘ（呼吸数推定部３４におけるデータ遅延時間。一例として約５秒）だけ経過した後の時刻Ｔ１３＋Ｘに、呼吸数推定部３４は、被験者Ｓに体動が生じたとの体動判定部３３の判定結果に基づき、トゥイッチ判定部３５の判定結果を参照して、その後の処理内容を決定する。
【０１１５】
ここでは、トゥイッチ判定部３５は、時刻Ｔ１３から時刻Ｔ１４の間に生じた体動はトゥイッチであると判定しているため、呼吸数推定部３４は、時刻Ｔ１３＋Ｘの後も、式２、式３を用いて新たに算出される推定値Ｒを、推定した呼吸数として逐次表示部５に出力する。推定値は、時刻Ｔ１３頃の被験者Ｓの状態を表わす値であるが、所定時間Ｘだけ遅延して時刻Ｔ１３＋Ｘ以降に出力される。
【０１１６】
なお、図１２（ａ）の上側のグラフに示すように、被験者Ｓにトゥイッチが生じている期間においても、被験者Ｓの呼吸波形に大きな乱れは生じない。そのため、推定用サンプリング期間にトゥイッチが生じている期間が含まれている場合であっても、式２、式３を用いて新たに算出される推定値Ｒの信頼性は高い。
【０１１７】
次いで、体動判定部３３が、時刻Ｔ１５に、被験者Ｓに体動が生じたと判定した場合、トゥイッチ判定部３５は、再度、当該時点から時間計測を開始する。そして、トゥイッチ判定部３５は、体動判定部３３の判定結果が体動有りである期間が閾値Ｄｔｈ（ここでは３秒）を越えた時点で、当該体動はトゥイッチではないと判定する。ここでは、時刻Ｔ１５＋Ｘよりも前に、体動判定部３３が、時刻Ｔ１５から継続して生じている体動はトゥイッチではないと判定したものとする。
【０１１８】
その後、時刻Ｔ１５から所定時間Ｘだけ経過した後の時刻Ｔ１５＋Ｘに、呼吸数推定部３４は、被験者Ｓに体動が生じたとの体動判定部３３の判定結果に基づき、トゥイッチ判定部３５の判定結果を再度参照して、その後の処理内容を決定する。
【０１１９】
ここでは、トゥイッチ判定部３５は、時刻Ｔ１５から継続して生じている体動はトゥイッチではないと判定しているため、呼吸数推定部３４は、時刻Ｔ１５＋Ｘの後は、時刻Ｔ１５＋Ｘの直前に算出した推定値Ｒを推定した呼吸数として表示部５に出力する。推定値は、時刻Ｔ１５以降の被験者Ｓの状態を表わす値であるが、所定時間Ｘだけ遅延して、時刻Ｔ１５＋Ｘ以降に出力される。
【０１２０】
第２実施形態の生体状態モニタリングシステム２００は、第１実施形態の生体状態モニタリングシステム１００が奏する効果に加えて、次の効果を奏する。
【０１２１】
第２実施形態の生体状態モニタリングシステム２００はトゥイッチ判定部３５を備え、被験者Ｓに体動が生じている場合であっても、当該体動が呼吸数の推定値の算出への影響が小さいトゥイッチである場合には、呼吸数推定部３４に、新たに算出される推定値Ｒの出力を続けさせる。このように、第２実施形態の生体状態モニタリングシステム２００は、体動がトゥイッチであるか否かを判定することにより、過去の安静期間に算出した推定値を代替的に出力する期間を必要最小限とすることができる。そのため、被験者の呼吸数の推定値をより一層高い信頼性をもって出力することができる。
【０１２２】
なお、上記の第２実施形態の生体状態モニタリングシステム２００の呼吸数推定部３４は、呼吸波形描画部３２から受け取った呼吸波形ＢＷ、及び体動判定部３３から受け取った体動判定の結果を、所定時間Ｘだけ遅延させて処理を行っていたが、これには限られない。例えば、システムにおける処理内容によっては、呼吸波形ＢＷの描画が、体動判定部３３における体動判定よりも、所定時間だけ遅れる場合もある。この場合は、呼吸数推定部３４において改めて遅延を行わなくても、トゥイッチ判定部３５における判定を、呼吸数推定部３４における処理よりも前に完了することが出来る。
【０１２３】
［変形例］
上記実施形態の生体状態モニタリングシステム１００、２００において、次の変形態様を採用することもできる。
【０１２４】
上記実施形態の生体状態モニタリングシステム１００、２００の呼吸数推定部３４は、被験者Ｓに体動が生じたと判定された後、体動期間内のみ、出力する呼吸数の推定値を、体動期間の直前に算出された推定値Ｒに維持するように構成されていてもよい。体動期間は、被験者Ｓに体動が生じている限り継続し得るため長期にわたり得る。一方で、体動終了後の所定期間は、長くても推定用サンプリング期間程度であり、上記第１、第２実施形態では２５秒程度である。したがって、この変形例のように、呼吸数推定部３４が出力する推定値の調整を体動期間においてのみ行う構成としても、被験者Ｓの体動の影響を相当程度抑制して、被験者の呼吸数の推定値を高い信頼性をもって出力することができる。
【０１２５】
上記実施形態の生体状態モニタリングシステム１００、２００は、被験者Ｓに体動が生じたと判定された後、体動期間、及び／又は体動期間直後の安静期間初期の所定期間において、被験者Ｓの呼吸数の推定値の出力を停止してもよい。この場合、生体状態モニタリングシステム１００、２００の使用者は、例えば最後に出力された呼吸数の推定値（即ち、体動期間の直前に算出された推定値Ｒ）を記憶することにより、被験者の現在の呼吸数の推定値を推測できる。この変形例の構成によっても、被験者Ｓの体動の影響を抑制して、被験者の呼吸数の推定値を高い信頼性をもって出力することができる。
【０１２６】
上記実施形態の生体状態モニタリングシステム１００、２００において、呼吸数推定部３４は、呼吸波形ＢＷを用いずに被験者Ｓの呼吸数を推定してもよい。
【０１２７】
具体的には例えば、所定のサンプリング期間に対応する荷重信号ｓ１〜ｓ４の少なくとも一つについてフーリエ解析を行い、呼吸に相当する周波数帯域（人間の呼吸は１分間に約１２〜２０回程度行われるため、約０．２Ｈｚ〜約０．３３Ｈｚである）に現れるピーク周波数を特定することにより行う。特定したピーク周波数より、その期間における被験者Ｓの呼吸数を算出（推定）できる。
【０１２８】
この場合も、サンプリング期間に体動期間が含まれる場合には、体動に起因する荷重信号ｓ１〜ｓ４の変動により、算出される被験者Ｓの呼吸数（推定値）の信頼性が低下し得る。したがって、上記実施形態の呼吸数推定部３４と同様の方法により出力する値を調整することで、呼吸数の推定値の信頼性を高めることができる。
【０１２９】
上記実施形態の生体状態モニタリングシステム１００、２００において、体動判定部３３は、標準偏差σ１〜σ４を呼吸波形ＢＷの平均振幅ＡＡＶｎで除することにより標準偏差σ１〜σ４の規格化を行い、被験者Ｓの体格や呼吸の深さが標準偏差σ１〜σ４の値に及ぼす影響を軽減（補償、補正）していたが、呼吸波形ＢＷの振幅を用いた補償（補正）の方法はこれには限られない。具体的には例えば、呼吸波形ＢＷの最新の振幅Ａｎで除することにより標準偏差σ１〜σ４の規格化を行ってもよい。
【０１３０】
なお、呼吸波形ＢＷの振幅を用いて標準偏差σ１〜σ４を補償した後に、補償後の値を閾値と比較する代わりに、呼吸波形ＢＷの振幅を用いて閾値を補償し、標準偏差σ１〜σ４と補償後の閾値とを比較することもできる。具体的には例えば、上記実施形態においては、標準偏差σ１〜σ４を平均振幅ＡＡＶｎで除する代わりに、所定の閾値にＡＡＶｎを乗じることにより、所望の補償（補正）を行うことが出来る。このように、両者は実質的に等価であり、比較に用いる標準偏差と閾値のいずれの値を呼吸波形の振幅により補償（補正）するかは適宜選択し得る。本明細書及び本発明においては、「標準偏差を呼吸波形の振幅により補償する」との文言は、閾値を呼吸波形の振幅により補償することも含むものとする。
【０１３１】
上記実施形態の生体状態モニタリングシステム１００、２００において、体動判定部３３は、規格化標準偏差σｓ１〜σｓ４の単純平均値σｓＡＶを閾値σｓＴＨと比較して体動判定を行っていたが、規格化標準偏差σｓ１〜σｓ４の少なくとも１つと閾値との比較、規格化標準偏差σｓ１〜σｓ４の少なくとも２つ以上の単純平均値と閾値との比較、規格化標準偏差σｓ１〜σｓ４の少なくとも２つ以上の合計値と閾値との比較等により体動判定を行うこともできる。
【０１３２】
また、上記実施形態の生体状態モニタリングシステム１００、２００において、体動判定部３３は、標準偏差に代えて、標準偏差の二乗である分散を用いることもできる。分散を呼吸波形の振幅で規格化する場合には、分散を呼吸波形の振幅の二乗で除してもよい。したがって本明細書及び本発明において、標準偏差は分散も含むものとする。
【０１３３】
体動判定部３３は、体動判定に用いる閾値にヒステリシスを持たせても良い。具体的には例えば、第１の閾値と、これよりも大きい第２の閾値を設定しておき、被験者Ｓに体動が生じていないと判定されている状況下においては、単純平均値σｓＡＶが第２の閾値以上となるまで、被験者Ｓに体動が生じたと判定しない。一方で、被験者Ｓに体動が生じていると判定されている状況下においては、単純平均値σｓＡＶが第２の閾値未満となっても被験者Ｓに体動が生じていないとは判定せず、単純平均値σｓＡＶが第１の閾値未満となった時点で、被験者Ｓに体動が生じていないと判定する。
【０１３４】
上記実施形態の生体状態モニタリングシステム１００、２００において、体動判定部３３は、規格化標準偏差σｓ１〜σｓ４を求めることなく、所定の閾値と、標準偏差σ１〜σ４の少なくとも１つ、又は標準偏差σ１〜σ４の少なくとも２つ以上の合計値若しくは単純平均値との比較に基づいて、被験者Ｓの体動の有無を判定しても良い。
【０１３５】
上記実施形態の生体状態モニタリングシステム１００、２００において、体動判定部３３は、標準偏差σ１〜σ４を用いることなく、被験者Ｓの体動の有無を判定することもできる。具体的には例えば、被験者Ｓの重心Ｇの移動に基づいて、被験者Ｓの体動の有無を判定する。
【０１３６】
前述の通り、被験者Ｓの重心Ｇは、被験者Ｓの呼吸に応じて、被験者Ｓの体軸ＳＡの方向に振動している（図５（ａ））。また、被験者Ｓの重心Ｇは、被験者Ｓに小さな体動や大きな体動が生じた場合には、これに応じて移動する。そして、所定期間における重心Ｇの移動距離は、被験者Ｓが呼吸のみを行っている期間、被験者Ｓが小さな体動を行っている期間、被験者Ｓが大きな体動を行っている期間の順に大きくなる。
【０１３７】
従って体動判定部３３は、所定期間における重心Ｇの移動距離を、所定の閾値と比較することにより、被験者Ｓに体動が生じているか否かを判定することができる。具体的には例えば、所定期間における重心Ｇの移動距離ｄが、所定の閾値ｄｔｈよりも大きい場合に、被験者Ｓに小さな体動が生じていると判定し得る。
【０１３８】
上記実施形態の生体状態モニタリングシステム１００、２００は、必ずしも荷重検出器１１〜１４の全てを備える必要はなく、このいずれか一つを備えるのみでもよい。例えば荷重検出器が３つの場合には、これを一直線に配置しなければ、ベッドＢＤ面上での被験者Ｓの重心位置Ｇを検出できる。また、荷重検出器は、必ずしもベッドの四隅に配置される必要はなく、ベッド上の被験者の荷重及びその変動を検出しうるように、任意の位置に配置し得る。また、荷重検出器１１〜１４は、ビーム形ロードセルを用いた荷重センサに限られず、例えばフォースセンサを使用することもできる。
【０１３９】
上記実施形態の生体状態モニタリングシステム１００、２００においては、荷重検出器１１〜１４の各々は、ベッドＢＤの脚の下端に取り付けられたキャスターＣの下に配置されていたがこれには限られない。荷重検出器１１〜１４の各々は、ベッドＢＤの４本の脚とベッドＢＤの床板との間に設けられてもよいし、ベッドＢＤの４本の脚が上下に分割可能であれば、上部脚と下部脚との間に設けられても良い。また、荷重検出器１１〜１４をベッドＢＤと一体に又は着脱可能に組み合わせて、ベッドＢＤと第１、第２実施形態の生体状態モニタリングシステム１００、２００とからなるベッドシステムＢＤＳを構成してもよい（図１３）。
【０１４０】
上記実施形態の生体状態モニタリングシステム１００、２００において、荷重検出部１とＡ／Ｄ変換部２との間に、荷重検出部１からの荷重信号を増幅する信号増幅部や、荷重信号からノイズを取り除くフィルタリング部を設けても良い。
【０１４１】
上記実施形態の生体状態モニタリングシステム１００、２００において、表示部５は、モニターに代えて、又はこれに加えて、生体状態（呼吸数、体動の有無等）を表わす情報を印字して出力するプリンタや、生体状態（呼吸数、体動の有無等）を表示するランプ等の簡易な視覚表示手段を備えてもよい。報知部６はスピーカーに代えて、又はこれに加えて、振動により報知を行う振動発生部を備えてもよい。
【０１４２】
本発明の特徴を維持する限り、本発明は上記実施の形態に限定されるものではなく、本発明の技術的思想の範囲内で考えられるその他の形態についても、本発明の範囲内に含まれる。
【産業上の利用可能性】
【０１４３】
本発明の生体状態モニタリングシステムによれば、被験者の体動の影響を抑制して、被験者の呼吸数を高い精度で推定することができる。したがって、本発明の生体状態モニタリングシステムを用いれば、高精度な呼吸数の推定に基づく、質の高い医療、介護を提供することができる。
【符号の説明】
【０１４４】
１  荷重検出部、１１，１２，１３，１４  荷重検出器、２  Ａ／Ｄ変換部、３  制御部、３１  標準偏差算出部、３２  呼吸波形描画部、３３  体動判定部、３４  呼吸数推定部、３５  トゥイッチ判定部、４  記憶部、５  表示部、６  報知部、７  入力部、１００  生体状態モニタリングシステム、ＢＤ  ベッド、ＢＤＳ  ベッドシステム、Ｓ  被験者
【図１】
2019187941.tif 000006
【図２】
2019187941.tif 000007
【図３】
2019187941.tif 000008
【図４】
2019187941.tif 000009
【図５】
2019187941.tif 000010
【図６】
2019187941.tif 000011
【図７】
2019187941.tif 000012
【図８】
2019187941.tif 000013
【図９】
2019187941.tif 000014
【図１０】
2019187941.tif 000015
【図１１】
2019187941.tif 000016
【図１２】
2019187941.tif 000017
【図１３】
2019187941.tif 000018