example6_2.m

<小波分析理论与MATLAB 7实现>是“MATLAB应用技术”系列丛书之一

t=0:pi/100:4*pi;
s=sin(t+3*pi/4);
plot(s);
title('原始信号s');
%====================================
%对s进行小波分解：db1  3层
[c,l]=wavedec(s,3,'db1');
%====================================
%提取小波分解的低频系数a3
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');
%====================================
%重构信号s
s1=waverec(c,l,'db1');
subplot(5,2,9);plot(s1);
Ylabel('s1');
%====================================
%下面用小波‘coif3’重复上述过程
[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');