load d1;
s = d1';
% 用db4对信号s进行7层小波分解
w = 'db4';
[c,l] = wavedec(s,7,w);
for i = 1:7
D(i,:) = wrcoef('d',c,l,w,i); %分解的7层结果分别顺序存在D的1到7行,即y1到y7
end
tt = 1:1000; % tt是x
subplot(9,1,1); plot(tt,s(tt),'r'); %画出原始数据
a7=s-sum(D); % 是s=a7+y7+y6+y5+y4+y3+y2+y1
subplot(9,1,2); plot(tt,a7(tt),'b');
title('Electrical Signal and Details');
for i = 1:7, subplot(9,1,i+2); plot(tt,D(7-i+1,tt),'g'); end % 画出y7到y1
% 对D中的数据,可以进行你想要的加工