(57)【要約】本発明は、複数の心拍数を決定し、複数の心拍数計算のうちのただ１つのみに基づいて複数の心拍数の間で選択を行なうパルス酸素濃度計に関する。第１の心拍数計算方法が選択され、その測定基準の精度が疑わしい場合を除き、第１の心拍数計算方法が用いられる。上記測定基準の精度が疑わしい場合、代替的な心拍数計算が利用可能であり、代わりに用いられる。一実施形態において、第１の心拍数計算方法は、集合平均された波形を用いないが、代替的な心拍数計算は、集合平均された波形を用いる。
2007527770.tif 000002
【特許請求の範囲】
【請求項１】
パルス酸素濃度計において心拍数を決定する方法であって、
第１の方法を用いることにより、パルス酸素測定信号から第１の心拍数を決定することと、
第２の方法を用いることにより、パルス酸素測定信号から第２の心拍数を決定することと、
該第１の方法に適用される測定基準を用いることにより、該第１の心拍数の信頼度を評価することと、
該第１の方法が信頼できるということを該測定基準が示唆するときに、該第１の心拍数を用いることと、
該第１の心拍数が信頼できないということを該測定基準が示唆するときに、該第２の心拍数を用いることと
を包含する、方法。
【請求項２】
直前のパルスが拒絶されたということを前記測定基準が示唆するときに、前記第１の心拍数は信頼できないということを決定することをさらに包含する、請求項１に記載の方法。
【請求項３】
前記第１の方法は、集合平均された波形を用いず、前記第２の方法は、集合平均された波形を用いる、請求項１に記載の方法。
【請求項４】
第１および第２の心拍数を各々決定することは、パルス周期を決定することを包含し、
用いられるパルス周期をパルス数に変換すること
をさらに包含する、請求項１に記載の方法。
【請求項５】
心拍数を決定するパルス酸素濃度計であって、
第１の方法を用いることにより、パルス酸素測定信号から第１の心拍数を決定する第１の心拍数計算器と、
第２の方法を用いることにより、パルス酸素測定信号から第２の心拍数を決定する第２の心拍数計算器と、
該第１の方法に適用される測定基準を用いることにより、該第１の心拍数の信頼度を決定するように構成された評価器と、
該第１の方法が信頼できるということを該測定基準が示唆するときに、該第１の心拍数を用い、該第１の心拍数が信頼できないということを該測定基準が示唆するときに、該第２の心拍数を用いるように構成された選択器と
を備える、パルス酸素濃度計。
【請求項６】
前記選択器は、直前のパルスが拒絶されたということを前記測定基準が示唆するときに、前記第１の心拍数は信頼できないということを決定する、請求項５に記載のパルス酸素濃度計。
【請求項７】
前記第１の心拍数計算器は、集合平均された波形を用いず、前記第２の心拍数計算器は、集合平均された波形を用いる、請求項５に記載のパルス酸素濃度計。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、酸素濃度計に関し、特に、パルス酸素濃度計から検出された波形における複数の機構によってパルス数を決定することに関する。
【背景技術】
【０００２】
パルス酸素濃度計は、典型的には、動脈血におけるヘモグロビンの血中酸素飽和度、組織に供給する個々の脈拍の血液量、および、患者の各心拍に対応する脈拍の速度を含む様々な血液化学特性を測定するために用いられるが、それらには限定されない。これら特性の測定は、非侵襲センサを用いることによって達成される。上記非侵襲センサは、血液が灌流する患者の組織の一部を介して光を散乱させ、そのような組織における光の吸収を光電子的に感知する。その後、様々な波長の吸収された量の光は、測定される血液成分の量を計算するために用いられる。
【０００３】
組織を介して散乱された光は、血液中に存在する血液成分の量を表す量において、血液によって吸収された１つ以上の波長として選択される。組織を介して散乱された透過光の量は、組織における血液成分の量の変化、および関連する光の吸収量の変化にしたがって、変動し得る。血中酸素レベルを測定するため、そのようなセンサは、血中酸素飽和度を測定するための周知の技術にしたがい、典型的には、少なくとも２つの異なる波長の光を生成することに適合した光源と、それら両波長を感知可能な光検出器とが提供される。
【０００４】
周知の非侵襲センサは、例えば指、耳、または頭皮のような体の一部に固定されるデバイスを含む。動物や人間において、これら体の一部の組織は血液によって灌流され、上記組織の表面は、センサに対して容易にアクセス可能である。
【０００５】
特許文献１、特許文献２、特許文献３は、「最良レート（ｂｅｓｔ  ｒａｔｅ）」モジュールを有するパルス酸素濃度計においてパルス数を計算する複数の方法を示している。「最良レート」モジュールは、複数のパルス数計算の間を取り持ち、各パルス数計算に関連する信頼度に基づいて最良レートを選択する。上記信頼度は、異なるパルス数計算の信頼度を決定する様々な測定基準（ｍｅｔｒｉｃ）を用いて計算される。また、特許文献４は、胎児心拍数モニタを示しており、上記胎児心拍数モニタは、複数の平行なフィルタ経路を用いることにより胎児心拍数を識別し、性能指数を用いることにより異なる心拍数評価に対する重み付け動作を行なう。
【０００６】
（Ｎ−１００）．１９８５年あたりまでのＮ−１００技術は、パルスのサイズ、パルスの形状、（周波数）発生の期待時間、および比Ｒ／ＩＲに関するパルス履歴に基づいて、パルスを受け入れるか拒絶するかしていた。
【０００７】
特に、Ｎ−１００は、後に負の最大傾斜が続き、その後に最小値を取る最大信号を探すことによって、パルスを見付けていた。処理は、「マンチ（ｍｕｎｃｈ）」と呼称される状態機械（ｓｔａｔｅ  ｍａｃｈｉｎｅ）において行なわれていた。各最大値は、ノイズゲート（ｎｏｉｓｅ  ｇａｇｅ）と呼称されるノイズ閾値を信号が下回るまで、適格とはされなかった。ノイズゲートのレベルは、期待される異なる信号振幅に適合するように、その後の処理ステップからのフィードバックによってセットされていたため、このことは、適合可能フィルタとしての役割を果たしていた。その後、パルスは、「レベル３（Ｌｅｖｅｌ３）」処理において受け入れまたは拒絶されていた。「レベル３」処理は、振幅、周期、および、履歴バッファにおける値の平均に対する新しいパルスの比の比率（ｒａｔｉｏ−ｏｆ−ｒａｔｉｏｓ；ＩＲに対するＲｅｄの比であり、ＲｅｄおよびＩＲは、ＤＣに対するＡＣの比として表されている）を比較し、その後、信頼レベル内に差が存在するか否かを決定することによって、変化する信号に適合するフィルタであった。新しいパルスが受け入れられると、履歴バッファは、新しいパルスに対する値を用いて更新されていた。ｌｅｖｅｌ３処理は、中心周波数と帯域幅（信頼限界）とが、フィルタ出力からのフィードバックに適合する、適合可能帯域通過フィルタとしての役割を果たしていた。
【０００８】
（Ｎ−２００）．Ｎ−２００は、ＥＣＧならびに内蔵のＥＣＧフィルタリングと同期化することができたという理由により、Ｎ−１００の改良であった。Ｎ−２００はまた、最大パルスの測定時間と最小パルスの測定時間との間でシフトするベースラインを補償するために、内挿（ｉｎｔｅｒｐｏｌａｔｉｏｎ）を追加していた。Ｎ−２００は、可変個数の信号サンプルの平均を計算する「ボックスカー（ｂｏｘｃａｒ）」フィルタのようなその他のフィルタリング特性もまた、同様に含んでいた。
【０００９】
Ｎ−２００は、様々なフィルタリングステップおよびスケーリングステップの後、デジタル化された信号を「ボックスカー」フィルタに印加し、「ボックスカー」フィルタは、Ｎ個のサンプルの平均を計算する。ここにＮは、フィルタリングされた心拍数にしたがい、その後の処理ステップからのフィードバックによってセットされる。新しいサンプルは、ボックスカーフィルタにおいて平均化されるが、古いサンプルは除外される。ボックスカー長（Ｎ）は、３つのパラメータ、すなわちパルス閾値、絶対最小パルス、および小パルスを決定するために用いられる。その後、集合平均（ｅｎｓｅｍｂｌｅ−ａｖｅｒａｇｉｎｇ；別名「スライダー（ｓｌｉｄｅｒ）」）フィルタは、新しいサンプルと、１パルス周期前の、先に集合平均されたサンプルとの加重平均を計算する。その後、サンプルは、Ｎ−１００のように、「マンチ」状態機械とノイズゲートとを通過する。Ｎ−１００処理に内挿特性が追加され、ベースラインレベルにおける変化を補償する。最小値および最大値は、異なる時間に発生する。変化するベースラインは、最小値を増加または減少させ、最大値を増加も減少もさせない。あるいはその逆に、変化するベースラインは、最大値を増加または減少させ、最小値を増加も減少もさせない。
【００１０】
「集合平均」は、複数のパルスからのサンプルを互いに平均化処理して合成パルスを形成する、ＮＥＬＬＣＯＲの登録商標であるＣ−Ｌｏｃｋの一体部分である。この処理はまた、「心臓ゲート平均（ｃａｒｄｉａｃ−ｇａｔｅｄ  ａｖｅｒａｇｉｎｇ）」としても知られている。上記処理は、各パルスの開始に印をつけるため、「トリガ（ｔｒｉｇｇｅｒ）」イベントを要求する。
【特許文献１】米国特許第６，０８３，１７２号明細書
【特許文献２】米国特許第５，８５３，３６４号明細書
【特許文献３】米国特許第６，４１１，８３３号明細書
【特許文献４】米国特許第５，５２４，６３１号明細書
【発明の開示】
【課題を解決するための手段】
【００１１】
（発明の簡単な要約）
本発明は、複数の心拍数を決定し、複数の心拍数計算のうちのただ１つに関する測定基準に基づいて、複数の心拍数の間から選択を行なう酸素濃度計に関する。第１の心拍数計算方法が選択され、その測定基準の精度が疑わしい場合を除き、第１の心拍数計算方法が用いられる。上記測定基準の精度が疑わしい場合、代替的な心拍数計算が利用可能であり、代わりに用いられる。
【００１２】
一実施形態において、第１の心拍数計算方法は、集合平均された波形を用いないが、代替的な心拍数計算は、集合平均された波形を用いる。第１の計算が、直前に検出された、第１の計算に関するパルスに適格性が付与されなかった場合に、代替的な心拍数計算が用いられる。
【発明を実施するための最良の形態】
【００１３】
（好ましい実施形態に関する詳細な記述）
図１は、本発明を組み込んだ酸素測定システムの一実施形態を図示している。センサ１０は、赤色ＬＥＤおよび赤外線ＬＥＤならびに光検出器を含んでいる。これらは、ケーブル１２によって基板１４に接続されている。ＬＥＤ駆動電流は、ＬＥＤ駆動インターフェース１６によって提供される。センサからの受信された光電流は、Ｉ−Ｖインターフェース１８に提供される。その後、ＩＲの電圧および赤色の電圧は、本発明を組み込んだシグマ・デルタインターフェース２０に提供される。シグマ・デルタインターフェース２０の出力は、マイクロコントローラ２２に提供される。マイクロコントローラ２２は、１０−ビットのＡ／Ｄコンバータを含む。コントローラ２２は、プログラム用のフラッシュメモリと、データ用のＥＥＰＲＯＭとを含む。プロセッサはまた、フラッシュメモリ２６に接続されたコントローラチップ２４を含む。最後に、クロック２８が用いられ、センサ１０におけるデジタル較正へのインターフェース３０が提供される。分離したホスト３２は、処理された情報を受信し、同時に、アナログ表示を提供するために、ライン３４上のアナログ信号を受信する。
【００１４】
（設計の要約）本発明の設計は、不要なノイズに対処することを意図している。信号の測定基準が測定され、フィルタの重み付けを決定するために用いられる。信号の測定基準は、例えば周波数（ヒトの心拍数の範囲に存在するような周波数）、形状（心臓のパルスに似た形状）、立ち上がり時間、等のように、パルスがプレチスモグラフ（ｐｌｅｔｈｙｓｍｏｇｒａｐｈ）またはノイズでありそうか否かを示唆するものである。類似する技術は、本明細書の背景におけるＮｅｌｌｃｏｒ  Ｎ２００において用いられていた。新しい設計は、例えば本発明において請求されている２つの集合平均のような、多くの異なる特徴と改変とを追加している。
【００１５】
アーキテクチャの詳細は、図２の図に示されている。この設計は、酸素飽和度とパルス数とを計算する。これらは、以下で別々に記述される。
【００１６】
（Ｉ．酸素飽和度の計算）
（Ａ．信号の調整）デジタルな赤色信号およびＩＲ信号は、以下により、このブロックにおいて受信および調整される。（１）１次微分を取り、ベースラインのシフトを取り除くこと、（２）一定の係数を用いて低域通過フィルタリングすること、（３）ＤＣ値で割ることにより、比を保つこと。信号調整サブシステムの機能は、ヒトのプレチスモグラフに生じる周波数を強調し、しばしば動きアーチファクトが集中する低周波数を減衰する。信号調整サブシステムは、初期化の間に識別されたハードウィア特性に基づいて、そのフィルタ係数を選択する。
入力−デジタルな赤色信号およびＩＲ信号
出力−前処理された赤色信号およびＩＲ信号。
【００１７】
（Ｂ．パルスの識別と適格性付与）低域通過フィルタリングされたデジタルな赤色信号およびＩＲ信号は、パルスを識別し、それらが動脈拍動に類似しているという適格性を付与するために、このブロックに提供される。このことは、あらかじめ訓練を受けた（ｐｒｅ−ｔｒａｉｎｅｄ）神経網を用いて行なわれ、主にＩＲ信号の上で行なわれる。パルスは、Ｎｅｌｌｃｏｒ  Ｎ−１００で行なわれたのと同様に、その振幅、形状および周波数を検討することによって、識別される。このブロックへの入力は、ブロックＤからの平均パルス周期である。この機能は、Ｎ−１００に類似しており、Ｎ−１００は、パルス数を用いて前もっての適格性付与を変更していた。出力は、不整脈の程度と個々のパルスの品質とを示唆する。
入力−（１）前処理された赤色信号およびＩＲ信号、（２）平均パルス周期、（３）低域通過フィルタからの低域通過波形。
出力−（１）不整脈の程度、（２）パルスの振幅の変動、（３）個々のパルスの品質、（４）パルスのビープ通知、（５）適格性付与されたパルスの周期およびエイジ（ａｇｅ）。
【００１８】
（Ｃ．信号品質の測定基準の計算）このブロックは、パルスの形状（微分スキュー；ｄｅｒｉｖａｔｉｖｅ  ｓｋｅｗ）、周期の可変性、パルスの振幅および可変性、比の比率の可変性、およびパルス数に関連する周波数コンテンツを決定する。
入力−（１）未加工のデジタルな赤色信号およびＩＲ信号、（２）不整脈の程度、個々のパルスの品質、パルスの振幅の変動、（３）前処理された赤色信号およびＩＲ信号、（４）平均パルス周期。
出力−（１）低域通過および集合平均のフィルタ重み、（２）センサ・オフ検出器に対する測定基準、（３）正規化された前処理波形、（４）パーセント変調（ｐｅｒｃｅｎｔ  ｍｏｄｕｌａｔｉｏｎ）。
【００１９】
（Ｄ．平均パルス周期）このブロックは、受信したパルスから平均パルス周期を計算する。
入力−適格性付与されたパルスの周期およびエイジ。
出力−平均パルス周期。
【００２０】
（Ｅ１．低域通過フィルタおよび集合平均）ブロックＥ１は、ブロックＡにより調整され、ブロックＣにより正規化された信号を、パルス数を識別するために、低域通過フィルタし、集合平均する。低域通過フィルタに対する重みは、信号測定基準ブロックＣによって決定される。信号はまた、集合平均され（このことは、パルス数付近の興味のある周波数以外の周波数とその高調波とを減衰する）、集合平均フィルタの重みはまた、信号測定基準ブロックＣによって決定される。信号が低下しているというフラグが付与されると、低い重みが割り当てられる。信号に不整脈であるというフラグが付与されると、集合平均は不整脈の間は適切ではないため、より高い重みが割り当てられる。赤色およびＩＲは、別々に処理されるが、同一のフィルタ重みで処理される。信号の測定基準が第１に計算されるようにするため、フィルタリングは、約１秒遅延される。
【００２１】
フィルタは、連続的に変動する重みを用いる。サンプルが集合平均されない場合、既にフィルタリングされたサンプルに対する重みは、加重平均において、ゼロにセットされ、さらに、新しいサンプルは、コードを介して処理される。このブロックは、信号のエイジ、すなわちフィルタリングの累積（処理における応答時間と遅延との和）を追跡する。あまりに古い結果は、（良いパルスがしばらく検出されなかった場合に）フラグが付与され得る。
入力−（１）正規化され前処理された赤色信号およびＩＲ信号、（２）平均パルス周期、（３）低域通過フィルタの重みおよび集合平均フィルタの重み、（４）利用可能な場合は、ＥＣＧトリガ、（５）ゼロ交差のトリガに対するＩＲ基本波。
出力−（１）フィルタリングされた赤色信号およびＩＲ信号、（２）エイジ。
【００２２】
（Ｆ．フィルタリングされた波形の相関に関する評価および平均重みの計算）これは、上述のＮ１００およびＮ２００に用いられるのと同様なノイズの測定基準を用いるが、フィードバックは用いない。フィルタに対する可変の重み付けは、比の比率の差異によって制御される。この可変重みのフィルタリングの効果は、アーチファクトが増大する際には比の比率がゆっくりと変化し、アーチファクトが減少する際には比の比率が迅速に変化するということである。サブシステムは、２つの応答モードを有する。高速モードにおけるフィルタリングは、エイジが３秒である測定基準を対象とする。正規モードでは、対象となるエイジは５秒である。高速モードにおいて、現在の値の最小の重み付けは、高いレベルに短縮される。言い換えると、ノイズが存在する場合は低い重みが新しい比の比率計算に割り当てられ、ノイズが存在しない場合は高い重みが割り当てられる。
入力−（１）フィルタリングされた赤色信号およびＩＲ信号ならびにエイジ、（２）較正係数、（３）応答モード（ユーザ速度設定）。
出力−比の比率計算に対する重みの平均化。
【００２３】
（Ｈ．飽和度の計算）飽和度は、較正係数と平均化された比の比率を有するアルゴリズムを用いて計算される。
入力−（１）平均化された比の比率、（２）較正係数。
出力−飽和度。
【００２４】
（ＩＩ．パルス数の計算）
（Ｅ２．低域通過フィルタおよび集合平均）ブロックＥ２は、ブロックＡにより調整された信号を、パルス数を識別するために、低域通過フィルタし、集合平均する。低域通過フィルタに対する重みは、信号測定基準ブロックＣによって決定される。信号はまた、集合平均され（このことは、パルス数付近の興味のある周波数以外の周波数と高調波とを減衰する）、集合平均フィルタの重みはまた、信号測定基準ブロックＣによって決定される。信号が低下しているというフラグが付与されると、低い重みが割り当てられる。信号に不整脈であるというフラグが付与されると、フィルタリングは不整脈の間は適切ではないため、より高い重みが割り当てられる。赤色およびＩＲは、別々に処理される。このブロックの処理は、信号の測定基準が第１の計算されるようにするため、約１秒遅延される。
【００２５】
フィルタは、連続的に変動する重みを用いる。サンプルが集合平均されない場合、既にフィルタリングされたサンプルに対する重みは、加重平均において、ゼロにセットされ、さらに、新しいサンプルは、コードを介して処理される。このブロックは、信号のエイジ、すなわちフィルタリングの累積（処理における応答時間と遅延との和）を追跡する。あまりに古い結果は、（良いパルスがしばらく検出されなかった場合に）フラグが付与され得る。
入力−（１）前処理された赤色信号およびＩＲ信号、（２）平均パルス周期、（３）低域通過フィルタの重みおよび集合平均フィルタの重み、（４）利用可能な場合は、ＥＣＧトリガ、（５）ゼロ交差のトリガに対するＩＲ基本波。
出力−（１）フィルタリングされた赤色信号およびＩＲ信号、（２）エイジ。
【００２６】
（Ｉ．フィルタリングされたパルスの識別と適格性付与）このブロックは、フィルタリングされた波形からのパルス周期を識別し、適格性付与を行なう。また、この結果は、パルスがブロックＢにより適格性を付与されなかった場合にのみ用いられる。
入力−（１）フィルタリングされた赤色信号およびＩＲ信号ならびにエイジ、（２）
平均パルス周期、（３）ハードウェアＩＤまたはノイズフロア。
出力−適格性付与されたパルスの周期およびエイジ。
【００２７】
（Ｊ．平均パルス周期およびパルス数計算）このブロックは、パルス数と平均パルス周期とを計算する。
入力−適格性付与されたパルスの周期およびエイジ。
出力−（１）平均パルス周期、（２）パルス数。
【００２８】
（ＩＩＩ．静脈拍動）
（Ｋ．静脈拍動の検出）ブロックＫは、ブロックＡからの前処理された赤色信号およびＩＲ信号ならびにエイジと、パルス数とを、入力として受信し、出力として、静脈拍動に関する指標を提供する。このサブシステムは、単一歯くし形フィルタ（ｓｉｎｇｌｅ−ｔｏｏｔｈ  ｃｏｍｂ  ｆｉｌｔｅｒ）を用いることにより、集合平均フィルタへの出力であるＩＲ基本波形を、時間領域に生成する。
入力−（１）フィルタリングされた赤色信号およびＩＲ信号ならびにエイジ、（２）パルス数。
出力−静脈拍動の指標、ＩＲ基本波。
【００２９】
（ＩＶ．センサ・オフ）
（Ｌ．センサ・オフの検出およびパルス振幅の損失）パルス損失およびセンサ・オフ検出サブシステムは、あらかじめ訓練を受けた神経網を用い、センサが患者から外れている（ｏｆｆ）か否かを決定する。神経網への入力は、数秒間にわたるＩＲおよび赤色の値の挙動に関するいくつかの局面を定量化する測定基準である。信号状態がパルスが存在する状態（Ｐｕｌｓｅ  Ｐｒｅｓｅｎｔ）またはセンサがオフであり得る状態（Ｓｅｎｓｏｒ  Ｍａｙｂｅ  Ｏｆｆ）のいずれでもない間、サンプルは、多数の酸素濃度計アルゴリズムのサブシステムによって無視される。信号状態の変数の値は、「パルスが存在する状態、不連続状態、パルス損失状態、センサがオフであり得る状態、およびセンサ・オフ状態」である。
入力−（１）測定基準、（２）フロントエンドのサーボ設定およびＩＤ。
出力−センサ・オフ指標を含む信号状態。
【００３０】
（パルス数計算サブシステム）サブシステムは、パルスの識別と適格性付与に関するサブシステムから適格性付与されたパルス周期を平均化する。上記サブシステムは、平均周期と対応するパルス数とを出力する。
【００３１】
酸素測定アルゴリズムは、このサブシステムの２つの例を含む。第１の例は、パルスの識別と適格性付与に関する例から、入力を受信する。その入力波形は、信号調整サブシステムによって処理され、その後低域通過フィルタリングされるが、集合平均サブシステムによって集合平均されない。パルス数計算サブシステムの第２の例は、パルス数の識別と適格性付与のサブシステムの２つの例から、入力を受信する。第１の例は、上記で記述されており、入力を受信する第２の例は、集合平均される。
【００３２】
（パルス周期のソースの選択）サブシステムの一例は、２つのソースから、適格性付与されたパルスを受信する。サブシステムは、１つのソース、すなわち「第１」ソースの解析だけに基づいて、これら２つのソースのうちのどちらをパルス数計算に用いるかを選択する。酸素測定アルゴリズムは、集合平均された波形を受信しないパルス識別と適格性付与の例を第１のソースとして指定し、パルスの識別と適格性付与に関するその他の例を適格性付与されたパルス周期の「代替的」ソースとして指定する。代替的なソースからの、適格性付与されたパルス周期は、第１ソースからの直前のパルスが拒絶されたか否かに対してのみ用いられる。適格性付与されたパルス周期が第１ソースによって受信されるとき、上記パルス周期は、常にパルス数計算を更新するために用いられ、第１ソースが再度パルス周期を拒絶する時点まで、代替的ソースからの適格性付与されたパルス周期が用いられるのを防ぎ得る。
【００３３】
（平均パルス周期の計算とパルス数の評価）サブシステムは、Ｑｕａｌｉｆｉｅｄ＿Ｐｕｌｓｅ＿Ｐｅｒｉｏｄを用いるとき、パルスベースの可変重みＩＩＲフィルタを用いることにより、平均パルス周期Ａｖｇ＿Ｐｅｒｉｏｄを更新し、その後、Ａｖｇ＿Ｐｅｒｉｏｄから、Ｒａｔｅ出力を計算する。このフィルタリング操作に対するステップは、
【００３４】
【数１】
2007527770.tif 000003
である。ここに、ＢＰＭにおいて、ｒｔはＱｕａｌｉｆｉｅｄ＿Ｐｕｌｓｅ＿Ｐｅｒｉｏｄに対応するパルス数であり、添え字ｔ−１は先に適格性付与されたパルスを示す。Δｔは、酸素測定アルゴリズムの秒単位でのサンプル間隔である。６０／Δｔは、１分あたりのサンプル数である。ｘは、典型的な成人のパルス数に対する７秒の応答時間を対象とするフィルタ重みである。
【００３５】
ｋは、最終フィルタの重みであり、ｘと連続的な値ｒｔの間の差分との両方に基づいている。はじめの数パルスの間、適格性付与された最初のパルスが等しく重み付けされ得るように、ｋは少なくとも１／Ｔｏｔａｌ＿Ｑｕａｌｉｆｉｅｄまで増加する。
【００３６】
Ｃｏｎｓｅｃｕｔｉｖｅ＿Ｑｕａｌｉｆｉｅｄは、連続的に適格性付与されたパルスの個数であり、Ｔｏｔａｌ＿Ｑｕａｌｉｆｉｅｄは、サブシステムが再初期化されてから適格性付与されたパルスの総数である。ｋが計算される前において、Ｃｏｎｓｅｃｕｔｉｖｅ＿ＱｕａｌｉｆｉｅｄとＴｏｔａｌ＿Ｑｕａｌｉｆｉｅｄとの両方は、Ｑｕａｌｉｆｉｅｄ＿Ｐｕｌｓｅ＿Ｐｅｒｉｏｄが用いられるたびに、増加させられる。Ｃｏｎｓｅｃｕｔｉｖｅ＿Ｑｕａｌｉｆｉｅｄは、現在使用中のパルス周期のソースによってパルスが拒絶されるときに、ゼロにセットされる。
【００３７】
上記のステップ６におけるＡｖｇ＿Ｐｅｒｉｏｄｔの更新に関する公式は、Ａｖｇ＿ＰｅｒｉｏｄｔとＱｕａｌｉｆｉｅｄ＿Ｐｕｌｓｅ＿Ｐｅｒｉｏｄとの幾何平均である。上記幾何平均は、サブシステムが、パルス対パルス（ｐｕｌｓｅｓ−ｔｏ−ｐｕｌｓｅ）の期間に関する大きな変動に対して、さらにはパルス数における持続性の大きな変化に対して、応答可能になるようにする。
【００３８】
一旦、Ｒａｔｅが非ゼロ値に初期化されると、Ｒａｔｅが更新されるか否かに関わらず、すべてのサンプルについてＲａｔｅ＿Ａｇｅは増加する。
【００３９】
（文脈図）図３は、パルス数計算サブシステムの文脈図である。上記サブシステムは、Ｑｕａｌｉｆｉｅｄ＿Ｐｕｌｓｅ＿ＰｅｒｉｏｄｓからのＡｖｇ＿ＰｅｒｉｏｄとＲａｔｅ出力とを更新する。上記サブシステムは、Ｎｏｔｉｆｙ＿Ｐｕｌｓｅ＿Ｒｅｊｅｃｔｅｄを第１ソースから最後に受信した場合に限り、Ａｌｔｅｒｎａｔｉｖｅ＿Ｐｅｒｉｏｄ＿ＳｏｕｒｃｅからのＱｕａｌｉｆｉｅｄ＿Ｐｕｌｓｅ＿Ｐｅｒｉｏｄｓを用いる。上記サブシステムは、Ｑｕａｌｉｆｉｅｄ＿Ｐｅｒｉｏｄ＿Ａｇｅに基づき、Ｒａｔｅ＿Ａｇｅ出力を更新する。Ｒａｔｅが更新されるとき、サブシステムはＰｕｌｓｅ＿Ｒａｔｅ＿Ｕｐｄａｔｅｄフラグをセットする。Ｒｅｉｎｉｔｉａｌｉｚｅは、サブシステムに対して自らを再初期化するように命令する。一旦Ｒａｔｅが初期化されると、Ｉｎｃｒｅｍｅｎｔ＿Ｒａｔｅ＿Ａｇｅは、サブシステムに対してすべてのサンプルについてＲａｔｅ＿Ａｇｅを増加させるように通知する。
【図面の簡単な説明】
【００４０】
【図１】図１は、本発明の実施形態を組み込んだ酸素濃度計のブロック図である。
【図２】図２は、本発明の実施形態を含む酸素濃度計のソフトウェア処理ブロックの図である。
【図３】図３は、パルス数計算サブシステムの文脈図である。
【図１】
2007527770.tif 000004
【図２】
2007527770.tif 000005
【図３】
2007527770.tif 000006
【国際調査報告】
2007527770.tif 000007
2007527770.tif 000008
2007527770.tif 000009