sym3_anly.m

来自「这是一个课程设计」· M 代码 · 共 26 行

ecg=fopen('ECG1.dat','r');
N=1201;
data=fread(ecg,N,'int16');
data=data/10;
save ECGdata data;
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='sym3';
level=5;
[c,l]=wavedec(x,level,wname);
a5=wrcoef('a',c,l,'sym3',5);
a4=wrcoef('a',c,l,'sym3',4);
a3=wrcoef('a',c,l,'sym3',3);
a2=wrcoef('a',c,l,'sym3',2);
a1=wrcoef('a',c,l,'sym3',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]);