matlab心电信号小波分解 去噪

clear; %读取ECG信号 fid=fopen('ecg.txt','r'); x=fscanf(fid,'%f %f %f',[600 3]); fclose(fid); N=length(x); for i=0:N/3-1 t(i+1)=x(3*i+1); sig1(i+1)=x(3*i+2); end s=sig1+0.6*randn(1,length(sig1)); figure(1); subplot(2,1,1); plot(t,s,'LineWidth',2); xlabel('时间'); axis([0 1.2 -inf inf]); ylabel('幅值'); title('原始ECG信号'); grid on %选用基波函数 wavelet='db5'; %分解级数 level=3; %压缩参数（去噪时alpha=3） alpha=1.5; %阈值类型（硬阀值） sorh='h'; %小波分解 [c,l]=wavedec(s,level,wavelet); %birge-massart算法计算阈值 [thr,nkeep]=wdcbm(c,l,alpha); %使用硬阈值压缩信号 %‘lvd’为允许设置各层阀值，‘gbl阀值’为固定阀值 [xc,cxc,lxc,perf0,perfl2]=wdencmp('lvd',c,l,wavelet,level,thr,sorh); subplot(2,1,2); plot(t,xc,'LineWidth',2); xlabel('时间'); axis([0 1.2 -inf inf]); ylabel('幅值'); title('去噪ECG信号'); grid on