example6_1.m

来自「电子工业出版社《matlab6.5辅助小波分析与应用》的源代码」· M 代码 · 共 24 行

%首先装载一原始信号  
load nearbrk;  
s=nearbrk;  
%===========================
%画出该信号的波形图  
subplot(4,2,1);  
plot(s);  
Ylabel('s');  
title('原始信号s和信号的近似a、细节d');  
%===========================
%用小波db2进行3层分解  
[c,l]=wavedec(s,3,'db2');
for i=1:3  %此处的循环可使编程简洁
   decmpa=wrcoef('a',c,l,'db2',4-i);  
   subplot(4,2,2*i+1);  
   plot(decmpa);  
   Ylabel(['a',num2str(4-i)]);  
end  
for i=1:3  
   decmpd=wrcoef('d',c,l,'db2',4-i);  
   subplot(4,2,2*(i+1));  
   plot(decmpd);  
   Ylabel(['d',num2str(4-i)]);  
end