x=0:0.002:2*pi; y=sin(x.^2); r=0.1*randn(size(x));
y1=y+r; subplot(211); plot(x,y1)   % 繪制原始信號
[cA,cD]=dwt(y1,'db4');             % 離散小波變換
subplot(223), plot(cA); subplot(224), plot(cD)  % 繪制近似和細節信號

y2=idwt(cA,cD,'db4'); norm(y1-y2)  % 檢驗反變換對信號的還原程度