%裝載leleccum信號
load leleccum;  
s = leleccum(1:3920); 
%用小波函數db1對信號進行三尺度小波分解
[C,L]=wavedec(s,6,'db1');   
subplot(3,1,1);
plot(s);
title('leleccum原始信號');
%用小波函數db1進行信號的重構
X = waverec(C,L,'db1');   
subplot(3,1,2);
plot(X);
title('小波重構');
diffence = s-X;
subplot(3,1,3);
plot(diffence);
title('小波重構誤差');