example6_2.asv

用MATLAB实现对一语音信号用不同小波进行分解

t=0:1:100*pi;
s=sin(3*t)+sin(0.3t)+sin(0.03t);
subplot(6,2,1);plot(s);
title('原始信号s');
%====================================
%对s进行小波分解：db3   5层
[c,l]=wavedec(s,5,'db3');
%====================================
%提取小波分解的低频系数
a5=appcoef(c,l,'db3',5);
a4=appcoef(c,l,'db3',4);
a3=appcoef(c,l,'db3',3);
a2=appcoef(c,l,'db3',2);
a1=appcoef(c,l,'db3',1);
%====================================
%提取小波分解的各层高频系数
d5=detcoef(c,l,5);
d4=detcoef(c,l,4);
d3=detcoef(c,l,3);
d2=detcoef(c,l,2);
d1=detcoef(c,l,1);	
%====================================
%绘出各系数的图形
subplot(6,2,3);plot(a5);
Ylabel('a5');
subplot(6,2,5);plot(a4);
Ylabel('a4');
subplot(6,2,7);plot(a3);
Ylabel('a3');
subplot(6,2,9);plot(a2);
Ylabel('a2');
subplot(6,2,11);plot(a1);
Ylabel('a1');
subplot(6,2,4);plot(d5);
Ylabel('d5');
subplot(6,2,6);plot(d4);
Ylabel('d4');
subplot(6,2,6);plot(d3);
Ylabel('d3');
subplot(6,2,8);plot(d2);
Ylabel('d2');
%====================================
%重构信号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');