wavelet_decompose.m

小波变换的条件限制的实现的编程

%信号的小波多分辨率分解与重构
%
%example of function=(s=cos(t*pi/50))
%
%信号的分解
t=0:pi/50:6*pi;
s=cos(t);
plot(s);
title('原始信号');
[c,l]=wavedec(s,3,'db1');%对S进行小波分解；用dbl小波进行3层分解
a3=appcoef(c,l,'db1',3);%提取小波低频系数；

d3=detcoef(c,l,3);
d2=detcoef(c,l,2);
d1=detcoef(c,l,1);

figure(2)
subplot(5,2,1);
plot(a3);
ylabel('a3');
subplot(5,2,3);
plot(d3);
ylabel('d3');
subplot(5,2,5);
plot(d2);
ylabel('d2');
subplot(5,2,7);
plot(d1);
ylabel('d1');


% 信号的重构
s1=waverec(c,l,'db1');
subplot(5,2,9);
plot(s1);
Ylabel('s1');

[c,l]=wavedec(s,3,'coif3');
a3=appcoef(c,l,'coif3',3);
d3=detcoef(c,l,3);
d2=detcoef(c,l,2);
d1=detcoef(c,l,1);
subplot(5,2,2);
plot(a3);
ylabel('a3');
subplot(5,2,4);
plot(d3);
ylabel('d3');
subplot(5,2,6);
plot(d2);
ylabel('d2');
subplot(5,2,8);
plot(d1);
ylabel('d1');

s2=waverec(c,l,'coif3');
subplot(5,2,10);
plot(s2);
ylabel('s2');