(57)【要約】【課題】胸部圧迫に起因する生理学的信号中の信号アーチファクトを低減すること。【解決手段】胸部圧迫適用中の生理学的信号を解析する装置は、胸部圧迫適用中の生理学的信号を取得する手段と、胸部圧迫の速度を検出するように構成された速度センサーの出力を取得する手段と、前記胸部圧迫に起因する前記生理学的信号中の少なくとも一つの信号アーチファクトを低減するために前記速度センサーの出力を使用する手段とを備える。【選択図】図１
2014195750.tif 000002
【特許請求の範囲】
【請求項１】
胸部圧迫適用中の生理学的信号を取得する手段と、
胸部圧迫の速度を検出するように構成された速度センサーの出力を取得する手段と、
前記胸部圧迫に起因する前記生理学的信号中の少なくとも一つの信号アーチファクトを低減するために前記速度センサーの出力を使用する手段とを備える、胸部圧迫適用中の生理学的信号を解析する装置。
【請求項２】
前記生理学的信号はＥＣＧ信号である請求項１に記載の装置。
【請求項３】
前記生理学的信号はＩＰＧである請求項１に記載の装置。
【請求項４】
前記生理学的信号はＩＣＧ信号である請求項１に記載の装置。
【請求項５】
前記生理学的信号はパルス酸素濃度信号である請求項１に記載の装置。
【請求項６】
前記生理学的信号中の少なくとも一つの信号アーチファクトを低減するために前記速度センサーの出力を使用する手段は、前記生理学的信号と前記速度との時間整合を行う請求項１または２に記載の装置。
【請求項７】
前記生理学的信号中の少なくとも一つの信号アーチファクトを低減するために前記速度センサーの出力を使用する手段は、胸部圧迫アーチファクトを低減するように調節された適応フィルタを含む請求項１または２に記載の装置。
【請求項８】
アーチファクトの低減された前記生理学的信号を処理して心室細動が存在するか否かを判定する心室細動検出アルゴリズムを更に備える請求項１または２に記載の装置。
【請求項９】
胸部圧迫が行われた時を検出し自動的に前記適応フィルタを始動する前処理手段を更に備える請求項７に記載の装置。
【請求項１０】
前記アルゴリズムが心室細動の存在を推定した場合に除細動ショックの伝達を可能にすることを特徴とする請求項８に記載の装置。
【請求項１１】
前記適応フィルタに供給した生理学的信号と、前記適応フィルタによるアーチファクト低減後の生理学的信号との差を表す差分信号が生成される、請求項７に記載の装置。
【請求項１２】
前記差分信号は前記生理学的信号中のアーチファクト量を表す測度となる請求項１１に記載の装置。
【請求項１３】
前記生理学的信号の以後の処理の修正に前記差分信号を使用する請求項１２に記載の装置。
【請求項１４】
前記アーチファクト量が第一の閾値を超えていることを前記差分信号が示す場合、前記心室細動検出アルゴリズムは、アーチファクトの影響に対する耐性を高めるように修正される請求項１３に記載の装置。
【請求項１５】
前記アーチファクト量が前記第一の閾値より高い第二の閾値を超えていることを前記差分信号が示す場合、前記心室細動検出アルゴリズムの使用は停止される請求項１４に記載の装置。
【請求項１６】
前記差分信号に対するスペクトル解析を行い、前記スペクトル解析の結果に基づいて前記生理学的信号のフィルタリングを調整する請求項１３に記載の装置。
【請求項１７】
適応フィルタに供給するに先立ち、前記速度センサーの出力は正規化前処理を受ける請求項７に記載の装置。
【請求項１８】
前記適応フィルタはＦＩＲフィルタを含む請求項７に記載の装置。
【請求項１９】
前記適応フィルタはゼロ次フィルタを含む請求項１８に記載の装置。
【請求項２０】
前記適応フィルタは前記生理学的信号の推定値によって動的に制御される係数を備える請求項７に記載の装置。
【請求項２１】
前記適応フィルタはフィルタ出力と測定された生理学的信号との差が閾値を超えているときに自動的にリセットされる機能を備える請求項７に記載の装置。
【請求項２２】
自動的なリセットはステップサイズを動的に変更してフィルタの収束と安定性との関係を改善する機能を備える請求項２１に記載の装置。
【請求項２３】
前記生理学的信号と前記速度センサーの出力とに対して整合させる時間整合処理手段を更に備える請求項１または２に記載の装置。
【請求項２４】
前記時間整合処理手段の出力に対して、前記生理学的信号と前記速度センサーの出力との間の誤差を低減させる適応フィルタ処理手段を更に備える請求項２３に記載の装置。
【請求項２５】
前記適応フィルタはカルマンフィルタを含む請求項７に記載の装置。
【請求項２６】
前記適応フィルタは適応等化を用いる請求項７に記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は心臓蘇生を支援する装置に関する。
【背景技術】
【０００２】
心拍停止した患者の蘇生処置には、患者の気道確保、患者に対する人工呼吸、患者の心臓、脳およびその他の生命維持に必要な器官への血流確保のための胸部圧迫が含まれる。患者がショックの効く心拍リズム（shockable heart rhythm）を有する場合には、蘇生として除細動治療が含まれることもある。一次救命処置（ＢＬＳ）という用語には次に挙げる要素のすべてが含まれる。初期診断、気道確保、呼気換気（人工呼吸）、胸部圧迫。これら３つの処置（気道による呼吸、および胸部圧迫を含む循環）の組み合わせには、心肺蘇生術（ＣＰＲ）という用語が使われる。
【発明の開示】
【発明が解決しようとする課題】
【０００３】
現在行われている自動化ＢＣＧ拍動分析法では、心肺蘇生（ＣＰＲ）を中断することにより、胸部の圧迫に由来するアーチファクトがＥＣＧに混入することを回避している。しかしＣＰＲを長時間中断すると蘇生成功率が低下することが明らかにされている。前胸部圧迫を中断すると自発循環回復率および２４時間生存率がいちじるしく低下するという研究が報告されている。Ｙ．Ｓａｔｏ、ＭＨ．Ｗｅｉｌ、Ｓ．Ｓｕｎ、Ｗ．Ｔａｎｇ、Ｊ．Ｘｉｅ、Ｍ．Ｎｏｃ、およびＪ．Ｂｉｓｅｒａ、心肺蘇生術時における前胸部圧迫中断の及ぼす悪影響について、Ｃｒｉｔｉｃａｌ  Ｃａｒｅ  Ｍｅｄｉｃｉｎｅ、２５巻（５）、７３３〜７３６（１９９７年）。Ｙｕほか、２００２年。Ｃｉｒｃｕｌａｔｉｏｎ、１０６、３６３３７２（２００２年）、Ｔ．Ｅｆｔｅｓｔｏｌ、Ｋ．Ｓｕｎｄｅ、およびＰ．Ａ．Ｓｔｅｅｎ、院外心拍停止時における除細動成功の確率試算に対し前胸部圧迫中断の与える影響、Ｃｉｒｃｕｌａｔｉｏｎ、１０５、２２７０〜２２７３、（２００２年）。
【０００４】
ＣＰＲ処置の別の重要問題に呼吸管理がある。呼吸モニタリングの典型的な方法としては、患者の経胸インピーダンスの変化を測定し追跡するインピーダンス呼吸曲線記録法が使用されている。しかし現在では、インピーダンスによる呼吸曲線記録法は、胸部を圧迫するとインピーダンス信号にいちじるしいアーチファクトが生じるため、胸部圧迫適用中の肺容量の指標としては信頼性に欠けている。
【０００５】
胸部圧迫に起因してＥＣＧに現れるアーチファクトを除去する手段として適応フィルタの使用が試みられている。ＳＯ．Ａｓｓｅ、Ｔ．Ｅｆｆｅｓｔｏｌ、ＪＨ．Ｈｕｓｏｙ、Ｋ．Ｓｕｎｄｅ、および、ＰＡ．Ｓｔｅｅｎ、管内蘇生中の心室細動時におけるＣＰＲアーチファクトの提言について、ＩＥＥＥ  Ｔｒａｎｓａｃｔｉｏｎｓ  ｏｎ  Ｂｉｏｍｅｄｉｃａｌ  Ｅｎｇｉｎｅｅｒｉｎｇ、４７巻、１４４０〜１４４９（２０００年）。Ａ．Ｌａｎｇｈｅｌｌｅ、Ｔ．Ｅｆｔｅｓｔｏｌ、Ｈ．Ｍｙｋｌｅｂｕｓｔ、Ｍ．Ｅｒｉｋｓｅｎ、ＭＴ．Ｈｏｌｔｅｎ、ＰＡ．Ｓｔｅｅｎ、最適多チャンネルフィルタを用いた心電図からのＣＰＲアーチファクト除去、Ｍａｒ；４８（３）：２７９９１（２００１年）。ＪＨ．Ｈｕｓｏｙ、Ｊ．Ｅｉｌｅｖｓｔｊｏｎｎ、Ｔ．Ｅｆｔｅｓｔｏｌ、ＳＯ．Ａｓｓｅ、Ｈ．Ｍｙｋｌｅｂｕｓｔ、および、Ｐ．Ａ．Ｓｔｅｅｎ、効率的に適合する追跡的アルゴリズムを用いた心電図からの心肺蘇生術アーチファクトの除去、ＩＥＥＥ  Ｔｒａｎｓａｃｔｉｏｎｓ  ｏｎ  Ｂｉｏｍｅｄｉｃａｌ  Ｅｎｇｉｎｅｅｒｉｎｇ、４９巻、１２８７〜１２９８（２００２年）。ＨＲ．Ｈａｌｐｅｒｉｉｎ、および、ＲＤ．Ｂｅｒｇｅｒ、ＣＰＲ  Ｃｈｅｓｔ  Ｃｏｍｐｒｅｓｓｉｏｎ  Ｍｏｎｉｔｏｒ、米国特許第６
，３９０，９９６号（２００２年）。Ａｓｓｅほか（２００２年）、および、Ｌａｎｇｈｅｌｌｅほか（２００１年）は適応フィルタに対する基準信号として圧迫深さおよび胸部インピーダンスを使用した。Ｈｕｓｏｙほか（２００２年）は、計算の複雑さを軽減するためマッチング追跡反復法を使ってこの研究を展開した。しかし、彼らが得た結果は、例えば高次の逆フィルタの計算を含むようなたいていが計算に集中したものである。Ｈａｌｐｅｒｉｎら（２００２年）は、信号の自動スペクトルおよび交差スペクトルを使う周波数ドメイン法と、ＥＣＧ信号の適応フィルタリングに対して回帰最小二条法を使う時間ドメイン法を提案した。両者とも大量の計算を必要とする。
【０００６】
適応フィルタに関しては多くの文献がある。例えば、Ｈｙｋｉｎ，Ａｄａｐｔｉｖｅ  Ｆｉｌｔｅｒ  Ｔｈｅｏｒｙ、第３版、Ｕｐｐｅｒ  Ｓａｄｄｌｅ  Ｒｉｖｅｒ、ＮＪ，ＵＳＡ．Ｐｒｅｎｔｉｃｅ−Ｈａｌｌ、１９９６年。
【課題を解決するための手段】
【０００７】
概して本発明は胸部圧迫適用中の生理学的信号（例えばＥＣＧ）を解析する方法に関する。本発明の一態様に従う方法は、胸部圧迫適用中の生理学的信号を取得し、そこから胸部圧迫の速度に関する情報を決定可能なセンサーの出力を取得し、胸部圧迫に起因する生理学的信号中の少なくとも一つの信号アーチファクトを低減するために速度に関する情報を使用することを含む。
【０００８】
本発明の好ましい実施例には以下の一つ以上を組み込むことができる。生理学的信号は、ＥＣＧ信号、ＩＰＧ信号、ＩＣＧ信号またはパルス酸素濃度信号などの各種生理学的信号のいずれであってもよい。センサーには速度センサーを使用することができ、速度に関する情報は速度センサーから決定することができる。センサーには加速度計を使用することもでき、速度情報は、加速度計の出力積分から決定してもよい。生理学的信号中の少なくとも一つの信号アーチファクトを低減するために速度に関する情報を使用することには、生理学的信号と速度との時間整合を含めることができる。生理学的信号中の少なくとも一つの信号アーチファクトを低減するために速度に関する情報を使用することは、胸部圧迫アーチファクトを低減するように調節された適応フィルタを使用することを含めることができる。その方法は、アーチファクトの低減された生理学的信号を処理して心室細動が存在するか否かを判定する心室細動検出アルゴリズムを含むことができる。その方法は、胸部圧迫が行われた時を検出し自動的に適応フィルタを始動する前処理工程を含むことができる。その方法は、アルゴリズムが心室細動の存在を推定した場合に除細動ショックの伝達を可能にすることを含むことができる。適応フィルタに供給した生理学的信号と、適応フィルタによるアーチファクト低減後の生理学的信号との差を表す差分信号を生成することができる。差分信号は生理学的信号中のアーチファクト量を表す測度とすることができる。生理学的信号の以後の処理の修正に差分信号を使用することができる。アーチファクト量が第一の閾値を超えていることを差分信号が示す場合、アーチファクトの影響に対する耐性を高めるように心室細動検出アルゴリズムを修正することができる。アーチファクト量が第一の閾値より高い第二の閾値を超えていることを差分信号が示す場合、心室細動検出アルゴリズムの使用を中止することができる。差分信号に対するスペクトル解析を行い、スペクトル解析の結果に基づいて生理学的信号のフィルタリングを調整することができる。適応フィルタに供給するに先立ち、速度信号に正規化前処理を行うことができる。適応フィルタはＦＩＲフィルタを含むことができる。適応フィルタはゼロ次フィルタを含むことができる。適応フィルタは生理学的信号の推定値によって動的に制御される係数を備えることができる。適応フィルタはフィルタ出力と測定された生理学的信号との差が閾値を超えているときに自動的にリセットされる機能を備えることができる。自動的なリセットはステップサイズを動的に変更してフィルタの収束と安定性との関係を改善する機能を備えることができる。生理学的信号と速度信号に対して行われ、これら二つの信号を圧迫に対して整合させる時間整合処理を更に備えてもよい。その方法は、時間整合処理の
出力に対して、生理学的信号と速度信号との間の誤差を低減させる適応フィルタ処理を行うことを含むことができる。適応フィルタはカルマンフィルタを含むことができる。適応フィルタは適応等化を用いることがえきる。
【０００９】
本発明の多くの利点（その一部は、本発明のさまざまな実施例のいくつかでのみ達成することができよう）のうちの一部は以下のとおりである。
本発明は、優れた技術、具体的には（ａ）蘇生術の実施によって誘発されるＥＣＧ信号中のアーチファクトを適応処理によって除去する技術、（ｂ）モニタリングするＥＣＧ信号を増幅する技術、および（ｃ）ＥＣＧのリズムアドバイザリアルゴリズムの信頼性を高める技術を提供する。
【００１０】
リズムアドバイザリアルゴリズムの一部として、本発明のさまざまな実施例をＥＣＧモニター、外部除細動器、ＥＣＧリズム分類器、または心室不整脈検出器の中に組み込むことができる。
【００１１】
本発明は、リズムアドバイザリアルゴリズムのためにＥＣＧデータを収集しながらＣＰＲを続行することができる。これによりＣＰＲの成績が向上し、ひいては例えば生存率が向上する。
【００１２】
さらに本発明は「浄化した」ＥＣＧ信号出力を除細動器の使用者への表示に提供することができる。
さらに本発明は、インピーダンス法によって胸部圧迫適用中の肺容積を測定する手段を初めて提供する。本発明は圧迫に由来するアーチファクトによって改悪された、例えばインピーダンス心拍記録やパルス酸素濃度計のような別の生理学的信号のフィルタリングに使用することができる
本発明は、ゼロ次ＦＩＲフィルタを使ったＣＰＲアーチファクトの除去に優れた性能を示す。これにより、本発明のいくつかの実施例は、従来技術で提案された適応フィルタの構造よりはるかに簡単かつ迅速なものになる。
【００１３】
本発明のいくつかの実施例では、基準信号の前処理と自動リセット機能の特徴によって、比較的大きなステップサイズを適応に使用することが可能であるため、収束がより速くより安定して実行される。
【００１４】
本発明のいくつかの実施例は、計算上のコストを低減しＣＰＲアーチファクトの除去に優れた性能を発揮する。
本発明のその他の特徴および利点は、本発明の詳細な説明、図面および特許請求の範囲に記載されている。
【図面の簡単な説明】
【００１５】
【図１】本発明の一実施例のブロック図。
【図２】正常洞調律に応答するＥＣＧ信号、ＣＰＲ基準信号、および適応フィルタの出力プロットを示す。
【図３】心室細動中におけるＥＣＧ信号、ＣＰＲ基準信号、および適応フィルタの出力プロットを示す。
【図４】フィルタＸ平均最小二乗（ＦＸＬＭＳＡＮＣ）アルゴリズムのブロック図。
【図５】図４のアルゴリズムを使った実施例のブロック図。
【図６】図５の実施例に関係する２つのスペクトル出力の分布を示す。
【発明を実施するための形態】
【００１６】
本発明の可能な実施例の数はあまりにも多く、そのすべてを本明細書に記載することが
できない。そこで、以下には現時点で好ましい若干の可能な実施例を記載する。しかし、これらは、本発明の実施例の記述であって、この章に記載した詳細な実施例に限定するものではなく、特許請求の範囲に広義の用語で記載した本発明を記載したものであることを強調しておく。
【００１７】
一つの可能な実施例を図１のフローチャートに示す。ＡＥＤのフロントエンドはＥＣＧ信号と、胸部圧迫速度に相当するＣＰＲ信号との両方を取得する。胸部の速度の代わりに胸部の変位や加速度が測定される場合には、測定信号から一以上の積分または微分操作を経て数学的に速度を求めてもよい。
【００１８】
速度信号は前処理されてから適応フィルタに供給される。好ましい実施例では、前処理は、適応フィルタに供給される信号を０と１の範囲内に限定するように測定信号を正規化することである。しかし、正規化を行わなくてもよい。別の実施例では、相互相関などの方法によってＥＣＧ信号と基準信号に対する時間整合処理が行われる。これによって、圧迫に関して２つの信号が整合され、その結果、適応フィルタの入力信号はよりよく整合される。しかしこの整合処理を行わなくてもよい。速度信号に別の前処理を行って適応フィルタの性能を上げることもできる。
【００１９】
図１において、ｘ（ｎ）およびｙ（ｎ）は適応フィルタの入力および出力であり、適応フィルタとしてはＦＩＲフィルタ、ＩＩＲフィルタ、その他のタイプのフィルタが使用できる。好ましい実施例では、フィルタの係数は、推定ＥＣＧ信号によって動的制御される。
ｈ（ｎ）＝ｈ（ｎ−１）＋ｍ×ｅ（ｎ）×Ｘ（ｎ）
式中、ｈ（ｎ）はフィルタ係数を含むベクトル、ｍは各フィルタのステップサイズを含むベクトル、ｃ（ｎ）は推定ＥＣＧ信号、そしてＸ（ｎ）は入力データを含むベクトルである。推定ＥＣＧは、（アーチファクトを含む）実測ＥＣＧ信号からフィルタ出力ｙ（ｎ）を減算して得られる。
【００２０】
いくつかの実施例では、自動化リセット機構が設けられる。フィルタ出力ｙ（ｎ）と測定ＥＣＧ  ｓ（ｎ）の間の差が閾値を超えると、適応フィルタが、その係数をリセットするため、系が不安定になることは避けられる。
【００２１】
図１に示した以外のフィルタ構造や、フィルタリングの別の数学的表現が可能である。
図２は、正常洞調律に応答する図１に示した適応フィルタの性能を示す例である。（ａ）の信号は、ＣＰＲアーチファクトを含むＥＣＧ信号である。（ｂ）の信号は基準信号として使用した圧迫速度である。（ｃ）の信号は適応フィルタの出力である。
【００２２】
図３は、心室細動中における図１に示した適応フィルタの性能を示す例である。（ａ）の信号は、ＣＰＲアーチファクトを含むＥＣＧ信号である。（ｂ）の信号は基準信号として使用した圧迫速度である。（ｃ）の信号は適応フィルタの出力である。
【００２３】
図２と図３の両方に示すように、図１の実施例は、実測ＥＣＧ信号（ａ）に混入したＣＰＲアーチファクトを低減することができる。推定ＥＣＧ信号（ｃ）では、ＣＰＲアーチファクトは完全ではないがほとんど除去されている。基準信号として使用した速度信号（ｂ）は、実測ＥＣＧ信号（ａ）に混入したＣＰＲアーチファクトと明らかに相関性を有している。
【００２４】
適応フィルタは、信号中のアーチファクトが基準信号と相関性を有するが、所要の信号（推定ＥＣＧ）とは相関性を有さないことを前提としており、それゆえ、基準信号を使ってアーチファクトを順応して推定し、実測ＥＣＧ信号からその推定したアーチファクトを
減算する。
【００２５】
図２に示す結果は、０次ＦＩＲフィルタに基づいたものであり、そのＦＩＲフィルタは単にＥＣＧ信号の電流サンプルのスケールを適応的に変更するものである。ＣＰＲアーチファクトは、完全ではないにしても、かなり除去された。このように、この実施例は簡便さと性能の効率とを両立している。
【００２６】
適応フィルタを使用した場合、適応収束の速度は、通常、ステップサイズ変数によって制御される。速い収束は大きなステップサイズを必要とし、それはフィルタの安定性を損なわせる傾向がある。いくつかの実施例の自動リセット機構は、ステップサイズを動的に変更することができ、ひいては収束と安定性の関係を改善することができる。
【００２７】
フィルタ係数はサンプルごとに更新される。係数の変化すなわちｈ（ｎ）−ｈ（ｎ−１）は、ステップサイズと基準信号の積に比例する。それゆえ、基準信号の振幅はフィルタの安定性と収束性に影響し得る。したがって、基準信号の前処理は、基準信号の調節によってフィルタの性能を高めることができる。
【００２８】
別の実施例では、相互相関法のような方法によってＥＣＧ信号と速度信号に時間整合処理が行われる。この処理によって、前記２つの信号を圧迫に対して整合させることができる。好ましくは、その後、例えば、ＥＣＧと速度との間の平均二乗誤差を最小にする方法に含まれる適応フィルタリング法が使用される。
【００２９】
圧迫が加えられているときを検出し、自動的に適応フィルタを作動させるための処理ユニットを設けることができる。適応フィルタの出力（すなわち、アーチファクトの低減されたＥＣＧ信号）は、自動体外式除細動器（ＡＥＤ）の心室細動（ＶＦ）検出アルゴリズム（例えば、ショックアドバイザリーアルゴリズム）に供給することができる。
【００３０】
適応フィルタのＥＣＧ入力とＥＣＧ出力の間の差を表すエラー信号を生成することができる。このエラー信号は信号中に含まれるＣＰＲアーチファクトの量を表す測度を与えることができ、ＥＣＧの以後の処理を修正する手段として有用である。例えば、アーチファクトのレベルがかなり高く（例えば第一の閾値より高く）なった場合、ＶＦ検出アルゴリズムの閾値を高くして、ＥＣＧ信号に依然として残留する任意のＣＰＲアーチファクトに対する耐性を高めてもよい。そのレベルがさらに高く（例えば第二の閾値より高く）なると、ＶＦ検出を完全に中止してもよい。
【００３１】
好ましい実施例では、フィルタ出力は、除細動器や心電図記録機能を備えたその他の医療装置の画面にグラフ表示される。フィルタ出力は、医療機器の帯状記録紙レコーダに印刷されてもよい。あるいは、フィルタ出力は、処理手段によって行われる以後の信号処理のための入力信号として使用されてもよい。このような信号処理の目的は、ＱＳＲ検出、ペース調整中におけるペース調整された心拍の検出、不整脈分析、および心室細動や他のショックの効くリズムの検出である。
【００３２】
エラー信号に対してスペクトル解析を行ってもよい。エラー信号の周波数成分の主要バンドに基づいてＥＣＧ信号の予備フィルタリングをＶＦ検出に先だって調整することができる。例えば、エラー信号が主として３〜５Ｈｚ帯に見られる場合、ＶＦ検出（または他のＥＣＧ処理）アルゴリズムへの入力に先だって、追加のフィルタリングをこの帯域に設けることができる。
【００３３】
上記以外に他の多くの本発明の実施例が特許請求の範囲によって定義される本発明の範囲に含まれる。
例えば、同期化および速度波形における位相エラーの両者を改善するには、適応チャンネル等化法を採用することができる。胸部圧迫による救命士の救命成績は時間とともに変化し、非定常処理としてモデル化する方がよい場合は、フィルタ性能を高めるために、カルマンフィルタ法を使用することもできる。
【００３４】
ＥＣＧ信号と速度信号の時間整合も、当業者によく知られている相互相関法のような方法によって実施可能である。これは圧迫に関する両信号の整合を実現する。好ましくは、適応フィルタリング法、例えばＥＣＧと速度との間の平均二乗誤差を最小化する場合に使用される方法が使われる。
【００３５】
さらに別の実施例では、ＣＰＲ胸部圧迫に起因して誘発されるＢＣＧアーチファクトをできるだけ抑制するために、より精密な信号処理法を使うことができる。例えば、フィードフォワード能動型雑音消去（ＦＡＮＣ）の名で知られる方法を使用できる。図４は、ＷｉｄｒｏｗとＢｕｒｇｅｓｓによって開発されたフィルタＸ平均最小二乗法（ＦＸＬＭＳ
ＡＮＣ）アルゴリズムのブロック図である。Ｐ（ｚ）は未知のプラントを表し、そこを通過することで信号ｘ（ｎ）はフィルタ処理される。デジタルフィルタＷ（ｚ）は、エラー信号ｅ（ｎ）をできるだけ小さくするように適応的に調整される。一実施例では、図５に示すように、ｘ（ｎ）はフィルタ処理前のＥＣＧ信号を表し、Ｐ（ｚ）は図から省略してあり、ｄ（ｎ）は圧迫速度信号ｖ（ｎ）で近似されている。ＬＭＳアルゴリズムでは、平均二乗コスト関数ξ（ｎ）＝Ｅ［ｅ２（ｎ）］としたとき、適応フィルタは、ステップサイズμで係数ベクトルを負の勾配方向に更新する最急降下アルゴリズムを使って、瞬間二乗誤差ξ（ｎ）＝ｅ２（ｎ）を最小化する。
ｗ（ｎ＋１）＝ｗ（ｎ）−μ／２＊Ｎ~ξ（ｎ）
式中、Ｎ~ξ（ｎ）は、ｎが−２ｖ（ｎ）ｅ（ｎ）に等しい時における平均二乗誤差（Ｍ
ＳＥ）勾配の瞬間推定値である。（上記のＮ~は、本願の国際公報パンフレット中のティ
ルデ付きＮを指す。）ＦＸＬＭＳ  ＡＮＣアルゴリズムの安定性と精度は、可変カットオフローパスフィルタＨ（ｚ）を追加し、胸部圧迫アーチファクトに関係しない、ＥＣＧ中の周波数成分を消去することで改善することができる。一般に、胸部圧迫アーチファクトのスペクトルエネルギーは、ＥＣＧのそれより圧倒的に低い。カットオフ周波数は、多くの場合約３Ｈｚで十分であるが、患者毎および胸部圧迫を行う救命士に応じて変化する。この難点を解消するため、ｖ（ｎ）と、ローパスフィルタに対する最適カットオフ周波数ｆＣを決めるカットオフ周波数推定（ＣＦＥ）法の入力に対してＦＦＴを行う。好ましい実施例では、その決定は、５Ｈｚを超えない、それ以下に波形エネルギーの８０％が存在する周波数の計算に基づいて行われるが、このパーセンテージは変更してもよく、追加の決定論理を採用してもよい。例えば、ｘ（ｎ）についてＦＦＴを計算してもよく、また、ＣＦＥ法に入力してもよい。圧迫アーチファクトの周波数スペクトルＸ（ｚ）の振幅ピークの振幅をまず初めに正規化し、その後、正規化した入力値Ｘ’（ｚ）から速度スペクトルＶ（ｚ）を減算することにより差スペクトルΔＸ’（ｚ）＝Ｘ’（ｚ）−Ｖ’（ｚ）が計算される。次に、スペクトルエネルギーの大部分（この実施例では９７％）がその中に入るように、Ｖ（ｚ）およびΔＸ’（ｚ）についての周波数すなわち図６にそれぞれｆＣＶおよびｆＣＸで示す周波数を決定する。次に、ｆＣＶとｆＣＸのいずれか低い方にｆＣを設定する。あるいはｆＣＶとｆＣＸの間の中間周波数にｆＣを設定することができる。
【００３６】
特に、センサーが、速度信号が取り出される運動センサーの近くの胸郭上にある場合、フィルタによって、他の生理学的信号、例えば当業者によく知られているインピーダンス心拍記録信号（ＩＣＧ）、インピーダンス呼吸曲線記録信号（ＩＰＧ）またはパルス酸素濃度信号の質を高めることもできる。インピーダンス呼吸曲線記録信号における圧迫アーチファクトの最小化は、前記方法のどれを使っても行うことができる。
【００３７】
適応フィルタは、適応フィルタ出力と基準信号との相互相関または適応フィルタ出力と
実測ＥＣＧ信号との相互相関を最小化するために使用できる。
【図１】
2014195750.tif 000003
【図２】
2014195750.tif 000004
【図３】
2014195750.tif 000005
【図４】
2014195750.tif 000006
【図５】
2014195750.tif 000007
【図６】
2014195750.tif 000008