A new LMS algorithm for analysis of atrial fibrillation signals

BackgroundA 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).MethodEquations 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.ResultsBased 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μV2/sample for the new LMS algorithm versus 0.72 ± 0.35μV2/sample for Widrow-Hoff LMS. The MSE for matching persistent CFAE averaged 0.55 ± 0.95μV2/sample for the new LMS algorithm versus 0.62 ± 0.55μV2/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.ConclusionsA 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.

[1]  Angelo B. Biviano,et al.  Different characteristics of complex fractionated atrial electrograms in acute paroxysmal versus long-standing persistent atrial fibrillation. , 2010, Heart rhythm.

[2]  B. Widrow,et al.  Adaptive noise cancelling: Principles and applications , 1975 .

[3]  Onkar Dabeer,et al.  Analysis of mean-square error and transient speed of the LMS Adaptive algorithm , 2002, IEEE Trans. Inf. Theory.

[4]  John J. Shynk,et al.  A successive cancellation algorithm for fetal heart-rate estimation using an intrauterine ECG signal , 2002, IEEE Trans. Biomed. Eng..

[5]  Bernard Widrow,et al.  The "Rubber-Mask" Technique I. Pattern Measurement and Analysis , 1973 .

[6]  K. Nademanee,et al.  A new approach for catheter ablation of atrial fibrillation: mapping of the electrophysiologic substrate. , 2004, Journal of the American College of Cardiology.

[7]  M. Akay,et al.  Biosignal pattern recognition and interpretation systems. 3. Methods of classification , 1994, IEEE Engineering in Medicine and Biology Magazine.

[8]  Luca Faes,et al.  A method for quantifying atrial fibrillation organization based on wave-morphology similarity , 2002, IEEE Transactions on Biomedical Engineering.

[9]  Bernard Widrow,et al.  The "rubber-mask" technique-II. Pattern storage and recognition , 1973, Pattern Recognit..

[10]  A. A. Goldstein,et al.  SENSOP: A derivative-free solver for nonlinear least squares with sensitivity scaling , 1993, Annals of Biomedical Engineering.

[11]  G. De Soete,et al.  Clustering and Classification , 2019, Data-Driven Science and Engineering.

[12]  Robert Plonsey,et al.  Bioelectricity: A Quantitative Approach Duke University’s First MOOC , 2013 .

[13]  E. Micheli-Tzanakou,et al.  The use of the differential steepest descent algorithm for adaptive template matching , 1992, Proceedings of the 1992 International Biomedical Engineering Days.

[14]  L. Faes,et al.  Principal component analysis and cluster analysis for measuring the local organisation of human atrial fibrillation , 2001, Medical and Biological Engineering and Computing.

[15]  Edward J. Ciaccio,et al.  Tonometric Arterial Pulse Sensor With Noise Cancellation , 2008, IEEE Transactions on Biomedical Engineering.

[16]  Bernard Widrow,et al.  A comparison of adaptive algorithms based on the methods of steepest descent and random search , 1976 .

[17]  C. Caraiscos,et al.  A roundoff error analysis of the LMS adaptive algorithm , 1984 .

[18]  K. Egiazarian,et al.  A transform domain LMS adaptive filter with variable step-size , 2002, IEEE Signal Processing Letters.

[19]  E. Micheli-Tzanakou,et al.  Characterization of spontaneous changes in electrogram morphology , 1994, Computers in Cardiology 1994.

[20]  B. Widrow,et al.  A variable leaky LMS adaptive algorithm , 2004, Conference Record of the Thirty-Eighth Asilomar Conference on Signals, Systems and Computers, 2004..

[21]  M. Akay,et al.  Biosignal pattern recognition and interpretation systems , 1993, IEEE Engineering in Medicine and Biology Magazine.

[22]  A. Gray,et al.  Unconstrained frequency-domain adaptive filter , 1982 .

[23]  P. van Gerwen,et al.  Digital echo cancellation for baseband data transmission , 1979 .

[24]  Peter J. McLane,et al.  A fast adaptive algorithm for MMSE receivers in DS-CDMA systems , 2004, IEEE Signal Processing Letters.

[25]  Edward J Ciaccio,et al.  Differences in Repeating Patterns of Complex Fractionated Left Atrial Electrograms in Longstanding Persistent Atrial Fibrillation as Compared With Paroxysmal Atrial Fibrillation , 2011, Circulation. Arrhythmia and electrophysiology.

[26]  V Barbaro,et al.  Mapping the Organization of Atrial Fibrillation with Basket Catheters Part II: Regional Patterns in Chronic Patients , 2001, Pacing and clinical electrophysiology : PACE.

[27]  D Durrer,et al.  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 , 1978, Circulation research.

[28]  Yang-yu Fan,et al.  A Variable Step-Size LMS Adaptive Filtering Algorithm , 2009, 2009 5th International Conference on Wireless Communications, Networking and Mobile Computing.

[29]  Harald Haas,et al.  Asilomar Conference on Signals, Systems, and Computers , 2006 .

[30]  Robert Plonsey,et al.  Bioelectric Potentials and Currents , 2000 .

[31]  Edward J. Ciaccio,et al.  Development of Gradient Descent Adaptive Algorithms to Remove Common Mode Artifact for Improvement of Cardiovascular Signal Quality , 2007, Annals of Biomedical Engineering.

[32]  R. Merletti,et al.  Advances in processing of surface myoelectric signals: Part 2 , 2006, Medical and Biological Engineering and Computing.

[33]  Khaled A. Mayyas,et al.  Performance analysis of the deficient length LMS adaptive algorithm , 2005, IEEE Transactions on Signal Processing.

[34]  V. Barbaro,et al.  Mapping the Organization of Atrial Fibrillation with Basket Catheters Part I: Validation of a Real‐Time Algorithm , 2001, Pacing and clinical electrophysiology : PACE.

[35]  R. Merletti,et al.  Advances in processing of surface myoelectric signals: Part 1 , 1995, Medical and Biological Engineering and Computing.

[36]  Bernard Widrow,et al.  The "Rubber-Mask" Technique II. Pattern Storage and Recognition , 1973 .

[37]  C. Henriquez,et al.  Study of Unipolar Electrogram Morphology in a Computer Model of Atrial Fibrillation , 2003, Journal of cardiovascular electrophysiology.

[38]  E. Micheli-Tzanakou,et al.  A least squares algorithm with one phase weight for adaptive cancelation of noise from biomedical signals , 1992, [1992] Proceedings of the Eighteenth IEEE Annual Northeast Bioengineering Conference.