wname='db4';
[Lo_D,Hi_D,Lo_R,Hi_R]=wfilters(wname);
% subplot(2,2,1);
% stem(Lo_D);
% title('分解低通滤波器');
% subplot(2,2,2); stem(Hi_D);
% title('分解高通滤波器');
% subplot(2,2,3); stem(Lo_R);
% title('重构低通滤波器');
% subplot(2,2,4); stem(Hi_R);
% title('重构高通滤波器');
file_t = fopen('波高数据3.txt','r');
A = fscanf(file_t,'%d');
subplot(2,2,1);
% figure(1);
plot(A(1:256));
x=A(1:256);
hold on;
% N=100;fs=1000; n=0:N-1; t=n/fs;
% x=sin(60*pi*t)+cos(100*pi*t); %原信号
%subplot(2,2,1); plot(x,'-r','LineWidth',1.5);
dim=4;
[cA1,cD1] = mydwt(x,Lo_D,Hi_D,dim);
[cA2,cD2]=mydwt(cA1,Lo_D,Hi_D,dim);
[cA3,cD3]=mydwt(cA2,Lo_D,Hi_D,dim);
[cA4,cD4]=mydwt(cA3,Lo_D,Hi_D,dim);
for i=1:length(cA)
if(abs(cA(i))<100) cA(i)=0;
end
cA(i)=floor (cA(i));
end
for i=1:length(cD)
if(abs(cD(i))<100) cD(i)=0;
end
cD(i)=floor (cD(i));
end
y1 = myidwt(cA1,cD1,Lo_R,Hi_R);
y2 = myidwt(cA2,cD2,Lo_R,Hi_R);
y3 = myidwt(cA3,cD3,Lo_R,Hi_R);
y4= myidwt(cA4,cD4,Lo_R,Hi_R);
subplot(2,2,2);
plot(y(40:200),'r','LineWidth',1.5);
subplot(2,2,3); stem(cA);
title('平均部分小波系数');
subplot(2,2,4); stem(cD);
title('细节部分小波系数');