One-Dimensional Signal Denoising Based on Kalman Filter (MATLAB)
clear all;load idealECG.mat;fs = 128;N = length(idealECG);t = (0:N-1)/fs; % SNR = 10;n_50 = 0.2 * sin(2 * pi * 50 * t);% nECG = awgn(idealECG,SNR,'measured')+n_50;% % save("nECG.mat","nECG");load nECG.mat; std(nECG-idealECG) std(nECG-idealECG-n_50) figure(1);subplot(2,1,1);plot(t,idealECG);xlim([5,10]);title('ideal ECG'); subplot(2,1,2);plot(t,nECG);xlim([5,10]);title("noisy ECG"); Scaler Kalman filter % Initialize the Kalman filter parametersA = 1; % State transition matrixH = 1; % Observation matrixR … Read more
