load noissin;%讀入白噪聲
s=noissin(1:1000);%去信號前1000個采樣點
[c,l]=wavedec(s,3,'db4');%對信號做層數為3的多尺度分解
[cd1,cd2,cd3]=detcoef(c,l,[1,2,3]);%得到三個尺度的細節數
ca3=appcoef(c,l,'db4',3);%得到尺度3的近似數
subplot(511);plot(1:1000,s);title('s');
subplot(512);plot(1:l(1),ca3);title('ca3');
subplot(513);plot(1:l(2),cd3);title('cd3');%l(i)存放的是相應參數的長度
subplot(514);plot(1:l(3),cd2);title('cd2');
subplot(515);plot(1:l(4),cd1);title('cd1');
%將原始信號和分解后得到的以一組近似系數和3組細節數的波形顯示出來