A new LMS algorithm for analysis of atrial fibrillation signals

BioMedical Engineering OnLine

Hasan Garan¹  William Whang¹  Angelo B Biviano¹  Edward J Ciaccio² 

[1]Department of Medicine – Division of Cardiology, Columbia University Medical Center, New York, USA

[2]Columbia University, Harkness Pavilion 804, 180 Fort Washington Avenue, New York, NY, 10032, USA

关键词: Mean-squared error;    LMS algorithm;    Fractionation;    F wave;    Electrocardiogram;    Atrial fibrillation;   

Others  :  798111
DOI  :  10.1186/1475-925X-11-15

 received in 2012-01-17, accepted in 2012-03-14,  发布年份 2012

【 摘 要 】

Background

A biomedical signal can be defined by its extrinsic features (x-axis and y-axis shift and scale) and intrinsic features (shape after normalization of extrinsic features). In this study, an LMS algorithm utilizing the method of differential steepest descent is developed, and is tested by normalization of extrinsic features in complex fractionated atrial electrograms (CFAE).

Method

Equations for normalization of x-axis and y-axis shift and scale are first derived. The algorithm is implemented for real-time analysis of CFAE acquired during atrial fibrillation (AF). Data was acquired at a 977 Hz sampling rate from 10 paroxysmal and 10 persistent AF patients undergoing clinical electrophysiologic study and catheter ablation therapy. Over 24 trials, normalization characteristics using the new algorithm with four weights were compared to the Widrow-Hoff LMS algorithm with four tapped delays. The time for convergence, and the mean squared error (MSE) after convergence, were compared. The new LMS algorithm was also applied to lead aVF of the electrocardiogram in one patient with longstanding persistent AF, to enhance the F wave and to monitor extrinsic changes in signal shape. The average waveform over a 25 s interval was used as a prototypical reference signal for matching with the aVF lead.

Results

Based on the derivation equations, the y-shift and y-scale adjustments of the new LMS algorithm were shown to be equivalent to the scalar form of the Widrow-Hoff LMS algorithm. For x-shift and x-scale adjustments, rather than implementing a long tapped delay as in Widrow-Hoff LMS, the new method uses only two weights. After convergence, the MSE for matching paroxysmal CFAE averaged 0.46 ± 0.49μV²/sample for the new LMS algorithm versus 0.72 ± 0.35μV²/sample for Widrow-Hoff LMS. The MSE for matching persistent CFAE averaged 0.55 ± 0.95μV²/sample for the new LMS algorithm versus 0.62 ± 0.55μV²/sample for Widrow-Hoff LMS. There were no significant differences in estimation error for paroxysmal versus persistent data. From all trials, the mean convergence time was approximately 1 second for both algorithms. The new LMS algorithm was useful to enhance the electrocardiogram F wave by subtraction of an adaptively weighted prototypical reference signal from the aVF lead. The extrinsic weighting over 25 s demonstrated that time-varying functions such as patient respiration could be identified and monitored.

Conclusions

A new LMS algorithm was derived and used for normalization of the extrinsic features in CFAE and for electrocardiogram monitoring. The weighting at convergence provides an estimate of the degree of similarity between two signals in terms of x-axis and y-axis shift and scale. The algorithm is computationally efficient with low estimation error. Based on the results, proposed applications include monitoring of extrinsic and intrinsic features of repetitive patterns in CFAE, enhancement of the electrocardiogram F wave and monitoring of time-varying signal properties, and to quantitatively characterize mechanistic differences in paroxysmal versus persistent AF.

【 授权许可】

   
2012 Ciaccio et al.; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

Figure 4.

Figure 5.

Figure 6.

Figure 7.

Figure 8.

【 参考文献 】

• [1]Widrow B: The Rubber Mask Technique (Part I). Pattern Measurement and Analysis. Pat Rec 1973, 5:175-197.
• [2]Widrow B: The Rubber Mask Technique (Part II). Pattern Storage and Recognition. Pat Rec 1973, 5:199-211.
• [3]Ciaccio EJ, Dunn SM, Akay M: Biosignal Pattern-Recognition and Interpretation Systems. 1. Fundamental Concepts. IEEE Eng Med Biol Mag 1993, 12:89-98.
• [4]Ciaccio EJ, Dunn SM, Akay M: Biosignal Pattern-Recognition and Interpretation Systems. 3. Clustering and Classification. IEEE Eng Med Biol Mag 1994, 13:129-135.
• [5]Widrow B, Glover JR, McCool JM, Kavnitz J, Williams CS, Hearn RH, Zeidler JR, Dong E: Adaptive noise cancelling: Principles and applications. Proc IEEE 1975, 63:1692-1716.
• [6]Widrow B, McCool JM, Larimore MG, Johnson CR: Stationary and nonstationary learning characteristics of the least mean squares adaptive filter. Proc IEEE 1976, 64:1151-1162.
• [7]Lai KC, Shynk JJ: A successive cancellation algorithm for fetal heart-rate estimation using an intrauterine ECG signal. IEEE Trans Biomed Engin 2002, 49:943-954.
• [8]Widrow B, McCool JM: A comparison of adaptive algorithms based on the method of steepest descent and random search. IEEE Trans Ant Prop 1976, 24:615-636.
• [9]Nademanee K, McKenzie J, Kosar E, Schwab M, Sunsaneewitayakul B, Vasavakul T, Khunnawat C, Ngarmukos T: A new approach for catheter ablation of atrial fibrillation: mapping of the electrophysiologic substrate. J Am Coll Cardiol 2004, 43:2044-2053.
• [10]Kléber AG, Janse MJ, van Capelle FJ, Durrer D: Mechanism and time course of S-T and T-Q segment changes during acute regional myocardial ischemia in the pig heart determined by extracellular and intracellular recordings. Circ Res 1978, 42:603-613.
• [11]Plonsey R, Barr RC: Plonsey R, Barr RC. In Bioelectricity: A Quantitative Approach. 2nd edition. Springer, New York; 2000:39-65.
• [12]Jacquemet V, Virag N, Ihara Z, Dang L, Blanc O, Zozor S, Vesin JM, Kappenberger L, Henriquez C: Study of unipolar electrogram morphology in a computer model of atrial fibrillation. Journal Cardiovasc Electrophys 2003, 14(10 Suppl):S172-179.
• [13]Lo Conte LR, Merletti R: Advances in processing of surface myoelectric signals: Part 2. Med Biol Eng Comput 1995, 33(3 Spec No):373-384.
• [14]Mansour D, Gray A: Unconstrained frequency domain adaptive filter. IEEE Trans Acoust Speech Sig Proc 1982, 30:726-734.
• [15]Ciaccio EJ, Micheli-Tzanakou E, Dunn SM, Wit AL: The use of the differential steepest descent algorithm for adaptive template matching. Proceedings International Biomedical Engineering Days 1992, 1:198-202.
• [16]Ciaccio EJ, Micheli-Tzanakou E: A least squares algorithm with one phase weight for adaptive cancelation of noise from biomedical signals. Proceedings of the Eighteenth IEEE Annual Northeast Biomedical Engineering Conference 1992, 1:41-42.
• [17]Ciaccio EJ, Dunn SM, Akay M, Micheli-Tzanakou E, Coromilas J, Costeas CA, Wit AL: Characterization of spontaneous changes in electrogram morphology. IEEE Computers in Cardiology 1994, 1:701-704.
• [18]Ciaccio EJ, Biviano AB, Whang W, Gambhir A, Garan H: Different characteristics of complex fractionated atrial electrograms in acute paroxysmal versus long-standing persistent atrial fibrillation. Heart Rhythm 2010, 7:1207-1215.
• [19]Ciaccio EJ, Biviano AB, Whang W, Vest JA, Gambhir A, Einstein AJ, Garan H: Differences in repeating patterns of complex fractionated left atrial electrograms in longstanding persistent atrial fibrillation as compared with paroxysmal atrial fibrillation. Circ Arrhythm Electrophysiol 2011, 4:470-477.
• [20]Chan IS, Goldstein AA, Bassingthwaighte JB: SENSOP: A Derivative-Free Solver for Nonlinear Least Squares with Sensitivity Scaling. Ann Biomed Eng 1993, 21:621-631.
• [21]Ciaccio EJ, Micheli-Tzanakou E: Development of gradient descent adaptive algorithms to remove common mode artifact for improvement of cardiovascular signal quality. Ann Biomed Eng 2007, 35:1146-1155.
• [22]Ciaccio EJ, Drzewiecki GM: Tonometric arterial pulse sensor with noise cancellation. IEEE Trans Biomed Eng 2008, 55:2388-2396.
• [23]Caraiscos C, Liu B: A roundoff error analysis of the LMS adaptive algorithm. IEEE Transactions on Acoustics, Speech and Signal Processing 1984, 32:34-41.
• [24]Bilcu RC, Kuosmanen P, Egiazarian K: A transform domain LMS adaptive filter with variable step-size. IEEE Signal Processing Letters 2002, 9:51-53.
• [25]Li XB, Fan YY, Peng K: A variable step-size LMS adaptive filtering algorithm. WiCom '09: 5th International Conference on Wireless Communications, Networking and Mobile Computing. Beijing; 2009:1-4.
• [26]Kamenetsky M, Widrow B: A variable leaky LMS adaptive algorithm'. Record of the Thirty-Eighth Asilomar Conference on Signals, Systems and Computers 2004, 1:125-128.
• [27]Dabeer O, Masry E: Analysis of mean-square error and transient speed of the LMS adaptive algorithm. IEEE Transactions Information Theory 2002, 48:1873-1894.
• [28]Faes L, Nollo G, Kirchner M, Olivetti E, Gaita F, Riccardi R, Antolini R: Principal component analysis and cluster analysis for measuring the local organisation of human atrial fibrillation. Med Biol Eng Comput 2001, 39:656-663.
• [29]Faes L, Nollo G, Antolini R, Gaita F, Ravelli F: A method for quantifying atrial fibrillation organization based on wave-morphology similarity. IEEE Trans Biomed Eng 2002, 49(12 Pt 2):1504-1513.
• [30]Barbaro V, Bartolini P, Calcagnini G, Censi F, Morelli S, Michelucci A: Mapping the organization of atrial fibrillation with basket catheters. Part I: Validation of a real-time algorithm. Pacing Clin Electrophysiol 2001, 24:1082-1088.
• [31]Hamouda WA, McLane PJ: A fast adaptive algorithm for MMSE receivers in DS-CDMA systems. IEEE Signal Processing Letters 2004, 11:86-89.
• [32]Verhueckx NAM, van der Elzen HC, Snijders FAM, van Gemin PJ: Distal echo cancellation for base band data transmission. IEEE Trans Acous Speech Sig Proc 1979, ASSP-27:768-781.
• [33]Mayyas K: Performance analysis of the deficient length LMS adaptive algorithm. IEEE Transactions on Signal Processing 2005, 53:2727-2734.