%裝載信號
load leleccum;  
%截取信號中的一段[2600，3100]
s=leleccum(2600:3100);  
%用小波db3對s進行三層分解
[c,l]=wavedec(s,3,'db3');  
%選用全局閾值進行信號壓縮處理
thr=40;  
[sd,csd,lsd,perf0,perfl2]=wdencmp('gbl',c,l,'db3',3,thr,'h',1);  
subplot(2,1,1);  
plot(s,'LineWidth',2);  
title('原始信號');  
xlabel('樣本序號 n');
ylabel('幅值 A');
subplot(2,1,2);  
plot(sd,'LineWidth',2);  
title('壓縮后的信號'); 
xlabel('樣本序號 n');
ylabel('幅值 A');