Hi all,
I've just started to do some programming on wavelet using MATLAB. I've noticed that there are alot of ways to write the program. However, I always hit errors when I try to write the programs. In my program, I need to read off an ECG file (.DAT or .DCD) and then perform the wavelet transform and get the baseline wander and noise free signal. I'd implement the similar process using IIR BPF. However, I really hit a dead end with the wavelet method. Any advise or a starter to guide me along?Thanks!
These are my codes:
-------------------------------------------------------
clear all
load Qon
load Qoff
[Lo_D,Hi_D,Lo_R,Hi_R] = wfilters(wname);
subplot(3,1,1); stem(Lo_D);
title('Decomposition low-pass filter');
subplot(3,1,2); stem(Hi_D);
title('Decomposition high-pass filter');
subplot(3,1,3); stem(Lo_R);
title('Reconstruction low-pass filter');
%subplot(3,1,4); stem(Hi_R);
title('Reconstruction high-pass filter');
xlabel('The four filters for db5')
for file_no = 111
close all;
filename = sprintf('d_%05d',file_no);
load(filename)
figure(1);
plot(ecg),grid on;hold on;
end
-------------------------------------------------------
I've just started to do some programming on wavelet using MATLAB. I've noticed that there are alot of ways to write the program. However, I always hit errors when I try to write the programs. In my program, I need to read off an ECG file (.DAT or .DCD) and then perform the wavelet transform and get the baseline wander and noise free signal. I'd implement the similar process using IIR BPF. However, I really hit a dead end with the wavelet method. Any advise or a starter to guide me along?Thanks!
These are my codes:
-------------------------------------------------------
clear all
load Qon
load Qoff
[Lo_D,Hi_D,Lo_R,Hi_R] = wfilters(wname);
subplot(3,1,1); stem(Lo_D);
title('Decomposition low-pass filter');
subplot(3,1,2); stem(Hi_D);
title('Decomposition high-pass filter');
subplot(3,1,3); stem(Lo_R);
title('Reconstruction low-pass filter');
%subplot(3,1,4); stem(Hi_R);
title('Reconstruction high-pass filter');
xlabel('The four filters for db5')
for file_no = 111
close all;
filename = sprintf('d_%05d',file_no);
load(filename)
figure(1);
plot(ecg),grid on;hold on;
end
-------------------------------------------------------