JPH0751236A - フィルタ処理されたｅｃｇ信号のリンギングを除去する装置 - Google Patents
フィルタ処理されたｅｃｇ信号のリンギングを除去する装置Info
- Publication number
- JPH0751236A JPH0751236A JP6164265A JP16426594A JPH0751236A JP H0751236 A JPH0751236 A JP H0751236A JP 6164265 A JP6164265 A JP 6164265A JP 16426594 A JP16426594 A JP 16426594A JP H0751236 A JPH0751236 A JP H0751236A
- Authority
- JP
- Japan
- Prior art keywords
- filter
- signal
- filtering
- ringing
- iir
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/30—Input circuits therefor
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Analogue/Digital Conversion (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Filters That Use Time-Delay Elements (AREA)
- Reverberation, Karaoke And Other Acoustics (AREA)
Abstract
(57)【要約】
【目的】 適切なオールパスフィルタを有するＥＣＧ信
号のリンギングを除去する装置を提供すること。 【構成】 ＩＩＲフィルタと、ＩＩＲフィルタの２倍の
位相シフトと等しい位相シフトを持つオールパスフィル
タとを有し、オールパスフィルタの零点は、ＩＩＲフィ
ルタの極と鏡像関係位置にある零点に置き換えた構成で
ある。
号のリンギングを除去する装置を提供すること。 【構成】 ＩＩＲフィルタと、ＩＩＲフィルタの２倍の
位相シフトと等しい位相シフトを持つオールパスフィル
タとを有し、オールパスフィルタの零点は、ＩＩＲフィ
ルタの極と鏡像関係位置にある零点に置き換えた構成で
ある。
Description
【０００１】
【産業上の利用分野】本発明は、アナログ入力信号を一
連のデジタル値に変換するＡ／Ｄコンバータと、Ａ／Ｄ
コンバータの出力側に接続されたフィルタユニットとを
有するフィルタ処理されたＥＣＧ信号のリンギングを除
去する装置に関する
連のデジタル値に変換するＡ／Ｄコンバータと、Ａ／Ｄ
コンバータの出力側に接続されたフィルタユニットとを
有するフィルタ処理されたＥＣＧ信号のリンギングを除
去する装置に関する
【０００２】
【従来の技術】信号の種々のフィルタ処理がＥＣＧ信号
の解析に必要である。
の解析に必要である。
【０００３】患者の体、たとえば腕と足とに配設された
２つの電極の間で記録されるＥＣＧ信号は、ＥＣＧ信号
に加えて、ＥＣＧ信号よりかなり大きなＤＣ電位も含ん
でいる。これらのＤＣ電位、または低周波信号は、たと
えば患者の呼吸運動の結果としての生体と電極の相対運
動のため変化する。
２つの電極の間で記録されるＥＣＧ信号は、ＥＣＧ信号
に加えて、ＥＣＧ信号よりかなり大きなＤＣ電位も含ん
でいる。これらのＤＣ電位、または低周波信号は、たと
えば患者の呼吸運動の結果としての生体と電極の相対運
動のため変化する。
【０００４】ツーパス（Ｔｗｏ−ｐａｓｓ）フィルタ処
理は、ＥＣＧ信号からＤＣ信号と低周波信号とを除去す
る１つの方法である。直線位相を持つフィルタは、非直
線位相を持つフィルタを用いて、信号を順方向（ｆｏｒ
ｗａｒｄ）と後進方向（ｂａｃｋｗａｒｄ）の両方のフ
ィルタ処理をすることにより得られ、英国特許第１，５
５６，５１２号明細書を参照すると、アナログフィルタ
を用いたこのフィルタ処理が記載されている。そのよう
なフィルタを用いて、位相ひずみが除去され、ＥＣＧ信
号形態のひずみを小さくできるので、フィルタ処理では
直線位相が望ましい。しかし、直線位相を持つフィルタ
の欠点は、リンギングの問題を完全に解決できないこと
である。リンギングは残存し、リンギングの形態と持続
時間はフィルタの伝達関数に依存する。いわゆるＦＩＲ
フィルタも直線位相を持つフィルタを実現するのに用い
ることができるが、ＦＩＲフィルタが制限されたメモリ
容量や制限された計算容量しか持たない装置で実行され
るものであるとすると、フィルタ曲線の傾斜の急峻さを
犠牲にしなければならない。
理は、ＥＣＧ信号からＤＣ信号と低周波信号とを除去す
る１つの方法である。直線位相を持つフィルタは、非直
線位相を持つフィルタを用いて、信号を順方向（ｆｏｒ
ｗａｒｄ）と後進方向（ｂａｃｋｗａｒｄ）の両方のフ
ィルタ処理をすることにより得られ、英国特許第１，５
５６，５１２号明細書を参照すると、アナログフィルタ
を用いたこのフィルタ処理が記載されている。そのよう
なフィルタを用いて、位相ひずみが除去され、ＥＣＧ信
号形態のひずみを小さくできるので、フィルタ処理では
直線位相が望ましい。しかし、直線位相を持つフィルタ
の欠点は、リンギングの問題を完全に解決できないこと
である。リンギングは残存し、リンギングの形態と持続
時間はフィルタの伝達関数に依存する。いわゆるＦＩＲ
フィルタも直線位相を持つフィルタを実現するのに用い
ることができるが、ＦＩＲフィルタが制限されたメモリ
容量や制限された計算容量しか持たない装置で実行され
るものであるとすると、フィルタ曲線の傾斜の急峻さを
犠牲にしなければならない。
【０００５】Ｅｕｒｏｐｅａｎ Ｈｅａｒｔ Ｊｏｕｒ
ｎａｌ（１９９１）、１２、４７３−４８０ページを参
照すると、梗塞後の患者のＱＲＳ複合信号の後部でのマ
イクロボルトレベルの小振幅を有する高周波信号の存在
は、将来生命の脅威となる心室の不整脈のリスクが高ま
っている明瞭な徴候であることが示されている。これら
のレイトポテンシャルを識別し、解析するためにハイパ
スフィルタ処理が必要である。したがって、１つの問題
は、ＱＲＳ信号自体によって惹起され、レイトポテンシ
ャルも埋もれてしまうリンギングである。この問題を解
決するために種々の提案がなされている。ＬｉｎｋとＴ
ｒａｈｍｓの”Ｈｉｇｈｐａｓｓ Ｆｉｌｔｅｒｓ ｆ
ｏｒ Ｄｅｔｅｃｔｉｎｇ Ｌａｔｅ Ｐｏｔｅｎｔｉ
ａｌｓ”，Ｐｒｏｃ．Ｃｏｍｐｕｔｅｒｓ ｉｎ Ｃａ
ｒｄｉｏｌｏｇｙ １９９２，１５９−１６２ページに
おいて、レイトポテンシャルを検出するために非巡回単
調２項ハイパスフィルタを利用することが提案されてい
る。
ｎａｌ（１９９１）、１２、４７３−４８０ページを参
照すると、梗塞後の患者のＱＲＳ複合信号の後部でのマ
イクロボルトレベルの小振幅を有する高周波信号の存在
は、将来生命の脅威となる心室の不整脈のリスクが高ま
っている明瞭な徴候であることが示されている。これら
のレイトポテンシャルを識別し、解析するためにハイパ
スフィルタ処理が必要である。したがって、１つの問題
は、ＱＲＳ信号自体によって惹起され、レイトポテンシ
ャルも埋もれてしまうリンギングである。この問題を解
決するために種々の提案がなされている。ＬｉｎｋとＴ
ｒａｈｍｓの”Ｈｉｇｈｐａｓｓ Ｆｉｌｔｅｒｓ ｆ
ｏｒ Ｄｅｔｅｃｔｉｎｇ Ｌａｔｅ Ｐｏｔｅｎｔｉ
ａｌｓ”，Ｐｒｏｃ．Ｃｏｍｐｕｔｅｒｓ ｉｎ Ｃａ
ｒｄｉｏｌｏｇｙ １９９２，１５９−１６２ページに
おいて、レイトポテンシャルを検出するために非巡回単
調２項ハイパスフィルタを利用することが提案されてい
る。
【０００６】米国特許明細書第４，４２２，４５９号明
細書に、患者のＱＲＳ複合信号の後方の部分に高周波信
号を含むタイムセグメントの有無を解明し、このタイム
セグメントの長さを測定する技術が記載されている。測
定されたアナログのＥＣＧ信号は、デジタル信号に変換
され、正常なＱＲＳ信号は、比較的ノイズの少ない複合
信号を得るために、数百程度の拍動にわたって平均化さ
れる。平均化された複合信号は４次のバタワースフィル
タでハイパスフィルタ処理される。フィルタ処理は、Ｑ
ＲＳの開始＋４０msまで順方向で実行され、また終わり
から開始＋４０msの時点まで後進方向で実行される。２
つの信号は結合され、合成信号を形成し、その合成信号
において、リンギングは、ＱＲＳ信号自体の”誤りとな
る”信号を犠牲にして、ＱＲＳ信号の直前、直後で避け
られる。フィルタ処理された信号は、いわゆるレイトポ
テンシャルの有無を確認するために測定される。この技
術の欠点は、フィルタ処理の曲線は不連続であり、ＱＲ
Ｓ複合信号は、誤りとなるエネルギー内容（ｃｏｎｔｅ
ｎｔ）を有することである。
細書に、患者のＱＲＳ複合信号の後方の部分に高周波信
号を含むタイムセグメントの有無を解明し、このタイム
セグメントの長さを測定する技術が記載されている。測
定されたアナログのＥＣＧ信号は、デジタル信号に変換
され、正常なＱＲＳ信号は、比較的ノイズの少ない複合
信号を得るために、数百程度の拍動にわたって平均化さ
れる。平均化された複合信号は４次のバタワースフィル
タでハイパスフィルタ処理される。フィルタ処理は、Ｑ
ＲＳの開始＋４０msまで順方向で実行され、また終わり
から開始＋４０msの時点まで後進方向で実行される。２
つの信号は結合され、合成信号を形成し、その合成信号
において、リンギングは、ＱＲＳ信号自体の”誤りとな
る”信号を犠牲にして、ＱＲＳ信号の直前、直後で避け
られる。フィルタ処理された信号は、いわゆるレイトポ
テンシャルの有無を確認するために測定される。この技
術の欠点は、フィルタ処理の曲線は不連続であり、ＱＲ
Ｓ複合信号は、誤りとなるエネルギー内容（ｃｏｎｔｅ
ｎｔ）を有することである。
【０００７】米国特許第４，４５８，６９１号明細書
に、レイトポテンシャルを検出するためのＥＣＧ信号を
順方向でハイパスフィルタ処理するシステムが記載され
ており、ＱＲＳ複合信号の種々のセグメントのフィルタ
処理を選択する適応ハイパスフィルタが、リンギングの
問題を解決するために用いられている。この技術もいく
つかの欠点を有する。したがって、リンギングは完全に
除去されず、できる限りリンギングの影響を軽減するた
めに、直線位相を有するフィルタ長の短いＦＩＲフィル
タが用いられており、周波数分割特性が悪くなる。
に、レイトポテンシャルを検出するためのＥＣＧ信号を
順方向でハイパスフィルタ処理するシステムが記載され
ており、ＱＲＳ複合信号の種々のセグメントのフィルタ
処理を選択する適応ハイパスフィルタが、リンギングの
問題を解決するために用いられている。この技術もいく
つかの欠点を有する。したがって、リンギングは完全に
除去されず、できる限りリンギングの影響を軽減するた
めに、直線位相を有するフィルタ長の短いＦＩＲフィル
タが用いられており、周波数分割特性が悪くなる。
【０００８】米国特許第４，４５８，６９２号明細書
に、ＥＣＧ信号の順方向のフィルタ処理を介して、レイ
トポテンシャルを検出する別の方法が記載されており、
リンギングの問題を解決するために、入力信号の大きさ
に依存させてフィルタの利得を調節する方法が用いられ
ている。利得を制御する信号はＱＲＳ信号を入力信号と
して用いるフィルタからの出力信号にすることができ、
入力信号が所定のレベルを超えた場合、フィルタの利得
を制限すると、たとえば出力信号に上のスパイク状のＲ
波の影響、つまりスパイク状のＲ波の後で発生するリン
ギングに関する影響を制限することができる。この方法
も全てのリンギングを除去できず、急峻なフィルタ曲線
の傾斜を犠牲にしなければならないという欠点を有す
る。
に、ＥＣＧ信号の順方向のフィルタ処理を介して、レイ
トポテンシャルを検出する別の方法が記載されており、
リンギングの問題を解決するために、入力信号の大きさ
に依存させてフィルタの利得を調節する方法が用いられ
ている。利得を制御する信号はＱＲＳ信号を入力信号と
して用いるフィルタからの出力信号にすることができ、
入力信号が所定のレベルを超えた場合、フィルタの利得
を制限すると、たとえば出力信号に上のスパイク状のＲ
波の影響、つまりスパイク状のＲ波の後で発生するリン
ギングに関する影響を制限することができる。この方法
も全てのリンギングを除去できず、急峻なフィルタ曲線
の傾斜を犠牲にしなければならないという欠点を有す
る。
【０００９】また、米国特許第４，４９２，２３５号明
細書は、順方向のＱＲＳ信号のフィルタ処理を通しての
レイトポテンシャルの検出に関し、リンギングの問題を
解決するために適応ハイパスフィルタ処理が用いられて
いる。
細書は、順方向のＱＲＳ信号のフィルタ処理を通しての
レイトポテンシャルの検出に関し、リンギングの問題を
解決するために適応ハイパスフィルタ処理が用いられて
いる。
【００１０】米国特許第５，０２５，７９４号明細書
に、レイトポテンシャルを順方向でフィルタ処理する別
の方法が記載されており、その方法において、ＱＲＳ信
号のサンプル値は、順方向と逆方向の両方のフィルタ処
理をされ、フィルタ処理された信号は一緒に加算され、
合成信号は、レイトポテンシャルと相応する信号成分を
生成するために平滑化される。前記明細書に記載の方法
は複雑である。
に、レイトポテンシャルを順方向でフィルタ処理する別
の方法が記載されており、その方法において、ＱＲＳ信
号のサンプル値は、順方向と逆方向の両方のフィルタ処
理をされ、フィルタ処理された信号は一緒に加算され、
合成信号は、レイトポテンシャルと相応する信号成分を
生成するために平滑化される。前記明細書に記載の方法
は複雑である。
【００１１】
【発明が解決しようとする課題】本発明の課題は、先技
術における前記の欠点を除去し、フィルタ処理したＥＣ
Ｇ信号のリンギングを除去する装置を提供することであ
る。
術における前記の欠点を除去し、フィルタ処理したＥＣ
Ｇ信号のリンギングを除去する装置を提供することであ
る。
【００１２】本発明の別の課題は、本発明による装置に
用いるのに適切なオールパスフィルタを製作する方法を
示すことである。
用いるのに適切なオールパスフィルタを製作する方法を
示すことである。
【００１３】
【課題を解決するための手段】本発明の課題は、フィル
タユニットが、巡回ＩＩＲフィルタと、ＩＩＲフィルタ
の２倍の位相シフトと等しい位相シフトを持つ後続のオ
ールパスフィルタとを有し、オールパスフィルタ処理
は、ＩＩＲフィルタ処理と関連して時間を反転して実行
され、およびオールパスフィルタがＩＩＲフィルタの零
点を、Ｚ平面の単位円でＩＩＲフィルタの極と鏡像関係
位置にある新しい零点に置き換えられた構成を有するこ
とにより達成される。
タユニットが、巡回ＩＩＲフィルタと、ＩＩＲフィルタ
の２倍の位相シフトと等しい位相シフトを持つ後続のオ
ールパスフィルタとを有し、オールパスフィルタ処理
は、ＩＩＲフィルタ処理と関連して時間を反転して実行
され、およびオールパスフィルタがＩＩＲフィルタの零
点を、Ｚ平面の単位円でＩＩＲフィルタの極と鏡像関係
位置にある新しい零点に置き換えられた構成を有するこ
とにより達成される。
【００１４】本発明による装置は、除去されなければな
らないリンギングを惹起するオリジナルフィルタの位相
シフトの２倍に正確に相当する位相シフトを持つオール
パスフィルタを有する。オールパスフィルタは、オリジ
ナルフィルタ処理の方向に対して時間を反転して、解析
すべき信号の一部に適用される。このように、問題とな
る装置に適用するために、非常に有利な方法でリンギン
グは除去される。したがって、本発明による装置によ
り、信号を巡回ハイパスフィルタを用いて順方向で、フ
ィルタの急峻さに何ら制限を加えることなくフィルタ処
理した後で、レイトポテンシャルを簡単に検出すること
ができる。また、本発明による装置は、ＱＲＳ複合信号
の波形の境界を定め、形を分類し、ＳＴレベルを求める
のに有利にも用いることができる。
らないリンギングを惹起するオリジナルフィルタの位相
シフトの２倍に正確に相当する位相シフトを持つオール
パスフィルタを有する。オールパスフィルタは、オリジ
ナルフィルタ処理の方向に対して時間を反転して、解析
すべき信号の一部に適用される。このように、問題とな
る装置に適用するために、非常に有利な方法でリンギン
グは除去される。したがって、本発明による装置によ
り、信号を巡回ハイパスフィルタを用いて順方向で、フ
ィルタの急峻さに何ら制限を加えることなくフィルタ処
理した後で、レイトポテンシャルを簡単に検出すること
ができる。また、本発明による装置は、ＱＲＳ複合信号
の波形の境界を定め、形を分類し、ＳＴレベルを求める
のに有利にも用いることができる。
【００１５】巡回フィルタでのフィルタ処理と、後続の
時間を反転したオールパスフィルタ処理の結果は、直線
位相を有するフィルタ処理ではないが、巡回フィルタの
位相応答を打ち消す位相応答を有するフィルタ処理であ
ることに注目すべきである。これは、問題となる用途に
おいて、リンギングを完全に除去するのに最適な特性で
ある。したがって、本発明による装置の全体の位相シフ
トは、巡回ＩＩＲフィルタで信号を後進方向にフィルタ
処理する場合に起こる位相シフトと同一である。このよ
うに、オリジナルフィルタ処理が従来の方法で順方向に
実行されるのであるが、後進方向のフィルタ処理と同一
の原理により、リンギングは重要な領域から”シフト”
される。
時間を反転したオールパスフィルタ処理の結果は、直線
位相を有するフィルタ処理ではないが、巡回フィルタの
位相応答を打ち消す位相応答を有するフィルタ処理であ
ることに注目すべきである。これは、問題となる用途に
おいて、リンギングを完全に除去するのに最適な特性で
ある。したがって、本発明による装置の全体の位相シフ
トは、巡回ＩＩＲフィルタで信号を後進方向にフィルタ
処理する場合に起こる位相シフトと同一である。このよ
うに、オリジナルフィルタ処理が従来の方法で順方向に
実行されるのであるが、後進方向のフィルタ処理と同一
の原理により、リンギングは重要な領域から”シフト”
される。
【００１６】
【実施例】次に、本発明による１つの実施例を本発明に
よるＥＣＧ信号を記録する装置のブロック図である添付
図面を参照して詳細に説明する。
よるＥＣＧ信号を記録する装置のブロック図である添付
図面を参照して詳細に説明する。
【００１７】ＥＣＧは、患者２の体表面、たとえば腕４
と足６とに配設される電極の間の電位として記録され
る。ＥＣＧ記録は、アナログ信号の調整と同様にユニッ
ト８で行われる。また、このユニット８において、ＥＣ
Ｇ信号に含まれる大きなＤＣ成分は、基線の周囲でＥＣ
Ｇ信号自体を測定できるように除去される。このように
して、後続のＡ／Ｄコンバータ１０に必要なダイナミッ
クレンジと分解能を低下することができる。
と足６とに配設される電極の間の電位として記録され
る。ＥＣＧ記録は、アナログ信号の調整と同様にユニッ
ト８で行われる。また、このユニット８において、ＥＣ
Ｇ信号に含まれる大きなＤＣ成分は、基線の周囲でＥＣ
Ｇ信号自体を測定できるように除去される。このように
して、後続のＡ／Ｄコンバータ１０に必要なダイナミッ
クレンジと分解能を低下することができる。
【００１８】Ａ／Ｄコンバータ１０の出力側は、多数の
ＱＲＳ複合信号にわたって平均信号を形成するように平
均化装置１２に接続され、前記平均信号は、バタワース
フィルタ等の巡回ＩＩＲフィルタ１４で順方向でフィル
タ処理される。
ＱＲＳ複合信号にわたって平均信号を形成するように平
均化装置１２に接続され、前記平均信号は、バタワース
フィルタ等の巡回ＩＩＲフィルタ１４で順方向でフィル
タ処理される。
【００１９】フィルタ１４の後に、平均化された複合信
号の後方の部分を、時間を反転してフィルタ処理する位
相シフト形のオールパスフィルタ１６がある。ＩＩＲフ
ィルタ１４で生ずるリンギングは、オールパスフィルタ
１６の出力信号から除去されるので、出力信号のより正
確な測定と解析ができる。
号の後方の部分を、時間を反転してフィルタ処理する位
相シフト形のオールパスフィルタ１６がある。ＩＩＲフ
ィルタ１４で生ずるリンギングは、オールパスフィルタ
１６の出力信号から除去されるので、出力信号のより正
確な測定と解析ができる。
【００２０】バタワースフィルタ１４の伝達関数は次の
とおり書くことができる。
とおり書くことができる。
【００２１】
【数１】
【００２２】オールパスフィルタ１６は、バタワースフ
ィルタ１４の零点を、単位円でフィルタの極と鏡像関係
に位置する新しい零点と置き換えることにより構成され
る。したがって、オールパスフィルタ１６の伝達関数は
次のようになる。
ィルタ１４の零点を、単位円でフィルタの極と鏡像関係
に位置する新しい零点と置き換えることにより構成され
る。したがって、オールパスフィルタ１６の伝達関数は
次のようになる。
【００２３】
【数２】
【００２４】オールパスフィルタ１６は、零点を、Ｚ平
面の単位円で極を鏡像関係位置にすることにより得られ
る補償しようとするフィルタと同一の極を有するバタワ
ースフィルタのものになる。
面の単位円で極を鏡像関係位置にすることにより得られ
る補償しようとするフィルタと同一の極を有するバタワ
ースフィルタのものになる。
【００２５】したがって、オールパスフィルタ１６の位
相シフトは次のようになる。
相シフトは次のようになる。
【００２６】ΦA（ω） ＝−４ω＋ΦD（ω） −Φ
D（−ω） ＝２ΦD（ω） −４ω ΦD（ω） ＝１つの零点を有するフィルタの位相シフト
である。すなわち、
D（−ω） ＝２ΦD（ω） −４ω ΦD（ω） ＝１つの零点を有するフィルタの位相シフト
である。すなわち、
【００２７】
【数３】
【００２８】であり、 ΦD（−ω） ＝−ΦD（ω） バタワースフィルタの位相シフトは、 ΦB（−ω） ＝−２ω＋ΦD（ω） したがって、本発明による装置のオールパスフィルタで
カスケード接続されるバタワースフィルタの全体の位相
シフトは、オールパスフィルタ処理が逆方向で行われる
ので、次のようになる。
カスケード接続されるバタワースフィルタの全体の位相
シフトは、オールパスフィルタ処理が逆方向で行われる
ので、次のようになる。
【００２９】Φtotal（ω） ＝ΦB（ω） −ΦA（ω）
＝−ΦD（ω） ＋２ω＝−ΦB（ω） 本発明による装置の後の位相シフトは、バタワースフィ
ルタで信号を後進方向にフィルタ処理する場合に得られ
る位相シフトと同一である。したがって、オリジナルフ
ィルタ処理は従来の方法で順方向に行われたが、リンギ
ングは重要な領域から”変移”される。
＝−ΦD（ω） ＋２ω＝−ΦB（ω） 本発明による装置の後の位相シフトは、バタワースフィ
ルタで信号を後進方向にフィルタ処理する場合に得られ
る位相シフトと同一である。したがって、オリジナルフ
ィルタ処理は従来の方法で順方向に行われたが、リンギ
ングは重要な領域から”変移”される。
【００３０】オールパスフィルタを製作する方法は、除
去されるべきリンギングを生成するバタワースフィルタ
から出発し、バタワースフィルタの零点は、Ｚ平面にお
ける単位円でフィルタの極と鏡像関係位置にある新しい
零点に置き換えられる。ＩＩＲフィルタの全ての零点は
Ｚ平面の単位円に位置するものと仮定する。このオール
パスフィルタに対する位相シフトは次のようになる。
去されるべきリンギングを生成するバタワースフィルタ
から出発し、バタワースフィルタの零点は、Ｚ平面にお
ける単位円でフィルタの極と鏡像関係位置にある新しい
零点に置き換えられる。ＩＩＲフィルタの全ての零点は
Ｚ平面の単位円に位置するものと仮定する。このオール
パスフィルタに対する位相シフトは次のようになる。
【００３１】 ΦA（ω） ＝２＊ΦH（ω） ＋（ｎH −ｎA ）ω ΦH（ω） ＝バタワースフィルタの位相シフトであり、 ｎH ＝バタワースフィルタの零点の数であり、 ｎA ＝オールパスフィルタの零点の数である。
【００３２】バタワースフィルタは極と同数の零点を有
するので、第２項はゼロとなる。他方で、補償は、オー
ルパスフィルタ処理の後で、一定の遅延でなされなけれ
ばならない。バタワースフィルタでのフィルタ処理と時
間を反転したのオールパスフィルタ処理との結果は、直
線位相を持つフィルタ処理ではないが、バタワースフィ
ルタの位相応答を打ち消す位相応答を持つフィルタ処理
であり、問題となる用途にとって最適のリンギングの除
去を達成することに注目すべきである。
するので、第２項はゼロとなる。他方で、補償は、オー
ルパスフィルタ処理の後で、一定の遅延でなされなけれ
ばならない。バタワースフィルタでのフィルタ処理と時
間を反転したのオールパスフィルタ処理との結果は、直
線位相を持つフィルタ処理ではないが、バタワースフィ
ルタの位相応答を打ち消す位相応答を持つフィルタ処理
であり、問題となる用途にとって最適のリンギングの除
去を達成することに注目すべきである。
【００３３】
【発明の効果】本発明により、ＥＣＧ信号のリンギング
を除去する装置を提供し、本発明に用いるのに適切なオ
ールパスフィルタを製作することができる。
を除去する装置を提供し、本発明に用いるのに適切なオ
ールパスフィルタを製作することができる。
【図１】本発明によるＥＣＧ信号を記録する装置のブロ
ック図である。
ック図である。
１４ ＩＩＲフィルタ １６ オールパスフィルタ
Claims (5)
- 【請求項１】 アナログ入力信号を一連のデジタル値に
変換するＡ／Ｄコンバータと、Ａ／Ｄコンバータの出力
側に接続されるフィルタユニットとを有するフィルタ処
理されたＥＣＧ信号のリンギングを除去する装置におい
て、フィルタユニットは、巡回ＩＩＲフィルタと、ＩＩ
Ｒフィルタの２倍の位相シフトと等しい位相シフトを持
つ後続のオールパスフィルタとを有し、オールパスフィ
ルタ処理は、ＩＩＲフィルタ処理と関連して時間を反転
して実行されることを特徴とするフィルタ処理されたＥ
ＣＧ信号のリンギングを除去する装置。 - 【請求項２】 平均化装置は、フィルタ処理の前に多数
のＱＲＳ複合信号に対して信号平均値を形成するため
に、Ａ／Ｄコンバータとフィルタユニットとの間に配設
されていることを特徴とする請求項１記載の装置。 - 【請求項３】 フィルタユニットは、ＱＲＳ複合信号の
後方の部分だけをオールパスフィルタでフィルタ処理す
るように構成されていることを特徴とする請求項１また
は２記載の装置。 - 【請求項４】 フィルタユニットは、ＱＲＳ複合信号の
前方の部分だけをオールパスフィルタでフィルタ処理す
るように構成されていることを特徴とする請求項１また
は２記載の装置。 - 【請求項５】 オールパスフィルタが、ＩＩＲフィルタ
の零点を、Ｚ平面の単位円でＩＩＲフィルタの極と鏡像
関係位置にある新しい零点に置き換えられた構成を有す
ることを特徴とする請求項１記載の装置。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| SE19939302433A SE9302433D0 (sv) | 1993-07-16 | 1993-07-16 | Anordning foer eliminering av ringningar hos filtrerade ekg-signaler |
| SE9302433-9 | 1993-07-16 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH0751236A true JPH0751236A (ja) | 1995-02-28 |
Family
ID=20390623
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP6164265A Pending JPH0751236A (ja) | 1993-07-16 | 1994-07-15 | フィルタ処理されたｅｃｇ信号のリンギングを除去する装置 |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US5532951A (ja) |
| EP (1) | EP0634137B1 (ja) |
| JP (1) | JPH0751236A (ja) |
| AT (1) | ATE148828T1 (ja) |
| DE (1) | DE69401727T2 (ja) |
| SE (1) | SE9302433D0 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109567788A (zh) * | 2018-11-29 | 2019-04-05 | 武汉中旗生物医疗电子有限公司 | 一种去除振铃的心电信号滤波方法 |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6249696B1 (en) | 1999-01-15 | 2001-06-19 | Medtronic Physio-Control Manufacturing Corp. | Method and apparatus for increasing the low frequency dynamic range of a digital ECG measuring system |
| US6280391B1 (en) | 1999-02-08 | 2001-08-28 | Physio-Control Manufacturing Corporation | Method and apparatus for removing baseline wander from an egg signal |
| US20070078353A1 (en) * | 2005-10-04 | 2007-04-05 | Welch Allyn, Inc. | Method and apparatus for removing baseline wander from an ECG signal |
| EP1909227B1 (en) * | 2006-10-03 | 2018-06-27 | Vestel Elektronik Sanayi ve Ticaret A.S. | Method of and apparatus for minimizing ringing artifacts in an input image |
Family Cites Families (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5271960A (en) * | 1975-12-11 | 1977-06-15 | Fukuda Denshi Kk | Realltime nonnphase filter using delay circuit |
| US4422459A (en) * | 1980-11-18 | 1983-12-27 | University Patents, Inc. | Electrocardiographic means and method for detecting potential ventricular tachycardia |
| US4458692A (en) * | 1982-02-11 | 1984-07-10 | Arrhythmia Research Technology, Inc. | System and method for predicting ventricular tachycardia with a gain controlled high pass filter |
| US4492235A (en) * | 1982-02-11 | 1985-01-08 | Arrhythmia Research Technology, Inc. | System and method for predicting ventricular tachycardia by derivative analysis |
| US4458691A (en) * | 1982-02-11 | 1984-07-10 | Arrhythmia Research Technology, Inc. | System and method for predicting ventricular tachycardia by adaptive high pass filter |
| US5025794A (en) * | 1988-08-30 | 1991-06-25 | Corazonix Corporation | Method for analysis of electrocardiographic signal QRS complex |
| DE3922469A1 (de) * | 1989-07-07 | 1991-01-17 | Nixdorf Computer Ag | Verfahren zum filtern digitalisierter signale |
| EP0407647B1 (de) * | 1989-07-14 | 1993-10-06 | Haberl, Ralph, Dr.med. | Einrichtung zur Bewertung ausgewählter Signalanteile in physiologischen Messsignalen, insbesondere von Spätpotentialen in Elektrokardiogrammen |
| US4939576A (en) * | 1989-07-24 | 1990-07-03 | Campbell Jack J | Adaptive ringing reducer for television signal processing |
| US5139027A (en) * | 1990-12-03 | 1992-08-18 | Cinventa Aktiebolag | Method of filtering an analog ecg signal |
| US5341811A (en) * | 1991-03-26 | 1994-08-30 | Allegheny-Singer Research Institute | Method and apparatus for observation of ventricular late potentials |
| US5269313A (en) * | 1991-09-09 | 1993-12-14 | Sherwood Medical Company | Filter and method for filtering baseline wander |
| US5297557A (en) * | 1992-10-14 | 1994-03-29 | Del Mar Avionics | Stress test system with bidirectional filter |
-
1993
- 1993-07-16 SE SE19939302433A patent/SE9302433D0/xx unknown
-
1994
- 1994-06-22 AT AT94109663T patent/ATE148828T1/de not_active IP Right Cessation
- 1994-06-22 DE DE69401727T patent/DE69401727T2/de not_active Expired - Fee Related
- 1994-06-22 EP EP94109663A patent/EP0634137B1/en not_active Expired - Lifetime
- 1994-06-27 US US08/266,233 patent/US5532951A/en not_active Expired - Fee Related
- 1994-07-15 JP JP6164265A patent/JPH0751236A/ja active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109567788A (zh) * | 2018-11-29 | 2019-04-05 | 武汉中旗生物医疗电子有限公司 | 一种去除振铃的心电信号滤波方法 |
| CN109567788B (zh) * | 2018-11-29 | 2021-08-20 | 武汉中旗生物医疗电子有限公司 | 一种去除振铃的心电信号滤波方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| SE9302433D0 (sv) | 1993-07-16 |
| EP0634137A1 (en) | 1995-01-18 |
| US5532951A (en) | 1996-07-02 |
| ATE148828T1 (de) | 1997-02-15 |
| DE69401727D1 (de) | 1997-03-27 |
| EP0634137B1 (en) | 1997-02-12 |
| DE69401727T2 (de) | 1997-08-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5772603A (en) | Device for filtering ECG signals | |
| US5297557A (en) | Stress test system with bidirectional filter | |
| US6280391B1 (en) | Method and apparatus for removing baseline wander from an egg signal | |
| US5025784A (en) | Apparatus and method for detecting and processing impedance rheogram | |
| US5318036A (en) | Method and apparatus for removing baseline wander from an ECG signal | |
| Shusterman et al. | Enhancing the precision of ECG baseline correction: selective filtering and removal of residual error | |
| EP1424637A1 (en) | Artifact removal from an electric signal | |
| US5908393A (en) | Reducing noise in a biological signal | |
| JP3171256B2 (ja) | 生体の測定信号から基線変動を除去する装置 | |
| Dotsinsky et al. | Optimization of bi-directional digital filtering for drift suppression in electrocardiogram signals | |
| JPH0852117A (ja) | 患者モニタ | |
| JPS63226338A (ja) | Ｅｃｇ波信号からペース・パルス信号を除去する方法とｅｃｇ波信号からペース・パルス信号を除去する装置 | |
| US6216031B1 (en) | Apparatus for enhancing signals in ECGs artefacts | |
| JPH0322770B2 (ja) | ||
| JPH0751236A (ja) | フィルタ処理されたｅｃｇ信号のリンギングを除去する装置 | |
| JP3136430B2 (ja) | 生体信号測定装置 | |
| JPH06197877A (ja) | 生体信号採取装置 | |
| Dotsinsky et al. | Tremor suppression in ECG | |
| JP2784046B2 (ja) | 心電図ｒ波検出装置 | |
| Philips | Adaptive base line correction of EEGs using warped polynomials | |
| JPS58223927A (ja) | デ−タ圧縮処理装置 | |
| JPS61206428A (ja) | 生体信号ア−チフアクト検出除去方式 | |
| EP3893723B1 (en) | Filtering unit for electrocardiography applications | |
| Dotsinsky et al. | Tremor suppression in the electrocardiogram | |
| Arzbaecher et al. | Computer Equalization of Low-Frequency Components in Tape Recorded Electrocardiographic Signals |