matlab例程

开发平台：

Matlab

db6chonggou.m：源码内容

							
t=0:0.01:20.5;  
s=5*sin(2*pi*5*t);
figure(1);
subplot(6,1,1); plot(s);
title('db6小波分解upcoef重构');
ylabel('s');
axis([0,2500, -5, 5]);
[c,l] = wavedec(s, 4, 'db6');
a1 = appcoef(c, l, 'db6', 1);
cd1 = detcoef(c,l,1);
scd1 = upcoef('d', cd1, 'db6', 1);
subplot(6,1,2); plot(scd1);
ylabel('d1');
axis([0,2500, -0.5, 0.5]);
a2 = appcoef(c, l, 'db6', 2);
cd2 = detcoef(c, l, 2);
scd2 = upcoef('d', cd2, 'db6', 2);
subplot(6,1,3); plot(scd2);
ylabel('d2');
axis([0,2500, -0.5, 0.5]);
a3 = appcoef(c, l, 'db6', 3);
cd3 = detcoef(c, l, 3);
scd3 = upcoef('d', cd3, 'db6', 3);
subplot(6,1,4); plot(scd3);
ylabel('d3');
axis([0,2500, -5, 5]);
cd4 = detcoef(c, l, 4);
scd4 = upcoef('d', cd4, 'db6', 4);
subplot(6,1,5); plot(scd4);
ylabel('d4');
axis([0,2500, -10, 10]);
ca4 = appcoef(c, l, 'db6', 4);
sca4 = upcoef('a', ca4, 'db6', 4);
subplot(6,1,6); plot(sca4);
ylabel('a4');
axis([0,2500, -5, 5]);
%------------------------------------------------------------------------
figure(2);
subplot(6,1,1); plot(s);
title('db6小波分解wrcoef重构');
ylabel('s');
axis([0,2500, -5, 5]);
[c1,l1]=wavedec(s,1,'db6');
scd1=wrcoef('d',c1,l1,'db6',1);
subplot(6,1,2); plot(scd1);
ylabel('d1');
axis([0,2500, -0.5, 0.5]);
[c2,l2]=wavedec(s,2,'db6');
scd2=wrcoef('d',c2,l2,'db6',2);
subplot(6,1,3); plot(scd2);
ylabel('d2');
axis([0,2500, -0.5, 0.5]);
[c3,l3]=wavedec(s,3,'db6');
scd3=wrcoef('d',c3,l3,'db6',3);
subplot(6,1,4); plot(scd3);
ylabel('d3');
axis([0,2500, -5, 5]);
[c4,l4]=wavedec(s,4,'db6');
scd4=wrcoef('d',c4,l4,'db6',4);
subplot(6,1,5); plot(scd4);
ylabel('d4');
axis([0,2500, -10, 10]);
sca4=wrcoef('a',c4,l4,'db6',4);
subplot(6,1,6); plot(sca4);
ylabel('a4');
axis([0,2500, -5, 5]);

						