ecg=fopen('ECG1.dat','r');
N=1201;
data=fread(ecg,N,'int16');
data=data/10;
fclose(ecg);
x=0:0.004:4.8;
figure(1);
subplot(321);
plot(x,data);
axis([0 4.8 -4000 4000]);
title('原始心電信號');
x=data;
wname='db6';
level=5;
[c,l]=wavedec(x,level,wname);
a5=wrcoef('a',c,l,'db6',5);
a4=wrcoef('a',c,l,'db6',4);
a3=wrcoef('a',c,l,'db6',3);
a2=wrcoef('a',c,l,'db6',2);
a1=wrcoef('a',c,l,'db6',1);
subplot(322);plot(a5);title('a5');axis([0 1201 -2000 2000]);
subplot(323);plot(a4);title('a4');axis([0 1201 -2000 2000]);
subplot(324);plot(a3);title('a3');axis([0 1201 -2000 2000]);
subplot(325);plot(a2);title('a2');axis([0 1201 -2000 2000]);
subplot(326);plot(a1);title('a1');axis([0 1201 -2000 2000]);