function myfilter(); 
clear all; 
% readfile = 'F:\project\experimentdata\basedrift\data0427lead9.dat';  
% readfile = 'F:\project\experimentdata\basedrift\data7.dat';  
readfile = 'F:\project\experimentdata\basedrift\data18.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('original'); 
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,15); 
plot(AA8); 
ylabel('A8'); 
axis([0,datalength,-500,500]); 
subplot(10,2,15); 
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,16); 
 plot(DD8); 
ylabel('D8'); 
axis([0,datalength,-200,200]); 
subplot(10,2,16); 
 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(A7)); 
 C2 = [A9,D9',D8',D7',D6',D5',D4',D3',D2',D1'];  % 
 s2 = waverec(C2,L,'coif3'); 
  
 subplot(311); 
 plot(s); 
 axis([0,datalength,yDn,yUp]); 
  
 %濾工頻 
b = [1,1,1,1,1,1,1,1,1,1]; 
a = [10]; 
nFactor = 1.5   %工頻需要適當放大濾波后的數據 
% b = [0.95,0,0,0,0,0,0,0,0,0,-0.95]; 
% a = [1,0,0,0,0,0,0,0,0,0,-0.9]; 
% nFactor = 1;   %工頻需要適當放大濾波后的數據 
 
yy = filter(b,a,s2); 
%修改幅值 
for i=1:datalength 
    yy(i) = yy(i)*nFactor; 
end 
 subplot(312); 
 plot(yy); 
 axis([0,datalength,yDn,yUp]); 
  
%噪聲信號 
for i = 1:datalength; 
    noisesig(i) = s(i) - yy(i); 
end 
subplot(313); 
 plot(noisesig); 
 axis([0,datalength,yDn,yUp]);