rwden.m

这是一个课程设计

%基于小波变换的心电信号的降噪
ecg=fopen('ECG1.dat','r');
N=1201;
data=fread(ecg,N,'int16');
data=data/10000;
fclose(ecg);
x=data;
wavename='db5';
level=4;
[c,l]=wavedec(x,level,wavename);
alpha=1.5;
sorh='h';
[thr,nkeep]=wdcbm(c,l,alpha);%使用硬阈值给信号降噪
[xc,cxc,lxc,perf0,perfl2]=wdencmp('lvd',c,l,wavename,level,thr,sorh);
t1=0:0.004:(length(x)-1)*0.004;
subplot(211);
plot(t1,x);
title('从人体采集的原始的ECG信号');
subplot(212);
plot(t1,xc);
title('Birge-Massart策略阈值降噪后的ECG信号（wname=db5 level=4）');