hwden.m

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

%基于小波变换的心电信号的压缩
ecg=fopen('ECG1.dat','r');
N=1201;
data=fread(ecg,N,'int16');
data=data/10000;
fclose(ecg);
x=data;
wavename='db5';
[c,l]=wavedec(x,4,wavename);
alpha=2;
sorh='h';
[thr,nkeep]=wdcbm(c,l,alpha);%使用硬阈值压缩信号
[xc,cxc,lxc,perf0,perfl2]=wdencmp('lvd',c,l,wavename,4,thr,sorh);
t=0:0.004:4.8;
subplot(511);
plot(t,x);
title('原始的ECG信号');
axis([0 5 -3 3]);
subplot(512);
plot(t,xc);
axis([0 5 -3 3]);
title('在第4层对高频系数量化压缩后的ECG信号（alpha=2 db5）');

[c,l]=wavedec(x,4,wavename);
alpha=4;
sorh='h';
[thr,nkeep]=wdcbm(c,l,alpha);%使用硬阈值压缩信号
[xc1,cxc1,lxc1,perf01,perfl21]=wdencmp('lvd',c,l,wavename,4,thr,sorh);
subplot(513);
plot(t,xc1);
axis([0 5 -3 3]);
title('在第4层对高频系数量化压缩后的ECG信号alpha=4 db5');


[c,l]=wavedec(x,2,wavename);
alpha=4;
sorh='h';
[thr,nkeep]=wdcbm(c,l,alpha);%使用硬阈值压缩信号
[xc1,cxc1,lxc1,perf01,perfl21]=wdencmp('lvd',c,l,wavename,2,thr,sorh);
subplot(514);
plot(t,xc1);
axis([0 5 -3 3]);
title('在第2层对高频系数量化压缩后的ECG信号alpha=4 db5');


[c,l]=wavedec(x,2,'db3');
alpha=4;
sorh='h';
[thr,nkeep]=wdcbm(c,l,alpha);%使用硬阈值压缩信号
[xc1,cxc1,lxc1,perf01,perfl21]=wdencmp('lvd',c,l,'db3',2,thr,sorh);
subplot(515);
plot(t,xc1);
axis([0 5 -3 3]);
title('在第2层对高频系数量化压缩后的ECG信号alpha=4 db3');