Adaptive augmented cubature Kalman filter/smoother for ECG denoising

Model-based Bayesian approaches have been widely applied in Electrocardiogram (ECG) signal processing, where their performances heavily rely on the accurate selection of model parameters, particularly the state and measurement noise covariance matrices. In this study, we introduce an adaptive augmented cubature Kalman filter/smoother (CKF/CKS) for ECG processing, which updates the noise covariance matrices at each time step to accommodate diverse noise types and input signal-to-noise ratios (SNRs). Additionally, we incorporate the dynamic time warping technique to enhance the filter's efficiency in the presence of heart rate variability. Furthermore, we propose a method to significantly reduce the computational complexity required for CKF/CKS implementation in ECG processing. The denoising performance of the proposed filter was compared to those of various nonlinear Kalman-based frameworks involving the Extended Kalman filter/smoother (EKF/EKS), the unscented Kalman filter/smoother (UKF/UKS), and the ensemble Kalman filter (EnKF) that was recently proposed for ECG enhancement. In this study, we conducted a comprehensive evaluation and comparison of the performance of various nonlinear Kalman-based frameworks for ECG signal processing, which have been proposed in recent years. Our assessment was carried out on multiple normal ECG segments extracted from different entries in the MIT-BIH Normal Sinus Rhythm Database (NSRDB). This database provides a diverse set of ECG recordings, allowing us to examine the filters' denoising capabilities across various scenarios. By comparing the performance of these filters on the same dataset, we aimed to provide a thorough analysis and identification of the most effective approach for ECG denoising. Two kinds of noises were introduced to such segments: 1-stationary white Gaussian noise and 2-non-stationary real muscle artifact noise. For evaluation, four comparable measures namely the SNR improvement, PRD, correlation coefficient and MSEWPRD were employed. The findings demonstrated that the suggested algorithm outperforms the EKF/EKS, EnKF/EnKS, UKF/UKS methods in both stationary and nonstationary environments regarding SNR improvement, PRD, correlation coefficient and MSEWPRD metrics.