JPH0751237A - Ecg信号の信号対雑音比を高める方法および装置 - Google Patents

Ecg信号の信号対雑音比を高める方法および装置

Info

Publication number
JPH0751237A
JPH0751237A JP6164269A JP16426994A JPH0751237A JP H0751237 A JPH0751237 A JP H0751237A JP 6164269 A JP6164269 A JP 6164269A JP 16426994 A JP16426994 A JP 16426994A JP H0751237 A JPH0751237 A JP H0751237A
Authority
JP
Japan
Prior art keywords
signals
signal
unit
correlation
time
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
Application number
JP6164269A
Other languages
English (en)
Inventor
Leif Soernmo
セルンモ ライフ
Thomas Ohlsson
オールソン トーマス
Roozbeh Atarius
アタリウス ロースベー
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens Elema AB
Original Assignee
Siemens Elema AB
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Siemens Elema AB filed Critical Siemens Elema AB
Publication of JPH0751237A publication Critical patent/JPH0751237A/ja
Pending legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/30Input circuits therefor
    • A61B5/307Input circuits therefor specially adapted for particular uses
    • A61B5/308Input circuits therefor specially adapted for particular uses for electrocardiography [ECG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/30Input circuits therefor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7246Details of waveform analysis using correlation, e.g. template matching or determination of similarity

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Signal Processing (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Cardiology (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Interface Circuits In Exchanges (AREA)
  • Stereo-Broadcasting Methods (AREA)
  • Noise Elimination (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Transforming Light Signals Into Electric Signals (AREA)

Abstract

(57)【要約】 【目的】 信号対雑音比を向上させ、小振幅信号成分を
推定するECG信号処理方法と装置を提供することであ
る。 【構成】 同期ユニットが、信号の集合を形成するため
に信号を時間同期させ、ウィンドウユニットに制御され
る計算ユニットが集合内の信号の間の一致を求め、最適
フィルタの決定に一致が用いられる構成である。

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は、ECG信号を記録し、
A/D変換し、処理する手段を有するECG信号の信号
対雑音比を高める方法および装置に関する。
【0002】
【従来の技術】生体表面で記録されたECG信号は、通
常筋肉の活動、記録されたアナログECG信号のデジタ
ル化、記録装置、および電源からの作用に起因する広帯
域のノイズを有する。しかし、現代の技術を用いること
により、記録装置と電源のノイズ源の作用をほとんど排
除することができる。しかし、患者の筋肉の活動とデジ
タル処理から生じるノイズは、しばしばECG信号内の
小振幅信号成分を検出できなくしてしまう。
【0003】従来ノイズを低減し、これらの小振幅信号
成分に対する信号対雑音比を向上させるために、ECG
信号の種々の平均化処理をしてきた。WeerdとMa
rtensの論文”Theory and Pract
ice of A Posteriori Wiene
r Filtering of Average Ev
oked Potentials”,Biol. Cy
bern.,30,ページ81−94,(1978)、
及びWeerdとKapの論文”A Posterio
ri Time Varying Filtering
of Averaged Evoked Poten
tials”,Biol. Cybern.,41,ペ
ージ223−234,(1981)を参照すると、時間
で変化するウィーナーフィルタ処理も、信号対雑音比を
高めるために用いられている。この技術は、種々の周波
数に対して求められた推定信号とノイズの作用に基づい
ている。この処理では、使用されるフーリエ変換処理に
おいて、いくらかのスペクトルの漏れが不可避である。
また、低い信号対雑音比では、この方法ではノイズ内の
信号を誤って推定することを避けることはできない。
【0004】
【発明が解決しようとする課題】本発明の課題は、全て
の前記欠点を有するフーリエ変換を用いずに、信号対雑
音比の向上を可能にし、次にいわゆる”レイトポテンシ
ャル”等のECG信号内の小振幅信号成分の推定を可能
にする新規の方法を提案し、新規の装置を提供すること
である。
【0005】
【課題を解決するための手段】本発明の課題は、複数の
心周期からのA/D変換した信号を、信号の集合を形成
するために時間同期し、集合内の信号の間の一致を、心
周期の少なくとも一部分の期間中に少なくとも2つのサ
ンプルを有する時間インターバルにおいて求める方法、
および同期ユニットが、信号の集合を形成するために、
複数の心周期からの信号を時間同期させるように構成さ
れ、ウィンドウユニットにより制御される計算ユニット
が、心周期の少なくとも一部分の期間中に少なくとも2
つのサンプルを有する時間インターバルにおいて、集合
内の信号の間の一致を求めるように構成されている装置
を用いて達成される。
【0006】本発明において、複数の心周期からのEC
G信号の再現部分は、有効な信号対雑音比と関連する特
性のフィルタでろ波して取り出される。
【0007】本発明によると、複数の心周期からの信号
は信号の集合を形成するために時間同期され、時間の相
関ではなく、信号の集合の間の相関が求められる。この
ようにして、結果として相互の共分散マトリクスの大き
さは小さくなり、あらゆる時点での最適フィルタの決定
が可能になる。
【0008】本発明の装置の実施例によると、信号の集
合内の信号の間の相関は、ウィンドウユニットを用いて
順次オーバーラップさせたインターバルで求められる。
いかなる時点でも最適フィルタは、短期間の解析を介し
て決定される。このように、時間ひずみは避けられる
(消失する)。
【0009】
【実施例】次に、本発明による装置の実施例を添付図面
を参照して詳細に説明する。
【0010】図1は、本発明による装置の概略のブロッ
ク図を示す。図1において、記録されたECG信号は、
A/Dコンバータ2でデジタルの形に変換され、QRS
群は検出器4で検出される。平均化された信号が平均化
装置6で形成され、各心周期の少なくとも一部分の期間
中のECG信号の信号対雑音比を高めるために、後続の
ユニット10で用いられるインターバルまたは時間窓
は、ウインドウユニット8で決定される。
【0011】ユニット6、8、および10を図2を参照
して詳細に説明する。
【0012】複数の心周期からの記録されたECG信号
は、時間同期され、信号の集合を形成する。低域遮断周
波数を有するハイパスフィルタでの低周波成分の除去
は、いかなる時点でも最適フィルタの決定において重要
なことである。他の全ての信号処理と同様にハイパスフ
ィルタによる処理は、信号処理ユニット5で行われる。
【0013】平均化装置6において、信号の集合の平均
値が形成される。平均化は集合内の全ての信号の”純
然”平均化として実行する必要はなく、平均値を”改善
し”、代表する信号周期を提供するために、異なる別の
方法で実行することができる。したがって、たとえば平
均化において、大きな歪みのある集合内の信号を省くこ
とができ、平均化は、たとえばウインドウユニット8に
より決定されるインターバルと順次オーバーラップして
実行することができる。また、平均化を中央値にて決定
することもできる。
【0014】相関決定ユニット12において、信号の集
合内の各信号の間の相関は、最尤推定法により決定され
る。それにより得られる集合の相関は、次の式により与
えられる。
【0015】
【数1】
【0016】上式において、Zi はi番目のろ波された
ECG信号を示し、Nは考察される心周期の数を示し、
Mはインターバルまたは時間窓のサンプル数を示す。
【0017】このように求められた集合の相関は、でき
る限り大きな信号対雑音比を得るために、平均化装置6
から得られる平均信号値をろ波する最適フィルタ14を
決定するのに用いられる。最適フィルタは、連続する各
インターバルまたは時間窓に対して決定され、結果とし
てフィルタは時間依存型になる。したがって、フィルタ
は各観察インターバルまたは各サンプルに対して所定の
重みを有する。ユニット14でのろ波は、相関決定ユニ
ット12から得られる重み関数hn による乗算になる。
これらの関連について図3に、信号の集合内のN個の信
号と、フィルタ14の特性を決定する重み関数hn とが
図示されている。
【0018】このようにして、実現される瞬時最適フィ
ルタは、周波数選択特性をもたないフィルタであり、信
号対雑音比の大きさはフィルタの重みを決定するのに用
いられる。信号対雑音比は、利用可能な心周期の関数で
あり、瞬時に最適化されるフィルタも利用可能な心周期
の関数であるので、その出力信号はこのパラメータの関
数になる。明らかにこのフィルタは、心周期の数が少な
い場合、相当変動する。しかし、この変動は、使用され
る心周期の数が増加すると低下する。
【0019】
【発明の効果】心周期の数Nが増加と共に、平均化され
た信号内のノイズは減少し、上述のとおり、フィルタが
瞬時に最適化されると同時に信号対雑音比の向上に寄与
する状況は、より多くの情報を得ることを介して改善さ
れる。したがって、これらの効果はともに、心周期の数
が増加すると、フィルタからの出力信号の変動を低減す
る。
【図面の簡単な説明】
【図1】本発明による装置の概略のブロック図である。
【図2】図1の装置の一部をより詳細に示すブロック図
である。
【図3】本発明の説明に供する線図である。
【符号の説明】
2 A/Dコンバータ 4 検出器 5 信号処理ユニット 6 平均化装置 8 ウィンドウユニット 12 相関決定ユニット 14 最適フィルタ

Claims (14)

    【特許請求の範囲】
  1. 【請求項1】 複数の心周期からのA/D変換した信号
    を、信号の集合を形成するために時間同期し、集合内の
    信号の間の一致を、心周期の少なくとも一部分の期間中
    に少なくとも2つのサンプルを有する時間インターバル
    において求めることを特徴とするECG信号の信号対雑
    音比を高める方法。
  2. 【請求項2】 一致を信号のろ波を制御する重みで表わ
    すことを特徴とする請求項1記載の方法。
  3. 【請求項3】 平均値を前記信号の集合内の信号を用い
    て形成することを特徴とする請求項1または2記載の方
    法。
  4. 【請求項4】 一致を求めるために、信号の集合内の信
    号の間の相関を求めることを特徴とする請求項1から3
    までのいずれか1項記載の方法。
  5. 【請求項5】 前記集合内の信号の間の相関を、順次オ
    ーバーラップさせたインターバルで求めることを特徴と
    する請求項4記載の方法。
  6. 【請求項6】 信号の平均値を、相関により求められる
    重み関数に応じてろ波して取り出すことを特徴とする請
    求項3または4記載の方法。
  7. 【請求項7】 重み関数により、信号の平均値をろ波し
    て取り出すための時間で変化する特性を有するフィルタ
    を決定することを特徴とする請求項6記載の方法。
  8. 【請求項8】 ECG信号を記録し、A/D変換し、処
    理をする手段(2、4、5)を有するECG信号の信号
    対雑音比を高める装置において、同期ユニットは、信号
    の集合を形成するために、複数の心周期からの信号を時
    間同期するように構成され、ウィンドウユニット(8)
    により制御される計算ユニット(10、12)は、心周
    期の少なくとも一部分の期間中に少なくとも2つのサン
    プルを有する時間インターバルにおいて集合内の信号の
    間の一致を求めるように構成されていることを特徴とす
    るECG信号の信号対雑音比を高める装置。
  9. 【請求項9】 計算ユニット(12)は、一致により求
    められる重み関数に応じて信号をろ波するために、フィ
    ルタユニット(14)に接続されていることを特徴とす
    る請求項8記載の装置。
  10. 【請求項10】 同期ユニットまたはウィンドウユニッ
    トにより制御される平均化装置(6)は、前記集合内の
    信号の平均値を形成するように構成されていることを特
    徴とする請求項8または9記載の装置。
  11. 【請求項11】 計算ユニット(12)は、一致を求め
    るために、信号の集合内の信号の間の相関を求めるよう
    に構成されていることを特徴とする請求項8から10ま
    でのいずれか1項記載の装置。
  12. 【請求項12】 前記ウィンドウユニット(8)は、前
    記集合内の信号の順次オーバーラップさせたインターバ
    ルにおいて、相関を求めるために、計算ユニット(1
    0、12)を制御するように構成されていることを特徴
    とする請求項11記載の装置。
  13. 【請求項13】 計算ユニット(12)は、信号平均値
    をろ波するフィルタユニットの時間で変化する特性を決
    定する重み関数を生成することを特徴とする請求項10
    から12までのいずれか1項記載の装置。
  14. 【請求項14】 計算ユニット(12)は最尤推定器を
    有することを特徴とする請求項11から13までのいず
    れか1項記載の装置。
JP6164269A 1993-07-16 1994-07-15 Ecg信号の信号対雑音比を高める方法および装置 Pending JPH0751237A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
SE9302436-2 1993-07-16
SE19939302436A SE9302436D0 (sv) 1993-07-16 1993-07-16 Foerfarande och anordning foer att hoeja signal-brus-foerhaallande hos ekg-signaler

Publications (1)

Publication Number Publication Date
JPH0751237A true JPH0751237A (ja) 1995-02-28

Family

ID=20390626

Family Applications (1)

Application Number Title Priority Date Filing Date
JP6164269A Pending JPH0751237A (ja) 1993-07-16 1994-07-15 Ecg信号の信号対雑音比を高める方法および装置

Country Status (6)

Country Link
US (1) US5564428A (ja)
EP (1) EP0634134B1 (ja)
JP (1) JPH0751237A (ja)
AT (1) ATE176064T1 (ja)
DE (1) DE69416046T2 (ja)
SE (1) SE9302436D0 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101426591B1 (ko) * 2012-03-13 2014-08-06 연세대학교 산학협력단 생체 신호의 노이즈 제거 장치 및 방법

Families Citing this family (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5697378A (en) * 1996-12-16 1997-12-16 Siemens Medical Systems, Inc. Method and apparatus for using multiple leads for QRS detection
US5827195A (en) * 1997-05-09 1998-10-27 Cambridge Heart, Inc. Electrocardiogram noise reduction using multi-dimensional filtering
DE19805530C1 (de) 1998-02-11 1999-07-22 Marquette Hellige Gmbh Vorrichtung zur Signalverbesserung bei artefaktbehaftetem EKG
US6004270A (en) * 1998-06-24 1999-12-21 Ecton, Inc. Ultrasound system for contrast agent imaging and quantification in echocardiography using template image for image alignment
US6056691A (en) 1998-06-24 2000-05-02 Ecton, Inc. System for collecting ultrasound imaging data at an adjustable collection image frame rate
US7171265B2 (en) * 1998-07-31 2007-01-30 Harbinger Medical, Inc. Apparatus and method for detecting lead adequacy and quality
US6468213B1 (en) 1999-09-14 2002-10-22 Ecton, Inc. Medical diagnostic ultrasound system and method
US7678048B1 (en) 1999-09-14 2010-03-16 Siemens Medical Solutions Usa, Inc. Medical diagnostic ultrasound system and method
US6497664B1 (en) 1999-09-14 2002-12-24 Ecton, Inc. Medical diagnostic ultrasound system and method
US6436039B1 (en) 1999-09-14 2002-08-20 Ecton, Inc. Medicial diagnostic ultrasound system and method
US6488625B1 (en) 1999-09-14 2002-12-03 Ecton, Inc. Medical diagnostic ultrasound system and method
US6561979B1 (en) 1999-09-14 2003-05-13 Acuson Corporation Medical diagnostic ultrasound system and method
USD445189S1 (en) 1999-09-14 2001-07-17 Ecton, Inc. Medical diagnostic ultrasound system
US6312381B1 (en) 1999-09-14 2001-11-06 Acuson Corporation Medical diagnostic ultrasound system and method
US6524244B1 (en) 1999-09-14 2003-02-25 Ecton Inc. Medical diagnostic ultrasound system and method
US6508763B1 (en) 1999-09-14 2003-01-21 Ecton, Inc. Medical diagnostic ultrasound system and method
US6339715B1 (en) 1999-09-30 2002-01-15 Ob Scientific Method and apparatus for processing a physiological signal
IL133780A (en) * 1999-12-29 2009-05-04 B S P Biolog Signal Proc Ltd Method and device for analyzing a periodic or semi-periodic signal
WO2002017782A2 (de) * 2000-08-26 2002-03-07 Squid International Ag Verfahren und vorrichtung zum adaptiven reduzieren des rauschens in einem signal, insbesondere in einem elektro- oder magnetokardiographischen signal
US20070078353A1 (en) * 2005-10-04 2007-04-05 Welch Allyn, Inc. Method and apparatus for removing baseline wander from an ECG signal
US7894885B2 (en) * 2007-05-02 2011-02-22 Biosense Webster, Inc. Coherent signal rejection in ECG
WO2014102653A1 (en) 2012-12-31 2014-07-03 Koninklijke Philips N.V. Method and apparatus for reducing motion artifacts in ecg signals
US10602944B2 (en) 2014-05-08 2020-03-31 Drägerwerk AG & Co. KGaA Detecting artifacts in a signal
CN113935378A (zh) * 2021-10-14 2022-01-14 山东省人工智能研究院 一种基于对抗性深度全卷积网络的ecg降噪方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL61465A (en) * 1980-11-12 1984-09-30 Univ Ramot Method and apparatus for monitoring electrocardiogram(ecg)signals
US4802491A (en) * 1986-07-30 1989-02-07 Massachusetts Institute Of Technology Method and apparatus for assessing myocardial electrical stability
US4769760A (en) * 1987-03-13 1988-09-06 Cherne Medical, Inc. Terrain biased dynamic multiple threshold synchronization method and apparatus
US5033473A (en) * 1989-10-24 1991-07-23 Hewlett-Packard Company Method for discriminating pace pulse tails
US5092341A (en) * 1990-06-18 1992-03-03 Del Mar Avionics Surface ecg frequency analysis system and method based upon spectral turbulence estimation
US5226424A (en) * 1991-05-17 1993-07-13 Caliber Medical Corporation Low energy consumptive device for acquisition of data relating to abnormal heart muscle activity
US5255186A (en) * 1991-08-06 1993-10-19 Telectronics Pacing Systems, Inc. Signal averaging of cardiac electrical signals using temporal data compression and scanning correlation

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101426591B1 (ko) * 2012-03-13 2014-08-06 연세대학교 산학협력단 생체 신호의 노이즈 제거 장치 및 방법

Also Published As

Publication number Publication date
US5564428A (en) 1996-10-15
EP0634134B1 (en) 1999-01-20
ATE176064T1 (de) 1999-02-15
DE69416046T2 (de) 1999-09-30
EP0634134A1 (en) 1995-01-18
SE9302436D0 (sv) 1993-07-16
DE69416046D1 (de) 1999-03-04

Similar Documents

Publication Publication Date Title
JPH0751237A (ja) Ecg信号の信号対雑音比を高める方法および装置
DE19703403C2 (de) Vorrichtung und Verfahren zum Reduzieren von elektromagnetischen Störungen
US5117833A (en) Bi-spectral filtering of electrocardiogram signals to determine selected QRS potentials
Iwata et al. Algorithm for detecting the first and the second heart sounds by spectral tracking
US6668189B2 (en) Method and system for measuring T-wave alternans by alignment of alternating median beats to a cubic spline
US8172764B2 (en) Method and system for high-resolution extraction of quasi-periodic signals
JPH0852117A (ja) 患者モニタ
US6216031B1 (en) Apparatus for enhancing signals in ECGs artefacts
US5048535A (en) Method and apparatus for detecting QRS complex of ECG
US5479933A (en) Method and apparatus for processing ECG signals
Christov Dynamic powerline interference subtraction from biosignals
CN111839501B (zh) 一种基于频域自相关的心电信号干扰去噪方法及系统
Abboud et al. An improved detection algorithm in fetal electrocardiography
Jane et al. Adaptive filtering of high-resolution ECG signals
EP0848335B1 (en) Method and apparatus for QRS detection
Tun et al. Analysis of computer aided identification system for ECG characteristic points
Alhussainy QRS complex detection and R–R interval computation based on discrete wavelet transform
CN104902809B (zh) 用于减少ecg信号中的运动伪影的方法和装置
US9826937B2 (en) Method and apparatus for reducing motion artifacts in ECG signals
JP3864390B2 (ja) 脳波解析方法
Cano et al. Enhancement of low-level ECG components in noise with time-sequenced adaptive filtering
Atarius et al. Maximum likelihood analysis of cardiac late potentials
CN113288169A (zh) 一种心电波形信号的波形识别方法、装置及设备
Atarius et al. Detection of cardiac late potentials in nonstationary noise
Cloutier et al. Spectral analysis of closing sounds produced by lonescu-Shiley bioprosthetic aortic heart valves: Part 2 Computer simulation of aortic closing sounds and estimation of their truncation level and signal-to-noise ratio