m0804.m

MATLAB仿真技术与实例应用教程

load nearbrk;  %载入一原始信号  
s=nearbrk;  

%画出该信号的波形图  
subplot(4,2,1);  
plot(s);  
Ylabel('s');  
title('原始信号s');  
 
[c,l]=wavedec(s,4,'haar');%用harr小波进行3层分解 
for i=1:4  %此处的循环可使编程简洁
   decmpa=wrcoef('a',c,l,'haar',5-i);  
   subplot(5,2,2*i+1);  
   plot(decmpa);  
   Ylabel(['a',num2str(5-i)]);  
end  
for i=1:4  
   decmpd=wrcoef('d',c,l,'haar',5-i);  
   subplot(5,2,2*(i+1));  
   plot(decmpd);  
   Ylabel(['d',num2str(5-i)]);  
end