小波变换的去基线程序

function myfilter(); clear all; % readfile = 'F:\project\experimentdata\basedrift\data0427lead9.dat'; % readfile = 'F:\project\experimentdata\basedrift\data7.dat'; readfile = 'F:\project\experimentdata\basedrift\data19.dat'; fid = fopen(readfile,'r'); %把filename文件名的内容读出，放f中。 [f,count] = fread(fid,'short'); fclose(fid); %关闭该文件; datalength = 3000; s = f(1000:4000); maxlev=9; %分解6尺度 [C,L] = wavedec(s,maxlev,'coif3'); figure(1); %绘出分解后重构的的逼近信息和细节细节 AA1=wrcoef('a',C,L,'coif3',1); AA2=wrcoef('a',C,L,'coif3',2); AA3=wrcoef('a',C,L,'coif3',3); AA4=wrcoef('a',C,L,'coif3',4); AA5=wrcoef('a',C,L,'coif3',5); AA6=wrcoef('a',C,L,'coif3',6); AA7=wrcoef('a',C,L,'coif3',7); AA8=wrcoef('a',C,L,'coif3',8); AA9=wrcoef('a',C,L,'coif3',9); DD1=wrcoef('d',C,L,'coif3',1); DD2=wrcoef('d',C,L,'coif3',2); DD3=wrcoef('d',C,L,'coif3',3); DD4=wrcoef('d',C,L,'coif3',4); DD5=wrcoef('d',C,L,'coif3',5); DD6=wrcoef('d',C,L,'coif3',6); DD7=wrcoef('d',C,L,'coif3',7); DD8=wrcoef('d',C,L,'coif3',8); DD9=wrcoef('d',C,L,'coif3',9); %提取分解后逼近系数和细节系数 A1=appcoef(C,L,'coif3',1); A2=appcoef(C,L,'coif3',2); A3=appcoef(C,L,'coif3',3); A4=appcoef(C,L,'coif3',4); A5=appcoef(C,L,'coif3',5); A6=appcoef(C,L,'coif3',6); A7=appcoef(C,L,'coif3',7); A8=appcoef(C,L,'coif3',8); A9=appcoef(C,L,'coif3',9); D1=detcoef(C,L,1); D2=detcoef(C,L,2); D3=detcoef(C,L,3); D4=detcoef(C,L,4); D5=detcoef(C,L,5); D6=detcoef(C,L,6); D7=detcoef(C,L,7); D8=detcoef(C,L,8); D9=detcoef(C,L,9); yUp = 2000; yDn = -2000; % yUp = 500; % yDn = -500; subplot(10,2,1); plot(s); title('sampled ECG signal'); axis([0,datalength,-500,500]); subplot(10,2,3); plot(AA1); ylabel('A1'); axis([0,datalength,-500,500]); subplot(10,2,5); plot(AA2); ylabel('A2'); axis([0,datalength,-500,500]); subplot(10,2,7); plot(AA3); ylabel('A3'); axis([0,datalength,-500,500]); subplot(10,2,9); plot(AA4); ylabel('A4'); axis([0,datalength,-500,500]); subplot(10,2,11); plot(AA5); ylabel('A5'); axis([0,datalength,-500,500]); subplot(10,2,13); plot(AA6); ylabel('A6'); axis([0,datalength,-500,500]); subplot(10,2,15); plot(AA7); ylabel('A7'); axis([0,datalength,-500,500]); subplot(10,2,17); plot(AA8); ylabel('A8'); axis([0,datalength,-500,500]); subplot(10,2,19); plot(AA9); ylabel('A9'); axis([0,datalength,-500,500]); subplot(10,2,4); plot(DD1); ylabel('D1'); axis([0,datalength,-20,20]); subplot(10,2,6); plot(DD2); ylabel('D2'); axis([0,datalength,-50,50]); subplot(10,2,8); plot(DD3); ylabel('D3'); axis([0,datalength,-200,200]); subplot(10,2,10); plot(DD4); ylabel('D4'); axis([0,datalength,-200,200]); subplot(10,2,12); plot(DD5); ylabel('D5'); axis([0,datalength,-200,200]); subplot(10,2,14); plot(DD6); ylabel('D6'); axis([0,datalength,-200,200]); subplot(10,2,16); plot(DD7); ylabel('D7'); axis([0,datalength,-200,200]); subplot(10,2,18); plot(DD8); ylabel('D8'); axis([0,datalength,-200,200]); subplot(10,2,20); plot(DD9); ylabel('D9'); axis([0,datalength,-200,200]); figure(2); %将第一尺度置零 D1= zeros(1,length(D1)); D2= zeros(1,length(D2)); A9=zeros(1,length(A9)); C2 = [A9,D9',D8',D7',D6',D5',D4',D3',D2,D1]; % s2 = waverec(C2,L,'coif3'); subplot(311); plot(s); axis([0,3000,-300,550]); subplot(312); plot(s2); axis([0,3000,-300,550]); %噪声信号 for i = 1:datalength; noisesig(i) = s(i) - s2(i); end subplot(313); plot(noisesig); axis([0,3000,-300,550]);