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